Introduction

Schizophrenia is a complex and disabling disorder characterized by psychotic and cognitive symptoms, motivational impairment, and affective dysregulation. Its pathophysiology still remains to be elucidated. There is a growing body of evidence that immune dysfunction is involved in the pathophysiology of schizophrenia. Genetic data, derived from the latest and largest genome-wide association study, identified significant associations with immune pathways in schizophrenia.1 Studies using information from population registries have found increased prevalence of autoimmune disorders in patients with schizophrenia and their family members.2 Furthermore, decreased prevalence of schizophrenia has been reported in men who have used glucocorticoids for somatic disorders,3 and some studies suggest efficacy of immunomodulatory drugs in patients with schizophrenia with the largest effect found in first-episode patients.46 Serum and cerebrospinal fluid (CSF) inflammatory markers, such as cytokines produced by microglial cells, have also been linked with schizophrenia and are increased predominantly in first-episode psychosis.7

Although there is evidence for altered immune function in schizophrenia, it is still unclear which immune processes and immune cells are affected. Various studies using post-mortem tissue of schizophrenia patients and controls described increased density and alterations of microglia.8,9 These findings could not be replicated in other studies, which may be due to small sample sizes and post-mortem confounders.8 Moreover, post-mortem brain tissue of patients in the early phase of the disease is very scarce.

Altered microglia density and phenotype is associated with increased 18 kDa translocator protein (TSPO) expression. Until now, the only method to investigate microglia cells in the brain in vivo is positron emission tomography (PET) using TSPO ligands.10 In classic central nervous system (CNS) autoimmune disorders such as multiple sclerosis increased binding of TSPO ligands was demonstrated.11 However, full-blown neuroinflammation is less evident in schizophrenia.12 PET studies examining TSPO binding in patients with schizophrenia have provided inconclusive results. Using the TSPO ligand (R)-[11C]PK11195, studies have reported increased binding potential in schizophrenia patients in total gray matter13 and hippocampus.14 Recently, significantly increased distribution volume ratios were reported in total gray matter, and frontal and temporal cortex gray matter of patients with schizophrenia, and ultra-high-risk individuals when compared with controls using the tracer [11C]PBR28.15 Other studies using the second-generation TSPO PET ligands [11C]DPA-713, [18F]FEPPA, and [11C]DAA1106, however, did not find significant differences between patients and controls.1618

Important reasons for these inconsistent findings using PET imaging could be differences in patients characteristics and PET methodology.19 For example, the disease duration of previous PET studies varied, as there was one ultra-high-risk assessment,15 two studies examined patients with schizophrenia within 5 years after disease onset,13,16 and other studies examined chronic patients.14,15,17 Clinical trials showed biggest effect of immunomodulatory drugs in the early stages of schizophrenia.5,6 The purpose of the present study was to examine (R)-[11C]PK11195 binding potential (BPND) in total gray matter and several gray matter regions of interest in recent onset psychosis. By using patients with recent onset psychotic disorder and matched controls the role of confounders was minimized.

Results

Demographics of patients with a recent onset psychotic disorder and healthy controls are shown in Table 1 and Supplementary Table 1. Patients displayed moderate symptoms at the time of PET scanning. The average disease duration as well as the average duration of treatment with antipsychotic medication was 1 year. Most patients were on antipsychotic medication at the time of the scanning, apart from 3 patients who were not taking medication and 1 patient who was antipsychotic naive. There were no significant differences in demographics between patients and controls. Although a significant difference in specific activity of (R)-[11C]PK11195 between groups was found, the injected mass of PK11195 did not differ between clinical groups, and no significant associations have been found between specific activity or injected mass and (R)-[11C]PK11195 BPND (data not shown). Moreover, all PET scans were performed at tracer level. Nicotine and cannabis use was not matched across groups, and was more prevalent in the patient group.

Table 1 Demographics

No significant difference in mean (R)-[11C]PK11195 BPND between groups was observed in total gray matter, either with or without age as covariate (Figure 1). Similarly, for regional (R)-[11C]PK11195 BPND no significant effect of clinical group was found (Figure 2, Table 2).

Figure 1
figure 1

Scatterplot showing individual binding potential (BPND) values of the total gray matter. No significant differences were found between patients with a recent onset psychotic disorder (black dots) compared with healthy controls (white dots). Bars indicate mean group BPND values.

Figure 2
figure 2

Scatterplot showing individual binding potential (BPND) values in each of the studied regions of interest (frontal cortex, temporal cortex, parietal cortex, striatum, and thalamus). No significant differences were found between patients with a recent onset psychotic disorder (black dots) compared with healthy controls (white dots). Bars indicate mean group BPND values.

Table 2 (R)-[11C]PK11195 binding potential (BPND)

A significant correlation was found between age and total gray matter (R)-[11C]PK11195 BPND when all data were analysed (N=36, r=0.38, P=0.02; Supplementary Figure 2), as well as age and regional (R)-[11C]PK11195 BPND in all ROIs (N=36, r=0.72–0.95, P<0.01). No significant correlations were observed between total and regional gray matter (R)-[11C]PK11195 BPND and other variables, such as alcohol, nicotine or cannabis use, disease duration, and disease severity (Supplementary Table 2; correlations between (R)-[11C]PK11195 BPND and PANSS subscale scores are available upon request). (R)-[11C]PK11195 BPND did not significantly differ between patients on and off antipsychotics, nor was there a significant difference between antipsychotics (Supplementary Table 3).

Discussion

In this study, microglia activation in gray matter ROIs of patients with a recent onset psychotic disorder were investigated in vivo. There was no significant difference in (R)-[11C]PK11195 BPND between patients and controls in total gray matter, nor in any of the ROIs investigated. In addition, no relationship was found between (R)-[11C]PK11195 BPND and clinical measurements, such as symptom severity.

This study is in line with previous negative findings of three PET studies examining patients with schizophrenia using the second-generation ligands [11C]DPA-713, [18F]FEPPA, and [11C]DAA1106, respectively,1618 but is in contrast to two (R)-[11C]PK11195 studies that have reported significant increases in (R)-[11C]PK11195 binding in total gray matter13 and hippocampus14 in patients with schizophrenia, and to one [11C]PBR28 study that reported increased distribution volume ratios in total gray matter and frontal and temporal cortex gray matter of patients with schizophrenia and ultra-high-risk individuals.15

In this study, whole-brain ROIs were selected because of the widespread TSPO availability in the brain. Combined ROIs were selected to avoid type I errors, which resulted in the inclusion of the hippocampus in the temporal region. When examining the hippocampus alone, there was also no significant difference (data not shown). Previous PET studies that examined widespread brain ROIs in schizophrenia,14 found a localized effect in temporal and frontal ROIs only.14,15

There can be several explanations for the discrepancy across the TSPO PET studies. First of all, this discrepancy in findings may be due to the differences in disease duration. Microglia can shift in function and morphology during the course of the illness20 and may be involved in (later stages of) schizophrenia rather than being associated with recent onset psychosis.

The patients in the present study have a short disease duration, which resulted in the inclusion of patients with a psychotic disorder not (yet) fulfilling criteria for schizophrenia. Our findings are in line with a recent study that reported no differences distribution volume (VT) values in recent onset schizophrenia.16 In addition, no differences in VT values were found in ultra-high-risk individuals using two-tissue compartmental model (2TCM).15 However, if analysed with a 2TCM accounting for endothelial vascular TSPO binding (2TCM-1K), this study did find an effect in ultra-high-risk individuals,15 and another (R)-[11C]PK11195 study found an effect in patients with an illness duration of less than first 5 years suggesting that microglia activation could be present in the early phase(s) of the disease.13

Second, a positive association between symptom severity and TSPO expression has been reported in some of the previous studies.15,17 This suggests that the level of symptomatology could be a reason for the discrepant results. However, as symptom severity in the present study was not different from that in previous PET studies reporting an effect, this seems unlikely.1317

Third, general concerns with the inclusion of schizophrenia patients are antipsychotic treatment and substance use, such as nicotine and cannabis use. Indeed, it has been shown in vitro that TSPO expression is reduced in cell lines treated with most atypical antipsychotics, apart from clozapine, which tended to increase TSPO expression.21 In the present study, no significant differences between patients on or off antipsychotics were found, but subgroups were small. Future studies may be able to address this issue by examining antipsychotic naive patients, which so far has only been performed in ultra-high-risk individuals.15 In addition, nicotine and cannabis use might reduce TSPO expression, as they have a suppressant effect on the immune system.22,23 Nicotine and cannabis use was more prevalent in the patient group. However, no significant association between nicotine or cannabis use and (R)-[11C]PK11195 BPND was observed.

Furthermore, the differences might be explained by differences in methodology.19 (R)-[11C]PK11195 has the advantage that it has been extensively investigated due to more than 30 years of research.24 A limitation is its lower signal-to-noise ratio in comparison with second-generation tracers, which has been demonstrated by animal blocking studies showing up to 80-fold higher specific binding in second-generation ligands, such as [11C]PBR28 and [18F]FEPPA, compared with (R)-[11C]PK11195.2527 In addition, biomathematical modeling of in silico/in vitro data indicated a higher predicted within-subject variability for (R)-[11C]PK11195 than second-generation ligands, such as [11C]PBR28 and [11C]DPA-713.28 Direct in vivo comparison between (R)-[11C]PK11195 and second-generation tracers in humans is limited, although the second-generation ligand [11C]DPA-713 showed higher VT values in comparison to (R)-[11C]PK11195 in a small group (N=2) of healthy volunteers.29 The negative findings of the present study seemed not to be tracer specific as previous (R)-[11C]PK11195 have reported an effect,13,14 and second-generation tracers have reported both negative and positive findings.1518 The two-tissue reversible plasma input model (2T4k) is the gold standard for analysis of dynamic (R)-[11C]PK11195 studies.30 However, plasma input models are more invasive than reference tissue methods. In general, there is no anatomical reference ROI available, as there is no brain region devoid of TSPO. If a reference region is not completely devoid of TSPO, this may result in underestimation of the signal. For instance, (R)-[11C]PK11195 may be affected by non-specific binding to endothelial tissue.31 However, no differences in time-activity curves of the reference region between patients and controls were found (data not shown). In addition, to avoid the effect of binding to venous sinuses, a vascular correction was applied in the present study.31,32

Lastly, a general problem is the large intersubject variability in TSPO expression, which is shown in this and other PET studies.1416 To correct this variability, normalization to whole-brain outcome was used in two studies that did report an effect, either by the use of the whole-brain gray matter binding potential as a covariate to correct for global uptake14 or by calculating the distribution volume ratio, which is defined as the ratio of VT in the ROI to VT in the whole-brain.15 Interestingly, the present results are in line with the (non-normalized) VT values of one study, as no difference in this outcome measure was found.15 It needs to be stressed that the intersubject variability in the present groups was probably too large to detect significant differences. Therefore, larger sample sizes may be needed to detect a possible effect between groups. The sample size was relatively large for a PET study and according to the sample size calculation of a previous data set,13 which was reanalysed using the same method as in the present study (TFvdD et al., unpublished), it should have been large enough to detect a clinically significant difference using the present method. However, microglia activation might be a relatively subtle phenomenon in psychosis, which is not easily detected.

In the present study, no corrections for partial volume effects were applied, as no reduction of gray matter volume in patients was found, nor a correlation with BPND (data not shown). In addition, a previous TSPO study did not report differences between partial volume corrected and uncorrected data in psychosis, indicating that partial volume effects are minor.17

The inclusion of patients with a short disease duration is the main strength of the present study. A limitation is the fact that most patients were on (atypical) antipsychotics. A general limitation for (R)-[11C]PK11195 and other TSPO ligands is that the actual cellular process underlying the PET signal is unclear.20 TSPO ligands do not exclusively bind to microglia, but may also bind to other cells in the CNS, such as astrocytes.33,34 Furthermore, for the microglia population, TSPO ligands cannot distinguish between pro-inflammatory and anti-inflammatory subpopulations. Future studies may develop more specific microglia markers that are also sensitive to the functional phenotype.19

In conclusion, no evidence was found with increased microglial activation in gray matter in recent onset psychosis. These findings suggest that microglia activation is not present in recent onset psychosis or that it is a subtle phenomenon that could not be detected due to other factors, such as patient characteristics and antipsychotic treatment. Further studies, preferably using more sensitive TSPO tracers, are warranted to ascertain the role of microglia in the pathophysiology of schizophrenia.

Materials and methods

Subjects

Overall, 19 patients with a recent onset psychotic disorder and 17 age and sex-matched healthy controls were included. Patients were recruited from the Department of Psychiatry of the University Medical Center Utrecht (UMCU), Utrecht, The Netherlands, and the Department of Psychiatry of the Academic Medical Center (AMC), Amsterdam, The Netherlands. Controls were recruited via advertisements. Eligible patients were aged 18–40 years, met diagnostic criteria for a psychotic disorder, and were studied within the first 2 years of duration of treatment or were antipsychotic naive. Duration of treatment was defined as first intake of antipsychotic medication. DSM-IV diagnoses were confirmed by two investigators using the Comprehensive Assessment of Symptoms and History (CASH) interview.35 Exclusion criteria were any neurological disorder, pregnancy, abnormal score on routine laboratory tests (blood tests and urine analysis), intake of non-cannabis addictive drugs or use of benzodiazepines and non-steroidal anti-inflammatory drugs (NSAIDs) (4 weeks before start of the study), and current and lifetime alcohol abuse.

Additional exclusion criteria for healthy controls were presence or history of psychiatric disorder confirmed by the CASH interview or a first-degree relative with a family history of schizophrenia spectrum disorder. The study complied with the Declaration of Helsinki and had been approved by the Medical Ethics Review Committee of UMCU. All subjects signed an informed consent form before inclusion.

Demographical data were collected within 1 week of performing the PET scan. Current and lifetime use of alcohol, nicotine, cannabis, and other illicit drugs was assessed using the Composite International Diagnostic Interview (CIDI)36 and disease severity with the positive and negative syndrome scale (PANSS).37 Complete blood count, serum chemistry, use of benzodiazepines and addictive drugs (3, 4-methylenedioxymethamphetamine, amphetamine, cannabis, cocaine, opiates), and pregnancy were tested with laboratory tests within 1 week of the PET scan.

PET procedures

PET scans were performed on an Philips Gemini TF-64 positron emission tomography/computed tomography (PET/CT) scanner (Philips Medical Systems, Best, The Netherlands).38 Tracer administration was standardized. After administration of a bolus injection of 370 MBq (R)-[11C]PK11195 using an infusion pump (Med-Rad, Beek, The Netherlands), a dynamic three-dimensional (3D) emission scan of 60.5 min was acquired as described previously.31 Before the emission scan, a low-dose computed tomography (CT) scan was acquired on the PET/CT scanner, which was used to correct the subsequent PET scan for tissue attenuations. Using laser beams, subject motion during scanning was checked visually at regular intervals and corrected immediately, whenever necessary. Moreover, motion was minimized by using a head immobilization device (head holder). PET data were corrected for dead time, decay, tissue attenuation, scatter, and randoms. PET images were reconstructed using a dedicated brain mode iterative algorithm (LOR-RAMLA)37 and had a voxel size of 2×2×2 mm3.

MRI scans

A structural T1-weighted magnetic resonance imaging (MRI) scan was acquired for each subject. Subjects were scanned on a Philips Achieva 3T scanner (Best, The Netherlands; voxel size 0.75×0.75×0.75 mm3). MRI scans were examined by a neuroradiologist to exclude major brain abnormalities in both patients and controls.

Image analysis

PET data were analysed using a reference tissue approach as described previously.31,32,39 In short, this procedure uses supervised cluster analysis algorithm with four kinetic classes (SVCA4) to extract the reference tissue input function. SCVA4 first applies an anatomical mask to include brain tissue only and then the algorithm clusters voxels in the following predefined kinetic classes: gray matter with specific (R)-[11C]PK11195 binding, gray matter without specific binding, white matter, and blood. Gray matter voxels without specific binding were used as reference tissue input. Next, BPND parametric images of (R)-[11C]PK11195 scans were generated using reference parametric mapping with a vascular correction (RPMVb),40 which is a basis function implementation of the simplified reference tissue model.41

T1-weighted MRI scans were co-registered with PET scans and used for segmentation of gray and white matter. Then, the co-registered MRI scan was used to delineate total gray matter, and to define five gray matter ROIs including the frontal cortex (volume weighted average of anterior cingulate cortex, gyrus rectus, orbitofrontal cortex, precentral gyrus, superior frontal cortex, middle frontal cortex, interior frontal cortex), temporal cortex (volume weighted average of hippocampus, anterior and posterior temporal cortex, parahippocampal, and ambient cortex, superior, middle and inferior temporal cortex, fusiform cortex), parietal cortex (volume weighted average of posterior cingulate cortex, postcentral gyrus, superior parietal cortex, and inferolateral parietal cortex), thalamus and striatum (volume weighted average of caudate nucleus and putamen) using probability map-based automatic brain delineation.42 This selection of ROIs was based on the widespread TSPO availability in the brain. The outcome measure was (R)-[11C]PK11195 binding potential (BPND; Supplementary Figure 1).41,43

Statistics

Analysis of variance (ANOVA) and multivariate analysis of variance (MANOVA) with ROI as within-subject factor, group as between-subject factor and age as covariate were carried out to test for differences between groups (two-sided) in total gray matter and regional (R)-[11C]PK11195 BPND, respectively. A subsequent MANOVA as a post-hoc analysis was implemented to identify regional group differences between patients and controls (two-sided). Homogeneity of variances was tested by the Levene's test for equality of variances and was not violated for (R)-[11C]PK11195 BPND in all regions tested.

Pearson’s correlation coefficient (r) was used to evaluate correlations between (R)-[11C]PK11195 BPND and other variables, such as demographics. The χ2 test was used to test between-group differences in gender. P values below 0.05 were considered to be significant and Bonferroni corrections were used for multiple comparisons. The power calculation was based on previous data,13 which were reanalysed using the same method as in the present study (using a reference tissue approach and a supervised cluster analysis algorithm to identify the reference region, TFvdD et al., unpublished). This resulted in the requirement of at least 16 subjects per group to detect a difference with 80% power to examine five ROIs (μ1=0.07, μ2=0.18, σ1=0.10, σ2=0.08, α=0.01).