i-fMRI: A virtual whiteboard discussion on multi-echo, simultaneous multi-slice EPI

Disclaimer: This isn't an April Fool!

I'd like to use the collective wisdom of the Internet to discuss the pros and cons of a general approach to simultaneous multislice (SMS) EPI that I've been thinking about recently, before anyone wastes time doing any actual programming or data acquisition.

Multi-echo EPI for de-noising fMRI data

There has been quite a lot of interest in using multi-echo EPI to characterize and de-noise time series data, e.g.

http://www.ncbi.nlm.nih.gov/pubmed/22209809

These methods rest on one critical aspect: they use in-plane parallel imaging (GRAPPA or SENSE, usually depending on the scanner vendor) to render the per slice acquisition time reasonable. For example, with R=2 acceleration it's possible to get three echo planar images per slice at TEs of around 15, 40 and 60 ms. The multiple echoes can then be used to characterize BOLD from non-BOLD signal variations, etc.

The immediate problem with this scheme is that the per slice acquisition time is still a lot longer than for normal EPI, meaning less brain coverage. The suggestion has been to use MB/SMS to regain speed in the slice dimension. This results in the combination of MB/SMS in the slice dimension and GRAPPA/SENSE in-plane, thereby complicating the reconstruction, possibly (probably) amplifying artifacts, enhancing motion sensitivity, etc. If we could eliminate the in-plane parallel imaging and do all the acceleration through MB/SMS then that would possibly reduce some of the artifact amplification, might simplify (slightly) the necessary reference data, etc.

A different approach? 

Physics for understanding fMRI artifacts: Part Thirteen

A tour through a real EPI pulse sequence

In some posts I've got planned it will be important for you to know something about all of the different functional modules that are included in a real EPI pulse sequence. So far in this PFUFA series I've used schematics of the particular segment of the sequence that I was writing about, e.g. the echo train that covers 2D k-space for single-shot EPI. Except that there comes a time when you need to know about the sequence in its entirety, as it is implemented on a scanner. Why? Because there are various events that I've given short shrift - fat saturation and N/2 ghost correction, for instance - that have significant temporal overheads in the sequence, and these additional delays obviously affect how quickly one can scan a brain.

So, without further ado, here is a pulse sequence for fat-suppressed, single-shot gradient echo EPI, as used for fMRI:

(Click to enlarge.)

Okay, so it's not the entire pulse sequence. The readout gradient echo train in this diagram has been curtailed after just nine of 64 total gradient echoes that will be acquired, for EPI with a matrix of 64x64 pixels. The omitted 55 echoes are simply clones of the nine echoes that you can see. (Note that there are no additional gradient episodes at the end of this particular EPI sequence; all the crusher gradients occur at the start of the sequence and these are visible in the above figure. More on crusher gradients below.) I should also point out that this is the timing diagram for acquisition of a single 64x64 matrix EPI slice. The pulse sequence as shown would be repeated n times for n slices within each TR. (See Note 1.)


Interpreting what you see

Let's first determine what information is being displayed on the figure above. There are five axes, all handily labeled on the far right-hand side of the figure. The top axis is the RF transmit channel; we've got two RF pulses in this sequence. The second axis down is the receiver, or analog-to-digital converter (ADC) channel. The scanner is receiving signals only when there's a rectangle specified on the second axis. Finally, the bottom three axes represent the pulsed field gradients, in the order X, Y, Z.

Just for fun, let's quickly determine what the scanner is doing in the logical frame of reference, before we delve into the nitty-gritty. The slice selection gradient will occur in concert with an RF excitation pulse, and we have two RF pulses to choose from. Slice selection can't be the first RF pulse because that pulse occurs without any concomitant gradients. Thus, the slice excitation pulse must be the second one and we can deduce that slice selection is along the Z axis, which is the magnet bore axis. We're doing axial slices.