Long-range functional connections mirror and link microarchitectural and cognitive hierarchies in the human brain

Transcriptomic contributions to a modern cytoarchitectonic parcellation of the human cerebral cortex

Transcriptomic contributions to the anatomical, functional, and network layout of the human cerebral cortex (HCC) have become a major interest in cognitive and systems neuroscience. Here, we tested if transcriptomic differences support a modern, algorithmic cytoarchitectonic parcellation of HCC. Using a data-driven approach, we identified a sparse subset of genes that differentially contributed to the cytoarchitectonic parcellation of HCC. A combined metric of cortical thickness and myelination (CT/M ratio), as well as cell density, correlated with gene expression. Enrichment analyses showed that genes specific to the cytoarchitectonic parcellation of the HCC were related to molecular functions such as transmembrane transport and ion channel activity. Together, the relationship between transcriptomics and cytoarchitecture bridges the gap among (i) gradients at the macro-scale (including thickness and myelination), (ii) areas at the meso-scale, and (iii) cell density at the microscale, as well as supports the recently proposed cortical spectrum theory and structural model.

Atypical connectome topography and signal flow in temporal lobe epilepsy

Temporal lobe epilepsy (TLE) is the most common pharmaco-resistant epilepsy in adults. While primarily associated with mesiotemporal pathology, recent evidence suggests that brain alterations in TLE extend beyond the paralimbic epicenter and impact macroscale function and cognitive functions, particularly memory. Using connectome-wide manifold learning and generative models of effective connectivity, we examined functional topography and directional signal flow patterns between large-scale neural circuits in TLE at rest. Studying a multisite cohort of 95 patients with TLE and 95 healthy controls, we observed atypical functional topographies in the former group, characterized by reduced differentiation between sensory and transmodal association cortices, with most marked effects in bilateral temporo-limbic and ventromedial prefrontal cortices. These findings were consistent across all study sites, present in left and right lateralized patients, and validated in a subgroup of patients with histopathological validation of mesiotemporal sclerosis and post-surgical seizure freedom. Moreover, they were replicated in an independent cohort of 30 TLE patients and 40 healthy controls. Further analyses demonstrated that reduced differentiation related to decreased functional signal flow into and out of temporolimbic cortical systems and other brain networks. Parallel analyses of structural and diffusion-weighted MRI data revealed that topographic alterations were independent of TLE-related cortical thinning but partially mediated by white matter microstructural changes that radiated away from paralimbic circuits. Finally, we found a strong association between the degree of functional alterations and behavioral markers of memory dysfunction. Our work illustrates the complex landscape of macroscale functional imbalances in TLE, which can serve as intermediate markers bridging microstructural changes and cognitive impairment.

Cytoarchitectonic gradients of laminar degeneration in behavioral variant frontotemporal dementia

Behavioral variant frontotemporal dementia (bvFTD) is a clinical syndrome primarily caused by either tau (bvFTD-tau) or TDP-43 (bvFTD-TDP) proteinopathies. We previously found lower cortical layers and dorsolateral regions accumulate greater tau than TDP-43 pathology; however, patterns of laminar neurodegeneration across diverse cytoarchitecture in bvFTD is understudied. We hypothesized that bvFTD-tau and bvFTD-TDP have distinct laminar distributions of pyramidal neurodegeneration along cortical gradients, a topologic order of cytoarchitectonic subregions based on increasing pyramidal density and laminar differentiation. Here, we tested this hypothesis in a frontal cortical gradient consisting of five cytoarchitectonic types (i.e., periallocortex, agranular mesocortex, dysgranular mesocortex, eulaminate-I isocortex, eulaminate-II isocortex) spanning anterior cingulate, paracingulate, orbitofrontal, and mid-frontal gyri in bvFTD-tau (n=27), bvFTD-TDP (n=47), and healthy controls (HC; n=32). We immunostained all tissue for total neurons (NeuN; neuronal-nuclear protein) and pyramidal neurons (SMI32; non-phosphorylated neurofilament) and digitally quantified NeuN-immunoreactivity (ir) and SMI32-ir in supragranular II-III, infragranular V-VI, and all I-VI layers in each cytoarchitectonic type. We used linear mixed-effects models adjusted for demographic and biologic variables to compare SMI32-ir between groups and examine relationships with the cortical gradient, long-range pathways, and clinical symptoms. We found regional and laminar distributions of SMI32-ir expected for HC, validating our measures within the cortical gradient framework. While SMI32-ir loss was not related to the cortical gradient in bvFTD-TDP, SMI32-ir progressively decreased along the cortical gradient of bvFTD-tau and included greater SMI32-ir loss in supragranular eulaminate-II isocortex in bvFTD-tau vs bvFTD-TDP ( p =0.039). In a structural model for long-range laminar connectivity between infragranular mesocortex and supragranular isocortex, we found a larger laminar ratio of mesocortex-to-isocortex SMI32-ir in bvFTD-tau vs bvFTD-TDP ( p =0.019), suggesting select long-projecting pathways may contribute to isocortical-predominant degeneration in bvFTD-tau. In cytoarchitectonic types with the highest NeuN-ir, we found lower SMI32-ir in bvFTD-tau vs bvFTD-TDP ( p =0.047), suggesting pyramidal neurodegeneration may occur earlier in bvFTD-tau. Lastly, we found that reduced SMI32-ir related to behavioral severity and frontal-mediated letter fluency, not temporal-mediated confrontation naming, demonstrating the clinical relevance and specificity of frontal pyramidal neurodegeneration to bvFTD-related symptoms. Our data suggest loss of neurofilament-rich pyramidal neurons is a clinically relevant feature of bvFTD that selectively worsens along a frontal cortical gradient in bvFTD-tau, not bvFTD-TDP. Therefore, tau-mediated degeneration may preferentially involve pyramidal-rich layers that connect more distant cytoarchitectonic types. Moreover, the hierarchical arrangement of cytoarchitecture along cortical gradients may be an important neuroanatomical framework for identifying which types of cells and pathways are differentially involved between proteinopathies.

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

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

Methods for decoding cortical gradients of functional connectivity

Macroscale gradients have emerged as a central principle for understanding functional brain organization. Previous studies have demonstrated that a principal gradient of connectivity in the human brain exists, with unimodal primary sensorimotor regions situated at one end and transmodal regions associated with the default mode network and representative of abstract functioning at the other. The functional significance and interpretation of macroscale gradients remains a central topic of discussion in the neuroimaging community, with some studies demonstrating that gradients may be described using meta-analytic functional decoding techniques. However, additional methodological development is necessary to fully leverage available meta-analytic methods and resources and quantitatively evaluate their relative performance. Here, we conducted a comprehensive series of analyses to investigate and improve the framework of data-driven, meta-analytic methods, thereby establishing a principled approach for gradient segmentation and functional decoding. We found that a two-segment solution determined by a k-means segmentation approach and an LDA-based meta-analysis combined with the NeuroQuery database was the optimal combination of methods for decoding functional connectivity gradients. Finally, we proposed a method for decoding additional components of the gradient decomposition. The current work aims to provide recommendations on best practices and flexible methods for gradient-based functional decoding of fMRI data.

Sex differences in intrinsic functional cortical organization reflect differences in network topology rather than cortical morphometry

Brain size robustly differs between sexes. However, the consequences of this anatomical dimorphism on sex differences in intrinsic brain function remain unclear. We investigated the extent to which sex differences in intrinsic cortical functional organization may be explained by differences in cortical morphometry, namely brain size, microstructure, and the geodesic distances of connectivity profiles. For this, we computed a low dimensional representation of functional cortical organization, the sensory-association axis, and identified widespread sex differences. Contrary to our expectations, observed sex differences in functional organization were not fundamentally associated with differences in brain size, microstructural organization, or geodesic distances, despite these morphometric properties being per se associated with functional organization and differing between sexes. Instead, functional sex differences in the sensory-association axis were associated with differences in functional connectivity profiles and network topology. Collectively, our findings suggest that sex differences in functional cortical organization extend beyond sex differences in cortical morphometry.

Teaser

Investigating sex differences in functional cortical organization and their association to differences in cortical morphometry.

C-SMB 2.0: Integrating over 25 years of motor sequencing research with the Discrete Sequence Production task

An exhaustive review is reported of over 25 years of research with the Discrete Sequence Production (DSP) task as reported in well over 100 articles. In line with the increasing call for theory development, this culminates into proposing the second version of the Cognitive framework of Sequential Motor Behavior (C-SMB 2.0), which brings together known models from cognitive psychology, cognitive neuroscience, and motor learning. This processing framework accounts for the many different behavioral results obtained with the DSP task and unveils important properties of the cognitive system. C-SMB 2.0 assumes that a versatile central processor (CP) develops multimodal, central-symbolic representations of short motor segments by repeatedly storing the elements of these segments in short-term memory (STM). Independently, the repeated processing by modality-specific perceptual and motor processors (PPs and MPs) and by the CP when executing sequences gradually associates successively used representations at each processing level. The high dependency of these representations on active context information allows for the rapid serial activation of the sequence elements as well as for the executive control of tasks as a whole. Speculations are eventually offered as to how the various cognitive processes could plausibly find their neural underpinnings within the intricate networks of the brain.

Diverging asymmetry of intrinsic functional organization in autism

Autism is a neurodevelopmental condition involving atypical sensory-perceptual functions together with language and socio-cognitive deficits. Previous work has reported subtle alterations in the asymmetry of brain structure and reduced laterality of functional activation in individuals with autism relative to non-autistic individuals (NAI). However, whether functional asymmetries show altered intrinsic systematic organization in autism remains unclear. Here, we examined inter- and intra-hemispheric asymmetry of intrinsic functional gradients capturing connectome organization along three axes, stretching between sensory-default, somatomotor-visual, and default-multiple demand networks, to study system-level hemispheric imbalances in autism. We observed decreased leftward functional asymmetry of language network organization in individuals with autism, relative to NAI. Whereas language network asymmetry varied across age groups in NAI, this was not the case in autism, suggesting atypical functional laterality in autism may result from altered developmental trajectories. Finally, we observed that intra- but not inter-hemispheric features were predictive of the severity of autistic traits. Our findings illustrate how regional and patterned functional lateralization is altered in autism at the system level. Such differences may be rooted in atypical developmental trajectories of functional organization asymmetry in autism.

Assortative mixing in micro-architecturally annotated brain connectomes

The wiring of the brain connects micro-architecturally diverse neuronal populations, but the conventional graph model, which encodes macroscale brain connectivity as a network of nodes and edges, abstracts away the rich biological detail of each regional node. Here, we annotate connectomes with multiple biological attributes and formally study assortative mixing in annotated connectomes. Namely, we quantify the tendency for regions to be connected based on the similarity of their micro-architectural attributes. We perform all experiments using four cortico-cortical connectome datasets from three different species, and consider a range of molecular, cellular, and laminar annotations. We show that mixing between micro-architecturally diverse neuronal populations is supported by long-distance connections and find that the arrangement of connections with respect to biological annotations is associated to patterns of regional functional specialization. By bridging scales of cortical organization, from microscale attributes to macroscale connectivity, this work lays the foundation for next-generation annotated connectomics.

Network evolution of regional brain volumes in young children reflects neurocognitive scores and mother's education

The maturation of regional brain volumes from birth to preadolescence is a critical developmental process that underlies emerging brain structural connectivity and function. Regulated by genes and environment, the coordinated growth of different brain regions plays an important role in cognitive development. Current knowledge about structural network evolution is limited, partly due to the sparse and irregular nature of most longitudinal neuroimaging data. In particular, it is unknown how factors such as mother's education or sex of the child impact the structural network evolution. To address this issue, we propose a method to construct evolving structural networks and study how the evolving connections among brain regions as reflected at the network level are related to maternal education and biological sex of the child and also how they are associated with cognitive development. Our methodology is based on applying local Fréchet regression to longitudinal neuroimaging data acquired from the RESONANCE cohort, a cohort of healthy children (245 females and 309 males) ranging in age from 9 weeks to 10 years. Our findings reveal that sustained highly coordinated volume growth across brain regions is associated with lower maternal education and lower cognitive development. This suggests that higher neurocognitive performance levels in children are associated with increased variability of regional growth patterns as children age.

BrainStat: A toolbox for brain-wide statistics and multimodal feature associations

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.

Micapipe: A pipeline for multimodal neuroimaging and connectome analysis

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.

Brain morphometric similarity and flexibility

Background

The cerebral cortex is represented through multiple multilayer morphometric similarity networks to study their modular structures. The approach introduces a novel way for studying brain networks' metrics across individuals, and can quantify network properties usually not revealed using conventional network analyses.

Methods

A total of 8 combinations or types of morphometric similarity networks were constructed – 4 combinations of the inter-regional cortical features on 2 brain atlases. The networks' modular structures were investigated by identifying those modular interactions that stay consistent across the combinations of inter-regional morphometric features and individuals.

Results

The results provide evidence of the community structures as the property of (i) cortical lobar divisions, and also as (ii) the product of different combinations of morphometric features used for the construction of the multilayer representations of the cortex. For the first time, this study has mapped out flexible and inflexible morphometric similarity hubs, and evidence has been provided about variations of the modular network topology across the multilayers with age and IQ.

Conclusions

The results contribute to understanding of intra-regional characteristics in cortical interactions, which potentially can be used to map heterogeneous neurodegeneration patterns in diseased brains.