Hostname: page-component-74656f9884-j7cmm Total loading time: 0 Render date: 2026-11-06T04:51:17.388Z Has data issue: false hasContentIssue false
  • English
  • Français

Disentangling heterogeneity of psychosis expression in the general population: sex-specific moderation effects of environmental risk factors on symptom networks

Published online by Cambridge University Press:  25 August 2021

Linda T. Betz Open the ORCID record for Linda T. Betz [Opens in a new window] ,
Nora Penzel Open the ORCID record for Nora Penzel [Opens in a new window] ,
Marlene Rosen ,
Kamaldeep Bhui Open the ORCID record for Kamaldeep Bhui [Opens in a new window] ,
Rachel Upthegrove Open the ORCID record for Rachel Upthegrove [Opens in a new window]  and
Joseph Kambeitz Open the ORCID record for Joseph Kambeitz [Opens in a new window]
Linda T. Betz*
Affiliation:
Department of Psychiatry and Psychotherapy, Faculty of Medicine and University Hospital Cologne, University of Cologne, Cologne, Germany
Nora Penzel
Affiliation:
Department of Psychiatry and Psychotherapy, Faculty of Medicine and University Hospital Cologne, University of Cologne, Cologne, Germany Department of Psychiatry and Psychotherapy, Ludwig-Maximilian-University, Munich, Germany
Marlene Rosen
Affiliation:
Department of Psychiatry and Psychotherapy, Faculty of Medicine and University Hospital Cologne, University of Cologne, Cologne, Germany
Kamaldeep Bhui
Affiliation:
Department of Psychiatry, University of Oxford, Oxford, UK Nuffield Department of Primary Care Health Sciences, University of Oxford, Oxford, UK
Rachel Upthegrove
Affiliation:
Institute for Mental Health and Centre for Human Brain Health, University of Birmingham, Birmingham, UK Birmingham Early Intervention Service, Birmingham Women's and Children's NHS Foundation Trust, Birmingham, UK
Joseph Kambeitz
Affiliation:
Department of Psychiatry and Psychotherapy, Faculty of Medicine and University Hospital Cologne, University of Cologne, Cologne, Germany
*
Author for correspondence: Linda T. Betz, E-mail: linda.betz@uk-koeln.de
Rights & Permissions [Opens in a new window]

Abstract

Background

Psychosis expression in the general population may reflect a behavioral manifestation of the risk for psychotic disorder. It can be conceptualized as an interconnected system of psychotic and affective experiences; a so-called ‘symptom network’. Differences in demographics, as well as exposure to adversities and risk factors, may produce substantial heterogeneity in symptom networks, highlighting potential etiological divergence in psychosis risk.

Methods

To explore this idea in a data-driven way, we employed a novel recursive partitioning approach in the 2007 English National Survey of Psychiatric Morbidity (N = 7242). We sought to identify ‘network phenotypes’ by explaining heterogeneity in symptom networks through potential moderators, including age, sex, ethnicity, deprivation, childhood abuse, separation from parents, bullying, domestic violence, cannabis use, and alcohol.

Results

Sex was the primary source of heterogeneity in symptom networks. Additional heterogeneity was explained by interpersonal trauma (childhood abuse and domestic violence) in women and domestic violence, cannabis use, ethnicity in men. Among women, especially those exposed to early interpersonal trauma, an affective loading within psychosis may have distinct relevance. Men, particularly those from minority ethnic groups, demonstrated a strong network connection between hallucinatory experiences and persecutory ideation.

Conclusion

Symptom networks of psychosis expression in the general population are highly heterogeneous. The structure of symptom networks seems to reflect distinct sex-related adversities, etiologies, and mechanisms of symptom-expression. Disentangling the complex interplay of sex, minority ethnic group status, and other risk factors may help optimize early intervention and prevention strategies in psychosis.

Keywords

Network analysispsychopathologyenvironmental exposuresex differencesaffective pathway to psychosisgeneral population

Information

Type
Original Article
Information
Psychological Medicine , Volume 53 , Issue 5 , April 2023 , pp. 1860 - 1869
Check for updates
Copyright
Copyright © The Author(s), 2021. Published by Cambridge University Press

Introduction

Recent research has advanced our understanding of psychosis through so-called ‘symptom networks’, i.e. causal systems of individual interacting experiences and symptoms (Betz et al., Reference Betz, Penzel, Kambeitz-Ilankovic, Rosen, Chisholm and Stainton2020; Hardy, O'Driscoll, Steel, van der Gaag, & van den Berg, Reference Hardy, O'Driscoll, Steel, van der Gaag and van den Berg2020; Isvoranu, Borsboom, van Os, & Guloksuz, Reference Isvoranu, Borsboom, van Os and Guloksuz2016; Isvoranu et al., Reference Isvoranu, Guloksuz, Epskamp, van Os and Borsboom2019, Reference Isvoranu, van Borkulo, Boyette, Wigman, Vinkers and Borsboom2017; Moffa et al., Reference Moffa, Catone, Kuipers, Kuipers, Freeman, Marwaha and Bebbington2017; Murphy, McBride, Fried, & Shevlin, Reference Murphy, McBride, Fried and Shevlin2018; Robinaugh, Hoekstra, Toner, & Borsboom, Reference Robinaugh, Hoekstra, Toner and Borsboom2020). Complex interactions between specific psychotic as well as non-psychotic experiences (e.g. depression and anxiety) in the general population may predate onset of psychosis in clinical settings (Guloksuz et al., Reference Guloksuz, van Nierop, Bak, de Graaf, Ten Have, van Dorsselaer and van Os2016, Reference Guloksuz, van Nierop, Lieb, van Winkel, Wittchen and van Os2015; Kelleher et al., Reference Kelleher, Keeley, Corcoran, Lynch, Fitzpatrick, Devlin and Cannon2012; Linscott & van Os, Reference Linscott and van Os2013; Murphy et al., Reference Murphy, McBride, Fried and Shevlin2018; van Os & Reininghaus, Reference van Os and Reininghaus2016). Additional lines of evidence indicate that there is considerable etiological continuity between subclinical and clinical levels of psychosis (Binbay et al., Reference Binbay, Drukker, Elbi, Tanık, Özkınay, Onay and Alptekin2012; DeRosse & Karlsgodt, Reference DeRosse and Karlsgodt2015; Kelleher & Cannon, Reference Kelleher and Cannon2011; Linscott & van Os, Reference Linscott and van Os2013). Thus, examining the symptom network structure of a transdiagnostic psychosis phenotype, reflecting a behavioral manifestation of risk for psychotic disorder in the general population that blends gradually into clinical syndromes, may help to better understand etiological mechanisms in psychosis and to develop prevention strategies (Bebbington, Reference Bebbington2015; Binbay et al., Reference Binbay, Drukker, Elbi, Tanık, Özkınay, Onay and Alptekin2012; DeRosse & Karlsgodt, Reference DeRosse and Karlsgodt2015; Isvoranu et al., Reference Isvoranu, Borsboom, van Os and Guloksuz2016; Kelleher & Cannon, Reference Kelleher and Cannon2011; Linscott & van Os, Reference Linscott and van Os2013; Robinaugh et al., Reference Robinaugh, Hoekstra, Toner and Borsboom2020; van Os & Reininghaus, Reference van Os and Reininghaus2016).

Importantly, symptomatology and involved etiological mechanisms in psychosis expression are highly variable by specific at risk groups (Bentall, Wickham, Shevlin, & Varese, Reference Bentall, Wickham, Shevlin and Varese2012; Isvoranu et al., Reference Isvoranu, Borsboom, van Os and Guloksuz2016; Linscott & van Os, Reference Linscott and van Os2013; van Os & Reininghaus, Reference van Os and Reininghaus2016). For example, in line with the theory of an affective pathway to psychosis, early traumatic events are strongly associated with connections between affective and psychotic symptomatology (Myin-Germeys & van Os, Reference Myin-Germeys and van Os2007; Upthegrove et al., Reference Upthegrove, Chard, Jones, Gordon-Smith, Forty, Jones and Craddock2015; van Nierop et al., Reference van Nierop, Viechtbauer, Gunther, van Zelst, de Graaf, Ten Have and van Winkel2015). In the presence of heterogeneity, averaged network models of psychosis may obscure important distinctions in relevant etiological mechanisms across specific risk groups (Jones, Mair, Simon, & Zeileis, Reference Jones, Mair, Simon and Zeileis2020; Moriarity, van Borkulo, & Alloy, Reference Moriarity, van Borkulo and Alloy2020). Thus far, however, heterogeneity in symptom networks of psychosis has been either overlooked or addressed in a partial way on a single candidate risk factor (such as sex, cannabis use, or socioeconomic background) at specific thresholds, or using summed environmental risk scores (Betz et al., Reference Betz, Penzel, Kambeitz-Ilankovic, Rosen, Chisholm and Stainton2020; Guloksuz et al., Reference Guloksuz, van Nierop, Bak, de Graaf, Ten Have, van Dorsselaer and van Os2016; Isvoranu et al., Reference Isvoranu, Borsboom, van Os and Guloksuz2016; Wüsten et al., Reference Wüsten, Schlier, Jaya, Fonseca-Pedrero, Peters and Lincoln2018), which lose specificity and relevance for real work prevention and intervention.

The characterization of ‘network phenotypes’ based on a comprehensive set of environmental and demographic factors may explain heterogeneity; that is, the structure of symptomatology is a function of types, combinations, and intensity of etiological loads in psychosis expression (Jones et al., Reference Jones, Mair, Simon and Zeileis2020; Moriarity et al., Reference Moriarity, van Borkulo and Alloy2020). With the goal of characterization of network phenotypes in mind, the current study uses novel work on recursive partitioning, a data-driven, explorative statistical technique that can sequentially extract isolated and combined moderation effects of a large set of environmental and demographic factors on symptom networks, without a priori specification of thresholds or combinations of risk factors (Jones et al., Reference Jones, Mair, Simon and Zeileis2020; Strobl, Malley, & Tutz, Reference Strobl, Malley and Tutz2009). Recursive partitioning identifies network phenotypes that are maximally distinct from each other (Jones et al., Reference Jones, Mair, Simon and Zeileis2020; Zeileis, Hothorn, & Hornik, Reference Zeileis, Hothorn and Hornik2008).

We used recursive partitioning to define meaningful network phenotypes of psychosis expression in the general population, using the 2007 Adult Psychiatric Morbidity in England Survey (APMS; National Centre for Social Research, University of Leicester, 2017). We hypothesized that exposure to environmental risk, if identified as defining a network phenotype, would be characteristically associated with more densely connected symptom networks when compared with samples not exposed to that specific environmental risk (Guloksuz et al., Reference Guloksuz, van Nierop, Bak, de Graaf, Ten Have, van Dorsselaer and van Os2016, Reference Guloksuz, van Nierop, Lieb, van Winkel, Wittchen and van Os2015; Isvoranu et al., Reference Isvoranu, Borsboom, van Os and Guloksuz2016; Lin, Fried, & Eaton, Reference Lin, Fried and Eaton2019; Russell, Keding, He, Li, & Herringa, Reference Russell, Keding, He, Li and Herringa2020). We also aimed to test whether the strength of connections between individual symptoms differed between network phenotypes.

Method

Data analytic strategy

We conducted all analyses in the R language for statistical computing, version 4.0.4. Throughout, we considered a significance level of α = 0.05. Data of the 2007 APMS (National Centre for Social Research, University of Leicester, 2017) used in the analyses are available from the UK Data Service (https://ukdataservice.ac.uk/). Code to reproduce the analyses can be accessed at www.github.com/LindaBetz/APMS_NetworkTree.

Sample

We present analyses based on the 2007 APMS of adults living in private households aged 16 and above who were recruited using a stratified multistage random probability sampling strategy (N = 7403) (McManus, Meltzer, Brugha, Bebbington, & Jenkins, Reference McManus, Meltzer, Brugha, Bebbington and Jenkins2009; Singleton, Bumpstead, O'Brien, Lee, & Meltzer, Reference Singleton, Bumpstead, O'Brien, Lee and Meltzer2003). Methods, procedures, and full details on sample characteristics have been described previously (McManus et al., Reference McManus, Meltzer, Brugha, Bebbington and Jenkins2009). For the present analyses, we excluded participants with missing values in the variables of interest, given that the methods employed do not allow missings. For comparing sample characteristics of included and excluded participants, we used permutation tests as implemented in the R package ‘coin’ (Hothorn, Hornik, van de Wiel, & Zeileis, Reference Hothorn, Hornik, van de Wiel and Zeileis2008).

Assessment of symptomatology

Selection and definition of symptom variables followed a previously published network analysis using data from the 2007 APMS (Moffa et al., Reference Moffa, Catone, Kuipers, Kuipers, Freeman, Marwaha and Bebbington2017), including measures from an affective domain (worry, sleep disturbance, generalized anxiety, and depression), and from a psychotic domain (persecutory ideation and hallucinatory experiences). All symptom variables in the network were coded in binary form (present or absent). For details on these assessments, see online Supplementary Method 1.

Assessment of environmental and demographic risk factors

Environmental risk factors comprised of psychosocial adversities in the form of physical abuse and sexual abuse before the age of 16, separation from parents until the age of 16 (local authority care and/or institutional care), lifetime experiences of bullying, and lifetime experiences of domestic violence. Additionally, we included sex, age, ethnic origin (White, Black, South Asian, and Mixed/Other), cannabis use in the past year, alcohol use, and socioeconomic deprivation. For details on these assessments, see online Supplementary Method 2.

Identification of network subgroups via recursive partitioning

In a first step, we estimated a partial correlation network (without regularization) based on the full sample, using the R package ‘qgraph’, version 1.6.5 (Epskamp, Cramer, Waldorp, Schmittmann, & Borsboom, Reference Epskamp, Cramer, Waldorp, Schmittmann and Borsboom2012). A partial correlation network depicts unique pairwise associations between variables (‘edges’ in network terminology), i.e. the share of the association between two variables that remains after controlling for all other variables in the network (Epskamp, Borsboom, & Fried, Reference Epskamp, Borsboom and Fried2018). We estimated the underlying zero-order correlations between the binary items using Pearson's φ, as recommended when employing recursive partitioning on binary data (Jones et al., Reference Jones, Mair, Simon and Zeileis2020). The stronger the partial correlation between two variables, the more likely it is that they co-occur, controlling for the other variables under consideration.

Second, we used a model-based recursive partitioning approach to identify meaningful subgroups of symptom networks given the included environmental and demographic factors, as implemented in the R package ‘networktree’, version 1.0.1 (Jones et al., Reference Jones, Mair, Simon and Zeileis2020). In brief, recursive partitioning sequentially creates a decision tree by either splitting or not splitting the sample along a set of potential moderating variables (Strobl et al., Reference Strobl, Malley and Tutz2009; Zeileis et al., Reference Zeileis, Hothorn and Hornik2008). The ‘networktree’ approach (Fig. 1) determines sample splits based on significant invariance in the correlation matrix of the network variables under consideration, yielding non-overlapping partitions of the sample with maximally heterogeneous symptom networks (Jones et al., Reference Jones, Mair, Simon and Zeileis2020). For a detailed account, we refer to online Supplementary Method 3 and available methodological articles (Jones et al., Reference Jones, Mair, Simon and Zeileis2020; Strobl et al., Reference Strobl, Malley and Tutz2009; Zeileis et al., Reference Zeileis, Hothorn and Hornik2008). For plotting, we transformed the correlation matrices to partial correlation matrices using the R package ‘qgraph’, such that edges reflect unique associations between two variables.

Fig. 1.

Recursive partitioning for symptom networks as applied to data from the 2007 Adult Psychiatric Morbidity Survey (APMS) study. The goal is to assess which of the included demographic and risk factors capture individual deviations from the correlation matrix of symptom scores, which underlies symptom networks. Starting with the whole sample, individual deviations from the correlation matrix of symptom scores are computed via a log-likelihood-based score function. The variable that explains these deviations best, as determined by a minimum p-value strategy at Bonferroni-corrected α, is selected (here: sex), and the sample split accordingly. Within the identified subgroups, the procedure is repeated recursively until no significant deviations, i.e. heterogeneity, is detected. We compared symptom networks of the identified subgroups in terms of global strength and individual edge weights. For a detailed account of the method, see online Supplementary Method 3 and Jones et al. (Reference Jones, Mair, Simon and Zeileis2020).

Comparison of identified subgroups

To delineate specific network differences between the identified subgroups (i.e. differences between subgroups as defined by a splitting factor in the recursive partitioning approach), we compared the overall strength of symptom connections, defined as the absolute sum of all individual partial correlation coefficients in the network (global strength; S), and differences in estimates of individual partial correlation coefficients (individual edge weights; ρ) within a Bayesian framework, using the R package ‘BGGM’, version 2.0.2 (Williams, Reference Williams2021; Williams & Mulder, Reference Williams and Mulder2019; Williams, Rast, Pericchi, & Mulder, Reference Williams, Rast, Pericchi and Mulder2020). Specifically, we used posterior predictive checks for assessing differences in overall connection strength (Williams et al., Reference Williams, Rast, Pericchi and Mulder2020), and evaluated the posterior distribution for each difference in partial correlation coefficients, where we deemed a difference significant if the 95% credible interval did not contain 0 (Williams, Reference Williams2021). The p-values derived from recursive partitioning are denoted as p RP, whereas p-values derived from post-hoc comparisons implemented in the package ‘BGGM’ are denoted as p BGGM.

Robustness analyses

We used the R package ‘bootnet’, version 1.4.3 (Epskamp et al., Reference Epskamp, Borsboom and Fried2018) to conduct robustness analyses to check stability and accuracy of the results. We investigated stability of symptom networks estimated in the full sample and identified subgroups by testing sensitivity to dropping cases. Specifically, we assessed the degree to which edge weights remained the same after re-estimating the networks with less cases via the correlation stability (CS) coefficient. The CS coefficient represents the maximum proportion of cases that can be dropped, such that the correlation between original edge weights and edge weights of networks based on subsets is 0.7 or higher (95% confidence). The CS coefficient should preferably be above 0.5 (good stability), and not be below 0.25 (acceptable stability) (Epskamp et al., Reference Epskamp, Borsboom and Fried2018). To investigate the accuracy of individual edge weights estimates across the networks in the full sample and identified subgroups, participants were randomly resampled 5000 times, and the bootstrapped confidence intervals (CIs) of the edge weights were estimated.

Results

Sample

Following removal of 161 participants (2.2% of the whole sample) with missing values in the variables of interest, the final sample comprised of 7242 participants, 56.8% of whom were women, with an average age of 50 (IQR = 30) years. Participants excluded due to missing data were on average older, less White and reported lower proportions of alcohol use and hallucinatory experiences, and higher proportions of depressive symptoms (online Supplementary Table S1).

Network variables and potential moderators

Table 1 presents positive endorsement and characteristics of the network variables and potential moderating risk factors in the sample. The most prevalent symptom was worry, and the most prevalent risk factor bullying.

Table 1.

Network variables and potential moderators with positive endorsement (%) or median (IQR) for the whole sample and disaggregated by sex

Overall symptom network structure and subgroups

The partial correlation network estimated in the full sample suggested positive relationships between all symptoms, with a mean edge weight of 0.11. Partial correlations within each symptom domain were, on average, stronger than between the domains. Recursive partitioning revealed that six of the tested demographic and environmental risk factors were linked to significant heterogeneity in symptom networks and split the sample accordingly in a hierarchical fashion: sex, childhood sexual abuse, childhood physical abuse, domestic violence, cannabis use, and ethnicity. Partial correlation matrices for the plotted networks are available at the linked GitHub repository. Sex was the primary source of heterogeneity (Fig. 2a, p RP < 0.001): the network of women was overall significantly less strongly connected (ΔS = −0.17, p BGGM = 0.002), and featured a significantly stronger connection between depression and hallucination (Δρ = 0.06), and a significantly weaker connection between sleep problems and persecutory ideation (Δρ = −0.07) than the network of men. This means that in women, depression and hallucination were more likely to co-occur than in men, whereas sleep problems and persecutory ideation were less likely to co-occur than in men. For networks of women and men, see Fig. 3.

Fig. 2.

Results from recursive partitioning, depicted as a decision tree of partial correlation networks. Numbers behind splitting factors give the sample size retained after the corresponding sample split. Symptom domains are differentiated by color. The thicker and less transparent the edge, the stronger the partial correlation between two symptoms. Blue (red) edges indicate positive (negative) relationships. To ensure visual comparability, edge weights were scaled identically across all networks. Only connections representing edge weights larger than 0.01 are depicted.

Fig. 3.

Partial correlation networks estimated in women (n = 4115) and men (n = 3127). Sex was identified as the first split in the recursive partitioning approach, suggesting that sex was the primary source of heterogeneity in symptom networks. Symptom domains are differentiated by color. The thicker and less transparent the edge, the stronger the partial correlation between two symptoms. Blue (red) edges indicate positive (negative) relationships. To ensure visual comparability, edge weights were scaled identically across both networks.

Distinct risk factors explained further heterogeneity in symptom networks of women and men, yielding eight different network phenotypes in total. Among women, experiences of childhood sexual abuse were the major source of heterogeneity in symptom networks (Fig. 2b, p RP = 0.016) linked to a stronger connection between anxiety and persecutory ideation (Δρ = 0.09). The difference in global strength of the symptom networks of women who reported sexual abuse and those who did not was not significant (ΔS = 0.08, p BGGM = 0.482). Among women who reported no childhood sexual abuse, exposure to childhood physical abuse explained further heterogeneity (Fig. 2c, p RP = 0.015), and was associated with a significantly stronger association between anxiety and hallucinations (Δρ = 0.26). Corresponding symptom networks did not differ significantly in global strength (ΔS = 0.70, p BGGM = 0.204). Finally, among those women that reported neither sexual nor physical abuse, exposure to domestic violence (Fig. 2d, p RP = 0.012) was linked to a stronger connection between worry and depression (Δρ = 0.12), as well as persecutory ideation and hallucinations (Δρ = 0.19). The difference in global strength of the corresponding symptom networks was not significant (ΔS = 0.32, p BGGM = 0.113).

Among men, in those who reported having experienced domestic violence (Fig. 2e, p RP = 0.007) the connection between sleep problems and anxiety was significantly stronger than that in men who did not report past domestic violence (Δρ = 0.17). Global strength was not significantly different between the corresponding networks (ΔS = 0.26, p BGGM = 0.376). Second, in men not reporting past domestic violence, cannabis use in the past year (Fig. 2f, p RP = 0.009) was associated with a significantly increased connection between worry and persecutory ideation (Δρ = 0.17), and a significantly weaker connection between hallucination and persecutory ideation (Δρ = −0.18). Global strength of the corresponding symptom networks did not differ significantly (ΔS = 0.33, p BGGM = 0.126). Finally, men reporting neither domestic violence nor cannabis use were further split by ethnic background (Fig. 2g, p RP = 0.011): the network of men with a Black or South Asian ethnic background was overall significantly more strongly connected (ΔS = 0.76, p BGGM = 0.003), and showed stronger connections between worry and depression (Δρ = 0.25), sleep problems and anxiety (Δρ = 0.20), anxiety and depression (Δρ = 0.16), depression and persecutory ideation (Δρ = 0.21), as well as persecutory ideation and hallucinatory experiences (Δρ = 0.20), and a weaker connection between sleep problems and depression (Δρ = −0.23) than the network of men from a White or Mixed ethnic background.

Age of the respondent, alcohol use, bullying, separation experiences, and socioeconomic deprivation were not identified as relevant sources of heterogeneity in symptom networks. Repeating analyses based on data from women and men separately yielded identical results regarding sex-specific moderators (online Supplementary Fig. S1).

Robustness analyses

The network estimated in the full sample, as well as all identified subgroup networks, showed good stability to dropping cases (online Supplementary Table S2). Accuracy analyses showed some relatively wide bootstrapped CIs in some of the identified subgroups with smaller sample sizes. In these cases, we recommend caution when interpreting the strength of weaker edges. However, the bootstrap mean was generally very close to the sample mean, indicating interpretable results (online Supplementary Figs S2–S16).

Discussion

In the current study, we employed a novel, data-driven recursive partitioning approach in a large national household survey to identify networks of psychotic and affective experiences in the population. Our findings point to considerable heterogeneity, which we explain with several phenotypic systems: six (out of 11) demographic and environmental risk factors yielded eight different network phenotypes, with sex being the primary source of heterogeneity in symptom networks. Among women and men, different risk factors were related to heterogeneity in symptom networks, suggesting potentially distinct relevance and mechanisms of these risk factors across the sexes, in line with a multidimensional model of sexual differentiation in psychosis risk (Riecher-Rössler, Butler, & Kulkarni, Reference Riecher-Rössler, Butler and Kulkarni2018). Overall, our findings on sex and other environmental differences illustrate that the multifactorial and heterogeneous nature of psychosis expression (Isvoranu et al., Reference Isvoranu, Borsboom, van Os and Guloksuz2016; Linscott & van Os, Reference Linscott and van Os2013; van Os & Reininghaus, Reference van Os and Reininghaus2016) appears to be reflected in symptom networks that differed considerably depending on the type, combination, and strength of demographic and environmental risk in a large general population sample.

Differences in symptom networks of women and men

The identification of multiple network phenotypes substantiates the notion that averaged symptom network models are likely not representative of psychosis expression in the general population (Jones et al., Reference Jones, Mair, Simon and Zeileis2020). Rather, observed differences in the strength of overall and specific symptom connections may point to diverse etiological mechanisms operating across different demographic and environmental risk factors. Corroborating a growing recognition that understanding variability by sex is central for the development of comprehensive etiological models of psychopathology (Hartung & Lefler, Reference Hartung and Lefler2019; Hodes & Epperson, Reference Hodes and Epperson2019; Riecher-Rössler et al., Reference Riecher-Rössler, Butler and Kulkarni2018; Rosen, Haidl, Ruhrmann, Vogeley, & Schultze-Lutter, Reference Rosen, Haidl, Ruhrmann, Vogeley and Schultze-Lutter2019), the primary source of heterogeneity in symptom networks of psychosis was sex.

Specifically, our results highlight how associations between affective and psychotic experiences may be differentially expressed in women and men. Prior research indicates that, following the theory of an affective pathway to psychosis, affective alterations, in particular depression and anxiety, may be fundamental driving forces of psychotic experiences (Betz et al., Reference Betz, Penzel, Kambeitz-Ilankovic, Rosen, Chisholm and Stainton2020; Isvoranu et al., Reference Isvoranu, van Borkulo, Boyette, Wigman, Vinkers and Borsboom2017; Myin-Germeys & van Os, Reference Myin-Germeys and van Os2007; Upthegrove et al., Reference Upthegrove, Lalousis, Mallikarjun, Chisholm, Griffiths and Iqbal2020; Upthegrove, Marwaha, & Birchwood, Reference Upthegrove, Marwaha and Birchwood2017; van Nierop et al., Reference van Nierop, Lecei, Myin-Germeys, Collip, Viechtbauer, Jacobs and van Winkel2018). Present findings suggest a particularly strong association between depression and hallucinatory experiences in the network of women compared to men, corroborating the idea that such an affective pathway to psychosis involving depression may be expressed to a greater degree in women, potentially funneled by increased emotional reactivity to life events and daily hassles (Davis, Matthews, & Twamley, Reference Davis, Matthews and Twamley1999; Hodes & Epperson, Reference Hodes and Epperson2019; Myin-Germeys & van Os, Reference Myin-Germeys and van Os2007; Stainton et al., Reference Stainton, Chisholm, Woodall, Hallett, Reniers, Lin and Wood2021). In the symptom network of men, by contrast, a previously identified link between sleep problems and persecutory ideation (Freeman et al., Reference Freeman, Brugha, Meltzer, Jenkins, Stahl and Bebbington2010) was stronger, and therefore, possibly more relevant, than in women. An intriguing potential clinical implication to be tested is that men may, on average, profit in particular from the use of interventions for sleep problems with demonstrated benefit for reducing persecutory ideation (Freeman et al., Reference Freeman, Sheaves, Goodwin, Yu, Nickless, Harrison and Espie2017).

Risk factors explaining heterogeneity in symptom networks of women and men

Among women, heterogeneity in symptom networks of psychosis expression was explained by exposure to interpersonal trauma, including childhood abuse and domestic violence. Specifically, exposure to childhood abuse was linked to stronger associations between anxiety and psychotic experiences. These findings are consistent with previous reports of increased proportions of mixed symptom expression following childhood trauma (Guloksuz et al., Reference Guloksuz, van Nierop, Lieb, van Winkel, Wittchen and van Os2015; Russell et al., Reference Russell, Keding, He, Li and Herringa2020; Upthegrove et al., Reference Upthegrove, Chard, Jones, Gordon-Smith, Forty, Jones and Craddock2015; van Nierop et al., Reference van Nierop, Viechtbauer, Gunther, van Zelst, de Graaf, Ten Have and van Winkel2015), but extend the literature by highlighting how sex may be an important determinant in this relationship. Following trauma, women are more likely to blame themselves, to view the world as dangerous, and to hold more negative views of themselves (Davis et al., Reference Davis, Matthews and Twamley1999; Tolin & Foa, Reference Tolin, Foa and Kimerling2002). This may facilitate a pathway from distressing interpretations of everyday events, including the experience of anxiety, to threat beliefs feeding into the formation of psychotic experiences, as proposed in cognitive models of psychosis (Freeman, Reference Freeman2007; Garety, Kuipers, Fowler, Freeman, & Bebbington, Reference Garety, Kuipers, Fowler, Freeman and Bebbington2001; Hardy et al., Reference Hardy, O'Driscoll, Steel, van der Gaag and van den Berg2020). Overall, the idea that a pathway from anxiety to psychotic experiences may be particularly relevant among women with a history of childhood abuse has potentially important repercussions for clinical practice and deserves further investigation (Bloomfield et al., Reference Bloomfield, Yusuf, Srinivasan, Kelleher, Bell and Pitman2020). Moreover, at a population level, it may well be that links between affective and psychotic experiences following childhood abuse are manifestations of personality function. The interplay between borderline personality functioning and affective instability, also involving psychosis, and subclinical and clinical levels of psychosis warrants further investigations (Barnow et al., Reference Barnow, Arens, Sieswerda, Dinu-Biringer, Spitzer and Lang2010).

Among men, cannabis use and minority ethnic group status were identified as potential sources of heterogeneity in network connections between psychotic and affective symptoms. Most striking differences were evident in the symptom network of men with a minority ethnic group status reporting no domestic violence or cannabis use. Documented variations in experience and reporting of hallucinations (Vanheusden et al., Reference Vanheusden, Mulder, van der Ende, Selten, van Lenthe, Verhulst and Mackenbach2008) and delusions (Berg et al., Reference Berg, Andreassen, Aminoff, Romm, Hauff and Melle2014) in minority ethnic groups seem to extend to the level of symptom networks. Here, they appear to be expressed as an increased co-occurrence of hallucinations and persecutory ideations in men from a minority ethnic background compared to men from the majority White or Mixed ethnic background. This finding agrees with the idea that, under the influence of risk factors, hallucinations and delusions can become connected, which has been linked to worse prognosis and symptom persistence (Binbay et al., Reference Binbay, Drukker, Elbi, Tanık, Özkınay, Onay and Alptekin2012; Smeets et al., Reference Smeets, Lataster, Dominguez, Hommes, Lieb, Wittchen and van Os2012; Smeets, Lataster, Viechtbauer, Delespaul, & G.R.O.U.P., Reference Smeets, Lataster, Viechtbauer and Delespaul2014; van Os & Reininghaus, Reference van Os and Reininghaus2016). Taken together with the present results, this may reinforce evidence that demonstrates that people from a minority ethnic background, particularly men, are at increased risk for poor mental health outcomes (Morgan et al., Reference Morgan, Fearon, Lappin, Heslin, Donoghue, Lomas and Dazzan2017; Singh et al., Reference Singh, Brown, Winsper, Gajwani, Islam, Jasani and Birchwood2015). With the present data, however, it cannot be excluded that ethnicity acts as a proxy measure for factors not covered by our analysis, such as specific forms of deprivation. Delineating how mental health outcomes in men from a minority ethnic background are determined is an outstanding task for future research and may help to design more effective interventions. Identifying potential commonalities underlying minority ethnic group status and domestic violence, both of which were associated with increased co-occurrence of psychotic experiences in men and women, respectively, may prove insightful in this context.

Except for domestic violence, which was a relevant moderating factor in women and men, different risk factors explained heterogeneity in symptom networks of women and men, suggesting a likely complex interplay between sex and risk factors in impacting psychosis risk. Childhood sexual and physical abuse, for instance, were sources of heterogeneity in symptom networks of women, but not men. This finding adds to previous research suggesting particularly detrimental effects of sexual and physical abuse on mental health of girls and women (Adams, Mrug, & Knight, Reference Adams, Mrug and Knight2018; Thompson, Kingree, & Desai, Reference Thompson, Kingree and Desai2004). One reason for the distinct role of adversities may lie in the sex-specific effects they have on the nervous system, against the backdrop of sex differences in maturation, structure, and functioning thereof (DeSantis et al., Reference DeSantis, Baker, Back, Spratt, Ciolino, Moran-Santa Maria and Brady2011; Dow-Edwards, Reference Dow-Edwards2020; Hodes & Epperson, Reference Hodes and Epperson2019; Popovic et al., Reference Popovic, Ruef, Dwyer, Antonucci, Eder and Sanfelici2020). Moreover, characteristics of some risk factors have been shown to differ by sex: men are more likely to engage in more escalating and chronic patterns of cannabis use than women, for example (Hawes, Trucco, Duperrouzel, Coxe, & Gonzalez, Reference Hawes, Trucco, Duperrouzel, Coxe and Gonzalez2019; Wagner & Anthony, Reference Wagner and Anthony2007). Girls, on the other hand, are more likely than boys to experience severe forms of sexual abuse within close victim–perpetrator relationships (Gold, Elhai, Lucenko, Swingle, & Hughes, Reference Gold, Elhai, Lucenko, Swingle and Hughes1998; Kendall-Tackett, Williams, & Finkelhor, Reference Kendall-Tackett, Williams and Finkelhor1993). Such variations may contribute to differing patterns of relationships between risk and symptom expression in women and men.

Overall, our results corroborate a growing realization that research should appraise that mechanisms contributing to psychosis expression may, at least in parts, differ by sex (Hodes & Epperson, Reference Hodes and Epperson2019; Riecher-Rössler et al., Reference Riecher-Rössler, Butler and Kulkarni2018; Rosen et al., Reference Rosen, Haidl, Ruhrmann, Vogeley and Schultze-Lutter2019; Stainton et al., Reference Stainton, Chisholm, Woodall, Hallett, Reniers, Lin and Wood2021). As clinical research works toward early identification and individually tailored preventive interventions, the complex interplay between sex and environmental factors in impacting psychosis risk needs to be better understood to optimize these efforts (Hartung & Lefler, Reference Hartung and Lefler2019; Riecher-Rössler et al., Reference Riecher-Rössler, Butler and Kulkarni2018; Rosen et al., Reference Rosen, Haidl, Ruhrmann, Vogeley and Schultze-Lutter2019; Stainton et al., Reference Stainton, Chisholm, Woodall, Hallett, Reniers, Lin and Wood2021). This includes disaggregating results by sex and gender in psychosis research more consistently (Hartung & Lefler, Reference Hartung and Lefler2019), for example by documenting differences and similarities in symptom networks of women and men.

Limitations

Results from the current study should be interpreted given several limitations. First, posterior predictive checks used for comparing the overall network connectivity tend to be conservative (Williams et al., Reference Williams, Rast, Pericchi and Mulder2020), which may have resulted in low sensitivity in post-hoc comparisons. This factor, and small sample sizes in some subgroups, may explain why we found no evidence that exposure to risk factors was associated with more densely connected symptom networks compared to non-exposure, contrary to our hypothesis. Effects of risk factors on symptom networks seem to be more specific, impacting single relations between symptoms rather than connectivity between all symptoms. Second, model-based recursive partitioning identifies those variables that reduce heterogeneity in symptom networks the most. Thus, age of the respondent, alcohol use, separation experiences, bullying, and socioeconomic deprivation may explain heterogeneity in symptom networks, but not to the same extent as the other risk factors tested. Related, differential relevance of risk factors, for example within ethnic groups, may have remained undetected due to small sample sizes in some subgroups. For a better understanding of the mechanisms relevant in different minority groups, targeted investigations in these populations with larger sample sizes are needed. Third, we did not incorporate complex design features of the APMS, such as weights to take non-response into account, due to the lack of established methods to do so for network models (Lin et al., Reference Lin, Fried and Eaton2019). Related, recursive partitioning currently only allows for complete case analyses. Even though the percentage of excluded participants was small, they differed from included participants in some important aspects, including hallucinatory and depressive symptoms, which may have biased our results. Although therefore not fully representative of the English population, our results are based on a large national household survey, with suitability for a data-intensive method, such as network-based recursive partitioning, unlike for smaller samples which would not offer the same opportunity. Fourth, the retrospective assessment of risk factors via self-report may be prone to memory biases and so directions of effect may be contested. Fifth, data used in the present analyses were gathered in a large epidemiological study; therefore, instruments and tools used were designed such that they were simple to understand and appropriate given their use in over 7000 people. This setting necessarily leads to less refined assessments of symptomatology and risk. Sixth, the analyses were based on cross-sectional data, meaning that the directions of interactions among the symptoms remain unknown. Longitudinal studies are therefore an important next step for this line of research, and extension of recursive partitioning methods to personalized network structures (e.g. derived from experience sampling methods) may allow for insights into how risk factors moderate dynamic associations between symptoms in individuals. Finally, some researchers have expressed concerns about stability and replicability of network models (e.g. estimates of edges; Forbes, Wright, Markon, & Krueger, Reference Forbes, Wright, Markon and Krueger2017; for a summary of the debate, see McNally, Reference McNally2021). Although our robustness analysis suggests that the networks and edge estimates are generally stable, especially weaker links in the networks of small subgroups should be interpreted with care. Given that recursive partitioning and network methodology are data-driven, replication of present findings in other samples is needed to establish generalizability (Fried et al., Reference Fried, Eidhof, Palic, Costantini, Huisman-van Dijk, Bockting and Karstoft2018).

Conclusion

Symptom networks of psychosis expression in the general population are highly heterogeneous. Sex was the primary source of heterogeneity, and different risk factors explained further variability in symptom networks of women and men, potentially reflecting distinct sex-specific mechanisms contributing to psychosis risk. Among women, an affective loading within psychosis, particularly following early interpersonal trauma, may have distinct importance. Among men, the symptom network of those from a minority ethnic background showed a particularly strong connection between hallucinatory experiences and persecutory ideation, which may reflect poorer outcomes including symptom resolution in this group. A better understanding and consideration of these sex differences provides an important opportunity to deliver high-quality prevention and patient-centered care in psychosis.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0033291721003470

Financial support

JK has received funding from the German Research Foundation (DFG; grant agreement no. KA 4413/1-1). MR has received funding from the Köln Fortune program (grant agreement no. 304/2020).

Conflict of interest

The authors declare no conflict of interests with relation to the work reported in this manuscript.

Ethical standards

The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008.

References

National Centre for Social Research, University of Leicester (2017). Adult psychiatric morbidity survey, 2007 (4th ed., Version 6379) [Data set]. UK Data Service. https://doi.org/10.5255/UKDA-SN-6379-2.Google Scholar
Adams, J., Mrug, S., & Knight, D. C. (2018). Characteristics of child physical and sexual abuse as predictors of psychopathology. Child Abuse & Neglect, 86, 167–177. https://doi.org/10.1016/j.chiabu.2018.09.019.CrossRefGoogle ScholarPubMed
Barnow, S., Arens, E. A., Sieswerda, S., Dinu-Biringer, R., Spitzer, C., & Lang, S. (2010). Borderline personality disorder and psychosis: A review. Current Psychiatry Reports, 12(3), 186–195. https://doi.org/10.1007/s11920-010-0107-9.CrossRefGoogle ScholarPubMed
Bebbington, P. (2015). Unravelling psychosis: Psychosocial epidemiology, mechanism, and meaning. Shanghai Archives of Psychiatry, 27(2), 70. https://doi.org/10.11919/j.issn.1002-0829.215027.Google ScholarPubMed
Bentall, R. P., Wickham, S., Shevlin, M., & Varese, F. (2012). Do specific early-life adversities lead to specific symptoms of psychosis? A study from the 2007 the adult psychiatric morbidity survey. Schizophrenia Bulletin, 38(4), 734–740. https://doi.org/10.1093/schbul/sbs049.CrossRefGoogle ScholarPubMed
Berg, A. O., Andreassen, O. A., Aminoff, S. R., Romm, K. L., Hauff, E., & Melle, I. (2014). The impact of immigration and visible minority status on psychosis symptom profile. Social Psychiatry and Psychiatric Epidemiology, 49(11), 1747–1757. https://doi.org/10.1007/s00127-014-0897-x.CrossRefGoogle ScholarPubMed
Betz, L. T., Penzel, N., Kambeitz-Ilankovic, L., Rosen, M., Chisholm, K., Stainton, A., … PRONIA Consortium (2020). General psychopathology links burden of recent life events and psychotic symptoms in a network approach. npj Schizophrenia, 6(1), 40. https://doi.org/10.1038/s41537-020-00129-w.CrossRefGoogle Scholar
Binbay, T., Drukker, M., Elbi, H., Tanık, F. A., Özkınay, F., Onay, H., … Alptekin, K. (2012). Testing the psychosis continuum: Differential impact of genetic and nongenetic risk factors and comorbid psychopathology across the entire spectrum of psychosis. Schizophrenia Bulletin, 38(5), 992–1002. https://doi.org/10.1093/schbul/sbr003.CrossRefGoogle ScholarPubMed
Bloomfield, M. A. P., Yusuf, F. N. I. B., Srinivasan, R., Kelleher, I., Bell, V., & Pitman, A. (2020). Trauma-informed care for adult survivors of developmental trauma with psychotic and dissociative symptoms: A systematic review of intervention studies. The Lancet Psychiatry, 7(5), 449–462. https://doi.org/10.1016/S2215-0366(20)30041-9.CrossRefGoogle ScholarPubMed
Davis, M. C., Matthews, K. A., & Twamley, E. W. (1999). Is life more difficult on Mars or Venus? A meta-analytic review of sex differences in major and minor life events. Annals of Behavioral Medicine, 21(1), 83–97. https://doi.org/10.1007/BF02895038.CrossRefGoogle ScholarPubMed
DeRosse, P., & Karlsgodt, K. H. (2015). Examining the psychosis continuum. Current Behavioral Neuroscience Reports, 2(2), 80–89. https://doi.org/10.1007/s40473-015-0040-7.CrossRefGoogle ScholarPubMed
DeSantis, S. M., Baker, N. L., Back, S. E., Spratt, E., Ciolino, J. D., Moran-Santa Maria, M., … Brady, K. T. (2011). Gender differences in the effect of early life trauma on hypothalamic-pituitary-adrenal axis functioning. Depression and Anxiety, 28(5), 383–392. https://doi.org/10.1002/da.20795.CrossRefGoogle ScholarPubMed
Dow-Edwards, D. (2020). Sex differences in the interactive effects of early life stress and the endocannabinoid system. Neurotoxicology and Teratology, 80, 106893. https://doi.org/10.1016/j.ntt.2020.106893.CrossRefGoogle ScholarPubMed
Epskamp, S., Borsboom, D., & Fried, E. I. (2018). Estimating psychological networks and their accuracy: A tutorial paper. Behavior Research Methods, 50(1), 195–212. https://doi.org/10.3758/s13428-017-0862-1.CrossRefGoogle ScholarPubMed
Epskamp, S., Cramer, A. O. J., Waldorp, L. J., Schmittmann, V. D., & Borsboom, D. (2012). qgraph: Network visualizations of relationships in psychometric data. Journal of Statistical Software, 48(4), 1–18. https://doi.org/10.18637/jss.v048.i04.CrossRefGoogle Scholar
Forbes, M. K., Wright, A. G. C., Markon, K. E., & Krueger, R. F. (2017). Evidence that psychopathology symptom networks have limited replicability. Journal of Abnormal Psychology, 126(7), 969–988. https://doi.org/10.1037/abn0000276.CrossRefGoogle ScholarPubMed
Freeman, D. (2007). Suspicious minds: The psychology of persecutory delusions. Clinical Psychology Review, 27(4), 425–457. https://doi.org/10.1016/j.cpr.2006.10.004.CrossRefGoogle ScholarPubMed
Freeman, D., Brugha, T., Meltzer, H., Jenkins, R., Stahl, D., & Bebbington, P. (2010). Persecutory ideation and insomnia: Findings from the second British national survey of psychiatric morbidity. Journal of Psychiatric Research, 44(15), 1021–1026. https://doi.org/10.1016/j.jpsychires.2010.03.018.CrossRefGoogle ScholarPubMed
Freeman, D., Sheaves, B., Goodwin, G. M., Yu, L.-M., Nickless, A., Harrison, P. J., … Espie, C. A. (2017). The effects of improving sleep on mental health (OASIS): A randomised controlled trial with mediation analysis. The Lancet. Psychiatry, 4(10), 749–758. https://doi.org/10.1016/S2215-0366(17)30328-0.CrossRefGoogle ScholarPubMed
Fried, E. I., Eidhof, M. B., Palic, S., Costantini, G., Huisman-van Dijk, H. M., Bockting, C. L. H., … Karstoft, K.-I. (2018). Replicability and generalizability of posttraumatic stress disorder (PTSD) networks: A cross-cultural multisite study of PTSD symptoms in four trauma patient samples. Clinical Psychological Science, 6(3), 335–351. https://doi.org/10.1177/2167702617745092.CrossRefGoogle ScholarPubMed
Garety, P. A., Kuipers, E., Fowler, D., Freeman, D., & Bebbington, P. E. (2001). A cognitive model of the positive symptoms of psychosis. Psychological Medicine, 31(2), 189–195. https://doi.org/10.1017/s0033291701003312.CrossRefGoogle ScholarPubMed
Gold, S. N., Elhai, J. D., Lucenko, B. A., Swingle, J. M., & Hughes, D. M. (1998). Abuse characteristics among childhood sexual abuse survivors in therapy: A gender comparison. Child Abuse & Neglect, 22(10), 1005–1012. https://doi.org/10.1016/S0145-2134(98)00085-4.CrossRefGoogle ScholarPubMed
Guloksuz, S., van Nierop, M., Bak, M., de Graaf, R., Ten Have, M., van Dorsselaer, S., … van Os, J. (2016). Exposure to environmental factors increases connectivity between symptom domains in the psychopathology network. BMC Psychiatry, 16, 223. https://doi.org/10.1186/s12888-016-0935-1.CrossRefGoogle ScholarPubMed
Guloksuz, S., van Nierop, M., Lieb, R., van Winkel, R., Wittchen, H.-U., & van Os, J. (2015). Evidence that the presence of psychosis in non-psychotic disorder is environment-dependent and mediated by severity of non-psychotic psychopathology. Psychological Medicine, 45(11), 2389–2401. https://doi.org/10.1017/S0033291715000380.CrossRefGoogle ScholarPubMed
Hardy, A., O'Driscoll, C., Steel, C., van der Gaag, M., & van den Berg, D. (2020). A network analysis of post-traumatic stress and psychosis symptoms. Psychological Medicine, 1–8. https://doi.org/10.1017/S0033291720001300.Google ScholarPubMed
Hartung, C. M., & Lefler, E. K. (2019). Sex and gender in psychopathology: DSM-5 and beyond. Psychological Bulletin, 145(4), 390–409. https://doi.org/10.1037/bul0000183.CrossRefGoogle ScholarPubMed
Hawes, S. W., Trucco, E. M., Duperrouzel, J. C., Coxe, S., & Gonzalez, R. (2019). Developmental pathways of adolescent cannabis use: Risk factors, outcomes and sex-specific differences. Substance Use & Misuse, 54(2), 271–281. https://doi.org/10.1080/10826084.2018.1517177.CrossRefGoogle ScholarPubMed
Hodes, G. E., & Epperson, C. N. (2019). Sex differences in vulnerability and resilience to stress across the life span. Biological Psychiatry, 86(6), 421–432. https://doi.org/10.1016/j.biopsych.2019.04.028.CrossRefGoogle ScholarPubMed
Hothorn, T., Hornik, K., van de Wiel, M., & Zeileis, A. (2008). Implementing a class of permutation tests: The coin package. Journal of Statistical Software, 28(8), 1–23. https://doi.org/10.18637/jss.v028.i08.CrossRefGoogle Scholar
Isvoranu, A.-M., Borsboom, D., van Os, J., & Guloksuz, S. (2016). A network approach to environmental impact in psychotic disorder: Brief theoretical framework. Schizophrenia Bulletin, 42(4), 870–873. https://doi.org/10.1093/schbul/sbw049.CrossRefGoogle ScholarPubMed
Isvoranu, A.-M., Guloksuz, S., Epskamp, S., van Os, J., Borsboom, D., & GROUP Investigators (2019). Toward incorporating genetic risk scores into symptom networks of psychosis. Psychological Medicine, 50(4), 1–8. https://doi.org/10.1017/S003329171900045X.Google ScholarPubMed
Isvoranu, A.-M., van Borkulo, C. D., Boyette, L.-L., Wigman, J. T. W., Vinkers, C. H., Borsboom, D. … Group Investigators (2017). A network approach to psychosis: Pathways between childhood trauma and psychotic symptoms. Schizophrenia Bulletin, 43(1), 187–196. https://doi.org/10.1093/schbul/sbw055.CrossRefGoogle ScholarPubMed
Jones, P. J., Mair, P., Simon, T., & Zeileis, A. (2020). Network trees: A method for recursively partitioning covariance structures. Psychometrika, 85(4), 926–945. https://doi.org/10.1007/s11336-020-09731-4.CrossRefGoogle Scholar
Kelleher, I., & Cannon, M. (2011). Psychotic-like experiences in the general population: Characterizing a high-risk group for psychosis. Psychological Medicine, 41(1), 1–6. https://doi.org/10.1017/S0033291710001005.CrossRefGoogle Scholar
Kelleher, I., Keeley, H., Corcoran, P., Lynch, F., Fitzpatrick, C., Devlin, N., … Cannon, M. (2012). Clinicopathological significance of psychotic experiences in non-psychotic young people: Evidence from four population-based studies. The British Journal of Psychiatry: The Journal of Mental Science, 201(1), 26–32. https://doi.org/10.1192/bjp.bp.111.101543.CrossRefGoogle ScholarPubMed
Kendall-Tackett, K. A., Williams, L. M., & Finkelhor, D. (1993). Impact of sexual abuse on children: A review and synthesis of recent empirical studies. Psychological Bulletin, 113(1), 164–180. https://doi.org/10.1037/0033-2909.113.1.164.CrossRefGoogle ScholarPubMed
Lin, S.-Y., Fried, E. I., & Eaton, N. R. (2019). The association of life stress with substance use symptoms: A network analysis and replication. Journal of Abnormal Psychology, 129(2), 204–214. https://doi.org/10.1037/abn0000485.CrossRefGoogle ScholarPubMed
Linscott, R. J., & van Os, J. (2013). An updated and conservative systematic review and meta-analysis of epidemiological evidence on psychotic experiences in children and adults: On the pathway from proneness to persistence to dimensional expression across mental disorders. Psychological Medicine, 43(6), 1133–1149. https://doi.org/10.1017/S0033291712001626.CrossRefGoogle ScholarPubMed
McManus, S., Meltzer, H., Brugha, T. S., Bebbington, P. E., & Jenkins, R. (2009). Adult psychiatric morbidity in England, 2007: Results of a household survey. Leeds: The NHS Information Centre for Health and Social Care. https://doi.org/10.13140/2.1.1563.5205.Google Scholar
McNally, R. J. (2021). Network analysis of psychopathology: Controversies and challenges. Annual Review of Clinical Psychology, 17, 31–53. https://doi.org/10.1146/annurev-clinpsy-081219-092850.CrossRefGoogle ScholarPubMed
Moffa, G., Catone, G., Kuipers, J., Kuipers, E., Freeman, D., Marwaha, S., … Bebbington, P. (2017). Using directed acyclic graphs in epidemiological research in psychosis: An analysis of the role of bullying in psychosis. Schizophrenia Bulletin, 43(6), 1273–1279. https://doi.org/10.1093/schbul/sbx013.CrossRefGoogle ScholarPubMed
Morgan, C., Fearon, P., Lappin, J., Heslin, M., Donoghue, K., Lomas, B., … Dazzan, P. (2017). Ethnicity and long-term course and outcome of psychotic disorders in a UK sample: The ÆSOP-10 study. The British Journal of Psychiatry: The Journal of Mental Science, 211(2), 88–94. https://doi.org/10.1192/bjp.bp.116.193342.CrossRefGoogle Scholar
Moriarity, D. P., van Borkulo, C., & Alloy, L. B. (2020). Inflammatory phenotype of depression symptom structure: A network perspective. Brain, Behavior, and Immunity, 93, 35–42. https://doi.org/10.1016/j.bbi.2020.12.005.CrossRefGoogle ScholarPubMed
Murphy, J., McBride, O., Fried, E., & Shevlin, M. (2018). Distress, impairment and the extended psychosis phenotype: A network analysis of psychotic experiences in an US general population sample. Schizophrenia Bulletin, 44(4), 768–777. https://doi.org/10.1093/schbul/sbx134.CrossRefGoogle Scholar
Myin-Germeys, I., & van Os, J. (2007). Stress-reactivity in psychosis: Evidence for an affective pathway to psychosis. Clinical Psychology Review, 27(4), 409–424. https://doi.org/10.1016/j.cpr.2006.09.005.CrossRefGoogle ScholarPubMed
Popovic, D., Ruef, A., Dwyer, D. B., Antonucci, L. A., Eder, J., Sanfelici, R., … PRONIA Consortium (2020). Traces of trauma: A multivariate pattern analysis of childhood trauma, brain structure, and clinical phenotypes. Biological Psychiatry, 88(11), 829–842. https://doi.org/10.1016/j.biopsych.2020.05.020.CrossRefGoogle ScholarPubMed
Riecher-Rössler, A., Butler, S., & Kulkarni, J. (2018). Sex and gender differences in schizophrenic psychoses-a critical review. Archives of Women's Mental Health, 21(6), 627–648. https://doi.org/10.1007/s00737-018-0847-9.CrossRefGoogle ScholarPubMed
Robinaugh, D. J., Hoekstra, R. H. A., Toner, E. R., & Borsboom, D. (2020). The network approach to psychopathology: A review of the literature 2008–2018 and an agenda for future research. Psychological Medicine, 50(3), 353–366. https://doi.org/10.1017/S0033291719003404.CrossRefGoogle Scholar
Rosen, M., Haidl, T. K., Ruhrmann, S., Vogeley, K., & Schultze-Lutter, F. (2019). Sex differences in symptomatology of psychosis-risk patients and in prediction of psychosis. Archives of Women's Mental Health, 23, 339–349. https://doi.org/10.1007/s00737-019-01000-3.CrossRefGoogle ScholarPubMed
Russell, J. D., Keding, T. J., He, Q., Li, J. J., & Herringa, R. J. (2020). Childhood exposure to interpersonal violence is associated with greater transdiagnostic integration of psychiatric symptoms. Psychological Medicine, 1–9. https://doi.org/10.1017/S0033291720003712.Google ScholarPubMed
Singh, S. P., Brown, L., Winsper, C., Gajwani, R., Islam, Z., Jasani, R., … Birchwood, M. (2015). Ethnicity and pathways to care during first episode psychosis: The role of cultural illness attributions. BMC Psychiatry, 15, 287. https://doi.org/10.1186/s12888-015-0665-9.CrossRefGoogle ScholarPubMed
Singleton, N., Bumpstead, R., O'Brien, M., Lee, A., & Meltzer, H. (2003). Psychiatric morbidity among adults living in private households, 2000. International Review of Psychiatry, 15(1–2), 65–73. https://doi.org/10.1080/0954026021000045967.CrossRefGoogle ScholarPubMed
Smeets, F., Lataster, T., Dominguez, M.-D.-G., Hommes, J., Lieb, R., Wittchen, H.-U., & van Os, J. (2012). Evidence that onset of psychosis in the population reflects early hallucinatory experiences that through environmental risks and affective dysregulation become complicated by delusions. Schizophrenia Bulletin, 38(3), 531–542. https://doi.org/10.1093/schbul/sbq117.CrossRefGoogle ScholarPubMed
Smeets, F., Lataster, T., Viechtbauer, W., Delespaul, P., & G.R.O.U.P. (2014). Evidence that environmental and genetic risks for psychotic disorder may operate by impacting on connections between core symptoms of perceptual alteration and delusional ideation. Schizophrenia Bulletin, 41(3), 687–697. https://doi.org/10.1093/schbul/sbu122.CrossRefGoogle ScholarPubMed
Stainton, A., Chisholm, K., Woodall, T., Hallett, D., Reniers, R. L. E. P., Lin, A., & Wood, S. J. (2021). Gender differences in the experience of psychotic-like experiences and their associated factors: A study of adolescents from the general population. Schizophrenia Research, 228, 410–416. https://doi.org/10.1016/j.schres.2021.01.008.CrossRefGoogle ScholarPubMed
Strobl, C., Malley, J., & Tutz, G. (2009). An introduction to recursive partitioning: Rationale, application, and characteristics of classification and regression trees, bagging, and random forests. Psychological Methods, 14(4), 323–348. https://doi.org/10.1037/a0016973.CrossRefGoogle ScholarPubMed
Thompson, M. P., Kingree, J. B., & Desai, S. (2004). Gender differences in long-term health consequences of physical abuse of children: Data from a nationally representative survey. American Journal of Public Health, 94(4), 599–604. https://doi.org/10.2105/ajph.94.4.599.CrossRefGoogle ScholarPubMed
Tolin, D. F., & Foa, E. B. (2002). Gender and PTSD: A cognitive model. In Kimerling, R. (Ed.), Gender and PTSD (Vol. 460, pp. 76–97). New York, NY, USA: The Guilford Press. Retrieved from https://psycnet.apa.org/fulltext/2003-04410-003.pdf.Google Scholar
Upthegrove, R., Chard, C., Jones, L., Gordon-Smith, K., Forty, L., Jones, I., & Craddock, N. (2015). Adverse childhood events and psychosis in bipolar affective disorder. The British Journal of Psychiatry: The Journal of Mental Science, 206(3), 191–197. https://doi.org/10.1192/bjp.bp.114.152611.CrossRefGoogle ScholarPubMed
Upthegrove, R., Lalousis, P., Mallikarjun, P., Chisholm, K., Griffiths, S. L., Iqbal, M., … PRONIA Consortium (2020). The psychopathology and neuroanatomical markers of depression in early psychosis. Schizophrenia Bulletin, 47(1), 249–258. https://doi.org/10.1093/schbul/sbaa094.CrossRefGoogle Scholar
Upthegrove, R., Marwaha, S., & Birchwood, M. (2017). Depression and schizophrenia: Cause, consequence, or trans-diagnostic issue? Schizophrenia Bulletin, 43(2), 240–244. https://doi.org/10.1093/schbul/sbw097.Google ScholarPubMed
Vanheusden, K., Mulder, C. L., van der Ende, J., Selten, J.-P., van Lenthe, F. J., Verhulst, F. C., & Mackenbach, J. P. (2008). Associations between ethnicity and self-reported hallucinations in a population sample of young adults in The Netherlands. Psychological Medicine, 38(8), 1095–1102. https://doi.org/10.1017/S0033291707002401.CrossRefGoogle Scholar
van Nierop, M., Lecei, A., Myin-Germeys, I., Collip, D., Viechtbauer, W., Jacobs, N., … van Winkel, R. (2018). Stress reactivity links childhood trauma exposure to an admixture of depressive, anxiety, and psychosis symptoms. Psychiatry Research, 260, 451–457. https://doi.org/10.1016/j.psychres.2017.12.012.CrossRefGoogle Scholar
van Nierop, M., Viechtbauer, W., Gunther, N., van Zelst, C., de Graaf, R., Ten Have, M., … van Winkel, R. (2015). Childhood trauma is associated with a specific admixture of affective, anxiety, and psychosis symptoms cutting across traditional diagnostic boundaries. Psychological Medicine, 45(6), 1277–1288. https://doi.org/10.1017/S0033291714002372.CrossRefGoogle ScholarPubMed
van Os, J., & Reininghaus, U. (2016). Psychosis as a transdiagnostic and extended phenotype in the general population. World Psychiatry, 15(2), 118–124. https://doi.org/10.1002/wps.20310.CrossRefGoogle ScholarPubMed
Wagner, F. A., & Anthony, J. C. (2007). Male-female differences in the risk of progression from first use to dependence upon cannabis, cocaine, and alcohol. Drug and Alcohol Dependence, 86(2–3), 191–198. https://doi.org/10.1016/j.drugalcdep.2006.06.003.CrossRefGoogle ScholarPubMed
Williams, D. R. (2021). Bayesian estimation for Gaussian graphical models: Structure learning, predictability, and network comparisons. Multivariate Behavioral Research, 56(2), 1–17. https://doi.org/10.1080/00273171.2021.1894412.CrossRefGoogle ScholarPubMed
Williams, D. R., & Mulder, J. (2019). BGGM: A R package for Bayesian Gaussian Graphical Models. Preprint. Retrieved from https://psyarxiv.com/3b5hf/download?format=pdf.Google Scholar
Williams, D. R., Rast, P., Pericchi, L. R., & Mulder, J. (2020). Comparing Gaussian graphical models with the posterior predictive distribution and Bayesian model selection. Psychological Methods, 25(5), 653–672. https://doi.org/10.1037/met0000254.CrossRefGoogle ScholarPubMed
Wüsten, C., Schlier, B., Jaya, E. S., Genetic Risk and Outcome of Psychosis (GROUP) Investigators, Fonseca-Pedrero, E., Peters, E., … Lincoln, T. M. (2018). Psychotic experiences and related distress: A cross-national comparison and network analysis based on 7141 participants from 13 countries. Schizophrenia Bulletin, 44(6), 1185–1194. https://doi.org/10.1093/schbul/sby087.CrossRefGoogle ScholarPubMed
Zeileis, A., Hothorn, T., & Hornik, K. (2008). Model-based recursive partitioning. Journal of Computational and Graphical Statistics, 17(2), 492–514. https://doi.org/10.1198/106186008X319331.CrossRefGoogle Scholar
Figure 0

Fig. 1. Recursive partitioning for symptom networks as applied to data from the 2007 Adult Psychiatric Morbidity Survey (APMS) study. The goal is to assess which of the included demographic and risk factors capture individual deviations from the correlation matrix of symptom scores, which underlies symptom networks. Starting with the whole sample, individual deviations from the correlation matrix of symptom scores are computed via a log-likelihood-based score function. The variable that explains these deviations best, as determined by a minimum p-value strategy at Bonferroni-corrected α, is selected (here: sex), and the sample split accordingly. Within the identified subgroups, the procedure is repeated recursively until no significant deviations, i.e. heterogeneity, is detected. We compared symptom networks of the identified subgroups in terms of global strength and individual edge weights. For a detailed account of the method, see online Supplementary Method 3 and Jones et al. (2020).

Figure 1

Table 1. Network variables and potential moderators with positive endorsement (%) or median (IQR) for the whole sample and disaggregated by sex

Figure 2

Fig. 2. Results from recursive partitioning, depicted as a decision tree of partial correlation networks. Numbers behind splitting factors give the sample size retained after the corresponding sample split. Symptom domains are differentiated by color. The thicker and less transparent the edge, the stronger the partial correlation between two symptoms. Blue (red) edges indicate positive (negative) relationships. To ensure visual comparability, edge weights were scaled identically across all networks. Only connections representing edge weights larger than 0.01 are depicted.

Figure 3

Fig. 3. Partial correlation networks estimated in women (n = 4115) and men (n = 3127). Sex was identified as the first split in the recursive partitioning approach, suggesting that sex was the primary source of heterogeneity in symptom networks. Symptom domains are differentiated by color. The thicker and less transparent the edge, the stronger the partial correlation between two symptoms. Blue (red) edges indicate positive (negative) relationships. To ensure visual comparability, edge weights were scaled identically across both networks.

Supplementary material: File

Betz et al. supplementary material

Betz et al. supplementary material

Download Betz et al. supplementary material(File)
File 5.4 MB