Lateral, medial, and ventral surfaces of the left hemisphere with average folding patterns of sulci (dark gray) and gyri (light gray), derived using the surface-based morphing procedure.20,31 Temporal-prefrontal regions of interest (ROIs) are shown in red. Lateral surface (A) includes dorsolateral prefrontal cortex (middle frontal gyrus and superior frontal sulcus), inferior frontal gyrus, and superior temporal gyrus. Medial surface (B) includes parahippocampal gyrus and anterior cingulate gyrus. Ventral surface (C) includes orbitofrontal cortex (H-shaped orbital sulcus, arcuate orbital sulcus, medial orbital sulcus, and fronto-orbital gyrus). Control ROIs are shown in blue. Lateral surface(A) includes central sulcus, precentral and postcentral gyri, the intraparietal sulcus, postcentral sulcus, and the superior parietal lobule. Medial surface(B) includes the calcarine fissure. Lowercase a indicates anterior; AG, angular gyrus; AOS, anterior occipital sulcus; ArOrS, arcuate orbital sulcus; AS, angular sulcus; CalcS, calcarine sulcus; CeS, central sulcus; CGa, anterior cingulate gyrus; CGp, posterior cingulate gyrus; CN, cuneus; CoS, collateral sulcus; CSa, anterior cingulate sulcus; CSI, circular sulcus of insula; CSp, posterior cingulate sulcus; FG, fusiform gyrus; FOrG, fronto-orbital gyrus; FP, frontal pole; H-OrS, H-shaped orbital sulcus; Hscl, Heschl's gyrus; IFG, inferior frontal gyrus; IFS, inferior frontal sulcus; Ins, insula; IPS, intraparietal sulcus; ITG, inferior temporal gyrus; ITS, inferior temporal sulcus; LG, lingual gyrus; LOG, lateral occipital gyrus; LOS, lateral occipital sulcus; MFG, middle frontal gyrus; MOrS, medial orbital sulcus; MTG, middle temporal gyrus; OLi, inferior occipital lateral gyrus; OP, occipital pole; OTS, occipitotemporal sulcus; lowercase p, posterior; PaC, paracingulate cortex; PCN, precuneus; PH, parahippocampal gyrus; PoCS, postcentral sulcus; PoG, postcentral gyrus; POS, parieto-occipital sulcus; PrCS, precentral sulcus; PrG, precentral gyrus; SFG, superior frontal gyrus; SFS, superior frontal sulcus; SMG, supramarginal gyrus; SOrS, superior orbital sulcus; SPL, superior parietal lobule; SPS, subparietal sulcus; STG, superior temporal gyrus; STS, superior temporal sulcus; SyF, sylvian fissure; and TP, temporal pole.
Cortical reconstruction procedures. A, Sagittal sections through individual brains (2 control subjects, left; 2 patients, right) with representations of the gray matter–white matter boundary (yellow line) and the pial surface (red line). B, Inflated cortical surfaces showing sulcal-gyral folding patterns of the same 4 individual brains. C, Each subject's reconstructed brain was morphed to a template guided by the mean and variance of the folding pattern at each point in the cortex,20,31 and regions of interest (ROIs) were painted onto this common surface template (the left inferior frontal gyrus, left superior temporal gyrus, and left central sulcus are illustrated in yellow). D, The transformation was reversed such that the ROIs were mapped back onto the cortical surfaces of each individual participant.
Statistical thickness difference maps from data sets 1 (A) and 2 (B) and the combined data sets (C). In mean thickness difference maps (D), yellow corresponds to mean thickness difference of at least 0.2 mm (percentage reduction in thickness of at least 8.7%). Cortical thickness differences are displayed on the average curvature (see Figure 1), with sulcal and gyral folds shown in light and dark gray, respectively. R indicates right; L, left.
Manual validation of automated thickness measures. A, Two small regions shown in red within the right orbitofrontal cortex (left panel) and the left calcarine sulcus (right panel) that did and did not show significant between-group differences, respectively, on the statistical thickness difference maps (Figure 3). B, These 2 regions were each mapped back onto the corresponding brain volumes of each of the 32 individual subjects in data set 2 using the spherical morphing algorithm applied in reverse (1 patient's brain volume is shown). Manual thickness measures of these 2 regions in each individual's volume were estimated by taking the shortest distance (over 3 cross-sectional views) between a visual estimate of the gray-white border and the pial surface, averaged across 10 points within each region. C, Mean and SDs of manual estimates of thickness were similar to the automated estimates of thickness in both patient and control groups. There was a significant difference between patients and control subjects in manual estimates of thickness in the right orbitofrontal cortex (t = 4.63; P<.0001) but not in the left calcarine fissure (t = 1.50; P =.14).
Kuperberg GR, Broome MR, McGuire PK, David AS, Eddy M, Ozawa F, Goff D, West WC, Williams SCR, van der Kouwe AJW, Salat DH, Dale AM, Fischl B. Regionally Localized Thinning of the Cerebral Cortex in Schizophrenia. Arch Gen Psychiatry. 2003;60(9):878–888. doi:10.1001/archpsyc.60.9.878
Schizophrenia is characterized by small reductions in cortical gray matter volume, particularly in the temporal and prefrontal cortices. The question of whether cortical thickness is reduced in schizophrenia has not been addressed using magnetic resonance imaging (MRI) techniques. Our objectives were to test the hypothesis that cortical thinning in patients with schizophrenia (relative to control subjects) is greater in temporal and prefrontal regions of interest (ROIs) than in control ROIs (superior parietal, calcarine, postcentral, central, and precentral cortices), and to obtain an unbiased estimate of the distribution of cortical thinning in patients (relative to controls) by constructing mean and statistical cortical thickness difference maps.
Participants included 33 right-handed outpatients receiving medication and meeting DSM-IV criteria for schizophrenia and 32 healthy volunteers, matched on age and parental socioeconomic status. After high-resolution MRI scans, models of the gray-white and pial surfaces were generated for each individual's cortex, and the distance between these 2 surfaces was used to compute cortical thickness. A surface-based averaging technique that aligned the main cortical folds across individuals allowed between-group comparisons of thickness within ROIs, and at multiple, uniformly sampled loci across the cortical ribbon.
Relative to controls, patients showed greater cortical thinning in temporal-prefrontal ROIs than in control ROIs, as revealed by a significant (P<.009) interaction between group and region type. Cortical thickness difference maps revealed significant (at P<.05, corrected) thinning within the orbitofrontal cortices bilaterally; the inferior frontal, inferior temporal, and occipitotemporal cortices on the left; and within the medial temporal and medial frontal cortices on the right. Superior parietal and primary somatosensory and motor cortices were relatively spared, even at subthreshold significance levels.
Patients with chronic schizophrenia showed widespread cortical thinning that particularly affected the prefrontal and temporal cortices. This thinning might reflect underlying neuropathological abnormalities in cortical structure.
SCHIZOPHRENIA IS characterized by subtle but significant gray matter deficits, particularly affecting the prefrontal and temporal cortices.1 Morphometric studies have focused on the measurement of cortical volume. The human cerebral cortex, however, is a 6-layered sheet of neural tissue that has been extensively folded, such that approximately two thirds is buried within cortical sulci. We identified sites of cortical thinning in schizophrenia using an automated morphometric method that yields measures of cortical gray matter thickness.
Morphometric studies examining cortical volume in schizophrenia have traditionally focused on regions of interest (ROIs) that are selected on the basis of a priori hypotheses. Although ROI studies have the advantage of being theoretically motivated, they are labor-intensive and time-consuming, and therefore usually do not examine the whole brain (with few exceptions, eg, Goldstein et al2). To complement the ROI approach, there have been attempts to automate the measurement of gray matter volume or density throughout the brain.3- 6 So that homologous regions can be averaged within and between subject groups, individual brains must be aligned, by registering to standardized volumetric space3,4,7 or by using computational matching strategies that align corresponding locations on the cortical surface.5,6 Maps showing mean and statistical differences in volume or density between groups can then be constructed. This approach (hereafter termed the morphometric brainmap approach) has the advantage of giving an unbiased and comprehensive assessment of morphometric changes at uniformly sampled loci throughout the brain. Taken together, morphometric brainmap studies in schizophrenia have confirmed subtle volumetric reductions in multiple anatomical regions within the prefrontal and temporal cortices1,8- 10 and within other regions such as the insula,8- 10 the occipitotemporal cortex,10 and, in patients with early-onset schizophrenia, the parietal cortex.5 This picture of subtle, widespread morphometric changes is consistent with a growing understanding of schizophrenia as a disorder characterized by cognitive abnormalities in a number of domains, with associated functional abnormalities that affect multiple neural systems, particularly in the prefrontal and temporal cortices.11
The question of whether cortical thickness is reduced in patients with schizophrenia has, thus far, not been addressed using magnetic resonance imaging (MRI) techniques. This is an important question because abnormalities in cortical thickness might reflect underlying neuropathological abnormalities in the intrinsic structure and integrity of cortical laminae.12,13 Some postmortem studies in schizophrenia report modest thinning of the cortical ribbon and/or layer-specific decreases in cell density in localized regions, such as within the dorsolateral prefrontal cortex.14- 16 Manual postmortem measurements of thickness across the entire cortical surface on large numbers of individuals, however, would be too time-consuming to be feasible. Similarly, it would also be extremely time-consuming and technically very difficult to measure thickness manually across the cortical surface from slice MR data.17
In the present study, we used automated, objective methods17- 21 to address the question of whether the thickness of the cortical gray matter is reduced in patients with schizophrenia and, if so, to determine the regional distribution of such thinning. Volume is a product of thickness and surface area. Hence, we predicted that cortical thinning might contribute to the volumetric changes in temporal and prefrontal cortices described in previous morphometric studies of schizophrenia. We computed the average cortical thickness in several ROIs in individual subjects to test our specific hypothesis that thinning would be greater in patients (relative to control subjects) in temporal and prefrontal ROIs than in control ROIs (the superior parietal cortex and primary visual, somatosensory, and motor cortices). In addition, we adopted a morphometric brainmap approach by constructing cortical thickness difference maps that gave an unbiased estimate of mean and statistical differences in thickness between patients and controls at multiple, uniformly sampled loci across the entire cortical ribbon.
Thirty-three subjects (17 patients and 16 matched controls [data set 1]) underwent scanning on a 1.5-T GE Signa NV/i system scanner (General Electric Company, Milwaukee, Wis), and an additional 32 subjects (16 patients and 16 matched controls [data set 2]) underwent scanning on a 1.5-T Siemens scanner (Siemens Medical Solutions, Iselin, NJ). Twenty-four subjects from data set 1 were recruited in London, England (12 patients from the Maudsley and Bethlem Royal Hospitals and 12 controls by advertisement); the remaining 9 subjects of data set 1 were recruited in Boston, Mass (5 patients from the Erich Lindemann Mental Health Center in Massachusetts, and 4 controls by advertisement). One subject from data set 1 was subsequently excluded, leaving 32 subjects (16 patients and 16 matched controls) in data set 1. All subjects in data set 2 were recruited in Boston (16 patients from the Erich Lindemann Mental Health Center and 16 controls by advertisement). All participants were right-handed as assessed using the modified Edinburgh Handedness Inventory22,23 and were aged 19 to 63 years. Exclusion criteria for all participants included neurological disease or damage, head trauma with documented cognitive sequelae or loss of consciousness of more than 5 minutes, medical disorders that can impair neurocognitive function, a history of substance abuse or dependence, and specific exclusions for MRI. Written consent was obtained from all subjects before participation according to the established guidelines of the Maudsley and Bethlem Royal Hospitals and Massachusetts General Hospital human subjects research committees.
Diagnoses were made by staff psychiatrists and were confirmed by 2 psychiatrists(G.R.K. and M.R.B.) using a structured clinical interview24 and examination of the case notes. All patients met DSM-IV criteria for schizophrenia25 and were receiving stable doses of antipsychotic medication. One patient was taking anticholinergic medication (psychopathological details are shown in Table 1). Healthy volunteers were not taking any medication and were screened24 to exclude the presence of psychiatric disorders. Patients and controls matched closely on demographic characteristics (Table 1).
Subjects in data set 1 underwent 2 structural scans on a 1.5-T GE Signa NV/i System at the Maudsley Hospital, London, or the Athinoula A. Martinos Center for Biomedical Imaging, Charlestown, Mass. Each scan consisted of a radiofrequency spoiled gradient refocused acquisition in the steady state[rf-spoiled GRASS] sequence with 124 sagittal slices; 1.3-mm thickness; repetition time (TR), 24 milliseconds; echo time (TE), minimum for symmetric echo; flip angle, 30°; bandwidth, 16 kHz and in-plane resolution, 1.3 × 1 mm. Scanning sequences and protocols were identical in the 2 sites.
Subjects in data set 2 underwent 2 structural scans on a 1.5-T Siemens scanner at the Athinoula A. Martinos Center for Biomedical Imaging, Boston. Each scan consisted of a 3-dimensional magnetization prepared rapid acquisition gradient-echo (3-D MPRAGE) sequence (128 sagittal slices; 1.3-mm thickness; TR, 7.25 milliseconds; TE, 3 milliseconds; flip angle, 7°; bandwidth, 195 Hz per pixel; in-plane resolution, 1.3 × 1 mm).
The 2 structural scans for each participant were averaged, after motion correction, to create a single high signal-to-noise average volume, thus avoiding confounds of poor scan quality or low signal. The resulting volume was used to segment cerebral white matter18,19 and to estimate the gray-white interface. Topological defects in the gray-white estimate were fixed,21 and this gray-white estimate was used as the starting point for a deformable surface algorithm designed to find the pial surface with submillimeter precision.17 The entire cortex in each individual subject was then visually inspected, and any inaccuracies in segmentation were manually corrected. One subject's data (from data set 1) was discarded owing to image artifacts leading to inaccurate surfaces.
For each subject, thickness measures across the cortex were computed by finding the point on the gray-white boundary surface that was closest to a given point on the estimated pial surface (and vice versa) and averaging between these 2 values.17 The accuracy of the thickness measures derived from this technique have been validated by direct comparisons with manual measures on postmortem brains.30
The surface representing the gray-white border was inflated,18,20 differences between individuals in the depth of gyri/sulci were normalized, and each subject's reconstructed brain was then morphed and registered to an average spherical surface representation that optimally aligned sulcal and gyral features across subjects.20,31 This spherical morphing procedure was used to compute the average thickness within ROIs and to construct the cortical thickness difference brain maps.
We selected the following 4 anatomical regions in the frontal lobe and 2 anatomical regions in the temporal lobe that we hypothesized would show group differences: the inferior frontal gyrus, the dorsolateral prefrontal cortex (broadly defined as the middle frontal gyrus plus superior frontal sulcus), the anterior cingulate gyrus, the orbitofrontal cortex (defined as the fronto-orbital gyrus plus orbital sulci), the parahippocampal gyrus, and the superior temporal gyrus. In addition, we selected the following 6 control ROIs that we predicted would not show significant difference in cortical thickness between groups: the central sulcus, the precentral gyrus, the postcentral gyrus, the postcentral sulcus and intraparietal lobule, the superior parietal lobule, and the calcarine fissure. These ROIs were drawn on maps of average folding patterns (derived using the surface-based morphing procedure as described by Fischl et al20,31) on the cortical surface of each hemisphere, with reference to an anatomical atlas32 (Figure 1).
Each ROI was mapped back on to each individual subject's unfolded surface by applying the same algorithm that morphed each subject's unfolded surface to the average spherical surface representation20,31 in reverse (Figure 2).
The average thickness within each ROI in each hemisphere in each individual subject was the dependent variable in a repeated-measures analysis of variance, with subjects as a random effect and group as a between-subject factor (2 levels: patients and controls). Region type (2 levels: temporal-prefrontal and control regions) and hemisphere (2 levels: right and left) were within-subject factors. Data set was entered into the model as a potentially confounding covariate. Analyses were repeated with the average thickness of each individual over the entire cortical surface as an additional covariate. We used a standard level of α = .05 to determine whether effects reached significance.
We used the spherical transformation to map the thickness measurements at each vertex on each subject's cortical surface into a common spherical coordinate system.20,31 The data were smoothed on the surface tessellation using an iterative nearest-neighbor averaging procedure (1000 iterations were applied, equivalent to applying a 2-dimensional gaussian smoothing kernel along the cortical surface with a full-width half-maximum of approximately 6 cm). Data were then averaged across participants in the spherical coordinate system31 to generate mean maps showing the mean differences in thickness (in millimeters) between the patient and control groups at multiple uniformly sampled loci across the cortical surface (Figure 3D).
Statistical thickness difference maps were constructed using a t statistic (Figure 3A-C). We used a general linear model in which the main effects of group (thickness differences between patients and controls) are shown, covarying for age, sex, and years of education. For the map in which data sets 1 and 2 were combined(Figure 3C), data set was an additional covariate. We had strong a priori hypotheses that differences would be observed within the prefrontal and temporal cortices. Given their surface area and the degree of surface smoothing applied, 42 effective tests were conducted within these temporal-prefrontal regions. Applying a Bonferroni correction, an uncorrected threshold of P<.001 is approximately equivalent to a corrected threshold of P<.05 (ie, 5% expected false-positive findings). The regions displayed in yellow on the brainmaps all met this criterion for significance. We also display thickness differences at subthreshold levels (an uncorrected threshold of P<.05) in red. This was to examine the full extent of cortical thickness differences across the cortical surface, particularly in data sets 1 and 2 that were collected independently on different scanners, and to determine whether regions outside temporal-prefrontal regions showed any hint of thinning in patients with schizophrenia.
Automated estimates of thickness in 2 regions that did and did not show significant (P<.001) between-group differences on the statistical thickness difference maps were validated to be reliable to less than 0.5 mm by direct comparisons with manual measures made on slice MR data (Figure 4).
Table 2 presents the mean and percentage differences in cortical thickness between patients and controls in all of these ROIs (averaged across the right and left hemispheres). An analysis of variance showed a significant interaction between group (patients vs controls) and region type (temporal-prefrontal regions vs control ROIs)(F1,61 = 7.23; P<.009). This interaction remained significant when we covaried for the average thickness across the entire cortex in each individual (F1,61 = 5.25; P<.02). Follow-up analyses confirmed that patients showed significant thinning (relative to controls) in the temporal-prefrontal cortices (F1,61 = 4.6; P<.03), but not in control regions(F1,61 = 2.5; P = .12). There were no significant interactions between data set and group.
Figure 3 (rows A and B) shows statistical difference maps of data sets 1 and 2 at an uncorrected threshold of P<.05. The overall pattern of thinning across the cortical surface was strikingly similar in both data sets, although they constituted independently recruited patients and controls who underwent scanning on different systems. This provided some internal validation of our main finding. Thinning was most apparent in the prefrontal cortex (particularly in the inferior frontal, orbitofrontal, and medial prefrontal cortices) and in the temporal cortex (superior, middle, and inferior temporal cortices, fusiform and parahippocampal gyri, and temporal poles). Thinning was also observed in the left occipitotemporal junction and the left insula. In both data sets, the primary visual cortex, central sulcus, precentral and postcentral gyri, and superior parietal cortex appeared to be relatively (although not entirely) spared. In data set 1, a small region of relative thickening was found in the patients around the central sulcus that is probably an artifact arising from suboptimal gray-white contrast in this region. Artifactual thickening was not observed in data set 2, where the acquisition sequence yielded better gray-white contrast in this region.
To maximize statistical power and to determine patterns of cortical thinning that were most robust across all patients and controls, we combined the 2 data sets to construct an overall statistical thickness difference map, covarying for the effect of data set (Figure 3C). Identical numbers of patients and controls were included in each data set, making it unlikely that their combination would systematically confound findings. This was confirmed by an analysis of variance that showed a significant main effect of group in global cortical thickness (ie, thickness across the entire cortical surface, collapsing across both hemispheres and data sets) (patients, 2.26 mm; controls, 2.32 mm; F1,60 = 5.30; P = .03), but no scanner × group interaction (F1,60 = 0.005; P = .95). Moreover, there were no significant interactions between group and data set in the ROI analyses described herein.
As shown in Figure 3C, regions of statistically thinner cortex in patients relative to controls (at P<.001 uncorrected; equivalent to P<.05 corrected, shown in yellow) included the orbitofrontal cortices bilaterally (approximate Brodmann areas [BAs] 10 and 11), the left inferior frontal gyrus (approximate BAs 44 and 45), the right medial frontal cortex (approximate BAs 32 and 24), the left occipitotemporal junction (approximate BAs 22, 21, and 37), the right medial temporal cortex (approximate BAs 28, 27, 34, and 35), and the left inferior temporal gyrus (approximate BA 20). All these regions except for the right orbitofrontal and right medial frontal cortices showed significant differences when the average thickness across the entire cortex in each individual was included as an additional covariate. At a subthreshold significance level (P<.05, uncorrected, shown in red) thinning in the patients was also observed in the right inferior frontal cortex (BAs 44, 45, and 47), the dorsolateral prefrontal cortex bilaterally (approximate BAs 8, 9, and 46), the left anterior medial frontal cortex (approximate BAs 32 and 24), the right inferior temporal cortex (approximate BA 20), the left medial temporal lobe (approximate BAs 27, 28, 34, and 35), the temporal poles bilaterally (approximate BAs 20 and 38), the fusiform gyri bilaterally (BA 37), and the left insula. These regions did not show significant differences in thickness when we covaried for the average thickness across the entire cortex in each individual. Even at P<.05 uncorrected, the superior parietal cortex and primary visual, somatosensory, and motor cortices appeared to be relatively spared. Moreover, although our tests were 2-tailed, the amount of thinning was overwhelmingly greater than thickening in the patients compared with the controls, ie, even at P<.05 uncorrected, virtually no false-positive findings were detected for thickening.
Figure 3D shows a map with mean differences in cortical thickness at uniformly sampled loci across the entire cortex, confirming an overall pattern of temporal-prefrontal thinning in the patients relative to controls.
In the patient group, age was significantly inversely correlated with cortical thickness (ie, thinner cortex with older age) in all temporal-prefrontal ROIs (all Pearson correlations, P<.004) except for the parahippocampal gyri (P = .68). In the control group, age showed significant inverse correlations with cortical thickness in all temporal-prefrontal ROIs (P<.04) except in the anterior cingulate cortices (P = .51) and the parahippocampal gyri (P = .58). There were no significant age × group interactions in any of the ROIs, ie, there was no evidence of different age-related effects between the 2 groups (P>.05 for all).
Years of education did not correlate with cortical thickness in any of the temporal-prefrontal ROIs in the patient or the control group. There were no significant differences in cortical thickness between men and women within the patient or the control group (P>.20 for all). Within the patient group, there were no significant correlations between cortical thickness in any of the temporal-prefrontal ROIs and medication dosage as indexed by chlorpromazine equivalents (P>.83 for all) or severity of illness as assessed using the global Positive and Negative Syndrome Scale (PANSS) score (P>.21 for all). There were significant correlations between cortical thickness and duration of illness in the inferior frontal cortex (P<.05), the dorsolateral prefrontal cortex (P<.02), and the anterior cingulate gyrus (P<.01). These correlations, however, did not remain significant when age was covaried as a potentially confounding variable (P>.35 for all).
We mapped the more focal regions that showed significant (P<.001) differences in the combined statistical thickness difference map (Figure 3C) back on each individual subject's cortical surface and calculated the average thickness within each of these focal regions in each individual. Within both groups, focal thickness within the left inferior frontal gyrus correlated significantly (P<.05) with the overall average thickness across all temporal-prefrontal ROIs. In the control group, but not the patient group, the other focal regions that correlated significantly with overall temporal-prefrontal thickness were the left occipitotemporal junction (P<.03), the left inferior temporal cortex (P<.03), and the right medial temporal cortex (P<.01).
When we covaried for overall temporal-prefrontal thickness in each individual, there remained significant (P<.01) differences between patients and controls in the thickness of all focal regions.
Using MR methods, we have demonstrated, for the first time to our knowledge, regionally selective thinning of the cortical ribbon in schizophrenia. The ROI analysis showed that cortical thinning was significantly greater in patients with schizophrenia (relative to matched controls) in temporal-prefrontal ROIs than in selected control ROIs (the superior parietal cortex and primary visual, somatosensory, and motor cortices). This was the case even when we covaried for global thickness across the entire cortex in each individual subject.
Both the mean thickness difference maps and subthreshold statistical thickness difference maps (for data sets 1 and 2 individually) confirmed an overall distribution of cortical thinning that primarily affected prefrontal and temporal cortices while relatively sparing the superior parietal cortex, primary visual regions, the central sulcus, or the precentral and postcentral gyri. The statistical thickness difference map that combined both data sets showed significant thinning within several more focal regions within the left lateral and inferior temporal-prefrontal cortex (left inferior frontal, orbitofrontal, inferior temporal, and occipitotemporal cortices) and within the right medial temporal-prefrontal cortex (right medial temporal, orbitofrontal, and medial frontal cortices). In comparing the ROI and brainmap findings more quantitatively, cortical thinning averaged across all temporal-prefrontal ROIs in patients relative to controls did not entirely explain thinning in these more focal regions, suggesting that thinning was nonhomogeneous within the temporal-prefrontal cortex.
We attempted to exclude several methodological factors as potential confounds of our main finding that cortical thinning in patients (relative to controls) was greater in temporal-prefrontal than control regions. First, we visually inspected the entire cortex in each individual subject, slice by slice, manually fixing any inaccuracies in segmentation that might have arisen from image artifact. Second, to exclude the possibility that our main finding arose as a result of systematic differences in the registration of the main cortical folds in the schizophrenic and control groups, we checked the accuracy of registration of the major gyri/sulci binding ROIs in each patient and control participant (Figure 2). Third, to exclude the possibility that our main finding was driven by differences between temporal-prefrontal and control regions in the between-subject variability of cortical thickness, we constructed mean thickness difference maps (that do not take into account variance between subjects). The overall pattern of thinning in these mean thickness difference maps (Figure 3D) was qualitatively similar to that seen in the statistical thickness maps (Figure 3C).
Volume is a product of cortical thickness and surface area. Volumetric deficits have been previously described in all of the regions in which we report cortical thinning in the present study,1,2,8- 10 suggesting that cortical thinning may contribute to volumetric abnormalities in schizophrenia. Reduced surface area has also been reported in some of the regions in which we demonstrated cortical thinning, such as the orbitofrontal cortex.33 Surface area and thickness are known to result from distinct stages during cortical development.34 If, as hypothesized, schizophrenia arises from abnormal neurodevelopment,35- 37 understanding the relationships between volume, thickness, and surface area will shed further light on its pathogenesis.
Because the methods we used probed gray matter thickness rather than volume, they do not index gray matter deficits that have been previously documented in subcortical structures such as the hippocampus38 and the thalamus.39 However, both of these subcortical regions have dense reciprocal projections with the cortical regions in which we observed thinning. Our findings are therefore consistent with models that have boldasized the role of hippocampocortical and thalamocortical circuitry in the pathoetiology of schizophrenia.38- 40
Abnormalities in cortical thickness in temporal and prefrontal regions may reflect underlying histopathological changes in these regions. Indeed, a recent study that used magnetization transfer imaging techniques suggested subtle or early neuropathological abnormalities in many of these temporal-prefrontal regions in patients with schizophrenia.41 What specific neuropathological abnormalities might underlie the cortical thinning we observed in the current study? Although there is evidence of volume loss and a decrease in the total number of cells in the hippocampus and thalamus, the precise nature of neuropathological abnormalities in the cerebral cortex in schizophrenia remains more controversial.14- 16 Almost all postmortem studies report only changes in cell density and do not report a reference volume, making it difficult to determine whether there is a loss in the total number of cells in the cortex.14 Some studies have reported cell-specific or layer-specific decreases in nerve cell density, particularly in interneurons in the prefrontal cortex.42 Others, however, report increases in cell density43,44 and, by inference, a decrease in the amount of interneuronal neuropil, in regions such as the dorsolateral prefrontal cortex.45 The precise relationships between these findings and overall cortical thickness remain undetermined. Combining the techniques presented herein with postmortem analysis, by targeting regions with maximal cortical thinning, may help clarify the neuropathological basis of morphological changes in schizophrenia.
The pattern of gray matter cortical deficits described in the present study can be linked with neuropsychological deficits in schizophrenia in a variety of cognitive domains.11 Temporal and inferior frontal atrophy may lead to abnormalites in semantic and language function; gray matter deficits in the dorsolateral prefrontal cortex may be linked with working memory and executive dysfunction; medial temporal lobe deficits may contribute to declarative memory dysfunction; and atrophy in the orbitofrontal, cingulate, and insula cortices may underlie abnormalities in emotional and social cognition.11 Abnormal patterns of activation in all of the regions that showed cortical thinning in this study have been described in functional neuroimaging studies in schizophrenia during cognitive activation paradigms of language and semantic processing,46 declarative memory,47 working memory,48,49 emotional processing,50 and social cognition.51 We are currently using the same cortical reconstruction and surface-based averaging techniques described herein in our analysis of functional MRI data. This will help discriminate between an anatomical vs a functional basis of activation abnormalities in schizophrenia, and allow us to address the fundamental question of how brain structure and function are linked in schizophrenia.
We have demonstrated selective thinning of the cerebral cortex within the temporal and prefrontal cortices in patients with schizophrenia. There are important considerations in the interpretation of our data. First, although we showed no significant correlations between cortical thickness in ROIs and overall disease severity as indexed using the global PANSS score, schizophrenia is a clinically heterogeneous disorder and the PANSS includes a wide variety of positive and negative symptoms. The findings reported herein could very well have arisen as expressions of effects in only parts of the clinical group. Evidence from volumetric studies suggests that some psychopathological variables, such as positive thought disorder,52 auditory hallucinations,53 negative symptoms,9 and global function,54 are correlated with morphometric changes in particular regions. We hope to determine whether these variables are also correlated with regional differences in cortical thinning in future studies in larger samples with a wider range of clinical symptoms and disease severity. It will also be interesting to determine whether regional cortical thinning is more marked in patients with a positive family history of schizophrenia.55
Second, both of the samples described herein included patients receiving long-term medication therapy. Volumetric reductions in some of the same regions in which we report thinning have been described in schizophrenic patients in their first episode of illness.56 Future studies using this technique will determine whether first-episode schizophrenic patients also show cortical thinning. In the present study, medication dosage and illness duration (covarying for age) did not predict cortical thinning in ROIs. However, we do not wish to overinterpret these nonsignificant correlations. Previous studies have suggested that patients may show a greater rate of decline with age in gray matter volume than controls.57,58 Future studies in a larger sample of patients will determine whether this is also true of cortical thickness changes. Indeed, the methods presented herein have potential for longitudinal studies of cortical thickness in schizophrenia. They are automated and objective and can be used to detect subtle changes in thickness over time, even in individual subjects. This will allow us to determine changes during the course of the schizophrenic illness and the effects of different pharmacological treatments.
Corresponding author and reprints: Gina R. Kuperberg, MD, PhD, MRCPsych, Department of Psychiatry, Massachusetts General Hospital (East), CNY-9, Bldg 149, Charlestown, MA 02129 (e-mail: firstname.lastname@example.org).
Submitted for publication November 1, 2002; final revision received January 27, 2003; accepted January 30, 2003.
Dr Kuperberg was supported by the Wellcome Trust, London, England, and is currently supported by grant K23 PA-00-004 from the National Institute of Mental Health, Rockville, Md. This work was also supported by the Institute for Mental Illness and Neuroscience Discovery, Albuquerque, NM. Dr Ozawa was supported by Hamamatsu Photonics, Shizuoka, Japan.
This research used a computer program, FreeSurfer, for data analysis. FreeSurfer, distributed freely for research purposes, is co-owned by an entity in which Dr Dale has a financial interest.
We thank Stephan Heckers, MD, for helpful discussion on the neuropathology of schizophrenia in relation to the current findings and Bruce Rosen, MD, PhD, for helpful discussion.