Functional near-infrared spectroscopy (fNIRS) is a noninvasive neuroimaging technique that monitors tissue oxygenation and cortical hemodynamics in the brain.1 Near-infrared light with wavelengths within the 650- to 900-nm range is projected onto the scalp by optical emitters, and is mainly absorbed by hemoglobin during its propagation in brain tissue. By collecting the scattered light with sensors positioned a few centimeters from the emitters, local changes in blood oxygenation on the cortex can be measured in terms of concentration changes in oxygenated hemoglobin (HbO) and deoxygenated hemoglobin (HbR).23.–4 By adding together the concentration variations in HbR and HbO, concentration changes in total hemoglobin (HbT) are obtained, which can be regarded as an indicator of regional cerebral blood volume (rCBV) variations under the hypothesis of constant hematocrit.5 Compared to alternatives such as single photon-emission computed tomography (SPECT), positron emission tomography (PET), or functional magnetic resonance imaging (fMRI), fNIRS offers the potential for long-term, noninvasive monitoring, and higher temporal resolution hemodynamic imaging.
fNIRS applications in the study of epilepsy have grown in recent years.6 Epilepsy is a common and diverse set of chronic neurological disorders characterized by seizures (or ictus), resulting from abnormal and excessive neuronal discharges.7,8 Seizures last from several seconds to a few minutes, and are usually accompanied by clinical manifestations such as altered mental state, tonic and/or clonic movements, and various other experiential symptoms.9,10 Brief paroxysmal discharges (20 to 200 ms) may occur in between these ictal events (i.e., seizures), and are called “interictal epileptiform discharges” (IEDs, a.k.a “spikes”). Both seizures and IEDs can usually be captured and observed with electroencephalography (EEG). Using fNIRS concurrently with EEG enables researchers not only to detect and monitor epileptic events, but also to study their impact on hemodynamics and oxygenation.
In past decades, different groups have assessed the clinical potential of EEG-fNIRS as an auxiliary method of epilepsy diagnosis1220.127.116.11.18.104.22.168.22.214.171.124.126.96.36.199.188.8.131.52.32.–33 (summarized in Table 1). In these works, the hemodynamic effects of epileptic events (mostly seizures) were preliminarily investigated in adults, children, and neonates. Early studies revealed heterogeneous patterns of oxygenation changes that might be specific to seizure types.11,12,1819.–20 For example, Villringer et al.20 reported cerebral blood volume increases concordant with seizure origin in three patients with presumed extratemporal epilepsy, while a reproducible deoxygenation was also observed over the frontal lobe by Steinhoff et al.19 on two patients with mesial temporal lobe epilepsy. Sokol et al.18 studied eight patients with drug refractory temporal epilepsy using fNIRS and showed distinct patterns of cerebral oxygenation for complex partial seizures and for rapidly secondarily generalized complex partial seizures. Haginoya et al.12 studied various types of pediatric epileptic seizures with fNIRS in 15 children and reported heterogeneous CBV changes depending on seizure types (convulsions, absence, tonic status epilepticus, and infantile spasms). These initial observations suggested the utility of fNIRS in epilepsy research but did not make direct progress toward evaluating its preclinical value for focus lateralization or localization because of the limited number of subjects, the use of a single fNIRS channel (frequently affixed to the hairless skin overlying the frontal cortex to avoid hair contamination) distant from the seizure focus, heterogeneity in seizure types, and/or lack of clear confirmation of the epileptogenic zone. One pioneer fNIRS study on focus lateralization was carried out by Watanabe et al.,23 where they performed bihemispheric sampling of brain hemodynamics and correctly lateralized the focus for 28 of the 29 patients. Their recent work on mesial temporal lobe epilepsy presented four successful focus lateralization cases when using HbO changes as the index.31 However, the behavior of other chromophores was not discussed in their reports and information on remote changes was not elaborated.
Summary of the literature on EEG-fNIRS in human epilepsy.
|Work||Number of patients||Number of channels||Epilepsy type||Major contributions|
|Villringer et al.20||3||1||Focal (extratemporal)||Increases in blood volume and HbO concentration were found over the presumed epileptic focus region|
|Steinhoff et al.19||2||1||Focal (temporal)||Reproducible cerebral deoxygenation was detected in the ipsilateral frontal cortex|
|Adelson et al.11||3||1||Focal (foci not specified)||A perfusion-metabolism mismatch was observed during seizures|
|Sokol et al.18||8||1||Focal (temporal)||Complex partial seizures and rapidly secondarily generalized complex partial seizures showed distinct pattern of cerebral oxygenation|
|Watanabe et al.21||12||8 to 24||Focal (10 temporal, 2 parietal)||In all 12 cases, rCBV increased rapidly on the focus side after partial seizure onsets|
|Haginoya et al.12||15||1||Cryptogenic or symptomatic, generalized or focal||Heterogeneous CBV changes were reported depending on seizure types|
|Watanabe et al.23||29||8 to 24||Focal (25 temporal, 2 frontal, 2 parietal)||Significant hyperperfusion on the side of seizure foci was observed in 96% of cases (28/29)|
|Buchheim et al.24||3||1||Generalized (absence)||Absence seizures were associated with a decrease in HbO and an increase in HbR over the frontal cortex|
|Munakata et al.13||3||24||Generalized (west syndrome)||An increase in rCBV occurred in multiple areas during clusters of spasms|
|Shuhaiber et al.16||1||2||Focal (frontal)||Seizures produced a progressive decline and wide fluctuations in regional oxygen saturation () on the side ipsilateral to the focus|
|Arca Diaz et al.25||1||2||Focal (frontal)||Fluctuations in were observed in bilateral frontal regions during concurrent seizures|
|Roche-Labarbe et al.14||6||1||Generalized||Generalized spike-and-wave discharges were related to complex oxygenation changes in frontal areas|
|Wallois et al.28||1||2||Focal (temporal, neonatal seizures)||Seizure-like discharges induced increases first in HbR, then in HbO and HbT; an undershoot was seen only for HbR|
|Cooper et al.26||4||Up to 60||Generalized (neonatal seizures)||Transient hemodynamic events identified were associated first with an increase followed by a rapid decrease of HbO concentration|
|Slone et al.17||5||2||Focal (temporal)||A decrease in HbO, HbT, and oxygen saturation () were observed in the frontal lobe up to 15 min before the seizure onset|
|Sato et al.34||1||48||Focal (frontal)||Increased cerebral blood flow was seen in the seizure onset zone after the onset of a SMA seizure, with rapid propagation to adjacent areas|
|Seyal15||6||1||Focal (temporal)||Preictal increase was observed in the frontal lobe, followed by periictal decrease and postictal increase|
|Singh et al.32||1||58||Neonatal seizures||High amplitude HbO, HbR, and HbT increased prior to extended decreases were seen to be associated with all seven electrographic neonatal seizures|
|Yücel et al.22||2||6 to 7||Focal||An increase in HbO, HbT, blood flow and metabolic rate of oxygen () was located over the epileptic focus region during seizures. A subsequent deoxygenation phase was also seen for one seizure|
|Jeppesen et al.29||15||2||Focal (10 temporal, 3 frontal, 1 parietal, 1 unspecified)||NIRS-measured hemoglobin concentration changes in the frontal lobe during seizures may not be suitable as biomarkers for seizure detection|
|Rizki et al.31||4||44||Focal (mesial temporal)||EEG-fNIRS captured HbO increases during seizures in mesial temporal lobe epilepsy, which could be used to lateralize the epileptic focus side|
|Sokoloff et al.30||11||2||Neonatal seizures||Neonatal seizures were associated with reduced on both hemispheres|
|Vinette et al.33||4||2||Focal (3 temporal, 1 central)||Slowly varying changes in hemodynamics may precede epileptic seizures|
In the past few years, the feasibility of applying an optical imaging technique in the study of the hemodynamic response to neonatal seizures has been explored by different groups26,28,30,32 (see also Ref. 35 for a review of early discoveries). Their recent results showed that the shape of the hemodynamic changes associated with neonate seizures might be distinct from the standard responses to functional stimulations on older children or adults. Since reading neonatal EEG can be extremely challenging, and performing SPECT, PET, or fMRI on neonates is usually not possible due to ethical and technical restrictions,28 their studies are remarkable in showing the potential of fNIRS for monitoring in neonatal intensive care units.
Another particular interest of recent fNIRS work was in the observation of oxygenation variations preceding the EEG onsets of epileptic seizures. With EEG-fNIRS, two different groups observed preictal oxygenation changes over the frontal lobe prior to temporal lobe seizures.15,17 Although a subsequent attempt of using fNIRS-measured hemoglobin concentration changes as biomarkers for seizure detection was unsuccessful,29 this endeavor suggests a potential for fNIRS as a home monitoring seizure detection or prediction device.
Some other studies focused on the application of recent fNIRS methodological improvements to epilepsy. For example, aiming at removing the movement artifact from clinical fNIRS recordings, Yücel et al.22 proposed to use collodion-fixed prism-based optical fibers, which were eventually proven to be helpful in recovering the hemodynamic response to three epileptic seizures in their study. Serving the same purpose, data-based filtering techniques were also adapted into the analysis of long-term fNIRS data.33 In addition, efforts have been made by Machado et al.36 to provide an optimal fNIRS emitter/detector montage in IED studies to maximize the sampling sensitivity over one or several brain regions.
In summary, preliminary work confirms the clinical potential and usefulness of fNIRS: (1) it was seen as a robust tool to explore hemodynamics associated with seizures on adults, children, and neonates; (2) studies of seizures indicated a potential for epileptic focus localization; (3) it can possibly be used to study preictal hemodynamic changes; and (4) fNIRS hardware and software have steadily improved over the last few years and will most likely continue to improve, which may in turn benefit future clinical studies.
Our group has been trying to complement the current research paradigm by recording from more channels (typically 100 to 150) to reconstruct two-dimensional (2-D) topographic images of activations with much larger brain coverage. The relatively large number of channels ascertained complex local and remote oxygenation changes during focal seizures in a fair number of patients.3738.–39 In addition, we were able to detect the hemodynamic responses associated with IEDs.3940.–41 We also assessed the degree of oxygenation demand due to IEDs in the epileptic focus based on EEG-fNIRS data collected from a total of 40 patients.42
Simultaneous Electroencephalography-Functional Near-Infrared Spectroscopy Recordings
Back in 2006, when one of us Dang Khoa Nguyen (DKN) initiated an EEG-fNIRS epilepsy research program, commercial fNIRS equipment did not incorporate EEG, and recordings were not long enough to study seizures. Helmets were often made with Velcro strips and were very sensitive to patient movement. A motorcycle helmet alternative was uncomfortably heavy. This motivated the design of a rigid and light helmet drilled with holes to mount NIRS optodes (sources and detectors) and EEG electrodes.43 Over the last 5 years, portable, long-term monitoring systems dedicated to epilepsy have been developed.4445.–46
In studies reviewed in the following sections,37184.108.40.206.–42,47 we used a home-modified commercial (nonportable) fNIRS imaging system (ISS Imagent, Champaign, Illinois) with helmets of different sizes to ensure the best fit to patients’ heads. Recording hardware included 19 carbon EEG electrodes, 64 fiber sources, and 16 fiber detectors mounted on the helmets. A detector is placed 3 to 5 cm away from several adjacent sources, forming multiple optical source–detector (SD) channels [Fig. 1(a)]. Optode positions were then coregistered onto a three-dimensional (3-D) high-resolution anatomical MRI (obtained prior to EEG-fNIRS) using BrainSight™ (Rogue Research, Montreal, Canada) [Fig. 1(b)]. 3-D coordinates of an SD channel were considered as the midpoint of its source and detector. In our studies, channel positions were intentionally arranged so that the covered area included the whole lobe with the preacknowledged epileptic focus, the corresponding contralateral lobe, and as much area from the other lobes as possible, especially on the suspected focus side. An alternative approach to generate the SD montage was recently proposed in Ref. 36, in which the authors aimed at maximizing the spatial sensitivity in one or several brain regions.
After installation of the helmet, multiple sessions of simultaneous EEG-fNIRS are recorded for each patient. During recordings, each patient is simply asked to sit comfortably in a chair and relax. Two wavelengths on either side of the 800-nm hemoglobin isosbestic point (i.e., the wavelength at which the absorption rates of HbR and HbO are equal) are used: 690 nm, which is more sensitive to HbR, and 830 nm, which is more sensitive to HbO. Optical intensities are recorded at 19.5 Hz by an oximeter. Hemoglobin concentration changes are then solved from the recorded optical intensities by using the modified Beer–Lambert law.48 Multiple SD channels deliver topographical images (i.e., activation maps on the surface of the brain), which are now commonly seen in the literature.
EEG data are recorded at a sampling rate of 500 Hz using a Neuroscan Synamps 2TM system (Compumedics, Charlotte, North Carolina), and are then band-pass filtered between 0.1 and 100 Hz to remove instrumental noise and other artificial disturbances. After acquisition, the EEG data are reviewed offline by certified clinical neurophysiologists. Seizure-onset and seizure-offset times, defined, respectively, as the earliest and latest clinical or electrographic evidence of seizure activity, are marked directly on the EEG along with IEDs using Analyzer 2.0 (Brain Products, Gilching, Germany). Physiological data are also marked on the electrocardiogram and are semiautomatically filtered.
Two software packages to analyze EEG-fNIRS data have been developed in house. The first one is an extension of HomER49 with helmet visualization and artifact detection. More explicitly, a topographic view onto the scalp is created with MATLAB (MathWorks) to display the evolution and general propagation of HbR and HbO activations. A student’s -test is subsequently performed to assess the significance of activation amplitude over baseline, assuming a confidence level of 95% corrected by Bonferroni thresholds. The calculated -values of all channels are then mapped to the topographic view without interpolation, according to the coordinates of the channels that have been acquired during the coregistration process. Overall, this software package allows quick inspection of fNIRS time courses and simplicity in statistical analysis.
The other data processing package (named nirs10) features more developed 2-D topographic reconstructions, which require a more elaborate statistical treatment for multiple comparisons. It is based on statistical parametric mapping (SPM)50 and allows automated preprocessing and standardized statistical analysis adapted to fNIRS data in a way similar to the treatment of fMRI data.51 A general linear model (GLM) is introduced to model the time-varying brain hemodynamics, as this can provide for more flexible tests of activations than simpler approaches.52,53 The detected hemoglobin concentration changes are decomposed by the GLM into a linear combination of the expected hemodynamic response of epileptic or physiological events (such as seizures, IEDs, cardiac oscillations, etc.) plus an error term. These responses are constructed by convolving the time and duration of the events (which have been previously marked on EEG data) with a basis function, normally a canonical hemodynamic response function (HRF).54 To obtain 2-D contrast maps from fNIRS data, the central difficulty in accounting for multiple comparisons is the sparseness of measurements (unlike fMRI data). Ye et al.5556.–57 recently addressed this problem in their near-infrared spectroscopy-statistical parametric mapping (NIRS-SPM) toolbox. Inhomogeneous interpolation was applied and temporal correlation of the residuals was targeted in a desired way by the “precoloring” process.51 -tests were then conducted at the interpolated measurements for each type of epileptic events. Their methodology was validated in controlled experiments, such as finger-tapping (where strong hemodynamic responses could be expected). The main NIRS-SPM statistical analysis developments were integrated into our nirs10 package.
Multichannel Electroencephalography-Functional Near-Infrared Spectroscopy for Focal Seizure Observation
Epileptic seizures are broadly classified into two different types, according to the source within the brain: focal seizures where the foci lie in only part of a cerebral hemisphere, and generalized seizures where the activation is distributed and involves both cerebral hemispheres.58 The most frequent focal seizures are temporal and frontal lobe seizures, while parietal, occipital, and insular seizures are less common and are thus less well studied in the literature. In this section, we review our work in studying hemodynamic responses associated with seizures during temporal lobe epilepsy,37 frontal lobe epilepsy,38 as well as posterior epilepsy39 with EEG-fNIRS. The data processing was carried out with our first analysis package (based on HomER) in the former two studies, and with the nirs10 toolbox in the latter study of posterior epilepsies. Our studies on seizures and IEDs were approved by the Ethics Committees of Sainte-Justine and Notre-Dame Hospitals. Informed consents were obtained from all subjects. On each subject, a comprehensive neurological evaluation was performed including ictal SPECT, PET, anatomical brain MRI, magnetoencephalography, and intracranial EEG (when needed). Focus localization was carried out by an epileptologist (DKN) based on multimodal analysis, and was confirmed by seizure-freedom following epilepsy surgery for most patients with a follow-up of more than a year.
Multichannel electroencephalography-functional near-infrared spectroscopy recording of temporal lobe seizures
Nine patients with refractory temporal lobe epilepsy were recruited for simultaneous EEG-fNIRS recordings at the Optical Imaging Laboratory of Hôpital Sainte-Justine (Montreal, Quebec, Canada). Eight temporal lobe seizures were successfully recorded in three patients while the other six patients had uneventful EEG-fNIRS sessions in terms of seizures (but not IEDs, analyzed later).
All detected seizures were accompanied by significant alterations in cerebral hemodynamics. Hemoglobin concentration changes occurred almost simultaneously with the suspected electrical onset. Variations from baseline to the highest peak reached 28% for HbR, 13% for HbO, and 11% for HbT. A detailed description of the hemodynamic responses to each seizure can be found in Ref. 37. Here, we first present the results from a group analysis where all eight seizures were aligned to seizure onset to assess overall variation of HbO and HbR (Fig. 2). In this figure, two distinct phases were identified in the focus area: in phase 1, a significant HbR decrease as well as an HbO increase occurred initially (consistent with a compensatory rise of rCBV to augment oxygen supply to discharging epileptic neurons59); in phase 2, after the initial decrease, a longer increase in HbR, which often outlasted the duration of the seizure, was seen for the three patients, while HbO continued to increase for patient 1 and even decreased for patients 2 and 3. This implies that during complex partial seizures, the brain might suffer from hypoxia. The reactive increase in rCBV might not be able to sufficiently meet the local oxygen demand.
Another notable finding in the study was the surprisingly similar hemodynamic profile seen over the region contralateral to the focus. Similar variations, albeit of lower amplitude, in hemoglobin concentration during and after temporal lobe seizures were observed early on in the opposite temporal lobe region. In other remote lobes of the ipsilateral hemisphere, hemoglobin concentrations tended to have changes of the opposite sign. While contralateral homologous and remote changes may be explained by seizure propagation, very early hemodynamic changes appeared to occur without any evidence of contralateral epileptic activity on EEG. It was then concluded that focal temporal lobe seizures had effects on remote regions both before and after they propagated to them.
The results were also evaluated to establish whether EEG-fNIRS could adequately lateralize the epileptic focus. Two laterality indices were calculated by comparing the extent of HbO activation as well as its peak -value between hemispheres. Among the eight temporal lobe seizures, six (from two out of the three patients) showed a larger extent of activation and/or higher peak amplitude on the ipsilateral side of the epileptic focus during the initial HbO increase. Lateralization of the focus to the correct side was possible for these two patients.
The following case study on seizure 4 in patient 1 illustrates these findings (Fig. 3).
Patient 1, seizure 4: This patient is a 49-year-old man with left medial temporal lobe epilepsy due to prior meningitis in childhood. One hundred and twenty optical channels were placed bilaterally over the frontal, temporal, and parietal lobes during EEG-fNIRS recordings. Five complex partial seizures with mean duration of 93 s were recorded during total recording time of 82 min. Figure 3(a) depicts the hemodynamic variations associated with seizure 4 over time in terms of HbO and HbR concentration changes measured by selected optodes over the left (ipsilateral) and right (contralateral) temporal lobes. Figure 3(b) presents uncorrected -maps at several specific times prior to obvious ictal rhythmic theta activity. In the temporal region ipsilateral to the focus side, an increase in HbO as well as an increase in HbT was first noticed (phase 1). The increase in HbO and in HbT mostly resulted from a local increase in rCBV, while HbR decreased. After phase 1, phase 2 followed, consisting of a significant long-duration increase in HbR as well as continued increases in HbO and in HbT. Hence, possible hypoxia over the seizing tissue was suggested by that phase. Ictal hemodynamic changes gradually returned to the baseline after the seizure offset. In the contralateral anterior temporal region, despite no evidence of epileptic contralateral activity on EEG at seizure onset, hemoglobin concentration changes behaved in a very similar way. Moreover, the results of seizure 4 helped to lateralize the epileptogenic focus to the correct hemisphere, as a larger spatial extent in HbO activation was clearly observed on the ipsilateral side.
Multichannel electroencephalography-functional near-infrared spectroscopy recording of frontal lobe seizures
A total of 20 patients with nonlesional, suspected refractory frontal lobe epilepsy also underwent EEG-fNIRS at Hôpital Sainte-Justine. The epileptic focus was eventually clearly identified for 11 patients following clinical investigation. Eighteen seizures from these 11 patients were recorded: five complex partial seizures, five simple partial seizures, and eight electrical seizures.
All detected seizures were observed to be associated with changes in the fNIRS optical signal over the epileptogenic zone. Details of ictal hemodynamic variations for each seizure can be found in Ref. 38. In summary, these variations usually occurred () after seizure onset, and persisted until seizure offset (seven patients) or even beyond (four patients). Over the seizure focus, an initial concentration increase for HbO and HbT from baseline to peak was noted, followed by a gradual decline to a plateau or to the baseline. The amplitude of the rise highly correlated seizure duration. Brief electrical seizures (mean duration: 9.4 s) were accompanied by small increases from the baseline (mean: 1.2%/0.9% for HbO/HbT), while larger increases (mean: 7.4%/5.7% for HbO/HbT) were seen with longer electroclinical seizures (mean duration: 25.3 s). HbR variation within the epilepticus focus was heterogeneous. In four patients, HbR behaved as expected, with an initial decrease to a nadir and a progressive return to the baseline or plateau. In the remaining patients, either no significant response or an opposite response was encountered. Hence, in contrast to what we observed with temporal lobe seizures,37 the local increase in rCBV appeared to be sufficient to meet the metabolic demand in the majority of recorded frontal lobe seizures (9 out of 11 patients) as most showed little change (suggesting sufficient oxygen supply) or a decrease (suggesting oxygen oversupply) in HbR. These observations may probably be due to the fact that recorded frontal lobe seizures were shorter in duration (mean duration: 18.2 s) compared to recorded temporal lobe seizures (mean duration: 66.5 s). Further recording of seizures is necessary to investigate the different factors that may affect the behavior of HbR.
In the homologous contralateral region, we again saw very similar alterations in HbO, HbR, and HbT as in our previous study on temporal lobe epilepsies.37 These closely mirrored activations occurred with no clear delay to ipsilateral hemoglobin variations, but their peak amplitudes were usually slightly lower.
To lateralize the epileptic focus, during the early phase of frontal lobe seizures, EEG-fNIRS discerned larger HbO activation extent on the focus side in 7 out of 8 patients, but higher peak amplitude of HbO activation on the focus side in only 4 out of 8 patients. The extent of HbO activation was more indicative of the epileptogenic hemisphere than peak HbO activation. Additional work is necessary to evaluate the EEG-fNIRS potential in seizure localization.
Multichannel electroencephalography-functional near-infrared spectroscopy recording of seizures from posterior epilepsies
Nine patients with refractory posterior epilepsies (such as occipital, parieto-occipital, or temporo-occipital epilepsies) participated in our EEG-fNIRS study.39 Seizures were seen in three patients, the first of whom had two simple partial seizures (lasting 34 and 52 s, respectively), the second suffered from numerous (34) brief auras of flashing lights but without clear electrographic correlates, and the third had a simple partial seizure (lasting 23 s), again with no obvious EEG change.
In most cases, GLM analysis revealed some responses to recorded epileptic seizures with various degrees of significance and agreement with the patients’ epileptic foci. In the first two patients, significant HbR decreases were located in the focus as well as in the corresponding contralateral areas. In the last patient, we did not observe any HbR decrease in the occipital lobes, but were able to locate an HbO increase broadly consistent with the suspected focus region.
In summary, our work in recent years confirms that continuous fNIRS-EEG can indeed detect hemodynamic changes during focal seizures and that these variations (occurring both locally over the epileptic focus and in remote areas) are more complex than previously thought.
Multichannel Electroencephalography-Functional Near-Infrared Spectroscopy for Interictal Epileptiform Discharge Observation
Results on seizures summarized earlier point to the clinical potential of continuous EEG-fNIRS in detecting focal seizures and monitoring their hemodynamic response before, during, and after their occurrence. Unfortunately, one may not always be able to record seizures during a prolonged EEG-fNIRS study as seizures are unpredictable. It is likely, however, that IEDs will be recorded as they are typically numerous in patients with drug-refractory focal epilepsy. Their analysis may prove useful in identifying epileptic foci or their networks.60,61
In Ref. 41, two analysis methods based on the GLM were used to process fNIRS data related to IEDs: a Bayesian general linear model (BGLM62) and wavelet generalized least square (WGLS63). Both methods were applied in simulation and to the analysis of responses to bilateral fronto-temporal IEDs in a patient with focal right frontal epilepsy. Both methods successfully detected a significant HbO increase and an HbR decrease (thus an increment of rCBV could be inferred) that were spatially concordant with the presumed epileptogenic focus. WGLS was overall more accurate, more specific, and more reliable than the BGLM, mainly because of its ability to model transient events for the drifts and its inherent noise whitening property. However, it was also noticed that both methods became inaccurate in case of low signal-to-noise ratio (SNR) and would lead to a high rate of false positives.
In Ref. 40, IEDs from three more patients were studied, comparing results between EEG-fNIRS and EEG-fMRI in the same subjects. These patients had IEDs that occurred very frequently, with an average of 32/26/26 IEDs/min, respectively. It thus appeared interesting to investigate the extent to which the hemodynamic response was linear when consecutive IEDs occurred within a brief timespan of one another. Clearly, when a large number of spikes occur in rapid succession, an inhibitory nonlinear effect is to be expected, as the hemodynamic response would grow arbitrarily large under a linear assumption. The method of Volterra kernels,64 an analog of the Taylor expansion of series as applied to functions, was used to quantify nonlinearities, implemented in the nirs10 software package. In these three patients, the first-order Volterra (linear response) and second-order Volterra (first nonlinear correction to the linear response) amplitudes were estimated using GLM methodology. Good concordance between fNIRS and fMRI in the three patients was observed for the amplitudes and ratio, and in the correspondence of a positive blood oxygenation level-dependent (BOLD) signal with a decrease of HbR and an increase of HbO. Adding this nonlinear contribution helped to obtain the correct localization of the epileptic focus with fMRI for one patient, which was missed in standard analyses without nonlinearities. For EEG-fNIRS, focus localization was less clear due to the presence of several responses in other areas of the brain, mainly mirror sites to the epileptic focus. Through this work, it was then concluded that including nonlinear terms to GLM might be beneficial when numerous IEDs take place in rapid succession.
As part of the study on posterior epilepsies,39 we applied GLM analysis to 7 out of 9 patients exhibiting focal IEDs during EEG-fNIRS. We located HbR responses and they were broadly concordant with the epileptic focus in 6 out of 7 patients. In five patients (5/7), the HbR response near the focus was recognized to be the most significant among all the activations with respect to -statistic scores and spatial extents, and was thus indicative of the focus. Nonlinearity in HbR response was also significant at a -statistics threshold of 3.0 in 3 out of 4 patients with frequent IEDs. Concentration changes in HbO in these patients were seen to be less sensitive () and specific () with the focus than HbR.
Sensitivity and specificity
In Ref. 42, we performed EEG-fNIRS on a larger sample set of patients (40) with focal epilepsy, to estimate the percentage of patients we can expect EEG-fNIRS to help in activation detection and focus localization. In this study, GLM analysis was performed in a standardized manner as described in Sec. 2.2 with the nirs10 toolbox, assuming an SPM8 canonical HRF50 as the basis function. In each patient, the anatomical MRI was segmented to obtain a gray matter layer, from which four 2-D cortical projections (frontal, dorsal, right, and left) were extracted. Optode positions on the standard SPM8 anatomical atlas were then mapped back to these topographical projections of the patients’ gray matter layer. One-tailed 2-D contrast maps for each IED type from a patient-level analysis (“group-of-sessions”) were finally generated by interpolation of the hemodynamic responses for the four views. As the Euler Characteristic65 and Bonferroni thresholds were a bit too severe at the patient-level, an attempt was made at calculating a peak false discovery rate threshold66 on 2-D maps, assuming a 95% confidence level. In the meantime, the 3-D position of the most plausible epileptic focus, previously marked by an epileptologist, was also projected onto the four 2-D views. Comparisons between the locations of detected activations and the focus were then made at each view.
Five patients were excluded from the sample because of too few IEDs or uncertain epileptic focus position. In the remaining 35 patients, 29 suffered from neocortical epilepsy, while the other six had mesial temporal lobe foci. In neocortical epilepsy, significant negative HbR concentration changes associated with IEDs were seen to be concordant with the epileptic focus in 18 out of 29 patients (sensitivity: 62%), while in eight patients, the HbR decreases near the focus were identified as being the most significant among all other clusters and thus would lead to successful identification of the epileptic focus (specificity: 28%). Similarly, significant HbO or HbT increases near the focus were observed in 38% of the 29 patients, and were the most significant in 21%. Specificity was further improved to 45%/24%/24% for HbR/HbO/HbT if analysis was restricted only to fNIRS activations that were topographically related on the basis of observed epileptiform activity. On the other hand, EEG-fNIRS was generally insensitive in the remaining six patients with mesial temporal lobe epilepsies (sensitivity: 0%/17%/17%, specificity: 0% for all chromophores). Although temporal IEDs detected on scalp EEG indicated that IEDs originating from mesial structures had projected to the temporal neocortex, corresponding hemodynamic changes were not observed in most of these six patients.
The relatively modest sensitivity and specificity rates obtained in this study were probably partially due to the fact that IEDs evoke a weaker neurovascular response compared to seizures. Moreover, we were still able to recognize some responses to IEDs in remote regions, e.g., similar HbR activations in corresponding areas of the contralateral lobe were observed on 61% of the 18 patients who already had significant HbR responses near the focus (50% for HbO or HbT). Further analysis is needed to better understand the pathophysiology (and thus, the metabolic demands) behind hemodynamic propagation.
Analysis of patient 39 is presented here as an individual case study (Fig. 4).
Patient 39, IEDs: This 18-year-old patient had an epileptic focus located in the right supplementary motor area (SMA) [Fig. 4(a)]. Four 15-min sessions followed by a 9.6-min final session were recorded with EEG-fNIRS. Some 121 channels were mounted on the helmet providing wide coverage of bilateral frontal lobes, parietal lobes, and central regions [Fig. 4(b)]; 550 central IEDs were captured during recordings [average: , Fig. 4(c)]. GLM analysis revealed a significant HbR decrease (with a minimum -value of ) located mainly in the right SMA at the time of IEDs [Fig. 4(d)], allowing correct localization of the epileptic focus. Increases in HbO or in HbT were not found in the right SMA, but were observed near corresponding areas contralateral to the focus region on the left side (a bit too posterior, maximum -values: for HbO/HbT).
In recent work,47 we sought to model patient-specific HRFs to IEDs for GLM analysis with fNIRS data (instead of assuming an HRF with a canonical shape as the basis function). For this purpose, we chose to follow the fMRI literature.6768.69.70.–71 Specifically, a second-order Volterra term was also added to the model to account for the nonlinearity in the hemodynamic response when frequent IEDs were present in the data. On five patients with refractory focal epilepsy, we consistently obtained higher statistical scores and/or larger spatial extents for activations using specific HRFs compared to the SPM8 canonical HRF. These preliminary results suggested that, in some cases, including specific HRFs in the analysis can potentially provide better sensitivity in the detection of neuronal activations associated with epileptic events.
In this paper, we briefly reviewed the basic principles of fNIRS before describing in some depth the methods that our group has been using for simultaneous EEG-fNIRS recordings. Through the studies of focal seizures, we reported that the hemodynamic changes can last much longer than the EEG evidence of seizure activity and that epileptic seizures were associated with increased oxygenation (in terms of an increase in HbO) over the epileptic focus in most cases. In addition, early complex hemodynamic changes outside the focus were frequently observed, notably in contralateral homologous areas but also in more remote regions as seizures evolve. HbR concentration behaved heterogeneously: it usually decreased, sometimes stayed unchanged, or even in some cases increased during some seizures, indicating a concomitant deoxygenation process. Along another research direction, using careful modeling, we have also revealed a clinical potential in studying the hemodynamic response to IEDs with EEG-fNIRS in focal epilepsy research.
In recent years, fNIRS studies of epilepsy have provided promising results, which in time may lead to clinical usefulness. However, much more work is necessary as fNIRS still has several limitations that prevents it from being a mature technique in clinical settings. To begin with, one problem concerns the interpretation of fNIRS results, as the relation between the measured hemodynamic change and the underlying neural activity is not fully understood. For example, evidence of a deoxygenation process can be detected during some seizures both in our study and in previous work.18,19,22,24 Whether the factors such as the type of seizures, their duration, their frequency, age, or other comorbid conditions are underlying the deoxygenation remains to be investigated. On the other hand, the HRF shape for IEDs was also reported to vary across patients and brain regions in previous BOLD-fMRI studies.72,73 Our recent work tried to estimate the HRF while accounting for the patient-specific variability and nonlinearity to some extent,47 but did not provide a standard solution for the analysis of IEDs with fNIRS. The other critical issue is the SNR of fNIRS signals. fNIRS data acquisition can be disturbed by movement artifacts especially during seizures. While several methods have been proposed to address this problem either with additional equipment22,74 or with signal denoising techniques,33,75,76 there is still no standard method for the correction of motion artifacts in fNIRS signals.77 Another factor that might lower the SNR of fNIRS signals is the unknown contribution from extracerebral tissue, as the averaged path of photon propagation is considered to follow a banana-shaped pattern from superficial layers (e.g., skin, scalp, skull, etc.) to cerebral cortex. The extracerebral contamination would further interfere with data interpretation by reducing fNIRS sensitivity to brain tissue, which gets worse if channels with shorter-than-typical separation (e.g., ) are employed in the SD layout.78,79 Recent studies concluded that signals recorded from very short SD separations (80) may be used to filter out the some of the physiological noise arising from superficial tissues.8182.–83 However, installing short channels together with long channels will probably lead to reduced head coverage, which is disadvantageous in epilepsy studies aimed at the localization of the focus region. Finally, it must be reminded that fNIRS can only reliably sample the superficial neocortex due to its limited penetration depth,6,19 a disadvantage shared by surface EEG.
While additional hardware/software development and clinical validation are necessary, results obtained so far indicate that fNIRS may become in time a useful complementary clinical tool. Current methods to visualize the hemodynamics associated with abnormal epileptiform discharges in patients with focal epilepsy include the application of ictal SPECT, PET, and combined EEG-fMRI. All the aforementioned techniques have their limitations. The localization yield of ictal SPECT is low (especially with extratemporal lobe epilepsies) as a single ictal perfusion brain image can rarely discriminate between the seizure onset zone and areas receiving propagated activity.84 Indeed, the ictal SPECT is a single snapshot in time that lacks temporal information, similarly for PET, which measures glucose metabolism in a single scan. EEG-fMRI studies have dealt quite successfully with IEDs, e.g., a recent study reported 88% sensitivity and 64% specificity in the focus localization of 21 patients,85 but this method has difficulties with seizures as they are unpredictable and rarely occur at the exact moment of scanning. The difficulty in recording seizures with fMRI is exacerbated by the fact that small patient movements during a seizure can have a severely negative effect on data quality.86 None of the above imaging tests can be used for continuous monitoring for several reasons including cost, availability, safety, and portability. Building on these unique qualities, potential fNIRS applications include the following: (a) rapid detection of the onset of a seizure to warn patients to protect themselves or nurses/physicians to intervene; (b) help in localizing the epileptic focus for eventual epilepsy surgery; (c) determination of the impact of seizures on brain oxygenation and cognition to allow more aggressive treatment (combination of antiepileptic drugs, earlier referral for surgery, oxygen therapy, etc.); and (d) lateralization of language function prior to epilepsy surgery. Other potential applications include the evaluation of critically ill patients in the intensive care unit: assessment of convulsive and nonconvulsive status epilepticus, appraisal of the impact of repeated seizures on tissue deoxygenation, continuous assessment of electrical and hemodynamic activity in acute stroke, examination of several ill-defined EEG patterns of unclear pathophysiology (periodic discharges, stimulus-induced rhythmic ictal discharges, etc.) but frequently seen in critically ill neurological patients, etc. Adding continuous NIRS to video-EEG telemetry has the potential to boost our capacity to monitor the brain in such situations, allowing early intervention before irreversible brain damage can ensue.
Work on fNIRS and epilepsy over the last two decades has provided cumulative promising results. To further validate the clinical utility of EEG-fNIRS in epilepsy research and treatment, the following challenges are proposed: a deeper understanding of hemodynamics associated with epileptic events; the development of a comfortable recording system with large spatial coverage combined with long term monitoring for days or even weeks, and portable use; and the improvements of analytical methods including appropriate modeling and artifact reduction. Randomized and blinded protocols with a larger patient number are also mandatory.
This work was supported by the Fonds de Recherche en Santé du Québec (FRSQ) Grant No. 14385, the Canadian Institutes of Health Research (CIHR), Institute of Circulatory and Respiratory Health (ICRH), and the Heart and Stroke Foundation of Canada (HSFC) Grant No. 62573, and the Savoy Foundation. The authors thank Parikshat Sirpal for suggested improvements to text quality.
Ke Peng is currently a PhD student in biomedical engineering at École Polytechnique de Montréal, Canada. He received his BSc degree in automation from Tsinghua University, China, in 2008, and his MScA degree in computer engineering from École Polytechnique de Montréal, Canada, in 2012.
Philippe Pouliot holds a PhD in physics from Rutgers University. He was a researcher in theoretical particle physics at the University of California at Santa Barbara, University of Texas at Austin, and Imperial College, London, UK. He is now a researcher in biomedical engineering at École Polytechnique de Montréal. His current research interests include neuroscience, epilepsy, aging, atherosclerosis, biophysical models, multimodal imaging (fMRI, fNIRS, OCT, 2-photon microscopy), and developing mathematical and software tools for data analysis.
Frédéric Lesage is a professor of electrical engineering at the École Polytechnique de Montréal and director of the Optical and Molecular Imaging Laboratory. He has pursued an interdisciplinary career in applied mathematics, physics, and imaging. His current research activities pertain to the development of novel imaging techniques for neuronal conditions, involving the analysis of optical signals impervious to physiological background noise, 3-D image reconstruction using multimodal instruments, time-domain optical parameter recovery, and multimodal imaging (fMRI-optical and EEG-optical).
Dang Khoa Nguyen completed his residency in neurology and a PhD in neuroscience at the Université de Montréal as well as a fellowship in epilepsy at the Yale University. He currently practices at Notre-Dame Hospital where he is the director of the Epilepsy Monitoring Unit. One of his research interests is the development of NIRS-EEG to monitor seizures in various clinical settings, including epilepsy monitoring units, acute stroke units, and intensive care units.