Human brain mapping: A systematic comparison of parcellation methods for the human cerebral cortex

Neuroimaging pattern interactions for suicide risk in depression captured by ensemble learning over transcriptome-defined parcellation

BACKGROUND

For suicide in major depression disorder, it is urgent to seek for a reliable neuroimaging biomarker with interpretable links to molecular tissue signatures. Accordingly, we developed an ensemble learning scheme over transcriptome-defined parcellations (TDP) to explore homogeneously parcellated brain patterns and their interactions.

METHODS

96 depressed patients without suicide attempt (SA), 86 with SA and 102 healthy controls were recruited for resting state fMRI scanning. Six genetic dimensions were created by homogenous transcriptomic delineations from Allen Human Brain Atlas. Spatially-continuous TDPs were generated according to expression-levels of each brain region along diverse dimensions. Subsequently, TDPs were integrated with a three-layer ensemble learning scheme, where brain dysfunction of each TDP related to suicide was quantified with a resting-state functional abnormality (RSFA) score. Then, personalized index of brain dysfunction was produced according to the interactive pattern across TDPs.

RESULTS

Ensemble learning over TDPs displayed higher suicide predictive performance, relative to that over the regions level, and over null model (95 % CI of accuracy: 73.23 ± 1.07 %; 64.59 ± 3.00 %; 65.41 ± 3.97 %, respectively). Empowered by specific parieto-occipital TDP (PO-TDP) pattern quantified with RSFA score in suicide risk prediction, its alternations of SA effects were spatially associated with transcriptional profiles of GRIN2A and GABRG2. Moreover, glutamatergic and GABAergic synapse were overrepresented in enrichment analysis.

CONCLUSION

Glutamatergic and GABAergic dysfunction in the visual cortex was suggested via the PO-TDP specific interaction pattern. The inherent excitatory/inhibitory imbalance could contribute to aberrant emotional processing and neurocognitive impairment, ultimately leading to suicide.

Direct Comparison of EEG Resting State and Task Functional Connectivity Patterns for Predicting Working Memory Performance Using Connectome-Based Predictive Modeling

Background: The integration of machine learning with advanced neuroimaging has emerged as a powerful approach for uncovering the relationship between neuronal activity patterns and behavioral traits. While resting-state neuroimaging has significantly contributed to understanding the neural basis of cognition, recent fMRI studies suggest that task-based paradigms may offer superior predictive power for cognitive outcomes. However, this hypothesis has never been tested using electroencephalography (EEG) data. Methods: We conducted the first experimental comparison of predictive models built on high-density EEG data recorded during both resting-state and an auditory working memory task. Multiple data processing pipelines were employed to ensure robustness and reliability. Model performance was evaluated by computing the Pearson correlation coefficient between predicted and observed behavioral scores, supplemented by mean absolute error and root mean square error metrics for each model configuration. Results: Consistent with prior fMRI findings, task-based EEG data yielded slightly better modeling performance than resting-state data. Both conditions demonstrated high predictive accuracy, with peak correlations between observed and predicted values reaching r = 0.5. Alpha and beta band functional connectivity were the strongest predictors of working memory performance, followed by theta and gamma bands. Additionally, the choice of parcellation atlas and connectivity method significantly influenced results, highlighting the importance of methodological considerations. Conclusion: Our findings support the advantage of task-based EEG over resting-state data in predicting cognitive performance, aligning with. The study underscores the critical role of frequency-specific functional connectivity and methodological choices in model performance. These insights should guide future experimental designs in cognitive neuroscience. Impact Statement This study provides the first direct comparison of EEG-based functional connectivity during rest and task conditions for predicting working memory performance using connectome-based predictive modeling (CPM). It demonstrates that task-based EEG data slightly outperforms resting-state data, with alpha and beta bands being the most predictive. The findings highlight the critical influence of methodological choices, such as parcellation atlases and connectivity metrics, on model outcomes. By bridging gaps in EEG research and validating CPM's applicability, this work advances the optimization of neuroimaging protocols for cognitive assessment, offering insights for future studies in cognitive neuroscience.

The generalizability of cortical area parcellations across early childhood

The cerebral cortex consists of distinct areas that develop through intrinsic embryonic patterning and postnatal experiences. Accurate parcellation of these areas in neuroimaging studies improves statistical power and cross-study comparability. Given significant brain changes in volume, microstructure, and connectivity during early life, we hypothesized that cortical areas in 1- to 3-year-olds would differ markedly from neonates and increasingly resemble adult patterns as development progresses. Here, we parcellated the cerebral cortex into putative areas using local functional connectivity (FC) gradients in 92 toddlers at 2 years old. We demonstrate high reproducibility of these cortical areas across 1- to 3-year-olds in two independent datasets. The area boundaries in 1- to 3-year-olds were more similar to those in adults than those in neonates. While the age-specific group area parcellation better fits the underlying FC in individuals during the first 3 years, adult area parcellations still have utility in developmental studies, especially in children older than 6 years. Additionally, we provide connectivity-based community assignments of the area parcels, showing fragmented anterior and posterior components based on the strongest connectivity, yet alignment with adult systems when weaker connectivity was included.

fMRI-based data-driven brain parcellation using independent component analysis

BACKGROUND

Studies using functional magnetic resonance imaging (fMRI) broadly require a method of parcellating the brain into regions of interest (ROIs). Parcellations can be based on standardized brain anatomy, such as the Montreal Neurological Institute's (MNI) 152 atlas, or an individual's functional activity patterns, such as the Personode software.

NEW METHOD

This work outlines and tests the independent component analysis (ICA)-based parcellation algorithm (IPA) when applied to a hypertension study (n = 48) that uses the independent components (ICs) output from group ICA (gICA) to build ROIs which are ideally spatially consistent and functionally homogeneous. After regression of ICs to all subjects, the IPA builds individualized parcellations while simultaneously obtaining a gICA-derived parcellation.

RESULTS

ROI spatial consistency quantified by dice similarity coefficients (DSCs) show individualized parcellations exhibit mean DSCs of 0.69 ± 0.14. Functional homogeneity, calculated as mean Pearson correlation value of all voxels comprising a ROI, shows individualized parcellations with a mean of 0.30 ± 0.14 and gICA-derived parcellations' mean of 0.38 ± 0.15.

COMPARISON WITH EXISTING METHOD(S)

Individualized Personode parcellations show decreased mean DSCs (0.43 ± 0.11) with the individualized parcellations, gICA-derived parcellations, and the MNI atlas having decreased homogeneity values of 0.28 ± 0.14, 0.31 ± 0.15, and 0.20 ± 0.11 respectively.

CONCLUSIONS

Results show that the IPA can more reliably define a ROI and does so with higher functional homogeneity. Given these findings, the IPA shows promise as a novel parcellation technique that could aid the analysis of fMRI data.

Functional-based parcellation of the mouse prefrontal cortex for network perturbation analysis

The prefrontal cortex (PFC) is a brain region involved in higher-order cognitive processes such as attention, emotional regulation, and social behavior. However, the delineation of distinct subdivisions within the mouse PFC and their contributions to the broader brain network function remain debated. This study utilizes resting-state functional magnetic resonance imaging (MRI) from a cohort of 100 C57BL/6J wild-type mice to derive the functional connectivity (FC)-based parcellation of the mouse PFC with voxel resolution. Our findings reveal clusters that deviate from the established anatomical subdivisions within the cingulate and prelimbic areas while aligning in infralimbic and orbital cortices. Upon the chemogenetic perturbation of one of the clusters, FC perturbations occur only within the functional network linked to the targeted cluster and do not spread to neighboring anatomical areas or functional clusters. We propose FC-based parcellation as a valuable approach for tracking the site of activation and network impact of neurostimulation strategies.

Precision neuroregulation combining liquid metal and magnetic stimulation

BACKGROUND

Electromagnetic field-based neuroregulation technology is a crucial technique for treating central nervous system and peripheral nervous system disorders. However, the use of invasive electrodes has unavoidable problems such as the risk of inflammation due to high hardness, electrical connections and the need for batteries. On the other hand, non-invasive magnetic stimulation has limitations such as centimeter-level focal areas and shallow stimulation depth.

METHODS

To enhance the precision and effectiveness of wireless magnetic stimulation, we employed a figure-8 magnetic stimulation coil (8-coil) to generate a magnetic field, combined with an injectable, highly conductive, and flexible liquid metal (LM) to produce a millimeter-scale focused electric field. A coaxial electric field measurement electrode was used to establish an agar phantom-based electric field measurement platform. The sciatic nerve of C57 mice was stimulated under acute anesthesia conditions, and electromyography (EMG) signals were collected to evaluate the enhancement of stimulation effects. Long-term safety was assessed through four weeks of implantation.

RESULTS

Theoretical analysis and finite element simulations demonstrated that the combination of LM and the 8-coil generated a millimeter-scale enhanced vector electric field within the tissue. Measured electric field distributions closely aligned with theoretical and simulation results. In the sciatic nerve experiments on mice, 1 µL of LM under a 0.45 T magnetic field significantly increased EMG signals and leg movement amplitude by approximately 500%. Long-term implantation under magnetic stimulation revealed no adverse effects.

CONCLUSIONS

This method utilizes focused electric fields to improve the precision and effectiveness of neuro-magnetic stimulation. It holds promise as a novel approach for precise stimulation. Preliminary evidence was provided for the safety of in vivo LM implantation under external magnetic fields.

Semantic-spatial feature-fused cortical surface parcellation: a scale-unified spatial learning network with boundary contrastive loss

The cortical surface parcellation provides prior guidance for studying mental disorders and human cognition. Graph neural networks (GNNs) have gained popularity in this task to preserve its spatial structure. However, previous GNNs struggled to effectively exploit the information contained in the complex spatial structure of the cortical surface and generally encountered an uneven node distribution issue. Meanwhile, labeling boundary nodes was also identified as a widespread problem in this task. Accordingly, this paper develops a scale-unified spatial learning network with a boundary contrastive loss (SSLNet) for cortical surface parcellation. Its core is the scale-unified spatial learning module. It devises neighbor feature extraction and aggregation strategies by fully integrating spatial coordinates and semantic structure to learn effective spatial features of local neighborhoods. More importantly, spatial scale unification is incorporated into this module to mitigate the negative effect on spatial learning caused by node distribution differences among local areas. Additionally, a universal boundary contrastive loss is constructed, enhancing the feature discriminability of boundary nodes by constraining them to be close to the same class nodes and apart from different class nodes in the feature space. It considerably improves boundary performance without increasing parameters or changing the network structure. Experiments regarding public Mindboggle demonstrate that the dice score and accuracy of SSLNet achieve 89.8 % and 90.89 % , respectively, surpassing existing methods.

Hippocampal network connectivity and episodic memory in individuals aging with traumatic brain injury

Aging is associated with marked declines in episodic memory corresponding with decreased volume in studies of morphology and reduced network response in studies of functional connectomics. Furthermore, recent research has demonstrated that reductions in resting state network connectivity are related to declines in episodic memory, specifically in the default mode and frontoparietal cortical networks. Additionally, the interactive effects of aging and traumatic brain injury (TBI) are associated with increased risk for neurodegeneration and episodic memory impairments. However, there is a gap in the literature examining episodic memory and hippocampal-subcortical resting state connectivity differences related to aging with and without TBI. The current work aims to investigate episodic memory differences between older adults with TBI (N = 45) and older adults with no history of TBI (N = 28) and how that relates to hippocampal-subcortical network differences at rest. We demonstrate a positive relationship between default mode and frontoparietal network connectivity and memory performance differentially between those aging with and without moderate-severe TBI (msTBI). Additionally, we demonstrate that reliability in the strength of resting state functional connectivity between parcellations is weakest among connections to the hippocampus compared to other cortical connections but is generally reliable across other connections.

Systematic bias in surface area asymmetry measurements from automatic cortical parcellations

Anatomical asymmetry is a hallmark of the human brain and may reflect hemispheric differences in its functional organization. Widely used software like FreeSurfer can automate neuroanatomical measurements and facilitate studies of hemispheric asymmetry. However, patterns of surface area lateralization measured using FreeSurfer are curiously consistent across diverse samples. Here, we demonstrate systematic biases in these measurements obtained from the default processing pipeline. We compared surface area asymmetry measured from reconstructions of original brains vs. the same scans after flipping their left-right orientation. The default pipeline returned implausible asymmetry patterns between the original and flipped brains: Many structures were always left- or right-lateralized. Notably, these biases occur prominently in key speech and language regions. In contrast, manual labeling and curvature-based parcellations of key structures both yielded the expected reversals of left/right lateralization in flipped brains. We determined that these biases result from discrepancies in how regional labels are defined in the left vs. right hemisphere in the default cortical parcellation atlases. These biases are carried into individual parcellations because the FreeSurfer parcellation algorithm prioritizes vertex correspondence to the template atlas relative to individual neuroanatomical variation. We further demonstrate several straightforward, bias-free approaches to measuring surface area asymmetry, including using symmetric registration templates and parcellation atlases, vertex-wise analyses, and within-subject curvature-based parcellations. These results highlight theoretical concerns about using only the default processing stream to make inferences about population-level brain asymmetry and underscore the need for validating bias-free neuroanatomical measurements, particularly when studying regions where structural lateralization may underlie functional lateralization.

A network correspondence toolbox for quantitative evaluation of novel neuroimaging results

The brain can be decomposed into large-scale functional networks, but the specific spatial topographies of these networks and the names used to describe them vary across studies. Such discordance has hampered interpretation and convergence of research findings across the field. We have developed the Network Correspondence Toolbox (NCT) to permit researchers to examine and report spatial correspondence between their novel neuroimaging results and multiple widely used functional brain atlases. We provide several exemplar demonstrations to illustrate how researchers can use the NCT to report their own findings. The NCT provides a convenient means for computing Dice coefficients with spin test permutations to determine the magnitude and statistical significance of correspondence among user-defined maps and existing atlas labels. The adoption of the NCT will make it easier for network neuroscience researchers to report their findings in a standardized manner, thus aiding reproducibility and facilitating comparisons between studies to produce interdisciplinary insights.

The Generalizability of Cortical Area Parcellations Across Early Childhood

The cerebral cortex consists of distinct areas that develop through intrinsic embryonic patterning and postnatal experiences. Accurate parcellation of these areas in neuroimaging studies improves statistical power and cross-study comparability. Given significant brain changes in volume, microstructure, and connectivity during early life, we hypothesized that cortical areas in 1- to 3-year-olds would differ markedly from neonates and increasingly resemble adult patterns as development progresses. Here, we parcellated the cerebral cortex into putative areas using local functional connectivity gradients in 92 toddlers at 2 years old. We demonstrate high reproducibility of these cortical regions across 1- to 3-year-olds in two independent datasets. The area boundaries in 1- to 3-year-olds were more similar to those in adults than those in neonates. While the age-specific group area parcellation better fit the underlying functional connectivity in individuals during the first 3 years, adult area parcellations might still have some utility in developmental studies, especially in children older than 6 years. Additionally, we provide connectivity-based community assignments of the parcels, showing fragmented anterior and posterior components based on the strongest connectivity, yet alignment with adult systems when weaker connectivity was included.

Individualized brain radiomics-based network tracks distinct subtypes and abnormal patterns in prodromal Parkinson's disease

Individuals in the prodromal phase of Parkinson's disease (PD) exhibit significant heterogeneity and can be divided into distinct subtypes based on clinical symptoms, pathological mechanisms, and brain network patterns. However, little has been done regarding the valid subtyping of prodromal PD, which hinders the early diagnosis of PD. Therefore, we aimed to identify the subtypes of prodromal PD using the brain radiomics-based network and examine the unique patterns linked to the clinical presentations of each subtype. Individualized brain radiomics-based network was constructed for normal controls (NC; N = 110), prodromal PD patients (N = 262), and PD patients (N = 108). A data-driven clustering approach using the radiomics-based network was carried out to cluster prodromal PD patients into higher-/lower-risk subtypes. Then, the dissociated patterns of clinical manifestations, anatomical structure alterations, and gene expression between these two subtypes were evaluated. Clustering findings indicated that one prodromal PD subtype closely resembled the pattern of NCs (N-P; N = 159), while the other was similar to the pattern of PD (P-P; N = 103). Significant differences were observed between the subtypes in terms of multiple clinical measurements, neuroimaging for morphological changes, and gene enrichment for synaptic transmission. Identification of prodromal PD subtypes based on brain connectomes and a full understanding of heterogeneity at this phase could inform early and accurate PD diagnosis and effective neuroprotective interventions.

Deep-learning based Embedding of Functional Connectivity Profiles for Precision Functional Mapping

Spatial correlation of functional connectivity profiles across matching anatomical locations in individuals is often calculated to delineate individual differences in functional networks. Likewise, spatial correlation is assessed across average functional connectivity profiles of groups to evaluate the maturity of functional networks during development. Despite its widespread use, spatial correlation is limited to comparing two samples at a time. In this study, we employed a variational autoencoder to embed functional connectivity profiles from various anatomical locations, individuals, and group averages for simultaneous comparison. We demonstrate that our variational autoencoder, with pre-trained weights, can project new functional connectivity profiles from the vertex space to a latent space with as few as two dimensions, yet still retain meaningful global and local structures in the data. Functional connectivity profiles from various functional networks occupy distinct compartments of the latent space. Moreover, the variability of functional connectivity profiles from the same anatomical location is readily captured in the latent space. We believe that this approach could be useful for visualization and exploratory analyses in precision functional mapping.

Near-critical corticothalamic eigenmodes: Effects of nonuniform connectivity on modes, activity, and communication channels

The effects of nonuniformities in axonal connectivity on natural modes of brain activity are explored to determine their contributions to modal eigenvalues, structure, and communication and to clarify the limits of validity of widely used uniform-connectivity approximations. Preferred channels of communication are demonstrated that are supported by natural modes of mean connectivity and resulting activity. The effects of axonal tracts on these modes are calculated using perturbation methods, and it is found that modes and their spectra are only moderately perturbed by even the largest white matter tracts. However, perturbations of activity are greatly magnified when modes are near-critical and realistic connectivity and gain perturbations can then enable rapid responses to stimuli on the observed timescales of evoked responses. It is thus argued that dynamic mode-mode communication channels complement ones based on white matter tracts and that both rely on near-criticality to have their observed effects.

Sensory and motor cortices parcellations estimated via distance-weighted sparse representation with application to autism spectrum disorder

BACKGROUND

Motor impairments and sensory processing abnormalities are prevalent in autism spectrum disorder (ASD), closely related to the core functions of the primary motor cortex (M1) and the primary somatosensory cortex (S1). Currently, there is limited knowledge about potential therapeutic targets in the subregions of M1 and S1 in ASD patients. This study aims to map clinically significant functional subregions of M1 and S1.

METHODS

Resting-state functional magnetic resonance imaging data (NTD = 266) from Autism Brain Imaging Data Exchange (ABIDE) were used for subregion modeling. We proposed a distance-weighted sparse representation algorithm to construct brain functional networks. Functional subregions of M1 and S1 were identified through consensus clustering at the group level. Differences in the characteristics of functional subregions were analyzed, along with their correlation with clinical scores.

RESULTS

We observed symmetrical and continuous subregion organization from dorsal to ventral aspects in M1 and S1, with M1 subregions conforming to the functional pattern of the motor homunculus. Significant intergroup differences and clinical correlations were found in the dorsal and ventral aspects of M1 (p < 0.05/3, Bonferroni correction) and the ventromedial BA3 of S1 (p < 0.05/5). These functional characteristics were positively correlated with autism severity. All subregions showed significant results in the ROI-to-ROI intergroup differential analysis (p < 0.05/80).

LIMITATIONS

The generalizability of the segmentation model requires further evaluation.

CONCLUSIONS

This study highlights the significance of M1 and S1 in ASD treatment and may provide new insights into brain parcellation and the identification of therapeutic targets for ASD.

Multilayer Network Analysis across Cortical Depths in Resting-State 7T fMRI

In graph theory, "multilayer networks" represent systems involving several interconnected topological levels. One example in neuroscience is the stratification of connections between different cortical depths or "laminae", which is becoming non-invasively accessible in humans using ultra-high-resolution functional MRI (fMRI). Here, we applied multilayer graph theory to examine functional connectivity across different cortical depths in humans, using 7T fMRI (1-mm3 voxels; 30 participants). Blood oxygenation level dependent (BOLD) signals were derived from five depths between the white matter and pial surface. We compared networks where the inter-regional connections were limited to a single cortical depth only ("layer-by-layer matrices") to those considering all possible connections between areas and cortical depths ("multilayer matrix"). We utilized global and local graph theory features that quantitatively characterize network attributes including network composition, nodal centrality, path-based measures, and hub segregation. Detecting functional differences between cortical depths was improved using multilayer connectomics compared to the layer-by-layer versions. Superficial depths of the cortex dominated information transfer and deeper depths drove clustering. These differences were largest in frontotemporal and limbic regions. fMRI functional connectivity across different cortical depths may contain neurophysiologically relevant information; thus, multilayer connectomics could provide a methodological framework for studies on how information flows across this stratification.

Discovering prominent differences in structural and functional connectomes using a multinomial stochastic block model

Understanding the differences between functional and structural human brain connectivity has been a focus of an extensive amount of neuroscience research. We employ a novel approach using the multinomial stochastic block model (MSBM) to explicitly extract components that characterize prominent differences across graphs. We analyze structural and functional connectomes derived from high-resolution diffusion-weighted MRI and fMRI scans of 250 Human Connectome Project subjects, analyzed at group connectivity level across 50 subjects. The inferred brain partitions revealed consistent, spatially homogeneous clustering patterns across inferred resolutions demonstrating the MSBM's reliability in identifying brain areas with prominent structure-function differences. Prominent differences in low-resolution brain maps (K = {3, 4} clusters) were attributed to weak functional connectivity in the bilateral anterior temporal lobes, while higher resolution results (K ≥ 25) revealed stronger interhemispheric functional than structural connectivity. Our findings emphasize significant differences in high-resolution functional and structural connectomes, revealing challenges in extracting meaningful connectivity measurements from both modalities, including tracking fibers through the corpus callosum and attenuated functional connectivity in anterior temporal lobe fMRI data, which we attribute to increased noise levels. The MSBM emerges as a valuable tool for understanding differences across graphs, with potential future applications and avenues beyond the current focus on characterizing modality-specific distinctions in connectomics data.

Prediction of Post Traumatic Epilepsy Using MR-Based Imaging Markers

Post-traumatic epilepsy (PTE) is a debilitating neurological disorder that develops after traumatic brain injury (TBI). Despite the high prevalence of PTE, current methods for predicting its occurrence remain limited. In this study, we aimed to identify imaging-based markers for the prediction of PTE using machine learning. Specifically, we examined three imaging features: Lesion volumes, resting-state fMRI-based measures of functional connectivity, and amplitude of low-frequency fluctuation (ALFF). We employed three machine-learning methods, namely, kernel support vector machine (KSVM), random forest, and an artificial neural network (NN), to develop predictive models. Our results showed that the KSVM classifier, with all three feature types as input, achieved the best prediction accuracy of 0.78 AUC (area under the receiver operating characteristic (ROC) curve) using nested cross-validation. Furthermore, we performed voxel-wise and lobe-wise group difference analyses to investigate the specific brain regions and features that the model found to be most helpful in distinguishing PTE from non-PTE populations. Our statistical analysis uncovered significant differences in bilateral temporal lobes and cerebellum between PTE and non-PTE groups. Overall, our findings demonstrate the complementary prognostic value of MR-based markers in PTE prediction and provide new insights into the underlying structural and functional alterations associated with PTE.

Breakthroughs and challenges for generating brain network-based biomarkers of treatment response in depression

Treatment outcomes widely vary for individuals diagnosed with major depressive disorder, implicating a need for deeper understanding of the biological mechanisms conferring a greater likelihood of response to a particular treatment. Our improved understanding of intrinsic brain networks underlying depression psychopathology via magnetic resonance imaging and other neuroimaging modalities has helped reveal novel and potentially clinically meaningful biological markers of response. And while we have made considerable progress in identifying such biomarkers over the last decade, particularly with larger, multisite trials, there are significant methodological and practical obstacles that need to be overcome to translate these markers into the clinic. The aim of this review is to review current literature on brain network structural and functional biomarkers of treatment response or selection in depression, with a specific focus on recent large, multisite trials reporting predictive accuracy of candidate biomarkers. Regarding pharmaco- and psychotherapy, we discuss candidate biomarkers, reporting that while we have identified candidate biomarkers of response to a single intervention, we need more trials that distinguish biomarkers between first-line treatments. Further, we discuss the ways prognostic neuroimaging may help to improve treatment outcomes to neuromodulation-based therapies, such as transcranial magnetic stimulation and deep brain stimulation. Lastly, we highlight obstacles and technical developments that may help to address the knowledge gaps in this area of research. Ultimately, integrating neuroimaging-derived biomarkers into clinical practice holds promise for enhancing treatment outcomes and advancing precision psychiatry strategies for depression management. By elucidating the neural predictors of treatment response and selection, we can move towards more individualized and effective depression interventions, ultimately improving patient outcomes and quality of life.

A quantitative comparison of atlas parcellations on the human superior temporal sulcus

The superior temporal sulcus (STS) has a functional topography that has been difficult to characterize through traditional approaches. Automated atlas parcellations may be one solution while also being beneficial for both dimensional reduction and standardizing regions of interest, but they yield very different boundary definitions along the STS. Here we evaluate how well machine learning classifiers can correctly identify six social cognitive tasks from STS activation patterns dimensionally reduced using four popular atlases (Glasser et al., 2016; Gordon et al., 2016; Power et al., 2011 as projected onto the surface by Arslan et al., 2018; Schaefer et al., 2018). Functional data was summarized within each STS parcel in one of four ways, then subjected to leave-one-subject-out cross-validation SVM classification. We found that the classifiers could readily label conditions when data was parcellated using any of the four atlases, evidence that dimensional reduction to parcels did not compromise functional fingerprints. Mean activation for the social conditions was the most effective metric for classification in the right STS, whereas all the metrics classified equally well in the left STS. Interestingly, even atlases constructed from random parcellation schemes (null atlases) classified the conditions with high accuracy. We therefore conclude that the complex activation maps on the STS are readily differentiated at a coarse granular level, despite a strict topography having not yet been identified. Further work is required to identify what features have greatest potential to improve the utility of atlases in replacing functional localizers.

OpenMAP-T1: A Rapid Deep-Learning Approach to Parcellate 280 Anatomical Regions to Cover the Whole Brain

This study introduces OpenMAP-T1, a deep-learning-based method for rapid and accurate whole-brain parcellation in T1- weighted brain MRI, which aims to overcome the limitations of conventional normalization-to-atlas-based approaches and multi-atlas label-fusion (MALF) techniques. Brain image parcellation is a fundamental process in neuroscientific and clinical research, enabling a detailed analysis of specific cerebral regions. Normalization-to-atlas-based methods have been employed for this task, but they face limitations due to variations in brain morphology, especially in pathological conditions. The MALF techniques improved the accuracy of the image parcellation and robustness to variations in brain morphology, but at the cost of high computational demand that requires a lengthy processing time. OpenMAP-T1 integrates several convolutional neural network models across six phases: preprocessing; cropping; skull-stripping; parcellation; hemisphere segmentation; and final merging. This process involves standardizing MRI images, isolating the brain tissue, and parcellating it into 280 anatomical structures that cover the whole brain, including detailed gray and white matter structures, while simplifying the parcellation processes and incorporating robust training to handle various scan types and conditions. The OpenMAP-T1 was validated on the Johns Hopkins University atlas library and eight available open resources, including real-world clinical images, and the demonstration of robustness across different datasets with variations in scanner types, magnetic field strengths, and image processing techniques, such as defacing. Compared with existing methods, OpenMAP-T1 significantly reduced the processing time per image from several hours to less than 90 s without compromising accuracy. It was particularly effective in handling images with intensity inhomogeneity and varying head positions, conditions commonly seen in clinical settings. The adaptability of OpenMAP-T1 to a wide range of MRI datasets and its robustness to various scan conditions highlight its potential as a versatile tool in neuroimaging.

Direct segmentation of cortical cytoarchitectonic domains using ultra-high-resolution whole-brain diffusion MRI

We assess the potential of detecting cortical laminar patterns and areal borders by directly clustering voxel values of microstructural parameters derived from high-resolution mean apparent propagator (MAP) magnetic resonance imaging (MRI), as an alternative to conventional template-warping-based cortical parcellation methods. We acquired MAP-MRI data with 200μm resolution in a fixed macaque monkey brain. To improve the sensitivity to cortical layers, we processed the data with a local anisotropic Gaussian filter determined voxel-wise by the plane tangent to the cortical surface. We directly clustered all cortical voxels using only the MAP-derived microstructural imaging biomarkers, with no information regarding their relative spatial location or dominant diffusion orientations. MAP-based 3D cytoarchitectonic segmentation revealed laminar patterns similar to those observed in the corresponding histological images. Moreover, transition regions between these laminar patterns agreed more accurately with histology than the borders between cortical areas estimated using conventional atlas/template-warping cortical parcellation. By cross-tabulating all cortical labels in the atlas- and MAP-based segmentations, we automatically matched the corresponding MAP-derived clusters (i.e., cytoarchitectonic domains) across the left and right hemispheres. Our results demonstrate that high-resolution MAP-MRI biomarkers can effectively delineate three-dimensional cortical cytoarchitectonic domains in single individuals. Their intrinsic tissue microstructural contrasts enable the construction of whole-brain mesoscopic cortical atlases.

Genetic fingerprinting with heritable phenotypes of the resting-state brain network topology

Cognitive, behavioral, and disease traits are influenced by both genetic and environmental factors. Individual differences in these traits have been associated with graph theoretical properties of resting-state networks, indicating that variations in connectome topology may be driven by genetics. In this study, we establish the heritability of global and local graph properties of resting-state networks derived from functional MRI (fMRI) and magnetoencephalography (MEG) using a large sample of twins and non-twin siblings from the Human Connectome Project. We examine the heritability of MEG in the source space, providing a more accurate estimate of genetic influences on electrophysiological networks. Our findings show that most graph measures are more heritable for MEG compared to fMRI and the heritability for MEG is greater for amplitude compared to phase synchrony in the delta, high beta, and gamma frequency bands. This suggests that the fast neuronal dynamics in MEG offer unique insights into the genetic basis of brain network organization. Furthermore, we demonstrate that brain network features can serve as genetic fingerprints to accurately identify pairs of identical twins within a cohort. These results highlight novel opportunities to relate individual connectome signatures to genetic mechanisms underlying brain function.

A next-generation, histological atlas of the human brain and its application to automated brain MRI segmentation

Magnetic resonance imaging (MRI) is the standard tool to image the human brain in vivo. In this domain, digital brain atlases are essential for subject-specific segmentation of anatomical regions of interest (ROIs) and spatial comparison of neuroanatomy from different subjects in a common coordinate frame. High-resolution, digital atlases derived from histology (e.g., Allen atlas [7], BigBrain [13], Julich [15]), are currently the state of the art and provide exquisite 3D cytoarchitectural maps, but lack probabilistic labels throughout the whole brain. Here we present NextBrain, a next-generation probabilistic atlas of human brain anatomy built from serial 3D histology and corresponding highly granular delineations of five whole brain hemispheres. We developed AI techniques to align and reconstruct ~10,000 histological sections into coherent 3D volumes with joint geometric constraints (no overlap or gaps between sections), as well as to semi-automatically trace the boundaries of 333 distinct anatomical ROIs on all these sections. Comprehensive delineation on multiple cases enabled us to build the first probabilistic histological atlas of the whole human brain. Further, we created a companion Bayesian tool for automated segmentation of the 333 ROIs in any in vivo or ex vivo brain MRI scan using the NextBrain atlas. We showcase two applications of the atlas: automated segmentation of ultra-high-resolution ex vivo MRI and volumetric analysis of Alzheimer's disease and healthy brain ageing based on ~4,000 publicly available in vivo MRI scans. We publicly release: the raw and aligned data (including an online visualisation tool); the probabilistic atlas; the segmentation tool; and ground truth delineations for a 100 μm isotropic ex vivo hemisphere (that we use for quantitative evaluation of our segmentation method in this paper). By enabling researchers worldwide to analyse brain MRI scans at a superior level of granularity without manual effort or highly specific neuroanatomical knowledge, NextBrain holds promise to increase the specificity of MRI findings and ultimately accelerate our quest to understand the human brain in health and disease.

Age-related differences in resting-state, task-related, and structural brain connectivity: graph theoretical analyses and visual search performance

Previous magnetic resonance imaging (MRI) research suggests that aging is associated with a decrease in the functional interconnections within and between groups of locally organized brain regions (modules). Further, this age-related decrease in the segregation of modules appears to be more pronounced for a task, relative to a resting state, reflecting the integration of functional modules and attentional allocation necessary to support task performance. Here, using graph-theoretical analyses, we investigated age-related differences in a whole-brain measure of module connectivity, system segregation, for 68 healthy, community-dwelling individuals 18-78 years of age. We obtained resting-state, task-related (visual search), and structural (diffusion-weighted) MRI data. Using a parcellation of modules derived from the participants' resting-state functional MRI data, we demonstrated that the decrease in system segregation from rest to task (i.e., reconfiguration) increased with age, suggesting an age-related increase in the integration of modules required by the attentional demands of visual search. Structural system segregation increased with age, reflecting weaker connectivity both within and between modules. Functional and structural system segregation had qualitatively different influences on age-related decline in visual search performance. Functional system segregation (and reconfiguration) influenced age-related decline in the rate of visual evidence accumulation (drift rate), whereas structural system segregation contributed to age-related slowing of encoding and response processes (nondecision time). The age-related differences in the functional system segregation measures, however, were relatively independent of those associated with structural connectivity.

A brain subcortical segmentation tool based on anatomy attentional fusion network for developing macaques

Magnetic Resonance Imaging (MRI) plays a pivotal role in the accurate measurement of brain subcortical structures in macaques, which is crucial for unraveling the complexities of brain structure and function, thereby enhancing our understanding of neurodegenerative diseases and brain development. However, due to significant differences in brain size, structure, and imaging characteristics between humans and macaques, computational tools developed for human neuroimaging studies often encounter obstacles when applied to macaques. In this context, we propose an Anatomy Attentional Fusion Network (AAF-Net), which integrates multimodal MRI data with anatomical constraints in a multi-scale framework to address the challenges posed by the dynamic development, regional heterogeneity, and age-related size variations of the juvenile macaque brain, thus achieving precise subcortical segmentation. Specifically, we generate a Signed Distance Map (SDM) based on the initial rough segmentation of the subcortical region by a network as an anatomical constraint, providing comprehensive information on positions, structures, and morphology. Then we construct AAF-Net to fully fuse the SDM anatomical constraints and multimodal images for refined segmentation. To thoroughly evaluate the performance of our proposed tool, over 700 macaque MRIs from 19 datasets were used in this study. Specifically, we employed two manually labeled longitudinal macaque datasets to develop the tool and complete four-fold cross-validations. Furthermore, we incorporated various external datasets to demonstrate the proposed tool's generalization capabilities and promise in brain development research. We have made this tool available as an open-source resource at https://github.com/TaoZhong11/Macaque_subcortical_segmentation for direct application.

Deep multimodal saliency parcellation of cerebellar pathways: Linking microstructure and individual function through explainable multitask learning

Parcellation of human cerebellar pathways is essential for advancing our understanding of the human brain. Existing diffusion magnetic resonance imaging tractography parcellation methods have been successful in defining major cerebellar fibre tracts, while relying solely on fibre tract structure. However, each fibre tract may relay information related to multiple cognitive and motor functions of the cerebellum. Hence, it may be beneficial for parcellation to consider the potential importance of the fibre tracts for individual motor and cognitive functional performance measures. In this work, we propose a multimodal data-driven method for cerebellar pathway parcellation, which incorporates both measures of microstructure and connectivity, and measures of individual functional performance. Our method involves first training a multitask deep network to predict various cognitive and motor measures from a set of fibre tract structural features. The importance of each structural feature for predicting each functional measure is then computed, resulting in a set of structure-function saliency values that are clustered to parcellate cerebellar pathways. We refer to our method as Deep Multimodal Saliency Parcellation (DeepMSP), as it computes the saliency of structural measures for predicting cognitive and motor functional performance, with these saliencies being applied to the task of parcellation. Applying DeepMSP to a large-scale dataset from the Human Connectome Project Young Adult study (n = 1065), we found that it was feasible to identify multiple cerebellar pathway parcels with unique structure-function saliency patterns that were stable across training folds. We thoroughly experimented with all stages of the DeepMSP pipeline, including network selection, structure-function saliency representation, clustering algorithm, and cluster count. We found that a 1D convolutional neural network architecture and a transformer network architecture both performed comparably for the multitask prediction of endurance, strength, reading decoding, and vocabulary comprehension, with both architectures outperforming a fully connected network architecture. Quantitative experiments demonstrated that a proposed low-dimensional saliency representation with an explicit measure of motor versus cognitive category bias achieved the best parcellation results, while a parcel count of four was most successful according to standard cluster quality metrics. Our results suggested that motor and cognitive saliencies are distributed across the cerebellar white matter pathways. Inspection of the final k = 4 parcellation revealed that the highest-saliency parcel was most salient for the prediction of both motor and cognitive performance scores and included parts of the middle and superior cerebellar peduncles. Our proposed saliency-based parcellation framework, DeepMSP, enables multimodal, data-driven tractography parcellation. Through utilising both structural features and functional performance measures, this parcellation strategy may have the potential to enhance the study of structure-function relationships of the cerebellar pathways.

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.

Brain Network Interconnectivity Dynamics Explain Metacognitive Differences in Listening Behavior

Complex auditory scenes pose a challenge to attentive listening, rendering listeners slower and more uncertain in their perceptual decisions. How can we explain such behaviors from the dynamics of cortical networks that pertain to the control of listening behavior? We here follow up on the hypothesis that human adaptive perception in challenging listening situations is supported by modular reconfiguration of auditory-control networks in a sample of N = 40 participants (13 males) who underwent resting-state and task functional magnetic resonance imaging (fMRI). Individual titration of a spatial selective auditory attention task maintained an average accuracy of ∼70% but yielded considerable interindividual differences in listeners' response speed and reported confidence in their own perceptual decisions. Whole-brain network modularity increased from rest to task by reconfiguring auditory, cinguloopercular, and dorsal attention networks. Specifically, interconnectivity between the auditory network and cinguloopercular network decreased during the task relative to the resting state. Additionally, interconnectivity between the dorsal attention network and cinguloopercular network increased. These interconnectivity dynamics were predictive of individual differences in response confidence, the degree of which was more pronounced after incorrect judgments. Our findings uncover the behavioral relevance of functional cross talk between auditory and attentional-control networks during metacognitive assessment of one's own perception in challenging listening situations and suggest two functionally dissociable cortical networked systems that shape the considerable metacognitive differences between individuals in adaptive listening behavior.

Fast connectivity gradient approximation: maintaining spatially fine-grained connectivity gradients while reducing computational costs

Brain connectome analysis suffers from the high dimensionality of connectivity data, often forcing a reduced representation of the brain at a lower spatial resolution or parcellation. This is particularly true for graph-based representations, which are increasingly used to characterize connectivity gradients, capturing patterns of systematic spatial variation in the functional connectivity structure. However, maintaining a high spatial resolution is crucial for enabling fine-grained topographical analysis and preserving subtle individual differences that might otherwise be lost. Here we introduce a computationally efficient approach to establish spatially fine-grained connectivity gradients. At its core, it leverages a set of landmarks to approximate the underlying connectivity structure at the full spatial resolution without requiring a full-scale vertex-by-vertex connectivity matrix. We show that this approach reduces computational time and memory usage while preserving informative individual features and demonstrate its application in improving brain-behavior predictions. Overall, its efficiency can remove computational barriers and enable the widespread application of connectivity gradients to capture spatial signatures of the connectome. Importantly, maintaining a spatially fine-grained resolution facilitates to characterize the spatial transitions inherent in the core concept of gradients of brain organization.

Systematic evaluation of fMRI data-processing pipelines for consistent functional connectomics

Functional interactions between brain regions can be viewed as a network, enabling neuroscientists to investigate brain function through network science. Here, we systematically evaluate 768 data-processing pipelines for network reconstruction from resting-state functional MRI, evaluating the effect of brain parcellation, connectivity definition, and global signal regression. Our criteria seek pipelines that minimise motion confounds and spurious test-retest discrepancies of network topology, while being sensitive to both inter-subject differences and experimental effects of interest. We reveal vast and systematic variability across pipelines' suitability for functional connectomics. Inappropriate choice of data-processing pipeline can produce results that are not only misleading, but systematically so, with the majority of pipelines failing at least one criterion. However, a set of optimal pipelines consistently satisfy all criteria across different datasets, spanning minutes, weeks, and months. We provide a full breakdown of each pipeline's performance across criteria and datasets, to inform future best practices in functional connectomics.

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.

An end-to-end infant brain parcellation pipeline

Objective

Accurate infant brain parcellation is crucial for understanding early brain development; however, it is challenging due to the inherent low tissue contrast, high noise, and severe partial volume effects in infant magnetic resonance images (MRIs). The aim of this study was to develop an end-to-end pipeline that enabled accurate parcellation of infant brain MRIs.

Methods

We proposed an end-to-end pipeline that employs a two-stage global-to-local approach for accurate parcellation of infant brain MRIs. Specifically, in the global regions of interest (ROIs) localization stage, a combination of transformer and convolution operations was employed to capture both global spatial features and fine texture features, enabling an approximate localization of the ROIs across the whole brain. In the local ROIs refinement stage, leveraging the position priors from the first stage along with the raw MRIs, the boundaries o the ROIs are refined for a more accurate parcellation.

Results

We utilized the Dice ratio to evaluate the accuracy of parcellation results. Results on 263 subjects from National Database for Autism Research (NDAR), Baby Connectome Project (BCP) and Cross-site datasets demonstrated the better accuracy and robustness of our method than other competing methods.

Conclusion

Our end-to-end pipeline may be capable of accurately parcellating 6-month-old infant brain MRIs.

Spectral Graph Neural Network-Based Multi-Atlas Brain Network Fusion for Major Depressive Disorder Diagnosis

Major Depressive Disorder (MDD) imposes a substantial burden within the healthcare domain, impacting millions of individuals worldwide. Functional Magnetic Resonance Imaging (fMRI) has emerged as a promising tool for the objective diagnosis of MDD, enabling the investigation of functional connectivity patterns in the brain associated with this disorder. However, most existing methods focus on a single brain atlas, which limits their ability to capture the complex, multi-scale nature of functional brain networks. To address these limitations, we propose a novel multi-atlas fusion method that incorporates early and late fusion in a unified framework. Our method introduces the concept of the holistic Functional Connectivity Network (FCN), which captures both intra-atlas relationships within individual atlases and inter-regional relationships between atlases with different brain parcellation scales. This comprehensive representation enables the identification of potential disease-related patterns associated with MDD in the early stage of our framework. Moreover, by decoding the holistic FCN from various perspectives through multiple spectral Graph Convolutional Neural Networks and fusing their results with decision-level ensembles, we further improve the performance of MDD diagnosis. Our approach is easily implemented with minimal modifications to existing model structures and demonstrates a robust performance across different baseline models. Our method, evaluated on public resting-state fMRI datasets, surpasses the current multi-atlas fusion methods, enhancing the accuracy of MDD diagnosis. The proposed novel multi-atlas fusion framework provides a more reliable MDD diagnostic technique. Experimental results show our approach outperforms both single- and multi-atlas-based methods, demonstrating its effectiveness in advancing MDD diagnosis.

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

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

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

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

Abnormal higher-order network interactions in Parkinson's disease visual hallucinations

Visual hallucinations in Parkinson's disease can be viewed from a systems-level perspective, whereby dysfunctional communication between brain networks responsible for perception predisposes a person to hallucinate. To this end, abnormal functional interactions between higher-order and primary sensory networks have been implicated in the pathophysiology of visual hallucinations in Parkinson's disease, however the precise signatures remain to be determined. Dimensionality reduction techniques offer a novel means for simplifying the interpretation of multidimensional brain imaging data, identifying hierarchical patterns in the data that are driven by both within- and between-functional network changes. Here, we applied two complementary non-linear dimensionality reduction techniques-diffusion-map embedding and t-distributed stochastic neighbour embedding (t-SNE)-to resting state functional MRI data, in order to characterize the altered functional hierarchy associated with susceptibility to visual hallucinations. Our study involved 77 people with Parkinson's disease (31 with hallucinations; 46 without hallucinations) and 19 age-matched healthy control subjects. In patients with visual hallucinations, we found compression of the unimodal-heteromodal gradient consistent with increased functional integration between sensory and higher order networks. This was mirrored in a traditional functional connectivity analysis, which showed increased connectivity between the visual and default mode networks in the hallucinating group. Together, these results suggest a route by which higher-order regions may have excessive influence over earlier sensory processes, as proposed by theoretical models of hallucinations across disorders. By contrast, the t-SNE analysis identified distinct alterations in prefrontal regions, suggesting an additional layer of complexity in the functional brain network abnormalities implicated in hallucinations, which was not apparent in traditional functional connectivity analyses. Together, the results confirm abnormal brain organization associated with the hallucinating phenotype in Parkinson's disease and highlight the utility of applying convergent dimensionality reduction techniques to investigate complex clinical symptoms. In addition, the patterns we describe in Parkinson's disease converge with those seen in other conditions, suggesting that reduced hierarchical differentiation across sensory-perceptual systems may be a common transdiagnostic vulnerability in neuropsychiatric disorders with perceptual disturbances.

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.

Multivariate patterns of brain functional connectome associated with COVID-19-related negative affect symptoms

Severe mental health problems with the representation of negative affect symptoms (NAS) have been increasingly reported during the coronavirus disease 2019 (COVID-19) pandemic. This study aimed to explore the multivariate patterns of brain functional connectome predicting COVID-19-related NAS. This cohort study encompassed a group of university students to undergo neuroimaging scans before the pandemic, and we re-contacted participants for 1-year follow-up COVID-related NAS evaluations during the pandemic. Regularized canonical correlation analysis was used to identify connectome-based dimensions of NAS to compute pairs of canonical variates. The predictive ability of identified functional connectome to NAS dimensional scores was examined with a nested cross-validation. Two dimensions (i.e. mode stress and mode anxiety) were related to distinct patterns of brain functional connectome (r2 = 0.911, PFDR = 0.048; r2 = 0.901, PFDR = 0.037, respectively). Mode anxiety was characterized by high loadings in connectivity between affective network (AFN) and visual network (VN), while connectivity of the default mode network with dorsal attention network (DAN) were remarkably prominent in mode stress. Connectivity patterns within the DAN and between DAN and VN, ventral attention network, and AFN was common for both dimensions. The identified functional connectome can reliably predict mode stress (r = 0.37, MAE = 5.1, p < 0.001) and mode anxiety (r = 0.28, MAE = 5.4, p = 0.005) in the cross-validation. Our findings provide new insight into multivariate dimensions of COVID-related NAS, which may have implications for developing network-based biomarkers in psychological interventions for vulnerable individuals in the pandemic.

Graph Learning for Cortical Parcellation from Tensor Decompositions of Resting-State fMRI

Cortical parcellation has long been a cornerstone in the field of neuroscience, enabling the cerebral cortex to be partitioned into distinct, non-overlapping regions that facilitate the interpretation and comparison of complex neuroscientific data. In recent years, these parcellations have frequently been based on the use of resting-state fMRI (rsfMRI) data. In parallel, methods such as independent components analysis have long been used to identify large-scale functional networks with significant spatial overlap between networks. Despite the fact that both forms of decomposition make use of the same spontaneous brain activity measured with rsfMRI, a gap persists in establishing a clear relationship between disjoint cortical parcellations and brain-wide networks. To address this, we introduce a novel parcellation framework that integrates NASCAR, a three-dimensional tensor decomposition method that identifies a series of functional brain networks, with state-of-the-art graph representation learning to produce cortical parcellations that represent near-homogeneous functional regions that are consistent with these brain networks. Further, through the use of the tensor decomposition, we avoid the limitations of traditional approaches that assume statistical independence or orthogonality in defining the underlying networks. Our findings demonstrate that these parcellations are comparable or superior to established atlases in terms of homogeneity of the functional connectivity across parcels, task contrast alignment, and architectonic map alignment. Our methodological pipeline is highly automated, allowing for rapid adaptation to new datasets and the generation of custom parcellations in just minutes, a significant advancement over methods that require extensive manual input. We describe this integrated approach, which we refer to as Untamed, as a tool for use in the fields of cognitive and clinical neuroscientific research. Parcellations created from the Human Connectome Project dataset using Untamed, along with the code to generate atlases with custom parcel numbers, are publicly available at https://untamed-atlas.github.io.

Effects of PRRT2 mutation on brain gray matter networks in paroxysmal kinesigenic dyskinesia

Although proline-rich transmembrane protein 2 is the primary causative gene of paroxysmal kinesigenic dyskinesia, its effects on the brain structure of paroxysmal kinesigenic dyskinesia patients are not yet clear. Here, we explored the influence of proline-rich transmembrane protein 2 mutations on similarity-based gray matter morphological networks in individuals with paroxysmal kinesigenic dyskinesia. A total of 51 paroxysmal kinesigenic dyskinesia patients possessing proline-rich transmembrane protein 2 mutations, 55 paroxysmal kinesigenic dyskinesia patients possessing proline-rich transmembrane protein 2 non-mutation, and 80 healthy controls participated in the study. We analyzed the structural connectome characteristics across groups by graph theory approaches. Relative to paroxysmal kinesigenic dyskinesia patients possessing proline-rich transmembrane protein 2 non-mutation and healthy controls, paroxysmal kinesigenic dyskinesia patients possessing proline-rich transmembrane protein 2 mutations exhibited a notable increase in characteristic path length and a reduction in both global and local efficiency. Relative to healthy controls, both patient groups showed reduced nodal metrics in right postcentral gyrus, right angular, and bilateral thalamus; Relative to healthy controls and paroxysmal kinesigenic dyskinesia patients possessing proline-rich transmembrane protein 2 non-mutation, paroxysmal kinesigenic dyskinesia patients possessing proline-rich transmembrane protein 2 mutations showed almost all reduced nodal centralities and structural connections in cortico-basal ganglia-thalamo-cortical circuit including bilateral supplementary motor area, bilateral pallidum, and right caudate nucleus. Finally, we used support vector machine by gray matter network matrices to classify paroxysmal kinesigenic dyskinesia patients possessing proline-rich transmembrane protein 2 mutations and paroxysmal kinesigenic dyskinesia patients possessing proline-rich transmembrane protein 2 non-mutation, achieving an accuracy of 73%. These results show that proline-rich transmembrane protein 2 related gray matter network deficits may contribute to paroxysmal kinesigenic dyskinesia, offering new insights into its pathophysiological mechanisms.

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.

Causal functional connectivity in Alzheimer's disease computed from time series fMRI data

Functional connectivity between brain regions is known to be altered in Alzheimer's disease and promises to be a biomarker for early diagnosis. Several approaches for functional connectivity obtain an un-directed network representing stochastic associations (correlations) between brain regions. However, association does not necessarily imply causation. In contrast, Causal Functional Connectivity (CFC) is more informative, providing a directed network representing causal relationships between brain regions. In this paper, we obtained the causal functional connectome for the whole brain from resting-state functional magnetic resonance imaging (rs-fMRI) recordings of subjects from three clinical groups: cognitively normal, mild cognitive impairment, and Alzheimer's disease. We applied the recently developed Time-aware PC (TPC) algorithm to infer the causal functional connectome for the whole brain. TPC supports model-free estimation of whole brain CFC based on directed graphical modeling in a time series setting. We compared the CFC outcome of TPC with that of other related approaches in the literature. Then, we used the CFC outcomes of TPC and performed an exploratory analysis of the difference in strengths of CFC edges between Alzheimer's and cognitively normal groups, based on edge-wise p-values obtained by Welch's t-test. The brain regions thus identified are found to be in agreement with literature on brain regions impacted by Alzheimer's disease, published by researchers from clinical/medical institutions.

Global Alterations of Whole Brain Structural Connectome in Parkinson's Disease: A Meta-analysis

Recent graph-theoretical studies of Parkinson's disease (PD) have examined alterations in the global properties of the brain structural connectome; however, reported alterations are not consistent. The present study aimed to identify the most robust global metric alterations in PD via a meta-analysis. A comprehensive literature search was conducted for all available diffusion MRI structural connectome studies that compared global graph metrics between PD patients and healthy controls (HC). Hedges' g effect sizes were calculated for each study and then pooled using a random-effects model in Comprehensive Meta-Analysis software, and the effects of potential moderator variables were tested. A total of 22 studies met the inclusion criteria for review. Of these, 16 studies reporting 10 global graph metrics (916 PD patients; 560 HC) were included in the meta-analysis. In the structural connectome of PD patients compared with HC, we found a significant decrease in clustering coefficient (g = -0.357, P = 0.005) and global efficiency (g = -0.359, P < 0.001), and a significant increase in characteristic path length (g = 0.250, P = 0.006). Dopaminergic medication, sex and age of patients were potential moderators of global brain network changes in PD. These findings provide evidence of decreased global segregation and integration of the structural connectome in PD, indicating a shift from a balanced small-world network to 'weaker small-worldization', which may provide useful markers of the pathophysiological mechanisms underlying PD.

Functional parcellation of the neonatal brain

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 adult- and older infant-derived parcels are a poor fit with neonatal data, emphasizing the need for neonatal-specific parcels. 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.

A personalized cortical atlas for functional regions of interest

Advances in functional MRI (fMRI) allow mapping an individual's brain function in vivo. Task fMRI can localize domain-specific regions of cognitive processing or functional regions of interest (fROIs) within an individual. Moreover, data from resting state (no task) fMRI can be used to define an individual's connectome, which can characterize that individual's functional organization via connectivity-based parcellations. However, can connectivity-based parcellations alone predict an individual's fROIs? Here, we describe an approach to compute individualized rs-fROIs (i.e., regions that correspond to given fROI constructed using only resting state data) for motor control, working memory, high-level vision, and language comprehension. The rs-fROIs were computed and validated using a large sample of young adults (n = 1,018) with resting state and task fMRI from the Human Connectome Project. First, resting state parcellations were defined across a sequence of resolutions from broadscale to fine-grained networks in a training group of 500 individuals. Second, 21 rs-fROIs were defined from the training group by identifying the rs network that most closely matched task-defined fROIs across all individuals. Third, the selectivity of rs-fROIs was investigated in a training set of the remaining 518 individuals. All computed rs-fROIs were indeed selective for their preferred category. Critically, the rs-fROIs had higher selectivity than probabilistic atlas parcels for nearly all fROIs. In conclusion, we present a potential approach to define selective fROIs on an individual-level circumventing the need for multiple task-based localizers.NEW & NOTEWORTHY We compute individualized resting state parcels that identify an individual's own functional regions of interest (fROIs) for high-level vision, language comprehension, motor control, and working memory, using only their functional connectome. This approach demonstrates a rapid and powerful alternative for finding a large set of fROIs in an individual, using only their unique connectivity pattern, which does not require the costly acquisition of multiple fMRI localizer tasks.

Hippocampal metabolic subregions and networks: Behavioral, molecular, and pathological aging profiles

INTRODUCTION

Hippocampal local and network dysfunction is the hallmark of Alzheimer's disease (AD).

METHODS

We characterized the spatial patterns of hippocampus differentiation based on brain co-metabolism in healthy elderly participants and demonstrated their relevance to study local metabolic changes and associated dysfunction in pathological aging.

RESULTS

The hippocampus can be differentiated into anterior/posterior and dorsal cornu ammonis (CA)/ventral (subiculum) subregions. While anterior/posterior CA show co-metabolism with different regions of the subcortical limbic networks, the anterior/posterior subiculum are parts of cortical networks supporting object-centered memory and higher cognitive demands, respectively. Both networks show relationships with the spatial patterns of gene expression pertaining to cell energy metabolism and AD's process. Finally, while local metabolism is generally lower in posterior regions, the anterior-posterior imbalance is maximal in late mild cognitive impairment with the anterior subiculum being relatively preserved.

DISCUSSION

Future studies should consider bidimensional hippocampal differentiation and in particular the posterior subicular region to better understand pathological aging.

Confounding Effects on the Performance of Machine Learning Analysis of Static Functional Connectivity Computed from rs-fMRI Multi-site Data

Resting-state functional magnetic resonance imaging (rs-fMRI) is a non-invasive imaging technique widely used in neuroscience to understand the functional connectivity of the human brain. While rs-fMRI multi-site data can help to understand the inner working of the brain, the data acquisition and processing of this data has many challenges. One of the challenges is the variability of the data associated with different acquisitions sites, and different MRI machines vendors. Other factors such as population heterogeneity among different sites, with variables such as age and gender of the subjects, must also be considered. Given that most of the machine-learning models are developed using these rs-fMRI multi-site data sets, the intrinsic confounding effects can adversely affect the generalizability and reliability of these computational methods, as well as the imposition of upper limits on the classification scores. This work aims to identify the phenotypic and imaging variables producing the confounding effects, as well as to control these effects. Our goal is to maximize the classification scores obtained from the machine learning analysis of the Autism Brain Imaging Data Exchange (ABIDE) rs-fMRI multi-site data. To achieve this goal, we propose novel methods of stratification to produce homogeneous sub-samples of the 17 ABIDE sites, as well as the generation of new features from the static functional connectivity values, using multiple linear regression models, ComBat harmonization models, and normalization methods. The main results obtained with our statistical models and methods are an accuracy of 76.4%, sensitivity of 82.9%, and specificity of 77.0%, which are 8.8%, 20.5%, and 7.5% above the baseline classification scores obtained from the machine learning analysis of the static functional connectivity computed from the ABIDE rs-fMRI multi-site data.

Effect of Extremely Preterm Birth on Adolescent Brain Network Organization

Introduction: Extremely preterm (EPT) birth, defined as birth at a gestational age (GA) <28 weeks, can have a lasting impact on cognition throughout the life span. Previous investigations reveal differences in brain structure and connectivity between infants born preterm and full-term (FT), but how does preterm birth impact the adolescent connectome? Methods: In this study, we investigate how EPT birth can alter broadscale network organization later in life by comparing resting-state functional magnetic resonance imaging connectome-based parcellations of the entire cortex in adolescents born EPT (N = 22) to age-matched adolescents born FT (GA ≥37 weeks, N = 28). We compare these parcellations to adult parcellations from previous studies and explore the relationship between an individual's network organization and behavior. Results: Primary (occipital and sensorimotor) and frontoparietal networks were observed in both groups. However, there existed notable differences in the limbic and insular networks. Surprisingly, the connectivity profile of the limbic network of EPT adolescents was more adultlike than the same network in FT adolescents. Finally, we found a relationship between adolescents' overall cognition score and their limbic network maturity. Discussion: Overall, preterm birth may contribute to the atypical development of broadscale network organization in adolescence and may partially explain the observed cognitive deficits.