MSM: a new flexible framework for Multimodal Surface Matching

JOSA: Joint surface-based registration and atlas construction of brain geometry and function

Surface-based cortical registration is an important topic in medical image analysis and facilitates many downstream applications. Current approaches for cortical registration are mainly driven by geometric features, such as sulcal depth and curvature, and often assume that registration of folding patterns leads to alignment of brain function. However, functional variability of anatomically corresponding areas across subjects has been widely reported, particularly in higher-order cognitive areas. In this work, we present JOSA, a novel cortical registration framework that jointly models the mismatch between geometry and function while simultaneously learning an unbiased population-specific atlas. Using a semi-supervised training strategy, JOSA achieves superior registration performance in both geometry and function to the state-of-the-art methods but without requiring functional data at inference. This learning framework can be extended to any auxiliary data to guide spherical registration that is available during training but is difficult or impossible to obtain during inference, such as parcellations, architectonic identity, transcriptomic information, and molecular profiles. By recognizing the mismatch between geometry and function, JOSA provides new insights into the future development of registration methods using joint analysis of brain structure and function.

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

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

Heritability of movie-evoked brain activity and connectivity

The neural bases of sensory processing are conserved across people, but no two individuals experience the same stimulus in exactly the same way. Recent work has established that the idiosyncratic nature of subjective experience is underpinned by individual variability in brain responses to sensory information. However, the fundamental origins of this individual variability have yet to be systematically investigated. Here, we establish a genetic basis for individual differences in sensory processing by quantifying (1) the heritability of high-dimensional brain responses to movies and (2) the extent to which this heritability is grounded in lower-level aspects of brain function. Specifically, we leverage 7T fMRI data collected from a twin sample to first show that movie-evoked brain activity and connectivity patterns are heritable across the cortex. Next, we use hyperalignment to decompose this heritability into genetic similarity in where vs. how sensory information is processed. Finally, we show that the heritability of brain activity patterns can be partially explained by the heritability of the neural timescale, a one-dimensional measure of local circuit functioning. These results demonstrate that brain responses to complex stimuli are heritable, and that this heritability is due, in part, to genetic control over stable aspects of brain function.

Dynamic layer-specific processing in the prefrontal cortex during working memory

The dorsolateral prefrontal cortex (dlPFC) is reliably engaged in working memory (WM) and comprises different cytoarchitectonic layers, yet their functional role in human WM is unclear. Here, participants completed a delayed-match-to-sample task while undergoing functional magnetic resonance imaging (fMRI) at ultra-high resolution. We examine layer-specific activity to manipulations in WM load and motor response. Superficial layers exhibit a preferential response to WM load during the delay and retrieval periods of a WM task, indicating a lamina-specific activation of the frontoparietal network. Multivariate patterns encoding WM load in the superficial layer dynamically change across the three periods of the task. Last, superficial and deep layers are non-differentially involved in the motor response, challenging earlier findings of a preferential deep layer activation. Taken together, our results provide new insights into the functional laminar circuitry of the dlPFC during WM and support a dynamic account of dlPFC coding.

β-Alanine supplementation improves fractional anisotropy scores in the hippocampus and amygdala in 60-80-year-old men and women

Recently, β-alanine (BA) supplementation was shown to improve cognitive function in older adults with decreased cognitive function. Mechanisms supporting these improvements have not been well defined. This study examined the effects of 10-weeks of BA supplementation on changes in circulating brain inflammatory markers, brain derived neurotrophic factor (BDNF), and brain morphology. Twenty participants were initially randomized into BA (2.4 g·d-1) or placebo (PL) groups. At each testing session, participants provided a resting blood sample and completed the Montreal cognitive assessment (MoCA) test and magnetic resonance imaging, which included diffusion tensor imaging to assess brain tissue integrity. Only participants that scored at or below normal for the MoCA assessment were analyzed (6 BA and 4 PL). The Mann-Whitney U test was used to examine Δ (POST-PRE) differences between the groups. No differences in Δ scores were noted in any blood marker (BDNF, CRP, TNF-α and GFAP). Changes in fractional anisotropy scores were significantly greater for BA than PL in the right hippocampus (p = 0.033) and the left amygdala (p = 0.05). No other differences were noted. The results provide a potential mechanism of how BA supplementation may improve cognitive function as reflected by improved tissue integrity within the hippocampus and amygdala.

Effects of connectivity hyperalignment (CHA) on estimated brain network properties: from coarse-scale to fine-scale

Recent gains in functional magnetic resonance imaging (fMRI) studies have been driven by increasingly sophisticated statistical and computational techniques and the ability to capture brain data at finer spatial and temporal resolution. These advances allow researchers to develop population-level models of the functional brain representations underlying behavior, performance, clinical status, and prognosis. However, even following conventional preprocessing pipelines, considerable inter-individual disparities in functional localization persist, posing a hurdle to performing compelling population-level inference. Persistent misalignment in functional topography after registration and spatial normalization will reduce power in developing predictive models and biomarkers, reduce the specificity of estimated brain responses and patterns, and provide misleading results on local neural representations and individual differences. This study aims to determine how connectivity hyperalignment (CHA)-an analytic approach for handling functional misalignment-can change estimated functional brain network topologies at various spatial scales from the coarsest set of parcels down to the vertex-level scale. The findings highlight the role of CHA in improving inter-subject similarities, while retaining individual-specific information and idiosyncrasies at finer spatial granularities. This highlights the potential for fine-grained connectivity analysis using this approach to reveal previously unexplored facets of brain structure and function.

A bimodal taxonomy of adult human brain sulcal morphology related to timing of fetal sulcation and trans-sulcal gene expression gradients

We developed a computational pipeline (now provided as a resource) for measuring morphological similarity between cortical surface sulci to construct a sulcal phenotype network (SPN) from each magnetic resonance imaging (MRI) scan in an adult cohort (n = 34,725; 45-82 years). Networks estimated from pairwise similarities of 40 sulci on 5 morphological metrics comprised two clusters of sulci, represented also by the bimodal distribution of sulci on a linear-to-complex dimension. Linear sulci were more heritable and typically located in unimodal cortex, and complex sulci were less heritable and typically located in heteromodal cortex. Aligning these results with an independent fetal brain MRI cohort (n = 228; 21-36 gestational weeks), we found that linear sulci formed earlier, and the earliest and latest-forming sulci had the least between-adult variation. Using high-resolution maps of cortical gene expression, we found that linear sulcation is mechanistically underpinned by trans-sulcal gene expression gradients enriched for developmental processes.

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

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

Dynamic adaptation to novelty in the brain is related to arousal and intelligence

How does the human brain respond to novelty? Here, we address this question using fMRI data wherein human participants watch the same movie scene four times. On the first viewing, this movie scene is novel, and on later viewings it is not. We find that brain activity is lower-dimensional in response to novelty. At a finer scale, we find that this reduction in the dimensionality of brain activity is the result of increased coupling in specific brain systems, most specifically within and between the control and dorsal attention systems. Additionally, we found that novelty induced an increase in between-subject synchronization of brain activity in the same brain systems. We also find evidence that adaptation to novelty, herein operationalized as the difference between baseline coupling and novelty-response coupling, is related to fluid intelligence. Finally, using separately collected out-of-sample data, we find that the above results may be linked to psychological arousal.

The heritability and structural correlates of resting-state fMRI complexity

The complexity of fMRI signals quantifies temporal dynamics of spontaneous neural activity, which has been increasingly recognized as providing important insights into cognitive functions and psychiatric disorders. However, its heritability and structural underpinnings are not well understood. Here, we utilize multi-scale sample entropy to extract resting-state fMRI complexity in a large healthy adult sample from the Human Connectome Project. We show that fMRI complexity at multiple time scales is heritable in broad brain regions. Heritability estimates are modest and regionally variable. We relate fMRI complexity to brain structure including surface area, cortical myelination, cortical thickness, subcortical volumes, and total brain volume. We find that surface area is negatively correlated with fine-scale complexity and positively correlated with coarse-scale complexity in most cortical regions, especially the association cortex. Most of these correlations are related to common genetic and environmental effects. We also find positive correlations between cortical myelination and fMRI complexity at fine scales and negative correlations at coarse scales in the prefrontal cortex, lateral temporal lobe, precuneus, lateral parietal cortex, and cingulate cortex, with these correlations mainly attributed to common environmental effects. We detect few significant associations between fMRI complexity and cortical thickness. Despite the non-significant association with total brain volume, fMRI complexity exhibits significant correlations with subcortical volumes in the hippocampus, cerebellum, putamen, and pallidum at certain scales. Collectively, our work establishes the genetic basis and structural correlates of resting-state fMRI complexity across multiple scales, supporting its potential application as an endophenotype for psychiatric disorders.

Macroscale brain states support the control of semantic cognition

A crucial aim in neuroscience is to understand how the human brain adapts to varying cognitive demands. This study investigates network reconfiguration during controlled semantic retrieval in differing contexts. We analyze brain responses to two semantic tasks of varying difficulty - global association and feature matching judgments - which are contrasted with non-semantic tasks on the cortical surface and within a whole-brain state space. Demanding semantic association tasks elicit activation in anterior prefrontal and temporal regions, while challenging semantic feature matching and non-semantic tasks predominantly activate posterior regions. Task difficulty also modulates activation along different dimensions of functional organization, suggesting different mechanisms of cognitive control. More demanding semantic association judgments engage cognitive control and default mode networks together, while feature matching and non-semantic tasks are skewed towards cognitive control networks. These findings highlight the brain's dynamic ability to tailor its networks to support diverse neurocognitive states, enriching our understanding of controlled cognition.

Longitudinally consistent registration and parcellation of cortical surfaces using semi-supervised learning

Temporally consistent and accurate registration and parcellation of longitudinal cortical surfaces is of great importance in studying longitudinal morphological and functional changes of human brains. However, most existing methods are developed for registration or parcellation of a single cortical surface. When applying to longitudinal studies, these methods independently register/parcellate each surface from longitudinal scans, thus often generating longitudinally inconsistent and inaccurate results, especially in small or ambiguous cortical regions. Essentially, longitudinal cortical surface registration and parcellation are highly correlated tasks with inherently shared constraints on both spatial and temporal feature representations, which are unfortunately ignored in existing methods. To this end, we unprecedentedly propose a novel semi-supervised learning framework to exploit these inherent relationships from limited labeled data and extensive unlabeled data for more robust and consistent registration and parcellation of longitudinal cortical surfaces. Our method utilizes the spherical topology characteristic of cortical surfaces. It employs a spherical network to function as an encoder, which extracts high-level cortical features. Subsequently, we build two specialized decoders dedicated to the tasks of registration and parcellation, respectively. To extract more meaningful spatial features, we design a novel parcellation map similarity loss to utilize the relationship between registration and parcellation tasks, i.e., the parcellation map warped by the deformation field in registration should match the atlas parcellation map, thereby providing extra supervision for the registration task and augmented data for parcellation task by warping the atlas parcellation map to unlabeled surfaces. To enable temporally more consistent feature representation, we additionally enforce longitudinal consistency among longitudinal surfaces after registering them together using their concatenated features. Experiments on two longitudinal datasets of infants and adults have shown that our method achieves significant improvements on both registration/parcellation accuracy and longitudinal consistency compared to existing methods, especially in small and challenging cortical regions.

A shifting role of thalamocortical connectivity in the emergence of cortical functional organization

The cortical patterning principle has been a long-standing question in neuroscience, yet how this translates to macroscale functional specialization in the human brain remains largely unknown. Here we examine age-dependent differences in resting-state thalamocortical connectivity to investigate its role in the emergence of large-scale functional networks during early life, using a primarily cross-sectional but also longitudinal approach. We show that thalamocortical connectivity during infancy reflects an early differentiation of sensorimotor networks and genetically influenced axonal projection. This pattern changes in childhood, when connectivity is established with the salience network, while decoupling externally and internally oriented functional systems. A developmental simulation using generative network models corroborated these findings, demonstrating that thalamic connectivity contributes to developing key features of the mature brain, such as functional segregation and the sensory-association axis, especially across 12-18 years of age. Our study suggests that the thalamus plays an important role in functional specialization during development, with potential implications for studying conditions with compromised internal and external processing.

A multidimensional investigation of sleep and biopsychosocial profiles with associated neural signatures

Sleep is essential for optimal functioning and health. Interconnected to multiple biological, psychological and socio-environmental factors (i.e., biopsychosocial factors), the multidimensional nature of sleep is rarely capitalized on in research. Here, we deployed a data-driven approach to identify sleep-biopsychosocial profiles that linked self-reported sleep patterns to inter-individual variability in health, cognition, and lifestyle factors in 770 healthy young adults. We uncovered five profiles, including two profiles reflecting general psychopathology associated with either reports of general poor sleep or an absence of sleep complaints (i.e., sleep resilience) respectively. The three other profiles were driven by the use of sleep aids and social satisfaction, sleep duration and cognitive performance, and sleep disturbance linked to cognition and mental health. Furthermore, identified sleep-biopsychosocial profiles displayed unique patterns of brain network organization. In particular, somatomotor network connectivity alterations were involved in the relationships between sleep and biopsychosocial factors. These profiles can potentially untangle the interplay between individuals' variability in sleep, health, cognition and lifestyle - equipping research and clinical settings to better support individual's well-being.

Profiling cortical morphometric similarity in perinatal brains: Insights from development, sex difference, and inter-individual variation

The topological organization of the macroscopic cortical networks important for the development of complex brain functions. However, how the cortical morphometric organization develops during the third trimester and whether it demonstrates sexual and individual differences at this particular stage remain unclear. Here, we constructed the morphometric similarity network (MSN) based on morphological and microstructural features derived from multimodal MRI of two independent cohorts (cross-sectional and longitudinal) scanned at 30-44 postmenstrual weeks (PMW). Sex difference and inter-individual variations of the MSN were also examined on these cohorts. The cross-sectional analysis revealed that both network integration and segregation changed in a nonlinear biphasic trajectory, which was supported by the results obtained from longitudinal analysis. The community structure showed remarkable consistency between bilateral hemispheres and maintained stability across PMWs. Connectivity within the primary cortex strengthened faster than that within high-order communities. Compared to females, male neonates showed a significant reduction in the participation coefficient within prefrontal and parietal cortices, while their overall network organization and community architecture remained comparable. Furthermore, by using the morphometric similarity as features, we achieved over 65 % accuracy in identifying an individual at term-equivalent age from images acquired after birth, and vice versa. These findings provide comprehensive insights into the development of morphometric similarity throughout the perinatal cortex, enhancing our understanding of the establishment of neuroanatomical organization during early life.

Spurious correlations in surface-based functional brain imaging

The study of functional MRI data is increasingly performed after mapping from volumetric voxels to surface vertices. Processing pipelines commonly used to achieve this mapping produce meshes with uneven vertex spacing, with closer neighbours in sulci compared to gyri. Consequently, correlations between the fMRI time series of neighbouring sulcal vertices are stronger than expected. However, the causes, extent, and impacts of this bias are not well understood or widely appreciated. We explain the origins of these biases, and using in-silico models of fMRI data, illustrate how they lead to spurious results. The bias leads to leakage of anatomical cortical folding information into fMRI time series. We show that many common analyses can be affected by this "gyral bias", including test-retest reliability, fingerprinting, functional parcellations, regional homogeneity, and brain-behaviour associations. Finally, we provide recommendations to avoid or remedy this spatial bias.

A multi-modal, asymmetric, weighted, and signed description of anatomical connectivity

The macroscale connectome is the network of physical, white-matter tracts between brain areas. The connections are generally weighted and their values interpreted as measures of communication efficacy. In most applications, weights are either assigned based on imaging features-e.g. diffusion parameters-or inferred using statistical models. In reality, the ground-truth weights are unknown, motivating the exploration of alternative edge weighting schemes. Here, we explore a multi-modal, regression-based model that endows reconstructed fiber tracts with directed and signed weights. We find that the model fits observed data well, outperforming a suite of null models. The estimated weights are subject-specific and highly reliable, even when fit using relatively few training samples, and the networks maintain a number of desirable features. In summary, we offer a simple framework for weighting connectome data, demonstrating both its ease of implementation while benchmarking its utility for typical connectome analyses, including graph theoretic modeling and brain-behavior associations.

A familiar face and person processing area in the human temporal pole

How does the brain process the faces of familiar people? Neuropsychological studies have argued for an area of the temporal pole (TP) linking faces with person identities, but magnetic susceptibility artifacts in this region have hampered its study with fMRI. Using data acquisition and analysis methods optimized to overcome this artifact, we identify a familiar face response in TP, reliably observed in individual brains. This area responds strongly to visual images of familiar faces over unfamiliar faces, objects, and scenes. However, TP did not just respond to images of faces, but also to a variety of high-level social cognitive tasks, including semantic, episodic, and theory of mind tasks. The response profile of TP contrasted with a nearby region of the perirhinal cortex that responded specifically to faces, but not to social cognition tasks. TP was functionally connected with a distributed network in the association cortex associated with social cognition, while PR was functionally connected with face-preferring areas of the ventral visual cortex. This work identifies a missing link in the human face processing system that specifically processes familiar faces, and is well placed to integrate visual information about faces with higher-order conceptual information about other people. The results suggest that separate streams for person and face processing reach anterior temporal areas positioned at the top of the cortical hierarchy.

Task and stimulus coding in the multiple-demand network

In the human brain, a multiple-demand (MD) network plays a key role in cognitive control, with core components in lateral frontal, dorsomedial frontal and lateral parietal cortex, and multivariate activity patterns that discriminate the contents of many cognitive activities. In prefrontal cortex of the behaving monkey, different cognitive operations are associated with very different patterns of neural activity, while details of a particular stimulus are encoded as small variations on these basic patterns (Sigala et al, 2008). Here, using the advanced fMRI methods of the Human Connectome Project and their 360-region cortical parcellation, we searched for a similar result in MD activation patterns. In each parcel, we compared multivertex patterns for every combination of three tasks (working memory, task-switching, and stop-signal) and two stimulus classes (faces and buildings). Though both task and stimulus category were discriminated in every cortical parcel, the strength of discrimination varied strongly across parcels. The different cognitive operations of the three tasks were strongly discriminated in MD regions. Stimulus categories, in contrast, were most strongly discriminated in a large region of primary and higher visual cortex, and intriguingly, in both parietal and frontal lobe regions adjacent to core MD regions. In the monkey, frontal neurons show a strong pattern of nonlinear mixed selectivity, with activity reflecting specific conjunctions of task events. In our data, however, there was limited evidence for mixed selectivity; throughout the brain, discriminations of task and stimulus combined largely linearly, with a small nonlinear component. In MD regions, human fMRI data recapitulate some but not all aspects of electrophysiological data from nonhuman primates.

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

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

The Multiscale Surface Vision Transformer

Surface meshes are a favoured domain for representing structural and functional information on the human cortex, but their complex topology and geometry pose significant challenges for deep learning analysis. While Transformers have excelled as domainagnostic architectures for sequence-to-sequence learning, the quadratic cost of the self-attention operation remains an obstacle for many dense prediction tasks. Inspired by some of the latest advances in hierarchical modelling with vision transformers, we introduce the Multiscale Surface Vision Transformer (MS-SiT) as a backbone architecture for surface deep learning. The self-attention mechanism is applied within local-mesh-windows to allow for high-resolution sampling of the underlying data, while a shifted-window strategy improves the sharing of information between windows. Neighbouring patches are successively merged, allowing the MS-SiT to learn hierarchical representations suitable for any prediction task. Results demonstrate that the MS-SiT outperforms existing surface deep learning methods for neonatal phenotyping prediction tasks using the Developing Human Connectome Project (dHCP) dataset. Furthermore, building the MS-SiT backbone into a U-shaped architecture for surface segmentation demonstrates competitive results on cortical parcellation using the UK Biobank (UKB) and manually-annotated MindBoggle datasets. Code and trained models are publicly available at https://github.com/metrics-lab/surface-vision-transformers.

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.

Mapping vascular network architecture in primate brain using ferumoxytol-weighted laminar MRI

Mapping the vascular organization of the brain is of great importance across various domains of basic neuroimaging research, diagnostic radiology, and neurology. However, the intricate task of precisely mapping vasculature across brain regions and cortical layers presents formidable challenges, resulting in a limited understanding of neurometabolic factors influencing the brain's microvasculature. Addressing this gap, our study investigates whole-brain vascular volume using ferumoxytol-weighted laminar-resolution multi-echo gradient-echo imaging in macaque monkeys. We validate the results with published data for vascular densities and compare them with cytoarchitecture, neuron and synaptic densities. The ferumoxytol-induced change in transverse relaxation rate (ΔR2*), an indirect proxy measure of cerebral blood volume (CBV), was mapped onto twelve equivolumetric laminar cortical surfaces. Our findings reveal that CBV varies 3-fold across the brain, with the highest vascular volume observed in the inferior colliculus and lowest in the corpus callosum. In the cerebral cortex, CBV is notably high in early primary sensory areas and low in association areas responsible for higher cognitive functions. Classification of CBV into distinct groups unveils extensive replication of translaminar vascular network motifs, suggesting distinct computational energy supply requirements in areas with varying cytoarchitecture types. Regionally, baseline R2* and CBV exhibit positive correlations with neuron density and negative correlations with receptor densities. Adjusting image resolution based on the critical sampling frequency of penetrating cortical vessels, allows us to delineate approximately 30% of the arterial-venous vessels. Collectively, these results mark significant methodological and conceptual advancements, contributing to the refinement of cerebrovascular MRI. Furthermore, our study establishes a linkage between neurometabolic factors and the vascular network architecture in the primate brain.

Multiscale Modes of Functional Brain Connectivity

Information processing in the brain spans from localised sensorimotor processes to higher-level cognition that integrates across multiple regions. Interactions between and within these subsystems enable multiscale information processing. Despite this multiscale characteristic, functional brain connectivity is often either estimated based on 10-30 distributed modes or parcellations with 100-1000 localised parcels, both missing across-scale functional interactions. We present Multiscale Probabilistic Functional Modes (mPFMs), a new mapping which comprises modes over various scales of granularity, thus enabling direct estimation of functional connectivity within- and across-scales. Crucially, mPFMs emerged from data-driven multilevel Bayesian modelling of large functional MRI (fMRI) populations. We demonstrate that mPFMs capture both distributed brain modes and their co-existing subcomponents. In addition to validating mPFMs using simulations and real data, we show that mPFMs can predict ~900 personalised traits from UK Biobank more accurately than current standard techniques. Therefore, mPFMs can offer a paradigm shift in functional connectivity modelling and yield enhanced fMRI biomarkers for traits and diseases.

BAI-Net: Individualized Anatomical Cerebral Cartography Using Graph Neural Network

Brain atlas is an important tool in the diagnosis and treatment of neurological disorders. However, due to large variations in the organizational principles of individual brains, many challenges remain in clinical applications. Brain atlas individualization network (BAI-Net) is an algorithm that subdivides individual cerebral cortex into segregated areas using brain morphology and connectomes. The presented method integrates group priors derived from a population atlas, adjusts areal probabilities using the context of connectivity fingerprints derived from the fiber-tract embedding of tractography, and provides reliable and explainable individualized brain areas across multiple sessions and scanners. We demonstrate that BAI-Net outperforms the conventional iterative clustering approach by capturing significantly heritable topographic variations in individualized cartographies. The topographic variability of BAI-Net cartographies has shown strong associations with individual variability in brain morphology, connectivity as well as higher relationship on individual cognitive behaviors and genetics. This study provides an explainable framework for individualized brain cartography that may be useful in the precise localization of neuromodulation and treatments on individual brains.

The Brain's Topographical Organization Shapes Dynamic Interaction Patterns That Support Flexible Behavior Based on Rules and Long-Term Knowledge

Adaptive behavior relies both on specific rules that vary across situations and stable long-term knowledge gained from experience. The frontoparietal control network (FPCN) is implicated in the brain's ability to balance these different influences on action. Here, we investigate how the topographical organization of the cortex supports behavioral flexibility within the FPCN. Functional properties of this network might reflect its juxtaposition between the dorsal attention network (DAN) and the default mode network (DMN), two large-scale systems implicated in top-down attention and memory-guided cognition, respectively. Our study tests whether subnetworks of FPCN are topographically proximal to the DAN and the DMN, respectively, and how these topographical differences relate to functional differences: the proximity of each subnetwork is anticipated to play a pivotal role in generating distinct cognitive modes relevant to working memory and long-term memory. We show that FPCN subsystems share multiple anatomical and functional similarities with their neighboring systems (DAN and DMN) and that this topographical architecture supports distinct interaction patterns that give rise to different patterns of functional behavior. The FPCN acts as a unified system when long-term knowledge supports behavior but becomes segregated into discrete subsystems with different patterns of interaction when long-term memory is less relevant. In this way, our study suggests that the topographical organization of the FPCN and the connections it forms with distant regions of cortex are important influences on how this system supports flexible behavior.

SUGAR: Spherical ultrafast graph attention framework for cortical surface registration

Cortical surface registration plays a crucial role in aligning cortical functional and anatomical features across individuals. However, conventional registration algorithms are computationally inefficient. Recently, learning-based registration algorithms have emerged as a promising solution, significantly improving processing efficiency. Nonetheless, there remains a gap in the development of a learning-based method that exceeds the state-of-the-art conventional methods simultaneously in computational efficiency, registration accuracy, and distortion control, despite the theoretically greater representational capabilities of deep learning approaches. To address the challenge, we present SUGAR, a unified unsupervised deep-learning framework for both rigid and non-rigid registration. SUGAR incorporates a U-Net-based spherical graph attention network and leverages the Euler angle representation for deformation. In addition to the similarity loss, we introduce fold and multiple distortion losses to preserve topology and minimize various types of distortions. Furthermore, we propose a data augmentation strategy specifically tailored for spherical surface registration to enhance the registration performance. Through extensive evaluation involving over 10,000 scans from 7 diverse datasets, we showed that our framework exhibits comparable or superior registration performance in accuracy, distortion, and test-retest reliability compared to conventional and learning-based methods. Additionally, SUGAR achieves remarkable sub-second processing times, offering a notable speed-up of approximately 12,000 times in registering 9,000 subjects from the UK Biobank dataset in just 32 min. This combination of high registration performance and accelerated processing time may greatly benefit large-scale neuroimaging studies.

Functional Localization of the Human Auditory and Visual Thalamus Using a Thalamic Localizer Functional Magnetic Resonance Imaging Task

Functional magnetic resonance imaging (fMRI) of the auditory and visual sensory systems of the human brain is an active area of investigation in the study of human health and disease. The medial geniculate nucleus (MGN) and lateral geniculate nucleus (LGN) are key thalamic nuclei involved in the processing and relay of auditory and visual information, respectively, and are the subject of blood-oxygen-level-dependent (BOLD) fMRI studies of neural activation and functional connectivity in human participants. However, localization of BOLD fMRI signal originating from neural activity in MGN and LGN remains a technical challenge, due in part to the poor definition of boundaries of these thalamic nuclei in standard T1-weighted and T2-weighted magnetic resonance imaging sequences. Here, we report the development and evaluation of an auditory and visual sensory thalamic localizer (TL) fMRI task that produces participant-specific functionally-defined regions of interest (fROIs) of both MGN and LGN, using 3 Tesla multiband fMRI and a clustered-sparse temporal acquisition sequence, in less than 16 minutes of scan time. We demonstrate the use of MGN and LGN fROIs obtained from the TL fMRI task in standard resting-state functional connectivity (RSFC) fMRI analyses in the same participants. In RSFC analyses, we validated the specificity of MGN and LGN fROIs for signals obtained from primary auditory and visual cortex, respectively, and benchmark their performance against alternative atlas- and segmentation-based localization methods. The TL fMRI task and analysis code (written in Presentation and MATLAB, respectively) have been made freely available to the wider research community.

Species -shared and -unique gyral peaks on human and macaque brains

Cortical folding is an important feature of primate brains that plays a crucial role in various cognitive and behavioral processes. Extensive research has revealed both similarities and differences in folding morphology and brain function among primates including macaque and human. The folding morphology is the basis of brain function, making cross-species studies on folding morphology important for understanding brain function and species evolution. However, prior studies on cross-species folding morphology mainly focused on partial regions of the cortex instead of the entire brain. Previously, our research defined a whole-brain landmark based on folding morphology: the gyral peak. It was found to exist stably across individuals and ages in both human and macaque brains. Shared and unique gyral peaks in human and macaque are identified in this study, and their similarities and differences in spatial distribution, anatomical morphology, and functional connectivity were also dicussed.

CiftiStorm pipeline: facilitating reproducible EEG/MEG source connectomics

We present CiftiStorm, an electrophysiological source imaging (ESI) pipeline incorporating recently developed methods to improve forward and inverse solutions. The CiftiStorm pipeline produces Human Connectome Project (HCP) and megconnectome-compliant outputs from dataset inputs with varying degrees of spatial resolution. The input data can range from low-sensor-density electroencephalogram (EEG) or magnetoencephalogram (MEG) recordings without structural magnetic resonance imaging (sMRI) to high-density EEG/MEG recordings with an HCP multimodal sMRI compliant protocol. CiftiStorm introduces a numerical quality control of the lead field and geometrical corrections to the head and source models for forward modeling. For the inverse modeling, we present a Bayesian estimation of the cross-spectrum of sources based on multiple priors. We facilitate ESI in the T1w/FSAverage32k high-resolution space obtained from individual sMRI. We validate this feature by comparing CiftiStorm outputs for EEG and MRI data from the Cuban Human Brain Mapping Project (CHBMP) acquired with technologies a decade before the HCP MEG and MRI standardized dataset.

A multidimensional investigation of sleep and biopsychosocial profiles with associated neural signatures

Sleep is essential for optimal functioning and health. Interconnected to multiple biological, psychological and socio-environmental factors (i.e., biopsychosocial factors), the multidimensional nature of sleep is rarely capitalized on in research. Here, we deployed a data-driven approach to identify sleep-biopsychosocial profiles that linked self-reported sleep patterns to inter-individual variability in health, cognition, and lifestyle factors in 770 healthy young adults. We uncovered five profiles, including two profiles reflecting general psychopathology associated with either reports of general poor sleep or an absence of sleep complaints (i.e., sleep resilience) respectively. The three other profiles were driven by sedative-hypnotics-use and social satisfaction, sleep duration and cognitive performance, and sleep disturbance linked to cognition and mental health. Furthermore, identified sleep-biopsychosocial profiles displayed unique patterns of brain network organization. In particular, somatomotor network connectivity alterations were involved in the relationships between sleep and biopsychosocial factors. These profiles can potentially untangle the interplay between individuals' variability in sleep, health, cognition and lifestyle - equipping research and clinical settings to better support individual's well-being.

Preterm Birth Alters the Regional Development and Structural Covariance of Cerebellum at Term-Equivalent Age

Preterm birth is associated with increased risk for a spectrum of neurodevelopmental disabilities. The cerebellum is implicated in a wide range of cognitive functions extending beyond sensorimotor control and plays an increasingly recognized role in brain development. Morphometric studies based on volume analyses have revealed impaired cerebellar development in preterm infants. However, the structural covariance between the cerebellum and cerebral cortex has not been studied during the neonatal period, and the extent to which structural covariance is affected by preterm birth remains unknown. In this study, using the structural MR images of 52 preterm infants scanned at term-equivalent age and 312 full-term controls from the Developing Human Connectome Project, we compared volumetric growth, local cerebellum shape development and cerebello-cerebral structural covariance between the two groups. We found that although there was no significant difference in the overall volume measurements between preterm and full-term infants, the shape measurements were different. Compared with the control infants, preterm infants had significantly larger thickness in the vermis and lower thickness in the lateral portions of the bilateral cerebral hemispheres. The structural covariance between the cerebellum and frontal and parietal lobes was significantly greater in preterm infants than in full-term controls. The findings in this study suggested that cerebellar development and cerebello-cerebral structural covariance may be affected by premature birth.

Whole-Brain Dynamics Disruptions in the Progression of Alzheimer's Disease: Understanding the Influence of Amyloid-Beta and Tau

INTRODUCTION

Alzheimer's disease (AD) affects brain structure and function along its evolution, but brain network dynamic changes remain largely unknown.

METHODS

To understand how AD shapes brain activity, we investigated the spatiotemporal dynamics and resting state functional networks using the intrinsic ignition framework, which characterizes how an area transmits neuronal activity to others, resulting in different degrees of integration. Healthy participants, MCI, and AD patients were scanned using resting state fMRI. Mixed effects models were used to assess the impact of ABeta and tau, at the regional and whole-brain levels.

RESULTS

Dynamic complexity is progressively reduced, with Healthy participants showing higher metastability (i.e., a more complex dynamical regime over time) than observed in the other stages, while AD subjects showed the lowest.

DISCUSSION

Our study provides further insight into how AD modulates brain network dynamics along its evolution, progressively disrupting the whole-brain and resting state network dynamics.

Bayesian Optimization of Neurostimulation (BOONStim)

BACKGROUND

Transcranial magnetic stimulation (TMS) treatment response is influenced by individual variability in brain structure and function. Sophisticated, user-friendly approaches, incorporating both established functional magnetic resonance imaging (fMRI) and TMS simulation tools, to identify TMS targets are needed.

OBJECTIVE

The current study presents the development and validation of the Bayesian Optimization of Neuro-Stimulation (BOONStim) pipeline.

METHODS

BOONStim uses Bayesian optimization for individualized TMS targeting, automating interoperability between surface-based fMRI analytic tools and TMS electric field modeling. Bayesian optimization performance was evaluated in a sample dataset (N=10) using standard circular and functional connectivity-defined targets, and compared to grid optimization.

RESULTS

Bayesian optimization converged to similar levels of total electric field stimulation across targets in under 30 iterations, converging within a 5% error of the maxima detected by grid optimization, and requiring less time.

CONCLUSIONS

BOONStim is a scalable and configurable user-friendly pipeline for individualized TMS targeting with quick turnaround.

Co-representation of Functional Brain Networks Is Shaped by Cortical Myeloarchitecture and Reveals Individual Behavioral Ability

Large-scale functional networks are spatially distributed in the human brain. Despite recent progress in differentiating their functional roles, how the brain navigates the spatial coordination among them and the biological relevance of this coordination is still not fully understood. Capitalizing on canonical individualized networks derived from functional MRI data, we proposed a new concept, that is, co-representation of functional brain networks, to delineate the spatial coordination among them. To further quantify the co-representation pattern, we defined two indexes, that is, the co-representation specificity (CoRS) and intensity (CoRI), for separately measuring the extent of specific and average expression of functional networks at each brain location by using the data from both sexes. We found that the identified pattern of co-representation was anchored by cortical regions with three types of cytoarchitectural classes along a sensory-fugal axis, including, at the first end, primary (idiotypic) regions showing high CoRS, at the second end, heteromodal regions showing low CoRS and high CoRI, at the third end, paralimbic regions showing low CoRI. Importantly, we demonstrated the critical role of myeloarchitecture in sculpting the spatial distribution of co-representation by assessing the association with the myelin-related neuroanatomical and transcriptomic profiles. Furthermore, the significance of manifesting the co-representation was revealed in its prediction of individual behavioral ability. Our findings indicated that the spatial coordination among functional networks was built upon an anatomically configured blueprint to facilitate neural information processing, while advancing our understanding of the topographical organization of the brain by emphasizing the assembly of functional networks.

Intra-V1 functional networks and classification of observed stimuli

Introduction

Previous studies suggest that co-fluctuations in neural activity within V1 (measured with fMRI) carry information about observed stimuli, potentially reflecting various cognitive mechanisms. This study explores the neural sources shaping this information by using different fMRI preprocessing methods. The common response to stimuli shared by all individuals can be emphasized by using inter-subject correlations or de-emphasized by deconvolving the fMRI with hemodynamic response functions (HRFs) before calculating the correlations. The latter approach shifts the balance towards participant-idiosyncratic activity.

Methods

Here, we used multivariate pattern analysis of intra-V1 correlation matrices to predict the Level or Shape of observed Navon letters employing the types of correlations described above. We assessed accuracy in inter-subject prediction of specific conjunctions of properties, and attempted intra-subject cross-classification of stimulus properties (i.e., prediction of one feature despite changes in the other). Weight maps from successful classifiers were projected onto the visual field. A control experiment investigated eye-movement patterns during stimuli presentation.

Results

All inter-subject classifiers accurately predicted the Level and Shape of specific observed stimuli. However, successful intra-subject cross-classification was achieved only for stimulus Level, but not Shape, regardless of preprocessing scheme. Weight maps for successful Level classification differed between inter-subject correlations and deconvolved correlations. The latter revealed asymmetries in visual field link strength that corresponded to known perceptual asymmetries. Post-hoc measurement of eyeball fMRI signals did not find differences in gaze between stimulus conditions, and a control experiment (with derived simulations) also suggested that eye movements do not explain the stimulus-related changes in V1 topology.

Discussion

Our findings indicate that both inter-subject common responses and participant-specific activity contribute to the information in intra-V1 co-fluctuations, albeit through distinct sub-networks. Deconvolution, that enhances subject-specific activity, highlighted interhemispheric links for Global stimuli. Further exploration of intra-V1 networks promises insights into the neural basis of attention and perceptual organization.

Hierarchical functional differences between gyri and sulci at different scales

Gyri and sulci are 2 fundamental cortical folding patterns of the human brain. Recent studies have suggested that gyri and sulci may play different functional roles given their structural and functional heterogeneity. However, our understanding of the functional differences between gyri and sulci remains limited due to several factors. Firstly, previous studies have typically focused on either the spatial or temporal domain, neglecting the inherently spatiotemporal nature of brain functions. Secondly, analyses have often been restricted to either local or global scales, leaving the question of hierarchical functional differences unresolved. Lastly, there has been a lack of appropriate analytical tools for interpreting the hierarchical spatiotemporal features that could provide insights into these differences. To overcome these limitations, in this paper, we proposed a novel hierarchical interpretable autoencoder (HIAE) to explore the hierarchical functional difference between gyri and sulci. Central to our approach is its capability to extract hierarchical features via a deep convolutional autoencoder and then to map these features into an embedding vector using a carefully designed feature interpreter. This process transforms the features into interpretable spatiotemporal patterns, which are pivotal in investigating the functional disparities between gyri and sulci. We evaluate the proposed framework on Human Connectome Project task functional magnetic resonance imaging dataset. The experiments demonstrate that the HIAE model can effectively extract and interpret hierarchical spatiotemporal features that are neuroscientifically meaningful. The analyses based on the interpreted features suggest that gyri are more globally activated, whereas sulci are more locally activated, demonstrating a distinct transition in activation patterns as the scale shifts from local to global. Overall, our study provides novel insights into the brain's anatomy-function relationship.

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.

Local gradient analysis of human brain function using the Vogt-Bailey Index

In this work, we take a closer look at the Vogt-Bailey (VB) index, proposed in Bajada et al. (NeuroImage 221:117140, 2020) as a tool for studying local functional homogeneity in the human cortex. We interpret the VB index in terms of the minimum ratio cut, a scaled cut-set weight that indicates whether a network can easily be disconnected into two parts having a comparable number of nodes. In our case, the nodes of the network consist of a brain vertex/voxel and its neighbours, and a given edge is weighted according to the affinity of the nodes it connects (as reflected by the modified Pearson correlation between their fMRI time series). Consequently, the minimum ratio cut quantifies the degree of small-scale similarity in brain activity: the greater the similarity, the 'heavier' the edges and the more difficult it is to disconnect the network, hence the higher the value of the minimum ratio cut. We compare the performance of the VB index with that of the Regional Homogeneity (ReHo) algorithm, commonly used to assess whether voxels in close proximity have synchronised fMRI signals, and find that the VB index is uniquely placed to detect sharp changes in the (local) functional organization of the human cortex.

Controlling the human connectome with spatially diffuse input signals

The human brain is never at "rest"; its activity is constantly fluctuating over time, transitioning from one brain state-a whole-brain pattern of activity-to another. Network control theory offers a framework for understanding the effort - energy - associated with these transitions. One branch of control theory that is especially useful in this context is "optimal control", in which input signals are used to selectively drive the brain into a target state. Typically, these inputs are introduced independently to the nodes of the network (each input signal is associated with exactly one node). Though convenient, this input strategy ignores the continuity of cerebral cortex - geometrically, each region is connected to its spatial neighbors, allowing control signals, both exogenous and endogenous, to spread from their foci to nearby regions. Additionally, the spatial specificity of brain stimulation techniques is limited, such that the effects of a perturbation are measurable in tissue surrounding the stimulation site. Here, we adapt the network control model so that input signals have a spatial extent that decays exponentially from the input site. We show that this more realistic strategy takes advantage of spatial dependencies in structural connectivity and activity to reduce the energy (effort) associated with brain state transitions. We further leverage these dependencies to explore near-optimal control strategies such that, on a per-transition basis, the number of input signals required for a given control task is reduced, in some cases by two orders of magnitude. This approximation yields network-wide maps of input site density, which we compare to an existing database of functional, metabolic, genetic, and neurochemical maps, finding a close correspondence. Ultimately, not only do we propose a more efficient framework that is also more adherent to well-established brain organizational principles, but we also posit neurobiologically grounded bases for optimal control.

On the stability of canonical correlation analysis and partial least squares with application to brain-behavior associations

Associations between datasets can be discovered through multivariate methods like Canonical Correlation Analysis (CCA) or Partial Least Squares (PLS). A requisite property for interpretability and generalizability of CCA/PLS associations is stability of their feature patterns. However, stability of CCA/PLS in high-dimensional datasets is questionable, as found in empirical characterizations. To study these issues systematically, we developed a generative modeling framework to simulate synthetic datasets. We found that when sample size is relatively small, but comparable to typical studies, CCA/PLS associations are highly unstable and inaccurate; both in their magnitude and importantly in the feature pattern underlying the association. We confirmed these trends across two neuroimaging modalities and in independent datasets with n ≈ 1000 and n = 20,000, and found that only the latter comprised sufficient observations for stable mappings between imaging-derived and behavioral features. We further developed a power calculator to provide sample sizes required for stability and reliability of multivariate analyses. Collectively, we characterize how to limit detrimental effects of overfitting on CCA/PLS stability, and provide recommendations for future studies.

Mapping the functional and structural connectivity of the scene network

The recognition and perception of places has been linked to a network of scene-selective regions in the human brain. While previous studies have focussed on functional connectivity between scene-selective regions themselves, less is known about their connectivity with other cortical and subcortical regions in the brain. Here, we determine the functional and structural connectivity profile of the scene network. We used fMRI to examine functional connectivity between scene regions and across the whole brain during rest and movie-watching. Connectivity within the scene network revealed a bias between posterior and anterior scene regions implicated in perceptual and mnemonic aspects of scene perception respectively. Differences between posterior and anterior scene regions were also evident in the connectivity with cortical and subcortical regions across the brain. For example, the Occipital Place Area (OPA) and posterior Parahippocampal Place Area (PPA) showed greater connectivity with visual and dorsal attention networks, while anterior PPA and Retrosplenial Complex showed preferential connectivity with default mode and frontoparietal control networks and the hippocampus. We further measured the structural connectivity of the scene network using diffusion tractography. This indicated both similarities and differences with the functional connectivity, highlighting biases between posterior and anterior regions, but also between ventral and dorsal scene regions. Finally, we quantified the structural connectivity between the scene network and major white matter tracts throughout the brain. These findings provide a map of the functional and structural connectivity of scene-selective regions to each other and the rest of the brain.

Molecular signatures of cortical expansion in the human fetal brain

The third trimester of human gestation is characterised by rapid increases in brain volume and cortical surface area. A growing catalogue of cells in the prenatal brain has revealed remarkable molecular diversity across cortical areas.1,2 Despite this, little is known about how this translates into the patterns of differential cortical expansion observed in humans during the latter stages of gestation. Here we present a new resource, μBrain, to facilitate knowledge translation between molecular and anatomical descriptions of the prenatal developing brain. Built using generative artificial intelligence, μBrain is a three-dimensional cellular-resolution digital atlas combining publicly-available serial sections of the postmortem human brain at 21 weeks gestation3 with bulk tissue microarray data, sampled across 29 cortical regions and 5 transient tissue zones.4 Using μBrain, we evaluate the molecular signatures of preferentially-expanded cortical regions during human gestation, quantified in utero using magnetic resonance imaging (MRI). We find that differences in the rates of expansion across cortical areas during gestation respect anatomical and evolutionary boundaries between cortical types5 and are founded upon extended periods of upper-layer cortical neuron migration that continue beyond mid-gestation. We identify a set of genes that are upregulated from mid-gestation and highly expressed in rapidly expanding neocortex, which are implicated in genetic disorders with cognitive sequelae. Our findings demonstrate a spatial coupling between areal differences in the timing of neurogenesis and rates of expansion across the neocortical sheet during the prenatal epoch. The μBrain atlas is available from: https://garedaba.github.io/micro-brain/ and provides a new tool to comprehensively map early brain development across domains, model systems and resolution scales.

Dissociation of Reliability, Heritability, and Predictivity in Coarse- and Fine-Scale Functional Connectomes during Development

The functional connectome supports information transmission through the brain at various spatial scales, from exchange between broad cortical regions to finer-scale, vertex-wise connections that underlie specific information processing mechanisms. In adults, while both the coarse- and fine-scale functional connectomes predict cognition, the fine scale can predict up to twice the variance as the coarse-scale functional connectome. Yet, past brain-wide association studies, particularly using large developmental samples, focus on the coarse connectome to understand the neural underpinnings of individual differences in cognition. Using a large cohort of children (age 9-10 years; n = 1,115 individuals; both sexes; 50% female, including 170 monozygotic and 219 dizygotic twin pairs and 337 unrelated individuals), we examine the reliability, heritability, and behavioral relevance of resting-state functional connectivity computed at different spatial scales. We use connectivity hyperalignment to improve access to reliable fine-scale (vertex-wise) connectivity information and compare the fine-scale connectome with the traditional parcel-wise (coarse scale) functional connectomes. Though individual differences in the fine-scale connectome are more reliable than those in the coarse-scale, they are less heritable. Further, the alignment and scale of connectomes influence their ability to predict behavior, whereby some cognitive traits are equally well predicted by both connectome scales, but other, less heritable cognitive traits are better predicted by the fine-scale connectome. Together, our findings suggest there are dissociable individual differences in information processing represented at different scales of the functional connectome which, in turn, have distinct implications for heritability and cognition.

Species -Shared and -Unique Gyral Peaks on Human and Macaque Brains

Cortical folding is an important feature of primate brains that plays a crucial role in various cognitive and behavioral processes. Extensive research has revealed both similarities and differences in folding morphology and brain function among primates including macaque and human. The folding morphology is the basis of brain function, making cross-species studies on folding morphology important for understanding brain function and species evolution. However, prior studies on cross-species folding morphology mainly focused on partial regions of the cortex instead of the entire brain. Previously, we defined a whole-brain landmark based on folding morphology: the gyral peak. It was found to exist stably across individuals and ages in both human and macaque brains. In this study, we identified shared and unique gyral peaks in human and macaque, and investigated the similarities and differences in the spatial distribution, anatomical morphology, and functional connectivity of them.

Basis of executive functions in fine-grained architecture of cortical and subcortical human brain networks

Theoretical models suggest that executive functions rely on both domain-general and domain-specific processes. Supporting this view, prior brain imaging studies have revealed that executive activations converge and diverge within broadly characterized brain networks. However, the lack of precise anatomical mappings has impeded our understanding of the interplay between domain-general and domain-specific processes. To address this challenge, we used the high-resolution multimodal magnetic resonance imaging approach of the Human Connectome Project to scan participants performing 3 canonical executive tasks: n-back, rule switching, and stop signal. The results reveal that, at the individual level, different executive activations converge within 9 domain-general territories distributed in frontal, parietal, and temporal cortices. Each task exhibits a unique topography characterized by finely detailed activation gradients within domain-general territory shifted toward adjacent resting-state networks; n-back activations shift toward the default mode, rule switching toward dorsal attention, and stop signal toward cingulo-opercular networks. Importantly, the strongest activations arise at multimodal neurobiological definitions of network borders. Matching results are seen in circumscribed regions of the caudate nucleus, thalamus, and cerebellum. The shifting peaks of local gradients at the intersection of task-specific networks provide a novel mechanistic insight into how partially-specialized networks interact with neighboring domain-general territories to generate distinct executive functions.

Modular subgraphs in large-scale connectomes underpin spontaneous co-fluctuation events in mouse and human brains

Previous studies have adopted an edge-centric framework to study fine-scale network dynamics in human fMRI. To date, however, no studies have applied this framework to data collected from model organisms. Here, we analyze structural and functional imaging data from lightly anesthetized mice through an edge-centric lens. We find evidence of "bursty" dynamics and events - brief periods of high-amplitude network connectivity. Further, we show that on a per-frame basis events best explain static FC and can be divided into a series of hierarchically-related clusters. The co-fluctuation patterns associated with each cluster centroid link distinct anatomical areas and largely adhere to the boundaries of algorithmically detected functional brain systems. We then investigate the anatomical connectivity undergirding high-amplitude co-fluctuation patterns. We find that events induce modular bipartitions of the anatomical network of inter-areal axonal projections. Finally, we replicate these same findings in a human imaging dataset. In summary, this report recapitulates in a model organism many of the same phenomena observed in previously edge-centric analyses of human imaging data. However, unlike human subjects, the murine nervous system is amenable to invasive experimental perturbations. Thus, this study sets the stage for future investigation into the causal origins of fine-scale brain dynamics and high-amplitude co-fluctuations. Moreover, the cross-species consistency of the reported findings enhances the likelihood of future translation.

Automating the Human Connectome Project's Temporal ICA Pipeline

Functional magnetic resonance imaging (fMRI) data are dominated by noise and artifacts, with only a small fraction of the variance relating to neural activity. Temporal independent component analysis (tICA) is a recently developed method that enables selective denoising of fMRI artifacts related to physiology such as respiration. However, an automated and easy to use pipeline for tICA has not previously been available; instead, two manual steps have been necessary: 1) setting the group spatial ICA dimensionality after MELODIC's Incremental Group-PCA (MIGP) and 2) labeling tICA components as artifacts versus signals. Moreover, guidance has been lacking as to how many subjects and timepoints are needed to adequately re-estimate the temporal ICA decomposition and what alternatives are available for smaller groups or even individual subjects. Here, we introduce a nine-step fully automated tICA pipeline which removes global artifacts from fMRI dense timeseries after sICA+FIX cleaning and MSMAll alignment driven by functionally relevant areal features. Additionally, we have developed an automated "reclean" Pipeline for improved spatial ICA (sICA) artifact removal. Two major automated components of the pipeline are 1) an automatic group spatial ICA (sICA) dimensionality selection for MIGP data enabled by fitting multiple Wishart distributions; 2) a hierarchical classifier to distinguish group tICA signal components from artifactual components, equipped with a combination of handcrafted features from domain expert knowledge and latent features obtained via self-supervised learning on spatial maps. We demonstrate that the dimensionality estimated for the MIGP data from HCP Young Adult 3T and 7T datasets is comparable to previous manual tICA estimates, and that the group sICA decomposition is highly reproducible. We also show that the tICA classifier achieved over 0.98 Precision-Recall Area Under Curve (PR-AUC) and that the correctly classified components account for over 95% of the tICA-represented variance on multiple held-out evaluation datasets including the HCP-Young Adult, HCP-Aging and HCP-Development datasets under various settings. Our automated tICA pipeline is now available as part of the HCP pipelines, providing a powerful and user-friendly tool for the neuroimaging community.

Static and dynamic fMRI-derived functional connectomes represent largely similar information

Functional connectivity (FC) of blood oxygen level-dependent (BOLD) fMRI time series can be estimated using methods that differ in sensitivity to the temporal order of time points (static vs. dynamic) and the number of regions considered in estimating a single edge (bivariate vs. multivariate). Previous research suggests that dynamic FC explains variability in FC fluctuations and behavior beyond static FC. Our aim was to systematically compare methods on both dimensions. We compared five FC methods: Pearson's/full correlation (static, bivariate), lagged correlation (dynamic, bivariate), partial correlation (static, multivariate), and multivariate AR model with and without self-connections (dynamic, multivariate). We compared these methods by (i) assessing similarities between FC matrices, (ii) by comparing node centrality measures, and (iii) by comparing the patterns of brain-behavior associations. Although FC estimates did not differ as a function of sensitivity to temporal order, we observed differences between the multivariate and bivariate FC methods. The dynamic FC estimates were highly correlated with the static FC estimates, especially when comparing group-level FC matrices. Similarly, there were high correlations between the patterns of brain-behavior associations obtained using the dynamic and static FC methods. We conclude that the dynamic FC estimates represent information largely similar to that of the static FC.