Faculty Bibliography

PURPOSE/OBJECTIVE:Rotational invariants (RIs) are at the root of many dMRI applications. Among others, they are presented as a sensible way of reducing the dimensionality of biophysical models. While thermal noise impact on diffusion metrics has been well studied, little is known on its effect on spherical harmonics-based RI (RISH) features and derived markers. In this work, we evaluate the effect of noise on RISH features and downstream Standard Model Imaging (SMI) estimates. THEORY AND METHODS/METHODS:Using simulated and test/retest multishell MRI data, we assess the accuracy and precision of RISH features and SMI parameters in the presence of thermal noise, as well as its robustness to variations in protocol design. We further propose and evaluate correction strategies that bypass the need of rotational invariant features as an intermediate step. RESULTS:Both RISH features and SMI estimates are impacted by SNR-dependent Rician biases. However, higher-order RISH features are susceptible to a secondary noise-related source of bias, which not only depends on SNR, but also protocol and underlying microstructure. Rician bias-correcting techniques are insufficient to maximize the accuracy of RISH and SMI features, or to ensure consistency across protocols. SMI estimators that avoid RISH features by fitting the model to the directional diffusion MRI data outperform RISH-based approaches in accuracy, repeatability, and reproducibility across acquisition protocols. CONCLUSIONS:RISH features are increasingly used in dMRI analysis, yet they are prone to various sources of noise that lower their accuracy and reproducibility. Understanding the impact of noise and mitigating such biases is critical to maximize the validity, repeatability, and reproducibility of dMRI studies.

The development of noninvasive MRI biomarkers as surrogates of histopathological features in kidney tissue requires detailed explorations of contrast. Therefore, we studied kidney diffusion kurtosis imaging (DKI) with a wide array of encodings, including flow compensation, variable directional sampling, and cardiac gating regimes. Twelve healthy volunteers underwent DKI at 5-10 diffusion weightings (b-values) ranging from 0 to 1200 smm^-2 with 12 or 30 directional samplings, bipolar or flow-compensated diffusion gradient waveforms, and at systolic or diastolic cardiac phases. DKI biomarkers, mean diffusivity (MD) and kurtosis (MK), were interrogated using a directionally robust fitting algorithm compared to conventional fits. The combination of flow compensation and cardiac triggering at the diastolic phase in the kidneys reduced flow effects on DKI. In systole, flow-compensated waveforms significantly reduced MD and MK for both cortex and medulla: cortex MD: 3.00 versus 2.55 μm² ms^-1, medulla MD: 2.80 versus 2.39 μm² ms^-1, cortex MK: 0.58 versus 0.45, and medulla MK: 0.60 versus 0.47 (all p < 0.05). Flow suppression alleviated requirements for processing the DKI at higher minimum b-values, as neither MD nor MK significantly differed at the diastolic phase for minimum b-values of 0 versus 200 smm^-2: cortex MD: 2.30 versus 2.28 μm² ms^-1, p = 0.278; medulla MD: 2.29 versus 2.28 μm² ms^-1, p = 0.437; cortex MK: 0.37 versus 0.36, p = 0.308; and medulla MK: 0.40 versus 0.40, p = 0.904. Flow-compensated waveforms mitigate cardiac and respiratory motion-related artifacts at higher diffusion encodings in addition to microcirculation effects. The robust fitting initially developed for brain DKI is highly applicable to the kidneys because it disentangles tissue-specific directional diffusion information from artifacts.

Ultra low-field (ULF) MRI is an accessible neuroimaging modality that can bridge healthcare disparities and advance population-level brain health research. However, the inherently low signal-to-noise ratio of ULF-MRI often necessitates reductions in spatial resolution and, combined with the field-dependency of MRI contrast, challenges the accurate extraction of clinically relevant brain morphology. We evaluate the current state of ULF-MRI brain volumetry utilizing techniques for enhancing spatial resolution and leveraging recent advancements in brain segmentation. This is based on the agreement between ULF and corresponding high-field (HF) MRI brain volumes, and test-retest repeatability for multiple ULF scans. In this study, we find that accurate brain volumes can be measured from ULF-MRIs when combining orthogonal imaging directions for T₂-weighted images to form a higher resolution image volume. We also demonstrate that not all orthogonal imaging directions contribute equally to volumetric accuracy and provide a recommended scan protocol given the constraints of the current technology.

Small-animal diffusion MRI (dMRI) has been used for methodological development and validation, characterizing the biological basis of diffusion phenomena, and comparative anatomy. The steps from animal setup and monitoring, to acquisition, analysis, and interpretation are complex, with many decisions that may ultimately affect what questions can be answered using the resultant data. This work aims to present selected considerations and recommendations from the diffusion community on best practices for preclinical dMRI of in vivo animals. We describe the general considerations and foundational knowledge that must be considered when designing experiments. We briefly describe differences in animal species and disease models and discuss why some may be more or less appropriate for different studies. We, then, give recommendations for in vivo acquisition protocols, including decisions on hardware, animal preparation, and imaging sequences, followed by advice for data processing including preprocessing, model-fitting, and tractography. Finally, we provide an online resource that lists publicly available preclinical dMRI datasets and software packages to promote responsible and reproducible research. In each section, we attempt to provide guides and recommendations, but also highlight areas for which no guidelines exist (and why), and where future work should focus. Although we mainly cover the central nervous system (on which most preclinical dMRI studies are focused), we also provide, where possible and applicable, recommendations for other organs of interest. An overarching goal is to enhance the rigor and reproducibility of small animal dMRI acquisitions and analyses, and thereby advance biomedical knowledge.

The value of preclinical diffusion MRI (dMRI) is substantial. While dMRI enables in vivo non-invasive characterization of tissue, ex vivo dMRI is increasingly being used to probe tissue microstructure and brain connectivity. Ex vivo dMRI has several experimental advantages including higher SNR and spatial resolution compared to in vivo studies, and enabling more advanced diffusion contrasts for improved microstructure and connectivity characterization. Another major advantage of ex vivo dMRI is the direct comparison with histological data, as a crucial methodological validation. However, there are a number of considerations that must be made when performing ex vivo experiments. The steps from tissue preparation, image acquisition and processing, and interpretation of results are complex, with many decisions that not only differ dramatically from in vivo imaging of small animals, but ultimately affect what questions can be answered using the data. This work represents "Part 2" of a three-part series of recommendations and considerations for preclinical dMRI. We describe best practices for dMRI of ex vivo tissue, with a focus on the value that ex vivo imaging adds to the field of dMRI and considerations in ex vivo image acquisition. We first give general considerations and foundational knowledge that must be considered when designing experiments. We briefly describe differences in specimens and models and discuss why some may be more or less appropriate for different studies. We then give guidelines for ex vivo protocols, including tissue fixation, sample preparation, and MR scanning. In each section, we attempt to provide guidelines and recommendations, but also highlight areas for which no guidelines exist (and why), and where future work should lie. An overarching goal herein is to enhance the rigor and reproducibility of ex vivo dMRI acquisitions and analyses, and thereby advance biomedical knowledge.

Preclinical diffusion MRI (dMRI) has proven value in methods development and validation, characterizing the biological basis of diffusion phenomena, and comparative anatomy. While dMRI enables in vivo non-invasive characterization of tissue, ex vivo dMRI is increasingly being used to probe tissue microstructure and brain connectivity. Ex vivo dMRI has several experimental advantages that facilitate high spatial resolution and high SNR images, cutting-edge diffusion contrasts, and direct comparison with histological data as a methodological validation. However, there are a number of considerations that must be made when performing ex vivo experiments. The steps from tissue preparation, image acquisition and processing, and interpretation of results are complex, with many decisions that not only differ dramatically from in vivo imaging of small animals, but ultimately affect what questions can be answered using the data. This work concludes a three-part series of recommendations and considerations for preclinical dMRI. Herein, we describe best practices for dMRI of ex vivo tissue, with a focus on image pre-processing, data processing, and comparisons with microscopy. In each section, we attempt to provide guidelines and recommendations but also highlight areas for which no guidelines exist (and why), and where future work should lie. We end by providing guidelines on code sharing and data sharing and point toward open-source software and databases specific to small animal and ex vivo imaging.

Tractography is a unique modality for the in vivo measurement of structural connectivity, crucial for understanding brain networks and neurological conditions. With increasing b-value, the diffusion-weighting signal becomes primarily sensitive to the intra-axonal signal. However, it remains unclear how tractography is affected by this observation. Here, using open-source datasets, we showed that at high b-values, DWI reduces the uncertainty in estimating fiber orientations. Specifically, we found the ratio of biologically-meaningful longer-range connections increases, accompanied with downstream impact of redistribution of connectome and network metrics. However, when going beyond b = 6000 s/mm², the loss of SNR imposed a penalty. Lastly, we showed that the data reaches satisfactory reproducibility with b-values above 1200 s/mm². Overall, the results suggest that using b-values above 2500 s/mm² is essential for more accurate connectome reconstruction by reducing uncertainty in fiber orientation estimation, supporting the use of higher b-value protocols in standard diffusion MRI scans and pipelines.

The clinical translation of diffusion magnetic resonance imaging (dMRI)-derived quantitative contrasts hinges on robust reproducibility, minimizing both same-scanner and cross-scanner variability. As multi-site data sets, including multi-shell dMRI, expand in scope, enhancing reproducibility across variable MRI systems and MRI protocols becomes crucial. This study evaluates the reproducibility of diffusion kurtosis imaging (DKI) metrics (beyond conventional diffusion tensor imaging (DTI)), at the voxel and region-of-interest (ROI) levels on magnitude and complex-valued dMRI data, using denoising with and without harmonization. We compared same-scanner, cross-scanner, and cross-protocol variability for a multi-shell dMRI protocol (2-mm isotropic resolution, b = 0, 1000, 2000 s/mm²) in 20 subjects. We first evaluated the effectiveness of Marchenko-Pastur Principal Component Analysis (MPPCA) based denoising strategies for both magnitude and complex data to mitigate noise-induced bias and variance, to improve dMRI parametric maps and reproducibility. Next, we examined the impact of denoising under different population analysis approaches, specifically comparing voxel-wise versus region of interest (ROI)-based methods. We also evaluated the role of denoising when harmonizing dMRI across scanners and protocols. The results indicate that DTI and DKI maps visually improve after MPPCA denoising, with noticeably fewer outliers in kurtosis maps. Denoising, either using magnitude or complex dMRI, enhances voxel-wise reproducibility, with test-retest variability of kurtosis indices reduced from 15%-20% without denoising to 5%-10% after denoising. Complex dMRI denoising reduces the noise floor by up to 60%. Denoising not only reduced variability across scans and protocols, but also increased statistical power for low SNR voxel-wise comparisons when comparing cross sectional groups. In conclusion, MPPCA denoising, either over magnitude or complex dMRI data, enhances the reproducibility and precision of higher-order diffusion metrics across same-scanner, cross-scanner, and cross-protocol assessments. The enhancement in data quality and precision facilitates the broader application and acceptance of these advanced imaging techniques in both clinical practice and large-scale neuroimaging studies.

The Adolescent Brain and Cognitive Development (ABCD) project is the largest study of adolescent brain development. ABCD longitudinally tracks 11,868 participants aged 9-10 years from 21 sites using standardized protocols for multi-site MRI data collection and analysis. While the multi-site and multi-scanner study design enhances the robustness and generalizability of analysis results, it may also introduce nonbiological variances including scanner-related variations, subject motion, and deviations from protocols. ABCD imaging data were collected biennially within a period of ongoing maturation in cortical thickness and integrity of cerebral white matter. These changes can bias the classical test-retest methodologies, such as intraclass correlation coefficients (ICC). We developed a site-wise adaptive ICC (AICC) to evaluate the reliability of imaging-derived phenotypes while accounting for ongoing brain development. AICC iteratively estimates the population-level age-related brain development trajectory using a weighted mixed model and updates age-corrected site-wise reliability until convergence. We evaluated the test-retest reliability of regional fractional anisotropy (FA) measures from diffusion tensor imaging and cortical thickness (CT) from structural MRI data for each site. The mean AICC for 20 FA tracts across sites was 0.61 ± 0.19, lower than the mean AICC for CT in 34 regions across sites, 0.76 ± 0.12. Remarkably, sites using Siemens scanners consistently showed significantly higher AICC values compared with those using GE/Philips scanners for both FA (AICC = 0.71 ± 0.12 vs. 0.46 ± 0.17, p < 0.001) and CT (AICC = 0.80 ± 0.10 vs. 0.69 ± 0.11, p < 0.001). These findings demonstrate site-and-scanner related variations in data quality and underscore the necessity for meticulous data curation in subsequent association analyses.

Joint modeling of diffusion and relaxation has seen growing interest due to its potential to provide complementary information about tissue microstructure. For brain white matter (WM), we designed an optimal diffusion-relaxometry MRI protocol that samples multiple b-values, B-tensor shapes, and echo times (TE). This variable-TE protocol (27 min) has as subsets a fixed-TE protocol (15 min) and a two-shell dMRI protocol (7 min), both characterizing diffusion only. We assessed the sensitivity, specificity, and reproducibility of these protocols with synthetic experiments and in six healthy volunteers. Compared with the fixed-TE protocol, the variable-TE protocol enables estimation of the free water fraction while also capturing compartmental