Estimation of Fiber Orientations Using Neighborhood Information

Data from diffusion magnetic resonance imaging (dMRI) can be used to reconstruct fiber tracts, for example, in muscle and white matter. Estimation of fiber orientations (FOs) is a crucial step in the reconstruction process and these estimates can be corrupted by noise. In this paper, a new method called Fiber Orientation Reconstruction using Neighborhood Information (FORNI) is described and shown to reduce the effects of noise and improve FO estimation performance by incorporating spatial consistency. FORNI uses a fixed tensor basis to model the diffusion weighted signals, which has the advantage of providing an explicit relationship between the basis vectors and the FOs. FO spatial coherence is encouraged using weighted l1-norm regularization terms, which contain the interaction of directional information between neighbor voxels. Data fidelity is encouraged using a squared error between the observed and reconstructed diffusion weighted signals. After appropriate weighting of these competing objectives, the resulting objective function is minimized using a block coordinate descent algorithm, and a straightforward parallelization strategy is used to speed up processing. Experiments were performed on a digital crossing phantom, ex vivo tongue dMRI data, and in vivo brain dMRI data for both qualitative and quantitative evaluation. The results demonstrate that FORNI improves the quality of FO estimation over other state of the art algorithms.Comment: Journal paper accepted in Medical Image Analysis. 35 pages and 16 figure

Diffeomorphic Metric Mapping and Probabilistic Atlas Generation of Hybrid Diffusion Imaging based on BFOR Signal Basis

We propose a large deformation diffeomorphic metric mapping algorithm to align multiple b-value diffusion weighted imaging (mDWI) data, specifically acquired via hybrid diffusion imaging (HYDI), denoted as LDDMM-HYDI. We then propose a Bayesian model for estimating the white matter atlas from HYDIs. We adopt the work given in Hosseinbor et al. (2012) and represent the q-space diffusion signal with the Bessel Fourier orientation reconstruction (BFOR) signal basis. The BFOR framework provides the representation of mDWI in the q-space and thus reduces memory requirement. In addition, since the BFOR signal basis is orthonormal, the L2 norm that quantifies the differences in the q-space signals of any two mDWI datasets can be easily computed as the sum of the squared differences in the BFOR expansion coefficients. In this work, we show that the reorientation of the $q$-space signal due to spatial transformation can be easily defined on the BFOR signal basis. We incorporate the BFOR signal basis into the LDDMM framework and derive the gradient descent algorithm for LDDMM-HYDI with explicit orientation optimization. Additionally, we extend the previous Bayesian atlas estimation framework for scalar-valued images to HYDIs and derive the expectation-maximization algorithm for solving the HYDI atlas estimation problem. Using real HYDI datasets, we show the Bayesian model generates the white matter atlas with anatomical details. Moreover, we show that it is important to consider the variation of mDWI reorientation due to a small change in diffeomorphic transformation in the LDDMM-HYDI optimization and to incorporate the full information of HYDI for aligning mDWI

Altered white matter microstructure is associated with social cognition and psychotic symptoms in 22q11.2 microdeletion syndrome.

22q11.2 Microdeletion Syndrome (22q11DS) is a highly penetrant genetic mutation associated with a significantly increased risk for psychosis. Aberrant neurodevelopment may lead to inappropriate neural circuit formation and cerebral dysconnectivity in 22q11DS, which may contribute to symptom development. Here we examined: (1) differences between 22q11DS participants and typically developing controls in diffusion tensor imaging (DTI) measures within white matter tracts; (2) whether there is an altered age-related trajectory of white matter pathways in 22q11DS; and (3) relationships between DTI measures, social cognition task performance, and positive symptoms of psychosis in 22q11DS and typically developing controls. Sixty-four direction diffusion weighted imaging data were acquired on 65 participants (36 22q11DS, 29 controls). We examined differences between 22q11DS vs. controls in measures of fractional anisotropy (FA), axial diffusivity (AD), and radial diffusivity (RD), using both a voxel-based and region of interest approach. Social cognition domains assessed were: Theory of Mind and emotion recognition. Positive symptoms were assessed using the Structured Interview for Prodromal Syndromes. Compared to typically developing controls, 22q11DS participants showed significantly lower AD and RD in multiple white matter tracts, with effects of greatest magnitude for AD in the superior longitudinal fasciculus. Additionally, 22q11DS participants failed to show typical age-associated changes in FA and RD in the left inferior longitudinal fasciculus. Higher AD in the left inferior fronto-occipital fasciculus (IFO) and left uncinate fasciculus was associated with better social cognition in 22q11DS and controls. In contrast, greater severity of positive symptoms was associated with lower AD in bilateral regions of the IFO in 22q11DS. White matter microstructure in tracts relevant to social cognition is disrupted in 22q11DS, and may contribute to psychosis risk

Tracking dynamic interactions between structural and functional connectivity : a TMS/EEG-dMRI study

Transcranial magnetic stimulation (TMS) in combination with neuroimaging techniques allows to measure the effects of a direct perturbation of the brain. When coupled with high-density electroencephalography (TMS/hd-EEG), TMS pulses revealed electrophysiological signatures of different cortical modules in health and disease. However, the neural underpinnings of these signatures remain unclear. Here, by applying multimodal analyses of cortical response to TMS recordings and diffusion magnetic resonance imaging (dMRI) tractography, we investigated the relationship between functional and structural features of different cortical modules in a cohort of awake healthy volunteers. For each subject, we computed directed functional connectivity interactions between cortical areas from the source-reconstructed TMS/hd-EEG recordings and correlated them with the correspondent structural connectivity matrix extracted from dMRI tractography, in three different frequency bands (alpha, beta, gamma) and two sites of stimulation (left precuneus and left premotor). Each stimulated area appeared to mainly respond to TMS by being functionally elicited in specific frequency bands, that is, beta for precuneus and gamma for premotor. We also observed a temporary decrease in the whole-brain correlation between directed functional connectivity and structural connectivity after TMS in all frequency bands. Notably, when focusing on the stimulated areas only, we found that the structure-function correlation significantly increases over time in the premotor area controlateral to TMS. Our study points out the importance of taking into account the major role played by different cortical oscillations when investigating the mechanisms for integration and segregation of information in the human brain

fMRI-Targeted High-Angular Resolution Diffusion MR Tractography to Identify Functional Language Tracts in Healthy Controls and Glioma Patients

Background MR Tractography enables non-invasive preoperative depiction of language subcortical tracts, which is crucial for the presurgical work-up of brain tumors; however, it cannot evaluate the exact function of the fibers. Purpose A systematic pipeline was developed to combine tractography reconstruction of language fiber bundles, based on anatomical landmarks (Anatomical-T), with language fMRI cortical activations. A fMRI-targeted Tractography (fMRI-T) was thus obtained, depicting the subsets of the anatomical tracts whose endpoints are located inside a fMRI activation. We hypothesized that fMRI-T could provide additional functional information regarding the subcortical structures, better reflecting the eloquent white matter structures identified intraoperatively. Methods Both Anatomical-T and fMRI-T of language fiber tracts were performed on 16 controls and preoperatively on 16 patients with left-hemisphere brain tumors, using a q-ball residual bootstrap algorithm based on High Angular Resolution Diffusion Imaging (HARDI) datasets (b = 3000 s/mm(2); 60 directions); fMRI ROIs were obtained using picture naming, verbal fluency, and auditory verb generation tasks. In healthy controls, normalized MNI atlases of fMRI-T and Anatomical-T were obtained. In patients, the surgical resection of the tumor was pursued by identifying eloquent structures with intraoperative direct electrical stimulation mapping and extending surgery to the functional boundaries. Post-surgical MRI allowed to identify Anatomical-T and fMRI-T non-eloquent portions removed during the procedure. Results MNI Atlases showed that fMRI-T is a subset of Anatomical-T, and that different task-specific fMRI-T involve both shared subsets and task-specific subsets - e.g., verbal fluency fMRI-T strongly involves dorsal frontal tracts, consistently with the phonogical-articulatory features of this task. A quantitative analysis in patients revealed that Anatomical-T removed portions of AF-SLF and IFOF were significantly greater than verbal fluency fMRI-T ones, suggesting that fMRI-T is a more specific approach. In addition, qualitative analyses showed that fMRI-T AF-SLF and IFOF predict the exact functional limits of resection with increased specificity when compared to Anatomical-T counterparts, especially the superior frontal portion of IFOF, in a subcohort of patients. Conclusion These results suggest that performing fMRI-T in addition to the 'classic' Anatomical-T may be useful in a preoperative setting to identify the 'high-risk subsets' that should be spared during the surgical procedure