ORIGINAL RESEARCH article

Front. Neuroinform., 07 October 2020

Volume 14 - 2020 | https://doi.org/10.3389/fninf.2020.00040

Classification of Schizophrenia by Functional Connectivity Strength Using Functional Near Infrared Spectroscopy

  • 1. China Academy of Information and Communications Technology, Beijing, China

  • 2. Institute of Electrical Engineering, Chinese Academy of Sciences, Beijing, China

  • 3. Peking 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

  • 4. School of Computer Science and Technology, Donghua University, Shanghai, China

  • 5. School of Computer Science and Technology, Xidian University, Xian, China

Abstract

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.

Introduction

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

Subject

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.

VFT Experiment

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

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 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:

Results

FC Matrices

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

Figure 4

Classification Results

Single-FCS Results

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

Table 1

LDAChannelBrain regionAccuracy
#41lSTC81.00%
#50lVLPFC78.50%
#35rVLPFC75.50%
#48mPFC74.00%
#36mPFC72.50%
KNN
#51lSTC82.00%
#35rVLPFC81.00%
#38mPFC79.50%
#50lVLPFC79.00%
#39lVLPFC78.00%
GPC
#5mPFC69.50%
#24rDLPFC68.50%
#39lVLPFC66.50%
#25rDLPFC67.00%
#19lDLPFC67.00%
SVM
#50lVLPFC83.50%
#41lSTC80.50%
#40lVLPFC79.00%
#44rSTC79.00%
#52lSTC77.50%

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).

52-FCS Results

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

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 (hein@cs.uni-saarland.de). 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

LOOCV10-fold20-fold
LDAAccuracy (PCA KernelPCA SparsePCA)80.50% 80.00% 82.00%79.00% 80.50% 80.50%81.50% 82.00% 82.50%
Sensitivity (PCA KernelPCA SparsePCA)82.00% 81.00% 83.00%81.00% 82.00% 83.00%83.00% 84.00% 83.00%
Specificity (PCA KernelPCA SparsePCA)79.00% 79.00% 81.00%77.00% 79.00% 78.00%80.00% 80.00% 82.00%
KNNAccuracy (PCA KernelPCA SparsePCA)79.00% 82.00% 79.50%80.50% 80.50% 82.00%80.00% 81.00% 80.50%
Sensitivity (PCA KernelPCA SparsePCA)82.00% 83.00% 81.00%81.00% 82.00% 84.00%82.00% 83.00% 81.00%
Specificity (PCA KernelPCA SparsePCA)76.00% 81.00% 78.00%80.00% 79.00% 80.00%78.00% 79.00% 80.00%
GPCAccuracy (PCA KernelPCA SparsePCA)72.00% 73.50% 75.00%68.50% 70.00% 70.50%70.00% 71.50% 70.50%
Sensitivity (PCA KernelPCA SparsePCA)73.00% 75.00% 77.00%70.00% 72.00% 72.00%72.00% 73.00% 72.00%
Specificity (PCA KernelPCA SparsePCA)71.00% 72.00% 73.00%67.00% 68.00% 70.00%68.00% 70.00% 69.00%
SVMAccuracy (PCA KernelPCA SparsePCA)87.50% 89.00% 88.00%85.50% 87.00% 88.50%87.50% 86.00% 89.50%
Sensitivity (PCA KernelPCA SparsePCA)89.00% 90.00% 89.00%86.00% 88.00% 89.00%86.00% 87.00% 90.00%
Specificity (PCA KernelPCA SparsePCA)86.00% 88.00% 87.00%85.00% 86.00% 88.00%89.00% 85.00% 89.00%

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

Selection of the channelsAccuracySensitivitySpecificity
#40#41#44#50#52
xx82.00%84.00%80%
xx81.50%82.00%81%
xx80.50%82.00%79%
xx83.00%85.00%81%
xx78.50%79.00%78%
xx82.00%83.00%79%
xx83.00%83.00%83%
xx79.00%80.00%78%
xx80.50%81.00%80%
xx83.50%85.00%82%
xxx81.50%83.00%80%
xxx85.00%87.00%83%
xxx84.00%86.00%82%
xxx82.00%82.00%80%
xxx84.50%85.00%84%
xxx83.00%85.00%81%
xxx80.00%82.00%78%
xxx79.50%81.00%78%
xxx80.00%81.00%79%
xxx81.00%84.00%78%
xxxx83.50%86.00%81%
xxxx83.00%84.00%82%
xxxx84.50%85.00%84%
xxxx82.00%85.00%79%
xxxx83.00%84.00%82%
xxxxx84.50%86.00%83%

Overall accuracy of different combinations of FCS using SVM (the optimized C and gamma are shown in Supplementary Figures 16–41).

Discussions

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.

Conclusion

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.

Statements

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 wutongning@caict.ac.cn.

Ethics statement

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.

Author contributions

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.

Funding

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fninf.2020.00040/full#supplementary-material.

References

  • 1

    AbdiH.WilliamsL. J. (2010). Principal component analysis. Wiley Interdiscip. Rev. Comput. Stat.2, 433450. 10.1002/wics.101

  • 2

    American Psychiatric Association. (1994). Diagnostic and Statistical Manual of Mental Disorders.Washington DC: American Psychiatric Association.

  • 3

    ArbabshiraniM. R.KiehlK. A.PearlsonG. D.CalhounV. D. (2013). Classification of schizophrenia patients based on resting-state functional network connectivity. Front. Neurosci.7:133. 10.3389/fnins.2013.00133

  • 4

    AzechiM.IwaseM.IkezawaK.TakahashiH.CanuetL.KurimotoR.et al. (2010). Discriminant analysis in schizophrenia and healthy subjects using prefrontal activation during frontal lobe tasks: a near-infrared spectroscopy. Schizophr. Res.117, 5260. 10.1016/j.schres.2009.10.003

  • 5

    BellaniM.CerrutiS.BrambillaP. (2010). Orbitofrontal cortex abnormalities in schizophrenia. Epidemiol. Psichiatr. Soc.19, 2325. 10.1017/s1121189x00001561

  • 6

    BuckleyP. F.MillerB. J.LehrerD. S.CastleD. J. (2009). Psychiatric comorbidities and schizophrenia. Schizophr. Bull.35, 383402. 10.1093/schbul/sbn135

  • 7

    BühlerT. (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/.

  • 8

    BullmoreE.SpornsO. (2009). Complex brain networks: graph theoretical analysis of structural and functional systems. Nat. Rev. Neurosci.10, 186198. 10.1038/nrn2575

  • 9

    ChuangC. C.NakagomeK.PuS.LanT. H.LeeC. Y.SunC. W. (2014). Discriminant analysis of functional optical topography for schizophrenia diagnosis. J. Biomed. Opt.19:011006. 10.1117/1.jbo.19.1.011006

  • 10

    DaiD.WangJ.HuaJ.HeH. (2012). Classification of ADHD children through multimodal magnetic resonance imaging. Front. Syst. Neurosci.6:63. 10.3389/fnsys.2012.00063

  • 11

    DemirciO.ClarkV. P.CalhounV. D. (2008). A projection pursuit algorithm to classify individuals using fMRI data: application to schizophrenia. NeuroImage39, 17741782. 10.1016/j.neuroimage.2007.10.012

  • 12

    DengS. P.HuW.CalhounV. D.WangY.-P. (2017). Integrating imaging genomic data in the quest for biomarkers for schizophrenia disease. IEEE/ACM Trans. Comput. Biol. Bioinform.15, 14801491. 10.1109/TCBB.2017.2748944

  • 13

    EhlisA. C.HerrmannM. J.PlichtaM. M.FallgatterA. 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, 113. 10.1016/j.pscychresns.2006.11.007

  • 14

    ElvevågB.GoldbergT. E. (2000). Cognitive impairment in schizophrenia is the core of the disorder. Crit. Rev. Neurobiol.14, 121. 10.1615/critrevneurobiol.v14.i1.10

  • 15

    EustonD.GruberA.McNaughtonB. (2012). The role of medial prefrontal cortex in memory and decision making. Neuron76, 10571070. 10.1016/j.neuron.2012.12.002

  • 16

    FristonK. J.FrithC. D. (1995). Schizophrenia: a disconnection syndrome. J. Clin. Neurosci.3, 8997.

  • 17

    FujitaY.TakebayashiM.HisaokaK.TsuchiokaM.MorinobuS.YamawakiS. (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, 132140. 10.1016/j.brainres.2011.06.052

  • 18

    GordonE.BarryR. J.AndersonJ.FawdryR.YongC.GrunewaldS.et al. (1994). Single photon emission computed tomography (SPECT) measures of brain function in schizophrenia. Aust. N. Z. J. Psychiatry28, 446452. 10.3109/00048679409075872

  • 19

    GreiciusM. (2008). Resting-state functional connectivity in neuro-psychiatric disorders. Curr. Opin. Neurol.21, 424430. 10.1097/WCO.0b013e328306f2c5

  • 20

    HahnT.MarquandA. F.PlichtaM. M.EhlisA. C.SchecklmannM. W.DreslerT.et al. (2013). A novel approach to probabilistic biomarker-based classification using functional near-infrared spectroscopy. Hum. Brain Mapp.34, 11021114. 10.1002/hbm.21497

  • 21

    HeZ.ZhaoJ.ShenJ.MuhlertN.ElliottR.ZhangD. (2020). The right VLPFC and downregulation of social pain: a TMS study. Hum. Brain Mapp.41, 13621371. 10.1002/hbm.24881

  • 22

    HeinM.BühlerT. (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, 847855.

  • 23

    HickokG.PoeppelD. (2004). Dorsal and ventral streams: a framework for understanding aspects of the functional anatomy of language. Cognition92, 6799. 10.1016/j.cognition.2003.10.011

  • 24

    HickokG.PoeppelD. (2007). The cortical organization of speech processing. Nat. Rev. Neurosci.8, 393402. 10.1038/nrn2113

  • 25

    HoshiY. (2003). Functional near-infrared optical imaging: utility and limitations in human brain mapping. Psychophysiology40, 511520. 10.1111/1469-8986.00053

  • 26

    HuppertT. J.DiamondS. G.FranceschiniM. A.BoasD. A. (2009). HomER: a review of time-series analysis methods for near-infrared spectroscopy of the brain. Appl. Opt.48, D280D298. 10.1364/ao.48.00d280

  • 27

    JiX. Y.QuanW. X.YangL.ChenJ.WangJ. J.WuT. N. (2020). Classification of schizophrenia by seed-based functional connectivity using prefronto-temporal functional near infrared spectroscopy. J. Neurosci. Meth.344. 10.1016/j.jneumeth.2020.108874

  • 28

    KanaharaN.SekineY.HaraguchiT.UchidaY.HashimotoK.ShimizuE.et al. (2013). Orbitofrontal cortex abnormality and deficit schizophrenia. Schizophr. Res.143, 246252. 10.1016/j.schres.2012.11.015

  • 29

    KarnathH. O. (2001). New insights into the functions of the superior temporal cortex. Nat. Rev. Neurosci.2, 568576. 10.1038/35086057

  • 30

    KitayamaM. (2020). MATLAB-Kernel-PCA.GitHub. Available online at: https://github.com/kitayama1234/MATLAB-Kernel-PCA. Accessed July 9, 2020.

  • 31

    KubotaY.ToichiM.ShimizuM.MasonR. A.CoconceaC. M.FindlingR. L.et al. (2005). Prefrontal activation during verbal fluency tests in schizophrenia—a near-infrared spectroscopy (NIRS) study. Schizophr. Res.77, 6573. 10.1016/j.schres.2005.01.007

  • 32

    LiZ.WangY.QuanW.WuT.LvB. (2015). Evaluation of different classification methods for the diagnosis of schizophrenia based on functional near-infrared spectroscopy. J. Neurosci. Methods241, 101110. 10.1016/j.jneumeth.2014.12.020

  • 33

    MaguireE. A.GadianD. G.JohnsrudeI. S.GoodC. D.AshburnerJ.FrackowiakR. S. J.et al. (2000). Navigation-related structural change in the hippocampi of taxi drivers. Proc. Natl. Acad. Sci. U S A97, 43984403. 10.1073/pnas.070039597

  • 34

    MearsD.PollardH. 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, 590605. 10.1002/jnr.23705

  • 35

    MiharaM.MiyaiI.HattoriN.HatakenakaM.YaguraH.KawanoT.et al. (2012). Neurofeedback using real-time near-infrared spectroscopy enhances motor imagery related cortical activation. PLoS One7:e32234. 10.1371/journal.pone.0032234

  • 36

    MitraP.MurthyC. A.PalS. K. (2002). Unsupervised feature selection using feature similarity. IEEE Trans. Pattern Anal. Mach. Intell.24, 301312. 10.1109/34.990133

  • 37

    MondenY.DanH.NagashimaM.DanI.KyutokuY.OkamotoM.et al. (2012). Clinically-oriented monitoring of acute effects of methylphenidate on cerebral hemodynamics in ADHD children using fNIRS. Clin. Neurophysiol.123, 11471157. 10.1016/j.clinph.2011.10.006

  • 38

    Mourão-MirandaJ.BokdeA. L.BornC.HampelH.StetterM. (2005). Classifying brain states and determining the discriminating activation patterns: support vector machine on functional MRI data. NeuroImage28, 980995. 10.1016/j.neuroimage.2005.06.070

  • 39

    OhtaniT.BouixS.HosokawaT.SaitoY.EckboR.BallingerT.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, 190197. 10.1016/j.schres.2014.05.016

  • 40

    OostenveldR.PraamstraP. (2001). The five percent electrode system for high-resolution EEG and ERP measurements. Clin. Neurophysiol.112, 713719. 10.1016/s1388-2457(00)00527-7

  • 41

    PengK.SteeleS. C.BecerraL.BorsookD. (2018). Brodmann area 10: collating, integrating and high level processing of nociception and pain. Prog. Neurobiol.161, 122. 10.1016/j.pneurobio.2017.11.004

  • 42

    Pina-CamachoL.PrietoJ. G.ParelladaM.Castro-FornielesJ.Gonzalez-PintoA. M.BombinI.et al. (2015). Predictors of schizophrenia spectrum disorders in early-onset first episodes of psychosis: a support vector machine model. J. Eur. Child Adolesc. Psychiatry24, 427440. 10.1007/s00787-014-0593-0

  • 43

    QuanW.WuT.LiZ.WangY.DongW.LvB. (2015). Reduced prefrontal activation during a verbal fluency task in Chinese-speaking patients with schizophrenia as measured by near-infrared spectroscopy. Prog. Neuropsychopharmacol. Biol. Psychiatry58, 5158. 10.1016/j.pnpbp.2014.12.005

  • 44

    RasmussenC. E.NickischH. (2010). Gaussian processes for machine learning (GPML) toolbox. J. Mach. Learn. Res.11, 30113015. 10.1115/1.4002474

  • 45

    RosenbaumD.HaiptA.FuhrK.HaeussingerF. B.MetzgerF. G.NuerkH.-C.et al. (2017). Aberrant functional connectivity in depression as an index of state and trait rumination. Sci. Rep.7:2174. 10.1038/s41598-017-02277-z

  • 46

    SchecklmannM.RomanosM.BretscherF. (2010). Prefrontal oxygenation during working memory in ADHD. J. Psychiatr. Res.44, 621628. 10.1016/j.jpsychires.2009.11.018

  • 47

    SongH.BogdanI. I. M.WangS.DongW.QuanW.DangW.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), 389394.

  • 48

    SongH.ChenL.GaoR. Q.Mihaita BogdanI. I.YangJ.WangS.et al. (2017). Automatic schizophrenic discrimination on fNIRS by using complex brain network analysis and SVM. BMC Med. Inform. Decis. Mak.17:166. 10.1186/s12911-017-0559-5

  • 49

    SugimuraY.WatanabeK.OgawaS.KodamaT.TakeshitaM.NodaT.et al. (2014). A discriminant model of mental disorders based on wave form patterns obtained by multi-channel near-infrared spectroscopy. Rinsho Byori62, 147152.

  • 50

    SutoT.FukudaM.ItoM.UeharaT.MikuniM. (2004). Multichannel near- infrared spectroscopy in depression and schizophrenia: cognitive brain activation study. Biol. Psychiatry55, 501511. 10.1016/j.biopsych.2003.09.008

  • 51

    TakizawaR.KasaiK.KawakuboY.MarumoK.KawasakiS.YamasueH.et al. (2008). Reduced frontopolar activation during verbal fluency task in schizophrenia: a multi-channel near-infrared spectroscopy study. Schizophr. Res.99, 250262. 10.1016/j.schres.2007.10.025

  • 52

    TanakaH.KaturaT. (2011). Classification of change detection and change blindness from near-infrared spectroscopy signals. J. Biomed. Opt.16:087001. 10.1117/1.3606494

  • 53

    TianJ.WangJ.QuanW.LvB.LiuJ.HobbieC.et al. (2019). The functional near-infrared spectroscopy in the diagnosis of schizophrenia. Eur. J. Psychiatry33, 97103. 10.1016/j.ejpsy.2019.05.001

  • 54

    TréhoutM.LerouxE.DelcroixN.DollfusS. (2017). Relationships between corpus callosum and language lateralization in patients with schizophrenia and bipolar disorders. Bipolar Disord.19, 496504. 10.1111/bdi.12526

  • 55

    van den HeuvelM. P.SpornsO. (2013). An anatomical substrate for integration among functional networks in human cortex. J. Neurosci.33, 1448914500. 10.1523/JNEUROSCI.2128-13.2013

  • 56

    WeissE. M.HoferA.GolaszewskiS.SiedentopfC.BrinkhoffC.KremserC.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, 287291. 10.1016/j.schres.2004.01.010

  • 57

    Whitfield-GabrieliS.ThermenosH. W.MilanovicS.TsuangM. T.FaraoneS. V.McCarleyR. 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 A106, 12791284. 10.1073/pnas.0809141106

  • 58

    YangH.LiuJ.SuiJ.PearlsonG.CalhounV. D. (2010). A hybrid machine learning method for fusing fMRI and genetic data: combining both improves classification of schizophrenia. Front. Hum. Neurosci.4:192. 10.3389/fnhum.2010.00192

  • 59

    YoonU.LeeJ.-M.ImK.ShinY.-W.ChoB. H.KimI. Y.et al. (2007). Pattern classification using principal components of cortical thickness and its discriminative pattern in schizophrenia. NeuroImage34, 14051415. 10.1016/j.neuroimage.2006.11.021

  • 60

    ZhengY. H.SongH. H.ZhangK. H.ZhangK.FanJ.LiuX. (2020). Dynamically spatio-temporal regularized correlation tracking. IEEE Trans. Neural Netw. Learn. Syst.31, 23362347. 10.1109/TNNLS.2019.2929407

  • 61

    ZhengY. H.SunL.WangS. F.ZhangJ.NingJ. (2019a). Spatially regularized structural support vector machine for robust visual tracking. IEEE Trans. Neural Netw. Learn. Syst.30, 30243034. 10.1109/TNNLS.2018.2855686

  • 62

    ZhengY. H.WangX. L.ZhangG. Q.XiaoB. H.XiaoF.ZhangJ. W. (2019b). Multi-kernel coupled projections for domain adaptive dictionary learning. IEEE Trans. Mult.21, 22922304. 10.1109/tmm.2019.2900166

  • 63

    ZhuY.LiuX.WangH.JiangT.FangY.HuH.et al. (2010). Reduced prefrontal activation during tower of London in first-episode schizophrenia: a multi-channel near-infrared spectroscopy study. Neurosci. Lett.478, 136140. 10.1016/j.neulet.2010.05.003

Summary

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

Volume

14 - 2020

Edited by

Heye Zhang, Sun Yat-sen University, China

Reviewed by

Zicheng Liu, Arctic University of Norway, Norway; Xi Cheng, Télécom ParisTech, France

Updates

Copyright

*Correspondence: Tongning Wu

These authors have contributed equally to this work

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics