An Immersive Virtual Reality Environment for Diagnostic Imaging

Lorenz Epprecht, Ahad Qureshi, Elliott D Kozin, Nicolas Vachicouras, Alexander M Huber, Ron Kikinis, Nikos Makris, Christian M Brown, Katherine L Reinshagen, and Daniel J Lee. 1/2020. “Human Cochlear Nucleus on 7 Tesla Diffusion Tensor Imaging: Insights Into Micro-anatomy and Function for Auditory Brainstem Implant Surgery.” Otol Neurotol.Abstract

OBJECTIVE: The cochlear nucleus (CN) is the target of the auditory brainstem implant (ABI). Most ABI candidates have Neurofibromatosis Type 2 (NF2) and distorted brainstem anatomy from bilateral vestibular schwannomas. The CN is difficult to characterize as routine structural MRI does not resolve detailed anatomy. We hypothesize that diffusion tensor imaging (DTI) enables both in vivo localization and quantitative measurements of CN morphology. STUDY DESIGN: We analyzed 7 Tesla (T) DTI images of 100 subjects (200 CN) and relevant anatomic structures using an MRI brainstem atlas with submillimetric (50 μm) resolution. SETTING: Tertiary referral center. PATIENTS: Young healthy normal hearing adults. INTERVENTION: Diagnostic. MAIN OUTCOME MEASURES: Diffusion scalar measures such as fractional anisotropy (FA), mean diffusivity (MD), mode of anisotropy (Mode), principal eigenvectors of the CN, and the adjacent inferior cerebellar peduncle (ICP). RESULTS: The CN had a lamellar structure and ventral-dorsal fiber orientation and could be localized lateral to the inferior cerebellar peduncle (ICP). This fiber orientation was orthogonal to tracts of the adjacent ICP where the fibers run mainly caudal-rostrally. The CN had lower FA compared to the medial aspect of the ICP (0.44 ± 0.09 vs. 0.64 ± 0.08, p < 0.001). CONCLUSIONS: 7T DTI enables characterization of human CN morphology and neuronal substructure. An ABI array insertion vector directed more caudally would better correspond to the main fiber axis of CN. State-of-the-art DTI has implications for ABI preoperative planning and future image guidance-assisted placement of the electrode array.

Guoqiang Xie, Fan Zhang, Laura Leung, Michael A Mooney, Lorenz Epprecht, Isaiah Norton, Yogesh Rathi, Ron Kikinis, Ossama Al-Mefty, Nikos Makris, Alexandra J Golby, and Lauren J O'Donnell. 1/2020. “Anatomical Assessment of Trigeminal Nerve Tractography Using Diffusion MRI: A Comparison of Acquisition B-Values and Single- and Multi-Fiber Tracking Strategies.” Neuroimage Clin, 25, Pp. 102160.Abstract

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.

Ørjan Bergmann, Rafael Henriques, Carl-Fredrik Westin, and Ofer Pasternak. 3/2020. “Fast and Accurate Initialization of the Free-water Imaging Model Parameters from Multi-shell Diffusion MRI.” NMR Biomed, 33, 3, Pp. e4219.Abstract

Cerebrospinal fluid partial volume effect is a known bias in the estimation of Diffusion Tensor Imaging (DTI) parameters from diffusion MRI data. The Free-Water Imaging model for diffusion MRI data adds a second compartment to the DTI model, which explicitly accounts for the signal contribution of extracellular free-water, such as cerebrospinal fluid. As a result the DTI parameters obtained through the free-water model are corrected for partial volume effects, and thus better represent tissue microstructure. In addition, the model estimates the fractional volume of free-water, and can be used to monitor changes in the extracellular space. Under certain assumptions, the model can be estimated from single-shell diffusion MRI data. However, by using data from multi-shell diffusion acquisitions, these assumptions can be relaxed, and the fit becomes more robust. Nevertheless, fitting the model to multi-shell data requires high computational cost, with a non-linear iterative minimization, which has to be initialized close enough to the global minimum to avoid local minima and to robustly estimate the model parameters. Here we investigate the properties of the main initialization approaches that are currently being used, and suggest new fast approaches to improve the initial estimates of the model parameters. We show that our proposed approaches provide a fast and accurate initial approximation of the model parameters, which is very close to the final solution. We demonstrate that the proposed initializations improve the final outcome of non-linear model fitting.

Laurent Chauvin, Kuldeep Kumar, Christian Wachinger, Marc Vangel, Jacques de Guise, Christian Desrosiers, William Wells, and Matthew Toews. 1/2020. “Neuroimage Signature from Salient Keypoints is Highly Specific to Individuals and Shared by Close Relatives.” Neuroimage, 204, Pp. 116208.Abstract

Neuroimaging studies typically adopt a common feature space for all data, which may obscure aspects of neuroanatomy only observable in subsets of a population, e.g. cortical folding patterns unique to individuals or shared by close relatives. Here, we propose to model individual variability using a distinctive keypoint signature: a set of unique, localized patterns, detected automatically in each image by a generic saliency operator. The similarity of an image pair is then quantified by the proportion of keypoints they share using a novel Jaccard-like measure of set overlap. Experiments demonstrate the keypoint method to be highly efficient and accurate, using a set of 7536 T1-weighted MRIs pooled from four public neuroimaging repositories, including twins, non-twin siblings, and 3334 unique subjects. All same-subject image pairs are identified by a similarity threshold despite confounds including aging and neurodegenerative disease progression. Outliers reveal previously unknown data labeling inconsistencies, demonstrating the usefulness of the keypoint signature as a computational tool for curating large neuroimage datasets.

Samo Lasič, Filip Szczepankiewicz, Erica Dall'Armellina, Arka Das, Christopher Kelly, Sven Plein, Jürgen E Schneider, Markus Nilsson, and Irvin Teh. 2/2020. “Motion-compensated b-tensor Encoding for in vivo Cardiac Diffusion-weighted Imaging.” NMR Biomed, 33, 2, Pp. e4213.Abstract

Motion is a major confound in diffusion-weighted imaging (DWI) in the body, and it is a common cause of image artefacts. The effects are particularly severe in cardiac applications, due to the nonrigid cyclical deformation of the myocardium. Spin echo-based DWI commonly employs gradient moment-nulling techniques to desensitise the acquisition to velocity and acceleration, ie, nulling gradient moments up to the 2nd order (M2-nulled). However, current M2-nulled DWI scans are limited to encode diffusion along a single direction at a time. We propose a method for designing b-tensors of arbitrary shapes, including planar, spherical, prolate and oblate tensors, while nulling gradient moments up to the 2nd order and beyond. The design strategy comprises initialising the diffusion encoding gradients in two encoding blocks about the refocusing pulse, followed by appropriate scaling and rotation, which further enables nulling undesired effects of concomitant gradients. Proof-of-concept assessment of in vivo mean diffusivity (MD) was performed using linear and spherical tensor encoding (LTE and STE, respectively) in the hearts of five healthy volunteers. The results of the M2-nulled STE showed that (a) the sequence was robust to cardiac motion, and (b) MD was higher than that acquired using standard M2-nulled LTE, where diffusion-weighting was applied in three orthogonal directions, which may be attributed to the presence of restricted diffusion and microscopic diffusion anisotropy. Provided adequate signal-to-noise ratio, STE could significantly shorten estimation of MD compared with the conventional LTE approach. Importantly, our theoretical analysis and the proposed gradient waveform design may be useful in microstructure imaging beyond diffusion tensor imaging where the effects of motion must be suppressed.

Neuroimaging Analysis Center

"understanding the human brain through imaging"

An Immersive Virtual Reality Environment for Diagnostic Imaging

Citation:

Recent Publications