Step A: Cortical thickness processing. Each image was aligned to stereotaxic space, corrected for nonuniformity artifacts, tissue classified, and masked, and inner and outer cortical surfaces were extracted. The intersection between surfaces was then determined to define cortical thickness. An atlas44 was used to parcellate the cortex, and a mean thickness for each region was calculated. Step B: Correlation matrix construction. Interregional correlations of regional thickness values were assessed in each group as a measure of interregional coupling. The strength of correlation between each pairwise combination of brain regions was color coded and represented in a matrix in which the rows and columns correspond to different brain regions. Step C: Correlation matrix comparison. Statistical comparisons were done on a correlation-by-correlation basis comparing cortexwide coupling between groups. Step D: Network construction. Matrices were thresholded to retain strong positive correlations and were converted to networks in which nodes represent brain regions and connections between regions represent the strong positive correlations that survived thresholding. Step E: Graph theory analyses. The networkwide property of density was assessed as a measure of overall connectivity. The centrality of each node was assessed based on the number of connections (degree) and positioning within the network (betweenness, eigenvector centrality, and closeness). This figure is adapted from Wheeler and Voineskos48 and Voineskos et al.49 MRI indicates magnetic resonance image.
Networks were constructed by thresholding interregional correlation matrices from eFigure 1 in the Supplement and retaining only the information from the strong positive correlations (in this case, r > 0.78). In the networks, the nodes (circles) represent brain regions and the strong positive correlations that survived thresholding are shown as connections between brain regions (red lines). The size of the nodes in the networks reflects the number of connections that are incident on it (degree). Large nodes have many connections (high degree) and small nodes have few connections (low degree). A, Networks formed from groups of patients with deficit schizophrenia and nondeficit schizophrenia, and from healthy controls. The deficit schizophrenia network has more connections and more highly connected nodes compared with the networks from the other groups. B, Networks formed from patients with bipolar I disorder and matched healthy controls. The bipolar disorder network has a similar number of connections and highly connected nodes as the matched control group. AG indicates angular gyrus; CG, cingulate; CUN, cuneus; FG, fusiform gyrus; GR, gyrus rectus; IC, insular cortex; IFG, inferior frontal gyrus; inf, inferior tier; IOG, inferior occipital gyrus; ITG, inferior temporal gyrus; L, left; LG, lingual gyrus; MFG, middle frontal gyrus tier; MOG, middle occipital gyrus; MTG, middle temporal gyrus; OFG, orbitofrontal gyrus; PCUN, precuneus; PHG, parahippocampal gyrus; PreCG, precentral gyrus; PostCG, postcentral gyrus; R, right; SFG, superior frontal gyrus; sup, superior; SMG, supramarginal gyrus; SOG, superior occipital gyrus; SPG, superior parietal gyrus; and STG, superior temporal gyrus.
Density is the number of connections in the network graph divided by the number of possible connections. Density values are shown for networks formed from a range of correlation strength thresholds in each group of subjects. These thresholds correspond to Pearson correlation coefficients (r) with significance levels that reflect P values that range from P = .0001 uncorrected to P = .0001 corrected (P = .000000075). At the far left of the graphs, r = 0.78 corresponds to P = .000000075, whereas at the far right end of the graphs, r = 0.61 corresponds to P = .0001. A, Network density values in patients with deficit schizophrenia and nondeficit schizophrenia, and in healthy controls. The network density values from the deficit schizophrenia groups are higher across all thresholds than are the nondeficit schizophrenia and control groups. B, Network density values in patients with bipolar I disorder and matched healthy controls. Density values are similar in the bipolar I disorder and matched healthy control groups across the range of threshold levels.
Four measures of centrality (betweenness, closeness, eigenvector centrality, and degree) were assessed. Centrality values were computed and ranked from lowest to highest for each brain region (y-axis) that participated in the largest connected component of each thresholded network (x-axis). These thresholds correspond to Pearson correlation coefficients (r) that range from 0.61 (P = .0001 uncorrected) to 0.78 (P = .0001 corrected). Centrality measure rankings are color coded, with warm colors representing high centrality rank and cool colors indicating low centrality rank. Black color-coded cells indicate that the node was not part of the largest connected component of the network at that threshold. High-centrality brain regions were those that ranked in the top 10 in more than 50% of the thresholded networks for 3 or more of the centrality measures and are identified in boldface type. Six of these regions were also identified as having high centrality in the analysis of each separate deficit groups from the Centre for Addiction and Mental Health and Zucker Hillside Hospital (eFigure 6 in the Supplement): the left and right supramarginal gyrus and middle temporal gyrus, the left inferior frontal gyrus, and the left superior temporal gyrus. Expansion of the brain region abbreviations is presented in the Figure 2 legend.
eMethods. Details About Study Participants, Their Characterization and MRI Processing
eTable 1. LONI Probabilistic Brain Atlas (LPBA40)
eTable 2. Anatomical Automatic Labeling (AAL) Brain Atlas
eTable 3. Demographic and Clinical Characteristics of Schizophrenia Subjects and Healthy Controls From CAMH
eTable 4. Demographic and Clinical Characteristics of Schizophrenia Subjects and Healthy Controls From ZHH
eTable 5. Demographic and Clinical Characteristics of Schizophrenia Subjects and Healthy Controls From CIMH
eFigure 1. Cortex-wide Structural Coupling Comparisons
eFigure 2. Cortex-wide Structural Coupling Comparisons Between Deficit and Nondeficit Patients After Regressing Out Potential Clinical Confounds
eFigure 3. Cortex-wide Structural Coupling Comparisons Between Deficit Patients, Nondeficit Patients and Healthy Controls From the ZHH and CAMH Samples
eFigure 4. Cortex-wide Structural Coupling Comparisons Using an Alternative Parcellation Atlas (AAL)
eFigure 5. Lobular Relationships That Are Significantly Stronger in Deficit Group Compared to the Healthy Group Using Two Different Parcellation Atlases
eFigure 6. Structural Coupling Network Graphs in Deficit Patients, Nondeficit Patients and Healthy Controls From the CAMH and ZHH Samples
eFigure 7. Network Graph Densities in Deficit Patients, Nondeficit Patients and Healthy Control Groups From the CAMH and ZHH Samples
eFigure 8. High Centrality Regions in Deficit Patient Networks in CAMH and ZHH Samples
Customize your JAMA Network experience by selecting one or more topics from the list below.
Wheeler AL, Wessa M, Szeszko PR, et al. Further Neuroimaging Evidence for the Deficit Subtype of SchizophreniaA Cortical Connectomics Analysis. JAMA Psychiatry. 2015;72(5):446–455. doi:10.1001/jamapsychiatry.2014.3020
Copyright 2015 American Medical Association. All Rights Reserved. Applicable FARS/DFARS Restrictions Apply to Government Use.
The clinical heterogeneity of schizophrenia has hindered neurobiological investigations aimed at identifying neural correlates of the disorder.
To identify network-based biomarkers across the spectrum of impairment present in schizophrenia by separately evaluating individuals with deficit and nondeficit subtypes of this disorder.
Design, Setting, and Participants
A university hospital network–based neuroimaging study was conducted between February 1, 2007, and February 28, 2012. Participants included patients with schizophrenia (n = 128) and matched healthy controls (n = 130) from two academic centers and patients with bipolar I disorder (n = 39) and matched healthy controls (n = 43) from a third site. Patients with schizophrenia at each site in the top quartile on the proxy scale for the deficit syndrome were classified as having deficit schizophrenia and those in the bottom quartile were classified as having nondeficit schizophrenia.
All participants underwent magnetic resonance brain imaging.
Main Outcomes and Measures
Network-level properties of cortical thickness were assessed in each group. Interregional cortexwide coupling was compared among the groups, and graph theoretical approaches were used to assess network density and node degree, betweenness, closeness, and eigenvector centrality.
Stronger frontoparietal and frontotemporal coupling was found in patients with deficit schizophrenia compared with those with nondeficit schizophrenia (17 of 1326 pairwise relationships were significantly different, P < .05; 5% false discovery rate) and in patients with deficit schizophrenia compared with healthy controls (49 of 1326 pairwise relationships were significantly different, P < .05; 5% false discovery rate). Participants with nondeficit schizophrenia and bipolar I disorder did not show significant differences in coupling relative to those in the control groups (for both comparisons, 0 of 1326 pairwise relationships were significantly different, P > .05; 5% false discovery rate). The networks formed from patients with deficit schizophrenia demonstrated increased density of connections relative to controls and nondeficit patients (range, 0.07-0.45 in controls, 0.09-0.43 in the nondeficit group, and 0.18-0.67 in the deficit group). High centrality nodes were identified in the supramarginal, middle, and superior temporal and inferior frontal regions in deficit schizophrenia networks based on ranking of 4 centrality metrics. High centrality regions were identified as those that ranked in the top 10 in 50% or more of the thresholded networks in 3 or more of the centrality measures. Network properties were similar in patients with deficit schizophrenia from both study sites.
Conclusions and Relevance
Patients with schizophrenia at one end of a spectrum show characteristic signatures of altered intracortical relationships compared with those at the other end of that spectrum, patients with bipolar I disorder, and healthy individuals. Cortical connectomic approaches can be used to reliably identify neural signatures in clinically heterogeneous groups of patients.
In patients with schizophrenia, heterogeneity of symptoms, of response to treatment, and of functional outcome1 have hindered neurobiological investigations aimed at identifying neural correlates of this debilitating disorder.2 Negative symptoms are important predictors of functional outcome,3 and there is considerable variability in negative symptom burden among patients with schizophrenia. Kirkpatrick et al4 demonstrated that patients with the deficit form of schizophrenia are characterized by negative symptoms, including restricted affect, diminished emotional range, poverty of speech, curbing of interests, diminished sense of purpose, and diminished social drive, that are both enduring and primary. Patients with deficit schizophrenia follow a course of illness characterized by persistent impairment and poor long-term prognosis with a low likelihood of experiencing recovery, as shown in a 20-year longitudinal follow-up study.5 Evidence6 suggests that deficit and nondeficit schizophrenia are not only distinguished by symptoms and course of illness but also have distinct etiologic factors and response to treatment. In addition, early functional and structural imaging studies7-9 pointed to alterations in frontoparietal brain circuits specific to deficit schizophrenia, whereas more recent evidence10,11 points to structural alterations in the temporal lobe.
Volumetric techniques have produced conflicting results, with some studies10,11 reporting more severe volume reductions in patients with deficit schizophreniaand others12,13 reporting more severe volume reductions in the nondeficit disease. In contrast, diffusion tensor imaging studies14-16 have identified replicable white matter tract impairments in patients with deficit schizophrenia in frontoparietal and frontotemporal circuitry that is important for emotion processing, emotion expression, and socioemotional functioning. These studies provide growing evidence that relationships among the cortical gray matter regions connected by these white matter tracts are altered in patients with deficit schizophrenia. Altered patterns of network connectivity in schizophrenia have been demonstrated,17 although the results remain inconsistent. The application of graph theoretical analysis to structural network data in schizophrenia has revealed broad network inefficiency, pointing to disturbed integration among brain regions.18-21 Topologically central nodes are critical for brain network integration, and network studies22 have identified node centrality alterations frequently in prefrontal regions and less consistently in temporal and parietal regions. In addition to methodologic differences, some of the inconsistent results of these network studies are likely attributable to the considerable clinical heterogeneity among people with a diagnosis of schizophrenia.
In the present study, we used an approach that leverages the heterogeneity of schizophrenia and includes a comparator disease group with bipolar I disorder. Patients with bipolar I disorder often experience the psychotic symptoms of schizophrenia, but they are typically free of negative symptoms.23 Through assessment of cortical network properties, using structural correlation and graph theory–based approaches, we first compared patients with deficit schizophrenia and patients with nondeficit schizophrenia with healthy individuals serving as controls. Cutting across the major psychoses, we then compared people with bipolar I disorder with matched controls. We hypothesized that (1) patients with deficit schizophrenia would demonstrate alterations in network properties compared with patients with nondeficit schizophrenia and controls, (2) characteristic network differences identified in patients with deficit schizophrenia would not be present in patients with bipolar I disorder, and (3) networks in patients with deficit schizophrenia would be characterized by node alterations localized primarily to frontal and parietal brain regions.
A total of 77 patients with schizophrenia and 79 matched healthy controls were recruited at the Centre for Addiction and Mental Health (CAMH) in Toronto, Ontario, Canada; 51 patients and 51 matched controls were recruited at Zucker Hillside Hospital (ZHH) in Glen Oaks, New York; and 39 patients with euthymic bipolar I disorder and 43 matched controls were recruited from the Central Institute of Mental Health (CIMH) in Mannheim, Germany. Individuals at all sites completed the Structured Clinical Interview for Diagnostic and Statistical Manual of Mental Disorders (Fourth Revision) for Axis I disorders (SCID-I).24 Participants recruited at CIMH also completed the SCID-II for Axis II personality disorder diagnoses.25 The samples were characterized for age at onset, medication history, years of education, IQ, handedness, current clinical symptoms, and current substance use or past dependence (detailed inclusion/exclusion criteria for each sample are provided in the eMethods in the Supplement). After complete description of the study to the participants, written informed consent was obtained as per local institutional review board regulations. Institutional review board approval was obtained from the CAMH, the ZHH, and the CIMH. Participants received financial compensation.
The Positive and Negative Syndrome Scale26 and the Brief Psychiatric Rating Scale27 were used to characterize clinical symptoms at CAMH and ZHH, respectively. The deficit and nondeficit schizophrenia groups were selected from the full schizophrenia samples based on characterization of the deficit syndrome according to the proxy for the deficit syndrome.28 Studies have examined the stability of the deficit syndrome classification using the proxy for the deficit syndrome and have shown the characterization to be stable over the short term (24 months)29 and across serial assessments over a 20-year follow-up period5 (eMethods in the Supplement). The proxy for the deficit syndrome has good specificity, sensitivity, and accuracy for identification of patients with deficit schizophrenia28 in both early-episode and chronic populations.30-35 Specifically, the proxy for the deficit syndrome is defined as the sum of the scores from either the Positive and Negative Syndrome Scale or the Brief Psychiatric Rating Scale of the anxiety, guilt feelings, depressive mood, and hostility items subtracted from the blunted affect item score. This calculation reflects the lack of negative affect and lack of dysphoria that are characteristic of patients with deficit schizophrenia along with their high degree of negative symptoms.4,28,36 The top quartile of proxy case identification scores was used to identify patients with deficit schizophrenia in each sample separately.37 To enhance the likelihood of correct classification, patients with nondeficit schizophrenia were identified as those who ranked in the bottom quartile of deficit scores in each sample. Healthy control participants from each site were selected based on matching to the deficit and nondeficit groups for age, sex, handedness, and parental education.
The neuroimaging aspect of the study was conducted between February 1, 2007, and February 28, 2012. Magnetic resonance (MR) images were acquired for all participants.
The MR images at CAMH were acquired using an 8-channel head coil on a 1.5-T system (EchoSpeed; General Electric Medical Systems). Axial inversion recovery–prepared spoiled gradient recall images (n = 124) were acquired using a 1.5-mm-thick slice acquisition with the following image parameters: echo time (TE), 5.3 milliseconds; repetition time (TR), 12.3 milliseconds; time to inversion, 300.0 milliseconds; and flip angle, 20°.
The MR imaging examinations at ZHH were conducted on a 3-T whole-body superconducting imaging system (Signa HDx; General Electric Medical Systems). The spoiled gradient recall images (n = 216) were acquired using a 1-mm-thick slice acquisition with the following image parameters: TE, 3 milliseconds; TR, 7.5 milliseconds; matrix size, 256 × 256; and field of view, 240 mm.
The MRI data at CIMH were acquired on a 3-T scanner (Magnetom Trio; Siemens Medical Solutions). The T1-weighted images (n = 160) were acquired using a 1.1-mm-thick slice acquisition with the following image parameters: TE, 2.98 milliseconds; TR, 2300 milliseconds; matrix size, 256 × 240 × 160; and field of view, 256 × 240 × 176 mm.
To calculate vertexwise cortical thickness, the T1-weighted images were submitted to the CIVET pipeline, version 1.1.10 (Montreal Neurological Institute, McGill University).38 Briefly, images were registered to a nonlinear template, with inhomogeneity corrected,39 skull stripped, and tissue classified.40,41 Deformable models were used to create gray and white matter surfaces,42 and white to gray matter surface distances were determined using the t-link metric.38 The thickness data were subsequently blurred using a 20-mm surface–based diffusion blurring kernel and nonlinearly aligned to a template43 (eMethods in the Supplement).
For each participant, vertexwise cortical thickness maps were parcellated using cortical regions defined in the Probabilistic Brain Atlas44 (http://www.loni.usc.edu/atlases/) with additional manual segmentation of the middle frontal gyri (eMethods in the Supplement). Mean cortical thickness was computed in a total of 52 distinct cortical regions (eTable 1 in the Supplement). For comparison, an alternative parcellation atlas with 74 cortical regions was also used (eTable 2 in the Supplement).45 Between-subject Pearson correlation coefficients (r) were computed for values in all pairwise combinations of brain regions to assess intracortical coupling. A linear regression was performed at every cortical region to remove the effects of site in the combined diagnostic groups from CAMH and ZHH (32 deficit, 32 nondeficit, and32 controls), and the correlations among resulting residuals were computed. Correlations in cortical thickness were examined in a sample of 32 patients with bipolar I disorder and 32 matched healthy controls from the CIMH sample. Correlation matrix construction and analyses were performed using MATLAB software (MathWorks Inc, version 126.96.36.1991 [R2013b]).
To evaluate intracortical coupling differences between diagnostic groups, all individual correlations were compared with no a priori predictions (1326 comparisons). Differences were assessed between the following pairwise combinations of groups: deficit schizophrenia vs schizophrenia-matched healthy controls; deficit schizophrenia vs nondeficit schizophrenia; nondeficit schizophrenia vs schizophrenia-matched healthy controls, and bipolar I disorder vs bipolar I disorder–matched healthy controls. Significant differences in individual correlations were determined with the Fisher r to z transformation. This procedure converts the r values to z values with the Fisher transform and compares the difference with the normal distribution.46 The false discovery rate was set at 5%.
Binary networks were created by thresholding cortexwide interregional cortical thickness correlations in each study group. The nodes in the networks represent brain regions, and the correlations that survived thresholding were considered connections. Graph theoretical measures (network density and node centrality) were assessed in networks formed at an array of correlation strength thresholds. The thresholds used ranged from a nominal value of P = .0001 to P = .000000075 (corresponding to a Bonferroni correction that accounts for all pairwise comparisons). This P value range corresponds to an r value range of 0.61 to 0.78. Networks were visualized with Cytoscape, version 3.1.10 (http://www.cytoscape.org.).
Network density is defined as the proportion of possible connections present in the network and is a measure of the overall level of connectivity. Network density was computed in all 5 groups at each network threshold level.
Node centrality was characterized in groups with coupling and density distinct from that in the control groups. To characterize centrality, 4 measures (degree, eigenvector centrality, closeness, and betweenness) were computed for each node that was part of the largest connected component of each network. Degree corresponds to the number of connections that are incident on a node. Betweenness is the number of all shortest paths in the network that pass through a given node. Closeness is the inverse of the mean shortest path length from one node to all other nodes in the network. Eigenvector centrality assigns relative scores to all nodes in the network based on the concept that connections to high-scoring nodes contribute more to the score than do equal connections to low-scoring nodes. All regions were ranked by degree, eigenvector centrality, closeness, and betweenness in each thresholded network. High-centrality brain regions were classified as those that ranked in the top 10 in more than 50% of the thresholded networks for 3 or more of the centrality measures. All graph theory measures were computed with functions from the Brain Connectivity Toolbox (http://www.brain-connectivity-toolbox.net/).47 Image processing and analysis steps are summarized in Figure 1.
All analyses were also performed in the individual samples of deficit schizophrenia, nondeficit schizophrenia, and matched controls from CAMH and ZHH. The additional evaluation was performed to determine whether findings were similar at both sites.
In the CAMH and ZHH samples, no significant differences were present among the healthy control, deficit schizophrenia, and nondeficit schizophrenia groups on any of the demographic variables, and no significant differences were detected between the deficit and nondeficit groups in age at onset or antipsychotic medication load (all P > .05) (eTable 3 and eTable 4 in the Supplement). Consistent with the classification of the deficit syndrome, patients with deficit schizophrenia displayed higher levels of blunted affect and lower levels of anxiety, feelings of guilt, depressive mood, and hostility as well as similar positive symptom scores in the CAMH and ZHH samples (eTable 3 and eTable 4 in the Supplement, respectively). Patients with bipolar I disorder did not differ significantly from the matched healthy controls in any demographic variables (eTable 5 in the Supplement). A summary of participant characteristics is presented in the Table.
Structural coupling comparisons revealed 49 relationships that were stronger in the patients with deficit schizophrenia than in the healthy controls and 17 relationships that were stronger in the deficit schizophrenia group than in the nondeficit schizophrenia group (all P < .05, corrected). Most of these differences were in frontoparietal and frontotemporal relationships. There were no significant differences between the patients with nondeficit schizophrenia and healthy controls (all P > .05, corrected) (eFigure 1A in the Supplement) or between patients with bipolar I disorder and matched healthy controls (all P > .05, corrected) (eFigure 1B in the Supplement). Similar frontoparietal and frontotemporal differences were also identified in the deficit schizophrenia group after regressing out potential clinical confounders (eFigure 2 in the Supplement) and in the separate CAMH and ZHH samples (eFigure 3 in the Supplement) as well as when using an alternative cortical parcellation atlas (eTable 2, eFigure 4, and eFigure 5 in the Supplement). Comparison of the intracortical coupling relationships between the control group matched to the schizophrenia groups and the control group matched to bipolar I disorder group revealed that only 1 correlation was significantly different between the 2 control groups.
Networks were constructed with brain regions as nodes and strong interregional cortical thickness correlations as connections (Figure 2). The density of network connections was greater in the deficit schizophrenia group than in the nondeficit schizophrenia group and healthy control group at the range of correlation thresholds explored (range, 0.07-0.45 in controls, 0.09-0.43 in the nondeficit group, and 0.18-0.67 in the deficit group) (Figure 3A). Networks constructed in the separate CAMH and ZHH samples (eFigure 6 in the Supplement) each demonstrated increased density of connections in the deficit groups compared with the nondeficit and healthy control groups (CAMH: range, 0.10-0.43 in controls, 0.10-0.48 in the nondeficit group, and 0.19-0.63 in the deficit group; ZHH: range, 0.09-0.33 in controls, 0.09-0.33 in the nondeficit group, and 0.31-0.68 in the deficit group) (eFigure 7 in the Supplement). Network density was similar in patients with bipolar I disorder and matched healthy controls (range, 0.08-0.43 in controls and 0.08-0.42 in the bipolar I disorder group) (Figure 3B) and in the schizophrenia- and bipolar I disorder–matched control groups.
Node centrality was characterized in the group of patients with deficit schizophrenia who had intracortical coupling and network density different from that of healthy controls. Assessing the 4 metrics of node centrality (degree, eigenvector centrality, closeness, and betweenness) in each brain region identified high centrality regions in networks of the deficit schizophrenia group . These regions were the inferior frontal gyrus, supramarginal gyrus, superior temporal gyrus, and middle temporal gyrus bilaterally and the caudal middle frontal gyrus, superior frontal gyrus, superior parietal gyrus, and fusiform gyrus in the right hemisphere (Figure 4). Six of these regions were also identified as having high centrality in each of the separate CAMH and ZHH samples of patients with deficit schizophrenia: the left and right supramarginal gyrus, left and right middle temporal gyrus, left inferior frontal gyrus, and left superior temporal gyrus (eFigure 8 in the Supplement).
Using connectomic approaches, we found a highly correlated structural network in patients with the deficit subtype of schizophrenia defined by high node centrality in inferior frontal, inferior parietal, and middle and superior temporal regions. Intracortical coupling among these regions was significantly stronger in people with deficit schizophrenia than in patients with nondeficit schizophrenia and controls. Network properties were similar in patients with deficit schizophrenia from both the CAMH and ZHH sites. When cortexwide structural coupling was compared in patients with bipolar I disorder and healthy controls, no significant differences were found. To our knowledge, this study provides the first graph theory–based results linking clinical heterogeneity with heterogeneity in neuroimaging findings in schizophrenia.
In contrast to the network-based approach applied in the present study, regional cortical thickness reductions are present in both the deficit and nondeficit schizophrenia groups compared with the control group14 and are present in bipolar I disorder in many of the same regions as in patients with schizophrenia.50 These findings suggest that conventional regional approaches may not be useful for identifying circuit- or network-based impairment.14 Previous network-based analyses17-20,51 in patients with chronic schizophrenia have reported deficits in frontotemporal and frontoparietal connectivity and altered centrality of regions. We found enhanced structural coupling associated with high regional centrality in the inferior frontal, inferior parietal, and middle and superior temporal cortices. Many of these same regions are connected by the white matter tracts (ie, inferior longitudinal fasciculus, arcuate fasciculus, and uncinate fasciculus) that have been identified14-16 as disrupted in patients with deficit schizophrenia. Published data52-54 support the hypothesis that regions that show strong structural coupling tend to be connected by white matter tracts. Although the exact relationship between disrupted white matter microstructure and enhance structural covariance cannot be elucidated here, both increased and decreased coupling between regions connected by white matter tracts has been reported55-57 in other neurodevelopmental disorders. When our present findings are considered with previous findings on white matter circuitry, the evidence is compelling that this subgroup of patients with schizophrenia is characterized by distinct impairments in neural circuitry. Furthermore, the possibility that recent network findings in schizophrenia18-21 are driven by this subset of patients deserves consideration.
Patients with deficit schizophrenia may represent a subgroup of people with schizophrenia who experience an early-onset, nonprogressive developmental process that interferes with basic cognitive and social skills beginning in early childhood.13,58 The increased density in structural coupling networks that we found may reflect altered early neurodevelopment in the deficit group. In support of this interpretation, a longitudinal developmental neuroimaging study59 in healthy individuals demonstrated that brain regions with synchronized developmental change show strong convergence with regions that demonstrate cross-sectional covariance of cortical thickness. Thus, increased density in the deficit schizophrenia networks found in the present study may reflect decreased differentiation of brain regions during development, which constitutes a disruption in healthy synchronized brain maturation. This diminished differentiation of brain regions may lead to reduced efficiency and increased cost for brain networks that strive for a balance between cost and efficiency.60 The high levels of centrality that we found in supramarginal, inferior frontal, superior, and middle temporal gyri in patients with deficit schizophrenia offer the interpretation that these regions may have undergone reduced differentiation more so than other brain regions. Reduced specialization of these regions may contribute to impairment in complex emergent properties that require module integration, such as social cognition, global cognition, and language, which are impaired in people with deficit schizophrenia.61
The impaired networks identified in patients with deficit schizophrenia are important for processes related to negative symptoms. For example, frontoparietal circuitry in the right hemisphere comprises the mirror neuron system, which is activated during basic emotion understanding and emotion experience sharing.62 These social cognitive domains are considerably impaired in people with schizophrenia (particularly those with the deficit form61), are an important determinant of functioning,63 and appear to be largely preserved in patients with bipolar I disorder.64 Future work characterizing cortical brain circuitry and social cognitive performance across patients with schizophrenia or bipolar I disorder can help further clarify whether the impairments observed in this study in patients with deficit schizophrenia are related to social cognitive impairments in the same group of patients.
One limitation of our study is that patient samples were collected at 3 academic centers and image acquisition parameters were different at the 3 sites. Combining schizophrenia groups from the CAMH and ZHH samples provided the power necessary to examine networkwide differences in deficit subtype patients, which are difficult to collect in large numbers. Additionally analyzing the separate samples revealed similar network properties in the deficit schizophrenia patient groups from both sites. A second limitation is that some have argued that regional reductions in brain structure may be an epiphenomenon of antipsychotic medication effects.65 Although it is possible that medication effects influenced the circuitry findings presented here, antipsychotic medication exposure was not different between the deficit and nondeficit schizophrenia groups at either site. Furthermore, although both approaches have been shown to accurately classify deficit and nondeficit subtypes, the Positive and Negative Syndrome Scale at CAMH and the Brief Psychiatric Rating Scale at ZHH were used to classify people with schizophrenia at each site. Finally, although use of an alternative cortical parcellation atlas produced similar results, some of the specific relationships that were found to be significantly different between groups were dependent on the template atlas used.
Our study may resolve some inconsistency in neuroimaging findings2,66 by showing that network properties of nodes within frontoparietal circuits are not characteristic of all patients with schizophrenia but rather are a neural signature of those with deficit schizophrenia. In addition, we showed that patients with bipolar I disorder do not demonstrate altered structural coupling among the thickness of cortical brain regions compared with control samples; to our knowledge, this is the first such assessment. Targeting these circuits should assist in novel therapeutic approaches aimed at treating negative symptoms and social impairment in people with schizophrenia.
Submitted for Publication: June 8, 2014; final revision received November 7, 2014; accepted November 9, 2014.
Corresponding Author: Aristotle N. Voineskos, MD, PhD, Research Imaging Centre, Centre for Addiction and Mental Health, 250 College St, Toronto, ON M5T 1R8, Canada (email@example.com).
Published Online: March 18, 2015. doi:10.1001/jamapsychiatry.2014.3020.
Author Contributions: Drs Wheeler and Voineskos had full access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis.
Study concept and design: Wheeler, Szeszko, Lerch, Voineskos.
Acquisition, analysis, or interpretation of data: All authors.
Drafting of the manuscript: Wheeler, Voineskos.
Critical revision of the manuscript for important intellectual content: All authors.
Statistical analysis: Wheeler, Chakravarty, Lerch, Linke, Voineskos.
Obtained funding: Wheeler, Wessa, Chakravarty, Mulsant, Malhotra, Voineskos.
Administrative, technical, or material support: Szeszko, Foussias, Chakravarty, DeRosse, Remington, Malhotra, Voineskos.
Study supervision: Szeszko, Malhotra, Voineskos.
Conflict of Interest Disclosures: Dr Foussias has served on advisory boards for Roche and has received speaker fees from Lundbeck, Novartis, and Roche. Dr Remington has received consultancy fees from Neurocrine, Novartis, and Synchroneuron as well as research and speaker fees from Novartis. Dr Mulsant has received research support in the form of medications for National Institutes of Health–funded clinical trials from Bristol-Myers Squibb, Eli Lilly and Company, and Pfizer. Dr Malhotra has served as a consultant for Genomind Inc. No other disclosures were reported.
Funding/Support: Work from Zucker Hillside Hospital was supported by grant P50MH080173 from the National Institute of Mental Health (NIMH). Work from the Centre for Addiction and Mental Health (CAMH) was supported in part by the CAMH Foundation thanks to the Kimel Family, Koerner New Scientist Award, and Paul E. Garfinkel New Investigator Catalyst Award, as well as the Brain and Behavior Research Foundation, Canadian Institutes of Health Research, Ontario Mental Health Foundation, and grant R01MH099167 from the NIMH. Work from Heidelberg was supported by the Deutsche Forschungsgemeinschaft (grant SFB636/C6, We3638/3-1).
Role of the Funder/Sponsor: The funding organizations had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Correction: This article was corrected on March 25, 2015, to fix a caption.