Non-negative data-driven mapping of structural connections with application to the neonatal brain.

Mapping connections in the neonatal brain can provide insight into the crucial early stages of neurodevelopment that shape brain organisation and lay the foundations for cognition and behaviour. Diffusion MRI and tractography provide unique opportunities for such explorations, through estimation of white matter bundles and brain connectivity. Atlas-based tractography protocols, i.e. a priori defined sets of masks and logical operations in a template space, have been commonly used in the adult brain to drive such explorations. However, rapid growth and maturation of the brain during early development make it challenging to ensure correspondence and validity of such atlas-based tractography approaches in the developing brain. An alternative can be provided by data-driven methods, which do not depend on predefined regions of interest. Here, we develop a novel data-driven framework to extract white matter bundles and their associated grey matter networks from neonatal tractography data, based on non-negative matrix factorisation that is inherently suited to the non-negative nature of structural connectivity data. We also develop a non-negative dual regression framework to map group-level components to individual subjects. Using in-silico simulations, we evaluate the accuracy of our approach in extracting connectivity components and compare with an alternative data-driven method, independent component analysis. We apply non-negative matrix factorisation to whole-brain connectivity obtained from publicly available datasets from the Developing Human Connectome Project, yielding grey matter components and their corresponding white matter bundles. We assess the validity and interpretability of these components against traditional tractography results and grey matter networks obtained from resting-state fMRI in the same subjects. We subsequently use them to generate a parcellation of the neonatal cortex using data from 323 new-born babies and we assess the robustness and reproducibility of this connectivity-driven parcellation.

Bayesian Optimisation of Large-Scale Biophysical Networks.

The relationship between structure and function in the human brain is well established, but not yet well characterised. Large-scale biophysical models allow us to investigate this relationship, by leveraging structural information (e.g. derived from diffusion tractography) in order to couple dynamical models of local neuronal activity into networks of interacting regions distributed across the cortex. In practice however, these models are difficult to parametrise, and their simulation is often delicate and computationally expensive. This undermines the experimental aspect of scientific modelling, and stands in the way of comparing different parametrisations, network architectures, or models in general, with confidence. Here, we advocate the use of Bayesian optimisation for assessing the capabilities of biophysical network models, given a set of desired properties (e.g. band-specific functional connectivity); and in turn the use of this assessment as a principled basis for incremental modelling and model comparison. We adapt an optimisation method designed to cope with costly, high-dimensional, non-convex problems, and demonstrate its use and effectiveness. Using five parameters controlling key aspects of our model, we find that this method is able to converge to regions of high functional similarity with real MEG data, with very few samples given the number of parameters, without getting stuck in local extrema, and while building and exploiting a map of uncertainty defined smoothly across the parameter space. We compare the results obtained using different methods of structural connectivity estimation from diffusion tractography, and find that one method leads to better simulations.

On the mechanical behaviour of PEEK and HA cranial implants under impact loading.

The human head can be subjected to numerous impact loadings such as those produced by a fall or during sport activities. These accidents can result in skull fracture and in some complex cases, part of the skull may need to be replaced by a biomedical implant. Even when the skull is not damaged, such accidents can result in brain swelling treated by decompressive craniectomy. Usually, after recovery, the part of the skull that has been removed is replaced by a prosthesis. In such situations, a computational tool able to analyse the choice of prosthesis material depending on the patient's specific activity has the potential to be extremely useful for clinicians. The work proposed here focusses on the development and use of a numerical model for the analysis of cranial implants under impact conditions. In particular, two main biomaterials commonly employed for this kind of prosthesis are polyether-ether-ketone (PEEK) and macroporous hydroxyapatite (HA). In order to study the suitability of these implants, a finite element head model comprising scalp, skull, cerebral falx, cerebrospinal fluid and brain tissues, with a cranial implant replacing part of the skull has been developed from magnetic resonance imaging data. The human tissues and these two biocompatible materials have been independently studied and their constitutive models are provided here. A computational model of the human head under impact loading is then implemented and validated, and a numerical comparison of the mechanical impact response of PEEK and HA implants is presented. This comparison was carried out in terms of the effectiveness of both implants in ensuring structural integrity and preventing traumatic brain injury. The results obtained in this work highlight the need to take into account environmental mechanical considerations to select the optimal implant depending on the specific patient: whereas HA implants present attractive biointegration properties, PEEK implant can potentially be a much more appropriate choice in a demanding mechanical life style. Finally, a novel methodology is proposed to assess the need for further clinical evaluation in case of impact with both implants over a large range of impact conditions.

A probabilistic atlas of the cerebellar white matter.

Imaging of the cerebellar cortex, deep cerebellar nuclei and their connectivity are gaining attraction, due to the important role the cerebellum plays in cognition and motor control. Atlases of the cerebellar cortex and nuclei are used to locate regions of interest in clinical and neuroscience studies. However, the white matter that connects these relay stations is of at least similar functional importance. Damage to these cerebellar white matter tracts may lead to serious language, cognitive and emotional disturbances, although the pathophysiological mechanism behind it is still debated. Differences in white matter integrity between patients and controls might shed light on structure-function correlations. A probabilistic parcellation atlas of the cerebellar white matter would help these studies by facilitating automatic segmentation of the cerebellar peduncles, the localization of lesions and the comparison of white matter integrity between patients and controls. In this work a digital three-dimensional probabilistic atlas of the cerebellar white matter is presented, based on high quality 3T, 1.25mm resolution diffusion MRI data from 90 subjects participating in the Human Connectome Project. The white matter tracts were estimated using probabilistic tractography. Results over 90 subjects were symmetrical and trajectories of superior, middle and inferior cerebellar peduncles resembled the anatomy as known from anatomical studies. This atlas will contribute to a better understanding of cerebellar white matter architecture. It may eventually aid in defining structure-function correlations in patients with cerebellar disorders.

High resolution whole brain diffusion imaging at 7T for the Human Connectome Project.

Mapping structural connectivity in healthy adults for the Human Connectome Project (HCP) benefits from high quality, high resolution, multiband (MB)-accelerated whole brain diffusion MRI (dMRI). Acquiring such data at ultrahigh fields (7T and above) can improve intrinsic signal-to-noise ratio (SNR), but suffers from shorter T2 and T2(⁎) relaxation times, increased B1(+) inhomogeneity (resulting in signal loss in cerebellar and temporal lobe regions), and increased power deposition (i.e. specific absorption rate (SAR)), thereby limiting our ability to reduce the repetition time (TR). Here, we present recent developments and optimizations in 7T image acquisitions for the HCP that allow us to efficiently obtain high quality, high resolution whole brain in-vivo dMRI data at 7T. These data show spatial details typically seen only in ex-vivo studies and complement already very high quality 3T HCP data in the same subjects. The advances are the result of intensive pilot studies aimed at mitigating the limitations of dMRI at 7T. The data quality and methods described here are representative of the datasets that will be made freely available to the community in 2015.

Effects of image reconstruction on fiber orientation mapping from multichannel diffusion MRI: reducing the noise floor using SENSE.

PURPOSE: To examine the effects of the reconstruction algorithm of magnitude images from multichannel diffusion MRI on fiber orientation estimation. THEORY AND METHODS: It is well established that the method used to combine signals from different coil elements in multichannel MRI can have an impact on the properties of the reconstructed magnitude image. Using a root-sum-of-squares approach results in a magnitude signal that follows an effective noncentral-χ distribution. As a result, the noise floor, the minimum measurable in the absence of any true signal, is elevated. This is particularly relevant for diffusion-weighted MRI, where the signal attenuation is of interest. RESULTS: In this study, we illustrate problems that such image reconstruction characteristics may cause in the estimation of fiber orientations, both for model-based and model-free approaches, when modern 32-channel coils are used. We further propose an alternative image reconstruction method that is based on sensitivity encoding (SENSE) and preserves the Rician nature of the single-channel, magnitude MR signal. We show that for the same k-space data, root-sum-of-squares can cause excessive overfitting and reduced precision in orientation estimation compared with the SENSE-based approach. CONCLUSION: These results highlight the importance of choosing the appropriate image reconstruction method for tractography studies that use multichannel receiver coils for diffusion MRI acquisition.