Group ICA for identifying biomarkers in schizophrenia: 'Adaptive' networks via spatially constrained ICA show more sensitivity to group differences than spatio-temporal regression

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.

Longitudinal resting-state network connectivity changes in electroconvulsive therapy patients compared to healthy controls

OBJECTIVE

Electroconvulsive therapy (ECT) is effective for major depressive episodes. Understanding of underlying mechanisms has been increased by examining changes of brain connectivity but studies often do not correct for test-retest variability in healthy controls (HC). In this study, we investigated changes in resting-state networks after ECT in a multicenter study.

METHODS

Functional resting-state magnetic resonance imaging data, acquired before start and within one week after ECT, from 90 depressed patients were analyzed, as well as longitudinal data of 24 HC. Group-information guided independent component analysis (GIG-ICA) was used to spatially restrict decomposition to twelve canonical resting-state networks. Selected networks of interest were the default mode network (DMN), salience network (SN), and left and right frontoparietal network (LFPN, and RFPN). Whole-brain voxel-wise analyses were used to assess group differences at baseline, group by time interactions, and correlations with treatment effectiveness. In addition, between-network connectivity and within-network strengths were computed.

RESULTS

Within-network strength of the DMN was lower at baseline in ECT patients which increased after ECT compared to HC, after which no differences were detected. At baseline, ECT patients showed lower whole-brain voxel-wise DMN connectivity in the precuneus. Increase of within-network strength of the LFPN was correlated with treatment effectiveness. We did not find whole-brain voxel-wise or between-network changes.

CONCLUSION

DMN within-network connectivity normalized after ECT. Within-network increase of the LFPN in ECT patients was correlated with higher treatment effectiveness. In contrast to earlier studies, we found no whole-brain voxel-wise changes, which highlights the necessity to account for test-retest effects.

Revisiting Functional Dysconnectivity: a Review of Three Model Frameworks in Schizophrenia

PURPOSE OF REVIEW

Over the last decade, evidence suggests that a combination of behavioral and neuroimaging findings can help illuminate changes in functional dysconnectivity in schizophrenia. We review the recent connectivity literature considering several vital models, considering connectivity findings, and relationships with clinical symptoms. We reviewed resting state fMRI studies from 2017 to 2023. We summarized the role of two sets of brain networks (cerebello-thalamo-cortical (CTCC) and the triple network set) across three hypothesized models of schizophrenia etiology (neurodevelopmental, vulnerability-stress, and neurotransmitter hypotheses).

RECENT FINDINGS

The neurotransmitter and neurodevelopmental models best explained CTCC-subcortical dysfunction, which was consistently connected to symptom severity and motor symptoms. Triple network dysconnectivity was linked to deficits in executive functioning, and the salience network (SN)-default mode network dysconnectivity was tied to disordered thought and attentional deficits. This paper links behavioral symptoms of schizophrenia (symptom severity, motor, executive functioning, and attentional deficits) to various hypothesized mechanisms.

Multi-study evaluation of neuroimaging-based prediction of medication class in mood disorders

Clinicians often face a dilemma in diagnosing bipolar disorder patients with complex symptoms who spend more time in a depressive state than a manic state. The current gold standard for such diagnosis, the Diagnostic and Statistical Manual (DSM), is not objectively grounded in pathophysiology. In such complex cases, relying solely on the DSM may result in misdiagnosis as major depressive disorder (MDD). A biologically-based classification algorithm that can accurately predict treatment response may help patients suffering from mood disorders. Here we used an algorithm to do so using neuroimaging data. We used the neuromark framework to learn a kernel function for support vector machine (SVM) on multiple feature subspaces. The neuromark framework achieves up to 95.45% accuracy, 0.90 sensitivity, and 0.92 specificity in predicting antidepressant (AD) vs. mood stabilizer (MS) response in patients. We incorporated two additional datasets to evaluate the generalizability of our approach. The trained algorithm achieved up to 89% accuracy, 0.88 sensitivity, and 0.89 specificity in predicting the DSM-based diagnosis on these datasets. We also translated the model to distinguish responders to treatment from nonresponders with up to 70% accuracy. This approach reveals multiple salient biomarkers of medication-class of response within mood disorders.

A Novel Method for Multi-subject fMRI Data Analysis: Independent Component Analysis with Clustering Embedded (ICA-CE)

The analysis of multi-subject functional magnetic resonance imaging (fMRI) data and the extraction of accurate brain functional networks (FNs) are of great importance. However, traditional independent component analysis (ICA) methods perform analysis on multi-subject fMRI data under the condition of known or assumed classes of subjects, which may decrease its ability to extract accurate individual brain FNs. Although a previous method named clusterwise ICA (C-ICA) clusters subjects and obtains shared FNs in group-level for each class, its clustering performance on complex data is not ideal. To address the issues, we propose a novel method called independent component analysis with clustering embedded (ICA-CE) that can achieve both the estimation of individual FNs and the clustering of subjects in an unsupervised or semi-supervised manner. Using the simulated data with different properties, ICA-CE achieved better clustering performance than group ICA followed by K-means and C-ICA, and the mean accuracy of extracted individual FNs obtained by ICA-CE was greater than 90%. Using the task-related fMRI data from Human Connectome Project (HCP), our method also achieved higher clustering accuracy, while extracting task-related class-specific FNs. In summary, ICA-CE is effective in estimating accurate brain FNs while achieving the clustering of multiple subjects.Clinical Relevance- Our method is promising in estimating accurate brain functional networks for patients with brain disorders and outputting related class label for each subject.

Altered large-scale dynamic connectivity patterns in Alzheimer's disease and mild cognitive impairment patients: A machine learning study

Alzheimer's disease (AD) is a common neurodegeneration disease associated with substantial disruptions in the brain network. However, most studies investigated static resting-state functional connections, while the alteration of dynamic functional connectivity in AD remains largely unknown. This study used group independent component analysis and the sliding-window method to estimate the subject-specific dynamic connectivity states in 1704 individuals from three data sets. Informative inherent states were identified by the multivariate pattern classification method, and classifiers were built to distinguish ADs from normal controls (NCs) and to classify mild cognitive impairment (MCI) patients with informative inherent states similar to ADs or not. In addition, MCI subgroups with heterogeneous functional states were examined in the context of different cognition decline trajectories. Five informative states were identified by feature selection, mainly involving functional connectivity belonging to the default mode network and working memory network. The classifiers discriminating AD and NC achieved the mean area under the receiver operating characteristic curve of 0.87 with leave-one-site-out cross-validation. Alterations in connectivity strength, fluctuation, and inter-synchronization were found in AD and MCIs. Moreover, individuals with MCI were clustered into two subgroups, which had different degrees of atrophy and different trajectories of cognition decline progression. The present study uncovered the alteration of dynamic functional connectivity in AD and highlighted that the dynamic states could be powerful features to discriminate patients from NCs. Furthermore, it demonstrated that these states help to identify MCIs with faster cognition decline and might contribute to the early prevention of AD.

A Scalable Approach to Independent Vector Analysis by Shared Subspace Separation for Multi-Subject fMRI Analysis

Joint blind source separation (JBSS) has wide applications in modeling latent structures across multiple related datasets. However, JBSS is computationally prohibitive with high-dimensional data, limiting the number of datasets that can be included in a tractable analysis. Furthermore, JBSS may not be effective if the data's true latent dimensionality is not adequately modeled, where severe overparameterization may lead to poor separation and time performance. In this paper, we propose a scalable JBSS method by modeling and separating the "shared" subspace from the data. The shared subspace is defined as the subset of latent sources that exists across all datasets, represented by groups of sources that collectively form a low-rank structure. Our method first provides the efficient initialization of the independent vector analysis (IVA) with a multivariate Gaussian source prior (IVA-G) specifically designed to estimate the shared sources. Estimated sources are then evaluated regarding whether they are shared, upon which further JBSS is applied separately to the shared and non-shared sources. This provides an effective means to reduce the dimensionality of the problem, improving analyses with larger numbers of datasets. We apply our method to resting-state fMRI datasets, demonstrating that our method can achieve an excellent estimation performance with significantly reduced computational costs.

Brain subnetworks most sensitive to alterations of functional connectivity in Schizophrenia: a data-driven approach

Functional connectivity (FC) of the brain changes in various brain disorders. Its complexity, however, makes it difficult to obtain a systematic understanding of these alterations, especially when they are found individually and through hypothesis-based methods. It would be easier if the variety of brain connectivity alterations is extracted through data-driven approaches and expressed as variation modules (subnetworks). In the present study, we modified a blind approach to determine inter-group brain variations at the network level and applied it specifically to schizophrenia (SZ) disorder. The analysis is based on the application of independent component analysis (ICA) over the subject's dimension of the FC matrices, obtained from resting-state functional magnetic resonance imaging (rs-fMRI). The dataset included 27 SZ people and 27 completely matched healthy controls (HC). This hypothesis-free approach led to the finding of three brain subnetworks significantly discriminating SZ from HC. The area associated with these subnetworks mostly covers regions in visual, ventral attention, and somatomotor areas, which are in line with previous studies. Moreover, from the graph perspective, significant differences were observed between SZ and HC for these subnetworks, while there was no significant difference when the same parameters (path length, network strength, global/local efficiency, and clustering coefficient) across the same limited data were calculated for the whole brain network. The increased sensitivity of those subnetworks to SZ-induced alterations of connectivity suggested whether an individual scoring method based on their connectivity values can be applied to classify subjects. A simple scoring classifier was then suggested based on two of these subnetworks and resulted in acceptable sensitivity and specificity with an area under the ROC curve of 77.5%. The third subnetwork was found to be a less specific building block (module) for describing SZ alterations. It projected a wider range of inter-individual variations and, therefore, had a lower chance to be considered as a SZ biomarker. These findings confirmed that investigating brain variations from a modular viewpoint can help to find subnetworks that are more sensitive to SZ-induced alterations. Altogether, our study results illustrated the developed method's ability to systematically find brain alterations caused by SZ disorder from a network perspective.

Multi-model order spatially constrained ICA reveals highly replicable group differences and consistent predictive results from resting data: A large N fMRI schizophrenia study

Brain functional networks identified from resting functional magnetic resonance imaging (fMRI) data have the potential to reveal biomarkers for brain disorders, but studies of complex mental illnesses such as schizophrenia (SZ) often yield mixed results across replication studies. This is likely due in part to the complexity of the disorder, the short data acquisition time, and the limited ability of the approaches for brain imaging data mining. Therefore, the use of analytic approaches which can both capture individual variability while offering comparability across analyses is highly preferred. Fully blind data-driven approaches such as independent component analysis (ICA) are hard to compare across studies, and approaches that use fixed atlas-based regions can have limited sensitivity to individual sensitivity. By contrast, spatially constrained ICA (scICA) provides a hybrid, fully automated solution that can incorporate spatial network priors while also adapting to new subjects. However, scICA has thus far only been used with a single spatial scale (ICA dimensionality, i.e., ICA model order). In this work, we present an approach using multi-objective optimization scICA with reference algorithm (MOO-ICAR) to extract subject-specific intrinsic connectivity networks (ICNs) from fMRI data at multiple spatial scales, which also enables us to study interactions across spatial scales. We evaluate this approach using a large N (N > 1,600) study of schizophrenia divided into separate validation and replication sets. A multi-scale ICN template was estimated and labeled, then used as input into scICA which was computed on an individual subject level. We then performed a subsequent analysis of multiscale functional network connectivity (msFNC) to evaluate the patient data, including group differences and classification. Results showed highly consistent group differences in msFNC in regions including cerebellum, thalamus, and motor/auditory networks. Importantly, multiple msFNC pairs linking different spatial scales were implicated. The classification model built on the msFNC features obtained up to 85% F1 score, 83% precision, and 88% recall, indicating the strength of the proposed framework in detecting group differences between schizophrenia and the control group. Finally, we evaluated the relationship of the identified patterns to positive symptoms and found consistent results across datasets. The results verified the robustness of our framework in evaluating brain functional connectivity of schizophrenia at multiple spatial scales, implicated consistent and replicable brain networks, and highlighted a promising approach for leveraging resting fMRI data for brain biomarker development.

A deep residual model for characterization of 5D spatiotemporal network dynamics reveals widespread spatiodynamic changes in schizophrenia

Schizophrenia is a severe brain disorder with serious symptoms including delusions, disorganized speech, and hallucinations that can have a long-term detrimental impact on different aspects of a patient's life. It is still unclear what the main cause of schizophrenia is, but a combination of altered brain connectivity and structure may play a role. Neuroimaging data has been useful in characterizing schizophrenia, but there has been very little work focused on voxel-wise changes in multiple brain networks over time, despite evidence that functional networks exhibit complex spatiotemporal changes over time within individual subjects. Recent studies have primarily focused on static (average) features of functional data or on temporal variations between fixed networks; however, such approaches are not able to capture multiple overlapping networks which change at the voxel level. In this work, we employ a deep residual convolutional neural network (CNN) model to extract 53 different spatiotemporal networks each of which captures dynamism within various domains including subcortical, cerebellar, visual, sensori-motor, auditory, cognitive control, and default mode. We apply this approach to study spatiotemporal brain dynamism at the voxel level within multiple functional networks extracted from a large functional magnetic resonance imaging (fMRI) dataset of individuals with schizophrenia (N = 708) and controls (N = 510). Our analysis reveals widespread group level differences across multiple networks and spatiotemporal features including voxel-wise variability, magnitude, and temporal functional network connectivity in widespread regions expected to be impacted by the disorder. We compare with typical average spatial amplitude and show highly structured and neuroanatomically relevant results are missed if one does not consider the voxel-wise spatial dynamics. Importantly, our approach can summarize static, temporal dynamic, spatial dynamic, and spatiotemporal dynamics features, thus proving a powerful approach to unify and compare these various perspectives. In sum, we show the proposed approach highlights the importance of accounting for both temporal and spatial dynamism in whole brain neuroimaging data generally, shows a high-level of sensitivity to schizophrenia highlighting global but spatially unique dynamics showing group differences, and may be especially important in studies focused on the development of brain-based biomarkers.

IABC: A Toolbox for Intelligent Analysis of Brain Connectivity

Brain functional networks and connectivity have played an important role in exploring brain function for understanding the brain and disclosing the mechanisms of brain disorders. Independent component analysis (ICA) is one of the most widely applied data-driven methods to extract brain functional networks/connectivity. However, it is hard to guarantee the reliability of networks/connectivity due to the randomness of component order and the difficulty in selecting an optimal component number in ICA. To facilitate the analysis of brain functional networks and connectivity using ICA, we developed a MATLAB toolbox called Intelligent Analysis of Brain Connectivity (IABC). IABC incorporates our previously proposed group information guided independent component analysis (GIG-ICA), NeuroMark, and splitting-merging assisted reliable ICA (SMART ICA) methods, which can estimate reliable individual-subject neuroimaging measures for further analysis. After user inputs functional magnetic resonance imaging (fMRI) data of multiple subjects that are regularly organized (e.g., in Brain Imaging Data Structure (BIDS)) and clicks a few buttons to set parameters, IABC automatically outputs brain functional networks, their related time courses, and functional network connectivity of each subject. All these neuroimaging measures are promising for providing clues in understanding brain function and differentiating brain disorders.

A method for capturing dynamic spectral coupling in resting fMRI reveals domain-specific patterns in schizophrenia

Introduction

Resting-state functional magnetic resonance imaging (rs-fMRI) is a powerful tool for assessing functional brain connectivity. Recent studies have focused on shorter-term connectivity and dynamics in the resting state. However, most of the prior work evaluates changes in time-series correlations. In this study, we propose a framework that focuses on time-resolved spectral coupling (assessed via the correlation between power spectra of the windowed time courses) among different brain circuits determined via independent component analysis (ICA).

Methods

Motivated by earlier work suggesting significant spectral differences in people with schizophrenia, we developed an approach to evaluate time-resolved spectral coupling (trSC). To do this, we first calculated the correlation between the power spectra of windowed time-courses pairs of brain components. Then, we subgrouped each correlation map into four subgroups based on the connectivity strength utilizing quartiles and clustering techniques. Lastly, we examined clinical group differences by regression analysis for each averaged count and average cluster size matrices in each quartile. We evaluated the method by applying it to resting-state data collected from 151 (114 males, 37 females) people with schizophrenia (SZ) and 163 (117 males, 46 females) healthy controls (HC).

Results

Our proposed approach enables us to observe the change of connectivity strength within each quartile for different subgroups. People with schizophrenia showed highly modularized and significant differences in multiple network domains, whereas males and females showed less modular differences. Both cell count and average cluster size analysis for subgroups indicate a higher connectivity rate in the fourth quartile for the visual network in the control group. This indicates increased trSC in visual networks in the controls. In other words, this shows that the visual networks in people with schizophrenia have less mutually consistent spectra. It is also the case that the visual networks are less spectrally correlated on short timescales with networks of all other functional domains.

Conclusions

The results of this study reveal significant differences in the degree to which spectral power profiles are coupled over time. Importantly, there are significant but distinct differences both between males and females and between people with schizophrenia and controls. We observed a more significant coupling rate in the visual network for the healthy controls and males in the upper quartile. Fluctuations over time are complex, and focusing on only time-resolved coupling among time-courses is likely to miss important information. Also, people with schizophrenia are known to have impairments in visual processing but the underlying reasons for the impairment are still unknown. Therefore, the trSC approach can be a useful tool to explore the reasons for the impairments.

Tracking functional network connectivity dynamics in the elderly

Introduction

Functional magnetic resonance imaging (fMRI) has shown that aging disturbs healthy brain organization and functional connectivity. However, how this age-induced alteration impacts dynamic brain function interaction has not yet been fully investigated. Dynamic function network connectivity (DFNC) analysis can produce a brain representation based on the time-varying network connectivity changes, which can be further used to study the brain aging mechanism for people at different age stages.

Method

This presented investigation examined the dynamic functional connectivity representation and its relationship with brain age for people at an elderly stage as well as in early adulthood. Specifically, the resting-state fMRI data from the University of North Carolina cohort of 34 young adults and 28 elderly participants were fed into a DFNC analysis pipeline. This DFNC pipeline forms an integrated dynamic functional connectivity (FC) analysis framework, which consists of brain functional network parcellation, dynamic FC feature extraction, and FC dynamics examination.

Results

The statistical analysis demonstrates that extensive dynamic connection changes in the elderly concerning the transient brain state and the method of functional interaction in the brain. In addition, various machine learning algorithms have been developed to verify the ability of dynamic FC features to distinguish the age stage. The fraction time of DFNC states has the highest performance, which can achieve a classification accuracy of over 88% by a decision tree.

Discussion

The results proved there are dynamic FC alterations in the elderly, and the alteration was found to be correlated with mnemonic discrimination ability and could have an impact on the balance of functional integration and segregation.

Multi-Subject Analysis for Brain Developmental Patterns Discovery via Tensor Decomposition of MEG Data

Identification of informative signatures from electrophysiological signals is important for understanding brain developmental patterns, where techniques such as magnetoencephalography (MEG) are particularly useful. However, less attention has been given to fully utilizing the multidimensional nature of MEG data for extracting components that describe these patterns. Tensor factorizations of MEG yield components that encapsulate the data's multidimensional nature, providing parsimonious models identifying latent brain patterns for meaningful summarization of neural processes. To address the need for meaningful MEG signatures for studies of pediatric cohorts, we propose a tensor-based approach for extracting developmental signatures of multi-subject MEG data. We employ the canonical polyadic (CP) decomposition for estimating latent spatiotemporal components of the data, and use these components for group level statistical inference. Using CP decomposition along with hierarchical clustering, we were able to extract typical early and late latency event-related field (ERF) components that were discriminative of high and low performance groups ([Formula: see text]) and significantly correlated with major cognitive domains such as attention, episodic memory, executive function, and language comprehension. We demonstrate that tensor-based group level statistical inference of MEG can produce signatures descriptive of the multidimensional MEG data. Furthermore, these features can be used to study group differences in brain patterns and cognitive function of healthy children. We provide an effective tool that may be useful for assessing child developmental status and brain function directly from electrophysiological measurements and facilitate the prospective assessment of cognitive processes.

Dynamic functional connectivity patterns associated with dementia risk

Background

This study assesses the relationships between dynamic functional network connectivity (DFNC) and dementia risk.

Methods

DFNC of the default mode (DMN), salience (SN), and executive control networks was assessed in 127 cognitively unimpaired older adults. Stepwise regressions were performed with dementia risk and protective factors and biomarkers as predictors of DFNC.

Results

Associations were found between times spent in (i) a “weakly connected” state and lower self-reported engagement in early- and mid-life cognitive activity and higher LDL cholesterol; (ii) a “SN-negatively connected” state and higher blood pressure, higher depression score, and lower body mass index (BMI); (iii) a “strongly connected” state and higher self-reported engagement in early-life cognitive activity, Preclinical Alzheimer’s cognitive composite-5 score, and BMI; and (iv) a “DMN-negatively connected” state and higher self-reported engagement in early- and mid-life stimulating activities and lower LDL cholesterol and blood pressure. The lower number of state transitions was associated with lower brain perfusion.

Conclusion

DFNC states are differentially associated with dementia risk and could underlie reserve.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13195-022-01006-7.

Ability of an altered functional coupling between resting-state networks to predict behavioral outcomes in subcortical ischemic stroke: A longitudinal study

Stroke can be viewed as an acute disruption of an individual's connectome caused by a focal or widespread loss of blood flow. Although individuals exhibit connectivity changes in multiple functional networks after stroke, the neural mechanisms that underlie the longitudinal reorganization of the connectivity patterns are still unclear. The study aimed to determine whether brain network connectivity patterns after stroke can predict longitudinal behavioral outcomes. Nineteen patients with stroke with subcortical lesions underwent two sessions of resting-state functional magnetic resonance imaging scanning at a 1-month interval. By independent component analysis, the functional connectivity within and between multiple brain networks (including the default mode network, the dorsal attention network, the limbic network, the visual network, and the frontoparietal network) was disrupted after stroke and partial recovery at the second time point. Additionally, regression analyses revealed that the connectivity between the limbic and dorsal attention networks at the first time point showed sufficient reliability in predicting the clinical scores (Fugl-Meyer Assessment and Neurological Deficit Scores) at the second time point. The overall findings suggest that functional coupling between the dorsal attention and limbic networks after stroke can be regarded as a biomarker to predict longitudinal clinical outcomes in motor function and the degree of neurological functional deficit. Overall, the present study provided a novel opportunity to improve prognostic ability after subcortical strokes.

Hierarchical Individual Naturalistic Functional Brain Networks with Group Consistency uncovered by a Two-Stage NAS-Volumetric Sparse DBN Framework

The functional magnetic resonance imaging under naturalistic paradigm (NfMRI) showed great advantages in identifying complex and interactive functional brain networks due to its dynamics and multimodal information. In recent years, various deep learning models, such as deep convolutional autoencoder (DCAE), deep belief network (DBN) and volumetric sparse deep belief network (vsDBN), can obtain hierarchical functional brain networks (FBN) and temporal features from fMRI data. Among them, the vsDBN model revealed a good capability in identifying hierarchical FBNs by modelling fMRI volume images. However, due to the high dimensionality of fMRI volumes and the diverse training parameters of deep learning methods, especially the network architecture that is the most critical parameter for uncovering the hierarchical organization of human brain function, researchers still face challenges in designing an appropriate deep learning framework with automatic network architecture optimization to model volumetric NfMRI. In addition, most of the existing deep learning models ignore the group-wise consistency and inter-subject variation properties embedded in NfMRI volumes. To solve these problems, we proposed a two-stage neural architecture search and vs DBN model (two-stage NAS-vsDBN model) to identify the hierarchical human brain spatio-temporal features possessing both group-consistency and individual-uniqueness under naturalistic condition. Moreover, our model defined reliable network structure for modelling volumetric NfMRI data via NAS framework, and the group-level and individual-level FBNs and associated temporal features exhibited great consistency. In general, our method well identified the hierarchical temporal and spatial features of the brain function and revealed the crucial properties of neural processes under natural viewing condition.Significance StatementIn this paper, we proposed and applied a novel analytical strategy - a two-stage NAS-vsDBN model to identify both group-level and individual-level spatio-temporal features at multi-scales from volumetric NfMRI data. The proposed PSO-based NAS framework can find optimal neural structure for both group-wise and individual-level vs-DBN models. Furthermore, with well-established correspondence between two stages of vsDBN models, our model can effectively detect group-level FBNs that reveal the consistency in neural processes across subjects and individual-level FBNs that maintain the subject specific variability, verifying the inherent property of brain function under naturalistic condition.

Reduced functional connectivity supports statistical learning of temporally distributed regularities

Statistical learning is a powerful ability that extracts regularities from our environment and makes predictions about future events. Using functional magnetic resonance imaging, we aimed to probe how a wide range of brain areas are intertwined to support statistical learning, characterising its architecture in the whole-brain functional connectivity (FC). Participants performed a statistical learning task of temporally distributed regularities. We used refined behavioural learning scores to associate individuals' learning performances with the FC changed by statistical learning. As a result, the learning performance was mediated by the activation strength in the lateral occipital cortex, angular gyrus, precuneus, anterior cingulate cortex, and superior frontal gyrus. Through a group independent component analysis, activations of the superior frontal network showed the largest correlation with the statistical learning performances. Seed-to-voxel whole-brain and seed-to-ROI FC analyses revealed that the FC between the superior frontal gyrus and the salience, language, and dorsal attention networks were reduced during statistical learning. We suggest that the weakened functional connections between the superior frontal gyrus and brain regions involved in top-down control processes serve a pivotal role in statistical learning, supporting better processing of novel information such as the extraction of new patterns from the environment.

Widespread cortical functional disconnection in gliomas: an individual network mapping approach

Assessment of impaired/preserved cortical regions in brain tumours is typically performed via intraoperative direct brain stimulation of eloquent areas or task-based functional MRI. One main limitation is that they overlook distal brain regions or networks that could be functionally impaired by the tumour.

This study aims: 1) to investigate the impact of brain tumours on the cortical synchronization of brain networks measured with resting-state functional magnetic resonance imaging (resting-state networks) both near the lesion and remotely; 2) to test whether potential changes in resting state networks correlate with cognitive status.

The sample included twenty-four glioma patients (mean age 58.1 ± 16.4y) with different pathological staging. We developed a new method for single subject localization of resting state networks abnormalities. First, we derived the spatial pattern of the main resting state networks by means of the group guided independent component analysis. This was informed by a high-resolution resting state networks template derived from an independent sample of healthy controls. Second, we developed a spatial similarity index to measure differences in network topography and strength between healthy controls and individual brain tumour patients. Next, we investigated the spatial relationship between altered networks and tumour location. Finally, multivariate analyses related cognitive scores across multiple cognitive domains (attention, language, memory, decision making) with patterns of multi-network abnormality.

We found that brain gliomas cause broad alterations of resting state networks topography that occurred mainly in structurally normal regions outside the tumour and oedema region. Cortical regions near the tumour often showed normal synchronization. Finally, multi-network abnormalities predicted attention deficits.

Overall, we present a novel method for the functional localization of resting state networks abnormalities in individual glioma patients. These abnormalities partially explain cognitive disabilities and shall be carefully navigated during surgery.

Individual prediction of trauma-focused psychotherapy response in youth with posttraumatic stress disorder using resting-state functional connectivity

•

ML and rs-fMRI have shown promise in predicting treatment-response in adults with PTSD.

•

Currently, no biomarkers for treatment-response are available in youth with PTSD.

•

FC between the FPN and SMN was stronger in treatment non-responders on the group-level.

•

A network within the bilateral STG predicted response for individual youth with 76% accuracy.

•

Future studies should test generalizability of these findings and test if larger cohorts increase accuracy.

Randomized controlled trials have shown efficacy of trauma-focused psychotherapies in youth with posttraumatic stress disorder (PTSD). However, response varies considerably among individuals. Currently, no biomarkers are available to assist clinicians in identifying youth who are most likely to benefit from treatment. In this study, we investigated whether resting-state functional magnetic resonance imaging (rs-fMRI) could distinguish between responders and non-responders on the group- and individual patient level. Pre-treatment rs-fMRI was recorded in 40 youth (ages 8–17 years) with (partial) PTSD before trauma-focused psychotherapy. Change in symptom severity from pre- to post-treatment was assessed using the Clinician-Administered PTSD scale for Children and Adolescents to divide participants into responders (≥30% symptom reduction) and non-responders. Functional networks were identified using meta-independent component analysis. Group-differences within- and between-network connectivity between responders and non-responders were tested using permutation testing. Individual predictions were made using multivariate, cross-validated support vector machine classification. A network centered on the bilateral superior temporal gyrus predicted treatment response for individual patients with 76% accuracy (p_FWE = 0.02, 87% sensitivity, 65% specificity, area-under-receiver-operator-curve of 0.82). Functional connectivity between the frontoparietal and sensorimotor network was significantly stronger in non-responders (t = 5.35, p_FWE = 0.01) on the group-level. Within-network connectivity was not significantly different between groups. This study provides proof-of-concept evidence for the feasibility to predict trauma-focused psychotherapy response in youth with PTSD at an individual-level. Future studies are required to test if larger cohorts could increase accuracy and to test further generalizability of the prediction models.

Multivariate patterns of brain-behavior associations across the adult lifespan

The nature of brain-behavior covariations with increasing age is poorly understood. In the current study, we used a multivariate approach to investigate the covariation between behavioral-health variables and brain features across adulthood. We recruited healthy adults aged 20–73 years-old (29 younger, mean age = 25.6 years; 30 older, mean age = 62.5 years), and collected structural and functional MRI (s/fMRI) during a resting-state and three tasks. From the sMRI, we extracted cortical thickness and subcortical volumes; from the fMRI, we extracted activation peaks and functional network connectivity (FNC) for each task. We conducted canonical correlation analyses between behavioral-health variables and the sMRI, or the fMRI variables, across all participants. We found significant covariations for both types of neuroimaging phenotypes (ps = 0.0004) across all individuals, with cognitive capacity and age being the largest opposite contributors. We further identified different variables contributing to the models across phenotypes and age groups. Particularly, we found behavior was associated with different neuroimaging patterns between the younger and older groups. Higher cognitive capacity was supported by activation and FNC within the executive networks in the younger adults, while it was supported by the visual networks’ FNC in the older adults. This study highlights how the brain-behavior covariations vary across adulthood and provides further support that cognitive performance relies on regional recruitment that differs between older and younger individuals.

BrainGNN: Interpretable Brain Graph Neural Network for fMRI Analysis

Understanding which brain regions are related to a specific neurological disorder or cognitive stimuli has been an important area of neuroimaging research. We propose BrainGNN, a graph neural network (GNN) framework to analyze functional magnetic resonance images (fMRI) and discover neurological biomarkers. Considering the special property of brain graphs, we design novel ROI-aware graph convolutional (Ra-GConv) layers that leverage the topological and functional information of fMRI. Motivated by the need for transparency in medical image analysis, our BrainGNN contains ROI-selection pooling layers (R-pool) that highlight salient ROIs (nodes in the graph), so that we can infer which ROIs are important for prediction. Furthermore, we propose regularization terms-unit loss, topK pooling (TPK) loss and group-level consistency (GLC) loss-on pooling results to encourage reasonable ROI-selection and provide flexibility to encourage either fully individual- or patterns that agree with group-level data. We apply the BrainGNN framework on two independent fMRI datasets: an Autism Spectrum Disorder (ASD) fMRI dataset and data from the Human Connectome Project (HCP) 900 Subject Release. We investigate different choices of the hyper-parameters and show that BrainGNN outperforms the alternative fMRI image analysis methods in terms of four different evaluation metrics. The obtained community clustering and salient ROI detection results show a high correspondence with the previous neuroimaging-derived evidence of biomarkers for ASD and specific task states decoded for HCP. Our code is available at https://github.com/xxlya/BrainGNN_Pytorch.

Developmental Changes in Dynamic Functional Connectivity From Childhood Into Adolescence

The longitudinal study of typical neurodevelopment is key for understanding deviations due to specific factors, such as psychopathology. However, research utilizing repeated measurements remains scarce. Resting-state functional magnetic resonance imaging (MRI) studies have traditionally examined connectivity as 'static' during the measurement period. In contrast, dynamic approaches offer a more comprehensive representation of functional connectivity by allowing for different connectivity configurations (time varying connectivity) throughout the scanning session. Our objective was to characterize the longitudinal developmental changes in dynamic functional connectivity in a population-based pediatric sample. Resting-state MRI data were acquired at the ages of 10 (range 8-to-12, n = 3,327) and 14 (range 13-to-15, n = 2,404) years old using a single, study-dedicated 3 Tesla scanner. A fully-automated spatially constrained group-independent component analysis (ICA) was applied to decompose multi-subject resting-state data into functionally homogeneous regions. Dynamic functional network connectivity (FNC) between all ICA time courses were computed using a tapered sliding window approach. We used a k-means algorithm to cluster the resulting dynamic FNC windows from each scan session into five dynamic states. We examined age and sex associations using linear mixed-effects models. First, independent from the dynamic states, we found a general increase in the temporal variability of the connections between intrinsic connectivity networks with increasing age. Second, when examining the clusters of dynamic FNC windows, we observed that the time spent in less modularized states, with low intra- and inter-network connectivity, decreased with age. Third, the number of transitions between states also decreased with age. Finally, compared to boys, girls showed a more mature pattern of dynamic brain connectivity, indicated by more time spent in a highly modularized state, less time spent in specific states that are frequently observed at a younger age, and a lower number of transitions between states. This longitudinal population-based study demonstrates age-related maturation in dynamic intrinsic neural activity from childhood into adolescence and offers a meaningful baseline for comparison with deviations from typical development. Given that several behavioral and cognitive processes also show marked changes through childhood and adolescence, dynamic functional connectivity should also be explored as a potential neurobiological determinant of such changes.

Dynamic connectivity predicts acute motor impairment and recovery post-stroke

Thorough assessment of cerebral dysfunction after acute lesions is paramount to optimize predicting clinical outcomes. We here built random forest classifier-based prediction models of acute motor impairment and recovery post-stroke. Predictions relied on structural and resting-state fMRI data from 54 stroke patients scanned within the first days of symptom onset. Functional connectivity was estimated via static and dynamic approaches. Motor performance was phenotyped in the acute phase and 6 months later. A model based on the time spent in specific dynamic connectivity configurations achieved the best discrimination between patients with and without motor impairments (out-of-sample area under the curve, 95% confidence interval: 0.67 ± 0.01). In contrast, patients with moderate-to-severe impairments could be differentiated from patients with mild deficits using a model based on the variability of dynamic connectivity (0.83 ± 0.01). Here, the variability of the connectivity between ipsilesional sensorimotor cortex and putamen discriminated the most between patients. Finally, motor recovery was best predicted by the time spent in specific connectivity configurations (0.89 ± 0.01) in combination with the initial impairment. Here, better recovery was linked to a shorter time spent in a functionally integrated configuration. Dynamic connectivity-derived parameters constitute potent predictors of acute impairment and recovery, which, in the future, might inform personalized therapy regimens to promote stroke recovery.

Fusing multimodal neuroimaging data with a variational autoencoder

Neuroimaging studies often collect multimodal data. These modalities contain both shared and mutually exclusive information about the brain. This work aims to find a scalable and interpretable method to fuse the information of multiple neuroimaging modalities into a lower-dimensional latent space using a variational autoencoder (VAE). To assess whether the encoder-decoder pair retains meaningful information, this work evaluates the representations using a schizophrenia classification task. The linear classifier, trained on the representations obtained through dimensionality reduction, achieves an area under the curve of the receiver operating characteristic (ROC-AUC) of 0.8609. Thus, training on a multimodal dataset with functional brain networks and a structural magnetic resonance imaging (sMRI) scan, leads to dimensionality reduction that retains meaningful information. The proposed dimensionality reduction outperforms both early and late fusion principal component analysis on the classification task.Clinical relevance - This work examines the interplay between neuroimaging modalities and their relation to mental disorders. This allows for more complex and rigorous analysis of multimodal neuroimaging data throughout clinical settings.

Dissemination in time and space in presymptomatic granulin mutation carriers: a GENFI spatial chronnectome study

The presymptomatic brain changes of granulin (GRN) disease, preceding by years frontotemporal dementia, has not been fully characterized. New approaches focus on the spatial chronnectome can capture both spatial network configurations and their dynamic changes over time. To investigate the spatial dynamics in 141 presymptomatic GRN mutation carriers and 282 noncarriers from the Genetic Frontotemporal dementia research Initiative cohort. We considered time-varying patterns of the default mode network, the language network, and the salience network, each summarized into 4 distinct recurring spatial configurations. Dwell time (DT) (the time each individual spends in each spatial state of each network), fractional occupacy (FO) (the total percentage of time spent by each individual in a state of a specific network) and total transition number (the total number of transitions performed by each individual in a specifict state) were considered. Correlations between DT, FO, and transition number and estimated years from expected symptom onset (EYO) and clinical performances were assessed. Presymptomatic GRN mutation carriers spent significantly more time in those spatial states characterised by greater activation of the insula and the parietal cortices, as compared to noncarriers (p < 0.05, FDR-corrected). A significant correlation between DT and FO of these spatial states and EYO was found, the longer the time spent in the spatial states, the closer the EYO. DT and FO significantly correlated with performances at tests tapping processing speed, with worse scores associated with increased spatial states' DT. Our results demonstrated that presymptomatic GRN disease presents a complex dynamic reorganization of brain connectivity. Change in both the spatial and temporal aspects of brain network connectivity could provide a unique glimpse into brain function and potentially allowing a more sophisticated evaluation of the earliest disease changes and the understanding of possible mechanisms in GRN disease.

Evidence of shared and distinct functional and structural brain signatures in schizophrenia and autism spectrum disorder

Schizophrenia (SZ) and autism spectrum disorder (ASD) share considerable clinical features and intertwined historical roots. It is greatly needed to explore their similarities and differences in pathophysiologic mechanisms. We assembled a large sample size of neuroimaging data (about 600 SZ patients, 1000 ASD patients, and 1700 healthy controls) to study the shared and unique brain abnormality of the two illnesses. We analyzed multi-scale brain functional connectivity among functional networks and brain regions, intra-network connectivity, and cerebral gray matter density and volume. Both SZ and ASD showed lower functional integration within default mode and sensorimotor domains, but increased interaction between cognitive control and default mode domains. The shared abnormalties in intra-network connectivity involved default mode, sensorimotor, and cognitive control networks. Reduced gray matter volume and density in the occipital gyrus and cerebellum were observed in both illnesses. Interestingly, ASD had overall weaker changes than SZ in the shared abnormalities. Interaction between visual and cognitive regions showed disorder-unique deficits. In summary, we provide strong neuroimaging evidence of the convergent and divergent changes in SZ and ASD that correlated with clinical features.

Identification of Methamphetamine Abstainers by Resting-State Functional Magnetic Resonance Imaging

Methamphetamine (MA) can cause brain structural and functional impairment, but there are few studies on whether this difference will sustain on MA abstainers. The purpose of this study is to investigate the correlation of brain networks in MA abstainers. In this study, 47 people detoxified for at least 14 months and 44 normal people took a resting-state functional magnetic resonance imaging (RS-fMRI) scan. A dynamic (i.e., time-varying) functional connectivity (FC) is obtained by applying sliding windows in the time courses on the independent components (ICs). The windowed correlation data for each IC were then clustered by k-means. The number of subjects in each cluster was used as a new feature for individual identification. The results show that the classifier achieved satisfactory performance (82.3% accuracy, 77.7% specificity, and 85.7% sensitivity). We find that there are significant differences in the brain networks of MA abstainers and normal people in the time domain, but the spatial differences are not obvious. Most of the altered functional connections (time-varying) are identified to be located at dorsal default mode network. These results have shown that changes in the correlation of the time domain may play an important role in identifying MA abstainers. Therefore, our findings provide valuable insights in the identification of MA and elucidate the pathological mechanism of MA from a resting-state functional integration point of view.

Construction of embedded fMRI resting-state functional connectivity networks using manifold learning

We construct embedded functional connectivity networks (FCN) from benchmark resting-state functional magnetic resonance imaging (rsfMRI) data acquired from patients with schizophrenia and healthy controls based on linear and nonlinear manifold learning algorithms, namely, Multidimensional Scaling, Isometric Feature Mapping, Diffusion Maps, Locally Linear Embedding and kernel PCA. Furthermore, based on key global graph-theoretic properties of the embedded FCN, we compare their classification potential using machine learning. We also assess the performance of two metrics that are widely used for the construction of FCN from fMRI, namely the Euclidean distance and the cross correlation metric. We show that diffusion maps with the cross correlation metric outperform the other combinations.

Altered Posterior Midline Activity in Patients with Jerky and Tremulous Functional Movement Disorders

Objective: To explore changes in resting-state networks in patients with jerky and tremulous functional movement disorders (JT-FMD). Methods: Resting-state functional magnetic resonance imaging data from seventeen patients with JT-FMD and seventeen age-, sex-, and education-matched healthy controls (HC) were investigated. Independent component analysis was used to examine the central executive network (CEN), salience network, and default mode network (DMN). Frequency distribution of network signal fluctuations and intra- and internetwork functional connectivity were investigated. Symptom severity was measured using the Clinical Global Impression-Severity scale. Beck Depression Inventory and Beck Anxiety Inventory scores were collected to measure depression and anxiety in FMD, respectively. Results: Compared with HC, patients with JT-FMD had significantly decreased power of lower range (0.01-0.10 Hz) frequency fluctuations in a precuneus and posterior cingulate cortex component of the DMN and in the dorsal attention network (DAN) component of the CEN (false discovery rate-corrected p < 0.05). No significant group differences were found for intra- and internetwork functional connectivity. In patients with JT-FMD, symptom severity was not significantly correlated with network measures. Depression scores were weakly correlated with intranetwork functional connectivity in the medial prefrontal cortex, while anxiety was not found to be related to network connectivity. Conclusions: Given the changes in the posterodorsal components of the DMN and DAN, we postulate that the JT-FMD-related functional alterations found in these regions could provide support for the concept that particularly attentional dysregulation is a fundamental disturbance in these patients.

Diagnostic power of resting-state fMRI for detection of network connectivity in Alzheimer's disease and mild cognitive impairment: A systematic review

Resting‐state fMRI (rs‐fMRI) detects functional connectivity (FC) abnormalities that occur in the brains of patients with Alzheimer's disease (AD) and mild cognitive impairment (MCI). FC of the default mode network (DMN) is commonly impaired in AD and MCI. We conducted a systematic review aimed at determining the diagnostic power of rs‐fMRI to identify FC abnormalities in the DMN of patients with AD or MCI compared with healthy controls (HCs) using machine learning (ML) methods. Multimodal support vector machine (SVM) algorithm was the commonest form of ML method utilized. Multiple kernel approach can be utilized to aid in the classification by incorporating various discriminating features, such as FC graphs based on “nodes” and “edges” together with structural MRI‐based regional cortical thickness and gray matter volume. Other multimodal features include neuropsychiatric testing scores, DTI features, and regional cerebral blood flow. Among AD patients, the posterior cingulate cortex (PCC)/Precuneus was noted to be a highly affected hub of the DMN that demonstrated overall reduced FC. Whereas reduced DMN FC between the PCC and anterior cingulate cortex (ACC) was observed in MCI patients. Evidence indicates that the nodes of the DMN can offer moderate to high diagnostic power to distinguish AD and MCI patients. Nevertheless, various concerns over the homogeneity of data based on patient selection, scanner effects, and the variable usage of classifiers and algorithms pose a challenge for ML‐based image interpretation of rs‐fMRI datasets to become a mainstream option for diagnosing AD and predicting the conversion of HC/MCI to AD.

Calhoun V, Grefkes C, Rost NS. Abnormal dynamic functional connectivity is linked to recovery after acute ischemic stroke

The aim of the current study was to explore the whole-brain dynamic functional connectivity patterns in acute ischemic stroke (AIS) patients and their relation to short and long-term stroke severity. We investigated resting-state functional MRI-based dynamic functional connectivity of 41 AIS patients two to five days after symptom onset. Re-occurring dynamic connectivity configurations were obtained using a sliding window approach and k-means clustering. We evaluated differences in dynamic patterns between three NIHSS-stroke severity defined groups (mildly, moderately, and severely affected patients). Furthermore, we built Bayesian hierarchical models to evaluate the predictive capacity of dynamic connectivity and examine the interrelation with clinical measures, such as white matter hyperintensity lesions. Finally, we established correlation analyses between dynamic connectivity and AIS severity as well as 90-day neurological recovery (ΔNIHSS). We identified three distinct dynamic connectivity configurations acutely post-stroke. More severely affected patients spent significantly more time in a configuration that was characterized by particularly strong connectivity and isolated processing of functional brain domains (three-level ANOVA: p < .05, post hoc t tests: p < .05, FDR-corrected). Configuration-specific time estimates possessed predictive capacity of stroke severity in addition to the one of clinical measures. Recovery, as indexed by the realized change of the NIHSS over time, was significantly linked to the dynamic connectivity between bilateral intraparietal lobule and left angular gyrus (Pearson's r = -.68, p = .003, FDR-corrected). Our findings demonstrate transiently increased isolated information processing in multiple functional domains in case of severe AIS. Dynamic connectivity involving default mode network components significantly correlated with recovery in the first 3 months poststroke.

A Deep Learning Model for Data-Driven Discovery of Functional Connectivity

Functional connectivity (FC) studies have demonstrated the overarching value of studying the brain and its disorders through the undirected weighted graph of functional magnetic resonance imaging (fMRI) correlation matrix. However, most of the work with the FC depends on the way the connectivity is computed, and it further depends on the manual post-hoc analysis of the FC matrices. In this work, we propose a deep learning architecture BrainGNN that learns the connectivity structure as part of learning to classify subjects. It simultaneously applies a graphical neural network to this learned graph and learns to select a sparse subset of brain regions important to the prediction task. We demonstrate that the model’s state-of-the-art classification performance on a schizophrenia fMRI dataset and demonstrate how introspection leads to disorder relevant findings. The graphs that are learned by the model exhibit strong class discrimination and the sparse subset of relevant regions are consistent with the schizophrenia literature.

ISOMAP and machine learning algorithms for the construction of embedded functional connectivity networks of anatomically separated brain regions from resting state fMRI data of patients with Schizophrenia

We construct Functional Connectivity Networks (FCN) from resting state fMRI (rsfMRI) recordings towards the classification of brain activity between healthy and schizophrenic subjects using a publicly available dataset (the COBRE dataset) of 145 subjects (74 healthy controls and 71 schizophrenic subjects). First, we match the anatomy of the brain of each individual to the Desikan-Killiany brain atlas. Then, we use the conventional approach of correlating the parcellated time series to construct FCN and ISOMAP, a nonlinear manifold learning algorithm to produce low-dimensional embeddings of the correlation matrices. For the classification analysis, we computed five key local graph-theoretic measures of the FCN and used the LASSO and Random Forest (RF) algorithms for feature selection. For the classification we used standard linear Support Vector Machines. The classification performance is tested by a double cross-validation scheme (consisting of an outer and an inner loop of "Leave one out" cross-validation (LOOCV)). The standard cross-correlation methodology produced a classification rate of 73.1%, while ISOMAP resulted in 79.3%, thus providing a simpler model with a smaller number of features as chosen from LASSO and RF, namely the participation coefficient of the right thalamus and the strength of the right lingual gyrus.

Atlas55+: Brain Functional Atlas of Resting-State Networks for Late Adulthood

Currently, several human brain functional atlases are used to define the spatial constituents of the resting-state networks (RSNs). However, the only brain atlases available are derived from samples of young adults. As brain networks are continuously reconfigured throughout life, the lack of brain atlases derived from older populations may influence RSN results in late adulthood. To address this gap, the aim of the study was to construct a reliable brain atlas derived only from older participants. We leveraged resting-state functional magnetic resonance imaging data from three cohorts of healthy older adults (total N = 563; age = 55-95 years) and a younger-adult cohort (N = 128; age = 18-35 years). We identified the major RSNs and their subdivisions across all older-adult cohorts. We demonstrated high spatial reproducibility of these RSNs with an average spatial overlap of 67%. Importantly, the RSNs derived from the older-adult cohorts were spatially different from those derived from the younger-adult cohort (P = 2.3 × 10-3). Lastly, we constructed a novel brain atlas, called Atlas55+, which includes the consensus of the major RSNs and their subdivisions across the older-adult cohorts. Thus, Atlas55+ provides a reliable age-appropriate template for RSNs in late adulthood and is publicly available. Our results confirm the need for age-appropriate functional atlases for studies investigating aging-related brain mechanisms.

Dynamic state with covarying brain activity-connectivity: On the pathophysiology of schizophrenia

The human brain is a dynamic system that incorporates the evolution of local activities and the reconfiguration of brain interactions. Reoccurring brain patterns, regarded as "brain states", have revealed new insights into the pathophysiology of brain disorders, particularly schizophrenia. However, previous studies only focus on the dynamics of either brain activity or connectivity, ignoring the temporal co-evolution between them. In this work, we propose to capture dynamic brain states with covarying activity-connectivity and probe schizophrenia-related brain abnormalities. We find that the state-based activity and connectivity show high correspondence, where strong and antagonistic connectivity is accompanied with strong low-frequency fluctuations across the whole brain while weak and sparse connectivity co-occurs with weak low-frequency fluctuations. In addition, graphical analysis shows that connectivity network efficiency is associated with the fluctuation of brain activities and such associations are different across brain states. Compared with healthy controls, schizophrenia patients spend more time in weakly-connected and -activated brain states but less time in strongly-connected and -activated brain states. schizophrenia patients also show lower efficiency in thalamic regions within the "strong" states. Interestingly, the atypical fractional occupancy of one brain state is correlated with individual attention performance. Our findings are replicated in another independent dataset and validated using different brain parcellation schemes. These converging results suggest that the brain spontaneously reconfigures with covarying activity and connectivity and such co-evolutionary property might provide meaningful information on the mechanism of brain disorders which cannot be observed by investigating either of them alone.

Dynamic Resting-State Connectivity Differences in Eyes Open Versus Eyes Closed Conditions

Introduction: Previous studies have shown significant conditional differences between eyes open, fixated at an image (EO) and eyes closed (EC) in the acquired resting-state functional magnetic resonance imaging (rs-fMRI) data. Aim: We recently showed significant functional network connectivity (FNC) differences between EO and EC across a variety of networks. In this study, we aim at further evaluating differences in dynamic FNC (dFNC) between EO and EC. Materials and Methods: Rs-fMRI were collected from adolescents aged 9-15 years old during both EO and EC conditions, and dFNC was calculated by using the independent component analysis framework. Results: We found that out of five states (clusters), state 1 was observed to be more dominant in the EO condition, whereas state 2 was observed to be more dominant in the EC condition. States 1 and 2 showed significant differences in the mean dwell time based on false discovery rate, and states 1, 2, 3, and 4 differed in the frequency of occurrences. These results are consistent with our previous study of static connectivity in suggesting that EO and EC differences not only are relatively strong but also importantly reveal that these differences vary over time, especially in one particularly transient connectivity pattern. Conclusion: Our results manifest as changes in the proportion of time spent in unique functional connectivity patterns, and they show unique transient functional connectivity patterns in a subset of identified states. Overall, our findings indicate that both static and dynamic rs-fMRI connectivity patterns are strongly impacted by basic conditional differences such as EO and EC. Impact statement Our findings not only suggest that eyes open, fixated at an image (EO) and eyes closed (EC) condition-related resting state functional magnetic resonance imaging differences are relatively strong, but they also reveal an important attribute of these conditions that these differences vary over time, especially in one particularly transient connectivity pattern. Our results manifest as changes in the proportion of time spent in unique functional connectivity patterns, and they show unique transient functional connectivity patterns in a subset of identified states. We believe there is benefit in having the EO/EC as a contrast of interest in future studies, if time allows.

Static and dynamic aspects of cerebro-cerebellar functional connectivity are associated with self-reported measures of impulsivity: A resting-state fMRI study

Human and animal brain studies bring converging evidence of a possible role for the cerebellum and the cerebro-cerebellar system in impulsivity. However, the precise nature of the relation between cerebro-cerebellar coupling and impulsivity is far from understood. Characterizing functional connectivity (FC) patterns between large-scale brain networks that mediate different forms of impulsivity, and the cerebellum may improve our understanding of this relation. Here, we analyzed static and dynamic features of cerebro-cerebellar FC using a highly sampled resting-state functional magnetic resonance imaging (rs-fMRI) dataset and tested their association with two widely used self-reports of impulsivity: the UPPS-P impulsive behavior scale and the behavioral inhibition/approach systems (BIS/BAS) in a large group of healthy subjects (N = 134, ≈ 1 hr of rs-fMRI/subject). We employed robust data-driven techniques to identify cerebral and cerebellar resting-state networks and extract descriptive summary measures of static and dynamic cerebro-cerebellar FC. We observed evidence linking BIS, BAS, sensation seeking, and lack of premeditation to the total strength and temporal variability of FC within networks connecting regions of the prefrontal cortex, precuneus, posterior cingulate cortex, basal ganglia, and thalamus with the cerebellum. Overall, our findings improve the existing knowledge of the neural correlates of impulsivity and the behavioral correlates of the cerebro-cerebellar system.

NeuroMark: An automated and adaptive ICA based pipeline to identify reproducible fMRI markers of brain disorders

•

Propose a new pipeline to link brain changes among different datasets, studies, and disorders.

•

Identify reproducible biomarkers in schizophrenia using independent data.

•

Find both common and unique brain impairments in schizophrenia and autism.

•

Reveal gradual changes from healthy controls to mild cognitive impairment to Alzheimer’s disease.

•

Obtain high classification accuracy (~90%) between bipolar disorder and major depressive disorder.

Many mental illnesses share overlapping or similar clinical symptoms, confounding the diagnosis. It is important to systematically characterize the degree to which unique and similar changing patterns are reflective of brain disorders. Increasing sharing initiatives on neuroimaging data have provided unprecedented opportunities to study brain disorders. However, it is still an open question on replicating and translating findings across studies. Standardized approaches for capturing reproducible and comparable imaging markers are greatly needed. Here, we propose a pipeline based on the priori-driven independent component analysis, NeuroMark, which is capable of estimating brain functional network measures from functional magnetic resonance imaging (fMRI) data that can be used to link brain network abnormalities among different datasets, studies, and disorders. NeuroMark automatically estimates features adaptable to each individual subject and comparable across datasets/studies/disorders by taking advantage of the reliable brain network templates extracted from 1828 healthy controls as guidance. Four studies including 2442 subjects were conducted spanning six brain disorders (schizophrenia, autism spectrum disorder, mild cognitive impairment, Alzheimer’s disease, bipolar disorder, and major depressive disorder) to evaluate validity of the proposed pipeline from different perspectives (replication of brain abnormalities, cross-study comparison, identification of subtle brain changes, and multi-disorder classification using identified biomarkers). Our results highlight that NeuroMark effectively identified replicated brain network abnormalities of schizophrenia across different datasets; revealed interesting neural clues on the overlap and specificity between autism and schizophrenia; demonstrated brain functional impairments present to varying degrees in mild cognitive impairments and Alzheimer's disease; and captured biomarkers that achieved good performance in classifying bipolar disorder and major depressive disorder.

Distinct thalamocortical network dynamics are associated with the pathophysiology of chronic low back pain

Thalamocortical dysrhythmia is a key pathology of chronic neuropathic pain, but few studies have investigated thalamocortical networks in chronic low back pain (cLBP) given its non-specific etiology and complexity. Using fMRI, we propose an analytical pipeline to identify abnormal thalamocortical network dynamics in cLBP patients and validate the findings in two independent cohorts. We first identify two reoccurring dynamic connectivity states and their associations with chronic and temporary pain. Further analyses show that cLBP patients have abnormal connectivity between the ventral lateral/posterolateral nucleus (VL/VPL) and postcentral gyrus (PoCG) and between the dorsal/ventral medial nucleus and insula in the less frequent connectivity state, and temporary pain exacerbation alters connectivity between the VL/VPL and PoCG and the default mode network in the more frequent connectivity state. These results extend current findings on thalamocortical dysfunction and dysrhythmia in chronic pain and demonstrate that cLBP pathophysiology and clinical pain intensity are associated with distinct thalamocortical network dynamics.

Thalamocortical dysrhythmia is a key pathology of chronic pain. Here, the authors propose an analytical pipeline to study dynamic fMRI brain networks and demonstrate that chronic low back pain pathophysiology and clinical pain intensity are associated with distinct thalamocortical network dynamics.

Testing Psychosis Phenotypes From Bipolar-Schizophrenia Network for Intermediate Phenotypes for Clinical Application: Biotype Characteristics and Targets

BACKGROUND

Psychiatry aspires to the molecular understanding of its disorders and, with that knowledge, to precision medicine. Research supporting such goals in the dimension of psychosis has been compromised, in part, by using phenomenology alone to estimate disease entities. To this end, we are proponents of a deep phenotyping approach in psychosis, using computational strategies to discover the most informative phenotypic fingerprint as a promising strategy to uncover mechanisms in psychosis.

METHODS

Doing this, the Bipolar-Schizophrenia Network for Intermediate Phenotypes (B-SNIP) has used biomarkers to identify distinct subtypes of psychosis with replicable biomarker characteristics. While we have presented these entities as relevant, their potential utility in clinical practice has not yet been demonstrated.

RESULTS

Here we carried out an analysis of clinical features that characterize biotypes. We found that biotypes have unique and defining clinical characteristics that could be used as initial screens in the clinical and research settings. Differences in these clinical features appear to be consistent with biotype biomarker profiles, indicating a link between biological features and clinical presentation. Clinical features associated with biotypes differ from those associated with DSM diagnoses, indicating that biotypes and DSM syndromes are not redundant and are likely to yield different treatment predictions. We highlight 3 predictions based on biotype that are derived from individual biomarker features and cannot be obtained from DSM psychosis syndromes.

CONCLUSIONS

In the future, biotypes may prove to be useful for targeting distinct molecular, circuit, cognitive, and psychosocial therapies for improved functional outcomes.

Trajectory Tracking of a Flexible Robot Manipulator by a New Optimized Fuzzy Adaptive Sliding Mode-Based Feedback Linearization Controller

This work presents a novel fuzzy adaptive sliding mode-based feedback linearization controller for trajectory tracking of a flexible robot manipulator. To reach this goal, after deriving the dynamical equations of the robot, the feedback linearization approach is utilized to change the nonlinear dynamics to a linear one and find the control law. Then, the sliding mode control strategy is implemented to design a stabilizer for trajectory tracking of the flexible robot. In order to adaptively tune the parameters of the designed controller, the gradient descent approach and the chain derivative rule are employed. Moreover, the Takagi–Sugeno–Kang fuzzy system is applied to regulate the controller gains. Finally, a multiobjective particle swarm optimization algorithm is used to find the optimum fuzzy rules. The conflicting objective functions considered as the integrals of the absolute values of the state error and the control effort should be minimized, simultaneously. The simulation results illustrate the effectiveness and capability of the introduced scenario in comparison with other methods.

Static and dynamic functional connectivity analysis of cerebrovascular reactivity: An fMRI study

BACKGROUND

Cerebrovascular reactivity (CVR) is an important aspect of brain function, and as such it is important to understand relationship between CVR and functional connectivity.

METHODS

This research studied the role of CVR, or the brain's ability to react to vasoactive stimuli on brain functional connectivity by scanning subjects with blood-oxygenation-level-dependent (BOLD) functional magnetic resonance imaging (fMRI) while they periodically inhale room air and a CO 2-enriched gas mixture. We developed a new metric to measure the effect of CVR on each intrinsic connectivity network (ICN), which contrasts to voxel-wise CVR. We also studied the changes in whole-brain connectivity patterns using both static functional network connectivity (sFNC) and dynamic FNC (dFNC).

RESULTS

We found that network connectivity is generally weaker during vascular dilation, which is supported by previous research. The dFNC analysis revealed that participants did not return to the pre-CO 2 inhalation state, suggesting that one-minute periods of room-air inhalation is not enough for the CO 2 effect to fully dissipate.

CONCLUSIONS

Cerebrovascular reactivity is one tool that the cerebrovascular system uses to ensure the constant, finely-tuned flow of oxygen to function properly. Understanding the relationship between CVR and brain dynamism can provide unique information about cerebrovascular diseases and general brain function. We observed that CVR has a wide, but consistent relationship to connectivity patterns between functional networks.

Identifying commonality and specificity across psychosis sub-groups via classification based on features from dynamic connectivity analysis

It is difficult to distinguish schizophrenia (SZ), schizoaffective disorder (SAD), and bipolar disorder with psychosis (BPP) as their clinical diagnoses rely on symptoms that overlap. In this paper, we investigate if there is biological evidence to support the symptom-based clinical categories by looking across the three disorders using dynamic connectivity measures, and provide meaningful characteristics on which brain functional connectivity measures are commonly or uniquely impaired. Large-sample functional magnetic resonance image (fMRI) datasets from 623 subjects including 238 healthy controls (HCs), 113 SZ patients, 132 SAD patients, and 140 BPP patients were analyzed. First, we computed whole-brain dynamic functional connectivity (DFC) using a sliding-window technique, and then extracted the individual connectivity states by applying our previously proposed decomposition-based DFC analysis method. Next, with the features from the dominant connectivity state, we assessed the clinical categories by performing both four-group (SZ, SAD, BPP and healthy control groups) and pair-wise classification using a support vector machine within cross-validation. Furthermore, we comprehensively summarized the shared and unique connectivity alterations among the disorders. In terms of the classification performance, our method achieved 69% in the four-group classification and >80% in the between-group classifications for the mean overall accuracy; and yielded 66% in the four-group classification and >80% in the between-group classifications for the mean balanced accuracy. Through summarizing the features that were automatically selected in the classifications, we found that among the three symptom-related disorders, their disorder-common impairments primarily included the decreased connectivity strength between thalamus and cerebellum and the increased strength between postcentral gyrus and thalamus. The disorder-unique changes included more various brain regions, mainly in the temporal and frontal gyrus. Our work demonstrates that dynamic functional connectivity provides biological evidence that both common and unique impairments exist in psychosis sub-groups.

Acute ischaemic stroke alters the brain's preference for distinct dynamic connectivity states

Acute ischaemic stroke disturbs healthy brain organization, prompting subsequent plasticity and reorganization to compensate for the loss of specialized neural tissue and function. Static resting state functional MRI studies have already furthered our understanding of cerebral reorganization by estimating stroke-induced changes in network connectivity aggregated over the duration of several minutes. In this study, we used dynamic resting state functional MRI analyses to increase temporal resolution to seconds and explore transient configurations of motor network connectivity in acute stroke. To this end, we collected resting state functional MRI data of 31 patients with acute ischaemic stroke and 17 age-matched healthy control subjects. Stroke patients presented with moderate to severe hand motor deficits. By estimating dynamic functional connectivity within a sliding window framework, we identified three distinct connectivity configurations of motor-related networks. Motor networks were organized into three regional domains, i.e. a cortical, subcortical and cerebellar domain. The dynamic connectivity patterns of stroke patients diverged from those of healthy controls depending on the severity of the initial motor impairment. Moderately affected patients (n = 18) spent significantly more time in a weakly connected configuration that was characterized by low levels of connectivity, both locally as well as between distant regions. In contrast, severely affected patients (n = 13) showed a significant preference for transitions into a spatially segregated connectivity configuration. This configuration featured particularly high levels of local connectivity within the three regional domains as well as anti-correlated connectivity between distant networks across domains. A third connectivity configuration represented an intermediate connectivity pattern compared to the preceding two, and predominantly encompassed decreased interhemispheric connectivity between cortical motor networks independent of individual deficit severity. Alterations within this third configuration thus closely resembled previously reported ones originating from static resting state functional MRI studies post-stroke. In summary, acute ischaemic stroke not only prompted changes in connectivity between distinct networks, but it also caused characteristic changes in temporal properties of large-scale network interactions depending on the severity of the individual deficit. These findings offer new vistas on the dynamic neural mechanisms underlying acute neurological symptoms, cortical reorganization and treatment effects in stroke patients.

Enhanced dynamic functional connectivity (whole-brain chronnectome) in chess experts

Multidisciplinary approaches have demonstrated that the brain is potentially modulated by the long-term acquisition and practice of specific skills. Chess playing can be considered a paradigm for shaping brain function, with complex interactions among brain networks possibly enhancing cognitive processing. Dynamic network analysis based on resting-state magnetic resonance imaging (rs-fMRI) can be useful to explore the effect of chess playing on whole-brain fluidity/dynamism (the chronnectome). Dynamic connectivity parameters of 18 professional chess players and 20 beginner chess players were evaluated applying spatial independent component analysis (sICA), sliding-time window correlation, and meta-state approaches to rs-fMRI data. Four indexes of meta-state dynamic fluidity were studied: i) the number of distinct meta-states a subject pass through, ii) the number of switches from one meta-state to another, iii) the span of the realized meta-states (the largest distance between two meta-states that subjects occupied), and iv) the total distance travelled in the state space. Professional chess players exhibited an increased dynamic fluidity, expressed as a higher number of occupied meta-states (meta-state numbers, 75.8 ± 7.9 vs 68.8 ± 12.0, p = 0.043 FDR-corrected) and changes from one meta-state to another (meta-state changes, 77.1 ± 7.3 vs 71.2 ± 11.0, p = 0.043 FDR-corrected) than beginner chess players. Furthermore, professional chess players exhibited an increased dynamic range, with increased traveling between successive meta-states (meta-state total distance, 131.7 ± 17.8 vs 108.7 ± 19.7, p = 0.0004 FDR-corrected). Chess playing may induce changes in brain activity through the modulation of the chronnectome. Future studies are warranted to evaluate if these potential effects lead to enhanced cognitive processing and if “gaming” might be used as a treatment in clinical practice.

Deep Temporal Organization of fMRI Phase Synchrony Modes Promotes Large-Scale Disconnection in Schizophrenia

Itinerant dynamics of the brain generates transient and recurrent spatiotemporal patterns in neuroimaging data. Characterizing metastable functional connectivity (FC) - particularly at rest and using functional magnetic resonance imaging (fMRI) - has shaped the field of dynamic functional connectivity (DFC). Mainstream DFC research relies on (sliding window) correlations to identify recurrent FC patterns. Recently, functional relevance of the instantaneous phase synchrony (IPS) of fMRI signals has been revealed using imaging studies and computational models. In the present paper, we identify the repertoire of whole-brain inter-network IPS states at rest. Moreover, we uncover a hierarchy in the temporal organization of IPS modes. We hypothesize that connectivity disorder in schizophrenia (SZ) is related to the (deep) temporal arrangement of large-scale IPS modes. Hence, we analyze resting-state fMRI data from 68 healthy controls (HC) and 51 SZ patients. Seven resting-state networks (and their sub-components) are identified using spatial independent component analysis. IPS is computed between subject-specific network time courses, using analytic signals. The resultant phase coupling patterns, across time and subjects, are clustered into eight IPS states. Statistical tests show that the relative expression and mean lifetime of certain IPS states have been altered in SZ. Namely, patients spend (45%) less time in a globally coherent state and a subcortical-centered state, and (40%) more time in states reflecting anticoupling within the cognitive control network, compared to the HC. Moreover, the transition profile (between states) reveals a deep temporal structure, shaping two metastates with distinct phase synchrony profiles. A metastate is a collection of states such that within-metastate transitions are more probable than across. Remarkably, metastate occupation balance is altered in SZ, in favor of the less synchronous metastate that promotes disconnection across networks. Furthermore, the trajectory of IPS patterns is less efficient, less smooth, and more restricted in SZ subjects, compared to the HC. Finally, a regression analysis confirms the diagnostic value of the defined IPS measures for SZ identification, highlighting the distinctive role of metastate proportion. Our results suggest that the proposed IPS features may be used for classification studies and for characterizing phase synchrony modes in other (clinical) populations.

Functional outcome is tied to dynamic brain states after mild to moderate traumatic brain injury

The current study set out to investigate the dynamic functional connectome in relation to long-term recovery after mild to moderate traumatic brain injury (TBI). Longitudinal resting-state functional MRI data were collected (at 1 and 3 months postinjury) from a prospectively enrolled cohort consisting of 68 patients with TBI (92% mild TBI) and 20 healthy subjects. Patients underwent a neuropsychological assessment at 3 months postinjury. Outcome was measured using the Glasgow Outcome Scale Extended (GOS-E) at 6 months postinjury. The 57 patients who completed the GOS-E were classified as recovered completely (GOS-E = 8; n = 37) or incompletely (GOS-E < 8; n = 20). Neuropsychological test scores were similar for all groups. Patients with incomplete recovery spent less time in a segregated brain state compared to recovered patients during the second visit. Also, these patients moved less frequently from one meta-state to another as compared to healthy controls and recovered patients. Furthermore, incomplete recovery was associated with disruptions in cyclic state transition patterns, called attractors, during both visits. This study demonstrates that poor long-term functional recovery is associated with alterations in dynamics between brain networks, which becomes more marked as a function of time. These results could be related to psychological processes rather than injury-effects, which is an interesting area for further work. Another natural progression of the current study is to examine whether these dynamic measures can be used to monitor treatment effects.