Next Article in Journal
Metabolic Effects of Dietary Glycerol Supplementation in Muscle and Liver of European Seabass and Rainbow Trout by 1H NMR Metabolomics
Previous Article in Journal
The metaRbolomics Toolbox in Bioconductor and beyond
Previous Article in Special Issue
The Biosynthetic Pathway of Major Avenanthramides in Oat
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of 1H-NMR Metabolomics for the Discovery of Blood Plasma Biomarkers of a Mediterranean Diet

1
Institute for Global Food Security, School of Biological Sciences, Queen’s University Belfast, Belfast BT9 5DL, UK
2
Beaumont Health, 3811 W. 13 Mile Road, Royal Oak, MI 48073, USA
3
Oakland University-William Beaumont School of Medicine, Rochester, MI 48309, USA
4
Centre for Public Health, Queen’s University Belfast, Belfast BT12 6BA, UK
5
School of Medicine, Dentistry and Biomedical Sciences, Queen’s University Belfast, Belfast BT9 7BL, UK
*
Author to whom correspondence should be addressed.
Metabolites 2019, 9(10), 201; https://doi.org/10.3390/metabo9100201
Submission received: 29 July 2019 / Revised: 22 September 2019 / Accepted: 23 September 2019 / Published: 27 September 2019
(This article belongs to the Special Issue Molecular Nutrition and Metabolism)

Abstract

:
The Mediterranean diet (MD) is a dietary pattern well-known for its benefits in disease prevention. Monitoring adherence to the MD could be improved by discovery of novel dietary biomarkers. The MEDiterranean Diet in Northern Ireland (MEDDINI) intervention study monitored the adherence of participants to the MD for up to 12 months. This investigation aimed to profile plasma metabolites, correlating each against the MD score of participants (n = 58). Based on an established 14-point scale MD score, subjects were classified into two groups (“low” and “high”). 1H-Nuclear Magnetic Resonance (1H-NMR) metabolomic analysis found that citric acid was the most significant metabolite (p = 5.99 × 10−4*; q = 0.03), differing between ‘low’ and ‘high’. Furthermore, five additional metabolites significantly differed (p < 0.05; q < 0.35) between the two groups. Discriminatory metabolites included: citric acid, pyruvic acid, betaine, mannose, acetic acid and myo-inositol. Additionally, the top five most influential metabolites in multivariate models were also citric acid, pyruvic acid, betaine, mannose and myo-inositol. Metabolites significantly correlated with the consumption of certain food types. For example, citric acid positively correlated fruit, fruit juice and vegetable constituents of the diet, and negatively correlated with sweet foods alone or when combined with carbonated drinks. Citric acid was the best performing biomarker and this was enhanced by paired ratio with pyruvic acid. The present study demonstrates the utility of metabolomic profiling for effectively assessing adherence to MD and the discovery of novel dietary biomarkers.

Graphical Abstract

1. Introduction

The Mediterranean diet (MD) is regarded as a healthy dietary pattern with identified health benefits [1], including protective effects from cardiovascular disease (CVD) [2]. The MD is characterised by a high intake of fruits and vegetables, whole grains, fish, nuts, legumes, olive oil and moderate wine consumption [3]. It is also characterised by low consumption of dairy products, and highly processed foods such as sugar sweetened beverages, processed meats and sweets, fast food and trans-fat. In the PREDIMED (Prevención con Dieta Mediterránea) study, compared with a low-fat diet, an MD supplemented with nuts or virgin olive oil decreased the incidence of CVD [4]. Furthermore, when supplemented with nuts, it improved blood lipoprotein profiles resulting in a lower risk of atherosclerosis [5]. An inverse relationship between MD adherence and the development of metabolic syndrome has also been demonstrated [6]. A cross-sectional investigation of 8821 Eastern Europeans found significant associations between Mediterranean diet score (MDS) and waist circumference, systolic blood pressure and triglycerides [7]. Individuals with an MDS in the highest quartile were significantly less likely to have metabolic syndrome compared with those in the lowest quartile [7]. Similarly, an observational study of 3232 healthy French adults reported that adherence to MD was associated with a reduced risk of metabolic syndrome [8]. Furthermore, the PREDIMED study examined 808 participants at high risk of CVD, finding that subjects in the highest MDS quartile were at the lowest risk of metabolic syndrome, and were less likely to have elevated triglycerides or decreased HDL [9]. Adherence to dietary patterns is currently assessed by self-reported food frequency questionnaires and by 24h food recall. The limitations of these methods include lack of accuracy and potential for bias [10]. For example, fruit and vegetable intake is shown to be overestimated when participants record their own intake [11]. Several dietary scores have been developed and applied to populations to evaluate the role of diet in health. However, to overcome the disadvantages of traditional methods, it may be helpful if these practices are combined with novel biomarkers of food exposure. Combining these could more accurately measure adherence to a given diet [12]. Biomarkers of food intake are any component in the body fluids that can be utilised to indicate dietary exposure. When complemented with traditional dietary assessment practices, they may provide a more accurate picture of individual food exposure [11]. The lack of robust biomarkers of compliance is also a major challenge in clinical nutrition research, hindering the interpretation of clinical trials and observational studies of free populations. Hence, there is a need for reliable, non-biased, convenient and cost-effective methods to overcome this. In a Northern European population, MD adherence over a period of 12 months tended to increase plasma levels of vitamin C, oleic acid and eicosapentaenoic acid [13], but clearly, there is potential for the discovery of biomarkers associated with an MD pattern.
In this regard, metabolomics could be utilised for discovering robust dietary biomarkers [10]. Metabolomics facilitates the accurate measurement of many metabolites in body fluids, making it a useful tool for assessing the nutritional status of an individual [14]. Common metabolomic platforms are based on either 1H-Nuclear Magnetic Resonance (1H-NMR) spectroscopy or mass spectrometry. While the latter is more sensitive, 1H-NMR can provide robust, reproducible and unbiased results with very low pre-treatment requirements [15], making it a suitable platform for large cohort studies.
The majority of previously reported studies investigating adherence to the MD have focused on urine. 1H-NMR metabolomic profiling of human urine found that MD adherence was associated with changes in the levels of 3-hydroxybutyrate, citric acid, and cis-aconitate, oleic acid, suberic acid, various amino acids and some microbial cometabolites [1]. A second study, also involving 1H-NMR profiling of urine, identified 34 metabolites discriminating between low and high adherence to MD [12]. A third NMR study compared urine from individuals consuming an MD (supplemented with Coenzyme Q10) with those consuming a Western diet rich in saturated fat [16]. This found higher hippurate in the MD group, which positively correlated with Coenzyme Q10 and beta-carotene plasma levels. Finally, a fourth study undertook liquid chromatography-mass spectrometry (LC-MS) analysis of urine using a subset of participants from the PREDIMED study [17]. Specifically, this examined metabolomic differences between non-walnut eaters (who never consumed walnuts) and walnut eaters (> 3 servings/week). A total of 18 metabolites were higher in walnut consumers and these included fatty acid metabolites, polyphenol microbial metabolites and compounds related to the tryptophan/serotonin pathway.
Thus far, only two studies have examined the blood metabolome of subjects consuming an MD. The first employed an LC-MS/MS methodology and indicated that the levels of several blood lipids are altered by an MD-based intervention [18]. The latter examined alternate MD correlation with serum metabolites and found that this diet was associated with 49 metabolites which positively correlated with different MD food components [19]. The overall aim of the current study was to discover novel blood-based metabolite biomarkers associated with the MD pattern. Therefore, we performed high-throughput 1H-NMR metabolomic analysis on blood plasma from a controlled study of MD adherence in a Northern European population.

2. Results

2.1. Univariate Analysis

Following metabolomic profiling of 137 plasma samples from 58 subjects collected at three different time points (baseline, 6 months, 12 months), we were able to identify and quantify 59 metabolites in each of the 1H-NMR spectra. Citric acid was the most significant metabolite differing between Low and High MDS (p = 5.99 × 10−4*; q = 0.03). Moreover, five additional metabolites significantly differed (p < 0.05; q < 0.35) between the two groups. These were: pyruvic acid (p = 0.005), betaine (p = 0.014), mannose (p = 0.021), acetic acid (p = 0.030) and myo-inositol (p = 0.035). Citric acid and betaine were significantly higher in the high MDS group by 20.5% and 15.5%, respectively. Mannose, pyruvic acid and myo-inositol were 28.7%, 20.9% and 24.3% lower, respectively (Table 1).
Table 1 shows the univariate statistical analysis of plasma metabolites in patients from Low (0–4) and High (5–10) MDS groups. Average concentrations (µM) of top 11 ranked metabolites and their respective standard deviations (SD). *p < 0.05 (citric acid, myo-inositol, pyruvic acid, mannose and betaine) with concentration averages, false discovery rate (FDR), area under the curve (AUC) of the receiver operating characteristic (ROC) curve and % change.

2.2. Multivariate Analysis

Figure 1A shows principal component analysis (PCA) of the 1H-NMR data. Three potential outliers were highlighted, which prompted us to check the pharmacological profiles of these individuals. The profiles were broadly similar to the other participants and thus none were excluded. Principal component 1 (PC1) explained 29.5% of the variance and component 2 (PC2) explained 8.9% of the variance. Two supervised multivariate methods were then applied, orthogonal partial least squares discriminant analysis (OPLS-DA; Figure 1B) and partial least squares discriminant analysis (PLS-DA; Figure 1C), both of which improved the separation of groups (Figure 1B). The PLS-DA model was then cross-validated on Metaboanalyst using the 10-fold CV cross validation method and using a maximum of three components. (Supplementary Materials Table S4), and by means of a Variable importance in projection (VIP) the top 15 ranking metabolites were identified. Citric acid, myo-inositol, pyruvic acid, mannose and betaine corresponded to the top five most discriminant metabolites in the model (Figure 1D).

2.3. Correlation of Metabolite Concentrations with MDS, and Consumption of Food Components

Citric acid, pyruvic acid, betaine, mannose and myo-inositol data were non-normally distributed and metabolite correlations were examined using Spearman’s rank correlation coefficient. These five metabolites significantly correlated with the 14-point scale MDS. Citric acid and betaine positively correlated with MDS while mannose, pyruvic acid and myo-inositol negatively correlated (Supplementary Materials Table S1) (Figure 2A). Application of multiple comparison using the Benjamini-Hochberg approach showed FDR to be higher than 0.05 (Supplementary Materials Table S1).
Correlations were examined between each of the top five metabolites and each food group (Supplementary Materials Table S1). These were: Fruit (Fr); Fruit Juice (J); Combined Fruit and Fruit Juice (FJ); Combined Fruit, Fruit Juice and Vegetables (FJV); Red meat (RM); combined Chicken and Turkey (CT); Fish (F); Nuts (N); Processed meat (PM), Legumes (L) and Alcohol (A); Olive Spreads (OS); Olive/Rapeseed oil (OR); Sweet Foods (S); Combined sweet foods and carbonated drinks (SD) and whole grain cereals (C). Citric acid positively correlated with J, FJ and FJV and negatively correlated with S and SD. Mannose significantly negatively correlated with J, FJ, FJV and C. Pyruvic acid significantly negatively correlated with J, FJ, FJV, F and OS and significantly positively correlated with PM. Myo-inositol was significantly negatively correlated with OS and positively correlated with PM, betaine significantly positively correlated with J and FJV. Application of multiple comparison using the Benjamini-Hochberg approach showed FDR to be higher than 0.05 (Supplementary Materials Table S1). T-test conducted to evaluate smoker and non-smoker differences showed no metabolites significantly differed between smokers and non-smokers. The top five shortlisted metabolites were citric acid (p = 0.07), betaine (p = 0.18), mannose (p = 0.34), pyruvic acid (p = 0.34) and myo-inositol (p = 0.34). Food group-to-Food group correlations within the study cohort and correlations of all 59 metabolites against each food group scored within the MDS can be seen on supplementary materials (Table S5 and Table S6 respectively)

2.4. Identification of Biomarkers

Table 1 shows the individual area under the curve (AUC) of the receiver operating characteristic (ROC) curve for each computed metabolite. Citric acid performed moderately well with an AUROC curve=0.67 which was superior to the other metabolites: pyruvic acid (0.64), betaine (0.62), mannose (0.61) and myo-inositol (0.60). Computed metabolite ratios showed improved ROC values. The citric acid/pyruvic acid was highest (0.74) (Figure 2C) followed by citric acid/lactic acid (0.73), citric acid/aspartate (0.72), citric acid/l-phenylalanine (0.71) and betaine/pyruvic acid (0.70) (Table 2).
Table 2 shows the best performing paired metabolite ratios with area under the curve (AUC) of the receiver operating characteristic (ROC) curve and % change. Features displayed in the table are ranked based on area under ROC curve with highest AUC (ROC) corresponding to the ratio Citric acid/Pyruvic acid, Citric acid/L-Lactic acid or Citric acid/Aspartate.

3. Discussion

Combining traditional methods such as food diaries or 24-h recall with biomarkers of food exposure may provide a more accurate picture for assessing adherence to a given diet. The present study used food diaries from 58 participants and quantified 59 blood metabolites by 1H-NMR and examined how their concentrations differed according to the reported Mediterranean diet score (MDS) of participants. Scoring systems are useful for measuring dietary patterns in a population but there is a need for continual improvement [20]. A previously validated 14-point item scale was adapted and adopted [3]. It was also employed in a UK population, demonstrating satisfactory accuracy for assessing MD adherence among individuals [21].
Adherence to the MD has often been assessed in Mediterranean populations where the MD patterns are regarded as familiar eating habits. Studying the efficacy of the MD in non-MD populations remains a challenge [22]. MD adherence in non-Mediterranean populations is generally lower [13]. In the MEDDINI cohort no participants achieved the highest MDS. The maximum MDS achieved in this non-Mediterranean population was 10, and “low” (0–4) and “high” (5–10) MDS groups were selected by dividing groups at the median. Unsupervised multivariate models incorporating all 59 metabolites only weakly separated these groups. This was modestly improved by the application of supervised modelling methods such as OPLS-DA and PLS-DA; however, it is clear that a degree of overlap exists between groups when the overall metabolomic profile is considered. This result is not surprising given the diverse range of components contributing to the MDS scoring system and the inherent heterogeneous character of the MD. Despite this, the multivariate statistical analysis identified five individual metabolites within the model, which were important for discriminating between “low” and “high” MDS groups. These were: citric acid, mannose, pyruvic acid, myo-inositol and betaine; therefore, we proceeded with univariate statistical analysis.
These same five metabolites significantly differed between “low” and “high” MDS groups. Citric acid and betaine were significantly higher in the high MDS group, by 20.5% and 15.5%, respectively. Myo-inositol, pyruvic acid and mannose were lower in the high MDS group, by 24.3%, 20.9% and 28.7%, respectively. Area-under-the-curve (AUC) values of the receiver operating characteristic (ROC) curves were compared for every single individual metabolite and all paired combinations of metabolites. Citric acid was the most influential metabolite in discriminating between ‘low and ‘high’ MD groups. As a single biomarker, citric acid was the best performing, it was the most capable of discriminating between “low” and “high” MDS and achieved the highest AUC
ROC curve (0.67). Promisingly, this was greatly enhanced through combination with pyruvic acid. A metabolite ratio of citric acid and pyruvic acid achieved an AUC ROC of 0.74.
In terms of the dietary components, citric acid positively correlated with the consumption of fruit juice, combined fruit and fruit juice and combined fruit, fruit juice and vegetables. It negatively correlated with the consumption of sweet foods and combined sweet foods and carbonated drinks. The strength of these relationships could perhaps be tested in larger cohorts with greater statistical power [23]; however, this is not the first time that citric acid has been associated with the MD. There is certainly evidence elsewhere to support a dietary link. Urine levels of citric acid were 26% higher in men and 35% higher in women in consuming a lactovegetarian diet compared with omnivorous control groups [23]. Similarly, in the PREDIMED study citric acid levels were significantly higher in urine from participants following an MD for 1 year compared with those following a low-fat diet [1]. This suggests that citric acid has a degree of specificity for the MD. There is evidence to suggest that where extra-virgin olive oil is consumed it could be a key source of dietary citric acid [1] and this specific association needs to be further explored. In contrast, there are metabolites reported in this study which were not significant in the present study, namely 3-hydroxybutyrate and cis-aconitate (these two, like citrate, were metabolites from the metabolism of carbohydrates, creatinine, amino acids (proline, N-acetylglutamine, glycine, branched-chain amino acids, and derived metabolites), lipids (oleic and suberic acids) and microbial cometabolites (phenylacetylglutamine and p-cresol)).
Other studies where adherence to MD was assessed explored the associations of an MDS with similar food groups. The alternate Mediterranean diet study found associations between MDS and fruit, vegetables, whole grains, fish and unsaturated fat. However, in this case, the highest point scale to assess adherence to MD was 8 and the highest score achieved was 4. This study found 21 identifiable metabolites associated with MD including four amino acids, one carbohydrate, two co-factors or vitamins, 11 lipids, and three xenobiotics [19].
Citric acid has previously been associated with resistant hypertension [24]. Our data showed a significant positive correlation between citric acid and both systolic and diastolic blood pressure. However, citric acid also proved to be significantly positively correlated to MDS. Correlations showed a stronger association and higher level of significance between citric acid and MDS (r = 0.280; p-value: 0.001) than those with blood pressure (r = 0.187; p-value: 0.032; r = 217; p-value: 0.012 for systolic and diastolic pressure, respectively). Systolic blood pressure did not significantly change throughout the intervention when doing an analysis of variance (ANOVA) comparing systolic blood pressure at baseline, 6 months and 12 months (p-value: 0.365.). Diastolic blood pressure showed significant differences between baseline and 12 months (p-value: 4.7 × 10−4) but not between baseline and 6 months (p-value: 0.138). The highest increase in citric acid was observed between baseline and 6 months. Citric acid showed a decrease between 6 and 12 months. Similarly, MD adherence was higher between baseline and 6 and decreased between 6 and 12 months, which strengthens a potential relationship between citric acid and MD adherence.
Pyruvic acid and citric acid are closely related metabolites in cellular metabolism and the results for pyruvic acid in the present study contrasted greatly with citric acid. Blood pyruvic acid was approximately 20% lower in individuals with a high MDS. Furthermore, pyruvic acid negatively correlated with the consumption of fruit, combined fruit and fruit juice, fish and olive spreads and positively correlated with processed meat. The negative relationship between pyruvic acid and fish consumption was the most statistically significant and strongest metabolite-food correlation observed. Furthermore, the positive correlation between pyruvic acid and processed meat was the second most statistically significant. It is worth remembering that all patients participating in the present study had a history of CVD or unstable angina, and that there was a significant difference in blood pressure between the two groups. Altered pyruvic acid has previously been associated with hypertension [25]. However, a correlation did not show any association between blood pressure and pyruvic acid in the present study. Excessive levels of pyruvic acid have also been linked with cancer [26], and even neurodegenerative diseases [27,28]. For example, pyruvic acid levels are around 20 times higher in highly invasive cancer cells compared with low invasive cancer cells [26], and pyruvic acid is two to four fold higher in diabetic rats than non-diabetic rats [29]. Pyruvic acid has not been previously linked to MD and further work should be completed to further elucidate this association. However, it is noteworthy that normal weight human subjects have significantly higher levels of both pyruvic acid and citric acid compared with severely obese patients [30]. In our dataset, patients’ weight was significantly positively associated with pyruvic acid but not citric acid (data not shown).
A number of other metabolite-food relationships were observed. Mannose significantly negatively correlated with fruit, combined fruit and fruit juice, combined fruit, fruit juice and vegetables and with whole grain cereals. We are not aware of any dietary-based studies that have been linked with mannose. However, the analysis of a large cohort of 2204 females found that after glucose, plasma mannose had the strongest metabolite association with type 2 diabetes [31].
We found that myo-inositol (the most abundant stereoisomer of inositol) negatively correlated with olive spreads and positively correlated with processed meat consumption. Neither fruit or vegetable consumption, nor red meat consumption were correlated with plasma myo-inositol levels. Initially, these observations seemed counter-intuitive given that myo-inositol occurs very abundantly in numerous plant-based foods associated with the MD, such as fruits, beans, grains and nuts [32]. However, is it well known that animal sources of inositol are more bioavailable than the phytate form, which is present in plant sources [33]. Furthermore, some of the richest, most bioavailable sources of inositol are organ meats, which may possibly explain the positive correlation observed for processed meat, but not for red meat.
Betaine (a derivative of choline) is naturally present in a variety of plant based foods such as cereals, grains and leafy vegetables and beets [34,35] and has previously been linked to MD [36]. It also may be a useful biomarker for identifying participants following the DASH (dietary approach to stop hypertension) diet [37]. Our investigation found that betaine was 15.5% higher in the “high” MD group, and it positively associated with the consumption of fruit juice and combined fruit, fruit juice and vegetables. This would indicate that diet is an important contributor. Such a link remains controversial. A cohort case study from the PREDIMED trial involving individuals with high CVD risk divided into an MD group and a control group found no significant associations between individual metabolites from the choline pathway and CVD after one year. However, the ratio betaine/choline showed an inverse association with CVD incidence although not with stroke alone [38]. Furthermore, in a cross-sectional subset of the Nutrition, Aging, and Memory in Elders cohort, betaine was also associated with lower levels of low-density lipoprotein cholesterol and triglycerides and reduced risk of diabetes mellitus [39]. However, further studies in different populations are necessary to confirm the link between betaine and CVD.
Citric acid and betaine significantly and positively correlated to MDS (r = 0.28, r = 0.19 respectively). Mannose, pyruvic acid and myo-inositol significantly and negatively correlated to MDS (r = −0.2, r = −0.02, r = −0.18 respectively). Citric acid showed to be the highest and most significant correlation with MDS, however overall these associations showed to be weak and when applying multiple comparison using the Benjamini-Hochberg approach, the FDR did not show to be less than 0.05 (Supplementary Materials Table S1). Further studies should be carried out to prove their association with MDS.
The limitations of the present study should be taken into consideration for future studies. A total of 58 patients participated in the study followed an MD during 12 months and plasma samples were collected at baseline, 6 months and 12 months. Samples were divided into two groups (low and high MDS) according to how well patients adhere to the 14-point scale MDS. Patients’ MDS was the lowest at baseline and it improved throughout the intervention; therefore, the majority of samples in the Low MDS group were from the baseline and the majority of samples in the High MDS group were collected at 6 and 12 months of the intervention. Despite the degree of dependency of the samples in both groups, this was not detrimental to the results. A parallel analysis was carried out dividing samples into timepoints: baseline, 6 and 12 months (paired samples) and results showed the same top five metabolites as the most significant ones: citric acid, betaine, mannose, myo-inositol and pyruvic acid. It is also worth noting that the concentrations of these five metabolites were not altered when applying paired or unpaired analysis.
Patients had a history of myocardial infarction or unstable angina previous to the intervention and were on a range of medication. However, drug intake diaries were revised to confirm that there was no association between outliers in the PCA model and drug intake. Patients were mostly overweight or obese; however, body mass index (BMI) did not show significant differences when comparing the low and high MDS groups. The present study focused on the discovery of biomarkers purely associated to diet and the effect of exercise was not recorded or taken into consideration. Blood pressure showed to be significant between the low and the high MDS groups, there were similar percentages of smokers in the low and high MDS groups and the percentage of males compared to females was higher in both groups. Gender did not show any significant differences in the results. Smoking status of individuals did not show any detectable influence on metabolite concentrations, although the p-value for citric acid was close to statistical significance. Further studies in larger cohorts could examine whether a relationship exists here. It should also be noted that none of the participants achieved the highest MDS on the 14-point scale and the highest scored achieved was 10.
The present study conducted untargeted metabolomics with the goal of profiling the plasma metabolome of two groups with high and low adherence to MD. 1H-NMR was used due to its high precision, reproducibility and quantitative power. This study found significant direct associations between MD and citric acid and betaine and significant inverse associations between MD and pyruvic acid, mannose and myo-inositol. These metabolites were significant following univariate analysis, but multivariate analysis corroborated these by indicating that these five metabolites were most responsible for the discrimination of the model (VIP scores of the PLS-DA model). Out of the five metabolites, citric acid and betaine have also been associated with MD on previous studies. Nevertheless, following further studies, targeted metabolomics should be carried out to confirm the validity of these five metabolites and their association with MD.

4. Materials and Methods

4.1. Patient Cohort

The Mediterranean Diet in Northern Ireland (MEDDINI) study was a pilot randomised controlled parallel group trial where 61 willing participants previously diagnosed with coronary heart disease (CHD) were recruited (samples from n = 58 participants were available for the present study) from the Cardiology Directorate, Royal Victoria Hospital, Belfast. Patients recruited were aged between 39 and 78 years and had a recent diagnosis (within four weeks prior to starting the study) of myocardial infarction (MI) or unstable angina. Patients with severe heart failure, pending on immediate inpatient coronary revascularisation, on warfarin therapy, on Omacor therapy, with cognitive impairment, those with records of extreme alcohol intake, taking multivitamin/fish oil supplements, unable to comply with the diet or to provide informed consent or those who were not expected to live longer than 6 months for any other causes were excluded from the study [13].
Plasma samples and dietary assessment measurements were collected at three-time points; baseline, 6 months and 12 months. Seven-day food diaries were used to collect food consumption data. Participants were asked to record the foods consumed over seven consecutive days, including an estimation of quantity consumed and information on preparation methods used. From seven-day food diaries, a database was created registering all food amounts eaten by all participants during the course of the intervention (baseline, 6 months and 12 months). 17 selected food groups from this database were used to correlate with shortlisted biomarkers (Supplementary Materials Table S1). Fasting blood samples were collected at the three time points to assess nutritional biomarkers. Ethical approval was obtained for this study from the Queen’s University Belfast Research Ethics Committee (ethical approval references: RGHT000049 and 15.42 for the Meddini original study and for the latest analysis respectively. Informed written consent was obtained from all participants [13].
The baseline characteristics of the sampled cohort (average BMI, blood pressure, age, sex and smoking characteristics) can be found in the Supporting information (Supplementary Materials Table S2). Mean age was 56.1 years for the low MDS group and 60.5 y for the high MDS group, and participants were on a range of medications. Mean BMI was 30.2 kg/m2 for the low MDS group and 29.8 for the high MDS group, which indicated that patients were mostly overweight/obese.

4.2. MD Scoring

Adherence to the MD was originally measured using a nine-point item scale, which has previously been validated [40] and consisted of a brief questionnaire for assessing adherence to MD and obtaining a rapid feedback. However, to include and control a wider range of food features of MD, adherence to an MD was rescored using the validated 14-point Mediterranean diet score (MDS) based on the PREDIMED score, where a score of 0 indicates lowest adherence to an MD and a score of 14 indicates highest adherence [3]. Recent MD advice was taken into consideration [41] and types and quantities of foods within the PREDIMED score were adapted to reflect the typical diet and dietary recommendations in Northern Ireland [42].

4.3. 1H-NMR Analysis of Blood

Plasma samples were collected from 58 participants (n = 137) at baseline (n = 58), 6 months (n = 44) and 12 months (n = 35). After blood collection, samples were left undisturbed at room temperature for 15 minutes for the blood to clot. Samples were centrifuged (2000× g; 10 min) and the supernatant (plasma) collected, aliquoted and stored at −80 °C. Plasma filtering was carried out using an adapted version of a previously used method [43]. Prior to filtration, centrifugal filter units (Amicon Ultracel; 0.5 mL; 3 kDa cut-off) were rinsed seven times each with H2O (0.5 mL) and centrifuged at 12,000× g for 20 min, to remove residual glycerol bound to the filter membranes. Each plasma sample (300 μL) was transferred to the centrifuge filter units. The samples were then centrifuged (13,000× g for 30 min at 4 °C) to remove macromolecules, which consist primarily of proteins and lipoproteins. Filtered plasma (285 µL) from each sample was transferred to a separate microcentrifuge tube. To each sample, D2O (35 μL) and buffer solution (30 μL; 11.7 mM DSS [disodium-2,2-dimethyl-2-silapentane-5-sulphonate], 1.75 M K2HPO4 (anhydrous), and 5.84 mM 2-chloro pyrimidine-5-carboxylic acid in H2O) was added. Samples (200 μL) were subsequently transferred to a standard Bruker 3mm NMR tube for analysis. Data collection was carried out on a 600 MHz Bruker ASCEND NMR spectrometer equipped with a five mm TCI cryoprobe using randomised order with two hundred and fifty-six transients acquired. Chemical shifts were reported in parts per million (ppm) of the operating frequency. DSS was used as the internal standard for chemical shift referencing and quantification. Using Bayesil (a web-based system that consists on a library of pure compounds specifically designed for the identification and quantification of 1D 1H-NMR metabolites) [44,45] all collected spectra were profiled using a custom library of 59 metabolites.

4.4. Statistical Analysis

To explore the metabolite associations, the 14-point MDS (Supplementary Materials Table S7) was used to divide the samples into two groups. The highest MDS achieved by any individual on the 14-point scale was 10; hence, the resulting groups were: Low Score: (0–4) and High Score (5–10) by splitting at the median (Supplementary Materials Table S3). Results obtained from these two groups were statistically analysed by performing univariate analysis (t-test with the software MetaboAnalyst (version 4.0) [46] to obtain the most significant features from identified metabolites. All data were normalised and auto-scaled (normalisation by sum and auto-scaling on Metaboanalyst) prior to univariate and multivariate analysis. False discovery rates (FDR, q-value) were calculated to account for multiple comparisons applying the Bonferroni correction on Metaboanalyst. p-value, FDR adjusted p-values (q-values) and fold-change were calculated. Principal Component Analysis (PCA), Partial Least Squares-Discriminant Analysis (PLSDA) and Orthogonal Projections to Latent Structures Discriminant Analysis (OPLS-DA) were carried out with both MetaboAnalyst (v4.0) and Simca (v14.0; Umetrics, Umea, Sweden) to observe group discrimination. PCA was also used to identify potential outliers, followed by PLS-DA to highlight significant metabolites which explained the maximum amount of variation between groups. An important issue with PLS-DA is deciding on the number of latent variables to be used to build the model. The default option offered by Metaboanalyst was used where the optimal number of latent variables was determined by Q2 (cross validated R2). Variable importance in projection (VIP) plots was used to identify the most influential metabolites responsible for the observed separation between groups. ROC curve analysis is widely used to describe the trade-off between sensitivity and specificity with regard to biomarker performance. The area under the ROC curve (AUROC) value is a robust measure used for biomarker discovery. MetaboAnalyst was used to develop ROC curves of individual metabolites (Table 1). In order to enhance this performance, pair-wise ratios of all possible metabolite concentrations were computed by Metaboanalyst showing those top-ranked ratios (Table 2). This procedure was applied for biomarker discovery rather than performance evaluation, due to the potential of over fitting. Data were tested for normality of distribution using the Shapiro Wilk test (IBM SPSS Statistics 25). All variables tested for correlations were nonparametric and therefore Spearman’s rank correlation coefficient was obtained. Correlations were deemed significant where p was ≤ 0.05. Multiple comparison correction was applied to p-values of the metabolite-food groups correlation using the Benjamini-Hochberg approach on SPSS. The effects of smoking on blood metabolites were assessed by conducting a t-test between smoker/non-smoker groups using Metaboanalyst.

5. Conclusions

In conclusion, this study reports for the first time that there is a potential association between blood levels of pyruvic acid, mannose and myo-inositol with MD consumption, and it corroborates previous associations with citric acid and betaine. It is only the third study to examine the relationship between blood metabolites levels and MD and it demonstrates the power of 1H-NMR metabolomics in dietary intervention studies. Further studies should explore the validity of these biomarkers/biomarker ratios. Given the heterogeneous nature of the MDS, the metabolites highlighted here are reasonably discriminative, and there is potential to incorporate these into currently available blood-based biomarkers panels of dietary intake.

Supplementary Materials

The following are available online at https://www.mdpi.com/2218-1989/9/10/201/s1, Table S1: Correlations of the top 5 metabolites correlations with MDS and each of the food groups assessed in the 14-item MD score, Table S2: Baseline characteristics of participants, Table S3: Assigned MDS scores listed by group, Table S4: Performance of multivariate modelling of metabolomic data, Table S5: Food group-to-Food group correlations within the study cohort, Table S6: Correlations of all 59 metabolites against each food group scored within the MDS, Table S7: 14-point criteria used to measure adherence to MD on the MEDDINI study.14 POINT MEDITERRANEAN DIET SCORE (MDS).

Author Contributions

B.D.G., J.V.W. and S.F.G. designed research and supervised the project. J.K. and A.Y. conducted research. S.E.M. and S.M. scored the food diaries of the MEDDINI patients. S.M. analysed data and wrote the manuscript with the support of B.D.G., J.V.W., S.F.G., J.K., A.Y. and S.E.M.; M.C.M. directed the dietary analysis. P.P.M. was involved in the original study design, advised on patient recruitment and helped the interpretation of the data. All authors provided critical feedback and helped shape the research, analysis and manuscript. All authors have read and approved the final manuscript

Funding

This research was funded by the Research and Development Office, Northern Ireland and the Department for Employment and Learning (DEL) Northern Ireland.

Acknowledgments

We thank the participants of the MEDDINI study and Karen Logan, who carried out the recruitment, delivered the interventions and completed all outcome assessment. A special thanks to Charlotte Neville for her assistance with governance procedures for the current analysis.

Conflicts of Interest

The authors have no conflict of interest to declare.

References

  1. Vázquez-Fresno, R.; Llorach, R.; Urpi-Sarda, M.; Lupianez-Barbero, A.; Estruch, R.; Corella, D.; Fitó, M.; Arós, F.; Ruiz-Canela, M.; Salas-Salvadó, J.; et al. Metabolomic pattern analysis after Mediterranean diet intervention in a nondiabetic population: A 1- and 3-year follow-up in the PREDIMED study. J. Proteome Res. 2015, 14, 531–540. [Google Scholar] [CrossRef] [PubMed]
  2. Mente, A.; de Koning, L.; Shannon, H.; Anand, S.S. A systematic review of the evidence supporting a causal link between dietary factors and coronary heart disease. Arch. Intern. Med. 2009, 169, 659–669. [Google Scholar] [CrossRef] [PubMed]
  3. Martínez-González, M.; García-Arellano, A.; Toledo, E.; Salas-Salvadó, J.; Buil-Cosiales, P.; Corella, D.; Covas, M.; Schröder, H.; Arós, F.; Gómez-Gracia, E.; et al. A 14-item Mediterranean diet assessment tool and obesity indexes among high-risk subjects: The PREDIMED trial. PLoS ONE 2012, 7, e43134. [Google Scholar] [CrossRef] [PubMed]
  4. Estruch, R.; Ros, E.; Salas-Salvadó, J.; Covas, M.I.; Corella, D.; Arós, F.; Gómez-Gracia, E.; Ruiz-Gutiérrez, V.; Fiol, M.; Lapetra, J.; et al. Primary Prevention of Cardiovascular Disease with a Mediterranean Diet Supplemented with Extra-Virgin Olive Oil or Nuts. N. Engl. J. Med. 2018, 378, e34. [Google Scholar] [CrossRef] [PubMed]
  5. Damasceno, N.; Sala-Vila, A.; Cofán, M.; Pérez-Heras, A.; Fitó, M.; Ruiz-Gutiérrez, V.; Martínez-González, M.; Corella, D.; Arós, F.; Estruch, R.; et al. Mediterranean diet supplemented with nuts reduces waist circumference and shifts lipoprotein subfractions to a less atherogenic pattern in subjects at high cardiovascular risk. Atherosclerosis 2013, 230, 347–353. [Google Scholar] [CrossRef] [PubMed]
  6. Grosso, G.; Mistretta, A.; Marventano, S.; Purrello, A.; Vitaglione, P.; Calabrese, G.; Drago, F.; Galvano, F. Beneficial effects of the Mediterranean diet on metabolic syndrome. Curr. Pharm. Des. 2014, 20, 5039–5044. [Google Scholar] [CrossRef] [PubMed]
  7. Grosso, G.; Stepaniak, U.; Micek, A.; Topor-Mądry, R.; Stefler, D.; Szafraniec, K.; Bobak, M.; Pająk, A. A Mediterranean-type diet is associated with better metabolic profile in urban Polish adults: Results from the HAPIEE study. Metabolism 2015, 64, 738–746. [Google Scholar] [CrossRef] [PubMed]
  8. Kesse-Guyot, E.; Ahluwalia, N.; Lassale, C.; Hercberg, S.; Fezeu, L.; Lairon, D. Adherence to Mediterranean diet reduces the risk of metabolic syndrome: A 6-year prospective study. Nutr. Metab. Cardiovasc. Dis. 2013, 23, 677–683. [Google Scholar] [CrossRef]
  9. Babio, N.; Bulló, M.; Basora, J.; Martínez-González, M.; Fernández-Ballart, J.; Márquez-Sandoval, F.; Molina, C.; Salas-Salvadó, J. Adherence to the Mediterranean diet and risk of metabolic syndrome and its components. Nutr. Metab. Cardiovasc. Dis. 2009, 19, 563–570. [Google Scholar] [CrossRef]
  10. Esko, T.; Hirschhorn, J.N.; Feldman, H.A.; Hsu, Y.H.; Deik, A.A.; Clish, C.B.; Ebbeling, C.B.; Ludwig, D.S. Metabolomic profiles as reliable biomarkers of dietary composition. Am. J. Clin. Nutr. 2017, 105, 547–554. [Google Scholar] [CrossRef]
  11. Woodside, J.V.; Draper, J.; Lloyd, A.; McKinley, M.C. Use of biomarkers to assess fruit and vegetable intake. Proc. Nutr. Soc. 2017, 76, 308–315. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Almanza-Aguilera, E.; Urpi-Sarda, M.; Llorach, R.; Vázquez-Fresno, R.; Garcia-Aloy, M.; Carmona, F.; Sanchez, A.; Madrid-Gambin, F.; Estruch, R.; Corella, D.; et al. Microbial metabolites are associated with a high adherence to a Mediterranean dietary pattern using a 1H-NMR-based untargeted metabolomics approach. J. Nutr. Biochem. 2017, 48, 36–43. [Google Scholar] [CrossRef] [PubMed]
  13. Logan, K.J.; Woodside, J.V.; Young, I.S.; McKinley, M.C.; Perkins-Porras, L.; McKeown, P.P. Adoption and maintenance of a Mediterranean diet in patients with coronary heart disease from a Northern European population: A pilot randomised trial of different methods of delivering Mediterranean diet advice. J. Hum. Nutr. Diet. 2010, 23, 30–37. [Google Scholar] [CrossRef] [PubMed]
  14. O’Gorman, A.; Gibbons, H.; Brennan, L. Metabolomics in the identification of biomarkers of dietary intake. Comput. Struct. Biotechnol. J. 2013, 4, e201301004. [Google Scholar] [CrossRef] [PubMed]
  15. Brennan, L. Progress in nuclear magnetic resonance spectroscopy. Prog. Nucl. Magn. Reson. Spectrosc. 2014, 83, 42–49. [Google Scholar] [CrossRef] [PubMed]
  16. González-Guardia, L.; Yubero-Serrano, E.; Delgado-Lista, J.; Perez-Martínez, P.; Garcia-Rios, A.; Marin, C.; Camargo, A.; Delgado-Casado, N.; Roche, H.; Perez-Jimenez, F.; et al. Effects of the Mediterranean diet supplemented with Coenzyme Q10 on metabolomic profiles in elderly men and women. J. Gerontol. A Biol. Sci. Med. Sci. 2014, 70, 78–84. [Google Scholar] [CrossRef]
  17. Garcia-Aloy, M.; Llorach, R.; Urpi-Sarda, M.; Tulipani, S.; Estruch, R.; Martínez-González, M.; Corella, D.; Fitó, M.; Ros, E.; Salas-Salvadó, J.; et al. Novel multimetabolite prediction of walnut consumption by a urinary biomarker model in a free-living population: the PREDIMED study. J. Proteome Res. 2014, 13, 3476–3483. [Google Scholar] [CrossRef] [PubMed]
  18. Bondia-Pons, I.; Martínez, J.; de la Iglesia, R.; Lopez-Legarrea, P.; Poutanen, K.; Hanhineva, K.; Zulet, M. Effects of short- and long-term Mediterranean-based dietary treatment on plasma LC-QTOF/MS metabolic profiling of subjects with metabolic syndrome features: The Metabolic Syndrome Reduction in Navarra (RESMENA) randomized controlled trial. Mol. Nutr. Food Res. 2015, 59, 711–728. [Google Scholar] [CrossRef]
  19. Playdon, M.; Moore, S.; Derkach, A.; Reedy, J.; Subar, A.; Sampson, J.; Albanaes, D.; Gu, F.; Kontto, J.; Lassale, C.; et al. Identifying biomarkers of dietary patterns by using metabolomics. Am. J. Clin. Nutr. 2017, 105, 450–465. [Google Scholar] [CrossRef]
  20. Schröder, H.; Fitó, M.; Estruch, R.; Martínez-González, M.A.; Corella, D.; Salas-Salvadó, J.; Lamuela-Raventós, R.; Ros, E.; Salaverría, I.; Fiol, M.; et al. A short screener is valid for assessing Mediterranean diet adherence among older Spanish men and women. J. Nutr. 2011, 141, 1140–1145. [Google Scholar] [CrossRef]
  21. Papadaki, A.; Johnson, L.; Toumpakari, Z.; England, C.; Rai, M.; Toms, S.; Penfold, C.; Zazpe, I.; Martínez-González, M.A.; Feder, G. Validation of the English Version of the 14-Item Mediterranean Diet Adherence Screener of the PREDIMED Study, in People at High Cardiovascular Risk in the UK. Nutrients 2018, 10, 138. [Google Scholar] [CrossRef] [PubMed]
  22. George, E.; Kucianski, T.; Mayr, H.; Moschonis, G.; Tierney, A.; Itsiopoulos, C. A Mediterranean Diet Model in Australia: Strategies for Translating the Traditional Mediterranean Diet into a Multicultural Setting. Nutrients 2018, 10, 465. [Google Scholar] [CrossRef] [PubMed]
  23. Xu, J.; Yang, S.; Cai, S.; Dong, J.; Li, X.; Chen, Z. Identification of biochemical changes in lactovegetarian urine using 1H NMR spectroscopy and pattern recognition. Anal. Bioanal. Chem. 2010, 396, 1451–1463. [Google Scholar] [CrossRef] [PubMed]
  24. Martin-Lorenzo, M.; Martinez, P.J.; Baldan-Martin, M.; Ruiz-Hurtado, G.; Prado, J.C.; Segura, J.; de la Cuesta, F.; Barderas, M.G.; Vivanco, F.; Ruilope, L.M.; et al. Citric Acid Metabolism in Resistant Hypertension Underlying Mechanisms and Metabolic Prediction of Treatment Response. Hypertension 2017, 70, 1049–1056. [Google Scholar] [CrossRef] [PubMed]
  25. Ameta, K.; Gupta, A.; Kumar, S.; Sethi, R.; Kumar, D.; Mahdi, A.A. Essential hypertension: A filtered serum based metabolomics study. Sci. Rep. 2017, 7, 2153. [Google Scholar] [CrossRef] [PubMed]
  26. Hu, S.; Wang, J.; Ji, E.; Christison, T.; Lopez, L.; Huang, Y. Targeted Metabolomic Analysis of Head and Neck Cancer Cells Using High Performance Ion Chromatography Coupled with a Q Exactive HF Mass Spectrometer. Anal. Chem. 2015, 87, 6371–6379. [Google Scholar] [CrossRef]
  27. Ahmed, S.; Santosh, W.; Kumar, S.; Thanka Christlet, H.T. Metabolic profiling of Parkinson’s disease: Evidence of biomarker from gene expression analysis and rapid neural network detection. J. Biomed. Sci. 2009, 16, 63–74. [Google Scholar] [CrossRef] [PubMed]
  28. Kori, M.; Aydın, B.; Unal, S.; Arga, K.; Kazan, D. Metabolic Biomarkers and Neurodegeneration: A Pathway Enrichment Analysis of Alzheimer’s Disease, Parkinson’s Disease, and Amyotrophic Lateral Sclerosis. OMICS 2016, 20, 645–661. [Google Scholar] [CrossRef]
  29. Zhang, S.; Nagana Gowda, G.A.; Asiago, V.; Shanaiah, N.; Barbas, C.; Raftery, D. Correlative and quantitative 1H NMR-based metabolomics reveals specific metabolic pathway disturbances in diabetic rats. Anal. Biochem. 2008, 383, 76–84. [Google Scholar] [CrossRef] [Green Version]
  30. Gralka, E.; Luchinat, C.; Tenori, L.; Ernst, B.; Thurnheer, M.; Schultes, B. Metabolomic fingerprint of severe obesity is dynamically affected by bariatric surgery in a procedure-dependent manner. Am. J. Clin. Nutr. 2015, 102, 1313–1322. [Google Scholar] [CrossRef]
  31. Menni, C.; Fauman, E.; Erte, I.; Perry, J.R.; Kastenmuller, G.; Shin, S.Y.; Petersen, A.K.; Hyde, C.; Psatha, M.; Ward, K.J.; et al. Biomarkers for type 2 diabetes and impaired fasting glucose using a nontargeted metabolomics approach. Diabetes 2013, 62, 4270–4276. [Google Scholar] [CrossRef] [PubMed]
  32. Clements, R.; Darnell, B. Myo-inositol content of common foods: Development of a high-myo-inositol diet. Am. J. Clin. Nutr. 1980, 33, 1954–1967. [Google Scholar] [CrossRef] [PubMed]
  33. Muscogiuri, G.; Palomba, S.; Laganà, A.; Orio, F. Current insights into inositol isoforms, Mediterranean and Ketogenic diets for polycystic ovary syndrome: From bench to bedside. Curr. Pharm. Des. 2016, 22, 5554–5557. [Google Scholar] [CrossRef] [PubMed]
  34. Zeisel, S.; Mar, M.; Howe, J.; Holden, J. Concentrations of choline-containing compounds and betaine in common foods. J. Nutr. 2003, 133, 1302–1307. [Google Scholar] [CrossRef] [PubMed]
  35. Graham, S.F.; Hollis, J.H.; Migaud, M.; Browne, R.A. Analysis of Betaine and Choline Contents of Aleurone, Bran, and Flour Fractions of Wheat (Triticum aestivum L.) Using 1H Nuclear Magnetic Resonance (NMR) Spectroscopy. J. Agric. Food Chem. 2009, 57, 1948–1951. [Google Scholar] [CrossRef] [PubMed]
  36. Detopoulou, P.; Panagiotakos, D.; Antonopoulou, S.; Pitsavos, C.; Stefanadis, C. Dietary choline and betaine intakes in relation to concentrations of inflammatory markers in healthy adults: The ATTICA study. Am. J. Clin. Nutr. 2008, 87, 424–430. [Google Scholar] [CrossRef] [PubMed]
  37. Rebholz, C.M.; Lichtenstein, A.H.; Zheng, Z.; Appel, L.J.; Coresh, J. Serum untargeted metabolomic profile of the Dietary Approaches to Stop Hypertension (DASH) dietary pattern. Am. J. Clin. Nutr. 2018, 108, 243–255. [Google Scholar] [CrossRef]
  38. Guasch-Ferré, M.; Hu, F.; Ruiz-Canela, M.; Bulló, M.; Toledo, E.; Wang, D.D.; Corella, D.; Gómez-Gracia, E.; Fiol, M.; Estruch, R.; et al. Plasma Metabolites from choline pathway and risk of cardiovascular disease in the PREDIMED (Prevention with Mediterranean Diet) study. J. Am. Heart Assoc. 2017, 6, e006524. [Google Scholar] [CrossRef]
  39. Roe, A.J.; Zhang, S.; Bhadelia, R.A.; Johnson, E.J.; Lichtenstein, A.H.; Rogers, G.T.; Rosenberg, I.H.; Smith, C.E.; Zeisel, S.H.; Scott, T.M. Choline and its metabolites are differently associated with cardiometabolic risk factors, history of cardiovascular disease, and MRI-documented cerebrovascular disease in older adults. Am. J. Clin. Nutr. 2017, 105, 1283–1290. [Google Scholar] [CrossRef]
  40. Martínez-González, M.; Fernández-Jarne, E.; Serrano-Martínez, M.; Wright, M.; Gomez-Gracia, E. Development of a short dietary intake questionnaire for the quantitative estimation of adherence to a cardioprotective Mediterranean diet. Eur. J. Clin. Nutr. 2004, 58, 1550–1552. [Google Scholar] [CrossRef] [Green Version]
  41. Bach-Faig, A.; Berry, E.; Lairon, D.; Reguant, J.; Trichopoulou, A.; Dernini, S.; Medina, F.; Battino, M.; Belahsen, R.; Miranda, G.; et al. Mediterranean diet pyramid today. Science and cultural updates. Public Health Nutr. 2011, 14, 2274–2284. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. McEvoy, C.T.; Moore, S.E.; Appleton, K.M.; Cupples, M.E.; Erwin, C.M.; Hunter, S.J.; Kee, F.; McCance, D.; Patterson, C.C.; Young, I.S.; et al. Trial to Encourage Adoption and Maintenance of a Mediterranean Diet (TEAM-MED): Protocol for a Randomised Feasibility Trial of a Peer Support Intervention for Dietary Behaviour Change in Adults at High Cardiovascular Disease Risk. Int. J. Environ. Res. Public Health 2018, 15, 1130. [Google Scholar] [CrossRef] [PubMed]
  43. Mercier, P.; Lewis, M.; Chang, D.; Baker, D.; Wishart, D. Towards automatic metabolomic profiling of high-resolution one-dimensional proton NMR spectra. J. Biomol. NMR 2011, 49, 307–323. [Google Scholar] [CrossRef] [PubMed]
  44. Ravanbakhsh, S.; Liu, P.; Bjordahl, T.C.; Mandal, R.; Grant, J.R.; Wilson, M.; Eisner, R.; Sinelnikov, I.; Hu, X.; Luchinat, C.; et al. Accurate, fully-automated NMR spectral profiling for metabolomics. PLoS ONE 2015, 10, e0124219. [Google Scholar] [CrossRef] [PubMed]
  45. Graham, S.F.; Holscher, C.; Green, B.D. Metabolic signatures of human Alzheimer’s disease (AD): 1H NMR analysis of the polar metabolome of post-mortem brain tissue. Metabolomics 2013, 10, 744–753. [Google Scholar] [CrossRef]
  46. Chong, J.; Soufan, O.; Li, C.; Caraus, I.; Li, S.; Bourque, G.; Wishart, D.S.; Xia, J. MetaboAnalyst 4.0: Towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 2018, 46, W486–W494. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Multivariate statistical modelling of 1H-Nuclear Magnetic Resonance 1H-NMR metabolomic data. Plots A–C show group separation achieved by principal component analysis (PCA), orthogonal partial least squares discriminant analysis (OPLS-DA) and partial least squares discriminant analysis (PLS-DA). Red circles (1) represent patients with high MDS and green circles (2) represent individuals with low MDS. D: Is the resulting variable importance in project (VIP) plot indicating the 15 most influential metabolites responsible for the observed separation in the PLS-DA model.
Figure 1. Multivariate statistical modelling of 1H-Nuclear Magnetic Resonance 1H-NMR metabolomic data. Plots A–C show group separation achieved by principal component analysis (PCA), orthogonal partial least squares discriminant analysis (OPLS-DA) and partial least squares discriminant analysis (PLS-DA). Red circles (1) represent patients with high MDS and green circles (2) represent individuals with low MDS. D: Is the resulting variable importance in project (VIP) plot indicating the 15 most influential metabolites responsible for the observed separation in the PLS-DA model.
Metabolites 09 00201 g001
Figure 2. Blood levels of citric acid and pyruvic acid and MDS. (A) Citric acid levels positively correlated with MDS, whereas Pyruvic acid levels negatively correlated. (B) A paired metabolite ratio of citric acid and pyruvic acid achieved the greatest area under the receiver operating characteristic (ROC) curve (AUC = 0.74) value which was the best performing biomarker for MDS. Optimal cutoff was represented on the curve with a red dot and with a red horizontal line on the box plot. The box-plot showed that the overall distribution profiles for low MDS (0–4) and high MDS (5–10) were broadly similar but the mean value was 53% higher in the high MDS group. Y-axis represents concentrations (µM). Data were median centred. The mean concentration of each group was indicated with a red diamond.
Figure 2. Blood levels of citric acid and pyruvic acid and MDS. (A) Citric acid levels positively correlated with MDS, whereas Pyruvic acid levels negatively correlated. (B) A paired metabolite ratio of citric acid and pyruvic acid achieved the greatest area under the receiver operating characteristic (ROC) curve (AUC = 0.74) value which was the best performing biomarker for MDS. Optimal cutoff was represented on the curve with a red dot and with a red horizontal line on the box plot. The box-plot showed that the overall distribution profiles for low MDS (0–4) and high MDS (5–10) were broadly similar but the mean value was 53% higher in the high MDS group. Y-axis represents concentrations (µM). Data were median centred. The mean concentration of each group was indicated with a red diamond.
Metabolites 09 00201 g002
Table 1. Metabolite levels differing according to low or high Mediterranean diet score MDS.
Table 1. Metabolite levels differing according to low or high Mediterranean diet score MDS.
NameMean [µM]
Low MDS
(n = 64)
SD. Low MDSMean [µM] High MDS
(n = 73)
SD. High MDSp-Valueq-Value (FDR)AUC (ROC)↑/↓% Change
1citric acid*76.6230.5092.3632.585.99 × 10−4*0.030.6720.5
2pyruvic acid*63.5034.8050.1934.360.005*0.160.64−20.9
3betaine*40.3717.8346.6418.420.01*0.280.6215.5
4mannose*55.0844.1439.2536.580.02*0.320.61−28.7
5acetic acid13.067.5414.946.550.03*0.340.6114.39
6myo-inositol*104.6774.4179.2247.250.03*0.340.60−24.3
7tyrosine4.582.745.482.860.060.520.5919.7
8glycerol1521.19817.491304.62875.720.080.540.59−14.23
9dimethylglycine3.011.752.371.560.090.540.58−21.3
10malonate8.357.038.994.350.090.540.587.66
11pyroglutamic acid87.2328.5795.3529.650.100.540.589.3
Table 2. Ratios of paired metabolites improve biomarker utility.
Table 2. Ratios of paired metabolites improve biomarker utility.
NameAUC (ROC)↑/↓% Change
1citric acid/pyruvic acid0.7453.3
2citric acid/l-lactic acid0.7329.2
3citric acid/aspartate0.7238.3
4citric acid/l-phenylalanine 0.71 27.1
5betaine/pyruvic acid 0.70 47.6

Share and Cite

MDPI and ACS Style

Macias, S.; Kirma, J.; Yilmaz, A.; Moore, S.E.; McKinley, M.C.; McKeown, P.P.; Woodside, J.V.; Graham, S.F.; Green, B.D. Application of 1H-NMR Metabolomics for the Discovery of Blood Plasma Biomarkers of a Mediterranean Diet. Metabolites 2019, 9, 201. https://doi.org/10.3390/metabo9100201

AMA Style

Macias S, Kirma J, Yilmaz A, Moore SE, McKinley MC, McKeown PP, Woodside JV, Graham SF, Green BD. Application of 1H-NMR Metabolomics for the Discovery of Blood Plasma Biomarkers of a Mediterranean Diet. Metabolites. 2019; 9(10):201. https://doi.org/10.3390/metabo9100201

Chicago/Turabian Style

Macias, Shirin, Joseph Kirma, Ali Yilmaz, Sarah E. Moore, Michelle C. McKinley, Pascal P. McKeown, Jayne V. Woodside, Stewart F. Graham, and Brian D. Green. 2019. "Application of 1H-NMR Metabolomics for the Discovery of Blood Plasma Biomarkers of a Mediterranean Diet" Metabolites 9, no. 10: 201. https://doi.org/10.3390/metabo9100201

APA Style

Macias, S., Kirma, J., Yilmaz, A., Moore, S. E., McKinley, M. C., McKeown, P. P., Woodside, J. V., Graham, S. F., & Green, B. D. (2019). Application of 1H-NMR Metabolomics for the Discovery of Blood Plasma Biomarkers of a Mediterranean Diet. Metabolites, 9(10), 201. https://doi.org/10.3390/metabo9100201

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop