Discriminant analysis of functional optical topography for schizophrenia diagnosis

Abstract. Abnormal prefrontal function plays a central role in the cognition deficits of schizophrenic patients; however, the character of the relationship between discriminant analysis and prefrontal activation remains undetermined. Recently, evidence of low prefrontal cortex (PFC) activation in individuals with schizophrenia has also been found during verbal fluency tests (VFT) and other cognitive tests with several neuroimaging methods. The purpose of this study is to assess the hemodynamic changes of the PFC and discriminant analysis between schizophrenia patients and healthy controls during VFT task by utilizing functional optical topography. A total of 99 subjects including 53 schizophrenic patients and 46 age- and gender-matched healthy controls were studied. The results showed that the healthy group had larger activation in the right and left PFC than in the middle PFC. Besides, the schizophrenic group showed weaker task performance and lower activation in the whole PFC than the healthy group. The result of the discriminant analysis showed a significant difference with P value <0.001 in six channels (CH 23, 29, 31, 40, 42, 52) between the schizophrenic and healthy groups. Finally, 68.69% and 71.72% of subjects are correctly classified as being schizophrenic or healthy with all 52 channels and six significantly different channels, respectively. Our findings suggest that the left PFC can be a feature region for discriminant analysis of schizophrenic diagnosis.


Introduction
Schizophrenia is an irreversible neurodevelopmental disorder with severe individual, family, and societal burdens. 1 Several studies have indicated that schizophrenic patients showed impaired performance in various aspects of social cognition, including theory of mind, emotion processing, and agency judgments. [2][3][4][5] The neuroimaging techniques such as functional magnetic resonance imaging (fMRI) and positron emission tomography (PET) have been widely used to image brain functions for many neurodegenerative diseases. [6][7][8][9][10][11][12] These techniques have greatly increased our knowledge about the neuropathological process of schizophrenia based on cognitive and emotional processes. Although these methods offer noninvasive imaging with excellent spatial resolution, the limitation of huge size of instruments and confinement of the participants to restricted positions inside the magnet cannot be applied for patientoriented diagnosis. Also, PET requires the injection of radioactive materials. These drawbacks make the imaging modalities difficult or impossible to apply in clinical diagnosis for neonates, children, old people, and claustrophobic patients.
Functional optical topography (fOT), imaging based on functional near-infrared spectroscopy measurement, is an another noninvasive neuroimaging technique that provides continuous recording of brain hemodynamic activity with detection of oxy-and deoxy-hemoglobin concentration changes. There are several benefits such as lower cost, nonionizing radiation imaging, real-time measurement, long time monitoring, and easy operation. 13 Although spatial resolution is limited, fOT offers more comprehensive information of brain activity than bloodoxygenation-level-dependent signal of fMRI. 14 Furthermore, the optical method can provide completely patient-oriented measurement. The fOT is also applicable for psychological test, because its temporal resolution is high enough to detect the changes of short duration such as the brain activation during cognitive task. In previous studies, fOT method has been shown to be sensitive enough for monitoring of physiological blood oxygenation changes during cognitive activation in schizophrenia. [15][16][17] Recently, fOT is being increasingly applied for investigation of prefrontal cortex (PFC) dysfunction in schizophrenia, because the social cognition of schizophrenic patients was associated with PFC activation. PFC dysfunction in individuals with schizophrenia has been detected and evidenced during verbal fluency test (VFT) and other cognitive tests by utilizing fOT. [18][19][20][21][22] However, the previous studies only show the group averaged data and correlational analysis between fOT and clinical symptoms. There are few studies reporting significantly different regions of PFC and discriminant analysis between schizophrenic and healthy subjects based on fOT measurement. The significantly different regions of PFC could be an indicator to evaluate the PFC function of schizophrenia in patient-oriented diagnosis. Thus, the goal of this study is to characterize/analyze fOT image between schizophrenia patients and healthy controls during VFT.
In this paper, we perform statistical analyses of PFC fOT image with VFT. The imaging results between schizophrenia and healthy controls are significantly different. Furthermore, the left PFC shows significant difference between schizophrenic and healthy subjects under the 0.1% significance level. In discriminant analysis, the 68.69% and 71.72% are correctly classified as being schizophrenic or healthy subjects with all 52 channels and six significantly different channels over the left PFC, respectively. Our results imply that the left PFC fOT imaging could be a feature region for schizophrenia-aided diagnosis.

Participants
The study includes 53 schizophrenic patients [female∕male: 27∕26; mean age: 34.4; standard deviation (SD): 11] and 46 age-and sex-matched healthy controls (female∕male: 25∕21; mean age 35.1; SD 11.2), who were recruited from the Tottori University Hospital in Yonago, Japan. The diagnosis for schizophrenia was based on the Structured Clinical Interview for Diagnostic and Statistical Manual of Mental Disorders (DSM)-IV. 23 Healthy controls had no history of psychiatric or neurologic disorders. All subjects were native Japanese speakers. There was no difference between schizophrenic and control groups with respect to age, gender, and body mass index (BMI). However, the education, estimated premorbid IQ, and number of words generated show significant differences. Schizophrenic symptoms were evaluated with the Positive and Negative Syndrome Scale (PANSS) and Global Assessment of Functioning. [24][25][26][27] Table 1 shows the clinical characteristics of the schizophrenic and healthy groups. All participants gave written informed consent. The study was in accordance with the latest version of the Declaration of Helsinki, and approved by the Ethics Committee of Tottori University Faculty of Medicine.

Functional Optical Topography
The fOT measurements were conducted with a 52-channel ETG-4000 Optical Topography System (Hitachi Medical Co., Tokyo, Japan). The near-infrared sources at dual-wavelength 695 and 830 nm were used, and the temporal resolution of detection is 0.1 s. According to the modified Beer-Lambert law, 28,29 the concentration changes in oxygenated hemoglobin (oxy-Hb) and deoxygenated hemoglobin (deoxy-Hb) could be obtained.
Previous studies indicate that [oxy-Hb] is the most sensitive indicator of changes in regional cerebral blood flow in fOT measurement. 30,31 In fOT system, 17 emitters and 16 detectors are plugged into a probing holder and arranged as 3 × 11 array for 52-channel detection with 30-mm source-detector separation. In order to cover the whole PFC, the array was positioned according to the international 10-20 system. 32 Specifically, the detector 13 was placed in F z ; the midline of channels nos. 46 and 49 were placed in the F p2 and F p1 , and the source 23 and 28 were fitted around T 4 and T 3 , respectively (as shown in Fig. 1).

Data Acquisition: VFT Tasks
The subjects were seated in a comfortable chair in a silent room and kept eyes open during the fOT monitoring. The VFT consists of a 60-s long VFT (letter version), a 30-s long pretask baseline, and a 70-s long posttask baseline. The three syllables (I: /to/, /se/ or /o/; II: /na/, /i/ or /ta/; III: /a/, /ki/ or /ha/) were announced in the order which was counterbalanced across participants. The fOT data were obtained by subtracting the baseline [oxy-Hb] while the subjects were simply repeating "a," "i," "u," "e," and "o" loudly. The VFT requires subject to say the words beginning with a specific syllable as many as possible. The schematic diagram of the task is shown in Fig. 2.   Figure 3 shows the flowchart of fOT data analyses. In each subject, the mean levels in [oxy-Hb] of all 52 channels are calculated from both baseline and activation task of VFT. In order to find the significantly different regions of PFC, the fOT data are analyzed based on statistical tests. We adopt a chi-square goodness-of-fit test for statistical hypothesis testing. The chi-square test-statistic can be expressed as 33,34

Data Analyses
where O i is an observed count and E i is an expected count; χ 2 is Pearson's cumulative test statistic; and n is the number of cells in the Due to physiological noise, the majority of the channels were found to be nonnormally distributed, and therefore, we used the nonparametric test on all the data. The two-sample Kolmogorov-Smirnov test (K-S test) is one of the most powerful nonparametric methods with highly sensitive comparison in both location and shape of the empirical cumulative distribution functions of the two samples. 35,36 The two-sample K-S test is expressed as D n;n 0 ¼ SUP x jG n ðxÞ − F n 0 ðxÞj; where SUP x is the supremum of the set of distances; and G n ðxÞ and F n 0 ðxÞ are the empirical distribution functions of the first and the second variables, respectively. The null hypothesis is rejected at significance level α if ffiffiffiffiffiffiffiffiffiffiffiffiffi nn 0 n þ n 0 r D n;n 0 > K α : (3) Equation (3) implies G n ðxÞ and F n 0 ðxÞ are significantly different. The test statistic K α is the maximum difference between the curves. After K-S testing, the significantly different regions of PFC are defined between healthy and schizophrenic groups. Then, the discriminant analysis is applied by leave-one-out cross-validation with k-means clustering. 37 In statistics and data mining, the k-means method is one of the most widely used clustering algorithms with nearest mean, drawing its popularity from its speed in practice. 37,38 Leave-one-out crossvalidation involves single observation from the original sample as the validation data, and the remaining observations as the training data. The validation step is looped until all observations are used once as the validation data. In this study, one subject is selected from all 99 subjects in sequence. Thus, the other 98 subjects are treated as training data with k-means clustering. The 99 times clustering of test data is then used for accuracy estimation of discriminant analysis. Figure 4 shows the [oxy-Hb] mapping of group-level statistics on averaging of healthy controls and schizophrenic patients. The [oxy-Hb] concentration of schizophrenic patients is significantly smaller than healthy controls during VFT.

Results
In case of healthy controls, the [oxy-Hb] response is stronger on both sides than in middle PFC. For quantitative analysis, the mean of [oxy-Hb] on right, middle, and left PFC during VFT are demonstrated in Fig. 5. There exists the significant difference between healthy and schizophrenic groups (paired t-test, *P < 0.005, **P < 0.001). Figure 6 shows the [oxy-Hb] mapping via 52 channels comparisons based on two-sample K-S test. The CH 23, 29, 31, 40, 42, and 52 are significantly different with P < 0.001. We can observe that the channel on the left regions of PFC is significantly different. This result implies that the left PFC may provide significant feature region to distinguish schizophrenia from healthy subjects in clinical diagnosis. The discriminant analysis reveals the significant difference with K-S test by utilizing leave-one-out cross-validation with k-means clustering. The four possible outcomes from a binary classifier of 52 and six Fig. 3 The flowchart of data analyses with fOT. The two-sample K-S test of nonparametric methods was used to find the significantly different channels because some channels in healthy group and schizophrenic group are not normally distributed by using the chi-square goodness-offit test. Finally, the leave-one-out cross-validation with k-means clustering was used for discriminant analysis from healthy to schizophrenia group. dimensions are shown in Table 2. The accuracy of correctly classified with 52 dimensions was 68.69% with a sensitivity of 85% and a specificity of 50% (i.e., true positive ¼ 45 of the 53 schizophrenic patients and true negative ¼ 23 of the 46 healthy controls). In the case of six significant difference dimensions, the accuracy of correctly classified was 71.72% with a sensitivity of 77% and a specificity of 65% (i.e., true positive ¼ 41 of the 53 schizophrenic patients and true negative ¼ 30 of the 46 healthy controls).

Discussion and Conclusions
The [oxy-Hb] responses on both sides of PFC are stronger than middle PFC from healthy controls, which implies deeper light penetration in the interhemispheric fissure that causes the fewer backscattered photons to be received. Furthermore, the [oxy-Hb] mapping of schizophrenic case is significantly weaker than normal on PFC and no significant difference is found among right, left, and middle forehead areas. It indicates the task-dependent profile of functional abnormalities can be observed from the brain metabolism on schizophrenic PFC. Also, the fOT measurement may offer a good potential to detect the structural change of brain. To our knowledge, we are the first to report the significant difference regions of PFC by utilizing discriminant analysis between schizophrenic and healthy subjects based on fOT monitoring. Although Azechi et al. used prefrontal activation with frontal lobe tasks for discriminant analyzing between schizophrenic and healthy subjects, 39 they did not provide the significant difference regions of PFC for discriminating schizophrenia. In this study, we pointed out the significant difference regions of PFC between schizophrenia and healthy controls. Figure 2 shows the group-level of [oxy-Hb] in right, left, and middle PFC that is significantly different between schizophrenia and healthy controls. According to two-sample K-S test, the channels of significant difference are analyzed on left PFC. Abnormal structural and/or functional hemisphere asymmetries in schizophrenia disorders have been assessed based on fOT. 18,40 In an fMRI study, reduced activation of the left frontal regions during word production was found in schizophrenia patients compared with normal. 41 The other study shows that schizophrenia produces a failure of attentional modulation that leads to a breakdown in the selective enhancement or inhibition of semantic/lexical representations, whose biological substrata are widely distributed across left (dominant) temporal and frontal lobes. 42 These findings can be interpreted as a sign of reduced specific lateralized PFC reactivity, possibly based on a left hemisphere functional deficit between schizophrenia and normal. Our result of the significant difference regions on left PFC is in agreement with previous studies. Therefore, the left PFC shows a great potential of feature region for discriminant analysis of schizophrenia. Table 2 demonstrates the accuracy of classification by applying six significant channels, which is better than 52 channels. Therefore, the significant channels in fOT can be used for schizophrenia diagnosis with discriminant analysis. However, the limitation of this study must be taken into account. Discriminant analysis with fOT can reveal significant differences between the probability distributions of signals in healthy and diseased groups, but cannot provide an individual diagnosis. Nevertheless, optical imaging with fOT can still be used as an auxiliary tool in schizophrenia diagnosis. Fig. 6 The result of two-sample K-S test. The CH 23, 29, 31, 40, 42, and 52 are significantly different with P < 0.001 between healthy and schizophrenia groups. This result further indicates that the activation of left PFC during VFT is significantly different between healthy and schizophrenia group.  In conclusion, the [oxy-Hb] PFC mapping of schizophrenia cases is weaker than healthy ones during VFT. The low [oxy-Hb] change of cognitive activation is concomitant with PFC dysfunction of schizophrenia. The significantly different channels of left PFC provide high-accuracy discrimination for schizophrenia. fOT offers both structural and functional information simultaneously with patient-oriented diagnosis. Furthermore, it could be used for the neural circuitry abnormality detection and estimation of pharmacological treatment and/ or cognitive rehabilitation.