In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI.

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency, ω $$ \omega $$ , in addition to the diffusion tensor, D $$ \mathbf{D} $$ , and relaxation, R 1 $$ {R}_1 $$ , R 2 $$ {R}_2 $$ , correlations. A D ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on their D ω - R 1 - R 2 $$ \mathbf{D}\left(\omega \right)-{R}_1-{R}_2 $$ distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.

"Pscore": A Novel Percentile-Based Metric to Accurately Assess Individual Deviations in Non-Gaussian Distributions of Quantitative MRI Metrics.

BACKGROUND: Quantitative magnetic resonance imaging (MRI) metrics could be used in personalized medicine to assess individuals against normative distributions. Conventional Zscore analysis is inadequate in the presence of non-Gaussian distributions. Therefore, if quantitative MRI metrics deviate from normality, an alternative is needed. PURPOSE: To confirm non-Gaussianity of diffusion MRI (dMRI) metrics on a publicly available dataset, and to propose a novel percentile-based method, "Pscore" to address this issue. STUDY TYPE: Retrospective cohort. POPULATION: Nine hundred and sixty-one healthy young adults (age: 22-35 years, females: 53%) from the Human Connectome Project. FIELD STRENGTH/SEQUENCE: 3-T, spin-echo diffusion echo-planar imaging, T1-weighted: MPRAGE. ASSESSMENT: The dMRI data were preprocessed using the TORTOISE pipeline. Forty-eight regions of interest (ROIs) from the JHU atlas were redrawn on a study-specific diffusion tensor (DT) template and average values were computed from various DT and mean apparent propagator (MAP) metrics. For each ROI, percentile ranks across participants were computed to generate "Pscores"-which normalized the difference between the median and a participant's value with the corresponding difference between the median and the 5th/95th percentile values. STATISTICAL TESTS: ROI-wise distributions were assessed using log transformations, Zscore, and the "Pscore" methods. The percentages of extreme values above-95th and below-5th percentile boundaries (PEV>95(%), PEV<5(%)) were also assessed in the overall white matter. Bootstrapping was performed to test the reliability of Pscores in small samples (N = 100) using 100 iterations. RESULTS: The dMRI metric distributions were systematically non-Gaussian, including positively skewed (eg, mean and radial diffusivity) and negatively skewed (eg, fractional and propagator anisotropy) metrics. This resulted in unbalanced tails in Zscore distributions (PEV>95 ≠ 5%, PEV<5 ≠ 5%) whereas "Pscore" distributions were symmetric and balanced (PEV>95 = PEV<5 = 5%); even for small bootstrapped samples (average PEV > 95 ¯ = PEV < 5 ¯ = 5 ± 0 % [SD]). DATA CONCLUSION: The inherent skewness observed for dMRI metrics may preclude the use of conventional Zscore analysis. The proposed "Pscore" method may help estimating individual deviations more accurately in skewed normative data, even from small datasets. LEVEL OF EVIDENCE: 1 TECHNICAL EFFICACY: Stage 1.

Neuroimaging Findings in US Government Personnel and Their Family Members Involved in Anomalous Health Incidents.

Importance: US government personnel stationed internationally have reported anomalous health incidents (AHIs), with some individuals experiencing persistent debilitating symptoms. Objective: To assess the potential presence of magnetic resonance imaging (MRI)-detectable brain lesions in participants with AHIs, with respect to a well-matched control group. Design, Setting, and Participants: This exploratory study was conducted at the National Institutes of Health (NIH) Clinical Center and the NIH MRI Research Facility between June 2018 and November 2022. Eighty-one participants with AHIs and 48 age- and sex-matched control participants, 29 of whom had similar employment as the AHI group, were assessed with clinical, volumetric, and functional MRI. A high-quality diffusion MRI scan and a second volumetric scan were also acquired during a different session. The structural MRI acquisition protocol was optimized to achieve high reproducibility. Forty-nine participants with AHIs had at least 1 additional imaging session approximately 6 to 12 months from the first visit. Exposure: AHIs. Main Outcomes and Measures: Group-level quantitative metrics obtained from multiple modalities: (1) volumetric measurement, voxel-wise and region of interest (ROI)-wise; (2) diffusion MRI-derived metrics, voxel-wise and ROI-wise; and (3) ROI-wise within-network resting-state functional connectivity using functional MRI. Exploratory data analyses used both standard, nonparametric tests and bayesian multilevel modeling. Results: Among the 81 participants with AHIs, the mean (SD) age was 42 (9) years and 49% were female; among the 48 control participants, the mean (SD) age was 43 (11) years and 42% were female. Imaging scans were performed as early as 14 days after experiencing AHIs with a median delay period of 80 (IQR, 36-544) days. After adjustment for multiple comparisons, no significant differences between participants with AHIs and control participants were found for any MRI modality. At an unadjusted threshold (P < .05), compared with control participants, participants with AHIs had lower intranetwork connectivity in the salience networks, a larger corpus callosum, and diffusion MRI differences in the corpus callosum, superior longitudinal fasciculus, cingulum, inferior cerebellar peduncle, and amygdala. The structural MRI measurements were highly reproducible (median coefficient of variation <1% across all global volumetric ROIs and <1.5% for all white matter ROIs for diffusion metrics). Even individuals with large differences from control participants exhibited stable longitudinal results (typically, <±1% across visits), suggesting the absence of evolving lesions. The relationships between the imaging and clinical variables were weak (median Spearman ρ = 0.10). The study did not replicate the results of a previously published investigation of AHIs. Conclusions and Relevance: In this exploratory neuroimaging study, there were no significant differences in imaging measures of brain structure or function between individuals reporting AHIs and matched control participants after adjustment for multiple comparisons.

Mapping gradient nonlinearity and miscalibration using diffusion-weighted MR images of a uniform isotropic phantom.

PURPOSE: To use diffusion measurements to map the spatial dependence of the magnetic field produced by the gradient coils of an MRI scanner with sufficient accuracy to correct errors in quantitative diffusion MRI (DMRI) caused by gradient nonlinearity and gradient amplifier miscalibration. THEORY AND METHODS: The field produced by the gradient coils is expanded in regular solid harmonics. The expansion coefficients are found by fitting a model to a minimum set of diffusion-weighted images of an isotropic diffusion phantom. The accuracy of the resulting gradient coil field maps is evaluated by using them to compute corrected b-matrices that are then used to process a multi-shell diffusion tensor imaging (DTI) dataset with 32 diffusion directions per shell. RESULTS: The method substantially reduces both the spatial inhomogeneity of the computed mean diffusivities (MD) and the computed values of the fractional anisotropy (FA), as well as virtually eliminating any artifactual directional bias in the tensor field secondary to gradient nonlinearity. When a small scaling miscalibration was purposely introduced in the x, y, and z, the method accurately detected the amount of miscalibration on each gradient axis. CONCLUSION: The method presented detects and corrects the effects of gradient nonlinearity and gradient gain miscalibration using a simple isotropic diffusion phantom. The correction would improve the accuracy of DMRI measurements in the brain and other organs for both DTI and higher order diffusion analysis. In particular, it would allow calibration of MRI systems, improving data harmony in multicenter studies.

Improved reproducibility of diffusion MRI of the human brain with a four-way blip-up and down phase-encoding acquisition approach.

PURPOSE: To assess the effects of blip-up and -down echo planar imaging (EPI) acquisition designs, with different choices of phase-encoding directions (PEDs) on the reproducibility of diffusion MRI (dMRI)-derived metrics in the human brain. METHODS: Diffusion MRI data in seven subjects were acquired five times, each with five different protocols. The base design included 64 diffusion directions acquired with anterior-posterior (AP) PED, the first and second protocols added reverse phase-encoded b=0s/mm2 posterior-anterior (PA) PED images. The third one included 32 directions all with PED acquisitions with opposite polarity (AP and PA). The fourth protocol, also with 32 unique directions used four PEDs (AP, PA, right-left (RL), and left-right (LR)). The scan time was virtually identical for all protocols. The variability of diffusion MRI metrics for each subject and each protocol was computed across the different sessions. RESULTS: The highest reproducibility for all dMRI metrics was obtained with protocol four (AP/PA-RL/LR, ie, four-way PED). Protocols that used only b=0s/mm2 for distortion correction, which are the most widely used designs, had the lowest reproducibility. CONCLUSIONS: An acquisition design with four PEDs, including all DWIs in addition to b=0s/mm2 images should be used to achieve high reproducibility in diffusion MRI studies.

Limits to anatomical accuracy of diffusion tractography using modern approaches.

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 tractography's 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.

Evaluating corrections for Eddy-currents and other EPI distortions in diffusion MRI: methodology and a dataset for benchmarking.

PURPOSE: To propose a methodology for assessment of algorithms that correct distortions due to motion, eddy-currents, and echo planar imaging in diffusion weighted images (DWIs). METHODS: The proposed method evaluates correction performance by measuring variability across datasets of the same object acquired with images having distortions in different directions, thereby overcoming the unavailability of ground-truth, undistorted DWIs. A comprehensive diffusion MRI dataset, collected using a suitable experimental design, is made available to the scientific community, consisting of three DWI shells (Bmax = 5000 s/mm2 ), 30 gradient directions, a replicate set of antipodal gradient directions, four phase-encoding directions, and three different head orientations. The proposed methodology was tested using the TORTOISE diffusion MRI processing pipeline. RESULTS: The median variability of the original distorted data was 123% higher for DWIs, 100-168% higher for tensor-derived metrics and 28-111% higher for MAPMRI metrics, than in the corrected versions. EPI distortions induced substantial variability, nearly comparable to the contribution of eddy-current distortions. CONCLUSIONS: The dataset and the evaluation strategy proposed herein enable quantitative comparison of different methods for correction of distortions due to motion, eddy-currents, and other EPI distortions, and can be useful in benchmarking newly developed algorithms.

Tensor-based morphometry using scalar and directional information of diffusion tensor MRI data (DTBM): Application to hereditary spastic paraplegia.

Tensor-based morphometry (TBM) performed using T1-weighted images (T1WIs) is a well-established method for analyzing local morphological changes occurring in the brain due to normal aging and disease. However, in white matter regions that appear homogeneous on T1WIs, T1W-TBM may be inadequate for detecting changes that affect specific pathways. In these regions, diffusion tensor MRI (DTI) can identify white matter pathways on the basis of their different anisotropy and orientation. In this study, we propose performing TBM using deformation fields constructed using all scalar and directional information provided by the diffusion tensor (DTBM) with the goal of increasing sensitivity in detecting morphological abnormalities of specific white matter pathways. Previously, mostly fractional anisotropy (FA) has been used to drive registration in diffusion MRI-based TBM (FA-TBM). However, FA does not have the directional information that the tensors contain, therefore, the registration based on tensors provides better alignment of brain structures and better localization of volume change. We compare our DTBM method to both T1W-TBM and FA-TBM in investigating differences in brain morphology between patients with complicated hereditary spastic paraplegia of type 11 (SPG11) and a group of healthy controls. Effect size maps of T1W-TBM of SPG11 patients showed diffuse atrophy of white matter. However, DTBM indicated that atrophy was more localized, predominantly affecting several long-range pathways. The results of our study suggest that DTBM could be a powerful tool for detecting morphological changes of specific white matter pathways in normal brain development and aging, as well as in degenerative disorders.

Analysis of the effects of noise, DWI sampling, and value of assumed parameters in diffusion MRI models.

PURPOSE: This study was a systematic evaluation across different and prominent diffusion MRI models to better understand the ways in which scalar metrics are influenced by experimental factors, including experimental design (diffusion-weighted imaging [DWI] sampling) and noise. METHODS: Four diffusion MRI models-diffusion tensor imaging (DTI), diffusion kurtosis imaging (DKI), mean apparent propagator MRI (MAP-MRI), and neurite orientation dispersion and density imaging (NODDI)-were evaluated by comparing maps and histogram values of the scalar metrics generated using DWI datasets obtained in fixed mouse brain with different noise levels and DWI sampling complexity. Additionally, models were fit with different input parameters or constraints to examine the consequences of model fitting procedures. RESULTS: Experimental factors affected all models and metrics to varying degrees. Model complexity influenced sensitivity to DWI sampling and noise, especially for metrics reporting non-Gaussian information. DKI metrics were highly susceptible to noise and experimental design. The influence of fixed parameter selection for the NODDI model was found to be considerable, as was the impact of initial tensor fitting in the MAP-MRI model. CONCLUSION: Across DTI, DKI, MAP-MRI, and NODDI, a wide range of dependence on experimental factors was observed that elucidate principles and practical implications for advanced diffusion MRI. Magn Reson Med 78:1767-1780, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Anatomical accuracy of brain connections derived from diffusion MRI tractography is inherently limited.

Tractography based on diffusion-weighted MRI (DWI) is widely used for mapping the structural connections of the human brain. Its accuracy is known to be limited by technical factors affecting in vivo data acquisition, such as noise, artifacts, and data undersampling resulting from scan time constraints. It generally is assumed that improvements in data quality and implementation of sophisticated tractography methods will lead to increasingly accurate maps of human anatomical connections. However, assessing the anatomical accuracy of DWI tractography is difficult because of the lack of independent knowledge of the true anatomical connections in humans. Here we investigate the future prospects of DWI-based connectional imaging by applying advanced tractography methods to an ex vivo DWI dataset of the macaque brain. The results of different tractography methods were compared with maps of known axonal projections from previous tracer studies in the macaque. Despite the exceptional quality of the DWI data, none of the methods demonstrated high anatomical accuracy. The methods that showed the highest sensitivity showed the lowest specificity, and vice versa. Additionally, anatomical accuracy was highly dependent upon parameters of the tractography algorithm, with different optimal values for mapping different pathways. These results suggest that there is an inherent limitation in determining long-range anatomical projections based on voxel-averaged estimates of local fiber orientation obtained from DWI data that is unlikely to be overcome by improvements in data acquisition and analysis alone.

DR-TAMAS: Diffeomorphic Registration for Tensor Accurate Alignment of Anatomical Structures.

In this work, we propose DR-TAMAS (Diffeomorphic Registration for Tensor Accurate alignMent of Anatomical Structures), a novel framework for intersubject registration of Diffusion Tensor Imaging (DTI) data sets. This framework is optimized for brain data and its main goal is to achieve an accurate alignment of all brain structures, including white matter (WM), gray matter (GM), and spaces containing cerebrospinal fluid (CSF). Currently most DTI-based spatial normalization algorithms emphasize alignment of anisotropic structures. While some diffusion-derived metrics, such as diffusion anisotropy and tensor eigenvector orientation, are highly informative for proper alignment of WM, other tensor metrics such as the trace or mean diffusivity (MD) are fundamental for a proper alignment of GM and CSF boundaries. Moreover, it is desirable to include information from structural MRI data, e.g., T1-weighted or T2-weighted images, which are usually available together with the diffusion data. The fundamental property of DR-TAMAS is to achieve global anatomical accuracy by incorporating in its cost function the most informative metrics locally. Another important feature of DR-TAMAS is a symmetric time-varying velocity-based transformation model, which enables it to account for potentially large anatomical variability in healthy subjects and patients. The performance of DR-TAMAS is evaluated with several data sets and compared with other widely-used diffeomorphic image registration techniques employing both full tensor information and/or DTI-derived scalar maps. Our results show that the proposed method has excellent overall performance in the entire brain, while being equivalent to the best existing methods in WM.

Clinical feasibility of using mean apparent propagator (MAP) MRI to characterize brain tissue microstructure.

Diffusion tensor imaging (DTI) is the most widely used method for characterizing noninvasively structural and architectural features of brain tissues. However, the assumption of a Gaussian spin displacement distribution intrinsic to DTI weakens its ability to describe intricate tissue microanatomy. Consequently, the biological interpretation of microstructural parameters, such as fractional anisotropy or mean diffusivity, is often equivocal. We evaluate the clinical feasibility of assessing brain tissue microstructure with mean apparent propagator (MAP) MRI, a powerful analytical framework that efficiently measures the probability density function (PDF) of spin displacements and quantifies useful metrics of this PDF indicative of diffusion in complex microstructure (e.g., restrictions, multiple compartments). Rotation invariant and scalar parameters computed from the MAP show consistent variation across neuroanatomical brain regions and increased ability to differentiate tissues with distinct structural and architectural features compared with DTI-derived parameters. The return-to-origin probability (RTOP) appears to reflect cellularity and restrictions better than MD, while the non-Gaussianity (NG) measures diffusion heterogeneity by comprehensively quantifying the deviation between the spin displacement PDF and its Gaussian approximation. Both RTOP and NG can be decomposed in the local anatomical frame for reference determined by the orientation of the diffusion tensor and reveal additional information complementary to DTI. The propagator anisotropy (PA) shows high tissue contrast even in deep brain nuclei and cortical gray matter and is more uniform in white matter than the FA, which drops significantly in regions containing crossing fibers. Orientational profiles of the propagator computed analytically from the MAP MRI series coefficients allow separation of different fiber populations in regions of crossing white matter pathways, which in turn improves our ability to perform whole-brain fiber tractography. Reconstructions from subsampled data sets suggest that MAP MRI parameters can be computed from a relatively small number of DWIs acquired with high b-value and good signal-to-noise ratio in clinically achievable scan durations of less than 10min. The neuroanatomical consistency across healthy subjects and reproducibility in test-retest experiments of MAP MRI microstructural parameters further substantiate the robustness and clinical feasibility of this technique. The MAP MRI metrics could potentially provide more sensitive clinical biomarkers with increased pathophysiological specificity compared to microstructural measures derived using conventional diffusion MRI techniques.

Tract Orientation and Angular Dispersion Deviation Indicator (TOADDI): A framework for single-subject analysis in diffusion tensor imaging.

The purpose of this work is to develop a framework for single-subject analysis of diffusion tensor imaging (DTI) data. This framework is termed Tract Orientation and Angular Dispersion Deviation Indicator (TOADDI) because it is capable of testing whether an individual tract as represented by the major eigenvector of the diffusion tensor and its corresponding angular dispersion are significantly different from a group of tracts on a voxel-by-voxel basis. This work develops two complementary statistical tests based on the elliptical cone of uncertainty, which is a model of uncertainty or dispersion of the major eigenvector of the diffusion tensor. The orientation deviation test examines whether the major eigenvector from a single subject is within the average elliptical cone of uncertainty formed by a collection of elliptical cones of uncertainty. The shape deviation test is based on the two-tailed Wilcoxon-Mann-Whitney two-sample test between the normalized shape measures (area and circumference) of the elliptical cones of uncertainty of the single subject against a group of controls. The False Discovery Rate (FDR) and False Non-discovery Rate (FNR) were incorporated in the orientation deviation test. The shape deviation test uses FDR only. TOADDI was found to be numerically accurate and statistically effective. Clinical data from two Traumatic Brain Injury (TBI) patients and one non-TBI subject were tested against the data obtained from a group of 45 non-TBI controls to illustrate the application of the proposed framework in single-subject analysis. The frontal portion of the superior longitudinal fasciculus seemed to be implicated in both tests (orientation and shape) as significantly different from that of the control group. The TBI patients and the single non-TBI subject were well separated under the shape deviation test at the chosen FDR level of 0.0005. TOADDI is a simple but novel geometrically based statistical framework for analyzing DTI data. TOADDI may be found useful in single-subject, graph-theoretic and group analyses of DTI data or DTI-based tractography techniques.

The diffusion tensor imaging (DTI) component of the NIH MRI study of normal brain development (PedsDTI).

The NIH MRI Study of normal brain development sought to characterize typical brain development in a population of infants, toddlers, children and adolescents/young adults, covering the socio-economic and ethnic diversity of the population of the United States. The study began in 1999 with data collection commencing in 2001 and concluding in 2007. The study was designed with the final goal of providing a controlled-access database; open to qualified researchers and clinicians, which could serve as a powerful tool for elucidating typical brain development and identifying deviations associated with brain-based disorders and diseases, and as a resource for developing computational methods and image processing tools. This paper focuses on the DTI component of the NIH MRI study of normal brain development. In this work, we describe the DTI data acquisition protocols, data processing steps, quality assessment procedures, and data included in the database, along with database access requirements. For more details, visit http://www.pediatricmri.nih.gov. This longitudinal DTI dataset includes raw and processed diffusion data from 498 low resolution (3 mm) DTI datasets from 274 unique subjects, and 193 high resolution (2.5 mm) DTI datasets from 152 unique subjects. Subjects range in age from 10 days (from date of birth) through 22 years. Additionally, a set of age-specific DTI templates are included. This forms one component of the larger NIH MRI study of normal brain development which also includes T1-, T2-, proton density-weighted, and proton magnetic resonance spectroscopy (MRS) imaging data, and demographic, clinical and behavioral data.

DR-BUDDI (Diffeomorphic Registration for Blip-Up blip-Down Diffusion Imaging) method for correcting echo planar imaging distortions.

We propose an echo planar imaging (EPI) distortion correction method (DR-BUDDI), specialized for diffusion MRI, which uses data acquired twice with reversed phase encoding directions, often referred to as blip-up blip-down acquisitions. DR-BUDDI can incorporate information from an undistorted structural MRI and also use diffusion-weighted images (DWI) to guide the registration, improving the quality of the registration in the presence of large deformations and in white matter regions. DR-BUDDI does not require the transformations for correcting blip-up and blip-down images to be the exact inverse of each other. Imposing the theoretical "blip-up blip-down distortion symmetry" may not be appropriate in the presence of common clinical scanning artifacts such as motion, ghosting, Gibbs ringing, vibrations, and low signal-to-noise. The performance of DR-BUDDI is evaluated with several data sets and compared to other existing blip-up blip-down correction approaches. The proposed method is robust and generally outperforms existing approaches. The inclusion of the DWIs in the correction process proves to be important to obtain a reliable correction of distortions in the brain stem. Methods that do not use DWIs may produce a visually appealing correction of the non-diffusion weighted images, but the directionally encoded color maps computed from the tensor reveal an abnormal anatomy of the white matter pathways.

Investigation of vibration-induced artifact in clinical diffusion-weighted imaging of pediatric subjects.

It has been reported that mechanical vibrations of the magnetic resonance imaging scanner could produce spurious signal dropouts in diffusion-weighted images resulting in artifactual anisotropy in certain regions of the brain with red appearance in the Directionally Encoded Color maps. We performed a review of the frequency of this artifact across pediatric studies, noting differences by scanner manufacturer, acquisition protocol, as well as weight and position of the subject. We also evaluated the ability of automated and quantitative methods to detect this artifact. We found that the artifact may be present in over 50% of data in certain protocols and is not limited to one scanner manufacturer. While a specific scanner had the highest incidence, low body weight and positioning were also associated with appearance of the artifact for both scanner types evaluated, making children potentially more susceptible than adults. Visual inspection remains the best method for artifact identification. Software for automated detection showed very low sensitivity (10%). The artifact may present inconsistently in longitudinal studies. We discuss a published case report that has been widely cited and used as evidence to set policy about diagnostic criteria for determining vegetative state. That report attributed longitudinal changes in anisotropy to white matter plasticity without considering the possibility that the changes were caused by this artifact. Our study underscores the need to check for the presence of this artifact in clinical studies, analyzes circumstances for when it may be more likely to occur, and suggests simple strategies to identify and potentially avoid its effects.

'Pscore' - A Novel Percentile-Based Metric to Accurately Assess Individual Deviations in Non-Gaussian Distributions of Quantitative MRI Metrics.

BACKGROUND: Quantitative MRI metrics could be used in personalized medicine to assess individuals against normative distributions. Conventional Zscore analysis is inadequate in the presence of non-Gaussian distributions. Therefore, if quantitative MRI metrics deviate from normality, an alternative is needed. PURPOSE: To confirm non-Gaussianity of diffusion MRI (dMRI) metrics on a publicly available dataset, and to propose a novel percentile-based method, 'Pscore' to address this issue. STUDY TYPE: Retrospective cohort. POPULATION: 961 healthy young-adults (age:22-35 years, Females:53%) from the Human Connectome Project. FIELD STRENGTH/SEQUENCE: 3-T, spin-echo diffusion echo-planar imaging, T1-weighted: MPRAGE. ASSESSMENT: The dMRI data were preprocessed using the TORTOISE pipeline. Forty-eight regions of interest (ROIs) from the JHU-atlas were redrawn on a study-specific diffusion tensor (DT) template and average values were computed from various DT and mean apparent propagator (MAP) metrics. For each ROI, percentile ranks across participants were computed to generate 'Pscores'- which normalized the difference between the median and a participant's value with the corresponding difference between the median and the 5th/95th percentile values. STATISTICAL TESTS: ROI-wise distributions were assessed using Log transformations, Zscore, and the 'Pscore' methods. The percentages of extreme values above-95th and below-5th percentile boundaries (PEV<5(%),PEV<5(%)) were also assessed in the overall white matter. Bootstrapping was performed to test the reliability of Pscores in small samples (n=100) using 100 iterations. RESULTS: The dMRI metric distributions were systematically non-Gaussian, including positively skewed (e.g., mean and radial distributions PEV>95≠5%,PEV<5≠5% whereas 'Pscore' distributions were symmetric and balanced PEV>95=PEV<5=5%; even for small bootstrapped samples (average PEV>95¯=PEV<5¯=5±0%SD). DATA CONCLUSION: The inherent skewness observed for dMRI metrics may preclude the use of conventional Zscore analysis. The proposed 'Pscore' method may help estimating individual deviations more accurately in skewed normative data, even from small datasets.

Quantifying brain development in the HEALthy Brain and Child Development (HBCD) Study: The magnetic resonance imaging and spectroscopy protocol.

The HEALthy Brain and Child Development (HBCD) Study, a multi-site prospective longitudinal cohort study, will examine human brain, cognitive, behavioral, social, and emotional development beginning prenatally and planned through early childhood. The acquisition of multimodal magnetic resonance-based brain development data is central to the study's core protocol. However, application of Magnetic Resonance Imaging (MRI) methods in this population is complicated by technical challenges and difficulties of imaging in early life. Overcoming these challenges requires an innovative and harmonized approach, combining age-appropriate acquisition protocols together with specialized pediatric neuroimaging strategies. The HBCD MRI Working Group aimed to establish a core acquisition protocol for all 27 HBCD Study recruitment sites to measure brain structure, function, microstructure, and metabolites. Acquisition parameters of individual modalities have been matched across MRI scanner platforms for harmonized acquisitions and state-of-the-art technologies are employed to enable faster and motion-robust imaging. Here, we provide an overview of the HBCD MRI protocol, including decisions of individual modalities and preliminary data. The result will be an unparalleled resource for examining early neurodevelopment which enables the larger scientific community to assess normative trajectories from birth through childhood and to examine the genetic, biological, and environmental factors that help shape the developing brain.

Mean apparent propagator (MAP) MRI: a novel diffusion imaging method for mapping tissue microstructure.

Diffusion-weighted magnetic resonance (MR) signals reflect information about underlying tissue microstructure and cytoarchitecture. We propose a quantitative, efficient, and robust mathematical and physical framework for representing diffusion-weighted MR imaging (MRI) data obtained in "q-space," and the corresponding "mean apparent propagator (MAP)" describing molecular displacements in "r-space." We also define and map novel quantitative descriptors of diffusion that can be computed robustly using this MAP-MRI framework. We describe efficient analytical representation of the three-dimensional q-space MR signal in a series expansion of basis functions that accurately describes diffusion in many complex geometries. The lowest order term in this expansion contains a diffusion tensor that characterizes the Gaussian displacement distribution, equivalent to diffusion tensor MRI (DTI). Inclusion of higher order terms enables the reconstruction of the true average propagator whose projection onto the unit "displacement" sphere provides an orientational distribution function (ODF) that contains only the orientational dependence of the diffusion process. The representation characterizes novel features of diffusion anisotropy and the non-Gaussian character of the three-dimensional diffusion process. Other important measures this representation provides include the return-to-the-origin probability (RTOP), and its variants for diffusion in one- and two-dimensions-the return-to-the-plane probability (RTPP), and the return-to-the-axis probability (RTAP), respectively. These zero net displacement probabilities measure the mean compartment (pore) volume and cross-sectional area in distributions of isolated pores irrespective of the pore shape. MAP-MRI represents a new comprehensive framework to model the three-dimensional q-space signal and transform it into diffusion propagators. Experiments on an excised marmoset brain specimen demonstrate that MAP-MRI provides several novel, quantifiable parameters that capture previously obscured intrinsic features of nervous tissue microstructure. This should prove helpful for investigating the functional organization of normal and pathologic nervous tissue.

In vivo disentanglement of diffusion frequency-dependence, tensor shape, and relaxation using multidimensional MRI.

Diffusion MRI with free gradient waveforms, combined with simultaneous relaxation encoding, referred to as multidimensional MRI (MD-MRI), offers microstructural specificity in complex biological tissue. This approach delivers intravoxel information about the microstructure, local chemical composition, and importantly, how these properties are coupled within heterogeneous tissue containing multiple microenvironments. Recent theoretical advances incorporated diffusion time dependency and integrated MD-MRI with concepts from oscillating gradients. This framework probes the diffusion frequency, ω, in addition to the diffusion tensor, D, and relaxation, R1, R2, correlations. A D(ω)-R1-R2 clinical imaging protocol was then introduced, with limited brain coverage and 3 mm3 voxel size, which hinder brain segmentation and future cohort studies. In this study, we introduce an efficient, sparse in vivo MD-MRI acquisition protocol providing whole brain coverage at 2 mm3 voxel size. We demonstrate its feasibility and robustness using a well-defined phantom and repeated scans of five healthy individuals. Additionally, we test different denoising strategies to address the sparse nature of this protocol, and show that efficient MD-MRI encoding design demands a nuanced denoising approach. The MD-MRI framework provides rich information that allows resolving the diffusion frequency dependence into intravoxel components based on their D(ω)-R1-R2 distribution, enabling the creation of microstructure-specific maps in the human brain. Our results encourage the broader adoption and use of this new imaging approach for characterizing healthy and pathological tissues.