A 2-year longitudinal investigation of insula subregional volumes in early psychosis ==================================================================================== * Andrew R. Kittleson * Maureen McHugo * Jinyuan Liu * Simon N. Vandekar * Kristan Armstrong * Baxter Rogers * Neil D. Woodward * Stephan Heckers * Julia M. Sheffield ## Abstract **Background** The insula is a heterogeneous cortical region with three cytoarchitectural subregions— agranular, dysgranular, and granular—that have distinct functional roles. Previous cross- sectional studies have shown smaller volume of all insula subregions in individuals with psychotic disorders. However, longitudinal trajectories of insula subregions in early psychosis, and the relationship between subregional volumes and relevant clinical phenomena, such as perceptual aberrations, have not been previously examined. **Methods** 66 early psychosis (EP) and 65 healthy comparison (HC) participants completed 2-4 study visits over 2 years. T1-weighted structural brain images were processed using longitudinal voxel- based morphometry in CAT12 and segmented into anatomic subregions. At baseline, participants completed the Perceptual Aberrations Scale (PAS) to capture bodily distortions. The EP group was further examined based on diagnostic trajectory over two years (stable schizophrenia, stable schizophreniform, and conversion from schizophreniform to schizophrenia). **Results** EP participants had smaller insula volumes in all subregions compared to HC participants, and these volumes were stable over two years. Compared to HC, insula volumes were significantly smaller in EP participants with a stable diagnosis of schizophrenia, but other diagnostic trajectory groups did not significantly differ from HC or the stable schizophrenia group. While perceptual aberrations were significantly elevated in EP participants, PAS scores were not significantly related to insula volume. **Conclusions** We find that all insula subregions are smaller in early psychosis and do not significantly decline over two years. These data suggest that all insula subregions are structurally impacted in schizophrenia-spectrum disorders and may be the result of abnormal neurodevelopment. Keywords * Insula * schizophrenia * early psychosis * structural neuroimaging * perceptual aberrations * longitudinal * diagnostic trajectory * clinical subgroups ## Introduction The human insula is a heterogeneous cortical brain region that is implicated in a wide variety of behaviors. Sometimes referred to as a fifth cerebral lobe, this region is among the first cortical areas to mature and continues to experience rapid growth into the early postnatal period1,2. In schizophrenia, insula abnormalities are widely reported, including structural alterations3, structural connectivity differences4, resting-state functional connectivity changes5,6, and abnormal task-based functional activation signatures7,8. One of the most common findings in this literature is the presence of smaller insula volumes in schizophrenia, which are observed both in first episode and chronic stages of illness9. Insula volumes are also smaller, albeit to a lesser extent, in those at clinical high-risk and with enriched genetic risk for psychosis10–12, suggesting that structural changes precede clinical diagnosis. Yet, the insula itself is not a homogenous region, and can be parsed into more precise subregions based on its cytoarchitecture, anatomy, and associated realms of function13. Structural parcellations of the insula often focus on three subregions, distinguished by the relative presence or absence or granular layer IV cortical neurons: agranular (anterior), dysgranular (central), and granular (posterior)14,15. These differing neural substrates have implications for associated function16. Von Economo neurons (VENs), spatially localized to the agranular insula in cortical layer V, are specialized bipolar neurons that allow for rapid long- range information transfer between brain regions17,18. They are thought to facilitate the agranular insula’s roles in complex socio-emotional processing like prosocial behaviors and social empathy19–22. Conversely, layer IV granular neurons found densely in the posterior granular insula are smaller and more superficial than VENs. These neurons receive neocortical, subcortical, and even visceral inputs23–25, and in synthesizing this sensory information create perceptions of exteroceptive and interoceptive stimuli, implicating the posterior insula in self- embodiment26,27. Distinct patterns of subregional volume change over time in early psychosis could reflect selective neuronal degradation in the insula, providing one potential mechanism through which neurobiology affects behavior. Prior work examining insula volumes in psychotic disorders have largely neglected these subregional distinctions, yet growing evidence points to their relevance. Findings from a large neurodevelopmental sample found smaller insula volumes for the agranular and dysgranular, but not the granular, subregions in psychosis-spectrum youth compared to typically developing youth and youth with other psychopathologies28. In addition, individuals with schizophrenia- spectrum disorders, including both those in early and chronic stages, have lower insula volumes compared with nonclinical participants in all subregions. Critically, this volume deficit was not observed in psychotic bipolar disorder, suggesting an additional consideration of diagnostic specificity28. Interestingly, the only prior study on longitudinal insula volume in early psychosis reported increases in anterior insula volumes and stable posterior insula volumes in psychosis participants scanned several years after their first-episode of psychosis, using broad anterior/posterior subdivisions and a highly heterogeneous definition of first-episode psychosis29. Together, this evidence suggests that insula volume changes in early psychosis may be specific to schizophrenia-spectrum disorders and has some subregional specificity. Yet, how these cytoarchitectural insula subregions change within the early course of illness, in a well-characterized cohort of individuals with schizophrenia-spectrum disorders, remains unknown. Finally, even within a longitudinal cohort that is robustly characterized and closely followed, schizophrenia-spectrum diagnoses exist along a spectrum of clinical and functional impairment. While some first-episode patients may never meet diagnostic criteria for schizophrenia (e.g., receive and maintain a diagnosis of schizophreniform disorder), others may progress to meet this diagnosis due to persistent or worsening symptoms (e.g., receive a schizophreniform diagnosis during the first episode, but subsequently meet criteria for schizophrenia)30,31. With this context, increasing evidence continues to underscore the importance of the insula in the phenomenology of psychotic illness28,32,33. With its roles in cognition, emotion, interoception, and salience processing34,35, the insula is crucial for a variety of functions that can go awry in psychosis, progressively and throughout the course of illness36–38. Indeed, it is possible that individuals with more persistent clinical and functional impairment may demonstrate smaller subregional volumes that decline over the early course of illness. Conversely, if insula alterations are similar across different diagnostic trajectories, this might suggest that insula volume deficits represent a developmental insult that precedes the earliest manifestations of a schizophrenia-spectrum disorder and do not clearly contribute to illness course. The current study aimed to elucidate insula subregional volume changes in early psychosis. A sample of individuals identified within two years of their first episode for a primary psychotic disorder was followed for two years with structural neuroimaging occurring every 6-8 months. Within this sample, we tested: 1) whether insula subregions demonstrated volumetric differences in early psychosis, 2) how subregional volumes changed over the course of two years in early psychosis participants, and 3) whether diagnostic trajectory subgroups had different patterns of insula volume changes over time. Finally, we explored how insula volumes related to perceptual bodily aberrations, which are linked with insula function, especially in the granular subregion in schizophrenia39. ## Methods and Materials ### Study participants One hundred forty-six (146) total participants (72 healthy comparison [HC] participants and 74 participants in the early stage of a psychotic disorder [EP]) were recruited for a prospective longitudinal study over the early course of psychosis. To target structural pathology unique to the early course of psychotic disorders, EP participants were recruited if it had been less than two years since their first episode of psychosis, with the majority of EP participants recruited within 6-8 months of illness presentation (average duration of psychosis at baseline visit, as measured using Symptom Onset in Schizophrenia Inventory [SOS]40, was 6.2 months). EP participants were recruited from outpatient and inpatient settings at Vanderbilt Psychiatric Hospital in Nashville, TN, while HC participants were recruited from the surrounding community. Inclusion and exclusion criteria for both groups can be found in the Supplement. See Table 1 for participant demographics and cognitive characteristics across the entire sample. View this table: [Table 1.](http://medrxiv.org/content/early/2024/11/27/2024.11.25.24317916/T1) Table 1. Participant demographics and clinical characteristics at baseline. Abbreviations: SCIP = Screen of Cognition in Psychiatry; WTAR = Wechsler Test of Adult Reading; PANSS = Positive and Negative Syndrome Scale; BPD-PF = bipolar disorder with psychotic features; SZA = schizoaffective disorder; SZF = schizophreniform; SZ = schizophrenia; APD = antipsychotic drug. From the total sample recruited, 65 HC (90.3%) and 66 EP (89.2%) participants completed scans that passed quality control thresholds. The majority (75% of HC participants and 62.1% of EP participants) of these study participants completed four study visits every 6-8 months over two years, each consisting of a structural MRI scan. Every participant included in analyses completed at least two study visits (Table 1). All participants provided informed consent to be enrolled in the study and were monetarily compensated for each visit they completed. This study was approved by the Vanderbilt University Institutional Review Board. Further details about participant attrition are included in the Supplement (see Figure S1). ### Clinical and cognitive characterization At baseline and final (i.e., 2-year follow-up) study visits, participants completed an in-person interview including the Structured Clinical Interview for DSM-IV, TR (SCID)41. Information from available medical records was reviewed and provided a complement to SCID results; all clinical data was reviewed for diagnostic consensus with S.H., practicing psychiatrist. For EP participants, psychotic symptomatology was assessed using the Positive and Negative Syndrome Scale (PANSS)42 at each visit. Current antipsychotic treatment and antipsychotic treatment dosages were assessed at every visit and converted into standardized chlorpromazine equivalents43,44. Diagnostic trajectory subgroupings were operationalized based on an EP participant’s diagnoses at baseline and 2-year follow-up visits. Three trajectories of illness were defined based on whether the participant (1) maintained a diagnosis of schizophreniform disorder over the course of the study (StableSZF; N = 17), (2) entered the study with a diagnosis of schizophreniform disorder at baseline but met criteria for schizophrenia or schizoaffective disorder at 2-year follow-up visit (ProgToSZ; N = 31), or (3) initially met criteria for schizophrenia or schizoaffective disorder upon entering the study (StableSZ; N = 18). The Wechsler Test of Adult Reading (WTAR) was used to estimate premorbid IQ at baseline45. Current cognitive function at each visit was assessed using the Screen for Cognitive Impairment in Psychiatry (SCIP)46. At baseline, participants completed the Perceptual Aberrations Scale (PAS), a 35-item True/False validated self-report measure of distortions in perception of the body and other objects47. The PAS is composed of statements about psychotic-like bodily distortions and perceptual experiences (e.g., “I can remember when it seemed as though one of my limbs took on an unusual shape.”); changes in body boundaries, parts, or appearance; and unreality (e.g., “I have sometimes felt confused as to whether my body was really my own.”). More information on clinical and cognitive characteristics across participant groups can be found in Table 1; demographic information for diagnostic trajectory subgroups is located in the Supplement. ### Data acquisition and processing Structural brain imaging data were collected on one of two 3 Tesla scanners (Philips Healthcare Inc., Best, Netherlands) with 32-channel head coils at the Vanderbilt University Institute of Imaging Sciences in Nashville, TN. At each scanning session, T1-weighted images were collected (voxel size = 1 mm3, field of view = 256 mm2, number of slices = 170, gap = 0 mm, echo time = 3.7 ms, repetition time = 8.0 ms). All images were visually inspected for artifacts prior to inclusion in analyses; no images were removed. In sum, 458 total images were collected (223 EP, 235 HC)—131 at scan 1, 131 at scan 2, 110 at scan 3, and 86 at scan 4. T1-weighted images were processed using the CAT12 longitudinal toolbox (Computational Anatomy Toolbox, Structural Brain Mapping group, Jena, Germany)48 for voxel-based morphometric analyses of normalized, modulated gray matter volumes. This pipeline is designed specifically for the analysis of longitudinal imaging data and includes spatial registration, tissue classification, intensity bias correction, and scaling. To test our hypotheses that insula subregional volumes change during the early course of psychosis, we used insular subregional masks developed by Farb and colleagues49 (see Figure S1) to estimate gray matter volume within the insula. These masks were based on the manual segmentation of insular cortex into agranular, dysgranular, and granular subregions, based on cytoarchitectural characteristics. Briefly, insula subregions of interest (ROIs) were manually drawn by Farb and colleagues to fit each visible gyrus of the template brain in MRIcron, yielding eight anatomically defined regions: five dorsal (four dysgranular and one granular) and three ventral (one agranular, one dysgranular, and one granular). ROIs were collapsed across dorsal and ventral divisions to create three nonoverlapping anatomic masks for insula subregions: agranular, dysgranular, and granular. A whole insula mask was produced by combining these three masks in MRIcron50. Subregional volumes were averaged across both hemispheres to reduce parameters in subsequent statistical modeling. ### Statistical analysis Statistical analyses were carried out using linear mixed effects models in R (R Core Team, 2019) with the lme4, rstatix, emmeans, car, and tidyverse packages51–55. A power analysis conducted prior to data analysis determined that with the present sample size, there would be 80% power to detect an effect size of *d*lJ=lJ0.46 in the primary longitudinal analysis with a type 1 error rate of 0.05. Linear mixed effects models testing group differences between HC and EP participants included age at baseline visit, sex assigned at birth, total intracranial volume (TIV), and scanner as fixed effects and included participant as a random effect. Models comparing longitudinal volume differences among diagnostic trajectory subgroups (HC, StableSZF, ProgToSZ, and StableSZ) were constructed using the same fixed and random effects. Cohen’s *d* effect sizes for follow-up contrasts were calculated from the F-statistics using the effectsize package in R56. First, to investigate how group influenced insula volumes, we constructed four models with insula Volume (whole, agranular, dysgranular, and granular) as the dependent variables and Group (EP, HC) as a fixed effect. Next, to determine trajectories of insula volume deficits over time in early psychosis, we constructed four separate models with insula Volume (whole, agranular, dysgranular, and granular) as the dependent variable, Group (EP, HC) and Time (months since baseline scan) as fixed effects, and a Group x Time interaction term. Models were tested with and without interaction terms; main effects reported are from models without interactions. We also tested our hypothesis that subregional insula volumes changed in group-specific ways over the early course of psychosis by constructing one model with Volume as the dependent variable, and Group (EP, HC), Subregion (agranular, dysgranular, and granular), and Time (months since baseline scan) as fixed effects. We then constructed similar models as above, but with Group defined by four diagnostic trajectory subgroups (HC; StableSZF; ProgToSZ; and StableSZ). Finally, to examine potential brain-behavior relationships between perceptual aberrations and insula volumes, we looked at PAS scores between Group (EP, HC) and tested associations between PAS score and both whole insula volumes and granular insula volumes. We tested these associations across all participants and then evaluated the Group difference. One EP outlier who scored 4+ standard deviations above the mean on the PAS was removed for these analyses. Relationships between insula subregional volumes and clinical and cognitive characteristics in the EP sample were also explored and are provided in the Supplement (see Table S2). ## Results ### Overall differences in insula volumes Volume was significantly lower in EP participants than HC participants for the whole insula (main effect of Group: F1,166.97 = 18.75, p < 0.001, *d* = 0.67; Figure 1) and within all three subregions: agranular (main effect of Group: F1,128.56 = 12.02, p < 0.001, *d* = 0.61), dysgranular (main effect of Group: F1,165.41 = 16.06, p < 0.001, *d* = 0.62) and granular (main effect of Group: F1,164.62 = 14.07, p < 0.001, *d* = 0.58; Figure 1). However, the magnitude of these group differences did not significantly differ across subregions (Group x Subregion interaction: F3,1691.84 = 0.49, p = 0.69). ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/11/27/2024.11.25.24317916/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2024/11/27/2024.11.25.24317916/F1) Figure 1. Overall whole, agranular, dysgranular, and granular insula volumes are all lower in EP participants than in HC participants. Error bars represent standard error above and below group means. Abbreviations: HC = healthy comparison; EP = early psychosis. Asterisk (*****) = p < 0.05. ### Longitudinal differences in insula volumes In longitudinal models of the whole insula, we found evidence of significantly decreasing insula volumes over two years (main effect of Time: F1,355.71 = 11.47, p < 0.001, *d* = 0.36), which differed between groups (Group x Time interaction: F1,327.87 = 6.23, p = 0.01, *d* = 0.28; Figure S3, Supplement). This was driven by HC insula volumes decreasing significantly over time (main effect of Time: whole insula: F1,179.19 = 23.46, p < 0.001, *d* = 0.72), which was not observed in EP participants (main effect of Time: whole insula: F1,171.32 = 0.11, p = 0.74, *d* = 0.05). Across all participants, we observed volume decline in all insula subregions: agranular (main effect of Time: F1,349.50 = 16.57, p < 0.001), dysgranular (main effect of Time: F1,354.66 = 3.99, p = 0.046), and granular (main effect of Time: F1,354.12 = 13.80, p < 0.001). However, there was no significant Time x Subregion interaction (F2,1235.58 = 0.11, p = 0.90). A significant Group x Time interaction was found for the dysgranular and granular subregions (dysgranular: F1,326.96 = 11.29, p < 0.001; granular: F1,326.41 = 4.02, p = 0.046), again driven by significantly decreasing subregional volumes in HC participants but not in EP participants (see Figure S3). However, the three-way interaction (Group x Subregion x Time) was not significant (F3,1684.96 = 0.58, p = 0.63), suggesting that the difference in volumetric trajectories between groups was not significantly different across the three subregions. ### Overall differences in insula volumes: effect of diagnostic trajectory We next examined how insula volumes differed amongst four diagnostic trajectory subgroups—1) HC (N = 65), 2) StableSZF (N = 17), 3) ProgToSZ (N = 31), and 4) StableSZ (N = 18). Whole insula and all subregions significantly differed across these four groups (whole insula main effect of Trajectory: F3,124.22 = 5.76, p = 0.001; agranular main effect of Trajectory: F3,124.82 = 4.75, p = 0.004; dysgranular main effect of Trajectory: F3,123.39 = 4.55, p = 0.005; granular main effect of Trajectory: F3,122.61 = 4.18, p = 0.007). We did not find subregional specificity in insula volume deficits that varied with diagnostic trajectory (Trajectory x Subregion interaction: F6,1231.44 = 1.12, p = 0.35). Post-hoc, Bonferroni-corrected pairwise comparisons were conducted to examine the pattern of volume differences across trajectory groups (Table 2). HC participants had significantly greater insula volumes across all measures compared only to the stable schizophrenia subgroup (main effect of Trajectory for HC vs. StableSZ: whole insula: F1,82.24 = 12.46; pcorrected < 0.001, *d* = 0.78; agranular insula: F1,81.54 = 10.32; pcorrected = 0.01, *d* = 0.71; dysgranular insula: F1,81.39 = 9.82; pcorrected = 0.01, *d* = 0.69; granular insula: F1,80.97 = 9.18; pcorrected = 0.02, *d* = 0.67; Figure 2). None of the diagnostic trajectory groups differed from one another. Of note, duration of illness did not significantly correlate with volume for any of the subregions (see Table S2). View this table: [Table 2.](http://medrxiv.org/content/early/2024/11/27/2024.11.25.24317916/T2) Table 2. Comparisons of insula volumes (whole, agranular, dysgranular, and granular) for each diagnostic trajectory subgroup pairing. Significance was determined at p < 0.008 for six comparisons; p values shown below are Bonferroni-uncorrected. Significant comparisons that survived Bonferroni correction for multiple comparisons are **bolded**. Abbreviations: SZF = schizophreniform; SZ = schizophrenia. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/11/27/2024.11.25.24317916/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2024/11/27/2024.11.25.24317916/F2) Figure 2. Whole insula volumes differed significantly among diagnostic trajectory subgroups (**A**), as did agranular, dysgranular, and granular subregional insula volumes (**B**). Significant pairwise differences (Bonferroni-corrected p < 0.008 for six comparisons) are indicated for each volume measure by *****; error bars represent standard error above and below subgroup estimated marginal means. Abbreviations: SZF = schizophreniform; SZ = schizophrenia. Asterisk (*****) = p < 0.05. ### Longitudinal differences in insula volumes: effect of diagnostic trajectory We did not find evidence of differences in whole insula volume change over time, based on diagnostic trajectory (Trajectory x Time interaction: F3,323.91 = 2.38, p = 0.07; Figure 3), nor did we find any subregional specificity in insula volume deficits that varied with diagnostic trajectory over time (Trajectory x Subregion x Time interaction: F6,1220.54 = 0.42, p = 0.86). Follow-up analyses revealed no evidence of significant changes in any insula volume (whole, agranular, dysgranular, or granular) over time in any of the three clinical trajectory subgroups (Table 3). ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/11/27/2024.11.25.24317916/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2024/11/27/2024.11.25.24317916/F3) Figure 3. Whole (**A**) insula volumes and some (dysgranular and granular) but not all subregional (**B**) insula volumes change differentially over time among diagnostic trajectory subgroups. Volumes shown are estimated marginal means at timepoints 0 months, 8 months, 16 months, and 24 months; error bars represent standard error above and below the estimated marginal mean. Abbreviations: SZF = schizophreniform disorder; SZ = schizophrenia. View this table: [Table 3.](http://medrxiv.org/content/early/2024/11/27/2024.11.25.24317916/T3) Table 3. Main effect of time for each diagnostic trajectory subgroup for each insula volume measure. Bonferroni-corrected p < 0.013 are considered significant (four comparisons per volume measure). Significant comparisons that survived Bonferroni correction for multiple comparisons are **bolded**. Abbreviations: SZF = schizophreniform; SZ = schizophrenia. ### Baseline associations of insula volumes with perceptual aberrations PAS scores were significantly elevated in EP compared to HC participants at baseline (*H* = 23.32, p < 0.001; Figure 4A). Across all participants, we did not find evidence for a significant association between PAS total and baseline whole insula volume (rS = -0.06, p = 0.50; Figure 4B) or baseline granular volume (rS = -0.05, p = 0.58; Figure 4C). Finally, based on our a priori interest in the granular insula, we did not find evidence of baseline PAS total scores significantly predicting granular subregional volume changes across EP participants in a linear mixed model (main effect of PAS: F1,58.69 = 0.026, p = 0.87). ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/11/27/2024.11.25.24317916/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2024/11/27/2024.11.25.24317916/F4) Figure 4. Self-reported perceptual aberrations (as measured by PAS total) are higher in EP participants at baseline (**A**). Neither whole insula volumes (**B**) nor granular insula subregional volumes (**C**) were significantly associated with self-reported perceptual aberrations at baseline among all participants. Abbreviations: EP = early psychosis; HC = healthy comparison; PAS = Perceptual Aberrations Scale. Asterisk (*****) = p < 0.05. ## Discussion With a longitudinal case-control study design we show that insula volumes are smaller by first episode of psychosis and do not change over two years. This baseline deficit is present in all granular subregions without apparent specificity, although individuals who already meet criteria for schizophrenia at the time of study entry show evidence of the largest deficits compared to nonclinical samples. Finally, despite higher self-reported perceptual aberrations in EP participants, no relationships were observed with insula volumes. Overall, our data demonstrate that individuals with schizophrenia-spectrum disorders have smaller insula volumes across all granular subregions at first episode of psychosis and that trajectory of illness may be important for understanding this process. Our results suggest that structural abnormalities of the insula occur prior to illness onset as neurodevelopmental insults. The development of the insula as a cortical brain region consists of several overlapping but temporally distinct processes, beginning with gyrification in prenatal development57,58, followed by sulcal deepening soon after birth59,60, cortical thickening in early childhood61,62, and subsequent cortical thinning in adolescence into early adulthood63. Prior work has shown perturbations to the insula in all of these structural metrics in individuals with psychosis64–68. Excessive synaptic pruning during adolescence has long been proposed as one mechanism by which individuals with schizophrenia develop disproportionally high gray matter loss69, and neuroimaging studies of adolescents with child-onset psychosis have shown evidence of earlier and more severe cortical gray matter loss compared to non-clinical samples70. Importantly, there seems to be some diagnostic specificity in these perturbations, as similar structural alterations to the insula have not been observed in individuals with other affective psychotic disorders such as psychotic bipolar disorder28. These morphological abnormalities in the insula, seemingly unique to schizophrenia, may underlie the earliest phenomenological signs of illness. Furthermore, it is important to note that our individuals with psychosis did not demonstrate gradually declining volumes over time as HC participants did. This may be due to prior gray matter loss in the prodrome. Indeed, across individuals at clinical high-risk for psychosis, cortical gray matter volume declines steepest for those who will convert to psychosis71, including within the insula10. In contrast, in studies of individuals with established schizophrenia-spectrum disorders, insula volumes did not decline progressively over time72. Taken together, this literature on insula structure in psychosis suggests that volume changes predate diagnosis onset, occur at least in part during the prodrome, and subsequently lessen in more chronic stages. The results of this study situate well into this existing literature by confirming that insula structural changes in psychosis occur prior to the onset of illness and may be pronounced for individuals with more severe illness trajectories, reflected by the smallest insula volumes among individuals who already met criteria for schizophrenia at the time of the study. In addition to the neurobiological mechanism underlying these structural brain changes in psychosis, it is crucial to better understand how the resulting aberrant volumes may affect clinical phenotypes, and how the insula may impact psychotic phenomenology more broadly. With its varied roles in somatosensation, salience detection, and social cognition, the insula is a key brain region for understanding the phenomenology of reality distortions, perceptual aberrations, and social cognitive impairments that characterize psychosis39,73,74. Prior literature has proposed that insula subregions are important for specific realms of function in normative populations and dysfunction in schizophrenia16,32. Thus, the extent of volume loss in specific insula subregions may also provide a neural basis for worsening psychotic symptomatology in specific domains, which has been shown across various realms of function and subregions of the insula75. Our evidence of comparable volume loss across all three subregions of the insula in early psychosis suggests that all associated symptom realms may be affected. With this context, we explored how insula gray matter volume might be associated with self- reported perceptual aberrations to investigate other brain-behavior relationships relevant to the insula. The insula is a key node within a proposed “interoceptive network” in schizophrenia76, so we investigated how volume changes might affect interoceptive deficits to extend our imaging findings into clinically meaningful brain-behavior relationships. A previous study by Brosey and Woodward of perceptual aberrations and insula volumes found a significant negative correlation between PAS scores and bilateral whole and granular insula volumes in participants with chronic psychosis39, which was not observed in our cohort. In fact, while our EP participants did score higher on the PAS than HC participants and recapitulate prior behavioral findings77, we did not find evidence of any significant relationships between self-reported perceptual aberrations at baseline and insula volumes. However, it may be that this association is not apparent until symptoms progress further; our cohort was captured soon after illness onset, while Brosey and Woodward’s cohort was older and more chronic in their illness course. Additionally, the phenomenology of psychosis is highly complex and multifaceted, requiring many brain regions and circuits to coordinate efforts78. Indeed, perceptual aberrations represent one such behavioral readout of psychotic symptomatology, and multiple brain regions, including but not limited to the insula, are involved in their production39. Probing the functional activation and connectivity of the insula, instead of investigating structural changes as here, may in fact be better for understanding links between this region and symptoms of schizophrenia-spectrum disorders. In sum, the current study is the first longitudinal investigation of both insula subregional volume and diagnostic trajectories in the early period of psychosis, and provides a valuable view into brain changes during this critical window. Our study suggests that insula volume changes in early psychosis differ from typical neurodevelopmental trajectories, and that illness-related volume changes likely precede clinical presentation and diagnosis. We also provide evidence here that among those with an emerging psychotic disorder, insula volumes are most different among patients with chronic schizophrenia compared to nonclinical comparisons. This study provides insight into the underlying neurobiological changes in the insula that may precede or occur alongside symptom progression in psychotic illness. Altered insula volumes, especially at the subregional level, may be an important developmental factor among those at high risk for developing schizophrenia. ## Supporting information Supplement [[supplements/317916_file06.docx]](pending:yes) ## Data Availability All data produced in the present study are available upon reasonable request to the authors. ## Disclosures The authors declare that they have no conflicts of interest or financial disclosures to report. ## Acknowledgments The work outlined in this manuscript was supported by the Charlotte and Donald Test Fund, NIMH grants R01-MH70560 (Heckers) R01-MH123563 (Vandekar), the Jack Martin, MD Research in Psychopharmacology endowment (Sheffield) and K23-MH126313. This work was conducted in part using the resources of the Center for Computational Imaging at the Vanderbilt University Institute of Imaging Sciences and the Advanced Computing Center for Research and Education (ACCRE) at Vanderbilt University in Nashville, TN. The authors would like to thank the participants of this study for their involvement in and commitment to research. * Received November 25, 2024. * Revision received November 25, 2024. * Accepted November 27, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Nishikuni, K. & Ribas, G. C. Study of fetal and postnatal morphological development of the brain sulci: Laboratory investigation. J. Neurosurg. Pediatr. 11, 1–11 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3171/2012.9.PEDS12122&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23140215&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 2. 2.Wai, M. S. M. et al. Development of the human insular cortex: differentiation, proliferation, cell death, and appearance of 5HT-2A receptors. Histochem. Cell Biol. 130, 1199–1204 (2008). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18762964&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 3. 3.Goodkind, M. et al. Identification of a Common Neurobiological Substrate for Mental Illness. JAMA Psychiatry 72, 305–315 (2015). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25651064&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 4. 4.Ellison-Wright, I. & Bullmore, E. Meta-analysis of diffusion tensor imaging studies in schizophrenia. Schizophr. Res. 108, 3–10 (2009). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.schres.2008.11.021&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19128945&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000264506500002&link_type=ISI) 5. 5.Sheffield, J. M., Rogers, B. P., Blackford, J. U., Heckers, S. & Woodward, N. D. Insula Functional Connectivity in Schizophrenia. Schizophr. Res. 220, 69–77 (2020). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32307263&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 6. 6.Tian, Y., Zalesky, A., Bousman, C., Everall, I. & Pantelis, C. Insula Functional Connectivity in Schizophrenia: Subregions, Gradients, and Symptoms. Biol. Psychiatry Cogn. Neurosci. Neuroimaging 4, 399–408 (2019). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30691966&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 7. 7.Moran, L. V. et al. Disruption of Anterior Insula Modulation of Large-Scale Brain Networks in Schizophrenia. Biol. Psychiatry 74, 467–474 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2013.02.029&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23623456&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000324051200013&link_type=ISI) 8. 8.Palaniyappan, L., Simmonite, M., White, T. P., Liddle, E. B. & Liddle, P. F. Neural Primacy of the Salience Processing System in Schizophrenia. Neuron 79, 814–828 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2013.06.027&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23972602&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000323587000018&link_type=ISI) 9. 9.Shepherd, A. M., Matheson, S. L., Laurens, K. R., Carr, V. J. & Green, M. J. Systematic Meta-Analysis of Insula Volume in Schizophrenia. Biol. Psychiatry 72, 775–784 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2012.04.020&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22621997&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000309610500010&link_type=ISI) 10. 10.Takahashi, T. et al. Insular cortex gray matter changes in individuals at ultra-high-risk of developing psychosis. Schizophr. Res. 111, 94–102 (2009). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.schres.2009.03.024&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19349150&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000266745900013&link_type=ISI) 11. 11.Zhu, X. et al. Polygenic Risk for Schizophrenia, Brain Structure, and Environmental Risk in UK Biobank. Schizophr. Bull. Open 2, sgab042 (2021). 12. 12.Fortea, A. et al. Cortical gray matter reduction precedes transition to psychosis in individuals at clinical high-risk for psychosis: A voxel-based meta-analysis. Schizophr. Res. 232, 98–106 (2021). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34029948&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 13. 13.Uddin, L. Q., Kinnison, J., Pessoa, L. & Anderson, M. L. Beyond the Tripartite Cognition–Emotion–Interoception Model of the Human Insular Cortex. J. Cogn. Neurosci. 26, 16–27 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access\_num=10.1162/jocn_a_00462&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23937691&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000327430900002&link_type=ISI) 14. 14.Quabs, J. et al. Cytoarchitecture, probability maps and segregation of the human insula. NeuroImage 260, 119453 (2022). 15. 15.Paxinos, G. & Watson, C. The Rat Brain in Stereotaxic Coordinates: Hard Cover Edition. (Elsevier, 2006). 16. 16.Kittleson, A. R., Woodward, N. D., Heckers, S. & Sheffield, J. M. The insula: Leveraging cellular and systems-level research to better understand its roles in health and schizophrenia. Neurosci. Biobehav. Rev. 160, 105643 (2024). 17. 17.Allman, J. M. et al. The von Economo neurons in fronto-insular and anterior cingulate cortex. Ann. N. Y. Acad. Sci. 1225, 59–71 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1749-6632.2011.06011.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21534993&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000291368900007&link_type=ISI) 18. 18.Allman, J. M. et al. The von Economo neurons in frontoinsular and anterior cingulate cortex in great apes and humans. Brain Struct. Funct. 214, 495–517 (2010). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00429-010-0254-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20512377&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000278898300009&link_type=ISI) 19. 19.Baez, S., García, A. M. & Ibanez, A. The Social Context Network Model in Psychiatric and Neurological Diseases. Curr. Top. Behav. Neurosci. 30, 379–396 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access\_num=10.1007/7854_2016_443&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27130326&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 20. 20.Masten, C. L., Morelli, S. A. & Eisenberger, N. I. An fMRI investigation of empathy for ‘social pain’ and subsequent prosocial behavior. NeuroImage 55, 381–388 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2010.11.060&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21122817&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000287008900038&link_type=ISI) 21. 21.Baez, S. & Ibanez, A. The effects of context processing on social cognition impairments in adults with Asperger’s syndrome. Front. Neurosci. 8, (2014). 22. 22.Singer, T., Critchley, H. D. & Preuschoff, K. A common role of insula in feelings, empathy and uncertainty. Trends Cogn. Sci. 13, 334–340 (2009). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tics.2009.05.001&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19643659&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000269411800005&link_type=ISI) 23. 23.Kurth, F. et al. Cytoarchitecture and Probabilistic Maps of the Human Posterior Insular Cortex. Cereb. Cortex N. Y. NY 20, 1448–1461 (2010). 24. 24.Harris, K. D. & Shepherd, G. M. G. The neocortical circuit: themes and variations. Nat. Neurosci. 18, 170–181 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nn.3917&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25622573&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 25. 25.Zhang, Z. & Oppenheimer, S. M. Characterization, distribution and lateralization of baroreceptor-related neurons in the rat insular cortex. Brain Res. 760, 243–250 (1997). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0006-8993(97)00284-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9237541&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997XK47700031&link_type=ISI) 26. 26.Craig, A. D. How do you feel? Interoception: the sense of the physiological condition of the body. Nat. Rev. Neurosci. 3, 655–666 (2002). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn894&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12154366&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000177447800016&link_type=ISI) 27. 27.Seth, A. K. Interoceptive inference, emotion, and the embodied self. Trends Cogn. Sci. 17, 565–573 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tics.2013.09.007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24126130&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000327230500008&link_type=ISI) 28. 28.Sheffield, J. M. et al. Insula sub-regions across the psychosis spectrum: morphology and clinical correlates. Transl. Psychiatry 11, 1–13 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41398-020-01158-w&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33414379&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 29. 29.Pigoni, A. et al. Insula volumes in first-episode and chronic psychosis: A longitudinal MRI study. Schizophr. Res. 241, 14–23 (2022). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35074528&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 30. 30.Polari, A. et al. Clinical trajectories in the ultra-high risk for psychosis population. Schizophr. Res. 197, 550–556 (2018). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29463457&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 31. 31.Allswede, D. M. et al. Characterizing Covariant Trajectories of Individuals at Clinical High Risk for Psychosis Across Symptomatic and Functional Domains. Am. J. Psychiatry 177, 164–171 (2020). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31509005&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 32. 32.Wylie, K. P. & Tregellas, J. R. The Role of the Insula in Schizophrenia. Schizophr. Res. 123, 93–104 (2010). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.schres.2010.08.027&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20832997&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 33. 33.Palaniyappan, L. & Liddle, P. F. Does the salience network play a cardinal role in psychosis? An emerging hypothesis of insular dysfunction. J. Psychiatry Neurosci. JPN 37, 17–27 (2012). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21693094&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 34. 34.Menon, V. & Uddin, L. Q. Saliency, switching, attention and control: a network model of insula function. Brain Struct. Funct. 214, 655–667 (2010). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00429-010-0262-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20512370&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000278898300021&link_type=ISI) 35. 35.Gasquoine, P. G. Contributions of the Insula to Cognition and Emotion. Neuropsychol. Rev. 24, 77–87 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11065-014-9246-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24442602&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 36. 36.Bora, E., Yücel, M. & Pantelis, C. Cognitive Impairment in Affective Psychoses: A Meta-analysis. Schizophr. Bull. 36, 112–125 (2010). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19767349&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 37. 37.Henry, J. D. et al. Emotion dysregulation in schizophrenia: Reduced amplification of emotional expression is associated with emotional blunting. Schizophr. Res. 95, 197–204 (2007). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.schres.2007.06.002&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17630254&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249836900026&link_type=ISI) 38. 38.Torregrossa, L. J., Amedy, A., Roig, J., Prada, A. & Park, S. Interoceptive functioning in schizophrenia and schizotypy. Schizophr. Res. 239, 151–159 (2022). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34891079&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 39. 39.Brosey, E. A. & Woodward, N. D. Neuroanatomical correlates of perceptual aberrations in psychosis. Schizophr. Res. 179, 125–131 (2017). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27729190&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 40. 40.Perkins, D. O. et al. Characterizing and dating the onset of symptoms in psychotic illness: the Symptom Onset in Schizophrenia (SOS) inventory. Schizophr. Res. 44, 1–10 (2000). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0920-9964(99)00161-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10867307&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000088092000001&link_type=ISI) 41. 41.Glasofer, D., Brown, A. & Riegel, M. Structured Clinical Interview for DSM-IV (SCID). in 1–4 (2015). doi:10.1007/978-981-287-087-2_80-1. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/978-981-287-087-2_80-1&link_type=DOI) 42. 42.Kay, S. R., Fiszbein, A. & Opler, L. A. The Positive and Negative Syndrome Scale (PANSS) for Schizophrenia. Schizophr. Bull. 13, 261–276 (1987). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/schbul/13.2.261&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3616518&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1987H977500007&link_type=ISI) 43. 43.Leucht, S. et al. Dose Equivalents for Second-Generation Antipsychotic Drugs: The Classical Mean Dose Method. Schizophr. Bull. 41, 1397–1402 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/schbul/sbv037&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25841041&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 44. 44.Gardner, D. M., Murphy, A. L., O’Donnell, H., Centorrino, F. & Baldessarini, R. J. International Consensus Study of Antipsychotic Dosing. Focus 12, 235–243 (2014). 45. 45.Ginsberg, J. P. Review of Wechsler Test of Adult Reading. Appl. Neuropsychol. 10, 182–184 (2003). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1207/S15324826AN1003_08&link_type=DOI) 46. 46.Purdon, S. *Purdon* (2005) *SCIP Manual*. (2005). 47. 47.Chapman, L. J., Chapman, J. P. & Raulin, M. L. Perceptual Aberration Scale. American Psychological Association doi:10.1037/t20156-000 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1037/t20156-000&link_type=DOI) 48. 48.Gaser, C. et al. CAT: a computational anatomy toolbox for the analysis of structural MRI data. GigaScience 13, giae049 (2024). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/gigascience/giae049&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=39102518&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 49. 49.Farb, N. A. S., Segal, Z. V. & Anderson, A. K. Attentional Modulation of Primary Interoceptive and Exteroceptive Cortices. Cereb. Cortex N. Y. NY 23, 114–126 (2013). 50. 50.NITRC: MRIcron: Tool/Resource Info. [https://www.nitrc.org/projects/mricron](https://www.nitrc.org/projects/mricron). 51. 51.Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting Linear Mixed-Effects Models using lme4. Preprint at doi:10.48550/arXiv.1406.5823 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.48550/arXiv.1406.5823&link_type=DOI) 52. 52.Lenth, R. V., et al. emmeans: Estimated Marginal Means, aka Least-Squares Means. (2024). 53. 53.Kassambara, A. rstatix: Pipe-Friendly Framework for Basic Statistical Tests. (2023). 54. 54.Fox, J. & Weisberg, S. An R Companion to Applied Regression. (SAGE Publications, 2018). 55. 55.car package - RDocumentation. [https://www.rdocumentation.org/packages/car/versions/3.1-2](https://www.rdocumentation.org/packages/car/versions/3.1-2). 56. 56.Package index. [https://easystats.github.io/effectsize/reference/index.html](https://easystats.github.io/effectsize/reference/index.html). 57. 57.Armstrong, E., Schleicher, A., Omran, H., Curtis, M. & Zilles, K. The Ontogeny of Human Gyrification. Cereb. Cortex 5, 56–63 (1995). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cercor/5.1.56&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7719130&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995QD03700005&link_type=ISI) 58. 58.Zilles, K., Palomero-Gallagher, N. & Amunts, K. Development of cortical folding during evolution and ontogeny. Trends Neurosci. 36, 275–284 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tins.2013.01.006&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23415112&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000319719900003&link_type=ISI) 59. 59.Meng, Y., Li, G., Lin, W., Gilmore, J. H. & Shen, D. Spatial distribution and longitudinal development of deep cortical sulcal landmarks in infants. NeuroImage 100, 206–218 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2014.06.004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24945660&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000344235800018&link_type=ISI) 60. 60.Im, K. & Grant, P. E. Sulcal pits and patterns in developing human brains. NeuroImage 185, 881–890 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2018.03.057&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29601953&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 61. 61.Ducharme, S. et al. Trajectories of cortical thickness maturation in normal brain development — The importance of quality control procedures. NeuroImage 125, 267–279 (2016). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26463175&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 62. 62.Remer, J. et al. Quantifying cortical development in typically developing toddlers and young children, 1–6 years of age. NeuroImage 153, 246–261 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2017.04.010&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28392489&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 63. 63.Huttenlocher, P. R. Morphometric study of human cerebral cortex development. Neuropsychologia 28, 517–527 (1990). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0028-3932(90)90031-I&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2203993&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990DT71500002&link_type=ISI) 64. 64.Rimol, L. M. et al. Cortical Volume, Surface Area, and Thickness in Schizophrenia and Bipolar Disorder. Biol. Psychiatry 71, 552–560 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2011.11.026&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22281121&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000300789100016&link_type=ISI) 65. 65.Nesvåg, R. et al. Regional thinning of the cerebral cortex in schizophrenia: Effects of diagnosis, age and antipsychotic medication. Schizophr. Res. 98, 16–28 (2008). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.schres.2007.09.015&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17933495&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000252494100003&link_type=ISI) 66. 66.Spalthoff, R., Gaser, C. & Nenadić, I. Altered gyrification in schizophrenia and its relation to other morphometric markers. Schizophr. Res. 202, 195–202 (2018). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30049600&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 67. 67.Yan, J. et al. Cortical thinning and flattening in schizophrenia and their unaffected parents. Neuropsychiatr. Dis. Treat. 15, 935–946 (2019). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31114205&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 68. 68.Takahashi, T. et al. Gross anatomical variations of the insular cortex in first-episode schizophrenia. Schizophr. Res. 260, 23–29 (2023). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=37549494&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 69. 69.Feinberg, I. Schizophrenia: Caused by a fault in programmed synaptic elimination during adolescence? J. Psychiatr. Res. 17, 319–334 (1982). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0022-3956(82)90038-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7187776&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 70. 70.Gogtay, N. Cortical Brain Development in Schizophrenia: Insights From Neuroimaging Studies in Childhood-Onset Schizophrenia. Schizophr. Bull. 34, 30– 36 (2008). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/schbul/sbm103&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17906336&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000252148800008&link_type=ISI) 71. 71.Chung, Y. & Cannon, T. D. Brain Imaging During the Transition from Psychosis Prodrome to Schizophrenia. J. Nerv. Ment. Dis. 203, 336 (2015). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25900551&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 72. 72.Takahashi, T. et al. Follow-up MRI study of the insular cortex in first-episode psychosis and chronic schizophrenia. Schizophr. Res. 108, 49–56 (2009). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.schres.2008.12.029&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19171466&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000264506500007&link_type=ISI) 73. 73.Huang, M.-X. et al. Somatosensory System Deficits in Schizophrenia Revealed by MEG during a Median-Nerve Oddball Task. Brain Topogr. 23, 82–104 (2010). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19943100&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 74. 74.Green, M. F., Horan, W. P. & Lee, J. Social cognition in schizophrenia. Nat. Rev. Neurosci. 16, 620–631 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn4005&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26373471&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) 75. 75.Sheffield, J. et al. Volume, Morphology, and Clinical Correlates of Insula Sub- Regions in Psychotic Disorders and Psychosis-Spectrum Youth. Biol. Psychiatry 89, S280 (2021). 76. 76.Yao, B. et al. Structural connectivity of an interoception network in schizophrenia. Psychiatry Res. Neuroimaging 331, 111636 (2023). 77. 77.Coleman, M. J., Levy, D. L., Lenzenweger, M. F. & Holzman, P. S. Thought disorder, perceptual aberrations, and schizotypy. J. Abnorm. Psychol. 105, 469–473 (1996). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1037/0021-843X.105.3.469&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8772019&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996VA92200019&link_type=ISI) 78. 78.Heckers, S. Neural models of schizophrenia. Dialogues Clin. Neurosci. 2, 267 (2000). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22033839&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F11%2F27%2F2024.11.25.24317916.atom)