Search

Diffu0sion-weighted magnetic resonance imaging (dMRI) is a non-invasive imaging technique that provides information about the barriers to the diffusion of water molecules in tissue. In the brain, this information can be used in several important ways, including to examine tissue abnormalities associated with brain disorders and to infer anatomical connectivity and the organization of white matter bundles through the use of tractography algorithms. However, dMRI also presents certain challenges. For example, historically, the biological validation of tractography models has shown only moderate correlations with anatomical connectivity as determined through invasive tract-tracing studies. Some of the factors contributing to such issues are low spatial resolution, low signal-to-noise ratios, and long scan times required for high-quality data, along with modeling challenges like complex fiber crossing patterns. Leveraging the capabilities provided by an ultra-high field scanner combined with denoising, we have acquired whole-brain, 0.58 mm isotropic resolution dMRI with a 2D-single shot echo planar imaging sequence on a 10.5 Tesla scanner in anesthetized macaques. These data produced high-quality tractograms and maps of scalar diffusion metrics in white matter. This work demonstrates the feasibility and motivation for in-vivo dMRI studies seeking to benefit from ultra-high fields.

The success of deep brain stimulation (DBS) surgeries for the treatment of movement disorders relies on the accurate placement of an electrode within the motor portion of subcortical brain targets. However, the high number of electrodes requiring relocation indicates that today's methods do not ensure sufficient accuracy for all patients. Here, with the goal of aiding DBS targeting, we use 7 Tesla (T) MRI data to identify the functional territories and parcellate the globus pallidus pars interna (GPi) into motor, associative and limbic regions in individual subjects. 7 T MRI scans were performed in seventeen patients (prior to DBS surgery) and one healthy control. Tractography-based parcellation of each patient's GPi was performed. The cortex was divided into four masks representing motor, limbic, associative and “other” regions. Given that no direct connections between the GPi and the cortex have been shown to exist, the parcellation was carried out in two steps: 1) The thalamus was parcellated based on the cortical targets, 2) The GPi was parcellated using the thalamus parcels derived from step 1. Reproducibility, via repeated scans of a healthy subject, and validity of the findings, using different anatomical pathways for parcellation, were assessed. Lastly, post-operative imaging data was used to validate and determine the clinical relevance of the parcellation. The organization of the functional territories of the GPi observed in our individual patient population agrees with that previously reported in the literature: the motor territory was located posterolaterally, followed anteriorly by the associative region, and further antero-ventrally by the limbic territory. While this organizational pattern was observed across patients, there was considerable variability among patients. The organization of the functional territories of the GPi was remarkably reproducible in intra-subject scans. Furthermore, the organizational pattern was observed consistently by performing the parcellation of the GPi via the thalamus and via a different pathway, going through the striatum. Finally, the active therapeutic contact of the DBS electrode, identified with a combination of post-operative imaging and post-surgery DBS programming, overlapped with the high-probability “motor” region of the GPi as defined by imaging-based methods. The consistency, validity, and clinical relevance of our findings have the potential for improving DBS targeting, by increasing patient-specific knowledge of subregions of the GPi to be targeted or avoided, at the stage of surgical planning, and later, at the stage when stimulation is adjusted.

We present a sparse Bayesian unmixing algorithm BusineX: Bayesian Unmixing for Sparse Inference-based Estimation of Fiber Crossings (X), for estimation of white matter fiber parameters from compressed (under-sampled) diffusion MRI (dMRI) data. BusineX combines compressive sensing with linear unmixing and introduces sparsity to the previously proposed multiresolution data fusion algorithm RubiX, resulting in a method for improved reconstruction, especially from data with lower number of diffusion gradients. We formulate the estimation of fiber parameters as a sparse signal recovery problem and propose a linear unmixing framework with sparse Bayesian learning for the recovery of sparse signals, the fiber orientations and volume fractions. The data is modeled using a parametric spherical deconvolution approach and represented using a dictionary created with the exponential decay components along different possible diffusion directions. Volume fractions of fibers along these directions define the dictionary weights. The proposed sparse inference, which is based on the dictionary representation, considers the sparsity of fiber populations and exploits the spatial redundancy in data representation, thereby facilitating inference from under-sampled q-space. The algorithm improves parameter estimation from dMRI through data-dependent local learning of hyperparameters, at each voxel and for each possible fiber orientation, that moderate the strength of priors governing the parameter variances. Experimental results on synthetic and in-vivo data show improved accuracy with a lower uncertainty in fiber parameter estimates. BusineX resolves a higher number of second and third fiber crossings. For under-sampled data, the algorithm is also shown to produce more reliable estimates.

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 (7 T 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 7 T image acquisitions for the HCP that allow us to efficiently obtain high-quality, high-resolution whole brain in-vivo dMRI data at 7 T. These data show spatial details typically seen only in ex-vivo studies and complement already very high quality 3 T HCP data in the same subjects. The advances are the result of intensive pilot studies aimed at mitigating the limitations of dMRI at 7 T. The data quality and methods described here are representative of the datasets that will be made freely available to the community in 2015.

High angular resolution diffusion imaging (HARDI) has become an important technique for imaging complex oriented structures in the brain and other anatomical tissues. This has motivated the recent development of several methods for computing the orientation probability density function (PDF) at each voxel. However, much less work has been done on developing techniques for filtering, interpolation, averaging and principal geodesic analysis of orientation PDF fields. In this paper, we present a Riemannian framework for performing such operations. The proposed framework does not require that the orientation PDFs be represented by any fixed parameterization, such as a mixture of von Mises-Fisher distributions or a spherical harmonic expansion. Instead, we use a nonparametric representation of the orientation PDF. We exploit the fact that under the square-root re-parameterization, the space of orientation PDFs forms a Riemannian manifold: the positive orthant of the unit Hilbert sphere. We show that various orientation PDF processing operations, such as filtering, interpolation, averaging and principal geodesic analysis, may be posed as optimization problems on the Hilbert sphere, and can be solved using Riemannian gradient descent. We illustrate these concepts with numerous experiments on synthetic, phantom and real datasets. We show their application to studying left/right brain asymmetries.

This paper presents a coupled thermo-hydro-mechanical (THM) model enriched with a buckling-type criterion for progressive spalling. In the first part of the paper, a general fully coupled multi-phase THM model describing the behaviour of concrete at moderate and high temperatures is presented. Then the spalling criterion and its numerical implementation in the framework of the finite element method are presented. Finally, a simple 1D numerical example will illustrate the effectiveness of the implemented numerical approach.