Accelerating fibre orientation estimation from diffusion weighted magnetic resonance imaging using GPUs

With the performance of central processing units (CPUs) having effectively reached a limit, parallel processing offers an alternative for applications with high computational demands. Modern graphics processing units (GPUs) are massively parallel processors that can execute simultaneously thousands of light-weight processes. In this study, we propose and implement a parallel GPU-based design of a popular method that is used for the analysis of brain magnetic resonance imaging (MRI). More specifically, we are concerned with a model-based approach for extracting tissue structural information from diffusion-weighted (DW) MRI data. DW-MRI offers, through tractography approaches, the only way to study brain structural connectivity, non-invasively and in-vivo. We parallelise the Bayesian inference framework for the ball & stick model, as it is implemented in the tractography toolbox of the popular FSL software package (University of Oxford). For our implementation, we utilise the Compute Unified Device Architecture (CUDA) programming model. We show that the parameter estimation, performed through Markov Chain Monte Carlo (MCMC), is accelerated by at least two orders of magnitude, when comparing a single GPU with the respective sequential single-core CPU version. We also illustrate similar speed-up factors (up to 120x) when comparing a multi-GPU with a multi-CPU implementation

Bayesian uncertainty quantification in linear models for diffusion MRI

Diffusion MRI (dMRI) is a valuable tool in the assessment of tissue microstructure. By fitting a model to the dMRI signal it is possible to derive various quantitative features. Several of the most popular dMRI signal models are expansions in an appropriately chosen basis, where the coefficients are determined using some variation of least-squares. However, such approaches lack any notion of uncertainty, which could be valuable in e.g. group analyses. In this work, we use a probabilistic interpretation of linear least-squares methods to recast popular dMRI models as Bayesian ones. This makes it possible to quantify the uncertainty of any derived quantity. In particular, for quantities that are affine functions of the coefficients, the posterior distribution can be expressed in closed-form. We simulated measurements from single- and double-tensor models where the correct values of several quantities are known, to validate that the theoretically derived quantiles agree with those observed empirically. We included results from residual bootstrap for comparison and found good agreement. The validation employed several different models: Diffusion Tensor Imaging (DTI), Mean Apparent Propagator MRI (MAP-MRI) and Constrained Spherical Deconvolution (CSD). We also used in vivo data to visualize maps of quantitative features and corresponding uncertainties, and to show how our approach can be used in a group analysis to downweight subjects with high uncertainty. In summary, we convert successful linear models for dMRI signal estimation to probabilistic models, capable of accurate uncertainty quantification.Comment: Added results from a group analysis and a comparison with residual bootstra

Advances in diffusion MRI acquisition and processing in the Human Connectome Project

The Human Connectome Project (HCP) is a collaborative 5-year effort to map human brain connections and their variability in healthy adults. A consortium of HCP investigators will study a population of 1200 healthy adults using multiple imaging modalities, along with extensive behavioral and genetic data. In this overview, we focus on diffusion MRI (dMRI) and the structural connectivity aspect of the project. We present recent advances in acquisition and processing that allow us to obtain very high-quality in-vivo MRI data, whilst enabling scanning of a very large number of subjects. These advances result from 2 years of intensive efforts in optimising many aspects of data acquisition and processing during the piloting phase of the project. The data quality and methods described here are representative of the datasets and processing pipelines that will be made freely available to the community at quarterly intervals, beginning in 2013

Improved Multi-Shot Diffusion-Weighted MRI with Zero-Shot Self-Supervised Learning Reconstruction

Diffusion MRI is commonly performed using echo-planar imaging (EPI) due to its rapid acquisition time. However, the resolution of diffusion-weighted images is often limited by magnetic field inhomogeneity-related artifacts and blurring induced by T2- and T2*-relaxation effects. To address these limitations, multi-shot EPI (msEPI) combined with parallel imaging techniques is frequently employed. Nevertheless, reconstructing msEPI can be challenging due to phase variation between multiple shots. In this study, we introduce a novel msEPI reconstruction approach called zero-MIRID (zero-shot self-supervised learning of Multi-shot Image Reconstruction for Improved Diffusion MRI). This method jointly reconstructs msEPI data by incorporating deep learning-based image regularization techniques. The network incorporates CNN denoisers in both k- and image-spaces, while leveraging virtual coils to enhance image reconstruction conditioning. By employing a self-supervised learning technique and dividing sampled data into three groups, the proposed approach achieves superior results compared to the state-of-the-art parallel imaging method, as demonstrated in an in-vivo experiment.Comment: 10 pages, 4 figure

A combined DTI-fMRI approach for optimizing the delineation of posteromedial vs. anterolateral entorhinal cortex

In the entorhinal cortex (EC), attempts have been made to identify the human homologue regions of the medial (MEC) and lateral (LEC) subdivision using either functional magnetic resonance imaging (fMRI) or diffusion tensor imaging (DTI). However, there are still discrepancies between entorhinal subdivisions depending on the choice of connectivity seed regions and the imaging modality used. While DTI can be used to follow the white matter tracts of the brain, fMRI can identify functionally connected brain regions. In this study, we used both DTI and resting-state fMRI in 103 healthy adults to investigate both structural and functional connectivity between the EC and associated cortical brain regions. Differential connectivity with these regions was then used to predict the locations of the human homologues of MEC and LEC. Our results from combining DTI and fMRI support a subdivision into posteromedial (pmEC) and anterolateral (alEC) EC and reveal a discrete border between the pmEC and alEC. Furthermore, the EC subregions obtained by either imaging modality showed similar distinct connectivity patterns: While pmEC showed increased connectivity preferentially with the default mode network, the alEC exhibited increased connectivity with regions in the dorsal attention and salience networks. Optimizing the delineation of the human homologues of MEC and LEC with a combined, cross-validated DTI-fMRI approach allows to define a likely border between the two subdivisions and has implications for both cognitive and translational neuroscience research

Ventralis intermedius nucleus anatomical variability assessment by MRI structural connectivity

The ventralis intermedius nucleus (Vim) is centrally placed in the dentato-thalamo-cortical pathway (DTCp) and is a key surgical target in the treatment of severe medically refractory tremor. It is not visible on conventional MRI sequences; consequently, stereotactic targeting currently relies on atlas-based coordinates. This fails to capture individual anatomical variability, which may lead to poor long-term clinical efficacy. Probabilistic tractography, combined with known anatomical connectivity, enables localisation of thalamic nuclei at an individual subject level. There are, however, a number of confounds associated with this technique that may influence results. Here we focused on an established method, using probabilistic tractography to reconstruct the DTCp, to identify the connectivity-defined Vim (cd-Vim) in vivo. Using 100 healthy individuals from the Human Connectome Project, our aim was to quantify cd-Vim variability across this population, measure the discrepancy with atlas-defined Vim (ad-Vim), and assess the influence of potential methodological confounds. We found no significant effect of any of the confounds. The mean cd-Vim coordinate was located within 1.88 mm (left) and 2.12 mm (right) of the average midpoint and 3.98 mm (left) and 5.41 mm (right) from the ad-Vim coordinates. cd-Vim location was more variable on the right, which reflects hemispheric asymmetries in the probabilistic DTC reconstructed. The method was reproducible, with no significant cd-Vim location differences in a separate test-retest cohort. The superior cerebellar peduncle was identified as a potential source of artificial variance. This work demonstrates significant individual anatomical variability of the cd-Vim that atlas-based coordinate targeting fails to capture. This variability was not related to any methodological confound tested. Lateralisation of cerebellar functions, such as speech, may contribute to the observed asymmetry. Tractography-based methods seem sensitive to individual anatomical variability that is missed by conventional neurosurgical targeting; these findings may form the basis for translational tools to improve efficacy and reduce side-effects of thalamic surgery for tremor

Diffusion MRI analysis:robust and efficient microstructure modeling

Diffusion MRI (dMRI) allows for investigating the structure of the human brain. This is useful for both scientific brain research as well as medical diagnosis. Since the raw dMRI data is not directly interpretable by humans, we use mathematical models to convert the raw dMRI data into something interpretable. These models can be computed using multiple different computational methods, each having a different trade-off in accuracy, robustness and efficiency. In this thesis we studied multiple different computational models for their usability and efficiency for dMRI modeling. In the end we provide the reader with methodological recommendations for dMRI modeling and provide a high performance GPU enabled dMRI computing platform containing all recommendations

Fusion in diffusion MRI for improved fibre orientation estimation: an application to the 3T and 7T data of the Human Connectome Project

Determining the acquisition parameters in diffusion magnetic resonance imaging (dMRI) is governed by a series of trade-offs. Images of lower resolution have less spatial specificity but higher signal to noise ratio (SNR). At the same time higher angular contrast, important for resolving complex fibre patterns, also yields lower SNR. Considering these trade-offs, the Human Connectome Project (HCP) acquires high quality dMRI data for the same subjects at different field strengths (3T and 7T), which are publically released. Due to differences in the signal behavior and in the underlying scanner hardware, the HCP 3T and 7T data have complementary features in k- and q-space. The 3T dMRI has higher angular contrast and resolution, while the 7T dMRI has higher spatial resolution. Given the availability of these datasets, we explore the idea of fusing them together with the aim of combining their benefits. We extend a previously proposed data-fusion framework and apply it to integrate both datasets from the same subject into a single joint analysis. We use a generative model for performing parametric spherical deconvolution and estimate fibre orientations by simultaneously using data acquired under different protocols. We illustrate unique features from each dataset and how they are retained after fusion. We further show that this allows us to complement benefits and improve brain connectivity analysis compared to analyzing each of the datasets individually

Subthalamic deep brain stimulation sweet spots and hyperdirect cortical connectivity in Parkinson’s disease

Objectives Firstly, to identify subthalamic region stimulation clusters that predict maximum improvement in rigidity, bradykinesia and tremor, or emergence of side-effects; and secondly, to map-out the cortical fingerprint, mediated by the hyperdirect pathways which predict maximum efficacy. Methods High angular resolution diffusion imaging in twenty patients with advanced Parkinson’s disease was acquired prior to bilateral subthalamic nucleus deep brain stimulation. All contacts were screened one-year from surgery for efficacy and side-effects at different amplitudes. Voxel-based statistical analysis of volumes of tissue activated models was used to identify significant treatment clusters. Probabilistic tractography was employed to identify cortical connectivity patterns associated with treatment efficacy. Results All patients responded well to treatment (46% mean improvement off medication UPDRS-III [p<0.0001]) without significant adverse events. Cluster corresponding to maximum improvement in tremor was in the posterior, superior and lateral portion of the nucleus. Clusters corresponding to improvement in bradykinesia and rigidity were nearer the superior border in a further medial and posterior location. The rigidity cluster extended beyond the superior border to the area of the zona incerta and Forel-H2 field. When the clusters where averaged, the coordinates of the area with maximum overall efficacy was X=-10(-9.5), Y=-13(-1) and Z=-7(-3) in MNI(AC-PC) space. Cortical connectivity to primary motor area was predictive of higher improvement in tremor; whilst that to supplementary motor area was predictive of improvement in bradykinesia and rigidity; and connectivity to prefrontal cortex was predictive of improvement in rigidity. Interpretation These findings support the presence of overlapping stimulation sites within the subthalamic nucleus and its superior border, with different cortical connectivity patterns, associated with maximum improvement in tremor, rigidity and bradykinesia