ABSTRACT
Following the first episode of psychosis, some patients develop poor social and occupational outcomes, while others display a pattern of preserved functioning. Several lines of evidence from preclinical, genetic and biochemical studies suggest a role for high oxidative stress in poor functional outcomes. The measurement of intracortical glutathione (GSH) using magnetic resonance spectroscopy (MRS) provides an opportunity to investigate the relationship between central antioxidant tone and functional outcomes at the time of first episode psychosis (FEP). A body of epidemiological studies indicates better functional outcomes in patients at early stages of schizophrenia compared to patients at a chronic, established phase of illness. We scanned 57 patients with FEP and 30 matched healthy controls and estimated GSH resonance using 7-Tesla MRS. We minimised the confounding effects of illness chronicity, long-term treatment exposure and metabolic complications by recruiting patients with <2 weeks of lifetime antipsychotic exposure on average and followed up this cohort for the next 1 year to determine functional outcomes. Patients with FEP who achieved employment/education or training status (EET) in the first year, had higher GSH at the baseline than healthy controls. Social and occupational functioning assessment scale (SOFAS) scores were also significantly higher in patients with higher GSH levels at the outset, after adjusting for various confounds including baseline SOFAS. Patients who were not in employment, education or training (NEET) did not differ from healthy subjects in their GSH levels. Our observations support a key role for the central antioxidant tone in the functional outcomes of early psychosis.
INTRODUCTION
For patients with schizophrenia, the probability of functional recovery is highest at the early stages of the illness, around the time of the first psychotic episode1–3; when a chronic pattern of the illness gets established, the chances of functional recovery greatly diminish, with only a small sub-group (∼13%) recovering at this stage 4,5. Currently, we do not know what mechanistic processes underlie these diminishing returns in recovery rates over time. Several clinical characteristics (e.g., the presence of negative, disorganised symptoms, cognitive deficits6) have been observed in association with poor functional outcomes; in particular, the degree of functioning at first presentation (baseline or premorbid) explains a significant amount of variance in long term functional outcomes7–9. Despite their explanatory power, these clinical associations do not offer an actionable mechanistic marker that can be harnessed for therapeutic purposes. There is an urgent need to understand the neurobiological factors that contribute to differences in functional outcomes in early stages of illness.
One promising approach to study variable outcomes in psychosis is quantifying the relative burden of oxidative stress experienced by patients during the first psychotic episode10. Fournier and colleagues utilised a data-driven stratification procedure on a cohort of patients with early psychosis and identified a subgroup with superior functional outcomes; this subgroup was characterised by lower levels of oxidative stress markers (especially glutathione peroxidase) in the blood11. Lower baseline blood levels of glutathione (GSH) predict later cognitive deficits12 as well as brain volume loss in early psychosis13. While peripheral markers of oxidative stress correlate with concentration of the primary intracortical antioxidant glutathione14, a direct link between central glutathione levels and functional outcome in first episode psychosis (FEP) is yet to be demonstrated. Wood and colleagues15 reported a 22% increase in medial temporal GSH levels in first episode psychosis; in a sub-sample from this study, treatment related increase in GSH was associated with a gain in global functioning scores16. In a small group of individuals with various mental health difficulties indicating a high-risk state for psychosis, we recently demonstrated higher GSH levels measured using magnetic resonance spectroscopy (MRS) from the dorsal anterior cingulate cortex (ACC) in those with better social and occupational functioning17. In the current study, we investigate if ACC GSH levels measured at the onset of psychosis, before establishing regular antipsychotic treatment, predict later functional outcomes in the first year of early intervention.
Early functional outcome status is a well-established indicator of long-term course of schizophrenia18,19. An exciting translational possibility of linking GSH levels with functional outcomes in psychosis, is the availability of targeted treatments that can improve intracortical GSH in patients. Several clinical trials have reported on the safety and efficacy of the glutathione precursor N-acetylcysteine (NAC) in patients with psychosis20. These trials (6 RCTs)21 indicate that NAC produces a modest, but significant improvement in cognitive deficits and negative symptoms (critical determinants of poor functional outcomes), when used as an adjunct to antipsychotics. Thus, in patients with psychosis, a deficit in intracortical GSH is likely to be a potentially modifiable pathway of poor outcomes.
METHODS
Participants
The sample for the present analysis consisted of 72 new referrals to the PEPP (Prevention and Early Intervention for Psychosis Program) at London Health Sciences Center. After exclusions were made due to missing/poor quality scan data our final sample consisted of 57 patients (48 males/9 females) (Table 1). All participants provided written, informed consent prior to participation as per approval provided by the Western University Health Sciences Research Ethics Board, London, Ontario. Inclusion criteria for study participation were as follows: individuals experiencing first episode psychosis, with not more than 14 days of cumulative lifetime antipsychotic exposure, no major head injuries (leading to a significant period of unconsciousness or seizures), or known neurological disorders, and no concurrent substance use disorder.
The mean lifetime total defined daily dose days (DDD × days on medication) for antipsychotic use was 1.80 days with 27 patients (47.4%) being completely antipsychotic naive at the time of scanning. Of those who had started antipsychotic treatment, (N=30; 52.6%), the median total defined daily dose days was 2.81 days (range of 0.4–14 DDD days). Patient consensus diagnosis was established using the best estimate procedure described by Leckman et al.22 following 6 months of treatment. Diagnoses were made based on the Structured Clinical Interview for DSM-5.
Healthy control subjects (n=30) were recruited through posters and word of mouth advertising. Healthy control subjects had no personal history of mental illness, no current use of medications, and no first-degree relatives with a history of psychotic disorders. Healthy controls were group matched to the FEP cohort for age and parental socio-economic status (the National Statistics Socioeconomic Classification: five-class version23). Similar to their FEP counterparts, those with a history of substance use disorders in the past 12 months, significant head injury or neurological disorders were excluded.
Clinical Measures
A clinical examination was conducted during the baseline assessment (the same day imaging was performed) by a research psychiatrist for patients or a trained rater for healthy controls. The clinical battery was used to assess patient symptom severity and to ensure that control subjects were free from current Axis I disorders and history of either psychotic illness or neurologic disorder. To assess symptoms of psychosis the Positive and Negative Syndrome Scale – 8 Item (PANSS-8) was used24. The PANSS-8 is an abbreviated version of the 30 Item PANSS clinical assessment of symptomology in schizophrenia and psychosis with acceptable internal consistency and highly correlated with the full PANSS24. Items are scored on a 1 (absent) to 7 (extreme) Likert type scale, assessing both positive (P1-delusions, P2 – Conceptual disorganization, & P3-hallucinations), and negative domains (N1-Blunted or flat affect, N4-passive social withdrawal, and N6-impoverishment of speech).
Additionally, the Social and Occupational Functioning Assessment Scale (SOFAS) was administered at baseline and follow up. The SOFAS is a single item measure of functioning scored between 1 (persistent inability to maintain minimum functioning without external support) and 100 (superior functioning in a wide range of activities). The SOFAS ratings of social and occupational functioning were made independent of symptom severity. In our study SOFAS was taken with consideration for current functioning (rather than highest level of functioning over the past year).
We obtained baseline SOFAS on the same day MRI data were acquired. The patients enrolled in the PEPP clinic were followed up over the next 12 months, and we ascertained their NEET status and follow-up SOFAS scores between 6 to 12 months. The functional assessment was based on multiple sources of information: patient interviews, information from the psychiatrist providing the clinical care, PEPP case managers and, where required, information from family members documented in clinical charts. Due to the need for multiple information sources, not all patients were assessed at the same time point after their illness onset, but the vocational (NEET) status of the cohort between the window of 6 to 12 months was captured, in addition to baseline and follow-up assessment of cross-sectional functioning using SOFAS. Patients were classified as NEET (vocationally inactive) if they were unemployed and not in any form of schooling/education for more than half of the time since the onset of treatment for psychosis. Individuals who were engaged in work or school for more than half of the duration of treatment were classified as EET (vocationally active). This definition considers a longer time frame (up to 6 months) than the 1-week period used by the Organisation for Economic Co-Operation and Development (OECD25), in line with its use in early intervention services for psychosis26,27. A consensus was reached within the research study team when discrepancies noted in reported functioning between patient’s accounts and those of care providers.
MRS assessment
A complete description of the MRS protocol used in the present study has been described in previous work by Dempster and colleagues28. Metabolite concentrations (glutamate and GSH) were estimated using single-voxel 1H-MRS data acquired with a Siemens MAGNETOM 7.0T head-only MRI scanner (Siemens, Erlangen, Germany) using an 8-channel transmit/32-channel receive head coil at the Center for Functional and Metabolic Mapping of Western University in London, Ontario. A 2.0 × 2.0 × 2.0 cm (8 cm3) voxel was placed on the bilateral dorsal ACC using a two-dimensional anatomical images acquired in the sagittal direction (37 slices, TR = 8000 ms, TE = 70 ms, flip-angle (α) = 120°, thickness = 3.5 mm, field of view = 240 × 191 mm). The posterior end of the voxel was set to coincide with the precentral gyrus and the caudal face of the voxel coincided with the most caudal location not part of the corpus callosum. The angulation of the voxel was determined to be tangential to the corpus callosum (Figure 1).
Statistical analyses
We used IBM SPSS Statistics version 24 for all analyses. First, the primary hypothesis of a relationship between NEET status in the first year and glutathione was tested using a one-way ANOVA, with healthy controls, NEET and EET patients as 3 groups of interest. Second, within the patient group, a bivariate correlation between follow-up SOFAS scores and GSH measurement was conducted, with bootstrapping for generating p-values and confidence intervals. We also tested if GSH levels retain their ability to predict follow-up SOFAS, after adjusting for the variance explained by baseline functioning (SOFAS at the time of MRS scanning), age and gender as covariates. Finally, we excluded patients without a clear diagnosis of schizophrenia by 6-12 months assessment and tested the relationship between GSH and SOFAS and NEET.
Results
When comparing healthy controls to FEPs, no statistically significant differences existed for age, parental socioeconomic status, or gender, although male participants were overrepresented in both groups (Table 1). As expected, there were several group-level differences in clinical scores between patient groups and healthy controls (Table 1). No significant differences in baseline demographic or clinical differences were identified between EET and NEET patients, with the exception of baseline SOFAS scores, which were higher among EET patients (t= 2.23, p= 0.031). Including subjects without useful MRS data did not affect this demographic profile (Supplementary Table 1).
Group differences in GSH levels
Glutathione levels were significantly different among the 3 groups (FEP-NEET, FEP-EET and HC) at the p<0.05 level [F(2,84)= 4.55, p= 0.01]. Post hoc comparisons using the Sidak test indicated that the mean GSH levels for FEP-EET subjects was significantly higher than healthy controls [[Mean(SD) of GSH: FEP-EET= 1.76(0.37), HC= 1.49(0.32); p=0.01]. However, FEP-NEET subjects [M(SD) =1.60(0.34)] did not significantly differ from the healthy controls (p= 0.25) or FEP-EET (p= 0.57). Of note, patients as a whole had higher GSH levels than healthy controls (FEP= 1.68 (0.36), HC= 1.49(0.32); t(85)= 2.46, p= 0.016). These results are shown in Figure 2.
Within the FEP group, GSH measure at baseline was positively correlated with follow-up SOFAS scores r(44) = 0.41 [bootstrap 95% CI: 0.2-0.61], p=0.006. GSH levels, in combination with age, gender and baseline SOFAS scores, predicted follow-up SOFAS (R2= 0.29, F(4,36)= 3.76, p= 0.013). GSH levels (β= 0.36, p= 0.016) showed the strongest association out of all independent predictors for follow-up SOFAS, after adjusting for the variance explained by baseline functioning, age and gender (β= 0.17-0.21, all p> 0.15). The correlation between adjusted GSH levels and follow-up SOFAS is displayed, stratified by NEET status, in Figure 3.
Analysis restricted to first episode schizophrenia only
Finally, we repeated the above analysis in a subset of patients with a consensus diagnosis of schizophrenia/schizoaffective disorder by 6-12 months, after excluding 5 subjects with major depressive disorder, schizophreniform or bipolar disorder. Glutathione levels continued to differ significantly among the 3 groups (FES-NEET, FES-EET and HC) at the p<0.05 level [F(2,79)= 4.16, p= 0.02]. Post hoc comparisons using the Sidak test indicated that the mean GSH levels for FES-EET subjects was significantly higher than healthy controls [[Mean(SD) of GSH in FES-EET = 1.75(0.36), HC= 1.49(0.32); p= 0.016]. However, FES-NEET subjects [M(SD) = 1.59(0.35)] did not significantly differ from the healthy controls (p= 0.26) or FES-EET (p= 0.67).
Within the FES group, GSH measure at baseline was positively correlated with follow-up SOFAS scores r(39)= 0.38 [bootstrap 95% CI: 0.15-0.61], p= 0.017. GSH levels, in combination with age, gender and baseline SOFAS scores, continued to predict follow-up SOFAS in the FES group at a trend level, (R2= 0.24, F(4,32)= 2.47, p= 0.06). GSH level (β= 0.32, p= 0.05) continued to be the most prominent independent predictor of the follow-up SOFAS, after adjusting for the variance explained by baseline functioning, age and gender (β= 0.10-0.23, all p> 0.16).
Prognostic relevance based on binarized GSH levels
To provide clinically relevant information to a patient with FEP whose GSH values have been quantified as described in this study, we performed a median split analysis of the patient group based on ACC GSH values of the entire sample (median= 1.586mM). We then compared the low-GSH and high-GSH FEP groups on the clinically relevant variables of NEET status and SOFAS at follow-up. The proportion of patients with FEP who were vocationally active (i.e., EET) significantly differed based on their baseline GSH status, with high-GSH FEP reporting 68% (23 of 34) in EET, while low-GSH FEP reporting 35% (8 of 23) in EET (Fisher’s exact p= 0.018). This translated to an odds ratio of 3.92 (95%CI: 1.34-11.2), and a relative risk of 1.95 (95%CI: 1.13-3.7) in favour of being in employment, education or training if an individual with FEP belonged to the high, instead of low-GSH group at the outset. The mean follow-up SOFAS scores for high-GSH FEP was also significantly higher than low-GSH FEP [Mean(SD) of SOFAS: high-GSH= 62.7(12.9), low-GSH= 53.5 (13.4); t (42)= 2.24; p= 0.03].
DISCUSSION
Using 7T MRS to measure GSH concentration from ACC in untreated FEP subjects followed up over 1 year period, we report 2 major findings: (1) GSH levels are higher in patients with FEP when compared to healthy controls; this difference was driven by those who later displayed superior vocational status and socio-occupational functioning and (2) the baseline GSH levels predict later socio-occupational functioning over and above what can be predicted from baseline functioning, indicating a specific mechanistic role for the central antioxidant status in the outcome trajectories following first psychotic episode. Patients with GSH levels that were lower than the median values in the sample were at approximately two times higher risk of being vocationally inactive, with approximately 10 points lower scores in SOFAS.
Our observed relationship between lower GSH levels and lower SOFAS scores and NEET status is in line with several prior reports. In an overlapping sample, we previously reported a predictive relationship between low GSH and delayed response to antipsychotics28. Lack of early response is a critical indicator of long-term poor outcomes in schizophrenia29–31. In line with studies linking lower GSH to higher residual symptom burden32, negative symptoms33 and cognitive deficits34 in schizophrenia, our observation highlights a prominent ‘pathoplastic’ role for antioxidant status in shaping the outcomes of this illness. In fact, a recent 7T-MRS study of the ACC reported that patients who fail first-line antipsychotic treatments are more likely to have an intracortical GSH-deficit35. In our study, the presence of higher ACC GSH in FEP may indicate the overall superior treatment responsiveness in this cohort, given their untreated status. Such higher GSH levels have been previously reported in other brain regions in FEP15,16. The small to medium effect size difference also underscores the considerable heterogeneity in intracortical GSH levels in psychosis, in line with several 3T MRS studies that did not report a notable group difference in the ACC36 as well as other regions37,38. This heterogeneity may relate to individual differences (i.e., subgroups of patients with high or low GSH levels) or stage-specific differences (acute excess, followed by later deficit). While several exogenous factors may also affect antioxidant levels, there is no evidence for a major role for these confounders in explaining aberrations in glutathione pathway in psychosis39.
An important strength of our study is its longitudinal nature; with the temporal information we can establish that higher GSH levels at the outset predict and thus influence superior socio-occupational functioning seen over the next 1 year among patients. But a major limitation is the availability of a single time point of GSH levels. To establish a more conclusive causal inference, we need longitudinal follow-up studies that capture multiple time points from untreated early stages of psychosis to a stable phase when functional outcome trajectories become more established. Long term follow-up of experimental studies with patients in early stages of psychosis that receive adjunctive NAC will also be illuminating in this regard. We also lacked peripheral and genotyping measures of antioxidant capacity in this cohort; while these measures are more accessible than 7T-MRS, they do not consistently reflect intracortical GSH40–42.
We conclude that in patients with first episode psychosis in whom intracortical GSH levels are higher during the acute phase of psychosis, functional outcomes are superior to those with lower levels, over the next 1 year. Thus, ACC GSH measure may be a useful indicator of resilience to oxidative stress and functional recovery after first episode of psychosis. Taken together with prior MRS studies from established cases of schizophrenia indicating a profile of treatment failures and residual symptoms in patients with intracortical GSH-deficit, our observation makes a compelling case to investigate the role of pre-emptive antioxidant interventions in early stages of psychosis.
Data Availability
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.
Institutional Review Board Statement
This study was conducted according to the guidelines of the declaration of Helsinki and approved by the Research Ethics Board of the University of Western Ontario (Project #108268; October 19, 2020).
Funding
This study was funded by CIHR Foundation Grant (375104/2017) to LP; AMOSO Opportunities fund to LP; partial salary support of PJ by NSERC Discovery Grant (RGPIN-2016-05055) to JT, Schulich School of Medicine and Dentistry Dean’s Scholarship to PJ; Parkwood Institute Studentship and the Jonathan and Joshua Memorial Scholarship to MM. Data acquisition was supported by the Canada First Excellence Research Fund to BrainSCAN, Western University (Imaging Core); Innovation fund for Academic Medical Organization of Southwest Ontario; Bucke Family Fund, The Chrysalis Foundation and The Arcangelo Rea Family Foundation (London, Ontario). LP acknowledges salary support from the Tanna Schulich Chair of Neuroscience and Mental Health.
Author Contributions
MM recruited patients, collected, and compiled demographic, functional outcome and treatment data and assisted in acquiring the scans. SF recruited patients, compiled demographic, functional outcome and treatment data and assisted in acquiring the scans. PJ acquired the MRS scans, performed the MRS quantification and analysis, supervised by JT. LP obtained the research funds, designed the project with JT and supervised MM, SF and PJ. LP conceived the analysis and wrote the first draft with MM and SF. All authors contributed to the writing and have approved the final version of the manuscript.
Conflict of Interest
LP reports personal fees from Otsuka Canada, SPMM Course Limited, UK, Canadian Psychiatric Association; book royalties from Oxford University Press; investigator-initiated educational grants from Janssen Canada, Sunovion and Otsuka Canada outside the submitted work. All other authors report no relevant conflicts.
Acknowledgements
We thank Mr. Trevor Szekeres, Mr. Scott Charlton, Mr. Joseph Gati for their assistance in data acquisition and archiving. We thank all research team members of the NIMI lab and all the staff members of the PEPP London team, particularly Drs. Kara Dempster, Julie Richard, Priya Subramanian and Hooman Ganjavi for their assistance in patient recruitment and supporting clinical care. We gratefully acknowledge the participants and their family members for their contributions.
Footnotes
↵* Shared first authors