Hostname: page-component-7c8c6479df-nwzlb Total loading time: 0 Render date: 2024-03-29T05:51:41.772Z Has data issue: false hasContentIssue false

The association between Mediterranean Diet Score and glucokinase regulatory protein gene variation on the markers of cardiometabolic risk: an analysis in the European Prospective Investigation into Cancer (EPIC)-Norfolk study

Published online by Cambridge University Press:  07 May 2014

Mercedes Sotos-Prieto
Affiliation:
Genetic and Molecular Epidemiology Unit, Department of Preventive Medicine and Public Health, School of Medicine, University of Valencia, Valencia, Spain. CIBER Fisiopatología de la Obesidad y Nutrición, Valencia, Spain Department of Epidemiology, Atherothrombosis and Cardiovascular Imaging, Centro Nacional de Investigaciones Cardiovasculares (CNIC), 28029Madrid, Spain MRC Epidemiology Unit, University of Cambridge School of Clinical Medicine, Institute of Metabolic Science, Cambridge Biomedical Campus, CambridgeCB2 0QQ, UK
Robert Luben
Affiliation:
Department of Public Health and Primary Care, University of Cambridge, Cambridge, UK
Kay-Tee Khaw
Affiliation:
Department of Public Health and Primary Care, University of Cambridge, Cambridge, UK
Nicholas J. Wareham
Affiliation:
MRC Epidemiology Unit, University of Cambridge School of Clinical Medicine, Institute of Metabolic Science, Cambridge Biomedical Campus, CambridgeCB2 0QQ, UK
Nita G. Forouhi*
Affiliation:
MRC Epidemiology Unit, University of Cambridge School of Clinical Medicine, Institute of Metabolic Science, Cambridge Biomedical Campus, CambridgeCB2 0QQ, UK
*
*Corresponding author: Dr N. G. Forouhi, fax +44 1223 330316, email nita.forouhi@mrc-epid.cam.ac.uk
Rights & Permissions [Opens in a new window]

Abstract

Consumption of a Mediterranean diet (MD) and genetic variation in the glucokinase regulatory protein (GCKR) gene have been reported to be associated with TAG and glucose metabolism. It is uncertain whether there is any interaction between these factors. Therefore, the aims of the present study were to test the association of adherence to a MD and rs780094 (G>A) SNP in the GCKR gene with the markers of cardiometabolic risk, and to investigate the interaction between genetic variation and MD adherence. We studied 20 986 individuals from the European Prospective Investigation into Cancer (EPIC)-Norfolk study. The relative Mediterranean Diet Score (rMED: range 0–18) was used to assess MD adherence. Linear regression was used to estimate the association between the rMED, genotype and cardiometabolic continuous traits, adjusting for potential confounders. In adjusted analyses, we observed independent associations of MD adherence and genotype with cardiometabolic risk, with the highest risk group (AA genotype; lowest rMED) having higher concentrations of TAG, total cholesterol and apoB (12·5, 2·3 and 3·1 %, respectively) v. those at the lowest risk (GG genotype; highest rMED). However, the associations of MD adherence with metabolic markers did not differ by genotype, with no significant gene–diet interactions for lipids or for glycated Hb. In conclusion, we found independent associations of the rMED and of the GCKR genotype with cardiometabolic profile, but found no evidence of interaction between them.

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

The Mediterranean diet (MD) pattern has been associated with reduced CHD, cancer and overall mortality( Reference Trichopoulou, Costacou and Bamia 1 ). The MD is characterised by a relatively greater proportion of quantity and diversity of plant-derived foods (whole-grain cereals, raw and cooked vegetables, fresh and dried fruits, legumes and nuts), fish, a relatively moderate intake of meat and dairy products, with olive oil as the added fat, and a moderate intake of wine during meals( Reference Trichopoulou, Costacou and Bamia 1 ). Although different levels of adherence to the MD have been shown by country (with lower adherence in Northern Europe, medium in Central Europe and higher in Southern Europe among the countries included in the European Prospective Investigation into Cancer (EPIC)-InterAct study( Reference Romaguera, Guevara and Norat 2 )), similar associations with CVD were observed across populations in different regions( Reference Truthmann, Richter and Thiele 3 ). Despite some variation in findings( Reference Panagiotakos, Tzima and Pitsavos 4 Reference Goulet, Lamarche and Nadeau 7 ), the overall evidence from epidemiological and interventional studies suggests a beneficial effect of adherence to the MD on continuous metabolic traits of cardiovascular risk such as lipids and glycaemia( Reference Panagiotakos, Tzima and Pitsavos 4 , Reference Estruch, Martínez-González and Corella 8 Reference Tzima, Pitsavos and Panagiotakos 11 ).

Recently, the potential importance of apolipoproteins in the prediction of cardiovascular risk has been appraised( Reference Di Angelantonio, Sarwar and Perry 12 , Reference Di Angelantonio, Gao and Pennells 13 ), and an interventional study has recently reported an improvement in apolipoprotein profile (higher apoA-1, and lower apoB and apoB:A-1 ratio) among individuals in the MD intervention group compared with a control diet-supplemented group( Reference Solá, Fitó and Estruch 14 ). In addition to the potential beneficial effects of diet, it is now also acknowledged that genetic factors may influence the regulation of plasma lipids and glucose metabolism( Reference Corella and Ordovas 15 , Reference Willer, Sanna and Jackson 16 ), and that gene–diet interactions may help explain the variation in metabolic risk( Reference Corella and Ordovas 15 ). Genome-wide association studies have identified a number of SNP that may be associated with these traits( Reference Willer, Sanna and Jackson 16 , Reference Saxena, Voight and Lyssenko 17 ).

The rs780094 (G>A) SNP in the glucokinase regulatory protein (GCKR) gene has been associated with higher TAG concentrations( Reference Willer, Sanna and Jackson 16 Reference Sparsø, Andersen and Nielsen 20 ). Conversely, the A allele has been associated with lower insulin levels, lower glucose levels and lower diabetes prevalence( Reference Saxena, Voight and Lyssenko 17 , Reference Orho-Melander, Melander and Guiducci 18 , Reference Dupuis, Langenberg and Prokopenko 21 , Reference Vaxillaire, Cavalcanti-Proença and Dechaume 22 ). In a high cardiovascular risk Spanish population (PREDIMED-Valencia (PREvencion con DIeta MEDiterranea-Valencia) study), we have recently reported that a MD modulated the association of GCKR gene variation on TAG concentrations( Reference Sotos-Prieto, Guillén and Vicente Sorli 23 ). However, it is currently unknown whether such effects may apply in the general British population, and whether they apply to other lipids, apolipoproteins and glycaemic markers, and whether there is an interaction between genetic and dietary factors on these metabolic parameters.

Therefore, the objectives of the present study were (1) to assess the association between adherence to the MD and continuous metabolic traits related to lipids and glucose metabolism in a large population-based study, (2) to investigate the association of GCKR rs780094 (G>A) with TAG concentrations and other metabolic markers and (3) to examine the joint effects and interaction of MD adherence and genetic effects on the metabolic markers of lipids and glycaemia.

Subjects and methods

Study participants and design

The EPIC-Norfolk Study recruited 25 639 men and women, aged 40–79 years at baseline (1993–7), who were resident in and around Norwich, England. The present study has been described in detail previously( Reference Day, Oakes and Luben 24 ), and it was conducted according to the guidelines laid down in the Declaration of Helsinki, and all procedures involving human participants were approved by the Norfolk District Health Authority Ethics Committee. Written informed consent was obtained from all participants. Since the baseline health-check visit, there were three follow-up assessments including two postal questionnaires and a repeat health-check visit, but this cross-sectional analysis is based on the baseline visit. Health and lifestyle information was collected using a baseline questionnaire, which asked about the participants' personal and family health, demography, lifestyle (including diet and physical activity) and social status (education and occupation). A standardised health check was performed by trained nurses, including measurement of height (cm), weight (kg) and waist circumference (WC; cm) as described previously( Reference Day, Oakes and Luben 24 ). Non-fasting blood samples were collected. A detailed description of the storage method has been described previously( Reference Riboli 25 ). Concentrations of serum lipids, total cholesterol (TC), HDL-cholesterol (HDL-C) and TAG concentrations were measured on an RA-1000 (Bayer Diagnostics). LDL-cholesterol (LDL-C) was calculated using the Friedewald formula( Reference Friedewald, Levy and Fredrickson 26 ); when serum TAG concentration exceeded 4·0 mmol/l, LDL-C concentration was not calculated. Concentrations of serum apoA-1 and apoB were measured using an Olympus AU640 Analyser (Olympus UK Limited). Measurement of glycated Hb (HbA1c) was added halfway through the baseline visit in 1995, and was available in approximately half of the cohort (n 10 780, with n 10 746 with HbA1c levels < 6·5 %). The level of HbA1c was measured with HPLC on a Bio-Rad Diamat (Bio-Rad).

We excluded the participants with missing baseline FFQ data or with ten missing FFQ lines as well as those with unavailable information on lipid measures. Thus, 20 986 participants from the baseline visit were eligible for inclusion in the present analysis.

Dietary assessment and Mediterranean Diet Score

Participants completed a validated 130-item semi-quantitative FFQ about their habitual diet in the past year( Reference Bingham, Gill and Welch 27 ). For all food items, respondents were asked to report the frequency of consumption on a 9-point scale for a ‘medium serving’ from ‘never or once per month’ to ‘more than six times per day’. Amounts of energy and individual nutrient intake were calculated from the frequency and amount (medium serving size) of each food reported in the FFQ and converted into g/d by using in-house software, the Compositional Analyses from Frequency Estimates (CAFE) program( Reference Welch, Luben and Khaw 28 ). The EPIC FFQ was validated against 16 d weighed food records (n 127)( Reference Kaaks and Riboli 29 ).

Adherence to the MD was assessed by using the relative Mediterranean Diet Score (rMED)( Reference Buckland, González and Agudo 30 , Reference Buckland, Agudo and Lujan 31 ), which is a variation of the original MD Score( Reference Trichopoulou, Costacou and Bamia 1 ). This variation consists mainly in the score criterion and in the olive oil component. Similarly to the MD Score, the rMED included nine nutritional components that are characteristic of the MD and included some presumed ‘beneficial’ components (vegetables, legumes, fruit and nuts, cereals, fish and seafood, olive oil (including olive oil consumption instead of the ratio of monounsaturated fat:saturated fat as in the original MD) and moderate alcohol consumption) and other presumed ‘detrimental’ components (meat and meat products and dairy products), as described previously( Reference Romaguera, Guevara and Norat 2 , Reference Buckland, González and Agudo 30 , Reference Buckland, Agudo and Lujan 31 ). While the original MD Score criteria assigned 1 point for intakes above the median among the study participants and 0 point for intakes below the median, with reversed scoring for meat and dairy intakes (score range 0–9), in the rMED, each component (apart from alcohol consumption) was calculated as a function of energy density (as g/4184 kJ (1000 kcal)), and was then divided into tertiles of intakes. We assigned a value of 0, 1 or 2 to the first, second and third tertiles, respectively, of the intakes of vegetables, legumes, fruit and nuts, cereals, and fish and seafood and positively scoring higher intakes for the beneficial components. The scoring was reversed for the two presumed detrimental components (meat and meat products and dairy products) by assigning a higher score for lower intakes. The scoring for olive oil was modified for the rMED( Reference Buckland, Agudo and Lujan 31 ) because of the relatively large number of non-consumers. Therefore, 0 was assigned to non-consumers, 1 was assigned to participants with an intake below the median olive oil consumption (calculated only within olive oil consumers) and 2 was assigned to participants whose intake was equal to or above this median. For alcohol, a value of 2 was given to men with moderate alcohol consumption (intakes from 10 to 50 g/d), and a value of 0 was assigned otherwise, whereas for women, the corresponding cut-off points were 5 and 25 g/d. Therefore, the rMED ranged from 0 (indicating the lowest adherence to the MD) to 18 (indicating the highest adherence to the MD), which is different from the original MD Score that ranged between 0 and 9. Finally, the rMED was further classified into categories to reflect low (0–6 points), medium (7–10 points) or high (11–18 points) adherence to the MD on the basis of previously published cut-off points( Reference Buckland, González and Agudo 30 ).

Determination of genotypes

DNA for genotyping was extracted from EDTA whole-blood aliquots collected at the first and second health checks using phenol–chloroform extraction. The call rate was 99·5 %. There were 19 800 participants with available data on genotype (rs780094). We also included data on rs1260326 SNP (in high linkage disequilibrium with the rs780094 SNP; r 2 0·94) in a subset (n 7273), to replicate findings. Genotype distributions did not deviate from Hardy–Weinberg expectations (P= 0·14, rs780094 and P= 0·80, rs1260326).

Statistical analysis

Continuous variables were assessed for normality of distribution, and skewed variables were normalised by log10 transformation. Data are presented as means and standard deviations or geometric means and 95 % CI for continuous variables, and frequencies and percentages for categorical variables. Baseline characteristics were described by rMED categories (low, medium and high). Differences in means among the groups were compared by using a one-factor ANOVA test, and comparisons of frequencies were conducted using the χ2 test. We fitted multiple linear regression models for the association between the rMED categories and continuous metabolic traits (TAG, TC, LDL-C, HDL-C, apoA-1, apoB and HbA1c) in the following way: model 1 adjusted for sex, age (in years, continuous). Model 2 additionally adjusted for potential confounders including physical activity (self-reported) that was derived into a four-scale index by combining levels of occupational and recreational physical activity( Reference Wareham, Jakes and Rennie 32 ) (inactive, moderately inactive, moderately active and active), smoking status (current smoker, former smoker and never smoker), total energy intake (EI, in kJ, continuous), lipid-lowering medication, social class (professional, managerial and technical, skilled non-manual, partly skilled and unskilled), and educational level (low education, O-level or equivalent to secondary school, A-level or equivalent including technical school, university degree or equivalent, including a higher vocational qualification). Model 3 additionally adjusted for BMI (continuous) and WC (continuous). We also examined the effects of per two-point increase in the rMED (as a continuous variable). The association between adherence to the MD and hypertriacylglycerolaemia status (defined according to the cut-off points ( ≥ 1·7 mmol TAG) was determined using multivariable-adjusted logistic regression models. We performed sensitivity analyses by excluding participants with chronic prevalent disease (heart disease, stroke, diabetes mellitus and/or cancer) and excluding misreporters of energy. Misreporting of EI was estimated by using the ratio of reported EI:predicted BMR (EI:BMR). Participants were classified as under-reporters (EI:BMR < 1·14), plausible reporters (EI:BMR = 1·14–2·1) or over-reporters (EI:BMR>2·1) of EI by using cut-off points proposed by Goldberg( Reference Goldberg, Black and Jebb 33 ).

Among the 19 800 individuals with available data, we analysed the rs780094 SNP using a co-dominant mode of inheritance as in previous studies( Reference Qi, Wu and Li 19 , Reference Sotos-Prieto, Guillén and Vicente Sorli 23 , Reference Bi, Kao and Boerwinkle 34 ). In multiple linear regression analysis, we adjusted for relevant covariates, with model 1 including age (continuous) and sex, and model 2 additionally including BMI and WC (both continuous) and use of lipid-lowering medication (no or yes). Multivariable logistic regression was used to estimate the OR of hypertriacylglycerolaemia associated with the polymorphism.

Joint effects of genotype (GG, GA and AA) and rMED (low, medium and high adherence to the MD) on continuous metabolic traits were examined using ANCOVA with nine possible combinations (high adherence+GG, high adherence+GA, high adherence+AA, medium adherence+GG, medium adherence+GA, medium adherence+AA, low adherence+GG, low adherence+GA and low adherence+AA).

To examine the gene × MD adherence interaction on the levels of cardiometabolic markers (TAG, TC, apoB and HbA1c), we used multivariable linear regression models including main effects and interaction terms.

Statistical analyses were conducted using SPSS version 15.0 for Windows (SPSS, Inc.). Statistical significance was set at the 0·05 level, and all tests were two-tailed. In addition, we also applied the Bonferroni method for multiple comparisons, and set a lower statistical significance accounting for the number of tests of significance( Reference Bender and Lange 35 ).

Results

Among the 20 986 individuals, 33 % (n 6924), 51 % (n 10 627), and 16 % (n 3435) were in the low, medium and high MD adherence categories (rMED).

Overall, a higher adherence to the MD was observed among women, individuals with a lower BMI and WC, older, physically active people, and individuals with a higher level of education (Table 1).

Table 1 Baseline characteristics according to the levels of adherence to Mediterranean dietary patterns (relative Mediterranean Diet Score (rMED) categories) in the European Prospective Investigation into Cancer (EPIC)-Norfolk study participants (Mean values and standard deviations; percentages; median values and interquartile ranges (IQR))

WC, waist circumference, PA, physical activity.

* Low adherence to the Mediterranean diet (MD; rMED 0–6); medium adherence to the MD (rMED 7–10); high adherence to the MD (rMED 11–18).

Total: n 20 986; low: n 6924 (33 %); medium: n 10 627 (51 %); high: n 3435 (16 %).

Association between Mediterranean diet adherence and metabolic markers

Table 2 shows that TAG concentrations and the apoB:apoA-1 ratio in the medium and high categories of the rMED were significantly lower, while HDL-C and apoA-1 concentrations were higher compared with the concentrations in the lowest rMED category. Further mutual adjustment for other lipids subtypes did not alter the findings. HbA1c concentration was lower with an increasing rMED. After Bonferroni correction for multiple comparisons, each of the concentrations of TAG, HDL-C and HbA1c remained significantly associated with the rMED.

Table 2 Association between adherence to the Mediterranean diet (MD) according to the relative Mediterranean Diet Score (rMED) categories and metabolic markers in the European Prospective Investigation into Cancer-Norfolk study participants (Mean or geometric mean values and 95 % confidence intervals)

TC, total cholesterol; HDL-C, HDL cholesterol; LDL-C, LDL cholesterol; HbA1c, glycated Hb.

* Low adherence to the MD (rMED 0–6); medium adherence to the MD (rMED 7–10); high adherence to the MD (rMED 11–18).

Low adherence: n 6924 (33 %); medium: n 10 627 (51 %); high: n 3435 (16 %).

P value is shown for the one-factor ANOVA test. A lower P value of 0·00 625 was applied after using the Bonferroni correction method for multiple comparisons.

§ All the variables are log-transformed.

Model 1: adjusted for age and sex.

Model 2: model 1+social class, educational level, physical activity, smoking status, energy intake and lipid-lowering medication.

** Model 3: model 2+BMI and waist circumference.

†† HbA1c n 10 746, with HbA1c levels < 6·5 %.

Each two-point increment in the rMED was significantly inversely associated with TAG (β − 0·005, 95 % CI − 0·007, − 0·002), LDL-C (β − 0·001, 95 % CI − 0·003, 0·000), apoB:apoA-1 ratio (β − 0·002, 95 % CI − 0·005, 0·000) and HbA1c (β − 0·001, 95 % CI − 0·002, 0·000), while it was directly associated with HDL-C (β 0·004, 95 % CI 0·002, 0·005) and apoA-1 (β 0·006, 95 % CI 0·002, 0·009) (β coefficients for log-transformed variables). The rMED was inversely associated with hypertriacylglycerolaemia (TAG>1·7 mmol/l), with OR 0·89, 95 % CI 0·81, 0·98; model 3. Sensitivity analyses did not alter the main findings.

Association between glucokinase regulatory protein SNP and metabolic markers

Baseline characteristics according to the GCKR rs780094 SNP showed that no significant differences were observed in relation to age, sex, BMI, WC, social class and physical activity. Significant differences were found for educational level and lipid medication (see online Supplementary Table S1). The AA genotype (frequency 15·6 %) was associated with 0·16 mmol/l higher TAG concentrations compared with GG individuals in adjusted analyses (Table 3). The AA genotype was also significantly positively associated with TC and apoB, but no significant differences were found for other lipid parameters (HDL-C and apoA-1) by genotype (Table 3). We found a tendency towards lower HbA1c concentration for the A allele, but this was not statistically significant. Application of the Bonferroni method for multiple comparisons did not alter these findings. Both GA and AA genotypes v. GG genotype were associated with a greater risk of hypertriacylglycerolaemia in adjusted analyses (GA: OR 1·23, 95 % CI 1·15, 1·33; AA: OR 1·47, 95 % CI 1·33, 1·62).

Table 3 Association between genotypes in the glucokinase regulatory protein (GCKR) gene (rs780094 (G>A)) and metabolic markers in the European Prospective Investigation into Cancer (EPIC)-Norfolk study participants (Mean or geometric mean values and 95 % confidence intervals)

TC, total cholesterol; HDL-C, HDL cholesterol; LDL-C, LDL cholesterol; HbA1c, glycated Hb.

* GG: n 7372 (37·2 %); GA: n 9337 (47·2 %), AA: n 3091 (15·6 %).

P value is shown for the one-factor ANOVA test. A lower P value of 0·00 625 was applied after using the Bonferroni correction method for multiple comparisons.

All the variables are log-transformed.

§ Model 1: adjusted for age and sex.

Model 2: model 1+BMI+waist circumference+lipid-lowering medication.

HbA1c n 10 746, with HbA1c levels < 6·5 %.

Joint results of genotype and Mediterranean diet adherence on metabolic markers

Table 4 demonstrates the joint results of genotype and MD adherence on metabolic parameters. Participants with the highest risk (AA genotype and with the lowest rMED) had 12·5 % higher TAG concentrations (mean 1·68 (95 % CI 1·63, 1·73) mmol/l) than those who had the lowest risk (mean 1·47 (95 % CI 1·45, 1·50) mmol/l; GG genotype with the highest rMED; Table 4). Similar results were found for TC concentrations (2·3 % higher) and apoB (3·1 % higher), respectively (Table 4). For HbA1c, the converse was the case with those who were GG homozygous with the lowest rMED having 3·3 % higher HbA1c concentration (Table 4). However, the association of MD adherence and metabolic markers did not differ significantly by genotype, with no significant gene–diet interactions for lipids or for HbA1c (TAG, P= 0·770; TC, P= 0·761; apoB, P= 0·860; HbA1c, P= 0·481). These findings were unchanged with correction for multiple comparisons using the Bonferroni method.

Table 4 Combined association between the relative Mediterranean Diet Score (rMED; low, medium or high adherence to Mediterranean dietary patterns) and the genotype (rs780094 (G>A) SNP in the glucokinase regulatory protein (GCKR) gene) on lipid concentrations and glycated Hb (HbA1c) (Mean or geometric mean values and 95 % confidence intervals)

TC, total cholesterol.

* GG: n 7372 (37·2 %); GA: n 9337 (47·2 %), AA: n 3091 (15·6 %).

P value for difference for the one-factor ANOVA test (model adjusted for age sex, BMI and lipid-lowering medication).

All the variables are log-transformed.

§ P value for comparison between each genotype according to the rMED. A lower P value of 0·0125 was applied after using the Bonferroni correction method for multiple comparisons.

Similar results were found when we analysed the rs1260326 SNP (in high linkage disequilibrium with the rs780094 SNP, r 2 0·94) in a small sample (n 7273; data not shown). Furthermore, we performed a subgroup analysis stratified by age ( < 59·1 years old v. >59·1 years old; cut-off points selected by median), physical activity (inactive or moderately inactive v. active or moderately active) and smoking status (smoker v. non-smoker) for the outcomes analysed in order to examine whether in any of these strata, a significant gene–diet interaction was obtained. We did not find any interaction, though there were differences by genotype in TAG levels in both strata for each of those three factors (data not shown).

Discussion

The findings from the present large population-based study in the UK suggest that higher adherence to the MD is associated with important metabolic benefits for lipid and glycaemic parameters, while the AA genotype of the GCKR rs780094 SNP is associated with an adverse distribution of lipid parameters. We also report that the TAG-raising A allele of rs780094 is associated with a trend towards a HbA1c-lowering effect. While there were significant differences in metabolic parameter distribution when comparing those with the highest and lowest risks (i.e. those with highest genetic susceptibility and lowest MD adherence v. those with lowest genetic susceptibility and highest MD adherence), the benefits of greater MD adherence were not differential by genotype.

The present study extends previous findings that were demonstrated for TAG concentrations and genotype in a high-cardiovascular risk Spanish population of 945 older adults( Reference Sotos-Prieto, Guillén and Vicente Sorli 23 ), to several metabolic parameters in a general British population in the EPIC-Norfolk study. The results of the present study indicate that a high adherence to the MD as defined by the rMED is associated with an improved metabolic profile (lower TAG concentrations, higher HDL-C concentrations and lower HbA1c levels). We provide evidence that the rs789004 SNP in the GCKR gene is associated with TAG concentrations, confirming previous findings( Reference Saxena, Voight and Lyssenko 17 , Reference Qi, Wu and Li 19 , Reference Vaxillaire, Cavalcanti-Proença and Dechaume 22 , Reference Bi, Kao and Boerwinkle 34 , Reference Chasman, Paré and Zee 36 ), and additionally report that carriers of the minor A allele also have higher TC and apoB concentrations. We further demonstrate that there are no gene–diet interactions on metabolic parameters, but that greater adherence to the MD confers a more favourable metabolic profile irrespective of genetic susceptibility.

The results of the present study are in general agreement with other cross-sectional, prospective and interventional studies of the association between MD adherence and markers of cardiometabolic risk( Reference Estruch, Martínez-González and Corella 8 , Reference Elhayany, Lustman and Abel 9 , Reference Tzima, Pitsavos and Panagiotakos 11 , Reference Chrysohoou, Panagiotakos and Pitsavos 37 Reference Athyros, Kakafika and Papageorgiou 42 ). However, although some studies have found inverse associations between MD adherence and TC( Reference Ambring, Friberg and Axelsen 5 , Reference Estruch, Martínez-González and Corella 8 ), we did not find such an association for TC. A weak association was shown for LDL-C when we considered rMED as continuous (per two points of increase in the rMED). These findings could be attributable to the differences in the definition of rMED in different studies. We also found an inverse association between the rMED and HbA1c levels, a finding that is in line with our recent report of an inverse association between adherence to the MD and incident diabetes( Reference Romaguera, Guevara and Norat 2 ).

Thus far, the association between MD adherence and apolipoproteins has not been widely studied. Notably, we found a direct association of MD adherence with apoA-1 concentration, but this was of borderline significance after Bonferroni correction for multiple comparisons, though it could be of potential clinical significance. This finding in a British general population is consistent with a randomised study in a high-cardiovascular risk Spanish population( Reference Solá, Fitó and Estruch 14 ). Our findings raise the possibility that MD adherence may be useful for achieving healthy apolipoprotein goals in the general British population, but require confirmation in future studies.

We confirmed the genetic association of the rare variant of the rs780094 SNP with higher TAG concentrations as described previously( Reference Willer, Sanna and Jackson 16 Reference Qi, Wu and Li 19 , Reference Vaxillaire, Cavalcanti-Proença and Dechaume 22 ). In the first study that identified this SNP in the GCKR gene as a target for TAG concentrations, this SNP was reported to explain 1 % of the residual variance in TAG levels( Reference Saxena, Voight and Lyssenko 17 ). In the present study, 0·5 % of the residual variance in TAG concentration is explained by this SNP, which is similar to that reported in Scandinavian samples (0·4–1·2 %)( Reference Orho-Melander, Melander and Guiducci 18 ) and higher than the variance reported in a French population (0·2 %)( Reference Vaxillaire, Cavalcanti-Proença and Dechaume 22 ). Past studies have found opposing associations of this SNP with TAG (elevated concentrations) and glucose (lowered concentrations)( Reference Saxena, Voight and Lyssenko 17 , Reference Orho-Melander, Melander and Guiducci 18 , Reference Dupuis, Langenberg and Prokopenko 21 , Reference Vaxillaire, Cavalcanti-Proença and Dechaume 22 , Reference Bi, Kao and Boerwinkle 34 ). We have extended this finding by examining the association with HbA1c concentration, not previously reported in the literature, and we report that the TAG-raising A allele of rs780094 is associated with an HbA1c-lowering effect. A mutational mechanism for the reported association of this variant with raised TAG and lower glucose levels has been proposed, such that the GCKR T-allele Pro446Leu (C>T; rs1260326) has reduced regulation by physiological concentrations of phosphate esters fructose 6, resulting indirectly in increased glucokinase activity( Reference Beer, Tribble and McCulloch 43 ). Altered glucokinase regulation in the liver is predicted to enhance glycolytic flux, promoting hepatic glucose metabolism and elevating concentrations of malonyl-CoA, a substrate for de novo lipogenesis, and thus relating to the reported association of this variant with raised TAG and lower glucose levels. Taken together, these data support an important role of GCKR in pathways regulating hepatic TAG as well as glucose metabolism in humans.

With respect to the study of genotype and other metabolic traits, we found an association between the rs780094 SNP and TC and apoB. Few studies have analysed other metabolic lipid or apolipoprotein traits in relation to this SNP( Reference Orho-Melander, Melander and Guiducci 18 Reference Sparsø, Andersen and Nielsen 20 , Reference Chasman, Paré and Zee 36 ). Sparsø et al. ( Reference Sparsø, Andersen and Nielsen 20 ) found a modest association between rs780094 and TC, and interestingly, in agreement with Chasman et al. ( Reference Chasman, Paré and Zee 36 ), apoB concentrations are also higher in carriers of the risk allele in the present study.

Although no significant gene–diet interaction was found, this is the first study to examine how a SNP (rs780094) in the GCKR gene previously shown to be consistently associated with TAG concentrations is associated not only with TAG levels, but also with other cardiometabolic risk phenotypes (TC, apoB and HbA1c) according to the groups of adherence to the MD and genetic risk allele carriers.

The strengths of the present study include a large population-based sample size of approximately 20 000 individuals, and a well-phenotyped dataset that enabled us to account for a number of relevant plausible confounders and to estimate differences and interactions with a high degree of precision. A further strength is that we have extended previous findings from nutrigenetic studies that have focused on associations between single diet or nutrient components and endpoints to examine dietary patterns by assessing the rMED score of MD adherence. The analysis of dietary patterns may provide a more accurate depiction of peoples' eating habits and yield better models for generating public health recommendations regarding healthy eating( Reference Hu 44 ). Limitations of the present study include the cross-sectional design of the study that does not allow temporal relationships to be examined. The FFQ for dietary assessment is prone to measurement error, and its use at a single time point in the present study does not permit the examination of any changes in diet over time, notwithstanding, the EPIC-Norfolk FFQ is a validated instrument( Reference Kaaks and Riboli 29 ). We had non-fasted samples in the present study, due to logistical issues relating to recruiting participants to a large epidemiological study. However, for glycaemia, we used HbA1c levels, where fasting is not required, and we do not believe that the absence of fasted lipid samples affects our overall findings and interpretation, in keeping with recent evidence( Reference Di Angelantonio, Sarwar and Perry 12 , Reference Sidhu and Naugler 45 ).

In conclusion, the present study provides evidence on the potential protective properties of the MD, evaluated as an a priori-defined score, on the markers of cardiometabolic risk, and confirms the association of genetic variation in GCKR with TAG concentrations, adding new information on its association with other lipids and apolipoproteins. The finding of both separate and joint associations of diet and genotype, but not an interaction between them, with the metabolic profile should be further investigated in other populations and in prospective studies.

Supplementary material

To view supplementary material for this article, please visit http://dx.doi.org/10.1017/S0007114514000580

Acknowledgements

We are grateful to the participants of the EPIC-Norfolk Study, and acknowledge the contribution of the EPIC-Norfolk study team as well as the general practitioners. We thank Amit Bhaniani for dataset preparation and Stephen Sharp for helpful discussions.

The EPIC-Norfolk study was supported by programme grants from the Medical Research Council UK and Cancer Research UK. We acknowledge core MRC Epidemiology Unit support through Programmes MC_UU_12015/1 and MC_UU_12015/5. M. S.-P. was further supported by AP2007-04663, FPU, Ministerio de Educación y Ciencia. The funders had no role in the design, analysis or writing of this article.

The contributions of the authors are as follows: M. S-.P. and N. G. F. conceptualised and designed the analysis; K-.T. K. and N. J. W. were the principal investigators of the EPIC-Norfolk study; R. L. provided information management for the study; M. S.-P. performed all the statistical analyses; M. S.-P. and N. G. F. wrote the manuscript; M. S.-P. and N. G. F. had responsibility for the final content; N. G. F. was the guarantor. All co-authors provided important intellectual contribution and edited the manuscript. All authors read and approved the final version.

None of the authors has any conflict of interest to declare.

References

1 Trichopoulou, A, Costacou, T, Bamia, C, et al. (2003) Adherence to a Mediterranean diet and survival in a Greek population. N Engl J Med 348, 25992608.CrossRefGoogle Scholar
2 Romaguera, D, Guevara, M, Norat, T, et al. (2011) Mediterranean diet and type 2 diabetes risk in the European Prospective Investigation into Cancer and Nutrition (EPIC) study: the InterAct project. Diabetes Care 34, 19131918.Google Scholar
3 Truthmann, J, Richter, A, Thiele, S, et al. (2012) Associations of dietary indices with biomarkers of dietary exposure and cardiovascular status among adolescents in Germany. Nutr Metab 9, 92.Google Scholar
4 Panagiotakos, DB, Tzima, N, Pitsavos, C, et al. (2007) The association between adherence to the Mediterranean diet and fasting indices of glucose homoeostasis: the ATTICA Study. J Am Coll Nutr 26, 3238.Google Scholar
5 Ambring, A, Friberg, P, Axelsen, M, et al. (2004) Effects of a Mediterranean-inspired diet on blood lipids, vascular function and oxidative stress in healthy subjects. Clin Sci (Lond) 106, 519525.Google Scholar
6 Andreoli, A, Lauro, S, Di Daniele, N, et al. (2008) Effect of a moderately hypoenergetic Mediterranean diet and exercise program on body cell mass and cardiovascular risk factors in obese women. Eur J Clin Nutr 62, 892897.Google Scholar
7 Goulet, J, Lamarche, B, Nadeau, G, et al. (2003) Effect of a nutritional intervention promoting the Mediterranean food pattern on plasma lipids, lipoproteins and body weight in healthy French-Canadian women. Atherosclerosis 170, 115124.Google Scholar
8 Estruch, R, Martínez-González, MA, Corella, D, et al. (2006) Effects of a Mediterranean-style diet on cardiovascular risk factors: a randomized trial. Ann Intern Med 145, 111.Google Scholar
9 Elhayany, A, Lustman, A, Abel, R, et al. (2010) A low carbohydrate Mediterranean diet improves cardiovascular risk factors and diabetes control among overweight patients with type 2 diabetes mellitus: a 1-year prospective randomized intervention study. Diabetes Obes Metab 12, 204209.CrossRefGoogle ScholarPubMed
10 Kastorini, CM, Milionis, HJ, Esposito, K, et al. (2011) 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 57, 12991313.Google Scholar
11 Tzima, N, Pitsavos, C, Panagiotakos, DB, et al. (2007) Mediterranean diet and insulin sensitivity, lipid profile and blood pressure levels, in overweight and obese people; the Attica study. Lipids Health Dis 6, 22.Google Scholar
12 Di Angelantonio, E, Sarwar, N, Perry, P, et al. (2009) Major lipids, apolipoproteins, and risk of vascular disease. JAMA 302, 19932000.Google Scholar
13 Emerging Risk Factor Collaboration Di Angelantonio, E, Gao, P, Pennells, L, et al. (2012) Lipid-related markers and cardiovascular disease prediction. JAMA 307, 24992506.Google Scholar
14 Solá, R, Fitó, M, Estruch, R, et al. (2011) Effect of a traditional Mediterranean diet on apolipoproteins B, A-I, and their ratio: a randomized, controlled trial. Atherosclerosis 218, 174180.Google Scholar
15 Corella, D & Ordovas, JM (2005) Single nucleotide polymorphisms that influence lipid metabolism: interaction with dietary factors. Annu Rev Nutr 25, 341390.Google Scholar
16 Willer, CJ, Sanna, S, Jackson, AU, et al. (2008) Newly identified loci that influence lipid concentrations and risk of coronary artery disease. Nat Genet 40, 161169.Google Scholar
17 Saxena, R, Voight, BF, Lyssenko, V, et al. (2007) Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels. Science 316, 13311336.Google Scholar
18 Orho-Melander, M, Melander, O, Guiducci, C, et al. (2008) Common missense variant in the glucokinase regulatory protein gene is associated with increased plasma triglyceride and C-reactive protein but lower fasting glucose concentrations. Diabetes 57, 31123121.Google Scholar
19 Qi, Q, Wu, Y, Li, H, et al. (2009) Association of GCKR rs780094, alone or in combination with GCK rs1799884, with type 2 diabetes and related traits in a Han Chinese population. Diabetologia 52, 834843.CrossRefGoogle ScholarPubMed
20 Sparsø, T, Andersen, G, Nielsen, T, et al. (2008) The GCKR rs780094 polymorphism is associated with elevated fasting serum triacylglycerol, reduced fasting and OGTT-related insulinaemia, and reduced risk of type 2 diabetes. Diabetologia 51, 7075.Google Scholar
21 Dupuis, J, Langenberg, C, Prokopenko, I, et al. (2010) New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk. Nat Genet 42, 105116.CrossRefGoogle ScholarPubMed
22 Vaxillaire, M, Cavalcanti-Proença, C, Dechaume, A, et al. (2008) The common P446L polymorphism in GCKR inversely modulates fasting glucose and triglyceride levels and reduces type 2 diabetes risk in the DESIR prospective general French population. Diabetes 57, 22532257.Google Scholar
23 Sotos-Prieto, M, Guillén, M, Vicente Sorli, J, et al. (2013) Relevant associations of the glucokinase regulatory protein/glucokinase gene variation with TAG concentrations in a high-cardiovascular risk population. Modulation by Mediterranean diet. Br J Nutr 109, 193201.CrossRefGoogle Scholar
24 Day, N, Oakes, S, Luben, R, et al. (1999) EPIC-Norfolk: study design and characteristics of the cohort. Br J Cancer 80, 95103.Google Scholar
25 Riboli, E (1992) Nutrition and cancer: background and rationale of the European Prospective Investigation into Cancer and Nutrition (EPIC). Ann Oncol 3, 783791.Google Scholar
26 Friedewald, WT, Levy, RI & Fredrickson, DS (1972) Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem 18, 499502.Google Scholar
27 Bingham, SA, Gill, C, Welch, A, et al. (1997) 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 26, Suppl. 1, S137S151.Google Scholar
28 Welch, AA, Luben, R, Khaw, KT, et al. (2005) The CAFE computer program for nutritional analysis of the EPIC-Norfolk food frequency questionnaire and identification of extreme nutrient values. J Hum Nutr Diet 18, 99116.Google Scholar
29 Kaaks, R & Riboli, E (1997) Validation and calibration of dietary intake measurements in the EPIC project: methodological considerations. Int J Epidemiol 26, Suppl. 1, S15S25.Google Scholar
30 Buckland, G, González, CA, Agudo, A, et al. (2009) Adherence to the Mediterranean diet and risk of coronary heart disease in the Spanish EPIC Cohort Study. Am J Epidemiol 170, 15181529.Google Scholar
31 Buckland, G, Agudo, A, Lujan, L, et al. (2010) Adherence to a Mediterranean diet and risk of gastric adenocarcinoma within the European Prospective Investigation into Cancer and Nutrition (EPIC) cohort study. Am J Clin Nutr 91, 381390.Google Scholar
32 Wareham, NJ, Jakes, RW, Rennie, KL, et al. (2003) Validity and repeatability of a simple index derived from the short physical activity questionnaire used in the European Prospective Investigation into Cancer and Nutrition (EPIC) study. Public Health Nutr 6, 407413.Google Scholar
33 Goldberg, GR, Black, AE, Jebb, SA, et al. (1991) Critical evaluation of energy intake data using fundamental principles of energy physiology: 1. Derivation of cut-off limits to identify under-recording. Eur J Clin Nutr 45, 569581.Google Scholar
34 Bi, M, Kao, WH, Boerwinkle, E, et al. (2010) Association of rs780094 in GCKR with metabolic traits and incident diabetes and cardiovascular disease: the ARIC Study. PLOS ONE 5, e11690.Google Scholar
35 Bender, R & Lange, S (2001) Adjusting for multiple testing – when and how? J Clin Epidemiol 54, 343349.Google Scholar
36 Chasman, DI, Paré, G, Zee, RY, et al. (2008) Genetic loci associated with plasma concentration of low-density lipoprotein cholesterol, high-density lipoprotein cholesterol, triglycerides, apolipoprotein A1, and apolipoprotein B among 6382 white women in genome-wide analysis with replication. Circ Cardiovasc Genet 1, 2130.Google Scholar
37 Chrysohoou, C, Panagiotakos, DB, Pitsavos, C, et al. (2004) Adherence to the Mediterranean diet attenuates inflammation and coagulation process in healthy adults: the ATTICA Study. J Am Coll Cardiol 44, 152158.Google Scholar
38 Shai, I, Schwarzfuchs, D, Henkin, Y, et al. (2008) Weight loss with a low-carbohydrate, Mediterranean, or low-fat diet. N Engl J Med 359, 229241.CrossRefGoogle ScholarPubMed
39 Esposito, K, Marfella, R, Ciotola, M, et al. (2004) Effect of a Mediterranean-style diet on endothelial dysfunction and markers of vascular inflammation in the metabolic syndrome: a randomized trial. JAMA 292, 14401446.Google Scholar
40 Esposito, K, Maiorino, MI, Di Palo, C, et al. (2009) Adherence to a Mediterranean diet and glycaemic control in type 2 diabetes mellitus. Diabet Med 26, 900907.Google Scholar
41 Papadaki, A & Scott, JA (2008) Follow-up of a web-based tailored intervention promoting the Mediterranean diet in Scotland. Patient Educ Couns 73, 256263.CrossRefGoogle ScholarPubMed
42 Athyros, VG, Kakafika, AI, Papageorgiou, AA, et al. (2011) Effect of a plant stanol ester-containing spread, placebo spread, or Mediterranean diet on estimated cardiovascular risk and lipid, inflammatory and haemostatic factors. Nutr Metab Cardiovasc Dis 21, 213221.Google Scholar
43 Beer, NL, Tribble, ND, McCulloch, LJ, et al. (2009) The P446L variant in GCKR associated with fasting plasma glucose and triglyceride levels exerts its effect through increased glucokinase activity in liver. Hum Mol Genet 18, 40814088.Google Scholar
44 Hu, FB (2002) Dietary pattern analysis: a new direction in nutritional epidemiology. Curr Opin Lipidol 13, 39.Google Scholar
45 Sidhu, D & Naugler, C (2012) Fasting time and lipid levels in a community-based population: a cross-sectional study. Arch Intern Med 172, 17071710.CrossRefGoogle Scholar
Figure 0

Table 1 Baseline characteristics according to the levels of adherence to Mediterranean dietary patterns (relative Mediterranean Diet Score (rMED) categories) in the European Prospective Investigation into Cancer (EPIC)-Norfolk study participants (Mean values and standard deviations; percentages; median values and interquartile ranges (IQR))

Figure 1

Table 2 Association between adherence to the Mediterranean diet (MD) according to the relative Mediterranean Diet Score (rMED) categories and metabolic markers in the European Prospective Investigation into Cancer-Norfolk study participants (Mean or geometric mean values and 95 % confidence intervals)

Figure 2

Table 3 Association between genotypes in the glucokinase regulatory protein (GCKR) gene (rs780094 (G>A)) and metabolic markers in the European Prospective Investigation into Cancer (EPIC)-Norfolk study participants (Mean or geometric mean values and 95 % confidence intervals)

Figure 3

Table 4 Combined association between the relative Mediterranean Diet Score (rMED; low, medium or high adherence to Mediterranean dietary patterns) and the genotype (rs780094 (G>A) SNP in the glucokinase regulatory protein (GCKR) gene) on lipid concentrations and glycated Hb (HbA1c) (Mean or geometric mean values and 95 % confidence intervals)

Supplementary material: PDF

Sotos-Prieto Supplementary Material

Table S1

Download Sotos-Prieto Supplementary Material(PDF)
PDF 32.2 KB