Shaoqiang Han1,2,3,4,5,6,7*
Ruiping Zheng1,2,3,4,5,6,7
Shuying Li8Bingqian Zhou1,2,3,4,5,6,7Yu Jiang1,2,3,4,5,6,7
Caihong Wang1,2,3,4,5,6,7
Yarui Wei1,2,3,4,5,6,7Jianyue Pang8Hengfen Li8
Yong Zhang1,2,3,4,5,6,7Yuan Chen1,2,3,4,5,6,7*
Jingliang Cheng1,2,3,4,5,6,7*
- 1Department of Magnetic Resonance Imaging, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
- 2Key Laboratory for Functional Magnetic Resonance Imaging and Molecular Imaging of Henan Province, Zhengzhou, China
- 3Engineering Technology Research Center for Detection and Application of Brain Function of Henan Province, Zhengzhou, China
- 4Engineering Research Center of Medical Imaging Intelligent Diagnosis and Treatment of Henan Province, Zhengzhou, China
- 5Key Laboratory of Magnetic Resonance and Brain Function of Henan Province, Zhengzhou, China
- 6Key Laboratory of Brain Function and Cognitive Magnetic Resonance Imaging of Zhengzhou, Zhengzhou, China
- 7Key Laboratory of Imaging Intelligence Research Medicine of Henan Province, Zhengzhou, China
- 8Department of Psychiatry, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, China
The pathophysiology and pharmacology of depression are hypothesized to be related to the imbalance of excitation–inhibition that gives rise to hierarchical dynamics (or intrinsic timescale gradient), further supporting a hierarchy of cortical functions. On this assumption, intrinsic timescale gradient is theoretically altered in depression. However, it remains unknown. We investigated altered intrinsic timescale gradient recently developed to measure hierarchical brain dynamics gradient and its underlying molecular architecture and brain-wide gene expression in depression. We first presented replicable intrinsic timescale gradient in two independent Chinese Han datasets and then investigated altered intrinsic timescale gradient and its possible underlying molecular and transcriptional bases in patients with depression. As a result, patients with depression showed stage-specifically shorter timescales compared with healthy controls according to illness duration. The shorter timescales were spatially correlated with monoamine receptor/transporter densities, suggesting the underlying molecular basis of timescale aberrance and providing clues to treatment. In addition, we identified that timescale aberrance-related genes ontologically enriched for synapse-related and neurotransmitter (receptor) terms, elaborating the underlying transcriptional basis of timescale aberrance. These findings revealed atypical timescale gradient in depression and built a link between neuroimaging, transcriptome, and neurotransmitter information, facilitating an integrative understanding of depression.
Introduction
As one of the leading disabling diseases worldwide (Murray et al., 2012), depression affects approximately 350 million people each year (Schmaal et al., 2017). Recent studies point out that imbalance of the excitation–inhibition (E/I) underlays the pathophysiology and pharmacology of the depression (Voineskos et al., 2019). Imbalance of E/I ratio hypothetically results in the aberrance of hierarchically organized intrinsic neural timescales (Kiebel et al., 2008) that support synchronizing large-scale brain networks usually measured with resting-state functional connectivity (rsFC; Buzsáki and Draguhn, 2004). Accordingly, the intrinsic timescale gradient is theoretically altered in depression that remains unknown yet.
The brain regions are hierarchically organized into increasing polyfunctional neural circuits embodied in topographic gradients of molecular, cellular, and anatomical properties (Huntenburg et al., 2018). Emerging through hierarchically organized feature (Burt et al., 2018) such as pyramidal cell dendritic spine density (Elston, 2003), long-range interactions (Wang, 2020), and gene expression gradients (Fulcher et al., 2019), intrinsic neural dynamics (or intrinsic timescale gradient) are also hierarchically organized, supporting a hierarchy of cortical functions (Kiebel et al., 2008; Hasson et al., 2015). Brain dynamics is also hierarchically organized along spatial gradients extending from sensorimotor regions to association cortex (Hasson et al., 2008) supporting functional communications between brain regions (Cocchi et al., 2016). Regions with longer “temporal receptive windows” are subsequently found to exhibit more slowly changing activity and vice versa (Hasson et al., 2015). In particular, regions such as prefrontal areas and parietal areas, densely interconnected central regions, have longer timescales compared to peripheral sensory areas (Chaudhuri et al., 2015) for the reason that prolonged neural timescale is needed to enable these high-order brain regions to integrate various information for robust sensory perception (Hasson et al., 2008), stable memory processing (Bernacchia et al., 2011), and decision-making (Cavanagh et al., 2016). By developing a large-scale biophysical model, Chaudhuri et al. (2015) elaborate that intrinsic timescale gradient depends crucially on recurrent network activity. Aberrance of neural timescales is supposed to the result of imbalance of the excitation–inhibition (E/I) ratio (Wengler et al., 2020). Wengler et al. (2020) find evidence for distinct hierarchical aberrance in timescale gradient as a function of hallucination and delusion, supporting glutamatergic and dopamine theories of psychosis (Corlett et al., 2011; Jardri et al., 2016). Imbalance of the E/I is also implicated in the pathophysiology and pharmacology of the depression (Covington et al., 2010; Voineskos et al., 2019) and the mechanism of fast-acting antidepressant is related to E/I rebalance (Li, 2020). Although the intrinsic timescale gradient should be altered in depression theoretically, it remains unknown yet.
Brain function such as rsFC is also modulated by genetic factors (Richiardi et al., 2015; Fornito et al., 2019; Richiardi et al., 2015) and coupled to neurotransmitters (Stagg et al., 2014; Kringelbach and Cruzat, 2020). Twin studies show that functional connectivity within the default-mode network and topological measures in the human brain are moderate to highly heritable (Glahn et al., 2010; Fornito et al., 2011; van den Heuvel et al., 2013). Recently, Allen Human Brain Atlas (AHBA) (Hawrylycz et al., 2012), a newly proposed brain-wide gene expression atlas, provides the possibility of bridging the gap between transcriptome and large-scale connectome organization (Fornito et al., 2019). Following the work of Richiardi et al. (2015) where they find that the transcriptome profile similarity within networks is higher than that between networks, a number of studies begin to explore the transcriptional basis of macroscopic neuroimaging phenotypes (Krienen et al., 2016; Vértes et al., 2016; Anderson et al., 2018; Li and Seidlitz, 2021). Recently, Zhu et al. (2021) find that spatial distribution of functional connectivity strength is modulated by genes enriched for terms such as synaptic transmission in health (Zhang et al., 2021). In depression, Li et al. (2021) identify that altered morphometric similarity network is correlated with transcriptional signatures (Li and Seidlitz, 2021). In addition, rsFC is found to be coupled to neurotransmitter transporters/receptors (Stagg et al., 2014; Kringelbach and Cruzat, 2020). Dysconnectivity in schizophrenia is linked to altered neurotransmission (Landek-Salgado et al., 2016; Limongi et al., 2020). Chen et al. (2021) find that abnormal functional topography of brain networks is associated with the dopaminergic and serotonergic systems underlying cognitive decline in schizophrenia by investigating the molecular architecture facilitating a link to treatment. The variation of timescales is hypothesized to arise from local biophysical properties of neurons across the cortical hierarchy, such as the density of glutamate receptors, calcium channels, and regulators of synaptic depression and facilitation (Zucker and Regehr, 2002; Wong and Wang, 2006). Investigating the molecular and transcriptional basis of altered intrinsic timescale gradient in depression helps to advance our understanding of how alterations at microscale architecture drive macroscale neuroimaging aberrance in depression.
In this study, we aimed to explore altered intrinsic timescale gradient and its underlying molecular and transcriptional signatures bridging the gap between molecular mechanism and macroscopic neuroimaging phenotypes in depression. First, we presented replicable landscape of intrinsic timescale gradient and its association with commonly used functional indicators including amplitude of low-frequency fluctuation (ALFF) and functional connectivity density (FCD) in two independent Chinese Han cohorts. Second, we investigated altered intrinsic timescale gradient in different stages of depression according to illness duration. Third, we inquired molecular and transcriptional basis of altered intrinsic timescale gradient in depression. Fourth, a functional enrichment analysis was performed to inquire ontological pathways of timescale aberrance-related genes in depression.
Materials and Methods
Datasets
Two independent Chinese Han datasets were used in this study. The first dataset come from the Southwest University Adult Lifespan Dataset (SALD) study including 494 healthy participants (female: male, 308:187, 19--80 years old). The second dataset included 121 HCs and 191 patients with depression. The resting-state functional MRI data were acquired and preprocessed using Data Processing & Analysis for Brain Imaging (DPABI)1 (Yan et al., 2016). The details about dataset description, scan acquisition, and preprocessing procedures were included in Supplementary Methods. The study was approved by the research ethical committee of the First Affiliated Hospital of Zhengzhou University.
Calculation of Timescales and Its Association With Functional Connectivity Density and Amplitude of Low-Frequency Fluctuation
Based on previous studies (Watanabe et al., 2019; Raut and Snyder, 2020), we calculated intrinsic neural timescales by calculating the magnitude of autocorrelation of the resting-state brain signals. There were two different definitions of timescales. First, the timescale was defined as the sum of positive autocorrelation function (ACF) values and then multiplied by the repetition time (TR) (Watanabe et al., 2019). The results reported in the next steps were based on this definition. The second was defined as the half of the full width at half maximum of the ACF (Raut and Snyder, 2020). To inquire the relationship between these two definitions, spatial correlation was obtained between the two mean timescale maps across healthy subjects in dataset 1 and healthy subjects in dataset 2, respectively. In addition, we also compared spatial correlation between the altered timescale maps of depression using these two definitions (see below).
As a newly proposed index, factors affecting landscape of timescales remained unclear. To explore these factors, we investigated whether factors such as gender (female vs. male), age, motion movement, and education level could affect intrinsic timescale gradient (Supplementary Methods). The timescales were calculated using custom MATLAB code publicly available at https://github.com/RaichleLab.
To intuitively elucidate the intrinsic timescale gradient measured, we inquired the association between intrinsic timescale gradient with other common functional indexes such as ALFF and FCD. For FCD maps, local, long-range, and global FCD maps were calculated (Tomasi and Volkow, 2010) where correlation threshold was determined by significance of single functional connection (p-value). A functional connection (correlation coefficient) was considered significant if its p < 0.05 (Bonferroni corrected). The obtained FCD maps were transformed to z-scores by subtracting the mean value and dividing by the standard deviation across gray matter voxels. The ALFF maps were calculated using resting-state functional magnetic resonance imaging data processing toolbox (REST) (Song et al., 2011), and normalized ALFF maps (dividing the mean value across gray matter voxels) were chosen for the following steps.
Dominance analysis was used to quantify the association between mean FCD (including local, long-range and global FCD), and mean ALFF maps to landscape of timescales across healthy subjects in each dataset (Budescu and David, 1993; Azen and Budescu, 2003).2 The details are provided in Supplementary Methods.
All of the above analysis steps were done in the discovery cohort (dataset 1) and validated in the replication cohort (HCs of dataset 2); results reported were based on the discovery cohort, unless stated otherwise.
Altered Intrinsic Timescale Gradient in Depression
Then, we explored whether intrinsic timescale gradient was altered in depression. The altered timescale was obtained by using two-tailed two-sample t-test equipped in SPM 12 where gender, age, mean FD, SNR0, and educational level were included as covariates. To explore whether aberrance of intrinsic timescale gradient was stage-dependent for the reason that mental disorders were found to present progressive brain structural alterations (Koutsouleris et al., 2014; Treadway et al., 2015; Cao et al., 2017; Zhang et al., 2017; Yüksel et al., 2018), patients with depression were further divided into three stages according to illness duration (Stage 1: 0 ≤ illness duration ≤12 months; Stage 2: 12 <illness duration ≤24 months; Stage 3: illness duration ≥24 months). Moreover, we also compared timescales in patients whose illness duration is less than 3/6 months to inquire whether timescale aberrance emerged from the beginning of the disease. Results reported in this study were corrected for multiple comparison (voxel-wise p < 0.001, cluster-level p < 0.05; GRF correction). To further explore the relationship between two definitions of timescales, the same statistical procedures were done in the second definition of timescales (half maximum of the ACF).
Spatial Correlation Between Altered Timescales of Depression With Receptor/Transporter Densities
To explore association between depression-induced changes in timescales and expression of a specific receptor/transporter, we evaluated the spatial relationship between altered timescales and the distribution of receptors/transporters. The timescale difference map was spatially correlated with PET/SPECT maps in JuSpace toolbox3 (Dukart et al., 2021). The default neuromorphometrics atlas excluding white matter and cerebrospinal fluid regions was used. Dopamine (D1 and D2), serotonin receptors (5-HT1a, 5-HT1b, and 5-HT2a), transporters (dopamine transporter and serotonin reuptake transporter 5-HTT), F-DOPA (a reflection of presynaptic dopamine synthesis capacity), the GABAergic receptor, and the noradrenaline transporter (NAT) were investigated. The correlation results were adjusted for spatial autocorrelation of local gray matter probabilities, and the significance of results was computed using permutation statistics (Dukart et al., 2021). To exclude the effect of atlas choice on our results, we used another two atlas (268 and 246 regions) (Shen et al., 2013; Fan et al., 2016) to validate these results.
What is more, to inquire whether the association with receptor/transporter densities was specific to altered timescales, we calculated spatial correlation between ALFF differences of depression with PET/SPECT maps using default neuromorphometrics atlas. The ALFF was chosen for the reason that it was widely used in resting-state fMRI studies.
Cortical Gene Expression Related to Altered Timescale of Depression
Inspired by a previous study (Reardon and Seidlitz, 2018), we ranked genes based on the spatial correlation between gene expression pattern and the voxel-wise unthresholded t-statistic map of timescale difference in depression. The gene expression data come from theAHBA4 (Hawrylycz et al., 2012), obtained from six adult human brains (Hawrylycz et al., 2012). Details and preprocessing procedures of AHBA were included in Supplementary Methods. The preprocessed AHBA used in this study comes from the Brain Annotation Toolbox (BAT)5 (Liu et al., 2019). Because only two right hemisphere data were included in the AHBA, we only considered the left hemisphere in our analysis (Arnatkeviciute et al., 2019). The correlation results were considered significant if |r| > 0.2 and p < 0.05 (FWE corrected). Finally, the positive and negative correlation gene lists were identified with timescale aberrance.
Enrichment Pathways Associated With Altered Timescales of Depression
We performed the gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways with the genes presenting significant spatial correlation with altered timescales of depression using Metascape (Zhou et al., 2019). Results reported here were corrected by the FDR (p < 0.05). This procedure was done in positive and negative correlation genes separately.
As done in a previous study (Li and Seidlitz, 2021), we further investigated shared enrichment terms between previously reported polygenic risk for depression and the timescale-related gene list (Wray et al., 2018; Howard and Adams, 2019). A multi-gene list meta-analysis was carried out between the timescale aberrance-related gene list and the gene list provided by these two studies.
Results
Clinical Demographics
The clinical demographics of subjects in dataset 1 and dataset 2 are included in Supplementary Tables 1, 2.
The Landscape of Resting-State Timescales and Its Association With Functional Connectivity Density and Amplitude of Low-Frequency Fluctuation
The mean timescale maps across healthy subjects in dataset 1 and dataset 2 are drawn in Figure 1. In accordance with a previous study (Watanabe et al., 2019), both dataset 1 and dataset 2 (only HCs) presented similar whole-brain patterns of timescales with longer timescales in frontal and parietal cortices and shorter timescales in other regions such as sensorimotor and visual areas (Figure 1). The spatial correlation between dataset 1 and dataset 2 (only HCs) was r = 0.783 (p < 0.05, permutation test). As there were two definitions of timescale, we calculated the spatial correlation between mean maps of the two definitions in dataset 1 and dataset 2 (only HCs). The landscapes of these two definitions were in good agreement (dataset 1, r = 0.961, p < 0.05 for permutation test; dataset 2, r = 0.893, p < 0.05 for permutation test) (Supplementary Figure 2). We observed significantly negative correlation between age and timescales only in dataset 1 (Supplementary Figure 3), suggesting that the timescales might be related to normal brain aging; results of HCs in dataset 2 were not significant, possibly due to the limited sample size. In addition, timescales in regions such as the left inferior temporal gyrus, sensorimotor cortex, and left middle frontal gyrus presented significantly negative correlation with educational level stating its potential role in the landscape of educational level in HCs of dataset 2 (Supplementary Figure 4).
Figure 1. The landscape of timescales and its association with ALFF and FCD. The “r” meant the spatial correlation between mean timescales maps in dataset 1 and dataset 2. The number in the dominance results meant the percentage of ALFF/FCD contributing to timescales where higher number meant higher association with timescales.
Then, we assessed the association between intrinsic timescale gradient with common resting-state functional indexes including FCD and ALFF using dominance analysis. Results revealed the relative importance of each predictor (collective R2 = 0.4528, long-range FCD = 29.19%, global FCD = 27.64%, local FCD = 16.05%, ALFF = 27.11%), suggesting that long-range and global FCD contributed the most to intrinsic neural timescales. These results were validated in HCs of dataset 2 (Supplementary Table 3).
Altered Timescales in Depression
Overall, there was no significant aberrance of timescales in patients with depression across stages according to illness duration. Whereafter, we investigated altered timescales in patients at different stages according to illness duration. Timescales in patients with depression presented stage-specific aberrance. In particular, patients presented shorter timescales in regions including right anterior insula extending to right putamen only at the first stage (<12 months). To further explore whether this aberrance occurred at disease onset, we investigated timescale aberrance in patients with shorter illness duration (<3/6 months). Regions such as ventral medial prefrontal cortex vmPFC/subgenual PFC, dorsal ACC, dorsal lateral PFC, the bilateral nucleus accumbens (NAcc), the striatum, and the bilateral insula presented decreased timescales in patients with an illness duration of less than 3/6 months. With the prolongation of the disease course, the timescale alterations gradually faded away (Figure 2 and Supplementary Table 4). In addition, to further explore the relationship between two definitions of timescales, the same statistical procedures were done in the second definition of timescales (half maximum of the ACF), and the results are included in Supplementary Figure 5. These results confirmed good consistency of altered intrinsic timescale gradient with different definitions.
Figure 2. Stage-specific aberrance of timescales in depression.
To exclude the possibility that the gradual reduction of timescales resulted from the samples used in the current study, we also inquired whether another functional index also presented stage-specific aberrance. For the reason mentioned above, ALFF was chosen. As a result, patients with depression did not present gradual ALFF aberrance in depression (Supplementary Figure 6).
Relationship to Receptor/Transporter Densities
Altered timescales were significantly correlated (p < 0.05 for permutation, FWE corrected) with seven receptor/transporter densities (Figure 3 and Supplementary Table 5) including 5-HT2a (5-HT subtype 2a), D1 (dopamine D1), DAT (dopamine transporter), F-DOPA (dopamine synthesis capacity), NAT (noradrenaline transporter), and SERT (serotonin transporter). These results were validated with a different brain atlas (Figure 3). In addition, to explore whether the correlation results were specific to altered timescales, we also calculated spatial correlation between ALFF differences of depression with receptor/transporter densities using default neuromorphometrics atlas. As a result, there was no significant spatial correlation between ALFF aberrance and receptor/transporter densities, hinting that the association was specific to timescale aberrance in depression.
Figure 3. The association between timescale aberrance with receptor/transporter densities. The “*” represented that the correlation was significant (p < 0.05 for permutation, FWE corrected).
Note: 5-HT1a, 5-HT subtype 1a; 5-HT1b, 5-HT subtype 1b; 5-HTaa, 5-HT subtype aa; D1, dopamine D1; D2, dopamine D2; DAT, dopamine transporter; F-DOPA, dopamine synthesis capacity; NAT, noradrenaline transporter; SERT, serotonin transporter.
Cortical Gene Expression Related to Altered Timescales in Depression
As timescales presented stage-specific aberrance and the differences faded away as the progression of illness in depression, we calculated the spatial correlation between gene expression with timescale differences in patients with an illness duration of less than 3 months (see before). As a result, 865/264 genes presented positive/negative correlation with the unthresholded timescale differences in patients with depression whose illness duration was less than 3 months (see section “Supplementary Methods”).
Enrichment Pathways Associated With Altered Timescales of Depression
We performed the GO biological processes and KEGG pathways with the associated gene lists obtained in the previous step using Metascape. The top 30 significant GO biological processes, such as “Trans-synaptic signaling,” “synapse signaling,” “synapse organization,” “dendrite development,” and “cognition,” and one KEGG pathway, “cGMP-PKG signaling pathway,” were identified (Figure 4) for the positive correlation gene list. Regulation of neurotransmitter receptor activity such as “glutamatergic synapse,” “dopaminergic synapse,” “glutamate receptor signaling pathway,” and “regulation of neurotransmitter secretion” was also identified (see section “Supplementary Methods”). These enrichment terms were clustered into clusters such as synaptic signaling, synapse organization (e.g., dendritic spine morphogenesis and dendrite development), regulation of transmembrane transport, and head development (e.g., hippocampus development and limbic system development) (Figure 4).
Figure 4. Functional enrichment of gene transcripts. (A) Top 30 enrichment terms of positive correlation genes. The size of the circle represented the number of genes enriched in a given term. The color bar represented the significance of a given term. (B) Metascape enrichment network visualization. Each term was represented as a circle node where its size is proportional to the number of genes enriched in the term and its color represented cluster identify.
Then, we investigated shared enrichment terms between the previously reported polygenic risk for depression and the positive correlation gene list by performing a multi-gene list meta-analysis (Zhou et al., 2019). As a result, we found 11 common pathways. The enrichment pathways included “synaptic signaling,” “synapse organization,” “cell–cell adhesion via plasma membrane adhesion molecules,” and “dendrite development” (Supplementary Figure 7).
Discussion
In this study, we investigated altered intrinsic neural timescale gradient in patients with depression and its possible underlying molecular and transcriptional signatures. Timescales presented stage-specific aberrance in depression. Specifically, patients at the beginning of illness (illness duration <3 months) presented shorter timescales in regions including vmPFC, ACC, the bilateral nucleus accumbens (NAcc), the striatum, and the bilateral insula. As the illness advanced, the difference faded away (disappeared when illness duration ≥12 months). Moreover, the shorter timescales at the beginning of depression were associated with receptor/transporter densities including 5-HT2a, D1/2, DAT, F-DOPA, NAT, and SERT, suggesting the underlying molecular basis of timescale aberrance and providing clues to treatment. Then, we identified timescale aberrance-related genes ontologically enriched for synapse-related and neurotransmitter (receptor) terms elaborating the underlying transcriptional basis of timescale aberrance. These findings revealed atypical intrinsic timescale gradient in depression and bridged the gap between neuroimaging, transcriptome, and neurotransmitter information facilitating an integrative understanding of depression.
Stage-Specifically Shorter Intrinsic Timescales in Depression
Patients with depression presented stage-specifically shorter timescales according to illness duration in regions including vmPFC, ACC, the bilateral NAc, the striatum, and the bilateral insula. The shorter timescales were only observed in patients with illness duration less than 12 months and then faded away as illness advanced. Converging lines of evidence confirmed that depression was a neuroprogressive illness (Kendler et al., 2001; Moylan et al., 2013); the morphometric alteration of critical brain regions was related to illness progression information (such as illness duration) (Frodl et al., 2003; McKinnon et al., 2009; van Eijndhoven et al., 2009; Alexander-Bloch et al., 2013; Chen V. C. et al., 2016). Consistent with this notion, we observed stage-specifically shorter timescales in patients with illness duration less than 12 months, and even 3 months, suggesting that shorter timescales occurred at the beginning of the disease. In our previous study, we identified that higher brain age was also stage-dependent (Han et al., 2021). This stage-specific aberrance might explain inconsistent findings in depression (Chen Z. et al., 2016). Note that the insignificant timescale aberrance in patients with longer illness duration did not necessarily mean the remission of depression for the reason that we did observe a significant difference in the total score of HAMD (p = 0.139, F = 2.00) across stages. Regions presenting shorter timescales were found to be related to blunted processing of incentive salience, weak reward source memory, and reinforcement learning underlying the anhedonia in depression (Whitton et al., 2015; Alloy et al., 2016; Han et al., 2020). The shorter timescales of these regions might be associated with inefficient responsiveness to rewards in depression (Whitton et al., 2015; Alloy et al., 2016). On the other hand, dorsal lateral PFC, subgenual PFC, and dorsal ACC belonged to dorsal systems inhibiting amygdala activity in the unstressed state (Phillips et al., 2008). The reduced neuronal size and diminished dendritic arborization in the dorsal system were found in depression (Jaako-Movits et al., 2006; Drevets et al., 2008). The chronic stress could affect the gene expression of monoamine (serotonin)-glutamate/GABA and subsequently affected the E/I balance (Dygalo et al., 2020). Consistent with some ideas, we identified the shorter timescales in the dorsal system resulting from imbalance of the excitation–inhibition (E/I) ratio (Wengler et al., 2020), suggesting the inability of the dorsal system to regulate stress response in depression (Phillips et al., 2008).
There were two possible interpretations for the stage-dependent timescale alteration in depression. First, this phenomenon might mirror the transition from a clinically unstable period, with large variability in functioning, to a relatively stable period, when patients have reached a plateau in functioning (Davidson and McGlashan, 1997; van Haren et al., 2003). Another possible interpretation was that shorter timescales might be associated with primal brain dysfunction in depression. Here, we preferred the latter one. Multimodal lines of evidence convergently indicated that depression was a neuroprogressive illness (Kendler et al., 2001; Moylan et al., 2013). Even patients with depression suffering from only one depressive episode also displayed characteristics of a progressive illness (Moylan et al., 2013). As the disease prolonged, brain tissue damage and physiological functioning gradually changed, which underpinned symptomatology and functional decline over time (Moylan et al., 2013). Our previous results (under review) revealed that progressive morphological alteration might originate from regions like vmPFC and then expand to other regions in depression. Similar progressive morphological alterations such as advanced illness were observed in schizophrenia (Jiang et al., 2018), epilepsy (Zhang et al., 2017), and generalized anxiety disorder (Chen et al., 2020). The original dysfunction might be of great significance to the pathogeny and the treatment of depression. Actually, early treatment of patients with depression is usually accompanied with better outcome of antidepressant treatment and remission, and the reverse was also true. For example, a longer duration of untreated illness was reported to have an unfavorable effect on the subsequent course of the illness (e.g., higher number of recurrences) (Altamura et al., 2007, 2008; Li et al., 2021). A shorter duration of untreated illness was related to better remission of depression and somatic symptoms (Bukh et al., 2013). On the other hand, a longer duration of untreated illness was found to be associated with a greater severity and a lower improvement percentage (Hung et al., 2015; Kraus and Kadriu, 2019). What is more, we found that the stage-specific aberrance might be specific to an intrinsic timescale and did not result from the sample used in the current study. In summary, the shorter timescales at the beginning of disease might reflect initial functional aberrance and mean a lot to subsequent treatment.
Molecular Architecture of the Shorter Timescales
To explore the potential neurophysiological mechanism underlying the shorter timescales observed in depression helping to facilitate a link to treatment (Chen et al., 2021), we calculated spatial correlation between maps of a variety of neurotransmitter systems with that of timescale aberrance (Dukart et al., 2021). In line with the monoamine hypothesis (Liu et al., 2018), shorter timescales were associated with monoamine neurotransmitters including serotonin, noradrenaline, and dopamine at the same time. It was not unexpected that the timescale differences were associated with serotonin and noradrenaline neurotransmissions because of their fatal roles in pathogenesis (Hamon and Blier, 2013) and the first-line treatment of depression by inhibiting the action of the serotonin/noradrenaline transporter to reduce reuptake of serotonin/noradrenaline (Pirker et al., 1995). Consistent with studies showing that 5-HT2a and SERT were decreased in patients with depression (Kambeitz and Howes, 2015; Steinberg et al., 2019), we observed that the substrate of timescale differences might be related to 5-HT2a and SERT. The reason might be that the 5-HT2A receptors have both excitatory and inhibitory roles underlying the potential biological mechanism of timescale hierarchies (Chaudhuri et al., 2015). The association between SERT binding and rsFC (Beliveau et al., 2015) and dysfunction of SERT binding could result in altered functional connectivity in depression (Han et al., 2019) followed by altered timescales in depression. In addition, we found that the timescale aberrance might be also related to dopaminergic neurotransmission. Dopaminergic neurotransmission playing an essential role by rewarding prediction error (Hollerman and Schultz, 1998; Bayer and Glimcher, 2005) and mediating motivational drive by the attribution of incentive salience to reward-related stimuli (Berridge, 2007) was also related to anhedonia and amotivation in depression (Mayberg et al., 2005). Reduced DAT density in the central and basal nuclei of the amygdala was found in a post-mortem study (Klimek et al., 2002). The association between timescale differences with dopaminergic neurotransmission suggested that shorter timescales of these regions might result in inefficient responsiveness to rewards (Whitton et al., 2015; Alloy et al., 2016). Engaging additional targets (e.g., DA) could help patients with residual symptoms and treatment-resistant depression (Blier, 2016). Combining with these findings, our results revealed the role of dopaminergic neurotransmission in timescale aberrance. What is more, the validation results confirmed the robustness (selection of different atlas) and specificity (Compared with ALFF) of association between timescale aberrance and neurotransmitter information. These results suggested a potential neurophysiological mechanism underlying the shorter timescales observed in depression, providing clues to treatment.
Altered Timescale-Related Gene Expressions Enriched for Functional Annotations
We identified genes whose expression pattern presented significantly spatial (positive/negative) correlation with timescale aberrance and their ontology terms elaborated the underlying transcriptional basis of timescale aberrance. Consistent with the monoamine hypothesis in depression (Liu et al., 2018), positive correlation genes related to timescale difference were significantly enriched in monoamine neurotransmitter-related GO biological processes/KEGG pathways including neurotransmitter secretion, transport, and receptor activity/complex. These results corresponded with the aforesaid findings about monoamine aberrance, suggesting transcriptional mechanisms of association between timescale difference and monoamine neurotransmitters. In recent years, the synaptic dysfunction hypothesis that depression was caused by disruption of homeostatic mechanisms controlling synaptic plasticity (Duman and Aghajanian, 2012) has been proposed in consideration of the moderate and delayed effectiveness of the widely prescribed serotonin selective reuptake inhibitors (SSRIs) (Trivedi et al., 2006) and rapid antidepressant actions of ketamine in treatment-resistant depressed patients (Berman et al., 2000; Zarate et al., 2006). Deficits of excitatory glutamate neurons and inhibitory GABA interneurons resulted in the vulnerability of these major neurotransmitter systems followed by dendritic atrophy and spine loss in neurons of the hippocampus and prefrontal cortex (Qiao et al., 2016; Duman et al., 2019). Dendrite complexity and synaptic density can also be increased after treatment with antidepressants (Li et al., 2010; Li et al., 2011; MacQueen and Frodl, 2011). In our study, the timescale difference-related genes were enriched in terms of the charge of the balance of excitation and inhibition including glutamatergic synapse, transmission, receptor signaling pathway, GABAergic synapse, regulation of NMDA receptor activity, and G protein-coupled receptor signaling pathway. In fact, the spatial correlation (p = 0.024 uncorrected) between timescale aberrance and GABAa (gamma-aminobutyric acid) was also observed in the current study. These results elaborated possible transcriptional basis underlying the altered intrinsic timescale gradient in depression and provided new lines of evidence supporting the synaptic dysfunction hypothesis. In addition to the overlapping ontology terms with that in GWAS in depression, the multi-gene list results stated that timescale difference-related genes were reliable and sensitive, providing additional function-related enrichment information for depression.
There were several limitations to be considered. First, the timescale differences were obtained on a single dataset. However, the stage-specific aberrance was also observed in accelerated brain aging GMV in our previous study and might not result from sample selection (Han et al., 2021). Second, there was a substantial variation across subjects reflecting the individual susceptibility of specific receptor systems (Dukart et al., 2021); future studies should use simultaneous PET and MRI to provide more direct evidence. Third, patients enrolled in our study were under a depressive state. Whether the timescale difference was differently altered in various mood states, such as remitted state (Rive et al., 2015), could be tested in further studies. Fourth, only cross-sectional data were included in this study, and future studies could explore whether altered intrinsic timescale gradient returned to normal with antidepressant treatment especially for fast-acting antidepressants (Li, 2020).
Conclusion
This study revealed atypical intrinsic timescale gradient for the first time. In virtue of brain-wide gene expression and molecular imaging atlases, we investigated possible underlying molecular and transcriptional basis of timescale aberrance linking transcriptome, neurotransmitter information, and neuroimaging findings in depression. These results consistently supported the synaptic dysfunction hypothesis and promoted an integrative understanding of hierarchical dynamics aberrance in depression.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by the Research Ethical Committee of The First Affiliated Hospital of Zhengzhou University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
SH and JC designed the study. RZ, SL, BZ, and YJ collected the data. SH and YC analyzed the data. SH and YZ drafted the work. YW, CW, JP, and HL revised the draft. All authors contributed to the article and approved the submitted version.
Funding
This research study was supported by the Natural Science Foundation of China (81601467, 81871327, and 62106229) and the Medical Science and Technology Research Project of Henan province (201701011, SBGJ202101013, and SBGJ202102103).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank all subjects who participated in this study.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2022.826609/full#supplementary-material
Footnotes
- ^ http://rfmri.org/dpabi
- ^ https://github.com/dominance-analysis/dominance-analysis
- ^ https://github.com/juryxy/JuSpace
- ^ http://human.brain-map.org/
- ^ https://istbi.fudan.edu.cn/lnen/info/1173/1788.htm
References
Alexander-Bloch, A., Giedd, J. N., and Bullmore, E. (2013). Imaging structural co-variance between human brain regions. Nat. Rev. Neurosci. 14, 322–336. doi: 10.1038/nrn3465
Alloy, L. B., Olino, T., Freed, R. D., and Nusslock, R. (2016). Role of reward sensitivity and processing in major depressive and bipolar spectrum disorders. Behav. Ther. 47, 600–621. doi: 10.1016/j.beth.2016.02.014
Altamura, A. C., Dell’Osso, B., Mundo, E., and Dell’Osso, L. (2007). Duration of untreated illness in major depressive disorder: a naturalistic study. Int. J. Clin. Pract. 61, 1697–1700. doi: 10.1111/j.1742-1241.2007.01450.x
Altamura, A. C., Dell’osso, B., Vismara, S., and Mundo, E. (2008). May duration of untreated illness influence the long-term course of major depressive disorder? Eur. Psychiatry 23, 92–96. doi: 10.1016/j.eurpsy.2007.11.004
Anderson, K. M., Krienen, F. M., Choi, E. Y., Reinen, J. M., Yeo, B. T. T., and Holmes, A. J. (2018). Gene expression links functional networks across cortex and striatum. Nat. Commun. 9:1428. doi: 10.1038/s41467-018-03811-x
Arnatkeviciute, A., Fulcher, B. D., and Fornito, A. (2019). A practical guide to linking brain-wide gene expression and neuroimaging data. Neuroimage 189, 353–367. doi: 10.1016/j.neuroimage.2019.01.011
Azen, R., and Budescu, D. V. (2003). The dominance analysis approach for comparing predictors in multiple regression. Psychol. Methods 8, 129–148. doi: 10.1037/1082-989x.8.2.129
Bayer, H. M., and Glimcher, P. W. (2005). Midbrain dopamine neurons encode a quantitative reward prediction error signal. Neuron 47, 129–141. doi: 10.1016/j.neuron.2005.05.020
Beliveau, V., Svarer, C., Frokjaer, V. G., Knudsen, G. M., Greve, D. N., and Fisher, P. M. (2015). Functional connectivity of the dorsal and median raphe nuclei at rest. Neuroimage 116, 187–195. doi: 10.1016/j.neuroimage.2015.04.065
Berman, R. M., Cappiello, A., Anand, A., Oren, D. A., Heninger, G. R., Charney, D. S., et al. (2000). Antidepressant effects of ketamine in depressed patients. Biol. Psychiatry 47, 351–354. doi: 10.1016/s0006-3223(99)00230-9
Bernacchia, A., Seo, H., Lee, D., and Wang, X. J. (2011). A reservoir of time constants for memory traces in cortical neurons. Nat. Neurosci. 14, 366–372. doi: 10.1038/nn.2752
Berridge, K. C. (2007). The debate over dopamine’s role in reward: the case for incentive salience. Psychopharmacology 191, 391–431. doi: 10.1007/s00213-006-0578-x
Blier, P. (2016). Neurobiology of depression and mechanism of action of depression treatments. J. Clin. Psychiatry 77:e319. doi: 10.4088/JCP.13097tx3c
Budescu David, V. (1993). Dominance analysis : a new approach to the problem of relative importance of predictors in multiple regression. Psychol. Bull. 114, 542–551.
Bukh, J. D., Bock, C., Vinberg, M., and Kessing, L. V. (2013). The effect of prolonged duration of untreated depression on antidepressant treatment outcome. J. Affect. Dis. 145, 42–48. doi: 10.1016/j.jad.2012.07.008
Burt, J. B., Demirtaş, M., Eckner, W. J., and Navejar, N. M. (2018). Hierarchy of transcriptomic specialization across human cortex captured by structural neuroimaging topography. Nat. Neurosci. 21, 1251–1259. doi: 10.1038/s41593-018-0195-0
Buzsáki, G., and Draguhn, A. (2004). Neuronal oscillations in cortical networks. Science 304, 1926–1929. doi: 10.1126/science.1099745
Cao, B., Passos, I. C., Mwangi, B., Amaral-Silva, H., and Tannous, J. (2017). Hippocampal subfield volumes in mood disorders. Mol. Psychiatry 22, 1352–1358. doi: 10.1038/mp.2016.262
Cavanagh, S. E., Wallis, J. D., Kennerley, S. W., and Hunt, L. T. (2016). Autocorrelation structure at rest predicts value correlates of single neurons during reward-guided choice. Elife 5:e18937. doi: 10.7554/eLife.18937
Chaudhuri, R., Knoblauch, K., Gariel, M. A., Kennedy, H., and Wang, X. J. (2015). A large-scale circuit mechanism for hierarchical dynamical processing in the primate cortex. Neuron 88, 419–431. doi: 10.1016/j.neuron.2015.09.008
Chen, J., Müller, V. I., Dukart, J., Hoffstaedter, F., Baker, J. T., Holmes, A. J., et al. (2021). Intrinsic connectivity patterns of task-defined brain networks allow individual prediction of cognitive symptom dimension of schizophrenia and are linked to molecular architecture. Biol. Psychiatry 89, 308–319. doi: 10.1016/j.biopsych.2020.09.024
Chen, V. C., Shen, C. Y., Liang, S. H., Li, Z. H., Tyan, Y. S., Liao, Y. T., et al. (2016). Assessment of abnormal brain structures and networks in major depressive disorder using morphometric and connectome analyses. J. Affect. Dis. 205, 103–111. doi: 10.1016/j.jad.2016.06.066
Chen, Y., Cui, Q., Fan, Y. S., Guo, X., Tang, Q., Sheng, W., et al. (2020). Progressive brain structural alterations assessed via causal analysis in patients with generalized anxiety disorder. Neuropsychopharmacology 45, 1689–1697. doi: 10.1038/s41386-020-0704-1
Chen, Z., Peng, W., Sun, H., Kuang, W., Li, W., Jia, Z., et al. (2016). High-field magnetic resonance imaging of structural alterations in first-episode, drug-naive patients with major depressive disorder. Trans. Psychiatry 6:e942. doi: 10.1038/tp.2016.209
Cocchi, L., Sale, M. V., Bell, P. T., Nguyen, V. T., Zalesky, A., and Mattingley, J. B. (2016). A hierarchy of timescales explains distinct effects of local inhibition of primary visual cortex and frontal eye fields. Elife 5:e15252. doi: 10.7554/eLife.15252
Corlett, P. R., Honey, G. D., Krystal, J. H., and Fletcher, P. C. (2011). Glutamatergic model psychoses: prediction error, learning, and inference. Neuropsychopharmacology 36, 294–315. doi: 10.1038/npp.2010.163
Covington, H. E. III, Lobo, M. K., Maze, I., Vialou, V., Hyman, J. M., Zaman, S., et al. (2010). Antidepressant effect of optogenetic stimulation of the medial prefrontal cortex. J. Neurosci. 30, 16082–16090. doi: 10.1523/jneurosci.1731-10.2010
Davidson, L., and McGlashan, T. H. (1997). The varied outcomes of schizophrenia. Can. J. Psychiatry 42, 34–43. doi: 10.1177/070674379704200105
Drevets, W. C., Price, J. L., and Furey, M. L. (2008). Brain structural and functional abnormalities in mood disorders: implications for neurocircuitry models of depression. Brain Struct. Funct. 213, 93–118. doi: 10.1007/s00429-008-0189-x
Dukart, J., Holiga, S., Rullmann, M., and Lanzenberger, R. (2021). JuSpace: a tool for spatial correlation analyses of magnetic resonance imaging data with nuclear imaging derived neurotransmitter maps. Hum. Brain Mapp. 42, 555–566. doi: 10.1002/hbm.25244
Duman, R. S., and Aghajanian, G. K. (2012). Synaptic dysfunction in depression: potential therapeutic targets. Science 338, 68–72. doi: 10.1126/science.1222939
Duman, R. S., Sanacora, G., and Krystal, J. H. (2019). Altered connectivity in depression: GABA and glutamate neurotransmitter deficits and reversal by novel treatments. Neuron 102, 75–90. doi: 10.1016/j.neuron.2019.03.013
Dygalo, N. N., Kalinina, T. S., and Shishkina, G. T. (2020). Stress-induced expression pattern of glutamate signaling genes associated with anhedonia. Stress 23, 700–707. doi: 10.1080/10253890.2020.1812574
Elston, G. N. (2003). Cortex, cognition and the cell: new insights into the pyramidal neuron and prefrontal function. Cereb Cortex 13, 1124–1138. doi: 10.1093/cercor/bhg093
Fan, L., Li, H., Zhuo, J., Zhang, Y., Wang, J., Chen, L., et al. (2016). The human brainnetome atlas: a new brain atlas based on connectional architecture. Cerebral Cortex 26, 3508–3526. doi: 10.1093/cercor/bhw157
Fornito, A., Arnatkevièiûtë, A., and Fulcher, B. D. (2019). Bridging the gap between connectome and transcriptome. Trends Cogn. Sci. 23, 34–50. doi: 10.1016/j.tics.2018.10.005
Fornito, A., Zalesky, A., Bassett, D. S., Meunier, D., Ellison-Wright, I., Yücel, M., et al. (2011). Genetic influences on cost-efficient organization of human cortical functional networks. J. Neurosci. 31, 3261–3270. doi: 10.1523/jneurosci.4858-10.2011
Frodl, T., Meisenzahl, E. M., Zetzsche, T., Born, C., Jäger, M., Groll, C., et al. (2003). Larger amygdala volumes in first depressive episode as compared to recurrent major depression and healthy control subjects. Biol. Psychiatry 53, 338–344. doi: 10.1016/s0006-3223(02)01474-9
Fulcher, B. D., Murray, J. D., Zerbi, V., and Wang, X. J. (2019). Multimodal gradients across mouse cortex. Proc. Natl. Acad. Sci. U.S.A. 116, 4689–4695. doi: 10.1073/pnas.1814144116
Glahn, D. C., Winkler, A. M., Kochunov, P., Almasy, L., Duggirala, R., Carless, M. A., et al. (2010). Genetic control over the resting brain. Proc. Natl. Acad. Sci. U.S.A. 107, 1223–1228. doi: 10.1073/pnas.0909969107
Hamon, M., and Blier, P. (2013). Monoamine neurocircuitry in depression and strategies for new treatments. Prog. Neuropsychopharmacol. Biol. Psychiatry 45, 54–63. doi: 10.1016/j.pnpbp.2013.04.009
Han, S., Chen, Y., Zheng, R., Li, S., Jiang, Y., Wang, C., et al. (2021). The stage-specifically accelerated brain aging in never-treated first-episode patients with depression. Hum. Brain Mapp. 42, 3656–3666. doi: 10.1002/hbm.25460
Han, S., Cui, Q., Wang, X., Chen, Y., Li, D., Li, L., et al. (2020). The anhedonia is differently modulated by structural covariance network of NAc in bipolar disorder and major depressive disorder. Prog. Neuropsychopharmacol. Biol. Psychiatry 99:109865. doi: 10.1016/j.pnpbp.2020.109865
Han, S., He, Z., Duan, X., Tang, Q., Chen, Y., Yang, Y., et al. (2019). Dysfunctional connectivity between raphe nucleus and subcortical regions presented opposite differences in bipolar disorder and major depressive disorder. Prog. Neuropsychopharmacol. Biol. Psychiatry 92, 76–82. doi: 10.1016/j.pnpbp.2018.12.017
Hasson, U., Chen, J., and Honey, C. J. (2015). Hierarchical process memory: memory as an integral component of information processing. Trends Cogn. Sci. 19, 304–313. doi: 10.1016/j.tics.2015.04.006
Hasson, U., Yang, E., Vallines, I., Heeger, D. J., and Rubin, N. (2008). A hierarchy of temporal receptive windows in human cortex. J. Neurosci. 28, 2539–2550. doi: 10.1523/jneurosci.5487-07.2008
Hawrylycz, M. J., Lein, E. S., Guillozet-Bongaarts, A. L., Shen, E. H., Ng, L., Miller, J. A., et al. (2012). An anatomically comprehensive atlas of the adult human brain transcriptome. Nature 489, 391–399. doi: 10.1038/nature11405
Hollerman, J. R., and Schultz, W. (1998). Dopamine neurons report an error in the temporal prediction of reward during learning. Nat. Neurosci. 1, 304–309. doi: 10.1038/1124
Howard, D. M., and Adams, M. J. (2019). Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions. Nat. Neurosci. 22, 343–352. doi: 10.1038/s41593-018-0326-7
Hung, C. I., Yu, N. W., Liu, C. Y., Wu, K. Y., and Yang, C. H. (2015). The impact of the duration of an untreated episode on improvement of depression and somatic symptoms. Neuropsychiatr. Dis. Treat. 11, 2245–2252. doi: 10.2147/ndt.s89498
Huntenburg, J. M., Bazin, P. L., and Margulies, D. S. (2018). Large-scale gradients in human cortical organization. Trends Cogn. Sci. 22, 21–31. doi: 10.1016/j.tics.2017.11.002
Jaako-Movits, K., Zharkovsky, T., Pedersen, M., and Zharkovsky, A. (2006). Decreased hippocampal neurogenesis following olfactory bulbectomy is reversed by repeated citalopram administration. Cell Mol. Neurobiol. 26, 1559–1570. doi: 10.1007/s10571-006-9090-4
Jardri, R., Hugdahl, K., Hughes, M., Brunelin, J., Waters, F., Alderson-Day, B., et al. (2016). Are hallucinations due to an imbalance between excitatory and inhibitory influences on the brain? Schizophr. Bull. 42, 1124–1134. doi: 10.1093/schbul/sbw075
Jiang, Y., Luo, C., Li, X., Duan, M., He, H., Chen, X., et al. (2018). Progressive reduction in gray matter in patients with schizophrenia assessed with MR imaging by using causal network analysis. Radiology 287, 633–642. doi: 10.1148/radiol.2017171832
Kambeitz, J. P., and Howes, O. D. (2015). The serotonin transporter in depression: meta-analysis of in vivo and post mortem findings and implications for understanding and treating depression. J. Affect. Dis. 186, 358–366. doi: 10.1016/j.jad.2015.07.034
Kendler, K. S., Thornton, L. M., and Gardner, C. O. (2001). Genetic risk, number of previous depressive episodes, and stressful life events in predicting onset of major depression. Am. J. Psychiatry 158, 582–586. doi: 10.1176/appi.ajp.158.4.582
Kiebel, S. J., Daunizeau, J., and Friston, K. J. (2008). A hierarchy of time-scales and the brain. PLoS Comput. Biol. 4:e1000209. doi: 10.1371/journal.pcbi.1000209
Klimek, V., Schenck, J. E., Han, H., Stockmeier, C. A., and Ordway, G. A. (2002). Dopaminergic abnormalities in amygdaloid nuclei in major depression: a postmortem study. Biol. Psychiatry 52, 740–748. doi: 10.1016/s0006-3223(02)01383-5
Koutsouleris, N., Davatzikos, C., Borgwardt, S., Gaser, C., Bottlender, R., Frodl, T., et al. (2014). Accelerated brain aging in schizophrenia and beyond: a neuroanatomical marker of psychiatric disorders. Schizophr. Bull. 40, 1140–1153. doi: 10.1093/schbul/sbt142
Kraus, C., and Kadriu, B. (2019). Prognosis and improved outcomes in major depression: a review. Trans. Psychiatry 9:127. doi: 10.1038/s41398-019-0460-3
Krienen, F. M., Yeo, B. T., Ge, T., Buckner, R. L., and Sherwood, C. C. (2016). Transcriptional profiles of supragranular-enriched genes associate with corticocortical network architecture in the human brain. Proc. Natl. Acad. Sci. U.S.A. 113, E469–E478. doi: 10.1073/pnas.1510903113
Kringelbach, M. L., and Cruzat, J. (2020). Dynamic coupling of whole-brain neuronal and neurotransmitter systems. Proc. Natl. Acad. Sci. U.S.A. 117, 9566–9576. doi: 10.1073/pnas.1921475117
Landek-Salgado, M. A., Faust, T. E., and Sawa, A. (2016). Molecular substrates of schizophrenia: homeostatic signaling to connectivity. Mol. Psychiatry 21, 10–28. doi: 10.1038/mp.2015.141
Li, J., and Seidlitz, J. (2021). Cortical structural differences in major depressive disorder correlate with cell type-specific transcriptional signatures. Nat. Commun. 12:1647. doi: 10.1038/s41467-021-21943-5
Li, N., Lee, B., Liu, R. J., Banasr, M., Dwyer, J. M., Iwata, M., et al. (2010). mTOR-dependent synapse formation underlies the rapid antidepressant effects of NMDA antagonists. Science 329, 959–964. doi: 10.1126/science.1190287
Li, N., Liu, R. J., Dwyer, J. M., Banasr, M., Lee, B., Son, H., et al. (2011). Glutamate N-methyl-D-aspartate receptor antagonists rapidly reverse behavioral and synaptic deficits caused by chronic stress exposure. Biol. Psychiatry 69, 754–761. doi: 10.1016/j.biopsych.2010.12.015
Li, X., Zhang, Y., Meng, C., Zhang, C., Zhao, W., Zhu, D. M., et al. (2021). Functional stability predicts depressive and cognitive improvement in major depressive disorder: a longitudinal functional MRI study. Prog. Neuropsychopharmacol. Biol. Psychiatry 111:110396. doi: 10.1016/j.pnpbp.2021.110396
Li, Y. F. (2020). A hypothesis of monoamine (5-HT) - glutamate/GABA long neural circuit: aiming for fast-onset antidepressant discovery. Pharmacol. Ther. 208:107494. doi: 10.1016/j.pharmthera.2020.107494
Limongi, R., Jeon, P., Mackinley, M., Das, T., Dempster, K., Théberge, J., et al. (2020). Glutamate and dysconnection in the salience network: neurochemical, effective connectivity, and computational evidence in schizophrenia. Biol. Psychiatry 88, 273–281. doi: 10.1016/j.biopsych.2020.01.021
Liu, J. J., Hezghia, A., Shaikh, S. R., and Cenido, J. F. (2018). Regulation of monoamine transporters and receptors by lipid microdomains: implications for depression. Neuropsychopharmacology 43, 2165–2179. doi: 10.1038/s41386-018-0133-6
Liu, Z., Rolls, E. T., Liu, Z., Zhang, K., Yang, M., Du, J., et al. (2019). Brain annotation toolbox: exploring the functional and genetic associations of neuroimaging results. Bioinformatics 35, 3771–3778. doi: 10.1093/bioinformatics/btz128
MacQueen, G., and Frodl, T. (2011). The hippocampus in major depression: evidence for the convergence of the bench and bedside in psychiatric research? Mol. Psychiatry 16, 252–264. doi: 10.1038/mp.2010.80
Mayberg, H. S., Lozano, A. M., Voon, V., McNeely, H. E., Seminowicz, D., Hamani, C., et al. (2005). Deep brain stimulation for treatment-resistant depression. Neuron 45, 651–660. doi: 10.1016/j.neuron.2005.02.014
McKinnon, M. C., Yucel, K., Nazarov, A., and MacQueen, G. M. (2009). A meta-analysis examining clinical predictors of hippocampal volume in patients with major depressive disorder. J. Psychiatry Neurosci. 34, 41–54.
Moylan, S., Maes, M., Wray, N. R., and Berk, M. (2013). The neuroprogressive nature of major depressive disorder: pathways to disease evolution and resistance, and therapeutic implications. Mol. Psychiatry 18, 595–606. doi: 10.1038/mp.2012.33
Murray, C. J., Vos, T., Lozano, R., Naghavi, M., Flaxman, A. D., Michaud, C., et al. (2012). Disability-adjusted life years (DALYs) for 291 diseases and injuries in 21 regions, 1990-2010: a systematic analysis for the global burden of disease study 2010. Lancet 380, 2197–2223. doi: 10.1016/s0140-6736(12)61689-4
Phillips, M. L., Ladouceur, C. D., and Drevets, W. C. (2008). A neural model of voluntary and automatic emotion regulation: implications for understanding the pathophysiology and neurodevelopment of bipolar disorder. Mol. Psychiatry 829, 833–857. doi: 10.1038/mp.2008.65
Pirker, W., Asenbaum, S., Kasper, S., Walter, H., Angelberger, P., Koch, G., et al. (1995). beta-CIT SPECT demonstrates blockade of 5HT-uptake sites by citalopram in the human brain in vivo. J. Neural. Trans. Gen. Sect. 100, 247–256. doi: 10.1007/bf01276462
Qiao, H., Li, M. X., Xu, C., Chen, H. B., An, S. C., and Ma, X. M. (2016). Dendritic spines in depression: what we learned from animal models. Neural. Plast 2016:8056370. doi: 10.1155/2016/8056370
Raut, R. V., and Snyder, A. Z. (2020). Hierarchical dynamics as a macroscopic organizing principle of the human brain. Proc. Natl. Acad. Sci. U.S.A. 117, 20890–20897. doi: 10.1073/pnas.2003383117
Reardon, P. K., and Seidlitz, J. (2018). Normative brain size variation and brain shape diversity in humans. Science 360, 1222–1227. doi: 10.1126/science.aar2578
Richiardi, J., Altmann, A., Milazzo, A. C., Chang, C., Chakravarty, M. M., Banaschewski, T., et al. (2015). Brain networks. correlated gene expression supports synchronous activity in brain networks. Science 348, 1241–1244. doi: 10.1126/science.1255905
Rive, M. M., Mocking, R. J. T., Koeter, M. W. J., Wingen, G. V., Wit, S. J. D., Heuvel, O. A. V. D., et al. (2015). State-dependent differences in emotion regulation between unmedicated bipolar disorder and major depressive disorder. JAMA Psychiatry 72:687. doi: 10.1001/jamapsychiatry.2015.0161
Schmaal, L., Hibar, D. P., Sämann, P. G., Hall, G. B., Baune, B. T., Jahanshad, N., et al. (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. Mol. Psychiatry 22, 900–909. doi: 10.1038/mp.2016.60
Shen, X., Tokoglu, F., Papademetris, X., and Constable, R. T. (2013). Groupwise whole-brain parcellation from resting-state fMRI data for network node identification. Neuroimage 82, 403–415. doi: 10.1016/j.neuroimage.2013.05.081
Song, X. W., Dong, Z. Y., Long, X. Y., Li, S. F., Zuo, X. N., Zhu, C. Z., et al. (2011). REST: a toolkit for resting-state functional magnetic resonance imaging data processing. PLoS One 6:e25031. doi: 10.1371/journal.pone.0025031
Stagg, C. J., Bachtiar, V., Amadi, U., Gudberg, C. A., Ilie, A. S., Sampaio-Baptista, C., et al. (2014). Local GABA concentration is related to network-level resting functional connectivity. Elife 3:e01465. doi: 10.7554/eLife.01465
Steinberg, L. J., Underwood, M. D., Bakalian, M. J., Kassir, S. A., Mann, J. J., and Arango, V. (2019). 5-HT1A receptor, 5-HT2A receptor and serotonin transporter binding in the human auditory cortex in depression. J. Psychiatry Neurosci. 44, 294–302. doi: 10.1503/jpn.180190
Tomasi, D., and Volkow, N. D. (2010). Functional connectivity density mapping. Proc. Natl. Acad. Sci. U.S.A. 107, 9885–9890. doi: 10.1073/pnas.1001414107
Treadway, M. T., Waskom, M. L., Dillon, D. G., Holmes, A. J., Park, M. T. M., Chakravarty, M. M., et al. (2015). Illness progression, recent stress, and morphometry of hippocampal subfields and medial prefrontal cortex in major depression. Biol. Psychiatry 77, 285–294. doi: 10.1016/j.biopsych.2014.06.018
Trivedi, M. H., Rush, A. J., Wisniewski, S. R., Nierenberg, A. A., Warden, D., Ritz, L., et al. (2006). Evaluation of outcomes with citalopram for depression using measurement-based care in STAR*D: implications for clinical practice. Am. J. Psychiatry 163, 28–40. doi: 10.1176/appi.ajp.163.1.28
van den Heuvel, M. P., van Soelen, I. L., Stam, C. J., Kahn, R. S., Boomsma, D. I., and Hulshoff Pol, H. E. (2013). Genetic control of functional brain network efficiency in children. Eur. Neuropsychopharmacol. 23, 19–23. doi: 10.1016/j.euroneuro.2012.06.007
van Eijndhoven, P., van Wingen, G., van Oijen, K., Rijpkema, M., Goraj, B., Jan Verkes, R., et al. (2009). Amygdala volume marks the acute state in the early course of depression. Biol. Psychiatry 65, 812–818. doi: 10.1016/j.biopsych.2008.10.027
van Haren, N. E., Cahn, W., Hulshoff Pol, H. E., Schnack, H. G., Caspers, E., Lemstra, A., et al. (2003). Brain volumes as predictor of outcome in recent-onset schizophrenia: a multi-center MRI study. Schizophr. Res. 64, 41–52. doi: 10.1016/s0920-9964(03)00018-5
Vértes, P. E., Rittman, T., Whitaker, K. J., Romero-Garcia, R., Váša, F., and Kitzbichler, M. G. (2016). Gene transcription profiles associated with inter-modular hubs and connection distance in human functional magnetic resonance imaging networks. Philos. Trans. R Soc. Lond. B Biol. Sci. 371:362. doi: 10.1098/rstb.2015.0362
Voineskos, D., Blumberger, D. M., Zomorrodi, R., Rogasch, N. C., Farzan, F., Foussias, G., et al. (2019). Altered transcranial magnetic stimulation-electroencephalographic markers of inhibition and excitation in the dorsolateral prefrontal cortex in major depressive disorder. Biol. Psychiatry 85, 477–486. doi: 10.1016/j.biopsych.2018.09.032
Wang, X. J. (2020). Macroscopic gradients of synaptic excitation and inhibition in the neocortex. Nat. Rev. Neurosci. 21, 169–178. doi: 10.1038/s41583-020-0262-x
Watanabe, T., Rees, G., and Masuda, N. (2019). Atypical intrinsic neural timescale in autism. Elife 8:e42256. doi: 10.7554/eLife.42256
Wengler, K., Goldberg, A. T., Chahine, G., and Horga, G. (2020). Distinct hierarchical alterations of intrinsic neural timescales account for different manifestations of psychosis. Elife 9, e56151. doi: 10.7554/eLife.56151
Whitton, A. E., Treadway, M. T., and Pizzagalli, D. A. (2015). Reward processing dysfunction in major depression, bipolar disorder and schizophrenia. Curr. Opin. Psychiatry 28, 7–12. doi: 10.1097/YCO.0000000000000122
Wong, K. F., and Wang, X. J. (2006). A recurrent network mechanism of time integration in perceptual decisions. J. Neurosci. 26, 1314–1328. doi: 10.1523/jneurosci.3733-05.2006
Wray, N. R., Ripke, S., Mattheisen, M., Trzaskowski, M., Byrne, E. M., Abdellaoui, A., et al. (2018). Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nat. Genet. 50, 668–681. doi: 10.1038/s41588-018-0090-3
Yan, C.-G., Wang, X.-D., Zuo, X.-N., and Zang, Y.-F. (2016). DPABI: data processing & analysis for (resting-state) brain imaging. Neuroinformatics 14, 339–351. doi: 10.1007/s12021-016-9299-4
Yüksel, D., Engelen, J., Schuster, V., Dietsche, B., Konrad, C., Jansen, A., et al. (2018). Longitudinal brain volume changes in major depressive disorder. J. Neural. Trans. 125, 1433–1447. doi: 10.1007/s00702-018-1919-8
Zarate, C. A. Jr., Singh, J. B., Carlson, P. J., Brutsche, N. E., Ameli, R., Luckenbaugh, D. A., et al. (2006). A randomized trial of an N-methyl-D-aspartate antagonist in treatment-resistant major depression. Arch. Gen. Psychiatry 63, 856–864. doi: 10.1001/archpsyc.63.8.856
Zhang, C., Cai, H., and Xu, X. (2021). Genetic architecture underlying differential resting-state functional connectivity of subregions within the human visual cortex. Cerebral Cortex. doi: 10.1093/cercor/bhab335 [Epub ahead of print].
Zhang, Z., Liao, W., Xu, Q., Wei, W., Zhou, H. J., Sun, K., et al. (2017). Hippocampus-associated causal network of structural covariance measuring structural damage progression in temporal lobe epilepsy. Hum. Brain Mapp. 38, 753–766. doi: 10.1002/hbm.23415
Zhou, Y., Zhou, B., Pache, L., and Chang, M. (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun. 10:1523. doi: 10.1038/s41467-019-09234-6
Zhu, D., Yuan, T., Gao, J., Xu, Q., Xue, K., Zhu, W., et al. (2021). Correlation between cortical gene expression and resting-state functional network centrality in healthy young adults. Hum. Brain Mapp. 42, 2236–2249. doi: 10.1002/hbm.25362
Keywords: first-episode depression, gene expression profiling, fMRI, intrinsic timescale gradient, neurotransmitter
Citation: Han S, Zheng R, Li S, Zhou B, Jiang Y, Wang C, Wei Y, Pang J, Li H, Zhang Y, Chen Y and Cheng J (2022) Integrative Functional, Molecular, and Transcriptomic Analyses of Altered Intrinsic Timescale Gradient in Depression. Front. Neurosci. 16:826609. doi: 10.3389/fnins.2022.826609
Received: 01 December 2021; Accepted: 10 January 2022;
Published: 17 February 2022.
Edited by:
Jiajia Zhu, First Affiliated Hospital of Anhui Medical University, ChinaReviewed by:
Le Gao, Yanshan University, ChinaFali Li, University of Electronic Science and Technology of China, China
Junjie Bu, Anhui Medical University, China
Copyright © 2022 Han, Zheng, Li, Zhou, Jiang, Wang, Wei, Pang, Li, Zhang, Chen and Cheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shaoqiang Han, c2hhcWlhbmdoYW5AMTYzLmNvbQ==; Yuan Chen, Y2hlbnl1YW5zaGl6dEAxNjMuY29t; Jingliang Cheng, ZmNjY2hlbmdqbEB6enUuZWR1LmNu