Education, tips and tricks to help you conduct better fMRI experiments.
Sure, you can try to fix it during data processing, but you're usually better off fixing the acquisition!

Tuesday, December 16, 2014

Updated checklist for fMRI acquisition methods reporting in the literature


This post updates the checklist that was presented back in January, 2013. The updated checklist is denoted Version 1.2. The main update is to include reporting for simultaneous multi-slice (SMS) (a.k.a. multi-band, MB) EPI.

Explanatory notes for parameter names appear in the lower portion of the post. Note that the present checklist was devised by considering typical fMRI experiments conducted on 1.5 T and 3 T scanners but the list should work reasonably well for 7 T sequences, too.

Please keep the comments and feedback coming. This is an ongoing, iterative process.




Release notes for Version 1.2

All changes from Version 1.1 have been highlighted in yellow, both on the list PDF and on the explanatory notes (below).

1. The "Spatial Encoding" parameter categories have been renamed "In-Plane Spatial Encoding" to better differentiate in-plane acceleration (e.g. GRAPPA) from slice dimension acceleration (SMS/MB).

2. When using slice dimension acceleration (i.e. SMS/MB), certain parameters that are listed as Supplemental for other EPI variants should be considered Essential. Specifically, it is suggested to report:
Matrix coil mode
Coil combination method

All the In-Plane Spatial Encoding parameters in the Supplemental category should be considered because there is a tendency to use SMS/MB to attain high spatial resolution, requiring long readout echo trains that can have higher distortion than found in typical EPI scans.

The In-plane reconstructed matrix parameter should be reported whenever partial Fourier sampling is used, as it often is for SMS/MB EPI.

All the RF & Contrast parameters in the Supplemental category should be reported because the shape, duration and amplitude of the excitation RF pulse are all integral components of the acceleration method.

The Shim routine should be reported if a non-standard shim is performed before SMS/MB EPI.

3. Pre-scan normalization has been added to the Supplemental section of RF & Contrast parameters. Large array coils produce strong receive field heterogeneity and the use of pre-scan normalization may improve the performance of post hoc motion correction.



Explanatory notes


Essential - Scanner

Magnetic field strength: In lieu of magnetic field strength the scanner operating frequency (in MHz) might be considered acceptable. I'm assuming we're all doing 1-H fMRI. If not, chances are the methods section's going to be very detailed anyway.


Essential - Hardware Options

Rx coil type: For standard coils provided by the scanner vendor a simple description consisting of the number of independent elements or channels should suffice, e.g. a 12-channel phased array coil, a 16-leg birdcage coil. Custom or third-party coils might warrant more detailed information, including the manufacturer. Most head-sized coils are Rx-only these days, but specifying Rx-only doesn't hurt if there could be any ambiguity, e.g. a birdcage coil could quite easily be Tx/Rx.


Essential - In-Plane Spatial Encoding

Pulse sequence type: A generic name/descriptor is preferred, e.g. single-shot EPI, or spiral in/out.

Number of shots (if > 1): Multi-shot EPI isn't a common pulse sequence; single-shot EPI is by far the most common variant, even if acquired using parallel imaging acceleration. I include this option to reinforce the importance of reporting the spatial encoding accurately.

PE acceleration factor (if > 1): This is usually called the acceleration factor, R in the literature. Vendors use their own notation, e.g. iPAT factor, SENSE factor, etc.

PE acceleration type (if > 1): It seems that most vendors use the generic, or published, names for parallel imaging methods such as SENSE, mSENSE and GRAPPA. I would think that trade names would also be acceptable provided that the actual (published) method can be deciphered from the scanner vendor and scanner type fields. But generic names/acroynms are to be preferred.

PE partial Fourier scheme (if used): Convention suggests listing the acquired portion/fraction of k-space rather than the omitted fraction. Any fraction that makes sense could be used, e.g. 6/8 or 48/64 are clearly equivalent.


Essential - Spatial Parameters

In-plane matrix: This should be the acquired matrix. If partial Fourier is used then I would suggest reporting the corresponding full k-space matrix and giving the partial Fourier scheme as listed above. But I wouldn't object to you reporting that you acquired a 64x48 partial Fourier matrix and later reported the reconstructed matrix size as 64x64. So long as everything is consistent it's all interpretable by a reader. (But see also the In-plane reconstructed matrix parameter in the Supplemental section.)

In-plane inline filtering (if any): Non-experts may be unaware that filtering might be applied to their "raw" images before they come off the scanner. It's imperative to check and report whether any spatial smoothing was applied on the scanner as well as during any "pre-processing" steps subsequent to porting the time series data offline.

Slice thickness and Inter-slice gap: For now I would use the numbers reported by the scanner, even though there may be some small variation across scanner vendors and across pulse sequences. For example, some vendors may use full width at half slice height while others may use a definition of width at or near the base, there may be pulse sequence options for RF excitation pulse shape and duration, etc. I see these details as secondary to the essential reporting improvements we're aiming for.

Slice acquisition order: Interleaved or contiguous would be sufficient, although explicit descending (e.g. head-to-foot) or ascending (e.g. foot-to-head) options for contiguous slices would be acceptable, too. Presumably, the subsequent use of slice timing correction will be reported under the post-processing steps (where most fMRIers call these "pre-processing" because they are applied before the statistical analysis).


Essential - Timing Parameters

TR: For single-shot EPI with no sparse sampling delay the TR also happens to be the acquisition time per volume of data. But if sparse sampling or multiple shot acquisition is being used then the TR should be clearly reported relative to these options. The conventional definition of TR is the time between successive RF excitations of the same slice. Thus, by this definition the reported TR would include any sparse sampling delay, but it would specify the time of each separate shot in a multi-shot acquisition and the per volume acquisition time would become TR x nshots.

No. of volumes in time series: Dummy scans (not saved to disk) should be reported separately. Likewise, the use or rejection of the first n volumes for task-related reasons, e.g. to allow a subject to acclimatize to the scanner sounds, should also be reported separately in the post-processing segment of the experiment. In this field we are only interested in the total quantity of data available to the experimenter.

No. of averages/volume (if > 1): I don't think I've ever seen anyone do anything but one average per TR for single-shot EPI/spiral (unless they've screwed something up) and I can't think of a reason why someone would want to do it for fMRI. But, if it happens for some reason then it's really, really important to specify it.


Essential - RF & Contrast

Fat suppression scheme: It's sufficient to state that fat saturation or fat suppression was used, for example. Further details aren't required unless the scheme was non-standard, e.g. a custom spatial-spectral excitation scheme.


Essential - Slice Acceleration

SMS/MB acceleration factor:  The number of slices excited simultaneously, equivalent to the spatial acceleration in the slice dimension. The total number of slices (Number of slices parameter) must be divisible by the SMS/MB acceleration factor. Note that the slice dimension acceleration (this parameter) is independent of any in-plane acceleration (e.g. SENSE, GRAPPA) and should be reported independently. For example, GRAPPA with R=2 should be reported as suggested in the In-Plane Spatial Encoding section, while slice acceleration should be reported here. It is not acceptable to combine the parameter reporting as one overall acceleration because the practical consequences of each are very different.

Sequence/recon name &/or version:  Unlike many single-shot EPI sequences which are product, most SMS/MB sequences are being actively developed and change frequently. It is imperative to report version numbers for both the acquisition and, if they aren’t coupled, the reconstruction. If reconstruction is performed offline, report the method used, and any version number, in full.


Essential - Customization

Sparse sampling delay (if used): Sometimes called "Delay in TR" on the scanner interface. Used most often for auditory stimulus or auditory response fMRI.

Prospective motion correction scheme (if used): PACE is one commercial option. These schemes fundamentally change the nature of the time series data that is available for subsequent processing and should be distinguished from retrospective (post-processing) corrections, e.g. affine registration such as MCFLIRT in FSL. It is also critical to know the difference between motion correction options on your scanner. On a Siemens Trio running VB15 or VB17, for instance, selecting the MoCo option enables PACE and a post hoc motion correction algorithm if you are using the sequence, ep2d_pace whereas only the post hoc motion correction algorithm - no PACE - is applied if you are using the ep2d_bold sequence. There's more detailed information on these options in my user training guide/FAQ.

Cardiac gating (if used): This isn't a common procedure for fMRI, and recording of cardiac information, e.g. using a pulse oximeter, is radically different from controlling the scanner acquisition via the subject's physiology. The recording of physiological information doesn't usually alter the MRI data acquisition, whereas gating does. Reporting of physio information is tangential to the reporting structure here, but if you are recording (and using) physio data then presumably you will report it accordingly somewhere in the experiment description.


Supplemental - Hardware Options

Gradient set type: It should be possible to infer the gradient coil from the scanner model. If not, e.g. because of a custom upgrade or use of a gradient insert set, then the specifications of the actual gradient coil should be reported independently.

Tx coil type (if non-standard): It should be possible to infer the Tx coil from the scanner model. If not, e.g. because of a custom upgrade or use of a combined Tx/Rx coil, then the new Tx coil should be reported independently. I would also advocate including the drive system used if the coil is used in anything but the typical quadrature mode.

Matrix coil mode (if used): There are typically default modes set on the scanner when one is using un-accelerated or accelerated (e.g. GRAPPA, SENSE) imaging.  If a non-standard coil element combination is used, e.g. acquisition of individual coil elements followed by an offline reconstruction using custom software, then that should be stated.

Coil combination method: Almost all fMRI studies using phased-array coils use root-sum-of-squares (rSOS) combination, but other methods exist. The image reconstruction is changed by the coil combination method (as for the matrix coil mode above), so anything non-standard should be reported.


Supplemental - In-Plane Spatial Encoding

PE direction: If you've shown any examples of EPI in your paper then the PE direction can usually be determined from the image. If N/2 ghosts or distortions aren't obvious, however, then it's rather important that the phase encode direction is stated, in concert with the readout echo spacing, so that a reader can infer your spatial accuracy.

Phase oversampling (if used): There's no reason to use phase oversampling for EPI - you're wasting acquisition time - but if you are using it for some reason then it should be reported consistent with the acquired matrix, acquired FOV, echo spacing and associated parameters.

Read gradient bandwith: Not an intrinsically useful parameter on its own, it does have value if reported in conjunction with the echo spacing. (An alternative to this single parameter would be the read gradient strength (in mT/m) and the digitizer bandwidth (in kHz).)

Readout echo spacing: Rarely reported but really useful! This number, in conjunction with the FOV and acquired matrix size, allows a reader to estimate the likely distortion in the phase encode direction.

Pulse sequence name: Could be invaluable for someone wanting to replicate a study. There may be multiple similar pulse sequences available, all capable of attaining the specifications given, but it is entirely feasible that only one of the sequences has a particular quirk in it!

k-space scheme: Readers will assume linear (monotonic) k-space steps in the PE direction unless indicated to the contrary. Centric ordering or other atypical schemes should be indicated, especially in concert with multiple shots if the Number of shots parameter is greater than one.

Read gradient strength: Could be useful in conjunction with information about the ramp sampling percentage and echo spacing time, otherwise probably of limited value to most readers.

(Ramp sampling percentage:) Ramp sampling can increase the N/2 ghost level considerably if there is appreciable gradient and digitization (data readout) mismatch. But determining the percentage of readout data points that are acquired on the flat vs. the ramp portions of each readout gradient episode can be involved. And for routine studies, as opposed to method development studies, there's probably not a whole lot of value here. Maybe remove it?

(Ghost correction method:) N/2 ghost correction usually happens invisibly to the user, but there are some options becoming available, especially useful for large array coils (e.g. 32-channel coils) where there may be local instabilities with some ghost correction methods. If known, and if non-standard, then it would be nice to report. But perhaps more overkill for fMRI methods?

PE partial Fourier reconstruction method: If the scanner offers more than one reconstruction option then the chosen option should be reported.


Supplemental - Spatial Parameters

In-plane resolution: This field is redundant if the reconstructed matrix (In-plane matrix parameter) and FOV are reported, but I for one wouldn't object to seeing the nominal in-plane pixel size given anyway. It may make the paper a faster read. Probably not worth arguing about. (Cue extensive debate...!)

In-plane reconstructed matrix: This is for reporting of zero filling (beyond the default zero filling that may have been done for a partial Fourier acquisition) to a larger matrix than acquired, prior to 2D FT. There may be modeling issues associated with the number of voxels in the image, not least of which is the size of the data set to be manipulated! It could save someone a lot of angst if she knows what you did to the data prior to uploading it to the public database.


Supplemental - Timing Parameters

No. of dummy scans: This is the number of dummy scans used to establish the T1 steady state. Many fMRI experiments also use acquired volumes subsequent to the (unsaved) dummy scans for neuroscientific reasons, e.g. to allow definition of a BOLD baseline or adjustment to the scanner noise. These two parameters should be reported separately.


Supplemental - RF & Contrast

Excitation RF pulse shape and Excitation RF pulse duration:  Not critical for standard pulse sequences on commercial scanners, but if atypical values are set in order to achieve very thin slices, for example, then reporting these parameters may be benenficial. These two parameters are essential, however, when reporting SMS/MB EPI.

Pre-scan normalization: Typically achieved with a low-resolution gradient echo scan, but the actual process may be invisible to the user. Instead, activating pre-scan normalization may simply require selection of the appropriate option. (This is the case on Siemens scanners.)


Supplemental - Slice Acceleration

SMS/MB reconstruction type:  Report details of the reconstruction type if possible, e.g. the kernel size for a GRAPPA-style reconstruction.

FOV shift:  The sequence may use blipped controlled aliasing along the slice direction to assist in the separation of the simultaneous slices. The blips are set to produce a fixed partial FOV shift, typically FOV/2, FOV/3 or FOV/4. The FOV shift may be user-controlled or set to a fixed default.

SIR/SER factor:  If simultaneous image (or echo) refocusing is used in addition to SMS/MB (as in Feinberg et al., PLoS One, 2010), report the SIR/SER acceleration factor as well as the SMS/MB acceleration. 


Supplemental - Customization

Image reconstruction type: Unless specified to the contrary, most readers will assume that magnitude images were taken from the 2D Fourier transform that yielded each individual EPI. If you use complex data - magnitude and phase - then that option should be specified along with the particular processing pipeline used to accommodate the atypical data type.

Shim routine: If manual shimming or an advanced phase map shimming routine is used, especially to improve the magnetic field over restricted brain volumes, then this information should be reported.

Receiver gain: Most scanners use some form of autogain to ensure that the dynamic range at the receiver is acceptable. If manual control over receiver gain is an option and is used then it should be reported because a mis-set gain could lead to artifacts that aren't typically seen in EPI, and a reader could subsequently attribute certain image features to other artifact sources.


_________________



Abbreviations:

FOV - Field-of-view
MB - Multi-band
N/2 - Half-FOV (for Nyquist ghosts)
PE - Phase encode
Rx - Radiofrequency (RF) receiver
SER - Simultaneous echo refocusing
SIR - Simultaeous image refocusing
SMS - Simultaneous multi-slice
Tx - Radiofrequency (RF) transmitter
TE - Echo time
TR - Repetition time


No comments:

Post a Comment