Pre-frontal parvalbumin interneurons in schizophrenia: a meta-analysis of post-mortem studies

1Schizophrenia is a severe mental disorder characterised by positive, negative and cognitive symptoms (Howes and Murray; Owen et al. ). It has ranked 12th in the leading global causes of disability for the last decade (Wang et al. ). Current pharmacological treatments, which are all dopamine receptor ligands (Howes et al. ), are ineffective for many patients, highlighting the need for a better understanding of the neurobiology underlying the disorder to develop more effective treatments (Howes and Kapur; Nutt and Need ). A number of neurotransmitters, including dopamine, glutamate and gamma-aminobutyric acid (GABA), are implicated in the pathophysiology of schizophrenia (Frohlich and Van Horn; Howes et al.; Krystal et al. ). GABAergic interneurons provide inhibitory control of cortical and subcortical circuits (Curley and Lewis; Konradi et al. ) and are, therefore, thought to lead to glutamatergic and dopaminergic dysfunction, which in turn leads to the symptoms of schizophrenia (Grace ). Abnormal GABAergic interneuron activity is thought to play a key role in working memory and other cognitive symptoms of schizophrenia (Frankle et al.; Gonzalez-Burgos et al.; Heckers and Konradi; Lewis et al.; Perry et al.; Stedehouder and Kushner; Taylor and Tso ).

2Although more than 20 different classes of GABAergic interneurons exist (DeFelipe et al. ), it has been suggested that in schizophrenia there is a loss of a particular type of GABAergic interneurons characterised by containing the calcium-binding protein parvalbumin (PV) (Beasley and Reynolds; Chung et al.; Enwright Iii et al.; Lisman et al.; Lodge and Grace; Perez and Lodge ). Parvalbumin interneurons are fast-spiking neurons that are critical to the generation of gamma oscillations, a type of high-frequency neuronal oscillation linked to working memory and other cognitive processes in healthy subjects and patients with schizophrenia (Gonzalez-Burgos et al.; Sohal et al.; Tallon-Baudry et al.; Uhlhaas and Singer ). So far, several post-mortem studies have examined parvalbumin interneurons in patients with schizophrenia, though the findings have not all been consistent (Chung et al.; Gonzalez-Burgos et al. ) with some early studies showing decreases in parvalbumin interneuron density in schizophrenia (Beasley and Reynolds; Reynolds et al. ) whilst others failed to find significant differences between patients and controls (Enwright et al.; Hashimoto et al.; Woo et al. ).

3Contemporary interpretations of the literature highlight the methodological limitations in measuring parvalbumin cell density; in particular, the danger of using immunoreactivity labelling techniques when protein expression is decreased (Enwright Iii et al. ). The prevailing view, therefore, remains that in schizophrenia parvalbumin interneuron abnormalities exist not at the level of neuronal morphology or density, but at the molecular level of gene expression and protein synthesis (Chung et al.; Enwright Iii et al.; Enwright et al.; Fung et al.; Hashimoto et al.; Joshi et al.; Volk et al. ). However, Toker et al. have challenged this view in a large transcriptomic study on bulk post-mortem tissue, using cell type-specific marker genes to indirectly show reduced parvalbumin cell density in schizophrenia. As far as the authors are aware, there has not been a previous meta-analysis of parvalbumin post-mortem studies in the pre-frontal cortex. In view of this, we aimed to systematically review and meta-analyse the available evidence on parvalbumin neuronal density and mRNA brain measures in the pre-frontal cortex in people with schizophrenia.

4A standardised search was conducted according to the PRISMA guidelines (Liberati et al. ). The Medline, EMBASE and PsycINFO electronic databases were searched. The electronic search using EMBASE and PsycINFO was carried out together using Ovid. The following keywords were used “(schizophrenia OR psychotic disorders OR psychosis) AND (parvalbumin)” until 8th February 2019 (for details of search strategy, see supplementary information).

5Inclusion criteria were (1) original studies reporting post-mortem findings of parvalbumin neuronal density or parvalbumin mRNA levels in the pre-frontal cortex (2) included patients with a confirmed diagnosis of schizophrenia (including schizoaffective disorder) (3) reporting mean and variance measures for patient and control groups. Exclusion criteria were pre-natal studies, animal studies, and articles in a language other than English, review articles not reporting original data, single cases reports, studies without a control group and articles that were not published in peer-reviewed journals. Studies that only measured GABA-related transcripts or markers but not parvalbumin protein or mRNA were excluded. We did not include studies that did not measure parvalbumin cell density and/or measured mRNA in regions other than the pre-frontal cortex, e.g. Konradi et al.. We also did not include micro-array studies looking at pre-frontal regions (n = 4) (Catts and Weickert; Enwright Iii et al.; Fung et al.; Mellios et al. ) as following exclusion for sample overlap or reuse of data, there were insufficient studies (n = 2) (Fung et al.; Mellios et al. ) for a quantitative meta-analysis; however, the results of these studies are discussed later in the paper.

6The study has two outcome measures: (1) the effect size for the difference in parvalbumin neuronal density between patients and controls, i.e. the density of cells defined by a presence or absence of detectable levels of parvalbumin protein or parvalbumin mRNA [either in neuron/mm2 or neuron/mm3 or neurons per identical volumes of sampled tissue in the case of Chung et al.]; 2) the effect size for the difference parvalbumin mRNA levels per neuron or in total grey matter between patients and controls. For all studies we extracted the mean and standard deviation for neuronal density and parvalbumin mRNA for each pre-frontal region in patients with schizophrenia and healthy controls using the available published data. WebPlotDigitizer (https://automeris.io/WebPlotDigitizer) was used to extract data from plots when necessary. If cell density or mRNA was measured across different laminar layers, these data were combined into one simple mean for the region. The online software package StatsToDo (StatsToDo 2014) was used to combine n, standard deviation (sd) and means for each laminar layer for each region into one overall n, mean and sd for each group (patient or control) in each region. In addition, we extracted the following variables: sample size, mean age, mean post-mortem interval, psychiatric medication, mode of death and methodology of the study. We took particular care to ensure studies which shared brain collections did not include the same brains and, in those studies, included in the quantitative meta-analysis that shared the same collection contact was made with the lead author to confirm no overlap. This led to the exclusion of a number of studies with overlapping brains (Beasley et al.; Bitanihirwe et al.; Catts and Weickert; Chung et al.; Enwright Iii et al.; Hashimoto et al.,; Knable et al.; Lewis; Torrey et al. ). In the case of the Hashimoto et al. and Volk et al. studies, means and standard deviations were calculated using raw data from the latter study that excluded the overlapping subjects. All studies included in the quantitative meta-analysis excluded patients with a history of or post-mortem findings suggestive of a neurological or neurodegenerative disease.

7RStudio statistical software running the “metafor” package was used to perform the meta-analysis (RStudio Team 2015). Group size, mean and standard deviation were used to determine standardised mean effect size (Hedges’ g). A random-effects meta-analytic model was used that did not assume homogeneity of effects. Heterogeneity was measured by calculating I2. Meta-regression was used to assess the effects of year of publication and PMI. Publication bias was tested for using funnel plot asymmetry and regression test. Where potential publication bias was suspected, trim and fill analysis was conducted to correct for putatively missing studies. A p value < 0.05 (two tailed) was taken as a significance level. Sub-analyses by fixation method (paraformaldehyde or formalin/paraffin) and by laminar layers III and IV were performed on the pre-frontal cortex cell density data (see Figs. 3 and 4).

8The searches identified 1188 publications, once duplicates were removed (see Fig. 1 for a PRISMA diagram of the literature search). Following screening, 24 post-mortem studies were identified (see Table 1 Supplementary Information for study characteristics). In total, 12 studies were conducted in pre-frontal regions, so were taken forward to the meta-analysis. Tables 1 and 2 summarise the clinico-demographic characteristics of the subjects included as well as the methodological approaches used, in the meta-analyses of parvalbumin cell density (n = 9) and parvalbumin mRNA (n = 4), respectively. This final sample (n = 13—Hashimoto et al. presents data for both parvalbumin cell density and mRNA measures) included 274 patients with schizophrenia (mean age: 49.8) and 275 healthy controls (mean age: 50.6).

9Nine studies (comprising 136 schizophrenia patients and 138 healthy controls) (Table 1) measured parvalbumin neuron density in the pre-frontal cortex. There was a significant reduction in parvalbumin cell density in the pre-frontal cortex regions of patients with schizophrenia relative to healthy controls (Hedges’ g = − 0.27; z = − 2.17; p = 0.03; 95% confidence interval (CI): − 0.51 to − 0.03) (Fig. 2). The I2 test revealed low heterogeneity (I2 = 0%; 95% CI: 0–57.53%). Visual inspection of the funnel plot analysis revealed one possible missing study on the left side (Supplementary Fig. 1). The funnel plot regression test was not significant (t = 0.28, df = 7, p = 0.8). After correction for potential publication bias the effect size reduction seen in parvalbumin cell density in patients with schizophrenia remained significant (Hedges’ g = − 0.29; p = 0.01). A sub-analysis by fixation method of the parvalbumin cell density studies in the pre-frontal cortex showed that those studies using paraformaldehyde, a more reliable fixation method, showed a non-significant reduction in patients with schizophrenia compared with healthy controls (Hedges’ g − 0.30; p = 0.05; CI: − 0.59–0.02) as did the formalin/paraffin-prepared studies (Hedges’ g − 0.19; p = 0.35; CI: − 0.60–0.22) (see Fig. 4). Leave-one-out analysis revealed that the results became non-significant if the Hashimoto et al. or the Bitanihirwe and Woo studies were taken out (see Supplementary Information).

10Four studies (comprising 138 schizophrenia patients and 137 healthy controls) (Table 2) measured parvalbumin mRNA in pre-frontal regions. The meta-analysis showed a non-significant reduction in parvalbumin mRNA in patients with schizophrenia when compared to healthy controls (Hedges’ g = − 0.44; z = − 1.56; p = 0.12; 95% CI − 0.99–0.12) (Fig. 3). The I2 test revealed high heterogeneity suggesting significant differences between studies (I2 = 79.56%; 95% CI: 36.8–98.49%). Visual inspection of the funnel plot revealed one missing studies on the right-hand side and the regression was not significant (t = − 0.15; df = 2, p = 0.89) (Supplementary Fig. 1). Following correction for publication bias, the effect size decrease in parvalbumin mRNA in patients with schizophrenia remained unchanged (Hedges’ g = − 0.44; p = 0.12).

11We found a significant reduction in parvalbumin cell density and a non-significant reduction in parvalbumin mRNA in patients with schizophrenia relative to healthy controls in pre-frontal cortical regions. Two studies of parvalbumin cell density in the frontal cortex showed large effect sizes, those by Hashimoto et al. and Bitanihirwe and Woo; the latter uniquely used a cRNA GAT-1 probe and a DIG-labelled PV mRNA probe to measure the density of neurons expressing both parvalbumin and GABA transporter 1 (GAT-1) mRNA in cortical layers 3 and 4 (data included in this meta-analysis) in patients with schizophrenia. Therefore, this study measured a specific molecular subtype of parvalbumin interneuron and suggests that decreases in parvalbumin interneuron density in schizophrenia could be larger within discrete molecular classes of parvalbumin expressing interneurons. However, the study by Hashimoto et al. used a parvalbumin mRNA probe alone to identify parvalbumin-positive cells, and showed a greater effect size in this random-effects meta-analytic model, suggesting that reductions in parvalbumin cell density are likely to be generalised across different molecular sub-types of parvalbumin interneuron.

12Our finding of a non-significant reduction in parvalbumin mRNA in frontal regions is unexpected, in particular given our finding of reduced parvalbumin cell density and the existence of micro-array studies which do show significant reductions in parvalbumin mRNA fontal regions (Enwright Iii et al.; Mellios et al. ). Unfortunately, there were insufficient micro-array studies measuring parvalbumin expression in frontal regions to be included in a meta-analysis.

13In terms of limitations, combining data from different studies in the context of meta-analyses lead to several potential sources of heterogeneity, such as methodological differences in labelling and microscopy techniques, and sample differences in illness duration, co-morbidity or prior treatment history of the patients included. All included studies took rigorous steps to ensure a reliable diagnosis of schizophrenia was confirmed according to established criteria, though several studies included patients with schizoaffective disorder as well as co-morbid substance misuse or dependence. Co-morbid alcohol dependence could be a significant confounding factor, in particular when looking at hippocampal pathology due to chronic alcohol-related changes in hippocampal size and volume (Agartz et al. ), though ultimately there were insufficient hippocampal studies to include in a quantitative meta-analysis. Co-morbidity within controls was also present in at least one study (Hashimoto et al. ). The method of tissue fixation, in particular the use of formalin and paraffin, may lead to a reduction in immunoreactive labelling of parvalbumin protein due to the effects of epitope masking, degradation and cross-linking (Ahram et al.; Hoetelmans et al. ). Consistent with this, the sub-analysis by fixation method in the frontal cortex, suggests that effects were larger in the paraformaldehyde group than the formalin group, suggesting future studies should use paraformaldehyde fixation (see Fig. 4). Furthermore, Enwright et al. suggest that studies which use monoclonal antibodies for parvalbumin protein to measure parvalbumin cell density may overestimate the decrease in parvalbumin cell density due to low protein levels rendering some parvalbumin cells undetectable, except at much higher levels of magnification. Furthermore, it is likely that the stereotactic three-dimensional counting techniques used by later studies provide a more reliable counting method because cell size can be more accurately measured when calculating density. However, Sakai et al. noted that accurate estimates of cell density can be made using two-dimensional counting methods if the Abercrombie (Abercrombie and Johnson ) correction is used to adjust for the confounding effect of cell size.

14In our calculation of mean differences between patients and controls we averaged findings across grey matter layers for each region of the frontal cortex to obtain an overall mean and SD for that particular region. Parvalbumin interneurons are predominantly found in the middle and lower cortical layers (Hof et al.; Tooney and Chahl ) and have been shown to be most affected in layers 3 and 4 in schizophrenia (Chung et al.; Hashimoto et al.; Sakai et al.; Tooney and Chahl ), and it maybe that averaging across layers underestimates the greater differences found between patients with schizophrenia and controls. We performed a sub-analysis for the studies that measured only layers 3 and 4 and found no significant difference between controls and patients for parvalbumin cell density in the frontal cortex (Hedges’ g = − 0.34 z = − 1.40 p = 0.16 95% CI: − 0.81–0.13—see Fig. 5).

15Cause of death did differ between patients and controls with nearly all studies having an excess of suicides within the patient sample, except Sakai et al. which excluded cases of suicide. Cause of death is known to be a significant potential confounding factor in post-mortem studies as agonal factors such as hypoxia, head injury and ingestion of neurotoxic substances are known to have significant effects on RNA integrity in the post-mortem brain (Tomita et al. ). However, it is not known whether parvalbumin interneurons are susceptible to this effect. Brain pH, an inverse correlate with agonal state, was controlled for within each study and was not found to be significantly different between patients and controls. Post-mortem interval was significantly different between patients and controls in one early study Reynolds et al.; however, the authors found no correlation between PMI and neuronal density. A meta-regression was performed with PMI as a co-variant in the studies measuring parvalbumin cell density and the relationship was not significant (z = − 0.05, p = 0.96) nor was year of publication (z =0.03, p = 0.97). However, in the meta-regression of parvalbumin mRNA studies post-mortem interval was significantly associated with the effect (Hedges’ g = 0.59, CI 0.29–0.89, p = <0.001). This suggests that the two studies which showed less of an effect in patients, Fung et al. and Joshi et al., may have done so because the difference in post-mortem interval was greater in patients versus controls. Again, year of publication was not significantly associated with the effect (z = 0.55, p = 0.58). It should also be noted that of the four studies included in the parvalbumin mRNA meta-analysis, the two which showed a difference in patients were from the Lewis lab, whereas the two which showed no difference were from the Weickert lab (see Fig. 3), so laboratory maybe a confounder for this meta-analysis.

16Antipsychotic medication exposure in patients is a potential source of heterogeneity in findings, although antipsychotic treatment has not been associated with parvalbumin interneuron loss or changes in parvalbumin transcript levels in pre-clinical models (Hashimoto et al.; Lewis et al. ).

17A final methodological consideration is the effect of statistical test. Hashimoto et al. found a reduction in parvalbumin mRNA-positive neurons that was not significant using an analysis of covariance (ANCOVA) model. However, when the same mean and sd data were inputted into the meta-analysis, it was found to exert a significant effect. This suggests that when other confounding factors are taken into account, the effect size showing a reduction in patients is reduced.

18Post-mortem studies of parvalbumin interneurons in schizophrenia initially focussed on interneuron morphology (Beasley and Reynolds; Benes et al.,; Bernstein et al.; Cotter et al.; Danos et al.; Falkai et al.; Kalus et al.; Konradi et al.; Pantazopoulos et al.,; Reynolds et al.; Wang et al.; Woo et al.; Zhang et al. ), while later studies began to measure parvalbumin-related proteins, mRNA and other molecular markers such as GAD67, GAT-1, ErB4q and PGC-1α (Bullock et al.; Byne et al.; Chung et al.; Fung et al.; Hashimoto et al.; Joshi et al.; Lucas et al.; McMeekin et al.; Mellios et al. ). Our findings extend other evidence of GABAergic dysfunction, such as reduced levels of GAD67 in the frontal cortex (Fung et al.; Volk et al. ), to identify a particular sub-class of GABAergic interneuron, parvalbumin-positive interneurons, as being affected in schizophrenia, although they do not exclude other sub-types also being involved.

19Although our meta-analysis of parvalbumin mRNA in the pre-frontal cortex found a non-significant reduction, a recent large micro-array study of parvalbumin mRNA per parvalbumin-containing neuron in layer three of the dorsolateral pre-frontal cortex found a significant reduction of 22.2% in patients with schizophrenia (Enwright Iii et al. ), although it should be noted that this study contains sample overlap with a number of the group’s other studies in frontal regions (see “Materials and methods”), including Volk et al. which is included in our PV mRNA meta-analysis. In addition, reduced interneuron PV mRNA appears to be associated with a concurrent loss of GAD67 mRNA (Chung et al.; Hashimoto et al. ), which could be due to defects within the PV gene (located 22q12–q13.1) or a factor that regulates gene expression. Toker et al. who used micro-array expression profiles of cell type marker genes in bulk post-mortem tissue as surrogate measures of cellular abundance suggest that their findings demonstrate the existence of a reduction in parvalbumin cell numbers in schizophrenia. Such findings suggest a deficiency secondary to either neuronal loss or altered parvalbumin transcription (Chung et al.; Glausier et al. ).

20Two main morphological sub-types of parvalbumin interneurons have been characterised within the neocortex: fast-spiking chandelier or axo-axonic cells, which target the axon of pyramidal cells, and fast-spiking basket cells that synapse onto pyramidal cells in the perisomatic basket region (Cobb et al.; Klausberger et al. ). Parvalbumin basket interneurons receive direct thalamic input (Jones ) and are thought to represent 80–90% of parvalbumin interneurons (Zaitsev et al. ). Such cells play a primary role in timing and synchronising pyramidal cell discharge through regulating rhythmic hyperpolarisation and local feedback loops (Bartos et al.; Hu et al.; Somogyi and Klausberger; Tremblay et al. ). In doing so, they contribute to global and regional neuronal oscillations, within the gamma range (Gonzalez-Burgos et al.; Senkowski and Gallinat ). Abnormalities in both gamma oscillation synchronicity and power are known to be present in schizophrenia (Gallinat et al.; Kwon et al.; Matthew et al.; Spencer et al. ). Therefore, the reduction in parvalbumin interneurons we found in frontal regions in schizophrenia could contribute to the gamma power and synchrony alterations seen in the disorder and the disrupted cortical–thalamic connections which are thought to underlie cognitive dysfunction found in schizophrenia (Cho et al. ).

21Our finding of reduced parvalbumin cell density in the pre-frontal cortex taken with other evidence implicating parvalbumin interneuron dysfunction (Enwright et al.; Sohal et al.; Volk et al.; Zhou et al. ), suggests that targeting parvalbumin interneurons may have potential as drug targets in the treatment of schizophrenia. One such possibility is agents that modulate Kv3.1/2 channels that are highly expressed on fast-spiking parvalbumin interneurons (Boddum et al.; Brown et al.; Rosato-Siri et al. ). Future post-mortem work on the molecular differentiation of parvalbumin interneuron subclasses and the regulational factors that act upon them will be needed to determine the true nature of parvalbumin interneuron pathology in schizophrenia.

22In summary, this meta-analysis finds that in schizophrenia parvalbumin interneurons in the frontal cortex show reduced neuronal density but the reduction in parvalbumin mRNA levels was not found to be significant. Both these findings should be viewed in the context of several methodological considerations and potential confounding factors. This highlights the need for further robust post-mortem studies in this and other brain regions, which take into account such factors. These findings add to the evidence that pathology of discrete neuronal classes exists in schizophrenia and that GABAergic interneurons play a significant role in the pathoetiology of the illness.

23Below is the link to the electronic supplementary material.