Hostname: page-component-586b7cd67f-dsjbd Total loading time: 0 Render date: 2024-12-08T21:53:31.973Z Has data issue: false hasContentIssue false

Identifying the metabolomic fingerprint of high and low flavonoid consumers

Published online by Cambridge University Press:  14 July 2017

Kerry L. Ivey*
Affiliation:
Department of Nutrition, Harvard T. H. Chan School of Public Health, Boston, MA, USA Department of Epidemiology, Harvard T. H. Chan School of Public Health, Boston, MA, USA
Eric B. Rimm
Affiliation:
Department of Nutrition, Harvard T. H. Chan School of Public Health, Boston, MA, USA Department of Epidemiology, Harvard T. H. Chan School of Public Health, Boston, MA, USA Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
Peter Kraft
Affiliation:
Department of Epidemiology, Harvard T. H. Chan School of Public Health, Boston, MA, USA Department of Biostatistics, Harvard T. H. Chan School of Public Health, Boston, MA, USA
Clary B. Clish
Affiliation:
Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA Broad Institute of the Massachusetts Institute of Technology, Boston, MA, USA
Aedin Cassidy
Affiliation:
Department of Nutrition, Norwich Medical School, University of East Anglia, Norwich, Norfolk, UK
Jonathan Hodgson
Affiliation:
University of Western Australia, School of Medicine and Pharmacology, Perth, WA, Australia
Kevin Croft
Affiliation:
University of Western Australia, School of Medicine and Pharmacology, Perth, WA, Australia
Brian Wolpin
Affiliation:
Channing Division of Network Medicine, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA Department of Medical Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, USA
Liming Liang
Affiliation:
Department of Epidemiology, Harvard T. H. Chan School of Public Health, Boston, MA, USA Department of Biostatistics, Harvard T. H. Chan School of Public Health, Boston, MA, USA
*
* Corresponding author: K. L. Ivey, email kivey@hsph.harvard.edu

Abstract

High flavonoid consumption can improve vascular health. Exploring flavonoid–metabolome relationships in population-based settings is challenging, as: (i) there are numerous confounders of the flavonoid–metabolome relationship; and (ii) the set of dependent metabolite variables are inter-related, highly variable and multidimensional. The Metabolite Fingerprint Score has been developed as a means of approaching such data. This study aims to compare its performance with that of more traditional methods, in identifying the metabolomic fingerprint of high and low flavonoid consumers. This study did not aim to identify biomarkers of intake, but rather to explore how systemic metabolism differs in high and low flavonoid consumers. Using liquid chromatography–tandem MS, 174 circulating plasma metabolites were profiled in 584 men and women who had complete flavonoid intake assessment. Participants were randomised to one of two datasets: (a) training dataset, to determine the models for the discrimination variables (n 399); and (b) validation dataset, to test the capacity of the variables to differentiate higher from lower total flavonoid consumers (n 185). The stepwise and full canonical variables did not discriminate in the validation dataset. The Metabolite Fingerprint Score successfully identified a unique pattern of metabolites that discriminated high from low flavonoid consumers in the validation dataset in a multivariate-adjusted setting, and provides insight into the relationship of flavonoids with systemic lipid metabolism. Given increasing use of metabolomics data in dietary association studies, and the difficulty in validating findings using untargeted metabolomics, this paper is of timely importance to the field of nutrition. However, further validation studies are required.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © The Author(s) 2017

Metabolomics is the study of the complement of metabolites present in biological samples( Reference Jenkins, Hardy and Beckmann 1 , Reference Fiehn 2 ). Application of metabolomic technology to targeted mechanistic investigations has elucidated many mechanisms underlying metabolic and disease pathways( Reference Bro, Kamstrup-Nielsen and Engelsen 3 Reference Mamas, Dunn and Neyses 6 ). With technological improvements, metabolomic analyses of human samples have now extended to exploratory, population-based, untargeted analyses, where the aim is to identify new or novel biomarkers of disease. Such analyses are typically conducted in well-designed, nested case–control and cohort studies( Reference Wang, Larson and Vasan 7 , Reference Kenny, Broadhurst and Dunn 8 ).

With increasing availability of metabolome data in large population-based studies with measures of habitual dietary intake, focus has shifted towards identifying metabolic profiles that are associated with different human exposures, such as diet and lifestyle( Reference Altmaier, Kastenmüller and Römisch-Margl 9 , Reference Kuo, Wang and Tian 10 ). Flavonoids are biologically active compounds found in many different foods such as blueberries, tea, wine and chocolate( Reference Harnly, Doherty and Beecher 11 , Reference Beecher 12 ). Given the short half-life of the parent flavonoid compounds( Reference Hollman, Gaag and Mengelers 13 , Reference Erlund, Meririnne and Alfthan 14 ), and their complex metabolism within the gastrointestinal tract( Reference Déprez, Brezillon and Rabot 15 Reference Stoupi, Williamson and Drynan 17 ), there are no validated and reliable biomarkers of total flavonoid intake. Furthermore, although having demonstrated vascular benefits, such as improved nitric oxide homeostasis and endothelial function( Reference Bondonno, Yang and Croft 18 Reference Schroeter, Heiss and Balzer 20 ), the mechanisms underlying the biological effects of flavonoids are complex and yet to be fully elucidated. As such, flavonoids represent an ideal candidate for application to untargeted exposure–metabolome analyses, where flavonoid intake is the independent variable, and the metabolites, the dependent variables, with the aim of identifying how systemic metabolism differs in high and low flavonoid consumers.

Unlike previous flavonoid–metabolome studies( Reference Manach, Hubert and Llorach 21 Reference Llorach, Urpi-Sarda and Jauregui 24 ), this project did not aim, nor did it have the capacity, to identify biomarkers of flavonoid intake or measure acute effects of flavonoid consumption. In fact, flavonoid compounds were not identified in the analytic platforms used in this study. Instead we aimed to describe the long-term metabolic profiles characteristic of those with habitual high or low total flavonoid intake. With this analysis comes the unique challenge of integrating multiple confounding variables into models where the dependent variables consist of many inter-related metabolite concentrations. The non-independence of the dependent metabolite variables makes the translation of results from metabolome-wide association studies challenging. Canonical discriminant analysis is a far more suitable modelling approach for such studies as it is capable of reducing the dimensionality of the rich dependent data matrix into a linear combination of the metabolite variables that summarise the metabolite variation between different levels of flavonoid intake.

Although reducing the dimensionality of the data, full canonical discriminant analysis incorporates information from all the dependent metabolite variables. Therefore, this analytical approach is not useful in identifying which metabolites are most important in discriminating differing levels of exposure in a population. By applying stepwise feature selection to canonical modelling we can identify the combination of metabolites that best discriminate people with differing flavonoid intake levels. However, the statistical properties of stepwise regression often result in biased parameter estimates.

There is large measurement error in metabolomic assessment, with peak area CV reaching as high as 25 %( Reference Townsend, Clish and Kraft 25 ). As such, metabolome data are better suited to ranking analyses and for discriminating samples with very low and very high metabolite concentrations. Conversely, due to the random error and the large number of parameters innately present in metabolomic datasets, metabolome data are less appropriate for complex model development where each metabolite is assigned a coefficient or weighting variable, such as canonical discriminant analysis. As a result of this over-fitting, we hypothesise that canonical variables created in one cohort (a training cohort) will perform poorly, in terms of discriminatory capacity, when applied to a validation cohort.

We instead put forth that a model based on identifying relative patterns of distinguishing metabolites, and collectively summarising the groups of metabolites that are relatively more or less abundant in individuals with different levels of total flavonoid intake. We term this summary variable the Metabolite Fingerprint Score, and we hypothesise that by identifying biological meaningful patterns that may be subtle, and ignoring the ‘noise’ in the data, the Metabolite Fingerprint Score developed in a training cohort will be able to differentiate high from low flavonoid consumers in a validation cohort. It is the aim of this study to compare the performance of metabolome-wide association studies, full canonical variables, step-wise canonical variables and the Metabolite Fingerprint Score in identifying the metabolomic fingerprint of high and low flavonoid consumers.

Methods

Participants

The Health Professionals Follow-Up Study (HPFS) was initiated in 1986 when 51 529 US men 40–75 years working in health professions completed a mailed biennial questionnaire( Reference Giovannucci, Ascherio and Rimm 26 ). The Nurses’ Health Study (NHS) was established in 1976 when 121 700 female nurses aged 30–55 years completed a mailed biennial questionnaire( Reference Colditz and Hankinson 27 ).

For this analysis we used data from a nested case–control study exploring the plasma metabolite profile associated with risk of pancreatic cancer( Reference Yuan, Clish and Wu 28 ): HPFS (n 237); NHS (n 370). Participants were selected for this study if they had incident pancreatic adenocarcinoma cases diagnosed after blood collection up to 2010 with available plasma and no prior history of cancer. For each case, we randomly selected two controls, matching on cohort, sex, year of birth, smoking status and fasting status (<8 h, ≥8 h). To reduce the influence of subclinical malignancy on plasma metabolite levels, we excluded cases diagnosed within 2 years of blood collection.

The final analysis included 584 participants (192 cases and 392 controls), as they had complete dietary flavonoid intake assessment in conjunction with plasma metabolic profiling, at the same time point. The study was approved by the Human Research Committee at Brigham and Women's Hospital (Boston, MA, USA), and all participants provided consent. For the study design, see Fig. 1.

Fig. 1. Study design.

Dietary intake assessment

In 1994 and 1990, the HPFS and NHS participants, respectively, completed a validated FFQ to determine habitual intake of foods and beverages over the preceding year. In particular, participants reported how often, on average, they consumed each food of a standard portion size. Based on the FFQ responses, as well as a comprehensive flavonoid food composition database, habitual daily total flavonoid intake was estimated using standard methods( Reference Cassidy, O'Reilly and Kay 29 ). Specifically, intake was computed by multiplying the frequency of consumption for a particular portion size by the flavonoid content in that particular food item and then summing the product across all food items. Total flavonoid intake was derived by summing up the five major flavonoid subclasses evaluated in the current analysis: flavonols; flavones; flavanones; flavan-3-ols (including both monomers and polymers); and anthocyanins. Isoflavones were not included in this analysis as they are almost exclusively associated with soya intake, with daily intakes less than 3 mg/d in this population of US adults.

Plasma samples

Blood samples in EDTA tubes were collected in the HPFS from 1993 to 1995, and in heparin tubes for the women in the NHS from 1989 to 1990. Blood samples were collected by participants, mailed overnight on cold packs, and then spun to collect and store plasma (delayed processing)( Reference Wei, Giovannucci and Fuchs 30 , Reference Hankinson, Willett and Manson 31 ).

Metabolite profiling

Plasma metabolites were measured as peak areas by a liquid chromatography–MS metabolomics platform directed by C. B. C. at the Broad Institute of the Massachusetts Institute of Technology and Harvard University (Cambridge, MA, USA). Metabolite profiling methods were developed using reference standards of metabolites to determine chromatographic retention times, MS multiple reaction monitoring transitions, declustering potentials and collision energies( Reference Wang, Larson and Vasan 7 ). We evaluated plasma from ten volunteers with plasma collected simultaneously in heparin and EDTA tubes. For the branched-chain amino acids, Spearman correlation coefficients between heparin and EDTA samples were 0·85 for isoleucine, 0·88 for leucine and 0·95 for valine.

In prior pilot work( Reference Townsend, Clish and Kraft 25 ), we determined that thirty-two metabolites had poor reproducibility in samples with delayed processing, so these metabolites were excluded as they could not be reliably measured in two of the participating cohorts. In the present study, three heparin plasma pools (fifty-seven total quality control (QC) samples) and three EDTA plasma pools (128 total QC samples) were randomly interspersed among participant samples as blinded QC samples.

This study included only identified peaks generated from each untargeted platform. We calculated mean CV for each metabolite across QC plasma pools and set an a priori threshold of ≤25 % for satisfactory reproducibility. Based on this criterion, twenty-two metabolites with mean CV > 25 % were excluded from our analyses. Five metabolites had undetectable levels for >10 % of participants and were also excluded (Supplementary Appendix S1). For other metabolites with ≤10 % undetectable, zero values were recoded to the minimum peak area value for each individual metabolite. A total of 174 metabolites passed QC and were included in the analysis (Supplementary Appendix S2).

Statistical analysis

Prior to analysis all metabolite peak areas were normalised by means of log transformation. Flavonoid intake was divided into tertiles, and metabolites were normalised (mean: 0, standard deviation: ±1) in both the HPFS and NHS cohorts. All data were then merged and participants were randomised to either the training dataset or the validation dataset in a 2 to 1 fashion. To ensure all metabolites have the opportunity to contribute equally to the canonical models, the values were also normalised after randomisation into the training and validation datasets. All analyses were performed with the SAS 9.2 statistical package (SAS Institute, Inc.).

The multivariate-adjusted models included case/control status, cohort, quintiles of energy intake, smoking status, age at blood collection, the Alternative Healthy Eating Index (minus alcohol) score and alcohol consumption.

Metabolome-wide association study

We conducted both unadjusted and multivariate-adjusted Spearman correlations between flavonoid intake and each of the 174 metabolites, in individual models. Using a Bonferroni correction for multiple hypothesis testing( Reference Bland and Altman 32 ), only those metabolites with P ≤ 0·0003 (0·05/174) were considered to be statistically significantly associated with total flavonoid intake.

Development of discriminatory variables in the training dataset

The stepwise canonical variable was computed from the canonical coefficients of a stepwise canonical correlation analysis. The independent variable was high (highest tertile) or low (lowest tertile) of flavonoid consumption (moderate flavonoid consumers were excluded from this analysis) and the dependent variable matrix consisted of all 174 metabolites. The significance level for adding variables in the forward selection mode, or removing them in the backward elimination mode, was 0·15. No metabolites were forced into the model.

In a full canonical discriminant analysis, the full canonical variable was computed from the canonical coefficients from the principal components analysis eigenvector encapsulating the greatest degree of discrimination. As such, all 174 metabolites contribute to the computation of the full canonical variable. The Metabolite Fingerprint Score is calculated as follows:

(1) $$\hbox{Metabolite Fingerprint Score} = \displaystyle{{\sum\nolimits_{i = 1}^{10} {h_i}} \over {10}} - \displaystyle{{\sum\nolimits_{i = 1}^{10} {l_i}} \over {10}}$$

where h is the peak area of metabolites with the highest canonical discriminant coefficient, and l is the peak area of metabolites with the lowest canonical discriminant coefficient. Variables included in the Flavonoid Metabolite Fingerprint Score computation are shown in Table 1.

Table 1. Variables included in the Flavonoid Metabolite Fingerprint Score computation

As outlined in Equation 1, the Metabolite Fingerprint Score incorporated the ten metabolites with the largest negative and positive canonical discriminant coefficients, identified by the full canonical discriminant analysis (described above). To explore the influence that fasting status and dietary fat intake have on the discriminatory capacity of the Metabolite Fingerprint Score, we repeated the multivariate-adjusted ANCOVA model, with additional adjustment for fasting status as well as intakes of trans-fat, cholesterol, saturated fat, monounsaturated fat and polyunsaturated fat. Model 1 excluded the Alternate Health Eating Index, and model 2 included the Alternate Health Eating Index. An additional sensitivity analysis was conducted to explore the relationship of fasting status and dietary fat intake with the Metabolite Fingerprint Score, via unadjusted and multivariate-adjusted (partial) Spearman correlation and ANCOVA models.

Evaluating the performance of the discriminatory variables in the validation cohort

Using the models developed in the training dataset, the three discriminatory variables (stepwise canonical variable, full canonical variable, and Metabolite Fingerprint Score) were computed in the validation dataset participants. We then performed separate unadjusted and multivariate-adjusted ANCOVA for each discriminatory variable, with tertiles of flavonoid intake as the independent variable.

Results

Cohort characteristics

The characteristics of the study population are presented in Table 2. The mean daily flavonoid intake was 344 (sd 310) mg/d, with high total flavonoid consumers having an over 5-fold greater mean total flavonoid intake than low total flavonoid consumers. Moderate and high total flavonoid consumers were less likely to have a history of smoking, and moderate consumers were older at the time of blood draw. Following random dichotomisation of the whole cohort, the test and validation cohorts were well matched in terms of total flavonoid intake, as well as other baseline characteristics.

Table 2. Cohort characteristics stratified by level of flavonoid consumption

(Mean values and standard deviations; numbers and percentages)

AHEI, Alternative Healthy Eating Index.

Metabolome-wide association study in the whole cohort

Using the significance cut-off of P ≤ 0·0003, in unadjusted analysis, the only metabolite significantly associated with flavonoid intake was cotinine; a metabolite of nicotine and a biomarker of tobacco smoke exposure( Reference Benowitz 33 ) (Spearman correlation coefficient −0·153; P = 0·0002). However, this relationship did not remain in the partial Spearman rank analysis which controlled for the variables in the multivariate models (Fig. 2).

Fig. 2. Multivariate-adjusted metabolome-wide association study of flavonoid intake and the twenty metabolites to which it is most strongly associated. -----, Bonferroni-corrected level of significance required, after accounting for 174 multiple comparisons. Multivariate-adjusted model includes case/control status, cohort, quintiles of energy intake, smoking status, age at blood collection, the Alternative Healthy Eating Index (minus alcohol) score and alcohol consumption.

We then repeated the partial Spearman rank analysis, with cigarette smoking as the independent variable, and controlled for flavonoid intake as well as the multivariate variables, and observed a strong positive association of cigarette smoking with cotinine concentration (P < 0·0001). This result highlights the importance of confounders in exposure–metabolome association studies, and, by observing the established causal association between cigarette exposure and cotinine concentrations, it validates our capacity to observe biologically meaningful exposure–metabolome associations in this cohort.

In the flavonoid–metabolome partial Spearman correlation analysis, no metabolites were significantly associated with flavonoid intake (Fig. 2).

Development of discriminatory variables in the training dataset

The training dataset was used to compute three discriminatory variables: (1) stepwise canonical variable; (2) full canonical variable; and (3) Metabolite Fingerprint Score.

Computing the stepwise canonical variable

`The forwards–backwards stepwise canonical analysis entered twenty-eight metabolites into the model, four of which were removed. The final model included the following metabolites: cotinine; 4-pyridoxate; thiamine; indoxylsulfate; isocitrate; β-hydoxybutyrate; aconitate; putrescine; sorbitol; glutamate; dimethylglycine; hydroxyphenylacetate; aminoisobutyric acid; quinolinate; sphingomyelin C16 : 0; TAG C48 : 1; TAG C56 : 2; TAG C54 : 8; TAG C50 : 4; TAG C50 : 0; cholesteryl ester C18 : 3; cholesteryl ester C20 : 5; phosphatidylcholine C32 : 2; diacylglycerol C34 : 2.

In the training dataset, the stepwise canonical variable was able to distinguish moderate and high flavonoid consumers from low flavonoid consumers (Table 3).

Table 3. ANCOVA of flavonoid discrimination variables by flavonoid intake group in the training dataset

(Least squared mean values with their standard errors)

* Significantly different from low consumers (P < 0·05).

Canonical value.

Multivariate-adjusted model includes case/control status, cohort, quintiles of energy intake, smoking status, age at blood collection, the Alternative Healthy Eating Index (minus alcohol) score and alcohol consumption.

Computing the full canonical variable

The greatest degree of discrimination between low and high flavonoid consumers was observed along eigenvector 1 of the full canonical discriminant analysis, and therefore, full canonical variable 1 was used for all subsequent analyses.

In the training dataset, the full canonical variable was able to distinguish moderate and high flavonoid consumers from low flavonoid consumers (Table 3).

Computing the Metabolite Fingerprint Score

The Metabolite Fingerprint Score was computed using the results of the full canonical model to dictate which metabolites are included in the computation (Equation 1). The ten ‘high’ and ten ‘low’ metabolites contributing to the Metabolite Fingerprint Score computation were all lipid metabolites. There were six TAG metabolites in the ‘high’ group, and five in the ‘low’ group. The other metabolite groups contributing to both the ‘high’ and ‘low’ metabolite groups were sphingomyelins, diacylglycerols and lysophosphatidylcholines. Phosphatidylcholine was unique to the ‘high’ metabolite group, and lysophosphatidylethanolamine was unique to the ‘low’ metabolite group.

In the training dataset, the Metabolite Fingerprint Score was able to distinguish high flavonoid consumers from low flavonoid consumers (P < 0·05) (Fig. 3).

Fig. 3. ANCOVA of the Metabolite Fingerprint Score by flavonoid intake group in the training dataset. Results are least squared mean values, with their standard errors represented by horizontal bars. * Significantly different from low consumers (P < 0·05). The multivariate-adjusted model includes case/control status, cohort, quintiles of energy intake, smoking status, age at blood collection, the Alternative Healthy Eating Index (minus alcohol) score and alcohol consumption. Low intake, n 123; moderate intake, n 147; high intake, n 129.

Assessing the performance of the discriminatory variables in the validation dataset

After developing the three discriminatory variables in the training dataset, we then sought to test their performance in the validation dataset. We were also interested in exploring the effect of multivariate adjustment on the capacity of the discriminatory variables to discriminate high from low flavonoid consumers.

There was no difference in the value of the stepwise canonical variable when comparing high flavonoid consumers with low flavonoid consumers, in either unadjusted or multivariate-adjusted models. Similarly, no discrimination was provided by the full canonical variable in either unadjusted or multivariate-adjusted models (Table 4).

Table 4. ANCOVA of stepwise canonical variable and full canonical variable by flavonoid intake group in the validation dataset*

(Least squared mean values with their standard errors)

* No values were significantly different from low consumers (P < 0·05).

Canonical value.

Multivariate-adjusted model includes case/control status, cohort, quintiles of energy intake, smoking status, age at blood collection, the Alternative Healthy Eating Index (minus alcohol) score and alcohol consumption.

Conversely, when compared with low flavonoid consumers, the value of the Metabolite Fingerprint Score was higher in high flavonoid consumers. This discriminant capacity was evident in both unadjusted and multivariate-adjusted models (P < 0·05) (Fig. 4).

Fig. 4. ANCOVA of the Metabolite Fingerprint Score by flavonoid intake group in the validation dataset. Results are least squared mean values, with their standard errors represented by horizontal bars. * Significantly different from low consumers (P < 0·05). The multivariate-adjusted model includes case/control status, cohort, quintiles of energy intake, smoking status, age at blood collection, the Alternative Healthy Eating Index (minus alcohol) score and alcohol consumption. Low intake, n 60; moderate intake, n 67; high intake, n 58.

Exploring the lipid nature of the Metabolite Fingerprint Score

As can be seen in Equation 1, the Metabolite Fingerprint Score is largely comprised of lipid metabolites, predominantly TAG, which appear on both the high and the low metabolite lists. To explore the potential role of dietary fat intake in contributing to this pattern, we repeated the validation analysis, adjusting for the variables in the multivariate-adjusted model (excluding Alternate Health Eating Index), with additional adjustment for intakes of trans-fat, cholesterol, saturated fat, monounsaturated fat and polyunsaturated fat. This did not influence the significance or strength of discrimination (multivariate adjusted least squared means of low and high consumers: −0·09 (sem 0·04) and +0·01 (sem 0·04), respectively; P = 0·0185). Furthermore, none of the dietary fats contributed significantly to the model (P > 0·05). Results were similar for the models which also included the Alternate Healthy Eating Index. Furthermore, discrimination capacity was not altered when fasting status was included in the multivariate-adjusted dietary fat model (multivariate adjusted least squared means of low and high consumers: −0·09 (sem 0·04) and +0·01 (sem 0·04), respectively; P = 0·0161).

We then sought to explore the independent association of dietary fat intake with the Metabolite Fingerprint Score. None of the dietary fats, nor fasting status, was correlated with the Metabolite Fingerprint Score in either unadjusted or partial (multivariate-adjusted) Spearman correlation models (multivariate-adjusted ANCOVA, P > 0·05). Furthermore, for each dietary fat, the Metabolite Fingerprint Score did not discriminate high from low consumers. This was also the case for fasting status, where the Metabolite Fingerprint Score was not significantly different in participants who did and did not fast (multivariate-adjusted ANCOVA, P > 0·05).

Exploring the Metabolite Fingerprint Score performance in the validation dataset

As the Metabolite Fingerprint Score is conceptually and statistically the dependent variable in these analyses, we are not concerned with whether resultant models explain the Metabolic Fingerprint Score, but rather whether the Metabolic Fingerprint Score is different in high and low flavonoid consumers, in the context of a multivariate setting. As such, to determine if the score consisting of ten high and ten low metabolites provides the greatest discrimination of high and low flavonoid consumers, we compared the F statistic assigned to the flavonoid consumption variable (high v. low intake) in a multivariate-adjusted ANCOVA. This was done for Metabolite Fingerprint Scores comprised of between five and fifteen metabolites, with the null hypothesis that in a multivariate-adjusted setting, the mean Metabolic Fingerprint Score of high and low flavonoid consumers would be the same.

The value of the F statistic ranged from 1·51 for the Metabolite Fingerprint Score computed from five high and five low metabolites, to 0·07 for the Metabolite Fingerprint Score computed from fifteen high and fifteen low metabolites. With an F statistic of 4·75, well above the critical threshold of 3·9, the Metabolite Fingerprint Score computed from ten high and ten low metabolites had the greatest discriminatory capacity.

Discussion

Using models created in a training dataset, when applied to a separate validation dataset containing different participants, the Metabolite Fingerprint Score was able to distinguish high from low flavonoid consumers, even after taking into account potential confounders of the flavonoid–metabolome relationship. Conversely, both the stepwise canonical variable and the full canonical variable did not distinguish high from low flavonoid consumers in the validation dataset. In addition, the multivariate-adjusted metabolome-wide association study failed to identify any significant metabolites.

The unadjusted metabolome-wide association study identified a significant, inverse, association between flavonoid intake and cotinine, a biomarker for cigarette smoking( Reference Benowitz 33 ). This spurious association manifested because of a strong positive association between cigarette smoking and cotinine peak area, and a strong inverse association between flavonoid intake and smoking status. This association disappeared following multivariate-adjustment, which included adjusting for smoking status, thus highlighting the importance of multivariate adjustment in diet–metabolome studies. Consistent with our hypothesis, the multivariate-adjusted metabolome-wide association study failed to identify any metabolites that were significantly associated with flavonoid consumption. As such, we pursued a pattern-based approach to analysing the data.

The Metabolite Fingerprint Score incorporated a set of ten ‘high’ metabolites that were in higher concentrations in high flavonoid consumers, and lower concentrations in low flavonoid consumers. An additional set of ten different metabolites were also incorporated into the Metabolite Fingerprint Score computation. These ‘low’ metabolites represented those that were in lower concentrations in high flavonoid consumers, and higher concentrations in low flavonoid consumers. This pattern consisting of ten ‘high’ and ten ‘low’ metabolites was found to provide the greatest degree of differentiation in the validation dataset.

This study utilised data from both lipid-soluble and polar metabolomic platforms. The stepwise and full canonical variables incorporated both lipid-soluble and polar metabolites, many of which included metabolites with large CV. Conversely, the Metabolite Fingerprint Score included only lipid metabolites that had a low degree of random measurement error (mean CV of metabolites included in the Metabolite Fingerprint Score: 9·6 %). If metabolites, or their pattern, are biologically related to flavonoid intake, these associations are more likely to be picked up in metabolites which can be assessed with a high degree of precision. As such, the preferential incorporation of low CV metabolites into the score was not by design, but rather because of the properties of the Metabolite Fingerprint Score itself.

Due to the dominance of lipid metabolites in the Metabolite Fingerprint Score, we thoroughly explored the role of dietary fat intake in influencing observed results. The inclusion of fasting status as well as intakes of cholesterol, trans-fat, saturated fat, unsaturated fat and polyunsaturated fat into the multivariate-adjusted model did not influence the discriminatory capacity of the Metabolite Fingerprint Score. Furthermore, there was no independent association of fasting status or any the dietary fats with the Metabolite Fingerprint Score. This is not surprising given the equal contribution of lipids to both the ‘high’, and the ‘low’, metabolite groups. These results suggest that dietary fat intake is not driving the observed results, and that the composition of the Metabolite Fingerprint Score may instead be reflecting the underlying differences in fatty acid metabolism amongst high and low flavonoid consumers. It is important to note that the lipid metabolite profile provides information on how lipids are metabolised and created systemically, which is quite distinct from the clinical relationship of flavonoids with total cholesterol, LDL-cholesterol, HDL-cholesterol and total TAG. Although there is a great degree of heterogeneity across studies, there is evidence to suggest possible effects of specific flavonoid compounds on some clinical lipid measurements( Reference Shrime, Bauer and McDonald 34 ). Further studies are warranted to fully elucidate the effect that flavonoids have on systemic lipid metabolism.

While it is beyond the scope of this paper to relate the lipid metabolite compounds to potential bioactivity, and the biological interpretation of results arising from lipid metabolomics platforms can be challenging, there are numerous reasons why one might expect the Metabolite Fingerprint Score to be predominated by lipids. It is hypothesised that due to humans’ capacity for fatty acid storage( Reference Ali, Koutsari and Mundi 35 ), lipid-soluble metabolites may be a better indicator of long-term alterations in metabolism, such as those explored in this study, when compared with the polar metabolites which often have a high systemic turnover rate. However, the most likely explanation comes from the nature of the metabolome data themselves. With a mean CV of 11·6 % for the lipid-soluble metabolites, and 19·4 % for the polar metabolites, our capacity to reproducibly assess metabolites was much greater for lipid-soluble metabolites in comparison with their water-soluble counterparts.

The Metabolite Fingerprint Score has identified a pattern of two sets of ten metabolites that, in relation to the metabolome as a whole, were in higher and lower concentrations in high flavonoid consumers, when compared with their low flavonoid-consuming counterparts. Furthermore, in a separate validation dataset, we were able to confirm that this pattern of metabolites characterised high from low flavonoid consumers, in individuals distinct from those upon which the score was created. This internal validation of the Metabolite Fingerprint Score is an important first step in understanding the global metabolic pathways affected by high flavonoid consumption. However, unlike polar metabolites which are well described in many common metabolic pathway analysis tools and databases, the role of fatty acids in systemic metabolism is less well characterised. As such, the results of population-based studies such as these are most useful in identifying candidate compounds which can then be further scrutinised in clinical trial settings, as well as by the use of more targeted spectroscopic methods, such as NMR, where the precise fatty acid structure can be scrutinised, and biological significance of this pattern can be unravelled.

Previous diet–plasma metabolome studies, using various study designs, have identified diet–metabolome relationships using a variety of analytical methods; ranging from correlation analysis to cluster analysis( Reference O'Sullivan, Gibney and Brennan 36 , Reference Menni, Zhai and MacGregor 37 ). However, it remains uncertain which is the optimal statistical method to apply to diet–metabolome analyses. Despite all three discrimination variables discriminating high from low flavonoid consumers in the dataset in which the models were developed, only the Metabolite Fingerprint Score discriminated high from low flavonoid consumers in the validation dataset. By weighting all included metabolites equally, and incorporating only directionality into the model, the Metabolite Fingerprint Score was able to avoid carrying forward artifacts, or unintended patterns, from the training dataset.

In conclusion, the Metabolite Fingerprint Score was the most valid means of identifying a pattern of inter-related metabolites that differentiated high from low flavonoid consumers, in a multivariate setting. However, to confirm the superior performance of the Metabolite Fingerprint Score, further validation studies using external cohorts, different exposures, and different metabolomics platforms are required. The composition of the Metabolite Fingerprint Score provides insight into how systemic lipid metabolism differs across high and low flavonoid consumers.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/jns.2017.27

Acknowledgements

This work was supported by a Nestlé Research Center award to the Broad Institute; and by the Robert T. and Judith B. Hale Fund for Pancreatic Cancer, Perry S. Levy Fund for Gastrointestinal Cancer Research, Pappas Family Research Fund for Pancreatic Cancer (National Institutes of Health (NIH) grant R01 CA124908), and NIH P50 CA127003 to CS Fuchs; and by the Department of Defense (CA130288), Howard Hughes Medical Institute, Lustgarten Foundation, and Promises for Purple to B. W. The HPFS is supported by NIH grant UM1 CA167552. The NHS is supported by NIH grants UM1 CA186107, P01 CA87969 and R01 CA49449. The salary of K. L. I. is supported by a National Health and Medical Research Council early career fellowship.

All authors made a substantial contribution to the protocol design, analysis and manuscript revisions. B. W. and C. B. C. also provided expertise on sample collection.

E. B. R., K. L. I. and A. C. reported receiving research support from the US Highbush Blueberry Research Council, for projects unrelated to this analysis. No other disclosures or conflicts of interest were reported.

References

1. Jenkins, H, Hardy, N, Beckmann, M, et al. (2004) A proposed framework for the description of plant metabolomics experiments and their results. Nat Biotechnol 22, 16011606.Google Scholar
2. Fiehn, O (2002) Metabolomics – the link between genotypes and phenotypes. Plant Mol Biol 48, 155171.Google Scholar
3. Bro, R, Kamstrup-Nielsen, MH, Engelsen, SB, et al. (2015) Forecasting individual breast cancer risk using plasma metabolomics and biocontours. Metabolomics 11, 13761380.Google Scholar
4. Kuhara, T (2007) Noninvasive human metabolome analysis for differential diagnosis of inborn errors of metabolism. J Chromatogr B Analyt Technol Biomed Life Sci 855, 4250.Google Scholar
5. Kind, T, Tolstikov, V, Fiehn, O, et al. (2007) A comprehensive urinary metabolomic approach for identifying kidney cancer. Anal Biochem 363, 185195.Google Scholar
6. Mamas, M, Dunn, WB, Neyses, L, et al. (2011) The role of metabolites and metabolomics in clinically applicable biomarkers of disease. Arch Toxicol 85, 517.Google Scholar
7. Wang, TJ, Larson, MG, Vasan, RS, et al. (2011) Metabolite profiles and the risk of developing diabetes. Nature Medicine 17, 448453.Google Scholar
8. Kenny, LC, Broadhurst, DI, Dunn, W, et al. (2010) Robust early pregnancy prediction of later preeclampsia using metabolomic biomarkers. Hypertension 56, 741749.Google Scholar
9. Altmaier, E, Kastenmüller, G, Römisch-Margl, W, et al. (2011) Questionnaire-based self-reported nutrition habits associate with serum metabolism as revealed by quantitative targeted metabolomics. Eur J Epidemiol 26, 145156.Google Scholar
10. Kuo, C-H, Wang, K-C, Tian, T-F, et al. (2012) Metabolomic characterization of laborers exposed to welding fumes. Chem Res Toxicol 25, 676686.CrossRefGoogle Scholar
11. Harnly, JM, Doherty, RF, Beecher, GR, et al. (2006) Flavonoid content of US fruits, vegetables, and nuts. J Agric Food Chem 54, 99669977.CrossRefGoogle ScholarPubMed
12. Beecher, GR (2003) Overview of dietary flavonoids: nomenclature, occurrence and intake. J Nutr 133, 3248S3254S.Google Scholar
13. Hollman, PC, Gaag, MV, Mengelers, MJ, et al. (1996) Absorption and disposition kinetics of the dietary antioxidant quercetin in man. Free Radic Biol Med 21, 703707.Google Scholar
14. Erlund, I, Meririnne, E, Alfthan, G, et al. (2001) Plasma kinetics and urinary excretion of the flavanones naringenin and hesperetin in humans after ingestion of orange juice and grapefruit juice. J Nutr 131, 235241.Google Scholar
15. Déprez, S, Brezillon, C, Rabot, S, et al. (2000) Polymeric proanthocyanidins are catabolized by human colonic microflora into low-molecular-weight phenolic acids. J Nutr 130, 27332738.CrossRefGoogle ScholarPubMed
16. Aura, AM, Martin-Lopez, P, O'Leary, KA, et al. (2005) In vitro metabolism of anthocyanins by human gut microflora. Eur J Nutr 44, 133142.Google Scholar
17. Stoupi, S, Williamson, G, Drynan, JW, et al. (2010) A comparison of the in vitro biotransformation of (–)-epicatechin and procyanidin B2 by human faecal microbiota. Mol Nutr Food Res 54, 747759.Google Scholar
18. Bondonno, CP, Yang, X, Croft, KD, et al. (2012) Flavonoid-rich apples and nitrate-rich spinach augment nitric oxide status and improve endothelial function in healthy men and women: a randomized controlled trial. Free Radic Biol Med 52, 95102.Google Scholar
19. Loke, WM, Hodgson, JM, Proudfoot, JM, et al. (2008) Pure dietary flavonoids quercetin and (−)-epicatechin augment nitric oxide products and reduce endothelin-1 acutely in healthy men. Am J Clin Nutr 88, 10181025.Google Scholar
20. Schroeter, H, Heiss, C, Balzer, J, et al. (2006) (–)-Epicatechin mediates beneficial effects of flavanol-rich cocoa on vascular function in humans. Proc Natl Acad Sci U S A 103, 10241029.Google Scholar
21. Manach, C, Hubert, J, Llorach, R, et al. (2009) The complex links between dietary phytochemicals and human health deciphered by metabolomics. Mol Nutr Food Res 53, 13031315.Google Scholar
22. Walsh, MC, Brennan, L, Pujos-Guillot, E, et al. (2007) Influence of acute phytochemical intake on human urinary metabolomic profiles. Am J Clin Nutr 86, 16871693.Google Scholar
23. Jacobs, DM, Fuhrmann, JC, van Dorsten, FA, et al. (2012) Impact of short-term intake of red wine and grape polyphenol extract on the human metabolome. J Agric Food Chem 60, 30783085.Google Scholar
24. Llorach, R, Urpi-Sarda, M, Jauregui, O, et al. (2009) An LC-MS-based metabolomics approach for exploring urinary metabolome modifications after cocoa consumption. J Proteome Res 8, 50605068.Google Scholar
25. Townsend, MK, Clish, CB, Kraft, P, et al. (2013) Reproducibility of metabolomic profiles among men and women in 2 large cohort studies. Clin Chem 59, 16571667.Google Scholar
26. Giovannucci, E, Ascherio, A, Rimm, EB, et al. (1995) Physical activity, obesity, and risk for colon cancer and adenoma in men. Ann Intern Med 122, 327334.Google Scholar
27. Colditz, GA & Hankinson, SE (2005) The Nurses’ Health Study: lifestyle and health among women. Nat Rev Cancer 5, 388396.Google Scholar
28. Yuan, C, Clish, CB, Wu, C, et al. (2016) Circulating metabolites and survival among patients with pancreatic cancer. J Natl Cancer Inst 108, djv409.Google Scholar
29. Cassidy, A, O'Reilly, ÉJ, Kay, C, et al. (2011) Habitual intake of flavonoid subclasses and incident hypertension in adults. Am J Clin Nutr 93, 338347.Google Scholar
30. Wei, EK, Giovannucci, E, Fuchs, CS, et al. (2005) Low plasma adiponectin levels and risk of colorectal cancer in men: a prospective study. J Natl Cancer Inst 97, 16881694.Google Scholar
31. Hankinson, SE, Willett, WC, Manson, JE, et al. (1995) Alcohol, height, and adiposity in relation to estrogen and prolactin levels in postmenopausal women. J Natl Cancer Inst 87, 12971302.Google Scholar
32. Bland, JM & Altman, DG (1995) Multiple significance tests: the Bonferroni method. BMJ 310, 170.Google Scholar
33. Benowitz, NL (1996) Cotinine as a biomarker of environmental tobacco smoke exposure. Epidemiol Rev 18, 188204.Google Scholar
34. Shrime, MG, Bauer, SR, McDonald, AC, et al. (2011) Flavonoid-rich cocoa consumption affects multiple cardiovascular risk factors in a meta-analysis of short-term studies. J Nutr 141, 19821988.Google Scholar
35. Ali, AH, Koutsari, C, Mundi, M, et al. (2011) Free fatty acid storage in human visceral and subcutaneous adipose tissue: role of adipocyte proteins. Diabetes 60, 23002307.Google Scholar
36. O'Sullivan, A, Gibney, MJ & Brennan, L (2011) Dietary intake patterns are reflected in metabolomic profiles: potential role in dietary assessment studies. Am J Clin Nutr 93, 314321.Google Scholar
37. Menni, C, Zhai, G, MacGregor, A, et al. (2013) Targeted metabolomics profiles are strongly correlated with nutritional patterns in women. Metabolomics 9, 506514.Google Scholar
Figure 0

Fig. 1. Study design.

Figure 1

Table 1. Variables included in the Flavonoid Metabolite Fingerprint Score computation

Figure 2

Table 2. Cohort characteristics stratified by level of flavonoid consumption(Mean values and standard deviations; numbers and percentages)

Figure 3

Fig. 2. Multivariate-adjusted metabolome-wide association study of flavonoid intake and the twenty metabolites to which it is most strongly associated. -----, Bonferroni-corrected level of significance required, after accounting for 174 multiple comparisons. Multivariate-adjusted model includes case/control status, cohort, quintiles of energy intake, smoking status, age at blood collection, the Alternative Healthy Eating Index (minus alcohol) score and alcohol consumption.

Figure 4

Table 3. ANCOVA of flavonoid discrimination variables by flavonoid intake group in the training dataset(Least squared mean values with their standard errors)

Figure 5

Fig. 3. ANCOVA of the Metabolite Fingerprint Score by flavonoid intake group in the training dataset. Results are least squared mean values, with their standard errors represented by horizontal bars. * Significantly different from low consumers (P < 0·05). The multivariate-adjusted model includes case/control status, cohort, quintiles of energy intake, smoking status, age at blood collection, the Alternative Healthy Eating Index (minus alcohol) score and alcohol consumption. Low intake, n 123; moderate intake, n 147; high intake, n 129.

Figure 6

Table 4. ANCOVA of stepwise canonical variable and full canonical variable by flavonoid intake group in the validation dataset*(Least squared mean values with their standard errors)

Figure 7

Fig. 4. ANCOVA of the Metabolite Fingerprint Score by flavonoid intake group in the validation dataset. Results are least squared mean values, with their standard errors represented by horizontal bars. * Significantly different from low consumers (P < 0·05). The multivariate-adjusted model includes case/control status, cohort, quintiles of energy intake, smoking status, age at blood collection, the Alternative Healthy Eating Index (minus alcohol) score and alcohol consumption. Low intake, n 60; moderate intake, n 67; high intake, n 58.

Supplementary material: File

Ivey supplementary material

Appendix

Download Ivey supplementary material(File)
File 54.8 KB