Learn more: PMC Disclaimer | PMC Copyright Notice
This is a preprint.
Sex differences in the functional network underpinnings of psychotic-like experiences in children
Associated Data
Abstract
Psychotic-like experiences (PLEs) include a range of sub-threshold symptoms that resemble aspects of psychosis but do not necessarily indicate the presence of psychiatric illness. These experiences are highly prevalent in youth and are associated with developmental disruptions across social, academic, and emotional domains. While not all youth who report PLEs develop psychosis, many develop other psychiatric illnesses during adolescence and adulthood. As such, PLEs are theorized to represent early markers of poor mental health. Here, we characterized the similarities and differences in the neurobiological underpinnings of childhood PLEs across the sexes using a large sample from the ABCD Study (n=5,260), revealing sex-specific associations between functional networks connectivity and PLEs. We find that although the networks associated with PLEs overlap to some extent across the sexes, there are also crucial differences. In females, PLEs are associated with dispersed cortical and non-cortical connections, whereas in males, they are primarily associated with functional connections within limbic, temporal parietal, somato/motor, and visual networks. These results suggest that early transdiagnostic markers of psychopathology may be distinct across the sexes, further emphasizing the need to consider sex in psychiatric research as well as clinical practice.
Introduction
Psychiatric illnesses are often preceded by premorbid and/or prodromal states that emerge months to years before the onset of a diagnosable illness[1]. These states include subtle changes in mood, behavior, thought, and cognition that do not meet the threshold for a formal diagnosis. Psychotic disorders manifest across a spectrum ranging from psychotic-like experiences (PLEs) to clinical high-risk states, first episode psychosis, and chronic psychotic disorder[2]. PLEs represent the earliest potential manifestations of psychosis and tend to emerge during childhood and adolescence. PLEs include a range of sub-clinical symptoms that resemble aspects of psychosis, such as perceptual distortions (i.e., hallucinations), false beliefs (i.e., delusions), and disorganized thinking or speech. As many as 60% of youth report PLEs[2, 3], and these experiences are associated with distress and can hinder social, academic, and emotional development[4]. While only a subset of children with PLEs develop psychosis, many develop other forms of psychiatric illnesses, including affective disorders (e.g., depression) and post-traumatic stress disorder[5]. This suggests that PLEs may represent early markers of not only psychosis, but transdiagnostic psychiatric illness more broadly.
Sex and gender differences exist in the development and expression of psychiatric illnesses[6], including psychosis[7]. Women (as per their gender, hereafter referred to as ‘women’) and individuals assigned female sex at birth (hereafter referred to as ‘females’) are more likely to exhibit internalizing problems (e.g., loneliness), while men (as per their gender, hereafter referred to as ‘men’) and individuals assigned male sex at birth (hereafter referred to as ‘males’) are more likely to exhibit externalizing problems (e.g., aggression) [8]. Psychosis tends to have an earlier onset in men and males along with more several clinical features relative to women and females[9]. PLEs also differ across sexes and genders such that delusions are more common in women and females while hallucinations are more common in men and males[10]. Some recent studies have examined the neurobiological underpinnings for these phenomenological differences. Functional magnetic resonance imaging is used to indirectly estimate neuronal activation based on changes in blood oxygenation. Regional functional activation signals can be correlated to assess the functional connectivity (or coupling) between brain regions. Disruptions in functional connectivity within the frontoparietal control and default mode networks have been implicated in a wide range of psychiatric illnesses and mental health traits[11–13]. Critically, these functional networks exhibit sex differences[14, 15] and they are uniquely related to transdiagnostic psychiatric symptom domains (e.g., internalizing and externalizing problems) across the sexes in youth[16]. However, it remains to established whether PLEs map onto shared or distinct functional networks across the sexes. An analysis of the sex-specific functional network underpinnings of PLEs in children will reveal early brain-based markers indicative of potential future onset of psychiatric illness that can be used to develop targeted preventative strategies.
Here, we sought to characterize the sex-specific functional network underpinnings of PLEs in children. To do so, we used multivariate analyses to quantify the functional networks associated with the PLEs in a large sample of children from the Adolescent Brain Cognitive Development (ABCD) Study (n=5620, 2671 females, ages 9–10) in a sex-specific manner. First, using a brain-based predictive modeling approach, we demonstrate that functional connectivity is significantly associated with the total number of PLEs and the resulting distress in females and males. Next, evaluating whether the functional networks underlying PLEs are shared or distinct across the sexes, we find that although there is some overlap, a distinct set of networks are implicated in PLEs in females and males. Finally, characterizing the functional network underpinnings of PLEs, we reveal that, in females, a diffuse set of cortical and non-cortical functional connections are associated with PLEs, whereas in males, PLEs predominantly map onto limbic, temporal parietal, somato/motor, and visual networks.
Results
Males report more PLEs and greater distress relative to females:
We assessed sex differences in PLEs in children. PLEs were assessed using the Prodromal Questionnaire-Brief Child Version, which is comprised of 21 questions. Total and Distress scores were computed, consistent with prior literature[3, 17, 18]. In our sample, 57% of females and 60% of males reported one or more PLEs (out of 21), with 18% of females and 20% of males reporting five of more PLEs. Across all females, an average Total score of 2.31±3.43 and Distress score of 5.74±10.40 was reported. Across all males, an average Total score of 2.55±3.54 and Distress score of 5.84±10.06 was reported. Although higher Total (Mann Whitney U statistic, U=3.31×106, pcorrected<0.01) and Distress (U=3.36×106, pcorrected<0.03) scores were reported in males relative to females, the observed differences in the distribution of the Total and Distress scores were minimal (Figure 1A). These results are in line with those previously reported in the ABCD Study[19], indicating that the subsample included here with complete resting-state functional MRI and behavioral data is representative of the larger ABCD cohort. Of note, these data indicate that PLEs are not uncommon in children.
Functional network connectivity is significantly associated with PLEs:
Using sex-specific cross-validated linear ridge regression models, we quantified the sex-specific associations between functional network connectivity and the Total and Distress scores for PLEs. As in similar studies using brain-based predictive models to study associations between functional connectivity and behavioral measures in the ABCD sample[13, 16, 20, 21], model performance was evaluated using measures of prediction accuracy (i.e., correlation between observed and predicted scores). Across both sexes, individual variability in functional connectivity was associated with the Total and Distress scores (Figure 1B). Although, these observed associations were numerically higher in males for both Total (prediction accuracy, rTotal= 0.10) and Distress (rDistress=0.11) scores compared to females (rTotal=0.09; rDistress=0.08), all brain-based prediction models yielded significant results (pcorrected<0.05 for all models relative to null models). These results demonstrate that individual variations in functional connectivity can predict both the endorsement of PLEs and the associated distress in youth. As such, functional network connectivity may serve as a potential early biomarker of transdiagnostic psychiatric illness risk.
Overlapping and distinct functional networks underlie PLEs across the sexes:
We applied the Haufe transformation[22] to the model feature weights to increase their interpretability and reliability[20]. As in prior work[15, 16], we averaged the transformed weights across the multiple train/test splits to obtain a mean feature importance score. To assess whether shared or distinct functional connections were associated with PLEs across the sexes, we computed the full correlation between the features extracted from each sex to predict Total and Distress scores (Figure 2A). Functional connections associated with Total and Distress scores were strongly correlated within each sex (rfemales=0.97, rmales=0.97), while functional connections associated with each score were moderately correlated across the sexes (rtotal=0.59, rdistress= 0.57). Finally, the associations between functional connectivity and the Total and Distress scores were summarized to a network-level by mapping the absolute regional pairwise feature weights onto 17 cortical networks[23] and one non-cortical network. In females, while the strongest associations were observed in temporal parietal, default, limbic, and ventral attention networks, dispersed functional connections throughout the brain were associated with PLEs (Figure 2Bfor Distress and S2 for Total scores). In males, limbic networks exhibited the strongest associations with PLEs followed by weaker associations in temporal parietal, somatomotor, and visual networks (Figure 2Bfor Distress and S2 for Total scores). Summarized network-level associations were also computed (Figure 3) and demonstrate similar results. Prior work in this area has identified associations between resting-state network connectivity in control, default mode, and cinguloparietal networks and PLEs[24]. Here, we find that while those networks are implicated to some extent, other networks (e.g., limbic, ventral attention, somatomotor, visual) also contribute to PLEs, and these associations between connectivity and PLEs differ across males and females. Taken together, this suggests that the functional network underpinnings of PLEs overlap to some extent across the sexes in children but also include distinct network contributions.
Discussion
PLEs represent one of the earliest markers of transdiagnostic psychopathology[25]. An understanding of the neural bases of PLEs is critical for the subsequent development of targeted prevention strategies to mitigate the potential negative sequelae of these experiences in at-risk youth. Here, we demonstrate that individual variations in functional network connectivity are significantly and uniquely associated with PLEs across the sexes in children. First, examining the prevalence in PLEs, we confirm prior work, showing that approximately 3 in 5 children report PLEs, and the number of experiences and the resulting distress is largely comparable across females and males. Next, quantifying the relationships between functional connectivity and PLEs, we demonstrate that individual differences in functional connectivity are significantly associated with PLEs in both sexes. Finally, comparing the network underpinnings of PLEs across the sexes, we find that they are moderately correlated but do exhibit notable differences. In females, dispersed functional connections throughout cortical and non-cortical structures are implicated in PLEs, whereas in males, PLEs are primarily associated with functional connections in limbic, temporal parietal, somatomotor, and visual networks. Collectively, our results suggest that the neurobiological underpinnings of PLEs may, in part, be influenced by assigned sex.
Extant literature has shown that PLEs are related to changes in neurobiological structure and function. PLEs have previously been linked to lower cortical gyrification in left fronto-temporal regions and higher cortical gyrification in right parietal cortex[26]. Reduced cortical and subcortical volumes have also been reported in youth with more distressing PLEs[27]. In the functional domain, connectivity between cerebellar and control networks, between cerebellar and cingulo-parietal networks, and within the control network have been implicated in distressing PLEs[27]. A separate study found that connections within the control, default, and cingulo-parietal networks were associated with PLEs[24]. Unfortunately, these studies only considered assigned sex as a covariate, if at all, restricting their ability to capture critical sex differences in the neurobiological underpinnings of PLEs[28]. Our present analyses demonstrate that the overlapping and unique individual variations in functional networks underlie PLEs across the sexes. Specifically, we observe significant associations between functional connectivity and Total and Distress PLE scores in the range of 0.08 and 0.11. PLEs are phenomenologically heterogeneous, so it is not surprising that these effect sizes, while significant, are relatively small in magnitude. Notably, our observed effect sizes are comparable to those reported between gender and risk-taking behavior (r=0.08) as well as combat exposure and PTSD (r=0.11)[29].
Self-reported PLEs are highly prevalent in children and are associated with parent-report psychopathology[3], delayed functional development[27], and future development of psychiatric illnesses[5]. Although not all youth who report PLEs develop psychosis, the risk of conversion is 3.5 to 4 times higher than individuals without PLEs[30, 31]. Moreover, many of them experience other mental health problems, including anxiety, depression, PTSD, substance dependence, and suicide attempts[5]. In fact, a New Zealand birth cohort study suggests more than 90% of children with strong PLEs have some form by clinically diagnosable psychiatric disorder by age 38[5]. This suggests that childhood PLEs may be predictive of poor mental health during adolescence and adulthood. As such, children who experience PLEs may benefit from early interventions directed at improving their mental well-being[32]. Here, we establish the sex-specific neural bases of PLEs in a large sample, providing a critical foundation for subsequent research on the neurobiological mechanisms driving the future emergence of transdiagnostic psychopathology in these individuals. Importantly, the network associations captured here are, for the most part, unique from those previously reported between functional connectivity and other psychiatric illness-linked behaviors (e.g., internalizing and externalizing) in this sample[16]. This suggests the brain-behavior associations captured here are specific to PLEs and do not represent general psychopathology in children. The neurobiological markers of PLEs identified here can subsequently be validated in smaller, disease-specific samples (e.g., NAPLS [33–35]). The present findings can also be used to guide the development of targeted intervention strategies that may modulate these functional connections, ultimately improving resilience.
Sex differences in the neurobiological underpinnings of PLEs, and psychosis more broadly, have been largely understudied to date. In schizophrenia, men and males have a greater risk of developing negative symptoms (e.g., avolition, anhedonia, alogia), while women and females tend to display more affective symptoms (e.g., depression, impulsivity, emotional instability)[9, 36]. Women and females also exhibit stronger cognitive functioning along executive, verbal, and processing domains, while men and males outperform females along memory and attention domains[36]. Women and females are at a particularly elevated risk of developing psychosis during critical hormonal transition periods (e.g., pregnancy, post-partum, menopause)[37], and earlier menarche is associated with a later onset of schizophrenia[38]. Here, we observe critical sex differences in the associations between functional brain networks and PLEs, providing an insight into the neural bases that may, in part, underlie sex differences in psychosis and schizophrenia. The sex differences observed here may be influenced by a number of sex- and gender-related factors, including differences in functional network maturation trajectories across the sexes[39–41]. Additional analyses using the longitudinal ABCD data and other datasets will inform whether the associations observed here remain consistent throughout adolescence and adulthood.
Importantly, the reported findings are subject to several limitations. First, these results were obtained using a single dataset. Although, our sample is relatively large (n=5,260) and we used a cross-validated, sex-specific brain-based predictive modeling approach known to generate reliable results[16], we cannot rule out the possibility that these results may be limited in their generalizability[42]. Future research in global open-access datasets with comparable neuroimaging and behavioral data could analyze whether the reported findings are specific to the ABCD sample considered here. Second, we only considered the influences of binary assigned sex in these analyses. Although sex is not binary, participants in this sample reported their assigned sex as Female or Male, and as such we were only able to compare the network underpinnings of PLEs between the binary sexes. Additionally, gender has been shown to influence functional network organization[15] and the manifestation of psychiatric illness[43]. Unfortunately, gender identity data was not acquired in this sample at the baseline time point and thus was not considered here. It is possible that our results are, in part, driven by gender rather than sex. Subsequent analyzes considering both sex and gender could disentangle their independent and intersecting influences on the neural bases of PLEs. Finally, we assessed the relationships between functional connectivity and PLEs at a single time-point from the baseline acquisition as we sought to identify the earliest neurobiological markers of transdiagnostic psychopathology. However, the networks implicated in these findings mature throughout adolescence[39]. Therefore, it is possible that the brain markers of PLEs, and poor mental health more broadly, shift throughout the transition from childhood to adolescence to adulthood. Future work incorporating multiple time points and a longitudinal analytical framework will be able to assess these potential changes.
Sex is a fundamental part of human identity and can influence both neurobiology and mental health. An understanding of how sex shapes the relationships between neurobiology and mental health is necessary for the advancement of mental health prevention and intervention strategies that are tailored to each sex. Here, we identify unique functional network markers of PLEs across the sexes, revealing sex-specific early markers that may be broadly indicative of poor mental health and forecast the development of transdiagnostic psychopathology.
Methods
The analyses described here use a similar framework as in our prior work [16, 44–46] to perform novel analyses to establish the functional network underpinnings of psychotic-like experiences in a sex-specific manner.
Dataset:
The Adolescent Brain Cognitive Development (ABCD) dataset is a large community-based sample of children and adolescents who were assessed on a comprehensive set of neuroimaging, behavioral, developmental, and psychiatric batteries[47]. Here, we used minimally preprocessed neuroimaging data and behavioral data acquired at the baseline from the NIMH Data Archive for ABCD Release 2.0.1. Magnetic resonance (MR) images were acquired across 21 sites in the United States of America using harmonized protocols for GE and Siemens scanners. In line with our prior work[13, 48], exclusion criteria were used to ensure quality control. As recommended by the ABCD consortium, we excluded individuals who were scanned using Philips scanners due to incorrect preprocessing (https://github.com/ABCD-STUDY/fMRI-cleanup). For the T1 data, individuals who did not pass recon-all quality control[49] were removed. For the functional connectivity data, functional runs with boundary-based registration (BBR) costs greater than 0.6 were excluded. Further, volumes with framewise displacement (FD)>0.3 mm or voxel-wise differentiated signal variance (DVARS)>50, along with one volume before and two volumes after, were marked as outliers and subsequently censored. Uncensored segments of data containing fewer than five contiguous volumes were also censored[50, 51]. Functional runs with over half of their volumes censored and/or max FD>5mm were removed. Individuals who did not have at least 4 minutes of data were also excluded. Individuals who did not have all of the relevant behavioral data were also excluded. Finally, we excluded siblings to prevent unintended biases due to inherent heritability in neurobiological and/or behavioral measures. Our final sample comprised 5260 children (2571 females, ages 9–10 years).
Sex:
We considered sex in terms of sex assigned at birth (referred to as ‘sex’)[52].
Psychotic-Like Experiences:
We considered the number and severity of psychotic-like experiences as measured by the self-report Prodromal Questionnaire–Brief Child Version. This questionnaire has previously been validated in the ABCD sample [3], and is considered to be a useful measure of early risk for psychosis. Moreover, psychotic-like experiences in the ABCD sample, as measured by this questionnaire, have been shown to be associated with a similar set of familial, cognitive, and emotional factors as in adults with psychosis[3]. Consistent with prior literature[3, 17, 18], we considered Total and Distress scores in these analyses. For each question, participants were first asked if they previously had that experience (1=Yes, 0=No). If they answered yes, they were then asked if the experience bothered them (1=Yes, 0=No). If they again answered yes, they were asked the extent to which the experience had bothered them (1=Not very bothered, 2=Slightly bothered, 3=Moderately bothered, 4=Very much bothered, 5=Extremely bothered). The sum of the endorsed questions corresponding to each PLE (range 0 to 21) was used to compute the Total score, and the total number of endorsed questions was weighted by the level of distress to compute the Distress score (range 0 to 126). Non-parametric Mann-Whitney U rank tests were used to evaluate sex differences in the Total and Distress scores. Resulting p-values were corrected for multiple comparisons across the two scores using the Benjamini-Hochberg False Discovery Rate (q=0.05) procedure[53].
Image Acquisition and Processing:
The minimally preprocessed MRI data were processed as previously described[13, 54]. Briefly, minimally preprocessed T1 data were further processed using FreeSurfer 5.3.0[55–58] to generate cortical surface meshes for each individual, which were then registered to a common spherical coordinate system[57, 58]. Minimally preprocessed fMRI data were further processed with the following steps: (1) removal of initial frames, with the number of frames removed depending on the type of scanner[49] and (2) alignment with the T1 images using boundary-based registration[59] with FsFast. Framewise displacement (FD)[60] and voxel-wise differentiated signal variance (DVARS)[61] were computed using fsl_motion_outliers. Respiratory pseudomotion was filtered out using a bandstop filter (0.31–0.43 Hz) before computing FD[62–64]. A total of 18 nuisance covariates were also regressed out of the fMRI time series: global signal, six motion correction parameters, averaged ventricular signal, averaged white matter signal, and their temporal derivatives. Regression coefficients were estimated from the non-censored volumes. Global signal regression was performed as it has been shown to improve behavioral prediction performance[65, 66]. Finally, the brain scans were interpolated across censored frames using least squares spectral estimation[67], band-pass filtered (0.009 Hz ≤ f ≤ 0.08 Hz), projected onto FreeSurfer fsaverage6 surface space, and smoothed using a 6 mm full-width half maximum kernel. Once processed, we extracted regional time series for 400 cortical[68] and 19 non-cortical[69] parcels. We computed full correlations between those time series yielding a 419×419 pairwise regional functional connectivity matrix. All processing as described was completed on a local server.
Predictive Modeling:
Linear ridge regression models capture robust, reliable, and interpretable associations between neuroimaging and phenotypic data while avoiding data leakage and minimizing overfitting[70, 71]. Here, we used sex-specific cross-validated linear ridge regression models to predict the Total and Distress scores based on functional connectivity in a sex-specific manner. For each sex, we randomly shuffled and split the data into 100 distinct train and test sets (at approximately a 2:1 ratio) without replacement. Imaging site was considered when splitting the data such that we placed all participants from a given site either in the train or test set but not split across the two. Within each train set, we optimized the regularization parameter using three-fold cross-validation while similarly accounting for imaging site as in the initial train-test split. Once optimized, we evaluated models on the corresponding test set. We repeated this process for each of 100 distinct train-test splits to obtain a distribution of prediction accuracy and explained variance. To evaluate model significance, for each set of predictive models, a corresponding set of null models was generated as follows: the output variable was randomly permuted 1000 times, and each permutation was used to train and test a null model using a randomly selected regularization parameter from the set of selected parameters from the original model. Prediction accuracy from each of the null models was then compared to the average accuracy from the corresponding distribution of model accuracies of the original models. The p-value for each model’s significance is defined as the proportion of null models with prediction accuracies or explained variances greater than or equal to those corresponding to the original distributions. All p-values were corrected for multiple comparisons across the two scores using the Benjamini-Hochberg False Discovery Rate (q=0.05) procedure[53].
Feature Weights:
We used the Haufe transformation[22] to transform feature weights obtained from the linear ridge regression models to increase their interpretability and reliability[13, 54, 72]. For each train split, we used feature weights obtained from the model, , the covariance of the input data (functional connectivity), , and the covariance of the output data (behavioral score), , to compute the Haufe-transformed feature weights, , as follows: . We then averaged the Haufe-transformed feature weights across the 100 splits to obtain a mean feature importance value. We computed full correlations between mean feature importance obtained from the sex-specific models to evaluate whether they relied on shared or unique features to predict PLEs. For all models, we also summarized pairwise regional feature importance at a network-level to support interpretability as previously described[45]. To do so, cortical parcels were assigned to one of 17 networks from the Yeo 17-network parcellation[73], and subcortical, brainstem, and cerebellar parcels were assigned to a single non-cortical network for convenience. Absolute regional pairwise feature weights were averaged to yield network-level estimates of associations between functional connectivity and psychotic-like experiences for each sex. Finally, pairwise network-level were summarized to a network-level by taking the sum across the pairwise network connections.
Acknowledgements
Data used in the preparation of this article were obtained from the Adolescent Brain Cognitive DevelopmentSM (ABCD) Study (https://abcdstudy.org), held in the NIMH Data Archive (NDA). This is a multisite, longitudinal study designed to recruit more than 10,000 children age 9–10 and follow them over 10 years into early adulthood. The ABCD Study® is supported by the National Institutes of Health and additional federal partners under award numbers U01DA041048, U01DA050989, U01DA051016, U01DA041022, U01DA051018, U01DA051037, U01DA050987, U01DA041174, U01DA041106, U01DA041117, U01DA041028, U01DA041134, U01DA050988, U01DA051039, U01DA041156, U01DA041025, U01DA041120, U01DA051038, U01DA041148, U01DA041093, U01DA041089, U24DA041123, U24DA041147. A full list of supporters is available at https://abcdstudy.org/federal-partners.html. A listing of participating sites and a complete listing of the study investigators can be found at https://abcdstudy.org/consortium_members/. ABCD consortium investigators designed and implemented the study and/or provided data but did not necessarily participate in the analysis or writing of this report. This manuscript reflects the views of the authors and may not reflect the opinions or views of the NIH or ABCD consortium investigators.
Funding Sources
This work was supported by the following funding sources: Northwell Health Advancing Women in Science and Medicine (Career Development Award to ED and Educational Achievement Award to ED), Feinstein Institutes for Medical Research (Emerging Scientist Award to ED), National Institute of Mental Health (K23MH127300 to JMR, R01MH123245 to AJH, R01MH120080 to AJH and BTTY, R01MH133334 to BTTY, and R01MH109508 to AKM), NUS Yong Loo Lin School of Medicine (NUHSRO/2020/124/TMR/LOA to BTTY), Singapore National Medical Research Council (NMRC) LCG (OFLCG19May-0035 to BTTY), NMRC CTG-IIT (CTGIIT23jan-0001 to BTTY), NMRC STaR (STaR20nov-0003 to BTTY), Singapore Ministry of Health (MOH) Centre Grant (CG21APR1009 to BTTY), Temasek Foundation (TF2223-IMH-01 to BTTY), Wellcome Trust (226716/Z/22/Z to AKM), and BD2 (Breakthrough Discoveries for thriving with Bipolar Disorder; Clinical Coordinating Center Grant to AKM). Any opinions, findings, and conclusions or recommendations expressed here are those of the authors and do not necessarily reflect the views of the funders.
Footnotes
Financial Disclosures
All authors reported no biomedical financial interests or potential conflicts of interest.