The present study acquired food diary data and concomitant blood samples from a Mediterranean diet intervention study in a Northern European population. Non-targeted metabolomic profiling was performed on plasma and this was combined with the calculated MDS to shortlist 73 features significantly differing between ‘low’ and ‘high’ consumers of a Mediterranean diet. Ultimately, this led to the putative identification of 7 high performing metabolite biomarkers (ROCAUC ≤ 0.79), which were also highly influential in multivariate modelling, and strongly correlated with MDS. Correlation of these biomarkers against each of the food groups involved in the calculation of MDS provided potential information on the possible origin of some of these. Using logistic regression analysis, we developed a model which accurately distinguished between the two dietary groups. We developed a predictive algorithm using the acquired data which had an AUC (95% CI) = 0.83 (0.76–0.89) with corresponding sensitivity and specificity equal to 0.79 (0.79–0.89) and 0.72 (0.62–0.83). To our knowledge, no logistic regression model currently exists for human plasma that distinguishes low and high MDS with such a high degree of accuracy.
The results appear impressive, given that adherence in the MEDDINI intervention study (from which samples and data were obtained) would be described as sub-maximal for MD (Maximum MDS achieved was 10 from the 14-point scale. Mean MDS in the High MDS group was 6.68). Participants were not supplemented with any foods, and as such, did not eat identical foods. We should be mindful of the potential complexities here, particularly that the participants had a history of CVD and were taking prescribed medications. However, a number of the metabolites uncovered as biomarkers have obvious connections with food consumption (especially those with ROCAUC values > 0.7). Top performing putative biomarkers included PC (P-18:1(9Z)/16:0), EPA (an n-3 long chain polyunsaturated fatty acid), and a closely related metabolite LysoPC(20:5(5Z,8Z,11Z,14Z,17Z)/0:0). The strong correlation of each of these, not only with MDS, but also with fish consumption underscores dietary relevance. MEDDINI participants did not take any fish oil supplements, which rules this out as a possible source. EPA has previously been reported as a validated biomarker of fish intake (20), but the lysophospholipid metabolite performed substantially better than the free fatty acid (spearman r = 0.524 vs 0.426). This novel biomarker ought to be closely examined in future dietary biomarker studies, but it is encouraging that one other study is supportive of it as a biomarker of MD (12). Unfortunately, due to time, budget and instrument availability it was not possible to check metabolite identifications against analytical standards. However, the intensity values for EPA and LysoPC(20:5(5Z,8Z,11Z,14Z,17Z)/0:0) were correlated against quantitated values acquired from a previously published investigation (13) which found strong associations (r = 0.595; p = 6.9E-14 and r = 0.589; p = 1.3E-13, respectively) thus providing reasonable confidence in their identity. Fish oil supplementation studies in human volunteers indicate that LysoPC (20:5(5Z,8Z,11Z,14Z,17Z)/0:0) is associated with EPA intake (21). It is thought that Lyso PC(20:5(5Z,8Z,11Z,14Z,17Z)/0:0) is more bioavailable than the free fatty acid, as gavage studies in mice increased the levels of EPA in the brain by > 100-fold(22).
Despite the potential utility of eicosapentaenoic acid metabolites, ultimately, none of these were incorporated into the optimised logistic regression model predicting MDS group. The three metabolites incorporated were MG(0:0/16:1(9Z)/0:0), PTX2SA, and PC(P-18:1(9Z)/16:0). The combined correlation of the 3 metabolites included in the model and MDS was (r = 0.64; p = 5.6E-17). MG(0:0/16:1(9Z)/0:0) was one of two monoacylglyceride biomarkers identified (the other being MG(0:0/20:3(5Z,8Z,11Z)/0:0). The respective fatty acid components of these MGs have potentially strong dietary relevance, and both markers significantly decreased with increasing MDS. Both of these metabolites were correlated with previous targeted measured blood triglycerides showing a highly significant association (r = 0.511; p = 1.2E-9) and (r = 0.348; p = 6.8E-5), respectively.
MG(0:0/20:3(5Z,8Z,11Z)/0:0) is the monoacylglyceride metabolite of mead acid (20:3(5Z,8Z,11Z). The levels of mead acid in plasma are known to be a marker for overall essential fatty acid (EFA) status (23, 24). Mead acid is not itself considered an essential fatty acid, but in the absence of adequate essential fatty acids in human tissues, the fatty acid MG(0:0/20:3(5Z,8Z,11Z)/0:0) is metabolised to mead acid (24–27). Given that MG(0:0/20:3(5Z,8Z,11Z)/0:0) is significantly negatively correlated with MDS potentially indicates that some individuals in the ‘low’ MDS group, particularly those with very high MG(0:0/20:3(5Z,8Z,11Z)/0:0) levels, may exhibit deficiency in essential fatty acids (28, 29). Adherence to MD has proved an enhancement in essential fatty acid levels (30, 31). This is consistent with our findings which showed a significant difference in EFA concentrations between low and high MD adherence groups (p = 0.001) and an 11.10% increase between low and high MD adherence.
Alternatively, it is also possible that MG (0:0/20:3(5Z,8Z,11Z)/0:0) is derived from the diet. For instance, mead acid is present in very high levels in animal cartilage, and it is noteworthy that this metabolite is strongly correlated with processed meat intake. The other MG biomarker identified is a metabolite of trans-palmitoleic acid (16:1(9Z)) which has previously been established as a marker of full fat dairy intake (32). Unfortunately, the correlation between this ion and dairy intake was not statistically significant in our study (p = 0.12).
Another novel high performing MDS biomarker which correlated strongly with fish consumption is a xenobiotic compound called Pectenotoxin-2 seco acid (PTX2SA). PTX2SA correlated closely with both MDS and fish intake. Produced by toxic dinoflagellates, pectenotoxins accumulate in shellfish, and humans can potentially be exposed to these through shellfish consumption. PTX2 is one of the family of pectenotoxin compounds, which are polyether macrolide toxins responsible for diarrheic shellfish poisoning (33). Reassuringly however, PTX2sa is a non-toxic metabolite of PTX2. Injection of mice with doses as high as 5mg/kg does not cause toxicity, and oral administration of PTX2sa is likely to be even less toxic (34). Intriguingly, PTX2sa and its epimer 7-epi-PTX2sa have previously been detected in Irish waters (35) however, we cannot find any evidence that PTX2sa has been detected in humans before. PTX2 is highly lipophilic and may not be released and absorbed during human digestion of shellfish. It also appears to be quite labile and, if it were to be liberated during digestion, acidity in the stomach would rapidly lead to metabolism to its non-toxic seco acid. PTX2sa has been detected in various marine samples, with mollusks and plankton being the most abundant sources (36). We are not aware that PTX2 or PTX2sa has been detected in fish specimens before, however, it seems likely given that trace amounts will occur, given that PTX2 has been measured at ≤ 8ng/l in seawater, ≤ 10 ng in suspended particular matter and ≤ 2ng/g in marine sediment (37). It is entirely plausible that this compound could originate from fish intake. Mollusk/shellfish consumption among MEDDINI participants during the surveyed period was, however, extremely rare, and the calculation of MDS was based on fish intake, and not shellfish intake. Alternatively, given the fact that PTX2/PTX2sa are highly lipophilic they may persist long after absorption and may reflect shellfish intake outside of the food diary data collection period.
The third metabolite to be incorporated in the model was a phosphatidylcholine identified PC(16:0/18:1(11Z)), comprised of palmitic acid and vaccenic acid. This metabolite was the best performing individual biomarker overall (AUCROC = 0.79) and it strongly positively correlated with MDS (r = 0.495; p = 1.04E-9). It is difficult to pinpoint the dietary origin of PC(16:0/18:1(11Z)) as it correlated with a number of food types, but the strongest association was with fish intake (r = 0.492, p-value = 1.3E-9)), which was almost equal in strength to its association with MDS.
One other noteworthy MDS biomarker detected was Xi-8-Hydroxyhexadecanedioic acid (Xi-8-HHDDA), also known as 8-hydroxyhexadecane dioic acid. Xi-8-HHDDA is a long-chain fatty acid, and here it was significantly associated with fruit, fruit juice and vegetable intake, and it also correlated with blood plasma levels of vitamin C (r = 0.30; p = 8E-5). It seems likely to originate from dietary plant intake, given that it is a cutin constituent of fruits and vegetables (38, 39). Its presence has been reported in fruit and tomatoes, which makes it a potential biomarker for the consumption of these food products (40, 41). It has also been identified as one of the major constituents of sweet cherries (42). Other dioic acids have been reported in other fruits, for example, 10,16-dihydroxyhexadecanoic was identified as major components of the of the cuticle of different apple varieties (43). To the best of our knowledge, this is the first time that the metabolite xi-8-Hydroxyhexadecanedioic has been identified as a potential biomarker of fruits and vegetable intake.
In conclusion, the present study is only the third to apply untargeted LC-MS metabolomics to a MD study and it provides a clear indication that this approach can be effective in a Northern European population with sub-maximal MD adherence. The findings further advance the ongoing search for a biomarker panel to determine adherence to a MD diet. Specifically, we propose a logistic regression model for accurately distinguishing low or high MDS, which will require careful validation using targeted and quantitative methods in other MD cohorts. There is clear evidence that the shortlisted metabolite biomarkers have statistically significant dietary associations (five for fish intake, one for fruit and vegetable intake, and two for processed meat), thus making the findings biologically plausible, and worthy of further investigation.