Original Research ARTICLE
Estimation of Symptom Severity Scores for Patients with Schizophrenia Using ERP Source Activations during a Facial Affect Discrimination Task
- 1Department of Biomedical Engineering, Chonnam National University, Yeosu, South Korea
- 2Psychiatry Department, Ilsan Paik Hospital, Inje University, Goyang, South Korea
- 3Department of Biomedical Engineering, Hanyang University, Seoul, South Korea
Precise diagnosis of psychiatric diseases and a comprehensive assessment of a patient's symptom severity are important in order to establish a successful treatment strategy for each patient. Although great efforts have been devoted to searching for diagnostic biomarkers of schizophrenia over the past several decades, no study has yet investigated how accurately these biomarkers are able to estimate an individual patient's symptom severity. In this study, we applied electrophysiological biomarkers obtained from electroencephalography (EEG) analyses to an estimation of symptom severity scores of patients with schizophrenia. EEG signals were recorded from 23 patients while they performed a facial affect discrimination task. Based on the source current density analysis results, we extracted voxels that showed a strong correlation between source activity and symptom scores. We then built a prediction model to estimate the symptom severity scores of each patient using the source activations of the selected voxels. The symptom scores of the Positive and Negative Syndrome Scale (PANSS) were estimated using the linear prediction model. The results of leave-one-out cross validation (LOOCV) showed that the mean errors of the estimated symptom scores were 3.34 ± 2.40 and 3.90 ± 3.01 for the Positive and Negative PANSS scores, respectively. The current pilot study is the first attempt to estimate symptom severity scores in schizophrenia using quantitative EEG features. It is expected that the present method can be extended to other cognitive paradigms or other psychological illnesses.
Since patients with schizophrenia have their own unique signs and symptoms that generally show dramatic changes over the progress of the illness or during treatment, it is of great importance to precisely evaluate the symptom severity of each patient. Although the assessment of symptom severity in schizophrenia is critical in establishing successful treatment strategies or evaluating the effectiveness of treatments, only few quantitative diagnosis tools exist to evaluate the symptom severity of individual patients with schizophrenia (Lakhan, 2006).
Thus far, symptom severity in schizophrenia has been generally diagnosed based on interview-based tests performed by trained psychiatrists. Although the inter-rater reliability or test-retest reliability of these criteria is reported to be high (Kay et al., 1987; Bell et al., 1994), it is still possible that the results might be highly influenced by initial assessment (Mortimer, 2007) or subjective opinions of the psychiatrists (Norman et al., 1996). Moreover, diagnosis of the negative symptoms of schizophrenia is relatively more difficult (Lindström et al., 1994; Norman et al., 1996) than that of positive symptoms, which are based on evident symptoms. Patients with negative symptoms might not show any evident signs or may often show general behaviors overlapping with other mental illnesses (Andreasen, 1982). Therefore, the development of quantitative methods to diagnose schizophrenia or evaluate symptom severity in individual patients can potentially help to enhance the diagnostic accuracy of schizophrenia.
Over the past several decades, researchers have put great effort into investigating predictable characteristics of schizophrenia. For instance, several studies have reported the possibility of diagnosing schizophrenia based on either functional (Sabeti et al., 2007, 2009) or structural characteristics (Takayanagi et al., 2011) of the brain. Various studies have also shown that electrophysiological measures such as P300 (Shenton et al., 1989; Strik et al., 1993; Kim et al., 2014), auditory steady-state response (ASSR) (Spencer et al., 2008; Shin et al., 2011), loudness-dependent auditory evoked potential (LDAEP) (Ostermann et al., 2012; Wyss et al., 2013), and mismatch negativity (MMN) (Hirayasu et al., 1998; Javitt et al., 2000; Youn et al., 2003) are strongly correlated with symptom severity in schizophrenia and thus might be used as trait or state markers of the disease. There is a growing consensus that the precise estimation of symptom severity is of great necessity as it can be used to predict quality of life (Bow-Thomas et al., 1999), long-term outcome (Harrison et al., 1996; Ho et al., 1998), relapse (Birchwood et al., 1989), and functional outcomes (Mueser et al., 1991; Bowie et al., 2008); however, to the best of our knowledge, no study has yet investigated how accurately such electrophysiological biomarkers can predict symptom severity of each individual patient.
In this preliminary study, we investigated whether it is possible to accurately predict the symptom severity of patients with schizophrenia from electrophysiological measurements. We used event-related potential (ERP) and its source imaging results as candidate variables for the estimation of symptom severity of patients with schizophrenia. Electroencephalography (EEG) was recorded while patients performed a facial affect discrimination task that was introduced in our previous study (Kim et al., 2013). We found that patients with schizophrenia not only showed significantly decreased neuronal activations during facial affective processing compared to normal controls, but also showed significant correlations between such neuronal activations and Positive and Negative Syndrome Scale (PANSS) scores (Kim et al., 2013). Based on these findings, we established a mathematical prediction model based on electrophysiological biomarkers to estimate the Positive/Negative PANSS scores of individual patients.
Materials and Methods
We recruited a total of 23 patients with schizophrenia for the current study. All participants were stable, right-handed, with normal or corrected-to-normal vision. The mean age of the participants was 32.2 ± 10.1 (mean ± SD) years, and 11 of them were female (Table 1). All participants had been diagnosed with schizophrenia based on the Structured Clinical Interview for Diagnostic and Statistical Manual of Mental Disorders, 4th Edition (DSM-IV), Axis I Psychiatric Disorders (Kay et al., 1987). All subjects were taking atypical antipsychotics (olanzapine, n = 11; risperidone, n = 12). This study was carried out in accordance with the recommendations of the Institutional Review Board (IRB) of Inje University Ilsan Paik Hospital with written informed consent from all subjects. All subjects gave written informed consent in accordance with the Declaration of Helsinki. The protocol was approved by the IRB of Inje University Ilsan Paik Hospital.
The assessment of the psychiatric symptom of the patients was done by a trained psychiatrist using the PANSS. The PANSS evaluates the severity of the two common types of symptoms (positive and negative) in schizophrenia as well as the general psychopathology of the patient based on the interview as well as reports of family members. The positive scale is rated by seven positive symptoms such as hallucination and delusion, which refer to an excess or distortion of normal functions. The negative scale is also rated by seven items of negative symptoms, which show loss or reduced functions compared to healthy subjects (e.g., emotional withdrawal). Each item is rated from 1 (absent) to 7 (extreme) and summed up for each category, resulting in a scale of 7–49 for each symptom scale. Additionally, the general scale of the patient is rated by 16 items that evaluates the general psychopathological symptoms, such as anxiety, tension, and poor attention. The average positive and negative score of the patients recruited for this study was 20.2 ± 7.8 and 18.7 ± 7.4, respectively (Figure 1).
Neural source activations were evaluated using EEG data recorded from 64 electrodes including two electrooculography channels (VEOG and HEOG). A total of 288 facial stimuli were presented to each patient in a randomized order, with an equal presentation probability with respect to three different emotional conditions (neutral, fear, and happy) (Lee et al., 2004). In this study, we only used EEG data from the neutral face condition in order to reduce the dimensionality of features because a greater number of candidate features does not always guarantee more precise prediction results.
As the first step of the analysis, we obtained source activation of distinct brain areas which has a significant correlation with the symptom score (Figure 2A) The overall analysis procedure is illustrated in Figure 2. The recorded EEG signals were preprocessed to eliminate unwanted artifacts and averaged to identify four ERP components known to be associated with facial emotion processing: P100, N170, N250, and P300 (Lee et al., 2010; Kim et al., 2013). The source image for each ERP component was reconstructed using sLORETA (Pascual-Marqui, 2002) for each patient. The sLORETA is one of the most widely-used source imaging methods that showed reliable source estimates even in the presence of noise or in low-density EEG configurations (Babiloni et al., 2010; Bae et al., 2011; Saavedra et al., 2012). In our previous study that used the same patient data, we have shown that source activation of N170 estimated using sLORETA was decreased in patients with schizophrenia compared to normal controls in multiple brain areas (Jung et al., 2012). After the execution of sLORETA software, voxels with a significant correlation with symptom scores were identified (Table 2). The significant voxels were then clustered based on anatomical location, and the voxel with the maximum correlation within each cluster was selected. The detailed procedures for data acquisition and data preprocessing are well described in our previous study (Kim et al., 2013).
Figure 2. An illustration of the overall analysis procedure: (A) Using the identified face-related ERP components (P100, N170, N250, P300), the source activation of each ERP was estimated using sLORETA. Multiple nearby voxels showing significant correlation between source activation and symptom severity scores were clustered. (B) The voxel activation with the highest correlation within each cluster was used as an independent variable of the general linear model to estimate the symptom severity score. The generated model was validated using leave-one-out cross validation (LOOCV).
Table 2. Brain regions showing significant correlation between PANSS scores and ERP source activation during the neutral face stimulus condition.
To estimate each patient's positive and negative symptom severity, we established a prediction model based on the general linear model (GLM) (Henderson et al., 2001; Coan and Allen, 2004). Each patient's ERP amplitudes and the source activation of the voxel that had the maximum correlation value in each cluster were selected as independent variables in the regression model (Figure 2B). The source activation was the standardized score obtained from sLORETA. The Positive and Negative scales of the PANSS were used as dependent variables in the regression model. To exclude unnecessary or relatively less relevant variables, a stepwise selection method was applied with entrance/exit tolerances of 0.05/0.10 (p-values).
The prediction model of the Positive scale of the PANSS was composed of two predictors [source activations at the medial frontal gyrus (N250) and the precentral gyrus (P100)]. The weighted linear combination of these predictors was significant [F(2, 22) = 35.402, p < 0.001] and explained approximately 78% of the variance in the Positive scale (R = 0.883, R2 = 0.780). A prediction model of the Negative scale of the PANSS was constructed using the source activations of the sub-gyral (P100) and middle frontal gyrus (N250), which was also significant [F(2, 22) = 17.507, p < 0.001] and explained approximately 63% of the variance in the Negative scale (R = 0.798, R2 = 0.636). The selected predictor variables of each model are listed in Table 3. Note that we also considered ERP amplitudes as candidate predictor variables, but none of them was included in the final prediction models as the ERP amplitudes were not selected by the stepwise selection method with the entrance tolerance of p = 0.05 at all. Even when we changed the entrance tolerance to p = 0.5 to construct a regression model only with ERP amplitudes, the best prediction models to estimate positive and negative PANSS scores could explain only 20.8 and 13.9% of the variances in each score, respectively. Likewise, we also considered age of patients as an independent variable of the model, but it was not included in the final prediction model.
To evaluate the performance of the prediction models, we adopted a leave-one-out cross-validation (LOOCV) strategy. LOOCV is the most extreme and accurate type of cross-validation strategy to estimate prediction errors where every data sample is used in turn for the validation of the model constructed with the remaining data (Hastie et al., 2003; Ko et al., 2011). LOOCV results showed that the symptom scores could be predicted with an average error of 3.34 ± 2.40 and 3.90 ± 3.01 for the Positive and Negative scales, respectively (Table 3).
For comparison, each patient's Positive and Negative scale scores were assumed to have either the average or median of the group symptom scores. Mean errors were 6.43 ± 4.18 (Positive) and 6.03 ± 4.05 (Negative) when the average symptom scores were used, and 6.09 ± 5.48 (Positive) and 5.96 ± 4.36 (Negative) when the median symptom scores were used, both of which were much higher than the scores estimated with our prediction models.
In this study, we investigated whether the symptom severity scores of individual patients with schizophrenia could be estimated using ERP current source activations during facial affect processing. We found that Positive and Negative scale scores can be estimated with a fairly high accuracy considering the wide range of each score (0~49), demonstrating the possibility of a quantitative assessment of psychiatric symptoms based on electrophysiological biomarkers. To the best of our knowledge, the current preliminary study is the first attempt to quantitatively estimate clinical symptom severity of schizophrenia from electrophysiological biomarkers.
The evaluation of a patient's psychiatric symptom scores using quantitative biomarkers provides a raw estimate of the patient's symptom severity without the need for an interview with psychiatrists or clinical psychologists. Specifically, this technique would be useful in estimating the symptom scores of patients who have low cognitive profiles or demonstrate unreliable behavior. Further, this technique can be used as an auxiliary tool for long-term longitudinal tracking of symptom changes in order to evaluate the responsiveness of psychopharmacological treatment and to establish an appropriate patient-specific treatment strategy.
In our study results, both Positive and Negative scales could be successfully estimated because abnormal facial affect processing is regarded as one of the key symptoms of schizophrenia associated with both positive and negative symptoms (Mandal et al., 1999; Kohler et al., 2003; Hofer et al., 2009). This suggests that the selection of proper symptom-related tasks is important in building reliable prediction models. For instance, we expected that an altered P300 amplitude and its source activation would be good candidate predictors of symptom scores because they have also been frequently reported to be strongly correlated with symptomatic scores of patients with schizophrenia (Shenton et al., 1989; Strik et al., 1993; Kim et al., 2014). In addition, other (electro-) physiological biomarkers (e.g., graph theory-based indices or behavioral outcomes) may potentially improve the overall estimation accuracy. Therefore, further studies are needed to compare the performance of various tasks and/or biomarkers in terms of estimation accuracy. We expect that our approach to building a prediction model to estimate symptom severity in schizophrenia from electrophysiological biomarkers can also be applied to other psychiatric diseases such as Alzheimer's disease, depression, and bipolar disorder, using specific paradigms reflecting each disease's unique signs or symptoms.
DK was responsible for the analyzing procedure wrote the manuscript. SL designed the study and wrote the protocol. MS analyzed the EEG data and produced ERP waveforms. CI supervised the study process and manuscript writing. All authors contributed to and have approved the final manuscript.
This work was supported by National Research Foundation of Korea (NRF) grants funded by the Korean government (MSIP) (2015R1A2A1A15054662 and NRF-2015M3C7A1031969).
Conflict of Interest Statement
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.
Babiloni, C., Vecchio, F., Buffo, P., Buttiglione, M., Cibelli, G., and Rossini, P. M. (2010). Cortical responses to consciousness of schematic emotional facial expressions: a highresolution EEG study. Hum. Brain Mapp. 31, 1556–1569. doi: 10.1002/hbm.20958
Bae, K.-Y., Kim, D.-W., Im, C.-H., and Lee, S.-H. (2011). Source imaging of P300 auditory evoked potentials and clinical correlations in patients with posttraumatic stress disorder. Prog. Neuro-psychopharmacol. Biol. Psychiatry 35, 1908–1917. doi: 10.1016/j.pnpbp.2011.08.002
Bell, M. D., Lysaker, P. H., Milstein, R. M., and Beamgoulet, J. L. (1994). Concurrent validity of the cognitive component of schizophrenia - relationship of panss scores to neuropsychological assessments. Psychiatry Res. 54, 51–58. doi: 10.1016/0165-1781(94)90064-7
Birchwood, M., Smith, J., Macmillan, F., Hogg, B., Prasad, R., Harvey, C., et al. (1989). Predicting relapse in schizophrenia - the development and implementation of an early signs monitoring-system using patients and families as observers, a preliminary investigation. Psychol. Med. 19, 649–656. doi: 10.1017/S0033291700024247
Bowie, C. R., Leung, W. W., Reichenberg, A., McClure, M. M., Patterson, T. L., Heaton, R. K., et al. (2008). Predicting schizophrenia patients' real-world behavior with specific neuropsychological and functional capacity measures. Biol. Psychiatry 63, 505–511. doi: 10.1016/j.biopsych.2007.05.022
Bow-Thomas, C. C., Velligan, D. I., Miller, A. L., and Olsen, J. (1999). Predicting quality of life from symptomatology in schizophrenia at exacerbation and stabilization. Psychiatry Res. 86, 131–142. doi: 10.1016/S0165-1781(99)00023-2
Henderson, H. A., Fox, N. A., and Rubin, K. H. (2001). Temperamental contributions to social behavior: the moderating roles of frontal EEG asymmetry and gender. J. Am. Acad. Child Adolesc. Psychiatry 40, 68–74. doi: 10.1097/00004583-200101000-00018
Hirayasu, Y., Potts, G. F., O'Donnell, B. F., Kwon, J. S., Arakaki, H., Akdag, S. J., et al. (1998). Auditory mismatch negativity in schizophrenia: topographic evaluation with a high-density recording montage. Am. J. Psychiatry 155, 1281–1284. doi: 10.1176/ajp.155.9.1281
Ho, B. C., Psych, M. R. C., Nopoulos, P., Flaum, M., Arndt, S., and Andreasen, N. C. (1998). Two-year outcome in first-episode schizophrenia: predictive value of symptoms for quality of life. Am. J. Psychiatry 155, 1196–1201. doi: 10.1176/ajp.155.9.1196
Hofer, A., Benecke, C., Edlinger, M., Huber, R., Kemmler, G., Rettenbacher, M. A., et al. (2009). Facial emotion recognition and its relationship to symptomatic, subjective, and functional outcomes in outpatients with chronic schizophrenia. Eur. Psychiatry 24, 27–32. doi: 10.1016/j.eurpsy.2008.06.008
Javitt, D. C., Shelley, A., and Ritter, W. (2000). Associated deficits in mismatch negativity generation and tone matching in schizophrenia. Clin. Neurophysiol. 111, 1733–1737. doi: 10.1016/S1388-2457(00)00377-1
Jung, H. T., Kim, D. W., Kim, S., Im, C. H., and Lee, S. H. (2012). Reduced source activity of event-related potentials for affective facial pictures in schizophrenia patients. Schizophr. Res. 136, 150–159. doi: 10.1016/j.schres.2011.10.023
Kim, D. W., Kim, H. S., Lee, S. H., and Im, C. H. (2013). Positive and negative symptom scores are correlated with activation in different brain regions during facial emotion perception in schizophrenia patients: a voxel-based sLORETA source activity study. Schizophr. Res. 151, 165–174. doi: 10.1016/j.schres.2013.10.025
Kim, D. W., Shim, M., Kim, J. I., Im, C. H., and Lee, S. H. (2014). Source activation of p300 correlates with negative symptom severity in patients with schizophrenia. Brain Topogr. 27, 307–317. doi: 10.1007/s10548-013-0306-x
Ko, L.-W., Wei, C.-S., Jung, T.-P., and Lin, C.-T. (2011). “Estimating the level of motion sickness based on EEG spectra,” in Foundations of Augmented Cognition. Directing the Future of Adaptive Systems, eds D. Schmorrow and C. Fidopiastis (Berlin; Heidelberg: Springer), 169–176.
Kohler, C. G., Turner, T. H., Bilker, W. B., Brensinger, C. M., Siegel, S. J., Kanes, S. J., et al. (2003). Facial emotion recognition in schizophrenia: intensity effects and error pattern. Am. J. Psychiatry 160, 1768–1774. doi: 10.1176/appi.ajp.160.10.1768
Lee, S. H., Kim, E. Y., Kim, S., and Bae, S. M. (2010). Event-related potential patterns and gender effects underlying facial affect processing in schizophrenia patients. Neurosci. Res. 67, 172–180. doi: 10.1016/j.neures.2010.03.001
Lee, W. H., Chae, J. H., Bahk, W. M., and Lee, K. U. (2004). Development and its preliminary standardization of pictures of facial expressions for affective neurosciences. J. Korean Neuropsychiatr Assoc. 43, 552–558.
Lindström, E., Wieselgren, I. M., and von Knorring, L. (1994). Interrater reliability of the structured clinical interview for the positive and negative syndrome scale for schizophrenia. Acta Psychiatr. Scand. 89, 192–195. doi: 10.1111/j.1600-0447.1994.tb08091.x
Mandal, M. K., Jain, A., Haque-Nizamie, S., Weiss, U., and Schneider, F. (1999). Generality and specificity of emotion-recognition deficit in schizophrenic patients with positive and negative symptoms. Psychiatry Res. 87, 39–46. doi: 10.1016/S0165-1781(99)00047-5
Mueser, K. T., Bellack, A. S., Douglas, M. S., and Wade, J. H. (1991). Prediction of social skill acquisition in schizophrenic and major affective-disorder patients from memory and symptomatology. Psychiatry Res. 37, 281–296. doi: 10.1016/0165-1781(91)90064-V
Norman, R. M. G., Malla, A. K., Cortese, L., and Diaz, F. (1996). A study of the interrelationship between and comparative interrater reliability of the SAPS, SANS and PANSS. Schizophr. Res. 19, 73–85. doi: 10.1016/0920-9964(95)00055-0
Ostermann, J., Uhl, I., Kohler, E., Juckel, G., and Norra, C. (2012). The loudness dependence of auditory evoked potentials and effects of psychopathology and psychopharmacotherapy in psychiatric inpatients. Hum. Psychopharmacol. 27, 595–604. doi: 10.1002/hup.2269
Saavedra, C., Iglesias, J., and Olivares, E. I. (2012). Event-related potentials elicited by face identity processing in elderly adults with cognitive impairment. Exp. Aging Res. 38, 220–245. doi: 10.1080/0361073X.2012.660057
Sabeti, M., Boostani, R., Katebi, S. D., and Price, G. W. (2007). Selection of relevant features for EEG signal classification of schizophrenic patients. Biomed. Signal Process. Control 2, 122–134. doi: 10.1016/j.bspc.2007.03.003
Sabeti, M., Katebi, S., and Boostani, R. (2009). Entropy and complexity measures for EEG signal classification of schizophrenic and control participants. Artif. Intell. Med. 47, 263–274. doi: 10.1016/j.artmed.2009.03.003
Shenton, M. E., Faux, S. F., McCarley, R. W., Ballinger, R., Coleman, M., Torello, M., et al. (1989). Correlations between abnormal auditory-P300 topography and positive symptoms in schizophrenia - a preliminary-report. Biol. Psychiatry 25, 710–716. doi: 10.1016/0006-3223(89)90242-4
Spencer, K. M., Salisbury, D. F., Shenton, M. E., and McCarley, R. W. (2008). Gamma-band auditory steady-state responses are impaired in first episode psychosis. Biol. Psychiatry 64, 369–375. doi: 10.1016/j.biopsych.2008.02.021
Strik, W. K., Dierks, T., and Maurer, K. (1993). Amplitudes of auditory P300 in remitted and residual schizophrenics - correlations with clinical-features. Neuropsychobiology 27, 54–60. doi: 10.1159/000118953
Takayanagi, Y., Takahashi, T., Orikabe, L., Mozue, Y., Kawasaki, Y., Nakamura, K., et al. (2011). Classification of first-episode schizophrenia patients and healthy subjects by automated MRI measures of regional brain volume and cortical thickness. PLoS ONE 6:e21047. doi: 10.1371/journal.pone.0021047
Wyss, C., Hitz, K., Hengartner, M. P., Theodoridou, A., Obermann, C., Uhl, I., et al. (2013). The loudness dependence of auditory evoked potentials (LDAEP) as an indicator of serotonergic dysfunction in patients with predominant schizophrenic negative symptoms. PLoS ONE 8:e68650. doi: 10.1371/journal.pone.0068650
Youn, T., Park, H. J., Kim, J. J., Kim, M. S., and Kwon, J. S. (2003). Altered hemispheric asymmetry and positive symptoms in schizophrenia: equivalent current dipole of auditory mismatch negativity. Schizophr. Res. 59, 253–260. doi: 10.1016/S0920-9964(02)00154-8
Keywords: schizophrenia, electroencephalogram (EEG), generalized linear model (GLM), electrophysiological biomarker, psychiatric diseases
Citation: Kim D-W, Lee S-H, Shim M and Im C-H (2017) Estimation of Symptom Severity Scores for Patients with Schizophrenia Using ERP Source Activations during a Facial Affect Discrimination Task. Front. Neurosci. 11:436. doi: 10.3389/fnins.2017.00436
Received: 20 April 2017; Accepted: 14 July 2017;
Published: 03 August 2017.
Edited by:Sung Chan Jun, Gwangju Institute of Science and Technology, South Korea
Reviewed by:Yingchun Zhang, University of Houston, United States
Elisabetta C. Del Re, Harvard Medical School, United States
Copyright © 2017 Kim, Lee, Shim and Im. 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) or licensor 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: Chang-Hwan Im, firstname.lastname@example.org