INTRODUCTION

Investigations of the structure of the human brain, both via in vivo imaging (Shenton et al, 2001; Niznikiewicz et al, 2003) and post-mortem studies (Harrison, 1999; Harrison and Lewis, 2003), are essential for understanding the neurobiology of schizophrenia. However, use of antipsychotic medications as treatment for schizophrenia represents a potential confounding factor in many of these studies. For example, MRI studies of structural brain changes across time are limited by the fact that—due to obvious ethical reasons—neither untreated subjects with schizophrenia nor control subjects exposed to antipsychotic medications can be used as comparison groups. Similarly, in post-mortem studies almost all subjects with schizophrenia have been treated with antipsychotic medications, whereas normal comparison subjects were not exposed to such medications.

We have advocated a three-fold strategy to address the problem of confounding due to chronic antipsychotic drug treatment in post-mortem studies of schizophrenia (Lewis, 2002). This strategy includes (1) the comparison of subjects with schizophrenia who were on or off medication at the time of death, (2) the comparison of subjects with other psychiatric disorders (eg major depression with psychosis) who were on or off antipsychotic medication at the time of death, and (3) the use of non-human primates exposed chronically to antipsychotic medications in a manner that mimics clinical practice. Although each of these approaches has its own intrinsic limitations, a convergence of findings across studies provides compelling evidence that antipsychotics do or do not account for the finding of interest. The third arm of this triad is particularly important because it is the only approach that permits fully controlled experimental conditions. In this regard, the non-human primate model may be informative not only for human post-mortem studies but also for interpreting the outcome of longitudinal in vivo imaging studies of subjects with schizophrenia. For example, such studies have found evidence of brain volume changes evolving with time (Thompson et al, 2001; Ho et al, 2003; Pantelis et al, 2003; Sporn et al, 2003) that are potentially related to disease progression, but that might also reflect an effect of treatment.

In addition to their effects on the brain in vivo (eg drug-induced changes in cell number in a given brain region), antipsychotic drugs might also induce changes in brain structure or composition that affect post-mortem tissue behavior. For example, knowing whether antipsychotic drug exposure in vivo causes differential tissue shrinkage during histological processing is essential for interpreting post-mortem volumetric studies.

Consequently, in this study, we developed a non-human primate model of oral antipsychotic treatment by chronically exposing groups of macaque monkeys to either haloperidol (a conventional antipsychotic), olanzapine (a second-generation antipsychotic) or sham (placebo). Using this model, we addressed the following questions: (1) Do antipsychotic medications induce changes in total or regional brain weight or volume? (2) Do antipsychotic medications affect the tissue shrinkage that occurs during post-mortem processing? (3) Do these effects differ as a function of type of antipsychotic drug?

METHODS

Animals and Training

Eighteen experimentally naïve male macaque monkeys (Macaca fascicularis), 4.5–5.3 years of age, were obtained from LABS of Virginia (Yemassee, SC). All animals were housed in single cages grouped together in the same room at the University of Pittsburgh Primate Research Laboratory. Six animals were randomly selected to form a group, with the remaining animals divided into two groups so that average body weights were approximately balanced across the three exposure groups (haloperidol, olanzapine, or sham). All studies were carried out in accordance with the NIH Guide for the Care and Use of Laboratory Animals and were approved by the University of Pittsburgh Institutional Animal Care and Use Committee.

All animals were trained to self-administer fruit punch-flavored sucrose pellets. Self-administration of the pellets was reinforced with raisins and was followed by drinking 60 ml of orange drink from a syringe. The orange drink was administrated in order to have an alternate vehicle for drug administration, should an animal become unreliable in ingesting the pellets. When animals achieved 100% compliance, sucrose pellets containing drug were introduced to the haloperidol group. Owing to delays in obtaining olanzapine, the animals in this group received sham pellets for the first 10 months. Two monkeys initially allocated to the sham group were switched with two initially allocated to the olanzapine group immediately before the initiation of the olanzapine treatment (see Animal Health section below).

Drug Preparation and Administration

Sucrose pellets (190 mg) containing haloperidol sulfate (Spectrum Chemicals and Laboratory Products, Gardena, CA), olanzapine (Eli Lilly, Indianapolis, IN), or no drug (sham pellets) were custom made (Research Diets Inc., New Brunswick, NJ). Two dosage levels per pellet were ordered for each drug. Quality control assays determined the content of haloperidol to be 1.0 or 2.0 mg/pellet and the content of olanzapine to be 0.55 or 1.1 mg/pellet. Methods for quality control were identical to those used for drug plasma level analysis (see Blood Sampling section below).

All medication and sham treatments were administered orally twice daily (between 08:30 and 09:30, and between 15:30–16:30), 5 days a week for the first 5 months and all 7 days a week for the rest of the study (22 months). Each monkey could receive the drug in one of three ways: (1) by ingesting intact pellets; (2) by ingesting the dissolved drug (haloperidol) or pulverized pellets in orange drink; and (3) by ingesting pellets pulverized in treats consisting of fresh or dehydrated fruit, or mixed with peanut butter or jelly. All monkeys were always given orange drink containing quinine sulfate as a taste mask to prevent the animals from knowing whether or not drug was present. During the study, animals in the haloperidol and the olanzapine groups required increasingly greater effort to maintain compliance with the dosing regimen. The animals seemed to develop an aversion to the taste and/or the subjective effects of the medication. Thus, dosing vehicles were adapted to ensure that the animals complied with drug ingestion.

Drug doses were increased over time based upon the monitored drug plasma levels. The goal was to achieve and maintain plasma concentrations of each drug that approximated those associated with a therapeutic response in humans. In a previous study (Akil et al, 1999; Pierri et al, 1999), we exposed monkeys to doses of haloperidol decanoate, resulting in mean (SD) trough haloperidol plasma levels of 4.3±1.1 ng/ml, which are similar to the mean (SD) plasma level of 4.9±2.9 ng/ml (total+reduced haloperidol) observed in humans treated according to the model of ‘neuroleptic threshold dosing’ (McEvoy et al, 1986, 1991). Consistent with this strategy, we observed notable extrapyramidal side effects in these animals. However, more recent imaging studies of dopamine D2 receptor occupancy suggest that antipsychotic efficacy, without pronounced extrapyramidal symptoms or hyperprolactinemia, occurs when 65–72% of striatal dopamine D2 receptor are occupied (Farde et al, 1992; Nordström et al, 1993; Kapur et al, 1998, 2000; Kapur and Seeman, 2001). In humans, this level of occupancy is associated with haloperidol plasma levels of 1–1.5 ng/ml (Nyberg et al, 1995; Kapur et al, 1997) or olanzapine plasma levels of 10–25 ng/ml (Kapur et al, 1998, 1999). Furthermore, clinical trials have documented that antipsychotic dosages corresponding to these plasma levels are associated with the desired antipsychotic effect (Perry et al, 1997, 2001; Apiquian et al, 2003; Keefe et al, 2004; Oosthuizen et al, 2004). Consequently, in this study, we sought to obtain plasma levels of 1–1.5 ng/ml for haloperidol and of 10–25 ng/ml for olanzapine. We chose to use oral administration of both drugs in order to mimic more closely the most common mode of drug administration in humans and because a depot form of olanzapine was not available.

Blood Sampling

Prior to the study, all animals were trained to enter a transfer box and then a restraining chute so that blood samples could be obtained without sedation, thus minimizing the exposure of animals to ketamine, an NMDA-receptor antagonist that in humans can induce clinical features reminiscent of schizophrenia. All samples were drawn from the femoral vein approximately 7–8 h after the morning dose in order to obtain trough plasma levels. Assays for haloperidol and olanzapine were performed using high-performance liquid chromatography, with the haloperidol assay based upon electrochemical detection (Gupta et al, 1995) and the olanzapine assay based upon ultraviolet detection (Olesen and Linnet, 1998). Interassay coefficients of variation, as estimated from the results of quality control samples, were 8.9% for haloperidol and 5.3% for olanzapine. The lower limit of quantitation for haloperidol was 0.25 ng/ml and the corresponding value for olanzapine was 2 ng/ml.

Animal Health and Behavior

No quantitative measures of behavior were obtained during the study, other than those associated with compliance with the drug exposure regimen. Qualitatively, animals in both the haloperidol and olanzapine groups showed signs of mild bradykinesia and sedation that appeared to worsen with increasing dose, and some animals seemed to become more aggressive with time. This aggressive behavior was noted to develop in four of the olanzapine-exposed animals and in two of the haloperidol-exposed animals.

Animals had free access to water, and were fed in the evening after receiving the afternoon dose of drugs to increase the likelihood that they would consume drug-containing food or drink. All animals were given between 20 and 30 monkey food pellets (Purina Monkey Chow) per day, which was standard for the population of macaque monkeys at the University of Pittsburgh Primate Research Laboratory. Throughout the study, all animals were provided access to other forms of enrichment including seeds for foraging, toys, and fruit and other food treats. Thus, although the animals were not given unlimited access to food, they were not food restricted.

Throughout the study, all animals were regularly monitored by the University of Pittsburgh Veterinary Staff. Over more than 2 years of study, only two animals—both in the olanzapine group—exhibited medical complications: One animal (2846) experienced recurrent episodes of diarrhea, which resolved with bowel rest or metronidazole, and which did not interfere with the animal's drug compliance. The other animal (2836) engaged in self-injurious behavior prior to starting drug treatment. This animal recurrently cut his hands and feet with the canines, and on several occasions required stitches to close wounds. During the study, this monkey also received NSAIDs and pain medication (buprenorphine) on several occasions. The monkey was initially in the sham group, but was switched into the olanzapine group (resulting in a shuffle of four monkeys total in order to keep the means of the body weights matched across groups), immediately before the initiation of olanzapine treatment, on the rationale that the olanzapine might decrease the frequency of self-injurious behavior. However, no olanzapine-related change in this behavior was seen. To address these behavioral problems, the animal was provided with increased human contact, doubling of cage space, and enhanced enrichment, and the cage was moved to an area of decreased stimulation. Under these circumstances, the animal's behavior stabilized. It should be noted that excluding the data from this animal did only change the results marginally and did not alter any conclusion.

Brain Removal and Measurements

Matched by terminal body weight, monkeys were euthanized in triads—one from each group. All three monkeys of each triad were euthanized on the same day, with the six triads euthanized over a period of 74 days. Each monkey was euthanized by an overdose of pentobarbital (50 mg/kg) and removal of the fresh brain was completed within approximately 25 min. The weight of the fresh, whole brain was obtained on a calibrated lab scale (PJ3600 DeltaRange, Mettler-Toledo, OH) with a precision of 0.01 g. The same scale was used for all subsequent measures (described in detail below). The cerebellum and brainstem were then separated from the cerebrum by a cut through the mesencephalon at the level of the superior colliculi. The cerebrum was split by a sagittal cut through the corpus callosum and the cerebellum and brainstem were similarly bisected. The right half of each brain was cut into slabs (5 mm) and frozen fresh in isopentane on dry ice and subsequently stored at −80°C. The left cerebellum/brainstem and the cerebrum were immediately placed in ice-cold 0.1 M phosphate buffer, pH=7.3.

After 5 min in the buffer, the left cerebrum was dissected into four pieces following the gyration pattern: (1) The frontal piece was removed by cutting (locally perpendicular to the pial surface) through the superior frontal gyrus and the anterior cingulate gyrus at the level of the central sulcus, within the corpus callosum along its dorsal and rostral edge, at the bottom of the central sulcus, along the rostral part of the circular sulcus of the insula, and around the temporal pole. During this procedure, the striatum was split. (2) The occipital piece was removed by a planar cut through the bottom of the lunate sulcus and extending to the medial surface of the hemisphere. (3) The parietal piece was separated from the remainder of the brain by cutting through the corpus callosum along its dorsal edge, at the bottom of the lateral fissure, and through the angular gyrus. (4) The remaining piece consisted of the temporal lobe, insula, and most of the basal ganglia, including the entire thalamus.

When the dissection was complete, the five pieces from the left hemisphere were individually blotted dry and weighed, and their volumes determined by Archimedes' principle of fluid displacement. We used a high-precision approach (Figure 1), in which the weight of the displaced fluid is measured (Scherle, 1970; Weibel, 1979). The pieces were then placed in ice-cold 4% paraformaldehyde in 0.1 M phosphate buffer (pH 7.4) for 48 h, washed in a graded series of sucrose solutions (12, 16, and 18%), and stored in an antifreeze solution (30 vol% ethylene glycol and 30 vol% glycerin in 0.024 M phosphate buffer) at −30°C. The post-mortem interval (PMI, time from death to immersion into fixative) for the left hemisphere pieces was kept at approximately 72 min for all 18 brains.

Figure 1
figure 1

Measurement of tissue weight and volume. Left: The tissue piece is weighed at the scale. Right: The principle of Archimedes ad modum Scherle (1970). A beaker containing a fluid is placed on the scale, the tissue supporter is inserted into the fluid to a mark on the handle, and the scale is reset to zero, preventing the volume of the supporter from biasing the measure. Then, the tissue piece on the supporter is lowered into the fluid to the mark, and the readout of the scale indicates the weight of the displaced fluid. The volume of the displaced fluid (ie the volume of the tissue piece) is then calculated using the density of the fluid. W and V are weight and volume of the tissue piece, and Wd and ρd are weight and density of the displaced fluid (storing solution in this example).

To assess the tissue deformation (predominantly shrinkage) introduced by these tissue-processing procedures, the weight and volume of each piece were determined, using the approach described above, at every major step in the tissue processing (ie fresh, after 48 h of fixation, after 12 days in sucrose solutions, after 11 days in storing solution at 4°C, and finally after approximately 1 year in storing solution at −30°C). Storing solution was used for the fluid displacement in the final measure, while all the previous measures were made in 0.1 M phosphate buffer.

Stereological Assessment of Parietal Piece Volumes

The volumes of gray and white matter in the parietal pieces were estimated using Cavalieri's principle (Gundersen and Jensen, 1987; Howard and Reed, 1998). The parietal pieces (equilibrated in storing solution) were embedded in agar—7% low-melt agarose (SeaPlaque Agarose, Cambrex, Rockland, ME) in storing solution—and cut into 2.5-mm thick slabs perpendicular to the intraparietal sulcus. The slicing procedure ensured that all cuts had a systematic, uniformly random position within each tissue piece. The slabs were returned to storing solution and coded as to subject number and group. The resulting cut surfaces of the slabs (12–15 slabs per brain) were observed using a Nikon SMZ1500 Stereoscopic Zoom Microscope, mounted on a Nikon C-DSDF diascopic stand F, illuminated from above by fiber optic light from a 150 W Tungsten Halogen Light Source (Schott-Fostec, Auburn, NY), and equipped with a 0.5 × HR Plan Apo Nikon objective, a Nikon DN100 Digital Net Camera (1280 × 960 pixels), and a Triniton Multiscan G420 color monitor. Each cut generates two cut surfaces, which obviously are mirror images of each other. Owing to the overall spindle shape of the parietal pieces, most slabs had a conical appearance where one cut surface was bigger than the other. It turned out to be very difficult to tell the exact outer boundary of the smaller cut surface due to overprojection of the larger cut surface through the transparent agar matrix. Therefore, the caudal surfaces of the rostral slabs were observed until the maximum total area was reached; the slab with maximal area was ignored; and the rostral surfaces of the caudal slabs were observed. In this way, each cut was observed using the cut surfaces that were the larger on its slab and therefore not biased by overprojection. All slabs were observed by one investigator (blind to exposure status) at a final magnification of 9 × at the monitor. A square point grid printed on a transparency was superimposed uniformly random over the monitor image of each cut surface. The areas of gray and white matter observed at each cut surface were estimated from the number of points hitting the gray or white matter, respectively. From these counts, the gray and white matter volumes for each parietal piece were estimated as:

where T=2.5 mm is the mean slab thickness, a=2.29 mm2 is the area per point of the point grid (corrected for magnification), Ai is the gray or white matter area of the ith cut surface, and Pi is the respective number of points hitting gray or white matter. For each parietal piece, on average a total of 430 points were counted over gray matter and 215 over white matter. The corresponding mean coefficients of error, CE, were estimated to be 0.02 for the gray matter volumes, 0.03 for the white matter volumes, and 0.01 for the combined total volume using the methods of Gundersen et al (1987, 1999).

Precision and Accuracy of the Archimedean Approach

In order to assess the precision of the weight and Archimedean measurements, we performed two sets of measurements for the 15 pieces from one triad of monkeys. The measurements were obtained by the same person, morning and afternoon of the same day, using samples from the same batch of storing solution for displacement. However, the tissue pieces were returned to the freezer, the scale recalibrated, and the glassware changed between the measurements. The concordance of the remeasured weights and volumes (Figure 2) indicated a very high level of precision corresponding to a CE≈0.004 for both weights and volumes.

Figure 2
figure 2

Precision of the Archimedean volume measures as demonstrated by 15 pairs of observations. Filled circles indicate sham; thin open circles, haloperidol; and thick open circles, olanzapine. Line is the identity line.

In order to assess the accuracy of the Archimedean measurements, we measured the volume of all 18 parietal pieces immediately before embedding in agar. In this way, we could compare the total volume measures obtained by the Archimedes' principle to corresponding volume estimates obtained by Cavalieri's principle. To measure the inter-rater reliability of the Archimedean measurements, two observers measured the same 15 tissue pieces at the same time, but blind to each other. The intraclass correlation coefficient of 0.9998 (95% confidence interval, 0.9993–0.9999) illustrated that our Archimedean measurements are as reliable as weighing the tissue. The weight and the Archimedean measures were performed by several observers who were not formally blind to the exposure status. Owing to the very high inter-rater reliability, precision, and accuracy of the methods, we judged there were virtually no room for any observer bias and no need for blinding at this initial stage.

Statistics

To examine the difference in initial body weight and age of the animals at the time of euthanasia among exposure groups (sham, haloperidol, and olanzapine), one-way analysis of variance (ANOVA) models were used, with the exposure group as the main effect. To explore the differences in terminal body weight among the three groups, an analysis of covariance (ANCOVA) model was used, with group as the main effect and initial body weight as a covariate.

One-way ANOVA models with the exposure group as the main effect were performed to examine the differences among the three groups for each primary response variable. For the multiple comparisons procedure, which compared both drug exposed groups with the sham group, we used Dunnett's procedure (Hochberg and Tamhane, 1987).

To examine the robustness of the ANOVA conclusions, the Kruskal–Wallis test, a nonparametric alternative to the one-way ANOVA, was carried out to test the group effect for each response variable. As a nonparametric analogue of the Dunnett's multiple comparisons procedure in one-way ANOVA, Steel's procedure (Hochberg and Tamhane, 1987) was used to confirm results of the comparisons of the drug-exposed groups with the sham group.

A two-way ANOVA model was employed to analyze the main group effect and the triad blocking effect on the change in the left brain volume as a result of tissue processing. Multiple comparisons among triads used Tukey's procedure (Neter et al, 1996). To examine the difference in shrinkage among the five regions of the brain, the multivariate analysis of variance (MANOVA) model based on repeated measures with a compound symmetric covariance structure (Littell et al, 1996) was used with exposure group as the main effect, triad as a blocking effect, and region as a within-subject effect. Regions were compared using the Bonferroni multiple comparisons procedure.

Analyses were implemented in SAS PROC GLM (SAS Institute Inc., 1990), SAS PROC MIXED (Littell et al, 1996), and SAS PROC NPAR1WAY (SAS Institute Inc., 1990). F-tests for the effect of diagnostic group are based on type III sum of squares. An alpha level of 0.05 was selected as indicator of significance for all tests and multiple comparisons.

RESULTS

Monkey Model of Chronic Antipsychotic Drug Exposure

The doses of haloperidol and olanzapine were gradually increased (Figure 3, top), until the monkeys were receiving 24–28 mg of haloperidol or 11.0–13.2 mg of olanzapine per day. The final trough drug plasma levels were 1.5 ng/ml for haloperidol and 15 ng/ml for olanzapine (Figure 3, bottom). Neither the initial (by design) nor the terminal mean body weight of the monkeys differed significantly across the groups. Mean (coefficient of variation CV=SD/mean) initial body weights were 4.40 (0.21) kg, 4.12 (0.14) kg, and 4.12 (0.29) kg for the sham, haloperidol, and olanzapine groups, respectively. The corresponding mean (CV) terminal body weights were 6.05 (0.20) kg, 6.13 (0.20) kg, and 6.05 (0.12) kg. Most animals did increase substantially in weight during the experimental period (Figure 4), consistent with the expected growth of male animals of this species and age. The mean (CV) age of the animals at time of euthanasia did not differ (F2,15=0.48, p=0.63) and were 84.0 (0.03) months, 85.2 (0.05) months, and 83.1 (0.04) months for the sham, haloperidol, and olanzapine groups, respectively.

Figure 3
figure 3

Top: Administrated daily drug dose as a function of time for each monkey in the haloperidol- (left) and olanzapine- (right) exposed groups. Numbers within each plot indicate how many monkeys received that drug dose. For example, at day 651, two monkeys were raised to a final haloperidol dosage of 28 mg/day, and 23 days later, additional three animals were raised to that level. The punctuated line and the light gray area at the beginning of the haloperidol graph indicate the initial period with drug administration limited to 5 days a week. Bottom: Mean trough plasma levels of the respective drug for the six monkeys in each group. The error bars indicate the standard error of the mean for each group at each time point. The plots at the bottom terminate earlier than the plots at the top because the monkeys were euthanized (in triads) over a 74-day period.

Figure 4
figure 4

Increase in mean body weights for each group (S, sham; H, haloperidol; O, olanzapine) across the course of the study. Vertical lines indicate the time points when each of the two drugs was introduced. The mean (CV) percentage gain in weight for each group was: sham, 38% (0.37); haloperidol, 53% (0.79); olanzapine, 53% (0.53).

Fresh Brain Weight and Volume

The mean fresh brain weights differed significantly among the three exposure groups (F2,15=3.99, p=0.041), with those of the haloperidol (H) and olanzapine (O) groups 8.1% and 9.6%, respectively, lower than that of the sham (S) group (Figure 5). The fresh weight and volume of the left cerebrum were calculated for each brain as the sum of the respective measures of the four cerebral tissue pieces, as this was the most anatomically well-defined unit in the study. Mean fresh left cerebral weights differed significantly (F2,15=4.85, p=0.024) among the groups; weights of the H and O groups were 8.9% and 10.9% lower, respectively, than the S group (Table 1). The corresponding mean fresh left cerebral volumes showed virtually the same pattern: H vs S decreased by 8.8% and O vs S decreased by 10.5% (F2,15=4.70, p=0.026).

Figure 5
figure 5

Fresh brain weights for the sham- (S), haloperidol- (H) and olanzapine- (O) exposed monkeys. The horizontal bars indicate group means.

Table 1 Mean (CV) Fresh Weight and Volume of the Left Cerebrum for Each Group of Animals

Comparison of mean fresh weights of individual tissue pieces suggested that the differences were more robust in the frontal and parietal regions (Figure 6). Specifically, the mean weight of the frontal piece was significantly (F2,15=4.55, p=0.029) lower by 10.1% and 10.4% in the H and O groups, respectively, relative to S; for the parietal piece, H was decreased by 11.2% and O was decreased by 13.6% (F2,15=5.20, p=0.019) relative to S. The other three regions followed the same general pattern, but the differences did not achieve statistical significance: temporal, H vs S decreased by 6.9% and O vs S decreased by 10.5% (F2,15=2.31, p=0.13); occipital, H vs S decreased by 9.7% and O vs S decreased by 9.8% (F2,15=0.60, p=0.56); and cerebellum and brainstem, H vs S decreased by 8.7% and O vs S decreased by 8.5% (F2,15=0.92, p=0.42). The pattern for the corresponding volumes was virtually identical, indicating only minor variability in the tissue density within the same brain region across the monkeys. Results of the multiple comparisons procedure are shown for all significant ANOVAs in Figure 7.

Figure 6
figure 6

Fresh weights of the individual tissue pieces from the left brain of the monkeys of the sham- (filled circles), haloperidol- (thin open circles) and olanzapine- (thick open circles) exposed groups. The horizontal bars indicate group means.

Figure 7
figure 7

The 95% confidence intervals for the relative difference of each exposure group (haloperidol, H; olanzapine, O) vs sham group (S) calculated using Dunnett's method for each variable that reached significance in the main ANOVAs. For those intervals that cross the origin, the difference between groups did not reach the 0.05 significance level.

Changes in Brain Weight and Volume during Tissue Processing

During the first four phases of tissue processing, brain weight and volume decreased in an identical manner across the three exposure groups (Figure 8). Over time in storing solution, brain weight and volume partially returned to initial values. The mean final tissue shrinkage due to tissue processing was 20.5% (CV=0.13) and was virtually identical (F2,10=1.02, p=0.39) across the three groups of animals (Table 2). However, we observed a significant difference in final shrinkage across the six triads (F5,10=13.31, p=0.0004), with the brains from the monkeys of the first triad shrinking significantly less than the five other triads, even though the experimental protocol did not differ across triads. Additionally, a marked difference in shrinkage pattern was observed across the five tissue regions (F4,68=71.31, p<0.0001). Pair-wise comparisons indicated significant differences in the amount of shrinkage across regions, with the greatest shrinkage occurring in the temporal–basal ganglia piece and the least in the occipital region (Table 2).

Figure 8
figure 8

Plots of the mean weight (left) and volume (right) of the entire left hemisphere at specific stages in tissue processing for the sham- (S), haloperidol- (H) and olanzapine- (O) exposed monkeys.

Table 2 Mean (CV) Final Shrinkage of Brain Regions and of the Whole Left Brain for Each Group

Gray and White Matter Volume of the Parietal Region

In order to determine whether the reductions in brain volume observed in the H and O animals reflected changes in gray and/or white matter volumes, we estimated the volumes of these compartments in the stored parietal pieces (Figure 9). The mean gray matter volume was significantly (F2,15=6.12, p=0.011) reduced in the H (11.8%) and O (15.2%) compared to S, while similar decreases in mean white matter volume (H, 13.3%; O, 12.7%) almost reached significance (F2,15=3.03, p=0.079).

Figure 9
figure 9

Gray and white matter volume in the parietal pieces determined by Cavalieri's principle for each of the sham- (S), haloperidol- (H) and olanzapine- (O) exposed monkeys. The horizontal bars indicate group means. Lettering at the top indicates the results of Dunnett's procedure; regions not sharing the same letter are significantly different at p<0.05.

Accuracy of the Archimedean Approach

In order to verify the accuracy of the Archimedean volume measurement, we assessed the volumes of the parietal pieces using two different methods (Archimedes' vs Cavalieri's principle). The comparison of Cavalieri estimates and Archimedean measures of total parietal volume indicated a high degree of concordance with a small systematic difference; the Archimedean measures were larger by an average of 2.9% (Figure 10), without a systematic difference in bias across the exposure groups.

Figure 10
figure 10

Comparison of the volume measures of the parietal pieces obtained by Archimedes' principle and Cavalieri's principle for each of the 18 animals. Filled circles indicate sham; thin open circles, haloperidol; and thick open circles, olanzapine. Line is the identity line.

DISCUSSION

Effects of Haloperidol and Olanzapine on Brain Weight and Volume

In the present study, we observed smaller brain weights and volumes in male macaque monkeys after 1.5–2.3 years of exposure to haloperidol or olanzapine at plasma drug levels comparable to those in treated humans. Exposure to either of the drugs was associated with a similar reduction in mean fresh brain weight, as well as mean fresh weight and volume of the left cerebrum, compared to matched, placebo-exposed animals. For both drugs, the magnitude of these effects was in the range of 8–11%. The reduction seemed to be global (ie including all brain regions), but was most robust in the frontal and parietal lobes. In addition, both gray and white matter volumes appeared to be reduced to a similar degree.

Several factors must be considered when interpreting the results. First, although macaque monkeys provide a more relevant brain model for human studies than do rodents, they metabolize antipsychotic drugs differently than humans. For example, approximately 55% of an oral dose of olanzapine is absorbed in monkeys vs at least 65% in humans; the terminal elimination half-life is 3 h in monkeys vs 27 h in humans; and due to different metabolic pathways, the plasma levels of the various metabolites of unknown bioactivity are quite likely very different (Kassahun et al, 1997; Mattiuz et al, 1997). Likewise, haloperidol is metabolized markedly faster by monkeys than humans (JM Perel and L Roizin, personal communication). Consistent with these differences in metabolism, the doses (in mg/kg) required in reaching comparable drug plasma levels were substantially higher in monkeys than those used in humans.

Second, it is not clear how an apparent drug effect on the brain of monkeys translates to humans with schizophrenia or other brain disorders. Specifically, we cannot conclude that normal human brains would show the same changes in weight and volume as monkeys, and if they do, whether the effect would be of the same, greater, or smaller magnitude in individuals with brain disorders. Finally, we cannot infer from this study whether such changes represent a beneficial, unwanted, or incidental effect.

To our knowledge, previous structural studies of primates exposed to antipsychotic medication have not assessed either total fresh or fixed volumes or weights (Akil et al, 1999; Pierri et al, 1999; Selemon et al, 1999), although all assume, based on cortical thickness measures, that no cortical volume loss occurred. However, a volume loss of 10% as seen here corresponds to a linear change of 3%, the detection of which requires more accurate and precise measurements than previously used. In this context, it should be mentioned that each of these previous studies based their structural analyses on density measures, which are known to be very sensitive to volume changes (Mendis-Handagama and Ewing, 1990). For example, both the increased 2-D densities of GAT-1-immunoreactive cartridges reported by Pierri et al (1999), and the increased glial 3-D densities reported by Selemon et al (1999) associated with antipsychotics, are very difficult to interpret unambiguously without measures of total cartridge or glia number and the corresponding reference volume. In light of our findings, the increased densities reported by these two studies could be a result of a reduction in the reference volume.

Potential Sources of the Observed Reductions in Brain Volume

What happened to the brain volume in our monkeys? Regions enriched in dopamine D2 receptors might be hypothesized to be preferentially affected. However, due to the tissue-blocking procedure (see next section), it is not possible from our data to determine whether regions enriched in dopamine D2 receptors were predominantly affected. As both the gray and white matter appeared to be reduced to a similar degree in the parietal piece, the most parsimonious explanation of the findings would be a change in a structural element common to both compartments, that is, involving neurons and their processes (axons and dendrites), glia cells, and/or the vascular bed. Owing to the limited number of non-human primate studies, the following discussion is to a large degree based upon what we know from human studies.

The observed differences could be due to a decrease in the number or size of one or several types of neurons. Modern stereological studies of total neuron number in a range of brain regions have not found evidence of a deficit in neuronal number in subjects with schizophrenia (most of whom had been treated with antipsychotics) (Heckers et al, 1991; Pakkenberg, 1993; Thune et al, 2001; Andersen and Pakkenberg, 2003). Even the seemingly robust findings of reduced neuronal number in the mediodorsal thalamic nucleus previously reported by several groups have recently been questioned (Cullen et al, 2003; Dorph-Petersen et al, 2004; Young et al, 2004). On the other hand, studies of various brain regions in (predominantly medically treated) subjects with schizophrenia have reported decreased somal cell volumes (Benes et al, 1991; Arnold et al, 1995; Rajkowska et al, 1998; Pierri et al, 2001; Sweet et al, 2003, 2004), suggesting that a treatment-related reduction in neuronal size (with associated reductions in dendritic and axonal arbors) may be present. However, although consistent, these neuronal somal size studies should be interpreted with some care as none of these studies fully applied unbiased stereological methods of sampling and measurement. In addition, a recent stereological study of brains from males chronically treated with antipsychotics did not find reductions in the total length or diameter of myelinated nerve fibers in prefrontal or global white matter (Marner and Pakkenberg, 2003). However, as unmyelinated fibers were not investigated, a change in these could not be ruled out.

Alternatively, the observed volume reduction could be due to changes in glia cells. Interestingly, several psychiatric disorders have been reported to have decreased glia cell number or density (Cotter et al, 2001, 2002; Rajkowska et al, 2001; Hof et al, 2003; Stark et al, 2004). However, additional studies are needed to determine whether these observed glial changes are specific to the respective diseases or related to treatment with antipsychotic drugs. A growing body of literature describes abnormalities in myelination in schizophrenia (Davis et al, 2003), although again it is not clear to what degree these findings reflect treatment with antipsychotics.

Another candidate for structural changes in both white and gray matter may be the vascular bed. However, as studies of different brain regions in different species have consistently found that the vascular bed comprises less than 5% of the total brain volume (Haug, 1986; Weis et al, 1996; Løkkegaard et al, 2001), a reduction in the vascular bed could contribute to, but not fully account for, the observed volume reductions.

Implication for Studies of Schizophrenia

The results of this study raise the possibility that reported findings of progressive reductions in the volume of certain brain regions in individuals with schizophrenia may be confounded by a treatment effect. A recent meta-analysis of post-mortem studies found a small (2%) but significant decrease in brain weight in schizophrenia (Harrison et al, 2003). Correspondingly, only 11 out of 50 imaging studies found reduced whole brain volume in schizophrenia (Shenton et al, 2001); however, all 11 of the positive studies were based upon chronic patients, while 10 of the negative studies were based upon subjects in their first episode of schizophrenia. In addition, it is not clear to what degree the inconsistency of the MRI findings reflects methodological limitations; assessments of accuracy of MRI-based volume estimates suggest that significant biases may be present (Gadeberg et al, 1999, 2002; Garcia-Fiñana et al, 2003). Assuming it exists, this small difference could to some degree be caused by medication effects like the one observed in this study; however, it appears that a medication-induced difference in whole brain size of a magnitude comparable to the 10% seen in our study does not occur, in general, in individuals with schizophrenia treated with antipsychotic drugs. It should also be noted, though, that a medication effect could be diluted by heterogeneities in subject populations as well as in treatment regimens. In addition, potential drug effects may show a more regional pattern in the human brain, compared to the monkey, again diluting the overall effects. Indeed, longitudinal imaging studies typically have reported progressive changes localized to specific brain regions such as reduced volume of middle temporal lobe (Okubo et al, 2001; Shenton et al, 2001; Niznikiewicz et al, 2003). Interestingly, a recent post-mortem schizophrenia study found a selective and significant 12% reduction in frontal gray matter volume and a 6% reduction in frontal white matter (Selemon et al, 2002), findings in the same range as those found in the current study.

Longitudinal imaging studies have reported increased volume of the caudate nucleus in humans treated with typical antipsychotics (Chakos et al, 1994; Keshavan et al, 1994; Doraiswamy et al, 1995), whereas reduced or unchanged caudate volume has been seen in subjects treated with atypical antipsychotics (Chakos et al, 1995; Gur et al, 1998; Corson et al, 1999; Scheepers et al, 2001). This pattern has been replicated in a recent stereological study of rats, where haloperidol exposure was associated with an increased volume of striatum and olanzapine exposure with a reduction (Andersson et al, 2002). Interestingly, a recent study of subjects with schizophrenia treated with low dosage haloperidol or atypical antipsychotics and normal comparison subjects found similar decrease in volumes of the caudate nucleus in all three groups after a 5-year follow-up (Tauscher-Wisniewski et al, 2002). Unfortunately, in our study, we have not been able to illuminate these findings. The caudate nucleus was divided into several pieces by the initial, cortex-focused dissection, and the largest part of the caudate was contained within the temporal piece. Furthermore, as the caudate comprised only a minor part of this piece and the reported caudate volume changes in humans were of limited magnitude (within ±15%), we do not have the necessary resolution to observe a potential drug-related caudate effect.

Tissue Shrinkage

Our findings also clearly illustrate how a standard histological protocol for cryoprotection and storage of brain tissue induces a significant shrinkage of 20%. This shrinkage occurs despite the gradual stepwise increment in sucrose concentration and the relative moderate final concentration of 18% sucrose; a single step of immersion in 30% sucrose as used by some laboratories may very well induce a greater degree of tissue shrinkage. The observed shrinkage in the current study was the same across all three groups of animals; we did not see any evidence of either global or regional differences in shrinkage due to drug exposure. It should be kept in mind, however, that our findings cannot rule out a distinct difference in post-mortem shrinkage of any minor subcompartment of the tissue (eg of a certain cell type), as such an effect—although pronounced—could be miniscule compared to the total volume and below the resolution of our study. The observation of similar post-mortem shrinkage across all three groups of animals is of importance when interpreting size or density data from human post-mortem studies of subjects treated with antipsychotics. For example, the finding of reduced global or regional brain volumes in subjects with schizophrenia (Pakkenberg, 1987; Selemon et al, 2002) has been discussed under the implicit assumption that the tissue volume was reduced in vivo, without considering the possibility that the brains from subjects with schizophrenia shrunk more during tissue processing. Finally, we observed that different brain regions shrunk differently. This might very well be due to regional differences in structure (such as white and gray matter content), volume, shape, and total external surface (natural and/or cut)—characteristics that are all of obvious importance for the exchange of fluids (predominantly water) that is likely the main cause of the processing-induced shrinkage. This indicates that the use of a global (whole brain) shrinkage correction factor for studies of a particular brain region should be avoided or only implemented and interpreted with very great care.

Precision and Accuracy of the Quantitative Approach

In the present study, we have documented that the implementation of Archimedes' principle is both precise and accurate. Although both Archimedes' and Cavalieri's methods for volume determination are unbiased in principle, they are of course potentially biased in implementation. The Archimedean approach was expected to have a small positive bias due to fluid or small air bubbles trapped in sulci at the surface of the tissue pieces, whereas the results based upon Cavalieri's principle in the present implementation could have a small negative bias due to under projection. The fact that potential biases of the two methods were of opposite directions is crucial for the comparison of the Archimedean results to those based upon Cavalieri's principle. Consistent with our expectations, we found that the combination of the systematic biases of the two methods is below 3%, indicating that both methods are very accurate. Furthermore, Figure 10 illustrates (as expected) that we did not register any observer biases across exposure groups due to lack of blinding during the last set of Archimedean measures.

The reported Cavalieri volumes of gray and white matter are influenced to an unknown (but limited) degree by over and under projections at the boundary between white/gray matter (eg a thin layer of white matter over gray matter may go unnoticed and vice versa). We sought to limit this bias as much as possible by cutting the parietal slabs perpendicular to the intraparietal sulcus, thereby minimizing the amount of cuts oblique to the pial surface. Perpendicular cuts are less prone to over- and underprojection problems due to better visualization of the white/gray matter boundary.

Conclusion

In summary, we found that chronic exposure of monkeys to haloperidol or olanzapine in a manner that mimics clinical use is associated with a significant reduction in brain volume that affects both gray and white matter. In contrast, although substantial and regionally specific reductions in tissue volume occur with histological processing, pre-mortem exposure to antipsychotics does not appear to affect this process.

Further studies are needed to confirm these observations of antipsychotic-related reductions in brain volume, to identify the affected neural elements, and to determine the mechanisms that produce these changes.