Connectivity-based parcellation: Critique and implications

fMRI-based spatio-temporal parcellations of the human brain

PURPOSE OF REVIEW

Human brain parcellation based on functional magnetic resonance imaging (fMRI) plays an essential role in neuroscience research. By segmenting vast and intricate fMRI data into functionally similar units, researchers can better decipher the brain's structure in both healthy and diseased states. This article reviews current methodologies and ideas in this field, while also outlining the obstacles and directions for future research.

RECENT FINDINGS

Traditional brain parcellation techniques, which often rely on cytoarchitectonic criteria, overlook the functional and temporal information accessible through fMRI. The adoption of machine learning techniques, notably deep learning, offers the potential to harness both spatial and temporal information for more nuanced brain segmentation. However, the search for a one-size-fits-all solution to brain segmentation is impractical, with the choice between group-level or individual-level models and the intended downstream analysis influencing the optimal parcellation strategy. Additionally, evaluating these models is complicated by our incomplete understanding of brain function and the absence of a definitive "ground truth".

SUMMARY

While recent methodological advancements have significantly enhanced our grasp of the brain's spatial and temporal dynamics, challenges persist in advancing fMRI-based spatio-temporal representations. Future efforts will likely focus on refining model evaluation and selection as well as developing methods that offer clear interpretability for clinical usage, thereby facilitating further breakthroughs in our comprehension of the brain.

The connectivity-based architecture of the human piriform cortex

The anatomy of the human piriform cortex (PC) is poorly understood. We used a bimodal connectivity-based-parcellation approach to investigate subregions of the PC and its connectional differentiation from the amygdala. One hundred (55 % female) genetically unrelated subjects from the Human Connectome Project were included. A region of interest (ROI) was delineated bilaterally covering PC and amygdala, and functional and structural connectivity of this ROI with the whole gray matter was computed. Spectral clustering was performed to obtain bilateral parcellations at granularities of k = 2-10 clusters and combined bimodal parcellations were computed. Validity of parcellations was assessed via their mean individual-to-group similarity per adjusted rand index (ARI). Individual-to-group similarity was higher than chance in both modalities and in all clustering solutions. The amygdala was clearly distinguished from PC in structural parcellations, and olfactory amygdala was connectionally more similar to amygdala than to PC. At higher granularities, an anterior and ventrotemporal and a posterior frontal cluster emerged within PC, as well as an additional temporal cluster at their boundary. Functional parcellations also showed a frontal piriform cluster, and similar temporal clusters were observed with less consistency. Results from bimodal parcellations were similar to the structural parcellations. Consistent results were obtained in a validation cohort. Distinction of the human PC from the amygdala, including its olfactory subregions, is possible based on its structural connectivity alone. The canonical fronto-temporal boundary within PC was reproduced in both modalities and with consistency. All obtained parcellations are freely available.

Comparative Analysis of Human-Chimpanzee Divergence in Brain Connectivity and its Genetic Correlates

Chimpanzees (Pan troglodytes) are humans' closest living relatives, making them the most directly relevant comparison point for understanding human brain evolution. Zeroing in on the differences in brain connectivity between humans and chimpanzees can provide key insights into the specific evolutionary changes that might have occured along the human lineage. However, conducting comparisons of brain connectivity between humans and chimpanzees remains challenging, as cross-species brain atlases established within the same framework are currently lacking. Without the availability of cross-species brain atlases, the region-wise connectivity patterns between humans and chimpanzees cannot be directly compared. To address this gap, we built the first Chimpanzee Brainnetome Atlas (ChimpBNA) by following a well-established connectivity-based parcellation framework. Leveraging this new resource, we found substantial divergence in connectivity patterns across most association cortices, notably in the lateral temporal and dorsolateral prefrontal cortex between the two species. Intriguingly, these patterns significantly deviate from the patterns of cortical expansion observed in humans compared to chimpanzees. Additionally, we identified regions displaying connectional asymmetries that differed between species, likely resulting from evolutionary divergence. Genes associated with these divergent connectivities were found to be enriched in cell types crucial for cortical projection circuits and synapse formation. These genes exhibited more pronounced differences in expression patterns in regions with higher connectivity divergence, suggesting a potential foundation for brain connectivity evolution. Therefore, our study not only provides a fine-scale brain atlas of chimpanzees but also highlights the connectivity divergence between humans and chimpanzees in a more rigorous and comparative manner and suggests potential genetic correlates for the observed divergence in brain connectivity patterns between the two species. This can help us better understand the origins and development of uniquely human cognitive capabilities.

Automated individual cortical parcellation via consensus graph representation learning

Cortical parcellation plays a pivotal role in elucidating the brain organization. Despite the growing efforts to develop parcellation algorithms using functional magnetic resonance imaging, achieving a balance between intra-individual specificity and inter-individual consistency proves challenging, making the generation of high-quality, subject-consistent cortical parcellations particularly elusive. To solve this problem, our paper proposes a fully automated individual cortical parcellation method based on consensus graph representation learning. The method integrates spectral embedding with low-rank tensor learning into a unified optimization model, which uses group-common connectivity patterns captured by low-rank tensor learning to optimize subjects' functional networks. This not only ensures consistency in brain representations across different subjects but also enhances the quality of each subject's representation matrix by eliminating spurious connections. More importantly, it achieves an adaptive balance between intra-individual specificity and inter-individual consistency during this process. Experiments conducted on a test-retest dataset from the Human Connectome Project (HCP) demonstrate that our method outperforms existing methods in terms of reproducibility, functional homogeneity, and alignment with task activation. Extensive network-based comparisons on the HCP S900 dataset reveal that the functional network derived from our cortical parcellation method exhibits greater capabilities in gender identification and behavior prediction than other approaches.

How the intrinsic functional connectivity patterns of the semantic network support semantic processing

Semantic processing, a core of language comprehension, involves the activation of brain regions dispersed extensively across the frontal, temporal, and parietal cortices that compose the semantic network. To comprehend the functional structure of this semantic network and how it prepares for semantic processing, we investigated its intrinsic functional connectivity (FC) and the relation between this pattern and semantic processing ability in a large sample from the Human Connectome Project (HCP) dataset. We first defined a well-studied brain network for semantic processing, and then we characterized the within-network connectivity (WNC) and the between-network connectivity (BNC) within this network using a voxel-based global brain connectivity (GBC) method based on resting-state functional magnetic resonance imaging (fMRI). The results showed that 97.73% of the voxels in the semantic network displayed considerably greater WNC than BNC, demonstrating that the semantic network is a fairly encapsulated network. Moreover, multiple connector hubs in the semantic network were identified after applying the criterion of WNC > 1 SD above the mean WNC of the semantic network. More importantly, three of these connector hubs (i.e., the left anterior temporal lobe, angular gyrus, and orbital part of the inferior frontal gyrus) were reliably associated with semantic processing ability. Our findings suggest that the three identified regions use WNC as the central mechanism for supporting semantic processing and that task-independent spontaneous connectivity in the semantic network is essential for semantic processing.

Principles of cortical areas and their implications for neuroimaging

Cortical organization should constrain the study of how the brain performs behavior and cognition. A fundamental concept in cortical organization is that of arealization: that the cortex is parceled into discrete areas. In part one of this report, we review how non-human animal studies have illuminated principles of cortical arealization by revealing: (1) what defines a cortical area, (2) how cortical areas are formed, (3) how cortical areas interact with one another, and (4) what "computations" or "functions" areas perform. In part two, we discuss how these principles apply to neuroimaging research. In doing so, we highlight several examples where the commonly accepted interpretation of neuroimaging observations requires assumptions that violate the principles of arealization, including nonstationary areas that move on short time scales, large-scale gradients as organizing features, and cortical areas with singular functionality that perfectly map psychological constructs. Our belief is that principles of neurobiology should strongly guide the nature of computational explanations.

Individual connectivity-based parcellations reflect functional properties of human auditory cortex

Neuroimaging studies of the functional organization of human auditory cortex have focused on group-level analyses to identify tendencies that represent the typical brain. Here, we mapped auditory areas of the human superior temporal cortex (STC) in 30 participants by combining functional network analysis and 1-mm isotropic resolution 7T functional magnetic resonance imaging (fMRI). Two resting-state fMRI sessions, and one or two auditory and audiovisual speech localizer sessions, were collected on 3-4 separate days. We generated a set of functional network-based parcellations from these data. Solutions with 4, 6, and 11 networks were selected for closer examination based on local maxima of Dice and Silhouette values. The resulting parcellation of auditory cortices showed high intraindividual reproducibility both between resting state sessions (Dice coefficient: 69-78%) and between resting state and task sessions (Dice coefficient: 62-73%). This demonstrates that auditory areas in STC can be reliably segmented into functional subareas. The interindividual variability was significantly larger than intraindividual variability (Dice coefficient: 57%-68%, p<0.001), indicating that the parcellations also captured meaningful interindividual variability. The individual-specific parcellations yielded the highest alignment with task response topographies, suggesting that individual variability in parcellations reflects individual variability in auditory function. Connectional homogeneity within networks was also highest for the individual-specific parcellations. Furthermore, the similarity in the functional parcellations was not explainable by the similarity of macroanatomical properties of auditory cortex. Our findings suggest that individual-level parcellations capture meaningful idiosyncrasies in auditory cortex organization.

Diurnal Variation of Brain Activity in the Human Suprachiasmatic Nucleus

The suprachiasmatic nucleus (SCN) is the central clock for circadian rhythms. Animal studies have revealed daily rhythms in the neuronal activity in the SCN. However, the circadian activity of the human SCN has remained elusive. In this study, to reveal the diurnal variation of the SCN activity in humans, we localized the SCN by employing an areal boundary mapping technique to resting-state functional images and investigated the SCN activity using perfusion imaging. In the first experiment (n = 27, including both sexes), we scanned each participant four times a day, every 6 h. Higher activity was observed at noon, while lower activity was recorded in the early morning. In the second experiment (n = 20, including both sexes), the SCN activity was measured every 30 min for 6 h from midnight to dawn. The results showed that the SCN activity gradually decreased and was not associated with the electroencephalography. Furthermore, the SCN activity was compatible with the rodent SCN activity after switching off the lights. These results suggest that the diurnal variation of the human SCN follows the zeitgeber cycles of nocturnal and diurnal mammals and is modulated by physical lights rather than the local time.

Intrinsic functional clustering of the macaque insular cortex

The functional organization of the primate insula has been studied using a variety of techniques focussing on regional differences in either architecture, connectivity, or function. These complementary methods offered insights into the complex organization of the insula and proposed distinct parcellation schemes at varying levels of detail and complexity. The advent of imaging techniques that allow non-invasive assessment of structural and functional connectivity, has popularized data-driven connectivity-based parcellation methods to investigate the organization of the human insula. Yet, it remains unclear if the subdivisions derived from these data-driven clustering methods reflect meaningful descriptions of the functional specialization of the insula. In this study, we employed hierarchical clustering to examine the cluster parcellations of the macaque insula. As our aim was exploratory, we examined parcellations consisting of two up to ten clusters. Three different cluster validation methods (fingerprinting, silhouette, elbow) converged on a four-cluster solution as the most optimal representation of our data. Examining functional response properties of these clusters, in addition to their brain-wide functional connectivity suggested a functional specialization related to processing gustatory, somato-motor, vestibular and social visual cues. However, a more detailed functional differentiation aligning with previous functional investigations of insula subfields became evident at higher cluster numbers beyond the proposed optimal four clusters. Overall, our findings demonstrate that resting-state-based hierarchical clustering can provide a meaningful description of the insula's functional organization at some level of detail. Nonetheless, cluster parcellations derived from this method are best combined with data obtained through other modalities, to provide a more comprehensive and detailed account of the insula's complex functional organization.

Can hubs of the human connectome be identified consistently with diffusion MRI?

Recent years have seen a surge in the use of diffusion MRI to map connectomes in humans, paralleled by a similar increase in processing and analysis choices. Yet these different steps and their effects are rarely compared systematically. Here, in a healthy young adult population (n = 294), we characterized the impact of a range of analysis pipelines on one widely studied property of the human connectome: its degree distribution. We evaluated the effects of 40 pipelines (comparing common choices of parcellation, streamline seeding, tractography algorithm, and streamline propagation constraint) and 44 group-representative connectome reconstruction schemes on highly connected hub regions. We found that hub location is highly variable between pipelines. The choice of parcellation has a major influence on hub architecture, and hub connectivity is highly correlated with regional surface area in most of the assessed pipelines (ρ > 0.70 in 69% of the pipelines), particularly when using weighted networks. Overall, our results demonstrate the need for prudent decision-making when processing diffusion MRI data, and for carefully considering how different processing choices can influence connectome organization.

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.

Inter-regional correlation estimators for functional magnetic resonance imaging

Functional magnetic resonance imaging (fMRI) functional connectivity between brain regions is often computed using parcellations defined by functional or structural atlases. Typically, some kind of voxel averaging is performed to obtain a single temporal correlation estimate per region pair. However, several estimators can be defined for this task, with various assumptions and degrees of robustness to local noise, global noise, and region size. In this paper, we systematically present and study the properties of 9 different functional connectivity estimators taking into account the spatial structure of fMRI data, based on a simple fMRI data spatial model. These include 3 existing estimators and 6 novel estimators. We demonstrate the empirical properties of the estimators using synthetic, animal, and human data, in terms of graph structure, repeatability and reproducibility, discriminability, dependence on region size, as well as local and global noise robustness. We prove analytically the link between regional intra-correlation and inter-region correlation, and show that the choice of estimator has a strong influence on inter-correlation values. Some estimators, including the commonly used correlation of averages (ca), are positively biased, and have more dependence to region size and intra-correlation than robust alternatives, resulting in spatially-dependent bias. We define the new local correlation of averages estimator with better theoretical guarantees, lower bias, significantly lower dependence on region size (Spearman correlation 0.40 vs 0.55, paired t-test T=27.2, p=1.1e-47), at negligible cost to discriminative power, compared to the ca estimator. The difference in connectivity pattern between the estimators is not distributed uniformly throughout the brain, but rather shows a clear ventral-dorsal gradient, suggesting that region size and intra-correlation plays an important role in shaping functional networks defined using the ca estimator, and leading to non-trivial differences in their connectivity structure. We provide an open source R package and equivalent Python implementation to facilitate the use of the new estimators, together with preprocessed rat time-series.

A parcellation scheme of mouse isocortex based on reversals in connectivity gradients

The brain is composed of several anatomically clearly separated structures. This parcellation is often extended into the isocortex, based on anatomical, physiological, or functional differences. Here, we derive a parcellation scheme based purely on the spatial structure of long-range synaptic connections within the cortex. To that end, we analyzed a publicly available dataset of average mouse brain connectivity, and split the isocortex into disjunct regions. Instead of clustering connectivity based on modularity, our scheme is inspired by methods that split sensory cortices into subregions where gradients of neuronal response properties, such as the location of the receptive field, reverse. We calculated comparable gradients from voxelized brain connectivity data and automatically detected reversals in them. This approach better respects the known presence of functional gradients within brain regions than clustering-based approaches. Placing borders at the reversals resulted in a parcellation into 41 subregions that differs significantly from an established scheme in nonrandom ways, but is comparable in terms of the modularity of connectivity between regions. It reveals unexpected trends of connectivity, such as a tripartite split of somatomotor regions along an anterior to posterior gradient. The method can be readily adapted to other organisms and data sources, such as human functional connectivity.

PPA: Principal parcellation analysis for brain connectomes and multiple traits

Our understanding of the structure of the brain and its relationships with human traits is largely determined by how we represent the structural connectome. Standard practice divides the brain into regions of interest (ROIs) and represents the connectome as an adjacency matrix having cells measuring connectivity between pairs of ROIs. Statistical analyses are then heavily driven by the (largely arbitrary) choice of ROIs. In this article, we propose a human trait prediction framework utilizing a tractography-based representation of the brain connectome, which clusters fiber endpoints to define a data-driven white matter parcellation targeted to explain variation among individuals and predict human traits. This leads to Principal Parcellation Analysis (PPA), representing individual brain connectomes by compositional vectors building on a basis system of fiber bundles that captures the connectivity at the population level. PPA eliminates the need to choose atlases and ROIs a priori, and provides a simpler, vector-valued representation that facilitates easier statistical analysis compared to the complex graph structures encountered in classical connectome analyses. We illustrate the proposed approach through applications to data from the Human Connectome Project (HCP) and show that PPA connectomes improve power in predicting human traits over state-of-the-art methods based on classical connectomes, while dramatically improving parsimony and maintaining interpretability. Our PPA package is publicly available on GitHub, and can be implemented routinely for diffusion image data.

Accurate Bayesian segmentation of thalamic nuclei using diffusion MRI and an improved histological atlas

The human thalamus is a highly connected brain structure, which is key for the control of numerous functions and is involved in several neurological disorders. Recently, neuroimaging studies have increasingly focused on the volume and connectivity of the specific nuclei comprising this structure, rather than looking at the thalamus as a whole. However, accurate identification of cytoarchitectonically designed histological nuclei on standard in vivo structural MRI is hampered by the lack of image contrast that can be used to distinguish nuclei from each other and from surrounding white matter tracts. While diffusion MRI may offer such contrast, it has lower resolution and lacks some boundaries visible in structural imaging. In this work, we present a Bayesian segmentation algorithm for the thalamus. This algorithm combines prior information from a probabilistic atlas with likelihood models for both structural and diffusion MRI, allowing segmentation of 25 thalamic labels per hemisphere informed by both modalities. We present an improved probabilistic atlas, incorporating thalamic nuclei identified from histology and 45 white matter tracts surrounding the thalamus identified in ultra-high gradient strength diffusion imaging. We present a family of likelihood models for diffusion tensor imaging, ensuring compatibility with the vast majority of neuroimaging datasets that include diffusion MRI data. The use of these diffusion likelihood models greatly improves identification of nuclear groups versus segmentation based solely on structural MRI. Dice comparison of 5 manually identifiable groups of nuclei to ground truth segmentations show improvements of up to 10 percentage points. Additionally, our chosen model shows a high degree of reliability, with median test-retest Dice scores above 0.85 for four out of five nuclei groups, whilst also offering improved detection of differential thalamic involvement in Alzheimer's disease (AUROC 81.98%). The probabilistic atlas and segmentation tool will be made publicly available as part of the neuroimaging package FreeSurfer (https://freesurfer.net/fswiki/ThalamicNucleiDTI).

Hippocampal anterior- posterior shift in childhood and adolescence

Hippocampal-cortical networks play an important role in neurocognitive development. Applying the method of Connectivity-Based Parcellation (CBP) on hippocampal-cortical structural covariance (SC) networks computed from T1-weighted magnetic resonance images, we examined how the hippocampus differentiates into subregions during childhood and adolescence (N = 1105, 6-18 years). In late childhood, the hippocampus mainly differentiated along the anterior-posterior axis similar to previous reported functional differentiation patterns of the hippocampus. In contrast, in adolescence a differentiation along the medial-lateral axis was evident, reminiscent of the cytoarchitectonic division into cornu ammonis and subiculum. Further meta-analytical characterization of hippocampal subregions in terms of related structural co-maturation networks, behavioural and gene profiling suggested that the hippocampal head is related to higher order functions (e.g. language, theory of mind, autobiographical memory) in late childhood morphologically co-varying with almost the whole brain. In early adolescence but not in childhood, posterior subicular SC networks were associated with action-oriented and reward systems. The findings point to late childhood as an important developmental period for hippocampal head morphology and to early adolescence as a crucial period for hippocampal integration into action- and reward-oriented cognition. The latter may constitute a developmental feature that conveys increased propensity for addictive disorders.

Homotopic local-global parcellation of the human cerebral cortex from resting-state functional connectivity

Resting-state fMRI is commonly used to derive brain parcellations, which are widely used for dimensionality reduction and interpreting human neuroscience studies. We previously developed a model that integrates local and global approaches for estimating areal-level cortical parcellations. The resulting local-global parcellations are often referred to as the Schaefer parcellations. However, the lack of homotopic correspondence between left and right Schaefer parcels has limited their use for brain lateralization studies. Here, we extend our previous model to derive homotopic areal-level parcellations. Using resting-fMRI and task-fMRI across diverse scanners, acquisition protocols, preprocessing and demographics, we show that the resulting homotopic parcellations are as homogeneous as the Schaefer parcellations, while being more homogeneous than five publicly available parcellations. Furthermore, weaker correlations between homotopic parcels are associated with greater lateralization in resting network organization, as well as lateralization in language and motor task activation. Finally, the homotopic parcellations agree with the boundaries of a number of cortical areas estimated from histology and visuotopic fMRI, while capturing sub-areal (e.g., somatotopic and visuotopic) features. Overall, these results suggest that the homotopic local-global parcellations represent neurobiologically meaningful subdivisions of the human cerebral cortex and will be a useful resource for future studies. Multi-resolution parcellations estimated from 1479 participants are publicly available (https://github.com/ThomasYeoLab/CBIG/tree/master/stable_projects/brain_parcellation/Yan2023_homotopic).

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.

Hypothalamic interaction with reward-related regions during subjective evaluation of foods

The reward system implemented in the midbrain, ventral striatum, orbitofrontal cortex, and ventromedial prefrontal cortex evaluates and compares various types of rewards given to the organisms. It has been suggested that autonomic factors influence reward-related processing via the hypothalamus, but how the hypothalamus modulates the reward system remains elusive. In this functional magnetic resonance imaging study, the hypothalamus was parcellated into individual hypothalamic nuclei performing different autonomic functions using boundary mapping parcellation analyses. The effective interaction during subjective evaluation of foods in a reward task was then investigated between the human hypothalamic nuclei and the reward-related regions. We found significant brain activity decrease in the paraventricular nucleus (PVH) and lateral nucleus in the hypothalamus in food evaluation compared with monetary evaluation. A psychophysiological interaction analysis revealed dual interactions between the PVH and (1) midbrain region and (2) ventromedial prefrontal cortex, with the former correlated with the stronger tendency of participants toward food-seeking. A dynamic causal modeling analysis further revealed unidirectional interactions from the PVH to the midbrain and ventromedial prefrontal cortex. These results suggest that the PVH in the human hypothalamus interacts with the reward-related regions in the cerebral cortex via multiple pathways (i.e., the midbrain pathway and ventromedial prefrontal pathway) to evaluate rewards for subsequent decision-making.

Differences in resting state functional connectivity underlie visuomotor performance declines in older adults with a genetic risk (APOE ε4) for Alzheimer’s disease

Introduction

Non-standard visuomotor integration requires the interaction of large networks in the brain. Previous findings have shown that non-standard visuomotor performance is impaired in individuals with specific dementia risk factors (family history of dementia and presence of the APOE ε4 allele) in advance of any cognitive impairments. These findings suggest that visuomotor impairments are associated with early dementia-related brain changes. The current study examined the underlying resting state functional connectivity (RSFC) associated with impaired non-standard visuomotor performance, as well as the impacts of dementia family history, sex, and APOE status.

Methods

Cognitively healthy older adults (n = 48) were tested on four visuomotor tasks where reach and gaze were increasingly spatially dissociated. Participants who had a family history of dementia or the APOE ε4 allele were considered to be at an increased risk for AD. To quantify RSFC within networks of interest, an EPI sequence sensitive to BOLD contrast was collected. The networks of interest were the default mode network (DMN), somatomotor network (SMN), dorsal attention network (DAN), ventral attention network (VAN), and frontoparietal control network (FPN).

Results

Individuals with the ε4 allele showed abnormalities in RSFC between posterior DMN nodes that predicted poorer non-standard visuomotor performance. Specifically, multiple linear regression analyses revealed lower RSFC between the precuneus/posterior cingulate cortex and the left inferior parietal lobule as well as the left parahippocampal cortex. Presence of the APOE ε4 allele also modified the relationship between mean DAN RSFC and visuomotor control, where lower mean RSFC in the DAN predicted worse non-standard visuomotor performance only in APOE ε4 carriers. There were otherwise no effects of family history, APOE ε4 status, or sex on the relationship between RSFC and visuomotor performance for any of the other resting networks.

Conclusion

The preliminary findings provide insight into the impact of APOE ε4-related genetic risk on neural networks underlying complex visuomotor transformations, and demonstrate that the non-standard visuomotor task paradigm discussed in this study may be used as a non-invasive, easily accessible assessment tool for dementia risk.

High-resolution cortical parcellation based on conserved brain landmarks for localization of multimodal data to the nearest centimeter

Precise cortical brain localization presents an important challenge in the literature. Brain atlases provide data-guided parcellation based on functional and structural brain metrics, and each atlas has its own unique benefits for localization. We offer a parcellation guided by intracranial electroencephalography, a technique which has historically provided pioneering advances in our understanding of brain structure-function relationships. We used a consensus boundary mapping approach combining anatomical designations in Duvernoy's Atlas of the Human Brain, a widely recognized textbook of human brain anatomy, with the anatomy of the MNI152 template and the magnetic resonance imaging scans of an epilepsy surgery cohort. The Yale Brain Atlas consists of 690 one-square centimeter parcels based around conserved anatomical features and each with a unique identifier to communicate anatomically unambiguous localization. We report on the methodology we used to create the Atlas along with the findings of a neuroimaging study assessing the accuracy and clinical usefulness of cortical localization using the Atlas. We also share our vision for the Atlas as a tool in the clinical and research neurosciences, where it may facilitate precise localization of data on the cortex, accurate description of anatomical locations, and modern data science approaches using standardized brain regions.

Disorganized functional architecture of amygdala subregional networks in obsessive-compulsive disorder

A precise understanding of amygdala-centered subtle networks may help refine neurocircuitry models of obsessive-compulsive disorder (OCD). We applied connectivity-based parcellation methodology to segment the amygdala based on resting-state fMRI data of 92 medication-free OCD patients without comorbidity and 90 matched healthy controls (HC). The amygdala was parcellated into two subregions corresponding to basolateral amygdala (BLA) and centromedial amygdala (CMA). Amygdala subregional functional connectivity (FC) maps were generated and group differences were evaluated with diagnosis-by-subregion flexible factorial ANOVA. We found significant diagnosis × subregion FC interactions in insula, supplementary motor area (SMA), midcingulate cortex (MCC), superior temporal gyrus (STG) and postcentral gyrus (PCG). In HC, the BLA demonstrated stronger connectivity with above regions compared to CMA, whereas in OCD, the connectivity pattern reversed to stronger CMA connectivity comparing to BLA. Relative to HC, OCD patients exhibited hypoconnectivity between left BLA and left insula, and hyperconnectivity between right CMA and SMA, MCC, insula, STG, and PCG. Moreover, OCD patients showed reduced volume of left BLA and right CMA compared to HC. Our findings characterized disorganized functional architecture of amygdala subregional networks in accordance with structural defects, providing direct evidence regarding the specific role of amygdala subregions in the neurocircuitry models of OCD.

Associations between early-life stress exposure and internalizing symptomatology during the COVID-19 pandemic: Assessing the role of neurobehavioral mediators

Background

The ongoing COVID-19 pandemic is a major stressor that has been associated with increased risk for psychiatric illness in the general population. Recent work has highlighted that experiences of early-life stress (ELS) may impact individuals' psychological functioning and vulnerability for developing internalizing psychopathology in response to pandemic-related stress. However, little is known about the neurobehavioral factors that may mediate the association between ELS exposure and COVID-related internalizing symptomatology. The current study sought to examine the mediating roles of pre-pandemic resting-state frontoamygdala connectivity and concurrent emotion regulation (ER) in the association between ELS and pandemic-related internalizing symptomatology.

Methods

Retrospective life-stress histories, concurrent self-reported ER strategies (i.e., reappraisal and suppression), concurrent self-reported internalizing symptomatology (i.e., depression- and anxiety-related symptomatology), and resting-state functional connectivity data from a sample of adults (N = 64, M age = 22.12, female = 68.75%) were utilized.

Results

There were no significant direct associations between ELS and COVID-related internalizing symptomatology. Neither frontoamygdala functional connectivity nor ER strategy use mediated an association between ELS and COVID-related internalizing symptomatology (ps > 0.05). Exploratory analyses identified a significant moderating effect of reappraisal use on the association between ELS and internalizing symptomatology (β = -0.818, p = 0.047), such that increased reappraisal use buffered the impact of ELS on psychopathology.

Conclusions

While frontoamygdala connectivity and ER do not appear to mediate the association between ELS and COVID-related internalizing symptomatology, our findings suggest that the use of reappraisal may buffer against the effect of ELS on mental health during the pandemic.

Connectivity-based thalamus parcellation and surgical targeting of somatosensory subnuclei

OBJECTIVE

The anatomy of the posterolateral thalamus varies substantially between individuals, presenting a challenge for surgical targeting. Patient-specific, connectivity-based parcellation of the thalamus may effectively approximate the ventrocaudal nucleus (Vc). This remains to be robustly validated or assessed as a method to guide surgical targeting. The authors assessed the validity of connectivity-based parcellation for targeting the Vc and its potential for improving clinical outcomes of pain surgery.

METHODS

A cohort of 19 patients with regional, chronic neuropathic pain underwent preoperative structural and diffusion MRI, then progressed to deep brain stimulation targeting the Vc based on traditional atlas coordinates. Surgical thalami were retrospectively segmented and then parcellated based on tractography estimates of thalamocortical connectivity. The location of each patient's electrode array was analyzed with respect to their primary somatosensory cortex (S1) parcel and compared across patients with reference to the thalamic homunculus.

RESULTS

Ten patients achieved long-term pain relief. Sixty-one percent of an average array (interquartile range 42%-74%) was located in the S1 parcel. In patients who achieved long-term benefit from surgery, array location in the individually generated S1 parcels was medial for face pain, centromedial for arm pain, and centrolateral for leg pain. Patients who did not benefit from surgery did not follow this pattern. Standard stereotactic coordinates of electrode locations diverged from this pattern.

CONCLUSIONS

Connectivity-based parcellation of the thalamus appears to be a reliable method for segmenting the Vc. Identifying the Vc in this way, and targeting mediolaterally as appropriate for the region of pain, merits exploration in an effort to increase the yield of successful surgical procedures.

Tasks activating the default mode network map multiple functional systems

Recent developments in network neuroscience suggest reconsidering what we thought we knew about the default mode network (DMN). Although this network has always been seen as unitary and associated with the resting state, a new deconstructive line of research is pointing out that the DMN could be divided into multiple subsystems supporting different functions. By now, it is well known that the DMN is not only deactivated by tasks, but also involved in affective, mnestic, and social paradigms, among others. Nonetheless, it is starting to become clear that the array of activities in which it is involved, might also be extended to more extrinsic functions. The present meta-analytic study is meant to push this boundary a bit further. The BrainMap database was searched for all experimental paradigms activating the DMN, and their activation likelihood estimation maps were then computed. An additional map of task-induced deactivations was also created. A multidimensional scaling indicated that such maps could be arranged along an anatomo-psychological gradient, which goes from midline core activations, associated with the most internal functions, to that of lateral cortices, involved in more external tasks. Further multivariate investigations suggested that such extrinsic mode is especially related to reward, semantic, and emotional functions. However, an important finding was that the various activation maps were often different from the canonical representation of the resting-state DMN, sometimes overlapping with it only in some peripheral nodes, and including external regions such as the insula. Altogether, our findings suggest that the intrinsic-extrinsic opposition may be better understood in the form of a continuous scale, rather than a dichotomy.

Cohesive parcellation of the human brain using resting-state fMRI

BACKGROUND

The data burden for resting-state fMRI analysis rises with increasing resolutions available at ultrahigh fields. Therefore, a fundamental preprocessing step in brain network analysis is to reduce the data, usually by performing some kind of data parcellation. Most functional parcellations based on rsfMRI connectivity are synthesized from the dense connectome. In contrast, most network analyses begin by reducing each parcel to a single exemplar time series. This disconnect between parcel formation and usage assumes that parcel exemplars adequately represent their member voxels, which is not always the case for commonly used parcellations.

Social cognitive network neuroscience

Over the past three decades, research from the field of social neuroscience has identified a constellation of brain regions that relate to social cognition. Although these studies have provided important insights into the specific neural regions underlying social behavior, they may overlook the broader neural context in which those regions and the interactions between them are embedded. Network neuroscience is an emerging discipline that focuses on modeling and analyzing brain networks-collections of interacting neural elements. Because human cognition requires integrating information across multiple brain regions and systems, we argue that a novel social cognitive network neuroscience approach-which leverages methods from the field of network neuroscience and graph theory-can advance our understanding of how brain systems give rise to social behavior. This review provides an overview of the field of network neuroscience, discusses studies that have leveraged this approach to advance social neuroscience research, highlights the potential contributions of social cognitive network neuroscience to understanding social behavior and provides suggested tools and resources for conducting network neuroscience research.

Imaging Genetic Based Mediation Analysis for Human Cognition

The brain connectome maps the structural and functional connectivity that forms an important neurobiological basis for the analysis of human cognitive traits while the genetic predisposition and our cognition ability are frequently found in close association. The issue of how genetic architecture and brain connectome causally affect human behaviors remains unknown. To seek for the potential causal relationship, in this paper, we carried out the causal pathway analysis from single nucleotide polymorphism (SNP) data to four common human cognitive traits, mediated by the brain connectome. Specifically, we selected 942 SNPs that are significantly associated with the brain connectome, and then estimated the direct and indirect effect on the human traits for each SNP. We found out that a majority of the selected SNPs have significant direct effects on human traits and discussed the trait-related brain regions and their implications.

Evaluating brain parcellations using the distance-controlled boundary coefficient

One important approach to human brain mapping is to define a set of distinct regions that can be linked to unique functions. Numerous brain parcellations have been proposed, using cytoarchitectonic, structural, or functional magnetic resonance imaging (fMRI) data. The intrinsic smoothness of brain data, however, poses a problem for current methods seeking to compare different parcellations. For example, criteria that simply compare within‐parcel to between‐parcel similarity provide even random parcellations with a high value. Furthermore, the evaluation is biased by the spatial scale of the parcellation. To address this problem, we propose the distance‐controlled boundary coefficient (DCBC), an unbiased criterion to evaluate discrete parcellations. We employ this new criterion to evaluate existing parcellations of the human neocortex in their power to predict functional boundaries for an fMRI data set with many different tasks, as well as for resting‐state data. We find that common anatomical parcellations do not perform better than chance, suggesting that task‐based functional boundaries do not align well with sulcal landmarks. Parcellations based on resting‐state fMRI data perform well; in some cases, as well as a parcellation defined on the evaluation data itself. Finally, multi‐modal parcellations that combine functional and anatomical criteria perform substantially worse than those based on functional data alone, indicating that functionally homogeneous regions often span major anatomical landmarks. Overall, the DCBC advances the field of functional brain mapping by providing an unbiased metric that compares the predictive ability of different brain parcellations to define brain regions that are functionally maximally distinct.

Sex Modulates the Pathological Aging Effect on Caudate Functional Connectivity in Mild Cognitive Impairment

Purpose

To assess the pathological aging effect on caudate functional connectivity among mild cognitive impairment (MCI) participants and examine whether and how sex and amyloid contribute to this process.

Materials and Methods

Two hundred and seventy-seven functional magnetic resonance imaging (fMRI) sessions from 163 cognitive normal (CN) older adults and 309 sessions from 139 participants with MCI were included as the main sample in our analysis. Pearson's correlation was used to characterize the functional connectivity (FC) between caudate nuclei and each brain region, then caudate nodal strength was computed to quantify the overall caudate FC strength. Association analysis between caudate nodal strength and age was carried out in MCI and CN separately using linear mixed effect (LME) model with covariates (education, handedness, sex, Apolipoprotein E4, and intra-subject effect). Analysis of covariance was conducted to investigate sex, amyloid status, and their interaction effects on aging with the fMRI data subset having amyloid status available. LME model was applied to women and men separately within MCI group to evaluate aging effects on caudate nodal strength and each region's connectivity with caudate nuclei. We then evaluated the roles of sex and amyloid status in the associations of neuropsychological scores with age or caudate nodal strength. An independent cohort was used to validate the sex-dependent aging effects in MCI.

Results

The MCI group had significantly stronger age-related increase of caudate nodal strength compared to the CN group. Analyzing women and men separately revealed that the aging effect on caudate nodal strength among MCI participants was significant only for women (left: P = 6.23 × 10-7, right: P = 3.37 × 10-8), but not for men (P > 0.3 for bilateral caudate nuclei). The aging effects on caudate nodal strength were not significantly mediated by brain amyloid burden. Caudate connectivity with ventral prefrontal cortex substantially contributed to the aging effect on caudate nodal strength in women with MCI. Higher caudate nodal strength is significantly related to worse cognitive performance in women but not in men with MCI.

Conclusion

Sex modulates the pathological aging effects on caudate nodal strength in MCI regardless of amyloid status. Caudate nodal strength may be a sensitive biomarker of pathological aging in women with MCI.

Meta-analytical transdiagnostic neural correlates in common pediatric psychiatric disorders

In the last decades, neuroimaging studies have attempted to unveil the neurobiological markers underlying pediatric psychiatric disorders. Yet, the vast majority of neuroimaging studies still focus on a single nosological category, which limit our understanding of the shared/specific neural correlates between these disorders. Therefore, we aimed to investigate the transdiagnostic neural correlates through a novel and data-driven meta-analytical method. A data-driven meta-analysis was carried out which grouped similar experiments’ topographic map together, irrespectively of nosological categories and task-characteristics. Then, activation likelihood estimation meta-analysis was performed on each group of experiments to extract spatially convergent brain regions. One hundred forty-seven experiments were retrieved (3124 cases compared to 3100 controls): 79 attention-deficit/hyperactivity disorder, 32 conduct/oppositional defiant disorder, 14 anxiety disorders, 22 major depressive disorders. Four significant groups of experiments were observed. Functional characterization suggested that these groups of aberrant brain regions may be implicated internally/externally directed processes, attentional control of affect, somato-motor and visual processes. Furthermore, despite that some differences in rates of studies involving major depressive disorders were noticed, nosological categories were evenly distributed between these four sets of regions. Our results may reflect transdiagnostic neural correlates of pediatric psychiatric disorders, but also underscore the importance of studying pediatric psychiatric disorders simultaneously rather than independently to examine differences between disorders.

A novel method for estimating connectivity-based parcellation of the human brain from diffusion MRI: Application to an aging cohort

Connectivity‐based parcellation (CBP) methods are used to define homogenous and biologically meaningful parcels or nodes—the foundations of brain network fingerprinting—by grouping voxels with similar patterns of brain connectivity. However, we still lack a gold standard method and the use of CBPs to study the aging brain remains scarce. Our study proposes a novel CBP method from diffusion MRI data and shows its potential to produce a more accurate characterization of the longitudinal alterations in brain network topology occurring in aging. For this, we constructed whole‐brain connectivity maps from diffusion MRI data of two datasets: an aging cohort evaluated at two timepoints (mean interval time: 52.8 ± 7.24 months) and a normative adult cohort—MGH‐HCP. State‐of‐the‐art clustering techniques were used to identify the best performing technique. Furthermore, we developed a new metric (connectivity homogeneity fingerprint [CHF]) to evaluate the success of the final CBP in improving regional/global structural connectivity homogeneity. Our results show that our method successfully generates highly homogeneous parcels, as described by the significantly larger CHF score of the resulting parcellation, when compared to the original. Additionally, we demonstrated that the developed parcellation provides a robust anatomical framework to assess longitudinal changes in the aging brain. Our results reveal that aging is characterized by a reorganization of the brain's structural network involving the decrease of intra‐hemispheric, increase of inter‐hemispheric connectivity, and topological rearrangement. Overall, this study proposes a new methodology to perform accurate and robust evaluations of CBP of the human brain.

Uncovering Cortical Units of Processing From Multi-Layered Connectomes

Modern diffusion and functional magnetic resonance imaging (dMRI/fMRI) provide non-invasive high-resolution images from which multi-layered networks of whole-brain structural and functional connectivity can be derived. Unfortunately, the lack of observed correspondence between the connectivity profiles of the two modalities challenges the understanding of the relationship between the functional and structural connectome. Rather than focusing on correspondence at the level of connections we presently investigate correspondence in terms of modular organization according to shared canonical processing units. We use a stochastic block-model (SBM) as a data-driven approach for clustering high-resolution multi-layer whole-brain connectivity networks and use prediction to quantify the extent to which a given clustering accounts for the connectome within a modality. The employed SBM assumes a single underlying parcellation exists across modalities whilst permitting each modality to possess an independent connectivity structure between parcels thereby imposing concurrent functional and structural units but different structural and functional connectivity profiles. We contrast the joint processing units to their modality specific counterparts and find that even though data-driven structural and functional parcellations exhibit substantial differences, attributed to modality specific biases, the joint model is able to achieve a consensus representation that well accounts for both the functional and structural connectome providing improved representations of functional connectivity compared to using functional data alone. This implies that a representation persists in the consensus model that is shared by the individual modalities. We find additional support for this viewpoint when the anatomical correspondence between modalities is removed from the joint modeling. The resultant drop in predictive performance is in general substantial, confirming that the anatomical correspondence of processing units is indeed present between the two modalities. Our findings illustrate how multi-modal integration admits consensus representations well-characterizing each individual modality despite their biases and points to the importance of multi-layered connectomes as providing supplementary information regarding the brain's canonical processing units.

Unrevealing Reliable Cortical Parcellation of Individual Brains Using Resting-State Functional Magnetic Resonance Imaging and Masked Graph Convolutions

Brain parcellation helps to understand the structural and functional organization of the cerebral cortex. Resting-state functional magnetic resonance imaging (fMRI) and connectivity analysis provide useful information to delineate individual brain parcels in vivo. We proposed an individualized cortical parcellation based on graph neural networks (GNN) to learn the reliable functional characteristics of each brain parcel on a large fMRI dataset and to infer the areal probability of each vertex on unseen subjects. A subject-specific confidence mask was implemented in the GNN model to account for the tradeoff between the topographic alignment across subjects and functional homogeneity of brain parcels on individual brains. The individualized brain parcellation achieved better functional homogeneity at rest and during cognitive tasks compared with the group-registered atlas (p-values < 0.05). In addition, highly reliable and replicable parcellation maps were generated on multiple sessions of the same subject (intrasubject similarity = 0.89), while notable variations in the topographic organization were captured across subjects (intersubject similarity = 0.81). Moreover, the intersubject variability of brain parcellation indicated large variations in the association cortices while keeping a stable parcellation on the primary cortex. Such topographic variability was strongly associated with the functional connectivity variability, significantly predicted cognitive behaviors, and generally followed the myelination, cytoarchitecture, and functional organization of the human brain. This study provides new avenues to the precise individualized mapping of the cortical areas through deep learning and shows high potentials in the personalized localization diagnosis and treatment of neurological disorders.

Browsing Multiple Subjects When the Atlas Adaptation Cannot Be Achieved via a Warping Strategy

Brain mapping studies often need to identify brain structures or functional circuits into a set of individual brains. To this end, multiple atlases have been published to represent such structures based on different modalities, subject sets, and techniques. The mainstream approach to exploit these atlases consists in spatially deforming each individual data onto a given atlas using dense deformation fields, which supposes the existence of a continuous mapping between atlases and individuals. However, this continuity is not always verified, and this "iconic" approach has limits. We present in this study an alternative, complementary, "structural" approach, which consists in extracting structures from the individual data, and comparing them without deformation. A "structural atlas" is thus a collection of annotated individual data with a common structure nomenclature. It may be used to characterize structure shape variability across individuals or species, or to train machine learning systems. This study exhibits Anatomist, a powerful structural 3D visualization software dedicated to building, exploring, and editing structural atlases involving a large number of subjects. It has been developed primarily to decipher the cortical folding variability; cortical sulci vary enormously in both size and shape, and some may be missing or have various topologies, which makes iconic approaches inefficient to study them. We, therefore, had to build structural atlases for cortical sulci, and use them to train sulci identification algorithms. Anatomist can display multiple subject data in multiple views, supports all kinds of neuroimaging data, including compound structural object graphs, handles arbitrary coordinate transformation chains between data, and has multiple display features. It is designed as a programming library in both C++ and Python languages, and may be extended or used to build dedicated custom applications. Its generic design makes all the display and structural aspects used to explore the variability of the cortical folding pattern work in other applications, for instance, to browse axonal fiber bundles, deep nuclei, functional activations, or other kinds of cortical parcellations. Multimodal, multi-individual, or inter-species display is supported, and adaptations to large scale screen walls have been developed. These very original features make it a unique viewer for structural atlas browsing.

A novel 5D brain parcellation approach based on spatio-temporal encoding of resting fMRI data from deep residual learning

OBJECTIVE

Brain parcellation is an essential aspect of computational neuroimaging research and deals with segmenting the brain into (possibly overlapping) sub-regions employed to study brain anatomy or function. In the context of functional parcellation, brain organization which is often measured via temporal metrics such as coherence, is highly dynamic. This dynamic aspect is ignored in most research, which typically applies anatomically based, fixed regions for each individual, and can produce misleading results.

METHODS

In this work, we propose a novel spatio-temporal-network (5D) brain parcellation scheme utilizing a deep residual network to predict the probability of each voxel belonging to a brain network at each point in time.

RESULTS

We trained 53 4D brain networks and evaluate the ability of these networks to capture spatial and temporal dynamics as well as to show sensitivity to individual or group-level variation (in our case with age).

CONCLUSION

The proposed system generates informative spatio-temporal networks that vary not only across individuals but also over time and space.

SIGNIFICANCE

The dynamic 5D nature of the developed approach provides a powerful framework that expands on existing work and has potential to identify novel and typically ignored findings when studying the healthy and disordered brain.

Spatial Stability of Functional Networks: A Measure to Assess the Robustness of Graph-Theoretical Metrics to Spatial Errors Related to Brain Parcellation

There is growing interest in studying human brain connectivity and in modelling the brain functional structure as a network. Brain network creation requires parcellation of the cerebral cortex to define nodes. Parcellation might be affected by possible errors due to inter- and intra-subject variability as a consequence of brain structural and physiological characteristics and shape variations related to ageing and diseases, acquisition noise, and misregistration. These errors could induce a knock-on effect on network measure variability. The aim of this study was to investigate spatial stability, a measure of functional connectivity variations induced by parcellation errors. We simulated parcellation variability with random small spatial changes and evaluated its effects on twenty-seven graph-theoretical measures. The study included subjects from three public online datasets. Two brain parcellations were performed using FreeSurfer with geometric atlases. Starting from these, 100 new parcellations were created by increasing the area of 30% of parcels, reducing the area of neighbour parcels, with a rearrangement of vertices. fMRI data were filtered with linear regression, CompCor, and motion correction. Adjacency matrices were constructed with 0.1, 0.2, 0.3, and 0.4 thresholds. Differences in spatial stability between datasets, atlases, and threshold were evaluated. The higher spatial stability resulted for Characteristic-path-length, Density, Transitivity, and Closeness-centrality, and the lower spatial stability resulted for Bonacich and Katz. Multivariate analysis showed a significant effect of atlas, datasets, and thresholds. Katz and Bonacich centrality, which was subject to larger variations, can be considered an unconventional graph measure, poorly implemented in the clinical field and not yet investigated for reliability assessment. Spatial stability (SS) is affected by threshold, and it decreases with increasing threshold for several measures. Moreover, SS seems to depend on atlas choice and scanning parameters. Our study highlights the importance of paying close attention to possible parcellation-related spatial errors, which may affect the reliability of functional connectivity measures.

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 many dimensions of human hippocampal organization and (dys)function

The internal organization of hippocampal formation has been studied for more than a century. Although early accounts emphasized its subfields along the medial-lateral axis, findings in recent decades have highlighted also the anterior-to-posterior (i.e., longitudinal) axis as a key contributor to this brain region's functional organization. Hence, understanding of hippocampal function likely demands characterizing both medial-to-lateral and anterior-to-posterior axes, an approach that has been concretized by recent advances in in vivo parcellation and gradient mapping techniques. Following a short historical overview, we review the evidence provided by these approaches in brain-mapping studies, as well as the perspectives they open for addressing the behavioral relevance of the interacting organizational axes in healthy and clinical populations.

Group-Wise Cortical Surface Parcellation Based on Inter-Subject Fiber Clustering. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2021;2021:2655-2659. [PMID: 34891798 DOI: 10.1109/embc46164.2021.9631099] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

We present an automatic algorithm for the group-wise parcellation of the cortical surface. The method is based on the structural connectivity obtained from representative brain fiber clusters, calculated via an inter-subject clustering scheme. Preliminary regions were defined from cluster-cortical mesh intersection points. The final parcellation was obtained using parcel probability maps to model and integrate the connectivity information of all subjects, and graphs to represent the overlap between parcels. Two inter-subject clustering schemes were tested, generating a total of 171 and 109 parcels, respectively. The resulting parcels were quantitatively compared with three state-of-the-art atlases. The best parcellation returned 69 parcels with a Dice similarity coefficient greater than 0.5. To the best of our knowledge, this is the first diffusion-based cortex parcellation method based on whole-brain inter-subject fiber clustering.

MULTIMODAL HUMAN BRAIN LONGITUDINAL PARCELLATION ACROSS LIFESPAN

Subdividing the human brain into several functionally distinct and spatially contiguous areas is important to understand the amazingly complex human cerebral cortex. However, adult aging is related to differences in the structure, function, and connectivity of brain areas, so that the single population subdivision does not apply to multiple age groups. Moreover, different modalities could provide affirmative and complementary information for the human brain subdivision. To obtain a more reasonable subdivision of the cerebral cortex, we make use of multimodal information to subdivide the human cerebral cortex across lifespan. Specifically, we first construct a population average functional connectivity matrix for each modality of each age group. Second, we separately calculate the population average similarity matrix for the cortical thickness and myelin modality of each age group. Finally, we fuse these population average matrixes to obtain the multimodal similarity matrix and feed it into the spectral clustering algorithm to generate the brain parcellation for each age group.

Connectivity alterations in autism reflect functional idiosyncrasy

Autism spectrum disorder (ASD) is commonly understood as an alteration of brain networks, yet case-control analyses against typically-developing controls (TD) have yielded inconsistent results. Here, we devised a novel approach to profile the inter-individual variability in functional network organization and tested whether such idiosyncrasy contributes to connectivity alterations in ASD. Studying a multi-centric dataset with 157 ASD and 172 TD, we obtained robust evidence for increased idiosyncrasy in ASD relative to TD in default mode, somatomotor and attention networks, but also reduced idiosyncrasy in lateral temporal cortices. Idiosyncrasy increased with age and significantly correlated with symptom severity in ASD. Furthermore, while patterns of functional idiosyncrasy were not correlated with ASD-related cortical thickness alterations, they co-localized with the expression patterns of ASD risk genes. Notably, we could demonstrate that patterns of atypical idiosyncrasy in ASD closely overlapped with connectivity alterations that are measurable with conventional case-control designs and may, thus, be a principal driver of inconsistency in the autism connectomics literature. These findings support important interactions between inter-individual heterogeneity in autism and functional signatures. Our findings provide novel biomarkers to study atypical brain development and may consolidate prior research findings on the variable nature of connectome level anomalies in autism.

Functional network topology in drug resistant and well-controlled idiopathic generalized epilepsy: a resting state functional MRI study

Despite an increasing number of drug treatment options for people with idiopathic generalized epilepsy (IGE), drug resistance remains a significant issue and the mechanisms underlying it remain poorly understood. Previous studies have largely focused on potential cellular or genetic explanations for drug resistance. However, epilepsy is understood to be a network disorder and there is a growing body of literature suggesting altered topology of large-scale resting networks in people with epilepsy compared with controls. We hypothesize that network alterations may also play a role in seizure control. The aim of this study was to compare resting state functional network structure between well-controlled IGE (WC-IGE), drug resistant IGE (DR-IGE) and healthy controls. Thirty-three participants with IGE (10 with WC-IGE and 23 with DR-IGE) and 34 controls were included. Resting state functional MRI networks were constructed using the Functional Connectivity Toolbox (CONN). Global graph theoretic network measures of average node strength (an equivalent measure to mean degree in a network that is fully connected), node strength distribution variance, characteristic path length, average clustering coefficient, small-world index and average betweenness centrality were computed. Graphs were constructed separately for positively weighted connections and for absolute values. Individual nodal values of strength and betweenness centrality were also measured and ‘hub nodes’ were compared between groups. Outcome measures were assessed across the three groups and between both groups with IGE and controls. The IGE group as a whole had a higher average node strength, characteristic path length and average betweenness centrality. There were no clear differences between groups according to seizure control. Outcome metrics were sensitive to whether negatively correlated connections were included in network construction. There were no clear differences in the location of ‘hub nodes’ between groups. The results suggest that, irrespective of seizure control, IGE interictal network topology is more regular and has a higher global connectivity compared to controls, with no alteration in hub node locations. These alterations may produce a resting state network that is more vulnerable to transitioning to the seizure state. It is possible that the lack of apparent influence of seizure control on network topology is limited by challenges in classifying drug response. It is also demonstrated that network topological features are influenced by the sign of connectivity weights and therefore future methodological work is warranted to account for anticorrelations in graph theoretic studies.

Hierarchical correction of p-values via an ultrametric tree running Ornstein-Uhlenbeck process

Statistical testing is classically used as an exploratory tool to search for association between a phenotype and many possible explanatory variables. This approach often leads to multiple testing under dependence. We assume a hierarchical structure between tests via an Ornstein-Uhlenbeck process on a tree. The process correlation structure is used for smoothing the p-values. We design a penalized estimation of the mean of the Ornstein-Uhlenbeck process for p-value computation. The performances of the algorithm are assessed via simulations. Its ability to discover new associations is demonstrated on a metagenomic dataset. The corresponding R package is available from https://github.com/abichat/zazou.

Ventral intermediate nucleus structural connectivity-derived segmentation: anatomical reliability and variability

The Ventral intermediate nucleus (Vim) of thalamus is the most targeted structure for the treatment of drug-refractory tremors. Since methodological differences across existing studies are remarkable and no gold-standard pipeline is available, in this study, we tested different parcellation pipelines for tractography-derived putative Vim identification. Thalamic parcellation was performed on a high quality, multi-shell dataset and a downsampled, clinical-like dataset using two different diffusion signal modeling techniques and two different voxel classification criteria, thus implementing a total of four parcellation pipelines. The most reliable pipeline in terms of inter-subject variability has been picked and parcels putatively corresponding to motor thalamic nuclei have been selected by calculating similarity with a histology-based mask of Vim. Then, spatial relations with optimal stimulation points for the treatment of essential tremor have been quantified. Finally, effect of data quality and parcellation pipelines on a volumetric index of connectivity clusters has been assessed. We found that the pipeline characterized by higher-order signal modeling and threshold-based voxel classification criteria was the most reliable in terms of inter-subject variability regardless data quality. The maps putatively corresponding to Vim were those derived by precentral and dentate nucleus-thalamic connectivity. However, tractography-derived functional targets showed remarkable differences in shape and sizes when compared to a ground truth model based on histochemical staining on seriate sections of human brain. Thalamic voxels connected to contralateral dentate nucleus resulted to be the closest to literature-derived stimulation points for essential tremor but at the same time showing the most remarkable inter-subject variability. Finally, the volume of connectivity parcels resulted to be significantly influenced by data quality and parcellation pipelines. Hence, caution is warranted when performing thalamic connectivity-based segmentation for stereotactic targeting.

Principles and open questions in functional brain network reconstruction

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network representation involves often covert theoretical assumptions and methodological choices which affect the way networks are reconstructed from experimental data, and ultimately the resulting network properties and their interpretation. Here, we review some fundamental conceptual underpinnings and technical issues associated with brain network reconstruction, and discuss how their mutual influence concurs in clarifying the organization of brain function.

In Vivo Super-Resolution Track-Density Imaging for Thalamic Nuclei Identification

The development of novel techniques for the in vivo, non-invasive visualization and identification of thalamic nuclei has represented a major challenge for human neuroimaging research in the last decades. Thalamic nuclei have important implications in various key aspects of brain physiology and many of them show selective alterations in various neurologic and psychiatric disorders. In addition, both surgical stimulation and ablation of specific thalamic nuclei have been proven to be useful for the treatment of different neuropsychiatric diseases. The present work aimed at describing a novel protocol for histologically guided delineation of thalamic nuclei based on short-tracks track-density imaging (stTDI), which is an advanced imaging technique exploiting high angular resolution diffusion tractography to obtain super-resolved white matter maps. We demonstrated that this approach can identify up to 13 distinct thalamic nuclei bilaterally with very high inter-subject (ICC: 0.996, 95% CI: 0.993-0.998) and inter-rater (ICC:0.981; 95% CI:0.963-0.989) reliability, and that both subject-based and group-level thalamic parcellation show a fair share of similarity to a recent standard-space histological thalamic atlas. Finally, we showed that stTDI-derived thalamic maps can be successfully employed to study structural and functional connectivity of the thalamus and may have potential implications both for basic and translational research, as well as for presurgical planning purposes.

Functional parcellation of human and macaque striatum reveals human-specific connectivity in the dorsal caudate

A wide homology between human and macaque striatum is often assumed as in both the striatum is involved in cognition, emotion and executive functions. However, differences in functional and structural organization between human and macaque striatum may reveal evolutionary divergence and shed light on human vulnerability to neuropsychiatric diseases. For instance, dopaminergic dysfunction of the human striatum is considered to be a pathophysiological underpinning of different disorders, such as Parkinson's disease (PD) and schizophrenia (SCZ). Previous investigations have found a wide similarity in structural connectivity of the striatum between human and macaque, leaving the cross-species comparison of its functional organization unknown. In this study, resting-state functional connectivity (RSFC) derived striatal parcels were compared based on their homologous cortico-striatal connectivity. The goal here was to identify striatal parcels whose connectivity is human-specific compared to macaque parcels. Functional parcellation revealed that the human striatum was split into dorsal, dorsomedial, and rostral caudate and ventral, central, and caudal putamen, while the macaque striatum was divided into dorsal, and rostral caudate and rostral, and caudal putamen. Cross-species comparison indicated dissimilar cortico-striatal RSFC of the topographically similar dorsal caudate. We probed clinical relevance of the striatal clusters by examining differences in their cortico-striatal RSFC and gray matter (GM) volume between patients (with PD and SCZ) and healthy controls. We found abnormal RSFC not only between dorsal caudate, but also between rostral caudate, ventral, central and caudal putamen and widespread cortical regions for both PD and SCZ patients. Also, we observed significant structural atrophy in rostral caudate, ventral and central putamen for both PD and SCZ while atrophy in the dorsal caudate was specific to PD. Taken together, our cross-species comparative results revealed shared and human-specific RSFC of different striatal clusters reinforcing the complex organization and function of the striatum. In addition, we provided a testable hypothesis that abnormalities in a region with human-specific connectivity, i.e., dorsal caudate, might be associated with neuropsychiatric disorders.

DarkASDNet: Classification of ASD on Functional MRI Using Deep Neural Network

Non-invasive whole-brain scans aid the diagnosis of neuropsychiatric disorder diseases such as autism, dementia, and brain cancer. The assessable analysis for autism spectrum disorders (ASD) is rationally challenging due to the limitations of publicly available datasets. For diagnostic or prognostic tools, functional Magnetic Resonance Imaging (fMRI) exposed affirmation to the biomarkers in neuroimaging research because of fMRI pickup inherent connectivity between the brain and regions. There are profound studies in ASD with introducing machine learning or deep learning methods that have manifested advanced steps for ASD predictions based on fMRI data. However, utmost antecedent models have an inadequacy in their capacity to manipulate performance metrics such as accuracy, precision, recall, and F1-score. To overcome these problems, we proposed an avant-garde DarkASDNet, which has the competence to extract features from a lower level to a higher level and bring out promising results. In this work, we considered 3D fMRI data to predict binary classification between ASD and typical control (TC). Firstly, we pre-processed the 3D fMRI data by adopting proper slice time correction and normalization. Then, we introduced a novel DarkASDNet which surpassed the benchmark accuracy for the classification of ASD. Our model's outcomes unveil that our proposed method established state-of-the-art accuracy of 94.70% to classify ASD vs. TC in ABIDE-I, NYU dataset. Finally, we contemplated our model by performing evaluation metrics including precision, recall, F1-score, ROC curve, and AUC score, and legitimize by distinguishing with recent literature descriptions to vindicate our outcomes. The proposed DarkASDNet architecture provides a novel benchmark approach for ASD classification using fMRI processed data.