Synopsis
Temporal signal fluctuations (tSNR) and local effective echo time (TElocal) are explored and their influence on BOLD sensitivity is investigated at 7T for healthy subjects and tumor patients, where prominent spatial variations in those two measures are expected. We show that tSNR may indicate sufficient sensitivity to detect activation but that BOLD sensitivity may be dramatically reduced by changes in TElocal close to pathologies and vital brain functions (motor, speech, auditory). Neglecting local TE variations can thus lead to false negative results in clinical fMRI. We thus suggest a new BOLD sensitivity metric based on TElocal∙tSNR.Purpose
To explore spatial variations in temporal signal-to-noise ratio (tSNR) and echo time (TE
local) and develop the most reliable metric for BOLD sensitivity estimation.
Introduction
Estimating BS is useful for fMRI studies in selecting the sequence
(simultaneous-multi-slice, 3D or slice-by-slice EPI) and optimizing
parameters
1–3, but also in the assessment of the reliability of the
functional results (considering false negatives and positives)
4. BS is
commonly estimated using two different metrics: 1) the product of the
local effective echo time and the image signal (TE
local∙S)
2,5 or 2) tSNR
3,6. We explore the pitfalls of each metric, suggest TE
local∙tSNR
as a more reliable measure of BS and investigate all three approaches
at 7T in healthy volunteers and tumor patients, where prominent regional
variations are expected.
Theory
Gradients in B
0, induced by tissue-air susceptibility differences, introduce intravoxel dephasing and signal loss, but also local modifications of the effective echo time through shifts in the k-space signal
7. If TE
local is reduced by the gradients, image signal (S) increases, but BS decreases, as BOLD contrast has less time to develop. In such regions BS metric 2 (tSNR) overestimates BS. If, on the other hand, the TE
local is increased, signal is reduced, but BOLD contrast increases. Here tSNR underestimates BS unless TE
local is so long that the signal is not acquired. At this point, signal, tSNR and BS suddenly fall to zero (type 2 signal loss). Metric 1 (TE
local∙S) considers both signal and TE
local changes, but does not identify dynamic physiological and technical noise sources, which are included by substituting S by tSNR. We, thus, suggest using a third metric: TE
local∙tSNR as a more reliable BS measure.
Methods and analysis
Measurements were carried out with a 7T Siemens scanner and a 32-channel head coil. Field maps (multi-echo GE, TEs=[5,10,16]ms, TR=400ms, 1.7x1.7x3.0mm
3, bandwidth=540Hz/pix) and a 5 min long resting state scan (EPI, TE=22ms, TR=2000ms, 1.7x1.7x3.0mm
3, bandwidth=1447Hz/pix) were acquired for 3 volunteers and 11 tumor patients. Field maps were calculated using the Hermitian inner product
8. TE
local maps were derived from field gradient according to Ref.
9 and normalized to TE=22ms. Voxels with type 2 signal loss (TE>47ms) were marked white in TE
local maps. EPI runs were coregistered to GE scans, distortion corrected, and mean signal S and tSNR were calculated.
Results
Fig. 1 shows normalized TElocal and two BS maps calculated according to metric 1 and 2 for a representative healthy subject. BS metric 1 shows high values in CSF and gray matter and low values in white matter (Fig.1c), whereas metric 2 shows rather the opposite: higher values in homogeneous white matter (Fig.1d). This does not, however, include variations due to changes in TElocal (Fig.1b). Since metric 1 does not consider temporal variations in signal, as tSNR does, it was omitted from further consideration.
Variations in TElocal and the differences in BS metric 2 and 3 representative for healthy volunteers are presented in Fig. 2. In a ventral slice (number 3) tSNR values are highest (about 70) close to sinuses and ear canals (Fig.2e, see arrows), where TElocal is strongly reduced (Fig.2d, by about 50%). The reduction in BS due to short TE is apparent in TElocal∙tSNR maps (reducing BS values from 70 to 35 in Fig.2f at the arrow positions). TElocal is also reduced by up to 40% in the basal ganglia (Fig.2d, slice 12, see arrows). Ventrally, close to the ear canals and dorsally in frontal lobes the TElocal increases strongly leading to type 2 signal loss (Fig.2d and f, white regions).
Results in the four tumor patients with the strongest gradients in B0 are presented in Fig. 3. In all cases large regions with strongly reduced TElocal (by about 50%) but relatively high tSNR were found in the proximity of pathologies (see arrows). For patient 1 and 2 this occurred close to primary motor cortex, for patient 3 near Broca’s area and patient 4 near auditory cortex. In BS metric 3 the tSNR values are readjusted accordingly (Fig.2f).
Discussion and conclusion
We have shown that tSNR maps may indicate sufficient BS to detect activation but that BS may be dramatically reduced by changes to TE
local. This effect was prominent at 7T in patients with large oedema, postoperative cavities close to vital brain functions (motor, speech, auditory). Neglecting changes to TE
local can lead to false negative results in clinical fMRI being hard to detect, because underlying EPI and tSNR indicate that there is adequate signal. We therefore advise that both, TE
local and tSNR, be considered in the assessment of bold sensitivity.
Acknowledgements
This study was funded by a DOC fellowship of the Austrian Academy of Science. Additional support was provided by the Austrian Science Fund (FWF KLI 264).References
1. Robinson S, Windischberger C, Rauscher A, Moser E. Optimized 3 T EPI of the amygdalae. NeuroImage 2004;22:203–210. doi: 10.1016/j.neuroimage.2003.12.048.
2. Weiskopf N, Hutton C, Josephs O, Deichmann R. Optimal EPI parameters for reduction of susceptibility-induced BOLD sensitivity losses: a whole-brain analysis at 3 T and 1.5 T. NeuroImage 2006;33:493–504. doi: 10.1016/j.neuroimage.2006.07.029.
3. Triantafyllou C, Hoge RD, Krueger G, Wiggins CJ, Potthast A, Wiggins GC, Wald LL. Comparison of physiological noise at 1.5 T, 3 T and 7 T and optimization of fMRI acquisition parameters. NeuroImage 2005;26:243–250. doi: 10.1016/j.neuroimage.2005.01.007.
4. Lipschutz B, Friston KJ, Ashburner J, Turner R, Price CJ. Assessing Study-Specific Regional Variations in fMRI Signal. NeuroImage 2001;13:392–398. doi: 10.1006/nimg.2000.0687.
5. Deichmann R, Josephs O, Hutton C, Corfield DR, Turner R. Compensation of Susceptibility-Induced BOLD Sensitivity Losses in Echo-Planar fMRI Imaging. NeuroImage 2002;15:120–135. doi: 10.1006/nimg.2001.0985.
6. Parrish TB, Gitelman DR, LaBar KS, Mesulam M-M. Impact of signal-to-noise on functional MRI. Magn. Reson. Med. 2000;44:925–932. doi: 10.1002/1522-2594(200012)44:6<925::AID-MRM14>3.0.CO;2-M.
7. Posse S. Direct imaging of magnetic field gradients by group spin-echoselection. Magn. Reson. Med. 1992;25:12–29. doi: 10.1002/mrm.1910250103.
8. Bernstein MA, Grgic M, Brosnan TJ, Pelc NJ. Reconstructions of phase contrast, phased array multicoil data. Magn. Reson. Med. 1994;32:330–334. doi: 10.1002/mrm.1910320308.
9. Chen N, Oshio K, Panych LP. Application of k-space energy spectrum analysis to susceptibility field mapping and distortion correction in gradient-echo EPI. NeuroImage 2006;31:609–622. doi: 10.1016/j.neuroimage.2005.12.022.