More than six decades have passed since the discovery of monoaminergic antidepressants. Yet, it remains a mystery why these drugs take weeks to months to achieve therapeutic effects, although their monoaminergic actions are present rapidly after treatment. In an attempt to solve this mystery, rather than studying the acute neurochemical effects of antidepressants, here we propose focusing on the early changes in the brain functional connectome using traditional statistics and machine learning approaches. Capitalizing on three independent datasets (n = 1,261) and recent developments in data and network science, we identified a specific connectome fingerprint that predates and predicts response to monoaminergic antidepressants. The discovered fingerprint appears to generalize to antidepressants with differing mechanism of action. We also established a consensus whole-brain hierarchical connectivity architecture and provided a set of model-based features engineering approaches suitable for identifying connectomic signatures of brain function in health and disease.

Overall, we aimed to determine whether early changes in brain functional connectivity predates and predicts response to slow acting antidepressants. To answer this question, we analyzed publicly available data from a relatively large neuroimaging clinical trial in which depressed patients (n = 202) were randomized to placebo or sertraline, a typical monoaminergic slow acting antidepressant (). We then leveraged a large dataset of high-quality fMRI from healthy volunteers (n = 1,003) who participated in the Human Connectome Project () to establish a whole-brain hierarchical parcellation of intrinsic connectivity networks (ICNs). The latter allowed us to extend the cortical findings to subcortical and cerebellar brain regions. Finally, encouraged by the predictive model findings, we conducted a pilot analysis to determine whether the identified connectomic signature related to sertraline response could predate and predict response to the rapid acting antidepressant ketamine, compared with both active and inactive control ().

Machine learning methods have been increasingly implemented in the study of psychopathology (). Among the strengths of machine learning algorithms are the fully data-driven approach and the cross-validation component often included in these predictive models, which could address the issue of over-fitting in interpretive models (e.g., regression) and may enhance generalizability to new data (). However, these machine learning approaches have two critical limitations in the study of psychopathology: (1) They are faced with large number of features (e.g., there are 1,799,970,000 unique edges in a cifti-based fMRI dense connectome) compared with the number of observations, which in psychiatry are often in the order of dozens to hundreds. (2) Considering the dimensionality reduction and weighting procedures involved, it is often not possible to back translate the selected/weighted features to the original space to visualize and understand the underlying neurobiology (). Connectome-based predictive model (CPM) is a linear machine learning approach, which retains the ability to back translate findings to the original feature space (). However, a limitation of CPM is the use of individual parcellated nodes, which leads to a relatively high number of features (e.g., there are 89,676 unique edges in the A424 atlas), which subsequently complicates the interpretability of the findings as the edges included in the final model are often in the order of thousands. Moreover, these nodes are based on anatomical location rather than on network affiliations. To address these limitations, we implemented an NRS-PM that reduces the connectome input features and facilitates the neurobiological interpretation of the final model, as all input features are network based. To assess the robustness of the approach and to rule out potential bias related to architecture selection, we investigated NRS-PM across all architecture levels—with and without subcortical structures, we deconstructed GBC into nodal internal and external NRS, and we determined the full connectome (FC-PM) results using a whole-brain multimodal atlas with 424 nodes (i.e., A424).

Borrowing from studies of complex systems (), NRS approaches address these limitations and have been successfully implemented in the study of psychopathology (). Briefly, the network-restricted approach used in this study is capable of comprehensively assessing all identifiable nodes in a given brain network (e.g., DM), to extract internal connectivity strength within that network and also calculate the external connectivity between networks. To ensure reproducibility in future studies, we opted to use a consensus hierarchical modularity atlas recently established in 1,003 subjects with high-quality functional magnetic resonance imaging (fMRI) data (). Based on the hypothesized hyperconnectivity in DM, we investigated whether reduction in internal DM connectivity would predict better response to sertraline. This was followed by exploratory analysis, with appropriate correction for multiple comparisons, to determine whether response to sertraline is predicted by connectivity within or between brain networks.

The DM comprises brain regions that are synchronously activated at rest and deactivated during external tasks (). Although not without inconsistency (), seed based, independent component, and meta-analysis results have suggested DM hyperconnectivity in MDD compared with controls (). To date, connectivity studies have primarily used seed-based approaches to identify abnormalities or changes in DM of patients with MDD. Among the limitations of this approach are (1) one (or a few) seeds do not cover all the nodes of the DM and the target clusters, hence potentially lowering sensitivity and complicating interpretability of the findings, and (2) it would be difficult to fully integrate findings across studies as they are often highly dependent on the seed location (). Alternatively, whole-brain topology measures address some of these limitations but fail to identify the specific interactions within and between networks (a.k.a., modules, systems, or communities; e.g., DM).

Global brain connectivity (GBC) is a measure of nodal strength within a network. Nodal strength is the most fundamental measure in a graph, as most other topology measures are affected by or are partially derived from nodal strength (). GBC is calculated as the average connectivity between a node (e.g., a gray matter voxel) and all other nodes within a network (e.g., all gray matter voxels) (). In major depressive disorder (MDD), widespread reduction in GBC and other comparable nodal strength measures has been repeatedly demonstrated across samples and by several research groups, particularly in regions within the prefrontal cortex (PFC) (). Moreover, increased GBC in bilateral caudate or lateral PFC following or during infusion with the rapid acting antidepressant ketamine was associated with enhanced response (). Compared with healthy controls, patients with MDD were also found to have high GBC in the posterior cingulate, a critical node within the default mode (DM) network (). This DM GBC dysconnectivity is also normalized following ketamine treatment (). The GBC approach has two essential strengths: (1) nodal strength is a fundamental topology measure within a network, and (2) whole-brain analysis can be conducted, without the limitation and the bias of a priori seed selection. However, a main limitation of GBC is the inability to determine which edges (i.e., connection between two regions) are driving the abnormalities. This limitation is addressed in the current study by the use of network-restricted strength (NRS) and NRS predictive model (NRS-PM) approaches, described in hypotheses 2 and 3. Another limitation of GBC is that it does not capture network changes in opposing directions, e.g., a balanced dynamic network shift from internal to external connectivity cannot be captured using GBC values. To address this issue, we here implemented two measures: nodal internal NRS (niNRS) and nodal external NRS (neNRS).

A consistent evidence in the field is that the behavioral severity of depression significantly improves following placebo treatment, at times leading to lack of behavioral difference between placebo treatment and well-established antidepressants (). The placebo response could be due to nonspecific effects or to the milieu effects of research studies (e.g., repeated visits and assessments). Although the placebo effect is a major impediment for efficacy studies, we here used the significant improvement in depression following placebo to our advantage, as it provided an optimal control for both the test-retest of fMRI and depression measures, as well as controlling for behavioral improvement due to nonspecific and milieu effects. This allowed us to identify the biological correlates of response to the studied antidepressants, controlling for test-retest and nonspecific response. Focusing on the interaction between treatment and response, we implemented data-driven approaches to test three hypotheses.

The serendipitous discovery of slow acting antidepressants in the 1950s has generated persistent interest in identifying the biological underpinnings of depression and unraveling the mechanism of action of antidepressants (). The clinical neuroscience field has since produced a wealth of knowledge related to the biological systems implicated in depression pathophysiology and to the neurochemical effects of these slow acting antidepressants, which tend to modulate monoaminergic neurotransmitters (). Yet, despite more than six decades of research, it remains a mystery as to why the therapeutic behavioral effects of these drugs are only evident following weeks to months of treatment, whereas the neurochemical effects are acutely present after administration (). Solving this mystery may be critical to developing novel efficacious and rapid acting treatments for the large population of patients who are treatment resistant to currently available antidepressants (). In recent years, accumulating evidence implicated brain circuitry and functionally connected networks in the pathology and treatment of depression (). Hence, rather than focusing on the acute synaptic neurochemical effects of monoaminergic antidepressants, it may be more revealing to examine the role of early brain network changes in the mechanisms of these slow acting antidepressants.

To further characterize the ketamine versus lanicemine findings, we conducted whole-brain NRS-PM at AA-24 and AA-50, as well as FC-PM. All three models significantly predicted treatment response to ketamine, compared with lanicemine (r = 0.41 to 0.48, CV = 10, iterations = 1,000, p < 0.05). As shown in Figure 6 , the identified models were largely comparable with those found in the sertraline versus placebo results.

(D) Using nodal strength within the FC-PM as cutoff to retain the highest top 2.5% negative predictive edges and top 2.5% positive predictive edges showed a pattern consistent with the NRS-PM findings, but at the expense of discarding 95% of the data.: The circular graphs in 6A, 6C, and 6D are labeled based on the Akiki-Abdallah (AA) whole-brain architecture at 24 modules (AA-24; see Figure 3 D; Table S1 ), whereas 6B is based on AA-50 (see Figure S3 D; Table S1 ). Modules are colored according to their AA-7 network affiliation (see Figure 3 C). Edges are colored based on the initiating module using a counter-clockwise path starting at 12 o'clock. Internal edges (i.e., within module NRS) are depicted as outer circles around the corresponding module. Thickness of edges reflect their corresponding weight in the predictive model. The module abbreviations of AA-24 and AA-50, along with further details about the affiliation of each node are reported in Table S1 . The predictive models will be made publicly available at https://github.com/emergelab

(C) The full connectome predictive model (FC-PM) predicted response to ketamine compared with active control, but it was not possible to visually discern the underlying signature considering the large number of edges retained in the PM.

(A and B) Network-restricted strength (NRS) predictive models (NRS-PMs; p = 0.003 and p = 0.002, respectively) revealed a comparable connectome fingerprint evident at 20 min during infusion, which predicts enhanced response to ketamine at 24 h, compared with active control. The overall pattern is consistent with findings in the sertraline models ( Figure 4 ), but at a significantly shorter scale (i.e., at 20 min compared with at 7 days).

Following FDR correction for multiple comparisons, the three whole-brain predictive models established in the sertraline study—here termed sertraline connectome fingerprint (CFP)—predicted ketamine response compared with lanicemine at AA-24-CFP (r = 0.52, n = 38, p = 0.0008), AA-50-CFP (r = 0.57, n = 38, p = 0.0002) and full connectome FC-CFP (r = 0.55, n = 38, p = 0.0003). All three sertraline CFPs failed to predict ketamine response, compared with placebo (p > 0.05).

Encouraged by the robust NRS-PM findings, we conducted a pilot follow-up analysis to investigate the generalizability of the identified whole-brain models in predicting the antidepressant response to ketamine, a well-established rapid acting antidepressant (). Here, we used data from a previously published pharmacoimaging study, which examined brain functional connectivity at a period that predates the antidepressant effects of ketamine (i.e., during infusion) (). In this study, 56 patients with MDD were randomized to ketamine, lanicemine, or normal saline. The antidepressant effects of ketamine often arise within hours of its administration. Thus, we examined whether connectivity changes during ketamine infusion predates and predicts response at 24h post treatment, using lanicemine and normal saline as active and inactive control arms, respectively. Ketamine and lanicemine are both N-methyl-D-aspartate (NMDA) receptor antagonists; therefore, the lanicemine arm provided control for both the milieu effect (assessments, infusion, scans, etc.) as well as the acute non-specific NMDA modulation of connectivity networks during infusion ().

Together, the network-restricted interpretative and predictive results supported the nS-PM GBC findings of increased connectivity in nodes within the CE and GPu SC but reduced connectivity in the SM and VI networks as predictor of response ( Figures 1 3 , and 4 ). However, the CE shift from internal to external connectivity cannot be captured by GBC measures, as the latter is an average of both internal and external connectivity. Therefore, to quantify the NRS shift, we deconstructed nS (i.e., nodal GBC) into two complementary measures: (1) nodal internal NRS (niNRS) calculated as the average connectivity between each node and all other nodes within the same ICN and (2) nodal external NRS (neNRS) calculated as the average connectivity between each node and all other nodes outside its ICN. Here, we used AA-7, which comprises the main brain networks while incorporating subcortical structures, including CE, DM, VS, DS, SC, SM, and VI networks ( Figure 3 C). Changes in neNRS at week 1 significantly predicted improvement at week 8 following sertraline treatment compared with placebo (r = 0.28, CV = 10, iterations = 1,000, p ≤ 0.001). Enhanced antidepressant response was associated with increased neNRS in brain regions primarily located within the CE and GPu SC networks ( Figures 5 A and 5B ). Similarly, changes in niNRS significantly predicted the antidepressant response (r = 0.24, CV = 10, iterations = 1,000, p = 0.005). Enhanced antidepressant response was associated with reduced niNRS in brain regions primarily located within the CE, GPu SC, SM, and VI networks ( Figures 5 C and 5D). These results quantitatively demonstrated the internal-to-external NRS shift as predictor of antidepressant response. The cortical and whole-brain neNRS-PMs and niNRS-PMs across all architectures are provided in Supplemental Information Figures S6 and S7 ). Similar to the NRS-PMs, including subcortical structures resulted in more consistent neNRS-PMs and niNRS-PMs across architectures ( Figures S6 and S7 ).

(C and D) A predictive model (p = 0.005) using nodal internal network-restricted strength (niNRS) as input, based on AA-7, predicted enhanced response to sertraline. The model showed reduced internal connectivity in brain regions within the CE, GPu SC, sensorimotor, and visual networks. There were no increases in internal connectivity. Red arrows point to regions that showed both reduced niNRS and increase neNRS, all of which located within the CE and GPu SC networks, consistent with the internal-to-external connectivity shift observed in previous analyses and quantitively supporting the observed pattern of a shift toward increased higher order control over primary cortices.

(A and B) A predictive model (p ≤ 0.001) using nodal external network-restricted strength (neNRS) as input, based on Akiki-Abdallah (AA) whole-brain atlas at the architecture with 7 modules (AA-7), predicted enhanced response to sertraline. The model showed increased global external connectivity in brain regions within the central executive (CE) and globus-pallidus-putamen subcortical (GPu SC) networks. There were no reductions in global external connectivity.

We next investigated whether a nodal based (i.e., 424x424 nodes) full connectome PM (FC-PM) would yield comparable or differing results, compared with whole-brain NRS-PM. The FC-PM significantly predicted the antidepressant response (r = 0.27, CV = 10, iterations = 1,000, p = 0.005). However, considering the large number of predictive edges, visualizing the PM is unable to discern the underlying connectomic signature ( Figure 4 C). Retaining the nodes with the highest degree (i.e., top 2.5% of each of positive and negative predictive edges) showed a pattern of reduced internal connectivity among nodes in the SM and VI network, but at the expense of discarding 95% of the model, and failed to reveal the internal-to-external shift within the CE network ( Figure 4 D).

To further assess the whole-brain NRS-PM results, the AA-24 and AA-50 were visualized in Figures 4 A and 4B , respectively. Predictive edges were comparable with findings with cortical atlas (see Supplemental Information and Figure S4 ), with richer connectomic signature and two highlights: (1) Modules within the CE and GPu SC ICNs showed reduced internal connectivity among each other but increased connectivity to the rest of the brain as predictors of enhanced antidepressant response. (2) Reduction in internal connectivity within the remaining ICNs, particularly SM and VI, also predicted better sertraline response. Modules containing the amygdala and insula also showed a shift from connection with primary cortices (e.g., sensory, motor, and visual) to increased connectivity with higher order association areas ( Figure 4 ).

(D) Using nodal strength within the FC-PM as cutoff to retain the highest top 2.5% negative predictive edges and top 2.5% positive predictive edges showed a pattern consistent with the NRS-PM findings, but at the expense of discarding 95% of the data, and did not fully depict the shift from internal to external strength within the CE and GPu SC networks.: The circular graphs in 4A, 4C, and 4D are labeled based on the Akiki-Abdallah (AA) whole-brain architecture at 24 modules (AA-24; see Figure 3 D; Table S1 ), whereas 4B is based on AA-50 (see Figure S3 D; Table S1 ). Modules are colored according to their AA-7 network affiliation (see Figure 3 C). Edges are colored based on the initiating module using a counter-clockwise path starting at 12 o'clock. Internal edges (i.e., within module NRS) are depicted as outer circles around the corresponding module. Thickness of edges reflects their corresponding weight in the predictive model. The module abbreviations of AA-24 and AA-50, along with further details about the affiliation of each node, are reported in Table S1 . The predictive models will be made publicly available at https://github.com/emergelab

(C) The full connectome predictive model (FC-PM; p = 0.005) predicted response to sertraline, but it was not possible to visually discern the underlying signature considering the large number of edges retained in the PM.

(A and B) Network-restricted strength (NRS) predictive models (NRS-PMs; p = 0.032 and p = 0.033, respectively) revealed a specific connectomic signature evident at week 1 post treatment, at a time of no clinically meaningful antidepressant effects, that predicts enhanced response to sertraline at week 8, compared with placebo. An overall pattern emerged consistent of reduced connectivity between modules within the central executive (CE) and globus-pallidus-putamen subcortical (GPu SC) networks, along with increased connectivity between these two networks (CE/GPu SC) and the rest of the brain as predictors of enhanced response to sertraline. The connectome fingerprint (CFP) also showed reduction in internal connectivity within the visual and sensorimotor networks, as well as reduced interference on these two networks from modules within the default mode and salience networks.

Following ICN hierarchical parcellation, we aimed to determine whether a whole-brain connectome fingerprint predates and predicts response to sertraline, compared with placebo. As shown in Figure S5 B, the whole-brain NRS-PM predicted the antidepressant response across AA-4 to AA-150 architectures (following FDR correction), with a peak at AA-58 (r = 0.27, CV = 10, iterations = 1,000, p = 0.003). Independently, the positive predictive edges peaked at AA-58 (r = 0.29, CV = 10, iterations = 1,000, p = 0.001) and the negative predictive edges peaked at AA-26 (r = 0.25, CV = 10, iterations = 1,000, p = 0.003). Compared with cortical NRS-PM, the inclusion of subcortical structure appears to have enhanced models, particularly positive predictions that were more consistent in independently predicting antidepressant response across architectures ( Figure S5 ).

To identify the ICN affiliation of the subcortical GBC findings and to extend the cortical NRS to subcortical and cerebellar regions, we next aimed to determine the hierarchical architecture of the whole brain using subject-level clustering of functional networks from 1,003 adult healthy subjects (). Using consistent definition of community across cortical, subcortical, and cerebellar structures, the whole-brain subject-level clustering and consensus community detection algorithm identified 136 significant architecture ( Figures 3 and S3 ). The hierarchical atlas, here termed Akiki-Abdallah (AA) atlas, delineated architectures ranging from 3 to 150 modules. The cortical modules remained largely consistent with the cortex-only atlases in our previous work () and in the literature (), as shown in Figures 2 A versus 2C. The AA architecture at 24 modules (AA-24) largely matches the cortical 22 modules architecture (see Supplemental Information and Figure S4 ) while adding the affiliation of subcortical and cerebellar nodes. AA-50 is the architecture with the highest similarity to subject-specific modularity ( Figures 3 B and S3 D). Notably, the ICN affiliations of the positive predictive GBC nodes findings were within the CE and global-pallidus-putamen subcortical (GPu SC) modules, whereas the negative predictive GBC nodes were within the SM and VI networks (see Figures 1 B versus 3 C).

(C and D) These maps show the whole-brain networks nodal affiliation at AA-7 (corresponds to the six main cortical networks) and AA-24 (corresponds to the peak subject-level similarity in the cortical atlas). The module abbreviations of AA-24, along with further details about the affiliation of each node, are reported in Table S1 . The figures were adapted with permission from the Emerge Research Program ( http://emerge.care ). The hierarchical atlas maps and codes will be made publicly available at https://github.com/emergelab

(B) Similarity plot showing the mean similarity between the partitioning in each consensus hierarchical level (i.e., from AA-3 to AA-150) and the subject-level clustering quantified by Z score of the Rand coefficient (blue line); the Rand score peaked at AA-50. The dashed red lines denote the brain architecture levels for Figures 3 D and S3 D.

These results confirmed the predictions of hypothesis 2. Yet, to better characterize the DM findings and inform future studies, we conducted a follow-up analysis, with FDR correction, examining internal and external connectivity of all cortical ICNs ( Figure 2 A, i.e., DM, CE, dorsal salience [DS], ventral salience [VS], SM, and VI). Following correction for multiple comparisons, there was significant treatment-by-response interactive effect on the CE-SM edge (F= 10.4, p = 0.001), reflecting increased connectivity predicts better response to sertraline compared with placebo. There was also significant treatment effect on the CE-SM edge (F= 11.0, p = 0.001), reflecting increased connectivity at week 1 following sertraline, compared with placebo ( Figure S2 ). There were no other significant treatment, response, or treatment-by-response effects (q > 0.05).

(B) The network-restricted strength (NRS) pentagon. Internal NRS is depicted as filled circles, whereas inter-networks external NRS is depicted as edges. * was used for p < 0.05, ** for p < 0.01, *** for p < 0.001. Interactions that survived FDR correction were denoted with squares (i.e., CE-SM edge). Filled circles and edges were colored gray for non-significant effects, blue for negative effects, and red for positive effects. The NRS interpretive analysis examining the interactive effects between treatment and response showed significant increase in connectivity between CE and SM at week 1 that predicts response to sertraline at week 8. This figure was adapted with permission from the Emerge Research Program ( http://emerge.care ).

Employing validated NRS methods (), we constructed a general linear model examining the effects of treatment, response, and treatment-by-response on changes in internal DM NRS (i.e., week 1 minus pretreatment). We found a significant treatment-by-response interaction (F= 5.0, p = 0.026), such that, compared with placebo, reduction in DM NRS at week 1 predicted better response to sertraline at week 8. There were no significant response or treatment effects (p > 0.1).

These results confirmed the predicted role of caudate GBC but failed to show an association between antidepressant response and increased lateral PFC GBC (hypothesis 1). We speculated that this failure may be due to the limitations (see Supplemental Information ) of the vertex-/voxel-wise interpretive analysis. To rule out this possibility, we implemented a nodal strength predictive model (nS-PM), with 1,000 iterations of 10 cross-validation (CV), to identify brain regions that significantly predict a continuous measure of improvement following sertraline compared with placebo. Changes in nodal GBC (i.e., nS) at week 1 significantly predicted improvement at week 8 following sertraline compared with placebo (r = 0.23, CV = 10, iterations = 1,000, p ≤ 0.001). Enhanced antidepressant response was associated with increased nodal GBC in bilateral caudate and left lateral PFC, along with other brain regions primarily located within the central executive (CE) ICN ( Figure 1 B). Reduced nodal GBC in several regions within the visual (VI) and sensorimotor (SM) ICNs also predicted better antidepressant response, compared with placebo ( Figure 1 B).

Guided by previous work with antidepressants (), we first conducted a whole-brain vertex-/voxel-wise GBC analysis, with false discovery rate (FDR) correction (q < 0.05), examining the interaction between treatment (sertraline versus placebo) and response (responders versus non-responders). We found that increased GBC in bilateral caudate and right rostral anterior cingulate (rACC) at week 1 were associated with better response to sertraline at week 8 compared with placebo ( Figure 1 A). Two clusters with reduced GBC in left Brodmann areas 4 (motor) and 5 (somatosensory) were associated with better response to sertraline ( Figure 1 A). The independent effects of response and treatment are provided in the Supplemental Information Figure S1 ).

(A) Whole-brain interpretive analysis showed a significant interaction between treatment and response, such that increased global brain connectivity (GBC) in the bilateral caudate and right rostral anterior cingulate (red circles) and decreased GBC in sensorimotor cortices (blue circles) at week 1 predicted response to sertraline at week 8 post treatment. The color bar represents the z values (threshold at p < 0.005, with circles denoting clusters that survived FDR correction at q < 0.05).

The sertraline dataset was acquired from the National Institute of Mental Health Data Archive (NDA), Establishing Moderators and Biosignatures of Antidepressant Response for Clinical Care (EMBARC). All participants (n = 202) with successful resting state fMRI at baseline and week 1, who completed depression assessment at week 8, were included in the current study. The protocol and results of the EMBARC clinical trial were reported elsewhere (). Briefly, unmedicated patients with chronic or recurrent MDD were randomized to 8 weeks of daily oral placebo or sertraline. Demographics and clinical features are described in Table 1 . Hamilton rating scale for depression (HAMD) was used to determine severity. Response was defined as at least 50% improvement in HAMD at week 8, the time point at which full clinical benefit is expected following slow acting antidepressant treatment.

Discussion

Using a set of traditional statistics and machine learning approaches, the study results provided strong evidence of a specific connectome fingerprint (CFP) that predates and predicts response to the slow acting antidepressant, sertraline. The study established a whole-brain hierarchical ICNs atlas and provided evidence of its relevance to the study of psychopathology and to clinical neuroscience discovery. It also presented specialized measures and tested the robustness of innovative approaches, i.e., NRS-PM, niNRS-PM, and neNRS-PM. The impact of these atlases, measures, and approaches is expected to go beyond the current findings, as similar approaches could be used to identify the CFPs of various brain functions in health and disease. Finally, the study provided pilot evidence about the potential generalizability, albeit at a much shorter timescale, of the identified CFP to the mechanisms of rapid acting antidepressants.

Partially supporting hypothesis 1, bilateral caudate clusters of increased global connectivity at week 1 post treatment predicted better response to sertraline compared with placebo at week 8 ( Figure 1 A). Treatment response was also associated with decreased global connectivity in the sensory motor area but increased connectivity in the rostral anterior cingulate cortex (ACC; Figure 1 A). However, the latter finding lost significance in the predictive model analysis ( Figure 1 B). Intriguingly, although GBC is a whole-brain based measure and not network restricted, the positive predictive GBC nodes were primarily located within the CE and GPu SC modules. Similarly, the negative predictive GBC nodes were primarily within the SM and VI networks (see Figures 1 B vs. 3 C). Although the primary analysis failed to show increased GBC in the lateral PFC as predicted by hypothesis 1, this may be due to limitations of the used vertex-/voxel-wise approach (see below) Supplemental Information.

The interpretive NRS model confirmed the predictions of hypothesis 2, showing significant association between better response to sertraline at week 8 and earlier reduction in internal DM connectivity at week 1 post treatment. Follow-up analyses revealed increased CE-SM connectivity as early predictor of enhanced response to sertraline ( Figure 2 B).

Compared with placebo, the NRS-PMs significantly predicted response to sertraline, supporting hypothesis 3, while revealing dynamic shifts in the brain circuits, critically underlying the negative and positive predictions of treatment effects ( Figures 4 and 5 ). In particular, three patterns of connectivity shifts have emerged as predictors of response to antidepressant treatment: (1) A reduction in internal connectivity among the CE and GPu SC modules, along with increased external connectivity between these modules and the rest of the brain (most evident in AA-50, but also in AA-24, and niNRS/neNRS). (2) Reduced internal connectivity among modules within the SM and VI networks. (3) In DM/VS modules containing the amygdala and insula (i.e., affective DM and inferior VS, respectively), there was reduced connectivity with perceptual and motor areas (i.e., SM and VI) but increased connectivity with higher order association regions, suggesting an early shift toward enhanced executive control Supplemental Information.

Yu et al., 2019 Yu M.

Linn K.A.

Shinohara R.T.

Oathes D.J.

Cook P.A.

Duprat R.

Moore T.M.

Oquendo M.A.

Phillips M.L.

McInnis M.

et al. Childhood trauma history is linked to abnormal brain connectivity in major depression. Yu et al., 2019 Yu M.

Linn K.A.

Shinohara R.T.

Oathes D.J.

Cook P.A.

Duprat R.

Moore T.M.

Oquendo M.A.

Phillips M.L.

McInnis M.

et al. Childhood trauma history is linked to abnormal brain connectivity in major depression. The current study focused on the changes overtime in brain functional connectivity and in psychopathology as reflected by antidepressant scores. This focused approach provided a powerful comprehensive assessment of brain connectivity in relation to treatment response. However, future complementary and follow-up studies are still required to determine whether pretreatment connectivity, as well as clinical features (e.g., early life stress or symptom clusters) would be associated with specific CFPs or would predict response to antidepressants. Notably, an outstanding recent study has investigated pretreatment connectivity and clinical features in the EMBARC cohort (). Before treatment, it was reported that various brain networks correlated with early life stressors and that MDD was characterized by increased NRS in intrinsic networks, including DM, but reduced NRS in task-positive networks, including CE (). Importantly, these pretreatment findings of altered DM and CE connectivity in MDD suggest that the antidepressant CFP identified in the current study may reflect an early pattern of normalization detectable biologically at week 1 but not evident behaviorally until week 8.

Abdallah et al., 2017b Abdallah C.G.

Averill L.A.

Collins K.A.

Geha P.

Schwartz J.

Averill C.

DeWilde K.E.

Wong E.

Anticevic A.

Tang C.Y.

et al. Ketamine treatment and global brain connectivity in major depression. Abdallah et al., 2018a Abdallah C.G.

Dutta A.

Averill C.L.

McKie S.

Akiki T.J.

Averill L.A.

William Deakin J. Ketamine, but not the NMDAR antagonist lanicemine, increases prefrontal global connectivity in depressed patients. Yeo et al., 2011 Yeo B.T.

Krienen F.M.

Sepulcre J.

Sabuncu M.R.

Lashkari D.

Hollinshead M.

Roffman J.L.

Smoller J.W.

Zollei L.

Polimeni J.R.

et al. The organization of the human cerebral cortex estimated by intrinsic functional connectivity. Choi et al., 2012 Choi E.Y.

Yeo B.T.

Buckner R.L. The organization of the human striatum estimated by intrinsic functional connectivity. Power et al., 2011 Power J.D.

Cohen A.L.

Nelson S.M.

Wig G.S.

Barnes K.A.

Church J.A.

Vogel A.C.

Laumann T.O.

Miezin F.M.

Schlaggar B.L.

Petersen S.E. Functional network organization of the human brain. In the context of inconsistent patterns in the literature, we did not a priori predict the widespread role of the CE network at multiple hierarchical architectures, even though previous results from our own studies have associated increased GBC in the lateral PFC and caudate with enhanced antidepressant response (). Better understanding of this oversight may be essential to enhance reproducibility in psychiatric neuroimaging studies. We believe that this oversight was primarily due to two common limitations of neuroimaging work. The first is that vertex-/voxel-wise traditional statistics (i.e., interpretive models, e.g., regression) are subject to overfitting owing to the lack of cross-validation and to type II error owing to the needed correction for multiple comparison. For example, an increase in vertex-wise GBC in the rostral ACC was identified in the interpretive analysis ( Figure 1 A) but did not hold when subjected to cross-validation procedures ( Figure 1 B). Moreover, the vertex-/voxel-wise interpretive model identifies only the peak effect (i.e., the vertices/voxels with the highest effect size), without the ability to fully map the behavioral brain effects (see Figures 1 A vs. 1B). The second limitation is that the majority of ICN atlases were not multiscale hierarchical (i.e., lacking the information about the upstream ICN architecture) and were limited to the cortex (e.g.,). The few ICN atlases that included subcortical structures were mostly post hoc (i.e., projecting the cortical ICNs on the subcortical regions without the use of a unified definition of ICN community; e.g.,) or were meta-analytically based missing important subcortical structures (e.g., the amygdala in). Thus, it was not commonly viewed that the caudate and frontoparietal cortex share the same upstream CE architecture. Furthermore, the biased focus of the literature on few key seeds and ICNs (e.g., DM) may have contributed to the apparent inconsistency across studies. In the current study, the full assessment of the connectome along with combined use of predictive models, hierarchical ICN architecture, and unified definition of ICN across cortical and subcortical nodes led to the robust identification of the CE ICN as critical predictor of treatment response using multiple data-driven analytical approaches.

In summary, based on the behavioral assessment of depression, there was no clinically meaningful behavioral response at week 1 (∼2 points reduction on HAMD) and there were no significant differences between the treatment arms. In contrast, the biological functional connectivity investigation at week 1 showed a robust connectomic signature that predates and predicts response to sertraline at week 8. The identified biosignature was evident using three distinct, yet overlapping and complimentary, approaches. A main impact of the findings is the identification of an antidepressant-related dynamic shift in brain networks connectivity, which is consistent with early increase in executive control weeks prior to the full therapeutic behavioral effects. Some features of the identified connectomic signature are reminiscent of findings with ketamine, and the predictive models at least partially generalized to ketamine treatment, raising the possibility that the discovered connectome fingerprint may generalize to antidepressants with differing mechanism of action.

Future studies can capitalize on the established set of approaches that combine the strengths of consensus hierarchical architecture of the brain, along with neuroscience-informed features engineering and machine learning to fully assess the functional connectome and successfully identify psychopathology-relevant reproducible fingerprint. These data-driven approaches are likely to have impact beyond answering the research question of the current study, by allowing the full assessment of functional connectivity in relation to behavior without the reliance on a limited number of seeds or networks. Together, the study findings enriched our understanding of the neurobiology of depression and revealed a previously unknown connectomic signature that may serve as a treatment target, as a biological indicator of response to optimize treatment regimen, or as a surrogate for early stages in drug development.