ORIGINAL RESEARCH article

Front. Nutr., 18 June 2025

Sec. Nutrition, Psychology and Brain Health

Volume 12 - 2025 | https://doi.org/10.3389/fnut.2025.1615978

Changes of functional brain network topology associated with nutritional indicator of patients with recurrent major depressive disorder

  • GL

    Guanghui Liu 1

  • YJ

    Yanli Jiang 2

  • CZ

    Chunpeng Zhou 1

  • XW

    Xiangxiang Wang 1

  • MD

    Min Dai 3

  • JZ

    Jinping Zhang 3

  • JY

    Jing Yu 4*

  • GZ

    Gangzhong Zhang 1*

  • 1. Department of Neurosurgery, Zhengzhou Seventh People’s Hospital, Zhengzhou, China

  • 2. Department of Radiology, Children’s Hospital of Nanjing Medical University, Nanjing, China

  • 3. Department of Neurology, Zhengzhou Seventh People’s Hospital, Zhengzhou, China

  • 4. Department of Radiology, Nanjing Brain Hospital, Affiliated Hospital of Nanjing Medical University, Nanjing, China

Article metrics

View details

1

Citations

1,8k

Views

412

Downloads

Abstract

Introduction:

Major depressive disorder (MDD) is generally categorized into first-episode MDD (fMDD) and recurrent MDD (rMDD). This study aimed to investigate the changes of brain network, as well as the relationships between relapses, brain regions, nutritional and metabolic indicators by resting-state functional magnetic resonance imaging (rs-fMRI) and graph theoretical analysis.

Methods:

Thirty-two fMDD patients, 32 rMDD patients, and 32 healthy controls (HCs) underwent rs-fMRI scanning. Graph theoretical analysis was applied to examine the nodal strength and nodal global efficiency of brain networks. In addition, the nutritional and metabolic indicators were acquired from all patients. The differences of demographic, clinical data and topological parameters between groups were compared. Moreover, the relationships between number of relapses, topological parameters of brain regions, nutritional and metabolic indicators were evaluated.

Results:

Patients with rMDD showed decreased level of brain-derived neurotrophic factor (BDNF) when compared with those with fMDD. Both fMDD and rMDD patients showed decreased nodal strength and global efficiency in the left amygdala compared to HCs. Additionally, rMDD patients exhibited more extensive network disruptions, including decreased nodal strength in the right middle frontal gyrus, left middle cingulate gyrus, left posterior cingulate gyrus, left hippocampus, and right amygdala, as well as decreased nodal global efficiency in multiple regions involved in emotional processing and cognitive control. Moreover, the number of episodes were negatively associated with the level of BDNF, nodal strength of right amygdala and nodal global efficiency of right amygdala of rMDD patients. The level of BDNF were positively related to nodal strength of right amygdala and nodal global efficiency of right amygdala of rMDD patients.

Conclusion:

Our findings revealed distinct patterns of brain network topology between fMDD and rMDD patients, with rMDD patients showing more widespread disruptions, which might be associated with greater number of relapses and worse level of neurological nutrition. These results suggested that recurrent depressive episodes might related to progressive disruptions in brain, particularly in regions involved in emotional processing and regulation.

1 Introduction

Major depressive disorder (MDD) represents one of the most prevalent and debilitating psychiatric conditions worldwide, characterized by persistent feelings of sadness, anhedonia, and various cognitive, emotional, and physical symptoms (1–3). Recent epidemiological studies indicate that MDD affects more than 264 million people globally, with a lifetime prevalence ranging from 3 to 21% in most developed countries (2, 4). The socioeconomic burden of MDD is substantial, with annual costs exceeding $326 billion in direct healthcare expenses and indirect costs such as lost productivity (5). The World Health Organization projects that by 2030, MDD will become the leading cause of disease burden globally, surpassing cardiovascular diseases and cancer (6). Moreover, approximately 35% of patients show inadequate response to first-line treatments, and up to 30% develop treatment-resistant depression, highlighting the urgent need for better understanding of disease mechanisms and more targeted therapeutic approaches (2, 7).

The clinical heterogeneity of MDD poses significant challenges for both research and treatment. MDD typically manifests in two distinct presentations: first-episode MDD (fMDD), representing the initial occurrence of depressive symptoms, and recurrent MDD (rMDD), characterized by multiple depressive episodes throughout the patient’s lifetime (8). The progression from first-episode to recurrent depression represents a critical transition in the disease course, with studies showing that approximately 60% of patients who recover from their first-episode experience recurrence, and this risk escalates to 90% after two episodes (9–11). This high recurrence rate suggests fundamental changes in brain function between first-episode and recurrent depression (12). Furthermore, each subsequent episode appears to increase the likelihood of future recurrences and may be associated with progressive cognitive decline and functional impairment. Despite the clinical significance of this progression, the underlying neural mechanisms driving the transition from first-episode to recurrent depression remain poorly understood (13, 14).

Recent advances in neuroscience have revolutionized our understanding of MDD pathophysiology, revealing that it extends far beyond simple monoamine deficiency to encompass complex alterations in neural circuits, neuroplasticity, and network dynamics (15–17). The emergence of network neuroscience has provided compelling evidence that psychiatric disorders, including MDD, can be conceptualized as disorders of brain network organization rather than isolated regional dysfunction (18). This network perspective has particular relevance for understanding the progression from fMDD to rMDD, as it may help explain how initial localized disruptions could propagate through neural circuits to create more widespread dysfunction over time (15, 18, 19). The network approach also aligns with the complex symptomatology of MDD, which includes disturbances in emotion, cognition, motivation, and physiological functions, suggesting dysfunction across multiple interacting neural systems rather than isolated brain regions.

Previous neuroimaging studies have identified widespread alterations in both structural and functional brain networks in MDD patients. These changes particularly affect regions involved in emotion processing (such as the amygdala and anterior cingulate cortex), cognitive control (such as the dorsolateral prefrontal cortex), and self-referential processing (such as the default mode network) (20–22). The limbic system, particularly the amygdala, has emerged as a critical hub in depression pathophysiology, showing consistent alterations across multiple studies (23). Neuroimaging studies have demonstrated that amygdala dysfunction correlates with symptom severity and treatment response, suggesting its potential role as a biomarker for disease progression. Additionally, altered connectivity between the amygdala and prefrontal regions has been implicated in emotional dysregulation, a core feature of MDD. However, most research has examined either fMDD or rMDD in isolation, leaving a significant gap in our understanding of how brain network organization differs between these two disease stages (24, 25).

The amygdala related pathways involved in emotion regulation mainly include the following (26–28): (1) thalamus-amygdala pathway, known as the short pathway: sensory information is transmitted from sensory organs to the thalamus, and is directly projected to the amygdala without going through fine processing in the cerebral cortex. (2) The cerebral cortex-thalamus-amygdala pathway, also known as the long pathway: sensory information is first transmitted from sensory organs to the thalamus, which then transmits the information to the cerebral cortex for detailed processing. (3) The amygdala-hippocampus pathway, which plays an important role in emotional learning and memory. The hippocampus also participates in the formation and storage of memory.

The occurrence of depression is also related to disturbances in nutrition of the central nervous system, as well as abnormal peripheral glucose and lipid metabolism (29). Brain neurogenesis is involved in the occurrence of depression through brain-derived neurotrophic factor (BDNF), an important neurotransmitter regulator (30). In the process of neurogenesis in depression, there are multiple regulatory signaling pathways, including the BDNF–TrkB signaling pathway. PI3K and Akt play important roles in neurogenesis under the action of BDNF. Almost all antidepressant drugs can enhance TrkB phosphorylation, activate this pathway to promote neurogenesis, and TrkB antagonists can block the antidepressant effect of drugs.

The purpose of this study was to investigate the differences in functional brain network topology between patients with fMDD and rMDD, and their relationships with nutritional, metabolic indicators using resting-state functional magnetic resonance imaging (rs-fMRI) and graph theory analysis. By examining network properties across the whole brain, we aimed to identify specific regional differences in nodal strength and global efficiency between these two patient groups and healthy controls (HCs). We hypothesized that both fMDD and rMDD patients would show altered network properties compared to HCs, and that rMDD patients might demonstrate more extensive network disruptions compared to fMDD patients, reflecting the cumulative impact of multiple depressive episodes on brain network organization. In addition, we hypothesized that the number of relapses might be associated with abnormal brain regions, impaired nutritional and metabolic indicators. Understanding the distinct patterns of brain network organization in fMDD versus rMDD might provide insights into the underlying mechanisms of disease progression and potentially guide the development of more targeted therapeutic strategies (such as neural regulation targeting bilateral amygdala or dorsolateral prefrontal cortex) for different stages of MDD. This knowledge could be particularly valuable for early intervention and prevention of recurrence in MDD patients, ultimately leading to improved treatment outcomes.

2 Materials and methods

2.1 Participants

This study was approved by the Medical Research Ethics Committee of Zhengzhou Seventh People’s Hospital and Nanjing Brain Hospital. Moreover, written informed consent was acquired from all subjects. A total of 32 fMDD patients and 32 rMDD patients were recruited. In addition, 32 right-handed HCs matched for age, gender and educational level were recruited.

The inclusion criteria were as follows: (1) 20–45 years of age; (2) at least 9 years of education; (3) right-handed; (4) Han Chinese. The diagnosis of MDD was made by two experienced psychiatrists according to the Diagnostic and Statistical Manual of Mental Disorders-fifth Edition (DSM-5) criteria. Both fMDD and rMDD patients had total scores of 17-item Hamilton Depression Rating Scale (HAMD-17) more than 17. The fMDD patients had no previous medication history. In addition, all HCs had total scores of HAMD-17 less than 7, and were screened through diagnostic interviews with the Structured Clinical Interview for DSM-5 Nonpatient Edition by two experienced psychiatrists, to exclude any current or history of psychiatric disorders.

The exclusion criteria were as follows: (1) uncontrolled or serious physical illness; (2) any other psychiatric, neurological, or somatic disorders; (3) history of any antidepressant medication 2 weeks before MRI scan; (4) history of head trauma or structural brain abnormalities identified by routine MRI; (5) abuse of alcohol or drugs; (6) female who were regnant, lactating or menstruating; (7) any contraindication for MRI scan.

The metabolic indicators including fasting blood glucose (FPG), glycated hemoglobin (HbA1c), total cholesterol (TC), triglycerides (TG), high-density lipoproteins (HDL) and low-density lipoproteins (LDL), as well as the nutritional indicators including BDNF and nerve growth factor (NGF), were acquired from all patients.

2.2 MRI data acquisition

MRI images were acquired using a 3.0 Tesla scanner (GE). All subjects were instructed to keep eyes closed, stay relaxed and awake, think of nothing in particular during MRI scanning. T1-weighted structural images were acquired with the following parameters: repetition time (TR) = 7.7 ms; echo time (TE) = 3.1 ms; slice thickness = 1 mm; field of view (FOV) = 256 × 256 mm2; matrix = 256 × 256; flip angle (FA) = 12°; slices = 160. The rs-MRI images were obtained with the following parameters: TR = 2000 ms; TE = 30 ms; slice thickness = 3.5 mm; FOV = 224 × 224 mm2; matrix = 80 × 80; FA = 90°; slices = 33; volumes = 240. Participants were excluded if structural abnormalities or artifacts were identified by structural MRI scans.

2.3 MRI data preprocessing

MRI data were preprocessed using the software of Data Processing Assistant for Resting-State fMRI (DPARSF) (31). The main steps were as follows: (1) removal of the first 10 time points; (2) slice timing; (3) head motion correction and subjects with translational motion exceeding 2 mm or rotational motion exceeding 2° were excluded; (4) spatial normalization into the standard Montreal Neurological Institute (MNI) space and re-sample into a voxel size of 3 mm × 3 mm × 3 mm; (5) spatial smoothing with a Gaussian kernel of 4 mm × 4 mm × 4 mm; (6) linear detrending; (7) temporal band-pass filtering at a frequency band of 0.01–0.1 Hz; (8) nuisance covariates regression (the head motion parameters, white matter and cerebrospinal fluid signals).

2.4 Construction of functional brain network

The graph theory analysis method defines a single brain region as a ‘node’ and defines the functional connectivity between brain regions as ‘edge’, simplifies the entire brain into a network diagram composed of nodes and edges. Firstly, the whole brain was divided into 90 nodes using the AAL template (including cortical and subcortical regions, excluding cerebellar regions) provided by the Montreal Institute (MNI) (32). Secondly, the Pearson correlation coefficients of the average time series between each pair of 90 brain regions were calculated to define the network edges.

2.5 Calculation of nodal topological parameters

Based on the GRETNA software package (33), the nodal topological parameters including nodal strength and nodal global efficiency were calculated (34). The sparsity was selected as 0.05–0.5, with an interval of 0.05. The area under curve (AUC) of these two metrics were calculated over all sparsity range.

2.6 Statistical analysis

The differences of demographic and clinical characteristics between patients and HCs were evaluated by one-way analysis of variance (ANOVA) and post hoc contrasts with the method of least-significant difference (LSD) (age and educational level), two-sample t-tests (HAMD sores, age of onset, metabolic and nutritional indicators), and Pearson chi-square test (gender). The statistical significances were set at p < 0.05.

In addition, the differences of nodal topological parameters including nodal strength and nodal global efficiency between patients and HCs were evaluated by ANOVA and post hoc contrasts with the method of LSD. The statistical significances were set at p < 0.05, and the multiple comparisons were corrected using the family-wise error (FWE) method with Bonferroni-corrected significance thresholds.

Moreover, the relationships between number of relapses, topological parameters of brain regions, nutritional and metabolic indicators were evaluated. The statistical significances were set at p < 0.05.

3 Results

3.1 Comparisons of demographic and clinical characteristics between patients and HCs

There were no significant differences in the gender, age and educational level among groups of fMDD, rMDD patients and HCs. In addition, no significant differences were found in the HAMD-17 scores, age of onset, metabolic indicators (FPG, HbA1c, TC, TG, HDL, LDL) and NGF level between fMDD and rMDD patients. However, patients with rMDD showed increased illness duration and decreased level of BDNF when compared with those with fMDD (Table 1).

Table 1

CharacteristicsfMDDrMDDHCsχ2/F/tp
Gender (male/female)15/1715/1714/180.0840.96
Age (years)33.06 ± 7.1634.06 ± 7.5931.47 ± 7.980.950.39
Educational level (years)13.66 ± 2.7513.84 ± 2.5213.19 ± 2.710.520.60
HAMD-17 (scores)25.09 ± 4.9226.47 ± 5.61-−1.040.30
Age of onset (years)33.06 ± 7.1632.19 ± 7.38-0.480.63
Number of episodes-2.94 ± 0.84---
Illness duration (months)3.84 ± 1.1413.47 ± 5.04-−10.55<0.01
Interval between episodes (months)-7.69 ± 3.91---
Metabolic indicators
FPG (mmol/L)5.08 ± 0.715.03 ± 0.69-0.240.81
HbA1c (%)5.12 ± 0.574.97 ± 0.63-1.030.31
TC (mmol/L)4.06 ± 0.763.94 ± 0.81-0.620.54
TG (mmol/L)1.08 ± 0.341.10 ± 0.37-−0.200.85
HDL (mmol/L)1.32 ± 0.121.36 ± 0.12-−1.190.24
LDL (mmol/L)2.61 ± 0.311.36 ± 0.32-1.230.23
Nutritional biomarkers
BDNF (ng/mL)30.15 ± 4.8826.97 ± 5.79-2.370.02
NGF (ng/L)9.49 ± 3.308.68 ± 2.93-1.030.31

The demographic and clinical characteristics of fMDD, rMDD patients and HCs.

MDD, major depressive disorder; fMDD, first-episode MDD patients; rMDD, recurrent MDD patients; HCs, healthy controls; HAMD-17, 17-item Hamilton Depression Rating Scale; FPG, fasting blood glucose; HbA1c, glycated hemoglobin; TC, total cholesterol; TG: triglycerides; HDL, high-density lipoproteins; LDL, low-density lipoproteins; BDNF, brain-derived neurotrophic factor; NGF, nerve growth factor. p values were obtained by one-way analysis of variance (ANOVA) and post hoc contrasts with the method of least-significant difference (LSD) (age and educational level), two-sample t-tests (HAMD sores, age of onset, metabolic and nutritional indicators), and Pearson chi-square test (gender). The statistical significances were set at p < 0.05.

3.2 Comparisons of nodal topological parameters between patients and HCs

Compared with HCs, rMDD patients showed decreased nodal strength in the right middle frontal gyrus, left middle cingulate gyrus, left posterior cingulate gyrus, left hippocampus, left amygdala and right amygdala. In addition, rMDD patients exhibited decreased nodal global efficiency in the right precental gyrus, right dorsolateral superior frontal gyrus, left orbital part of superior frontal gyrus, right middle frontal gyrus, left middle cingulate gyrus, right middle cingulate gyrus, left posterior cingulate gyrus, right posterior cingulate gyrus, left hippocampus, left amygdala and right amygdala (Table 2; Figures 1, 2).

Table 2

Brain regionsfMDDrMDDHCsFp
Nodal strength
Right middle frontal gyrus6.16 ± 2.797.37 ± 4.4510.29 ± 4.908.370.00045a
Left middle cingulate gyrus4.86 ± 2.347.48 ± 3.9310.09 ± 4.9514.450.000003a
Left posterior cingulate gyrus5.39 ± 2.418.08 ± 4.1710.37 ± 5.2311.800.000027a
Left hippocampus5.19 ± 2.976.28 ± 4.219.26 ± 4.778.640.00036a
Left amygdala6.61 ± 2.218.93 ± 1.8111.95 ± 3.3635.383.75E-12a,b,c
Right amygdala6.61 ± 2.639.25 ± 3.2711.69 ± 3.3321.651.91E-8a,c
Nodal global efficiency
Right precental gyrus0.18 ± 0.030.20 ± 0.040.21 ± 0.048.780.00032a
Right dorsolateral superior frontal gyrus0.17 ± 0.020.18 ± 0.040.21 ± 0.048.490.00041a
Left orbital part of superior frontal gyrus0.14 ± 0.030.15 ± 0.040.17 ± 0.048.580.00038a
Right middle frontal gyrus0.16 ± 0.030.17 ± 0.050.21 ± 0.059.080.00025a
Left middle cingulate gyrus0.15 ± 0.030.18 ± 0.040.21 ± 0.0516.686.47E-7a
Right middle cingulate gyrus0.14 ± 0.030.16 ± 0.040.18 ± 0.048.570.00038a
Left posterior cingulate gyrus0.15 ± 0.030.18 ± 0.050.21 ± 0.0513.730.000006a
Right posterior cingulate gyrus0.16 ± 0.030.17 ± 0.040.19 ± 0.049.520.00017a
Left hippocampus0.15 ± 0.040.17 ± 0.050.20 ± 0.058.770.00032a
Left amygdala0.16 ± 0.030.19 ± 0.030.22 ± 0.0429.231.42E-10a,b,c
Right amygdala0.16 ± 0.030.20 ± 0.040.22 ± 0.0422.997.70E-9a,c

Changed nodal topological parameters of functional brain networks in fMDD and rMDD patients.

MDD, major depressive disorder; fMDD, first-episode MDD patients; rMDD, recurrent MDD patients; HCs, healthy controls. adecreased nodal topological parameters in rMDD patients when compared with HCs; bdecreased nodal topological parameters in fMDD patients when compared with HCs; cdecreased nodal topological parameters in rMDD patients when compared with fMDD patients. p values were obtained by one-way analysis of variance (ANOVA) and post hoc contrasts with the method of least-significant difference (LSD). The statistical significances were set at P < 0.05, and the multiple comparisons were corrected using the family-wise error (FWE) method with Bonferroni-corrected significance thresholds.

Figure 1

Figure 2

Compared with HCs, fMDD patients demonstrated decreased nodal strength in the left amygdala and decreased nodal global efficiency in the left amygdala (Table 2; Figures 1, 2).

Compared with fMDD patients, rMDD patients demonstrated decreased nodal strength in the right amygdala and decreased nodal global efficiency in the right amygdala (Table 2; Figures 1, 2).

3.3 Relationships between abnormal brain regions and clinical characteristics

The number of episodes of MDD were negatively associated with the level of BDNF (r = −0.64; p < 0.01), nodal strength of the right amygdala (r = −0.63; p < 0.01) and nodal global efficiency of the right amygdala (r = −0.56; p < 0.01) of rMDD patients. In addition, the level of BDNF were positively related to the nodal strength of the right amygdala (r = 0.48; p < 0.01) and nodal global efficiency of the right amygdala (r = 0.42; p < 0.01) in rMDD patients (Figure 3).

Figure 3

4 Discussion

In this study, we investigated the differences in functional brain network topology between patients with fMDD and rMDD using rs-fMRI and graph theory analysis. Our findings revealed that both fMDD and rMDD patients exhibited altered brain network properties when compared to HCs, with rMDD patients showing more extensive disruptions in network organization, particularly in the amygdala and other emotion-related regions. In addition, the number of episodes were negatively associated with BDNF level, nodal parameters of right amygdala while BDNF level were positively related to nodal parameters of right amygdala in rMDD patients. These findings suggested that rMDD patients had more widespread disruptions, which might be related to greater number of relapses and worse level of neurological nutrition. All these demonstrated that that recurrent depressive episodes might related to progressive disruptions in brain.

Firstly, in comparison with HCs, both fMDD and rMDD patients showed decreased nodal strength and global efficiency in the left amygdala (the affected side of the amygdala), a critical region in emotional processing (35). The decreased nodal strength indicated reduced functional connectivity between the amygdala and other brain regions, while reduced global efficiency suggested compromised information processing capability within the amygdala-centered network (36, 37). These alterations aligned with previous studies demonstrating abnormal amygdala function in MDD, suggesting that disrupted amygdala connectivity might be a core feature of depression regardless of episode status (38, 39).

Notably, rMDD patients exhibited more extensive network disruptions compared to both HCs and fMDD patients. Beyond the amygdala, rMDD patients showed decreased nodal strength in multiple regions including the right middle frontal gyrus, left middle cingulate gyrus, left posterior cingulate gyrus, and left hippocampus. These regions are key components of the emotional processing and cognitive control networks (40, 41). The middle frontal gyrus is involved in emotional regulation and executive function, while the cingulate cortex plays a crucial role in emotional awareness and cognitive control (40, 41). Additionally, the reduced global efficiency in the right precentral gyrus and right dorsolateral superior frontal gyrus suggested compromised top-down control circuits essential for both motor control and executive function (3). The decreased efficiency in the left orbital part of superior frontal gyrus, a region crucial for emotional processing and decision-making, further indicated disruption of frontal-limbic circuits that were fundamental for mood regulation (42). The bilateral involvement of middle and posterior cingulate regions indicated a more systemic disruption of the emotion regulation network, potentially explaining the increased difficulty in emotional regulation observed in rMDD patients (43). The hippocampus, working in concert with the amygdala, is essential for emotional memory processing (44). The widespread network disruptions in rMDD patients suggested that recurrent episodes might be associated with more pervasive alterations in brain network organization.

Our findings highlighted a particularly interesting hemispheric difference in amygdala dysfunction between fMDD and rMDD patients. While fMDD patients showed only left amygdala (the affected side of the amygdala) alterations, rMDD patients demonstrated bilateral amygdala abnormalities. This hemispheric difference might reflect the progression of network dysfunction from first-episode to recurrence. The bilateral nature of amygdala dysfunction in rMDD patients might explain their potentially higher risk of recurrence and more severe clinical presentations (45, 46). This pattern of progressive network disruption suggested that the cumulative impact of multiple depressive episodes might be associated with more severe disruptions in the brain’s emotional processing networks.

The reduced nodal global efficiency observed in rMDD patients across multiple regions suggested compromised network integration and information processing. These more extensive disruptions in rMDD patients indicated that the chronic nature of the illness might exacerbate these deficits, which might be associated with a more compromised brain network organization. This finding had important implications for understanding disease progression and treatment approaches, as it suggested that early intervention might be crucial in preventing the development of more widespread network dysfunction.

The ‘Neurotrophic Hypothesis of Depression’ suggests that the interruption of neurotrophic support is a key mechanism for MDD related synaptic and brain related functional changes (16). Neurotrophic factors are responsible for the formation, support, and plasticity of neural networks. Numerous studies have shown that the levels of neurotrophic factors in the blood are reduced in patients with persistent and recurrent depression, as well as in animal models of depression. In this study, BDNF level were negatively associated with the number of episodes, and positively related to nodal parameters of right amygdala in rMDD patients. These findings suggested that rMDD patients had worse level of neurological nutrition, which might related to progressive disruptions in brain.

These findings had several potential clinical implications. The distinct patterns of network disruption between fMDD and rMDD suggested that different therapeutic approaches might be needed for these two groups. For example, interventions targeting bilateral amygdala function might be more relevant for rMDD patients, while focused left amygdala (the affected side of the amygdala) interventions might be more appropriate for fMDD patients. Treatment strategies that addressed both the emotional and cognitive aspects of the disorder, such as cognitive behavioral therapy (CBT) or neuromodulation techniques like transcranial magnetic stimulation (TMS), might be particularly beneficial for patients with recurrent depression. These therapies could be tailored to target the specific brain regions with altered network properties, offering a more effective approach for managing both acute and recurrent depressive episodes (47–49).

Several limitations of this study should be noted. Firstly, the cross-sectional design prevents us from determining whether the observed network differences were causes or consequences of recurrence. Longitudinal studies were needed to clarify the temporal evolution of these network alterations. Secondly, while the current study focused on rs-fMRI (49), future research could incorporate task-based fMRI or other neuroimaging techniques such as diffusion tensor imaging (DTI) to further explore the structural connectivity changes that might underlie the observed functional disruptions. Finally, exploring the relationship between brain network alterations and clinical factors such as different symptoms and response to treatment would provide a more comprehensive understanding of the functional significance of these changes.

5 Conclusion

In conclusion, this study revealed distinct patterns of functional brain network topology between patients with fMDD and rMDD using rs-fMRI and graph theoretical analysis. Our findings demonstrated that both fMDD and rMDD patients showed altered network properties in the left amygdala (the affected side of the amygdala) compared to HCs. More importantly, rMDD patients exhibited more extensive network disruptions, including bilateral amygdala abnormalities and decreased nodal strength and global efficiency in multiple emotion-related brain regions. Greater number of relapses were associated with decreased BDNF level and impaired function of right amygdala in rMDD patients while decreased BDNF level were related to impaired function of right amygdala in rMDD patients. These findings suggested that recurrent depressive episodes might be associated with progressive disruptions in brain, particularly in regions involved in emotional processing and regulation, which provided important implications for early intervention and personalized treatment strategies in MDD patients.

Statements

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 studies involving humans were approved by the Medical Research Ethics Committee of Zhengzhou Seventh People’s Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

GL: Software, Validation, Data curation, Investigation, Resources, Supervision, Writing – review & editing, Project administration, Visualization, Conceptualization, Writing – original draft, Methodology, Formal analysis, Funding acquisition. YJ: Software, Investigation, Writing – original draft, Resources, Funding acquisition, Visualization, Validation, Formal analysis, Data curation, Supervision, Project administration, Writing – review & editing, Methodology, Conceptualization. CZ: Investigation, Writing – original draft, Resources, Visualization, Funding acquisition, Software, Formal analysis, Validation, Conceptualization, Data curation, Supervision, Writing – review & editing, Project administration, Methodology. XW: Methodology, Data curation, Conceptualization, Investigation, Validation, Supervision, Writing – review & editing, Visualization, Software, Writing – original draft, Formal analysis, Resources, Project administration, Funding acquisition. MD: Formal analysis, Visualization, Data curation, Resources, Validation, Project administration, Writing – review & editing, Software, Methodology, Writing – original draft, Funding acquisition, Conceptualization, Investigation, Supervision. JZ: Writing – original draft, Formal analysis, Software, Visualization, Funding acquisition, Resources, Methodology, Project administration, Supervision, Conceptualization, Investigation, Data curation, Writing – review & editing, Validation. JY: Project administration, Methodology, Validation, Formal analysis, Supervision, Data curation, Writing – original draft, Software, Conceptualization, Writing – review & editing, Resources, Investigation, Funding acquisition, Visualization. GZ: Writing – original draft, Formal analysis, Visualization, Resources, Funding acquisition, Writing – review & editing, Project administration, Methodology, Investigation, Validation, Conceptualization, Software, Data curation, Supervision.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

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.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

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.

References

  • 1.

    LiuZLiuYZhaoXZhangHFengTPangJet al. Correlation between aβ(1-42), Dnmt3a2, urinary Ad7c-Ntp and cognitive dysfunction in first-episode and recurrent Mdd: a case-control study. Indian J Psychiatry. (2022) 64:5606. doi: 10.4103/indianjpsychiatry.indianjpsychiatry_111_22

  • 2.

    GuanJSunYFanYLiangJLiuCYuHet al. Effects and neural mechanisms of different physical activity on major depressive disorder based on cerebral multimodality monitoring: a narrative review. Front Hum Neurosci. (2024) 18:1406670. doi: 10.3389/fnhum.2024.1406670

  • 3.

    ZhuYHuangTLiRYangQZhaoCYangMet al. Distinct resting-state effective connectivity of large-scale networks in first-episode and recurrent major depression disorder: evidence from the rest-meta-MDD consortium. Front Neurosci. (2023) 17:1308551. doi: 10.3389/fnins.2023.1308551

  • 4.

    Gutiérrez-RojasLPorras-SegoviaADunneHAndrade-GonzálezNCervillaJA. Prevalence and correlates of major depressive disorder: a systematic review. Braz J Psychiatry. (2020) 42:65772. doi: 10.1590/1516-4446-2020-0650

  • 5.

    GreenbergPFournierASisitskyTGreenbergPEFournierA-ASimesMet al. The economic burden of adults with major depressive disorder in the United States (2010 and 2018). PharmacoEconomics. (2021) 39:65365. doi: 10.1007/s40273-021-01019-4

  • 6.

    World Health Organization. Global burden of mental disorders and the need for a comprehensive, coordinated response from health and social sectors at the country level: report by the secretariat. Geneva: World Health Organization (2012).

  • 7.

    ChenYJiaLGaoWWuCMuQFangZet al. Alterations of brainstem volume in patients with first-episode and recurrent major depressive disorder. BMC Psychiatry. (2023) 23:687. doi: 10.1186/s12888-023-05146-4

  • 8.

    WangKHuYYanCLiMWuYQiuJet al. Brain structural abnormalities in adult major depressive disorder revealed by voxel- and source-based morphometry: evidence from the rest-meta-MDD consortium. Psychol Med. (2023) 53:367282. doi: 10.1017/S0033291722000320

  • 9.

    DienerCKuehnerCBrusniakWUblBWessaMFlorH. A meta-analysis of neurofunctional imaging studies of emotion and cognition in major depression. NeuroImage. (2012) 61:67785. doi: 10.1016/j.neuroimage.2012.04.005

  • 10.

    LyeMTeyYTorYLyeMSTeyYYTorYSet al. Predictors of recurrence of major depressive disorder. PLoS One. (2020) 15:e0230363. doi: 10.1371/journal.pone.0230363

  • 11.

    BurcusaSIaconoW. Risk for recurrence in depression. Clin Psychol Rev. (2007) 27:95985. doi: 10.1016/j.cpr.2007.02.005

  • 12.

    SunFYanJPangJSongMWangMHuangTet al. Distinct effects of first-episode and recurrent major depressive disorder on brain changes related to suicidal ideation: evidence from the rest-meta-mdd project. J Affect Disord. (2024) 351:47280. doi: 10.1016/j.jad.2024.01.213

  • 13.

    JappeLKlimes-DouganBCrowS. An fMRI study of emotional face processing in adolescent major depression in partial fulfillment of the requirements for the degree of master of arts emotion face processing in adolescent depression emotion face processing in adolescent depression. Minneapolis, MN: University of Minnesota (2013).

  • 14.

    ChenDLiZLuoY. Alterations in cognitive and emotional brain networks related to depression. Chin Sci Bull. (2017) 62:24929. doi: 10.1360/N972016-00892

  • 15.

    YeMYangTQingPLeiXQiuJLiuG. Changes of functional brain networks in major depressive disorder: a graph theoretical analysis of resting-state fMRI. PLoS One. (2015) 10:e0133775. doi: 10.1371/journal.pone.0133775

  • 16.

    FriesGSaldanaVFinnsteinJFriesGRSaldanaVAReinT. Molecular pathways of major depressive disorder converge on the synapse. Mol Psychiatry. (2023) 28:28497. doi: 10.1038/s41380-022-01806-1

  • 17.

    CuiLLiSWangSWuXLiuYYuWet al. Major depressive disorder: hypothesis, mechanism, prevention and treatment. Signal Transduct Target Ther. (2024) 9:30. doi: 10.1038/s41392-024-01738-y

  • 18.

    AgziyartEAbbasianKMakoueiSAgziyartEAMohammadiSB. Investigating changes of functional brain networks in major depressive disorder by graph theoretical analysis of resting-state fMRI. Psychiatry Res Neuroimaging. (2024) 344:111880. doi: 10.1016/j.pscychresns.2024.111880

  • 19.

    ZhiDMaXLvLKeQYangYYangXet al. Abnormal dynamic functional network connectivity and graph theoretical analysis in major depressive disorder. Annu Int Conf IEEE Eng Med Biol Soc. (2018) 2018:55861. doi: 10.1109/EMBC.2018.8512340

  • 20.

    ZhangFPengWSweeneyJSweeneyJAJiaZ‐YGongQ‐Y. Brain structure alterations in depression: psychoradiological evidence. CNS Neurosci Ther. (2018) 24:9941003. doi: 10.1111/cns.12835

  • 21.

    IbrahimHKulikovaALyHIbrahimHMRushAJSherwood BrownE. Anterior cingulate cortex in individuals with depressive symptoms: a structural MRI study. Psychiatry Res Neuroimaging. (2022) 319:111420. doi: 10.1016/j.pscychresns.2021.111420

  • 22.

    ShelineYBarchDPriceJShelineYIBarchDMPriceJLet al. The default mode network and self-referential processes in depression. Proc Natl Acad Sci USA. (2009) 106:19427. doi: 10.1073/pnas.0812686106

  • 23.

    ConnollyCHoTBlomEConnollyCGHoTCBlomEHet al. Resting-state functional connectivity of the amygdala and longitudinal changes in depression severity in adolescent depression. J Affect Disord. (2017) 207:8694. doi: 10.1016/j.jad.2016.09.026

  • 24.

    WangLZhaoYEdmistonEKWomerFYZhangRZhaoPet al. Structural and functional abnormities of amygdala and prefrontal cortex in major depressive disorder with suicide attempts. Front Psych. (2019) 10:923. doi: 10.3389/fpsyt.2019.00923

  • 25.

    SunJDuZMaYChenLWangZGuoCet al. Altered functional connectivity in first-episode and recurrent depression: a resting-state functional magnetic resonance imaging study[J]. Front Neurol. (2022) 13:922207. doi: 10.3389/fneur.2022.922207

  • 26.

    SergerieKChocholCArmonyJ. The role of the amygdala in emotional processing: a quantitative meta-analysis of functional neuroimaging studies. Neurosci Biobehav Rev. (2008) 32:81130. doi: 10.1016/j.neubiorev.2007.12.002

  • 27.

    LabuschagneIJfDGraceSMizziSHenryJDPetersCet al. Specialization of amygdala subregions in emotion processing. Hum Brain Mapp. (2024) 45:e26673. doi: 10.1002/hbm.26673

  • 28.

    SeoHOemischM. Decoding emotion: the amygdala-prefrontal cortex pathway for emotion regulation of children. Biol Psychiatry. (2020) 88:5179. doi: 10.1016/j.biopsych.2020.07.017

  • 29.

    MehdiSMACostaAPSvobCPanLDartoraWJTalatiAet al. Depression and cognition are associated with lipid dysregulation in both a multigenerational study of depression and the National Health and nutrition examination survey. Transl Psychiatry. (2024) 14:142. doi: 10.1038/s41398-024-02847-6

  • 30.

    HanHYaoJWuJMaoSPanHQvLet al. Implications of neurogenesis in depression through BDNF: rodent models, regulatory pathways, gut microbiota, and potential therapy. Mol Psychiatry. (2025). doi: 10.1038/s41380-025-03044-7 [Epub ahead of print].

  • 31.

    Chao-GanYYu-FengZ. DPARSF: a Matlab toolbox for "pipeline" data analysis of resting-state fMRI. Front Syst Neurosci. (2010) 4:–13. doi: 10.3389/fnsys.2010.00013

  • 32.

    Tzourio-MazoyerNLandeauBPapathanassiouDCrivelloFEtardODelcroixNet al. Automated anatomical labeling of activations in Spm using a macroscopic anatomical parcellation of the Mni Mri single-subject brain. NeuroImage. (2002) 15:27389. doi: 10.1006/nimg.2001.0978

  • 33.

    WangJWangXXiaMLiaoXEvansAHeY. Gretna: a graph theoretical network analysis toolbox for imaging connectomics. Front Hum Neurosci. (2015) 9:386. doi: 10.3389/fnhum.2015.00386

  • 34.

    RubinovMSpornsO. Complex network measures of brain connectivity: uses and interpretations. NeuroImage. (2010) 52:105969. doi: 10.1016/j.neuroimage.2009.10.003

  • 35.

    WuFTuZSunJet al. Abnormal functional and structural connectivity of amygdala-prefrontal circuit in first-episode adolescent depression: a combined fMRI and DTI study. Front Psych. (2019) 10:983. doi: 10.3389/fpsyt.2019.00983

  • 36.

    NachshonOFarahRHorowitz-KrausT. Decreased functional connectivity between the left amygdala and frontal regions interferes with reading, emotional, and executive functions in children with reading difficulties. Front Hum Neurosci. (2020) 14:104. doi: 10.3389/fnhum.2020.00104

  • 37.

    AloiJMcCuskerMCLewBJSchantellMEastmanJAFrenzelMRet al. Altered amygdala-cortical connectivity in individuals with cannabis use disorder. J Psychopharmacol. (2021) 35:136574. doi: 10.1177/02698811211054163

  • 38.

    ShengFWangYLiRLiXChenXZhangZet al. Altered effective connectivity among face-processing systems in major depressive disorder. J Psychiatry Neurosci. (2024) 49:E145e156. doi: 10.1503/jpn.230123

  • 39.

    UedaSTakemoto-KimuraS. Exploring the molecular and neuronal bases involved in central amygdala-dependent control of emotional behaviors. Nihon Yakurigaku Zasshi. (2024) 159:31620. doi: 10.1254/fpj.23052

  • 40.

    OuYNiXGaoXYuYZhangYWangYet al. Structural and functional changes of anterior cingulate cortex subregions in migraine without aura: relationships with pain sensation and pain emotion. Cereb Cortex. (2024) 34:bhae040. doi: 10.1093/cercor/bhae040

  • 41.

    WangYZhuWXiaoMZhangQZhaoYZhangHet al. Hostile attribution bias mediates the relationship between structural variations in the left middle frontal gyrus and trait angry rumination. Front Psychol. (2018) 9:526. doi: 10.3389/fpsyg.2018.00526

  • 42.

    DvorakJHilkeMTrettinMWenzlerSHagenMGhirmaiNet al. Aberrant brain network topology in fronto-limbic circuitry differentiates euthymic bipolar disorder from recurrent major depressive disorder. Brain Behav. (2019) 9:e01257. doi: 10.1002/brb3.1257

  • 43.

    ShengWCuiQJiangKChenYTangQWangCet al. Individual variation in brain network topology is linked to course of illness in major depressive disorder. Cereb Cortex. (2022) 32:530110. doi: 10.1093/cercor/bhac015

  • 44.

    KhatriUKwonG. Alzheimer's disease diagnosis and biomarker analysis using resting-state functional MRI functional brain network with multi-measures features and hippocampal subfield and amygdala volume of structural MRI. Front Aging Neurosci. (2022) 14:818871. doi: 10.3389/fnagi.2022.818871

  • 45.

    LiuHWangCLanXLiWZhangFFuLet al. Functional connectivity of the amygdala and the antidepressant and antisuicidal effects of repeated ketamine infusions in major depressive disorder. Front Neurosci. (2023) 17:1123797. doi: 10.3389/fnins.2023.1123797

  • 46.

    JacobsRHBarbaAGowinsJRKlumppHJenkinsLMMickeyBJet al. Decoupling of the amygdala to other salience network regions in adolescent-onset recurrent major depressive disorder. Psychol Med. (2016) 46:105567. doi: 10.1017/S0033291715002615

  • 47.

    AustSBrakemeierESpiesJBrakemeierE-LHerrera-MelendezALKaiserTet al. Efficacy of augmentation of cognitive behavioral therapy with transcranial direct current stimulation for depression: a randomized clinical trial. JAMA Psychiatr. (2022) 79:52837. doi: 10.1001/jamapsychiatry.2022.0696

  • 48.

    KönigPZwikyEKüttnerAUhligMRedlichR. Brain functional effects of cognitive behavioral therapy for depression: A systematic review of task-based fMRI studies. J Affect Disord. (2025) 368:87287. doi: 10.1016/j.jad.2024.09.084

  • 49.

    KhameneASPourSBHeidariAFarahNParvinE. Comparison of the efficacy of trans-cranial magnetic stimulation and cognitive-behavioral therapy on depression. Iran J Educ Sociol. (2019) 2:1827. doi: 10.29252/ijes.2.2.182

Summary

Keywords

major depressive disorder, resting-state functional magnetic resonance imaging, functional brain network, recurrence, brain-derived neurotrophic factor

Citation

Liu G, Jiang Y, Zhou C, Wang X, Dai M, Zhang J, Yu J and Zhang G (2025) Changes of functional brain network topology associated with nutritional indicator of patients with recurrent major depressive disorder. Front. Nutr. 12:1615978. doi: 10.3389/fnut.2025.1615978

Received

22 April 2025

Accepted

04 June 2025

Published

18 June 2025

Volume

12 - 2025

Edited by

Chuanfeng Tang, Nanjing University of Chinese Medicine, China

Reviewed by

Yifan Wang, Nanjing University of Chinese Medicine, China

Hong Ding, Anhui Medical University, China

Updates

Copyright

*Correspondence: Jing Yu, Gangzhong Zhang,

†These authors share first authorship

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics