Abstract
Introduction
Approximately a third of patients with schizophrenia do not respond to standard antipsychotic medications and are hence categorised as having treatment-resistant schizophrenia (TRS) (Howes et al., 2017). Clozapine is considered the gold standard treatment for TRS due to its superior efficacy in this patient population (The Lancet Psychiatry, 2024). However, the neurobiological basis of TRS and the mechanism of action of clozapine remain unclear.
Schizophrenia is associated with alterations in both brain glutamate concentrations (Egerton et al., 2020) and cerebral blood flow (CBF; Goozée et al., 2014). Recent systematic reviews and meta-analyses have indicated higher Glx (glutamate plus glutamine) in the basal ganglia across the schizophrenia-spectrum disorders (Merritt et al., 2016, 2023; Nakahara et al., 2022), lower medial frontal cortex glutamatergic levels with individual differences in the degree of glutamate dysfunction (Merritt et al., 2023), and cortical hypoperfusion and subcortical hyperperfusion (Percie du Sert et al., 2023). Neurobiological alterations in brain glutamate and blood flow might be more prominent in TRS. For example, in the meta-analysis by Nakahara et al. (2022), while there were no significant alterations in frontal cortical glutamate metabolites overall in patients compared to controls, the TRS group exhibited higher glutamate and Glx in the mid-cingulate cortex, though this effect was not seen in the basal ganglia. We have also observed widespread reductions in CBF in TRS compared to control, across cortical (anterior cingulate cortex (ACC)) and subcortical (hippocampus, thalamus, & striatum) regions (Sun et al., 2024).
Alterations in glutamate concentration and local CBF are thought to be functionally related via neurovascular coupling. As the primary excitatory neurotransmitter, glutamate plays a crucial role in regulating local CBF to meet metabolic demands (Attwell et al., 2010). Glutamate release and activation of N-methyl-D-aspartate (NMDA) receptors on astrocytes evoke increases in intracellular Ca2+ and release of vasoactive substances to cause local vasodilation (Attwell et al., 2010; Fergus and Lee, 1997; Zhang et al., 2024). While there are few human brain
Emerging evidence suggests that antipsychotic treatment dynamically affects the relationships between neurotransmitter activity and CBF. For example, risperidone treatment reversed the direction of the association between ACC Glx and the blood oxygen level-dependent (BOLD) response in schizophrenia after 6 weeks (Cadena et al., 2018). There is separate evidence that glutamate and CBF are modulated by antipsychotic treatment and that, potentially, this may contribute to symptomatic response. Longitudinal studies have suggested a pattern of decreasing glutamate levels during antipsychotic administration in schizophrenia, particularly in frontal regions (Egerton et al., 2017; Kubota et al., 2020). Decreases in glutamate following antipsychotic treatment are also supported by a recent meta-analysis, which shows a negative association between antipsychotic dose and glutamate or Glx concentration in the medial frontal cortex, including the ACC (Merritt et al., 2021). A separate meta-analysis indicated higher Glx in the dorsolateral prefrontal cortex and hippocampus only in unmedicated patients compared to healthy controls (Nakahara et al., 2022).
The efficacy of clozapine in TRS may also be partially related to glutamatergic effects, including positive modulation of NMDA receptor-mediated glutamate transmission (Arvanov et al., 1997; Arvanov and Wang, 1999; Chen and Yang, 2002; Gemperle et al., 2003; Gray et al., 2009; Homayoun and Moghaddam, 2007; Malhotra et al., 1997), suppression of hyperglutamatergic activity via activation of presynaptic group metabotropic glutamate receptors (e.g. III-mGluR) (Fukuyama et al., 2019), and attenuating MK-801-induced increases in glutamate release in the frontal cortex (López-Gil et al., 2007). We previously detected a decrease in striatal, but not ACC, glutamate levels during clozapine treatment in schizophrenia, which correlated with improvement in symptom severity (McQueen et al., 2021).
The changes in CBF observed during antipsychotic treatment appear to be regionally and potentially treatment-specific (Goozée et al., 2014). For example, haloperidol, a prototypical antipsychotic, decreased CBF in frontal and cerebellar regions and increased CBF in the basal ganglia (Lahti et al., 2003; Miller et al., 2001); while risperidone, an atypical antipsychotic, decreased CBF only in the cerebellar region (Miller et al., 2001). With regard to clozapine, mixed findings have been reported, including CBF decreases (Sun et al., 2024), increases (Ertugrul et al., 2009; Lahti et al., 2003; Sharafi, 2005), and no change in cortical regions (Zhao et al., 2006). Using the magnetic resonance imaging (MRI) approach of pseudo-continuous arterial spin labelling (pCASL), we recently reported decreases in CBF in the ACC over 12 weeks of clozapine treatment, consistent with some previous radiotracer studies (Lahti et al., 2003; Zhao et al., 2006), and that higher striatal and hippocampal CBF prior to commencing clozapine was associated with greater subsequent improvement in symptom severity (Sun et al., 2024).
The aim of this study was to examine the local relationships between glutamate metabolites and CBF in TRS and their change over 12 weeks of clozapine treatment. We hypothesised that, within the ACC and striatum, glutamate and CBF would be positively correlated, and that their changes over 12 weeks of clozapine treatment would also show a positive correlation. Exploratory analyses examined associations between glutamate, CBF, and symptom improvement. In addition, as glutamate may influence neural activity and thereby potentially CBF in distal brain regions (Duncan et al., 2014; Kiemes et al., 2021; Zahid et al., 2023), we also examined the relationships between ACC or striatal glutamate metabolites and CBF using a whole-brain voxel-wise approach.
Methods
Study design and participants
Ethical approval for this study was obtained from the London South East National Health Service Research Ethics Committee (Reference 13/LO/1857), and informed consent to participate was obtained from participants with the mental capacity to do so. For participants without the capacity to provide consent, a carer provided advice on their behalf.
Participants were recruited from inpatient and outpatient services at the South London and Maudsley and Oxleas NHS Foundation Trusts. Inclusion required that participants were ≥18 years old, met ICD-10 criteria for schizophrenia (F20) or schizoaffective disorder (F25), and were either clozapine naïve or had not taken clozapine for ≥3 months prior to the consent visit. Criteria for TRS were inferred from medical records and the treating psychiatrist, including having ≥2 previous trials of non-clozapine antipsychotics, each within the recommended dose for ≥6 weeks, as according to clinical judgement, and were being referred for clozapine initiation. Exclusion criteria for all participants included substance dependency as defined using ICD-10 criteria, pregnancy, or contraindication to MRI. The observation period was from baseline (after referral for clozapine but prior to clozapine initiation/reinitiation) to 12 weeks after clozapine initiation. The 12-week period was selected based on consensus guidelines, recommended for evaluating clozapine response (Howes et al., 2017).
Clinical and demographic measures
Demographic variables included age and sex. Medical history was obtained through clinical interview and review of medical records and included age of onset (also used to calculate duration of illness), diagnosis, number of hospital admissions, and previous antipsychotic trials, name and dose of current antipsychotic prior to clozapine initiation, and any co-medications. At the time of each MRI visit, symptom severity was assessed using the Positive and Negative Syndrome Scale (PANSS; Kay et al., 1987), and social and occupational functioning was assessed with the Global Assessment of Functioning Scale (GAF) (Hall, 1995). Blood samples were taken to monitor clozapine plasma levels, to determine whether clozapine levels met the indicated therapeutic threshold of 350 ng/ml (Siskind et al., 2021) at week 6, week 12, and additionally at week 8 (if ≤350 ng/ml at week 6) (Rosenheck et al., 1999). Change in total symptom severity was calculated as the percentage change in PANSS total from baseline to week 12, calculated using the formula (Leucht et al., 2007):
1H-MRS Acquisition and Quantification
All data were acquired on a GE Healthcare MR750 3T scanner (Waukesha, USA) at the Centre for Neuroimaging Sciences, Institute of Psychiatry, Psychology & Neuroscience (King’s College London), across two sessions (‘baseline’, ‘week 12’). Acquisition of parameters for structural (T1-weighted) data was acquired using a 3D T1-weighted inversion recovery prepared spoiled gradient recalled sequence (Inversion time (TI) = 400 mseconds, Repetition time (TR) = 7.31 mseconds, Echo time (TE) = 3.02 mseconds, voxel size = .05 × 1.05 × 1.2 mm3, Field of view (FoV) = 270 mm, flip angle = 11°).
The proton magnetic resonance spectroscopy (1H-MRS) data was handled and described in the same manner, including parameters for data acquisition, analysis, and quality control, as in the study with an overlapping sample (McQueen et al., 2021). 1H-MRS glutamate and Glx levels were obtained from the conventional Point-Resolved Spectroscopy (PRESS) sequence, incorporating a CHEmically Selective Suppression (CHESS) water suppression routine (TE = 30 mseconds in the ACC and TE = 35 mseconds in the striatum, TR = 3000 mseconds, 96 averages). For each acquisition, unsuppressed water reference spectra (16 averages) were also acquired for eddy current correction and water scaling. The region of interest (ROI) was prescribed from the T1-weighted image, to be: (1) in the centre of the ACC (20 × 20 × 20 mm) (Egerton et al., 2012, 2018); (2) in the right striatum (20 × 20 × 20 mm, mainly covering the right head of striatum with the lower end of the voxel located dorsal to the anterior commissure, to include the maximum amount of grey matter and minimum amount of cerebrospinal fluid (CSF) (de la Fuente-Sandoval et al., 2013; McQueen et al., 2018). Images of voxel placement are provided in Figure S1.
Spectra were analysed using Linear Combination of Model Spectra (LCModel) version 6.3-0I (Provencher, 1993), using LCModel basis sets matched to the TEs for each voxel. Values for glutamate metabolites were corrected for voxel tissue composition by using the formula:
where WM, GM, and CSF indicate the fraction of white and grey matter and cerebrospinal fluid content in the voxel. This formula assumes a CSF water concentration of 55,556 mol/m3 (Gasparovic et al., 2006). Voxel GM ratio was calculated as
pCASL Acquisition and Quantification
The data was handled and described in the same manner, including parameters for data acquisition, computation of CBF maps, and procedures for spatial normalisation of these maps to the Montreal Neurological Institute (MNI) reference space, as in the study with an overlapping sample (Sun et al., 2024). pCASL images were acquired using a 3D Fast Spin Echo spiral multi-shot readout. The parameters of the arterial labelling part of the sequence were: label duration = 1800 mseconds, post labelling delay = 2025 mseconds, TR = 4968 mseconds. For the readout part, 3D images were acquired with a Fast Spin Echo ‘stack of spirals’ routine, using 8 spiral arms per slice location, 512 points per arm, 55 slice locations of 3 mm thickness, TE = 11.09 mseconds, FoV = 240 mm. A proton-density image with the same parameters was collected within the same sequence, which was used to derive a whole brain CBF map in physiological units (ml blood/100 g tissue/minute) following the guidelines recently published for CBF computation (Alsop et al., 2015). The total acquisition time for the complete ASL pulse sequence and the proton-density image was in 6:28 minutes.
CBF maps were pre-processed and spatially normalised to the standard space of the MNI using Automatic Software for ASL Processing (ASAP) toolbox 26 (Mato Abad et al., 2016) running in Statistical Parametric Mapping Software (SPM-12, Functional Imaging Laboratory, Wellcome Centre for Human Neuroimaging, Institute of Neurology, University College London, UK, http://www.fil.ion.ucl.ac.uk/spm/software/spm12/). These spatially normalised images were used for ROI and voxel-wise analysis. For primary analyses, ROI masks in the bilateral ACC and striatum were defined using WFU PickAtlas in SPM-12 (‘ACC’ and ‘striatum’). For confirmatory analyses, as 1H-MRS data were acquired in the right striatum only, a separate ROI mask for the right striatum only was generated from the bilateral mask (anatomical ROI right striatum), and, to extract CBF values from the 1H-MRS voxel specifically, two additional masks were generated in SPM-12 from the ACC and right striatal 1H-MRS voxels. Mean CBF values were extracted from each mask using the ASAP toolbox. Mean global CBF values were extracted from the group grey matter mask to investigate global differences and to account for inter-individual differences in global perfusion.
Statistical Analyses
Clinical data
The within-subject change in symptom severity and functioning following clozapine treatment was determined using paired t-tests with a two-tailed alpha of 0.05, in Python version 3.12.6.
Integration of 1H-MRS and pCASL data
ROI analysis
Shapiro-Wilk tests checked CBF and 1H-MRS variables for normality of distribution. Correlations between variables were analysed using Pearson’s tests. Associations between ACC or striatal glutamate, Glx and CBF within the same regions were determined at baseline, week 12, as well as their changes over 12 weeks. We gave a
Exploratory analyses of the relationships between glutamate concentrations, regional CBF, and symptoms were conducted using hierarchical linear regression. Separate hierarchical models examined the relationships between (a) change in glutamate (or Glx) and CBF and improvement in total symptom severity, and (b) the baseline glutamate (or Glx) and CBF and subsequent improvement in symptom severity. For analysis (a), change in global CBF was entered in step 1 (Model 1, control variable); change in glutamate metabolites and change in CBF values as predictors in step 2, covarying for change in global CBF (Model 2); and the interaction term change in glutamate metabolites × CBF values was entered in step 3 (Model 3). Similarly, for analysis (b), baseline global CBF was entered in step 1 (Model 1, control variable); baseline glutamate metabolites and baseline CBF values as predictors in step 2, controlling for baseline global CBF (Model 2); and the interaction term baseline glutamate metabolites x CBF values was entered in step 3 (Model 3). In addition, the week 12 clozapine dose was entered as a covariate (Model 4) in both analyses. The threshold for statistical significance was
Whole-brain voxel-wise analyses
Supplementary voxel-wise analysis was conducted in SPM-12. To examine relationships between the change in glutamate (or Glx) and the change in CBF during clozapine treatment, using ImCalc, we first calculated one individual CBF difference image in MNI space for each participant (follow-up image – baseline image), representing the change in CBF during clozapine treatment. Changes in individual glutamate levels were then entered as regressors in a voxel-wise ANOVA, with or without week 12 clozapine dose entered as a covariate. To examine glutamate (or Glx)-CBF associations at the baseline and week 12 timepoints separately, individual participant glutamate concentration values were entered as regressors in a voxel-wise ANOVA, with or without global CBF as a covariate of no interest. The threshold for statistical significance was defined as
Results
Both 1H-MRS and pCASL data were available for 30 individuals with TRS before clozapine initiation (baseline). Of these, data were available for 20 individuals at follow-up, 12 weeks after clozapine treatment. This study integrated a subsample of participants from previous publications reporting 1H-MRS (McQueen et al., 2021) and CBF (Sun et al., 2024) separately. Characteristics of the participants are shown in Table 1. Six of the 30 TRS subjects had previously taken clozapine. Symptom severity and functioning improved over 12 weeks of clozapine treatment. No significant demographic or clinical differences existed between those who did or did not complete the 12-week follow-up visit (Table S1). As previously reported (McQueen et al., 2021; Sun et al., 2024), ACC CBF (
Descriptive statistics of the patient sample at baseline and follow-up.
Data is presented as mean (standard deviation) unless otherwise specified.
ACC: Anterior cingulate cortex; CBF: Cerebral blood flow; CPZE: Chlorpromazine equivalent dose; GAF: Global assessment of functioning; Glx: Glutamate + Glutamine; GABA: γ-aminobutyric acid; GABAergic medication including benzodiazepines, valproate, lithium, zopiclone, pregabalin, lamotrigine; PANSS: Positive and Negative Syndrome Scale.
Relationships between CBF and glutamate metabolites in the ACC and striatum
At the baseline (prior to clozapine), there were no significant correlations between CBF and glutamate (or Glx) within the ACC or striatum (all
Pearson’s correlations between CBF and glutamate metabolites in the ACC and striatum.
ACC: Anterior cingulate cortex; CBF: Cerebral blood flow; Glu: Glutamate; Glx: Glutamate + Glutamine; MRS: Magnetic resonance spectroscopy.

Scatterplot of change in glutamate and change in CBF during clozapine treatment in TRS (
Relationships with clinical variables
To investigate how CBF and glutamate metabolites may change in relation to clinical symptoms, we performed hierarchical multiple regression analyses. In these analyses, we first assessed regional CBF and glutamate or Glx in relation to total symptom improvement via multiple regression while considering the effect of global CBF. We then assessed the same relationship with interaction via moderation analyses. Finally, the clozapine dose at week 12 was included in an additional model.
The change over time in glutamate (or Glx) and CBF in the ACC and striatum did not significantly contribute to explaining variance in total symptom improvement during clozapine treatment (all
ACC CBF and ACC Glx prior to commencing clozapine significantly contributed to explaining variance in subsequent total symptom improvement during clozapine treatment (controlling for baseline global CBF: overall model:
Voxel-wise analyses
Voxel-wise analysis additionally detected a positive relationship between change in striatal glutamate and change in CBF in the parietal cortex during clozapine treatment (one cluster: MNI (

Relationship between change in striatal glutamate and change in cortical CBF during clozapine (

Relationship between ACC glutamate and cerebellar CBF after 12 weeks of clozapine (
Discussion
This study examined the relationships between changes in brain glutamate and CBF over 12 weeks of clozapine treatment in TRS. Our main finding was that longitudinal changes in striatal glutamate and blood flow following clozapine treatment were positively correlated. Second, patients with higher Glx and CBF in the ACC prior to treatment tended to have better clinical improvement in total symptom severity during subsequent clozapine treatment, with this effect being significantly driven by ACC Glx, although these effects fell below the threshold for significance when the clozapine dose at week 12 was added to the model.
A positive association between reductions in striatal glutamate and striatal CBF during clozapine treatment is in keeping with a mechanistic association between neural glutamatergic activity and local CBF to support metabolism (Attwell et al., 2010), and indicates that reductions in striatal glutamate during clozapine treatment (McQueen et al., 2021) were correlated with functional consequences. According to established theories (Attwell et al., 2010), glutamate release increases intracellular Ca2+, which then triggers the release of vasoactive molecules, dilating blood vessel diameter (Filosa et al., 2004; Mishra et al., 2016; Zonta et al., 2003). Neurovascular coupling predicts that an increase in inhibitory activity will result in a decreased/negative BOLD functional MRI response (Donahue et al., 2010; Krause et al., 2014; Logothetis, 2008). Within our study, while it should be noted that 1H-MRS does not measure glutamate involved in neurotransmission specifically, the decrease in glutamate levels during clozapine treatment may indicate reduced excitatory activity, which could, in turn, lower metabolic demand and local CBF. Clozapine may act to normalise glutamatergic synaptic function (Hribkova et al., 2022), enhance NMDA receptor activity while blocking NMDA antagonist effects, and inhibit glutamate transmission (Fukuyama et al., 2019; Hribkova et al., 2022). Consequently, a decrease in striatal glutamate-mediated signalling may reduce neurovascular responses and decrease local blood flow (Shaw et al., 2021; Sun et al., 2013; Takano et al., 2006).
It is also possible that the observed positive relationship between the reduction in glutamate and the reduction in CBF in the striatum during clozapine treatment may reflect adaptive mechanisms related to discontinuing non-clozapine antipsychotics, either in addition to or instead of reflecting the effects of clozapine treatment specifically. For example, treatment with non-clozapine antipsychotics appears to be associated with increases in striatal blood flow (Goozée et al., 2014) and grey matter volume (Ebdrup et al., 2013; Jørgensen et al., 2016; Yang et al., 2021), whereas decreases in striatal volume are observed during clozapine treatment (Krajner et al., 2022; Sun et al., 2024). The relative contributions of non-clozapine antipsychotic discontinuation and clozapine treatment to striatal glutamate and CBF changes could be further explored in pre-clinical models. It is also possible that structural changes in the striatum could affect both striatal glutamate and CBF, although previous findings in this cohort suggest, changes in CBF and glutamate levels seem independent of the decrease in grey matter volume (Krajner et al., 2022; Sun et al., 2024).
In contrast to the striatum, although CBF was decreased in the ACC during clozapine treatment, this was not significantly associated with a change in ACC glutamate metabolites. This may indicate a lesser extent of neurovascular coupling in the ACC, and/or that the decrease in CBF in this region reflects non-glutamatergic processes, for example, changes in GABAergic activity, neuroinflammation, or vascular tone. Interestingly, a recent bimodal meta-analysis of regional CBF and glucose metabolism in schizophrenia indicates conjoint findings in the bilateral striatum, but disjoint findings, indicating neurovascular uncoupling, in dorsomedial frontal areas partially overlapping with the ACC region in our analysis (Sukumar et al., 2020).
Previous cross-sectional studies detected a positive relationship between glutamate and white matter CBF in individuals with schizophrenia (Wijtenburg et al., 2017) and between glutamate and grey matter CBF in healthy controls (Bojesen et al., 2018) in the ACC. We were unable to investigate the relationship between glutamate and white matter CBF due to the low CBF signal in white matter in our data. We did not observe any significant relationship between ACC glutamate and ACC grey matter CBF in participants with TRS prior to starting clozapine treatment, which may reflect illness or medication effects. Furthermore, ACC Glx predicted subsequent total symptom improvement during clozapine treatment, whereas ACC CBF did not contribute. This effect became non-significant when additional variance in symptom improvement was accounted for by including clozapine dose in the model. This exploratory finding warrants investigation in a larger, independent sample but may tentatively suggest that neurochemical signals can better predict clozapine response than blood flow.
The findings of this study should be considered in the context of its limitations, including the relatively small sample size of 20 participants at both time points. With 80% power, alpha level of 0.05, this sample size provides sensitivity to detect a correlation of approximately
Conclusions
In conclusion, our study provides evidence of a positive relationship between changes in glutamate and CBF in the striatum after switching from a non-clozapine antipsychotic to clozapine treatment. This indicates that glutamatergic changes, as observed with 1H-MRS, are accompanied by changes in brain activity. In contrast, this association was not observed in the ACC, which may reflect regional differences in neurovascular coupling in schizophrenia.
Supplemental Material
sj-docx-1-jop-10.1177_02698811251409141 – Supplemental material for Associations between glutamate and cerebral blood flow in treatment-resistant schizophrenia during clozapine treatment
Supplemental material, sj-docx-1-jop-10.1177_02698811251409141 for Associations between glutamate and cerebral blood flow in treatment-resistant schizophrenia during clozapine treatment by Junyu Sun, Fernando Zelaya, Kyra-Verena Sendt, Grant McQueen, Amy L. Gillespie, John Lally, Owen O’Daly, David J. Lythgoe, Oliver D. Howes, Gareth J. Barker, Philip McGuire, James H. MacCabe and Alice Egerton in Journal of Psychopharmacology
Footnotes
Funding
Declaration of conflicting interests
Supplemental material
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
