Altered Brain Function in First-Episode and Recurrent Depression: A Resting-State Functional Magnetic Resonance Imaging Study

Background: Studies on differences in brain function activity between the first depressive episode (FDE) and recurrent depressive episodes (RDE) are scarce. In this study, we used regional homogeneity (ReHo) and amplitude of low-frequency fluctuations (ALFF) as indices of abnormal brain function activity. We aimed to determine the differences in these indices between patients with FDE and those with RDE, and to investigate the correlation between areas of abnormal brain function and clinical symptoms.

Methods: A total of 29 patients with RDE, 28 patients with FDE, and 29 healthy controls (HCs) who underwent resting-state functional magnetic resonance imaging were included in this study. The ReHo and ALFF measurements were used for image analysis and further analysis of the correlation between different brain regions and clinical symptoms.

Results: Analysis of variance showed significant differences among the three groups in ReHo and ALFF in the frontal, parietal, temporal, and occipital lobes. ReHo was higher in the right inferior frontal triangular gyrus and lower in the left inferior temporal gyrus in the RDE group than in the FDE group. Meanwhile, ALFF was higher in the right inferior frontal triangular gyrus, left anterior cingulate gyrus, orbital part of the left middle frontal gyrus, orbital part of the left superior frontal gyrus, and right angular gyrus, but was lower in the right lingual gyrus in the RDE group than in the FDE group. ReHo and ALFF were lower in the left angular gyrus in the RDE and FDE groups than in the HC group. Pearson correlation analysis showed a positive correlation between the ReHo and ALFF values in these abnormal areas in the frontal lobe and the severity of depressive symptoms (P < 0.05). Abnormal areas in the temporal and occipital lobes were negatively correlated with the severity of depressive symptoms (P < 0.05).

Conclusion: The RDE and FDE groups had abnormal neural function activity in some of the same brain regions. ReHo and ALFF were more widely distributed in different brain regions and had more complex neuropathological mechanisms in the RDE group than in the FDE group, especially in the right inferior frontal triangular gyrus of the frontal lobe.

Introduction

Major depressive disorder (MDD) is a common clinical psychiatric disorder characterized by affective, cognitive, and somatic symptoms. Most of its clinical manifestations are common symptoms such as depressed mood, diminished interest, slowed thinking, loss of appetite, and insomnia (Wolfers et al., 2015). The World Health Organization reports that an estimated 350 million people worldwide have depression and that depression is the leading cause of suicide, with approximately 800,000 people dying by suicide each year (Zhang et al., 2018). By 2030, depression is expected to become the disease with the largest burden worldwide (Friedrich, 2017). Currently, the diagnosis of MDD relies on clinical scale assessment of the patient and the experience of the psychiatrist; however, knowledge about precise neurobiological biomarkers is lacking.

According to the classification criteria of the 10th revision of the International Statistical Classification of Diseases and Related Health Problems, MDD is classified into the first depressive episode (FDE) and recurrent depressive episodes (RDE) (Yuksel et al., 2018). Previous studies have demonstrated differences in depressive and somatic symptoms (Roca et al., 2011), cognitive functioning (Roca et al., 2015; Zu et al., 2021; Varghese et al., 2022), and quality of life (Zu et al., 2021) between patients with RDE and those with FDE. Epidemiological surveys have shown that depression has a high recurrence rate and that, once the first episode has occurred, it relapses within 5 years. In addition, approximately 80% of patients with depression have a history of two recurrences (Burcusa and Iacono, 2007). Moreover, the severity of MDD increases with the number of relapses (Harlev et al., 2021), and it is increasingly being recognized that the challenge in patients with depression is preventing relapse rather than promoting recovery (Fava et al., 2017). Therefore, the differences in pathogenesis between RDE and FDE need to be elucidated from a neuropathological perspective.

In recent years, with rapid developments in neuroimaging technology, resting-state functional magnetic resonance imaging (rs-fMRI) has gradually been applied to the study of insomnia (Marques et al., 2018), schizophrenia (Sheffield and Barch, 2016), autism (Bathelt and Geurts, 2021), and other psychiatric systemic disorders. It has also been applied to the study of MDD subtypes (Brown et al., 2019; Yue et al., 2020). rs-fMRI indirectly reflects neuronal spontaneous activity by measuring the blood oxygen level-dependent (BOLD) signal, which enables the investigation of brain function abnormalities in the early disease stage and has advantages of simplicity, non-invasiveness, and reproducibility (Lee et al., 2013). Regional homogeneity (ReHo) and amplitude of low-frequency fluctuations (ALFF) are two of the commonly used metrics in rs-fMRI analysis. ReHo mainly assesses the synchronous reflection of the sequence between a given voxel and its neighboring voxels, reflecting the temporal homogeneity of signals related to regional blood oxygen levels (Zang et al., 2004). ALFF is the sum of the spectral amplitude of each voxel signal in the low-frequency range (usually 0.01–0.08 Hz), reflecting the amplitude of low-frequency fluctuations caused by spontaneous neuronal activity (Zang et al., 2007). Combining the ReHo and ALFF methods can contribute to a better understanding of the abnormalities of brain function in patients with MDD (Ni et al., 2016).

However, the combination of ALFF and ReHo has been mainly applied in studies on mild cognitive impairment (Ni et al., 2016), schizophrenia (Zhao et al., 2018), and anxiety disorders (Shen et al., 2020), whereas the differences between FDE and RDE have not been investigated. A study in patients with FDE and remitting MDD (rMDD) showed differences in ALFF and ReHo in the temporal lobes, although patients with rMDD had been cured with antidepressant medication (Yang et al., 2018). Another study observed alterations in bilateral frontal BOLD signals in patients with FDE or RDE; however, the study used a task-state observation approach (Yuksel et al., 2018). The frontal lobe is an important brain region in the pathogenesis of MDD, and previous studies have also found differences in frontal lobe function between FDE and RDE (Talarowska et al., 2015).

In this study, we aimed to determine the differences in brain function activity among patients with FDE, patients with RDE, and healthy controls (HCs) by using the ReHo and ALFF methods based on rs-fMRI techniques. Further, we also aimed to analyze the correlation between different brain areas and clinical characteristics. We hypothesized that differences would be observed in the brain neural circuits of patients with FDE and those with RDE, especially in areas closely related to the frontal lobe. This study will provide a neuroimaging basis for the differences in neuropathological mechanisms between FDE and RDE and some insights for clinical research.

Materials and Methods

Participants

A total of 59 outpatients with MDD from Guang’anmen Hospital, Chinese Academy of Traditional Chinese Medicine, Beijing First Hospital of Integrative Medicine, were recruited for this study. All patients met the Diagnostic and Statistical Manual of Mental Disorders Fifth Edition criteria for MDD. We used the 17-item Hamilton Rating Scale for Depression (HAMD-17) (Hamilton, 1960) to assess the severity of depression in all patients and classified these patients into those with RDE [n = 31; mean frequency of recurrence, 2.27 (standard deviation = 0.79)] and those with FDE (n = 28, 0 recurrence). The inclusion criteria were (1) age 18–55 years; (2) HAMD-17 score > 17; and (3) The FDE group all had their FDE prior to enrollment and were not receiving any antidepressant medication. The RDE group had a previous history of depression, cured by antidepressant medication, now recurring and a history of antidepressant withdrawal for at least 4 weeks prior to enrollment. We also included 29 sex- and age-matched HCs (21 women and 8 men) who (1) were aged 18–55 years, (2) had a HAMD–17 score of < 7, (3) had right-handedness, and (4) had no history of any mental illness in first-degree relatives.

The exclusion criteria for patients and HCs were as follows: (1) serious mental illness and other diseases such as cardiovascular and cerebrovascular disorders; (2) history of drug and alcohol abuse; (3) any contraindications to MRI, such as presence of a heart pacemaker, metal fixed false teeth, or severe claustrophobia; (4) pregnant or lactating status; and (5) bipolar disorder or suicidal ideation.

All patients were required to sign an informed consent form before enrollment. This study was approved by the ethics committee of Guang’anmen Hospital, Chinese Academy of Traditional Chinese Medicine.

Scan Acquisition

All patients in this study underwent MRI using a Magnetom Skyra 3.0-T scanner (Siemens, Erlangen, Germany). Before the scanning procedure, the patients were instructed to remain awake and avoid active thinking. During the scanning process, the patients were required to wear earplugs and noise-canceling headphones, to use a hood to immobilize the head, and to lie flat on the examination bed. The scanning procedure involved a localizer scan, high-resolution three-dimensional T1-weighted imaging, and BOLD-fMRI.

The scanning parameters were as follows: for three-dimensional T1-weighted imaging, time repetition/time echo = 2,500/2.98 ms, flip angle = 7°, matrix = 64 × 64, field of view = 256 mm × 256 mm, slice thickness = 1 mm, slice number = 48, slices = 192, scanning time = 6 min 3 s; for BOLD-fMRI, time repetition/time echo = 2,000/30 ms, flip angle = 90°, matrix = 64 × 64, field of view = 240 mm × 240 mm, slice number = 43, slice thickness/spacing = 3.0/1.0 mm, number of obtained volumes = 200, and scanning time = 6 min 40 s.

Image Processing

Functional Magnetic Resonance Imaging Data Preprocessing

The rs-fMRI data were preprocessed using DPARSF (Data Processing Assistant for rs-fMRI) software (DPABI5.0)¹ (Yan and Zang, 2010) in MATLAB (Mathworks Inc., Natick, MA, United States), according to the following process: (1) conversion of DICOM raw data to NIFTI format; (2) removal of the first 10 time points to place the data in a stable state; (3) slice timing; (4) realignment of head motion (removal of patients with head movements > 2 mm in any direction and motor rotation > 2°); (5) regression of covariates, including brain white matter signal, cerebrospinal fluid signal, and head movement parameters; (6) spatial normalization (the functional images of all patients were converted to Montreal Neurological Institute standard space using the DARTEL method); and (7) linear detrending and filtering (0.01–0.08 Hz).

Regional Homogeneity Analysis

The DPARSF (Data Processing Assistant for Resting-State fMRI Advanced Edition_V5.1_201001, see text footnote 1) software was used to analyze the ReHo and ALFF of the pre-processed data. The similarity of the time series of each voxel to its neighboring voxels (26 neighboring voxels) was assessed using the Kendall’s coefficient of concordance (KCC) (Kendall, 1990), i.e., ReHo values. The whole-brain ReHo images of the subjects were obtained by calculating the KCC values of the whole-brain voxels. To improve the signal-to-noise ratio, the Re Ho images were spatially smoothed using a 6 mm × 6 mm × 6 mm full-width half-height Gaussian kernel.

Amplitude of Low-Frequency Fluctuations Analysis

Data were spatially normalized and smoothed, and a fast Fourier transform was performed to switch the time series to the frequency domain to obtain the power spectrum. The square root of the power spectrum at each frequency was calculated to obtain the average square root of the ALFF measurement for each voxel in the range of 0.01–0.08 Hz. Finally, time bandpass filtering (0.01–0.08 Hz) was performed. To reduce the inter-individual variability, ALFF was transformed to zALFF using Fisher’s z transformation before statistical analysis.

Statistical Analyses

Clinical Data Analysis

Clinical data were analyzed using the SPSS 23.0 statistical software (IBM Corporation, Somers, NY, United States). One-way analysis of variance was used to compare age and educational level among the three groups, and the chi-square test was used to compare sex differences. A two-sample t-test was used to compare the duration of disease, HAMD-17 scores, and frequency of recurrence between the two patient groups, with a threshold of P < 0.05 (two-tailed) set as statistically significant.

Functional Magnetic Resonance Imaging Data Analysis

Imaging data were analyzed using the DPARSF toolbox, and a voxel-based one-way analysis of variance was performed to compare the whole-brain ReHo/ALFF map among the three groups. Sex, age, educational level, and framewise displacement (a metric derived from Jenkinson’s formula) were used as covariates, and brain areas with ReHo/ALFF differences among the three groups were corrected for Gaussian random fields. The corrected cluster level was set at P < 0.05 (two-tailed), and threshold voxel levels of P < 0.005 were defined as statistically different. We performed post-hoc t-test analysis using DPARSF 5.1 software for two-by-two comparisons between groups, and Bonferroni correction was applied to the results, setting a threshold of P < 0.016 (0.05/3) for statistical significance. The threshold was set to clusters > 10 voxels.

To verify the relationship between ReHo/ALFF values and clinical symptoms, we extracted the mean ReHo/ALFF values of three different brain regions and performed Pearson correlation analysis of the clinical scale scores of each group. Significance was set at a statistical threshold of P < 0.05 (two-tailed).

Results

Characteristics of Research Samples

Two patients with RDE were excluded because of excessive head movement displacement. Therefore, a total of 29 patients with RDE, 28 patients with FDE, and 29 HCs met the inclusion criteria. No statistical differences among the three groups were found in terms of sex, age, and years of education. The HAMD-17 scores were not statistically different between the RDE and FDE groups, whereas a statistical difference was observed in the duration of illness and frequency of recurrence (Table 1).

TABLE 1

Table 1. Demographic and clinical characteristics of the study participants.

Differences in Regional Homogeneity/Amplitude of Low-Frequency Fluctuations Between the Recurrent Depressive Episodes, First Depressive Episode, and Healthy Controls

One-way analysis of variance showed significant differences in ReHo and ALFF among the three groups in the right inferior frontal triangular gyrus, left anterior cingulate cortex, left middle temporal gyrus/left angular gyrus, and left inferior temporal gyrus/left inferior occipital gyrus. Meanwhile, ReHo was significantly different in the left superior temporal gyrus and right Rolandic operculum gyrus. ALFF was also significantly different in the right angular gyrus, right lingual gyrus, and Vermis_3 (Tables 2, 3 and Figures 1, 2).

TABLE 2

Table 2. ReHo differences in RDE, FDE, and HCs.

TABLE 3

Table 3. ALFF differences in RDE, FDE, and HCs.

FIGURE 1

Figure 1. Brain regions with abnormal ReHo among the three groups based on one-way analysis of variance. The color bars indicate the F-value.

FIGURE 2

Figure 2. Brain regions with abnormal ALFF among the three groups based on one-way analysis of variance. The color bars indicate the F-value.

Compared with the FDE group, the RDE group had elevated ReHo in the right inferior frontal triangular gyrus and decreased ReHo in the left inferior temporal gyrus. Meanwhile, the RDE group had elevated ALFF in the right inferior frontal triangular gyrus, left anterior cingulate gyrus, orbital part of the left middle frontal gyrus, orbital part of the left superior frontal gyrus, and right angular gyrus, but had decreased ALFF in the right lingual gyrus (Tables 2, 3 and Figure 3).

FIGURE 3

Figure 3. Brain regions with abnormal ReHo (left) and ALFF (right) between the RDE group and the FDE group based on a post hoc t-test. The color bars indicate the T-value.

Compared with the HC group, the RDE group had elevated ReHo in the right inferior frontal triangular gyrus, left anterior cingulate gyrus, and left superior temporal gyrus, but had reduced ReHo in the left angular gyrus. Meanwhile, the RDE group had elevated ALFF in the left anterior cingulate gyrus and right angular gyrus, but had decreased ALFF in the left angular gyrus (Tables 2, 3 and Figure 4).

FIGURE 4

Figure 4. Brain regions with abnormal ReHo (left) and ALFF (right) between the RDE group and the HC group based on a post hoc t-test. The color bars indicate the T-value.

Compared with the HC group, the FDE group had elevated ReHo in the left inferior temporal gyrus and decreased ReHo in the left middle temporal gyrus/left angular gyrus. Meanwhile, the FDE group had elevated ALFF in the left inferior temporal gyrus and Vermis_3, but had decreased ALFF in the left middle temporal gyrus/left angular gyrus (Tables 2, 3 and Figure 5).

FIGURE 5

Figure 5. Brain regions with abnormal ReHo (left) and ALFF (right) between the FDE group and the HC group based on a post hoc t-test. The color bars indicate the T-value.

Significant Correlation Between Functional Image and Clinical Feature

To test the correlation between areas of abnormal brain activity and the severity of clinical depressive symptoms, we further performed Pearson correlation analysis. We found a positive correlation between the ReHo/ALFF values of the right inferior frontal triangular gyrus and the HAMD-17 scores in the RDE group (r = 0.436, P = 0.018; r = 0.394, P = 0.034). Meanwhile, we observed a positive correlation between the ReHo values in the left anterior cingulate cortex/orbital part of the left middle frontal gyrus and the HAMD-17 scores in the FDE group (r = 0.488, P = 0.008). Furthermore, we found a negative correlation between the ReHo values in the left inferior temporal gyrus/left inferior occipital gyrus and right lingual gyrus and the HAMD-17 scores in the FDE group (r = −0.412, P = 0.029; r = −0.408, P = 0.030) (Table 4 and Figure 6).

TABLE 4

Table 4. Correlation of abnormal brain areas with clinical symptoms.

FIGURE 6

Figure 6. Positive correlation between the ReHo/ALFF values of abnormal brain regions and the HAMD-17 scores. (A) ReHo values in the recurrent depressive episode (RDE) group. (C) ALFF values in the RDE group. (D) ALFF values in the first depressive episode (FDE) group. Negative correlation between the ReHo/ALFF values of abnormal brain regions and the HAMD-17 scores. (B) ReHo values in the FDE group. (E) ALFF values in the FDE group. Frontal_Inf_Tri_R, right inferior frontal triangular gyrus; Cingulum_Ant_L/Frontal_Med_Orb_L, left anterior cingulate cortex/orbital part of the left middle frontal gyrus; Temporal_Inf_L/Occipital_Inf_L, left inferior temporal gyrus/left inferior occipital gyrus; Lingual_R, right lingual gyrus; ReHo, regional homogeneity; ALFF, amplitude of low-frequency fluctuations; HAMD-17, 17-item Hamilton Rating Scale for Depression.

Discussion

To our knowledge, this is the first study to investigate the changes in local brain function activity in patients with RDE or FDE using the ReHo and ALFF methods. The results of this study showed that the RDE and FDE groups had abnormal neural function activity in some of the same brain regions. ReHo and ALFF were more widely distributed in different brain regions and had more complex neuropathological mechanisms in the RDE group than in the FDE group. This study provides a reference for the differences in brain function activity between RDE and FDE.

We found increased ReHo and ALFF in the right inferior frontal triangular gyrus in the RDE group compared with the FDE group. The right inferior frontal triangular gyrus is located in the dorsolateral prefrontal cortex (DLPFC), an important component of executive function that is closely related to working memory, thought activity, and cognitive control (Wang et al., 2020; Brosch et al., 2021; Nejati et al., 2022). Previous studies have also shown that depressed mood and cognitive behavioral impairment in patients with depression are associated with abnormalities in executive function (Sun et al., 2018). Patients with mild impairment in the DLPFC often show symptoms of depression, such as loss of interest, poor memory, slow thinking, lack of motivation, and insomnia (Elliott, 2003; Li et al., 2021). Previous studies have found that the right inferior frontal triangle gyrus plays an important role in the early pathogenesis of MDD (Zhang et al., 2021) and that the gray matter volume of the DLPFC is reduced in patients with RDE after 6 weeks of antidepressant therapy (Li et al., 2010). Another study showed that stimulation of the right and left DLPFC with transcranial direct current stimulation was effective in reducing the risk of recurrence in patients with depression (Aparicio et al., 2019). The BOLD signals in the DLPFC have also been found to differ between RDE and FDE (Yuksel et al., 2018). This suggests that patients with RDE and those with FDE have differences in executive function. The correlation analysis in this study showed that the ReHo and ALFF values in the right inferior frontal triangular gyrus were positively correlated with the HAMD-17 scores in the RDE group, whereas no such correlation was found in the FDE group. These results suggest that the right inferior frontal triangular gyrus is an important neurobiological imaging marker for RDE and an important brain region for differentiating RDE from FDE.

We also found differences in ALFF in the left anterior cingulate gyrus and left orbitofrontal gyrus between the RDE and FDE groups. The anterior cingulate gyrus is an important component of the limbic system. It has extensive fibrous connections to many cortical and subcortical structures; is involved in the regulation of a wide range of functions, such as emotion, cognition, and motivation; and is closely associated with the onset of depression (Xiao and Zhang, 2018; Zheng et al., 2018; Rolls, 2019). Previous studies have found higher ALFF in the right ventral anterior cingulate gyrus in patients with anxious depression than in patients with rMDD and in HCs. In addition, abnormal activation of the anterior cingulate gyrus at rest in patients with MDD may be related to the failure of emotional control, which is a central factor in negative rumination and persistent self-focus in these patients (Amodio and Frith, 2006; Liu et al., 2015). Therefore, the results of this study suggest that patients with RDE show abnormal activation of the left side of the anterior cingulate gyrus. The orbitofrontal cortex is an important part of the reward network and is closely associated with emotional information and sensory stimuli (O’Doherty, 2004; Liu et al., 2017). A previous study found that patients with refractory depression, which is characterized by persistence and recurrence, had significantly higher ALFF in the orbitofrontal cortex than HCs (Liu et al., 2014). This suggests that abnormal activation of the orbitofrontal cortex is a cause of the complexity of RDE. Another study also showed that patients with FDE had reduced ALFF in the left and right orbitofrontal cortex compared with HCs, suggesting that patients with FDE have reduced regulation of the reward network (Zhang et al., 2014). Although these results were from different studies, they all support the involvement of the anterior cingulate gyrus and orbitofrontal cortex in regulating the pathophysiology of MDD. The correlation analysis in the current study showed that the ALFF values in the left anterior cingulate cortex/orbital part of the left middle frontal gyrus were positively correlated with the HAMD-17 scores in the FDE group, thus identifying an important brain region for differentiating between RDE and FDE.

In addition, we found that the RDE group had higher ALFF values in the right lingual gyrus than the FDE group. The lingual gyrus is part of the occipital lobe, and previous studies have shown that the lingual gyrus is also involved in activities related to visual memory processing and is closely associated with the development of MDD (Jung et al., 2014; Le et al., 2017; Palejwala et al., 2021). One study found that patients with rMDD had lower ALFF values in the right lingual gyrus than HCs, suggesting that reduced ALFF in the right lingual gyrus is a marker for the remission of depression (Yang et al., 2018). Another study also showed that patients with FDE had lower ALFF values in the right lingual gyrus than HCs (Wang et al., 2012). Our correlation analysis showed that the ALFF values in the right lingual gyrus were negatively correlated with the HAMD-17 scores in the FDE group. This suggests that the right lingual gyrus can be used as a status marker for FDE.

Meanwhile, we found that the ALFF in the right angular gyrus was higher in the RDE group than in the FDE group. A previous study reported that the right angular gyrus is associated with the self-localization function in humans and is potentially associated with psychiatric disorders, and that an abnormality in this function affects both sensory and perceptual functions, with a causal relationship between them (de Boer et al., 2020). This may also be informative in elucidating the differences in neuropathological mechanisms between RDE and FDE. The inferior temporal gyrus is involved in functions such as social cognition, emotional stimulus processing, self-referential processing, and semantic processing and is closely associated with MDD development (Chao et al., 1999; Cabeza and Nyberg, 2000; Herath et al., 2001; Hu et al., 2017; Kocsis et al., 2021). Previous studies have found that ReHo is significantly higher in the left middle temporal gyrus in patients with FDE than in those with rMDD, suggesting that enhanced metabolism in the left middle temporal gyrus may be one of the pathogenic mechanisms of FDE (Yang et al., 2018). Therefore, we suggest that differences exist between RDE and FDE in terms of abnormalities in the frontal, temporal, parietal, and occipital lobes, particularly in the right inferior frontal triangular gyrus of the frontal lobe.

Interestingly, we also found that the RDE and FDE groups had lower ReHo and ALFF in the left angular gyrus than the HC group. The angular gyrus is located in the posterior part of the inferior parietal lobe and is mainly involved in human semantic and numerical processing, memory retrieval, spatial cognition, word reading and comprehension, reasoning, and social cognition, with the left angular gyrus playing a more important role in situational simulation and memory (Seghier, 2013; Thakral et al., 2017; Ramanan and Bellana, 2019). The angular gyrus is also an important part of the default mode network (Raichle et al., 2001; Raichle, 2015). Previous studies have found differences in BOLD signals in the angular gyrus between patients with RDE or FDE and HCs, consistent with the present study (Yuksel et al., 2018). Other studies have also reported that the ReHo and ALFF values in the left angular gyrus were lower in patients with MDD than in HCs (Wang et al., 2012; Liu et al., 2021). Abnormalities in angular gyrus function can lead to cognitive impairment, a common clinical manifestation of MDD (Lee et al., 2012). Therefore, this suggests that (1) both RDE and FDE are characterized by default mode network dysfunction and (2) ReHo and ALFF abnormalities in the left angular gyrus are important markers for differentiating patients with MDD from HCs.

Some limitations of this study should be considered. First, the patients with RDE enrolled in this study might have still been affected by an underlying antidepressant action despite having stopped medication for 4 weeks. Second, the number of recurrences in the RDE group was inconsistent, and a first-recurrence study of patients with RDE seem to have greater research value. Third, this study focused on only one scale (HAMD-17). To enhance the scientific value of this study, more scales need to be used in the future to focus on the detailed correlation of cognitive, somatic, anxiety, and insomnia symptoms with RDE and FDE. Finally, this study did not find a difference between RDE and FDE in depressive symptoms, which may be related to the specificity and small size of the study population. Further studies with a larger sample size are needed to confirm or refute the findings of this study.

Conclusion

This study used ReHo and ALFF, which are indices based on the rs-fMRI technique, to preliminarily analyze the differences between RDE and FDE in terms of neural activity in different brain regions. The RDE and FDE groups showed abnormal changes in neural function activity in some of the same brain regions, with ReHo and ALFF being more widely distributed in different brain regions and the neuropathological mechanisms being more complex in the RDE group than in the FDE group, particularly in the right inferior frontal triangular gyrus of the frontal lobe.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The experimental protocol was approved by the Ethics Committee of Guang’anmen Hospital, Chinese Academy of Traditional Chinese Medicine (No. 2017-021-SQ), Trial registration: China Clinical Trials Registry, chiCTR1800014277. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

JF conceived and designed this experiment. JS wrote and revised the manuscript and participated in the collection of cases and statistical analysis of the data. LC drew diagrams and made statistical analysis of data and revised the manuscript. ZD analyzed the data and revised the manuscript. JH, DG, YM, ZW, CG, and YL participated in data analysis, case collection, and manuscript writing. YH, LZ, and JC performed fMRI on the patient. FX, XH, XX, JT, and XY involved in case collection and symptom assessment of patients. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by the Science and Technology Innovation Project of Chinese Academy of Traditional Chinese Medicine (CI2021A03301), the National Natural Science Foundation of China (82174282 and 81774433), and the National Key Research and Development Program of China (2018YFC1705800).

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.

Acknowledgments

We thank all the patients, healthy people who participated in the experiment for their support.

Footnotes

1. ^ http://www.rfmri.org/DPARSF

References

Amodio, D. M., and Frith, C. D. (2006). Meeting of minds: the medial frontal cortex and social cognition. Nat. Rev. Neurosci. 7, 268–277. doi: 10.1038/nrn1884

PubMed Abstract | CrossRef Full Text | Google Scholar

Aparicio, L. V. M., Rosa, V., Razza, L. M., Sampaio-Junior, B., Borrione, L., Valiengo, L., et al. (2019). Transcranial direct current stimulation (tDCS) for preventing major depressive disorder relapse: Results of a 6-month follow-up. Depress Anxiety 36, 262–268. doi: 10.1002/da.22878

PubMed Abstract | CrossRef Full Text | Google Scholar

Bathelt, J., and Geurts, H. M. (2021). Difference in default mode network subsystems in autism across childhood and adolescence. Autism 25, 556–565. doi: 10.1177/1362361320969258

PubMed Abstract | CrossRef Full Text | Google Scholar

Brosch, K., Stein, F., Meller, T., Schmitt, S., Yuksel, D., Ringwald, K. G., et al. (2021). DLPFC volume is a neural correlate of resilience in healthy high-risk individuals with both childhood maltreatment and familial risk for depression. Psychol. Med. 1–7. doi: 10.1017/S0033291721001094 [Epub online ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, E., Clark, D., Hassel, S., MacQueen, G., and Ramasubbu, R. J. N. d (2019). Intrinsic thalamocortical connectivity varies in the age of onset subtypes in major depressive disorder. Neuropsychiatr. Dis. Treat. 15, 75–82. doi: 10.2147/ndt.S184425

PubMed Abstract | CrossRef Full Text | Google Scholar

Burcusa, S. L., and Iacono, W. G. (2007). Risk for recurrence in depression. Clin. Psychol. Rev. 27, 959–985. doi: 10.1016/j.cpr.2007.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Cabeza, R., and Nyberg, L. (2000). Imaging cognition II: An empirical review of 275 PET and fMRI studies. J. Cogn. Neurosci. 12, 1–47. doi: 10.1162/08989290051137585

PubMed Abstract | CrossRef Full Text | Google Scholar

Chao, L. L., Haxby, J. V., and Martin, A. (1999). Attribute-based neural substrates in temporal cortex for perceiving and knowing about objects. Nat. Neurosci. 2, 913–919. doi: 10.1038/13217

PubMed Abstract | CrossRef Full Text | Google Scholar

de Boer, D. M. L., Johnston, P. J., Kerr, G., Meinzer, M., and Cleeremans, A. (2020). A causal role for the right angular gyrus in self-location mediated perspective taking. Sci. Rep. 10:19229. doi: 10.1038/s41598-020-76235-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Elliott, R. (2003). Executive functions and their disorders. Br. Med. Bull. 65, 49–59. doi: 10.1093/bmb/65.1.49

PubMed Abstract | CrossRef Full Text | Google Scholar

Fava, G. A., Cosci, F., Guidi, J., and Tomba, E. (2017). Well-being therapy in depression: New insights into the role of psychological well-being in the clinical process. Depress Anxiety 34, 801–808. doi: 10.1002/da.22629

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedrich, M. J. (2017). Depression Is the Leading Cause of Disability Around the World. JAMA 317:1517. doi: 10.1001/jama.2017.3826

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamilton, M. (1960). A rating scale for depression. J. Neurol. Neurosurg. Psychiatr. 23, 56–62. doi: 10.1136/jnnp.23.1.56

PubMed Abstract | CrossRef Full Text | Google Scholar

Harlev, D., Ravona-Springer, R., Nuriel, Y., and Fruchter, E. (2021). Sleep Monitoring Using WatchPAT Device to Predict Recurrence of Major Depression in Patients at High Risk for Major Depression Disorder Recurrence: A Case Report. Front. Psychiatr. 12:572660. doi: 10.3389/fpsyt.2021.572660

PubMed Abstract | CrossRef Full Text | Google Scholar

Herath, P., Kinomura, S., and Roland, P. E. (2001). Visual recognition: evidence for two distinctive mechanisms from a PET study. Hum. Brain Mapp. 12, 110–119. doi: 10.1002/1097-0193(200102)12:2<110::aid-hbm1008<3.0.co;2-0

CrossRef Full Text | Google Scholar

Hu, B., Rao, J., Li, X., Cao, T., Li, J., Majoe, D., et al. (2017). Emotion Regulating Attentional Control Abnormalities In Major Depressive Disorder: An Event-Related Potential Study. Sci. Rep. 7:13530. doi: 10.1038/s41598-017-13626-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Jung, J., Kang, J., Won, E., Nam, K., Lee, M. S., Tae, W. S., et al. (2014). Impact of lingual gyrus volume on antidepressant response and neurocognitive functions in Major Depressive Disorder: a voxel-based morphometry study. J. Affect. Disord. 169, 179–187. doi: 10.1016/j.jad.2014.08.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Kendall, M. G. (1990). Rank Correlation Methods (5th ed ed.). Oxford: Oxford University Press.

Google Scholar

Kocsis, K., Holczer, A., Kazinczi, C., Boross, K., Horvath, R., Nemeth, L. V., et al. (2021). Voxel-based asymmetry of the regional gray matter over the inferior temporal gyrus correlates with depressive symptoms in medicated patients with major depressive disorder. Psychiatr. Res. Neuroimaging 317:111378. doi: 10.1016/j.pscychresns.2021.111378

PubMed Abstract | CrossRef Full Text | Google Scholar

Le, T. M., Borghi, J. A., Kujawa, A. J., Klein, D. N., and Leung, H. C. (2017). Alterations in visual cortical activation and connectivity with prefrontal cortex during working memory updating in major depressive disorder. Neuroimage Clin. 14, 43–53. doi: 10.1016/j.nicl.2017.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, M. H., Smyser, C. D., and Shimony, J. S. (2013). Resting-state fMRI: a review of methods and clinical applications. Am. J. Neuroradiol. 34, 1866–1872. doi: 10.3174/ajnr.A3263

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, R. S., Hermens, D. F., Porter, M. A., and Redoblado-Hodge, M. A. (2012). A meta-analysis of cognitive deficits in first-episode Major Depressive Disorder. J. Affect. Disord. 140, 113–124. doi: 10.1016/j.jad.2011.10.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C. T., Lin, C. P., Chou, K. H., Chen, I. Y., Hsieh, J. C., Wu, C. L., et al. (2010). Structural and cognitive deficits in remitting and non-remitting recurrent depression: a voxel-based morphometric study. Neuroimage 50, 347–356. doi: 10.1016/j.neuroimage.2009.11.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, D., Zhang, H., Liu, Y., Liang, X., Chen, Y., Zheng, Y., et al. (2021). Abnormal Functional Connectivity of Posterior Cingulate Cortex Correlates With Phonemic Verbal Fluency Deficits in Major Depressive Disorder. Front. Neurol. 12:724874. doi: 10.3389/fneur.2021.724874

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Ma, X., Song, L., Fan, J., Wang, W., and Lv, X. (2015). Abnormal spontaneous neural activity in the anterior insular and anterior cingulate cortices in anxious depression. Behav. Brain Res. 281, 339–347. doi: 10.1016/j.bbr.2014.11.047

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Ren, L., Womer, F. Y., Wang, J., Fan, G., Jiang, W., et al. (2014). Alterations in amplitude of low frequency fluctuation in treatment-naive major depressive disorder measured with resting-state fMRI. Hum. Brain Mapp. 35, 4979–4988. doi: 10.1002/hbm.22526

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, S., Ma, R., Luo, Y., Liu, P., Zhao, K., Guo, H., et al. (2021). Facial Expression Recognition and ReHo Analysis in Major Depressive Disorder. Front. Psychol. 12:688376. doi: 10.3389/fpsyg.2021.688376

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Zhao, X., Cheng, Z., Zhang, F., Chang, J., Wang, H., et al. (2017). Regional homogeneity associated with overgeneral autobiographical memory of first-episode treatment-naive patients with major depressive disorder in the orbitofrontal cortex: A resting-state fMRI study. J. Affect. Disord. 209, 163–168. doi: 10.1016/j.jad.2016.11.044

PubMed Abstract | CrossRef Full Text | Google Scholar

Marques, D. R., Gomes, A. A., Caetano, G., and Castelo-Branco, M. (2018). Insomnia Disorder and Brain’s Default-Mode Network. Curr. Neurol. Neurosci. Rep. 18:45. doi: 10.1007/s11910-018-0861-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Nejati, V., Majidinezhad, M., and Nitsche, M. (2022). The role of the dorsolateral and ventromedial prefrontal cortex in emotion regulation in females with major depressive disorder (MDD): A tDCS study. J. Psychiatr. Res. 148, 149–158. doi: 10.1016/j.jpsychires.2022.01.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Ni, L., Liu, R., Yin, Z., Zhao, H., Nedelska, Z., Hort, J., et al. (2016). Aberrant Spontaneous Brain Activity in Patients with Mild Cognitive Impairment and concomitant Lacunar Infarction: A Resting-State Functional MRI Study. J. Alzheimers. Dis. 50, 1243–1254. doi: 10.3233/JAD-150622

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Doherty, J. P. (2004). Reward representations and reward-related learning in the human brain: insights from neuroimaging. Curr. Opin. Neurobiol. 14, 769–776. doi: 10.1016/j.conb.2004.10.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Palejwala, A. H., Dadario, N. B., Young, I. M., O’Connor, K., Briggs, R. G., Conner, A. K., et al. (2021). Anatomy and White Matter Connections of the Lingual Gyrus and Cuneus. World Neurosurg. 151, e426–e437. doi: 10.1016/j.wneu.2021.04.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Raichle, M. E. (2015). The brain’s default mode network. Annu. Rev. Neurosci. 38, 433–447. doi: 10.1146/annurev-neuro-071013-014030

PubMed Abstract | CrossRef Full Text | Google Scholar

Raichle, M. E., MacLeod, A. M., Snyder, A. Z., Powers, W. J., Gusnard, D. A., and Shulman, G. L. (2001). A default mode of brain function. Proc. Natl. Acad. Sci. U S A 98, 676–682. doi: 10.1073/pnas.98.2.676

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramanan, S., and Bellana, B. (2019). A Domain-General Role for the Angular Gyrus in Retrieving Internal Representations of the External World. J. Neurosci. 39, 2978–2980. doi: 10.1523/JNEUROSCI.3231-18.2019

PubMed Abstract | CrossRef Full Text | Google Scholar

Roca, M., Armengol, S., Garcia-Garcia, M., Rodriguez-Bayon, A., Ballesta, I., Serrano, M. J., et al. (2011). Clinical differences between first and recurrent episodes in depressive patients. Compr. Psychiatr. 52, 26–32. doi: 10.1016/j.comppsych.2010.04.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Roca, M., Lopez-Navarro, E., Monzon, S., Vives, M., Garcia-Toro, M., and Garcia-Campayo, J. (2015). Cognitive impairment in remitted and non-remitted depressive patients: A follow-up comparison between first and recurrent episodes. Eur. Neuropsychopharmacol. 25, 1991–1998. doi: 10.1016/j.euroneuro.2015.07.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Rolls, E. T. (2019). The cingulate cortex and limbic systems for action, emotion, and memory. Handb. Clin. Neurol. 166, 23–37. doi: 10.1016/B978-0-444-64196-0.00002-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Seghier, M. L. (2013). The angular gyrus: multiple functions and multiple subdivisions. Neuroscientist 19, 43–61. doi: 10.1177/1073858412440596

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheffield, J. M., and Barch, D. M. (2016). Cognition and resting-state functional connectivity in schizophrenia. Neurosci. Biobehav. Rev. 61, 108–120. doi: 10.1016/j.neubiorev.2015.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, Z., Zhu, J., Ren, L., Qian, M., Shao, Y., Yuan, Y., et al. (2020). Aberrant amplitude low-frequency fluctuation (ALFF) and regional homogeneity (ReHo) in generalized anxiety disorder (GAD) and their roles in predicting treatment remission. Ann. Transl. Med. 8:1319. doi: 10.21037/atm-20-6448

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, H., Luo, L., Yuan, X., Zhang, L., He, Y., Yao, S., et al. (2018). Regional homogeneity and functional connectivity patterns in major depressive disorder, cognitive vulnerability to depression and healthy subjects. J. Affect. Disord. 235, 229–235. doi: 10.1016/j.jad.2018.04.061

PubMed Abstract | CrossRef Full Text | Google Scholar

Talarowska, M., Zajączkowska, M., and Gałecki, P. J. P. D. (2015). Cognitive functions in first-episode depression and recurrent depressive disorder. Psychiatr. Danub. 27, 38–43.

PubMed Abstract | Google Scholar

Thakral, P. P., Madore, K. P., and Schacter, D. L. (2017). A Role for the Left Angular Gyrus in Episodic Simulation and Memory. J. Neurosci. 37, 8142–8149. doi: 10.1523/JNEUROSCI.1319-17.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

Varghese, S., Frey, B. N., Schneider, M. A., Kapczinski, F., and de Azevedo Cardoso, T. (2022). Functional and cognitive impairment in the first episode of depression: A systematic review. Acta. Psychiatr. Scand. 145, 156–185. doi: 10.1111/acps.13385

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Dai, W., Su, Y., Wang, G., Tan, Y., Jin, Z., et al. (2012). Amplitude of low-frequency oscillations in first-episode, treatment-naive patients with major depressive disorder: a resting-state functional MRI study. PLoS One 7:e48658. doi: 10.1371/journal.pone.0048658

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Zhang, A., Yang, C., Li, G., Sun, N., Liu, P., et al. (2020). Enhanced Functional Connectivity Within Executive Function Network in Remitted or Partially Remitted MDD Patients. Front. Psychiatr. 11:538333. doi: 10.3389/fpsyt.2020.538333

PubMed Abstract | CrossRef Full Text | Google Scholar

Wolfers, T., Buitelaar, J. K., Beckmann, C. F., Franke, B., and Marquand, A. F. (2015). From estimating activation locality to predicting disorder: A review of pattern recognition for neuroimaging-based psychiatric diagnostics. Neurosci. Biobehav. Rev. 57, 328–349. doi: 10.1016/j.neubiorev.2015.08.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, X., and Zhang, Y. Q. (2018). A new perspective on the anterior cingulate cortex and affective pain. Neurosci. Biobehav. Rev. 90, 200–211. doi: 10.1016/j.neubiorev.2018.03.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, C. G., and Zang, Y. F. (2010). DPARSF: A MATLAB Toolbox for “Pipeline” Data Analysis of Resting-State fMRI. Front. Syst. Neurosci. 4:13. doi: 10.3389/fnsys.2010.00013

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, C., Zhang, A., Jia, A., Ma, J. X., Sun, N., Wang, Y., et al. (2018). Identify abnormalities in resting-state brain function between first-episode, drug-naive major depressive disorder and remitted individuals: a 3-year retrospective study. Neuroreport 29, 907–916. doi: 10.1097/WNR.0000000000001054

PubMed Abstract | CrossRef Full Text | Google Scholar

Yue, Y., Jiang, Y., Shen, T., Pu, J., Lai, H., and Zhang, B. J. F. i. n (2020). ALFF and ReHo Mapping Reveals Different Functional Patterns in Early- and Late-Onset Parkinson’s Disease. Front. Neurosci. 14:141. doi: 10.3389/fnins.2020.00141

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuksel, 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. Psychiatr. 84, 39–49. doi: 10.1016/j.pnpbp.2018.02.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Zang, Y., Jiang, T., Lu, Y., He, Y., and Tian, L. (2004). Regional homogeneity approach to fMRI data analysis. Neuroimage. 22, 394–400. doi: 10.1016/j.neuroimage.2003.12.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Zang, Y. F., He, Y., Zhu, C. Z., Cao, Q. J., Sui, M. Q., Liang, M., et al. (2007). Altered baseline brain activity in children with ADHD revealed by resting-state functional MRI. Brain Dev. 29, 83–91. doi: 10.1016/j.braindev.2006.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, F. F., Peng, W., Sweeney, J. A., Jia, Z. Y., and Gong, Q. Y. (2018). Brain structure alterations in depression: Psychoradiological evidence. CNS Neurosci. Ther. 24, 994–1003. doi: 10.1111/cns.12835

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Zhu, X., Wang, X., Zhu, X., Zhong, M., Yi, J., et al. (2014). First-episode medication-naive major depressive disorder is associated with altered resting brain function in the affective network. PLoS One 9:e85241. doi: 10.1371/journal.pone.0085241

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Chen, Y., Wei, W., Yang, X., Meng, Y., Yu, H., et al. (2021). Changes in Regional Homogeneity of Medication-Free Major Depressive Disorder Patients With Different Onset Ages. Front. Psychiatr. 12:713614. doi: 10.3389/fpsyt.2021.713614

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, C., Zhu, J., Liu, X., Pu, C., Lai, Y., Chen, L., et al. (2018). Structural and functional brain abnormalities in schizophrenia: A cross-sectional study at different stages of the disease. Prog. Neuropsychopharmacol. Biol. Psychiatr. 83, 27–32. doi: 10.1016/j.pnpbp.2017.12.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, H., Li, F., Bo, Q., Li, X., Yao, L., Yao, Z., et al. (2018). The dynamic characteristics of the anterior cingulate cortex in resting-state fMRI of patients with depression. J. Affect. Disord. 227, 391–397. doi: 10.1016/j.jad.2017.11.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Zu, S., Wang, D., Fang, J., Xiao, L., Zhu, X., Wu, W., et al. (2021). Comparison of Residual Depressive Symptoms. Functioning, and Quality of Life Between Patients with Recurrent Depression and First Episode Depression After Acute Treatment in China. Neuropsychiatr Dis Treat 17, 3039–3051. doi: 10.2147/NDT.S317770

PubMed Abstract | CrossRef Full Text | Google Scholar