ORIGINAL RESEARCH article

Front. Neurol., 15 April 2025

Sec. Applied Neuroimaging

Volume 16 - 2025 | https://doi.org/10.3389/fneur.2025.1552774

Quantitative SISCOM assessment for epileptogenic zone localization: insights from a multicenter study comparing two software platforms in temporal lobe epilepsy

Carla Oliveira Young
Carla Oliveira Young1*Brunno Machado de CamposBrunno Machado de Campos2Edna Marina de SouzaEdna Marina de Souza3Sergio Querino BrunettoSergio Querino Brunetto3Maria Julia de Oliveira Santos GualbertoMaria Julia de Oliveira Santos Gualberto4Leonardo Alexandre-SantosLeonardo Alexandre-Santos4Marcos Geraldo Merichelo de OliveiraMarcos Geraldo Merichelo de Oliveira4Marina Koutsodontis Machado AlvimMarina Koutsodontis Machado Alvim2Fernando CendesFernando Cendes2Elba EtchebehereElba Etchebehere1Lauro Wichert-AnaLauro Wichert-Ana4Brbara Juarez Amorim
Bárbara Juarez Amorim1*
  • 1Division of Nuclear Medicine, Department of Radiology and Oncology, University of Campinas (UNICAMP), Campinas, Brazil
  • 2Division of Epilepsy, Department of Neurology, University of Campinas (UNICAMP), Campinas, Brazil
  • 3Centre of Biomedical Engineering, University of Campinas (UNICAMP), Campinas, Brazil
  • 4Nuclear Medicine and PET/CT Laboratory, Department of Medical Imaging, Haematology and Clinical Oncology, Ribeirão Preto Medical School, University of São Paulo, Ribeirão Preto, Brazil

Introduction: Pharmacoresistant epilepsy affects around one-third of individuals with epilepsy, requiring precise diagnosis, particularly in cases where surgical resection of the epileptogenic zone (EZ) is an option. Functional imaging techniques, such as ictal-interictal subtraction SPECT coregistered to MRI (SISCOM), have proven useful in pre-surgical evaluation by improving EZ localization accuracy. However, the widespread use of SISCOM is limited by the high costs and technical complexity of commercial software. Statistical Parametric Mapping (SPM) has been demonstrated to be a viable alternative for SISCOM analysis, displaying the potential for cost-effective EZ localization.

Materials and methods: In this retrospective study, we evaluated patients with pharmacoresistant temporal lobe epilepsy from two reference centers of epilepsy in Brazil, who underwent ictal and interictal SPECT imaging as part of their pre-surgical evaluation, achieving favorable outcomes (Engel I or II) after surgical resection. The EZ reference standard was determined according to anatomopathological findings and good clinical outcomes. SISCOM was performed using a semi-automated approach with Statistical Parametric Mapping (SPM) and a proprietary software – Analyze. Data from each method were compared to the EZ reference standard and classified as concordant, partially concordant, or discordant.

Results: We included 20 patients, 14 (70%) with left temporal lobe epilepsy and six (30%) with right temporal lobe epilepsy. Hippocampal sclerosis was the most common pathology (80%). Both SPM and Analyze were concordant with the EZ reference standard in 14 cases (70%), showing no difference in sensitivity between the methods. However, SPM generated smaller, more localized clusters, while Analyze produced larger clusters with broader spatial coverage. Concordance between the two methods was poor (Kappa = 0.0179), reflecting methodological differences.

Conclusion: This study evidences technical differences between SISCOM performed with SPM and Analyze, but with similar sensitivity (70%) for EZ localization. Further studies with larger sample sizes are required to confirm these findings. The data presented here suggest that SISCOM-SPM, due to its rapid and semi-automated workflow, may offer a practical and accessible alternative to proprietary software for epilepsy surgical planning.

1 Introduction

Pharmacoresistant epilepsy affects around one-third of individuals with epilepsy (1). This condition is defined by the failure to achieve seizure control, even after the use of at least two appropriately selected and well-tolerated antiepileptic drugs (1). For these patients, surgical resection of the epileptogenic zone (EZ) is the most effective therapeutic option to achieve seizure freedom (2). However, the precise localization of the EZ remains significantly challenging, particularly in patients with normal or inconclusive magnetic resonance imaging (MRI) (3). Thus, pharmacoresistant epilepsy requires considerable clinical attention, demanding advanced diagnostic and therapeutic strategies (1).

In this context, functional neuroimaging techniques, such as single-photon emission computed tomography (SPECT), have proven to be valuable tools in the pre-surgical evaluation of patients with pharmacoresistant epilepsy, complementing structural MRI findings (3, 4). SPECT, particularly when performed during a seizure (ictal SPECT), assists in identifying the EZ by highlighting increased cerebral perfusion. Generally, cerebral SPECT images are analyzed visually, presenting a diagnostic sensitivity for temporal lobe epilepsy of approximately 97% for ictal SPECT and 44% for interictal SPECT (5). However, visual analysis can often be challenging (6), and an accurate definition of the EZ is crucial (7).

The subtraction ictal-interictal SPECT coregistered to MRI (SISCOM) technique optimizes this process by allowing a direct comparison between ictal (hyperperfused) and interictal (hypoperfused) images, thereby increasing sensitivity and specificity in detecting the EZ (711). SISCOM images have higher spatial resolution than ictal or interictal images individually and can accurately localize the EZ, even in patients with focal cortical dysplasia and a normal MRI (12, 13). Studies such as Foiadelli et al. (4) have demonstrated the efficacy of SISCOM, particularly in pediatric epilepsy, assisting in surgical interventions even in cases with inconclusive MRI findings. Additionally, some studies reinforce the role of SISCOM in epilepsy surgery, highlighting improved postoperative outcomes (14).

Despite its clinical value, the widespread use of SISCOM is limited by the high costs and technical complexity of the commercial software required for its analysis (15). Few studies have developed accessible tools to replicate SISCOM’s functionalities at a lower cost, highlighting the need for different solutions (8, 10).

As an alternative, Statistical Parametric Mapping (SPM) is a free software well-documented in the literature for performing cerebral quantification (1618). Integrating SISCOM with Statistical Parametric Mapping (SPM) has improved image quantification and provided a cost-effective platform (4, 19).

The present authors have already used SPM for performing SISCOM in a previous scientific study, and this approach was designated as SISCOM-SPM (20). The aim of the present study was to apply SISCOM-SPM and compare it with another well-established proprietary software for localizing the EZ. We conducted a multicenter study across two major national reference centers in epilepsy.

2 Materials and methods

This is a retrospective, multicentre study involving patients from the University of Campinas Clinical Hospital and the Ribeirão Preto Medical School Clinical Hospital who had previously undergone ictal and interictal cerebral SPECT imaging as part of epilepsy follow-up. Both participating institutions have expertise in these quantitative analyses, including SISCOM-SPM and a proprietary software named Analyze.

This study was approved by the Research Ethics Committee of the Faculty of Medical Sciences of the University of Campinas (UNICAMP). The protocol was registered under the CAAE: 36549220.9.1001.5404.

2.1 Inclusion and exclusion criteria

Inclusion criteria for this study were as follows: patients with pharmacoresistant epilepsy, unifocal temporal epilepsy; those who previously undergone ictal and interictal cerebral SPECT studies as part of the routine investigation for the EZ; those who agreed to provide informed consent; submitted to surgical resection; seizure-free or almost seizure-free outcome after surgery (Engel I and II).

Exclusion criteria consisted of extratemporal or multifocal temporal epilepsy, not having both ictal and interictal cerebral SPECT studies available, and refusal to participate in the study or to provide informed consent.

As all patients underwent surgery for EZ resection, the reference standard for EZ localization was defined by the location of anatomopathological findings indicating areas consistent with epileptogenic lesions alongside favorable patient outcomes (Engel I or II). This is considered the gold standard for EZ localization (13).

2.2 Patients selection

We analyzed the medical records of patients who underwent cerebral SPECT for epileptogenic focus investigation between 2015 and 2019, who proceeded to surgical resection of the EZ and demonstrated a favorable clinical outcome, remaining seizure-free or with rare seizures (Engel I and II). Patients were selected from two national reference institutions, the University of Campinas and the Ribeirão Preto Medical School, in Brazil.

2.3 Patients preparation for interictal SPECT imaging

All patients remained in a low-light, quiet room for 15 min, with permanent intravenous access through a butterfly connected to a catheter with saline solution. While at rest, 1110 MBq (30 mCi) of 99mTc-ECD were injected. The patients rested for another 10 min before the SPECT acquisition. EEG monitoring was not performed during interictal SPECT acquisition, but no seizures were observed clinically.

2.4 Patients preparation for ictal SPECT imaging

All patients were monitored with a long-term video-EEG. The antiseizure medication was reduced in selected cases to increase the chance of epileptic seizures. Patients rested while continuous video-EEG was recorded. All patients remained with permanent intravenous access through a butterfly connected to a catheter with a saline solution. A syringe with the radiotracer was attached to the catheter and protected with a lead shield to ensure a fast injection. Upon seizure onset, around 1,110 MBq (30 mCi) of 99mTc-ECD were quickly injected. EEG and video recordings confirmed seizures. SPECT images were acquired 30–90 min after seizure cessation and patients’ symptoms. For the five patients from the Campinas cohort, radiotracer injection occurred between 13 and 24 s after seizure onset (mean: 16.6 s), consistent with literature recommendations for optimal EZ localization (13). Injection time data were not available for the Ribeirão Preto cohort due to retrospective data limitations.

2.5 Brain SPECT acquisition

SPECT was performed using a Symbia T2 SPECT/CT system (Siemens, Erlangen, Bayern, Germany) with a high resolution, low energy, two-head collimator. The SPECT images were acquired, with photopeak centered at 140 keV and 15% window, 128 × 128 matrix, 2.67 zoom (which could be variable), and 64 views for each head (37 s per view). Raw data were reconstructed with 3D OSEM (17 intersections, 16 subsets), 7.65 mm Gaussian filter and CT attenuation correction. Images were displayed in transaxial, coronal, and sagittal planes for visual analysis.

2.6 MRI acquisition

The MRI epilepsy protocol was conducted on a 3 Tesla Philips Achieva scanner and included: (1) Coronal images perpendicular to the hippocampus: T1-weighted “inversion recovery” (voxel of 0.64×0.64×3 mm3, gap of 0.3 mm, 45 slices, TR 3,550 ms, TE 15 ms, IR delay of 400 ms); (2) Axial images parallel to the hippocampus: DWI (voxel of 1.19×1.19×4 mm3, 28 slices, TR 3,474 ms, TE 70 ms, max b-factor of 1,000 s/mm2); (3) Additional 3D sequences aligned to the corpus callosum: 3D T1-weighted image (1 mm isotropic voxels, matrix of 240×240, 180 sagittal slices, TR of 7 ms, TE of 3.2 ms), 3D-FLAIR (voxel of 1.2×1.2×1.2 mm3, matrix of 180×180 and 300 slices, TR of 1,000 ms, TI of 2,400 ms), T2WI 3D (1 mm isotropic voxels, matrix of 250×250, 327 slices, TR of 2000 ms and TE of 364 ms).

2.7 SPECT visual analysis

Two experienced nuclear medicine physicians with over 20 years of experience in nuclear medicine brain images performed the analysis. These two nuclear medicine physicians evaluated the images in consensus, looking for focal areas of hyper or hypoperfusion and comparing both cerebral hemispheres. They also compared brain perfusion with the cerebellar perfusion. Nuclear medicine physicians were aware of clinical and electroencephalographic findings. The EZ was defined as a focal area of hyperperfusion in the ictal SPECT images, which became hypoperfused or normally perfused in interictal SPECT images. When there was more than one hyperperfused area, we correlated the findings with the time injection and ictal semiology to define the EZ. The other areas were considered as propagation areas.

2.8 Subtraction ictal-interictal SPECT co-registered to MRI with SPM (SISCOM-SPM)

A trained physician performed the ictal and interictal subtractions using the Statistical Parametric Mapping (SPM) software, version 12 (SPM12) (Wellcome Department of Imaging Neuroscience, University College London, UK)1 in Matlab (MathWorks, Natick, MA, USA) (21).

After acquiring and reconstructing interictal and ictal images, they were converted from DICOM to Analyze format using MRIcro software. Both ictal and interictal images in the Analyze format were loaded into SPM12, which runs in Matlab software.

The images were subsequently realigned using the anterior commissure as a reference. In the registration step, ictal and interictal images were registered using routines based on Mutual Information maximization algorithms. Both images were taken as target and source in two different registrations. The registered images were realigned, and the mean intensity value calculated, normalizing the uptake levels toward the whole brain. The registered images’ positional correspondence and mean were checked in SPM12. After establishing that all areas were correspondent, the subtraction was done. To obtain the cerebral perfusion differences, the transformed and normalized interictal SPECT image was subtracted voxel by voxel from the ictal SPECT image using the LCN12 subtract routine (21). The difference image was transformed into a z-score using the mean and standard deviation (SD) of the differences in all brain voxels (22). All voxels exceeding two z-scores were considered significant. Usually, a hyperperfused region appears as a cluster, which is a set of voxels.

The subtraction process with SPM took approximately 15 min.

2.9 Subtraction ictal-interictal SPECT co-registered to MRI with Analyze (SISCOM-Analyze)

SISCOM was executed using the ANALYZE © 10.0 Software (AnalyzeDirect, Inc., Overland Park, Kansas, USA), applying a methodology adapted from prior research (23). Subtraction was carried out by subtracting the interictal SPECT scan from the ictal SPECT, transforming the difference in the signal into z-score maps, calculated from the mean and standard deviation across all brain voxels, producing images that reflected changes in regional cerebral blood flow (rCBF). This quantitative data matrix was then fused with the patient’s MRI to integrate functional and anatomical insights within a single image. After functional overlay, only regions with significant changes in rCBF—defined as greater than two standard deviations above or below the mean—were displayed.

The subtraction process with Analyze typically takes approximately 60 min when performed by an expert, whereas it can exceed 2 h for less experienced individuals. While this reflects general practice, future studies are needed to accurately quantify the time required under different levels of expertise.

2.10 EZ analysis with SISCOM-SPM and SISCOM-Analyze

Specific parameters were adopted to investigate the localization of the epileptogenic zone, including the voxel with the highest intensity and the largest clusters acquired through SISCOM performed using both methods: (a) EZ would be the highest intensity cluster (Z score); (b) EZ would be the largest cluster (number of voxels); (c) EZ would be the second largest cluster (number of voxels); and (d) EZ would be the third-largest cluster (number of voxels).

2.11 SISCOM analyses vs. EZ reference standard

The location of the EZ identified by both methods, SISCOM-SPM and SISCOM-Analyze, was compared to the reference standard and categorized as follows:

• Concordant: when there was an overlap between the hyperperfused area and the EZ (Figure 1);

• Partially concordant: when the hyperperfused area was adjacent or very close to the EZ;

• Discordant: when the hyperperfused area was distant from the EZ (Figure 2).

Figure 1
www.frontiersin.org

Figure 1. Patient 4: 50-year-old female with left hippocampal sclerosis. Both methods (SPM and Analyze) demonstrated concordance with the gold standard, identifying the EZ in the left temporal lobe. (A) Axial slices from the SISCOM analysis using SPM, showing the largest cluster and the voxel with the highest intensity located in the left temporal lobe, concordant with the EZ. (B) Axial slices from the Analyze software, also identifying the largest cluster and the voxel with the highest intensity in the left temporal lobe, both concordant with the EZ.

Figure 2
www.frontiersin.org

Figure 2. Patient 17: 45-year-old female with focal cortical dysplasia (FCD) IB in the right temporal lobe. Neither method demonstrated concordance with the gold standard. (A) Axial slices from the SISCOM analysis using SPM, showing multiple clusters discordant with the EZ, which was located in the right temporal lobe. (B) Axial slices from the Analyze software, also displaying multiple clusters inconsistent with the EZ in the right temporal lobe.

2.12 Statistical analysis

Descriptive analysis was performed for categorical variables using frequency tables and quantitative variables using measures of central tendency and dispersion. Agreement between methods (SISCOM-Analyze and SISCOM-SPM with the EZ reference standard) was evaluated using the Kappa coefficient, with values classified as poor (<0.40), moderate (0.40–0.75), or excellent (>0.75). Mann–Whitney and Fisher’s exact tests were applied to compare groups, and Spearman’s rank correlation coefficient was used to calculate the numerical variables. The level of significance adopted was 5%.

In some analyses, cases classified as discordant and partially concordant were grouped into a non-concordant category to facilitate statistical evaluation and simplify interpretation.

3 Results

We analyzed the medical records of 97 individuals. Among them, we included five patients from the University of Campinas and 15 patients from Ribeirão Preto School, who had ictal and interictal SPECT scans, and proceeded to surgical resection of the EZ and had a favorable postoperative outcome (Engel I and II). Therefore, we included 20 patients from both institutions in total. Their mean age was 42.5 ± 12.2 years (median: 44.5 years; range: 10.0–58.0 years). Eleven patients (55.0%) were female and nine (45.0%) were male. All patients had temporal lobe epilepsy, and based on the reference standard, six (30.0%) had right temporal lobe epilepsy, and 14 (70.0%) had left temporal lobe epilepsy.

According to anatomopathological findings, 16 (80.0%) had hippocampal sclerosis, one (5.0%) focal cortical dysplasia and three (15.0%) other pathologies (astrogliosis, grade 1 ganglioglioma and fibrillar gliosis). Table 1 summarizes the demographic data, gold standard localization, anatomopathological findings, and post-surgical outcomes (Engel classification) of the 20 patients included in the study.

Table 1
www.frontiersin.org

Table 1. Patients’ clinical and anatomopathological characteristics.

3.1 SISCOM results

The mean total cluster volume was 29,778.1 ± 17,920.2 mm3 using SPM and 17,386.8 ± 3,969.9 voxels using Analyze.

The mean number of clusters in hippocampal sclerosis patients was 11.8 (SD = 12.3) and 34.0 in the other pathologies (SD = 20.7), and this difference had statistical significance (p = 0.0252).

Table 2 presents the highest voxel intensity, total number of clusters, and the volumes of the three largest clusters (in mm3) identified by the SPM and Analyze methods.

Table 2
www.frontiersin.org

Table 2. Quantitative SISCOM results.

Highest intensity cluster – The average intensity was 3.7 ± 0.7 in SPM and 252.6 ± 2.8 in Analyze. The mean number of clusters was 16.2 ± 16.4 in SPM and 58.8 ± 31.0 in Analyze.

Largest cluster – The mean volume of the largest cluster was 21,855.6 ± 16,798.5 mm3 in SPM and 26,121.8 ± 16,032.4 mm3 in Analyze, with mean intensities of 2.7 ± 0.5 and 49.7 ± 9.2, respectively.

Second largest cluster – The mean volumes were 4,422.2 ± 4,471.7 mm3 in SPM and 8,706.3 ± 4,241.7 mm3 in Analyze, with mean intensities of 2.5 ± 0.4 and 39.1 ± 11.0, respectively.

Third largest cluster – The mean volumes were 1,741.2 ± 1,815.0 mm3 in SPM and 4,335.2 ± 2,594.9 mm3 in Analyze, with mean intensities of 2.4 ± 0.4 and 42.0 ± 10.2, respectively.

3.2 Patients based analysis

In the 20 patients, SPM identified 14 patients as concordant with EZ, seven by highest intensity and 14 by highest cluster volume with a sensitivity of 70%. The largest cluster also identified all seven patients with the highest intensity (Table 3). Analyze identified EZ in 14 patients, either by the highest intensity or largest cluster (Table 3), with a 70% sensitivity.

Table 3
www.frontiersin.org

Table 3. Concordance of SISCOM analyses with the epileptogenic zone (EZ) reference standard.

4 Discussion

The present study analyzed the free software Statistical Parametric Mapping (SPM) and the proprietary software Analyze for EZ localization in patients with drug-resistant temporal lobe epilepsy from two national reference centers for epilepsy studies. We demonstrated that both software tools used to perform SISCOM have similar capabilities to localize the EZ, with a sensitivity of 70%.

The primary institution has experience with SPM (24, 25) and previously published results from a study featuring 23 patients with pharmacoresistant epilepsy, demonstrating that SISCOM-SPM improved EZ localization by 21.8% compared to visual analysis alone (20). Meanwhile, the participating institution has scientific experience with Analyze software, which has been applied by the team in epilepsy studies (3, 23, 26) as well as for detecting reperfusion in systemic lupus erythematosus (27).

We noticed a low concordance between SISCOM performed with SPM and Analyze regarding the highest intensity voxel (Kappa: 0.0179). Matsuda et al. (10) also recognized significant variability in the performance of different software tools in clinical and research involving multicenter studies. It can be attributed to differences in how each algorithm processes the data, potentially impacting its sensitivity and spatial resolution. In our study, the poor concordance (Kappa = 0.0179) likely reflects intrinsic differences between the software algorithms, including variations in image preprocessing, intensity normalization, z-score thresholding, and spatial smoothing. SPM tends to identify smaller, focal hyperperfused clusters, while Analyze detects broader activation patterns that may also include propagation zones. Future studies incorporating voxel-wise overlap analysis and algorithmic benchmarking could provide deeper insights into these discrepancies and support standardization across platforms.

A significant difference was observed between the methods compared to applicability in clinical practice. Analyze depends on the technician who performs it, with processing steps lasting over 2 h, depending on the professional’s experience. In contrast, the method developed by the researchers of this study using SPM offers a practically automated flow, with processing completed in a few minutes, making it more accessible and efficient in clinical routines. From a clinical perspective, the streamlined workflow of SISCOM-SPM, requiring approximately 15 min for image processing, offers a considerable advantage over the more time-consuming and operator-dependent Analyze software. This difference has practical implications for routine epilepsy centers, where timely decision-making and reproducibility are essential for surgical planning.

In the present study, SISCOM performed with SPM revealed clusters with smaller volumes (mean volume of the largest cluster: 21,855.6 mm3 for SPM versus 26,121.8 mm3 for Analyze) and lower absolute values of maximum voxel intensity (3.7 ± 0.7 for SPM compared to 252.6 ± 2.8 for Analyze). These differences in cluster characteristics likely contributed to the discordance observed between the methods, as SPM tends to delineate more compact regions while Analyze identifies larger regions. Despite this, SPM demonstrated superior anatomical discrimination, particularly in cases of hippocampal sclerosis, further supporting its clinical applicability in patients with mesial temporal lobe epilepsy, where the epileptogenic zone occupies a smaller area.

Visual analyses of the images obtained from each method revealed that SISCOM performed using Analyze generates clusters with larger volumes and a higher number of clusters overall (an average of 14 clusters in Analyze versus 9 clusters in SPM). However, it was observed that these clusters often also include areas unrelated to the EZ, corresponding to regions of seizure propagation. This emphasizes the importance of establishing appropriate z-score thresholds, as discussed by De Coster et al. (22), who noted that lower thresholds increase the method’s sensitivity but may also include activation areas not directly related to the EZ.

There were some limitations in the present study. First of all, among the challenges of conducting a multicenter study is the difficulty in easily accessing some data. In particular, the injection time for patients from the Ribeirão Preto cohort was unavailable, limiting our ability to fully assess timing-related variability in SISCOM performance. As discussed by O’Brien et al., it is known that shorter intervals between the onset of the seizure and the administration of the radiotracer increase the accuracy in localizing the EZ. We found an unideal sensitivity in EZ localization with both tools, achieving 70%. This is slightly less than the literature which is around 92–93% in temporal lobe epilepsy (5, 13) and between 76–85% for refractory epilepsy in general (7, 17). This lower sensitivity may, in part, reflect longer delays between seizure onset and radiotracer injection in some cases. An additional limitation is the absence of EEG monitoring during interictal SPECT acquisitions, which may affect the ability to detect electrographic seizures that do not manifest clinically, especially in temporal lobe epilepsy. Furthermore, visual interpretation of the SPECT images was performed with access to clinical and EEG data. While this reflects standard clinical practice, it may have introduced bias in visual assessment.

Moreover, the small number of participants in the study limits definite conclusions. Only a few patients with postoperative Engel I-II outcomes from both institutions had the ictal and interictal SPECT available, resulting in a small number of participants in this study, consequently reducing the robustness of the statistical analyses. The numerical distribution from the statistical analysis indicates similar results for both methods, with comparable levels of agreement and disagreement with the reference standard. Consequently, the performance of the methods is similar; however, it is not statistically possible to determine which method is more accurate. Future prospective multicenter studies with standardized acquisition protocols, harmonized SISCOM analysis pipelines, and voxel-level comparison metrics are warranted to improve reproducibility, strengthen clinical reliability, and support broader adoption in surgical decision-making.

5 Conclusion

This study detected some technical differences between SISCOM using SPM and Analyze for the localization of the EZ in temporal lobe epilepsy. While Analyze presents larger and more numerous clusters, potentially increasing the possibility of encompassing the correct epileptogenic focus, it also includes seizure propagation regions and is lengthier. In contrast, SPM provides a refined anatomical delineation, particularly in cases with focal alterations such as hippocampal sclerosis, being faster and more feasible in clinical practice. However, despite these technical differences, both tools could localize EZ with similar sensitivity, suggesting that either can be effectively used to evaluate epilepsy. In this context, SISCOM-SPM offers the advantage of a semi-automated and examiner-independent workflow, reducing not only processing time but also improving reproducibility and usability in clinical settings. These qualities make it a particularly suitable option for pre-surgical evaluation in epilepsy centers, especially those seeking practical and accessible imaging solutions.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement

The studies involving humans were approved by Comitê de Ética em Pesquisa da Universidade Estadual de Campinas (CEP-Unicamp), affiliated with State University of Campinas (Unicamp), Brazil. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

CO: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. BC: Conceptualization, Data curation, Formal analysis, Methodology, Software, Validation, Visualization, Writing – review & editing. ES: Formal analysis, Writing – review & editing, Investigation, Methodology. SB: Methodology, Writing – review & editing. MG: Data curation, Investigation, Writing – review & editing. LA-S: Conceptualization, Data curation, Methodology, Software, Validation, Visualization, Writing – review & editing. MO: Data curation, Investigation, Software, Writing – review & editing. MA: Writing – review & editing. FC: Conceptualization, Funding acquisition, Methodology, Resources, Supervision, Validation, Visualization, Writing – review & editing. EE: Resources, Writing – review & editing. LW-A: Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – review & editing. BA: Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by the Brazilian Federal Agency for Support and Evaluation of Graduate Education (CAPES) – Finance Code 88882.434897/2019-01.

Acknowledgments

The authors thank Patrick Dupont and Karolien Goffin from the University of Leuven (Belgium) and Eduardo Marçal for their support and help.

Conflict of interest

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.

Generative AI statement

The authors declare that Gen AI was used in the creation of this manuscript. Artificial intelligence assisted in the article by providing specific suggestions for synonyms and technical terms, without any involvement in purely intellectual content or text development.

Publisher’s note

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.

Footnotes

References

1. Kalilani, L, Sun, X, Pelgrims, B, Noack-Rink, M, and Villanueva, V. The epidemiology of drug-resistant epilepsy: a systematic review and meta-analysis. Epilepsia. (2018) 59:2179–93. doi: 10.1111/epi.14596

PubMed Abstract | Crossref Full Text | Google Scholar

2. Sukprakun, C, and Tepmongkol, S. Nuclear imaging for localization and surgical outcome prediction in epilepsy: a review of latest discoveries and future perspectives. Front Neurol. (2022) 13:1083775. doi: 10.3389/fneur.2022.1083775

PubMed Abstract | Crossref Full Text | Google Scholar

3. Bianchin, MM, Wichert-Ana, L, Velasco, TR, Martins, AP, and Sakamoto, AC. Imaging epilepsy with SISCOM. Nat Rev Neurol. (2011) 7:240–2. doi: 10.1038/nrneurol.2010.131-c1

PubMed Abstract | Crossref Full Text | Google Scholar

4. Foiadelli, T, Lagae, L, Goffin, K, Theys, T, de Amici, M, Sacchi, L, et al. Subtraction ictal SPECT coregistered to MRI (SISCOM) as a guide in localizing childhood epilepsy. Epilepsia Open. (2020) 5:61–72. doi: 10.1002/epi4.12373

PubMed Abstract | Crossref Full Text | Google Scholar

5. Devous, MD Sr, Thisted, RA, Morgan, GF, Leroy, RF, and Rowe, CC. SPECT brain imaging in epilepsy: a meta-analysis. J Nucl Med. (1998) 39:285–93.

PubMed Abstract | Google Scholar

6. Kumar, A, and Chugani, HT. The role of radionuclide imaging in epilepsy, part 1: sporadic temporal and Extratemporal lobe epilepsy. J Nucl Med Technol. (2017) 45:14–21. doi: 10.2967/jnumed.112.114397

PubMed Abstract | Crossref Full Text | Google Scholar

7. Chen, T, and Guo, L. The role of SISCOM in preoperative evaluation for patients with epilepsy surgery: a meta-analysis. Seizure. (2016) 41:43–50. doi: 10.1016/j.seizure.2016.06.024

PubMed Abstract | Crossref Full Text | Google Scholar

8. Moreau, JT, Saint-Martin, C, Baillet, S, and Dudley, RWR. MNI SISCOM: an open-source tool for computing subtraction ictal single-photon emission CT Coregistered to MRI. J Digit Imaging. (2021) 34:357–61. doi: 10.1007/s10278-021-00422-9

PubMed Abstract | Crossref Full Text | Google Scholar

9. Suárez-Piñera, M, Mestre-Fusco, A, Ley, M, González, S, Medrano, S, Principe, A, et al. Perfusion SPECT, SISCOM and PET (18)F-FDG in the assessment of drug-refractory epilepsy patients candidates for epilepsy surgery. Rev Esp Med Nucl Imagen Mol. (2015) 34:350–7. doi: 10.1016/j.remn.2015.05.002

PubMed Abstract | Crossref Full Text | Google Scholar

10. Matsuda, H, Matsuda, K, Nakamura, F, Kameyama, S, Masuda, H, Otsuki, T, et al. Contribution of subtraction ictal SPECT coregistered to MRI to epilepsy surgery: a multicenter study. Ann Nucl Med. (2009) 23:283–91. doi: 10.1007/s12149-009-0236-6

PubMed Abstract | Crossref Full Text | Google Scholar

11. O'brien, TJ, So, EL, Mullan, BP, Cascino, GD, Hauser, MF, Brinkmann, BH, et al. Subtraction peri-ictal SPECT is predictive of extratemporal epilepsy surgery outcome. Neurology. (2000) 55:1668–77. doi: 10.1212/WNL.55.11.1668

PubMed Abstract | Crossref Full Text | Google Scholar

12. Stamoulis, C, Verma, N, Kaulas, H, Halford, JJ, Duffy, FH, Pearl, PL, et al. The promise of subtraction ictal SPECT co-registered to MRI for improved seizure localization in pediatric epilepsies: affecting factors and relationship to the surgical outcome. Epilepsy Res. (2017) 129:59–66. doi: 10.1016/j.eplepsyres.2016.11.020

PubMed Abstract | Crossref Full Text | Google Scholar

13. O'Brien, TJ, So, EL, Mullan, BP, Hauser, MF, Brinkmann, BH, Bohnen, NI, et al. Subtraction ictal SPECT co-registered to MRI improves clinical usefulness of SPECT in localizing the surgical seizure focus. Neurology. (1998) 50:445–54. doi: 10.1212/WNL.50.2.445

PubMed Abstract | Crossref Full Text | Google Scholar

14. Delev, D, Quesada, CM, Grote, A, Boström, JP, Elger, C, Vatter, H, et al. A multimodal concept for invasive diagnostics and surgery based on neuronavigated voxel-based morphometric MRI postprocessing data in previously nonlesional epilepsy. J Neurosurg. (2018) 128:1178–86. doi: 10.3171/2016.12.JNS161676

PubMed Abstract | Crossref Full Text | Google Scholar

15. Fisher, RS, Acevedo, C, Arzimanoglou, A, Bogacz, A, Cross, JH, Elger, CE, et al. ILAE official report: a practical clinical definition of epilepsy. Epilepsia. (2014) 55:475–82. doi: 10.1111/epi.12550

PubMed Abstract | Crossref Full Text | Google Scholar

16. Tousseyn, S, Dupont, P, Goffin, K, Sunaert, S, and Paesschen, WV. Correspondence between large-scale ictal and interictal epileptic networks revealed by single photon emission computed tomography (SPECT) and electroencephalography (EEG)-functional magnetic resonance imaging (fMRI). Epilepsia. (2015) 56:382–92. doi: 10.1111/epi.12910

PubMed Abstract | Crossref Full Text | Google Scholar

17. Sulc, V, Stykel, S, Hanson, DP, Brinkmann, BH, Jones, DT, Holmes, DR, et al. Statistical SPECT processing in MRI-negative epilepsy surgery. Neurology. (2014) 82:932–9. doi: 10.1212/WNL.0000000000000209

PubMed Abstract | Crossref Full Text | Google Scholar

18. Knowlton, RC, Lawn, ND, Mountz, JM, Buddhiwardhan, O, Miller, S, Burneo, JG, et al. Ictal single-photon emission computed tomography imaging in extra temporal lobe epilepsy using statistical parametric mapping. J Neuroimaging. (2004) 14:324–30. doi: 10.1177/1051228404268570

PubMed Abstract | Crossref Full Text | Google Scholar

19. Penny, WD, Friston, KJ, Ashburner, JT, Kiebel, SJ, and Nichols, TE. Statistical parametric mapping: the analysis of functional brain images. London: Elsevier (2011).

Google Scholar

20. Oliveira Young, C, Etchbehere, ECSC, Souza, EM, Brunetto, SQ, Santos, AO, Lima, MCL, et al. Clinical usefulness of SISCOM-SPM compared to visual analysis to locate the epileptogenic zone. Front Neurol. (2020) 11:467. doi: 10.3389/fneur.2020.00467

PubMed Abstract | Crossref Full Text | Google Scholar

21. Dupont, P. Multi-modal imaging in epilepsy: SPECT and PET. JBR-BTR. (2008) 91:258–61.

Google Scholar

22. De Coster, L, Van Laere, K, Cleeren, E, Baete, K, Dupont, P, Van Paesschen, W, et al. On the optimal z-score threshold for SISCOM analysis to localize the ictal onset zone. EJNMMI. (2018) 8:34. doi: 10.1186/s13550-018-0381-9

PubMed Abstract | Crossref Full Text | Google Scholar

23. Wichert-Ana, L, Ferruzzi, EH, Alexandre, V Jr, Velasco, TR, Bianchin, MM, Araújo, WM, et al. Sphenoid sinus bleeding during generalized seizure: a rare SISCOM (subtraction ictal SPECT coregistered with MRI) finding mimicking skull base tumor. Clin Nucl Med. (2007) 32:45–6. doi: 10.1097/01.rlu.0000249761.32486.5e

PubMed Abstract | Crossref Full Text | Google Scholar

24. Amorim, BJ, Etchebehere, EC, Camargo, EE, Rio, PA, Bonilha, L, Rorden, C, et al. Statistical voxel-wise analysis of ictal SPECT reveals pattern of abnormal perfusion in patients with temporal lobe epilepsy. Arq Neuropsiquiatr. (2005) 63:977–83. doi: 10.1590/s0004-282x2005000600014

PubMed Abstract | Crossref Full Text | Google Scholar

25. Amorim, BJ, Ramos, CD, dos Santos, AO, da Cunha Lopes de Lima, M, Min, LL, Camargo, EE, et al. Brain SPECT in mesial temporal lobe epilepsy: comparison between visual analysis and SPM. Arq Neuropsiquiatr. (2010) 68:153–60. doi: 10.1590/s0004-282x2010000200001

PubMed Abstract | Crossref Full Text | Google Scholar

26. Wichert-Ana, L, de Azevedo-Marques, PM, Oliveira, LF, Fernandes, RM, Velasco, TR, Santos, AC, et al. Ictal technetium-99 m ethyl cysteinate dimer single-photon emission tomographic findings in epileptic patients with polymicrogyria syndromes: a subtraction of ictal-interictal SPECT coregistered to MRI study. Eur J Nucl Med Mol Imaging. (2008) 35:1159–70. doi: 10.1007/s00259-007-0655-3

PubMed Abstract | Crossref Full Text | Google Scholar

27. Trevisan, AC, Louzada, P Jr, Pitella, FA, Kato, M, Santos, AC, and Wichert-Ana, L. The reperfusion of ischemic vasculitis in a patient with neuropsychiatric lupus: a serial single-photon emission computed tomography study. Indian J Nucl Med. (2019) 34:167–8. doi: 10.4103/ijnm.IJNM_135_18

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: SISCOM, SPECT, quantitative analysis, temporal lobe epilepsy, pharmacoresistant epilepsy, epileptogenic zone

Citation: Oliveira Young C, de Campos BM, de Souza EM, Brunetto SQ, Gualberto MJdOS, Alexandre-Santos L, de Oliveira MGM, Alvim MKM, Cendes F, Etchebehere E, Wichert-Ana L and Amorim BJ (2025) Quantitative SISCOM assessment for epileptogenic zone localization: insights from a multicenter study comparing two software platforms in temporal lobe epilepsy. Front. Neurol. 16:1552774. doi: 10.3389/fneur.2025.1552774

Received: 29 December 2024; Accepted: 31 March 2025;
Published: 15 April 2025.

Edited by:

Marino M. Bianchin, Federal University of Rio Grande do Sul, Brazil

Reviewed by:

Carolina Torres, Clinical Hospital of Porto Alegre, Brazil
Dikhra Khan, All India Institute of Medical Sciences, India

Copyright © 2025 Oliveira Young, de Campos, de Souza, Brunetto, Gualberto, Alexandre-Santos, de Oliveira, Alvim, Cendes, Etchebehere, Wichert-Ana and Amorim. 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) and the copyright owner(s) 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: Carla Oliveira Young, Y2FybGEuc29saXlvdW5nQGdtYWlsLmNvbQ==; Bárbara Juarez Amorim, YmFyYmp1QHVuaWNhbXAuYnI=

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.