Minimal gradient encoding for robust estimation of diffusion anisotropy

Relationship between regional volume changes and water diffusion in fixed marmoset brains: an in vivo and ex vivo comparison

Ex vivo studies of the brain are often employed as experimental systems in neuroscience. In general, brains for ex vivo MRI studies are usually fixed with paraformaldehyde to preserve molecular structure and prevent tissue destruction during long-term storage. As a result, fixing brain tissue causes microstructural changes and a decrease in brain volume. Therefore, the purpose of this study was to investigate the regional effect of brain volume and microstructural changes on the restricted diffusion of water molecules in the common marmoset brain using in vivo and ex vivo brains from the same individual. We used 9.4T magnetic resonance imaging and also compared the T2-weighted images and diffusion-weighted imaging (DWI) data between in vivo and ex vivo brains to investigate changes in brain volume and diffusion of water molecules in 12 common marmosets. We compared fractional anisotropy, mean diffusivity, AD (axial diffusivity), and radial diffusivity values in white matter and gray matter between in vivo and ex vivo brains. We observed that AD showed the strongest correlation with regional volume changes in gray matter. The results showed a strong correlation between AD and changes in brain volume. By comparing the in vivo and ex vivo brains of the same individual, we identified significant correlations between the local effects of perfusion fixation on microstructural and volumetric changes of the brain and alterations in the restricted diffusion of water molecules within the brain. These findings provide valuable insights into the complex relationships between tissue fixation, brain structure, and water diffusion properties in the marmoset brain.

Impact of truncating diffusion MRI scans on diffusional kurtosis imaging

OBJECTIVE

Diffusional kurtosis imaging (DKI) extends diffusion tensor imaging (DTI), characterizing non-Gaussian diffusion effects but requires longer acquisition times. To ensure the robustness of DKI parameters, data acquisition ordering should be optimized allowing for scan interruptions or shortening. Three methodologies were used to examine how reduced diffusion MRI scans impact DKI histogram-metrics: 1) the electrostatic repulsion model (OptEEM); 2) spherical codes (OptSC); 3) random (RandomTRUNC).

MATERIALS AND METHODS

Pre-acquired diffusion multi-shell data from 14 female healthy volunteers (29±5 years) were used to generate reordered data. For each strategy, subsets containing different amounts of the full dataset were generated. The subsampling effects were assessed on histogram-based DKI metrics from tract-based spatial statistics (TBSS) skeletonized maps. To evaluate each subsampling method on simulated data at different SNRs and the influence of subsampling on in vivo data, we used a 3-way and 2-way repeated measures ANOVA, respectively.

RESULTS

Simulations showed that subsampling had different effects depending on DKI parameter, with fractional anisotropy the most stable (up to 5% error) and radial kurtosis the least stable (up to 26% error). RandomTRUNC performed the worst while the others showed comparable results. Furthermore, the impact of subsampling varied across distinct histogram characteristics, the peak value the least affected (OptEEM: up to 5% error; OptSC: up to 7% error) and peak height (OptEEM: up to 8% error; OptSC: up to 11% error) the most affected.

CONCLUSION

The impact of truncation depends on specific histogram-based DKI metrics. The use of a strategy for optimizing the acquisition order is advisable to improve DKI robustness to exam interruptions.

The effect of elimination of gibbs ringing, noise and systematic errors on the DTI metrics and tractography in a rat brain

Diffusion tensor imaging (DTI) metrics and tractography can be biased due to low signal-to-noise ratio (SNR) and systematic errors resulting from image artifacts and imperfections in magnetic field gradients. The imperfections include non-uniformity and nonlinearity, effects caused by eddy currents, and the influence of background and imaging gradients. We investigated the impact of systematic errors on DTI metrics of an isotropic phantom and DTI metrics and tractography of a rat brain measured at high resolution. We tested denoising and Gibbs ringing removal methods combined with the B matrix spatial distribution (BSD) method for magnetic field gradient calibration. The results showed that the performance of the BSD method depends on whether Gibbs ringing is removed and the effectiveness of stochastic error removal. Region of interest (ROI)-based analysis of the DTI metrics showed that, depending on the size of the ROI and its location in space, correction methods can remove systematic bias to varying degrees. The preprocessing pipeline proposed and dedicated to this type of data together with the BSD method resulted in an even - 90% decrease in fractional anisotropy (FA) (globally and locally) in the isotropic phantom and - 45% in the rat brain. The largest global changes in the rat brain tractogram compared to the standard method without preprocessing (sDTI) were noticed after denoising. The direction of the first eigenvector obtained from DTI after denoising, Gibbs ringing removal and BSD differed by an average of 56 and 10 degrees in the ROI from sDTI and from sDTI after denoising and Gibbs ringing removal, respectively. The latter can be identified with the amount of improvement in tractography due to the elimination of systematic errors related to imperfect magnetic field gradients. Based on the results, the systematic bias for high resolution data mainly depended on SNR, but the influence of non-uniform gradients could also be seen. After denoising, the BSD method was able to further correct both the metrics and tractography of the diffusion tensor in the rat brain by taking into account the actual distribution of magnetic field gradients independent of the examined object and uniquely dependent on the scanner and sequence. This means that in vivo studies are also subject to this type of errors, which should be taken into account when processing such data.

Effect of number of diffusion encoding directions in neonatal diffusion tensor imaging using Tract-Based Spatial Statistical analysis

Diffusion tensor imaging (DTI) has been used to study the developing brain in early childhood, infants and in utero studies. In infants, number of used diffusion encoding directions has traditionally been smaller in earlier studies down to the minimum of 6 orthogonal directions. Whereas the more recent studies often involve more directions, number of used directions remain an issue when acquisition time is optimized without compromising on data quality and in retrospective studies. Variability in the number of used directions may introduce bias and uncertainties to the DTI scalar estimates that affect cross-sectional and longitudinal study of the brain. We analysed DTI images of 133 neonates, each data having 54 directions after quality control, to evaluate the effect of number of diffusion weighting directions from 6 to 54 with interval of 6 to the DTI scalars with Tract-Based Spatial Statistics (TBSS) analysis. The TBSS analysis was applied to DTI scalar maps, and the mean region of interest (ROI) values were extracted using JHU atlas. We found significant bias in ROI mean values when only 6 directions were used (positive in fractional anisotropy [FA] and negative in fractional anisotropy [MD], axial diffusivity [AD] and fractional anisotropy [RD]), while when using 24 directions and above, the difference to scalar values calculated from 54 direction DTI was negligible. In repeated measures voxel-wise analysis, notable differences to 54 direction DTI were observed with 6, 12 and 18 directions. DTI measurements from data with at least 24 directions may be used in comparisons with DTI measurements from data with higher numbers of directions.

Cognitive-motor telerehabilitation in multiple sclerosis (CoMoTeMS): study protocol for a randomised controlled trial

Background

The management of cognitive impairment is an important goal in the treatment of multiple sclerosis (MS). While cognitive rehabilitation has been proven to be effective in improving cognitive performance in MS, research in the elderly indicates a higher effectiveness of combined cognitive-motor rehabilitation. Here, we present the protocol of a randomised controlled clinical trial to assess whether a combined cognitive-motor telerehabilitation programme is more effective in improving working memory than only cognitive or motor training.

Methods/design

The CoMoTeMS-trial is a two-centre, randomised, controlled and blinded clinical trial. A total of 90 patients with MS will receive 12 weeks of either a combined cognitive-motor telerehabilitation programme or only cognitive or motor training. The primary outcome is a change in the digit span backwards. Secondary outcomes are other cognitive changes (Brief International Cognitive Assessment for Multiple Sclerosis and Backward Corsi), Expanded Disability Status Scale (EDSS), 6-Min Walk Test, 25-Foot Walk Test, 9-Hole Peg Test, anxiety and depression, fatigue, quality of life, cognitive and physical activity level, electroencephalography and magnetic resonance imaging of the brain.

Discussion

We hypothesise that the improvement in digit span backwards after 12 weeks of treatment will be significantly higher in the group treated with the combined cognitive-motor telerehabilitation programme, compared to the groups receiving only cognitive and only motor training.

Trial registration

ClinicalTrials.gov NCT05355389. Registered on 2 May 2022.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13063-022-06697-9.

Effect of number of diffusion-encoding directions in diffusion metrics of 5-year-olds using tract-based spatial statistical analysis

Methodological aspects and effects of different imaging parameters on DTI (diffusion tensor imaging) results and their reproducibility have been recently studied comprehensively in adult populations. Although MR imaging of children's brains has become common, less interest has been focussed on researching whether adult-based optimised parameters and pre-processing protocols can be reliably applied to paediatric populations. Furthermore, DTI scalar values of preschool aged children are rarely reported. We gathered a DTI dataset from 5-year-old children (N = 49) to study the effect of the number of diffusion-encoding directions on the reliability of resultant scalar values with TBSS (tract-based spatial statistics) method. Additionally, the potential effect of within-scan head motion on DTI scalars was evaluated. Reducing the number of diffusion-encoding directions deteriorated both the accuracy and the precision of all DTI scalar values. To obtain reliable scalar values, a minimum of 18 directions for TBSS was required. For TBSS fractional anisotropy values, the intraclass correlation coefficient with two-way random-effects model (ICC[2,1]) for the subsets of 6 to 66 directions ranged between 0.136 [95%CI 0.0767;0.227] and 0.639 [0.542;0.740], whereas the corresponding values for subsets of 18 to 66 directions were 0.868 [0.815;0.913] and 0.995 [0.993;0.997]. Following the exclusion of motion-corrupted volumes, minor residual motion did not associate with the scalar values. A minimum of 18 diffusion directions is recommended to result in reliable DTI scalar results with TBSS. We suggest gathering extra directions in paediatric DTI to enable exclusion of volumes with motion artefacts and simultaneously preserve the overall data quality.

A groupwise registration and tractography framework for cardiac myofiber architecture description by diffusion MRI: An application to the ventricular junctions

Background

Knowledge of the normal myocardial–myocyte orientation could theoretically allow the definition of relevant quantitative biomarkers in clinical routine to diagnose heart pathologies. A whole heart diffusion tensor template representative of the global myofiber organization over species is therefore crucial for comparisons across populations. In this study, we developed a groupwise registration and tractography framework to resolve the global myofiber arrangement of large mammalian sheep hearts. To demonstrate the potential application of the proposed method, a novel description of sub-regions in the intraventricular septum is presented.

Methods

Three explanted sheep (ovine) hearts (size ~12×8×6 cm3, heart weight ~ 150 g) were perfused with contrast agent and fixative and imaged in a 9.4T magnet. A group-wise registration of high-resolution anatomical and diffusion-weighted images were performed to generate anatomical and diffusion tensor templates. Diffusion tensor metrics (eigenvalues, eigenvectors, fractional anisotropy …) were computed to provide a quantitative and spatially-resolved analysis of cardiac microstructure. Then tractography was performed using deterministic and probabilistic algorithms and used for different purposes: i) Visualization of myofiber architecture, ii) Segmentation of sub-area depicting the same fiber organization, iii) Seeding and Tract Editing. Finally, dissection was performed to confirm the existence of macroscopic structures identified in the diffusion tensor template.

Results

The template creation takes advantage of high-resolution anatomical and diffusion-weighted images obtained at an isotropic resolution of 150 μm and 600 μm respectively, covering ventricles and atria and providing information on the normal myocardial architecture. The diffusion metric distributions from the template were found close to the one of the individual samples validating the registration procedure. Small new sub-regions exhibiting spatially sharp variations in fiber orientation close to the junctions of the septum and ventricles were identified. Each substructure was defined and represented using streamlines. The existence of a fiber-bundles in the posterior junction was validated by anatomical dissection. A complex structural organization of the anterior junction in comparison to the posterior junction was evidenced by the high-resolution acquisition.

Conclusions

A new framework combining cardiac template generation and tractography was applied on the whole sheep heart. The framework can be used for anatomical investigation, characterization of microstructure and visualization of myofiber orientation across samples. Finally, a novel description of the ventricular junction in large mammalian sheep hearts was proposed.

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 T₂-DKI-FWE model that exploits the T₂ relaxation properties of both compartments, thereby better conditioning the parameter estimation problem and providing, at the same time, an additional potential biomarker (the T₂ of tissue). In our approach, the T₂ of tissue is estimated as an unknown parameter, whereas the T₂ 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 T₂ of free water is studied. Next, the improved conditioning of T₂-DKI-FWE compared to DKI-FWE is illustrated using the Cramér-Rao lower bound matrix. Finally, the performance of the T₂-DKI-FWE model is compared to that of the DKI-FWE and T₂-DKI models on both simulated and real datasets. The error due to a biased approximation of the T₂ 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 T₂-DKI-FWE model is beneficial for the identifiability of the model parameters. Our results suggest that the T₂-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 T₂ 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.

Renal Damages in Deoxycorticosterone Acetate-Salt Hypertensive Rats: Assessment with Diffusion Tensor Imaging and T2-mapping

PURPOSE

This study aimed to investigate the feasibility of diffusion tensor imaging (DTI) and T2-mapping to assess temporal renal damage in deoxycorticosterone acetate-salt (DOCA-salt) hypertensive rats and compare the results with histopathologic and immunohistochemical findings.

PROCEDURES

After baseline renal magnetic resonance imaging (MRI), 24 out of 30 uninephrectomized Sprague-Dawley rats with DOCA-salt-induced hypertension were divided equally into four groups. Group 1 had renal MRI at weeks 2, 4, 6, and 8, and groups 2, 3, and 4 had MRI at weeks 2, 4, and 6, respectively. The remaining 6 rats were used as sham controls. The renal cortex and outer and inner stripes of the outer medulla were examined over time using fractional anisotropy (FA), apparent diffusion coefficient (ADC), and T2-mapping, and the results were compared with baseline values. The degree of glomerular and tubular injury, endothelial cell thickening, hyaline arteriolosclerosis, macrophage infiltration, microcyst formation, and fibrosis in different zones at different time points in the DOCA-salt rats were compared with controls.

RESULTS

Compared with baseline values, DOCA-salt rats demonstrated a significant decrease in renal cortical FA from week 4 to week 8 (0.244 ± 0.015 vs 0.172 ± 0.014-0.150 ± 0.016, P = 0.018-0.002), corresponding to significantly more glomerular damage, arteriolosclerosis, macrophage infiltration, and fibrosis. The DOCA-salt rats had significantly increased cortical ADC and T2 values at weeks 6 and 8 (1.778 ± 0.051 × 10^-3 mm²/s vs 1.872 ± 0.058-1.917 ± 0.066 × 10^-3 mm²/s; 93.7 ± 4.9 ms vs 98.0 ± 2.9-100.7 ± 4.0 ms, respectively, all P < 0.05), consistent with excessively fluid-filled microcysts (aquaporin-2+). Despite DOCA-salt rats harbored markedly increased fibrosis in outer and inner stripes of the outer medulla at weeks 6 and 8, only nonsignificant decreases in FA were observed in comparison with the controls suggesting that only limited microstructural changes were present.

CONCLUSIONS

Renal cortical FA is useful for the early detection and monitoring of renal damage in DOCA-salt hypertensive rats.

Accelerating in vivo fast spin echo high angular resolution diffusion imaging with an isotropic resolution in mice through compressed sensing

PURPOSE

Echo planar imaging (EPI) is commonly used to acquire the many volumes needed for high angular resolution diffusion Imaging (HARDI), posing a higher risk for artifacts, such as distortion and deformation. An alternative to EPI is fast spin echo (FSE) imaging, which has fewer artifacts but is inherently slower. The aim is to accelerate FSE such that a HARDI data set can be acquired in a time comparable to EPI using compressed sensing.

METHODS

Compressed sensing was applied in either q-space or simultaneously in k-space and q-space, by undersampling the k-space in the phase-encoding direction or retrospectively eliminating diffusion directions for different degrees of undersampling. To test the replicability of the acquisition and reconstruction, brain data were acquired from six mice, and a numerical phantom experiment was performed. All HARDI data were analyzed individually using constrained spherical deconvolution, and the apparent fiber density and complexity metric were evaluated, together with whole-brain tractography.

RESULTS

The apparent fiber density and complexity metric showed relatively minor differences when only q-space undersampling was used, but deteriorate when k-space undersampling was applied. Likewise, the tract density weighted image showed good results when only q-space undersampling was applied using 15 directions or more, but information was lost when fewer volumes or k-space undersampling were used.

CONCLUSION

It was found that acquiring 15 to 20 diffusion directions with a full k-space and reconstructed using compressed sensing could suffice for a replicable measurement of quantitative measures in mice, where areas near the sinuses and ear cavities are untainted by signal loss.

A data-driven approach to optimising the encoding for multi-shell diffusion MRI with application to neonatal imaging

Diffusion MRI has the potential to provide important information about the connectivity and microstructure of the human brain during normal and abnormal development, noninvasively and in vivo. Recent developments in MRI hardware and reconstruction methods now permit the acquisition of large amounts of data within relatively short scan times. This makes it possible to acquire more informative multi-shell data, with diffusion sensitisation applied along many directions over multiple b-value shells. Such schemes are characterised by the number of shells acquired, and the specific b-value and number of directions sampled for each shell. However, there is currently no clear consensus as to how to optimise these parameters. In this work, we propose a means of optimising multi-shell acquisition schemes by estimating the information content of the diffusion MRI signal, and optimising the acquisition parameters for sensitivity to the observed effects, in a manner agnostic to any particular diffusion analysis method that might subsequently be applied to the data. This method was used to design the acquisition scheme for the neonatal diffusion MRI sequence used in the developing Human Connectome Project (dHCP), which aims to acquire high quality data and make it freely available to the research community. The final protocol selected by the algorithm, and currently in use within the dHCP, consists of 20 b=0 images and diffusion-weighted images at b = 400, 1000 and 2600 s/mm² with 64, 88 and 128 directions per shell, respectively.

Brain MRI Diffusion Encoding Direction Number Affects Tract-Based Spatial Statistics Results in Multiple Sclerosis

BACKGROUND AND PURPOSE

Diffusion tensor imaging (DTI) is a promising approach to detect the underlying brain pathology. These alterations can be seen in several diseases such as multiple sclerosis. Tract-based spatial statistics (TBSS) is an easy to use and robust way for analyzing diffusion data. The effect of acquisition parameters of DTI on TBSS has not been evaluated, especially the number of diffusion encoding directions (NDED), which is directly proportional with scan time.

METHODS

We analyzed a large set of DTI data of healthy controls (N = 126) and multiple sclerosis patients (N = 78). The highest NDED (60 directions) was reduced and a tensor calculation was done separately for every subset. We calculated the mean and standard deviation of DTI parameters under the white matter mask. Moreover, the FMRIB Software Library TBSS pipeline was used on DTI images with 15, 30, 45, and 60 directions to compare differences between groups. Mean DTI parameters were compared between groups as a function of NDED.

RESULTS

The mean value of FA and AD decreased with increasing number of directions. This was more pronounced in areas with smaller FA values. RD and MD were constant. The skeleton size reduced with elevating NDED along with the number of significant voxels. The TBSS analysis showed significant differences between groups throughout the majority of the skeleton and the group difference was associated with NDED.

CONCLUSION

Our results suggested that results of TBSS depended on the NDED, which should be considered when comparing DTI data with varying protocols.

Effects of signal averaging, gradient encoding scheme, and spatial resolution on diffusion kurtosis imaging: An empirical study using 7T MRI

BACKGROUND

Although diffusion gradient directions and b-values have been optimized for diffusion kurtosis imaging (DKI), little is known about the effect of signal averaging on DKI reliability.

PURPOSE

To evaluate how signal averaging influences the reliability of DKI indices using two gradient encoding schemes with three spatial resolutions.

STUDY TYPE

Prospective.

ANIMAL MODEL

Fifteen naïve Sprague-Dawley rats.

FIELD STRENGTH/SEQUENCE

DKI was performed at 7T using two schemes (30 directions with three b-values [30d-3b] and six directions with 15 b-values [6d-15b]), three resolutions, and eight repetitions.

ASSESSMENT

DKI reliability was assessed using voxelwise relative error (σ) and test-retest error of fractional anisotropy (FA), mean diffusivity (MD), and mean kurtosis (MK) within gray matter (GM) and white matter (WM). The number of excitations (NEX) was optimized by considering DKI reliability. The influence of the partial volume effect (PVE) was also assessed.

STATISTICAL TEST

One-way analysis of variance.

RESULTS

The 30d-3b scheme, compared with the 6d-15b scheme, exhibited apparently smaller σ_FA and σ_MK (eg, at NEX 1, in GM, for three resolutions, σ_FA : 19.9-38.2% vs. 34.2-61.4%, σ_MK : 6.9-11.4% vs. 14.1-15.4%) and similar σ_MD (all differences between two schemes <1.6%). The optimal NEX was determined as 2 for enabling a reliable measurement of DKI-derived indices. The PVE at the lowest resolution apparently increased σ_FA for both schemes (19.9% for 30d-3b and 34.2% for 6d-15b) and σ_MK for the 6d-15b scheme (14.7%) in GM, and exerted lower effects on MK values for the 30d-3b scheme (P > 0.05).

DATA CONCLUSION

A higher number of diffusion directions would benefit FA and MK estimation. A higher spatial resolution helps to reduce PVE. By using the 30d-3b scheme, MK is considered a robust index to reflect microstructural changes in GM and WM. We propose a systematic approach to determine the optimal DKI protocols for appropriate preclinical settings.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;50:1593-1603.

Diffusion tensor imaging of the sciatic nerve in Charcot-Marie-Tooth disease type I patients: a prospective case-control study

OBJECTIVES

This study aimed to evaluate whether diffusion tensor imaging (DTI) parameters and cross-sectional area (CSA) can differentiate between the sciatic nerve of Charcot-Marie-Tooth (CMT) disease type I (demyelinating form) patients and that of controls.

METHODS

This prospective comparison study included 18 CMT type I patients and 18 age/sex-matched volunteers. Magnetic resonance imaging including DTI and axial T2-weighted Dixon sequence was performed for each subject. Region of interest analysis was independently performed by two radiologists on each side of the sciatic nerve at four levels: hamstring tendon origin (level 1), lesser trochanter of the femur (level 2), gluteus maximus tendon insertion (level 3), and mid-femur (level 4). Fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), and radial diffusivity (RD) were calculated. The CSA of the sciatic nerve bundle was measured using axial water-only image at each level. Comparisons of DTI parameters between the two groups were performed using the two-sample t test and Mann-Whitney U test. Interobserver agreement analysis was also conducted.

RESULTS

Interobserver agreement was excellent for all DTI parameter analyses. FA was significantly lower at all four levels in CMT patients than controls. RD, MD, and CSA were significantly higher at all four levels in CMT patients. AD was significantly higher at level 2 in CMT patients.

CONCLUSION

DTI assessment of the sciatic nerve is reproducible and can discriminate the demyelinating nerve pathology of CMT type I patients from normal nerves. The CSA of the sciatic nerve is also a potential parameter for diagnosing nerve abnormality in CMT type I patients.

KEY POINTS

• Diffusion tensor imaging parameters of the sciatic nerve at proximal to mid-femur level revealed significant differences between the Charcot-Marie-Tooth disease patients and controls. • The cross-sectional area of the sciatic nerve was significantly larger in the Charcot-Marie-Tooth disease patients. • Interobserver agreement was excellent (intraclass coefficient > 0.8) for all diffusion tensor imaging parameter analyses.

Minimal number of gradient directions for robust measurement of spherical mean diffusion weighted signal

PURPOSE

Determination of the minimum number of gradient directions (N_min) for robust measurement of spherical mean diffusion weighted signal (S¯).

METHODS

Computer simulations were employed to characterize the relative standard deviation (RSD) of the measured spherical mean signal as a function of the number of gradient directions (N). The effects of diffusion weighting b-value and signal-to-noise ratio (SNR) were investigated. Multi-shell high angular resolution Human Connectome Project diffusion data were analyzed to support the simulation results.

RESULTS

RSD decreases with increasing N, and the minimum number of N needed for RSD ≤ 5% is referred to as N_min. At high SNRs, N_min increases with increasing b-value to achieve sufficient sampling. Simulations showed that N_min is linearly dependent on the b-value. At low SNRs, N_min increases with increasing b-value to reduce the noise. RSD can be estimated as σS¯N, where σ = 1/SNR is the noise level. The experimental results were in good agreement with the simulation results. The spherical mean signal can be measured accurately with a subset of gradient directions.

CONCLUSION

As N_min is affected by b-value and SNR, we recommend using 10 × b / b₁ (b₁ = 1 ms/μm²) uniformly distributed gradient directions for typical human diffusion studies with SNR ~ 20 for robust spherical mean signal measurement.

Diffusion kurtosis imaging with free water elimination: A bayesian estimation approach

PURPOSE

Diffusion kurtosis imaging (DKI) is an advanced magnetic resonance imaging modality that is known to be sensitive to changes in the underlying microstructure of the brain. Image voxels in diffusion weighted images, however, are typically relatively large making them susceptible to partial volume effects, especially when part of the voxel contains cerebrospinal fluid. In this work, we introduce the "Diffusion Kurtosis Imaging with Free Water Elimination" (DKI-FWE) model that separates the signal contributions of free water and tissue, where the latter is modeled using DKI.

THEORY AND METHODS

A theoretical study of the DKI-FWE model, including an optimal experiment design and an evaluation of the relative goodness of fit, is carried out. To stabilize the ill-conditioned estimation process, a Bayesian approach with a shrinkage prior (BSP) is proposed. In subsequent steps, the DKI-FWE model and the BSP estimation approach are evaluated in terms of estimation error, both in simulation and real data experiments.

RESULTS

Although it is shown that the DKI-FWE model parameter estimation problem is ill-conditioned, DKI-FWE was found to describe the data significantly better compared to the standard DKI model for a large range of free water fractions. The acquisition protocol was optimized in terms of the maximally attainable precision of the DKI-FWE model parameters. The BSP estimator is shown to provide reliable DKI-FWE model parameter estimates.

CONCLUSION

The combination of the DKI-FWE model with BSP is shown to be a feasible approach to estimate DKI parameters, while simultaneously eliminating free water partial volume effects. Magn Reson Med 80:802-813, 2018. © 2018 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution NonCommercial License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited and is not used for commercial purposes.

Single- and Multiple-Shell Uniform Sampling Schemes for Diffusion MRI Using Spherical Codes

In diffusion MRI (dMRI), a good sampling scheme is important for efficient acquisition and robust reconstruction. Diffusion weighted signal is normally acquired on single or multiple shells in q-space. Signal samples are typically distributed uniformly on different shells to make them invariant to the orientation of structures within tissue, or the laboratory coordinate frame. The Electrostatic Energy Minimization (EEM) method, originally proposed for single shell sampling scheme in dMRI, was recently generalized to multi-shell schemes, called Generalized EEM (GEEM). GEEM has been successfully used in the Human Connectome Project (HCP). However, EEM does not directly address the goal of optimal sampling, i.e., achieving large angular separation between sampling points. In this paper, we propose a more natural formulation, called Spherical Code (SC), to directly maximize the minimal angle between different samples in single or multiple shells. We consider not only continuous problems to design single or multiple shell sampling schemes, but also discrete problems to uniformly extract sub-sampled schemes from an existing single or multiple shell scheme, and to order samples in an existing scheme. We propose five algorithms to solve the above problems, including an incremental SC (ISC), a sophisticated greedy algorithm called Iterative Maximum Overlap Construction (IMOC), an 1-Opt greedy method, a Mixed Integer Linear Programming (MILP) method, and a Constrained Non-Linear Optimization (CNLO) method. To our knowledge, this is the first work to use the SC formulation for single or multiple shell sampling schemes in dMRI. Experimental results indicate that SC methods obtain larger angular separation and better rotational invariance than the state-of-the-art EEM and GEEM. The related codes and a tutorial have been released in DMRITool.

Effects of Orientation and Anisometry of Magnetic Resonance Imaging Acquisitions on Diffusion Tensor Imaging and Structural Connectomes

Diffusion-weighted imaging (DWI) quantifies water molecule diffusion within tissues and is becoming an increasingly used technique. However, it is very challenging as correct quantification depends on many different factors, ranging from acquisition parameters to a long pipeline of image processing. In this work, we investigated the influence of voxel geometry on diffusion analysis, comparing different acquisition orientations as well as isometric and anisometric voxels. Diffusion-weighted images of one rat brain were acquired with four different voxel geometries (one isometric and three anisometric in different directions) and three different encoding orientations (coronal, axial and sagittal). Diffusion tensor scalar measurements, tractography and the brain structural connectome were analyzed for each of the 12 acquisitions. The acquisition direction with respect to the main magnetic field orientation affected the diffusion results. When the acquisition slice-encoding direction was not aligned with the main magnetic field, there were more artifacts and a lower signal-to-noise ratio that led to less anisotropic tensors (lower fractional anisotropic values), producing poorer quality results. The use of anisometric voxels generated statistically significant differences in the values of diffusion metrics in specific regions. It also elicited differences in tract reconstruction and in different graph metric values describing the brain networks. Our results highlight the importance of taking into account the geometric aspects of acquisitions, especially when comparing diffusion data acquired using different geometries.

The effect of diffusion gradient direction number on corticospinal tractography in the human brain: an along-tract analysis

OBJECTIVES

We evaluated diffusion imaging measures of the corticospinal tract obtained with a probabilistic tractography algorithm applied to data of two acquisition protocols based on different numbers of diffusion gradient directions (NDGDs).

MATERIALS AND METHODS

The corticospinal tracts (CST) of 18 healthy subjects were delineated using 22 and 66-NDGD data. An along-tract analysis of diffusion metrics was performed to detect possible local differences due to NDGD.

RESULTS

FA values at 22-NDGD showed an increase along the central portion of the CST. The mean of partial volume fraction of the orientation of the second fiber (f2) was higher at 66-NDGD bilaterally, because for 66-NDGD data the algorithm more readily detects dominant fiber directions beyond the first, thus the increase in FA at 22-NDGD is due to a substantially reduced detection of crossing fiber volume. However, the good spatial correlation between the tracts drawn at 22 and 66 NDGD shows that the extent of the tract can be successfully defined even at lower NDGD.

CONCLUSIONS

Given the spatial tract localization obtained even at 22-NDGD, local analysis of CST can be performed using a NDGD compatible with clinical protocols. The probabilistic approach was particularly powerful in evaluating crossing fibers when present.

A Novel Technique for Region and Linguistic Specific nTMS-based DTI Fiber Tracking of Language Pathways in Brain Tumor Patients

Navigated transcranial magnetic stimulation (nTMS) has recently been introduced as a non-invasive tool for functional mapping of cortical language areas prior to surgery. It correlates well with intraoperative neurophysiological monitoring (IONM) findings, allowing defining the best surgical strategy to preserve cortical language areas during surgery for language-eloquent tumors. Nevertheless, nTMS allows only for cortical mapping and postoperative language deficits are often caused by injury to subcortical language pathways. Nowadays, the only way to preoperatively visualize language subcortical white matter tracts consists in DTI fiber tracking (DTI-FT). However, standard DTI-FT is based on anatomical landmarks that vary interindividually and can be obscured by the presence of the tumor itself. It has been demonstrated that combining nTMS with DTI-FT allows for a more reliable visualization of the motor pathway in brain tumor patients. Nevertheless, no description about such a combination has been reported for the language network. The aim of the present study is to describe and assess the feasibility and reliability of using cortical seeding areas defined by error type-specific nTMS language mapping (nTMS-positive spots) to perform DTI-FT in patients affected by language-eloquent brain tumors. We describe a novel technique for a nTMS-based DTI-FT to visualize the complex cortico-subcortical connections of the language network. We analyzed quantitative findings, such as fractional anisotropy values and ratios, and the number of visualized connections of nTMS-positive spots with subcortical pathways, and we compared them with results obtained by using the standard DTI-FT technique. We also analyzed the functional concordance between connected cortical nTMS-positive spots and subcortical pathways, and the likelihood of connection for nTMS-positive vs. nTMS-negative cortical spots. We demonstrated, that the nTMS-based approach, especially what we call the “single-spot” strategy, is able to provide a reliable and more detailed reconstruction of the complex cortico-subcortical language network as compared to the standard DTI-FT. We believe this technique represents a beneficial new strategy for customized preoperative planning in patients affected by tumors in presumed language eloquent location, providing anatomo-functional information to plan language-preserving surgery.

Reduced Field-of-View Diffusion-Weighted Imaging of the Lumbosacral Enlargement: A Pilot In Vivo Study of the Healthy Spinal Cord at 3T

Diffusion tensor imaging (DTI) has recently started to be adopted into clinical investigations of spinal cord (SC) diseases. However, DTI applications to the lower SC are limited due to a number of technical challenges, related mainly to the even smaller size of the SC structure at this level, its position relative to the receiver coil elements and the effects of motion during data acquisition. Developing methods to overcome these problems would offer new means to gain further insights into microstructural changes of neurological conditions involving the lower SC, and in turn could help explain symptoms such as bladder and sexual dysfunction. In this work, the feasibility of obtaining grey and white matter (GM/WM) DTI indices such as axial/radial/mean diffusivity (AD/RD/MD) and fractional anisotropy (FA) within the lumbosacral enlargement (LSE) was investigated using a reduced field-of-view (rFOV) single-shot echo-planar imaging (ss-EPI) acquisition in 14 healthy participants using a clinical 3T MR system. The scan-rescan reproducibility of the measurements was assessed by calculating the percentage coefficient of variation (%COV). Mean FA was higher in WM compared to GM (0.58 and 0.4 in WM and GM respectively), AD and MD were higher in WM compared to GM (1.66 μm2ms-1 and 0.94 μm2ms-1 in WM and 1.2 μm2ms-1 and 0.82 μm2ms-1 in GM for AD and MD respectively) and RD was lower in WM compared to GM (0.58 μm2ms-1 and 0.63 μm2ms-1 respectively). The scan-rescan %COV was lower than 10% in all cases with the highest values observed for FA and the lowest for MD. This pilot study demonstrates that it is possible to obtain reliable tissue-specific estimation of DTI indices within the LSE using a rFOV ss-EPI acquisition. The DTI acquisition and analysis protocol presented here is clinically feasible and may be used in future investigations of neurological conditions implicating the lower SC.

REDUCING CSF PARTIAL VOLUME EFFECTS TO ENHANCE DIFFUSION TENSOR IMAGING METRICS OF BRAIN MICROSTRUCTURE

Technological advances over recent decades now allow for in vivo observation of human brain tissue through the use of neuroimaging methods. While this field originated with techniques capable of capturing macrostructural details of brain anatomy, modern methods such as diffusion tensor imaging (DTI) that are now regularly implemented in research protocols have the ability to characterize brain microstructure. DTI has been used to reveal subtle micro-anatomical abnormalities in the prodromal phase ofº various diseases and also to delineate "normal" age-related changes in brain tissue across the lifespan. Nevertheless, imaging artifact in DTI remains a significant limitation for identifying true neural signatures of disease and brain-behavior relationships. Cerebrospinal fluid (CSF) contamination of brain voxels is a main source of error on DTI scans that causes partial volume effects and reduces the accuracy of tissue characterization. Several methods have been proposed to correct for CSF artifact though many of these methods introduce new limitations that may preclude certain applications. The purpose of this review is to discuss the complexity of signal acquisition as it relates to CSF artifact on DTI scans and review methods of CSF suppression in DTI. We will then discuss a technique that has been recently shown to effectively suppress the CSF signal in DTI data, resulting in fewer errors and improved measurement of brain tissue. This approach and related techniques have the potential to significantly improve our understanding of "normal" brain aging and neuropsychiatric and neurodegenerative diseases. Considerations for next-level applications are discussed.

K-Optimal Gradient Encoding Scheme for Fourth-Order Tensor-Based Diffusion Profile Imaging

The design of an optimal gradient encoding scheme (GES) is a fundamental problem in diffusion MRI. It is well studied for the case of second-order tensor imaging (Gaussian diffusion). However, it has not been investigated for the wide range of non-Gaussian diffusion models. The optimal GES is the one that minimizes the variance of the estimated parameters. Such a GES can be realized by minimizing the condition number of the design matrix (K-optimal design). In this paper, we propose a new approach to solve the K-optimal GES design problem for fourth-order tensor-based diffusion profile imaging. The problem is a nonconvex experiment design problem. Using convex relaxation, we reformulate it as a tractable semidefinite programming problem. Solving this problem leads to several theoretical properties of K-optimal design: (i) the odd moments of the K-optimal design must be zero; (ii) the even moments of the K-optimal design are proportional to the total number of measurements; (iii) the K-optimal design is not unique, in general; and (iv) the proposed method can be used to compute the K-optimal design for an arbitrary number of measurements. Our Monte Carlo simulations support the theoretical results and show that, in comparison with existing designs, the K-optimal design leads to the minimum signal deviation.

A model-based reconstruction for undersampled radial spin-echo DTI with variational penalties on the diffusion tensor

Radial spin-echo diffusion imaging allows motion-robust imaging of tissues with very low T2 values like articular cartilage with high spatial resolution and signal-to-noise ratio (SNR). However, in vivo measurements are challenging, due to the significantly slower data acquisition speed of spin-echo sequences and the less efficient k-space coverage of radial sampling, which raises the demand for accelerated protocols by means of undersampling. This work introduces a new reconstruction approach for undersampled diffusion-tensor imaging (DTI). A model-based reconstruction implicitly exploits redundancies in the diffusion-weighted images by reducing the number of unknowns in the optimization problem and compressed sensing is performed directly in the target quantitative domain by imposing a total variation (TV) constraint on the elements of the diffusion tensor. Experiments were performed for an anisotropic phantom and the knee and brain of healthy volunteers (three and two volunteers, respectively). Evaluation of the new approach was conducted by comparing the results with reconstructions performed with gridding, combined parallel imaging and compressed sensing and a recently proposed model-based approach. The experiments demonstrated improvements in terms of reduction of noise and streaking artifacts in the quantitative parameter maps, as well as a reduction of angular dispersion of the primary eigenvector when using the proposed method, without introducing systematic errors into the maps. This may enable an essential reduction of the acquisition time in radial spin-echo diffusion-tensor imaging without degrading parameter quantification and/or SNR.

Effect of increasing diffusion gradient direction number on diffusion tensor imaging fiber tracking in the human brain

Objective

To assess the effects of varying the number of diffusion gradient directions (NDGDs) on diffusion tensor fiber tracking (FT) in human brain white matter using tract characteristics.

Materials and Methods

Twelve normal volunteers underwent diffusion tensor imaging (DTI) scanning with NDGDs of 6, 11, 15, 21, and 31 orientations. Three fiber tract groups, including the splenium of the corpus callosum (CC), the entire CC, and the full brain tract, were reconstructed by deterministic DTI-FT. Tract architecture was first qualitatively evaluated by visual observation. Six quantitative tract characteristics, including the number of fibers (NF), average length (AL), fractional anisotropy (FA), relative anisotropy (RA), mean diffusivity (MD), and volume ratio (VR) were measured for the splenium of the CC at the tract branch level, for the entire CC at tract level, and for the full brain tract at the whole brain level. Visual results and those of NF, AL, FA, RA, MD, and VR were compared among the five different NDGDs.

Results

The DTI-FT with NDGD of 11, 15, 21, and 31 orientations gave better tracking results compared with NDGD of 6 after the visual evaluation. NF, FA, RA, MD, and VR values with NDGD of six were significantly greater (smallest p = 0.001 to largest p = 0.042) than those with four other NDGDs (11, 15, 21, or 31 orientations), whereas AL measured with NDGD of six was significantly smaller (smallest p = 0.001 to largest p = 0.041) than with four other NDGDs (11, 15, 21, or 31 orientations). No significant differences were observed in the results among the four NDGD groups of 11, 15, 21, and 31 directions (smallest p = 0.059 to largest p = 1.000).

Conclusion

The main fiber tracts were detected with NDGD of six orientations; however, the use of larger NDGD (≥ 11 orientations) could provide improved tract characteristics at the expense of longer scanning time.

Effects of rejecting diffusion directions on tensor-derived parameters

Diffusion Tensor Imaging (DTI) is adversely affected by subject motion. It is necessary to discard the corrupted images before diffusion parameter estimation. However, the consequences of rejecting those images are not well understood. In this study, we investigated the effects of excluding one or more volumes of diffusion weighted images by analyzing the changes in fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (AD), radial diffusivity (RD) and the primary eigenvector (V1). Based on the full set of diffusion images acquired by the Jones30 diffusion scheme, we generated incomplete sets of at least six in three different ways: random, uniform and clustered rejections. The results showed that MD was not significantly affected by rejecting diffusion directions. In the cases of random rejections, FA, AD, RD and V1 were overestimated more greatly with increasing number of rejections and the overestimations were worse in low FA regions than high FA regions. For uniform rejections, at which the remaining diffusion directions are evenly distributed on a sphere, little change was observed in FA and in V1. Clustered rejections, on the other hand, displayed the most significant overestimation of the parameters, and the resulting accuracy depended on the relative orientation of the underlying fibers with respect to the excluded directions. In practice, if diffusion direction data is excluded, it is important to note the number and location of directions rejected, in order to make a more precise analysis of the data.

Linked alterations in gray and white matter morphology in adults with high-functioning autism spectrum disorder: a multimodal brain imaging study

Growing evidence suggests that a broad range of behavioral anomalies in people with autism spectrum disorder (ASD) can be linked with morphological and functional alterations in the brain. However, the neuroanatomical underpinnings of ASD have been investigated using either structural magnetic resonance imaging (MRI) or diffusion tensor imaging (DTI), and the relationships between abnormalities revealed by these two modalities remain unclear. This study applied a multimodal data-fusion method, known as linked independent component analysis (ICA), to a set of structural MRI and DTI data acquired from 46 adult males with ASD and 46 matched controls in order to elucidate associations between different aspects of atypical neuroanatomy of ASD. Linked ICA identified two composite components that showed significant between-group differences, one of which was significantly correlated with age. In the other component, participants with ASD showed decreased gray matter (GM) volumes in multiple regions, including the bilateral fusiform gyri, bilateral orbitofrontal cortices, and bilateral pre- and post-central gyri. These GM changes were linked with a pattern of decreased fractional anisotropy (FA) in several white matter tracts, such as the bilateral inferior longitudinal fasciculi, bilateral inferior fronto-occipital fasciculi, and bilateral corticospinal tracts. Furthermore, unimodal analysis for DTI data revealed significant reductions of FA along with increased mean diffusivity in those tracts for ASD, providing further evidence of disrupted anatomical connectivity. Taken together, our findings suggest that, in ASD, alterations in different aspects of brain morphology may co-occur in specific brain networks, providing a comprehensive view for understanding the neuroanatomy of this disorder.

•

Structural alterations of gray (GM) and white matter (WM) in ASD were investigated.

•

Linked independent component analysis was used for multimodal data analysis.

•

Alterations of GM and WM in ASD co-occurred in cognitive and affective networks.

•

Results reveal an integrative view of multiple aspects of structural changes in ASD.

Improving the accuracy of PGSE DTI experiments using the spatial distribution of b matrix

A novel method for improving the accuracy of diffusion tensor imaging (DTI) is proposed. It takes into account the b matrix spatial variations, which can be easily determined using a simple anisotropic diffusion phantom. In opposite to standard numerical procedure of the b matrix calculation that requires the exact knowledge of amplitudes, shapes and time dependencies of diffusion gradients, the new method, which we call BSD-DTI (B-matrix spatial distribution in DTI), relies on direct measurements of its space-dependent components. The proposed technique was demonstrated on the Bruker Biospec 94/20USR system, using the spin echo diffusion sequence to image an isotropic water phantom and an anisotropic capillary phantom. The accuracy of the diffusion tensor determination was improved by an overall factor of about 8 for the isotropic water phantom.

Diffusion tensor imaging in Alzheimer's disease: insights into the limbic-diencephalic network and methodological considerations

Glucose hypometabolism and gray matter atrophy are well known consequences of Alzheimer's disease (AD). Studies using these measures have shown that the earliest clinical stages, in which memory impairment is a relatively isolated feature, are associated with degeneration in an apparently remote group of areas—mesial temporal lobe (MTL), diencephalic structures such as anterior thalamus and mammillary bodies, and posterior cingulate. These sites are thought to be strongly anatomically inter-connected via a limbic-diencephalic network. Diffusion tensor imaging or DTI—an imaging technique capable of probing white matter tissue microstructure—has recently confirmed degeneration of the white matter connections of the limbic-diencephalic network in AD by way of an unbiased analysis strategy known as tract-based spatial statistics (TBSS). The present review contextualizes the relevance of these findings, in which the fornix is likely to play a fundamental role in linking MTL and diencephalon. An interesting by-product of this work has been in showing that alterations in diffusion behavior are complex in AD—while early studies tended to focus on fractional anisotropy, recent work has highlighted that this measure is not the most sensitive to early changes. Finally, this review will discuss in detail several technical aspects of DTI both in terms of image acquisition and TBSS analysis as both of these factors have important implications to ensure reliable observations are made that inform understanding of neurodegenerative diseases.

Diffusion tensor imaging of the sural nerve in normal controls

OBJECTIVE

To develop a diffusion tensor imaging (DTI) protocol for assessing the sural nerve in healthy subjects.

METHODS

Sural nerves in 25 controls were imaged using DTI at 3T with 6, 15, and 32 gradient directions. Fractional anisotropy (FA) and apparent diffusion coefficient (ADC) were computed from nerve regions of interest co-registered with T2-weighted images.

RESULTS

Coronal images with 0.5(RL)× 2.0(FH)× 0.5(AP)mm(3) resolution successfully localized the sural nerve. FA maps showed less variability with 32 directions (0.559 ± 0.071) compared to 15(0.590 ± 0.080) and 6(0.659 ± 0.109).

CONCLUSIONS

Our DTI protocol was effective in imaging sural nerves in controls to establish normative FA/ADC, with potential to be used non-invasively in diseased nerves of patients.

Reproducibility of diffusion tensor imaging in normal subjects: an evaluation of different gradient sampling schemes and registration algorithm

INTRODUCTION

Diffusion tensor imaging (DTI) is very useful for investigating white matter integrity in ageing and neurological disorders; thus, evaluating its reproducibility under different acquisition protocols and analysis methods may assist in the design of clinical studies.

METHODS

To measure the reproducibility of DTI in normal subjects, this study include (1) depicting the reproducibility of DTI measurements in commonly used regions-of-interest analysis by intraclass correlation coefficient (ICC) and coefficient of variation (CV), (2) evaluating and comparing inter and intrasession test-retest reproducibility, and (3) illustrating the effect of the number of diffusion-encoding directions (NDED) and registration algorithms on measurement reproducibility.

RESULTS

DTI measurements exhibit high reproducibility, with overall (430/480) ICC ≥ 0.70, (478/480) within-subject CV (CVws) ≤10.00 % and between-subject CV (CVbs) ranging from 1.32 to 13.63 %. Repeated measures ANOVAs and paired t tests were conducted to compare inter and intrasession reproducibility with different diffusion sampling schemes and registration algorithms. Our results also confirmed that increasing the NDED could improve the accuracy and reproducibility of DTI measurements. In addition, we compared reproducibility indices that were derived using different registration algorithms, and a tensor-based deformable registration yielded the most reproducible results. Finally, we found that increasing the NDED could reduce the difference between the reproducibility of measurement derived using different registration algorithms and between the reproducibility of intersession and intrasession.

CONCLUSION

Our results suggest that the choice of DTI acquisition protocol and post-processing methods can influence the accurate estimation and reproducibility of DTI measurements and should be considered carefully for clinical applications.

Development of a high angular resolution diffusion imaging human brain template

Brain diffusion templates contain rich information about the microstructure of the brain, and are used as references in spatial normalization or in the development of brain atlases. The accuracy of diffusion templates constructed based on the diffusion tensor (DT) model is limited in regions with complex neuronal micro-architecture. High angular resolution diffusion imaging (HARDI) overcomes limitations of the DT model and is capable of resolving intravoxel heterogeneity. However, when HARDI is combined with multiple-shot sequences to minimize image artifacts, the scan time becomes inappropriate for human brain imaging. In this work, an artifact-free HARDI template of the human brain was developed from low angular resolution multiple-shot diffusion data. The resulting HARDI template was produced in ICBM-152 space based on Turboprop diffusion data, was shown to resolve complex neuronal micro-architecture in regions with intravoxel heterogeneity, and contained fiber orientation information consistent with known human brain anatomy.

Determination of the appropriate b value and number of gradient directions for high-angular-resolution diffusion-weighted imaging

High-angular-resolution diffusion-weighted imaging (HARDI) is one of the most common MRI acquisition schemes for use with higher order models of diffusion. However, the optimal b value and number of diffusion-weighted (DW) directions for HARDI are still undetermined, primarily as a result of the large number of available reconstruction methods and corresponding parameters, making it impossible to identify a single criterion by which to assess performance. In this study, we estimate the minimum number of DW directions and optimal b values required for HARDI by focusing on the angular frequency content of the DW signal itself. The spherical harmonic (SH) series provides the spherical analogue of the Fourier series, and can hence be used to examine the angular frequency content of the DW signal. Using high-quality data acquired along 500 directions over a range of b values, we estimate that SH terms above l = 8 are negligible in practice for b values up to 5000 s/mm(2), implying that a minimum of 45 DW directions is sufficient to fully characterise the DW signal. l > 0 SH terms were found to increase as a function of b value, levelling off at b = 3000 s/mm(2), suggesting that this value already provides the highest achievable angular resolution. In practice, it is recommended to acquire more than the minimum of 45 DW directions to avoid issues with imperfections in the uniformity of the DW gradient directions and to meet signal-to-noise requirements of the intended reconstruction method.

DTI of human skeletal muscle: the effects of diffusion encoding parameters, signal-to-noise ratio and T2 on tensor indices and fiber tracts

In this study, we have performed simulations to address the effects of diffusion encoding parameters, signal-to-noise ratio (SNR) and T2 on skeletal muscle diffusion tensor indices and fiber tracts. Where appropriate, simulations were corroborated and validated by in vivo diffusion tensor imaging (DTI) of human skeletal muscle. Specifically, we have addressed: (i) the accuracy and precision of the diffusion parameters and eigenvectors at different SNR levels; (ii) the effects of the diffusion gradient direction encoding scheme; (iii) the optimal b value for diffusion tensor estimation; (iv) the effects of changes in skeletal muscle T2; and, finally, the influence of SNR on fiber tractography and derived (v) fiber lengths, (vi) pennation angles and (vii) fiber curvatures. We conclude that accurate DTI of skeletal muscle requires an SNR of at least 25, a b value of between 400 and 500 s/mm(2), and data acquired with at least 12 diffusion gradient directions homogeneously distributed on half a sphere. Furthermore, for DTI studies focusing on skeletal muscle injury or pathology, apparent changes in the diffusion parameters need to be interpreted with great care in view of the confounding effects of T2, particularly for moderate to low SNR values.

A connectome-based comparison of diffusion MRI schemes

Diffusion MRI has evolved towards an important clinical diagnostic and research tool. Though clinical routine is using mainly diffusion weighted and tensor imaging approaches, Q-ball imaging and diffusion spectrum imaging techniques have become more widely available. They are frequently used in research-oriented investigations in particular those aiming at measuring brain network connectivity. In this work, we aim at assessing the dependency of connectivity measurements on various diffusion encoding schemes in combination with appropriate data modeling. We process and compare the structural connection matrices computed from several diffusion encoding schemes, including diffusion tensor imaging, q-ball imaging and high angular resolution schemes, such as diffusion spectrum imaging with a publically available processing pipeline for data reconstruction, tracking and visualization of diffusion MR imaging. The results indicate that the high angular resolution schemes maximize the number of obtained connections when applying identical processing strategies to the different diffusion schemes. Compared to the conventional diffusion tensor imaging, the added connectivity is mainly found for pathways in the 50-100mm range, corresponding to neighboring association fibers and long-range associative, striatal and commissural fiber pathways. The analysis of the major associative fiber tracts of the brain reveals striking differences between the applied diffusion schemes. More complex data modeling techniques (beyond tensor model) are recommended 1) if the tracts of interest run through large fiber crossings such as the centrum semi-ovale, or 2) if non-dominant fiber populations, e.g. the neighboring association fibers are the subject of investigation. An important finding of the study is that since the ground truth sensitivity and specificity is not known, the comparability between results arising from different strategies in data reconstruction and/or tracking becomes implausible to understand.

Design of multishell sampling schemes with uniform coverage in diffusion MRI

PURPOSE

In diffusion MRI, a technique known as diffusion spectrum imaging reconstructs the propagator with a discrete Fourier transform, from a Cartesian sampling of the diffusion signal. Alternatively, it is possible to directly reconstruct the orientation distribution function in q-ball imaging, providing so-called high angular resolution diffusion imaging. In between these two techniques, acquisitions on several spheres in q-space offer an interesting trade-off between the angular resolution and the radial information gathered in diffusion MRI. A careful design is central in the success of multishell acquisition and reconstruction techniques.

METHODS

The design of acquisition in multishell is still an open and active field of research, however. In this work, we provide a general method to design multishell acquisition with uniform angular coverage. This method is based on a generalization of electrostatic repulsion to multishell.

RESULTS

We evaluate the impact of our method using simulations, on the angular resolution in one and two bundles of fiber configurations. Compared to more commonly used radial sampling, we show that our method improves the angular resolution, as well as fiber crossing discrimination.

DISCUSSION

We propose a novel method to design sampling schemes with optimal angular coverage and show the positive impact on angular resolution in diffusion MRI.

Motion Detection in Diffusion MRI via Online ODF Estimation

The acquisition of high angular resolution diffusion MRI is particularly long and subject motion can become an issue. The orientation distribution function (ODF) can be reconstructed online incrementally from diffusion-weighted MRI with a Kalman filtering framework. This online reconstruction provides real-time feedback throughout the acquisition process. In this article, the Kalman filter is first adapted to the reconstruction of the ODF in constant solid angle. Then, a method called STAR (STatistical Analysis of Residuals) is presented and applied to the online detection of motion in high angular resolution diffusion images. Compared to existing techniques, this method is image based and is built on top of a Kalman filter. Therefore, it introduces no additional scan time and does not require additional hardware. The performance of STAR is tested on simulated and real data and compared to the classical generalized likelihood ratio test. Successful detection of small motion is reported (rotation under 2°) with no delay and robustness to noise.

Tracking cerebral white matter changes across the lifespan: insights from diffusion tensor imaging studies

Delineating the normal development of brain white matter (WM) over the human lifespan is crucial to improved understanding of underlying WM pathology in neuropsychiatric and neurological conditions. We review the extant literature concerning diffusion tensor imaging studies of brain WM development in healthy individuals available until October 2012, summarise trends of normal development of human brain WM and suggest possible future research directions. Temporally, brain WM maturation follows a curvilinear pattern with an increase in fractional anisotropy (FA) from newborn to adolescence, decelerating in adulthood till a plateau around mid-adulthood, and a more rapid decrease of FA from old age onwards. Spatially, brain WM tracts develop from central to peripheral regions, with evidence of anterior-to-posterior maturation in commissural and projection fibres. The corpus callosum and fornix develop first and decline earlier, whilst fronto-temporal WM tracts like cingulum and uncinate fasciculus have protracted maturation and decline later. Prefrontal WM is most vulnerable with greater age-related FA reduction compared with posterior WM. Future large scale studies adopting longitudinal design will better clarify human brain WM changes over time.

Diffusion MR imaging: How to get the maximum from the experimental time

Diffusion-based MR imaging is the only non-invasive method for characterising the microstructural organization of brain tissue in vivo. Diffusion tensor MRI (DT-MRI) is currently routinely used in both research and clinical practice. However, other diffusion approaches are gaining more and more popularity and an increasing number of researchers express interest in using them concomitantly with DT-MRI. While non tensor-based methods hold great promises for increasing the specificity of diffusion MR imaging, including them in the experimental routine inevitably leads to longer experimental times. In most cases, this may preclude the translation of the full protocol to clinical practice, especially when these methods are to be used with subjects that are not compatible with long scanning sessions (e.g., with elderly and pediatric subjects who have difficulties in maintaining a fixed head position during a long imaging session).

The aim of this review is to guide the end-users on obtaining the maximum from the experimental time allocated to collecting diffusion MRI data. This is done by: (i) briefly reviewing non tensor-based approaches; (ii) reviewing the optimal protocols for both tensor and non tensor-based imaging; and (iii) drawing the conclusions for different experimental times.

Optimal real-time estimation in diffusion tensor imaging

Diffusion tensor imaging (DTI) constitutes the most used paradigm among the diffusion-weighted magnetic resonance imaging (DW-MRI) techniques due to its simplicity and application potential. Recently, real-time estimation in DW-MRI has deserved special attention, with several proposals aiming at the estimation of meaningful diffusion parameters during the repetition time of the acquisition sequence. Specifically focusing on DTI, the underlying model of the noise present in the acquired data is not taken into account, leading to a suboptimal estimation of the diffusion tensor. In this paper, we propose an optimal real-time estimation framework for DTI reconstruction in single-coil acquisitions. By including an online estimation of the time-changing noise variance associated to the acquisition process, the proposed method achieves the sequential best linear unbiased estimator. Results on both synthetic and real data show that our method outperforms those so far proposed, reaching the best performance of the existing proposals by processing a substantially lower number of diffusion images.