Improved detection of fMRI activity in ventromedial prefrontal cortex using multi-echo EPI
Brice Fernandez1, Laura Leuchs2, Phillip G. Sämann 2, Michael Czisch2, and Victor I. Spoormaker2

1Applications and Workflow, GE Healthcare, Munich, Germany, 2Neuroimaging Unit, Max Planck Institute of Psychiatry, Munich, Germany

Synopsis

Standard fMRI suffers from signal loss in the ventromedial and orbital prefrontal cortex, a region of special interest in affective neuroimaging. Multi-echo EPI (MEPI) is known to have several advantages over EPI. In this work, we test if MEPI is able to reach better performance in detecting task-induced activation in the ventromedial prefrontal cortex (vmPFC) during fear conditioning, known for eliciting activity in this area. We demonstrate that MEPI (by means of the weighted sum combination approach) outperforms standard EPI in vmPFC, which is highly relevant for affective neuroscience and psychiatry given its critical role in emotion regulation.

Introduction

A limitation of fMRI in the investigation of affective symptomatology is signal dropout in EPI, especially pronounced in ventromedial prefrontal cortex (vmPFC), which is a region of special interest in affective neuroscience and psychiatry1. Standard fMRI at 3 Tesla leads to significant signal and BOLD sensitivity loss in vmPFC and other orbitofrontal and temporal regions affected by signal dropout2. Multi-echo EPI (MEPI) is known to have several advantages3-9 over EPI. To date, studies have indicated3-8 tends better activation detection, also in vmPFC. However, to our very best knowledge, no study has reported a statistical comparison of MEPI with the standard fMRI using an appropriate task that specifically elicits activity in vmPFC.

Purpose

To test if MEPI, with the weighted sum combination approach5, outperforms standard EPI in detecting task induced activation in vmPFC using a task that robustly activates this area.

Methods

Thirty-six young and healthy subjects (age 25.0±3.2 years, 17 female) were included, four data sets had to be excluded (e.g. incomplete data).

In the fear-conditioning task, three visual stimuli (geometric shapes) were presented 10 times each with 12-14s inter-stimulus-intervals. Two of the stimuli were paired with mild electric shocks at stimulus offset, administered to the wrist at an individually titrated level that was uncomfortable but not painful. The two conditioned stimuli were paired with shocks in 100% (100% reinforcement) and 60% (60% reinforcement) of all stimulus presentations and can be summarized as CS+. The third stimulus was never paired with the shock (safety stimulus – CS-). The comparison CS->CS+ generally shows the vmPFC10. The reverse contrast shows an increase of activity in dorsal anterior cingulate (dACC) and bilateral insula10.

The data were acquired on a GE Discovery MR750 3T scanner with a 32-channel head coil. The fMRI data were acquired using a MEPI pulse sequence with 3 echoes (E1, E2, E3) in an orientation aligned with the AC-PC. Acquisition parameters were: 64x64 matrix, 3.4mm isotropic resolution, 36 slices, acceleration-factor 2, TR=2.56s, TE1/TE2/TE3=12/29/46ms. Additionally, a 3D high-resolution T1w image was acquired.

The fMRI analysis was performed with SPM8. Data were first corrected for slice timing differences and realigned. At this stage, the weighted sum combination5 was calculated resulting in a weighted echo Ec. Then, the MEPI data were co-registered to the 3D-T1, the 3D-T1 was normalized to the MNI space and the spatial transformation applied to the fMRI data. To be able to statistically compare the methods, a linear detrending followed by a z-score transformation was applied to all voxel time courses of MEPI (E1, E2, E3 and Ec) before spatial smoothing (Gaussian FWHM 6x6x6 mm3).

First level models comprised hemodynamic response function convolved stimulus regressors, motion coefficients as nuisance regressors were estimated, and the contrast CS ->CS+ was computed. At the second level a full factorial analysis of variance was used to compare the different methods (i.e. E1, E2, E3, Ec) with a repeated measures factor of four levels (method). The primary contrast of interest was the comparison of E2 (standard fMRI) with Ec. To control for multiple tests, familywise error (FWE) correction were used.

Results

All methods showed the expected task activation but at different levels of statistical significance (see Figure 1). Figure 2 shows the main effect of method to highlight the regions in which the methods differ, together with the contrasts estimates for three exemplary differentially activated regions. Figure 3 depicts the Ec>E2 (standard EPI) comparison where both the vmPFC and the bilateral insula demonstrate significant differences. Figure 4 illustrates the contribution (the weights) of each echo in the brain regions of interest. Figure 5 shows the extension of the clusters in Ec compared to the others while showing the regions affected by the dropout.

Discussion

Our results show that MEPI outperforms standard EPI in detecting task-induced vmPFC activation, even with a basic weighted sum. Further, results show the contribution of each echo in different brain regions (Figure 4, 5). It is worth noting that MEPI does not only improve the contrast in brain regions affected by signal dropout where E1 contributes significantly, but also in homogeneous brain regions where E3 contributes significantly because it is closer to the optimal TE than E2. This is in accordance with previous studies2-10 and a strong advantage of MEPI. Recent methodologies like ME-ICA7-8, multiband and coil development are expected to take these results further.

Conclusion

We have demonstrated that, in a fear-conditioning task, MEPI allow an overall better detection of fMRI activity in vmPFC and other involved brain regions. Future work will include advanced analyses7-8, other tasks and resting state fMRI.

Acknowledgements

No acknowledgement found.

References

1. Price and Drevets, Neuropsychopharmacology, 2010;35(1):192-216.

2. Deichmann R, Josephs O, Hutton C, et al. Compensation of susceptibility-induced BOLD sensitivity losses in echo-planar fMRI imaging. Neuroimage. 2002;15:120-135.

3. Posse S. Multi-echo acquisition. Neuroimage. 2012;62:665-671.

4. Posse S, Wiese S, Gembris D, et al. Enhancement of BOLD-contrast sensitivity by single-shot multi-echo functional MR imaging. Magn Reson Med. 1999;42:87-97

5. Poser, B. A.; Versluis, M. J.; Hoogduin, et al. BOLD contrast sensitivity enhancement and artifact reduction with multiecho EPI: parallel-acquired inhomogeneity-desensitized fMRI. Magn Reson Med. 2006;55:1227-1235.

6. Bhavsar S, Zvyagintsev M, Mathiak K. BOLD sensitivity and SNR characteristics of parallel imaging-accelerated single-shot multi-echo EPI for fMRI. Neuroimage. 2013;84C:65-75

7. Kundu P, Inati S, Evans J, et al. Differentiating BOLD and non-BOLD signals in fMRI time series using multi-echo EPI. Neuroimage. 2012;60:1759-1770.

8. Kundu P, Brenowitz N, Voon V, et al. Integrated strategy for improving functional connectivity mapping using multiecho fMRI. Proc Natl Acad Sci USA. 2013;110:16187-16192

9. Fernandez B, Czish M. Simulation of BOLD Sensitivity of Single-Shot Multi-Echo EPI versus Sample-Induced Susceptibility Gradients. ISMRM, 2014;2982

10. Fullana M, Harrison B, Soriano-Mas C, et al. Neural signatures of human fear conditioning: an updated and extended meta-analysis of fMRI studies. Mol Psychiatry. 2015; 1476-5578.

Figures

Figure 1: The contrast of the safety versus fear conditioned stimuli CS->CS+ (hot colors) and its reverse (cold colors) in MNI space for the four methods E1, E2, E3 and Ec. T-contrast p(uncorrected)<10-4 (T>4.22) and cluster extent>160. Note that all cluster shown are p(FWE clusterwise)<0.05.

Figure 2: Main effect of the factor method in MNI space. (a) F-contrast p(FWE voxelwise)<0.05, cluster extent>20 and the contrast estimates of the peak voxel in vmPFC (b), dACC (c) and right insula (d).

Figure 3: Comparison Ec >E2 for the contrast CS->CS+ (hot color) and its reverse (cold color) in MNI space. T-contrast with p(FWE clusterwise)<0.05, cluster extent>100. Note that both effect directions indicate improved detection of activation/deactivation by Ec compared with E2.

Figure 4: For illustrative purposes, color-coded representation in MNI space of the contribution (the weights) of E1 (line 1), E2 (line 2) and E3 (line 3) used to create Ec. Contour of the clusters shown on figure 2 are also displayed.

Figure 5: For illustration purpose, activation shown on figure 1 overlaid on their respective mean (across subjects) normalized EPI images (MNI coordinate [12 0 -18]). Note the extension of the clusters in Ec compared to E1, E2, E3 and in regards of the brain regions affected by signal dropout.



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