Altered Regional Cerebral Blood Flow Associated with Mood and Sleep in Shift Workers: Cerebral Perfusion Magnetic Resonance Imaging Study.

BACKGROUND AND
PURPOSE: Shift work disrupts the body's circadian rhythms and increases the risk of health problems. Despite evidence of neuropsychological disturbances in shift workers (SW), the brain functional status as measured by brain perfusion in chronic shift work has not been evaluated previously. We investigated the regional cerebral blood flow (rCBF) in SW using perfusion MRI (pMRI) and evaluated the relationships between altered rCBF and sleep, mood, psychometric measures, and quality of life.
METHODS: Fifteen rotational SW and 15 day workers (DW) were enrolled. The participants were all female nurses working at a university-affiliated hospital. During 2 weeks of actigraphy they underwent pMRI scanning and psychometric testing on the last day immediately after working. Demographic characteristics, insomnia, daytime sleepiness, and mood were compared between the groups.
RESULTS: The participants were aged 35.3±2.9 years (mean±SD) and had been performing their current work for more than 2 years. The demographic characteristics did not differ between SW and DW, but the levels of insomnia, anxiety, depression, and hyperactivity-restlessness in psychometric measures were higher in SW than in DW. Cerebral perfusion in SW was significantly decreased in the cuneus, fusiform/parahippocampal gyri, and cerebellum of the right hemisphere, while it was increased in the inferior occipital gyrus of the left hemisphere. Perfusion changes in SW were significantly correlated with depression and insomnia severity. The onset and duration irregularity of sleep among SW were related to insomnia, mood, hyperactivity/ restlessness, and quality of life.
CONCLUSIONS: SW experience considerably more insomnia and mood disturbances than do DW, and this is significantly related to perfusion changes in multiple brain areas.

INTRODUCTION

Shift work disrupts the physiological circadian rhythms and acutely deteriorates sleep and other biological functions, which in turn alters both physical and psychological well-being.1 Chronic shift work can increase the risks of cognitive impairment and depression.2 The circadian system not only governs the sleep-wake cycle but also regulates cognition and affective behavior via multiple mechanisms, including noradrenergic innervation of the cortex, clock-gene expression, dendritic spine density, and serotonin transport.3 It is therefore expected that circadian dysfunction results in multiple neuroanatomical, neurophysiological, molecular and neurochemical disruptions during shift work. Circadianrhythm-disrupted mice show a loss of dendritic length, decreased complexity of neurons in the prelimbic prefrontal cortex, and reduced cognitive flexibility or emotional changes consistent with the changes seen in neural architecture.4 The localized brain damage found in that study suggest that certain areas of the brain are selectively vulnerable to circadian-rhythm disruption and that these regions are responsible for the symptoms seen in shift-work disorders.

Despite experimental evidence that disrupted circadian rhythms can induce nerve damage and cognitive/affective system changes, which suggest alterations of the central nervous system in certain brain regions, there have been few neuroimaging studies of shift-work disorders. The first structural MRI study performed in shift workers (SW) showed a reduced volume of the pontomesencephalic tegmentum and related poor sleep quality compared to day workers (DW).5 That preliminary study suggested a chronic disruption of circadian rhythms or decreased exposure to bright light in SW. However, it had limitations in that the participants had been performing shift work for less than 2 years, and that the depression and sleepiness observed in SW were not explained by regional volume changes of the brainstem.

Functional neuroimaging studies are useful for the early detection of metabolic and hemodynamic dysfunction compared to structural MRI.6 Perfusion MRI (pMRI) can be used to estimate various hemodynamic parameters such as the cerebral blood volume (CBV), cerebral blood flow (CBF), and mean transit time (MTT). pMRI has been applied to investigate pathomechanisms of acute sleep restriction7 and obstructive sleep apnea syndrome.8 Those studies showed that regional CBF (rCBF) changes were associated with the relevant symptoms, even without accompanying changes in the gray-matter volume.78 Arterial spin labeling (ASL) is the most recent noninvasive pMRI technology, but it may require a longer imaging time, underestimate CBF changes, and a show poor signal-to-noise ratio.9 The dynamic susceptibility contrast (DSC) technique utilizes exogenous contrast to measure regional blood flow, but it has low sensitivity and spatial resolution.10

To the best of our knowledge, no previous study has investigated brain perfusion in SW. Our hypothesis was that chronic disruption of the sleep-wake pattern by shift work causes localized neuronal damage and associated changes in regional cerebral perfusion. Observations of the regional blood flow alteration might be more sensitive than measuring changes in the gray-matter volume in brain sites earlier in the condition.

The aims of this study were to identify differences in sleep duration and quality between SW and DW and also rCBF changes between them using DSC MRI, while considering their sleep-wake schedules, cognitive function, and mood scores.

METHODS

Participants

Fifteen permanent day-fixed nurses (working during 08:00–17:00) and 15 fast-rotating 3-shift nurses (performing day work during 07:00–15:30, evening work during 15:00–23:30, and night work during 23:00–07:30) were recruited from a metropolitan hospital. Inclusion criteria were 1) aged 25–65 years, 2) duration of current work of ≥2 years, 3) no psychiatric or serious medical problem, 4) not pregnant or planning pregnancy, and 5) no definite obstructive sleep apnea syndrome [apnea-hypopnea index <5/h, as determined in a home sleep study using the Embletta MPR device (Natus Medical, Pleasanton, CA, USA)]. All participants gave written informed consent, and the study protocol was approved by the Ethics Committee of Samsung Medical Center (IRB No. 2014-12-004).

Procedure and questionnaires

Each subject wore an actigraph and completed a sleep- and work-schedule diary for 2 weeks. The Joggle test and MRI scanning were conducted on the 14th day after finishing work on that day (Fig. 1).

Fig. 1

Diagrammatic representation of work schedule and timeline. The work schedule of DW was constant while that of SW involved rotational three shifts. Participants wore an actigraph for 2 weeks and underwent MRI scanning immediately after working. In the case of SW, MRI scanning was performed on the day after two consecutive night shifts. AW: awakening, D: day work, DW: day workers, E: evening work, N: night work, SW: shift workers.

Considering diurnal fluctuations of the brain at rest over the course of a single day and the different effects of fatigue from workload according to shifts, MRI scanning was performed on the day following two consecutive night shifts in SW and at the same time in DW.

The following self-reported questionnaires were applied: Epworth Sleepiness Scale, with a score ≥11 indicating clinically significant daytime sleepiness11; Horne-Östberg Morning-Eveningness Questionnaire12; Insomnia Severity Index (ISI), with >14 indicating clinical insomnia)13; World Health Organization Quality of Life Scale Abbreviated Version (WHOQOL-BREF)14 and Hospital Anxiety and Depression Scale (HADS), with scores ≥8 being clinically significant15; and the Korean version of Conners' Adult Attention Deficit Hyperactivity Disorder Rating Scales (CAARS-K) for measuring inattention/memory problems, hyperactivity/restlessness, impulsivity/emotional lability, and problems with self-concept.16

Actigraphy

Participants were instructed to continuously wear an Actiwatch Spectrum device (Respironics, Murrysville, PA, USA) on the nondominant wrist. This device can detect accelerations down to 0.025 g with a bandwidth between 0.35 Hz and 7.5 Hz and a sampling frequency of 32 Hz. Wrist activity data were analyzed using Actiware (version 5.70; Philips Respironics, Murrysville, PA, USA) sleep software, with the following sleep parameters derived: total sleep time (TST), sleep onset latency (SL), wakefulness after sleep onset (WASO), and sleep efficiency (SE). Variations in the sleep onset time and duration as recorded by the actigraph are regarded as indicators of irregular sleep patterns.17 The irregularity in the sleep onset time was assessed by measuring the standard deviation from the daily sleep onset time and converting this into the accumulated minutes from 0 h on the 24-h system.

Joggle test

Cognitive testing was conducted using automated psychometric measures that are commercially available from Joggle Research (https://admin.joggleresearch.com). The selected battery of tests was as follows: motor praxis task for measuring sensorimotor speed, visual object learning task for visual learning and spatial working memory, fractal two-back test for working memory, abstract matching task for abstraction, line orientation task for spatial orientation, digit symbol substitution task for complex scanning and visual tracking, balloon analog risk task for risk decision-making, and the psychomotor vigilance test for vigilant attention.18

Perfusion MRI

DSC MRI

All MRI images were acquired using a 3.0-T MRI scanner (Intera Achieva 3T, Philips Medical Systems, Best, the Netherlands) with an eight-channel sensitivity-encoding head coil. High-spatial-resolution axial T2*-weighted echo planar images were acquired using the following sequence parameters: repetition time (TR)=1,720 ms, echo time (TE)=35 ms, field of view (FOV)=240×240 mm, matrix=256×256, slices= 50, slice thickness=5.0 mm, gap=1.5 mm, number of dynamic scans=50, and temporal resolution=1.8 s. Baseline images were acquired for 10 dynamic scans, followed by the automatic injection of gadolinium-based contrast agent (0.1 mmoL/kg; Gadovist, Bayer HealthCare Pharmaceuticals, Berlin, Germany) at a rate of 3 mL/s through an 18G intravenous catheter, followed by a 40-mL normal saline flush with additional data acquisition. High-spatial-resolution sagittal T1-weighted images were acquired using the following scan parameters before injecting the contrast agent: TR=10 ms, TE=4.6 ms, FOV=240×240 mm, matrix=480×480, slices=360, slice thickness=0.5 mm.

Quantification of DSC MRI

The flowchart for the quantification of DSC-MRI data is shown in Fig. 2. The MRI signals did not exhibit a linear relationship with the concentration of the contrast agent, and thus the curves of the MRI signal versus time were converted into curves of contrast-agent concentration versus time. The concentration of contrast agent in a given tissue voxel was estimated as C(t)= , where C(t) is the concentration of contrast agent, S(t) is the signal intensity at time t, and S is the baseline signal intensity (Fig. 2A and B). For quantification of DSC-MRI data, the concentration of contrast agent in tissue [C(t] can be written as C(t=CBF(t)R(t−τ)dτ=C(t)⊗[CBF×R(t)], where C(t) is the concentration of contrast agent in a large vessel. In model-independent techniques, CBF×R(t) can be estimated using a deconvolution method. CBF×R(t) was computed using singular value decomposition (SVD) with block circulant matrices.19 SVD was implemented using the svd function in MATLAB (The Mathworks, Natick, MA, USA). The threshold for the diagonal matrix generated using SVD was set to 0.15 of the maximum eigenvalue in order to reduce oscillations of the derived tissue residue function R(t).20 The arterial input function [C(t)] was manually measured by visually inspecting the concentration–time curves near the M1 segment of the middle cerebral artery (Fig. 2C). rCBF was computed as the maximum of the deconvolved tissue impulse response (CBF×R(t) at t=0) using the property R(0)=1, and the regional CBV was calculated as the ratio of the area under the concentration of contrast agent in a given voxel to the area under the curve of the contrast-agent concentration versus time of the arterial input function (Fig. 2D). MTT corresponded to the ratio of CBV to CBF. To obtain the absolute values of CBF and CBV, their regional values were scaled by assuming a normal white-matter CBF of 22 mL/100 g/min and a CBV of 2 mL/100 g.19 For spatial registration, the DSC-MRI image before injecting contrast agent was registered to the T1-weighted image using the affine transformation (from DSC space to T1 space), and then the T1-weighted image was registered to the Montreal Neurological Institute (MNI) space (the MNI152 template in AFNI software) using the affine transformation (from T1 space to MNI space) for each subject. The affine transformation from DSC space to MNI space was calculated using both the affine transformation matrices from DSC space to T1 space and from T1 space to MNI space. Using the obtained affine transformation from DSC space to MNI space, the computed parametric maps were registered to MNI space for group analysis. We performed a voxelby-voxel comparison of perfusion maps between two groups (3dttest++ function in AFNI), and the statistically significant regions were selected for the further analysis of the correlations of clinical scores.

Fig. 2

Flowchart for quantifying the DSC-MRI data. A: DSC-MRI data. B: Concentration of contrast agent. DSC-MRI data were converted into the concentration of contrast agent. C: Concentration–time curves in a large vessel (red box in B) and tissue. The arterial input function was manually determined in large voxels. D: The tissue impulse response function [R(t)] multiplied by the CBF was computed using deconvolution. The regional CBF was computed as the maximum of the deconvolved tissue impulse response, the regional CBV was computed as the ratio of the AUC in a tissue to the AUC of the arterial input function, and the MTT was computed as the ratio of CBV to CBF. E: Observed C(t) (blue dots) and fitted curve using the estimated tissue impulse response function and CBF (red line). AUC: area under the curve, CBF: cerebral blood flow, CBV: cerebral blood volume, Cp: concentration of contrast agent in a large vessel, Ct: concentration of contrast agent in tissue, DSC: dynamic susceptibility contrast, MTT: mean transit time, TR: repetition time.

Statistical analysis

Statistical analyses were performed with SPSS software (version 18.0, SPSS Inc., Chicago, IL, USA). A p value of less than 0.05 was considered to indicate statistical significance. All continuous variables were analyzed using the independent t-test or Mann-Whitney U test, and categorical variables were analyzed using Pearson's chi-square tests. Spearman partial correlations between sleep irregularity and the questionnaire variables in each group were assessed after adjusting for age and the number of years performing the current work. A two-sample t-test was performed to detect significant differences in the mean kurtosis of rCBF between groups. The statistical threshold was set at p<0.001, with a spatial-extent threshold of a cluster size of 50 mm3 (about eight voxels in the original space). For each statistically significant cluster, the averaged values across voxels within each cluster were used in the correlation analysis.

RESULTS

Demographic and clinical characteristics

The demographic characteristics and the findings from the self-reported questionnaires are summarized in Table 1. SW reported a much higher level of insomnia than did DW, with 40% of the SW suffering from clinically significant insomnia. The depression and anxiety levels were also significantly higher in SW, with 80% of them showing a clinically significant depressive mood. In terms of psychometric measures, hyperactivity and restlessness scores were also higher in SW. No significant differences between SW and DW were found in most of the cognitive functions measured by the Joggle test (Supplementary Table 1 in the online-only Data Supplement).

Table 1

Demographic and clinical characteristics

Data are mean±standard-deviation, median (interquartile range), or n (%) values.

*p<0.05, t-test, †Mann-Whitney U test.

BMI: body mass index, CAARS-K: Korean version of Conners' Adult Attention Deficit Hyperactivity Disorder Rating Scales, DW: day workers, ESS: Epworth Sleepiness Scale, HADS: Hospital Anxiety and Depression Scale, ISI: Insomnia Severity Index, MEQ: Horne-Östberg Morning-Eveningness Questionnaire, SW: shift workers, WHOQOL-BREF: World Health Organization Quality of Life Scale Abbreviated Version.

Actigraph-derived sleep parameters

The sleep-wake patterns differed between SW and DW according to their work schedule (Fig. 3). The sleep parameters (TST, SL, WASO, and SE) did not differ between the groups for 2 weeks prior to and on the day before pMRI scanning (Table 2).

Fig. 3

Plots of actigraphy data. Light blue boxes indicate sleeping periods, dark blue boxes indicate where actigraphy recording was started, and arrows indicate working periods.

Table 2

Actigraph-derived sleep parameters

Data are mean±standard-deviation values.

*p<0.05, t-test.

DW: day workers, pMRI: perfusion MRI, SE: sleep efficiency, SL: sleep onset latency, SW: shift workers, WASO: wakefulness after sleep onset.

SW showed a significantly irregular sleep duration [variability in daily sleep duration: 2.07±0.40 vs. 1.37±0.50 h (mean±standard deviation), p<0.001] and sleep onset time (variability in daily sleep onset time: 4.5±0.7 vs. 1.36±0.70 h, p<0.001) compared to DW (Table 2).

The sleep irregularity of SW was positively correlated with scores for hyperactivity and restlessness in the psychometric test, insomnia, and HADS. The WHOQOL-BREF score was negatively correlated with sleep duration and onset irregularity. Sleep onset irregularity was positively correlated with inattention and memory problems in the psychometric test (Table 3).

Table 3

Sleep irregularity and its relationships with clinical factors

Data are r (p) values.

*p<0.05.

CAARS-K: Korean version of Conners' Adult Attention Deficit Hyperactivity Disorder Rating Scales, HADS: Hospital Anxiety and Depression Scale, ISI: Insomnia Severity Index, r: spearman partial correlation coefficient, WHOQOL-BREF: World Health Organization Quality of Life Scale Abbreviated Version.

pMRI findings

rCBF in the cuneus, fusiform/parahippocampal gyri, and cerebellum in the right hemisphere was significantly lower in SW than in DW, while that in the left inferior occipital gyrus was increased in SW (Fig. 4, Table 4).

Fig. 4

Comparison of regional perfusion MRI data between SW and DW. Left inferior occipital gyrus (A), right cerebellum (B), right fusiform/parahippocampal gyri (C), and right cuneus (D). Red color indicates areas with lower CBF in SW compared to DW, and blue color indicates areas with higher CBF. Statistical images are displayed with a cluster probability threshold of p<0.05, corrected for multiple comparisons. CBF: cerebral blood flow, DW: day workers, SW: shift workers.

Table 4

Clusters of brain regions showing cerebral perfusion changes

DW: day workers, MNI: Montreal Neurological Institute, rCBF: regional cerebral blood flow, SW: shift workers.

rCBF in the right fusiform gyrus cluster was negatively correlated with the HADS depression score (r=−0.496, p<0.01) and ISI (r=−0.498, p<0.01). rCBF in the right cerebellum showed a negative correlation with the HADS depression score (r=−0.466, p<0.01), while rCBF in the left inferior occipital gyrus showed a positive correlation with the HADS depression score (r=0.472, p<0.01). The CAARS-K score and actigraph-derived sleep parameters were not correlated with rCBF alterations.

DISCUSSION

The SW in this study were more depressed and anxious as well as more hyperactive and restless than the DW. pMRI showed that rCBF of SW was considerably altered in temporo-occipital and cerebellar areas compared to DW, and those areas were significantly correlated with depression and insomnia severity. rCBF was significantly decreased in the right fusiform gyrus, cuneus, and cerebellum, while the perfusion was somewhat increased in the left inferior occipital area. rCBF in the right fusiform gyrus and cerebellum was negatively correlated with depression severity, while it was positively correlated with the perfusion in the left inferior occipital lobe. Patients with depression show structural or functional derangement in visual regions such as the fusiform/lingual gyri and middle occipital cortex,2122 and derangement of the fusiform gyrus in memory processing is responsible for cognitive vulnerability and depression.23

Functional MRI studies have revealed that the cerebellum is involved in cognition and emotion in addition to motor control.242526 Our findings are in line with previous research showing that emotional memory and the severity of depressive symptoms are associated with structural changes in anterior gray-matter regions of the cerebellum in patients with depression.26

A particularly interesting finding of the present study was the increased rCBF in the left inferior occipital gyrus in SW and the perfusion in this region being positively correlated with depression scores. A previous ASL MRI study in patients with depression found increased perfusion in the temporo-occipital area.27 Another paper reported that hyperactivity in that area might be associated with experiencing sad events or the recognition of unpleasant stimuli.28 The causal relationship between perfusion change and mood disturbance in SW is not clear. All of the present participants denied any mood problems before starting work. Together with findings from previous studies, the present results suggest that perfusion changes in the temporo-occipito-cerebellar regions are associated with the occurrence of mood disturbances in SW.

rCBF in the right fusiform gyrus of SW was negatively correlated with insomnia severity. Previous functional neuroimaging studies have produced mixed findings. The fusiform gyrus is frequently affected in insomnia, exhibiting both increased and decreased activity.2930 One study found that the activity was dependent on the sleep-wake state, with the activity of the fusiform area being increased during sleep but decreased while awake.31 The fusiform gyrus is involved in multimodal combined and facial identification and has been linked to depression.32 These findings together with previous ones may indicate that dysfunction in the region contributes to affective difficulties with insomnia in SW.

The SW in this study exhibited irregular sleep onset and sleep duration, with these symptoms being significantly related to a higher level of depression and insomnia as well as a lower quality of life. Variations in the sleep onset time as recorded by the actigraph was used as an indicator of irregular sleep patterns, and the standard deviation of sleep onset for each participant was regarded as an index of irregularity of the sleep-wake pattern.33 The sleep irregularity of SW was highly significant for the onset time (4.5 h in SW vs. 1.36 h in DW), which may indicate that keeping a regular sleep-wake pattern is essential to prevent sleep and psychological health problems.

This study had several strengths that should be considered. The experience of shift work was long enough (4–12 years) to reflect the characteristics of SW. Most prior studies have relied on self-reported questionnaires to assess the sleep behaviors of SW, whereas we used actigraphy to obtain objective information about sleep and behavior patterns and extensively investigated mood, psychometric features, and cognitive performance. Polysomnography was not conducted in this study because the different sleep schedules of SW may lead to different sleep structures or arousal indices according to particular shifts. Instead, participants were evaluated using home sleep tests, which revealed no evidence of sleep-disordered breathing in any of them. To the best of our knowledge, there has been no previous cerebral perfusion study using DSC-MRI data in SW. Our analysis of DSC-MRI data revealed specific brain regions showing alteration of cerebral perfusion in SW compared to DW, and these brain regions showed significant correlations with depressive mood.

Several limitations of this study should be addressed. The participants were all females in their mid-30 s, and their cognitive changes induced by circadian disruption may differ from those in males.34 The Joggle test in SW indicated that their cognitive function was not disturbed, which is probably related to a lower level of difficulty in test items considering the high innate capacities of younger participants. Perfusion differences between SW and DW might have been attenuated by the relative youngness of the participants and the smallness of the sample. In this study, DW consumed more alcohol than SW. Chronic alcohol consumption exceeding 15 drinks/week is known to reduce the global CBF.35 Although the present DW consumed significantly less alcohol than those in the previous study, it is still possible that the CBF was reduced due to alcohol consumption, which would have underestimated the CBF difference between the two groups. The present MRI studies were conducted on the day following two consecutive night shifts in SW and at the same time in DW, and the schedules of the SW might have differed before their night shifts. Also, the time of day is important when interpreting functional connectivity or rCBF, and circadian and homeostatic influences are important when understanding the results of resting-state data. However, in shift-work studies, the effects of fatigue and workload according to different shifts might overwhelm the time-of-day effects. A previous investigation found a significant correlation between rCBF and the cortisol level in the morning in the anterior cingulate cortex.36 In the present study, the regional perfusion in the anterior cingulate cortex did not differ between SW and DW. This might reduce the need to consider the time of day, but the different MRI study times impair the ability to interpret our results. There could be potential errors in computer perfusion parameters due to susceptibility artifacts of DSC-MRI data, registration errors using an affine transformation, and image distortion due to the phase-encoding polarity.

This study has shown that disrupted circadian rhythm is associated with multiregional cerebral dysfunction; however, the underlying mechanisms remain unknown and so require further investigation.

In summary, we found that female SW were more depressed, anxious, and hyperactive-restlessness than DW, and that the severity of these features was significantly correlated with the irregularity of sleep. Mood and insomnia symptoms were substantially correlated with perfusion changes in SW. This study suggests that SW have altered functional changes in their brains that may partly explain the impact of desynchronization of the circadian rhythm on the development of affective disorder in SW.