Stability metrics for optic radiation tractography: Towards damage prediction after resective surgery

Subject-Specific Automatic Reconstruction of White Matter Tracts

MRI-based tractography is still underexploited and unsuited for routine use in brain tumor surgery due to heterogeneity of methods and functional-anatomical definitions and above all, the lack of a turn-key system. Standardization of methods is therefore desirable, whereby an objective and reliable approach is a prerequisite before the results of any automated procedure can subsequently be validated and used in neurosurgical practice. In this work, we evaluated these preliminary but necessary steps in healthy volunteers. Specifically, we evaluated the robustness and reliability (i.e., test-retest reproducibility) of tractography results of six clinically relevant white matter tracts by using healthy volunteer data (N = 136) from the Human Connectome Project consortium. A deep learning convolutional network-based approach was used for individualized segmentation of regions of interest, combined with an evidence-based tractography protocol and appropriate post-tractography filtering. Robustness was evaluated by estimating the consistency of tractography probability maps, i.e., averaged tractograms in normalized space, through the use of a hold-out cross-validation approach. No major outliers were found, indicating a high robustness of the tractography results. Reliability was evaluated at the individual level. First by examining the overlap of tractograms that resulted from repeatedly processed identical MRI scans (N = 10, 10 iterations) to establish an upper limit of reliability of the pipeline. Second, by examining the overlap for subjects that were scanned twice at different time points (N = 40). Both analyses indicated high reliability, with the second analysis showing a reliability near the upper limit. The robust and reliable subject-specific generation of white matter tracts in healthy subjects holds promise for future validation of our pipeline in a clinical population and subsequent implementation in brain tumor surgery.

An atlas of white matter anatomy, its variability, and reproducibility based on constrained spherical deconvolution of diffusion MRI

Virtual dissection of white matter (WM) using diffusion MRI tractography is confounded by its poor reproducibility. Despite the increased adoption of advanced reconstruction models, early region-of-interest driven protocols based on diffusion tensor imaging (DTI) remain the dominant reference for virtual dissection protocols. Here we bridge this gap by providing a comprehensive description of typical WM anatomy reconstructed using a reproducible automated subject-specific parcellation-based approach based on probabilistic constrained-spherical deconvolution (CSD) tractography. We complement this with a WM template in MNI space comprising 68 bundles, including all associated anatomical tract selection labels and associated automated workflows. Additionally, we demonstrate bundle inter- and intra-subject variability using 40 (20 test-retest) datasets from the human connectome project (HCP) and 5 sessions with varying b-values and number of b-shells from the single-subject Multiple Acquisitions for Standardization of Structural Imaging Validation and Evaluation (MASSIVE) dataset. The most reliably reconstructed bundles were the whole pyramidal tracts, primary corticospinal tracts, whole superior longitudinal fasciculi, frontal, parietal and occipital segments of the corpus callosum and middle cerebellar peduncles. More variability was found in less dense bundles, e.g., the fornix, dentato-rubro-thalamic tract (DRTT), and premotor pyramidal tract. Using the DRTT as an example, we show that this variability can be reduced by using a higher number of seeding attempts. Overall inter-session similarity was high for HCP test-retest data (median weighted-dice = 0.963, stdev = 0.201 and IQR = 0.099). Compared to the HCP-template bundles there was a high level of agreement for the HCP test-retest data (median weighted-dice = 0.747, stdev = 0.220 and IQR = 0.277) and for the MASSIVE data (median weighted-dice = 0.767, stdev = 0.255 and IQR = 0.338). In summary, this WM atlas provides an overview of the capabilities and limitations of automated subject-specific probabilistic CSD tractography for mapping white matter fasciculi in healthy adults. It will be most useful in applications requiring a reproducible parcellation-based dissection protocol, and as an educational resource for applied neuroimaging and clinical professionals.

Fully automated delineation of the optic radiation for surgical planning using clinically feasible sequences

Quadrantanopia caused by inadvertent severing of Meyer's Loop of the optic radiation is a well-recognised complication of temporal lobectomy for conditions such as epilepsy. Dissection studies indicate that the anterior extent of Meyer's Loop varies considerably between individuals. Quantifying this for individual patients is thus an important step to improve the safety profile of temporal lobectomies. Previous attempts to delineate Meyer's Loop using diffusion MRI tractography have had difficulty estimating its full anterior extent, required manual ROI placement, and/or relied on advanced diffusion sequences that cannot be acquired routinely in most clinics. Here we present CONSULT: a pipeline that can delineate the optic radiation from raw DICOM data in a completely automated way via a combination of robust pre-processing, segmentation, and alignment stages, plus simple improvements that bolster the efficiency and reliability of standard tractography. We tested CONSULT on 696 scans of predominantly healthy participants (539 unique brains), including both advanced acquisitions and simpler acquisitions that could be acquired in clinically acceptable timeframes. Delineations completed without error in 99.4% of the scans. The distance between Meyer's Loop and the temporal pole closely matched both averages and ranges reported in dissection studies for all tested sequences. Median scan-rescan error of this distance was 1 mm. When tested on two participants with considerable pathology, delineations were successful and realistic. Through this, we demonstrate not only how to identify Meyer's Loop with clinically feasible sequences, but also that this can be achieved without fundamental changes to tractography algorithms or complex post-processing methods.

Fourier Transform on the Homogeneous Space of 3D Positions and Orientations for Exact Solutions to Linear PDEs

Fokker–Planck PDEs (including diffusions) for stable Lévy processes (including Wiener processes) on the joint space of positions and orientations play a major role in mechanics, robotics, image analysis, directional statistics and probability theory. Exact analytic designs and solutions are known in the 2D case, where they have been obtained using Fourier transform on SE(2). Here, we extend these approaches to 3D using Fourier transform on the Lie group SE(3) of rigid body motions. More precisely, we define the homogeneous space of 3D positions and orientations R3⋊S2:=SE(3)/({0}×SO(2)) as the quotient in SE(3). In our construction, two group elements are equivalent if they are equal up to a rotation around the reference axis. On this quotient, we design a specific Fourier transform. We apply this Fourier transform to derive new exact solutions to Fokker–Planck PDEs of α-stable Lévy processes on R3⋊S2. This reduces classical analysis computations and provides an explicit algebraic spectral decomposition of the solutions. We compare the exact probability kernel for α=1 (the diffusion kernel) to the kernel for α=12 (the Poisson kernel). We set up stochastic differential equations (SDEs) for the Lévy processes on the quotient and derive corresponding Monte-Carlo methods. We verified that the exact probability kernels arise as the limit of the Monte-Carlo approximations.

Topological principles and developmental algorithms might refine diffusion tractography

The identification and reconstruction of axonal pathways in the living brain or "ex-vivo" is promising a revolution in connectivity studies bridging the gap from animal to human neuroanatomy with extensions to brain structural-functional correlates. Unfortunately, the methods suffer from juvenile drawbacks. In this perspective paper we mention several computational and developmental principles, which might stimulate a new generation of algorithms and a discussion bridging the neuroimaging and neuroanatomy communities.

Meyer's loop tractography for image-guided surgery depends on imaging protocol and hardware

Introduction

Surgical resection is an effective treatment for temporal lobe epilepsy but can result in visual field defects. This could be minimized if surgeons knew the exact location of the anterior part of the optic radiation (OR), the Meyer's loop. To this end, there is increasing prevalence of image-guided surgery using diffusion MRI tractography. Despite considerable effort in developing analysis methods, a wide discrepancy in Meyer's loop reconstructions is observed in the literature. Moreover, the impact of differences in image acquisition on Meyer's loop tractography remains unclear. Here, while employing the same state-of-the-art analysis protocol, we explored the extent to which variance in data acquisition leads to variance in OR reconstruction.

Methods

Diffusion MRI data were acquired for the same thirteen healthy subjects using standard and state-of-the-art protocols on three scanners with different maximum gradient amplitudes (MGA): Siemens Connectom (MGA = 300 mT/m); Siemens Prisma (MGA = 80 mT/m) and GE Excite-HD (MGA = 40 mT/m). Meyer's loop was reconstructed on all subjects and its distance to the temporal pole (ML-TP) was compared across protocols.

Results

A significant effect of data acquisition on the ML-TP distance was observed between protocols (p < .01 to 0.0001). The biggest inter-acquisition discrepancy for the same subject across different protocols was 16.5 mm (mean: 9.4 mm, range: 3.7-16.5 mm).

Conclusion

We showed that variance in data acquisition leads to substantive variance in OR tractography. This has direct implications for neurosurgical planning, where part of the OR is at risk due to an under-estimation of its location using conventional acquisition protocols.