Keywords
fMRI, Schizophrenia, Time-series, Classification
This article is included in the Brainhack Global collection.
fMRI, Schizophrenia, Time-series, Classification
This version of the manuscript is solely updated in accordance with the suggestions given by the respected referees. This version mainly elaborates the discussion section. It describes each of the identified brain regions showing changes in functional activation and compares with the existing studies to validate the finding of the paper. Here, some demographic details of the subjects are added to the dataset table. This version also states the scope of future work of the study.
See the author's detailed response to the review by Sagarika Bhattacharjee
Schizophrenia is a severe mental disorder that affects different regions of the brain, often involving hallucinations and delusions. Functional magnetic resonance imaging (fMRI) data comprising 3D brain scans acquired over time (thus resulting in a 4D set) is often used to study brain regions affected by schizophrenia. Each voxel of the 3D brain volume is associated with a time series of signal intensity values. General linear model (GLM)1 and independent component analysis (ICA)2 are often employed to study the voxel activity by transforming the 4D time-series data to a 3D spatial map.
The present work involves a novel application of mean deviation on time-series fMRI data to identify the distinct voxels that show high functional activation during a task. The work aims to identify the relevant brain regions that are affected in schizophrenia. Further, the identified voxels (features) are used to distinguish between schizophrenia patients and healthy subjects.
The time-series fMRI data having 1.5T strength was taken from the FBIRN phase – II data repository3 available at site 0009 and site 0010. From the dataset, four different runs of auditory oddball task data of 34 schizophrenic patients (group G1) and 34 healthy controls (group G2) were extracted. Every run of each subject’s data contains 140 brain volumes acquired in 280 seconds time (TR = 2 seconds). Table 1 shows the dataset details.
Pre-processing of the fMRI data was done using SPM8 toolbox in Matlab2014b. The temporal variation was corrected using slice timing correction, followed by the motion correction using realignment. Each of the fMRI scans was spatially normalized into standard Montreal Neurological Institute (MNI) space using an EPI template yielding voxel dimension of 3×3×3 mm3. Finally smoothing was done using a 9×9×9 mm3 full width at half maximum (FWHM) Gaussian kernel, resulting in a 3D brain volume containing 53×63×46, i.e., 1,53,594 voxels.
The activation pattern of the voxels was analysed in two phases.
Phase I. In the first phase, identification of voxels exhibiting high activation pattern (anytime during its time-course) is carried out for each subject. As the study focused on the variation in the signal intensity of the voxels (V) over time, absolute mean deviation () for each of the 140 time points was computed for each voxel, and the median (M) of the 140 values of was found. Mean deviation () values were compared with α times M (α was chosen to be 3, based on experimentation) to identify whether a voxel exhibited high level of activation at any time during the 140 units of time. This voxel-wise analysis was performed for all the voxels of a given subject. Thus, a set of relevant voxels showing high degree of activation was obtained for each subject.
Phase II. In the second phase, a common subset of voxels exhibiting high degree of activation across all the subjects within a group was obtained. Finally, both the subsets belonging to groups G1 (schizophrenia patients) and G2 (healthy controls) were merged to get the set S. The voxels in set S were backtracked to MNI brain space and finally mapped into Talairach’s space4 to identify the brain regions. This procedure has been described in Algorithm 1.
Classification. The set S was used to distinguish between schizophrenia patients and healthy subjects using two classifiers, viz., support vector machine (SVM) with sigmoid kernel5 and extreme learning machine (ELM) classifier6.
Experimental settings. All the implementations were done in MATLAB2014b. Parameter α was varied in the range of 1 to 7 in steps of 1 to identify the number of voxels that exhibited a high level of activations during the task. When the value of α was taken as 1 and 2, a large number of voxels showed activation level higher than α times M, resulting in set S having voxels that represents almost the entire brain. However, for α = 3, it was found that set S contained only 1580 distinct voxels that mapped to the brain regions which are generally affected in schizophrenia. When α was taken as more than 3, the number of voxels in the set S were close to zero rendering it too small for any meaningful analysis. Thus, α = 3 was found to be the most suitable value.
Further, the set S of voxels obtained α = 3 was used to fine-tune the classifiers. The SVM classifier gave the best results for the regularization parameter C = 1.09, and sigmoid kernel based ELM classifier gave best the results with 503 hidden neurons.
To evaluate the distinguishing capability of the voxels/features in set S, a comparison was done between the classification accuracy obtained using S and the accuracy obtained using the voxels set given by the GLM based approach. In this case, GLM was applied using SPM8 toolbox to convert the 4D time-series fMRI data to 3D contrast map for each subject. The GLM yielded an activation map comprising around 60000 voxels out of 153594 which were activated during the task.
Notations:
m (=34): the number of subjects in each group
n (=140): the number of observations in a run
Vi : time-series of ith voxel
i.e. Vi = [vi,1 vi,2 vi,3 ⋯ vi,n];
Steps:
1. Calculate absolute mean deviation for each voxel using
2. Find median Mi of
3. For each subject k ∈{1,2,...,m}, select the set Vsk of voxels that show deviation higher than αMi.
4. Find the group wise intersection of the voxels selected in step 3 for groups G1 and G2
5. Merge the two sets, obtained in step 4 to obtain set S
6. Map S into the brain space to identify affected regions.
A comparison of the results of the classification accuracies obtained using feature sets given by the GLM and the proposed approach is shown in Table 2. The features selected by the proposed approach when backtracked to Talairach’s space revealed the brain regions that are generally affected in schizophrenia7–9, which validates the efficacy of the approach. The distribution of the selected voxels that distinguish the schizophrenia patients from the healthy subjects is shown in Figure 1 (a–d). The results show the increased changes in functional activation in the regions such as occipital lobe, frontal lobe, posterior lobe, and temporal lobe. When looking into the level of gyri, certain changes in activation pattern are seen in superior temporal gyrus, lingual gyrus, cuneus, declive, medial frontal gyrus, and middle occipital gyrus. Some regions in Brodmann areas (BA 18, 10, 9, 17, 19, 32, 6, 37, 21, 22, 46, and 47) also show distinct changes in functional activation in schizophrenics when compared to healthy controls. Figure 2 (a–c) show the activated voxels when plotted on a sample fMRI image for an axial, coronal and sagittal view of the brain.
GLM | Proposed approach | |
---|---|---|
Number of voxels | ~ 60,000 | ~ 1580 |
SVM with Sigmoid kernel | 32.45% | 76.47% |
ELM with Sigmoid kernel | 57.35% | 61.46% |
Identified brain regions at different levels of hierarchy, namely, hemisphere level (a), lobes level (b), gyrus level (c), and Brodmann’s area level (d).
Unlike other conventional methods such as GLM to select the voxels showing a statistically significant response to the experimental conditions10, the proposed approach identifies the neural activity in a particular voxel with the help of bold signal over time, irrespective of any experimental condition. The proposed approach does not require any details for the task and conditions. It works on the temporal values of each voxel for each subject's data one by one. Like other multi-voxel pattern analysis (MVPA) methods10–12, this approach also tries to find the participation of multiple voxels when selecting the final set of relevant voxels across a particular group of the subjects.
The classification accuracies, as shown in Table 2, demonstrate the efficacy of the proposed methodology. The reduced set of 1580 voxels achieved a much higher accuracy when compared to the GLM approach. The approach gives a better result for both of the SVM and ELM classifiers when compared to the GLM approach.
Figures 1 (a–d) show the distribution of the selected voxels for each level of brain regions. These regions show distinct changes in functional activation in schizophrenia patients when compared to healthy controls, and thereby distinguish between schizophrenia and healthy subjects with high classification accuracy. Most of the regions identified in the study comply with the existing literature13–16.
From the Figure 1 (b), change in functional activation can be seen in the frontal lobe which is responsible for motor function, executive functions and attention17,18. Literature19,20 suggest frontal lobe functional dysfunction in schizophrenia. Significant changes are observed in the temporal lobe and occipital lobe as seen in Figure 1 (b). The temporal lobe is basically responsible for holding primary auditory perception such as hearing, and occipital lobe is responsible for visual perception. As schizophrenia patients suffer from auditory and visual hallucinations, the functional deficit in these regions are responsible for the symptoms, studies21,22 also suggest functional changes in these areas in schizophrenia.
As seen in Figure 1 (c), significant changes in functional activations are found in regions such as superior frontal gyrus, superior temporal gyrus, lingual gyrus, medial frontal gyrus, middle occipital gyrus, anterior cingulate, cuneus, and declive. Superior temporal gyrus contains the primary auditory cortex which is responsible for processing sound, sending sensory information to auditory cortex and also to specify the sound frequencies precisely. Previous studies23–25 also showed that the superior temporal gyrus gets affected in schizophrenia. Figure 1 (c) shows functional changes in superior frontal gyrus, which is mainly involved in self-awareness26. Literature13,27 suggests changes in superior frontal gyrus. The literature also suggests functional changes in middle occipital gyrus28. Since schizophrenia patients also suffer from visual hallucinations and deficiency in visual attention, the dysfunctioning of the areas such as declive and cuneus (BA 17) may play role in the disorder. Studies found that cuneus15,29 and declive30 show functional changes in schizophrenia. Lingual gyrus is basically linked to function for visual processing31. The result of this study also shows subtle functional changes in lingual gyrus indicating difficulties in visual abilities in schizophrenia32,33. Even functional abnormality in anterior cingulate was found in several studies16,34.
In the level of Brodmann’s area, as seen in Figure 1 (d), BA 18 and 19 show functional changes in schizophrenia patients when compared to healthy controls. These regions lie in the occipital cortex, mainly responsible for the interpretation of images35. Studies36,37 show that these regions are commonly affected in schizophrenia. BA 10, lies in the prefrontal cortex, is responsible for executive functions such as attention, processing of working memory and taking decision for future actions38. Similar to the previous studies39,40, this study also identifies the functional changes in BA 10. BA 9 lies in the frontal lobe, mainly responsible for short-term memory, auditory verbal attention. This region may play a vital role in auditory hallucination and in retrieving the short-term memory in schizophrenia patients41,42. BA 37, responsible for visual fixation43 and recognizing true-false memory, is also mentioned in the previous studies44,45. BA 21 and 22 lie in the temporal cortex, believed to play a role in auditory processing are found to be affected in schizophrenia19,20. Other than these regions, the result shows significant changes in functional activations in the areas such as BA 6, BA 37, BA 46, and BA 11, which are also reported in the literature14,46,47.
This paper identifies the affected brain regions in schizophrenia and compares them with the previous studies. As the study focused on the statistical measures derived from the voxel values across the time course, the effect of covariates such as level of education, duration of disease, and medication history could not be incorporated. However, during the group-wise analysis (mentioned in Phase II of the methodology), a grouping of subjects’ data was also performed on the basis of gender and age. The results obtained were quite uniform across all ages and genders. Although this study was performed on the auditory oddball task fMRI data, it would be interesting to explore the applicability of the approach on the resting-state fMRI data and other performance-based tasks.
This work describes a simple and fast feature selection algorithm based on mean deviation for time-series fMRI data to identify the activated brain voxels that are generally affected in schizophrenia. The proposed approach was found to be efficient in selecting a minimal set of relevant voxels directly from time-series 4D fMRI data. The obtained voxel set was capable of distinguishing between healthy and schizophrenic subjects. One may explore the possibility of applying this approach to fMRI data of other psychological disorders.
The Matlab source codes, a text file containing dataset details including subject ID and their age, and the instructions for the study can be found at: https://github.com/IndraChatterjee/AnomalyDetection_TimeSeries_fMRI_Schizophrenia.
The complete source codes are archived in a publicly accessible record at: https://doi.org/10.5281/zenodo.143853948
License: CC0
The four runs of auditory oddball task fMRI data from the FBIRN phase II repository can be downloaded from http://schizconnect.org/ querying 1.5T fMRI data for healthy and schizophrenia subjects available at site 0009 and 0010. The list of subjects chosen for this study is mentioned in the ‘DataDetails_FBIRN15T.txt’ file available at the GitHub repository. Users are required to sign-up to SchizConnect to download data and conditions of use are as written in the data use agreement of the FBIRN project.
Cameron Craddock confirms that the author has an appropriate level of expertise to conduct this research, and confirms that the submission is of an acceptable scientific standard. Cameron Craddock declares the following competing interests: I am the Chair of Brainhack, and this organisation awarded this paper this year's Brainhack poster prize. Affiliation: Associate Professor of Diagnostic Medicine, Dell Medical School, The University of Texas at Austin, Austin, TX, USA.
The author would like to thank the organizers and all the attendees of 2018 OHBM Brainhack Singapore.
Data used for this study were hosted in the Function BIRN Data Repository (http://fbirnbdr.birncommunity.org:8080/BDR/) using Project Accession Number 2007-BDR-6UHZ1, supported by grants to the Function BIRN (U24-RR021992) Testbed funded by the National Center for Research Resources at the National Institutes of Health, U.S.A.
Views | Downloads | |
---|---|---|
F1000Research | - | - |
PubMed Central
Data from PMC are received and updated monthly.
|
- | - |
Competing Interests: I know the author personally from a conference.
Is the work clearly and accurately presented and does it cite the current literature?
Yes
Is the study design appropriate and is the work technically sound?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Yes
If applicable, is the statistical analysis and its interpretation appropriate?
I cannot comment. A qualified statistician is required.
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Partly
Competing Interests: I know the author personally from a conference
Is the work clearly and accurately presented and does it cite the current literature?
Yes
Is the study design appropriate and is the work technically sound?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Yes
If applicable, is the statistical analysis and its interpretation appropriate?
Yes
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Yes
Competing Interests: No competing interests were disclosed.
Alongside their report, reviewers assign a status to the article:
Invited Reviewers | ||
---|---|---|
1 | 2 | |
Version 2 (revision) 20 Dec 18 |
read | |
Version 1 08 Oct 18 |
read | read |
Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list:
Sign up for content alerts and receive a weekly or monthly email with all newly published articles
Already registered? Sign in
The email address should be the one you originally registered with F1000.
You registered with F1000 via Google, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Google account password, please click here.
You registered with F1000 via Facebook, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Facebook account password, please click here.
If your email address is registered with us, we will email you instructions to reset your password.
If you think you should have received this email but it has not arrived, please check your spam filters and/or contact for further assistance.
Comments on this article Comments (0)