Zijing Deng1,2†
Xiaowei Jiang1,3†
Wen Liu1,2Wenhui Zhao1,2Linna Jia1,2Qikun Sun4Yu Xie1Yifang Zhou1,2Ting Sun1,2Feng Wu1,2
Lingtao Kong1,2
Yanqing Tang2,5*- 1Brain Function Research Section, The First Affiliated Hospital of China Medical University, Shenyang, China
- 2Department of Psychiatry, The First Affiliated Hospital of China Medical University, Shenyang, China
- 3Department of Radiology, The First Affiliated Hospital of China Medical University, Shenyang, China
- 4Department of Radiation Oncology, The First Affiliated Hospital of China Medical University, Shenyang, China
- 5Department of Gerontology, The First Affiliated Hospital of China Medical University, Shenyang, China
Background: Insomnia is considered one of the manifestations of sleep disorders, and its intensity is linked to the treatment effect or suicidal thoughts. Major depressive disorder (MDD) is classified into various subtypes due to heterogeneous symptoms. Melancholic MDD has been considered one of the most common subtypes with special sleep features. However, the brain functional mechanisms in melancholic MDD with insomnia remain unclear.
Materials and methods: Melancholic MDD and healthy controls (HCs, n = 46) were recruited for the study. Patients were divided into patients with melancholic MDD with low insomnia (mMDD-LI, n = 23) and patients with melancholic MDD with high insomnia (mMDD-HI, n = 30), according to the sleep disturbance subscale of the 17-item Hamilton Depression Rating Scale. The dynamic amplitude of low-frequency fluctuation was employed to investigate the alterations of brain activity among the three groups. Then, the correlations between abnormal dALFF values of brain regions and the severity of symptoms were investigated.
Results: Lower dALFF values were found in the mMDD-HI group in the right middle temporal gyrus (MTG)/superior temporal gyrus (STG) than in the mMDD-LI (p = 0.014) and HC groups (p < 0.001). Melancholic MDD groups showed decreased dALFF values than HC in the right middle occipital gyri (MOG)/superior occipital gyri (SOG), the right cuneus, the bilateral lingual gyrus, and the bilateral calcarine (p < 0.05). Lower dALFF values than HC in the left MOG/SOG and the left cuneus in melancholic MDD groups were found, but no significant difference was found between the mMDD-LI group and HC group (p = 0.079). Positive correlations between the dALFF values in the right MTG/STG and HAMD-SD scores (the sleep disturbance subscale of the HAMD-17) in the mMDD-HI group (r = 0.41, p = 0.042) were found. In the pooled melancholic MDD, the dALFF values in the right MOG/SOG and the right cuneus (r = 0.338, p = 0.019), the left MOG/SOG and the left cuneus (r = 0.299, p = 0.039), and the bilateral lingual gyrus and the bilateral calcarine (r = 0.288, p = 0.047) were positively correlated with adjusted HAMD scores.
Conclusion: The occipital cortex may be related to depressive symptoms in melancholic MDD. Importantly, the right MTG/STG may play a critical role in patients with melancholic MDD with more severe insomnia.
Introduction
Major depressive disorder (MDD) is a widely distributed, heterogeneous, and disabling psychiatric disease (1, 2), characterized by different clinical manifestations, such as persistently depressed mood, anhedonia, low self-esteem, somatization, weight change, cognitive dysfunction, retardation, and sleep disturbances (1). Sleep disturbances, particularly insomnia, are prevalent and prodromal clinical characteristics in MDD (3), and 92% of patients with major depressive episodes reported substantial sleep complaints (4). Insomnia, as one of the main manifestations of sleep disturbances, is not only an important risk factor for developing depression (5, 6) but also for recurrence (7), and its severity is associated with the quality of life (8), the severity of depression (8), the therapeutic effect (9, 10), and suicidal thoughts (11, 12). In addition, the link between MDD and insomnia showed a strong bond based on substantial shared genetic liability (13, 14).
Previous research has shown a significant relationship between sleep disturbance (insomnia especially) and depression (15, 16). Furthermore, several neuroimaging studies have also found objective differences in brain regions and changes in neural activity in patients with MDD (17–21) and insomnia individuals (22–25) by using distinctive imaging approaches, such as the amplitude of low-frequency fluctuation (ALFF), functional connectivity (FC), and regional homogeneity (ReHo). These studies implied that neuroimaging characteristics can correctly reflect individual sleep disturbances and performance and could be used as diagnostic biomarkers of MDD (26). There is also some research into the underlying brain mechanisms of patients with MDD with sleep disturbances using resting-state functional MRI (rs-fMRI). Carried out using various methods of rs-fMRI, such as ALFF or FC in patients with MDD with insomnia, researchers noticed changes in the activity of some brain regions, such as the salience network, the suprachiasmatic nuclei, and the default mode network (27–29). A study using a multicenter dataset found that the combination of gray matter (GM) density and fractional ALFF can accurately predict individual total sleep disturbance scores of the 17-Item Hamilton Depression Rating Scale (HAMD-17) which includes items 4 (sleep initiation disorder), 5 (sleep maintenance disorder), and 6 (early awakening) of HAMD-17 in patients with MDD (26). According to certain findings, patients with MDD with more severe insomnia had smaller cortical surface areas in several frontoparietal cortical areas (30), and the interaction between depression and insomnia was related to reduced GM volume in the right orbitofrontal cortex (31). Besides, abnormal global FC density in the visual system was considered a biomarker in patients with MDD with insomnia (32). As a result, the underlying biological mechanisms of sleep disruption in MDD are being investigated constantly.
However, these studies mentioned above did not go deep into each subtype of MDD, and the heterogeneity of symptoms may bring inconsistent experimental results. In addition, newer imaging methods can also be used in future research.
As a heterogeneous mental disorder, dividing MDD into subcategories based on diagnostic criteria is essential for identifying underlying pathophysiological mechanisms, and this form of individualized diagnosis is required for precise MDD treatment (26). Melancholic MDD has been viewed as the most serious subtype of MDD (33) with specific clinical symptoms, such as persistent anhedonia, psychomotor disturbances, cognitive impairment, early morning awakening, excessive guilt, and anorexia (34). In addition, there are several biological indicators, such as disturbances in sleep architecture, occurring more frequently in melancholia than in other types of depression (35, 36). Studies also found that patients with melancholic MDD had a higher rate of nightmares, middle, and terminal insomnia than patients with non-melancholic MDD (37, 38). Previous research revealed that patients with melancholic MDD displayed evident sleep-electroencephalogram (EEG) changes, which involved low quantities of slow-wave sleep, disrupted sleep, a short rapid eye movement (REM) latency, and a high REM density compared to non-melancholic depression (36, 39–41). However, limited studies use imaging methods to investigate underlying neuroimaging mechanisms of melancholic MDD with insomnia.
Resting-state functional magnetic resonance imaging (rs-fMRI) is a method to measure spontaneous activity in the brain that has been generally utilized to explore functional changes in the human brain (42). There are also several studies about melancholic MDD. These studies found that the melancholic MDD had distinct fractional ALFF values in the right middle temporal gyrus (MTG)/and bilateral superior occipital gyrus (SOG) (43), lower ReHo values in the right SOG/middle occipital gyrus (MOG) (44), different network homogeneity values in the right posterior cingulate cortex (PCC)/precuneus, right angular gyrus, and the right MTG (45), decreased voxel-mirrored homotopic connectivity values in the fusiform gyrus, PCC, and SOG (46). As a promising approach to rs-fMRI, ALFF has been adopted widely because it can reflect the spontaneous neuronal activity of specific regions (47). Simultaneously, ALFF is also considered to be momentarily stationary during a typical rs-MRI session (48). However, the activity of brain areas is inherently dynamic (49). Therefore, to capture the dynamic properties of distinct brain regions’ activity effectively, the dynamic sliding window method such as dALFF has been adopted in several types of research in contrast to the “static” approaches of rs-fMRI. Furthermore, a previous study demonstrated that dALFF changes are related to changes in EEG band power (50). Recently, many psychiatric illnesses have discovered a significant change in dALFF, such as MDD (51), bipolar disorder (52), Schizophrenia (53), primary insomnia (54), and generalized anxiety disorder (55). For an instance, Staner found that temporal dALFF can predict suicidal thoughts in depressed patients, while static ALFF cannot (16). Nevertheless, the dynamics of brain activity in patients with melancholic MDD with insomnia have not yet been quantified.
Therefore, our goal is to investigate the dynamic spontaneous brain activity in patients with melancholic MDD along with insomnia compared to HCs by employing the dALFF. We assume that melancholic MDD with variable degrees of insomnia may show certain particular dynamic brain functional changes and these features are associated with insomnia symptoms of melancholic MDD.
Materials and methods
Participants
In the current research, a total of 53 melancholic MDD individuals and 46 healthy controls (HCs) matched for age, gender, and education were recruited from the Shenyang Mental Health Center and The First Affiliated Hospital of China Medical University. Participants ranged in age from 18 to 40. The Medical Research Ethical Committee of the First Affiliated Hospital of China Medical University has approved the study. Before starting the study, all participants were informed of the purpose of the study and signed an informed consent form.
The diagnosis of melancholic MDD was independently performed by two seasoned psychiatrists, according to the Diagnostic and Statistical Manual of Mental Disorders, fourth edition (DSM-IV) criteria. The inclusion criteria for melancholic MDD were: (1) continuously anhedonia; (2) at least three of the symptoms as follows are required: depressive mood, feeling worse in the morning, early morning awakening, remarkable psychomotor disturbances, severe anorexia or weight loss, and excessive or inappropriate guilt. Patients were enrolled in the research while experiencing a depressive episode with a total of HAMD-17 scores ≥ 17 on the day of the MRI scan. These patients were excluded if diagnosed with other Axis I or Axis II disorder.
HCs with HAMD-17 scores ≤ 7 were recruited from the local community and they had no Axis I or Axis II disorder. Meanwhile, their first-degree relatives could not have any family history of psychiatric disorders.
The following exclusion criteria applied to all participants: history of neurological such as head injury, stroke, seizures, or transient ischemic attack; caffeine, drug, or alcohol abuse; and pregnancy.
Assessments
The depression severity was evaluated by HAMD-17 (56) and the anxiety severity was measured based on the Hamilton Anxiety Rating Scale (HAMA) (57). Insomnia symptoms were evaluated using the insomnia subscale of the HAMD-17 (HAMD-SD) (22, 26, 58), which involved items 4 (sleep initiation disorder), 5 (sleep maintenance disorder), and 6 (early awakening). The scores of the HAMD insomnia subscale and sleep diary data are widely considered to be well correlated, which has been a global measure of insomnia severity in depressive disorders (59). The adjusted HAMD scores meant the HAMD-17 scores after the omission of the sleep subscale scores (27, 32). Patients with melancholic MDD were separated into patients with melancholic MDD with low insomnia group (mMDD-LI, HAMD-SD ≤ 3, n = 23) and patients with melancholic MDD with high insomnia (mMDD-HI, HAMD-SD ≥ 4, n = 30), according to the HAMD-SD.
Image acquisition
The rs-MRI functional images were acquired by adopting a 3.0 T GE SIGNA MRI system at the Image Institute of The First Affiliated Hospital of China Medical University, Shenyang, China. The following are the parameters that were obtained using a gradient echo-planar imaging (EPI) sequence: repetition time (TR) = 2,000 ms, echo time (TE) = 40 ms, field of view (FOV) = 240 × 240 mm2, flip angle (FA) = 90°, image matrix size = 64 × 64, slices = 35, slice thickness = 3 mm, spacing between slices = 3 mm. Participants were required to wear earplugs and use foam pads to reduce scanner noise and head motion, and their head was fixedly positioned during scan time. The resting state meant that participants were not doing any cognitive tasks during the scan time and they were asked to keep relaxed and awake with their eyes closed while simultaneously keeping their minds blank. After completing scanning, we obtained a total of 200 volumes of images.
Image pre-processing
Images were processed by applying the Statistical Parametric Mapping 12 (SPM12)1 and the Data Processing Assistant for Resting-State fMRI (DPABI 4.1, Advanced edition) (60) based on the custom code written in MATLAB. To maintain the stability of the initial signal, the first 10 volumes of each participant’s scanned data were eliminated. Then, the remaining 190 images were adjusted for slice-timing and head motion. Participants with the translation of more than 2 mm or rotation of more than 2° of head motion in each direction were not included in the present study. The mean framewise displacement (FD) was used to measure the scrubbing-related micro-head motion of each participant. Additionally, no statistical difference existed in mean FD when comparing the three groups. After realignment, the imaging data corrected were spatially normalized into a standard EPI template in the Montreal Neurological Institute (MNI) space and resampled to a voxel size of 3 mm × 3 mm × 3 mm. The spatial smoothing of the EPI images adopted a 4 mm full width at half maximum (FWHM) Gaussian kernel. The BOLD signals were then detrended to correct a linear trend. Finally, the linear regression of the nuisance covariates was performed to remove the effects, including head motion parameters, cerebrospinal fluid signal, and white matter signal. Normalization of T1 images was also conducted, and details and results are shown in Supplementary Table 3 and Supplementary Figure 3, and pre-processing part.
Dynamic amplitude of low-frequency fluctuation acquisition
The dynamic ALFF was computed using temporal dynamic analysis (TDA) toolkits in DPABI 4.1 (60) based on a sliding window analysis. The window length which was considered a crucial parameter of resting-state dynamic computation should be larger than 1/fmin. The fmin referred to the minimum frequency of the time series. As a shorter window length is more likely to introduce misleading fluctuations in the observed dALFF, a longer window length may fail to discover the potential dALFF. Therefore, a window length of 50 TRs (100 s) was selected to compute the temporal variability of ALFF according to prior research (51, 61). We calculated the ALFF of each subject with a sliding window and obtained the ALFF values of each given voxel, the time series of which were then transformed into a frequency domain with a fast Fourier transform. The square root of the power spectrum of each voxel was computed and summed from 0.01 to 0.08 Hz. The dALFF values were obtained through computing the standard deviation (SD) of ALFF values at each voxel across the sliding window dynamics. SD is a useful and quantified indicator that is generally applied to describe the degree of change in dALFF in the research of temporal dynamic brain activity. Therefore, we used SD as dALFF for the next analysis.
Statistical analysis
The demographic features and clinical symptoms were analyzed by several statistical methods, such as one-way analysis of variance (ANOVA), two-sample t-test, and chi-square test, and p < 0.05 was set as significance. The specific use of these statistical methods is shown in Table 1. The relationships between the total HAMD-17 scores, adjusted HAMD scores, HAMD-SD scores, and HAMA scores were explored by Pearson’s correlation in the pooled melancholic MDD groups, and the statistical significance was also set at p < 0.05. Images data were analyzed by two software, and two types of correction are used. Altered dynamic ALFF values across the three groups were investigated by ANOVA in DPABI 4.1 and the results were corrected with the Gaussian random field (GRF) one-tailed (62), the threshold of voxel-wise was p < 0.001, and the cluster-level was p < 0.05. The extracted dALFF values were evaluated by SPSS (Statistical Product and Service Solutions) and the one-way analysis of covariance (ANCOVA) was used to compare differences among three groups when age, gender, years of education, and FD were covariates. Significance was set at p < 0.05 after correcting by Bonferroni correction. The partial correlation analysis was conducted to examine the relationship between the altered dALFF values and clinical features (HAMD-SD scores, HAMA scores, and adjusted HAMD scores) in melancholic MDD groups when age, gender, medication, HAMA scores, adjusted HAMD scores, or HAMD-SD scores were used as covariates. Significance was set at p < 0.05 after correcting by Bonferroni correction.
Table 1. Demographic and clinical measurement among three groups.
Validation analysis
To verify the results of dALFF variability obtained from a sliding window length of 50 TR, we adopted using two other window lengths (30 and 70 TR).
Results
Demographic characteristics and clinical symptoms
As Table 1 illustrates age, sex, years of education, and mean FD were not significantly different among the mMDD-LI, mMDD-HI, and HC groups (p > 0.05). And there were no significant differences in medication, duration of disease, or number of episodes between the two melancholic MDD groups (p > 0.05). Compared with the mMDD-LI, mMDD-HI showed higher HAMD (p = 0.012) and HAMA scores (p = 0.039). However, the adjusted HAMD scores were not statistically different between the two melancholic MDD groups (p = 0.50). In the pooled melancholic MDD, the HAMD-SD scores were significantly positively associated with the HAMD scores (r = 0.452, p = 0.001) and HAMA scores (r = 0.377, p = 0.005), but not with the adjusted HAMD scores (r = 0.178, p = 0.203). The HAMA scores were also positively associated with the HAMD scores (r = 0.68, p < 0.001) and the adjusted HAMD scores (r = 0.629, p < 0.001) in the pooled melancholic MDD.
Differences in dynamic amplitude of low-frequency fluctuation values among three groups
Significant group differences of dALFF were found in the right MTG/superior temporal gyrus (STG), the bilateral MOG/SOG, the bilateral cuneus, the bilateral calcarine, and the bilateral lingual gyrus (p < 0.001, GRF correction). The ANCOVA revealed that the mMDD-HI group displayed lower dALFF values in the right MTG/STG than the mMDD-LI group (p = 0.014, Bonferroni correction) and HC group (p < 0.001, Bonferroni correction). Meanwhile, two melancholic MDD groups showed lower dALFF values than HC in the right MOG/SOG, the right cuneus, the bilateral lingual gyrus, and the bilateral calcarine (p < 0.05, Bonferroni correction). The two melancholic MDD groups had lower dALFF values than HC in the left MOG/SOG and the left cuneus, although the difference between the mMDD-LI group and HC group was not statistically significant (p = 0.079, Bonferroni correction). Table 2 and Figure 1 illustrate these findings.
Table 2. Brain regions showing significant group differences in dALFF values.
Figure 1. Differences in dALFF values and post hoc analysis of dALFF values among three groups. Cluster 1: right MTG/STG, Cluster 2: bilateral lingual gyrus/calcarine; Cluster 3: right middle occipital gyrus/superior occipital gyrus/cuneus; Cluster 4: left middle occipital gyrus/superior occipital gyrus/cuneus. (A) Significant difference in dALFF values among the three groups. Significant at p < 0.001 after GRF correction. R, right side; L, left side. (B) Post hoc analysis of dALFF values with significant variations across the three groups. ***p < 0.001 level, **p < 0.01 level, *p < 0.05 level, Bonferroni correction.
Correlations between dynamic amplitude of low-frequency fluctuation values and clinical characteristics
The partial correlation analysis revealed that the dALFF values in the right MTG/STG were positively correlated with HAMD-SD scores in the mMDD-HI group when age, gender, medication, and adjusted HAMD scores, and HAMA scores were covariates (r = 0.41, p = 0.042). In the pooled melancholic MDD, the dALFF values in the right MOG/SOG and the right cuneus were significantly positively correlated with adjusted HAMD scores (r = 0.338, p = 0.019), the dALFF values in the left MOG/SOG and the left cuneus were positively correlated with adjusted HAMD scores (r = 0.299, p = 0.039), and the dALFF values in the bilateral lingual gyrus and the bilateral calcarine were positively correlated with adjusted HAMD scores (r = 0.288, p = 0.047) when age, gender, medication, HAMD-SD scores, HAMA scores were as covariates. The above results are shown in Tables 3, 4. No other correlation was found.
Table 3. The partial correlation analysis between the dALFF values and adjusted HAMD scores in the pooled melancholic MDD.
Table 4. The partial correlation analysis between the dALFF values and HAMD-SD scores in patients with melancholic MDD with high insomnia.
Validation results
In the present study, we validated our results by using different sliding window lengths (30 and 70 TR). Finally, the relevant results were presented in Supplementary Tables 1, 2 and Supplementary Figures 1, 2.
Discussion
In the present study, we investigated the neuroimaging features of melancholic MDD with varying degrees of insomnia by adopting the dynamic sliding window method. Our research found that the mMDD-HI group displayed significantly decreased dALFF values in the right MTG/STG than the mMDD-LI and HC groups. Furthermore, the dALFF values in the right MTG/STG were significantly positively correlated with HAMD-SD scores in the mMDD-HI group. We also found two melancholic MDD groups showed decreased dALFF values than HC in the right MOG/SOG, the right cuneus, the bilateral lingual gyrus, and the bilateral calcarine. The two melancholic MDD groups had lower dALFF values than HC in the left MOG/SOG and the left cuneus, but the difference between the mMDD-LI group and HC group was not statistically significant. Meanwhile, the dALFF values in the bilateral MOG/SOG, the bilateral cuneus, the bilateral lingual gyrus, and the bilateral calcarine were significantly positively correlated with adjusted HAMD scores in the pooled melancholic MDD. These findings suggested that the bilateral occipital cortex may be associated with the neuroimaging mechanisms of melancholic MDD depressive symptoms other than insomnia, while the temporal cortex may play an important role in more severe insomnia symptoms in patients with melancholic MDD.
In the present study, the mMDD-HI group revealed alterations of dALFF in the right temporal cortex. Meanwhile, the right temporal cortex may be associated with more severe insomnia symptoms in patients with melancholic MDD. Previous studies also speculated that functional alterations in the right MTG/STG may be a neurobiological characteristic of melancholic MDD (44, 45). Some studies found GM hypertrophies (63) and diminished ALFF in the right MTG (64) in patients with insomnia disorder. The altered fALFF brain areas of patients with insomnia were found in the right STG, which was unaffected by age and gender factors (65). This evidence suggests that the right MTG/STG may play an important role in insomnia symptoms of melancholic depression. Previous studies have proposed the hyperarousal hypothesis to explain the mechanism of insomnia, which is the result of hyperarousal in emotional or cognitive networks (66, 67). Furthermore, emotional and cognitive functions, such as emotional regulation, social cognition, and the process of memory, have been considered to be associated with the temporal lobe (55, 68), which may explain the hyperarousal hypothesis. Hence, we speculated that the abnormal dALFF values in the right temporal cortex may be one of the neuroimaging mechanisms of high insomnia symptoms in melancholic MDD.
However, there are also studies that are inconsistent with our results. They found several brain regions, such as the left inferior temporal gyrus (32), the posterior parahippocampal/hippocampal gyrus (32), frontoparietal cortical areas (30), caudate nuclei (69), and the anterior insular subregions (70), were associated with insomnia symptoms in patients with MDD. We went on to think about the reasons for these differences, which may be caused by distinct rs-fMRI methods or the heterogeneous symptoms of patients with MDD. This suggests that we need to expand the sample size, pay more attention to the heterogeneity of MDD, and adopt newer methods for analysis in future.
In the current study, we also discovered that melancholic MDD caused changes in dALFF values in the occipital cortex. Furthermore, these findings revealed that the occipital cortex may be associated with depressive characteristics. In recent years, multiple studies reported that patients with melancholic MDD manifested considerably lower values in the MOG/SOG than patients with non-melancholic MDD or HCs by applying rs-fMRI neuroimaging methods (44–46). However, there is not much literature on melancholic MDD, and we found similar results in other studies on MDD. Another research also found lower ALFF in the bilateral occipital cortex in treatment-free MDD subjects (71). The results indicate that decreased dALFF values of the calcarine in MDD are essentially consistent with findings from another study, where lower ALFF values were investigated in the left calcarine (27). Apart from functional changes, melancholic depression involves structural abnormalities, according to growing neuroimaging research evidence. Previous studies have found that patients with melancholic MDD showed asymmetrical cerebrospinal fluid (CSF) volume enlargement in the Sylvian fissure region (72, 73). The increased CSF volume may cause occipital bending, which may result in functional changes in the occipital cortex (74). These present findings are in line with the results of the previous studies. Therefore, these results revealed that alterations in the occipital cortex may be specific in melancholic MDD.
Besides, correlation analysis results show that the occipital cortex may be associated with depressive traits. The lingual gyrus and cuneus, located in the medial occipital lobe, are necessary for both basic and advanced visual processing (75). The function of the occipital cortex is to process initial visual stimuli, consolidate information into visual working memory (76), and revolve in the consciousness of facial expressions (77), and its aberrant changes may explain working memory impairment in melancholic MDD (78). A study demonstrated functional changes in the occipital cortex when depressed adolescents are presented with sad faces, and the occipital cortex is important to the interaction of depression status, attention, and mood processing (79). In addition, the functional stability in the middle occipital gyri (MOG) was negatively linked with depressive symptoms, which could predict long-term recovery in depressive symptoms in patients with MDD (80). As a result, we hypothesized that the occipital cortex maybe plays a major role in the neuroimaging mechanisms of other depressive symptoms more than insomnia symptoms in melancholic.
In addition, we also observed that the mMDD-HI group showed lower dALFF values in all brain regions found in the current study compared to the mMDD-LI group and HC group, although certain differences did not achieve statistical significance. However, previous similar investigations using different imaging methods have yielded inconsistent results (27, 32). As a result, future studies will require a larger sample size to verify these results.
From our research, we observed several limitations. First, HAMD was adopted to evaluate the severity of insomnia, which may cause an incomplete assessment. In further studies, we could apply more objective and subjective measurements, such as Pittsburgh Insomnia Rating Scale or polysomnography, to improve the precision of the severity of insomnia assessments. Second, most patients had been on medicine before the scan and the drug-related effects cannot be excluded. In subsequent statistical analysis, we attempted to limit the impact of medication on the results by taking medication as a covariance. More medication-free patients should be enrolled in our research. Finally, our study was a cross-sectional study with a small sample size for each subgroup, and the sample size should be expanded in future, and relevant longitudinal work should be conducted.
In conclusion, the present study investigated brain functional alterations in melancholic MDD with insomnia. We found that the occipital cortex may be related to depressive symptoms in those with melancholic MDD, which further supports the findings of previous imaging studies on melancholic MDD. Moreover, the major finding of the present study revealed that the right temporal cortex was associated with more severe insomnia symptoms in melancholic MDD, which is consistent with our earlier hypothesis, and it suggested that the right temporal cortex could be a biomarker for melancholic MDD with high insomnia.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the corresponding author.
Ethics statement
The studies involving human participants were reviewed and approved by the Medical Research Ethical Committee of The First Affiliated Hospital of China Medical University. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author contributions
ZD performed investigation, data processing, statistical analyses, and visualization, and wrote the original draft. XJ was responsible for conceptualization, data curation, writing of original draft, and funding acquisition. WL wrote the original draft. WZ and LJ performed investigation and data processing. QS acquired the data. YX was responsible for statistical analyses. YZ and TS validated the results. FW and LK recruited patients, confirmed the diagnosis, and acquired funding. YT was responsible for conceptualization, project administration, and funding acquisition. All authors reviewed and approved the manuscript.
Funding
This work was supported by the National Key R&D Program of China (Grants #2018YFC1311600 and 2016YFC1306900 to YT), the Liaoning Revitalization Talents Program (Grant #XLYC1808036 to YT), the Natural Science Foundation of Liaoning Province (2020-MS-176 to XJ), the National Key R&D Program “Science and Technology Winter Olympics” (2021YFF0306503 to FW), the Joint Fund of National Natural Science Foundation of China (U1808204 to FW), and the Natural Science Foundation of Liaoning Province (2019-MS-05 to FW).
Acknowledgments
We thank all participants who contributed to this research.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2022.958994/full#supplementary-material
Footnotes
References
1. American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders (DSM-5). Arlington, V A: American Psychiatric Association (2013).
2. Lu J, Xu X, Huang Y, Li T, Ma C, Xu G, et al. Prevalence of depressive disorders and treatment in China: a cross-sectional epidemiological study. Lancet Psychiatry. (2021) 8:981–90. doi: 10.1016/s2215-0366(21)00251-0
3. Perlis ML, Giles DE, Buysse DJ, Tu X, Kupfer DJ. Self-reported sleep disturbance as a prodromal symptom in recurrent depression. J Affect Disord. (1997) 42:209–12. doi: 10.1016/s0165-0327(96)01411-5
4. Geoffroy PA, Hoertel N, Etain B, Bellivier F, Delorme R, Limosin F, et al. Insomnia and hypersomnia in major depressive episode: Prevalence, sociodemographic characteristics and psychiatric comorbidity in a population-based study. J Affect Disord. (2018) 226:132–41. doi: 10.1016/j.jad.2017.09.032
5. Riemann D, Voderholzer U. Primary insomnia: a risk factor to develop depression? J Affect Disord. (2003) 76:255–9. doi: 10.1016/s0165-0327(02)00072-1
6. Baglioni C, Battagliese G, Feige B, Spiegelhalder K, Nissen C, Voderholzer U, et al. Insomnia as a predictor of depression: a meta-analytic evaluation of longitudinal epidemiological studies. J Affect Disord. (2011) 135:10–9. doi: 10.1016/j.jad.2011.01.011
7. Li SX, Lam SP, Chan JW, Yu MW, Wing YK. Residual sleep disturbances in patients remitted from major depressive disorder: a 4-year naturalistic follow-up study. Sleep. (2012) 35:1153–61. doi: 10.5665/sleep.2008
8. McCall WV, Reboussin BA, Cohen W. Subjective measurement of insomnia and quality of life in depressed inpatients. J Sleep Res. (2000) 9:43–8. doi: 10.1046/j.1365-2869.2000.00186.x
9. Jindal RD, Thase ME. Treatment of insomnia associated with clinical depression. Sleep Med Rev. (2004) 8:19–30. doi: 10.1016/s1087-0792(03)00025-x
10. Howland RH. Sleep interventions for the treatment of depression. J Psychos Nurs Mental Health Serv. (2011) 49:17–20. doi: 10.3928/02793695-20101208-01
11. Richardson JD, Thompson A, King L, Corbett B, Shnaider P, St Cyr K, et al. Insomnia, psychiatric disorders and suicidal ideation in a National Representative Sample of active Canadian Forces members. BMC Psychiatry. (2017) 17:211. doi: 10.1186/s12888-017-1372-5
12. Pompili M, Innamorati M, Forte A, Longo L, Mazzetta C, Erbuto D, et al. Insomnia as a predictor of high-lethality suicide attempts. Int J Clin Pract. (2013) 67:1311–6. doi: 10.1111/ijcp.12211
13. Baranova A, Cao H, Zhang F. Shared genetic liability and causal effects between major depressive disorder and insomnia. Hum Mol Genet. (2021) 31:1336–45. doi: 10.1093/hmg/ddab328
14. Melhuish Beaupre LM, Tiwari AK, Gonçalves VF, Zai CC, Marshe VS, Lewis CM, et al. Potential genetic overlap between insomnia and sleep symptoms in major depressive disorder: a polygenic risk score analysis. Front Psychiatry. (2021) 12:734077. doi: 10.3389/fpsyt.2021.734077
15. Riemann D, Krone LB, Wulff K, Nissen C. Sleep, insomnia, and depression. Neuropsychopharmacology. (2020) 45:74–89. doi: 10.1038/s41386-019-0411-y
16. Staner L. Comorbidity of insomnia and depression. Sleep Med Rev. (2010) 14:35–46. doi: 10.1016/j.smrv.2009.09.003
17. Yao X, Yin Z, Liu F, Wei S, Zhou Y, Jiang X, et al. Shared and distinct regional homogeneity changes in bipolar and unipolar depression. Neurosci Lett. (2018) 673:28–32. doi: 10.1016/j.neulet.2018.02.033
18. Sun H, Luo L, Yuan X, Zhang L, He Y, Yao S, et al. Regional homogeneity and functional connectivity patterns in major depressive disorder, cognitive vulnerability to depression and healthy subjects. J Affect Disord. (2018) 235:229–35. doi: 10.1016/j.jad.2018.04.061
19. Wang M, Ju Y, Lu X, Sun J, Dong Q, Liu J, et al. Longitudinal changes of amplitude of low-frequency fluctuations in MDD patients: A 6-month follow-up resting-state functional magnetic resonance imaging study. J Affect Disord. (2020) 276:411–7. doi: 10.1016/j.jad.2020.07.067
20. Kaiser RH, Andrews-Hanna JR, Wager TD, Pizzagalli DA. Large-Scale network dysfunction in major depressive disorder: a meta-analysis of resting-state functional connectivity. JAMA Psychiatry. (2015) 72:603–11. doi: 10.1001/jamapsychiatry.2015.0071
21. Wang L, Kong Q, Li K, Su Y, Zeng Y, Zhang Q, et al. Frequency-dependent changes in amplitude of low-frequency oscillations in depression: a resting-state fMRI study. Neurosci Lett. (2016) 614:105–11. doi: 10.1016/j.neulet.2016.01.012
22. Liu CH, Liu CZ, Zhang J, Yuan Z, Tang LR, Tie CL, et al. Reduced spontaneous neuronal activity in the insular cortex and thalamus in healthy adults with insomnia symptoms. Brain research (2016) 1648(Pt A):317–24. doi: 10.1016/j.brainres.2016.07.024
23. Zhu Y, Zhao X, Yin H, Zhang M. Functional connectivity density abnormalities and anxiety in primary insomnia patients. Brain Imaging Behav. (2021) 15:114–21. doi: 10.1007/s11682-019-00238-w
24. Zhang Y, Zhang Z, Wang Y, Zhu F, Liu X, Chen W, et al. Dysfunctional beliefs and attitudes about sleep are associated with regional homogeneity of left inferior occidental gyrus in primary insomnia patients: a preliminary resting state functional magnetic resonance imaging study. Sleep Med. (2021) 81:188–93. doi: 10.1016/j.sleep.2021.02.039
25. Ran Q, Chen J, Li C, Wen L, Yue F, Shu T, et al. Abnormal amplitude of low-frequency fluctuations associated with rapid-eye movement in chronic primary insomnia patients. Oncotarget. (2017) 8:84877–88. doi: 10.18632/oncotarget.17921
26. Shi Y, Zhang L, He C, Yin Y, Song R, Chen S, et al. Sleep disturbance-related neuroimaging features as potential biomarkers for the diagnosis of major depressive disorder: a multicenter study based on machine learning. J Affect Disord. (2021) 295:148–55. doi: 10.1016/j.jad.2021.08.027
27. Liu CH, Guo J, Lu SL, Tang LR, Fan J, Wang CY, et al. Increased Salience Network Activity in Patients With Insomnia Complaints in Major Depressive Disorder. Front Psychiatry. (2018) 9:93. doi: 10.3389/fpsyt.2018.00093
28. McKinnon AC, Hickie IB, Scott J, Duffy SL, Norrie L, Terpening Z, et al. Current sleep disturbance in older people with a lifetime history of depression is associated with increased connectivity in the Default Mode Network. J Affect Disord. (2018) 229:85–94. doi: 10.1016/j.jad.2017.12.052
29. Gong L, Yu SY, Xu RH, Liu D, Dai XJ, Wang ZY, et al. The abnormal reward network associated with insomnia severity and depression in chronic insomnia disorder. Brain Imaging Behav. (2021) 15:1033–42. doi: 10.1007/s11682-020-00310-w
30. Leerssen J, Blanken TF, Pozzi E, Jahanshad N, Aftanas L, Andreassen OA, et al. Brain structural correlates of insomnia severity in 1053 individuals with major depressive disorder: results from the ENIGMA MDD Working Group. Transl Psychiatry. (2020) 10:425. doi: 10.1038/s41398-020-01109-5
31. Yu S, Shen Z, Lai R, Feng F, Guo B, Wang Z, et al. The orbitofrontal cortex gray matter is associated with the interaction between insomnia and depression. Front Psychiatry. (2018) 9:651. doi: 10.3389/fpsyt.2018.00651
32. Gong L, Xu R, Liu D, Zhang C, Huang Q, Zhang B, et al. Abnormal functional connectivity density in patients with major depressive disorder with comorbid insomnia. J Affect Disord. (2020) 266:417–23. doi: 10.1016/j.jad.2020.01.088
33. Cui X, Guo W, Wang Y, Yang TX, Yang XH, Wang Y, et al. Aberrant default mode network homogeneity in patients with first-episode treatment-naive melancholic depression. Int J Psychophysiol. (2017) 112:46–51. doi: 10.1016/j.ijpsycho.2016.12.005
34. Parker G, Fink M, Shorter E, Taylor MA, Akiskal H, Berrios G, et al. Issues for DSM-5: whither melancholia? The case for its classification as a distinct mood disorder. Am J Psychiatry. (2010) 167:745–7. doi: 10.1176/appi.ajp.2010.09101525
35. Armitage R. Sleep and circadian rhythms in mood disorders. Acta Psychiatr Scand. (2007) 115:104–15. doi: 10.1111/j.1600-0447.2007.00968.x
36. Antonijevic I. HPA axis and sleep: identifying subtypes of major depression. Stress Int J Biol Stress. (2008) 11:15–27. doi: 10.1080/10253890701378967
37. Agargun MY, Besiroglu L, Cilli AS, Gulec M, Aydin A, Incl R, et al. Nightmares, suicide attempts, and melancholic features in patients with unipolar major depression. J Affect Disord. (2007) 98:267–70. doi: 10.1016/j.jad.2006.08.005
38. Besiroglu L, Agargun MY, Inci R. Nightmares and terminal insomnia in depressed patients with and without melancholic features. Psychiatry Res. (2005) 133:285–7. doi: 10.1016/j.psychres.2004.12.001
39. Wichniak A, Wierzbicka A, Jernajczyk W. Sleep as a biomarker for depression. Int Rev Psychiatry (Abingdon, England). (2013) 25:632–45. doi: 10.3109/09540261.2013.812067
40. Harald B, Gordon P. Meta-review of depressive subtyping models. J Affect Disord. (2012) 139:126–40. doi: 10.1016/j.jad.2011.07.015
41. Bruun CF, Arnbjerg CJ, Kessing LV. Electroencephalographic parameters differentiating melancholic depression, non-melancholic depression, and healthy controls. a systematic review. Front Psychiatry. (2021) 12:648713. doi: 10.3389/fpsyt.2021.648713
42. Biswal B, Yetkin FZ, Haughton VM, Hyde JS. Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Mag Reson Med. (1995) 34:537–41. doi: 10.1002/mrm.1910340409
43. Zhang Y, Cui X, Ou Y, Liu F, Li H, Chen J, et al. Differentiating melancholic and non-melancholic major depressive disorder using fractional amplitude of low-frequency fluctuations. Front Psychiatry. (2021) 12:763770. doi: 10.3389/fpsyt.2021.763770
44. Yan M, He Y, Cui X, Liu F, Li H, Huang R, et al. Disrupted regional homogeneity in melancholic and non-melancholic major depressive disorder at rest. Front Psychiatry. (2021) 12:618805. doi: 10.3389/fpsyt.2021.618805
45. Yan M, Cui X, Liu F, Li H, Huang R, Tang Y, et al. Abnormal default-mode network homogeneity in melancholic and nonmelancholic major depressive disorder at rest. Neural Plasticity. (2021) 2021:6653309. doi: 10.1155/2021/6653309
46. Shan X, Cui X, Liu F, Li H, Huang R, Tang Y, et al. Shared and distinct homotopic connectivity changes in melancholic and non-melancholic depression. J Affect Disord. (2021) 287:268–75. doi: 10.1016/j.jad.2021.03.038
47. Tian S, Zhu R, Chattun MR, Wang H, Chen Z, Zhang S, et al. Temporal dynamics alterations of spontaneous neuronal activity in anterior cingulate cortex predict suicidal risk in bipolar II patients. Brain Imag Behav. (2021) 15:2481–91. doi: 10.1007/s11682-020-00448-7
48. Zhang L, Zhang R, Han S, Womer FY, Wei Y, Duan J, et al. Three major psychiatric disorders share specific dynamic alterations of intrinsic brain activity. Schizophr Res. (2021) 243:322–9. doi: 10.1016/j.schres.2021.06.014
49. Calhoun VD, Miller R, Pearlson G, Adalı T. The chronnectome: time-varying connectivity networks as the next frontier in fMRI data discovery. Neuron. (2014) 84:262–74. doi: 10.1016/j.neuron.2014.10.015
50. Liao W, Zhang Z, Mantini D, Xu Q, Ji GJ, Zhang H, et al. Dynamical intrinsic functional architecture of the brain during absence seizures. Brain Struct Funct. (2014) 219:2001–15. doi: 10.1007/s00429-013-0619-2
51. Zheng R, Chen Y, Jiang Y, Wen M, Zhou B, Li S, et al. Dynamic altered amplitude of low-frequency fluctuations in patients with major depressive disorder. Front Psychiatry. (2021) 12:683610. doi: 10.3389/fpsyt.2021.683610
52. Liang Y, Jiang X, Zhu W, Shen Y, Xue F, Li Y, et al. Disturbances of dynamic function in patients with bipolar disorder i and its relationship with executive-function deficit. Front Psychiatry. (2020) 11:537981. doi: 10.3389/fpsyt.2020.537981
53. Li Q, Cao X, Liu S, Li Z, Wang Y, Cheng L, et al. Dynamic alterations of amplitude of low-frequency fluctuations in patients with drug-naïve first-episode early onset schizophrenia. Front Neurosci. (2020) 14:901. doi: 10.3389/fnins.2020.00901
54. van den Noort M, Bosch P, Yeo S, Lim S. Emotional memory processing: which comes first - depression or poor sleep? Sleep Med. (2016) 22:100. doi: 10.1016/j.sleep.2015.07.029
55. Beauregard M, Paquette V, Lévesque J. Dysfunction in the neural circuitry of emotional self-regulation in major depressive disorder. Neuroreport. (2006) 17:843–6. doi: 10.1097/01.wnr.0000220132.32091.9f
56. Hamilton M. A rating scale for depression. J Neurol Neurosurg Psychiatry. (1960) 23:56–62. doi: 10.1136/jnnp.23.1.56
57. Hamilton M. The assessment of anxiety states by rating. Br J Med Psychol. (1959) 32:50–5. doi: 10.1111/j.2044-8341.1959.tb00467.x
58. Park SC, Kim JM, Jun TY, Lee MS, Kim JB, Jeong SH, et al. Prevalence and clinical correlates of insomnia in depressive disorders: the crescend study. Psychiatry Invest. (2013) 10:373–81. doi: 10.4306/pi.2013.10.4.373
59. Manber R, Blasey C, Arnow B, Markowitz JC, Thase ME, Rush AJ, et al. Assessing insomnia severity in depression: comparison of depression rating scales and sleep diaries. J Psychiatr Res. (2005) 39:481–8. doi: 10.1016/j.jpsychires.2004.12.003
60. Yan CG, Wang XD, Zuo XN, Zang YFDPABI. Data processing & analysis for (resting-state) brain imaging. Neuroinformatics. (2016) 14:339–51. doi: 10.1007/s12021-016-9299-4
61. Li J, Duan X, Cui Q, Chen H, Liao W. More than just statics: temporal dynamics of intrinsic brain activity predicts the suicidal ideation in depressed patients. Psychol Med. (2019) 49:852–60. doi: 10.1017/s0033291718001502
62. Chen X, Lu B, Yan CG. Reproducibility of R-fMRI metrics on the impact of different strategies for multiple comparison correction and sample sizes. Hum Brain Mapp. (2018) 39:300–18. doi: 10.1002/hbm.23843
63. Li S, Wang BA, Li C, Feng Y, Li M, Wang T, et al. Progressive gray matter hypertrophy with severity stages of insomnia disorder and its relevance for mood symptoms. Eur Radiol. (2021) 31:6312–22. doi: 10.1007/s00330-021-07701-7
64. Wu Y, Zhuang Y, Qi J. Explore structural and functional brain changes in insomnia disorder: A PRISMA-compliant whole brain ALE meta-analysis for multimodal MRI. Medicine. (2020) 99:e19151. doi: 10.1097/md.0000000000019151
65. Wang YK, Shi XH, Wang YY, Zhang X, Liu HY, Wang XT, et al. Evaluation of the age-related and gender-related differences in patients with primary insomnia by fractional amplitude of low-frequency fluctuation: a resting-state functional magnetic resonance imaging study. Medicine. (2020) 99:e18786. doi: 10.1097/md.0000000000018786
66. Riemann D, Spiegelhalder K, Feige B, Voderholzer U, Berger M, Perlis M, et al. The hyperarousal model of insomnia: a review of the concept and its evidence. Sleep Med Rev. (2010) 14:19–31. doi: 10.1016/j.smrv.2009.04.002
67. Levenson JC, Kay DB, Buysse DJ. The pathophysiology of insomnia. Chest. (2015) 147:1179–92. doi: 10.1378/chest.14-1617
68. Gallagher HL, Frith CD. Functional imaging of ‘theory of mind’. Trends Cogn Sci. (2003) 7:77–83. doi: 10.1016/s1364-6613(02)00025-6
69. Blank TS, Meyer BM, Wieser MK, Rabl U, Schögl P, Pezawas L. Brain morphometry and connectivity differs between adolescent- and adult-onset major depressive disorder. Depress Anxiety. (2022) 39:387–96. doi: 10.1002/da.23254
70. Wu H, Zheng Y, Zhan Q, Dong J, Peng H, Zhai J, et al. Covariation between spontaneous neural activity in the insula and affective temperaments is related to sleep disturbance in individuals with major depressive disorder. Psychol Med. (2021) 51:731–40. doi: 10.1017/s0033291719003647
71. Liu J, Ren L, Womer FY, Wang J, Fan G, Jiang W, et al. Alterations in amplitude of low frequency fluctuation in treatment-naïve major depressive disorder measured with resting-state fMRI. Hum Brain Mapp. (2014) 35:4979–88. doi: 10.1002/hbm.22526
72. Via E, Cardoner N, Pujol J, Martínez-Zalacaín I, Hernández-Ribas R, Urretavizacaya M, et al. Cerebrospinal fluid space alterations in melancholic depression. PLoS One. (2012) 7:e38299. doi: 10.1371/journal.pone.0038299
73. Greenberg DL, Payne ME, MacFall JR, Steffens DC, Krishnan RR. Hippocampal volumes and depression subtypes. Psychiatry Res. (2008) 163:126–32. doi: 10.1016/j.pscychresns.2007.12.009
74. Maller JJ, Thomson RH, Rosenfeld JV, Anderson R, Daskalakis ZJ, Fitzgerald PB. Occipital bending in depression. Brain J Neurol. (2014) 137(Pt 6):1830–7. doi: 10.1093/brain/awu072
75. Palejwala AH, Dadario NB, Young IM, O’Connor K, Briggs RG, Conner AK, et al. Anatomy and white matter connections of the lingual gyrus and cuneus. World Neurosurg. (2021) 151:e426–37. doi: 10.1016/j.wneu.2021.04.050
76. Garrett A, Kelly R, Gomez R, Keller J, Schatzberg AF, Reiss AL. Aberrant brain activation during a working memory task in psychotic major depression. Am J Psychiatry. (2011) 168:173–82. doi: 10.1176/appi.ajp.2010.09121718
77. Tao H, Guo S, Ge T, Kendrick KM, Xue Z, Liu Z, et al. Depression uncouples brain hate circuit. Mol Psychiatry. (2013) 18:101–11. doi: 10.1038/mp.2011.127
78. Makovski T, Lavidor M. Stimulating occipital cortex enhances visual working memory consolidation. Behav Brain Res. (2014) 275:84–7. doi: 10.1016/j.bbr.2014.09.004
79. Colich NL, Foland-Ross LC, Eggleston C, Singh MK, Gotlib IH. Neural aspects of inhibition following emotional primes in depressed adolescents. J Clin child Adolesc Psychol. (2016) 45:21–30. doi: 10.1080/15374416.2014.982281
Keywords: melancholic depression, resting-state, magnetic resonance imaging, the dynamic amplitude of low-frequency fluctuation, sleep disturbance, insomnia
Citation: Deng Z, Jiang X, Liu W, Zhao W, Jia L, Sun Q, Xie Y, Zhou Y, Sun T, Wu F, Kong L and Tang Y (2022) The aberrant dynamic amplitude of low-frequency fluctuations in melancholic major depressive disorder with insomnia. Front. Psychiatry 13:958994. doi: 10.3389/fpsyt.2022.958994
Received: 01 June 2022; Accepted: 29 July 2022;
Published: 22 August 2022.
Edited by:
Ning Hua Wang, Fourth Military Medical University, ChinaCopyright © 2022 Deng, Jiang, Liu, Zhao, Jia, Sun, Xie, Zhou, Sun, Wu, Kong and Tang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yanqing Tang, tangyanqing@cmu.edu.cn
†These authors have contributed equally to this work and share first authorship