Psychological states influence our happiness and productivity; however, estimates of their impact have historically been assumed to be limited by the accuracy with which introspection can quantify them. Over the last two decades, studies have shown that introspective descriptions of psychological states correlate with objective indicators of cognition, including task performance and metrics of brain function, using techniques like functional magnetic resonance imaging (fMRI). Such evidence suggests it may be possible to quantify the mapping between self-reports of experience and objective representations of those states (e.g., those inferred from measures of brain activity). Here, we used machine learning to show that self-reported descriptions of experiences across tasks can reliably map the objective landscape of task states derived from brain activity. In our study, 194 participants provided descriptions of their psychological states while performing tasks for which the contribution of different brain systems was available from prior fMRI studies. We used machine learning to combine these reports with descriptions of brain function to form a ‘state-space’ that reliably predicted patterns of brain activity based solely on unseen descriptions of experience (N = 101). Our study demonstrates that introspective reports can share information with the objective task landscape inferred from brain activity. Using machine learning, this study shows that self-reported experience of task states maps onto ‘brain space’, i.e., features of different task states identified in fMRI.
Complex macro-scale patterns of brain activity that emerge during periods of wakeful rest provide insight into the organisation of neural function, how these differentiate individuals based on their traits, and the neural basis of different types of self-generated thoughts. Although brain activity during wakeful rest is valuable for understanding important features of human cognition, its unconstrained nature makes it difficult to disentangle neural features related to personality traits from those related to the thoughts occurring at rest. Our study builds on recent perspectives from work on ongoing conscious thought that highlight the interactions between three brain networks - ventral and dorsal attention networks, as well as the default mode network. We combined measures of personality with state-of-the-art indices of ongoing thoughts at rest and brain imaging analysis and explored whether this ‘tri-partite’ view can provide a framework within which to understand the contribution of states and traits to observed patterns of neural activity at rest. To capture macro-scale relationships between different brain systems, we calculated cortical gradients to describe brain organisation in a low dimensional space. Our analysis established that for more introverted individuals, regions of the ventral attention network were functionally more aligned to regions of the somatomotor system and the default mode network. At the same time, a pattern of detailed self-generated thought was associated with a decoupling of regions of dorsal attention from regions in the default mode network. Our study, therefore, establishes that interactions between attention systems and the default mode network are important influences on ongoing thought at rest and highlights the value of integrating contemporary perspectives on conscious experience when understanding patterns of brain activity at rest.
Recent efforts to chart human brain growth across the lifespan using large-scale MRI data have provided reference standards for human brain development. However, similar models for nonhuman primate (NHP) growth are lacking. The rhesus macaque, a widely used NHP in translational neuroscience due to its similarities in brain anatomy, phylogenetics, cognitive, and social behaviors to humans, serves as an ideal NHP model. This study aimed to create normative growth charts for brain structure across the macaque lifespan, enhancing our understanding of neurodevelopment and aging, and facilitating cross-species translational research. Leveraging data from the PRIMatE Data Exchange (PRIME-DE) and other sources, we aggregated 1,522 MRI scans from 1,024 rhesus macaques. We mapped non-linear developmental trajectories for global and regional brain structural changes in volume, cortical thickness, and surface area over the lifespan. Our findings provided normative charts with centile scores for macaque brain structures and revealed key developmental milestones from prenatal stages to aging, highlighting both species-specific and comparable brain maturation patterns between macaques and humans. The charts offer a valuable resource for future NHP studies, particularly those with small sample sizes. Furthermore, the interactive open resource (https://interspeciesmap.childmind.org) supports cross-species comparisons to advance translational neuroscience research.
Alzheimer's disease (AD) is a brain network disorder where pathological proteins accumulate through networks and drive cognitive decline. Yet, the role of network connectivity in facilitating this accumulation remains unclear. Using in-vivo multimodal imaging, we show that the distribution of tau and reactive microglia in humans follows spatial patterns of connectivity variation, the so-called gradients of brain organization. Notably, less distinct connectivity patterns ("gradient contraction") are associated with cognitive decline in regions with greater tau, suggesting an interaction between reduced network differentiation and tau on cognition. Furthermore, by modeling tau in subject-specific gradient space, we demonstrate that tau accumulation in the frontoparietal and temporo-occipital cortices is associated with greater baseline tau within their functionally and structurally connected hubs, respectively. Our work unveils a role for both functional and structural brain organization in pathology accumulation in AD, and supports subject-specific gradient space as a promising tool to map disease progression.
Background Higher-order cognition is hypothesized to be implemented via distributed cortical networks that are linked via long-range connections. However, it is unknown how computational advantages of long-range connections reflect cortical microstructure and microcircuitry. Methods We investigated this question by (i) profiling long-range cortical connectivity using resting-state functional magnetic resonance imaging (MRI) and cortico-cortical geodesic distance mapping, (ii) assessing how long-range connections reflect local brain microarchitecture, and (iii) examining the microarchitectural similarity of regions connected through long-range connections. Results Analysis of 2 independent datasets indicated that sensory/motor areas had more clustered short-range connections, while transmodal association systems hosted distributed, long-range connections. Meta-analytical decoding suggested that this topographical difference mirrored shifts in cognitive function, from perception/action towards emotional/social processing. Analysis of myelin-sensitive in vivo MRI as well as postmortem histology and transcriptomics datasets established that gradients in functional connectivity distance are paralleled by those present in cortical microarchitecture. Notably, long-range connections were found to link spatially remote regions of association cortex with an unexpectedly similar microarchitecture. Conclusions By mapping covarying topographies of long-range functional connections and cortical microcircuits, the current work provides insights into structure-function relations in human neocortex.
A central goal in neuroscience is the development of a comprehensive mapping between structural and functional brain features, which facilitates mechanistic interpretation of brain function. However, the interpretability of structure-function brain models remains limited by a lack of biological detail. Here, we characterize human structural brain networks weighted by multiple white matter microstructural features including total intra-axonal cross-sectional area and myelin content. We report edge-weight-dependent spatial distributions, variance, small-worldness, rich club, hubs, as well as relationships with function, edge length, and myelin. Contrasting networks weighted by the total intra-axonal cross-sectional area and myelin content of white matter tracts, we find opposite relationships with functional connectivity, an edge-length-independent inverse relationship with each other, and the lack of a canonical rich club in myelin-weighted networks. When controlling for edge length, networks weighted by either fractional anisotropy, radial diffusivity, or neurite density show no relationship with whole-brain functional connectivity. We conclude that the co-utilization of structural networks weighted by total intra-axonal cross-sectional area and myelin content could improve our understanding of the mechanisms mediating the structure-function brain relationship.
The wiring of the brain connects micro-architecturally diverse neuronal populations. The conventional graph model encodes macroscale brain connectivity as a network of nodes and edges, but abstracts away the rich biological detail of each regional node. Regions are different in terms of their microscale attributes, many of which are readily available through modern technological advances and data-sharing initiatives. How is macroscale connectivity related to nodal attributes? Here we investigate the systematic arrangement of white-matter connectivity with respect to multiple biological annotations. Namely, we formally study assortative mixing in annotated connectomes by quantifying the tendency for regions to be connected with each other based on the similarity of their micro-architectural attributes. We perform all experiments using four cortico-cortical connectome datasets from three different species (human, macaque and mouse), and consider a range of molecular, cellular and laminar annotations, including gene expression, neurotransmitter receptors, neuron density, laminar thickness and intracortical myelin. Importantly, we disentangle the relationship between neural wiring, regional heterogeneity and spatial embedding using spatial autocorrelation-preserving null models. We show that mixing between micro-architecturally diverse neuronal populations is supported by long-distance connections. Using meta-analytic decoding, we find that the arrangement of connectivity patterns with respect to biological annotations shape patterns of regional functional specialization. Specifically, regions that connect to biologically similar regions are associated with executive function; conversely, regions that connect with biologically dissimilar regions are associated with memory function. By bridging scales of cortical organization, from microscale attributes to macroscale connectivity, this work lays the foundation for next-generation annotated connectomics.
Background Autism spectrum disorder (ASD) is a common neurodevelopmental diagnosis showing substantial phenotypic heterogeneity. A leading example can be found in verbal and nonverbal cognitive skills, which vary from elevated to impaired compared with neurotypical individuals. Moreover, deficits in verbal profiles often coexist with normal or superior performance in the nonverbal domain. Methods To study brain substrates underlying cognitive imbalance in ASD, we capitalized categorical and dimensional IQ profiling as well as multimodal neuroimaging. Results IQ analyses revealed a marked verbal to nonverbal IQ imbalance in ASD across 2 datasets (Dataset-1: 155 ASD, 151 controls; Dataset-2: 270 ASD, 490 controls). Neuroimaging analysis in Dataset-1 revealed a structure-function substrate of cognitive imbalance, characterized by atypical cortical thickening and altered functional integration of language networks alongside sensory and higher cognitive areas. Conclusion Although verbal and nonverbal intelligence have been considered as specifiers unrelated to autism diagnosis, our results indicate that intelligence disparities are accentuated in ASD and reflected by a consistent structure-function substrate affecting multiple brain networks. Our findings motivate the incorporation of cognitive imbalances in future autism research, which may help to parse the phenotypic heterogeneity and inform intervention-oriented subtyping in ASD.
Recent theories of cortical organisation suggest features of function emerge from the spatial arrangement of brain regions. For example, association cortex is located furthest from systems involved in action and perception. Association cortex is also ‘interdigitated’ with adjacent regions having different patterns of functional connectivity. It is assumed that topographic properties, such as distance between regions, constrains their functions, however, we lack a formal description of how this occurs. Here we use variograms, a quantification of spatial autocorrelation, to profile how function changes with the distance between cortical regions. We find function changes with distance more gradually within sensory-motor cortex than association cortex. Importantly, systems within the same type of cortex (e.g., fronto-parietal and default mode networks) have similar profiles. Primary and association cortex, therefore, are differentiated by how function changes over space, emphasising the value of topographical features of a region when estimating its contribution to cognition and behaviour.
Objective Temporal lobe epilepsy (TLE) is the most common drug-resistant epilepsy in adults. Here, we aimed to profile local neural function in TLE in vivo , building on prior evidence that has identified widespread structural alterations. Using multimodal MRI, we mapped intrinsic neural timescales (INT) at rest, examined associations to TLE-related structural compromise, and evaluated the clinical utility of INT. Methods We studied 46 TLE patients and 44 healthy controls from two independent sites, and mapped INT changes in patients relative to controls across hippocampal, subcortical, and neocortical regions. We examined region-specific associations to structural alterations and explored effects of age and epilepsy duration. A supervised machine learning paradigm assessed utility of INT for classifying patients- vs -controls and seizure focus lateralization. Results Relative to controls, TLE showed marked INT reductions across multiple regions bilaterally, indexing faster changing resting activity, with strongest effects in ipsilateral medial and lateral temporal regions, and sensorimotor cortices. Findings were consistent in each site and robust, albeit with reduced effect sizes, when correcting for structural alterations. TLE-related INT reductions increased with advancing disease duration, yet findings differed from aging effects seen in controls. Classifiers based on INT distinguished patients- vs -controls (balanced accuracy, 5-fold: 76±2.65%; cross-site, 72-83%) and lateralized the focus in TLE (balanced accuracy, 5-fold: 96±2.10%; cross-site, 95-97%) with high accuracies and generalization. Conclusions Our findings robustly demonstrate atypical macroscale function in TLE in a topography that extends beyond mesiotemporal epicenters. INT measurements can assist in TLE diagnosis, seizure focus lateralization, and monitoring of disease progression, which suggests clinical utility.
Tau pathology can spread through connectivity-based networks, with certain regions (or epicenters) accumulating more tau than others. Such spatial vulnerabilities may be due to their unique apical position in the cortical hierarchy, which can be elucidated through ‘gradients of connectivity’ (Margulies 2016 PNAS). Prior work showed that the primary gradient of functional connectivity unveils a uni-to-transmodal topography of the healthy neocortex which highly correlates with a cognitive gradient of perception-to-abstraction. Here, we hypothesized that the gradients of functional/structural connectivity interact with tau to affect cognitive functions in Alzheimer’s disease (AD). We included 213 participants from TRIAD (103 CN Aß-, 103 CN Aß+, and 75 CI Aß+) with diffusion-weighted MRI, resting-state functional MRI, 18F-MK6240 tau-PET, and an extensive cognitive battery. We performed non-linear dimensionality reduction on the individual functional and structural connectomes, and extracted the first components (‘gradients’) -explaining most variance (G1FC and G1SC). First, we compared G1FC_or_SC between diagnostic groups. Second, we investigated the interaction effect of G1FC_or_SC*tauSUVR on cognition (across 9 cognitive domains). Last, we investigated whether the tau-cognition relationship changed in a topography-specific manner along the cortical hierarchy, within (equally-sized) gradient-derived meta-ROIs along G1FC_or_SC. Results were compared to Braak-derived regional associations. Analyses were adjusted for age, sex, APOE, education, and multiple comparisons. We observed reduced segregation of functional networks in AD compared to controls, with unimodal (lower-order cognitive) and transmodal (higher-order cognitive) regions moving closer on G1FC. This may indicate loss of network specialization in AD. Participants who had both higher tau and G1FC alterations had more cognitive impairment (Fig.1A; shown for MMSE/language). This interaction-effect was less pronounced with G1SC (Fig.1B). Last, tau correlated with cognition in a topography-specific progressive manner (i.e., along the transmodal-unimodal G1FC axis and anterior-posterior G1SC axis) (Fig.1C). Notably, tau correlated with delayed memory progressively along the posterior-anterior G1SC axis (R2 = 0.93 in all and R2 = 0.95 in A+) and BraakI-VI axis (R2 = 0.77 in all and R2 = 0.28 in A+). Our work supports the contribution of connectome-driven tau distribution on cognitive impairment in AD. Connectome gradients may provide a spatial framework to study tau spreading along the major axes of brain organization underlying specific cognitive domains.
Brain imaging research enjoys increasing adoption of supervised machine learning for single-participant disease classification. Yet, the success of these algorithms likely depends on population diversity, including demographic differences and other factors that may be outside of primary scientific interest. Here, we capitalize on propensity scores as a composite confound index to quantify diversity due to major sources of population variation. We delineate the impact of population heterogeneity on the predictive accuracy and pattern stability in 2 separate clinical cohorts: the Autism Brain Imaging Data Exchange (ABIDE, n = 297) and the Healthy Brain Network (HBN, n = 551). Across various analysis scenarios, our results uncover the extent to which cross-validated prediction performances are interlocked with diversity. The instability of extracted brain patterns attributable to diversity is located preferentially in regions part of the default mode network. Collectively, our findings highlight the limitations of prevailing deconfounding practices in mitigating the full consequences of population diversity.
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.
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.
While reading, our mind can wander to unrelated autobiographical information, creating a perceptually decoupled state detrimental to narrative comprehension. To understand how this mind-wandering state emerges, we asked whether retrieving autobiographical content necessitates functional disengagement from visual input. In Experiment 1, brain activity was recorded using functional magnetic resonance imaging (fMRI) in an experimental situation mimicking naturally occurring mind-wandering, allowing us to precisely delineate neural regions involved in memory and reading. Individuals read expository texts and ignored personally relevant autobiographical memories, as well as the opposite situation. Medial regions of the default mode network (DMN) were recruited during memory retrieval. In contrast, left temporal and lateral prefrontal regions of the DMN, as well as ventral visual cortex, were recruited when reading for comprehension. Experiment two used functional connectivity both at rest and during tasks to establish that (i) DMN regions linked to memory are more functionally decoupled from regions of ventral visual cortex than regions in the same network engaged when reading; and (ii) individuals with more self-generated mental contents and poorer comprehension, while reading in the lab, showed more decoupling between visually connected DMN sites important for reading and primary visual cortex. A similar pattern of connectivity was found in Experiment 1, with greater coupling between this DMN site and visual cortex when participants reported greater focus on reading in the face of conflict from autobiographical memory cues; moreover, the retrieval of personally relevant memories increased the decoupling of these sites. These converging data suggest we lose track of the narrative when our minds wander because generating autobiographical mental content relies on cortical regions within the DMN which are functionally decoupled from ventral visual regions engaged during reading.
The hippocampus is a uniquely infolded allocortical structure in the medial temporal lobe that consists of the microstructurally and functionally distinct subregions: subiculum, cornu ammonis, and dentate gyrus. The hippocampus is a remarkably plastic region that is implicated in learning and memory. At the same time it has been shown that hippocampal subregion volumes are heritable, and that genetic expression varies along a posterior to anterior axis. Here, we studied how a heritable, stable, hippocampal organisation may support its flexible function in healthy adults. Leveraging the twin set-up of the Human Connectome Project with multimodal neuroimaging, we observed that the functional connectivity between hippocampus and cortex was heritable and that microstructure of the hippocampus genetically correlated with cortical microstructure. Moreover, both functional and microstructural organisation could be consistently captured by anterior-to-posterior and medial-to-lateral axes across individuals. However, heritability of functional, relative to microstructural, organisation was found reduced, suggesting individual variation in functional organisation may be explained by experience-driven factors. Last, we demonstrate that structure and function couple along an inherited macroscale organisation, suggesting an interplay of stability and plasticity within the hippocampus. Our study provides new insights on the heritability of the hippocampal of the structure and function within the hippocampal organisation.
Multimodal neuroimaging grants a powerful window into the structure and function of the human brain at multiple scales. Recent methodological and conceptual advances have enabled investigations of the interplay between large-scale spatial trends (also referred to as gradients) in brain microstructure and connectivity, offering an integrative framework to study multiscale brain organization. Here, we share a multimodal MRI dataset for Microstructure-Informed Connectomics (MICA-MICs) acquired in 50 healthy adults (23 women; 29.54 ± 5.62 years) who underwent high-resolution T1-weighted MRI, myelin-sensitive quantitative T1 relaxometry, diffusion-weighted MRI, and resting-state functional MRI at 3 Tesla. In addition to raw anonymized MRI data, this release includes brain-wide connectomes derived from (i) resting-state functional imaging, (ii) diffusion tractography, (iii) microstructure covariance analysis, and (iv) geodesic cortical distance, gathered across multiple parcellation scales. Alongside, we share large-scale gradients estimated from each modality and parcellation scale. Our dataset will facilitate future research examining the coupling between brain microstructure, connectivity, and function. MICA-MICs is available on the Canadian Open Neuroscience Platform data portal ( https://portal.conp.ca ) and the Open Science Framework ( https://osf.io/j532r/ ).
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.
The spread of abnormal tau along connectivity-based networks is a hallmark of Alzheimer’s disease (AD). In addition, activated microglia was shown to affect tau spread across networks. However, traditional network approaches using atlas-defined regions suffer from assuming within-region homogenous connectivity, as well as abrupt and linear connectivity changes between the regions; while, in fact, brain structure and function exist as complex overlapping axes of connectivity variation. Here, we studied smooth/continuous spatial transitions or ‘gradients’ of connectivity and their link to gradients of pathology in AD. This work provides novel insights into brain architecture driving pathology spread. We enrolled 220 individuals including 106 cognitively normal Aβ-negative (A-) and 37 A+, and 77 MCI/AD A+ (TRIAD). From diffusion-weighted MRI and fMRI data, we estimated region-to-region fiber count and correlated timeseries, respectively, and averaged subject-wise connectivity matrices to create structural and functional ‘template connectomes’. From 18 F-MK6240 and 11 C-PBR28 PET data, we extracted regional tau and neuroinflammation SUVR values, respectively, and created covariance matrices for each disease stage. Next, the template connectomes/covariance matrices were subjected to non-linear dimensionality reduction to extract embedded components (‘gradients’). We correlated the primary gradients of functional or structural connectivity (G1 FC /G1 SC ) with the primary and secondary gradients of tau or neuroinflammation (G1 TAU /G2 TAU ; G1 INFLAM /G2 INFLAM ). We employed three brain atlases for validation: high-resolution customized Glasser-atlas (1332 equally-sized subregions), DKT, and Schaefer-100. G1 FC replicated the previously reported ‘uni-to-transmodal’ functional gradient ( Margulies 2016) and correlated strongly with G1 TAU in cognitively impaired A+ ( Fig.1 A); illustrating a common topographic variation in tau spread and functional connectivity in symptomatic stages. Conversely, G1 SC ran between medial-temporal-lobe (MTL)/posterior and frontal/anterior regions, and correlated strongly with both G2 TAU ( Fig.1 B) and G1 INFLAM ( Fig.1 C) in cognitively (un)impaired A+; illustrating a common topographic variation in tau/microglia and structural connectivity already present at early disease stages. Different brain atlases led to similar results. We identified two gradients of continuous topographic variation in tau and neuroinflammation, potentially reflecting heterogeneity of AD-subgroups. While microglial activation and early-stage tau propagation may be largely dominated by structural organisation, later-stage tau propagation may be more heterogeneously associated with variation in structural and particularly functional connectivity.
Semantic cognition allows us to make sense of our varied experiences, including the words we hear and the objects we see. Contemporary accounts identify multiple interacting components that underpin semantic cognition, including diverse unimodal "spoke" systems that are integrated by a heteromodal "hub", and control processes that allow us to access weakly-encoded as well as dominant aspects of knowledge to suit the circumstances. The current study examined how these dimensions of semantic cognition might be related to whole-brain-derived components (or gradients) of connectivity. A nonlinear dimensionality reduction technique was applied to resting-state functional magnetic resonance imaging from 176 participants to characterise the strength of two key connectivity gradients in each individual: the principal gradient captured the separation between unimodal and heteromodal cortex, while the second gradient corresponded to the distinction between motor and visual cortex. We then examined whether the magnitude of these gradients within the semantic network was related to specific aspects of semantic cognition by examining individual differences in semantic and non-semantic tasks. Participants whose intrinsic connectivity showed a better fit with Gradient 1 had faster identification of weak semantic associations. Furthermore, a better fit with Gradient 2 was linked to faster performance on picture semantic judgements. These findings show that individual differences in aspects of semantic cognition can be related to components of connectivity within the semantic network.