Participants
The ARMS data for the present study were collected as a part of the Minds in Transition (MinT) project [32], which is a longitudinal study focused on the transition from ARMS to schizophrenia. The research was conducted in collaboration with early psychosis services located in metropolitan, regional, and rural centres across New South Wales, Australia. Participant referrals were obtained from a variety of sources, including the national Headspace initiative (https://headspace.org.au), mental health workers, general practitioners, school counsellors, and self-referrals.
The original MinT study recruited 102 ARMS individuals and 61 healthy control (HC) participants. In the present study, we analysed structural brain imaging data available from a subset of the MinT study, specifically 44 ARMS individuals aged 16 years and older (mean age 19.7, SD 2.1, range 16.2 – 24.1 years; 21 males and 23 females). Additionally, we included data from 19 schizophrenia patients (SCZ) who met the DSM-IV diagnostic criteria and were younger than 25 years of age (mean age 22.6, SD 1.5, range 19.9 – 24.8 years; 12 males and 7 females; mean age ARMS < SCZ: p < .001), which was obtained from the Australian Schizophrenia Research Bank (ASRB [33]). Lastly, we also analysed data from 36 HC participants (mean age 21.1, SD 2.0, range 16.6 – 24.8 years; 16 males and 20 females) pooled from both the MinT study (n = 17) and the ASRB (n = 19). The selection of participants from the ASRB was based on their age being less than 25 years.
ARMS was assessed with the Comprehensive Assessment of At-risk Mental State (CAARMS; version Yung et al. [34]). CAARMS defines ARMS as a significant decline of functioning over a one-year period, indicated by a drop of at least 30 points on the Global Assessment of Function (GAF) rating scale [35]. This decline of functioning is accompanied by (i) emerging, attenuated psychotic symptoms and/or, (ii) brief limited intermittent psychotic symptoms and/or (iii) an immediate family history of schizophrenia. In this study, the ARMS group was further divided into two subgroups based on the median split derived from the CAARMS composite score. The composite score was developed through expert consultation aiming to capture at-risk mental state as broadly as possible, thus avoiding over-reliance on individual symptom domains. Therefore, the composite score provides a better representation of the psychopathology observed in ARMS. The composite score was calculated by summing the intensity rating scores for unusual thought content, non-bizarre ideas, perceptual abnormalities, disorganized speech, alogia, avolition/apathy, anhedonia, social isolation, impaired role function, disorganising/odd/stigmatising behaviour, aggression/dangerous behaviour, mania, depression, mood swings/liability, and anxiety. The two subgroups consisted of 22 subjects with low and 22 subjects with high at-risk symptom ratings (Table 4). The Social and Occupational Functioning Assessment Scale (SOFAS [35]) was also employed in the assessment.
For gender and age matching, we selected the best matched HC participants for comparison with the two respective clinical groups. We included 29 HC individuals (mean age 20.4, SD 1.5, range 16.6 – 22.8 years; 13 males and 16 females) for comparison with the ARMS groups and 26 HC participants (mean age 22.0, SD 1.5, range 19.7 – 24.8 years; 12 males and 14 females) for comparison with the SCZ group (Tables 4 and 5).
Insert Tables 4 and 5 about here
Exclusion criteria for the ARMS participants included pre-existing psychosis individuals whose symptoms exceeded the CAARMS psychosis threshold and individuals receiving antipsychotic pharmacotherapy. Substance use was assessed with the Alcohol Use Disorders Identification Test (AUDIT [36] , the Cannabis Use Disorders Identification Test (CUDIT [37]), and the Opiate Treatment Index: drug use all types (OTI [38]). Participants diagnosed with drug dependence, as assessed by either the Structured Clinical Interview for DSM-IV Axis I Disorders (Clinical Version; SCID-CV) or the Kiddie Schedule for Affective Disorders and Schizophrenia for School-aged Children, Present and Lifetime Version (K-SADS-PL), were also excluded. Additionally, participants with a history of head injury causing loss of consciousness for more than 15 min, organic brain impairment, estimated pre-morbid IQ lower than 70, impaired hearing (>20 dB [SPL]), history of nasal trauma, or those meeting MRI exclusion criteria were also excluded from the study.
Study protocol
Upon study entry, all ARMS participants undertook a battery of clinical and neuropsychological tests and electroencephalographic recordings over the course of 2 to 3 days (reported in [32]). Participants were also given the opportunity to participate to undergo MRI brain scans. For the first year of the study, ARMS participants were contacted every three months to assess their clinical status. At the 12-month follow-up, potential transition to psychosis was assessed by applying a DSM-IV diagnosis, using either the SCID-CV or the K-SADS-PL.
MRI data acquisition
All MRI data used in this study were collected with 1.5 T Siemens Avanto MRI scanners. The acquisition protocol was consistent across the five participating sites, including the MinT and ASRB projects, which were conducted concurrently. The T1-weighted magnetisation-prepared rapid-acquisition gradient echo sequence used by all five sites employed the following parameters: a repetition time of 1980 ms, an echo time of 4.3 ms, a voxel size of .9765625 x .9765625 x 1mm3, and a flip angle of 15º.
Image processing
The software Freesurfer 5.1 [39-42] was used to estimate the cortical thickness, surface area of the grey/white matter interface, and intracranial volume (ICV). In order to ensure data quality, we implemented a quality control process, which involved an iterative process of visual inspection, editing and re-running of Freesurfer 5.1 as required, following the recommended protocols (http://surfer.nmr/mgh.harvard.edu/fswiki/Edits). The rigorous quality control process allowed the achievement of accurate representations of the pial and white matter boundary.
Statistical analyses
IBM SPSS Statistics for macOS, Version 25.0 (IBM Corp. Released 2017, Armonk, NY) was used to conduct statistical analysis, including tests to examine the effects of the MRI scanner site, ICV, age and gender on the average left and right grey matter thickness and surface area of the HC participants.
The software applications mris_preproc, mri_surf2surf and mri_glmfit (Freesurfer 6.0) were used to perform group analyses and correlations at the vertex level with the cortical measures (grey matter thickness and surface area, respectively). The correlation analysis of grey matter thickness and surface area included symptom (CAARMS composite scores) and functional ratings (GAF, SOFAS), as well as AUDIT and CUDIT scores at the vertex level with both cortical measures.
For all surface area analyses, the ICV [43] was included as a nuisance variable. Freesurfer 6.0 (http://surfer.nmr.mgh.harvard.edu/) was used for these analyses because the Freesurfer 5.1 version of mris_preproc software does not apply a Jacobian correction for the surface area by default when transforming to the average space (target atlas, fsaverage). When using Freesurfer 6.0 , the total quantity of surface area for each subject is conserved across the transformation to the target atlas (fsaverge [44] .
A full-width half-maximum (FWHM) kernel of 20 mm was used for grey matter thickness and surface area analyses, together with a frontal-temporal mask to define the region of interest (Figure 2). The frontal-temporal mask was derived from merging the frontal and temporal regions as described by the Desikan-Killiany Atlas [45]. The merged parcellations from the frontal lobe included the superior frontal, rostral middle frontal, caudal middle frontal, pars opercularis, pars orbitalis, pars triangularis, lateral orbitofrontal, medial orbitofrontal, precentral, paracentral, and frontal pole. The temporal lobe regions included were the superior temporal, middle temporal, inferior temporal, banks of the superior temporal sulcus, fusiform, transverse temporal, entorhinal, temporal pole, and parahippocampal.
Insert Figure 2 about here
Multiple comparison correction was performed using permutation testing with 10,000 tests, a cluster-forming threshold of .05 and a cluster-wise threshold of .05. Permutation testing was used to control for false positives that may occur with the specified parameter settings. Monte Carlo simulations were used to perform the permutation testing [44]. In addition, a Bonferroni correction was applied to take both hemispheres into account. Finally, in order to optimise the computational time and to improve the accuracy of the p value when clusters were found with p < .1 after 10,000 tests, permutation testing with 100,000 tests were performed.