Spatial transformations of diffusion tensor magnetic resonance images

The authors address the problem of applying spatial transformations (or “image warps”) to diffusion tensor magnetic resonance images. The orientational information that these images contain must be handled appropriately when they are transformed spatially during image registration. The authors present solutions for global transformations of three-dimensional images up to 12-parameter affine complexity and indicate how their methods can be extended for higher order transformations. Several approaches are presented and tested using synthetic data. One method, the preservation of principal direction algorithm, which takes into account shearing, stretching and rigid rotation, is shown to be the most effective. Additional registration experiments are performed on human brain data obtained from a single subject, whose head was imaged in three different orientations within the scanner. All of the authors' methods improve the consistency between registered and target images over naive warping algorithms

Reward circuitry is perturbed in the absence of the serotonin transporter

The serotonin transporter (SERT) modulates the entire serotonergic system in the brain and influences both the dopaminergic and norepinephrinergic systems. These three systems are intimately involved in normal physiological functioning of the brain and implicated in numerous pathological conditions. Here we use high-resolution magnetic resonance imaging (MRI) and spectroscopy to elucidate the effects of disruption of the serotonin transporter in an animal model system: the SERT knock-out mouse. Employing manganese-enhanced MRI, we injected Mn^(2+) into the prefrontal cortex and obtained 3D MR images at specific time points in cohorts of SERT and normal mice. Statistical analysis of co-registered datasets demonstrated that active circuitry originating in the prefrontal cortex in the SERT knock-out is dramatically altered, with a bias towards more posterior areas (substantia nigra, ventral tegmental area, and Raphé nuclei) directly involved in the reward circuit. Injection site and tracing were confirmed with traditional track tracers by optical microscopy. In contrast, metabolite levels were essentially normal in the SERT knock-out by in vivo magnetic resonance spectroscopy and little or no anatomical differences between SERT knock-out and normal mice were detected by MRI. These findings point to modulation of the limbic cortical–ventral striatopallidal by disruption of SERT function. Thus, molecular disruptions of SERT that produce behavioral changes also alter the functional anatomy of the reward circuitry in which all the monoamine systems are involved

Bayesian Estimation of White Matter Atlas from High Angular Resolution Diffusion Imaging

We present a Bayesian probabilistic model to estimate the brain white matter atlas from high angular resolution diffusion imaging (HARDI) data. This model incorporates a shape prior of the white matter anatomy and the likelihood of individual observed HARDI datasets. We first assume that the atlas is generated from a known hyperatlas through a flow of diffeomorphisms and its shape prior can be constructed based on the framework of large deformation diffeomorphic metric mapping (LDDMM). LDDMM characterizes a nonlinear diffeomorphic shape space in a linear space of initial momentum uniquely determining diffeomorphic geodesic flows from the hyperatlas. Therefore, the shape prior of the HARDI atlas can be modeled using a centered Gaussian random field (GRF) model of the initial momentum. In order to construct the likelihood of observed HARDI datasets, it is necessary to study the diffeomorphic transformation of individual observations relative to the atlas and the probabilistic distribution of orientation distribution functions (ODFs). To this end, we construct the likelihood related to the transformation using the same construction as discussed for the shape prior of the atlas. The probabilistic distribution of ODFs is then constructed based on the ODF Riemannian manifold. We assume that the observed ODFs are generated by an exponential map of random tangent vectors at the deformed atlas ODF. Hence, the likelihood of the ODFs can be modeled using a GRF of their tangent vectors in the ODF Riemannian manifold. We solve for the maximum a posteriori using the Expectation-Maximization algorithm and derive the corresponding update equations. Finally, we illustrate the HARDI atlas constructed based on a Chinese aging cohort of 94 adults and compare it with that generated by averaging the coefficients of spherical harmonics of the ODF across subjects

Scanner Invariant Representations for Diffusion MRI Harmonization

Purpose: In the present work we describe the correction of diffusion-weighted MRI for site and scanner biases using a novel method based on invariant representation. Theory and Methods: Pooled imaging data from multiple sources are subject to variation between the sources. Correcting for these biases has become very important as imaging studies increase in size and multi-site cases become more common. We propose learning an intermediate representation invariant to site/protocol variables, a technique adapted from information theory-based algorithmic fairness; by leveraging the data processing inequality, such a representation can then be used to create an image reconstruction that is uninformative of its original source, yet still faithful to underlying structures. To implement this, we use a deep learning method based on variational auto-encoders (VAE) to construct scanner invariant encodings of the imaging data. Results: To evaluate our method, we use training data from the 2018 MICCAI Computational Diffusion MRI (CDMRI) Challenge Harmonization dataset. Our proposed method shows improvements on independent test data relative to a recently published baseline method on each subtask, mapping data from three different scanning contexts to and from one separate target scanning context. Conclusion: As imaging studies continue to grow, the use of pooled multi-site imaging will similarly increase. Invariant representation presents a strong candidate for the harmonization of these data

Random noise in Diffusion Tensor Imaging, its Destructive Impact and Some Corrections

The empirical origin of random noise is described, its influence on DTI variables is illustrated by a review of numerical and in vivo studies supplemented by new simulations investigating high noise levels. A stochastic model of noise propagation is presented to structure noise impact in DTI. Finally, basics of voxelwise and spatial denoising procedures are presented. Recent denoising procedures are reviewed and consequences of the stochastic model for convenient denoising strategies are discussed

Diffeomorphic Metric Mapping of High Angular Resolution Diffusion Imaging based on Riemannian Structure of Orientation Distribution Functions

In this paper, we propose a novel large deformation diffeomorphic registration algorithm to align high angular resolution diffusion images (HARDI) characterized by orientation distribution functions (ODFs). Our proposed algorithm seeks an optimal diffeomorphism of large deformation between two ODF fields in a spatial volume domain and at the same time, locally reorients an ODF in a manner such that it remains consistent with the surrounding anatomical structure. To this end, we first review the Riemannian manifold of ODFs. We then define the reorientation of an ODF when an affine transformation is applied and subsequently, define the diffeomorphic group action to be applied on the ODF based on this reorientation. We incorporate the Riemannian metric of ODFs for quantifying the similarity of two HARDI images into a variational problem defined under the large deformation diffeomorphic metric mapping (LDDMM) framework. We finally derive the gradient of the cost function in both Riemannian spaces of diffeomorphisms and the ODFs, and present its numerical implementation. Both synthetic and real brain HARDI data are used to illustrate the performance of our registration algorithm

Altered Neurocircuitry in the Dopamine Transporter Knockout Mouse Brain

The plasma membrane transporters for the monoamine neurotransmitters dopamine, serotonin, and norepinephrine modulate the dynamics of these monoamine neurotransmitters. Thus, activity of these transporters has significant consequences for monoamine activity throughout the brain and for a number of neurological and psychiatric disorders. Gene knockout (KO) mice that reduce or eliminate expression of each of these monoamine transporters have provided a wealth of new information about the function of these proteins at molecular, physiological and behavioral levels. In the present work we use the unique properties of magnetic resonance imaging (MRI) to probe the effects of altered dopaminergic dynamics on meso-scale neuronal circuitry and overall brain morphology, since changes at these levels of organization might help to account for some of the extensive pharmacological and behavioral differences observed in dopamine transporter (DAT) KO mice. Despite the smaller size of these animals, voxel-wise statistical comparison of high resolution structural MR images indicated little morphological change as a consequence of DAT KO. Likewise, proton magnetic resonance spectra recorded in the striatum indicated no significant changes in detectable metabolite concentrations between DAT KO and wild-type (WT) mice. In contrast, alterations in the circuitry from the prefrontal cortex to the mesocortical limbic system, an important brain component intimately tied to function of mesolimbic/mesocortical dopamine reward pathways, were revealed by manganese-enhanced MRI (MEMRI). Analysis of co-registered MEMRI images taken over the 26 hours after introduction of Mn^(2+) into the prefrontal cortex indicated that DAT KO mice have a truncated Mn^(2+) distribution within this circuitry with little accumulation beyond the thalamus or contralateral to the injection site. By contrast, WT littermates exhibit Mn^(2+) transport into more posterior midbrain nuclei and contralateral mesolimbic structures at 26 hr post-injection. Thus, DAT KO mice appear, at this level of anatomic resolution, to have preserved cortico-striatal-thalamic connectivity but diminished robustness of reward-modulating circuitry distal to the thalamus. This is in contradistinction to the state of this circuitry in serotonin transporter KO mice where we observed more robust connectivity in more posterior brain regions using methods identical to those employed here

Connectivity-based parcellation of the thalamus explains specific cognitive and behavioural symptoms in patients with bilateral thalamic infarct

A novel approach based on diffusion tractography was used here to characterise the cortico-thalamic connectivity in two patients, both presenting with an isolated bilateral infarct in the thalamus, but exhibiting partially different cognitive and behavioural profiles. Both patients (G.P. and R.F.) had a pervasive deficit in episodic memory, but only one of them (R.F.) suffered also from a dysexecutive syndrome. Both patients had an MRI scan at 3T, including a T1-weighted volume. Their lesions were manually segmented. T1-volumes were normalised to standard space, and the same transformations were applied to the lesion masks. Nineteen healthy controls underwent a diffusion-tensor imaging (DTI) scan. Their DTI data were normalised to standard space and averaged. An atlas of Brodmann areas was used to parcellate the prefrontal cortex. Probabilistic tractography was used to assess the probability of connection between each voxel of the thalamus and a set of prefrontal areas. The resulting map of corticothalamic connections was superimposed onto the patients' lesion masks, to assess whether the location of the thalamic lesions in R.F. (but not in G. P.) implied connections with prefrontal areas involved in dysexecutive syndromes. In G.P., the lesion fell within areas of the thalamus poorly connected with prefrontal areas, showing only a modest probability of connection with the anterior cingulate cortex (ACC). Conversely, R.F.'s lesion fell within thalamic areas extensively connected with the ACC bilaterally, with the right dorsolateral prefrontal cortex, and with the left supplementary motor area. Despite a similar, bilateral involvement of the thalamus, the use of connectivity-based segmentation clarified that R.F.'s lesions only were located within nuclei highly connected with the prefrontal cortical areas, thus explaining the patient's frontal syndrome. This study confirms that DTI tractography is a useful tool to examine in vivo the effect of focal lesions on interconnectivity brain patterns