More Than Just Statics: Temporal Dynamic Changes in Inter- and Intrahemispheric Functional Connectivity in First-Episode, Drug-Naive Patients With Major Depressive Disorder
Yu Jiang1,2,3,4,5,6,7† 
Yuan Chen1,2,3,4,5,6,7† 
Ruiping Zheng1,2,3,4,5,6,7 Bingqian Zhou1,2,3,4,5,6,7 
Ying Wei1,2,3,4,5,6,7 
Ankang Gao1,2,3,4,5,6,7 
Yarui Wei1,2,3,4,5,6,7 Shuying Li8 
Jinxia Guo9 
Shaoqiang Han1,2,3,4,5,6,7* 
Yong Zhang1,2,3,4,5,6,7* 
Jingliang Cheng1,2,3,4,5,6,7*
- 1Department of Magnetic Resonance Imaging, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
- 2Key Laboratory for Functional Magnetic Resonance Imaging and Molecular Imaging of Henan Province, Zhengzhou, China
- 3Engineering Technology Research Center for Detection and Application of Brain Function of Henan Province, Zhengzhou, China
- 4Engineering Research Center of Medical Imaging Intelligent Diagnosis and Treatment of Henan Province, Zhengzhou, China
- 5Key Laboratory of Magnetic Resonance and Brain Function of Henan Province, Zhengzhou, China
- 6Key Laboratory of Brain Function and Cognitive Magnetic Resonance Imaging of Zhengzhou, Zhengzhou, China
- 7Key Laboratory of Imaging Intelligence Research Medicine of Henan Province, Zhengzhou, China
- 8Department of Psychiatry, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
- 9GE Healthcare, Beijing, China
Several functional magnetic resonance imaging (fMRI) studies have demonstrated abnormalities in static intra- and interhemispheric functional connectivity among diverse brain regions in patients with major depressive disorder (MDD). However, the dynamic changes in intra- and interhemispheric functional connectivity patterns in patients with MDD remain unclear. Fifty-eight first-episode, drug-naive patients with MDD and 48 age-, sex-, and education level-matched healthy controls (HCs) underwent resting-state fMRI. Whole-brain functional connectivity, analyzed using the functional connectivity density (FCD) approach, was decomposed into ipsilateral and contralateral functional connectivity. We computed the intra- and interhemispheric dynamic FCD (dFCD) using a sliding window analysis to capture the dynamic patterns of functional connectivity. The temporal variability in functional connectivity was quantified as the variance of the dFCD over time. In addition, intra- and interhemispheric static FCD (sFCD) patterns were calculated. Associations between the dFCD variance and sFCD in abnormal brain regions and the severity of depressive symptoms were analyzed. Compared to HCs, patients with MDD showed lower interhemispheric dFCD variability in the inferior/middle frontal gyrus and decreased sFCD in the medial prefrontal cortex/anterior cingulate cortex and posterior cingulate cortex/precuneus in both intra- and interhemispheric comparisons. No significant correlations were found between any abnormal dFCD variance or sFCD at the intra- and interhemispheric levels and the severity of depressive symptoms. Our results suggest intra- and interhemispheric functional connectivity alterations in the dorsolateral prefrontal cortex (DLPFC) and default mode network regions involved in cognition, execution and emotion. Furthermore, our study emphasizes the essential role of altered interhemispheric communication dynamics in the DLPFC in patients with MDD. These findings contribute to our understanding of the pathophysiology of MDD.
Introduction
Major depressive disorder (MDD), characterized by anhedonia, a persistent negative mood, and cognitive impairments, is a common mental disorder (Lu et al., 2013). It has been reported that MDD is the second leading cause of disability in China (Lu et al., 2021). The high recurrence of MDD and the severe distress experienced by the affected individuals leads to a high risk of suicide, which increases the societal and economic burden of the disease (Nock et al., 2009; Richards, 2011). A better understanding of the pathophysiology of MDD is therefore essential to improve therapeutic outcomes and reduce the social burden.
Functional magnetic resonance imaging (fMRI) provides a new direction for investigating resting-state functional connectivity (FC) and offers effective information for evaluating the integration of functional neural networks in vivo. Altered whole-brain FC patterns in MDD have been observed in some resting-state fMRI studies (Zhuo et al., 2017, 2020; Shi et al., 2020), and the FC abnormalities were found to be closely related to cognitive impairments and emotional disorders in patients with MDD (Shi et al., 2020). Moreover, researchers have observed decreased long-range connectivity in the visual cortex, local connectivity patterns in the frontal and temporal cortex in patients with MDD (Zou et al., 2016), and associations between altered local connectivity in MDD before and after antidepressant treatment and the improvements in clinical severity (Wang et al., 2018). According to these fMRI studies, the characteristics of under-connectivity and over-connectivity are effective markers for MDD.
Intra- and interhemispheric interactions play an important role in ensuring fast and efficient computation and information processing in many perceptual and cognitive functions including facial recognition, language understanding, and visuospatial processing (Ringo et al., 1994; Kong et al., 2018). Some fMRI studies have reported high levels of interhemispheric FC in the sensorimotor and auditory cortex in human brain (Biswal et al., 1995; Cordes et al., 2000). Intrahemispheric FC was also reported in earlier studies (Quigley et al., 2003), reflecting the phenomenon of intrahemispheric integration. Some studies have suggested that the left hemisphere may be dominant for cognitive processing, while the right hemisphere may be more important for affective processing (Wittling and Roschmann, 1993; Gainotti, 2012). Some mental disorders, including schizophrenia (Zhu et al., 2018) and autism spectrum disorders (Lee et al., 2016), have been related to altered intra- and interhemispheric interactions. Abnormal hemispheric communication patterns have also been observed in patients with MDD (Guo et al., 2018; Ran et al., 2020). Earlier studies have mainly focused on interhemispheric connectivity, namely, voxel-mirrored homotopic connectivity, which is the correlation between each voxel and the corresponding voxel in the opposite hemisphere (Zuo et al., 2010). Using this measure, individuals with MDD have been found to show disruptions in interhemispheric FC, suggesting impairments in interhemispheric communication during cognitive and emotional processing (Guo et al., 2013, 2018; Lai and Wu, 2014; Hermesdorf et al., 2016; Hou et al., 2016). A few resting-state fMRI studies have observed intra- and interhemispheric FC dysfunction in some cognitive and emotional areas in MDD (Jiang et al., 2019; Ran et al., 2020; Ding et al., 2021). These studies indicate that such impairments in intra- and interhemispheric FC are strongly linked to deficits in perceptual, cognitive, and emotional processing in MDD, and that they may thus be a powerful internal brain activity metric that can help our understanding of MDD as well as aid in the discovery of the underlying neural resting-state pathway.
However, these earlier studies have relied on the assumption that functional connections are static during resting-state fMRI scans and have ignored the presence and potential of FC dynamics in patients with MDD. Recent studies on FC have considered the variability in spatial organization of the brain over time, which can be quantified by measuring the temporal variability of neural signals (Zhu et al., 2020; Tang et al., 2022). FC dynamics may be even more evident when mental activity is unconstrained in the resting-state. Some researchers believe that individuals freely engage in some types of mental activity during the conscious resting state (Delamillieure et al., 2010), and that the predominant mental activity influences FC and modular organization throughout the brain (Doucet et al., 2012). Early studies emphasized that dynamic changes in FC may be related to changes in the state of consciousness (Hutchison et al., 2013), in arousal (Chang et al., 2013), in emotional state (Cribben et al., 2012), and in vigilance (Patanaik et al., 2018). The altered temporal characteristics of dynamic FC may be potential biomarkers of various neurological and psychiatric disorders, including Alzheimer’s disease (Jones et al., 2012), schizophrenia (Sakoğlu et al., 2010), and MDD (Cui et al., 2020). By virtue of the significant dynamic FC abnormalities of specific brain regions and networks in MDD which accurately capture changes in neural activity (Zhang et al., 2020; Zhu et al., 2020), characterizing dynamic FC at intra- and interhemispheric levels is necessary to acquire more worthwhile information about changes in brain states and network properties. Altered temporal dynamics in intra- and interhemispheric FC patterns have been reported in autistic children, suggesting that this may be behind the impairments in social communication characteristic of autism spectrum disorders (Guo et al., 2020). To date, the dynamic changes in intra- and interhemispheric FC patterns in MDD remain unclear. Further investigation of FC dynamics may help to offer a novel understanding and explanation of abnormal brain communication in MDD.
Functional connectivity density (FCD) mapping is a data-driven graph theory method that describes the number of connections of each voxel with all other voxels and identifies the distribution of major connected hubs in the whole brain (Tomasi and Volkow, 2010). In comparison with independent component analysis and seed-based FC methods that rely on specific regions or networks, FCD mapping can evaluate the functional properties of the connectomics of the entire brain, thus evaluating the brain’s information communication capability (Zuo et al., 2012). Static FCD (sFCD) using the entire time series reflects the characteristics of FC within whole-brain networks, while dynamic FCD (dFCD) can be utilized to detect abnormalities in brain connectivity over time by integrating the FCD with a sliding window correlation approach (Li et al., 2019). After computing the FCD within each window, the temporal variability in the dFCD could characterize the time-varying FC patterns within the brain networks. To the best of our knowledge, alterations in the dFCD patterns in MDD remain unclear. Since dynamic and static FC patterns have been reported to offer overlapping or complementary information (Qiao et al., 2020), we combined dFCD and sFCD methods to obtain a more comprehensive explanation of the intra- and interhemispheric FC patterns observed in MDD.
In the present study, we thus decomposed the whole-brain FCD into ipsilateral and contralateral parts, which represent intra- and interhemispheric FC, respectively. Similar to global FCD, intra- and interhemispheric FCD can be used to identify the distribution of major connected hubs at hemispheric levels in brain networks. The associations between dFCD variance and sFCD in abnormal brain regions and the severity of depressive symptoms were also analyzed. We hypothesized that patients with MDD would exhibit abnormal dFCD and sFCD patterns of intra- and interhemispheric connections in some brain regions and specific dFCD variance and sFCD changes at intra- and interhemispheric levels related to the severity of their depressive symptoms.
Materials and Methods
Participants
A total of 106 individuals were enrolled in this study, including 58 first-episode drug-naive patients with MDD and 48 age-, sex-, and education level-matched healthy controls (HCs). MDD was diagnosed by two trained psychiatrists independently, according to the Structured Clinical Interview of the Diagnostic and Statistical Manual of Mental Disorders (SCID-I/P, Chinese version). To rule out the presence of current or past psychiatric illness, participants were carefully screened through a SCID-I/P. The severity of depressive symptoms was assessed using the 24-items Hamilton Depression Scale (HAMD-24), an auxiliary tool for evaluating clinical status (Hamilton, 1960). Higher HAMD-24 scores correspond to more severe the symptoms of depression. Patients with MDD were experiencing a depressive episode at the time of scanning. The inclusion criteria for patients with MDD were as follows: (1) fulfilled the SCID-I/P diagnostic criteria for MDD; (2) first episode, and never received any antidepressant treatment; (3) a cut-off score ≥ 20 on the HAMD-24; and (4) Han ethnicity and right-handed. HCs were recruited from the local community via poster advertisements. The exclusion criteria for all participants included in the current study were as follows: (1) a history of other mental diseases; (2) organic brain disease or serious physical disease; (3) drug abuse or alcohol dependence; and (4) contraindications to MRI examination, such as metal implants or claustrophobia.
Table 1 shows the demographic and clinical characteristics of the MDD and HCs groups. There were no statistical differences between the two groups in terms of age (t = −0.950, P = 0.345), sex (χ2 = 2.045, P = 0.153), years of education (t = −1.133, P = 0.261), or mean framewise displacement (FD) (t = −1.728, P = 0.088). The range of HAMD-24 scores in the MDD group was 20–69.
Table 1. Demographics and clinical characteristics of all participants.
Written informed consent was obtained from all participants. The study was approved by the Medical Ethics Committee of The First Affiliated Hospital of Zhengzhou University.
Image Acquisition
MRI data were acquired using a GE Discovery MR750 3.0T scanner (General Electric, Milwaukee, WI, United States), with an eight-channel prototype quadrature birdcage head coil. All participants lay in a supine position and had their head fixed with foam pads to reduce head movement. Earplugs were used to reduce noise stimulation. During scanning, all participants were instructed to keep their eyes closed, stay awake, and avoid intentional thinking. Resting-state fMRI data were obtained using a gradient-echo echo-planar imaging sequence with the following scanning parameters: repetition time/echo time (TR/TE) = 2,000/40 ms, field of view (FOV) = 220 mm × 220 mm, slices = 32, matrix size = 64 × 64, section thickness = 4 mm, gap = 0.5 mm, flip angle (FA) = 90°, and 180 volumes (360 s in total). Structural T1-weighted images were obtained with 3D spoiled gradient echo scan sequence: TR/TE = 8,164/3.18 ms, inversion time = 900 ms, FA = 7°, resolution matrix = 256 × 256, thickness = 1.0 mm, and slices = 188.
Data Preprocessing
All preprocessing was conducted using the Data Processing and Analysis for Brain Imaging (DPABI)1 toolbox. The first 10 volumes were discarded to ensure signal stability, followed by slice timing and realignment. Head movement exceeding 2.5 mm in maximum displacement or 2.5° rotation in angular motion were excluded. FD was computed, and data with mean FD values exceeding 0.5 mm were excluded (Power et al., 2012). Images were normalized to the standard Montreal Neurological Institute (MNI) space, and each voxel was resampled to 3 mm × 3 mm × 3 mm. The normalized images were detrended to correct the linear trend. Then several spurious variances including 24 head motion parameters, cerebrospinal fluid signal and white matter signal were regressed out from the data. The global signal was not regressed, since previous studies have reported that global signal regression may introduce distortions in time series data (Anderson et al., 2011), and that global signal variance and topography can show alterations in mental disorders, including MDD (Yang et al., 2014; Han et al., 2019; Zhang et al., 2019). Furthermore, band pass filtering was carried out at a frequency range of 0.01–0.08 Hz. Finally, scrubbing with cubic spline interpolation was performed to eliminate the influence of head motion and ensure continuous time points. The mean FD was entered as a covariate into the group-level analysis.
Inter- and Intrahemispheric Functional Connectivity Density Analysis
An sFCD analysis was performed to evaluate global, contralateral, and ipsilateral FC. Automated anatomical labeling (AAL) atlas excluding the cerebellum was used as a cerebral gray matter mask. The data-driven FCD mapping method was applied to calculate the global sFCD, which was defined as the number of connections of each voxel with all other voxels in the brain mask. Then, the global sFCD was decomposed into two parts according to the relative positions of the seed and target voxels: First, into the contralateral (or interhemispheric) sFCD, which refers to the number of FC of a given voxel with all voxels in the opposite hemispheric brain regions, and second, into the ipsilateral (or intrahemispheric) sFCD, which represents the number of FC with all voxels in the same hemispheric brain regions. To evaluate the connectivity between two voxels, linear Pearson correlations were performed, and the correlation threshold was set in accordance with the criterion of p < 0.05 (uncorrected) (Zou et al., 2016; Guo et al., 2020; Yang et al., 2021). When the correlation coefficient between two voxels is larger than the correlation threshold, connectivity is considered to be present.
In the dFCD analysis, we applied a sliding window approach to calculate the temporal variability of global, contralateral, and ipsilateral FC. We chose a sliding window ranging from 10 to 180 s, based on a previous study that suggested this size to be appropriate for capturing the dynamic characteristics of intrinsic brain activities (Gonzalez-Castillo et al., 2015). Prior studies have also demonstrated that to avoid introducing spurious fluctuations, the minimum window length should be above 1/fmin, where fmin is the minimum frequency of the time series (Leonardi and Van De Ville, 2015). Based on previous research (Guo et al., 2020; Luo et al., 2021; Zheng et al., 2021), a window length of 50 TRs (100s) and a step size of 2 TRs (4s) were selected to compute the temporal variability of FCD. The full-length time course was segmented into 61 windows for each subject. The correlation threshold of dFCD was set to be identical to that in the sFCD analysis (Guo et al., 2020). The global FCD calculated within each window was decomposed into contralateral and ipsilateral FCD, and then the variance of the time-varying FCD patterns over time was calculated to estimate temporal variability in global, contralateral, and ipsilateral dFCD.
Statistical and Correlation Analysis
For each subject, the dFCD variance and sFCD maps were converted to z-values using Fisher z-transformation to enhance the normality of the data distribution. Then, all normalized images were smoothed with a 6 mm × 6 mm × 6 mm full-width at half-maximum (FWHM) Gaussian kernel. A two-sample t-test was employed on the FCD maps to analyze the differences between the MDD and HCs groups. Age, sex, years of education, and mean FD were used as covariates. In order to control the false positive rate, multiple comparisons were corrected using a Gaussian random field theory (GRF) approach at a voxel-wise threshold of p < 0.005 with the cluster-level threshold set to p < 0.05 (the cluster-level correction refers to family wise error), and cluster-extent threshold was set at 20 voxels.
The brain clusters showing significant global, intra-, and interhemispheric FCD differences between patients with MDD and HCs were defined as regions of interest (ROIs). The values of dFCD variance and sFCD at global and intra- and interhemispheric levels in all ROIs were extracted separately. A Pearson correlation analysis was performed to identify correlations of abnormal dFCD variance and sFCD with depressive symptom severity (p < 0.05, uncorrected).
Validation Analysis
Two additional window lengths (30 and 80 TRs) were applied to test the robustness of the sliding window analysis. Moreover, dFCD variability with 1 TR as the shifting step and 50 TRs as the window length was calculated to estimate the confounding influence of different step sizes on our main dFCD findings.
Results
Inter- and Intrahemispheric Dynamic Functional Connectivity Density Variance
Compared with the HCs group, the MDD group showed decreased global dFCD variance in the inferior frontal gyrus (IFG) and the middle frontal gyrus (MFG) (Figure 1 and Table 2). Furthermore, the comparisons of the contralateral dFCD variance maps showed patterns that were similar to the global dFCD variance (Figure 1 and Table 2). The MDD group showed decreased contralateral dFCD variance in the IFG/MFG compared to the HCs group (Figure 1 and Table 2). No significant difference in ipsilateral dFCD variability was found between the MDD and HCs groups.
Figure 1. The group difference for (A) global and (B) contralateral dFCD variance between the MDD and HCs groups. The color bars indicate the T value based on two-sample t-test. The statistical significance level was set at Pvoxel < 0.005, Pcluster < 0.05 under Gaussian random field theory (GRF) corrected. The histogram shows the average (C) global and (D) contralateral dFCD variance in the IFG/MFG. *p < 0.05. These results are mean ± standard error. Abbreviations: dFCD, dynamic functional connectivity density; IFG, inferior frontal gyrus; MFG, middle frontal gyrus; MDD, major depressive disorder; HCs, healthy controls; L, left; R, right.
Table 2. Altered global and contralateral dFCD variance in MDD group.
Inter- and Intrahemispheric Static Functional Connectivity Density
As shown in Figure 2 and Table 3, patients with MDD showed decreased global sFCD in the anterior default mode network (DMN) regions, including the medial prefrontal cortex (mPFC) and anterior cingulate cortex (ACC), as well as in posterior DMN regions including the posterior cingulate cortex (PCC) and precuneus (PCu). The group differences in the contralateral sFCD and ipsilateral sFCD revealed patterns similar to those of the global sFCD. Both decreased contralateral sFCD and ipsilateral sFCD in the mPFC/ACC and PCC/PCu were found in patients with MDD (Figure 2 and Table 3).
Figure 2. The group difference for (A) global, (B) contralateral and (C) ipsilateral sFCD between the MDD and HCs groups. The color bars indicate the T value based on two-sample t-test. The statistical significance level was set at Pvoxel < 0.005, Pcluster < 0.05 under Gaussian random field theory (GRF) corrected. The histogram shows the average (D) global, (E) contralateral and (F) ipsilateral sFCD values in the mPFC/ACC and PCC/PCu. *p < 0.05. These results are mean ± standard error. Abbreviations: sFCD, static functional connectivity density; mPFC, medial prefrontal cortex; ACC, anterior cingulate cortex; PCC, posterior cingulate cortex; PCu, precuneus; MDD, major depressive disorder; HCs, healthy controls; L, left; R, right.
Table 3. Altered global, contralateral, and ipsilateral sFCD in MDD group.
Correlation Analysis
There were no significant correlations between dFCD variance or sFCD at global, intra-, and interhemispheric levels in abnormal brain regions and the total HAMD-24 scores (p > 0.05, uncorrected).
Validation Analysis
The results of the contralateral dFCD variance analysis were consistent with our main dFCD results. For the global dFCD analysis, the findings in the IFG/MFG regions with a shifting step of 1 TR were successfully repeated, while the results of the two additional window lengths of 30 and 80 TRs were not preserved after GRF correction. The detailed results of the validation analysis are presented in the Supplementary Material.
Discussion
Since dynamic and static FC patterns have been reported to offer overlapping or complementary information (Qiao et al., 2020), the current study employed both dFCD and sFCD to compare intra- and interhemispheric FC patterns between patients with MDD and HCs. We found that patients with MDD exhibited decreased interhemispheric dFCD variability in the dorsolateral prefrontal cortex (DLPFC) regions, including the IFG/MFG, which is associated with attention, emotion, and cognitive control (Aron et al., 2014; Song et al., 2019). Both intra- and interhemispheric comparisons showed reduced sFCD in regions related to attention focus and cognitive control, namely the DMN regions (mPFC/ACC and PCC/PCu) in patients with MDD.
To the best of our knowledge, the present study is the first to explore the dynamic abnormalities of intra- and interhemispheric FC in the intrinsic neural networks of patients with MDD. Previous studies have focused on static FC abnormalities at intra- and interhemispheric levels in MDD (Jiang et al., 2019; Ran et al., 2020; Ding et al., 2021), and have suggested the significance of exploring intra- and interhemispheric information transfer pathways in patients with MDD. Our results thus expand these previous static FC findings into the realm of dynamic FC, by using the whole-brain dFCD approach at intra- and interhemispheric levels. Notably, our study reveals abnormal FC patterns of time-varying intrinsic brain functional networks in MDD, which is in accordance with prior studies on abnormal interaction dynamics in the DLPFC and other brain regions in patients with MDD (Cui et al., 2020; Zhang et al., 2020). The combination of dynamic and static functional connectome analyses employed here reflects abnormal integration of intra- and interhemispheric brain regions and therefore point the way for future research in patients with MDD.
Patients with MDD showed decreased interhemispheric dFCD variability in the IFG/MFG relative to the HCs. The IFG and MFG are parts of the DLPFC, with the right IFG playing an important role in emotion and cognitive control (Aron et al., 2014). The right MFG is associated with sustained attention, working memory, and cognitive control functions, and may be the node that links the dorsal and ventral attention networks (Koyama et al., 2017; Yüksel et al., 2018; Song et al., 2019). A previous study showed decreased global sFCD in the IFG and MFG in patients with MDD (Zhuo et al., 2020). Another study by Shi et al. (2020) applied an FC strength (FCS) method to measure all the functional connection values of each voxel and demonstrated decreased global static FCS in the IFG in first-episode drug-naive patients with MDD, indicating that static FCS within DLPFC regions changes during the early stages of MDD (Shi et al., 2020). These studies indicate that the IFG and MFG are crucial hubs of MDD in the brain. Our results extend these findings using the dFCD method and offer nuanced information on the temporal changes in the major FC hubs, which may be masked in traditional static studies. Interhemispheric under-connectivity patterns, rather than intrahemispheric under-connectivity, may contribute to global under-connectivity in the MDD group. The invariable interhemispheric dFCD in the IFG/MFG region exhibited in our study further highlights abnormal dynamic functional coordination in MDD, which points to dysfunctions in attention, emotion, and cognitive control.
Our results show both intra- and interhemispheric sFCD alterations in the mPFC/ACC and PCC/PCu in patients with MDD. The mPFC/ACC is centrally located in the anterior DMN regions and is involved in social perception and cognitive control (Adolphs, 2009). A previous resting-state fMRI study showed decreased intra- and interhemispheric static FC in the mPFC and ACC in MDD, which is in line with our findings (Jiang et al., 2019). Decreased interhemispheric static FC between anatomically symmetrical voxels has also been reported in the mPFC (Guo et al., 2013; Wang et al., 2013; Lai and Wu, 2014) and ACC (Lai and Wu, 2014) in MDD. The PCC/PCu, belonging to a core region of the posterior DMN, plays an important role in episodic memory, attention focus, and internally directed cognition (Maddock et al., 2001; Leech and Sharp, 2014). A previous study showed decreased intra- and interhemispheric static FC in the PCC as well as decreased interhemispheric static FC in PCu in MDD (Ran et al., 2020). However, this study investigated intra- and interhemispheric static FC only in the DMN areas and based their analyses on a ROI-wise correlation approach, which most likely led to the differences in findings compared to our study. In addition, reduced interhemispheric static FC in the PCC (Guo et al., 2013, 2018; Wang et al., 2015) and PCu (Guo et al., 2013; Hermesdorf et al., 2016) have been found in MDD in previous studies. Our study demonstrated that there is no overlap between the results from dFCD and sFCD. The results of intra- and interhemispheric sFCD present complementary information that differs from the information provided by dFCD. By measuring sFCD, the present results demonstrate the impact of impairments in intra- and interhemispheric functional coordination and deepen our understanding of FC abnormalities in MDD.
Surprisingly, we did not find significant correlations between abnormal dFCD variance or sFCD and the severity of depressive symptoms in our patients. We hypothesize that this is attributable to the small sample size and narrow scoring range of HAMD-24. This suggests that while dFCD and sFCD can reflect abnormal intra- and interhemispheric functional connections, they may not be useful as biological indicators for the quantitative analysis of depression severity.
The present study has some limitations. First, while we selected 50 TRs as the window length for the main dFCD analysis to capture the dynamic patterns of intrinsic brain activities (Guo et al., 2020; Luo et al., 2021; Zheng et al., 2021), the optimal window length remains controversial, and the effects of different window lengths are unknown. Second, the sample size of our study is relatively small, which may restrict the generalizability of our findings. Third, future longitudinal investigations are necessary to further verify the findings of this cross-sectional design study. Fourth, our study included teenagers and adults, and the influence of developmental factors on the experimental results should be considered. Future research should focus on age-related changes in inter- and intrahemispheric FC in MDD. Fifth, the duration of the scan from which resting-state fMRI data was derived was 6 min, which is a relatively short period and may give rise to less reliable results for FC analysis (Birn et al., 2013).
Conclusion
The results of our dFCD analysis present complementary information that differs from the information that sFCD provides at the intra- and interhemispheric level. The present study demonstrates altered interhemispheric dFCD variability in IFG/MFG in patients with MDD. Additionally, we found global alterations in sFCD in DMN regions, including the mPFC/ACC and the PCC/PCu in our patients with MDD. Our findings reveal the impact that impairments in intra- and interhemispheric functional coordination have in MDD, which adds to our understanding of the pathophysiology of this disease.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Ethics Statement
The studies involving human participants were reviewed and approved by the Medical Ethics Committee of The First Affiliated Hospital of Zhengzhou University. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
YJ, YC, and SH conceived and designed the study. YJ, RZ, and BZ supervised the conduct of the study. YC, RZ, BZ, SL, YW, and AG are responsible for the data acquisition. YJ and YC drafted the initial manuscript, analyzed the data, and took responsibility for the manuscript. YJ, YC, BZ, and YRW assisted with the literature review. JC, YZ, SH, and JG reviewed and revised the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (81601467, 81871327, and 62106229) and Medical Science and Technology research project of Henan province (201701011).
Conflict of Interest
JG is an employee of GE Healthcare.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We would like to express our thanks to all subjects for their participation.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnhum.2022.868135/full#supplementary-material
Footnotes
References
Adolphs, R. (2009). The social brain: neural basis of social knowledge. Annu. Rev. Psychol. 60, 693–716. doi: 10.1146/annurev.psych.60.110707.163514
Anderson, J. S., Druzgal, T. J., Lopez-Larson, M., Jeong, E. K., Desai, K., and Yurgelun-Todd, D. (2011). Network anticorrelations, global regression, and phase-shifted soft tissue correction. Hum. Brain Mapp. 32, 919–934. doi: 10.1002/hbm.21079
Aron, A. R., Robbins, T. W., and Poldrack, R. A. (2014). Inhibition and the right inferior frontal cortex: one decade on. Trends Cogn. Sci. 18, 177–185. doi: 10.1016/j.tics.2013.12.003
Birn, R. M., Molloy, E. K., Patriat, R., Parker, T., Meier, T. B., Kirk, G. R., et al. (2013). The effect of scan length on the reliability of resting-state fMRI connectivity estimates. Neuroimage 83, 550–558. doi: 10.1016/j.neuroimage.2013.05.099
Biswal, B., Yetkin, F. Z., Haughton, V. M., and Hyde, J. S. (1995). Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn. Reson. Med. 34, 537–541. doi: 10.1002/mrm.1910340409
Chang, C., Metzger, C. D., Glover, G. H., Duyn, J. H., Heinze, H. J., and Walter, M. (2013). Association between heart rate variability and fluctuations in resting-state functional connectivity. Neuroimage 68, 93–104. doi: 10.1016/j.neuroimage.2012.11.038
Cordes, D., Haughton, V. M., Arfanakis, K., Wendt, G. J., Turski, P. A., Moritz, C. H., et al. (2000). Mapping functionally related regions of brain with functional connectivity MR imaging. AJNR Am. J. Neuroradiol. 21, 1636–1644.
Cribben, I., Haraldsdottir, R., Atlas, L. Y., Wager, T. D., and Lindquist, M. A. (2012). Dynamic connectivity regression: determining state-related changes in brain connectivity. Neuroimage 61, 907–920. doi: 10.1016/j.neuroimage.2012.03.070
Cui, G., Wang, Y., Wang, X., Zheng, L., Li, L., Li, P., et al. (2020). Static and dynamic functional connectivity of the prefrontal cortex during resting-state predicts self-serving bias in depression. Behav. Brain Res. 379:112335. doi: 10.1016/j.bbr.2019.112335
Delamillieure, P., Doucet, G., Mazoyer, B., Turbelin, M. R., Delcroix, N., Mellet, E., et al. (2010). The resting state questionnaire: an introspective questionnaire for evaluation of inner experience during the conscious resting state. Brain Res. Bull. 81, 565–573. doi: 10.1016/j.brainresbull.2009.11.014
Ding, Y. D., Yang, R., Yan, C. G., Chen, X., Bai, T. J., Bo, Q. J., et al. (2021). Disrupted hemispheric connectivity specialization in patients with major depressive disorder: evidence from the REST-meta-MDD Project. J. Affect. Disord. 284, 217–228. doi: 10.1016/j.jad.2021.02.030
Doucet, G., Naveau, M., Petit, L., Zago, L., Crivello, F., Jobard, G., et al. (2012). Patterns of hemodynamic low-frequency oscillations in the brain are modulated by the nature of free thought during rest. Neuroimage 59, 3194–3200. doi: 10.1016/j.neuroimage.2011.11.059
Gainotti, G. (2012). Unconscious processing of emotions and the right hemisphere. Neuropsychologia 50, 205–218. doi: 10.1016/j.neuropsychologia.2011.12.005
Gonzalez-Castillo, J., Hoy, C. W., Handwerker, D. A., Robinson, M. E., Buchanan, L. C., Saad, Z. S., et al. (2015). Tracking ongoing cognition in individuals using brief, whole-brain functional connectivity patterns. Proc. Natl. Acad. Sci. USA 112, 8762–8767. doi: 10.1073/pnas.1501242112
Guo, W., Cui, X., Liu, F., Chen, J., Xie, G., Wu, R., et al. (2018). Decreased interhemispheric coordination in the posterior default-mode network and visual regions as trait alterations in first-episode, drug-naive major depressive disorder. Brain Imag. Behav. 12, 1251–1258. doi: 10.1007/s11682-017-9794-8
Guo, W., Liu, F., Dai, Y., Jiang, M., Zhang, J., Yu, L., et al. (2013). Decreased interhemispheric resting-state functional connectivity in first-episode, drug-naive major depressive disorder. Prog. Neuropsychopharmacol. Biol. Psychiatry 41, 24–29. doi: 10.1016/j.pnpbp.2012.11.003
Guo, X., Duan, X., Chen, H., He, C., Xiao, J., Han, S., et al. (2020). Altered inter- and intrahemispheric functional connectivity dynamics in autistic children. Hum. Brain Mapp. 41, 419–428. doi: 10.1002/hbm.24812
Hamilton, M. (1960). A rating scale for depression. J. Neurol. Neurosurg. Psychiatry 23, 56–62. doi: 10.1136/jnnp.23.1.56
Han, S., Wang, X., He, Z., Sheng, W., Zou, Q., Li, L., et al. (2019). Decreased static and increased dynamic global signal topography in major depressive disorder. Prog. Neuropsychopharmacol. Biol. Psychiatry 94:109665. doi: 10.1016/j.pnpbp.2019.109665
Hermesdorf, M., Sundermann, B., Feder, S., Schwindt, W., Minnerup, J., Arolt, V., et al. (2016). Major depressive disorder: findings of reduced homotopic connectivity and investigation of underlying structural mechanisms. Hum. Brain Mapp. 37, 1209–1217. doi: 10.1002/hbm.23097
Hou, Z., Sui, Y., Song, X., and Yuan, Y. (2016). Disrupted interhemispheric synchrony in default mode network underlying the impairment of cognitive flexibility in late-onset depression. Front. Aging Neurosci. 8:230. doi: 10.3389/fnagi.2016.00230
Hutchison, R. M., Womelsdorf, T., Gati, J. S., Everling, S., and Menon, R. S. (2013). Resting-state networks show dynamic functional connectivity in awake humans and anesthetized macaques. Hum. Brain Mapp. 34, 2154–2177. doi: 10.1002/hbm.22058
Jiang, X., Shen, Y., Yao, J., Zhang, L., Xu, L., Feng, R., et al. (2019). Connectome analysis of functional and structural hemispheric brain networks in major depressive disorder. Transl. Psychiatry 9:136. doi: 10.1038/s41398-019-0467-9
Jones, D. T., Vemuri, P., Murphy, M. C., Gunter, J. L., Senjem, M. L., Machulda, M. M., et al. (2012). Non-stationarity in the “resting brain’s” modular architecture. PLoS One 7:e39731. doi: 10.1371/journal.pone.0039731
Kong, X. Z., Mathias, S. R., Guadalupe, T., Glahn, D. C., Franke, B., Crivello, F., et al. (2018). Mapping cortical brain asymmetry in 17,141 healthy individuals worldwide via the ENIGMA Consortium. Proc. Natl. Acad. Sci. USA 115, E5154–E5163. doi: 10.1073/pnas.1718418115
Koyama, M. S., O’Connor, D., Shehzad, Z., and Milham, M. P. (2017). Differential contributions of the middle frontal gyrus functional connectivity to literacy and numeracy. Sci. Rep. 7:17548. doi: 10.1038/s41598-017-17702-6
Lai, C. H., and Wu, Y. T. (2014). Decreased inter-hemispheric connectivity in anterior sub-network of default mode network and cerebellum: significant findings in major depressive disorder. Int. J. Neuropsychopharmacol. 17, 1935–1942. doi: 10.1017/s1461145714000947
Lee, J. M., Kyeong, S., Kim, E., and Cheon, K. A. (2016). Abnormalities of inter- and intra-hemispheric functional connectivity in autism spectrum disorders: a study using the autism brain imaging data exchange database. Front. Neurosci. 10:191. doi: 10.3389/fnins.2016.00191
Leech, R., and Sharp, D. J. (2014). The role of the posterior cingulate cortex in cognition and disease. Brain 137, 12–32. doi: 10.1093/brain/awt162
Leonardi, N., and Van De Ville, D. (2015). On spurious and real fluctuations of dynamic functional connectivity during rest. Neuroimage 104, 430–436. doi: 10.1016/j.neuroimage.2014.09.007
Li, R., Wang, L., Chen, H., Guo, X., Liao, W., Tang, Y. L., et al. (2019). Abnormal dynamics of functional connectivity density in children with benign epilepsy with centrotemporal spikes. Brain Imag. Behav. 13, 985–994. doi: 10.1007/s11682-018-9914-0
Lu, J., Xu, X., Huang, Y., Li, T., Ma, C., Xu, G., et al. (2021). Prevalence of depressive disorders and treatment in China: a cross-sectional epidemiological study. Lancet Psychiatry 8, 981–990. doi: 10.1016/s2215-0366(21)00251-0
Lu, S., Peng, H., Wang, L., Vasish, S., Zhang, Y., Gao, W., et al. (2013). Elevated specific peripheral cytokines found in major depressive disorder patients with childhood trauma exposure: a cytokine antibody array analysis. Compr. Psychiatry 54, 953–961. doi: 10.1016/j.comppsych.2013.03.026
Luo, Z., Chen, G., Jia, Y., Zhong, S., Gong, J., Chen, F., et al. (2021). Shared and specific dynamics of brain segregation and integration in bipolar disorder and major depressive disorder: a resting-state functional magnetic resonance imaging study. J. Affect. Disord. 280, 279–286. doi: 10.1016/j.jad.2020.11.012
Maddock, R. J., Garrett, A. S., and Buonocore, M. H. (2001). Remembering familiar people: the posterior cingulate cortex and autobiographical memory retrieval. Neuroscience 104, 667–676. doi: 10.1016/s0306-4522(01)00108-7
Nock, M. K., Hwang, I., Sampson, N., Kessler, R. C., Angermeyer, M., Beautrais, A., et al. (2009). Cross-national analysis of the associations among mental disorders and suicidal behavior: findings from the WHO World Mental Health Surveys. PLoS Med. 6:e1000123. doi: 10.1371/journal.pmed.1000123
Patanaik, A., Tandi, J., Ong, J. L., Wang, C., Zhou, J., and Chee, M. W. L. (2018). Dynamic functional connectivity and its behavioral correlates beyond vigilance. Neuroimage 177, 1–10. doi: 10.1016/j.neuroimage.2018.04.049
Power, J. D., Barnes, K. A., Snyder, A. Z., Schlaggar, B. L., and Petersen, S. E. (2012). Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 59, 2142–2154. doi: 10.1016/j.neuroimage.2011.10.018
Qiao, D., Zhang, A., Sun, N., Yang, C., Li, J., Zhao, T., et al. (2020). Altered static and dynamic functional connectivity of habenula associated with suicidal ideation in first-episode. drug-naïve patients with major depressive disorder. Front. Psychiatry 11:608197. doi: 10.3389/fpsyt.2020.608197
Quigley, M., Cordes, D., Turski, P., Moritz, C., Haughton, V., Seth, R., et al. (2003). Role of the corpus callosum in functional connectivity. AJNR Am. J. Neuroradiol. 24, 208–212.
Ran, S., Zuo, Z., Li, C., Yin, X., Qu, W., Tang, Q., et al. (2020). Atrophic corpus callosum associated with altered functional asymmetry in major depressive disorder. Neuropsychiatr. Dis. Treat. 16, 1473–1482. doi: 10.2147/ndt.S245078
Richards, D. (2011). Prevalence and clinical course of depression: a review. Clin. Psychol. Rev. 31, 1117–1125. doi: 10.1016/j.cpr.2011.07.004
Ringo, J. L., Doty, R. W., Demeter, S., and Simard, P. Y. (1994). Time is of the essence: a conjecture that hemispheric specialization arises from interhemispheric conduction delay. Cereb. Cortex. 4, 331–343. doi: 10.1093/cercor/4.4.331
Sakoğlu, U., Pearlson, G. D., Kiehl, K. A., Wang, Y. M., Michael, A. M., and Calhoun, V. D. (2010). A method for evaluating dynamic functional network connectivity and task-modulation: application to schizophrenia. Magma 23, 351–366. doi: 10.1007/s10334-010-0197-8
Shi, Y., Li, J., Feng, Z., Xie, H., Duan, J., Chen, F., et al. (2020). Abnormal functional connectivity strength in first-episode, drug-naïve adult patients with major depressive disorder. Prog. Neuropsychopharmacol. Biol. Psychiatry 97:109759. doi: 10.1016/j.pnpbp.2019.109759
Song, P., Lin, H., Liu, C., Jiang, Y., Lin, Y., Xue, Q., et al. (2019). Transcranial magnetic stimulation to the middle frontal gyrus during attention modes induced dynamic module reconfiguration in brain networks. Front. Neuroinform. 13:22. doi: 10.3389/fninf.2019.00022
Tang, Q., Cui, Q., Chen, Y., Deng, J., Sheng, W., Yang, Y., et al. (2022). Shared and distinct changes in local dynamic functional connectivity patterns in major depressive and bipolar depressive disorders. J. Affect. Disord. 298, 43–50. doi: 10.1016/j.jad.2021.10.109
Tomasi, D., and Volkow, N. D. (2010). Functional connectivity density mapping. Proc. Natl. Acad. Sci. USA 107, 9885–9890. doi: 10.1073/pnas.1001414107
Wang, J., Wei, Q., Yuan, X., Jiang, X., Xu, J., Zhou, X., et al. (2018). Local functional connectivity density is closely associated with the response of electroconvulsive therapy in major depressive disorder. J. Affect. Disord. 225, 658–664. doi: 10.1016/j.jad.2017.09.001
Wang, L., Li, K., Zhang, Q. E., Zeng, Y. W., Jin, Z., Dai, W. J., et al. (2013). Interhemispheric functional connectivity and its relationships with clinical characteristics in major depressive disorder: a resting state fMRI study. PLoS One 8:e60191. doi: 10.1371/journal.pone.0060191
Wang, Y., Zhong, S., Jia, Y., Zhou, Z., Wang, B., Pan, J., et al. (2015). Interhemispheric resting state functional connectivity abnormalities in unipolar depression and bipolar depression. Bipolar. Disord. 17, 486–495. doi: 10.1111/bdi.12315
Wittling, W., and Roschmann, R. (1993). Emotion-related hemisphere asymmetry: subjective emotional responses to laterally presented films. Cortex 29, 431–448. doi: 10.1016/s0010-9452(13)80252-3
Yang, G. J., Murray, J. D., Repovs, G., Cole, M. W., Savic, A., Glasser, M. F., et al. (2014). Altered global brain signal in schizophrenia. Proc. Natl. Acad. Sci. USA 111, 7438–7443. doi: 10.1073/pnas.1405289111
Yang, Y., Cui, Q., Pang, Y., Chen, Y., Tang, Q., Guo, X., et al. (2021). Frequency-specific alteration of functional connectivity density in bipolar disorder depression. Prog. Neuropsychopharmacol. Biol. Psychiatry 104:110026. doi: 10.1016/j.pnpbp.2020.110026
Yüksel, D., Dietsche, B., Konrad, C., Dannlowski, U., Kircher, T., and Krug, A. (2018). Neural correlates of working memory in first episode and recurrent depression: an fMRI study. Prog. Neuropsychopharmacol. Biol. Psychiatry 84, 39–49. doi: 10.1016/j.pnpbp.2018.02.003
Zhang, H., Palaniyappan, L., Wu, Y., Cong, E., Wu, C., Ding, L., et al. (2020). The concurrent disturbance of dynamic functional and structural brain connectome in major depressive disorder: the prefronto-insular pathway. J. Affect. Disord. 274, 1084–1090. doi: 10.1016/j.jad.2020.05.148
Zhang, J., Magioncalda, P., Huang, Z., Tan, Z., Hu, X., Hu, Z., et al. (2019). Altered global signal topography and its different regional localization in motor cortex and hippocampus in mania and depression. Schizophr Bull. 45, 902–910. doi: 10.1093/schbul/sby138
Zheng, R., Chen, Y., Jiang, Y., Wen, M., Zhou, B., Li, S., et al. (2021). Dynamic altered amplitude of low-frequency fluctuations in patients with major depressive disorder. Front. Psychiatry 12:683610. doi: 10.3389/fpsyt.2021.683610
Zhu, D. M., Yang, Y., Zhang, Y., Wang, C., Wang, Y., Zhang, C., et al. (2020). Cerebellar-cerebral dynamic functional connectivity alterations in major depressive disorder. J. Affect. Disord. 275, 319–328. doi: 10.1016/j.jad.2020.06.062
Zhu, F., Liu, F., Guo, W., Chen, J., Su, Q., Zhang, Z., et al. (2018). Disrupted asymmetry of inter- and intra-hemispheric functional connectivity in patients with drug-naive, first-episode schizophrenia and their unaffected siblings. EBioMedicine 36, 429–435. doi: 10.1016/j.ebiom.2018.09.012
Zhuo, C., Zhou, C., Lin, X., Tian, H., Wang, L., Chen, C., et al. (2020). Common and distinct global functional connectivity density alterations in drug-naïve patients with first-episode major depressive disorder with and without auditory verbal hallucination. Prog. Neuropsychopharmacol. Biol. Psychiatry 96:109738. doi: 10.1016/j.pnpbp.2019.109738
Zhuo, C., Zhu, J., Wang, C., Qu, H., Ma, X., and Qin, W. (2017). Different spatial patterns of brain atrophy and global functional connectivity impairments in major depressive disorder. Brain Imag. Behav. 11, 1678–1689. doi: 10.1007/s11682-016-9645-z
Zou, K., Gao, Q., Long, Z., Xu, F., Sun, X., Chen, H., et al. (2016). Abnormal functional connectivity density in first-episode, drug-naive adult patients with major depressive disorder. J. Affect. Disord. 194, 153–158. doi: 10.1016/j.jad.2015.12.081
Zuo, X. N., Ehmke, R., Mennes, M., Imperati, D., Castellanos, F. X., Sporns, O., et al. (2012). Network centrality in the human functional connectome. Cereb. Cortex. 22, 1862–1875. doi: 10.1093/cercor/bhr269
Keywords: major depressive disorder, resting-state functional magnetic resonance imaging, interhemisphere, intrahemisphere, dynamic functional connectivity, static functional connectivity
Citation: Jiang Y, Chen Y, Zheng R, Zhou B, Wei Y, Gao A, Wei Y, Li S, Guo J, Han S, Zhang Y and Cheng J (2022) More Than Just Statics: Temporal Dynamic Changes in Inter- and Intrahemispheric Functional Connectivity in First-Episode, Drug-Naive Patients With Major Depressive Disorder. Front. Hum. Neurosci. 16:868135. doi: 10.3389/fnhum.2022.868135
Received: 02 February 2022; Accepted: 14 March 2022;
Published: 08 April 2022.
Edited by:
Valeria Mondini, Graz University of Technology, AustriaReviewed by:
Jiajia Zhu, First Affiliated Hospital of Anhui Medical University, ChinaPeng Wang, Chinese Academy of Medical Sciences and Peking Union Medical College, China
Copyright © 2022 Jiang, Chen, Zheng, Zhou, Wei, Gao, Wei, Li, Guo, Han, Zhang and Cheng. 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: Shaoqiang Han, ShaoqiangHan@163.com; Yong Zhang, zzuzhangyong2013@163.com; Jingliang Cheng, fccchengjl@zzu.edu.cn
†These authors have contributed equally to this work and share first authorship