Individual Parcellation Research Articles

A hierarchical Bayesian brain parcellation framework for fusion of functional imaging datasets

Abstract Different task-based and resting-state imaging datasets provide complementary information about the organization of the human brain. Brain parcellations based on single datasets will therefore be biased towards the particular type of information present in each dataset. To overcome this limitation, we propose here a hierarchical Bayesian framework that can learn a probabilistic brain parcellation across numerous task-based and resting-state datasets, exploiting their combined strengths. The framework is partitioned into a spatial arrangement model that defines the probability of each voxel belonging to a specific parcel (the probabilistic group atlas), and a set of dataset-specific emission models that define the probability of the observed data given the parcel of the voxel. Using the human cerebellum as an example, we show that the framework optimally combines information from different datasets to achieve a new population-based atlas that outperforms atlases based on single datasets. Furthermore, we demonstrate that using only 10 min of individual data, the framework is able to generate individual brain parcellations that outperform group atlases.

D-MHGCN: An End-to-End Individual Behavioral Prediction Model Using Dual Multi-Hop Graph Convolutional Network.

Predicting individual behavior is a crucial area of research in neuroscience. Graph Neural Networks (GNNs), as powerful tools for extracting graph-structured features, are increasingly being utilized in various functional connectivity (FC) based behavioral prediction tasks. However, current predictive models primarily focus on enhancing GNNs' ability to extract features from FC networks while neglecting the importance of upstream individual network construction quality. This oversight results in constructed functional networks that fail to adequately represent individual behavioral capacity, thereby affecting the subsequent prediction accuracy. To address this issue, we proposed a new GNN-based behavioral prediction framework, named Dual Multi-Hop Graph Convolutional Network (D-MHGCN). Through the joint training of two GCNs, this framework integrates individual functional network construction and behavioral prediction into a unified optimization model. It allows the model to dynamically adjust the individual functional cortical parcellation according to the downstream tasks, thus creating task-aware, individual-specific FCNs that largely enhance its ability to predict behavior scores. Additionally, we employed multi-hop graph convolution layers instead of traditional single-hop methods in GCN to capture complex hierarchical connectivity patterns in brain networks. Our experimental evaluations, conducted on the large, public Human Connectome Project dataset, demonstrate that our proposed method outperforms existing methods in various behavioral prediction tasks. Moreover, it produces more functionally homogeneous cortical parcellation, showcasing its practical utility and effectiveness. Our work not only enhances the accuracy of individual behavioral prediction but also provides deeper insights into the neural mechanisms underlying individual differences in behavior.

Brain network dynamics in patients with single- and multiple-domain amnestic mild cognitive impairment.

Brain network dynamics have been extensively explored in patients with amnestic mild cognitive impairment (aMCI); however, differences in single- and multiple-domain aMCI (SD-aMCI and MD-aMCI) remain unclear. Using multicenter datasets, coactivation patterns (CAPs) were constructed and compared among normal control (NC), SD-aMCI, MD-aMCI, and Alzheimer's disease (AD) patients based on individual high-order cognitive network (HOCN) and primary sensory network (PSN) parcellations. Correlations between spatiotemporal characteristics and neuropsychological scores were analyzed. Compared to NC, SD-aMCI showed temporal alterations in HOCN-dominant CAPs, while MD-aMCI showed alterations in PSN-dominant CAPs. In addition, transitions from SD-aMCI to AD may involve PSN, while MD-aMCI to AD involves both PSN and HOCN. Results were generally consistent across datasets from Chinese and White populations. The HOCN and PSN are distinctively involved in aMCI subtypes and in the transformation between aMCI subtypes and AD, highlighting the necessity of aMCI subtype classification in AD studies. Individual functional network parcellations and coactivation pattern (CAP) analysis were performed to characterize spatiotemporal differences between single- and multiple-domain amnestic mild cognitive impairment (SD-aMCI and MD-aMCI), and between distinct aMCI subtypes and Alzheimer's disease (AD). The analysis of multicenter datasets converged on four pairs of recurrent CAPs, including primary sensory networks (PSN)-dominant CAPs, high-order cognitive networks (HOCN)-dominant CAPs, and PSN-HOCN-interacting CAPs. The HOCN and PSN are distinctively involved in aMCI subtypes and in the transformation between distinct aMCI subtypes and AD.

Multimodal Connectivity-Based Individual Parcellation and Analysis for Humans and Rhesus Monkeys.

Individual brains vary greatly in morphology, connectivity and organization. Individualized brain parcellation is capable of precisely localizing subject-specific functional regions. However, most individualization approaches have examined single modalities of data and have not generalized to nonhuman primates. The present study proposed a novel multimodal connectivity-based individual parcellation (MCIP) method, which optimizes within-region homogeneity, spatial continuity and similarity to a reference atlas with the fusion of personal functional and anatomical connectivity. Comprehensive evaluation demonstrated that MCIP outperformed state-of-the-art multimodal individualization methods in terms of functional and anatomical homogeneity, predictability of cognitive measures, heritability, reproducibility and generalizability across species. Comparative investigation showed a higher topographic variability in humans than that in macaques. Therefore, MCIP provides improved accurate and reliable mapping of brain functional regions over existing methods at an individual level across species, and could facilitate comparative and translational neuroscience research.

TS-AI: A deep learning pipeline for multimodal subject-specific parcellation with task contrasts synthesis

Accurate mapping of brain functional subregions at an individual level is crucial. Task-based functional MRI (tfMRI) captures subject-specific activation patterns during various functions and behaviors, facilitating the individual localization of functionally distinct subregions. However, acquiring high-quality tfMRI is time-consuming and resource-intensive in both scientific and clinical settings. The present study proposes a two-stage network model, TS-AI, to individualize an atlas on cortical surfaces through the prediction of tfMRI data. TS-AI first synthesizes a battery of task contrast maps for each individual by leveraging tract-wise anatomical connectivity and resting-state networks. These synthesized maps, along with feature maps of tract-wise anatomical connectivity and resting-state networks, are then fed into an end-to-end deep neural network to individualize an atlas. TS-AI enables the synthesized task contrast maps to be used in individual parcellation without the acquisition of actual task fMRI scans. In addition, a novel feature consistency loss is designed to assign vertices with similar features to the same parcel, which increases individual specificity and mitigates overfitting risks caused by the absence of individual parcellation ground truth. The individualized parcellations were validated by assessing test-retest reliability, homogeneity, and cognitive behavior prediction using diverse reference atlases and datasets, demonstrating the superior performance and generalizability of TS-AI. Sensitivity analysis yielded insights into region-specific features influencing individual variation in functional regionalization. Additionally, TS-AI identified accelerated shrinkage in the medial temporal and cingulate parcels during the progression of Alzheimer's disease, suggesting its potential in clinical research and applications.

Individual patterns of functional connectivity in neonates as revealed by surface-based Bayesian modeling.

Resting-state functional connectivity is a widely used approach to study the functional brain network organization during early brain development. However, the estimation of functional connectivity networks in individual infants has been rather elusive due to the unique challenges involved with functional magnetic resonance imaging (fMRI) data from young populations. Here, we use fMRI data from the developing Human Connectome Project (dHCP) database to characterize individual variability in a large cohort of term-born infants (N = 289) using a novel data-driven Bayesian framework. To enhance alignment across individuals, the analysis was conducted exclusively on the cortical surface, employing surface-based registration guided by age-matched neonatal atlases. Using 10 minutes of resting-state fMRI data, we successfully estimated subject-level maps for fourteen brain networks/subnetworks along with individual functional parcellation maps that revealed differences between subjects. We also found a significant relationship between age and mean connectivity strength in all brain regions, including previously unreported findings in higher-order networks. These results illustrate the advantages of surface-based methods and Bayesian statistical approaches in uncovering individual variability within very young populations.

Predicting Cognitive Functioning for Patients with a High-Grade Glioma: Evaluating Different Representations of Tumor Location in a Common Space

Cognitive functioning is increasingly considered when making treatment decisions for patients with a brain tumor in view of a personalized onco-functional balance. Ideally, one can predict cognitive functioning of individual patients to make treatment decisions considering this balance. To make accurate predictions, an informative representation of tumor location is pivotal, yet comparisons of representations are lacking. Therefore, this study compares brain atlases and principal component analysis (PCA) to represent voxel-wise tumor location. Pre-operative cognitive functioning was predicted for 246 patients with a high-grade glioma across eight cognitive tests while using different representations of voxel-wise tumor location as predictors. Voxel-wise tumor location was represented using 13 different frequently-used population average atlases, 13 randomly generated atlases, and 13 representations based on PCA. ElasticNet predictions were compared between representations and against a model solely using tumor volume. Preoperative cognitive functioning could only partly be predicted from tumor location. Performances of different representations were largely similar. Population average atlases did not result in better predictions compared to random atlases. PCA-based representation did not clearly outperform other representations, although summary metrics indicated that PCA-based representations performed somewhat better in our sample. Representations with more regions or components resulted in less accurate predictions. Population average atlases possibly cannot distinguish between functionally distinct areas when applied to patients with a glioma. This stresses the need to develop and validate methods for individual parcellations in the presence of lesions. Future studies may test if the observed small advantage of PCA-based representations generalizes to other data.

Reduced lateralization of multiple functional brain networks in autistic males

BackgroundAutism spectrum disorder has been linked to a variety of organizational and developmental deviations in the brain. One such organizational difference involves hemispheric lateralization, which may be localized to language-relevant regions of the brain or distributed more broadly.Methods In the present study, we estimated brain hemispheric lateralization in autism based on each participant’s unique functional neuroanatomy rather than relying on group-averaged data. Additionally, we explored potential relationships between the lateralization of the language network and behavioral phenotypes including verbal ability, language delay, and autism symptom severity. We hypothesized that differences in hemispheric asymmetries in autism would be limited to the language network, with the alternative hypothesis of pervasive differences in lateralization. We tested this and other hypotheses by employing a cross-sectional dataset of 118 individuals (48 autistic, 70 neurotypical). Using resting-state fMRI, we generated individual network parcellations and estimated network asymmetries using a surface area-based approach. A series of multiple regressions were then used to compare network asymmetries for eight significantly lateralized networks between groups.ResultsWe found significant group differences in lateralization for the left-lateralized Language (d = -0.89), right-lateralized Salience/Ventral Attention-A (d = 0.55), and right-lateralized Control-B (d = 0.51) networks, with the direction of these group differences indicating less asymmetry in autistic males. These differences were robust across different datasets from the same participants. Furthermore, we found that language delay stratified language lateralization, with the greatest group differences in language lateralization occurring between autistic males with language delay and neurotypical individuals.ConclusionsThese findings evidence a complex pattern of functional lateralization differences in autism, extending beyond the Language network to the Salience/Ventral Attention-A and Control-B networks, yet not encompassing all networks, indicating a selective divergence rather than a pervasive one. Moreover, we observed an association between Language network lateralization and language delay in autistic males.

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.

Prior-guided individualized thalamic parcellation based on local diffusion characteristics.

Comprising numerous subnuclei, the thalamus intricately interconnects the cortex and subcortex, orchestrating various facets of brain functions. Extracting personalized parcellation patterns for these subnuclei is crucial, as different thalamic nuclei play varying roles in cognition and serve as therapeutic targets for neuromodulation. However, accurately delineating the thalamic nuclei boundary at the individual level is challenging due to intersubject variability. In this study, we proposed a prior-guided parcellation (PG-par) method to achieve robust individualized thalamic parcellation based on a central-boundary prior. We first constructed probabilistic atlas of thalamic nuclei using high-quality diffusion MRI datasets based on the local diffusion characteristics. Subsequently, high-probability voxels in the probabilistic atlas were utilized as prior guidance to train unique multiple classification models for each subject based on a multilayer perceptron. Finally, we employed the trained model to predict the parcellation labels for thalamic voxels and construct individualized thalamic parcellation. Through a test-retest assessment, the proposed prior-guided individualized thalamic parcellation exhibited excellent reproducibility and the capacity to detect individual variability. Compared with group atlas registration and individual clustering parcellation, the proposed PG-par demonstrated superior parcellation performance under different scanning protocols and clinic settings. Furthermore, the prior-guided individualized parcellation exhibited better correspondence with the histological staining atlas. The proposed prior-guided individualized thalamic parcellation method contributes to the personalized modeling of brain parcellation.

Functional connectivity analyses of individual hippocampal subregions in major depressive disorder with electroconvulsive therapy.

The hippocampus has been widely reported to be involved in the neuropathology of major depressive disorder (MDD). All the previous researches adopted group-level hippocampus subregions atlas to investigate abnormal functional connectivities in MDD in absence of capturing individual variability. In addition, the molecular basis of functional impairments of hippocampal subregions in MDD remains elusive. We aimed to reveal functional disruptions and recovery of individual hippocampal subregions in MDD patients before and after ECT and linked these functional connectivity differences to transcriptomic profiles to reveal molecular mechanism. we used group guided individual functional parcellation approach to define individual subregions of hippocampus for each participant. Resting-state functional connectivity (FC) analysis of individual hippocampal subregions was conducted to investigate functional disruptions and recovery in MDD patients before and after ECT. Spatial association between functional connectivity differences and transcriptomic profiles was employed to reveal molecular mechanism. MDD patients showed increased FCs of the left tail part of hippocampus with dorsolateral prefrontal cortex and middle temporal gyrus while decreased FC with primary visual cortex. These abnormal FCs in MDD patients were normalized after ECT. In addition, we found that functional disruptions of the left tail part of hippocampus in MDD were mainly related to synaptic signaling and transmission, ion transport, cell-cell signaling and neurogenesis. Our findings provide initial evidence for functional connectome disruption of individual hippocampal subregions and their molecular basis in MDD.

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

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

Functional connectivity reveals different brain networks underlying the idiopathic foreign accent syndrome.

Foreign accent syndrome (FAS) is characterized by new onset speech that is perceived as foreign. Available data from acquired cases suggests focal brain damage in language and sensorimotor brain networks, but little remains known about abnormal functional connectivity in idiopathic cases of FAS without structural damage. Here, connectomic analyses were completed on three patients with idiopathic FAS to investigate unique functional connectivity abnormalities underlying accent change for the first time. Machine learning (ML)-based algorithms generated personalized brain connectomes based on a validated parcellation scheme from theHuman Connectome Project (HCP). Diffusion tractography was performed on each patient to rule out structural fiber damage to the language system. Resting-state-fMRI was assessed with ML-based software to examine functional connectivity between individual parcellations within language and sensorimotor networks and subcortical structures. Functional connectivity matrices were created and compared against a dataset of 200 healthy subjects to identify abnormally connected parcellations. Three female patients (28-42 years) who presented with accent changes from Australian English to Irish (n = 2) or American English to British English (n = 1) demonstrated fully intact language system structural connectivity. All patients demonstrated functional connectivity anomalies within language and sensorimotor networks in numerous left frontal regions and between subcortical structures in one patient. Few commonalities in functional connectivity anomalies were identified between all three patients, specifically 3 internal-network parcellation pairs. No common inter-network functional connectivity anomalies were identified between all patients. The current study demonstrates specific language, and sensorimotor functional connectivity abnormalities can exist and be quantitatively shown in the absence of structural damage for future study.

HOA2.0-ComPaRe: A next generation Harvard-Oxford Atlas comparative parcellation reasoning method for human and macaque individual brain parcellation and atlases of the cerebral cortex.

Comparative structural neuroanatomy is a cornerstone for understanding human brain structure and function. A parcellation framework that relates systematically to fundamental principles of histological organization is an essential step in generating structural comparisons between species. In the present investigation, we developed a comparative parcellation reasoning system (ComPaRe), which is a formal ontological system in human and non-human primate brains based on the cortical cytoarchitectonic mapping used for both species as detailed by Brodmann. ComPaRe provides a theoretical foundation for mapping neural systems in humans and other species using neuroimaging. Based on this approach, we revised the methodology of the original Harvard-Oxford Atlas (HOA) system of brain parcellation to produce a comparative framework for the human (hHOA) and the rhesus monkey (mHOA) brains, which we refer to as HOA2.0-ComPaRe. In addition, we used dedicated segmentation software in the publicly available 3D Slicer platform to parcellate an individual human and rhesus monkey brain. This method produces quantitative morphometric parcellations in the individual brains. Based on these parcellations we created a representative template and 3D brain atlas for the two species, each based on a single subject. Thus, HOA2.0-ComPaRe provides a theoretical foundation for mapping neural systems in humans and other species using neuroimaging, while also representing a significant revision of the original human and macaque monkey HOA parcellation schemas. The methodology and atlases presented here can be used in basic and clinical neuroimaging for morphometric (volumetric) analysis, further generation of atlases, as well as localization of function and structural lesions.

Individual functional parcellation revealed compensation of dynamic limbic network organization in healthy ageing.

Using group-level functional parcellations and constant-length sliding window analysis, dynamic functional connectivity studies have revealed network-specific impairment and compensation in healthy ageing. However, functional parcellation and dynamic time windows vary across individuals; individual-level ageing-related brain dynamics are uncertain. Here, we performed individual parcellation and individual-length sliding window clustering to characterize ageing-related dynamic network changes. Healthy participants (n= 637, 18-88 years) from the Cambridge Centre for Ageing and Neuroscience dataset were included. An individual seven-network parcellation, varied from group-level parcellation, was mapped for each participant. For each network, strong and weak cognitive brain states were revealed by individual-length sliding window clustering and canonical correlation analysis. The results showed negative linear correlations between age and change ratios of sizes in the default mode, frontoparietal, and salience networks and a positive linear correlation between age and change ratios of size in the limbic network (LN). With increasing age, the occurrence and dwell time of strong states showed inverted U-shaped patterns or a linear decreasing pattern in most networks but showed a linear increasing pattern in the LN. Overall, this study reveals a compensative increase in emotional networks (i.e., the LN) and a decline in cognitive and primary sensory networks in healthy ageing. These findings may provide insights into network-specific and individual-level targeting during neuromodulation in ageing and ageing-related diseases.

Atlas-guided parcellation: Individualized functionally-homogenous parcellation in cerebral cortex

Resting-state Magnetic resonance imaging-based parcellation aims to group the voxels/vertices non-invasively based on their connectivity profiles, which has achieved great success in understanding the fundamental organizational principles of the human brain. Given the substantial inter-individual variability, the increasing number of studies focus on individual parcellation. However, current methods perform individual parcellations independently or are based on the group prior, requiring expensive computational costs, precise parcel alignment, and extra group information. In this work, an efficient and flexible parcellation framework of individual cerebral cortex was proposed based on a region growing algorithm by merging the unassigned and neighbor vertex with the highest-correlated parcel iteratively. It considered both consistency with prior atlases and individualized functional homogeneity of parcels, which can be applied to a single individual without parcel alignment and group information. The proposed framework was leveraged to 100 unrelated subjects for functional homogeneity comparison and individual identification, and 186 patients with Parkison's disease for symptom prediction. Results demonstrated our framework outperformed other methods in functional homogeneity, and the generated parcellations provided 100% individual identification accuracy. Moreover, the default mode network (DMN) exhibited higher functional homogeneity, intra-subject parcel reproducibility and fingerprinting accuracy, while the sensorimotor network did the opposite, reflecting that the DMN is the most representative, stable, and individual-identifiable network in the resting state. The correlation analysis showed that the severity of the disease symptoms was related negatively to the similarity of individual parcellation and the atlases of healthy populations. The disease severity can be correctly predicted using machine learning models based on individual topographic features such as parcel similarity and parcel size. In summary, the proposed framework not only significantly improves the functional homogeneity but also captures individualized and disease-related brain topography, serving as a potential tool to explore brain function and disease in the future.

Automatic method for individual parcellation of manganese-enhanced magnetic resonance imaging of rat brain.

AimsTo construct an automatic method for individual parcellation of manganese-enhanced magnetic resonance imaging (MEMRI) of rat brain with high accuracy, which could preserve the inherent voxel intensity and Regions of interest (ROI) morphological characteristics simultaneously.Methods and resultsThe transformation relationship from standardized space to individual space was obtained by firstly normalizing individual image to the Paxinos space and then inversely transformed. On the other hand, all the regions defined in the atlas image were separated and resaved as binary mask images. Then, transforming the mask images into individual space via the inverse transformations and reslicing using the 4th B-spline interpolation algorithm. The boundary of these transformed regions was further refined by image erosion and expansion operator, and finally combined together to generate the individual parcellations. Moreover, two groups of MEMRI images were used for evaluation. We found that the individual parcellations were satisfied, and the inherent image intensity was preserved. The statistical significance of case-control comparisons was further optimized.ConclusionsWe have constructed a new automatic method for individual parcellation of rat brain MEMRI images, which could preserve the inherent voxel intensity and further be beneficial in case-control statistical analyses. This method could also be extended to other imaging modalities, even other experiments species. It would facilitate the accuracy and significance of ROI-based imaging analyses.

Source imaging of high-density visual evoked potentials with multi-scale brain parcellations and connectomes

We describe the multimodal neuroimaging dataset VEPCON (OpenNeuro Dataset ds003505). It includes raw data and derivatives of high-density EEG, structural MRI, diffusion weighted images (DWI) and single-trial behavior (accuracy, reaction time). Visual evoked potentials (VEPs) were recorded while participants (n = 20) discriminated briefly presented faces from scrambled faces, or coherently moving stimuli from incoherent ones. EEG and MRI were recorded separately from the same participants. The dataset contains raw EEG and behavioral data, pre-processed EEG of single trials in each condition, structural MRIs, individual brain parcellations at 5 spatial resolutions (83 to 1015 regions), and the corresponding structural connectomes computed from fiber count, fiber density, average fractional anisotropy and mean diffusivity maps. For source imaging, VEPCON provides EEG inverse solutions based on individual anatomy, with Python and Matlab scripts to derive activity time-series in each brain region, for each parcellation level. The BIDS-compatible dataset can contribute to multimodal methods development, studying structure-function relations, and to unimodal optimization of source imaging and graph analyses, among many other possibilities.

Group‐guided individual functional parcellation of the hippocampus and application to normal aging

Aging is closely associated with cognitive decline affecting attention, memory and executive functions. The hippocampus is the core brain area for human memory, learning, and cognition processing. To delineate the individual functional patterns of hippocampus is pivotal to reveal the neural basis of aging. In this study, we developed a group‐guided individual parcellation approach based on semisupervised affinity propagation clustering using the resting‐state functional magnetic resonance imaging to identify individual functional subregions of hippocampus and to identify the functional patterns of each subregion during aging. A three‐way group parcellation was yielded and was taken as prior information to guide individual parcellation of hippocampus into head, body, and tail in each subject. The superiority of individual parcellation of hippocampus is validated by higher intraregional functional similarities by compared to group‐level parcellation results. The individual variations of hippocampus were associated with coactivation patterns of three typical functions of hippocampus. Moreover, the individual functional connectivities of hippocampus subregions with predefined target regions could better predict age than group‐level functional connectivities. Our study provides a novel framework for individual brain functional parcellations, which may facilitate the future individual researches for brain cognitions and brain disorders and directing accurate neuromodulation.

Tractography density affects whole-brain structural architecture and resting-state dynamical modeling

Dynamical modeling of the resting-state brain dynamics essentially relies on the empirical neuroimaging data utilized for the model derivation and validation. There is however still no standardized data processing for magnetic resonance imaging pipelines and the structural and functional connectomes involved in the models. In this study, we thus address how the parameters of diffusion-weighted data processing for structural connectivity (SC) can influence the validation results of the whole-brain mathematical models informed by SC. For this, we introduce a set of simulation conditions including the varying number of total streamlines of the whole-brain tractography (WBT) used for extraction of SC, cortical parcellations based on functional and anatomical brain properties and distinct model fitting modalities. The main objective of this study is to explore how the quality of the model validation can vary across the considered simulation conditions. We observed that the graph-theoretical network properties of structural connectome can be affected by varying tractography density and strongly relate to the model performance. We also found that the optimal number of the total streamlines of WBT can vary for different brain atlases. Consequently, we suggest a way how to improve the model performance based on the network properties and the optimal parameter configurations from multiple WBT conditions. Furthermore, the population of subjects can be stratified into subgroups with divergent behaviors induced by the varying WBT density such that different recommendations can be made with respect to the data processing for individual subjects and brain parcellations.