Dissection of the Temporofrontal Extreme Capsule Fasciculus Using Diffusion MRI Tractography and Association with Lexical Retrieval.

The well-known arcuate fasciculus that connects the posterior superior temporal region with the language production region in the ventrolateral frontal cortex constitutes the classic peri-Sylvian dorsal stream of language. A second temporofrontal white matter tract connects ventrally the anterior to intermediate lateral temporal cortex with frontal areas via the extreme capsule. This temporofrontal extreme capsule fasciculus (TFexcF) constitutes the ventral stream of language processing. The precise origin, course, and termination of this pathway has been examined in invasive tract tracing studies in macaque monkeys, but there have been no standard protocols for its reconstruction in the human brain using diffusion imaging tractography. Here we provide a protocol for the dissection of the TFexcF in vivo in the human brain using diffusion magnetic resonance imaging (MRI) tractography which provides a solid basis for exploring its functional role. A key finding of the current dissection protocol is the demonstration that the TFexcF is left hemisphere lateralized. Furthermore, using the present dissection protocol, we demonstrate that the TFexcF is related to lexical retrieval scores measured with the category fluency test, in contrast to the classical arcuate fasciculus (the dorsal language pathway) that was also dissected and was related to sentence repetition.

Neurophysiological basis of contrast dependent BOLD orientation tuning.

Recent work in early visual cortex of humans has shown that the BOLD signal exhibits contrast dependent orientation tuning, with an inverse oblique effect (oblique > cardinal) at high contrast and a horizontal effect (vertical > horizontal) at low contrast. This finding is at odds with decades of neurophysiological research demonstrating contrast invariant orientation tuning in primate visual cortex, yet the source of this discrepancy is unclear. We hypothesized that contrast dependent BOLD orientation tuning may arise due to contrast dependent influences of feedforward (FF) and feedback (FB) synaptic activity, indexed through gamma and alpha rhythms, respectively. To quantify this, we acquired EEG and BOLD in healthy humans to generate and compare orientation tuning curves across all neural frequency bands with BOLD. As expected, BOLD orientation selectivity in V1 was contrast dependent, preferring oblique orientations at high contrast and vertical at low contrast. On the other hand, EEG orientation tuning was contrast invariant, though frequency-specific, with an inverse-oblique effect in the gamma band (FF) and a horizontal effect in the alpha band (FB). Therefore, high-contrast BOLD orientation tuning closely matched FF activity, while at low contrast, BOLD best resembled FB orientation tuning. These results suggest that contrast dependent BOLD orientation tuning arises due to the reduced contribution of FF input to overall neurophysiological activity at low contrast, shifting BOLD orientation tuning towards the orientation preferences of FB at low contrast.

3D-SSF: A bio-inspired approach for dynamic multi-subject clustering of white matter tracts.

There is growing interest in the study of white matter (WM) variation across subjects, and in particular the analysis of specific WM bundles, to better understand brain development and aging, as well as to improve early detection of some diseases. Several WM multi-subject clustering methods have been proposed to study WM bundles. These methods aim to overcome the complexity of the problem, which includes the huge size of the WM tractography datasets generated from multiple subjects, the existence of various streamlines with different positions, lengths and geometric forms, as well as the presence of outliers. However, the current methods are not sufficiently flexible to address all of these constraints. Here we introduce a novel dynamic multi-subject clustering framework based on a distributed multiagent implementation of the Multiple Species Flocking model, that we name 3D-Streamlines Stream Flocking (3D-SSF). Specifically, we consider streamlines from different subjects as data streams, and each streamline is assigned to a mobile agent. Agents work together following flocking rules in order to form a flock. Thanks to a similarity function, the agents that are associated with similar streamlines form a flock, whereas the agents that are associated with dissimilar streamlines are considered outliers. We use various experiments performed on noisy synthetic and real human brain data to validate 3D-SSF and demonstrate that it is more efficient and robust to outliers compared to other classical approaches. 3D-SSF is able to extract WM bundles at a population level, while considering WM variation across subjects and eliminating outlier streamlines.

Denoising and fast diffusion imaging with physically constrained sparse dictionary learning.

Diffusion-weighted imaging (DWI) allows imaging the geometry of water diffusion in biological tissues. However, DW images are noisy at high b-values and acquisitions are slow when using a large number of measurements, such as in Diffusion Spectrum Imaging (DSI). This work aims to denoise DWI and reduce the number of required measurements, while maintaining data quality. To capture the structure of DWI data, we use sparse dictionary learning constrained by the physical properties of the signal: symmetry and positivity. The method learns a dictionary of diffusion profiles on all the DW images at the same time and then scales to full brain data. Its performance is investigated with simulations and two real DSI datasets. We obtain better signal estimates from noisy measurements than by applying mirror symmetry through the q-space origin, Gaussian denoising or state-of-the-art non-local means denoising. Using a high-resolution dictionary learnt on another subject, we show that we can reduce the number of images acquired while still generating high resolution DSI data. Using dictionary learning, one can denoise DW images effectively and perform faster acquisitions. Higher b-value acquisitions and DSI techniques are possible with approximately 40 measurements. This opens important perspectives for the connectomics community using DSI.

Robust clustering of massive tractography datasets.

This paper presents a clustering method that detects the fiber bundles embedded in any MR-diffusion based tractography dataset. Our method can be seen as a compressing operation, capturing the most meaningful information enclosed in the fiber dataset. For the sake of efficiency, part of the analysis is based on clustering the white matter (WM) voxels rather than the fibers. The resulting regions of interest are used to define subset of fibers that are subdivided further into consistent bundles using a clustering of the fiber extremities. The dataset is reduced from more than one million fiber tracts to about two thousand fiber bundles. Validations are provided using simulated data and a physical phantom. We see our approach as a crucial preprocessing step before further analysis of huge fiber datasets. An important application will be the inference of detailed models of the subdivisions of white matter pathways and the mapping of the main U-fiber bundles.

Spherical wavelet transform for ODF sharpening.

The choice of local HARDI reconstruction technique is crucial for discerning multiple fiber orientations, which is itself of substantial importance for tractography, and reliable and accurate assessment of white matter fiber geometry. Due to the complexity of the diffusion process and its milieu, distinct diffusion compartments can have different frequency signatures, making the HARDI signal spread over multiple frequency bands. Therefore, we put forth the idea of multiscale analysis with localized basis functions, ensuring that different frequency ranges are probed. With the aim of truthful recovery of fiber orientations, we reconstruct the orientation distribution function (ODF), by incorporating a spherical wavelet transform (SWT) into the Funk-Radon transform. First, we apply and validate our proposed SWT method on real physical phantoms emulating fiber bundle crossings. Then, we apply the SWT method to a real brain data set. The analysis of the real data set suggests that different angular frequencies may capture different information, thus stressing the importance of multiscale analysis. For both phantom and real data, we compare the SWT reconstruction with state-of-the-art q-ball imaging and spherical deconvolution reconstruction methods. We demonstrate the algorithm efficiency in diffusion ODF denoising and sharpening that is of particular importance for applications to fiber tracking (especially for probabilistic approaches), and brain connectome mapping. Also, the algorithm results in considerable data compression that could prove beneficial in applications to fiber bundle segmentation, and for HARDI based white matter morphometry methods.

Mathematical methods for diffusion MRI processing.

In this article, we review recent mathematical models and computational methods for the processing of diffusion Magnetic Resonance Images, including state-of-the-art reconstruction of diffusion models, cerebral white matter connectivity analysis, and segmentation techniques. We focus on Diffusion Tensor Images (DTI) and Q-Ball Images (QBI).

Detection of multiple pathways in the spinal cord using q-ball imaging.

Magnetic resonance diffusion tensor imaging (DTI) has been extensively applied to the spinal cord for depicting its architecture and for assessing its integrity following spinal lesions. However, DTI is limited in representing complex white matter architecture, notably in the presence of crossing fibres. Recently, q-ball imaging (QBI) has been proposed as a new method for recovering complex white matter architecture. We applied this technique to both ex vivo and in vivo spinal cords of cats using a 3T scanner. For the purpose of comparison, gradients have been applied in 55 and 100 encoding directions and b-values varied from 800 to 3000 s/mm(2). As a result, QBI was able to retrieve crossing fibre information, where the DTI approach was constrained in a unique diffusion direction. To our knowledge, this is the first study demonstrating the benefits of QBI for detecting the presence of longitudinal, commissural and dorso-ventral fibres in the spinal cord. It is a first step towards in vivo characterization of the healthy and injured human spinal cord using high angular resolution diffusion imaging and QBI.

Validation of vessel-based registration for correction of brain shift.

The displacement and deformation of brain tissue is a major source of error in image-guided neurosurgery systems. We have designed and implemented a method to detect and correct brain shift using pre-operative MR images and intraoperative Doppler ultrasound data and present its validation with both real and simulated data. The algorithm uses segmented vessels from both modalities, and estimates the deformation using a modified version of the iterative closest point (ICP) algorithm. We use the least trimmed squares (LTS) to reduce the number of outliers in the point matching procedure. These points are used to drive a thin-plate spline transform to achieve non-linear registration. Validation was completed in two parts. First, the technique was tested and validated using realistic simulations where the results were compared to the known deformation. The registration technique recovered 75% of the deformation in the region of interest accounting for deformations as large as 20 mm. Second, we performed a PVA-cryogel phantom study where both MR and ultrasound images of the phantom were obtained for three different deformations. The registration results based on MR data were used as a gold standard to evaluate the performance of the ultrasound based registration. On average, deformations of 7.5 mm magnitude were corrected to within 1.6 mm for the ultrasound based registration and 1.07 mm for the MR based registration.