Estimation of field inhomogeneity map following magnitude-based ambiguity-resolved water-fat separation.

Magnitude-based PDFF (Proton Density Fat Fraction) and R2∗ mapping with resolved water-fat ambiguity is extended to calculate field inhomogeneity (field map) using the phase images. The estimation is formulated in matrix form, resolving the field map in a least-squares sense. PDFF and R2∗ from magnitude fitting may be updated using the estimated field maps. The limits of quantification of our voxel-independent implementation were assessed. Bland-Altman was used to compare PDFF and field maps from our method against a reference complex-based method on 152 UK Biobank subjects (1.5 T Siemens). A separate acquisition (3 T Siemens) presenting field inhomogeneities was also used. The proposed field mapping was accurate beyond double the complex-based limit range. High agreement was obtained between the proposed method and the reference in UK. Robust field mapping was observed at 3 T, for inhomogeneities over 400 Hz including rapid variation across edges. Field mapping following unambiguous magnitude-based water-fat separation was demonstrated in-vivo and showed potential at 3 T.

Assessing Myocardial Microstructure With Biophysical Models of Diffusion MRI.

Biophysical models are a promising means for interpreting diffusion weighted magnetic resonance imaging (DW-MRI) data, as they can provide estimates of physiologically relevant parameters of microstructure including cell size, volume fraction, or dispersion. However, their application in cardiac microstructure mapping (CMM) has been limited. This study proposes seven new two-compartment models with combination of restricted cylinder models and a diffusion tensor to represent intra- and extracellular spaces, respectively. Three extended versions of the cylinder model are studied here: cylinder with elliptical cross section (ECS), cylinder with Gamma distributed radii (GDR), and cylinder with Bingham distributed axes (BDA). The proposed models were applied to data in two fixed mouse hearts, acquired with multiple diffusion times, q-shells and diffusion encoding directions. The cylinderGDR-pancake model provided the best performance in terms of root mean squared error (RMSE) reducing it by 25% compared to diffusion tensor imaging (DTI). The cylinderBDA-pancake model represented anatomical findings closest as it also allows for modelling dispersion. High-resolution 3D synchrotron X-ray imaging (SRI) data from the same specimen was utilized to evaluate the biophysical models. A novel tensor-based registration method is proposed to align SRI structure tensors to the MR diffusion tensors. The consistency between SRI and DW-MRI parameters demonstrates the potential of compartment models in assessing physiologically relevant parameters.

Improved compressed sensing and super-resolution of cardiac diffusion MRI with structure-guided total variation.

PURPOSE: Structure-guided total variation is a recently introduced prior that allows reconstruction of images using knowledge of the location and orientation of edges in a reference image. In this work, we demonstrate the advantages of a variant of structure-guided total variation known as directional total variation (DTV), over traditional total variation (TV), in the context of compressed-sensing reconstruction and super-resolution. METHODS: We compared TV and DTV in retrospectively undersampled ex vivo diffusion tensor imaging and diffusion spectrum imaging data from healthy, sham, and hypertrophic rat hearts. RESULTS: In compressed sensing at an undersampling factor of 8, the RMS error of mean diffusivity and fractional anisotropy relative to the fully sampled ground truth were 44% and 20% lower in DTV compared with TV. In super-resolution, these values were 29% and 14%, respectively. Similarly, we observed improvements in helix angle, transverse angle, sheetlet elevation, and sheetlet azimuth. The RMS error of the diffusion kurtosis in the undersampled data relative to the ground truth was uniformly lower (22% on average) with DTV compared to TV. CONCLUSION: Acquiring one fully sampled non-diffusion-weighted image and 10 diffusion-weighted images at 8× undersampling would result in an 80% net reduction in data needed. We demonstrate efficacy of the DTV algorithm over TV in reducing data sampling requirements, which can be translated into higher apparent resolution and potentially shorter scan times. This method would be equally applicable in diffusion MRI applications outside the heart.

Automatic lesion detection, segmentation and characterization via 3D multiscale morphological sifting in breast MRI.

Previous studies on computer aided detection/diagnosis (CAD) in 4D breast magnetic resonance imaging (MRI) usually regard lesion detection, segmentation and characterization as separate tasks, and typically require users to manually select 2D MRI slices or regions of interest as the input. In this work, we present a breast MRI CAD system that can handle 4D multimodal breast MRI data, and integrate lesion detection, segmentation and characterization with no user intervention. The proposed CAD system consists of three major stages: region candidate generation, feature extraction and region candidate classification. Breast lesions are firstly extracted as region candidates using the novel 3D multiscale morphological sifting (MMS). The 3D MMS, which uses linear structuring elements to extract lesion-like patterns, can segment lesions from breast images accurately and efficiently. Analytical features are then extracted from all available 4D multimodal breast MRI sequences, including T1-, T2-weighted and DCE sequences, to represent the signal intensity, texture, morphological and enhancement kinetic characteristics of the region candidates. The region candidates are lastly classified as lesion or normal tissue by the random under-sampling boost (RUSboost), and as malignant or benign lesion by the random forest. Evaluated on a breast MRI dataset which contains a total of 117 cases with 141 biopsy-proven lesions (95 malignant and 46 benign lesions), the proposed system achieves a true positive rate (TPR) of 0.90 at 3.19 false positives per patient (FPP) for lesion detection and a TPR of 0.91 at a FPP of 2.95 for identifying malignant lesions without any user intervention. The average dice similarity index (DSI) is0.72±0.15for lesion segmentation. Compared with previously proposed lesion detection, detection-segmentation and detection-characterization systems evaluated on the same breast MRI dataset, the proposed CAD system achieves a favourable performance in breast lesion detection and characterization.

The impact of signal-to-noise ratio, diffusion-weighted directions and image resolution in cardiac diffusion tensor imaging - insights from the ex-vivo rat heart.

BACKGROUND: Cardiac diffusion tensor imaging (DTI) is limited by scan time and signal-to-noise (SNR) restrictions. This invariably leads to a trade-off between the number of averages, diffusion-weighted directions (ND), and image resolution. Systematic evaluation of these parameters is therefore important for adoption of cardiac DTI in clinical routine where time is a key constraint. METHODS: High quality reference DTI data were acquired in five ex-vivo rat hearts. We then retrospectively set 2 ≤ SNR ≤ 97, 7 ≤ ND ≤ 61, varied the voxel volume by up to 192-fold and investigated the impact on the accuracy and precision of commonly derived parameters. RESULTS: For maximal scan efficiency, the accuracy and precision of the mean diffusivity is optimised when SNR is maximised at the expense of ND. With typical parameter settings used clinically, we estimate that fractional anisotropy may be overestimated by up to 13% with an uncertainty of ±30%, while the precision of the sheetlet angles may be as poor as ±31°. Although the helix angle has better precision of ±14°, the transmural range of helix angles may be under-estimated by up to 30° in apical and basal slices, due to partial volume and tapering myocardial geometry. CONCLUSIONS: These findings inform a baseline of understanding upon which further issues inherent to in-vivo cardiac DTI, such as motion, strain and perfusion, can be considered. Furthermore, the reported bias and reproducibility provides a context in which to assess cardiac DTI biomarkers.

High-Resolution Magnetic Resonance Imaging of the Regenerating Adult Zebrafish Heart.

The adult zebrafish is a well-established model for studying heart regeneration, but due to its tissue opaqueness, repair has been primarily assessed using destructive histology, precluding repeated investigations of the same animal. We present a high-resolution, non-invasive in vivo magnetic resonance imaging (MRI) method incorporating a miniature respiratory and anaesthetic perfusion set-up for live adult zebrafish, allowing for visualization of scar formation and heart regeneration in the same animal over time at an isotropic 31 µm voxel resolution. To test the method, we compared well and poorly healing cardiac ventricles using a transgenic fish model that exhibits heat-shock (HS) inducible impaired heart regeneration. HS-treated groups revealed persistent scar tissue for 10 weeks, while control groups were healed after 4 weeks. Application of the advanced MRI technique allowed clear discrimination of levels of repair following cryo- and resection injury for several months. It further provides a novel tool for in vivo time-lapse imaging of adult fish for non-cardiac studies, as the method can be readily applied to image wound healing in other injured or diseased tissues, or to monitor tissue changes over time, thus expanding the range of questions that can be addressed in adult zebrafish and other small aquatic species.

Validation of diffusion tensor MRI measurements of cardiac microstructure with structure tensor synchrotron radiation imaging.

BACKGROUND: Diffusion tensor imaging (DTI) is widely used to assess tissue microstructure non-invasively. Cardiac DTI enables inference of cell and sheetlet orientations, which are altered under pathological conditions. However, DTI is affected by many factors, therefore robust validation is critical. Existing histological validation is intrinsically flawed, since it requires further tissue processing leading to sample distortion, is routinely limited in field-of-view and requires reconstruction of three-dimensional volumes from two-dimensional images. In contrast, synchrotron radiation imaging (SRI) data enables imaging of the heart in 3D without further preparation following DTI. The objective of the study was to validate DTI measurements based on structure tensor analysis of SRI data. METHODS: One isolated, fixed rat heart was imaged ex vivo with DTI and X-ray phase contrast SRI, and reconstructed at 100 µm and 3.6 µm isotropic resolution respectively. Structure tensors were determined from the SRI data and registered to the DTI data. RESULTS: Excellent agreement in helix angles (HA) and transverse angles (TA) was observed between the DTI and structure tensor synchrotron radiation imaging (STSRI) data, where HADTI-STSRI = -1.4° ± 23.2° and TADTI-STSRI = -1.4° ± 35.0° (mean ± 1.96 standard deviation across all voxels in the left ventricle). STSRI confirmed that the primary eigenvector of the diffusion tensor corresponds with the cardiomyocyte long-axis across the whole myocardium. CONCLUSIONS: We have used STSRI as a novel and high-resolution gold standard for the validation of DTI, allowing like-with-like comparison of three-dimensional tissue structures in the same intact heart free of distortion. This represents a critical step forward in independently verifying the structural basis and informing the interpretation of cardiac DTI data, thereby supporting the further development and adoption of DTI in structure-based electro-mechanical modelling and routine clinical applications.

Monte Carlo Simulations of Diffusion Weighted MRI in Myocardium: Validation and Sensitivity Analysis.

A model of cardiac microstructure and diffusion MRI is presented, and compared with experimental data from ex vivo rat hearts. The model includes a simplified representation of individual cells, with physiologically correct cell size and orientation, as well as intra- to extracellular volume ratio. Diffusion MRI is simulated using a Monte Carlo model and realistic MRI sequences. The results show good correspondence between the simulated and experimental MRI signals. Similar patterns are observed in the eigenvalues of the diffusion tensor, the mean diffusivity (MD), and the fractional anisotropy (FA). A sensitivity analysis shows that the diffusivity is the dominant influence on all three eigenvalues of the diffusion tensor, the MD, and the FA. The area and aspect ratio of the cell cross-section affect the secondary and tertiary eigenvalues, and hence the FA. Within biological norms, the cell length, volume fraction of cells, and rate of change of helix angle play a relatively small role in influencing tissue diffusion. Results suggest that the model could be used to improve understanding of the relationship between cardiac microstructure and diffusion MRI measurements, as well as in testing and refinement of cardiac diffusion MRI protocols.

Evaluation of non-Gaussian diffusion in cardiac MRI.

PURPOSE: The diffusion tensor model assumes Gaussian diffusion and is widely applied in cardiac diffusion MRI. However, diffusion in biological tissue deviates from a Gaussian profile as a result of hindrance and restriction from cell and tissue microstructure, and may be quantified better by non-Gaussian modeling. The aim of this study was to investigate non-Gaussian diffusion in healthy and hypertrophic hearts. METHODS: Thirteen rat hearts (five healthy, four sham, four hypertrophic) were imaged ex vivo. Diffusion-weighted images were acquired at b-values up to 10,000 s/mm2 . Models of diffusion were fit to the data and ranked based on the Akaike information criterion. RESULTS: The diffusion tensor was ranked best at b-values up to 2000 s/mm2 but reflected the signal poorly in the high b-value regime, in which the best model was a non-Gaussian "beta distribution" model. Although there was considerable overlap in apparent diffusivities between the healthy, sham, and hypertrophic hearts, diffusion kurtosis and skewness in the hypertrophic hearts were more than 20% higher in the sheetlet and sheetlet-normal directions. CONCLUSION: Non-Gaussian diffusion models have a higher sensitivity for the detection of hypertrophy compared with the Gaussian model. In particular, diffusion kurtosis may serve as a useful biomarker for characterization of disease and remodeling in the heart. Magn Reson Med 78:1174-1186, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

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.

Prospective acceleration of diffusion tensor imaging with compressed sensing using adaptive dictionaries.

PURPOSE: Diffusion MRI requires acquisition of multiple diffusion-weighted images, resulting in long scan times. Here, we investigate combining compressed sensing and a fast imaging sequence to dramatically reduce acquisition times in cardiac diffusion MRI. METHODS: Fully sampled and prospectively undersampled diffusion tensor imaging data were acquired in five rat hearts at acceleration factors of between two and six using a fast spin echo (FSE) sequence. Images were reconstructed using a compressed sensing framework, enforcing sparsity by means of decomposition by adaptive dictionaries. A tensor was fit to the reconstructed images and fiber tractography was performed. RESULTS: Acceleration factors of up to six were achieved, with a modest increase in root mean square error of mean apparent diffusion coefficient (ADC), fractional anisotropy (FA), and helix angle. At an acceleration factor of six, mean values of ADC and FA were within 2.5% and 5% of the ground truth, respectively. Marginal differences were observed in the fiber tracts. CONCLUSION: We developed a new k-space sampling strategy for acquiring prospectively undersampled diffusion-weighted data, and validated a novel compressed sensing reconstruction algorithm based on adaptive dictionaries. The k-space undersampling and FSE acquisition each reduced acquisition times by up to 6× and 8×, respectively, as compared to fully sampled spin echo imaging. Magn Reson Med 76:248-258, 2016. © 2015 Wiley Periodicals, Inc.

Fully automatic lesion segmentation in breast MRI using mean-shift and graph-cuts on a region adjacency graph.

PURPOSE: To present and evaluate a fully automatic method for segmentation (i.e., detection and delineation) of suspicious tissue in breast MRI. MATERIALS AND METHODS: The method, based on mean-shift clustering and graph-cuts on a region adjacency graph, was developed and its parameters tuned using multimodal (T1, T2, DCE-MRI) clinical breast MRI data from 35 subjects (training data). It was then tested using two data sets. Test set 1 comprises data for 85 subjects (93 lesions) acquired using the same protocol and scanner system used to acquire the training data. Test set 2 comprises data for eight subjects (nine lesions) acquired using a similar protocol but a different vendor's scanner system. Each lesion was manually delineated in three-dimensions by an experienced breast radiographer to establish segmentation ground truth. The regions of interest identified by the method were compared with the ground truth and the detection and delineation accuracies quantitatively evaluated. RESULTS: One hundred percent of the lesions were detected with a mean of 4.5 ± 1.2 false positives per subject. This false-positive rate is nearly 50% better than previously reported for a fully automatic breast lesion detection system. The median Dice coefficient for Test set 1 was 0.76 (interquartile range, 0.17), and 0.75 (interquartile range, 0.16) for Test set 2. CONCLUSION: The results demonstrate the efficacy and accuracy of the proposed method as well as its potential for direct application across different MRI systems. It is (to the authors' knowledge) the first fully automatic method for breast lesion detection and delineation in breast MRI.