Connectopic mapping with resting-state fMRI

Functional gradients in the human lateral prefrontal cortex revealed by a comprehensive coordinate-based meta-analysis

The lateral prefrontal cortex (LPFC) of humans enables flexible goal-directed behavior. However, its functional organization remains actively debated after decades of research. Moreover, recent efforts aiming to map the LPFC through meta-analysis are limited, either in scope or in the inferred specificity of structure-function associations. These limitations are in part due to the limited expressiveness of commonly-used data analysis tools, which restricts the breadth and complexity of questions that can be expressed in a meta-analysis. Here, we adopt NeuroLang, a novel approach to more expressive meta-analysis based on probabilistic first-order logic programming, to infer the organizing principles of the LPFC from 14,371 neuroimaging studies. Our findings reveal a rostrocaudal and a dorsoventral gradient, respectively explaining the most and second most variance in meta-analytic connectivity across the LPFC. Moreover, we identify a unimodal-to-transmodal spectrum of coactivation patterns along with a concrete-to-abstract axis of structure-function associations extending from caudal to rostral regions of the LPFC. Finally, we infer inter-hemispheric asymmetries along the principal rostrocaudal gradient, identifying hemisphere-specific associations with topics of language, memory, response inhibition, and sensory processing. Overall, this study provides a comprehensive meta-analytic mapping of the LPFC, grounding future hypothesis generation on a quantitative overview of past findings.

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.

An Open MRI Dataset For Multiscale Neuroscience

Multimodal neuroimaging grants a powerful window into the structure and function of the human brain at multiple scales. Recent methodological and conceptual advances have enabled investigations of the interplay between large-scale spatial trends (also referred to as gradients) in brain microstructure and connectivity, offering an integrative framework to study multiscale brain organization. Here, we share a multimodal MRI dataset for Microstructure-Informed Connectomics (MICA-MICs) acquired in 50 healthy adults (23 women; 29.54 ± 5.62 years) who underwent high-resolution T1-weighted MRI, myelin-sensitive quantitative T1 relaxometry, diffusion-weighted MRI, and resting-state functional MRI at 3 Tesla. In addition to raw anonymized MRI data, this release includes brain-wide connectomes derived from (i) resting-state functional imaging, (ii) diffusion tractography, (iii) microstructure covariance analysis, and (iv) geodesic cortical distance, gathered across multiple parcellation scales. Alongside, we share large-scale gradients estimated from each modality and parcellation scale. Our dataset will facilitate future research examining the coupling between brain microstructure, connectivity, and function. MICA-MICs is available on the Canadian Open Neuroscience Platform data portal ( https://portal.conp.ca ) and the Open Science Framework ( https://osf.io/j532r/ ).

In-vivo data-driven parcellation of Heschl's gyrus using structural connectivity

The human auditory cortex around Heschl’s gyrus (HG) exhibits diverging patterns across individuals owing to the heterogeneity of its substructures. In this study, we investigated the subregions of the human auditory cortex using data-driven machine-learning techniques at the individual level and assessed their structural and functional profiles. We studied an openly accessible large dataset of the Human Connectome Project and identified the subregions of the HG in humans using data-driven clustering techniques with individually calculated imaging features of cortical folding and structural connectivity information obtained via diffusion magnetic resonance imaging tractography. We characterized the structural and functional profiles of each HG subregion according to the cortical morphology, microstructure, and functional connectivity at rest. We found three subregions. The first subregion (HG₁) occupied the central portion of HG, the second subregion (HG₂) occupied the medial-posterior-superior part of HG, and the third subregion (HG₃) occupied the lateral-anterior-inferior part of HG. The HG₃ exhibited strong structural and functional connectivity to the association and paralimbic areas, and the HG₁ exhibited a higher myelin density and larger cortical thickness than other subregions. A functional gradient analysis revealed a gradual axis expanding from the HG₂ to the HG₃. Our findings clarify the individually varying structural and functional organization of human HG subregions and provide insights into the substructures of the human auditory cortex.

Disorganization of language and working memory systems in frontal versus temporal lobe epilepsy

Cognitive impairment is a common comorbidity of epilepsy and adversely impacts people with both frontal lobe (FLE) and temporal lobe (TLE) epilepsy. While its neural substrates have been investigated extensively in TLE, functional imaging studies in FLE are scarce. In this study, we profiled the neural processes underlying cognitive impairment in FLE and directly compared FLE and TLE to establish commonalities and differences. We investigated 172 adult participants (56 with FLE, 64 with TLE and 52 controls) using neuropsychological tests and four functional MRI tasks probing expressive language (verbal fluency, verb generation) and working memory (verbal and visuo-spatial). Patient groups were comparable in disease duration and anti-seizure medication load. We devised a multiscale approach to map brain activation and deactivation during cognition and track reorganization in FLE and TLE. Voxel-based analyses were complemented with profiling of task effects across established motifs of functional brain organization: (i) canonical resting-state functional systems; and (ii) the principal functional connectivity gradient, which encodes a continuous transition of regional connectivity profiles, anchoring lower-level sensory and transmodal brain areas at the opposite ends of a spectrum. We show that cognitive impairment in FLE is associated with reduced activation across attentional and executive systems, as well as reduced deactivation of the default mode system, indicative of a large-scale disorganization of task-related recruitment. The imaging signatures of dysfunction in FLE are broadly similar to those in TLE, but some patterns are syndrome-specific: altered default-mode deactivation is more prominent in FLE, while impaired recruitment of posterior language areas during a task with semantic demands is more marked in TLE. Functional abnormalities in FLE and TLE appear overall modulated by disease load. On balance, our study elucidates neural processes underlying language and working memory impairment in FLE, identifies shared and syndrome-specific alterations in the two most common focal epilepsies and sheds light on system behaviour that may be amenable to future remediation strategies.

Serotonergic psychedelic drugs LSD and psilocybin reduce the hierarchical differentiation of unimodal and transmodal cortex

Lysergic acid diethylamide (LSD) and psilocybin are serotonergic psychedelic compounds with potential in the treatment of mental health disorders. Past neuroimaging investigations have revealed that both compounds can elicit significant changes to whole-brain functional organization and dynamics. A recent proposal linked past findings into a unified model and hypothesized reduced whole-brain hierarchical organization as a key mechanism underlying the psychedelic state, but this has yet to be directly tested. We applied a non-linear dimensionality reduction technique previously used to map hierarchical connectivity gradients to assess cortical organization in the LSD and psilocybin state from two previously published pharmacological resting-state fMRI datasets (N = 15 and 9, respectively). Results supported our primary hypothesis: The principal gradient of cortical connectivity, describing a hierarchy from unimodal to transmodal cortex, was significantly flattened under both drugs relative to their respective placebo conditions. Between-condition contrasts revealed that this was driven by a reduction of functional differentiation at both hierarchical extremes - default and frontoparietal networks at the upper end, and somatomotor at the lower. Gradient-based connectivity mapping indicated that this was underpinned by a disruption of modular unimodal connectivity and increased unimodal-transmodal crosstalk. Results involving the second and third gradient, which, respectively represent axes of sensory and executive differentiation, also showed significant alterations across both drugs. These findings provide support for a recent mechanistic model of the psychedelic state relevant to therapeutic applications of psychedelics. More fundamentally, we provide the first evidence that macroscale connectivity gradients are sensitive to an acute pharmacological manipulation, supporting a role for psychedelics as scientific tools to perturb cortical functional organization.

Multimodal Gradient Mapping of Rodent Hippocampus

The hippocampus plays a central role in supporting our coherent and enduring sense of self and our place in the world. Understanding its functional organisation is central to understanding this complex role. Previous studies suggest function varies along a long hippocampal axis, but there is disagreement about the presence of sharp discontinuities or gradual change along that axis. Other open questions relate to the underlying drivers of this variation and the conservation of organisational principles across species. Here, we delineate the primary organisational principles underlying patterns of hippocampal functional connectivity (FC) in the mouse using gradient analysis on resting state fMRI data. We further applied gradient analysis to mouse gene co-expression data to examine the relationship between variation in genomic anatomy and functional organisation. Two principal FC gradients along a hippocampal axis were revealed. The principal gradient exhibited a sharp discontinuity that divided the hippocampus into dorsal and ventral compartments. The second, more continuous, gradient followed the long axis of the ventral compartment. Dorsal regions were more strongly connected to areas involved in spatial navigation while ventral regions were more strongly connected to areas involved in emotion, recapitulating patterns seen in humans. In contrast, gene co-expression gradients showed a more segregated and discrete organisation. Our findings suggest that hippocampal functional organisation exhibits both sharp and gradual transitions and that hippocampal genomic anatomy exerts only a subtle influence on this organisation.

Intrinsic timescales as an organizational principle of neural processing across the whole rhesus macaque brain

Hierarchical temporal dynamics are a fundamental computational property of the brain; however, there are no whole brain, noninvasive investigations into timescales of neural processing in animal models. To that end, we used the spatial resolution and sensitivity of ultrahigh field functional magnetic resonance imaging (fMRI) performed at 10.5 T to probe timescales across the whole macaque brain. We uncovered within-species consistency between timescales estimated from fMRI and electrophysiology. Crucially, we extended existing electrophysiological hierarchies to whole-brain topographies. Our results validate the complementary use of hemodynamic and electrophysiological intrinsic timescales, establishing a basis for future translational work. Further, with these results in hand, we were able to show that one facet of the high-dimensional functional connectivity (FC) topography of any region in the brain is closely related to hierarchical temporal dynamics. We demonstrated that intrinsic timescales are organized along spatial gradients that closely match FC gradient topographies across the whole brain. We conclude that intrinsic timescales are a unifying organizational principle of neural processing across the whole brain.

Mapping dopaminergic projections in the human brain with resting-state fMRI

The striatum receives dense dopaminergic projections, making it a key region of the dopaminergic system. Its dysfunction has been implicated in various conditions including Parkinson’s disease (PD) and substance use disorder. However, the investigation of dopamine-specific functioning in humans is problematic as current MRI approaches are unable to differentiate between dopaminergic and other projections. Here, we demonstrate that ‘connectopic mapping’ – a novel approach for characterizing fine-grained, overlapping modes of functional connectivity – can be used to map dopaminergic projections in striatum. We applied connectopic mapping to resting-state functional MRI data of the Human Connectome Project (population cohort; N = 839) and selected the second-order striatal connectivity mode for further analyses. We first validated its specificity to dopaminergic projections by demonstrating a high spatial correlation (r = 0.884) with dopamine transporter availability – a marker of dopaminergic projections – derived from DaT SPECT scans of 209 healthy controls. Next, we obtained the subject-specific second-order modes from 20 controls and 39 PD patients scanned under placebo and under dopamine replacement therapy (L-DOPA), and show that our proposed dopaminergic marker tracks PD diagnosis, symptom severity, and sensitivity to L-DOPA. Finally, across 30 daily alcohol users and 38 daily smokers, we establish strong associations with self-reported alcohol and nicotine use. Our findings provide evidence that the second-order mode of functional connectivity in striatum maps onto dopaminergic projections, tracks inter-individual differences in PD symptom severity and L-DOPA sensitivity, and exhibits strong associations with levels of nicotine and alcohol use, thereby offering a new biomarker for dopamine-related (dys)function in the human brain.

Connectivity gradients on tractography data: Pipeline and example applications

Gray matter connectivity can be described in terms of its topographical organization, but the differential role of white matter connections underlying that organization is often unknown. In this study, we propose a method for unveiling principles of organization of both gray and white matter based on white matter connectivity as assessed using diffusion magnetic ressonance imaging (MRI) tractography with spectral embedding gradient mapping. A key feature of the proposed approach is its capacity to project the individual connectivity gradients it reveals back onto its input data in the form of projection images, allowing one to assess the contributions of specific white matter tracts to the observed gradients. We demonstrate the ability of our proposed pipeline to identify connectivity gradients in prefrontal and occipital gray matter. Finally, leveraging the use of tractography, we demonstrate that it is possible to observe gradients within the white matter bundles themselves. Together, the proposed framework presents a generalized way to assess both the topographical organization of structural brain connectivity and the anatomical features driving it.

The relation of integrated psychological therapy to resting state functional brain connectivity networks in patients with schizophrenia

Functional brain dysconnectivity measured with resting state functional magnetic resonance imaging (rsfMRI) has been linked to cognitive impairment in schizophrenia. This study investigated the effects on functional brain connectivity of Integrated Psychological Therapy (IPT), a cognitive behavioral oriented group intervention program, in 31 patients with schizophrenia. Patients received IPT or an equal intensity non-specific psychological treatment in a non-randomized design. Evidence of improvement in executive and social functions, psychopathology and overall level of functioning was observed after treatment completion at six months only in the IPT treatment group and was partially sustained at one-year follow up. Independent Component Analysis and Isometric Mapping (ISOMAP), a non-linear manifold learning algorithm, were used to construct functional connectivity networks from the rsfMRI data. Functional brain dysconnectivity was observed in patients compared to a group of 17 healthy controls, both globally and specifically including the default mode (DMN) and frontoparietal network (FPN). DMN and FPN connectivity were reversed towards healthy control patterns only in the IPT treatment group and these effects were sustained at follow up for DMN but not FPN. These data suggest the use of rsfMRI as a biomarker for accessing and monitoring the therapeutic effects of cognitive remediation therapy in schizophrenia.

Neuroimaging brain growth charts: A road to mental health

Mental disorders are common health concerns and contribute to a heavy global burden on our modern society. It is challenging to identify and treat them timely. Neuroimaging evidence suggests the incidence of various psychiatric and behavioral disorders is closely related to the atypical development of brain structure and function. The identification and understanding of atypical brain development provide chances for clinicians to detect mental disorders earlier, perhaps even prior to onset, and treat them more precisely. An invaluable and necessary method in identifying and monitoring atypical brain development are growth charts of typically developing individuals in the population. The brain growth charts can offer a series of standard references on typical neurodevelopment, representing an important resource for the scientific and medical communities. In the present paper, we review the relationship between mental disorders and atypical brain development from a perspective of normative brain development by surveying the recent progress in the development of brain growth charts, including four aspects on growth chart utility: 1) cohorts, 2) measures, 3) mechanisms, and 4) clinical translations. In doing so, we seek to clarify the challenges and opportunities in charting brain growth, and to promote the application of brain growth charts in clinical practice.

Resting-state functional heterogeneity of the right insula contributes to pain sensitivity

Previous studies have described the structure and function of the insular cortex in terms of spatially continuous gradients. Here we assess how spatial features of insular resting state functional organization correspond to individual pain sensitivity. From a previous multicenter study, we included 107 healthy participants, who underwent resting state functional MRI scans, T1-weighted scans and quantitative sensory testing on the left forearm. Thermal and mechanical pain thresholds were determined. Connectopic mapping, a technique using non-linear representations of functional organization was employed to describe functional connectivity gradients in both insulae. Partial coefficients of determination were calculated between trend surface model parameters summarizing spatial features of gradients, modal and modality-independent pain sensitivity. The dominant connectopy captured the previously reported posteroanterior shift in connectivity profiles. Spatial features of dominant connectopies in the right insula explained significant amounts of variance in thermal (R² = 0.076; p < 0.001 and R² = 0.031; p < 0.029) and composite pain sensitivity (R² = 0.072; p < 0.002). The left insular gradient was not significantly associated with pain thresholds. Our results highlight the functional relevance of gradient-like insular organization in pain processing. Considering individual variations in insular connectopy might contribute to understanding neural mechanisms behind pain and improve objective brain-based characterization of individual pain sensitivity.

Recent developments in representations of the connectome

Research into the human connectome (i.e., all connections in the human brain) with the use of resting state functional MRI has rapidly increased in popularity in recent years, especially with the growing availability of large-scale neuroimaging datasets. The goal of this review article is to describe innovations in functional connectome representations that have come about in the past 8 years, since the 2013 NeuroImage special issue on 'Mapping the Connectome'. In the period, research has shifted from group-level brain parcellations towards the characterization of the individualized connectome and of relationships between individual connectomic differences and behavioral/clinical variation. Achieving subject-specific accuracy in parcel boundaries while retaining cross-subject correspondence is challenging, and a variety of different approaches are being developed to meet this challenge, including improved alignment, improved noise reduction, and robust group-to-subject mapping approaches. Beyond the interest in the individualized connectome, new representations of the data are being studied to complement the traditional parcellated connectome representation (i.e., pairwise connections between distinct brain regions), such as methods that capture overlapping and smoothly varying patterns of connectivity ('gradients'). These different connectome representations offer complimentary insights into the inherent functional organization of the brain, but challenges for functional connectome research remain. Interpretability will be improved by future research towards gaining insights into the neural mechanisms underlying connectome observations obtained from functional MRI. Validation studies comparing different connectome representations are also needed to build consensus and confidence to proceed with clinical trials that may produce meaningful clinical translation of connectome insights.

Striatal topographical organization: Bridging the gap between molecules, connectivity and behavior

The striatum represents the major hub of the basal ganglia, receiving projections from the entire cerebral cortex and it is assumed to play a key role in a wide array of complex behavioral tasks. Despite being extensively investigated during the last decades, the topographical organization of the striatum is not well understood yet. Ongoing efforts in neuroscience are focused on analyzing striatal anatomy at different spatial scales, to understand how structure relates to function and how derangements of this organization are involved in various neuropsychiatric diseases. While being subdivided at the macroscale level into dorsal and ventral divisions, at a mesoscale level the striatum represents an anatomical continuum sharing the same cellular makeup. At the same time, it is now increasingly ascertained that different striatal compartments show subtle histochemical differences, and their neurons exhibit peculiar patterns of gene expression, supporting functional diversity across the whole basal ganglia circuitry. Such diversity is further supported by afferent connections which are heterogenous both anatomically, as they originate from distributed cortical areas and subcortical structures, and biochemically, as they involve a variety of neurotransmitters. Specifically, the cortico-striatal projection system is topographically organized delineating a functional organization which is maintained throughout the basal ganglia, subserving motor, cognitive and affective behavioral functions. While such functional heterogeneity has been firstly conceptualized as a tripartite organization, with sharply defined limbic, associative and sensorimotor territories within the striatum, it has been proposed that such territories are more likely to fade into one another, delineating a gradient-like organization along medio-lateral and ventro-dorsal axes. However, the molecular and cellular underpinnings of such organization are less understood, and their relations to behavior remains an open question, especially in humans. In this review we aimed at summarizing the available knowledge on striatal organization, especially focusing on how it links structure to function and its alterations in neuropsychiatric diseases. We examined studies conducted on different species, covering a wide array of different methodologies: from tract-tracing and immunohistochemistry to neuroimaging and transcriptomic experiments, aimed at bridging the gap between macroscopic and molecular levels.

Striatal Subdivisions Estimated via Deep Embedded Clustering With Application to Parkinson's Disease

Recent fMRI connectivity-based parcellation (CBP) methods have been developed to obtain homogeneous and functionally coherent brain parcels. However, most of these studies utilize traditional clustering methods that neglect hidden nonlinear features. To enhance parcellation performance, here we propose a deep embedded connectivity-based parcellation (DECBP) framework and apply it to determine functional subdivisions of the striatum in public resting state fMRI data sets. This framework integrates fMRI connectivity features into deep embedded clustering (DEC), a deep neural network based on a stacked autoencoder. Compared to three prevalent clustering methods and their combinations with principal component analysis (PCA), the DECBP exhibited a significantly higher similarity between scans, individuals, and groups, indicating enhanced reproducibility. The generated reliable parcellations were also largely consistent with other public atlases. We further explored the functional subunits in the striatum in a data set from 23 Parkinson's disease (PD) subjects and 27 age-matched healthy controls (HC). All putaminal subregions of PD demonstrated lower interhemispheric connectivity than those of HC, which might reflect imbalance in the pathological progression of PD. Such hypo-connectivity was also observed between putaminal subregions and other brain regions, reflecting neuroimaging manifestations of the altered cortico-striato-thalamo-cortical circuit. These observed weaker couplings were associated with PD severity and duration. Our results support the utilization of the DECBP framework and suggest that abnormal connectivity in putaminal subregions may be a potential indicator of PD.

Hierarchical modelling of functional brain networks in population and individuals from big fMRI data

•

We introduce stochastic PROFUMO (sPROFUMO) for inferring functional brain networks from big data.

•

sPROFUMO hierarchically estimates fMRI networks for the population and every individual.

•

We characterised high dimensional resting state fMRI networks from UK Biobank.

•

Model outperforms ICA and dual regression for estimation of individual-specific network topography.

•

We demonstrate the model's utility for predicting cognitive traits, and capturing subject variability in network topographies versus connectivity.

A major goal of large-scale brain imaging datasets is to provide resources for investigating heterogeneous populations. Characterisation of functional brain networks for individual subjects from these datasets will have an enormous potential for prediction of cognitive or clinical traits. We propose for the first time a technique, Stochastic Probabilistic Functional Modes (sPROFUMO), that is scalable to UK Biobank (UKB) with expected 100,000 participants, and hierarchically estimates functional brain networks in individuals and the population, while allowing for bidirectional flow of information between the two. Using simulations, we show the model's utility, especially in scenarios that involve significant cross-subject variability, or require delineation of fine-grained differences between the networks. Subsequently, by applying the model to resting-state fMRI from 4999 UKB subjects, we mapped resting state networks (RSNs) in single subjects with greater detail than has been possible previously in UKB (>100 RSNs), and demonstrate that these RSNs can predict a range of sensorimotor and higher-level cognitive functions. Furthermore, we demonstrate several advantages of the model over independent component analysis combined with dual-regression (ICA-DR), particularly with respect to the estimation of the spatial configuration of the RSNs and the predictive power for cognitive traits. The proposed model and results can open a new door for future investigations into individualised profiles of brain function from big data.

Functional harmonics reveal multi-dimensional basis functions underlying cortical organization

The human brain consists of specialized areas that flexibly interact to form a multitude of functional networks. Complementary to this notion of modular organization, brain function has been shown to vary along a smooth continuum across the whole cortex. We demonstrate a mathematical framework that accounts for both of these perspectives: harmonic modes. We calculate the harmonic modes of the brain's functional connectivity graph, called "functional harmonics," revealing a multi-dimensional, frequency-ordered set of basis functions. Functional harmonics link characteristics of cortical organization across several spatial scales, capturing aspects of intra-areal organizational features (retinotopy, somatotopy), delineating brain areas, and explaining macroscopic functional networks as well as global cortical gradients. Furthermore, we show how the activity patterns elicited by seven different tasks are reconstructed from a very small subset of functional harmonics. Our results suggest that the principle of harmonicity, ubiquitous in nature, also underlies functional cortical organization in the human brain.

Atypical Integration of Sensory-to-Transmodal Functional Systems Mediates Symptom Severity in Autism

A notable characteristic of autism spectrum disorder (ASD) is co-occurring deficits in low-level sensory processing and high-order social interaction. While there is evidence indicating detrimental cascading effects of sensory anomalies on the high-order cognitive functions in ASD, the exact pathological mechanism underlying their atypical functional interaction across the cortical hierarchy has not been systematically investigated. To address this gap, here we assessed the functional organisation of sensory and motor areas in ASD, and their relationship with subcortical and high-order trandmodal systems. In a resting-state fMRI data of 107 ASD and 113 neurotypical individuals, we applied advanced connectopic mapping to probe functional organization of primary sensory/motor areas, together with targeted seed-based intrinsic functional connectivity (iFC) analyses. In ASD, the connectopic mapping revealed topological anomalies (i.e., excessively more segregated iFC) in the motor and visual areas, the former of which patterns showed association with the symptom severity of restricted and repetitive behaviors. Moreover, the seed-based analysis found diverging patterns of ASD-related connectopathies: decreased iFCs within the sensory/motor areas but increased iFCs between sensory and subcortical structures. While decreased iFCs were also found within the higher-order functional systems, the overall proportion of this anomaly tends to increase along the level of cortical hierarchy, suggesting more dysconnectivity in the higher-order functional networks. Finally, we demonstrated that the association between low-level sensory/motor iFCs and clinical symptoms in ASD was mediated by the high-order transmodal systems, suggesting pathogenic functional interactions along the cortical hierarchy. Findings were largely replicated in the independent dataset. These results highlight that atypical integration of sensory-to-high-order systems contributes to the complex ASD symptomatology.

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

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

Shifting gradients of macroscale cortical organization mark the transition from childhood to adolescence

The transition from childhood to adolescence is marked by pronounced shifts in brain structure and function that coincide with the development of physical, cognitive, and social abilities. Prior work in adult populations has characterized the topographical organization of the cortex, revealing macroscale functional gradients that extend from unimodal (somatosensory/motor and visual) regions through the cortical association areas that underpin complex cognition in humans. However, the presence of these core functional gradients across development as well as their maturational course have yet to be established. Here, leveraging 378 resting-state functional MRI scans from 190 healthy individuals aged 6 to 17 y old, we demonstrate that the transition from childhood to adolescence is reflected in the gradual maturation of gradient patterns across the cortical sheet. In children, the overarching organizational gradient is anchored within the unimodal cortex, between somatosensory/motor and visual territories. Conversely, in adolescence, the principal gradient of connectivity transitions into an adult-like spatial framework, with the default network at the opposite end of a spectrum from primary sensory and motor regions. The observed gradient transitions are gradually refined with age, reaching a sharp inflection point in 13 and 14 y olds. Functional maturation was nonuniformly distributed across cortical networks. Unimodal networks reached their mature positions early in development, while association regions, in particular the medial prefrontal cortex, reached a later peak during adolescence. These data reveal age-dependent changes in the macroscale organization of the cortex and suggest the scheduled maturation of functional gradient patterns may be critically important for understanding how cognitive and behavioral capabilities are refined across development.

Motor Dysfunction of Autonomic Nervous System Based on Resting State Functional Magnetic Resonance Imaging

Epilepsy is caused by highly synchronized abnormal discharges of neurons in the brain. Resting State FMRI is widely used as a non-invasive checking method, clinical and basic research in the field of epilepsy has a huge influence, to study the plant nerve system movement function of brain activity and the neural network function of the complex contact provides broad prospects, greatly promote the development of clinical neurology and imaging. Therefore, this paper studies the functional connection and cognitive function of verbal working memory (VWM) in patients with epilepsy by applying resting state FMRI. The experimental results showed that VWM and cognitive functions of epileptic patients showed a trend of decline or even disappearance, to realize the detection of autonomic nervous system motor dysfunction of patients.

Structural Connectivity Gradients of the Temporal Lobe Serve as Multiscale Axes of Brain Organization and Cortical Evolution

The temporal lobe is implicated in higher cognitive processes and is one of the regions that underwent substantial reorganization during primate evolution. Its functions are instantiated, in part, by the complex layout of its structural connections. Here, we identified low-dimensional representations of structural connectivity variations in human temporal cortex and explored their microstructural underpinnings and associations to macroscale function. We identified three eigenmodes which described gradients in structural connectivity. These gradients reflected inter-regional variations in cortical microstructure derived from quantitative magnetic resonance imaging and postmortem histology. Gradient-informed models accurately predicted macroscale measures of temporal lobe function. Furthermore, the identified gradients aligned closely with established measures of functional reconfiguration and areal expansion between macaques and humans, highlighting their potential role in shaping temporal lobe function throughout primate evolution. Findings were replicated in several datasets. Our results provide robust evidence for three axes of structural connectivity in human temporal cortex with consistent microstructural underpinnings and contributions to large-scale brain network function.

Inter-individual body mass variations relate to fractionated functional brain hierarchies

Variations in body mass index (BMI) have been suggested to relate to atypical brain organization, yet connectome-level substrates of BMI and their neurobiological underpinnings remain unclear. Studying 325 healthy young adults, we examined associations between functional connectivity and inter-individual BMI variations. We utilized non-linear connectome manifold learning techniques to represent macroscale functional organization along continuous hierarchical axes that dissociate low level and higher order brain systems. We observed an increased differentiation between unimodal and heteromodal association networks in individuals with higher BMI, indicative of a disrupted modular architecture and hierarchy of the brain. Transcriptomic decoding and gene enrichment analyses identified genes previously implicated in genome-wide associations to BMI and specific cortical, striatal, and cerebellar cell types. These findings illustrate functional connectome substrates of BMI variations in healthy young adults and point to potential molecular associations.

Gradients of connectivity as graph Fourier bases of brain activity

The application of graph theory to model the complex structure and function of the brain has shed new light on its organization, prompting the emergence of network neuroscience. Despite the tremendous progress that has been achieved in this field, still relatively few methods exploit the topology of brain networks to analyze brain activity. Recent attempts in this direction have leveraged on the one hand graph spectral analysis (to decompose brain connectivity into eigenmodes or gradients) and the other graph signal processing (to decompose brain activity "coupled to" an underlying network in graph Fourier modes). These studies have used a variety of imaging techniques (e.g., fMRI, electroencephalography, diffusion-weighted and myelin-sensitive imaging) and connectivity estimators to model brain networks. Results are promising in terms of interpretability and functional relevance, but methodologies and terminology are variable. The goals of this paper are twofold. First, we summarize recent contributions related to connectivity gradients and graph signal processing, and attempt a clarification of the terminology and methods used in the field, while pointing out current methodological limitations. Second, we discuss the perspective that the functional relevance of connectivity gradients could be fruitfully exploited by considering them as graph Fourier bases of brain activity.

Differences in subcortico-cortical interactions identified from connectome and microcircuit models in autism

The pathophysiology of autism has been suggested to involve a combination of both macroscale connectome miswiring and microcircuit anomalies. Here, we combine connectome-wide manifold learning with biophysical simulation models to understand associations between global network perturbations and microcircuit dysfunctions in autism. We studied neuroimaging and phenotypic data in 47 individuals with autism and 37 typically developing controls obtained from the Autism Brain Imaging Data Exchange initiative. Our analysis establishes significant differences in structural connectome organization in individuals with autism relative to controls, with strong between-group effects in low-level somatosensory regions and moderate effects in high-level association cortices. Computational models reveal that the degree of macroscale anomalies is related to atypical increases of recurrent excitation/inhibition, as well as subcortical inputs into cortical microcircuits, especially in sensory and motor areas. Transcriptomic association analysis based on postmortem datasets identifies genes expressed in cortical and thalamic areas from childhood to young adulthood. Finally, supervised machine learning finds that the macroscale perturbations are associated with symptom severity scores on the Autism Diagnostic Observation Schedule. Together, our analyses suggest that atypical subcortico-cortical interactions are associated with both microcircuit and macroscale connectome differences in autism.

Mesoscale hierarchical organization of primary somatosensory cortex captured by resting-state-fMRI in humans

The primary somatosensory cortex (S1) plays a key role in the processing and integration of afferent somatosensory inputs along an anterior-to-posterior axis, contributing towards necessary human function. It is believed that anatomical connectivity can be used to probe hierarchical organization, however direct characterization of this principle in-vivo within humans remains elusive. Here, we use resting-state functional connectivity as a complement to anatomical connectivity to investigate topographical principles of human S1. We employ a novel approach to examine mesoscopic variations of functional connectivity, and demonstrate a topographic organisation spanning the region's hierarchical axis that strongly correlates with underlying microstructure while tracing along architectonic Brodmann areas. Our findings characterize anatomical hierarchy of S1 as a 'continuous spectrum' with evidence supporting a functional boundary between areas 3b and 1. The identification of this topography bridges the gap between structure and connectivity, and may be used to help further current understanding of sensorimotor deficits.

Functional connectivity of the hippocampus and its subfields in resting-state networks

Many neuroimaging studies have shown that the hippocampus participates in a resting‐state network called the default mode network. However, how the hippocampus connects to the default mode network, whether the hippocampus connects to other resting‐state networks and how the different hippocampal subfields take part in resting‐state networks remains poorly understood. Here, we examined these issues using the high spatial‐resolution 7T resting‐state fMRI dataset from the Human Connectome Project. We used data‐driven techniques that relied on spatially‐restricted Independent Component Analysis, Dual Regression and linear mixed‐effect group‐analyses based on participant‐specific brain morphology. The results revealed two main activity hotspots inside the hippocampus. The first hotspot was located in an anterior location and was correlated with the somatomotor network. This network was subserved by co‐activity in the CA1, CA3, CA4 and Dentate Gyrus fields. In addition, there was an activity hotspot that extended from middle to posterior locations along the hippocampal long‐axis and correlated with the default mode network. This network reflected activity in the Subiculum, CA4 and Dentate Gyrus fields. These results show how different sections of the hippocampus participate in two known resting‐state networks and how these two resting‐state networks depend on different configurations of hippocampal subfield co‐activity.

An expanding manifold in transmodal regions characterizes adolescent reconfiguration of structural connectome organization

Adolescence is a critical time for the continued maturation of brain networks. Here, we assessed structural connectome development in a large longitudinal sample ranging from childhood to young adulthood. By projecting high-dimensional connectomes into compact manifold spaces, we identified a marked expansion of structural connectomes, with strongest effects in transmodal regions during adolescence. Findings reflected increased within-module connectivity together with increased segregation, indicating increasing differentiation of higher-order association networks from the rest of the brain. Projection of subcortico-cortical connectivity patterns into these manifolds showed parallel alterations in pathways centered on the caudate and thalamus. Connectome findings were contextualized via spatial transcriptome association analysis, highlighting genes enriched in cortex, thalamus, and striatum. Statistical learning of cortical and subcortical manifold features at baseline and their maturational change predicted measures of intelligence at follow-up. Our findings demonstrate that connectome manifold learning can bridge the conceptual and empirical gaps between macroscale network reconfigurations, microscale processes, and cognitive outcomes in adolescent development.

The Graded Change in Connectivity across the Ventromedial Prefrontal Cortex Reveals Distinct Subregions

The functional heterogeneity of the ventromedial prefrontal cortex (vmPFC) suggests it may include distinct functional subregions. To date these have not been well elucidated. Regions with differentiable connectivity (and as a result likely dissociable functions) may be identified using emergent data-driven approaches. However, prior parcellations of the vmPFC have only considered hard splits between distinct regions, although both hard and graded connectivity changes may exist. Here we determine the full pattern of change in structural and functional connectivity across the vmPFC for the first time and extract core distinct regions. Both structural and functional connectivity varied along a dorsomedial to ventrolateral axis from relatively dorsal medial wall regions to relatively lateral basal orbitofrontal cortex. The pattern of connectivity shifted from default mode network to sensorimotor and multimodal semantic connections. This finding extends the classical distinction between primate medial and orbital regions by demonstrating a similar gradient in humans for the first time. Additionally, core distinct regions in the medial wall and orbitofrontal cortex were identified that may show greater correspondence to functional differences than prior hard parcellations. The possible functional roles of the orbitofrontal cortex and medial wall are discussed.

Atypical neural topographies underpin dysfunctional pattern separation in temporal lobe epilepsy

Episodic memory is the ability to accurately remember events from our past. The process of pattern separation is hypothesized to underpin this ability and is defined as the ability to orthogonalize memory traces, to maximize the features that make them unique. Contemporary cognitive neuroscience suggests that pattern separation entails complex interactions between the hippocampus and the neocortex, where specific hippocampal subregions shape neural reinstatement in the neocortex. To test this hypothesis, the current work studied both healthy controls and patients with temporal lobe epilepsy (TLE) who present with hippocampal structural anomalies. In all participants, we measured neural activity using functional magnetic resonance imaging (fMRI) while they retrieved memorized items compared to lure items which share features with the target. Behaviorally, TLE patients were less able to exclude lures than controls, and showed a reduction in pattern separation. To assess the hypothesized relationship between neural patterns in the hippocampus and the neocortex, we identified topographic gradients of intrinsic connectivity along neocortical and hippocampal subfield surfaces and identified the topographic profile of the neural activity accompanying pattern separation. In healthy controls, pattern separation followed a graded pattern of neural activity, both along the hippocampal long axis (and peaked in anterior segments that are more heavily engaged in transmodal processing) and along the neocortical hierarchy running from unimodal to transmodal regions (peaking in transmodal default mode regions). In TLE patients, however, this concordance between task-based functional activations and topographic gradients was markedly reduced. Furthermore, person specific measures of concordance between task-related activity and connectivity gradients in patients and controls related to inter-individual differences in behavioral measures of pattern separation and episodic memory, highlighting the functional relevance of the observed topographic motifs. Our work is consistent with an emerging understanding that successful discrimination between memories with similar features entails a shift in the locus of neural activity away from sensory systems, a pattern that is mirrored along the hippocampal long axis and with respect to neocortical hierarchies. More broadly, our study establishes topographic profiling using intrinsic connectivity gradients captures the functional underpinnings of episodic memory processes in manner that is sensitive to their reorganization in pathology.

SLC40A1 Mediates Ferroptosis and Cognitive Dysfunction in Type 1 Diabetes

Cognitive dysfunction often accompanies diabetes. Both hypoglycemia and hyperglycemia cause cognitive dysfunctions. However, the underlying pathophysiology remains unclear. Recent evidence show that ferroptosis primarily triggers nerve cell death, Alzheimer's disease (AD), Huntington's disease (HD), and Parkinson's disease (PD). The present study aimed to investigate whether ferroptosis is a vital pathogenic pathway in diabetes-induced cognitive dysfunction. Type 1 diabetic rat model was created by intraperitoneal injection of streptozotocin (STZ). Significant cognitive dysfunction was observed in the diabetic rats as evidenced by increase in latency period to find a hidden platform and decreased cumulative time spent in the target quadrant (TQ) in the Morris water maze test. We detected the amplitude of low-frequency fluctuation (ALFF) of the BOLD (Blood Oxygenation Level-Dependent) signal using resting-state functional magnetic resonance imaging (rs-fMRI). Consequently, we found that the ALFF values, as well as the T2 relaxation time of the bilateral hippocampus, were reduced in Type 1 diabetic rats. We detected Fe²⁺ level and lipid peroxidation products (malondialdehyde (MDA) and 4-Hydroxynonenal (4-HNE)) in the hippocampus. Mitochondria and neuron injury in the STZ-induced diabetic rats were determined using a Transmission Electron Microscope and Nissl body staining. Iron overload and ferroptosis were detected in the hippocampus. Furthermore, mRNA microarray analysis revealed 201 dysregulated mRNAs in STZ-induced type 1 diabetes (T1D). Pathway enrichment analyses indicated that differentially expressed mRNAs associated-coding genes were associated with ferroptosis. Among ferroptosis signaling pathway genes, Slc40a1 gene (ferroportin) was downregulated. We show that ferroptosis is associated with diabetic cognitive dysfunction and Slc40a1 mediates ferroptosis in T1D.

A covariate-constraint method to map brain feature space into lower dimensional manifolds

Human brain connectome studies aim to both explore healthy brains, and extract and analyze relevant features associated with pathologies of interest. Usually this consists of modeling the brain connectome as a graph and using graph metrics as features. A fine brain description requires graph metrics computation at the node level. Given the relatively reduced number of patients in standard cohorts, such data analysis problems fall in the high-dimension, low-sample-size framework. In this context, our goal is to provide a machine learning technique that exhibits flexibility, gives the investigator an understanding of the features and covariates, allows visualization and exploration, and yields insight into the data and the biological phenomena at stake. The retained approach is dimension reduction in a manifold learning methodology; the originality is that the investigator chooses one (or several) reduced variables. The proposed method is illustrated in two studies. The first one addresses comatose patients; the second one compares young and elderly populations. The method sheds light on the differences between brain connectivity graphs using graph metrics and potential clinical interpretations of these differences.

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

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

Propagating patterns of intrinsic activity along macroscale gradients coordinate functional connections across the whole brain

The intrinsic activity of the human brain, observed with resting-state fMRI (rsfMRI) and functional connectivity, exhibits macroscale spatial organization such as functional networks and gradients. Dynamic analysis techniques have shown that functional connectivity is a mere summary of time-varying patterns with distinct spatial and temporal characteristics. A better understanding of these patterns might provide insight into aspects of the brain's intrinsic activity that cannot be inferred by functional connectivity or the spatial maps derived from it, such as functional networks and gradients. Here, we describe three spatiotemporal patterns of coordinated activity across the whole brain obtained by averaging similar ~20-second-long segments of rsfMRI timeseries. In each of these patterns, activity propagates along a particular macroscale functional gradient, simultaneously across the cerebral cortex and in most other brain regions. In some regions, like the thalamus, the propagation suggests previously-undescribed gradients. The coordinated activity across areas is consistent with known tract-based connections, and nuanced differences in the timing of peak activity between regions point to plausible driving mechanisms. The magnitude of correlation within and particularly between functional networks is remarkably diminished when these patterns are regressed from the rsfMRI timeseries, a quantitative demonstration of their significant role in functional connectivity. Taken together, our results suggest that a few recurring patterns of propagating intrinsic activity along macroscale gradients give rise to and coordinate functional connections across the whole brain.

Systematically disrupted functional gradient of the cortical connectome in generalized epilepsy: Initial discovery and independent sample replication

Genetic generalized epilepsy is a network disorder typically involving distributed areas identified by classical neuroanatomy. However, the finer topological relationships in terms of continuous spatial arrangement between these systems are still ambiguous. Connectome gradients provide the topological representations of human macroscale hierarchy in an abstract low-dimensional space by embedding the functional connectome into a set of axes. Leveraging connectome gradients, we systematically scrutinized abnormalities of functional connectome gradient in patients with genetic generalized epilepsy with tonic-clonic seizure (GGE-GTCS, n = 78) compared to healthy controls (HC, n = 85), and further examined the reproducibility across multiple processing configurations and in an independent validation sample (patients with GGE-GTCS, n = 28; HC, n = 31). Our findings demonstrated an extended principal gradient at different spatial scales, network-level and vertex-level, in patients with GGE-GTCS. We found consistent results across processing parameters and in validation sample. The extended principal gradient revealed the excessive functional segregation between unimodal and transmodal systems associated with duration of epilepsy and age at seizure onset in patients. Furthermore, the connectivity profile of regions with abnormal principal gradients verified the disrupted functional hierarchy revealed by gradients. Together, our findings provided a novel view of functional system hierarchy alterations, which facilitated a continuous spatial arrangement of macroscale networks, to increase our understanding of the functional connectome hierarchy in generalized epilepsy.

A Common Space Approach to Comparative Neuroscience

Comparative neuroscience is entering the era of big data. New high-throughput methods and data-sharing initiatives have resulted in the availability of large, digital data sets containing many types of data from ever more species. Here, we present a framework for exploiting the new possibilities offered. The multimodality of the data allows vertical translations, which are comparisons of different aspects of brain organization within a single species and across scales. Horizontal translations compare particular aspects of brain organization across species, often by building abstract feature spaces. Combining vertical and horizontal translations allows for more sophisticated comparisons, including relating principles of brain organization across species by contrasting horizontal translations, and for making formal predictions of unobtainable data based on observed results in a model species.

Cross-species neuroscience: closing the explanatory gap

Neuroscience has seen substantial development in non-invasive methods available for investigating the living human brain. However, these tools are limited to coarse macroscopic measures of neural activity that aggregate the diverse responses of thousands of cells. To access neural activity at the cellular and circuit level, researchers instead rely on invasive recordings in animals. Recent advances in invasive methods now permit large-scale recording and circuit-level manipulations with exquisite spatio-temporal precision. Yet, there has been limited progress in relating these microcircuit measures to complex cognition and behaviour observed in humans. Contemporary neuroscience thus faces an explanatory gap between macroscopic descriptions of the human brain and microscopic descriptions in animal models. To close the explanatory gap, we propose adopting a cross-species approach. Despite dramatic differences in the size of mammalian brains, this approach is broadly justified by preserved homology. Here, we outline a three-armed approach for effective cross-species investigation that highlights the need to translate different measures of neural activity into a common space. We discuss how a cross-species approach has the potential to transform basic neuroscience while also benefiting neuropsychiatric drug development where clinical translation has, to date, seen minimal success. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Estimation of simultaneous BOLD and dynamic FDG metabolic brain activations using a multimodality concatenated ICA (mcICA) method

Simultaneous magnetic resonance and positron emission tomography provides an opportunity to measure brain haemodynamics and metabolism in a single scan session, and to identify brain activations from multimodal measurements in response to external stimulation. However, there are few analysis methods available for jointly analysing the simultaneously acquired blood-oxygen-level dependant functional MRI (fMRI) and 18-F-fluorodeoxyglucose functional PET (fPET) datasets. In this work, we propose a new multimodality concatenated ICA (mcICA) method to identify joint fMRI-fPET brain activations in response to a visual stimulation task. The mcICA method produces a fused map from the multimodal datasets with equal contributions of information from both modalities, measured by entropy. We validated the method in silico, and applied it to an in vivo visual stimulation experiment. The mcICA method estimated the activated brain regions in the visual cortex modulated by both BOLD and FDG signals. The mcICA provides a fully data-driven analysis approach to analyse cerebral haemodynamic response and glucose uptake signals arising from exogenously induced neuronal activity.

Toward a connectivity gradient-based framework for reproducible biomarker discovery

Despite myriad demonstrations of feasibility, the high dimensionality of fMRI data remains a critical barrier to its utility for reproducible biomarker discovery. Recent efforts to address this challenge have capitalized on dimensionality reduction techniques applied to resting-state fMRI, identifying principal components of intrinsic connectivity which describe smooth transitions across different cortical systems, so called "connectivity gradients". These gradients recapitulate neurocognitively meaningful organizational principles that are present in both human and primate brains, and also appear to differ among individuals and clinical populations. Here, we provide a critical assessment of the suitability of connectivity gradients for biomarker discovery. Using the Human Connectome Project (discovery subsample=209; two replication subsamples= 209 × 2) and the Midnight scan club (n = 9), we tested the following key biomarker traits - reliability, reproducibility and predictive validity - of functional gradients. In doing so, we systematically assessed the effects of three analytical settings, including i) dimensionality reduction algorithms (i.e., linear vs. non-linear methods), ii) input data types (i.e., raw time series, [un-]thresholded functional connectivity), and iii) amount of the data (resting-state fMRI time-series lengths). We found that the reproducibility of functional gradients across algorithms and subsamples is generally higher for those explaining more variances of whole-brain connectivity data, as well as those having higher reliability. Notably, among different analytical settings, a linear dimensionality reduction (principal component analysis in our study), more conservatively thresholded functional connectivity (e.g., 95-97%) and longer time-series data (at least ≥20mins) was found to be preferential conditions to obtain higher reliability. Those gradients with higher reliability were able to predict unseen phenotypic scores with a higher accuracy, highlighting reliability as a critical prerequisite for validity. Importantly, prediction accuracy with connectivity gradients exceeded that observed with more traditional edge-based connectivity measures, suggesting the added value of a low-dimensional and multivariate gradient approach. Finally, the present work highlights the importance and benefits of systematically exploring the parameter space for new imaging methods before widespread deployment.

Cross-species functional alignment reveals evolutionary hierarchy within the connectome

Evolution provides an important window into how cortical organization shapes function and vice versa. The complex mosaic of changes in brain morphology and functional organization that have shaped the mammalian cortex during evolution, complicates attempts to chart cortical differences across species. It limits our ability to fully appreciate how evolution has shaped our brain, especially in systems associated with unique human cognitive capabilities that lack anatomical homologues in other species. Here, we develop a function-based method for cross-species alignment that enables the quantification of homologous regions between humans and rhesus macaques, even when their location is decoupled from anatomical landmarks. Critically, we find cross-species similarity in functional organization reflects a gradient of evolutionary change that decreases from unimodal systems and culminates with the most pronounced changes in posterior regions of the default mode network (angular gyrus, posterior cingulate and middle temporal cortices). Our findings suggest that the establishment of the default mode network, as the apex of a cognitive hierarchy, has changed in a complex manner during human evolution - even within subnetworks.

Joint embedding: A scalable alignment to compare individuals in a connectivity space

A common coordinate space enabling comparison across individuals is vital to understanding human brain organization and individual differences. By leveraging dimensionality reduction algorithms, high-dimensional fMRI data can be represented in a low-dimensional space to characterize individual features. Such a representative space encodes the functional architecture of individuals and enables the observation of functional changes across time. However, determining comparable functional features across individuals in resting-state fMRI in a way that simultaneously preserves individual-specific connectivity structure can be challenging. In this work we propose scalable joint embedding to simultaneously embed multiple individual brain connectomes within a common space that allows individual representations across datasets to be aligned. Using Human Connectome Project data, we evaluated the joint embedding approach by comparing it to the previously established orthonormal alignment model. Alignment using joint embedding substantially increased the similarity of functional representations across individuals while simultaneously capturing their distinct profiles, allowing individuals to be more discriminable from each other. Additionally, we demonstrated that the common space established using resting-state fMRI provides a better overlap of task-activation across participants. Finally, in a more challenging scenario - alignment across a lifespan cohort aged from 6 to 85 - joint embedding provided a better prediction of age (r2 = 0.65) than the prior alignment model. It facilitated the characterization of functional trajectories across lifespan. Overall, these analyses establish that joint embedding can simultaneously capture individual neural representations in a common connectivity space aligning functional data across participants and populations and preserve individual specificity.

Dispersion of functional gradients across the adult lifespan

Ageing is commonly associated with changes to segregation and integration of functional brain networks, but, in isolation, current network-based approaches struggle to elucidate changes across the many axes of functional organisation. However, the advent of gradient mapping techniques in neuroimaging provides a new means of studying functional organisation in a multi-dimensional connectivity space. Here, we studied ageing and behaviourally-relevant differences in a three-dimensional connectivity space using the Cambridge Centre for Ageing Neuroscience cohort (n = 643). Building on gradient mapping techniques, we developed a set of measures to quantify the dispersion within and between functional communities. We detected a strong shift of the visual network across the adult lifespan from an extreme to a more central position in the 3D gradient space. In contrast, the dispersion distance between transmodal communities (dorsal attention, ventral attention, frontoparietal and default mode) did not change. However, these communities themselves were increasingly dispersed with increasing age, reflecting more dissimilar functional connectivity profiles within each community. Increasing dispersion of frontoparietal, attention and default mode networks, in particular, were associated negatively with cognition, measured by fluid intelligence. By using a technique that explicitly captures the ordering of functional systems in a multi-dimensional hierarchical framework, we identified behaviorally-relevant age-related differences of within and between network organisation. We propose that the study of functional gradients across the adult lifespan could provide insights that may facilitate the development of new strategies to maintain cognitive ability across the lifespan in health and disease.

Cortical gradients of functional connectivity are robust to state-dependent changes following sleep deprivation

Sleep deprivation leads to significant impairments in cognitive performance and changes to the interactions between large scale cortical networks, yet the hierarchical organization of cortical activity across states is still being explored. We used functional magnetic resonance imaging to assess activations and connectivity during cognitive tasks in 20 healthy young adults, during three states: (i) following a normal night of sleep, (ii) following 24hr of total sleep deprivation, and (iii) after a morning recovery nap. Situating cortical activity during cognitive tasks along hierarchical organizing gradients based upon similarity of functional connectivity patterns, we found that regional variations in task-activations were captured by an axis differentiating areas involved in executive control from default mode regions and paralimbic cortex. After global signal regression, the range of functional differentiation along this axis at baseline was significantly related to decline in working memory performance (2-back task) following sleep deprivation, as well as the extent of recovery in performance following a nap. The relative positions of cortical regions within gradients did not significantly change across states, except for a lesser differentiation of the visual system and increased coupling of the posterior cingulate cortex with executive control areas after sleep deprivation. This was despite a widespread increase in the magnitude of functional connectivity across the cortex following sleep deprivation. Cortical gradients of functional differentiation thus appear relatively insensitive to state-dependent changes following sleep deprivation and recovery, suggesting that there are no large-scale changes in cortical functional organization across vigilance states. Certain features of particular gradient axes may be informative for the extent of decline in performance on more complex tasks following sleep deprivation, and could be beneficial over traditional voxel- or parcel-based approaches in identifying realtionships between state-dependent brain activity and behavior.

High-dimensional brain-wide functional connectivity mapping in magnetoencephalography

BACKGROUND

Brain functional connectivity (FC) analyses based on magneto/electroencephalography (M/EEG) signals have yet to exploit the intrinsic high-dimensional information. Typically, these analyses are constrained to regions of interest to avoid the curse of dimensionality, with the latter leading to conservative hypothesis testing.

NEW METHOD

We removed such constraint by estimating high-dimensional source-based M/EEG-FC using cluster-permutation statistic (CPS) and demonstrated the feasibility of this approach by identifying resting-state changes in mild cognitive impairment (MCI), a prodromal stage of Alzheimer's disease. Particularly, we proposed a unified framework for CPS analysis together with a novel neighbourhood measure to estimate more compact and neurophysiological plausible neural communication. As clusters could more confidently reveal interregional communication, we proposed and tested a cluster-strength index to demonstrate other advantages of CPS analysis.

RESULTS

We found clusters of increased communication or hypersynchronization in MCI compared to healthy controls in delta (1-4 Hz) and higher-theta (6-8 Hz) bands oscillations. These mainly consisted of interactions between occipitofrontal and occipitotemporal regions in the left hemisphere, which may be critically affected in the early stages of Alzheimer's disease.

CONCLUSIONS

Our approach could be important to create high-resolution FC maps from neuroimaging studies in general, allowing the multimodal analysis of neural communication across multiple spatial scales. Particularly, FC clusters more robustly represent the interregional communication by identifying dense bundles of connections that are less sensitive to inter-individual anatomical and functional variability. Overall, this approach could help to better understand neural information processing in healthy and disease conditions as needed for developing biomarker research.

A multi-scale cortical wiring space links cellular architecture and functional dynamics in the human brain

The vast net of fibres within and underneath the cortex is optimised to support the convergence of different levels of brain organisation. Here, we propose a novel coordinate system of the human cortex based on an advanced model of its connectivity. Our approach is inspired by seminal, but so far largely neglected models of cortico-cortical wiring established by postmortem anatomical studies and capitalises on cutting-edge in vivo neuroimaging and machine learning. The new model expands the currently prevailing diffusion magnetic resonance imaging (MRI) tractography approach by incorporation of additional features of cortical microstructure and cortico-cortical proximity. Studying several datasets and different parcellation schemes, we could show that our coordinate system robustly recapitulates established sensory-limbic and anterior-posterior dimensions of brain organisation. A series of validation experiments showed that the new wiring space reflects cortical microcircuit features (including pyramidal neuron depth and glial expression) and allowed for competitive simulations of functional connectivity and dynamics based on resting-state functional magnetic resonance imaging (rs-fMRI) and human intracranial electroencephalography (EEG) coherence. Our results advance our understanding of how cell-specific neurobiological gradients produce a hierarchical cortical wiring scheme that is concordant with increasing functional sophistication of human brain organisation. Our evaluations demonstrate the cortical wiring space bridges across scales of neural organisation and can be easily translated to single individuals.

Signal diffusion along connectome gradients and inter-hub routing differentially contribute to dynamic human brain function

Human cognition is dynamic, alternating over time between externally-focused states and more abstract, often self-generated, patterns of thought. Although cognitive neuroscience has documented how networks anchor particular modes of brain function, mechanisms that describe transitions between distinct functional states remain poorly understood. Here, we examined how time-varying changes in brain function emerge within the constraints imposed by macroscale structural network organization. Studying a large cohort of healthy adults (n = 326), we capitalized on manifold learning techniques that identify low dimensional representations of structural connectome organization and we decomposed neurophysiological activity into distinct functional states and their transition patterns using Hidden Markov Models. Structural connectome organization predicted dynamic transitions anchored in sensorimotor systems and those between sensorimotor and transmodal states. Connectome topology analyses revealed that transitions involving sensorimotor states traversed short and intermediary distances and adhered strongly to communication mechanisms of network diffusion. Conversely, transitions between transmodal states involved spatially distributed hubs and increasingly engaged long-range routing. These findings establish that the structure of the cortex is optimized to allow neural states the freedom to vary between distinct modes of processing, and so provides a key insight into the neural mechanisms that give rise to the flexibility of human cognition.