The main findings of this meta-analysis are that schizophrenia is associated with increased variability in the concentrations of glutamatergic metabolites in the brain, together with regional differences in mean glutamatergic metabolite concentrations. The greatest amount of data were available for the MFC, where glutamate levels were lower and the variability of all glutamate metabolites (glutamate, glutamine and Glx) were increased in schizophrenia compared to controls. While the emphasis has been on glutamate in the MFC in schizophrenia, our results additionally indicate glutamatergic dysregulation in other brain regions, finding increased variability in glutamine and Glx in the DLPFC, higher levels of Glx in the basal ganglia, and higher levels of glutamine alongside increased variability of glutamate and Glx in the thalamus. The increased variability in glutamatergic metabolites tended to be most apparent in studies examining patients who were younger or more symptomatic.
Glutamate is tightly regulated in the brain through complex feedback mechanisms which may be disrupted in schizophrenia (65). The increased variability in glutamatergic metabolites in schizophrenia indicate a spectrum of disturbances in glutamate homeostatic control between individual patients, resulting in a wide range of concentration values. Regionally, our analysis of the available data in schizophrenia indicates dysfunctional regulation (i.e. increased variability) of glutamatergic metabolite concentrations in the MFC, DLPFC and thalamus, and mean value decreases in the MFC and increases in the thalamus and basal ganglia. Increased glutamatergic variability in schizophrenia, as observed at the macro scale with 1H-MRS, could relate to individual differences in the nature or in the extent of the underlying molecular pathophysiological mechanisms. We did not find evidence for a bimodal distribution within individual participant data, as would be expected if between-subject glutamatergic differences were driven by the presence or absence of a simple, discrete variable. In contrast, the observed unimodal distributions are consistent with the view that glutamatergic pathology in schizophrenia arises secondary to a range of factors that vary amongst patients (e.g polygenic and multiple environmental vulnerability variables)(7–9). Therefore, findings of higher glutamate levels in treatment resistant patients in comparison to treatment responders(15–24) may represent a continuum, rather than the presence of discrete subtypes.
There was evidence in some regions that the variability of glutamate and glutamine was greater in studies including younger or more symptomatic patients. The variability of MFC glutamine and basal ganglia glutamate in patients with lower symptom severity appeared to be more like that seen in controls. This may indicate corrective regulatory mechanisms over time, with the distribution becoming more uniform with age and its co-segregates (e.g. medication or illness duration) or possibly as symptoms improve. In the case of basal ganglia glutamate this is supported by the finding that variability is also lower in older patients. Our results indicate that in older patients, glutamate in the temporal lobe and basal ganglia may be even more homogenous than that seen in controls, which could potentially reflect over-compensation. Alternatively, the greater glutamatergic variability in studies examining younger, more symptomatic patients may reflect greater clinical heterogeneity within first episode cohorts, for example in respect to diagnosis or severity of symptoms. It should be noted that the relationship between basal ganglia glutamate and symptom severity is largely influenced by a single study, and thus further studies in highly symptomatic patients are needed to confirm this finding.
As increased variability was of comparable magnitudes in antipsychotic-naive and medicated patient cohorts, this suggests that variability in glutamate does not result from differential effects of antipsychotic medication on glutamate levels between individuals(14, 21). Although a recent meta-analysis of DLPFC glutamate found higher variability in antipsychotic medicated patients and lower variability in medication-naïve patients(29). Our analysis, including 12 more recent studies and investigating DLPFC glutamate and Glx separately, did not find any effect of medication status in this region. In fact, in the basal ganglia, it appeared that antipsychotic-treated cohorts displayed reduced variability compared to controls. Potentially this could be due to a regulatory effect of antipsychotics on basal ganglia Glx(66, 67).
The meta-analysis of standardised mean differences found lower MFC glutamate and higher thalamic glutamine and basal ganglia Glx in patients compared to controls. These results are consistent with recent meta-analyses(2–6), but substantially extend them by including 14 new datasets for MFC glutamate and 2 new datasets for thalamic glutamine. Stratified analyses found that lower MFC glutamate was observed across studies examining medicated patients but not across studies examining antipsychotic-naïve patients. Thus, antipsychotic medication could lower MFC glutamate levels, as indicated by longitudinal studies(14, 21) and a mega-analysis(12), although the meta-regression with CPZ dose was not significant. Our analysis also revealed relationships between the proportion of males in the study and MFC glutamate and frontal white matter Glx effect sizes, such that a higher proportion of males was associated with lower glutamate levels in patients compared to controls. The effect of sex should be further investigated on the individual level through large studies or mega-analyses. Although there was no difference in glutamate metabolite variability in the temporal lobe, higher glutamate levels were found in studies which included more highly symptomatic patients. This is also consistent with our recent mega-analysis(12). Finally, other than in the basal ganglia and frontal white matter, meta-regressions did not find an accelerated loss of glutamatergic metabolites in patients with age for the majority of brain regions, also consistent with our recent mega-analysis(12).
A limitation of the meta-analysis is the high between-study inconsistency, as measured by I2, for most brain regions studied. This was highest in the temporal lobe and may relate to the difficulty of obtaining good quality 1H-MRS imaging in this region. There is a possibility that case-control differences in variability result from greater movement artifacts in patient populations(68). A recent meta-analysis emphasises the importance of using strict Cramér–Rao lower bound criteria (≤ 7%) and short echo times (≤ 20 ms) to improve 1H-MRS consistency(5). Furthermore, the glutamine signal cannot be accurately resolved from glutamate below 3T, although the majority of studies reporting glutamine were conducted above 3T. As voxel placement varied between studies, broad categories of brain regions were used, limiting the regional specificity of our results. Meta-regression analyses of clinical and demographic variables are limited to the study level and are not sensitive to variation within individual studies (although meta-regressions with the SD of clinical and demographic variables were carried out). Lastly, the number of included studies is low for some brain regions, such as the thalamus, and there are a small number of studies examining antipsychotic-naïve patients in all regions except the MFC, and so these sensitivity analyses should be considered preliminary.
In summary, this meta-analysis demonstrates increased regional variability in glutamatergic metabolites in schizophrenia in addition to mean differences compared to controls. Increased inter-individual differences in glutamatergic metabolites in schizophrenia are likely to have a complex mechanistic basis. Further work is also required to determine the clinical consequences along the spectrum of glutamate dysregulation. Both glutamatergic metabolite levels(12) and interindividual variability appear to be greater in younger and more symptomatic patients. Neurobiological heterogeneity may also relate to heterogeneity in antipsychotic response in schizophrenia, and some studies have shown that glutamatergic metabolite levels in the MFC, thalamus, DLPFC and striatum associate with the degree of antipsychotic response(15, 16, 67, 69, 70, 17–24). Our findings are relevant to the on-going effort to develop novel drug therapies to target glutamate dysfunction in schizophrenia, as the presence of glutamatergic heterogeneity may indicate the importance of targeting more specific patient subgroups.