The association between adherence to the Mediterranean diet and hepatic steatosis: cross-sectional analysis of two independent studies, the UK Fenland Study and the Swiss CoLaus Study
BMC Medicine volume 17, Article number: 19 (2019)
Background and aims
The risk of hepatic steatosis may be reduced through changes to dietary intakes, but evidence is sparse, especially for dietary patterns including the Mediterranean diet. We investigated the association between adherence to the Mediterranean diet and prevalence of hepatic steatosis.
Cross-sectional analysis of data from two population-based adult cohorts: the Fenland Study (England, n = 9645, 2005–2015) and CoLaus Study (Switzerland, n = 3957, 2009–2013). Habitual diet was assessed using cohort-specific food frequency questionnaires. Mediterranean diet scores (MDSs) were calculated in three ways based on adherence to the Mediterranean dietary pyramid, dietary cut-points derived from a published review, and cohort-specific tertiles of dietary consumption. Hepatic steatosis was assessed by abdominal ultrasound and fatty liver index (FLI) in Fenland and by FLI and non-alcoholic fatty liver disease (NAFLD) score in CoLaus. FLI includes body mass index (BMI), waist circumference, gamma-glutamyl transferase, and triglyceride; NAFLD includes diabetes, fasting insulin level, fasting aspartate-aminotransferase (AST), and AST/alanine transaminase ratio. Associations were assessed using Poisson regression.
In Fenland, the prevalence of hepatic steatosis was 23.9% and 27.1% based on ultrasound and FLI, respectively, and in CoLaus, 25.3% and 25.7% based on FLI and NAFLD score, respectively. In Fenland, higher adherence to pyramid-based MDS was associated with lower prevalence of hepatic steatosis assessed by ultrasound (prevalence ratio (95% confidence interval), 0.86 (0.81, 0.90) per one standard deviation of MDS). This association was attenuated [0.95 (0.90, 1.00)] after adjustment for body mass index (BMI). Associations of similar magnitude were found for hepatic steatosis assessed by FLI in Fenland [0.82 (0.78, 0.86)] and in CoLaus [0.85 (0.80, 0.91)], and these were also attenuated after adjustment for BMI. Findings were similar when the other two MDS definitions were used.
Greater adherence to the Mediterranean diet was associated with lower prevalence of hepatic steatosis, largely explained by adiposity. These findings suggest that an intervention promoting a Mediterranean diet may reduce the risk of hepatic steatosis.
Hepatic steatosis is a major cause of chronic liver disease worldwide, with prevalence estimates ranging from 25 to 45% in the general population . Hepatic steatosis, usually defined as fat accumulation > 5% in hepatocytes, is the first recognisable stage for both alcoholic and non-alcoholic fatty liver disease (NAFLD). Hepatic steatosis, especially NAFLD, may lead to progressive liver fibrosis and cirrhosis and increased risk of type 2 diabetes and cardiovascular diseases . NAFLD is the primary hepatic outcome of metabolic syndrome and further cardiometabolic diseases, which include insulin resistance as well as dyslipidaemia and obesity as key pathologic mechanisms according to recent scientific advances [3, 4]. Previous dietary studies reported that, compared with healthy controls, individuals with hepatic steatosis had a higher consumption of carbohydrates, saturated fat, meat, and soft drinks, along with a tendency towards a lower intake of omega-3 polyunsaturated fatty acids (PUFA) and fish [5,6,7,8]. However, there are few data concerning associations between dietary factors and overall dietary pattern with hepatic steatosis among healthy adults.
Recently, the Mediterranean diet has been recommended for the management of NAFLD by the EASL-EASD-EASO Clinical Practice Guideline . Meta-analyses and systematic reviews of published evidence have supported a beneficial impact of adherence to the Mediterranean diet on overall mortality and on risks of metabolic syndrome, type 2 diabetes, and cardiovascular disease [9,10,11]. However, whether the Mediterranean diet is associated with risk of hepatic steatosis remains unclear. Higher adherence to the Mediterranean diet has beneficial effects on the progression of hepatic steatosis, but this evidence was derived from small-scale trials (n < 90 followed up for < 6 months) in patients with existing hepatic steatosis [12,13,14,15,16,17]. Among adults free from clinically manifest hepatic steatosis, evidence from three studies is available, but still inconclusive. Two cross-sectional studies reported conflicting results: one among obese Spanish adults with high cardiovascular risk (n = 794) reported an inverse association , and the other among apparently healthy Chinese adults (n = 332) reported a null association . The Framingham Heart Study, only one longitudinal study thus far, has recently reported a significant inverse association of greater adherence to the Mediterranean diet with incident hepatic steatosis (n = 1521 adults over 6 years follow-up) .
Given the limited and inconsistent evidence, which is mainly restricted to a single study or country, we aimed to investigate the cross-sectional association of adherence to the Mediterranean diet with hepatic steatosis among middle-aged healthy adults in two independent population-based cohorts: Fenland (East England, UK) and CoLaus (Lausanne, Switzerland).
Study design and population
The Fenland Study is an ongoing population-based cohort of adults from general practice lists in Cambridgeshire (Cambridge, Ely, Wisbech, and surrounding villages) in the UK . Overall, 12,435 adults born between 1950 and 1975 (aged 30 to 65 years at recruitment) attended baseline clinical assessments in 2005–2015. This cohort was established to investigate environmental and genetic risk factors for the development of type 2 diabetes and related metabolic disorders. Exclusion criteria were the presence of known diabetes, being pregnant, being unable to walk unaided, or having psychosis or a terminal illness. The Cambridge Local Research Ethics Committee approved the study (04/Q0108/19), and all participants provided written informed consent.
The CoLaus Study is an ongoing population-based cohort including a random sample of 6733 individuals aged 35 to 75 years in the city of Lausanne, Switzerland. Details of the study have been described previously . This cohort was established in 2003 to investigate the clinical, biological, and genetic determinants of cardiovascular disease, and it included participants of European origin. For the present analysis of CoLaus, we used data from the first follow-up (study period 2009–2013; n = 5064) when dietary assessment was initiated. The Institutional Ethics Committee of the University of Lausanne approved the study (reference 16/03), and all participants provided written informed consent.
In this analysis of the two cohorts, we used standardised exclusion criteria to remove participants with diabetes [defined as glycated haemoglobin ≥ 48 mmol/mol, fasting plasma glucose ≥ 7.0 mmol/L, 2-h glucose ≥ 11.1 mmol/L, or use of glucose-lowering drugs], those who are pregnant, and those missing dietary data, outcome data, or key covariates. We excluded those with implausible energy intake based on sex-specific thresholds (< 500 or > 3500 kcal/day in women; and < 800 or > 4000 kcal/day in men) . Participants with missing marital status data (n = 1506) in Fenland were retained for analyses and coded with a missing indicator variable.
In Fenland, a self-administered, 130-item, semi-quantitative food frequency questionnaire (FFQ) was used to assess habitual dietary intake over the previous year . The validity of the FFQ was previously assessed against 16-day weighed dietary record, 24-h recall, and biomarkers [25,26,27]. In this FFQ, the average consumption frequencies of each food item ranging from ‘never or less than once per month’ to ‘six or more times per day’ (nine categories) were provided. Total energy and macronutrient intake were estimated based on the UK food composition database .
In CoLaus, a self-administered, 97-item, semi-quantitative FFQ was used to assess dietary intake over the preceding 4 weeks , the validity of which had been assessed in nearby Canton of Geneva against 24-h recalls [30, 31]. For each item, consumption frequencies ranging from ‘less than once during the last 4 weeks’ to ‘two or more times per day’ (seven categories) were provided, in addition to the serving size (smaller, equal, or bigger) in comparison to a reference size.
Mediterranean diet scores
We calculated three Mediterranean diet scores (MDSs). We evaluated the pyramid-based MDS (pyrMDS)  based on the Mediterranean dietary pyramid  as our primary exposure. The Mediterranean dietary pyramid was proposed by the Mediterranean Diet Foundation to be applied to both Mediterranean and non-Mediterranean countries . We previously confirmed the content validity of pyrMDS in a non-Mediterranean setting, with its higher scores being associated with lower incidence of cardiovascular disease and cardiovascular and all-cause mortality in a UK population . The detailed scoring method has been reported previously . Briefly, a continuous score of 0 to 1 was assigned for each recommended consumption level of the 15 components of the pyramid (possible range 0–15): vegetables, legumes, and fish as healthy food groups; red meat, processed meat, potato, and sweets as unhealthy food groups; and fruits, cereals, nuts, eggs, dairy, white meat, and alcoholic beverages as items for which a moderate consumption was recommended (Additional file 1: Table S1). The second MDS was based on an algorithm proposed by Sofi and Casini  based on a systematic review of the published literature (literature-based MDS, LitMDS; possible range 0–18). The LitMDS accounted for 9 dietary items: vegetables, legumes, fruits and nuts, cereals, dairy, fish, meat, alcohol, and olive oil. The third MDS was based on each cohort’s sex-specific tertiles (tMDS; possible range 0–18) and accounted for the same 9 dietary items as LitMDS (Additional file 1: Table S1). We also previously tested these scores in a UK population . The MDS calculation was adjusted to an energy intake of 2000 kcal/day (8.37 MJ/day) based on the residual method [23, 32].
Ascertainment of hepatic steatosis
In Fenland, hepatic steatosis was ascertained by abdominal ultrasound, which is considered the first-line diagnostic procedure for hepatic steatosis . A semi-quantitative grading system was used to define normal hepatic echotexture or mild, moderate, and severe steatosis. The images were scored retrospectively according to standardised criteria by two trained operators who were unaware of other study measures. Hepatic steatosis scoring criteria were (a) increased echotexture of the liver parenchyma (bright liver in comparison with the kidney), (b) decreased visualisation of the intra-hepatic vasculature, and (c) attenuation of the ultrasound beam. Each criterion was scored on a 4-point scale, and a cumulative liver fat score based on the sum of the scores was created (possible range 3–12) . A score of ≤ 4 was classified as normal liver, 5–7 as mild steatosis, 8–10 as moderate steatosis, and ≥ 11 as severe steatosis. In our primary analysis, the mild, moderate, and severe categories were grouped together. The diagnostic accuracy of ultrasound was previously assessed against proton magnetic resonance spectroscopy, with sensitivity and specificity of 96% and 94%, respectively .
In both cohorts, we also evaluated fatty liver index (FLI) as an outcome, using anthropometry measures and fasting blood markers, calculated based on a logistic function including body mass index (BMI), waist circumference (WC), triglyceride (TG), and gamma-glutamyl transferase (GGT) levels as follows:
FLI × 100 ranged from 0 to 100, and the presence of hepatic steatosis was defined by FLI ≥ 60, a value with sensitivity and specificity of 61% and 86%, respectively . The diagnostic accuracy of FLI in comparison to ultrasonography has been reported to have an area under the receiver operating characteristic curve of 0.813 (95% CI 0.797, 0.830) .
In the CoLaus Study, abdominal ultrasound measurements were not available and hepatic steatosis was assessed by FLI (as above) and additionally by the ‘NAFLD liver fat score’ . The NAFLD liver fat score was based on an algorithm including the presence of metabolic syndrome defined by the International Diabetes Federation criteria , presence of type 2 diabetes, and fasting concentrations of insulin, aspartate-aminotransferase (AST), and the AST/alanine transaminase (ALT) ratio:
Compared to proton magnetic resonance imaging, the presence of hepatic steatosis defined by a NAFLD liver fat score greater than or equal to − 0.640 had a sensitivity of 86% and a specificity of 71% . In Fenland, AST levels were not available to calculate NAFLD liver fat score.
Assessment of covariates
In both studies, socio-demographic, lifestyle, and health characteristics were collected by self-administered questionnaires. Socio-demographic data included age, sex, marital status (single, married/cohabiting, and widowed/separated/divorced), occupational social class (routine/manual and administrative/professional in Fenland; working and not working in CoLaus), and educational level (compulsory, secondary, and university). Health characteristics included the presence of metabolic syndrome and family history of diabetes. In the Fenland Study, test site (Cambridge, Ely, and Wisbech) and household income (< £20,000, £20,000–40,000, and > £40,000) were also used as covariates. Smoking status was classified as ‘never’, ‘former’, and ‘current’.
In Fenland, physical activity was assessed objectively using combined heart rate and movement sensing for over 6 days (Actiheart, CamNTech, Cambridge, UK) with individual calibration for heart rate using a treadmill test . To estimate intensity time series, free-living data were pre-processed and modelled using a branched equation framework then summarised over time as daily physical activity energy expenditure (kcal/day). In CoLaus, physical activity was assessed with a validated self-administered quantitative physical activity frequency questionnaire [41, 42].
In both cohorts, body weight and height were measured with participants barefoot and in light indoor clothes; BMI was calculated as weight (kg) divided by height squared (m2). Waist circumference was measured by tape mid-way between the lowest rib and the iliac crest. In Fenland, body fat mass was also measured with dual-energy X-ray absorptiometry. Blood pressure was measured three times using an automated oscillometric sphygmomanometer in both cohorts, and the average of the two last measurements was used to define systolic and diastolic blood pressure.
Fasting venous blood samples were collected. In Fenland, blood samples were placed on ice, centrifuged, and stored at − 70 °C until analysis. In CoLaus, all assays were performed on the blood samples within 2 h of blood collection. In both cohorts, plasma TG, high-density lipoprotein cholesterol, and glucose were measured using standard enzymatic methods and ALT, AST (only in CoLaus), and GGT were measured using reference method as standardised by the International Federation of Clinical Chemistry.
In both cohorts, alcohol consumption was assessed by self-report number of alcoholic beverage units consumed in the preceding week, categorised as ‘abstainers’ (0 unit/week), ‘moderate’ (1–21 units/week for men, 1–14 for women), and ‘heavy’ (> 21 units/week for men, > 14 for women) drinkers. We undertook two approaches regarding alcohol consumption: one as a covariate and the second as a component of the MDS to evaluate adherence to the Mediterranean diet from which alcohol consumption was separated out as it is an established risk factor for hepatic steatosis .
Statistical analyses were performed using Stata (version 14; StataCorp, College Station, TX). Descriptive statistics were presented as mean ± standard deviation (SD) for continuous variables and proportions for categorical variables. Cohen’s kappa statistics were calculated to assess the agreement between the FLI and ultrasound liver fat score in Fenland and between the FLI and NAFLD liver fat score in CoLaus.
Each of the three MDSs was modelled both categorically (quintiles) and continuously (per SD). Multivariable-adjusted Poisson regression was used to estimate the prevalence ratios (PR) and 95% confidence intervals (CI)  and to examine the association between hepatic steatosis (ultrasound liver fat score and FLI in Fenland; FLI and NAFLD liver fat score in CoLaus) as dependent variables and the different MDSs (pyrMDS, litMDS, and tMDS) as independent variables.
Analyses were adjusted for potential confounders including age, sex, marital status, occupational status, educational level, smoking status, energy intake, and physical activity energy expenditure (Fenland) or estimated total energy expenditure (CoLaus). Further adjustments for BMI as potential confounder were also conducted. In Fenland, the FFQ aimed at assessing habitual dietary intake across the previous year, whereas in CoLaus, it was over the past 4 weeks. Hence, to adjust for possible seasonal variation in CoLaus, the dates of dietary intake assessment were included in regression models.
A priori, we examined whether the association between adherence to the pyrMDS and hepatic steatosis varied by alcohol consumption, testing for statistical interaction by alcohol consumption and adherence to the Mediterranean diet, and we also conducted analysis stratified by alcohol consumption (abstainers, moderate, and heavy drinkers). We also assessed the influence of adjustment for both BMI and waist circumference; for body fat mass (Fenland only); for alcohol consumption (units/week); for clinical variables, including blood pressure > 130/85 mmHg (yes/no), TG level > 1.7 mmol/L (yes/no), high-density lipoprotein level < 1.29 mmol/L for men and < 1.03 mmol/L for women (yes/no), and glucose level ≥ 5.6 mmol/L (yes/no); and for family history of diabetes and metabolic syndrome (except for NAFLD liver fat score as metabolic syndrome is one of its components).
We conducted sensitivity analyses excluding the alcohol component from the pyrMDS to rule out the possible impact of alcohol on the observed association; excluding participants with BMI ≥ 30 kg/m2; excluding participants with excessive alcohol consumption; including participants with probable implausible energy intake; excluding participants with probable secondary causes of hepatic steatosis such as hepatitis B or C, HIV, or hepatotoxic medications; or including participants with diabetes (only for NAFLD liver fat score in CoLaus). The British National Formulary codes were used to identify hepatitis B and C or HIV in the Fenland cohort and the Anatomical Therapeutic Chemical classification of the World Health Organization to identify hepatotoxic medications in the CoLaus Study. Finally, we assessed the robustness of the results when modelling ultrasound liver fat score (1) using 7 as the cut-point to define hepatic steatosis (normal/mild vs. moderate/severe), and (2) continuously, we also examined the association of adherence to the Mediterranean diet with ALT and GGT levels as crude markers of hepatic steatosis [1, 44]. The latter were natural log transformed prior to regression analysis.
Possible interactions between the different MDSs with age, sex, and BMI in the main model were examined using the Wald test. We prespecified stratified analyses if significant interaction was identified. Statistical significance was considered for a two-sided test with a p value < 0.05.
Characteristics of included participants
Of the 12,435 participants in the Fenland Study, 2790 (22.4%) were excluded, leaving 9645 participants (54.4% women; 48.9 ± 7.4 years) for analysis. Of the 5064 participants in CoLaus, 1107 (21.8%) were excluded, leaving 3957 participants (56.8% women; 57.0 ± 10.4 years) for analysis. The reasons for exclusion are summarised in Additional file 1: Figure S1, and the characteristics of included and excluded participants are shown in Additional file 1: Table S2. The mean pyrMDS score was 9.07 ± 1.43 and 8.45 ± 1.24 in Fenland and CoLaus, respectively. In both cohorts, adherence to the Mediterranean diet was higher among women compared to men, positively correlated with socio-economic status and negatively correlated with TG, liver enzyme levels, BMI, waist circumference, prevalence of current smokers, and prevalence of the metabolic syndrome (Table 1).
Adherence to the Mediterranean diet and prevalence of hepatic steatosis
In Fenland, the prevalence of hepatic steatosis was 23.9% and 27.1% based on ultrasound and FLI, respectively. Diagnoses using FLI and ultrasound liver fat score were concordant (kappa = 0.79). In CoLaus, the prevalence of hepatic steatosis was 25.3% and 25.7% based on FLI and NAFLD score, respectively, with high concordance (kappa = 0.83) between the measures. Hepatic steatosis was more prevalent among men; those with higher BMI, waist circumference, and liver enzymes; and those with metabolic syndrome (Additional file 1: Table S2).
In Fenland, the multivariable-adjusted analysis showed an inverse association between pyrMDS quintiles and hepatic steatosis based on ultrasound (ptrend < 0.001) (Table 2), but the association was attenuated after further adjustment for BMI (ptrend < 0.043). An inverse association was also seen per one SD difference in pyrMDS [PR of 0.86 (95% CI 0.81–0.90)] but was attenuated after adjustment for BMI [0.95 (0.90, 1.00)].
In both cohorts, when FLI was the outcome, there was an inverse association between pyrMDS quintiles and hepatic steatosis (ptrend < 0.001) (Table 3), but the association attenuated after adjustment for BMI (ptrend = 0.001 and 0.009 for Fenland and CoLaus, respectively). Significant inverse association was seen per one SD difference in the pyrMDS in both cohorts; the observed association was attenuated but remained significant in Fenland after adjustment for BMI (Table 3).
In CoLaus, when NAFLD liver fat score was the outcome, there was a non-significant inverse association between pyrMDS quintiles and hepatic steatosis (ptrend = 0.022). Using pyrMDS as a continuous variable, there was an inverse association per one SD difference [0.93 (0.87–0.99)], which attenuated after adjustment for BMI (Table 3).
Associations using LitMDS or tMDS as exposure variables were largely consistent with those of pyrMDS (Additional file 1: Table S3).
In both cohorts, analyses stratified by alcohol consumption showed that the inverse association between pyrMDS quintiles and hepatic steatosis was non-significant among ‘abstainers’ but significant among ‘moderate’ drinkers, except for NAFLD liver fat score in CoLaus. Among ‘heavy’ drinkers, the results were inconsistent across the two cohorts (Additional file 1: Table S4). Further adjustment for BMI attenuated the inverse associations between pyrMDS and hepatic steatosis within alcohol consumption strata.
Adjusting for both BMI and waist circumference attenuated the associations (ptrend > 0.05), except for FLI score in both studies (Additional file 1: Table S5). In Fenland, the inverse association between pyrMDS and hepatic steatosis as assessed by ultrasound remained significant after adjusting for body fat mass (Additional file 1: Table S5). Further adjustment for other potential confounders or mediators including alcohol consumption with and without adiposity measures, clinical variables, and presence of metabolic syndrome partly attenuated the associations (Additional file 1: Table S5). The results were not materially different when excluding the alcohol component from the pyrMDS; when excluding participants based on high BMI ≥ 30 kg/m2, excessive alcohol consumption, or secondary causes of hepatic steatosis; or when including participants with implausible energy intakes or with diabetes (Additional file 1: Table S5).
In Fenland, the use of a different dichotomisation of the ultrasound liver fat score weakened the association between pyrMDS quintiles and hepatic steatosis (Additional file 1: Table S6). The association per SD of pyrMDS was not altered, but CIs were widened (Additional file 1: Table S6). Associations modelling ultrasound liver fat score as a continuous outcome variable attenuated the inverse association, but it remained statistically significant (ptrend < 0.001) (Additional file 1: Table S6). Adherence to the Mediterranean diet was inversely but not significantly associated with ALT and GGT (Additional file 1: Table S7).
Finally, there were no significant interactions between adherence to the Mediterranean diet and age or sex in either cohort. A significant interaction between pyrMDS and BMI was found in both cohorts (pinteraction < 0.05), except when considering ultrasound liver fat score as an outcome in Fenland (pinteraction = 0.33); hence, due to inconsistent results, we conducted an analysis stratified by BMI categories (underweight/normal, overweight, and obese) (Additional file 1: Table S8).
In Fenland, the inverse association between pyrMDS quintiles and hepatic steatosis was non-significant (ptrend > 0.05) in the analyses stratified by BMI, except among overweight participants when considering FLI as an outcome (ptrend < 0.001). In CoLaus, considering FLI as an outcome, a non-significant inverse association between pyrMDS quintiles and hepatic steatosis was found among underweight/normal (ptrend = 0.006) and overweight participants (ptrend = 0.023). In contrast, considering NAFLD liver fat score as an outcome in CoLaus, stratified analyses showed no significant association between pyrMDS quintiles and hepatic steatosis (ptrend > 0.05).
In two independent population-based cohorts in the UK and Switzerland, we found an inverse association between greater adherence to the Mediterranean diet and the prevalence of hepatic steatosis. This association was consistent across three different definitions of adherence to the Mediterranean diet and with different non-invasive criteria for assessing hepatic steatosis but was attenuated when we included measures of adiposity as a covariate.
Comparison with other studies
Our findings add to the limited number of studies that evaluated the associations between single food groups or other dietary patterns and hepatic steatosis. A cross-sectional study previously reported that an eating pattern high in alcohol and meat intakes and with a low tea consumption was associated with higher liver fat content among 354 individuals from the PopGen cohort in Germany . Another cross-sectional study among Israeli adults reported a positive association between soft drink consumption and risk of NAFLD, whereas higher intake of omega-3 polyunsaturated fatty acids was negatively associated with the risk of NAFLD . Additionally, one Australian prospective study among adolescents indicated that a Western dietary pattern characterised by a high consumption of takeaway foods, confectionary, red or processed meat, refined grain, chips, full-fat dairy, and soft drinks at age 14 was associated with the presence of NAFLD at age 17 years . Our findings for the Mediterranean diet are consistent with a recent secondary cross-sectional analysis of the PREDIMED trial in Spain, although that study focused on overweight and obese participants at high cardiovascular risk, and assessed a dietary inflammatory index in conjunction with adherence to the Mediterranean diet . Moreover, our results are in agreement with the only prospective study to date reported a significant inverse association between adherence to the Mediterranean diet and incident hepatic steatosis among healthy individuals without hepatic steatosis .
Nevertheless, our findings of an inverse association between Mediterranean diet and hepatic steatosis were not consistent with those from other studies. Possible explanations may include differences in population and methods. For instance, one cross-sectional study of healthy Chinese adults in Hong Kong reported no significant association between higher adherence to the Mediterranean diet and NAFLD prevalence even in the analyses without BMI adjustment . This could be because the authors defined MDS based on population-specific medians of the study population, which might be too crude [19, 32]. Additionally, a small case-control study conducted in Athens compared 73 patients with NAFLD with 58 healthy controls and found no association between greater MDS and lower likelihood of having NAFLD .
In our study, the inverse association between adherence to the Mediterranean diet and hepatic steatosis was no longer significant after adjustment for adiposity (BMI and waist circumference) suggesting that adiposity may be a confounder or a mediator of the association. While our study was not designed to distinguish between these two phenomena, our findings were consistent with previous reports of non-significant associations between adherence to the Mediterranean diet and NAFLD after adjustment for BMI or abdominal fat [19, 47]. Our findings are also in agreement with a previous prospective study suggesting that the Western dietary pattern acts predominantly on hepatic steatosis incidence via the obesity pathway .
Previous studies have reported positive  or null  associations between the Mediterranean diet and liver enzymes. The null associations reported previously are in line with our results; the possible explanation could be the normal levels of liver enzymes in the majority of individuals with hepatic steatosis .
Possible mechanisms and implications
The potential mechanisms explaining the inverse association of the individual components of the Mediterranean diet on hepatic steatosis have been reviewed previously [51, 52]. Adherence to the Mediterranean diet is associated with high antioxidant capacity and anti-inflammatory and anti-fibrotic effects due to high fruit and vegetable consumption . The Mediterranean diet is also associated with high levels of monounsaturated fatty acids (MUFA), mostly from olive oil intake. In animal models, high MUFA from olive oil and/or PUFA improved lipid profile, decreased cytokine expression in visceral adipocytes, and decreased liver enzymes and hepatic TGs, thus attenuating hepatic steatosis [54,55,56]. These results have been corroborated in a meta-analysis of human studies, which confirmed that omega-3 fatty acids were negatively associated with hepatic steatosis .
Although epidemiological studies linking Mediterranean diet to hepatic steatosis are scarce, there is strong evidence regarding the inverse association between adherence to the Mediterranean diet and cardiovascular disease, the metabolic syndrome, and its components [9, 10, 58]. Few observational studies and randomised clinical trials have evaluated the impact of adherence to the Mediterranean diet on hepatic steatosis progression [12, 48, 49, 59]. All these studies, despite their difference in study design, sample population, outcome measure assessment methods, and various Mediterranean diet-scoring methods, reported the beneficial impact of greater adherence to the Mediterranean diet on reduced progression of hepatic steatosis. Thus, the Mediterranean diet may represent an alternative therapy for hepatic steatosis, given the well-recognised problems of achieving sustainable weight loss, which is currently the primary treatment for hepatic steatosis through lifestyle therapy involving diet and exercise. Therefore, the clinical importance of the Mediterranean diet for the prevention of hepatic steatosis deserves further discussion.
Strengths and limitations
Our study in two population-based cohorts in two independent settings of the UK and Switzerland, together with the large sample size of nearly 14,000 participants, allowed us to conduct stratified analysis with adequate statistical power. The application of three different scoring algorithms for the Mediterranean diet enhanced the utility of representing adherence to the diet in different epidemiological settings. Finally, the consistency of the findings in two independent cohorts suggests that the results may be reasonably generalisable to other populations.
Our study also has some limitations. First, its cross-sectional design precludes any causal interpretation of the association between adherence to the Mediterranean diet and hepatic steatosis. The study design was also limited by its inability to distinguish whether adiposity mediated or confounded the observed association. Second, hepatic steatosis was based on different indicators and not on liver biopsy, which is considered to be the gold standard technique. However, liver biopsy is an invasive method, which cannot be implemented in large-scale epidemiological studies with predominantly healthy participants. In the Fenland Study, hepatic steatosis was based on ultrasound, a valid method for detecting hepatic steatosis , supplemented with an algorithm-based FLI, while in the CoLaus Study, two hepatic steatosis indices were used. These validated markers of hepatic steatosis are based on non-invasive and easily ascertained measurements which make them suitable for large epidemiological studies [38, 39, 60, 61], and the results were consistent irrespective of the marker used. Third, excluding participants with diabetes might have led to an underestimation of hepatic steatosis prevalence as assessed by NAFLD liver fat score and its association with the Mediterranean diet. However, findings were similar upon including participants with diabetes in the sensitivity analysis. Fourth, the methods for measuring variables (e.g. physical activity) differed between cohorts; however, the consistency of our findings suggests that the impact of the differing methods on results was minimal. Fifth, self-reported dietary intake and covariates are prone to measurement errors but the instruments used were validated against reference methods, and we found consistency of our estimates across the levels of adjustments.
Greater adherence to the Mediterranean diet was associated with lower likelihood of hepatic steatosis, and this association was largely related to markers of adiposity. These cross-sectional findings were strengthened by their consistency in two independent settings in the UK and Switzerland and are hypothesis generating for the role of the Mediterranean diet in hepatic steatosis in non-Mediterranean settings. The cross-sectional findings from our population-based cohorts warrant further interventional or observational studies to test whether the improvement of adherence to the Mediterranean diet is efficacious and effective in diverse settings for the primary prevention of hepatic steatosis.
Body mass index
Food frequency questionnaire
Fatty liver index
Literature-based Mediterranean diet scores
Mediterranean diet scores
Monounsaturated fatty acids
Non-alcoholic fatty liver disease
Polyunsaturated fatty acids
Pyramid-based Mediterranean diet scores
Tertile-based Mediterranean diet scores
Sattar N, Forrest E, Preiss D. Non-alcoholic fatty liver disease. BMJ. 2014;349:g4596.
Preiss D, Sattar N. Non-alcoholic fatty liver disease: an overview of prevalence, diagnosis, pathogenesis and treatment considerations. Clin Sci (Lond). 2008;115(5):141–50.
EASL-EASD-EASO Clinical Practice Guidelines for the management of non-alcoholic fatty liver disease. J Hepatol. 2016;64(6):1388-1402. https://doi.org/10.1016/j.jhep.2015.11.004.
Abenavoli L, Milic N, Di Renzo L, Preveden T, Medic-Stojanoska M, De Lorenzo A. Metabolic aspects of adult patients with nonalcoholic fatty liver disease. World J Gastroenterol. 2016;22(31):7006–16.
Zelber-Sagi S, Nitzan-Kaluski D, Goldsmith R, Webb M, Blendis L, Halpern Z, Oren R. Long term nutritional intake and the risk for non-alcoholic fatty liver disease (NAFLD): a population based study. J Hepatol. 2007;47(5):711–7.
Ricci G, Canducci E, Pasini V, Rossi A, Bersani G, Ricci E, Alvisi V. Nutrient intake in Italian obese patients: relationships with insulin resistance and markers of non-alcoholic fatty liver disease. Nutrition. 2011;27(6):672–6.
Kim CH, Kallman JB, Bai C, Pawloski L, Gewa C, Arsalla A, Sabatella ME, Younossi ZM. Nutritional assessments of patients with non-alcoholic fatty liver disease. Obes Surg. 2010;20(2):154–60.
Musso G, Gambino R, De Michieli F, Cassader M, Rizzetto M, Durazzo M, Faga E, Silli B, Pagano G. Dietary habits and their relations to insulin resistance and postprandial lipemia in nonalcoholic steatohepatitis. Hepatology. 2003;37(4):909–16.
Rees K, Hartley L, Flowers N, Clarke A, Hooper L, Thorogood M, Stranges S. 'Mediterranean' dietary pattern for the primary prevention of cardiovascular disease. Cochrane Database Syst Rev. 2013;(8). https://doi.org/10.1002/14651858.CD009825.pub2.
Kastorini CM, Milionis HJ, Esposito K, Giugliano D, Goudevenos JA, Panagiotakos DB. The effect of Mediterranean diet on metabolic syndrome and its components: a meta-analysis of 50 studies and 534,906 individuals. J Am Coll Cardiol. 2011;57(11):1299–313.
Sofi F, Cesari F, Abbate R, Gensini GF, Casini A. Adherence to Mediterranean diet and health status: meta-analysis. BMJ. 2008;337:a1344.
Ryan MC, Itsiopoulos C, Thodis T, Ward G, Trost N, Hofferberth S, O’Dea K, Desmond PV, Johnson NA, Wilson AM. The Mediterranean diet improves hepatic steatosis and insulin sensitivity in individuals with non-alcoholic fatty liver disease. J Hepatol. 2013;59(1):138–43.
Abenavoli L, Greco M, Nazionale I, Peta V, Milic N, Accattato F, Foti D, Gulletta E, Luzza F. Effects of Mediterranean diet supplemented with silybin-vitamin E-phospholipid complex in overweight patients with non-alcoholic fatty liver disease. Expert Rev Gastroenterol Hepatol. 2015;9(4):519–27.
Trovato FM, Catalano D, Martines GF, Pace P, Trovato GM. Mediterranean diet and non-alcoholic fatty liver disease: the need of extended and comprehensive interventions. Clin Nutr. 2015;34(1):86–8.
Gelli C, Tarocchi M, Abenavoli L, Di Renzo L, Galli A, De Lorenzo A. Effect of a counseling-supported treatment with the Mediterranean diet and physical activity on the severity of the non-alcoholic fatty liver disease. World J Gastroenterol. 2017;23(17):3150–62.
Katsagoni CN, Papatheodoridis GV, Ioannidou P, Deutsch M, Alexopoulou A, Papadopoulos N, Papageorgiou M-V, Fragopoulou E, Kontogianni MD. Improvements in clinical characteristics of patients with non-alcoholic fatty liver disease, after an intervention based on the Mediterranean lifestyle: a randomised controlled clinical trial. Br J Nutr. 2018;120(2):164–75.
Properzi C, O'Sullivan TA, Sherriff JL, Ching HL, Jeffrey GP, Buckley RF, Tibballs J, MacQuillan GC, Garas G, Adams LA. Ad libitum Mediterranean and Low Fat Diets both Significantly Reduce Hepatic Steatosis: a Randomized Controlled Trial. Hepatology. 2018;68(5):1741-1754.
Cantero I, Abete I, Babio N, Aros F, Corella D, Estruch R, Fito M, Hebert JR, Martinez-Gonzalez MA, Pinto X, et al. Dietary Inflammatory Index and liver status in subjects with different adiposity levels within the PREDIMED trial. Clin Nutr. 2018;37(5):1736-1743.
Chan R, Wong VW, Chu WC, Wong GL, Li LS, Leung J, Chim AM, Yeung DK, Sea MM, Woo J, et al. Diet-quality scores and prevalence of nonalcoholic fatty liver disease: a population study using proton-magnetic resonance spectroscopy. PLoS One. 2015;10(9):e0139310.
Ma J, Hennein R, Liu C, Long MT, Hoffmann U, Jacques PF, Lichtenstein AH, Hu FB, Levy D. Improved diet quality associates with reduction in liver fat, particularly in individuals with high genetic risk scores for nonalcoholic fatty liver disease. Gastroenterology. 2018;155(1):107–17.
Rolfe Ede L, Loos RJ, Druet C, Stolk RP, Ekelund U, Griffin SJ, Forouhi NG, Wareham NJ, Ong KK. Association between birth weight and visceral fat in adults. Am J Clin Nutr. 2010;92(2):347–52.
Firmann M, Mayor V, Vidal PM, Bochud M, Pecoud A, Hayoz D, Paccaud F, Preisig M, Song KS, Yuan X, et al. The CoLaus Study: a population-based study to investigate the epidemiology and genetic determinants of cardiovascular risk factors and metabolic syndrome. BMC Cardiovasc Disord. 2008;8:6.
Willett WC, Howe GR, Kushi LH. Adjustment for total energy intake in epidemiologic studies. Am J Clin Nutr. 1997;65(4 Suppl):1220S–8S discussion 1229S-1231S.
Bingham SA, Welch AA, McTaggart A, Mulligan AA, Runswick SA, Luben R, Oakes S, Khaw KT, Wareham N, Day NE. Nutritional methods in the European Prospective Investigation of Cancer in Norfolk. Public Health Nutr. 2001;4(3):847–58.
Bingham SA, Gill C, Welch A, Day K, Cassidy A, Khaw KT, Sneyd MJ, Key TJ, Roe L, Day NE. Comparison of dietary assessment methods in nutritional epidemiology: weighed records v. 24 h recalls, food-frequency questionnaires and estimated-diet records. Br J Nutr. 1994;72(4):619–43.
Bingham SA, Gill C, Welch A, Cassidy A, Runswick SA, Oakes S, Lubin R, Thurnham DI, Key TJ, Roe L, et al. Validation of dietary assessment methods in the UK arm of EPIC using weighed records, and 24-hour urinary nitrogen and potassium and serum vitamin C and carotenoids as biomarkers. Int J Epidemiol. 1997;26(Suppl 1):S137–51.
Bingham SA, Cassidy A, Cole TJ, Welch A, Runswick SA, Black AE, Thurnham D, Bates C, Khaw KT, Key TJ, et al. Validation of weighed records and other methods of dietary assessment using the 24 h urine nitrogen technique and other biological markers. Br J Nutr. 1995;73(4):531–50.
Mulligan AA, Luben RN, Bhaniani A, Parry-Smith DJ, O’Connor L, Khawaja AP, Forouhi NG, Khaw KT. A new tool for converting food frequency questionnaire data into nutrient and food group values: FETA research methods and availability. BMJ Open. 2014;4(3):e004503.
Morabia A, Bernstein M, Kumanyika S, Sorenson A, Mabiala I, Prodolliet B, Rolfo I, Luong BL. Développement et validation d’un questionnaire alimentaire semi-quantitatif à partir d’une enquête de population. Soz Praventivmed. 1994;39(6):345–69.
Bernstein L, Huot I, Morabia A. Amélioration des performances d’un questionnaire alimentaire semi-quantitatif comparé à un rappel des 24 heures. Sante Publique (Bucur). 1995;7:403–13.
Beer-Borst S, Costanza MC, Pechere-Bertschi A, Morabia A. Twelve-year trends and correlates of dietary salt intakes for the general adult population of Geneva, Switzerland. Eur J Clin Nutr. 2009;63(2):155–64.
Tong TY, Wareham NJ, Khaw KT, Imamura F, Forouhi NG. Prospective association of the Mediterranean diet with cardiovascular disease incidence and mortality and its population impact in a non-Mediterranean population: the EPIC-Norfolk study. BMC Med. 2016;14(1):135.
Bach-Faig A, Berry EM, Lairon D, Reguant J, Trichopoulou A, Dernini S, Medina FX, Battino M, Belahsen R, Miranda G, et al. Mediterranean diet pyramid today. Science and cultural updates. Public Health Nutr. 2011;14(12a):2274–84.
Sofi F, Casini A. Mediterranean diet and non-alcoholic fatty liver disease: new therapeutic option around the corner? World J Gastroenterol. 2014;20(23):7339–46.
Mehta SR, Thomas EL, Bell JD, Johnston DG, Taylor-Robinson SD. Non-invasive means of measuring hepatic fat content. World J Gastroenterol. 2008;14(22):3476–83.
De Lucia Rolfe E, Brage S, Sleigh A, Finucane F, Griffin SJ, Wareham NJ, Ong KK, Forouhi NG. Validity of ultrasonography to assess hepatic steatosis compared to magnetic resonance spectroscopy as a criterion method in older adults. PLoS One. 2018;13(11):e0207923.
Bedogni G, Bellentani S, Miglioli L, Masutti F, Passalacqua M, Castiglione A, Tiribelli C. The fatty liver index: a simple and accurate predictor of hepatic steatosis in the general population. BMC Gastroenterol. 2006;6:33.
Koehler EM, Schouten JN, Hansen BE, Hofman A, Stricker BH, Janssen HL. External validation of the fatty liver index for identifying nonalcoholic fatty liver disease in a population-based study. Clin Gastroenterol Hepatol. 2013;11(9):1201–4.
Kotronen A, Peltonen M, Hakkarainen A, Sevastianova K, Bergholm R, Johansson LM, Lundbom N, Rissanen A, Ridderstrale M, Groop L, et al. Prediction of non-alcoholic fatty liver disease and liver fat using metabolic and genetic factors. Gastroenterology. 2009;137(3):865–72.
A new international diabetes federation worldwide definition of the metabolic syndrome: the rationale and the results. Rev Esp Cardiol 2005, 58(12):1371-1375.
Brage S, Ekelund U, Brage N, Hennings MA, Froberg K, Franks PW, Wareham NJ. Hierarchy of individual calibration levels for heart rate and accelerometry to measure physical activity. J Appl Physiol. 2007;103(2):682–92.
Bernstein M, Sloutskis D, Kumanyika S, Sparti A, Schutz Y, Morabia A. Data-based approach for developing a physical activity frequency questionnaire. Am J Epidemiol. 1998;147(2):147–54.
Spiegelman D, Hertzmark E. Easy SAS calculations for risk or prevalence ratios and differences. Am J Epidemiol. 2005;162(3):199–200.
Rinella ME. Nonalcoholic fatty liver disease: a systematic review. JAMA. 2015;313(22):2263–73.
Koch M, Borggrefe J, Schlesinger S, Barbaresko J, Groth G, Jacobs G, Lieb W, Laudes M, Muller MJ, Bosy-Westphal A, et al. Association of a lifestyle index with MRI-determined liver fat content in a general population study. J Epidemiol Community Health. 2015;69(8):732–7.
Oddy WH, Herbison CE, Jacoby P, Ambrosini GL, O’Sullivan TA, Ayonrinde OT, Olynyk JK, Black LJ, Beilin LJ, Mori TA, et al. The Western dietary pattern is prospectively associated with nonalcoholic fatty liver disease in adolescence. Am J Gastroenterol. 2013;108(5):778–85.
Kontogianni MD, Tileli N, Margariti A, Georgoulis M, Deutsch M, Tiniakos D, Fragopoulou E, Zafiropoulou R, Manios Y, Papatheodoridis G. Adherence to the Mediterranean diet is associated with the severity of non-alcoholic fatty liver disease. Clin Nutr. 2014;33(4):678–83.
Abenavoli L, Greco M, Milic N, Accattato F, Foti D, Gulletta E, Luzza F. Effect of Mediterranean Diet and Antioxidant Formulation in Non-Alcoholic Fatty Liver Disease: A Randomized Study. Nutrients. 2017;9(8). https://doi.org/10.3390/nu9080870.
Fraser A, Abel R, Lawlor DA, Fraser D, Elhayany A. A modified Mediterranean diet is associated with the greatest reduction in alanine aminotransferase levels in obese type 2 diabetes patients: results of a quasi-randomised controlled trial. Diabetologia. 2008;51(9):1616–22.
Clark JM, Brancati FL, Diehl AM. The prevalence and etiology of elevated aminotransferase levels in the United States. Am J Gastroenterol. 2003;98(5):960–7.
Zelber-Sagi S, Salomone F, Mlynarsky L. The Mediterranean dietary pattern as the diet of choice for non-alcoholic fatty liver disease: evidence and plausible mechanisms. Liver Int. 2017;37(7):936–49.
Mahady SE, George J. Exercise and diet in the management of nonalcoholic fatty liver disease. Metabolism. 2016;65(8):1172–82.
Salomone F, Godos J, Zelber-Sagi S. Natural antioxidants for non-alcoholic fatty liver disease: molecular targets and clinical perspectives. Liver Int. 2016;36(1):5–20.
Jurado-Ruiz E, Varela LM, Luque A, Berna G, Cahuana G, Martinez-Force E, Gallego-Duran R, Soria B, de Roos B, Romero Gomez M, et al. An extra virgin olive oil rich diet intervention ameliorates the nonalcoholic steatohepatitis induced by a high-fat "Western-type" diet in mice. Mol Nutr Food Res. 2017;61(3):1600549.
Depner CM, Torres-Gonzalez M, Tripathy S, Milne G, Jump DB. Menhaden oil decreases high-fat diet-induced markers of hepatic damage, steatosis, inflammation, and fibrosis in obese Ldlr-/- mice. J Nutr. 2012;142(8):1495–503.
Byrne CD. Fatty liver: role of inflammation and fatty acid nutrition. Prostaglandins Leukot Essent Fatty Acids. 2010;82(4-6):265–71.
Parker HM, Johnson NA, Burdon CA, Cohn JS, O’Connor HT, George J. Omega-3 supplementation and non-alcoholic fatty liver disease: a systematic review and meta-analysis. J Hepatol. 2012;56(4):944–51.
Estruch R, Ros E, Salas-Salvado J, Covas MI, Corella D, Aros F, Gomez-Gracia E, Ruiz-Gutierrez 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(25):e34.
Aller R, Izaola O, de la Fuente B, De Luis Roman DA. Mediterranean diet is associated with liver histology in patients with nonalcoholic fatty liver disease. Nutr Hosp. 2015;32(6):2518–24.
Ruhl CE, Everhart JE. Fatty liver indices in the multiethnic United States National Health and Nutrition Examination Survey. Aliment Pharmacol Ther. 2015;41(1):65–76.
Machado MV, Cortez-Pinto H. Non-invasive diagnosis of non-alcoholic fatty liver disease. A critical appraisal. J Hepatol. 2013;58(5):1007–19.
The authors are grateful to all the participants, General Practitioners, and staff of Fenland and CoLaus studies (http://www.mrc-epid.cam.ac.uk/research/studies/fenland/ and https://www.colaus-psycolaus.ch/). The authors thank the Fenland Study co-ordination team, the Field Epidemiology team, the data management team and the laboratory team of the MRC Epidemiology Unit, and are grateful to Dr Sara Upponi at Addenbrooke’s Hospital, Cambridge, for reviewing and advising on ultrasound scans with potential abnormalities.
The Fenland Study is funded by the Medical Research Council Epidemiology Unit (SJAG/004 RG71153) and the Wellcome Trust. We acknowledge support from Medical Research Council programmes MC_UU_12015/1 and MC_UU_12015/5 and from NIHR Biomedical Research Centre Cambridge: Nutrition, Diet, and Lifestyle Research Theme (IS-BRC-1215-20014). The CoLaus study is supported by research grants from GlaxoSmithKline, the Faculty of Biology and Medicine of Lausanne, and the Swiss National Science Foundation (grants 33CSCO-122661, 33CS30-139468, and 33CS30-148401). SKS supported by the Swiss National Science Foundation (Doc.Mobility number P1LAP3-171805). NJW is an NIHR Senior Investigator.
Availability of data and materials
For both Fenland and CoLaus studies, non-identifiable individual-level data are available for researchers who seek to answer questions related to health and disease in the context of research projects who meet the criteria for data sharing by research committees. Please follow the instructions at http://www.mrc-epid.cam.ac.uk/research/data-sharing/ and https://www.colaus-psycolaus.ch/ for information on how to submit an application for gaining access to Fenland and CoLaus data, respectively.
Ethics approval and consent to participate
The Fenland Study was approved by the Cambridge Local Research Ethics Committee (NRES Committee East of England-Cambridge Central – 04/Q0108/19). The institutional Ethics Committee of the University of Lausanne, which later became the Ethics Commission of Canton Vaud (www.cer-vd.ch) approved the CoLaus study (reference 16/03). Both studies were performed in agreement with the Helsinki declaration and its former amendments. Informed consent was obtained from all participants in both studies.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Flow chart for the Fenland and CoLaus sample selection. Table S1. Mediterranean diet score specified by three definitions and their components. Table S2. Characteristics of participants included and excluded from the analysis, Fenland and CoLaus studies. Table S3. Association between adherence to the literature-based and tertile-based Mediterranean diet and prevalence of hepatic steatosis, Fenland and CoLaus studies. Table S4. Association between adherence to the Mediterranean diet and prevalence of hepatic steatosis within alcohol consumption strata, Fenland and CoLaus studies. Table S5. Sensitivity analyses for the association between adherence to the Mediterranean diet and prevalence of hepatic steatosis, Fenland and CoLaus studies. Table S6. Association between adherence to the Mediterranean diet and prevalence of hepatic steatosis, Fenland Study. Table S7. Association between adherence to the Mediterranean diet and ALT and GGT, Fenland and CoLaus studies. Table S8. Association between adherence to the Mediterranean diet and prevalence of hepatic steatosis within BMI strata, Fenland and CoLaus studies. (DOCX 124 kb)
About this article
Cite this article
Khalatbari-Soltani, S., Imamura, F., Brage, S. et al. The association between adherence to the Mediterranean diet and hepatic steatosis: cross-sectional analysis of two independent studies, the UK Fenland Study and the Swiss CoLaus Study. BMC Med 17, 19 (2019). https://doi.org/10.1186/s12916-019-1251-7