Corticospinal tract modeling for neurosurgical planning by tracking through regions of peritumoral edema and crossing fibers using two-tensor unscented Kalman filter tractography

A systematic review of automated methods to perform white matter tract segmentation

White matter tract segmentation is a pivotal research area that leverages diffusion-weighted magnetic resonance imaging (dMRI) for the identification and mapping of individual white matter tracts and their trajectories. This study aims to provide a comprehensive systematic literature review on automated methods for white matter tract segmentation in brain dMRI scans. Articles on PubMed, ScienceDirect [NeuroImage, NeuroImage (Clinical), Medical Image Analysis], Scopus and IEEEXplore databases and Conference proceedings of Medical Imaging Computing and Computer Assisted Intervention Society (MICCAI) and International Symposium on Biomedical Imaging (ISBI), were searched in the range from January 2013 until September 2023. This systematic search and review identified 619 articles. Adhering to the specified search criteria using the query, "white matter tract segmentation OR fiber tract identification OR fiber bundle segmentation OR tractography dissection OR white matter parcellation OR tract segmentation," 59 published studies were selected. Among these, 27% employed direct voxel-based methods, 25% applied streamline-based clustering methods, 20% used streamline-based classification methods, 14% implemented atlas-based methods, and 14% utilized hybrid approaches. The paper delves into the research gaps and challenges associated with each of these categories. Additionally, this review paper illuminates the most frequently utilized public datasets for tract segmentation along with their specific characteristics. Furthermore, it presents evaluation strategies and their key attributes. The review concludes with a detailed discussion of the challenges and future directions in this field.

Reconstructing the somatotopic organization of the corticospinal tract remains a challenge for modern tractography methods

The corticospinal tract (CST) is a critically important white matter fiber tract in the human brain that enables control of voluntary movements of the body. The CST exhibits a somatotopic organization, which means that the motor neurons that control specific body parts are arranged in order within the CST. Diffusion magnetic resonance imaging (MRI) tractography is increasingly used to study the anatomy of the CST. However, despite many advances in tractography algorithms over the past decade, modern, state-of-the-art methods still face challenges. In this study, we compare the performance of six widely used tractography methods for reconstructing the CST and its somatotopic organization. These methods include constrained spherical deconvolution (CSD) based probabilistic (iFOD1) and deterministic (SD-Stream) methods, unscented Kalman filter (UKF) tractography methods including multi-fiber (UKF2T) and single-fiber (UKF1T) models, the generalized q-sampling imaging (GQI) based deterministic tractography method, and the TractSeg method. We investigate CST somatotopy by dividing the CST into four subdivisions per hemisphere that originate in the leg, trunk, hand, and face areas of the primary motor cortex. A quantitative and visual comparison is performed using diffusion MRI data (N = 100 subjects) from the Human Connectome Project. Quantitative evaluations include the reconstruction rate of the eight anatomical subdivisions, the percentage of streamlines in each subdivision, and the coverage of the white matter-gray matter (WM-GM) interface. CST somatotopy is further evaluated by comparing the percentage of streamlines in each subdivision to the cortical volumes for the leg, trunk, hand, and face areas. Overall, UKF2T has the highest reconstruction rate and cortical coverage. It is the only method with a significant positive correlation between the percentage of streamlines in each subdivision and the volume of the corresponding motor cortex. However, our experimental results show that all compared tractography methods are biased toward generating many trunk streamlines (ranging from 35.10% to 71.66% of total streamlines across methods). Furthermore, the coverage of the WM-GM interface in the largest motor area (face) is generally low (under 40%) for all compared tractography methods. Different tractography methods give conflicting results regarding the percentage of streamlines in each subdivision and the volume of the corresponding motor cortex, indicating that there is generally no clear relationship, and that reconstruction of CST somatotopy is still a large challenge. Overall, we conclude that while current tractography methods have made progress toward the well-known challenge of improving the reconstruction of the lateral projections of the CST, the overall problem of performing a comprehensive CST reconstruction, including clinically important projections in the lateral (hand and face areas) and medial portions (leg area), remains an important challenge for diffusion MRI tractography.

Quantitative Validation of the Correlation Between Optimized Pyramidal Tract Delineation After Brain Shift Compensation and Direct Electrical Subcortical Stimulation During Brain Tumor Surgery

It remains unclear whether tractography of pyramidal tracts is correlated with the intraoperative direct electrical subcortical stimulation (DESS), and brain shift further complicates the issue. The objective of this research is to quantitatively verify the correlation between optimized tractography (OT) of pyramidal tracts after brain shift compensation and DESS during brain tumor surgery. OT was performed for 20 patients with lesions in proximity to the pyramidal tracts based on preoperative diffusion-weighted magnetic resonance imaging. During surgery, tumor resection was guided by DESS. A total of 168 positive stimulation points and their corresponding stimulation intensity thresholds were recorded. Using the brain shift compensation algorithm based on hierarchical B-spline grids combined with a Gaussian resolution pyramid, we warped the preoperative pyramidal tract models and used receiver operating characteristic (ROC) curves to investigate the reliability of our brain shift compensation method based on anatomic landmarks. Additionally, the minimum distance between the DESS points and warped OT (wOT) model was measured and correlated with DESS intensity threshold. Brain shift compensation was achieved in all cases, and the area under the ROC curve was 0.96 in the registration accuracy analysis. The minimum distance between the DESS points and the wOT model was found to have a significantly high correlation with the DESS stimulation intensity threshold (r = 0.87, P < 0.001), with a linear regression coefficient of 0.96. Our OT method can provide comprehensive and accurate visualization of the pyramidal tracts for neurosurgical navigation and was quantitatively verified by intraoperative DESS after brain shift compensation.

The organization of frontostriatal brain wiring in non-affective early psychosis compared with healthy subjects using a novel diffusion imaging fiber cluster analysis

BACKGROUND

Alterations in brain connectivity may underlie neuropsychiatric conditions such as schizophrenia. We here assessed the degree of convergence of frontostriatal fiber projections in 56 young adult healthy controls (HCs) and 108 matched Early Psychosis-Non-Affective patients (EP-NAs) using our novel fiber cluster analysis of whole brain diffusion magnetic resonance imaging tractography.

METHODS

Using whole brain tractography and our fiber clustering methodology on harmonized diffusion magnetic resonance imaging data from the Human Connectome Project for Early Psychosis we identified 17 white matter fiber clusters that connect frontal cortex (FCtx) and caudate (Cd) per hemisphere in each group. To quantify the degree of convergence and, hence, topographical relationship of these fiber clusters, we measured the inter-cluster mean distances between the endpoints of the fiber clusters at the level of the FCtx and of the Cd, respectively.

RESULTS

We found (1) in both groups, bilaterally, a non-linear relationship, yielding convex curves, between FCtx and Cd distances for FCtx-Cd connecting fiber clusters, driven by a cluster projecting from inferior frontal gyrus; however, in the right hemisphere, the convex curve was more flattened in EP-NAs; (2) that cluster pairs in the right (p = 0.03), but not left (p = 0.13), hemisphere were significantly more convergent in HCs vs EP-NAs; (3) in both groups, bilaterally, similar clusters projected significantly convergently to the Cd; and, (4) a significant group by fiber cluster pair interaction for 2 right hemisphere fiber clusters (numbers 5, 11; p = .00023; p = .00023) originating in selective PFC subregions.

CONCLUSIONS

In both groups, we found the FCtx-Cd wiring pattern deviated from a strictly topographic relationship and that similar clusters projected significantly more convergently to the Cd. Interestingly, we also found a significantly more convergent pattern of connectivity in HCs in the right hemisphere and that 2 clusters from PFC subregions in the right hemisphere significantly differed in their pattern of connectivity between groups.

Directionally encoded color track density imaging in brain tumor patients: A potential application to neuro-oncology surgical planning

BACKGROUND

Diffusion magnetic resonance imaging white matter tractography, an increasingly popular preoperative planning modality used for pre-surgical planning in brain tumor patients, is employed with the goal of maximizing tumor resection while sparing postoperative neurological function. Clinical translation of white matter tractography has been limited by several shortcomings of standard diffusion tensor imaging (DTI), including poor modeling of fibers crossing through regions of peritumoral edema and low spatial resolution for typical clinical diffusion MRI (dMRI) sequences. Track density imaging (TDI) is a post-tractography technique that uses the number of tractography streamlines and their long-range continuity to map the white matter connections of the brain with enhanced image resolution relative to the acquired dMRI data, potentially offering improved white matter visualization in patients with brain tumors. The aim of this study was to assess the utility of TDI-based white matter maps in a neurosurgical planning context compared to the current clinical standard of DTI-based white matter maps.

METHODS

Fourteen consecutive brain tumor patients from a single institution were retrospectively selected for the study. Each patient underwent 3-Tesla dMRI scanning with 30 gradient directions and a b-value of 1000 s/mm2. For each patient, two directionally encoded color (DEC) maps were produced as follows. DTI-based DEC-fractional anisotropy maps (DEC-FA) were generated on the scanner, while DEC-track density images (DEC-TDI) were generated using constrained spherical deconvolution based tractography. The potential clinical utility of each map was assessed by five practicing neurosurgeons, who rated the maps according to four clinical utility statements regarding different clinical aspects of pre-surgical planning. The neurosurgeons rated each map according to their agreement with four clinical utility statements regarding if the map 1 identified clinically relevant tracts, (2) helped establish a goal resection margin, (3) influenced a planned surgical route, and (4) was useful overall. Cumulative link mixed effect modeling and analysis of variance were performed to test the primary effect of map type (DEC-TDI vs. DEC-FA) on rater score. Pairwise comparisons using estimated marginal means were then calculated to determine the magnitude and directionality of differences in rater scores by map type.

RESULTS

A majority of rater responses agreed with the four clinical utility statements, indicating that neurosurgeons found both DEC maps to be useful. Across all four investigated clinical utility statements, the DEC map type significantly influenced rater score. Rater scores were significantly higher for DEC-TDI maps compared to DEC-FA maps. The largest effect size in rater scores in favor of DEC-TDI maps was observed for clinical utility statement 2, which assessed establishing a goal resection margin.

CONCLUSION

We observed a significant neurosurgeon preference for DEC-TDI maps, indicating their potential utility for neurosurgical planning.

Spatially regularized low-rank tensor approximation for accurate and fast tractography

Tractography based on diffusion Magnetic Resonance Imaging (dMRI) is the prevalent approach to the in vivo delineation of white matter tracts in the human brain. Many tractography methods rely on models of multiple fiber compartments, but the local dMRI information is not always sufficient to reliably estimate the directions of secondary fibers. Therefore, we introduce two novel approaches that use spatial regularization to make multi-fiber tractography more stable. Both represent the fiber Orientation Distribution Function (fODF) as a symmetric fourth-order tensor, and recover multiple fiber orientations via low-rank approximation. Our first approach computes a joint approximation over suitably weighted local neighborhoods with an efficient alternating optimization. The second approach integrates the low-rank approximation into a current state-of-the-art tractography algorithm based on the unscented Kalman filter (UKF). These methods were applied in three different scenarios. First, we demonstrate that they improve tractography even in high-quality data from the Human Connectome Project, and that they maintain useful results with a small fraction of the measurements. Second, on the 2015 ISMRM tractography challenge, they increase overlap, while reducing overreach, compared to low-rank approximation without joint optimization or the traditional UKF, respectively. Finally, our methods permit a more comprehensive reconstruction of tracts surrounding a tumor in a clinical dataset. Overall, both approaches improve reconstruction quality. At the same time, our modified UKF significantly reduces the computational effort compared to its traditional counterpart, and to our joint approximation. However, when used with ROI-based seeding, joint approximation more fully recovers fiber spread.

Quantitative mapping of the brain's structural connectivity using diffusion MRI tractography: A review

Diffusion magnetic resonance imaging (dMRI) tractography is an advanced imaging technique that enables in vivo reconstruction of the brain's white matter connections at macro scale. It provides an important tool for quantitative mapping of the brain's structural connectivity using measures of connectivity or tissue microstructure. Over the last two decades, the study of brain connectivity using dMRI tractography has played a prominent role in the neuroimaging research landscape. In this paper, we provide a high-level overview of how tractography is used to enable quantitative analysis of the brain's structural connectivity in health and disease. We focus on two types of quantitative analyses of tractography, including: 1) tract-specific analysis that refers to research that is typically hypothesis-driven and studies particular anatomical fiber tracts, and 2) connectome-based analysis that refers to research that is more data-driven and generally studies the structural connectivity of the entire brain. We first provide a review of methodology involved in three main processing steps that are common across most approaches for quantitative analysis of tractography, including methods for tractography correction, segmentation and quantification. For each step, we aim to describe methodological choices, their popularity, and potential pros and cons. We then review studies that have used quantitative tractography approaches to study the brain's white matter, focusing on applications in neurodevelopment, aging, neurological disorders, mental disorders, and neurosurgery. We conclude that, while there have been considerable advancements in methodological technologies and breadth of applications, there nevertheless remains no consensus about the "best" methodology in quantitative analysis of tractography, and researchers should remain cautious when interpreting results in research and clinical applications.

Effects of Multi-Shell Free Water Correction on Glioma Characterization

Diffusion MRI is a useful tool to investigate the microstructure of brain tumors. However, the presence of fast diffusing isotropic signals originating from non-restricted edematous fluids, within and surrounding tumors, may obscure estimation of the underlying tissue characteristics, complicating the radiological interpretation and quantitative evaluation of diffusion MRI. A multi-shell regularized free water (FW) elimination model was therefore applied to separate free water from tissue-related diffusion components from the diffusion MRI of 26 treatment-naïve glioma patients. We then investigated the diagnostic value of the derived measures of FW maps as well as FW-corrected tensor-derived maps of fractional anisotropy (FA). Presumed necrotic tumor regions display greater mean and variance of FW content than other parts of the tumor. On average, the area under the receiver operating characteristic (ROC) for the classification of necrotic and enhancing tumor volumes increased by 5% in corrected data compared to non-corrected data. FW elimination shifts the FA distribution in non-enhancing tumor parts toward higher values and significantly increases its entropy (p ≤ 0.003), whereas skewness is decreased (p ≤ 0.004). Kurtosis is significantly decreased (p < 0.001) in high-grade tumors. In conclusion, eliminating FW contributions improved quantitative estimations of FA, which helps to disentangle the cancer heterogeneity.

Clinical applications of magnetic resonance imaging based functional and structural connectivity

Advances in computational neuroimaging techniques have expanded the armamentarium of imaging tools available for clinical applications in clinical neuroscience. Non-invasive, in vivo brain MRI structural and functional network mapping has been used to identify therapeutic targets, define eloquent brain regions to preserve, and gain insight into pathological processes and treatments as well as prognostic biomarkers. These tools have the real potential to inform patient-specific treatment strategies. Nevertheless, a realistic appraisal of clinical utility is needed that balances the growing excitement and interest in the field with important limitations associated with these techniques. Quality of the raw data, minutiae of the processing methodology, and the statistical models applied can all impact on the results and their interpretation. A lack of standardization in data acquisition and processing has also resulted in issues with reproducibility. This limitation has had a direct impact on the reliability of these tools and ultimately, confidence in their clinical use. Advances in MRI technology and computational power as well as automation and standardization of processing methods, including machine learning approaches, may help address some of these issues and make these tools more reliable in clinical use. In this review, we will highlight the current clinical uses of MRI connectomics in the diagnosis and treatment of neurological disorders; balancing emerging applications and technologies with limitations of connectivity analytic approaches to present an encompassing and appropriate perspective.

Exploring Peritumoral Neural Tracts by Using Neurite Orientation Dispersion and Density Imaging

Diffusion Tensor Imaging (DTI) tractography has been widely used in brain tumor surgery to ensure thorough resection and minimize functional damage. However, due to enhanced anisotropic uncertainty in the area with peritumoral edema, diffusion tractography is generally not practicable leading to high false-negative results in neural tracking. In this study, we evaluated the usefulness of the neurite orientation dispersion and density imaging (NODDI) derived tractography for investigating structural heterogeneity of the brain in patients with brain tumor. A total of 24 patients with brain tumors, characterized by peritumoral edema, and 10 healthy counterparts were recruited from 2014 to 2021. All participants underwent magnetic resonance imaging. Moreover, we used the images obtained from the healthy participants for calibrating the orientation dispersion threshold for NODDI-derived corticospinal tract (CST) reconstruction. Compared to DTI, NODDI-derived tractography has a great potential to improve the reconstruction of fiber tracking through regions of vasogenic edema. The regions with edematous CST in NODDI-derived tractography demonstrated a significant decrease in the intracellular volume fraction (VF_ic, p < 0.000) and an increase in the isotropic volume fraction (VF_iso, p < 0.014). Notably, the percentage of the involved volume of the concealed CST and lesion-to-tract distance could reflect the motor function of the patients. After the tumor resection, four patients with 1–5 years follow-up were showed subsidence of the vasogenic edema and normal CST on DTI tractography. NODDI-derived tractography revealed tracts within the edematous area and could assist neurosurgeons to locate the neural tracts that are otherwise not visualized by conventional DTI tractography.

Comparison of multiple tractography methods for reconstruction of the retinogeniculate visual pathway using diffusion MRI

The retinogeniculate visual pathway (RGVP) conveys visual information from the retina to the lateral geniculate nucleus. The RGVP has four subdivisions, including two decussating and two nondecussating pathways that cannot be identified on conventional structural magnetic resonance imaging (MRI). Diffusion MRI tractography has the potential to trace these subdivisions and is increasingly used to study the RGVP. However, it is not yet known which fiber tracking strategy is most suitable for RGVP reconstruction. In this study, four tractography methods are compared, including constrained spherical deconvolution (CSD) based probabilistic (iFOD1) and deterministic (SD‐Stream) methods, and multi‐fiber (UKF‐2T) and single‐fiber (UKF‐1T) unscented Kalman filter (UKF) methods. Experiments use diffusion MRI data from 57 subjects in the Human Connectome Project. The RGVP is identified using regions of interest created by two clinical experts. Quantitative anatomical measurements and expert anatomical judgment are used to assess the advantages and limitations of the four tractography methods. Overall, we conclude that UKF‐2T and iFOD1 produce the best RGVP reconstruction results. The iFOD1 method can better quantitatively estimate the percentage of decussating fibers, while the UKF‐2T method produces reconstructed RGVPs that are judged to better correspond to the known anatomy and have the highest spatial overlap across subjects. Overall, we find that it is challenging for current tractography methods to both accurately track RGVP fibers that correspond to known anatomy and produce an approximately correct percentage of decussating fibers. We suggest that future algorithm development for RGVP tractography should take consideration of both of these two points.

Spatial correspondence of spinal cord white matter tracts using diffusion tensor imaging, fibre tractography, and atlas-based segmentation

PURPOSE

Neuroimaging provides great utility in complex spinal surgeries, particularly when anatomical geometry is distorted by pathology (tumour, degeneration, etc.). Spinal cord MRI diffusion tractography can be used to generate streamlines; however, it is unclear how well they correspond with white matter tract locations along the cord microstructure. The goal of this work was to evaluate the spatial correspondence of DTI tractography with anatomical MRI in healthy anatomy (where anatomical locations can be well defined in T1-weighted images).

METHODS

Ten healthy volunteers were scanned on a 3T system. T1-weighted (1 × 1 × 1 mm) and diffusion-weighted images (EPI readout, 2 × 2 × 2 mm, 30 gradient directions) were acquired and subsequently registered (Spinal Cord Toolbox (SCT)). Atlas-based (SCT) anatomic label maps of the left and right lateral corticospinal tracts were identified for each vertebral region (C2-C6) from T1 images. Tractography streamlines were generated with a customized approach, enabling seeding of specific spinal tract regions corresponding to individual vertebral levels. Spatial correspondence of generated fibre streamlines with anatomic tract segmentations was compared in unseeded regions of interest (ROIs).

RESULTS

Spatial correspondence of the lateral corticospinal tract streamlines was good over a single vertebral ROI (Dice's similarity coefficient (DSC) = 0.75 ± 0.08, Hausdorff distance = 1.08 ± 0.17 mm). Over larger ROI, fair agreement between tractography and anatomical labels was achieved (two levels: DSC = 0.67 ± 0.13, three levels: DSC = 0.52 ± 0.19).

CONCLUSION

DTI tractography produced good spatial correspondence with anatomic white matter tracts, superior to the agreement between multiple manual tract segmentations (DSC ~ 0.5). This supports further development of spinal cord tractography for computer-assisted neurosurgery.

Creation of a novel trigeminal tractography atlas for automated trigeminal nerve identification

Diffusion MRI (dMRI) tractography has been successfully used to study the trigeminal nerves (TGNs) in many clinical and research applications. Currently, identification of the TGN in tractography data requires expert nerve selection using manually drawn regions of interest (ROIs), which is prone to inter-observer variability, time-consuming and carries high clinical and labor costs. To overcome these issues, we propose to create a novel anatomically curated TGN tractography atlas that enables automated identification of the TGN from dMRI tractography. In this paper, we first illustrate the creation of a trigeminal tractography atlas. Leveraging a well-established computational pipeline and expert neuroanatomical knowledge, we generate a data-driven TGN fiber clustering atlas using tractography data from 50 subjects from the Human Connectome Project. Then, we demonstrate the application of the proposed atlas for automated TGN identification in new subjects, without relying on expert ROI placement. Quantitative and visual experiments are performed with comparison to expert TGN identification using dMRI data from two different acquisition sites. We show highly comparable results between the automatically and manually identified TGNs in terms of spatial overlap and visualization, while our proposed method has several advantages. First, our method performs automated TGN identification, and thus it provides an efficient tool to reduce expert labor costs and inter-operator bias relative to expert manual selection. Second, our method is robust to potential imaging artifacts and/or noise that can prevent successful manual ROI placement for TGN selection and hence yields a higher successful TGN identification rate.

Investigation of local white matter abnormality in Parkinson's disease by using an automatic fiber tract parcellation

The deficits of white matter (WM) microstructure are involved during Parkinson's disease (PD) progression. Most current methods identify key WM tracts relying on cortical regions of interest (ROIs). However, such ROI methods can be challenged due to low diffusion anisotropy near the gray matter (GM), which could result in a low sensitivity of tract identification. This work proposes an automatic WM parcellation method to improve the accuracy of WM tract identification and locate abnormal tracts by using sensitive features. The proposed method consists of 1) whole brain WM parcellation using an established fiber clustering method, without using any ROIs, 2) features of fasciculus were calculated to quantify diffusion measures at each equal cross-section along the whole cluster. Then, we use the proposed features to investigate the WM difference in PD compared with healthy controls (HC). We also use these features to investigate the relationship of clinical symptoms and specific fiber tracts. The novelty of the proposed method is that it automatically identifies the abnormal WM fibers in cluster degree. Experiment results indicated that the proposed method had advantage in detecting the local WM abnormality by performing between-group statistical analysis in 30 patients with PD and 28 HC. We found 13 hemisphere clusters and 8 commissural clusters had significant group difference (p < 0.05, corrected by FDR method) in local regions, which belonged to multiple fiber tracts including cingulum bundle (CB), inferior occipito-frontal fasciculus (IoFF), corpus callosum (CC), external capsule (EC), uncinate fasciculus (UF), superior longitudinal fasciculus (SLF) and thalamo front (TF). We also found clusters that had relevance with clinical indices of cognitive function (2 clusters), athletic function (6 clusters), and depressive state (2 clusters) in these significant clusters. From the experiment results, it confirmed the ability of the proposed method to identify potential WM microstructure abnormality.

Individual variations of the human corticospinal tract and its hand-related motor fibers using diffusion MRI tractography

The corticospinal tract (CST) is one of the most well studied tracts in human neuroanatomy. Its clinical significance can be demonstrated in many notable traumatic conditions and diseases such as stroke, spinal cord injury (SCI) or amyotrophic lateral sclerosis (ALS). With the advent of diffusion MRI and tractography the computational representation of the human CST in a 3D model became available. However, the representation of the entire CST and, specifically, the hand motor area has remained elusive. In this paper we propose a novel method, using manually drawn ROIs based on robustly identifiable neuroanatomic structures to delineate the entire CST and isolate its hand motor representation as well as to estimate their variability and generate a database of their volume, length and biophysical parameters. Using 37 healthy human subjects we performed a qualitative and quantitative analysis of the CST and the hand-related motor fiber tracts (HMFTs). Finally, we have created variability heat maps from 37 subjects for both the aforementioned tracts, which could be utilized as a reference for future studies with clinical focus to explore neuropathology in both trauma and disease states.

Implementation of clinical tractography for pre-surgical planning of space occupying lesions: An investigation of common acquisition and post-processing methods compared to dissection studies

Background and purpose

There is limited standardization of acquisition and processing methods in diffusion tractography for pre-surgical planning, leading to a range of approaches. In this study, a number of representative acquisition variants and post processing methods are considered, to assess their importance when implementing a clinical tractography program.

Methods

Diffusion MRI was undertaken in ten healthy volunteers, using protocols typical of clinical and research scanning: a 32-direction diffusion acquisition with and without peripheral gating, and a non-gated 64 diffusion direction acquisition. All datasets were post-processed using diffusion tensor reconstruction with streamline tractography, and with constrained spherical deconvolution (CSD) with both streamline and probabilistic tractography, to delineate the cortico-spinal tract (CST) and optic radiation (OR). The accuracy of tractography results was assessed against a histological atlas using a novel probabilistic Dice overlap technique, together with direct comparison to tract volumes and distance of Meyer’s loop to temporal pole (ML-TP) from dissections studies. Three clinical case studies of patients with space occupying lesions were also investigated.

Results

Tracts produced by CSD with probabilistic tractography provided the greatest overlap with the histological atlas (overlap scores of 44% and 52% for the CST and OR, respectively) and best matched tract volume and ML-TP distance from dissection studies. The acquisition protocols investigated had limited impact on the accuracy of the tractography. In all patients, the CSD based probabilistic tractography created tracts with greatest anatomical plausibility, although in one case anatomically plausible pathways could not be reconstructed without reducing the probabilistic threshold, leading to an increase in false positive tracts.

Conclusions

Advanced post processing techniques such as CSD with probabilistic tractography are vital for pre-surgical planning. However, overall accuracy relative to dissection studies remains limited.

Tractography in Neurosurgery: A Systematic Review of Current Applications

The ability to visualize the brain's fiber connections noninvasively in vivo is relatively young compared with other possibilities of functional magnetic resonance imaging. Although many studies showed tractography to be of promising value for neurosurgical care, the implications remain inconclusive. An overview of current applications is presented in this systematic review. A search was conducted for (("tractography" or "fiber tracking" or "fibre tracking") and "neurosurgery") that produced 751 results. We identified 260 relevant articles and added 20 more from other sources. Most publications concerned surgical planning for resection of tumors (n = 193) and vascular lesions (n = 15). Preoperative use of transcranial magnetic stimulation was discussed in 22 of these articles. Tractography in skull base surgery presents a special challenge (n = 29). Fewer publications evaluated traumatic brain injury (TBI) (n = 25) and spontaneous intracranial bleeding (n = 22). Twenty-three articles focused on tractography in pediatric neurosurgery. Most authors found tractography to be a valuable addition in neurosurgical care. The accuracy of the technique has increased over time. There are articles suggesting that tractography improves patient outcome after tumor resection. However, no reliable biomarkers have yet been described. The better rehabilitation potential after TBI and spontaneous intracranial bleeding compared with brain tumors offers an insight into the process of neurorehabilitation. Tractography and diffusion measurements in some studies showed a correlation with patient outcome that might help uncover the neuroanatomical principles of rehabilitation itself. Alternative corticofugal and cortico-cortical networks have been implicated in motor recovery after ischemic stroke, suggesting more complex mechanisms in neurorehabilitation that go beyond current models. Hence tractography may potentially be able to predict clinical deficits and rehabilitation potential, as well as finding possible explanations for neurologic disorders in retrospect. However, large variations of the results indicate a lack of data to establish robust diagnostical concepts at this point. Therefore, in vivo tractography should still be interpreted with caution and by experienced surgeons.

SlicerDMRI: Diffusion MRI and Tractography Research Software for Brain Cancer Surgery Planning and Visualization

PURPOSE

We present SlicerDMRI, an open-source software suite that enables research using diffusion magnetic resonance imaging (dMRI), the only modality that can map the white matter connections of the living human brain. SlicerDMRI enables analysis and visualization of dMRI data and is aimed at the needs of clinical research users. SlicerDMRI is built upon and deeply integrated with 3D Slicer, a National Institutes of Health-supported open-source platform for medical image informatics, image processing, and three-dimensional visualization. Integration with 3D Slicer provides many features of interest to cancer researchers, such as real-time integration with neuronavigation equipment, intraoperative imaging modalities, and multimodal data fusion. One key application of SlicerDMRI is in neurosurgery research, where brain mapping using dMRI can provide patient-specific maps of critical brain connections as well as insight into the tissue microstructure that surrounds brain tumors.

PATIENTS AND METHODS

In this article, we focus on a demonstration of SlicerDMRI as an informatics tool to enable end-to-end dMRI analyses in two retrospective imaging data sets from patients with high-grade glioma. Analyses demonstrated here include conventional diffusion tensor analysis, advanced multifiber tractography, automated identification of critical fiber tracts, and integration of multimodal imagery with dMRI.

RESULTS

We illustrate the ability of SlicerDMRI to perform both conventional and advanced dMRI analyses as well as to enable multimodal image analysis and visualization. We provide an overview of the clinical rationale for each analysis along with pointers to the SlicerDMRI tools used in each.

CONCLUSION

SlicerDMRI provides open-source and clinician-accessible research software tools for dMRI analysis. SlicerDMRI is available for easy automated installation through the 3D Slicer Extension Manager.

Anatomical assessment of trigeminal nerve tractography using diffusion MRI: A comparison of acquisition b-values and single- and multi-fiber tracking strategies

•

Investigation of the performance of multiple dMRI acquisitions and fiber models for trigeminal nerve (TGN) identification.

•

Expert rating study of over 1000 TGN visualizations using seven proposed expert rating anatomical criteria.

•

The two-tensor tractography method had better performance on identifying true positive structures, while generating more false positive streamlines in comparison to the single-tensor tractography method.

•

TGN tracking performance was significantly different across the three b-values for almost all structures studied.

Background

The trigeminal nerve (TGN) is the largest cranial nerve and can be involved in multiple inflammatory, compressive, ischemic or other pathologies. Currently, imaging-based approaches to identify the TGN mostly rely on T2-weighted magnetic resonance imaging (MRI), which provides localization of the cisternal portion of the TGN where the contrast between nerve and cerebrospinal fluid (CSF) is high enough to allow differentiation. The course of the TGN within the brainstem as well as anterior to the cisternal portion, however, is more difficult to display on traditional imaging sequences. An advanced imaging technique, diffusion MRI (dMRI), enables tracking of the trajectory of TGN fibers and has the potential to visualize anatomical regions of the TGN not seen on T2-weighted imaging. This may allow a more comprehensive assessment of the nerve in the context of pathology. To date, most work in TGN tracking has used clinical dMRI acquisitions with a b-value of 1000 s/mm² and conventional diffusion tensor MRI (DTI) tractography methods. Though higher b-value acquisitions and multi-tensor tractography methods are known to be beneficial for tracking brain white matter fiber tracts, there have been no studies conducted to evaluate the performance of these advanced approaches on nerve tracking of the TGN, in particular on tracking different anatomical regions of the TGN.

Objective

We compare TGN tracking performance using dMRI data with different b-values, in combination with both single- and multi-tensor tractography methods. Our goal is to assess the advantages and limitations of these different strategies for identifying the anatomical regions of the TGN.

Methods

We proposed seven anatomical rating criteria including true and false positive structures, and we performed an expert rating study of over 1000 TGN visualizations, as follows. We tracked the TGN using high-quality dMRI data from 100 healthy adult subjects from the Human Connectome Project (HCP). TGN tracking performance was compared across dMRI acquisitions with b = 1000 s/mm², b = 2000 s/mm² and b = 3000 s/mm², using single-tensor (1T) and two-tensor (2T) unscented Kalman filter (UKF) tractography. This resulted in a total of six tracking strategies. The TGN was identified using an anatomical region-of-interest (ROI) selection approach. First, in a subset of the dataset we identified ROIs that provided good TGN tracking performance across all tracking strategies. Using these ROIs, the TGN was then tracked in all subjects using the six tracking strategies. An expert rater (GX) visually assessed and scored each TGN based on seven anatomical judgment criteria. These criteria included the presence of multiple expected anatomical segments of the TGN (true positive structures), specifically branch-like structures, cisternal portion, mesencephalic trigeminal tract, and spinal cord tract of the TGN. False positive criteria included the presence of any fibers entering the temporal lobe, the inferior cerebellar peduncle, or the middle cerebellar peduncle. Expert rating scores were analyzed to compare TGN tracking performance across the six tracking strategies. Intra- and inter-rater validation was performed to assess the reliability of the expert TGN rating result.

Results

The TGN was selected using two anatomical ROIs (Meckel's Cave and cisternal portion of the TGN). The two-tensor tractography method had significantly better performance on identifying true positive structures, while generating more false positive streamlines in comparison to the single-tensor tractography method. TGN tracking performance was significantly different across the three b-values for almost all structures studied. Tracking performance was reported in terms of the percentage of subjects achieving each anatomical rating criterion. Tracking of the cisternal portion and branching structure of the TGN was generally successful, with the highest performance of over 98% using two-tensor tractography and b = 1000 or b = 2000. However, tracking the smaller mesencephalic and spinal cord tracts of the TGN was quite challenging (highest performance of 37.5% and 57.07%, using two-tensor tractography with b = 1000 and b = 2000, respectively). False positive connections to the temporal lobe (over 38% of subjects for all strategies) and cerebellar peduncles (100% of subjects for all strategies) were prevalent. High joint probability of agreement was obtained in the inter-rater (on average 83%) and intra-rater validation (on average 90%), showing a highly reliable expert rating result.

Conclusions

Overall, the results of the study suggest that researchers and clinicians may benefit from tailoring their acquisition and tracking methodology to the specific anatomical portion of the TGN that is of the greatest interest. For example, tracking of branching structures and TGN-T2 overlap can be best achieved with a two-tensor model and an acquisition using b = 1000 or b = 2000. In general, b = 1000 and b = 2000 acquisitions provided the best-rated tracking results. Further research is needed to improve both sensitivity and specificity of the depiction of the TGN anatomy using dMRI.

Test-retest reproducibility of white matter parcellation using diffusion MRI tractography fiber clustering

There are two popular approaches for automated white matter parcellation using diffusion MRI tractography, including fiber clustering strategies that group white matter fibers according to their geometric trajectories and cortical-parcellation-based strategies that focus on the structural connectivity among different brain regions of interest. While multiple studies have assessed test-retest reproducibility of automated white matter parcellations using cortical-parcellation-based strategies, there are no existing studies of test-retest reproducibility of fiber clustering parcellation. In this work, we perform what we believe is the first study of fiber clustering white matter parcellation test-retest reproducibility. The assessment is performed on three test-retest diffusion MRI datasets including a total of 255 subjects across genders, a broad age range (5-82 years), health conditions (autism, Parkinson's disease and healthy subjects), and imaging acquisition protocols (three different sites). A comprehensive evaluation is conducted for a fiber clustering method that leverages an anatomically curated fiber clustering white matter atlas, with comparison to a popular cortical-parcellation-based method. The two methods are compared for the two main white matter parcellation applications of dividing the entire white matter into parcels (i.e., whole brain white matter parcellation) and identifying particular anatomical fiber tracts (i.e., anatomical fiber tract parcellation). Test-retest reproducibility is measured using both geometric and diffusion features, including volumetric overlap (wDice) and relative difference of fractional anisotropy. Our experimental results in general indicate that the fiber clustering method produced more reproducible white matter parcellations than the cortical-parcellation-based method.

Advances in computational and statistical diffusion MRI

Computational methods are crucial for the analysis of diffusion magnetic resonance imaging (MRI) of the brain. Computational diffusion MRI can provide rich information at many size scales, including local microstructure measures such as diffusion anisotropies or apparent axon diameters, whole-brain connectivity information that describes the brain's wiring diagram and population-based studies in health and disease. Many of the diffusion MRI analyses performed today were not possible five, ten or twenty years ago, due to the requirements for large amounts of computer memory or processor time. In addition, mathematical frameworks had to be developed or adapted from other fields to create new ways to analyze diffusion MRI data. The purpose of this review is to highlight recent computational and statistical advances in diffusion MRI and to put these advances into context by comparison with the more traditional computational methods that are in popular clinical and scientific use. We aim to provide a high-level overview of interest to diffusion MRI researchers, with a more in-depth treatment to illustrate selected computational advances.

Investigation into local white matter abnormality in emotional processing and sensorimotor areas using an automatically annotated fiber clustering in major depressive disorder

This work presents an automatically annotated fiber cluster (AAFC) method to enable identification of anatomically meaningful white matter structures from the whole brain tractography. The proposed method consists of 1) a study-specific whole brain white matter parcellation using a well-established data-driven groupwise fiber clustering pipeline to segment tractography into multiple fiber clusters, and 2) a novel cluster annotation method to automatically assign an anatomical tract annotation to each fiber cluster by employing cortical parcellation information across multiple subjects. The novelty of the AAFC method is that it leverages group-wise information about the fiber clusters, including their fiber geometry and cortical terminations, to compute a tract anatomical label for each cluster in an automated fashion. We demonstrate the proposed AAFC method in an application of investigating white matter abnormality in emotional processing and sensorimotor areas in major depressive disorder (MDD). Seven tracts of interest related to emotional processing and sensorimotor functions are automatically identified using the proposed AAFC method as well as a comparable method that uses a cortical parcellation alone. Experimental results indicate that our proposed method is more consistent in identifying the tracts across subjects and across hemispheres in terms of the number of fibers. In addition, we perform a between-group statistical analysis in 31 MDD patients and 62 healthy subjects on the identified tracts using our AAFC method. We find statistical differences in diffusion measures in local regions within a fiber tract (e.g. 4 fiber clusters within the identified left hemisphere cingulum bundle (consisting of 14 clusters) are significantly different between the two groups), suggesting the ability of our method in identifying potential abnormality specific to subdivisions of a white matter structure.

An anatomically curated fiber clustering white matter atlas for consistent white matter tract parcellation across the lifespan

This work presents an anatomically curated white matter atlas to enable consistent white matter tract parcellation across different populations. Leveraging a well-established computational pipeline for fiber clustering, we create a tract-based white matter atlas including information from 100 subjects. A novel anatomical annotation method is proposed that leverages population-based brain anatomical information and expert neuroanatomical knowledge to annotate and categorize the fiber clusters. A total of 256 white matter structures are annotated in the proposed atlas, which provides one of the most comprehensive tract-based white matter atlases covering the entire brain to date. These structures are composed of 58 deep white matter tracts including major long range association and projection tracts, commissural tracts, and tracts related to the brainstem and cerebellar connections, plus 198 short and medium range superficial fiber clusters organized into 16 categories according to the brain lobes they connect. Potential false positive connections are annotated in the atlas to enable their exclusion from analysis or visualization. In addition, the proposed atlas allows for a whole brain white matter parcellation into 800 fiber clusters to enable whole brain connectivity analyses. The atlas and related computational tools are open-source and publicly available. We evaluate the proposed atlas using a testing dataset of 584 diffusion MRI scans from multiple independently acquired populations, across genders, the lifespan (1 day-82 years), and different health conditions (healthy control, neuropsychiatric disorders, and brain tumor patients). Experimental results show successful white matter parcellation across subjects from different populations acquired on multiple scanners, irrespective of age, gender or disease indications. Over 99% of the fiber tracts annotated in the atlas were detected in all subjects on average. One advantage in terms of robustness is that the tract-based pipeline does not require any cortical or subcortical segmentations, which can have limited success in young children and patients with brain tumors or other structural lesions. We believe this is the first demonstration of consistent automated white matter tract parcellation across the full lifespan from birth to advanced age.

Anatomical accuracy of standard-practice tractography algorithms in the motor system - A histological validation in the squirrel monkey brain

For two decades diffusion fiber tractography has been used to probe both the spatial extent of white matter pathways and the region to region connectivity of the brain. In both cases, anatomical accuracy of tractography is critical for sound scientific conclusions. Here we assess and validate the algorithms and tractography implementations that have been most widely used - often because of ease of use, algorithm simplicity, or availability offered in open source software. Comparing forty tractography results to a ground truth defined by histological tracers in the primary motor cortex on the same squirrel monkey brains, we assess tract fidelity on the scale of voxels as well as over larger spatial domains or regional connectivity. No algorithms are successful in all metrics, and, in fact, some implementations fail to reconstruct large portions of pathways or identify major points of connectivity. The accuracy is most dependent on reconstruction method and tracking algorithm, as well as the seed region and how this region is utilized. We also note a tremendous variability in the results, even though the same MR images act as inputs to all algorithms. In addition, anatomical accuracy is significantly decreased at increased distances from the seed. An analysis of the spatial errors in tractography reveals that many techniques have trouble properly leaving the gray matter, and many only reveal connectivity to adjacent regions of interest. These results show that the most commonly implemented algorithms have several shortcomings and limitations, and choices in implementations lead to very different results. This study should provide guidance for algorithm choices based on study requirements for sensitivity, specificity, or the need to identify particular connections, and should serve as a heuristic for future developments in tractography.

Differentiation of Edematous, Tumoral and Normal Areas of Brain Using Diffusion Tensor and Neurite Orientation Dispersion and Density Imaging

BACKGROUND

Presurigical planning for glioma tumor resection and radiotherapy treatment require proper delineation of tumoral and peritumoral areas of brain. Diffusion tensor imaging (DTI) is the most common mathematical model applied for diffusion weighted MRI data. Neurite orientation dispersion and density imaging (NODDI) is another mathematical model for DWI data modeling.

OBJECTIVE

We studied whether extracted parameters of DTI, and NODDI models can be used to differentiate between edematous, tumoral, and normal areas in brain white matter (WM).

MATERIAL AND METHODS

12 patients with peritumoral edema underwent 3T multi-shell diffusion imaging with b-values of 1000 and 2000 smm-2 in 30 and 64 gradient directions, respectively. We fitted DTI and NODDI to data in manually drawn regions of interest and used their derived parameters to characterize edematous, tumoral and normal brain areas.

RESULTS

We found that DTI parameters fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD) and radial diffusivity (RD) all significantly differentiated edematous from contralateral normal brain WM (p<0.005). However, only FA was found to distinguish between edematous WM fibers and tumor invaded fibers (p = 0.001). Among NODDI parameters, the intracellular volume fraction (ficvf) had the best distinguishing power with (p = 0.001) compared with the isotropic volume fraction (fiso), the orientation dispersion index (odi), and the concentration parameter of Watson distribution (κ), while comparing fibers inside normal, tumoral, and edematous areas.

CONCLUSION

The combination of two diffusion based methods, i.e. DTI and NODDI parameters can distinguish and characterize WM fibers involved in edematus, tumoral, and normal brain areas with reasonable confidence. Further studies will be required to improve the detectability of WM fibers inside the solid tumor if they hypothetically exist in tumoral parenchyma.

Free water modeling of peritumoral edema using multi-fiber tractography: Application to tracking the arcuate fasciculus for neurosurgical planning

Purpose

Peritumoral edema impedes the full delineation of fiber tracts due to partial volume effects in image voxels that contain a mixture of cerebral parenchyma and extracellular water. The purpose of this study is to investigate the effect of incorporating a free water (FW) model of edema for white matter tractography in the presence of edema.

Materials and methods

We retrospectively evaluated 26 consecutive brain tumor patients with diffusion MRI and T2-weighted images acquired presurgically. Tractography of the arcuate fasciculus (AF) was performed using the two-tensor unscented Kalman filter tractography (UKFt) method, the UKFt method with a reduced fiber tracking stopping fractional anisotropy (FA) threshold (UKFt+rFA), and the UKFt method with the addition of a FW compartment (UKFt+FW). An automated white matter fiber tract identification approach was applied to delineate the AF. Quantitative measurements included tract volume, edema volume, and mean FW fraction. Visual comparisons were performed by three experts to evaluate the quality of the detected AF tracts.

Results

The AF volume in edematous brain hemispheres was significantly larger using the UKFt+FW method (p<0.0001) compared to UKFt, but not significantly larger (p = 0.0996) in hemispheres without edema. The AF size increase depended on the volume of edema: a significant correlation was found between AF volume affected by (intersecting) edema and AF volume change with the FW model (Pearson r = 0.806, p<0.0001). The mean FW fraction was significantly larger in tracts intersecting edema (p = 0.0271). Compared to the UKFt+rFA method, there was a significant increase of the volume of the AF tract that intersected the edema using the UKFt+FW method, while the whole AF volumes were similar. Expert judgment results, based on the five patients with the smallest AF volumes, indicated that the expert readers generally preferred the AF tract obtained by using the FW model, according to their anatomical knowledge and considering the potential influence of the final results on the surgical route.

Conclusion

Our results indicate that incorporating biophysical models of edema can increase the sensitivity of tractography in regions of peritumoral edema, allowing better tract visualization in patients with high grade gliomas and metastases.

SlicerDMRI: Open Source Diffusion MRI Software for Brain Cancer Research

Diffusion MRI (dMRI) is the only noninvasive method for mapping white matter connections in the brain. We describe SlicerDMRI, a software suite that enables visualization and analysis of dMRI for neuroscientific studies and patient-specific anatomic assessment. SlicerDMRI has been successfully applied in multiple studies of the human brain in health and disease, and here, we especially focus on its cancer research applications. As an extension module of the 3D Slicer medical image computing platform, the SlicerDMRI suite enables dMRI analysis in a clinically relevant multimodal imaging workflow. Core SlicerDMRI functionality includes diffusion tensor estimation, white matter tractography with single and multi-fiber models, and dMRI quantification. SlicerDMRI supports clinical DICOM and research file formats, is open-source and cross-platform, and can be installed as an extension to 3D Slicer (www.slicer.org). More information, videos, tutorials, and sample data are available at dmri.slicer.org Cancer Res; 77(21); e101-3. ©2017 AACR.

Whole brain white matter connectivity analysis using machine learning: An application to autism

In this paper, we propose an automated white matter connectivity analysis method for machine learning classification and characterization of white matter abnormality via identification of discriminative fiber tracts. The proposed method uses diffusion MRI tractography and a data-driven approach to find fiber clusters corresponding to subdivisions of the white matter anatomy. Features extracted from each fiber cluster describe its diffusion properties and are used for machine learning. The method is demonstrated by application to a pediatric neuroimaging dataset from 149 individuals, including 70 children with autism spectrum disorder (ASD) and 79 typically developing controls (TDC). A classification accuracy of 78.33% is achieved in this cross-validation study. We investigate the discriminative diffusion features based on a two-tensor fiber tracking model. We observe that the mean fractional anisotropy from the second tensor (associated with crossing fibers) is most affected in ASD. We also find that local along-tract (central cores and endpoint regions) differences between ASD and TDC are helpful in differentiating the two groups. These altered diffusion properties in ASD are associated with multiple robustly discriminative fiber clusters, which belong to several major white matter tracts including the corpus callosum, arcuate fasciculus, uncinate fasciculus and aslant tract; and the white matter structures related to the cerebellum, brain stem, and ventral diencephalon. These discriminative fiber clusters, a small part of the whole brain tractography, represent the white matter connections that could be most affected in ASD. Our results indicate the potential of a machine learning pipeline based on white matter fiber clustering.

Comparison of Diffusion-Weighted MRI Reconstruction Methods for Visualization of Cranial Nerves in Posterior Fossa Surgery

Diffusion-weighted imaging (DWI)-based tractography has gained increasing popularity as a method for detailed visualization of white matter (WM) tracts. Different imaging techniques, and more novel, advanced imaging methods provide significant WM structural detail. While there has been greater focus on improving tract visualization for larger WM pathways, the relative value of each method for cranial nerve reconstruction and how this methodology can assist surgical decision-making is still understudied. Images from 10 patients with posterior fossa tumors (4 male, mean age: 63.5), affecting either the trigeminal nerve (CN V) or the facial/vestibular complex (CN VII/VIII), were employed. Three distinct reconstruction methods [two tensor-based methods: single diffusion tensor tractography (SDT) (3D Slicer), eXtended streamline tractography (XST), and one fiber orientation distribution (FOD)-based method: streamline tractography using constrained spherical deconvolution (CSD)-derived estimates (MRtrix3)], were compared to determine which of these was best suited for use in a neurosurgical setting in terms of processing speed, anatomical accuracy, and accurate depiction of the relationship between the tumor and affected CN. Computation of the tensor map was faster when compared to the implementation of CSD to provide estimates of FOD. Both XST and CSD-based reconstruction methods tended to give more detailed representations of the projections of CN V and CN VII/VIII compared to SDT. These reconstruction methods were able to more accurately delineate the course of CN V and CN VII/VIII, differentiate CN V from the cerebellar peduncle, and delineate compression of CN VII/VIII in situations where SDT could not. However, CSD-based reconstruction methods tended to generate more invalid streamlines. XST offers the best combination of anatomical accuracy and speed of reconstruction of cranial nerves within this patient population. Given the possible anatomical limitations of single tensor models, supplementation with more advanced tensor-based reconstruction methods might be beneficial.

Performance of unscented Kalman filter tractography in edema: Analysis of the two-tensor model

Diffusion MRI tractography is increasingly used in pre-operative neurosurgical planning to visualize critical fiber tracts. However, a major challenge for conventional tractography, especially in patients with brain tumors, is tracing fiber tracts that are affected by vasogenic edema, which increases water content in the tissue and lowers diffusion anisotropy. One strategy for improving fiber tracking is to use a tractography method that is more sensitive than the traditional single-tensor streamline tractography.

We performed experiments to assess the performance of two-tensor unscented Kalman filter (UKF) tractography in edema. UKF tractography fits a diffusion model to the data during fiber tracking, taking advantage of prior information from the previous step along the fiber. We studied UKF performance in a synthetic diffusion MRI digital phantom with simulated edema and in retrospective data from two neurosurgical patients with edema affecting the arcuate fasciculus and corticospinal tracts. We compared the performance of several tractography methods including traditional streamline, UKF single-tensor, and UKF two-tensor. To provide practical guidance on how the UKF method could be employed, we evaluated the impact of using various seed regions both inside and outside the edematous regions, as well as the impact of parameter settings on the tractography sensitivity. We quantified the sensitivity of different methods by measuring the percentage of the patient-specific fMRI activation that was reached by the tractography.

We expected that diffusion anisotropy threshold parameters, as well as the inclusion of a free water model, would significantly influence the reconstruction of edematous WM fiber tracts, because edema increases water content in the tissue and lowers anisotropy. Contrary to our initial expectations, varying the fractional anisotropy threshold and including a free water model did not affect the UKF two-tensor tractography output appreciably in these two patient datasets. The most effective parameter for increasing tracking sensitivity was the generalized anisotropy (GA) threshold, which increased the length of tracked fibers when reduced to 0.075. In addition, the most effective seeding strategy was seeding in the whole brain or in a large region outside of the edema.

Overall, the main contribution of this study is to provide insight into how UKF tractography can work, using a two-tensor model, to begin to address the challenge of fiber tract reconstruction in edematous regions near brain tumors.

•

Reconstruction of edematous white matter from diffusion MRI is investigated.

•

The performance of two–tensor unscented Kalman filter (UKF) tractography is assessed.

•

The two–tensor model in UKF is analyzed in phantom and patient data experiments.

•

Practical guidance on employing the UKF method in neurosurgical patients is provided

White matter tractography for neurosurgical planning: A topography-based review of the current state of the art

We perform a review of the literature in the field of white matter tractography for neurosurgical planning, focusing on those works where tractography was correlated with clinical information such as patient outcome, clinical functional testing, or electro-cortical stimulation. We organize the review by anatomical location in the brain and by surgical procedure, including both supratentorial and infratentorial pathologies, and excluding spinal cord applications. Where possible, we discuss implications of tractography for clinical care, as well as clinically relevant technical considerations regarding the tractography methods. We find that tractography is a valuable tool in variable situations in modern neurosurgery. Our survey of recent reports demonstrates multiple potentially successful applications of white matter tractography in neurosurgery, with progress towards overcoming clinical challenges of standardization and interpretation.

Automated white matter fiber tract identification in patients with brain tumors

We propose a method for the automated identification of key white matter fiber tracts for neurosurgical planning, and we apply the method in a retrospective study of 18 consecutive neurosurgical patients with brain tumors. Our method is designed to be relatively robust to challenges in neurosurgical tractography, which include peritumoral edema, displacement, and mass effect caused by mass lesions. The proposed method has two parts. First, we learn a data-driven white matter parcellation or fiber cluster atlas using groupwise registration and spectral clustering of multi-fiber tractography from healthy controls. Key fiber tract clusters are identified in the atlas. Next, patient-specific fiber tracts are automatically identified using tractography-based registration to the atlas and spectral embedding of patient tractography.

Results indicate good generalization of the data-driven atlas to patients: 80% of the 800 fiber clusters were identified in all 18 patients, and 94% of the 800 fiber clusters were found in 16 or more of the 18 patients. Automated subject-specific tract identification was evaluated by quantitative comparison to subject-specific motor and language functional MRI, focusing on the arcuate fasciculus (language) and corticospinal tracts (motor), which were identified in all patients. Results indicate good colocalization: 89 of 95, or 94%, of patient-specific language and motor activations were intersected by the corresponding identified tract. All patient-specific activations were within 3mm of the corresponding language or motor tract. Overall, our results indicate the potential of an automated method for identifying fiber tracts of interest for neurosurgical planning, even in patients with mass lesions.

•

Spectral clustering machine learning approach for white matter tract identification

•

Data-driven white matter parcellation learned from healthy subjects tractography

•

White matter parcellation applied to 18 consecutive patients with brain tumors

•

Arcuate fasciculus and corticospinal tracts identified in all patients

•

All tracts within 3 mm of corresponding patient-specific functional activations