Investigating apparent differences between standard DKI and axisymmetric DKI and its consequences for biophysical parameter estimates.

PURPOSE: The purpose of the study is to identify differences between axisymmetric diffusion kurtosis imaging (DKI) and standard DKI, their consequences for biophysical parameter estimates, and the protocol choice influence on parameter estimation. METHODS: Noise-free and noisy, synthetic diffusion MRI human brain data is simulated using standard DKI for a standard and the fast "199" acquisition protocol. First the noise-free "baseline" difference between both DKI models is estimated and the influence of fiber complexity is investigated. Noisy data is used to establish the signal-to-noise ratio at which the baseline difference exceeds noise variability. The influence of protocol choices and denoising is investigated. The five axisymmetric DKI tensor metrics (AxTM), the parallel and perpendicular diffusivity and kurtosis and mean of the kurtosis tensor are used to compare both DKI models. Additionally, the baseline difference is also estimated for the five parameters of the WMTI-Watson model. RESULTS: The parallel and perpendicular kurtosis and all of the WMTI-Watson parameters had large baseline differences. Using a Westin or FA mask reduced the number of voxels with large baseline difference, that is, by selecting voxels with less complex fibers. For the noisy data, precision was worsened by the fast "199" protocol but adaptive denoising can help counteract these effects. CONCLUSION: For the diffusivities and mean of the kurtosis tensor, axisymmetric DKI with a standard protocol delivers similar results as standard DKI. Fiber complexity is one main driver of the baseline differences. Using the "199" protocol worsens precision in noisy data but adaptive denoising mitigates these effects.

Axisymmetric diffusion kurtosis imaging with Rician bias correction: A simulation study.

PURPOSE: To compare the estimation accuracy of axisymmetric diffusion kurtosis imaging (DKI) and standard DKI in combination with Rician bias correction (RBC). METHODS: Axisymmetric DKI is more robust against noise-induced variation in the measured signal than standard DKI because of its reduced parameter space. However, its susceptibility to Rician noise bias at low signal-to-noise ratios (SNR) is unknown. Here, we investigate two main questions: first, does RBC improve estimation accuracy of axisymmetric DKI?; second, is estimation accuracy of axisymmetric DKI increased compared to standard DKI? Estimation accuracy was investigated on the five axisymmetric DKI tensor metrics (AxTM): the parallel and perpendicular diffusivity and kurtosis and mean of the kurtosis tensor, using a noise simulation study based on synthetic data of tissues with varying fiber alignment and in-vivo data focusing on white matter. RESULTS: RBC mainly increased accuracy for the parallel AxTM in tissues with highly to moderately aligned fibers. For the perpendicular AxTM, axisymmetric DKI without RBC performed slightly better than with RBC. However, the combination of axisymmetric DKI with RBC was the overall best performing algorithm across all five AxTM in white matter and axisymmetric DKI itself substantially improved accuracy in axisymmetric tissues with low fiber alignment. CONCLUSION: Combining axisymmetric DKI with RBC facilitates accurate DKI parameter estimation at unprecedented low SNRs ( ≈ 15 $$ \approx 15 $$ ) in white matter, possibly making it a valuable tool for neuroscience and clinical research studies where scan time is a limited resource. The tools used here are available in the open-source ACID toolbox for SPM.

Error quantification in multi-parameter mapping facilitates robust estimation and enhanced group level sensitivity.

Multi-Parameter Mapping (MPM) is a comprehensive quantitative neuroimaging protocol that enables estimation of four physical parameters (longitudinal and effective transverse relaxation rates R1 and R2*, proton density PD, and magnetization transfer saturation MTsat) that are sensitive to microstructural tissue properties such as iron and myelin content. Their capability to reveal microstructural brain differences, however, is tightly bound to controlling random noise and artefacts (e.g. caused by head motion) in the signal. Here, we introduced a method to estimate the local error of PD, R1, and MTsat maps that captures both noise and artefacts on a routine basis without requiring additional data. To investigate the method's sensitivity to random noise, we calculated the model-based signal-to-noise ratio (mSNR) and showed in measurements and simulations that it correlated linearly with an experimental raw-image-based SNR map. We found that the mSNR varied with MPM protocols, magnetic field strength (3T vs. 7T) and MPM parameters: it halved from PD to R1 and decreased from PD to MTsat by a factor of 3-4. Exploring the artefact-sensitivity of the error maps, we generated robust MPM parameters using two successive acquisitions of each contrast and the acquisition-specific errors to down-weight erroneous regions. The resulting robust MPM parameters showed reduced variability at the group level as compared to their single-repeat or averaged counterparts. The error and mSNR maps may better inform power-calculations by accounting for local data quality variations across measurements. Code to compute the mSNR maps and robustly combined MPM maps is available in the open-source hMRI toolbox.

Denoising for Improved Parametric MRI of the Kidney: Protocol for Nonlocal Means Filtering.

In order to tackle the challenges caused by the variability in estimated MRI parameters (e.g., T2* and T2) due to low SNR a number of strategies can be followed. One approach is postprocessing of the acquired data with a filter. The basic idea is that MR images possess a local spatial structure that is characterized by equal, or at least similar, noise-free signal values in vicinities of a location. Then, local averaging of the signal reduces the noise component of the signal. In contrast, nonlocal means filtering defines the weights for averaging not only within the local vicinity, bur it compares the image intensities between all voxels to define "nonlocal" weights. Furthermore, it generally compares not only single-voxel intensities but small spatial patches of the data to better account for extended similar patterns. Here we describe how to use an open source NLM filter tool to denoise 2D MR image series of the kidney used for parametric mapping of the relaxation times T2* and T2.This chapter is based upon work from the COST Action PARENCHIMA, a community-driven network funded by the European Cooperation in Science and Technology (COST) program of the European Union, which aims to improve the reproducibility and standardization of renal MRI biomarkers.

Example dataset for the hMRI toolbox.

The hMRI toolbox is an open-source toolbox for the calculation of quantitative MRI parameter maps from a series of weighted imaging data, and optionally additional calibration data. The multi-parameter mapping (MPM) protocol, incorporating calibration data to correct for spatial variation in the scanner's transmit and receive fields, is the most complete protocol that can be handled by the toolbox. Here we present a dataset acquired with such a full MPM protocol, which is made freely available to be used as a tutorial by following instructions provided on the associated toolbox wiki pages, which can be found at http://hMRI.info, and following the theory described in: hMRI - A toolbox for quantitative MRI in neuroscience and clinical research [1].

hMRI - A toolbox for quantitative MRI in neuroscience and clinical research.

Neuroscience and clinical researchers are increasingly interested in quantitative magnetic resonance imaging (qMRI) due to its sensitivity to micro-structural properties of brain tissue such as axon, myelin, iron and water concentration. We introduce the hMRI-toolbox, an open-source, easy-to-use tool available on GitHub, for qMRI data handling and processing, presented together with a tutorial and example dataset. This toolbox allows the estimation of high-quality multi-parameter qMRI maps (longitudinal and effective transverse relaxation rates R1 and R2⋆, proton density PD and magnetisation transfer MT saturation) that can be used for quantitative parameter analysis and accurate delineation of subcortical brain structures. The qMRI maps generated by the toolbox are key input parameters for biophysical models designed to estimate tissue microstructure properties such as the MR g-ratio and to derive standard and novel MRI biomarkers. Thus, the current version of the toolbox is a first step towards in vivo histology using MRI (hMRI) and is being extended further in this direction. Embedded in the Statistical Parametric Mapping (SPM) framework, it benefits from the extensive range of established SPM tools for high-accuracy spatial registration and statistical inferences and can be readily combined with existing SPM toolboxes for estimating diffusion MRI parameter maps. From a user's perspective, the hMRI-toolbox is an efficient, robust and simple framework for investigating qMRI data in neuroscience and clinical research.

Magnetic resonance advection imaging of cerebrovascular pulse dynamics.

We analyze the pulsatile signal component of dynamic echo planar imaging data from the brain by modeling the dependence between local temporal and spatial signal variability. The resulting magnetic resonance advection imaging maps depict the location of major arteries. Color direction maps allow for visualization of the direction of blood vessels. The potential significance of magnetic resonance advection imaging maps is demonstrated on a functional magnetic resonance imaging data set of 19 healthy subjects. A comparison with the here introduced pulse coherence maps, in which the echo planar imaging signal is correlated with a cardiac pulse signal, shows that the magnetic resonance advection imaging approach results in a better spatial definition without the need for a pulse reference. In addition, it is shown that magnetic resonance advection imaging velocities can be estimates of pulse wave velocities if certain requirements are met, which are specified. Although for this application magnetic resonance advection imaging velocities are not quantitative estimates of pulse wave velocities, they clearly depict local pulsatile dynamics. Magnetic resonance advection imaging can be applied to existing dynamic echo planar imaging data sets with sufficient spatiotemporal resolution. It is discussed whether magnetic resonance advection imaging might have the potential to evolve into a biomarker for the health of the cerebrovascular system.

Improving Accuracy and Temporal Resolution of Learning Curve Estimation for within- and across-Session Analysis.

Estimation of learning curves is ubiquitously based on proportions of correct responses within moving trial windows. Thereby, it is tacitly assumed that learning performance is constant within the moving windows, which, however, is often not the case. In the present study we demonstrate that violations of this assumption lead to systematic errors in the analysis of learning curves, and we explored the dependency of these errors on window size, different statistical models, and learning phase. To reduce these errors in the analysis of single-subject data as well as on the population level, we propose adequate statistical methods for the estimation of learning curves and the construction of confidence intervals, trial by trial. Applied to data from an avoidance learning experiment with rodents, these methods revealed performance changes occurring at multiple time scales within and across training sessions which were otherwise obscured in the conventional analysis. Our work shows that the proper assessment of the behavioral dynamics of learning at high temporal resolution can shed new light on specific learning processes, and, thus, allows to refine existing learning concepts. It further disambiguates the interpretation of neurophysiological signal changes recorded during training in relation to learning.

More Specific Signal Detection in Functional Magnetic Resonance Imaging by False Discovery Rate Control for Hierarchically Structured Systems of Hypotheses.

Signal detection in functional magnetic resonance imaging (fMRI) inherently involves the problem of testing a large number of hypotheses. A popular strategy to address this multiplicity is the control of the false discovery rate (FDR). In this work we consider the case where prior knowledge is available to partition the set of all hypotheses into disjoint subsets or families, e. g., by a-priori knowledge on the functionality of certain regions of interest. If the proportion of true null hypotheses differs between families, this structural information can be used to increase statistical power. We propose a two-stage multiple test procedure which first excludes those families from the analysis for which there is no strong evidence for containing true alternatives. We show control of the family-wise error rate at this first stage of testing. Then, at the second stage, we proceed to test the hypotheses within each non-excluded family and obtain asymptotic control of the FDR within each family at this second stage. Our main mathematical result is that this two-stage strategy implies asymptotic control of the FDR with respect to all hypotheses. In simulations we demonstrate the increased power of this new procedure in comparison with established procedures in situations with highly unbalanced families. Finally, we apply the proposed method to simulated and to real fMRI data.

Evidence for early, non-lesional cerebellar damage in patients with multiple sclerosis: DTI measures correlate with disability, atrophy, and disease duration.

BACKGROUND: Common symptoms of multiple sclerosis (MS) such as gait ataxia, poor coordination of the hands, and intention tremor are usually the result of dysfunctionality in the cerebellum. Magnetic resonance imaging (MRI) has frequently failed to detect cerebellar damage in the form of inflammatory lesions in patients presenting with symptoms of cerebellar dysfunction. OBJECTIVE: To detect microstructural cerebellar tissue alterations in early MS patients with a "normal appearing" cerebellum using diffusion tensor imaging (DTI). METHODS: A total of 68 patients with relapsing-remitting MS (RRMS) and without cerebellar lesions and 26 age-matched healthy controls were admitted to high-resolution MRI and DTI to assess microstructure and volume of the cerebellar white matter (CBWM). RESULTS: We found cerebellar fractional anisotropy (FA) and CBWM volume reductions in the group of 68 patients. Interestingly, a subgroup of these patients that was derived by including only patients with early and mild MS (N=23, median age 30 years, median Expanded Disability Status Scale =1.5, median duration 28 months) showed already cerebellar FA but no CBWM volume reductions. FA reductions were correlated with disability, atrophy, and disease duration. CONCLUSION: "Normal appearing" cerebellar white matter can be damaged in a very early stage of RRMS. DTI seems to be a sensitive tool for detecting this hidden cerebellar damage.

Recovery of thalamic microstructural damage after Shiga toxin 2-associated hemolytic-uremic syndrome.

INTRODUCTION: The underlying pathophysiology of neurological complications in patients with hemolytic-uremic syndrome (HUS) remains unclear. It was recently attributed to a direct cytotoxic effect of Shiga toxin 2 (Stx2) in the thalamus. Conventional MRI of patients with Stx2-caused HUS revealed - despite severe neurological symptoms - only mild alterations if any, mostly in the thalamus. Against this background, we questioned: Does diffusion tensor imaging (DTI) capture the thalamic damage better than conventional MRI? Are neurological symptoms and disease course better reflected by thalamic alterations as detected by DTI? Are other brain regions also affected? METHODS: Three women with serious neurological deficits due to Stx2-associated HUS were admitted to MRI/DTI at disease onset. Two of them were longitudinally examined. Fractional anisotropy (FA) and mean diffusivity were computed to assess Stx2-caused microstructural damage. RESULTS: Compared to 90 healthy women, all three patients had significantly reduced thalamic FA. Thalamic mean diffusivity was only reduced in two patients. DTI of the longitudinally examined women demonstrated slow normalization of thalamic FA, which was paralleled by clinical improvement. CONCLUSION: Whereas conventional MRI only shows slight alterations based on subjective evaluation, DTI permits quantitative, objective, and longitudinal assessment of cytotoxic cerebral damage in individual patients.

Local estimation of the noise level in MRI using structural adaptation.

We present a method for local estimation of the signal-dependent noise level in magnetic resonance images. The procedure uses a multi-scale approach to adaptively infer on local neighborhoods with similar data distribution. It exploits a maximum-likelihood estimator for the local noise level. The validity of the method was evaluated on repeated diffusion data of a phantom and simulated data using T1-data corrupted with artificial noise. Simulation results were compared with a recently proposed estimate. The method was also applied to a high-resolution diffusion dataset to obtain improved diffusion model estimation results and to demonstrate its usefulness in methods for enhancing diffusion data.

POAS4SPM: a toolbox for SPM to denoise diffusion MRI data.

We present an implementation of a recently developed noise reduction algorithm for dMRI data, called multi-shell position orientation adaptive smoothing (msPOAS), as a toolbox for SPM. The method intrinsically adapts to the structures of different size and shape in dMRI and hence avoids blurring typically observed in non-adaptive smoothing. We give examples for the usage of the toolbox and explain the determination of experiment-dependent parameters for an optimal performance of msPOAS.

High-resolution diffusion kurtosis imaging at 3T enabled by advanced post-processing.

Diffusion Kurtosis Imaging (DKI) is more sensitive to microstructural differences and can be related to more specific micro-scale metrics (e.g., intra-axonal volume fraction) than diffusion tensor imaging (DTI), offering exceptional potential for clinical diagnosis and research into the white and gray matter. Currently DKI is acquired only at low spatial resolution (2-3 mm isotropic), because of the lower signal-to-noise ratio (SNR) and higher artifact level associated with the technically more demanding DKI. Higher spatial resolution of about 1 mm is required for the characterization of fine white matter pathways or cortical microstructure. We used restricted-field-of-view (rFoV) imaging in combination with advanced post-processing methods to enable unprecedented high-quality, high-resolution DKI (1.2 mm isotropic) on a clinical 3T scanner. Post-processing was advanced by developing a novel method for Retrospective Eddy current and Motion ArtifacT Correction in High-resolution, multi-shell diffusion data (REMATCH). Furthermore, we applied a powerful edge preserving denoising method, denoted as multi-shell orientation-position-adaptive smoothing (msPOAS). We demonstrated the feasibility of high-quality, high-resolution DKI and its potential for delineating highly myelinated fiber pathways in the motor cortex. REMATCH performs robustly even at the low SNR level of high-resolution DKI, where standard EC and motion correction failed (i.e., produced incorrectly aligned images) and thus biased the diffusion model fit. We showed that the combination of REMATCH and msPOAS increased the contrast between gray and white matter in mean kurtosis (MK) maps by about 35% and at the same time preserves the original distribution of MK values, whereas standard Gaussian smoothing strongly biases the distribution.

Adaptive smoothing as inference strategy: more specificity for unequally sized or neighbouring regions.

Although spatial smoothing of fMRI data can serve multiple purposes, increasing the sensitivity of activation detection is probably its greatest benefit. However, this increased detection power comes with a loss of specificity when non-adaptive smoothing (i.e. the standard in most software packages) is used. Simulation studies and analysis of experimental data was performed using the R packages neuRosim and fmri. In these studies, we systematically investigated the effect of spatial smoothing on the power and number of false positives in two particular cases that are often encountered in fMRI research: (1) Single condition activation detection for regions that differ in size, and (2) multiple condition activation detection for neighbouring regions. Our results demonstrate that adaptive smoothing is superior in both cases because less false positives are introduced by the spatial smoothing process compared to standard Gaussian smoothing or FDR inference of unsmoothed data.

Dissociations between behavioural and functional magnetic resonance imaging-based evaluations of cognitive function after brain injury.

Functional neuroimaging methods hold promise for the identification of cognitive function and communication capacity in some severely brain-injured patients who may not retain sufficient motor function to demonstrate their abilities. We studied seven severely brain-injured patients and a control group of 14 subjects using a novel hierarchical functional magnetic resonance imaging assessment utilizing mental imagery responses. Whereas the control group showed consistent and accurate (for communication) blood-oxygen-level-dependent responses without exception, the brain-injured subjects showed a wide variation in the correlation of blood-oxygen-level-dependent responses and overt behavioural responses. Specifically, the brain-injured subjects dissociated bedside and functional magnetic resonance imaging-based command following and communication capabilities. These observations reveal significant challenges in developing validated functional magnetic resonance imaging-based methods for clinical use and raise interesting questions about underlying brain function assayed using these methods in brain-injured subjects.

Structural adaptive segmentation for statistical parametric mapping.

Functional Magnetic Resonance Imaging inherently involves noisy measurements and a severe multiple test problem. Smoothing is usually used to reduce the effective number of multiple comparisons and to locally integrate the signal and hence increase the signal-to-noise ratio. Here, we provide a new structural adaptive segmentation algorithm (AS) that naturally combines the signal detection with noise reduction in one procedure. Moreover, the new method is closely related to a recently proposed structural adaptive smoothing algorithm and preserves shape and spatial extent of activation areas without blurring their borders.

High-resolution fMRI: overcoming the signal-to-noise problem.

Increasing the spatial resolution in functional Magnetic Resonance Imaging (fMRI) inherently lowers the signal-to-noise ratio (SNR). In order to still detect functionally significant activations in high-resolution images, spatial smoothing of the data is required. However, conventional non-adaptive smoothing comes with a reduced effective resolution, foiling the benefit of the higher acquisition resolution. We show how our recently proposed structural adaptive smoothing procedure for functional MRI data can improve signal detection of high-resolution fMRI experiments regardless of the lower SNR. The procedure is evaluated on human visual and sensory-motor mapping experiments. In these applications, the higher resolution could be fully utilized and high-resolution experiments were outperforming normal resolution experiments by means of both statistical significance and information content.

Diffusion tensor imaging: structural adaptive smoothing.

Diffusion Tensor Imaging (DTI) data is characterized by a high noise level. Thus, estimation errors of quantities like anisotropy indices or the main diffusion direction used for fiber tracking are relatively large and may significantly confound the accuracy of DTI in clinical or neuroscience applications. Besides pulse sequence optimization, noise reduction by smoothing the data can be pursued as a complementary approach to increase the accuracy of DTI. Here, we suggest an anisotropic structural adaptive smoothing procedure, which is based on the Propagation-Separation method and preserves the structures seen in DTI and their different sizes and shapes. It is applied to artificial phantom data and a brain scan. We show that this method significantly improves the quality of the estimate of the diffusion tensor, by means of both bias and variance reduction, and hence enables one either to reduce the number of scans or to enhance the input for subsequent analysis such as fiber tracking.

Functional MRI of the zebra finch brain during song stimulation suggests a lateralized response topography.

Electrophysiological and activity-dependent gene expression studies of birdsong have contributed to the understanding of the neural representation of natural sounds. However, we have limited knowledge about the overall spatial topography of song representation in the avian brain. Here, we adapt the noninvasive functional MRI method in mildly sedated zebra finches (Taeniopygia guttata) to localize and characterize song driven brain activation. Based on the blood oxygenation level-dependent signal, we observed a differential topographic responsiveness to playback of bird's own song, tutor song, conspecific song, and a pure tone as a nonsong stimulus. The bird's own song caused a stronger response than the tutor song or tone in higher auditory areas. This effect was more pronounced in the medial parts of the forebrain. We found left-right hemispheric asymmetry in sensory responses to songs, with significant discrimination between stimuli observed only in the right hemisphere. This finding suggests that perceptual responses might be lateralized in zebra finches. In addition to establishing the feasibility of functional MRI in sedated songbirds, our results demonstrate spatial coding of song in the zebra finch forebrain, based on developmental familiarity and experience.