Sleep deprivation leads to significant impairments in cognitive performance and changes to the interactions between large scale cortical networks, yet the hierarchical organization of cortical activity across states is still being explored. We used functional magnetic resonance imaging to assess activations and connectivity during cognitive tasks in 20 healthy young adults, during three states: (i) following a normal night of sleep, (ii) following 24hr of total sleep deprivation, and (iii) after a morning recovery nap. Situating cortical activity during cognitive tasks along hierarchical organizing gradients based upon similarity of functional connectivity patterns, we found that regional variations in task-activations were captured by an axis differentiating areas involved in executive control from default mode regions and paralimbic cortex. After global signal regression, the range of functional differentiation along this axis at baseline was significantly related to decline in working memory performance (2-back task) following sleep deprivation, as well as the extent of recovery in performance following a nap. The relative positions of cortical regions within gradients did not significantly change across states, except for a lesser differentiation of the visual system and increased coupling of the posterior cingulate cortex with executive control areas after sleep deprivation. This was despite a widespread increase in the magnitude of functional connectivity across the cortex following sleep deprivation. Cortical gradients of functional differentiation thus appear relatively insensitive to state-dependent changes following sleep deprivation and recovery, suggesting that there are no large-scale changes in cortical functional organization across vigilance states. Certain features of particular gradient axes may be informative for the extent of decline in performance on more complex tasks following sleep deprivation, and could be beneficial over traditional voxel- or parcel-based approaches in identifying realtionships between state-dependent brain activity and behavior.

Analysis and interpretation of neuroimaging datasets has become a multidisciplinary endeavor, relying not only on statistical methods, but increasingly on associations with respect to other brain-derived features such as gene expression, histological data, and functional as well as cognitive architectures. Here, we introduce BrainStat - a toolbox for (i) univariate and multivariate linear models in volumetric and surface-based brain imaging datasets, and (ii) multidomain feature association of results with respect to spatial maps of post-mortem gene expression and histology, task-based fMRI meta-analysis, as well as resting-state fMRI motifs across several common surface templates. The combination of statistics and feature associations into a turnkey toolbox streamlines analytical processes and accelerates cross-modal research. The toolbox is implemented in both Python and MATLAB, two widely used programming languages in the neuroimaging and neuroinformatics communities. BrainStat is openly available and complemented by an expandable documentation.

Insular cortex is a core hub involved in multiple cognitive and socio-affective processes. Yet, the anatomical mechanisms that explain how it is involved in such a diverse array of functions remain incompletely understood. Here, we tested the hypothesis that changes in myeloarchitecture across the insular cortex explain how it can be involved in many different facets of cognitive function. Detailed intracortical profiling, performed across hundreds of insular locations on the basis of myelin-sensitive magnetic resonance imaging (MRI), was compressed into a lower-dimensional space uncovering principal axes of myeloarchitectonic variation. Leveraging two datasets with different high-resolution MRI contrasts, we obtained robust support for two principal dimensions of insular myeloarchitectonic differentiation in vivo, one running from ventral anterior to posterior banks and one radiating from dorsal anterior towards both ventral anterior and posterior subregions. Analyses of post mortem 3D histological data showed that the antero-posterior axis was mirrored in cytoarchitectural markers, even when controlling for sulco-gyral folding. Resting-state functional connectomics in the same individuals and ad hoc meta-analyses showed that myelin gradients in the insula relate to diverse affiliation to macroscale intrinsic functional systems, showing differential shifts in functional network embedding across each myelin-derived gradient. Collectively, our findings offer a novel approach to capture structure-function interactions of a key node of the limbic system, and suggest a multidimensional structural basis underlying the diverse functional roles of the insula.

Multimodal magnetic resonance imaging (MRI) has accelerated human neuroscience by fostering the analysis of brain microstructure, geometry, function, and connectivity across multiple scales and in living brains. The richness and complexity of multimodal neuroimaging, however, demands processing methods to integrate information across modalities and to consolidate findings across different spatial scales. Here, we present micapipe, an open processing pipeline for multimodal MRI datasets. Based on BIDS-conform input data, micapipe can generate i) structural connectomes derived from diffusion tractography, ii) functional connectomes derived from resting-state signal correlations, iii) geodesic distance matrices that quantify cortico-cortical proximity, and iv) microstructural profile covariance matrices that assess inter-regional similarity in cortical myelin proxies. The above matrices can be automatically generated across established 18 cortical parcellations (100-1000 parcels), in addition to subcortical and cerebellar parcellations, allowing researchers to replicate findings easily across different spatial scales. Results are represented on three different surface spaces (native, conte69, fsaverage5), and outputs are BIDS-conform. Processed outputs can be quality controlled at the individual and group level. micapipe was tested on several datasets and is available at https://github.com/MICA-MNI/micapipe, documented at https://micapipe.readthedocs.io/, and containerized as a BIDS App http://bids-apps.neuroimaging.io/apps/. We hope that micapipe will foster robust and integrative studies of human brain microstructure, morphology, function, cand connectivity.

Ongoing brain function is largely determined by the underlying wiring of the brain, but the specific rules governing this relationship remain unknown. Emerging literature has suggested that functional interactions between brain regions emerge from the structural connections through mono- as well as polysynaptic mechanisms. Here, we propose a novel approach based on diffusion maps and Riemannian optimization to emulate this dynamic mechanism in the form of random walks on the structural connectome and predict functional interactions as a weighted combination of these random walks. Our proposed approach was evaluated in two different cohorts of healthy adults (Human Connectome Project, HCP; Microstructure-Informed Connectomics, MICs). Our approach outperformed existing approaches and showed that performance plateaus approximately around the third random walk. At macroscale, we found that the largest number of walks was required in nodes of the default mode and frontoparietal networks, underscoring an increasing relevance of polysynaptic communication mechanisms in transmodal cortical networks compared to primary and unimodal systems.

Human cognition is dynamic, alternating over time between externally-focused states and more abstract, often self-generated, patterns of thought. Although cognitive neuroscience has documented how networks anchor particular modes of brain function, mechanisms that describe transitions between distinct functional states remain poorly understood. Here, we examined how time-varying changes in brain function emerge within the constraints imposed by macroscale structural network organization. Studying a large cohort of healthy adults (n = 326), we capitalized on manifold learning techniques that identify low dimensional representations of structural connectome organization and we decomposed neurophysiological activity into distinct functional states and their transition patterns using Hidden Markov Models. Structural connectome organization predicted dynamic transitions anchored in sensorimotor systems and those between sensorimotor and transmodal states. Connectome topology analyses revealed that transitions involving sensorimotor states traversed short and intermediary distances and adhered strongly to communication mechanisms of network diffusion. Conversely, transitions between transmodal states involved spatially distributed hubs and increasingly engaged long-range routing. These findings establish that the structure of the cortex is optimized to allow neural states the freedom to vary between distinct modes of processing, and so provides a key insight into the neural mechanisms that give rise to the flexibility of human cognition.

Ageing is commonly associated with changes to segregation and integration of functional brain networks, but, in isolation, current network-based approaches struggle to elucidate changes across the many axes of functional organisation. However, the advent of gradient mapping techniques in neuroimaging provides a new means of studying functional organisation in a multi-dimensional connectivity space. Here, we studied ageing and behaviourally-relevant differences in a three-dimensional connectivity space using the Cambridge Centre for Ageing Neuroscience cohort (n = 643). Building on gradient mapping techniques, we developed a set of measures to quantify the dispersion within and between functional communities. We detected a strong shift of the visual network across the adult lifespan from an extreme to a more central position in the 3D gradient space. In contrast, the dispersion distance between transmodal communities (dorsal attention, ventral attention, frontoparietal and default mode) did not change. However, these communities themselves were increasingly dispersed with increasing age, reflecting more dissimilar functional connectivity profiles within each community. Increasing dispersion of frontoparietal, attention and default mode networks, in particular, were associated negatively with cognition, measured by fluid intelligence. By using a technique that explicitly captures the ordering of functional systems in a multi-dimensional hierarchical framework, we identified behaviorally-relevant age-related differences of within and between network organisation. We propose that the study of functional gradients across the adult lifespan could provide insights that may facilitate the development of new strategies to maintain cognitive ability across the lifespan in health and disease.