Miseon Shim
Chang-Hwan Im
Seung-Hwan Lee
Han-Jeong Hwang- 1Department of Electronics and Information, Korea University, Sejong, South Korea
- 2Industry Development Institute, Korea University, Sejong, South Korea
- 3Department of Biomedical Engineering, Hanyang University, Seoul, South Korea
- 4Department of Psychiatry, Ilsan Paik Hospital, Inje University, Goyang, South Korea
- 5Clinical Emotion and Cognition Research Laboratory, Goyang, South Korea
- 6Interdisciplinary Graduate Program for Artificial Intelligence Smart Convergence Technology, Korea University, Sejong, South Korea
Electroencephalography (EEG)-based diagnosis of psychiatric diseases using machine-learning approaches has made possible the objective diagnosis of various psychiatric diseases. The objective of this study was to improve the performance of a resting-state EEG-based computer-aided diagnosis (CAD) system to diagnose post-traumatic stress disorder (PTSD), by optimizing the frequency bands used to extract EEG features. We used eyes-closed resting-state EEG data recorded from 77 PTSD patients and 58 healthy controls (HC). Source-level power spectrum densities (PSDs) of the resting-state EEG data were extracted from 6 frequency bands (delta, theta, alpha, low-beta, high-beta, and gamma), and the PSD features of each frequency band and their combinations were independently used to discriminate PTSD and HC. The classification performance was evaluated using support vector machine with leave-one-out cross validation. The PSD features extracted from slower-frequency bands (delta and theta) showed significantly higher classification performance than those of relatively higher-frequency bands. The best classification performance was achieved when using delta PSD features (86.61%), which was significantly higher than that reported in a recent study by about 13%. The PSD features selected to obtain better classification performances could be explained from a neurophysiological point of view, demonstrating the promising potential to develop a clinically reliable EEG-based CAD system for PTSD diagnosis.
Introduction
Post-traumatic stress disorder (PTSD) is a psychiatric disorder caused by experiencing or witnessing traumatic events (Yehuda, 2002), and PTSD patients are diagnosed through interview with clinical experts based on the Diagnostic and Statistical Manual of Mental Disorders, 5th Edition (DSM-5, note that all abbreviations were summarized in Supplementary Table 1) (Friedman et al., 2011). Post-traumatic stress disorder patients generally show a high rate of comorbidity with other mental illnesses, which leads to confusion of diagnosis (Ginzburg et al., 2010; Elhai et al., 2011; Gros et al., 2012). According to a previous study, 73.3% of veteran PTSD patients had comorbid anxiety disorders (e.g., general, panic, and social anxiety disorder) (Magruder et al., 2005). Another study reported that 68% of veteran PTSD patients met the criteria for major depressive disorder (MDD) (Grubaugh et al., 2010). These results imply that when interviewing with clinical experts, if PTSD patients hide their traumatic histories or symptoms, they could be misdiagnosed as other mental diseases with high probability. Thus, it is necessary to introduce a PTSD diagnosis tool to compliment the diagnostic failure rate of traditional diagnosis (Sumpter and McMillan, 2005).
In recent years, the machine-learning-based computer-aided diagnosis (CAD) system has received increased attention, due to its ability to predict the state of neuropsychiatric diseases using objective neurophysiological biomarkers (McBride et al., 2011; Mitra et al., 2016; Vergara et al., 2017, 2018). Early CAD systems developed for those suffering from traumatic events have used neuroimaging-based features for its diagnosis, but they have focused on patients with traumatic brain injury (TBI), rather than PTSD (McBride et al., 2011; Mitra et al., 2016; Vergara et al., 2017, 2018). For example, Mitra et al., achieved 68% classification accuracy when differentiating TBI from healthy controls (HC) using diffusion tensor imaging (DTI) features (Mitra et al., 2016), meanwhile, Vergara et al., reported a classification accuracy of 84% when differentiating TBI patients from HC using DTI and functional magnetic resonance imaging (fMRI)-based features (Vergara et al., 2017). However, although both TBI and PTSD are developed by traumatic events, TBI and PTSD should be studied independently because the cause and characteristics of TBI and PTSD are totally different, in that TBI patients are troubled with physical brain damage, while PTSD patients suffer from mental problems (Bryant, 2011).
A few studies have attempted to differentiate PTSD patients from HC using neuroimaging-based features. Rangaprakash et al. (2017) achieved 81% classification accuracy when classifying PTSD patients and HC using effective connectivity network features extracted from resting-state fMRI data. In addition, Zhang et al., introduced a classification model for PTSD diagnosis using magnetoencephalographic (MEG) connectomes, and reported the performance of an area-under-the-curve (AUC) value of 0.9 (Zhang et al., 2020). However, while fMRI and MEG have shown the potential to differentiate PTSD patients from HC, due to their high cost and low portability, they are not usable in practice for both clinicians and patients (Wienbruch et al., 2006). To overcome the mentioned limitations, electroencephalogram (EEG) could be an adequate alternative neuroimaging tool for the diagnosis of PTSD patients (Chen, 2001).
Electroencephalogram (EEG) is more portable than other neuroimaging tools, such as fMRI and MEG (Chen, 2001), and is also suitable to investigate dynamic neuronal changes due to its high-temporal resolution (Saletu et al., 1991; Burle et al., 2015). In particular, since abnormal neuronal changes in psychiatric patients reflect their pathophysiology (Newson and Thiagarajan, 2019), they could be used for the diagnosis of various psychiatric diseases (Orru et al., 2012). Therefore, many researchers have introduced EEG-based CAD systems for the diagnosis of various psychiatric disorders, and achieved promising classification accuracies when differentiating psychiatric disorders from HC (Orru et al., 2012). Recently, we introduced 2 EEG-based CAD systems to assist the accurate diagnosis of PTSD patients (Shim et al., 2019; Kim Y. W. et al., 2020).
In our previous study, we attempted to classify PTSD patients and HC using P300 event-related potential (ERP) based on machine-learning technique (Shim et al., 2019), and obtained a classification accuracy of 80% in differentiating 2 groups. Although we achieved acceptable classification accuracy, the previous study was limited in terms of its usability; PTSD patients were required to perform an auditory attention task to evoke P300 activation even though they generally have difficulty in concentrating on an attention-related task. To overcome this limitation, resting-state EEG could be an appropriate alternative because no effort is required to record resting-state EEG from PTSD patients; also, resting-state EEG reflects the pathophysiological traits of PTSD patients (McFarlane et al., 2005; Veltmeyer et al., 2006). Recently, we investigated the possibility of using resting-state EEG to distinguish PTSD patients from HC, where we employed 2 types of source-level features: (i) power spectrum densities (PSDs) and (ii) network indices based on graph theory (Kim Y. W. et al., 2020). We confirmed that PSD features showed significantly higher classification accuracy than network features, and obtained a maximum classification accuracy of 73.09%.
The objective of this study was to enhance the performance of classifying PTSD patients from HC by optimizing the frequency band used for extracting PSD features from resting-state EEG. Our previous study (Kim Y. W. et al., 2020) did not consider a delta-frequency band of 1–4 Hz even though slow EEG waves were closely related to the typical endophenotypes of PTSD patients (McFarlane et al., 2005; Veltmeyer et al., 2006; Newson and Thiagarajan, 2019), and an optimization of frequency bands was also not performed using traditional 6 frequency bands (delta, theta, alpha, low-beta, high-beta and gamma) in terms of classification performance. Therefore, in this study, we aimed to optimize PSD frequency bands by introducing the delta band and investigating all possible combinations of the 6 frequency bands in order to increase the classification performance between PTSD patients and HC. Furthermore, we investigated the neurophysiological meanings of most discriminable features selected to obtain the best classification performance, thereby contributing to the development of a reliable CAD system to assist the diagnosis of PTSD patients. The main highlights of our contributions are briefly listed below:
i) Optimization of resting-state EEG-based features for assisting diagnosis of PTSD patients.
ii) Significant improvement of classification performance as compared to a recent study by about 13%.
iii) Interpretation of useful features selected to attain the best classification performance from a neurophysiological point of view, thereby providing the basis to develop a reliable CAD system for PTSD diagnosis.
Materials and Methods
Participants
Seventy-seven PTSD patients aged between 20 and 60 years and 58 HC aged between 23 and 60 years were recruited for this study from the Psychiatry Department of Inje University Ilsan Paik Hospital (see Table 1 for detailed demographic data). The patients’ diagnosis by a board-certified psychiatrist was based on the Diagnostic and Statistical Manual of Mental Disorders, 4th edition (DSM-IV) Axis I Psychiatric Disorders. Patients were excluded if they accorded with the following criteria: (1) abnormality of the central nervous system, (2) medical histories of alcohol or drug abuse, (3) intellectual disability, (4) history of head injuries with loss of consciousness and experience with electrical therapy (e.g., electroconvulsive therapy, ECT), and (5) psychotic symptoms lasting for at least 24 h. HC was recruited from the local community through local newspapers and posters. Individuals without any psychiatric medical history were recruited for HC. If the HC was taking or had taken any kinds of psychotropic medication, they were excluded from the study. The study protocol was approved by the Institutional Review Board of Inje University Ilsan Paik Hospital (2015-09-018), and all participants provided written informed consent.
Table 1. Demographic data of post-traumatic stress disorder (PTSD) patients and healthy controls (HC). The p-values represent significant differences between the two groups.
Three psychiatric symptoms of PTSD patients were evaluated by clinical experts. Impact of Event Scale-Revised (IES-R) (Weiss, 2007) was used to determine whether the patients had PTSD or not by evaluating the response severity of traumatic events. Beck Anxiety Inventory (BAI) (Beck and Steer, 1990) and Beck Depression Inventory (BDI) (Beck et al., 1996) were used to check anxiety symptom and depression symptom, respectively (Table 1). Three types of psychotropic medications were prescribed to the patients based on patients’ clinical symptoms examined by clinical experts: antidepressants for depressive symptoms [selective serotonin reuptake inhibitors (n = 67), venlafaxine (n = 10)], antipsychotics for psychotic symptoms [aripiprazole (n = 5), quetiapine (n = 17)], and sedative-hypnotics for anxiety symptoms [lorazepam (n = 37), clonazepam (n = 27), diazepam (n = 15), and alprazolam (n = 35)]. All patients received combined psychotropic medications as follows: (1) antidepressant, antipsychotics, and sedative-hypnotics (n = 32); (2) antidepressant and antipsychotics (n = 12); (3) antidepressant and sedative-hypnotics (n = 33).
Electroencephalogram Recording and Preprocessing
Resting-state EEG data used in this study were the same as presented in our previous study (Shim et al., 2017), where we only investigated disrupted brain networks and the relationships between network indices and symptoms of PTSD patients. Resting-state EEG data were recorded with a band-pass filter of 1–100 Hz for 5 min with eyes closed (sampling rate: 1,000 Hz), for which 64 Ag/AgCl scalp electrodes were evenly mounted on the scalp according to the extended international 10-20 system [NeuroScan SynAmps2 (Compumedics United States, El Paso, TX, United States); references: M1 and M2]. To reduce both external and internal artifacts of EEG data, a series of pre-processing approaches were applied to raw EEG data. Firstly, external artifacts, such as electrocardiography (ECG) and eye-related artifacts (e.g., eyes blinks and horizontal movements), were removed using established mathematical procedures based on regression approach (Semlitsch et al., 1986), and other gross artifacts (e.g., if electrodes showed amplitudes higher than 200 μV) were rejected by visual inspection. To correct the baseline, the DC offset of EEG channels was removed by subtracting the average values of EEG data from each time point for each channel. The baseline corrected EEG data were then band-pass filtered between 1 and 55 Hz using a third order Butterworth IIR filter to remove high-frequency external artifacts, and they were epoched with a length of 4.096 s. Epochs were rejected if they contained significant physiological artifacts (± 100 μV) at any electrode (Kim et al., 2019), and 10 artifact-free epochs extracted for each subject were used for further analysis, as in previous studies (Gudmundsson et al., 2007; Shim et al., 2017). Since a previous study reported that a total of 40 s epoch is sufficient to obtain reliable results of quantifying resting-state EEG data, we used 10 epochs (40.96 s) to extract classification features from resting-state EEG (Gudmundsson et al., 2007).
Source-Level Power Spectrum Densities Feature Extraction
Because related studies reported the superiority of source-level features as compared to sensor-level features (Shim et al., 2016, 2019; Kim J. et al., 2020) as well as of PSD features as compared to network features in terms of classification performance (Kim Y. W. et al., 2020), we used source-level PSD features for the classification of PTSD patients from HC. To estimate the source-level time series, a lead-field matrix was computed using a three-layer (inner skull, outer skull, and scalp) boundary element method (BEM) model, which was constructed from the standard head model (Colin 27) using the OpenMEEG toolbox (Gramfort et al., 2010). An inverse operator was created using a weighted minimum-norm estimation (wMNE) algorithm implemented in Brainstorm toolbox (Tadel et al., 2011). A time-series of source activities at 15,000 cortical vertices was estimated for every time point using the EEG data created by concatenating 10 artifact-free EEG epochs to improve computational efficacy (Kang et al., 2018). After computing current source densities, representative source signals at 68 regions of interests (ROIs) based on the Desikan–Killiany atlas were estimated using the 1st component of principal component analysis (PCA) (Dimitriadis et al., 2018). We excluded 18 ROIs that showed statistical difference between 2 groups (PTSD vs. HC) in terms of the variance explained by the 1st PC because in this case the 1st PCs of 2 groups did not similarly explain the variance of original data (independent t-test; Bonferroni corrected p < 0.05). Therefore, the source signals at the remaining 50 ROIs were used for further data analysis, and the source signals of each ROI was epoched into 4.069 s. Supplementary Table 2 provides the name of all 68 ROIs, their corresponding variances explained by the 1st PCs, and statistical test results, respectively. Time-varying source-level PSDs were then estimated by a complex Morlet-Wavelet method using “ft_freqanalysis” Matlab function implemented in Fieldtrip toolbox (Oostenveld et al., 2011). A crucial input parameter “cgf.width” in “ft_freqanalysis” was set as 3 to appropriately determine a wavelet width (= cfg.width/frequency/pi) according to the recommendation of the Fieldtrip guideline, thereby guaranteeing the accurate estimation of PSDs for all frequencies with the EEG epoch of 4.096 s. The Morlet-wavelet transform with a sinusoidal wave modified by a Gaussian shape was applied to the source-level time series of each ROI. Source-level PSDs of each ROI were independently quantified by averaging time-varying PSDs in 6 frequency bands, i.e., delta [1–4 Hz], theta [4–8 Hz], alpha [8–12 Hz], low-beta [12–22 Hz], high-beta [22–30 Hz], and gamma [30–55Hz]. Note that the beta-frequency band was divided into 2 sub-bands: low-beta [12–22 Hz] and high-beta [22–30 Hz] bands (Kim et al., 2018; Shim et al., 2018). 30 PSD feature sets were finally constructed by integrating the PSD features of different frequency bands (Table 2).
Table 2. Thirty power spectrum density (PSD) feature sets constructed by combining different frequency bands and the number of features for each feature set. Fifty features were extracted for each frequency band.
Classification
30 PSD feature sets were independently tested to evaluate the performance of classifying PTSD patients and HC, thereby finding an optimal combination of PSD frequency bands with respect to classification performance. To this end, classification performances were evaluated using the features by sequentially eliminating recursive features from all features for each feature set using sequential backward selection (SBS) method (García-Laencina et al., 2014). The classification accuracy was evaluated using a linear support vector machine (SVM) classifier (Orru et al., 2012; Alimardani et al., 2018) with a 10-fold cross validation method to prevent overestimate classification performance and improve computation efficacy for each of the 30 feature sets (Kim Y. W. et al., 2020). Note that we also tested 2 other machine-learning classifiers, Random Forest and AdaBoost, as well as 3 deep learning classifiers based on convolutional neural network (CNN), shallow ConvNet (Schirrmeister et al., 2017), EEGNet (Lawhern et al., 2018), and a 13 layers-based deep CNN (Acharya et al., 2018), but we only report SVM results due to its superior classification performance as compared to the other classifiers. In this study, the number of subjects (sample size) in each group was imbalanced (PTSD – 77 and HC – 58), which could lead to biased classification performances. To complement the effect of the imbalance sample size on classification, we employed 2 strategies: (1) a cost-sensitive SVM classifier that modifies the weight of margin penalty proportional to sample size and (2) a balanced classification accuracy [sensitivity + specificity)/2]. Moreover, we evaluated the receiver operating characteristics (ROC) curve and the area under the ROC (AUC) as another performance measure (Long et al., 2017). In addition, to investigate what features (ROIs) were most importantly used for classification, SVM coefficients were evaluated for each ROI (Bosilj et al., 2018) when maximum classification accuracies were obtained for 2 best feature sets extracted from delta and theta frequency bands. Since different features were selected within each cross-validation loop, the absolute values of SVM coefficients were averaged for each of the ROIs selected across cross validation, and they were normalized between 0 and 1. Then, we visualized the ROIs with different sizes and colors along with PSD features where the size of ROIs was proportional to SVM coefficients. A higher SVM coefficient (lager circle) means higher importance in terms of classification based on which neurophysiological interpretation is possible (Bosilj et al., 2018). Figure 1 represents the flowchart of this study.
Figure 1. Flowchart of machine-learning-based classification approach.
Results
Classification Accuracy
Figures 2A,B show the maximum classification accuracies of each frequency band and the corresponding ROC curves with AUC values, respectively: delta: 86.61% with 0.93 (AUC); theta: 82.06% with 0.86; alpha: 72.79% with 0.80; low-beta: 69.36% with 0.79; high-beta: 54.50% with 0.56; and gamma: 50.00% with 0.56 (see Supplementary Table 3 for the balanced classification accuracies of all feature sets). The delta and theta feature sets only showed acceptable classification performance for a practical binary classification system (> 70%) (Perelmouter and Birbaumer, 2000; Hwang et al., 2014). The combinatory feature set of delta and theta features showed almost same classification performance with the delta feature set (86.39%). Other combinatory feature sets, including the delta feature set, also showed comparable classification performance with the delta feature set, e.g., 86.17% for delta + alpha, delta + low-beta, and delta + high-beta, and 85.96% for delta + theta + alpha + low-beta. Note that no feature sets outperformed the delta feature set in terms of classification performance (Supplementary Table 3).
Figure 2. (A) Maximum classification accuracies of each frequency band, and (B) Receiver operating characteristic (ROC) curves with corresponding areas under the curves (AUC) for each frequency band. *indicates a maximum classification accuracy.
Spatial Power Spectrum Density Distribution and Important Features
Figure 3 presents the spatial PSD distributions of delta and theta bands, respectively, and ROIs selected when achieving the maximum classification accuracies. Overall, the PSDs of PTSD patients were considerably reduced, as compared to those of HC (first and second columns of Figure 3). Red circles represent the important ROIs that have SVM coefficients over the upper bound of 95% confidence interval (mean + 2 standard deviation) and blue circles indicates the other selected ROIs. Most of the selected ROIs were overlapped between the 2 frequency bands and they were located in the fronto-temporal area.
Figure 3. Spatial PSD distributions for each group (first and second column) with respect to the frequency band, and the ROIs selected when achieving the maximum classification accuracy for each frequency band (third column). Red circles represent the important ROIs that have SVM coefficients over the upper bound of 95% confidence interval and blue circles indicate the other selected ROIs. The size of circles is proportional to SVM coefficients.
Discussion
In the present study, we investigated the optimal PSD frequency bands to improve the performance of a resting-state EEG-based CAD system to assist the diagnosis of PTSD patients using machine-learning technique. The classification accuracies of the lower frequency PSD feature sets (delta and theta) were significantly higher than those of the relatively higher frequency PSD feature sets (alpha, low-beta, high-beta, and gamma). The best classification performance was obtained when using delta PSDs (86.61% and AUC – 0.93). The features that were selected to attain the best classification accuracy were closely related to the neurophysiological characteristics of PTSD patients, which will be discussed in detail.
Source-Level Power Spectrum Density Features for the Classification of Post-traumatic Stress Disorder Patients
It has been well documented that the results of sensor-level analysis would be distorted and smeared by volume conduction effects due to different tissue conductivities (van den Broek et al., 1998), and thereby sensor-level analysis could not be used to accurately extract neuronal information (Babiloni et al., 1997). The introduction of source-level analysis could redeem the weaknesses of sensor-level analysis (Michel et al., 2004). In fact, many researchers have utilized source-level features to develop CAD systems to assist the diagnosis of psychiatric disorders. For example, previous EEG and MEG studies showed higher classification accuracies when using source-level features than sensor-level features to differentiate schizophrenia patients from HC (Shim et al., 2016; Kim J. et al., 2020), PTSD patients from HC (Shim et al., 2019; Kim Y. W. et al., 2020), and PTSD patients from major depressive disorder patients (Shim et al., 2019). Therefore, in this study, we only considered source-level features to classify PTSD patients from HC rather than using sensor-level analysis, and achieved reasonable classification performance.
Power spectrum densities (PSDs) have been mainly used as a classification feature type to develop a CAD system to assist in the diagnosis of psychiatric patients (Cassani et al., 2018; Radenković and Lopez, 2019) because distinct and abnormal PSD patterns were shown in psychiatric disorders, as compared to HC (Newson and Thiagarajan, 2019); these abnormal PSD patterns were also related to genetic traits of psychiatric disorders (Venables et al., 2009). PTSD patients also showed altered PSD patterns compared to HC; PTSD patients revealed significantly diminished PSDs in delta and theta bands (McFarlane et al., 2005; Veltmeyer et al., 2006; Newson and Thiagarajan, 2019). Moreover, the diminished PSDs in slow-frequency bands significantly correlated with patients’ symptom scores (Clinician-Administered PTSD scale, CAPS) (Veltmeyer et al., 2006), indicating that the altered PSD patterns would reflect the pathophysiology of PTSD patients, such as re-experiencing and arousal (Veltmeyer et al., 2006). Abnormal PSD patterns of PTSD patients were also clearly observed in this study (Figure 3), and they played an important role in classifying PTSD patients from HC. However, unfortunately, we did not find any significant relationship between the PSD features used to obtain a relatively high performance (i.e., 18 features for delta band and 35 features for theta band) and symptom scores (i.e., IES-R, BAI, and BDI) in this study. This would be because these features were selected from a machine learning perspective for better classification; low frequency PSD features for each group (PTSD and HC) were separately grouped with a discriminable distance (i.e., relatively high PSD values for HC vs. low PSD values for PTSD, as shown in Figure 3), but which does not guarantee that the low frequency PSD features of PTSD patients are necessarily correlated with neuropsychological estimates (e.g., symptom scores). Note that the discriminable low frequency PSD features were selected by comparing those of PTSD patients and HC, but symptom scores were acquired only from PTSD patients. Because of the mentioned reason, the low frequency PSD features of PTSD patients did not seem to be proportional to symptom scores even though the low frequency PSD features can be used for the discrimination of PTSD patients from HC. In order to find a certain feature set that is directly correlated with neuropsychological estimates, regression should be used instead of classification by focusing on only PTSD data, but which would be beyond the scope of this study, focusing on the discrimination of 2 groups. Nevertheless, it is important to investigate neuropsychological traits using objective neurophysiological biomarkers in order to improve the understanding of neural mechanism in psychiatric disorders. Therefore, we will keep going to develop a new feature-based CAD system that can provide both the accurate diagnosis of psychiatric patients in terms of machine learning and the understanding of neurophysiological traits in terms of neuroscience by adopting other EEG-based metrics, such as effective connectivity and complexity measures.
Crucial Role of Slow Brain Waves
In the present study, the best classification performance was obtained when using delta PSD features. The best classification performance obtained in this study (86.61%) was considerably higher than that of our previous PTSD study by approximately 13% (Kim Y. W. et al., 2020). Unlike the previous study, we used delta PSD features when differentiating PTSD patients and HC, which mainly led to the relatively higher classification accuracy. Moreover, when using theta PSD features, a reasonable classification accuracy of 82.06% was also achieved, which was also higher than that of the previous study (Kim Y. W. et al., 2020). The performance improvement would be caused by the distinct characteristics of slow EEG waves of PTSD patients in the delta and theta frequency bands. Many previous studies reported that slow EEG waves were a typical endophenotype of PTSD patients (Franke et al., 2016; Sheerin et al., 2018; Newson and Thiagarajan, 2019). As mentioned above, PTSD patients showed significantly decreased delta and theta PSDs compared to HC, and the decreased slow-frequency PSDs were closely associated with their altered symptoms, such as arousal and numbing (Franke et al., 2016; Sheerin et al., 2018). That is, due to the distinct pathological traits of slow EEG waves in PTSD patients, the present study could significantly improve the classification performance for the diagnosis of PTSD patients, compared to that of our previous study (Kim Y. W. et al., 2020).
Spatial Distribution of Power Spectrum Density Features
All PSD features selected when attaining maximum classification accuracies for delta and theta bands were extracted from fronto-temporal areas, such as frontal pole, opercular, anterior cingulate gyrus, superior temporal gyrus, and temporal pole (third column of Figure 3). The mean PSD values of each ROI of the mentioned fronto-temporal areas were significantly smaller in PTSD patients than those in HC. According to previous studies, the frontal regions, including the frontal pole and opercular part of the inferior frontal area, are involved in emotion regulation processing (Grecucci et al., 2013; Bramson et al., 2019) and temporal areas are significantly related to rumination symptom of PTSD patients (Ferdek et al., 2016); PTSD patients showed significantly reduced brain activation in both brain areas in resting-state (Rabe et al., 2006), as shown in this study. Therefore, it is neurophysiologically plausible to obtain better classification performance for PTSD diagnosis when using PSD features extracted from fronto-temporal areas. This also indicates that employing only fronto-temporal areas might be sufficient to implement a reliable CAD system for PTSD diagnosis, thereby facilitating the development of a more clinically practical CAD system.
In this study, we used EEG data measured in resting state during which default mode network (DMN) or salience network is more active than the task period (Choi et al., 2021), and thus DMN might be usefully used for PTSD diagnosis. We investigated the ROIs selected when the maximum classification accuracy was obtained using delta PSD features, and found that 16 of 18 DMN ROIs were selected together with 2 other ROIs, i.e., entorhinal and para hippocampal areas. Both entorhinal and para hippocampal areas are known to be related to clinical symptoms in PTSD patients, such as altered cognitive function and disrupted memory (Kerr et al., 2007). This result means that it is necessary to use the features extracted from the ROIs related to neuropsychological traits of patients together with those related to current brain state (resting state in this study) to obtain a relatively high diagnosis accuracy.
Machine Learning Approach With Interpretable Electroencephalogram Biomarkers
In recent years, machine learning approaches have received increasing attention in the development of an EEG-based CAD system to assist the accurate diagnosis of psychiatric disorders (Bzdok and Meyer-Lindenberg, 2018; Fazel and O’Reilly, 2020; Rutherford, 2020; Koutsouleris et al., 2021). In particular, it is important to use neurophysiologically interpretable EEG biomarkers in the development of an EEG-based CAD system to improve its reliability, as in deep learning (Sturm et al., 2016). In this study, we introduced low frequency EEG features that were closely associated with the neurophysiological characteristics of PTSD patients (Franke et al., 2016; Vergara et al., 2017; Newson and Thiagarajan, 2019), and thereby not only the proposed low frequency EEG biomarkers could significantly improve the performance of the CAD system, but also they were neurophysiologically interpretable. To evaluate the potential of low frequency EEG features for the performance improvement of the CAD system, we used a traditional machine learning approach using PSDs as features and SVM as a classifier, and the improved classification accuracy of 86.61% would not be still enough to be used in clinics. Even though we also used 3 different CNN-based deep learning algorithms for the diagnosis of PTSD patients as mentioned in the method section, we could not obtain a comparable classification accuracy to that obtained using the SVM approach (73.13% for Shallow ConvNet, 76.15% for EEGNet, and 77.03% for 13-Layers CNN). This result would be derived due to a relatively small number of samples (data amount) used in deep learning algorithms; traditional machine-learning models show better performance than that of deep learning models with relatively small samples (Kumar and Manash, 2019). Therefore, we will keep trying to improve the diagnosis accuracy by introducing more advanced algorithms, such as data augmentation to compensate a small number of samples in our future studies (Wang et al., 2018), thereby developing a clinically usable CAD system for PTSD patients.
Limitations
First, since all recruited PTSD patients were on medication, we could not control the compounding effects from medications. It has been reported that EEG characteristics can be changed by medications. For example, antipsychotics increased slow waves, such as delta and theta (Amann et al., 2003), antidepressants modified alpha patterns (Bruder et al., 2008), and sedative-hypnotics enhanced low frequency power (< 15 Hz) (Ferri et al., 2017). All psychotic medications tend to increase PSD values in particular for relatively low frequency bands, but delta and theta PSDs for PTSD patients were still significantly lower than those of HC, as shown in Figure 3. Therefore, it could be reasonably thought that delta and theta PSDs of PTSD patients could be increased by medications, but which did not reach those of HC. Thereby, significant difference between two groups was still kept in terms of low frequency PSD values, and they were used as key features for accurate diagnosis. Our previous study also used same types of medications used in this study, i.e., antipsychotics, antidepressants, and sedative-hypnotics (Kim Y. W. et al., 2020). Because patients recruited in the previous study were independent from those recruited in the present study, it is impossible to directly compare the results of the two studies. However, assuming that similar medication impacts occurred on EEG characteristics due to using same types of medications for both groups of PTSD patients, it can be reasonably thought that the enhanced classification performance in this study would be caused by using optimal PSD features as compared to the previous study, not medication effects. Further studies should follow with drug-naïve PTSD patients to accurately investigate the medication impact on the development of an EEG-based CAD system, which can also allow for PTSD diagnosis at initial screening stages. Second, even though we controlled other psychiatric illness, we did not control comorbid depression. Third, we used a total of 40.96-s EEG data measured using 64 electrodes for data analysis (e.g., source estimation). It was reported that reliable source signal estimation can be possible using the EEG data acquired from more than 60 channels (Michel et al., 2004), and also reliable resting-state EEG data analysis with a data length of more than 40 s (Gudmundsson et al., 2007). However, since it is obvious that use of longer EEG data measured from a high-density EEG system with more channels (e.g., > 128) can allow for more reliable EEG data analysis, including source imaging (Haartsen et al., 2020), using such the high-density EEG data is necessary to improve the reliability and accuracy of EEG data analysis in future studies. Fourth, PTSD patients showed different EEG characteristics according to gender. For example, female patients showed enhanced alpha asymmetry than female HC (Metzger et al., 2004), while male patients showed decreased source activities in theta band and decreased alpha PSD as compared to male HC (Jokić-begić and Begić, 2003; Todder et al., 2012). Therefore, diagnostic performance could be further improved by considering gender-specific EEG features, which would be one of the interesting future topics.
Conclusion
We investigated optimal PSD frequency bands to improve the classification performances of a resting-state EEG-based CAD system for precise PTSD diagnosis. Low-frequency PSD features in delta and theta frequency bands showed significantly higher classification performances than relatively high-frequency PSD features (alpha, low-beta, high-beta, and gamma), and the best classification performance (86.61% and AUC of 0.93) was obtained when using delta PSD features. In addition, most meaningful features were extracted from fronto-tempora areas, which coincided with the neurophysiological findings of previous EEG-based PTSD studies. Although the present study showed relatively high classification performances between PTSD patients and HC, there is still room to improve the classification performance by introducing novel feature extraction and classification methods based on deep learning algorithms, which could be interesting future research topics.
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 Inje University Ilsan Paik Hospital (2015-09-018). The patients/participants provided their written informed consent to participate in this study.
Author Contributions
MS, S-HL, and H-JH: conception and design of the work and interpretation of data. MS and S-HL: data acquisition and analysis. MS, C-HI, S-HL, and H-JH: writing (original draft preparation and editing). All authors read and contributed to the final version of the manuscript.
Funding
This work was supported by the Institute for Information & Communications Technology Planning & Evaluation (IITP) grant funded by the Korea government (No. 2017-0-00451; Development of BCI based Brain and Cognitive Computing Technology for Recognizing User’s Intentions using Deep Learning), and by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the MSIT (2020R1A4A1017775) and by the Ministry of Education (NRF-2019R1I1A1A01063313), and the Korea Science and Engineering Foundation (KOSEF), funded by the Korean government (NRF-2018R1A2A2A05018505).
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fninf.2022.811756/full#supplementary-material
References
Acharya, U. R., Oh, S. L., Hagiwara, Y., Tan, J. H., Adeli, H., and Subha, D. P. (2018). Automated EEG-based screening of depression using deep convolutional neural network. Comput. Methods Programs Biomed. 161, 103–113. doi: 10.1016/j.cmpb.2018.04.012
Alimardani, F., Cho, J.-H., Boostani, R., and Hwang, H.-J. (2018). Classification of bipolar disorder and schizophrenia using steady-state visual evoked potential based features. IEEE Access 6, 40379–40388. doi: 10.1109/access.2018.2854555
Amann, B. L., Pogarell, O., Mergl, R., Juckel, G., Grunze, H., Mulert, C., et al. (2003). EEG abnormalities associated with antipsychotics: a comparison of quetiapine, olanzapine, haloperidol and healthy subjects. Hum. Psychopharmacol. 18, 641–646. doi: 10.1002/hup.537
Babiloni, F., Babiloni, C., Carducci, F., Fattorini, L., Anello, C., Onorati, P., et al. (1997). High resolution EEG: a new model-dependent spatial deblurring method using a realistically-shaped MR-constructed subject’s head model. Electroencephalogr. Clin. Neurophysiol. 102, 69–80. doi: 10.1016/s0921-884x(96)96508-x
Beck, A. T., and Steer, R. A. (1990). Manual for the Beck Anxiety Inventory. San Antonio, TX: Psychological Corporation.
Beck, A. T., Steer, R. A., and Brown, G. K. (1996). Manual for the Beck Depression Inventory-ii, Vol. 1. San Antonio, TX: Psychological Corporation, 82.
Bosilj, P., Duckett, T., and Cielniak, G. (2018). Connected attribute morphology for unified vegetation segmentation and classification in precision agriculture. Comput. Ind. 98, 226–240. doi: 10.1016/j.compind.2018.02.003
Bramson, B., Folloni, D., Verhagen, L., Hartogsveld, B., Mars, R. B., Toni, I., et al. (2019). Human lateral frontal pole contributes to control over social-emotional behaviour. bioRxiv [Preprint] doi: 10.1101/584896
Bruder, G. E., Sedoruk, J. P., Stewart, J. W., McGrath, P. J., Quitkin, F. M., and Tenke, C. E. (2008). EEG alpha measures predict therapeutic response to an SSRI antidepressant: pre and post treatment findings. Biol. Psychiatry 63:1171. doi: 10.1016/j.biopsych.2007.10.009
Bryant, R. (2011). Post-traumatic stress disorder vs traumatic brain injury. Dialogues Clin. Neurosci. 13:251. doi: 10.31887/dcns.2011.13.2/rbryant
Burle, B., Spieser, L., Roger, C., Casini, L., Hasbroucq, T., and Vidal, F. (2015). Spatial and temporal resolutions of EEG: is it really black and white? A scalp current density view. Int. J. Psychophysiol. 97, 210–220. doi: 10.1016/j.ijpsycho.2015.05.004
Bzdok, D., and Meyer-Lindenberg, A. (2018). Machine learning for precision psychiatry: opportunities and challenges. Biol. Psychiatry Cogn. Neurosci. Neuroimaging 3, 223–230. doi: 10.1016/j.bpsc.2017.11.007
Cassani, R., Estarellas, M., San-Martin, R., Fraga, F. J., and Falk, T. H. (2018). Systematic review on resting-state EEG for Alzheimer’s disease diagnosis and progression assessment. Dis. Markers 2018:5174815. doi: 10.1155/2018/5174815
Chen, A. C. (2001). New perspectives in EEG/MEG brain mapping and PET/fMRI neuroimaging of human pain. Int. J. Psychophysiol. 42, 147–159. doi: 10.1016/s0167-8760(01)00163-5
Choi, K. M., Kim, J. Y., Kim, Y. W., Han, J. W., Im, C. H., and Lee, S. H. (2021). Comparative analysis of default mode networks in major psychiatric disorders using resting-state EEG. Sci. Rep. 11, 1–10. doi: 10.1038/s41598-021-03736-4
Dimitriadis, S. I., López, M. E., Bruña, R., Cuesta, P., Marcos, A., Maestú, F., et al. (2018). How to build a functional connectomic biomarker for mild cognitive impairment from source reconstructed meg resting-state activity: the combination of ROI representation and connectivity estimator matters. Front. Neurosci. 12:306. doi: 10.3389/fnins.2018.00306
Elhai, J. D., de Francisco Carvalho, L., Miguel, F. K., Palmieri, P. A., Primi, R., and Frueh, B. C. (2011). Testing whether posttraumatic stress disorder and major depressive disorder are similar or unique constructs. J. Anxiety Disord. 25, 404–410. doi: 10.1016/j.janxdis.2010.11.003
Fazel, S., and O’Reilly, L. (2020). Machine learning for suicide research–can it improve risk factor identification? JAMA Psychiatry 77, 13–14. doi: 10.1001/jamapsychiatry.2019.2896
Ferdek, M. A., van Rijn, C. M., and Wyczesany, M. (2016). Depressive rumination and the emotional control circuit: an EEG localization and effective connectivity study. Cogn. Affect. Behav. Neurosci. 16, 1099–1113. doi: 10.3758/s13415-016-0456-x
Ferri, R., Rundo, F., Silvani, A., Zucconi, M., Bruni, O., Ferini-Strambi, L., et al. (2017). REM sleep EEG instability in REM sleep behavior disorder and clonazepam effects. Sleep 40, zsx080. doi: 10.1093/sleep/zsx080
Franke, L. M., Walker, W. C., Hoke, K. W., and Wares, J. R. (2016). Distinction in EEg slow oscillations between chronic mild traumatic brain injury and ptsd. Int. J. Psychophysiol. 106, 21–29. doi: 10.1016/j.ijpsycho.2016.05.010
Friedman, M. J., Resick, P. A., Bryant, R. A., and Brewin, C. R. (2011). Considering PTSD for DSM-5. Depress. Anxiety 28, 750–769. doi: 10.1002/da.20767
García-Laencina, P. J., Rodríguez-Bermudez, G., and Roca-Dorda, J. (2014). Exploring dimensionality reduction of EEG features in motor imagery task classification. Expert Syst. Appl. 41, 5285–5295. doi: 10.1016/j.eswa.2014.02.043
Ginzburg, K., Ein-Dor, T., and Solomon, Z. (2010). Comorbidity of posttraumatic stress disorder, anxiety and depression: a 20-year longitudinal study of war veterans. J. Affect. Disord. 123, 249–257. doi: 10.1016/j.jad.2009.08.006
Gramfort, A., Papadopoulo, T., Olivi, E., and Clerc, M. (2010). OpenMEEG: opensource software for quasistatic bioelectromagnetics. Biomed. Eng. Online 9:45. doi: 10.1186/1475-925X-9-45
Grecucci, A., Giorgetta, C., Bonini, N., and Sanfey, A. G. (2013). Reappraising social emotions: the role of inferior frontal gyrus, temporo-parietal junction and insula in interpersonal emotion regulation. Front. Hum. Neurosci. 7:523. doi: 10.3389/fnhum.2013.00523
Gros, D. F., Price, M., Magruder, K. M., and Frueh, B. C. (2012). Symptom overlap in posttraumatic stress disorder and major depression. Psychiatry Res. 196, 267–270. doi: 10.1016/j.psychres.2011.10.022
Grubaugh, A. L., Long, M. E., Elhai, J. D., Frueh, B. C., and Magruder, K. M. (2010). An examination of the construct validity of posttraumatic stress disorder with veterans using a revised criterion set. Behav. Res. Ther. 48, 909–914. doi: 10.1016/j.brat.2010.05.019
Gudmundsson, S., Runarsson, T. P., Sigurdsson, S., Eiriksdottir, G., and Johnsen, K. (2007). Reliability of quantitative EEG features. Clin. Neurophysiol. 118, 2162–2171. doi: 10.1016/j.clinph.2007.06.018
Haartsen, R., van der Velde, B., Jones, E. J., Johnson, M. H., and Kemner, C. (2020). Using multiple short epochs optimises the stability of infant EEG connectivity parameters. Sci. Rep. 10, 1–13. doi: 10.1038/s41598-020-68981-5
Hwang, H.-J., Lim, J.-H., Kim, D.-W., and Im, C.-H. (2014). Evaluation of various mental task combinations for near-infrared spectroscopy-based brain-computer interfaces. J. Biomed. Opt. 19:077005. doi: 10.1117/1.JBO.19.7.077005
Jokić-begić, N., and Begić, D. (2003). Quantitative electroencephalogram (qEEG) in combat veterans with post-traumatic stress disorder (PTSD). Nord. J. Psychiatry 57, 351–355. doi: 10.1080/08039480310002688
Kang, S. S., MacDonald, A. W. III, Chafee, M. V., Im, C. H., Bernat, E. M., Davenport, N. D., et al. (2018). Abnormal cortical neural synchrony during working memory in schizophrenia. Clin. Neurophysiol. 129, 210–221. doi: 10.1016/j.clinph.2017.10.024
Kerr, K. M., Agster, K. L., Furtak, S. C., and Burwell, R. D. (2007). Functional neuroanatomy of the parahippocampal region: the lateral and medial entorhinal areas. Hippocampus 17, 697–708. doi: 10.1002/hipo.20315
Kim, J., Kim, M.-Y., Kwon, H., Kim, J.-W., Im, W.-Y., Lee, S. M., et al. (2020). Feature optimization method for machine learning-based diagnosis of schizophrenia using magnetoencephalography. J. Neurosci. Methods 338:108688. doi: 10.1016/j.jneumeth.2020.108688
Kim, J.-Y., Son, J.-B., Leem, H.-S., and Lee, S.-H. (2019). Psychophysiological alteration after virtual reality experiences using smartphone-assisted head mount displays: an EEG-based source localization study. Appl. Sci. 9:2501. doi: 10.3390/app9122501
Kim, S., Kim, J. S., Shim, M., Im, C.-H., and Lee, S.-H. (2018). Altered cortical functional network during behavioral inhibition in individuals with childhood trauma. Sci. Rep. 8, 1–10. doi: 10.1038/s41598-018-28329-6
Kim, Y.-W., Kim, S., Shim, M., Jin, M. J., Jeon, H., Lee, S.-H., et al. (2020). Riemannian classifier enhances the accuracy of machine-learning-based diagnosis of PTSD using resting EEG. Prog. Neuropsychopharmacol. Biol. Psychiatry 102:109960. doi: 10.1016/j.pnpbp.2020.109960
Koutsouleris, N., Dwyer, D. B., Degenhardt, F., Maj, C., Urquijo-Castro, M. F., Sanfelici, R., et al. (2021). Multimodal machine learning workflows for prediction of psychosis in patients with clinical high-risk syndromes and recent-onset depression. JAMA Psychiatry 78, 195–209. doi: 10.1001/jamapsychiatry.2020.3604
Kumar, P. R., and Manash, E. B. K. (2019). “Deep learning: a branch of machine learning,” in Paper Presented of the Journal of Physics: Conference Series, (Andhra Pradesh, India), Vol. 1228:012045. doi: 10.1016/j.jad.2020.11.104
Lawhern, V. J., Solon, A. J., Waytowich, N. R., Gordon, S. M., Hung, C. P., and Lance, B. J. (2018). EEGNet: a compact convolutional neural network for EEG-based brain—computer interfaces. J. Neural Eng. 15:056013. doi: 10.1088/1741-2552/aace8c
Long, X., Chen, L., Jiang, C., Zhang, L., and Alzheimer’s Disease Neuroimaging Initiative. (2017). Prediction and classification of Alzheimer disease based on quantification of MRI deformation. PLoS One 12:e0173372. doi: 10.1371/journal.pone.0173372
Magruder, K. M., Frueh, B. C., Knapp, R. G., Davis, L., Hamner, M. B., Martin, R. H., et al. (2005). Prevalence of posttraumatic stress disorder in veterans affairs primary care clinics. Gen. Hosp. Psychiatry 27, 169–179. doi: 10.1016/j.genhosppsych.2004.11.001
McBride, J., Zhao, X., Nichols, T., Abdul-Ahad, T., Wilson, M., Vagnini, V., et al. (2011). “Classification of traumatic brain injury using support vector machine analysis of event-related tsallis entropy,” in Proceedings of the 2011 Biomedical Sciences and Engineering Conference: Image Informatics and Analytics in Biomedicine, (Piscataway, NJ: IEEE), 1–4. doi: 10.29284/ijasis.2.1.2016.1-4
McFarlane, A., Clark, C. R., Bryant, R. A., Williams, L. M., Niaura, R., Paul, R. H., et al. (2005). The impact of early life stress on psychophysiological, personality and behavioral measures in 740 non-clinical subjects. J. Integr. Neurosci. 4, 27–40. doi: 10.1142/s0219635205000689
Metzger, L. J., Paige, S. R., Carson, M. A., Lasko, N. B., Paulus, L. A., Pitman, R. K., et al. (2004). PTSD arousal and depression symptoms associated with increased right-sided parietal EEG asymmetry. J. Abnorm. Psychol. 113:324. doi: 10.1037/0021-843X.113.2.324
Michel, C. M., Murray, M. M., Lantz, G., Gonzalez, S., Spinelli, L., and de Peralta, R. G. (2004). EEG source imaging. Clin. Neurophysiol. 115, 2195–2222. doi: 10.1016/j.clinph.2004.06.001
Mitra, J., Shen, K-k, Ghose, S., Bourgeat, P., Fripp, J., Salvado, O., et al. (2016). Statistical machine learning to identify traumatic brain injury (TBI) from structural disconnections of white matter networks. Neuroimage 129, 247–259. doi: 10.1016/j.neuroimage.2016.01.056
Newson, J. J., and Thiagarajan, T. C. (2019). EEG frequency bands in psychiatric disorders: a review of resting state studies. Front. Hum. Neurosci. 12:521. doi: 10.3389/fnhum.2018.00521
Oostenveld, R., Fries, P., Maris, E., and Schoffelen, J.-M. (2011). Fieldtrip: open source software for advanced analysis of MEG,EEG, and invasive electrophysiological data. Comput. Intell. Neurosci. 2011:156869. doi: 10.1155/2011/156869
Orru, G., Pettersson-Yeo, W., Marquand, A. F., Sartori, G., and Mechelli, A. (2012). Using support vector machine to identify imaging biomarkers of neurological and psychiatric disease: a critical review. Neurosci. Biobehav. Rev. 36, 1140–1152. doi: 10.1016/j.neubiorev.2012.01.004
Perelmouter, J., and Birbaumer, N. (2000). A binary spelling interface with random errors. IEEE Trans. Rehabil. Eng. 8, 227–232. doi: 10.1109/86.847824
Rabe, S., Beauducel, A., Zöllner, T., Maercker, A., and Karl, A. (2006). Regional brain electrical activity in posttraumatic stress disorder after motor vehicle accident. J. Abnorm. Psychol. 115:687. doi: 10.1037/0021-843X.115.4.687
Radenković, M. Č, and Lopez, V. L. (2019). Machine learning approaches for detecting the depression from resting-state electroencephalogram (EEG): a review study. arXiv [Preprint]. arXiv:190903115. doi: 10.2196/19548
Rangaprakash, D., Deshpande, G., Daniel, T. A., Goodman, A. M., Robinson, J. L., Salibi, N., et al. (2017). Compromised hippocampus-striatum pathway as a potential imaging biomarker of mild-traumatic brain injury and posttraumatic stress disorder. Hum. Brain Mapp. 38, 2843–2864. doi: 10.1002/hbm.23551
Rutherford, S. (2020). The promise of machine learning for psychiatry. Biol. Psychiatry 88, e53–e55.
Saletu, B., Anderer, P., Paulus, E., Grünberger, J., Wicke, L., Neuhold, A., et al. (1991). EEG brain mapping in diagnostic and therapeutic assessment of dementia. Alzheimer Dis. Assoc. Disord. 5, S57–S75. doi: 10.1097/00002093-199100051-00010
Schirrmeister, R. T., Springenberg, J. T., Fiederer, L. D. J., Glasstetter, M., Eggensperger, K., Tangermann, M., et al. (2017). Deep learning with convolutional neural networks for EEG decoding and visualization. Hum. Brain. Mapp. 38, 5391–5420. doi: 10.1002/hbm.23730
Semlitsch, H. V., Anderer, P., Schuster, P., and Presslich, O. (1986). A solution for reliable and valid reduction of ocular artifacts, applied to the p300 ERP. Psychophysiology 23, 695–703. doi: 10.1111/j.1469-8986.1986.tb00696.x
Sheerin, C. M., Franke, L. M., Aggen, S. H., Amstadter, A. B., and Walker, W. C. (2018). Evaluating the contribution of EEG power profiles to characterize and discriminate posttraumatic stress symptom factors in a combat-exposed population. Clin. EEG Neurosci. 49, 379–387. doi: 10.1177/1550059418767583
Shim, M., Hwang, H.-J., Kim, D.-W., Lee, S.-H., and Im, C.-H. (2016). Machine-learning-based diagnosis of schizophrenia using combined sensor-level and source-level EEG features. Schizophr. Res. 176, 314–319. doi: 10.1016/j.schres.2016.05.007
Shim, M., Im, C.-H., Kim, Y.-W., and Lee, S.-H. (2018). Altered cortical functional network in major depressive disorder: a resting-state electroencephalogram study. NeuroImage Clin. 19, 1000–1007. doi: 10.1016/j.nicl.2018.06.012
Shim, M., Im, C.-H., and Lee, S. (2017). Disrupted cortical brain network in post-traumatic stress disorder patients: a resting-state electroencephalographic study. Transl. Psychiatry 7, e1231–e1231. doi: 10.1038/tp.2017.200
Shim, M., Jin, M. J., Im, C.-H., and Lee, S.-H. (2019). Machine-learning-based classification between post-traumatic stress disorder and major depressive disorder using P300 features. NeuroImage Clin. 24:102001. doi: 10.1016/j.nicl.2019.102001
Sturm, I., Lapuschkin, S., Samek, W., and Müller, K.-R. (2016). Interpretable deep neural networks for single-trial EEG classification. J. Neurosci. Methods 274, 141–145. doi: 10.1016/j.jneumeth.2016.10.008
Sumpter, R. E., and McMillan, T. M. (2005). Misdiagnosis of post-traumatic stress disorder following severe traumatic brain injury. Br. J. Psychiatry 186, 423–426. doi: 10.1192/bjp.186.5.423
Tadel, F., Baillet, S., Mosher, J. C., Pantazis, D., and Leahy, R. M. (2011). Brainstorm: a user-friendly application for MEG/EEG analysis. Comput. Intell. Neurosci. 2011:879716. doi: 10.1155/2011/879716
Todder, D., Levine, J., Abujumah, A., Mater, M., Cohen, H., and Kaplan, Z. (2012). The quantitative electroencephalogram and the low-resolution electrical tomographic analysis in posttraumatic stress disorder. Clin. EEG Neurosci. 43, 48–53. doi: 10.1177/1550059411428716
van den Broek, S. P., Reinders, F., Donderwinkel, M., and Peters, M. (1998). Volume conduction effects in EEG and MEG. Electroencephalogr. Clin. Neurophysiol 106, 522–534. doi: 10.1016/s0013-4694(97)00147-8
Veltmeyer, M. D., McFarlane, A. C., Bryant, R. A., Mayo, T., Gordon, E., and Clark, C. R. (2006). Integrative assessment of brain function in PTSD: brain stability and working memory. J. Integr. Neurosci. 5, 123–138. doi: 10.1142/s0219635206001057
Venables, N. C., Bernat, E. M., and Sponheim, S. R. (2009). Genetic and disorder-specific aspects of resting state EEG abnormalities in schizophrenia. Schizophr. Bull. 35, 826–839. doi: 10.1093/schbul/sbn021
Vergara, V. M., Mayer, A. R., Damaraju, E., Kiehl, K. A., and Calhoun, V. (2017). Detection of mild traumatic brain injury by machine learning classification using resting state functional network connectivity and fractional anisotropy. J. Neurotrauma 34, 1045–1053. doi: 10.1089/neu.2016.4526
Vergara, V. M., Mayer, A. R., Kiehl, K. A., and Calhoun, V. D. (2018). Dynamic functional network connectivity discriminates mild traumatic brain injury through machine learning. NeuroImage Clin. 19, 30–37. doi: 10.1016/j.nicl.2018.03.017
Wang, F., Zhong, S. H., Peng, J., Jiang, J., and Liu, Y. (2018). “Data augmentation for eeg-based emotion recognition with deep convolutional neural networks,” in Proceedings of the International Conference on Multimedia Modeling, eds K. Schoeffmann, T. H. Chalidabhongse, C. W. Ngo, S. Aramvith, N. E. O’Connor, Y.-S. Ho, et al. (Cham: Springer International Publishing), 82–93. doi: 10.1155/2021/2520394
Weiss, D. S. (2007). “The impact of event scale: revised,” in Cross-Cultural Assessment of Psychological Trauma and PTSD, eds J. P. Wilson and C. S.-k Tang (Heidelberg: Springer), 219–238.
Wienbruch, C., Candia, V., Svensson, J., Kleiser, R., and Kollias, S. S. (2006). A portable and low-cost fMRI compatible pneumatic system for the investigation of the somatosensensory system in clinical and research environments. Neurosci. Lett. 398, 183–188. doi: 10.1016/j.neulet.2006.01.025
Keywords: machine-learning technique, classification, computer-aided diagnosis, resting-state electroencephalogram (EEG), slow-frequency EEG oscillation, post-traumatic stress disorder (PTSD)
Citation: Shim M, Im C-H, Lee S-H and Hwang H-J (2022) Enhanced Performance by Interpretable Low-Frequency Electroencephalogram Oscillations in the Machine Learning-Based Diagnosis of Post-traumatic Stress Disorder. Front. Neuroinform. 16:811756. doi: 10.3389/fninf.2022.811756
Received: 09 November 2021; Accepted: 28 March 2022;
Published: 26 April 2022.
Edited by:
Liana Fattore, CNR Neuroscience Institute (IN), ItalyReviewed by:
Stavros I. Dimitriadis, Greek Association of Alzheimer’s Disease and Related Disorders, GreeceAbdulhakim AL-Ezzi, University of Technology Petronas, Malaysia
Görsev Yener, İzmir University of Economics, Turkey
Copyright © 2022 Shim, Im, Lee and Hwang. 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: Seung-Hwan Lee, [email protected]; Han-Jeong Hwang, [email protected]