Abstract
Keywords
Introduction
In humans, sleep-wake activity follows a diurnal pattern known as circadian rest-activity rhythm (RAR). 1 Proper RAR is fundamental to optimal health, 2 whereas RAR disruption is associated with an increased risk of chronic diseases, such as cancer and diabetes.1,3 RAR characteristics, including regularity, fragmentation, amplitude, and fitting to 24-h models, can be objectively quantified by applying various methods to several days of continuous accelerometer recordings. 4
As one of the most common causes of permanent functional disability and mortality in elderly people,5,6 cerebral small vessel disease (CSVD) affects deep perforating vessels in the brain and, thus, can be seen on magnetic resonance imaging (MRI). 7 White matter hyperintensity (WMH), lacunar infarct (LI), and cerebral microbleeds (CMBs) are the most acknowledged MRI features.8,9 Furthermore, since these MRI markers do not always occur separately, several scoring systems were invented to represent the overall burden of CSVD in clinical practice.10,11
Besides common features of CSVD, such as cognitive impairment, gait abnormalities, and emotional changes, recent research has elucidated a potential relationship between CSVD and RAR characteristics as well as sleep quality. 12 At present, the identified RAR feature discovered in CSVD patients is still far from being consistent. For example, fragmentation (reflected by intradaily variability [IV]) RAR might be related to MRI markers for CSVD, such as white matter lesions and cerebral microbleeds. 13 In contrast, Oosterman et al. reported that stability (reflected by interdaily stability [IS]) and amplitude (reflected by relative amplitude [RA]) rather than fragmentation were related to deep WMH. 14 In this study, we sought to identify the RAR pattern of elderly persons with CSVD by analyzing actigraphy data using the extended cosinor model and non-parametric methods and clarify the relationship between RAR and CSVD MRI markers.
Methods
Participants
This cross-sectional observational study was conducted between November 1, 2020, and June 20, 2023. It included 115 consecutive elder CSVD cases at the inpatient of Department of Neurology of Seventh Medical Center of PLA General Hospital (Beijing, China). The current research protocol was approved by the Academic Ethics Committee of the Biological Sciences Division of the PLA General Hospital in Beijing, China. Most of the participants showed complaints about slight dizziness, mild headache, mood disturbances, cognitive disorders, etc., with strong will of MRI scan. Patients with the following conditions were excluded: major tremor, cerebral bleeding episodes, or stroke; additional factors, inducing leukoencephalopathy (caused by immune or genetic factors or demyelination); multisystem disorders, such as polyarteritis nodosa, connective tissue diseases–related nervous system vasculitis, arthritis, infection-associated vasculitis; psychotropic medications; musculoskeletal diseases; contraindications for MRI together with additional neurodegenerative diseases such as Alzheimer's disease (clinical diagnostic criteria < NIA-AA>), Parkinson's disease, and frontotemporal dementia. Age, sex, educational level, body mass index (BMI), and a series of vascular risk factors, including hypertension, diabetes, coronary artery disease, hyperlipidemia, smoking, and alcohol use, were considered possible covariates. Hypertension was defined as systolic blood pressure ≥140 mmHg, diastolic blood pressure ≥90 mmHg and/or use of antihypertensive medications. 15 Diabetes was defined as ongoing or newly initiated therapy with antidiabetic drugs or hemoglobin A1c of ≥6.5%. 16 Hyperlipidemia was defined as one of the following points in fasting venous plasma test: total cholesterol ≥6.2 mmol/L; low density lipoprotein cholesterol ≥4.1 mmol/L; triglycerides ≥2.3 mmol/L. 17 All the vascular risk factors were in binomial form.
MRI measurements
The 3.0T MRI brain scan (Discovery MR750; GE Healthcare, USA) revealed white matter lesions that reflected the severity of CSVD. This study used brain MRI (inter-slice and slice thicknesses of 1.5 and 5 mm, separately) with T1 fluid-attenuated inversion recovery (TR, 1750 ms; TE, 23 ms; T1, 780 ms; FOV, 24 cm) and T2-weighted imaging (TR, 7498 ms; TE, 105 ms; FOV, 24 cm). The evaluators had no knowledge of imaging results.
Simple MRI score
The CSVD burden was calculated according to Olama et al. 18 and the simple CSVD MRI score in the current study ranged from 0 to 3 points. In detail, 1 point was awarded for each of the following: WMH graded above 2 points according to Fazekas scale and the presence of LIs (≥3) and CMBs (present). Simple MRI score was related to cognitive deficits, mortality, and neurologic dysfunction in CSVD patients.19–21
Wrist actigraphy
In accordance with our protocol, all cases were asked to wear the ActiGraph GT3X + device (ActiGraph, Pensacola, USA) on the nondominant wrist all day (except when swimming or bathing) for 4–7 days.
ActiLife software (ActiGraph, Pensacola, FL, USA) was used to download data after each wear period. Data with sufficient wear time were chosen and analyzed.
Actigraphic data capture
This work obtained actigraphic data at the epoch of 60 s using the corresponding software. Due to the triaxial accelerometer nature of the ActiGraph GT3X + device, this work sampled diurnal VM using the formula below:
RAR data analysis
First, we used the extended cosinor model, which is derived from R code implementing the sigmoidally transformed cosine (Figure 1). RAR features were exhibited in the following variables.

Illustration of extended cosinor–modeled RAR parameters from a study participant. Time (x-axis) is set in days. RAR α indicates the relative width of the active period, with higher values corresponding to more narrow active periods. RAR β indicates the steepness of the modeled rhythm, with higher values corresponding to more steep curves. Amplitude and mesor are height parameters, and up-mesor, acrophase, and down-mesor are time parameters, as described in the text.

Scatter graphs demonstrating the relationship between WMH and RAR β (A), LI and RAR α (B), as well as CMBs and IV (C).
RAR α: a parameter indicating the relative width of active to rest periods (higher alpha indicates narrower active relative to resting periods). RAR β: a parameter indicating the steepness of RARs (higher values indicate steeper or “square-like” RARs). Acrophase: a parameter indicating the time when activity peaks in the modeled rhythm. Amplitude: a parameter reflecting the difference in magnitude of activity between active and rest phases. Up-mesor: the time when activity passes up through mesor, approximating the time the participant “gets going” in the morning. Down-mesor: the time when activity passes down through the mesor, approximating the time of day the participant “settles down” for the night pseudo-F: a parameter reflecting how well the modeled rhythm fits the observed data (lower values indicate poorer model fit, which suggests that the rhythm might be erratic and/or variable).
Second, non-parametric properties were determined using ActiGraph enumeration data (vector magnitude data) as follows: (1) IS; (2) IV; (3) time occurrence with associated activity counts for a 10-h period with the highest activity (M10) and 5-h period with the lowest activity (L5); (4) RA. 22
IS provides data on RAR synchronization for environmental stimulation with claimed stability, which is calculated using formula 1 below based on an average 24-h profile:
M10 is the maximum sum for a 10-h continuous activity log, and L5 represents the smallest sum for a 5-h continuous activity log. The ratio of (M10-L5)/(M10 + L5) can be used to calculate RA. 23
Sleep quality
Aside from RAR features, sleep variables such as sleep latency (SL), wake after sleep onset (WASO), sleep efficiency (SE), average duration of awakenings (ADA), and times of awakenings (TA) could be automatically analyzed and exported. Daytime sleepiness was evaluated by the Epworth sleepiness scale (ESS). The ESS consists of 8 self-rated items, each scoring from 0 to 3, which measure a subject's habitual “likelihood of dozing or falling asleep” in common situations of daily living. 24
Statistics
For analyzing the differences in demographic and clinical data between groups, Dunnett for one-way ANOVA and the Bonferroni method for the chi-square test were selected to compare the demographic data. Analysis of Covariance (ANCOVA), adjusting for ages was used to compare the RAR parameters and sleep quality. A linear regression analysis adjusted for age, gender, and educational level was used to find the association between RAR variables and CSVD severity. A p < 0.05 was used to indicate statistical significance. SPSS22.0 was used for statistical analysis.
Results
Patients’ demographic characteristics are shown in Table 1. Among groups, differences in demographic features were not significant except for age (67.47 ± 5.18 years versus 69.20 ± 7.03 years versus 72.93 ± 9.32 years versus 73.09 ± 9.70 years; p = 0.014) and hypertension frequency (15.79% versus 50.00% versus 60.60% versus 55.17%; p = 0.009). For the extended cosinor model, elderly patients with simple MRI score of 2–3 points showed statistically lower amplitude compared with individuals with simple MRI score of 0 points (110.57 ± 106.56 versus 68.20 ± 49.72 versus 52.72 ± 52.85, p = 0.044). For the non-parametric method, elderly people with a simple MRI score of 1–3 points exhibited higher fragmentation (higher IV) compared to those with a simple MRI score of 0 points (0.77 ± 0.37 versus 0.97 ± 0.33 versus 1.05 ± 0.32 versus 1.11 ± 0.33, p = 0.020). It seemed that elder patients with a simple MRI score of 2–3 points showed lower RA compared with individuals with a simple MRI score of 0 points (0.78 ± 0.13 versus 0.70 ± 0.12 versus 0.67 ± 0.13, p = 0.057), whereas the trends were only marginal. We did not find obvious differences regarding other RAR parameters and sleep quality. Further details are given in Table 2.
Clinical and demographic characteristics of participants.
SD: standard deviation; CAD: coronary artery disease; MRI: magnetic resonance imaging.
*p < 0.05 simple MRI score = 0 versus simple MRI score = 3.
p < 0.05 simple MRI score = 0 versus simple MRI score = 2.
p < 0.05 simple MRI score = 0 versus simple MRI score = 1.
RAR parameters and sleep quality of participants.
RAR: rest-activity rhythm; IS: interdaily stability; IV: intradaily variability; RA: relative amplitude; MRI: magnetic resonance imaging; SL: sleep latency; SD: standard deviation; SE: sleep efficiency; TTB: total time in bed; TST: total sleep time; WASO: wake after sleep onset; TA: times of awakenings; ESS: Epworth sleepiness scale.
*p < 0.05 simple MRI score = 0 versus simple MRI score = 3.
p < 0.05 simple MRI score = 0 versus simple MRI score = 2.
p < 0.05 simple MRI score = 0 versus simple MRI score = 1.
Furthermore, the association between RAR parameters and CSVD was explored using linear regression analysis. As given in Table 3, for the extended cosinor model, WMH was significantly related to higher RAR β (standardized β = 0.387, p = 0.002), and LI was significantly related to higher RAR α (standardized β = 0.297, p = 0.045). The trends were still significant after adjusting for age, gender, education, BMI, and vascular risk factors. As shown in Table 4, for the non-parametric method, elderly persons with CMBs were more likely to have higher fragmentation reflected by higher IV (standardized β = 0.359, p = 0.015), and the statistical association was unchanged after controlling for age, gender, education, BMI, and vascular risk factors.
Association between CSVD and extended cosinor model RAR parameters.
RAR: rest-activity rhythm; WMH: white matter hyperintensity; LI: lacunar infarct; MRI: magnetic resonance imaging; CSVD: cerebral small vessel disease.
Model 1 represents the association between CSVD and extended cosinor model RAR parameters without adjustment.
Model 2 represents the association between CSVD and extended cosinor model RAR parameters adjusted for age, gender, education, BMI, and vascular risk factors.
*p < 0.05
Association between CSVD and non-parametric RAR parameters.
IS: interdaily stability; IV: intradaily variability; RA: relative amplitude; RAR: rest-activity rhythm; WMH: white matter hyperintensity; LI: lacunar infarct; MRI: magnetic resonance imaging; CSVD: cerebral small vessel disease.
Model 1 represents the association between CSVD and non-parametric RAR parameters without adjustment.
Model 2 represents the association between CSVD and non-parametric RAR parameters adjusted for age, gender, education, BMI, and vascular risk factors.
*p < 0.05
Moreover, we explored the relationship between sleep quality and CSVD. WMH, LI, and CMBs were not related to SE, TTB, WASO, and ESS. Corresponding data are given in Table 5.
Association between CSVD and sleep quality parameters.
SE: sleep efficiency; TTB: total time in bed; WASO: wake after sleep onset; ESS: Epworth sleepiness scale; CSVD: cerebral small vessel disease; WMH: white matter hyperintensity; LI: lacunar infarct.
Model 1 represents the association between CSVD and sleep quality parameters without adjustment.
Model 2 represents the association between CSVD and sleep quality parameters adjusted for age, gender, education, BMI, and vascular risk factors.
Discussion
We found that elderly adults with CSVD displayed a disturbed RAR pattern. Using the extended cosinor model, the amplitude was the most obvious parameter disturbed in CSVD patients. Additionally, using a non-parametric method, RA, which is similar to amplitude derived from the former model, was the most disturbed parameter captured. All these findings might be caused by the diurnal low physical activity in CSVD patients, which was reported in our previous research. 25
Interestingly, we did not find any markedly affected parameter related to RAR shape and time in elderly patients with CSVD. However, IV, a non-parametric parameter, was marginally disrupted in our study, which was according to the results of Zuurbier et al. 13 Being different from the extended cosinor model, the non-parametric approach is characterized as “without making any assumptions about the shape of RAR and time schedule”. 26 Thus, we proposed that CSVD patients showed a fragmented RAR feature; however, this kind of RAR disturbance found in CSVD patients was not in a characteristic shape or exact time schedule.
Furthermore, different CSVD MRI markers were associated with circadian RAR. WMH was negatively associated with RAR β, which was beyond our proposal. Most previous studies reported that lower RAR β (decreased steepness of the RAR curve) was found in patients with higher depressive symptom severity27 or higher risk of stroke. 28 However, there is opposite evidence. For example, higher RAR β was found in elderly people with more severe perceived physical fatigability compared to healthy controls. 29 This could be explained by the fact that “the demands of caregiving may have been the predominant driver determining the beta of the RAR”. 27 Subjects in our study were collected from a ward of neurology, where most severe CSVD patients received high-level caregiving, which might be a possible reason. Research with a larger sample is still needed in the future to clarify the relationship between WMH and RAR patterns.
By assessing the autopsy of CSVD, Sommers et al. 30 failed to find the association between microinfarct and non-parametric RAR parameters. In another study on MRI markers of CSVD, Zuurbier and colleagues did not report the association between LI and non-parametric RAR parameters either. 13 In the current study, LI was positively and independently associated with RAR α. Thus, elderly people with more LI showed narrower active periods of RAR. According to the results of Zhao et al., 31 LI was associated with excessive daytime sleepiness (reflected by ESS score). Although their study did not include the sleep quality of participants with CSVD, the authors assumed that LI might impair the sleep-wake cycle. However, our results did not show an association between LI and objective sleep quality or ESS score. We assumed that other factors, such as apathy and depression, might mediate higher RAR α found with more LI, because both apathy and depression were commonly encountered symptoms in LI.25,32
In line with previous findings, 13 CMBs were positively associated with IV. Besides, there was a lack of an association between objective sleep quality and CSVD. CMBs, especially those located in the cortex, are related to amyloid angiopathy. 33 Our results could suggest that amyloid metabolism in the brain was a potential cause of fragmentation in CSVD patients. A similar conclusion was drawn in a study of Alzheimer's disease. 3 Considering that CSVD can interrupt frontal-subcortical connections, we agreed with the hypothesis that subclinical brain damage might affect circadian rhythms.
As a cost-effective, objective, and comfortable device, wearable actimetry is popularly used in clinical research and large-scale cohorts such as UK Biobank.34,35 Circadian RAR measurement based on actigraphyic device is thus becoming more and more widely recognized. In recent years, abnormal circadian RAR has been evidenced as a new emerging predictor in aged adult, for prodromal Alzheimer's disease, Parkinson disease, even long-term all-cause mortality.3,34,36 More in detail, low amplitude (derived from extended cosinor model) or relative amplitude (derived from non-parametric model) was reported to be a risk factor for adverse health outcomes. Other parameters, reflecting high fragmentation (IV derived from non-parametric model), delayed acrophase (derived from non-parametric model), as well as low robustness (pseudo-F derived from extended cosinor model), were also found to be contributors for delirium and dementia in elderly. 35 Based on the actigraphyic data, interventions should be done by researchers of chronobiology, chrono-exercise, and chrono-nutrition, to clarify the effect of improving RAR on health and longevity in future.
There are two major strengths in the present work. First, the two most commonly used RAR analyzing methods, the extended cosinor model and the non-parametric method, were used to describe RAR features. Second, a simple CSVD score was selected to discover whether RAR was disrupted more severely in elderly adults with a higher CSVD burden. This work also had some limitations. First, the sample size in this study was small. Therefore, a larger sample size is required. Second, we cannot elucidate the temporal relationship between RAR disruption and CSVD since this was a cross-sectional study. Third, enlarged perivascular space, which is correlated with sleep disorder and was not evaluated in this study, should be added in future research.
In summary, the RAR pattern was changed in elderly persons with CSVD. WMH, LI, and CMBs were independently associated with RAR β, RAR α, and IV, respectively.
Footnotes
Acknowledgments
We thank Prof Stephen Smagula for technological support.
Author contributions
Hongyi Zhao (Formal analysis; Writing—original draft); Haiyang Zhang (Formal analysis; Software); Yu Ding (Data curation); Hong Li (Writing—review & editing); Yonghua Huang (Conceptualization; Funding acquisition; Project administration)
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The present work was funded by Wu Jieping Foundation (grant no.: 320.6750.18456).
Declaration of conflicting interests
Yonghua Huang is an Editorial Board Member of this journal but was not involved in the peer-review process of this article nor had access to any information regarding its peer-review.
The remaining authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Data availability
The data that support the findings of this study are available on request from the corresponding author Yonghua Huang.
