Connectopic mapping with resting-state fMRI

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.

Shifts in structural connectome organization in the limbic and sensory systems of patients with episodic migraine

Migraine is a complex neurological condition characterized by recurrent headaches, which is often accompanied by various neurological symptoms. Magnetic resonance imaging (MRI) is a powerful tool for investigating whole-brain connectivity patterns; however, systematic assessment of structural connectome organization has rarely been performed. In the present study, we aimed to examine the changes in structural connectivity in patients with episodic migraines using diffusion MRI. First, we computed structural connectivity using diffusion MRI tractography, after which we applied dimensionality reduction techniques to the structural connectivity and generated three low-dimensional eigenvectors. We subsequently calculated the manifold eccentricity, defined as the Euclidean distance between each data point and the center of the data in the manifold space. We then compared the manifold eccentricity between patients with migraines and healthy controls, revealing significant between-group differences in the orbitofrontal cortex, temporal pole, and sensory/motor regions. Between-group differences in subcortico-cortical connectivity further revealed significant changes in the amygdala, accumbens, and caudate nuclei. Finally, supervised machine learning effectively classified patients with migraines and healthy controls using cortical and subcortical structural connectivity features, highlighting the importance of the orbitofrontal and sensory cortices, in addition to the caudate, in distinguishing between the groups. Our findings confirmed that episodic migraine is related to the structural connectome changes in the limbic and sensory systems, suggesting its potential utility as a diagnostic marker for migraine.

A shifting role of thalamocortical connectivity in the emergence of cortical functional organization

The cortical patterning principle has been a long-standing question in neuroscience, yet how this translates to macroscale functional specialization in the human brain remains largely unknown. Here we examine age-dependent differences in resting-state thalamocortical connectivity to investigate its role in the emergence of large-scale functional networks during early life, using a primarily cross-sectional but also longitudinal approach. We show that thalamocortical connectivity during infancy reflects an early differentiation of sensorimotor networks and genetically influenced axonal projection. This pattern changes in childhood, when connectivity is established with the salience network, while decoupling externally and internally oriented functional systems. A developmental simulation using generative network models corroborated these findings, demonstrating that thalamic connectivity contributes to developing key features of the mature brain, such as functional segregation and the sensory-association axis, especially across 12-18 years of age. Our study suggests that the thalamus plays an important role in functional specialization during development, with potential implications for studying conditions with compromised internal and external processing.

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.

Characterizing Large-Scale Human Circuit Development with In Vivo Neuroimaging

Large-scale coordinated patterns of neural activity are crucial for the integration of information in the human brain and to enable complex and flexible human behavior across the life span. Through recent advances in noninvasive functional magnetic resonance imaging (fMRI) methods, it is now possible to study this activity and how it emerges in the living fetal brain across the second half of human gestation. This work has demonstrated that functional activity in the fetal brain has several features in keeping with highly organized networks of activity, which are undergoing a highly programmed and rapid sequence of development before birth, in which long-range connections emerge and core features of the mature functional connectome (such as hub regions and a gradient organization) are established. In this review, the findings of these studies are summarized, their relationship to the known changes in developmental neurobiology is considered, and considerations for future work in the context of limitations to the fMRI approach are presented.

Systematically altered connectome gradient in benign childhood epilepsy with centrotemporal spikes: Potential effect on cognitive function

OBJECTIVE

Benign childhood epilepsy with centrotemporal spikes (BECTS) affects brain network hierarchy and cognitive function; however, itremainsunclearhowhierarchical changeaffectscognition in patients with BECTS. A major aim of this study was to examine changes in the macro-network function hierarchy in BECTS and its potential contribution to cognitive function.

METHODS

Overall, the study included 50 children with BECTS and 69 healthy controls. Connectome gradient analysis was used to determine the brain network hierarchy of each group. By comparing gradient scores at each voxel level and network between groups, we assessed changes in whole-brain voxel-level and network hierarchy. Functional connectivity was used to detect the functional reorganization of epilepsy caused by these abnormal brain regions based on these aberrant gradients. Lastly, we explored the relationships between the change gradient and functional connectivity values and clinical variables and further predicted the cognitive function associated with BECTS gradient changes.

RESULTS

In children with BECTS, the gradient was extended at different network and voxel levels. The gradient scores frontoparietal network was increased in the principal gradient of patients with BECTS. The left precentral gyrus (PCG) and right angular gyrus gradient scores were significantly increased in the principal gradient of children with BECTS. Moreover, in regions of the brain with abnormal principal gradients, functional connectivity was disrupted. The left PCG gradient score of children with BECTS was correlated with the verbal intelligence quotient (VIQ), and the disruption of functional connectivity in brain regions with abnormal principal gradients was closely related to cognitive function. VIQ was significantly predicted by the principal gradient map of patients.

SIGNIFICANCE

The results indicate connectome gradient disruption in children with BECTS and its relationship to cognitive function, thereby increasing our understanding of the functional connectome hierarchy and providing potential biomarkers for cognitive function of children with BECTS.

Cortical gradient perturbation in attention deficit hyperactivity disorder correlates with neurotransmitter-, cell type-specific and chromosome- transcriptomic signatures

AIMS

This study aimed to illuminate the neuropathological landscape of attention deficit hyperactivity disorder (ADHD) by a multiscale macro-micro-molecular perspective from in vivo neuroimaging data.

METHODS

The "ADHD-200 initiative" repository provided multi-site high-quality resting-state functional connectivity (rsfc-) neuroimaging for ADHD children and matched typically developing (TD) cohort. Diffusion mapping embedding model to derive the functional connectome gradient detecting biologically plausible neural pattern was built, and the multivariate partial least square method to uncover the enrichment of neurotransmitomic, cellular and chromosomal gradient-transcriptional signatures of AHBA enrichment and meta-analytic decoding.

RESULTS

Compared to TD, ADHD children presented connectopic cortical gradient perturbations in almost all the cognition-involved brain macroscale networks (all pBH <0.001), but not in the brain global topology. As an intermediate phenotypic variant, such gradient perturbation was spatially enriched into distributions of GABAA/BZ and 5-HT2A receptors (all pBH <0.01) and co-varied with genetic transcriptional expressions (e.g. DYDC2, ATOH7, all pBH <0.01), associated with phenotypic variants in episodic memory and emotional regulations. Enrichment models demonstrated such gradient-transcriptional variants indicated the risk of both cell-specific and chromosome- dysfunctions, especially in enriched expression of oligodendrocyte precursors and endothelial cells (all pperm <0.05) as well enrichment into chromosome 18, 19 and X (pperm <0.05).

CONCLUSIONS

Our findings bridged brain macroscale neuropathological patterns to microscale/cellular biological architectures for ADHD children, demonstrating the neurobiologically pathological mechanism of ADHD into the genetic and molecular variants in GABA and 5-HT systems as well brain-derived enrichment of specific cellular/chromosomal expressions.

The flattening of spacetime hierarchy of the N,N-dimethyltryptamine brain state is characterized by harmonic decomposition of spacetime (HADES) framework

The human brain is a complex system, whose activity exhibits flexible and continuous reorganization across space and time. The decomposition of whole-brain recordings into harmonic modes has revealed a repertoire of gradient-like activity patterns associated with distinct brain functions. However, the way these activity patterns are expressed over time with their changes in various brain states remains unclear. Here, we investigate healthy participants taking the serotonergic psychedelic N,N-dimethyltryptamine (DMT) with the Harmonic Decomposition of Spacetime (HADES) framework that can characterize how different harmonic modes defined in space are expressed over time. HADES demonstrates significant decreases in contributions across most low-frequency harmonic modes in the DMT-induced brain state. When normalizing the contributions by condition (DMT and non-DMT), we detect a decrease specifically in the second functional harmonic, which represents the uni- to transmodal functional hierarchy of the brain, supporting the leading hypothesis that functional hierarchy is changed in psychedelics. Moreover, HADES' dynamic spacetime measures of fractional occupancy, life time and latent space provide a precise description of the significant changes of the spacetime hierarchical organization of brain activity in the psychedelic state.

Predicting treatment outcome based on resting-state functional connectivity in internalizing mental disorders: A systematic review and meta-analysis

Predicting treatment outcome in internalizing mental disorders prior to treatment initiation is pivotal for precision mental healthcare. In this regard, resting-state functional connectivity (rs-FC) and machine learning have often shown promising prediction accuracies. This systematic review and meta-analysis evaluates these studies, considering their risk of bias through the Prediction Model Study Risk of Bias Assessment Tool (PROBAST). We examined the predictive performance of features derived from rs-FC, identified features with the highest predictive value, and assessed the employed machine learning pipelines. We searched the electronic databases Scopus, PubMed and PsycINFO on the 12th of December 2022, which resulted in 13 included studies. The mean balanced accuracy for predicting treatment outcome was 77% (95% CI: [72%- 83%]). rs-FC of the dorsolateral prefrontal cortex had high predictive value in most studies. However, a high risk of bias was identified in all studies, compromising interpretability. Methodological recommendations are provided based on a comprehensive exploration of the studies' machine learning pipelines, and potential fruitful developments are discussed.

Body mass index-dependent shifts along large-scale gradients in human cortical organization explain dietary regulatory success

Making healthy dietary choices is essential for keeping weight within a normal range. Yet many people struggle with dietary self-control despite good intentions. What distinguishes neural processing in those who succeed or fail to implement healthy eating goals? Does this vary by weight status? To examine these questions, we utilized an analytical framework of gradients that characterize systematic spatial patterns of large-scale neural activity, which have the advantage of considering the entire suite of processes subserving self-control and potential regulatory tactics at the whole-brain level. Using an established laboratory food task capturing brain responses in natural and regulatory conditions (N = 123), we demonstrate that regulatory changes of dietary brain states in the gradient space predict individual differences in dietary success. Better regulators required smaller shifts in brain states to achieve larger goal-consistent changes in dietary behaviors, pointing toward efficient network organization. This pattern was most pronounced in individuals with lower weight status (low-BMI, body mass index) but absent in high-BMI individuals. Consistent with prior work, regulatory goals increased activity in frontoparietal brain circuits. However, this shift in brain states alone did not predict variance in dietary success. Instead, regulatory success emerged from combined changes along multiple gradients, showcasing the interplay of different large-scale brain networks subserving dietary control and possible regulatory strategies. Our results provide insights into how the brain might solve the problem of dietary control: Dietary success may be easier for people who adopt modes of large-scale brain activation that do not require significant reconfigurations across contexts and goals.

Repetition Suppression Reveals Cue-Specific Spatial Representations for Landmarks and Self-Motion Cues in the Human Retrosplenial Cortex

The efficient use of various spatial cues within a setting is crucial for successful navigation. Two fundamental forms of spatial navigation, landmark-based and self-motion-based, engage distinct cognitive mechanisms. The question of whether these modes invoke shared or separate spatial representations in the brain remains unresolved. While nonhuman animal studies have yielded inconsistent results, human investigation is limited. In our previous work (Chen et al., 2019), we introduced a novel spatial navigation paradigm utilizing ultra-high field fMRI to explore neural coding of positional information. We found that different entorhinal subregions in the right hemisphere encode positional information for landmarks and self-motion cues. The present study tested the generalizability of our previous finding with a modified navigation paradigm. Although we did not replicate our previous finding in the entorhinal cortex, we identified adaptation-based allocentric positional codes for both cue types in the retrosplenial cortex (RSC), which were not confounded by the path to the spatial location. Crucially, the multi-voxel patterns of these spatial codes differed between the cue types, suggesting cue-specific positional coding. The parahippocampal cortex exhibited positional coding for self-motion cues, which was not dissociable from path length. Finally, the brain regions involved in successful navigation differed from our previous study, indicating overall distinct neural mechanisms recruited in our two studies. Taken together, the current findings demonstrate cue-specific allocentric positional coding in the human RSC in the same navigation task for the first time and that spatial representations in the brain are contingent on specific experimental conditions.

Dopamine D1-Receptor Organization Contributes to Functional Brain Architecture

Recent work has recognized a gradient-like organization in cortical function, spanning from primary sensory to transmodal cortices. It has been suggested that this axis is aligned with regional differences in neurotransmitter expression. Given the abundance of dopamine D1-receptors (D1DR), and its importance for modulation and neural gain, we tested the hypothesis that D1DR organization is aligned with functional architecture, and that inter-regional relationships in D1DR co-expression modulate functional cross talk. Using the world's largest dopamine D1DR-PET and MRI database (N = 180%, 50% female), we demonstrate that D1DR organization follows a unimodal-transmodal hierarchy, expressing a high spatial correspondence to the principal gradient of functional connectivity. We also demonstrate that individual differences in D1DR density between unimodal and transmodal regions are associated with functional differentiation of the apices in the cortical hierarchy. Finally, we show that spatial co-expression of D1DR primarily modulates couplings within, but not between, functional networks. Together, our results show that D1DR co-expression provides a biomolecular layer to the functional organization of the brain.

Whole-brain structural connectome asymmetry in autism

Autism spectrum disorder is a common neurodevelopmental condition that manifests as a disruption in sensory and social skills. Although it has been shown that the brain morphology of individuals with autism is asymmetric, how this differentially affects the structural connectome organization of each hemisphere remains under-investigated. We studied whole-brain structural connectivity-based brain asymmetry in individuals with autism using diffusion magnetic resonance imaging obtained from the Autism Brain Imaging Data Exchange initiative. By leveraging dimensionality reduction techniques, we constructed low-dimensional representations of structural connectivity and calculated their asymmetry index. Comparing the asymmetry index between individuals with autism and neurotypical controls, we found atypical structural connectome asymmetry in the sensory and default-mode regions, particularly showing weaker asymmetry towards the right hemisphere in autism. Network communication provided topological underpinnings by demonstrating that the inferior temporal cortex and limbic and frontoparietal regions showed reduced global network communication efficiency and decreased send-receive network navigation in the inferior temporal and lateral visual cortices in individuals with autism. Finally, supervised machine learning revealed that structural connectome asymmetry could be used as a measure for predicting communication-related autistic symptoms and nonverbal intelligence. Our findings provide insights into macroscale structural connectome alterations in autism and their topological underpinnings.

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.

Abnormal higher-order network interactions in Parkinson's disease visual hallucinations

Visual hallucinations in Parkinson's disease can be viewed from a systems-level perspective, whereby dysfunctional communication between brain networks responsible for perception predisposes a person to hallucinate. To this end, abnormal functional interactions between higher-order and primary sensory networks have been implicated in the pathophysiology of visual hallucinations in Parkinson's disease, however the precise signatures remain to be determined. Dimensionality reduction techniques offer a novel means for simplifying the interpretation of multidimensional brain imaging data, identifying hierarchical patterns in the data that are driven by both within- and between-functional network changes. Here, we applied two complementary non-linear dimensionality reduction techniques-diffusion-map embedding and t-distributed stochastic neighbour embedding (t-SNE)-to resting state functional MRI data, in order to characterize the altered functional hierarchy associated with susceptibility to visual hallucinations. Our study involved 77 people with Parkinson's disease (31 with hallucinations; 46 without hallucinations) and 19 age-matched healthy control subjects. In patients with visual hallucinations, we found compression of the unimodal-heteromodal gradient consistent with increased functional integration between sensory and higher order networks. This was mirrored in a traditional functional connectivity analysis, which showed increased connectivity between the visual and default mode networks in the hallucinating group. Together, these results suggest a route by which higher-order regions may have excessive influence over earlier sensory processes, as proposed by theoretical models of hallucinations across disorders. By contrast, the t-SNE analysis identified distinct alterations in prefrontal regions, suggesting an additional layer of complexity in the functional brain network abnormalities implicated in hallucinations, which was not apparent in traditional functional connectivity analyses. Together, the results confirm abnormal brain organization associated with the hallucinating phenotype in Parkinson's disease and highlight the utility of applying convergent dimensionality reduction techniques to investigate complex clinical symptoms. In addition, the patterns we describe in Parkinson's disease converge with those seen in other conditions, suggesting that reduced hierarchical differentiation across sensory-perceptual systems may be a common transdiagnostic vulnerability in neuropsychiatric disorders with perceptual disturbances.

Structural connectome alterations between individuals with autism and neurotypical controls using feature representation learning

Autism spectrum disorder is one of the most common neurodevelopmental conditions associated with sensory and social communication impairments. Previous neuroimaging studies reported that atypical nodal- or network-level functional brain organization in individuals with autism was associated with autistic behaviors. Although dimensionality reduction techniques have the potential to uncover new biomarkers, the analysis of whole-brain structural connectome abnormalities in a low-dimensional latent space is underinvestigated. In this study, we utilized autoencoder-based feature representation learning for diffusion magnetic resonance imaging-based structural connectivity in 80 individuals with autism and 61 neurotypical controls that passed strict quality controls. We generated low-dimensional latent features using the autoencoder model for each group and adopted an integrated gradient approach to assess the contribution of the input data for predicting latent features during the encoding process. Subsequently, we compared the integrated gradient values between individuals with autism and neurotypical controls and observed differences within the transmodal regions and between the sensory and limbic systems. Finally, we identified significant associations between integrated gradient values and communication abilities in individuals with autism. Our findings provide insights into the whole-brain structural connectome in autism and may help identify potential biomarkers for autistic connectopathy.

Reorganization of the cortical connectome functional gradient in age-related hearing loss

Age-related hearing loss (ARHL), one of the most common sensory deficits in elderly individuals, is a risk factor for dementia; however, it is unclear how ARHL affects the decline in cognitive function. To address this issue, a connectome gradient framework was used to identify critical features of information integration between sensory and cognitive processing centers using resting-state functional magnetic resonance imaging (rs-fMRI) data from 40 individuals with ARHL and 36 healthy controls (HCs). The first three functional gradient alterations associated with ARHL were investigated at the global, network and regional levels. Using a support vector machine (SVM) model, our analysis distinguished individuals with ARHL with normal cognitive function from those with cognitive decline. Compared to HCs, individuals with ARHL had a contracted principal primary-to-transmodal gradient axis, especially in the visual and default mode networks, with an altered gradient explained ratio and variance. Among individuals with ARHL, cognitive decline was detected in the visual network in the principal gradient as well as in the limbic, salience and default mode networks in the third gradient (salience to frontoparietal/default mode). These results suggest that ARHL is associated with disrupted information processing from the primary sensory networks to higher-order cognitive networks and highlight the key nodes closely associated with cognitive decline during cognitive processing in ARHL, providing new insights into the mechanism of cognitive impairment and suggesting potential treatments related to ARHL.

Graded functional organization in the left inferior frontal gyrus: evidence from task-free and task-based functional connectivity

The left inferior frontal gyrus has been ascribed key roles in numerous cognitive domains, such as language and executive function. However, its functional organization is unclear. Possibilities include a singular domain-general function, or multiple functions that can be mapped onto distinct subregions. Furthermore, spatial transition in function may be either abrupt or graded. The present study explored the topographical organization of the left inferior frontal gyrus using a bimodal data-driven approach. We extracted functional connectivity gradients from (i) resting-state fMRI time-series and (ii) coactivation patterns derived meta-analytically from heterogenous sets of task data. We then sought to characterize the functional connectivity differences underpinning these gradients with seed-based resting-state functional connectivity, meta-analytic coactivation modeling and functional decoding analyses. Both analytic approaches converged on graded functional connectivity changes along 2 main organizational axes. An anterior-posterior gradient shifted from being preferentially associated with high-level control networks (anterior functional connectivity) to being more tightly coupled with perceptually driven networks (posterior). A second dorsal-ventral axis was characterized by higher connectivity with domain-general control networks on one hand (dorsal functional connectivity), and with the semantic network, on the other (ventral). These results provide novel insights into an overarching graded functional organization of the functional connectivity that explains its role in multiple cognitive domains.

Early visual experience refines the retinotopic organization within and across visual cortical regions

Early visual areas are retinotopically organized in human and non-human primates. Population receptive field (pRF) size increases with eccentricity and from lower- to higher-level visual areas. Furthermore, the cortical magnification factor (CMF), a measure of how much cortical space is devoted to each degree of visual angle, is typically larger for foveal as opposed to peripheral regions of the visual field. Whether this fine-scale organization within and across visual areas depends on early visual experience has yet been unknown. Here, we employed 7T functional magnetic resonance imaging pRF mapping to assess the retinotopic organization of early visual regions (i.e., V1, V2, and V3) in eight sight recovery individuals with a history of congenital blindness until a maximum of 4 years of age. Compared with sighted controls, foveal pRF sizes in these individuals were larger, and pRF sizes did not show the typical increase with eccentricity and down the visual cortical processing stream (V1-V2-V3). Cortical magnification was overall diminished and decreased less from foveal to parafoveal visual field locations. Furthermore, cortical magnification correlated with visual acuity in sight recovery individuals. The results of this study suggest that early visual experience is essential for refining a presumably innate prototypical retinotopic organization in humans within and across visual areas, which seems to be crucial for acquiring full visual capabilities.

Probing neurodynamics of experienced emotions-a Hitchhiker's guide to film fMRI

Film functional magnetic resonance imaging (fMRI) has gained tremendous popularity in many areas of neuroscience. However, affective neuroscience remains somewhat behind in embracing this approach, even though films lend themselves to study how brain function gives rise to complex, dynamic and multivariate emotions. Here, we discuss the unique capabilities of film fMRI for emotion research, while providing a general guide of conducting such research. We first give a brief overview of emotion theories as these inform important design choices. Next, we discuss films as experimental paradigms for emotion elicitation and address the process of annotating them. We then situate film fMRI in the context of other fMRI approaches, and present an overview of results from extant studies so far with regard to advantages of film fMRI. We also give an overview of state-of-the-art analysis techniques including methods that probe neurodynamics. Finally, we convey limitations of using film fMRI to study emotion. In sum, this review offers a practitioners' guide to the emerging field of film fMRI and underscores how it can advance affective neuroscience.

Functional gradients reveal altered functional segregation in patients with amnestic mild cognitive impairment and Alzheimer's disease

Alzheimer's disease and amnestic mild cognitive impairment are associated with disrupted functional organization in brain networks, involved with alteration of functional segregation. Connectome gradients are a new tool representing brain functional topological organization to smoothly capture the human macroscale hierarchy. Here, we examined altered topological organization in amnestic mild cognitive impairment and Alzheimer's disease by connectome gradient mapping. We further quantified functional segregation by gradient dispersion. Then, we systematically compared the alterations observed in amnestic mild cognitive impairment and Alzheimer's disease patients with those in normal controls in a two-dimensional functional gradient space from both the whole-brain level and module level. Compared with normal controls, the first gradient, which described the neocortical hierarchy from unimodal to transmodal regions, showed a more distributed and significant suppression in Alzheimer's disease than amnestic mild cognitive impairment patients. Furthermore, gradient dispersion showed significant decreases in Alzheimer's disease at both the global level and module level, whereas this alteration was limited only to limbic areas in amnestic mild cognitive impairment. Notably, we demonstrated that suppressed gradient dispersion in amnestic mild cognitive impairment and Alzheimer's disease was associated with cognitive scores. These findings provide new evidence for altered brain hierarchy in amnestic mild cognitive impairment and Alzheimer's disease, which strengthens our understanding of the progressive mechanism of cognitive decline.

Evaluating functional brain organization in individuals and identifying contributions to network overlap

Individual differences in the spatial organization of resting state networks have received increased attention in recent years. Measures of individual-specific spatial organization of brain networks and overlapping network organization have been linked to important behavioral and clinical traits and are therefore potential biomarker targets for personalized psychiatry approaches. To better understand individual-specific spatial brain organization, this paper addressed three key goals. First, we determined whether it is possible to reliably estimate weighted (non-binarized) resting state network maps using data from only a single individual, while also maintaining maximum spatial correspondence across individuals. Second, we determined the degree of spatial overlap between distinct networks, using test-retest and twin data. Third, we systematically tested multiple hypotheses (spatial mixing, temporal switching, and coupling) as candidate explanations for why networks overlap spatially. To estimate weighted network organization, we adopt the Probabilistic Functional Modes (PROFUMO) algorithm, which implements a Bayesian framework with hemodynamic and connectivity priors to supplement optimization for spatial sparsity/independence. Our findings showed that replicable individual-specific estimates of weighted resting state networks can be derived using high quality fMRI data within individual subjects. Network organization estimates using only data from each individual subject closely resembled group-informed network estimates (which was not explicitly modeled in our individual-specific analyses), suggesting that cross-subject correspondence was largely maintained. Furthermore, our results confirmed the presence of spatial overlap in network organization, which was replicable across sessions within individuals and in monozygotic twin pairs. Intriguingly, our findings provide evidence that network overlap is indicative of linear additive coupling. These results suggest that regions of network overlap concurrently process information from all contributing networks, potentially pointing to the role of overlapping network organization in the integration of information across multiple brain systems.

Age-related differences in the functional topography of the locus coeruleus and their implications for cognitive and affective functions

The locus coeruleus (LC) is an important noradrenergic nucleus that has recently attracted a lot of attention because of its emerging role in cognitive and psychiatric disorders. Although previous histological studies have shown that the LC has heterogeneous connections and cellular features, no studies have yet assessed its functional topography in vivo, how this heterogeneity changes over aging, and whether it is associated with cognition and mood. Here, we employ a gradient-based approach to characterize the functional heterogeneity in the organization of the LC over aging using 3T resting-state fMRI in a population-based cohort aged from 18 to 88 years of age (Cambridge Centre for Ageing and Neuroscience cohort, n=618). We show that the LC exhibits a rostro-caudal functional gradient along its longitudinal axis, which was replicated in an independent dataset (Human Connectome Project [HCP] 7T dataset, n=184). Although the main rostro-caudal direction of this gradient was consistent across age groups, its spatial features varied with increasing age, emotional memory, and emotion regulation. More specifically, a loss of rostral-like connectivity, more clustered functional topography, and greater asymmetry between right and left LC gradients was associated with higher age and worse behavioral performance. Furthermore, participants with higher-than-normal Hospital Anxiety and Depression Scale (HADS) ratings exhibited alterations in the gradient as well, which manifested in greater asymmetry. These results provide an in vivo account of how the functional topography of the LC changes over aging, and imply that spatial features of this organization are relevant markers of LC-related behavioral measures and psychopathology.

Fine-grained topographic organization within somatosensory cortex during resting-state and emotional face-matching task and its association with ASD traits

Sensory atypicalities are particularly common in autism spectrum disorders (ASD). Nevertheless, our knowledge about the divergent functioning of the underlying somatosensory region and its association with ASD phenotype features is limited. We applied a data-driven approach to map the fine-grained variations in functional connectivity of the primary somatosensory cortex (S1) to the rest of the brain in 240 autistic and 164 neurotypical individuals from the EU-AIMS LEAP dataset, aged between 7 and 30. We estimated the S1 connection topography ('connectopy') at rest and during the emotional face-matching (Hariri) task, an established measure of emotion reactivity, and accessed its association with a set of clinical and behavioral variables. We first demonstrated that the S1 connectopy is organized along a dorsoventral axis, mapping onto the S1 somatotopic organization. We then found that its spatial characteristics were linked to the individuals' adaptive functioning skills, as measured by the Vineland Adaptive Behavior Scales, across the whole sample. Higher functional differentiation characterized the S1 connectopies of individuals with higher daily life adaptive skills. Notably, we detected significant differences between rest and the Hariri task in the S1 connectopies, as well as their projection maps onto the rest of the brain suggesting a task-modulating effect on S1 due to emotion processing. All in all, variation of adaptive skills appears to be reflected in the brain's mesoscale neural circuitry, as shown by the S1 connectivity profile, which is also differentially modulated during rest and emotional processing.

The Cellular Underpinnings of the Human Cortical Connectome

The functional properties of the human brain arise, in part, from the vast assortment of cell types that pattern the cortex. The cortical sheet can be broadly divided into distinct networks, which are further embedded into processing streams, or gradients, that extend from unimodal systems through higher-order association territories. Here, using transcriptional data from the Allen Human Brain Atlas, we demonstrate that imputed cell type distributions are spatially coupled to the functional organization of cortex, as estimated through fMRI. Cortical cellular profiles follow the macro-scale organization of the functional gradients as well as the associated large-scale networks. Distinct cellular fingerprints were evident across networks, and a classifier trained on post-mortem cell-type distributions was able to predict the functional network allegiance of cortical tissue samples. These data indicate that the in vivo organization of the cortical sheet is reflected in the spatial variability of its cellular composition.

Connectopic mapping techniques do not reflect functional gradients in the brain

Functional gradients, in which response properties change gradually across a brain region, have been proposed as a key organising principle of the brain. Recent studies using both resting-state and natural viewing paradigms have indicated that these gradients may be reconstructed from functional connectivity patterns via "connectopic mapping" analyses. However, local connectivity patterns may be confounded by spatial autocorrelations artificially introduced during data analysis, for instance by spatial smoothing or interpolation between coordinate spaces. Here, we investigate whether such confounds can produce illusory connectopic gradients. We generated datasets comprising random white noise in subjects' functional volume spaces, then optionally applied spatial smoothing and/or interpolated the data to a different volume or surface space. Both smoothing and interpolation induced spatial autocorrelations sufficient for connectopic mapping to produce both volume- and surface-based local gradients in numerous brain regions. Furthermore, these gradients appeared highly similar to those obtained from real natural viewing data, although gradients generated from real and random data were statistically different in certain scenarios. We also reconstructed global gradients across the whole-brain - while these appeared less susceptible to artificial spatial autocorrelations, the ability to reproduce previously reported gradients was closely linked to specific features of the analysis pipeline. These results indicate that previously reported gradients identified by connectopic mapping techniques may be confounded by artificial spatial autocorrelations introduced during the analysis, and in some cases may reproduce poorly across different analysis pipelines. These findings imply that connectopic gradients need to be interpreted with caution.

Age-related differences in the functional topography of the locus coeruleus: implications for cognitive and affective functions

The locus coeruleus (LC) is an important noradrenergic nucleus that has recently attracted a lot of attention because of its emerging role in cognitive and psychiatric disorders. Although previous histological studies have shown that the LC has heterogeneous connections and cellular features, no studies have yet assessed its functional topography in vivo, how this heterogeneity changes over aging and whether it is associated with cognition and mood. Here we employ a gradient-based approach to characterize the functional heterogeneity in the organization of the LC over aging using 3T resting-state fMRI in a population-based cohort aged from 18 to 88 years old (Cambridge Centre for Ageing and Neuroscience cohort, n=618). We show that the LC exhibits a rostro-caudal functional gradient along its longitudinal axis, which was replicated in an independent dataset (Human Connectome Project 7T dataset, n=184). Although the main rostro-caudal direction of this gradient was consistent across age groups, its spatial features varied with increasing age, emotional memory and emotion regulation. More specifically, a loss of rostral-like connectivity, more clustered functional topography and greater asymmetry between right and left LC gradients was associated with higher age and worse behavioral performance. Furthermore, participants with higher-than-normal Hospital Anxiety and Depression Scale ratings exhibited alterations in the gradient as well, which manifested in greater asymmetry. These results provide an in vivo account of how the functional topography of the LC changes over aging, and imply that spatial features of this organization are relevant markers of LC-related behavioral measures and psychopathology.

Geometric constraints on human brain function

The anatomy of the brain necessarily constrains its function, but precisely how remains unclear. The classical and dominant paradigm in neuroscience is that neuronal dynamics are driven by interactions between discrete, functionally specialized cell populations connected by a complex array of axonal fibres1-3. However, predictions from neural field theory, an established mathematical framework for modelling large-scale brain activity4-6, suggest that the geometry of the brain may represent a more fundamental constraint on dynamics than complex interregional connectivity7,8. Here, we confirm these theoretical predictions by analysing human magnetic resonance imaging data acquired under spontaneous and diverse task-evoked conditions. Specifically, we show that cortical and subcortical activity can be parsimoniously understood as resulting from excitations of fundamental, resonant modes of the brain's geometry (that is, its shape) rather than from modes of complex interregional connectivity, as classically assumed. We then use these geometric modes to show that task-evoked activations across over 10,000 brain maps are not confined to focal areas, as widely believed, but instead excite brain-wide modes with wavelengths spanning over 60 mm. Finally, we confirm predictions that the close link between geometry and function is explained by a dominant role for wave-like activity, showing that wave dynamics can reproduce numerous canonical spatiotemporal properties of spontaneous and evoked recordings. Our findings challenge prevailing views and identify a previously underappreciated role of geometry in shaping function, as predicted by a unifying and physically principled model of brain-wide dynamics.

Hi-Fi fMRI: High-resolution, fast-sampled and sub-second whole-brain functional MRI at 3T in humans

Functional magnetic resonance imaging (fMRI) is a methodological cornerstone of neuroscience. Most studies measure blood-oxygen-level-dependent (BOLD) signal using echo-planar imaging (EPI), Cartesian sampling, and image reconstruction with a one-to-one correspondence between the number of acquired volumes and reconstructed images. However, EPI schemes are subject to trade-offs between spatial and temporal resolutions. We overcome these limitations by measuring BOLD with a gradient recalled echo (GRE) with 3D radial-spiral phyllotaxis trajectory at a high sampling rate (28.24ms) on standard 3T field-strength. The framework enables the reconstruction of 3D signal time courses with whole-brain coverage at simultaneously higher spatial (1mm ³ ) and temporal (up to 250ms) resolutions, as compared to optimized EPI schemes. Additionally, artifacts are corrected before image reconstruction; the desired temporal resolution is chosen after scanning and without assumptions on the shape of the hemodynamic response. By showing activation in the calcarine sulcus of 20 participants performing an ON-OFF visual paradigm, we demonstrate the reliability of our method for cognitive neuroscience research.

Subcortical functional connectivity gradients in temporal lobe epilepsy

BACKGROUND AND MOTIVATION

Functional gradients have been used to study differences in connectivity between healthy and diseased brain states, however this work has largely focused on the cortex. Because the subcortex plays a key role in seizure initiation in temporal lobe epilepsy (TLE), subcortical functional-connectivity gradients may help further elucidate differences between healthy brains and TLE, as well as differences between left (L)-TLE and right (R)-TLE.

METHODS

In this work, we calculated subcortical functional-connectivity gradients (SFGs) from resting-state functional MRI (rs-fMRI) by measuring the similarity in connectivity profiles of subcortical voxels to cortical gray matter voxels. We performed this analysis in 24 R-TLE patients and 31 L-TLE patients (who were otherwise matched for age, gender, disease specific characteristics, and other clinical variables), and 16 controls. To measure differences in SFGs between L-TLE and R-TLE, we quantified deviations in the average functional gradient distributions, as well as their variance, across subcortical structures.

RESULTS

We found an expansion, measured by increased variance, in the principal SFG of TLE relative to controls. When comparing the gradient across subcortical structures between L-TLE and R-TLE, we found that abnormalities in the ipsilateral hippocampal gradient distributions were significantly different between L-TLE and R-TLE.

CONCLUSION

Our results suggest that expansion of the SFG is characteristic of TLE. Subcortical functional gradient differences exist between left and right TLE and are driven by connectivity changes in the hippocampus ipsilateral to the seizure onset zone.

Functional gradients of the medial parietal cortex in a healthy cohort with family history of sporadic Alzheimer's disease

BACKGROUND

The medial parietal cortex is an early site of pathological protein deposition in Alzheimer's disease (AD). Previous studies have identified different subregions within this area; however, these subregions are often heterogeneous and disregard individual differences or subtle pathological alterations in the underlying functional architecture. To address this limitation, here we measured the continuous connectivity gradients of the medial parietal cortex and assessed their relationship with cerebrospinal fluid (CSF) biomarkers, ApoE ε4 carriership and memory in asymptomatic individuals at risk to develop AD.

METHODS

Two hundred sixty-three cognitively normal participants with a family history of sporadic AD who underwent resting-state and task-based functional MRI using encoding and retrieval tasks were included from the PREVENT-AD cohort. A novel method for characterizing spatially continuous patterns of functional connectivity was applied to estimate functional gradients in the medial parietal cortex during the resting-state and task-based conditions. This resulted in a set of nine parameters that described the appearance of the gradient across different spatial directions. We performed correlation analyses to assess whether these parameters were associated with CSF biomarkers of phosphorylated tau181 (p-tau), total tau (t-tau), and amyloid-ß1-42 (Aß). Then, we compared the spatial parameters between ApoE ε4 carriers and noncarriers, and evaluated the relationship between these parameters and memory.

RESULTS

Alterations involving the superior part of the medial parietal cortex, which was connected to regions of the default mode network, were associated with higher p-tau, t-tau levels as well as lower Aß/p-tau levels during the resting-state condition (p < 0.01). Similar alterations were found in ApoE ε4 carriers compared to non-carriers (p < 0.003). In contrast, lower immediate memory scores were associated with changes in the middle part of the medial parietal cortex, which was connected to inferior temporal and posterior parietal regions, during the encoding task (p = 0.001). No results were found when using conventional connectivity measures.

CONCLUSIONS

Functional alterations in the medial parietal gradients are associated with CSF AD biomarkers, ApoE ε4 carriership, and lower memory in an asymptomatic cohort with a family history of sporadic AD, suggesting that functional gradients are sensitive to subtle changes associated with early AD stages.

Gradients of striatal function in antipsychotic-free first-episode psychosis and schizotypy

Both psychotic illness and subclinical psychosis-like experiences (PLEs) have been associated with cortico-striatal dysfunction. This work has largely relied on a discrete parcellation of the striatum into distinct functional areas, but recent evidence suggests that the striatum comprises multiple overlapping and smoothly varying gradients (i.e., modes) of functional organization. Here, we investigated two of these functional connectivity modes, previously associated with variations in the topographic patterning of cortico-striatal connectivity (first-order gradient), and dopaminergic innervation of the striatum (second-order gradient), and assessed continuities in striatal function from subclinical to clinical domains. We applied connectopic mapping to resting-state fMRI data to obtain the first-order and second-order striatal connectivity modes in two distinct samples: (1) 56 antipsychotic-free patients (26 females) with first-episode psychosis (FEP) and 27 healthy controls (17 females); and (2) a community-based cohort of 377 healthy individuals (213 females) comprehensively assessed for subclinical PLEs and schizotypy. The first-order "cortico-striatal" and second-order "dopaminergic" connectivity gradients were significantly different in FEP patients compared to controls bilaterally. In the independent sample of healthy individuals, variations in the left first-order "cortico-striatal" connectivity gradient were associated with inter-individual differences in a factor capturing general schizotypy and PLE severity. The presumed cortico-striatal connectivity gradient was implicated in both subclinical and clinical cohorts, suggesting that variations in its organization may represent a neurobiological trait marker across the psychosis continuum. Disruption of the presumed dopaminergic gradient was only noticeable in patients, suggesting that neurotransmitter dysfunction may be more apparent to clinical illness.

Evaluate the efficacy and reliability of functional gradients in within-subject designs

The cerebral cortex is characterized as the integration of distinct functional principles that correspond to basic primary functions, such as vision and movement, and domain-general functions, such as attention and cognition. Diffusion embedding approach is a novel tool to describe transitions between different functional principles, and has been successively applied to investigate pathological conditions in between-group designs. What still lacking and urgently needed is the efficacy of this method to differentiate within-subject circumstances. In this study, we applied the diffusion embedding to eyes closed (EC) and eyes on (EO) resting-state conditions from 145 participants. We found significantly lower within-network dispersion of visual network (VN) (p = 7.3 × 10^-4 ) as well as sensorimotor network (SMN) (p = 1 × 10^-5 ) and between-network dispersion of VN (p = 2.3 × 10^-4 ) under EC than EO, while frontoparietal network (p = 9.2 × 10^-4 ) showed significantly higher between-network dispersion during EC than EO. Test-retest reliability analysis further displayed fair reliability (intraclass correlation coefficient [ICC] < 0.4) of the network dispersions (mean ICC = 0.116 ± 0.143 [standard deviation]) except for the within-network dispersion of SMN under EO (ICC = 0.407). And the reliability under EO was higher but not significantly higher than reliability under EC. Our study demonstrated that the diffusion embedding approach that shows fair reliability is capable of distinguishing EC and EO resting-state conditions, such that this method could be generalized to other within-subject designs.

Whole-brain functional gradients reveal cortical and subcortical alterations in patients with episodic migraine

Migraine is a type of headache with multiple neurological symptoms. Prior neuroimaging studies in patients with migraine based on functional magnetic resonance imaging have found regional as well as network-level alterations in brain function. Here, we expand on prior studies by establishing whole-brain functional connectivity patterns in patients with migraine using dimensionality reduction techniques. We studied functional brain connectivity in 50 patients with episodic migraine and sex- and age-matched healthy controls. Using dimensionality reduction techniques that project high-dimensional functional connectivity onto low-dimensional representations (i.e., eigenvectors), we found significant between-group differences in the eigenvectors between patients with migraine and healthy controls, particularly in the sensory/motor and limbic cortices. Furthermore, we assessed between-group differences in subcortical connectivity with subcortical weighted manifolds defined by subcortico-cortical connectivity multiplied by cortical eigenvectors and revealed significant alterations in the amygdala. Finally, leveraging supervised machine learning, we moderately predicted headache frequency using cortical and subcortical functional connectivity features, again indicating that sensory and limbic regions play a particularly important role in predicting migraine frequency. Our study confirmed that migraine is a hierarchical disease of the brain that shows alterations along the sensory-limbic axis, and therefore, the functional connectivity in these areas could be a useful marker to investigate migraine symptomatology.

Correspondence of functional connectivity gradients across human isocortex, cerebellum, and hippocampus

Gradient mapping is an important technique to summarize high dimensional biological features as low dimensional manifold representations in exploring brain structure-function relationships at various levels of the cerebral cortex. While recent studies have characterized the major gradients of functional connectivity in several brain structures using this technique, very few have systematically examined the correspondence of such gradients across structures under a common systems-level framework. Using resting-state functional magnetic resonance imaging, here we show that the organizing principles of the isocortex, and those of the cerebellum and hippocampus in relation to the isocortex, can be described using two common functional gradients. We suggest that the similarity in functional connectivity gradients across these structures can be meaningfully interpreted within a common computational framework based on the principles of predictive processing. The present results, and the specific hypotheses that they suggest, represent an important step toward an integrative account of brain function.

Joint-embeddings reveal functional differences in default-mode network architecture between marmosets and humans

The default-mode network (DMN) is a distributed functional brain system integral for social and higher-order cognition in humans with implications in a myriad of neuropsychological disorders. In this study, we compared the functional architecture of the DMN between humans and marmosets to assess their similarities and differences using joint gradients. This approach permits simultaneous large-scale mapping of functional systems across the cortex of humans and marmosets, revealing evidence of putative homologies between them. In doing so, we find that the DMN architecture of the marmoset exhibits differences along its anterolateral-posterior axis. Specifically, the anterolateral node of the DMN (dorsolateral prefrontal cortex) displayed weak connections and inconsistent connection topographies as compared to its posterior DMN-nodes (posterior cingulate and posterior parietal cortices). We also present evidence that the marmoset medial prefrontal cortex and temporal lobe areas correspond to other macroscopical distributed functional systems that are not part of the DMN. Given the importance of the marmoset as a pre-clinical primate model for higher-order cognitive functioning and the DMN's relevance to cognition, our results suggest that the marmoset may lack the capacity to integrate neural information to subserve cortical dynamics that is necessary for supporting diverse cognitive demands.

Comparison Between Gradients and Parcellations for Functional Connectivity Prediction of Behavior

Resting-state functional connectivity (RSFC) is widely used to predict behavioral measures. To predict behavioral measures, representing RSFC with parcellations and gradients are the two most popular approaches. Here, we compare parcellation and gradient approaches for RSFC-based prediction of a broad range of behavioral measures in the Human Connectome Project (HCP) and Adolescent Brain Cognitive Development (ABCD) datasets. Among the parcellation approaches, we consider group-average "hard" parcellations (Schaefer et al., 2018), individual-specific "hard" parcellations (Kong et al., 2021a), and an individual-specific "soft" parcellation (spatial independent component analysis with dual regression; Beckmann et al., 2009). For gradient approaches, we consider the well-known principal gradients (Margulies et al., 2016) and the local gradient approach that detects local RSFC changes (Laumann et al., 2015). Across two regression algorithms, individual-specific hard-parcellation performs the best in the HCP dataset, while the principal gradients, spatial independent component analysis and group-average "hard" parcellations exhibit similar performance. On the other hand, principal gradients and all parcellation approaches perform similarly in the ABCD dataset. Across both datasets, local gradients perform the worst. Finally, we find that the principal gradient approach requires at least 40 to 60 gradients to perform as well as parcellation approaches. While most principal gradient studies utilize a single gradient, our results suggest that incorporating higher order gradients can provide significant behaviorally relevant information. Future work will consider the inclusion of additional parcellation and gradient approaches for comparison.

Omnipresence of the sensorimotor-association axis topography in the human connectome

Low-dimensional representations are increasingly used to study meaningful organizational principles within the human brain. Most notably, the sensorimotor-association axis consistently explains the most variance in the human connectome as its so-called principal gradient, suggesting that it represents a fundamental organizational principle. While recent work indicates these low dimensional representations are relatively robust, they are limited by modeling only certain aspects of the functional connectivity structure. To date, the majority of studies have restricted these approaches to the strongest connections in the brain, treating weaker or negative connections as noise despite evidence of meaningful structure among them. The present work examines connectivity gradients of the human connectome across a full range of connectivity strengths and explores the implications for outcomes of individual differences, identifying potential dependencies on thresholds and opportunities to improve prediction tasks. Interestingly, the sensorimotor-association axis emerged as the principal gradient of the human connectome across the entire range of connectivity levels. Moreover, the principal gradient of connections at intermediate strengths encoded individual differences, better followed individual-specific anatomical features, and was also more predictive of intelligence. Taken together, our results add to evidence of the sensorimotor-association axis as a fundamental principle of the brain's functional organization, since it is evident even in the connectivity structure of more lenient connectivity thresholds. These more loosely coupled connections further appear to contain valuable and potentially important information that could be used to improve our understanding of individual differences, diagnosis, and the prediction of treatment outcomes.

Graded functional organisation in the left inferior frontal gyrus: evidence from task-free and task-based functional connectivity

The left inferior frontal gyrus (LIFG) has been ascribed key roles in numerous cognitive domains, including language, executive function and social cognition. However, its functional organisation, and how the specific areas implicated in these cognitive domains relate to each other, is unclear. Possibilities include that the LIFG underpins a domain-general function or, alternatively, that it is characterized by functional differentiation, which might occur in either a discrete or a graded pattern. The aim of the present study was to explore the topographical organisation of the LIFG using a bimodal data-driven approach. To this end, we extracted functional connectivity (FC) gradients from 1) the resting-state fMRI time-series of 150 participants (77 female), and 2) patterns of co-activation derived meta-analytically from task data across a diverse set of cognitive domains. We then sought to characterize the FC differences driving these gradients with seed-based resting-state FC and meta-analytic co-activation modelling analyses. Both analytic approaches converged on an FC profile that shifted in a graded fashion along two main organisational axes. An anterior-posterior gradient shifted from being preferentially associated with high-level control networks (anterior LIFG) to being more tightly coupled with perceptually-driven networks (posterior). A second dorsal-ventral axis was characterized by higher connectivity with domain-general control networks on one hand (dorsal LIFG), and with the semantic network, on the other (ventral). These results provide novel insights into a graded functional organisation of the LIFG underpinning both task-free and task-constrained mental states, and suggest that the LIFG is an interface between distinct large-scale functional networks.

Subcortical Functional Connectivity Gradients in Temporal Lobe Epilepsy

Background and Motivation

Functional gradients have been used to study differences in connectivity between healthy and diseased brain states, however this work has largely focused on the cortex. Because the subcortex plays a key role in seizure initiation in temporal lobe epilepsy (TLE), subcortical functional-connectivity gradients may help further elucidate differences between healthy brains and TLE, as well as differences between left (L)-TLE and right (R)-TLE.

Methods

In this work, we calculated subcortical functional-connectivity gradients (SFGs) from resting-state functional MRI (rs-fMRI) by measuring the similarity in connectivity profiles of subcortical voxels to cortical gray matter voxels. We performed this analysis in 23 R-TLE patients and 32 L-TLE patients (who were otherwise matched for age, gender, disease specific characteristics, and other clinical variables), and 16 controls. To measure differences in SFGs between L-TLE and R-TLE, we quantified deviations in the average functional gradient distributions, as well as their variance, across subcortical structures.

Results

We found an expansion, measured by increased variance, in the principal SFG of TLE relative to controls. When comparing the gradient across subcortical structures between L-TLE and R-TLE, we found that abnormalities in the ipsilateral hippocampal gradient distributions were significantly different between L-TLE and R-TLE.

Conclusion

Our results suggest that expansion of the SFG is characteristic of TLE. Subcortical functional gradient differences exist between left and right TLE and are driven by connectivity changes in the hippocampus ipsilateral to the seizure onset zone.

Does the current state of biomarker discovery in autism reflect the limits of reductionism in precision medicine? Suggestions for an integrative approach that considers dynamic mechanisms between brain, body, and the social environment

Over the past decade, precision medicine has become one of the most influential approaches in biomedical research to improve early detection, diagnosis, and prognosis of clinical conditions and develop mechanism-based therapies tailored to individual characteristics using biomarkers. This perspective article first reviews the origins and concept of precision medicine approaches to autism and summarises recent findings from the first "generation" of biomarker studies. Multi-disciplinary research initiatives created substantially larger, comprehensively characterised cohorts, shifted the focus from group-comparisons to individual variability and subgroups, increased methodological rigour and advanced analytic innovations. However, although several candidate markers with probabilistic value have been identified, separate efforts to divide autism by molecular, brain structural/functional or cognitive markers have not identified a validated diagnostic subgroup. Conversely, studies of specific monogenic subgroups revealed substantial variability in biology and behaviour. The second part discusses both conceptual and methodological factors in these findings. It is argued that the predominant reductionist approach, which seeks to parse complex issues into simpler, more tractable units, let us to neglect the interactions between brain and body, and divorce individuals from their social environment. The third part draws on insights from systems biology, developmental psychology and neurodiversity approaches to outline an integrative approach that considers the dynamic interaction between biological (brain, body) and social mechanisms (stress, stigma) to understanding the origins of autistic features in particular conditions and contexts. This requires 1) closer collaboration with autistic people to increase face validity of concepts and methodologies; (2) development of measures/technologies that enable repeat assessment of social and biological factors in different (naturalistic) conditions and contexts, (3) new analytic methods to study (simulate) these interactions (including emergent properties), and (4) cross-condition designs to understand which mechanisms are transdiagnostic or specific for particular autistic sub-populations. Tailored support may entail both creating more favourable conditions in the social environment and interventions for some autistic people to increase well-being.

The thalamus in psychosis spectrum disorder

Psychosis spectrum disorder (PSD) affects 1% of the world population and results in a lifetime of chronic disability, causing devastating personal and economic consequences. Developing new treatments for PSD remains a challenge, particularly those that target its core cognitive deficits. A key barrier to progress is the tenuous link between the basic neurobiological understanding of PSD and its clinical phenomenology. In this perspective, we focus on a key opportunity that combines innovations in non-invasive human neuroimaging with basic insights into thalamic regulation of functional cortical connectivity. The thalamus is an evolutionary conserved region that forms forebrain-wide functional loops critical for the transmission of external inputs as well as the construction and update of internal models. We discuss our perspective across four lines of evidence: First, we articulate how PSD symptomatology may arise from a faulty network organization at the macroscopic circuit level with the thalamus playing a central coordinating role. Second, we discuss how recent animal work has mechanistically clarified the properties of thalamic circuits relevant to regulating cortical dynamics and cognitive function more generally. Third, we present human neuroimaging evidence in support of thalamic alterations in PSD, and propose that a similar "thalamocortical dysconnectivity" seen in pharmacological imaging (under ketamine, LSD and THC) in healthy individuals may link this circuit phenotype to the common set of symptoms in idiopathic and drug-induced psychosis. Lastly, we synthesize animal and human work, and lay out a translational path for biomarker and therapeutic development.

Striatal connectopic maps link to functional domains across psychiatric disorders

Transdiagnostic approaches to psychiatry have significant potential in overcoming the limitations of conventional diagnostic paradigms. However, while frameworks such as the Research Domain Criteria have garnered significant enthusiasm among researchers and clinicians from a theoretical angle, examples of how such an approach might translate in practice to understand the biological mechanisms underlying complex patterns of behaviors in realistic and heterogeneous populations have been sparse. In a richly phenotyped clinical sample (n = 186) specifically designed to capture the complex nature of heterogeneity and comorbidity within- and between stress- and neurodevelopmental disorders, we use exploratory factor analysis on a wide range of clinical questionnaires to identify four stable functional domains that transcend diagnosis and relate to negative valence, cognition, social functioning and inhibition/arousal before replicating them in an independent dataset (n = 188). We then use connectopic mapping to map inter-individual variation in fine-grained topographical organization of functional connectivity in the striatum-a central hub in motor, cognitive, affective and reward-related brain circuits-and use multivariate machine learning (canonical correlation analysis) to show that these individualized topographic representations predict transdiagnostic functional domains out of sample (r = 0.20, p = 0.026). We propose that investigating psychiatric symptoms across disorders is a promising path to linking them to underlying biology, and can help bridge the gap between neuroscience and clinical psychiatry.

Functional Gradients of the Cerebellum: a Review of Practical Applications

Gradient-based analyses have contributed to the description of cerebellar functional neuroanatomy. More recently, functional gradients of the cerebellum have been used as a multi-purpose tool for neuroimaging research. Here, we provide an overview of the many practical applications of cerebellar functional gradient analyses. These practical applications include examination of intra-cerebellar and cerebellar-extracerebellar organization; transformation of functional gradients into parcellations with discrete borders; projection of functional gradients calculated within cerebellar structures to other extracerebellar structures; interpretation of cerebellar neuroimaging findings using qualitative and quantitative methods; detection of differences in patient populations; and other more complex practical applications of cerebellar gradient-based analyses. This review may serve as an introduction and catalog of options for neuroscientists who wish to design and analyze imaging studies using functional gradients of the cerebellum.

The human egomotion network

All volitional movement in a three-dimensional space requires multisensory integration, in particular of visual and vestibular signals. Where and how the human brain processes and integrates self-motion signals remains enigmatic. Here, we applied visual and vestibular self-motion stimulation using fast and precise whole-brain neuroimaging to delineate and characterize the entire cortical and subcortical egomotion network in a substantial cohort (n=131). Our results identify a core egomotion network consisting of areas in the cingulate sulcus (CSv, PcM/pCi), the cerebellum (uvula), and the temporo-parietal cortex including area VPS and an unnamed region in the supramarginal gyrus. Based on its cerebral connectivity pattern and anatomical localization, we propose that this region represents the human homologue of macaque area 7a. Whole-brain connectivity and gradient analyses imply an essential role of the connections between the cingulate sulcus and the cerebellar uvula in egomotion perception. This could be via feedback loops involved updating visuo-spatial and vestibular information. The unique functional connectivity patterns of PcM/pCi hint at central role in multisensory integration essential for the perception of self-referential spatial awareness. All cortical egomotion hubs showed modular functional connectivity with other visual, vestibular, somatosensory and higher order motor areas, underlining their mutual function in general sensorimotor integration.

A dynamic gradient architecture generates brain activity states

The human brain exhibits a diverse yet constrained range of activity states. While these states can be faithfully represented in a low-dimensional latent space, our understanding of the constitutive functional anatomy is still evolving. Here we applied dimensionality reduction to task-free and task fMRI data to address whether latent dimensions reflect intrinsic systems and if so, how these systems may interact to generate different activity states. We find that each dimension represents a dynamic activity gradient, including a primary unipolar sensory-association gradient underlying the global signal. The gradients appear stable across individuals and cognitive states, while recapitulating key functional connectivity properties including anticorrelation, modularity, and regional hubness. We then use dynamical systems modeling to show that gradients causally interact via state-specific coupling parameters to create distinct brain activity patterns. Together, these findings indicate that a set of dynamic, intrinsic spatial gradients interact to determine the repertoire of possible brain activity states.

Functional Parcellation of Human Brain Using Localized Topo-Connectivity Mapping

The analysis of connectivity between parcellated regions of cortex provides insights into the functional architecture of the brain at a systems level. However, the derivation of functional structures from voxel-wise analyses at finer scales remains a challenge. We propose a novel method, called localized topo-connectivity mapping with singular-value-decomposition-informed filtering (or filtered LTM), to identify and characterize voxel-wise functional structures in the human brain from resting-state fMRI data. Here we describe its mathematical formulation and provide a proof-of-concept using simulated data that allow an intuitive interpretation of the results of filtered LTM. The algorithm has also been applied to 7T fMRI data acquired as part of the Human Connectome Project to generate group-average LTM images. Generally, most of the functional structures revealed by LTM images agree in the boundaries with anatomical structures identified by T1-weighted images and fractional anisotropy maps derived from diffusion MRI. In addition, the LTM images also reveal subtle functional variations that are not apparent in the anatomical structures. To assess the performance of LTM images, the subcortical region and occipital white matter were separately parcellated. Statistical tests were performed to demonstrate that the synchronies of fMRI signals in LTM-derived functional parcels are significantly larger than those with geometric perturbations. Overall, the filtered LTM approach can serve as a tool to investigate the functional organization of the brain at the scale of individual voxels as measured in fMRI.

An evaluation of how connectopic mapping reveals visual field maps in V1

Abstract

Functional gradients, in which response properties change gradually across the cortical surface, have been proposed as a key organising principle of the brain. However, the presence of these gradients remains undetermined in many brain regions. Resting-state neuroimaging studies have suggested these gradients can be reconstructed from patterns of functional connectivity. Here we investigate the accuracy of these reconstructions and establish whether it is connectivity or the functional properties within a region that determine these “connectopic maps”. Different manifold learning techniques were used to recover visual field maps while participants were at rest or engaged in natural viewing. We benchmarked these reconstructions against maps measured by traditional visual field mapping. We report an initial exploratory experiment of a publicly available naturalistic imaging dataset, followed by a preregistered replication using larger resting-state and naturalistic imaging datasets from the Human Connectome Project. Connectopic mapping accurately predicted visual field maps in primary visual cortex, with better predictions for eccentricity than polar angle maps. Non-linear manifold learning methods outperformed simpler linear embeddings. We also found more accurate predictions during natural viewing compared to resting-state. Varying the source of the connectivity estimates had minimal impact on the connectopic maps, suggesting the key factor is the functional topography within a brain region. The application of these standardised methods for connectopic mapping will allow the discovery of functional gradients across the brain.

Protocol registration

The stage 1 protocol for this Registered Report was accepted in principle on 19 April 2022. The protocol, as accepted by the journal, can be found at 10.6084/m9.figshare.19771717.

Multiscale neural gradients reflect transdiagnostic effects of major psychiatric conditions on cortical morphology

It is increasingly recognized that multiple psychiatric conditions are underpinned by shared neural pathways, affecting similar brain systems. Here, we carried out a multiscale neural contextualization of shared alterations of cortical morphology across six major psychiatric conditions (autism spectrum disorder, attention deficit/hyperactivity disorder, major depression disorder, obsessive-compulsive disorder, bipolar disorder, and schizophrenia). Our framework cross-referenced shared morphological anomalies with respect to cortical myeloarchitecture and cytoarchitecture, as well as connectome and neurotransmitter organization. Pooling disease-related effects on MRI-based cortical thickness measures across six ENIGMA working groups, including a total of 28,546 participants (12,876 patients and 15,670 controls), we identified a cortex-wide dimension of morphological changes that described a sensory-fugal pattern, with paralimbic regions showing the most consistent alterations across conditions. The shared disease dimension was closely related to cortical gradients of microstructure as well as neurotransmitter axes, specifically cortex-wide variations in serotonin and dopamine. Multiple sensitivity analyses confirmed robustness with respect to slight variations in analytical choices. Our findings embed shared effects of common psychiatric conditions on brain structure in multiple scales of brain organization, and may provide insights into neural mechanisms of transdiagnostic vulnerability.

Topographic organization of the human caudate functional connectivity and age-related changes with resting-state fMRI

The striatum is postulated to play a central role in gating cortical processing during goal-oriented behavior. While many human neuroimaging studies have treated the striatum as an undivided whole or several homogeneous compartments, some recent studies showed that its circuitry is topographically organized and has more complex relations with the cortical networks than previously assumed. Here, we took a gradient functional connectivity mapping approach that utilizes the entire anatomical space of the caudate nucleus to examine the organization of its functional relationship with the rest of the brain and how its topographic mapping changes with age. We defined the topography of the caudate functional connectivity using three publicly available resting-state fMRI datasets. We replicated and extended previous findings. First, we found two stable gradients of caudate connectivity patterns along its medial-lateral (M-L) and anterior-posterior (A-P) axes, supporting findings from previous tract-tracing studies of non-human primates that there are at least two main organizational principles within the caudate nucleus. Second, unlike previous emphasis of the A-P topology, we showed that the differential connectivity patterns along the M-L gradient of caudate are more clearly organized with the large-scale neural networks; such that brain networks associated with internal vs. external orienting behavior are respectively more closely linked to the medial vs. lateral extent of the caudate. Third, the caudate's M-L organization showed greater age-related reduction in integrity, which was further associated with age-related changes in behavioral measures of executive functions. In sum, our analysis confirmed a sometimes overlooked M-L functional connectivity gradient within the caudate nucleus, with its lateral longitudinal zone more closely linked to the frontoparietal cortical circuits and age-related changes in cognitive control. These findings provide a more precise mapping of the human caudate functional connectivity, both in terms of the gradient organization with cortical networks and age-related changes in such organization.