Intended for healthcare professionals

CCBY Open access
Research

Association of plasma biomarkers of fruit and vegetable intake with incident type 2 diabetes: EPIC-InterAct case-cohort study in eight European countries

BMJ 2020; 370 doi: https://doi.org/10.1136/bmj.m2194 (Published 08 July 2020) Cite this as: BMJ 2020;370:m2194
  1. Ju-Sheng Zheng, principal investigator1 2,
  2. Stephen J Sharp, senior statistician1,
  3. Fumiaki Imamura, senior investigator scientist1,
  4. Rajiv Chowdhury, principal researcher3,
  5. Thomas E Gundersen, chief scientific officer4,
  6. Marinka Steur, career development fellow1,
  7. Ivonne Sluijs, assistant professor5,
  8. Yvonne T van der Schouw, professor5,
  9. Antonio Agudo, unit head6,
  10. Dagfinn Aune, associate professor7 8 9,
  11. Aurelio Barricarte, epidemiologist10,
  12. Heiner Boeing, professor11,
  13. María-Dolores Chirlaque, specialist in preventive medicine and public health12 13,
  14. Miren Dorronsoro, unit chief14,
  15. Heinz Freisling, scientist15,
  16. Douae El-Fatouhi, biostatistician16,
  17. Paul W Franks, professor17,
  18. Guy Fagherazzi, research leader16 18,
  19. Sara Grioni, research assistant19,
  20. Marc J Gunter, section head15,
  21. Cecilie Kyrø, senior researcher20,
  22. Verena Katzke, postdoctoral scientist21,
  23. Tilman Kühn, group head21,
  24. Kay-Tee Khaw, professor22,
  25. Nasser Laouali, postdoctoral researcher16,
  26. Giovanna Masala, medical epidemiologist23,
  27. Peter M Nilsson, professor17,
  28. Kim Overvad, professor24 25,
  29. Salvatore Panico, professor26,
  30. Keren Papier, nutritional epidemiologist27,
  31. J Ramón Quirós, unit head28,
  32. Olov Rolandsson, professor29,
  33. Daniel Redondo-Sánchez, research assistant13 30 31,
  34. Fulvio Ricceri, epidemiologist32 33,
  35. Matthias B Schulze, professor34 35 36,
  36. Annemieke M W Spijkerman, senior scientist37,
  37. Anne Tjønneland, research head20,
  38. Tammy Y N Tong, nutritional epidemiologist27,
  39. Rosario Tumino, senior pathologist38,
  40. Elisabete Weiderpass, director15,
  41. John Danesh, professor3 39 40 41,
  42. Adam S Butterworth, reader3 39,
  43. Elio Riboli, professor7,
  44. Nita G Forouhi, professor1,
  45. Nicholas J Wareham, professor1
  1. 1MRC Epidemiology Unit, University of Cambridge School of Clinical Medicine, Institute of Metabolic Science, Cambridge Biomedical Campus, Cambridge CB2 0QQ, UK
  2. 2Key Laboratory of Growth Regulation and Translation Research of Zhejiang Province, School of Life Sciences, Westlake University, Hangzhou, China
  3. 3MRC/BHF Cardiovascular Epidemiology Unit, Department of Public Health and Primary Care, University of Cambridge, Cambridge, UK
  4. 4Vitas AS, Oslo, Norway
  5. 5Julius Centre for Health Sciences and Primary Care, University Medical Center Utrecht, Utrecht University, Utrecht, Netherlands
  6. 6Unit of Nutrition and Cancer, Cancer Epidemiology Research Programme, Catalan Institute of Oncology, Group of Research on Nutrition and Cancer, Bellvitge Biomedical Research Institute, L'Hospitalet of Llobregat, Barcelona, Spain
  7. 7School of Public Health, Imperial College London, London, UK
  8. 8Department of Nutrition, Bjørknes University College, Oslo, Norway
  9. 9Department of Endocrinology, Morbid Obesity and Preventive Medicine, Oslo University Hospital, Oslo, Norway
  10. 10Navarra Public Health Institute, Pamplona, Spain
  11. 11Department of Epidemiology, German Institute of Human Nutrition Potsdam-Rehbruecke, Nuthetal, Germany
  12. 12Department of Epidemiology, Regional Health Council, IMIB-Arrixaca, Murcia University, Murcia, Spain
  13. 13CIBER in Epidemiology and Public Health (CIBERESP), Madrid, Spain
  14. 14Public Health Division of Gipuzkoa, San Sebastian, Spain
  15. 15International Agency for Research on Cancer, Lyon, France
  16. 16Centre of Research in Epidemiology and Population Health, UMR 1018 Inserm, Institut Gustave Roussy, Paris-Sud Paris-Saclay University, Villejuif, France
  17. 17Department of Clinical Sciences, Lund University, Malmö, Sweden
  18. 18Luxembourg Institute of Health (LIH), Department of Population Health, Strassen, Luxembourg
  19. 19Epidemiology and Prevention Unit, Milan, Italy
  20. 20Danish Cancer Society Research Centre, Copenhagen, Denmark
  21. 21Division of Cancer Epidemiology, German Cancer Research Centre, Heidelberg, Germany
  22. 22Department of Public Health and Primary Care, University of Cambridge, Cambridge, UK
  23. 23Cancer Risk Factors and Life-Style Epidemiology Unit, Institute for Cancer Research, Prevention and Clinical Network, Florence, Italy
  24. 24Department of Public Health, Aarhus University, Aarhus, Denmark
  25. 25Department of Cardiology, Aalborg University Hospital, Aarhus, Denmark
  26. 26Department of Clinical Medicine and Surgery, Federico II University, Naples, Italy
  27. 27Cancer Epidemiology Unit, Nuffield Department of Population Health, University of Oxford, Oxford, UK
  28. 28Public Health Directorate, Asturias, Spain
  29. 29Department of Public Health and Clinical Medicine, Family Medicine, Umeå University, Umeå, Sweden
  30. 30Andalusian School of Public Health, Granada, Spain
  31. 31Institute of Biosanitary Research of Granada, Granada, Spain
  32. 32Department of Clinical and Biological Sciences, University of Turin, Turin, Italy
  33. 33Unit of Epidemiology, Regional Health Service ASL TO3, Grugliasco, Turin, Italy
  34. 34Department of Molecular Epidemiology, German Institute of Human Nutrition Potsdam-Rehbruecke, Nuthetal, Germany
  35. 35German Centre for Diabetes Research, Neuherberg, Germany
  36. 36University of Potsdam, Institute of Nutritional Sciences, Potsdam, Germany
  37. 37National Institute for Public Health and the Environment, Bilthoven, Netherlands
  38. 38Cancer Registry and Histopathology Department, Azienda Sanitaria Provinciale, Ragusa, Italy
  39. 39NIHR Blood and Transplant Research Unit in Donor Health and Genomics, Department of Public Health and Primary Care, University of Cambridge, Cambridge, UK
  40. 40British Heart Foundation Cambridge Centre of Excellence, Division of Cardiovascular Medicine, Addenbrooke's Hospital, Cambridge, UK
  41. 41Department of Human Genetics, Wellcome Trust Sanger Institute, Hinxton, Cambridge, UK
  1. Correspondence to: N G Forouhi nita.forouhi{at}mrc-epid.cam.ac.uk
  • Accepted 2 April 2020

Abstract

Objective To investigate the association of plasma vitamin C and carotenoids, as indicators of fruit and vegetable intake, with the risk of type 2 diabetes.

Design Prospective case-cohort study.

Setting Populations from eight European countries.

Participants 9754 participants with incident type 2 diabetes, and a subcohort of 13 662 individuals from the European Prospective Investigation into Cancer and Nutrition (EPIC) cohort of 340 234 participants: EPIC-InterAct case-cohort study.

Main outcome measure Incident type 2 diabetes.

Results In a multivariable adjusted model, higher plasma vitamin C was associated with a lower risk of developing type 2 diabetes (hazard ratio per standard deviation 0.82, 95% confidence interval 0.76 to 0.89). A similar inverse association was shown for total carotenoids (hazard ratio per standard deviation 0.75, 0.68 to 0.82). A composite biomarker score (split into five equal groups), comprising vitamin C and individual carotenoids, was inversely associated with type 2 diabetes with hazard ratios 0.77, 0.66, 0.59, and 0.50 for groups 2-5 compared with group 1 (the lowest group). Self-reported median fruit and vegetable intake was 274 g/day, 396 g/day, and 508 g/day for participants in categories defined by groups 1, 3, and 5 of the composite biomarker score, respectively. One standard deviation difference in the composite biomarker score, equivalent to a 66 (95% confidence interval 61 to 71) g/day difference in total fruit and vegetable intake, was associated with a hazard ratio of 0.75 (0.67 to 0.83). This would be equivalent to an absolute risk reduction of 0.95 per 1000 person years of follow up if achieved across an entire population with the characteristics of the eight European countries included in this analysis.

Conclusions These findings indicate an inverse association between plasma vitamin C, carotenoids, and their composite biomarker score, and incident type 2 diabetes in different European countries. These biomarkers are objective indicators of fruit and vegetable consumption, and suggest that diets rich in even modestly higher fruit and vegetable consumption could help to prevent development of type 2 diabetes.

Introduction

The global burden of type 2 diabetes has risen over the past decades and its prevention is a public health priority.1 High fruit and vegetable intake has been suggested to have an important role in prevention of this disorder.2 Evidence from prospective studies linking fruit and vegetable intake with type 2 diabetes is inconsistent and weak,345, and evidence from randomised controlled trials is sparse.6 Previous research studies have typically used dietary food frequency questionnaires to assess fruit and vegetable intake, which are subject to measurement error and recall bias.7 Circulating plasma vitamin C and carotenoids have been proposed as objective biomarkers of fruit and vegetable intake, with evidence for their validity from observational and experimental studies.89101112 In meta-analyses of randomised controlled trials, groups provided with more fruits and vegetables had increased blood concentrations of a panel of fruit and vegetable related biomarkers.101112 A meta-analysis of up to 96 intervention studies found that blood vitamin C and several carotenoids were the most consistently responsive biomarkers for fruit and vegetable intake.11 In individual participant meta-analysis of controlled intervention studies, evidence was found of a positive dose-response association between fruit and vegetable consumption and biomarker concentrations.12 Moreover, in comparative analyses, compliance at group level with fruit and vegetable interventions was indicated equally well by blood vitamin C or individual carotenoids.10 Thus, investigation of the association between these biomarkers and type 2 diabetes could provide insight into the association of fruit and vegetable intake with this disorder.

Results from previous small studies describing the association between circulating plasma carotenoids and the incidence of type 2 diabetes are inconclusive, with an inverse association being reported in some studies9131415 and others finding no significant association.16 Only one previously published prospective cohort study has examined circulating vitamin C and incident type 2 diabetes17 reporting an inverse association between the concentration of plasma vitamin C and incident type 2 diabetes in a population in the United Kingdom, but evidence is lacking in other countries or populations with different lifestyles and dietary behaviours.

Our study aimed at examining the association of baseline levels of circulating vitamin C and carotenoids with incident type 2 diabetes in the European Prospective Investigation into Cancer and Nutrition (EPIC)-InterAct study, which is based on more than 340 000 community based adults from eight European countries. We also aimed to construct a composite biomarker score to examine the association of the combination of biomarkers with incident type 2 diabetes.

Methods

Study design and population

The EPIC-InterAct study is a prospective case-cohort study, nested within the European EPIC study.18 In brief, cases of incident type 2 diabetes occurring in eight of the 10 EPIC countries (Denmark, France, Germany, Italy, Netherlands, Spain, Sweden, and UK) between 1991 and 2007 were ascertained and verified by the InterAct consortium across 26 study centres. We ascertained and verified 12 403 individuals with incident type 2 diabetes over 3.99 million person years of follow-up from a cohort of 340 234 participants with stored blood and buffy coat. From this cohort, a centre stratified subcohort was assembled by randomly selecting 16 835 individuals. A total of 16 154 participants remained in the subcohort after exclusion of 548 with prevalent diabetes, 129 with uncertain diabetes status, and 4 with diabetes after censoring. The case-cohort design has the advantages of temporal sequence and power of a cohort study (in that it involves the complete number of incident cases) with the measurement efficiency of a case-control study. The random subcohort included 778 participants with verified incident type 2 diabetes, which is a design feature of the case-cohort approach.18 We further excluded 6543 participants with no plasma samples available for the measurement of plasma vitamin C (including all Swedish samples (n=5401)) and 4952 with no carotenoids (including one Swedish centre, Malmo (n=3556)). We subsequently excluded those with samples of low volume or samples that failed in biochemical analysis or were haemolysed (n=218 for vitamin C; n=85 for carotenoids). We, therefore, included a total of 22 833 participants, with 9754 participants with incident type 2 diabetes and 13 662 subcohort participants, with an average follow up of 9.7 years (supplemental figure S1). As a design feature of the case-cohort study the final eligible subcohort included 583 participants with incident type 2 diabetes. For the analysis of different biomarkers, the sample size varied according to the number of missing data for the biomarkers or covariates, and the final sample size for each biomarker analysis is presented.

All participants gave written informed consent and the study was approved by the local ethics committee in the participating centres, and the internal review board of the International Agency for Research on Cancer.

Ascertainment and verification of cases of type 2 diabetes

Ascertainment of cases of incident type 2 diabetes up until 31 December 2007 involved a review of self-reported data, linkage to primary and secondary care registers, use of medication (drug registers), hospital admissions, and mortality data, which we have described in detail previously.18 No cases were ascertained solely by self-report because we confirmed any self-reported case with at least one other independent source of information. Cases of type 2 diabetes in Denmark and Sweden were not ascertained by self-report but were identified through local and national diabetes and pharmaceutical registers, and were considered verified.

Measurement of plasma vitamin C and carotenoids

Non-fasting blood samples were collected at the baseline visit of the EPIC study. Plasma samples stored at −196°C (or at −150°C for samples from Denmark) were shipped to Vitas Analytical Services (Oslo, Norway) for the measurement of plasma vitamin C and six individual carotenoids (α carotene, β carotene, lycopene, lutein, zeaxanthin, β cryptoxanthin) by high performance liquid chromatography-ultraviolet methods, as described previously, by staff masked to participant information.1920 Assays for plasma vitamin C and carotenoids were performed on separate plasma aliquots, with treatment methods and chromatographic procedures.1920 For the vitamin C assay, meta-phosphoric acid was added to stabilise the vitamin C and prevent oxidation. For carotenoids, internal standard β apo-8-carotenal and butylated hydroxytoluene (as an antioxidant) were added before treatment. Coefficients of variation were 4.2-4.5% for vitamin C and 2.7-6.7% for carotenoids. The lower limit of detection was 0.7 µmol/L for vitamin C, 0.009 µmol/L for α carotene and β carotene, 0.018 µmol/L for lycopene, 0.009 µmol/L for lutein and zeaxanthin, and 0.007 µmol/L for β cryptoxanthin. We imputed random values between zero and the lower limit of detection for those below the lower limit (1.0% for vitamin C, 1.7% for α carotene, 6.7% for zeaxanthin, and <0.5% for others).

To assess the reproducibility of the assay of vitamin C and carotenoids after long term storage, we compared these biomarkers in the period of measurement for this study (2013-14) with those previously measured in 2004917 in a subset of Norfolk (UK) participants of the EPIC study (n=1582 for vitamin C; n between 520 and 678 for individual carotenoids). The reproducibility for the measurement of these biomarkers was high for most biomarkers (r=0.89-0.93), and moderate for zeaxanthin (r=0.63).

Measurement of other baseline characteristics

Baseline weight, height, and waist circumference were collected by trained health professionals during a visit to a study centre, except in France and in the Oxford cohort.21 Baseline dietary information was collected using a dietary questionnaire that was self-administered or administered by an interviewer, which was developed within each country to estimate the usual food intake of participants.212223 Information on baseline physical activity, smoking status, sociodemographic factors, and medical history was obtained from self-administered questionnaires.21

Serum lipid biomarkers (except for Umeå centre in Sweden, where plasma was used), including total cholesterol, high density lipoprotein cholesterol, and triglycerides were measured at SHL (Stichting Ingenhousz Laboratory, Etten-Leur, Netherlands). Low density lipoprotein cholesterol was calculated based on the Friedewald formula.24 Both high and low density lipoprotein cholesterol were considered as confounders because both can influence carotenoid levels (as fat soluble vitamins) and could influence the association between nutritional biomarker and disease.925

Statistical analysis

All the statistical analyses were performed in Stata 14 (Statacorp, College Station, TX). The concentration of total carotenoids was calculated as the sum of the six individual carotenoids. Plasma vitamin C, total, and six individual carotenoids (α carotene, β carotene, lycopene, lutein, zeaxanthin, β cryptoxanthin) were the exposures for this analysis. A composite biomarker score was generated by calculating the average of the standardised values (each value was standardised using the mean and standard deviation from the subcohort) of vitamin C and six individual carotenoids, as previously described.9 The Spearman correlations between plasma vitamin C and six individual carotenoids were calculated in the subcohort. In cross sectional analyses, to identify potential correlates of plasma vitamin C and carotenoids, we log transformed the values of the total and individual carotenoids, owing to their skewed distributions. In the subcohort, we used linear regression to estimate country-specific β coefficients and 95% confidence intervals, representing the associations of demographic, lifestyle, and dietary factors with plasma vitamin C, total and six individual carotenoids, and the composite biomarker score, with mutual adjustment for all the demographic, lifestyle, and dietary factors in the model. We then pooled the estimated associations across countries using random effects meta-analysis.

To help interpret the results of the composite biomarker score in relation to fruit and vegetable intake, we calculated the difference in total fruit and vegetable intake for each one standard deviation higher composite biomarker score in the subcohort using a linear regression. These analyses included adjustment for age, sex, centre, physical activity, smoking status, employment, marital status, education, alcohol intake, total energy intake, body mass index, waist circumference, and dietary covariates (intake of potatoes, cereal and cereal products, milk and dairy products, fruit and vegetable juice, soft drinks, fish, red meat, legumes, egg and egg products, nuts and seeds, offal, and vitamin supplement use). In addition, we used a regression model with restricted cubic splines (three knots) to assess the shape of the association of the composite biomarker score with fruit and vegetable intake in the subcohort, adjusted for the same covariates as above.

We used Prentice weighted Cox regression26 to estimate country-specific hazard ratio and 95% confidence interval for incident type 2 diabetes, comparing biomarkers divided into five equal groups (from the lowest, group 1 to the highest, group 5) and for each one standard deviation (groups 1-5 and the standard deviation were both calculated in the subcohort). We then pooled the estimates using random effects meta-analysis. We fitted three statistical models, adjusting for potential sociodemographic, lifestyle behavioural, biochemical, and anthropometric confounders: model 1a adjusted for age (as underlying time scale, continuous), sex, and centre; model 1b, as model 1a, plus physical activity (inactive, moderately inactive, moderately active, active), smoking status (never, former, current), employment (no, yes), marital status (single, married, separated/divorced, widowed), education (low, middle, high), alcohol drinking (never, 0 to <6, 6 to <12, 12 to <24, and ≥24 g/day), total energy intake (continuous), high density lipoprotein cholesterol (continuous; only for carotenoids analyses) and low density lipoprotein cholesterol (continuous; only for carotenoids analyses); and model 2, as model 1b, plus body mass index (continuous) and waist circumference (continuous). Based on model 2, we further modelled restricted cubic spline terms (three knots) for analysis of each biomarker to assess the shape of their associations with incident type 2 diabetes.2728 To investigate the association with type 2 diabetes of an estimated daily fruit and vegetable intake that meets or exceeds the “five-a-day” recommendation (≥400 g/day) compared with not meeting the recommendation (<400 g/day), we categorised participants into two groups: biomarker predicted five or more portions a day and fewer than five portions a day, using the value of the composite biomarker score, which equated to 400 grams a day of fruit and vegetable intake in the prediction model, as the cut-off point.

In sensitivity analyses, we imputed missing values of covariates using multiple imputation with chained equations, 29 to assess the effect of missing data on the biomarker-disease associations; we performed analysis in each of 10 imputed datasets and combined the estimates using Rubin’s rules. The imputation model included all covariates from the analysis model, the event variable, and the Nelson-Aalen estimate of cumulative hazard.29 To examine the influence of each individual biomarker on the composite biomarker score results, we amended the composite biomarker score by excluding one biomarker at a time, and then examined its association with type 2 diabetes. We also assessed robustness of the results from model 2 to take accountof other potential confounders, including mutual adjustment for the other individual nutritional biomarkers (vitamin C and six individual carotenoids), use of vitamin supplements, season of blood draw, family history of diabetes, baseline prevalence of stroke, coronary heart diseases or cancers (excluding non-melanoma skin cancer), hormone use and menopausal status in women, genetic risk score for body mass index and for insulin resistance, and a diet quality score (Mediterranean diet score). The genetic risk scores for body mass index and insulin resistance were both unweighted, and were calculated as the sum of the number of risk alleles for body mass index (97 genetic variants)30 or insulin resistance (53 genetic variants).31 We also accounted for potential reverse causality by excluding people with haemoglobin A1c greater than or equal to 6.5% (48 mmol/mol) at baseline or those confirmed as having type 2 diabetes within the first two years or first four years after baseline. We performed further sensitivity analysis by excluding participants with baseline cancer or cardiovascular disease or by restricting the analysis to the first eight years of follow-up. We tested for multiplicative interaction between each nutritional biomarker and the following prespecified variables: age, sex, body mass index, physical activity, smoking status, season of blood draw, and vitamin supplement use. We used a lower threshold for significance based on the Bonferroni correction. Analyses within subgroups defined by the variables were performed if the P value for interaction was <0.05.

Patient and public involvement

No patients were involved in setting the research question, nor were they involved in the design and implementation of the study. No plans exist to involve patients in dissemination of the results.

Results

Population characteristics

Baseline mean (standard deviation) concentration of plasma vitamin C and total carotenoids was 36.3 (18.3) and 1.3 (0.7) µmol/L, respectively, among individuals who developed type 2 diabetes compared with 42.3 (19.2) µmol/L and 1.7 (0.8) µmol/L in subcohort participants (table 1). In the whole subcohort, participants in Germany had the highest mean levels of plasma vitamin C (48.4 (standard deviation 19.2) µmol/L), whereas those from Italy had the lowest levels (36.6 (18.6) µmol/L). Participants from France had the highest mean plasma total carotenoid levels (2.3 (1.0) µmol/L), whereas those from Denmark had the lowest levels (1.3 (0.8) µmol/L; table 1). Plasma vitamin C, total and individual carotenoids were all positively correlated with each other (supplemental table S1).

Table 1

Distribution of plasma vitamin C and carotenoids by case status and by population characteristics in the subcohort of EPIC-InterAct study

View this table:

Association of plasma vitamin C, total and individual carotenoids with potential correlates

After adjustment for demographic, lifestyle, and dietary factors, positive associations were found for sex (women compared with men), physical activity, and education level with both plasma vitamin C and total carotenoids, while waist circumference and current smoking status (compared with never smoking) were inversely associated (supplemental figure S2, supplemental table S2).

Dietary intake of total fruits and vegetables, citrus fruits, non-citrus fruits, fruiting vegetables, and fruit and vegetable juice were all positively associated with both plasma vitamin C, total carotenoids, and the composite biomarker score (fig 1). Each 100 g/day higher intake of fruits and vegetables was associated with 0.10 higher level of the composite biomarker score (fig 1); conversely, every one standard deviation higher composite biomarker score was associated with a 66 (95% confidence interval 61 to 71) g/day higher intake of fruits and vegetables. The median fruit and vegetable intake was 274, 357, 396, 452, and 508 g/day for levels 1 (lowest) to 5 (highest) of the composite biomarker score, respectively (table 2). Forest plots showing the association of fruit and vegetable intake subgroups with each plasma biomarker by country and overall are presented in supplemental figure S3; the directions of these associations were generally consistent across countries.

Fig 1
Fig 1

Association of dietary factors with plasma vitamin C and carotenoids in the subcohort of the European Prospective Investigation into Cancer and Nutrition (EPIC)-InterAct study. The values in the box represent the differences in vitamin C, carotenoids, or the composite biomarker score (in standard deviaiton units) for each one standardised unit per category difference in dietary factors. The mean (standard deviation) or percentage of participants for the examined dietary factors are presented in the second column. Linear regression was used to obtain the country specific estimate of an association, adjusting for demographic, lifestyle, and other examined dietary factors if available in that country, except the overlapping food groups. The country-specific estimates were then combined using random effects meta-analysis. *P<0.05; †P<0.001; ‡P<0.0001. Sample size was 10 584 for vitamin C and 11 537 for the carotenoid estimation. All the values are expressed in a red scale for different levels of positive associations and blue scale for different levels of negative associations. All the carotenoids variables were natural log transformed before statistical analysis in the linear regression. Standard deviation was 19.2 µmol/L for vitamin C, 0.50 for total carotenoids (log transformed), 0.91 for α carotene (log transformed), 0.73 for β carotene (log transformed), 0.68 for lycopene (log transformed), 0.51 for lutein (log transformed), 0.91 for zeaxanthin (log transformed), 0.94 for β cryptoxanthin (log transformed), and 0.57 for the composite biomarker score. A β coefficient of 0.10 represents the increase in the composite biomarker score per 100 g/day of fruit and vegetable intake. Performing this cross sectional analysis the other way round, every change in one standard deviation in the composite biomarker score was associated with a 66 g/day increase in fruit and vegetable intake

Table 2

Prospective associations between plasma vitamin C and carotenoids and incident type 2 diabetes in the EPIC-InterAct study*

View this table:

Association of plasma vitamin C, total and individual carotenoids with incident type 2 diabetes

Higher levels of plasma vitamin C were associated with a lower hazard of type 2 diabetes; hazard ratio for each standard deviation 0.71 (95% confidence interval 0.65 to 0.78) in a model adjusted for age, sex, and study centre (model 1a). Hazard ratios comparing groups 2, 3, 4, and 5 of vitamin C with group 1 (the lowest group) were 0.84 (0.77 to 0.92), 0.61 (0.53 to 0.71), 0.53 (0.46 to 0.61), and 0.39 (0.31 to 0.50). These inverse associations were attenuated after further adjustment for socioeconomic factors, lifestyle factors, and adiposity (model 2), but the trend across the groups remained: hazard ratios 0.94 (0.83 to 1.08), 0.72 (0.60 to 0.87), 0.68 (0.56 to 0.82), and 0.58 (0.47 to 0.72) comparing groups 2, 3, 4, and 5, with group 1, respectively (table 2).

For plasma total carotenoids, the inverse association with type 2 diabetes estimated in models 1a and 1b was attenuated after further adjustment for adiposity (model 2), but remained significant. The hazard ratio for each standard deviation of total carotenoids was 0.75 (95% confidence interval 0.68 to 0.82). The hazard ratios comparing groups 2, 3, 4, and 5 of total carotenoids with group 1 were 0.83 (0.74 to 0.92), 0.64 (0.54 to 0.75), 0.61 (0.53 to 0.70), and 0.51 (0.43 to 0.60), respectively. In model 2, all plasma individual carotenoids were inversely associated with type 2 diabetes, except for zeaxanthin (table 2). For plasma total and individual carotenoids, the heterogeneity (I2) across countries ranged from 48.0% for β cryptoxanthin to 75.2% for β carotene (supplemental figure S4), but the directions of the associations were consistent across different countries.

The composite biomarker score, which included contributions from all the examined biomarkers, was inversely associated with type 2 diabetes. The hazard ratios comparing groups 2, 3, 4, and 5 of the combined biomarker score with group 1 were 0.77 (0.68 to 0.87), 0.66 (0.54 to 0.80), 0.59 (0.48 to 0.72), and 0.50 (0.40 to 0.62), respectively. One standard deviation difference in the composite biomarker score, equivalent to a 66 (95% confidence interval 61 to 71) g/day difference in total fruit and vegetable intake, was associated with a hazard ratio of 0.75 (95% confidence interval 0.67 to 0.83). This would be equivalent to an absolute risk reduction of 0.95 per 1000 person years of follow-up if achieved across an entire population with the characteristics of the eight European countries included in this analysis. Thus the crude overall incidence rate of type 2 diabetes would be estimated to be reduced from 3.8 to 2.85 per 1000 person years of follow-up in this population.

We found evidence of non-linear associations (P value for non-linearity <0.001) for plasma total carotenoids, α carotene, β carotene, and lutein, with a strong inverse association with type 2 diabetes at low to moderate concentrations, but weaker at higher concentrations (fig 2). The shape of the associations of the composite biomarker score with total fruit and vegetable intake and with type 2 diabetes is presented in figure 3. In an analysis of the association with type 2 diabetes using this biomarker score as a cut-off level to define consumption of five or more portions a day of fruits and vegetables, the hazard ratio compared with consuming an estimated fewer than five portions a day was 0.69 (95% confidence interval 0.63 to 0.76).

Fig 2
Fig 2

Associations of plasma vitamin C and carotenoids with incident type 2 diabetes in the European Prospective Investigation into Cancer and Nutrition (EPIC)-InterAct Study. The purple solid line and the shaded area represent estimates of hazard ratios and the 95% confidence intervals, respectively, for each biomarker (median in the subcohort as a reference). Covariates included age (as underlying time scale), sex, centre, physical activity, smoking status, employment, marital status, education, alcohol intake, total energy intake, high density lipoprotein cholesterol (for analyses of carotenoids only), low density lipoprotein cholesterol (for analyses of carotenoids only), body mass index, and waist circumference. For total carotenoids, α carotene, β carotene, and lutein, we found evidence of a non-linear association (P value for non-linearity <0.001)

Fig 3
Fig 3

Associations of plasma composite biomarker score with incident type 2 diabetes and total fruit and vegetable intake in the European Prospective Investigation into Cancer and Nutrition (EPIC)-InterAct Study. The purple solid line and shaded areas represent hazard ratios and 95% confidence intervals for the association of the composite biomarker score with type 2 diabetes (median as a reference), adjusted for age (as underlying time scale), sex, centre, physical activity, smoking status, employment, marital status, education, alcohol intake, total energy intake, high density lipoprotein cholesterol, low density lipoprotein cholesterol, body mass index, and waist circumference. The red dashed lines represent the association of the composite biomarker score with fruit and vegetable intake (best estimates and their 95% confidence intervals), adjusted for age, sex, centre, physical activity, smoking status, employment, marital status, education, alcohol intake, total energy intake, body mass index, waist circumference, and dietary intake of potatoes, cereal and cereal products, milk and dairy products, fruit and vegetable juice, soft drinks, fish, red meat, legumes, egg and egg products, nuts and seeds, offals, and vitamin supplement use. P values for non-linearity were 0.056 for the hazard ratios, and <0.001 for estimates for fruit and vegetable intake. The vertical lines within the figure represent the quintiles (20th, 40th, 60th, and 80th centiles) of the composite biomarker score in the subcohort

In a sensitivity analysis, after mutual adjustment for the other individual biomarkers, the inverse association with type 2 diabetes for each of the biomarkers was attenuated (supplemental table S3). Lycopene and β cryptoxanthin were no longer significantly associated with type 2 diabetes. Other sensitivity analyses had little effect on the findings (supplemental tables S3 and S4).

We identified evidence of a significant interaction only between baseline age and both plasma vitamin C and α carotene (P value for interaction <0.001), with the inverse associations with type 2 diabetes being weaker in the oldest age group (supplemental figure S5).

Discussion

In this large study with dietary diversity across eight European countries, higher concentrations of plasma vitamin C and carotenoids were associated with a lower incidence of type 2 diabetes. A composite biomarker score, comprising the seven examined biomarkers, was associated with consumption of fruit and vegetables and was also inversely related to the risk of type 2 diabetes. These findings provide strong evidence from objectively measured biomarkers for the recommendation that fruit and vegetable intake should be increased to prevent type 2 diabetes.

Strengths and limitations of this study

The strengths of our study include our report of the association of plasma vitamin C and individual carotenoids with type 2 diabetes in a large case-cohort study. The study included a subcohort of 13 662 individuals and 9754 participants with incident type 2 diabetes, from an original cohort of 340 234 participants with 3.99 million person years of follow-up. This analysis, based on predominantly white participants from eight European countries, indicates the generalisability of the study findings to populations of European descent, though further research in other populations is warranted. The number of participants with incident type 2 diabetes in this study was 10-fold greater than in previous studies of the association of biomarkers of fruit and vegetable intake with risk of type 2 diabetes, which not only provided high statistical power for the main effect analyses but also permitted a dose-response analysis and other sensitivity analyses.

This study has several limitations. The biomarkers could degrade during storage over years, which would make it difficult to quantify the absolute concentration of vitamin C and carotenoids. The likelihood of degradation in this study was minimised, however, by storing the samples in liquid nitrogen and if it did occur, it should be the same in cases and non-cases. The vitamin C and carotenoid biomarkers were measured only at baseline, and intra-individual variation over time is likely. However, when follow-up was restricted to the first 8 years, the associations remained similar. Another limitation is the observational nature of this study and the potential for residual confounding due to mis-measured or unmeasured confounders. The inverse associations we observed could, in part, be explained by the confounding effect of a healthy lifestyle and favourable socioeconomic conditions, which might not have been fully accounted for in the analyses. Moreover, differential misclassification in ascertainment of type 2 diabetes could have occurred if those with high fruit and vegetable consumption were more health conscious, sought medical assessment, and hence were diagnosed with the disorder. This possibility could cause bias towards the null and an underestimation of the inverse association. The results could be explained by reverse causality, but our sensitivity analysis indicated that this was unlikely. An assessment of causal associations of these nutritional biomarkers with type 2 diabetes was beyond the scope of our study. In addition, a common limitation for nutritional biomarkers is that their relationship with diet is often influenced by both genetic variation and nutrient metabolism.32 Further challenges in interpretation arise. Thus the presence of these biomarkers might not be exclusively due to the consumption of fruit and vegetables, but some carotenoids, for instance, could also originate from the food production process, such as food colouring. The biomarkers could reflect underlying dietary patterns or food substitution effects when fruit or vegetables replace meat in the diet, for instance.

Comparison with other studies

Evidence from previous studies for the magnitude of the association between circulating vitamin C and carotenoids and the risk of type 2 diabetes is limited, mainly because of the small number of participants in previous studies.13141617 Additionally, those studies did not investigate differences in those associations by country. Biomarkers can be used not only to investigate the relationship between total fruit and vegetables and the risk of diabetes but can also help to examine the effect of specific fruit and vegetable subgroups. For example, α carotene and β carotene were highly correlated with root vegetables, in this and a previous study.8 Therefore, the biomarker finding of α carotene and β carotene could suggest a potential inverse association with type 2 diabetes for root vegetables; this deserves further study. In this study, the inverse association between some biomarkers and type 2 diabetes was attenuated when mutually adjusted for the other individual biomarkers. These results suggest that a combination of nutrients in fruit and vegetables explains the association with type 2 diabetes rather than a single nutrient. Thus our finding of an inverse association between the biomarker composite score and risk of type 2 diabetes could be capturing this combined effect.

Our results could also provide a possible explanation of the previous null findings of randomised trials for single vitamin supplements. In the Women’s Antioxidant Cardiovascular Study, dietary supplements of vitamin C, vitamin E, or β carotene had no significant effects on the risk of developing type 2 diabetes during a follow-up of 9.2 years.33 Similarly, long term β carotene supplementation did not affect the risk of incident type 2 diabetes among healthy men in a randomised trial over a 12 year follow-up period,34 or in male smokers in a randomised trial with 6.1 years of follow-up.35 The absence of a significant intervention effect on the risk of developing type 2 diabetes in these trials could have been because they focused on supplementation of these individual nutrients themselves rather than the foods or a dietary pattern for which these nutrients are objective biomarkers. Fruits and vegetables contain many other components besides vitamin C and carotenoids, including fibre, phytochemicals such as flavonoids, and other antioxidants, all of which could have beneficial health effects.36 Alternatively, the null results for vitamin supplementation trials could be explained by the focus of those trials on individual rather than composite nutrients, or their focus on increasing average levels rather than reducing deficiency. We observed non-linear associations for carotenoids, and thus the effect of supplementation on type 2 diabetes could be greater in those who had lower levels at baseline. Finally, the distribution of vitamin levels differed in those trials from that in our observational study. We found a mean concentration of plasma β carotene of 0.47 µmol/L. This result contrasts with a mean of 2.24 µmol/L in the Women’s Antioxidant Cardiovascular Study33 and 5.6 µmol/L in the Alpha-Tocopherol, Beta-Carotene Cancer Prevention Study,35 after β carotene intervention of 50 mg every other day or 20 mg each day, respectively.

The potential mechanisms of association could be that consumption of fruit and vegetables helps to regulate weight and adiposity, as well as glucose-insulin homoeostasis and the inflammatory status of the participants,3738 thus leading to a reduced risk of type 2 diabetes. Another possible mechanism is the influence of the fruit and vegetable intake on the gut microbiota, given their high concentration of fibre or other gut microbiota-related nutrients.3940 A sensitivity analysis of additional adjustment for a diet quality score did not alter the findings for any of the biomarkers, which further confirms that higher fruit and vegetable intake is likely to be beneficial for prevention of type 2 diabetes regardless of the overall diet quality. We could not rule out the possibility, however, that the observed association partly reflected metabolic effects independent of dietary factors.

What this study adds

Various dietary guidelines have recommended increasing fruit and vegetable intake as an important component of a healthy diet. However, evidence derived from a food frequency questionnaire for the specific role of fruit and vegetables and their subtypes in the prevention of type 2 diabetes has previously been weak and inconsistent.3441 The potential overall benefits of fruits and vegetables have also been questioned within certain popular dietary regimens that favour low carbohydrate intake, including advice to limit the consumption of many fruits and vegetables.4243 Although five portions a day of fruit and vegetables have been recommended for decades, in 2014-15, 69% of UK adults ate fewer than this number,44 and this proportion is even higher in European (EU) adults (86%).45 The low population level concordance with the “five a day” recommendation provides an incentive to quantify the benefits of making small changes in consumption of fruits and vegetables even below the threshold of the widely recommended guideline level, as suggested by our findings from the dose-response relationship of the biomarkers with fruit and vegetable intake presented in figure 3. Using biomarkers to quantify fruit and vegetable intake is an adjunct to self-reported questionnaires as objective assessment methods are less affected by the measurement errors of subjective dietary assessment tools. The use of biomarkers in this context is at present a research tool used to advance scientific understanding rather than a suggestion for their use in clinical practice.

This study also extends our understanding of the correlates of plasma vitamin C, total and individual carotenoids in different European countries. Previous studies,84647 based on a subsample of about 3000 people in the EPIC cohort, reported that the strongest predictors of individual carotenoids were fruits for β cryptoxanthin, total carrots and root vegetable for α carotene, and tomato products for lycopene. Our findings with a larger subcohort of people from the EPIC study are broadly consistent with the previous work, and provide additional information about the correlates of vitamin C, another strong biomarker of fruit and vegetable intake. In addition to the association of plasma vitamin C and carotenoids with fruit and vegetable intake, we have also described their association with other food groups, and demographic and lifestyle factors. We found a positive association between consumption of fruit and vegetable juice and the composite biomarker score, which could reflect a combination of factors, including naturally occurring vitamin C and the content of carotenoids in juice as well as fortified products. We could not distinguish between these, but our analyses of association between the biomarkers and incident type 2 diabetes were adjusted for vitamin supplement use, and suggest that as biomarkers of fruit and vegetable intake these findings endorse the consumption of fruit and vegetables, not that of supplements.

Conclusions

Higher levels of plasma vitamin C, total and individual carotenoid biomarkers, and their composite biomarker score were associated with a lower incidence of type 2 diabetes in diverse European populations in eight countries. These biomarkers are indicators of fruit and vegetable consumption, reducing the measurement error and bias of dietary self-reports. Our findings suggest that higher fruit and vegetable consumption is inversely associated with the incidence of type 2 diabetes, regardless of whether this increase in consumption is from a level below or above the recommended five-a-day threshold. The public health implication of this observation is that the consumption of even a moderately increased amount of fruit and vegetables among populations who typically consume low levels could help to prevent type 2 diabetes. It should be noted that these findings and other available evidence suggest that fruit and vegetable intake, rather than vitamin supplements, is potentially beneficial for the prevention of type 2 diabetes.

What is already known on this topic

  • Investigation of a link between fruit and vegetable intake and the risk of type 2 diabetes has relied on self-reported dietary questionnaires, with inconsistent findings

  • Evidence from objective markers of fruit and vegetable intake is sparse but would be complementary to self-report

What this study adds

  • Results from the assessment of individual blood biomarkers and a composite biomarker score comprising plasma vitamin C and carotenoids indicate an inverse association with new onset type 2 diabetes in a pan-European population

  • This study suggests that even a modest increase in fruit and vegetable intake could help to prevent type 2 diabetes, indicated by objective biomarkers of consumption, regardless of whether the increase is among people with initially low or high intake

Acknowledgments

We thank all EPIC participants and staff for their contribution to the study. We thank Nicola Kerrison (MRC Epidemiology Unit, Cambridge) for managing the data for the InterAct Project. We thank the technical and functional operational teams of the MRC Epidemiology Unit and laboratory team at VITAS AS, Norway, for the measurements of plasma vitamin C and carotenoids. We thank staff from the EPIC-CVD and EPIC-InterAct Coordinating Centres for carrying out sample preparation and data-handling work, particularly Sarah Spackman (EPIC-CVD Data Manager), and Cambridge Genomic Services for genotyping.

Footnotes

  • Contributors: NGF and NJW contributed equally to this work. NJW was chief investigator and NGF and SJS were investigators of the InterAct project. J-SZ and NGF had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis; J-SZ performed the statistical analyses and drafted the manuscript with input from NGF. The working group (J-SZ, SJS, FI, TEG, MS, RC, YTvdS, IS, AB, NGF) and all authors contributed to interpretation of data, revised the article critically for important intellectual content, and approved the final version of the manuscript. NGF and NJW are the guarantors, and they contributed equally to the work. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted.

  • Funding: The InterAct project was funded by the EU FP6 programme (grant number LSHM_CT_2006_037197). Biomarker measurements for vitamin C and carotenoids were funded jointly by the InterAct project, the EPIC-CVD project, and the MRC Cambridge Initiative (RG71466, SJAH/004). EPIC-CVD has been supported by the UK Medical Research Council (MR/L003120/1), the British Heart Foundation (RG/13/13/30194; RG/18/13/33946), the European Commission Framework Programme 7 (HEALTH-F2-2012-279233), the European Research Council (268834), and the National Institute for Health Research (NIHR; Cambridge Biomedical Research Centre at the Cambridge University Hospitals NHS Foundation Trust). JD is funded by the National Institute for Health Research (senior investigator award). NJW and NGF acknowledge funding from the following agencies: Medical Research Council Epidemiology Unit MC_UU_12015/1 and MC_UU_12015/5, and NIHR Biomedical Research Centre Cambridge: nutrition, diet, and lifestyle research theme (IS-BRC-1215-20014). JSZ has received funding from Westlake University (No YSYY0209) and European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 701708. TYNT has received funding from the Cancer Research UK (C570/A16491 and C8221/A19170), UK Medical Research Council (MR/M012190/1). PWF has received funding from Swedish Research Council, Novo Nordisk, Swedish Diabetes Association, Swedish Heart-Lung Foundation, European Research Council. K-TK has received funding from Medical Research Council UK, Cancer Research UK. TK has received funding from German Cancer Aid, German Cancer Research Center (DKFZ), German Federal Ministry of Education and Research (BMBF). PMN has received funding from Swedish Research Council. KO has received funding from Danish Cancer Society. SP has received funding from Compagnia di San Paolo. JRQ has received funding from regional government of Asturias. OR has received funding from the Västerboten County Council. ER has received funding from Imperial College Biomedical Research Centre. EPIC Bilthoven and Utrecht acknowledge the Dutch Ministry of Public Health, Welfare and Sports, Netherlands Cancer Registry, LK Research Funds, Dutch Prevention Funds, Dutch ZON (Zorg Onderzoek Nederland), World Cancer Research Fund (WCRF). The coordination of EPIC-Spain was financially supported by the Health Research Fund (FIS) - Instituto de Salud Carlos III (ISCIII), regional governments of Andalucía, Asturias, Basque Country, Murcia, and Navarra, and the Catalan Institute of Oncology (ICO) (Spain). The funders of the studies had no role in the study design, data collection, data analysis, data interpretation, or report preparation.

  • Competing interests: All authors have completed the ICMJE uniform disclosure form at www.icmje.org/coi_disclosure.pdf and declare: support from European commission, UK Medical Research Council, British Heart Foundation, and National Institute for Health Research Cambridge Biomedical Research Centre for the submitted work; no financial relationships with any organisations that might have an interest in the submitted work in the previous three years; no other relationships or activities that could appear to have influenced the submitted work. Where authors are identified as personnel of the International Agency for Research on Cancer/World Health Organization, the authors alone are responsible for the views expressed in this article and they do not necessarily represent the decisions, policy, or views of the International Agency for Research on Cancer/WHO.

  • Ethical approval: All included studies were approved by local review boards, and all participants gave written informed consent to participate in the study.

  • Data sharing: EPIC-InterAct Study data cannot be deposited publicly as these collaborative data originate from multiple research institutions across eight European countries with different legal frameworks. The authors confirm that researchers seeking the analysis dataset for this work can submit a data request to the EPIC-InterAct study central contact point by emailing interact@mrc-epid.cam.ac.uk.

  • The corresponding author affirms that the manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as planned (and, if relevant, registered) have been explained.

  • Dissemination to participants and related patient and public communities: Findings will be disseminated via the communications team of the corresponding author’s institute.

http://creativecommons.org/licenses/by/4.0/

This is an Open Access article distributed in accordance with the terms of the Creative Commons Attribution (CC BY 4.0) license, which permits others to distribute, remix, adapt and build upon this work, for commercial use, provided the original work is properly cited. See: http://creativecommons.org/licenses/by/4.0/.

References