ABSTRACT
Background Evidence suggests dysfunction of GABAergic interneurons in psychosis, and prior research has linked GABAergic function with a tendency toward negative affective states. Magnetic resonance spectroscopy (MRS) studies measuring GABA have yielded inconsistent findings. We investigate GABA concentrations in young adults with attenuated psychosis syndrome (APS) and first episode psychosis (FEP), as well as testing the hypothesis that negative affect is a clinical phenotype that is associated with reduced GABA.
Materials MRS data were obtained from 14 patients with FEP, 7 patients with APS and 15 healthy controls (HC), using a MEGA-PRESS sequence on a 3T Philips Ingenia scanner. Voxels were placed in rostral MFC and midline-occipital cortex. Gannet 3.1 was used to determine GABA+ and Glx (glutamate and glutamine combined) concentrations.
Results We found a trend towards increased rostral MFC GABA+ concentrations in FEP, but no group differences in occipital GABA+ concentrations. When covarying for scores on the Psychological Stress Index, rostral MFC GABA+ levels in FEP were significantly greater than APS and HC. Planned comparisons revealed a trend towards increased rostral MFC GABA+ in APS relative to HC. No group differences in Glx or occipital GABA+ were found.
Conclusion These results, considered alongside previously published findings, suggest multiple factors influencing GABA+ in psychosis. We conclude a process exists which drives up GABA+ in early psychosis, alongside a separate process in which reduced GABA+ is associated with increased negative affect. These multiple processes have resulted in contradictory findings, and their untangling is critical to understanding of GABA+ in psychosis.
1. INTRODUCTION
Considerable evidence implicates dysfunction of GABAergic systems in schizophrenia, alongside other psychosis spectrum disorders. Post-mortem findings of reductions in the expression of glutamic acid decarboxylase-67 (GAD67), an enzyme critical for GABA synthesis, have been well replicated in both schizophrenia (Curley et al., 2011; Guidotti et al., 2000; Lewis et al., 2012; Schmidt & Mirnics, 2015) and bipolar disorder (Guidotti et al., 2000; Torrey et al., 2005; Woo et al., 2008), and observed in both the cerebral cortex and hippocampus. GAD67 reductions also have been found in rodent models of schizophrenia (Behrens et al., 2007; Liu et al., 2001; Takao et al., 2013). Using a pharmacological challenge with lorazepam, a non-subtype selective benzodiazepine that potentiates GABA function, in combination with blood oxygen-level dependent (BOLD) fMRI, our group has demonstrated GABAergic dysfunction in schizophrenia both while patients passively viewed emotionally salient images (Taylor et al., 2014) and while performing a face processing task (Tso et al., 2015). Further, the altered BOLD response to lorazepam in schizophrenia was positively correlated with increased negative affect.
While traditionally there has been a focus on the association of GABA with cognitive deficits, emerging links between GABAergic systems, affect, and emotion processing indicate an important line of inquiry. Affective dysregulation is often exhibited by patients with psychosis spectrum disorders (Myin-Germeys & van Os, 2007; Phillips et al., 2007). We conceptualise this dysregulation as a high propensity to experience negative affect states such as anxiety and dysthymia, and to exhibit a low tolerance for even minor, everyday stressors. Patients with schizophrenia report more negative experiences in daily life (Docherty, 1996; Horan et al., 2005; Jones & Fernyhough, 2007; Myin-Germeys et al., 2000) and perceive emotional stimuli with a more negative bias in laboratory settings (Burbridge & Barch, 2002; A. S. Cohen & Minor, 2010; Holt et al., 2006). Importantly, measures of negative affect have significant prognostic value in schizophrenia, predicting hospitalisation and poor functional outcome (Blanchard et al., 1998; Grove et al., 2016; Horan et al., 2008).
Here, we aimed to build upon findings from our group, which demonstrated negative affect associated GABAergic abnormalities in schizophrenia using in-vivo functional imaging, by exploring these questions using magnetic resonance spectroscopy. Thus far, studies of GABA concentrations in schizophrenia using magnetic resonance spectroscopy have reported mixed findings, with evidence for increased (de la Fuente-Sandoval et al., 2018; Öngür et al., 2010; Yang et al., 2015), decreased (Goto et al., 2009; Menschikov et al., 2016; Thakkar et al., 2017; Yoon et al., 2010) and as normal GABA concentrations when compared with healthy controls (C. M. A. Chen et al., 2014; Tayoshi et al., 2010). These divergent findings likely reflect the influence or multiple factors or processes upon GABAergic systems, such as illness stage and symptom profile, as well as the location of the MRS voxel in the brain. Indeed, when examining rostral medial frontal cortex (MFC) GABA concentrations, several groups studying early, unmedicated patients have revealed elevated GABA concentrations (T. Chen et al., 2017; de la Fuente-Sandoval et al., 2018; Kegeles et al., 2012), although not all groups have replicated this finding (e.g. Wang et al., 2016). By contrast, studies of older medicated patients have revealed no abnormalities (Chiu et al., 2018; Kegeles et al., 2012; Rowland et al., 2013), or reductions when compared with healthy controls (Rowland et al., 2013; Thakkar et al., 2017).
In this study, we used MRS to compare GABA+ concentrations in the rostral MFC and midline-occipital cortex in patients with first episode psychosis (FEP) and individuals meeting criteria for attenuated psychosis syndrome (APS) with healthy, age-matched controls. Additionally, we explore links between GABA+ concentrations and levels of negative affect (NA) as measured by the Psychological Stress Index (PSI9: Tso et al., 2012). Based on prior reports of increased rostral MFC GABA in early schizophrenia, we hypothesized that GABA+ concentrations in this voxel would be increased relative to controls. We also predicted, based on findings linking increased NA with altered BOLD responses, that GABA+ concentrations would be associated with NA. As glutamate has been linked with emotion processing in healthy controls (Stan et al., 2014), as well as reported abnormalities in the excitation/inhibition balance in psychosis (for a review see Liu et al., 2021), we also performed an exploratory analysis examining relationships between NA and Glx (the combined resonance of glutamate and glutamine).
2. METHODS
2.1 Participants
Study participants were 14 patients with FEP, 7 individuals meeting APS criteria, and 15 healthy controls (HC) who were group matched to the psychosis groups in terms of age, gender, and years of parental education (Table 1). Participants were recruited from the local community via advertisements and referrals from clinicians at local mental health clinics. The Structured Interview for Psychosis Syndromes (SIPS: Miller et al., 2003) and Structured Clinical Interview for DSM-IV Axis I Disorders (SCID-IV: First et al., 1996) were used to determine APS status and psychiatric diagnoses. Participants completed the 9-item PSI9 (Tso et al., 2012), a validated self-report scale which was modelled on Perceived Stress Scale (S. Cohen et al., 1983). The PSI9 was designed to measure negative affect and stress sensitivity over the last month, but can also capture trait-like tendencies, and has demonstrated predictive power for both functional and clinical outcomes (Tso et al., 2012).
All participants were aged 15 – 27 years old, with normal, or corrected to normal vision. Exclusion criteria included contraindications to MRI (e.g., pacemakers, metallic implants), alcohol or substance use in the past month or dependence in the last 6 months, history of head injury or concussions with unconsciousness lasting > 5 minutes. Healthy control participants were free from histories of psychiatric disorders, and history of psychotic or bipolar disorder in their first-degree relatives.
Study procedures were approved by the Institutional Review Board of the University of Michigan Medical School. Following full explanation of the study, written informed consent was obtained from participants aged 18 years old and over. For those under 18 years, written parental consent was obtained, alongside written assent from the participant.
2.2 MR acquisition
MR data were acquired using a 3T Phillips Ingenia system with a 32-channel head coil. Head motion was minimized using padding and Velcro straps. For anatomical reference and tissue classification, we obtained a structural scan for each participant (3D-GRE, slice thickness = 1.4 mm, TR = 25 ms, TE = minimum, flip angle = 35, field of view = 24 cm).
For the assessment of GABA+ and Glx, MRS data were collected using a MEGA-PRESS sequence (Mescher et al., 1998), with the following parameters: TE = 68 ms (TE1 = 15 ms, TE2 = 53 ms); TR = 1.8 seconds; 256 transients of 2000 data points; spectral width 2 kHz; frequency selective editing pulses (14 ms) applied at 1.9 ppm (ON) and 7.46 ppm (OFF).
Data were collected from two voxels (Figure 1): a 34 × 36 × 20mm rostral MFC voxel placed anterior to the genu of the corpus collosum, and a 35 × 25 × 25mm occipital cortex voxel centered on the calcarine sulcus, parallel to the cerebellar tentorium. Acquisition times were 11 mins 24 sec for the MFC voxel and 11 mins 12 sec for the occipital voxel.
Recommended minimum reporting details for the MRS data are included in Supplemental Table 3, set out as per Minimum Reporting Standard for in-vivo Magnetic Resonance Spectroscopy (MRSinMRS): Experts’ Consensus Recommendations (Lin et al., 2021).
2.3 Data Preprocessing and Analysis
Preprocessing and quantification of the edited MRS spectra were performed using Gannet 3.1 (Edden et al., 2014), an open-source MATLAB based toolbox. Standard processing steps including frequency and phase corrections of time-resolved data using spectral registration were used. GABA and Glx levels at 3.00 and 3.75 ppm respectively were estimated using nonlinear, least-squares fitting. As the GABA signal detected at 3.00 pm using our experimental parameters is expected to contain contributions from macromolecules and homocarnosine, it is referred to as GABA+ throughout this manuscript.
GABA+ and Glx concentrations were evaluated relative to the unsuppressed water signal. While there is debate regarding the best choice of reference metabolite (Mullins et al., 2014), evidence suggests that Cr is decreased in schizophrenia (Öngür et al., 2009), thus we opted against this as our primary measure. Nevertheless, many studies use Cr as a baseline, and we report GABA+/Cr and Glx/Cr in our supplementary materials, to allow comparison with these reports.
Since psychosis, even in its early stages, is associated with morphological changes in the brain such as frontal and temporo-limbic reductions in grey matter (Takahashi & Suzuki, 2018), we calculated estimates of GABA+ and Glx that are corrected for the tissue composition of the MRS voxel. Specifically, Gannet 3.1 coregisters the MRS voxels on to the T1-weighted structural images which we acquired during the same scanning session, then uses SPM12 segmentation functions to calculate the fractions of grey matter (GM), white matter (WM) and cerebrospinal fluid (CSF) in each voxel. As CSF contains negligible amounts of GABA or Glx, Gannet 3.1 uses these fractions to remove the impact of CSF fraction on the voxel, as well as applying tissue specific relaxation and water visibility constraints for each of the three tissue types (Harris et al., 2015). We present the alpha-correction in the main manuscript as our primary measure, which assumes the concentration of GABA+ and Glx in GM is twice that of WM (Ganji et al., 2014) and applies an alpha of 0.5 to correct for between-subject differences in voxel tissue composition. Uncorrected measures referenced to water and creatine included in the supplementary materials.
Spectra were visually inspected for artefacts across Gannet 3.1’s preprocessing steps and as well as goodness of model fit. Gannet 3.1 provides a measure of fit error and to ensure the robustness of our findings, only estimates of GABA+ and Glx with a fit error of less than 15% were included in our analysis. Gannet 3.1 also provides spectra quality metrics including signal to noise ratio (SNR) and linewidth (FWHM) which were also assessed, and spectra with values greater than 3 SD of the sample mean for these metrics were excluded. Any extreme GABA+ values of greater than 5 SD of the sample mean were excluded, to prevent these values having a high influence on statistical analyses.
2.4 Statistical analyses
All statistical analyses were performed in R. Group comparisons of demographic and clinical measures were performed using ANOVA, or chi-square tests where measures were categorical. To investigate group differences in GABA+ and Glx, we performed ANOVA, followed by ANCOVA which included PSI9 scores as a covariate. Since our primary measure of interest was GABA+ concentrations in the two voxels, we controlled for multiple comparisons by setting our significance threshold at 0.025 (representing a Bonferroni correction of 0.05/5). Following significant main effects of group, planned comparisons were used to compare GABA+ concentration differences between the FEP group and the APS and HC participants (FEP vs HC + APS), and between APS and HC.
3. RESULTS
3.1 Sample characteristics
Summaries of demographic and clinical measures are presented in Table 1. Groups did not differ on age, gender, handedness, parental education, or socioeconomic status. APS and FEP groups did not differ on PANSS positive scores, however the FEP group had higher PANSS total scores (t17.11 = -4.14, p<.001), as well as higher PANSS negative (t16.60 = 3.31, p = .004) and PANSS general psychopathology (t18.92 = -3.16, p = .005) scores. There was a significant group difference in PSI9 scores (F2,33, = 12.14, p<.001), with post-hoc Tukey tests showing PSI9 scores differed significantly between the FEP and HC groups (p < .0001), but not between the APS and HC groups (p = .08) or between the FEP and APS groups (p = .21).
Figure 2 presents the degree of spatial overlap for both voxels, alongside their edited spectra. In the rostral MFC voxel, GABA+ data was excluded from three FEP participants and one HC due to GABA+ fit errors of greater than 15%. Extreme GABA+ values (greater than 5 SDs from the mean) led to the exclusion of rostral MFC GABA+ and Glx from one FEP participant, and Occ GABA+ and Glx from one APS participant.
3.2 Voxel tissue composition and spectral quality
Table 2 includes a breakdown of voxel tissue compositions. No group tissue composition differences were found in the rostral MFC voxel (GM: F2,32 = .92, p = .41; WM: F2,32 = .96, p = .39; CSF: F2,32 = 0.20, p = .82). In the occipital voxel, there was a significant main effect of group on GM fraction (F2,32 = 4.16, p = .02), with significantly greater GM fraction in HC participants compared with FEP (p =.01). Group differences in WM fraction trended towards significance (F2,32 = 2.59, p = .09). There were no group differences in CSF fraction (F2,32 = .40, p = .68).
GABA+ and Glx fit errors are presented in Table 3. There were no group differences in GABA+ fit error % (MFC: F2,28 = 1.57, p = .27; Occ: F2,32 = 0.22, p = .81) or Glx fit error % (MFC: F2,32 = .71, p = .50; Occ: F2,32 = 0.34, p = .72). MRS spectral quality measures are presented in Supplemental Table 1. There were no group differences in linewidth (FWHM) or SNR.
3.3 Group Differences in GABA+ and Glx
Group average metabolite concentrations are presented in Table 3 and Figure 3. When we conducted ANOVAs to examine group differences in GABA+, we found a trend towards a significant effect of group in the MFC (F2,28 = 3.04, p = .06), but no differences in the occipital cortex (F2,32 = 0.49, p = .62). There were no main effects of group on Glx concentration in either voxel (MFC: F2,28 = 0.08, p = .93; Occ: F2,32 = 0.34, p = .72).
When PSI9 was included in the model to adjust for group differences in negative affect, we found a significant main effect of group on MFC GABA+ concentrations (F2,28 = 4.26, p = .02). Planned comparisons revealed elevated MFC GABA+ in the FEP group in comparison with the HC and APS groups (t27= 2.90, p = .01). There was a trend towards elevated MFC GABA+ in the APS group when compared with healthy controls (t27= 2.00, p = .0548). We did not find significant effects of group on rostral MFC Glx (F2,31 = 0.26, p = .77), or occipital GABA+ (F2,31 = 0.48, p = .62) or Glx (F2,32 = 1.03, p = .23) when PSI9 was included in the model.
3.4 Relationships between GABA+ and anti-psychotic medication dose
In the FEP group, correlations between rostral MFC GABA+ and CPZeq were r = .53, p = .11. Correlations between midline occipital GABA+ and CPZeq were r = .07, p = .81. Since only one participant with APS was receiving anti-psychotic medication, we did not calculate correlations between GABA+ and CPZeq in this group.
4. DISCUSSION
Using MRS, we investigated rostral MFC and midline occipital cortex GABA+ and Glx concentrations in participants with FEP and APS. With negative affect (NA), measured by PSI9, as a covariate, we found significantly elevated rostral MFC GABA+ in patients with FEP as well as a trend towards elevated rostral MFC GABA+ in individuals meeting APS criteria. This elevation of GABA +concentrations is consistent with prior investigations of rostral GABA in acute schizophrenia, in unmedicated individuals (T. Chen et al., 2017; de la Fuente-Sandoval et al., 2018; Kegeles et al., 2012), although some groups have also reported decreased GABA concentrations in the region of this rostral voxel (Chiu et al., 2018; Wang et al., 2016).
Prior divergent MRS findings and the relatively small number of studies leads us to believe that there are multiple factors influencing GABA levels in psychosis. Synthesizing the results of the current study with previously published reports from our group in which a benzodiazepine challenge was used to demonstrate an abnormal BOLD response in schizophrenia which correlated with levels of negative affect (Taylor et al., 2014; Tso et al., 2015), we suggest there is developing evidence to support at least two separate pathological processes involving rostral MFC GABA+ in psychosis.
4.1 Increased GABA+ concentrations in the rostral MFC in early psychosis
The first process is distinguished by elevated GABA+ concentrations in the rostral MFC, which occur early in the illness and may be driven down by treatment or other factors which has led to contradictory findings in the literature. This process may reflect a primary pathophysiological process or a compensatory response – for example, to environmental stressors. A chronic unpredictable stress paradigm in rodents has been shown to increase ACC GABA levels, as measured by MRS (Perrine et al., 2014). Alternatively, increased GABA may be a secondary response to glutamatergic dysfunction. For example, ketamine infusion, which blocks excitatory NMDA receptors, has been shown to increase medial prefrontal cortex GABA levels in humans (Rodriguez et al., 2015). We examined Glx concentrations in our sample and found no evidence for abnormal concentrations in psychosis. Glx, however, is the combined signal of glutamate and its precursor glutamine, and a balanced cycling between these two amino acids is essential for normal brain functioning. Investigating Glx concentrations may miss subtle shifts in the balances between glutamate (Glu) and glutamine (Gln), and therefore fail to detect glutamatergic dysfunction. A recent meta-analysis of glutamatergic metabolites does provide some evidence for elevated Glx and Glu in treatment resistant schizophrenia (Nakahara et al., 2021), lending support to the hypothesis of glutamatergic dysfunction in psychosis.
When comparing the APS group with controls, we found a trend suggesting increased rostral MFC GABA+ in APS, controlling for levels of negative affect. While this effect did not reach significance and should be interpreted with caution, our APS sample was small, and this analysis was undoubtedly underpowered. This trend is in line with previous evidence indicating elevated rostral MFC GABA+ in individuals at an ultra-high risk of developing psychosis (De La Fuente-Sandoval et al., 2015), and indicates that processes resulting in abnormal GABA concentrations in the rostral MFC may be present prior to the conversion to full psychosis.
4.2 Increased Negative Affect is associated with reduced GABA+ concentrations
As we note, multiple factors may drive GABA+ levels down and obscure group differences, and we propose there is evidence for a second process in which increased NA is associated with reduced rostral MFC GABA+, as well as the altered BOLD response observed in benzodiazepine challenge studies published by our group (Taylor et al., 2014; Tso et al., 2015). In the current study, we observed significantly higher levels of NA in our FEP group when compared with controls, and once we controlled for levels of NA in our statistical analysis, we found significant group differences in rostral MFC GABA+ concentrations. These differences were at a trend level without the inclusion of levels of NA in our model, indicating that NA may have been driving GABA+ levels down and obscuring group differences. The existence of two processes affecting GABA+ levels could explain the otherwise confusing finding from de la Fuente-Sandoval and colleagues, who found an inverse correlation between rostral MFC GABA levels and PANSS general psychopathology scores – which also capture NA - in first episode psychosis patients, who otherwise had higher GABA levels (de la Fuente-Sandoval et al., 2018).
Evidence from benzodiazepine challenge studies has linked greater levels of NA with altered GABAergic function in the dorsomedial prefrontal cortex while viewing salient stimuli (Taylor et al., 2014) and in the visual cortex during a face processing task (Tso et al., 2015). These findings indicate that this process may be more widespread throughout the cortex rather than confined to frontal regions, although we found no evidence for a relationship with NA in our midline occipital MRS voxel. In addition, Taylor and colleagues (2014) found that negative affect correlated positively with the altered BOLD signal in both patients and controls, suggesting that this process may be present on a continuum from psychosis to disease.
Reviewing the literature provides further support for a link between reduced GABA and NA. Considerable evidence has supported the involvement of reduced GABAergic activity in patients with major depressive disorder (MDD), who have an abundance of NA (Duman et al., 2019). GABA MRS studies have shown reduced GABA in MDD patients, including the MFC where we placed our voxels (Gabbey et al., 2013; Kantrowitz et al., 2021) .The demonstration that neuropeptides which potentiate GABAergic transmission have rapid antidepressant effects provides additional evidence of a role for GABA in NA (Fee et al., 2017) in experimental paradigms, inducing NA with a fear-of shock challenge acutely decreased rostral MFC GABA in healthy controls (Hasler et al., 2010). GAD1 (the gene coding for the GAD67 protein) knock-down mouse models have demonstrated impaired fear extinction (Brown et al., 2015) and emotion regulation but preserved working memory function (Kolata et al., 2018). Although as the major inhibitory neurotransmitter, GABA activity is tightly coupled in an excitatory-inhibitory balance which is difficult to easily characterize in pathological states, multiple lines of evidence show signs of reduced GABA activity, linked to various NA states.
4.3 Limitations and Future Directions
A key limitation of this study was its modest sample size, which means the results we describe should be interpreted with the appropriate caution, and replication of these effects is vital. As we note above, the small sample of participants with APS perhaps meant that our exploration of GABA+ in these individuals was underpowered. In addition, the majority of our FEP group and one member of our APS group were receiving anti-psychotic medication, which meant our sample was not free of medication confounds. Since our sample size was relatively small, we lacked the statistical power to probe the influence of medication on GABA+ concentrations while investigating linkages between GABA+ and negative affect. Prior research indicates that the relationship between GABA+ concentrations and anti-psychotic medications is perhaps complex, with evidence showing a negative correlation between GABA+ levels in the ACC and anti-psychotic dose (Tayoshi et al., 2010) and reductions in mPFC GABA levels following a 4 week course of atypical anti-psychotics (de la Fuente-Sandoval et al., 2018), contrasted with findings that a 6 month treatment with atypical anti-psychotics did not affect GABAergic transmission in the frontal lobe, left basal ganglia or parietal-occipital lobe (Goto et al., 2010). Further investigations of linkages between GABA concentrations and negative affect would benefit from a larger sample size which allows the impact of medication to be carefully explored. Additionally, future work should extend our MRS findings, and the findings from our group’s pharmacological probe of the impact of lorazepam on the BOLD response in the MFC, to investigate if reduced GABA levels in the rostral MFC are associated with an altered BOLD response in the same individuals.
4.4 Conclusions
In conclusion, we found evidence for significantly elevated rostral MFC GABA+ concentrations in psychosis, when controlling for levels of negative affect. When considered alongside previously published demonstrations of altered BOLD signal in response to a benzodiazepine challenge, as well as linkages between GABA and negative affective states, we suggest the presence of two pathological processes involving GABA. The first is distinct from negative affect and is characterized by an increase of GABA in the rostral MFC which occurs early in the illness. The second process is reflected by reduced GABA which is associated with increased negative affect and an altered BOLD response to benzodiazepine. While more data is required to untangle these processes, and clarify GABAergic dysfunction in the psychosis spectrum, NA is a strong, transdiagnostic predictor of functional outcome, and these findings establish potential therapeutic leverage points linking GABAergic treatments, with NA as the clinical outcome.
Data Availability
All data produced in the present study are available upon reasonable request to the authors