Multisubject independent component analysis of fMRI: a decade of intrinsic networks, default mode, and neurodiagnostic discovery

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

BACKGROUND

The onset of anorexia nervosa (AN) frequently occurs during adolescence and is associated with preoccupation with body weight and shape and extreme underweight. Altered resting state functional connectivity in the brain has been described in individuals with AN, but only from a static perspective. The current study investigated the temporal dynamics of functional connectivity in adolescents with AN and how it relates to clinical features.

METHOD

99 female patients acutely ill with AN and 99 pairwise age-matched female healthy control (HC) participants were included in the study. Using resting-state functional MRI data and an established sliding-window analytic approach, we identified dynamic resting-state functional connectivity states and extracted dynamic indices such as dwell time (the duration spent in a state), fraction time (the proportion of the total time occupied by a state), and number of transitions (number of switches) from one state to another, to test for group differences.

RESULTS

Individuals with AN had relatively reduced fraction time in a mildly connected state with pronounced connectivity within the default mode network (DMN) and an overall reduced number of transitions between states.

CONCLUSIONS

These findings revealed by a dynamic, but not static analytic approach might hint towards a more "rigid" connectivity, a phenomenon commonly observed in internalizing mental disorders, and in AN possibly related to a reduction in energetic costs as a result of nutritional deprivation.

A Telescopic Independent Component Analysis on Functional Magnetic Resonance Imaging Data Set

Brain function can be modeled as the dynamic interactions between functional sources at different spatial scales, and each spatial scale can contain its functional sources with unique information, thus using a single scale may provide an incomplete view of brain function. This paper introduces a novel approach, termed "telescopic independent component analysis (TICA)," designed to construct spatial functional hierarchies and estimate functional sources across multiple spatial scales using fMRI data. The method employs a recursive ICA strategy, leveraging information from a larger network to guide the extraction of information about smaller networks. We apply our model to the default mode network (DMN), visual network (VN), and right frontoparietal network (RFPN). We investigate further on DMN by evaluating the difference between healthy people and individuals with schizophrenia. We show that the TICA approach can detect the spatial hierarchy of DMN, VS, and RFPN. In addition, TICA revealed DMN-associated group differences between cohorts that may not be captured if we focus on a single-scale ICA. In sum, our proposed approach represents a promising new tool for studying functional sources.

Abnormal large-scale resting-state functional networks in anti-N-methyl-D-aspartate receptor encephalitis

Background

Patients with anti-N-methyl-D-aspartate receptor (anti-NMDAR) encephalitis often experience severe symptoms. Resting-state functional MRI (rs-fMRI) has revealed widespread impairment of functional networks in patients. However, the changes in information flow remain unclear. This study aims to investigate the intrinsic functional connectivity (FC) both within and between resting-state networks (RSNs), as well as the alterations in effective connectivity (EC) between these networks.

Methods

Resting-state functional MRI (rs-fMRI) data were collected from 25 patients with anti-NMDAR encephalitis and 30 healthy controls (HCs) matched for age, sex, and educational level. Changes in the intrinsic functional connectivity (FC) within and between RSNs were analyzed using independent component analysis (ICA). The functional interaction between RSNs was identified by granger causality analysis (GCA).

Results

Compared to HCs, patients with anti-NMDAR encephalitis exhibited lower performance on the Wisconsin Card Sorting Test (WCST), both in terms of correct numbers and correct categories. Additionally, these patients demonstrated decreased scores on the Montreal Cognitive Assessment (MoCA). Neuroimaging studies revealed abnormal intra-FC within the default mode network (DMN), increased intra-FC within the visual network (VN) and dorsal attention network (DAN), as well as increased inter-FC between VN and the frontoparietal network (FPN). Furthermore, aberrant effective connectivity (EC) was observed among the DMN, DAN, FPN, VN, and somatomotor network (SMN).

Conclusion

Patients with anti-NMDAR encephalitis displayed noticeable deficits in both memory and executive function. Notably, these patients exhibited widespread impairments in intra-FC, inter-FC, and EC. These results may help to explain the pathophysiological mechanism of anti-NMDAR encephalitis.

Unraveling the Neural Landscape of Mental Disorders using Double Functional Independent Primitives (dFIPs)

Mental illnesses extract a high personal and societal cost, and thus explorations of the links between mental illness and functional connectivity in the brain are critical. Investigating major mental illnesses, believed to arise from disruptions in sophisticated neural connections, allows us to comprehend how these neural network disruptions may be linked to altered cognition, emotional regulation, and social interactions. Although neuroimaging has opened new avenues to explore neural alterations linked to mental illnesses, the field still requires precise and sensitive methodologies to inspect these neural substrates of various psychological disorders. In this study, we employ a hierarchical methodology to derive double functionally independent primitives (dFIPs) from resting state functional magnetic resonance neuroimaging data (rs-fMRI). These dFIPs encapsulate canonical overlapping patterns of functional network connectivity (FNC) within the brain. Our investigation focuses on the examination of how combinations of these dFIPs relate to different mental disorder diagnoses. The central aim is to unravel the complex patterns of FNC that correspond to the diverse manifestations of mental illnesses. To achieve this objective, we used a large brain imaging dataset from multiple sites, comprising 5805 total individuals diagnosed with schizophrenia (SCZ), autism spectrum disorder (ASD), bipolar disorder (BPD), major depressive disorder (MDD), and controls. The key revelations of our study unveil distinct patterns associated with each mental disorder through the combination of dFIPs. Notably, certain individual dFIPs exhibit disorder-specific characteristics, while others demonstrate commonalities across disorders. This approach offers a novel, data-driven synthesis of intricate neuroimaging data, thereby illuminating the functional changes intertwined with various mental illnesses. Our results show distinct signatures associated with psychiatric disorders, revealing unique connectivity patterns such as heightened cerebellar connectivity in SCZ and sensory domain hyperconnectivity in ASD, both contrasted with reduced cerebellar-subcortical connectivity. Utilizing the dFIP concept, we pinpoint specific functional connections that differentiate healthy controls from individuals with mental illness, underscoring its utility in identifying neurobiological markers. In summary, our findings delineate how dFIPs serve as unique fingerprints for different mental disorders.

Association between the oral microbiome and brain resting state connectivity in schizophrenia

Recent microbiome-brain axis findings have shown evidence of the modulation of microbiome community as an environmental mediator in brain function and psychiatric illness. This work is focused on the role of the microbiome in understanding a rarely investigated environmental involvement in schizophrenia (SZ), especially in relation to brain circuit dysfunction. We leveraged high throughput microbial 16s rRNA sequencing and functional neuroimaging techniques to enable the delineation of microbiome-brain network links in SZ. N = 213 SZ and healthy control subjects were assessed for the oral microbiome. Among them, 139 subjects were scanned by resting-state functional magnetic resonance imaging (rsfMRI) to derive brain functional connectivity. We found a significant microbiome compositional shift in SZ beta diversity (weighted UniFrac distance, p = 6 × 10-3; Bray-Curtis distance p = 0.021). Fourteen microbial species involving pro-inflammatory and neurotransmitter signaling and H2S production, showed significant abundance alterations in SZ. Multivariate analysis revealed one pair of microbial and functional connectivity components showing a significant correlation of 0.46. Thirty five percent of microbial species and 87.8 % of brain functional network connectivity from each component also showed significant differences between SZ and healthy controls with strong performance in classifying SZ from healthy controls, with an area under curve (AUC) = 0.84 and 0.87, respectively. The results suggest a potential link between oral microbiome dysbiosis and brain functional connectivity alteration in relation to SZ, possibly through immunological and neurotransmitter signaling pathways and the hypothalamic-pituitary-adrenal axis, supporting for future work in characterizing the role of oral microbiome in mediating effects on SZ brain functional activity.

ADHD diagnosis guided by functional brain networks combined with domain knowledge

Utilizing functional magnetic resonance imaging (fMRI) to model functional brain networks (FBNs) is increasingly prominent in attention-deficit/hyperactivity disorder (ADHD) research, revealing neural impact and mechanisms through the exploration of activated brain regions. However, current FBNs-based methods face two major challenges. The primary challenge stems from the limitations of existing modeling methods in accurately capturing both regional correlations and long-distance dependencies (LDDs) within the dynamic brain, thereby affecting the diagnostic accuracy of FBNs as biomarkers. Additionally, limited sample size and class imbalance also pose a challenge to the learning performance of the model. To address the issues, we propose an automated diagnostic framework, which integrates modeling, multimodal fusion, and classification into a unified process. It aims to extract representative FBNs and efficiently incorporate domain knowledge to guide ADHD classification. Our work mainly includes three-fold: 1) A multi-head attention-based region-enhancement module (MAREM) is designed to simultaneously capture regional correlations and LDDs across the entire sequence of brain activity, which facilitates the construction of representative FBNs. 2) The multimodal supplementary learning module (MSLM) is proposed to integrate domain knowledge from phenotype data with FBNs from neuroimaging data, achieving information complementarity and alleviating the problems of insufficient medical data and unbalanced sample categories. 3) An ADHD automatic diagnosis framework guided by FBNs and domain knowledge (ADF-FAD) is proposed to help doctors make more accurate decisions, which is applied to the ADHD-200 dataset to confirm its effectiveness. The results indicate that the FBNs extracted by MAREM perform well in modeling and classification. After with MSLM, the model achieves accuracy of 92.4%, 74.4%, and 80% at NYU, PU, and KKI, respectively, demonstrating its ability to effectively capture crucial information related to ADHD diagnosis. Codes are available at https://github.com/zhuimengxuebao/ADF-FAD.

Learning Spatiotemporal Brain Dynamics in Adolescents via Multimodal MEG and fMRI Data Fusion Using Joint Tensor/Matrix Decomposition

OBJECTIVE

Brain function is understood to be regulated by complex spatiotemporal dynamics, and can be characterized by a combination of observed brain response patterns in time and space. Magnetoencephalography (MEG), with its high temporal resolution, and functional magnetic resonance imaging (fMRI), with its high spatial resolution, are complementary imaging techniques with great potential to reveal information about spatiotemporal brain dynamics. Hence, the complementary nature of these imaging techniques holds much promise to study brain function in time and space, especially when the two data types are allowed to fully interact.

METHODS

We employed coupled tensor/matrix factorization (CMTF) to extract joint latent components in the form of unique spatiotemporal brain patterns that can be used to study brain development and function on a millisecond scale.

RESULTS

Using the CMTF model, we extracted distinct brain patterns that revealed fine-grained spatiotemporal brain dynamics and typical sensory processing pathways informative of high-level cognitive functions in healthy adolescents. The components extracted from multimodal tensor fusion possessed better discriminative ability between high- and low-performance subjects than single-modality data-driven models.

CONCLUSION

Multimodal tensor fusion successfully identified spatiotemporal brain dynamics of brain function and produced unique components with high discriminatory power.

SIGNIFICANCE

The CMTF model is a promising tool for high-order, multimodal data fusion that exploits the functional resolution of MEG and fMRI, and provides a comprehensive picture of the developing brain in time and space.

Static and Dynamic Cross-Network Functional Connectivity Shows Elevated Entropy in Schizophrenia Patients

Schizophrenia (SZ) patients exhibit abnormal static and dynamic functional connectivity across various brain domains. We present a novel approach based on static and dynamic inter-network connectivity entropy (ICE), which represents the entropy of a given network's connectivity to all the other brain networks. This novel approach enables the investigation of how connectivity strength is heterogeneously distributed across available targets in both SZ patients and healthy controls. We analyzed fMRI data from 151 schizophrenia patients and demographically matched 160 healthy controls. Our assessment encompassed both static and dynamic ICE, revealing significant differences in the heterogeneity of connectivity levels across available brain networks between SZ patients and healthy controls (HC). These networks are associated with subcortical (SC), auditory (AUD), sensorimotor (SM), visual (VIS), cognitive control (CC), default mode network (DMN) and cerebellar (CB) functional brain domains. Elevated ICE observed in individuals with SZ suggests that patients exhibit significantly higher randomness in the distribution of time-varying connectivity strength across functional regions from each source network, compared to healthy control group. C-means fuzzy clustering analysis of functional ICE correlation matrices revealed that SZ patients exhibit significantly higher occupancy weights in clusters with weak, low-scale functional entropy correlation, while the control group shows greater occupancy weights in clusters with strong, large-scale functional entropy correlation. k-means clustering analysis on time-indexed ICE vectors revealed that cluster with highest ICE have higher occupancy rates in SZ patients whereas clusters characterized by lowest ICE have larger occupancy rates for control group. Furthermore, our dynamic ICE approach revealed that it appears healthy for a brain to primarily circulate through complex, less structured connectivity patterns, with occasional transitions into more focused patterns. However, individuals with SZ seem to struggle with transiently attaining these more focused and structured connectivity patterns. Proposed ICE measure presents a novel framework for gaining deeper insights into understanding mechanisms of healthy and disease brain states and a substantial step forward in the developing advanced methods of diagnostics of mental health conditions.

Multiscale Modes of Functional Brain Connectivity

Information processing in the brain spans from localised sensorimotor processes to higher-level cognition that integrates across multiple regions. Interactions between and within these subsystems enable multiscale information processing. Despite this multiscale characteristic, functional brain connectivity is often either estimated based on 10-30 distributed modes or parcellations with 100-1000 localised parcels, both missing across-scale functional interactions. We present Multiscale Probabilistic Functional Modes (mPFMs), a new mapping which comprises modes over various scales of granularity, thus enabling direct estimation of functional connectivity within- and across-scales. Crucially, mPFMs emerged from data-driven multilevel Bayesian modelling of large functional MRI (fMRI) populations. We demonstrate that mPFMs capture both distributed brain modes and their co-existing subcomponents. In addition to validating mPFMs using simulations and real data, we show that mPFMs can predict ~900 personalised traits from UK Biobank more accurately than current standard techniques. Therefore, mPFMs can offer a paradigm shift in functional connectivity modelling and yield enhanced fMRI biomarkers for traits and diseases.

A spatially constrained independent component analysis jointly informed by structural and functional network connectivity

There are a growing number of neuroimaging studies motivating joint structural and functional brain connectivity. Brain connectivity of different modalities provides insight into brain functional organization by leveraging complementary information, especially for brain disorders such as schizophrenia. In this paper, we propose a multi-modal independent component analysis (ICA) model that utilizes information from both structural and functional brain connectivity guided by spatial maps to estimate intrinsic connectivity networks (ICNs). Structural connectivity is estimated through whole-brain tractography on diffusion-weighted MRI (dMRI), while functional connectivity is derived from resting-state functional MRI (rs-fMRI). The proposed structural-functional connectivity and spatially constrained ICA (sfCICA) model estimates ICNs at the subject level using a multi-objective optimization framework. We evaluated our model using synthetic and real datasets (including dMRI and rs-fMRI from 149 schizophrenia patients and 162 controls). Multi-modal ICNs revealed enhanced functional coupling between ICNs with higher structural connectivity, improved modularity, and network distinction, particularly in schizophrenia. Statistical analysis of group differences showed more significant differences in the proposed model compared to the unimodal model. In summary, the sfCICA model showed benefits from being jointly informed by structural and functional connectivity. These findings suggest advantages in simultaneously learning effectively and enhancing connectivity estimates using structural connectivity.

Coactivation pattern analysis reveals altered whole-brain functional transient dynamics in autism spectrum disorder

Recent studies on autism spectrum disorder (ASD) have identified recurring states dominated by similar coactivation pattern (CAP) and revealed associations between dysfunction in seed-based large-scale brain networks and clinical symptoms. However, the presence of abnormalities in moment-to-moment whole-brain dynamics in ASD remains uncertain. In this study, we employed seed-free CAP analysis to identify transient brain activity configurations and investigate dynamic abnormalities in ASD. We utilized a substantial multisite resting-state fMRI dataset consisting of 354 individuals with ASD and 446 healthy controls (HCs, from HC groups and 2). CAP were generated from a subgroup of all HC subjects (HC group 1) through temporal K-means clustering, identifying four CAPs. These four CAPs exhibited either the activation or inhibition of the default mode network (DMN) and were grouped into two pairs with opposing spatial CAPs. CAPs for HC group 2 and ASD were identified by their spatial similarity to those for HC group 1. Compared with individuals in HC group 2, those with ASD spent more time in CAPs involving the ventral attention network but less time in CAPs related to executive control and the dorsal attention network. Support vector machine analysis demonstrated that the aberrant dynamic characteristics of CAPs achieved an accuracy of 74.87% in multisite classification. In addition, we used whole-brain dynamics to predict symptom severity in ASD. Our findings revealed whole-brain dynamic functional abnormalities in ASD from a single transient perspective, emphasizing the importance of the DMN in abnormal dynamic functional activity in ASD and suggesting that temporally dynamic techniques offer novel insights into time-varying neural processes.

Orthogonal extended infomax algorithm

Objective.The extended infomax algorithm for independent component analysis (ICA) can separate sub- and super-Gaussian signals but converges slowly as it uses stochastic gradient optimization. In this paper, an improved extended infomax algorithm is presented that converges much faster.Approach.Accelerated convergence is achieved by replacing the natural gradient learning rule of extended infomax by a fully-multiplicative orthogonal-group based update scheme of the ICA unmixing matrix, leading to an orthogonal extended infomax algorithm (OgExtInf). The computational performance of OgExtInf was compared with original extended infomax and with two fast ICA algorithms: the popular FastICA and Picard, a preconditioned limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) algorithm belonging to the family of quasi-Newton methods.Main results.OgExtInf converges much faster than original extended infomax. For small-size electroencephalogram (EEG) data segments, as used for example in online EEG processing, OgExtInf is also faster than FastICA and Picard.Significance.OgExtInf may be useful for fast and reliable ICA, e.g. in online systems for epileptic spike and seizure detection or brain-computer interfaces.

Mapping dynamic spatial patterns of brain function with spatial-wise attention

Objective: Using functional magnetic resonance imaging (fMRI) and deep learning to discover the spatial pattern of brain function, or functional brain networks (FBNs) has been attracted many reseachers. Most existing works focus on static FBNs or dynamic functional connectivity among fixed spatial network nodes, but ignore the potential dynamic/time-varying characteristics of the spatial networks themselves. And most of works based on the assumption of linearity and independence, that oversimplify the relationship between blood-oxygen level dependence signal changes and the heterogeneity of neuronal activity within voxels.Approach: To overcome these problems, we proposed a novel spatial-wise attention (SA) based method called Spatial and Channel-wise Attention Autoencoder (SCAAE) to discover the dynamic FBNs without the assumptions of linearity or independence. The core idea of SCAAE is to apply the SA to generate FBNs directly, relying solely on the spatial information present in fMRI volumes. Specifically, we trained the SCAAE in a self-supervised manner, using the autoencoder to guide the SA to focus on the activation regions. Experimental results show that the SA can generate multiple meaningful FBNs at each fMRI time point, which spatial similarity are close to the FBNs derived by known classical methods, such as independent component analysis.Main results: To validate the generalization of the method, we evaluate the approach on HCP-rest, HCP-task and ADHD-200 dataset. The results demonstrate that SA mechanism can be used to discover time-varying FBNs, and the identified dynamic FBNs over time clearly show the process of time-varying spatial patterns fading in and out.Significance: Thus we provide a novel method to understand human brain better. Code is available athttps://github.com/WhatAboutMyStar/SCAAE.

Spatial-temporal convolutional attention for discovering and characterizing functional brain networks in task fMRI

Functional brain networks (FBNs) are spatial patterns of brain function that play a critical role in understanding human brain function. There are many proposed methods for mapping the spatial patterns of brain function, however they oversimplify the underlying assumptions of brain function and have various limitations such as linearity and independence. Additionally, current methods fail to account for the dynamic nature of FBNs, which limits their effectiveness in accurately characterizing these networks. To address these limitations, we present a novel deep learning and spatial-wise attention based model called Spatial-Temporal Convolutional Attention (STCA) to accurately model dynamic FBNs. Specifically, we train STCA in a self-supervised manner by utilizing a Convolutional Autoencoder to guide the STCA module in assigning higher attention weights to regions of functional activity. To validate the reliability of the results, we evaluate our approach on the HCP-task motor behavior dataset, the experimental results demonstrate that the STCA derived FBNs have higher spatial similarity with the templates and that the spatial similarity between the templates and the FBNs derived by STCA fluctuates with the task design over time, suggesting that STCA can reflect the dynamic changes of brain function, providing a powerful tool to better understand human brain function. Code is available at https://github.com/SNNUBIAI/STCAE.

Dysfunction of the triple-network model is associated with cognitive impairment in patients with cerebral small vessel disease

Purpose

This study aimed to demonstrate the correlations between the altered functional connectivity patterns in the triple-network model and cognitive impairment in patients with cerebral small vascular disease (CSVD).

Methods

Resting-state functional magnetic resonance imaging data were obtained from 22 patients with CSVD and 20 healthy controls. The resting-state data were analyzed using independent component analysis and functional network connectivity (FNC) analysis to explore the functional alterations in the intrinsic triple-network model including the salience network (SN), default mode network (DMN), and central executive network (CEN), and their correlations with the cognitive deficits and clinical observations in the patients with CSVD.

Results

Compared to the healthy controls, the patients with CSVD exhibited increased connectivity patterns in the CEN-DMN and decreased connectivity patterns in the DMN-SN, CEN-SN, intra-SN, and intra-DMN. Significant negative correlations were detected between the intra-DMN connectivity pattern and the Montreal Cognitive Assessment (MoCA) total scores (r = -0.460, p = 0.048) and MoCA abstraction scores (r = -0.565, p = 0.012), and a positive correlation was determined between the intra-SN connectivity pattern and the MoCA abstraction scores (r = 0.491, p = 0.033).

Conclusions

Our study findings suggest that the functional alterations in the triple-network model are associated with the cognitive deficits in patients with CSVD and shed light on the importance of the triple-network model in the pathogenesis of CSVD.

Association between the oral microbiome and brain resting state connectivity in schizophrenia

Recent microbiome-brain axis findings have shown evidence of the modulation of microbiome community as an environmental mediator in brain function and psychiatric illness. This work is focused on the role of the microbiome in understanding a rarely investigated environmental involvement in schizophrenia (SZ), especially in relation to brain circuit dysfunction. We leveraged high throughput microbial 16s rRNA sequencing and functional neuroimaging techniques to enable the delineation of microbiome-brain network links in SZ. N=213 SZ and healthy control (HC) subjects were assessed for the oral microbiome. Among them, 139 subjects were scanned by resting-state functional magnetic resonance imaging (rsfMRI) to derive brain functional connectivity. We found a significant microbiome compositional shift in SZ beta diversity (weighted UniFrac distance, p= 6×10 -3 ; Bray-Curtis distance p = 0.021). Fourteen microbial species involving pro-inflammatory and neurotransmitter signaling and H 2 S production, showed significant abundance alterations in SZ. Multivariate analysis revealed one pair of microbial and functional connectivity components showing a significant correlation of 0.46. Thirty five percent of microbial species and 87.8% of brain functional network connectivity from each component also showed significant differences between SZ and HC with strong performance in classifying SZ from HC, with an area under curve (AUC) = 0.84 and 0.87, respectively. The results suggest a potential link between oral microbiome dysbiosis and brain functional connectivity alteration in relation to SZ, possibly through immunological and neurotransmitter signaling pathways and the hypothalamic-pituitary-adrenal axis, supporting for future work in characterizing the role of oral microbiome in mediating effects on SZ brain functional activity.

Functional magnetic resonance imaging study of group independent components underpinning item responses to paranoid-depressive scale

BACKGROUND

Our study expand upon a large body of evidence in the field of neuropsychiatric imaging with cognitive, affective and behavioral tasks, adapted for the functional magnetic resonance imaging (MRI) (fMRI) experimental environment. There is sufficient evidence that common networks underpin activations in task-based fMRI across different mental disorders.

AIM

To investigate whether there exist specific neural circuits which underpin differential item responses to depressive, paranoid and neutral items (DN) in patients respectively with schizophrenia (SCZ) and major depressive disorder (MDD).

METHODS

60 patients were recruited with SCZ and MDD. All patients have been scanned on 3T magnetic resonance tomography platform with functional MRI paradigm, comprised of block design, including blocks with items from diagnostic paranoid (DP), depression specific (DS) and DN from general interest scale. We performed a two-sample t-test between the two groups-SCZ patients and depressive patients. Our purpose was to observe different brain networks which were activated during a specific condition of the task, respectively DS, DP, DN.

RESULTS

Several significant results are demonstrated in the comparison between SCZ and depressive groups while performing this task. We identified one component that is task-related and independent of condition (shared between all three conditions), composed by regions within the temporal (right superior and middle temporal gyri), frontal (left middle and inferior frontal gyri) and limbic/salience system (right anterior insula). Another component is related to both diagnostic specific conditions (DS and DP) e.g. It is shared between DEP and SCZ, and includes frontal motor/language and parietal areas. One specific component is modulated preferentially by to the DP condition, and is related mainly to prefrontal regions, whereas other two components are significantly modulated with the DS condition and include clusters within the default mode network such as posterior cingulate and precuneus, several occipital areas, including lingual and fusiform gyrus, as well as parahippocampal gyrus. Finally, component 12 appeared to be unique for the neutral condition. In addition, there have been determined circuits across components, which are either common, or distinct in the preferential processing of the sub-scales of the task.

CONCLUSION

This study has delivers further evidence in support of the model of trans-disciplinary cross-validation in psychiatry.

Functional Connectivity Between Default Mode and Ventral Attention Networks Mediates the Effects of Chronotype on Daily Physical Activity

Daily physical activity (dPA) is closely related to circadian rhythm and chronotype. The functional connectivity (FC) within or between the default mode (DMN) and ventral attention network (vAN) were associated with dPA and chronotype. DMN-vAN FC was investigated for its role in chronotype and dPA. 153 participants completed the reduced version of the Morningness-Eveningness Questionnaire (rMEQ), dPA was measured via actigraphy (5-day), and then resting-state fMRI scans were performed. rMEQ scores and steps recorded by the actigraphic devices (with each hour as the time window to calculate steps for five consecutive days per hour, subsequently yielding the maximum number of steps and its corresponding time, ie, SM and SMT) represent chronotype and dPA respectively. The results found that the rMEQ scores were significantly negatively correlated with SMT. The positive correlation between the rMEQ scores and the DMN-vAN FC was significant. There were also significant positive correlations between SMT and DMN-vAN FC. Further analysis revealed that DMN-vAN mediates the relationship between chronotype and SMT. The FC of DMN-vAN may be the underlying neural mechanism through which chronotype influences dPA. These findings could support the development of reasonable activity schedules or specific intervention programs to improve physical health.

Modulatory effects of fMRI acquisition time of day, week and year on adolescent functional connectomes across spatial scales: Implications for inference

Metabolic, hormonal, autonomic and physiological rhythms may have a significant impact on cerebral hemodynamics and intrinsic brain synchronization measured with fMRI (the resting-state connectome). The impact of their characteristic time scales (hourly, circadian, seasonal), and consequently scan timing effects, on brain topology in inherently heterogeneous developing connectomes remains elusive. In a cohort of 4102 early adolescents with resting-state fMRI (median age = 120.0 months; 53.1 % females) from the Adolescent Brain Cognitive Development Study, this study investigated associations between scan time-of-day, time-of-week (school day vs weekend) and time-of-year (school year vs summer vacation) and topological properties of resting-state connectomes at multiple spatial scales. On average, participants were scanned around 2 pm, primarily during school days (60.9 %), and during the school year (74.6 %). Scan time-of-day was negatively correlated with multiple whole-brain, network-specific and regional topological properties (with the exception of a positive correlation with modularity), primarily of visual, dorsal attention, salience, frontoparietal control networks, and the basal ganglia. Being scanned during the weekend (vs a school day) was correlated with topological differences in the hippocampus and temporoparietal networks. Being scanned during the summer vacation (vs the school year) was consistently positively associated with multiple topological properties of bilateral visual, and to a lesser extent somatomotor, dorsal attention and temporoparietal networks. Time parameter interactions suggested that being scanned during the weekend and summer vacation enhanced the positive effects of being scanned in the morning. Time-of-day effects were overall small but spatially extensive, and time-of-week and time-of-year effects varied from small to large (Cohen's f ≤ 0.1, Cohen's d<0.82, p < 0.05). Together, these parameters were also positively correlated with temporal fMRI signal variability but only in the left hemisphere. Finally, confounding effects of scan time parameters on relationships between connectome properties and cognitive task performance were assessed using the ABCD neurocognitive battery. Although most relationships were unaffected by scan time parameters, their combined inclusion eliminated associations between properties of visual and somatomotor networks and performance in the Matrix Reasoning and Pattern Comparison Processing Speed tasks. Thus, scan time of day, week and year may impact measurements of adolescent brain's functional circuits, and should be accounted for in studies on their associations with cognitive performance, in order to reduce the probability of incorrect inference.

Resting state network connectivity is associated with cognitive flexibility performance in youth in the Adolescent Brain Cognitive Development Study

Cognitive flexibility is an executive functioning skill that develops in childhood, and when impaired, has transdiagnostic implications for psychiatric disorders. To identify how intrinsic neural architecture at rest is linked to cognitive flexibility performance, we used the data-driven method of independent component analysis (ICA) to investigate resting state networks (RSNs) and their whole-brain connectivity associated with levels of cognitive flexibility performance in children. We hypothesized differences by cognitive flexibility performance in RSN connectivity strength in cortico-striatal circuitry, which would manifest via the executive control network, right and left frontoparietal networks (FPN), salience network, default mode network (DMN), and basal ganglia network. We selected participants from the Adolescent Brain Cognitive Development (ABCD) Study who scored at the 25th, ("CF-Low"), 50th ("CF-Average"), or 75th percentiles ("CF-High") on a cognitive flexibility task, were early to middle puberty, and did not exhibit significant psychopathology (n = 967, 47.9% female; ages 9-10). We conducted whole-brain ICA, identifying 14 well-characterized RSNs. Groups differed in connectivity strength in the right FPN, anterior DMN, and posterior DMN. Planned comparisons indicated CF-High had stronger connectivity between right FPN and supplementary motor/anterior cingulate than CF-Low. CF-High had more anti-correlated connectivity between anterior DMN and precuneus than CF-Average. CF-Low had stronger connectivity between posterior DMN and supplementary motor/anterior cingulate than CF-Average. Post-hoc correlations with reaction time by trial type demonstrated significant associations with connectivity. In sum, our results suggest childhood cognitive flexibility performance is associated with DMN and FPN connectivity strength at rest, and that there may be optimal levels of connectivity associated with task performance that vary by network.

Identifying canonical and replicable multi-scale intrinsic connectivity networks in 100k+ resting-state fMRI datasets

Despite the known benefits of data-driven approaches, the lack of approaches for identifying functional neuroimaging patterns that capture both individual variations and inter-subject correspondence limits the clinical utility of rsfMRI and its application to single-subject analyses. Here, using rsfMRI data from over 100k individuals across private and public datasets, we identify replicable multi-spatial-scale canonical intrinsic connectivity network (ICN) templates via the use of multi-model-order independent component analysis (ICA). We also study the feasibility of estimating subject-specific ICNs via spatially constrained ICA. The results show that the subject-level ICN estimations vary as a function of the ICN itself, the data length, and the spatial resolution. In general, large-scale ICNs require less data to achieve specific levels of (within- and between-subject) spatial similarity with their templates. Importantly, increasing data length can reduce an ICN's subject-level specificity, suggesting longer scans may not always be desirable. We also find a positive linear relationship between data length and spatial smoothness (possibly due to averaging over intrinsic dynamics), suggesting studies examining optimized data length should consider spatial smoothness. Finally, consistency in spatial similarity between ICNs estimated using the full data and subsets across different data lengths suggests lower within-subject spatial similarity in shorter data is not wholly defined by lower reliability in ICN estimates, but may be an indication of meaningful brain dynamics which average out as data length increases.

Aberrant dynamic functional network connectivity in type 2 diabetes mellitus individuals

An increasing number of recent brain imaging studies are dedicated to understanding the neuro mechanism of cognitive impairment in type 2 diabetes mellitus (T2DM) individuals. In contrast to efforts to date that are limited to static functional connectivity, here we investigate abnormal connectivity in T2DM individuals by characterizing the time-varying properties of brain functional networks. Using group independent component analysis (GICA), sliding-window analysis, and k-means clustering, we extracted thirty-one intrinsic connectivity networks (ICNs) and estimated four recurring brain states. We observed significant group differences in fraction time (FT) and mean dwell time (MDT), and significant negative correlation between the Montreal Cognitive Assessment (MoCA) scores and FT/MDT. We found that in the T2DM group the inter- and intra-network connectivity decreases and increases respectively for the default mode network (DMN) and task-positive network (TPN). We also found alteration in the precuneus network (PCUN) and enhanced connectivity between the salience network (SN) and the TPN. Our study provides evidence of alterations of large-scale resting networks in T2DM individuals and shed light on the fundamental mechanisms of neurocognitive deficits in T2DM.

Denoising brain networks using a fixed mathematical phase change in independent component analysis of magnitude-only fMRI data

Brain networks extracted by independent component analysis (ICA) from magnitude-only fMRI data are usually denoised using various amplitude-based thresholds. By contrast, spatial source phase (SSP) or the phase information of ICA brain networks extracted from complex-valued fMRI data, has provided a simple yet effective way to perform the denoising using a fixed phase change. In this work, we extend the approach to magnitude-only fMRI data to avoid testing various amplitude thresholds for denoising magnitude maps extracted by ICA, as most studies do not save the complex-valued data. The main idea is to generate a mathematical SSP map for a magnitude map using a mapping framework, and the mapping framework is built using complex-valued fMRI data with a known SSP map. Here we leverage the fact that the phase map derived from phase fMRI data has similar phase information to the SSP map. After verifying the use of the magnitude data of complex-valued fMRI, this framework is generalized to work with magnitude-only data, allowing use of our approach even without the availability of the corresponding phase fMRI datasets. We test the proposed method using both simulated and experimental fMRI data including complex-valued data from University of New Mexico and magnitude-only data from Human Connectome Project. The results provide evidence that the mathematical SSP denoising with a fixed phase change is effective for denoising spatial maps from magnitude-only fMRI data in terms of retaining more BOLD-related activity and fewer unwanted voxels, compared with amplitude-based thresholding. The proposed method provides a unified and efficient SSP approach to denoise ICA brain networks in fMRI data.

A novel subject-wise dictionary learning approach using multi-subject fMRI spatial and temporal components

The conventional dictionary learning (DL) algorithms aim to adapt the dictionary/sparse code to individual functional magnetic resonance imaging (fMRI) data. Thus, lacking the capability to consolidate the spatiotemporal diversities offered by other subjects. Considering that subject-wise (sw) data matrix can be decomposed into the sparse linear combination of multi-subject (MS) time courses and MS spatial maps, two new algorithms, sw sequential DL (swsDL) and sw block DL (swbDL), have been proposed. They are based on the novel framework, defined by the mixing model, where base matrices prepared by operating a computationally fast sparse spatiotemporal blind source separation method over multiple subjects are employed to adapt the mixing matrices to sw training data. They solve the optimization models formulated using [Formula: see text]/[Formula: see text]-norm penalization/constraints through dictionary/sparse code pair update and alternating minimization approach. They are unique because no existing sparse DL method can incorporate MS spatiotemporal components while updating sw atoms/sparse codes, which can eventually be assembled using neuroscience knowledge to extract group-level dynamics. Various fMRI datasets are used to evaluate and compare the performance of the proposed algorithms with existing state-of-the-art algorithms. Specifically, overall, a [Formula: see text] increase in the mean correlation value and [Formula: see text] reduction in the mean computation time exhibited by swsDL and swbDL, respectively, over the adaptive consistent sequential dictionary algorithm.

Regression models for partially localized fMRI connectivity analyses

Background

Brain functional connectivity analysis of resting-state functional magnetic resonance imaging (fMRI) data is typically performed in a standardized template space assuming consistency of connections across subjects. Analysis methods can come in the form of one-edge-at-a-time analyses or dimension reduction/decomposition methods. Common to these approaches is an assumption that brain regions are functionally aligned across subjects; however, it is known that this functional alignment assumption is often violated.

Methods

In this paper, we use subject-level regression models to explain intra-subject variability in connectivity. Covariates can include factors such as geographic distance between two pairs of brain regions, whether the two regions are symmetrically opposite (homotopic), and whether the two regions are members of the same functional network. Additionally, a covariate for each brain region can be included, to account for the possibility that some regions have consistently higher or lower connectivity. This style of analysis allows us to characterize the fraction of variation explained by each type of covariate. Additionally, comparisons across subjects can then be made using the fitted connectivity regression models, offering a more parsimonious alternative to edge-at-a-time approaches.

Results

We apply our approach to Human Connectome Project data on 268 regions of interest (ROIs), grouped into eight functional networks. We find that a high proportion of variation is explained by region covariates and network membership covariates, while geographic distance and homotopy have high relative importance after adjusting for the number of predictors. We also find that the degree of data repeatability using our connectivity regression model-which uses only partial location information about pairs of ROI's-is comparably as high as the repeatability obtained using full location information.

Discussion

While our analysis uses data that have been transformed into a common template-space, we also envision the method being useful in multi-atlas registration settings, where subject data remains in its own geometry and templates are warped instead. These results suggest the tantalizing possibility that fMRI connectivity analysis can be performed in subject-space, using less aggressive registration, such as simple affine transformations, multi-atlas subject-space registration, or perhaps even no registration whatsoever.

Dynamic functional connectivity changes associated with decreased memory performance in betel quid dependence

The temporal variability of the dynamic functional connectivity (dFC) has been suggested as a useful metric for studying abnormal cognitive function. This study aimed to explore the associations between the temporal properties of dFC and memory performance in betel quid dependence (BQD). Sixty-four BQD individuals and 47 gender- and age-matched healthy controls (HCs) underwent functional magnetic resonance imaging and a series of neuropsychological assessments. The dFC was constructed by calculating the Pearson correlation coefficients within a sliding window and was clustered into three functional connectivity states using k-means clustering. The dFC temporal properties derived from the cluster results were compared between the BQD and HC groups. The results showed that States 1 and 3 featured more frequent and weak connectivity, and State 2 featured less frequent and strong connectivity. There were significant differences for mean dwell time (MDT) in State 3 (p = 0.022) and fraction of time in State 2 (p = 0.018) between the BQD and HC groups. Pearson correlation analyses showed that the MDT in State 1 was negatively correlated with long delay free recall and short delay free recall, and the MDT in State 3 was positively correlated with false positive of long delay recall. Our findings provide strong evidence that MDT match the memory performance and suggest new insights into the pathophysiological mechanism of memory disorders in BQD individuals.

Functional brain network identification and fMRI augmentation using a VAE-GAN framework

Recently, deep learning models have achieved superior performance for mapping functional brain networks from functional magnetic resonance imaging (fMRI) data compared with traditional methods. However, due to the lack of sufficient data and the high dimensionality of brain volume, deep learning models of fMRI tend to suffer from overfitting. In addition, existing methods rarely studied fMRI data augmentation and its application. To address these issues, we developed a VAE-GAN framework that combined a VAE (variational auto-encoder) with a GAN (generative adversarial net) for functional brain network identification and fMRI augmentation. As a generative model, the VAE-GAN models the distribution of fMRI so that it enables the extraction of more generalized features, and thus relieve the overfitting issue. The VAE-GAN is easier to train on fMRI than a standard GAN since it uses latent variables from VAE to generate fake data rather than relying on random noise that is used in a GAN, and it can generate higher quality of fake data than VAE since the discriminator can promote the training of the generator. In other words, the VAE-GAN inherits the advantages of VAE and GAN and avoids their limitations in modeling of fMRI data. Extensive experiments on task fMRI datasets from HCP have proved the effectiveness and superiority of the proposed VAE-GAN framework for identifying both temporal features and functional brain networks compared with existing models, and the quality of fake data is higher than those from VAE and GAN. The results on resting state fMRI of Attention Deficit Hyperactivity Disorder (ADHD)-200 dataset further demonstrated that the fake data generated by the VAE-GAN can help improve the performance of brain network modeling and ADHD classification.

Extracting functional connectivity brain networks at the resting state from pulsed arterial spin labeling data

Introduction

Functional connectivity in the brain is often studied with blood oxygenation level dependent (BOLD) resting state functional magnetic resonance imaging (rsfMRI), but the BOLD signal is several steps removed from neuronal activity. Arterial spin labeling (ASL), particularly pulsed ASL (PASL), has also the capacity to measure the blood-flow changes in response to activity. In this paper, we investigated the feasibility of extracting major brain networks from PASL data, in contrast with rsfMRI analsyis.

Materials and methods

In this retrospective study, we analyzed a cohort dataset that consists of 21 mild traumatic brain injury (mTBI) patients and 29 healthy controls, which was collected in a previous study. By extracting 10 major brain networks from the data of both PASL and rsfMRI, we contrasted their similarities and differences in the 10 networks extracted from both modalities.

Results

Our data demonstrated that PASL could be used to extract all 10 major brain networks. Eight out of 10 networks demonstrated over 60 % similarity to rsfMRI data. Meanwhile, there are similar but not identical changes in networks detected between mTBI patients and healthy controls with both modalities. Notably, the PASL-extracted default mode network (DMN), other than the rsfMRI-extracted DMN, includes some regions known to be associated with the DMN in other studies. It demonstrated that PASL data can be analyzed to identify resting state networks with reasonable reliability, even without rsfMRI data.

Conclusion

Our analysis provides an opportunity to extract functional connectivity information in heritage datasets in which ASL but not BOLD was collected.

Leading basic modes of spontaneous activity drive individual functional connectivity organization in the resting human brain

Spontaneous activity of the human brain provides a window to explore intrinsic principles of functional organization. However, most studies have focused on interregional functional connectivity. The principles underlying rich repertoires of instantaneous activity remain largely unknown. We apply a recently proposed eigen-microstate analysis to three resting-state functional MRI datasets to identify basic modes that represent fundamental activity patterns that coexist over time. We identify five leading basic modes that dominate activity fluctuations. Each mode exhibits a distinct functional system-dependent coactivation pattern and corresponds to specific cognitive profiles. In particular, the spatial pattern of the first leading basis mode shows the separation of activity between the default-mode and primary and attention regions. Based on theoretical modelling, we further reconstruct individual functional connectivity as the weighted superposition of coactivation patterns corresponding to these leading basic modes. Moreover, these leading basic modes capture sleep deprivation-induced changes in brain activity and interregional connectivity, primarily involving the default-mode and task-positive regions. Our findings reveal a dominant set of basic modes of spontaneous activity that reflect multiplexed interregional coordination and drive conventional functional connectivity, furthering the understanding of the functional significance of spontaneous brain activity.

A deep learning method for autism spectrum disorder identification based on interactions of hierarchical brain networks

BACKGROUND

It has been recently shown that deep learning models exhibited remarkable performance of representing functional Magnetic Resonance Imaging (fMRI) data for the understanding of brain functional activities. With hierarchical structure, deep learning models can infer hierarchical functional brain networks (FBN) from fMRI. However, the applications of the hierarchical FBNs have been rarely studied.

METHODS

In this work, we proposed a hierarchical recurrent variational auto-encoder (HRVAE) to unsupervisedly model the fMRI data. The trained HRVAE encoder can predict hierarchical temporal features from its three hidden layers, and thus can be regarded as a hierarchical feature extractor. Then LASSO (least absolute shrinkage and selection operator) regression was applied to estimate the corresponding hierarchical FBNs. Based on the hierarchical FBNs from each subject, we constructed a novel classification framework for brain disorder identification and test it on the Autism Brain Imaging Data Exchange (ABIDE) dataset, a world-wide multi-site database of autism spectrum disorder (ASD). We analyzed the hierarchy organization of FBNs, and finally used the overlaps of hierarchical FBNs as features to differentiate ASD from typically developing controls (TDC).

RESULTS

The experimental results on 871 subjects from ABIDE dataset showed that the HRVAE model can effectively derive hierarchical FBNs including many well-known resting state networks (RSN). Moreover, the classification result improved the state-of-the-art by achieving a very high accuracy of 82.1 %.

CONCLUSIONS

This work presents a novel data-driven deep learning method using fMRI data for ASD identification, which could provide valuable reference for clinical diagnosis. The classification results suggest that the interactions of hierarchical FBNs have association with brain disorder, which promotes the understanding of FBN hierarchy and could be applied to other brain disorder analysis.

Transdiagnostic indicators predict developmental changes in cognitive control resting-state networks

Over the past decade, transdiagnostic indicators in relation to neurobiological processes have provided extensive insight into youth's risk for psychopathology. During development, exposure to childhood trauma and dysregulation (i.e., so-called AAA symptomology: anxiety, aggression, and attention problems) puts individuals at a disproportionate risk for developing psychopathology and altered network-level neural functioning. Evidence for the latter has emerged from resting-state fMRI studies linking mental health symptoms and aberrations in functional networks (e.g., cognitive control (CCN), default mode networks (DMN)) in youth, although few of these investigations have used longitudinal designs. Herein, we leveraged a three-year longitudinal study to identify whether traumatic exposures and concomitant dysregulation trigger changes in the developmental trajectories of resting-state functional networks involved in cognitive control (N = 190; 91 females; time 1 Mage = 11.81). Findings from latent growth curve analyses revealed that greater trauma exposure predicted increasing connectivity between the CCN and DMN across time. Greater levels of dysregulation predicted reductions in within-network connectivity in the CCN. These findings presented in typically developing youth corroborate connectivity patterns reported in clinical populations, suggesting there is predictive utility in using transdiagnostic indicators to forecast alterations in resting-state networks implicated in psychopathology.

Evaluation of boundaries between mood and psychosis disorder using dynamic functional network connectivity (dFNC) via deep learning classification

The validity and reliability of diagnoses in psychiatry is a challenging topic in mental health. The current mental health categorization is based primarily on symptoms and clinical course and is not biologically validated. Among multiple ongoing efforts, neurological observations alongside clinical evaluations are considered to be potential solutions to address diagnostic problems. The Bipolar-Schizophrenia Network on Intermediate Phenotypes (B-SNIP) has published multiple papers attempting to reclassify psychotic illnesses based on biological rather than symptomatic measures. However, the effort to investigate the relationship between this new categorization approach and other neuroimaging techniques, including resting-state fMRI data, is still limited. This study focused on investigating the relationship between different psychotic disorders categorization methods and resting-state fMRI-based measures called dynamic functional network connectivity (dFNC) using state-of-the-art artificial intelligence (AI) approaches. We applied our method to 613 subjects, including individuals with psychosis and healthy controls, which were classified using both the Diagnostic and Statistical Manual of Mental Disorders (DSM-IV) and the B-SNIP biomarker-based (Biotype) approach. Statistical group differences and cross-validated classifiers were performed within each framework to assess how different categories. Results highlight interesting differences in occupancy in both DSM-IV and Biotype categorizations compared to healthy individuals, which are distributed across specific transient connectivity states. Biotypes tended to show less distinctiveness in occupancy level and included fewer cellwise differences. Classification accuracy obtained by DSM-IV and Biotype categories were both well above chance. Results provided new insights and highlighted the benefits of both DSM-IV and biology-based categories while also emphasizing the importance of future work in this direction, including employing further data types.

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.

Fine temporal brain network structure modularizes and localizes differently in men and women: insights from a novel explainability framework

Deep learning has become an effective tool for classifying biological sex based on functional magnetic resonance imaging (fMRI). However, research on what features within the brain are most relevant to this classification is still lacking. Model interpretability has become a powerful way to understand "black box" deep-learning models, and select features within the input data that are most relevant to the correct classification. However, very little work has been done employing these methods to understand the relationship between the temporal dimension of functional imaging signals and the classification of biological sex. Consequently, less attention has been paid to rectifying problems and limitations associated with feature explanation models, e.g. underspecification and instability. In this work, we first provide a methodology to limit the impact of underspecification on the stability of the measured feature importance. Then, using intrinsic connectivity networks from fMRI data, we provide a deep exploration of sex differences among functional brain networks. We report numerous conclusions, including activity differences in the visual and cognitive domains and major connectivity differences.

Multi-head attention-based masked sequence model for mapping functional brain networks

The investigation of functional brain networks (FBNs) using task-based functional magnetic resonance imaging (tfMRI) has gained significant attention in the field of neuroimaging. Despite the availability of several methods for constructing FBNs, including traditional methods like GLM and deep learning methods such as spatiotemporal self-attention mechanism (STAAE), these methods have design and training limitations. Specifically, they do not consider the intrinsic characteristics of fMRI data, such as the possibility that the same signal value at different time points could represent different brain states and meanings. Furthermore, they overlook prior knowledge, such as task designs, during training. This study aims to overcome these limitations and develop a more efficient model by drawing inspiration from techniques in the field of natural language processing (NLP). The proposed model, called the Multi-head Attention-based Masked Sequence Model (MAMSM), uses a multi-headed attention mechanism and mask training approach to learn different states corresponding to the same voxel values. Additionally, it combines cosine similarity and task design curves to construct a novel loss function. The MAMSM was applied to seven task state datasets from the Human Connectome Project (HCP) tfMRI dataset. Experimental results showed that the features acquired by the MAMSM model exhibit a Pearson correlation coefficient with the task design curves above 0.95 on average. Moreover, the model can extract more meaningful networks beyond the known task-related brain networks. The experimental results demonstrated that MAMSM has great potential in advancing the understanding of functional brain networks.

Developmental and aging resting functional magnetic resonance imaging brain state adaptations in adolescents and adults: A large N (>47K) study

The brain's functional architecture and organization undergo continual development and modification throughout adolescence. While it is well known that multiple factors govern brain maturation, the constantly evolving patterns of time-resolved functional connectivity are still unclear and understudied. We systematically evaluated over 47,000 youth and adult brains to bridge this gap, highlighting replicable time-resolved developmental and aging functional brain patterns. The largest difference between the two life stages was captured in a brain state that indicated coherent strengthening and modularization of functional coupling within the auditory, visual, and motor subdomains, supplemented by anticorrelation with other subdomains in adults. This distinctive pattern, which we replicated in independent data, was consistently less modular or absent in children and presented a negative association with age in adults, thus indicating an overall inverted U-shaped trajectory. This indicates greater synchrony, strengthening, modularization, and integration of the brain's functional connections beyond adolescence, and gradual decline of this pattern during the healthy aging process. We also found evidence that the developmental changes may also bring along a departure from the canonical static functional connectivity pattern in favor of more efficient and modularized utilization of the vast brain interconnections. State-based statistical summary measures presented robust and significant group differences that also showed significant age-related associations. The findings reported in this article support the idea of gradual developmental and aging brain state adaptation processes in different phases of life and warrant future research via lifespan studies to further authenticate the projected time-resolved brain state trajectories.

Identification of Homogeneous Subgroups from Resting-State fMRI Data

The identification of homogeneous subgroups of patients with psychiatric disorders can play an important role in achieving personalized medicine and is essential to provide insights for understanding neuropsychological mechanisms of various mental disorders. The functional connectivity profiles obtained from functional magnetic resonance imaging (fMRI) data have been shown to be unique to each individual, similar to fingerprints; however, their use in characterizing psychiatric disorders in a clinically useful way is still being studied. In this work, we propose a framework that makes use of functional activity maps for subgroup identification using the Gershgorin disc theorem. The proposed pipeline is designed to analyze a large-scale multi-subject fMRI dataset with a fully data-driven method, a new constrained independent component analysis algorithm based on entropy bound minimization (c-EBM), followed by an eigenspectrum analysis approach. A set of resting-state network (RSN) templates is generated from an independent dataset and used as constraints for c-EBM. The constraints present a foundation for subgroup identification by establishing a connection across the subjects and aligning subject-wise separate ICA analyses. The proposed pipeline was applied to a dataset comprising 464 psychiatric patients and discovered meaningful subgroups. Subjects within the identified subgroups share similar activation patterns in certain brain areas. The identified subgroups show significant group differences in multiple meaningful brain areas including dorsolateral prefrontal cortex and anterior cingulate cortex. Three sets of cognitive test scores were used to verify the identified subgroups, and most of them showed significant differences across subgroups, which provides further confirmation of the identified subgroups. In summary, this work represents an important step forward in using neuroimaging data to characterize mental disorders.

Exploring the longitudinal associations of functional network connectivity and psychiatric symptom changes in youth

BACKGROUND

Functional connectivity has been associated with psychiatric problems, both in children and adults, but inconsistencies are present across studies. Prior research has mostly focused on small clinical samples with cross-sectional designs.

METHODS

We adopted a longitudinal design with repeated assessments to investigate associations between functional network connectivity (FNC) and psychiatric problems in youth (9- to 17-year-olds, two time points) from the general population. The largest single-site study of pediatric neurodevelopment was used: Generation R (N = 3,131 with data at either time point). Psychiatric symptoms were measured with the Child Behavioral Checklist as broadband internalizing and externalizing problems, and its eight specific syndrome scales (e.g., anxious-depressed). FNC was assessed with two complementary approaches. First, static FNC (sFNC) was measured with graph theory-based metrics. Second, dynamic FNC (dFNC), where connectivity is allowed to vary over time, was summarized into 5 states that participants spent time in. Cross-lagged panel models were used to investigate the longitudinal bidirectional relationships of sFNC with internalizing and externalizing problems. Similar cross-lagged panel models were run for dFNC.

RESULTS

Small longitudinal relationships between dFNC and certain syndrome scales were observed, especially for baseline syndrome scales (i.e., rule-breaking, somatic complaints, thought problems, and attention problems) predicting connectivity changes. However, no association between any of the psychiatric problems (broadband and syndrome scales) with either measure of FNC survived correction for multiple testing.

CONCLUSION

We found no or very modest evidence for longitudinal associations between psychiatric problems with dynamic and static FNC in this population-based sample. Differences in findings may stem from the population drawn, study design, developmental timing, and sample sizes.

Joint connectivity matrix independent component analysis: Auto-linking of structural and functional connectivities

The study of human brain connectivity, including structural connectivity (SC) and functional connectivity (FC), provides insights into the neurophysiological mechanism of brain function and its relationship to human behavior and cognition. Both types of connectivity measurements provide crucial yet complementary information. However, integrating these two modalities into a single framework remains a challenge, because of the differences in their quantitative interdependencies as well as their anatomical representations due to distinctive imaging mechanisms. In this study, we introduced a new method, joint connectivity matrix independent component analysis (cmICA), which provides a data-driven parcellation and automated-linking of SC and FC information simultaneously using a joint analysis of functional magnetic resonance imaging (MRI) and diffusion-weighted MRI data. We showed that these two connectivity modalities produce common cortical segregation, though with various degrees of (dis)similarity. Moreover, we show conjoint FC networks and structural white matter tracts that directly link these cortical parcellations/sources, within one analysis. Overall, data-driven joint cmICA provides a new approach for integrating or fusing structural connectivity and FC systematically and conveniently, and provides an effective tool for connectivity-based multimodal data fusion in brain.

Multimodal fusion of multiple rest fMRI networks and MRI gray matter via multilink joint ICA reveals highly significant function/structure coupling in Alzheimer's disease

In this paper we focus on estimating the joint relationship between structural MRI (sMRI) gray matter (GM) and multiple functional MRI (fMRI) intrinsic connectivity networks (ICN) using a novel approach called multi-link joint independent component analysis (ml-jICA). The proposed model offers several improvements over the existing joint independent component analysis (jICA) model. We assume a shared mixing matrix for both the sMRI and fMRI modalities, while allowing for different mixing matrices linking the sMRI data to the different ICNs. We introduce the model and then apply this approach to study the differences in resting fMRI and sMRI data from patients with Alzheimer's disease (AD) versus controls. The results yield significant differences with large effect sizes that include regions in overlapping portions of default mode network, and also hippocampus and thalamus. Importantly, we identify two joint components with partially overlapping regions which show opposite effects for Alzheimer's disease versus controls, but were able to be separated due to being linked to distinct functional and structural patterns. This highlights the unique strength of our approach and multimodal fusion approaches generally in revealing potentially biomarkers of brain disorders that would likely be missed by a unimodal approach. These results represent the first work linking multiple fMRI ICNs to gray matter components within a multimodal data fusion model and challenges the typical view that brain structure is more sensitive to AD than fMRI.

Determining four confounding factors in individual cognitive traits prediction with functional connectivity: an exploratory study

Resting-state functional connectivity (RSFC) has been widely adopted for individualized trait prediction. However, multiple confounding factors may impact the predicted brain-behavior relationships. In this study, we investigated the impact of 4 confounding factors including time series length, functional connectivity (FC) type, brain parcellation choice, and variance of the predicted target. The data from Human Connectome Project including 1,206 healthy subjects were employed, with 3 cognitive traits including fluid intelligence, working memory, and picture vocabulary ability as the prediction targets. We compared the prediction performance under different settings of these 4 factors using partial least square regression. Results demonstrated appropriate time series length (300 time points) and brain parcellation (independent component analysis, ICA100/200) can achieve better prediction performance without too much time consumption. FC calculated by Pearson, Spearman, and Partial correlation achieves higher accuracy and lower time cost than mutual information and coherence. Cognitive traits with larger variance among subjects can be better predicted due to the well elaboration of individual variability. In addition, the beneficial effects of increasing scan duration to prediction partially arise from the improved test-retest reliability of RSFC. Taken together, the study highlights the importance of determining these factors in RSFC-based prediction, which can facilitate standardization of RSFC-based prediction pipelines going forward.

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.

Adaptive working memory training does not produce transfer effects in cognition and neuroimaging

Despite growing interest in cognitive interventions from academia and industry, it remains unclear if working memory (WM) training, one of the most popular cognitive interventions, produces transfer effects. Transfer effects are training-induced gains in performance in untrained cognitive tasks, while practice effects are improvements in trained task. The goal of this study was to evaluate potential transfer effects by comprehensive cognitive testing and neuroimaging. In this prospective, randomized-controlled, and single-blind study, we administered an 8-week n-back training to 55 healthy middle-aged (50-64 years) participants. State-of-the-art multimodal neuroimaging was used to examine potential anatomic and functional changes. Relative to control subjects, who performed non-adaptive WM training, no near or far transfer effects were detected in experimental subjects, who performed adaptive WM training. Equivalently, no training-related changes were observed in white matter integrity, amplitude of low frequency fluctuations, glucose metabolism, functional and metabolic connectivity. Exploratory within-group comparisons revealed some gains in transfer tasks, which, however, cannot be attributed to an increased WM capacity. In conclusion, WM training produces transfer effects neither at the cognitive level nor in terms of neural structure or function. These results speak against a common view that training-related gains reflect an increase in underlying WM capacity. Instead, the presently observed practice effects may be a result of optimized task processing strategies, which do not necessarily engage neural plasticity.

A precise language network revealed by the independent component-based lesion mapping in post-stroke aphasia

Brain lesion mapping studies have provided the strongest evidence regarding the neural basis of cognition. However, it remained a problem to identify symptom-specific brain networks accounting for observed clinical and neuroanatomical heterogeneity. Independent component analysis (ICA) is a statistical method that decomposes mixed signals into multiple independent components. We aimed to solve this issue by proposing an independent component-based lesion mapping (ICLM) method to identify the language network in patients with moderate to severe post-stroke aphasia. Lesions were first extracted from 49 patients with post-stroke aphasia as masks applied to fMRI data in a cohort of healthy participants to calculate the functional connectivity (FC) within the masks and non-mask brain voxels. ICA was further performed on a reformatted FC matrix to extract multiple independent networks. Specifically, we found that one of the lesion-related independent components (ICs) highly resembled classical language networks. Moreover, the damaged level within the language-related lesioned network is strongly associated with language deficits, including aphasia quotient, naming, and auditory comprehension scores. In comparison, none of the other two traditional lesion mapping methods found any regions responsible for language dysfunction. The language-related lesioned network extracted with the ICLM method showed high specificity in detecting aphasia symptoms compared with the performance of resting ICs and classical language networks. In total, we detected a precise language network in patients with aphasia and proved its efficiency in the relationship with language symptoms. In general, our ICLM could successfully identify multiple lesion-related networks from complicated brain diseases, and be used as an effective tool to study brain-behavior relationships and provide potential biomarkers of particular clinical behavioral deficits.

Constrained CPD of Complex-Valued Multi-Subject fMRI Data via Alternating Rank-R and Rank-1 Least Squares

Complex-valued shift-invariant canonical polyadic decomposition (CPD) under a spatial phase sparsity constraint (pcsCPD) shows excellent separation performance when applied to band-pass filtered complex-valued multi-subject fMRI data. However, some useful information may also be eliminated when using a band-pass filter to suppress unwanted noise. As such, we propose an alternating rank- R and rank-1 least squares optimization to relax the CPD model. Based upon this optimization method, we present a novel constrained CPD algorithm with temporal shift-invariance and spatial sparsity and orthonormality constraints. More specifically, four steps are conducted until convergence for each iteration of the proposed algorithm: 1) use rank- R least-squares fit under spatial phase sparsity constraint to update shared spatial maps after phase de-ambiguity; 2) use orthonormality constraint to minimize the cross-talk between shared spatial maps; 3) update the aggregating mixing matrix using rank- R least-squares fit; 4) utilize shift-invariant rank-1 least-squares on a series of rank-1 matrices reconstructed by each column of the aggregating mixing matrix to update shared time courses, and subject-specific time delays and intensities. The experimental results of simulated and actual complex-valued fMRI data show that the proposed algorithm improves the estimates for task-related sensorimotor and auditory networks, compared to pcsCPD and tensorial spatial ICA. The proposed alternating rank- R and rank-1 least squares optimization is also flexible to improve CPD-related algorithm using alternating least squares.

Multimodel Order Independent Component Analysis: A Data-Driven Method for Evaluating Brain Functional Network Connectivity Within and Between Multiple Spatial Scales

Background: While functional connectivity is widely studied, there has been little work studying functional connectivity at different spatial scales. Likewise, the relationship of functional connectivity between spatial scales is unknown. Methods: We proposed an independent component analysis (ICA)-based approach to capture information at multiple-model orders (component numbers), and to evaluate functional network connectivity (FNC) both within and between model orders. We evaluated the approach by studying group differences in the context of a study of resting-state functional magnetic resonance imaging (rsfMRI) data collected from schizophrenia (SZ) individuals and healthy controls (HC). The predictive ability of FNC at multiple spatial scales was assessed using support vector machine-based classification. Results: In addition to consistent predictive patterns at both multiple-model orders and single-model orders, unique predictive information was seen at multiple-model orders and in the interaction between model orders. We observed that the FNC between model orders 25 and 50 maintained the highest predictive information between HC and SZ. Results highlighted the predictive ability of the somatomotor and visual domains both within and between model orders compared with other functional domains. Also, subcortical-somatomotor, temporal-somatomotor, and temporal-subcortical FNCs had relatively high weights in predicting SZ. Conclusions: In sum, multimodel order ICA provides a more comprehensive way to study FNC, produces meaningful and interesting results, which are applicable to future studies. We shared the spatial templates from this work at different model orders to provide a reference for the community, which can be leveraged in regression-based or fully automated (spatially constrained) ICA approaches. Impact statement Multimodel order independent component analysis (ICA) provides a comprehensive way to study brain functional network connectivity within and between multiple spatial scales, highlighting findings that would have been ignored in single-model order analysis. This work expands upon and adds to the relatively new literature on resting functional magnetic resonance imaging-based classification and prediction. Results highlighted the differentiating power of specific intrinsic connectivity networks on classifying brain disorders of schizophrenia patients and healthy participants, at different spatial scales. The spatial templates from this work provide a reference for the community, which can be leveraged in regression-based or fully automated ICA approaches.