Subjects
A total of 173 subjects (93 patients with CI and 80 HCs) diagnosed by the Department of Psychology and Psychiatry of Guangdong Second Provincial General Hospital or recruited from the local community between 2021 and 2023 were included. All subjects were right-handed and had not taken any psychotropic drugs for at least 2 weeks before and during the study to eliminate the effects of the drugs. This study was approved by the Ethics Committee of Guangdong Second Provincial General Hospital. Informed consent was obtained from all participants.
All patients with CI were assessed by two senior psychiatrists based on criteria from the Diagnostic and Statistical Manual of Mental Disorders, version 5. The inclusion criteria for patients with CI were as follows: (a) complaint of difficulty falling asleep, sleep maintenance, and/or early awakening; (b) insomnia occurred at least three times a week for at least 3 months; (c) related daytime dysfunction (e.g., fatigue, emotional disorder, or cognitive impairment); (d) Insomnia Severity Index (ISI) score at least 7 points; and (e) no other sleep disorders (e.g. obstructive sleep apnea or sleep-related movement disorders), severe organic diseases, mental disorders, substance use, or head trauma.
We also recruited 80 age, education, and handedness-matched HCs. The inclusion criteria for the HC group were as follows: (a) good sleep quality; (b) no history of head trauma and no brain lesions, confirmed by conventional T1- or T2-weighted fluid-attenuated inversion-recovery (FLAIR) magnetic resonance imaging; (c) no neuropsychiatric diseases; (d) no pregnant, lactating, or menstruating if female.
Clinical assessments
All patients underwent a series of neuropsychological assessments, including sleep, mood, and cognition-related scales. These scales included the Pittsburgh Sleep Quality Index (PSQI), Insomnia Severity Index (ISI), Self-Rating Depression Scale (SDS), Self-Rating Anxiety Scale (SAS), Montreal Cognitive Assessment (MoCA; total score, visual-spatial and executive functions, naming, memory, attention, language, abstraction, delay recall, and orientation), Digit Symbol Substitution Test (DSST), and Digit Span Test (DST). HCs completed cognitive-related scales consisting of MoCA, DSST, and DST.
MRI data acquisition
All subjects were scanned using a Philips 3.0T MR imaging system (Ingenia, Philips, The Netherlands) at the Guangdong Second Provincial General Hospital. All subjects were placed in an advanced supine position, wearing earplugs to minimize noise, and secured with tight straps and foam pads to reduce head movement. They were banned from taking caffeine, alcohol, or other psychotropic drugs for 48 hours before the scan and instructed to stay awake and try not to think about anything during the scan. Conventional sequences (T1-weighted images and T2-FLAIR images) were obtained to screen for whether lesion in the brain. The rs-fMRI data were performed using a gradient-echo planar imaging sequence with the following parameters: repetition time (TR) = 2000 ms, echo time (TE) = 30 ms, matrix = 64 × 61, field-of-view (FOV) = 224 mm × 224 mm2, flip angle (FA) = 90◦, slice thickness = 3.5 mm with a 1.0 mm gap, interleaved scanning, and 240 volumes in total. Next, high-resolution anatomical images were acquired using a T1-weighted 3D gradient-echo sequence with the following parameters: 185 axial slices, TR = 7.9 ms, TE = 3.6 ms, FA = 8◦, slice thickness = 1.0 mm, no gap, matrix = 256 × 256, and FOV = 256 mm × 256 mm.
Data preprocessing
All images were preprocessed in the Data Processing Assistant for Resting-State fMRI (DPARSF). The preprocessing steps were as follows: (a) remove the first ten volumes to ensure signal stabilization; (b) slice timing; (c) realignment: subjects whose head motion > 2 mm or rotation > 2°in any direction were excluded. Six patients with CI and two healthy controls were excluded due to excessive head motion. (d) spatial normalization to Montreal Neurological Institute space with an isotropic voxel size 3 mm × 3 mm × 3 mm; (e) nuisance covariate regression (Friston’s 24 head motion parameters, white matter signal, and cerebrospinal fluid signal); and (f) Temporal band-pass (0.01–0.10 Hz) filtering.
Static ReHo and dynamic ReHo Calculation
The static ReHo was estimated by computing Kendall’s coefficient of concordance of the times series for each of the 27 nearest neighbor voxels. Then, the images were Z-transformed and spatially smoothed with a Gaussian kernel of 6 mm full-width half-maximum for statistical analysis.
The dynamic ReHo was obtained by using the sliding window method based on the Temporal Dynamic Analysis (TDA) toolkit. The window length and step size are the key factors that affect the algorithm. If the window length is short, the risk of spurious fluctuations increases. The window length is too long to capture its dynamics. In this study, a moderate window length of 30 TR and a shift step size of 1 TR were used and a total of 201 ReHo brain maps were obtained for each participant by calculating the ReHo of the whole brain voxels under each time window. The temporal variability of dynamic ReHo was defined as the standard deviation (SD) of the window-based ReHo maps.
Statistical analysis
The demographic and clinical data between the two groups were compared using the two-sample t-tests or Mann-Whitney U test and chi-square test with SPSS (version 25.0; IBM, Armonk, NY). The threshold for statistical significance was set at p < 0.05, and all hypothesis tests were two-tailed.
To examine the differences between the CI group and HC group in static and dynamic ReHo, a two-sample t-test was employed, with age, education, gender, and mean framewise displacement (FD) as covariates. Multiple comparison correction was performed based on the Gaussian random field theory (GRF, voxel-wise p < 0.005, cluster-wise p < 0.05, two-tailed).
To explore the relationship between altered brain regions and clinical scales in CI, we used the Kolmogorov-Smirnov test to verify whether the data were normally distributed and then selected Spearman correlation or Pearson correlation analysis. P < 0.05 was considered statistically significant.
Finally, binary logistic regression was used to calculate the prediction probabilities of altered static ReHo and dynamic ReHo brain regions, separately and in combination, for diagnosing CI. The receiver operating characteristic (ROC) curves were constructed using the prediction probabilities. The area under the ROC curve (AUC) was calculated to evaluate the diagnostic power of these indicators.