- 1Department of Neurology, Zhengzhou University People’s Hospital, Zhengzhou, Henan, China
- 2Department of Neurology, Henan Provincial People’s Hospital, People’s Hospital of Zhengzhou University, Zhengzhou, Henan, China
- 3Department of Rehabilitation, The First Hospital of Shanxi Medical University, Taiyuan, Shanxi, China
- 4School of Medical Engineering, Xinxiang Medical University, Xinxiang, Henan, China
- 5School of Life Sciences and Technology, Xi’an Jiaotong University, Xi’an, Shaanxi, China
- 6Department of Neurology, People’s Hospital of Henan University, Zhengzhou, Henan, China
Objective: Cognitive impairment (CI) is a common disorder in patients with epilepsy (PWEs). Objective assessment method for diagnosing CI in PWEs would be beneficial in reality. This study proposed to construct a diagnostic model for CI in PWEs using the clinical and the phase locking value (PLV) functional connectivity features of the electroencephalogram (EEG).
Methods: PWEs who met the inclusion and exclusion criteria were divided into a cognitively normal (CON) group (n = 55) and a CI group (n = 76). The 23 clinical features and 684 PLVEEG features at the time of patient visit were screened and ranked using the Fisher score. Adaptive Boosting (AdaBoost) and Gradient Boosting Decision Tree (GBDT) were used as algorithms to construct diagnostic models of CI in PWEs either with pure clinical features, pure PLVEEG features, or combined clinical and PLVEEG features. The performance of these models was assessed using a five-fold cross-validation method.
Results: GBDT-built model with combined clinical and PLVEEG features performed the best with accuracy, precision, recall, F1-score, and an area under the curve (AUC) of 90.11, 93.40, 89.50, 91.39, and 0.95%. The top 5 features found to influence the model performance based on the Fisher scores were the magnetic resonance imaging (MRI) findings of the head for abnormalities, educational attainment, PLVEEG in the beta (β)-band C3-F4, seizure frequency, and PLVEEG in theta (θ)-band Fp1-Fz. A total of 12 of the top 5% of features exhibited statistically different PLVEEG features, while eight of which were PLVEEG features in the θ band.
Conclusion: The model constructed from the combined clinical and PLVEEG features could effectively identify CI in PWEs and possess the potential as a useful objective evaluation method. The PLVEEG in the θ band could be a potential biomarker for the complementary diagnosis of CI comorbid with epilepsy.
Introduction
Cognitive impairment (CI) is one of the very common comorbidities occurring in 70–80% of patients with epilepsy (PWEs) (Helmstaedter and Witt, 2017). Previous studies have revealed several factors that may induce CI in PWEs, including age at onset, duration of illness, surgical head trauma, perinatal injury, temporal lobe epilepsy, hippocampal abnormalities, seizures, status epilepticus, medications, and psychiatric factors (Black et al., 2010; Titiz et al., 2014; Vrinda et al., 2019; Jarcuskova et al., 2020; Wang et al., 2020; Novak et al., 2022a). Furthermore, interictal epileptiform discharges (IEDs) in electroencephalogram (EEG) recordings are an important indicator of CI in PWEs (Ung et al., 2017; Gavrilovic et al., 2019; Balcik et al., 2020), but the exact role of EEG in diagnosing CI in such patients has rarely been studied.
Cognitive scales serve as the primary method for diagnosing CI, with the Montreal Cognitive Assessment (MoCA) scale considered the most appropriate and more sensitive than the Mini-Mental State Examination (MMSE) scale for screening cognitive impairment in epileptic individuals (Montano-Lozada et al., 2021; Huang et al., 2022; Novak et al., 2022b). Notably, the MoCA-30 point scale is superior to the MoCA-20 scale for CI assessment in clinical practices (Bergeron et al., 2017; Del Brutto et al., 2019; Rodrigues et al., 2020; Melikyan et al., 2021). However, the scale has some shortcomings, most notably its susceptibility to subjective factors from both patients and physicians, which may lead to errors in the test. Although the MoCA scale is well suited to screening for CI in epileptic patients, however, it is a generic neurological screening tool for cognitive assessments. Therefore, there is an urgent need for developing an efficient objective assessment indicator for cognitive functions, specifically for individuals with epileptic symptoms.
Electroencephalogram plays a vital role in the diagnosis and management of epilepsy, as it provides an objective and accurate response to functional changes in the brain, thus avoiding the influence of subjective factors in the patient. A growing body of research has demonstrated a strong correlation between altered cognitive functions and the neural connectivity of different brain regions (He et al., 2018; Fadaie et al., 2021; Duma et al., 2022). Functional connectivity is a type of neural connectivity that mediates the temporal correlation between neurophysiological events at different brain regions and is primarily used to measure the degree of dependency and correlation between the signals. The phase locking value (PLV) is one of the quantitative indicators for functional connectivity (Elahian et al., 2017; Duma et al., 2021). Furthermore, EEG-based functional connectivity is employed to predict vagus nerve stimulation (VNS) responsiveness in children with refractory epilepsies (Ma et al., 2022), as well as to diagnose CI in patients comorbid with Parkinson’s disease (PD) (Cai et al., 2021). However, this approach has not been applied to the diagnosis of cognitive dysfunctions in PWEs. The Adaptive Boosting (AdaBoost) and Gradient-Boosted Decision Trees (GBDT) are classic algorithms for ensemble learning (EL) and have been widely used in areas of neurologic disorders such as epilepsy, Alzheimer’s disease (AD), PD, etc. (Peng et al., 2020; Wenbo et al., 2021; Zhang S. et al., 2021; Edeh et al., 2022). These follow the models constructed based on the clinical and PLVEEG functional connectivity features of EL algorithms and have shown the potential of an efficient objective evaluation tool for diagnosing CI in PWEs.
Here, we used EL algorithms to construct three distinct models for the diagnosis of CI in PWEs, purely based on the clinical and PLVEEG features. Additionally, we investigated to identify potential biomarkers for the diagnosis of cognitive functions in PWEs.
Materials and methods
Selection of the participants
A total of 131 PWEs from the outpatient clinic of the Epilepsy Center of Henan Provincial People’s Hospital between June 2018 and May 2022 were retrospectively screened and enrolled in the study. The inclusion criteria were: (1) the patient must meet the criteria of the International League Against Epilepsy (ILAE) for the diagnosis of epilepsy, seizures, and other epileptic syndromes (Fisher et al., 2014); (2) the age range at the time of consultation must be 12–60 years; (3) the patient must had a MoCA test at the time of consultation and should not have any history of MoCA scale testing in the last year; (4) at least 20 min of outpatient scalp EEG at the time of consultation, along with the availability of retrospective EEG data; and (5) the patient must have a complete clinical history and previous cranial MRI findings. Subjects were excluded if: (1) the patient’s age was less than 12 years or more than 60 years at the time of consultation; (2) the patient was diagnosed with psychogenic non-epileptic seizures, or epilepsy syndrome; (3) the patient was treated with drugs other than antiseizures medications that affect cognitive functions, such as benzodiazepines, anti-psychotics, and memory-enhancing drugs, at the time of consultation; and (4) the patient was missing the 20-min EEG recording data at the time of the enrollment.
Based on the patients’ MoCA scores during their visits to the epilepsy clinic, 131 PWEs were recruited for the study and were subsequently divided into the control (CON) group (MoCA ≥ 26; n = 55) and the CI group (MoCA < 26; n = 76) (Figure 1 and Table 1). The study was approved by the Ethics Committee of Henan Provincial People’s Hospital and all eligible subjects signed the written informed consent before their final recruitment to the study.
Figure 1. Flow chart. PWEs, patients with epilepsy; MoCA, Montreal Cognitive Assessment; EEG, electroencephalogram; CI, cognitive impairment; CON, cognitively normal; EL, ensemble learning.
Clinical features
Based on the patients’ medical history and clinical investigations at the time of the current visit to the epilepsy clinic, 23 clinical features were identified, in conjunction with previous studies: (1) age; (2) age at the first onset; (3) time from the first onset to current visit (Black et al., 2010); (4) gender; (5) family history of epilepsy (defined as whether a first or second degree relative had epilepsy); (6) history of previous head surgery or trauma; (7) history of previous the central nervous system (CNS) infections; (8) history of perinatal injuries due to premature birth, obstructed labor, hypoxia, and/or intracranial hematoma; (9) TLE; (10) MRI of the head for abnormalities; (11) hippocampal atrophy, or sclerosis (Titiz et al., 2014); (12) different types of seizures like generalized, focal, or both; (13) status epilepticus; (14) generalized tonic-clonic seizures (GTCS); (15) seizure frequency in the last year (Wang et al., 2020) (rare: ≤1 event; occasional: 2–3 events; frequent: ≥4 events); (16) class of antiseizures medications (Wang et al., 2020); (17) valproate (VPA) therapy in the last year; (18) phenytoin (PHT) therapy in the last year; (19) topiramate (TPM) therapy in the last year; (20) aura of epilepsy; (21) anxiety [according to the Hamilton Anxiety Inventory (HAI) scale rating: none, possible, definitely, or definitely obvious]; (22) depression [according to the Hamilton Depression Inventory (HDI) scale rating: none, possible, or definite]; and (23) educational attainment (≤6 years, 7–9 years, 10–12 years, or ≥13 years) (Table 2).
EEG acquisition and preprocessing
All patients in both CON and CI groups had scalp EEG recordings monitored for at least 20 min during this visit. All tests were performed in the awake closed-eye state, while EEG recordings performed during the sleep and awake open-eye states were excluded. The EEG-1200°C machine (Nihon Kohden, Tokyo, Japan), with a sampling frequency of 256 Hz, an amplification multiplier of 1000×, a low-pass filter of 70 Hz, and a high-pass filter of 0.5 Hz, was used for this study. This system uniformly used the international 10–20 lead system for placing the scalp electrodes, including 19 recording leads, namely Fp1, Fp2, Fz, Cz, Pz, C3, C4, T3, T4, T5, T6, F3, F4, F7, F8, O1, O2, P3, and P4, and 2 reference leads A1 and A2.
Preprocessing of EEG data was performed using the EEGLAB toolbox in MATLAB software (Mathworks Inc., USA) (Delorme and Makeig, 2004). Briefly, the EEG recordings were first filtered to extract only the 0.5–30 Hz recordings. Afterward, the artifacts of eye movements in electromyogram (EMG) were removed using independent component analysis. Finally, the 20-min EEG recording of each patient was intercepted into 6 s segments, and PLVEEG features were extracted.
Parameters setting for AdaBoost and GBDT
AdaBoost and GBDT are typical methods of boosting algorithm. In the AdaBoost model, the number and learning rate of base classifiers were also determined by grid search, ranging from 50 to 150 and 0 to 1, respectively and the algorithm of AdaBoost set to SAMME.R. The base classifier of AdaBoost was SVM, the kernel was RBF and the C and gamma of which were also determined by grid search, ranging from 2–10 to 210 and 0.0001 to 10, respectively. Other parameters were set to default values. In the GBDT model, the number, learning rate, and subsample of base classifiers were also determined by grid search, ranging from 50 to 150, 0 to 1 and 0.5 to 0.8, respectively. The base classifier of GBDT was CART, the max depth and the max leaf nodes of which were also determined by grid, search ranging from 10 to 15 and 10 to 30, respectively. Other parameters were set to default values. In order to reduce the contingency and improve the generalization ability, the five-fold cross-validation method was used to evaluate the performance of the model and select the best model. All of the above algorithms were programmed and realized by sklearn in PyCharm IDE using Python 3.7. The computer system is windows 10 professional, the CPU is Inter Core i7-10700K Processor @3.9 GHz, and the RAM is 32 GB. The final parameters of the model are shown in Table 3.
PLV-based functional connectivity features
Phase locking value is a type of connection characteristic, which quantifies the degree of phase synchronization between the two EEG signals (Aydore et al., 2013; Leguia et al., 2021). The Hilbert transform was first applied to the preprocessed EEG data to calculate the instantaneous amplitude and instantaneous phase for each lead site. The PLV indicator was then calculated using the following formula:
Where N denoted the number of EEG segments per subject, θ(t,n) presented the instantaneous phase difference between different leads of the same segment, exp(jθ(t,n)) represented the complex signal obtained with the help of Euler’s formula using phase, and represented the superimposed value of the complex signals of all segments of a patient, which was averaged to obtain the PLV feature value of a subject.
The PLV feature was then quantized into a value in the range [0,1]. When PLV = 1, the phase difference between the two signals was constant, i.e., perfectly synchronized. When PLV = 0, the phase difference was uniformly distributed over the complex plane unit circle according to time, indicating that there was no synchronization. Between 0 and 1, the signal difference exhibited an “overall convergence” nature, such that as PLV tended to 1, two close signals exhibited better synchronization.
Since it would be more accurate to calculate the instantaneous phase of narrowband signals using the Hilbert transform, the preprocessed EEG segments were divided into four narrow bands according to different frequency ranges, namely delta (δ) (1–4 Hz), θ (4–7), alpha (α) (8–13 Hz), and β (14–30 Hz) bands. The PLVEEG values of these four frequency bands were calculated separately for 200 windows (6 s) of each subject’s 20-min EEG recording. Finally, 200 PLVEEG feature matrices of 19 × 19 in each of the four frequency bands were obtained for each subject and averaged into a single matrix for each frequeny band, so that each subject ended up with a total of four feature matrices for four frequency bands. These PLVEEG feature matrices would be further filtered and sorted characterized (Figure 2).
Figure 2. Mean PLVEEG features in the four frequency bands for the CON and the CI groups of PWEs. PWEs, patients with epilepsy; PLV, phase locking value; CI, cognitive impairment; CON, cognitively normal.
Feature extraction
As shown earlier, 23 clinical features were selected based on the previous studies and contents of available medical records. The EEG records of all subjects were divided into four different frequency bands. For each subject’s 200 6 s segments in any of the frequency bands, 19 leads were paired as two by two, and a 19 × 19 PLVEEG functional connectivity matrix was calculated for each segment’s EEG, excluding duplicate PLVEEG features that made comparisons with the leads themselves, to obtain a total of 171 PLVEEG features for the EEG recordings of a given subject. The PLVEEG features from 80 segments were then averaged. A total of 707 clinical-PLVEEG features, including 684 PLVEEG and 23 clinical features, were obtained in the four frequency bands for each subject. However, it was unknown which features were valid for a particular learning algorithm, and for this reason, we needed to filter all the features to select those that were beneficial to the learning algorithm. Filtering features not only optimized the algorithm to make the model more generalized but also reduced the running time of the algorithm resolving overfitting issues and the difficulty of the learning task, thereby improving the efficiency and the interpretability of the model.
Fisher score is a common feature filtering method (Zhang J. et al., 2021). Features with a strong discriminatory performance exhibit the smallest possible intra-class distance and the largest possible inter-class distance. The higher the inter-class variance and the lower the intra-class variance of PLVEEG features in the same frequency band from different patients, the higher the Fisher score value is. We ranked the features from the largest to the smallest, based on their Fisher score values, with the higher ranked features being theoretically more discriminative.
Modeling process
The classification models were trained using AdaBoost and GBDT platforms as classifiers. Models were constructed based on the pure clinical features, PLVEEG features, and combined clinical- PLVEEG features, as well. To improve the classification performance, generalization skills, and speed of each model, Fisher scores were used to filter the features. Five-fold cross-validation was used to construct the classification model, using 80% of the two sets of data each time, and the remaining 20% of the data was used for model validation.
Statistical analysis
To compare the variability of clinical and normalized PLVEEG features between the CON and CI groups, the quantitative data were first tested for normality using the Shapiro–Wilk test, followed by a comparison of the data with a normal distribution expressed as mean ± standard deviation (SD) using the independent samples t-test, and the Mann–Whitney U-test was applied for data with an abnormal distribution expressed as median ± interquartile range (IQR). For qualitative information, the chi-squared (χ2) test or Fisher’s exact test was used to assess the variability between the two data sets. A p- or p’- value of < 0.05 was considered statistically significant, where p’ referred to a p-value that was corrected by the false discovery rate (FDR) correction. We used SPSS v26.0 for all kinds of statistical analyses.
Results
Clinical feature-based model construction
Of the 23 clinical features, we used Fisher scores to filter the top 15 clinical features in terms of weightage to construct the diagnostic model (Table 4A). The selected features were educational attainment, seizure frequency, VPA, class of antiseizures medications, hippocampal atrophy and sclerosis, age, status epilepticus, MRI of the head for abnormalities, time from the first onset to the current visit, history of previous CNS infections, TLE, anxiety, age at the first onset, history of previous head surgeries or trauma, and gender. The features that showed significant statistical differences between the two groups were educational attainment, seizure frequency, VPA, class of antiseizures medications, hippocampal atrophy and sclerosis, and time from the first onset to the current visit. In the classification model, constructed based on the clinical features using AdaBoost, the model performances after a five-fold cross-validation for accuracy, precision, recall, F1-score, and AUC were 67.89, 66.69, 91.57, 76.71, and 0.75%, respectively. While, in case of the classification model built by GBDT, the final performances after cross-validation for accuracy, precision, recall, F1-score, and AUC were, respectively, 68.09, 70.80, 75.84, 72.62, and 0.76% (Figure 3 and Figure 4A). Therefore, these two algorithms were found to differ slightly in the construction of a model for identifying impaired consciousness in epilepsy patients using the clinical features only.
Figure 3. The evaluation indexes after five-fold cross-validation. GBDT, Gradient Boosting Decision Tree; AdaBoost, Adaptive Boosting.
Figure 4. The performance of six models. (A) Pure clinical features. (B) Pure electroencephalogram (EEG) features. (C) Combined clinical and PLVEEG features. GBDT, Gradient Boosting Decision Tree; AdaBoost, Adaptive Boosting; AUC, area under the curve; ROC, receiver operating-characteristic curve; std. dev, standard deviation.
PLVEEG feature-based model construction
A total of 171 PLVEEG features were extracted for each of the 4 bands of the 20-min EEG recording for each patient, accounting for a total of 684 features (Table 4B). Then the model was constructed using those features with Fisher scores in the top 150 ranks. In the AdaBoost-based classification model, the model performance after a five-fold cross-validation for accuracy, precision, recall, F1-score, and AUC were 83.93, 84.76, 88.08, 86.30, and 0.91%, respectively. Likewise, for the GBDT-based classification model, the final performances after the cross-validation for accuracy, precision, recall, F1-score, and AUC were 88.58, 92.17, 88.17, 90.05, and 0.94%, respectively (Figure 3 and Figure 4B). Importantly, the GBDT was found to outperform AdaBoost in classification model construction using PLVEEG features, demonstrating that the GBDT-based model could be more accurate in identifying epilepsy patients suffering from cognitive dysfunctions. It was also found that PLVEEG features in θ band T5-T6, θ band Fp1-Pz, δ band Fp1-Pz, β band P3-F4, and β band C3-F4 were the top 5 most important ones that might influence the model.
A combined clinical-PLVEEG feature-based model construction
The combined clinical-PLVEEG features were found the most appropriate for constructing the best performing classification models, using either AdaBoost or GBDT algorithm. A total of 707 features were screened using Fisher scores for 23 clinical features and 684 PLVEEG features. A total of 4 clinical features were selected within the top 10 weighted features, namely MRI of the head for abnormalities in the first rank, educational attainment in the second rank, seizure frequency in the fourth rank, and hippocampal atrophy or sclerosis in the sixth rank; all of which were significantly differed between the two groups. Between the two groups, the remaining PLVEEG features with significant differences were C3-F4 in the β-band, Fp1-Pz in the θ-band, F3-F8 in the β-band, C3-P4 in the β-band, C3-P4 in the θ-band, and T5-T6 in the β-band, with only Fp1-Pz in the θ-band, and T5-T6 in the β-band. Although many features were not statistically different between the two groups, they exhibited a very strong impact on the model after the Fisher score screening. Whereas a total of 12 PLVEEG features in the top 5% of features affecting the model performance were significantly different between the two groups, including eight features in the θ band and three PLVEEG features in the β band. We suspected that PLVEEG in the θ band might be the biomarker that could distinguish between these two groups (Table 4C and Figure 5).
Figure 5. In the combined clinical-PLVEEG model, there were statistically significant differences in 12 PLVEEG features between the CON and CI groups of PWEs. The higher the fisher score, the tighter the connection between the leads. P < 0.05 is considered statistically significant. (A) Alpha band; (B) beta band; (C) theta band; (D) delta band. PWEs, patients with epilepsy; PLV, phase locking value; CI, cognitive impairment; CON, cognitively normal.
For AdaBoost, the top 150 Fisher scores were selected to build the classification model, and the final performances after five-fold cross-validation were 87.78, 85.95, 93.17, 89.35, and 0.92% for accuracy, precision, recall, F1-score, and AUC, respectively. While for GBDT, the top 250 Fisher scores were selected to build the classification model, and the model performances after five-fold cross-validation were 90.11, 93.40, 89.50, 91.39, and 0.95% for accuracy, precision, recall, F1-score, and AUC, respectively (Figure 3 and Figure 4C). The recall performance of the AdaBoost model was found to be slightly higher than that of the GDBT, while GDBT outperformed AdaBoost in terms of other metrics.
Comparison between different models
Six models, based on the clinical features only, PLVEEG features only, and combined clinical-PLVEEG features, were constructed for 55 CON and 76 epilepsy patients suffering from cognitive dysfunctions, using the ensemble algorithms like AdaBoost and GBDT. We found that the models constructed with combined clinical-PLVEEG features outperformed those developed with either pure clinical or pure PLVEEG features for both the AdaBoost and GBDT algorithms. Notably, the models constructed solely with clinical features performed the worst. The cross-sectional comparisons also revealed that GBDT-built models outperformed the AdaBoost-based ones in both classification models constructed with PLVEEG features. Furthermore, GBDT also outperformed AdaBoost in cases of both pure clinical features and combined clinical-PLVEEG features, with an exception for recall performance (Table 5).
Not only that, but we could also identify potential biomarkers like EEG indicators using the combined clinical-PLVEEG feature-based models that might be able to detect CI in epilepsy patients, which could be highly useful in the diagnosis of epilepsy in clinical settings. Additionally, many of the clinical features used have also been reported in previous studies suggesting their strong association with CI symptoms in epilepsy patients, but have not been ranked to the extent to which these clinical features might affect cognition. Therefore, we ranked these clinical features by their respective Fisher scores. Our findings suggest that EEG could be of great interest to subjects with cognitive deficits, especially those with epileptic symptoms. Previously, technical limitations were the main obstacle in improving the application of EEG for epilepsy diagnosis and treatment. By estimating the combined effects of clinical and PLVEEG features, we could predict the current cognitive status in epilepsy patients, providing clinicians with more options for precise diagnosis and effective treatment plans.
Discussion
To the best of our knowledge, the present study is the first of its kind to use an integrated algorithm for the construction of a classification model for facilitating the diagnosis of CI in PWE by combined clinical and PLVEEG functional connectivity features.
Advantages of combined clinical-PLVEEG features for classification model building
Although several risk factors affecting cognitive functions in epilepsy have been identified, however, only a few studies have used these clinical features to predict whether PWEs have a comorbid CI situation. Importantly, it’s been difficult to determine the extent to which these clinical features might affect cognition with a background of epilepsy. A meta-analysis (Novak et al., 2022a) has found that duration of epilepsy, frequency of seizures, and use of antiseizures medications are important clinical features that can affect cognition. Moreover, some studies suggest that education, history of surgical head trauma, anxiety and depression, hippocampal abnormalities, TLE, and seizure types may influence cognitive functions in PWEs (Piazzini et al., 2006; Bell et al., 2011; Vrinda et al., 2019; Jarcuskova et al., 2020; Wang et al., 2020; Phuong et al., 2021; Elsherif and Esmael, 2022). A previous study (Lin et al., 2021) collected 12 clinical features from outpatients with epilepsy to construct a model for diagnosing CI with a performance accuracy, recall, precision, and AUC of 60, 51, 88, and 0.71%, respectively, and concluded that status epilepticus, history of previous surgical head trauma, and seizure frequency were the top three clinical features affecting cognition. However, the clinical features considered in this study were not comprehensive enough, for example, it did not take into account important factors affecting PWEs such as education level and the classes of antiseizures medications taken (Wang et al., 2020). It was previously thought that VPA, PHT, and TPM could cause cognitive dysfunctions in PWEs (Brunbech and Sabers, 2002; Dang et al., 2021; Lozano-Garcia et al., 2021), and for this reason, the presence or absence of these three drugs was used as a clinical feature. The study showed that only VPA had significant weightage for this model, while PHT and TPM, probably due to insufficient data, were not statistically significant, and did not contribute to the construction of the model.
Of the models constructed using pure clinical features, the performance accuracy, recall, precision, and AUC for the AdaBoost/GBDT models were 67.89/68.03%, 91.57/75.84%, 66.69/70.80%, and 0.75/0.76%, respectively. Using Fisher scores, we selected 23 clinical features. Of these, education level, seizure frequency, and VPA therapy ranked the top three clinical characteristics affecting cognition in PWEs. Among the models constructed with combined clinical and PLVEEG features, the accuracy, recall, precision, and AUC of the AdaBoost/GBDT models were 87.78/90.11%, 93.17/89.50%, 85.95/93.40%, and 0.92/0.95%, respectively. We applied the Fisher scoring method for the 23 clinical features and 684 PLVEEG features to jointly screen and rank. Among these features, MRI abnormalities, education level, and seizure frequency were the top 3 most influential clinical features. The performance of the models constructed using clinical features alone was better than that shown in previous studies for all metrics, except for the performance accuracy. While the performance of the models constructed using combined clinical and PLVEEG features was significantly improved than that reported previously. Thus, we concluded that combined clinical and PLVEEG features were more appropriate for PWEs and that a combination of different types of features would be an optimal choice for constructing diagnostic prediction models.
PLVEEG features are valid indicators for diagnosing CI in PWEs
PLVEEG is used to remotely examine the task-induced changes in neural activities, synchronized in EEG recordings, which is a classic metric for computing functional brain connectivity features. Jones et al. (2022) have used PLVEEG functional connectivity features as an evaluation metric for assessing the efficacy of transcranial alternating current stimulation (tACS) on age-associated cognitive decline. Li et al. (2022) have constructed a model combining the clinical and PLVEEG features to diagnose Alzheimer’s disease (AD), which exhibits satisfactory performance and robustness. Another study (Lanzone et al., 2021) has found that PLVEEG in the α band of patients who were effective on treatment with perampanel as an add-on drug could be used as a biomarker to predict the responsiveness to perampanel drugs. Cho et al. (2017) have reported that PLVEEG in the γ band may be a potential biomarker for predicting seizures. In this study, the accuracy, recall, precision, and AUC of the AdaBoost/GBDT models were 83.93/88.58%, 84.76/92.17%, 88.08/88.17%, 86.30/90.05%, and 0.91/0.94%, respectively, when only the PLVEEG features were used for the model construction. The θ-band T5- T6, θ-band Fp1-Pz, and δ-band Fp1-Pz were the top three PLVEEG features affecting the model weightage, indicating that the PLVEEG functional connectivity features might be valid indicators for the diagnosis of cognitive dysfunctions comorbid with epilepsy.
PLVEEG features in the θ band may be a potential biomarker for diagnosing CI in PWEs
Here, we calculated the PLVEEG features of the four frequency bands (α, β, θ, δ), and found that the PLVEEG features, especially of the θ band, might be potential biomarkers to distinguish between epilepsy patients with or without comorbid CI. In our constructed model of the combined clinical and PLVEEG features, we employed Fisher scoring to rank individual features, which revealed 12 PLV features that ranked in the top 30 were significantly different between the CON and CI groups. Notably, eight of these features were related to the θ band and three to the β band.
The θ band has been found to have an important relationship with epilepsy and cognitive function in previous studies. One study (Douw et al., 2010) has demonstrated that functional connectivity features in the θ band could be used to aid in the diagnosis of epilepsy with a recall of 62% and a specificity of 72%. Other studies (Jun et al., 2020) have also suggested that stimulation of the hippocampus may increase the release of θ rhythms, thereby improving the associative memory function. These studies suggest that increasing the θ rhythm in the hippocampus may provide a theoretical basis for the neural mechanisms of memory enhancement. Moreover, Gupta et al. (2012) have identified that θ rhythms in the hippocampus of rats are associated with visuospatial abilities and executive abilities related to memory and cognition. Another study (Braithwaite et al., 2020) has revealed that increased power of the θ rhythm in children can be a valid biomarker for predicting non-verbal cognitive abilities. Furthermore, it (Ahmadlou et al., 2014) has been concluded that functional connectivity features in the θ band could be used to differentiate between patients with mild CI and healthy elderly populations. Briels et al. (2020) have found that functional connectivity indicators in the θ and β frequency bands in AD patients may help diagnose the disease severity. Other studies (Singh et al., 2018) have shown that a reduction in midfrontal θ wave frequency responds to the degree of effective control of cognitive functions in PD patients. The θ rhythms in the frontal lobe are highly correlated with cognitive function (Cavanagh and Frank, 2014), with Fp1-Fz being within the frontal lobe. Our results showed that the PLVEEG features of Fp1-Fz in the θ band were significantly different between the CON and CI groups of epilepsy patients, accounting for a high weightage in the diagnostic model. In this context, one study (Cao et al., 2022) has reported an important relationship between the θ rhythm and cognition in patients with schizophrenia, indicating that superior cognitive performance may be significantly associated with a smaller θ wave power, and altered θ rhythm and cognition are highly correlated mainly in the parieto-occipital lobe. The P4 and T5 were close to the occipital region in our investigation. The PLVEEG for P4-T5 were also significantly different between the two groups and accounted for a higher weightage in the model. Furthermore, it is shown (Usami et al., 2019) that β oscillations can enhance the responsiveness of the cerebral cortex to inputs from distant cortices, suggesting that β frequencies may have an important role in functional connectivity. Interestingly, α frequency is significantly increased in AD patients presenting with mild cognitive dysfunctions (Moretti, 2015). The α frequency was found to be less influential in our study, in terms of statistical significance and the weightage of the model, possibly due to the exclusion of AD patients’ data.
Previous studies have amply demonstrated the significance of functional connectivity features in the θ band in the diagnosis of epilepsy and cognitive dysfunctions. Therefore, our study demonstrated that PLVEEG features in the θ band might be reliable biomarkers for diagnosing CI in PWEs, especially those with high Fisher scores.
Limitations and future directions
Despite these excellent results, there are still certain limitations to this study. First, this was a single-center retrospective study with data from only one institutional epilepsy center and a small sample population. Although the combined clinical and PLVEEG features and advanced algorithms ensured the accuracy of our results, multi-center prospective studies are warranted for the generalization of our results. Here, we provided a theoretical basis and demonstrated the possibilities of further improving the diagnostic methods for PWEs comorbid with CI. Second, this study was based on the MoCA scale. However, we classified the features based on the total MoCA scores rather than the subtest scores. Although our model could address the issue of differentiating PWEs with or without cognitive deficits, the content of each subtest should be investigated more carefully in the future. Finally, the potential biomarkers that we extracted were mainly functional connectivity features of the EEG and a subset of clinical features. The future brain network features extracted from MRI examinations can be useful in improving the accuracy and superiority of the model. We propose to validate the performance of our models with larger datasets from multiple epilepsy centers in the future, as well as add new features to improve the accuracy of the model.
Conclusion
In this study, we constructed a diagnostic model for CI in PWEs based on the combined clinical and PLVEEG features. Besides, we found that PLVEEG functional connectivity features in the θ band might be potential biomarkers for the diagnosis of CI in PWEs.
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 Institutional Review Board of the Henan Provincial People’s Hospital. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author contributions
ZR and MY designed the study. XH obtained funding. ZR, YZ, YH, YXH, and QW acquired the data. PZ, TZ, and NW analyzed EEG recordings. ZR, ZZ, and BWe worked on EEG preprocessing and machine learning process. ZR, ZZ, and BWa conducted the statistical analysis. ZR, ZZ, YH, YXH, and QW analyzed and interpreted the data. ZR and XH drafted and revised the manuscript. All authors revised this draft and read and approved the final manuscript.
Funding
This study was sponsored by Henan Province’s Gong Jian Program (Authorization number: SB201901074) and 23456 Talent Engineering (Authorization number: ZC20200371).
Acknowledgments
We would like to thank the staff members at the Department of Neurology at the Henan Provincial People’s Hospital for their assistance with this research study and we would like to express our most profound gratitude to our patients and their families.
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.
Abbreviations
PWEs, patients with epilepsy; MoCA, Montreal Cognitive Assessment; CI, cognitive impairment; CON, cognitively normal; EEG, electroencephalogram; PLV, phase locking value; AdaBoost, Adaptive Boosting; GBDT, Gradient Boosting Decision Tree; MRI, magnetic resonance imaging; AUC, area under the curve.
References
Ahmadlou, M., Adeli, A., Bajo, R., and Adeli, H. (2014). Complexity of functional connectivity networks in mild cognitive impairment subjects during a working memory task. Clin. Neurophysiol. 125, 694–702. doi: 10.1016/j.clinph.2013.08.033
Aydore, S., Pantazis, D., and Leahy, R. M. A. (2013). Note on the Phase locking value and its properties. Neuroimage 74, 231–244. doi: 10.1016/j.neuroimage.2013.02.008
Balcik, Z. E., Senadim, S., Tekin, B., Ceyhan Dirican, A., Eren, F., Karahan, M. G., et al. (2020). Do interictal EEG findings reflect cognitive function in juvenile myoclonic epilepsy? Epilepsy Behav. 111:107281. doi: 10.1016/j.yebeh.2020.107281
Bell, B., Lin, J. J., Seidenberg, M., and Hermann, B. (2011). The neurobiology of cognitive disorders in temporal lobe epilepsy. Nat. Rev. Neurol. 7, 154–164. doi: 10.1038/nrneurol.2011.3
Bergeron, D., Flynn, K., Verret, L., Poulin, S., Bouchard, R. W., Bocti, C., et al. (2017). Multicenter validation of an MMSE-MoCA conversion table. J. Am. Geriatr. Soc. 65, 1067–1072. doi: 10.1111/jgs.14779
Black, L. C., Schefft, B. K., Howe, S. R., Szaflarski, J. P., Yeh, H. S., and Privitera, M. D. (2010). The effect of seizures on working memory and executive functioning performance. Epilepsy Behav. 17, 412–419. doi: 10.1016/j.yebeh.2010.01.006
Braithwaite, E. K., Jones, E. J. H., Johnson, M. H., and Holmboe, K. (2020). Dynamic modulation of frontal theta power predicts cognitive ability in infancy. Dev. Cogn. Neurosci. 45:100818. doi: 10.1016/j.dcn.2020.100818
Briels, C. T., Schoonhoven, D. N., Stam, C. J., de Waal, H., Scheltens, P., and Gouw, A. A. (2020). Reproducibility of EEG functional connectivity in Alzheimer’s disease. Alzheimers Res. Ther. 12:68. doi: 10.1186/s13195-020-00632-3
Brunbech, L., and Sabers, A. (2002). Effect of antiepileptic drugs on cognitive function in individuals with epilepsy: A comparative review of newer versus older agents. Drugs 62, 593–604. doi: 10.2165/00003495-200262040-00004
Cai, M., Dang, G., Su, X., Zhu, L., Shi, X., Che, S., et al. (2021). Identifying mild cognitive impairment in Parkinson’s disease with electroencephalogram functional connectivity. Front. Aging Neurosci. 13:701499. doi: 10.3389/fnagi.2021.701499
Cao, Y., Han, C., Peng, X., Su, Z., Liu, G., Xie, Y., et al. (2022). Correlation between resting theta power and cognitive performance in patients with schizophrenia. Front. Hum. Neurosci. 16:853994. doi: 10.3389/fnhum.2022.853994
Cavanagh, J. F., and Frank, M. J. (2014). Frontal theta as a mechanism for cognitive control. Trends Cogn. Sci. 18, 414–421. doi: 10.1016/j.tics.2014.04.012
Cho, D., Min, B., Kim, J., and Lee, B. (2017). EEG-based prediction of epileptic seizures using phase synchronization elicited from noise-assisted multivariate empirical mode decomposition. IEEE Trans. Neural Syst. Rehabil. Eng. 25, 1309–1318. doi: 10.1109/TNSRE.2016.2618937
Dang, Y. L., Foster, E., Lloyd, M., Rayner, G., Rychkova, M., Ali, R., et al. (2021). Adverse events related to antiepileptic drugs. Epilepsy Behav. 115:107657. doi: 10.1016/j.yebeh.2020.107657
Del Brutto, O. H., Mera, R. M., Del Brutto, V. J., Zambrano, M., Wright, C. B., and Rundek, T. (2019). Clinical and neuroimaging risk factors for cognitive decline in community-dwelling older adults living in rural ecuador. A population-based prospective cohort study. Int. J. Geriatr. Psychiatry 34, 447–452. doi: 10.1002/gps.5037
Delorme, A., and Makeig, S. (2004). EEGLAB: An open source toolbox for analysis of single-trial EEG dynamics including independent component analysis. J. Neurosci. Methods 134, 9–21. doi: 10.1016/j.jneumeth.2003.10.009
Douw, L., de Groot, M., van Dellen, E., Heimans, J. J., Ronner, H. E., Stam, C. J., et al. (2010). ’Functional connectivity’ is a sensitive predictor of epilepsy diagnosis after the first seizure. PLoS One 5:e10839. doi: 10.1371/journal.pone.0010839
Duma, G. M., Danieli, A., Mattar, M. G., Baggio, M., Vettorel, A., Bonanni, P., et al. (2022). Resting state network dynamic reconfiguration and neuropsychological functioning in temporal lobe epilepsy: An HD-EEG investigation. Cortex 157, 1–13. doi: 10.1016/j.cortex.2022.08.010
Duma, G. M., Danieli, A., Vettorel, A., Antoniazzi, L., Mento, G., and Bonanni, P. (2021). Investigation of dynamic functional connectivity of the source reconstructed epileptiform discharges in focal epilepsy: A graph theory approach. Epilepsy Res. 176:106745. doi: 10.1016/j.eplepsyres.2021.106745
Edeh, M. O., Dalal, S., Dhaou, I. B., Agubosim, C. C., Umoke, C. C., Richard-Nnabu, N. E., et al. (2022). Artificial intelligence-based ensemble learning model for prediction of hepatitis C disease. Front. Public Health 10:892371. doi: 10.3389/fpubh.2022.892371
Elahian, B., Yeasin, M., Mudigoudar, B., Wheless, J. W., and Babajani-Feremi, A. (2017). Identifying seizure onset zone from electrocorticographic recordings: A machine learning approach based on phase locking value. Seizure 51, 35–42. doi: 10.1016/j.seizure.2017.07.010
Elsherif, M., and Esmael, A. (2022). Hippocampal atrophy and quantitative EEG markers in mild cognitive impairment in temporal lobe epilepsy versus extra-temporal lobe epilepsy. Neurol. Sci. 43, 1975–1986. doi: 10.1007/s10072-021-05540-4
Fadaie, F., Lee, H. M., Caldairou, B., Gill, R. S., Sziklas, V., Crane, J., et al. (2021). Atypical functional connectome hierarchy impacts cognition in temporal lobe epilepsy. Epilepsia 62, 2589–2603. doi: 10.1111/epi.17032
Fisher, R. S., Acevedo, C., Arzimanoglou, A., Bogacz, A., Cross, J. H., Elger, C. E., et al. (2014). Ilae official report: A practical clinical definition of epilepsy. Epilepsia 55, 475–482. doi: 10.1111/epi.12550
Gavrilovic, A., Toncev, G., Boskovic Matic, T., Vesic, K., Ilic Zivojinovic, J., and Gavrilovic, J. (2019). Impact of epilepsy duration, seizure control and EEG abnormalities on cognitive impairment in drug-resistant epilepsy patients. Acta Neurol. Belg. 119, 403–410. doi: 10.1007/s13760-019-01090-x
Gupta, A. S., van der Meer, M. A., Touretzky, D. S., and Redish, A. D. (2012). Segmentation of spatial experience by hippocampal theta sequences. Nat. Neurosci. 15, 1032–1039. doi: 10.1038/nn.3138
He, X., Bassett, D. S., Chaitanya, G., Sperling, M. R., Kozlowski, L., and Tracy, J. I. (2018). Disrupted dynamic network reconfiguration of the language system in temporal lobe epilepsy. Brain 141, 1375–1389. doi: 10.1093/brain/awy042
Helmstaedter, C., and Witt, J. A. (2017). Epilepsy and cognition–a bidirectional relationship? Seizure 49, 83–89. doi: 10.1016/j.seizure.2017.02.017
Huang, H., Cui, G., Tang, H., Kong, L., Wang, X., Cui, C., et al. (2022). Relationships between plasma expression levels of microrna-146a and microrna-132 in epileptic patients and their cognitive, mental and psychological disorders. Bioengineered 13, 941–949. doi: 10.1080/21655979.2021.2015528
Jarcuskova, D., Palusna, M., Gazda, J., Feketeova, E., and Gdovinova, Z. (2020). Which clinical and neuropsychological factors are responsible for cognitive impairment in patients with epilepsy? Int. J. Public Health 65, 947–956. doi: 10.1007/s00038-020-01401-7
Jones, K. T., Johnson, E. L., Gazzaley, A., and Zanto, T. P. (2022). Structural and functional network mechanisms of rescuing cognitive control in aging. Neuroimage 262:119547. doi: 10.1016/j.neuroimage.2022.119547
Jun, S., Lee, S. A., Kim, J. S., Jeong, W., and Chung, C. K. (2020). Task-dependent effects of intracranial hippocampal stimulation on human memory and hippocampal theta power. Brain Stimul. 13, 603–613. doi: 10.1016/j.brs.2020.01.013
Lanzone, J., Ricci, L., Tombini, M., Boscarino, M., Mecarelli, O., Pulitano, P., et al. (2021). The effect of perampanel on EEG spectral power and connectivity in patients with focal epilepsy. Clin. Neurophysiol. 132, 2176–2183. doi: 10.1016/j.clinph.2021.05.026
Leguia, M. G., Andrzejak, R. G., Rummel, C., Fan, J. M., Mirro, E. A., Tcheng, T. K., et al. (2021). Seizure cycles in focal epilepsy. JAMA Neurol. 78, 454–463. doi: 10.1001/jamaneurol.2020.5370
Li, X., Zhou, T., and Qiu, S. (2022). Alzheimer’s disease analysis algorithm based on no-threshold recurrence plot convolution network. Front. Aging Neurosci. 14:888577. doi: 10.3389/fnagi.2022.888577
Lin, F., Han, J., Xue, T., Lin, J., Chen, S., Zhu, C., et al. (2021). Predicting cognitive impairment in outpatients with epilepsy using machine learning techniques. Sci. Rep. 11:20002. doi: 10.1038/s41598-021-99506-3
Lozano-Garcia, A., Hampel, K. G., Villanueva, V., Gonzalez-Bono, E., and Cano-Lopez, I. (2021). The number of anti-seizure medications mediates the relationship between cognitive performance and quality of life in temporal lobe epilepsy. Epilepsy Behav. 115:107699. doi: 10.1016/j.yebeh.2020.107699
Ma, J., Wang, Z., Cheng, T., Hu, Y., Qin, X., Wang, W., et al. (2022). A prediction model integrating synchronization biomarkers and clinical features to identify responders to vagus nerve stimulation among pediatric patients with drug-resistant epilepsy. CNS Neurosci. Ther. 28, 1838–1848. doi: 10.1111/cns.13923
Melikyan, Z. A., Malek-Ahmadi, M., O’Connor, K., Atri, A., Kawas, C. H., and Corrada, M. M. (2021). Norms and equivalences for MoCA-30, MoCA-22, and MMSE in the oldest-old. Aging Clin. Exp. Res. 33, 3303–3311. doi: 10.1007/s40520-021-01886-z
Montano-Lozada, J. M., Lopez, N., Espejo-Zapata, L. M., Soto-Anari, M., Ramos-Henderson, M., Caldichoury-Obando, N., et al. (2021). Cognitive changes in patients with epilepsy identified through the moca test during neurology outpatient consultation. Epilepsy Behav. 122:108158. doi: 10.1016/j.yebeh.2021.108158
Moretti, D. V. (2015). Theta and alpha EEG frequency interplay in subjects with mild cognitive impairment: Evidence from EEG, MRI, and SPECT brain modifications. Front. Aging Neurosci. 7:31. doi: 10.3389/fnagi.2015.00031
Novak, A., Vizjak, K., and Rakusa, M. (2022a). Cognitive impairment in people with epilepsy. J. Clin. Med. 11:267. doi: 10.3390/jcm11010267
Novak, A., Vizjak, K., Gacnik, A., and Rakusa, M. (2022b). Cognitive impairment in people with epilepsy: Montreal cognitive assessment (MoCA) as a screening tool. Acta. Neurol. Belg. 8:4. doi: 10.1007/s13760-022-02046-4
Peng, T., Chen, X., Wan, M., Jin, L., Wang, X., Du, X., et al. (2020). The prediction of hepatitis E through ensemble learning. Int. J. Environ. Res. Public Health 18:159. doi: 10.3390/ijerph18010159
Phuong, T. H., Houot, M., Mere, M., Denos, M., Samson, S., and Dupont, S. (2021). Cognitive impairment in temporal lobe epilepsy: Contributions of lesion, localization and lateralization. J. Neurol. 268, 1443–1452. doi: 10.1007/s00415-020-10307-6
Piazzini, A., Canevini, M. P., Turner, K., Chifari, R., and Canger, R. (2006). Elderly people and epilepsy: Cognitive function. Epilepsia 47(Suppl. 5) 82–84. doi: 10.1111/j.1528-1167.2006.00884.x
Rodrigues, S. G., Gouveia, R. G., and Bentes, C. (2020). Moca as a cognitive assessment tool for absence status epilepticus. Epileptic Disord. 22, 229–232. doi: 10.1684/epd.2020.1149
Singh, A., Richardson, S. P., Narayanan, N., and Cavanagh, J. F. (2018). Mid-frontal theta activity is diminished during cognitive control in Parkinson’s disease. Neuropsychologia 117, 113–122. doi: 10.1016/j.neuropsychologia.2018.05.020
Titiz, A. S., Mahoney, J. M., Testorf, M. E., Holmes, G. L., and Scott, R. C. (2014). Cognitive impairment in temporal lobe epilepsy: Role of online and offline processing of single cell information. Hippocampus 24, 1129–1145. doi: 10.1002/hipo.22297
Ung, H., Cazares, C., Nanivadekar, A., Kini, L., Wagenaar, J., Becker, D., et al. (2017). Interictal epileptiform activity outside the seizure onset zone impacts cognition. Brain 140, 2157–2168. doi: 10.1093/brain/awx143
Usami, K., Milsap, G. W., Korzeniewska, A., Collard, M. J., Wang, Y., Lesser, R. P., et al. (2019). Cortical responses to input from distant areas are modulated by local spontaneous alpha/beta oscillations. Cereb. Cortex 29, 777–787. doi: 10.1093/cercor/bhx361
Vrinda, M., Arun, S., Srikumar, B. N., Kutty, B. M., and Shankaranarayana Rao, B. S. (2019). Temporal lobe epilepsy-induced neurodegeneration and cognitive deficits: Implications for aging. J. Chem. Neuroanat. 95, 146–153. doi: 10.1016/j.jchemneu.2018.02.005
Wang, L., Chen, S., Liu, C., Lin, W., and Huang, H. (2020). Factors for cognitive impairment in adult epileptic patients. Brain Behav. 10:e01475. doi: 10.1002/brb3.1475
Wenbo, W., Yang, S., and Guici, C. (2021). Blood glucose concentration prediction based on VMD-KELM-adaboost. Med. Biol. Eng. Comput. 59, 2219–2235. doi: 10.1007/s11517-021-02430-x
Zhang, J., Xu, D., Hao, K., Zhang, Y., Chen, W., Liu, J., et al. (2021). FS-GBDT: Identification multicancer-risk module via a feature selection algorithm by integrating fisher score and Gbdt. Brief Bioinform. 22:bbaa189. doi: 10.1093/bib/bbaa189
Keywords: epilepsy, cognitive impairment, EEG, phase locking value, GBDT, AdaBoost, diagnostic model, Fisher score
Citation: Ren Z, Zhao Y, Han X, Yue M, Wang B, Zhao Z, Wen B, Hong Y, Wang Q, Hong Y, Zhao T, Wang N and Zhao P (2023) An objective model for diagnosing comorbid cognitive impairment in patients with epilepsy based on the clinical-EEG functional connectivity features. Front. Neurosci. 16:1060814. doi: 10.3389/fnins.2022.1060814
Received: 15 October 2022; Accepted: 28 December 2022;
Published: 12 January 2023.
Edited by:
Reza Rastmanesh, The Nutrition Society, United KingdomReviewed by:
Gian Marco Duma, Eugenio Medea (IRCCS), ItalyAleksandar Gavrilovic, University of Kragujevac, Serbia
Copyright © 2023 Ren, Zhao, Han, Yue, Wang, Zhao, Wen, Hong, Wang, Hong, Zhao, Wang and Zhao. 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: Xiong Han, aGFueGlvbmdAenp1LmVkdS5jbg==
†These authors have contributed equally to this work
‡These authors share first authorship