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

Review of weighted exponential random graph models frameworks applied to neuroimaging

Neuro-imaging data can often be represented as statistical networks, especially for functional magnetic resonance imaging (fMRI) data, where brain regions are defined as nodes and the functional interactions between those regions are taken as edges. Such networks are commonly divided into classes depending on the type of edges, namely binary or weighted. A binary network means edges can either be present or absent. Whereas the edges of a weighted network are associated with weight values, and fMRI networks belong to weighted networks. Statistical methods are often adopted to analyse such networks, among which, the exponential random graph model (ERGM) is an important network analysis approach. Typically ERGMs are applied to binary networks, and weighted networks often need to be binarised by arbitrarily selecting a threshold value to define the presence of the edges, which can lead to non-robustness and loss of valuable edge weight information representing the strength of fMRI interaction in fMRI networks. While it is therefore important to gain deeper insight in adopting ERGM on weighted networks, there only exists a few different ERGM frameworks for weighted networks; some of these are not directly implementable on fMRI networks based on their original proposal. We systematically review, implement, analyse and compare five such frameworks via a simulation study and provide guidelines on each modelling framework as well as conclude the suitability of them on fMRI networks based on a range of criteria. We concluded that Multi-Layered ERGM is currently the most suitable framework.

fMRI-based spatio-temporal parcellations of the human brain

PURPOSE OF REVIEW

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

RECENT FINDINGS

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

SUMMARY

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

The Transdiagnostic Connectome Project: a richly phenotyped open dataset for advancing the study of brain-behavior relationships in psychiatry

An important aim in psychiatry is the establishment of valid and reliable associations linking profiles of brain functioning to clinically relevant symptoms and behaviors across patient populations. To advance progress in this area, we introduce an open dataset containing behavioral and neuroimaging data from 241 individuals aged 18 to 70, comprising 148 individuals meeting diagnostic criteria for a broad range of psychiatric illnesses and a healthy comparison group of 93 individuals. These data include high-resolution anatomical scans, multiple resting-state, and task-based functional MRI runs. Additionally, participants completed over 50 psychological and cognitive assessments. Here, we detail available behavioral data as well as raw and processed MRI derivatives. Associations between data processing and quality metrics, such as head motion, are reported. Processed data exhibit classic task activation effects and canonical functional network organization. Overall, we provide a comprehensive and analysis-ready transdiagnostic dataset, which we hope will accelerate the identification of illness-relevant features of brain functioning, enabling future discoveries in basic and clinical neuroscience.

In vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here, we noninvasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the gamma-aminobutyric acid (GABA) agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in the association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 y old) and Asian (7.2 to 7.9 y old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

Altered asymmetry of functional connectome gradients in major depressive disorder

Introduction

Major depressive disorder (MDD) is a debilitating disease involving sensory and higher-order cognitive dysfunction. Previous work has shown altered asymmetry in MDD, including abnormal lateralized activation and disrupted hemispheric connectivity. However, it remains unclear whether and how MDD affects functional asymmetries in the context of intrinsic hierarchical organization.

Methods

Here, we evaluate intra- and inter-hemispheric asymmetries of the first three functional gradients, characterizing unimodal-transmodal, visual-somatosensory, and somatomotor/default mode-multiple demand hierarchies, to study MDD-related alterations in overarching system-level architecture.

Results

We find that, relative to the healthy controls, MDD patients exhibit alterations in both primary sensory regions (e.g., visual areas) and transmodal association regions (e.g., default mode areas). We further find these abnormalities are woven in heterogeneous alterations along multiple functional gradients, associated with cognitive terms involving mind, memory, and visual processing. Moreover, through an elastic net model, we observe that both intra- and inter-asymmetric features are predictive of depressive traits measured by BDI-II scores.

Discussion

Altogether, these findings highlight a broad and mixed effect of MDD on functional gradient asymmetry, contributing to a richer understanding of the neurobiological underpinnings in MDD.

Long-range inhibitory neurons mediate cortical neurovascular coupling

To meet the high energy demands of brain function, cerebral blood flow (CBF) parallels changes in neuronal activity by a mechanism known as neurovascular coupling (NVC). However, which neurons play a role in mediating NVC is not well understood. Here, we identify in mice and humans a specific population of cortical GABAergic neurons that co-express neuronal nitric oxide synthase and tachykinin receptor 1 (Tacr1). Through whole-tissue clearing, we demonstrate that Tacr1 neurons extend local and long-range projections across functionally connected cortical areas. We show that whisker stimulation elicited Tacr1 neuron activity in the barrel cortex through feedforward excitatory pathways. Additionally, through optogenetic experiments, we demonstrate that Tacr1 neurons are instrumental in mediating CBF through the relaxation of mural cells in a similar fashion to whisker stimulation. Finally, by electron microscopy, we observe that Tacr1 processes contact astrocytic endfeet. These findings suggest that Tacr1 neurons integrate cortical activity to mediate NVC.

Using a deep generation network reveals neuroanatomical specificity in hemispheres

Asymmetry is an important property of brain organization, but its nature is still poorly understood. Capturing the neuroanatomical components specific to each hemisphere facilitates the understanding of the establishment of brain asymmetry. Since deep generative networks (DGNs) have powerful inference and recovery capabilities, we use one hemisphere to predict the opposite hemisphere by training the DGNs, which automatically fit the built-in dependencies between the left and right hemispheres. After training, the reconstructed images approximate the homologous components in the hemisphere. We use the difference between the actual and reconstructed hemispheres to measure hemisphere-specific components due to asymmetric expression of environmental and genetic factors. The results show that our model is biologically plausible and that our proposed metric of hemispheric specialization is reliable, representing a wide range of individual variation. Together, this work provides promising tools for exploring brain asymmetry and new insights into self-supervised DGNs for representing the brain.

In-vivo whole-cortex marker of excitation-inhibition ratio indexes cortical maturation and cognitive ability in youth

A balanced excitation-inhibition ratio (E/I ratio) is critical for healthy brain function. Normative development of cortex-wide E/I ratio remains unknown. Here we non-invasively estimate a putative marker of whole-cortex E/I ratio by fitting a large-scale biophysically-plausible circuit model to resting-state functional MRI (fMRI) data. We first confirm that our model generates realistic brain dynamics in the Human Connectome Project. Next, we show that the estimated E/I ratio marker is sensitive to the GABA-agonist benzodiazepine alprazolam during fMRI. Alprazolam-induced E/I changes are spatially consistent with positron emission tomography measurement of benzodiazepine receptor density. We then investigate the relationship between the E/I ratio marker and neurodevelopment. We find that the E/I ratio marker declines heterogeneously across the cerebral cortex during youth, with the greatest reduction occurring in sensorimotor systems relative to association systems. Importantly, among children with the same chronological age, a lower E/I ratio marker (especially in association cortex) is linked to better cognitive performance. This result is replicated across North American (8.2 to 23.0 years old) and Asian (7.2 to 7.9 years old) cohorts, suggesting that a more mature E/I ratio indexes improved cognition during normative development. Overall, our findings open the door to studying how disrupted E/I trajectories may lead to cognitive dysfunction in psychopathology that emerges during youth.

Age-related decline in cognitive flexibility and inadequate preparation: evidence from task-state network analysis

Behavioral evidence showed decreased cognitive flexibility in older adults. However, task-based network mechanisms of cognitive flexibility in aging (CFA) remain unclear. Here, we provided the first task-state network evidence that CFA was associated with inadequate preparation for switching trials by revealing age-related changes in functional integration. We examined functional integration in a letter-number switch task that distinguished between the cue and target stages. Both young and older adults showed decreased functional integration from the cue stage to the target stage, indicating that control-related processes were executed as the task progressed. However, compared to young adults, older adults showed less cue-to-target reduction in functional integration, which was primarily driven by higher network integration in the target stage. Moreover, less cue-to-target reductions were correlated with age-related decreases in task performance in the switch task. To sum up, compared to young adults, older adults pre-executed less control-related processes in the cue stage and more control-related processes in the target stage. Therefore, the decline in cognitive flexibility in older adults was associated with inadequate preparation for the impending demands of cognitive switching. This study offered novel insights into network mechanisms underlying CFA. Furthermore, we highlighted that training the function of brain networks, in conjunction with providing more preparation time for older adults, may be beneficial to their cognitive flexibility.

Whole-brain structural connectome asymmetry in autism

Autism spectrum disorder is a common neurodevelopmental condition that manifests as a disruption in sensory and social skills. Although it has been shown that the brain morphology of individuals with autism is asymmetric, how this differentially affects the structural connectome organization of each hemisphere remains under-investigated. We studied whole-brain structural connectivity-based brain asymmetry in individuals with autism using diffusion magnetic resonance imaging obtained from the Autism Brain Imaging Data Exchange initiative. By leveraging dimensionality reduction techniques, we constructed low-dimensional representations of structural connectivity and calculated their asymmetry index. Comparing the asymmetry index between individuals with autism and neurotypical controls, we found atypical structural connectome asymmetry in the sensory and default-mode regions, particularly showing weaker asymmetry towards the right hemisphere in autism. Network communication provided topological underpinnings by demonstrating that the inferior temporal cortex and limbic and frontoparietal regions showed reduced global network communication efficiency and decreased send-receive network navigation in the inferior temporal and lateral visual cortices in individuals with autism. Finally, supervised machine learning revealed that structural connectome asymmetry could be used as a measure for predicting communication-related autistic symptoms and nonverbal intelligence. Our findings provide insights into macroscale structural connectome alterations in autism and their topological underpinnings.

Individual variability in functional connectivity of human auditory cortex

Individual variability in functional connectivity underlies individual differences in cognition and behaviors, yet its association with functional specialization in the auditory cortex remains elusive. Using resting-state functional magnetic resonance imaging data from the Human Connectome Project, this study was designed to investigate the spatial distribution of auditory cortex individual variability in its whole-brain functional network architecture. An inherent hierarchical axis of the variability was discerned, which radiates from the medial to lateral orientation, with the left auditory cortex demonstrating more pronounced variations than the right. This variability exhibited a significant correlation with the variations in structural and functional metrics in the auditory cortex. Four auditory cortex subregions, which were identified from a clustering analysis based on this variability, exhibited unique connectional fingerprints and cognitive maps, with certain subregions showing specificity to speech perception functional activation. Moreover, the lateralization of the connectional fingerprint exhibited a U-shaped trajectory across the subregions. These findings emphasize the role of individual variability in functional connectivity in understanding cortical functional organization, as well as in revealing its association with functional specialization from the activation, connectome, and cognition perspectives.

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.

Modular brain network in volitional eyes closing: enhanced integration with a marked impact on hubs

Volitional eyes closing would shift brain's information processing modes from the "exteroceptive" to "interoceptive" state. This transition induced by the eyes closing is underpinned by a large-scale reconfiguration of brain network, which is still not fully comprehended. Here, we investigated the eyes-closing-relevant network reconfiguration by examining the functional integration among intrinsic modules. Our investigation utilized a publicly available dataset with 48 subjects being scanned in both eyes closed and eyes open conditions. It was found that the modular integration was significantly enhanced during the eyes closing, including lower modularity index, higher participation coefficient, less provincial hubs, and more connector hubs. Moreover, the eyes-closing-enhanced integration was particularly noticeable in the hubs of network, mainly located in the default-mode network. Finally, the hub-dominant modular enhancement was positively correlated to the eyes-closing-reduced entropy of BOLD signal, suggesting a close connection to the diminished consciousness of individuals. Collectively, our findings strongly suggested that the enhanced modular integration with substantially reorganized hubs characterized the large-scale cortical underpinning of the volitional eyes closing.

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.

Evidence for a Compensatory Relationship between Left- and Right-Lateralized Brain Networks

The two hemispheres of the human brain are functionally asymmetric. At the network level, the language network exhibits left-hemisphere lateralization. While this asymmetry is widely replicated, the extent to which other functional networks demonstrate lateralization remains a subject of Investigation. Additionally, it is unknown how the lateralization of one functional network may affect the lateralization of other networks within individuals. We quantified lateralization for each of 17 networks by computing the relative surface area on the left and right cerebral hemispheres. After examining the ecological, convergent, and external validity and test-retest reliability of this surface area-based measure of lateralization, we addressed two hypotheses across multiple datasets (Human Connectome Project = 553, Human Connectome Project-Development = 343, Natural Scenes Dataset = 8). First, we hypothesized that networks associated with language, visuospatial attention, and executive control would show the greatest lateralization. Second, we hypothesized that relationships between lateralized networks would follow a dependent relationship such that greater left-lateralization of a network would be associated with greater right-lateralization of a different network within individuals, and that this pattern would be systematic across individuals. A language network was among the three networks identified as being significantly left-lateralized, and attention and executive control networks were among the five networks identified as being significantly right-lateralized. Next, correlation matrices, an exploratory factor analysis, and confirmatory factor analyses were used to test the second hypothesis and examine the organization of lateralized networks. We found general support for a dependent relationship between highly left- and right-lateralized networks, meaning that across subjects, greater left lateralization of a given network (such as a language network) was linked to greater right lateralization of another network (such as a ventral attention/salience network) and vice versa. These results further our understanding of brain organization at the macro-scale network level in individuals, carrying specific relevance for neurodevelopmental conditions characterized by disruptions in lateralization such as autism and schizophrenia.

Brain decoding of the Human Connectome Project tasks in a dense individual fMRI dataset

Brain decoding aims to infer cognitive states from patterns of brain activity. Substantial inter-individual variations in functional brain organization challenge accurate decoding performed at the group level. In this paper, we tested whether accurate brain decoding models can be trained entirely at the individual level. We trained several classifiers on a dense individual functional magnetic resonance imaging (fMRI) dataset for which six participants completed the entire Human Connectome Project (HCP) task battery >13 times over ten separate fMRI sessions. We evaluated nine decoding methods, from Support Vector Machines (SVM) and Multi-Layer Perceptron (MLP) to Graph Convolutional Neural Networks (GCN). All decoders were trained to classify single fMRI volumes into 21 experimental conditions simultaneously, using ∼7 h of fMRI data per participant. The best prediction accuracies were achieved with GCN and MLP models, whose performance (57-67 % accuracy) approached state-of-the-art accuracy (76 %) with models trained at the group level on >1 K hours of data from the original HCP sample. Our SVM model also performed very well (54-62 % accuracy). Feature importance maps derived from MLP -our best-performing model- revealed informative features in regions relevant to particular cognitive domains, notably in the motor cortex. We also observed that inter-subject classification achieved substantially lower accuracy than subject-specific models, indicating that our decoders learned individual-specific features. This work demonstrates that densely-sampled neuroimaging datasets can be used to train accurate brain decoding models at the individual level. We expect this work to become a useful benchmark for techniques that improve model generalization across multiple subjects and acquisition conditions.

CAMBA framework: Unveiling the brain asymmetry alterations and longitudinal changes after stroke using resting-state EEG

Hemispheric asymmetry or lateralization is a fundamental principle of brain organization. However, it is poorly understood to what extent the brain asymmetries across different levels of functional organizations are evident in health or altered in brain diseases. Here, we propose a framework that integrates three degrees of brain interactions (isolated nodes, node-node, and edge-edge) into a unified analysis pipeline to capture the sliding window-based asymmetry dynamics at both the node and hemisphere levels. We apply this framework to resting-state EEG in healthy and stroke populations and investigate the stroke-induced abnormal alterations in brain asymmetries and longitudinal asymmetry changes during poststroke rehabilitation. We observe that the mean asymmetry in patients was abnormally enhanced across different frequency bands and levels of brain interactions, with these abnormal patterns strongly associated with the side of the stroke lesion. Compared to healthy controls, patients displayed significant alterations in asymmetry fluctuations, disrupting and reconfiguring the balance of inter-hemispheric integration and segregation. Additionally, analyses reveal that specific abnormal asymmetry metrics in patients tend to move towards those observed in healthy controls after short-term brain-computer interface rehabilitation. Furthermore, preliminary evidence suggests that baseline clinical and asymmetry features can predict poststroke improvements in the Fugl-Meyer assessment of the lower extremity (mean absolute error of about 2). Overall, these findings advance our understanding of hemispheric asymmetry. Our framework offers new insights into the mechanisms underlying brain alterations and recovery after a brain lesion, may help identify prognostic biomarkers, and can be easily extended to different functional modalities.

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.

The regional variation of laminar thickness in the human isocortex is related to cortical hierarchy and interregional connectivity

The human isocortex consists of tangentially organized layers with unique cytoarchitectural properties. These layers show spatial variations in thickness and cytoarchitecture across the neocortex, which is thought to support function through enabling targeted corticocortical connections. Here, leveraging maps of the 6 cortical layers based on 3D human brain histology, we aimed to quantitatively characterize the systematic covariation of laminar structure in the cortex and its functional consequences. After correcting for the effect of cortical curvature, we identified a spatial pattern of changes in laminar thickness covariance from lateral frontal to posterior occipital regions, which differentiated the dominance of infra- versus supragranular layer thickness. Corresponding to the laminar regularities of cortical connections along cortical hierarchy, the infragranular-dominant pattern of laminar thickness was associated with higher hierarchical positions of regions, mapped based on resting-state effective connectivity in humans and tract-tracing of structural connections in macaques. Moreover, we show that regions with similar laminar thickness patterns have a higher likelihood of structural connections and strength of functional connections. In sum, here we characterize the organization of laminar thickness in the human isocortex and its association with cortico-cortical connectivity, illustrating how laminar organization may provide a foundational principle of cortical function.