A Common Space Approach to Comparative Neuroscience.

Comparative neuroscience is entering the era of big data. New high-throughput methods and data-sharing initiatives have resulted in the availability of large, digital data sets containing many types of data from ever more species. Here, we present a framework for exploiting the new possibilities offered. The multimodality of the data allows vertical translations, which are comparisons of different aspects of brain organization within a single species and across scales. Horizontal translations compare particular aspects of brain organization across species, often by building abstract feature spaces. Combining vertical and horizontal translations allows for more sophisticated comparisons, including relating principles of brain organization across species by contrasting horizontal translations, and for making formal predictions of unobtainable data based on observed results in a model species.

SARS-CoV-2 is associated with changes in brain structure in UK Biobank.

There is strong evidence of brain-related abnormalities in COVID-191-13. However, it remains unknown whether the impact of SARS-CoV-2 infection can be detected in milder cases, and whether this can reveal possible mechanisms contributing to brain pathology. Here we investigated brain changes in 785 participants of UK Biobank (aged 51-81 years) who were imaged twice using magnetic resonance imaging, including 401 cases who tested positive for infection with SARS-CoV-2 between their two scans-with 141 days on average separating their diagnosis and the second scan-as well as 384 controls. The availability of pre-infection imaging data reduces the likelihood of pre-existing risk factors being misinterpreted as disease effects. We identified significant longitudinal effects when comparing the two groups, including (1) a greater reduction in grey matter thickness and tissue contrast in the orbitofrontal cortex and parahippocampal gyrus; (2) greater changes in markers of tissue damage in regions that are functionally connected to the primary olfactory cortex; and (3) a greater reduction in global brain size in the SARS-CoV-2 cases. The participants who were infected with SARS-CoV-2 also showed on average a greater cognitive decline between the two time points. Importantly, these imaging and cognitive longitudinal effects were still observed after excluding the 15 patients who had been hospitalised. These mainly limbic brain imaging results may be the in vivo hallmarks of a degenerative spread of the disease through olfactory pathways, of neuroinflammatory events, or of the loss of sensory input due to anosmia. Whether this deleterious effect can be partially reversed, or whether these effects will persist in the long term, remains to be investigated with additional follow-up.

Universal dynamic fitting of magnetic resonance spectroscopy.

PURPOSE: Dynamic (2D) MRS is a collection of techniques where acquisitions of spectra are repeated under varying experimental or physiological conditions. Dynamic MRS comprises a rich set of contrasts, including diffusion-weighted, relaxation-weighted, functional, edited, or hyperpolarized spectroscopy, leading to quantitative insights into multiple physiological or microstructural processes. Conventional approaches to dynamic MRS analysis ignore the shared information between spectra, and instead proceed by independently fitting noisy individual spectra before modeling temporal changes in the parameters. Here, we propose a universal dynamic MRS toolbox which allows simultaneous fitting of dynamic spectra of arbitrary type. METHODS: A simple user-interface allows information to be shared and precisely modeled across spectra to make inferences on both spectral and dynamic processes. We demonstrate and thoroughly evaluate our approach in three types of dynamic MRS techniques. Simulations of functional and edited MRS are used to demonstrate the advantages of dynamic fitting. RESULTS: Analysis of synthetic functional 1H-MRS data shows a marked decrease in parameter uncertainty as predicted by prior work. Analysis with our tool replicates the results of two previously published studies using the original in vivo functional and diffusion-weighted data. Finally, joint spectral fitting with diffusion orientation models is demonstrated in synthetic data. CONCLUSION: A toolbox for generalized and universal fitting of dynamic, interrelated MR spectra has been released and validated. The toolbox is shared as a fully open-source software with comprehensive documentation, example data, and tutorials.

Diffusion-weighted MR spectroscopy: Consensus, recommendations, and resources from acquisition to modeling.

Brain cell structure and function reflect neurodevelopment, plasticity, and aging; and changes can help flag pathological processes such as neurodegeneration and neuroinflammation. Accurate and quantitative methods to noninvasively disentangle cellular structural features are needed and are a substantial focus of brain research. Diffusion-weighted MRS (dMRS) gives access to diffusion properties of endogenous intracellular brain metabolites that are preferentially located inside specific brain cell populations. Despite its great potential, dMRS remains a challenging technique on all levels: from the data acquisition to the analysis, quantification, modeling, and interpretation of results. These challenges were the motivation behind the organization of the Lorentz Center workshop on "Best Practices & Tools for Diffusion MR Spectroscopy" held in Leiden, the Netherlands, in September 2021. During the workshop, the dMRS community established a set of recommendations to execute robust dMRS studies. This paper provides a description of the steps needed for acquiring, processing, fitting, and modeling dMRS data, and provides links to useful resources.

GABA and Glutamate in hMT+ Link to Individual Differences in Residual Visual Function After Occipital Stroke.

BACKGROUND: Damage to the primary visual cortex following an occipital stroke causes loss of conscious vision in the contralateral hemifield. Yet, some patients retain the ability to detect moving visual stimuli within their blind field. The present study asked whether such individual differences in blind field perception following loss of primary visual cortex could be explained by the concentration of neurotransmitters γ-aminobutyric acid (GABA) and glutamate or activity of the visual motion processing, human middle temporal complex (hMT+). METHODS: We used magnetic resonance imaging in 19 patients with chronic occipital stroke to measure the concentration of neurotransmitters GABA and glutamate (proton magnetic resonance spectroscopy) and functional activity in hMT+ (functional magnetic resonance imaging). We also tested each participant on a 2-interval forced choice detection task using high-contrast, moving Gabor patches. We then measured and assessed the strength of relationships between participants' residual vision in their blind field and in vivo neurotransmitter concentrations, as well as visually evoked functional magnetic resonance imaging activity in their hMT+. Levels of GABA and glutamate were also measured in a sensorimotor region, which served as a control. RESULTS: Magnetic resonance spectroscopy-derived GABA and glutamate concentrations in hMT+ (but not sensorimotor cortex) strongly predicted blind-field visual detection abilities. Performance was inversely related to levels of both inhibitory and excitatory neurotransmitters in hMT+ but, surprisingly, did not correlate with visually evoked blood oxygenation level-dependent signal change in this motion-sensitive region. CONCLUSIONS: Levels of GABA and glutamate in hMT+ appear to provide superior information about motion detection capabilities inside perimetrically defined blind fields compared to blood oxygenation level-dependent signal changes-in essence, serving as biomarkers for the quality of residual visual processing in the blind-field. Whether they also reflect a potential for successful rehabilitation of visual function remains to be determined.

Tensor image registration library: Deformable registration of stand-alone histology images to whole-brain post-mortem MRI data.

BACKGROUND: Accurate registration between microscopy and MRI data is necessary for validating imaging biomarkers against neuropathology, and to disentangle complex signal dependencies in microstructural MRI. Existing registration methods often rely on serial histological sampling or significant manual input, providing limited scope to work with a large number of stand-alone histology sections. Here we present a customisable pipeline to assist the registration of stand-alone histology sections to whole-brain MRI data. METHODS: Our pipeline registers stained histology sections to whole-brain post-mortem MRI in 4 stages, with the help of two photographic intermediaries: a block face image (to undistort histology sections) and coronal brain slab photographs (to insert them into MRI space). Each registration stage is implemented as a configurable stand-alone Python script using our novel platform, Tensor Image Registration Library (TIRL), which provides flexibility for wider adaptation. We report our experience of registering 87 PLP-stained histology sections from 14 subjects and perform various experiments to assess the accuracy and robustness of each stage of the pipeline. RESULTS: All 87 histology sections were successfully registered to MRI. Histology-to-block registration (Stage 1) achieved 0.2-0.4 mm accuracy, better than commonly used existing methods. Block-to-slice matching (Stage 2) showed great robustness in automatically identifying and inserting small tissue blocks into whole brain slices with 0.2 mm accuracy. Simulations demonstrated sub-voxel level accuracy (0.13 mm) of the slice-to-volume registration (Stage 3) algorithm, which was observed in over 200 actual brain slice registrations, compensating 3D slice deformations up to 6.5 mm. Stage 4 combined the previous stages and generated refined pixelwise aligned multi-modal histology-MRI stacks. CONCLUSIONS: Our open-source pipeline provides robust automation tools for registering stand-alone histology sections to MRI data with sub-voxel level precision, and the underlying framework makes it readily adaptable to a diverse range of microscopy-MRI studies.

Functional connectivity between interoceptive brain regions is associated with distinct health-related domains: A population-based neuroimaging study.

Interoception is the sensation, perception, and integration of signals from within the body. It has been associated with a broad range of physiological and psychological processes. Further, interoceptive variables are related to specific regions and networks in the human brain. However, it is not clear whether or how these networks relate empirically to different domains of physiological and psychological health at the population level. We analysed a data set of 19,020 individuals (10,055 females, 8965 males; mean age: 63 years, age range: 45-81 years), who have participated in the UK Biobank Study, a very large-scale prospective epidemiological health study. Using canonical correlation analysis (CCA), allowing for the examination of associations between two sets of variables, we related the functional connectome of brain regions implicated in interoception to a selection of nonimaging health and lifestyle related phenotypes, exploring their relationship within modes of population co-variation. In one integrated and data driven analysis, we obtained four statistically significant modes. Modes could be categorised into domains of arousal and affect and cardiovascular health, respiratory health, body mass, and subjective health (all p < .0001) and were meaningfully associated with distinct neural circuits. Circuits represent specific neural "fingerprints" of functional domains and set the scope for future studies on the neurobiology of interoceptive involvement in different lifestyle and health-related phenotypes. Therefore, our research contributes to the conceptualisation of interoception and may lead to a better understanding of co-morbid conditions in the light of shared interoceptive structures.

Longitudinal connections and the organization of the temporal cortex in macaques, great apes, and humans.

The temporal association cortex is considered a primate specialization and is involved in complex behaviors, with some, such as language, particularly characteristic of humans. The emergence of these behaviors has been linked to major differences in temporal lobe white matter in humans compared with monkeys. It is unknown, however, how the organization of the temporal lobe differs across several anthropoid primates. Therefore, we systematically compared the organization of the major temporal lobe white matter tracts in the human, gorilla, and chimpanzee great apes and in the macaque monkey. We show that humans and great apes, in particular the chimpanzee, exhibit an expanded and more complex occipital-temporal white matter system; additionally, in humans, the invasion of dorsal tracts into the temporal lobe provides a further specialization. We demonstrate the reorganization of different tracts along the primate evolutionary tree, including distinctive connectivity of human temporal gray matter.

Human decisions about when to act originate within a basal forebrain-nigral circuit.

Decisions about when to act are critical for survival in humans as in animals, but how a desire is translated into the decision that an action is worth taking at any particular point in time is incompletely understood. Here we show that a simple model developed to explain when animals decide it is worth taking an action also explains a significant portion of the variance in timing observed when humans take voluntary actions. The model focuses on the current environment's potential for reward, the timing of the individual's own recent actions, and the outcomes of those actions. We show, by using ultrahigh-field MRI scanning, that in addition to anterior cingulate cortex within medial frontal cortex, a group of subcortical structures including striatum, substantia nigra, basal forebrain (BF), pedunculopontine nucleus (PPN), and habenula (HB) encode trial-by-trial variation in action time. Further analysis of the activity patterns found in each area together with psychophysiological interaction analysis and structural equation modeling suggested a model in which BF integrates contextual information that will influence the decision about when to act and communicates this information, in parallel with PPN and HB influences, to nigrostriatal circuits. It is then in the nigrostriatal circuit that action initiation per se begins.

Identifying microstructural changes in diffusion MRI; How to circumvent parameter degeneracy.

Biophysical models that attempt to infer real-world quantities from data usually have many free parameters. This over-parameterisation can result in degeneracies in model inversion and render parameter estimation ill-posed. However, in many applications, we are not interested in quantifying the parameters per se, but rather in identifying changes in parameters between experimental conditions (e.g. patients vs controls). Here we present a Bayesian framework to make inference on changes in the parameters of biophysical models even when model inversion is degenerate, which we refer to as Bayesian EstimatioN of CHange (BENCH). We infer the parameter changes in two steps; First, we train models that can estimate the pattern of change in the measurements given any hypothetical direction of change in the parameters using simulations. Next, for any pair of real data sets, we use these pre-trained models to estimate the probability that an observed difference in the data can be explained by each model of change. BENCH is applicable to any type of data and models and particularly useful for biophysical models with parameter degeneracies, where we can assume the change is sparse. In this paper, we apply the approach in the context of microstructural modelling of diffusion MRI data, where the models are usually over-parameterised and not invertible without injecting strong assumptions. Using simulations, we show that in the context of the standard model of white matter our approach is able to identify changes in microstructural parameters from conventional multi-shell diffusion MRI data. We also apply our approach to a subset of subjects from the UK-Biobank Imaging to identify the dominant standard model parameter change in areas of white matter hyperintensities under the assumption that the standard model holds in white matter hyperintensities.

Accurate predictions of individual differences in task-evoked brain activity from resting-state fMRI using a sparse ensemble learner.

Modelling and predicting individual differences in task-fMRI activity can have a wide range of applications from basic to clinical neuroscience. It has been shown that models based on resting-state activity can have high predictive accuracy. Here we propose several improvements to such models. Using a sparse ensemble learner, we show that (i) features extracted using Stochastic Probabilistic Functional Modes (sPROFUMO) outperform the previously proposed dual-regression approach, (ii) that the shape and overall intensity of individualised task activations can be modelled separately and explicitly, (iii) training the model on predicting residual differences in brain activity further boosts individualised predictions. These results hold for both surface-based analyses of the Human Connectome Project data as well as volumetric analyses of UK-biobank data. Overall, our model achieves state of the art prediction accuracy on par with the test-retest reliability of task-fMRI scans, suggesting that it has potential to supplement traditional task localisers.

An automated pipeline for extracting histological stain area fraction for voxelwise quantitative MRI-histology comparisons.

The acquisition of MRI and histology in the same post-mortem tissue sample enables direct correlation between MRI and histologically-derived parameters. However, there still lacks a standardised automated pipeline to process histology data, with most studies relying on manual intervention. Here, we introduce an automated pipeline to extract a quantitative histological measure for staining density (stain area fraction, SAF) from multiple immunohistochemical (IHC) stains. The pipeline is designed to directly address key IHC artefacts related to tissue staining and slide digitisation. Here, the pipeline was applied to post-mortem human brain data from multiple subjects, relating MRI parameters (FA, MD, RD, AD, R2*, R1) to IHC slides stained for myelin, neurofilaments, microglia and activated microglia. Utilising high-quality MRI-histology co-registrations, we then performed whole-slide voxelwise comparisons (simple correlations, partial correlations and multiple regression analyses) between multimodal MRI- and IHC-derived parameters. The pipeline was found to be reproducible, robust to artefacts and generalisable across multiple IHC stains. Our partial correlation results suggest that some simple MRI-SAF correlations should be interpreted with caution, due to the co-localisation of other tissue features (e.g., myelin and neurofilaments). Further, we find activated microglia-a generic biomarker of inflammation-to consistently be the strongest predictor of high DTI FA and low RD, which may suggest sensitivity of diffusion MRI to aspects of neuroinflammation related to microglial activation, even after accounting for other microstructural changes (demyelination, axonal loss and general microglia infiltration). Together, these results show the utility of this approach in carefully curating IHC data and performing multimodal analyses to better understand microstructural relationships with MRI.

Estimating axial diffusivity in the NODDI model.

To estimate microstructure-related parameters from diffusion MRI data, biophysical models make strong, simplifying assumptions about the underlying tissue. The extent to which many of these assumptions are valid remains an open research question. This study was inspired by the disparity between the estimated intra-axonal axial diffusivity from literature and that typically assumed by the Neurite Orientation Dispersion and Density Imaging (NODDI) model (d∥=1.7µm2/ms). We first demonstrate how changing the assumed axial diffusivity results in considerably different NODDI parameter estimates. Second, we illustrate the ability to estimate axial diffusivity as a free parameter of the model using high b-value data and an adapted NODDI framework. Using both simulated and in vivo data we investigate the impact of fitting to either real-valued or magnitude data, with Gaussian and Rician noise characteristics respectively, and what happens if we get the noise assumptions wrong in this high b-value and thus low SNR regime. Our results from real-valued human data estimate intra-axonal axial diffusivities of ∼2-2.5µm2/ms, in line with current literature. Crucially, our results demonstrate the importance of accounting for both a rectified noise floor and/or a signal offset to avoid biased parameter estimates when dealing with low SNR data.

Modelling white matter in gyral blades as a continuous vector field.

Many brain imaging studies aim to measure structural connectivity with diffusion tractography. However, biases in tractography data, particularly near the boundary between white matter and cortical grey matter can limit the accuracy of such studies. When seeding from the white matter, streamlines tend to travel parallel to the convoluted cortical surface, largely avoiding sulcal fundi and terminating preferentially on gyral crowns. When seeding from the cortical grey matter, streamlines generally run near the cortical surface until reaching deep white matter. These so-called "gyral biases" limit the accuracy and effective resolution of cortical structural connectivity profiles estimated by tractography algorithms, and they do not reflect the expected distributions of axonal densities seen in invasive tracer studies or stains of myelinated fibres. We propose an algorithm that concurrently models fibre density and orientation using a divergence-free vector field within gyral blades to encourage an anatomically-justified streamline density distribution along the cortical white/grey-matter boundary while maintaining alignment with the diffusion MRI estimated fibre orientations. Using in vivo data from the Human Connectome Project, we show that this algorithm reduces tractography biases. We compare the structural connectomes to functional connectomes from resting-state fMRI, showing that our model improves cross-modal agreement. Finally, we find that after parcellation the changes in the structural connectome are very minor with slightly improved interhemispheric connections (i.e, more homotopic connectivity) and slightly worse intrahemispheric connections when compared to tracers.

Hierarchical modelling of functional brain networks in population and individuals from big fMRI data.

A major goal of large-scale brain imaging datasets is to provide resources for investigating heterogeneous populations. Characterisation of functional brain networks for individual subjects from these datasets will have an enormous potential for prediction of cognitive or clinical traits. We propose for the first time a technique, Stochastic Probabilistic Functional Modes (sPROFUMO), that is scalable to UK Biobank (UKB) with expected 100,000 participants, and hierarchically estimates functional brain networks in individuals and the population, while allowing for bidirectional flow of information between the two. Using simulations, we show the model's utility, especially in scenarios that involve significant cross-subject variability, or require delineation of fine-grained differences between the networks. Subsequently, by applying the model to resting-state fMRI from 4999 UKB subjects, we mapped resting state networks (RSNs) in single subjects with greater detail than has been possible previously in UKB (>100 RSNs), and demonstrate that these RSNs can predict a range of sensorimotor and higher-level cognitive functions. Furthermore, we demonstrate several advantages of the model over independent component analysis combined with dual-regression (ICA-DR), particularly with respect to the estimation of the spatial configuration of the RSNs and the predictive power for cognitive traits. The proposed model and results can open a new door for future investigations into individualised profiles of brain function from big data.

Clinical applications of magnetic resonance imaging based functional and structural connectivity.

Advances in computational neuroimaging techniques have expanded the armamentarium of imaging tools available for clinical applications in clinical neuroscience. Non-invasive, in vivo brain MRI structural and functional network mapping has been used to identify therapeutic targets, define eloquent brain regions to preserve, and gain insight into pathological processes and treatments as well as prognostic biomarkers. These tools have the real potential to inform patient-specific treatment strategies. Nevertheless, a realistic appraisal of clinical utility is needed that balances the growing excitement and interest in the field with important limitations associated with these techniques. Quality of the raw data, minutiae of the processing methodology, and the statistical models applied can all impact on the results and their interpretation. A lack of standardization in data acquisition and processing has also resulted in issues with reproducibility. This limitation has had a direct impact on the reliability of these tools and ultimately, confidence in their clinical use. Advances in MRI technology and computational power as well as automation and standardization of processing methods, including machine learning approaches, may help address some of these issues and make these tools more reliable in clinical use. In this review, we will highlight the current clinical uses of MRI connectomics in the diagnosis and treatment of neurological disorders; balancing emerging applications and technologies with limitations of connectivity analytic approaches to present an encompassing and appropriate perspective.

FSL-MRS: An end-to-end spectroscopy analysis package.

PURPOSE: We introduce FSL-MRS, an end-to-end, modular, open-source MRS analysis toolbox. It provides spectroscopic data conversion, preprocessing, spectral simulation, fitting, quantitation, and visualization. METHODS: The FSL-MRS package is modular. Its programs operate on data in a standard format (Neuroimaging Informatics Technology Initiative [NIfTI]) capable of storing single-voxel and multivoxel spectroscopy, including spatial orientation information. The FSL-MRS toolbox includes tools for preprocessing of raw spectroscopy data, including coil combination, frequency and phase alignment, and filtering. A density matrix simulation program is supplied for generation of basis spectra from simple text-based descriptions of pulse sequences. Fitting is based on linear combination of basis spectra and implements Markov chain Monte Carlo optimization for the estimation of the full posterior distribution of metabolite concentrations. Validation of the fitting is carried out on independently created simulated data, phantom data, and three in vivo human data sets (257 single-voxel spectroscopy and 8 MRSI data sets) at 3 T and 7 T. Interactive HTML reports are automatically generated by processing and fitting stages of the toolbox. The FSL-MRS package can be used on the command line or interactively in the Python language. RESULTS: Validation of the fitting shows low error in simulation (median error of 11.9%) and in phantom (3.4%). Average correlation between a third-party toolbox (LCModel) and FSL-MRS was high (0.53-0.81) in all three in vivo data sets. CONCLUSION: The FSL-MRS toolbox is designed to be flexible and extensible to new forms of spectroscopic acquisitions. Custom fitting models can be specified within the framework for dynamic or multivoxel spectroscopy. It is available as part of the FMRIB Software Library.

Quantifying myelin in crossing fibers using diffusion-prepared phase imaging: Theory and simulations.

PURPOSE: Myelin has long been the target of neuroimaging research. However, most available techniques can only provide a voxel-averaged estimate of myelin content. In the human brain, white matter fiber pathways connecting different brain areas and carrying different functions often cross each other in the same voxel. A measure that can differentiate the degree of myelination of crossing fibers would provide a more specific marker of myelination. THEORY AND METHODS: One MRI signal property that is sensitive to myelin is the phase accumulation. This sensitivity is used by measuring the phase accumulation of the signal remaining after diffusion-weighting, which is called diffusion-prepared phase imaging (DIPPI). Including diffusion-weighting before estimating the phase accumulation has two distinct advantages for estimating the degree of myelination: (1) It increases the relative contribution of intra-axonal water, whose phase is related linearly to the thickness of the surrounding myelin (in particular the log g-ratio); and (2) it gives directional information, which can be used to distinguish between crossing fibers. Here the DIPPI sequence is described, an approach is proposed to estimate the log g-ratio, and simulations are used and DIPPI data acquired in an isotropic phantom to quantify other sources of phase accumulation. RESULTS: The expected bias is estimated in the log g-ratio for reasonable in vivo acquisition parameters caused by eddy currents (~4%-10%), remaining extra-axonal signal (~15%), and gradients in the bulk off-resonance field (<10% for most of the brain). CONCLUSION: This new sequence may provide a g-ratio estimate per fiber population crossing within a voxel.

Structural Variability in the Human Brain Reflects Fine-Grained Functional Architecture at the Population Level.

Human brain structure topography is thought to be related in part to functional specialization. However, the extent of such relationships is unclear. Here, using a data-driven, multimodal approach for studying brain structure across the lifespan (N = 484, n = 260 females), we demonstrate that numerous structural networks, covering the entire brain, follow a functionally meaningful architecture. These gray matter networks (GMNs) emerge from the covariation of gray matter volume and cortical area at the population level. We further reveal fine-grained anatomical signatures of functional connectivity. For example, within the cerebellum, a structural separation emerges between lobules that are functionally connected to distinct, mainly sensorimotor, cognitive and limbic regions of the cerebral cortex and subcortex. Structural modes of variation also replicate the fine-grained functional architecture seen in eight well defined visual areas in both task and resting-state fMRI. Furthermore, our study shows a structural distinction corresponding to the established segregation between anterior and posterior default-mode networks (DMNs). These fine-grained GMNs further cluster together to form functionally meaningful larger-scale organization. In particular, we identify a structural architecture bringing together the functional posterior DMN and its anticorrelated counterpart. In summary, our results demonstrate that the relationship between structural and functional connectivity is fine-grained, widespread across the entire brain, and driven by covariation in cortical area, i.e. likely differences in shape, depth, or number of foldings. These results suggest that neurotrophic events occur during development to dictate that the size and folding pattern of distant, functionally connected brain regions should vary together across subjects.SIGNIFICANCE STATEMENT Questions about the relationship between structure and function in the human brain have engaged neuroscientists for centuries in a debate that continues to this day. Here, by investigating intersubject variation in brain structure across a large number of individuals, we reveal modes of structural variation that map onto fine-grained functional organization across the entire brain, and specifically in the cerebellum, visual areas, and default-mode network. This functionally meaningful structural architecture emerges from the covariation of gray matter volume and cortical folding. These results suggest that the neurotrophic events at play during development, and possibly evolution, which dictate that the size and folding pattern of distant brain regions should vary together across subjects, might also play a role in functional cortical specialization.

Improved fibre dispersion estimation using b-tensor encoding.

Measuring fibre dispersion in white matter with diffusion magnetic resonance imaging (MRI) is limited by an inherent degeneracy between fibre dispersion and microscopic diffusion anisotropy (i.e., the diffusion anisotropy expected for a single fibre orientation). This means that estimates of fibre dispersion rely on strong assumptions, such as constant microscopic anisotropy throughout the white matter or specific biophysical models. Here we present a simple approach for resolving this degeneracy using measurements that combine linear (conventional) and spherical tensor diffusion encoding. To test the accuracy of the fibre dispersion when our microstructural model is only an approximation of the true tissue structure, we simulate multi-compartment data and fit this with a single-compartment model. For such overly simplistic tissue assumptions, we show that the bias in fibre dispersion is greatly reduced (~5x) for single-shell linear and spherical tensor encoding data compared with single-shell or multi-shell conventional data. In in-vivo data we find a consistent estimate of fibre dispersion as we reduce the b-value from 3 to 1.5 ms/µm2, increase the repetition time, increase the echo time, or increase the diffusion time. We conclude that the addition of spherical tensor encoded data to conventional linear tensor encoding data greatly reduces the sensitivity of the estimated fibre dispersion to the model assumptions of the tissue microstructure.