A graded tractographic parcellation of the temporal lobe

Graded and sharp transitions in semantic function in left temporal lobe

Recent work has focussed on how patterns of functional change within the temporal lobe relate to whole-brain dimensions of intrinsic connectivity variation (Margulies et al., 2016). We examined two such 'connectivity gradients' reflecting the separation of (i) unimodal versus heteromodal and (ii) visual versus auditory-motor cortex, examining visually presented verbal associative and feature judgments, plus picture-based context and emotion generation. Functional responses along the first dimension sometimes showed graded change between modality-tuned and heteromodal cortex (in the verbal matching task), and other times showed sharp functional transitions, with deactivation at the extremes and activation in the middle of this gradient (internal generation). The second gradient revealed more visual than auditory-motor activation, regardless of content (associative, feature, context, emotion) or task process (matching/generation). We also uncovered subtle differences across each gradient for content type, which predominantly manifested as differences in relative magnitude of activation or deactivation.

Exploring graded profiles of hippocampal atrophy along the anterior-posterior axis in semantic dementia and Alzheimer's disease

Mounting evidence indicates marked hippocampal degeneration in semantic dementia (SD) however, the spatial distribution of hippocampal atrophy profiles in this syndrome remains unclear. Using a recently developed parcellation approach, we extracted hippocampal volumes from four distinct subregions running from anterior to posterior along the longitudinal axis (anterior, intermediate rostral, intermediate caudal, and posterior). Volumetric differences in hippocampal subregions were compared between 21 SD, 24 matched Alzheimer's disease (AD), and 27 healthy older Control participants. Despite comparable overall hippocampal volume loss, SD and AD groups diverged in terms of the magnitude of atrophy along the anterior-posterior axis of the hippocampus. Global hippocampal atrophy was observed in AD, with no discernible gradation or lateralisation. In contrast, SD patients displayed graded bilateral hippocampal atrophy, most pronounced on the left-hand side, and concentrated in anterior relative to posterior subregions. Finally, we found preliminary evidence that disease-specific vulnerability along the anterior-posterior axis of the hippocampus was associated with canonical clinical features of these syndromes.

Feature similarity gradients detect alterations in the neonatal cortex associated with preterm birth

The early life environment programmes cortical architecture and cognition across the life course. A measure of cortical organisation that integrates information from multimodal MRI and is unbound by arbitrary parcellations has proven elusive, which hampers efforts to uncover the perinatal origins of cortical health. Here, we use the Vogt-Bailey index to provide a fine-grained description of regional homogeneities and sharp variations in cortical microstructure based on feature gradients, and we investigate the impact of being born preterm on cortical development at term-equivalent age. Compared with term-born controls, preterm infants have a homogeneous microstructure in temporal and occipital lobes, and the medial parietal, cingulate, and frontal cortices, compared with term infants. These observations replicated across two independent datasets and were robust to differences that remain in the data after matching samples and alignment of processing and quality control strategies. We conclude that cortical microstructural architecture is altered in preterm infants in a spatially distributed rather than localised fashion.

Impact of brain parcellation on prediction performance in models of cognition and demographics

Brain connectivity analysis begins with the selection of a parcellation scheme that will define brain regions as nodes of a network whose connections will be studied. Brain connectivity has already been used in predictive modelling of cognition, but it remains unclear if the resolution of the parcellation used can systematically impact the predictive model performance. In this work, structural, functional and combined connectivity were each defined with five different parcellation schemes. The resolution and modality of the parcellation schemes were varied. Each connectivity defined with each parcellation was used to predict individual differences in age, education, sex, executive function, self-regulation, language, encoding and sequence processing. It was found that low-resolution functional parcellation consistently performed above chance at producing generalisable models of both demographics and cognition. However, no single parcellation scheme showed a superior predictive performance across all cognitive domains and demographics. In addition, although parcellation schemes impacted the graph theory measures of each connectivity type (structural, functional and combined), these differences did not account for the out-of-sample predictive performance of the models. Taken together, these findings demonstrate that while high-resolution parcellations may be beneficial for modelling specific individual differences, partial voluming of signals produced by the higher resolution of the parcellation likely disrupts model generalisability.

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.

The hierarchical organization of the precuneus captured by functional gradients

The precuneus shows considerable heterogeneity in multiple dimensions including anatomy, function, and involvement in brain disorders. Leveraging the state-of-the-art functional gradient approach, we aimed to investigate the hierarchical organization of the precuneus, which may hold promise for a unified understanding of precuneus heterogeneity. Resting-state functional MRI data from 793 healthy individuals were used to discover and validate functional gradients of the precuneus, which were calculated based on the voxel-wise precuneus-to-cerebrum functional connectivity patterns. Then, we further explored the potential relationships of the precuneus functional gradients with cortical morphology, intrinsic geometry, canonical functional networks, and behavioral domains. We found that the precuneus principal and secondary gradients showed dorsoanterior-ventral and ventroposterior-dorsal organizations, respectively. Concurrently, the principal gradient was associated with cortical morphology, and both the principal and secondary gradients showed geometric distance dependence. Importantly, precuneus functional subdivisions corresponding to canonical functional networks (behavioral domains) were distributed along both gradients in a hierarchical manner, i.e., from the sensorimotor network (somatic movement and sensation) at one extreme to the default mode network (abstract cognitive functions) at the other extreme for the principal gradient and from the visual network (vision) at one end to the dorsal attention network (top-down control of attention) at the other end for the secondary gradient. These findings suggest that the precuneus functional gradients may provide mechanistic insights into the multifaceted nature of precuneus heterogeneity.

Functional connectivity gradients of the cingulate cortex

Heterogeneity of the cingulate cortex is evident in multiple dimensions including anatomy, function, connectivity, and involvement in networks and diseases. Using the recently developed functional connectivity gradient approach and resting-state functional MRI data, we found three functional connectivity gradients that captured distinct dimensions of cingulate hierarchical organization. The principal gradient exhibited a radiating organization with transitions from the middle toward both anterior and posterior parts of the cingulate cortex and was related to canonical functional networks and corresponding behavioral domains. The second gradient showed an anterior-posterior axis across the cingulate cortex and had prominent geometric distance dependence. The third gradient displayed a marked differentiation of subgenual and caudal middle with other parts of the cingulate cortex and was associated with cortical morphology. Aside from providing an updated framework for understanding the multifaceted nature of cingulate heterogeneity, the observed hierarchical organization of the cingulate cortex may constitute a novel research agenda with potential applications in basic and clinical neuroscience.

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.

The convergence and divergence of episodic and semantic functions across lateral parietal cortex

Decades of research have highlighted the importance of lateral parietal cortex (LPC) across a myriad of cognitive domains. Yet, the underlying function of LPC remains unclear. Two domains that have emphasized LPC involvement are semantic memory and episodic memory retrieval. From each domain, sophisticated functional models have been proposed, as well as the more domain-general assumption that LPC is engaged by any form of internally directed cognition (episodic/semantic retrieval being examples). Here we used a combination of functional magnetic resonance imaging, functional connectivity, and diffusion tensor imaging white-matter connectivity to show that (i) ventral LPC (angular gyrus [AG]) was positively engaged during episodic retrieval but disengaged during semantic memory retrieval and (ii) activity negatively varied with task difficulty in the semantic task whereas episodic activation was independent of difficulty. In contrast, dorsal LPC (intraparietal sulcus) showed domain general activation that was positively correlated with task difficulty. Finally, (iii) a dorsal-ventral and anterior-posterior gradient of functional and structural connectivity was found across the AG (e.g. mid-AG connected with episodic retrieval). We propose a unifying model in which LPC as a whole might share a common underlying neurocomputation (multimodal buffering) with variations in the emergent cognitive functions across subregions arising from differences in the underlying connectivity.

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 cross-domain functional organization of posterior lateral temporal cortex: insights from ALE meta-analyses of 7 cognitive domains spanning 12,000 participants

The posterior lateral temporal cortex is implicated in many verbal, nonverbal, and social cognitive domains and processes. Yet without directly comparing these disparate domains, the region's organization remains unclear; do distinct processes engage discrete subregions, or could different domains engage shared neural correlates and processes? Here, using activation likelihood estimation meta-analyses, the bilateral posterior lateral temporal cortex subregions engaged in 7 domains were directly compared. These domains comprised semantics, semantic control, phonology, biological motion, face processing, theory of mind, and representation of tools. Although phonology and biological motion were predominantly associated with distinct regions, other domains implicated overlapping areas, perhaps due to shared underlying processes. Theory of mind recruited regions implicated in semantic representation, tools engaged semantic control areas, and faces engaged subregions for biological motion and theory of mind. This cross-domain approach provides insight into how posterior lateral temporal cortex is organized and why.

Functional connectivity gradients of the insula to different cerebral systems

The diverse functional roles of the insula may emerge from its heavy connectivity to an extensive network of cortical and subcortical areas. Despite several previous attempts to investigate the hierarchical organization of the insula by applying the recently developed gradient approach to insula-to-whole brain connectivity data, little is known about whether and how there is variability across connectivity gradients of the insula to different cerebral systems. Resting-state functional MRI data from 793 healthy subjects were used to discover and validate functional connectivity gradients of the insula, which were computed based on its voxel-wise functional connectivity profiles to distinct cerebral systems. We identified three primary patterns of functional connectivity gradients of the insula to distinct cerebral systems. The connectivity gradients to the higher-order transmodal associative systems, including the prefrontal, posterior parietal, temporal cortices, and limbic lobule, showed a ventroanterior-dorsal axis across the insula; those to the lower-order unimodal primary systems, including the motor, somatosensory, and occipital cortices, displayed radiating transitions from dorsoanterior toward both ventroanterior and dorsoposterior parts of the insula; the connectivity gradient to the subcortical nuclei exhibited an organization along the anterior-posterior axis of the insula. Apart from complementing and extending previous literature on the heterogeneous connectivity patterns of insula subregions, the presented framework may offer ample opportunities to refine our understanding of the role of the insula in many brain disorders.

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/).

Subregions of DLPFC Display Graded yet Distinct Structural and Functional Connectivity

The human dorsolateral prefrontal cortex (DLPFC; approximately corresponding to Brodmann areas 9 and 46) has demonstrable roles in diverse executive functions such as working memory, cognitive flexibility, planning, inhibition, and abstract reasoning. However, it remains unclear whether this is the result of one functionally homogeneous region or whether there are functional subdivisions within the DLPFC. Here, we divided the DLPFC into seven areas along rostral-caudal and dorsal-ventral axes anatomically and explored their respective patterns of structural and functional connectivity. In vivo probabilistic tractography (11 females and 13 males) and resting-state functional magnetic resonance imaging (fMRI; 57 females and 21 males) were employed to map out the patterns of connectivity from each DLPFC subregion. Structural connectivity demonstrated graded intraregional connectivity within the DLPFC. The patterns of structural connectivity between the DLPFC subregions and other cortical areas revealed that the dorsal-rostral subregions connections were restricted to other frontal and limbic areas, whereas the ventral-caudal region was widely connected to frontal, parietal, and limbic cortex. Functional connectivity analyses demonstrated that subregions of DLPFC were strongly interconnected to each other. The dorsal subregions were associated with the default mode network (DMN), while middle dorsal-rostral subregions were linked with the multiple demand network (MDN). The ventral-caudal subregion showed increased functional coupling with both DMN and MDN. Our results suggest that the connectivity of the DLPFC may be subdivided along a dorsorostral-ventrocaudal axis with differing (albeit graded) patterns of connectivity reflecting the integrative executive function of the DLPFC.SIGNIFICANCE STATEMENT Research has shown that the dorsolateral prefrontal cortex (DLPFC) plays a role in various executive functions. By dividing the DLPFC into seven areas along rostral-caudal and dorsal-ventral axes anatomically, we explored their patterns of structural and functional connectivity. The patterns of connectivity within DLPFC subregions demonstrated graded intraregional connectivity. There were distinctive patterns of connectivity with other cortical areas in dorsal-rostral and ventral-caudal DLPFC subregions. Divisions across DLPFC subregions seem to align with their structural and functional connectivity. Our results suggest that DLPFC may be subdivided by the diagonal axis of the dorsal-ventral axis and rostral-caudal axis, supporting the framework of a functional organization along the anterior-posterior axis in the lateral PFC.

Learning Cortical Parcellations Using Graph Neural Networks

Deep learning has been applied to magnetic resonance imaging (MRI) for a variety of purposes, ranging from the acceleration of image acquisition and image denoising to tissue segmentation and disease diagnosis. Convolutional neural networks have been particularly useful for analyzing MRI data due to the regularly sampled spatial and temporal nature of the data. However, advances in the field of brain imaging have led to network- and surface-based analyses that are often better represented in the graph domain. In this analysis, we propose a general purpose cortical segmentation method that, given resting-state connectivity features readily computed during conventional MRI pre-processing and a set of corresponding training labels, can generate cortical parcellations for new MRI data. We applied recent advances in the field of graph neural networks to the problem of cortical surface segmentation, using resting-state connectivity to learn discrete maps of the human neocortex. We found that graph neural networks accurately learn low-dimensional representations of functional brain connectivity that can be naturally extended to map the cortices of new datasets. After optimizing over algorithm type, network architecture, and training features, our approach yielded mean classification accuracies of 79.91% relative to a previously published parcellation. We describe how some hyperparameter choices including training and testing data duration, network architecture, and algorithm choice affect model performance.

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.

A tale of two gradients: differences between the left and right hemispheres predict semantic cognition

Decomposition of whole-brain functional connectivity patterns reveals a principal gradient that captures the separation of sensorimotor cortex from heteromodal regions in the default mode network (DMN). Functional homotopy is strongest in sensorimotor areas, and weakest in heteromodal cortices, suggesting there may be differences between the left and right hemispheres (LH/RH) in the principal gradient, especially towards its apex. This study characterised hemispheric differences in the position of large-scale cortical networks along the principal gradient, and their functional significance. We collected resting-state fMRI and semantic, working memory and non-verbal reasoning performance in 175 + healthy volunteers. We then extracted the principal gradient of connectivity for each participant, tested which networks showed significant hemispheric differences on the gradient, and regressed participants’ behavioural efficiency in tasks outside the scanner against interhemispheric gradient differences for each network. LH showed a higher overall principal gradient value, consistent with its role in heteromodal semantic cognition. One frontotemporal control subnetwork was linked to individual differences in semantic cognition: when it was nearer heteromodal DMN on the principal gradient in LH, participants showed more efficient semantic retrieval—and this network also showed a strong hemispheric difference in response to semantic demands but not working memory load in a separate study. In contrast, when a dorsal attention subnetwork was closer to the heteromodal end of the principal gradient in RH, participants showed better visual reasoning. Lateralization of function may reflect differences in connectivity between control and heteromodal regions in LH, and attention and visual regions in RH.

Does the temporal cortex make us human? A review of structural and functional diversity of the primate temporal lobe

Temporal cortex is a primate specialization that shows considerable variation in size, morphology, and connectivity across species. Human temporal cortex is involved in many behaviors that are considered especially well developed in humans, including semantic processing, language, and theory of mind. Here, we ask whether the involvement of temporal cortex in these behaviors can be explained in the context of the 'general' primate organization of the temporal lobe or whether the human temporal lobe contains unique specializations indicative of a 'step change' in the lineage leading to modern humans. We propose that many human behaviors can be explained as elaborations of temporal cortex functions observed in other primates. However, changes in temporal lobe white matter suggest increased integration of information within temporal cortex and between posterior temporal cortex and other association areas, which likely enable behaviors not possible in other species.

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.

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.

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.

A tutorial and tool for exploring feature similarity gradients with MRI data

There has been an increasing interest in examining organisational principles of the cerebral cortex (and subcortical regions) using different MRI features such as structural or functional connectivity. Despite the widespread interest, introductory tutorials on the underlying technique targeted for the novice neuroimager are sparse in the literature. Articles that investigate various "neural gradients" (for example based on region studied "cortical gradients," "cerebellar gradients," "hippocampal gradients" etc … or feature of interest "functional gradients," "cytoarchitectural gradients," "myeloarchitectural gradients" etc …) have increased in popularity. Thus, we believe that it is opportune to discuss what is generally meant by "gradient analysis". We introduce basics concepts in graph theory, such as graphs themselves, the degree matrix, and the adjacency matrix. We discuss how one can think about gradients of feature similarity (the similarity between timeseries in fMRI, or streamline in tractography) using graph theory and we extend this to explore such gradients across the whole MRI scale; from the voxel level to the whole brain level. We proceed to introduce a measure for quantifying the level of similarity in regions of interest. We propose the term "the Vogt-Bailey index" for such quantification to pay homage to our history as a brain mapping community. We run through the techniques on sample datasets including a brain MRI as an example of the application of the techniques on real data and we provide several appendices that expand upon details. To maximise intuition, the appendices contain a didactic example describing how one could use these techniques to solve a particularly pernicious problem that one may encounter at a wedding. Accompanying the article is a tool, available in both MATLAB and Python, that enables readers to perform the analysis described in this article on their own data. We refer readers to the graphical abstract as an overview of the analysis pipeline presented in this work.

The anterior and medial thalamic nuclei and the human limbic system: tracing the structural connectivity using diffusion-weighted imaging

The limbic system is a phylogenetically old, behaviorally defined system that serves as a center for emotions. It controls the expression of anger, fear, and joy and also influences sexual behavior, vegetative functions, and memory. The system comprises a collection of tel-, di-, and mesencephalic structures whose components have evolved and increased over time. Previous animal research indicates that the anterior nuclear group of the thalamus (ANT), as well as the habenula (Hb) and the adjacent mediodorsal nucleus (MD) each play a vital role in the limbic circuitry. Accordingly, diffusion imaging data of 730 subjects obtained from the Human Connectome Project and the masks of six nuclei (anterodorsal, anteromedial, anteroventral, lateral dorsal, Hb, and MD) served as seed regions for a direct probabilistic tracking to the rest of the brain using diffusion-weighted imaging. The results revealed that the ANT nuclei are part of the limbic and the memory system as they mainly connect via the mammillary tract, mammillary body, anterior commissure, fornix, and retrosplenial cortices to the hippocampus, amygdala, medio-temporal, orbito-frontal and occipital cortices. Furthermore, the ANT nuclei showed connections to the mesencephalon and brainstem to varying extents, a pattern rarely described in experimental findings. The habenula—usually defined as part of the epithalamus—was closely connected to the tectum opticum and seems to serve as a neuroanatomical hub between the visual and the limbic system, brainstem, and cerebellum. Finally, in contrast to experimental findings with tracer studies, directly determined connections of MD were mainly confined to the brainstem, while indirect MD fibers form a broad pathway connecting the hippocampus and medio-temporal areas with the mediofrontal cortex.

Connectivity Gradient in the Human Left Inferior Frontal Gyrus: Intraoperative Cortico-Cortical Evoked Potential Study

In the dual-stream model of language processing, the exact connectivity of the ventral stream to the anterior temporal lobe remains elusive. To investigate the connectivity between the inferior frontal gyrus (IFG) and the lateral part of the temporal and parietal lobes, we integrated spatiotemporal profiles of cortico-cortical evoked potentials (CCEPs) recorded intraoperatively in 14 patients who had undergone surgical resection for a brain tumor or epileptic focus. Four-dimensional visualization of the combined CCEP data showed that the pars opercularis (Broca’s area) is connected to the posterior temporal cortices and the supramarginal gyrus, whereas the pars orbitalis is connected to the anterior lateral temporal cortices and angular gyrus. Quantitative topographical analysis of CCEP connectivity confirmed an anterior–posterior gradient of connectivity from IFG stimulus sites to the temporal response sites. Reciprocality analysis indicated that the anterior part of the IFG is bidirectionally connected to the temporal or parietal area. This study shows that each IFG subdivision has different connectivity to the temporal lobe with an anterior–posterior gradient and supports the classical connectivity concept of Dejerine; that is, the frontal lobe is connected to the temporal lobe through the arcuate fasciculus and also a double fan-shaped structure anchored at the limen insulae.

Linking individual differences in semantic cognition to white matter microstructure

Semantic cognition is thought to involve the interaction of heteromodal conceptual representations with control processes that (i) focus retrieval on currently-relevant information, and (ii) suppress dominant yet irrelevant features and associations. Research suggests that semantic control demands are higher when retrieving a link between weakly-associated word pairs, since there is a mismatch between the pattern of semantic retrieval required by the task and the dominant associations of individual words. In addition, given that heteromodal concepts are thought to reflect the integration of vision, audition, valence and other features, the control demands of semantic tasks should be higher when there is less consistency between these features. In the present study, 62 volunteers completed a semantic decision task, where association strength and semantic-affective congruence were manipulated. We used diffusion tensor magnetic resonance imaging to obtain fractional anisotropy (FA) measures of white matter tracts hypothesized to be part of the semantic network. The behavioural data revealed an interaction between semantic-affective congruence and strength of association, suggesting these manipulations both contribute to semantic control demands. Next we considered how individual differences in these markers of semantic control relate to the microstructure of canonical white matter tracts, complementing previous studies that have largely focused on measures of intrinsic functional connectivity. Repeated-measures analysis of covariance showed opposing interactions between semantic control markers and FA of two tracts: left inferior longitudinal fasciculus (ILF) and right inferior fronto-occipital fasciculus (IFOF). Participants with higher FA in left ILF showed more efficient retrieval of weak associations, and more accurate performance for weak associations when meaning and valence were incongruent, consistent with the hypothesis that this left hemisphere tract supports semantic control. In contrast, participants with higher FA in right IFOF were more accurate for trials in which meaning and valence were congruent, and consequently when semantic control demands were minimised. These findings are consistent with recent studies showing that semantic control processes are strongly left-lateralised. In contrast, long-range connections from vision to semantic regions in the right hemisphere might support relatively automatic patterns of semantic retrieval.

Principles of temporal association cortex organisation as revealed by connectivity gradients

To establish the link between structure and function of any large area of the neocortex, it is helpful to identify its principles of organisation. One way to establish such principles is to investigate how differences in whole-brain connectivity are structured across the area. Here, we use Laplacian eigenmaps on diffusion MRI tractography data to investigate the organisational principles of the human temporal association cortex. We identify three overlapping gradients of connectivity that are, for the most part, consistent across hemispheres. The first gradient reveals an inferior–superior organisation of predominantly longitudinal tracts and separates visual and auditory unimodal and multimodal cortices. The second gradient radiates outward from the posterior middle temporal cortex with the arcuate fascicle as a distinguishing feature; the third gradient is concentrated in the anterior temporal lobe and emanates towards its posterior end. We describe the functional relevance of each of these gradients through the meta-analysis of data from the neuroimaging literature. Together, these results unravel the overlapping dimensions of structural organization of the human temporal cortex and provide a framework underlying its functional multiplicity.

CBPtools: a Python package for regional connectivity-based parcellation

Regional connectivity-based parcellation (rCBP) is a widely used procedure for investigating the structural and functional differentiation within a region of interest (ROI) based on its long-range connectivity. No standardized software or guidelines currently exist for applying rCBP, making the method only accessible to those who develop their own tools. As such, there exists a discrepancy between the laboratories applying the procedure each with their own software solutions, making it difficult to compare and interpret the results. Here, we outline an rCBP procedure accompanied by an open source software package called CBPtools. CBPtools is a Python (version 3.5+) package that allows users to run an extensively evaluated rCBP analysis workflow on a given ROI. It currently supports two modalities: resting-state functional connectivity and structural connectivity based on diffusion-weighted imaging, along with support for custom connectivity matrices. Analysis parameters are customizable and the workflow can be scaled to a large number of subjects using a parallel processing environment. Parcellation results with corresponding validity metrics are provided as textual and graphical output. Thus, CBPtools provides a simple plug-and-play, yet customizable way to conduct rCBP analyses. By providing an open-source software we hope to promote reproducible and comparable rCBP analyses and, importantly, make the rCBP procedure readily available. Here, we demonstrate the utility of CBPtools using a voluminous data set on an average compute-cluster infrastructure by performing rCBP on three ROIs prominently featured in parcellation literature.

BrainSpace: a toolbox for the analysis of macroscale gradients in neuroimaging and connectomics datasets

Understanding how cognitive functions emerge from brain structure depends on quantifying how discrete regions are integrated within the broader cortical landscape. Recent work established that macroscale brain organization and function can be described in a compact manner with multivariate machine learning approaches that identify manifolds often described as cortical gradients. By quantifying topographic principles of macroscale organization, cortical gradients lend an analytical framework to study structural and functional brain organization across species, throughout development and aging, and its perturbations in disease. Here, we present BrainSpace, a Python/Matlab toolbox for (i) the identification of gradients, (ii) their alignment, and (iii) their visualization. Our toolbox furthermore allows for controlled association studies between gradients with other brain-level features, adjusted with respect to null models that account for spatial autocorrelation. Validation experiments demonstrate the usage and consistency of our tools for the analysis of functional and microstructural gradients across different spatial scales.

A structural connectivity convergence zone in the ventral and anterior temporal lobes: Data-driven evidence from structural imaging

The hub-and-spoke model of semantic cognition seeks to reconcile embodied views of a fully distributed semantic network with patient evidence, primarily from semantic dementia, who demonstrate modality-independent conceptual deficits associated with atrophy centred on the ventrolateral anterior temporal lobe. The proponents of this model have recently suggested that the temporal cortex is a graded representational space where concepts become less linked to a specific modality as they are processed farther away from primary and secondary sensory cortices and towards the ventral anterior temporal lobe. To explore whether there is evidence that the connectivity patterns of the temporal lobe converge in its ventral anterior end the current study uses three dimensional Laplacian eigenmapping, a technique that allows visualisation of similarity in a low dimensional space. In this space similarity is encoded in terms of distances between data points. We found that the ventral and anterior temporal lobe is in a unique position of being at the centre of mass of the data points within the connective similarity space. This can be interpreted as the area where the connectivity profiles of all other temporal cortex voxels converge. This study is the first to explicitly investigate the pattern of connectivity and thus provides the missing link in the evidence that the ventral anterior temporal lobe can be considered a multi-modal graded hub.

LittleBrain: A gradient-based tool for the topographical interpretation of cerebellar neuroimaging findings

Gradient-based approaches to brain function have recently unmasked fundamental properties of brain organization. Diffusion map embedding analysis of resting-state fMRI data revealed a primary-to-transmodal axis of cerebral cortical macroscale functional organization. The same method was recently used to analyze resting-state data within the cerebellum, revealing for the first time a sensorimotor-fugal macroscale organization principle of cerebellar function. Cerebellar gradient 1 extended from motor to non-motor task-unfocused (default-mode network) areas, and cerebellar gradient 2 isolated task-focused processing regions. Here we present a freely available and easily accessible tool that applies this new knowledge to the topographical interpretation of cerebellar neuroimaging findings. LittleBrain illustrates the relationship between cerebellar data (e.g., volumetric patient study clusters, task activation maps, etc.) and cerebellar gradients 1 and 2. Specifically, LittleBrain plots all voxels of the cerebellum in a two-dimensional scatterplot, with each axis corresponding to one of the two principal functional gradients of the cerebellum, and indicates the position of cerebellar neuroimaging data within these two dimensions. This novel method of data mapping provides alternative, gradual visualizations that complement discrete parcellation maps of cerebellar functional neuroanatomy. We present application examples to show that LittleBrain can also capture subtle, progressive aspects of cerebellar functional neuroanatomy that would be difficult to visualize using conventional mapping techniques. Download and use instructions can be found at https://xaviergp.github.io/littlebrain.

Fiber length profiling: A novel approach to structural brain organization

There has been a recent increased interest in the structural connectivity of the cortex. However, an important feature of connectivity remains relatively unexplored; tract length. In this article, we develop an approach to characterize fiber length distributions across the human cerebral cortex. We used data from 76 participants of the Adult WU-Minn Human Connectome Project using probabilistic tractography. We found that connections of different lengths are not evenly distributed across the cortex. They form patterns where certain areas have a high density of fibers of a specific length while other areas have very low density. To assess the relevance of these new maps in relation to established characteristics, we compared them to structural indices such as cortical myelin content and cortical thickness. Additionally, we assessed their relation to resting state network organization. We noted that areas with very short fibers have relatively more myelin and lower cortical thickness while the pattern is inverted for longer fibers. Furthermore, the cortical fiber length distributions produce specific correlation patterns with functional resting state networks. Specifically, we find evidence that as resting state networks increase in complexity, their length profiles change. The functionally more complex networks correlate with maps of varying lengths while primary networks have more restricted correlations. We posit that these maps are a novel way of differentiating between 'local modules' that have restricted connections to 'neighboring' areas and 'functional integrators' that have more far reaching connectivity.

From percept to concept in the ventral temporal lobes: Graded hemispheric specialisation based on stimulus and task

The left and right ventral anterior temporal lobes (vATL) have been implicated as key regions for the representation of conceptual knowledge. However, the nature and degree of hemispheric specialisation in their function is unclear. To address this issue, we investigated hemispheric specialisation in the ventral temporal lobes using a distortion-corrected spin-echo fMRI protocol that enhanced signal in vATL. We employed an orthogonal manipulation of stimulus (written words vs pictured objects) and task (naming vs recognition). Words elicited left-lateralised vATL activation while objects elicited bilateral activation with no hemispheric bias. In contrast, posterior ventral temporal cortex exhibited a rightward bias for objects as well as a leftward bias for words. Naming tasks produced left-lateralised activation in vATL while activity for recognition was equal in left and right vATLs. These findings are incompatible with proposals that left and right ATLs are strongly modular in function, since these predict rightward as well as leftward biases. Instead, they support an alternative model in which (a) left and right ATL together form a bilateral, integrated system for the representation of concepts and (b) within this system, graded hemispheric specialisation emerges as a consequence of differential connectivity with other neural systems. On this view, greater left vATL activation for written word processing develops as a consequence of the inputs this region receives from left-lateralised visual word processing system in posterior temporal cortex. Greater left vATL activation during naming tasks is most likely due to connectivity with left-lateralised speech output systems in prefrontal and motor cortices.

Cortical parcellation based on structural connectivity: A case for generative models

One of the major challenges in systems neuroscience is to identify brain networks and unravel their significance for brain function -this has led to the concept of the 'connectome'. Connectomes are currently extensively studied in large-scale international efforts at multiple scales, and follow different definitions with respect to their connections as well as their elements. Perhaps the most promising avenue for defining the elements of connectomes originates from the notion that individual brain areas maintain distinct (long-range) connection profiles. These connectivity patterns determine the areas' functional properties and also allow for their anatomical delineation and mapping. This rationale has motivated the concept of connectivity-based cortex parcellation. In the past ten years, non-invasive mapping of human brain connectivity has led to immense advances in the development of parcellation techniques and their applications. Unfortunately, many of these approaches primarily aim for confirmation of well-known, existing architectonic maps and, to that end, unsuitably incorporate prior knowledge and frequently build on circular argumentation. Often, current approaches also tend to disregard the specific apertures of connectivity measurements, as well as the anatomical specificities of cortical areas, such as spatial compactness, regional heterogeneity, inter-subject variability, the multi-scaling nature of connectivity information, and potential hierarchical organisation. From a methodological perspective, however, a useful framework that regards all of these aspects in an unbiased way is technically demanding. In this commentary, we first outline the concept of connectivity-based cortex parcellation and discuss its prospects and limitations in particular with respect to structural connectivity. To improve reliability and efficiency, we then strongly advocate for connectivity-based cortex parcellation as a modelling approach; that is, an approximation of the data based on (model) parameter inference. As such, a parcellation algorithm can be formally tested for robustness -the precision of its predictions can be quantified and statistics about potential generalization of the results can be derived. Such a framework also allows the question of model constraints to be reformulated in terms of hypothesis testing through model selection and offers a formative way to integrate anatomical knowledge in terms of prior distributions.

An emergent functional parcellation of the temporal cortex

The temporal lobe has been associated with various cognitive functions which include memory, auditory cognition and semantics. However, at a higher level of conceptualisation, all of the functions associated with the temporal lobe can be considered as lying along one major axis; from modality-specific to modality-general processing. This paper used a spectral reordering technique on resting-state and task-based functional data to extract the major organisational axis of the temporal lobe in a bottom-up, data-driven fashion. Independent parcellations were performed on resting-state scans from 71 participants and active semantic task scans from 23 participants acquired using dual echo gradient echo planar imaging in order to preserve signal in inferior temporal cortex. The resulting organisational axis was consistent (over dataset and hemisphere) and progressed from superior temporal gyrus and posterior inferior temporal cortex to ventrolateral anterior temporal cortex. A hard parcellation separated a posterior (superior temporal and posterior fusiform and inferior temporal gyri) and an anterior cluster (ventrolateral anterior temporal lobe). The functional connectivity of the hard clusters supported the hypothesis that the connectivity gradient separated modality-specific and modality-general regions. This hypothesis was then directly tested by performing a VOI analysis upon an independent semantic task-based data set including auditory and visually presented stimuli. This confirmed that the ventrolateral anterior aspects of the temporal lobe are associated with modality-general processes whilst posterior and superior aspects are specific to certain modalities, with the posterior inferior subregions involved in visual processes and superior regions involved in audition.