ORIGINAL RESEARCH article
Classification of Schizophrenia by Functional Connectivity Strength Using Functional Near Infrared Spectroscopy
- 1China Academy of Information and Communications Technology, Beijing, China
- 2Institute of Electrical Engineering, Chinese Academy of Sciences, Beijing, China
- 3Peking University Sixth Hospital, Peking University Institute of Mental Health, NHC Key Laboratory of Mental Health (Peking University), National Clinical Research Center for Mental Disorders (Peking University Sixth Hospital), Beijing, China
- 4School of Computer Science and Technology, Donghua University, Shanghai, China
- 5School of Computer Science and Technology, Xidian University, Xian, China
Functional near-infrared spectroscopy (fNIRS) has been widely employed in the objective diagnosis of patients with schizophrenia during a verbal fluency task (VFT). Most of the available methods depended on the time-domain features extracted from the data of single or multiple channels. The present study proposed an alternative method based on the functional connectivity strength (FCS) derived from an individual channel. The data measured 100 patients with schizophrenia and 100 healthy controls, who were used to train the classifiers and to evaluate their performance. Different classifiers were evaluated, and support machine vector achieved the best performance. In order to reduce the dimensional complexity of the feature domain, principal component analysis (PCA) was applied. The classification results by using an individual channel, a combination of several channels, and 52 ensemble channels with and without the dimensional reduced technique were compared. It provided a new approach to identify schizophrenia, improving the objective diagnosis of this mental disorder. FCS from three channels on the medial prefrontal and left ventrolateral prefrontal cortices rendered accuracy as high as 84.67%, sensitivity at 92.00%, and specificity at 70%. The neurophysiological significance of the change at these regions was consistence with the major syndromes of schizophrenia.
Schizophrenia is a kind of psychiatric disorder characterized by a series of positive/psychotic (e.g., hallucinations and delusions), negative/deficit (e.g., insufficiency of thought and loss of motivation), and cognitive (e.g., impairment of memory and attention) symptoms (Buckley et al., 2009). Conventionally, clinical diagnostic criteria are predominately based on the relative subjective approaches, for example, according to the diagnostic manuals (American Psychiatric Association, 1994). With the development of neuroimaging, a number of objective methods to identify schizophrenia patients have emerged, e.g., single photon emission computed tomography (SPECT; Gordon et al., 1994), diffusion tensor imaging (DTI; Ohtani et al., 2014), functional magnetic resonance imaging (fMRI; Weiss et al., 2004; Deng et al., 2017; Tréhout et al., 2017), and functional near infrared spectroscopy (fNIRS; Kubota et al., 2005; Rosenbaum et al., 2017).
fNIRS is a noninvasive hemodynamic imaging technique used to assess functional activities in the human brain (Hoshi, 2003). It detects the concentration of oxygenated hemoglobin (Oxy-Hb) and deoxygenated hemoglobin (Deoxy-Hb) by measuring the absorption and reflection of specific near infrared spectrums in the cortices during tasks. Compared to other neuroimaging instruments, fNIRS has the benefit of being low cost with a high portability. These advantages have enabled its application in the diagnosis of schizophrenia, which was mainly based on the effect of hypofrontality (reduced activation around the bilateral prefrontal cortex) during various verbal fluency tasks (VFTs; Suto et al., 2004; Ehlis et al., 2007; Takizawa et al., 2008; Ji et al., 2020). In practice, the majority of these studies extracted the time-domain features from single or multiple channels of healthy subjects and patients with schizophrenia (Suto et al., 2004; Kanahara et al., 2013; Sugimura et al., 2014; Tian et al., 2019). As a consequence, diverse machine learning classifiers (Li et al., 2015) were trained and distinguished the patients with schizophrenia from the healthy subjects.
In contrast, substantial neuroimaging studies of other modalities have found abnormal dysconnectivity between the prefrontal cortex and temporal cortex in schizophrenia patients (Friston and Frith, 1995; Maguire et al., 2000; Greicius, 2008; Bullmore and Sporns, 2009; Whitfield-Gabrieli et al., 2009), and these experiments using EEG and fMRI have proposed classification methods based on brain network properties (Demirci et al., 2008; Yang et al., 2010; Arbabshirani et al., 2013). Nevertheless, the method on whole-brain network properties cannot be directly applied to fNIRS analysis because the conventional clinical fNIRS only measures signals from the frontotemporal cortex. One fNIRS study discriminated patients with schizophrenia using four global network properties (Song et al., 2016). The achieved overall accuracy was 85.5%, but the local changes could not be investigated with the approach. Hence, the analysis on regional functional connectivity (FC), integrating both the spatial and temporal relation of brain activities, is hypothesized to provide new insights on classifying schizophrenia.
In this article, we provided an FC-based method to identify schizophrenia patients. Oxy-Hb data from 100 schizophrenia patients and 100 healthy subjects during VFT were used in the experiment. functional connectivity strength (FCS) from single channel, from the ensemble 52 channels, from the dimensional reduced 52 channels, and from different combinations of 2–5 channels were used to trained four popular classifiers (Linear Discriminant Analysis: LDA, k-Nearest Neighbor: KNN, Gaussian Processes classifier: GPC, and Support Vector Machine: SVM), respectively. The best accuracy was 85.00% (LOOCV), with sensitivity as 87.00% (LOOCV) and specificity as 83.00% (LOOCV), by using FCS from three channels. Theneurophysiological significance was discussed. The FCS-based method provided a new and effective approach for schizophrenic identification.
Materials and Methods
The Oxy-Hb dataset included 100 schizophrenic (male/female: 50/50, 33.81 ± 11.52 years old and ranging from 18 to 53 years old) and 100 healthy subjects (male/female: 47/53, 34.43 ± 12.36 years old and ranging by 18–78 years old) who were recruited from Peking University Sixth Hospital. The diagnosis for schizophrenia was based on DSM-IV and conducted by two clinical doctors. All subjects were native Chinese speakers and right-handed. This study was carried out in conformity with the Declaration of Helsinki and was sustained by the ethics committee of Peking University Sixth Hospital. All subjects provided written consent after being fully informed of the procedures in the study.
The experiment was conducted in a quiet room and no entry was permitted during the experiment. The Chinese VFT (Quan et al., 2015) was initiated by a 30-s pre-task baseline period, followed by a 60-s task period and a 30-s post-task baseline period (Figure 1). There was a screen 1 m in front of the participants. During the pre-task and post-task baseline periods, the participants were asked to stare at the center of the screen and count from 1 to 5. During the 60-s task period, three Chinese characters (“,” “,” and “,” indicating white, sky, and big, respectively) were displayed on the screen and changed every 20 s. The participants were instructed to produce as many phrases or four-character idioms starting with these characters as they could.
Figure 1. Experimental flowchart. The experiment has three procedures: 30-s pre-task, 60-s verbal fluency task (VFT) task, and 30-s post-task.
fNIRS Measurement and Data Preprocessing
The measurement was conducted using a 52-channel near infrared spectrometer (ETG-4000, Hitachi Medical Co., Japan). The instrument had 33 probes (17 emitters and 16 detectors; Figure 2). The positioning of the receivers emitters was referred to an international 10-20 system of Electroencephalography (Oostenveld and Praamstra, 2001). To specify, the detector between Channel #5 and #6 was located at Fz, the emitters close to #43 and #52 were fitted around T4 and T3, and #46 and #49 were placed in Fp2 and Fp1, respectively. The measurement area covered the bilateral prefrontal and temporal cortices (Figure 2). The separation between the channels was 3 cm. In the experiments, each subject was measured with 120 s (30 s pre-task baseline, 60 s VFT and 30 s post-task baseline) at a sampling rate of 10 Hz. Hence, there were 1,200 signal points for each channel per subject. The measured Oxy-Hb signal was organized as a matrix with 300 × 1,200 × 52 (number of subjects × signal points × amount of channels).
Figure 2. The setting of functional near infrared spectroscopy (fNIRS) probe and channel. Hexagons stand for near-infrared light emitters, diamonds stand for near-infrared light receptors, and cycles stand for fNIRS channels. T3, T4, Fp1, Fp2, and Fz are the electrode positions in the international 10-20 system. rSFC, right superior frontal cortex; rSTC, right superior temporal cortex; rDLPFC, right dorsolateral prefrontal cortex; rVLPFC, right ventrolateral prefrontal cortex; mPFC, medial prefrontal cortex; lDLPFC, left dorsolateral prefrontal cortex; lVLPFC, left ventrolateral prefrontal cortex; lSFC, left superior frontal cortex; lSTC, left superior temporal cortex. The placement is in line with the configuration of Schecklmann et al. (2010).
The raw Oxy-Hb data were preprocessed through a band-pass filter of 0.009–0.08 Hz to remove the motion artifacts. The least square method was used to eliminate and remove the linear trend from the Oxy-Hb signals. MATLAB toolkit HomER2 (Huppert et al., 2009) was used to preprocess the original data.
Feature Extraction for Classification
The conventional classification methods usually utilized the time-domain features (for example, mean amplitude of Oxy-Hb during VFT). Frontal functional dysconnectivity is a salient feature of schizophrenia but it has not yet been applied for identifying schizophrenia. FCS was selected to characterize the effect and the following steps were used to obtain the value:
(1) Pearson’s correlation among the data from 52 channels was calculated by:
where , are the mean, and Sx, Sy are the standard deviations of the measured data xi and yi, respectively; n is the number of the data.
(2) FCS was calculated by:
As a consequence, 52 FCSs were derived per subject. We assessed the results from three kinds of approaches to identify schizophrenia.
- Single-feature: FCS from single channel
- Ensemble 52 FCSs: FCSs from 52 channels with or without dimensional reduction techniques
- Combine-features: FCS from the combination of several single-features
Four widely used classifiers in schizophrenic identification, LDA, GPC, KNN, and SVM, were trained and applied in the study. Since the details of these classifiers were extensively discussed previously (Mourão-Miranda et al., 2005; Yoon et al., 2007; Azechi et al., 2010; Tanaka and Katura, 2011; Dai et al., 2012; Arbabshirani et al., 2013; Hahn et al., 2013), we will not repeat them again. The MATLAB toolkits, GPML (v3.41; Rasmussen and Nickisch, 2010), LIBSVM (v3.1.22, and LDA (V1.0.0.03) were used in the analysis. KNN was realized by MATLAB function KNNCLASSIFY.
Note that the major parameters of these toolkits and functions used default or empirical values, with the exceptions of:
- KNN: k = (100)1/2 = 10; Euclidean distance is adopted to calculate the distances between the unlabeled sample and the labeled training samples. Traditionally, the Euclidean distance is appropriate when the issue included mutually correlated observations. As such, this distance needs to consider every variable and does not remove redundancies. The situation is very similar to our situation: the data from 52 channels are highly correlated and none can be simply removed. k is to set The method has been proposed by Mitra et al. (2002);
- SVM: RBF kennel; C and gamma were optimized by automated grid search and evaluated via 10-fold cross-validation. The optimization was conducted per case and the best RBF factors were provided along with the results.
Evaluation of the Classification Performance
To evaluate the performance of the individual classifier, both leave-one-out cross-validation (LOOCV) and 10-fold/ 20-fold CV were used to estimate the performance of the classifier. The subjects involved in the experiments were schizophrenia patients (positive, P) and the healthy controls (negative, N). The true positive (TP) and the true negative (TN) are the number of patients and healthy people being correctly classified, respectively. The false positive (FP) is the number of healthy people being classified as patients. The false negative (FN) is the number of patients being classified as healthy people. The performance of the classification method was assessed in terms of accuracy, sensitivity, and specificity as shown in:
Figure 3 shows the waveform of Oxy-Hb from 52 channels. The results were averaged across the healthy control and schizophrenic group. The reduced Oxy-Hb during VFT was obvious in patients with schizophrenia, being consistent with previous literature. The derived FCs were mapped in Figure 4. It was revealed that functional connections with high intensity were observed in the healthy controls, indicating dysconnectivity of schizophrenia.
Figure 3. The concentration curve of Oxy-Hb from 52 channels averaged from the healthy controls (blue lines) and the schizophrenic group (red lines). SZ, schizophrenic patients; HC, healthy controls.
Figure 4. Functional connectivity (FC) matrices averaged over schizophrenia group and healthy controls. SZ, schizophrenic patients (A); HC, healthy controls (B).
The overall accuracy from the top five channels to identify schizophrenia is shown in Figure 5 and is summarized in Table 1. LDA, KNN, and SVM demonstrated similar overall levels of accuracy, i.e., LDA: 72.50–81.00% (LOOCV), KNN: 78.00–82.00% (LOOCV), and SVM: 77.50–83.50% (LOOCV). GPC had the lowest accuracy at 67.00–69.50% (LOOCV). In terms of spatial distribution of the channels, although LDA, KNN, and SVM demonstrated laterality (left or right sidedness), in general, the best channels identified by these three classifiers were on the ventral part of the frontal cortices. In contrast, GPC utilized the FCS from the dorsal channels.
Figure 5. Overall accuracy achieved by single-channel feature using four classifiers [LDA (A), KNN (B), GPC (C), SVM (D)]. The top five channels were indicated by blue (the optimized C and gamma of the top five channels of SVM are shown in the Supplementary Figures 1–5). LDA, Linear Discriminant Analysis; KNN, k-Nearest Neighbor; GPC, Gaussian Processes classifier.
Table 1. Top five channels with functional connectivity strength (FCS) representing the best overall accuracy from the four classifiers (the optimized C and gamma of five channels of Support Vector Machine (SVM) are shown in Supplementary Figures 1–5).
Figure 6 shows the 52-feature results from different classifiers. Again, GPC achieved the lowest results (accuracy at 51.00% (LOOCV), with sensitivity at 55.00% (LOOCV) and specificity at 47.0% (LOOCV)). The other three classifiers had similar performances whilst SVM slightly outperformed the other two. In summary, SVM achieved the best accuracy at 86.50% (LOOCV), sensitivity at 91.00% (LOOCV), and specificity at 82% (LOOCV). LDA had the best accuracy at 83.00% (LOOCV), sensitivity at 85.00% (LOOCV), and specificity at 81.00% (LOOCV). KNN yielded the best accuracy at 77.00% (LOOCV), sensitivity at 84.00% (LOOCV), and specificity at 70.00% (LOOCV). The performance of the classifiers initially increased with the numbers of channels but stabilized when more channels were taken into consideration. It may indicate the existence of redundancy in this feature space.
Figure 6. Fifty-two functional connectivity strength (FCS) results [accuracy (A), sensitivity (B), specificity (C); the optimized C and gamma of best accuracy as 86.50% of SVM are shown in Supplementary Figure 6].
Dimensional Reduced 52-FCS Results
Principal component analysis (PCA) can convert multiple observations of potentially correlated variables into a set of linearly independent components. It implemented singular value decomposition to reduce the dimensionality of a dataset that consisted of a large number of interrelated variables, while retaining as much variation present in the dataset as possible (Abdi and Williams, 2010). For comparison, we also used two other variants of PCA in the analysis: Kernel PCA and Sparse PCA. Kernel PCA uses various kernel functions to project datasets into a higher dimensional feature space, where it is linearly separable. We selected Gaussian kernel in this case. It was realized by MATLAB function KernelPca.m (Kitayama, 2020). Sparse PCA is implemented on the basis of the inverse power method for nonlinear eigenproblems, which is introduced in detail by Hein and Bühler (2010). Moreover, the deflation scheme proposed by Bühler (2015) is adopted to compute multiple principal components. It was realized by the free software sparsePCA developed by Matthias Hein and Thomas Bühler [Copyright 2010–2020 Thomas Bühler and Matthias Hein (firstname.lastname@example.org). Machine Learning Group, Saarland University, Germany4].
We derived the first 21 principal components representing cumulative rates exceeding 93.4%. Classification based on the selected principal components is shown in Table 2. The results of 10-fold and 20-fold cross validation were provided for comparison.
Table 2. Classification performance of leave-one-out cross-validation (LOOCV), 10-fold, and 20-fold by using three dimensionality reduction methods (the optimized C and gamma of SVM are shown in Supplementary Figures 7–15).
Combined FCS Results
Further effort was made to assess the capability of schizophrenic identification using a certain combination of the channels. Since SVM yielded the best overall accuracy, the experiments were conducted only using this classifier. FCSs from 2, 3, 4, and 5 channels were selected from the five channels presenting the top capability on schizophrenic identification (presented in Figure 5). The results are shown in Table 3. Classification using FCS from three channels can achieve accuracy at 85.00% (LOOCV), sensitivity at 87.00% (LOOCV), and specificity at 83.00% (LOOCV).
Table 3. Overall accuracy of different combinations of FCS using SVM (the optimized C and gamma are shown in Supplementary Figures 16–41).
Schizophrenia has been considered a disorder of connectivity between various brain units (Elvevåg and Goldberg, 2000). The connections were found to be reduced by schizophrenia, as shown in Figure 4. This finding was consistent with studies using other imaging modalities (Bellani et al., 2010; Deng et al., 2017).
FCS measures the connectivity across different brain units, so as to identify the hubs playing important roles in information processing and communication during cognitive tasks (van den Heuvel and Sporns, 2013; Mears and Pollard, 2016). As shown in Table 1, the capability of discriminating schizophrenia was evident for the FCS at VLPFC and mPFC. mPFC relates to decision making and short- and long-term memory (Euston et al., 2012), and coordinates VLPFC and DLPFC functions (Peng et al., 2018). The neurophysiological functions of this cortex are associated with the symptoms of schizophrenia. The left VLPFC associates with the production of articulate language and in nonlinguistic tasks (Hickok and Poeppel, 2004, 2007), while the right VLPFC plays a role in linking working memory with episodic memory and in a series of complicated social behaviors (He et al., 2020). The reduced FCS of VLPFC in patients with schizophrenia may relate to the impairment of both verbal skills and social functions, which are the major symptoms of schizophrenia. In contrast, some channels were at STC, which mediates spatial awareness and exploration (Karnath, 2001). To summarize, these changes during VFT corresponded to the perturbed performance of schizophrenia patients (difficulty or incapability to produce four-character idioms).
PCA reduced the dimension of the feature space and saved the computational cost, while achieving satisfactory accuracy. The disadvantage of PCA was that the principle components could not be attributed to the data from the specific channel, thus concealing the regional neurophysiological changes. Using the FCS from three channels, the achieved performance was comparable to the current results: accuracy at 70–86%, sensitivity at 70–84%, and specificity at 65–93% (Arbabshirani et al., 2013; Chuang et al., 2014; Li et al., 2015; Pina-Camacho et al., 2015; Song et al., 2017). The method was not calculated from the time-domain values on single or multiple channels. It means that reliable results could be provided when integrated with the time-domain approaches.
There are some limitations to the present study. First, the individual schizophrenic episode was not identified and taken into analysis. However, it may have implications on the effected sites (Zhu et al., 2010). Second, the patients receiving medications and physical treatment were not ruled out from the study. Although previous studies have revealed a negligible medication effect on fNIRS signals, investigations of drug-free patients or of those receiving physical treatment (e.g., transcranial magnetic stimulation, electroconvulsive therapy, and neurofeedback) will be needed to allow further clinical applications of fNIRS (Fujita et al., 2011; Mihara et al., 2012; Monden et al., 2012). Third, educational background may have an impact on language ability. In our study, we did not categorize the subjects into more educational background groups because the number of subjects in each group would be sparse. But the two groups matched their educational backgrounds (the schizophrenia group included 14 graduate degrees, 20 undergraduate degrees, 20 college degrees, 31 senior high school degrees, and 15 junior high school degrees, and the healthy group included 15 graduate degrees, 20 undergraduate degrees, 20 college degrees, 30 high school degrees, and 15 junior high school degrees). Lastly, only three machine-learning classifiers, LDA, KNN, and SVM, were used in the study because they were the most popular machine-learning classifiers in discriminating patients with schizophrenia. The comparison of their performance was a topic being widely discussed while the individual performance seemed to be signal- and feature-dependent (Mourão-Miranda et al., 2005; Hahn et al., 2013; Li et al., 2015). In the present study, we conducted a similar comparison. Other classifiers, such as artificial neural networking, has not yet been applied but theoretically could be utilized in the identification of schizophrenia (Zheng et al., 2019a,b, 2020). We will try it in our future study.
The study proposed an FCS-based method to identify patients with schizophrenia. 52-channel Oxy-Hb data of fronto-temporal fNIRS were obtained during VFT from healthy and schizophrenic subjects. The FCS of each channel was calculated as features for classification. We investigated the performance of different classifiers, from FCS of all the 52 channels or from several channels. The method was in sharp contrast to most previous studies using the time-average data obtained from multiple channels. The classification results were comparable to the existing results. In addition, the method can detect the changes of hubs during VFT, which was in consistency with the symptoms of schizophrenia.
Data Availability Statement
The datasets presented in this article are not readily available because dataset are strictly restricted to non-commercial uses. Requests to access the datasets should be directed to email@example.com.
The studies involving human participants were reviewed and approved by the ethics committee of Peking University Sixth Hospital. The patients/participants provided their written informed consent to participate in this study.
TW conceptualized the experiments and supervised their implementation. JY and XJ contributed code to the project. WQ contributed to data acquisition. YL and BW validated the results and visualized them. XJ wrote the first version of the manuscript with input from JY and TN. TW and JY wrote the final version of the manuscript.
This study was supported by grants from the National Natural Science Foundation Project (No. 61971445), the National Key Research and Development Program of China (Nos. 2018YFC1314200 and 2019YFF0216302), and the National Science and Technology Major Project (No. 2018ZX10301201).
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.
- ^ www.gaussianprocess.org/gpml
- ^ www.csie.ntu.edu.tw/ cjlin/libsvm/
- ^ https://ww2.mathworks.cn/matlabcentral/fileexchange/29673-lda-linear-discriminant-analysis
- ^ www.ml.uni-saarland.de
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fninf.2020.00040/full#supplementary-material.
Arbabshirani, M. R., Kiehl, K. A., Pearlson, G. D., and Calhoun, V. D. (2013). Classification of schizophrenia patients based on resting-state functional network connectivity. Front. Neurosci. 7:133. doi: 10.3389/fnins.2013.00133
Azechi, M., Iwase, M., Ikezawa, K., Takahashi, H., Canuet, L., Kurimoto, R., et al. (2010). Discriminant analysis in schizophrenia and healthy subjects using prefrontal activation during frontal lobe tasks: a near-infrared spectroscopy. Schizophr. Res. 117, 52–60. doi: 10.1016/j.schres.2009.10.003
Bühler, T. (2015). A Flexible Framework for Solving Constrained Ratio Problems in Machine Learning. Germany: Saarland University. Ph.D. Thesis. Available online at: http://scidok.sulb.uni-saarland.de/volltexte/2015/6159/.
Chuang, C. C., Nakagome, K., Pu, S., Lan, T. H., Lee, C. Y., and Sun, C. W. (2014). Discriminant analysis of functional optical topography for schizophrenia diagnosis. J. Biomed. Opt. 19:011006. doi: 10.1117/1.jbo.19.1.011006
Demirci, O., Clark, V. P., and Calhoun, V. D. (2008). A projection pursuit algorithm to classify individuals using fMRI data: application to schizophrenia. NeuroImage 39, 1774–1782. doi: 10.1016/j.neuroimage.2007.10.012
Deng, S. P., Hu, W., Calhoun, V. D., and Wang, Y.-P. (2017). Integrating imaging genomic data in the quest for biomarkers for schizophrenia disease. IEEE/ACM Trans. Comput. Biol. Bioinform. 15, 1480–1491. doi: 10.1109/TCBB.2017.2748944
Ehlis, A. C., Herrmann, M. J., Plichta, M. M., and Fallgatter, A. J. (2007). Cortical activation during two verbal fluency tasks in schizophrenic patients and healthy controls as assessed by multi-channel near-infrared spectroscopy. Psychiatry Res. 156, 1–13. doi: 10.1016/j.pscychresns.2006.11.007
Fujita, Y., Takebayashi, M., Hisaoka, K., Tsuchioka, M., Morinobu, S., and Yamawaki, S. (2011). Asymmetric alternation of the hemodynamic response at the prefrontal cortex in patients with schizophrenia during electroconvulsive therapy: a near-infrared spectroscopy study. Brain Res. 1410, 132–140. doi: 10.1016/j.brainres.2011.06.052
Gordon, E., Barry, R. J., Anderson, J., Fawdry, R., Yong, C., Grunewald, S., et al. (1994). Single photon emission computed tomography (SPECT) measures of brain function in schizophrenia. Aust. N. Z. J. Psychiatry 28, 446–452. doi: 10.3109/00048679409075872
Hahn, T., Marquand, A. F., Plichta, M. M., Ehlis, A. C., Schecklmann, M. W., Dresler, T., et al. (2013). A novel approach to probabilistic biomarker-based classification using functional near-infrared spectroscopy. Hum. Brain Mapp. 34, 1102–1114. doi: 10.1002/hbm.21497
Hein, M., and Bühler, T. (2010). “An inverse power method for nonlinear eigenproblems with applications in 1-spectral clustering and sparse PCA,” in Advances in Neural Information Processing Systems 23 (NIPS 2010), Neural Information Processing Systems Conference, 847–855.
Huppert, T. J., Diamond, S. G., Franceschini, M. A., and Boas, D. A. (2009). HomER: a review of time-series analysis methods for near-infrared spectroscopy of the brain. Appl. Opt. 48, D280–D298. doi: 10.1364/ao.48.00d280
Ji, X. Y., Quan, W. X., Yang, L., Chen, J., Wang, J. J., and Wu, T. N. (2020). Classification of schizophrenia by seed-based functional connectivity using prefronto-temporal functional near infrared spectroscopy. J. Neurosci. Meth. 344. doi: 10.1016/j.jneumeth.2020.108874
Kanahara, N., Sekine, Y., Haraguchi, T., Uchida, Y., Hashimoto, K., Shimizu, E., et al. (2013). Orbitofrontal cortex abnormality and deficit schizophrenia. Schizophr. Res. 143, 246–252. doi: 10.1016/j.schres.2012.11.015
Kitayama, M. (2020). MATLAB-Kernel-PCA. GitHub. Available online at: https://github.com/kitayama1234/MATLAB-Kernel-PCA. Accessed July 9, 2020.
Kubota, Y., Toichi, M., Shimizu, M., Mason, R. A., Coconcea, C. M., Findling, R. L., et al. (2005). Prefrontal activation during verbal fluency tests in schizophrenia—a near-infrared spectroscopy (NIRS) study. Schizophr. Res. 77, 65–73. doi: 10.1016/j.schres.2005.01.007
Li, Z., Wang, Y., Quan, W., Wu, T., and Lv, B. (2015). Evaluation of different classification methods for the diagnosis of schizophrenia based on functional near-infrared spectroscopy. J. Neurosci. Methods 241, 101–110. doi: 10.1016/j.jneumeth.2014.12.020
Maguire, E. A., Gadian, D. G., Johnsrude, I. S., Good, C. D., Ashburner, J., Frackowiak, R. S. J., et al. (2000). Navigation-related structural change in the hippocampi of taxi drivers. Proc. Natl. Acad. Sci. U S A 97, 4398–4403. doi: 10.1073/pnas.070039597
Mears, D., and Pollard, H. B. (2016). Network science and the human brain: using graph theory to understand the brain and one of its hubs, the amygdala, in health and disease. J. Neurosci. Res. 94, 590–605. doi: 10.1002/jnr.23705
Mihara, M., Miyai, I., Hattori, N., Hatakenaka, M., Yagura, H., Kawano, T., et al. (2012). Neurofeedback using real-time near-infrared spectroscopy enhances motor imagery related cortical activation. PLoS One 7:e32234. doi: 10.1371/journal.pone.0032234
Monden, Y., Dan, H., Nagashima, M., Dan, I., Kyutoku, Y., Okamoto, M., et al. (2012). Clinically-oriented monitoring of acute effects of methylphenidate on cerebral hemodynamics in ADHD children using fNIRS. Clin. Neurophysiol. 123, 1147–1157. doi: 10.1016/j.clinph.2011.10.006
Mourão-Miranda, J., Bokde, A. L., Born, C., Hampel, H., and Stetter, M. (2005). Classifying brain states and determining the discriminating activation patterns: support vector machine on functional MRI data. NeuroImage 28, 980–995. doi: 10.1016/j.neuroimage.2005.06.070
Ohtani, T., Bouix, S., Hosokawa, T., Saito, Y., Eckbo, R., Ballinger, T., et al. (2014). Abnormalities in white matter connections between orbitofrontal cortex and anterior cingulate cortex and their associations with negative symptoms in schizophrenia: a DTI study. Schizophr. Res. 157, 190–197. doi: 10.1016/j.schres.2014.05.016
Peng, K., Steele, S. C., Becerra, L., and Borsook, D. (2018). Brodmann area 10: collating, integrating and high level processing of nociception and pain. Prog. Neurobiol. 161, 1–22. doi: 10.1016/j.pneurobio.2017.11.004
Pina-Camacho, L., Prieto, J. G., Parellada, M., Castro-Fornieles, J., Gonzalez-Pinto, A. M., Bombin, I., et al. (2015). Predictors of schizophrenia spectrum disorders in early-onset first episodes of psychosis: a support vector machine model. J. Eur. Child Adolesc. Psychiatry 24, 427–440. doi: 10.1007/s00787-014-0593-0
Quan, W., Wu, T., Li, Z., Wang, Y., Dong, W., and Lv, B. (2015). Reduced prefrontal activation during a verbal fluency task in Chinese-speaking patients with schizophrenia as measured by near-infrared spectroscopy. Prog. Neuropsychopharmacol. Biol. Psychiatry 58, 51–58. doi: 10.1016/j.pnpbp.2014.12.005
Rosenbaum, D., Haipt, A., Fuhr, K., Haeussinger, F. B., Metzger, F. G., Nuerk, H.-C., et al. (2017). Aberrant functional connectivity in depression as an index of state and trait rumination. Sci. Rep. 7:2174. doi: 10.1038/s41598-017-02277-z
Song, H., Bogdan, I. I. M., Wang, S., Dong, W., Quan, W., Dang, W., et al. (2016). “Automatic schizophrenia discrimination on fNIRS by using PCA and SVM,” in Proceedings of the 2016 IEEE International Conference on Bioinformatics and Biomedicine (BIBM), (Shenzhen, China: IEEE), 389–394.
Song, H., Chen, L., Gao, R. Q., Mihaita Bogdan, I. I., Yang, J., Wang, S., et al. (2017). Automatic schizophrenic discrimination on fNIRS by using complex brain network analysis and SVM. BMC Med. Inform. Decis. Mak. 17:166. doi: 10.1186/s12911-017-0559-5
Sugimura, Y., Watanabe, K., Ogawa, S., Kodama, T., Takeshita, M., Noda, T., et al. (2014). A discriminant model of mental disorders based on wave form patterns obtained by multi-channel near-infrared spectroscopy. Rinsho Byori 62, 147–152.
Suto, T., Fukuda, M., Ito, M., Uehara, T., and Mikuni, M. (2004). Multichannel near- infrared spectroscopy in depression and schizophrenia: cognitive brain activation study. Biol. Psychiatry 55, 501–511. doi: 10.1016/j.biopsych.2003.09.008
Takizawa, R., Kasai, K., Kawakubo, Y., Marumo, K., Kawasaki, S., Yamasue, H., et al. (2008). Reduced frontopolar activation during verbal fluency task in schizophrenia: a multi-channel near-infrared spectroscopy study. Schizophr. Res. 99, 250–262. doi: 10.1016/j.schres.2007.10.025
Tian, J., Wang, J., Quan, W., Lv, B., Liu, J., Hobbie, C., et al. (2019). The functional near-infrared spectroscopy in the diagnosis of schizophrenia. Eur. J. Psychiatry 33, 97–103. doi: 10.1016/j.ejpsy.2019.05.001
Tréhout, M., Leroux, E., Delcroix, N., and Dollfus, S. (2017). Relationships between corpus callosum and language lateralization in patients with schizophrenia and bipolar disorders. Bipolar Disord. 19, 496–504. doi: 10.1111/bdi.12526
Weiss, E. M., Hofer, A., Golaszewski, S., Siedentopf, C., Brinkhoff, C., Kremser, C., et al. (2004). Brain activation patterns during a verbal fluency test—a functional MRI study in healthy volunteers and patients with schizophrenia. Schizophr. Res. 70, 287–291. doi: 10.1016/j.schres.2004.01.010
Whitfield-Gabrieli, S., Thermenos, H. W., Milanovic, S., Tsuang, M. T., Faraone, S. V., McCarley, R. W., et al. (2009). Hyperactivity and hyperconnectivity of the default network in schizophrenia and in first-degree relatives of persons with schizophrenia. Proc. Natl. Acad. Sci. U S A 106, 1279–1284. doi: 10.1073/pnas.0809141106
Yang, H., Liu, J., Sui, J., Pearlson, G., and Calhoun, V. D. (2010). A hybrid machine learning method for fusing fMRI and genetic data: combining both improves classification of schizophrenia. Front. Hum. Neurosci. 4:192. doi: 10.3389/fnhum.2010.00192
Yoon, U., Lee, J.-M., Im, K., Shin, Y.-W., Cho, B. H., Kim, I. Y., et al. (2007). Pattern classification using principal components of cortical thickness and its discriminative pattern in schizophrenia. NeuroImage 34, 1405–1415. doi: 10.1016/j.neuroimage.2006.11.021
Zheng, Y. H., Song, H. H., Zhang, K. H., Zhang, K., Fan, J., and Liu, X. (2020). Dynamically spatio-temporal regularized correlation tracking. IEEE Trans. Neural Netw. Learn. Syst. 31, 2336–2347. doi: 10.1109/TNNLS.2019.2929407
Zheng, Y. H., Sun, L., Wang, S. F., Zhang, J., and Ning, J. (2019a). Spatially regularized structural support vector machine for robust visual tracking. IEEE Trans. Neural Netw. Learn. Syst. 30, 3024–3034. doi: 10.1109/TNNLS.2018.2855686
Zheng, Y. H., Wang, X. L., Zhang, G. Q., Xiao, B. H., Xiao, F., and Zhang, J. W. (2019b). Multi-kernel coupled projections for domain adaptive dictionary learning. IEEE Trans. Mult. 21, 2292–2304. doi: 10.1109/tmm.2019.2900166
Zhu, Y., Liu, X., Wang, H., Jiang, T., Fang, Y., Hu, H., et al. (2010). Reduced prefrontal activation during tower of London in first-episode schizophrenia: a multi-channel near-infrared spectroscopy study. Neurosci. Lett. 478, 136–140. doi: 10.1016/j.neulet.2010.05.003
Keywords: functional near infrared spectroscopy (fNIRS), schizophrenia, functional connectivity strength (FCS), support machine vector, classification
Citation: Yang J, Ji X, Quan W, Liu Y, Wei B and Wu T (2020) Classification of Schizophrenia by Functional Connectivity Strength Using Functional Near Infrared Spectroscopy. Front. Neuroinform. 14:40. doi: 10.3389/fninf.2020.00040
Received: 02 June 2020; Accepted: 22 July 2020;
Published: 07 October 2020.
Edited by:Heye Zhang, Sun Yat-sen University, China
Copyright © 2020 Yang, Ji, Quan, Liu, Wei and Wu. 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: Tongning Wu, firstname.lastname@example.org
† These authors have contributed equally to this work