Power spectrum detects corpus callosum directionality using T2-weighted MRI in secondary progressive MS patients and controls
Shrushrita Sharma1 and Yunyan Zhang2

1Biomedical Engineering Program, University of Calgary, Calgary, AB, Canada, 2Departments of Radiology and Clinical Neurosciences, University of Calgary, Calgary, AB, Canada

Synopsis

Standard MRI is routinely collected in patient care but is limited in assessing changes in tissue microstructure. We developed a new method to assess tissue directionality using the power spectrum of T2-weighted MRI and validated it using the highly coherent structure, corpus callosum. In controls, power spectrum-derived angles corresponded exactly with the predicted aligning directions of the corpus callosum, and such aligning patterns were interrupted in advanced MS patients with increased variability and angular entropy. Fourier-based power spectrum may provide advanced measures of tissue directionality following myelin and axonal pathology using clinical scans.

PURPOSE

To determine whether analysis of the power spectrum generated from Fourier transform (FT) can detect white matter alignment using clinical scans and to evaluate how tissue coherency changes in patients with advanced multiple sclerosis (MS) using this new method.

METHODS

Brain MR images were acquired using a 3T scanner (GE Healthcare, DISCOVERY MR750, Milwaukee, USA) from a clinical study focusing on the structure and function of corpus callosum. Recruitment included 8 Secondary Progressive MS patients, whose disability score was all equal or above 6 out of 10, and 10 age- and gender-matched healthy controls. MRI protocol included a T2-weighted MR sequence with the following parameters: TE/TR = 82/60 field of view = 240 X 240 mm2, matrix size = 256 X 256. Recently we developed a new technique to evaluate the coherency and anisotropy of tissue structure based on the power spectrum of the FT. This method involves the following steps: 1) Calculating the FT and power spectra of an image; 2) performing logarithmic normalization and intensity thresholding of the power spectrum to preserve information within the upper 80 to 100 percentile; 3) conducting polar conversion of the processed power spectrum and plotting the orientation profile of the image (Fig. 1). To verify the utility of our method, we examined the directionality of corpus callosum as it is the largest interhemispheric white matter structure with highly organized nerve fibers. Regions of interest (ROIs) were drawn in the left, centra and right aspects of both genu and splenium of the structure in axial T2-weighted MRI< with an optimized size of 36 to 49 pixels each (Fig. 2). The distribution of aligning directions (degrees) in each ROI was plotted as a histogram, from which the angular entropy and dominant orientation of individual ROIs were derived and that were used to compare with the visually observed directions of the regions in healthy and to detect changes in anisotropy in patients.

RESULTS

Visually, the estimated directions of alignment were approximately 135 (right), 0 (center) and 45 (left) degrees in the genu and 45, 0 and 135 degrees in the splenium of corpus callosum in the image domain (Fig. 2). Using our assessing method, we found that the highest peak of power spectra in the right, center and left areas of genu was at 45, 90 and 135 degrees, which was equal to 135, 0 and 45 degrees respectively in the image domain at the same region. Similarly, in the right, center and left ROIs of the splenium, the spectral peak was at 135, 90 and 45 degrees, equivalent to 45, 0 and 135 degrees in the image (Fig. 3). In patients, corresponding regions showed a similar aligning profile but with reduced orientation strength and increased complexity as measured by angular entropy, suggesting loss of tissue or alignment (p<0.05 as per multi-comparison statistical analysis; Fig. 4). This is most severe in the left and right regions than in the central regions. All ROIs showed a similar trend; however, there was larger variability in the central regions of the genu and the splenium than the other areas.

DISCUSSION

Alterations in tissue microstructures can lead to changes in signal intensity and alignment of image voxels [1]. Assessing tissue anisotropy has been conducted mostly using advanced MRI such as diffusion-weighted imaging, which are not routinely acquired in a clinical setting [2]. In the study, we showed the possibility of obtaining similar structural information using standard T2-weighted MRI. The orientational property in different areas of the corpus callosum calculated using our method showed an extremely high match with the visualized angles of alignment in these regions. Moreover, the loss of white matter anisotropy known in patients with SPMS in detected using this method, even in the lesion free structure such as the corpus callosum [3]. While central ROIs demonstrate relativity large variability, likely due to the position and size of the ROI and needs further confirmation, our overall findings suggest the sensitivity and potential of this image post-processing method. The underlying pathologies of such MRI results may include demyelination, axonal loss and gliosis [4]. Future studies seek to validate current results in large sample size and compared MRI metrics with tissue histopathology.

CONCLUSION

This study suggests the potential of using FT power spectrum in assessing tissue directionality using clinical MRI. This may become an alternative method for determining white matter coherency and can have important clinical implications in improving the utility of clinical scans and in determining injury or repair in patients with MS and other neurological disorder.

Acknowledgements

We gratefully acknowledge the funding from University of Calgary Queen Elizabeth II Scholarship, NSERC I3T CREATE Scholarship and MS Society of Canada for their support.

References

[1] Ge Y, Multiple Sclerosis: The Role of MR Imaging. Am J Neuroradiol. 2006; 27(6):1165-76.

[2] Bayan C, Levitt JM, Miller E, Kaplan D, Georgekoudi I. Fully automated, quantitative, noninvasive assessment of collagen fiber content and organization in thick collagen gels. J Appl Phys. 2009; 105(10):1-11.

[3] Begin S, Belanger E, Laffray S, Aube B, Chamma E, Belise J, et al. Local assessment of myelin health in a multiple sclerosis mouse model using 2D Fourier transform approach. Biomed Opt Express. 2013; 4(10):2003-14.

[4] Kim A, Lakshman N, Petrol WM. Quantitative assessment of local collagen matrix remodeling in 3-D culture: The role of Rho kinase. Exp Cell Res. 2006; 312(18):3683-92.

Figures

Method demonstration. A: a sample T2-weighted brain MRI with a region of interest (ROI) in right genu(135° angle). B & C: zoomed view and Fourier Transform of the ROI. D-F: Power spectrum after logarithmic normalization (D) and thresholding (80-100 percentile, E), and angular distribution of the ROI (45° dominance, F).

Image (A) shows a T2-weighted MR image with regions of interest (ROIs) highlighted in the genu (1:right; 2:center; 3:left) and the splenium (1:right; 2:center; 3:left). Panel (B) shows the zoomed views of these ROIs in (A) showing their alignment.

Orientation plots of representative corpus callosum ROIs derived from the power spectrum. The degrees represent alignment directions in the ROI, and peak height reflects the strength of an alignment. These aligning degrees derived from the frequency domain are perpendicular to those seen in image domain based on the Fourier theory.

Plots of dominant alignment and angular entropy of all genu (G) and splenium (S) areas in patient and controls. The orientation of the filled lines in the bar represents the dominant aligning direction of the area including 45 ( right-G, left-S); 90 (center-G, center-S), and 135 degrees (right-S, left-S).



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
1924