Skip to main content

METHODS article

Front. Neurosci., 19 September 2014
Sec. Brain Imaging Methods
This article is part of the Research Topic Functional brain mapping of epilepsy networks: methods and applications View all 24 articles

De-noising with a SOCK can improve the performance of event-related ICA

  • 1The Florey Institute of Neuroscience and Mental Health, Austin Hospital, The University of Melbourne, Melbourne, VIC, Australia
  • 2Department of Medicine, The University of Melbourne, Melbourne, VIC, Australia
  • 3Department of Radiology, The University of Melbourne, Melbourne, VIC, Australia

Event-related ICA (eICA) is a partially data-driven analysis method for event-related fMRI that is particularly suited to analysis of simultaneous EEG-fMRI of patients with epilepsy. EEG-fMRI studies in epileptic patients are typically analyzed using the general linear model (GLM), often with assumption that the onset and offset of neuronal activity match EEG event onset and offset, the neuronal activation is sustained at a constant level throughout the epileptiform event and that associated fMRI signal changes follow the canonical HRF. The eICA method allows for less constrained analyses capable of detecting early, non-canonical responses. A key step of eICA is the initial deconvolution which can be confounded by various sources of structured noise present in the fMRI signal. To help overcome this, we have extend the eICA procedure by utilizing a fully standalone and automated fMRI de-noising procedure to process the fMRI data from an EEG-fMRI acquisition prior to running eICA. Specifically we first apply ICA to the entire fMRI time-series and use a classifier to remove noise-related components. The automated objective de-noiser, “Spatially Organized Component Klassificator” (SOCK) is used; it has previously been shown to distinguish a substantial fraction of noise from true activation, without rejecting the latter, in resting-state fMRI. A second ICA is then performed, this time on the event-related response estimates derived from the denoised data (according to the usual eICA procedure). We hypothesize that SOCK + eICA has the potential to be more sensitive than eICA alone. We test the effectiveness of SOCK by comparing activation obtained in an eICA analysis of EEG-fMRI data with and without the use of SOCK for 14 patients with rolandic epilepsy who exhibited stereotypical IEDs arising from a focus in the rolandic fissure.

1. Introduction

Event-related functional magnetic resonance imaging (fMRI) is an MRI technique that can be used to detect changes in the Blood Oxygen Level Dependent (BOLD) hemodynamic response to neural activity in response to certain events. The conventional method for detecting event-related responses in fMRI consists of modeling the expected fMRI response to an event by convolving a stimulus presentation time-course with an assumed canonical Haemodynamic Response Function (HRF) and using linear regression to identify voxels with a significant correlation to this expected response (Josephs et al., 1997). One typically assumes that the onset and offset of neuronal activity match stimuli onset and offset, the neuronal activation is sustained at a constant level throughout the stimulus and that evoked fMRI signal changes follow the canonical HRF.

There are instances, however, when these assumptions may not be satisfied. An example is interictal epileptiform discharges (IEDs), which are pathological patterns of activity generated by the brain of patients with epilepsy between seizures (de Curtis et al., 2012). IEDs produce marked and stereotyped trace deviations on electroencephalography (EEG) recordings and can be studied using fMRI by using a simultaneous acquisition of EEG (EEG-fMRI) in order to identify the event timings (Lemieux et al., 2001; Bnar et al., 2002). Studies have shown that the onset of the neuronal activity underlying the EEG discharge may not always coincide with the EEG onset (Bai et al., 2010; Carney et al., 2010; Masterton et al., 2010). For example, Carney et al. (2010) identified changes in BOLD signal which precede the onset of epileptiform activity. In addition, it is also reported that the use of the same HRF in all patients may not be appropriate and that individual-based HRF models provide increases in extent and degree of activation (Masterton et al., 2010; Storti et al., 2013).

To address the above issues, we developed an algorithm, dubbed event-related independent components analysis (eICA), which allowed for less constrained analyses capable of detecting early, non-canonical responses (Masterton et al., 2013a,b). Event-related ICA is a technique that provides an estimate of the underlying components that give rise to the observed event-related fMRI signal changes throughout the brain, and importantly, does not rely upon the specification of an HRF model or predefined Regions of Interest (ROIs). Unlike a standard independent components analysis (ICA), which is applied to the entire fMRI time series, the eICA method is applied only to the event-related time courses at each voxel (an estimate of the event-related signal at each voxel is first obtained by deconvolution of the observed fMRI signal with the observed EEG event timing), which means that only a small number of components are generated that are all explicitly related to the event of interest. Event-related ICA can be applied to data from individual subjects and also to group data using a temporal concatenation approach. We previously demonstrated that the eICA method, when applied to EEG-fMRI data acquired from a group of patients with Benign epilepsy with centro-temporal spikes (BECTS), provided better performance than a standard event-related analysis and a linear deconvolution approach, with a better detection rate in single-subject analyses (73 vs. 53%) and only event-related ICA finding significant group-level activation (Masterton et al., 2013b).

A key element of the eICA is the initial deconvolution. However, the stability of the deconvolution can be compromised by various sources of structured noise (Biswal et al., 1996; Friston et al., 1996; Glover et al., 2000) present in the fMRI signal. These include rapid and slow head movements, physiological activity (breathing and heartbeat) and potential acquisition artifacts. Data driven techniques, especially ICA, are increasingly being employed to separate signal and noise in conventional fMRI data (Thomas et al., 2002; Kochiyama et al., 2005; McKeown et al., 2006; Perlbarg et al., 2007; Stevens et al., 2007; Calhoun et al., 2008; Tohka et al., 2008; Sui et al., 2009; Beckmann, 2012; Kundu et al., 2012; Bhaganagarapu et al., 2013; Salimi-Khorshidi et al., 2014). However, in the context of EEG-fMRI studies in epilepsy, the interpretation of the results from an application of ICA can be difficult as it may produce more than a hundred different components per subject with the majority of these likely having no relationship to the EEG event of interest (Rodionov et al., 2007; LeVan et al., 2010).

To address this we developed a strategy for the automated objective identification of artifactual components from an ICA, that we have dubbed a Spatially Organized Component Klassificator (SOCK) (Bhaganagarapu et al., 2013). The primary objective of SOCK is to distinguish noise from true activation without rejecting the latter. SOCK automatically classifies ICs into one of two categories; artifact or unlikely artifact. It does so using spatial measures likely to indicate motion, physiological noise, or machine or undetermined noise. SOCK was shown to successfully remove artifactual components, without rejecting true activation in resting state data (Bhaganagarapu et al., 2013). Unlike existing automatic classifier methods which are primarly dependent on training data to inform classification (De Martino et al., 2007; Tohka et al., 2008; Salimi-Khorshidi et al., 2014) or require querying a public database (Sochat et al., 2014), SOCK is a standalone, automated and objective method that does not require the user to train the algorithm. It is able to identify a high proportion of artifact-related ICs without removing components that are likely to be of neuronal origin (Bhaganagarapu et al., 2013).

In this paper, we extend the eICA procedure by utilizing SOCK to automatically de-noise fMRI data from an EEG-fMRI acquisition prior to running eICA. As we are de-noising the entire fMRI time series prior to the eICA, we hypothesize that this approach has the potential to be more sensitive than eICA alone. The use of an automated de-noising procedure like SOCK in the context of eICA is a novel methodology and to our knowledge has not been investigated previously. We demonstrate the effectiveness of SOCK by comparing the extent of activation obtained in a standard eICA analysis of EEG-fMRI data with and without the use of SOCK for 14 patients with rolandic epilepsy who exhibited stereotypical IEDs arising from a focus in the rolandic fissure.

2. Methods

2.1. Methods Overview

SOCK is applied to de-noise fMRI data prior to event-related ICA. An overview of the de-noising process is given below (see also Figure 1) and more detail is provided in the sections that follow.

1. ICA was applied to the pre-processed fMRI data (see Section 2.6) using MELODIC (Beckmann and Smith, 2004), yielding both thresholded and unthresholded ICs and associated time courses and power spectra1.

2. ICs were classified into one of two categories using SOCK: artifact or unlikely artifact.

3. All ICs classified into the artifact category were discarded and a de-noised fMRI data series was constructed with only the unlikely artifact ICs.

4. An event-related ICA (eICA) was preformed using the de-noised fMRI data (along with EEG).

The above process was performed for both group and individual studies.

FIGURE 1
www.frontiersin.org

Figure 1. We assess the performance of SOCK by comparing the activation obtained in an eICA analysis with and without the use of SOCK (for both group and individual studies). In the with-SOCK processing stream, ICA was applied to pre-processed fMRI data yielding spatial component maps with associated time courses and power spectra, SOCK automatically classified ICs into one of two categories; artifact or unlikely artifact. Rejecting all ICs classified into the artifact category, a de-noised fMRI data series is formed that is then processed with eICA.

2.2. ICA Decomposition

In the with-SOCK processing stream, ICA is employed to decompose the 4D fMRI time series into a linear combination of spatially independent component maps with an associated time-course (McKeown et al., 1998; Hyvärinen, 1999). In practice this decomposition is usually too computationally expensive to perform on raw fMRI data, so a preliminary data reduction step using principal components analysis is applied prior to ICA. Several freely available software packages are available to perform this preprocessing and decomposition; we used MELODIC which is part of the FSL package (Beckmann and Smith, 2004). The output is a set of spatial maps with associated time courses and power spectra. These then form the input for the automatic classifier, SOCK.

2.3. Classification of ICs Using SOCK

SOCK classifies ICs using features likely to indicate motion, physiological noise, or machine or undetermined noise. The algorithm is described in detail elsewhere (Bhaganagarapu et al., 2013) and our implementation is freely available at www.brain.org.au/software. Briefly, individual slices in each IC are assessed for:

1. Smoothness: contributions of low and high spatial frequency content, to detect components with a large number of isolated very small clusters or isolated voxels (i.e., a spotty appearance).

2. Edge activity: extent of activity in an edge mask.

3. Ventricular activity: extent of activity in a Cerebrospinal fluid (CSF) mask.

4. Temporal Frequency Noise (TFN): the power in temporal frequency beyond 0.08 Hz.

Based on these measures and with the assistance of k-means clustering, ICs dominated by artifact are classified into an Artifact category and all other ICs (i.e., those containing possible neuronal signal) into an Unlikely Artifact category.

The SOCK procedure was implemented in MATLAB (R2010b, The MathWorks, Natick, MA, USA). Source code of our implementation of the method is available at http://www.brain.org.au/software.

2.4. Constructing De-Noised Data

After SOCK classification, all ICs classified in the artifact category are discarded from the fMRI data set and a de-noised fMRI data set is assembled from the remaining components. This is done via the FSL function, fsl_regfilt (with the ‘-a’ aggressive filtering option). The de-noised fMRI data along with the original EEG timings are then input to an eICA analysis.

2.5. Event-Related Independent Components Analysis (eICA)

The eICA method, described in detail elsewhere (Masterton et al., 2013b), can be applied at either an individual or group level. In brief, eICA uses two separate steps to identify events observed in the EEG: firstly, a linear deconvolution (via GLM) provides an estimate of the event-related BOLD response at each voxel in the brain in a time-window spanning from 30 s before to 30 s after the event onset. The deconvolution does not assume any particular response shape and allows for changes occurring before the event onset. ICA is then used to separate the estimated event-related fMRI signal changes into a small number of spatial maps and associated time-courses that summarize the timing of activity within different spatial sources.

To estimate the event-related response across the group, the ICA decomposition was performed upon temporally concatenated data (Calhoun et al., 2009); note that in this context the event related responses (rather then the original fMRI time series) were concatenated. This provided a common set of spatial maps for each group with subject specific time courses. Components of interest were identified as those exhibiting activity in the vicinity of the ipsilateral rolandic region (Masterton et al., 2013b).

The eICA procedure was implemented in MATLAB (R2010b, The MathWorks, Natick, MA, USA) using the SPM8 software (Wellcome Department of Cognitive Neurology, http://www.fil.ion.ucl.ac.uk/spm) to perform the GLM parameter estimation, and the FastICA and ICASSO (Hyvärinen, 1999; Himberg et al., 2004) toolboxes to perform the ICA decomposition. The resulting spatial maps were transformed into z-statistics maps by fitting a mixture model to the data (see Masterton et al., 2013b for more details). This eICA procedure was applied to EEG-fMRI data in separate analyses with and without de-noising the fMRI data via SOCK.

2.6. fMRI Data

2.6.1. Subjects

The fMRI data used for this study was the same as that previously studied with eICA (without SOCK) and described in detail by Masterton et al. (2013b). We summarize key subject details below and in Table 1. Data from fourteen patients with typical BECTS, recruited for EEG-fMRI from the Royal Childrens Hospital, Monash Medical Centre and Austin Hospital in Melbourne, Australia, are included in the analysis. One patient had independent left and right-sided CTS; the remainder had unilateral discharges—this provided a total of fifteen different events for study. More detail on the patient cohort is provided in Lillywhite et al. (2009). A representative EEG recording of CTS discharges in the MRI scanner is also provided in Masterton et al. (2010). This cohort was chosen because the previously published eICA could be used as a gold standard when assessing the performance of SOCK to de-noise fMRI data. This study had approval from the Human Research Ethics Committee at each recruiting hospital and all subjects (or their parents) provided written informed consent.

TABLE 1
www.frontiersin.org

Table 1. Patient details.

2.6.2. Data acquisition

The patients underwent 30 min of simultaneous EEG and fMRI scanning. fMRI images were acquired in a 3T GE Signa LX scanner (General Electric, Milwaukee, WI, USA) using a BOLD-weighted gradient-recalled echo-planar imaging sequence (TR = 3 s; TE = 40 ms; FOV = 24 × 24 cm; 128 × 128 matrix; 25 interleaved 4 mm slices with 1 mm gap). In three studies (Subjects 3, 7, and 9 in Table 1) a slightly different fMRI acquisition was used (TR = 3 s; TE = 40 ms; FOV = 22 × 22 cm; 64 × 64 matrix; 35 interleaved 3.2 mm slices with 0.2 mm gap).

Simultaneous EEG was acquired during fMRI scanning using an MR-compatible EEG system (developed in-house) with scalp electrodes positioned in the standard 10–20 locations and filtering to remove the effect of cardioballistic and motion artifacts (Masterton et al., 2007). The patients' EEG was reviewed by experienced electroencephalographers according to the guidelines developed in our group (Flanagan et al., 2009) and the timing of all identified CTS was recorded.

2.6.3. Data analysis

Image conversion was performed using iBrain (Abbott and Jackson, 2001), preprocessing and statistical analysis utilized SPM8 with the aid of the iBrain Analysis Toolbox for SPM (Abbott et al., 2011; www.brain.org.au/software). Preprocessing included temporal alignment of slices within each volume to the first slice, rigid-body spatial realignment to correct for subject motion, spatial normalization to a symmetric template and spatial smoothing with a Gaussian kernel (FWHM = 8 mm). The symmetric template was created specifically for this patient group using SPM8 software by normalizing each subject's brain to MNI space, averaging these images together (along with a left-right flipped version of each image), and then smoothing with an 8 mm Gaussian filter (Wilke et al., 2002). To enable grouping of data between subjects with left and right-sided CTS, the data from subjects with right-sided CTS were flipped in the left-right direction prior to group analysis.

3. Results

3.1. ICA Analysis and SOCK Classification

MELODIC ICA yielded an average of 114 components per subject (range: 57–285). SOCK classified between 27 and 53% of each subject's components as artifact (mean 41%). These ICs were discarded to construct a de-noised fMRI data set for each subject. See Table 2 for summary details of the ICA decomposition and the SOCK classification for all 14 subjects.

TABLE 2
www.frontiersin.org

Table 2. ICA decomposition and the SOCK classification for 14 patients who underwent an EEG-fMRI study as described in Section 2.6.2.

3.2. Individual Subject Event-Related Analysis

Fifteen individual analyses were performed for eICA, each with and without the use of SOCK to de-noise fMRI data. This included two analyses for the patient (Subject 5) that had independent left and right-sided CTS, which were analyzed as separate events. The results are summarized in Table 3.

TABLE 3
www.frontiersin.org

Table 3. A summary of the results of the individual analyses comparing the number of ICs yielded from an eICA with and without the use of SOCK prior to eICA.

3.3. eICA (without SOCK)

eICA estimated an average of 7 (range: 2–11) different components for each analysis. In 11 out of 15 analyses (73%) at least one component was identified showing activity in the vicinity of the ipsilateral rolandic region.

3.4. eICA (with SOCK)

After de-noising the fMRI data with SOCK, eICA estimated an average of 6 (range: 1–10) different components for each analysis. In 13 out of 15 analyses (87%) at least one component was identified showing activity in the vicinity of the ipsilateral rolandic region.

Figure 2 displays sample slices and time-courses of the rolandic component derived from an eICA of subject 11, with and without SOCK. Activation is overlaid onto a mean functional image for that subject. Comparison of the left (without SOCK) and right (with SOCK) columns shows more robust activation in the area of interest after de-noising the data with SOCK and additional activation on the contralateral side (see green arrows). Furthermore, the shape of the time-course after applying SOCK is also qualitatively smoother.

FIGURE 2
www.frontiersin.org

Figure 2. Spatial maps (thresholded at p < 0.05) with time courses for subject 11 indicating the differences in activation with and without the use of SOCK (left and right columns respectively). Activation is overlaid onto a mean functional image for this subject. Warm and cool colors indicate respectively a positive or negative correlation with the component time course. Arrows in green show areas of increased activation within the region of interest when SOCK was used. Furthermore, the shape of the time-course after applying SOCK is also qualitatively smoother then prior to using SOCK. The zero time-point, indicated by the vertical dotted line in the center of the time-course plot, represents the onset time of the EEG discharge.

The associated spatial maps (with and without SOCK) and time courses for all remaining subjects are provided in Figure 3.

FIGURE 3
www.frontiersin.org

Figure 3. Spatial maps (thresholded at p < 0.05) with time courses for all subjects (except subject 11 already shown in Figure 2) indicating the differences in activation with and without the use of SOCK (left and right columns respectively). Activation is overlaid onto a mean functional image for each subject. Warm and cool colors indicate respectively a positive or negative correlation with the component time course. The time course axes are similar to Figure 2. The zero time-point, indicated by the vertical dotted line in the center of each time-course plot, represents the onset time of the EEG discharge. Two subjects (8 and 10) yielded activation in the vicinity of the ipsilateral rolandic region when analyzed with eICA after de-noising with SOCK but not when analyzed by eICA alone. Furthermore, the shape of the time-course after applying SOCK is also qualitatively smoother then prior to using SOCK.

In two of these analyses (Subjects 8 and 10) activation in the vicinity of the ipsilateral rolandic region was identified when analyzed with eICA after de-noising with SOCK but not when analyzed by eICA alone (see Figure 3). The shape of the time-courses after applying SOCK for these subjects was also qualitatively consistent with the other subjects' peri-rolandic component time courses.

Furthermore, using SOCK prior to running an eICA has qualitatively decreased the noise in both the spatial maps and time courses. For example, the spatial maps for Subjects 9 and 13 (Figure 3) are observed to have little or no activation on the edge of the brain and in the CSF after applying SOCK. In addition, the time courses are observed to follow a BOLD response more consistent with the other subjects.

3.5. Group Event-Related Analysis

3.6. eICA (without SOCK)

eICA group analysis (without SOCK) estimated 14 different components out of which a single CTS-related component of interest was identified low in the ipsilateral post-central gyrus, extending along its opercular surface into the lateral fissure. A much smaller activation was also seen on the contralateral side (see “Without SOCK” panel in Figure 4). It is important to note that the term “activation” used here is defined based upon the direction of signal change near the EEG event onset at time 0; however if the haemodynamic contribution to the response is assumed canonical then the larger than canonical post event dip suggests there is substantial neuronal deactivation following an initial smaller positive neuronal activation event in this component.

FIGURE 4
www.frontiersin.org

Figure 4. The group eICA result. Left panel: Only a single group component of interest was found for the eICA only analysis (Without SOCK), which is dominated by activity low in the ipsilateral post-central gyrus, extending along its opercular surface into the lateral fissure. A much smaller region was also seen on the contralateral side. Right panel: The SOCK+eICA analysis separated these regions into distinct components with different time-courses and revealed a substantially larger extent of de-activation on the contralateral side (see far right component) that appears to have a time course somewhat delayed from the ipsilateral-only component. Arrows in green highlight areas of substantial improvement when SOCK was used. The components are displayed as z-statistic maps thresholded at p < 0.05 corrected for multiple comparisons and overlaid upon the group-mean fMRI image. The time-course at the bottom represents the average modulation of this network across all the subjects i.e., the event-related impulse response function. The zero time-point, indicated by the vertical dotted line in the center of the plot, represents the onset time of the EEG discharge. Warm and cool colors indicate respectively a positive or negative correlation with the component time course.

3.7. eICA (with SOCK)

The SOCK+eICA analysis yielded 15 different components out of which two CTS-related components of interest were identified; (1) low in the ipsilateral post-central gyrus, extending along its opercular surface into the lateral fissure (also observed in eICA without SOCK analysis) and (2) a bilateral component containing de-activation in both the ipsilateral region above and contralaterally (see Figure 4).

4. Discussion

In this paper we have demonstrated the superiority of SOCK+eICA compared with eICA alone for mapping functional brain activity associated with epileptic spikes. It has previously been demonstrated that that eICA is superior to conventional event-related analyses when the BOLD response does not closely match the canonical haemodynamic response function (HRF) (Masterton et al., 2013b). Taken together, our results suggest that SOCK+eICA should replace eICA alone as the preferred method for such analyses.

The centro-temporal spikes of Rolandic epilepsy served as a good test case for our analysis methodology, as it is known that the BOLD response does not well match the standard HRF. This is in part due to neuronal activity associated with the spikes being detectable with fMRI before activity becomes sufficiently widespread and synchronized to manifest as a spike visible on the EEG, and in part due to a larger post-spike undershoot (Masterton et al., 2010, 2013b). The eICA procedure provides a less constrained approach than a GLM incorporating a conventional HRF model, however this flexibility comes at the cost of lower power (increased susceptibility to noise). Whilst the event-related nature of the eICA approach provides a much stronger constraint than conventional ICA on the full fMRI time series, the eICA method is still susceptible to noise, particularly in the initial deconvolution step. Temporally non-stationary noise would be expected to increase the heterogeneity of the raw signal response associated with events, making deconvolution more challenging, and spatial non-stationary of the noise would be expected to increase heterogeneity of the derived event-related responses across voxels. This would then deleteriously affect the performance of the subsequent ICA of the event-related responses. Thus using a procedure which removes a substantial quantity of noise from the input data may improve the end result. The SOCK procedure is a standalone, automated and objective method which is able to remove a substantial fraction of noise without removing biologically interesting signal (Bhaganagarapu et al., 2013). The results of the present study indicate in practice the improvement can be substantial when SOCK is used to de-noise fMRI data prior to eICA.

Applying SOCK+eICA to the existing EEG-fMRI of our BECTS cohort has improved the confidence in the initial results of Masterton et al. (2010) and Masterton et al. (2013b) (i.e., that the centro-temporal spikes arise from low in the ipsilateral post-central gyrus, extending along its opercular surface into the lateral fissure), with two of the previously negative-result individuals now showing activity in this region. There now remain just two individuals with negative results: These two subjects (3 and 6 in Table 1) along with subject 14 registered the smallest number of events compared to all other subjects. The lower the event count, the less power one has to detect an effect (Flanagan et al., 2009).

In these BECTS subjects the epileptic spikes were unilateral during the imaging session. Homologous regions of cortex are connected via fibers projecting through the corpus callosum and typically inhibit neural activity in the contralateral hemisphere. The new SOCK+eICA group analysis reveals a more complex response than previously evident, with initial ipsilateral activity, followed by a more extensive bilateral pattern of deactivation (i.e., the time-course of the deactivation component displays a later rise and peak compared to the ipsilateral-only component). We would interpret the new finding as distinguishing the activation of the ipsilateral cortex during epileptiform events, and a later bilateral decrease in activity in response to this activation.

5. Conclusion

We have demonstrated a novel application of our ICA classifier, SOCK, to de-noise fMRI prior to an event-related ICA in patients with rolandic epilepsy. The procedure outlined in this paper harnesses the advantage of both techniques: (1) SOCK de-noises fMRI in an objective and automated manner utilizing the entire fMRI time-series. (2) eICA utilizes the EEG information to derive event-related responses which are input into an ICA, thus constraining the final eICA decomposition to a small number of components time-locked to the events of interest. The use of SOCK increased power to detect activity of interest in both individual and group analyses.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This study was supported by the National Health and Medical Research Council of Australia (Project grants 368650 and 318900, Program Grant 628952, and a practitioner fellowship 527800 to Graeme D. Jackson), the Austin Hospital Medical Research Foundation, and the Operational Infrastructure Support Program of the State Government of Victoria, Australia.

Footnotes

1. ^Temporal information expressed in the frequency domain. This is done mathematically by taking the discrete Fourier Transform of the time course.

References

Abbott, D., Masterton, R., Waites, A., Bhaganagarapu, K., Pell, G., Harvey, M., et al. (2011). “The iBrain™ analysis toolbox for SPM,” in Proceedings of the 17th Annual Meeting of the Organisation for Human Brain Mapping, Quebec, QC.

Abbott, D. F., and Jackson, G. D. (2001). iBrain software for analysis of visualization of functional MR images. Neuroimage 13, s59. doi: 10.1016/S1053-8119(01)91402-8

CrossRef Full Text

Bai, X., Vestal, M., Berman, R., Negishi, M., Spann, M., Vega, C., et al. (2010). Dynamic time course of typical childhood absence seizures: EEG, behavior, and functional magnetic resonance imaging. J. Neurosci. 30, 5884–5893. doi: 10.1523/JNEUROSCI.5101-09.2010

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Beckmann, C., and Smith, S. (2004). Probabilistic independent component analysis for functional magnetic resonance imaging. IEEE Trans. Med. Imaging 23, 137–152. doi: 10.1109/TMI.2003.822821

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Beckmann, C. F. (2012). Modelling with independent components. Neuroimage 62, 891–901. doi: 10.1016/j.neuroimage.2012.02.020

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Bhaganagarapu, K., Jackson, G. D., and Abbott, D. F. (2013). An automated method for identifying artifact in independent component analysis of resting-state fMRI. Front. Hum. Neurosci. 7:343. doi: 10.3389/fnhum.2013.00343

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Biswal, B., Deyoe, E. A., and Hyde, J. S. (1996). Reduction of physiological fluctuations in fMRI using digital filters. Magn. Reson. Med. 35, 107–113. doi: 10.1002/mrm.1910350114

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Bnar, C.-G., Gross, D. W., Wang, Y., Petre, V., Pike, B., Dubeau, F., et al. (2002). The BOLD response to interictal epileptiform discharges. Neuroimage 17, 1182–1192. doi: 10.1006/nimg.2002.1164

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Calhoun, V. D., Liu, J., and Adali, T. (2009). A review of group ICA for fMRI data and ICA for joint inference of imaging, genetic, and ERP data. Neuroimage 45(1 Suppl):S163–S172. doi: 10.1016/j.neuroimage.2008.10.057

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Calhoun, V. D., Maciejewski, P. K., Pearlson, G. D., and Kiehl, K. A. (2008). Temporal lobe and default hemodynamic brain modes discriminate between schizophrenia and bipolar disorder. Hum. Brain Mapp. 29, 1265–1275. doi: 10.1002/hbm.20463

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Carney, P. W., Masterton, R. A. J., Harvey, A. S., Scheffer, I. E., Berkovic, S. F., and Jackson, G. D. (2010). The core network in absence epilepsy. differences in cortical and thalamic BOLD response. Neurology 75, 904–911. doi: 10.1212/WNL.0b013e3181f11c06

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

de Curtis, M., Jefferys, J. G. R., and Avoli, M. (2012). “Interictal epileptiform discharges in partial epilepsy: complex neurobiological mechanisms based on experimental and clinical evidence,” in Jasper's Basic Mechanisms of the Epilepsies, 4th Edn National Center for Biotechnology Information (US), eds J. L. Noebels, M. Avoli, M. A. Rogawski, R. W. Olsen, and A. V. Delgado-Escueta (Bethesda, MD: Oxford University Press), 213–223.

De Martino, F., Gentile, F., Esposito, F., Balsi, M., Di Salle, F., Goebel, R., et al. (2007). Classification of fMRI independent components using IC-fingerprints and support vector machine classifiers. Neuroimage 34, 177–194. doi: 10.1016/j.neuroimage.2006.08.041

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Flanagan, D., Abbott, D. F., and Jackson, G. D. (2009). How wrong can we be? the effect of inaccurate mark-up of EEG/fMRI studies in epilepsy. Clin. Neurophysiol. 120, 1637–1647. doi: 10.1016/j.clinph.2009.04.025

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Friston, K. J., Williams, S., Howard, R., Frackowiak, R. S., and Turner, R. (1996). Movement-related effects in fMRI time-series. Magn. Reson. Med. 35, 346–355. doi: 10.1002/mrm.1910350312

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Glover, G. H., Li, T. Q., and Ress, D. (2000). Image-based method for retrospective correction of physiological motion effects in fMRI: RETROICOR. Magn. Reson. Med. 44, 162–167. doi: 10.1002/1522-2594(200007)44:1<162::AID-MRM23>3.0.CO;2-E

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Himberg, J., Hyvärinen, A., and Esposito, F. (2004). Validating the independent components of neuroimaging time series via clustering and visualization. Neuroimage 22, 1214–1222. doi: 10.1016/j.neuroimage.2004.03.027

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Hyvärinen, A. (1999). Fast and robust fixed-point algorithms for independent component analysis. IEEE Trans. Neural Netw. 10, 626–634. doi: 10.1109/72.761722

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Josephs, O., Turner, R., and Friston, K. (1997). Event-related fMRI. Hum. Brain Mapp. 5, 243–248. doi: 10.1002/(SICI)1097-0193(1997)5:4<243::AID-HBM7>3.0.CO;2-3

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kochiyama, T., Morita, T., Okada, T., Yonekura, Y., Matsumura, M., and Sadato, N. (2005). Removing the effects of task-related motion using independent-component analysis. Neuroimage 25, 802–814. doi: 10.1016/j.neuroimage.2004.12.027

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Kundu, P., Inati, S. J., Evans, J. W., Luh, W.-M., and Bandettini, P. A. (2012). Differentiating BOLD and non-BOLD signals in fMRI time series using multi-echo EPI. Neuroimage 60, 1759–1770. doi: 10.1016/j.neuroimage.2011.12.028

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lemieux, L., Salek-Haddadi, A., Josephs, O., Allen, P., Toms, N., Scott, C., et al. (2001). Event-related fMRI with simultaneous and continuous EEG: description of the method and initial case report. Neuroimage 14, 780–787. doi: 10.1006/nimg.2001.0853

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

LeVan, P., Tyvaert, L., Moeller, F., and Gotman, J. (2010). Independent component analysis reveals dynamic ictal BOLD responses in EEG-fMRI data from focal epilepsy patients. Neuroimage 49, 366–378. doi: 10.1016/j.neuroimage.2009.07.064

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Lillywhite, L. M., Saling, M. M., Simon Harvey, A., Abbott, D. F., Archer, J. S., Vears, D. F., et al. (2009). Neuropsychological and functional MRI studies provide converging evidence of anterior language dysfunction in BECTS. Epilepsia 50, 2276–2284. doi: 10.1111/j.1528-1167.2009.02065.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Masterton, R. A. J., Abbott, D. F., Fleming, S. W., and Jackson, G. D. (2007). Measurement and reduction of motion and ballistocardiogram artefacts from simultaneous EEG and fMRI recordings. Neuroimage 37, 202–211. doi: 10.1016/j.neuroimage.2007.02.060

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Masterton, R. A. J., Carney, P. W., Abbott, D. F., and Jackson, G. D. (2013a). Absence epilepsy subnetworks revealed by event-related independent components analysis of functional magnetic resonance imaging. Epilepsia 54, 801–808. doi: 10.1111/epi.12163

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Masterton, R. A. J., Harvey, A. S., Archer, J. S., Lillywhite, L. M., Abbott, D. F., Scheffer, I. E., et al. (2010). Focal epileptiform spikes do not show a canonical BOLD response in patients with benign rolandic epilepsy (BECTS). Neuroimage 51, 252–260. doi: 10.1016/j.neuroimage.2010.01.109

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Masterton, R. A. J., Jackson, G. D., and Abbott, D. F. (2013b). Mapping brain activity using event-related independent components analysis (eICA): specific advantages for EEG-fMRI. Neuroimage 70, 164–174. doi: 10.1016/j.neuroimage.2012.12.025

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

McKeown, M., Hu, Y.-J., and Wang, Z. (2006). “ICA denoising for event-related fMRI studies,” in In Engineering in Medicine and Biology Society, 2005. IEEE-EMBS 2005. 27th Annual International Conference (Shanghai), 157–161.

McKeown, M. J., Makeig, S., Brown, G. G., Jung, T. P., Kindermann, S. S., Bell, A. J., et al. (1998). Analysis of fMRI data by blind separation into independent spatial components. Hum. Brain Mapp. 6, 160–188. doi: 10.1002/(SICI)1097-0193(1998)6:3<160::AID-HBM5>3.0.CO;2-1

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Perlbarg, V., Bellec, P., Anton, J.-L., Plgrini-Issac, M., Doyon, J., and Benali, H. (2007). CORSICA: correction of structured noise in fMRI by automatic identification of ICA components. Magn. Reson. Imaging 25, 35–46. doi: 10.1016/j.mri.2006.09.042

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Rodionov, R., De Martino, F., Laufs, H., Carmichael, D. W., Formisano, E., Walker, M., et al. (2007). Independent component analysis of interictal fMRI in focal epilepsy: comparison with general linear model-based EEG-correlated fMRI. Neuroimage 38, 488–500. doi: 10.1016/j.neuroimage.2007.08.003

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Salimi-Khorshidi, G., Douaud, G., Beckmann, C. F., Glasser, M. F., Griffanti, L., and Smith, S. M. (2014). Automatic denoising of functional MRI data: combining independent component analysis and hierarchical fusion of classifiers. Neuroimage 90, 449–468. doi: 10.1016/j.neuroimage.2013.11.046

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sochat, V., Supekar, K., Bustillo, J., Calhoun, V., Turner, J. A., and Rubin, D. L. (2014). A robust classifier to distinguish noise from fMRI independent components. PLoS ONE 9:e95493. doi: 10.1371/journal.pone.0095493

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Stevens, M. C., Kiehl, K. A., Pearlson, G., and Calhoun, V. D. (2007). Functional neural circuits for mental timekeeping. Hum. Brain Mapp. 28, 394–408. doi: 10.1002/hbm.20285

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Storti, S., Formaggio, E., Bertoldo, A., Manganotti, P., Fiaschi, A., and Toffolo, G. (2013). Modelling hemodynamic response function in epilepsy. Clin. Neurophysiol. 124, 2108–2118. doi: 10.1016/j.clinph.2013.05.024

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Sui, J., Adali, T., Pearlson, G. D., and Calhoun, V. D. (2009). An ICA-based method for the identification of optimal FMRI features and components using combined group-discriminative techniques. Neuroimage 46, 73–86. doi: 10.1016/j.neuroimage.2009.01.026

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Thomas, C. G., Harshman, R. A., and Menon, R. S. (2002). Noise reduction in BOLD-Based fMRI using component analysis. Neuroimage 17, 1521–1537. doi: 10.1006/nimg.2002.1200

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Tohka, J., Foerde, K., Aron, A. R., Tom, S. M., Toga, A. W., and Poldrack, R. A. (2008). Automatic independent component labeling for artifact removal in fMRI. Neuroimage 39, 1227–1245. doi: 10.1016/j.neuroimage.2007.10.013

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Wilke, M., Schmithorst, V. J., and Holland, S. K. (2002). Assessment of spatial normalization of whole-brain magnetic resonance images in children. Hum. Brain Mapp. 17, 48–60. doi: 10.1002/hbm.10053

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Keywords: functional magnetic resonance imaging (fMRI), independent component analysis (ICA), automated classification, artifacts, denoising, filter, event related ICA, Benign epilepsy with centro-temporal spikes (BECTS)

Citation: Bhaganagarapu K, Jackson GD and Abbott DF (2014) De-noising with a SOCK can improve the performance of event-related ICA. Front. Neurosci. 8:285. doi: 10.3389/fnins.2014.00285

Received: 30 June 2014; Accepted: 22 August 2014;
Published online: 19 September 2014.

Edited by:

Stephen C. Strother, University of Toronto, Canada

Reviewed by:

Michael Hanke, Otto-von-Guericke-University, Germany
Christian Windischberger, Medical University of Vienna, Austria

Copyright © 2014 Bhaganagarapu, Jackson and Abbott. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: David F. Abbott, The Florey Institute of Neuroscience and Mental Health, Melbourne Brain Centre - Austin Campus, 245 Burgundy Street, Heidelberg, VIC 3084, Australia e-mail: david.abbott@florey.edu.au

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.