Next Article in Journal
Spatiotemporal Dynamics of Covert Versus Overt Processing of Happy, Fearful and Sad Facial Expressions
Next Article in Special Issue
Role of Vitamin E and the Orexin System in Neuroprotection
Previous Article in Journal
Disrupted Functional Rich-Club Organization of the Brain Networks in Children with Attention-Deficit/Hyperactivity Disorder, a Resting-State EEG Study
Previous Article in Special Issue
The Current Role of Dexmedetomidine as Neuroprotective Agent: An Updated Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Glutathione as a Molecular Marker of Functional Impairment in Patients with At-Risk Mental State: 7-Tesla 1H-MRS Study

1
Department of Medical Biophysics, Western University, London, ON N6A 3K7, Canada
2
Lawson Health Research Institute, Imaging Division, London, ON N6A 4V2, Canada
3
Robarts Research Institute, Western University, London, ON N6A 3K7, Canada
4
Department of Psychiatry, Western University, London, ON N6A 3K7, Canada
5
Department of Neuroscience, Western University, London, ON N6A 3K7, Canada
6
Department of Psychology, Western University, London, ON N6A 3K7, Canada
7
St. Joseph’s Health Care, Diagnostic Imaging, London, ON N6A 4V2, Canada
8
Department of Medical Imaging, Western University, London, ON N6A 3K7, Canada
*
Author to whom correspondence should be addressed.
Submission received: 6 April 2021 / Revised: 2 July 2021 / Accepted: 13 July 2021 / Published: 17 July 2021
(This article belongs to the Special Issue Neuroprotection and Oxidative Distress)

Abstract

:
A substantial number of individuals with clinical high-risk (CHR) mental state do not transition to psychosis. However, regardless of future diagnostic trajectories, many of these individuals develop poor social and occupational functional outcomes. The levels of glutathione, a crucial cortical antioxidant, may track variations in functional outcomes in early psychosis and prodromal states. Thirteen clinical high-risk and 30 healthy control volunteers were recruited for a 7-Tesla magnetic resonance spectroscopy scan with a voxel positioned within the dorsal anterior cingulate cortex (ACC). Clinical assessment scores were collected to determine if any association was observable with glutathione levels. The Bayesian Spearman’s test revealed a positive association between the Social and Occupational Functioning Assessment Scale (SOFAS) and the glutathione concentration in the clinical high-risk group but not in the healthy control group. After accounting for variations in the SOFAS scores, the CHR group had higher GSH levels than the healthy subjects. This study is the first to use 7-Tesla magnetic resonance spectroscopy to test whether ACC glutathione levels relate to social and occupational functioning in a clinically high-risk group and offers preliminary support for glutathione levels as a clinically actionable marker of prognosis in emerging adults presenting with risk features for various severe mental illnesses.

1. Introduction

Emerging adults with attenuated or brief and limited psychotic symptoms are said to be in a clinical high-risk (CHR) state (or “at-risk or ultrahigh-risk” mental state) that later develops into multiple diagnostic outcomes including schizophrenia, mood disorders, such as bipolar disorder, or a major depressive disorder [1,2]. A substantial number of individuals with CHR develop poor long-term functional (i.e., social and occupational) outcomes irrespective of diagnostic transitions. Longitudinal studies indicate that a large proportion of individuals with CHR do not transition to psychosis (65–89% not psychotic over 2–10 years [3,4,5]) but have poor social and occupational outcomes (48% functionally impaired at 3–10 years [3,6]). While functional outcomes improve over time in CHR patients who have good functioning at the baseline, persistent deficits are seen in those who start with lower levels of functioning [4]. In other words, lower levels of functioning in CHR state before the onset of diagnosable psychiatric disorders such as schizophrenia and bipolar disorders reliably predict the trajectory of continued poor functioning over a long time period. The molecular bases of such pervasive functional deficits continue to be unknown [7], proving to be a major hurdle in developing meaningful treatments aimed at CHR state.
Oxidative stress has emerged as a key mechanism underlying the pathophysiology of many psychiatric disorders including psychosis [8]. Destructive free radicals that damage brain tissue are by-products of oxidative metabolism but are effectively scavenged by antioxidants. Glutathione (GSH), the cardinal antioxidant in brain cells, shows 27–52% reduction [9,10,11] in established schizophrenia. Genetic [12,13] and cell biology studies [14,15,16] indicate that in a subset of patients, GSH production on demand is likely to be reduced [17]. We have recently demonstrated the prognostic importance of low GSH in predicting early clinical response to antipsychotics in first-episode schizophrenia. In this study, we observed that for every 10% baseline difference in anterior cingulate cortex (ACC) GSH among the patients, seven additional days of delay in response occurred after treatment initiation [18]. Lack of early response is a critical indicator of long-term poor outcomes in schizophrenia [19,20,21]. We and others have also related lower GSH to various determinants of functional outcomes including residual symptom burden [22], negative symptoms [23] and cognitive deficits [24], supporting the notion that the “hub of oxidative stress” indexed by GSH [25,26] is likely a critical determinant of functioning.
To date, various imaging modalities have been employed to predict both conversion to psychosis and functional outcomes in youth at CHR [27]. Structural MRI studies report pronounced cortical thinning [28,29] and longitudinal reduction in gray matter volume [30] in the ACC of CHR patients. These ACC abnormalities precede the psychosis onset [28,31]. Furthermore, aberrant activation of ACC in fMRI studies using tasks that demand executive controls has also been reported in CHR subjects who develop poor symptomatic and functional outcomes [32,33,34]. While several other distributed structural and functional changes precede and predict psychosis (see the review by Andreou and Borgwardt [35]), a recent meta-analysis has established that dorsal ACC (dACC) is the only brain region whose structural gray matter changes precede psychosis (i.e., in CHR state) and persist even after long periods of treatment (i.e., in both acute and chronic stages) [36], highlighting the crucial role of ACC for long-term outcomes.
We recently synthesized in vivo MRS studies on ACC and demonstrated a significant GSH reduction in the established cases of schizophrenia but an elevation in bipolar disorders [37], indicating that GSH levels may track the variations in functional outcomes that typify the prognostic course of psychiatric disorders. Such divergence between disorders may mean that in the “pluripotent” CHR state that includes patients with varying levels of functioning as a single group, GSH levels may not differ from healthy controls but relate to variations in levels of functioning. In fact, in the only previous MRS study of cortical glutathione in clinical high-risk state [38], da Silva et al. reported no difference between healthy controls and CHR subjects in anterior cingulate glutathione [39]. Functional outcomes were not evaluated in this study; thus, the role of GSH as a transdiagnostic prognostic marker in CHR state remains unknown.
In the current study, we use ultrahigh field 7T MRS for the first time to test if ACC GSH levels relate to social and occupational functioning in the CHR group. We expected GSH levels to be reduced among patients with poor functioning. Furthermore, we aimed at establishing the difference in GSH levels between the CHR and healthy control subjects after accounting for variations in functioning. We evaluated these hypotheses using Bayesian and frequentist statistical approaches.

2. Materials and Methods

2.1. Participants

We recruited 13 clinical high-risk (CHR) volunteers along with 30 healthy control (HC) volunteers group-matched for age, gender and parental socio-economic status. Patient volunteers were recruited from the referrals received by the PROSPECT (Prodromal Symptoms of Psychosis—Early Clinical Identification and Treatment) program at London Health Sciences Center, London, Ontario. The patients were help-seeking individuals referred to the clinic by community physicians, healthcare workers or friends/family. All the referrals were reviewed by an intake coordinator via telephone using a validated instrument (PRIME Screen—Revised). If found eligible for further assessment, the patients were evaluated within 2 weeks of referral using the Structured Interview for Psychosis-risk Syndromes (SIPS) [40]. Patients with medical conditions, pervasive developmental disorders or intellectual disability underlying the reported symptoms, those who received treatment with antipsychotic medications to treat presenting symptoms (minimal effective dose for a period of at least 2 weeks) and those with psychotic symptoms secondary to active substance use (intoxication effects) were excluded. Based on the SIPS, the patients satisfying attenuated psychotic syndrome (APS) or brief and limited intermittent psychosis (BLIPS) were both included in the CHR group. The healthy volunteers had no personal history of mental illness with no family history of psychotic disorder. All the participants were screened to exclude significant head injury, major medical illness or MRI contraindications and provided written informed consent according to the guidelines of the Human Research Ethics Board for Health Sciences at Western University, London, Ontario.

2.2. MRS Acquisition and Analysis

A Siemens MAGNETOM 7T head-only MRI scanner (Siemens, Erlangen, Germany) was used for all MRS acquisition along with a site-built head coil (8-channel transmit, 32-channel receive coil array) at the Centre for Functional and Metabolic Mapping of Western University (London, ON, Canada). A two-dimensional sagittal anatomical image (37 slices, TR = 8000 ms, TE = 70 ms, flip-angle (α) = 120°, thickness = 3.5 mm, field of view = 240 × 191 mm2) was used as reference to prescribe a 2.0 × 2.0 × 2.0 cm3 (8 cm3) 1H-MRS voxel on the bilateral dorsal ACC (Figure 1). The voxel position was prescribed by setting the posterior face of the voxel to coincide with the precentral gyrus and setting the position of the inferior face of the voxel to the most caudal point not part of the corpus callosum. The voxel angle was set to be tangential to the corpus callosum. A semi-LASER 1H-MRS sequence (TR = 7500 ms, TE = 100 ms, bandwidth = 6000 Hz, N = 2048) was used to acquire 32 channel-combined, VAPOR [41] water-suppressed spectra as well as a water-unsuppressed spectrum to be used for spectral editing and quantification. All the participants were asked to fix their gaze on a white cross (50% gray background) during MRS acquisition.
Using the techniques outlined by Near et al. [42], the 32 spectra were phase- and frequency-corrected before being averaged into a single spectrum to be used for all subsequent analyses. QUECC [43] and HSVD [44] were applied to the spectrum for lineshape deconvolution and removal of the residual water signal, respectively. Spectral fitting was done using fitMAN [45], a time-domain fitting algorithm that uses a nonlinear iterative Levenberg–Marquardt minimization algorithm to echo time-specific prior knowledge templates. The metabolite fitting template included 17 brain metabolites: alanine, aspartate, choline, creatine, γ-aminobutyric acid (GABA), glucose, glutamate, glutamine, glutathione, glycine, lactate, myo-inositol, N-acetylaspartate, N-acetylaspartylglutamate, phosphorylethanolamine, scyllo-inositol and taurine. Due to the long echo time used, no significant macromolecular contribution was expected. Metabolite quantification was then performed using Barstool [46] with corrections made for tissue-specific (gray matter, white matter, CSF) T1 and T2 relaxation through partial volume segmentation calculations of voxels mapped onto T1-weighted images acquired using a 0.75-mm isotropic MP2RAGE sequence (TR = 6000 ms, TI1 = 800 ms, TI2 = 2700 ms, flip-angle 1 (α1) = 4°, flip-angle 2 (α2) = 5°, FOV = 350 mm × 263 mm × 350 mm, Tacq = 9 min 38 s, iPATPE = 3 and 6/8 partial k-space). All spectral fit underwent visual quality inspection as well as the Cramer–Rao lower bounds (CRLB) assessment for each metabolite.
The quality of metabolite quantification was measured using CRLB percentages for both groups using a CRLB threshold < 30% for glutathione to determine inclusion toward further analyses, in line with our prior study [18]. Notably, the mean CRLB for these metabolites were over two times lower than the individual threshold percentages. There was no significant difference in CRLB between the clinical high-risk group and the healthy controls for both metabolites reported in this study. We present the concentration and CRLB of other metabolites in our fitting template, along with the two presently mentioned, in the Table S1. A sample of fitted spectrum for a single participant is presented in Figure 1.

2.3. Clinical Assessments

Symptom severity was measured using the scale of prodromal symptoms (SOPS) on the same day of the scan. We also quantified the overall social and occupational functioning at the time of first presentation using the SOFAS [47] administered on the same day of the scanning. To determine cannabis use in the previous six months, the Cannabis Abuse Screening Test (CAST) was used [48]. The CAST is a six-item Likert-scale self-report questionnaire which asks the participant about cannabis use and how it affects their daily activities and relationships. Scores range from 6 to 30, with higher scores indicating more cannabis use. To determine alcohol use in the previous 6 months, the Alcohol Use Disorders Identification Test—Concise (AUDIT—C) [49] was used. The AUDIT—C is a three-item Likert-scale self-report questionnaire which asks the participant about alcohol use frequency and quantity. Scores range from 0 to 12, with higher scores indicating more alcohol use. Alcohol users and nonusers were classified by AUDIT—C scores of four or more and less than four, respectively. Lastly, nicotine use in the previous six months was determined using the single-item Fagerström Test for Nicotine Dependence and the smoking index [50]. The Fagerström test indicates time to the first cigarette after waking, and the smoking index is calculated by multiplying the number of years regularly smoking by the number of cigarettes per day divided by 20 cigarettes per pack. A lower Fagerström test value indicates more nicotine dependence, and a higher smoking index indicates more nicotine use. The 10-item Drug Abuse Screening Test (DAST-10) [51] was also employed for substances other than cannabis, alcohol and nicotine, though our cohort did not endorse any such use.

2.4. Bayesian Analysis

We evaluated the association between GSH and the SOFAS, the SOPS and the CAST scores in the CHR group and the relationship between GSH and the SOFAS scores in the HC group by using a Bayesian Spearman’s test [52]. This approach relies on data augmentation via the Metropolis-within-Gibbs sampling algorithm. Briefly, we assumed the rank data as the reflection of a latent (truncated) normal distribution which allowed us to use a conventional likelihood function. That is to say, the latent continuous scores would manifest as “degraded” rank values. Following this assumption, the data augmentation algorithm would yield samples from a truncated posterior distribution. Here, we tested the null hypothesis that ρ = 0 versus the alternative hypothesis that ρ ~ Uniform [−1,1] (i.e., following a uniform prior distribution). We drew 11,000 samples using a Markov chain Monte Carlo (MCMC) method using the “spearmanCorrelation.R” function in R as specified by van Doorn et al. [52]. We reported the Bayes factor relative to the null model (BF10). BF10 > 1.0 suggests evidence in support of the (alternative) hypothesis and vice versa. We also reported the mode and the proportion of the posterior distribution (i.e., posterior proportion, PP) of the estimated ρ (rho) values differing from zero along with the 95% highest density interval of the most credible values (HDI).
We estimated the posterior distribution of the (estimated) between-groups differences in the CAST, the AUDIT—C and the SOFAS scores by means of a generalized linear model (GLM) within the context of hierarchical Bayesian parameter estimation as follows:
s c o r e s i = β 0 + g r o u p β g r o u p x g r o u p ( i )
where the data conformed to normal distribution around the predicted value ( s c o r e ) with (wide) data-scaled uniform prior distribution for the standard deviation (σi). The baseline parameter ( β 0 ) had a data-scaled normal prior distribution with the mean equal to the data mean and the (wide) standard deviation relative to the standard deviation (SDdata) of the data (1/(SDdata × 5)2). Group deflection parameters ( β g r o u p ) had normal prior distributions with mean zero and gamma prior distribution for the standard deviation σβ with data-scaled shape and rate parameters (SDdata/2 and 2 × SDdata, respectively). This means that σβ provided informed priors on each group (deflection) parameter. In other words, groups would act as priors between each other. In total, we estimated posterior distributions of five free parameters (σi, β 0 ,   β H C ,   β C H R and σβ ) . The posteriors were estimated in RJAGS using MCMC, drawing 11,000 samples (thinning = 10). We reported the PP of the between-groups difference in scores.
To evaluate the group effect after accounting for the effect of the SOFAS scores, we included those scores as a covariate in the GLM,
[ G S H ] i = β 0 + g r o u p β g r o u p x g r o u p ( i ) + β s o f a s x s o f a s ( i )
in which we added a normal prior distribution of the covariate parameter ( β s o f a s ) which had zero mean and the data-scaled standard deviation equal to 1/(2 × SDGSG_data/SDSOFAS_data)2. In total, we estimated posterior distributions of six free parameters (σi, β 0 ,   β H C ,   β C H R , σβ and β s o f a s ) .
The posteriors were estimated in RJAGS using MCMC, drawing 11,000 samples (thinning = 10). We reported the PP of the between-groups difference in GSH along with the 95% HDI. The posterior distribution of the effect size of this difference is also reported. There were two reasons for employing a Bayesian approach. First, we were interested in recovering the whole posterior distribution of possible values of the parameters given the data, constraining our discussion to the 95% highest density interval. Second, a Bayesian approach allowed us to identify the support not only for the alternative hypothesis, but also for the null hypothesis.

2.5. Frequentist Analysis

All frequentist statistical tests were computed using IBM SPSS Statistics version 26 [53]. Group demographic differences were calculated using t-tests and chi-squared tests for continuous and dichotomous variables, respectively. Hierarchical regression was used to assess the effect of the SOFAS scores and diagnosis (dummy coded CHR = 0, healthy controls = 1), with parameter estimates examined to test individual variable effects. Lastly, Spearman’s correlation was used to determine the correlation between metabolite levels and clinical scores (SOFAS, CAST, AUDIT—C, SOPS).

3. Results

3.1. Demographic Data

Subject demographic and clinical data are summarized in Table 1. A small number within the clinical high-risk group were being administered antidepressants (N = 3) or benzodiazepine (N = 2) at the time of the scan. Mean percent CRLB values for CHR and HC GSH were 10 ± 4% and 11 ± 4%, respectively. The CHR patients had substantially high levels of functional impairment. The groups differed in the CAST scores, which were higher in the CHR group than in the HC group (mode = 4.19, posterior proportion (PP) = 1.0), but not in the AUDIT—C scores (mode = −0.01, PP = 0.6). The SOFAS scores were higher in the HC group than in the CHR group (mode = 16.7, PP = 1.0).

3.2. GSH, CHR Status and Social and Occupational Functioning

The Spearman’s test revealed a positive association between the Social and Occupational Functioning Assessment Scale (SOFAS) scores and GSH in the CHR group (mode ρ = 0.58, posterior proportion (PP) = 0.98, Bayesian factor in favor of H1 over null H0 (BF10) = 2.1), whereas in the HC, the test indicated “absence of association” (mode ρ = 0.11, PP = 0.44, BF10 = 0.23). In the CHR group, there was neither an effect of the SOPS (mode ρ = −0.17, PP = 0.78, BF10 = 0.22) nor an effect of the CAST (mode ρ = 0.32, PP = 0.87, BF10 = 0.46) on GSH (Table 2). Note that most healthy controls reported no cannabis use, leading to limited variance in the CAST scores that could not be meaningfully related to GSH.
After accounting for the SOFAS scores, the metabolite level in the HC group was smaller than in the CHR group (mode difference = −0.27, PP = 0.97; effect size −1.06, PP = 0.97). Summary statistics of the posterior distributions of the model’s parameter estimates are reported in Table 2. Figure 2 shows the posterior distributions of the estimated between-groups difference in GSH.

3.3. Frequentist Analysis

With the positive association between the SOFAS scores and GSH already established in CHR using Bayesian analysis, our frequentist analysis was performed as a replication rather than as an addition. Median splitting of glutathione concentrations in the patient group (Figure 3) revealed a significant difference between the SOFAS scores of the low-glutathione (< 1.60 mM) and high-glutathione (> 1.60 mM) subgroups (t(11) = −2.49, p = 0.03), consistent with the Bayesian results relating GSH to the SOFAS scores in the CHR group.
In a hierarchical regression model with the GSH level as the dependent variable, with the SOFAS score entered as a predictor for all the subjects, the adjusted R2 of the model was −0.024, F(1, 41) = 0.013, p = 0.9, an insignificant effect. When the CHR status was included in the model, the R2 of the model increased to 0.07. This R2 increase was statistically significant, F(1, 40) = 5.18, p = 0.028. The regression coefficient for the CHR status was significantly negative (B = −0.56, t = −2.28), indicating that the GSH level was significantly higher in CHR subjects than in healthy controls after controlling for variance due to SOFAS. In this model, the SOFAS had a trend level of significance as a predictor, indicating that higher SOFAS scores are seen in the presence of higher GSH levels (B = 0.43, t = 1.75, p = 0.09). These results are in keeping with the Bayesian analysis reported in the manuscript.

4. Discussion

Our data provide evidence in support of a relationship between GSH levels and social and occupational functioning in clinical high-risk state and the presence of higher GSH in CHR subjects when the variance related to functional impairment is accounted for. In this study, we observed no relationship between GSH levels in the ACC and the SOFAS in the healthy volunteers, especially as the functional variability was within a narrow range among the healthy subjects. Furthermore, we did not observe any correlations between GSH and prodromal positive symptom severity. Taken together, these results support our hypothesis that GSH is a key molecular substrate underlying the functional deficits seen in the CHR state.
We report higher GSH after accounting for differences in functioning in CHR compared to healthy subjects, with reduced GSH being associated with reduced functioning among the CHR subjects. In the presence of CHR, higher levels of GSH may occur as a compensatory process across the spectrum (thus leading to a group level increase in GSH), but those who do not mount a sufficient GSH response may have poorer functional outcomes. In healthy subjects, this compensatory drive does not occur. In line with our prior work demonstrating the potentially counteracting effect of GSH on glutamatergic excitotoxicity in early stages of psychosis [54,55], we only examined the dACC in this study. GSH levels in the dACC also relate to core symptoms of psychosis such as disorganization [56], which strongly predict poor outcomes when present in CHR [57]. Nevertheless, antioxidant aberrations related to psychosis are unlikely to be limited only to the dACC, as striatal and thalamic reductions of GSH were also demonstrated in prior MRS studies [24,58].
An exciting translational utility of identifying GSH deficit in low-functioning patients is the therapeutic possibility of correcting it. A number of compounds with the potential to correct the effects of GSH deficit are in the pipeline [59,60,61]. Of these, N-acetylcysteine has been shown to improve cognition and negative symptoms in schizophrenia (six RCTs) [62] and global functioning in mood disorders [63,64]. Given the persistent nature of functional deficits, reversing them will likely require longer trials that are substantially difficult to complete. Antioxidants that increase GSH levels are more likely to benefit patients whose GSH levels are lower to begin with [65]. Our results support stratifying antioxidant trials on the basis of the baseline functional impairment or GSH levels in the future.
Our study has a number of strengths. We used the 7T-MRS sequence with improved specificity to detect GSH resonance with reduced macromolecular interference [18,66]. Among the MRS studies specifically optimized for GSH detection, 7T studies [22,67] report higher effect sizes for GSH reduction in schizophrenia compared to 3T [12,68]. We also recruited patients who had not been treated with antipsychotics and evaluated an age-, gender- and parental socio-economic status-matched control group. Nevertheless, our sample size was limited compared to the prior study addressing this issue using a 3T-MRS sequence. Furthermore, we lacked the follow-up data necessary to identify transition to psychosis among the CHR groups. From the published meta-analytical data, we expect from one to three converters in the next 2–5 years of observation [69]. The final limitation of our work was the imbalance of certain demographic characteristics, namely, CAST score differences between groups and the imbalance of gender within the CHR group.

5. Conclusions

In summary, our data offer preliminary support for the GSH level as a clinically actionable marker of prognosis in emerging adults presenting with risk features for various severe mental illnesses. The use of a longitudinal approach to track GSH levels in future CHR studies may help establish the mechanistic primacy of the antioxidant status in determining long-term outcomes.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/brainsci11070941/s1, Table S1: Mean metabolite concentration (SD) and mean CRLB (SD).

Author Contributions

Conceptualization, L.P. (Lena Palaniyappan); methodology, P.J.; validation, P.J. and R.L.; formal analysis, P.J. and R.L.; investigation, P.J., R.L., S.D.F. and L.P. (Lena Palaniyappan); resources, S.D.F., C.B., M.M., M.G. and L.P. (Laura Powe); data curation, P.J., S.D.F. and C.B.; writing—original draft preparation, P.J., R.L. and L.P. (Lena Palaniyappan); writing—review and editing, P.J., R.L., S.D.F., J.T. and L.P. (Lena Palaniyappan); visualization, P.J. and R.L.; supervision, J.T. and L.P. (Lena Palaniyappan); funding acquisition, L.P. (Lena Palaniyappan). All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by a CIHR Foundation grant (375104/2017) to L.P. (Lena Palaniyappan); an AMOSO Opportunities Fund grant to L.P. (Lena Palaniyappan); a BrainSCAN grant to R.L.; a Schulich School of Medicine and Dentistry Dean’s Scholarship to P.J.; and a Parkwood Institute Studentship to M.M. Data acquisition was supported by the Canada First Excellence Research Fund (BrainSCAN, Western University (Imaging Core)); the Innovation Fund for Academic Medical Organization of Southwest Ontario; the Bucke Family Fund; The Chrysalis Foundation; and The Arcangelo Rea Family Foundation (London, ON, Canada). The PROSPECT program is supported by a grant from Children’s Health Foundation, London Health Sciences Centre, London, ON, Canada.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Western University Research Ethics Board of University of Western Ontario (Project # 107028; October 19, 2020).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. Written informed consent has been obtained from the patients to publish this paper.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We thank Trevor Szekeres, Scott Charlton and Joseph Gati for their assistance in data acquisition and archiving. We thank Rob Bartha and Dickson Wong for consultation provided on MRS analysis and William Pavlovsky for clinical radiology consultations. We thank all the research team members of the NIMI laboratory and all the staff members of the PEPP London team for their assistance in patient recruitment and supporting clinical care. We gratefully acknowledge the participants and their family members for their contributions.

Conflicts of Interest

L.P. (Lena Palaniyappan) 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 conflict.

References

  1. McGorry, P.D.; Yung, A.R.; Phillips, L.J. The “Close-in” or Ultra High-Risk Model: A Safe and Effective Strategy for Research and Clinical Intervention in Prepsychotic Mental Disorder. Schizophr. Bull. 2003, 29, 771–790. [Google Scholar] [CrossRef] [Green Version]
  2. Yung, A.R.; McGorry, P.D. The Prodromal Phase of First-episode Psychosis: Past and Current Conceptualizations. Schizophr. Bull. 1996, 22, 353–370. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Beck, K.; Studerus, E.; Andreou, C.; Egloff, L.; Leanza, L.; Simon, A.E.; Borgwardt, S.; Riecher-Rössler, A. Clinical and functional ultra-long-term outcome of patients with a clinical high risk (CHR) for psychosis. Eur. Psychiatry 2019, 62, 30–37. [Google Scholar] [CrossRef] [PubMed]
  4. Carrión, R.E.; McLaughlin, D.; Goldberg, T.E.; Auther, A.M.; Olsen, R.H.; Olvet, D.M.; Correll, C.U.; Cornblatt, B.A. Prediction of Functional Outcome in Individuals at Clinical High Risk for Psychosis. JAMA Psychiatry 2013, 70, 1133. [Google Scholar] [CrossRef] [Green Version]
  5. Simon, A.E.; Velthorst, E.; Nieman, D.H.; Linszen, D.; Umbricht, D.; de Haan, L. Ultra high-risk state for psychosis and non-transition: A systematic review. Schizophr. Res. 2011, 132, 8–17. [Google Scholar] [CrossRef] [PubMed]
  6. Rutigliano, G.; Valmaggia, L.; Landi, P.; Frascarelli, M.; Cappucciati, M.; Sear, V.; Rocchetti, M.; De Micheli, A.; Jones, C.; Palombini, E.; et al. Persistence or recurrence of non-psychotic comorbid mental disorders associated with 6-year poor functional outcomes in patients at ultra high risk for psychosis. J. Affect. Disord. 2016, 203, 101–110. [Google Scholar] [CrossRef] [Green Version]
  7. Cannon, T.D. The Current State of the Clinical High Risk for Psychosis Research Paradigm. Biol. Psychiatry 2020, 88, 284–286. [Google Scholar] [CrossRef] [PubMed]
  8. Flatow, J.; Buckley, P.; Miller, B.J. Meta-analysis of oxidative stress in schizophrenia. Biol. Psychiatry 2013, 74. [Google Scholar] [CrossRef] [Green Version]
  9. Do, K.Q.; Trabesinger, A.H.; Kirsten-Krüger, M.; Lauer, C.J.; Dydak, U.; Hell, D.; Holsboer, F.; Boesiger, P.; Cuénod, M. Schizophrenia: Glutathione deficit in cerebrospinal fluid and prefrontal cortex in vivo. Eur. J. Neurosci. 2000, 12, 3721–3728. [Google Scholar] [CrossRef]
  10. Yao, J.K.; Leonard, S.; Reddy, R. Altered Glutathione Redox State in Schizophrenia. Dis. Markers 2006, 22, 83–93. [Google Scholar] [CrossRef] [Green Version]
  11. Gawryluk, J.W.; Wang, J.-F.; Andreazza, A.C.; Shao, L.; Young, L.T. Decreased levels of glutathione, the major brain antioxidant, in post-mortem prefrontal cortex from patients with psychiatric disorders. Int. J. Neuropsychopharmacol. 2011, 14, 123–130. [Google Scholar] [CrossRef] [Green Version]
  12. Xin, L.; Mekle, R.; Fournier, M.; Baumann, P.S.; Ferrari, C.; Alameda, L.; Jenni, R.; Lu, H.; Schaller, B.; Cuenod, M.; et al. Genetic Polymorphism Associated Prefrontal Glutathione and Its Coupling With Brain Glutamate and Peripheral Redox Status in Early Psychosis. Schizophr. Bull. 2016, 42, 1185–1196. [Google Scholar] [CrossRef] [PubMed]
  13. Kim, S.; Kang, S.; Chung, J.-H.; Park, H.; Cho, K.; Park, M.-S. Genetic Polymorphisms of Glutathione-Related Enzymes (GSTM1, GSTT1, and GSTP1) and Schizophrenia Risk: A Meta-Analysis. Int. J. Mol. Sci. 2015, 16, 19602–19611. [Google Scholar] [CrossRef] [Green Version]
  14. Gysin, R.; Kraftsik, R.; Sandell, J.; Bovet, P.; Chappuis, C.; Conus, P.; Deppen, P.; Preisig, M.; Ruiz, V.; Steullet, P.; et al. Impaired glutathione synthesis in schizophrenia: Convergent genetic and functional evidence. Proc. Natl. Acad. Sci. USA 2007, 104, 16621–16626. [Google Scholar] [CrossRef] [Green Version]
  15. Tosic, M.; Ott, J.; Barral, S.; Bovet, P.; Deppen, P.; Gheorghita, F.; Matthey, M.-L.; Parnas, J.; Preisig, M.; Saraga, M.; et al. Schizophrenia and Oxidative Stress: Glutamate Cysteine Ligase Modifier as a Susceptibility Gene. Am. J. Hum. Genet. 2006, 79, 586–592. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Matigian, N.; Abrahamsen, G.; Sutharsan, R.; Cook, A.L.; Vitale, A.M.; Nouwens, A.; Bellette, B.; An, J.; Anderson, M.; Beckhouse, A.G.; et al. Disease-specific, neurosphere-derived cells as models for brain disorders. Dis. Model. Mech. 2010, 3, 785–798. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Perkins, D.O.; Jeffries, C.D.; Do, K.Q. Potential Roles of Redox Dysregulation in the Development of Schizophrenia. Biol. Psychiatry 2020, 88, 326–336. [Google Scholar] [CrossRef]
  18. Dempster, K.; Jeon, P.; MacKinley, M.; Williamson, P.; Théberge, J.; Palaniyappan, L. Early treatment response in first episode psychosis: A 7-T magnetic resonance spectroscopic study of glutathione and glutamate. Mol. Psychiatry 2020, 25, 1640–1650. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Lambert, M.; Naber, D.; Schacht, A.; Wagner, T.; Hundemer, H.-P.; Karow, A.; Huber, C.G.; Suarez, D.; Haro, J.M.; Novick, D.; et al. Rates and predictors of remission and recovery during 3 years in 392 never-treated patients with schizophrenia. Acta Psychiatr. Scand. 2008, 118, 220–229. [Google Scholar] [CrossRef]
  20. Derks, E.M.; Fleischhacker, W.W.; Boter, H.; Peuskens, J.; Kahn, R.S. Antipsychotic Drug Treatment in First-Episode Psychosis should patients be switched to a different antipsychotic drug after 2, 4, or 6 weeks of nonresponse? J. Clin. Psychopharmacol. 2010, 30, 176–180. [Google Scholar] [CrossRef]
  21. Carbon, M.; Correll, C.U. Clinical predictors of therapeutic response to antipsychotics in schizophrenia. Dialogues Clin. Neurosci. 2014, 16, 505–524. [Google Scholar] [CrossRef] [PubMed]
  22. Kumar, J.; Liddle, E.B.; Fernandes, C.C.; Palaniyappan, L.; Hall, E.L.; Robson, S.E.; Simmonite, M.; Fiesal, J.; Katshu, M.Z.; Qureshi, A.; et al. Glutathione and glutamate in schizophrenia: A 7T MRS study. Mol. Psychiatry 2020, 25, 873–882. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Matsuzawa, D.; Obata, T.; Shirayama, Y.; Nonaka, H.; Kanazawa, Y.; Yoshitome, E.; Takanashi, J.; Matsuda, T.; Shimizu, E.; Ikehira, H.; et al. Negative Correlation between Brain Glutathione Level and Negative Symptoms in Schizophrenia: A 3T 1H-MRS Study. PLoS ONE 2008, 3, e1944. [Google Scholar] [CrossRef]
  24. Wang, A.M.; Pradhan, S.; Coughlin, J.M.; Trivedi, A.; DuBois, S.L.; Crawford, J.L.; Sedlak, T.W.; Nucifora, F.C.; Nestadt, G.; Nucifora, L.G.; et al. Assessing Brain Metabolism With 7-T Proton Magnetic Resonance Spectroscopy in Patients With First-Episode Psychosis. JAMA Psychiatry 2019, 76, 314. [Google Scholar] [CrossRef] [PubMed]
  25. Steullet, P.; Cabungcal, J.H.; Monin, A.; Dwir, D.; O’Donnell, P.; Cuenod, M.; Do, K.Q. Redox dysregulation, neuroinflammation, and NMDA receptor hypofunction: A “central hub” in schizophrenia pathophysiology? Schizophr. Res. 2016, 176, 41–51. [Google Scholar] [CrossRef] [Green Version]
  26. Maas, D.A.; Vallès, A.; Martens, G.J.M. Oxidative stress, prefrontal cortex hypomyelination and cognitive symptoms in schizophrenia. Transl. Psychiatry 2017, 7, e1171. [Google Scholar] [CrossRef] [Green Version]
  27. Worthington, M.A.; Cao, H.; Cannon, T.D. Discovery and Validation of Prediction Algorithms for Psychosis in Youths at Clinical High Risk. Biol. Psychiatry Cogn. Neurosci. Neuroimaging 2020, 5, 738–747. [Google Scholar] [CrossRef]
  28. Fornito, A.; Yung, A.R.; Wood, S.J.; Phillips, L.J.; Nelson, B.; Cotton, S.; Velakoulis, D.; McGorry, P.D.; Pantelis, C.; Yücel, M. Anatomic Abnormalities of the Anterior Cingulate Cortex Before Psychosis Onset: An MRI Study of Ultra-High-Risk Individuals. Biol. Psychiatry 2008, 64, 758–765. [Google Scholar] [CrossRef]
  29. Takayanagi, Y.; Kulason, S.; Sasabayashi, D.; Takahashi, T.; Katagiri, N.; Sakuma, A.; Obara, C.; Nakamura, M.; Kido, M.; Furuichi, A.; et al. Reduced Thickness of the Anterior Cingulate Cortex in Individuals With an At-Risk Mental State Who Later Develop Psychosis. Schizophr. Bull. 2017, 43, 907–913. [Google Scholar] [CrossRef]
  30. Pantelis, C.; Velakoulis, D.; McGorry, P.D.; Wood, S.J.; Suckling, J.; Phillips, L.J.; Yung, A.R.; Bullmore, E.T.; Brewer, W.; Soulsby, B.; et al. Neuroanatomical abnormalities before and after onset of psychosis: A cross-sectional and longitudinal MRI comparison. Lancet 2003, 361, 281–288. [Google Scholar] [CrossRef] [Green Version]
  31. Witthaus, H.; Kaufmann, C.; Bohner, G.; Özgürdal, S.; Gudlowski, Y.; Gallinat, J.; Ruhrmann, S.; Brüne, M.; Heinz, A.; Klingebiel, R.; et al. Gray matter abnormalities in subjects at ultra-high risk for schizophrenia and first-episode schizophrenic patients compared to healthy controls. Psychiatry Res. Neuroimaging 2009, 173, 163–169. [Google Scholar] [CrossRef]
  32. Whalley, H.C.; Simonotto, E.; Moorhead, W.; McIntosh, A.; Marshall, I.; Ebmeier, K.P.; Owens, D.G.C.; Goddard, N.H.; Johnstone, E.C.; Lawrie, S.M. Functional Imaging as a Predictor of Schizophrenia. Biol. Psychiatry 2006, 60, 454–462. [Google Scholar] [CrossRef]
  33. Sabb, F.W.; van Erp, T.G.M.; Hardt, M.E.; Dapretto, M.; Caplan, R.; Cannon, T.D.; Bearden, C.E. Language network dysfunction as a predictor of outcome in youth at clinical high risk for psychosis. Schizophr. Res. 2010, 116, 173–183. [Google Scholar] [CrossRef] [Green Version]
  34. Morey, R.A.; Inan, S.; Mitchell, T.V.; Perkins, D.O.; Lieberman, J.A.; Belger, A. Imaging Frontostriatal Function in Ultra-High-Risk, Early, and Chronic Schizophrenia During Executive Processing. Arch. Gen. Psychiatry 2005, 62, 254. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Andreou, C.; Borgwardt, S. Structural and functional imaging markers for susceptibility to psychosis. Mol. Psychiatry 2020, 25, 2773–2785. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Liloia, D.; Brasso, C.; Cauda, F.; Mancuso, L.; Nani, A.; Manuello, J.; Costa, T.; Duca, S.; Rocca, P. Updating and characterizing neuroanatomical markers in high-risk subjects, recently diagnosed and chronic patients with schizophrenia: A revised coordinate-based meta-analysis. Neurosci. Biobehav. Rev. 2021, 123, 83–103. [Google Scholar] [CrossRef] [PubMed]
  37. Das, T.K.; Javadzadeh, A.; Dey, A.; Sabesan, P.; Théberge, J.; Radua, J.; Palaniyappan, L. Antioxidant defense in schizophrenia and bipolar disorder: A meta-analysis of MRS studies of anterior cingulate glutathione. Prog. Neuro-Psychopharmacol. Biol. Psychiatry 2019, 91, 94–102. [Google Scholar] [CrossRef]
  38. Romeo, B.; Petillion, A.; Martelli, C.; Benyamina, A. Magnetic resonance spectroscopy studies in subjects with high risk for psychosis: A meta-analysis and review. J. Psychiatr. Res. 2020, 125, 52–65. [Google Scholar] [CrossRef]
  39. Da Silva, T.; Hafizi, S.; Andreazza, A.C.; Kiang, M.; Bagby, R.M.; Navas, E.; Laksono, I.; Truong, P.; Gerritsen, C.; Prce, I.; et al. Glutathione, the Major Redox Regulator, in the Prefrontal Cortex of Individuals at Clinical High Risk for Psychosis. Int. J. Neuropsychopharmacol. 2018, 21, 311–318. [Google Scholar] [CrossRef] [Green Version]
  40. Miller, T.J.; McGlashan, T.H.; Rosen, J.L.; Cadenhead, K.; Ventura, J.; McFarlane, W.; Perkins, D.O.; Pearlson, G.D.; Woods, S.W. Prodromal Assessment with the Structured Interview for Prodromal Syndromes and the Scale of Prodromal Symptoms: Predictive Validity, Interrater Reliability, and Training to Reliability. Schizophr. Bull. 2003, 29, 703–715. [Google Scholar] [CrossRef] [Green Version]
  41. Tkáć, I.; Gruetter, R. Methodology of1H NMR spectroscopy of the human brain at very high magnetic fields. Appl. Magn. Reson. 2005, 29, 139–157. [Google Scholar] [CrossRef]
  42. Near, J.; Edden, R.; Evans, C.J.; Paquin, R.; Harris, A.; Jezzard, P. Frequency and phase drift correction of magnetic resonance spectroscopy data by spectral registration in the time domain. Magn. Reson. Med. 2015, 73, 44–50. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Bartha, R.; Drost, D.J.; Menon, R.S.; Williamson, P.C. Spectroscopic lineshape correction by QUECC: Combined QUALITY deconvolution and eddy current correction. Magn. Reson. Med. 2000, 44, 641–645. [Google Scholar] [CrossRef]
  44. van den Boogaart, A.; Ala-Korpela, M.; Jokisaari, J.; Griffiths, J.R. Time and frequency domain analysis of NMR data compared: An application to 1D1H spectra of lipoproteins. Magn. Reson. Med. 1994, 31, 347–358. [Google Scholar] [CrossRef]
  45. Bartha, R.; Drost, D.J.; Williamson, P.C. Factors affecting the quantification of short echoin-vivo1H MR spectra: Prior knowledge, peak elimination, and filtering. NMR Biomed. 1999, 12, 205–216. [Google Scholar] [CrossRef]
  46. Wong, D. MRI Investigations of Metabolic and Structural Brain Changes in Alzheimer’s Disease and Vitamin D Deprivation. Ph.D. Thesis, The University of Western Ontario, London, ON, Canada, 2019. [Google Scholar]
  47. Rybarczyk, B. Social and Occupational Functioning Assessment Scale (SOFAS). In Encyclopedia of Clinical Neuropsychology; Springer: Cham, Switzerland, 2018; p. 3206. [Google Scholar]
  48. Stentebjerg-Olesen, M.; Jeppesen, P.; Pagsberg, A.K.; Fink-Jensen, A.; Kapoor, S.; Chekuri, R.; Carbon, M.; Al-Jadiri, A.; Kishimoto, T.; Kane, J.M.; et al. Early Nonresponse Determined by the Clinical Global Impressions Scale Predicts Poorer Outcomes in Youth with Schizophrenia Spectrum Disorders Naturalistically Treated with Second-Generation Antipsychotics. J. Child Adolesc. Psychopharmacol. 2013, 23, 665–675. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Bush, K.; Kivlahan, D.R.; McDonell, M.B.; Fihn, S.D.; Bradley, K.A. The AUDIT Alcohol Consumption Questions (AUDIT-C): An Effective Brief Screening Test for Problem Drinking. Arch. Intern. Med. 1998, 158, 1789. [Google Scholar] [CrossRef] [Green Version]
  50. Heatherton, T.F.; Kozlowski, L.T.; Frecker, R.C.; Fagerstrom, K.-O. The Fagerstrom Test for Nicotine Dependence: A revision of the Fagerstrom Tolerance Questionnaire. Addiction 1991, 86, 1119–1127. [Google Scholar] [CrossRef]
  51. Skinner, H.A. The drug abuse screening test. Addict. Behav. 1982, 7, 363–371. [Google Scholar] [CrossRef]
  52. van Doorn, J.; Ly, A.; Marsman, M.; Wagenmakers, E.-J. Bayesian rank-based hypothesis testing for the rank sum test, the signed rank test, and Spearman’s ρ . J. Appl. Stat. 2020, 47, 2984–3006. [Google Scholar] [CrossRef] [Green Version]
  53. IBM Corp. IBM SPSS Statistics for Windows, Version 26.0; IBM: Armonk, NY, USA, 2019. [Google Scholar]
  54. Limongi, R.; Jeon, P.; Théberge, J.; Palaniyappan, L. Counteracting Effects of Glutathione on the Glutamate-Driven Excitation/Inhibition Imbalance in First-Episode Schizophrenia: A 7T MRS and Dynamic Causal Modeling Study. Antioxidants 2021, 10, 75. [Google Scholar] [CrossRef]
  55. Limongi, R.; Jeon, P.; Mackinley, M.; Das, T.; Dempster, K.; Théberge, J.; Bartha, R.; Wong, D.; Palaniyappan, L. Glutamate and Dysconnection in the Salience Network: Neurochemical, Effective Connectivity, and Computational Evidence in Schizophrenia. Biol. Psychiatry 2020, 88, 273–281. [Google Scholar] [CrossRef]
  56. Pan, Y.; Dempster, K.; Jeon, P.; Théberge, J.; Khan, A.R.; Palaniyappan, L. Acute conceptual disorganization in untreated first-episode psychosis: A combined magnetic resonance spectroscopy and diffusion imaging study of the cingulum. J. Psychiatry Neurosci. 2021, 46, E337–E346. [Google Scholar] [CrossRef] [PubMed]
  57. Perkins, D.O.; Jeffries, C.D.; Cornblatt, B.A.; Woods, S.W.; Addington, J.; Bearden, C.E.; Cadenhead, K.S.; Cannon, T.D.; Heinssen, R.; Mathalon, D.H.; et al. Severity of thought disorder predicts psychosis in persons at clinical high-risk. Schizophr. Res. 2015, 169, 169–177. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Reyes-Madrigal, F.; León-Ortiz, P.; Mao, X.; Mora-Durán, R.; Shungu, D.C.; de la Fuente-Sandoval, C. Striatal Glutathione in First-episode Psychosis Patients Measured In Vivo with Proton Magnetic Resonance Spectroscopy. Arch. Med. Res. 2019, 50, 207–213. [Google Scholar] [CrossRef] [PubMed]
  59. Pereira, R.; Sousa, C.; Costa, A.; Andrade, P.; Valentão, P. Glutathione and the Antioxidant Potential of Binary Mixtures with Flavonoids: Synergisms and Antagonisms. Molecules 2013, 18, 8858–8872. [Google Scholar] [CrossRef] [Green Version]
  60. Yagishita, Y.; Gatbonton-Schwager, T.N.; McCallum, M.L.; Kensler, T.W. Current Landscape of NRF2 Biomarkers in Clinical Trials. Antioxidants 2020, 9, 716. [Google Scholar] [CrossRef]
  61. Zalachoras, I.; Hollis, F.; Ramos-Fernández, E.; Trovo, L.; Sonnay, S.; Geiser, E.; Preitner, N.; Steiner, P.; Sandi, C.; Morató, L. Therapeutic potential of glutathione-enhancers in stress-related psychopathologies. Neurosci. Biobehav. Rev. 2020, 114, 134–155. [Google Scholar] [CrossRef] [PubMed]
  62. Yolland, C.O.B.; Hanratty, D.; Neill, E.; Rossell, S.L.; Berk, M.; Dean, O.M.; Castle, D.J.; Tan, E.J.; Phillipou, A.; Harris, A.W.F.; et al. Meta-analysis of randomised controlled trials with N -acetylcysteine in the treatment of schizophrenia. Aust. New Zeal. J. Psychiatry 2020, 54, 453–466. [Google Scholar] [CrossRef]
  63. Kishi, T.; Miyake, N.; Okuya, M.; Sakuma, K.; Iwata, N. N-acetylcysteine as an adjunctive treatment for bipolar depression and major depressive disorder: A systematic review and meta-analysis of double-blind, randomized placebo-controlled trials. Psychopharmacology (Berl) 2020, 237. [Google Scholar] [CrossRef]
  64. Fernandes, B.S.; Dean, O.M.; Dodd, S.; Malhi, G.S.; Berk, M. N-acetylcysteine in depressive symptoms and functionality: A systematic review and meta-analysis. J. Clin. Psychiatry 2016, 77, e457–e466. [Google Scholar] [CrossRef] [PubMed]
  65. Conus, P.; Seidman, L.J.; Fournier, M.; Xin, L.; Cleusix, M.; Baumann, P.S.; Ferrari, C.; Cousins, A.; Alameda, L.; Gholam-Rezaee, M.; et al. N-acetylcysteine in a Double-Blind Randomized Placebo-Controlled Trial: Toward Biomarker-Guided Treatment in Early Psychosis. Schizophr. Bull. 2018, 44, 317–327. [Google Scholar] [CrossRef] [Green Version]
  66. Wong, D.; Schranz, A.L.; Bartha, R. Optimized in vivo brain glutamate measurement using long-echo-time semi-LASER at 7 T. NMR Biomed. 2018, 31, e4002. [Google Scholar] [CrossRef] [PubMed]
  67. Reid, M.A.; Salibi, N.; White, D.M.; Gawne, T.J.; Denney, T.S.; Lahti, A.C. 7T Proton Magnetic Resonance Spectroscopy of the Anterior Cingulate Cortex in First-Episode Schizophrenia. Schizophr. Bull. 2019, 45, 180–189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Hermens, D.F.; Lagopoulos, J.; Naismith, S.L.; Tobias-Webb, J.; Hickie, I.B. Distinct neurometabolic profiles are evident in the anterior cingulate of young people with major psychiatric disorders. Transl. Psychiatry 2012, 2, e110. [Google Scholar] [CrossRef] [Green Version]
  69. Fusar-Poli, P.; Bonoldi, I.; Yung, A.R.; Borgwardt, S.; Kempton, M.J.; Valmaggia, L.; Barale, F.; Caverzasi, E.; McGuire, P. Predicting psychosis: Meta-analysis of transition outcomes in individuals at high clinical risk. Arch. Gen. Psychiatry 2012, 69, 220. [Google Scholar] [CrossRef] [Green Version]
Figure 1. MRS voxel and spectra. (A) Sagittal, (B) axial and (C) coronal views of voxel positioning on the dorsal anterior cingulate cortex. (D) Sample spectra obtained from a single healthy participant. The bolded black line represents the fitted spectra with the residuals above and each individual metabolite contributions below.
Figure 1. MRS voxel and spectra. (A) Sagittal, (B) axial and (C) coronal views of voxel positioning on the dorsal anterior cingulate cortex. (D) Sample spectra obtained from a single healthy participant. The bolded black line represents the fitted spectra with the residuals above and each individual metabolite contributions below.
Brainsci 11 00941 g001
Figure 2. Bayesian analysis. Posterior distributions of the estimated between-groups difference in glutathione.
Figure 2. Bayesian analysis. Posterior distributions of the estimated between-groups difference in glutathione.
Brainsci 11 00941 g002
Figure 3. Median split analysis of glutathione (GSH) on the SOFAS. (A) Mean (± SD) glutathione concentrations (mM) of the low-GSH (< 1.60 mM; N = 7) and high-GSH (> 1.60 mM; N = 6) subgroups of clinical high-risk (CHR) cohorts as well as of the healthy controls (HC; N = 30). (B) Mean (± SD) SOFAS scores of the same three groups. Asterisks denote significant differences between the groups based on one-way ANOVA with the Tukey’s multiple comparison test, with ** denoting p < 0.05 and **** denoting p < 0.001.
Figure 3. Median split analysis of glutathione (GSH) on the SOFAS. (A) Mean (± SD) glutathione concentrations (mM) of the low-GSH (< 1.60 mM; N = 7) and high-GSH (> 1.60 mM; N = 6) subgroups of clinical high-risk (CHR) cohorts as well as of the healthy controls (HC; N = 30). (B) Mean (± SD) SOFAS scores of the same three groups. Asterisks denote significant differences between the groups based on one-way ANOVA with the Tukey’s multiple comparison test, with ** denoting p < 0.05 and **** denoting p < 0.001.
Brainsci 11 00941 g003
Table 1. Demographics and clinical characteristics.
Table 1. Demographics and clinical characteristics.
CharacteristicsClinical High Risk
(N = 13)
Healthy Controls
(N = 30)
t/χ2p
Gender (male/female)11/219/111.9480.163
Marital status (Mar/S)1/122/280.0090.926
Age (M/SD)22.2/4.021.8/3.7−0.3300.745
Parental NS-SEC (M/SD)3.54/0.883.03/1.38−1.4440.158
Cannabis use (Y/N)6/71/294.9260.026
Antidepressant use (Y/N)3/10-
Benzodiazepine use (Y/N)2/11-
Antipsychotic use (Y/N)1/12-
SOPS, total (M/SD)8.4/5.0-
CAST, total (M/SD)11.4/7.86.2/0.8−2.200.052
Audit—C, total (M/SD)5.5/3.85.9/2.50.2810.783
SOFAS (M/SD)
APS/BLIPS
Converted (Y/N)
67.7/9.5
13/0
3/10
82.7/3.8

6.267

0.000

Chi-squared analyses for categorical variables and independent t-tests for continuous variables were used to calculate p-values for differences between groups. Mar—married, S—single, M—mean, SD—standard deviation, NS-SEC—National Statistics Socio-economic Classification, Y—yes, N—no, SOPS—scale of prodromal symptoms, CAST—Cannabis Abuse Screening Test, AUDIT—C—Alcohol Use Disorders Identification Test—Concise, SOFAS—Social and Occupational Functioning Assessment Scale, APS—attenuated psychotic syndrome, BLIPS—brief and limited intermittent psychosis.
Table 2. Parameter estimates (posteriors) of the hierarchical Bayesian linear model and the Spearman’s correlation test.
Table 2. Parameter estimates (posteriors) of the hierarchical Bayesian linear model and the Spearman’s correlation test.
ParameterMeanMedianModeHDIlowHDIhigh
β 0 0.8490.8420.823−0.1651.875
β[CHR]0.1350.1350.135−0.0130.274
β[HC]−0.135−0.135−0.135−0.2740.013
βSOFAS0.0090.0100.010−0.0040.023
σβ0.4400.3350.1990.0001.155
σi0.2700.2670.2640.2120.332
ρ [HC]SOFAS0.1540.1430.114−01260.378
ρ [CHR]SOFAS0.4780.5090.5860.0450.769
ρ [CHR]SOPS−0.073−0.076−0.171−0.3350.196
ρ [CHR]CAST0.2010.2070.315−0.1240.500
Note: HDI represents 95% of the most credible values. HDI—highest density interval, β 0 —intercept, β[CHR]—deflection parameter for the CHR group, β[HC]—deflection parameter for the HC group, βSOFAS—deflection parameter for the SOFAS covariate, σβ—standard deviation of the baseline parameter, σi—standard deviation of the predicted value, ρ [HC]SOFAS—Spearman’s correlation between the SOFAS score and GSH of the control group, ρ [CHR]SOFAS—Spearman’s correlation between the SOFAS score and GSH of the CHR group, ρ [CHR]SOPS—Spearman’s correlation between the SOPS score and GSH of the CHR group, ρ [CHR]CAST—Spearman’s correlation between the CAST score and GSH of the CHR group.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jeon, P.; Limongi, R.; Ford, S.D.; Branco, C.; Mackinley, M.; Gupta, M.; Powe, L.; Théberge, J.; Palaniyappan, L. Glutathione as a Molecular Marker of Functional Impairment in Patients with At-Risk Mental State: 7-Tesla 1H-MRS Study. Brain Sci. 2021, 11, 941. https://0-doi-org.brum.beds.ac.uk/10.3390/brainsci11070941

AMA Style

Jeon P, Limongi R, Ford SD, Branco C, Mackinley M, Gupta M, Powe L, Théberge J, Palaniyappan L. Glutathione as a Molecular Marker of Functional Impairment in Patients with At-Risk Mental State: 7-Tesla 1H-MRS Study. Brain Sciences. 2021; 11(7):941. https://0-doi-org.brum.beds.ac.uk/10.3390/brainsci11070941

Chicago/Turabian Style

Jeon, Peter, Roberto Limongi, Sabrina D. Ford, Cassandra Branco, Michael Mackinley, Maya Gupta, Laura Powe, Jean Théberge, and Lena Palaniyappan. 2021. "Glutathione as a Molecular Marker of Functional Impairment in Patients with At-Risk Mental State: 7-Tesla 1H-MRS Study" Brain Sciences 11, no. 7: 941. https://0-doi-org.brum.beds.ac.uk/10.3390/brainsci11070941

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop