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

Tell me why: A scoping review on the fundamental building blocks of fMRI-based network analysis

Understanding complex brain-behaviour relationships in psychiatric and neurological conditions is crucial for advancing clinical insights. This review explores the current landscape of network estimation methods in the context of functional MRI (fMRI) based network neuroscience, focusing on static undirected network analysis. We focused on papers published in a single year (2022) and characterised what we consider the fundamental building blocks of network analysis: sample size, network size, association type, edge inclusion strategy, edge weights, modelling level, and confounding factors. We found that the most common methods across all included studies (n = 191) were the use of pairwise correlations to estimate the associations between brain regions (79.6 %), estimation of weighted networks (95.3 %), and estimation of the network at the individual level (86.9 %). Importantly, a substantial number of studies lacked comprehensive reporting on their methodological choices, hindering the synthesis of research findings within the field. This review underscores the critical need for careful consideration and transparent reporting of fMRI network estimation methodologies to advance our understanding of complex brain-behaviour relationships. By facilitating the integration between network neuroscience and network psychometrics, we aim to significantly enhance our clinical understanding of these intricate connections.

Mapping the Psychosis Spectrum - Imaging Neurosubtypes from Multi-Scale Functional Network Connectivity

This study aims to identify Psychosis Imaging Neurosubtypes (PINs)-homogeneous subgroups of individuals with psychosis characterized by distinct neurobiology derived from imaging features. Specifically, we utilized resting-state fMRI data from 2103 B-SNIP 1&2 participants (1127 with psychosis, 350 relatives, 626 controls) to compute subject-specific multiscale functional network connectivity (msFNC). We then derived a low-dimensional neurobiological subspace, termed Latent Network Connectivity (LNC), which captured system-wide interconnected multiscale information across three components (cognitive-related, typical, psychosis-related). Projections of psychosis participants' msFNC onto this subspace revealed three PINs through unsupervised learning, each with distinct cognitive, clinical, and connectivity profiles, spanning all DSM diagnoses (Schizophrenia, Bipolar, Schizoaffective). PIN-1, the most cognitively impaired, showed Cerebellar-Subcortical and Visual-Sensorimotor hypoconnectivity, alongside Visual-Subcortical hyperconnectivity. Most cognitively preserved PIN-2 showed Visual-Subcortical, Subcortical-Sensorimotor, and Subcortical-Higher Cognition hypoconnectivity. PIN-3 exhibited intermediate cognitive function, showing Cerebellar-Subcortical hypoconnectivity alongside Cerebellar-Sensorimotor and Subcortical-Sensorimotor hyperconnectivity. Notably, 55% of relatives aligned with the same neurosubtype as their affected family members-a significantly higher rate than random chance (p-valueRelatives-to-PIN-1 < 0.001, p-valueRelatives-to-PIN-2 < 0.05, p-valueRelatives-to-PIN-3 < 0.001) compared to a non-significant 37% DSM-based classification, supporting a biological basis of these neurosubtypes. Cognitive performance reliably aligns with distinct brain connectivity patterns, which are also evident in relatives, supporting their construct validity. Our PINs differed from original B-SNIP Biotypes, which were determined from electrophysiological, cognitive, and oculomotor data. These findings underscore the limitations of DSM-based classifications in capturing the biological complexity of psychotic disorders and highlight the potential of imaging-based neurosubtypes to enhance our understanding of the psychosis spectrum.

Neurobiology-based cognitive biotypes using multi-scale intrinsic connectivity networks in psychotic disorders

Understanding neurobiology and developing effective interventions for cognitive dysfunction in psychotic disorders remain elusive. Insufficient knowledge about the biological heterogeneity of cognitive dysfunction hinders progress. We aimed to identify subgroups of patients with psychosis and distinct patterns of functional brain alterations related to cognition (cognitive biotypes). We analyzed B-SNIP consortium data (2 270 participants including participants with psychotic disorders, relatives, and controls, 55% females). We used reference-informed independent component analysis with the standardized and fully automated framework NeuroMark and the 100k multi-scale intrinsic connectivity networks (ICN) template to obtain subject-specific ICNs and whole-brain functional network connectivity (FNC). FNC features associated with cognitive performance were identified using multivariate joint analysis. K-means clustering identified patient subgroups based on these features. Two biotypes with different functional brain alteration patterns were identified. Relative to controls, biotype 1 exhibited hypoconnectivity in cerebellar-subcortical and somatomotor-visual networks and worse cognitive performance. Biotype 2 exhibited hyperconnectivity in somatomotor-subcortical networks, hypoconnectivity in somatomotor-high cognitive processing networks, and better-preserved cognitive performance. Demographic, clinical, cognitive, and FNC characteristics of biotypes were consistent in discovery and replication sets and in relatives. 76.56% of relatives were assigned to a psychosis biotype, of those, 70.12% were to the same biotype as their affected family members. These findings suggest two distinctive psychosis-related cognitive biotypes with differing functional brain patterns shared with their relatives. Instead of traditional diagnosis, patient stratification based on these biotypes may help optimize future research and identify biological targets for the treatment of cognitive dysfunction in psychosis.

CNN AND LSTM MODELS FOR FMRI-BASED SCHIZOPHRENIA CLASSIFICATION USING C-ICA OF DFNC

Resting-state fMRI (rs-fMRI) captures brain activity at rest, it demonstrates information on how different regions interact without explicity task-based influences. This provides insights into both healthy and disordered brain states. However, clinical application of rs-fMRI remains challenging due to the wide variability in functional connectivity across individuals. Traditional data-driven methods like independent component analysis (ICA) struggle to balance these individual differences with broader patterns. Constrained methods, such as constrained ICA (cICA), have been introduced to address this by integrating templates from multiple external datasets to enhance accuracy and consistency. In our study, we analyzed rs-fMRI data from 100,517 individuals from diverse datasets, processed through a robust quality-control dynamic connectivity pipeline established in previous work. Using the resulting brain state templates as cICA priors, we examined the effectiveness of cICA for schizophrenia classification using a combined CNN and LSTM architecture. Results showed stable classification accuracy (87.6% to 86.43%) for the CNN model, while the LSTM model performed less optimally, likely due to sequence processing, yet still yielded comparable results. These findings underscore the potential of group-informed methods and prior data templates in constrained dynamic ICA, offering improved reliability and clinical relevance in rs-fMRI analysis and advancing our understanding of brain function.

A telescopic independent component analysis on functional magnetic resonance imaging dataset

Brain function can be modeled as 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 independent component analysis (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 the DMN by evaluating the difference between healthy people and individuals with schizophrenia. We show that the TICA approach can detect the spatial hierarchy of the DMN, VN, and RFPN. In addition, the 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.

Topologically Optimized Intrinsic Brain Networks

The estimation of brain networks is instrumental in quantifying and evaluating brain function. Nevertheless, achieving precise estimations of subject-level networks has proven to be a formidable task. In response to this challenge, researchers have developed group-inference frameworks that leverage robust group-level estimations as a common reference point to infer corresponding subject-level networks. Generally, existing approaches either leverage the common reference as a strict, voxel-wise spatial constraint (i.e., strong constraints at the voxel level) or impose no constraints. Here, we propose a targeted approach that harnesses the topological information of group-level networks to encode a high-level representation of spatial properties to be used as constraints, which we refer to as Topologically Optimized Intrinsic Brain Networks (TOIBN). Consequently, our method inherits the significant advantages of constraint-based approaches, such as enhancing estimation efficacy in noisy data or small sample sizes. On the other hand, our method provides a softer constraint than voxel-wise penalties, which can result in the loss of individual variation, increased susceptibility to model biases, and potentially missing important subject-specific information. Our analyses show that the subject maps from our method are less noisy and true to the group networks while promoting subject variability that can be lost from strict constraints. We also find that the topological properties resulting from the TOIBN maps are more expressive of differences between individuals with schizophrenia and controls in the default mode, subcortical, and visual networks.

Building Multivariate Molecular Imaging Brain Atlases Using the NeuroMark PET Independent Component Analysis Framework

Molecular imaging analyses using positron emission tomography (PET) data often rely on macro-anatomical regions of interest (ROI), which may not align with chemo-architectural boundaries and obscure functional distinctions. While methods such as independent component analysis (ICA) have been useful to address this limitation, the fully data-driven nature can make it challenging to compare results across studies. Here, we introduce the NeuroMark PET approach, utilizing spatially constrained independent component analysis to define overlapping regions that may reflect the brain's molecular architecture. We first generate an ICA template for the PET radiotracer florbetapir (FBP), targeting amyloid-β (Aβ) accumulation in the brain, using blind ICA on large datasets to identify replicable independent components. Only components that targeted Aβ were included in this study, defined as Aβ networks (AβNs), by omitting components targeting myelin or other non-Aβ targets. Next, we use the AβNs as priors for spatially constrained ICA, resulting in a fully automated ICA pipeline called NeuroMark PET. This NeuroMark pipeline, including its AβNs, was validated against a standard neuroanatomical PET atlas, using data from the Alzheimer's Disease Neuroimaging Initiative (ADNI). The study included 296 cognitively normal participants with FBP PET scans and 173 with florbetaben (FBB) PET scans, an analogue radiotracer also targeting Aβ accumulation. Our results show that NeuroMark PET captures biologically meaningful, participant-specific features, such as subject specific loading values, consistent across individuals, and also shows higher sensitivity and power for detecting age-related changes compared to traditional atlas-based ROIs. Using this framework, we also highlight some of the advantages of using ICA analysis for PET data. In this study, an AβN consists of weighted voxels and forms a pattern throughout the entire brain. For example, components may have weighted values at every voxel and can overlap with one another, enabling the separation of artifacts which may coincide with the AβNs of interest. In addition, this approach allows for the differentiation, separating white matter components, which may overlap in complex ways with the AβNs, mainly residing in the neighboring gray matter. Results also showed that the most age associated AβN (representing the cognitive control network, CC1) exhibited a stronger association with age compared with macro-anatomical regions of interest. This may suggest that each NeuroMark FBP AβN represents a spatial network following chemo-architectural uptake with greater biological relevance compared with anatomical ROIs. In summary, the proposed NeuroMark PET approach offers a fully automated framework, providing accurate and reproducible brain AβNs. This approach enhances our ability to investigate the molecular underpinnings of brain function and pathology, offering an alternative to traditional ROI-based analyses.

NEUROMARK DFNC PATTERNS: A FULLY AUTOMATED PIPELINE TO ESTIMATE SUBJECT-SPECIFIC STATES FROM RS-FMRI DATA VIA CONSTRAINED ICA OF DFNC IN +100K SUBJECTS

Resting-state functional MRI (rs-fMRI) provides valuable insights into brain function during rest, but faces challenges in clinical applications due to individual differences in functional connectivity. While Independent Component Analysis (ICA) is commonly used, it struggles to balance individual variations with inter-subject information. To address this, constrained ICA (cICA) approaches have been developed using templates from multiple datasets to improve accuracy and comparability. In this study, we collected rs-fMRI data from 100,517 individuals across diverse datasets. Data were preprocessed through a standard fMRI pipeline. Our method first used replicable fMRI component templates as priors in constrained ICA (the NeuroMark pipeline), then estimated dynamic functional network connectivity (dFNC). Through clustering analysis, we generated replicable dFNC states, which were then used as priors in constrained ICA to automatically estimate subject-specific states from new subjects.This approach provides a robust framework for analyzing individual rs-fMRI data while maintaining consistency across large datasets, potentially advancing clinical applications of rs-fMRI.

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 SZ patients and 160 demographically matched healthy controls (HC). Our assessment encompassed both static and dynamic ICE, revealing significant differences in the heterogeneity of connectivity levels across available functional brain networks between SZ patients and 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 HC. 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 in HC, the brain primarily communicates through complex, less structured connectivity patterns, with occasional transitions into more focused patterns. Individuals with SZ are significantly less likely to attain these more focused and structured transient connectivity patterns. The proposed ICE measure presents a novel framework for gaining deeper insight into mechanisms of healthy and diseased brain states and represents a useful step forward in developing advanced methods to help diagnose mental health conditions.

Integrating fMRI spatial network dynamics and EEG spectral power: insights into resting state connectivity

Introduction

The Integration of functional magnetic resonance imaging (fMRI) and electroencephalography (EEG) has allowed for a novel exploration of the brain's spatial-temporal resolution. While functional brain networks show variations in both spatial and temporal dimensions, most studies focus on fixed spatial networks that change together over time.

Methods

In this study, for the first time, we link spatially dynamic brain networks with EEG spectral properties recorded simultaneously, which allows us to concurrently capture high spatial and temporal resolutions offered by these complementary imaging modalities. We estimated time-resolved brain networks using sliding window-based spatially constrained independent component analysis (scICA), producing resting brain networks that evolved over time at the voxel level. Next, we assessed their coupling with four time-varying EEG spectral power (delta, theta, alpha, and beta).

Results

Our analysis demonstrated how the networks' volumes and their voxel-level activities vary over time and revealed significant correlations with time-varying EEG spectral power. For instance, we found a strong association between increasing volume of the primary visual network and alpha band power, consistent with our hypothesis for eyes open resting state scan. Similarly, the alpha, theta, and delta power of the Pz electrode were localized to voxel-level activities of primary visual, cerebellum, and temporal networks, respectively. We also identified a strong correlation between the primary motor network and alpha (mu rhythm) and beta activity. This is consistent with motor tasks during rest, though this remains to be tested directly.

Discussion

These association between space and frequency observed during rest offer insights into the brain's spatial-temporal characteristics and enhance our understanding of both spatially varying fMRI networks and EEG band power.

Spatiotemporal Complexity in the Psychotic Brain

Psychotic disorders, such as schizophrenia and bipolar disorder, pose significant diagnostic challenges with major implications on mental health. The measures of resting-state fMRI spatiotemporal complexity offer a powerful tool for identifying irregularities in brain activity. To capture global brain connectivity, we employed information-theoretic metrics, overcoming the limitations of pairwise correlation analysis approaches. This enables a more comprehensive exploration of higher-order interactions and multiscale intrinsic connectivity networks (ICNs) in the psychotic brain. In this study, we provide converging evidence suggesting that the psychotic brain exhibits states of randomness across both spatial and temporal dimensions. To further investigate these disruptions, we estimated brain network connectivity using redundancy and synergy measures, aiming to assess the integration and segregation of topological information in the psychotic brain. Our findings reveal a disruption in the balance between redundant and synergistic information, a phenomenon we term brainquake in this study, which highlights the instability and disorganization of brain networks in psychosis. Moreover, our exploration of higher-order topological functional connectivity reveals profound disruptions in brain information integration. Aberrant information interactions were observed across both cortical and subcortical ICNs. We specifically identified the most easily affected irregularities in the sensorimotor, visual, temporal, default mode, and fronto-parietal networks, as well as in the hippocampal and amygdalar regions, all of which showed disruptions. These findings underscore the severe impact of psychotic states on multiscale critical brain networks, suggesting a profound alteration in the brain's complexity and organizational states.

Addressing Inconsistency in Functional Neuroimaging: A Replicable Data-Driven Multi-Scale Functional Atlas for Canonical Brain Networks

The advent of multiple neuroimaging methodologies has greatly aided in the conceptualization of large-scale functional brain networks in the field of cognitive neuroscience. However, there is inconsistency across studies in both nomenclature and the functional entities being described. There is a need for a unifying framework that standardizes terminology across studies while also bringing analyses and results into the same reference space. Here we present a whole-brain atlas of canonical functional brain networks derived from more than 100,000 resting-state fMRI datasets. These data-driven functional networks are highly replicable across datasets and capture information from multiple spatial scales. We have organized, labeled, and described the networks with terms familiar to the fields of cognitive and affective neuroscience in order to optimize their utility in future neuroimaging analyses and enhance the accessibility of new findings. The benefits of this atlas are not limited to future template-based or reference-guided analyses, but also extend to other data-driven neuroimaging approaches across modalities, such as those using blind independent component analysis (ICA). Future studies utilizing this atlas will contribute to greater harmonization and standardization in functional neuroimaging research.

Born Connected: Do Infants Already Have Adult-Like Multi-Scale Connectivity Networks?

The human brain undergoes remarkable development with the first six postnatal months witnessing the most dramatic structural and functional changes, making this period critical for in-depth research. rsfMRI studies have identified intrinsic connectivity networks (ICNs), including the default mode network, in infants. Although early formation of these networks has been suggested, the specific identification and number of ICNs reported in infants vary widely, leading to inconclusive findings. In adults, ICNs have provided valuable insights into brain function, spanning various mental states and disorders. A recent study analyzed data from over 100,000 subjects and generated a template of 105 multi-scale ICNs enhancing replicability and generalizability across studies. Yet, the presence of these ICNs in infants has not been investigated. This study addresses this significant gap by evaluating the presence of these multi-scale ICNs in infants, offering critical insight into the early stages of brain development and establishing a baseline for longitudinal studies. To accomplish this goal, we employ two sets of analyses. First, we employ a fully data-driven approach to investigate the presence of these ICNs from infant data itself. Towards this aim, we also introduce burst independent component analysis (bICA), which provides reliable and unbiased network identification. The results reveal the presence of these multi-scale ICNs in infants, showing a high correlation with the template (rho > 0.5), highlighting the potential for longitudinal continuity in such studies. We next demonstrate that reference-informed ICA-based techniques can reliably estimate these ICNs in infants, highlighting the feasibility of leveraging the NeuroMark framework for robust brain network extraction. This approach not only enhances cross-study comparisons across lifespans but also facilitates the study of brain changes across different age ranges. In summary, our study highlights the novel discovery that the infant brain already possesses ICNs that are widely observed in older cohorts.

Networks extracted from nonlinear fMRI connectivity exhibit unique spatial variation and enhanced sensitivity to differences between individuals with schizophrenia and controls

Schizophrenia is a chronic brain disorder associated with widespread alterations in functional brain connectivity. Although data-driven approaches such as independent component analysis are often used to study how schizophrenia impacts linearly connected networks, alterations within the underlying nonlinear functional connectivity structure remain largely unknown. Here we report the analysis of networks from explicitly nonlinear functional magnetic resonance imaging connectivity in a case-control dataset. We found systematic spatial variation, with higher nonlinear weight within core regions, suggesting that linear analyses underestimate functional connectivity within network centers. We also found that a unique nonlinear network incorporating default-mode, cingulo-opercular and central executive regions exhibits hypoconnectivity in schizophrenia, indicating that typically hidden connectivity patterns may reflect inefficient network integration in psychosis. Moreover, nonlinear networks including those previously implicated in auditory, linguistic and self-referential cognition exhibit heightened statistical sensitivity to schizophrenia diagnosis, collectively underscoring the potential of our methodology to resolve complex brain phenomena and transform clinical connectivity analysis.

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.

4D dynamic spatial brain networks at rest linked to cognition show atypical variability and coupling in schizophrenia

Despite increasing interest in the dynamics of functional brain networks, most studies focus on the changing relationships over time between spatially static networks or regions. Here we propose an approach to study dynamic spatial brain networks in human resting state functional magnetic resonance imaging (rsfMRI) data and evaluate the temporal changes in the volumes of these 4D networks. Our results show significant volumetric coupling (i.e., synchronized shrinkage and growth) between networks during the scan, that we refer to as dynamic spatial network connectivity (dSNC). We find that several features of such dynamic spatial brain networks are associated with cognition, with higher dynamic variability in these networks and higher volumetric coupling between network pairs positively associated with cognitive performance. We show that these networks are modulated differently in individuals with schizophrenia versus typical controls, resulting in network growth or shrinkage, as well as altered focus of activity within a network. Schizophrenia also shows lower spatial dynamical variability in several networks, and lower volumetric coupling between pairs of networks, thus upholding the role of dynamic spatial brain networks in cognitive impairment seen in schizophrenia. Our data show evidence for the importance of studying the typically overlooked voxel-wise changes within and between brain networks.

Beyond Artifacts: Rethinking Motion-Related Signals in Resting-State fMRI Analysis

Resting-state functional magnetic resonance imaging (rsfMRI) plays a pivotal role in estimating intrinsic brain functional connectivity within healthy and clinical populations. However, the pervasive impact of head motion confounds the interpretation of rsfMRI data and is typically addressed through preprocessing without further exploration. This investigation aims to scrutinize the intricate interplay between head motion and neurobiologically relevant BOLD signal as well as its potential clinical implications. Here, we use independent component analysis (ICA) to extract large-scale brain networks from BOLD fMRI and modeled head motion time series for 508 subjects sourced from three major psychosis projects. Our approach uncovers the presence of latent network information within modeled head motion data. Moreover, we find altered functional network connectivity (FNC) between healthy controls (HC) and individuals with schizophrenia (SZ) for BOLD and motion networks, revealing that projections of BOLD time series onto network features extracted from head motion data reflect cohort-specific information. Our approach challenges conventional perspectives by treating motion-related signals not as mere noise, but as potential repositories of valuable insights into functional brain connectivity across diverse populations.

Identifying Canonical multi-scale Intrinsic Connectivity Networks in Infant resting-state fMRI and their Association with Age

Intrinsic Connectivity Networks (ICNs) reflect functional brain organization responsible for various cognitive processes, including sensory perception, motor control, memory, and attention. In this study, we used the Multivariate-Objective Optimization Independent Component Analysis with Reference (MOO-ICAR) and the NeuroMark 2.1 (adult) template to estimate subject-specific ICNs in resting-state functional magnetic resonance imaging (rsfMRI) data of infants. The NeuroMark 2.1 template contains 105 multi-scale canonical ICNs derived from 100k+ adults across multiple datasets. The multi-scale ICNs capture functional segregation across various levels of granularity across brain, revealing functional sources and their interactions. The results showed that the 105 ICNs in infants were spatially aligned with those in the template and revealed age-related distinctive patterns in static Functional Network Connectivity (sFNC), particularly in the sub-cortical and high-level cognitive domains. This study is the first to investigate the presence and development of these multi-scale ICNs in infant rsfMRI data. Our findings confirmed the presence of identifiable canonical ICNs in infants as young as six months, showcasing a strong association between these networks and age and suggesting potential biomarkers for early identification of neurodevelopmental disability.

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.

Systematic evaluation of fMRI data-processing pipelines for consistent functional connectomics

Functional interactions between brain regions can be viewed as a network, enabling neuroscientists to investigate brain function through network science. Here, we systematically evaluate 768 data-processing pipelines for network reconstruction from resting-state functional MRI, evaluating the effect of brain parcellation, connectivity definition, and global signal regression. Our criteria seek pipelines that minimise motion confounds and spurious test-retest discrepancies of network topology, while being sensitive to both inter-subject differences and experimental effects of interest. We reveal vast and systematic variability across pipelines' suitability for functional connectomics. Inappropriate choice of data-processing pipeline can produce results that are not only misleading, but systematically so, with the majority of pipelines failing at least one criterion. However, a set of optimal pipelines consistently satisfy all criteria across different datasets, spanning minutes, weeks, and months. We provide a full breakdown of each pipeline's performance across criteria and datasets, to inform future best practices in functional connectomics.

Neurobiology-based Cognitive Biotypes Using Multi-scale Intrinsic Connectivity Networks in Psychotic Disorders

Objective

Understanding the neurobiology of cognitive dysfunction in psychotic disorders remains elusive, as does developing effective interventions. Limited knowledge about the biological heterogeneity of cognitive dysfunction hinders progress. This study aimed to identify subgroups of patients with psychosis with distinct patterns of functional brain alterations related to cognition (cognitive biotypes).

Methods

B-SNIP consortium data (2,270 participants including participants with psychotic disorders, relatives, and controls) was analyzed. Researchers used reference-informed independent component analysis and the NeuroMark 100k multi-scale intrinsic connectivity networks (ICN) template to obtain subject-specific ICNs and whole-brain functional network connectivity (FNC). FNC features associated with cognitive performance were identified through multivariate joint analysis. K-means clustering identified subgroups of patients based on these features in a discovery set. Subgroups were further evaluated in a replication set and in relatives.

Results

Two biotypes with different functional brain alteration patterns were identified. Biotype 1 exhibited brain-wide alterations, involving hypoconnectivity in cerebellar-subcortical and somatomotor-visual networks and worse cognitive performance. Biotype 2 exhibited hyperconnectivity in somatomotor-subcortical networks and hypoconnectivity in somatomotor-high cognitive processing networks, and better preserved cognitive performance. Demographic, clinical, cognitive, and FNC characteristics of biotypes were consistent in discovery and replication sets, and in relatives. 70.12% of relatives belonged to the same biotype as their affected family members.

Conclusions

These findings suggest two distinctive psychosis-related cognitive biotypes with differing functional brain patterns shared with their relatives. Patient stratification based on these biotypes instead of traditional diagnosis may help to optimize future research and clinical trials addressing cognitive dysfunction in psychotic disorders.

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.

A whole-brain neuromark resting-state fMRI analysis of first-episode and early psychosis: Evidence of aberrant cortical-subcortical-cerebellar functional circuitry

Psychosis (including symptoms of delusions, hallucinations, and disorganized conduct/speech) is a main feature of schizophrenia and is frequently present in other major psychiatric illnesses. Studies in individuals with first-episode (FEP) and early psychosis (EP) have the potential to interpret aberrant connectivity associated with psychosis during a period with minimal influence from medication and other confounds. The current study uses a data-driven whole-brain approach to examine patterns of aberrant functional network connectivity (FNC) in a multi-site dataset comprising resting-state functional magnetic resonance images (rs-fMRI) from 117 individuals with FEP or EP and 130 individuals without a psychiatric disorder, as controls. Accounting for age, sex, race, head motion, and multiple imaging sites, differences in FNC were identified between psychosis and control participants in cortical (namely the inferior frontal gyrus, superior medial frontal gyrus, postcentral gyrus, supplementary motor area, posterior cingulate cortex, and superior and middle temporal gyri), subcortical (the caudate, thalamus, subthalamus, and hippocampus), and cerebellar regions. The prominent pattern of reduced cerebellar connectivity in psychosis is especially noteworthy, as most studies focus on cortical and subcortical regions, neglecting the cerebellum. The dysconnectivity reported here may indicate disruptions in cortical-subcortical-cerebellar circuitry involved in rudimentary cognitive functions which may serve as reliable correlates of psychosis.

Networks extracted from nonlinear fMRI connectivity exhibit unique spatial variation and enhanced sensitivity to differences between individuals with schizophrenia and controls

Functional magnetic resonance imaging (fMRI) studies often estimate brain intrinsic connectivity networks (ICNs) from temporal relationships between hemodynamic signals using approaches such as independent component analysis (ICA). While ICNs are thought to represent functional sources that play important roles in various psychological phenomena, current approaches have been tailored to identify ICNs that mainly reflect linear statistical relationships. However, the elements comprising neural systems often exhibit remarkably complex nonlinear interactions that may be involved in cognitive operations and altered in psychiatric conditions such as schizophrenia. Consequently, there is a need to develop methods capable of effectively capturing ICNs from measures that are sensitive to nonlinear relationships. Here, we advance a novel approach to estimate ICNs from explicitly nonlinear whole-brain functional connectivity (ENL-wFC) by transforming resting-state fMRI (rsfMRI) data into the connectivity domain, allowing us to capture unique information from distance correlation patterns that would be missed by linear whole-brain functional connectivity (LIN-wFC) analysis. Our findings provide evidence that ICNs commonly extracted from linear (LIN) relationships are also reflected in explicitly nonlinear (ENL) connectivity patterns. ENL ICN estimates exhibit higher reliability and stability, highlighting our approach's ability to effectively quantify ICNs from rsfMRI data. Additionally, we observed a consistent spatial gradient pattern between LIN and ENL ICNs with higher ENL weight in core ICN regions, suggesting that ICN function may be subserved by nonlinear processes concentrated within network centers. We also found that a uniquely identified ENL ICN distinguished individuals with schizophrenia from healthy controls while a uniquely identified LIN ICN did not, emphasizing the valuable complementary information that can be gained by incorporating measures that are sensitive to nonlinearity in future analyses. Moreover, the ENL estimates of ICNs associated with auditory, linguistic, sensorimotor, and self-referential processes exhibit heightened sensitivity towards differentiating between individuals with schizophrenia and controls compared to LIN counterparts, demonstrating the translational value of our approach and of the ENL estimates of ICNs that are frequently reported as disrupted in schizophrenia. In summary, our findings underscore the tremendous potential of connectivity domain ICA and nonlinear information in resolving complex brain phenomena and revolutionizing the landscape of clinical FC analysis.