Fuzzy Fibers: Uncertainty in dMRI Tractography

Fiber tracking based on diffusion weighted Magnetic Resonance Imaging (dMRI) allows for noninvasive reconstruction of fiber bundles in the human brain. In this chapter, we discuss sources of error and uncertainty in this technique, and review strategies that afford a more reliable interpretation of the results. This includes methods for computing and rendering probabilistic tractograms, which estimate precision in the face of measurement noise and artifacts. However, we also address aspects that have received less attention so far, such as model selection, partial voluming, and the impact of parameters, both in preprocessing and in fiber tracking itself. We conclude by giving impulses for future research

An Open-Source Tool for Anisotropic Radiation Therapy Planning in Neuro-oncology Using DW-MRI Tractography.

There is evidence from histopathological studies that glioma tumor cells migrate preferentially along large white matter bundles. If the peritumoral white matter structures can be used to predict the likely trajectory of migrating tumor cells outside of the surgical margin, then this information could be used to inform the delineation of radiation therapy (RT) targets. In theory, an anisotropic expansion that takes large white matter bundle anatomy into account may maximize the chances of treating migrating cancer cells and minimize the amount of brain tissue exposed to high doses of ionizing radiation. Diffusion-weighted MRI (DW-MRI) can be used in combination with fiber tracking algorithms to model the trajectory of large white matter pathways using the direction and magnitude of water movement in tissue. The method presented here is a tool for translating a DW-MRI fiber tracking (tractography) dataset into a white matter path length (WMPL) map that assigns each voxel the shortest distance along a streamline back to a specified region of interest (ROI). We present an open-source WMPL tool, implemented in the package Diffusion Imaging in Python (DIPY), and code to convert the resulting WMPL map to anisotropic contours for RT in a commercial treatment planning system. This proof-of-concept lays the groundwork for future studies to evaluate the clinical value of incorporating tractography modeling into treatment planning

Improving Fiber Alignment in HARDI by Combining Contextual PDE Flow with Constrained Spherical Deconvolution

We propose two strategies to improve the quality of tractography results computed from diffusion weighted magnetic resonance imaging (DW-MRI) data. Both methods are based on the same PDE framework, defined in the coupled space of positions and orientations, associated with a stochastic process describing the enhancement of elongated structures while preserving crossing structures. In the first method we use the enhancement PDE for contextual regularization of a fiber orientation distribution (FOD) that is obtained on individual voxels from high angular resolution diffusion imaging (HARDI) data via constrained spherical deconvolution (CSD). Thereby we improve the FOD as input for subsequent tractography. Secondly, we introduce the fiber to bundle coherence (FBC), a measure for quantification of fiber alignment. The FBC is computed from a tractography result using the same PDE framework and provides a criterion for removing the spurious fibers. We validate the proposed combination of CSD and enhancement on phantom data and on human data, acquired with different scanning protocols. On the phantom data we find that PDE enhancements improve both local metrics and global metrics of tractography results, compared to CSD without enhancements. On the human data we show that the enhancements allow for a better reconstruction of crossing fiber bundles and they reduce the variability of the tractography output with respect to the acquisition parameters. Finally, we show that both the enhancement of the FODs and the use of the FBC measure on the tractography improve the stability with respect to different stochastic realizations of probabilistic tractography. This is shown in a clinical application: the reconstruction of the optic radiation for epilepsy surgery planning

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

Probing white-matter microstructure with higher-order diffusion tensors and susceptibility tensor MRI.

Diffusion MRI has become an invaluable tool for studying white matter microstructure and brain connectivity. The emergence of quantitative susceptibility mapping and susceptibility tensor imaging (STI) has provided another unique tool for assessing the structure of white matter. In the highly ordered white matter structure, diffusion MRI measures hindered water mobility induced by various tissue and cell membranes, while susceptibility sensitizes to the molecular composition and axonal arrangement. Integrating these two methods may produce new insights into the complex physiology of white matter. In this study, we investigated the relationship between diffusion and magnetic susceptibility in the white matter. Experiments were conducted on phantoms and human brains in vivo. Diffusion properties were quantified with the diffusion tensor model and also with the higher order tensor model based on the cumulant expansion. Frequency shift and susceptibility tensor were measured with quantitative susceptibility mapping and susceptibility tensor imaging. These diffusion and susceptibility quantities were compared and correlated in regions of single fiber bundles and regions of multiple fiber orientations. Relationships were established with similarities and differences identified. It is believed that diffusion MRI and susceptibility MRI provide complementary information of the microstructure of white matter. Together, they allow a more complete assessment of healthy and diseased brains

Ensemble tractography

Fiber tractography uses diffusion MRI to estimate the trajectory and cortical projection zones of white matter fascicles in the living human brain. There are many different tractography algorithms and each requires the user to set several parameters, such as curvature threshold. Choosing a single algorithm with a specific parameters sets poses two challenges. First, different algorithms and parameter values produce different results. Second, the optimal choice of algorithm and parameter value may differ between different white matter regions or different fascicles, subjects, and acquisition parameters. We propose using ensemble methods to reduce algorithm and parameter dependencies. To do so we separate the processes of fascicle generation and evaluation. Specifically, we analyze the value of creating optimized connectomes by systematically combining candidate fascicles from an ensemble of algorithms (deterministic and probabilistic) and sweeping through key parameters (curvature and stopping criterion). The ensemble approach leads to optimized connectomes that provide better cross-validatedprediction error of the diffusion MRI data than optimized connectomes generated using the singlealgorithms or parameter set. Furthermore, the ensemble approach produces connectomes that contain both short- and long-range fascicles, whereas single-parameter connectomes are biased towards one or the other. In summary, a systematic ensemble tractography approach can produce connectomes that are superior to standard single parameter estimates both for predicting the diffusion measurements and estimating white matter fascicles.Fil: Takemura, Hiromasa. University of Stanford; Estados Unidos. Osaka University; JapónFil: Caiafa, César Federico. Provincia de Buenos Aires. Gobernación. Comisión de Investigaciones Científicas. Instituto Argentino de Radioastronomía. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto Argentino de Radioastronomía; ArgentinaFil: Wandell, Brian A.. University of Stanford; Estados UnidosFil: Pestilli, Franco. Indiana University; Estados Unido

Connectivity of the Superficial Muscles of the Human Perineum: A Diffusion Tensor Imaging-Based Global Tractography Study.

Despite the importance of pelvic floor muscles, significant controversy still exists about the true structural details of these muscles. We provide an objective analysis of the architecture and orientation of the superficial muscles of the perineum using a novel approach. Magnetic Resonance Diffusion Tensor Images (MR-DTI) were acquired in 10 healthy asymptomatic nulliparous women, and 4 healthy males. Global tractography was then used to generate the architecture of the muscles. Micro-CT imaging of a male cadaver was performed for validation of the fiber tracking results. Results show that muscles fibers of the external anal sphincter, from the right and left side, cross midline in the region of the perineal body to continue as transverse perinea and bulbospongiosus muscles of the opposite side. The morphology of the external anal sphincter resembles that of the number '8' or a "purse string". The crossing of muscle fascicles in the perineal body was supported by micro-CT imaging in the male subject. The superficial muscles of the perineum, and external anal sphincter are frequently damaged during child birth related injuries to the pelvic floor; we propose the use of MR-DTI based global tractography as a non-invasive imaging technique to assess damage to these muscles