Efficient gradient calibration based on diffusion MRI

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

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

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

PURPOSE

To use diffusion measurements to map the spatial dependence of the magnetic field produced by the gradient coils of an MRI scanner with sufficient accuracy to correct errors in quantitative diffusion MRI (DMRI) caused by gradient nonlinearity and gradient amplifier miscalibration.

THEORY AND METHODS

The field produced by the gradient coils is expanded in regular solid harmonics. The expansion coefficients are found by fitting a model to a minimum set of diffusion-weighted images of an isotropic diffusion phantom. The accuracy of the resulting gradient coil field maps is evaluated by using them to compute corrected b-matrices that are then used to process a multi-shell diffusion tensor imaging (DTI) dataset with 32 diffusion directions per shell.

RESULTS

The method substantially reduces both the spatial inhomogeneity of the computed mean diffusivities (MD) and the computed values of the fractional anisotropy (FA), as well as virtually eliminating any artifactual directional bias in the tensor field secondary to gradient nonlinearity. When a small scaling miscalibration was purposely introduced in the x, y, and z, the method accurately detected the amount of miscalibration on each gradient axis.

CONCLUSION

The method presented detects and corrects the effects of gradient nonlinearity and gradient gain miscalibration using a simple isotropic diffusion phantom. The correction would improve the accuracy of DMRI measurements in the brain and other organs for both DTI and higher order diffusion analysis. In particular, it would allow calibration of MRI systems, improving data harmony in multicenter studies.

Regional variations in ex-vivo diffusion tensor anisotropy are associated with cardiomyocyte remodeling in rats after left ventricular pressure overload

BACKGROUND

Pressure overload left ventricular (LV) hypertrophy is characterized by increased cardiomyocyte width and ventricle wall thickness, however the regional variation of this remodeling is unclear. Cardiovascular magnetic resonance (CMR) diffusion tensor imaging (DTI) may provide a non-invasive, comprehensive, and geometrically accurate method to detect regional differences in structural remodeling in hypertrophy. We hypothesized that DTI parameters, such as fractional and planar anisotropy, would reflect myocyte remodeling due to pressure overload in a regionally-dependent manner.

METHODS

We investigated the regional distributions of myocyte remodeling in rats with or without transverse aortic constriction (TAC) via direct measurement of myocyte dimensions with confocal imaging of thick tissue sections, and correlated myocyte cross-sectional area and other geometric features with parameters of diffusivity from ex-vivo DTI in the same regions of the same hearts.

RESULTS

We observed regional differences in several parameters from DTI between TAC hearts and SHAM controls. Consistent with previous studies, helix angles from DTI correlated strongly with those measured directly from histological sections (p < 0.001, R2 = 0.71). There was a transmural gradient in myocyte cross-sectional area in SHAM hearts that was diminished in the TAC group. We also found several regions of significantly altered DTI parameters in TAC LV compared to SHAM, especially in myocyte sheet angle dispersion and planar anisotropy. Among others, these parameters correlated significantly with directly measured myocyte aspect ratios.

CONCLUSIONS

These results show that structural remodeling in pressure overload LV hypertrophy is regionally heterogeneous, especially transmurally, with a greater degree of remodeling in the sub-endocardium compared to the sub-epicardium. Additionally, several parameters derived from DTI correlated significantly with measurements of myocyte geometry from direct measurement in histological sections. We suggest that DTI may provide a non-invasive, comprehensive method to detect regional structural myocyte LV remodeling during disease.

Assessment of the systematic errors caused by diffusion gradient inhomogeneity in DTI-computer simulations

Diffusion tensor imaging (DTI) is a powerful MRI modality that allows the investigation of the microstructure of tissues both in vivo and noninvasively. Its reliability is strictly dependent on the performance of diffusion-sensitizing gradients, of which spatial nonuniformity is a known issue in the case of virtually all clinical MRI scanners. The influence of diffusion gradient inhomogeneity on the accuracy of the diffusion tensor imaging was investigated by means of computer simulations supported by an MRI experiment performed at the isocenter and 15 cm away. The DTI measurements of two diffusion phantoms were simulated assuming a nonuniform diffusion-sensitizing gradient and various levels of noise. Thereafter, the tensors were calculated by two methods: (i) assuming a spatially constant b-matrix (standard DTI) and (ii) applying the b-matrix spatial distribution in the DTI (BSD-DTI) technique, a method of indicating the b-matrix for each voxel separately using an anisotropic phantom as a standard of diffusion. The average eigenvalues and fractional anisotropy across the homogeneous region of interest were calculated and compared with the expected values. Diffusion gradient inhomogeneity leads to overestimation of the largest eigenvalue, underestimation of the smallest one and thus overestimation of fractional anisotropy. The effect is similar to that caused by noise; however, it could not be corrected by increasing SNR. The MRI measurements, performed using a 3 T clinical scanner, revealed that the split of the eigenvalues measured 15 cm away from the isocenter is significant (up to 25%). The BSD-DTI calibration allowed the reduction of the measured fractional anisotropy of the isotropic medium from 0.174 to 0.031, suggesting that gradient inhomogeneity was the main cause of this error. For the phantom measured at the isocenter, however, the split was almost not observed; the average eigenvalues were shifted from the expected value by ~ 5%.

Multicenter Evaluation of Geometric Accuracy of MRI Protocols Used in Experimental Stroke

It has recently been suggested that multicenter preclinical stroke studies should be carried out to improve translation from bench to bedside, but the accuracy of magnetic resonance imaging (MRI) scanners routinely used in experimental stroke has not yet been evaluated. We aimed to assess and compare geometric accuracy of preclinical scanners and examine the longitudinal stability of one scanner using a simple quality assurance (QA) protocol. Six 7 Tesla animal scanners across six different preclinical imaging centers throughout Europe were used to scan a small structural phantom and estimate linear scaling errors in all orthogonal directions and volumetric errors. Between-scanner imaging consisted of a standard sequence and each center's preferred sequence for the assessment of infarct size in rat models of stroke. The standard sequence was also used to evaluate the drift in accuracy of the worst performing scanner over a period of six months following basic gradient calibration. Scaling and volumetric errors using the standard sequence were less variable than corresponding errors using different stroke sequences. The errors for one scanner, estimated using the standard sequence, were very high (above 4% scaling errors for each orthogonal direction, 18.73% volumetric error). Calibration of the gradient coils in this system reduced scaling errors to within ±1.0%; these remained stable during the subsequent 6-month assessment. In conclusion, despite decades of use in experimental studies, preclinical MRI still suffers from poor and variable geometric accuracy, influenced by the use of miscalibrated systems and various types of sequences for the same purpose. For effective pooling of data in multicenter studies, centers should adopt standardized procedures for system QA and in vivo imaging.

Resolving Fine Cardiac Structures in Rats with High-Resolution Diffusion Tensor Imaging

Cardiac architecture is fundamental to cardiac function and can be assessed non-invasively with diffusion tensor imaging (DTI). Here, we aimed to overcome technical challenges in ex vivo DTI in order to extract fine anatomical details and to provide novel insights in the 3D structure of the heart. An integrated set of methods was implemented in ex vivo rat hearts, including dynamic receiver gain adjustment, gradient system scaling calibration, prospective adjustment of diffusion gradients, and interleaving of diffusion-weighted and non-diffusion-weighted scans. Together, these methods enhanced SNR and spatial resolution, minimised orientation bias in diffusion-weighting, and reduced temperature variation, enabling detection of tissue structures such as cell alignment in atria, valves and vessels at an unprecedented level of detail. Improved confidence in eigenvector reproducibility enabled tracking of myolaminar structures as a basis for segmentation of functional groups of cardiomyocytes. Ex vivo DTI facilitates acquisition of high quality structural data that complements readily available in vivo cardiac functional and anatomical MRI. The improvements presented here will facilitate next generation virtual models integrating micro-structural and electro-mechanical properties of the heart.

Mapping cardiac microstructure of rabbit heart in different mechanical states by high resolution diffusion tensor imaging: A proof-of-principle study

Myocardial microstructure and its macroscopic materialisation are fundamental to the function of the heart. Despite this importance, characterisation of cellular features at the organ level remains challenging, and a unifying description of the structure of the heart is still outstanding. Here, we optimised diffusion tensor imaging data to acquire high quality data in ex vivo rabbit hearts in slack and contractured states, approximating diastolic and systolic conditions. The data were analysed with a suite of methods that focused on different aspects of the myocardium. In the slack heart, we observed a similar transmural gradient in helix angle of the primary eigenvector of up to 23.6°/mm in the left ventricle and 24.2°/mm in the right ventricle. In the contractured heart, the same transmural gradient remained largely linear, but was offset by up to +49.9° in the left ventricle. In the right ventricle, there was an increase in the transmural gradient to 31.2°/mm and an offset of up to +39.0°. The application of tractography based on each eigenvector enabled visualisation of streamlines that depict cardiomyocyte and sheetlet organisation over large distances. We observed multiple V- and N-shaped sheetlet arrangements throughout the myocardium, and insertion of sheetlets at the intersection of the left and right ventricle. This study integrates several complementary techniques to visualise and quantify the heart's microstructure, projecting parameter representations across different length scales. This represents a step towards a more comprehensive characterisation of myocardial microstructure at the whole organ level.