Abstract
Background: Disorganized thinking is a core feature of acute psychotic episodes that is linked to social and vocational functioning. Several lines of evidence implicate disrupted cognitive control, excitatory overdrive and oxidative stress relating to the anterior cingulate cortex as mechanisms of conceptual disorganization (CD). We examined 3 candidate mechanistic markers related to CD in first-episode psychosis: glutamate excess, cortical antioxidant (glutathione) status and the integrity of the cingulum bundle that connects regions implicated in cognitive control.
Methods: We used fractional anisotropy maps from 7 T diffusion-weighted imaging to investigate the bilateral cingulum based on a probabilistic white matter atlas. We compared high CD, low CD and healthy control groups and performed probabilistic fibre tracking from the identified clusters (regions of interest within the cingulum) to the rest of the brain. We quantified glutamate and glutathione using magnetic resonance spectroscopy (MRS) in the dorsal anterior cingulate cortex.
Results: We found a significant fractional anisotropy reduction in a cluster in the left cingulum in the high CD group compared to the low CD group (Cohen’s d = 1.39; p < 0.001) and controls (Cohen’s d = 0.86; p = 0.009). Glutamate levels did not vary among groups, but glutathione levels were higher in the high CD group than in the low CD group. We also found higher glutathione related to lower fractional anisotropy in the cingulum cluster in the high CD group.
Limitations: The MRS measures of glutamine were highly uncertain, and MRS was acquired from a single voxel only.
Conclusion: Acute CD relates to indicators of oxidative stress, as well as reduced white matter integrity of the cingulum, but not to MRI-based glutamatergic excess. We propose that both oxidative imbalance and structural dysconnectivity underlie acute disorganization.
Introduction
Conceptual disorganization (CD) refers to difficulties in the goal-directed sequencing of thoughts that manifest as circumstantial, illogical or tangential speech or weakened goal of thinking (loose associations) in patients with psychotic disorders such as schizophrenia.1 Among patients with schizophrenia, CD is an important contributor to reduced real-world functioning.2 It is closely linked to the construct of formal thought disorder (FTD), but unlike FTD, CD is a clinical interpretation that is not restricted to the speech domain and refers specifically to the positive (rather than the negative) dimension of disorganized thinking.
Based on the idea that CD relates to the aberrant shunting of cognitive processes,3,4 structural disconnection among brain regions has been investigated as the neural basis using diffusion tensor imaging (DTI). Viher and colleagues5 and Cavelti and colleagues6 found no association between white matter integrity and the severity of disorganization or FTD, but Bopp and colleagues7 noted more severe FTD in patients with reduced integrity of the cingulum tract. More recently, de Boer and colleagues8 found an association between computational linguistic indices relating to negative symptoms and white matter integrity in schizophrenia. This effect was stronger for the entire white matter structure, rather than for language-related tracts.8 These previous DTI studies were limited to patients with established illness (> 2 years) and with long exposure to antipsychotic medications that affected both DTI indices9 (but see Voineskos and colleagues10) and symptom severity scores.6,8 At present, the structural basis of acute CD in minimally treated patients in first-episode psychosis (FEP) is unknown.11,12
An experimental approach that provokes symptoms of acute CD is subanesthetic ketamine infusion in both healthy volunteers and patients with schizophrenia.13 Ketamine disrupts inhibitory interneuron function, resulting in the disinhibition of the excitatory glutamatergic outputs of the prefrontal cortex. Sustained disinhibition increases extracellular glutamate, resulting in oxidative stress.14 Glutathione (GSH), a neuronal antioxidant, is the primary defence against excitotoxicity in the human brain.15 Subgroups of patients with psychoses and higher levels of CD show elevated peripheral markers of oxidative stress.16–18 Reduced GSH and glutamate have been reported in the anterior cingulate cortex (ACC) in several magnetic resonance spectroscopy (MRS) studies in schizophrenia.19,20 However, despite the relevance of the ketamine model to CD,11 the role of brain glutamate and GSH in CD remains unstudied to date.
In this study, we focused on the cingulum (a frontal association bundle) and the dorsal ACC related to acute CD for 3 important reasons. First, the structural integrity of the cingulum has repeatedly been observed to be abnormal in psychosis,7,21,22 and this disruption has been linked to the positive symptoms of schizophrenia.23–26 Second, cognitive control function is presumed to be disrupted in patients with high CD.27 Cingulum integrity underlies several aspects of the cognitive control dysfunction seen in schizophrenia.28,29 The cingulum is the primary association tract that connects the dorsal ACC with other sites that have been implicated in cognitive control, such as the precuneus.30 Third, in previous structural studies we reported volumetric abnormalities in networks anchored on the dorsal ACC in relation to thought disorder.31 In particular, disorganized thinking was associated with structural variations in the cognitive control regions, more than in language regions.32 Abnormalities in the structural integrity of the cingulum have been linked to CD by Bopp and colleagues.7 Finally, given its protracted development compared to other long bundles,33 the cingulum is considered to be particularly vulnerable to maturational aberrations relevant for schizophrenia. Markers of oxidative stress, such as reduced GSH in the ACC, relate to the maturational integrity of the cingulum bundle.34
In the current study, we bring together 2 candidate mechanistic markers related to acute CD in FEP: glutamate excess and antioxidant (GSH) status in the ACC using MRS; and the structural integrity of the cingulum based on fractional anisotropy (FA), using DTI. We also investigated the relationship between these mechanistic features in patients with untreated FEP and a median antipsychotic medication exposure of less than 3 days. This study is part of an ongoing preregistered observational trial (NCT02882204) investigating the pathophysiology of FTD in the early stages of schizophrenia with a focus on medial prefrontal connectivity, glutamatergic excess and GSH abnormalities.
Methods
Participants
Patients with FEP (31 with schizophrenia; 9 with schizoaffective/schizophreniform psychosis) between the ages of 15 and 39 years were recruited from the Prevention and Early Intervention for Psychosis Program at London Health Sciences Centre in London, Ontario, Canada, between April 2017 and July 2019. The program uses an assertive case-management model to provide assessment and treatment to individuals experiencing FEP.
In total, all patients and 25 healthy controls had scans of adequate quality (without obvious artifact and visible structural damage) for DTI analysis. Patients included in the study were experiencing FEP and had received antipsychotic treatment for less than 14 days. Both inpatients and outpatients were eligible to participate if they were able to provide informed consent and safely participate in the MRI protocol. All participants received a consensus diagnosis from 3 psychiatrists (L.P./K.D. and the primary treatment provider) after approximately 6 months, on the basis of the best estimate procedure (as described by Leckman and colleagues35) and the Structured Clinical Interview for DSM-5.36 After the 6-month consensus diagnosis, participants who met the criteria for bipolar disorder with psychotic features, major depressive disorder with psychotic features or suspected drug-induced psychoses were excluded from further analyses. Care was provided as usual to the study participants through their psychiatrist and other allied health members of the Prevention and Early Intervention for Psychosis Program. Antipsychotic medications were chosen by the treating psychiatrist in collaboration with the patient and/or their substitute decision-maker. Patients were offered the option of treatment with a long-acting injectable at the earliest opportunity,37 in accordance with current national guidelines for the treatment of FEP.
Recruitment of healthy controls (n = 25) was carried out through posters that advertised the opportunity to participate in a neuroimaging study to track the outcomes of FEP. Healthy controls had no personal history of mental illness and no first-degree relatives with a history of psychotic disorders. Group matching with the FEP cohort for age, sex and parental education was maintained.
Exclusion criteria for both the FEP and healthy control groups were as follows: meeting the criteria for a substance use disorder in the past year, according to the DSM-5;36 a history of a major head injury (leading to a substantial period of unconsciousness or seizures); a substantial, uncontrolled medical illness; or any MRI scan contraindications.
All potential study participants provided written, informed consent before participation, as outlined by the Western University Health Sciences Research Ethics Board. The sample reported here partially overlaps with those of 2 of our previous studies.38,39
Symptom assessment
We assessed the severity of symptoms in patients using the 8-item Positive and Negative Syndrome Scale (PANSS-8),40 which measures the severity of positive and negative symptoms. We also used the Social and Occupational Functioning Assessment Scale (SOFAS) to assess the overall level of functioning. 41 We used a modified digit symbol substitution task (DSST) to quantify processing speed, a cognitive function that shows the most prominent reduction among patients with psychosis.42 The written and oral items on the DSST were scored separately, and we used the mean scores from these items for assessment, as in our previous study.43
We used ultra-high-resolution (7 T) diffusion-weighted imaging in 40 patients with FEP and 25 controls to investigate the bilateral cingulum, defined according to the Johns Hopkins University white matter atlas.44 Based on the CD scores from the PANSS-8, we divided the patients into 2 groups: low CD (P2 = 1–3) and high CD (P2 = 4–6; see Appendix 1, available at jpn.ca/200167-a1, for more on this rationale). We chose this cut-off to distinguish between those with clinically severe symptoms and those with minimal symptoms, as recommended by the Remission Working Group.45 This cut-off has also recently been employed to distinguish patients with or without FTD and language dysfunction. 46 These 2 groups were not different in terms of the severity of other positive or negative symptoms. We applied tract-based spatial statistics (TBSS) to compare FA values using FSL software47 with 5000-permutation correction within the selected tract. We used F tests for group comparisons.
MRI data acquisition
We acquired all magnetic resonance images on a 7.0 T Siemens Magnetom Terra MRI scanner using an 8-channel transmit/32-channel receive head coil at the Centre for Functional and Metabolic Mapping, Robarts Research Institute, Western University, London, Ontario, Canada. We acquired DTI using 64 diffusion-weighted volumes (b value = 1000 s/mm2) and an unweighted volume, with a 2-dimensional, spin-echo, echo-planar imaging sequence (matrix size = 104 × 104; 72 slices; voxel resolution 2 × 2 × 2 mm; repetition time 5100 ms; echo time 50.2 ms; generalized autocalibrating partial parallel acquisition 3; multi-band 3; phase partial Fourier 6/8).
MRI preprocessing
We processed diffusion-weighted images using an open-source, containerized application called “prepdwi” (www.github.com/khanlab/prepdwi), which uses the BIDS48 and BIDS Apps49 standards to perform standardized preprocessing and fitting. We first preprocessed diffusion-weighted MRI data with denoising using a local PCA method with “dwidenoise” from mrtrix50 and correction of ringing artifacts with the “unring” tool.51 We corrected eddy current and susceptibility-based distortions using “topup”52 and “eddy”53 based on FSL, with the “–repol” option enabled for outlier replacement.54 We performed within-subject rigid registration of the corrected diffusion-weighted imaging volume and the T1-weighted volume using block-matching55 to bring the diffusion-weighted images into the same space as the T1-weighted images. We corrected gradient nonlinearities using vendor-provided spherical harmonic coefficients for the gradient system (AC84), applied with the “gradient_unwarp” tool56 to generate a nonlinear transformation, which was composed with the T1-weighted linear transformation to resample the diffusion-weighted images into the corrected T1-weighted space in a single step. We used modulation with the determinant of the Jacobian matrix of the unwarping to correct for intensity differences in the magnitude images because of gradient nonlinearities. We then used preprocessed diffusion-weighted images in the T1-weighted space to estimate diffusion tensor metrics using “dtifit.”57 All participants’ FA images were first nonlinearly aligned to an FMRIB58_FA standard brain image using nonlinear registration and then used to create a mean FA skeleton representing the centres of all tracts common to all participants (threshold of 0.2). After that, each participant’s aligned FA was projected onto this skeleton and fed into voxel-wise comparison. The voxels with significant group differences were labelled using the Johns Hopkins University white matter atlas44 and enlarged using the “TBSS_fill” tool (https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/TBSS/UserGuide?highlight=%28TBSS_fill%29) for illustrative purposes. In addition to FA, we undertook post hoc comparisons of radial diffusivity (reflecting myelin content), axial diffusivity (reflecting axonal integrity) and mode of anisotropy (reflecting fibre organization) for the regions that showed significant FA abnormalities.
We performed fibre tracking analysis using Bedpostx (based on FSL) for local modelling of diffusion parameters and the Probtrackx (https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/FDT) for tractography and connectivity-based segmentation (Appendix 1, Figure S1). Probabilistic fibre tracking is performed from the seed region that shows a difference in the group comparison to the whole brain of fibres in the native space, with parameters as follows: way total number = 70 000; number of samples (number of individual pathways that are drawn through the probability distributions on principle fibre direction) = 5000; curvature threshold = 0.2; maximum number of steps = 2000; step length = 0.5 mm; fibthresh (subsidiary fibre volume threshold) = 0.01. We performed 2-stage registration, initially with FLIRT for linear registration and then with FNIRT for nonlinear registration, both using FSL. We extracted the mean FA value and the volume of the traced fibres from each participant for further analysis.
Metabolites assessment
As described in our previous study,58 we measured metabolites (glutamate, GSH, glutamine) using a Siemens Magnetom 7.0 T MRI (Siemens Healthineers) using an 8-channel transmit/32-channel receive, head-only radiofrequency coil at the Centre for Metabolic Mapping at Western University. A 2.0 × 2.0 × 2.0 cm (8 cm3) proton magnetic resonance spectroscopy (1H-MRS) voxel was placed in the bilateral dorsal ACC using a 2-dimensional anatomic imaging sequence in the sagittal direction (37 slices; repetition time 8000 ms; echo time 70 ms; flip angle [α] 120°; thickness 3.5 mm; field of view 240 × 191 mm). We acquired a total of 32 channel-combined, water-suppressed spectral averages using the semi-LASER 1H-MRS pulse sequence (repetition time 7500 ms, echo time 100 ms) during each scan session; participants were at rest and asked to stare at a white cross on a screen for 4 minutes. Water suppression was achieved using the VAPOUR preparation sequence, and water-unsuppressed spectra were acquired for spectral postprocessing (lineshape deconvolution) and quantification (water amplitude as concentration reference). Each of the 32 spectral acquisitions was corrected for frequency and phase drifts as described in Near et al59 before averaging. Metabolite concentrations were corrected for grey and white matter volumes using previously described methods60 in Barstool.61 Water-subtracted spectra were modelled using the fitMAN a priori–based minimization algorithm and a quantification template that included 17 metabolite spectral signatures derived from simulation.61 Importantly, we used a long echo time acquisition, optimized for GSH and glutamate signals (in line with the trial preregistration) with minimal macromolecular interference, as per Wong and colleagues.61 We inspected all spectra visually for quality, and included only metabolites with Cramer–Rao lower bounds (CRLB) of less than 10% in the analyses. An example of representation of fit is provided in Figure 1. For glutamine, the mean CRLB of over 20% for both FEP and healthy control groups did not meet our rigorous CRLB cut-off of less than 10%, indicating higher uncertainty in its quantification compared with glutamate and GSH spectra. We have reported the results for glutamine in Appendix 1.
Example of representation of fit for metabolites. GABA = γ-aminobutyric acid.
We previously examined the overall resting functional MRI (fMRI) correlates of CD from an overlapping sample; this has been reported elsewhere.62
Statistical analyses
We carried out statistical analyses of participant demographic and clinical characteristics and parameters of imaging using SPSS 22.0 (IBM). We made group comparisons of the variables using independent-samples t tests or F tests for continuous variables; Kruskal–Wallis H tests for group comparison of ordinal items; or χ2 tests for proportions, as appropriate. We performed voxel-wise group comparisons and within-group comparisons on the mean FA skeleton in the framework of the generalized linear model using the Randomize tool in FSL.63 All analyses were conducted with 5000 permutations. We set the significance level at p < 0.05 with threshold-free cluster enhancement correction64 and family-wise error rate correction for multiple comparisons.
Results
Demographic and clinical characteristics
All patients had acute psychosis, were minimally treated and had significant social and occupational dysfunction at the time of scanning. The demographic and clinical characteristics of the sample are presented in Table 1. The 2 patient groups did not differ in terms of age, sex, socioeconomic status, duration of untreated illness in days, exposure to antipsychotic medication in days, total daily dose of antipsychotic medication or number of patients with schizophrenia in either group (diagnosed 6 months after the initial presentation). Across the 3 groups, we observed significant differences in SOFAS and DSST scores, confirmed by analysis of variance.
Demographic and clinical characteristics*
Group differences in white matter integrity
In the voxel-wise F test, we observed a cluster (MNI coordinates x, y, z = −8, −1, 34) with a significant difference in FA (p = 0.036, corrected) among the 3 groups in the left cingulum (Figure 2A). We extracted mean FA, axial diffusivity, radial diffusivity and mode of anisotropy value in this cluster (Figure 2B). The high CD group showed lower FA than the low CD group (Cohen’s d = 1.39, p < 0.001) and the healthy controls (Cohen’s d = 0.86, p = 0.009). Moreover, the high CD group had lower axial diffusivity than the low CD group (Cohen’s d = 0.84, p = 0.021). Interestingly, the low CD group had significantly higher axial diffusivity than the healthy controls (Cohen’s d = 0.63, p = 0.025), indicating an aberrant increase in the axonal integrity of the cingulum in patients with lower levels of CD. We did not find any significant difference in radial diffusivity or mode of anisotropy value among the 3 groups.
A small part of the cingulum tract showed reduced FA and AD in the group with high conceptual disorganization. (A) Results of F tests of tract-based spatial statistics among 3 groups in the bilateral cingulum tract with permutation testing (n = 5000). (B) Values for FA, RD, AD, MO, total volume of fibres and FA of fibres tracked from the centroid of the seed region, identified in tract-based spatial statistics to the rest of the brain. AD = axial diffusivity; FA = fractional anisotropy; HC = healthy controls; MO = mode of anisotropy; RD = radial diffusivity.
Upon fibre tracking from the seed region to the whole brain, the high CD group showed significantly less volume of all fibres than the low CD group (Cohen’s d = 0.69, p = 0.030) and healthy control group (Cohen’s d = 0.81, p = 0.019). The mean FA of the tracked fibres did not differ among the 3 groups, indicating localized changes to the cingulum bundle (independent t tests, p = 0.11 to 0.76), although the mean FA of the tracked fibres was lower in patients with lower DSST scores (r = 0.381, p = 0.015). We found no significant relationship between FA and DSST in the healthy control group (r = −0.152, p = 0.47). The volume of all fibres was negatively correlated with P2 value (r = −0.264, p = 0.034).
Group differences in MRS metabolite measures
The F tests showed no significant difference in glutamate levels among the 3 groups, but a significant difference in GSH levels (F = 3.894, p= 0.027; Figure 3). The high CD group had higher GSH levels than the healthy controls (Cohen’s d = 1.096, p = 0.008). Given the suspected differences in the neuronal versus astrocytic origin of the glutamate and GSH signals, respectively (see Duarte and Xin65 and Harris and colleagues66 for a review), we calculated the glutamate–GSH ratio as a proxy indicator of relative neuronal density. We found that the high CD group had a significantly lower ratio (driven by higher GSH concentration) than the low CD group (Cohen’s d = 0.912, p = 0.011) and the healthy controls (Cohen’s d = 1.663, p = 0.001) indicating a possible increase in metabolite signals of astrocytic origin (GSH) compared to those of neuronal origin (glutamate) in relation to disorganization.
Probabilistic map of magnetic resonance spectroscopy cingulate cortex voxel placement. Colour bar indicates the probability of the placed voxel across all participants (0 to 1).
Relationship between cingulum microstructure and MRS glutamate and GSH
Axial diffusivity of the cingulum cluster was negatively related to glutamate (r16 = −0.586, p = 0.009) and GSH (r16 = −0.523, p = 0.026) in healthy controls, but only to GSH (r9 = −0.797, p = 0.003) in the high CD group. We found no DTI-MRS correlations in the low CD group (all p > 0.1).
To determine if the observed relationship between DTI metrics, GSH and P2 score on the PANSS-8 scale (conceptual disorganization) was driven by other symptom differences or antipsychotic exposure status, we related GSH concentration, axial diffusivity and fibre volume to negative symptom burden, positive symptom burden, total defined daily dose and total medication days. None of these correlations was significant (r = −0.16 to 0.2; p = 0.27 to 0.98). Cumulative medication exposure (defined daily dose × days) did not differ between the low CD and high CD groups (t = −0.178, p = 0.86); we also found no differences between these groups in terms of the total number of antipsychotic exposure days (t = 1.486, p = 0.15). To determine if the observed relationship between DTI metrics, GSH and P2 score was an artifact of binarization, we ran a Spearman’s correlation between GSH concentration, axial diffusivity, and fibre volume and P2 scores in all participants. We noted that P2 scores were negatively correlated to the volume of fibres from the seed region to whole brain (r = −0.264, p = 0.034), but we found no significant correlation between P2 score, axial diffusivity and GSH (p > 0.1).
Discussion
To our knowledge, this is the first study to link CD to markers of excitotoxicity (glutamate), oxidative stress (GSH) and white matter integrity in the cingulum in FEP. We report 3 major observations: the presence of severe acute CD in FEP is associated with reduced white matter integrity of the cingulum, characterized by reduced axial diffusivity and decreased volume of fibres in the cingulum tract; it is also associated with increased GSH concentration in the dorsal ACC; and it is not associated with increased glutamate concentration. In addition, higher cingulum integrity was related to lower glutamate and GSH levels in healthy controls. Our findings suggest that the severity of CD in untreated FEP may be related to markers of oxidative stress and white matter damage, but not to glutamate deficits as indexed using MRS.
Our observation of reduced white matter integrity of the cingulum in patients with high CD was consistent with the findings of Bopp and colleagues,7 who systematically studied white matter integrity in relation to persistent disorganization. Reduced white matter integrity in the cingulum is a well-replicated finding in psychotic disorders,22,67–69 although heterogeneous changes are often reported at the whole brain level.70 The cingulum is a large body of fibre connecting multiple midline brain regions that play a crucial role in self-related processing (precuneus, ventromedial prefrontal cortex and medial temporal cortex).71 Individual differences in cingulum microstructure be linked to cognitive control measures in healthy individuals,30 as well as in patients with schizophrenia.72,73 Disrupted connectivity of the cingulum reported here supports the notion that the severity of CD relates to disrupted domain-general control aspects of information processing in schizophrenia.32
We found increased GSH levels in patients with a high degree of CD in acute FEP, despite the normal levels of glutamate in the ACC. This association was specific to CD, as neither negative symptom burden nor positive symptom severity was related to GSH concentrations in this sample. In an overlapping sample, we previously noted that people with FEP and higher GSH were likely to respond more quickly to treatment, despite their higher symptom burden. Interestingly, among adolescents at risk for psychosis, increased GSH levels in the ACC have been linked to a higher degree of disorganization.74 Taken together, these results suggest that an acute oxidative stress response associated with disrupted white matter integrity may contribute to CD in the early stages of psychosis, but higher GSH may indicate an intact antioxidant response mounted to neutralize the oxidative stress. Patients who mount this GSH response may ultimately experience a better treatment response, and those with a persistent illness profile may show a significant, albeit small, effect-size reduction in GSH.19,75 An alternative possibility is that higher GSH itself may be an index of excitotoxicity, because GSH is known to act as a reservoir for glutamate when normal routes of glutamate synthesis are impaired. 76 Our recent fMRI-based causal modelling did not support this idea, instead favouring a protective role for GSH against glutamatergic excitotoxicity.77 Finally, a reduction in glutamate–GSH ratio reported here may indicate acute astrogliosis and/or rarefication of glutamatergic synapses. These mechanistic underpinnings could not be resolved with our current data; pharmaco-fMRI studies manipulating GSH levels will be critical for future insights.
Limitations
Our study had notable strengths, including the use of highly resolved 7 T spectroscopy (which allowed for concurrent GSH and glutamate quantification) a combined DTI-MRS protocol and the recruitment of minimally medicated patients with FEP.
Several limitations must be considered as well. First, because we optimized our acquisition for concurrent measurement of glutamate and GSH, glutamine spectra were estimated with higher uncertainty (Appendix 1). Second, we relied on a single voxel placed at the dorsal ACC, in line with other 7 T MRS studies in psychosis;78,79 the status of GSH and glutamate in relation to CD may differ in other brain regions. We also restricted our DTI analysis to the cingulum, given the previous work of Bopp and colleagues.7 A TBSS analysis across the entire white matter skeleton (reported in Appendix 1) revealed no significant group differences. Although we had sufficient power to detect a group difference between the high and low CD groups in a hypothesis-driven search restricted to a single tract, our sample was likely of insufficient size for a whole brain mass univariate search with necessary corrections for multiple testing.
Conclusion
The neural correlates of acute conceptual disorganization likely relate to cortical oxidative stress and disrupted white matter integrity. Given the functional importance of persistent disorganization, we call for clinical trials of antioxidants and agents that improve neuroglial interactions to focus on addressing the burden of disorganization in psychosis.
Footnotes
Competing interests: L. Palaniyappan reports grants and personal fees from Otsuka Canada, grants and personal fees from Janssen Canada, grants from Sunovion, personal fees from SPMM Course (United Kingdom), other fees from Oxford University Press and personal fees from the Canadian Psychiatric Association, all outside the submitted work. He is an associate editor of JPN, but was not involved in the decision to accept this manuscript for publication. No other competing interests were declared.
Contributors: L. Palaniyappan designed the study. K. Dempster and J. Théberge acquired the data, which Y. Pan, P. Jeon and A. Khan analyzed. Y. Pan and L. Palaniyappan wrote the article, which K. Dempster, P. Jeon, J. Théberge and A. Khan reviewed. All authors approved the final version to be published and can certify that no other individuals not listed as authors have made substantial contributions to the paper.
Funding: This study was funded by CIHR Foundation Grant No. 375104/2017 (to LP); Schulich School of Medicine Clinical Investigator Fellowship (to KD); graduate student salary support of PJ from NSERC Discovery Grant No. RGPIN2016-05055 (to JT); Canada Graduate Scholarship (to KD); the National Natural Science Foundation of China, grant numbers 81671335, 81701325, 81801353 (to YP). 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).
- Received August 29, 2020.
- Revision received December 23, 2020.
- Accepted January 19, 2021.
This is an Open Access article distributed in accordance with the terms of the Creative Commons Attribution (CC BY-NC-ND 4.0) licence, which permits use, distribution and reproduction in any medium, provided that the original publication is properly cited, the use is noncommercial (i.e., research or educational use), and no modifications or adaptations are made. See: https://creativecommons.org/licenses/by-nc-nd/4.0/