Optimizing SMS-BOLD image reconstruction for resting state analysis and reconstruction time
Ross W. Mair1,2, R. Matthew Hutchison1,3, Stephanie McMains1, and Steven Cauley2

1Center for Brain Science, Harvard University, Cambridge, MA, United States, 2A.A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States, 3Department of Psychology, Harvard University, Cambridge, MA, United States

Synopsis

The computational processes required for slice-unaliasing in SMS-BOLD scans are taxing on the scanner reconstruction computers, so data is sometimes observed far from real-time, and reconstruction may lag up to tens of minutes behind the acquisition. A channel compression algorithm has been proposed to counter computational demands of these reconstruction processes. We studied functional networks derived from resting-state scans as a function of slice acceleration, slice-GRAPPA kernel size and channel compression to find an optimal solution for an existing, conventional 3.0 T scanner. Slice-GRAPPA kernel size played little effect in functional network definition and two-fold channel compression was beneficial to reconstruction time without impacting functional network data quality.

Purpose

Slice-accelerated EPI using multiband RF pulses that allow simultaneous multi-slice (SMS) acquisition of BOLD contrast images1,2 can significantly enhance the temporal and spatial resolution of fMRI. This advance enables whole-brain fMRI at sub-second volume-TRs,3,4 which is advantageous for resting-state fMRI, which attempts to capture sub-second frequency fluctuations.5,6 However, the computational processes required for slice-unaliasing can be taxing for the reconstruction computer, so data can be observed far from real-time, and reconstruction may lag up to tens of minutes behind the acquisition. A channel compression algorithm has been proposed to counter these computational demands.7 Here, functional networks derived from resting-state scans are studied as a function of slice acceleration, slice-GRAPPA kernel size,2 and channel compression to find an optimal solution in terms of data quality and reconstruction time for a conventional 3.0 T scanner.

Methods

All measurements were performed using a Siemens Tim Trio 3.0 T MRI scanner equipped with the fastest Step4 image reconstruction computer, and using the vendor-supplied 32-channel head coil. Slice-accelerated BOLD images were acquired using Siemens WIP770A sequence with modified reconstruction code that improved computational speed and enabled channel compression. All resting-state scans employed ~6 min of data acquisition while subjects rested with eyes open. A 3mm isotropic, unaccelerated scan was acquired with the conventional BOLD sequence as a reference point. Slice-accelerated scans were acquired at 2mm nominal resolution (6/8 partial-fourier), with SMS=1, 3, 4, 6 or 8, enabling TR=5.90, 2.00, 1.44, 0.93 and 0.76 sec respectively. TE=30 ms, BW~1600 Hz/px, echo spacing~0.71 ms for all 2mm scans. The raw data was reconstructed at the scanner multiple times with varying slice-GRAPPA kernel sizes (3×3, 5×5, 7×7) and channel compression factors (CCR=1, 2, 4). Notation throughout indicates kernel size and CCR as: kX×kY×CCR. Seed-based functional connectivity analysis was performed using in-house scripts.8 Slice-time correction was performed for the 3mm and MB1 scans, but omitted for all slice-accelerated scans.

Results

Fig. 1 shows an example for the Default Network (posterior cingulate seed8) where SMS=4 shows a larger number of voxels detected than the SMS=8 case, using the same reconstruction parameters. For a given SMS factor, larger kernel sizes with CCR=1 or 2 yielded improved tSNR over the default 3×3×1 reconstruction, however this generally did not improve network definition. Fig. 2 shows example maps of the Default Network for SMS=8, with different levels of channel compression. CCR=2 has minimal effect on network definition, but CCR=4 resulted in lower tSNR for all kernel sizes, and significant loss of network definition as well as introducing significant noise and false positives. Similar results are seen in the motor and visual networks. Fig. 3 highlights differences between pairs of reconstructions (from the same acquisition). z scores for major network nodes are slightly lower when CCR>1, usually ~0.05 – these have been masked out. Fig. 4 shows voxel counts for Default and motor networks as a function of SMS and reconstruction parameters. While SMS=8 (and maximum number of time points) aids motor network definition, the default network is optimally shown with SMS=4.

Discussion

Despite the increase in the number of time-points acquired, increasing the SMS factor to 6 or 8 was not always beneficial for resting-state analysis. While slice-accelerated scans always resulted in better network definition over unaccelerated scans, for some networks, SMS=4 provided the maximum number of voxels above threshold. Motor and visual networks benefited the most from SMS=8 acceleration. For a given SMS-acquired dataset, 5×5 kernels had been commonly used for SMS>6 accelerations,2,3 but this is unnecessary for fcMRI analysis. Although the tSNR is increased by the larger kernel size (after normalizing out other effects), this had minimal effect on network definition, while reconstruction times were heavily impacted. CCR=2 had minimal effect on data from SMS=4 and SMS=8 scans reconstructed with any kernel size, however CCR=4 introduced unacceptable levels of false activations and loss of network definition. Without channel compression, increasing the kernel size lengthened the reconstruction time by ~40 %. CCR=2 results in a 10% time reduction for a 3×3 kernel, but ~40% time reduction for 5×5. One acceptable compromise for maximum tSNR without impacting reconstruction time may be to use 5×5×2, although similar fcMRI maps will likely be achievable with 3×3×1. These conclusions relate specifically to the use of the Siemens 32-channel head coil on a Tim Trio scanner – higher kernel sizes and compression factors may be optimal with higher-array head coils or other scanner platforms, especially those equipped with GPU-assisted reconstruction computers. Future work will consider the effect of physiological noise effects over the different TR values on the results above.

Acknowledgements

Harvard Center for Brain Science; NIH Grant P41-RR14075. Himanshu Bhat and Thomas Benner of Siemens for sequence modifications.

References

1. Moeller S, et al, Multiband multislice GE-EPI at 7 Tesla with 16-fold acceleration using Partial Parallel Imaging with application to high spatial and temporal whole-brain FMRI. Magn. Reson. Med., 2010;63:1144–1153.

2. Setsompop K, et al, Blipped-controlled aliasing in parallel imaging for simultaneous multislice echo planar imaging with reduced g-factor penalty. Magn. Reson. Med., 2012;67:1210-1224.

3. Feinberg DA, Moeller S, et al, Multiplexed echo planar imaging for sub-second whole brain FMRI and fast diffusion imaging. PLoS One, 2010;5:e15710.

4. Feinberg DA, Setsompop K, Ultra-fast MRI of the human brain with simultaneous multi-slice imaging. J. Magn. Reson., 2013;229:90–100.

5. Smith SM, Miller KL, Moeller S, et al, Temporally-independent functional modes of spontaneous brain activity. Proc. Natl. Acad. Sci., 2012;109:3131–3136.

6. Smith SM, Beckmann CF, Andersson J, Auerbach EJ, Bijsterbosch J, et al, Resting-state fMRI in the Human Connectome Project. NeuroImage 2013;80:144-168.

7. Cauley SF, et al, Geometric-decomposition Coil Compression for Real-time Simultaneous MultiSlice EPI reconstruction at high MultiBand factors. Proc ISMRM 2014;22:4404.

8. van Dijk KR, et al, Intrinsic Functional Connectivity as a tool for Human Connectomics: Theory, Properties, and Optimization. J. Neurophysiol., 2010:103: 297–321.

Figures

Maps of the Default Network (posterior cingulate seed), for SMS 4 or 8, with different Slice-GRAPPA kernels. Variation in network definition is minimal for 3×3 and 5×5 kernels, at both SMS 4 or 8, except when CCR=2 is used, in which case 3×3 kernels show slightly better network extent.

Maps of the Default Network from a posterior cingulate seed, for SMS8, with different levels of channel compression. CCR=2 has minimal effect on network definition, but CCR=4 degrades the network significantly, while increasing noise. Similar results are seen in other networks

Differences in network maps between pairs of datasets with different reconstruction parameters. Minimal difference is observed between 3×3×1 and 3×3×2, or 5×5×1 and 3×3×1. However, noticeable differences are apparent in 3×3×4 vs 3×3×1. Pixel count analysis showed CCR=4 resulted in lower network node definition but higher false positive voxels.

Number of voxels with z>0.4 for the Default and Motor networks, as a % of voxels in the 3mm unaccelerated scan. Increased slice acceleration generally increases network definition, although SMS=4 gave maximal extent for the Default network. Despite SNR increases, larger kernels don’t improve network extent for either network.



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