Search

Reproducible identification of white matter pathways across subjects is essential for the study of structural connectivity of the human brain. One of the key challenges is anatomical differences between subjects and human rater subjectivity in labeling. Labeling white matter regions of interest presents many challenges due to the need to integrate both local and global information. Clearly communicating the manual processes to capture this information is cumbersome, yet essential to lay a solid foundation for comprehensive atlases. Segmentation protocols must be designed so the interpretation of the requested tasks as well as locating structural landmarks is anatomically accurate, intuitive and reproducible. In this work, we quantified the reproducibility of a first iteration of an open/public multi-bundle segmentation protocol. This allowed us to establish a baseline for its reproducibility as well as to identify the limitations for future iterations. The protocol was tested/evaluated on both typical 3 T research acquisition Baltimore Longitudinal Study of Aging (BLSA) and high-acquisition quality Human Connectome Project (HCP) datasets. The results show that a rudimentary protocol can produce acceptable intra-rater and inter-rater reproducibility. However, this work highlights the difficulty in generalizing reproducible results and the importance of reaching consensus on anatomical description of white matter pathways. The protocol has been made available in open source to improve generalizability and reliability in collaboration. The goal is to improve upon the first iteration and initiate a discussion on the anatomical validity (or lack thereof) of some bundle definitions and the importance of reproducibility of tractography segmentation., (Copyright © 2021 Elsevier Inc. All rights reserved.)

Background: Achieving inter-site / inter-scanner reproducibility of diffusion weighted magnetic resonance imaging (DW-MRI) metrics has been challenging given differences in acquisition protocols, analysis models, and hardware factors., Purpose: Magnetic field gradients impart scanner-dependent spatial variations in the applied diffusion weighting that can be corrected if the gradient nonlinearities are known. However, retrieving manufacturer nonlinearity specifications is not well supported and may introduce errors in interpretation of units or coordinate systems. We propose an empirical approach to mapping the gradient nonlinearities with sequences that are supported across the major scanner vendors., Study Type: Prospective observational study., Subjects: A spherical isotropic diffusion phantom, and a single human control volunteer., Field Strength/sequence: 3 T (two scanners). Stejskal-Tanner spin echo sequence with b-values of 1000, 2000 s/mm 2 with 12, 32, and 384 diffusion gradient directions per shell., Assessment: We compare the proposed correction with the prior approach using manufacturer specifications against typical diffusion pre-processing pipelines (i.e., ignoring spatial gradient nonlinearities). In phantom data, we evaluate metrics against the ground truth. In human and phantom data, we evaluate reproducibility across scans, sessions, and hardware., Statistical Tests: Wilcoxon rank-sum test between uncorrected and corrected data., Results: In phantom data, our correction method reduces variation in mean diffusivity across sessions over uncorrected data (p < 0.05). In human data, we show that this method can also reduce variation in mean diffusivity across scanners (p < 0.05)., Conclusion: Our method is relatively simple, fast, and can be applied retroactively. We advocate incorporating voxel-specific b-value and b-vector maps should be incorporated in DW-MRI harmonization preprocessing pipelines to improve quantitative accuracy of measured diffusion parameters., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Purpose: Diffusion-weighted magnetic resonance imaging (DW-MRI) is of critical importance for characterizing in-vivo white matter. Models relating microarchitecture to observed DW-MRI signals as a function of diffusion sensitization are the lens through which DW-MRI data are interpreted. Numerous modern approaches offer opportunities to assess more complex intra-voxel structures. Nevertheless, there remains a substantial gap between intra-voxel estimated structures and ground truth captured by 3-D histology., Methods: Herein, we propose a novel data-driven approach to model the non-linear mapping between observed DW-MRI signals and ground truth structures using a sequential deep neural network regression using residual block deep neural network (ResDNN). Training was performed on two 3-D histology datasets of squirrel monkey brains and validated on a third. A second validation was performed using scan-rescan datasets of 12 subjects from Human Connectome Project. The ResDNN was compared with multiple micro-structure reconstruction methods and super resolved-constrained spherical deconvolution (sCSD) in particular as baseline for both the validations., Results: Angular correlation coefficient (ACC) is a correlation/similarity measure and can be interpreted as accuracy when compared with a ground truth. The median ACC of ResDNN is 0.82 and median ACC's of different variants of CSD are 0.75, 0.77, 0.79. The mean, median and std. of ResDNN & sCSD ACC across 12 subjects from HCP are 0.74, 0.88, 0.31 and 0.61, 0.71, 0.31 respectively., Conclusion: This work highlights the ability of deep learning to capture linkages between ex-vivo ground truth data with feasible MRI sequences. The data-driven approach is applicable to human in-vivo data and results in intriguingly high reproducibility of orientation structure., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Neuroimaging often involves acquiring high-resolution anatomical images along with other low-resolution image modalities, like diffusion and functional magnetic resonance imaging. Performing gray matter statistics with low-resolution image modalities is a challenge due to registration artifacts and partial volume effects. Gray matter surface based spatial statistics (GS-BSS) has been shown to provide higher sensitivity using gray matter surfaces compared to that of skeletonization approach of gray matter based spatial statistics which is adapted from tract based spatial statistics in diffusion studies. In this study, we improve upon GS-BSS incorporating neurite orientation dispersion and density imaging (NODDI) based search (denoted N-GSBSS) by 1) enhancing metrics mapping from native space, 2) incorporating maximum orientation dispersion index (ODI) search along surface normal, and 3) proposing applicability to other modalities, such as functional MRI (fMRI). We evaluated the performance of N-GSBSS against three baseline pipelines: volume-based registration, FreeSurfer's surface registration and ciftify pipeline for fMRI and simulation studies. First, qualitative mean ODI results are shown for N-GSBSS with and without NODDI based search in comparison with ciftify pipeline. Second, we conducted one-sample t-tests on working memory activations in fMRI to show that the proposed method can aid in the analysis of low resolution fMRI data. Finally we performed a sensitivity test in a simulation study by varying percentage change of intensity values within a region of interest in gray matter probability maps. N-GSBSS showed higher sensitivity in the simulation test compared to the other methods capturing difference between the groups starting at 10% change in the intensity values. The computational time of N-GSBSS is 68 times faster than that of traditional surface-based or 86 times faster than that of ciftify pipeline analysis., (Copyright © 2019 Elsevier Inc. All rights reserved.)

The ability to evaluate empirical diffusion MRI acquisitions for quality and to correct the resulting imaging metrics allows for improved inference and increased replicability. Previous work has shown promise for estimation of bias and variance of generalized fractional anisotropy (GFA) but comes at the price of computational complexity. This paper aims to provide methods for estimating GFA, bias of GFA and standard deviation of GFA quickly and accurately. In order to provide a method for bias and variance estimation that can return results faster than the previously studied statistical techniques, three deep, fully-connected neural networks are developed for GFA, bias of GFA, and standard deviation of GFA. The results of these networks are compared to the observed values of the metrics as well as those fit from the statistical techniques (i.e. Simulation Extrapolation (SIMEX) for bias estimation and wild bootstrap for variance estimation). Our GFA network provides predictions that are closer to the true GFA values than a Q-ball fit of the observed data (root-mean-square error (RMSE) 0.0077 vs 0.0082, p < .001). The bias network also shows statistically significant improvement in comparison to the SIMEX-estimated error of GFA (RMSE 0.0071 vs. 0.01, p < .001)., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Diffusion weighted MRI (DWMRI) and the myriad of analysis approaches (from tensors to spherical harmonics and brain tractography to body multi-compartment models) depend on accurate quantification of the apparent diffusion coefficient (ADC). Signal drift during imaging (e.g., due to b0 drift associated with heating) can cause systematic non-linearities that manifest as ADC changes if not corrected. Herein, we present a case study on two phantoms on one scanner. Different scan protocols exhibit different degrees of drift during similar scans and may be sensitive to the order of scans within an exam. Vos et al. recently reviewed the effects of signal drift in DWMRI acquisitions and proposed a temporal model for correction. We propose a novel spatial-temporal model to correct for higher order aspects of the signal drift and derive a statistically robust variant. We evaluate the Vos model and propose a method using two phantoms that mimic the ADC of the relevant brain tissue (0.36-2.2 × 10-3 mm 2 /s) on a single 3 T scanner. The phantoms are (1) a spherical isotropic sphere consisting of a single concentration of polyvinylpyrrolidone (PVP) and (2) an ice-water phantom with 13 vials of varying PVP concentrations. To characterize the impact of interspersed minimally weighted volumes ("b0's"), image volumes with b-value equal to 0.1 s/mm 2 are interspersed every 8, 16, 32, 48, and 96 diffusion weighted volumes in different trials. Signal drift is found to have spatially varying effects that are not accounted for with temporal-only models. The novel model captures drift more accurately (i.e., reduces the overall change per-voxel over the course of a scan) and results in more consistent ADC metrics., (Copyright © 2018 Elsevier Inc. All rights reserved.)

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/mm 2 ) 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/mm 2 ). 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/mm 2 . 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 (V ic ) 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., (Copyright © 2018 Elsevier Inc. All rights reserved.)

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/mm 2 ) 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., (Copyright © 2017 Elsevier Inc. All rights reserved.)