Hostname: page-component-7c8c6479df-ws8qp Total loading time: 0 Render date: 2024-03-28T18:54:04.267Z Has data issue: false hasContentIssue false

Does a better adherence to dietary guidelines reduce mortality risk and environmental impact in the Dutch sub-cohort of the European Prospective Investigation into Cancer and Nutrition?

Published online by Cambridge University Press:  03 August 2017

Sander Biesbroek*
Affiliation:
Centre for Nutrition, Prevention and Health Services, National Institutefor Public Health and the Environment (RIVM), Antonie van Leeuwenhoeklaan 9, Bilthoven 3721 MA, The Netherlands
W. M. Monique Verschuren
Affiliation:
Centre for Nutrition, Prevention and Health Services, National Institutefor Public Health and the Environment (RIVM), Antonie van Leeuwenhoeklaan 9, Bilthoven 3721 MA, The Netherlands Julius Centre for Health Sciences and Primary Care, University Medical Centre Utrecht, Universiteitsweg 100, 3584 CG Utrecht, The Netherlands
Jolanda M. A. Boer
Affiliation:
Centre for Nutrition, Prevention and Health Services, National Institutefor Public Health and the Environment (RIVM), Antonie van Leeuwenhoeklaan 9, Bilthoven 3721 MA, The Netherlands
Mirjam E. van de Kamp
Affiliation:
Centre for Nutrition, Prevention and Health Services, National Institutefor Public Health and the Environment (RIVM), Antonie van Leeuwenhoeklaan 9, Bilthoven 3721 MA, The Netherlands
Yvonne T. van der Schouw
Affiliation:
Julius Centre for Health Sciences and Primary Care, University Medical Centre Utrecht, Universiteitsweg 100, 3584 CG Utrecht, The Netherlands
Anouk Geelen
Affiliation:
Division of Human Nutrition, Wageningen University & Research, PO Box 17, 6700 AA Wageningen, The Netherlands
Moniek Looman
Affiliation:
Division of Human Nutrition, Wageningen University & Research, PO Box 17, 6700 AA Wageningen, The Netherlands
Elisabeth H. M. Temme
Affiliation:
Centre for Nutrition, Prevention and Health Services, National Institutefor Public Health and the Environment (RIVM), Antonie van Leeuwenhoeklaan 9, Bilthoven 3721 MA, The Netherlands
*
*Corresponding author: S. Biesbroek, email sander.biesbroek@rivm.nl
Rights & Permissions [Opens in a new window]

Abstract

Guidelines for a healthy diet aim to decrease the risk of chronic diseases. It is unclear as to what extent a healthy diet is also an environmentally friendly diet. In the Dutch sub-cohort of the European Prospective Investigation into Cancer and Nutrition, the diet was assessed with a 178-item FFQ of 40 011 participants aged 20–70 years between 1993 and 1997. The WHO’s Healthy Diet Indicator (HDI), the Dietary Approaches to Stop Hypertension (DASH) score and the Dutch Healthy Diet index 2015 (DHD15-index) were investigated in relation to greenhouse gas (GHG) emissions, land use and all-cause mortality risk. GHG emissions were associated with HDI scores (−3·7 % per sd increase (95 % CI −3·4, −4·0) for men and −1·9 % (95 % CI −0·4, −3·4) for women), with DASH scores in women only (1·1 % per sd increase, 95 % CI 0·9, 1·3) and with DHD15-index scores (−2·5 % per sd increase (95 % CI −2·2, −2·8) for men and −2·0 % (95 % CI −1·9, −2·2) for women). For all indices, higher scores were associated with less land use (ranging from −1·3 to −3·1 %). Mortality risk decreased with increasing scores for all indices. Per sd increase of the indices, hazard ratios for mortality ranged from 0·88 (95 % CI 0·82, 0·95) to 0·96 (95 % CI 0·92, 0·99). Our results showed that adhering to the WHO and Dutch dietary guidelines will lower the risk of all-cause mortality and moderately lower the environmental impact. The DASH diet was associated with lower mortality and land use, but because of high dairy product consumption in the Netherlands it was also associated with higher GHG emissions.

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/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.
Copyright
Copyright © The Authors 2017

Global warming has led to an increased interest in environmentally friendly dietary patterns. At the end of 2015, the Paris Climate Agreement and the United Nations Sustainable Development Goals were initiated( 1 , 2 ). Both agreements reflect the world’s recognition that action is needed. In the European Union, the food sector is responsible for 20–30 % of the total greenhouse gas (GHG) emissions( Reference Tukker, Huppes and Guinée 3 ). A typical Western dietary pattern high in animal products, soft drinks and processed foods is reported to have a large environmental impact( Reference Pimentel and Pimentel 4 ) and is associated with a higher risk of diseases compared with diets rich in vegetables, fruit and fibre-rich cereals( Reference Lock, Pomerleau and Causer 5 Reference Cordain, Eaton and Sebastian 7 ). Correspondingly, a Mediterranean dietary pattern high in fruit and vegetables, legumes, wine, fish and oils and low in meats, dairy products and processed foods is found to be healthier( Reference Serra‐Majem, Roman and Estruch 8 ), and adherence to the Mediterranean guidelines is more environmentally friendly( Reference Sáez-Almendros, Obrador and Bach-Faig 9 ). In a Dutch setting, substituting meat with vegetables, fruit–nuts–seeds, fish or pasta–rice–couscous was associated with both a lower mortality risk (6 to 19 %) and a reduced environmental burden, measured as GHG emissions (4–11 %) and land use (10–12 %)( Reference Biesbroek, Bueno-de-Mesquita and Peeters 10 ). Shifts in dietary patterns can therefore potentially benefit both the environment and health.

As early as 1986, Gussow & Clancy( Reference Gussow and Clancy 11 ) proposed that dietary guidelines should take into account the impact of dietary patterns on global natural resources. Yet, current guidelines are still primarily based on health outcomes. Examples of such guidelines are the World Health Organization( 12 ) and Dutch( Reference Kromhout, Spaaij and de Goede 13 ) dietary guidelines and the Dietary Approaches to Stop Hypertension (DASH) diet( Reference Appel, Moore and Obarzanek 14 ). For research purposes, the levels of adherence to the WHO and DASH dietary guidelines have been operationalised in the Healthy Diet Indicator (HDI) and DASH score, respectively( Reference Huijbregts, Feskens and Rasanen 15 , Reference Fung, Chiuve and McCullough 16 ). Recently, the dietary guidelines of the Netherlands were reviewed( Reference Kromhout, Spaaij and de Goede 13 ), and therefore the Dutch Healthy Diet index, reflecting adherence to the Dutch dietary guidelines of 2006( Reference van Lee, Geelen and van Huysduynen 17 ), was updated to measure adherence to these new guidelines: the Dutch Healthy Diet index 2015 (DHD15-index)( Reference Looman, Feskens and de Rijk 18 ).

Greater adherence to the WHO guidelines for a healthy diet (HDI) has been associated with increased longevity in European and American elderly( Reference Jankovic, Geelen and Streppel 19 ), whereas the DASH diet but not the HDI was significantly associated with a lower risk of developing CVD, CHD and stroke in the European Prospective Investigation into Cancer and Nutrition – Dutch cohort (EPIC-NL) study( Reference Struijk, May and Wezenbeek 20 ). The previous DHD-index was not associated with disease burden (disability-adjusted life years) or CVD risk in one Dutch cohort( Reference Struijk, May and Wezenbeek 20 , Reference Struijk, Beulens and May 21 ) but was associated with decreased risk of all-cause mortality and CVD mortality in another( Reference van Lee, Geelen and Kiefte-de Jong 22 ).

To the best of our knowledge, no prospective study has compared the effect of adherence to dietary guidelines on indicators of environmental impact and health outcomes. Therefore, in this study, we first quantified the association between better adherence to these guidelines and dietary environmental impact. Second, we studied the association between adherence to the dietary guidelines and risk of all-cause mortality. We used data from the Dutch contribution to the EPIC-NL study.

Methods

Study population

EPIC-NL( Reference Beulens, Monninkhof and Verschuren 23 ) consists of 40 011 subjects of EPIC-Prospect( Reference Boker, Van Noord and Van Der Schouw 24 ) and EPIC-MORGEN( Reference Verschuren, Blokstra and Picavet 25 , Reference Blokstra, Smit and Bueno de Mesquita 26 ), both carried out between 1993 and 1997.The EPIC-Prospect cohort included 17 357 women aged 49–70 years living in the city of Utrecht and its vicinity. The EPIC-MORGEN cohort included 22 654 men and women aged 20–65 years, living in Amsterdam, Maastricht and Doetinchem. This study was conducted according to the guidelines laid down in the Declaration of Helsinki, and all procedures involving human participants were approved by the Institutional Review Board of the University Medical Centre Utrecht and the Medical Ethical Committee of TNO Nutrition and Food Research. Written informed consent was obtained from all participants. The design, cohort profile and rationale of EPIC-NL are described elsewhere( Reference Beulens, Monninkhof and Verschuren 23 ). On average, the EPIC-NL cohort had a participation rate of 40 % (35 % in EPIC-Prospect and 45 % in EPIC-MORGEN). The non-response was previously found to have an impact on prevalence estimates of, for example, smoking but not examined associations( Reference Van Loon, Tijhuis and Picavet 27 ).

Diet and environmental impact assessment

Usual daily dietary intake was estimated by a 178-item FFQ, which has been validated against twelve 24-h recalls and biomarkers in 24-h urine and blood( Reference Ocke, Bueno-de-Mesquita and Goddijn 28 , Reference Ocke and Bueno-de-Mesquita 29 ). Spearman’s rank correlation coefficients based on estimates of the FFQ and 24-h recalls in men were 0·51 for potatoes, 0·36 for vegetables, 0·68 for fruits, 0·39 for meat, 0·69 for dairy products, 0·76 for sugar and sweet products and 0·52 for biscuits and pastry. Results for women were similar. Energy intake and daily nutrient intakes were estimated using the 1996 Dutch Food Composition table( 30 ). Blonk Consultants assessed the environmental impact of food items consumed by the Dutch population( Reference Marinussen, Kramer and Pluimers 31 ). Environmental impact was calculated based on life cycle assessments (LCA). The LCA were cradle to grave and included all steps from production, transport, preparation, to waste. The impact value of a food item is a weighted average of different subtypes of the product, for example, by country of origin, which is based on the Dutch production (for the Dutch market) and Dutch import data. GHG emissions are expressed as kg CO2 equivalents/d. Land use is expressed as m2×year/d. These LCA data were combined with the EPIC-NL FFQ data to calculate daily GHG emissions and land use associated with the usual diet. For a more elaborate description of the calculation of the environmental impact of the diet, see our previous paper( Reference Biesbroek, Bueno-de-Mesquita and Peeters 10 ).

Dietary indices

Healthy Diet Indicator

The HDI is based on the 2002 WHO guidelines( 12 ) and has previously been used in other EPIC-NL papers( Reference Struijk, May and Wezenbeek 20 , Reference Struijk, Beulens and May 21 , Reference Fransen, Beulens and May 32 , Reference Berentzen, Beulens and Hoevenaar-Blom 33 ). The HDI score consists of six nutrients (SFA, PUFA, cholesterol, protein, dietary fibre and free sugars) and one food group (fruits and vegetables) (Table 1). When the intake was within the recommended range according to WHO’s guidelines, a score of 1 was assigned to that component; otherwise, 0 points were given. The final HDI score was the sum of all these components, ranging from 0 (minimal adherence) to 7 (maximal adherence).

Table 1 Components and scoring criteria of the indices measuring adherence to dietary guidelines

En%, percentage of total energy intake (excluding alcohol); DASH, Dietary Approaches to Stop Hypertension.

* Higher quintile represents higher intake. Scoring of the components of the DASH score depends on the sex-specific quintile (1, 2, 3, 4 or 5 points).

A score above the recommended intake is 10 points, whereas an intake below is given a proportional score between 0 and 10 points.

A maximum of 40 g cheese/d could be included.

§ A maximum of 4 g lean fish/d could be included.

Dietary Approaches To Stop Hypertension score

The DASH score is based on eight criteria from the DASH clinical trial from 1997( Reference Appel, Moore and Obarzanek 14 ) and has been used in other analyses in EPIC-NL( Reference Struijk, May and Wezenbeek 20 , Reference Joosten, Grobbee and Verschuren 34 ). Included components are fruit, vegetables, nuts and legumes, whole grains, low-fat dairy products, Na, red and processed meats and sweetened beverages. For each of the components, participants were classified into sex-specific quintiles according to their intake. A score ranging from 1 to 5 was given to each quintile. For all components except Na, red and processed meats and sweetened beverages, higher intakes were given higher scores, whereas for the latter components higher intakes were given lower scores (Table 1). In our population, only Na incorporated in food products is accounted for and added salt during cooking or at the table is not included. The overall DASH score is the sum of all components and can range from 8 (lowest adherence) to 40 (highest adherence).

Dutch Healthy Diet index 2015

The DHD15-index is based on the Dutch dietary guidelines of 2015( Reference Kromhout, Spaaij and de Goede 13 ) and is an updated version of the previous DHD-index presented by Van Lee et al. ( Reference van Lee, Geelen and van Huysduynen 17 ). The DHD15-index consists of fifteen components (Table 1). A proportional score between 0 and 10 was assigned to the components. There are five component types included, which are adequacy, moderation, optimum, ratio and quality components. For adequacy components, the recommendation is to consume at least the mentioned quantity. The moderation components represent the foods for which the intake should be lowered. The optimum component has an optimal range of intake. The ratio components are based on replacement of one group by another food group. The quality component is based on the type of food group. The components vegetables, fruit, legumes, nuts, fish and tea are adequacy components, and the components red meat, processed meat, sweetened beverages and fruit juices, alcohol and Naare moderation components. The component dairy product is an optimum component with an optimal range of intake, whereas the fats and oils component is defined as a ratio component. The coffee component is defined as a qualitative component based on type of coffee (filtered or unfiltered). The whole-grain component is scored with two subcomponents: an adequacy component for whole-grain consumption and a ratio component for the ratio whole-grain products and refined grain products.

More detailed information on the calculation of the DHD15-index will be published in another paper( Reference Looman, Feskens and de Rijk 18 ). In short, for the adequacy components, a lower limit of intake was given. For example, for fruits, it is recommended to eat at least 200 g/d. This level of intake received the maximum score of 10 points, gradually decreasing to an intake level of 0, which received 0 points. For the moderation components, an upper limit was given in the guidelines. For example, for red meat it was recommended to eat a maximum of 45 g/d. This level of intake received 10 points, decreasing to a score of 0 points at intakes of 100 g/d or more. The optimum component (dairy product) was calculated by assigning 10 points when the intake was within the optimum range (300–450 g/d). When intake was below the optimum range, the scores decreased linearly with lower intakes. When intake was above the optimum range, higher intakes were given linearly fewer points with a score of 0 points assigned to intakes of 750 g/d or more. The ratio components were scored by calculating the ratio between the recommended food group and the food group that needed to be replaced, and dividing this ratio by the difference between threshold and cut-off value. The maximum score of 10 points was assigned when the ratio was higher than the cut-off value, and gradually decreased to 0 points at the threshold value. In addition, for some foods a quality aspect was defined (coffee); the maximum score of 10 points was assigned if all coffee consumed was filtered or if there was no coffee consumption. If all coffee consumed was unfiltered, the score would be 0. The final DHD15-index was the sum of all components and ranged from 0 (minimal adherence) to 140 (maximal adherence).

All-cause mortality assessment

The vital status of all EPIC-NL participants was obtained through linkage with the municipal population registries. Participants were followed up over time until death by any cause, loss to follow-up or were censored on 1 January 2015. During a mean follow-up of 19·2 (sd 3·3) years, 3845 deaths were documented.

Lifestyle and anthropometric variables

At baseline, study participants completed a standardised structured general questionnaire on the presence of chronic diseases, related potential risk factors and lifestyle factors. Blood pressure, weight and height were measured by trained staff according to standardised protocols( Reference Beulens, Monninkhof and Verschuren 23 ). BMI was calculated by dividing weight by height squared (kg/m2). Physical activity was assessed with a validated questionnaire( Reference Haftenberger, Schuit and Tormo 35 ) and classified according to the Cambridge Physical Activity Index (CPAI) with imputed data for missing values (n 4930)( Reference Pols, Peeters and Ocke 36 ). The CPAI is categorised into inactive, moderately inactive, moderately active and active. Smoking was operationalised as current, former and never smoker. Educational level was coded as low (lower vocational training or primary school), medium (intermediate vocational training or secondary school) or high (higher vocational training or university).

Statistical analysis

Participants without dietary information at baseline were excluded (n 218) from this study. Participants with implausible dietary intake – that is, those in the highest and lowest 0·5 % of the ratio of reported energy intake:BMR – were also excluded (n 400). Participants without informed consent for linkage to municipal registries were excluded (n 1034). Participants with a self-reported history of cancer (n 1616), diabetes (n 759), myocardial infarction (n 514) or stroke (n 451) at baseline were excluded because their reported usual diet may not reflect their diet before diagnosis. Participants with missing information on BMI (n 17), educational level (n 191) or smoking status (n 23) were also excluded. After these exclusions, 35 031 participants remained for analysis (EPIC-Prospect=14 770 and EPIC-MORGEN=20 261). As these are secondary analyses based on an existing large cohort with a long follow-up time, the justification for the sample size is not required.

Because the diets of men and women differ with respect to total energy intake and environmental impact, all analyses were stratified by sex. For the DASH score and DHD15-index, tertiles of adherence were created. Because the HDI ranged from 0 to 7 only, the HDI was categorised into three groups with a score of 0–2, 3 and 4–7, respectively. Tabulations of sociodemographic data by HDI categories and tertiles of the DASH score and DHD15-index were made.

General linear models were used to calculate differences in mean GHG emission and land use in the different categories of dietary indices. The first category/tertile was used as reference. These analyses were adjusted for age at baseline, total energy intake and physical activity to compare the environmental impact based on dietary choices independently of total amount of foods consumed. In addition, we added educational level to the model in a sensitivity analysis. We assessed multicollinearity using the variance inflation factor, homoscedasticity and independence of the residuals, the mean and normality of the residuals and checked for linearity and found that all criteria were met.

Cox proportional hazard models were used to estimate hazard ratios (HR) with 95 % CI for the associations between dietary indices and all-cause mortality. The first category/tertile was used as reference. The model was pooled for sub-cohort (EPIC-MORGEN or EPIC-Prospect) and adjusted for confounding by age, BMI, educational level, smoking status, total energy intake and physical activity. The HDI and DASH score did not include alcohol consumption; these models were also adjusted for this variable. The proportional hazards assumption was checked using the Schoenfeld residuals test and showed that none of the P values were significant. P values for the linear trend across the categories were calculated by including the mean score of each tertile as continuous variable in the model.

All analyses were repeated for a (continuous) change of 1 sd in the dietary score to enable a better comparison of the associations between the different dietary scores.

All statistical analyses were performed using SAS software (version 9.4; SAS Institute Inc.). A two-sided P<0·05 was considered statistically significant.

Results

Participants with higher DASH and DHD15-index scores tend to be older, whereas those with a higher HDI score are younger than participants with lower scores (Table 2). The BMI and prevalence of current smoking is consistently lower with higher dietary index scores, whereas physical activity levels are higher. With higher dietary index scores, energy intake is lower for the DASH score and DHD15-index but higher for the HDI. Alcohol intake is lower with higher dietary index scores, with the exception of women for whom alcohol intake remains constant across tertiles of the DASH score. Mean dietary index component scores and percentages of participants meeting the guidelines are presented in the online Supplementary Table S1.

Table 2 Baseline characteristics of the Dutch sub-cohort of the European Investigation into Cancer and Nutrition according to tertiles of the Healthy Diet Indicator (HDI), Dietary Approaches to Stop Hypertension (DASH) score and Dutch Healthy Diet index 2015 (DHD15-index) (Mean values and standard deviations; percentages)

Adherence to dietary guidelines and environmental impact

On average, men’s diets have a higher total environmental impact than women’s diet (respectively, 4·6 kg CO2-eq/d and 4·4 m2×year/d v. 3·7 kg CO2-eq/d 3·5 m2×year/d). However, when expressed per 4184 kJ (1000 kcal), women’s diets have higher GHG emissions and land use (1·8 kg CO2-eq and 1·7 m2×year v. 2·0 kg CO2-eq and 1·9 m2 year).

Mean GHG emissions and land uses according to the three categories of adherence to the three dietary indices are presented in Fig. 1. In men, comparing the highest category with the lowest, dietary GHG emissions are significantly lower for the HDI (−9·1 %; 95 % CI −8·4, −9·9) and for the DHD15-index (−5·5 %; 95 % CI −4·7, −6·3) but not for the DASH score (0·6; 95 % CI −0·2, 1·4) after adjusting for age, energy intake and physical activity levels (Table 3). In women, better adherence to the guidelines is associated with statistically lower GHG emissions for the HDI and DHD15-index. However, higher scores on the DASH diet were associated with significantly higher GHG emissions (2·3 %; 95 % CI 1·8, 2·9). Analysing the results continuously per difference of 1 sdof the score gave very similar results. Including educational level as a possible additional confounder did not change the results (not shown).

Fig. 1 Greenhouse gas (GHG) emissions (a) and land use (b) according to tertiles of the Healthy Diet Indicator (HDI), Dietary Approaches to Stop Hypertension (DASH) diet and Dutch Healthy Diet index 2015 (DHD15-index). Values are adjusted means with their standard errors. , Category/tertile 1; , category/tertile 2; , category/tertile 3. All values are adjusted for age at baseline, energy intake and physical activity level. Significance compared with category/tertile 1. * P<0·05, ** P<0·0001.

Table 3 Greenhouse gas (GHG) emissions and land use according to the Healthy Diet Indicator (HDI), Dietary Approaches to Stop Hypertension (DASH) diet and Dutch Healthy Diet index 2015 (DHD15-index)Footnote * (Adjusted differences and 95 % confidence intervals; mean values and standard deviations)

Ref., referent values.

* All values adjusted for age at baseline, energy intake and physical activity level.

All dietary guideline indices showed that an increase in the score is associated with lower land use (Fig. 1 and Table 3). For men, comparing category 3 with category 1, land use is significantly lower by −7·7 % (95 % CI −6·9, −8·5) for the HDI, −2·7 % (95 % CI −1·9, −3·5) for the DASH score and −7·1 % (95 % CI −6·3, −7·9) for the DHD15-index. In women, land use of category 3 v. category 1 is−5·0 % (95 % CI −4·5, −5·6) lower for the HDI, −3·5 % (95 % CI −2·9, −4·0) for the DASH score and −9·6 % (95 % CI −9·1, −10·1) for the DHD15-index.

Adherence to dietary guidelines and all-cause mortality

Mortality risk is significantly lower in the highest compared with the lowest category of the HDI for both men (HRC3–C1 0·82; 95 % CI 0·70, 0·97) and women (HRC3–C1 0·83; 95 % CI 0·76, 0·91) (Table 4). Adherence to the DASH diet is not associated with all-cause mortality in men and women when analysed in tertiles. However, analysing the DASH score continuously per sd increase of the score, a better adherence is significantly associated with lower risk of all-cause mortality (Table 4). Mortality risk is significantly lower with better adherence to the DHD15-index in men (HRT3–T1 0·84; 95 % CI 0·69, 0·98) and women (0·85; 95 % CI 0·78, 0·96).

Table 4 Associations between the Healthy Diet Indicator (HDI), Dietary Approaches to Stop Hypertension (DASH) diet and Dutch Healthy Diet index 2015 (DHD15-index), and all-cause mortality among 35 031 European Prospective Investigation into Cancer and Nutrition – Dutch cohort (EPIC-NL) participantsFootnote * (Adjusted hazard ratios (HR) and 95 % confidence intervals; medians, mean values and standard deviations)

Ref., referent values.

Model is pooled for cohort (EPIC-Prospect or EPIC-MORGEN).

* Mortality risk: adjusted for age at baseline, BMI, educational level, smoking status, total daily energy intake and physical activity level. In addition, the HDI and DASH models are adjusted for alcohol intake.

Discussion

Our study shows that better adherence to the dietary guidelines from the WHO, DASH (in men only) and Dutch Health Council is associated with lower environmental impact and lower risk of all-cause mortality. In men, the largest difference in environmental impact is observed for higher scores on the HDI (−3·7 % per sd for GHG emissions and −3·3 % per sd for land use). For women, the largest environmental impact differences are observed for the DHD15-index (−2·0 % per sd for GHG emission and −3·1 % per sd for land use). Higher DHD15-index scores are associated with the lowest relative all-cause mortality risk of the indices (HR sd of 0·88 for men and 0·92 for women).

The different guidelines – and therefore the indices – have conceptual differences. Despite these differences, the HDI, DASH and DHD15-index have quite comparable associations with all-cause mortality, and the HDI and DHD15-index show comparable possible reductions in GHG emissions and land use of the diet in our study.

Government and health organisations promote dietary patterns linked to a broad range of positive health effects. Considering the significant environmental impact of our current diet, these diets should ideally be accompanied by lower GHG emissions and land use to meet the Food and Agriculture Organization’s definition of a sustainable diet( 37 ). In a literature review of sixteen studies by Payne et al., it is stated that dietary patterns that primarily aim to reduce GHG emissions may not always improve nutritional quality or health outcomes compared with the average dietary patterns( Reference Payne, Scarborough and Cobiac 38 ). Similarly, in a previous paper, we showed that total GHG emissions and land use of the diet were not associated with all-cause mortality( Reference Biesbroek, Bueno-de-Mesquita and Peeters 10 ). On the other hand, our current results suggest that adhering to some of the dietary guidelines will both improve health and moderately reduce environmental impact. Although diets according to guidelines most likely ensure a higher dietary quality, their links with environmental impact are less clear. A review of modelling studies showed that meeting dietary guidelines may reduce GHG emissions by 0–35 % and land use by 15–50 % compared with the average observed food consumption in a population( Reference Hallström, Carlsson-Kanyama and Börjesson 39 ). However, in five of the fourteen scenarios the reduction potential was <10 %. Although the increased consumption of some food groups, such as fruits and vegetables, legumes and nuts, increases the environmental impact, this is usually outweighed by the lower consumption of meat and products, such as snacks, sweets and pastries( Reference Van Dooren, Marinussen and Blonk 40 ). In addition, fish consumption would have to increase to meet the recommendation, which will not only increase GHG emissions and land use (farmed fish), but it will also put pressure on wild fish stocks. However, total GHG emissions will not necessarily increase when fish replaces other protein-rich foods with a higher environmental impact, such as beef( Reference Seves, Temme and Brosens 41 ). This underlines the importance of looking at dietary patterns from both a health and an environmental perspective together to create environmentally friendly and healthy diets.

Although limited to two studies, environmental impact has been studied for the HDI( Reference Green, Milner and Dangour 42 ) and the DASH diet( Reference Monsivais, Scarborough and Lloyd 43 ) before. Green et al. ( Reference Green, Milner and Dangour 42 ) modelled and optimised the current UK consumption to completely match the WHO guidelines and showed a possible 17 % reduction in GHG emissions. The difference found by the authors was larger than we observed in our study (category 3 (4–7 points) v. category 1 (0–2 points): 4 % lower GHG emissions in women and 9 % in men), but we compared categories of adherence and did not model the impact when all recommendations are completely met (mean score was 4·5 out of a possible 7 in the highest category). In addition, this was an optimisation study in which specific foods with a low environmental impact were selected to replace unhealthy foods instead of observing actual food intakes as in our study. Monsivais et al. ( Reference Monsivais, Scarborough and Lloyd 43 ) studied the GHG emissions associated with the DASH diet within the UK population of the EPIC cohort. GHG emissions in the highest quintile of adherence were 17 % lower than in the lowest quintile. We observed a small increase in GHG emissions (a not significant increase of 0·6 % for men and a significant 2·3 % increase for women) between the lowest to the highest tertile of the DASH score. There is a noticeable difference between our studies in the calculation of the DASH score. We analyse the foods in g/d, whereas Monsivais et al. used food groups expressed as energy percentage. Low-fat dairy product consumption is promoted in the DASH diet, but the mean consumption and the variation in consumption in the UK cohort are much lower than in our population (174·9 (sd 104·0)( Reference O’Connor, Lentjes and Luben 44 ) v. 250·8 (sd 209·1) g/d). As the DASH score is based on quintiles of intake, the difference between the consumption in Q5 and Q1 of low-fat dairy product is much larger in our population than in the UK one. Consequently, the difference in environmental impact of this food group is also much larger because dairy product is a major contributor to total GHG emissions. Total dairy product is on average responsible for 25 % of the diet-related GHG emissions in our population( Reference Biesbroek, Bueno-de-Mesquita and Peeters 10 ). Also, differences between quintiles of the total DASH score (Monsivais et al. ( Reference Monsivais, Scarborough and Lloyd 43 )) are expected to be larger than between tertiles (our study). Combined, these factors may account for the small increase in GHG emissions in our cohort compared with the decrease observed in the UK cohort.

Previously, higher levels of adherence to the HDI were inversely associated with all-cause mortality risk in two European population studies with similar risk estimates as we present here( Reference Jankovic, Geelen and Streppel 19 , Reference Huijbregts, Feskens and Räsänen 45 ). The DASH diet was associated with lower all-cause mortality in older adults( Reference Reedy, Krebs-Smith and Miller 46 ) and adults with hypertension( Reference Parikh, Lipsitz and Natarajan 47 ), similar to our results. For the previous DHD-index, based on the Dutch guidelines from 2006, supporting evidence of the overall association with all-cause and cause-specific mortality is inconclusive( Reference Struijk, May and Wezenbeek 20 Reference van Lee, Geelen and Kiefte-de Jong 22 ), whereas the new DHD15-index is clearly associated with mortality in our population.

Energy intake is correlated with GHG emissions. In a UK population, for every 4184 kJ (1000 kcal) GHG emissions increased by, on average, 3 kg CO2-eq( Reference Monsivais, Scarborough and Lloyd 43 ). Therefore, we adjusted our analysis for energy intake to independently study the effect of dietary quality on GHG emissions and land use. We can conclude that at equal energy intake adhering to the HDI, DASH (only for land use) and the DHD15-index is better for the environment. Taking the obesity trend in the Netherlands into consideration, not only dietary quality should be improved, but also the limiting of energy intake should be crucial. Besides reducing the burden of disease associated with overweight and obesity( Reference Guh, Zhang and Bansback 48 ), this would also reduce the environmental impact of our diet by less food being eaten. A modelling scenario in the US in which the energy intake was reduced to maintain a healthy body weight without changing the actual food mix resulted in an approximate reduction in GHG emissions of 9 %( Reference Tom, Fischbeck and Hendrickson 49 ).

The reductions in GHG emissions that are observed in our study are only moderate, but are accompanied by a clear reduction in mortality risk. The added value of sustainability aspects in nutritional guidelines would be that consumers would learn what the healthy foods are and could combine this with the sustainable choice between food groups and also within each food group( Reference Gonzalez Fischer and Garnett 50 ). However, to achieve a more sustainable food system, change should not be limited to the choices made by consumers. Within each food group, producers, retailers and transport businesses may invest in new technologies that retain crop yield while protecting biodiversity, reuse materials, provide better storage and produce less waste to facilitate a more sustainable production of foods. This is why both aspects are combined in the UN Sustainable Development Goals under the header ‘responsible consumption and production’( 2 ). However, if the consumers’ demand for sustainable and healthy foods becomes more eminent, this might force companies to follow.

Major strengths of this study are that we used both dietary and environmental impact data of the same population and have linked these with registered mortality data. In addition, we have a large population-based cohort with a long follow-up of 19 years. We compared three different indices for healthy diets and found similar and thus robust results. Some limitations need to be addressed. Our study assessed dietary intake and its environmental impact only at baseline and only in adults. We assume in our current analyses that both intake and impact are stable. However, a previous study compared the environmental impact of the Dutch diet of 2007/2010 with that of 1997/1998 and observed a 4·9 % lower GHG emission in men and 7·0 % for women because of changing dietary intakes( Reference Marinussen, Kramer and Pluimers 31 ). According to the Dutch National Food Consumption Survey, between 2012 and 2014 and 2007 and 2010, the average diet of Dutch adults showed a decrease in the consumption of potatoes, fats and oils, alcoholic beverages, dairy products, cakes and biscuits and meat v. an increase in the intake of non-alcoholic drinks and condiments and sauces. Adolescents showed similar dietary changes, but in addition fruit intake increased by 20 %. Not taking such dietary changes into account can result in an underestimation or overestimation of the association between the dietary patterns and mortality because participants can be misclassified. If we were to apply the Dutch guidelines to the current food consumption data, most Dutch people would not meet the guidelines at this moment( Reference Rossum, Buurma-Rethans and Vennemann 51 ) and thus our message that dietary change is needed to increase sustainability of the diet and health remains a priority.

Under- and over-reporting of dietary energy intake might affect the associations of the dietary patterns with environmental impact. Therefore, we excluded participants in the highest and lowest 0·5 % of the ratio of reported energy intake:BMR as proxy for over- and under-reporters. The indicators for environmental impact are based on Dutch LCA data and apply to a Dutch setting only. Our results for environmental impact of the dietary patterns may therefore not be directly extrapolated to other countries in which production methods, productivity, fossil energy use, import and export and ways of transport may differ.

In conclusion, national and international guidelines for a healthy diet are aimed at decreasing the risk of chronic diseases. If these guidelines were adhered to by a larger proportion of the Dutch population, lower risk of all-cause mortality, as well as reductions in GHG emissions and land use, could be achieved. The possible reductions in GHG emissions and land use seem to be moderate. Eating more plant-based instead of animal-based products and, according to the guidelines for a healthy diet, limiting energy intake to match energy requirement are all strategies that need to be combined and applied to maximise the health potential while limiting the environmental impact of our diet.

Acknowledgements

The authors would like to thank Rianne Weggemans and Caroline Spaaij from the Dutch Health Council for their help in interpreting the Dutch dietary guidelines in order to create the DHD15-index. The authors would also like to thank Edith Feskens from Wageningen University for her help and expertise in creating and calculating the DHD15-index.

This project was funded by a research grant from the Strategic Program of the RIVM. Project code: S133006.

S. B., W. M. M. W. and E. H. M. T. designed the research; S. B. conducted research and analysed the data; S. B., A. G., M. L. and E. H. M. T. designed the DHD15-index. J. M. A. B., M. E. v. d. K. and Y. T. v. d. S. provided valuable input to the manuscript. All authors read and approved the final manuscript.

The authors declare that there are no conflicts of interest.

Supplementary material

For supplementary material/s referred to in this article, please visit https://doi.org/10.1017/S0007114517001878

References

1. United Nations (2015) Paris agreement. https://treaties.un.org/doc/Treaties/2016/02/20160215%2006-03%20PM/Ch_XXVII-7-d.pdf (accessed September 2016).Google Scholar
2. United Nations General Assembly (2015) Resolution 70/1 transforming our world: the 2030 agenda for sustainable development. https://documents-dds-ny.un.org/doc/UNDOC/GEN/N15/291/89/pdf/N1529189.pdf?OpenElement (accessed September 2016).Google Scholar
3. Tukker, A, Huppes, G, Guinée, J, et al. (2006) Environmental Impact of Products (EIPRO) analysis of the life cycle environmental impacts related to the final consumption of the EU-25, Report EUR 22284 EN. Brussels, Belgium.Google Scholar
4. Pimentel, D & Pimentel, M (2003) Sustainability of meat-based and plant-based diets and the environment. Am J Clin Nutr 78, 660S663S.Google Scholar
5. Lock, K, Pomerleau, J, Causer, L, et al. (2005) The global burden of disease attributable to low consumption of fruit and vegetables: implications for the global strategy on diet. Bull World Health Organ 83, 100108.Google Scholar
6. Sinha, R, Cross, AJ, Graubard, BI, et al. (2009) Meat intake and mortality: a prospective study of over half a million people. Arch Intern Med 169, 562571.Google Scholar
7. Cordain, L, Eaton, SB, Sebastian, A, et al. (2005) Origins and evolution of the Western diet: health implications for the 21st century. Am J Clin Nutr 81, 341354.Google Scholar
8. Serra‐Majem, L, Roman, B & Estruch, R (2006) Scientific evidence of interventions using the Mediterranean diet: a systematic review. Nutr Rev 64, s1, S27S47.CrossRefGoogle ScholarPubMed
9. Sáez-Almendros, S, Obrador, B, Bach-Faig, A, et al. (2013) Environmental footprints of Mediterranean versus Western dietary patterns: beyond the health benefits of the Mediterranean diet. Environ Health 12, 118.Google Scholar
10. Biesbroek, S, Bueno-de-Mesquita, HB, Peeters, PH, et al. (2014) Reducing our environmental footprint and improving our health: greenhouse gas emission and land use of usual diet and mortality in EPIC-NL – a prospective cohort study. Environ Health 13, 27.Google Scholar
11. Gussow, JD & Clancy, KL (1986) Dietary guidelines for sustainability. J Nutr Educ 118, 15.Google Scholar
12. World Health Organization (2003) Diet, Nutrition and the Prevention of Chronic Diseases. Geneva: WHO.Google Scholar
13. Kromhout, D, Spaaij, C, de Goede, J, et al. (2016) The 2015 Dutch food-based dietary guidelines. Eur J Clin Nutr 70, 869878.CrossRefGoogle ScholarPubMed
14. Appel, LJ, Moore, TJ, Obarzanek, E, et al. (1997) A clinical trial of the effects of dietary patterns on blood pressure. N Eng J Med 336, 11171124.Google Scholar
15. Huijbregts, P, Feskens, E, Rasanen, L, et al. (1997) Dietary pattern and 20-year mortality in elderly men in Finland, Italy and The Netherlands: longitudinal cohort study. BMJ 315, 1317.Google Scholar
16. Fung, TT, Chiuve, SE, McCullough, ML, et al. (2008) Adherence to a DASH-style diet and risk of coronary heart disease and stroke in women. Arch Intern Med 168, 713720.Google Scholar
17. van Lee, L, Geelen, A, van Huysduynen, E, et al. (2012) The Dutch Healthy Diet index (DHD-index): an instrument to measure adherence to the Dutch Guidelines for a Healthy Diet. Nutr J 11, 19.Google Scholar
18. Looman, M, Feskens, EJ, de Rijk, M, et al. (2017) Development and evaluation of the Dutch Healthy Diet index 2015. Public Health Nutr (Epublication ahead of print version 19 June 2017).Google Scholar
19. Jankovic, N, Geelen, A, Streppel, MT, et al. (2014) Adherence to a healthy diet according to the World Health Organisation guidelines and all-cause mortality in elderly adults from Europe and the United States. Am J Epidemiol 180, 978988.Google Scholar
20. Struijk, EA, May, AM, Wezenbeek, NL, et al. (2014) Adherence to dietary guidelines and cardiovascular disease risk in the EPIC-NL cohort. Int J Cardiol 176, 354359.Google Scholar
21. Struijk, EA, Beulens, JW, May, AM, et al. (2014) Dietary patterns in relation to disease burden expressed in disability-adjusted life years. Am J Clin Nutr 100, 11581165.Google Scholar
22. van Lee, L, Geelen, A, Kiefte-de Jong, JC, et al. (2016) Adherence to the Dutch dietary guidelines is inversely associated with 20-year mortality in a large prospective cohort study. Eur J Clin Nutr 70, 262268.CrossRefGoogle Scholar
23. Beulens, JWJ, Monninkhof, EM, Verschuren, WMM, et al. (2010) Cohort profile: the EPIC-NL study. Int J Epidemiol 39, 11701178.Google Scholar
24. Boker, LK, Van Noord, P, Van Der Schouw, Y, et al. (2001) Prospect-EPIC Utrecht: study design and characteristics of the cohort population. Eur J Epidemiol 17, 10471053.Google Scholar
25. Verschuren, W, Blokstra, A, Picavet, H, et al. (2008) Cohort profile: the Doetinchem cohort study. Int J Epidemiol 37, 12361241.Google Scholar
26. Blokstra, A, Smit, H, Bueno de Mesquita, H, et al. (2005) Monitoring of risk factors and health in the Netherlands (MORGEN-cohort) 1993–1997. Lifestyle- and risk factors: prevalences and trends, Report no. 263200008. Bilthoven: RIVM (in Dutch).Google Scholar
27. Van Loon, AJM, Tijhuis, M, Picavet, HSJ, et al. (2003) Survey non-response in the Netherlands: effects on prevalence estimates and associations. Ann Epidemiol 13, 105110.Google Scholar
28. Ocke, MC, Bueno-de-Mesquita, HB, Goddijn, HE, et al. (1997) The Dutch EPIC food frequency questionnaire. I. Description of the questionnaire, and relative validity and reproducibility for food groups. Int J Epidemiol 26, Suppl. 1, S37.Google Scholar
29. Ocke, MC, Bueno-de-Mesquita, HB, et al. (1997) The Dutch EPIC food frequency questionnaire. II. Relative validity and reproducibility for nutrients. Int J Epidemiol 26, Suppl. 1, S49.Google Scholar
30. NEVO-table (1996) Dutch Food Composition Table. Zeist: NEVO (Nederlandse Voedingsstoffenbestand) Foundation (in Dutch).Google Scholar
31. Marinussen, M, Kramer, G, Pluimers, J, et al. (2012) The Environmental Impact of Our Diet: An Analysis Based on de Nutitional Consumption Survey of 2007–2010, (summary in English). Breda: Blonk Milieu Advies (in Dutch, summary in English).Google Scholar
32. Fransen, HP, Beulens, JW, May, AM, et al. (2015) Dietary patterns in relation to quality-adjusted life years in the EPIC-NL cohort. Prev Med 77, 119124.CrossRefGoogle ScholarPubMed
33. Berentzen, NE, Beulens, JW, Hoevenaar-Blom, MP, et al. (2013) Adherence to the WHO’s Healthy Diet Indicator and overall cancer risk in the EPIC-NL cohort. PLOS ONE 8, e70535.Google Scholar
34. Joosten, MM, Grobbee, DE, Verschuren, WM, et al. (2010) Combined effect of alcohol consumption and lifestyle behaviors on risk of type 2 diabetes. Am J Clin Nutr 91, 17771783.Google Scholar
35. Haftenberger, M, Schuit, A, Tormo, M, et al. (2002) Physical activity of subjects aged 50–64 years involved in the European Prospective Investigation into Cancer and Nutrition (EPIC). Public Health Nutr 5, 11631178.Google Scholar
36. Pols, MA, Peeters, P, Ocke, MC, et al. (1997) Estimation of reproducibility and relative validity of the questions included in the EPIC Physical Activity Questionnaire. Int J Epidemiol 26, Suppl. 1, S181.Google Scholar
37. Food and Agriculture Organization (2010) Biodiversity and Sustainable Diets United Against Hunger. International Scientific Symposium. Rome: FAO.Google Scholar
38. Payne, CL, Scarborough, P & Cobiac, L (2016) Do low-carbon-emission diets lead to higher nutritional quality and positive health outcomes? A systematic review of the literature. Public Health Nutr 19, 26542661.CrossRefGoogle ScholarPubMed
39. Hallström, E, Carlsson-Kanyama, A & Börjesson, P (2015) Environmental impact of dietary change: a systematic review. J Clean Prod 91, 111.Google Scholar
40. Van Dooren, C, Marinussen, M, Blonk, H, et al. (2014) Exploring dietary guidelines based on ecological and nutritional values: a comparison of six dietary patterns. Food Policy 44, 3646.Google Scholar
41. Seves, SM, Temme, EHM, Brosens, MCC, et al. (2016) Sustainability aspects and nutritional composition of fish: evaluation of wild and cultivated fish species consumed in the Netherlands. Climatic Change 135, 597610.Google Scholar
42. Green, R, Milner, J, Dangour, A, et al. (2015) The potential to reduce greenhouse gas emissions in the UK through healthy and realistic dietary change. Clim Change 129, 253265.Google Scholar
43. Monsivais, P, Scarborough, P, Lloyd, T, et al. (2015) Greater accordance with the DASH dietary pattern is associated with lower diet-related greenhouse gas production but higher dietary costs in the United Kingdom. Am J Clin Nutr 102, 138145.Google Scholar
44. O’Connor, LM, Lentjes, MAH, Luben, RN, et al. (2014) Dietary dairy product intake and incident type 2 diabetes: a prospective study using dietary data from a 7-day food diary. Diabetologia 57, 909917.Google Scholar
45. Huijbregts, P, Feskens, E, Räsänen, L, et al. (1997) Dietary pattern and 20 year mortality in elderly men in Finland, Italy, and the Netherlands: longitudinal cohort study. BMJ 315, 1317.Google Scholar
46. Reedy, J, Krebs-Smith, SM, Miller, PE, et al. (2014) Higher diet quality is associated with decreased risk of all-cause, cardiovascular disease, and cancer mortality among older adults. J Nutr 144, 881889.Google Scholar
47. Parikh, A, Lipsitz, SR & Natarajan, S (2009) Association between a DASH-like diet and mortality in adults with hypertension: findings from a population-based follow-up study. Am J Hypertens 22, 409416.Google Scholar
48. Guh, DP, Zhang, W, Bansback, N, et al. (2009) The incidence of co-morbidities related to obesity and overweight: a systematic review and meta-analysis. BMC Public Health 9, 88.Google Scholar
49. Tom, MS, Fischbeck, PS & Hendrickson, CT (2016) Energy use, blue water footprint, and greenhouse gas emissions for current food consumption patterns and dietary recommendations in the US. Environ Syst Decis 36, 92103.Google Scholar
50. Gonzalez Fischer, C & Garnett, T (2016) Plates, Pyramids, Planet: Developments in National Healthy and Sustainable Dietary Guidelines: A State of Play Assessment . Rome: The Food and Agriculture Organisation of the United Nations and the Food Climate Research Network at the University of Oxford.Google Scholar
51. Rossum, CTMv, Buurma-Rethans, EJM, Vennemann, FBC, et al. (2016) The diet of the Dutch: results of the first two years of the Dutch National Food Consumption Survey 2012–2016, Letter Report no. 2016-0082. Bilthoven: RIVM.Google Scholar
Figure 0

Table 1 Components and scoring criteria of the indices measuring adherence to dietary guidelines

Figure 1

Table 2 Baseline characteristics of the Dutch sub-cohort of the European Investigation into Cancer and Nutrition according to tertiles of the Healthy Diet Indicator (HDI), Dietary Approaches to Stop Hypertension (DASH) score and Dutch Healthy Diet index 2015 (DHD15-index) (Mean values and standard deviations; percentages)

Figure 2

Fig. 1 Greenhouse gas (GHG) emissions (a) and land use (b) according to tertiles of the Healthy Diet Indicator (HDI), Dietary Approaches to Stop Hypertension (DASH) diet and Dutch Healthy Diet index 2015 (DHD15-index). Values are adjusted means with their standard errors. , Category/tertile 1; , category/tertile 2; , category/tertile 3. All values are adjusted for age at baseline, energy intake and physical activity level. Significance compared with category/tertile 1. * P<0·05, ** P<0·0001.

Figure 3

Table 3 Greenhouse gas (GHG) emissions and land use according to the Healthy Diet Indicator (HDI), Dietary Approaches to Stop Hypertension (DASH) diet and Dutch Healthy Diet index 2015 (DHD15-index)* (Adjusted differences and 95 % confidence intervals; mean values and standard deviations)

Figure 4

Table 4 Associations between the Healthy Diet Indicator (HDI), Dietary Approaches to Stop Hypertension (DASH) diet and Dutch Healthy Diet index 2015 (DHD15-index), and all-cause mortality among 35 031 European Prospective Investigation into Cancer and Nutrition – Dutch cohort (EPIC-NL) participants* (Adjusted hazard ratios (HR) and 95 % confidence intervals; medians, mean values and standard deviations)

Supplementary material: PDF

Biesbroek supplementary material

Biesbroek supplementary material 1

Download Biesbroek supplementary material(PDF)
PDF 271.9 KB
Supplementary material: PDF

Biesbroek supplementary material

Biesbroek supplementary material 2

Download Biesbroek supplementary material(PDF)
PDF 150.3 KB