Abstract

Background. The role of childhood maltreatment (CM) is believed to be crucial in the aberrant function of reward circuit in adolescent-onset major depressive disorder (AO-MDD). Nevertheless, the impact of abnormalities in the GP-based reward circuit on the association between CM and the severity of AO-MDD remains largely unknown. Methods. The GP-based resting-state functional connectivity (RSFC) was analyzed in a sample of 75 patients with AO-MDD and 80 healthy controls in order to identify potential abnormalities in the GP-based reward circuit in AO-MDD patients. Furthermore, we investigated the possible associations between aberrant GP-based reward circuit functioning, CM and its subtypes (namely, childhood abuse and childhood neglect), and the severity of AO-MDD. Results. Compared to the healthy control group, patients with AO-MDD demonstrated a reduction in RSFC between the left posterior GP and the right dorsomedial prefrontal cortex (DMPFC). Our moderation analysis revealed that the abnormal RSFC between the posterior GP and DMPFC had a moderating effect on the relationship between CM and the severity of AO-MDD. Furthermore, upon further interaction decomposition, we observed a positive correlation between CM and AO-MDD severity exclusively in patients with AO-MDD who exhibited lower RSFC between the posterior GP and DMPFC. For AO-MDD patients with higher RSFC between posterior GP and DMPFC, the relationship between CM and AO-MDD severity was not discernible. Conclusions. Our findings underscore the crucial role of anomalies in the reward circuit in AO-MDD and furnish novel leads for probing the relationship among CM, malfunctioning of the reward circuit, and AO-MDD.

1. Introduction

In the past few years, there has been a significant surge in the incidence of adolescent-onset major depressive disorder (AO-MDD), with a 12-month depression prevalence of about 13% [1]. Adolescence is a crucial stage of quick emotional and cognitive growth, and AO-MDD can significantly impede neurodevelopment and social functioning (such as academic progress and interpersonal relationships) [2]. Given its strong correlation with dropping out of school, nonsuicidal self-harm, persistent ailment, and recurrent depressive episodes in adulthood, AO-MDD has emerged as a growing societal issue [3, 4].

Childhood maltreatment (CM), which includes emotional, physical, and sexual abuse, as well as emotional and physical neglect, is a prevalent and potent risk element in the development of AO-MDD. Previous research has demonstrated that CM impacts over 35% of the worldwide population [5] and 45.6% of depressed patients have a history of CM [6]. AO-MDD patients exposed to CM often display an earlier onset of symptoms, more cooccurring disorders, and an inadequate response to clinical treatment [6, 7]. However, meta-analytic evidence suggests that traumatic childhood experiences are preventable and that AO-MDD patients could benefit from active treatment [8]. Consequently, as a disruptive yet positively modifiable aspect, CM should be accorded high priority in the intervention strategies of AO-MDD [8].

Clinical investigations have unveiled variegated responses to traumatic experiences among patients with AO-MDD [9, 10], indicating the existence of a subgroup of AO-MDD patients who are highly susceptible to the adverse effects of traumatic experiences, i.e., highly sensitive to CM. The neurobiological mechanisms that underlie hypersensitivity to CM in AO-MDD remain unclear, but certain clues suggest that it may be associated with aberrant resting-state functional connectivity (RSFC) of the reward circuit [11, 12]. The reward circuit is involved in the processing of rewarding effects (e.g., pleasure and euphoria), and the upregulation of dopamine system function in the reward circuit increases euphoric activity [13]. The core of the reward circuit is comprised of the globus pallidus (GP) and the nucleus accumbens, which receive signal input from the prefrontal cortex, hippocampus, and amygdala [14, 15]. While measuring the electrophysiological activity of the GP and nucleus accumbens in rats performing reward-based tasks, Ottenheimer et al. observed a greater number of reward-sensitive neurons within the GP compared to those in the nucleus accumbens [16]. Furthermore, a preference-related signal that flexibly reports the relative value of reward outcomes across multiple conditions was discovered in the GP [16]. These findings underscore the pivotal role of the GP in the reward circuit, suggesting that the GP is not simply a relay station for downstream signals from the nucleus accumbens but rather a core brain region within the reward circuit [16].

Previous research has demonstrated the adverse effects of CM on the reward circuit, specifically in the GP region. For example, Dillon et al. have shown that individuals who have been maltreated exhibit increased anhedonia and decreased responsivity of the GP towards rewarding stimuli [17]. However, it is important to note that the GP is not a homogenous structure, rather it is comprised of functionally and structurally heterogeneous nuclei. In a recent study, Tian et al. have constructed a novel set of brain atlases that apply to subcortical nuclei, based on a sample of over 1,000 healthy adults [18]. This pioneering work deviates from previous approaches that relied on structural magnetic resonance imaging (MRI) to map brain anatomy, as it utilizes functional MRI to elucidate the intricate topographical organization of the human subcortex, culminating in one of the most detailed gradient atlases of the subcortex to date [18]. According to the brain atlases derived by Tian et al., the GP is segmented into four subregions, which are the bilateral anterior globus pallidus (aGP) and the bilateral posterior globus pallidus (pGP) [18]. The aGP primarily serves motivational and cognitive functions [19], while the pGP is mainly involved in regulating motor functions, particularly the precision grip movement parameters [20].

The RSFC analysis represents a fundamental method for characterizing interregional brain activity, as it allows the examination of temporal correlations of blood oxygen level-dependent signals among different brain regions [21]. RSFC analysis provides more objective evidence, shedding light on significant pathophysiological processes in depression [22]. In recent years, RSFC has been employed to investigate functional abnormalities in the brains of individuals with MDD with CM [23, 24]. In the present study, we selected four subregions of the GP in the bilateral hemispheres as regions of interest (ROI) for RSFC analysis, with the aim of exploring abnormalities in the reward circuit of patients with AO-MDD. Specifically, we aimed to investigate the relationship between CM, altered function of the reward circuit, and the outcomes of AO-MDD. We generated the following hypothesis: (i) compared to healthy adolescents, patients with AO-MDD exhibit anomalous reward circuit function that is centered on the GP, and (ii) in patients with AO-MDD, altered reward circuit function may play a critical role in explaining the relationship between CM and the severity of depression. Our study contributes new evidence to confirm the presence of abnormal reward circuit in AO-MDD and provides additional neuroimaging clues that may aid in identifying AO-MDD patients who are highly sensitive to traumatic experiences.

2. Methods

2.1. Participants

Seventy-five patients diagnosed with AO-MDD and 80 healthy controls (HCs) were included in the baseline data of a registered clinical trial (ChiCTR2100042346). The study received approval from the Ethics Committee of the Affiliated Brain Hospital of Guangzhou Medical University, and all participants and their parents provided written informed consent. Both AO-MDD patients and HCs underwent a structured clinical interview based on the diagnostic and statistical manual of mental disorders, fifth edition (DSM-5), and the diagnosis of major depressive disorder was made by a trained clinician. All patients met the following inclusion criteria: (i) age between 12 and 18 years old, (ii) first episode of AO-MDD, (iii) 17-item Hamilton Rating Scale for Depression (HAMD-17) , and (iv) had been medication-free for at least 4 weeks prior to inclusion in the trial. First episode and unmedicated AO-MDD patients were recruited from the outpatient department of the Affiliated Brain Hospital of Guangzhou Medical University. Healthy controls matched for age and education level were recruited through advertising. Participants with a history of developmental disorders, tic disorders, or attention deficit hyperactivity disorder were excluded, as were those with any other mental illness (e.g., bipolar disorder, posttraumatic stress disorder, and schizophrenia), neurological disorders or other major physical illnesses, and contraindications to MRI scanning.

The Hamilton Rating Scale for Depression (HAMD-17) was employed to evaluate the severity of depression in AO-MDD [25]. Meanwhile, the Childhood Trauma Questionnaire (CTQ) was utilized for measuring CM, encompassing two major categories of traumatic experiences, namely, childhood neglect and childhood abuse [26, 27]. The average score of all 25 items of the CTQ served as the comprehensive CM score (). The score for childhood abuse was determined by the average score of 15 items that assessed abuse (), while the score for childhood neglect was calculated based on the average score of 10 items that evaluated neglect ().

2.2. Definition of Regions of Interest

In the seed-based RSFC analysis, four subregions of the GP, comprising the bilateral aGP and bilateral pGP, were designated as the ROI based on the brain atlas mapped by Tian et al. The anatomical position of the GP and its subregions are delineated in Figure S1.

2.3. MRI Data Acquisition

All images were acquired using a Siemens Magnetom Prisma 3.0T MRI Scanner, equipped with a 64-channel head coil, at the Magnetic Resonance Center of Affiliated Brain Hospital of Guangzhou Medical University. Participants were provided with instructions to maintain eye closure and wakefulness and minimize cognitive activity during the scanning procedure. The whole-brain functional echo-planar images were obtained transversely with a gradient-recalled echo-planar imaging pulse sequence utilizing the following parameters: repetition  ms, echo  ms, flip angle =56°, slice  mm, number of , and field of  mm. The anatomical image was obtained utilizing a sagittal three-dimensional gradient echo T1-weighted sequence with the following settings: repetition  ms, echo  ms, inversion  ms, flip angle =8°, number of , slice  mm, and .

2.4. MRI Data Preprocessing

The functional images underwent preprocessing using Data Processing and Analysis for Brain Imaging (DPABI_V6.1, http://rfmri.org/dpabi), which was implemented in MATLAB (version R2019b). Signal stabilization was achieved by removing the first 25 functional volumes, and the remaining 425 volumes were corrected for timing differences between slices. Motion correction of the functional images was accomplished using a six motion parameter (rigid body). It is noteworthy that the mean framewise displacement (FD) based on the Jenkinson model (FD-Jenkinson) was computed by averaging the FD from every time point for each participant [28]. No participants with head movement exceeding 2 mm or rotation exceeding 2° were excluded. Following motion correction, the individual structural image (T1-weighted image) was coregistered to the mean functional image. The transformed structural images were then segmented into gray matter, white matter, and cerebrospinal fluid. The rigid-body six model, the white matter signal, and the cerebrospinal fluid signal were subsequently treated as nuisance covariates and regressed out. Finally, the images were temporally band-pass filtered between 0.01 Hz and 0.08 Hz and smoothed using a 4 mm full width at half maximum Gaussian kernel.

2.5. Calculation of RSFC with Globus Pallidus as Regions of Interest

The RSFC analysis was executed utilizing the DPABI software. The RSFC analysis with the four subregions of the GP as ROI comprised of a three-step approach. Firstly, the time series within the selected ROIs were extracted. Subsequently, a voxel-wise correlation analysis was carried out between the subregions of GP and other brain areas to derive RSFC maps. Finally, the obtained RSFC maps underwent a Z-transformation to enhance the normality of the data distribution.

2.6. Statistical Analysis

In this study, the statistical analysis can be segregated into four distinct components as outlined below. (i)Differences in demographic characteristics between groups, such as age, gender, education level, and duration of the disease, as well as clinical characteristics such as HAMD-17 score and CTQ score, were analyzed using -tests and tests with SPSS 23.0 (IBM Corporation, Armonk, NY, USA)(ii)We utilized two-sample -tests to examine differences between patients with AO-MDD and HCs in RSFC analysis with the GP serving as the ROI. To control for potential confounding variables, we included gender, age, education level, and mean FD as covariates. Additionally, given that the majority of participants were female (59%), we conducted a validation analysis for females to verify the robustness of our results (refer to supplementary materials). Gaussian random field (GRF) correction was performed with a significance threshold of voxel- and cluster-. Additional analysis was conducted to investigate the relationship between altered reward circuit function and the severity of anhedonia, as measured by the score of item 7 of the HAMD-17. In specificity analyses, we explored whether analogous results would emerge for the accumbens (comprising the nucleus accumbens shell and nucleus accumbens core)(iii)Correlational analyses were performed separately on AO-MDD patients and HCs to investigate the relationship between the severity of depression (including dimensions such as depression severity, anxiety/somatization, retardation, cognitive disturbance, sleep disruption, and weight) and traumatic experiences (such as the CTQ total score, childhood abuse score, and childhood neglect score)(iv)Multiple linear regression was further employed to investigate whether the abnormal RSFC between dorsomedial prefrontal cortex (DMPFC) and pGP (serving as moderating variable ) moderated the relationship between traumatic experiences (serving as independent variable , encompassing childhood abuse, childhood neglect, and CTQ total score) and AO-MDD severity (serving as dependent variable ) in AO-MDD patients and HCs. Covariates such as age, education, and gender were treated as controlling variables. Traumatic experiences and abnormal RSFC between DMPFC and pGP were considered continuous variables in all models. In the post hoc specificity analysis, anxiety severity (represented by HAMA score) was utilized as an alternative outcome variable, with traumatic experiences serving as the independent variable and abnormal RSFC between DMPFC and pGP as the moderating variable

3. Results

3.1. Demographic and Clinical Characteristics

The demographic and clinical characteristics of all participants are delineated in Table 1. Notably, AO-MDD patients did not exhibit any significant difference from HCs in relation to age, education level, and mean FD. However, there was a significant difference between the two groups with regard to gender, which warranted consideration as a covariate to mitigate potential confounding effects. Furthermore, significant differences were observed between AO-MDD patients and HCs in relation to HAMD-17 total score, HAMA, CTQ total score, childhood abuse score, and childhood neglect score.

3.2. Correlation Analysis between Depression Severity and CM

Correlation analyses were conducted to investigate the relationship between depression severity and CM in AO-MDD patients and HCs. Specifically, the HAMD-17 total score and its dimensions (namely, depression severity, anxiety/somatization, retardation, cognitive disturbance, sleep disruption, and weight) were analyzed alongside the CTQ total score and its subscales (including childhood abuse and childhood neglect). As demonstrated in Figure S2, the HAMD-17 total score exhibited a significant association with childhood abuse, while sleep disruption was significantly associated with the CTQ total score in AO-MDD patients. However, in HCs, we only observed a correlation between the HAMD-17 total score and childhood abuse.

3.3. Abnormal RSFC with Globus Pallidus as Regions of Interest in AO-MDD Patient

Analyses of two-sample -test were conducted to examine variations in ROI-based RSFC between patients with AO-MDD and HCs. As illustrated in Figure 1 and Table S1, there was a statistically significant decrease in RSFC between the left pGP and the DMPFC (BA8) in AO-MDD patients compared to HCs. However, there was no significant difference between AO-MDD patients and HCs in the RSFC analyses involving the bilateral aGP and the right pGP as ROIs. We examined the RSFC of DMPFC-pGP in relation to anhedonia severity, which was assessed using the score of item 7 of the HAMD-17. Our analysis revealed a close-to-significant association with anhedonia severity, with an -value of 0.199 and a value of 0.086. Additionally, we did not observe any significant differences between AO-MDD patients and HCs in the RSFC analyses involving the four subregions of the nucleus accumbens as ROIs.

3.4. Moderation Analysis

Table 2 presents the results of the multiple linear regression analysis conducted on AO-MDD patients to investigate the impact of traumatic experiences on AO-MDD severity. Model 1 highlights that the interaction between RSFC of DMPFC-pGP and CTQ total score had a significant effect on AO-MDD severity. In model 2, the interaction between RSFC of DMPFC-pGP and childhood abuse was also found to be significant. However, the interaction between RSFC of DMPFC-pGP and childhood neglect in model 3 did not reach statistical significance. Additionally, Table S2 of the supplementary material shows that no significant interaction effects were detected in the HCs. Therefore, the results presented below pertain exclusively to the dataset of AO-MDD patients.

To further decompose the interaction effects of RSFC between the DMPFC and the pGP, as well as traumatic experiences (given that model 3 failed to attain significance, only the CTQ total score and childhood abuse were incorporated) on the severity of AO-MDD, patients with AO-MDD were stratified into high, medium, and low reward circuit function groups. Additional multiple linear regressions were executed individually in each subgroup, wherein high, medium, and low reward circuit functions were denoted by high, medium, and low RSFC between DMPFC and pGP, respectively. Specifically, the low reward circuit function group encompassed AO-MDD patients exhibiting RSFC between DMPFC and (below mean -0.5 standard deviation (SD)), the medium reward circuit function group comprised AO-MDD patients displaying RSFC between DMPFC and and < -0.0040 (mean -0.5 SD to mean +0.5 SD), and the high reward circuit function group included AO-MDD patients possessing RSFC between DMPFC and (above mean +0.5 SD). The mean ±0.5 SD threshold was selected to ensure similarity of participant numbers in each subgroup. These multiple linear regressions treated CTQ total score and childhood abuse as independent variables, AO-MDD severity as the dependent variable, and gender, age, and education level as covariates. As evinced in Tables 3 and 4 and Figure 2, CTQ total score and childhood abuse manifested a significant association with AO-MDD severity solely in the low reward circuit function group, and the more severe the traumatic experiences, the more severe the AO-MDD.

To investigate whether the moderating impact of aberrant RSFC between the DMPFC and pGP was exclusive to the severity of AO-MDD, we conducted an additional moderation analysis. In this analysis, we utilized anxiety severity, as determined by HAMA scores, as the dependent variable, traumatic experience as the independent variable, and aberrant RSFC between the DMPFC and pGP as the moderating variable. Subsequently, a post hoc specificity analysis revealed that the moderating impact of aberrant RSFC between the DMPFC and pGP was only observed with respect to the severity of AO-MDD and not in relation to anxiety severity. This indicates that the moderating effect of abnormal RSFC between the DMPFC and pGP is specific to the severity of AO-MDD (refer to Table S3 for details).

3.5. Validation

Furthermore, given that a majority of our study participants were of the female sex (59%), we conducted the aforementioned analyses solely on female subjects to ascertain the robustness of our results. As depicted in Table S3 of the supplementary materials, we observed that our findings remained unaffected by any potential gender bias.

4. Discussion

In the present study, we examined the aberrant RSFC with the GP as the ROI in patients with AO-MDD. Our aim was to explore the relationship between CM, AO-MDD severity, and the abnormal functioning of the reward circuit centered on the GP. The following findings were obtained: (i) A decrease in the RSFC between the pGP and the DMPFC was observed in AO-MDD patients compared to HCs, suggesting that patients with AO-MDD had an aberrant reward circuit function with the GP as its core. (ii) The association between CM and AO-MDD severity was moderated by the abnormal functioning of the reward circuit centered on the GP. In AO-MDD patients with reduced reward circuit function, greater CM was associated with increased severity of AO-MDD. (iii) Similarly, reward circuit abnormalities moderated the association between childhood abuse and AO-MDD severity. However, this moderating effect of reward circuit abnormalities was not observed between childhood neglect and AO-MDD severity. (iv) Post hoc specificity analysis revealed that the moderating effect of reward circuit abnormalities was observed only between CM/childhood abuse and AO-MDD severity, but not between CM and anxiety severity. This indicated that the moderating effect of reward circuit abnormalities was specific to depression severity in AO-MDD patients. Furthermore, validation analysis demonstrated that our findings were not influenced by gender bias, indicating that our results were robust. Our examination of reward circuits in individuals with AO-MDD enhances the understanding of the diminished positive emotions observed in AO-MDD and offers novel insights into the neural mechanisms underlying mood and subjective experience in this population.

Most prior investigations concerning reward anticipation have focused on abnormalities in the reward circuit associated with the nucleus accumbens, while disregarding the significant role of the GP in this circuit. Building on Ottenheimer et al.’s innovative proposal for a reward circuit centered on the GP [16], the present study provides the first confirmation that patients with AO-MDD demonstrate reduced RSFC between the pGP and DMPFC, offering novel evidence of abnormal reward circuit in AO-MDD patients. There is prior research that has shown that within the reward circuitry, information flow from the ventral striatum (where the GP is a component) can project to the frontal cortex, including the prefrontal regions [29]. The pGP is recognized as a critical region for motor and motivational functions, as it receives inputs from various reward-related neural circuits, effectively linking goal-driven behaviors with the anticipation of rewards [30]. Animal models of adverse life events have revealed alterations in reward-seeking behaviors [31] and increased self-administration of various psychoactive substances [32]. In addition, Pechtel and Pizzagalli found that females who had experienced childhood sexual abuse exhibited lower behavioral accuracy when utilizing learned information in reward tasks [33]. Additionally, we observed that reduced reward circuit function moderated the relationship between CM and AO-MDD severity, lending support to the notion that not all AO-MDD patients who experience CM suffer significant negative effects from traumatic experiences and that AO-MDD patients with diminished reward circuit function are more sensitive to CM. These findings are consistent with prior research; for example, a prospective cohort study by Vidal-Ribas et al. found that exposure to stressful life events during childhood predicted reduced neural responses to reward anticipation in basal ganglia regions, particularly the nucleus accumbens and GP [34]. Additionally, diminished responses to reward anticipation in the medial prefrontal and anterior cingulate cortex predicted greater stress reactivity in the future [34]. This offers evidence for a bidirectional association between stress and reward circuit function, whereby early stress may impact reward anticipation, while reduced reward anticipation may increase susceptibility to stress [34]. Furthermore, Ding et al. discovered that recurrent MDD patients exhibit reduced nucleus accumbens RSFC in both the reward network and default mode network [35]. Taken together with prior research, our findings confirm that dysregulation of corticosubcortical connectivity may be a plausible mechanism for abnormal reward processing in AO-MDD [36].

We have observed that abnormalities in the reward circuit moderated the association between childhood abuse and the severity of AO-MDD. Specifically, we found that the impact of childhood abuse on AO-MDD severity was attenuated in patients exhibiting higher reward circuit function. Conversely, we did not identify a moderating effect of reward circuit abnormalities on the relationship between childhood neglect and AO-MDD severity, indicating that childhood abuse and childhood neglect have differential effects on patients with AO-MDD. Prior investigations have reported similar findings, namely, that childhood abuse is strongly associated with affective symptoms, while childhood neglect is highly correlated with cognitive symptoms [37, 38]. For example, Kilian et al. have demonstrated that a history of childhood neglect is a significant predictor of impaired social cognition in adulthood, including attention/agility and verbal learning ability [37]. Furthermore, children who experience childhood abuse, particularly physical and sexual abuse, are more susceptible to developing symptoms of emotional dysregulation in adulthood, such as impulsive aggression, repeated self-harm, and chronic suicidal tendencies [39]. An important and as-yet-unanswered question is why childhood abuse and neglect have differing effects on AO-MDD patients. Existing research suggests that the neurobiological mechanisms underlying the effects of childhood abuse and neglect on AO-MDD patients are distinct [40, 41], with childhood abuse exerting a more significant impact on reward circuit function in patients with AO-MDD [42]. Our discovery that abnormal reward circuits modulate the relationship between childhood abuse and AO-MDD severity provides novel insights into exploring the underlying neural mechanisms in AO-MDD patients who have suffered childhood abuse.

Upon performing a post hoc specificity analysis, it was discovered that the moderating effect of abnormalities within the reward circuit was exclusively evident in the association between CM and the severity of AO-MDD, while no such effect was observed between CM and anxiety severity. This suggests that the moderating effect of reward circuit abnormalities is specific to the severity of AO-MDD, which may be due to the distinct neurobiological mechanisms underlying depressive and anxiety symptoms. Using a mouse model, Fan et al. discovered a strong link between depressive-like states associated with social status loss and abnormal function of the reward circuit [43]. Soares-Cunha et al. found that subtypes of spiny neurons located in the nucleus accumbens, a critical area of the reward circuit, signal both reward and aversion to regulate depression-like symptoms [44]. Additionally, Morel et al. utilized the chronic social defeat stress paradigm to reveal that the midbrain projection to the basolateral amygdala encodes anxiety-like behaviors, but not depression-like behaviors [45]. In conjunction with prior research, we suggest that the decreased function of the reward circuit is related to the underlying neurobiological mechanisms of depression severity, while overactivation of the amygdala-related pathway is associated with anxiety symptoms.

The present study is subject to certain limitations. Firstly, our study’s cross-sectional design does not permit us to draw inferences regarding the causal relationships among variables. Secondly, due to the limited sample size, we were unable to categorize patients with AO-MDD into AO-MDD with and without CM groups based on their exposure to traumatic events. Thirdly, the evaluation of CM is based on retrospective questionnaires, which may be subject to personal bias and memory distortion. However, it is pertinent to note that self-reported accounts of CM are considered valid in the absence of objective evidence, as they represent an individual’s perception of early-life adversity [46]. Fourthly, due to the gender mismatch between the AO-MDD patients and the HCs, gender was accounted for as a covariate. To assess whether gender bias influenced our findings, we performed a validation analysis exclusively on female participants. Our results, as demonstrated in the Supplementary Material, were not impacted by gender bias.

In sum, we have observed a decline in RSFC between the left pGP and the DMPFC among patients with AO-MDD. This suggests that patients with AO-MDD have a weakened reward circuit function centered on the GP. Furthermore, the analysis of moderation and interaction decomposition indicated that AO-MDD patients with lower reward circuit function exhibited higher susceptibility to CM, i.e., hypersensitivity to traumatic experiences. Our findings emphasize the crucial role of reward circuit abnormalities in the onset of AO-MDD and offer new leads to explore the relationship between CM, reward circuit dysfunction, and AO-MDD.

Data Availability

The data that support the findings of this study are available from the corresponding authors upon reasonable request.

Disclosure

The funding source had no role in the study design, analysis, or interpretation of data or in the preparation of the report or decision to publish.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

WC acquired the data, analyzed and interpreted the data, and drafted the manuscript. ZB designed and conceptualized the study, analyzed and interpreted the data, and critically revised the manuscript. CW, XF, MZ, FZ, ZR, YX, HY, ZJ, YX, YY, YF, KW, GH, JC, and GX acquired the data and critically revised the manuscript. YL and YP critically revised the manuscript. All authors read and approved the final manuscript. Weicheng Li and Zhibo Hu contributed equally to this work.

Acknowledgments

The authors declare the following financial interests/personal relationships which may be considered as potential competing interests. Author Yuping Ning has received research grant support from the Science and Technology Plan Project of Guangdong Province (grant number 2019B030316001) and Science and Technology Program of Guangzhou Liwan District (grant number 202201003). Author Yanling Zhou has received research grant support from the Guangdong Basic and Applied Basic Research Foundation (grant number 2022A1515011567) and Guangzhou Science and Technology Planning Project (grant numbers 202103000032 and 202201010714). Author Chengyu Wang has received research grant support from the Science and Technology Plan Project of Guangzhou (2023A03J0842). Author Xiaofeng Lan has received research grant support from the Guangzhou Traditional Chinese Medicine and Integrated Traditional Chinese and Western Medicine Technology Project (grant number 20232A010013). This work was supported by the National Natural Science Foundation of China (grant number 81801343), Guangdong Basic and Applied Basic Research Foundation (grant number 2019A151501136), Science and Technology Plan Project of Guangdong Province (grant number 2019B030316001), Guangdong Medical Science and Technology Research Project (grant number B2020178), Guangdong Basic and Applied Basic Research Foundation Outstanding Youth Project (2021B1515020064), Science and Technology Plan Project of Guangzhou (grant number 202102020557), and Guangzhou Municipal Key Discipline in Medicine (2017-2019).

Supplementary Materials

Figure S1: visualization of the subregions of globus pallidus. Four seeds of the globus pallidus in the bilateral hemisphere. Figure S2: correlation analysis between depression severity and CM in the AO-MDD and HC groups. Table S1: regions showing significant differences in functional connectivity among the AO-MDD and HC groups (). Table S2: multiple linear regression of HAMA total by CTQ total and abuse and neglect with RSFC as a moderator. Table S3: regions showing significant differences in functional connectivity among the AO-MDD and HC groups in female participants (). (Supplementary Materials)