A Comparative Study of Regional Homogeneity of Resting-State fMRI Between the Early-Onset and Late-Onset Recurrent Depression in Adults

Background: Neurobiological mechanisms underlying the recurrence of major depressive disorder (MDD) at different ages are unclear, and this study used the regional homogeneity (ReHo) index to compare whether there are differences between early onset recurrent depression (EORD) and late onset recurrent depression (LORD).

Methods: Eighteen EORD patients, 18 LORD patients, 18 young healthy controls (HCs), and 18 older HCs were included in the rs-fMRI scans. ReHo observational metrics were used for image analysis and further correlation of differential brain regions with clinical symptoms was analyzed.

Results: ANOVA analysis revealed significant differences between the four groups in ReHo values in the prefrontal, parietal, temporal lobes, and insula. Compared with EORD, the LORD had higher ReHo in the right fusiform gyrus/right middle temporal gyrus, left middle temporal gyrus/left angular gyrus, and right middle temporal gyrus/right angular gyrus, and lower ReHo in the right inferior frontal gyrus/right insula and left superior temporal gyrus/left insula. Compared with young HCs, the EORD had higher ReHo in the right inferior frontal gyrus/right insula, left superior temporal gyrus/left insula, and left rolandic operculum gyrus/left superior temporal gyrus, and lower ReHo in the left inferior parietal lobule, right inferior parietal lobule, and left middle temporal gyrus/left angular gyrus. Compared with old HCs, the LORD had higher ReHo in the right fusiform gyrus/right middle temporal gyrus, right middle temporal gyrus/right angular gyrus, and left rolandic operculum gyrus/left superior temporal gyrus, and lower ReHo in the right inferior frontal gyrus/right insula. ReHo in the right inferior frontal gyrus/right insula of patients with LORD was negatively correlated with the severity of 17-item Hamilton Rating Scale for Depression (HAMD-17) scores (r = −0.5778, p = 0.0120).

Conclusion: Adult EORD and LORD patients of different ages have abnormal neuronal functional activity in some brain regions, with differences closely related to the default mode network (DMN) and the salience network (SN), and patients of each age group exhibit ReHo abnormalities relative to matched HCs.

Clinical Trial Registration: [http://www.chictr.org.cn/], [ChiCTR1800014277].

Introduction

Major depressive disorder (MDD) is a common clinical psychiatric disorder that is characterized by depressed mood, slowed thinking, reduced interest, and diminished cognitive function as the main clinical manifestations, and is one of the main causes of disability (Rotenstein et al., 2016; Cladder-Micus et al., 2018). It is estimated that MDD will become the second leading cause of disability in the world (Ferrari et al., 2013). Approximately 25% of MDD patients relapse within 6 months of discharge (Harlev et al., 2021), approximately one-third of MDD patients are prone to relapse within a year (Cosci et al., 2020), and the severity of MDD increases with the number of depressive relapses (Harlev et al., 2021). There is a growing awareness that the challenge of depression is relapse prevention rather than recovery (Fava et al., 2017). Therefore, understanding the pathogenesis of recurrent depression is of great importance for clinical work.

Previous studies have found that patients with recurrent depression at different ages have different clinical features (Driscoll et al., 2005; Hammen et al., 2008; Sung et al., 2013). Early onset recurrent depression (EORD) is associated with a higher risk of life-long persistent depression, a higher risk of suicide, worse social cognitive function, and greater anxiety than late onset recurrent depression (LORD) (Hammen et al., 2008; Sung et al., 2013). And LORD is also associated with more cognitive dysfunction than EORD (Driscoll et al., 2005). A study has shown that EORD has a greater risk of relapse than LORD (Cassano et al., 1993). Therefore, recurrent depression in different age groups may be associated with different neuropathological mechanisms.

Moreover, recurrent depression is more heritable and increases the risk of depression in offspring (Burcusa and Iacono, 2007; Benjet et al., 2020; Jaffee et al., 2021). For example, one study found that adolescent children whose parents had recurrent depression had a 4.21-fold greater probability of reaching MDD than adolescents whose parents were never depressed (Jaffee et al., 2021). A meta analysis showed a degree of specificity in different subtypes of depression, and understanding their specificity can be a guide to clinical treatment (Harald and Gordon, 2012).

However, the clinical characteristics of EORD and LORD depend on the different cut-off ages (Cassano et al., 1993; Driscoll et al., 2005; Hammen et al., 2008). An age-specific study of recurrent depression in adolescents showed that EORD recurred between 15 and 20 years of age, and LORD after 20 years of age (Hammen et al., 2008). A study of recurrent depression in older adults showed that the age of EORD was defined as before 59 years of age, while LORD was defined as after 60 years of age (Driscoll et al., 2005). Another study defined the age of EORD before 44 years of age and LORD after 45 years of age (Cassano et al., 1993). The above studies suggest that the biological evidence for the age limits of EORD and LORD is unclear and there is no unified consensus on age classification criteria. In this study, the age of EORD was defined as 18 to 29 years and LORD was defined from 30 to 45 years based on previous studies to avoid the effect of (>45 years) cerebrovascular disease (Shen et al., 2017).

In recent years, with the development of neuroimaging technology, resting-state functional MRI (rs-MRI) has been widely used in insomnia (Li et al., 2019), schizophrenia (Huang et al., 2020), autism (Shi et al., 2020), and other psychiatric disorders. Regional homogeneity (ReHo) is an important indicator for the study of rs-fMRI. ReHo can reflect the synchronization of brain functional activity state in localized regions of whole-brain voxels and is used to assess the level of coordination of neural activity in local brain regions (Zang et al., 2004). In addition, the ReHo indicator has been used to study changes in neuronal functional activity in depressive subtypes of the brain, and abnormalities have been found in brain regions such as the frontal lobe, precuneus, and insula (Chen et al., 2012; Shen et al., 2017; Zhang et al., 2021). The ReHo values of these abnormal brain regions are closely correlated with the default mode network (DMN) and the salience network (SN). The correlation of ReHo values with abnormal brain regions in subtypes of depression is useful for understanding the degree of severity of the disease and neuroimaging markers (Shen et al., 2017; Liu et al., 2021). Previous studies have also found that differences in brain function across subtypes of depression are associated with abnormalities in DMN and SN (Li et al., 2017; Cai et al., 2021). However, there is a lack of neuroimaging studies on the differences between EORD and LORD, and the differences in regional neural activity between EROD and LORD are unclear.

In this study, we focused on the differences in local brain functional activity between EORD (18–29 years) and LORD (30–45 years) and whether there was a correlation between the depression group and clinical symptoms. Based on previous studies of age differences in subtypes of depression, we hypothesized that (1) differences in local brain functional activity in EORD and LORD may be closely related to the DMN and SN. (2) Increased or decreased local brain function of ReHo in EORD and LORD may be associated with clinical depressive symptoms.

Materials and Methods

Subjects

A total of 36 patients (18 with EORD and 18 with LORD) were diagnosed with recurrent depression from the Department of Psychosomatic Medicine, Guang’anmen Hospital, China Academy of Chinese Medical Sciences, the Department of Psychiatry, Beijing First Hospital of Integrative Medicine, the Department of Psychiatry, Yuquan Hospital, Qinghua University, and the Department of Psychiatry, Xuanwu Hospital, Capital Medical University. The inclusion criteria were as follows: (1) According to the diagnostic criteria for depression in the Diagnostic and Statistical Manual of Mental Disorders (DSM-V) of the American Psychiatric Association, all patients were diagnosed by two experienced psychiatrists; (2) All patients recurred after antidepressant treatment, relapsed between the ages of 18 and 45 years, and were off medication for at least 4 weeks prior to admission; (3) 17-item Hamilton Rating Scale for Depression (HAMD-17) score > 17; and (4) right-handedness. 36 gender- and age-matched healthy controls (HCs; 25 females and 11 males) were also included in the HCs: (1) HAMD-17 score < 7; (2) right-handedness; and (3) no history of mental illness.

Exclusion criteria of patients and HCs include the following items: (1) suffering from serious mental illness and other diseases such as cardiovascular and cerebrovascular; (2) with a history of drug and alcohol abuse; (3) had any MRI contraindications, such as heart pacemaker, metal fixed false teeth, or severe claustrophobia; and (4) were pregnant or lactating.

Clinical Materials and Subgroups

Data were collected for the study population including gender, age, education level, frequency of recurrence, and duration of disease. Subjects were assessed for depressive symptoms by an experienced psychiatrist using the HAMD-17 scale. Based on previous studies on age division of first-episode depression to avoid the influence of cerebrovascular diseases (Shen et al., 2016, 2017), we initially divided all recurrent depression patients with recurrent depression into EORD group (18–29 years) and LORD group (30–45 years). The HCs were also matched corresponding to two subgroups: young HCs (18–29 years) and old HCs (30–45 years).

Scan Acquisition

All subjects in this study underwent data acquisition using a Magneton Skyra 3.0 T MRI scanner (Siemens, Germany). Subjects were informed prior to the scan to remain awake during the scan and to avoid active thinking activities. During the scan, the subject is required to use earplugs and wear noise-canceling headphones, use a hood to immobilize the head, and lie flat on the examination bed. The scanning procedure contains a localizer, a high resolution three-dimensional T1-weighted imaging (3D-T1WI), and a blood oxygenation level-dependent fMRI (BOLD-fMRI).

The scanning parameters were as follows: 3D-T1WI: time repetition (TR)/time echo (TE) = 2500/2.98 ms, flip angle = 7°, matrix = 64 × 64, field of view (FOV) = 256 mm × 256 mm, slice thickness = 1 mm, slice number = 48, slices =192, scanning time 6 min 3 s; BOLD-fMRI: TR/TE = 2000/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, and 200 volumes were obtained, scanning time 6 min 40 s.

fMRI Data Analysis

fMRI Data Preprocessing

The rs-fMRI data were pre-processed using Data Processing Assistant for Resting-State fMRI (DPARSF) software (DPABI5.0; Chao-Gan and Yu-Feng, 2010) in MATLAB (Mathworks, Inc., Natick, MA, United States).¹ First, image format conversion was performed to convert the raw data DICOM format to NIFTI format. Next, the first 10 time points were excluded. Then, slice timing, realign, was performed to remove subjects with head movement translations exceeding 2.0 mm and rotations exceeding 2.0°. Spatial normalization was performed by normalizing the functional images to Montreal Neurological Institute (MNI) space. Linear regression was performed on the covariates of head movement, white matter (WM), and cerebrospinal fluid (CSF) signals to reduce the effects. Finally, we removed the linear drift and set the filter to 0.01–0.08 Hz to reduce the effect of noise.

ReHo Analysis

The pre-processed data images were analyzed using DPARSF software, and ReHo was calculated by Kendall correlation coefficient (Kendall & Gibbons, 1990), based on voxels to calculate the synchronization of the time series of a given voxel with the time series changes of its 26 adjacent voxels, and a ReHo map was obtained for each subject. The ReHo map of each subject was divided by the whole-brain mean ReHo value to obtain a normalized ReHo map. Finally, the smReHo maps were obtained by smoothing using a 6 mm × 6 mm × 6 mm Gaussian smoothing kernel for subsequent statistical analysis.

Statistical Analysis

Clinical Data Analysis

The clinical data were analyzed with SPSS 23.0 statistics software (IBM Corp, Somers, New York, United States). One-way analysis of variance (ANOVA) was used to compare age and education level between the four groups, and the chi-square test was used to compare gender. The duration of illness, frequency of recurrence, and HAMD-17 scores were compared between patient groups using a two-sample t-test. p < 0.05 was statistically significant.

fMRI Data Analysis

Image data statistics were analyzed using the DPARSF toolbox, and a voxel-based one-way ANOVA was performed to compare whole-brain ReHo maps of the four groups. Gender, age, education level, and framewise displacement (FD) metric (derived from Jenkinson’s formula) of the four groups of subjects were used as covariates, brain areas with ReHo differences between the four groups were corrected for Gaussian random fields (GRF), and corrected cluster levels were settled as p < 0.05 and threshold voxel levels p < 0.005 were defined as statistically different. The time series means of the peak voxel ReHo values of the four contrasting brain regions were extracted, and post-hoc two-sample t-tests were performed on the ReHo values between the paired groups (EORD vs. young HCs, LORD vs. old HCs, LORD vs. EORD) using SPSS 23.0 software, and the results were Bonferroni corrected (p < 0.01).

In order to verify the relationship between ReHo values and clinical symptoms, Pearson correlation analysis was performed by extracting ReHo values of abnormal brain regions and HAMD-17 scores in patient groups, separately, and p < 0.05 was statistically significant.

Results

Characteristics of Research Samples

A total of 18 patients with EROD, 18 patients with LORD, 18 young HCs, and 18 old HCs met the study criteria. There were no statistical differences between the MDD group and the HCs and in terms of gender and years of education, and there were statistical differences in age comparisons between the four groups. There were no statistical differences in disease duration, frequency of recurrence, and HAMD-17 scores between the two MDD groups (Table 1).

TABLE 1

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

ReHo: Group Differences

A one-way ANOVA was performed on the ReHo values of the four groups using age, gender, years of education, and mean frame displacement as covariates. ReHo was found to be significantly significant difference in the right inferior frontal gyrus/right insula, left superior temporal gyrus/left insula, left inferior parietal lobule, right inferior parietal lobule, right fusiform gyrus/right middle temporal gyrus, left middle temporal gyrus/left angular gyrus, right middle temporal gyrus/right angular gyrus, and left rolandic operculum gyrus/left superior temporal gyrus significant differences between superior gyrus groups (Table 2; Figure 1).

TABLE 2

Table 2. One-way ANOVA between the four groups for differences in ReHo values Brain areas.

FIGURE 1

Figure 1. Statistical maps showing ANOVA results for ReHo abnormalities between EORD or LORD patients and HCs (GRF corrected). Color bar indicates red to yellow showing enhancement of ReHo values.

Compared with young HCs, the EORD had higher ReHo in the right inferior insular frontal gyrus/right insula, left superior temporal gyrus/left insula, and left rolandic operculum gyrus/superior temporal gyrus, and lower ReHo values in the left inferior parietal lobule, right inferior parietal lobule, and left middle temporal gyrus/left angular gyrus (Figure 2).

FIGURE 2

Figure 2. Post-hoc two-sample t-tests (Bonferroni corrected) comparison showing ReHo values differences at peak voxel between each pair group (EORD vs. young HCs, LORD vs. old HCs, and LORD vs. EORD). EORD, early onset recurrent depression; LORD, late onset recurrent depression; young HCs, young healthy controls; and old HCs, old healthy controls. Frontal_Inf_Oper_R/Insula_R = right inferior frontal gyrus/right insula; Temporal_Sup_L/Insula_L = left superior temporal gyrus/left insula; Parietal_Inf_L = left inferior parietal lobule; Parietal_Inf_R = right inferior parietal lobule; Fusiform_R/Temporal_Mid_R = right fusiform gyrus/right middle temporal gyrus; Temporal_Mid_L/Angular_L = left middle temporal gyrus/left angular gyrus; Temporal_Mid_R/Angular_R = right middle temporal gyrus/right angular gyrus; and Rolandic_Oper_L/Temporal_Sup_L = left rolandic operculum gyrus/left superior temporal gyrus.

Abnormal ReHo in LORD Patients vs Old HC

Compared with old HCs, the LORD had higher ReHo in the right fusiform gyrus/right middle temporal gyrus, right middle temporal gyrus/right angular gyrus, left rolandic operculum gyrus/left superior temporal gyrus, and lower ReHo values in the right inferior frontal gyrus/right insula (Figure 2).

Abnormal ReHo in LORD Patients vs EORD

Compared with EORD, the LORD group had higher ReHo in the right fusiform gyrus/right middle temporal gyrus, left middle temporal gyrus/left angular gyrus, and right middle temporal gyrus/right angular gyrus, and lower ReHo in the right inferior frontal gyrus/right insula and left superior temporal gyrus/left insula (Figure 2).

Correlations Between the Abnormal ReHo and Depressive Symptoms in the Patients

To test whether there was a relationship between the severity of depressive symptoms and abnormal ReHo in the depressed group in this study, we performed further correlation analyses. We found a negative correlation between HAMD-17 scores and ReHo values in the right inferior frontal gyrus/right insula in the LORD group (r = −0.5778, p = 0.0120; Figure 3).

FIGURE 3

Figure 3. The correlation of the ReHo and HDRS symptom in LORD patients. LORD, late onset recurrent depression.

Discussion

To our knowledge, this is the first to analyze and compare the regional neural activity of LORD and EORD based on rs-MRI technique using the ReHo method. The results of this study found that adult EORD and LORD have abnormal neuronal activity in some of the same brain regions at different ages. And relative to matched controls, different ReHo abnormalities were demonstrated at each age. Correlation analysis showed that ReHo values in the right inferior frontal gyrus/right insula were negatively correlated with clinical scales and could be used as a neuroimaging marker to distinguish LORD from EORD. This study provides new insights into functional brain activity in recurrent depression at different ages.

The results of this study showed that LORD had higher ReHo in the right fusiform gyrus/right middle temporal gyrus, left middle temporal gyrus/left angular gyrus, and right middle temporal gyrus/right angular gyrus compared with EORD. The concept of DMN was introduced by (Raichle et al., 2001) in 2001, which suggests that some brain areas are more active than others in a quiet, resting state without a task and are “activated,” while functional activity decreases when performing a task and becomes “negatively activated.” This region consists of the ventral medial prefrontal cortex, posterior cingulate/retosplenial cortex, lateral temporal cortex, dorsal medial prefrontal cortex, and inferior parietal lobules (including the angular gyrus) and is involved in the extraction of situational memory, detection of surroundings, and introspective states (Raichle et al., 2001; Raichle and Snyder, 2007; Buckner et al., 2008). The results of this study showed that most of the brain regions belonged to DMN. Most studies have regarded abnormal DMN function as a neurobiological marker of MDD (Gao et al., 2018; Goldstein-Piekarski et al., 2018; Sendi et al., 2021), DMN also changes in different stages of disease and at different ages, which is closely related to changes in clinical depression symptoms, and different therapeutic effects will also produce different responses to DMN (Huang et al., 2020; Pan et al., 2020). The middle temporal gyrus and angular gyrus are important constituent parts of the DMN (Buckner et al., 2008). The middle temporal gyrus not only participates in self-referential processing, but also plays an important role in autobiographical memory and regulates human emotional and mental activities (Sheldon et al., 2016). The angular gyrus is located at the junction of the parietal, temporal, and occipital lobes and is primarily involved in ontological functions such as semantic processing, number processing, attention, and memory (Seghier, 2013). Previous studies have shown abnormal DMN functional activity in patients with recurrent depression (Marchetti et al., 2012; Yan et al., 2019; Berwian et al., 2020). It was found that the recurrent depression group had hyperactivation in the right middle temporal gyrus compared to the HCs (Liu et al., 2017, 2019). Another study also found that the recurrent depression had abnormal bold signaling in the angular gyrus compared to the HCs (Yüksel et al., 2018). All of the above studies suggest abnormalities in the middle temporal gyrus and angular gyrus in patients with recurrent depression. The fusiform gyrus belongs to the visual cognitive network and together with the inferior temporal gyrus forms the ventral visual pathway, which is involved in the recognition of objects, words, and faces. The neural response of the fusiform gyrus increases with the intensity of the sad stimulus (Wu et al., 2011). ReHo was higher in the left superior temporal gyrus (Chen et al., 2012) and right fusiform gyrus (Shen et al., 2017) in late onset depression (LOD) compared with early onset depression (EOD), suggesting that compensatory elevation of the left superior temporal gyrus and right fusiform gyrus may have a pathogenetic mechanism for LOD. This is broadly consistent with the results of the present study. Therefore, the results of the present study suggest that ReHo is elevated in the right fusiform/right middle temporal gyrus, left middle temporal/left angular gyrus, and right middle temporal/right angular gyrus in LORD compared with EORD, and this partial compensatory elevation of DMN may be one of the differences in the brain mechanisms of recurrent depression in these two subtypes.

We also found that LORD had lower ReHo than EORD in the right inferior frontal gyrus/right insula and left superior temporal gyrus/left insula. The right inferior frontal gyrus, part of the ventral lateral prefrontal cortex, plays an important role in response inhibition and is primarily responsible for downregulating negative emotional responses and inhibiting unwanted information or inappropriate behavior (Kravitz et al., 2011). The insula is a cortical structure located deep in the brain and is involved in the processing of information related to emotion, attention, visceral sensation, etc., that is transmitted to the brain, as well as being involved in taste and vision (Sprengelmeyer et al., 2011). The insula cortex is connected to the frontal limbic area (Guo et al., 2015). Previous studies have found that metabolism in the insula cortex is enhanced in HCs when they recall sad events (Reiman et al., 1997). These two brain areas are closely related to the SN, which has been shown to be a functional brain network involved in perception, cognition, emotion, and social awareness (Menon and Uddin, 2010; Ham et al., 2013; Uddin, 2015; Chand and Dhamala, 2016; Touroutoglou et al., 2016). The SN has an important regulatory role for the DMN and the executive control network (ECN) (Touroutoglou et al., 2016). Previous, it was found that patients with recurrent depression had lower amplitude of low frequency fluctuations (ALFF) (Liu et al., 2017) and gray matter volume (Stratmann et al., 2014) in the right insula compared to HCs. It has also been found that middle-aged and older women with a history of depression have reduced FC between SN (right insula seed) and ECN (Vega et al., 2020). All of the above studies found that abnormalities in the function of the insula were closely associated with recurrent depression. Previous studies have found significantly enhanced ReHo in the right insula of LOD compared to old HCs, suggesting SN abnormalities in patients with LOD depression (Shen et al., 2017). Therefore, the results of this study suggest that there are differences in SN between LORD and EORD, which may also be the neuropathological pathogenesis of these two subtypes of recurrent depression. Further correlation analysis in this study showed that ReHo values in the right inferior insular frontal gyrus/right insula of the LORD group were negatively correlated with HAMD-17 scale scores. It is suggested that this region can be used as a neuroimaging marker with high sensitivity and specificity to distinguish between the two subtypes of recurrent depression. More research is needed in the future to clarify this point.

On the other hand, we also found that the ReHo in the left superior temporal gyrus/left insula were lower in the LORD compared to the EORD. The superior temporal gyrus is part of the auditory-verbal center and plays an important role in the processing of emotions, memory, and mental activity (Yao et al., 2018). Previous studies found that the LOD had higher ReHo values in the left superior temporal gyrus compared to the EOD, suggesting that compensatory elevation of ReHo in the left superior temporal gyrus may be the pathogenesis of LOD patients (Chen et al., 2012). Therefore, the results of this study suggest that the left superior temporal gyrus/left insula is also an important brain region in which differences exist between the two subgroups of patients with recurrent depression.

Interestingly, we observed that both the LORD and EROD groups together had increased ReHo in the left rolandic operculum gyrus/left superior temporal gyrus compared to HCs. The rolandic operculum gyrus and superior temporal gyrus belong to the auditory network (Oh et al., 2018; Koshiyama et al., 2020), which has the function of integrating auditory stimuli and processing them, and is closely related to human emotions (Phillips et al., 2001; Wang et al., 2019). Previous studies have found that depression scale scores are positively correlated with gray matter in the left rolandic operculum gyrus and left superior temporal gyrus, suggesting that these two brain regions also play a role in mood regulation (Besteher et al., 2017). One study showed that the insula is closely related to the operculum gyrus and that the insula/operculum gyrus is involved in interoception and interoceptive awareness, processes signals critical for self-awareness (Blefari et al., 2017). A meta-analysis showed that depressed patients had significantly increased activation of the right superior temporal gyrus during working memory (WM) compared to the HCs (Wang et al., 2015). In addition, abnormal fMRI connections between the right superior temporal lobe and subgenual cingulate cortex were also found in previous studies in patients with recurrent depression (Lythe et al., 2015). The above findings favorably support our observation of ReHo anomalies in the left rolandic operculum gyrus and left superior temporal gyrus. Therefore, these findings suggest that (Rotenstein et al., 2016) patients with recurrent depression may have features of auditory network dysfunction. (Rotenstein et al., 2016) ReHo abnormalities in the left rolandic operculum gyrus and left superior temporal gyrus may be important brain regions in distinguishing the recurrent depression group from the HCs, independent of age of onset.

However, there are still several limitations to this study. First, the small sample size of this study will affect the statistical validity to some extent, and a larger sample size is needed to confirm or overturn the current results. Second, only patients with recurrent depression were recruited in this study. Although there was a discontinuation period of at least 4 weeks, there may still be an influence of potential factors such as antidepressants. Finally, although there was no statistical difference in the frequency of recurrence between the two subtypes of depression groups, the study of first recurrence of depression in different age groups is more clinically relevant and further research in this area will be enhanced in the future.

Conclusion

In conclusion, this study was based on rs-fMRI technique and used ReHo index analysis to initially explore the differences in regional neural activity between LORD and EORD. We found that although LORD and EORD had similar clinical symptoms, abnormal changes in neurological functional activity existed in some of the same brain regions, suggesting that different pathogenesis may exist in patients with different age onset. In addition, patients of each age group also exhibited different ReHo abnormalities relative to HCs.

Data Availability Statement

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

Ethics Statement

This study was reviewed and approved by the Ethics Committee of Guang’anmen Hospital, Chinese Academy of Traditional Chinese Medicine, China (NO. 2017-021-SQ). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

J-fS and L-mC drafted the manuscript and participated in data collection and analysis. J-kH, D-qG, C-lG, YM, YL, YH, and F-qX involved in data analysis and project design work. J-lF involved in the design of the experimental study and the revision of the manuscript. 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), and National Natural Science Foundation of China (82174282, 81774433), 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

The authors thank all the patients and healthy people who participated in the experiment for their support.

Footnotes

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

References

Benjet, C., Albor, Y. C., Bocanegra, E. S., Borges, G., Méndez, E., Casanova, L., et al. (2020). Incidence and recurrence of depression from adolescence to early adulthood: a longitudinal follow-up of the Mexican adolescent mental health survey. J. Affect. Disord. 263, 540–546. doi: 10.1016/j.jad.2019.11.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Berwian, I. M., Wenzel, J. G., Kuehn, L., et al. (2020). The relationship between resting-state functional connectivity, antidepressant discontinuation and depression relapse. Sci. Rep. 10:22346. doi: 10.1038/s41598-020-79170-9

CrossRef Full Text | Google Scholar

Besteher, B., Gaser, C., Langbein, K., Dietzek, M., Sauer, H., and Nenadić, I. (2017). Effects of subclinical depression, anxiety and somatization on brain structure in healthy subjects. J. Affect. Disord. 215, 111–117. doi: 10.1016/j.jad.2017.03.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Blefari, M. L., Martuzzi, R., Salomon, R., Bello-Ruiz, J., Herbelin, B., Serino, A., et al. (2017). Bilateral Rolandic operculum processing underlying heartbeat awareness reflects changes in bodily self-consciousness. Eur. J. Neurosci. 45, 1300–1312. doi: 10.1111/ejn.13567

PubMed Abstract | CrossRef Full Text | Google Scholar

Buckner, R. L., Andrews-Hanna, J. R., and Schacter, D. L. (2008). The brain’s default network: anatomy, function, and relevance to disease. Ann. N. Y. Acad. Sci. 1124, 1–38. doi: 10.1196/annals.1440.011

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

Cai, Y., Elsayed, N. M., and Barch, D. M. (2021). Contributions from resting state functional connectivity and familial risk to early adolescent-onset MDD: results from the adolescent brain cognitive development study. J. Affect. Disord. 287, 229–239. doi: 10.1016/j.jad.2021.03.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Cassano, G. B., Akiskal, H. S., Savino, M., Soriani, A., Musetti, L., and Perugi, G. (1993). Single episode of major depressive disorder. First episode of recurrent mood disorder or distinct subtype of late-onset depression? Eur. Arch. Psychiatry Clin. Neurosci. 242, 373–380. doi: 10.1007/BF02190251

CrossRef Full Text | Google Scholar

Chand, G. B., and Dhamala, M. (2016). The salience network dynamics in perceptual decision-making. NeuroImage 134, 85–93. doi: 10.1016/j.neuroimage.2016.04.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Chao-Gan, Y., and Yu-Feng, Z. (2010). DPARSF: A MATLAB Toolbox for “Pipeline” Data Analysis of Resting-State fMRI. Front. Syst. Neurosci. 4:13. doi: 10.3389/fnsys.2010.00013

CrossRef Full Text | Google Scholar

Chen, J. D., Liu, F., Xun, G. L., Chen, H. F., Hu, M. R., Guo, X. F., et al. (2012). Early and late onset, first-episode, treatment-naive depression: same clinical symptoms, different regional neural activities. J. Affect. Disord. 143, 56–63. doi: 10.1016/j.jad.2012.05.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Cladder-Micus, M. B., Speckens, A. E. M., Vrijsen, J. N., Donders, A. R. T., Becker, E. S., and Spijker, J. (2018). Mindfulness-based cognitive therapy for patients with chronic, treatment-resistant depression: a pragmatic randomized controlled trial. Depress. Anxiety 35, 914–924. doi: 10.1002/da.22788

PubMed Abstract | CrossRef Full Text | Google Scholar

Cosci, F., Guidi, J., Mansueto, G., and Fava, G. A. (2020). Psychotherapy in recurrent depression: efficacy, pitfalls, and recommendations. Expert. Rev. Neurother. 20, 1169–1175. doi: 10.1080/14737175.2020.1804870

PubMed Abstract | CrossRef Full Text | Google Scholar

Driscoll, H. C., Basinski, J., Mulsant, B. H., Butters, M. A., Dew, M. A., Houck, P. R., et al. (2005). Late-onset major depression: clinical and treatment-response variability. Int. J. Geriatr. Psychiatry 20, 661–667. doi: 10.1002/gps.1334

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

Ferrari, A. J., Charlson, F. J., Norman, R. E., Patten, S. B., Freedman, G., Murray, C. J. L., et al. (2013). Burden of depressive disorders by country, sex, age, and year: findings from the global burden of disease study 2010. PLoS Med. 10:e1001547. doi: 10.1371/journal.pmed.1001547

CrossRef Full Text | Google Scholar

Gao, W., Chen, S., Biswal, B., Lei, X., and Yuan, J. (2018). Temporal dynamics of spontaneous default-mode network activity mediate the association between reappraisal and depression. Soc. Cogn. Affect. Neurosci. 13, 1235–1247. doi: 10.1093/scan/nsy092

PubMed Abstract | CrossRef Full Text | Google Scholar

Goldstein-Piekarski, A. N., Staveland, B. R., Ball, T. M., Yesavage, J., Korgaonkar, M. S., and Williams, L. M. (2018). Intrinsic functional connectivity predicts remission on antidepressants: a randomized controlled trial to identify clinically applicable imaging biomarkers. Transl. Psychiatry 8:57. doi: 10.1038/s41398-018-0100-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, W., Liu, F., Xiao, C., Zhang, Z., Liu, J., Yu, M., et al. (2015). Decreased insular connectivity in drug-naive major depressive disorder at rest. J. Affect. Disord. 179, 31–37. doi: 10.1016/j.jad.2015.03.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Ham, T., Leff, A., de Boissezon, X., Joffe, A., and Sharp, D. J. (2013). Cognitive control and the salience network: an investigation of error processing and effective connectivity. J. Neurosci. 33, 7091–7098. doi: 10.1523/JNEUROSCI.4692-12.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Hammen, C., Brennan, P. A., Keenan-Miller, D., and Herr, N. R. (2008). Early onset recurrent subtype of adolescent depression: clinical and psychosocial correlates. J. Child Psychol. Psychiatry 49, 433–440. doi: 10.1111/j.1469-7610.2007.01850.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Harald, B., and Gordon, P. (2012). Meta-review of depressive subtyping models. J. Affect. Disord. 139, 126–140. doi: 10.1016/j.jad.2011.07.015

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. Psychiatry 12:572660. doi: 10.3389/fpsyt.2021.572660

CrossRef Full Text | Google Scholar

Huang, H., Botao, Z., Jiang, Y., Tang, Y., Zhang, T., Tang, X., et al. (2020). Aberrant resting-state functional connectivity of salience network in first-episode schizophrenia. Brain Imaging Behav. 14, 1350–1360. doi: 10.1007/s11682-019-00040-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Jaffee, S. R., Sligo, J. L., McAnally, H. M., Bolton, A. E., Baxter, J. M., and Hancox, R. J. (2021). Early-onset and recurrent depression in parents increases risk of intergenerational transmission to adolescent offspring. J. Child Psychol. Psychiatry 62, 979–988. doi: 10.1111/jcpp.13356

PubMed Abstract | CrossRef Full Text | Google Scholar

Kendall, M. G., and Gibbons, J. D. (1990). in Rank Correlation Methods. 5th Edn. ed. E. Arnold (London New York, NY: Oxford University Press).

Google Scholar

Koshiyama, D., Miyakoshi, M., Tanaka-Koshiyama, K., Joshi, Y. B., Molina, J. L., Sprock, J., et al. (2020). Neurophysiologic Characterization of Resting State Connectivity Abnormalities in Schizophrenia Patients. Front. Psychiatry 11:608154. doi: 10.3389/fpsyt.2020.608154

PubMed Abstract | CrossRef Full Text | Google Scholar

Kravitz, D. J., Saleem, K. S., Baker, C. I., and Mishkin, M. (2011). A new neural framework for visuospatial processing. Nat. Rev. Neurosci. 12, 217–230. doi: 10.1038/nrn3008

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, W., Wang, Y., Ward, B. D., Antuono, P. G., Li, S. J., and Goveas, J. S. (2017). Intrinsic inter-network brain dysfunction correlates with symptom dimensions in late-life depression. J. Psychiatr. Res. 87, 71–80. doi: 10.1016/j.jpsychires.2016.12.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M., Wang, R., Zhao, M., Zhai, J., Liu, B., Yu, D., et al. (2019). Abnormalities of thalamus volume and resting state functional connectivity in primary insomnia patients. Brain Imaging Behav. 13, 1193–1201. doi: 10.1007/s11682-018-9932-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C. H., Ma, X., Yuan, Z., Song, L. P., Jing, B., Lu, H. Y., et al. (2017). Decreased resting-state activity in the Precuneus is associated With depressive episodes in recurrent depression. J. Clin. Psychiatry 78, e372–e382. doi: 10.4088/JCP.15m10022

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C. H., Tang, L. R., Gao, Y., Zhang, G. Z., Li, B., Li, M., et al. (2019). Resting-state mapping of neural signatures of vulnerability to depression relapse. J. Affect. Disord. 250, 371–379. doi: 10.1016/j.jad.2019.03.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, P., Tu, H., Zhang, A., Yang, C., Liu, Z., Lei, L., et al. (2021). Brain functional alterations in MDD patients with somatic symptoms: a resting-state fMRI study. J. Affect. Disord. 295, 788–796. doi: 10.1016/j.jad.2021.08.143

PubMed Abstract | CrossRef Full Text | Google Scholar

Lythe, K. E., Moll, J., Gethin, J. A., Workman, C. I., Green, S., Lambon Ralph, M. A., et al. (2015). Self-blame-selective hyperconnectivity between anterior temporal and subgenual cortices and prediction of recurrent depressive episodes. JAMA Psychiat. 72, 1119–1126. doi: 10.1001/jamapsychiatry.2015.1813

PubMed Abstract | CrossRef Full Text | Google Scholar

Marchetti, I., Koster, E. H., Sonuga-Barke, E. J., and De Raedt, R. (2012). The default mode network and recurrent depression: a neurobiological model of cognitive risk factors. Neuropsychol. Rev. 22, 229–251. doi: 10.1007/s11065-012-9199-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Menon, V., and Uddin, L. Q. (2010). Saliency, switching, attention and control: a network model of insula function. Brain Struct. Funct. 214, 655–667. doi: 10.1007/s00429-010-0262-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Oh, S. Y., Boegle, R., Ertl, M., Stephan, T., and Dieterich, M. (2018). Multisensory vestibular, vestibular-auditory, and auditory network effects revealed by parametric sound pressure stimulation. NeuroImage 176, 354–363. doi: 10.1016/j.neuroimage.2018.04.057

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, F., Xu, Y., Zhou, W., Chen, J., Wei, N., Lu, S., et al. (2020). Disrupted intrinsic functional connectivity of the cognitive control network underlies disease severity and executive dysfunction in first-episode, treatment-naive adolescent depression. J. Affect. Disord. 264, 455–463. doi: 10.1016/j.jad.2019.11.076

PubMed Abstract | CrossRef Full Text | Google Scholar

Phillips, M. D., Lowe, M. J., Lurito, J. T., Dzemidzic, M., and Mathews, V. P. (2001). Temporal lobe activation demonstrates sex-based differences during passive listening. Radiology 220, 202–207. doi: 10.1148/radiology.220.1.r01jl34202

PubMed Abstract | CrossRef Full Text | Google Scholar

Raichle, M. E., Mac Leod, 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

Raichle, M. E., and Snyder, A. Z. (2007). A default mode of brain function: a brief history of an evolving idea. NeuroImage 37, 1083–1090. doi: 10.1016/j.neuroimage.2007.02.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Reiman, E. M., Lane, R. D., Ahern, G. L., Schwartz, G. E., Davidson, R. J., Friston, K. J., et al. (1997). Neuroanatomical correlates of externally and internally generated human emotion. Am. J. Psychiatry 154, 918–925. doi: 10.1176/ajp.154.7.918

PubMed Abstract | CrossRef Full Text | Google Scholar

Rotenstein, L. S., Ramos, M. A., Torre, M., Segal, J. B., Peluso, M. J., Guille, C., et al. (2016). Prevalence of depression, depressive symptoms, and suicidal ideation Among medical students: a systematic review and meta-analysis. JAMA 316, 2214–2236. doi: 10.1001/jama.2016.17324

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

Sendi, M. S. E., Zendehrouh, E., Sui, J., Fu, Z., Zhi, D., Lv, L., et al. (2021). Abnormal dynamic functional network connectivity estimated from default mode network predicts symptom severity in major depressive disorder. Brain Connect. 11, 838–849. doi: 10.1089/brain.2020.0748

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheldon, S., Farb, N., Palombo, D. J., and Levine, B. (2016). Intrinsic medial temporal lobe connectivity relates to individual differences in episodic autobiographical remembering. Cortex 74, 206–216. doi: 10.1016/j.cortex.2015.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, Z., Cheng, Y., Yang, S., et al. (2016). Changes of grey matter volume in first-episode drug-naive adult major depressive disorder patients with different age-onset. Neuroimage Clin. 12, 492–498. doi: 10.1016/j.nicl.2016.08.016

CrossRef Full Text | Google Scholar

Shen, Z., Jiang, L., Yang, S., et al. (2017). Identify changes of brain regional homogeneity in early and later adult onset patients with first-episode depression using resting-state fMRI. PLoS One 12:e0184712. doi: 10.1371/journal.pone.0184712

CrossRef Full Text | Google Scholar

Shi, L. J., Zhou, H. Y., Wang, Y., Shen, Y. M., Fang, Y. M., He, Y. Q., et al. (2020). Altered empathy-related resting-state functional connectivity in adolescents with early-onset schizophrenia and autism spectrum disorders. Asian J. Psychiatr. 53:102167. doi: 10.1016/j.ajp.2020.102167

PubMed Abstract | CrossRef Full Text | Google Scholar

Sprengelmeyer, R., Steele, J. D., Mwangi, B., Kumar, P., Christmas, D., Milders, M., et al. (2011). The insular cortex and the neuroanatomy of major depression. J. Affect. Disord. 133, 120–127. doi: 10.1016/j.jad.2011.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Stratmann, M., Konrad, C., Kugel, H., et al. (2014). Insular and hippocampal gray matter volume reductions in patients with major depressive disorder. PLoS One 9:e102692. doi: 10.1371/journal.pone.0102692

CrossRef Full Text | Google Scholar

Sung, S. C., Wisniewski, S. R., Balasubramani, G. K., Zisook, S., Kurian, B., Warden, D., et al. (2013). Does early-onset chronic or recurrent major depression impact outcomes with antidepressant medications? A CO-MED trial report. Psychol. Med. 43, 945–960. doi: 10.1017/S0033291712001742

CrossRef Full Text | Google Scholar

Touroutoglou, A., Bliss-Moreau, E., Zhang, J., Mantini, D., Vanduffel, W., Dickerson, B. C., et al. (2016). A ventral salience network in the macaque brain. NeuroImage 132, 190–197. doi: 10.1016/j.neuroimage.2016.02.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Uddin, L. Q. (2015). Salience processing and insular cortical function and dysfunction. Nat. Rev. Neurosci. 16, 55–61. doi: 10.1038/nrn3857

CrossRef Full Text | Google Scholar

Vega, J. N., Taylor, W. D., Gandelman, J. A., Boyd, B. D., Newhouse, P. A., Shokouhi, S., et al. (2020). Persistent intrinsic functional network connectivity alterations in middle-aged and older women with remitted depression. Front. Psychiatry 11:62. doi: 10.3389/fpsyt.2020.00062

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X. L., Du, M. Y., Chen, T. L., et al. (2015). Neural correlates during working memory processing in major depressive disorder. Prog. Neuro-Psychopharmacol. Biol. Psychiatry 56, 101–108. doi: 10.1016/j.pnpbp.2014.08.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, D., Xue, S. W., Tan, Z., Wang, Y., Lian, Z., and Sun, Y. (2019). Altered hypothalamic functional connectivity patterns in major depressive disorder. Neuroreport 30, 1115–1120. doi: 10.1097/WNR.0000000000001335

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Q. Z., Li, D. M., Kuang, W. H., Zhang, T. J., Lui, S., Huang, X. Q., et al. (2011). Abnormal regional spontaneous neural activity in treatment-refractory depression revealed by resting-state fMRI. Hum. Brain Mapp. 32, 1290–1299. doi: 10.1002/hbm.21108

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, C. G., Chen, X., Li, L., Castellanos, F. X., Bai, T. J., Bo, Q. J., et al. (2019). Reduced default mode network functional connectivity in patients with recurrent major depressive disorder. Proc. Natl. Acad. Sci. U. S. A. 116, 9078–9083. doi: 10.1073/pnas.1900390116

PubMed Abstract | CrossRef Full Text | Google Scholar

Yao, X., Yin, Z., Liu, F., Wei, S., Zhou, Y., Jiang, X., et al. (2018). Shared and distinct regional homogeneity changes in bipolar and unipolar depression. Neurosci. Lett. 673, 28–32. doi: 10.1016/j.neulet.2018.02.033

PubMed Abstract | CrossRef Full Text | Google Scholar

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. Neuro-Psychopharmacol. Biol. Psychiatry 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

CrossRef Full Text | Google Scholar

Zhang, Z., Chen, Y., Wei, W., et al. (2021). Changes in regional homogeneity of medication-free major depressive disorder patients with different onset ages. Front. Psychiatry 12:713614. doi: 10.3389/fpsyt.2021.713614

CrossRef Full Text | Google Scholar