Genetic Variability of the Vitamin D Receptor Affects Susceptibility to Parkinson’s Disease and Dopaminergic Treatment Adverse Events

Vitamin D is a lipid-soluble molecule and an important transcriptional regulator in many tissues and organs, including the brain. Its role has been demonstrated also in Parkinson’s disease (PD) pathogenesis. Vitamin D receptor (VDR) is responsible for the initiation of vitamin D signaling cascade. The aim of this study was to assess the associations of VDR genetic variability with PD risk and different PD-related phenotypes. We genotyped 231 well characterized PD patients and 161 healthy blood donors for six VDR single nucleotide polymorphisms, namely rs739837, rs4516035, rs11568820, rs731236, rs2228570, and rs1544410. We observed that VDR rs2228570 is associated with PD risk (p < 0.001). Additionally, we observed associations of specific VDR genotypes with adverse events of dopaminergic treatment. VDR rs1544410 (GG vs. GA + AA: p = 0.005; GG vs. GA: p = 0.009) was associated with the occurrence of visual hallucinations and VDR rs739837 (TT vs. GG: p = 0.036), rs731236 (TT vs. TC + CC: p = 0.011; TT vs. TC: p = 0.028; TT vs. CC: p = 0.035), and rs1544410 (GG vs. GA: p = 0.014) with the occurrence of orthostatic hypotension. We believe that the reported study may support personalized approach to PD treatment, especially in terms of monitoring vitamin D level and vitamin D supplementation in patients with high risk VDR genotypes.

Introduction

Vitamin D is a lipid-soluble secosteroid obtained through dermal synthesis by exposure to sunlight or through dietary intake (Fullard and Duda, 2020). The active compound acts as a hormone and binds to the vitamin D receptor (VDR), which regulates expression of more than 200 different genes (Samuel and Sitrin, 2008; Fullard and Duda, 2020; Lv et al., 2020). VDR is present in several brain areas, such as hippocampus, hypothalamus, thalamus, cortex, and substantia nigra (Moretti et al., 2018). Vitamin D thus has an important role in the central nervous system as it is involved in neurotransmission, neuroprotection, and neuroplasticity (DeLuca et al., 2013). All this prompted many studies to investigate the relation between vitamin D homeostasis and neurological disorders. Several disorders with diverse etiologies have thus been associated with vitamin D levels already: multiple sclerosis, Alzheimer’s disease, amyotrophic lateral sclerosis, and Parkinson’s disease (PD) (DeLuca et al., 2013).

Parkinson’s disease is the second most prevalent neurodegenerative brain disorder and is presenting with a variety of motor and non-motor signs and symptoms. Its hallmarks are aggregation of the misfolded α-synuclein in the intracellular inclusions called Lewy bodies and degeneration of dopaminergic neurons in the nigrostriatal pathway (Kalia and Lang, 2015; Poewe et al., 2017; Tolosa et al., 2021). In later stages of PD, motor and non-motor complications due to prolonged dopaminergic treatment develop and they severely impact patients’ and caregivers’ health related quality of life (Zhao et al., 2021). Vitamin D deficiency is common in PD patients. However, it is still unknown how is it related to PD, whether it is a cause or a consequence of PD pathology (Luo et al., 2018). Nevertheless, it was shown in rodent models of PD that VDR affects expression of tyrosine hydroxylase as the rate-limiting step in dopamine synthesis, and also the expression of glial cell line-derived neurotrophic factor (GDNF), which facilitates neuronal regrowth and protects dopaminergic nerve terminals (Lv et al., 2020).

Studies have already shown associations between vitamin D levels and PD risk. A 29-year long prospective study conducted in Finland showed that people with increased vitamin D serum concentrations had lower risk for PD development (Knekt et al., 2010). Several studies also showed that higher exposure to sunlight, which increases vitamin D production in the skin, protects against PD development (Kenborg et al., 2011; Zhu et al., 2014; Wang et al., 2016; Kravietz et al., 2017). Additionally, vitamin D levels were also shown to be correlated with severity of motor symptoms in PD, scored by the MDS-Unified Parkinson’s Disease Rating Scale (MDS-UPDRS) and Hoehn and Yahr (H&Y) scales. Both scales were inversely correlated with serum and plasma vitamin D levels (Chitsaz et al., 2013; Ding et al., 2013; Fullard and Duda, 2020). There are also some reports available stating that vitamin D levels are associated with non-motor symptoms, although these studies are fairly inconsistent. Higher plasma vitamin D levels were associated with better verbal fluency and verbal memory in PD without dementia (Peterson et al., 2013). Several studies have shown an association between vitamin D levels and depression in the general elderly population (Wilkins et al., 2006; Peterson et al., 2013), similarly as orthostatic hypotension (Duval et al., 2015). Finally, vitamin D levels also showed association with olfactory dysfunction (Kim et al., 2018).

Genetic variability in VDR has already been investigated in several different settings. More specifically, VDR rs2228570 was associated with the risk of PD development (Han et al., 2012; Török et al., 2013; Li et al., 2015; Niu et al., 2015; Wang et al., 2019; Gao et al., 2020; Hu et al., 2020; Agliardi et al., 2021). This single nucleotide polymorphism (SNP) was also associated with the disease progression (Suzuki et al., 2012) and with cognitive decline in PD patients (Gatto et al., 2015). A few other studies demonstrated the association of VDR rs1544410 with PD risk (Kim et al., 2005; Li et al., 2015).

However, none of the studies investigated the association of VDR genetic variability with the occurrence of adverse events (AEs) of treatment with levodopa and dopamine agonists (ADs). Additionally, no studies have assessed whether VDR genetic variability influences the daily dose requirements of dopaminergic treatment necessary to adequately control PD symptoms in an individual patient. The aim of our study was to address these two scientific questions and additionally also to assess the association between VDR SNPs and PD risk in a Slovenian cohort of PD patients.

Materials and Methods

Participants and Clinical Data

In total, 231 consecutive PD patients were included in this study. Patients were recruited at the Department of Neurology, University Medical Centre Ljubljana, Slovenia, between October 2016 and April 2018 (Redenšek et al., 2019). Inclusion criteria are listed elsewhere already (Redenšek et al., 2019). The criteria were the following: (1) diagnosis of PD as defined in the United Kingdom Parkinson’s Disease Society Brain Bank criteria (Goetz et al., 2008) by an experienced movement disorders specialist, (2) available clinical data, (3) at least 3 months of levodopa and/or DAs treatment duration, (4) ongoing dopaminergic treatment with levodopa and/or DAs. Patients with atypical and secondary forms of parkinsonisms were not included in the study.

Patients and their caregivers underwent a structured interview to obtain demographic and clinical data of interest. Additional clinical data was obtained from medical records. Data on gender, side of disease initiation, presence of tremor, treatment with DAs, age at diagnosis, disease duration, and levodopa treatment duration was collected. Levodopa equivalent dose (LED) at the inclusion in the study was used as a measure of the disease severity level. We also collected data on the occurrence of AEs during the course of dopaminergic treatment. The dopaminergic treatment related AEs of interest were the following: motor fluctuations (MFs), dyskinesia, excessive daytime sleepiness and sleep attacks, visual hallucinations (VHs), nausea/vomiting, orthostatic hypotension, peripheral edema, and impulse control disorders (ICDs).

A control group of 161 unrelated Slovenian blood donors without neurological and oncological disorders were included in the study for the purpose of PD susceptibility evaluation. We collected data on their age and gender.

The study protocol was approved by the Slovenian Ethics Committee for Research in Medicine (KME 42/05/16). All subjects gave written informed consent in accordance with the Declaration of Helsinki.

Single Nucleotide Polymorphism Selection

In total, six functional SNPs from regulatory regions of the VDR gene were selected for the analysis, namely rs739837, rs4516035, rs11568820, rs731236, rs2228570, and rs1544410. All of the SNPs had a function already reported in the literature or predicted by the SNP function prediction tool (Xu and Taylor, 2009). We selected SNPs based on their previous association with PD and related phenotypes (Kim et al., 2005; Han et al., 2012; Suzuki et al., 2012; Török et al., 2013; Gatto et al., 2015; Li et al., 2015; Niu et al., 2015; Wang et al., 2019; Gao et al., 2020; Hu et al., 2020; Agliardi et al., 2021) and based on their location in the gene and minor allele frequency (MAF). Only very common SNPs with the MAF above 10% were selected. SNPs from coding region, 5′UTR, and 3′UTR were prioritized. Two selected SNPs are located in the coding region (rs731236 and rs2228570), two SNPs in the 5′UTR (rs4516035 and rs11568820), one in the 3′UTR (rs739837), and one in the intronic region (rs1544410). The intronic SNP was included into the investigation due to the previously reported association with PD risk (Kim et al., 2005; Li et al., 2015). The characteristics of the investigated SNPs are presented in the Table 1.

TABLE 1

Table 1. Genotype frequencies and characteristics of the genetic variants studied in our PD cohort.

DNA Isolation and Genotyping

Peripheral blood samples were obtained for DNA extraction and genomic DNA was isolated using the FlexiGene DNA Kit (Qiagen, Hilden, Germany) in the course of our previous study (Redenšek et al., 2019). Six SNPs (VDR rs739837, rs4516035, rs11568820, rs731236, rs2228570, and rs1544410) were genotyped with KASPar assays (KBiosciences, Herts, United Kingdom and LGC Genomics, United Kingdom) according to manufacturer’s instructions. In total, 10% of samples were genotyped in duplicate as quality control and all the results were concordant.

Statistical Analysis

Median and 25th to 75th percentile range were used to describe central tendency and variability of continuous variables, while frequencies were used to describe the distribution of categorical variables. The agreement of genotype frequencies with Hardy–Weinberg equilibrium was assessed by Chi-squared test in the control cohort.

Logistic regression was used to calculate odds ratios (ORs) and 95% confidence intervals (95% CIs) to examine the associations of selected SNPs and clinical data with the risk for AEs. Logistic regression was used for the risk analysis as well. Dominant, recessive, and additive genetic models were used according to the genotype frequencies. Non-parametric Mann−Whitney U test and Kruskal−Wallis test were used to assess the effects of genotypes on the LED. ANCOVA was used to adjust for the clinical variables.

All statistical tests were two-sided. Bonferroni correction was used to account for multiple comparisons to prevent false positive results. p-values up to 0.0083 (0.05/6) were considered statistically significant, while p-values between 0.0083 and 0.0500 were considered nominally significant. For an allelic variant with an average MAF of 0.39 and with a 32% prevalence of an AE, this study had 80% or more power to detect OR of 0.40 or less and OR of 2.23 or more. More detailed power calculations are presented in the Supplementary Table 1. Power calculations were conducted by the PS Power and sample size calculations, version 3.0. All statistical analyses were carried out by IBM SPSS Statistics, version 21.0 (IBM Corporation, Armonk, NY, United States).

Results

Patients’ Characteristics

A total of 132 male and 99 female PD patients were included in the study. Patients’ median age at enrollment was 72.5 years (65.7–78.0). Their median age at diagnosis was 62.1 years (54.8–71.7) and their median disease duration was 7.6 years (3.8–13.6). The median dopaminergic treatment duration was 7.3 years (3.6–13.5) and the median LED at enrollment was 975 mg/day (600–1,363.5).

The control cohort included 125 men and 36 women with a median age of 55 years (52–58.5).

The main patient and control cohort characteristics are presented in the Table 2. Other characteristics of the patient cohort and the data on the frequencies of AEs are summarized in the Supplementary Table 2, which is adjusted according to Redenšek et al. (2019).

TABLE 2

Table 2. Patient and control cohort characteristics.

Single Nucleotide Polymorphism Genotyping Analysis

In total, six SNPs were analyzed. Data on genotype frequencies of these variants along with their locations, MAF, and predicted functions are presented in the Table 1. None of the genotype distributions deviated from the Hardy–Weinberg equilibrium as seen in Table 1.

Influence of Genetic Variability on the Risk for Parkinson’s Disease Development

Risk analysis showed a significant association between VDR rs2228570 and PD development. Carriers of at least one VDR rs2228570 T allele had increased odds for development of the disease (OR = 2.75; 95% CI = 1.73–4.38: p < 0.001). The association remained significant after adjustment for age and gender (OR = 3.07; 95% CI = 1.86–5.08: p < 0.001). The other investigated polymorphisms were not associated with PD risk (Table 3).

TABLE 3

Table 3. Associations between VDR genotypes and the risk for Parkinson’s disease.

Influence of VDR Genetic Variability on the Symptomatology of Parkinson’s Disease

LED was used as a measure of PD severity at the point of enrollment. The VDR rs4516035 was nominally significantly associated with the dose of dopaminergic drugs needed to adequately manage PD symptoms. Carriers of at least one VDR rs4516035 C allele required higher LED than carriers of the TT genotype (p = 0.023). There were no significant or nominally significant associations detected in the additive model. After adjustment of the nominally significant association from the dominant model for gender and disease duration, only disease duration showed significant effect on the LED (p < 0.001). Complete data is presented in the Supplementary Table 3.

Influence of VDR Genetic Variability on the Occurrence of Adverse Events of Dopaminergic Treatment

We checked the association of selected SNPs with the occurrence of eight different AEs of dopaminergic treatment. We found significant and nominally significant associations with four AEs, namely VHs, orthostatic hypotension, peripheral edema, and MFs.

Carriers of the VDR rs1544410 GA genotype (OR = 0.43; 95% CI = 0.22–0.84; p = 0.013) or carriers of at least one VDR rs1544410 A allele (OR = 0.46; 95% CI = 0.25–0.86; p = 0.014) had decreased odds for development of VHs. In regards to orthostatic hypotension, carriers of the VDR rs739837 GG genotype had increased odds (OR = 2.23; 95% CI = 1.05–4.70; p = 0.036) for development of this AE, whereas carriers of the VDR rs731236 TC (OR = 0.52; 95% CI = 0.29–0.93; p = 0.026) and CC genotype (OR = 0.38; 95% CI = 0.16–0.91; p = 0.029), as well as carriers of at least one VDR rs731236 C allele (OR = 0.48; 95% CI = 0.28–0.84; p = 0.010) had decreased odds for development of orthostatic hypotension. VDR rs1544410 was also associated with orthostatic hypotension. Carriers of the VDR rs1544410 GA genotype (OR = 0.48; 95% CI = 0.26–0.86; p = 0.013), AA genotype (OR = 0.42; 95% CI = 0.18–0.99; p = 0.046), and carriers of at least one A allele (OR = 0.46; 95% CI = 0.26–0.81; p = 0.007) had decreased odds for development of this AE. Finally, carriers of the VDR rs731236 TC genotype had decreased odds for development of peripheral edema (OR = 0.43; 95% CI = 0.20–0.90; p = 0.026). The results for all the investigated SNPs and non-motor AEs are presented in Supplementary Tables 4–9.

In regards to the motor AEs, only MFs presented with a nominally significant association with the VDR genetic variability. Carriers of at least one VDR rs2228570 T allele had increased odds for development of MFs (OR = 2.01; 95% CI = 1.01–4.02; p = 0.047). All the data for motor AEs is presented in Supplementary Tables 10, 11.

Significant and nominally significant associations obtained in univariate analysis were then adjusted for significant clinical parameters in multivariate analysis. Clinical parameters associated with the occurrence of non-motor and motor AEs in our patient cohort were already reported in our previous study (Redenšek et al., 2019). These significant clinical parameters were used for adjustments of associations with the corresponding AEs. The results of the multivariate analysis are presented in the Table 4. All of the nominally significant or significant associations from the univariate analysis retained significance levels also after adjustments, except for the association of the VDR rs2228570 with MFs. Additionally, the association between VDR rs731236 and peripheral edema was not adjusted for any clinical parameter because none of them were significantly associated with this AE. Both associations between the occurrence of VHs and the VDR rs1544410 genotype remained significant with p-values dropping even lower compared to the univariate analysis. Additionally, the majority of associations between VDR genetic variability and the occurrence of orthostatic hypotension remained nominally significant. Carriers of the VDR rs739837 GG genotype (OR = 2.23; 95% CI = 1.05–4.71; p = 0.036) had increased odds for development of this AE. Carriers of the VDR rs731236 TC genotype (OR = 0.52; 95% CI = 0.29–0.93; p = 0.028), CC genotype (OR = 0.39; 95% CI = 0.16–0.94; p = 0.035), or at least one C allele (OR = 0.49; 95% CI = 0.28–0.85; p = 0.011) had decreased odds for development of orthostatic hypotension. Finally, carriers of the VDR rs1544410 GA genotype (OR = 0.48; 95% CI = 0.27–0.86; p = 0.014) had decreased chances for this AE. The remaining two associations showed p-values above 0.05.

TABLE 4

Table 4. Significant associations from the univariate analysis adjusted for significant clinical parameters.

Discussion

The presented study assessed the associations between selected functional SNPs of VDR and symptomatology of PD. This is the first study to report strong associations between VDR polymorphisms and non-motor AEs of dopaminergic treatment, such as VHs and orthostatic hypotension. Additionally, it confirms the association of VDR rs2228570 with the risk of PD development. However, we observed no statistically significant associations with LED and consequently disease severity.

In our cohort, VDR rs2228570 conferred almost threefold increase in the risk for PD development in carriers of at least one T allele. The product of VDR rs2228570 wild-type C allele is a three amino-acids shorter version of the protein and is more transcriptionally potent than the longer version (Uitterlinden et al., 2004; Lv et al., 2020). This is even more important since VDR regulates the expression of several neurotrophic factors, such as nerve growth factor (Wion et al., 1991), neurotrophin 3 (Neveu et al., 1994), and GDNFs (Naveilhan et al., 1993). It has also been reported that VDR rs2228570 C allele carriers have higher capacity for intestinal calcium absorption, which may lead to higher vitamin D levels (Arai et al., 1997; Uitterlinden et al., 2004). Additionally, vitamin D signaling via VDR is involved in protection of neurons against oxidative stress (Fullard and Duda, 2020). Since lower vitamin D levels were reported as a risk factor for PD development several times (Knekt et al., 2010), we can expect that higher vitamin D levels could present a protective influence. Similarly, it has been shown that a significantly lower (50%) dose of 1,25-dihydroxyvitamin D₃ is necessary for the effect in the VDR rs2228570 CC genotype than in the VDR rs2228570 CT genotype carriers (Colin et al., 2000). A randomized control trial showed that vitamin D₃ supplementation may slow the progression of PD in patients with VDR rs2228570 CT and TT genotypes (Suzuki et al., 2013). The TT genotype was already associated with the more severe form of PD (Suzuki et al., 2012), with cognitive decline in PD (Gatto et al., 2015), and with PD risk (Hu et al., 2020; Agliardi et al., 2021). Some reported studies also showed an opposite effect in regards to PD risk, however most of these studies were done in Asian patient cohorts (Li et al., 2015; Wang et al., 2019), which may be related to the ethnic differences in genetic predisposition to diseases. All of the above indicates and supports a very high importance of the finding that VDR rs2228570 genotype partly defines the development and course of PD.

We observed an association between the VDR rs4516035 and LED as a measure of the general disease severity, which lost its statistical significance after adjustment for the disease duration. This variant modifies recognition sequence of the GATA-3 transcription factor (Fibla and Caruz, 2010), which might decrease expression of VDR and thus lead to worse PD symptomatology and consequently higher LED. However, we can conclude that disease duration is a much more important parameter in determining LED than VDR rs4516035 genotype, which is in concordance with the expectations. We have to emphasize that LED is not the most suitable parameter for assessment of PD severity in a particular point of the disease course. However, PD clinical rating scales were not available at the point of patient enrollment. Nevertheless, VDR genotype has been assessed in association with different PD-related phenotypes already (Suzuki et al., 2012; Gatto et al., 2015), which indicates that further research in independent patient cohorts is warranted.

We also detected some interesting novel associations of VDR SNPs with the occurrence of AEs of dopaminergic treatment. We observed that carriers of the VDR rs1544410 GA genotype and carriers of at least one A allele had decreased odds for development of VHs. It has been reported that the G allele is associated with lower mRNA stability, which indicates that there could be less VDR available for vitamin D signaling (Fibla and Caruz, 2010). Knowing this, higher VDR expression might protect against VHs. Additionally, knowing the genotype of patients in advance could enable compensation for lower VDR availability with vitamin D supplementation to avoid VHs development. For example, it has already been shown in Alzheimer’s disease that vitamin D supplementation protects against psychosis (Fan et al., 2020). Furthermore, VDR rs739837, rs731236, and rs1544410 genotypes were strongly associated with the occurrence of orthostatic hypotension. A study by Jang et al. (2015) showed that serum 25-hydroxyvitamin D₃ levels were significantly lower in PD patients with orthostatic hypotension in comparison to the patients without it, indicating an importance of vitamin D signaling in orthostatic hypotension occurrence. The exact mechanism is not known, but it was shown that vitamin D signaling plays a role in vascular endothelial function, blood pressure, renal function and cardiac autonomic activity (Soysal et al., 2014). We report that carriers of the VDR rs739837 GG genotype have more than twofold higher odds for development of orthostatic hypotension in comparison to the TT genotype. It was postulated that GG genotype leads to lower expression of VDR (Küchler et al., 2021), which might explain the increased risk for orthostatic hypotension. Furthermore, carriers of at least one VDR rs731236 C allele and VDR rs1544410 GA genotype had decreased odds for the development of orthostatic hypotension. Both genetic changes increase mRNA stability (Fibla and Caruz, 2010), which might be the reason for protective effect of the variants. Similarly, carriers of the VDR rs731236 TC genotype had decreased odds for the development of peripheral edema. The presented study is to the best of our knowledge the first one to report the association between VDR genetic variability and peripheral edema. However, the mechanism remains to be elucidated. The associations shown in this study all present with a clinical potential of pre-emptive identification of patients at risk for AE development and of patients that might benefit from vitamin D supplementation.

Our study presents important novel findings, but it should be emphasized that there are some limitations that we have to bear in mind when interpreting our results. A study cohort is of moderate size, but still comparable to other pharmacogenetic studies performed in the cohorts of PD patients and of a uniform genetic background as well. As explained above, LED is not the most suitable measure of disease severity. Clinical scales would be more relevant for this assessment. Unfortunately, this data was not available at the point of patient enrollment. The control group within the PD susceptibility analysis is in average almost 20 years younger than the patient group, which means that we might have missed some PD cases. However, the prevalence of PD is rather low in the general population, which is why we speculate that the age difference between the groups does not affect the results to a significant extent. Finally, we only considered the AEs as binary categorical variables since clinical scales are not routinely used to assess AEs in PD. Quantified severity of AEs would likely give us an opportunity to analyze associations more in depth.

Despite its limitations, the presented study provided us with important findings, which should nevertheless be assessed in independent populations as well. The association between VDR genetic variability and occurrence of orthostatic hypotension is a highly relevant finding which has not been elucidated yet. This indicates that patients with a specific VDR genotype might benefit from vitamin D supplementation to prevent orthostatic hypotension. VDR genotype could serve as a valuable predictive biomarker of this AE as well. Furthermore, we confirmed a very strong association between VDR rs2228570 and increased risk for PD development. This has recently been pointed out by several studies conducted in the Caucasian population from which a cohort of Slovenian patients included in the reported study originate as well. The latter confirms the reliability of this finding, which should be accounted for in future polygenic risk scores for prediction of PD since none of them so far take into account the VDR genotype (Ibanez et al., 2017; Dehestani et al., 2021). Additionally, patients with at risk VDR genotypes are suitable candidates for close monitoring of vitamin D levels.

Conclusion

This is the first study that evaluated a comprehensive spectrum of PD-related phenotypes in one cohort of PD patients in relation to VDR genetic variability. Our results may help identify people at risk to develop PD as well as PD patients at risk to develop certain AEs of dopaminergic treatment that would benefit from close monitoring of vitamin D levels and vitamin D supplementation. Results of the present study should be validated in an independent cohort. Nevertheless, we believe that these results will contribute to a more personalized treatment strategy in individual PD patients.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author/s.

Ethics Statement

The studies involving human participants were reviewed and approved by the Medical Ethics Committee of the Republic of Slovenia. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

SR, MT, and VD: conceptualization. SR and VD: methodology. SR, TK, and TB: formal analysis. SR: statistical analysis. SR and TK: writing—original draft preparation. SR, TK, TB, MT, and VD: writing—review and editing. SR and TB: visualization. MT and VD: supervision. VD: funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Slovenian Research Agency (ARRS) (grant no. P1-0170).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank all the patients for their participation in the study. We also thank Dušan Flisar, MD, Maja Kojovič, MD, Ph.D., Milica Gregorič-Kramberger, MD, Ph.D., Zvezdan Pirtošek, MD, Ph.D., Dejan Georgiev, MD, Ph.D., and Nina Zupančič Križnar, MD, for their contribution to patient enrollment.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2022.853277/full#supplementary-material

Abbreviations

VDR, vitamin D receptor; PD, Parkinson’s disease; GDNF, glial cell line-derived neurotrophic factor; MDS-UPDRS, MDS-Unified Parkinson’s Disease Rating Scale; H&Y, Hoehn and Yahr; SNP, single nucleotide polymorphism; AE, adverse event; DA, dopamine agonist; LED, levodopa equivalent dose; MFs, motor fluctuations; VHs, visual hallucinations; ICDs, impulse control disorders; ORs, odds ratios; CIs, confidence intervals.

References

Agliardi, C., Guerini, F. R., Zanzottera, M., Bolognesi, E., Meloni, M., Riboldazzi, G., et al. (2021). The VDR FokI (rs2228570) polymorphism is involved in Parkinson’s disease. J. Neurol. Sci. 428:117606.

Google Scholar

Arai, H., Miyamoto, K., Taketani, Y., Yamamoto, H., Iemori, Y., Morita, K., et al. (1997). A vitamin D receptor gene polymorphism in the translation initiation codon: effect on protein activity and relation to bone mineral density in Japanese women. J. Bone Miner. Res. 12, 915–921. doi: 10.1359/jbmr.1997.12.6.915

PubMed Abstract | CrossRef Full Text | Google Scholar

Chitsaz, A., Maracy, M., Basiri, K., Izadi Boroujeni, M., Tanhaei, A. P., Rahimi, M., et al. (2013). 25-hydroxyvitamin d and severity of Parkinson’s disease. Int. J. Endocrinol. 2013:689149.

Google Scholar

Colin, E. M., Weel, A. E., Uitterlinden, A. G., Buurman, C. J., Birkenhäger, J. C., Pols, H. A., et al. (2000). Consequences of vitamin D receptor gene polymorphisms for growth inhibition of cultured human peripheral blood mononuclear cells by 1, 25-dihydroxyvitamin D3. Clin. Endocrinol. 52, 211–216. doi: 10.1046/j.1365-2265.2000.00909.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Dehestani, M., Liu, H., and Gasser, T. (2021). Polygenic Risk scores contribute to personalized medicine of Parkinson’s disease. J. Pers. Med. 11:1030. doi: 10.3390/jpm11101030

PubMed Abstract | CrossRef Full Text | Google Scholar

DeLuca, G. C., Kimball, S. M., Kolasinski, J., Ramagopalan, S. V., and Ebers, G. C. (2013). Review: the role of vitamin D in nervous system health and disease. Neuropathol. Appl. Neurobiol. 39, 458–484. doi: 10.1111/nan.12020

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, H., Dhima, K., Lockhart, K. C., Locascio, J. J., Hoesing, A. N., Duong, K., et al. (2013). Unrecognized vitamin D3 deficiency is common in Parkinson disease: Harvard Biomarker Study. Neurology 81, 1531–1537. doi: 10.1212/wnl.0b013e3182a95818

PubMed Abstract | CrossRef Full Text | Google Scholar

Duval, G. T., Brangier, A., Barré, J., Launay, C. P., Beauchet, O., and Annweiler, C. (2015). Vitamin D Deficiency and Incident onset of orthostatic hypotension in older adults: preliminary results from the ‘MERE’, Study. J. Am. Geriatr. Soc. 63, 1245–1247. doi: 10.1111/jgs.13504

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, P., Qi, X., Sweet, R. A., and Wang, L. (2020). Network Systems pharmacology-based mechanism study on the beneficial effects of Vitamin D against Psychosis in Alzheimer’s Disease. Sci. Rep. 10:6136. doi: 10.1038/s41598-020-63021-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Fibla, J., and Caruz, A. (2010). “Vitamin D and HIV Infection,” in Soluble Factors Mediating Innate Immune Responses to HIV Infection, ed. M. Alfano (Sharjah: Bentham Science Publishers), 111–144. doi: 10.2174/978160805006211001010111

CrossRef Full Text | Google Scholar

Fullard, M. E., and Duda, J. E. (2020). A review of the relationship between Vitamin D and Parkinson disease symptoms. Front. Neurol. 11:454. doi: 10.3389/fneur.2020.00454

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, J., Teng, J., Liu, Z., Cai, M., and Xie, A. (2020). Association between vitamin D receptor polymorphisms and susceptibility to Parkinson’s disease: an updated meta-analysis. Neurosci. Lett. 720:134778. doi: 10.1016/j.neulet.2020.134778

PubMed Abstract | CrossRef Full Text | Google Scholar

Gatto, N. M., Sinsheimer, J. S., Cockburn, M., Escobedo, L. A., Bordelon, Y., and Ritz, B. (2015). Vitamin D receptor gene polymorphisms and Parkinson’s disease in a population with high ultraviolet radiation exposure. J. Neurol. Sci. 352, 88–93. doi: 10.1016/j.jns.2015.03.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Goetz, C. G., Poewe, W., Dubois, B., Schrag, A., Stern, M. B., Lang, A. E., et al. (2008). MDS-UPDRS. Available online at: http://www.movementdisorders.org/MDS-Files1/PDFs/MDS-UPDRS-Rating-Scales/NewUPDRS7308final.pdf (accessed September 18, 2017).

Google Scholar

Han, X., Xue, L., Li, Y., Chen, B., and Xie, A. (2012). Vitamin D receptor gene polymorphism and its association with Parkinson’s disease in Chinese Han population. Neurosci. Lett. 525, 29–33. doi: 10.1016/j.neulet.2012.07.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, W., Wang, L., Chen, B., and Wang, X. (2020). Vitamin D receptor rs2228570 polymorphism and Parkinson’s disease risk in a Chinese population. Neurosci. Lett. 717:134722. doi: 10.1016/j.neulet.2019.134722

PubMed Abstract | CrossRef Full Text | Google Scholar

Ibanez, L., Dube, U., Saef, B., Budde, J., Black, K., Medvedeva, A., et al. (2017). Parkinson disease polygenic risk score is associated with Parkinson disease status and age at onset but not with alpha-synuclein cerebrospinal fluid levels. BMC Neurol. 17:198. doi: 10.1186/s12883-017-0978-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Jang, W., Park, J., Kim, J. S., Youn, J., Oh, E., Kwon, K. Y., et al. (2015). Vitamin D deficiency in Parkinson’s disease patients with orthostatic hypotension. Acta Neurol. Scand. 132, 242–250. doi: 10.1111/ane.12390

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalia, L. V., and Lang, A. E. (2015). Parkinson’s disease. Lancet 386, 896–912.

Google Scholar

Kenborg, L., Lassen, C. F., Ritz, B., Schernhammer, E. S., Hansen, J., Gatto, N. M., et al. (2011). Outdoor work and risk for Parkinson’s disease: a population-based case-control study. Occup. Environ. Med. 68, 273–278. doi: 10.1136/oem.2010.057448

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, J. E., Oh, E., Park, J., Youn, J., Kim, J. S., and Jang, W. (2018). Serum 25-hydroxyvitamin D3 level may be associated with olfactory dysfunction in de novo Parkinson’s disease. J. Clin. Neurosci. 57, 131–135. doi: 10.1016/j.jocn.2018.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, J. S., Kim, Y. I., Song, C., Yoon, I., Park, J. W., Choi, Y. B., et al. (2005). Association of vitamin D receptor gene polymorphism and Parkinson’s disease in Koreans. J. Korean Med. Sci. 20, 495–498. doi: 10.3346/jkms.2005.20.3.495

PubMed Abstract | CrossRef Full Text | Google Scholar

Knekt, P., Kilkkinen, A., Rissanen, H., Marniemi, J., Sääksjärvi, K., and Heliövaara, M. (2010). Serum vitamin D and the risk of Parkinson disease. Arch. Neurol. 67, 808–811.

Google Scholar

Kravietz, A., Kab, S., Wald, L., Dugravot, A., Singh-Manoux, A., Moisan, F., et al. (2017). Association of UV radiation with Parkinson disease incidence: a nationwide French ecologic study. Environ. Res. 154, 50–56. doi: 10.1016/j.envres.2016.12.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Küchler, E. C., Schröder, A., Teodoro, V. B., Nazet, U., Scariot, R., Spanier, G., et al. (2021). The role of 25-hydroxyvitamin-D3 and vitamin D receptor gene in human periodontal ligament fibroblasts as response to orthodontic compressive strain: an in vitro study. BMC Oral Health 21:386. doi: 10.1186/s12903-021-01740-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Qi, H., Wei, S., Wang, L., Fan, X., Duan, S., et al. (2015). Vitamin D receptor gene polymorphisms and the risk of Parkinson’s disease. Neurol. Sci. 36, 247–255.

Google Scholar

Luo, X., Ou, R., Dutta, R., Tian, Y., Xiong, H., and Shang, H. (2018). Association between Serum Vitamin D levels and Parkinson’s Disease: a systematic review and meta-analysis. Front. Neurol. 9:909. doi: 10.3389/fneur.2018.00909

PubMed Abstract | CrossRef Full Text | Google Scholar

Lv, L., Tan, X., Peng, X., Bai, R., Xiao, Q., Zou, T., et al. (2020). The relationships of vitamin D, vitamin D receptor gene polymorphisms, and vitamin D supplementation with Parkinson’s disease. Transl. Neurodegener. 9:34. doi: 10.1186/s40035-020-00213-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Moretti, R., Morelli, M. E., and Caruso, P. (2018). Vitamin D in Neurological diseases: a rationale for a pathogenic impact. Int. J. Mol. Sci. 19:2245. doi: 10.3390/ijms19082245

PubMed Abstract | CrossRef Full Text | Google Scholar

Naveilhan, P., Neveu, I., Baudet, C., Ohyama, K. Y., Brachet, P., and Wion, D. (1993). Expression of 25(OH) vitamin D3 24-hydroxylase gene in glial cells. Neuroreport 5, 255–257. doi: 10.1097/00001756-199312000-00018

PubMed Abstract | CrossRef Full Text | Google Scholar

Neveu, I., Naveilhan, P., Jehan, F., Baudet, C., Wion, D., De Luca, H. F., et al. (1994). 1,25-dihydroxyvitamin D3 regulates the synthesis of nerve growth factor in primary cultures of glial cells. Brain Res. Mol. Brain Res. 24, 70–76. doi: 10.1016/0169-328x(94)90119-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Niu, M. Y., Wang, L., and Xie, A. M. (2015). ApaI, BsmI, FokI, and TaqI Polymorphisms in the Vitamin D Receptor Gene and Parkinson’s Disease. Chin. Med. J. 128, 1809–1814. doi: 10.4103/0366-6999.159358

PubMed Abstract | CrossRef Full Text | Google Scholar

Peterson, A. L., Murchison, C., Zabetian, C., Leverenz, J. B., Watson, G. S., Montine, T., et al. (2013). Memory, mood, and vitamin D in persons with Parkinson’s disease. J. Parkinsons Dis. 3, 547–555. doi: 10.3233/JPD-130206

PubMed Abstract | CrossRef Full Text | Google Scholar

Poewe, W., Seppi, K., Tanner, C. M., Halliday, G. M., Brundin, P., Volkmann, J., et al. (2017). Parkinson disease. Nat. Rev. Dis. Primers 3:17013.

Google Scholar

Redenšek, S., Flisar, D., Kojoviæ, M., Gregoriè Kramberger, M., Georgiev, D., Pirtošek, Z., et al. (2019). Dopaminergic Pathway genes influence adverse events related to dopaminergic treatment in Parkinson’s Disease. Front. Pharmacol. 10:8. doi: 10.3389/fphar.2019.00008

PubMed Abstract | CrossRef Full Text | Google Scholar

Samuel, S., and Sitrin, M. D. (2008). Vitamin D’s role in cell proliferation and differentiation. Nutr. Rev. 66(10 Suppl. 2), 1753–4887. doi: 10.1111/j.1753-4887.2008.00094.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Soysal, P., Yay, A., and Isik, A. T. (2014). Does vitamin D deficiency increase orthostatic hypotension risk in the elderly patients? Arch. Gerontol. Geriatr. 59, 74–77. doi: 10.1016/j.archger.2014.03.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Suzuki, M., Yoshioka, M., Hashimoto, M., Murakami, M., Kawasaki, K., Noya, M., et al. (2012). 25-hydroxyvitamin D, vitamin D receptor gene polymorphisms, and severity of Parkinson’s disease. Mov. Disord. 27, 264- -271. doi: 10.1002/mds.24016

PubMed Abstract | CrossRef Full Text | Google Scholar

Suzuki, M., Yoshioka, M., Hashimoto, M., Murakami, M., Noya, M., Takahashi, D., et al. (2013). Randomized, double-blind, placebo-controlled trial of vitamin D supplementation in Parkinson disease. Am. J. Clin. Nutr. 97, 1004–1013. doi: 10.3945/ajcn.112.051664

PubMed Abstract | CrossRef Full Text | Google Scholar

Tolosa, E., Garrido, A., Scholz, S. W., and Poewe, W. (2021). Challenges in the diagnosis of Parkinson’s disease. Lancet Neurol. 20, 385–397.

Google Scholar

Tomlinson, C. L., Stowe, R., Patel, S., Rick, C., Gray, R., and Clarke, C. E. (2010). Systematic review of levodopa dose equivalency reporting in Parkinson’s disease. Mov. Disord. 25, 2649–2653. doi: 10.1002/mds.23429

PubMed Abstract | CrossRef Full Text | Google Scholar

Török, R., Török, N., Szalardy, L., Plangar, I., Szolnoki, Z., Somogyvari, F., et al. (2013). Association of vitamin D receptor gene polymorphisms and Parkinson’s disease in Hungarians. Neurosci. Lett. 551, 70–74. doi: 10.1016/j.neulet.2013.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Uitterlinden, A. G., Fang, Y., Van Meurs, J. B., Pols, H. A., and Van Leeuwen, J. P. (2004). Genetics and biology of vitamin D receptor polymorphisms. Gene 338, 143–156. doi: 10.1016/j.gene.2004.05.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Yang, D., Yu, Y., Shao, G., and Wang, Q. (2016). Vitamin D and sunlight exposure in newly-diagnosed Parkinson’s Disease. Nutrients 8:142. doi: 10.3390/nu8030142

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Shen, N., Lu, Y., and Tan, K. (2019). Vitamin D receptor polymorphisms and the susceptibility of Parkinson’s disease. Neurosci. Lett. 699, 206–211. doi: 10.1016/j.neulet.2019.02.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilkins, C. H., Sheline, Y. I., Roe, C. M., Birge, S. J., and Morris, J. C. (2006). Vitamin D deficiency is associated with low mood and worse cognitive performance in older adults. Am. J. Geriatr. Psychiatry 14, 1032–1040. doi: 10.1097/01.JGP.0000240986.74642.7c

PubMed Abstract | CrossRef Full Text | Google Scholar

Wion, D., MacGrogan, D., Neveu, I., Jehan, F., Houlgatte, R., and Brachet, P. (1991). 1,25-Dihydroxyvitamin D3 is a potent inducer of nerve growth factor synthesis. J. Neurosci. Res. 28, 110–114. doi: 10.1002/jnr.490280111

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Z., and Taylor, J. A. (2009). SNPinfo: integrating GWAS and candidate gene information into functional SNP selection for genetic association studies. Nucleic Acids Res. 37, W600–W6005. doi: 10.1093/nar/gkp290

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, N., Yang, Y., Zhang, L., Zhang, Q., Balbuena, L., Ungvari, G. S., et al. (2021). Quality of life in Parkinson’s disease: a systematic review and meta-analysis of comparative studies. CNS Neurosci. Ther. 27, 270–279. doi: 10.1111/cns.13549

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, D., Liu, G. Y., Lv, Z., Wen, S. R., Bi, S., and Wang, W. Z. (2014). Inverse associations of outdoor activity and vitamin D intake with the risk of Parkinson’s disease. J. Zhejiang Univ. Sci. B 15, 923–927. doi: 10.1631/jzus.B1400005

PubMed Abstract | CrossRef Full Text | Google Scholar