A confounder controlled machine learning approach: Group analysis and classification of schizophrenia and Alzheimer's disease using resting-state functional network connectivity.

Resting-state functional magnetic resonance imaging (rs-fMRI) has increasingly been used to study both Alzheimer's disease (AD) and schizophrenia (SZ). While most rs-fMRI studies being conducted in AD and SZ compare patients to healthy controls, it is also of interest to directly compare AD and SZ patients with each other to identify potential biomarkers shared between the disorders. However, comparing patient groups collected in different studies can be challenging due to potential confounds, such as differences in the patient's age, scan protocols, etc. In this study, we compared and contrasted resting-state functional network connectivity (rs-FNC) of 162 patients with AD and late mild cognitive impairment (LMCI), 181 schizophrenia patients, and 315 cognitively normal (CN) subjects. We used confounder-controlled rs-FNC and applied machine learning algorithms (including support vector machine, logistic regression, random forest, and k-nearest neighbor) and deep learning models (i.e., fully-connected neural networks) to classify subjects in binary and three-class categories according to their diagnosis labels (e.g., AD, SZ, and CN). Our statistical analysis revealed that FNC between the following network pairs is stronger in AD compared to SZ: subcortical-cerebellum, subcortical-cognitive control, cognitive control-cerebellum, and visual-sensory motor networks. On the other hand, FNC is stronger in SZ than AD for the following network pairs: subcortical-visual, subcortical-auditory, subcortical-sensory motor, cerebellum-visual, sensory motor-cognitive control, and within the cerebellum networks. Furthermore, we observed that while AD and SZ disorders each have unique FNC abnormalities, they also share some common functional abnormalities that can be due to similar neurobiological mechanisms or genetic factors contributing to these disorders' development. Moreover, we achieved an accuracy of 85% in classifying subjects into AD and SZ where default mode, visual, and subcortical networks contributed the most to the classification and accuracy of 68% in classifying subjects into AD, SZ, and CN with the subcortical domain appearing as the most contributing features to the three-way classification. Finally, our findings indicated that for all classification tasks, except AD vs. SZ, males are more predictable than females.

Topological state-space estimation of functional human brain networks.

We introduce an innovative, data-driven topological data analysis (TDA) technique for estimating the state spaces of dynamically changing functional human brain networks at rest. Our method utilizes the Wasserstein distance to measure topological differences, enabling the clustering of brain networks into distinct topological states. This technique outperforms the commonly used k-means clustering in identifying brain network state spaces by effectively incorporating the temporal dynamics of the data without the need for explicit model specification. We further investigate the genetic underpinnings of these topological features using a twin study design, examining the heritability of such state changes. Our findings suggest that the topology of brain networks, particularly in their dynamic state changes, may hold significant hidden genetic information.

The risk of cannabis use disorder is mediated by altered brain connectivity: A chronnectome study.

The brain mechanisms underlying the risk of cannabis use disorder (CUD) are poorly understood. Several studies have reported changes in functional connectivity (FC) in CUD, although none have focused on the study of time-varying patterns of FC. To fill this important gap of knowledge, 39 individuals at risk for CUD and 55 controls, stratified by their score on a self-screening questionnaire for cannabis-related problems (CUDIT-R), underwent resting-state functional magnetic resonance imaging. Dynamic functional connectivity (dFNC) was estimated using independent component analysis, sliding-time window correlations, cluster states and meta-state indices of global dynamics and were compared among groups. At-risk individuals stayed longer in a cluster state with higher within and reduced between network dFNC for the subcortical, sensory-motor, visual, cognitive-control and default-mode networks, relative to controls. More globally, at-risk individuals had a greater number of meta-states and transitions between them and a longer state span and total distance between meta-states in the state space. Our findings suggest that the risk of CUD is associated with an increased dynamic fluidity and dynamic range of FC. This may result in altered stability and engagement of the brain networks, which can ultimately translate into altered cortical and subcortical function conveying CUD risk. Identifying these changes in brain function can pave the way for early pharmacological and neurostimulation treatment of CUD, as much as they could facilitate the stratification of high-risk individuals.

Cross-cohort replicable resting-state functional connectivity in predicting symptoms and cognition of schizophrenia.

Schizophrenia (SZ) is a debilitating mental illness characterized by adolescence or early adulthood onset of psychosis, positive and negative symptoms, as well as cognitive impairments. Despite a plethora of studies leveraging functional connectivity (FC) from functional magnetic resonance imaging (fMRI) to predict symptoms and cognitive impairments of SZ, the findings have exhibited great heterogeneity. We aimed to identify congruous and replicable connectivity patterns capable of predicting positive and negative symptoms as well as cognitive impairments in SZ. Predictable functional connections (FCs) were identified by employing an individualized prediction model, whose replicability was further evaluated across three independent cohorts (BSNIP, SZ = 174; COBRE, SZ = 100; FBIRN, SZ = 161). Across cohorts, we observed that altered FCs in frontal-temporal-cingulate-thalamic network were replicable in prediction of positive symptoms, while sensorimotor network was predictive of negative symptoms. Temporal-parahippocampal network was consistently identified to be associated with reduced cognitive function. These replicable 23 FCs effectively distinguished SZ from healthy controls (HC) across three cohorts (82.7%, 90.2%, and 86.1%). Furthermore, models built using these replicable FCs showed comparable accuracies to those built using the whole-brain features in predicting symptoms/cognition of SZ across the three cohorts (r = .17-.33, p < .05). Overall, our findings provide new insights into the neural underpinnings of SZ symptoms/cognition and offer potential targets for further research and possible clinical interventions.

Oral ketamine effects on dynamics of functional network connectivity in patients treated for chronic suicidality.

The underlying brain mechanisms of ketamine in treating chronic suicidality and the characteristics of patients who will benefit from ketamine treatment remain unclear. To address these gaps, we investigated temporal variations of brain functional synchronisation in patients with suicidality treated with ketamine in a 6-week open-label oral ketamine trial. The trial's primary endpoint was the Beck Scale for Suicide Ideation (BSS). Patients who experienced greater than 50% improvement in BSS scores or had a BSS score less than 6 at the post-treatment and follow-up (10 weeks) visits were considered responders and persistent responders, respectively. The reoccurring and transient connectivity pattern (termed brain state) from 29 patients (45.6 years ± 14.5, 15 females) were investigated by dynamic functional connectivity analysis of resting-state functional MRI at the baseline, post-treatment, and follow-up. Post-treatment patients showed significantly more (FDR-Q = 0.03) transitions among whole brain states than at baseline. We also observed increased dwelling time (FDR-Q = 0.04) and frequency (FDR-Q = 0.04) of highly synchronised brain state at follow-up, which were significantly correlated with BSS scores (both FDR-Q = 0.008). At baseline, persistent responders had higher fractions (FDR-Q = 0.03, Cohen's d = 1.39) of a cognitive control network state with high connectivities than non-responders. These findings suggested that ketamine enhanced brain changes among different synchronisation patterns and enabled high synchronisation patterns in the long term, providing a possible biological pathway for its suicide-prevention effects. Moreover, differences in cognitive control states at baseline may be used for precise ketamine treatment planning.

MRI morphometry of the anterior and posterior cerebellar vermis and its relationship to sensorimotor and cognitive functions in children.

INTRODUCTION: The human cerebellum emerges as a posterior brain structure integrating neural networks for sensorimotor, cognitive, and emotional processing across the lifespan. Developmental studies of the cerebellar anatomy and function are scant. We examine age-dependent MRI morphometry of the anterior cerebellar vermis, lobules I-V and posterior neocortical lobules VI-VII and their relationship to sensorimotor and cognitive functions. METHODS: Typically developing children (TDC; n=38; age 9-15) and healthy adults (HAC; n=31; 18-40) participated in high-resolution MRI. Rigorous anatomically informed morphometry of the vermis lobules I-V and VI-VII and total brain volume (TBV) employed manual segmentation computer-assisted FreeSurfer Image Analysis Program [http://surfer.nmr.mgh.harvard.edu]. The neuropsychological scores (WASI-II) were normalized and related to volumes of anterior, posterior vermis, and TBV. RESULTS: TBVs were age independent. Volumes of I-V and VI-VII were significantly reduced in TDC. The ratio of VI-VII to I-V (∼60%) was stable across age-groups; I-V correlated with visual-spatial-motor skills; VI-VII with verbal, visual-abstract and FSIQ. CONCLUSIONS: In TDC neither anterior I-V nor posterior VI-VII vermis attained adult volumes. The "inverted U" developmental trajectory of gray matter peaking in adolescence does not explain this finding. The hypothesis of protracted development of oligodendrocyte/myelination is suggested as a contributor to TDC's lower cerebellar vermis volumes.

Uncovering Effects of Schizophrenia upon a Maximally Significant, Minimally Complex Subset of Default Mode Network Connectivity Features.

A common analysis approach for resting state functional magnetic resonance imaging (rs-fMRI) dynamic functional network connectivity (dFNC) data involves clustering windowed correlation time-series and assigning time windows to clusters (i.e., states) that can be quantified to summarize aspects of the dFNC dynamics. However, those methods can be dominated by a select few features and obscure key dynamics related to less dominant features. This study presents an iterative feature learning approach to identify a maximally significant and minimally complex subset of dFNC features within the default mode network (DMN) in schizophrenia (SZ). Utilizing dFNC data from individuals with SZ and healthy controls (HC), our approach uncovers a subset of features that has a greater number of dFNC states with disorder-related dynamics than is found when all features are present in the clustering. We find that anterior cingulate cortex/posterior cingulate cortex (ACC/PCC) interactions are consistently related to SZ across the most significant iterations of the feature learning analysis and that individuals with SZ tend to spend more time in states with greater intra-ACC anticorrelation and almost no time in a state of high intra-ACC correlation that HCs periodically enter. Our findings highlight the need for nuanced analyses to reveal disorder-related dynamics and advance our understanding of neuropsychiatric disorders.

Identifying EEG Biomarkers of Depression with Novel Explainable Deep Learning Architectures.

Deep learning methods are increasingly being applied to raw electroencephalogram (EEG) data. However, if these models are to be used in clinical or research contexts, methods to explain them must be developed, and if these models are to be used in research contexts, methods for combining explanations across large numbers of models must be developed to counteract the inherent randomness of existing training approaches. Model visualization-based explainability methods for EEG involve structuring a model architecture such that its extracted features can be characterized and have the potential to offer highly useful insights into the patterns that they uncover. Nevertheless, model visualization-based explainability methods have been underexplored within the context of multichannel EEG, and methods to combine their explanations across folds have not yet been developed. In this study, we present two novel convolutional neural network-based architectures and apply them for automated major depressive disorder diagnosis. Our models obtain slightly lower classification performance than a baseline architecture. However, across 50 training folds, they find that individuals with MDD exhibit higher ß power, potentially higher δ power, and higher brain-wide correlation that is most strongly represented within the right hemisphere. This study provides multiple key insights into MDD and represents a significant step forward for the domain of explainable deep learning applied to raw EEG. We hope that it will inspire future efforts that will eventually enable the development of explainable EEG deep learning models that can contribute both to clinical care and novel medical research discoveries.

Parallel Multilink Group Joint ICA: Fusion of 3D Structural and 4D Functional Data Across Multiple Resting fMRI Networks.

Multimodal neuroimaging research plays a pivotal role in understanding the complexities of the human brain and its disorders. Independent component analysis (ICA) has emerged as a widely used and powerful tool for disentangling mixed independent sources, particularly in the analysis of functional magnetic resonance imaging (fMRI) data. This paper extends the use of ICA as a unifying framework for multimodal fusion, introducing a novel approach termed parallel multilink group joint ICA (pmg-jICA). The method allows for the fusion of gray matter maps from structural MRI (sMRI) data to multiple fMRI intrinsic networks, addressing the limitations of previous models. The effectiveness of pmg-jICA is demonstrated through its application to an Alzheimer's dataset, yielding linked structure-function outputs for 53 brain networks. Our approach leverages the complementary information from various imaging modalities, providing a unique perspective on brain alterations in Alzheimer's disease. The pmg-jICA identifies several components with significant differences between HC and AD groups including thalamus, caudate, putamen with in the subcortical (SC) domain, insula, parahippocampal gyrus within the cognitive control (CC) domain, and the lingual gyrus within the visual (VS) domain, providing localized insights into the links between AD and specific brain regions. In addition, because we link across multiple brain networks, we can also compute functional network connectivity (FNC) from spatial maps and subject loadings, providing a detailed exploration of the relationships between different brain regions and allowing us to visualize spatial patterns and loading parameters in sMRI along with intrinsic networks and FNC from the fMRI data. In essence, developed approach combines concepts from joint ICA and group ICA to provide a rich set of output characterizing data-driven links between covarying gray matter networks, and a (potentially large number of) resting fMRI networks allowing further study in the context of structure/function links. We demonstrate the utility of the approach by highlighting key structure/function disruptions in Alzheimer's individuals.

Neuroimaging alterations and relapse in early-stage psychosis.

BACKGROUND: Recent reports have indicated that symptom exacerbation after a period of improvement, referred to as relapse, in early-stage psychosis could result in brain changes and poor disease outcomes. We hypothesized that substantial neuroimaging alterations may exist among patients who experience relapse in early-stage psychosis. METHODS: We studied patients with psychosis within 2 years after the first psychotic event and healthy controls. We divided patients into 2 groups, namely those who did not experience relapse between disease onset and the magnetic resonance imaging (MRI) scan (no-relapse group) and those who did experience relapse between these 2 timings (relapse group). We analyzed 3003 functional connectivity estimates between 78 regions of interest (ROIs) derived from resting-state functional MRI data by adjusting for demographic and clinical confounding factors. RESULTS: We studied 85 patients, incuding 54 in the relapse group and 31 in the no-relapse group, along with 94 healthy controls. We observed significant differences in 47 functional connectivity estimates between the relapse and control groups after multiple comparison corrections, whereas no differences were found between the no-relapse and control groups. Most of these pathological signatures (64%) involved the thalamus. The Jonckheere-Terpstra test indicated that all 47 functional connectivity changes had a significant cross-group progression from controls to patients in the no-relapse group to patients in the relapse group. LIMITATIONS: Longitudinal studies are needed to further validate the involvement and pathological importance of the thalamus in relapse. CONCLUSION: We observed pathological differences in neuronal connectivity associated with relapse in early-stage psychosis, which are more specifically associated with the thalamus. Our study implies the importance of considering neurobiological mechanisms associated with relapse in the trajectory of psychotic disorders.

A method for estimating dynamic functional network connectivity gradients (dFNG) from ICA captures smooth inter-network modulation.

Dynamic functional network connectivity (dFNC) analysis is a widely used approach for studying brain function and offering insight into how brain networks evolve over time. Typically, dFNC studies utilized fixed spatial maps and evaluate transient changes in coupling among time courses estimated from independent component analysis (ICA). This manuscript presents a complementary approach that relaxes this assumption by spatially reordering the components dynamically at each timepoint to optimize for a smooth gradient in the FNC (i.e., a smooth gradient among ICA connectivity values). Several methods are presented to summarize dynamic FNC gradients (dFNGs) over time, starting with static FNC gradients (sFNGs), then exploring the reordering properties as well as the dynamics of the gradients themselves. We then apply this approach to a dataset of schizophrenia (SZ) patients and healthy controls (HC). Functional dysconnectivity between different brain regions has been reported in schizophrenia, yet the neural mechanisms behind it remain elusive. Using resting state fMRI and ICA on a dataset consisting of 151 schizophrenia patients and 160 age and gender-matched healthy controls, we extracted 53 intrinsic connectivity networks (ICNs) for each subject using a fully automated spatially constrained ICA approach. We develop several summaries of our functional network connectivity gradient analysis, both in a static sense, computed as the Pearson correlation coefficient between full time series, and a dynamic sense, computed using a sliding window approach followed by reordering based on the computed gradient, and evaluate group differences. Static connectivity analysis revealed significantly stronger connectivity between subcortical (SC), auditory (AUD) and visual (VIS) networks in patients, as well as hypoconnectivity in sensorimotor (SM) network relative to controls. sFNG analysis highlighted distinctive clustering patterns in patients and HCs along cognitive control (CC)/ default mode network (DMN), SC/ AUD/ SM/ cerebellar (CB), and VIS gradients. Furthermore, we observed significant differences in the sFNGs between groups in SC and CB domains. dFNG analysis suggested that SZ patients spend significantly more time in a SC/ CB state based on the first gradient, while HCs favor the DMN state. For the second gradient, however, patients exhibited significantly higher activity in CB/ VIS domains, contrasting with HCs' DMN engagement. The gradient synchrony analysis conveyed more shifts between SM/ SC networks and transmodal CC/ DMN networks in patients. In addition, the dFNG coupling revealed distinct connectivity patterns between SC, SM and CB centroids in SZ patients compared to HCs. To recap, our results advance our understanding of brain network modulation by examining smooth connectivity trajectories. This provides a more complete spatiotemporal summary of the data, contributing to the growing body of current literature regarding the functional dysconnectivity in schizophrenia patients. By employing dFNG, we highlight a new perspective to capture large scale fluctuations across the brain while maintaining the convenience of brain networks and low dimensional summary measures.

Towards a multimodal neuroimaging-based risk score for mild cognitive impairment by combining clinical studies with a large (N>37000) population-based study.

Alzheimer's disease (AD) is the most common form of age-related dementia, leading to a decline in memory, reasoning, and social skills. While numerous studies have investigated the genetic risk factors associated with AD, less attention has been given to identifying a brain imaging-based measure of AD risk. This study introduces a novel approach to assess mild cognitive impairment MCI, as a stage before AD, risk using neuroimaging data, referred to as a brain-wide risk score (BRS), which incorporates multimodal brain imaging. To begin, we first categorized participants from the Open Access Series of Imaging Studies (OASIS)-3 cohort into two groups: controls (CN) and individuals with MCI. Next, we computed structure and functional imaging features from all the OASIS data as well as all the UK Biobank data. For resting functional magnetic resonance imaging (fMRI) data, we computed functional network connectivity (FNC) matrices using fully automated spatially constrained independent component analysis. For structural MRI data we computed gray matter (GM) segmentation maps. We then evaluated the similarity between each participant's neuroimaging features from the UK Biobank and the difference in the average of those features between CN individuals and those with MCI, which we refer to as the brain-wide risk score (BRS). Both GM and FNC features were utilized in determining the BRS. We first evaluated the differences in the distribution of the BRS for CN vs MCI within the OASIS-3 (using OASIS-3 as the reference group). Next, we evaluated the BRS in the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort (using OASIS-3 as the reference group), showing that the BRS can differentiate MCI from CN in an independent data set. Subsequently, using the sMRI BRS, we identified 10 distinct subgroups and similarly, we identified another set of 10 subgroups using the FNC BRS. For sMRI and FNC we observed results that mutually validate each other, with certain aspects being complementary. For the unimodal analysis, sMRI provides greater differentiation between MCI and CN individuals than the fMRI data, consistent with prior work. Additionally, by utilizing a multimodal BRS approach, which combines both GM and FNC assessments, we identified two groups of subjects using the multimodal BRS scores. One group exhibits high MCI risk with both negative GM and FNC BRS, while the other shows low MCI risk with both positive GM and FNC BRS. Moreover, in the UKBB we have 46 participants diagnosed with AD showed FNC and GM patterns similar to those in high-risk groups, defined in both unimodal and multimodal BRS. Finally, to ensure the reproducibility of our findings, we conducted a validation analysis using the ADNI as an additional reference dataset and repeated the above analysis. The results were consistently replicated across different reference groups, highlighting the potential of FNC and sMRI-based BRS in early Alzheimer's detection.

Developmentally sensitive multispectral cortical connectivity profiles serving visual selective attention.

Throughout childhood and adolescence, the brain undergoes significant structural and functional changes that contribute to the maturation of multiple cognitive domains, including selective attention. Selective attention is crucial for healthy executive functioning and while key brain regions serving selective attention have been identified, their age-related changes in neural oscillatory dynamics and connectivity remain largely unknown. We examined the developmental sensitivity of selective attention circuitry in 91 typically developing youth aged 6 - 13 years old. Participants completed a number-based Simon task while undergoing magnetoencephalography (MEG) and the resulting data were preprocessed and transformed into the time-frequency domain. Significant oscillatory brain responses were imaged using a beamforming approach, and task-related peak voxels in the occipital, parietal, and cerebellar cortices were used as seeds for subsequent whole-brain connectivity analyses in the alpha and gamma range. Our key findings revealed developmentally sensitive connectivity profiles in multiple regions crucial for selective attention, including the temporoparietal junction (alpha) and prefrontal cortex (gamma). Overall, these findings suggest that brain regions serving selective attention are highly sensitive to developmental changes during the pubertal transition period.

Analysis of Longitudinal Change Patterns in Developing Brain Using Functional and Structural Magnetic Resonance Imaging via Multimodal Fusion.

Functional and structural magnetic resonance imaging (fMRI and sMRI) are complementary approaches that can be used to study longitudinal brain changes in adolescents. Each individual modality offers distinct insights into the brain. Each individual modality may overlook crucial aspects of brain analysis. By combining them, we can uncover hidden brain connections and gain a more comprehensive understanding. In previous work, we identified multivariate patterns of change in whole-brain function during adolescence. In this work, we focus on linking functional change patterns (FCPs) to brain structure. We introduce two approaches and applied them to data from the Adolescent Brain and Cognitive Development (ABCD) dataset. First, we evaluate voxelwise sMRI-FCP coupling to identify structural patterns linked to our previously identified FCPs. Our approach revealed multiple interesting patterns in functional network connectivity (FNC) and gray matter volume (GMV) data that were linked to subject level variation. FCP components 2 and 4 exhibit extensive associations between their loadings and voxel-wise GMV data. Secondly, we leveraged a symmetric multimodal fusion technique called multiset canonical correlation analysis (mCCA) + joint independent component analysis (jICA). Using this approach, we identify structured FCPs such as one showing increased connectivity between visual and sensorimotor domains and decreased connectivity between sensorimotor and cognitive control domains, linked to structural change patterns (SCPs) including alterations in the bilateral sensorimotor cortex. Interestingly, females exhibit stronger coupling between brain functional and structural changes than males, highlighting sex-related differences. The combined results from both asymmetric and symmetric multimodal fusion methods underscore the intricate sex-specific nuances in neural dynamics. By utilizing two complementary multimodal approaches, our study enhances our understanding of the dynamic nature of brain connectivity and structure during the adolescent period, shedding light on the nuanced processes underlying adolescent brain development.

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.

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).

Dynamic functional connectivity in anorexia nervosa: alterations in states of low connectivity and state transitions.

BACKGROUND: The onset of anorexia nervosa (AN) frequently occurs during adolescence and is associated with preoccupation with body weight and shape and extreme underweight. Altered resting state functional connectivity in the brain has been described in individuals with AN, but only from a static perspective. The current study investigated the temporal dynamics of functional connectivity in adolescents with AN and how it relates to clinical features. METHOD: 99 female patients acutely ill with AN and 99 pairwise age-matched female healthy control (HC) participants were included in the study. Using resting-state functional MRI data and an established sliding-window analytic approach, we identified dynamic resting-state functional connectivity states and extracted dynamic indices such as dwell time (the duration spent in a state), fraction time (the proportion of the total time occupied by a state), and number of transitions (number of switches) from one state to another, to test for group differences. RESULTS: Individuals with AN had relatively reduced fraction time in a mildly connected state with pronounced connectivity within the default mode network (DMN) and an overall reduced number of transitions between states. CONCLUSIONS: These findings revealed by a dynamic, but not static analytic approach might hint towards a more "rigid" connectivity, a phenomenon commonly observed in internalizing mental disorders, and in AN possibly related to a reduction in energetic costs as a result of nutritional deprivation.

SMART (Splitting-Merging Assisted Reliable) Independent Component Analysis for Extracting Accurate Brain Functional Networks.

Functional networks (FNs) hold significant promise in understanding brain function. Independent component analysis (ICA) has been applied in estimating FNs from functional magnetic resonance imaging (fMRI). However, determining an optimal model order for ICA remains challenging, leading to criticism about the reliability of FN estimation. Here, we propose a SMART (splitting-merging assisted reliable) ICA method that automatically extracts reliable FNs by clustering independent components (ICs) obtained from multi-model-order ICA using a simplified graph while providing linkages among FNs deduced from different-model orders. We extend SMART ICA to multi-subject fMRI analysis, validating its effectiveness using simulated and real fMRI data. Based on simulated data, the method accurately estimates both group-common and group-unique components and demonstrates robustness to parameters. Using two age-matched cohorts of resting fMRI data comprising 1,950 healthy subjects, the resulting reliable group-level FNs are greatly similar between the two cohorts, and interestingly the subject-specific FNs show progressive changes while age increases. Furthermore, both small-scale and large-scale brain FN templates are provided as benchmarks for future studies. Taken together, SMART ICA can automatically obtain reliable FNs in analyzing multi-subject fMRI data, while also providing linkages between different FNs.

Multiview hyperedge-aware hypergraph embedding learning for multisite, multiatlas fMRI based functional connectivity network analysis.

Recently, functional magnetic resonance imaging (fMRI) based functional connectivity network (FCN) analysis via graph convolutional networks (GCNs) has shown promise for automated diagnosis of brain diseases by regarding the FCNs as irregular graph-structured data. However, multiview information and site influences of the FCNs in a multisite, multiatlas fMRI scenario have been understudied. In this paper, we propose a Class-consistency and Site-independence Multiview Hyperedge-Aware HyperGraph Embedding Learning (CcSi-MHAHGEL) framework to integrate FCNs constructed on multiple brain atlases in a multisite fMRI study. Specifically, for each subject, we first model brain network as a hypergraph for every brain atlas to characterize high-order relations among multiple vertexes, and then introduce a multiview hyperedge-aware hypergraph convolutional network (HGCN) to extract a multiatlas-based FCN embedding where hyperedge weights are adaptively learned rather than employing the fixed weights precalculated in traditional HGCNs. In addition, we formulate two modules to jointly learn the multiatlas-based FCN embeddings by considering the between-subject associations across classes and sites, respectively, i.e., a class-consistency module to encourage both compactness within every class and separation between classes for promoting discrimination in the embedding space, and a site-independence module to minimize the site dependence of the embeddings for mitigating undesired site influences due to differences in scanning platforms and/or protocols at multiple sites. Finally, the multiatlas-based FCN embeddings are fed into a few fully connected layers followed by the soft-max classifier for diagnosis decision. Extensive experiments on the ABIDE demonstrate the effectiveness of our method for autism spectrum disorder (ASD) identification. Furthermore, our method is interpretable by revealing ASD-relevant brain regions that are biologically significant.

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.