Next Article in Journal
The Effects of 6 Weeks of Tribulus terrestris L. Supplementation on Body Composition, Hormonal Response, Perceived Exertion, and CrossFit® Performance: A Randomized, Single-Blind, Placebo-Controlled Study
Next Article in Special Issue
Are We What We Eat? Impact of Diet on the Gut–Brain Axis in Parkinson’s Disease
Previous Article in Journal
Grifola frondosa (Maitake) Extract Reduces Fat Accumulation and Improves Health Span in C. elegans through the DAF-16/FOXO and SKN-1/NRF2 Signalling Pathways
Previous Article in Special Issue
Tocotrienols Ameliorate Neurodegeneration and Motor Deficits in the 6-OHDA-Induced Rat Model of Parkinsonism: Behavioural and Immunohistochemistry Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diet Quality and Risk of Parkinson’s Disease: The Rotterdam Study

by
Anne J. Strikwerda
1,
Lisanne J. Dommershuijsen
1,
M. Kamran Ikram
1,2 and
Trudy Voortman
1,3,*
1
Department of Epidemiology, Erasmus MC University Medical Center, P.O. Box 2040, 3000 CA Rotterdam, The Netherlands
2
Department of Neurology, Erasmus MC University Medical Center, P.O. Box 2040, 3000 CA Rotterdam, The Netherlands
3
Division of Human Nutrition and Health, Wageningen University & Research, P.O. Box 9101, 6700 HB Wageningen, The Netherlands
*
Author to whom correspondence should be addressed.
Nutrients 2021, 13(11), 3970; https://doi.org/10.3390/nu13113970
Submission received: 27 September 2021 / Revised: 29 October 2021 / Accepted: 3 November 2021 / Published: 7 November 2021
(This article belongs to the Special Issue Nutrition and Parkinson's Disease)

Abstract

:
The Mediterranean diet has been associated with the risk of Parkinson’s disease (PD), but limited research has been performed on other dietary patterns. We studied the relationship between overall diet quality and PD risk in the general population. We included 9414 participants from the Rotterdam Study, a prospective population-based study in the Netherlands. Diet was defined using a Dutch diet quality score, a Mediterranean diet score and data-driven dietary patterns constructed with principal component analysis (PCA). During an average follow-up of 14.1 years, PD was diagnosed in 129 participants. We identified a ‘Prudent’, ‘Unhealthy’ and ‘Traditional Dutch’ pattern from the PCA. We found a possible association between the Mediterranean diet (Hazard ratio (HR) per standard deviation (SD) 0.89 (95% confidence interval (CI) 0.74–1.07)), the ‘Prudent’ pattern (HR per SD 0.81 (95% CI 0.61–1.08)) and the risk of PD. However, no associations with PD risk were found for the Dutch diet quality score (HR per SD 0.93 (95% CI 0.77–1.12)), the ‘Unhealthy’ pattern (HR per SD 1.05 (95% CI 0.85–1.29)) or the ‘Traditional Dutch’ pattern (HR per SD 0.90 (95% CI 0.69–1.17)). In conclusion, our results corroborate previous findings of a possible protective effect of the Mediterranean diet. Further research is warranted to study the effect of other dietary patterns on PD risk.

1. Introduction

Parkinson’s disease (PD) is estimated to affect over 6 million individuals worldwide; the disease burden is expected to increase even further [1]. The cause of PD is multifactorial, as genetic and environmental factors play an important role [2,3]. Previously, several lifestyle factors had been associated with PD risk [4]. Given the hypothesis that the origin of PD might be found in the gut, and the important influence of diet on the gut microbiome, modification of diet could potentially affect PD risk [5,6,7,8].
Indeed, several studies have linked intake of specific foods or nutrients to PD risk. Frequent intake of certain dairy products has been be associated with an increased risk of PD, whereas higher coffee and tea intake has been associated with a decreased risk of PD [9,10,11,12,13,14]. For other foods, inconsistent results have been described [15,16,17,18]. In addition to food groups, complete dietary patterns have been identified as a comprehensive way to study the role of diet in chronic diseases [19]. Several studies have reported a reduced risk of PD with better adherence to the Mediterranean diet, but inconsistent associations were found between the American Healthy Diet Index and PD risk [20,21,22,23]. The question thus remains whether it is specifically the Mediterranean diet or a ‘healthy’ dietary pattern in general that is associated with PD risk.
In this study, we determined the association between diet and the risk of PD in the general population, through pre-defined diet scores and data-driven dietary patterns. Additionally, we determined the effects of specific food groups on the risk of PD to further unravel the association between diet and PD.

2. Materials and Methods

2.1. Study Population

The Rotterdam Study is a prospective population-based cohort study that aims to examine etiology, natural history and targets for intervention in chronic diseases. The Rotterdam Study was initiated in 1990, recruiting all inhabitants aged 55 years or older from a defined suburban district in Rotterdam (cohort Rotterdam Study-I, RS-I). The cohort was extended in 2000 (recruiting participants ≥ 55 years) and again in 2006 (recruiting participants ≥ 45 years), forming cohort RS-II and RS-III, respectively. At baseline, and at every three to five year follow-up visit, participants took part in a home interview and a wide set of examinations at the research center [24].
The baseline cohort of the Rotterdam Study consisted of 14,926 participants (RS-I: 7983, RS-II: 3011, RS-III: 3932). From this group, dietary data were available for 9788 participants. Based on previous studies, we determined a reported energy intake of less than 500 kcal/day or more than 5000 kcal/day as unreliable and, therefore, excluded 88 participants [25]. From the remaining 9700 participants with reliable and available dietary data, participants with prevalent parkinsonism or PD (n = 51), without a screening for parkinsonism at baseline (n = 174) or who did not give consent for parkinsonism follow-up (n = 61), were excluded. This resulted in a final eligible cohort of 9414 participants for the current study. The flowchart of study participants can be found in Figure 1.

2.2. Dietary Intake Assessment and Dietary Pattern Construction

Diet was assessed at baseline with semi-quantitative food frequency questionnaires (FFQ), which were validated for the ranking of food and nutrient intake in the Rotterdam Study and other populations [26,27]. The FFQ evaluated the intake of 170 food items in the preceding year for RS-I and RS-II and was upgraded to evaluate 389 food items in the preceding month for RS-III. Energy and nutrient intake were calculated using Dutch food composition tables (NEVO) [28].

2.3. Dietary Patterns

Adherence to the Dutch dietary guidelines was determined through a Dutch diet quality score, which was previously described in more detail [25]. In short, the score is based on adherence to fourteen components of the 2015 Dutch food-based dietary guidelines [29]. A score of zero was given to participants who did not meet the criterion and a score of one for participants who met the criterion, resulting in a maximum sum score of fourteen points. Adherence to the Mediterranean diet was determined with the Mediterranean diet score that was created by Panagiotakos et al. [23]. This sum score is based on eleven main components of the Mediterranean diet. Participants were assigned zero to five points for adherence to each component, resulting in a maximum sum score of 55 points [23].
In addition to pre-defined diet scores, we created dietary patterns derived from actual food intake in the population in order to determine differences in effects as compared to the pre-defined diet scores. For this, food items from the FFQ were combined into 23 food groups based on nutrient composition and everyday use. The food groups are described in Table S1. These food groups were used for the construction of dietary patterns with principal component analysis (PCA) and Varimax rotation for better interpretation. Based on scree plots, we included components with an eigenvalue above 1.5 for further analyses with PD risk. A component is defined by a cluster of food groups, where each food group has a factor loading that indicates the strength of the association between the food group and the component. To characterize dietary patterns from the components, we focused on factor loadings below −0.3 or higher than 0.3. For each participant and for each of the dietary patterns standardized adherence scores (Z-scores) were calculated by summing up the intake per food group and weighing this by the factor loading corresponding to that food group. To determine the effects of the individual food groups on PD risk, all 23 food groups constructed for the PCA were also analyzed separately.

2.4. Parkinson’s Disease Ascertainment

Incident PD during follow-up was ascertained through continuous linkage with clinical records from general practitioners, nursing home physicians, hospitals and pharmaceutical records. In addition, participants could self-report a new diagnosis of PD during the home interview and were screened for parkinsonian symptoms during follow-up visits. All case reports of incident events were evaluated by a panel led by an experienced neurologist, where the final diagnosis was confirmed. This process has previously been described more thoroughly [30]. Incident parkinsonism included: PD, drug-induced parkinsonism, vascular parkinsonism, Lewy body disease, parkinsonism with dementia other than Lewy body disease, multi-system atrophy, progressive supranuclear palsy, corticobasal degeneration or parkinsonism resulting from a tumor. Follow-up continued until the first of January 2018, or until ascertained parkinsonism, death, or loss to follow-up, whichever came first.

2.5. Covariate Assessment

Covariate data were collected from the study round in which the FFQ was completed. Weight and height were measured at the research center, from which body mass index (BMI) was calculated by dividing weight by height squared (kg/m²). Further covariate data were self-reported using questionnaires. Smoking behavior was determined as all forms of tobacco smoking and was classified into never, former or current smoking. Levels of education were low (primary education), lower-middle (lower or intermediate general or lower vocational education), middle (intermediate vocational or higher general education) or high (higher vocational education or university). Genetic ancestry based on genome-wide association studies (GWAS), was estimated using Admixture [31]. Ancestral groups were created including individuals having at least 50% genetic material from that group, and individuals with less than 50% genetic material in any group were classified as ‘Admixed’. The groups of genetic ancestry were European, East Asian, African and Admixed.

2.6. Statistical Analysis

Missing values of genetic ancestry (10.2%), education (0.6%), BMI (0.7%) and smoking behavior (0.5%) were multiply imputed using all covariates, the Dutch diet quality score and the outcome as predictors (5 imputations). We used Cox proportional hazards models to assess the association between the pre-defined diet scores, the data-driven dietary patterns, the food groups and PD risk. Hazard ratios (HR) and 95% confidence intervals (CI) were computed and calculated per standard deviation (SD) increase for the diet scores, the dietary patterns and the food groups. Additionally, we calculated the HRs for the diet scores and dietary patterns per tertile, with the lowest tertile as reference. The basic model was adjusted for sex, age and cohort and the covariate model was additionally adjusted for education, smoking behavior, BMI and energy intake. The proportional hazard assumption was checked with Schoenfeld residuals; there were no clear violations [32].
We performed several sensitivity analyses to examine the robustness of our results. First, we reran the analysis for different episodes of follow-up to compare short- versus long-term effects. Second, we studied the association between the dietary patterns and PD for the cohorts RS-I, RS-II and RS-III separately to assess whether the difference in FFQ or other diet-related cohort differences affected our results. Finally, we evaluated differences in effects between males and females and evaluated whether using all-cause parkinsonism as the outcome changed our results.
Instead of dichotomizing the found effects into ‘significant’ and ‘non-significant’, we will show effect estimates and their corresponding confidence intervals in order to describe the range of estimates compatible with the data. All analyses were performed in R.4.0.3 (R Core Team, Vienna, Austria) [33].

3. Results

3.1. Population Characteristics

The characteristics of the study population are shown in Table 1. The median age of the study population was 62.2 years (interquartile range, IQR (58–70)) and 57.8% were female. Most participants were current or former smokers (67.8%) and had a low or lower-middle level of education (56.6%). The median energy intake was 2020 kcal/day (IQR (1682–2410)), the median Dutch diet quality score was 7.0 out of 14 (IQR (5–8)) and the median Mediterranean diet score was 37 out of 55 (IQR (34–39)). After a mean follow-up of 14.1 years, 129 participants were diagnosed with PD. The mean age at diagnosis of PD was 76.6 years (standard deviation, SD ± 7.4) and 47.3% of the participants with incident PD were female.

3.2. Dietary Patterns and Parkinson’s Disease Risk

Using PCA, we identified three components with an eigenvalue above 1.5 for which the Varimax rotated factor loadings are presented in Table 2. We named the dietary patterns based on the factor loadings of the clustered food groups. The pattern we named ‘Prudent’ was defined by a high intake of vegetables, legumes, white meat, fish and nuts. The ‘Unhealthy’ pattern was defined by high intake of alcoholic beverages, coffee, eggs, red or processed meat, animal-based fats and by low intake of fruits and tea. The ‘Traditional Dutch’ pattern was defined by high intake of cheese, vegetable oils and spreads, potatoes and fries, wholegrain products and sweet snacks. The three patterns explained 8.6%, 7.5% and 7.7% of the total variance, respectively, resulting in a total explained variance of 23.8%.
Associations between the diet scores, data-driven dietary patterns and the risk of PD are shown in Table 3. No associations with the risk of PD were observed for the Dutch diet quality score, the ‘Unhealthy’ pattern and the ‘Traditional Dutch’ pattern, neither in the continuous nor the categorical analysis (adjusted HRs per SD: Dutch diet quality score 0.93 (95% CI 0.77–1.12), ‘Unhealthy’ pattern 1.05 (95% CI 0.85–1.29), ‘Traditional Dutch’ pattern 0.90 (95% CI 0.69–1.17)). Estimates of the unadjusted and adjusted models for the Mediterranean diet score (adjusted HR per SD 0.89 (95% CI 0.74–1.07)) and the ‘Prudent’ pattern (adjusted HR per SD 0.81 (95% CI 0.61–1.08)) suggested a lower hazard of PD with higher scores, although confidence intervals were wide.
Figure 2 shows the associations between the individual food groups from the PCA and the risk of PD. A higher intake of yoghurt and fermented milk (adjusted HR per SD 1.13 (95% CI 0.98–1.31)), animal-based fats (adjusted HR per SD 1.19 (95% CI 1.00–1.43)) and sweet snacks (adjusted HR per SD 1.21 (95% CI 1.03–1.41)) was associated with a higher risk of PD. Although confidence intervals were wide, we observed a lower risk of PD with a higher intake of fruits (adjusted HR per SD 0.82 (95% CI 0.63–1.07)), vegetables (adjusted HR per SD 0.84 (95% CI 0.65–1.09)), cheese (adjusted HR per SD 0.79 (95% CI 0.62–1.01)) and savory snacks (adjusted HR per SD 0.73 (95% CI 0.45–1.19)).

3.3. Sensitivity Analyses

We observed no clear difference in associations between the dietary patterns and the risk of PD with different follow-up times or for the three cohorts separately (Tables S2 and S3). When stratifying for sex, stronger effect estimates were found in the association of dietary patterns with the risk of PD for males, although confidence intervals overlapped between the two strata (Table S4). Finally, changing the outcome to all-cause parkinsonism increased the number of total incident events in the eligible cohort to 254. In this analysis, the association between the ‘Prudent’ dietary pattern and the Mediterranean diet score with risk of all-cause parkinsonism became more pronounced (Table S5: ‘Prudent’ pattern adjusted HR per SD 0.76 (95% CI 0.61–0.95), Mediterranean diet adjusted HR per SD 0.86 (95% CI 0.76–0.98)).

4. Discussion

In this population-based cohort study, we corroborated previous findings of a relation between the Mediterranean diet and the risk of PD. However, we did not find clear effects of adherence to the Dutch dietary guidelines on the risk of PD. In addition, our study provides new evidence for the possible protective and detrimental effects of separate food groups on the risk of PD.
A previous Finnish cohort study found no significant association between the Alternate Healthy Eating Index (AHEI)—a diet score reflecting US guidelines—and risk of PD, while an American cohort study found a lower risk of PD with a higher AHEI score [20,22,34]. However, in the same cohort studies, and also more recently in a Swedish cohort study, a higher Mediterranean diet score was consistently associated with a lower risk of PD [20,21,22,23,35]. These previous studies are in line with our findings of a possible association between the Mediterranean diet and the ‘Prudent’ dietary pattern and the risk of PD, without clear effects of the Dutch diet quality score.
Our more detailed food group analyses revealed that higher intake of yoghurt and fermented milk, animal-based fats and sweet snacks was associated with a higher risk of PD. The association between intake of sweet snacks and PD is novel; it is yet unclear what underlies this association. Previous research on impulsive behavior suggested higher impulsivity in PD patients, and possibly, the impulsive intake of sweet snacks could explain the found association [36]. This might also imply that snacking is a consequence of prodromal PD rather than a cause, which needs further investigation. Based on previous research, we also hypothesized that a frequent intake of dairy products, mostly milk intake, would be associated with an increased PD risk, but we found inconsistent results concerning this association [12,14]. We found no association between higher milk and cream intake and PD, but we did find a possible higher risk of PD for yoghurt and fermented milk intake and animal-based fat intake, which mostly consists of dairy butter. The positive association between dietary animal fat intake and PD is in line with results from a previous case–control study by Logroscino et al. [15]. However, the association remains inconsistent in the literature, as Chen et al. found no significant association between animal fat intake and PD in the Nurses’ Health Study and the Health Professionals Follow-up Study [16]. The possible associations we found between cheese and savory snacks and the risk of PD are remarkable and have not been reported by previous studies, and thus, require further study.
Our study has several strengths. First, the population-based study design of the Rotterdam Study with a continuous follow-up of participants through clinical records and follow-up visits reinforces the validity of the incident PD events in this study. Second, by using pre-defined diet scores, data-driven dietary patterns and separate food groups, we were able to secure a complete image of diet quality. In addition, with both pre-defined diet scores and data-driven dietary patterns derived from the PCA, we managed to explore and compare PD risk based on both desired intake reflecting guidelines and the actual diet patterns from the data. Using this two-faced approach, we observed similar results for the Mediterranean diet score and the comparable ‘Prudent’ data-driven dietary pattern, which strengthens our conclusions. We found only one study that has used the combination of data-driven patterns and pre-defined diet scores in association with PD [22].
Nonetheless, there are also limitations in our study. Firstly, we used only baseline diet information. As dietary habits are subject to changes over time, diet could have altered during the follow-up period, which was not taken into account in the current study [37]. Secondly, the low number of incident PD events in our study gave rise to wide confidence intervals around our estimates, which indeed became remarkably smaller when more events were included by using all-cause parkinsonism as the outcome. Thirdly, our study sample included solely Dutch individuals living in the Rotterdam district Ommoord; the dietary patterns identified in this study and their associations with PD are thus not necessarily generalizable to other study populations. Lastly, two different FFQs were used at baseline to assess dietary intake, which could have resulted in differences in the determinant between the study cohorts. Nevertheless, both FFQs included the necessary items to calculate data-driven dietary patterns and the specific food groups. Moreover, repeating the analyses for each cohort separately revealed similar results, implying that this difference in questionnaires did not impact the findings much.
Taking into account both our findings and the findings of previous studies, the Mediterranean diet seems to be a promising focus for future research aiming to study the preventive potential of lifestyle interventions for PD [20,22,35]. Nevertheless, the mechanisms underlying the association between the Mediterranean diet and PD are still to be resolved. Potentially, its anti-inflammatory character and ability to reduce oxidative stress could contribute to the beneficial effect on PD risk [7,38]. These mechanisms, as well as potential beneficial effects of other healthy dietary patterns, remain subject to further investigation.

5. Conclusions

In summary, our results corroborate previous findings of an association between the Mediterranean diet and the risk of PD. In addition, our results suggest a possible link between dairy products, fruits, vegetables and sweet and savory snacks and the risk of PD, which warrants replication by future studies.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/nu13113970/s1, Table S1. Description of the food groups; Table S2: Time-stratified analyses for association between dietary patterns and Parkinson’s disease; Table S3: Analyses per cohort for association between dietary patterns and Parkinson’s disease; Table S4: Sex-stratified analyses for association between dietary patterns and Parkinson’s disease; Table S5: Associations between dietary patterns and risk of parkinsonism (n = 9414).

Author Contributions

Conceptualization, A.J.S., L.J.D. and T.V.; Methodology, A.J.S. and L.J.D.; Software, A.J.S.; Validation, L.J.D.; Formal Analysis, A.J.S. and L.J.D.; Investigation, A.J.S. and L.J.D.; Resources, T.V. and M.K.I.; Data Curation, T.V., M.K.I. and L.J.D.; Writing—Original Draft Preparation, A.J.S.; Writing—Review and Editing, L.J.D., M.K.I. and T.V.; Visualization, A.J.S.; Supervision, L.J.D., T.V. and M.K.I.; Project Administration, A.J.S. and L.J.D.; Funding Acquisition, T.V. and M.K.I. All authors have read and agreed to the published version of the manuscript.

Funding

The Rotterdam Study is funded by Erasmus Medical Center and Erasmus University, Rotterdam, Netherlands Organization for the Health Research and Development (ZonMw), the Research Institute for Diseases in the Elderly (RIDE), the Ministry of Education, Culture and Science, the Ministry for Health, Welfare and Sports, the European Commission (DG XII), and the Municipality of Rotterdam. This study received further support from Stiching ParkinsonFonds.

Institutional Review Board Statement

The Rotterdam Study has been approved by the Medical Ethics Committee of the Erasmus MC (registration number MEC 02.1015) and by the Dutch Ministry of Health, Welfare and Sport (Population Screening Act WBO, license number 1071272-159521-PG). The Rotterdam Study Personal Registration Data collection is filed with the Erasmus MC Data Protection Officer under registration number EMC1712001. The Rotterdam Study has been entered into the Netherlands National Trial Register (NTR; www.trialregister.nl) and into the WHO International Clinical Trials Registry Platform (ICTRP; www.who.int/ictrp/network/primary/en/) under shared catalogue number NTR6831.

Informed Consent Statement

All participants provided written informed consent to participate in the study and to have their information obtained from treating physicians.

Data Availability Statement

Data can be obtained upon request. Requests should be directed towards the management team of the Rotterdam Study ([email protected]), which has a protocol for approving data requests. Because of restrictions based on privacy regulations and informed consent of the participants, data cannot be made freely available in a public repository.

Acknowledgments

We are thankful for the contribution of the participants of the Rotterdam study, the staff of the Rotterdam Study and the general practitioners and pharmacists.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dorsey, E.R.; Elbaz, A.; Nichols, E.; Abd-Allah, F.; Abdelalim, A.; Adsuar, J.C.; Ansha, M.G.; Brayne, C.; Choi, J.-Y.J.; Collado-Mateo, D.; et al. Global, regional, and national burden of Parkinson’s disease, 1990–2016: A systematic analysis for the Global Burden of Disease Study 2016. Lancet Neurol. 2018, 17, 939–953. [Google Scholar] [CrossRef] [Green Version]
  2. De Lau, L.M.; Breteler, M.M. Epidemiology of Parkinson’s Disease. Lancet Neurol. 2006, 5, 525–535. [Google Scholar] [CrossRef]
  3. Warner, T.T.; Schapira, A.H. Genetic and environmental factors in the cause of Parkinson’s Disease. Ann. Neurol. 2003, 53, 16–23. [Google Scholar] [CrossRef]
  4. Ascherio, A.; Schwarzschild, M.A. The epidemiology of Parkinson’s disease: Risk factors and prevention. Lancet Neurol. 2016, 15, 1257–1272. [Google Scholar] [CrossRef]
  5. Meng, L.; Yuan, X.; Cao, X.; Zhang, Z. The gut-brain axis in the pathogenesis of Parkinson’s disease. Brain Sci. Adv. 2019, 5, 73–81. [Google Scholar] [CrossRef] [Green Version]
  6. Wu, G.D.; Chen, J.; Hoffmann, C.; Bittinger, K. Linking long-term dietary patterns with gut microbiotal enterotypes. Science 2011, 334, 105–108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Jackson, A.; Forsyth, C.B.; Shaikh, M.; Voigt, R.M.; Engen, P.A.; Ramirez, V.; Keshavarzian, A. Diet in Parkinson’s Disease: Critical role for the Microbiome. Front. Neurol. 2019, 10, 1245. [Google Scholar] [CrossRef] [PubMed]
  8. Horsager, J.; Andersen, K.B.; Knudsen, K.; Skjærbæk, C.; Fedorova, T.D.; Okkels, N.; Schaeffer, E.; Bonkat, S.K.; Geday, J.; Otto, M.; et al. Brain-first versus body-first Parkinson’s disease: A multimodal imaging case-control study. Brain 2020, 143, 3077–3088. [Google Scholar] [CrossRef] [PubMed]
  9. Palacios, N.; Gao, X.; McCullough, M.L.; Schwarzschild, M.A.; Shah, R.; Gapstur, S.; Ascherio, A. Caffeine and risk of Parkinson’s disease in a large cohort of men and women. Mov. Disord. 2012, 27, 1276–1282. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Qi, H.; Li, S. Dose-response meta-analysis on coffee, tea and caffeine consumption with risk of Parkinson’s disease. Geriatr. Gerontol. Int. 2014, 14, 430–439. [Google Scholar] [CrossRef]
  11. Chen, H.; O’Reilly, E.; McCullough, M.L.; Rodriquez, C.; Schwarzschild, M.A.; Calle, E.E.; Thun, M.J.; Ascherio, A. Consumption of dairy products and risk of Parkinson’s disease. Am. J. Epidemiol. 2007, 165, 998–1006. [Google Scholar] [CrossRef]
  12. Jiang, W.; Ju, C.; Jiang, H.; Zhang, D. Dairy foods intake and risk of Parkinson’s disease: A dose-response meta-analysis of prospective cohort studies. Eur. J. Epidemiol. 2014, 29, 613–619. [Google Scholar] [CrossRef]
  13. Hughes, K.C.; Gao, X.; Kim, I.Y.; Wang, M.; Weisskopf, M.G.; Schwarzschild, M.A.; Ascherio, A. Intake of dairy foods and risk of Parkinson Disease. Neurology 2017, 89, 46–52. [Google Scholar] [CrossRef] [PubMed]
  14. Olsson, E.; Byberg, L.; Höijer, J.; Kilander, L.; Larsson, S.C. Milk and Fermented Milk Intake and Parkinson’s Disease: Cohort Study. Nutrients 2020, 12, 2763. [Google Scholar] [CrossRef] [PubMed]
  15. Logroscino, G.; Marder, K.; Cote, L.; Tang, M.X.; Shea, S.; Mayeux, R. Dietary lipids and antioxidants in Parkinson’s disease: A population-based, case-control study. Ann. Neurol. 1996, 39, 89–94. [Google Scholar] [CrossRef] [PubMed]
  16. Chen, H.; Zhang, S.M.; Hernán, M.A.; Willett, W.C.; Ascherio, A. Dietary intakes of fat and risk of Parkinson’s Disease. Am. J. Epidemiol. 2003, 157, 1007–1014. [Google Scholar] [CrossRef] [Green Version]
  17. Palacios, N.; Gao, X.; O’Reilly, E.; Schwarzschild, M.A. Alcohol an Risk of Parkinson Disease in a Large Prospective Cohort of Men and Women. Mov. Disord. 2012, 27, 980–987. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Zhang, D.; Jiang, H.; Xie, J. Alcohol intake and risk of Parkinson’s disease: A meta-analysis of observational studies. Mov. Disord. 2014, 29, 819–822. [Google Scholar] [CrossRef]
  19. Jacques, P.F.; Tucker, K.L. Are dietary patterns useful for understanding the role of diet in chronic disease? Am. J. Clin. Nutr. 2001, 73, 1–2. [Google Scholar] [CrossRef] [Green Version]
  20. Sääksjärvi, K.; Knekt, P.; Lundqvist, A.; Männistö, S.; Heliövaara, H.; Rissanen, H.; Järvinen, R. A cohort study on diet and the risk of Parkinson’s disease: The role of food groups and diet quality. Br. J. Nutr. 2013, 109, 329–337. [Google Scholar] [CrossRef] [Green Version]
  21. Maraki, M.I.; Yannakoulia, M.; Stamelou, M.; Stefanis, L.; Xiromerisiou, G.; Kosmidis, M.H.; Dardiotis, E.; Hadjigeorgiou, G.M.; Sakka, P.; Anastasiou, C.A.; et al. Mediterranean diet adherence is related to reduced probability of prodromal Parkinson’s disease. Mov. Disord. 2019, 34, 48–57. [Google Scholar] [CrossRef]
  22. Gao, X.; Chen, H.; Fung, T.T.; Logroscino, G.; Schwarzschild, M.A.; Hu, F.B.; Ascherio, A. Prospective study of dietary pattern and risk of Parkinson Disease. Am. J. Clin. Nutr. 2007, 86, 1486–1494. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Panagiotakos, D.B.; Pitsavos, C.; Stefanadis, C. Dietary patterns: A Mediteranean diet score and its relation to clinical and biological markers of cardiovascular disease risk. Nutr. Metab. Cardiovasc. Dis. 2006, 16, 559–568. [Google Scholar] [CrossRef] [PubMed]
  24. Ikram, M.A.; Brusselle, G.; Ghanbari, M.; Goedgebure, A.; Ikram, M.K.; Kavousi, M.; Kieboom, B.C.T.; Klaver, C.C.W.; De Knegt, R.J.; Luik, A.I.; et al. Objectives, design and main findings until 2020 from the Rotterdam Study. Eur. J. Epidemiol. 2020, 35, 483–517. [Google Scholar] [CrossRef]
  25. Voortman, T.; Kiefte-de Jong, J.C.; Ikram, M.A.; Stricker, B.H.; Van Rooij, F.J.A.; Lahousse, L.; Tiemeier, H.; Brusselle, G.; Franco, O.H.; Schoufour, J.D. Adherence to the 2015 Dutch dietary guidelines and risk of non-communicable diseases and mortality in the Rotterdam Study. Eur. J. Epidemiol. 2017, 32, 993–1005. [Google Scholar] [CrossRef] [PubMed]
  26. Klipstein-Grobusch, K.; Den Breeijen, J.H.; Goldbohm, R.A.; Geleijnse, J.M.; Hofman, A.; Grobbee, D.E.; Witteman, J.C. Dietary assessment in the elderly: Validation of a semiquantitative food frequency questionnaire. Eur. J. Clin. Nutr. 1998, 52, 588–596. [Google Scholar] [CrossRef] [Green Version]
  27. Goldbohm, R.A.; Van den Brandt, P.A.; Brants, H.A.; Van’t Veer, P.; Al, M.; Sturmans, F.; Hermus, R.J. Validation of a dietary questionnaire used in a large-scale prospective cohort study on diet and cancer. Eur. J. Clin. Nutr. 1994, 48, 253–265. [Google Scholar]
  28. RIVM. NEVO-Table. Updated 2019. Available online: https://nevo-online.rivm.nl/ (accessed on 5 November 2021).
  29. Kromhout, D.; Spaaij, C.J.K.; De Goede, J.; Weggemans, R.M. The 2015 Dutch food-based dietary guidelines. Eur. J. Clin. Nutr. 2016, 70, 869–878. [Google Scholar] [CrossRef]
  30. Darweesh, S.K.L.; Koudstaal, P.J.; Stricker, B.H.; Hofman, A.; Ikram, M.A. Incidence of Parkinson Disease in the General Population: The Rotterdam Study. Am. J. Epidemiol. 2016, 183, 1018–1026. [Google Scholar] [CrossRef] [Green Version]
  31. Alexander, D.H.; Novembre, J.; Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009, 19, 1655–1664. [Google Scholar] [CrossRef] [Green Version]
  32. Grambsch, P.M.; Therneau, T.M. Proportional hazards tests and diagnostics based on weighted residuals. Biometrika 1994, 81, 515–526. [Google Scholar] [CrossRef]
  33. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  34. McCullough, M.L.; Feskanich, D.; Stampfer, M.J.; Giovannucci, E.L.; Rimm, E.B.; Hu, F.B.; Spiegelman, D.; Hunter, D.J.; Colditz, G.A.; Willett, W.C. Diet quality and major chronic disease risk in men and women: Moving toward improved dietary guidance. Am. J. Clin. Nutr. 2002, 76, 1261–1271. [Google Scholar] [CrossRef] [Green Version]
  35. Yin, W.; Löf, M.; Pedersen, N.L.; Sandin, S.; Fang, F. Mediterranean Dietary Pattern at Middle Age and Risk of Parkinson’s Disease: A Swedish Cohort Study. Mov. Disord. 2021, 36, 255–260. [Google Scholar] [CrossRef]
  36. Aumann, M.A.; Stark, A.J.; Hughes, S.B.; Lin, Y.C.; Kang, H.; Bradly, E.; Zald, D.H.; Claassen, D.O. Self-reported rates of impulsivity in Parkinson’s Disease. Ann. Clin. Transl. Neurol. 2020, 7, 437–448. [Google Scholar] [CrossRef] [PubMed]
  37. Mishra, G.D.; McNaughton, S.A.; Bramwell, G.D.; Wadsworth, M.E.J. Longitudinal changes in dietary patterns during adult life. Br. J. Nutr. 2006, 96, 735–744. [Google Scholar] [PubMed]
  38. Alcalay, R.; Gu, Y.; Meija-Santana, H.; Cote, L.; Marder, K.S.; Scarmeas, N. The Association between Mediterranean Diet Adherence and Parkinson’s Disease. Mov. Disord. 2012, 27, 771–774. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Flowchart for selection of study participants.
Figure 1. Flowchart for selection of study participants.
Nutrients 13 03970 g001
Figure 2. Associations between the food groups and Parkinson’s disease. The hazard ratios (HR), shown as diamonds, and the 95% confidence interval (CI), shown as lines, were obtained using Cox proportional hazard models. Hazard ratios are shown per standard deviation increase in intake for the 23 food groups, which were also used in the principal component analysis.
Figure 2. Associations between the food groups and Parkinson’s disease. The hazard ratios (HR), shown as diamonds, and the 95% confidence interval (CI), shown as lines, were obtained using Cox proportional hazard models. Hazard ratios are shown per standard deviation increase in intake for the 23 food groups, which were also used in the principal component analysis.
Nutrients 13 03970 g002
Table 1. Baseline characteristics.
Table 1. Baseline characteristics.
n = 9414
Sex, female5439 (57.8%)
Age, years62.2 (58.0–70.0)
Genetic ancestry, European a9251 (98.3%)
Education a
Low1456 (15.5%)
Lower-middle3872 (41.1%)
Middle2631 (27.9%)
High1455 (15.5%)
Smoking behavior a
Never3023 (32.1%)
Former4178 (44.4%)
Current2213 (23.5%)
Alcohol intake, drinks/day0.5 (0.0–1.8)
BMI, kg/m² a26.8 (±4.0)
Energy intake, kcal/day2020 (1,682–2,410)
Dutch diet quality score b7.0 (5.0–8.0)
Mediterranean diet score c37.0 (34.0–39.0)
Values are median (interquartile range, IQR), mean (standard deviation, ±SD) or n (%). a Genetic ancestry, education, smoking behavior and body mass index (BMI) are based on imputed values. The number of missing values were 969 for genetic ancestry, 55 for education, 45 for smoking behavior and 70 for BMI. b Theoretical range 0–14, where a higher score is better adherence to the Dutch dietary guidelines. c Theoretical range 0–55, where a higher score is better adherence to the Mediterranean diet.
Table 2. Factor loadings for each group per dietary pattern.
Table 2. Factor loadings for each group per dietary pattern.
Food GroupsLoadings Per Dietary Pattern
PrudentUnhealthyTraditional Dutch
Vegetables0.56−0.230.15
Fruits0.41−0.44-
Legumes0.46--
Milk and cream−0.11-0.18
Yoghurt and fermented milk-−0.38-
Eggs0.280.320.12
Cheese--0.39
Unprocessed white meat0.41--
Processed and red meat-0.430.44
Fish and seafood0.51-−0.13
Vegetable oils and spreads-−0.160.63
Animal-based fats-0.370.12
Wholegrain products0.14−0.320.54
Refined grain products0.350.25-
Potatoes and fries−0.17-0.62
Sugar-containing beverages0.270.24-
Tea-−0.47-
Coffee-0.340.34
Alcoholic beverages0.250.46-
Savory snacks0.370.21-
Sweet snacks--0.32
Nuts and seeds0.45--
Soy products and meat replacers0.32−0.18-
% of explained variance8.67.57.7
Dietary patterns and factor loadings were obtained from the principal component analysis. The three components with an eigenvalue >1.5 are presented in the table as the ‘Prudent’, ‘Unhealthy’ and ‘Traditional Dutch’ dietary pattern, which were named based on the factor loadings belonging to the clustered food groups in each component. The presented loadings were rotated with Varimax rotation for better interpretation. All factor loadings that are higher than 0.3 and lower than −0.3 are printed in bold. Details on the food groups are presented in Table S1.
Table 3. Associations between dietary patterns and Parkinson’s disease (n=9414).
Table 3. Associations between dietary patterns and Parkinson’s disease (n=9414).
Dietary PatternNo. of Incident Parkinson’s Disease Basic Model
HR (95% CI)
Covariate Model
HR (95% CI)
Dutch diet quality score per SD1290.97 (0.81–1.17)0.93 (0.77–1.12)
Tertiles
Low (reference)3711
Medium481.02 (0.66–1.59)0.97 (0.62–1.52)
High441.07 (0.68–1.69)0.98 (0.62–1.55)
Mediterranean diet score per SD1290.93 (0.78–1.12)0.89 (0.74–1.07)
Tertiles
Low3811
Medium571.39 (0.92–2.11)1.33 (0.87–2.02)
High340.88 (0.55–1.40)0.80 (0.50–1.29)
Prudent pattern per SD1290.87 (0.67–1.14)0.81 (0.61–1.08)
Tertiles
Low (reference)5311
Medium470.97 (0.65–1.46)0.94 (0.63–1.42)
High290.95 (0.57–1.56)0.89 (0.52–1.50)
Unhealthy pattern per SD1290.98 (0.80–1.21)1.05 (0.85–1.29)
Tertiles
Low (reference)4411
Medium481.07 (0.70–1.62)1.15 (0.75–1.74)
High370.96 (0.60–1.52)1.08 (0.67–1.74)
Traditional Dutch pattern per SD1290.99 (0.82–1.20)0.90 (0.69–1.17)
Tertiles
Low (reference)2811
Medium511.50 (0.94–2.41)1.44 (0.88–2.35)
High501.27 (0.77–2.10)1.16 (0.64–2.10)
The hazard ratios (HR) and 95% confidence interval (CI), obtained using Cox proportional hazard models, are shown per standard deviation (SD) increase for the dietary patterns and scores. The basic model was adjusted for sex, age at baseline and Rotterdam Study cohort. The covariate model was adjusted for all items in the basic model and additionally for body mass index (BMI), education, smoking behavior and energy intake.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Strikwerda, A.J.; Dommershuijsen, L.J.; Ikram, M.K.; Voortman, T. Diet Quality and Risk of Parkinson’s Disease: The Rotterdam Study. Nutrients 2021, 13, 3970. https://doi.org/10.3390/nu13113970

AMA Style

Strikwerda AJ, Dommershuijsen LJ, Ikram MK, Voortman T. Diet Quality and Risk of Parkinson’s Disease: The Rotterdam Study. Nutrients. 2021; 13(11):3970. https://doi.org/10.3390/nu13113970

Chicago/Turabian Style

Strikwerda, Anne J., Lisanne J. Dommershuijsen, M. Kamran Ikram, and Trudy Voortman. 2021. "Diet Quality and Risk of Parkinson’s Disease: The Rotterdam Study" Nutrients 13, no. 11: 3970. https://doi.org/10.3390/nu13113970

APA Style

Strikwerda, A. J., Dommershuijsen, L. J., Ikram, M. K., & Voortman, T. (2021). Diet Quality and Risk of Parkinson’s Disease: The Rotterdam Study. Nutrients, 13(11), 3970. https://doi.org/10.3390/nu13113970

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop