Hostname: page-component-7b9c58cd5d-9klzr Total loading time: 0 Render date: 2025-03-16T06:15:43.381Z Has data issue: false hasContentIssue false

Local dynamic spontaneous brain activity changes in first-episode, treatment-naïve patients with major depressive disorder and their associated gene expression profiles

Published online by Cambridge University Press:  30 October 2020

Kaizhong Xue
Affiliation:
Department of Radiology and Tianjin Key Laboratory of Functional Imaging, Tianjin Medical University General Hospital, Tianjin 300052, China
Sixiang Liang
Affiliation:
Tianjin Anding Hospital, Tianjin 300222, China The National Clinical Research Center for Mental Disorders & Beijing Key Laboratory of Mental Disorders, Beijing Anding Hospital, Capital Medical University & the Advanced Innovation Center for Human Brain Protection, Capital Medical University, Beijing 100088, China
Bingbing Yang
Affiliation:
Department of Radiology and Tianjin Key Laboratory of Functional Imaging, Tianjin Medical University General Hospital, Tianjin 300052, China
Dan Zhu
Affiliation:
Department of Radiology and Tianjin Key Laboratory of Functional Imaging, Tianjin Medical University General Hospital, Tianjin 300052, China
Yingying Xie
Affiliation:
Department of Radiology and Tianjin Key Laboratory of Functional Imaging, Tianjin Medical University General Hospital, Tianjin 300052, China
Wen Qin
Affiliation:
Department of Radiology and Tianjin Key Laboratory of Functional Imaging, Tianjin Medical University General Hospital, Tianjin 300052, China
Feng Liu*
Affiliation:
Department of Radiology and Tianjin Key Laboratory of Functional Imaging, Tianjin Medical University General Hospital, Tianjin 300052, China
Yong Zhang*
Affiliation:
Tianjin Anding Hospital, Tianjin 300222, China
Chunshui Yu*
Affiliation:
Department of Radiology and Tianjin Key Laboratory of Functional Imaging, Tianjin Medical University General Hospital, Tianjin 300052, China CAS Center for Excellence in Brain Science and Intelligence Technology, Chinese Academy of Sciences, Shanghai 200031, China
*
Author for correspondence: Feng Liu, E-mail: fengliu@tijmu.edu.cn; Yong Zhang, E-mail: zhangyong@tjmhc.com; Chunshui Yu, E-mail: chunshuiyu@tijmu.edu.cn
Author for correspondence: Feng Liu, E-mail: fengliu@tijmu.edu.cn; Yong Zhang, E-mail: zhangyong@tjmhc.com; Chunshui Yu, E-mail: chunshuiyu@tijmu.edu.cn
Author for correspondence: Feng Liu, E-mail: fengliu@tijmu.edu.cn; Yong Zhang, E-mail: zhangyong@tjmhc.com; Chunshui Yu, E-mail: chunshuiyu@tijmu.edu.cn
Rights & Permissions [Opens in a new window]

Abstract

Background

Major depressive disorder (MDD) is a common debilitating disorder characterized by impaired spontaneous brain activity, yet little is known about its alterations in dynamic properties and the molecular mechanisms associated with these changes.

Methods

Based on the resting-state functional MRI data of 65 first-episode, treatment-naïve patients with MDD and 66 healthy controls, we compared dynamic regional homogeneity (dReHo) of spontaneous brain activity between the two groups, and we investigated gene expression profiles associated with dReHo alterations in MDD by leveraging transcriptional data from the Allen Human Brain Atlas and weighted gene co-expression network analysis.

Results

Compared with healthy controls, patients with MDD consistently showed reduced dReHo in both fusiform gyri and in the right temporal pole and hippocampus. The expression profiles of 16 gene modules were correlated with dReHo alterations in MDD. These gene modules were enriched for various biological process terms, including immune, synaptic signalling, ion channels, mitochondrial function and protein metabolism, and were preferentially expressed in different cell types.

Conclusions

Patients with MDD have reduced dReHo in brain areas associated with emotional and cognitive regulation, and these changes may be related to complex polygenetic and polypathway mechanisms.

Type
Original Article
Copyright
Copyright © The Author(s) 2020. Published by Cambridge University Press

Introduction

Major depressive disorder (MDD) is a leading cause of disability worldwide, affecting 350 million people each year (World Health Organization, 2008). The conditions of more than half of patients with MDD become chronic or recurrent after the first depressive episode, resulting in huge economic and medical burdens (Cuijpers, Beekman, & Reynolds, Reference Cuijpers, Beekman and Reynolds2012). Although substantial efforts have been made in the past decade, the pathophysiological mechanism of MDD remains largely unknown.

Resting-state functional magnetic resonance imaging (rs-fMRI) is increasingly emerging as a promising tool to delineate spontaneous brain activity of the human brain, which has prominent advantages in clinical research due to the minimal need for cooperation from patients. Previous studies have demonstrated altered resting-state spontaneous brain activity in MDD via various approaches (Fitzgerald, Laird, Maller, & Daskalakis, Reference Fitzgerald, Laird, Maller and Daskalakis2008; Geng et al., Reference Geng, Yan, Shi, Chen, Mo, Shao and Lu2019; Guo et al., Reference Guo, Liu, Xue, Yu, Ma, Tan and Zhao2011a, Reference Guo, Sun, Liu, Xu, Wu, Liu and Zhao2011b; Xia et al., Reference Xia, Si, Sun, Ma, Liu, Wang and Group2019), among which regional homogeneity (ReHo) shows high reproducibility (Zang, Jiang, Lu, He, & Tian, Reference Zang, Jiang, Lu, He and Tian2004) and has been repeatedly reported to be changed in MDD (Lai, Reference Lai2018; Sun et al., Reference Sun, Luo, Yuan, Zhang, He, Yao and Xiao2018; Xia et al., Reference Xia, Si, Sun, Ma, Liu, Wang and Group2019). However, ReHo can reflect only static temporal synchronization of spontaneous brain activity among nearby voxels, which is inconsistent with the concept that spontaneous brain activity is a time-dependent dynamic process (Hutchison et al., Reference Hutchison, Womelsdorf, Allen, Bandettini, Calhoun, Corbetta and Chang2013). In recent years, dynamic functional connectivity has been used to assess large-scale dynamic patterns of temporal synchronization of spontaneous brain activity between spatially remote brain regions/voxels (Friston, Reference Friston2011; Hutchison et al., Reference Hutchison, Womelsdorf, Allen, Bandettini, Calhoun, Corbetta and Chang2013), and its changes have been reported in neuropsychiatric diseases, including schizophrenia (Sakoglu et al., Reference Sakoglu, Pearlson, Kiehl, Wang, Michael and Calhoun2010), Alzheimer's disease (de Vos et al., Reference de Vos, Koini, Schouten, Seiler, van der Grond, Lechner and Rombouts2018), epilepsy (Liu et al., Reference Liu, Wang, Li, Wang, Li, Zhang and Chen2017), and MDD (Demirtaş et al., Reference Demirtaş, Tornador, Falcón, López-Solà, Hernández-Ribas, Pujol and Deco2016; Han et al., Reference Han, Wang, He, Sheng, Zou, Li and Chen2019; Kaiser et al., Reference Kaiser, Whitfield-Gabrieli, Dillon, Goer, Beltzer, Minkel and Pizzagalli2016; Qiu et al., Reference Qiu, Xia, Cheng, Yuan, Kuang, Bi and Gong2018). However, few studies have focused on the dynamic properties of regional temporal synchronization of spontaneous brain activity among nearby voxels, which can be measured by dynamic ReHo (dReHo). Furthermore, no studies have investigated dReHo changes in MDD.

MDD has an estimated heritability of approximately 37% (Sullivan, Neale, & Kendler, Reference Sullivan, Neale and Kendler2000); however, reproducible risk loci of MDD rarely have been identified in previous genome-wide association studies (GWASs) (Howard et al., Reference Howard, Adams, Shirali, Clarke, Marioni, Davies and McIntosh2018, Reference Howard, Adams, Clarke, Hafferty, Gibson, Shirali and McIntosh2019, Sullivan et al., Reference Sullivan, de Geus, Willemsen, James, Smit, Zandbelt and Penninx2009; Wray et al., Reference Wray, Ripke, Mattheisen, Trzaskowski, Byrne and Abdellaoui2018), possibly due to the inadequate sample size. As intermediate phenotypes, neuroimaging features of MDD are theoretically closer to genetic substrates of MDD (Gottesman & Gould, Reference Gottesman and Gould2003). In addition, gene expression is theoretically closer to neuroimaging features of MDD than genetic variants. Therefore, several studies have used cross-sample or cross-region spatial correlation analysis to investigate associations between brain gene expression and neuroimaging features of major psychiatric disorders and have achieved great success (Liu, Tian, Li, Li, & Zhuo, Reference Liu, Tian, Li, Li and Zhuo2019; Morgan et al., Reference Morgan, Seidlitz, Whitaker, Romero-Garcia, Clifton, Scarpazza and Bullmore2019; Romero-Garcia, Warrier, Bullmore, Baron-Cohen, & Bethlehem, Reference Romero-Garcia, Warrier, Bullmore, Baron-Cohen and Bethlehem2019; Romme, de Reus, Ophoff, Kahn, & van den Heuvel, Reference Romme, de Reus, Ophoff, Kahn and van den Heuvel2017). Nevertheless, no transcription-neuroimaging association studies have been carried out to identify genes whose brain expression profiles are correlated with dReHo changes in MDD.

We hypothesize that brain gene expression is related to regional spontaneous brain activity changes in MDD. Here, gene expression data were derived from six postmortem brains provided by the Allen Human Brain Atlas (AHBA) (Hawrylycz et al., Reference Hawrylycz, Miller, Menon, Feng, Dolbeare, Guillozet-Bongaarts and Lein2015), and regional spontaneous brain activity was assessed by dReHo since it is more consistent with the fact that brain activity is dynamic (Hutchison et al., Reference Hutchison, Womelsdorf, Allen, Bandettini, Calhoun, Corbetta and Chang2013). We then identified genes associated with dReHo changes in MDD by performing cross-sample spatial correlations between gene expression data from the postmortem AHBA brains and case-control dReHo differences from living human brains. Finally, the function of the identified genes associated with dReHo changes in MDD was explored by several enrichment analyses. An overview of the analytical framework is summarized in Fig. 1.

Fig. 1. Overview of the analysis pipeline. (a) a sliding window method was utilized to compute dReHo, and between-group comparisons of dReHo maps were performed; (b) the association between dReHo differences and gene expression from the AHBA was explored; (c) stable cortical and subcortical modules related to MDD were identified; (d) enrichment analysis was applied for these modules. dReHo, dynamic regional homogeneity; AHBA, Allen Human Brain Atlas.

Materials and methods

Subjects

Sixty-seven first-episode, treatment-naïve patients with MDD were recruited from Tianjin Medical University General Hospital. MDD was diagnosed based on the diagnostic criteria of the Structure Clinical Interview of DSM-IV (SCID). Following previous studies (Khan, Bhat, Kolts, Thase, & Brown, Reference Khan, Bhat, Kolts, Thase and Brown2010), we recruited only MDD patients with Hamilton Depression Rating Scale-17 (HDRS-17) scores ⩾18. The Hamilton Anxiety Rating Scale-14 (HARS-14) was used to evaluate anxiety symptoms in these patients, but we did not set a cut-off value. Moreover, patients with other psychiatric diseases, including schizophrenia, bipolar disorder, mental retardation and personality disorder, were excluded. A total of 69 age- and sex-matched healthy controls were recruited from the local community, and the SCID Non-Patient Edition was used to confirm that they were free of any psychiatric disorders. Common exclusion criteria for all participants were (1) age younger than 18 years or older than 65 years; (2) left-handedness; (3) the presence of organic intracranial lesions; (4) any history of comorbid alcohol or drug abuse; (5) poor image quality; and (6) MRI contradictions. The Ethics Institution of Tianjin Medical University General Hospital approved the study protocol and written informed consent was obtained from participants.

Image acquisition

All MRI images were obtained using a 3.0-Tesla Discovery MR750 scanner (General Electric, Milwaukee, WI). Foam paddings and earplugs were utilized to minimize head movement and reduce noise. The rs-fMRI data were obtained using a single-shot gradient-recalled-echo echo-planar-imaging (SS-GRE-EPI) sequence with the following parameters: repetition time (TR) = 2000 ms; echo time (TE) = 30 ms; field of view (FOV) = 220 mm × 220 mm; matrix = 64 × 64; flip angle (FA) = 90°; slice thickness = 3 mm, no gap; 36 interleaved transverse slices; and 180 volumes. During the fMRI scanning, subjects were asked to close their eyes, hold still, think of nothing in particular, and not fall asleep. Sagittal 3D T1-weighted images were acquired by a brain volume sequence using the following parameters: TR = 8.14 ms; TE = 3.17 ms; inversion time = 450 ms; FOV = 256 mm × 256 mm; matrix = 256 × 256; FA = 12°; slice thickness = 1 mm, no gap; and 188 slices.

Image preprocessing

The imaging quality of each subject was carefully examined, and subjects with poor imaging quality were excluded. The rs-fMRI data were preprocessed using the Data Processing Assistant for Resting-State fMRI (DPARSF) (Chao-Gan & Yu-Feng, Reference Chao-Gan and Yu-Feng2010) based on Statistical Parametric Mapping 12 (http://www.fil.ion.ucl.ac.uk/spm). Briefly, the first ten volumes were discarded to allow the signal to reach equilibrium and the remaining 170 volumes were corrected for temporal differences between slices and head motion. Participants with maximum displacement >2.0 mm and maximum rotation >2.0° were excluded from the subsequent analyses. In addition, the mean framewise displacement (FD) (Van Dijk, Sabuncu, & Buckner, Reference Van Dijk, Sabuncu and Buckner2012) was calculated. Subsequently, individual structural images were co-registered to the mean functional image, and then the transformed structural images were segmented into grey matter, white matter and cerebrospinal fluid. Based on these segmented images, the normalization parameters from individual native space to Montreal Neurological Institute (MNI) space were estimated by the Diffeomorphic Anatomical Registration Through Exponentiated Lie algebra (DARTEL) tool (Ashburner, Reference Ashburner2007). The motion-corrected functional imaging data were normalized to MNI space using these normalization parameters and resampled to 3-mm cubic voxels. Spurious variances were removed by regression, including linear drift, Friston-24 head motion parameters (Friston, Williams, Howard, Frackowiak, & Turner, Reference Friston, Williams, Howard, Frackowiak and Turner1996), and signals from white matter and cerebrospinal fluid. Finally, temporal bandpass filtering (0.01–0.08 Hz) was applied to reduce the effects of low-frequency drift and high-frequency noise.

dReHo calculation

dReHo was calculated using the sliding window approach with the Data Processing & Analysis of Brain Imaging (DPABI) toolbox (Yan, Wang, Zuo, & Zang, Reference Yan, Wang, Zuo and Zang2016). Since there is no consensus on the selection of window length and step size, we calculated individual dReHo maps (Fig. 1a and online Supplementary Table S1) based on 25 combinations of the two parameters (window lengths: 18, 34, 50, 66, 82 TRs; step sizes: 1, 2, 3, 4, 5 TRs) that have been used in prior studies (Hutchison & Morton, Reference Hutchison and Morton2015; Leonardi & Van De Ville, Reference Leonardi and Van De Ville2015; Li, Duan, Cui, Chen, & Liao, Reference Li, Duan, Cui, Chen and Liao2019; Li, Lu, & Yan, Reference Li, Lu and Yan2020; Liao et al., Reference Liao, Zhang, Mantini, Xu, Ji, Zhang and Lu2014). For each combination, ReHo was computed for each sliding window and then the standard deviation was calculated across these ReHo maps to characterize the dReHo. The dReHo maps were finally smoothed with a Gaussian kernel of 8-mm full width at half-maximum.

dReHo comparison

Voxel-wise two-sample t tests were performed to identify dReHo differences between the two groups in each combination (five window lengths and five step sizes), and 25 t-maps were generated while adjusting for age, sex and mean FD. Multiple comparisons were corrected by the topological false discovery rate (FDR) method (the initial height threshold: uncorrected p < 0.001; topological FDR: q < 0.05). The topological FDR is an extension of the Benjamini-Hochberg FDR (BH-FDR) and is specially designed to correct multiple comparisons in neuroimaging analyses (Chumbley, Worsley, Flandin, & Friston, Reference Chumbley, Worsley, Flandin and Friston2010). To reduce the influence of the selection of window lengths and step sizes on the dReHo differences between groups, we created a probability map to show highly consistent brain regions with significant differences in dReHo calculated by different parametric combinations. In addition, mean dReHo values of clusters with 100% probability were extracted to perform partial correlations between dReHo and clinical variables (HDRS, HARS and illness duration) in patients with MDD while adjusting for age, sex, and mean FD. Given that whole-brain coverage including the cerebellum and brain stem was not available for all subjects, we restricted between-group comparisons to within the cerebrum.

Gene expression data preprocessing

Publicly available gene expression data of six postmortem neurotypical adult brains were obtained from the AHBA database (Hawrylycz et al., Reference Hawrylycz, Lein, Guillozet-Bongaarts, Shen, Ng, Miller and Jones2012), which, in total, consists of normalized microarray expression data of more than 20 000 genes measured in 3702 brain samples. Among them, four donors had expression data only in the left hemisphere, and two donors had expression data in both hemispheres. Based on the newly proposed pipeline (Arnatkevic Iūtė, Fulcher, & Fornito, Reference Arnatkevic Iūtė, Fulcher and Fornito2019), we separately processed the two datasets and obtained two sample × gene expression matrices of 1196 × 10 185 from the four donors and 1209 × 10 185 from the two donors (Fig. 1b and online Supplementary Material). Here, the samples represent brain tissue samples, and the screening process of tissue samples is provided in online Supplementary Fig. S1.

Creation of gene modules

Weighted gene co-expression network analysis (WGCNA) is an effective tool that clusters genes into network modules and detects biologically meaningful results (Zhang & Horvath, Reference Zhang and Horvath2005). Therefore, we first constructed a signed network for each sample × gene expression matrix. Then, a WGCNA consensus network was created from two networks to identify common expression patterns across the two aforementioned expression matrices, and 19 gene modules were finally obtained (online Supplementary material, Supplementary Fig. S2 and Table S2). For each module, the expression pattern (a matrix represented by sample × number of genes in this module) was summarized by the first principal component, defined as the module eigengene (ME, number of samples × 1 vector).

Transcription-neuroimaging association analysis

To establish the sample correspondence between gene expression and neuroimaging data, we defined a 6-mm radius sphere centred at the MNI coordinate of each tissue sample and extracted the mean value on each t-statistical map of dReHo difference (25 t-statistical maps). Then, for each module, we calculated Pearson's correlation between the ME and the t-value across brain tissue samples. The correlation analysis was performed separately for the cortical and subcortical samples because of great differences in gene expression profiles between cortical and subcortical regions (Hawrylycz et al., Reference Hawrylycz, Lein, Guillozet-Bongaarts, Shen, Ng, Miller and Jones2012). Multiple comparisons were corrected using the BH-FDR method (q < 0.05). Gene modules with significant correlations between the ME and the t-value in cortical samples were defined as cortical modules, and subcortical modules were defined by the same approach. As shown in Fig. 1c, based on the results of correlation analyses, we first intersected the results of 25 combinations of dynamic parameters in cortical and subcortical modules, respectively. The cortical/subcortical modules that were correlated with all t-statistical maps derived from different combinations were considered stable modules. Then, the stable cortical and subcortical gene modules were also intersected to identify shared and distinct gene modules for cortical and subcortical regions.

Enrichment analysis

Enrichment analysis (Fig. 1d) for stable modules was mainly performed by the web server g: Profiler (https://biit.cs.ut.ee/gprofiler), including the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases and putative transcription factor binding sites (TFBS) from the TRANSFAC database (Reimand et al., Reference Reimand, Arak, Adler, Kolberg, Reisberg, Peterson and Vilo2016). All the 10 185 included genes were defined as a background gene set, and the AutoAnnotate Cytoscape application (Shannon et al., Reference Shannon, Markiel, Ozier, Baliga, Wang, Ramage and Ideker2003) was used to show the enrichment results. In addition, cell-type enrichment analysis was conducted using Fisher's exact test with a specificity index probability (pSI) (Dougherty, Schmidt, Nakajima, & Heintz, Reference Dougherty, Schmidt, Nakajima and Heintz2010) of 0.05 (BH-FDR correction) based on previously published gene expression data for different cell types of the human brain (Zhang et al., Reference Zhang, Sloan, Clarke, Caneda, Plaza, Blumenthal and Barres2016).

Results

Demographic information and clinical characteristics

After excluding two subjects (one patient and one control) due to poor image quality and three subjects (one patient and two controls) due to excessive head motion, 65 MDD patients and 66 healthy controls were finally included in this study. The demographic and clinical data of these subjects are displayed in online Supplementary Table S3. The two groups did not show significant differences in sex (chi-squared test, p = 0.78), age (two-sample t test, p = 0.65), or mean FD (two-sample t test, p = 0.31).

Case-control difference in dReHo

Based on individual dReHo maps from each of the 25 parametric combinations, we compared voxel-wise dReHo differences between patients with MDD and healthy controls (q < 0.05, topological FDR corrected), from which we created a probability map (Fig. 2) to show the probability of a given voxel with a significant intergroup difference in these comparisons. Compared with healthy controls, patients with MDD consistently (100% probability) showed significantly decreased dReHo in both fusiform gyri and the right hippocampus and temporal pole (superior part) (Fig. 2 and online Supplementary Table S4). We did not find significantly increased dReHo in any comparisons in patients with MDD relative to healthy controls. We also tested the impact of HARS score on our results, but we did not find significant effects (online Supplementary Results and Supplementary Fig. S3).

Fig. 2. A probability map showing brain voxels with significant dReHo differences (q < 0.05, topological FDR corrected) between the MDD and HC groups. The colour bar represents the probability of a given voxel showing a significant intergroup difference in the 25 comparisons. Each comparison is performed on individual dReHo maps calculated by a distinct combination of dynamic parameters. FDR, false discovery rate; HC, healthy control; L, left; MDD, major depressive disorder; R, right.

Correlations between dReHo and clinical variables

The mean dReHo values of clusters with consistent (100% probability) reduction in patients with MDD were extracted. Partial correlation analysis was conducted between the mean dReHo values and clinical variables, including HARS and HDRS scores and illness duration, while adjusting for age, sex, and FD. No significant correlations (all p > 0.05) were found between dReHo and any clinical variables in patients with MDD.

Transcription-neuroimaging associations

Nineteen consensus modules were identified from WGCNA analysis based on expression data from four left-hemispheric and two bi-hemispheric brain samples (Fig. 3). Cross-sample Pearson's correlation was calculated between MEs of each module and dReHo differences to show the consistent correlation pattern of modules across the two expression datasets and the 25 combinations. The results showed that there were 16 consistently significant MDD-related modules, 14 of which belonged to cortical modules and six of which belonged to subcortical modules (overlapping four modules) (Figs. 4 and 5a). The positive correlation means that brain regions with significant dReHo reduction in MDD show higher gene expression, and the negative correlation means that those regions show lower gene expression (online Supplementary Figs. S4 and S5).

Fig. 3. Consensus analysis of expression data. (a) WGCNA dendrogram showing consensus modules based on the topological co-expression of 10 185 genes from four left hemispheres and two whole brains. (b) To avoid the scenario that co-expression networks constructed accidentally based on the set of WGCNA parameters used, we rebuilt networks by applying different parameter sets. The final consensus modules used in this study were conserved among the majority of the tests with different parameter sets. WGCNA, weighted gene co-expression network analysis.

Fig. 4. Correlations between gene expression and dReHo changes in MDD. (a) The heatmap shows Pearson's correlation coefficients between MEs of gene modules and t values of case-control dReHo differences in cortical samples. Each column represents a gene module (n = 18), and each row represents a case-control dReHo difference map derived from dReHo maps calculated with a given combination (n = 25) of window size and step size. In the rows of the upper half, gene expression data were obtained from two AHBA donors with whole-brain coverage; in the rows of the lower half, however, they were obtained from four AHBA donors only with left-brain coverage. (b) The heatmap shows correlations in subcortical samples. * indicates that the correlation reaches the significance of q < 0.05 (BH-FDR corrected) in all analyses. The colour bar represents Pearson's correlation coefficients between MEs of gene modules and t values of case-control dReHo differences.

Fig. 5. Enrichment of genes in stable modules and hub genes. (a) Counts of MDD-related cortical and subcortical modules. (b) KEGG enrichment for genes in all modules (top two pathways shown for each module). All obtained results could pass q < 0.05 (BH-FDR correction) and showed odds ratios >1. (c) GO enrichment for genes in all modules (top two pathways shown for each module marked by different colours). (d) Cell-type marker enrichment of genes in cortical modules. The bubble size represents significance for – log10FDR (BH-FDR correction), and the bubble colour represents whether the module is significant or not. (e) The top 10 hub genes are shown for all MDD-related modules. Edges are weighted by the strength of correlation between nodes (genes). KEGG, Kyoto Encyclopedia of Genes and Genomes; GO, Gene Ontology.

Enrichment analysis

GO enrichment analysis was conducted for each preserved module (see detailed results in online Supplementary Data Table S1). Genes in these preserved modules were enriched for pathways related to synaptic signalling, adenosine triphosphate synthesis, epithelial regulation, immune response, the regulation of transcription and translation, protein localization and metabolic regulation, ion channel activities, mitochondrial function, axon ensheathment and central nervous system development. KEGG enrichment analysis revealed that genes in these preserved modules were related to the glutamatergic system, serotonergic system, neurodegenerative diseases (such as Alzheimer's disease, Parkinson disease and Huntington disease), autoimmune diseases and circadian rhythms (Figs. 5b, c, online Supplementary Tables S5 and S6). As illustrated in Fig. 5d and online Supplementary Table S7, there were seven preserved modules specifically enriched for neurons, two for astrocytes, two for microglia, one for endothelial cells and one for oligodendrocytes. Moreover, TFBS enrichment showed that 11 out of 16 preserved modules were related to transcriptional control (online Supplementary Table S8 and Data Table S1).

Discussion

To our knowledge, this is the first attempt to apply dReHo to investigate regional temporal synchronization changes in spontaneous brain activity in first-episode, treatment-naïve patients with MDD and, in parallel, to directly correlate brain gene expression with case-control dReHo differences to explore gene modules associated with dReHo changes in MDD and to characterize the involved biological pathways.

In the present study, we observed that first-episode, treatment-naïve patients with MDD consistently showed significantly decreased dReHo in the fusiform gyrus, hippocampus and temporal pole, which were present in all 25 combinations of window lengths and step sizes. The fusiform gyrus is involved in the processing of emotional faces (Kanwisher, McDermott, & Chun, Reference Kanwisher, McDermott and Chun1997; Weiner & Grill-Spector, Reference Weiner and Grill-Spector2010). Thus, our finding of reduced dReHo in the fusiform gyrus supports and extends previous reports of the structural and functional impairments in this area in MDD (Schmaal et al., Reference Schmaal, Hibar, Samann, Hall, Baune, Jahanshad and Veltman2017; Truong et al., Reference Truong, Minuzzi, Soares, Frey, Evans, MacQueen and Hall2013). As an association area, the temporal pole is involved in various high-level cognitive functions, such as facial recognition (Olsson & Ochsner, Reference Olsson and Ochsner2008), memory (Munoz-Lopez, Mohedano-Moriano, & Insausti, Reference Munoz-Lopez, Mohedano-Moriano and Insausti2010), and language processing (Hickok & Poeppel, Reference Hickok and Poeppel2007). This structure also plays a role in emotion regulation and the theory of mind mainly by interacting with the amygdala and orbital frontal cortex (Frith & Frith, Reference Frith and Frith2010), which also showed a high probability of having decreased dReHo in patients with MDD. Moreover, this result extends the finding of the reduced static ReHo in MDD (Guo et al., Reference Guo, Sun, Liu, Xu, Wu, Liu and Zhao2011b; Liu et al., Reference Liu, Xu, Xu, Wang, Zhao, Lv and Du2010) to the impairment of regional dynamics of temporal synchronization of spontaneous brain activity in the disorder. It has been widely accepted that exposure to life stress is a major risk factor for MDD (McEwen, Reference McEwen2008), which is related to the fact that overwhelming stress can impair normal hippocampal function. An increasing number of studies have demonstrated that the hippocampus is highly sensitive to stress (Pittenger & Duman, Reference Pittenger and Duman2008), the most important environmental risk factor for the development of MDD (Kessler, Reference Kessler1997) and it is also the main target associated with the genetic risk of MDD (Xu et al., Reference Xu, Li, Qin, Jun Li, Zhuo, Liu and Yu2018). This is why our study and others always show hippocampal impairment in patients with MDD (Gold, Reference Gold2015; Liu et al., Reference Liu, Guo, Liu, Long, Ma, Xue and Chen2013). Taking previous and our findings together, we propose that functional impairment in the cognition- and emotion-regulatory networks is an important feature of MDD.

Unexpectedly, we did not find any significant correlations between dReHo values in the consistently reduced dReHo regions and clinical variables (HARS and HDRS scores and illness duration) in patients with MDD. The reasons for the non-significant findings are elusive, but they may be related to the clinical features of patients (first-episode, treatment-naïve), the lack of clinical assessments (such as rumination) reflecting other aspects of MDD, the relatively short illness duration, and others. Thus, the clinical value of dReHo changes in MDD needs to be clarified in future studies.

In contrast to previous GWASs that failed to identify reproducible risk loci of MDD (Howard et al., Reference Howard, Adams, Shirali, Clarke, Marioni, Davies and McIntosh2018, Reference Howard, Adams, Clarke, Hafferty, Gibson, Shirali and McIntosh2019; Sullivan et al., Reference Sullivan, de Geus, Willemsen, James, Smit, Zandbelt and Penninx2009; Wray et al., Reference Wray, Ripke, Mattheisen, Trzaskowski, Byrne and Abdellaoui2018), by combining an endophenotype (dReHo) of MDD and brain gene expression data, we found a dozen gene modules with various functions associated with dReHo changes in MDD, indicating that MDD may develop via complex polygenetic and polypathway mechanisms. Gene expression profiles in both cortical and subcortical areas, in multiple cell types, in different biological functions and pathways might be related to the development of MDD. Gene expression in both cortical and subcortical areas was associated with dReHo changes in MDD, mimicking the pattern of the brain structural and functional changes in this disorder (Chen et al., Reference Chen, Liu, Xun, Chen, Hu, Guo and Zhao2012; Liu et al., Reference Liu, Guo, Liu, Long, Ma, Xue and Chen2013; Ma et al., Reference Ma, Ding, Li, Guo, Long, Liu and Chen2012; Vasic, Walter, Hose, & Wolf, Reference Vasic, Walter, Hose and Wolf2008). The co-existence of the shared and specific gene modules between cortical and subcortical areas reflects the complexity of the molecular substrates of MDD. This complexity was also supported by the involvement of the related gene modules in various cell types, biological functions and pathways. For example, several modules enriched for neurons were mainly implicated in synaptic transmission and neuronal development, indicating that neuronal impairment and maldevelopment may be possible mechanisms of MDD (Krishnan & Nestler, Reference Krishnan and Nestler2008; Willner, Scheel-Kruger, & Belzung, Reference Willner, Scheel-Kruger and Belzung2013). Several modules enriched for microglia are mainly implicated in the immune system, which is consistent with the stress model of MDD via microglial activation and the inflammatory response (Felger & Lotrich, Reference Felger and Lotrich2013; Hodes, Kana, Menard, Merad, & Russo, Reference Hodes, Kana, Menard, Merad and Russo2015; Setiawan et al., Reference Setiawan, Wilson, Mizrahi, Rusjan, Miler, Rajkowska and Meyer2015). Since dReHo is calculated based on blood-oxygen-level-dependent (BOLD) signals, the association of dReHo changes in MDD with the expression profiles of gene modules enriched for neurons, astrocytes and endothelial cells is consistent with the neurovascular coupling theory of BOLD signals because these cells are the main components of the neurovascular unit (Iadecola, Reference Iadecola2017; Muoio, Persson, & Sendeski, Reference Muoio, Persson and Sendeski2014). These findings indicate that impaired neurovascular coupling may also contribute to depressive symptoms. In addition, the enrichment of MDD-related genes for neurodegenerative diseases may be a reflection of the shared assaulting target (such as hippocampus) by the risk factors for these disorders (Xu et al., Reference Xu, Li, Qin, Jun Li, Zhuo, Liu and Yu2018). Moreover, hub gene detection can provide candidate genes for further studies. For example, several immune-related genes were identified by hub gene detection, which may be used to explore the immune mechanism of MDD.

At least two limitations should be considered when interpreting our findings. The first limitation is that gene expression data and brain imaging data were not obtained from the same subjects. The identified correlations between gene expression data from six postmortem brains and brain imaging data from other living human brains may be correct only if gene expression profiles are highly conserved across individuals. That is, this study would inevitably miss genes or gene modules with great individual differences. Since gene expression profiles are dependent on age, sex, and ethnicity (Mogil et al., Reference Mogil, Andaleon, Badalamenti, Dickinson, Guo, Rotter and Wheeler2018; Somel et al., Reference Somel, Guo, Fu, Yan, Hu, Xu and Khaitovich2010), the mismatch of age, sex, and ethnicity between expression and imaging data may introduce bias to this study. Another limitation is the impact of the selection of dynamic parameters on the calculation of dReHo. Although we used 25 combinations of dynamic parameters to identify consistent results, the development of strategies to select the optimal combination of parameters is still critically important.

In this study, we found that first-episode, treatment-naïve patients with MDD showed reduced dReHo in the regulatory networks of cognition and emotion, indicating that functional abnormalities in these networks are important features of MDD. The case-control differences in dReHo showed replicable spatial correlations with the expression profiles of gene modules associated with various cell types, biological functions and pathways, indicating that MDD may develop via complex polygenetic and polypathway mechanisms.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0033291720003876.

Acknowledgements

We thank the Allen Institute for Brain Science founders and staff who supplied the brain expression data. This study was funded by the National Key Research and Development Program of China (2018YFC1314300), the National Natural Science Foundation of China (82030053 and 82072001), the Tianjin Key Technology R&D Program (17ZXMFSY00090), and the Natural Science Foundation of Tianjin Municipal Science and Technology Commission (18JCQNJC10900).

Author contributions

Kaizhong Xue and Sixiang Liang made substantial contributions to the conception and design of the work and to the acquisition, analysis, and interpretation of data, and they drafted the manuscript. Bingbing Yang, Dan Zhu and Yingying Xie helped draft the manuscript. Wen Qin and Feng Liu critically revised the manuscript for important intellectual content. Yong Zhang and Chunshui Yu approved the final version of the manuscript to be published and agreed to be accountable for all aspects of the work, ensuring that questions related to the accuracy or integrity of any part of the work were appropriately investigated and resolved. All authors read, critically revised, and approved the manuscript.

Conflict of interest

None.

Ethical standards

The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008.

Footnotes

*

These authors contributed equally to this work.

References

Arnatkevic Iūtė, A., Fulcher, B. D., & Fornito, A. (2019). A practical guide to linking brain-wide gene expression and neuroimaging data. Neuroimage, 189, 353367. doi: 10.1016/j.neuroimage.2019.01.011.CrossRefGoogle Scholar
Ashburner, J. (2007). A fast diffeomorphic image registration algorithm. Neuroimage, 38(1), 95113. doi: 10.1016/j.neuroimage.2007.07.007.CrossRefGoogle ScholarPubMed
Chao-Gan, Y., & Yu-Feng, Z. (2010). DPARSF: A MATLAB toolbox for “pipeline” data analysis of resting-state fMRI. Frontiers in Systems Neuroscience, 4, 13. doi: 10.3389/fnsys.2010.00013.Google ScholarPubMed
Chen, J. D., Liu, F., Xun, G. L., Chen, H. F., Hu, M. R., Guo, X. F., … Zhao, J. P. (2012). Early and late onset, first-episode, treatment-naive depression: Same clinical symptoms, different regional neural activities. Journal of Affective Disorders, 143(1–3), 5663. doi: 10.1016/j.jad.2012.05.025.CrossRefGoogle ScholarPubMed
Chumbley, J., Worsley, K., Flandin, G., & Friston, K. (2010). Topological FDR for neuroimaging. Neuroimage, 49(4), 30573064. doi: 10.1016/j.neuroimage.2009.10.090.CrossRefGoogle ScholarPubMed
Cuijpers, P., Beekman, A. T. F., & Reynolds, C. F., 3rd. (2012). Preventing depression: A global priority. JAMA, 307(10), 10331034. doi: 10.1001/jama.2012.271.CrossRefGoogle ScholarPubMed
Demirtaş, M., Tornador, C., Falcón, C., López-Solà, M., Hernández-Ribas, R., Pujol, J., … Deco, G. (2016). Dynamic functional connectivity reveals altered variability in functional connectivity among patients with major depressive disorder. Human Brain Mapping, 37(8), 29182930. doi: 10.1002/hbm.23215.CrossRefGoogle ScholarPubMed
de Vos, F., Koini, M., Schouten, T. M., Seiler, S., van der Grond, J., Lechner, A., … Rombouts, S. (2018). A comprehensive analysis of resting state fMRI measures to classify individual patients with Alzheimer's disease. Neuroimage, 167, 6272. doi: 10.1016/j.neuroimage.2017.11.025.CrossRefGoogle ScholarPubMed
Dougherty, J. D., Schmidt, E. F., Nakajima, M., & Heintz, N. (2010). Analytical approaches to RNA profiling data for the identification of genes enriched in specific cells. Nucleic Acids Research, 38(13), 42184230. doi: 10.1093/nar/gkq130.CrossRefGoogle ScholarPubMed
Felger, J. C., & Lotrich, F. E. (2013). Inflammatory cytokines in depression: Neurobiological mechanisms and therapeutic implications. Neuroscience, 246, 199229. doi: 10.1016/j.neuroscience.2013.04.060.CrossRefGoogle ScholarPubMed
Fitzgerald, P. B., Laird, A. R., Maller, J., & Daskalakis, Z. J. (2008). A meta-analytic study of changes in brain activation in depression. Human Brain Mapping, 29(6), 683695. doi: 10.1002/hbm.20426.CrossRefGoogle ScholarPubMed
Friston, K. J. (2011). Functional and effective connectivity: A review. Brain Connectivity, 1(1), 1336. doi: 10.1089/brain.2011.0008.CrossRefGoogle ScholarPubMed
Friston, K. J., Williams, S., Howard, R., Frackowiak, R. S., & Turner, R. (1996). Movement-related effects in fMRI time-series. Magnetic Resonance in Medicine, 35(3), 346355. doi: 10.1002/mrm.1910350312.CrossRefGoogle ScholarPubMed
Frith, U., & Frith, C. (2010). The social brain: Allowing humans to boldly go where no other species has been. Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences, 365(1537), 165176. doi: 10.1098/rstb.2009.0160.CrossRefGoogle ScholarPubMed
Geng, J., Yan, R., Shi, J., Chen, Y., Mo, Z., Shao, J., … Lu, Q. (2019). Altered regional homogeneity in patients with somatic depression: A resting-state fMRI study. Journal of Affective Disorders, 246, 498505. doi: 10.1016/j.jad.2018.12.066.CrossRefGoogle ScholarPubMed
Gold, P. W. (2015). The organization of the stress system and its dysregulation in depressive illness. Molecular Psychiatry, 20(1), 3247. doi: 10.1038/mp.2014.163.CrossRefGoogle ScholarPubMed
Gottesman, I. I., & Gould, T. D. (2003). The endophenotype concept in psychiatry: Etymology and strategic intentions. The American Journal of Psychiatry, 160(4), 636645. https://doi.org/10.1176/appi.ajp.160.4.636.CrossRefGoogle ScholarPubMed
Guo, W. B., Liu, F., Xue, Z. M., Yu, Y., Ma, C. Q., Tan, C. L., … Zhao, J. P. (2011a). Abnormal neural activities in first-episode, treatment-naïve, short-illness-duration, and treatment-response patients with major depressive disorder: A resting-state fMRI study. Journal of Affective Disorders, 135(1–3), 326331. doi: 10.1016/j.jad.2011.06.048.CrossRefGoogle Scholar
Guo, W. B., Sun, X. L., Liu, L., Xu, Q., Wu, R. R., Liu, Z. N., … Zhao, J. P. (2011b). Disrupted regional homogeneity in treatment-resistant depression: A resting-state fMRI study. Progress in Neuro-Psychopharmacology & Biological Psychiatry, 35(5), 12971302. doi: 10.1016/j.pnpbp.2011.02.006.CrossRefGoogle Scholar
Han, S., Wang, X., He, Z., Sheng, W., Zou, Q., Li, L., … Chen, H. (2019). Decreased static and increased dynamic global signal topography in major depressive disorder. Progress in Neuro-Psychopharmacology & Biological Psychiatry, 94, 109665. doi: 10.1016/j.pnpbp.2019.109665.CrossRefGoogle ScholarPubMed
Hawrylycz, M. J., Lein, E. S., Guillozet-Bongaarts, A. L., Shen, E. H., Ng, L., Miller, J. A., … Jones, A. R. (2012). An anatomically comprehensive atlas of the adult human brain transcriptome. Nature, 489(7416), 391399. doi: 10.1038/nature11405.CrossRefGoogle ScholarPubMed
Hawrylycz, M., Miller, J. A., Menon, V., Feng, D., Dolbeare, T., Guillozet-Bongaarts, A. L., … Lein, E. (2015). Canonical genetic signatures of the adult human brain. Nature Neuroscience, 18(12), 18321844. doi: 10.1038/nn.4171.CrossRefGoogle ScholarPubMed
Hickok, G., & Poeppel, D. (2007). The cortical organization of speech processing. Nature Reviews. Neuroscience, 8(5), 393402. doi: 10.1038/nrn2113.CrossRefGoogle ScholarPubMed
Hodes, G. E., Kana, V., Menard, C., Merad, M., & Russo, S. J. (2015). Neuroimmune mechanisms of depression. Nature Neuroscience, 18(10), 13861393. doi: 10.1038/nn.4113.CrossRefGoogle ScholarPubMed
Howard, D. M., Adams, M. J., Clarke, T.-K., Hafferty, J. D., Gibson, J., Shirali, M., … McIntosh, A. M. (2019). Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions. Nature Neuroscience, 22(3), 343352. doi: 10.1038/s41593-018-0326-7.CrossRefGoogle ScholarPubMed
Howard, D. M., Adams, M. J., Shirali, M., Clarke, T. K., Marioni, R. E., Davies, G., … McIntosh, A. M. (2018). Genome-wide association study of depression phenotypes in UK biobank identifies variants in excitatory synaptic pathways. Nature Communications, 9(1), 1470. doi: 10.1038/s41467-018-03819-3.CrossRefGoogle ScholarPubMed
Hutchison, R. M., & Morton, J. B. (2015). Tracking the brain's functional coupling dynamics over development. The Journal of Neuroscience: The Official Journal of the Society for Neuroscience, 35(17), 68496859. doi: 10.1523/jneurosci.4638-14.2015.CrossRefGoogle ScholarPubMed
Hutchison, R. M., Womelsdorf, T., Allen, E. A., Bandettini, P. A., Calhoun, V. D., Corbetta, M., … Chang, C. (2013). Dynamic functional connectivity: Promise, issues, and interpretations. Neuroimage, 80, 360378. doi: 10.1016/j.neuroimage.2013.05.079.CrossRefGoogle ScholarPubMed
Iadecola, C. (2017). The neurovascular unit coming of age: A journey through neurovascular coupling in health and disease. Neuron, 96(1), 1742. doi: 10.1016/j.neuron.2017.07.030.CrossRefGoogle ScholarPubMed
Kaiser, R. H., Whitfield-Gabrieli, S., Dillon, D. G., Goer, F., Beltzer, M., Minkel, J., … Pizzagalli, D. A. (2016). Dynamic resting-state functional connectivity in major depression. Neuropsychopharmacology, 41(7), 18221830. doi: 10.1038/npp.2015.352.CrossRefGoogle ScholarPubMed
Kanwisher, N., McDermott, J., & Chun, M. M. (1997). The fusiform face area: A module in human extrastriate cortex specialized for face perception. The Journal of Neuroscience: The Official Journal of the Society for Neuroscience, 17(11), 43024311. doi: 10.1523/JNEUROSCI.17-11-04302.1997.CrossRefGoogle ScholarPubMed
Kessler, R. C. (1997). The effects of stressful life events on depression. Annual Review of Psychology, 48, 191214. doi: 10.1146/annurev.psych.48.1.191.CrossRefGoogle ScholarPubMed
Khan, A., Bhat, A., Kolts, R., Thase, M. E., & Brown, W. (2010). Why has the antidepressant-placebo difference in antidepressant clinical trials diminished over the past three decades? CNS Neuroscience & Therapeutics, 16(4), 217226. doi: 10.1111/j.1755-5949.2010.00151.x.CrossRefGoogle ScholarPubMed
Krishnan, V., & Nestler, E. J. (2008). The molecular neurobiology of depression. Nature, 455(7215), 894902. doi: 10.1038/nature07455.CrossRefGoogle ScholarPubMed
Lai, C.-H. (2018). The regional homogeneity of cingulate-precuneus regions: The putative biomarker for depression and anxiety. Journal of Affective Disorders, 229, 171176. doi: 10.1016/j.jad.2017.12.086.CrossRefGoogle ScholarPubMed
Leonardi, N., & Van De Ville, D. (2015). On spurious and real fluctuations of dynamic functional connectivity during rest. Neuroimage, 104, 430436. doi: 10.1016/j.neuroimage.2014.09.007.CrossRefGoogle ScholarPubMed
Li, J., Duan, X., Cui, Q., Chen, H., & Liao, W. (2019). More than just statics: Temporal dynamics of intrinsic brain activity predicts the suicidal ideation in depressed patients. Psychological Medicine, 49(5), 852860. doi: 10.1017/s0033291718001502.CrossRefGoogle ScholarPubMed
Li, L., Lu, B., & Yan, C. G. (2020). Stability of dynamic functional architecture differs between brain networks and states. Neuroimage, 216, 116230. doi: 10.1016/j.neuroimage.2019.116230.CrossRefGoogle ScholarPubMed
Liao, W., Zhang, Z., Mantini, D., Xu, Q., Ji, G. J., Zhang, H., … Lu, G. (2014). Dynamical intrinsic functional architecture of the brain during absence seizures. Brain Structure & Function, 219(6), 20012015. doi: 10.1007/s00429-013-0619-2.CrossRefGoogle ScholarPubMed
Liu, F., Guo, W., Liu, L., Long, Z., Ma, C., Xue, Z., … Chen, H. (2013). Abnormal amplitude low-frequency oscillations in medication-naive, first-episode patients with major depressive disorder: A resting-state fMRI study. Journal of Affective Disorders, 146(3), 401406. doi: 10.1016/j.jad.2012.10.001.CrossRefGoogle ScholarPubMed
Liu, F., Tian, H., Li, J., Li, S., & Zhuo, C. (2019). Altered voxel-wise gray matter structural brain networks in schizophrenia: Association with brain genetic expression pattern. Brain Imaging and Behavior, 13(2), 493502. doi: 10.1007/s11682-018-9880-6.CrossRefGoogle ScholarPubMed
Liu, F., Wang, Y., Li, M., Wang, W., Li, R., Zhang, Z., … Chen, H. (2017). Dynamic functional network connectivity in idiopathic generalized epilepsy with generalized tonic-clonic seizure. Human Brain Mapping, 38(2), 957973. doi: 10.1002/hbm.23430.CrossRefGoogle ScholarPubMed
Liu, Z., Xu, C., Xu, Y., Wang, Y., Zhao, B., Lv, Y., … Du, C. (2010). Decreased regional homogeneity in insula and cerebellum: A resting-state fMRI study in patients with major depression and subjects at high risk for major depression. Psychiatry Research, 182(3), 211215. doi: 10.1016/j.pscychresns.2010.03.004.CrossRefGoogle ScholarPubMed
Ma, C., Ding, J., Li, J., Guo, W., Long, Z., Liu, F., … Chen, H. (2012). Resting-state functional connectivity bias of middle temporal gyrus and caudate with altered gray matter volume in major depression. PLoS ONE, 7(9), e45263. doi: 10.1371/journal.pone.0045263.CrossRefGoogle ScholarPubMed
McEwen, B. S. (2008). Understanding the potency of stressful early life experiences on brain and body function. Metabolism: Clinical and Experimental, 57(Suppl 2), S1115. doi: 10.1016/j.metabol.2008.07.006.CrossRefGoogle ScholarPubMed
Mogil, L. S., Andaleon, A., Badalamenti, A., Dickinson, S. P., Guo, X., Rotter, J. I., … Wheeler, H. E. (2018). Genetic architecture of gene expression traits across diverse populations. PLoS Genetics, 14(8), e1007586. doi: 10.1371/journal.pgen.1007586.CrossRefGoogle ScholarPubMed
Morgan, S. E., Seidlitz, J., Whitaker, K. J., Romero-Garcia, R., Clifton, N. E., Scarpazza, C., … Bullmore, E. T. (2019). Cortical patterning of abnormal morphometric similarity in psychosis is associated with brain expression of schizophrenia-related genes. Proceedings of the National Academy of Sciences of the United States of America, 116(19), 96049609. doi: 10.1073/pnas.1820754116CrossRefGoogle ScholarPubMed
Munoz-Lopez, M. M., Mohedano-Moriano, A., & Insausti, R. (2010). Anatomical pathways for auditory memory in primates. Frontiers in Neuroanatomy, 4, 129. doi: 10.3389/fnana.2010.00129.CrossRefGoogle ScholarPubMed
Muoio, V., Persson, P. B., & Sendeski, M. M. (2014). The neurovascular unit – concept review. Acta Physiologica (Oxford, England), 210(4), 790798. doi: 10.1111/apha.12250.CrossRefGoogle ScholarPubMed
Olsson, A., & Ochsner, K. N. (2008). The role of social cognition in emotion. Trends in Cognitive Scences, 12(2), 6571. doi: 10.1016/j.tics.2007.11.010.CrossRefGoogle ScholarPubMed
Pittenger, C., & Duman, R. S. (2008). Stress, depression, and neuroplasticity: A convergence of mechanisms. Neuropsychopharmacology: Official Publication of the American College of Neuropsychopharmacology, 33(1), 88109. doi: 10.1038/sj.npp.1301574.CrossRefGoogle Scholar
Qiu, L., Xia, M., Cheng, B., Yuan, L., Kuang, W., Bi, F., … Gong, Q. (2018). Abnormal dynamic functional connectivity of amygdalar subregions in untreated patients with first-episode major depressive disorder. Journal of Psychiatry & Neuroscience: JPN, 43(4), 262272. doi: 10.1503/jpn.170112.CrossRefGoogle ScholarPubMed
Reimand, J., Arak, T., Adler, P., Kolberg, L., Reisberg, S., Peterson, H., & Vilo, J. (2016). g:Profiler – a web server for functional interpretation of gene lists (2016 update). Nucleic Acids Research, 44(W1), W83W89. doi: 10.1093/nar/gkw199.CrossRefGoogle Scholar
Romero-Garcia, R., Warrier, V., Bullmore, E. T., Baron-Cohen, S., & Bethlehem, R. A. I. (2019). Synaptic and transcriptionally downregulated genes are associated with cortical thickness differences in autism. Molecular Psychiatry, 24(7), 10531064. doi: 10.1038/s41380-018-0023-7.CrossRefGoogle ScholarPubMed
Romme, I. A., de Reus, M. A., Ophoff, R. A., Kahn, R. S., & van den Heuvel, M. P. (2017). Connectome disconnectivity and cortical gene expression in patients with schizophrenia. Biological Psychiatry, 81(6), 495502. doi: 10.1016/j.biopsych.2016.07.012.CrossRefGoogle ScholarPubMed
Sakoglu, U., Pearlson, G. D., Kiehl, K. A., Wang, Y. M., Michael, A. M., & Calhoun, V. D. (2010). A method for evaluating dynamic functional network connectivity and task-modulation: Application to schizophrenia. Magma (New York. N.Y.), 23(5–6), 351366. doi: 10.1007/s10334-010-0197-8.Google Scholar
Schmaal, L., Hibar, D. P., Samann, P. G., Hall, G. B., Baune, B. T., Jahanshad, N., … Veltman, D. J. (2017). Cortical abnormalities in adults and adolescents with major depression based on brain scans from 20 cohorts worldwide in the ENIGMA Major Depressive Disorder Working Group. Molecular Psychiatry, 22(6), 900909. doi: 10.1038/mp.2016.60.CrossRefGoogle ScholarPubMed
Setiawan, E., Wilson, A. A., Mizrahi, R., Rusjan, P. M., Miler, L., Rajkowska, G., … Meyer, J. H. (2015). Role of translocator protein density, a marker of neuroinflammation, in the brain during major depressive episodes. JAMA Psychiatry, 72(3), 268275. doi: 10.1001/jamapsychiatry.2014.2427.CrossRefGoogle ScholarPubMed
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., … Ideker, T. (2003). Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Research, 13(11), 24982504. doi: 10.1101/gr.1239303.CrossRefGoogle ScholarPubMed
Somel, M., Guo, S., Fu, N., Yan, Z., Hu, H. Y., Xu, Y., … Khaitovich, P. (2010). MicroRNA, mRNA, and protein expression link development and aging in human and macaque brain. Genome Research, 20(9), 12071218. doi: 10.1101/gr.106849.110.CrossRefGoogle ScholarPubMed
Sullivan, P. F., de Geus, E. J., Willemsen, G., James, M. R., Smit, J. H., Zandbelt, T., … Penninx, B. W. (2009). Genome-wide association for major depressive disorder: A possible role for the presynaptic protein piccolo. Molecular Psychiatry, 14(4), 359375. doi: 10.1038/mp.2008.125.CrossRefGoogle ScholarPubMed
Sullivan, P. F., Neale, M. C., & Kendler, K. S. (2000). Genetic epidemiology of major depression: Review and meta-analysis. The American Journal of Psychiatry, 157(10), 15521562. doi: 10.1176/appi.ajp.157.10.1552.CrossRefGoogle ScholarPubMed
Sun, H., Luo, L., Yuan, X., Zhang, L., He, Y., Yao, S., … Xiao, J. (2018). Regional homogeneity and functional connectivity patterns in major depressive disorder, cognitive vulnerability to depression and healthy subjects. Journal of Affective Disorders, 235, 229235. doi: 10.1016/j.jad.2018.04.061.CrossRefGoogle ScholarPubMed
Truong, W., Minuzzi, L., Soares, C. N., Frey, B. N., Evans, A. C., MacQueen, G. M., & Hall, G. B. (2013). Changes in cortical thickness across the lifespan in major depressive disorder. Psychiatry Research, 214(3), 204211. doi: 10.1016/j.pscychresns.2013.09.003.CrossRefGoogle ScholarPubMed
Van Dijk, K. R., Sabuncu, M. R., & Buckner, R. L. (2012). The influence of head motion on intrinsic functional connectivity MRI. Neuroimage, 59(1), 431438. doi: 10.1016/j.neuroimage.2011.07.044.CrossRefGoogle ScholarPubMed
Vasic, N., Walter, H., Hose, A., & Wolf, R. C. (2008). Gray matter reduction associated with psychopathology and cognitive dysfunction in unipolar depression: A voxel-based morphometry study. Journal of Affective Disorders, 109(1–2), 107116. doi: 10.1016/j.jad.2007.11.011.CrossRefGoogle ScholarPubMed
Weiner, K. S., & Grill-Spector, K. (2010). Sparsely-distributed organization of face and limb activations in human ventral temporal cortex. Neuroimage, 52(4), 15591573. doi: 10.1016/j.neuroimage.2010.04.262.CrossRefGoogle ScholarPubMed
Willner, P., Scheel-Kruger, J., & Belzung, C. (2013). The neurobiology of depression and antidepressant action. Neuroscience and Biobehavioral Reviews, 37(10 Pt 1), 23312371. doi: 10.1016/j.neubiorev.2012.12.007CrossRefGoogle ScholarPubMed
World Health Organization (2008). The Global Burden of Disease: 2004 Update. 2004 Update. Geneva, Switzerland: World Health Organization, p. 146.Google Scholar
Wray, N. R., Ripke, S., Mattheisen, M., Trzaskowski, M., Byrne, E. M., & Abdellaoui, A., … Major Depressive Disorder Working Group of the Psychiatric Genomics, C. (2018). Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nature Genetics, 50(5), 668681. doi: 10.1038/s41588-018-0090-3.CrossRefGoogle ScholarPubMed
Xia, M., Si, T., Sun, X., Ma, Q., Liu, B., Wang, L., … Group, D. I.-M. D. D. W. (2019). Reproducibility of functional brain alterations in major depressive disorder: Evidence from a multisite resting-state functional MRI study with 1434 individuals. Neuroimage, 189, 700714. doi: 10.1016/j.neuroimage.2019.01.074.CrossRefGoogle ScholarPubMed
Xu, J., Li, Q., Qin, W., Jun Li, M., Zhuo, C., Liu, H., … Yu, C. (2018). Neurobiological substrates underlying the effect of genomic risk for depression on the conversion of amnestic mild cognitive impairment. Brain: A Journal of Neurology, 141(12), 34573471. doi: 10.1093/brain/awy277.CrossRefGoogle ScholarPubMed
Yan, C. G., Wang, X. D., Zuo, X. N., & Zang, Y. F. (2016). DPABI: Data processing & analysis for (resting-state) brain imaging. Neuroinformatics, 14(3), 339351. doi: 10.1007/s12021-016-9299-4.CrossRefGoogle ScholarPubMed
Zang, Y., Jiang, T., Lu, Y., He, Y., & Tian, L. (2004). Regional homogeneity approach to fMRI data analysis. Neuroimage, 22(1), 394400. doi: 10.1016/j.neuroimage.2003.12.030.CrossRefGoogle ScholarPubMed
Zhang, B., & Horvath, S. (2005). A general framework for weighted gene co-expression network analysis. Statistical Applications in Genetics and Molecular Biology, 4, Article17. doi: 10.2202/1544-6115.1128CrossRefGoogle ScholarPubMed
Zhang, Y., Sloan, S. A., Clarke, L. E., Caneda, C., Plaza, C. A., Blumenthal, P. D., … Barres, B. A. (2016). Purification and characterization of progenitor and mature human astrocytes reveals transcriptional and functional differences with mouse. Neuron, 89(1), 3753. doi: 10.1016/j.neuron.2015.11.013.CrossRefGoogle ScholarPubMed
Figure 0

Fig. 1. Overview of the analysis pipeline. (a) a sliding window method was utilized to compute dReHo, and between-group comparisons of dReHo maps were performed; (b) the association between dReHo differences and gene expression from the AHBA was explored; (c) stable cortical and subcortical modules related to MDD were identified; (d) enrichment analysis was applied for these modules. dReHo, dynamic regional homogeneity; AHBA, Allen Human Brain Atlas.

Figure 1

Fig. 2. A probability map showing brain voxels with significant dReHo differences (q < 0.05, topological FDR corrected) between the MDD and HC groups. The colour bar represents the probability of a given voxel showing a significant intergroup difference in the 25 comparisons. Each comparison is performed on individual dReHo maps calculated by a distinct combination of dynamic parameters. FDR, false discovery rate; HC, healthy control; L, left; MDD, major depressive disorder; R, right.

Figure 2

Fig. 3. Consensus analysis of expression data. (a) WGCNA dendrogram showing consensus modules based on the topological co-expression of 10 185 genes from four left hemispheres and two whole brains. (b) To avoid the scenario that co-expression networks constructed accidentally based on the set of WGCNA parameters used, we rebuilt networks by applying different parameter sets. The final consensus modules used in this study were conserved among the majority of the tests with different parameter sets. WGCNA, weighted gene co-expression network analysis.

Figure 3

Fig. 4. Correlations between gene expression and dReHo changes in MDD. (a) The heatmap shows Pearson's correlation coefficients between MEs of gene modules and t values of case-control dReHo differences in cortical samples. Each column represents a gene module (n = 18), and each row represents a case-control dReHo difference map derived from dReHo maps calculated with a given combination (n = 25) of window size and step size. In the rows of the upper half, gene expression data were obtained from two AHBA donors with whole-brain coverage; in the rows of the lower half, however, they were obtained from four AHBA donors only with left-brain coverage. (b) The heatmap shows correlations in subcortical samples. * indicates that the correlation reaches the significance of q < 0.05 (BH-FDR corrected) in all analyses. The colour bar represents Pearson's correlation coefficients between MEs of gene modules and t values of case-control dReHo differences.

Figure 4

Fig. 5. Enrichment of genes in stable modules and hub genes. (a) Counts of MDD-related cortical and subcortical modules. (b) KEGG enrichment for genes in all modules (top two pathways shown for each module). All obtained results could pass q < 0.05 (BH-FDR correction) and showed odds ratios >1. (c) GO enrichment for genes in all modules (top two pathways shown for each module marked by different colours). (d) Cell-type marker enrichment of genes in cortical modules. The bubble size represents significance for – log10FDR (BH-FDR correction), and the bubble colour represents whether the module is significant or not. (e) The top 10 hub genes are shown for all MDD-related modules. Edges are weighted by the strength of correlation between nodes (genes). KEGG, Kyoto Encyclopedia of Genes and Genomes; GO, Gene Ontology.

Supplementary material: File

Xue et al. Supplementary Materials

Xue et al. Supplementary Materials

Download Xue et al. Supplementary Materials(File)
File 41.6 MB