Spinal cord 1H MR Spectroscopy (1H-MRS) is a promising method for musculoskeletal research. However, due to the spine’s anatomical location there is a significant degradation of signal quality due to magnetic field inhomogeneities, rendering most MRS approaches inaccurate. Although there has been measurement of ΔB0 in spinal cord MRS, there are no comprehensive assessments of temporal changes in B0 and B1+ relating physiological disturbances with MRS accuracy. Thus our goal was to continually measure temporal changes in B0 and B1+ during the length of a typical MEGA-PRESS scan (10min).
Experiments were performed using a 3T GE MR750 scanner and a 32 channel receive-only
head/neck array RF coil (General Electric Healthcare, Milwaukee, WI). Measurements were made on a static home
designed/built spinal cord phantom and a healthy human subject. In separate experiments B0 and B1+ field maps were repeatedly acquired over a 10 minute timespan in the phantom and
cervical spinal cord (C4) of a human subject (Fig.1). The B0 and B1+ sequences,
each taking 6s
and 11s respectively, were
auto-shimmed only for the first measurements so the same shim currents were
maintained for the duration of the 10 minutes. Acquisition of heart rate and respiration were
performed using a pulsed oximeter and thoracic respiratory bellows, both temporally
synchronized with the beginning of the first scan. These were continually recorded over the
length of the 10 minute scan at a sampling rate of 10Hz and 25Hz, respectively.
ROIs were selected in B0 and B1+ field maps of the spinal cord phantom and human C-spine (Fig.2) and means and standard deviations were calculated. There was slight variation in the B0 (mean±SD) 60.7±0.51Hz and B1+ 28.4±2.72 degrees/flip angle fields within the spinal cord, over the 10minutes of scanning. Whereas there was significant variation in the B0 (mean±SD) 48.5±2.72Hz and B1+ 21.4±4.61 degrees/flip angle fields within the spinal cord, over the 10minutes of scanning. Hence, physiological variation was more the dominant contribution to variability. The heart rate data was Fourier transformed (FT) and the power spectral density (PSD) was used to identify the dominant frequency contribution of motion for each field map. A standard deviation window for each field map was found for the respiratory motion contribution. The physiological data was further analyzed using principal component analysis (PCA) (Fig.3). The dominant principal component in B0 variability, heart rate, yielded a correlation coefficient of 0.71 whereas for B1+, the dominant principal component, respiratory motion, yielded a correlation coefficient of 0.52.
[1] Henning A., et al. (2008) Magn Reson Med 59(6):1250-1258.
[2] Cooke FJ. et al (2004) Magn Reson Med 51(6):1122-1128.
[3] Sacolick Li. et al. (2010) Magn Reson Med 63(5):1315-1322.