Vishwesh Nath

Empirical consideration of the effects of acquisition parameters and analysis model on clinically feasible q-ball imaging

Q-ball imaging (QBI) is a popular high angular resolution diffusion imaging (HARDI) technique used to study brain architecture in vivo. Simulation and phantom-based studies suggest that QBI results are affected by the b-value, the number of diffusion weighting directions, and the signal-to-noise ratio (SNR). However, optimal acquisition schemes for QBI in clinical settings are largely undetermined given empirical (observed) imaging considerations. In this study, we acquire a HARDI dataset at five b-values with 11 repetitions on a single subject to investigate the effects of acquisition scheme and subsequent analysis models on the accuracy and precision of measures of tissue composition and fiber orientation derived from clinically feasible QBI at 3T. Clinical feasibility entails short scan protocols - less than 5minutes in the current study - resulting in lower SNR, lower b-values, and fewer diffusion directions than are typical in most QBI protocols with research applications, where time constraints are less prevalent. In agreement with previous studies, we find that the b-value and number of diffusion directions impact the magnitude and variation of QBI indices in both white matter and gray matter regions; however, QBI indices are most heavily dependent on the maximum order of the spherical harmonic (SH) series used to represent the diffusion orientation distribution function (ODF). Specifically, to ensure numerical stability and reduce the occurrence of false peaks and inflated anisotropy, we recommend oversampling by at least 8-12 more diffusion directions than the number of estimated coefficients for a given SH order. In addition, in an equal scan time comparison of QBI accuracy, we find that increasing the directional resolution of the HARDI dataset is preferable to repeating observations; however, our results indicate that as few as 32 directions at a low b-value (1000s/mm2) captures most of the angular information in the q-ball ODF. Our findings provide guidance for determining an optimal acquisition scheme for QBI in the low SNR and low scan time regime, and suggest that care must be taken when choosing the basis functions used to represent the QBI ODF.

Effects of b-value and number of gradient directions on diffusion MRI measures obtained with Q-ball imaging

High-angular-resolution diffusion-weighted imaging (HARDI) MRI acquisitions have become common for use with higher order models of diffusion. Despite successes in resolving complex fiber configurations and probing microstructural properties of brain tissue, there is no common consensus on the optimal b-value and number of diffusion directions to use for these HARDI methods. While this question has been addressed by analysis of the diffusion-weighted signal directly, it is unclear how this translates to the information and metrics derived from the HARDI models themselves. Using a high angular resolution data set acquired at a range of b-values, and repeated 11 times on a single subject, we study how the b-value and number of diffusion directions impacts the reproducibility and precision of metrics derived from Q-ball imaging, a popular HARDI technique. We find that Q-ball metrics associated with tissue microstructure and white matter fiber orientation are sensitive to both the number of diffusion directions and the spherical harmonic representation of the Q-ball, and often are biased when under sampled. These results can advise researchers on appropriate acquisition and processing schemes, particularly when it comes to optimizing the number of diffusion directions needed for metrics derived from Q-ball imaging.

Limits to anatomical accuracy of diffusion tractography using modern approaches

&NA; Diffusion MRI fiber tractography is widely used to probe the structural connectivity of the brain, with a range of applications in both clinical and basic neuroscience. Despite widespread use, tractography has well‐known pitfalls that limits the anatomical accuracy of this technique. Numerous modern methods have been developed to address these shortcomings through advances in acquisition, modeling, and computation. To test whether these advances improve tractography accuracy, we organized the 3‐D Validation of Tractography with Experimental MRI (3D‐VoTEM) challenge at the ISBI 2018 conference. We made available three unique independent tractography validation datasets – a physical phantom and two ex vivo brain specimens ‐ resulting in 176 distinct submissions from 9 research groups. By comparing results over a wide range of fiber complexities and algorithmic strategies, this challenge provides a more comprehensive assessment of tractographys inherent limitations than has been reported previously. The central results were consistent across all sub‐challenges in that, despite advances in tractography methods, the anatomical accuracy of tractography has not dramatically improved in recent years. Taken together, our results independently confirm findings from decades of tractography validation studies, demonstrate inherent limitations in reconstructing white matter pathways using diffusion MRI data alone, and highlight the need for alternative or combinatorial strategies to accurately map the fiber pathways of the brain. HighlightsOrganized international tractography challenge utilizing three validation datasets.Anatomical accuracy of modern diffusion tractography techniques is limited.Advancements are needed to overcome limited sensitivity/specificity of reconstructions.

SHARD: spherical harmonic-based robust outlier detection for HARDI methods

High Angular Resolution Diffusion Imaging (HARDI) models are used to capture complex intra-voxel microarchitectures. The magnetic resonance imaging sequences that are sensitized to diffusion are often highly accelerated and prone to motion, physiologic, and imaging artifacts. In diffusion tensor imaging, robust statistical approaches have been shown to greatly reduce these adverse factors without human intervention. Similar approaches would be possible with HARDI methods, but robust versions of each distinct HARDI approach would be necessary. To avoid the computational and pragmatic burdens of creating individual robust HARDI analysis variants, we propose a robust outlier imputation model to mitigate outliers prior to traditional HARDI analysis. This model uses a weighted spherical harmonic fit of diffusion weighted magnetic resonance imaging scans to estimate the values which had been corrupted during acquisition to restore them. Briefly, spherical harmonics of 6th order were used to generate basis function which were weighted by diffusion signal for detection of outliers. For validation, a single healthy volunteer was scanned for a single session comprising of two scans one without head movement and the other with deliberate head movement at a b-value of 3000 s/mm2 with 64 diffusion weighted directions with a single b0 (5 averages) per scan. The deliberate motion from the volunteer created natural artifacts in the acquisition of one of the scans. The imputation model shows reduction in root mean squared error of the raw signal intensities and improvement for the HARDI method Q-ball in terms of the Angular Correlation Coefficient. The results reveal that there is quantitative and qualitative improvement. The proposed model can be used as general pre-processing model before implementing any HARDI model in general to restore the artifacts which are created because of the outlier diffusion signal in certain gradient volumes.

Evaluation of inter-site bias and variance in diffusion-weighted MRI

An understanding of the bias and variance of diffusion weighted magnetic resonance imaging (DW-MRI) acquisitions across scanners, study sites, or over time is essential for the incorporation of multiple data sources into a single clinical study. Studies that combine samples from various sites may be introducing confounding factors due to site-specific artifacts and patterns. Differences in bias and variance across sites may render the scans incomparable, and, without correction, inferences obtained from these data may be misleading. We present an analysis of the bias and variance of scans of the same subjects across different sites and evaluate their impact on statistical analyses. In previous work, we presented a simulation extrapolation (SIMEX) technique for bias estimation as well as a wild bootstrap technique for variance estimation in metrics obtained from a Q-ball imaging (QBI) reconstruction of empirical high angular resolution diffusion imaging (HARDI) data. We now apply those techniques to data acquired from 5 healthy volunteers on 3 independent scanners under closely matched acquisition protocols. The bias and variance of GFA measurements were estimated on a voxel-wise basis for each scan and compared across study sites to identify site-specific differences. Further, we provide model recommendations that can be used to determine the extent of the impact of bias and variance as well as aspects of the analysis to account for these differences. We include a decision tree to help researchers determine if model adjustments are necessary based on the bias and variance results.

Empirical single sample quantification of bias and variance in Q‐ball imaging

The bias and variance of high angular resolution diffusion imaging methods have not been thoroughly explored in the literature and may benefit from the simulation extrapolation (SIMEX) and bootstrap techniques to estimate bias and variance of high angular resolution diffusion imaging metrics.

Empirical reproducibility, sensitivity, and optimization of acquisition protocol, for Neurite Orientation Dispersion and Density Imaging using AMICO

Neurite Orientation Dispersion and Density Imaging (NODDI) has been gaining prominence for estimating multiple diffusion compartments from MRI data acquired in a clinically feasible time. To establish a pathway for adoption of NODDI in clinical studies, it is important to understand the sensitivity and reproducibility of NODDI metrics on empirical data in the context of acquisition protocol and brain anatomy. Previous studies addressed reproducibility across the 3 T scanners and within session and between subject reproducibility at 1.5 T and 3 T. However, empirical reproducibility on the performance of NODDI metrics based on b-value and diffusion-sensitized directions has not yet been addressed. In this study, we investigate a high angular resolution dataset with 11 repeats of a study with five b-values shells (1000, 1500, 2000, 2500 and 3000 s/mm2) and 96 directions per shell on a single subject. We validated the findings with a dataset from second subject with 10 repeats and 3 b-value shells (1000, 2000, 3000 s/mm2). The NODDI model was estimated using Accelerated Microstructure Imaging via Convex Optimization (AMICO) for different b-values and gradient directions on two-shell High Angular Resolution Density Imaging (HARDI) data fixing the lower shell at b = 1000 s/mm2. NODDI model applied to all acquired imaging data was used as a baseline gold standard for comparison. Additionally, we characterize orientation dispersion index (ODI) reproducibility using single-shell data. The experimental findings confirmed the sensitivity of intracellular volume fraction (Vic) with the choice of outer shell b-value more than with the choice of gradient directions. On the other hand, ODI is more sensitive to the number of gradient directions compared to b-value selection. Single-shell results for ODI are more comparable to 2-shell data at lower b-values than higher b-values. Recommended settings by region of interest and acquisition time are reported for the researchers considering using NODDI in human studies and/or comparing results across acquisition protocols.

Towards Portable Large-Scale Image Processing with High-Performance Computing

High-throughput, large-scale medical image computing demands tight integration of high-performance computing (HPC) infrastructure for data storage, job distribution, and image processing. The Vanderbilt University Institute for Imaging Science (VUIIS) Center for Computational Imaging (CCI) has constructed a large-scale image storage and processing infrastructure that is composed of (1) a large-scale image database using the eXtensible Neuroimaging Archive Toolkit (XNAT), (2) a content-aware job scheduling platform using the Distributed Automation for XNAT pipeline automation tool (DAX), and (3) a wide variety of encapsulated image processing pipelines called “spiders.” The VUIIS CCI medical image data storage and processing infrastructure have housed and processed nearly half-million medical image volumes with Vanderbilt Advanced Computing Center for Research and Education (ACCRE), which is the HPC facility at the Vanderbilt University. The initial deployment was natively deployed (i.e., direct installations on a bare-metal server) within the ACCRE hardware and software environments, which lead to issues of portability and sustainability. First, it could be laborious to deploy the entire VUIIS CCI medical image data storage and processing infrastructure to another HPC center with varying hardware infrastructure, library availability, and software permission policies. Second, the spiders were not developed in an isolated manner, which has led to software dependency issues during system upgrades or remote software installation. To address such issues, herein, we describe recent innovations using containerization techniques with XNAT/DAX which are used to isolate the VUIIS CCI medical image data storage and processing infrastructure from the underlying hardware and software environments. The newly presented XNAT/DAX solution has the following new features: (1) multi-level portability from system level to the application level, (2) flexible and dynamic software development and expansion, and (3) scalable spider deployment compatible with HPC clusters and local workstations.

Comparison of multi-fiber reproducibility of PAS-MRI and Q-ball with empirical multiple b-value HARDI

Crossing fibers are prevalent in human brains and a subject of intense interest for neuroscience. Diffusion tensor imaging (DTI) can resolve tissue orientation but is blind to crossing fibers. Many advanced diffusion-weighted magnetic resolution imaging (MRI) approaches have been presented to extract crossing-fibers from high angular resolution diffusion imaging (HARDI), but the relative sensitivity and specificity of approaches remains unclear. Here, we examine two leading approaches (PAS and q-ball) in the context of a large-scale, single subject reproducibility study. A single healthy individual was scanned 11 times with 96 diffusion weighted directions and 10 reference volumes for each of five b-values (1000, 1500, 2000, 2500, 3000 s/mm2) for a total of 5830 volumes (over the course of three sessions). We examined the reproducibility of the number of fibers per voxel, volume fraction, and crossing-fiber angles. For each method, we determined the minimum resolvable angle for each acquisition. Reproducibility of fiber counts per voxel was generally high (~80% consensus for PAS and ~70% for q-ball), but there was substantial bias between individual repetitions and model estimated with all data (~10% lower consensus for PAS and ~15% lower for q-ball). Both PAS and q-ball predominantly discovered fibers crossing at near 90 degrees, but reproducibility was higher for PAS across most measures. Within voxels with low anisotropy, q-ball finds more intra-voxel structure; meanwhile, PAS resolves multiple fibers at greater than 75 degrees for more voxels. These results can inform researchers when deciding between HARDI approaches or interpreting findings across studies.