Comparison Between Gradients and Parcellations for Functional Connectivity Prediction of Behavior

Heritability of functional gradients in the human subcortico-cortical connectivity

The human subcortex plays a pivotal role in cognition and is widely implicated in the pathophysiology of many psychiatric disorders. However, the heritability of functional gradients based on subcortico-cortical functional connectivity remains elusive. Here, leveraging twin functional MRI (fMRI) data from both the Human Connectome Project (n = 1023) and the Adolescent Brain Cognitive Development study (n = 936) datasets, we construct large-scale subcortical functional gradients and delineate an increased principal functional gradient pattern from unimodal sensory/motor networks to transmodal association networks. We observed that this principal functional gradient is heritable, and the strength of heritability exhibits a heterogeneous pattern along a hierarchical unimodal-transmodal axis in subcortex for both young adults and children. Furthermore, employing a machine learning framework, we show that this heterogeneous pattern of the principal functional gradient in subcortex can accurately discern the relationship between monozygotic twin pairs and dizygotic twin pairs with an accuracy of 76.2% (P < 0.001). The heritability of functional gradients is associated with the anatomical myelin proxied by MRI-derived T1-weighted/T2-weighted (T1w/T2w) ratio mapping in subcortex. This study provides new insights into the biological basis of subcortical functional hierarchy by revealing the structural and genetic properties of the subcortical functional gradients.

Tau follows principal axes of functional and structural brain organization in Alzheimer's disease

Alzheimer's disease (AD) is a brain network disorder where pathological proteins accumulate through networks and drive cognitive decline. Yet, the role of network connectivity in facilitating this accumulation remains unclear. Using in-vivo multimodal imaging, we show that the distribution of tau and reactive microglia in humans follows spatial patterns of connectivity variation, the so-called gradients of brain organization. Notably, less distinct connectivity patterns ("gradient contraction") are associated with cognitive decline in regions with greater tau, suggesting an interaction between reduced network differentiation and tau on cognition. Furthermore, by modeling tau in subject-specific gradient space, we demonstrate that tau accumulation in the frontoparietal and temporo-occipital cortices is associated with greater baseline tau within their functionally and structurally connected hubs, respectively. Our work unveils a role for both functional and structural brain organization in pathology accumulation in AD, and supports subject-specific gradient space as a promising tool to map disease progression.

Fast connectivity gradient approximation: maintaining spatially fine-grained connectivity gradients while reducing computational costs

Brain connectome analysis suffers from the high dimensionality of connectivity data, often forcing a reduced representation of the brain at a lower spatial resolution or parcellation. This is particularly true for graph-based representations, which are increasingly used to characterize connectivity gradients, capturing patterns of systematic spatial variation in the functional connectivity structure. However, maintaining a high spatial resolution is crucial for enabling fine-grained topographical analysis and preserving subtle individual differences that might otherwise be lost. Here we introduce a computationally efficient approach to establish spatially fine-grained connectivity gradients. At its core, it leverages a set of landmarks to approximate the underlying connectivity structure at the full spatial resolution without requiring a full-scale vertex-by-vertex connectivity matrix. We show that this approach reduces computational time and memory usage while preserving informative individual features and demonstrate its application in improving brain-behavior predictions. Overall, its efficiency can remove computational barriers and enable the widespread application of connectivity gradients to capture spatial signatures of the connectome. Importantly, maintaining a spatially fine-grained resolution facilitates to characterize the spatial transitions inherent in the core concept of gradients of brain organization.

Two common and distinct forms of variation in human functional brain networks

The cortex has a characteristic layout with specialized functional areas forming distributed large-scale networks. However, substantial work shows striking variation in this organization across people, which relates to differences in behavior. While most previous work treats individual differences as linked to boundary shifts between the borders of regions, here we show that cortical 'variants' also occur at a distance from their typical position, forming ectopic intrusions. Both 'border' and 'ectopic' variants are common across individuals, but differ in their location, network associations, properties of subgroups of individuals, activations during tasks, and prediction of behavioral phenotypes. Border variants also track significantly more with shared genetics than ectopic variants, suggesting a closer link between ectopic variants and environmental influences. This work argues that these two dissociable forms of variation-border shifts and ectopic intrusions-must be separately accounted for in the analysis of individual differences in cortical systems across people.

Functional gradients reveal cortical hierarchy changes in multiple sclerosis

Functional gradient (FG) analysis represents an increasingly popular methodological perspective for investigating brain hierarchical organization but whether and how network hierarchy changes concomitant with functional connectivity alterations in multiple sclerosis (MS) has remained elusive. Here, we analyzed FG components to uncover possible alterations in cortical hierarchy using resting-state functional MRI (rs-fMRI) data acquired in 122 MS patients and 97 healthy control (HC) subjects. Cortical hierarchy was assessed by deriving regional FG scores from rs-fMRI connectivity matrices using a functional parcellation of the cerebral cortex. The FG analysis identified a primary (visual-to-sensorimotor) and a secondary (sensory-to-transmodal) component. Results showed a significant alteration in cortical hierarchy as indexed by regional changes in FG scores in MS patients within the sensorimotor network and a compression (i.e., a reduced standard deviation across all cortical parcels) of the sensory-transmodal gradient axis, suggesting disrupted segregation between sensory and cognitive processing. Moreover, FG scores within limbic and default mode networks were significantly correlated (ρ = 0.30 $$ \rho =0.30 $$ , p < .005 after Bonferroni correction for both) with the symbol digit modality test (SDMT) score, a measure of information processing speed commonly used in MS neuropsychological assessments. Finally, leveraging supervised machine learning, we tested the predictive value of network-level FG features, highlighting the prominent role of the FG scores within the default mode network in the accurate prediction of SDMT scores in MS patients (average mean absolute error of 1.22 ± 0.07 points on a hold-out set of 24 patients). Our work provides a comprehensive evaluation of FG alterations in MS, shedding light on the hierarchical organization of the MS brain and suggesting that FG connectivity analysis can be regarded as a valuable approach in rs-fMRI studies across different MS populations.

Unveiling the Core Functional Networks of Cognition: An Ontology-Guided Machine Learning Approach

Deciphering the functional architecture that underpins diverse cognitive functions is fundamental quest in neuroscience. In this study, we employed an innovative machine learning framework that integrated cognitive ontology with functional connectivity analysis to identify brain networks essential for cognition. We identified a core assembly of functional connectomes, primarily located within the association cortex, which showed superior predictive performance compared to two conventional methods widely employed in previous research across various cognitive domains. Our approach achieved a mean prediction accuracy of 0.13 across 16 cognitive tasks, including working memory, reading comprehension, and sustained attention, outperforming the traditional methods' accuracy of 0.08. In contrast, our method showed limited predictive power for sensory, motor, and emotional functions, with a mean prediction accuracy of 0.03 across 9 relevant tasks, slightly lower than the traditional methods' accuracy of 0.04. These cognitive connectomes were further characterized by distinctive patterns of resting-state functional connectivity, structural connectivity via white matter tracts, and gene expression, highlighting their neurogenetic underpinnings. Our findings reveal a domain-general functional network fingerprint that pivotal to cognition, offering a novel computational approach to explore the neural foundations of cognitive abilities.

Dissociation of Reliability, Heritability, and Predictivity in Coarse- and Fine-Scale Functional Connectomes during Development

The functional connectome supports information transmission through the brain at various spatial scales, from exchange between broad cortical regions to finer-scale, vertex-wise connections that underlie specific information processing mechanisms. In adults, while both the coarse- and fine-scale functional connectomes predict cognition, the fine scale can predict up to twice the variance as the coarse-scale functional connectome. Yet, past brain-wide association studies, particularly using large developmental samples, focus on the coarse connectome to understand the neural underpinnings of individual differences in cognition. Using a large cohort of children (age 9-10 years; n = 1,115 individuals; both sexes; 50% female, including 170 monozygotic and 219 dizygotic twin pairs and 337 unrelated individuals), we examine the reliability, heritability, and behavioral relevance of resting-state functional connectivity computed at different spatial scales. We use connectivity hyperalignment to improve access to reliable fine-scale (vertex-wise) connectivity information and compare the fine-scale connectome with the traditional parcel-wise (coarse scale) functional connectomes. Though individual differences in the fine-scale connectome are more reliable than those in the coarse-scale, they are less heritable. Further, the alignment and scale of connectomes influence their ability to predict behavior, whereby some cognitive traits are equally well predicted by both connectome scales, but other, less heritable cognitive traits are better predicted by the fine-scale connectome. Together, our findings suggest there are dissociable individual differences in information processing represented at different scales of the functional connectome which, in turn, have distinct implications for heritability and cognition.

Should one go for individual- or group-level brain parcellations? A deep-phenotyping benchmark

The analysis and understanding of brain characteristics often require considering region-level information rather than voxel-sampled data. Subject-specific parcellations have been put forward in recent years, as they can adapt to individual brain organization and thus offer more accurate individual summaries than standard atlases. However, the price to pay for adaptability is the lack of group-level consistency of the data representation. Here, we investigate whether the good representations brought by individualized models are merely an effect of circular analysis, in which individual brain features are better represented by subject-specific summaries, or whether this carries over to new individuals, i.e., whether one can actually adapt an existing parcellation to new individuals and still obtain good summaries in these individuals. For this, we adapt a dictionary-learning method to produce brain parcellations. We use it on a deep-phenotyping dataset to assess quantitatively the patterns of activity obtained under naturalistic and controlled-task-based settings. We show that the benefits of individual parcellations are substantial, but that they vary a lot across brain systems.

Multilayer meta-matching: translating phenotypic prediction models from multiple datasets to small data

Resting-state functional connectivity (RSFC) is widely used to predict phenotypic traits in individuals. Large sample sizes can significantly improve prediction accuracies. However, for studies of certain clinical populations or focused neuroscience inquiries, small-scale datasets often remain a necessity. We have previously proposed a "meta-matching" approach to translate prediction models from large datasets to predict new phenotypes in small datasets. We demonstrated large improvement of meta-matching over classical kernel ridge regression (KRR) when translating models from a single source dataset (UK Biobank) to the Human Connectome Project Young Adults (HCP-YA) dataset. In the current study, we propose two meta-matching variants ("meta-matching with dataset stacking" and "multilayer meta-matching") to translate models from multiple source datasets across disparate sample sizes to predict new phenotypes in small target datasets. We evaluate both approaches by translating models trained from five source datasets (with sample sizes ranging from 862 participants to 36,834 participants) to predict phenotypes in the HCP-YA and HCP-Aging datasets. We find that multilayer meta-matching modestly outperforms meta-matching with dataset stacking. Both meta-matching variants perform better than the original "meta-matching with stacking" approach trained only on the UK Biobank. All meta-matching variants outperform classical KRR and transfer learning by a large margin. In fact, KRR is better than classical transfer learning when less than 50 participants are available for finetuning, suggesting the difficulty of classical transfer learning in the very small sample regime. The multilayer meta-matching model is publicly available at GITHUB_LINK.