Original Research ARTICLE
More Consistently Altered Connectivity Patterns for Cerebellum and Medial Temporal Lobes than for Amygdala and Striatum in Schizophrenia
- 1Department of Psychiatry, Technische Universität München, München, Germany
- 2TUM-Neuroimaging Center, Klinikum Rechts der Isar, Technische Universität München, München, Germany
- 3Klinik für Psychiatrie und Psychotherapie, Klinikum der Universität München, München, Germany
- 4School of Computer Science and Engineering, University of Electronic Science and Technology of China, Chengdu, China
- 5Big Data Research Center, University of Electronic Science and Technology of China, Chengdu, China
- 6Center for Information in BioMedicine, University of Electronic Science and Technology of China, Chengdu, China
- 7Department of Neuroradiology, Technische Universität München, München, Germany
- 8Department of Nuclear Medicine, Technische Universität München, München, Germany
Background: Brain architecture can be divided into a cortico-thalamic system and modulatory “subcortical-cerebellar” systems containing key structures such as striatum, medial temporal lobes (MTLs), amygdala, and cerebellum. Subcortical-cerebellar systems are known to be altered in schizophrenia. In particular, intrinsic functional brain connectivity (iFC) between these systems has been consistently demonstrated in patients. While altered connectivity is known for each subcortical-cerebellar system separately, it is unknown whether subcortical-cerebellar systems’ connectivity patterns with the cortico-thalamic system are comparably altered across systems, i.e., if separate subcortical-cerebellar systems’ connectivity patterns are consistent across patients.
Methods: To investigate this question, 18 patients with schizophrenia (3 unmedicated, 15 medicated with atypical antipsychotics) and 18 healthy controls were assessed by resting-state functional magnetic resonance imaging (fMRI). Independent component analysis of fMRI data revealed cortical intrinsic brain networks (NWs) with time courses representing proxies for cortico-thalamic system activity. Subcortical-cerebellar systems’ activity was represented by fMRI-based time courses of selected regions-of-interest (ROIs; i.e., striatum, MTL, amygdala, cerebellum). Correlation analysis among ROI- and NWs-time courses yielded individual connectivity matrices [i.e., connectivity between NW and ROIs (allROIs-NW, separateROI-NW), only NWs (NWs-NWs), and only ROIs (allROIs-allROIs)] as main outcome measures, which were classified by support-vector-machine-based (SVM) leave-one-out cross-validation. Differences in classification accuracy were statistically evaluated for consistency across subjects and systems.
Results: Correlation matrices based on allROIs-NWs yielded 91% classification accuracy, which was significantly superior to allROIs-allROIs and NWs-NWs (56 and 74%, respectively). Considering separate subcortical-cerebellar systems, cerebellum-NWs and MTL-NWs reached highest accuracy values with 91 and 85%, respectively, while those of striatum-NW and amygdala-NW were significantly lower with about 65% classification accuracy.
Conclusion: Results provide initial evidence for differential consistency of altered intrinsic connectivity patterns between subcortical-cerebellar systems and the cortico-thalamic system. Data suggest that differential dysconnectivity patterns between subcortical-cerebellar and cortical systems might reflect different disease states or patient subgroups.
Widespread brain dysconnectivity is an essential element in the pathophysiology of schizophrenia (Stephan et al., 2009). The cortico-thalamic system is particularly affected by aberrant functional and structural connectivity (Flynn et al., 2003; Begré and Koenig, 2008; Lynall et al., 2010; van den Heuvel et al., 2010; Welsh et al., 2010; Woodward et al., 2012; Palaniyappan et al., 2013; Wagner et al., 2013, 2015; Manoliu et al., 2014). For example, the coherence of ongoing, slowly fluctuating neural activity (<0.1 Hz), which reflects a basic form of brain organization called intrinsic functional connectivity (iFC; Fox and Raichle, 2007), is consistently changed in patients relative to healthy controls for both cortico-thalamic (Welsh et al., 2010; Woodward et al., 2012) and cortico-cortical connectivity (Lynall et al., 2010; Palaniyappan et al., 2013; Manoliu et al., 2014; van den Heuvel and Fornito, 2014). One should note that, rather than overall reduced or increased connectivity, alterations reflect rather complex reorganization of network interaction.
Beyond the cortico-thalamic system, several other non-cortical/non-thalamic brain systems are involved in schizophrenia. For example, the striatum shows increased pre-synaptic dopamine concentrations that are crucial for schizophrenia’s psychotic symptoms (Howes et al., 2012) and increased activity is linked with symptom dimensions, for putamen with positive symptoms during psychosis, for ventral striatum with negative symptoms during remission (Sorg et al., 2013). Moreover, substantial alterations within the medial temporal lobes (MTLs) include altered microstructure, transmitters, perfusion, and activity, particularly in the hippocampus (Tamminga et al., 2010). In more detail, structural deficits were recently attributed to hippocampal subfields in large samples of schizophrenia patients (Mathew et al., 2014; Arnold et al., 2015; Haukvik et al., 2015) and were associated with decreased default-mode-network suppression in healthy subjects at genetic risk (Seidman et al., 2014). Patients’ memory deficits, particularly with regard to relational information (vs. item-specific), have been linked with diminished hippocampal activation (Ragland et al., 2015). Furthermore, abnormalities in amygdala and cerebellum, ranging from molecular to microcircuit alterations, are consistent features of schizophrenia and are linked with emotional and cognitive deficits, respectively (Andreasen and Pierson, 2008; Benes, 2010; Anticevic et al., 2012; Mier et al., 2014; Millan et al., 2014; Vai et al., 2015). In more detail, amygdala alterations are heterogeneous and appear to differ depending on the clinical syndrome and disorder state (Rasetti et al., 2009). For example, amygdala hyperactivity has been shown to relate to paranoid ideation in schizophrenia, suggesting symptom-based subgroups (Pinkham et al., 2015). Findings across subjects at high-risk, early and chronic states of schizophrenia, suggest differentially affected amygdala iFC from increased to reduced levels, also reflecting symptom severity (Anticevic et al., 2014; Mukherjee et al., 2014). Considering altered cerebellar functional connectivity, mainly reduced connectivity has been reported consistently across different tasks and disorder states (Lungu et al., 2013). Interestingly, the cerebellum has been identified as one of the key regions of hypoconnectivity with the thalamus related to conversion of subjects at ultra high risk to psychosis (Anticevic et al., 2015).
Intriguingly, all these brain systems—striatum, MTL, amygdala, and cerebellum—are massively connected with the cortico-thalamic system (Swanson, 2003). They are assumed to play key modulatory roles in cortico-thalamic processes (Buzsaki, 2006) and beyond structural aberrations (Bois et al., 2015; van Erp et al., 2015), several studies demonstrated their aberrant iFC with the cortico-thalamic system in patients with schizophrenia (Collin et al., 2011; Anticevic et al., 2012; Fornito et al., 2013; Kraguljac et al., 2014; Peters et al., in press). We use the term “subcortical-cerebellar system” to refer to these systems. We are aware that it is a matter of debate to separate these systems from the cortico-thalamic system; for example, the striatum is massively integrated in cortico-basal ganglia-thalamo-cortical loops. However, differences in anatomy and connectivity, physiology and function, as well as their distinct role in schizophrenia pathophysiology suggest that treating these systems separately from thalamo-cortical connectivity is plausible (Buzsaki, 2006; Andreasen and Pierson, 2008; Benes, 2010; Tamminga et al., 2010; Howes et al., 2012). While previous studies investigated each “subcortical-cerebellar” system and its connectivity separately, it is unknown how the whole pattern of connectivity changes between “subcortical-cerebellar” systems and the cortico-thalamic system is configured in schizophrenia, particularly whether abnormal connectivity is consistent across patients for different systems. This might be an important question, since graded patterns of aberrant connectivity across subcortical-cerebellar systems might vary across disease states and patient subgroups. To obtain initial evidence for this question, the current study investigates connectivity patterns in the cortico-thalamic system across mentioned subcortical-cerebellar systems, and—most importantly—between these two types of systems in psychotic patients with schizophrenia to evaluate the inter-individual consistency of potential changes across systems.
To reach the study’s objective, the following steps have been performed (for an overview, see Figure S1): (i) We assessed slowly fluctuating ongoing brain activity (Leopold and Maier, 2012) in patients with schizophrenia and healthy controls with resting-state functional MRI (rs-fMRI). (ii) Ongoing brain activity is organized in intrinsic brain networks (NWs) characterized by specific iFC patterns (Fox and Raichle, 2007). We used intrinsic networks, which cover preferentially the cerebral cortex, as proxy for cortico-thalamic sub-systems (Figure 1). Corresponding time courses, reflecting networks’ activity (NWs) were used for further analysis. (iii) Concerning subcortical-cerebellar systems, we focused on striatum, amygdala, MTL, and cerebellum (Figure 2). For each system, we chose typical sub-regions, which have been used previously to study systems’ iFC (details in the method part). For these regions-of-interest, we extracted representative rs-fMRI signal time courses (ROIs). (iv) Then, we calculated iFC connectivity matrices (i.e., Pearson correlation of pairs of signal time courses of NWs or ROIs, respectively) for combinations of allROIs-NWs, NWs-NWs, allROIs-allROIs as well as separateROI-NWs for each subject. (v) Finally, we used multivariate pattern classification (i.e., support vector machine) and leave-one-out cross validation to estimate group differences and particularly inter-individual consistency of iFC patterns for different subcortical-cerebellar systems.
Figure 1. Cortical networks. Displayed are spatial maps of cortical networks (NWs) derived from independent component analysis (ICA) of patients’ and healthy controls’ resting state fMRI data, representing the cortico-thalamic sub-systems. The 22 NWs shown were selected following an automated spatial multiple-regression on reference-templates provided by Allen et al. (2011).
Figure 2. Subcortical-cerebellar systems. Displayed are regions of interest (ROI) based on coordinates derived from previous studies in case of striatum, medial temporal lobe and cerebellum (Kahn et al., 2008; Etkin et al., 2009; Krienen and Buckner, 2009; Peters et al., in press). Definition of amygdala ROI was derived from Anatomy toolbox for SPM (http://www.fz-juelich.de/inb/inb-3/spm_anatomy_toolbox).
We hypothesized that: 1. Connectivity matrices based on allROIs-NWs reflect more consistent connectivity patterns within groups of patients and healthy controls than allROIs-allROIs and NWs-NWs connectivity matrices. Therefore, they classify subjects more accurately. 2. Separate-ROI-NWs connectivity matrices describe differentially consistent connectivity patterns. Therefore, separate ROIs classify subjects with different accuracy.
Materials and Methods
Eighteen healthy controls and Eighteen patients participated in the study (Table 1). All participants had been investigated in a previous study (Manoliu et al., 2014), which addressed cortical iFC changes for three selected intrinsic networks i.e., default mode, salience, and central executive network; here, we focus on a larger number of cortical networks and their iFC with subcortical-cerebellar systems. Written informed consent in accordance with the Human Research Committee guidelines of the Klinikum Rechts der Isar, Technische Universität München was obtained from all participants. Patients were recruited from the Department of Psychiatry, Klinikum Rechts der Isar, TU München, controls by word-of-mouth advertising. Participants’ examination included medical history, psychiatric interview, psychometric assessment, and additionally blood tests for patients. Psychiatric diagnoses relied on the Diagnostic and Statistical Manual of Mental Disorders (DSM-IV). To assess psychiatric diagnoses, the Structured Clinical Interview for DSM-IV (SCID-I) was used (Spitzer et al., 1992). The global level of social, occupational, and psychological functioning was measured with the Global Assessment of Functioning Scale (GAF; Spitzer et al., 1992). For rating severity of clinical symptoms on the day of scanning, the Positive and Negative Syndrome Scale (PANSS) was applied (Kay et al., 1987). Clinical-psychometric assessment was completed by psychiatrists DS and MS who have been professionally trained for SCID and PANSS-based interviews with inter-rater reliability of more than 95%. Inclusion criteria for the study were diagnosis of schizophrenia and age between 18 and 60 years. Exclusion criteria were current or past neurological or internal systemic disorder, current depressive or manic episode, substance abuse (except for nicotine), schizoaffective disorder, and cerebral pathology in MRI.
All patients included were diagnosed with paranoid schizophrenia during acute psychosis as indicated by clinical exacerbation and increased positive symptom scores on the PANSS (Table 1). Seven out of eighteen had significant hallucinations (PANSS P3 ≥ 3), 15 delusions (P1 ≥ 3). The mean duration of illness was 7.15 years (SD = 6.89 years), the mean number of hospital stays was 2.98 (SD = 2.48). Concerning medication, three patients were free of any antipsychotic medication, the others received mono- or dual therapy with atypical antipsychotic medication including Amisulpride (n = 2), Olanzapine (n = 11), Clozapine (n = 4), Quetiapine (n = 2), Ziprasidone (n = 1), Risperidone (n = 5), Aripiprazole (n = 2), Paliperidone (n = 3; cp. Table S2 for individual medication protocols and dosage; see Table S1 for mean chlorpromazine equivalent dose). All healthy subjects were free of psychotropic medication and any neurological and psychiatric disorder, current and in history.
All participants underwent 10 min of rs-fMRI with the instruction to keep their eyes closed and not to fall asleep. We verified that subjects stayed awake by interrogating via intercom immediately after the rs-fMRI scan. Before and after scanning, a medical examination of patients validated their stable condition and investigated whether they had feelings of odd situations during the scanning. No patient dropped out during the scanning session.
MRI was carried out on a 3 T whole body scanner (Achieva, Philips Healthcare). FMRI was based on gradient echo EPI sequence (TE = 35 ms, TR = 2000 ms, flip angle = 82°, FoV = 220 × 220 mm2, matrix = 80 × 80, 32 slices, slice thickness = 4 mm, and 0 mm interslice gap; 300 volumes). T1-weighted anatomical MRI was based on magnetization-prepared rapid acquisition gradient echo sequence (TE = 4 ms, TR = 9 ms, TI = 100 ms, flip angle = 5°, FoV = 240 × 240 mm2, matrix = 240 × 240, 170 slices, voxel size = 1 × 1 × 1 mm3).
Functional MRI Data Preprocessing and Analysis
For each participant, the first three rs-fMRI scans were discarded due to magnetization effects. SPM8 (Wellcome Department of Cognitive Neurology, London) was used for motion correction, spatial normalization into the stereotactic space of the Montreal Neurological Institute (MNI) and spatial smoothing with an 8 mm × 8 mm × 8 mm Gaussian kernel. To ensure data quality, particularly concerning motion-induced artifacts, temporal signal-to-noise ratio (tSNR) and point-to-point head motion were estimated for each subject (Murphy et al., 2007; Van Dijk et al., 2012). Point-to-point motion was defined as the absolute displacement of each brain volume compared to its previous volume. Moreover, root mean square (RMS) of the translational head movement parameters was calculated for each subject. Excessive head motion (cumulative motion translation >3 mm and mean point-to-point translation or rotation >0.15 mm or 0.1°) was applied as exclusion criterion. None of the participants had to be excluded. Two-sample t-tests yielded no significant differences between groups regarding mean point-to-point translation or rotation of any direction (p > 0.18), RMS (p > 0.25), or tSNR (p > 0.35).
Subcortical-Cerebellar Systems: ROIs and Preprocessing
For each subcortical-cerebellar system, representative ROIs were defined based on previous studies, which investigated systems’ cortical iFC (Kahn et al., 2008; Etkin et al., 2009; Krienen and Buckner, 2009; Peters et al., in press; Figure 2). For the amygdala, left and right basolateral amygdala ROI were derived from the Anatomy toolbox for SPM1 following Etkin et al. (2009) and converted to corresponding ROIs via Marsbar2. For other subcortical-cerebellar systems, center coordinates in striatum, cerebellum, and MTL, respectively, were derived from the literature and spherical ROIs (6 mm radius) were created via Marsbar: left and right ventral and dorsal striatum following Peters et al. (in press) with MNI coordinates [XYZ]: (±12, 9, −9) and (±24, 12, 0), left and right Crus 1 and 2 following Krienen and Buckner (2009) with MNI coordinates [XYZ]: (±12.0, −80.0, −24.0) and (±22.0, −86.0, −40.0), left and right hippocampus and parahippocampus ROIs following Kahn et al. (2008) with MNI coordinates [XYZ]: (±24.0, −18.0, −18.0) and (±26.0, −40.0, −12.0).
Using in-house scripts employing SPM 8 and matlab3 routines, for each subject and ROI, fMRI-time courses were extracted, Butterworth bandpass-filtered for the frequency range from 0.009 to 0.08 Hz, and reduced to ROI-representative time courses by singular value decomposition. To account for partial volume and movement effects, linear effects of global gray matter (GM), white matter (WM), cerebrospinal fluid (CSF) fMRI-signal, and six movement parameters were regressed out. We included the global GM signal as a nuisance regressor since it is thought to reflect a combination of physiological processes (such as cardiac and respiratory fluctuations) and scanner drift. To extract the nuisance covariate time series for GM, WM and CSF, each individual’s high-resolution T1-weighted structural image was segmented. Mean images of the study sample’s T1-segmentation were used to create ROIs for the extraction of GM, WM, and CSF nuisance signals.
Cortico-Thalamic Networks: NW Identification and Selection
NWs were identified by ICA of pre-processed fMRI data within an identical framework as defined by Allen et al. (2011) and as applied previously (Manoliu et al., 2014). Briefly, selection of the optimal ICA model-order to analyze rs-fMRI data is subject of ongoing debate. However, it has been demonstrated that a model-order around 70 components may represent an optimal level to detect between-group differences and to avoid false positive results. Therefore, we followed the approach of Allen et al. (2011), who used ICA model order of 75 components. Further reasons to follow Allen’s approach were both greater comparability of results across studies and most importantly, reduction of subjective bias for network selection by using Allen’s templates to identify networks-of-interest. Preprocessed data were decomposed into 75 spatially independent components within a group-ICA framework based on the infomax-algorithm as implemented in the GIFT-software4. FMRI data were concatenated and reduced by two-step principal component analysis, followed by independent component estimation with the infomax-algorithm. We subsequently ran 40 ICA (ICASSO) to ensure stability of the estimated components. This resulted in a set of average group components, which are then back-reconstructed into single subject space. Each back-reconstructed component consists of a spatial z-map reflecting the component’s functional connectivity pattern across space and an associated time course reflecting component’s activity across time.
To select ICs which covered mainly the cortico-thalamic system, we chose templates as provided by Allen et al. (2011), i.e., t-maps of 28 components that reflect canonical intrinsic networks. We chose components mainly covering the cerebral cortex (22 of 28 maps, see Supplementary Material) to reflect networks-of-interest in an automated and objective way, and performed multiple spatial regression analyses of our 75 independent components’ spatial maps on these templates. Components of highest correlation coefficient with the templates were selected, resulting in 22 ICs of interest. In the end, this approach yielded a component’s z-map and time course for each subject and cortical intrinsic network, which reflect network’s coherent activity. Time courses of NW were used as proxies for cortico-thalamic system activity for further analysis.
To generate individual feature vectors as input parameter for pattern classification via SVM, respective Pearson-correlations of NWs’ and ROIs’ time courses and subsequent Fisher’s r-to-z transformation were performed for each subject (allROIs-NWs) and represented as typical correlation matrices. Correlation matrices for NWs-NWs and allROIs-allROIs connectivity were analogously created. Besides a global analysis including all ROIs and NWs as input for classification, also separateROIs-NWs correlation matrices were analyzed, including only one of the four subcortical-cerebellar system-ROIs at a time (i.e., cerebellum, amygdala, striatum, or MTL), thus facilitating comparison of classification accuracy obtained by each of the subcortical-cerebellar systems’ component.
The resulting correlation matrices of the two groups (feature vectors) were used for pattern recognition based classification via the WEKA toolbox5. Specifically, a support-vector machine based algorithm (SVM) was applied together with leave-one-out cross validation. The basic idea of SVM procedures (Vapnik, 1999) is to construct a separating hyperplane between multi-dimensional training instances (i.e., correlation matrices) of both classes (i.e., groups). Among all possible hyperplanes, that one with the maximum margin between classes is selected and subsequently applied to classify test instances. In order to evaluate classification effectiveness, SVM procedure has been realized within a leave-one-out validation framework, in which in each round of classification one subject is used as test instance while the others represent training instances. The specific algorithm implemented in WEKA applies sequential minimal optimization and a linear kernel (Platt, 1999). Measures-of-interest for SVM-based classification were classification accuracy, sensitivity and specificity, which together reflect the power of the particular feature (i.e., NWs-NWs, allROIs-allROIs, allROIs-NWs, and separateROIs-NMs) to discriminate groups. These measures reflect how consistent patients/controls are correctly classified due to the corresponding iFC pattern. We intended to examine different subcortical-cerebellar systems for such consistency of connectivity pattern changes. In addition, statistical significance of differences in classification measures was determined by reading out the individual classification results (i.e., correct or misclassified) and pairwise identification of diverging classification accuracy for pairs of features between NWs-NWs, allROIs-allROIs, allROIs-NWs, and separateROIs-NWs (Japkowicz and Shah, 2011). In more detail, let (X1,Y1), (X2,Y2), …, (X35,Y35) indicate the individual classification results, where Xi = 0 or Yi = 0 (misclassified), Xi = 1 or Yi = 1 (correct). We counted the number of pairs Xi > Yi. Let W be the number of pairs (Xi > Yi), N is the number of pairs that Xi and Yi are different, then W follows a binomial distribution W ~ B(N, 0.5). The exact p-value was obtained by checking binomial distribution.
Typical Cortical Intrinsic Networks have been Identified in Patients and Controls
ICA analysis of rs-fMRI data from patients and controls identified 22 intrinsic networks, which cover mainly the cerebral cortex (Figure 1; for more details, Figures S2–6). Networks were consistent across groups, and matched with previously defined networks of an identical methodological approach in 603 healthy controls (Allen et al., 2011).
Superior Classification Accuracy for allROIs-NWs than allROIs-allROIs or NWs-NWs
When separating individual patients from healthy controls via SVM classification accuracy for iFC patterns between all subcortical-cerebellar systems’ ROIs and cortical NWs (allROIs-NWs) was 91% with high sensitivity and specificity levels of about 90% (Table 2). Classification accuracies for iFC patterns among subcortical-cerebellar systems allROIs-allROIs and cortical networks NWs-NWs, were significantly lower with 56 and 73%, respectively (Tables 2, 4).
Table 2. Classification results based on intrinsic functional connectivity (iFC) among subcortical-cerebellar systems and cortical networks, respectively.
Superior Classification Accuracy for MTL-NWs and Cerebellum-NWs than Striatum-NWs and Amygdala-NWs
When separating individual patients from healthy controls via SVM for different subcortical-cerebellar systems separately, classification accuracies were significantly higher for iFC patterns of cerebellum-NWs (91%, sensitivity 82%, specificity 100%) and MTL-NWs (85%, sensitivity 77%, specificity 94%) than for striatum-NWs (65%, sensitivity 65%, specificity 65%) and amygdala-NW (68%, sensitivity 71%, specificity 65%; Table 3 and S3, Figure 3).
Table 3. Classification results based on intrinsic functional connectivity (iFC) between single subcortical-cerebellar systems and cortical networks.
Table 4. Statistical evaluation of differences in classification accuracy for different connectivity matrices.
Figure 3. ROI-network connectivity. Connectivity matrices between cortical networks (NWs, cp. Figure 1) and different ROI, i.e., cerebellum, medial temporal lobe, amygdala and striatum (cp. Figure 2), served as input features for support-vector-machine-based classification of schizophrenia patients and healthy subjects. Displayed classification accuracies [%] indicate differential discriminatory power of separate ROI-NWs-systems and the degree of inter-subject consistency of ROI-specific connectivity pattern changes in schizophrenia.
The current study focused on iFC between different subcortical-cerebellar systems (ROIs) and cortical intrinsic networks (NWs) in patients with schizophrenia. (i) Taking all subcortical-cerebellar systems of interest together, i.e., MTL, cerebellum, striatum, and amygdala (allROIs-NWs), iFC separated patients from controls with an accuracy of about 90%, which was significantly superior in comparison to allROIs-allROIs and NWs-NWs connectivity. This result highlights the role of iFC between subcortical-cerebellar systems and cortical networks in contrast to both cortico-cortical and subcortical-cerebellar-subcortical-cerebellar systems’ connectivity in schizophrenia. (ii) Across subcortical-cerebellar systems, MTL-NWs and cerebellum-NWs connectivity showed highest classification accuracies of about 90%, which were significantly higher in comparison to iFC patterns of striatum-NWs and amygdala-NWs with accuracies of about 65%. This finding provides first evidence for differential across-system consistency of changes in intrinsic connectivity between different subcortical-cerebellar systems and the cortico-thalamic system. Data suggest that the whole pattern of cortical-subcortical-cerebellar system dysconnectivity might be differential for patients, maybe indicating different patient subgroups or disease states.
iFC Between Subcortical-Cerebellar Systems and Cortical Networks is more Consistently Changed than iFC within Modulatory Systems and Cortical Networks, Respectively
Changes in iFC patterns between all subcortical-cerebellar and cortical networks are highly consistent in psychotic patients with schizophrenia indicated by classification accuracy, sensitivity, and specificity of about 90% (Table 2). Interestingly, such iFC changes are more consistent than those among only subcortical-cerebellar (allROIs-allROIs) or only cortical networks (NWs-NWs), respectively, highlighting the prominent role of subcortical-cerebellar system connectivity into the cortico-thalamic system in schizophrenia (Tables 2, 5). This result is in line with several previous studies, which demonstrated separately impaired iFC between a selected subcortical-cerebellar system and cortical regions (Collin et al., 2011; Anticevic et al., 2012; Fornito et al., 2013; Kraguljac et al., 2014; Wagner et al., 2015; Peters et al., in press). For instance, in a study by Wagner et al. (2015) an abnormal effective connectivity was observed between thalamus, anterior cingulate and dorsolateral prefrontal cortex in patients with schizophrenia compared to healthy controls. Moreover, patients have been found to show a decreased functional connectivity between the putamen with right anterior insula and dorsal prefrontal cortex and the ventral striatum with left anterior insula compared to healthy control subjects (Peters et al., in press). The present result goes beyond these studies by demonstrating that the whole pattern of iFC between several subcortical-cerebellar systems and cortical systems is more consistently aberrant than iFC among subcortical-cerebellar systems or among cortical systems alone.
Table 5. Statistical evaluation of differences in classification accuracy for different connectivity matrices focused on single subcortical-cerebellar systems.
Across Subcortical-Cerebellar Systems, iFC Changes with Cortical Networks are more Consistent for MTLs and Cerebellum than Striatum and Amygdala
Intrinsic connectivity changes between cerebellum and MTL respectively, and cortical networks were highly consistent in patients with accuracies of about 90% (Table 3, Figure 3). This result is well in line with previous findings reporting aberrant seed-based iFC of hippocampus and cerebellum, respectively, with cortical areas (Collin et al., 2011; Kraguljac et al., 2014). In a longitudinal study over 2 years in more than 60 patients with schizophrenia, Duan et al. (2015) demonstrated lasting iFC decrease between hippocampus and both, temporal and parietal areas, suggesting aberrant hippocampus-cortex iFC as a permanent feature in schizophrenia. Collin et al. (2011) showed consistently decreased cerebellar iFC into widespread parts of the cortex, in patients and their unaffected siblings, suggesting impaired cerebellum iFC as putative risk endophenotype for schizophrenia and therefore, as stable feature in the disease. In line with the studies from Collin et al. (2011) and Duan et al. (2015), which suggest longer lasting MTL- and cerebellum-iFC changes, we found highly consistent changes of both MTL- and cerebellum-cortex iFC in patients. This finding suggests that aberrant cortex-cerebellum iFC and cortex-MTL iFC together are stable features in schizophrenia.
In contrast, iFC between striatum and amygdala, respectively, and cortical networks separated patients from controls with significantly lower accuracy levels of about 65% (Table 3 and S3, Figure 3). Aberrant amygdala and striatum iFC with the cortex is in line with several previous studies. Anticevic et al. (2012) showed decreased iFC between amygdala and the orbitofrontal cortex, which was only present in patients of early and chronic disease stages but not in unaffected high-risk persons. This finding suggests that iFC between cortex and amygdala is modulated by disease state. Concerning striatum, previous studies demonstrated a complex pattern of in- and decreased iFC with the cortex along the ventromedial to dorsolateral axis of the striatum; in particular, this pattern of change seems to be modulated by the disease state and to be intimately linked with psychotic symptoms (Fornito et al., 2013; Manoliu et al., 2013; Peters et al., in press). However, psychotic symptoms seem to be not necessarily linked with pathophysiological striatum changes. For example, Demjaha et al. (2012) demonstrated in psychotic non-responders of anti-dopaminergic treatment, that presynaptic dopamine synthesis levels are unchanged in comparison to healthy controls, suggesting striatal dopamine synthesis being independent from psychosis. In line with the studies from Anticevic and Demjaha, which suggest inconsistent amygdala and striatum changes in patients with schizophrenia, we found that changes of both amygdala- and striatum-cortex iFC were of moderate accuracy and sensitivity, and of significant lower accuracy than iFC changes of cerebellum and MTL. Considering that cerebellum and MTL iFC changes appear as rather permanent features in schizophrenia, this finding suggests that aberrant cortex-amygdala iFC and cortex-striatum iFC are less consistent in schizophrenia and may reflect different disease states or heterogeneity across patients.
Differential Patterns of Aberrant iFC Between Subcortical-Cerebellar Systems and Cortical Networks as Candidates to Separate Sub-Groups of Patients and/or Disease States
Widespread brain dysconnectivity is central to schizophrenia (Stephan et al., 2009). Several models proposed specific patterns of dysconnectivity being associated with specific syndromes, symptoms, states, and subgroups of patients in the disease (Aleman and Kahn, 2005; Andreasen and Pierson, 2008; Stephan et al., 2009; Tamminga et al., 2010; Williamson and Allman, 2012). Due to both their essential role in modulating the cortico-thalamic system (Buzsaki, 2006) and the overwhelming evidence of their multi-level changes in schizophrenia (Andreasen and Pierson, 2008; Benes, 2010; Tamminga et al., 2010; Howes et al., 2012), striatum, cerebellum, MTL, and amygdala are promising candidates to specify schizophrenia’s dysconnectivity. The current study provides initial evidence that these main modulatory systems of the cortico-thalamic system and their iFC with cortical networks might specify distinct dysconnectivity patterns across patients despite limitations such as medication effects in particular as laid out further below. Future studies are necessary to identify how patterns of subcortical-cerebellar systems’ dysconnectivity vary across disease states, unaffected at-risk persons, and patients, and whether specific patterns may define patient subgroups.
Limitations and Methodological Issues
(i) Most patients of our study were treated with antipsychotic drugs (Table S1). Antipsychotic medication has been demonstrated to alter functional connectivity in patients with schizophrenia (Sambataro et al., 2010). Therefore, we cannot exclude that observed changes in iFC in comparison to healthy controls are confounded by medication effects. However, one should note that the main focus of the study was differential iFC changes across subcortical-cerebellar systems in patients i.e., differences in changes of patients instead of differences in patients relative to controls, suggesting that—due to consistent confound across patients—findings might be more robust against medication influences. Nevertheless, studies in non-medicated patients would be favorable. (ii) Sample size of the study is small, limiting the power of our study results to some degree. Therefore and due to confounding medication effects, we have to categorize our study as a preliminary study about differential iFC changes across subcortical-cerebellar systems in schizophrenia. (iii) One might ask why we used SVM instead of canonical two-sample t-tests to compare iFC across subjects. Typical two-sample t-tests for seed-based iFC maps would provide information about spatially specific group differences in cortical iFC for distinct subcortical-cerebellar systems; but such group difference maps are difficult to compare across subjects due to both, varying spatial extent of cortical changes and most importantly, the lack of direct measures about how consistent iFC changes are across systems and subjects. Combining measures of classification and their statistical comparison via binominal testing yields exactly such information. (iv) Why did we select intrinsic networks as proxies for cortico-thalamic sub-systems instead of brain-atlas-based anatomically defined regions? Concerning brain activity, we focused on slowly fluctuating activity as measured by rs-fMRI and correspondingly on iFC of such ongoing activity. Overwhelming evidence indicates that iFC is organized by intrinsic networks particularly in the human cortex and exceeds anatomically defined cortex parcellations (Fox and Raichle, 2007).
Results provide preliminary evidence for more consistent iFC changes with cortical networks for MTL and cerebellum than for striatum and amygdala in schizophrenia. Differential iFC changes might reflect distinct disease states or patient subgroups.
CS and HP designed the study; JB, MS, DS and HF recruited participants and acquired data; VR, CZ and AW acquired data; JS and HP analyzed data; CS, KK, VR, AW and HP interpreted data; CS, KK and HP drafted the article; all authors critically revised and approved the final version of the article.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by the German Federal Ministry of Education and Research (BMBF 01EV0710 to AW) and the Kommission für Klinische Forschung of the Klinikum Rechts der Isar der Technischen Universität München (KKF 8765162 to CS). This work was partially supported by the National Natural Science Foundation of China (Grant Nos. 61403062 and 61433014 to JS) and Fundamental Research Funds for the Central Universities (Grant No. ZYGX2014J053 to JS).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fnhum.2016.00055/abstract
- ^ www.fz-juelich.de/inb/inb-3/spm_anatomy_toolbox
- ^ http://marsbar.sourceforge.net
- ^ www.mathworks.com
- ^ http://icatb.sourceforge.net
- ^ www.cs.waikato.ac.nz/ml/weka/index.html
Allen, E. A., Erhardt, E. B., Damaraju, E., Gruner, W., Segall, J. M., Silva, R. F., et al. (2011). A baseline for the multivariate comparison of resting-state networks. Front. Syst. Neurosci. 5:2. doi: 10.3389/fnsys.2011.00002
Anticevic, A., Haut, K., Murray, J. D., Repovs, G., Yang, G. J., Diehl, C., et al. (2015). Association of thalamic dysconnectivity and conversion to psychosis in youth and young adults at elevated clinical risk. JAMA Psychiatry 72, 882–891. doi: 10.1001/jamapsychiatry.2015.0566
Anticevic, A., Repovs, G., and Barch, D. M. (2012). Emotion effects on attention, amygdala activation and functional connectivity in schizophrenia. Schizophr. Bull. 38, 967–980. doi: 10.1093/schbul/sbq168
Anticevic, A., Tang, Y., Cho, Y. T., Repovs, G., Cole, M. W., Savic, A., et al. (2014). Amygdala connectivity differs among chronic, early course and individuals at risk for developing schizophrenia. Schizophr. Bull. 40, 1105–1116. doi: 10.1093/schbul/sbt165
Arnold, S. J., Ivleva, E. I., Gopal, T. A., Reddy, A. P., Jeon-Slaughter, H., Sacco, C. B., et al. (2015). Hippocampal volume is reduced in schizophrenia and schizoaffective disorder but not in psychotic bipolar I disorder demonstrated by both manual tracing and automated parcellation (FreeSurfer). Schizophr. Bull. 41, 233–249. doi: 10.1093/schbul/sbu009
Bois, C., Levita, L., Ripp, I., Owens, D. C. G., Johnstone, E. C., Whalley, H. C., et al. (2015). Hippocampal, amygdala and nucleus accumbens volume in first-episode schizophrenia patients and individuals at high familial risk: a cross-sectional comparison. Schizophr. Res. 165, 45–51. doi: 10.1016/j.schres.2015.03.024
Collin, G., Hulshoff Pol, H. E., Haijma, S. V., Cahn, W., Kahn, R. S., and van den Heuvel, M. P. (2011). Impaired cerebellar functional connectivity in schizophrenia patients and their healthy siblings. Front. Psychiatry 2:73. doi: 10.3389/fpsyt.2011.00073
Demjaha, A., Murray, R. M., McGuire, P. K., Kapur, S., and Howes, O. D. (2012). Dopamine synthesis capacity in patients with treatment-resistant schizophrenia. Am. J. Psychiatry 169, 1203–1210. doi: 10.1176/appi.ajp.2012.12010144
Duan, H. F., Gan, J. L., Yang, J. M., Cheng, Z. X., Gao, C. Y., Shi, Z. J., et al. (2015). A longitudinal study on intrinsic connectivity of hippocampus associated with positive symptom in first-episode schizophrenia. Behav. Brain Res. 283, 78–86. doi: 10.1016/j.bbr.2015.01.022
Etkin, A., Prater, K. E., Schatzberg, A. F., Menon, V., and Greicius, M. D. (2009). Disrupted amygdalar subregion functional connectivity and evidence of a compensatory network in generalized anxiety disorder. Arch. Gen. Psychiatry 66, 1361–1372. doi: 10.1001/archgenpsychiatry.2009.104
Flynn, S. W., Lang, D. J., Mackay, A. L., Goghari, V., Vavasour, I. M., Whittall, K. P., et al. (2003). Abnormalities of myelination in schizophrenia detected in vivo with MRI and post-mortem with analysis of oligodendrocyte proteins. Mol. Psychiatry 8, 811–820. doi: 10.1038/sj.mp.4001337
Fornito, A., Harrison, B. J., Goodby, E., Dean, A., Ooi, C., Nathan, P. J., et al. (2013). Functional dysconnectivity of corticostriatal circuitry as a risk phenotype for psychosis. JAMA Psychiatry 70, 1143–1151. doi: 10.1001/jamapsychiatry.2013.1976
Haukvik, U. K., Westlye, L. T., Mørch-Johnsen, L., Jørgensen, K. N., Lange, E. H., Dale, A. M., et al. (2015). In vivo hippocampal subfield volumes in schizophrenia and bipolar disorder. Biol. Psychiatry 77, 581–588. doi: 10.1016/j.biopsych.2014.06.020
Howes, O. D., Kambeitz, J., Kim, E., Stahl, D., Slifstein, M., Abi-Dargham, A., et al. (2012). The nature of dopamine dysfunction in schizophrenia and what this means for treatment: meta-analysis of imaging studies. Arch. Gen. Psychiatry 69, 776–786. doi: 10.1001/archgenpsychiatry.2012.169
Kahn, I., Andrews-Hanna, J. R., Vincent, J. L., Snyder, A. Z., and Buckner, R. L. (2008). Distinct cortical anatomy linked to subregions of the medial temporal lobe revealed by intrinsic functional connectivity. J. Neurophysiol. 100, 129–139. doi: 10.1152/jn.00077.2008
Kraguljac, N. V., White, D. M., Hadley, J., Reid, M. A., and Lahti, A. C. (2014). Hippocampal-parietal dysconnectivity and glutamate abnormalities in unmedicated patients with schizophrenia. Hippocampus 24, 1524–1532. doi: 10.1002/hipo.22332
Lungu, O., Barakat, M., Laventure, S., Debas, K., Proulx, S., Luck, D., et al. (2013). The incidence and nature of cerebellar findings in schizophrenia: a quantitative review of fMRI literature. Schizophr. Bull. 39, 797–806. doi: 10.1093/schbul/sbr193
Lynall, M. E., Bassett, D. S., Kerwin, R., McKenna, P. J., Kitzbichler, M., Muller, U., et al. (2010). Functional connectivity and brain networks in schizophrenia. J. Neurosci. 30, 9477–9487. doi: 10.1523/JNEUROSCI.0333-10.2010
Manoliu, A., Meng, C., Brandl, F., Doll, A., Tahmasian, M., Scherr, M., et al. (2013). Insular dysfunction within the salience network is associated with severity of symptoms and aberrant inter-network connectivity in major depressive disorder. Front. Hum. Neurosci. 7:930. doi: 10.3389/fnhum.2013.00930
Manoliu, A., Riedl, V., Zherdin, A., Mühlau, M., Schwerthöffer, D., Scherr, M., et al. (2014). Aberrant dependence of default mode/central executive network interactions on anterior insular salience network activity in schizophrenia. Schizophr. Bull. 40, 428–437. doi: 10.1093/schbul/sbt037
Mathew, I., Gardin, T. M., Tandon, N., Eack, S., Francis, A. N., Seidman, L. J., et al. (2014). Medial temporal lobe structures and hippocampal subfields in psychotic disorders: findings from the bipolar-schizophrenia network on intermediate phenotypes (B-SNIP) study. JAMA Psychiatry 71, 769–777. doi: 10.1001/jamapsychiatry.2014.453
Mier, D., Lis, S., Zygrodnik, K., Sauer, C., Ulferts, J., Gallhofer, B., et al. (2014). Evidence for altered amygdala activation in schizophrenia in an adaptive emotion recognition task. Psychiatry Res. 221, 195–203. doi: 10.1016/j.pscychresns.2013.12.001
Millan, M. J., Fone, K., Steckler, T., and Horan, W. P. (2014). Negative symptoms of schizophrenia: clinical characteristics, pathophysiological substrates, experimental models and prospects for improved treatment. Eur. Neuropsychopharmacol. 24, 645–692. doi: 10.1016/j.euroneuro.2014.03.008
Mukherjee, P., Whalley, H. C., McKirdy, J. W., Sprengelmeyer, R., Young, A. W., McIntosh, A. M., et al. (2014). Altered amygdala connectivity within the social brain in schizophrenia. Schizophr. Bull. 40, 152–160. doi: 10.1093/schbul/sbt086
Murphy, K., Bodurka, J., and Bandettini, P. A. (2007). How long to scan? The relationship between fMRI temporal signal to noise ratio and necessary scan duration. Neuroimage 34, 565–574. doi: 10.1016/j.neuroimage.2006.09.032
Palaniyappan, L., Simmonite, M., White, T. P., Liddle, E. B., and Liddle, P. F. (2013). Neural primacy of the salience processing system in schizophrenia. Neuron 79, 814–828. doi: 10.1016/j.neuron.2013.06.027
Peters, H., Riedl, V., Manoliu, A., Scherr, M., Schwerthöffer, D., Zimmer, C., et al. (in press). Reduced functional connectivity between putamen and right anterior insula in psychotic patients with schizophrenia. Br. J. Psychiatry.
Pinkham, A. E., Liu, P., Lu, H., Kriegsman, M., Simpson, C., and Tamminga, C. (2015). Amygdala hyperactivity at rest in paranoid individuals with schizophrenia. Am. J. Psychiatry 172, 784–792. doi: 10.1176/appi.ajp.2014.14081000
Platt, J. C. (1999). “Fast training of support vector machines using sequential minimal optimization,” in Advances in Kernel Methods, eds S. Bernhard, J. Lkopf, C. B. Christopher, and J. S. Alexander (Cambridge, MA: MIT Press), 185–208.
Ragland, J. D., Ranganath, C., Harms, M. P., Barch, D. M., Gold, J. M., Layher, E., et al. (2015). Functional and neuroanatomic specificity of episodic memory dysfunction in schizophrenia: a functional magnetic resonance imaging study of the relational and item-specific encoding task. JAMA Psychiatry 72, 909–916. doi: 10.1001/jamapsychiatry.2015.0276
Rasetti, R., Mattay, V. S., Wiedholz, L. M., Kolachana, B. S., Hariri, A. R., Callicott, J. H., et al. (2009). Evidence that altered amygdala activity in schizophrenia is related to clinical. Am. J. Psychiatry 166, 216–225. doi: 10.1176/appi.ajp.2008.08020261
Sambataro, F., Blasi, G., Fazio, L., Caforio, G., Taurisano, P., Romano, R., et al. (2010). Treatment with olanzapine is associated with modulation of the default mode network in patients with Schizophrenia. Neuropsychopharmacology 35, 904–912. doi: 10.1038/npp.2009.192
Seidman, L. J., Rosso, I. M., Thermenos, H. W., Makris, N., Juelich, R., Gabrieli, J. D., et al. (2014). Medial temporal lobe default mode functioning and hippocampal structure as vulnerability indicators for schizophrenia: a MRI study of non-psychotic adolescent first-degree relatives. Schizophr. Res. 159, 426–434. doi: 10.1016/j.schres.2014.09.011
Sorg, C., Manoliu, A., Neufang, S., Myers, N., Peters, H., Schwerthöffer, D., et al. (2013). Increased intrinsic brain activity in the striatum reflects symptom dimensions in schizophrenia. Schizophr. Bull. 39, 387–395. doi: 10.1093/schbul/sbr184
Spitzer, R. L., Williams, J. B., Gibbon, M., and First, M. B. (1992). The Structured Clinical Interview for DSM-III-R (SCID). I: history, rationale and description. Arch. Gen. Psychiatry 49, 624–629. doi: 10.1001/archpsyc.1992.01820080032005
Stephan, K. E., Friston, K. J., and Frith, C. D. (2009). Dysconnection in schizophrenia: from abnormal synaptic plasticity to failures of self-monitoring. Schizophr. Bull. 35, 509–527. doi: 10.1093/schbul/sbn176
Vai, B., Sferrazza Papa, G., Poletti, S., Radaelli, D., Donnici, E., Bollettini, I., et al. (2015). Abnormal cortico-limbic connectivity during emotional processing correlates with symptom severity in schizophrenia. Eur. Psychiatry 30, 590–597. doi: 10.1016/j.eurpsy.2015.01.002
van den Heuvel, M. P., Mandl, R. C., Stam, C. J., Kahn, R. S., and Hulshoff Pol, H. E. (2010). Aberrant frontal and temporal complex network structure in schizophrenia: a graph theoretical analysis. J. Neurosci. 30, 15915–15926. doi: 10.1523/jneurosci.2874-10.2010
van Erp, T. G. M., Hibar, D. P., Rasmussen, J. M., Glahn, D. C., Pearlson, G. D., Andreassen, O. A., et al. (2015). Subcortical brain volume abnormalities in 2028 individuals with schizophrenia and 2540 healthy controls via the ENIGMA consortium. Mol. Psychiatry doi: 10.1038/mp.2015.118 [Epub ahead of print].
Wagner, G., De la Cruz, F., Schachtzabel, C., Güllmar, D., Schultz, C. C., Schlösser, R. G., et al. (2015). Structural and functional dysconnectivity of the fronto-thalamic system in schizophrenia: a DCM-DTI study. Cortex 66, 35–45. doi: 10.1016/j.cortex.2015.02.004
Wagner, G., Koch, K., Schachtzabel, C., Schultz, C. C., Gaser, C., Reichenbach, J. R., et al. (2013). Structural basis of the fronto-thalamic dysconnectivity in schizophrenia: a combined DCM-VBM study. Neuroimage Clin. 3, 95–105. doi: 10.1016/j.nicl.2013.07.010
Welsh, R. C., Chen, A. C., and Taylor, S. F. (2010). Low-frequency BOLD fluctuations demonstrate altered thalamocortical connectivity in schizophrenia. Schizophr. Bull. 36, 713–722. doi: 10.1093/schbul/sbn145
Keywords: fMRI, support vector machine, schizophrenia, functional connectivity (FC), multivariate pattern analysis, subcortical
Citation: Peters H, Shao J, Scherr M, Schwerthöffer D, Zimmer C, Förstl H, Bäuml J, Wohlschläger A, Riedl V, Koch K and Sorg C (2016) More Consistently Altered Connectivity Patterns for Cerebellum and Medial Temporal Lobes than for Amygdala and Striatum in Schizophrenia. Front. Hum. Neurosci. 10:55. doi: 10.3389/fnhum.2016.00055
Received: 16 November 2015; Accepted: 05 February 2016;
Published: 17 February 2016.
Edited by:Rachael D. Seidler, University of Michigan, USA
Reviewed by:Stephan F. Taylor, University of Michigan, USA
Jessica A. Bernard, Texas A&M University, USA
Copyright © 2016 Peters, Shao, Scherr, Schwerthöffer, Zimmer, Förstl, Bäuml, Wohlschläger, Riedl, Koch and Sorg. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution and reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kathrin Koch, firstname.lastname@example.org
† These authors have contributed equally to this work.