
==== Front
Mol Psychiatry
Mol Psychiatry
Molecular Psychiatry
1359-4184
1476-5578
Nature Publishing Group UK London

38351174
2450
10.1038/s41380-024-02450-7
Article
Divergent suicidal symptomatic activations converge on somato-cognitive action network in depression
http://orcid.org/0000-0003-3517-6420
Li Jiao jiaoli@uestc.edu.cn

12
Wang Dajing 12
Xia Jie 12
Zhang Chao 12
Meng Yao 12
Xu Shuo 12
http://orcid.org/0000-0002-4062-4753
Chen Huafu chenhf@uestc.edu.cn

12
http://orcid.org/0000-0001-7406-7193
Liao Wei weiliao.wl@gmail.com

12
1 https://ror.org/04qr3zq92 grid.54549.39 0000 0004 0369 4060 The Clinical Hospital of Chengdu Brain Science Institute, School of Life Science and Technology, University of Electronic Science and Technology of China, Chengdu, 611731 P.R. China
2 https://ror.org/04qr3zq92 grid.54549.39 0000 0004 0369 4060 MOE Key Lab for Neuroinformation, High-Field Magnetic Resonance Brain Imaging Key Laboratory of Sichuan Province, University of Electronic Science and Technology of China, Chengdu, 611731 P.R. China
14 2 2024
14 2 2024
2024
29 7 19801989
25 9 2023
19 1 2024
23 1 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Individuals with depression have the highest lifetime prevalence of suicide attempts (SA) among mental illnesses. Numerous neuroimaging studies have developed biomarkers from task-related neural activation in depressive patients with SA, but the findings are inconsistent. Empowered by the contemporary interconnected view of depression as a neural system disorder, we sought to identify a specific brain circuit utilizing published heterogeneous neural activations. We systematically reviewed all published cognitive and emotional task-related functional MRI studies that investigated differences in the location of neural activations between depressive patients with and without SA. We subsequently mapped an underlying brain circuit functionally connecting to each experimental activation using a large normative connectome database (n = 1000). The identified SA-related functional network was compared to the network derived from the disease control group. Finally, we decoded this convergent functional connectivity network using microscale transcriptomic and chemo-architectures, and macroscale psychological processes. We enrolled 11 experimental tasks from eight studies, including depressive patients with SA (n = 147) and without SA (n = 196). The heterogeneous SA-related neural activations localized to the somato-cognitive action network (SCAN), exhibiting robustness to little perturbations and specificity for depression. Furthermore, the SA-related functional network was colocalized with brain-wide gene expression involved in inflammatory and immunity-related biological processes and aligned with the distribution of the GABA and noradrenaline neurotransmitter systems. The findings demonstrate that the SA-related functional network of depression is predominantly located at the SCAN, which is an essential implication for understanding depressive patients with SA.

Subject terms

Depression
Diagnostic markers
https://doi.org/10.13039/501100001809 National Natural Science Foundation of China (National Science Foundation of China) 82202250 62036003 Li Jiao Chen Huafu https://doi.org/10.13039/501100002858 China Postdoctoral Science Foundation BX2021057 2022M710615 Li Jiao The National Key R&amp;D Program of China (2022YFC2009906 and 2022YFC2009900).The Fundamental Research Funds for the Central Universities (ZYGX2022YGRH008)issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Depression is a common mental illness, with approximately 280 million people affected around the world [1]. A systematic review indicated that depressive patients have a higher risk of suicidality than individuals without depression [2]. Furthermore, the lifetime mortality rate due to suicide in patients with depression reaches 15% [3]. Suicide attempts (SA) are an important factor contributing to mortality [4]. Thus, identifying the underlying brain circuits in depressive patients with SA could help us understand the relevant neurological basis and ultimately prevent suicide.

Considerable efforts have been dedicated to elucidating neural mechanisms in depressive patients with SA during various cognitive and emotional tasks [5, 6], including the Balloon Analog Risk Task [7], the facial emotion processing task [8–10], the Go/NoGo task [11], the Iowa Gambling task [12], the discounting task [13], and the Cyberball Game [14]. However, loci of neural activation are plagued with low reproducibility. For instance, Pan et al. found that depressive patients with SA exhibited significantly increased activation in the right anterior cingulate gyrus, left dorsolateral prefrontal cortex, primary sensory cortex, and right middle temporal gyrus during a facial emotional task compared with depressive patients without SA [10]. However, Chase et al. did not find different activations between depressive patients with SA and without SA during the same task [15].

The conventional activation likelihood estimation (ALE) meta-analysis technique [16] has been used to identify commonality amidst the heterogeneities to solve the “reproducibility crisis” in neuroimaging [17, 18]. However, reproducibility in ALE meta-analysis is mainly focused on whether the same brain regions are activated, ignoring the connectivity among activated regions. Recently, there has been a shift in the localization of psychiatric diseases from isolated brain regions to shared brain networks [19, 20]. Cash and colleagues [20] identified a robust and clinically significant distributed brain network of depression that was irreconcilable with the ALE method using a task-related database. Considering depression as a neural system disorder by the contemporary interconnected view [21] and the functional integration between different brain regions during psychological processes [22], we hypothesized that heterogeneous abnormal neuroimaging activations in depression with SA could map to a common brain network.

To uncover the common brain network underpinning SA in depressive patients by the abnormal activation coordinates during task performance, we chose a recently proposed method termed “activation network mapping” (ANM) [23], which was developed from the “lesion network mapping” technique [24]. The ANM method replaces brain lesions with reported activation coordinates as input to obtain the brain network of specific symptoms based on a large cohort of resting-state normative functional connectivity. The ANM technique has been used to derive brain networks for schizophrenia with the symptoms of formal thought disorder and auditory verbal hallucinations [25]. Similar to ANM, previous studies have used the coordinate-based network mapping technique to identify networks associated with depression and late-life depression using brain structure differences [19] and a shared brain network associated with Parkinson’s disease with dementia from neuroimaging findings [26]. Overall, the ANM method integrates the heterogeneous neuroimaging data to identify a common network in depression with SA from task-based functional MRI (fMRI) studies.

In this study, we aimed to identify the corresponding brain circuit in depressive patients with SA using the ANM method. We first tested whether the heterogeneous abnormal brain activations in depressive patients with SA across multiple experimental tasks would localize to a common network. Then, we tested whether this identified network was robust to perturbations of analyzed parameters and experimental tasks and was specific in depression compared with the disease control group. Finally, we decoded the identified brain network by microscale transcriptomic and chemo-architectures and macroscale psychological processes.

Methods

Search strategy

This systematic meta-analysis was conducted according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) standard [27].

The search for neuroimaging studies investigating functional alterations in depressive patients with SA compared to those without SA was performed up to 1 October 2022 on PubMed and Web of Science. Keywords restricted to human studies were: (suicide OR “suicide*”) AND (fMRI OR “Functional MRI” OR “Functional MRI” OR “Functional Magnetic Resonance Imaging”) AND (depression OR MDD OR “major depression” OR “major depressive disorder” OR “depressive disorder” OR depressed). Other sources included reference lists of previous reviews and meta-analysis articles that summarized the studies about task-related alterations in depressive patients with SA.

Data extraction

The inclusion criteria were as follows: (i) specific tasks were performed during fMRI scanning; (ii) depressive patients without SA were included; (iii) whole-brain activations were compared between depressive patients with and without SA; (iv) the Talairach or Montreal Neurological Institute (MNI) coordinates of significant brain regions were reported; (v) studies did not include other mental illness patients other than depressive patients. The following exclusion criteria were used: (i) specific regions of interest were analyzed, or the brain was only partially covered; (ii) no significant differences between depressive patients with SA and without SA were reported; (iii) the coordinates were not reported; (iv) differences in functional connectivity rather than activation were examined.

Information was extracted from each identified study, including author, publication year, experimental tasks, sample size, demographics, clinical scale score, and coordinates. The Talairach coordinates were converted into the MNI coordinates.

Activation network mapping

The ANM method [23] was used to test whether heterogeneous abnormal activation regions in depression with SA localize to a common brain network. This method is similar to coordinate-based network mapping [28]. Resting-state functional connectivity (RSFC) analysis section: First, a 4-mm-radius sphere at each reported coordinate was created. Spheres from the same experimental task were merged to produce a combined seed. Eleven combined seeds were generated for subsequent analyses. Second, functional connectivity was calculated between each combined seed and voxels at the whole-brain level for each subject based on a publicly available resting-state connectome database [29, 30], i.e., Brain Genomics Superstruct Project (https://dataverse.harvard.edu/dataverse/GSP). Third, each subject-level correlation r-map (above 1000) was transformed to a Fisher z-map. Then, the experiment-level t-map was generated by comparing all the 1000 subject-level Fisher z-maps from the same experimental task against zero using the one-sample t test. Overlapping section: First, each experiment-level t-map was thresholded at t > 3 (corresponding to voxel-wise false discovery rate [FDR] correction with P < 0.01). Second, all binarized maps were overlapped and thresholded at 60% [23] to create a group-level network overlap map, named SA-related ANM in depression. Several control analyses were performed to test the robustness of the identified SA-related ANM in depression (Supplementary Methods).

Specificity analysis of the SA-related ANM in depression

To assess the specificity of the SA-related ANM in depression, we explored SA-related ANM in patients with bipolar disorder (BD) and schizophrenia (SCZ) (hereafter termed “disease control group”) using the same methods and parameters as in depression patients. Four experimental tasks were identified from three articles, which explored the brain activity differences between BD or SCZ patients with SA and without SA. Considering the limited number of studies in BD and SCZ patients with SA, the two disorders were defined as one disease control group. A priori-defined region of interest (ROI) in the insula, obtained from a meta-analysis that integrated the task-related brain activation abnormalities in depression with SA [31], was used to test the specificity of the SA-related ANM in depression. All experiment-level t-statistics were converted to z-statistics using the t2z function in SPM12. Group mean z-maps were separately calculated for patients with depression and the disease control group. Next, the voxels’ z-values within the insula were extracted from mean z-map for each group. Finally, differences of z-value distribution were evaluated between the depression and disease control groups using Cohen’s d test. Meanwhile, the precentral and postcentral gyri were used as the post-hoc ROI to test the specificity of SA-related ANM in depression. The priori-defined and post-hoc ROIs were generated from the AAL3 brain atlas [32].

Transcriptomic analysis

The Allen Human Brain Atlas dataset [33] (https://human.brain-map.org/) allows us to investigate the relationship between SA-related ANM in depression and transcriptome. Gene expression data of brain tissue samples were collected from six postmortem human brains. The microarray expression data were prerocessed by the abagen toolbox [34], which provides reproducible workflows to process gene expression data following an established practical guide [35].

After preprocessing gene expression data using the abagen toolbox with the Human Brainnetome atlas with 246 regions [36], we obtained a regional expression matrix for each donor with 246 rows corresponding to brain regions, and 15,633 columns corresponding to the retained genes. To get stably expressed genes, we selected the genes with a mean similarity across brain regions above 0.1 between donors, resulting in 12,824 genes. We then averaged the retained genes across all donors. Because only two donors have gene expression in the right hemisphere, we extracted the left hemisphere gene expression data for the following analyses. Consequently, we generated a 123 × 12,824 matrix describing the transcript levels of the left hemisphere.

Partial least squares (PLS) correlation analysis [37] was used to determine the relationship between SA-related ANM in depression and the gene expression patterns of all 12,824 genes. The SA-related ANM in depression was firstly downsampled to 246 regions, i.e., the averaged value of all voxels within the region as the region’s value. Gene expression data were used as predictor variables, and SA-related ANM in depression were used as the response variables in PLS correlation. The first component of PLS (PLS1) was strongly correlated with SA-related ANM in depression across 123 regions. To correct the spatial autocorrelation, we test the statistical significance of variance explained by randomizing the response variables 5,000 times using the Moran spectral randomization [38, 39]. A bootstrapping method was used to estimate each gene’s weighting coefficient in the PLS analysis. The ratio of the weight of each gene to its bootstrap standard error (5000 times) was utilized to calculate z scores of each gene weight on PLS1 [40]. Genes were ranked according to their contribution to PLS1.

The gene enrichment analysis was conducted by Metascape (https://metascape.org/), which can perform gene meta-analysis with over 40 independent knowledge bases [41]. Genes with PLS1 weights of z > 3.69 or z < −3.69 (all PFDR < 0.001) were input into the Metascape website to perform Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) biological processes analyses with the background of 16,831 brain-expressed genes [42]. The background genes were obtained by excluding probes if they did not exceed the background in at least 20% of samples after probe to gene re-annotation. The obtained biological processes were corrected by FDR correction (P < 0.05).

Neurotransmitter receptors and transporters analysis

We assessed the relationships between SA-related ANM in depression and neurotransmitter systems across all brain regions. To this end, we obtained whole-brain volumetric PET receptor images from a cohort of over 1200 healthy individuals to construct an overview of 19 distinct neurotransmitter receptors and transporters across nine different neurotransmitter systems [43] (https://github.com/netneurolab/hansen_receptors/tree/main/data/PET_nifti_images). Based on the previous study [44], five neurotransmitters associated with depression were analyzed in this study, including serotonin (5-HT1a, 5-HT1b, 5-HT2a, 5-HT4, 5-HT6, and 5-HTT), noradrenaline (NAT), dopamine (D1, D2, and DAT), glutamate (mGluR5 and NMDAR), and gamma-aminobutyric acid (GABAA). Similar to the procedures in Hansen et al., receptors and transporters with more than one mean image of the same tracer (i.e., 5-HT1b, D2, and mGluR5) were averaged together in a manner that weights each image by the number of participants in the cohort [43]. Different neurotransmitter receptors and transporters from the same category were averaged to generate a corresponding neurotransmitter system map. The five neurotransmitter system maps were assigned to 246 regions, and the averaged value of all voxels within the region was defined as the region’s neurotransmitter value. After each neurotransmitter map z-normalized across 246 regions, we used a multivariate linear regression model to explore the contributions of neurotransmitter systems to the SA-related ANM in depression. The relaimpo (relative importance of regressor in linear models, version 2.2-5) package in R was used. Relative importance metrics can be used to address linear regression with multiple collinear regression [45]. The model was defined as follows:Ov=β0+β1×Dopamine+β2×GABA+β3×Glutamate+β4×Noradrenaline+β5×Serotonin+ϵ,

where Ov (overlap value) is the values of 246 ROIs from the SA-related ANM in depression. After determining the percentage of explained variance from these predictors, we assessed the relative contribution of each neurotransmitter system using bootstrapping. The multiple comparisons were corrected by FDR with P < 0.05.

NeuroSynth meta-analysis

The SA-related ANM in depression was decoded by the NeuroSynth database (https://www.neurosynth.org) to identify the relevant cognitive functions [46, 47]. A public python code (https://github.com/NeuroanatomyAndConnectivity/gradient_analysis/blob/master/05_metaanalysis_neurosynth.ipynb) was used to perform this analysis [48]. Relevant cognitive functions were identified using FDR correction (P < 0.05).

Results

The SA-related ANM in depression

We identified eight studies with 11 experimental tasks that reported abnormal activations in depressive patients with SA (total n = 147) compared to depressive patients without SA (total n = 196) (Fig. 1; Table 1). Brain regions with abnormal activation from these studies were highly heterogeneous, but the majority (7/11, 64%) of combined seeds, generated by the abnormal activation coordinates (Supplementary Table 1) from the same experimental task, connected to the somato-cognitive action network (SCAN). The amygdala, left hippocampus, and right orbitofrontal and left temporal cortices were also included in this network map (Fig. 2). The identified SA-related ANM in depression was robust to parameter selection of threshold values in the Overlapping section (Dice coefficient = 0.62; Supplementary Fig. 1A) and seed sizes in the RSFC analysis section (6 mm: Dice coefficient = 0.42; 8 mm: Dice coefficient = 0.47; 10 mm: Dice coefficient = 0.49; Supplementary Fig. 1B). Using the leave-one-experiment-out approach, we verified that the identified SA-related ANM in depression is stable between all 11 experimental tasks and the maps with one experimental task removed, showing moderate to high similarity (Dice coefficient = 0.49–0.80; Supplementary Fig. 1C).Fig. 1 The PRISMA selection diagram of task-related neuroimaging studies in depression with SA.

Table 1 Experimental tasks and clinical information of subjects in included studies.

First author (year)	Experimental task	Number of subjects	Mean age (years)	BDI	HDRS	
		MDD+SA N=147	MDD N=196	HC N=164	MDD+SA	MDD	HC	MDD+SA	MDD	HC	MDD+SA	MDD	HC	
Jollant et al. (2008) [9]	Happy face;

Angry face

	13	14	16	40.3	43.9	32.4	7.2	8.1	2.6	2.4	2.9	1.1	
Pan et al. (2011) [11]	Go/NoGo	15	15	14	16.2	15.9	15.2	15.5	4.4	2.1	N.A.	N.A.	N.A.	
Pan et al. (2013a) [10]	Happy face;

Angry face

	14	15	15	16.2	15.9	15.3	15.1	4.4	1.9	N.A.	N.A.	N.A.	
Pan et al. (2013b) [12]	Iowa Gambling Task (IGT)	15	14	13	16.2	15.8	15.2	15.5	4.6	2.2	N.A.	N.A.	N.A.	
Vanyukov et al. (2016) [13]	Discounting task of value difference;

Discounting task of tracking value of choices with longer versus shorter delays

	13	13	22	70.4	73.4	72.1	N.A.	N.A.	N.A.	14.5	13.2	2.5	
Olie et al. (2017) [14]	Cyberball Game	36	41	28	39.5	37.6	38.9	3.5	3	0	3.5	2	2	
Ai et al. (2018) [8]	Faces task	18	54	26	37.7	37.6	39.0	N.A.	N.A.	N.A.	N.A.	N.A.	N.A.	
Ji et al. (2021) [7]	Balloon Analog Risk Task (BART)	23	30	30	21.4	23.4	20.6	35.2	33.2	3.7	N.A.	N.A.	N.A.	
BDI Beck Depression Inventory score, HDRS Hamilton Depression Rating Scale (17-item) score, MDD+SA major depressive disorder patients with suicide attempts, MDD major depressive disorder patients without suicidal attempts, HC healthy controls, NA not available.

Fig. 2 Network localization of heterogeneous activations in depression with SA.

A Four-millimeter spheres centered on abnormal coordinates reported in each experimental task were created and combined to create a combined seed. B Brain regions are significantly connected to each combined seed across a normative connectome (n = 1000). All 1000 maps for each experimental task were compared against zero using a voxel-wise one-sample t test and with a threshold of t > 3. C Binarized (t = 3) t maps from all experimental tasks were overlapped to identify the regions connected to most of the combined seeds.

Specificity of the SA-related ANM in depression

To determine the specificity of SA-related ANM in depression, we identified three studies with four experimental tasks in the disease control group (Supplementary Fig. 2 and Fig. 3; Supplementary Table 2), which reported the abnormal activation coordinates in BD and SCZ patients with SA (total n = 52) compared to patients without SA (total n = 81) (Fig. 3A; Supplementary Table 3). We found that the SA-related ANM in depression exhibited stronger connectivity in the prior-defined insula and post-hoc ROIs of precentral and postcentral gyri than the disease control group (Fig. 3B).Fig. 3 The specificity of SA-related ANM in depression.

A Left: Combined abnormal activation coordinates across all experimental tasks of the disease control group. Coordinates from the same experimental task are shown with the same color. Right: SA-related ANM in the disease control group. Brain regions connected to over 60% of the combined seeds are shown in two specific brain slices. B Differences were calculated with Cohen’s d test in the priori insula and the post hoc precentral and postcentral gyri. At the upper right corner of the two density difference maps, the AAL3 mask of the corresponding brain regions is shown.

Molecular mechanisms of the SA-related ANM in depression

We next performed transcriptomic analysis to determine which genes were associated with the SA-related ANM in depression. We found that the PLS1 explained 17.63% variance of the SA-related ANM in depression (Fig. 4A), which was significantly more than expected by chance (Pmoran = 0.03). The spatial distribution of PLS1 was positively correlated with the SA-related ANM in depression (Fig. 4B; r(121) = 0.42, Pmoran = 0.001). This positive correlation means that genes positively (negatively) weighted on PLS1 are overexpressed in the brain regions with higher (lower) overlap [42, 49]. After correcting the normalized weights of PLS1, we found that 1458 positive weighted genes (PLS1+ gene set) and 1444 negative weighted genes (PLS1− gene set) (Fig. 4C; PFDR < 0.001) remained. The PLS1+ gene set is mainly enriched in the inflammatory and immune biological processes, such as “leukocyte activation”, “immune response-activating signaling pathway”, and “primary immunodeficiency” (Fig. 4D). The PLS1− gene set is mainly enriched in the neuro projection- and synapse-related biological processes, such as “modulation of chemical synaptic transmission”, “regulation of synapse organization”, “motor neuron axon guidance”, and “dopaminergic synapse” (Fig. 4E).Fig. 4 Gene PLS analysis and enrichment analyses.

A Spatial map of PLS1 score in the left hemisphere. B The spatial correlation between the PLS1 score and the SA-related ANM in depression. C The ranked gene weights of PLS1 (PFDR < 0.001). D Left: The significantly enriched biological processes (PFDR < 0.05) of the PLS1+ gene set (q-value: the value from the Benjamini–Hochberg procedure for multiple comparison correction; color of the circle: different biological processes corresponding to the nodes of the enrichment network in the right). Right: Metascape enrichment network representing the similarities within each cluster and between clusters (circle nodes represent enriched biological processes; their size depends on the proportion of genes related to the enriched biological process among all input genes; nodes with the same color belong to the same cluster; nodes within the same cluster are typically similar to each other; biological processes with similarity >0.3 are connected through edges). E Left: The significantly enriched biological processes (PFDR < 0.05) of the PLS1− gene set. Right: Metascape enrichment network representing the similarities within each cluster and between clusters.

To further decode the SA-related ANM in depression by chemo-architectures, we investigated the relationship between neurotransmitter systems and the SA-related ANM in depression using a multiple linear regression model (Fig. 5A). The model explained 18% of the variance in the SA-related ANM in depression (F(5240) = 11.5, P < 0.001, adjusted R2 = 0.18; Fig. 5B, C). The noradrenaline and GABA systems significantly contributed to the SA-related ANM in depression (noradrenaline weight = 0.33, PFDR < 0.001; GABA weight = −0.24, PFDR = 0.03). Notably, the noradrenaline system had the highest relative contribution (relative contribution = 63.5%; Fig. 5D).Fig. 5 The relationships between different neurotransmitter systems and the SA-related ANM in depression.

A A multiple linear regression model was used to determine the relationships between different neurotransmitter systems and the SA-related ANM in depression. B The predicted SA-related ANM in depression. C A scatterplot was generated to display the similarity between the fitted network and the observed SA-related ANM in depression. D The relative contributions of each neurotransmitter system were assessed using the relaimpo R package. Error bar, 95% bootstrap confidence intervals. An asterisk represents that P values survived after FDR correction with P < 0.05.

Cognitive functions related to SA in depression

We binarized the SA-related ANM in depression with a 60% threshold to generate a mask for NeuroSynth meta-analysis. We found fifteen significantly related cognitive functions (Fig. 6; PFDR < 0.05). The three most relevant cognitive functions were “execution imagine and somatosensory” (z = 21.1), “observation and goal execution” (z = 17.3), and “performing function” (z = 8.9).Fig. 6 The NeuroSynth meta-analysis of the SA-related ANM in depression.

The terms shown here are significantly associated with the SA-related ANM in depression.

Discussion

Using the recently proposed ANM technique, we found that heterogeneous neural activations were functionally connected to the SCAN in depression with SA. The identified SA-related ANM in depression was robust to the parameter perturbations and leave-one-experiment-out analyses, and was specific for depression. In the transcriptomic and chemo-architectures analyses, we revealed that genes significantly related to the SA-related ANM in depression were mainly enriched in inflammatory and immunity-related biological processes. The spatial distributions of noradrenaline and GABA systems recapitulated the robust SA-related ANM in depression. These findings provide new insight into unifying heterogeneous neuroimaging studies, further advancing our integrative understanding of psychopathology in depression with SA.

The heterogeneous brain regions with abnormal activation in depressive patients with SA across a series of experimental tasks were localized to a common brain network using the ANM approach. This finding might help us to explain the low reproducibility of functional neuroimaging studies. We found that the ANM result exhibited a larger experimental contribution (73%) than the ALE result (27%) to the consistent finding in the left insula (Supplementary Materials; Supplementary Fig. 4). This finding is consistent with numerous studies on network localization of neuropsychiatric symptoms [50], neurodegenerative diseases [28], and emotion processing [23]. The identified SA-related ANM shows robustness against parameter perturbations and specificity for other spectrum disorders (depression patients with SA vs. SCZ and BD patients with SA), suggesting that ANM may be one of the methods to resolve the low reproducibility across functional neuroimaging studies. Specifically, our findings show that regions located outside the abnormal activation loci, but functionally connected to these areas, can also be involved in the depressive process. This finding is consistent with a growing literature suggesting that symptoms localize to connected brain circuits, not to single brain abnormalities [19, 26, 51]. Furthermore, brain lesions caused depression and depression-relieving neuromodulation sites functionally link to a brain network that is comparable across 14 different datasets, suggesting therapeutic utility [52].

Heterogeneous activation findings in depressive patients with SA mainly converged on the SCAN which includes the motor area (M1), supplementary motor area (SMA), thalamus, posterior putamen, and postural cerebellum, and functionally connected to the dorsal anterior cingulate cortex (dACC), parietal cortex, and insula [53]. The SCAN plays an important role in the integration of mind and body, which has been argued to be crucial in psychology. The SMA, M1, insula, and parietal cortex in the SCAN are most often involved in response inhibition, which is the ability to suppress unneeded and unsafe actions or thoughts [54]. Abnormalities of these brain regions can lead to deficient inhibitory control (impulsivity), which may be an important factor in patients with SA [55]. Structural and functional abnormalities in the M1, insula, and parietal cortex were found in depressive patients with SA [10, 56, 57]. In addition, reminiscing on painful memories might harm response inhibition [58]. Overwhelming sad and painful memories can encourage rumination and mind-wandering to compete for limited cognitive resources [58], influencing the ability to stop or suppress suicidal behavior. Ji et al. found that individuals with depression were affected by a solid motivation to evade psychological pain [7]. This motivation can potentially result in suicide behavior by occupying brain resources, thereby disrupting the normal response inhibition function. In conclusion, dysfunction of the SCAN may play an important role in modulating emotion and cognition in depressive patients with SA.

The strongest peak of the SA-related ANM in depression was observed in the thalamus that integrates sensory input from the external environment [59]. The thalamus is crucial for processing information after stress exposure and may act as an interface for the stress response in the body and mind [60]. Suicidal people were shown to have an enlarged thalamus [61] and to suffer from stress dysregulation of the hypothalamus–pituitary–adrenal (HPA) axis [62]. In addition, the SA-related ANM in depression included the somatosensory cortex, amygdala, left hippocampus, and right orbitofrontal and left temporal cortices, which also have great significance in depression with SA. The amygdala, left hippocampus, and right orbitofrontal cortex are essential areas of the limbic system, which is generally known to be involved in emotion regulation and is implicated in depression and suicidal behavior [63, 64]. The somatosensory cortex is involved in body awareness and complex experience of pain, and it is less activated in response to pain stimuli [65]. Pain symptoms in patients with depression occupied a high rate (50%) based on a large-scale community study [66], which is significantly associated with lifetime suicide attempts [67]. Furthermore, a mega-analysis revealed that the organization of cortical networks involved in processing sensory information might be a more stable neuroimaging marker for depression than previously assumed alterations in higher-order neural networks like default mode and frontoparietal control networks [68].

Our transcriptional analysis revealed that genes associated with the SA-related ANM in depression were enriched in several biological processes, such as inflammatory and immune processes. A growing body of evidence indicates that inflammation plays an essential role in individuals with depression and SA [69]. The pro-inflammatory cytokines IL-1, IL-6, and TNF-α were found to be associated with behavioral changes in depressive patients [70]. Compared to the healthy control group, SA individuals displayed significantly different leukocyte levels, which are part of our immune system [71]. Our neurotransmitter system analysis revealed that the GABA and noradrenaline neurotransmitter systems were significantly correlated with the SA-related ANM in depression. GABA serves as the major inhibitory neurotransmitter in the brain, which improves mood, increases relaxation, and alleviates pain. Reduced GABA levels have been observed in postmortem tissues of depressive patients with SA [72]. Noradrenaline is associated with executive functioning and governs cognition, motivation, and intelligence. These three factors are crucial to social relationships, and their dysfunction affects the quality of life of depressive patients [73]. Noradrenaline may underly the pathophysiology of depression, and may be a biomarker for the treatment of depressive patients [74]. Furthermore, the relationships between GABA and noradrenaline neurotransmitters with SA-related ANM in depression were opposite, which may be related to the distinct roles in HPA axis activation. The GABA and noradrenaline can separately inhibit and excite the hypothalamic activation of the HPA axis [62]. The HPA axis is abnormally hyperactivated in depressive patients with SA [62]. These transcriptomic and chemo-architecture analyses possibly demonstrated the molecular mechanism of the SA-related ANM in depression. In addition to micro- and meso-associations with SA-related ANM in depression, we also found behavioral relevance of SA-related ANM, such as the execution function, memory, and attention. These cognitive functions were found to be impaired in depressive patients [75, 76]. In particular, dysfunction of execution [77] and attention [78] was found to be more severe in depression patients with SA.

Several limitations need to be considered in our study. First, given that we used brain activation abnormalities as inputs rather than lesions, our results provide only correlational but not causational information [79]. The causation can be tested with prospective studies and brain stimulation studies [80]. Second, because of the limited brain activation investigation in depressive patients with SA, the number of studies included is relatively low, which may cause negative results. In addition, this study included various cognitive and emotional experimental tasks, which may increase the heterogeneity and impact the reproducibility and generalizability of the findings. However, we conducted the leave-one-experiment-out analysis to validate the stability of SA-related ANM in depression. Similar issues existed in the disease control group (i.e., schizophrenia and bipolar disorders). Thus, the specificity of SA-related ANM in depression should be described cautiously. Future studies can enroll more participants in the specific task to investigate the SA-related ANM in depression [23]. Finally, we explored the relationships between the SA-related ANM in depression and gene expression levels and neurotransmitter distributions of healthy controls. Future studies should collect the patients’ gene expression and neurotransmitter data in the brain, which can reveal the correlated directionality between micro- and meso-scale distributions and the SA-related ANM in depression.

In summary, we integrated the heterogeneous abnormal neural activations into the SCAN in depressive patients with SA using the ANM framework. The findings support the notion that body–mind integration is crucial to depressive patients with SA. The robust SA-related ANM can be decoded by transcriptomic and molecular profiles, possibly suggesting several important multiscale interactions in depression patients with SA.

Supplementary information

Supplemental information

Supplementary information

The online version contains supplementary material available at 10.1038/s41380-024-02450-7.

Acknowledgements

This work was supported by the National Key R&D Program of China (2022YFC2009906 and 2022YFC2009900), the National Natural Science Foundation of China (82202250 and 62036003), the China Postdoctoral Science Foundation (BX2021057 and 2022M710615), the Fundamental Research Funds for the Central Universities (ZYGX2022YGRH008), and the Medical-Engineering Cooperation Funds from University of Electronic Science and Technology of China (ZYGX2021YGLH201).

Author contributions

JL, WL and HC concepted and designed this study. JL and DW analyzed data, interpreted data and created the figures and table. JX, CZ, YM and SX help to perform the analyses and interprete the results. JL, DW and WL wrote the manuscript. All authors reviewed and commented on the manuscript.

Data availability

Human gene expression data that support the findings of this study are available in the Allen Brain Atlas (https://human.brain-map.org/static/download). A publicly available resting-state connectome database was obtained from Brain Genomics Superstruct Project (https://dataverse.harvard.edu/dataverse/GSP). Neurotransmitter receptors and transporters were obtained from https://github.com/netneurolab/hansen_receptors/tree/main/data/PET_nifti_images.

Code availability

Transcriptomic data was preprocessed by abagen toolbox (https://github.com/rmarkello/abagen). Human gene The gene enrichment analysis was conducted by Metascape (https://metascape.org/). Multiple linear-regression model was performed using the relaimpo in R package (http://prof.beuth-hochschule.de/groemping/software/relaimpo/).

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Jiao Li, Dajing Wang.
==== Refs
References

1. Evans-Lacko S Aguilar-Gaxiola S Al-Hamzawi A Alonso J Benjet C Bruffaerts R Socio-economic variations in the mental health treatment gap for people with anxiety, mood, and substance use disorders: results from the WHO World Mental Health (WMH) surveys Psychol Med 2018 48 1560 71 10.1017/S0033291717003336 29173244
Evans-Lacko S, Aguilar-Gaxiola S, Al-Hamzawi A, Alonso J, Benjet C, Bruffaerts R, et al. Socio-economic variations in the mental health treatment gap for people with anxiety, mood, and substance use disorders: results from the WHO World Mental Health (WMH) surveys. Psychol Med. 2018;48:1560–71.29173244
2. Cai H Xie XM Zhang Q Cui X Lin JX Sim K Prevalence of suicidality in major depressive disorder: a systematic review and meta-analysis of comparative studies Front Psychiatry 2021 12 690130 10.3389/fpsyt.2021.690130 34603096
Cai H, Xie XM, Zhang Q, Cui X, Lin JX, Sim K, et al. Prevalence of suicidality in major depressive disorder: a systematic review and meta-analysis of comparative studies. Front Psychiatry. 2021;12:690130.34603096
3. Pantazatos SP Huang YY Rosoklija GB Dwork AJ Arango V Mann JJ Whole-transcriptome brain expression and exon-usage profiling in major depression and suicide: evidence for altered glial, endothelial and ATPase activity Mol Psychiatry 2017 22 760 73 10.1038/mp.2016.130 27528462
Pantazatos SP, Huang YY, Rosoklija GB, Dwork AJ, Arango V, Mann JJ. Whole-transcriptome brain expression and exon-usage profiling in major depression and suicide: evidence for altered glial, endothelial and ATPase activity. Mol Psychiatry. 2017;22:760–73.27528462
4. Marx W Penninx B Solmi M Furukawa TA Firth J Carvalho AF Major depressive disorder Nat Rev Dis Prim 2023 9 44 10.1038/s41572-023-00454-1 37620370
Marx W, Penninx B, Solmi M, Furukawa TA, Firth J, Carvalho AF, et al. Major depressive disorder. Nat Rev Dis Prim. 2023;9:44.37620370
5. Auerbach RP Pagliaccio D Allison GO Alqueza KL Alonso MF Neural correlates associated with suicide and nonsuicidal self-injury in youth Biol Psychiat 2021 89 119 33 10.1016/j.biopsych.2020.06.002 32782140
Auerbach RP, Pagliaccio D, Allison GO, Alqueza KL, Alonso MF. Neural correlates associated with suicide and nonsuicidal self-injury in youth. Biol Psychiat. 2021;89:119–33.32782140
6. Schmaal L van Harmelen AL Chatzi V Lippard ETC Toenders YJ Averill LA Imaging suicidal thoughts and behaviors: a comprehensive review of 2 decades of neuroimaging studies Mol Psychiatr 2020 25 408 27 10.1038/s41380-019-0587-x
Schmaal L, van Harmelen AL, Chatzi V, Lippard ETC, Toenders YJ, Averill LA, et al. Imaging suicidal thoughts and behaviors: a comprehensive review of 2 decades of neuroimaging studies. Mol Psychiatr. 2020;25:408–27.
7. Ji X Zhao J Li H Pizzagalli DA Law S Lin P From motivation, decision-making to action: An fMRI study on suicidal behavior in patients with major depressive disorder J Psychiatr Res 2021 139 14 24 10.1016/j.jpsychires.2021.05.007 34004553
Ji X, Zhao J, Li H, Pizzagalli DA, Law S, Lin P, et al. From motivation, decision-making to action: An fMRI study on suicidal behavior in patients with major depressive disorder. J Psychiatr Res. 2021;139:14–24.34004553
8. Ai H van Tol MJ Marsman JBC Veltman DJ Ruhe HG van der Wee NJA Differential relations of suicidality in depression to brain activation during emotional and executive processing J Psychiatr Res 2018 105 78 85 10.1016/j.jpsychires.2018.08.018 30212727
Ai H, van Tol MJ, Marsman JBC, Veltman DJ, Ruhe HG, van der Wee NJA, et al. Differential relations of suicidality in depression to brain activation during emotional and executive processing. J Psychiatr Res. 2018;105:78–85.30212727
9. Jollant F Lawrence NS Giampietro V Brammer MJ Fullana MA Drapier D Orbitofrontal cortex response to angry faces in men with histories of suicide attempts Am J Psychiat 2008 165 740 8 10.1176/appi.ajp.2008.07081239 18346998
Jollant F, Lawrence NS, Giampietro V, Brammer MJ, Fullana MA, Drapier D, et al. Orbitofrontal cortex response to angry faces in men with histories of suicide attempts. Am J Psychiat. 2008;165:740–8.18346998
10. Pan LA Hassel S Segreti AM Nau SA Brent DA Phillips ML Differential patterns of activity and functional connectivity in emotion processing neural circuitry to angry and happy faces in adolescents with and without suicide attempt Psychol Med 2013 43 2129 42 10.1017/S0033291712002966 23298821
Pan LA, Hassel S, Segreti AM, Nau SA, Brent DA, Phillips ML. Differential patterns of activity and functional connectivity in emotion processing neural circuitry to angry and happy faces in adolescents with and without suicide attempt. Psychol Med. 2013a;43:2129–42.23298821
11. Pan LA Batezati-Alves SC Almeida JRC Segreti A Akkal D Hassel S Dissociable patterns of neural activity during response inhibition in depressed adolescents with and without suicidal behavior J Am Acad Child Psy 2011 50 602 11 10.1016/j.jaac.2011.03.018
Pan LA, Batezati-Alves SC, Almeida JRC, Segreti A, Akkal D, Hassel S, et al. Dissociable patterns of neural activity during response inhibition in depressed adolescents with and without suicidal behavior. J Am Acad Child Psy. 2011;50:602–11.
12. Pan L Segreti A Almeida J Jollant F Lawrence N Brent D Preserved hippocampal function during learning in the context of risk in adolescent suicide attempt Psychiat Res-Neuroim 2013 211 112 8 10.1016/j.pscychresns.2012.07.008
Pan L, Segreti A, Almeida J, Jollant F, Lawrence N, Brent D, et al. Preserved hippocampal function during learning in the context of risk in adolescent suicide attempt. Psychiat Res-Neuroim. 2013b;211:112–8.
13. Vanyukov PM Szanto K Hallquist MN Siegle GJ Reynolds CF Forman SD Paralimbic and lateral prefrontal encoding of reward value during intertemporal choice in attempted suicide Psychol Med 2016 46 381 91 10.1017/S0033291715001890 26446615
Vanyukov PM, Szanto K, Hallquist MN, Siegle GJ, Reynolds CF, Forman SD, et al. Paralimbic and lateral prefrontal encoding of reward value during intertemporal choice in attempted suicide. Psychol Med. 2016;46:381–91.26446615
14. Olié E Jollant F Deverdun J de Champfleur NM Cyprien F Le Bars E The experience of social exclusion in women with a history of suicidal acts: a neuroimaging study Sci Rep. 2017 7 89 10.1038/s41598-017-00211-x 28273888
Olié E, Jollant F, Deverdun J, de Champfleur NM, Cyprien F, Le Bars E, et al. The experience of social exclusion in women with a history of suicidal acts: a neuroimaging study. Sci Rep. 2017;7:89.28273888
15. Chase HW Segreti AM Fournier JC Phillips ML Brent D Pan L Prefrontal BOLD responses coupled to changing emotional faces in adolescents with and without a history of suicide attempt J Med Psychol 2020 22 45 54 10.3233/JMP-180018
Chase HW, Segreti AM, Fournier JC, Phillips ML, Brent D, Pan L. Prefrontal BOLD responses coupled to changing emotional faces in adolescents with and without a history of suicide attempt. J Med Psychol. 2020;22:45–54.
16. Eickhoff SB Bzdok D Laird AR Kurth F Fox PT Activation likelihood estimation meta-analysis revisited Neuroimage 2012 59 2349 61 10.1016/j.neuroimage.2011.09.017 21963913
Eickhoff SB, Bzdok D, Laird AR, Kurth F, Fox PT. Activation likelihood estimation meta-analysis revisited. Neuroimage. 2012;59:2349–61.21963913
17. Muller VI Cieslik EC Serbanescu I Laird AR Fox PT Eickhoff SB Altered brain activity in unipolar depression revisited: meta-analyses of neuroimaging studies JAMA Psychiatry 2017 74 47 55 10.1001/jamapsychiatry.2016.2783 27829086
Muller VI, Cieslik EC, Serbanescu I, Laird AR, Fox PT, Eickhoff SB. Altered brain activity in unipolar depression revisited: meta-analyses of neuroimaging studies. JAMA Psychiatry. 2017;74:47–55.27829086
18. van Heeringen K Bijttebier S Desmyter S Vervaet M Baeken C Is there a neuroanatomical basis of the vulnerability to suicidal behavior? A coordinate-based meta-analysis of structural and functional MRI studies Front Hum Neurosci 2014 8 824 25374525
van Heeringen K, Bijttebier S, Desmyter S, Vervaet M, Baeken C. Is there a neuroanatomical basis of the vulnerability to suicidal behavior? A coordinate-based meta-analysis of structural and functional MRI studies. Front Hum Neurosci. 2014;8:824.25374525
19. Zhukovsky P Anderson JAE Coughlan G Mulsant BH Cipriani A Voineskos AN Coordinate-based network mapping of brain structure in major depressive disorder in younger and older adults: a systematic review and meta-analysis Am J Psychiat 2021 178 1119 28 10.1176/appi.ajp.2021.21010088 34645274
Zhukovsky P, Anderson JAE, Coughlan G, Mulsant BH, Cipriani A, Voineskos AN. Coordinate-based network mapping of brain structure in major depressive disorder in younger and older adults: a systematic review and meta-analysis. Am J Psychiat. 2021;178:1119–28.34645274
20. Cash RFH Müller VI Fitzgerald PB Eickhoff SB Zalesky A Altered brain activity in unipolar depression unveiled using connectomics Nat Ment Health 2023 1 174 85 10.1038/s44220-023-00038-8
Cash RFH, Müller VI, Fitzgerald PB, Eickhoff SB, Zalesky A. Altered brain activity in unipolar depression unveiled using connectomics. Nat Ment Health. 2023;1:174–85.
21. Malhi GS Mann JJ Depression Lancet 2018 392 2299 312 10.1016/S0140-6736(18)31948-2 30396512
Malhi GS, Mann JJ. Depression. Lancet. 2018;392:2299–312.30396512
22. Hohwy J Functional integration and the mind Synthese 2007 159 315 28 10.1007/s11229-007-9240-3
Hohwy J. Functional integration and the mind. Synthese. 2007;159:315–28.
23. Peng SL Xu PF Jiang YY Gong GL Activation network mapping for integration of heterogeneous fMRI findings Nat Hum Behav 2022 6 1417 29 10.1038/s41562-022-01371-1 35654963
Peng SL, Xu PF, Jiang YY, Gong GL. Activation network mapping for integration of heterogeneous fMRI findings. Nat Hum Behav. 2022;6:1417–29.35654963
24. Fox MD Mapping symptoms to brain networks with the human connectome N. Engl J Med 2018 379 2237 45 10.1056/NEJMra1706158 30575457
Fox MD. Mapping symptoms to brain networks with the human connectome. N. Engl J Med. 2018;379:2237–45.30575457
25. Chang X Zhao W Kang JJ Xiang ST Xie C Corona-Hernandez H Language abnormalities in schizophrenia: binding core symptoms through contemporary empirical evidence Schizophrenia-Uk 2022 8 95 10.1038/s41537-022-00308-x
Chang X, Zhao W, Kang JJ, Xiang ST, Xie C, Corona-Hernandez H, et al. Language abnormalities in schizophrenia: binding core symptoms through contemporary empirical evidence. Schizophrenia-Uk. 2022;8:95.
26. Weil RS Hsu JK Darby RR Soussand L Fox MD Neuroimaging in Parkinson’s disease dementia: connecting the dots Brain Commun 2019 1 fcz006 10.1093/braincomms/fcz006 31608325
Weil RS, Hsu JK, Darby RR, Soussand L, Fox MD. Neuroimaging in Parkinson’s disease dementia: connecting the dots. Brain Commun. 2019;1:fcz006.31608325
27. Page MJ Moher D Bossuyt PM Boutron I Hoffmann TC Mulrow CD PRISMA 2020 explanation and elaboration: updated guidance and exemplars for reporting systematic reviews BMJ 2021 372 n160 10.1136/bmj.n160 33781993
Page MJ, Moher D, Bossuyt PM, Boutron I, Hoffmann TC, Mulrow CD, et al. PRISMA 2020 explanation and elaboration: updated guidance and exemplars for reporting systematic reviews. BMJ. 2021;372:n160.33781993
28. Darby RR Joutsa J Fox MD Network localization of heterogeneous neuroimaging findings Brain 2019 142 70 79 10.1093/brain/awy292 30551186
Darby RR, Joutsa J, Fox MD. Network localization of heterogeneous neuroimaging findings. Brain. 2019;142:70–79.30551186
29. Holmes AJ Hollinshead MO O’Keefe TM Petrov VI Fariello GR Wald LL Brain genomics superstruct project initial data release with structural, functional, and behavioral measures Sci Data 2015 2 150031 10.1038/sdata.2015.31 26175908
Holmes AJ, Hollinshead MO, O’Keefe TM, Petrov VI, Fariello GR, Wald LL, et al. Brain genomics superstruct project initial data release with structural, functional, and behavioral measures. Sci Data. 2015;2:150031.26175908
30. Yeo BTT Krienen FM Sepulcre J Sabuncu MR Lashkari D Hollinshead M The organization of the human cerebral cortex estimated by intrinsic functional connectivity J Neurophysiol 2011 106 1125 65 10.1152/jn.00338.2011 21653723
Yeo BTT, Krienen FM, Sepulcre J, Sabuncu MR, Lashkari D, Hollinshead M, et al. The organization of the human cerebral cortex estimated by intrinsic functional connectivity. J Neurophysiol. 2011;106:1125–65.21653723
31. Li H Chen Z Gong Q Jia Z Voxel-wise meta-analysis of task-related brain activation abnormalities in major depressive disorder with suicide behavior Brain Imaging Behav 2020 14 1298 308 10.1007/s11682-019-00045-3 30790165
Li H, Chen Z, Gong Q, Jia Z. Voxel-wise meta-analysis of task-related brain activation abnormalities in major depressive disorder with suicide behavior. Brain Imaging Behav. 2020;14:1298–308.30790165
32. Rolls ET Huang CC Lin CP Feng JF Joliot M Automated anatomical labelling atlas 3 Neuroimage 2020 206 116189 10.1016/j.neuroimage.2019.116189 31521825
Rolls ET, Huang CC, Lin CP, Feng JF, Joliot M. Automated anatomical labelling atlas 3. Neuroimage. 2020;206:116189.31521825
33. Hawrylycz MJ Lein ES Guillozet-Bongaarts AL Shen EH Ng L Miller JA An anatomically comprehensive atlas of the adult human brain transcriptome Nature 2012 489 391 9 10.1038/nature11405 22996553
Hawrylycz MJ, Lein ES, Guillozet-Bongaarts AL, Shen EH, Ng L, Miller JA, et al. An anatomically comprehensive atlas of the adult human brain transcriptome. Nature. 2012;489:391–9.22996553
34. Markello RD Arnatkeviciute A Poline JB Fulcher B Fornito A Misic B Standardizing workflows in imaging transcriptomics with the abagen toolbox Elife 2021 10 e72129 10.7554/eLife.72129 34783653
Markello RD, Arnatkeviciute A, Poline JB, Fulcher B, Fornito A, Misic B. Standardizing workflows in imaging transcriptomics with the abagen toolbox. Elife. 2021;10:e72129.34783653
35. Arnatkeviciute A Fulcher BD Fornito A A practical guide to linking brain-wide gene expression and neuroimaging data Neuroimage 2019 189 353 67 10.1016/j.neuroimage.2019.01.011 30648605
Arnatkeviciute A, Fulcher BD, Fornito A. A practical guide to linking brain-wide gene expression and neuroimaging data. Neuroimage. 2019;189:353–67.30648605
36. Fan LZ Li H Zhuo JJ Zhang Y Wang JJ Chen LF The human brainnetome atlas: a new brain atlas based on connectional architecture Cereb Cortex 2016 26 3508 26 10.1093/cercor/bhw157 27230218
Fan LZ, Li H, Zhuo JJ, Zhang Y, Wang JJ, Chen LF, et al. The human brainnetome atlas: a new brain atlas based on connectional architecture. Cereb Cortex. 2016;26:3508–26.27230218
37. Abdi H Partial least squares regression and projection on latent structure regression (PLS Regression) WIREs Comput Stat 2010 2 97 106 10.1002/wics.51
Abdi H. Partial least squares regression and projection on latent structure regression (PLS Regression). WIREs Comput Stat. 2010;2:97–106.
38. Wagner HH Dray S Generating spatially constrained null models for irregularly spaced data using Moran spectral randomization methods Methods Ecol Evol 2015 6 1169 78 10.1111/2041-210X.12407
Wagner HH, Dray S. Generating spatially constrained null models for irregularly spaced data using Moran spectral randomization methods. Methods Ecol Evol. 2015;6:1169–78.
39. de Wael RV Benkarim O Paquola C Lariviere S Royer J Tavakol S BrainSpace: a toolbox for the analysis of macroscale gradients in neuroimaging and connectomics datasets Commun Biol 2020 3 103 10.1038/s42003-020-0794-7 32139786
de Wael RV, Benkarim O, Paquola C, Lariviere S, Royer J, Tavakol S, et al. BrainSpace: a toolbox for the analysis of macroscale gradients in neuroimaging and connectomics datasets. Commun Biol. 2020;3:103.32139786
40. Whitaker KJ Vertes PE Romero-Garcia R Vasa F Moutoussis M Prabhu G Adolescence is associated with genomically patterned consolidation of the hubs of the human brain connectome P Natl Acad Sci USA 2016 113 9105 10 10.1073/pnas.1601745113
Whitaker KJ, Vertes PE, Romero-Garcia R, Vasa F, Moutoussis M, Prabhu G, et al. Adolescence is associated with genomically patterned consolidation of the hubs of the human brain connectome. P Natl Acad Sci USA. 2016;113:9105–10.
41. Zhou YY Zhou B Pache L Chang M Khodabakhshi AH Tanaseichuk O Metascape provides a biologist-oriented resource for the analysis of systems-level datasets Nat Commun 2019 10 1523 10.1038/s41467-019-09234-6 30944313
Zhou YY, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10:1523.30944313
42. Li J Seidlitz J Suckling J Fan FY Ji GJ Meng Y Cortical structural differences in major depressive disorder correlate with cell type-specific transcriptional signatures Nat Commun 2021 12 1647 10.1038/s41467-021-21943-5 33712584
Li J, Seidlitz J, Suckling J, Fan FY, Ji GJ, Meng Y, et al. Cortical structural differences in major depressive disorder correlate with cell type-specific transcriptional signatures. Nat Commun. 2021;12:1647.33712584
43. Hansen JY Shafiei G Markello RD Smart K Cox SML Norgaard M Mapping neurotransmitter systems to the structural and functional organization of the human neocortex Nat Neurosci 2022 25 1569 81 10.1038/s41593-022-01186-3 36303070
Hansen JY, Shafiei G, Markello RD, Smart K, Cox SML, Norgaard M, et al. Mapping neurotransmitter systems to the structural and functional organization of the human neocortex. Nat Neurosci. 2022;25:1569–81.36303070
44. Werner FM Covenas R Classical neurotransmitters and neuropeptides involved in major depression: a review Int J Neurosci 2010 120 455 70 10.3109/00207454.2010.483651 20583898
Werner FM, Covenas R. Classical neurotransmitters and neuropeptides involved in major depression: a review. Int J Neurosci. 2010;120:455–70.20583898
45. Gromping U Relative importance for linear regression in R: The package relaimpo J Stat Softw 2006 17 1 27 10.18637/jss.v017.i01
Gromping U. Relative importance for linear regression in R: The package relaimpo. J Stat Softw. 2006;17:1–27.
46. Lu QY Towne JM Lock M Jiang C Cheng ZX Habes M Toward coordinate-based cognition dictionaries: a brainMap and neurosynth demo Neuroscience 2022 493 109 18 10.1016/j.neuroscience.2022.02.016 35569642
Lu QY, Towne JM, Lock M, Jiang C, Cheng ZX, Habes M, et al. Toward coordinate-based cognition dictionaries: a brainMap and neurosynth demo. Neuroscience. 2022;493:109–18.35569642
47. Yarkoni T Poldrack RA Nichols TE Van Essen DC Wager TD Large-scale automated synthesis of human functional neuroimaging data Nat Methods 2011 8 665 U695 10.1038/nmeth.1635 21706013
Yarkoni T, Poldrack RA, Nichols TE, Van Essen DC, Wager TD. Large-scale automated synthesis of human functional neuroimaging data. Nat Methods. 2011;8:665–U695.21706013
48. Margulies DS Ghosh SS Goulas A Falkiewicz M Huntenburg JM Langs G Situating the default-mode network along a principal gradient of macroscale cortical organization P Natl Acad Sci USA 2016 113 12574 9 10.1073/pnas.1608282113
Margulies DS, Ghosh SS, Goulas A, Falkiewicz M, Huntenburg JM, Langs G, et al. Situating the default-mode network along a principal gradient of macroscale cortical organization. P Natl Acad Sci USA. 2016;113:12574–9.
49. Morgan SE Seidlitz J Whitaker KJ Romero-Garcia R Clifton NE Scarpazza C Cortical patterning of abnormal morphometric similarity in psychosis is associated with brain expression of schizophrenia-related genes Proc Natl Acad Sci USA 2019 116 9604 9 10.1073/pnas.1820754116 31004051
Morgan SE, Seidlitz J, Whitaker KJ, Romero-Garcia R, Clifton NE, Scarpazza C, et al. Cortical patterning of abnormal morphometric similarity in psychosis is associated with brain expression of schizophrenia-related genes. Proc Natl Acad Sci USA. 2019;116:9604–9.31004051
50. Taylor JJ Lin C Talmasov D Ferguson MA Schaper F Jiang J A transdiagnostic network for psychiatric illness derived from atrophy and lesions Nat Hum Behav 2023 7 420 9 10.1038/s41562-022-01501-9 36635585
Taylor JJ, Lin C, Talmasov D, Ferguson MA, Schaper F, Jiang J, et al. A transdiagnostic network for psychiatric illness derived from atrophy and lesions. Nat Hum Behav. 2023;7:420–9.36635585
51. Ji GJ Zalesky A Wang YR He KL Wang L Du RR Linking personalized brain atrophy to Schizophrenia network and treatment response Schizophrenia Bull 2023 49 43 52 10.1093/schbul/sbac162
Ji GJ, Zalesky A, Wang YR, He KL, Wang L, Du RR, et al. Linking personalized brain atrophy to Schizophrenia network and treatment response. Schizophrenia Bull. 2023;49:43–52.
52. Siddiqi SH Schaper FLWVJ Horn A Hsu J Padmanabhan JL Brodtmann A Brain stimulation and brain lesions converge on common causal circuits in neuropsychiatric disease Nat Hum Behav 2021 5 1707 16 10.1038/s41562-021-01161-1 34239076
Siddiqi SH, Schaper FLWVJ, Horn A, Hsu J, Padmanabhan JL, Brodtmann A, et al. Brain stimulation and brain lesions converge on common causal circuits in neuropsychiatric disease. Nat Hum Behav. 2021;5:1707–16.34239076
53. Gordon EM Chauvin RJ Van AN Rajesh A Nielsen A Newbold DJ A somato-cognitive action network alternates with effector regions in motor cortex Nature 2023 617 351 9 10.1038/s41586-023-05964-2 37076628
Gordon EM, Chauvin RJ, Van AN, Rajesh A, Nielsen A, Newbold DJ, et al. A somato-cognitive action network alternates with effector regions in motor cortex. Nature. 2023;617:351–9.37076628
54. Bari A Robbins TW Inhibition and impulsivity: behavioral and neural basis of response control Prog Neurobiol 2013 108 44 79 10.1016/j.pneurobio.2013.06.005 23856628
Bari A, Robbins TW. Inhibition and impulsivity: behavioral and neural basis of response control. Prog Neurobiol. 2013;108:44–79.23856628
55. Albanese BJ Macatee RJ Gallyer AJ Stanley IH Joiner TE Schmidt NB Impaired conflict detection differentiates suicide attempters from ideating nonattempters: evidence from event-related potentials Biol Psychiatry Cogn Neurosci Neuroimaging 2019 4 902 12 31399392
Albanese BJ, Macatee RJ, Gallyer AJ, Stanley IH, Joiner TE, Schmidt NB. Impaired conflict detection differentiates suicide attempters from ideating nonattempters: evidence from event-related potentials. Biol Psychiatry Cogn Neurosci Neuroimaging. 2019;4:902–12.31399392
56. Hwang JP Lee TW Tsai SJ Chen TJ Yang CH Lirng JF Cortical and subcortical abnormalities in late-onset depression with history of suicide attempts investigated with MRI and voxel-based morphometry J Geriatr Psychiatry Neurol 2010 23 171 84 10.1177/0891988710363713 20430976
Hwang JP, Lee TW, Tsai SJ, Chen TJ, Yang CH, Lirng JF, et al. Cortical and subcortical abnormalities in late-onset depression with history of suicide attempts investigated with MRI and voxel-based morphometry. J Geriatr Psychiatry Neurol. 2010;23:171–84.20430976
57. Tsujii N Mikawa W Tsujimoto E Adachi T Niwa A Ono H Reduced left precentral regional responses in patients with major depressive disorder and history of suicide attempts Plos One 2017 12 e0175249 10.1371/journal.pone.0175249 28380030
Tsujii N, Mikawa W, Tsujimoto E, Adachi T, Niwa A, Ono H, et al. Reduced left precentral regional responses in patients with major depressive disorder and history of suicide attempts. Plos One. 2017;12:e0175249.28380030
58. King JS Sad mood, emotion regulation, and response inhibition J Cogn Psychol 2020 32 573 9 10.1080/20445911.2020.1777418
King JS. Sad mood, emotion regulation, and response inhibition. J Cogn Psychol. 2020;32:573–9.
59. Roxo MR Franceschini PR Zubaran C Kleber FD Sander JW The limbic system conception and its historical evolution ScientificWorldJournal 2011 11 2427 40 10.1100/2011/157150
Roxo MR, Franceschini PR, Zubaran C, Kleber FD, Sander JW. The limbic system conception and its historical evolution. ScientificWorldJournal. 2011;11:2427–40.
60. Reinelt J Uhlig M Muller K Lauckner ME Kumral D Schaare HL Acute psychosocial stress alters thalamic network centrality Neuroimage 2019 199 680 90 10.1016/j.neuroimage.2019.06.005 31173902
Reinelt J, Uhlig M, Muller K, Lauckner ME, Kumral D, Schaare HL, et al. Acute psychosocial stress alters thalamic network centrality. Neuroimage. 2019;199:680–90.31173902
61. Young KA Bonkale WL Holcomb LA Hicks PB German DC Major depression, 5HTTLPR genotype, suicide and antidepressant influences on thalamic volume Br J Psychiatry 2008 192 285 9 10.1192/bjp.bp.107.039180 18378990
Young KA, Bonkale WL, Holcomb LA, Hicks PB, German DC. Major depression, 5HTTLPR genotype, suicide and antidepressant influences on thalamic volume. Br J Psychiatry. 2008;192:285–9.18378990
62. Berardelli I Serafini G Cortese N Fiasche F O’Connor RC Pompili M The involvement of Hypothalamus-Pituitary-Adrenal (HPA) axis in suicide risk Brain Sci. 2020 10 653 10.3390/brainsci10090653 32967089
Berardelli I, Serafini G, Cortese N, Fiasche F, O’Connor RC, Pompili M. The involvement of Hypothalamus-Pituitary-Adrenal (HPA) axis in suicide risk. Brain Sci. 2020;10:653.32967089
63. Lemke H Probst S Warneke A Waltemate L Winter A Thiel K The course of disease in major depressive disorder is associated with altered activity of the limbic system during negative emotion processing Biol Psychiatry-Cogn N. 2022 7 323 32
Lemke H, Probst S, Warneke A, Waltemate L, Winter A, Thiel K, et al. The course of disease in major depressive disorder is associated with altered activity of the limbic system during negative emotion processing. Biol Psychiatry-Cogn N. 2022;7:323–32.
64. Balcioglu YH Kose S Neural substrates of suicide and suicidal behaviour: from a neuroimaging perspective Psychiatry Clin Psych 2018 28 314 28
Balcioglu YH, Kose S. Neural substrates of suicide and suicidal behaviour: from a neuroimaging perspective. Psychiatry Clin Psych. 2018;28:314–28.
65. Malejko K Brown RC Plener PL Bonenberger M Graf H Abler B Differential neural processing of unpleasant sensory stimulation in patients with major depression Eur Arch Psy Clin N. 2021 271 557 65
Malejko K, Brown RC, Plener PL, Bonenberger M, Graf H, Abler B. Differential neural processing of unpleasant sensory stimulation in patients with major depression. Eur Arch Psy Clin N. 2021;271:557–65.
66. Ohayon MM Schatzberg AF Chronic pain and major depressive disorder in the general population J Psychiatr Res 2010 44 454 61 10.1016/j.jpsychires.2009.10.013 20149391
Ohayon MM, Schatzberg AF. Chronic pain and major depressive disorder in the general population. J Psychiatr Res. 2010;44:454–61.20149391
67. Park MJ Choi KW Na EJ Hong JP Cho MJ Fava M Multiple types of somatic pain increase suicide attempts and depression: A nationwide community sample of Korean adults Compr Psychiatry 2019 90 43 48 10.1016/j.comppsych.2018.12.006 30684832
Park MJ, Choi KW, Na EJ, Hong JP, Cho MJ, Fava M, et al. Multiple types of somatic pain increase suicide attempts and depression: A nationwide community sample of Korean adults. Compr Psychiatry. 2019;90:43–48.30684832
68. Javaheripour N Li M Chand T Krug A Kircher T Dannlowski U Altered resting-state functional connectome in major depressive disorder: a mega-analysis from the PsyMRI consortium Transl Psychiatry 2021 11 511 10.1038/s41398-021-01619-w 34620830
Javaheripour N, Li M, Chand T, Krug A, Kircher T, Dannlowski U, et al. Altered resting-state functional connectome in major depressive disorder: a mega-analysis from the PsyMRI consortium. Transl Psychiatry. 2021;11:511.34620830
69. Serafini G Pompili M Seretti ME Stefani H Palermo M Coryell W The role of inflammatory cytokines in suicidal behavior: A systematic review Eur Neuropsychopharm 2013 23 1672 86 10.1016/j.euroneuro.2013.06.002
Serafini G, Pompili M, Seretti ME, Stefani H, Palermo M, Coryell W, et al. The role of inflammatory cytokines in suicidal behavior: A systematic review. Eur Neuropsychopharm. 2013;23:1672–86.
70. Leonard BE Changes in the immune system in depression and dementia: causal or co-incidental effects? Int J Dev Neurosci 2001 19 305 12 10.1016/S0736-5748(01)00014-4 11337199
Leonard BE. Changes in the immune system in depression and dementia: causal or co-incidental effects? Int J Dev Neurosci. 2001;19:305–12.11337199
71. Yagci I Avci S Biochemical predictors in presentations to the emergency department after a suicide attempt Bratisl Med J 2021 122 224 9 10.4149/BLL_2021_012
Yagci I, Avci S. Biochemical predictors in presentations to the emergency department after a suicide attempt. Bratisl Med J. 2021;122:224–9.
72. Lu K Hong Y Tao M Shen L Zheng Z Fang K Depressive patient-derived GABA interneurons reveal abnormal neural activity associated with HTR2C EMBO Mol Med 2023 15 e16364 10.15252/emmm.202216364 36373384
Lu K, Hong Y, Tao M, Shen L, Zheng Z, Fang K, et al. Depressive patient-derived GABA interneurons reveal abnormal neural activity associated with HTR2C. EMBO Mol Med. 2023;15:e16364.36373384
73. Moret C Briley M The importance of norepinephrine in depression Neuropsychiatr Dis Treat 2011 7 9 13 21750623
Moret C, Briley M. The importance of norepinephrine in depression. Neuropsychiatr Dis Treat. 2011;7:9–13.21750623
74. Stahl SM. Stahl’s essential psychopharmacology neuroscientific basis and practical applications. 5th Edition, 2021.
75. Baune BT Miller R McAfoose J Johnson M Quirk F Mitchell D The role of cognitive impairment in general functioning in major depression Psychiatry Res 2010 176 183 9 10.1016/j.psychres.2008.12.001 20138370
Baune BT, Miller R, McAfoose J, Johnson M, Quirk F, Mitchell D. The role of cognitive impairment in general functioning in major depression. Psychiatry Res. 2010;176:183–9.20138370
76. Rock PL Roiser JP Riedel WJ Blackwell AD Cognitive impairment in depression: a systematic review and meta-analysis Psychol Med 2014 44 2029 40 10.1017/S0033291713002535 24168753
Rock PL, Roiser JP, Riedel WJ, Blackwell AD. Cognitive impairment in depression: a systematic review and meta-analysis. Psychol Med. 2014;44:2029–40.24168753
77. Gujral S Dombrovski AY Butters M Clark L Reynolds CF Szanto K Impaired executive function in contemplated and attempted suicide in late life Am J Geriat Psychiatry 2014 22 811 9 10.1016/j.jagp.2013.01.025
Gujral S, Dombrovski AY, Butters M, Clark L, Reynolds CF, Szanto K. Impaired executive function in contemplated and attempted suicide in late life. Am J Geriat Psychiatry. 2014;22:811–9.
78. Keilp JG Gorlyn M Oquendo MA Burke AK Mann JJ Attention deficit in depressed suicide attempters Psychiatry Res 2008 159 7 17 10.1016/j.psychres.2007.08.020 18329724
Keilp JG, Gorlyn M, Oquendo MA, Burke AK, Mann JJ. Attention deficit in depressed suicide attempters. Psychiatry Res. 2008;159:7–17.18329724
79. Siddiqi SH Kording KP Parvizi J Fox MD Causal mapping of human brain function Nat Rev Neurosci 2022 23 361 75 10.1038/s41583-022-00583-8 35444305
Siddiqi SH, Kording KP, Parvizi J, Fox MD. Causal mapping of human brain function. Nat Rev Neurosci. 2022;23:361–75.35444305
80. Siddiqi SH, Khosravani S, Rolston JD, Fox MD. The future of brain circuit-targeted therapeutics. Neuropsychopharmacology. 2024;49:179–88.
