Dynamic Functional Connectivity Predicts Treatment Response to Electroconvulsive Therapy in Major Depressive Disorder

The overlap across psychotic disorders: A functional network connectivity analysis

Functional network connectivity (FNC) has previously been shown to distinguish patient groups from healthy controls (HC). However, the overlap across psychiatric disorders such as schizophrenia (SZ), bipolar (BP), and schizoaffective disorder (SAD) is not evident yet. This study focuses on studying the overlap across these three psychotic disorders in both dynamic and static FNC (dFNC/sFNC). We used resting-state fMRI, demographics, and clinical information from the Bipolar-Schizophrenia Network on Intermediate Phenotypes cohort (BSNIP). The data includes three groups of patients with schizophrenia (SZ, N = 181), bipolar (BP, N = 163), and schizoaffective (SAD, N = 130) and HC (N = 238) groups. After estimating each individual's dFNC, we group them into three distinct states. We evaluated two dFNC features, including occupancy rate (OCR) and distance travelled over time. Finally, the extracted features, including both sFNC and dFNC, are tested statistically across patients and HC groups. In addition, we explored the link between the clinical scores and the extracted features. We evaluated the connectivity patterns and their overlap among SZ, BP, and SAD disorders (false discovery rate or FDR corrected p < 0.05). Results showed dFNC captured unique information about overlap across disorders where all disorder groups showed similar pattern of activity in state 2. Moreover, the results showed similar patterns between SZ and SAD in state 1 which was different than BP. Finally, the distance travelled feature of SZ (average R = 0.245, p < 0.01) and combined distance travelled from all disorders was predictive of the PANSS symptoms scores (average R = 0.147, p < 0.01).

Searching Reproducible Brain Features using NeuroMark: Templates for Different Age Populations and Imaging Modalities

A primary challenge to the data-driven analysis is the balance between poor generalizability of population-based research and characterizing more subject-, study- and population-specific variability. We previously introduced a fully automated spatially constrained independent component analysis (ICA) framework called NeuroMark and its functional MRI (fMRI) template. NeuroMark has been successfully applied in numerous studies, identifying brain markers reproducible across datasets and disorders. The first NeuroMark template was constructed based on young adult cohorts. We recently expanded on this initiative by creating a standardized normative multi-spatial-scale functional template using over 100,000 subjects, aiming to improve generalizability and comparability across studies involving diverse cohorts. While a unified template across the lifespan is desirable, a comprehensive investigation of the similarities and differences between components from different age populations might help systematically transform our understanding of the human brain by revealing the most well-replicated and variable network features throughout the lifespan. In this work, we introduced two significant expansions of NeuroMark templates first by generating replicable fMRI templates for infants, adolescents, and aging cohorts, and second by incorporating structural MRI (sMRI) and diffusion MRI (dMRI) modalities. Specifically, we built spatiotemporal fMRI templates based on 6,000 resting-state scans from four datasets. This is the first attempt to create robust ICA templates covering dynamic brain development across the lifespan. For the sMRI and dMRI data, we used two large publicly available datasets including more than 30,000 scans to build reliable templates. We employed a spatial similarity analysis to identify replicable templates and investigate the degree to which unique and similar patterns are reflective in different age populations. Our results suggest remarkably high similarity of the resulting adapted components, even across extreme age differences. With the new templates, the NeuroMark framework allows us to perform age-specific adaptations and to capture features adaptable to each modality, therefore facilitating biomarker identification across brain disorders. In sum, the present work demonstrates the generalizability of NeuroMark templates and suggests the potential of new templates to boost accuracy in mental health research and advance our understanding of lifespan and cross-modal alterations.

Improved clinical outcome prediction in depression using neurodynamics in an emotional face-matching functional MRI task

Introduction

Approximately one in six people will experience an episode of major depressive disorder (MDD) in their lifetime. Effective treatment is hindered by subjective clinical decision-making and a lack of objective prognostic biomarkers. Functional MRI (fMRI) could provide such an objective measure but the majority of MDD studies has focused on static approaches, disregarding the rapidly changing nature of the brain. In this study, we aim to predict depression severity changes at 3 and 6 months using dynamic fMRI features.

Methods

For our research, we acquired a longitudinal dataset of 32 MDD patients with fMRI scans acquired at baseline and clinical follow-ups 3 and 6 months later. Several measures were derived from an emotion face-matching fMRI dataset: activity in brain regions, static and dynamic functional connectivity between functional brain networks (FBNs) and two measures from a wavelet coherence analysis approach. All fMRI features were evaluated independently, with and without demographic and clinical parameters. Patients were divided into two classes based on changes in depression severity at both follow-ups.

Results

The number of coherence clusters (nCC) between FBNs, reflecting the total number of interactions (either synchronous, anti-synchronous or causal), resulted in the highest predictive performance. The nCC-based classifier achieved 87.5% and 77.4% accuracy for the 3- and 6-months change in severity, respectively. Furthermore, regression analyses supported the potential of nCC for predicting depression severity on a continuous scale. The posterior default mode network (DMN), dorsal attention network (DAN) and two visual networks were the most important networks in the optimal nCC models. Reduced nCC was associated with a poorer depression course, suggesting deficits in sustained attention to and coping with emotion-related faces. An ensemble of classifiers with demographic, clinical and lead coherence features, a measure of dynamic causality, resulted in a 3-months clinical outcome prediction accuracy of 81.2%.

Discussion

The dynamic wavelet features demonstrated high accuracy in predicting individual depression severity change. Features describing brain dynamics could enhance understanding of depression and support clinical decision-making. Further studies are required to evaluate their robustness and replicability in larger cohorts.

Brain dynamics reflecting an intra-network brain state is associated with increased posttraumatic stress symptoms in the early aftermath of trauma

This study examines the association between brain dynamic functional network connectivity (dFNC) and current/future posttraumatic stress (PTS) symptom severity, and the impact of sex on this relationship. By analyzing 275 participants' dFNC data obtained ~2 weeks after trauma exposure, we noted that brain dynamics of an inter-network brain state link negatively with current (r=-0.179, pcorrected= 0.021) and future (r=-0.166, pcorrected= 0.029) PTS symptom severity. Also, dynamics of an intra-network brain state correlated with future symptom intensity (r = 0.192, pcorrected = 0.021). We additionally observed that the association between the network dynamics of the inter-network brain state with symptom severity is more pronounced in females (r=-0.244, pcorrected = 0.014). Our findings highlight a potential link between brain network dynamics in the aftermath of trauma with current and future PTSD outcomes, with a stronger protective effect of inter-network brain states against symptom severity in females, underscoring the importance of sex differences.

Explainable fuzzy clustering framework reveals divergent default mode network connectivity dynamics in schizophrenia

Introduction

Dynamic functional network connectivity (dFNC) analysis of resting state functional magnetic resonance imaging data has yielded insights into many neurological and neuropsychiatric disorders. A common dFNC analysis approach uses hard clustering methods like k-means clustering to assign samples to states that summarize network dynamics. However, hard clustering methods obscure network dynamics by assuming (1) that all samples within a cluster are equally like their assigned centroids and (2) that samples closer to one another in the data space than to their centroids are well-represented by their centroids. In addition, it can be hard to compare subjects, as in some cases an individual may not manifest a state strongly enough to enter a hard cluster. Approaches that allow a dimensional approach to connectivity patterns (e.g., fuzzy clustering) can mitigate these issues. In this study, we present an explainable fuzzy clustering framework by combining fuzzy c-means clustering with several explainability metrics and novel summary features.

Methods

We apply our framework for schizophrenia (SZ) default mode network analysis. Namely, we extract dFNC from individuals with SZ and controls, identify 5 dFNC states, and characterize the dFNC features most crucial to those states with a new perturbation-based clustering explainability approach. We then extract several features typically used in hard clustering and further present a variety of unique features specially designed for use with fuzzy clustering to quantify state dynamics. We examine differences in those features between individuals with SZ and controls and further search for relationships between those features and SZ symptom severity.

Results

Importantly, we find that individuals with SZ spend more time in states of moderate anticorrelation between the anterior and posterior cingulate cortices and strong anticorrelation between the precuneus and anterior cingulate cortex. We further find that individuals with SZ tend to transition more rapidly than controls between low-magnitude and high-magnitude dFNC states.

Conclusion

We present a novel dFNC analysis framework and use it to identify effects of SZ upon network dynamics. Given the ease of implementing our framework and its enhanced insight into network dynamics, it has great potential for use in future dFNC studies.

Manifold learning for fMRI time-varying functional connectivity

Whole-brain functional connectivity (FC) measured with functional MRI (fMRI) evolves over time in meaningful ways at temporal scales going from years (e.g., development) to seconds [e.g., within-scan time-varying FC (tvFC)]. Yet, our ability to explore tvFC is severely constrained by its large dimensionality (several thousands). To overcome this difficulty, researchers often seek to generate low dimensional representations (e.g., 2D and 3D scatter plots) hoping those will retain important aspects of the data (e.g., relationships to behavior and disease progression). Limited prior empirical work suggests that manifold learning techniques (MLTs)-namely those seeking to infer a low dimensional non-linear surface (i.e., the manifold) where most of the data lies-are good candidates for accomplishing this task. Here we explore this possibility in detail. First, we discuss why one should expect tvFC data to lie on a low dimensional manifold. Second, we estimate what is the intrinsic dimension (ID; i.e., minimum number of latent dimensions) of tvFC data manifolds. Third, we describe the inner workings of three state-of-the-art MLTs: Laplacian Eigenmaps (LEs), T-distributed Stochastic Neighbor Embedding (T-SNE), and Uniform Manifold Approximation and Projection (UMAP). For each method, we empirically evaluate its ability to generate neuro-biologically meaningful representations of tvFC data, as well as their robustness against hyper-parameter selection. Our results show that tvFC data has an ID that ranges between 4 and 26, and that ID varies significantly between rest and task states. We also show how all three methods can effectively capture subject identity and task being performed: UMAP and T-SNE can capture these two levels of detail concurrently, but LE could only capture one at a time. We observed substantial variability in embedding quality across MLTs, and within-MLT as a function of hyper-parameter selection. To help alleviate this issue, we provide heuristics that can inform future studies. Finally, we also demonstrate the importance of feature normalization when combining data across subjects and the role that temporal autocorrelation plays in the application of MLTs to tvFC data. Overall, we conclude that while MLTs can be useful to generate summary views of labeled tvFC data, their application to unlabeled data such as resting-state remains challenging.

Exploring the Neural Processes behind Narrative Engagement: An EEG Study

Past cognitive neuroscience studies using naturalistic stimuli have considered narratives holistically and focused on cognitive processes. In this study, we incorporated the narrative structure, the dramatic arc, as an object of investigation, to examine how engagement levels fluctuate across a narrative-aligned dramatic arc. We explored the possibility of predicting self-reported engagement ratings from neural activity and investigated the idiosyncratic effects of each phase of the dramatic arc on brain responses as well as the relationship between engagement and brain responses. We presented a movie excerpt following the six-phase narrative arc structure to female and male participants while collecting EEG signals. We then asked this group of participants to recall the excerpt, another group to segment the video based on the dramatic arc model, and a third to rate their engagement levels while watching the movie. The results showed that the self-reported engagement ratings followed the pattern of the narrative dramatic arc. Moreover, while EEG amplitude could not predict group-averaged engagement ratings, other features comprising dynamic intersubject correlation (dISC), including certain frequency bands, dynamic functional connectivity patterns and graph features were able to achieve this. Furthermore, neural activity in the last two phases of the dramatic arc significantly predicted engagement patterns. This study is the first to explore the cognitive processes behind the dramatic arc and its phases. By demonstrating how neural activity predicts self-reported engagement, which itself aligns with the narrative structure, this study provides insights on the interrelationships between narrative structure, neural responses, and viewer engagement.

Information flow and dynamic functional connectivity during electroconvulsive therapy in patients with depression

BACKGROUND

Electroconvulsive therapy is effectively used for treatment-resistant depression; however, its neural mechanism is largely unknown. Resting-state functional magnetic resonance imaging is promising for monitoring outcomes of electroconvulsive therapy for depression. This study aimed to explore the imaging correlates of the electroconvulsive therapy effects on depression using Granger causality analysis and dynamic functional connectivity analyses.

METHODS

We performed advanced analyses of resting-state functional magnetic resonance imaging data at the beginning and intermediate stages and end of the therapeutic course to identify neural markers that reflect or predict the therapeutic effects of electroconvulsive therapy on depression.

RESULTS

We demonstrated that information flow between the functional networks analyzed by Granger causality changes during electroconvulsive therapy, and this change was correlated with the therapeutic outcome. Information flow and the dwell time (an index reflecting the temporal stability of functional connectivity) before electroconvulsive therapy are correlated with depressive symptoms during and after treatment.

LIMITATIONS

First, the sample size was small. A larger group is needed to confirm our findings. Second, the influence of concomitant pharmacotherapy on our results was not fully addressed, although we expected it to be minimal because only minor changes in pharmacotherapy occurred during electroconvulsive therapy. Third, different scanners were used the groups, although the acquisition parameters were the same; a direct comparison between patient and healthy participant data was not possible. Thus, we presented the data of the healthy participants separately from that of the patients as a reference.

CONCLUSIONS

These results show the specific properties of functional brain connectivity.

The link between static and dynamic brain functional network connectivity and genetic risk of Alzheimer's disease

Apolipoprotein E (APOE) polymorphic alleles are genetic factors associated with Alzheimer's disease (AD) risk. Although previous studies have explored the link between AD genetic risk and static functional network connectivity (sFNC), to the best of our knowledge, no previous studies have evaluated the association between dynamic FNC (dFNC) and AD genetic risk. Here, we examined the link between sFNC, dFNC, and AD genetic risk with a data-driven approach. We used rs-fMRI, demographic, and APOE data from cognitively normal individuals (N = 886) between 42 and 95 years of age (mean = 70 years). We separated individuals into low, moderate, and high-risk groups. Using Pearson correlation, we calculated sFNC across seven brain networks. We also calculated dFNC with a sliding window and Pearson correlation. The dFNC windows were partitioned into three distinct states with k-means clustering. Next, we calculated the proportion of time each subject spent in each state, called occupancy rate or OCR and frequency of visits. We compared both sFNC and dFNC features across individuals with different genetic risks and found that both sFNC and dFNC are related to AD genetic risk. We found that higher AD risk reduces within-visual sensory network (VSN) sFNC and that individuals with higher AD risk spend more time in a state with lower within-VSN dFNC. We also found that AD genetic risk affects whole-brain sFNC and dFNC in women but not men. In conclusion, we presented novel insights into the links between sFNC, dFNC, and AD genetic risk.

Predictive signature of static and dynamic functional connectivity for ECT clinical outcomes

Introduction: Electroconvulsive therapy (ECT) remains one of the most effective approaches for treatment-resistant depressive episodes, despite the potential cognitive impairment associated with this treatment. As a potent stimulator of neuroplasticity, ECT might normalize aberrant depression-related brain function via the brain's reconstruction by forming new neural connections. Multiple lines of evidence have demonstrated that functional connectivity (FC) changes are reliable indicators of antidepressant efficacy and cognitive changes from static and dynamic perspectives. However, no previous studies have directly ascertained whether and how different aspects of FC provide complementary information in terms of neuroimaging-based prediction of clinical outcomes. Methods: In this study, we implemented a fully automated independent component analysis framework to an ECT dataset with subjects (n = 50, age = 65.54 ± 8.92) randomized to three treatment amplitudes (600, 700, or 800 milliamperes [mA]). We extracted the static functional network connectivity (sFNC) and dynamic FNC (dFNC) features and employed a partial least square regression to build predictive models for antidepressant outcomes and cognitive changes. Results: We found that both antidepressant outcomes and memory changes can be robustly predicted by the changes in sFNC (permutation test p < 5.0 × 10^-3). More interestingly, by adding dFNC information, the model achieved higher accuracy for predicting changes in the Hamilton Depression Rating Scale 24-item (HDRS₂₄, t = 9.6434, p = 1.5 × 10^-21). The predictive maps of clinical outcomes show a weakly negative correlation, indicating that the ECT-induced antidepressant outcomes and cognitive changes might be associated with different functional brain neuroplasticity. Discussion: The overall results reveal that dynamic FC is not redundant but reflects mechanisms of ECT that cannot be captured by its static counterpart, especially for the prediction of antidepressant efficacy. Tracking the predictive signatures of static and dynamic FC will help maximize antidepressant outcomes and cognitive safety with individualized ECT dosing.

Manifold Learning for fMRI time-varying FC

Whole-brain functional connectivity ( FC ) measured with functional MRI (fMRI) evolve over time in meaningful ways at temporal scales going from years (e.g., development) to seconds (e.g., within-scan time-varying FC ( tvFC )). Yet, our ability to explore tvFC is severely constrained by its large dimensionality (several thousands). To overcome this difficulty, researchers seek to generate low dimensional representations (e.g., 2D and 3D scatter plots) expected to retain its most informative aspects (e.g., relationships to behavior, disease progression). Limited prior empirical work suggests that manifold learning techniques ( MLTs )-namely those seeking to infer a low dimensional non-linear surface (i.e., the manifold) where most of the data lies-are good candidates for accomplishing this task. Here we explore this possibility in detail. First, we discuss why one should expect tv FC data to lie on a low dimensional manifold. Second, we estimate what is the intrinsic dimension (i.e., minimum number of latent dimensions; ID ) of tvFC data manifolds. Third, we describe the inner workings of three state-of-the-art MLTs : Laplacian Eigenmaps ( LE ), T-distributed Stochastic Neighbor Embedding ( T-SNE ), and Uniform Manifold Approximation and Projection ( UMAP ). For each method, we empirically evaluate its ability to generate neuro-biologically meaningful representations of tvFC data, as well as their robustness against hyper-parameter selection. Our results show that tvFC data has an ID that ranges between 4 and 26, and that ID varies significantly between rest and task states. We also show how all three methods can effectively capture subject identity and task being performed: UMAP and T-SNE can capture these two levels of detail concurrently, but L E could only capture one at a time. We observed substantial variability in embedding quality across MLTs , and within- MLT as a function of hyper-parameter selection. To help alleviate this issue, we provide heuristics that can inform future studies. Finally, we also demonstrate the importance of feature normalization when combining data across subjects and the role that temporal autocorrelation plays in the application of MLTs to tvFC data. Overall, we conclude that while MLTs can be useful to generate summary views of labeled tvFC data, their application to unlabeled data such as resting-state remains challenging.

Mapping intrinsic functional network topological architecture in major depression disorder after electroconvulsive therapy

Disrupted topological organization of functional brain networks has been well documented in major depressive disorder (MDD). However, there is no report about how electroconvulsive therapy (ECT), a rapid way for depression remission, affects whole-brain functional network topological architecture to improve clinical symptoms in individuals with MDD. In this study, resting-state functional magnetic resonance imaging (rs-fMRI) data were collected for twenty-four MDD patients before and after receiving ECT and 25 gender-, age- and education-matched healthy controls (HC). The functional brain network for each subject was mapped using Brainnetome Atlas and graph-theory was applied to measure topological properties for both binary and weighted network. The results showed that ECT can significantly increase shortest path length and decrease global efficiency in MDD patients. In addition, significant alterations in nodal degree, nodal efficiency as well as between nodal functional connectivity strength were found in MDD patients after ECT. The network nodes showing changed degree, efficiency and connectivity were primarily distributed in default mode network (DMN), fronto-parietal network (FPN), and limbic system. Our findings demonstrates that ECT improves depressive symptoms by reorganizing disrupted network topological architecture in MDD patients and highlights the important role of functional reorganization of DMN, FPN, and limbic network contributing to depression remission.