Characterizing brain mechanics through 7 tesla magnetic resonance elastography.

Magnetic resonance elastography (MRE) is a non-invasive method for determining the mechanical response of tissues using applied harmonic deformation and motion-sensitive MRI. MRE studies of the human brain are typically performed at conventional field strengths, with a few attempts at the ultra-high field strength, 7T, reporting increased spatial resolution with partial brain coverage. Achieving high-resolution human brain scans using 7T MRE presents unique challenges of decreased octahedral shear strain-based signal-to-noise ratio (OSS-SNR) and lower shear wave motion sensitivity. In this study, we establish high resolution MRE at 7T with a custom 2D multi-slice single-shot spin-echo echo-planar imaging sequence, using the Gadgetron advanced image reconstruction framework, applying Marchenko-Pastur Principal component analysis denoising, and using nonlinear viscoelastic inversion. These techniques allowed us to calculate the viscoelastic properties of the whole human brain at 1.1 mm isotropic imaging resolution with high OSS-SNR and repeatability. Using phantom models and 7T MRE data of eighteen healthy volunteers, we demonstrate the robustness and accuracy of our method at high-resolution while quantifying the feasible tradeoff between resolution, OSS-SNR, and scan time. Using these post-processing techniques, we significantly increased OSS-SNR at 1.1 mm resolution with whole-brain coverage by approximately 4-fold and generated elastograms with high anatomical detail. Performing high-resolution MRE at 7T on the human brain can provide information on different substructures within brain tissue based on their mechanical properties, which can then be used to diagnose pathologies (e.g. Alzheimer's disease), indicate disease progression, or better investigate neurodegeneration effects or other relevant brain disorders,in vivo.

A site-wise reliability analysis of the ABCD diffusion fractional anisotropy and cortical thickness: impact of scanner platforms.

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 non-biological 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 to 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.

brainlife.io: a decentralized and open-source cloud platform to support neuroscience research.

Neuroscience is advancing standardization and tool development to support rigor and transparency. Consequently, data pipeline complexity has increased, hindering FAIR (findable, accessible, interoperable and reusable) access. brainlife.io was developed to democratize neuroimaging research. The platform provides data standardization, management, visualization and processing and automatically tracks the provenance history of thousands of data objects. Here, brainlife.io is described and evaluated for validity, reliability, reproducibility, replicability and scientific utility using four data modalities and 3,200 participants.

Assessment of Precision and Accuracy of Brain White Matter Microstructure using Combined Diffusion MRI and Relaxometry.

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, 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 2-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 free water fractions while also capturing compartmental T2 relaxation times. Jointly measuring diffusion and relaxation offers increased sensitivity and specificity to microstructure parameters in brain white matter with voxelwise coefficients of variation below 10%.

Optimization and Validation of the DESIGNER dMRI preprocessing pipeline in white matter aging.

Various diffusion MRI (dMRI) preprocessing pipelines are currently available to yield more accurate diffusion parameters. Here, we evaluated accuracy and robustness of the optimized Diffusion parameter EStImation with Gibbs and NoisE Removal (DESIGNER) pipeline in a large clinical dMRI dataset and using ground truth phantoms. DESIGNER has been modified to improve denoising and target Gibbs ringing for partial Fourier acquisitions. We compared the revisited DESIGNER (Dv2) (including denoising, Gibbs removal, correction for motion, EPI distortion, and eddy currents) against the original DESIGNER (Dv1) pipeline, minimal preprocessing (including correction for motion, EPI distortion, and eddy currents only), and no preprocessing on a large clinical dMRI dataset of 524 control subjects with ages between 25 and 75 years old. We evaluated the effect of specific processing steps on age correlations in white matter with DTI and DKI metrics. We also evaluated the added effect of minimal Gaussian smoothing to deal with noise and to reduce outliers in parameter maps compared to DESIGNER (Dv2)'s noise removal method. Moreover, DESIGNER (Dv2)'s updated noise and Gibbs removal methods were assessed using ground truth dMRI phantom to evaluate accuracy. Results show age correlation in white matter with DTI and DKI metrics were affected by the preprocessing pipeline, causing systematic differences in absolute parameter values and loss or gain of statistical significance. Both in clinical dMRI and ground truth phantoms, DESIGNER (Dv2) pipeline resulted in the smallest number of outlier voxels and improved accuracy in DTI and DKI metrics as noise was reduced and Gibbs removal was improved. Thus, DESIGNER (Dv2) provides more accurate and robust DTI and DKI parameter maps as compared to no preprocessing or minimal preprocessing.

White Matter Cerebrovascular Reactivity: Effects of Microangiopathy and Proximal Occlusions on the Dynamic BOLD Response.

Introduction: Cerebral microangiopathy often manifests as white matter hyperintensities (WMH) on T2-weighted MR images and is associated with elevated stroke risk. Large vessel steno-occlusive disease (SOD) is also independently associated with stroke risk, however, the interaction of microangiopathy and SOD is not well understood. Cerebrovascular reactivity (CVR) describes the capacity of cerebral circulation to adapt to changes in perfusion pressure and neurovascular demand, and its impairment portends future infarctions. CVR can be measured with blood oxygen level dependent (BOLD) imaging following acetazolamide stimulus (ACZ-BOLD). We studied CVR differences between WMH and normal-appearing white matter (NAWM) in patients with chronic SOD, hypothesizing additive influences upon CVR measured by novel, fully dynamic CVR maxima ( CVR max ). Methods: A cross sectional study was conducted to measure per-voxel, per-TR maximal CVR ( CVR max ) using a custom computational pipeline in 23 subjects with angiographically-proven unilateral SOD. WMH and NAWM masks were applied to CVR max maps. White matter was subclassified with respect to the SOD-affected hemisphere, including: i. contralateral NAWM; ii. contralateral WMH iii. ipsilateral NAWM; iv. ipsilateral WMH. CVR max was compared between these groups with a Kruskal-Wallis test followed by a Dunn-Sidak post-hoc test for multiple comparisons. Results: 19 subjects (age 50±12 years, 53% female) undergoing 25 examinations met criteria. WMH volume was asymmetric in 16/19 subjects with 13/16 exhibiting higher volumes ipsilateral to SOD. Pairwise comparisons of CVR max between groups was significant with ipsilateral WMH CVR max lower than contralateral NAWM (p=0.015) and contralateral WMH (p=0.003) when comparing in-subject medians and lower than all groups when comparing pooled voxelwise values across all subjects (p<0.0001). No significant relationship between WMH lesion size and CVR max was detected. Conclusion: Our results suggest additive effects of microvascular and macrovascular disease upon white matter CVR, but with greater overall effects relating to macrovascular SOD than to apparent microangiopathy. Dynamic ACZ-BOLD presents a promising path towards a quantitative stroke risk imaging biomarker. BACKGROUND: Cerebral white matter (WM) microangiopathy manifests as sporadic or sometimes confluent high intensity lesions in MR imaging with T2-weighting, and bears known associations with stroke, cognitive disability, depression and other neurological disorders 1-5 . Deep white matter is particularly susceptible to ischemic injury owing to the deprivation of collateral flow between penetrating arterial territories, and hence deep white matter hyperintensities (WMH) may portend future infarctions 6-8 . The pathophysiology of WMH is variable but commonly includes a cascade of microvascular lipohyalinosis and atherosclerosis together with impaired vascular endothelial and neurogliovascular integrity, leading to blood brain barrier dysfunction, interstitial fluid accumulation, and eventually tissue damage 9-14 . Independent of the microcirculation, cervical and intracranial large vessel steno-occlusive disease (SOD) often results from atheromatous disease and is associated with increased risk of stroke owing to thromboembolic phenomena, hypoperfusion, or combinations thereof 15-17 . White matter disease is more common in the affected hemisphere of patients with asymmetric or unilateral SOD, producing both macroscopic WMH detectable by routine structural MRI, as well as microstructural changes and altered structural connectivity detected by advanced diffusion microstructural imaging 18, 19 . An improved understanding of the interaction of microvascular disease (i.e., WMH) and macrovascular steno-occlusion could better inform stroke risk stratification and guide treatment strategies when coexistent. Cerebrovascular reactivity (CVR) is an autoregulatory adaptation characterized by the capacity of the cerebral circulation to respond to physiological or pharmacological vasodilatory stimuli 20-22 . CVR may be heterogeneous and varies across tissue type and pathological states 1, 16 . Alterations in CVR are associated with elevated stroke risk in SOD patients, although white matter CVR, and in particular the CVR profiles of WMH, are only sparsely studied and not fully understood 1, 23-26 . We have previously employed blood oxygen level dependent (BOLD) imaging following a hemodynamic stimulus with acetazolamide (ACZ) in order to measure CVR (i.e. ACZ-BOLD) 21, 27, 28 . Despite the emergence of ACZ-BOLD as a technique for clinical and experimental use, poor signal-to-noise characteristics of the BOLD effect have generally limited its interpretation to coarse, time-averaged assessment of the terminal ACZ response at arbitrarily prescribed delays following ACZ administration (e.g. 10-20 minutes) 29 . More recently, we have introduced a dedicated computational pipeline to overcome historically intractable signal-to-noise ratio (SNR) limitations of BOLD, enabling fully dynamic characterization of the cerebrovascular response, including identification of previously unreported, unsustained or transient CVR maxima ( CVR max ) following hemodynamic provocation 27, 30 . In this study, we compared such dynamic interrogation of true CVR maxima between WMH and normal appearing white matter (NAWM) among patients with chronic, unilateral SOD in order to quantify their interaction and to assess the hypothesized additive effects of angiographically-evident macrovascular stenoses when intersecting microangiopathic WMH.

In vivo evidence of microstructural hypo-connectivity of brain white matter in 22q11.2 deletion syndrome.

22q11.2 deletion syndrome, or 22q11.2DS, is a genetic syndrome associated with high rates of schizophrenia and autism spectrum disorders, in addition to widespread structural and functional abnormalities throughout the brain. Experimental animal models have identified neuronal connectivity deficits, e.g., decreased axonal length and complexity of axonal branching, as a primary mechanism underlying atypical brain development in 22q11.2DS. However, it is still unclear whether deficits in axonal morphology can also be observed in people with 22q11.2DS. Here, we provide an unparalleled in vivo characterization of white matter microstructure in participants with 22q11.2DS (12-15 years) and those undergoing typical development (8-18 years) using a customized magnetic resonance imaging scanner which is sensitive to axonal morphology. A rich array of diffusion MRI metrics are extracted to present microstructural profiles of typical and atypical white matter development, and provide new evidence of connectivity differences in individuals with 22q11.2DS. A recent, large-scale consortium study of 22q11.2DS identified higher diffusion anisotropy and reduced overall diffusion mobility of water as hallmark microstructural alterations of white matter in individuals across a wide age range (6-52 years). We observed similar findings across the white matter tracts included in this study, in addition to identifying deficits in axonal morphology. This, in combination with reduced tract volume measurements, supports the hypothesis that abnormal microstructural connectivity in 22q11.2DS may be mediated by densely packed axons with disproportionately small diameters. Our findings provide insight into the in vivo white matter phenotype of 22q11.2DS, and promote the continued investigation of shared features in neurodevelopmental and psychiatric disorders.

brainlife.io: A decentralized and open source cloud platform to support neuroscience research.

Neuroscience research has expanded dramatically over the past 30 years by advancing standardization and tool development to support rigor and transparency. Consequently, the complexity of the data pipeline has also increased, hindering access to FAIR data analysis to portions of the worldwide research community. brainlife.io was developed to reduce these burdens and democratize modern neuroscience research across institutions and career levels. Using community software and hardware infrastructure, the platform provides open-source data standardization, management, visualization, and processing and simplifies the data pipeline. brainlife.io automatically tracks the provenance history of thousands of data objects, supporting simplicity, efficiency, and transparency in neuroscience research. Here brainlife.io's technology and data services are described and evaluated for validity, reliability, reproducibility, replicability, and scientific utility. Using data from 4 modalities and 3,200 participants, we demonstrate that brainlife.io's services produce outputs that adhere to best practices in modern neuroscience research.

Diaschisis Profiles in the Cerebellar Response to Hemodynamic Stimuli: Insights From Dynamic Measurement of Cerebrovascular Reactivity to Identify Occult and Transient Maxima.

BACKGROUND: Crossed cerebellar diaschisis (CCD) refers to depressions in perfusion and metabolism within the cerebellar hemisphere contralateral to supratentorial disease. Prior investigation into CCD in cerebrovascular reactivity (CVR) has been limited to terminal CVR estimations (CVRend ). We recently have demonstrated the presence of unsustained CVR maxima (CVRmax ) using dynamic CVR analysis, offering a fully dynamic characterization of CVR to hemodynamic stimuli. PURPOSE: To investigate CCD in CVRmax from dynamic blood oxygen level-dependent (BOLD) MRI, by comparison with conventional CVRend estimation. STUDY TYPE: Retrospective. POPULATION: A total of 23 patients (median age: 51 years, 10 females) with unilateral chronic steno-occlusive cerebrovascular disease, without prior knowledge of CCD status. FIELD STRENGTH/SEQUENCE: A 3-T, T1-weighted magnetization-prepared rapid gradient-echo (MPRAGE) and acetazolamide-augmented BOLD imaging performed with a gradient-echo echo-planar imaging (EPI) sequence. ASSESSMENT: A custom denoising pipeline was used to generate BOLD-CVR time signals. CVRend was established using the last minute of the BOLD response relative to the first-minute baseline. Following classification of healthy versus diseased cerebral hemispheres, CVRmax and CVRend were calculated for bilateral cerebral and cerebellar hemispheres. Three independent observers evaluated all data for the presence of CCD. STATISTICAL TESTS: Pearson correlations for comparing CVR across hemispheres, two-proportion Z-tests for comparing CCD prevalence, and Wilcoxon signed-rank tests for comparing median CVR. The level of statistical significance was set at P ≤ 0.05. RESULTS: CCD-related changes were observed on both CVRend and CVRmax maps, with all CCD+ cases identifiable by inspection of either map. Diseased cerebral and contralateral cerebellar hemispheric CVR correlations in CCD+ patients were stronger when using CVRend (r = 0.728) as compared to CVRmax (r = 0.676). CVR correlations between healthy cerebral hemispheres and contralateral cerebellar hemispheres were stronger for CVRmax (r = 0.739) than for CVRend (r = 0.705). DATA CONCLUSION: CCD-related alterations could be observed in CVR examinations. Conventional CVRend may underestimate CVR and could exaggerate CCD. EVIDENCE LEVEL: 4. TECHNICAL EFFICACY: Stage 3.

Acetazolamide-augmented BOLD MRI to Assess Whole-Brain Cerebrovascular Reactivity in Chronic Steno-occlusive Disease Using Principal Component Analysis.

Background Exhaustion of cerebrovascular reactivity (CVR) portends increased stroke risk. Acetazolamide-augmented blood oxygenation level-dependent (BOLD) MRI has been used to estimate CVR, but low signal-to-noise conditions relegate its use to terminal CVR (CVRend) measurements that neglect dynamic features of CVR. Purpose To demonstrate comprehensive characterization of acetazolamide-augmented BOLD MRI response in chronic steno-occlusive disease using a computational framework to precondition signal time courses for dynamic whole-brain CVR analysis. Materials and Methods This study focused on retrospective analysis of consecutive patients with unilateral chronic steno-occlusive disease who underwent acetazolamide-augmented BOLD imaging for recurrent minor stroke or transient ischemic attack at an academic medical center between May 2017 and October 2020. A custom principal component analysis-based denoising pipeline was used to correct spatially varying non-signal-bearing contributions obtained by a local principal component analysis of the MRI time series. Standard voxelwise CVRend maps representing terminal responses were produced and compared with maximal CVR (CVRmax) as isolated from binned (per-repetition time) denoised BOLD time course. A linear mixed-effects model was used to compare CVRmax and CVRend in healthy and diseased hemispheres. Results A total of 23 patients (median age, 51 years; IQR, 42-61, 13 men) who underwent 32 BOLD examinations were included. Processed MRI data showed twofold improvement in signal-to-noise ratio, allowing improved isolation of dynamic characteristics in signal time course for sliding window CVRmax analysis to the level of each BOLD repetition time (approximately 2 seconds). Mean CVRmax was significantly higher than mean CVRend in diseased (5.2% vs 3.8%, P < .01) and healthy (5.5% vs 4.0%, P < .01) hemispheres. Several distinct time-signal signatures were observed, including nonresponsive; delayed/blunted; brisk; and occasionally nonmonotonic time courses with paradoxical features in normal and abnormal tissues (ie, steal and reverse-steal patterns). Conclusion A principal component analysis-based computational framework for analysis of acetazolamide-augmented BOLD imaging can be used to measure unsustained CVRmax through twofold improvements in signal-to-noise ratio. © RSNA, 2023 Supplemental material is available for this article.

Cardiac Phase and Flow Compensation Effects on REnal Flow and Microstructure AnisotroPy MRI in Healthy Human Kidney.

BACKGROUND: Renal diffusion-weighted imaging (DWI) involves microstructure and microcirculation, quantified with diffusion tensor imaging (DTI), intravoxel incoherent motion (IVIM), and hybrid models. A better understanding of their contrast may increase specificity. PURPOSE: To measure modulation of DWI with cardiac phase and flow-compensated (FC) diffusion gradient waveforms. STUDY TYPE: Prospective. POPULATION: Six healthy volunteers (ages: 22-48 years, five females), water phantom. FIELD STRENGTH/SEQUENCE: 3-T, prototype DWI sequence with 2D echo-planar imaging, and bipolar (BP) or FC gradients. 2D Half-Fourier Single-shot Turbo-spin-Echo (HASTE). Multiple-phase 2D spoiled gradient-echo phase contrast (PC) MRI. ASSESSMENT: BP and FC water signal decays were qualitatively compared. Renal arteries and velocities were visualized on PC-MRI. Systolic (peak velocity), diastolic (end stable velocity), and pre-systolic (before peak velocity) phases were identified. Following mutual information-based retrospective self-registration of DWI within each kidney, and Marchenko-Pastur Principal Component Analysis (MPPCA) denoising, combined IVIM-DTI analysis estimated mean diffusivity (MD), fractional anisotropy (FA), and eigenvalues (λi) from tissue diffusivity (Dt ), perfusion fraction (fp ), and pseudodiffusivity (Dp , Dp,axial , Dp,radial ), for each tissue (cortex/medulla, segmented on b0/FA respectively), phase, and waveform (BP, FC). Monte Carlo water diffusion simulations aided data interpretation. STATISTICAL TESTS: Mixed model regression probed differences between tissue types and pulse sequences. Univariate general linear model analysis probed variations among cardiac phases. Spearman correlations were measured between diffusion metrics and renal artery velocities. Statistical significance level was set at P < 0.05. RESULTS: Water BP and FC signal decays showed no differences. Significant pulse sequence dependence occurred for λ1 , λ3 , FA, Dp , fp , Dp,axial , Dp,radial in cortex and medulla, and medullary λ2 . Significant cortex/medulla differences occurred with BP for all metrics except MD (systole [P = 0.224]; diastole [P = 0.556]). Significant phase dependence occurred for Dp , Dp,axial , Dp,radial for BP and medullary λ1 , λ2 , λ3 , MD for FC. FA correlated significantly with velocity. Monte Carlo simulations indicated medullary measurements were consistent with a 34 µm tubule diameter. DATA CONCLUSION: Cardiac gating and flow compensation modulate of measurements of renal diffusion. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY STAGE: 2.

Reproducibility of the Standard Model of diffusion in white matter on clinical MRI systems.

Estimating intra- and extra-axonal microstructure parameters, such as volume fractions and diffusivities, has been one of the major efforts in brain microstructure imaging with MRI. The Standard Model (SM) of diffusion in white matter has unified various modeling approaches based on impermeable narrow cylinders embedded in locally anisotropic extra-axonal space. However, estimating the SM parameters from a set of conventional diffusion MRI (dMRI) measurements is ill-conditioned. Multidimensional dMRI helps resolve the estimation degeneracies, but there remains a need for clinically feasible acquisitions that yield robust parameter maps. Here we find optimal multidimensional protocols by minimizing the mean-squared error of machine learning-based SM parameter estimates for two 3T scanners with corresponding gradient strengths of 40and80mT/m. We assess intra-scanner and inter-scanner repeatability for 15-minute optimal protocols by scanning 20 healthy volunteers twice on both scanners. The coefficients of variation all SM parameters except free water fraction are ≲10% voxelwise and 1-4% for their region-averaged values. As the achieved SM reproducibility outcomes are similar to those of conventional diffusion tensor imaging, our results enable robust in vivo mapping of white matter microstructure in neuroscience research and in the clinic.

Improved diffusion parameter estimation by incorporating T2 relaxation properties into the DKI-FWE model.

The free water elimination (FWE) model and its kurtosis variant (DKI-FWE) can separate tissue and free water signal contributions, thus providing tissue-specific diffusional information. However, a downside of these models is that the associated parameter estimation problem is ill-conditioned, necessitating the use of advanced estimation techniques that can potentially bias the parameter estimates. In this work, we propose the T2-DKI-FWE model that exploits the T2 relaxation properties of both compartments, thereby better conditioning the parameter estimation problem and providing, at the same time, an additional potential biomarker (the T2 of tissue). In our approach, the T2 of tissue is estimated as an unknown parameter, whereas the T2 of free water is assumed known a priori and fixed to a literature value (1573 ms). First, the error propagation of an erroneous assumption on the T2 of free water is studied. Next, the improved conditioning of T2-DKI-FWE compared to DKI-FWE is illustrated using the Cramér-Rao lower bound matrix. Finally, the performance of the T2-DKI-FWE model is compared to that of the DKI-FWE and T2-DKI models on both simulated and real datasets. The error due to a biased approximation of the T2 of free water was found to be relatively small in various diffusion metrics and for a broad range of erroneous assumptions on its underlying ground truth value. Compared to DKI-FWE, using the T2-DKI-FWE model is beneficial for the identifiability of the model parameters. Our results suggest that the T2-DKI-FWE model can achieve precise and accurate diffusion parameter estimates, through effective reduction of free water partial volume effects and by using a standard nonlinear least squares approach. In conclusion, incorporating T2 relaxation properties into the DKI-FWE model improves the conditioning of the model fitting, while only requiring an acquisition scheme with at least two different echo times.

What's new and what's next in diffusion MRI preprocessing.

Diffusion MRI (dMRI) provides invaluable information for the study of tissue microstructure and brain connectivity, but suffers from a range of imaging artifacts that greatly challenge the analysis of results and their interpretability if not appropriately accounted for. This review will cover dMRI artifacts and preprocessing steps, some of which have not typically been considered in existing pipelines or reviews, or have only gained attention in recent years: brain/skull extraction, B-matrix incompatibilities w.r.t the imaging data, signal drift, Gibbs ringing, noise distribution bias, denoising, between- and within-volumes motion, eddy currents, outliers, susceptibility distortions, EPI Nyquist ghosts, gradient deviations, B1 bias fields, and spatial normalization. The focus will be on "what's new" since the notable advances prior to and brought by the Human Connectome Project (HCP), as presented in the predecessing issue on "Mapping the Connectome" in 2013. In addition to the development of novel strategies for dMRI preprocessing, exciting progress has been made in the availability of open source tools and reproducible pipelines, databases and simulation tools for the evaluation of preprocessing steps, and automated quality control frameworks, amongst others. Finally, this review will consider practical considerations and our view on "what's next" in dMRI preprocessing.

Nanostructure-specific X-ray tomography reveals myelin levels, integrity and axon orientations in mouse and human nervous tissue.

Myelin insulates neuronal axons and enables fast signal transmission, constituting a key component of brain development, aging and disease. Yet, myelin-specific imaging of macroscopic samples remains a challenge. Here, we exploit myelin's nanostructural periodicity, and use small-angle X-ray scattering tensor tomography (SAXS-TT) to simultaneously quantify myelin levels, nanostructural integrity and axon orientations in nervous tissue. Proof-of-principle is demonstrated in whole mouse brain, mouse spinal cord and human white and gray matter samples. Outcomes are validated by 2D/3D histology and compared to MRI measurements sensitive to myelin and axon orientations. Specificity to nanostructure is exemplified by concomitantly imaging different myelin types with distinct periodicities. Finally, we illustrate the method's sensitivity towards myelin-related diseases by quantifying myelin alterations in dysmyelinated mouse brain. This non-destructive, stain-free molecular imaging approach enables quantitative studies of myelination within and across samples during development, aging, disease and treatment, and is applicable to other ordered biomolecules or nanostructures.

Toward more robust and reproducible diffusion kurtosis imaging.

PURPOSE: The general utility of diffusion kurtosis imaging (DKI) is challenged by its poor robustness to imaging artifacts and thermal noise that often lead to implausible kurtosis values. THEORY AND METHODS: A robust scalar kurtosis index can be estimated from powder-averaged diffusion-weighted data. We introduce a novel DKI estimator that uses this scalar kurtosis index as a proxy for the mean kurtosis to regularize the fit. RESULTS: The regularized DKI estimator improves the robustness and reproducibility of the kurtosis metrics and results in parameter maps with enhanced quality and contrast. CONCLUSION: Our novel DKI estimator promotes the wider use of DKI in clinical research and potentially diagnostics by improving the reproducibility and precision of DKI fitting and, as such, enabling enhanced visual, quantitative, and statistical analyses of DKI parameters.

The variability of MR axon radii estimates in the human white matter.

The noninvasive quantification of axonal morphology is an exciting avenue for gaining understanding of the function and structure of the central nervous system. Accurate non-invasive mapping of micron-sized axon radii using commonly applied neuroimaging techniques, that is, diffusion-weighted MRI, has been bolstered by recent hardware developments, specifically MR gradient design. Here the whole brain characterization of the effective MR axon radius is presented and the inter- and intra-scanner test-retest repeatability and reproducibility are evaluated to promote the further development of the effective MR axon radius as a neuroimaging biomarker. A coefficient-of-variability of approximately 10% in the voxelwise estimation of the effective MR radius is observed in the test-retest analysis, but it is shown that the performance can be improved fourfold using a customized along-tract analysis.

Improved Task-based Functional MRI Language Mapping in Patients with Brain Tumors through Marchenko-Pastur Principal Component Analysis Denoising.

Background Functional MRI improves preoperative planning in patients with brain tumors, but task-correlated signal intensity changes are only 2%-3% above baseline. This makes accurate functional mapping challenging. Marchenko-Pastur principal component analysis (MP-PCA) provides a novel strategy to separate functional MRI signal from noise without requiring user input or prior data representation. Purpose To determine whether MP-PCA denoising improves activation magnitude for task-based functional MRI language mapping in patients with brain tumors. Materials and Methods In this Health Insurance Portability and Accountability Act-compliant study, MP-PCA performance was first evaluated by using simulated functional MRI data with a known ground truth. Right-handed, left-language-dominant patients with brain tumors who successfully performed verb generation, sentence completion, and finger tapping functional MRI tasks were retrospectively identified between January 2017 and August 2018. On the group level, for each task, histograms of z scores for original and MP-PCA denoised data were extracted from relevant regions and contralateral homologs were seeded by a neuroradiologist blinded to functional MRI findings. Z scores were compared with paired two-sided t tests, and distributions were compared with effect size measurements and the Kolmogorov-Smirnov test. The number of voxels with a z score greater than 3 was used to measure task sensitivity relative to task duration. Results Twenty-three patients (mean age ± standard deviation, 43 years ± 18; 13 women) were evaluated. MP-PCA denoising led to a higher median z score of task-based functional MRI voxel activation in left hemisphere cortical regions for verb generation (from 3.8 ± 1.0 to 4.5 ± 1.4; P < .001), sentence completion (from 3.7 ± 1.0 to 4.3 ± 1.4; P < .001), and finger tapping (from 6.9 ± 2.4 to 7.9 ± 2.9; P < .001). Median z scores did not improve in contralateral homolog regions for verb generation (from -2.7 ± 0.54 to -2.5 ± 0.40; P = .90), sentence completion (from -2.3 ± 0.21 to -2.4 ± 0.37; P = .39), or finger tapping (from -2.3 ± 1.20 to -2.7 ± 1.40; P = .07). Individual functional MRI task durations could be truncated by at least 40% after MP-PCA without degradation of clinically relevant correlations between functional cortex and functional MRI tasks. Conclusion Denoising with Marchenko-Pastur principal component analysis led to higher task correlations in relevant cortical regions during functional MRI language mapping in patients with brain tumors. © RSNA, 2020 Online supplemental material is available for this article.

Cross-scanner and cross-protocol multi-shell diffusion MRI data harmonization: Algorithms and results.

Cross-scanner and cross-protocol variability of diffusion magnetic resonance imaging (dMRI) data are known to be major obstacles in multi-site clinical studies since they limit the ability to aggregate dMRI data and derived measures. Computational algorithms that harmonize the data and minimize such variability are critical to reliably combine datasets acquired from different scanners and/or protocols, thus improving the statistical power and sensitivity of multi-site studies. Different computational approaches have been proposed to harmonize diffusion MRI data or remove scanner-specific differences. To date, these methods have mostly been developed for or evaluated on single b-value diffusion MRI data. In this work, we present the evaluation results of 19 algorithms that are developed to harmonize the cross-scanner and cross-protocol variability of multi-shell diffusion MRI using a benchmark database. The proposed algorithms rely on various signal representation approaches and computational tools, such as rotational invariant spherical harmonics, deep neural networks and hybrid biophysical and statistical approaches. The benchmark database consists of data acquired from the same subjects on two scanners with different maximum gradient strength (80 and 300 â€‹mT/m) and with two protocols. We evaluated the performance of these algorithms for mapping multi-shell diffusion MRI data across scanners and across protocols using several state-of-the-art imaging measures. The results show that data harmonization algorithms can reduce the cross-scanner and cross-protocol variabilities to a similar level as scan-rescan variability using the same scanner and protocol. In particular, the LinearRISH algorithm based on adaptive linear mapping of rotational invariant spherical harmonics features yields the lowest variability for our data in predicting the fractional anisotropy (FA), mean diffusivity (MD), mean kurtosis (MK) and the rotationally invariant spherical harmonic (RISH) features. But other algorithms, such as DIAMOND, SHResNet, DIQT, CMResNet show further improvement in harmonizing the return-to-origin probability (RTOP). The performance of different approaches provides useful guidelines on data harmonization in future multi-site studies.