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.

Heritability of hippocampal functional and microstructural organisation.

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.

Automated individual-level parcellation of Broca's region based on functional connectivity.

Broca's region can be subdivided into its constituent areas 44 and 45 based on established differences in connectivity to superior temporal and inferior parietal regions. The current study builds on our previous work manually parcellating Broca's area on the individual-level by applying these anatomical criteria to functional connectivity data. Here we present an automated observer-independent and anatomy-informed parcellation pipeline with comparable precision to the manual labels at the individual-level. The method first extracts individualized connectivity templates of areas 44 and 45 by assigning to each surface vertex within the ventrolateral frontal cortex the partial correlation value of its functional connectivity to group-level templates of areas 44 and 45, accounting for other template connectivity patterns. To account for cross-subject variability in connectivity, the partial correlation procedure is then repeated using individual-level network templates, including individual-level connectivity from areas 44 and 45. Each node is finally labeled as area 44, 45, or neither, using a winner-take-all approach. The method also incorporates prior knowledge of anatomical location by weighting the results using spatial probability maps. The resulting area labels show a high degree of spatial overlap with the gold-standard manual labels, and group-average area maps are consistent with cytoarchitectonic probability maps of areas 44 and 45. To facilitate reproducibility and to demonstrate that the method can be applied to resting-state fMRI datasets with varying acquisition and preprocessing parameters, the labeling procedure is applied to two open-source datasets from the Human Connectome Project and the Nathan Kline Institute Rockland Sample. While the current study focuses on Broca's region, the method is adaptable to parcellate other cortical regions with distinct connectivity profiles.

Relating sex-bias in human cortical and hippocampal microstructure to sex hormones.

Determining sex-bias in brain structure is of great societal interest to improve diagnostics and treatment of brain-related disorders. So far, studies on sex-bias in brain structure predominantly focus on macro-scale measures, and often ignore factors determining this bias. Here we study sex-bias in cortical and hippocampal microstructure in relation to sex hormones. Investigating quantitative intracortical profiling in-vivo using the T1w/T2w ratio in 1093 healthy females and males of the cross-sectional Human Connectome Project young adult sample, we find that regional cortical and hippocampal microstructure differs between males and females and that the effect size of this sex-bias varies depending on self-reported hormonal status in females. Microstructural sex-bias and expression of sex hormone genes, based on an independent post-mortem sample, are spatially coupled. Lastly, sex-bias is most pronounced in paralimbic areas, with low laminar complexity, which are predicted to be most plastic based on their cytoarchitectural properties. Albeit correlative, our study underscores the importance of incorporating sex hormone variables into the investigation of brain structure and plasticity.

Macroscale Thalamic Functional Organization Disturbances and Underlying Core Cytoarchitecture in Early-Onset Schizophrenia.

BACKGROUND AND HYPOTHESIS: Schizophrenia is a polygenetic mental disorder with heterogeneous positive and negative symptom constellations, and is associated with abnormal cortical connectivity. The thalamus has a coordinative role in cortical function and is key to the development of the cerebral cortex. Conversely, altered functional organization of the thalamus might relate to overarching cortical disruptions in schizophrenia, anchored in development. STUDY DESIGN: Here, we contrasted resting-state fMRI in 86 antipsychotic-naive first-episode early-onset schizophrenia (EOS) patients and 91 typically developing controls to study whether macroscale thalamic organization is altered in EOS. Employing dimensional reduction techniques on thalamocortical functional connectome (FC), we derived lateral-medial and anterior-posterior thalamic functional axes. STUDY RESULTS: We observed increased segregation of macroscale thalamic functional organization in EOS patients, which was related to altered thalamocortical interactions both in unimodal and transmodal networks. Using an ex vivo approximation of core-matrix cell distribution, we found that core cells particularly underlie the macroscale abnormalities in EOS patients. Moreover, the disruptions were associated with schizophrenia-related gene expression maps. Behavioral and disorder decoding analyses indicated that the macroscale hierarchy disturbances might perturb both perceptual and abstract cognitive functions and contribute to negative syndromes in patients. CONCLUSIONS: These findings provide mechanistic evidence for disrupted thalamocortical system in schizophrenia, suggesting a unitary pathophysiological framework.

Heritability and cross-species comparisons of human cortical functional organization asymmetry.

The human cerebral cortex is symmetrically organized along large-scale axes but also presents inter-hemispheric differences in structure and function. The quantified contralateral homologous difference, that is asymmetry, is a key feature of the human brain left-right axis supporting functional processes, such as language. Here, we assessed whether the asymmetry of cortical functional organization is heritable and phylogenetically conserved between humans and macaques. Our findings indicate asymmetric organization along an axis describing a functional trajectory from perceptual/action to abstract cognition. Whereas language network showed leftward asymmetric organization, frontoparietal network showed rightward asymmetric organization in humans. These asymmetries were heritable in humans and showed a similar spatial distribution with macaques, in the case of intra-hemispheric asymmetry of functional hierarchy. This suggests (phylo)genetic conservation. However, both language and frontoparietal networks showed a qualitatively larger asymmetry in humans relative to macaques. Overall, our findings suggest a genetic basis for asymmetry in intrinsic functional organization, linked to higher order cognitive functions uniquely developed in humans.

Genetic and phylogenetic uncoupling of structure and function in human transmodal cortex.

Brain structure scaffolds intrinsic function, supporting cognition and ultimately behavioral flexibility. However, it remains unclear how a static, genetically controlled architecture supports flexible cognition and behavior. Here, we synthesize genetic, phylogenetic and cognitive analyses to understand how the macroscale organization of structure-function coupling across the cortex can inform its role in cognition. In humans, structure-function coupling was highest in regions of unimodal cortex and lowest in transmodal cortex, a pattern that was mirrored by a reduced alignment with heritable connectivity profiles. Structure-function uncoupling in macaques had a similar spatial distribution, but we observed an increased coupling between structure and function in association cortices relative to humans. Meta-analysis suggested regions with the least genetic control (low heritable correspondence and different across primates) are linked to social-cognition and autobiographical memory. Our findings suggest that genetic and evolutionary uncoupling of structure and function in different transmodal systems may support the emergence of complex forms of cognition.

The impact of ischemic stroke on connectivity gradients.

The functional organization of the brain can be represented as a low-dimensional space that reflects its macroscale hierarchy. The dimensions of this space, described as connectivity gradients, capture the similarity of areas' connections along a continuous space. Studying how pathological perturbations with known effects on functional connectivity affect these connectivity gradients provides support for their biological relevance. Previous work has shown that localized lesions cause widespread functional connectivity alterations in structurally intact areas, affecting a network of interconnected regions. By using acute stroke as a model of the effects of focal lesions on the connectome, we apply the connectivity gradient framework to depict how functional reorganization occurs throughout the brain, unrestricted by traditional definitions of functional network boundaries. We define a three-dimensional connectivity space template based on functional connectivity data from healthy controls. By projecting lesion locations into this space, we demonstrate that ischemic strokes result in dimension-specific alterations in functional connectivity over the first week after symptom onset. Specifically, changes in functional connectivity were captured along connectivity Gradients 1 and 3. The degree of functional connectivity change was associated with the distance from the lesion along these connectivity gradients (a measure of functional similarity) regardless of the anatomical distance from the lesion. Together, these results provide support for the biological validity of connectivity gradients and suggest a novel framework to characterize connectivity alterations after stroke.