Matteo Bastiani

Human cortical connectome reconstruction from diffusion weighted MRI: The effect of tractography algorithm

Reconstructing the macroscopic human cortical connectome by Diffusion Weighted Imaging (DWI) is a challenging research topic that has recently gained a lot of attention. In the present work, we investigate the effects of intra-voxel fiber direction modeling and tractography algorithm on derived structural network indices (e.g. density, small-worldness and global efficiency). The investigation is centered on three semi-independent distinctions within the large set of available diffusion models and tractography methods: i) single fiber direction versus multiple directions in the intra-voxel diffusion model, ii) deterministic versus probabilistic tractography and iii) local versus global measure-of-fit of the reconstructed fiber trajectories. The effect of algorithm and parameter choice has two components. First, there is the large effect of tractography algorithm and parameters on global network density, which is known to strongly affect graph indices. Second, and more importantly, there are remaining effects on graph indices which range in the tens of percent even when global density is controlled for. This is crucial for the sensitivity of any human structural network study and for the validity of study comparisons. We then investigate the effect of the choice of tractography algorithm on sensitivity and specificity of the resulting connections with a connectome dissection quality control (QC) approach. In this approach, evaluation of Tract Specific Density Coefficients (TSDCs) measures sensitivity while careful inspection of tractography path results assesses specificity. We use this to discuss interactions in the combined effects of these methods and implications for future studies.

Histological validation of high-resolution DTI in human post mortem tissue

Diffusion tensor imaging (DTI) is amongst the simplest mathematical models available for diffusion magnetic resonance imaging, yet still by far the most used one. Despite the success of DTI as an imaging tool for white matter fibers, its anatomical underpinnings on a microstructural basis remain unclear. In this study, we used 65 myelin-stained sections of human premotor cortex to validate modeled fiber orientations and oft used microstructure-sensitive scalar measures of DTI on the level of individual voxels. We performed this validation on high spatial resolution diffusion MRI acquisitions investigating both white and gray matter. We found a very good agreement between DTI and myelin orientations with the majority of voxels showing angular differences less than 10°. The agreement was strongest in white matter, particularly in unidirectional fiber pathways. In gray matter, the agreement was good in the deeper layers highlighting radial fiber directions even at lower fractional anisotropy (FA) compared to white matter. This result has potentially important implications for tractography algorithms applied to high resolution diffusion MRI data if the aim is to move across the gray/white matter boundary. We found strong relationships between myelin microstructure and DTI-based microstructure-sensitive measures. High FA values were linked to high myelin density and a sharply tuned histological orientation profile. Conversely, high values of mean diffusivity (MD) were linked to bimodal or diffuse orientation distributions and low myelin density. At high spatial resolution, DTI-based measures can be highly sensitive to white and gray matter microstructure despite being relatively unspecific to concrete microarchitectural aspects.

Comparative analysis of the macroscale structural connectivity in the macaque and human brain.

The macaque brain serves as a model for the human brain, but its suitability is challenged by unique human features, including connectivity reconfigurations, which emerged during primate evolution. We perform a quantitative comparative analysis of the whole brain macroscale structural connectivity of the two species. Our findings suggest that the human and macaque brain as a whole are similarly wired. A region-wise analysis reveals many interspecies similarities of connectivity patterns, but also lack thereof, primarily involving cingulate regions. We unravel a common structural backbone in both species involving a highly overlapping set of regions. This structural backbone, important for mediating information across the brain, seems to constitute a feature of the primate brain persevering evolution. Our findings illustrate novel evolutionary aspects at the macroscale connectivity level and offer a quantitative translational bridge between macaque and human research.

The developing human connectome project: A minimal processing pipeline for neonatal cortical surface reconstruction

The Developing Human Connectome Project (dHCP) seeks to create the first 4-dimensional connectome of early life. Understanding this connectome in detail may provide insights into normal as well as abnormal patterns of brain development. Following established best practices adopted by the WU-MINN Human Connectome Project (HCP), and pioneered by FreeSurfer, the project utilises cortical surface-based processing pipelines. In this paper, we propose a fully automated processing pipeline for the structural Magnetic Resonance Imaging (MRI) of the developing neonatal brain. This proposed pipeline consists of a refined framework for cortical and sub-cortical volume segmentation, cortical surface extraction, and cortical surface inflation, which has been specifically designed to address considerable differences between adult and neonatal brains, as imaged using MRI. Using the proposed pipeline our results demonstrate that images collected from 465 subjects ranging from 28 to 45 weeks post-menstrual age (PMA) can be processed fully automatically; generating cortical surface models that are topologically correct, and correspond well with manual evaluations of tissue boundaries in 85% of cases. Results improve on state-of-the-art neonatal tissue segmentation models and significant errors were found in only 2% of cases, where these corresponded to subjects with high motion. Downstream, these surfaces will enhance comparisons of functional and diffusion MRI datasets, supporting the modelling of emerging patterns of brain connectivity.

Towards a comprehensive framework for movement and distortion correction of diffusion MR images: Within volume movement

ABSTRACT Most motion correction methods work by aligning a set of volumes together, or to a volume that represents a reference location. These are based on an implicit assumption that the subject remains motionless during the several seconds it takes to acquire all slices in a volume, and that any movement occurs in the brief moment between acquiring the last slice of one volume and the first slice of the next. This is clearly an approximation that can be more or less good depending on how long it takes to acquire one volume and in how rapidly the subject moves. In this paper we present a method that increases the temporal resolution of the motion correction by modelling movement as a piecewise continous function over time. This intra‐volume movement correction is implemented within a previously presented framework that simultaneously estimates distortions, movement and movement‐induced signal dropout. We validate the method on highly realistic simulated data containing all of these effects. It is demonstrated that we can estimate the true movement with high accuracy, and that scalar parameters derived from the data, such as fractional anisotropy, are estimated with greater fidelity when data has been corrected for intra‐volume movement. Importantly, we also show that the difference in fidelity between data affected by different amounts of movement is much reduced when taking intra‐volume movement into account. Additional validation was performed on data from a healthy volunteer scanned when lying still and when performing deliberate movements. We show an increased correspondence between the “still” and the “movement” data when the latter is corrected for intra‐volume movement. Finally we demonstrate a big reduction in the telltale signs of intra‐volume movement in data acquired on elderly subjects. HighlightsWe introduce a method to correct for intra‐volume movement into an existing framework for movement and distortion correction.It has been validated on realistic simulated data and on experimental data with deliberate movement.The results indicate that one can reliably reverse the adverse effects of intra‐volume movement.

Unraveling the multiscale structural organization and connectivity of the human brain: the role of diffusion MRI

The structural architecture and the anatomical connectivity of the human brain show different organizational principles at distinct spatial scales. Histological staining and light microscopy techniques have been widely used in classical neuroanatomical studies to unravel brain organization. Using such techniques is a laborious task performed on 2-dimensional histological sections by skilled anatomists possibly aided by semi-automated algorithms. With the recent advent of modern magnetic resonance imaging (MRI) contrast mechanisms, cortical layers and columns can now be reliably identified and their structural properties quantified post-mortem. These developments are allowing the investigation of neuroanatomical features of the brain at a spatial resolution that could be interfaced with that of histology. Diffusion MRI and tractography techniques, in particular, have been used to probe the architecture of both white and gray matter in three dimensions. Combined with mathematical network analysis, these techniques are increasingly influential in the investigation of the macro-, meso-, and microscopic organization of brain connectivity and anatomy, both in vivo and ex vivo. Diffusion MRI-based techniques in combination with histology approaches can therefore support the endeavor of creating multimodal atlases that take into account the different spatial scales or levels on which the brain is organized. The aim of this review is to illustrate and discuss the structural architecture and the anatomical connectivity of the human brain at different spatial scales and how recently developed diffusion MRI techniques can help investigate these.

White matter microstructural changes in adolescent anorexia nervosa including an exploratory longitudinal study

Background Anorexia nervosa (AN) often begins in adolescence, however, the understanding of the underlying pathophysiology at this developmentally important age is scarce, impeding early interventions. We used diffusion tensor imaging (DTI) to investigate microstructural white matter (WM) brain changes including an experimental longitudinal follow-up. Methods We acquired whole brain diffusion-weighted brain scans of 22 adolescent female hospitalized patients with AN at admission and nine patients longitudinally at discharge after weight rehabilitation. Patients (10–18 years) were compared to 21 typically developing controls (TD). Tract-based spatial statistics (TBSS) were applied to compare fractional anisotropy (FA) across groups and time points. Associations between average FA values of the global WM skeleton and weight as well as illness duration parameters were analyzed by multiple linear regression. Results We observed increased FA in bilateral frontal, parietal and temporal areas in AN patients at admission compared to TD. Higher FA of the global WM skeleton at admission was associated with faster weight loss prior to admission. Exploratory longitudinal analysis showed this FA increase to be partially normalized after weight rehabilitation. Conclusions Our findings reveal a markedly different pattern of WM microstructural changes in adolescent AN compared to most previous results in adult AN. This could signify a different susceptibility and reaction to semi-starvation in the still developing brain of adolescents or a time-dependent pathomechanism differing with extend of chronicity. Higher FA at admission in adolescents with AN could point to WM fibers being packed together more closely.

Automatic Segmentation of Human Cortical Layer-Complexes and Architectural Areas Using Ex vivo Diffusion MRI and Its Validation

Recently, several magnetic resonance imaging contrast mechanisms have been shown to distinguish cortical substructure corresponding to selected cortical layers. Here, we investigate cortical layer and area differentiation by automatized unsupervised clustering of high-resolution diffusion MRI data. Several groups of adjacent layers could be distinguished in human primary motor and premotor cortex. We then used the signature of diffusion MRI signals along cortical depth as a criterion to detect area boundaries and find borders at which the signature changes abruptly. We validate our clustering results by histological analysis of the same tissue. These results confirm earlier studies which show that diffusion MRI can probe layer-specific intracortical fiber organization and, moreover, suggests that it contains enough information to automatically classify architecturally distinct cortical areas. We discuss the strengths and weaknesses of the automatic clustering approach and its appeal for MR-based cortical histology.

Optimal short-time acquisition schemes in high angular resolution diffusion-weighted imaging.

This work investigates the possibilities of applying high-angular-resolution-diffusion-imaging- (HARDI-) based methods in a clinical setting by investigating the performance of non-Gaussian diffusion probability density function (PDF) estimation for a range of b-values and diffusion gradient direction tables. It does so at realistic SNR levels achievable in limited time on a high-performance 3T system for the whole human brain in vivo. We use both computational simulations and in vivo brain scans to quantify the angular resolution of two selected reconstruction methods: Q-ball imaging and the diffusion orientation transform. We propose a new analytical solution to the ODF derived from the DOT. Both techniques are analytical decomposition approaches that require identical acquisition and modest postprocessing times and, given the proposed modifications of the DOT, can be analyzed in a similar fashion. We find that an optimal HARDI protocol given a stringent time constraint (<10 min) combines a moderate b-value (around 2000 s/mm2) with a relatively low number of acquired directions (>48). Our findings generalize to other methods and additional improvements in MR acquisition techniques.

Improved tractography using asymmetric fibre orientation distributions

Diffusion MRI allows us to make inferences on the structural organisation of the brain by mapping water diffusion to white matter microstructure. However, such a mapping is generally ill-defined; for instance, diffusion measurements are antipodally symmetric (diffusion along x and –x are equal), whereas the distribution of fibre orientations within a voxel is generally not symmetric. Therefore, different sub-voxel patterns such as crossing, fanning, or sharp bending, cannot be distinguished by fitting a voxel-wise model to the signal. However, asymmetric fibre patterns can potentially be distinguished once spatial information from neighbouring voxels is taken into account. We propose a neighbourhood-constrained spherical deconvolution approach that is capable of inferring asymmetric fibre orientation distributions (A-fods). Importantly, we further design and implement a tractography algorithm that utilises the estimated A-fods, since the commonly used streamline tractography paradigm cannot directly take advantage of the new information. We assess performance using ultra-high resolution histology data where we can compare true orientation distributions against sub-voxel fibre patterns estimated from down-sampled data. Finally, we explore the benefits of A-fods-based tractography using in vivo data by evaluating agreement of tractography predictions with connectivity estimates made using different in-vivo modalities. The proposed approach can reliably estimate complex fibre patterns such as sharp bending and fanning, which voxel-wise approaches cannot estimate. Moreover, histology-based and in-vivo results show that the new framework allows more accurate tractography and reconstruction of maps quantifying (symmetric and asymmetric) fibre complexity.