The Allen Mouse Brain Common Coordinate Framework: A 3D Reference Atlas.

Recent large-scale collaborations are generating major surveys of cell types and connections in the mouse brain, collecting large amounts of data across modalities, spatial scales, and brain areas. Successful integration of these data requires a standard 3D reference atlas. Here, we present the Allen Mouse Brain Common Coordinate Framework (CCFv3) as such a resource. We constructed an average template brain at 10 µm voxel resolution by interpolating high resolution in-plane serial two-photon tomography images with 100 µm z-sampling from 1,675 young adult C57BL/6J mice. Then, using multimodal reference data, we parcellated the entire brain directly in 3D, labeling every voxel with a brain structure spanning 43 isocortical areas and their layers, 329 subcortical gray matter structures, 81 fiber tracts, and 8 ventricular structures. CCFv3 can be used to analyze, visualize, and integrate multimodal and multiscale datasets in 3D and is openly accessible (https://atlas.brain-map.org/).

Functional parcellation of the neonatal cortical surface.

The cerebral cortex is organized into distinct but interconnected cortical areas, which can be defined by abrupt differences in patterns of resting state functional connectivity (FC) across the cortical surface. Such parcellations of the cortex have been derived in adults and older infants, but there is no widely used surface parcellation available for the neonatal brain. Here, we first demonstrate that existing parcellations, including surface-based parcels derived from older samples as well as volume-based neonatal parcels, are a poor fit for neonatal surface data. We next derive a set of 283 cortical surface parcels from a sample of n = 261 neonates. These parcels have highly homogenous FC patterns and are validated using three external neonatal datasets. The Infomap algorithm is used to assign functional network identities to each parcel, and derived networks are consistent with prior work in neonates. The proposed parcellation may represent neonatal cortical areas and provides a powerful tool for neonatal neuroimaging studies.

The cytoarchitectonic landscape revealed by deep learning method facilitated precise positioning in mouse neocortex.

Neocortex is a complex structure with different cortical sublayers and regions. However, the precise positioning of cortical regions can be challenging due to the absence of distinct landmarks without special preparation. To address this challenge, we developed a cytoarchitectonic landmark identification pipeline. The fluorescence micro-optical sectioning tomography method was employed to image the whole mouse brain stained by general fluorescent nucleotide dye. A fast 3D convolution network was subsequently utilized to segment neuronal somas in entire neocortex. By approach, the cortical cytoarchitectonic profile and the neuronal morphology were analyzed in 3D, eliminating the influence of section angle. And the distribution maps were generated that visualized the number of neurons across diverse morphological types, revealing the cytoarchitectonic landscape which characterizes the landmarks of cortical regions, especially the typical signal pattern of barrel cortex. Furthermore, the cortical regions of various ages were aligned using the generated cytoarchitectonic landmarks suggesting the structural changes of barrel cortex during the aging process. Moreover, we observed the spatiotemporally gradient distributions of spindly neurons, concentrated in the deep layer of primary visual area, with their proportion decreased over time. These findings could improve structural understanding of neocortex, paving the way for further exploration with this method.

A framework to improve the alignment of individual cytoarchitectonic maps of the Julich-Brain atlas using cortical folding landmarks.

The segregation of the cortical mantle into cytoarchitectonic areas provides a structural basis for the specialization of different brain regions. In vivo neuroimaging experiments can be linked to this postmortem cytoarchitectonic parcellation via Julich-Brain. This atlas embeds probabilistic maps that account for inter-individual variability in the localization of cytoarchitectonic areas in the reference spaces targeted by spatial normalization. We built a framework to improve the alignment of architectural areas across brains using cortical folding landmarks. This framework, initially designed for in vivo imaging, was adapted to postmortem histological data. We applied this to the first 14 brains used to establish the Julich-Brain atlas to infer a refined atlas with more focal probabilistic maps. The improvement achieved is significant in the primary regions and some of the associative areas. This framework also provides a tool for exploring the relationship between cortical folding patterns and cytoarchitectonic areas in different cortical regions to establish new landmarks in the remainder of the cortex.

Distinct Lateral Prefrontal Regions Are Organized in an Anterior-Posterior Functional Gradient.

The dorsolateral prefrontal cortex (dlPFC) is composed of multiple anatomically defined regions involved in higher-order cognitive processes, including working memory and selective attention. It is organized in an anterior-posterior global gradient where posterior regions track changes in the environment, whereas anterior regions support abstract neural representations. However, it remains unknown if such a global gradient results from a smooth gradient that spans regions or an emergent property arising from functionally distinct regions, that is, an areal gradient. Here, we recorded single neurons in the dlPFC of nonhuman primates trained to perform a memory-guided saccade task with an interfering distractor and analyzed their physiological properties along the anterior-posterior axis. We found that these physiological properties were best described by an areal gradient. Further, population analyses revealed that there is a distributed representation of spatial information across the dlPFC. Our results validate the functional boundaries between anatomically defined dlPFC regions and highlight the distributed nature of computations underlying working memory across the dlPFC.SIGNIFICANCE STATEMENT Activity of frontal lobe regions is known to possess an anterior-posterior functional gradient. However, it is not known whether this gradient is the result of individual brain regions organized in a gradient (like a staircase), or a smooth gradient that spans regions (like a slide). Analysis of physiological properties of individual neurons in the primate frontal regions suggest that individual regions are organized as a gradient, rather than a smooth gradient. At the population level, working memory was more prominent in posterior regions, although it was also present in anterior regions. This is consistent with the functional segregation of brain regions that is also observed in other systems (i.e., the visual system).

The 3D Structural Architecture of the Human Hand Area Is Nontopographic.

The functional topography of the human primary somatosensory cortex hand area is a widely studied model system to understand sensory organization and plasticity. It is so far unclear whether the underlying 3D structural architecture also shows a topographic organization. We used 7 Tesla (7T) magnetic resonance imaging (MRI) data to quantify layer-specific myelin, iron, and mineralization in relation to population receptive field maps of individual finger representations in Brodman area 3b (BA 3b) of human S1 in female and male younger adults. This 3D description allowed us to identify a characteristic profile of layer-specific myelin and iron deposition in the BA 3b hand area, but revealed an absence of structural differences, an absence of low-myelin borders, and high similarity of 3D microstructure profiles between individual fingers. However, structural differences and borders were detected between the hand and face areas. We conclude that the 3D structural architecture of the human hand area is nontopographic, unlike in some monkey species, which suggests a high degree of flexibility for functional finger organization and a new perspective on human topographic plasticity.SIGNIFICANCE STATEMENT Using ultra-high-field MRI, we provide the first comprehensive in vivo description of the 3D structural architecture of the human BA 3b hand area in relation to functional population receptive field maps. High similarity of precise finger-specific 3D profiles, together with an absence of structural differences and an absence of low-myelin borders between individual fingers, reveals the 3D structural architecture of the human hand area to be nontopographic. This suggests reduced structural limitations to cortical plasticity and reorganization and allows for shared representational features across fingers.

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.

Organization of thalamocortical structural covariance and a corresponding 3D atlas of the mouse thalamus.

For information from sensory organs to be processed by the brain, it is usually passed to appropriate areas of the cerebral cortex. Almost all of this information passes through the thalamus, a relay structure that reciprocally connects to the vast majority of the cortex. The thalamus facilitates this information transfer through a set of thalamocortical connections that vary in cellular structure, molecular profiles, innervation patterns, and firing rates. Additionally, corticothalamic connections allow for intracortical information transfer through the thalamus. These efferent and afferent connections between the thalamus and cortex have been the focus of many studies, and the importance of cortical connectivity in defining thalamus anatomy is demonstrated by multiple studies that parcellate the thalamus based on cortical connectivity profiles. Here, we examine correlated morphological variation between the thalamus and cortex, or thalamocortical structural covariance. For each voxel in the thalamus as a seed, we construct a cortical structural covariance map that represents correlated cortical volume variation, and examine whether high structural covariance is observed in cortical areas that are functionally relevant to the seed. Then, using these cortical structural covariance maps as features, we subdivide the thalamus into six non-overlapping regions (clusters of voxels), and assess whether cortical structural covariance is associated with cortical connectivity that specifically originates from these regions. We show that cortical structural covariance is high in areas of the cortex that are functionally related to the seed voxel, cortical structural covariance varies along cortical depth, and sharp transitions in cortical structural covariance profiles are observed when varying seed locations in the thalamus. Subdividing the thalamus based on structural covariance, we additionally demonstrate that the six thalamic clusters of voxels stratify cortical structural covariance along the dorsal-ventral, medial-lateral, and anterior-posterior axes. These cluster-associated structural covariance patterns are prominently detected in cortical regions innervated by fibers projecting out of their related thalamic subdivisions. Together, these results advance our understanding of how the thalamus and the cortex couple in their volumes. Our results indicate that these volume correlations reflect functional organization and structural connectivity, and further provides a novel segmentation of the mouse thalamus that can be used to examine thalamic structural variation and thalamocortical structural covariation in disease models.

Development of the next-generation functional neuro-cognitive imaging protocol - Part 1: A 3D sliding-window convolutional neural net for automated brain parcellation.

Functional MRI has emerged as a powerful tool to assess the severity of Post-concussion syndrome (PCS) and to provide guidance for neuro-cognitive therapists during treatment. The next-generation functional neuro-cognitive imaging protocol (fNCI2) has been developed to provide this assessment. This paper covers the first step in the analysis process, the development of a rapidly re-trainable, machine-learning, brain parcellation tool. The use of a sufficiently deep U-Net architecture encompassing a small (39 × 39 × 39 voxel input, 27 × 27 × 27 voxel output) sliding window to sample the entirety of the 3D image allows for the prediction of the entire image using only a single trained network. A large number of training, validating, and testing windows are thus generated from the 101 manually-labeled Mindboggle images, and full-image prediction is provided via a voxel-vote method using overlapping windows. Our method produces parcellated images that are highly consistent with standard atlas-based methods in under 3 min on a modern GPU, and the single network architecture allows for rapid retraining (<36 hr) as needed.

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.

Organization of the human cerebral cortex estimated within individuals: networks, global topography, and function.

The cerebral cortex is populated by specialized regions that are organized into networks. Here we estimated networks from functional MRI (fMRI) data in intensively sampled participants. The procedure was developed in two participants (scanned 31 times) and then prospectively applied to 15 participants (scanned 8-11 times). Analysis of the networks revealed a global organization. Locally organized first-order sensory and motor networks were surrounded by spatially adjacent second-order networks that linked to distant regions. Third-order networks possessed regions distributed widely throughout association cortex. Regions of distinct third-order networks displayed side-by-side juxtapositions with a pattern that repeated across multiple cortical zones. We refer to these as supra-areal association megaclusters (SAAMs). Within each SAAM, two candidate control regions were adjacent to three separate domain-specialized regions. Response properties were explored with task data. The somatomotor and visual networks responded to body movements and visual stimulation, respectively. Second-order networks responded to transients in an oddball detection task, consistent with a role in orienting to salient events. The third-order networks, including distinct regions within each SAAM, showed two levels of functional specialization. Regions linked to candidate control networks responded to working memory load across multiple stimulus domains. The remaining regions dissociated across language, social, and spatial/episodic processing domains. These results suggest that progressively higher-order networks nest outward from primary sensory and motor cortices. Within the apex zones of association cortex, there is specialization that repeatedly divides domain-flexible from domain-specialized regions. We discuss implications of these findings, including how repeating organizational motifs may emerge during development.NEW & NOTEWORTHY The organization of cerebral networks was estimated within individuals with intensive, repeat sampling of fMRI data. A hierarchical organization emerged in each individual that delineated first-, second-, and third-order cortical networks. Regions of distinct third-order association networks consistently exhibited side-by-side juxtapositions that repeated across multiple cortical zones, with clear and robust functional specialization among the embedded regions.

Altered thalamocortical tract trajectory growth with undisrupted thalamic parcellation pattern in human lissencephaly brain at mid-gestational stage.

Proper topographically organized neural connections between the thalamus and the cerebral cortex are mandatory for thalamus function. Thalamocortical (TC) fiber growth begins during the embryonic period and completes by the third trimester of gestation, so that human neonates at birth have a thalamus with a near-facsimile of adult functional parcellation. Whether congenital neocortical anomaly (e.g., lissencephaly) affects TC connection in humans is unknown. Here, via diffusion MRI fiber-tractography analysis of long-term formalin-fixed postmortem fetal brain diagnosed as lissencephaly in comparison with an age-matched normal one, we found similar topological patterns of thalamic subregions and of internal capsule parcellated by TC fibers. However, lissencephaly fetal brain showed white matter structural changes, including fewer/less organized TC fibers and optic radiations, and much less cortical plate invasion by TC fibers - particularly around the shallow central sulcus. Diffusion MRI fiber tractography of normal fetal brains at 15, 23, and 26 gestational weeks (GW) revealed dynamic volumetric change of each parcellated thalamic subregion, suggesting coupled developmental progress of the thalamus with the corresponding cortex. Moreover, from GW23 and GW26 normal fetal brains, TC endings in the cortical plate could be delineated to reflect cumulative progressive TC invasion of cortical plate. By contrast, lissencephaly brain showed a dramatic decrease in TC invasion of the cortical plate. Our study thus shows the feasibility of diffusion MRI fiber tractography in postmortem long-term formalin-fixed fetal brains to disclose the developmental progress of TC tracts coordinating with thalamic and neocortical growth both in normal and lissencephaly fetal brains at mid-gestational stage.

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

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

Diffusion-tensor-imaging 1-year-old and 2-year-old infant brain atlases with comprehensive gray and white matter labels.

Human infancy is marked by fastest postnatal brain structural changes. It also coincides with the onset of many neurodevelopmental disorders. Atlas-based automated structure labeling has been widely used for analyzing various neuroimaging data. However, the relatively large and nonlinear neuroanatomical differences between infant and adult brains can lead to significant offsets of the labeled structures in infant brains when adult brain atlas is used. Age-specific 1- and 2-year-old brain atlases covering all major gray and white matter (GM and WM) structures with diffusion tensor imaging (DTI) and structural MRI are critical for precision medicine for infant population yet have not been established. In this study, high-quality DTI and structural MRI data were obtained from 50 healthy children to build up three-dimensional age-specific 1- and 2-year-old brain templates and atlases. Age-specific templates include a single-subject template as well as two population-averaged templates from linear and nonlinear transformation, respectively. Each age-specific atlas consists of 124 comprehensively labeled major GM and WM structures, including 52 cerebral cortical, 10 deep GM, 40 WM, and 22 brainstem and cerebellar structures. When combined with appropriate registration methods, the established atlases can be used for highly accurate automatic labeling of any given infant brain MRI. We demonstrated that one can automatically and effectively delineate deep WM microstructural development from 3 to 38 months by using these age-specific atlases. These established 1- and 2-year-old infant brain DTI atlases can advance our understanding of typical brain development and serve as clinical anatomical references for brain disorders during infancy.

Graph alignment exploiting the spatial organization improves the similarity of brain networks.

Every brain is unique, having its structural and functional organization shaped by both genetic and environmental factors over the course of its development. Brain image studies tend to produce results by averaging across a group of subjects, under the common assumption that it is possible to subdivide the cortex into homogeneous areas while maintaining a correspondence across subjects. We investigate this assumption: can the structural properties of a specific region of an atlas be assumed to be the same across subjects? This question is addressed by looking at the network representation of the brain, with nodes corresponding to brain regions and edges to their structural relationships. Using an unsupervised graph matching strategy, we align the structural connectomes of a set of healthy subjects, considering parcellations of different granularity, to understand the connectivity misalignment between regions. First, we compare the obtained permutations with four different algorithm initializations: Spatial Adjacency, Identity, Barycenter, and Random. Our results suggest that applying an alignment strategy improves the similarity across subjects when the number of parcels is above 100 and when using Spatial Adjacency and Identity initialization (the most plausible priors). Second, we characterize the obtained permutations, revealing that the majority of permutations happens between neighbors parcels. Lastly, we study the spatial distribution of the permutations. By visualizing the results on the cortex, we observe no clear spatial patterns on the permutations and all the regions across the context are mostly permuted with first and second order neighbors.

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.

Measuring functional connectivity in frequency-domain helps to better characterize brain function.

Resting-state functional connectivity (FC) is widely used in multivariate pattern analysis of functional magnetic resonance imaging (fMRI), including identifying the locations of putative brain functional borders, predicting individual phenotypes, and diagnosing clinical mental diseases. However, limited attention has been paid to the analysis of functional interactions from a frequency perspective. In this study, by contrasting coherence-based and correlation-based FC with two machine learning tasks, we observed that measuring FC in the frequency domain helped to identify finer functional subregions and achieve better pattern discrimination capability relative to the temporal correlation. This study has proven the feasibility of coherence in the analysis of fMRI, and the results indicate that modeling functional interactions in the frequency domain may provide richer information than that in the time domain, which may provide a new perspective on the analysis of functional neuroimaging.

Global diversity in individualized cortical network topography.

Individualized cortical network topography (ICNT) varies between people and exhibits great variability in the association networks in the human brain. However, these findings were mainly discovered in Western populations. It remains unclear whether and how ICNT is shaped by the non-Western populations. Here, we leveraged a multisession hierarchical Bayesian model to define individualized functional networks in White American and Han Chinese populations with data from both US and Chinese Human Connectome Projects. We found that both the size and spatial topography of individualized functional networks differed between White American and Han Chinese groups, especially in the heteromodal association cortex (including the ventral attention, control, language, dorsal attention, and default mode networks). Employing a support vector machine, we then demonstrated that ethnicity-related ICNT diversity can be used to identify an individual's ethnicity with high accuracy (74%, pperm < 0.0001), with heteromodal networks contributing most to the classification. This finding was further validated through mass-univariate analyses with generalized additive models. Moreover, we reveal that the spatial heterogeneity of ethnic diversity in ICNT correlated with fundamental properties of cortical organization, including evolutionary cortical expansion, brain myelination, and cerebral blood flow. Altogether, this case study highlights a need for more globally diverse and publicly available neuroimaging datasets.

The parcellation of cingulate cortex in neonatal period based on resting-state functional MRI.

The human cingulate cortex (CC) is a complex region that is characterized by heterogeneous cytoarchitecture, connectivity, and function, and it is associated with various cognitive functions. The adult CC has been divided into various subregions, and this subdivision is highly consistent with its functional differentiation. However, only a few studies have focused on the function of neonatal CC. The aim of this study was to describe the cingulate segregation and the functional connectivity of each subdivision in full-term neonates (n = 60) based on resting-state functional magnetic resonance imaging. The neonatal CC was divided into three subregions, and each subregion showed specific connectivity patterns. The anterior cingulate cortex was mainly correlated with brain regions related to the salience (affected) network and default mode network (DMN), the midcingulate cortex was related to motor areas, and the posterior cingulate cortex was coupled with DMN. Moreover, we found that the cingulate subregions showed distinct functional profiles with major brain networks, which were defined using independent component analysis, and exhibited functional lateralization. This study provided new insights into the understanding of the functional specialization of neonatal CC, and these findings may have significant clinical implications, especially in predicting neurological disorder.

Parcellating the human brain using resting-state dynamic functional connectivity.

Brain cartography has expanded substantially over the past decade. In this regard, resting-state functional connectivity (FC) plays a key role in identifying the locations of putative functional borders. However, scant attention has been paid to the dynamic nature of functional interactions in the human brain. Indeed, FC is typically assumed to be stationary across time, which may obscure potential or subtle functional boundaries, particularly in regions with high flexibility and adaptability. In this study, we developed a dynamic FC (dFC)-based parcellation framework, established a new functional human brain atlas termed D-BFA (DFC-based Brain Functional Atlas), and verified its neurophysiological plausibility by stereo-EEG data. As the first dFC-based whole-brain atlas, the proposed D-BFA delineates finer functional boundaries that cannot be captured by static FC, and is further supported by good correspondence with cytoarchitectonic areas and task activation maps. Moreover, the D-BFA reveals the spatial distribution of dynamic variability across the brain and generates more homogenous parcels compared with most alternative parcellations. Our results demonstrate the superiority and practicability of dFC in brain parcellation, providing a new template to exploit brain topographic organization from a dynamic perspective. The D-BFA will be publicly available for download at https://github.com/sliderplm/D-BFA-618.