Skip to main content

Olive oil consumption, plasma metabolites, and risk of type 2 diabetes and cardiovascular disease

Abstract

Background

Olive oil consumption has been inversely associated with the risk of type 2 diabetes (T2D) and cardiovascular disease (CVD). However, the impact of olive oil consumption on plasma metabolites remains poorly understood. This study aims to identify plasma metabolites related to total and specific types of olive oil consumption, and to assess the prospective associations of the identified multi-metabolite profiles with the risk of T2D and CVD.

Methods

The discovery population included 1837 participants at high cardiovascular risk from the PREvención con DIeta MEDiterránea (PREDIMED) trial with available metabolomics data at baseline. Olive oil consumption was determined through food-frequency questionnaires (FFQ) and adjusted for total energy. A total of 1522 participants also had available metabolomics data at year 1 and were used as the internal validation sample. Plasma metabolomics analyses were performed using LC–MS. Cross-sectional associations between 385 known candidate metabolites and olive oil consumption were assessed using elastic net regression analysis. A 10-cross-validation (CV) procedure was used, and Pearson correlation coefficients were assessed between metabolite-weighted models and FFQ-derived olive oil consumption in each pair of training–validation data sets within the discovery sample. We further estimated the prospective associations of the identified plasma multi-metabolite profile with incident T2D and CVD using multivariable Cox regression models.

Results

We identified a metabolomic signature for the consumption of total olive oil (with 74 metabolites), VOO (with 78 metabolites), and COO (with 17 metabolites), including several lipids, acylcarnitines, and amino acids. 10-CV Pearson correlation coefficients between total olive oil consumption derived from FFQs and the multi-metabolite profile were 0.40 (95% CI 0.37, 0.44) and 0.27 (95% CI 0.22, 0.31) for the discovery and validation sample, respectively. We identified several overlapping and distinct metabolites according to the type of olive oil consumed. The baseline metabolite profiles of total and extra virgin olive oil were inversely associated with CVD incidence (HR per 1SD: 0.79; 95% CI 0.67, 0.92 for total olive oil and 0.70; 0.59, 0.83 for extra virgin olive oil) after adjustment for confounders. However, no significant associations were observed between these metabolite profiles and T2D incidence.

Conclusions

This study reveals a panel of plasma metabolites linked to the consumption of total and specific types of olive oil. The metabolite profiles of total olive oil consumption and extra virgin olive oil were associated with a decreased risk of incident CVD in a high cardiovascular-risk Mediterranean population, though no associations were observed with T2D incidence.

Trial registration: The PREDIMED trial was registered at ISRCTN (http://www.isrctn.com/, ISRCTN35739639).

Introduction

Olive oil has been traditionally used as the main culinary and dressing fat in Mediterranean regions and is gaining global popularity due to its remarkable nutritional profile and health benefits. It has been proposed as one of the key components of the Mediterranean Diet (MedDiet), which makes it cardio-protective [1]. The best nutritional and organoleptic quality olive oil, extra-virgin olive oil, and virgin olive oil (VOO) varieties are obtained by mechanically pressing olives. They contain an exceptional matrix of lipids rich in monounsaturated fatty acids (mainly oleic acid) and high amounts of bioactive compounds, including polyphenols (hydroxytyrosol and oleuropein), lipid derivatives (squalene, tocopherols), and vitamin E, and have a richer taste, color, and aroma than other common varieties. Because of its processing, the refined or common variety of olive oil (COO) has a similar lipid profile but a lower content of phytochemicals [2]. Evidence from epidemiological studies and clinical trials support that olive oil consumption is associated with a lower risk of type 2 diabetes (T2D), cardiovascular diseases (CVD), and mortality [3,4,5], and has demonstrated protective effects on lipid metabolism, inflammation, endothelial function, and oxidative stress [6,7,8]. Despite the mounting evidence supporting olive oil’s beneficial impact on human health, the lack of detailed differentiation between olive oil types in existing studies is a significant limitation. This oversight hinders a comprehensive understanding of the distinctive health properties offered by VOO.

Despite the wide range of health benefits attributed to olive oil consumption, the biological mechanisms underlying these salutary effects have not been well defined. Nutritional metabolomics, an evolving approach, holds great promise in enhancing our understanding of the biological effects of nutritional factors and may also help to identify potential novel biomarkers related to dietary intake or the subsequent metabolic response to this intake. Along these lines, the distinct phytochemical concentrations found in olive oil varieties may confer specific health advantages [7]. For example, VOO consumption has been associated with better HDL functionality, indicating a potential mediating role between metabolic pathways associated with cardiovascular health and VOO consumption [9, 10]. We hypothesized that the amount of total olive oil consumed, and specifically if the consumption is mainly from VOO or common olive oil, is associated with unique metabolite profiles and that these multi-metabolite profiles are related to a lower incidence of T2D and CVD.

In the current study, leveraging dietary and metabolomics data in the PREvención con DIeta MEDiterránea (PREDIMED) study, we used an agnostic machine learning approach to identify plasma metabolite profiles associated with total olive oil, VOO, and COO. We then evaluated whether the identified multi-metabolite profiles are associated with T2D and CVD incidence risk independently of known risk factors and diet.

Methods

Study population

The present analysis was conducted in the context of the PREDIMED study. This study was a Spanish multicenter randomized controlled nutritional intervention trial conducted between 2003 and 2010. The main objective was to examine the effect of the MedDiet on the primary prevention of CVD in a population with several risk factors for CVD. The full protocol of the PREDIMED study has been previously published and can be found on the study website (http://www.predimed.es/) [11, 12]. Study participants provided written informed consent and all the study centers have approved the protocol by their Institutional Review Boards. The PREDIMED trial is registered at ISRCTN (http://www.isrctn.com/ ISRCTN35739639).

The participants involved in this analysis are from three nested study samples: the first study for T2D (the PREDIMED-T2D study), the second study for CVD (the PREDIMED-CVD study), and a third subset of PREDIMED participants who completed an oral glucose tolerance test (OGTT) at baseline. The first case-cohort study consisted of 251 participants with incident T2D cases and 694 participants without T2D at baseline (overlapping participants n = 53 between cases and cohort participants) [13, 14]. The second case-cohort study consisted of 229 participants with incident CVD cases (a composite of myocardial infarction, stroke, and CVD mortality) and 788 participants without CVD at baseline (overlapping n = 37 between cases and cohort participants). The third study consisted of 132 participants. More information about these studies is available elsewhere [15, 16].

Participants who had complete metabolomics data and nutritional data at baseline (from semi-quantitative food frequency questionnaires (FFQs)) and were not duplicated (overlapping participants between databases from different metabolomics sub-studies, n = 122) were selected (n = 1,882). Additionally, participants with missing values in FFQs at baseline (n = 11), a daily energy intake lower than 500 kcal for women or 800 kcal for men and higher than 3500 kcal for women and 4000 kcal for men (n = 30) [17], or with missing values in ≥ 20% metabolites (n = 4) were excluded. Therefore, this analysis included 1,837 participants at baseline (634 participants allocated to the MedDiet supplemented with VOO group, 630 to the MedDiet supplemented with nuts, and 573 participants to the control group) (Additional file 1: Figure S1).

Further, an internal validation in the same population was conducted using dietary and metabolomics data from the 1-year visit in 1522 study participants (Additional file 1: Figure S1).

Dietary assessment

Dietary data were obtained using a validated 137-item semi-quantitative FFQ that trained dietitians collected from the participants at baseline and 1-year visits in face-to-face interviews [18]. Food, nutrients, and energy intake were estimated using Spanish food composition tables [19, 20]. Total olive oil, VOO, and COO consumption were derived from these FFQs. Total olive oil consumption was considered the sum of VOO and COO.

Anthropometric measurements and other covariates

At baseline and 1-year visits, blood pressure (in triplicate) and anthropometric measurements such as weight, height, and waist circumference were measured according to the study protocol by trained staff. Additionally, physical activity was assessed with the validated Spanish version of the Minnesota Leisure-Time Physical Activity questionnaire [21] and other information about lifestyle, medical conditions, or medication use was also collected.

Metabolite profiling

At baseline and 1-year visits, overnight fasting plasma EDTA samples (> 8 h) were collected, processed, and stored in -80◦C freezers at each recruiting center. Before metabolomics assays, case-cohort participant samples were randomized in pairs (baseline plus 1-year visit) and sent to the Broad Institute of Harvard University and the Massachusetts Institute of Technology for analysis. Metabolic profiling of the plasma samples was performed using high-throughput liquid chromatography-tandem mass spectrometry (LC–MS) techniques [22]. After quality filtration and standardization, 400 known metabolites were quantified, of which 19 metabolites were removed from the analyses (3 metabolites that were considered as internal standards (1,2-didodecanoyl-sn-glycerol-3-phosphocholine, valine-d8, and phenylalanine-d8), 7 metabolites that were drugs (acetaminophen, metronidazole, metformin, valsartan, warfarin, verapamil, atenolol) and 9 metabolites due to > 20% of missing values). The analyses were conducted with 381 known metabolites.

To quantitatively profile polar metabolites and plasma lipids, LC–MS was used as previously described [23,24,25]. Amino acids (AAs) and other polar metabolites were profiled with a Nexera X2 U-HPLC (Shimadzu Corp., Marlborough, MA) coupled to a Q-Exactive mass spectrometer (ThermoFisher Scientific, Waltham, MA). Metabolites were extracted from 10 μL plasma and 90 μL of acetonitrile/methanol/formic acid (74.9:24.9:0.2 vol:vol:vol) that contained stable isotope-labeled internal standards [valine-d8 (Sigma-Aldrich) and phenylalanine-d8 (Cambridge Isotope Laboratories)]. After centrifuging at 9,000 × g for 10 min at 4◦C, the samples supernatants were injected directly onto a 150 × 2-mm, 3-μm Atlantis HILIC column (Waters). The column was eluted isocratically at a flow rate of 250 L/min with 5% mobile phase A (10 mmol ammonium formate/L and 0.1% formic acid in water) for 0.5 min followed by a linear gradient to 40% mobile phase B (acetonitrile with 0.1% formic acid) over 10 min. MS analyses were carried out using electrospray ionization in the positive-ion mode. Full-scan spectra were acquired over 70–800 m/z. Fatty acids and other lipids were also profiled using a Nexera X2 U-HPLC (Shimadzu Corp., Marlborough, MA) coupled to an Exactive Plus orbitrap MS (Thermo Fisher Scientific) and were extracted from 10 μL plasma using 190 μL of isopropanol containing 1,2-didodecanoyl-sn-glycerol-3-phosphocholine (Avanti Polar Lipids) as an internal standard. The lipid extraction (2 μL) was injected into a 100 × 2.1-mm, 1.7-μm ACQUITY BEH C8 column (Waters). The column was eluted isocratically with 80% mobile-phase A of (95:5:0.1 vol:vol:vol 10 mM ammonium acetate/methanol/formic acid) for 1 min followed by a linear gradient to 80% mobile-phase B (99.9:0.1 vol:vol methanol/formic acid) over 2 min, a linear gradient to 100% mobile-phase B over 7 min, and then 3 min at 100% mobile-phase B. For the AAs, MS analyses were carried out using electrospray ionization in the positive-ion mode using full-scan analysis over 200–1100 m/z. Raw data were processed using Trace Finder version 3.1 and 3.3 (Thermo Fisher Scientific) and Progenesis QI (Nonlinear Dynamics). Polar metabolite identities were confirmed using authentic reference standards. Lipids were identified using the head group, total acyl carbon numbers, and total acyl double bond content. Pairs of pooled plasma reference samples were analyzed in intervals of 20 participant samples to assess data quality and to facilitate data standardization across the analytical queue and sample batches. One sample of each pair of the pooled references functioned as a passive quality control to assess the analytical measurement reproducibility of each metabolite. The other pooled sample was used to standardize using a “nearest neighbor” approach, i.e., standardized values were calculated using the ratio of the value in each sample over the nearest pooled plasma reference multiplied by the median value measured across the pooled references.

Statistical analysis

Baseline characteristics were presented as means and standard deviations (SD) when variables were quantitative and as percentages (n) when variables were categorical. The metabolites associated with total olive oil, VOO, and COO consumption were selected using plasma baseline metabolomics data (i.e., discovery sample). Plasma metabolomics data at 1 year were used as validation samples. Total olive oil, VOO, and COO consumption at baseline and 1 year were adjusted for total energy intake using the residual method [17].

The statistical quality controls used with the metabolomics data were as follows: When metabolites presented missing values of less than 20% (i.e., not detectable/quantifiable concentrations or not present metabolite), they were imputed using a random forest approach ("missForest" function from the "missForest" R package) as previous publications have recommended [26,27,28]. Metabolites were normalized and scaled using Blom's rank-based inverse normal transformation [29].

Linear regression models were used to assess associations between plasma baseline metabolites and total olive oil, VOO, and COO consumption. Models were adjusted for recruiting center, age, sex, smoking status (former smoker, never smoker, smoker), BMI, education, and physical activity (METs-min/day). Multiple testing correction was performed using the Benjamini–Hochberg procedure (FDR) and reported findings with FDR p-Adjusted value < 0.05.

To determine the metabolomic profile associated with each exposure (total olive oil, VOO, and COO) using an agnostic metabolomics approach including knowing metabolites, Gaussian linear regression models were used with the elastic net penalty (ENR) (“caret” v 6.0–84 and “glmnet” R Package). A tenfold cross-validation (CV) approach was performed with the discovery population (PREDIMED baseline data). First, the sample was split into 10 training-validation sets (90–10%, respectively) and we performed a tenfold CV to find the optimal value of the tuning parameter (λ) that results in a mean squared error within 1 SD of the minimum (minMSE). Additionally, the α parameter was evaluated from 0 (i.e., a Ridge regression) to 1 (i.e., a Lasso regression) in 0.1 increments to test the best parameters for these analyses. The best-predicting accuracy in the validation sets was obtained with α = 0.4 for total olive oil, 0.2 for VOO, and 0.3 for COO models. After evaluating the tuning parameters, a tenfold CV was performed again with the discovery population, the coefficients from each tenfold CV iteration were extracted and were constructed weighted models using the regression coefficients of the selected metabolites from each ENR in the training set to validate with its validation data set pair. Pearson correlation coefficients were determined between total olive oil, VOO, and COO consumption and the metabolomics profile in the pair of training validation data sets in both the discovery and internal validation samples (i.e., PREDIMED baseline data and PREDIMED 1-year data). For reproducibility, regression coefficients were reported using 10 iterations of the tenfold CV elastic net regression in the entire data set. These analyses were based on consistency among CV runs; therefore, any P-value was derived.

The associations between the identified metabolite profiles of total olive oil, VOO, and COO consumption (1 SD) and T2D risk (245 events at baseline and 161 incident events at 1 year) within the T2D nested case-cohort study and CVD risk (222 events at baseline and 159 incident events at 1 year) within the CVD nested case-cohort study were run with Cox regressions with Barlow weights and robust variance estimator. Four multivariate models were assessed. For the baseline analysis, the first model was adjusted for age (years), sex, and propensity scores, as previously described [12], and was stratified by intervention group and recruitment center. The second model was additionally adjusted for BMI (kg/m2), smoking status (current, former, or never), alcohol consumption in g/day (and adding a quadratic term), educational level (primary, secondary, or college), physical activity (METs/min/day), family history of CVD (yes/no), the baseline prevalence of dyslipidemia and hypertension, and dyslipidemia and hypertension medication. The third model included all covariables from model 2 in addition to the consumption of vegetables, fruits, cereals, nuts, eggs, legumes, fish, meat, and dairy products (all g/day). The fourth model included all covariates from model 3 in addition to the baseline consumption of total olive oil, VOO, or COO from which the metabolite set was derived, respectively. For the 1-year analyses, we used the same models as at baseline excluding those T2D or CVD cases diagnosed during the first year. Interactions between olive oil profiles and intervention groups were evaluated using the likelihood ratio test including the interaction product terms as covariables.

Several sensitivity analyses were conducted. First, we performed a principal component analysis (PCA) using the mean ENR coefficients from the metabolites consistently selected (i.e., 10 times) in each olive oil profile. A zero value was assigned whenever a particular metabolite was not found by a specific approach. Coefficients were centered and scaled before PCA analysis. Second, we evaluated the specificity of each metabolite profile using Pearson correlations between the consumption of each type of olive oil with each metabolomics profile. Third, we conducted stratified analyses by intervention group. Fourth, we further adjusted the multi-variable Cox regression models for coffee and tea intake.

All analyses were considered statistically significant when P < 0.05 and were performed using R version 4.2.2 statistical software (R Foundation for Statistical Computing).

Results

Characteristics of the study participants

The baseline characteristics of the study participants according to tertiles of energy-adjusted total olive oil consumption are shown in Table 1. The mean ± SD consumption of total olive oil at baseline was 22 ± 8 g/day in the lowest tertile compared to 56 ± 10 g/day in the highest tertile. For VOO and COO, respectively, the mean consumption was 9 ± 11 g/day and 12 ± 12 g/day in the first tertile, and 36 ± 27 g/day and 19 ± 25 g/day in the third tertile. In the top tertile, a higher percentage of participants had T2D, were less likely to smoke, and consumed lower amounts of fruits, vegetables, legumes, cereals, meat, and wine, as compared to participants in the lowest tertile.

Table 1 Baseline characteristics of study participants according to tertiles of energy-adjusted total olive oil consumption

Identification of metabolites associated with olive oil consumption

Cross-sectional associations between baseline plasma metabolites and baseline total olive oil and VOO consumption are shown in Figs. 1 and 2, and Additional file 1: Table S1. Associations between baseline plasma metabolites and baseline common olive oil consumption are shown in Additional file 1: Figure S2 and Additional file 1: Table S1. 78 metabolites were individually significantly associated with total olive oil consumption, 28 metabolites were significantly associated with VOO consumption, and 5 metabolites were significantly associated with COO consumption (FDR < 0.05).

Fig. 1
figure 1

Volcano plot showing the associations between plasma metabolites and total olive oil consumption at baseline. The models were adjusted by recruiting center, smoking status (former smoker, never smoker, smoker), sex, BMI, age, education, and physical activity (METs/day). An FDR < 0.05 was considered statistically significant (up dotted line)

Fig. 2
figure 2

Volcano plot showing the associations between plasma metabolites and virgin olive oil consumption at baseline. The models were adjusted by recruiting center, smoking status (former smoker, never smoker, smoker), sex, BMI, age, education, and physical activity (METs/day). An FDR < 0.05 was considered statistically significant (up dotted line)

Figures 3 and 4 and Additional file 1: Figure S3 show the coefficients (mean and SD) for the ten times selected metabolites in the 10-cross validation of the continuous elastic regression for total olive oil, VOO, and COO. For total olive oil, the metabolites with the strongest inverse associations were C14:2 carnitine, C16:1 cholesterol ester (CE), and cotinine. For VOO, lactose, malate, and sphinganine exhibited the highest negative correlation coefficients; for COO, the strongest inverse associations were found for piperine, pantothenate, and adenosine diphosphate (ADP). For total olive oil, the strongest positive associations were observed for C18:1 carnitine, C34:1 PC plasmalogen A, and C34:2 PC plasmalogen; for VOO, piperine, C34:1 PC plasmalogen A, and ADP showed the strongest direct associations; and for COO, malate, sphingosine, and ornithine were the metabolites with the strongest direct associations. Only piperine was selected in the three metabolomic profiles. This metabolite was positively associated with total olive oil and VOO, but not with COO. Figure 5 shows Venn diagrams showing the number of overlapping and different metabolites for total and subtypes of olive oil identified using elastic net continuous regressions. Thirty metabolites were selected only in the total olive oil profile, thirty-four were only in the VOO profile, and six were only selected in the COO profile.

Fig. 3
figure 3

Coefficients (mean and SD) for the metabolites were selected ten times in the 10-cross validation of the continuous elastic regression for energy-adjusted total olive oil consumption. The sets of metabolites were selected using elastic continuous regression models (with lambda.min) employing the whole dataset of subjects (n = 1833). Negative coefficients are plotted on the left, whereas positive coefficients are shown on the right

Fig. 4
figure 4

Coefficients (mean and SD) for the metabolites were selected ten times in the 10-cross validation of the continuous elastic regression for energy-adjusted virgin oil. The sets of metabolites were selected using elastic continuous regression models (with lambda.min) employing the whole dataset of subjects (n = 1833). Negative coefficients are plotted on the left, whereas positive coefficients are shown on the right

Fig. 5
figure 5

Venn diagram showing the overlapping and different metabolites for total and subtypes of olive oil consumption were identified using elastic net continuous regressions. TOO total olive oil, VOO virgin olive oil, COO common olive oil

Furthermore, to identify the principal components consisting of metabolites most associated with total olive oil, VOO, and/or COO, we also performed a PCA using the coefficients of the metabolites selected by the different olive oil consumption profiles with ENR (Additional file 1: Figure S4). The first principal component (1PC) accounted for 53% of the variability, while the second principal component (2PC) accounted for 36.3% of the variability. 1PC differentiated total olive oil and VOO profiles of the COO profile while 2PC showed differences between total olive oil and COO profiles versus VOO profile. In the PCA biplot (Additional file 1: Figure S4), we observed clusters of metabolites clustered close to the three different profiles.

Table 2 shows the Pearson correlation coefficients between each consistently selected metabolite and the consumption of total olive oil, VOO, and COO in the PREDIMED baseline data (discovery population) and 1-year data (validation sample). The Pearson correlations between the metabolite profiles and energy-adjusted olive oil consumption derived from the FFQ at baseline were 0.40 (95% CI: 0.37, 0.44) for total olive oil consumption, 0.37 (95% CI: 0.33, 0.41) for VOO consumption, and 0.25 (95% CI: 0.21, 0.29) for COO consumption. At 1 year, the Pearson correlations were 0.27 (95% CI: 0.22, 0.31) for total olive oil, 0.23 (95% CI: 0.18, 0.28) for VOO, and 0.16 (95% CI: 0.12, 0.21) for COO. Total olive oil consumption was associated with 74 metabolites, VOO with 78, and COO with 17. Additional file 1: Figure S5 shows a correlation plot between the FFQ-derived consumption of total and subtypes of olive oil and metabolite profiles.

Table 2 Pearson correlation coefficients between metabolomics signatures and olive oil consumption

Associations between the identified metabolomic profiles of olive oil consumption and the risk of T2D and CVD

Table 3 shows the prospective associations between the identified metabolomic profiles of total and subtypes of olive oil consumption and the incidence of T2D and CVD. No significant associations were observed for baseline and 1-year olive oil metabolite profiles and T2D risk. The identified baseline metabolomic profiles of total and subtypes of olive oil showed significant associations with CVD incidence. After adjusting for lifestyle and dietary risk factors, the HR for CVD and 95%CI (for every 1 SD increase) was 0.79 (0.67, 0.92; P-value = 0.003) for the total olive oil metabolite profile, and 0.70 (0.59, 0.83; P-value =  < 0.001) for the VOO, but 1.37 (1.15, 1.63; P-value =  < 0.001) for the COO. In sensitivity analysis, additionally adjusting for coffee and tea consumption, the results remained consistent. Only the 1-year VOO metabolite profile was inversely associated with CVD risk (0.81; 0.65, 1.00; P-value = 0.049).

Table 3 Hazard ratios (95% CIs) for incident type 2 diabetes and cardiovascular diseases according to multi-metabolite associated with 1 SD of olive oil consumption in the PREDIMED study

Interactions between olive oil metabolomic profiles and intervention groups were significant in the CVD models (P-values < 0.05) and non-significant in the T2D models (P-values > 0.05). In the stratified analysis by intervention group (Additional file 1: Tables S2), no significant associations were observed between the olive oil metabolomic profiles and T2D incidence for any of the study intervention groups. However, inverse associations between the metabolite profile of total olive oil and VOO at baseline and the risk of developing CVD were stronger in participants allocated to the MedDiet + VOO group (Additional file 1: Table S3). A direct association between the metabolite signature of COO and a higher risk of CVD was observed. In those participants allocated to the MedDiet + nuts group, only the COO profile showed a positive significant association with CVD incidence. In participants allocated to the control group, an inverse association was observed for the VOO metabolite profile and risk of CVD. At 1 year, only in those participants allocated to the MedDiet + VOO group, a significant inverse association between the total olive oil, and VOO metabolomics signature and CVD was observed, and the COO metabolomic profile was directly associated with CVD risk (Additional file 1: Table S3).

Discussion

Leveraging the integrated dietary and metabolomics data in a Mediterranean population at high cardiovascular disease risk, using an agnostic machine-learning approach, we identified a metabolomic signature for the consumption of total olive oil (with 74 metabolites), VOO (with 78 metabolites), and COO (with 17 metabolites). The metabolomics signature included several lipids, acylcarnitines, and amino acids. Notably, we identified several overlapping and distinct metabolites according to the type of olive oil consumed. Among the significant findings, certain lipids, including plasmalogens, triacylglycerol (TAG), and the organic acid ADP, exhibited positive associations with total olive oil and virgin olive oil (VOO), while no such correlations were observed with common olive oil (COO). Furthermore, we observed a strong and consistent positive correlation between the multi-metabolite profiles and olive oil consumption as assessed through food-frequency questionnaires (FFQ) at both baseline and one year. Linking these metabolites to disease risk, the metabolite profiles associated with total olive oil and VOO showed inverse associations with CVD risk after adjusting for sociodemographic and dietary factors including FFQ-derived olive oil consumption. However, no significant associations between the multi-metabolite profiles and T2D were observed. To the best of our knowledge, this is the first study evaluating the association between plasma metabolite profiles of total and specific types of olive oil and the risk of cardiometabolic diseases. Overall, our findings provide novel insights into the health-promoting benefits of olive oil consumption and highlight the potential relevance of olive oil-related metabolites in relation to chronic disease risk. These results pave the way for further research on the specific metabolic pathways impacted by olive oil consumption and its implications for preventing cardiometabolic diseases.

Olive oil contains a wide variety of minor phytochemicals, such as tocopherols, carotenoids, and phenolic compounds, with recognized biological activity. Compared to common olive oil (a mixture of refined with a minor quantity of virgin olive oil), VOO is obtained exclusively by physical procedures, such as first-pressing and centrifugation, preserving a large part of its phenolic molecules, mainly hydroxytyrosol, oleocanthal, and oleuropein [30, 31]. Several studies identified a dose-dependent relationship between hydroxytyrosol in plasma or urine and VOO consumption [32, 33]. This metabolite was used as a biomarker of compliance with the intervention (MedDiet + VOO) in the PREDIMED study [12] and its biological metabolite (homovanillyl alcohol) was associated with a lower risk of CVD and total mortality [34]. In the EPIC cohort, urinary hydroxytyrosol was also correlated with total olive oil intake [35]. Unfortunately, our metabolomic platform does not allow us to determine plasma phenolic compounds. In addition, hydroxytyrosol can also be obtained from the endogenous hydroxylation of tyrosol, another common phenolic compound present in beer and wine [36, 37]. The determination of metabolite profiles of total olive oil, VOO, or COO consumption may reflect more specifically the variety of olive oil consumed and the metabolic pathways implicated after its consumption.

In previous analysis within the context of the PREDIMED Study, which evaluated the overall Mediterranean dietary pattern [38], several metabolites selected in the signature were found to be associated with olive oil consumption. Specifically, the consumption of olive oil, measured by the Mediterranean Diet Adherence Score (MEDAS), showed positive correlations with various lipids (C20:5 CE, C22:5 CE, some plasmalogens, and C24:1 SM), C18:1 carnitine and 4-pyridoxate, and negative correlations with other lipid molecules including TAGs, PC, CE or DAG, C4 and C14:2 carnitines, glycine, and pyroglutamic acid. In the present study, we also observed some of the previously mentioned associations in relation to different types of reported olive oil consumption. For instance, C20:5 CE, C18:1 carnitine, and various PCs were positively associated with VOO and/or total olive oil consumption, while C14:2 carnitine or C20:4 LPE were associated negatively. Interestingly, our findings indicate that there was only one metabolite on this list that showed an association with COO consumption. This could be attributed to the relatively low overall consumption of this type of olive oil in our study population when compared to VOO.

While the majority of the metabolites identified in this analysis are involved in internal metabolism, it is noteworthy that some of them originate from food sources or result from microbial activity in the gastrointestinal tract. For instance, piperine showed a positive association with total olive oil and VOO consumption, while exhibiting a negative association with COO consumption. Piperine is an alkaloid found in high concentrations in black pepper [39], which happens to be one of the most commonly used spices in Mediterranean cuisine [40]. Likewise, another interesting finding was the negative association between VOO consumption and ectoine. Ectoine is a secondary metabolite [41, 42] produced by certain bacterial genera, such as Streptomyces spp., and has been identified in small concentrations in the human intestinal microbiota [43]. Its presence in the context of olive oil consumption suggests a potential interaction between the gut microbiome and dietary patterns. Other metabolites derived from intestinal microbiota such as 4-hydroxyhippuric acid, indole-3-propionate, and TMAO were positively related to VOO consumption, and have been associated with the metabolism of polyphenols and the consumption of foods and beverages rich in polyphenols in other studies [44, 45] In previous studies, indole-3-propionate and TMAO were found to have negative associations with VOO and COO consumption, respectively [46]. Additionally, hydroxy cotinine and cotinine (both nicotine derivatives [47], as well as caffeine, AAMU, and N1-metil-2-piridona-5-carboxamide (caffeine derivatives [48, 49]) showed negative associations with total olive oil consumption and/or VOO. These findings may reflect the lower prevalence of smokers and coffee consumers among participants who consume VOO.

Olive oil contains mainly MUFA in the form of oleic acid [50]. Many of the identified metabolites are related to MUFA and lipid metabolic pathways, reflecting the composition of olive oil. Only five lipid metabolites (C14:0 CE, C40:6 PS, C18:2 LPC, C34:5 PC plasmalogen, and C54:5 TAG) were exclusively related to the consumption of COO. C5, C9, and C18:1 carnitines were positively associated with VOO consumption, while C10:2, C12, and C14:2 carnitines were negatively associated. Elevated concentrations of acylcarnitines may be a product of the dysregulation of fatty acid oxidation and mitochondrial function. In the PREDIMED study, medium- and long-chain acylcarnitines have been previously associated with an increased risk of CVD [15], but our results indicated that MedDiet interventions may mitigate the adverse associations shown between higher concentrations of acylcarnitines and CVD.

Among the metabolites related to energy and carbohydrate metabolism, fructose 6-phosphate and inositol were found in the VOO profile, citrate in the total olive oil profile, lactate in both COO and total olive oil profiles, and malate in both VOO and COO profiles. Given that olive oil primarily consists of fatty acids, it is not surprising that only a few carbohydrate-related metabolites were associated with olive oil consumption. Similar observations have been reported in animal studies (51) which showed that diets rich in fats, like the Mediterranean Diet, can influence the Randle cycle, leading to increased malonyl-CoA production from β-oxidation, which serves as a substrate for the TCA cycle and/or gluconeogenesis [52]. Plasma glycocholate, glycodeoxycholic acid, bilirubin, and biliverdin levels are synthesized conjugated bile acids that have been positively associated with total olive oil and/or VOO in our study. Olive oil acts on the gallbladder providing its complete emptying, stimulating the synthesis of bile salts in the liver, and increasing the hepatic secretion of cholesterol [53], thus potentially explaining the association found in our study.

Our signatures have also identified several metabolites related to purine pathways. In our study, 1-methylguanine, urate, indoxyl sulfate (a uremic solute), inosine, adenosine, and ADP were positively associated with VOO or total olive oil consumption. It has been suggested that higher levels of urate in subjects with CVD may represent a compensatory response to counteract oxidative stress [54]. Another oxidative stress marker such as allantoin was negatively associated with VOO but positively associated with common olive oil [55], which suggests that VOO has a key role in these pathways with an antioxidant effect. Some metabolites were identified only in the VOO signature as adipate, acetylcholine, γ-butyrobetaine, 2-aminoadipate, DMGV, pantothenol or trimethyl benzene, and other metabolites only with total olive oil consumption as cytosine or pipecolic acid. We did not find an explanation in the existing literature for why these metabolites showed associations with VOO but not with COO.

These unexplained associations between specific metabolites and VOO or total olive oil consumption could potentially be attributed to the complex interactions between dietary components and individual metabolic responses.

Some metabolites that are part of the identified metabolite profiles have previously been associated with CVD and T2D [56, 57], potentially explaining the positive effects that have been reported for olive oil consumption, especially VOO, on cardiometabolic health [1]. For example, C24:0 ceramide or α-hydroxybutyrate have previously been associated with insulin resistance and increased risk of T2D, and C16:1 CE, several acylcarnitines, cortisol, or deoxycortisol (intermediate of cortisol) with increased CVD risk [56,57,58,59,60]. All these metabolites have been inversely associated with VOO consumption in the signature. However, other lipids such as C54:2 TAG, and C36:2 PC plasmalogen, which have also been positively associated with the consumption of total olive oil or VOO, have been associated with an increased risk of CVD and T2D. After adjusting for potential confounders, we found that the total olive oil and VOO metabolic profiles were associated with a 21% and 26% lower CVD risk, respectively, whereas the COO metabolic profile was associated with a 26% higher risk of CVD risk. However, after 1-year, no statistically significant associations between the total olive oil and COO metabolomic profiles and CVD were observed, similar to the previous analysis of the PREDIMED [61]. In the context of a MedDiet, VOO has been demonstrated to improve lipid profile, markers of glucose and inflammation, and decreased blood pressure, all considered CVD risk factors [1, 62,63,64]. VOO is known to have a higher concentration of polyphenols compared to COO, and these bioactive compounds are associated with various health benefits, including antioxidant and anti-inflammatory effects or improvements in lipid profile [2]. These observed differences in varying polyphenol content between olive oil subtypes may potentially explain the differences between the metabolomic profiles of VOO and COO and the risk of CVD. Additionally, the statistical differences between VOO and COO profiles could be influenced by the consumption patterns of the study participants, where the lower consumption of COO compared to the consumption of VOO implies an opposite distribution of CVD events. Curiously, similar significant results were seen in the MedDiet groups but not in the control group when stratified by intervention groups. On the other hand, no significant associations were observed between the metabolomic signatures and T2D, and some controversial findings have been reported in the literature for the associations between olive oil consumption and T2D. In the PREDIMED study, while participants allocated to the MedDiet + VOO had a lower risk of developing T2D, the consumption of olive oil alone has not been associated with T2D [65, 66]. Reverse causality might explain these findings. It is possible that participants with higher olive oil consumption, leading to higher T2DM prevalence, had a better overall lifestyle and diet compared to others in the study, influencing the observed outcomes. This highlights the complexity of dietary research and the need for careful consideration of confounding factors.

These findings need to be interpreted in the context of some limitations. First, the identified metabolite profiles are not an objective biomarker of olive oil per se, but they reflect the overall homeostasis associated with olive oil consumption, the substitution of other food by the consumption of olive oil, and the individual biological responses to diet. In addition, since olive oil consumption is usually accompanied by other foods, some of the selected metabolites may be associated with the consumption of other foods. Second, the metabolite profiles were derived from a pool of 385 annotated metabolites, while thousands of unique metabolites have been identified to date. We cannot exclude that more biologically relevant metabolites regarding olive oil intake were absent from our data set. For example, the metabolomics approach used for quantifying lipids did not identify the specific fatty acids for each molecule (we can only provide the number of carbons and double bonds of each lipid), and measures of polyphenols and phytochemicals are not available; consequently, some relevant olive oil biomarkers may have been missed. Therefore, the specificity of the identified metabolite profiles of olive oil intake remains uncertain. Future studies are warranted to identify additional objective biomarkers of olive oil intake, including urinary metabolites from dietary intake of phenolic compounds that will be assessed in the future.

Third, due to the use of FFQs for collecting dietary data, measurement errors may be present compared to the use of short-term biomarkers of intake. However, the validity and reproducibility of the FFQ have been reported previously. Of note, the correlation between total olive oil intake assessed by our FFQ and 3-d dietary records was relatively high (r = 0.60) [18]. Because of our study's observational design, we cannot establish the causality of the association between the metabolomic signatures and cardiometabolic diseases. Nevertheless, we performed a rigorous multivariable adjustment to minimize residual confounding. Further, although we evaluated the cross-population reproducibility of the metabolite profiles, it should be validated in independent populations. Our population was a Mediterranean population with a high consumption of VOO even before being randomized in the study. Therefore, we cannot be ruled out reverse causation in participants with diabetes at baseline.

The current study also has several strengths. The PREDIMED study is an ideal setting to identify metabolite profiles associated with olive oil consumption because the consumption of olive oil at baseline (mean total olive oil at baseline in this population is 40 g/d) is much higher than in other populations. In addition, we were able to differentiate between subtypes of olive oil, which have different nutritional compositions, this is a limitation of other existing studies where detailed data on specific subtypes of olive oil and comprehensive metabolomics data are not available. The present study has a large sample size and detailed covariate data to control for confounding and well-defined outcomes. In addition, we used agnostic machine learning models using more than > 350 well-annotated metabolites. We cross-validated our results internally in the discovery population using baseline data and conducted replication analysis using data at year 1.

Conclusions

In summary, our study analyzed 385 candidate metabolites and identified distinct panels associated with the consumption of total olive oil, VOO, and COO. Specifically, we identified a metabolomic signature associated with the consumption of total olive oil (with 74 metabolites), VOO (with 78 metabolites), and COO (with 17 metabolites). Moreover, our findings highlight the significance of the VOO metabolite profile in reducing the risk of CVD in Mediterranean individuals with high cardiovascular risk. These novel insights provide potential biomarkers of olive oil consumption (especially VOO) and offer valuable information on the mechanisms underlying the relationship between olive oil consumption and cardiometabolic diseases, including CVD and type 2 diabetes (T2D). By advancing our understanding of the metabolic responses to olive oil consumption, our research contributes to the broader field of nutrition and health, emphasizing the health benefits of incorporating extra virgin olive oil into dietary patterns.

Fundings

The PREDIMED study was funded by NIH grants R01 HL118264, R01 DK102896 and R01 3R01 HL118264, the Spanish Ministry of Health (Instituto de Salud Carlos III, The PREDIMED Network grant RD 06/0045, 2006–2013, coordinated by M.A. Martínez-González; and a previous network grant RTIC-G03/140, 2003–2005, coordinated R. Estruch). Additional grants were received from the Ministerio de Economía y Competitividad-Fondo Europeo de Desarrollo Regional (Projects CNIC-06/2007, CIBER 06/03, PI06-1326, PI07-0954, PI11/02505, SAF2009-12304, and AGL2010–22319-C03-03) and by the Generalitat Valenciana (PROMETEO17/2017 and PROMETEO 21/2021). Dr. Guasch-Ferré acknowledges that The Novo Nordisk Foundation Center for Basic Metabolic Research is supported by and unrestricted grant from the Novo Nordisk Foundation (grant no. NNF18CC0034900). Dr. Salas-Salvadó gratefully acknowledges the financial support by ICREA under the ICREA Academia program.

Availability of data and materials

The dataset generated and/or analyzed during the current study is not publicly available due to the lack of authorization from PREDIMED participants. Those wishing to access the PREDIMED trial data used in this study can request the corresponding author and it will then be passed on to members of the PREDIMED Steering Committee for deliberation.

Abbreviations

95% CI:

95% Confidence interval

AAMU:

5-Acetamido-6-amino-3-methyl-uracil

ADP:

Adenosine diphosphate

BMI:

Body mass index

CE:

Cholesterol ester

COO:

Common olive oil

CV:

Cross-validation

CVD:

Cardiovascular diseases

DAG:

Diacylglycerol

DMGV:

Dimethylguanidino valeric acid

VOO:

Virgin olive oil

FFQ:

Food frequency questionnaires

HR:

Hazard ratio

LC-MS:

Liquid chromatography-tandem mass spectrometry

LPC:

Lysophosphatidylcholine

LPE:

Lysophosphatidylethanolamine

MAG:

Monoacylglycerol

MEDAS:

Mediterranean Diet Adherence Screener

MedDiet:

Mediterranean Diet

MUFA:

Monounsaturated fatty acid

PC:

Phosphatidylcholine

PREDIMED:

PREvención con DIeta MEDiterránea

PS:

Phosphatidylserine

SM:

Sphingomyelin

T2D:

Type 2 diabetes

TAG:

Triacylglycerol

TCA:

Tricarboxylic acid cycle

TMAO:

Trimethylamine N-oxide

References

  1. Gaforio, Visioli, Alarcón-de-la-Lastra, Castañer, Delgado-Rodríguez, Fitó, et al. Virgin Olive Oil and Health: Summary of the III International Conference on Virgin Olive Oil and Health Consensus Report, JAEN (Spain) 2018. Nutrients. 2019 Sep 1;11(9):2039.

  2. Ros E. Olive oil and CVD: accruing evidence of a protective effect. Br J Nutr. 2012;108(11):1931–3.

    Article  CAS  PubMed  Google Scholar 

  3. Martínez-González MA, Sayón-Orea C, Bullón-Vela V, Bes-Rastrollo M, Rodríguez-Artalejo F, Yusta-Boyo MJ, et al. Effect of olive oil consumption on cardiovascular disease, cancer, type 2 diabetes, and all-cause mortality: a systematic review and meta-analysis. Clin Nutr. 2022;41(12):2659–82.

    Article  PubMed  Google Scholar 

  4. Schwingshackl L, Lampousi AM, Portillo MP, Romaguera D, Hoffmann G, Boeing H. Olive oil in the prevention and management of type 2 diabetes mellitus: a systematic review and meta-analysis of cohort studies and intervention trials. Nutr Diabetes. 2017;7(4):e262–e262.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Guasch-Ferré M, Li Y, Willett WC, Sun Q, Sampson L, Salas-Salvadó J, et al. Consumption of olive oil and risk of total and cause-specific mortality among U.S. Adults J Am Coll Cardiol. 2022;79(2):101–12.

    Article  PubMed  Google Scholar 

  6. Markellos C, Ourailidou ME, Gavriatopoulou M, Halvatsiotis P, Sergentanis TN, Psaltopoulou T. Olive oil intake and cancer risk: A systematic review and meta-analysis. PLoS ONE. 2022;17(1):e0261649.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Schwingshackl L, Krause M, Schmucker C, Hoffmann G, Rücker G, Meerpohl JJ. Impact of different types of olive oil on cardiovascular risk factors: a systematic review and network meta-analysis. Nutr Metab Cardiovasc Dis. 2019;29(10):1030–9.

    Article  CAS  PubMed  Google Scholar 

  8. Yubero-Serrano EM, Lopez-Moreno J, Gomez-Delgado F, Lopez-Miranda J. Extra virgin olive oil: more than a healthy fat. Eur J Clin Nutr. 2019;72(S1):8–17.

    Article  PubMed  Google Scholar 

  9. Hernáez Á, Fernández-Castillejo S, Farràs M, Catalán Ú, Subirana I, Montes R, et al. Olive oil polyphenols enhance high-density lipoprotein function in humans. Arterioscler Thromb Vasc Biol. 2014;34(9):2115–9.

    Article  PubMed  Google Scholar 

  10. Fernández-Castillejo S, Pedret A, Catalán Ú, Valls R, Farràs M, Rubió L, et al. Virgin olive oil phenolic compounds modulate the HDL lipidome in hypercholesterolaemic subjects: a lipidomic analysis of the VOHF study. Mol Nutr Food Res. 2021;65(9):2001192.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Martinez-Gonzalez MA, Corella D, Salas-Salvado J, Ros E, Covas MI, Fiol M, et al. Cohort profile: design and methods of the PREDIMED study. Int J Epidemiol. 2012;41(2):377–85.

    Article  PubMed  Google Scholar 

  12. Estruch R, Ros E, Salas-Salvadó J, Covas MI, Corella D, Arós F, 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(25):e34.

    Article  CAS  PubMed  Google Scholar 

  13. Guasch-Ferré M, Ruiz-Canela M, Li J, Zheng Y, Bulló M, Wang DD, et al. Plasma acylcarnitines and risk of type 2 diabetes in a mediterranean population at high cardiovascular risk. J Clin Endocrinol Metab. 2019;104(5):1508–19.

    Article  PubMed  Google Scholar 

  14. Ruiz-Canela M, Guasch-Ferré M, Toledo E, Clish CB, Razquin C, Liang L, et al. Plasma branched chain/aromatic amino acids, enriched Mediterranean diet and risk of type 2 diabetes: case-cohort study within the PREDIMED trial. Diabetologia. 2018;61(7):1560–71.

    Article  CAS  PubMed Central  Google Scholar 

  15. Guasch-Ferré M, Zheng Y, Ruiz-Canela M, Hruby A, Martínez-González MA, Clish CB, et al. Plasma acylcarnitines and risk of cardiovascular disease: effect of Mediterranean diet interventions. Am J Clin Nutr. 2016;103(6):1408–16.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Ruiz-Canela M, Toledo E, Clish CB, Hruby A, Liang L, Salas-Salvadó J, et al. Plasma branched-chain amino acids and incident cardiovascular disease in the PREDIMED trial. Clin Chem. 2016;62(4):582–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Willet W. Nutritional epidemiology. 3rd ed. New York: Oxford University Press; 2012.

    Book  Google Scholar 

  18. Fernández-Ballart JD, Piñol JL, Zazpe I, Corella D, Carrasco P, Toledo E, et al. Relative validity of a semi-quantitative food-frequency questionnaire in an elderly Mediterranean population of Spain. Br J Nutr. 2010;103(12):1808–16.

    Article  PubMed  Google Scholar 

  19. Mataix J. Tablas de composición de alimentos. 4th ed. Granada, Spain: Universidad de Granada; 2003.

    Google Scholar 

  20. Moreiras O, Carvajal A, Cabrera L. Tablas de Composición de Alimentos [Food Composition Tables]. 9th ed. Moreiras O, Carvajal A, Cabrera L, editors. Madrid, Spain: Ediciones Pirámide; 2005.

  21. Elosua R, Marrugat J, Molina L, Pons S, Pujol E. Validation of the minnesota leisure time physical activity questionnaire in spanish men. Am J Epidemiol. 1994;139(12):1197–209.

    Article  CAS  PubMed  Google Scholar 

  22. Hernández-Alonso P, Papandreou C, Bulló M, Ruiz-Canela M, Dennis C, Deik A, et al. Plasma metabolites associated with frequent red wine consumption: a metabolomics approach within the PREDIMED study. Mol Nutr Food Res. 2019;63(17):1900140.

    Article  Google Scholar 

  23. O’Sullivan JF, Morningstar JE, Yang Q, Zheng B, Gao Y, Jeanfavre S, et al. Dimethylguanidino valeric acid is a marker of liver fat and predicts diabetes. J Clin Investig. 2017;127(12):4394–402.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Paynter NP, Balasubramanian R, Giulianini F, Wang DD, Tinker LF, Gopal S, et al. Metabolic predictors of incident coronary heart disease in women. Circulation. 2018;137(8):841–53.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Wang TJ, Larson MG, Vasan RS, Cheng S, Rhee EP, McCabe E, et al. Metabolite profiles and the risk of developing diabetes. Nat Med. 2011;17(4):448–53.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Stekhoven DJ, Buhlmann P. MissForest–non-parametric missing value imputation for mixed-type data. Bioinformatics. 2012;28(1):112–8.

    Article  CAS  PubMed  Google Scholar 

  27. Gromski P, Xu Y, Kotze H, Correa E, Ellis D, Armitage E, et al. Influence of missing values substitutes on multivariate analysis of metabolomics data. Metabolites. 2014;4(2):433–52.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Wei R, Wang J, Su M, Jia E, Chen S, Chen T, et al. Missing value imputation approach for mass spectrometry-based metabolomics data. Sci Rep. 2018;8(1):663.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Blom G. Statistical estimates and transformed beta-variables. New York: Wiley; 1958.

    Google Scholar 

  30. Tripoli E, Giammanco M, Tabacchi G, Di Majo D, Giammanco S, La Guardia M. The phenolic compounds of olive oil: structure, biological activity and beneficial effects on human health. Nutr Res Rev. 2005;18(1):98–112.

    Article  CAS  PubMed  Google Scholar 

  31. Lozano-Castellón J, López-Yerena A, Rinaldi de Alvarenga JF, Romero del Castillo-Alba J, Vallverdú-Queralt A, Escribano-Ferrer E, et al. Health-promoting properties of oleocanthal and oleacein: two secoiridoids from extra-virgin olive oil. Crit Rev Food Sci Nutr. 2020;60(15):2532–48.

    Article  PubMed  Google Scholar 

  32. Garcia-Aloy M, Hulshof PJM, Estruel-Amades S, Osté MCJ, Lankinen M, Geleijnse JM, et al. Biomarkers of food intake for nuts and vegetable oils: an extensive literature search. Genes Nutr. 2019;14(1):7.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Nikou T, Sakavitsi ME, Kalampokis E, Halabalaki M. Metabolism and bioavailability of olive bioactive constituents based on in vitro, in vivo and human studies. Nutrients. 2022;14(18):3773.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. De la Torre R, Corella D, Castañer O, Martínez-González MA, Salas-Salvador J, Vila J, et al. Protective effect of homovanillyl alcohol on cardiovascular disease and total mortality: virgin olive oil, wine, and catechol-methylathion. Am J Clin Nutr. 2017;105(6):1297–304.

    Article  PubMed  Google Scholar 

  35. Almanza-Aguilera E, Davila-Cordova E, Guiñón-Fort D, Farràs M, Masala G, Santucci de Magistris M, et al. Correlation analysis between dietary intake of tyrosols and their food sources and urinary excretion of tyrosol and hydroxytyrosol in a european population. Antioxidants. 2023;12(3):715.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Soldevila-Domenech N, Boronat A, Mateus J, Diaz-Pellicer P, Matilla I, Pérez-Otero M, et al. Generation of the antioxidant hydroxytyrosol from tyrosol present in beer and red wine in a randomized clinical trial. Nutrients. 2019;11(9):2241.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Boronat A, Mateus J, Soldevila-Domenech N, Guerra M, Rodríguez-Morató J, Varon C, et al. Cardiovascular benefits of tyrosol and its endogenous conversion into hydroxytyrosol in humans. A randomized, controlled trial. Free Radic Biol Med. 2019;143:471–81.

    Article  CAS  PubMed  Google Scholar 

  38. Li J, Guasch-Ferré M, Chung W, Ruiz-Canela M, Toledo E, Corella D, et al. The Mediterranean diet, plasma metabolome, and cardiovascular disease risk. Eur Heart J. 2020;41(28):2645–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Meghwal M, Goswami TK. Piper nigrum and Piperine: an Update. Phytother Res. 2013;27(8):1121–30.

    Article  CAS  PubMed  Google Scholar 

  40. Talib WH, AlHur MJ, Al Naimat S, Ahmad RE, Al-Yasari AH, Al-Dalaeen A, et al. Anticancer effect of spices used in mediterranean diet: preventive and therapeutic potentials. Front Nutr. 2022. https://doi.org/10.3389/fnut.2022.905658.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Bethlehem L, van Echten-Deckert G. Ectoines as novel anti-inflammatory and tissue protective lead compounds with special focus on inflammatory bowel disease and lung inflammation. Pharmacol Res. 2021;164:105389.

    Article  CAS  PubMed  Google Scholar 

  42. Bownik A, Stępniewska Z. Ectoine as a promising protective agent in humans and animals. Arch Ind Hyg Toxicol. 2016;67(4):260–5.

    CAS  Google Scholar 

  43. Cuozzo S, de LeBlanc AD, LeBlanc JG, Hoffmann N, Tortella GR. Streptomyces genus as a source of probiotics and its potential for its use in health. Microbiol Res. 2023;266:127248.

    Article  CAS  PubMed  Google Scholar 

  44. van Dorsten FA, Grün CH, van Velzen EJJ, Jacobs DM, Draijer R, van Duynhoven JPM. The metabolic fate of red wine and grape juice polyphenols in humans assessed by metabolomics. Mol Nutr Food Res. 2010;54(7):897–908.

    Article  PubMed  Google Scholar 

  45. Pallister T, Jackson MA, Martin TC, Zierer J, Jennings A, Mohney RP, et al. Hippurate as a metabolomic marker of gut microbiome diversity: modulation by diet and relationship to metabolic syndrome. Sci Rep. 2017;7(1):13670.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Tang WHW, Wang Z, Levison BS, Koeth RA, Britt EB, Fu X, et al. Intestinal Microbial Metabolism of Phosphatidylcholine and Cardiovascular Risk. N Engl J Med. 2013;368(17):1575–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Grizzell JA, Echeverria V. New Insights into the Mechanisms of action of cotinine and its distinctive effects from nicotine. Neurochem Res. 2015;40(10):2032–46.

    Article  CAS  PubMed  Google Scholar 

  48. Weimann A, Sabroe M, Poulsen HE. Measurement of caffeine and five of the major metabolites in urine by high-performance liquid chromatography/tandem mass spectrometry. J Mass Spectrom. 2005;40(3):307–16.

    Article  CAS  PubMed  Google Scholar 

  49. Wong P, Bachki A, Banerjee K, Leyland-Jones B. Identification of N1-methyl-2-pyridone-5-carboxamide and N1-methyl-4-pyridone-5-carboxamide as components in urine extracts of individuals consuming coffee. J Pharm Biomed Anal. 2002;30(3):773–80.

    Article  CAS  PubMed  Google Scholar 

  50. Jimenez-Lopez C, Carpena M, Lourenço-Lopes C, Gallardo-Gomez M, Lorenzo JM, Barba FJ, et al. Bioactive compounds and quality of extra virgin olive oil. Foods. 2020;9(8):1014.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. An Y, Xu W, Li H, Lei H, Zhang L, Hao F, et al. High-fat diet induces dynamic metabolic alterations in multiple biological matrices of rats. J Proteome Res. 2013;12(8):3755–68.

    Article  CAS  PubMed  Google Scholar 

  52. Randle PJ, Garland PB, Hales CN, Newsholme EA. The glucose fatty-acid cycle its role in insulin sensitivity and the metabolic disturbances of Diabetes mellitus. The Lancet. 1963;281(7285):785–9.

    Article  Google Scholar 

  53. Krzeminski R, Gorinstein S, Leontowicz H, Leontowicz M, Gralak M, Czerwinski J, et al. Effect of different olive oils on bile excretion in rats fed cholesterol-containing and cholesterol-free diets. J Agric Food Chem. 2003;51(19):5774–9.

    Article  CAS  PubMed  Google Scholar 

  54. Bae JS, Shin DH, Park PS, Choi BY, Kim MK, Shin MH, et al. The impact of serum uric acid level on arterial stiffness and carotid atherosclerosis: The Korean Multi-Rural Communities Cohort study. Atherosclerosis. 2013;231(1):145–51.

    Article  CAS  PubMed  Google Scholar 

  55. Lotito S, Frei B. Consumption of flavonoid-rich foods and increased plasma antioxidant capacity in humans: cause, consequence, or epiphenomenon? Free Radic Biol Med. 2006;41(12):1727–46.

    Article  CAS  PubMed  Google Scholar 

  56. Morze J, Wittenbecher C, Schwingshackl L, Danielewicz A, Rynkiewicz A, Hu FB, et al. Metabolomics and type 2 diabetes risk: an updated systematic review and meta-analysis of prospective cohort studies. Diabetes Care. 2022;45(4):1013–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. McGranaghan P, Saxena A, Rubens M, Radenkovic J, Bach D, Schleußner L, et al. Predictive value of metabolomic biomarkers for cardiovascular disease risk: a systematic review and meta-analysis. Biomarkers. 2020;25(2):101–11.

    Article  PubMed  Google Scholar 

  58. Ruiz-Canela M, Guasch-Ferré M, Razquin C, Toledo E, Hernández-Alonso P, Clish CB, et al. Plasma acylcarnitines and risk of incident heart failure and atrial fibrillation: the Prevención con dieta mediterránea study. Revista Española de Cardiología (English Edition). 2022;75(8):649–58.

    Article  Google Scholar 

  59. Iob E, Steptoe A. Cardiovascular disease and hair cortisol: a novel biomarker of chronic stress. Curr Cardiol Rep. 2019;21(10):116.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Gall WE, Beebe K, Lawton KA, Adam KP, Mitchell MW, Nakhle PJ, et al. α-Hydroxybutyrate Is an early biomarker of insulin resistance and glucose intolerance in a nondiabetic population. PLoS ONE. 2010;5(5):e10883.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Guasch-Ferré M, Hu FB, Martínez-González MA, Fitó M, Bulló M, Estruch R, et al. Olive oil intake and risk of cardiovascular disease and mortality in the PREDIMED Study. BMC Med. 2014;12(1):78.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Urpi-Sarda M, Casas R, Chiva-Blanch G, Romero-Mamani ES, Valderas-Martínez P, Arranz S, et al. Virgin olive oil and nuts as key foods of the Mediterranean diet effects on inflammatory biomakers related to atherosclerosis. Pharmacol Res. 2012;65(6):577–83.

    Article  CAS  PubMed  Google Scholar 

  63. Gimeno E, de la Torre-Carbot K, Lamuela-Raventós RM, Castellote AI, Fitó M, de la Torre R, et al. Changes in the phenolic content of low density lipoprotein after olive oil consumption in men. A randomized crossover controlled trial. Br J Nutr. 2007;98(6):1243–50.

    Article  CAS  PubMed  Google Scholar 

  64. Bullo M, Lamuela-Raventos R, Salas-Salvado J. Mediterranean diet and oxidation: nuts and olive oil as important sources of fat and antioxidants. Curr Top Med Chem. 2011;11(14):1797–810.

    Article  CAS  PubMed  Google Scholar 

  65. Salas-Salvadó J, Bulló M, Babio N, Martínez-González MÁ, Ibarrola-Jurado N, Basora J, et al. Reduction in the incidence of type 2 diabetes with the mediterranean diet. Diabet Care. 2011;34(1):14–9.

    Article  Google Scholar 

  66. Salas-Salvadó J, Bulló M, Estruch R, Ros E, Covas MI, Ibarrola-Jurado N, et al. Prevention of diabetes with mediterranean diets. Ann Intern Med. 2014;160(1):1–10.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors are grateful to the PREDIMED participants and personnel, as well as the affiliated primary care center staff, for their altruistic collaboration and excellent assistance.

Author information

Authors and Affiliations

Authors

Contributions

JGG, MGF, MR-C, CR, ET, MAM-G, JSS, FBH designed the research; DC, RE, MF, LSM, MAM-G, JSS coordinated the subject recruitment at the outpatient clinics and clinical data collection in Prevención con Dieta Mediterránea (PREDIMED); CBC conducted the metabolomics data analysis; JGG conducted the statistical analysis; JGG, MGF, and JSS drafted the manuscript; JGG, and JSS, are the guarantors, and, as such, take responsibility for the integrity of the data and the accuracy of the data analysis; JGG, MAM-G, JSS, FBH had access to all the data in the study; JGG, NB, MF, J-PD-C, ET, CW, MAM-G, JSS, FBH interpreted the data; and all authors made critical revisions to the manuscript for key intellectual content and read and approved the final manuscript.

Corresponding authors

Correspondence to Jesús F. García-Gavilán or Marta Guasch-Ferré.

Ethics declarations

Ethics approval and consent to participate

The Research Ethics Committee of the Institut d’Investigacions Sanitàries Pere i Virgili approved the PREDIMED study protocol and all participants provided written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare no conflict of interest in this article. The California Walnut Commission, Borges SA, and Morella Nuts SA donated nuts for the PREDIMED study. None of the funding sources played a role in the design, collection, analysis, or interpretation of the data or in the decision to submit the manuscript for publication.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file1

: Table S1. Associations between plasma metabolites and total olive oil, extra virgin olive oil, and common olive oil consumption at baseline. Table S2. Hazard ratios (95% CIs) of type 2 diabetes incidence according to metabolomics profiles of energy-adjusted olive oil and its subtypes in the PREDIMED study groups stratified by intervention group. Table S3. Hazard ratios (95% CIs) of cardiovascular disease incidence according to metabolomics profiles of energy-adjusted olive oil and its subtypes in the PREDIMED study groups stratified by intervention group. Figure S1. Flowchart of participants and analysis. Figure S2. Volcano plot showing the associations between plasma metabolites and common olive oil consumption at baseline. Figure S3. Metabolite coefficients (mean and SD) selected ten times in the 10-cross validation of the continuous elastic regression for energy-adjusted common olive oil. Figure S4. Biplot of the principal component analysis using metabolites’ coefficients derived by the elastic net continuous regression of each olive oil consumption approach. Figure S5. Correlation plot showing the correlation between each self-reported olive oil intake and olive oil metabolomics profiles and each food group consumption.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

García-Gavilán, J.F., Babio, N., Toledo, E. et al. Olive oil consumption, plasma metabolites, and risk of type 2 diabetes and cardiovascular disease. Cardiovasc Diabetol 22, 340 (2023). https://doi.org/10.1186/s12933-023-02066-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12933-023-02066-1

Keywords