Peripheral Clock System Abnormalities in Patients With Parkinson’s Disease

Abstract

Objective: To evaluate the altered expression of peripheral clock genes, circulating melatonin levels, and their correlations with sleep-wake phenotypes including probable rapid eye movement sleep behavior disorder (pRBD) symptoms in a relatively large sample of Parkinson’s disease (PD) patients. Methods: We determined the expression profiles of five principal clock genes, BMAL1, CLOCK, CRY1, PER1, and PER2, in the peripheral blood mononuclear cells (PBMCs) of PD patients (n=326), and healthy controls (HC, n=314) using real-time quantitative PCR. Melatonin concentrations in the plasma of two groups were evaluated by enzyme-linked immunosorbent assay. Then we performed comprehensive association analyses on the PBMCs clock gene expression, plasma melatonin levels and sleep characteristics. Results: Our data showed that the expression levels of BMAL1, CLOCK, CRY1, PER1, and PER2 were significantly decreased in the PBMCs of PD as compared with that of HC (P< 0.05). PD patients had reduced plasma melatonin levels compared with HC (P< 0.0001). pRBD and excessive daytime sleepiness are common in these PD patients and are associated with the expression levels of all five clock genes (r= -0.344~-0.789, P< 0.01) and melatonin concentrations (r= -0.509~-0.753, P< 0.01). Statistical analyses also revealed that a combination of five clock genes and melatonin could reach a high diagnostic performance (areas under the curves, 97%) for PD comorbid pRBD. Conclusions: This case-control study demonstrates that peripheral BMAL1, CLOCK, CRY1, PER1, PER2 and melatonin levels are altered in PD patients and may serve as endogenous markers for sleep and wakefulness disturbances of PD.

Introduction

Parkinson’s disease (PD) is a common and progressive neurodegenerative disorder that manifests with a broad range of motor and non-motor symptoms (Kalia and Lang, 2015; Schapira et al., 2017). Among a variety of non-motor manifestations of PD, disturbance of sleep and wakefulness is recognized as the most common symptoms affecting as many as 90% of PD patients (Tandberg et al., 1998; Videnovic et al., 2014). Sleep-wake disturbances mainly manifest as poor night-time sleep quality, excessive daytime sleepiness (EDS) and rapid eye movement sleep behavior disorder (RBD) (Li et al., 2021), which ultimately lead to poor quality of life and increased mortality in PD patients (Tandberg et al., 1998; Mattis and Sehgal, 2016).

The sleep-wake cycle is regulated by circadian rhythms (Mattis and Sehgal, 2016). Circadian rhythm is a physiological and behavioral process with a periodicity approximately every 24 h, driven by the endogenous circadian clock in the suprachiasmatic nucleus (SCN), the central circadian pacemaker (Edgar et al., 2012). In addition to sleep-wake disturbances, other circadian-related symptoms including motor activity, visual performance, autonomic dysfunction, and responsiveness to PD medications have also been shown daily fluctuations in PD patients (Videnovic and Willis, 2016; De et al., 2018). Besides, a growing body of research has reported that there exist complicated interactions between the circadian rhythm and dopaminergic systems (Korshunov et al., 2017; Sundaram et al., 2019). All this evidence highlights the complex characteristics of circadian disruption in PD patients, suggesting the need for a reliable and accurate evaluation to identify the main cause of circadian disruption in PD patients.

The circadian rhythm is regulated by a group of clock genes that coordinated their transcription and translation negative feedback loops at a cellular level (Rijo-Ferreira and Takahashi, 2019; Liu et al., 2021). Two of the most important transcription factors are called brain muscle aryl hydrocarbon receptor nuclear translocator-like protein 1 (BMAL1) and circadian locomotor output cycles kaput (CLOCK), which together activate the transcription of cryptochrome (CRY) and period (PER) genes (Hergenhan et al., 2020). Clock genes are not only expressed in the central circadian pacemaker but also in most of the peripheral tissues and cells (Hergenhan et al., 2020). It has been shown that the clock machinery is rhythmically expressed in the human peripheral blood mononuclear cells (PBMCs) (Johansson et al., 2016; Zhang S. et al., 2020). Besides, the SCN also transmits its circadian signal to the pineal gland, the organ that secretes melatonin during the dark phase (Wu and Swaab, 2005). Circulating melatonin is a well-accepted marker of endogenous circadian rhythm and its level is altered in patients with PD (Videnovic et al., 2014; Videnovic and Willis, 2016). Therefore, expression levels of clock genes in the PBMCs along with circulating melatonin may provide new avenues to evaluate the circadian abnormalities in PD patients.

There are only a few case studies of circadian rhythm abnormalities in a small number of PD patients in recent years (Fifel, 2017). The objective of this study is therefore to determine whether the expression levels of five principal clock genes, BMAL1, CLOCK, CRY1, PER1, and PER2 in the PBMCs and circulating melatonin levels are altered in PD patients compared to healthy controls (HC). Comprehensive analyses on the associations between the expression levels of PBMCs clock genes, plasma melatonin concentration and the sleep disturbances phenotypes of PD patients have also been performed in this study.

Materials and Methods

Participants

A total of 640 participants (326 PD patients and 314 HC) were enrolled between January 2017 and June 2020 (Table 1). Patients with PD in this study were recruited from those who visited the Neurology Department of the First Affiliated Hospital of Dalian Medical University. The PD patients were diagnosed as idiopathic PD according to the Movement Disorder Society Clinical Diagnostic Criteria for PD (Postuma et al., 2015). Among 326 PD patients, 107 were of recent-onset without any PD treatment, the other 219 patients were treated with PD medications. PD disease severity was assessed by Modified Hoehn and Yahr (H-Y) staging (Goetz et al., 2004). Participants of HC were enrolled from the Health Examination Center of Dalian Medical University, showing that they did not have any obvious neurological disorders or non-neurological disorders. Severe depression and psychosis, and other conditions which may affect the circadian rhythms such as melatonin treatment, shiftwork traveling to different time zones within 30 days before the enrollment, were excluded from all the participants. This study has been granted ethical approval by the Ethics Committee of the First Affiliated Hospital of Dalian Medical University (approval number: LCKY2014-29). Written informed consent was obtained from all participants.

TABLE 1

Table 1. Demographics, clinical and sleep characteristics of HC and patients with PD.

Sleep-Related Measurement Scales

The severity of sleep and wakefulness disturbances in PD patients was assessed by the sleep-related measurement scales including the RBD questionnaire–Hong Kong (RBDQ-HK), Epworth Sleepiness Scale (ESS), and Pittsburgh Sleep Quality Index (PSQI). The RBDQ-HK is a Chinese instrument for screening and modifying the severity of RBD (Dauvilliers et al., 2008). Since patients in this study had not received a polysomnography diagnosis, those who have the RBDQ-HK score of more than 18 (the cutoff point) were determined as the probable RBD (pRBD) (Shen et al., 2014; Wong et al., 2016). Daytime sleepiness was measured by the Chinese version of the ESS. A total score ≥ 10 indicates EDS (Videnovic et al., 2014). PSQI was used in the assessment of the night sleep quality of participants (Mollayeva et al., 2016).

Blood Sampling, Plasma, and Peripheral Blood Mononuclear Cells Separation

When the blood sampling, all the participants (PD patients and HC) were not strictly shielded from external light with the constant ambient temperature at 21 ± 2°C. Fasting peripheral blood samples (2 mL) were collected by direct venipuncture at 6–7 a.m. The plasma was aliquoted (800 μL) into a sterile tube and stored at −80°C. PBMCs were separated from the Peripheral Lymphocyte Separation Medium (HAOYANG, Tianjin, China) by centrifugation at 450 g for 20 min and were stored at −80°C.

Peripheral Blood Mononuclear Cells mRNA Extraction and Quantification

Total RNAs from the PBMCs were extracted utilizing the mirVana mRNA Isolation Kit (Ambion, Carlsbad, CA, United States). The target five clock genes mRNA levels were then measured by quantitative real-time PCR (qRT-PCR). The PCR reactions were executed using ABI 7500 fast real-time PCR system (Applied Biosystems, Foster City, CA, United States) in a total volume of 20 μL for each reaction. GAPDH gene was used as the housekeeping gene. The specific primers targeting PBMCs BMAL1, CLOCK, CRY1, PER1, and PER2 are presented in Supplementary Table 1. After 94°C for 30 s, the experimental reaction consisted of 40 cycles of 94°C for 5 s and 60°C for 34 s, and the target genes were detected by the fluorescent dye SYBR Green I (TransGen, Beijing, China). The target genes expression was determined using the 2^{–delta Ct} method.

Plasma Melatonin Measurement

A total of 309 plasma samples from the participants in this study (156 from HC and 153 from PD) were separated simultaneously with the extraction of PBMCs. Plasma melatonin concentration was evaluated using the Human Melatonin enzyme-linked immunosorbent assay (ELISA) Kit (Elabscience, Wuhan, China) according to the manufacturer’s instructions. The concentration of melatonin in the plasma was finally determined by comparing the OD values of the samples to the standard curves. The analytical sensitivity of the assay was 9.38 pg/mL. The intra-assay coefficients of variation for the kits were 9.4% (PD) and 9.5% (HC), and the inter-assay variation coefficients were 9.2% (PD) and 9.5% (HC).

Statistical Analysis

Data were examined for normality using the Shapiro-Wilk test. Quantitative data were expressed as mean ± standard error (SE) for each group. We analyzed group differences using the Kruskal-Wallis and non-parametric Mann-Whitney U-test. Spearman’s coefficients were calculated to evaluate bivariate associations between the demographic, disease severity, sleep characteristics, plasma melatonin concentration and the expression levels of clock genes. Partial correlation analyses were utilized to adjust for age, sex, and pharmacotherapy. Receiver operating characteristic (ROC) curves and areas under the curves (AUC) were used to evaluate the diagnostic performance of the potential biomarkers. The alterations of the clock genes levels and melatonin concentration were evaluated using two-way ANOVA for the effect of sunlight. Statistical analyses were performed with the SPSS software version 22.0 (SPSS Inc., Chicago, IL, United States). A P-value < 0.05 was considered as statistical significance.

Results

Characteristics of the Study Population

In this study, we collected a total of 640 PBMCs samples: 326 from patients with idiopathic PD and 314 from HC. All participants enrolled in this study were Chinese. Demographics, disease status, and sleep characteristics of the study cohort are outlined in Table 1. No significant difference in sex or age was found between the PD and HC groups. In the 264 PD patients (80.98%) who returned the sleep-related questionnaires, RBDQ-HK scores ranged from 2 to 91, with 88 PD patients (33.33%) classified as patients with pRBD (RBDQ-HK score greater than 18). 160 patients (60.6%) had EDS, as defined by an ESS score of more than 9. PD patients with pRBD presented more serious EDS and poorer sleep quality (P < 0.0001). Similarly, patients with EDS also had more severe RBD symptoms and poorer sleep quality (P < 0.0001). Demographics, disease duration, and H-Y scores were not significantly different in the subgroups of PD patients (Table 1).

Altered Expressions of BMAL1, CLOCK, CRY1, PER1, and PER2 Genes in the Peripheral Blood Mononuclear Cells of Parkinson’s Disease Compared With Healthy Controls

We determined five clock genes (BMAL1, CLOCK, CRY1, PER1, and PER2) mRNA levels in the PBMCs of PD patients (n = 326) and HC (n = 314) by qRT-PCR. Our data showed significantly lower levels of BMAL1 (P < 0.0001), CLOCK (P < 0.01), CRY1 (P < 0.01), PER1 (P < 0.05), and PRY2 (P < 0.01) in PD than those in HC (Figures 1A–E).

FIGURE 1

Figure 1. Scatter plots of BMAL1, CLOCK, CRY1, PER1, and PER2 relative mRNA expression levels in the PBMCs of HC (n = 314) and PD (n = 326). (A–E) Significantly lower levels of BMAL1 (A, P < 0.0001), CLOCK (B, P < 0.01), CRY1 (C, P < 0.01), PER1 (D, P < 0.05) and PRY2 (E, P < 0.01) were seen in PD patients compared to those in HC. Horizontal bars represent mean and SE values. *P < 0.05, **P < 0.01, and ****P < 0.0001.

Among PD patients, those with pRBD (n = 88) had significantly lower expression levels of all five clock genes compared with those in HC (n = 314, P < 0.05, Figures 2A–E), in particular, they had lower expression levels of BMAL1 when comparing to the PD patients without pRBD (n = 176, P < 0.05, Figure 2A). PBMCs BMAL1 (P < 0.05), CLOCK (P < 0.05), and PER2 (P < 0.05) levels in PD patients without pRBD were also decreased markedly compared with those in HC (Figures 2A,B,E). Similarly, expression levels of all five clock genes in PD patients with EDS (n = 160) were significantly lower than those in HC (n = 314, P < 0.05, Figures 2F–J). Significant differences in BMAL1 (P < 0.01) and CLOCK (P < 0.05) expression levels were seen between the subgroups of PD patients with EDS and without EDS (Figures 2F,G). However, when comparing the five clock genes expression levels in the subgroups of patients without EDS and HC, no significant difference was observed.

FIGURE 2

Figure 2. Scatter plots of BMAL1, CLOCK, CRY1, PER1, and PER2 relative mRNA expression levels in the PBMCs of HC and different subgroups of PD. (A–E) PD patients with pRBD (n = 88) had significantly lower expression levels of five clock genes compared with those in HC (n = 314, BMAL1, P < 0.0001; CLOCK, P < 0.01; CRY1, P < 0.05; PER1, P < 0.05; PER2, P < 0.01, respectively), and had decreased level of BMAL1 compared with that of PD patients without pRBD (n = 176, P < 0.05). (F–J) Expression levels of all five clock genes in PD patients with EDS (n = 160) were significantly lower than that of HC (n = 314, BMAL1, P < 0.0001; CLOCK, P < 0.001; CRY1, P < 0.01; PER1, P < 0.05; PER2, P < 0.01, respectively). Significantly differences of BMAL1 (P < 0.01) and CLOCK (P < 0.05) expression levels were seen between the subgroups of PD patients with EDS and without EDS. Horizontal bars represent mean and SE values. *P < 0.05, **P < 0.01, ***P < 0.001, and ****P < 0.0001.

We then analyzed associations between the expression levels of clock genes and the clinical characteristics of PD patients. Correlation coefficients showed significant associations among the expression levels of five clock genes (Table 2). Expression levels of CLOCK (r = −0.143, P < 0.05) and CRY1 (r = −0.12, P < 0.05) were inversely associated with the age of the PD patients (Table 2). However, disease duration (years after the onset of disease symptoms) and disease severity (the H-Y score) in the PD group were not significantly correlated with the expression levels of the five clock genes (correcting for age, sex, and pharmacotherapy, Table 2). In the analysis of the relationships between PD pharmacotherapy and clock genes expressions, no significant difference of the five clock genes levels was found among the four different pharmacotherapy subgroups of PD (correcting for age and sex, Supplementary Figure 1A).

TABLE 2

Table 2. The correlation coefficients between the expression levels of clock genes in the PBMCs, plasma melatonin concentration, the age and sex of the PD patients, and HC.

Altered Concentration of Melatonin in the Plasma of Parkinson’s Disease Compared With Healthy Controls

Results from ELISA show that plasma melatonin concentration was significantly decreased in patients with PD than in HC (P < 0.0001, Figure 3A). PD patients with pRBD (n = 56) had significantly lower melatonin concentration compared with those in HC (n = 314, P < 0.0001) and PD patients without pRBD (n = 97, P < 0.0001). Besides, there is a significant difference between melatonin concentration in HC and PD patients without pRBD (P < 0.0001, Figure 3B). Similar alterations of plasma melatonin concentration were also found between EDS group (n = 92) and negative EDS group (n = 64) in PD patients (P < 0.0001), and both of two groups have significant lower melatonin concentration when comparing to HC (PD with EDS, P < 0.0001; PD without EDS, P < 0.01; Figure 3C). Spearman and partial correlation analyses show that there are no significant associations between melatonin concentration, age, sex, disease duration or disease severity (Tables 2, 3). Furthermore, plasma melatonin concentration in four different pharmacotherapy subgroups of PD didn’t reach a statistical difference (Supplementary Figure 1B).

FIGURE 3

Figure 3. Scatter plots of melatonin concentration in the plasma of HC (n = 156) and PD (n = 153). (A) Plasma melatonin concentration was significantly decreased in patients with PD than in HC (P < 0.0001). (B) PD patients with pRBD (n = 56) had significantly lower melatonin concentration compared with those in HC (n = 156, P < 0.0001) and PD patients those without pRBD (n = 97, P < 0.0001). A significant difference in melatonin concentration was found between HC and PD without pRBD (P < 0.0001). (C) Alterations of plasma melatonin concentration were found between PD with EDS (n = 92) and PD without EDS (n = 64, P < 0.0001), and both of two groups have significant lower melatonin concentration when comparing to HC (PD with EDS, P < 0.0001; PD without EDS, P < 0.01). ∗∗P < 0.01 and ****P < 0.0001.

TABLE 3

Table 3. The correlation coefficients between the expression levels of clock genes in the PBMCs, sleep, and disease characteristics of the PD patients.

Associations Between the Expression Levels of Clock Genes, Plasma Melatonin, and Sleep-Wake Disturbances in Parkinson’s Disease Patients

We then determined the correlations between the severity of pRBD (RBDQ-HK score), the severity of daytime sleepiness (ESS score), self-reported sleep quality (PSQI score), PBMCs levels of five clock genes, and plasma melatonin levels of PD patients (correcting for the age, sex, and pharmacotherapy). Our results showed that the RBDQ-HK, ESS, and PSQI scores completed by the PD patients (n = 264) were inversely associated with the expression levels of all five clock genes (r = −0.344∼−0.789, P < 0.01) and plasma melatonin (r = −0.509∼−0.753, P < 0.01, Table 3).

The performance of PBMCs clock genes and plasma melatonin for the diagnosis of PD comorbid pRBD and EDS were evaluated by the AUC values based on the ROC curve analysis. The results showed that AUCs of the combination of five clock genes and melatonin for pRBD were 0.97 (95% CI, 0.94–0.99, P < 0.01), which outperformed those of five clock genes and melatonin alone (Figure 4A). Similarly, the AUC value of the combined five clock genes with melatonin was 0.91 (95% CI, 0.86–0.96, P < 0.01) for EDS, which reached a higher performance compared to these factors alone (Figure 4B).

FIGURE 4

Figure 4. Receiver operating characteristic (ROC) curves of BMAL1, CLOCK, CRY1, PER1, PER2 and melatonin levels for PD with pRBD vs. PD without pRBD, and those for PD with EDS vs. PD without EDS. (A) The AUCs values of BMAL1, CLOCK, CRY1, PER1, PER2, melatonin, and the combination of five clock genes and melatonin for pRBD were 0.94 (BMAL1, 95% CI, 0.9–0.97, P < 0.05), 0.8 (CLOCK, 95% CI, 0.73–0.88, P < 0.05), 0.81 (CRY1, 95% CI, 0.74–0.88, P < 0.05), 0.76 (PER1, 95% CI, 0.68–0.84, P < 0.05), 0.8 (PER2, 95% CI, 0.73–0.87, P < 0.05), 0.82 (melatonin 95% CI, 0.75–0.89, P < 0.05), and 0.97 (combined five clock genes and melatonin, 95% CI, 0.94–0.99, P < 0.05), respectively. (B) The AUCs values for EDS were 0.86 (BMAL1, 95% CI, 0.8–0.92, P < 0.05), 0.82 (CLOCK, 95% CI, 0.75–0.89, P < 0.05), 0.78 (CRY1, 95% CI, 0.71–0.86, P < 0.05), 0.73 (PER1, 95% CI, 0.65–0.81, P < 0.05), 0.7 (PER2, 95% CI, 0.61–0.78, P < 0.05), 0.83 (melatonin, 95% CI, 0.76–0.9, P < 0.05), and 0.91 (combined five clock genes and melatonin, 95% CI, 0.86–0.96, P < 0.05), respectively.

Altered Clock Genes and Melatonin With the Changes of Sunlight Intensity

Since we collected the peripheral blood samples from the subjects in this study at a specific time which is 6--7 a.m. but in different months from January 1, 2018, to March 31, 2020, which means that there are variations of the environmental sunlight intensity when sampling. To clarify the alterations of the clock genes as the sunlight changes in different months, subjects were further divided into 4 subgroups according to the annual sunrise schedule of Dalian city, China,1 including I: Sampling before the sunrise (blood collection in January, February, and December), II: Sampling around the sunrise (blood collection in March and November), III: Sampling 0.5 h after the sunrise (blood collection in April, September, and October), IV: Sampling 1 h after the sunrise (blood collection in May, June, July, and August) (Figure 5A). The results showed that all these four subgroups in PD presented lower levels of BMAL1, CLOCK, CRY1, PER1, PER2, and melatonin than those in HC (Figures 5B–G). In the group of HC, PBMCs levels of BMAL1 (P < 0.01), PER1 (P < 0.05) and plasma melatonin (P < 0.0001) exhibited significant variations with the enhancement of the sunlight intensity. However, there was a lack of light-dependent variation in BMAL1 (P = 0.41, Figure 5B), PER1 expression (P = 0.19, Figure 5E), and melatonin concentration (P = 0.48, Figure 5F) in PD patients. Down-regulation of CLOCK, CRY1, and PER2 expression with higher light intensity was also found in HC and PD patients; however, no statistically significance was found (Figures 5C,D,F).

FIGURE 5

Figure 5. Expression levels of the five clock genes and melatonin in different sunlight intensity of HC and patients with PD. (A) A schematic diagram of the sunlight conditions when sampling in different groups of months. The time of blood sampling from the subjects was fixed at 6–7 a.m. Subjects were divided into 4 groups according to the local sunrise schedule. I: Sampling before the sunrise (January, February, and December, n (HC: PD) = 79: 80); II: Sampling around the sunrise (March and November, n (HC: PD) = 58: 64); III: Sampling 0.5 h after the sunrise (April, September, and October, n (HC: PD) = 57: 63); IV: Sampling 1 h after the sunrise (May, June, July, and August; n (HC: PD) = 120: 119). (B–F) PBMCs levels of BMAL1 (B), CLOCK (C), CRY1 (D), PER1 (E), and PER2 (F) expression in the group of PD and HC, with different intensities of sunlight. (G) Plasma melatonin concentration in PD and HC with different intensities of sunlight. The results are the mean ± SEM values. *P < 0.05, **P < 0.01, and ****P < 0.0001.

Discussion

In this study, we measured the expression profiles of five principal clock genes (BMAL1, CLOCK, CRY1, PER1, and PER2) in the PBMCs, as well as the melatonin concentration in the plasma of PD patients and HC simultaneously. To our knowledge, it was the first study to analyze the expression levels of five clock genes from a relatively large population of PD patients and controls. The main finding from this study is that the peripheral molecular clock is altered in PD, which could be potential markers for the sleep and wakefulness disturbances of PD. Furthermore, this study reported for the first time the loss of light-dependent variations of BMAL1 and PER1 levels in PD patients.

The circadian rhythm in mammals is cell-autonomous depending on transcriptional autoregulatory feedback loops organized by a core set of clock genes (Zhang S. et al., 2020). BMAL1 is a central component of the mammalian circadian clock system (Yang et al., 2018). CLOCK as the first discovered mammalian clock gene, along with BMAL1 can form a heterodimer (BMAL1: CLOCK) that binds to regulatory elements (E-boxes) sites and induces the expression of other clock genes, including PER and CRY (Trott et al., 2018; Cox and Takahashi, 2019). The heterodimerized PER and CRY in the cytoplasm translocate to the nucleus, where they inhibit BMAL1: CLOCK-mediated transcription (Cox and Takahashi, 2019). Overall, these interlocked transcriptional feedback loops participate in the molecular foundation of circadian rhythms.

Prior studies of clock genes in PD patients are scarce. Similar to our findings, Cai et al. (2010) reported lower PBMCs levels of BMAL1 in PD (n = 68) compared with HC (n = 54). While Breen et al. (2014). found no significant difference between the BMAL1 levels in PD and HC at 6 a.m., there was an increased expression of PER2 at 4 a.m. in PD. This study also found that plasma melatonin levels were decreased in PD patients compared with controls, which is consistent with the previous reports (Bordet et al., 2003; Breen et al., 2014; Videnovic et al., 2014). Besides, expression levels of CLOCK and CRY1 in PD patients showed a significant association with age in PD patients but not in HC. It has been reported that aging is associated with the disturbance of the circadian system resulting in an internal desynchronization of rhythm (De et al., 2018). We suspect that circadian abnormalities could be more sensitive to aging in PD patients than that in HC. Although previous studies reported that signaling mediated by the dopamine activity could promote the transcriptional capacity of the BMAL1-CLOCK complex (Landgraf and Oster, 2012), we found no significant difference among the four different pharmacotherapy groups of PD patients in this study. Our results suggest that the molecular basis of altered circadian rhythm gene expression in PD is more likely to be the consequence of PD pathology. One recent brain bank case-control study found the disease-related neuropathological changes in the SCN of PD patients, indicating the pathophysiological mechanisms underlying the PD pathology and circadian rhythm abnormalities (Pablo-Fenandez et al., 2018). On the other hand, it has been proposed that circadian disruption could also influence the neurodegenerative process in PD (Musiek and Holtzman, 2016). Most notably, the alterations of clock genes and melatonin levels can be observed in the patients who were in the early stage of PD or recent-onset without any PD treatment, and these alterations were not significantly related to the disease severity and duration of PD. Our result suggests that the dysregulation of the peripheral molecular clock could be an early sign for PD, which helps to explain why the circadian-related symptoms, for instance, RBD and autonomic dysfunction, occur in the prodromal phase of PD (Liu and Le, 2020).

One of the most prominent outputs of the circadian system is the synchronization of the sleep-wake cycle (Landgraf and Oster, 2012). Our data further demonstrate that the disturbances of sleep and wakefulness, which are mainly manifested as RBD symptoms, EDS, and poor night sleep quality, are common in PD patients. Numerous studies have revealed RBD to be an original symptom and anticipant biomarker of PD (Li and Le, 2019; Zhang F. et al., 2020). Until now, the association between RBD and PD in the regulation of the circadian rhythm is still not fully understood. Given the cost and facilities required for the gold standard-polysomnography (PSG) for the diagnosis of RBD, questionnaires that adequately screen for pRBD are useful for clinical studies (Dauvilliers et al., 2008). This study using RBDQ-HK evaluation of pRBD documents suggested that there may be interactions between circadian rhythm disruption and the pathology of PD even without RBD. However, PD patients with comorbid RBD may have a severer disturbance of molecular circadian rhythm.

This study also provides important evidence on the potential of peripheral clock genes and melatonin for the identification of sleep and wakefulness disturbances in PD patients. This is probably because the altered clock genes and melatonin level may cause a disruption to the neural circuitry controlling circadian rhythms, which in turn leads to disturbances in the sleep-wake cycle of PD patients. There’s plenty of evidence showing that clock genes and the sleep-wake cycle are tightly linked (Wulff et al., 2009). Study on the BMAL1-deficient mice illustrated reduced sleep-wake rhythmicity and increased sleep fragmentation (Laposky et al., 2005). Besides, the mice lacking CRY1 also illustrated sleep structure alteration, including an increase in non-REM time and EEG delta power (Wisor et al., 2002). Gene association studies in two independent populations reported an association between sequence variants of CLOCK and sleep duration (Allebrandt et al., 2010). Polymorphisms of the human PER1 and PER2 gene has also been found in association with abnormal sleep pattern (Toh et al., 2001; Carpen et al., 2006). In addition to the clock genes, Breen et al. also found that alterations in melatonin output were significantly associated with reduced slow-wave and REM sleep (Maggio et al., 2019). Our results indicate that combined BMAL1, CLOCK, CRY1, PER1, PER2, and melatonin could be potential biomarkers for evaluating the sleep-wake rhythm disturbances of PD patients, which may help improve the accuracy of the clinical assessment and subjective sleep questionnaires. Future studies are needed to further explore our observations using objective measurements of RBD, EDS, and sleep quality.

Although circadian rhythm is endogenous, it is also adapted to the local environmental factors called “zeitgeber,” and light is generally regarded as the most important zeitgeber (Aumann et al., 2016). This study demonstrated a loss of rhythmic expressions of BMAL1, PER1, and melatonin levels in PD with the intensity of environmental light changes. To our knowledge, only one study has analyzed the effect of photoperiod on the tyrosine hydroxylase and dopamine transporter immunoreactivities in the midbrain of people who died in summer vs. winter (Aumann et al., 2016). The density of tyrosine hydroxylase and dopamine transporter immunostaining positive neurons was higher in summer compared with winter, indicating the relationships between environmental photoperiod and the activity of midbrain dopaminergic neurons (Aumann et al., 2016). The reason for the clock gene alterations caused by changing seasonal patterns is still not clear and warrants further investigations. Our results indicate that the greater intensity of the light when sampling, the closer levels of clock genes and melatonin are seen between PD patients and HC. From this point of view, light exposure may have therapeutic effects on PD via regulating melatonin output and clock gene expressions and improving the consolidation of sleep-wake cycles.

Conclusion

In a relatively large clinical sample, this study found that the expression levels of BMAL1, CLOCK, CRY1, PER1, and PER2 in the PBMCs, as well as plasma melatonin levels of PD patients, are significantly decreased compared with HC. Sleep and wakefulness disturbances are common in these PD patients and are associated with all five clock genes and melatonin levels. PD patients with comorbid pRBD had a more severe disturbance of molecular circadian rhythm. We also discovered that a combination of BMAL1, CLOCK, CRY1, PER1, PER2, and melatonin could be a potential biomarker for PD comorbid RBD and EDS. Further, our findings indicate the necessity of strengthening and stabilize circadian rhythmicity through therapeutic intervention on clock genes and melatonin to improve the life quality of PD patients.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Committee of The First Affiliated Hospital of Dalian Medical University (approval number: LCKY2014-29). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

TL and WL designed the project of this manuscript. TL, CC, YLi, NW, YY, and MA-N carried out all the experiments. TL, CC, CJ, and HY contributed to statistical analyses and results interpretation. TL, CJ, and WL contributed to the drafting of the manuscript. TL, YLe, JQ, and WL revised the manuscript. WL contributed to the research concept, research administration. All authors edited and approved the final version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (NSFC 81771521 and 82001483), the National Key R&D Program (2016YFC1306600), and the Guangdong Provincial Key R&D Program (2018B030337001).

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 the patients and controls for participation and nurses in the Department of Neurology of The First Affiliated Hospital, Dalian Medical University for their caring for the patients and contribution of blood sampling.

Supplementary Material

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

Footnotes

1. ^ https://www.timeanddate.com/sun/china/luda

References

Allebrandt, K. V., Teder-Laving, M., and Akyol, M. (2010). CLOCK gene variants associate with sleep duration in two independent populations. Biol. Psychiatry 67, 1040–1047. doi: 10.1016/j.biopsych.2009.12.026

Aumann, T. D., Raabus, M., Tomas, D., Prijanto, A., Churilov, L., Spitzer, N. C., et al. (2016). Differences in number of midbrain dopamine neurons associated with summer and winter photoperiods in humans. PLoS One 11:e0158847. doi: 10.1371/journal.pone.0158847

Bordet, R., Devos, D., Brique, S., Touitou, Y., Guieu, J. D., Libersa, C., et al. (2003). Study of circadian melatonin secretion pattern at different stages of Parkinson’s disease. Clin. Neuropharmacol. 26, 65–72. doi: 10.1097/00002826-200303000-00005

Breen, D. P., Vuono, R., Nawarathna, U., Fisher, K., Shneerson, J. M., Reddy, A. B., et al. (2014). Sleep and circadian rhythm regulation in early Parkinson disease. JAMA Neurol. 71, 589–595. doi: 10.1001/jamaneurol.2014.65

Cai, Y., Liu, S., Sothern, R., Xu, S., and Chan, P. (2010). Expression of clock genes Per1 and Bmal1 in total leukocytes in health and Parkinson’s disease. Eur. J. Neurol. 17, 550–554. doi: 10.1111/j.1468-1331.2009.02848.x

Carpen, J. D., Schantz, M., and Smits, M. (2006). A silent polymorphism in the PER1 gene associates with extreme diurnal preference in humans. J. Hum. Genet. 51, 1122–1125. doi: 10.1007/s10038-006-0060-y

Cox, K. H., and Takahashi, J. S. (2019). Circadian clock genes and the transcriptional architecture of the clock mechanism. J. Mol. Endocrinol. 63, R93–R102. doi: 10.1530/JME-19-0153

Dauvilliers, Y., Schenck, C., Postuma, R. B., Iranzo, A., Luppi, P. H., Plazzi, G., et al. (2008). REM sleep behavior disorder. Nat. Rev. Dis. Primers 4:19. doi: 10.1038/s41572-018-0016-5

De, F., Bisaglia, M., Zordan, M. A., and Sandrelli, F. (2018). Circadian rhythm abnormalities in Parkinson’s disease from humans to flies and back. Int. J. Mol. Sci. 19:3911. doi: 10.3390/ijms19123911

Edgar, R. S., Green, E. W., Zhao, Y., Ooijen, G., Olmedo, M., Qin, X., et al. (2012). Peroxiredoxins are conserved markers of circadian rhythms. Nature 485, 459–464. doi: 10.1038/nature11088

Fifel, K. (2017). Alterations of the circadian system in Parkinson’s disease patients. Mov. Disord. 5, 682–692. doi: 10.1002/mds.26865

Goetz, C. G., Poewe, W., Rascol, O., Sampaio, C., Stebbins, G. T., Counsell, C., et al. (2004). Movement disorder society task Force report on the Hoehn and Yahr staging scale: status and recommendations. Mov. Disord. 19, 1020–1028. doi: 10.1002/mds.20213

Hergenhan, S., Holtkamp, S., and Scheiermann, C. (2020). Molecular interactions between components of the circadian clock and the immune system. J. Mol. Biol. 432, 3700–3713. doi: 10.1016/j.jmb.2019.12.044

Johansson, A. S., Owe-Larsson, B. R., Hetta, J., and Lundkvist, G. B. (2016). Altered circadian clock gene expression in patients with schizophrenia. Schizophr. Res. 174, 17–23. doi: 10.1016/j.schres.2016.04.029

Kalia, L. V., and Lang, A. E. (2015). Parkinson’s disease. Lancet 386, 896–912. doi: 10.1016/S0140-6736(14)61393-3

Korshunov, K. S., Blakemore, L. J., and Trombley, P. Q. (2017). Dopamine: A modulator of circadian rhythms in the central nervous system. Front. Cell Neurosci. 11:91. doi: 10.3389/fncel.2017.00091

Landgraf, D., and Oster, A. S. (2012). Clock genes and sleep. Pflugers Arch. 463, 3–14. doi: 10.1007/s00424-011-1003-9

Laposky, A., Easton, A., Dugovic, C., Walisser, J., Bradfield, C., and Turek, F. (2005). Deletion of the mammalian circadian clock gene BMAL1/Mop3 alters baseline sleep architecture and the response to sleep deprivation. Sleep 28, 395–409. doi: 10.1093/sleep/28.4.395

Li, S., Jia, C., Li, T., and Le, W. (2021). Hot topics in recent Parkinson’s disease research: where we are and where we should go. Neurosci. Bull. Online Head of Print. doi: 10.1007/s12264-021-00749-x

Li, T., and Le, W. (2019). Biomarkers for Parkinson’s disease: how good are they? Neurosci. Bull. 36, 183–194. doi: 10.1007/s12264-019-00433-1

Liu, X., and Le, W. (2020). Profiling non-motor symptoms in monogenic Parkinson’s disease. Front. Aging Neurosci. 12:591183. doi: 10.3389/fnagi.2020.591183

Liu, Y., Niu, L., Liu, X., Cheng, C., and Le, W. (2021). Recent progress in non-motor features of Parkinson’s disease with a focus on circadian rhythm disregulation. Neurosci. Bull. 38, 1010–1024. doi: 10.1007/s12264-021-00711-x

Maggio, R., Vaglini, F., Rossi, M., Fasciani, I., Pietrantoni, I., Marampon, F., et al. (2019). Parkinson’s disease and light: the bright and the dark sides. Brain Res. Bull. 150, 290–296. doi: 10.1016/j.brainresbull.2019.06.013

Mattis, J., and Sehgal, A. (2016). Circadian rhythms, sleep, and disorders of aging. Trends Endocrinol. Metab. 27, 192–203. doi: 10.1016/j.tem.2016.02.003

Mollayeva, T., Turairajah, P., Burton, K., Mollayeva, S., Shapiro, C. M., and Colantonio, A. (2016). The Pittsburgh sleep quality index as a screening tool for sleep dysfunction in clinical and non-clinical samples: a systematic review and meta-analysis. Sleep Med. Rev. 25, 52–73. doi: 10.1016/j.smrv.2015.01.009

Musiek, E. S., and Holtzman, D. M. (2016). Mechanisms linking circadian clocks, sleep, and neurodegeneration. Science 354, 1004–1008. doi: 10.1126/science.aah4968

Pablo-Fenandez, E., Courtney, R., Warner, T., and Holton, J. (2018). A histologic study of the circadian system in Parkinson’s disease, multiple system atrophy, and progressive supranuclear palsy. JAMA Neurol. 75, 1008–1012. doi: 10.1001/jamaneurol.2018.0640

Postuma, R., Berg, D., Stern, M., Poewe, W., Olanow, C. W., Oertel, W., et al. (2015). MDS clinical diagnostic criteria for Parkinson’s disease. Mov. Disord. 30, 1591–1601. doi: 10.1002/mds.26424

Rijo-Ferreira, F., and Takahashi, J. S. (2019). Genomics of circadian rhythms in health and disease. Genome Med. 11, 1–16. doi: 10.1186/s13073-019-0704-0

Schapira, A. H., Chaudhuri, K. R., and Jenner, P. (2017). Non-motor features of Parkinson disease. Nat. Rev. Neurosci. 18, 435–450. doi: 10.1038/nrn.2017.62

Shen, S., Shen, Y., Xiong, K. P., Chen, J., Mao, C. J., Huang, J. Y., et al. (2014). Validation study of REM sleep behavior disorder questionnaire-Hong Kong (RBDQ-HK) in east China. Sleep Med. 15, 952–958. doi: 10.1016/j.sleep.2014.03.020

Sundaram, S., Hughes, R. L., Peterson, E., Muller-Oehring, E. M., Bronte-Stewart, H. M., Poston, K. L., et al. (2019). Establishing a framework for neuropathological correlates and glymphatic system functioning in Parkinson’s disease. Neurosci. Biobehav. Rev. 103, 305–315. doi: 10.1016/j.neubiorev.2019.05.016

Tandberg, E., Larsen, J. P., and Karlsen, K. (1998). A community-based study of sleep disorders in patients with Parkinson’s disease. Mov. Disord. 13, 895–899. doi: 10.1002/mds.870130606

Toh, K. L., Jones, C. R., He, Y., Eide, E. J., Hinz, W. A., Virshup, D. M., et al. (2001). An hPer2 phosphorylation site mutation in familial advanced sleep phase syndrome. Science 291, 1040–1043. doi: 10.1126/science.1057499

Trott, A. J., Menet, J. S., and Achim, K. (2018). Regulation of circadian clock transcriptional output by CLOCK: BMAL1. PLoS Genet. 14:e1007156. doi: 10.1371/journal.pgen.1007156

Videnovic, A., and Willis, G. L. (2016). Circadian system—a novel diagnostic and therapeutic target in Parkinson’s disease? Mov. Disord. 31, 260–269. doi: 10.1002/mds.26509

Videnovic, A., Noble, C., Reid, K. J., Peng, J., Turek, F. W., Marconi, A., et al. (2014). Circadian melatonin rhythm and excessive daytime sleepiness in Parkinson disease. JAMA Neurol. 71, 463–469. doi: 10.1001/jamaneurol.2013.6239

Wisor, J. P., O’Hara, B. F., Terao, A., Selby, C. P., Kilduff, T. S., Sancar, A., et al. (2002). A role for cryptochromes in sleep regulation. BMC Neurosci. 3:20. doi: 10.1186/1471-2202-3-20

Wong, J., Li, J., Pavlova, M., Chen, S., Wu, A., Wu, S., et al. (2016). Risk factors for probable REM sleep behavior disorder: a community-based study. Neurology 86, 1306–1312. doi: 10.1212/WNL.0000000000002414

Wu, Y. H., and Swaab, D. F. (2005). The human pineal gland and melatonin in aging and Alzheimer’s disease. J. Pineal Res. 38, 145–152. doi: 10.1111/j.1600-079X.2004.00196.x

Wulff, K., Porcheret, K., Cussans, E., and Foster, R. G. (2009). Sleep and circadian rhythm disturbances: multiple genes and multiple phenotypes. Curr. Opin. Genet. Dev. 19, 237–246. doi: 10.1016/j.gde.2009.03.007

Yang, G., Chen, L., Grant, G. R., Paschos, G., Song, W. L., Musiek, E. S., et al. (2018). Timing of expression of the core clock gene Bmal1 influences its effects on aging and survival. Sci. Transl. Med. 8:324ra16. doi: 10.1126/scitranslmed.aad3305

Zhang, F., Niu, L., Liu, X. Y., Liu, Y. F., Li, S., Yu, H., et al. (2020). Rapid eye movement sleep behavior disorder and neurodegenerative diseases: an update. Aging Dis. 11, 315–326. doi: 10.14336/AD.2019.0324

Zhang, S., Dai, M., Wang, X., Jiang, S. H., Hu, L. P., Zhang, X. L., et al. (2020). Signalling entrains the peripheral circadian clock. Cell. Signal. 69:109433. doi: 10.1016/j.cellsig.2019.109433