Apparent Diffusion Coefficient Approximation and Diffusion Anisotropy Characterization in DWI

Characterizing brain structures and remodeling after TBI based on information content, diffusion entropy

BACKGROUND

To overcome the limitations of conventional diffusion tensor magnetic resonance imaging resulting from the assumption of a Gaussian diffusion model for characterizing voxels containing multiple axonal orientations, Shannon's entropy was employed to evaluate white matter structure in human brain and in brain remodeling after traumatic brain injury (TBI) in a rat.

METHODS

Thirteen healthy subjects were investigated using a Q-ball based DTI data sampling scheme. FA and entropy values were measured in white matter bundles, white matter fiber crossing areas, different gray matter (GM) regions and cerebrospinal fluid (CSF). Axonal densities' from the same regions of interest (ROIs) were evaluated in Bielschowsky and Luxol fast blue stained autopsy (n = 30) brain sections by light microscopy. As a case demonstration, a Wistar rat subjected to TBI and treated with bone marrow stromal cells (MSC) 1 week after TBI was employed to illustrate the superior ability of entropy over FA in detecting reorganized crossing axonal bundles as confirmed by histological analysis with Bielschowsky and Luxol fast blue staining.

RESULTS

Unlike FA, entropy was less affected by axonal orientation and more affected by axonal density. A significant agreement (r = 0.91) was detected between entropy values from in vivo human brain and histologically measured axonal density from post mortum from the same brain structures. The MSC treated TBI rat demonstrated that the entropy approach is superior to FA in detecting axonal remodeling after injury. Compared with FA, entropy detected new axonal remodeling regions with crossing axons, confirmed with immunohistological staining.

CONCLUSIONS

Entropy measurement is more effective in distinguishing axonal remodeling after injury, when compared with FA. Entropy is also more sensitive to axonal density than axonal orientation, and thus may provide a more accurate reflection of axonal changes that occur in neurological injury and disease.

Opportunities and pitfalls in the quantification of fiber integrity: what can we gain from Q-ball imaging?

The quantification of fiber integrity is central to the clinical application of diffusion imaging. Compared to diffusion tensor imaging (DTI), Q-ball imaging (QBI) allows for the depiction of multiple fiber directions within a voxel. However, this advantage has not yet been shown to translate directly to superior quantification of fiber integrity. Furthermore, recent developments in QBI reconstruction with solid angle consideration have led to sharper and intrinsically normalized orientation distribution functions. The implications of this technique on quantification are also unknown. To investigate this, the generalized fractional anisotropy (GFA) from the original and the more recent QBI reconstruction scheme and the DTI derived fractional anisotropy (FA) were evaluated comparatively using Monte Carlo simulations and real MRI measurements of crossing fiber phantoms. Contrast-to-noise ratio, accuracy, independence of the acquisition setup and the relation of single fiber anisotropies to measured anisotropy in crossings were assessed. In homogeneous single-fiber regions at b-values around 1000 s/mm2, the FA performed best. While the original QBI reconstruction does not show a clear advantage even at higher b-values and in crossing regions, the new reconstruction scheme yields superior properties and is recommended for quantification at higher b-values and especially in regions of heterogeneous fiber configuration.

Information-theoretic analysis of brain white matter fiber orientation distribution functions

We propose a new information-theoretic metric, the symmetric Kullback-Leibler divergence (sKL-divergence), to measure the difference between two water diffusivity profiles in high angular resolution diffusion imaging (HARDI). Water diffusivity profiles are modeled as probability density functions on the unit sphere, and the sKL-divergence is computed from a spherical harmonic series, which greatly reduces computational complexity. Adjustment of the orientation of diffusivity functions is essential when the image is being warped, so we propose a fast algorithm to determine the principal direction of diffusivity functions using principal component analysis (PCA). We compare sKL-divergence with other inner-product based cost functions using synthetic samples and real HARDI data, and show that the sKL-divergence is highly sensitive in detecting small differences between two diffusivity profiles and therefore shows promise for applications in the nonlinear registration and multisubject statistical analysis of HARDI data.

Automatic identification and classification of characteristic kinetic curves of breast lesions on DCE-MRI

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) of the breast is being used increasingly in the detection and diagnosis of breast cancer as a complementary modality to mammography and sonography. Although the potential diagnostic value of kinetic curves in DCE-MRI is established, the method for generating kinetic curves is not standardized. The inherent reason that curve identification is needed is that the uptake of contrast agent in a breast lesion is often heterogeneous, especially in malignant lesions. It is accepted that manual region of interest selection in 4D breast magnetic resonance (MR) images to generate the kinetic curve is a time-consuming process and suffers from significant inter- and intraobserver variability. We investigated and developed a fuzzy c-means (FCM) clustering-based technique for automatically identifying characteristic kinetic curves from breast lesions in DCE-MRI of the breast. Dynamic contrast-enhanced MR images were obtained using a T1-weighted 3D spoiled gradient echo sequence with Gd-DTPA dose of 0.2 mmol/kg and temporal resolution of 69 s. FCM clustering was applied to automatically partition the signal-time curves in a segmented 3D breast lesion into a number of classes (i.e., prototypic curves). The prototypic curve with the highest initial enhancement was selected as the representative characteristic kinetic curve (CKC) of the lesion. Four features were then extracted from each characteristic kinetic curve to depict the maximum contrast enhancement, time to peak, uptake rate, and washout rate of the lesion kinetics. The performance of the kinetic features in the task of distinguishing between benign and malignant lesions was assessed by receiver operating characteristic analysis. With a database of 121 breast lesions (77 malignant and 44 benign cases), the classification performance of the FCM-identified CKCs was found to be better than that from the curves obtained by averaging over the entire lesion and similar to kinetic curves generated from regions drawn within the lesion by a radiologist experienced in breast MRI.