- 1Department of Neurosurgery, The Translational Research Institute for Neurological Disorders of Wannan Medical College, The First Affiliated Hospital of Wannan Medical College, Yijishan Hospital of Wannan Medical College, Wuhu, China
- 2Department of Radiology, The First Affiliated Hospital of Wannan Medical College, Yijishan Hospital of Wannan Medical College, Wuhu, China
Background: The tumor invasion of the frontal lobe induces changes in the executive control network (ECN). It remains unclear whether epileptic seizures in frontal glioma patients exacerbate the structural and functional alterations within the ECN, and whether these changes can be used to identify glioma-related seizures at an early stage. This study aimed to investigate the altered structural and functional patterns of ECN in frontal gliomas without epilepsy (non-FGep) and frontal gliomas with epilepsy (FGep) and to evaluate whether the patterns can accurately distinguish glioma-related epilepsy.
Methods: We measured gray matter (GM) volume, regional homogeneity (ReHo), and functional connectivity (FC) within the ECN to identify the structural and functional changes in 50 patients with frontal gliomas (29 non-FGep and 21 FGep) and 39 healthy controls (CN). We assessed the relationships between the structural and functional changes and cognitive function using partial correlation analysis. Finally, we applied a pattern classification approach to test whether structural and functional abnormalities within the ECN can distinguish non-FGep and FGep from CN subjects.
Results: Within the ECN, non-FGep and FGep showed increased local structure (GM) and function (ReHo), and decreased FC between brain regions compared to CN. Also, non-FGep and FGep showed differential patterns of structural and functional abnormalities within the ECN, and these abnormalities are more severe in FGep than in non-FGep. Lastly, FC between the right superior frontal gyrus and right dorsolateral prefrontal cortex was positively correlated with episodic memory scores in non-FGep and FGep. In particular, the support vector machine (SVM) classifier based on structural and functional abnormalities within ECN could accurately distinguish non-FGep and FGep from CN, and FGep from non-FGep on an individual basis with very high accuracy, area under the curve (AUC), sensitivity, and specificity.
Conclusion: Tumor invasion of the frontal lobe induces local structural and functional reorganization within the ECN, exacerbated by the accompanying epileptic seizures. The ECN abnormalities can accurately distinguish the presence or absence of epileptic seizures in frontal glioma patients. These findings suggest that differential ECN patterns can assist in the early identification and intervention of epileptic seizures in frontal glioma patients.
Introduction
Glioma is the most common malignant tumor of the central nervous system (Ostrom et al., 2019), and the incidence of epilepsy in glioma can be as high as 40–90% (Pallud et al., 2014; Fisher et al., 2017; Akeret et al., 2020; You et al., 2020). Tumor-related epilepsy profoundly impacts patients’ prognosis and quality of life (Maschio and Dinapoli, 2012). However, how tumors induce epilepsy is still unclear. This poses a significant challenge for antiepileptic therapy in patients with newly diagnosed tumors (Siomin et al., 2005; Tremont-Lukats et al., 2008). Tumor invasion of the frontal lobe induces changes in the executive control network (ECN) (Liu Y. et al., 2020). Similarly, the onset of epilepsy also leads to functional network alterations (Cao et al., 2014; Fang et al., 2022). Therefore, an in-depth understanding of the structural and functional patterns of the ECN in glioma-associated epilepsy is required for the early differentiation and diagnosis of glioma seizures and the development of individualized intervention strategies.
Executive control network is a frontoparietal circuit involved in executive functions (Crone et al., 2006; Sakai and Passingham, 2006; Dosenbach et al., 2007; Seeley et al., 2007, 2009; Brier et al., 2012; Chen et al., 2016b). Tumor invasion of the frontal lobe (located in ECN) affects the brain’s executive functions (Seeley et al., 2009; Smith et al., 2009; Park et al., 2015; Zhang et al., 2016; Liu et al., 2019). Recent studies have consistently indicated that when the tumor invades the frontal lobe, the frontal glioma shows structural [gray matter (GM) reorganization and reduction in cortical thickness] (Kinno et al., 2020; Liu Y. et al., 2020) and functional [functional reorganization and decreased functional connectivity (FC)] abnormalities within the ECN (Liu et al., 2019; Liu Y. et al., 2020; Tordjman et al., 2021). Furthermore, idiopathic epilepsy reduces FC in networks (such as the salience network, default network, and ECN) (Cao et al., 2014; Li et al., 2017). In particular, some studies indicate that tumor-related epilepsy induces topological changes in white matter (WM) networks (Zhou et al., 2022) and functional networks (Fang et al., 2022). Based on these findings, we hypothesize that tumor invasion of the frontal lobe, accompanied by epileptic seizures, will lead to disruption and reorganization of the ECN. This damage will be aggravated compared to glioma patients without epilepsy.
In glioma surgery, after tumor control, treatment of glioma-related epilepsy is the second primary purpose (You et al., 2020). When combined with preoperative individualized and standardized antiepileptic therapies, glioma surgery can be highly effective in epilepsy control (You et al., 2020). However, identifying epileptic seizures early in glioma and administering individualized anti-seizure medicine is a significant challenge in clinical practice. As the clinical importance of radiomics and artificial intelligence (AI) has increased in recent years, deep learning and machine learning (ML) combined with multimodal imaging features have been widely used to predict glioma grades (Vamvakas et al., 2019; Abdelaziz Ismael et al., 2020), survival (Sanghani et al., 2018; Peeken et al., 2019; Choi et al., 2020), and molecular phenotype (Tan et al., 2019; Bangalore Yogananda et al., 2020). Furthermore, AI-based ML methods have been applied to various areas of epilepsy research, such as structural connectome and diffusion tensor imaging features (Munsell et al., 2015; Gleichgerrcht et al., 2018; Abbasi and Goldenholz, 2019). With this knowledge, we theorize that the combination of ML and multimodal neuroimaging data may be helpful in the early identification of glioma seizures. In particular, the structural and functional abnormalities of the ECN may be used for the early identification of frontal glioma-associated epileptic seizures.
This study aimed to investigate the altered structural and functional patterns of ECN in non-FGep and FGep subjects and to evaluate whether these patterns can accurately distinguish glioma-associated epileptic seizures. Based above-mentions studies (Liu et al., 2019; Liu Y. et al., 2020; Tordjman et al., 2021), we hypothesized that non-FGep and FGep would show structural (i.e., GM) and functional [i.e., regional homogeneity (ReHo) and FC] abnormalities of ECN and display differential damage patterns. Furthermore, we hypothesize that the altered patterns can accurately distinguish tumor invasion of the frontal lobe with or without an epileptic seizure. To test the hypothesis, we identified structural and functional changes within the ECN in non-FGep and FGep patients. Next, we applied a pattern classification approach to test how well structural and functional abnormalities within ECN could distinguish non-FGep and FGep from CN subjects.
Materials and Methods
Subjects
We recruited 50 patients with frontal gliomas from inpatients at the Department of Neurosurgery, Yijishan Hospital of Wannan Medical College (China). Patients were divided into 29 non-FGep (no clear history of seizures) and 21 FGep (clear history of seizures). The diagnosis of glioma-related epilepsy was based on clinical markers and electroencephalography findings. The CN were recruited through advertisements and matched to patients with frontal gliomas according to age, gender, and education (Table 1). The subjects signed written informed consent forms approved by the Human Participants Ethics Committee of Yijishan Hospital of Wannan Medical College (China).
The inclusion criteria for the patient groups were as follows: (1) tumor pathology was confirmed as primary glioma by surgery, (2) the tumor invaded the frontal lobe, (3) the extension of the tumor had not reached the central sulcus, (4) glioma had unilateral tumor invasion, (5) glioma had no brain injury, and (6) no history of biopsy, radiotherapy, or chemotherapy. The exclusion criteria were: (1) multiple lesion foci, (2) history of substance abuse, (3) magnetic resonance imaging (MRI) contraindications, and (4) gliomas involving bilateral prefrontal lobes.
Neuropsychological Assessments
A neuropsychological battery was utilized to assess episodic memory (auditory verbal learning test- 20-min delay recall, AVLT-20-DR), information processing speed (trail making test-A, TMT), visuospatial function (clock drawing test, CDT), and executive function (digit span test, DST). All subjects underwent a standardized clinical interview and comprehensive neuropsychological assessments that were evaluated by two neuropsychologists. These are senior doctors with more than 15 years of experience. Now we have added this information in section “Materials and Methods.” Specific characteristics for all subjects are provided in Table 1.
Magnetic Resonance Imaging Data Acquisition
The MRI data from 2019 to 2021 were scanned preoperatively. Pre-surgery MRI images were acquired with a 3.0 GE Signa HDxt scanner with an 8-channel head-coil in the department of radiology, Yijishan Hospital of Wannan Medical College. High-resolution T1-weighted MR images were obtained by a 3D magnetization-prepared rapid gradient-echo (MPRAGE) with the following parameters: repeat time (TR) = 1900 milliseconds (ms), echo time (TE) = 2.49 ms, time inversion (TI) = 900 ms, matrix = 256 × 256, flip angle (FA) = 90°, thickness = 1 millimeter (mm), gap = 0.5 mm, slices = 176. Resting-state functional images, including 240 volumes, were obtained using a gradient-recalled echo-planar imaging (GRE-EPI) sequence, with TR = 2000 ms, TE = 30 ms, FA = 90°, acquisition matrix = 64 × 64, field of view (FOV) = 220 mm × 220 mm, thickness = 4.0 mm, gap = 0 mm, number of slices = 36, and voxel size = 3.4 mm × 3.4 mm × 4 mm.
Image Pre-processing
Magnetic resonance imaging data were pre-processed using Matlab (Math Works Inc., Natick, MA, United States) and SPM8. The image processing procedure was conducted per published studies (Liu Y. et al., 2020; Chen et al., 2022). Briefly, the functional MRI (fMRI) image pre-processing steps included discarding the first 10 volumes, slice-timing, head motion corrections, spatial normalization, spatial smoothing, denoising, and temporal filtering (Liu Y. et al., 2020; Chen et al., 2022). A total of 2 non-FGep patents and 1 FGep patent were excluded from the subsequent analysis due to the excessive head movement (cumulative translation or rotation > 3.0 mm or 3.0°). We found no significant differences in head motion parameters between the groups (Power et al., 2012; Van Dijk et al., 2012).
During fMRI image pre-processing, the Diffeomorphic Anatomical Registration Through Exponentiated Lie Algebra (DARTEL) algorithm was used to normalize and segment the structural images into GM, WM, and cerebrospinal fluid (CSF) partitions (Ashburner and Friston, 2009). The native and DARTEL versions were used to compute total intracranial volumes (TIV). TIV was obtained by calculating the sum of GM, WM, and CSF tissues of all voxels in native space by using internal code. The structural and functional image processing procedure details are provided in the Supplementary Methods 1.
Executive Control Network Functional Connectivity Analysis
Referring to previously published studies (Zhang et al., 2018), we made patients’ group GM mask without tumor area before performing FC analysis. Subsequent FC analysis was constrained within this GM mask. The details regarding the construction of group GM mask without tumor area are provided in the Supplementary Methods 2.
Several studies have agreed that the dorsolateral prefrontal cortex (DLPFC) is the critical brain region within the ECN (Seeley et al., 2007; Fransson and Marrelec, 2008; Smith et al., 2009). According to converging evidence (Seeley et al., 2007; Fransson and Marrelec, 2008; Smith et al., 2009; Dong et al., 2020; Liu et al., 2021), this study selected the right DLPFC (MNI space: 48, 12, 34) to draw 6-mm spheres as seed region of interest (ROI) (Seeley et al., 2007; Fransson and Marrelec, 2008; Smith et al., 2009; Dong et al., 2020; Liu et al., 2021). We performed a voxel-wise cross-correlation analysis between the individual averaged time series for all voxel within DLPFC and the whole brain within the patients’ group GM mask (Chen et al., 2016a,2022; Liu et al., 2021).
Regional Homogeneity Analysis
As described in previous studies (Zang et al., 2004; Liu Y. et al., 2020), we measured ReHo to characterize the homogeneity of the time series in a local neighborhood of voxels within the ECN. The details of the ReHo analysis are provided in Supplementary Methods 3.
Pattern Classification Based on the Altered Gray Matter, Regional Homogeneity, and Functional Connectivity Within Executive Control Network
As described in previous studies (Chen et al., 2022), an support vector machine (SVM) approach was applied to test how well GM, ReHo, and FC within ECN could distinguish non-FGep and FGep from CN subjects. We performed a linear SVM classifier using the LIBSVM software.1 We assess the generalization of this SVM classifier and its accuracy, sensitivity, and specificity using a leave-one-out cross-validation (LOOCV) strategy. Accuracies obtained for each tested subject were averaged to obtain the LOOCV accuracy. The AUC was used to assess the classification power of the SVM classifier. The details of the LOOCV strategy are provided in Supplementary Methods 4.
Statistical Analysis
Demographic, Neuropsychological Data, and Head Rotation Parameters
We performed a one-way analysis of variance (ANOVA) and chi-square test (only applied in the comparison of gender) to compare differences in demographic data (age, gender, and education level), cognitive measures, and head rotation parameters among CN, non-FGep, and FGep subjects (p < 0.05). For each neuropsychological test, the individual raw scores were transformed to Z scores to improve the statistical power for cognitive measures, according to the mean and standard deviation of the scores for all subjects. Notably, for tests measured in time (TMT), the raw scores were defined as the reciprocal of the time required for the test. Post hoc analysis was performed to investigate the source of ANOVA differences.
Group-Level Intrinsic Connectivity Analysis
To determine the ECN patterns for the CN, non-FGep, and FGep subjects, the FC spatial maps of each group were submitted to a random-effect analysis using a one-sample t-test after controlling for age, sex, education, and TIV. A stringent threshold of p < 0.001 was set using the permutation test with Threshold-Free Cluster Enhancement (TFCE) (Smith and Nichols, 2009) and the family wise error (FWE) correction for the whole brain to reveal regions that are most robustly correlated with each seed.
We used two-sample t-tests to assess the between-group differences of FCs of ECN after controlling for the effects of age, sex, education, and TIV. A stringent threshold of p < 0.001 was set using the permutation test with TFCE and the false discovery rate (FDR) correction.
Clinical Behavioral Significance of Altered Gray Matter, Regional Homogeneity, and Functional Connectivity
To assess the clinical behavioral significance of altered GM, ReHo, and FC, we extracted the values of the regions showing differential patterns among CN, non-FGep, and FGep subjects. We performed a correlation analysis to assess the relationships between the extracted GM, ReHo, and FC values and cognition after controlling for age, sex, and education. The significance level was set to p < 0.05 using FDR correction for multiple comparisons.
Results
Demographic and Neuropsychological Data
As shown in Table 1, no significant differences in age, gender, and education level were found between CN, non-FGep, and FGep groups (p > 0.05). The non-FGep showed significant deficits in information processing indicated by lower TMT scores. Similarly, FGep showed a substantial decrease in episodic memory (lower AVLT scores) and information processing speed (lower TMT scores) relative to CN (p-values < 0.05).
Functional Connectivity Patterns and Inter-Group Differences Within the Executive Control Network
Functional connectivity patterns within the ECN are shown in Figure 1. CN, non-FGep, and FGep showed differential FC patterns within ECN (p < 0.001 corrected by TFCE-FWE). Compared with CN, non-FGep had reduced FC between the right DLPFC and left precentral gyrus (PreCG), insula (INS), middle cingulate and paracingulate gyri (MCC), postcentral gyrus (PoCG), inferior parietal gyrus (IPG), and bilateral supramarginal gyrus (SMG) within the ECN (Figure 2A and Table 2, p < 0.001 corrected by TFCE-FDR). The FGep group had decreased FC between the right DLPFC and left PreCG, bilateral superior frontal gyrus (SFGdor), bilateral middle frontal gyrus (MFG), right inferior frontal gyrus, opercular part (IFGoperc), right inferior frontal gyrus, and triangular part (IFGtriang) within the ECN, compared to CN (Figure 2B, p < 0.001 corrected by TFCE-FDR). Comparing non-FGep and FGep, within the ECN, FGep had diminished FC between the right DLPFC and the right ANG and bilateral SFGdor, and increased FC between the right DLPFC and left PoCG (Figure 2C, p < 0.001 corrected by TFCE-FDR).
Figure 1. Resting-state functional connectivity patterns of the executive control network within group maps in CN (A), non-FGep (B), and FGep (C) subjects. CN, healthy controls; non-FGep, frontal glioma without epilepsy; FGep, frontal glioma with epilepsy; PreCG, Precentral gyrus; L, left hemisphere; R, right hemisphere.
Figure 2. Comparison of FC within the executive control network between CN, non-FGep, and FGep subjects. (A) Altered regions of FC within the executive control network in non-FGep patients compared with CN. The bar chart on the right shows the quantitative comparison of FC in these altered regions. ***p < 0.001. (B) Altered regions of FC within the executive control network in FGep patients compared with CN. The bar chart on the right shows the quantitative comparison of FC in these altered regions. ***p < 0.001. (C) Different regions of FC within the executive control network between FGep and non-FGep patients. The bar chart on the right shows the quantitative comparison of FC in these different regions. ***p < 0.001. CN, healthy controls; non-FGep, frontal glioma without epilepsy; FGep, frontal glioma with epilepsy; L, left hemisphere; R, right hemisphere; PreCG, Precentral gyrus; INS, Insula; MCC, Middle Cingulate and paracingulate gyri; PoCG, Postcentral gyrus; IPG, Inferior parietal gyrus; SMG, Supramarginal gyrus; SFGdor, Superior frontal gyrus, dorsolateral; MFG, Middle frontal gyrus; IFGoperc, Inferior frontal gyrus, opercular part; IFGtriang, Inferior frontal gyrus, triangular part; ANG, Angular gyrus; FC, functional connectivity.
Differences in Gray Matter and Regional Homogeneity Within the Executive Control Network
Compared with CN, non-FGep showed increased GM in left PoCG and left IPG (Figure 3A and Table 2, p < 0.05 corrected by FDR) and increased ReHo in left INS and left PreCG (Figure 3C and Table 2, p < 0.05 corrected by FDR). The FGep showed increased GM in left PreCG and right SFGdor than CN (Figure 3B and Table 2, p < 0.05 corrected by FDR).
Figure 3. Comparison of ReHo and GM within the executive control network among CN, non-FGep, and FGep subjects. (A) Altered regions of GM within the executive control network in non-FGep patients compared with CN. The bar chart on the right shows the quantitative comparison of GM in these altered regions. ***p < 0.001. (B) Altered regions of GM within the executive control network in FGep patients compared with CN. The bar chart on the right shows the quantitative comparison of FC in these altered regions. ***p < 0.001. (C) Altered regions of ReHo within the executive control network in non-FGep patients compared to CN. The bar chart on the right shows the quantitative comparison of ReHo in these different regions. ***p < 0.001. CN, healthy controls; non-FGep, frontal glioma without epilepsy; FGep, frontal glioma with epilepsy; PreCG, Precentral gyrus; INS, Insula; PoCG, Postcentral gyrus; IPG, Inferior parietal gyrus; SFGdor, Superior frontal gyrus, dorsolateral; GM, gray matter; ReHo, regional homogeneity.
Clinical Behavioral Significance of Altered Gray Matter, Regional Homogeneity, and Functional Connectivity
The correlation analysis demonstrated that FC between right DLPFC and left PoCG within the ECN was positively correlated with the AVLT score in FGep and non-FGep (Figure 4, r = 0.390, p = 0.007, FDR correction). Moreover, no significant correlations were found between GM, ReHo, FC, and other cognitive functions (p-values < 0.05).
Figure 4. Clinical behavioral significance of brain regions with differential FC between non-FGep and FGep subjects. non-FGep, frontal glioma without epilepsy; FGep, frontal glioma with epilepsy; SFG, superior frontal gyrus; FC, functional connectivity; AVLT, Auditory Verbal Learning Test- 20-min delayed recall.
Classification of CN, Non-FGep, and FGep Patients Based on the Altered Gray Matter Volumes, Regional Homogeneity, and Functional Connectivity
As shown in Figure 5, the SVM classifier’s classification accuracy was 94.12% for non-FGep from CN, 86.67% for FGep from CN, and 76% for FGep from non-FGep. The ROC curve of the SVM classifier showed a high power for distinguishing non-FGep from CN, FGep from CN, and FGep from non-FGep on an individual subject basis, with high AUC (99.0, 89.0, and 83.0%, respectively), sensitivity (97.0, 82.0, and 76.0%, respectively), and specificity (97.0, 86.0, and 86.0%, respectively).
Figure 5. Classification power of the MRI-based “classifier” in distinguishing individuals as non-FGep vs. CN, FGep vs. CN, and non-FGep vs. FGep. (A) The schematic diagram indicated a support vector machine “classifier” based on altered structural (GM) and functional (ReHo and FC) characteristics within the executive control network. (B) ROC curve showing the classification power of the MRI-based “classifier.” The values of ACC, AUC, sensitivity, and specificity are marked in the lower right of the figure. CN, healthy controls; non-FGep, frontal glioma without epilepsy; FGep, frontal glioma with epilepsy; ACC, accuracy; AUC, area under the curve; GM, gray matter; ReHo, regional homogeneity; FC, functional connectivity.
Discussion
To our best knowledge, this was the first study to evaluate structural and functional abnormalities of ECN in patients with frontal glioma-related epilepsy and to investigate whether the abnormalities can accurately identify glioma-related seizures. Three main results of the study should be emphasized. First, compared with CN, non-FGep and FGep showed convergent and divergent altered patterns of ECN based on the structural and functional level of the brain. Second, FC between right SFG and right DLPFC was positively correlated with the AVLT scores in non-FGep and FGep. Lastly, the SVM classifier based on structural and functional abnormalities within the ECN could accurately differentiate non-FGep and FGep from CN, and FGep from non-FGep on an individual basis with very high accuracy, AUC, sensitivity, and specificity. These findings indicate that specific functional changes can assist in the early identification and intervention of epileptic seizures in frontal glioma patients.
This study found that non-FGep had reduced FC within ECN, mainly distributed in PreCG, INS, MCC, PoCG, IPG, and SMG. FGep also showed decreased FC within ECN, distributed primarily in PreCG, SFGdor, MFG, IFGoperc, and IFGtriang. These results suggest that non-FGep and FGep showed convergent and divergent altered patterns of ECN. These ECN network findings are supported by behavioral results with significantly lower executive function scores (TMT scores) in both non-FGep and FGep. Our results are inconsistent with previous studies (Liu et al., 2019; Liu Y. et al., 2020; Tordjman et al., 2021), which showed that frontal glioma patients have the clinical compensation phenomenon of intact executive function and increased FC within ECN. Their interpretation is that tumor invasion induces FC reorganization of the ECN. This inconsistency could be because compensation and decompensation coexist dynamically in disease progression, and our study subjects may have been in the decompensation stage (Belleville et al., 2011; Chen et al., 2015, 2016a; Liu Y. et al., 2020). Another reason could be that the method for driving FC analysis of ECN in this study is inconsistent with their approach. Our research adopts a seed-based, data-driven approach in the whole brain for calculating FC. In contrast, their research used an independent component analysis and FC between regions within the ECN mask. Furthermore, our findings suggest that tumor invasion in non-FGep was more likely to diminish FC in remote brain regions, while tumor invasion in FGep was more likely to decrease FC in nearby brain regions. These findings support the idea that tumor invasion of cortical lesions leads to the disruption of connections between local and remote brain regions in the form of functional networks (Vigneau et al., 2006; Carter et al., 2010; Xiang et al., 2010; Briganti et al., 2012).
Notably, FGep showed decreased FC in ANG, PoCG, and SFG.dor within ECN compared to non-FGep and CN. Contrastingly, non-FGep showed no difference in FC in SFG.dor. FC between right SFG and right DLPFC was also positively correlated with AVLT scores in non-FGep and FGep. These findings suggest that tumor invasion of the frontal lobe accompanied by epileptic seizures exacerbates functional abnormalities within ECN, and the reduced FC contributes to episodic memory impairment in FGep. Clinical behavioral results support this theory, wherein non-FGep showed intact episodic memory indicated by AVLT scores, while FGep showed significant episodic memory impairment. Indeed, previous studies have confirmed the interaction between the ECN and the episodic memory network with the frontal lobe as the interactive hub (Yuan et al., 2016; Xu et al., 2020). Therefore, we propose that the frontal lobe (SFG) may be an early neuroimaging marker of episodic memory decline during epileptic episodes in patients with frontal glioma.
Another unique contribution of this study is that tumor invasion of the frontal lobe induced local structural (GM volume) and functional (ReHo) reorganization only, regardless of an epileptic seizure. These results align with previous studies that found that tumor invasion induces the brain’s structural and functional reorganization (Liu et al., 2019; Hu et al., 2020; Liu D. et al., 2020; Liu Y. et al., 2020). However, unlike these studies, our results found reorganization occurred at perilesional and remote recruitment (Duffau, 2014; Almairac et al., 2018). It is also possible that different isocitrate dehydrogenase 1 (IDH1) mutation statuses lead to different patterns of network reorganization (Qi et al., 2022). This study did not identify the IDH1 mutation status. Interestingly, reorganization of tumor invasion with seizure occurred in perilesional recruitment, whereas reorganization without seizure occurred in remote recruitment. This suggests that the seizures may have more to do with the tumor itself. Therefore, seizure-induced brain reorganization is more likely to occur around the lesion.
This study also showed that non-FGep and FGep showed compensatory activity in the same brain region: PreCG. The PreCG is associated with sensorimotor and visuo-spatial functions (Davion et al., 2021; Li et al., 2022). Our behavioral results showed that visuospatial cognitive function was intact in both groups. Therefore, we suggest that tumor invasion induced structural and functional reorganization of the PreCG compensated for the visuospatial cognitive function of the patients. Although the SFG.dor was compensated in FGep, its FC decreased significantly, and episodic memory was impaired. This indicates that seizures in glioma patients may accelerate the progression from the compensatory to the decompensated stage.
Our most fascinating finding was that by combining structural and functional abnormalities within the ECN, the SVM classifier could accurately differentiate non-FGep and FGep from CN and FGep from non-FGep on an individual basis with very high accuracy, AUC, sensitivity, and specificity. Structural and functional features within ECN may reflect epilepsy susceptibility in patients with frontal glioma and can effectively detect frontal glioma-related epilepsy. Although, in recent years, some studies have constructed prediction models of epileptic seizures in glioma by using radiomic features and ML (Liu et al., 2018; Gao et al., 2021), they only use the information of segmented and masked tumor lesions. These models have some limitations. First, there are many manual factors in the segmentation and selection of tumors. Second, most glioma-related seizures are not located in the lesion but the peripheral or remote brain region. The radiomic features based on tumor lesions cannot truly reflect the characteristics of epileptic seizures.
Furthermore, some studies suggest that tumor invasion affects brain function over the whole network rather than isolated regions (Vigneau et al., 2006; Xiang et al., 2010; Briganti et al., 2012). In particular, some studies indicated that tumor-related epilepsy induces topological changes in WM networks (Zhou et al., 2022) and functional networks (Fang et al., 2022). Therefore, the network-based features and the structural features outside the tumor lesion can reflect more information about glioma-related epilepsy. The effective neuroimaging features will be highly conducive to the prediction of seizures. In addition, radiomic analysis can incorporate too many features in a small sample of data, which can easily result in over-fitting or bias of the prediction model. Our study included the structural and functional abnormalities within the ECN that have the highest association with the impact of tumor invasion of the frontal lobe as features of the model, which improved the predictive power and generalization of the model.
However, this study has certain limitations: First, a relatively small sample size, especially for patients with glioma with seizures (n = 21). Larger samples will be needed in the future to verify our results. Second, we only evaluated the GM volume, and future studies need more indicators to assess the structure, such as fiber integrity measured by diffusion tensor technology. Finally, this study was a single-center study that could not evaluate the model’s generalization. The stability of the model needs to be verified in independent samples. Future multi-center studies are needed to verify the generalization of our model in distinguishing glioma epileptic seizures.
Conclusion
Tumor invasion of the frontal lobe induces local structural and functional reorganization within the ECN. These structural and functional abnormalities are exacerbated by accompanying epileptic seizures. ECN abnormalities, which are related to cognitive impairment, can be used to accurately distinguish the presence or absence of epileptic seizures in patients with frontal glioma at an early stage. Our findings raise the possibility that specific functional changes can assist in the early identification and intervention of epileptic seizures in frontal glioma patients.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by the Human Participants Ethics Committee of Yijishan Hospital of Wannan Medical College (China). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
GD and XJ designed the study. MT, RX, WZ, KD, and ZH collected the data. GD, RX, XC, and HZ analyzed the data and prepared the manuscript. All authors read and approved the final manuscript.
Funding
This study was provided by the grant from the Funding of “Peak” Training Program for Scientific Research of Yijishan Hospital (GF2019G01 and KPF2019003).
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.
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.
Acknowledgments
We thank all neurosurgical nurses and doctors for their help in this study.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2022.916771/full#supplementary-material
Footnotes
References
Abbasi, B., and Goldenholz, D. M. (2019). Machine learning applications in epilepsy. Epilepsia 60, 2037–2047.
Abdelaziz Ismael, S. A., Mohammed, A., and Hefny, H. (2020). An enhanced deep learning approach for brain cancer MRI images classification using residual networks. Artif. Intell. Med. 102:101779. doi: 10.1016/j.artmed.2019.101779
Akeret, K., Stumpo, V., Staartjes, V. E., Vasella, F., Velz, J., Marinoni, F., et al. (2020). Topographic brain tumor anatomy drives seizure risk and enables machine learning based prediction. NeuroImage Clin. 28:102506. doi: 10.1016/j.nicl.2020.102506
Almairac, F., Duffau, H., and Herbet, G. (2018). Contralesional macrostructural plasticity of the insular cortex in patients with glioma: a VBM study. Neurology 91, e1902–e1908. doi: 10.1212/WNL.0000000000006517
Ashburner, J., and Friston, K. J. (2009). Computing average shaped tissue probability templates. NeuroImage 45, 333–341. doi: 10.1016/j.neuroimage.2008.12.008
Bangalore Yogananda, C. G., Shah, B. R., Vejdani-Jahromi, M., Nalawade, S. S., Murugesan, G. K., Yu, F. F., et al. (2020). A novel fully automated MRI-based deep-learning method for classification of IDH mutation status in brain gliomas. Neuro Oncol. 22, 402–411. doi: 10.1093/neuonc/noz199
Belleville, S., Clement, F., Mellah, S., Gilbert, B., Fontaine, F., and Gauthier, S. (2011). Training-related brain plasticity in subjects at risk of developing Alzheimer’s disease. Brain 134, 1623–1634. doi: 10.1093/brain/awr037
Brier, M. R., Thomas, J. B., Snyder, A. Z., Benzinger, T. L., Zhang, D., Raichle, M. E., et al. (2012). Loss of intranetwork and internetwork resting state functional connections with Alzheimer’s disease progression. J. Neurosci. 32, 8890–8899. doi: 10.1523/JNEUROSCI.5698-11.2012
Briganti, C., Sestieri, C., Mattei, P. A., Esposito, R., Galzio, R. J., Tartaro, A., et al. (2012). Reorganization of functional connectivity of the language network in patients with brain gliomas. AJNR Am. J. Neuroradiol. 33, 1983–1990. doi: 10.3174/ajnr.A3064
Cao, X., Qian, Z., Xu, Q., Shen, J., Zhang, Z., and Lu, G. (2014). Altered intrinsic connectivity networks in frontal lobe epilepsy: a resting-state fMRI study. Comput. Math. Methods Med. 2014:864979. doi: 10.1155/2014/864979
Carter, A. R., Astafiev, S. V., Lang, C. E., Connor, L. T., Rengachary, J., Strube, M. J., et al. (2010). Resting interhemispheric functional magnetic resonance imaging connectivity predicts performance after stroke. Ann. Neurol. 67, 365–375. doi: 10.1002/ana.21905
Chen, J., Chen, R., Xue, C., Qi, W., Hu, G., Xu, W., et al. (2022). Hippocampal-Subregion mechanisms of repetitive transcranial magnetic stimulation causally associated with amelioration of episodic memory in amnestic mild cognitive impairment. J. Alzheimers Dis. JAD 85, 1329–1342. doi: 10.3233/jad-210661
Chen, J., Duan, X., Shu, H., Wang, Z., Long, Z., Liu, D., et al. (2016a). Differential contributions of subregions of medial temporal lobe to memory system in amnestic mild cognitive impairment: insights from fMRI study. Sci. Rep. 6:26148. doi: 10.1038/srep26148
Chen, J., Shu, H., Wang, Z., Zhan, Y., Liu, D., Liao, W., et al. (2016b). Convergent and divergent intranetwork and internetwork connectivity patterns in patients with remitted late-life depression and amnestic mild cognitive impairment. Cortex 83, 194–211. doi: 10.1016/j.cortex.2016.08.001
Chen, J., Zhang, Z., and Li, S. (2015). Can multi-modal neuroimaging evidence from hippocampus provide biomarkers for the progression of amnestic mild cognitive impairment? Neurosci. Bull. 31, 128–140. doi: 10.1007/s12264-014-1490-8
Choi, Y. S., Ahn, S. S., Chang, J. H., Kang, S. G., Kim, E. H., Kim, S. H., et al. (2020). Machine learning and radiomic phenotyping of lower grade gliomas: improving survival prediction. Eur. Radiol. 30, 3834–3842. doi: 10.1007/s00330-020-06737-5
Crone, E. A., Wendelken, C., Donohue, S. E., and Bunge, S. A. (2006). Neural evidence for dissociable components of task-switching. Cereb. Cortex 16, 475–486. doi: 10.1093/cercor/bhi127
Davion, J. B., Lopes, R., Jougleux, C., Viard, R., Dumont, J., Leclerc, X., et al. (2021). Brief International Cognitive Assessment for Multiple Sclerosis scores are associated with the cortical thickness of specific cortical areas in relapsing-remitting patients. Rev. Neurol. 178, 326–336. doi: 10.1016/j.neurol.2021.06.014
Dong, W., Qiu, C., Jiang, X., Shen, B., Zhang, L., Liu, W., et al. (2020). Can the executive control network be used to diagnose Parkinson’s disease and as an efficacy indicator of deep brain stimulation? Parkinsons Dis. 2020:6348102. doi: 10.1155/2020/6348102
Dosenbach, N. U., Fair, D. A., Miezin, F. M., Cohen, A. L., Wenger, K. K., Dosenbach, R. A., et al. (2007). Distinct brain networks for adaptive and stable task control in humans. Proc. Natl. Acad. Sci. U.S.A. 104, 11073–11078. doi: 10.1073/pnas.0704320104
Duffau, H. (2014). Diffuse low-grade gliomas and neuroplasticity. Diagn. Interv. Imaging 95, 945–955. doi: 10.1016/j.diii.2014.08.001
Fang, S., Li, L., Weng, S., Guo, Y., Zhang, Z., Wang, L., et al. (2022). Decreasing shortest path length of the sensorimotor network induces frontal Glioma-related epilepsy. Front. Oncol. 12:840871. doi: 10.3389/fonc.2022.840871
Fisher, R. S., Cross, J. H., French, J. A., Higurashi, N., Hirsch, E., Jansen, F. E., et al. (2017). Operational classification of seizure types by the international league against epilepsy: position paper of the ILAE commission for classification and terminology. Epilepsia 58, 522–530. doi: 10.1111/epi.13670
Fransson, P., and Marrelec, G. (2008). The precuneus/posterior cingulate cortex plays a pivotal role in the default mode network: evidence from a partial correlation network analysis. NeuroImage 42, 1178–1184. doi: 10.1016/j.neuroimage.2008.05.059
Gao, A., Yang, H., Wang, Y., Zhao, G., Wang, C., Wang, H., et al. (2021). Radiomics for the prediction of epilepsy in patients with frontal glioma. Front. Oncol. 11:725926. doi: 10.3389/fonc.2021.725926
Gleichgerrcht, E., Munsell, B., Bhatia, S., Vandergrift, W. A. III, Rorden, C., McDonald, C., et al. (2018). Deep learning applied to whole-brain connectome to determine seizure control after epilepsy surgery. Epilepsia 59, 1643–1654. doi: 10.1111/epi.14528
Hu, G., Hu, X., Yang, K., Liu, D., Xue, C., Liu, Y., et al. (2020). Restructuring of contralateral gray matter volume associated with cognition in patients with unilateral temporal lobe glioma before and after surgery. Hum. Brain Mapp. 41, 1786–1796. doi: 10.1002/hbm.24911
Kinno, R., Muragaki, Y., Maruyama, T., Tamura, M., Tanaka, K., Ono, K., et al. (2020). Differential effects of a left frontal glioma on the cortical thickness and complexity of both Hemispheres. Cereb. Cortex Commun. 1:tgaa027. doi: 10.1093/texcom/tgaa027
Li, G., Jiang, S., Meng, J., Chai, G., Wu, Z., Fan, Z., et al. (2022). Assessing differential representation of hand movements in multiple domains using stereo-electroencephalographic recordings. NeuroImage 250:118969. doi: 10.1016/j.neuroimage.2022.118969
Li, Q., Chen, Y., Wei, Y., Chen, S., Ma, L., He, Z., et al. (2017). Functional network connectivity patterns between idiopathic generalized epilepsy with myoclonic and absence seizures. Front. Comput. Neurosci. 11:38. doi: 10.3389/fncom.2017.00038
Liu, D., Chen, J., Hu, X., Hu, G., Liu, Y., Yang, K., et al. (2020). Contralesional homotopic functional plasticity in patients with temporal glioma. J. Neurosurg. 134, 417–425. doi: 10.3171/2019.11.JNS191982
Liu, D., Hu, X., Liu, Y., Yang, K., Xiao, C., Hu, J., et al. (2019). Potential intra- or cross-network functional reorganization of the triple unifying networks in patients with frontal glioma. World Neurosurg. 128, e732–e743. doi: 10.1016/j.wneu.2019.04.248
Liu, W., Liu, L., Cheng, X., Ge, H., Hu, G., Xue, C., et al. (2021). Functional integrity of executive control network contributed to retained executive abilities in mild cognitive impairment. Front. Aging Neurosci. 13:710172. doi: 10.3389/fnagi.2021.710172
Liu, Y., Hu, G., Yu, Y., Jiang, Z., Yang, K., Hu, X., et al. (2020). Structural and functional reorganization within cognitive control network associated with protection of executive function in patients with unilateral frontal gliomas. Front. Oncol. 10:794. doi: 10.3389/fonc.2020.00794
Liu, Z., Wang, Y., Liu, X., Du, Y., Tang, Z., Wang, K., et al. (2018). Radiomics analysis allows for precise prediction of epilepsy in patients with low-grade gliomas. NeuroImage Clin. 19, 271–278. doi: 10.1016/j.nicl.2018.04.024
Maschio, M., and Dinapoli, L. (2012). Patients with brain tumor-related epilepsy. J. Neuro Oncol. 109, 1–6. doi: 10.1016/b978-0-12-417043-8.00001-8
Munsell, B. C., Wee, C. Y., Keller, S. S., Weber, B., Elger, C., da Silva, L. A., et al. (2015). Evaluation of machine learning algorithms for treatment outcome prediction in patients with epilepsy based on structural connectome data. NeuroImage 118, 219–230. doi: 10.1016/j.neuroimage.2015.06.008
Ostrom, Q. T., Cioffi, G., Gittleman, H., Patil, N., Waite, K., Kruchko, C., et al. (2019). CBTRUS statistical report: primary brain and other central nervous system tumors diagnosed in the United States in 2012-2016. Neuro Oncol. 21, v1–v100. doi: 10.1093/neuonc/noz150
Pallud, J., Audureau, E., Blonski, M., Sanai, N., Bauchet, L., Fontaine, D., et al. (2014). Epileptic seizures in diffuse low-grade gliomas in adults. Brain 137, 449–462. doi: 10.1093/brain/awt345
Park, J. E., Kim, H. S., Kim, S. J., Kim, J. H., and Shim, W. H. (2015). Alteration of long-distance functional connectivity and network topology in patients with supratentorial gliomas. Neuroradiology 58, 311–320. doi: 10.1007/s00234-015-1621-6
Peeken, J. C., Goldberg, T., Pyka, T., Bernhofer, M., Wiestler, B., Kessel, K. A., et al. (2019). Combining multimodal imaging and treatment features improves machine learning-based prognostic assessment in patients with glioblastoma multiforme. Cancer Med. 8, 128–136. doi: 10.1002/cam4.1908
Power, J. D., Barnes, K. A., Snyder, A. Z., Schlaggar, B. L., and Petersen, S. E. (2012). Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. NeuroImage 59, 2142–2154. doi: 10.1016/j.neuroimage.2011.10.018
Qi, C., Wang, R., Meng, L., Li, S., and Li, Y. (2022). Functional reorganization of contralesional networks varies according to isocitrate dehydrogenase 1 mutation status in patients with left frontal lobe glioma. Neuroradiology doi: 10.1007/s00234-022-02932-x
Sakai, K., and Passingham, R. E. (2006). Prefrontal set activity predicts rule-specific neural processing during subsequent cognitive performance. J. Neurosci. 26, 1211–1218. doi: 10.1523/JNEUROSCI.3887-05.2006
Sanghani, P., Ang, B. T., King, N. K. K., and Ren, H. (2018). Overall survival prediction in glioblastoma multiforme patients from volumetric, shape and texture features using machine learning. Surg. Oncol. 27, 709–714. doi: 10.1016/j.suronc.2018.09.002
Seeley, W. W., Crawford, R. K., Zhou, J., Miller, B. L., and Greicius, M. D. (2009). Neurodegenerative diseases target large-scale human brain networks. Neuron 62, 42–52. doi: 10.1016/j.neuron.2009.03.024
Seeley, W. W., Menon, V., Schatzberg, A. F., Keller, J., Glover, G. H., Kenna, H., et al. (2007). Dissociable intrinsic connectivity networks for salience processing and executive control. J. Neurosci. 27, 2349–2356. doi: 10.1523/JNEUROSCI.5587-06.2007
Siomin, V., Angelov, L., Li, L., and Vogelbaum, M. A. (2005). Results of a survey of neurosurgical practice patterns regarding the prophylactic use of anti-epilepsy drugs in patients with brain tumors. J. Neuro Oncol. 74, 211–215. doi: 10.1007/s11060-004-6912-4
Smith, S. M., and Nichols, T. E. (2009). Threshold-free cluster enhancement: addressing problems of smoothing, threshold dependence and localisation in cluster inference. NeuroImage 44, 83–98. doi: 10.1016/j.neuroimage.2008.03.061
Smith, S. M., Fox, P. T., Miller, K. L., Glahn, D. C., Fox, P. M., Mackay, C. E., et al. (2009). Correspondence of the brain’s functional architecture during activation and rest. Proc. Natl. Acad. Sci. U.S.A. 106, 13040–13045. doi: 10.1073/pnas.0905267106
Tan, Y., Zhang, S. T., Wei, J. W., Dong, D., Wang, X. C., Yang, G. Q., et al. (2019). A radiomics nomogram may improve the prediction of IDH genotype for astrocytoma before surgery. Eur. Radiol. 29, 3325–3337. doi: 10.1007/s00330-019-06056-4
Tordjman, M., Madelin, G., Gupta, P. K., Cordova, C., Kurz, S. C., Orringer, D., et al. (2021). Functional connectivity of the default mode, dorsal attention and fronto-parietal executive control networks in glial tumor patients. J. Neuro Oncol. 152, 347–355. doi: 10.1007/s11060-021-03706-w
Tremont-Lukats, I. W., Ratilal, B. O., Armstrong, T., and Gilbert, M. R. (2008). Antiepileptic drugs for preventing seizures in people with brain tumors. Cochrane Database Syst. Rev. 2008:CD004424. doi: 10.1002/14651858.CD004424.pub2
Vamvakas, A., Williams, S. C., Theodorou, K., Kapsalaki, E., Fountas, K., Kappas, C., et al. (2019). Imaging biomarker analysis of advanced multiparametric MRI for glioma grading. Phys. Med. 60, 188–198. doi: 10.1016/j.ejmp.2019.03.014
Van Dijk, K. R., Sabuncu, M. R., and Buckner, R. L. (2012). The influence of head motion on intrinsic functional connectivity MRI. NeuroImage 59, 431–438. doi: 10.1016/j.neuroimage.2011.07.044
Vigneau, M., Beaucousin, V., Herve, P. Y., Duffau, H., Crivello, F., Houde, O., et al. (2006). Meta-analyzing left hemisphere language areas: phonology, semantics, and sentence processing. NeuroImage 30, 1414–1432. doi: 10.1016/j.neuroimage.2005.11.002
Xiang, H. D., Fonteijn, H. M., Norris, D. G., and Hagoort, P. (2010). Topographical functional connectivity pattern in the perisylvian language networks. Cereb. Cortex 20, 549–560. doi: 10.1093/cercor/bhp119
Xu, W., Chen, S., Xue, C., Hu, G., Ma, W., Qi, W., et al. (2020). Functional MRI-specific alterations in executive control network in mild cognitive impairment: an ALE meta-analysis. Front. Aging Neurosci. 12:578863. doi: 10.3389/fnagi.2020.578863
You, G., Sha, Z., and Jiang, T. (2020). Clinical diagnosis and perioperative management of glioma-related epilepsy. Front. Oncol. 10:550353. doi: 10.3389/fonc.2020.550353
Yuan, B., Chen, J., Gong, L., Shu, H., Liao, W., Wang, Z., et al. (2016). Mediation of episodic memory performance by the executive function network in patients with amnestic mild cognitive impairment: a resting-state functional MRI study. Oncotarget 7, 64711–64725. doi: 10.18632/oncotarget.11775
Zang, Y., Jiang, T., Lu, Y., He, Y., and Tian, L. (2004). Regional homogeneity approach to fMRI data analysis. NeuroImage 22, 394–400. doi: 10.1016/j.neuroimage.2003.12.030
Zhang, H., Shi, Y., Yao, C., Tang, W., Yao, D., Zhang, C., et al. (2016). Alteration of the intra- and cross- hemisphere posterior default mode network in frontal lobe glioma patients. Sci. Rep. 6:26972. doi: 10.1038/srep26972
Zhang, N., Xia, M., Qiu, T., Wang, X., Lin, C. P., Guo, Q., et al. (2018). Reorganization of cerebro-cerebellar circuit in patients with left hemispheric gliomas involving language network: a combined structural and resting-state functional MRI study. Hum. Brain Mapp. 39, 4802–4819. doi: 10.1002/hbm.24324
Keywords: epileptic seizure, executive control network, functional connectivity, gray matter, frontal glioma, pattern classification, regional homogeneity
Citation: Di G, Tan M, Xu R, Zhou W, Duan K, Hu Z, Cao X, Zhang H and Jiang X (2022) Altered Structural and Functional Patterns Within Executive Control Network Distinguish Frontal Glioma-Related Epilepsy. Front. Neurosci. 16:916771. doi: 10.3389/fnins.2022.916771
Received: 10 April 2022; Accepted: 05 May 2022;
Published: 26 May 2022.
Edited by:
Jiaojian Wang, University of Electronic Science and Technology of China, ChinaReviewed by:
Yu-Chen Chen, Nanjing Medical University, ChinaRui Liu, Capital Medical University, China
Copyright © 2022 Di, Tan, Xu, Zhou, Duan, Hu, Cao, Zhang and Jiang. 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: Xiaochun Jiang, jiangxiaochun2019@hotmail.com