Hostname: page-component-586b7cd67f-t7fkt Total loading time: 0 Render date: 2024-11-22T13:52:32.910Z Has data issue: false hasContentIssue false
  • English
  • Français

Psychotic experiences and negative symptoms from adolescence to emerging adulthood: developmental trajectories and associations with polygenic scores and childhood characteristics

Published online by Cambridge University Press:  03 October 2022

Laura Havers Open the ORCID record for Laura Havers [Opens in a new window] ,
Sophie von Stumm ,
Alastair G. Cardno ,
Daniel Freeman  and
Angelica Ronald
Laura Havers*
Affiliation:
Department of Psychological Sciences, Birkbeck, University of London, London, UK
Sophie von Stumm
Affiliation:
Department of Education, University of York, Heslington, UK
Alastair G. Cardno
Affiliation:
Division of Psychological and Social Medicine, University of Leeds, Leeds, UK
Daniel Freeman
Affiliation:
Department of Psychiatry, University of Oxford, Oxford, UK Oxford Health NHS Foundation Trust, Oxford, UK
Angelica Ronald
Affiliation:
Department of Psychological Sciences, Birkbeck, University of London, London, UK
*
Author for correspondence: Laura Havers, E-mail: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Background

Psychotic experiences and negative symptoms (PENS) are common in non-clinical populations. PENS are associated with adverse outcomes, particularly when they persist. Little is known about the trajectories of PENS dimensions in young people, nor about the precursory factors associated with these trajectories.

Methods

We conducted growth mixture modelling of paranoia, hallucinations, and negative symptoms across ages 16, 17, and 22 in a community sample (N = 12 049–12 652). We then described the emergent trajectory classes through their associations with genome-wide polygenic scores (GPS) for psychiatric and educational phenotypes, and earlier childhood characteristics.

Results

Three trajectory classes emerged for paranoia, two for hallucinations, and two for negative symptoms. Across PENS, GPS for clinical help-seeking, major depressive disorder, and attention deficit hyperactivity disorder were associated with increased odds of being in the most elevated trajectory class (OR 1.07–1.23). Lower education GPS was associated with the most elevated trajectory class for hallucinations and negative symptoms (OR 0.77–0.91). Conversely for paranoia, higher education GPS was associated with the most elevated trajectory class (OR 1.25). Trajectory class associations were not significant for schizophrenia, obsessive-compulsive disorder, bipolar disorder, or anorexia GPS. Emotional/behaviour problems and life events in childhood were associated with increased odds of being in the most elevated trajectory class across PENS.

Conclusions

Our results suggest latent heterogeneity in the development of paranoia, hallucinations, and negative symptoms in young people that is associated with specific polygenic scores and childhood characteristics.

Keywords

Adolescencecommunitydevelopmental trajectoriesemerging adulthoodgrowth mixture modellingnegative symptomspolygenic scorespsychotic experiences
Type
Original Article
Information
Psychological Medicine , Volume 53 , Issue 12 , September 2023 , pp. 5685 - 5697
Check for updates
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2022. Published by Cambridge University Press

Introduction

Experiences that are characteristic of psychosis are frequently observed in community samples (Healy et al., Reference Healy, Brannigan, Dooley, Coughlan, Clarke, Kelleher and Cannon2019; McGrath et al., Reference McGrath, Saha, Al-Hamzawi, Alonso, Bromet, Bruffaerts and Kessler2015). Psychotic experiences such as delusions, paranoia and hallucinations, and negative symptoms, including flat affect and poverty of speech, have been reported in childhood and adolescence, and across the lifespan (Barragan, Laurens, Navarro, & Obiols, Reference Barragan, Laurens, Navarro and Obiols2011; Dhossche, Ferdinand, Ende, Hofstra, & Verhulst, Reference Dhossche, Ferdinand, Ende, Hofstra and Verhulst2002; Dominguez, Saka, Lieb, Wittchen, & van Os, Reference Dominguez, Saka, Lieb, Wittchen and van Os2010; Kelleher et al., Reference Kelleher, Connor, Clarke, Devlin, Harley and Cannon2012; Maric, Krabbendam, Vollebergh, de Graaf, & van Os, Reference Maric, Krabbendam, Vollebergh, de Graaf and van Os2003; Ronald et al., Reference Ronald, Sieradzka, Cardno, Haworth, McGuire and Freeman2014). Whilst psychotic experiences reported by individuals in community samples generally abate (Linscott & van Os, Reference Linscott and van Os2013; van Os, Linscott, Myin-Germeys, Delespaul, & Krabbendam, Reference van Os, Linscott, Myin-Germeys, Delespaul and Krabbendam2009), persistence compared to the transience of psychotic experiences and negative symptoms (PENS) has been found to predict poor clinical and functional outcomes (e.g. Hielscher et al. Reference Hielscher, DeVylder, Hasking, Connell, Martin and Scott2021; Kaymaz et al. Reference Kaymaz, Drukker, Lieb, Wittchen, Werbeloff, Weiser and van Os2012; Mackie, Castellanos-Ryan, & Conrod, Reference Mackie, Castellanos-Ryan and Conrod2011; Wigman et al. Reference Wigman, van Winkel, Raaijmakers, Ormel, Verhulst, Reijneveld and Vollebergh2011b). For example, psychotic experiences reported in an adolescent community sample were associated with increased odds for psychotic impairment in a dose–response manner, with odds ratios of 1.5, 5, and 9.9 for the presence of psychotic experiences at one, two, and three time-points, respectively (Dominguez, Wichers, Lieb, Wittchen, & van Os, Reference Dominguez, Wichers, Lieb, Wittchen and van Os2011). Less is known about the developmental course of negative symptoms in the community, but findings from a small number of studies suggest that persistence is associated with adverse outcomes including psychotic and functional impairment (Dominguez et al., Reference Dominguez, Saka, Lieb, Wittchen and van Os2010; Janssens et al., Reference Janssens, Boyette, Heering, Bartels-Velthuis, Lataster, Kahn and Myin-Germeys2016).

Previous longitudinal studies have often focussed on total scores or aggregated measures of psychotic experiences. Persistence of these experiences has been found to be associated with family background characteristics such as familial psychiatric history and socioeconomic status (SES); and with characteristics reported in childhood such as lower educational attainment, childhood trauma and other adverse life events, and emotional and behavioural problems (Bourque, Afzali, O'Leary-Barrett, & Conrod, Reference Bourque, Afzali, O'Leary-Barrett and Conrod2017; Cougnard et al., Reference Cougnard, Marcelis, Myin-Germeys, Graaf, Vollebergh, Krabbendam and van Os2007; DeVylder, Lehmann, & Chen, Reference DeVylder, Lehmann and Chen2015; Janssens et al., Reference Janssens, Boyette, Heering, Bartels-Velthuis, Lataster, Kahn and Myin-Germeys2016; Kalman, Bresnahan, Schulze, & Susser, Reference Kalman, Bresnahan, Schulze and Susser2019; Rammos et al., Reference Rammos, Sullivan, Kounali, Jones, Hammerton, Hines and Zammit2021).

However, PENS that are reported in the community have a multidimensional psychometric structure (Ronald et al., Reference Ronald, Sieradzka, Cardno, Haworth, McGuire and Freeman2014; Stefanis et al., Reference Stefanis, Delespaul, Henquet, Bakoula, Stefanis and van Os2004; Yung et al., Reference Yung, Nelson, Baker, Buckby, Baksheev and Cosgrave2009), and specific dimensions show distinct associations with other types of psychopathology (Armando et al., Reference Armando, Nelson, Yung, Ross, Birchwood, Girardi and Nastro2010; Ronald et al., Reference Ronald, Sieradzka, Cardno, Haworth, McGuire and Freeman2014; Wigman et al., Reference Wigman, Vollebergh, Raaijmakers, Iedema, van Dorsselaer, Ormel and van Os2011c; Yung et al., Reference Yung, Nelson, Baker, Buckby, Baksheev and Cosgrave2009), environmental exposures (Cosgrave et al., Reference Cosgrave, Purple, Haines, Porcheret, van Heugten-van der Kloet, Johns and Wulff2021; Shakoor et al., Reference Shakoor, McGuire, Cardno, Freeman, Plomin and Ronald2015), and genetic factors (Zavos et al., Reference Zavos, Freeman, Haworth, McGuire, Plomin, Cardno and Ronald2014). Gaining a better understanding of both the development of separate PENS dimensions, and the precursory factors associated with their development is important. It will allow for the delineation of dimension-specific theoretical models, which may be implemented to identify and help individuals at risk for both concurrent psychopathology and later poor outcomes (Armando et al., Reference Armando, Nelson, Yung, Ross, Birchwood, Girardi and Nastro2010; Cosgrave et al., Reference Cosgrave, Purple, Haines, Porcheret, van Heugten-van der Kloet, Johns and Wulff2021; Steenkamp et al., Reference Steenkamp, Tiemeier, Blanken, Hillegers, Kushner and Bolhuis2021; Yung et al., Reference Yung, Nelson, Baker, Buckby, Baksheev and Cosgrave2009).

The handful of studies that have investigated the development of paranoia/delusions and hallucinations separately have mainly relied on measures of a few items (Connell et al., Reference Connell, Betts, McGrath, Alati, Najman, Clavarino and Scott2016; De Loore et al., Reference De Loore, Gunther, Drukker, Feron, Sabbe, Deboutte and Myin-Germeys2011; Hielscher et al., Reference Hielscher, DeVylder, Hasking, Connell, Martin and Scott2021; Sheaves et al., Reference Sheaves, Bebbington, Goodwin, Harrison, Espie, Foster and Freeman2016; Steenkamp et al., Reference Steenkamp, Tiemeier, Blanken, Hillegers, Kushner and Bolhuis2021). These studies found the persistence of paranoia/delusions and hallucinations to be associated with a range of maladies. Building on the findings of these studies that used brief measurement tools, measures that assess quantitative variation across a broader range of experiences (e.g. Bartels-Velthuis, van de Willige, Jenner, van Os, & Wiersma, Reference Bartels-Velthuis, van de Willige, Jenner, van Os and Wiersma2011; Freeman et al. Reference Freeman, Stahl, McManus, Meltzer, Brugha, Wiles and Bebbington2012) have the potential to further enhance understanding of PENS by capturing experiences that may otherwise go undetected (Mitchell et al., Reference Mitchell, Maybery, Russell-Smith, Collerton, Gignac and Waters2017).

In terms of genetic factors, several studies have investigated their influence on PENS at single time-points or assessments (Ronald & Pain, Reference Ronald and Pain2018), and findings from a small number of family studies further suggest that genetic factors are associated with the development of PENS (Ericson, Tuvblad, Raine, Young-Wolff, & Baker, Reference Ericson, Tuvblad, Raine, Young-Wolff and Baker2011; Havers, Taylor, & Ronald, Reference Havers, Taylor and Ronald2019; Janssens et al., Reference Janssens, Boyette, Heering, Bartels-Velthuis, Lataster, Kahn and Myin-Germeys2016; Wigman et al., Reference Wigman, van Winkel, Jacobs, Wichers, Derom, Thiery and van Os2011a). The prior study with the largest sample size (N = 1448 twin pairs) found that 38–62% of the covariance in separate PENS dimensions measured across two time-points in adolescence was accounted for by genetic influences (Havers et al., Reference Havers, Taylor and Ronald2019). Genome-wide polygenic scores (GPS) can also be used as an index of an individual's polygenic propensity to a given outcome. Only one study to date has employed polygenic score methods in the context of the persistence of PENS. This study reported a null association between schizophrenia GPS and the persistence of aggregated psychotic experiences, measured across ages 12–24 in the Avon Longitudinal Study of Parents and Children (Rammos et al., Reference Rammos, Sullivan, Kounali, Jones, Hammerton, Hines and Zammit2021). PENS likely reflect vulnerability for poor functional and clinical outcomes, broadly, rather than solely for psychosis (Healy et al., Reference Healy, Brannigan, Dooley, Coughlan, Clarke, Kelleher and Cannon2019; van Os & Reininghaus, Reference van Os and Reininghaus2016; Yung et al., Reference Yung, Nelson, Baker, Buckby, Baksheev and Cosgrave2009): the extent to which the development of PENS dimensions is associated with polygenic liability across a broad range of phenotypes, including psychiatric disorders, clinical help-seeking, intelligence and educational attainment, is therefore of interest, but has not previously been tested.

Previous studies that have investigated trajectories of PENS have done so primarily by manually grouping individuals according to the presence or absence of these experiences across time-points (e.g. Cougnard et al. Reference Cougnard, Marcelis, Myin-Germeys, Graaf, Vollebergh, Krabbendam and van Os2007; Dominguez et al. Reference Dominguez, Wichers, Lieb, Wittchen and van Os2011; Hafeez & Yung, Reference Hafeez and Yung2021; Hielscher et al. Reference Hielscher, DeVylder, Hasking, Connell, Martin and Scott2021; Janssens et al. Reference Janssens, Boyette, Heering, Bartels-Velthuis, Lataster, Kahn and Myin-Germeys2016; Rammos et al. Reference Rammos, Sullivan, Kounali, Jones, Hammerton, Hines and Zammit2021; Steenkamp et al. Reference Steenkamp, Tiemeier, Blanken, Hillegers, Kushner and Bolhuis2021; van Rossum, Dominguez, Lieb, Wittchen, & van Os, Reference van Rossum, Dominguez, Lieb, Wittchen and van Os2011). In contrast, latent variable modelling [specifically here, growth mixture modelling (GMM)] can be used, figuratively, to investigate whether individuals can be classified according to similarities in their underlying, latent trajectories (Herle et al., Reference Herle, Micali, Abdulkadir, Loos, Bryant-Waugh, Hübel and De Stavola2020). Of the studies that have used GMM, multiple latent trajectory classes including a persistent or increasing class have been identified for broadly-defined psychotic experiences measured across adolescence (Bourque et al., Reference Bourque, Afzali, O'Leary-Barrett and Conrod2017; Lin et al., Reference Lin, Wigman, Nelson, Vollebergh, van Os, Baksheev and Yung2011; Mackie et al., Reference Mackie, Castellanos-Ryan and Conrod2011, Reference Mackie, O'Leary-Barrett, Al-Khudhairy, Castellanos-Ryan, Struve, Topper and Conrod2013; Thapar et al., Reference Thapar, Heron, Jones, Owen, Lewis and Zammit2012; Wigman et al., Reference Wigman, Vollebergh, Raaijmakers, Iedema, van Dorsselaer, Ormel and van Os2011c), and adulthood (Wigman et al., Reference Wigman, van Winkel, Jacobs, Wichers, Derom, Thiery and van Os2011a). Yet, there are currently no published findings that estimate the latent development of paranoia and hallucinations separately, or of negative symptoms, reported in the community.

The current study builds on prior research by modelling latent heterogeneity in the development of separate PENS dimensions. We focused on the period of mid-late adolescence to emerging adulthood, which is a common period of onset for a range of mental health problems, including psychosis (Kessler et al., Reference Kessler, Amminger, Aguilar-Gaxiola, Alonso, Lee and Ustun2007; Kim-Cohen et al., Reference Kim-Cohen, Caspi, Moffitt, Harrington, Milne and Poulton2003; Maibing et al., Reference Maibing, Pedersen, Benros, Mortensen, Dalsgaard and Nordentoft2015). As pre-registered (https://osf.io/pax6k), we hypothesised that multiple classes, including a persistent class, would be identified for each of paranoia, hallucinations, and negative symptoms, and that persistence would be associated with the following: (i) lower SES and family psychiatric history, (ii) more emotional and behavioural difficulties, more life events, and lower educational attainment (both in childhood and adulthood), and (iii) higher GPS of psychiatric and clinical help-seeking outcomes, and lower GPS of intelligence and educational attainment. We also predicted that male sex would be associated with persistent negative symptoms (e.g. Dominguez et al. Reference Dominguez, Saka, Lieb, Wittchen and van Os2010; Roy, Maziade, Labbé, & Mérette, Reference Roy, Maziade, Labbé and Mérette2001).

Methods

Participants

Participants were part of the Twins Early Development Study (TEDS). Families of twins born 1994–1996 in England and Wales were invited by the Office for National Statistics to take part on behalf of TEDS. Sixteen thousand eight hundred and ten (16 810) families responded to the invitation. Online Supplementary Table S1 shows participation rates and exclusion details (see Rimfeld et al., Reference Rimfeld, Malanchini, Spargo, Spickernell, Selzam, McMillan and Plomin2019 for a recent overview of TEDS). Paranoia and hallucinations (as well as the additional measures at age 22) were self-reported, and negative symptoms were reported by parents (N = 12 049–12 652) at mean ages 16.32 (s.d. 0.69), 17.06 (s.d. 0.88), and 22.85 (s.d. 0.88) (online Supplementary Table S2). Online Supplementary Tables S4–S6 show demographic comparisons between individuals with complete and incomplete PENS data. Earlier in the study, parents completed assessments of their children's behaviour at mean age 7.06 (s.d. 0.25), and teachers reported educational achievement at mean age 7.20 (s.d. 0.27).

Measures

PENS

Paranoia, hallucinations, and negative symptoms were assessed using subscales of the Specific Psychotic Experiences Questionnaire (SPEQ; Ronald et al., Reference Ronald, Sieradzka, Cardno, Haworth, McGuire and Freeman2014), described in online Supplementary Information S1. Paranoia was measured by 15 items adapted from the Paranoia Checklist (Freeman et al., Reference Freeman, Garety, Bebbington, Smith, Rollinson, Fowler and Dunn2005), and hallucinations by nine items adapted from the Cardiff Anomalous Perceptions Scale (Bell, Halligan, & Ellis, Reference Bell, Halligan and Ellis2006), both on a 6-point scale and both adapted for use in adolescents by clinical experts (Ronald et al., Reference Ronald, Sieradzka, Cardno, Haworth, McGuire and Freeman2014). Negative symptoms were measured by eight items on a 4-point scale, adapted from the Scale for the Assessment of Negative Symptoms (Andreasen, Reference Andreasen1982). Descriptive statistics are reported in online Supplementary Table S2. Reliability and validity information regarding use of the SPEQ in the current sample is reported in Ronald et al. (Reference Ronald, Sieradzka, Cardno, Haworth, McGuire and Freeman2014).

Additional measures

At ages 7 and 22, the Strengths and Difficulties Questionnaire (SDQ; Goodman, Reference Goodman1997) was used to assess emotional and behavioural problems, and the Short Mood and Feeling Questionnaire (MFQ; Angold, Costello, Messer, and Pickles, Reference Angold, Costello, Messer and Pickles1995) was used at age 22 to assess depressive symptoms. The assessment measures for SES, family psychiatric history, educational attainment, and life events are described in online Supplementary Information S2.

GPS

Genotyping of participants is described in online Supplementary Information S3. GPSs were calculated by other TEDS researchers (Selzam et al., Reference Selzam, McAdams, Coleman, Carnell, O'Reilly, Plomin and Llewellyn2018, Reference Selzam, Ritchie, Pingault, Reynolds, O'Reilly and Plomin2019). Descriptions of the GPSs and their calculation are detailed in online Supplementary Information S4. Standardized residuals of the GPS regressed on the first 10 principal components of ancestry, batch, and chip were used. GPS corresponding to the most predictive fraction (f) of causal markers were used in the association analyses (online Supplementary Information S4 and Tables S14, S24, S34).

Statistical analyses

A structural equation modelling framework was used for all analyses, using Mplus (version 8.6). Observed total score data at each age was modelled using full information maximum likelihood to accommodate missing data under the assumption that data were missing at random, and robust estimation was used to accommodate multivariate nonnormality of residuals. Data from related individuals were accommodated by adjustment to the standard errors, with the family ID specified as the unit of clustering.

Prior to conducting GMM, longitudinal measurement invariance was assessed, and a series of latent growth curve models (LGCM) were run to determine the optimal functional form of growth. GMM was used to identify latent trajectory classes using the growth form suggested by the LGCM. GMM extends LGCM by parameterising the trajectories of a prespecified number of heterogenous latent classes (online Supplementary Fig. S1) (Ram & Grimm, Reference Ram and Grimm2009).

A 1-class model was first estimated, followed by an increasing number of classes, up to the point where there was consistent nonconvergence (Jung & Wickrama, Reference Jung and Wickrama2008). Within each k-class (where k refers to the number of classes), two models were initially run: A latent class growth analysis model (Model LCGA, with no variance on the latent growth factors), and an unconstrained GMM (Model 0). The following parameters were freely estimated in Model 0: (i) growth factor means, (ii) growth factor variances and covariances, (iii) residual variances. For k > 1 models, these parameters were estimated for each class. A series of constrained models were run where there were convergence issues with these models (Sijbrandij et al., Reference Sijbrandij, Hoekstra, Almansa, Peeters, Bültmann and Reijneveld2020; van de Schoot, Sijbrandij, Winter, Depaoli, & Vermunt, Reference van de Schoot, Sijbrandij, Winter, Depaoli and Vermunt2017).

For the subsequently tested models, variance parameter constraints were as follows: Model 1A: Within-class residual variances. Model 1B: Between-class residual variances. Model 1C: Between-class growth factor variances. Model 2A: Within-class and between-class residual variances. Model 2B: Within-class residual variances and between-class growth factor variances. Model 2C: Between-class residual variances and between-class growth factor variances.

Online Supplementary Fig. S2 shows a decision-making flowchart regarding estimation specifications and model adjustments. Variation in time-scores (ages) was modelled by allowing individual (random) slope factor loadings (Mehta & West, Reference Mehta and West2000), specified using TSCORES in Mplus.

For each k-class, the model with the lowest BIC value was selected. These models were rerun using the two seed values corresponding to the highest replicated loglikelihood value (Jung & Wickrama, Reference Jung and Wickrama2008; Shireman, Steinley, & Brusco, Reference Shireman, Steinley and Brusco2016). Where BIC values were indistinguishable (difference <10), the AIC was referred to and the model with the lowest AIC value (difference >2) was selected. The best fitting overall model was determined by jointly considering, (i) BIC (and AIC where necessary), (ii) entropy values (with a value of one reflecting perfect classification accuracy), (iii) empirical and theoretical plausibility of the within-class parameter estimates and overall solution. Two sets of post-hoc sensitivity tests were conducted, described in online Supplementary Information S5. Likelihood ratio tests for comparing k-class with k-1-class GMMs are not available for models with individual time-scores.

Multinomial logistic regression analyses were conducted to assess the relationship between the family background, age 7, and GPS variables (specified as auxiliary variables), and the latent class variable, using the automatically implemented ‘3-step’ procedure in Mplus (Asparouhov & Muthén, Reference Asparouhov and Muthén2014; Vermunt, Reference Vermunt2017). The automatically implemented ‘BCH’ procedure was used to estimate the class-specific means of all auxiliary variables (Asparouhov & Muthén, Reference Asparouhov and Muthén2014; Bakk & Vermunt, Reference Bakk and Vermunt2016; Bolck, Croon, & Hagenaars, Reference Bolck, Croon and Hagenaars2004): Both methods estimate the GMM and the most likely class membership values, and then make estimates of the relations between the auxiliary variables and the latent class variable, adjusted for classification error. The False Discovery Rate method was used to correct for multiple testing at q < 0.05 (Benjamini & Hochberg, Reference Benjamini and Hochberg1995). Single-predictor regressions were first run. Significant predictors (at q < 0.05) were entered into multiple-predictor regressions.

Results

Descriptive statistics

Descriptive statistics for PENS at each age are shown in online Supplementary Table S2. Cross-age correlation coefficients are shown in online Supplementary Table S3.

Measurement invariance

Partial scalar-level invariance was found for paranoia and hallucinations, and partial strict-level invariance was found for negative symptoms (online Supplementary Tables S7, S18, S28).

LGCM

Linear growth models provided an acceptable approximation of the data across PENS (online Supplementary Tables S8, S19, S29; parameter estimates shown in online Supplementary Tables S9, S20, S30). Quadratic and latent basis models provided a better fit than the linear model for hallucinations, though these models were highly constrained (to achieve just- and over-identification). Linear models were taken forward in the interest of parsimony and consistency, to reduce the risk of overfitting, and for simpler estimation using individual time-scores (Sterba, Reference Sterba2014). Online Supplementary Figs S3–S5 show spaghetti plots for observed individual trajectories and mean trajectories estimated in the LGCMs.

GMM

Model fitting results for converged models are shown in Table 1. Estimated parameters and trajectory plots from the best fitting models are shown in Table 2 and Fig. 1, respectively.

Table 1. Growth mixture model fit results for converged models of paranoia, hallucinations, and negative symptoms

Note. k, number of classes; Par., number of estimated parameters for final model; LL, log-likelihood value; AIC, Akaike's information criterion; BIC, Bayesian information criterion. 1 = No calculation for 1-class model. Bold typeset indicates lowest BIC value for each dimension. a = slope factor variance (and covariance) fixed to zero for class #3. Only converged models that were run are shown in this table. All models that were run are shown in online Supplementary Tables S10, S21, S31.

Table 2. Parameter estimates of best-fitting model for paranoia, hallucinations, and negative symptoms

Note. k, number of classes; W1–W3, data collection waves 1–3.

Variance of W1–W3 represents residual variance at data collection waves 1–3. Covariance represents covariance between the intercept and slope factors. Model 0: Unconstrained model. Model 1A: Model with within-class residual variances constrained.

a Parameter fixed to zero.

Fig. 1. Trajectory plots of best-fitting models for paranoia, hallucinations, and negative symptoms.

Note. Lines represent mean trajectories; bands represent 95% confidence intervals. Parameter estimates for the trajectories are reported in Table 2.

Paranoia (online Supplementary Tables S10, S12, S13)

The 3-class Model 0 had the lowest BIC of all best-fitting k-class models. Sensitivity analyses suggested better fit of this model compared to the 2-class Models 1A-2C (online Supplementary Table S11), and 3-class homoscedastic model (df = 16, loglikelihood = −77 087.757, BIC = 154 325.861, entropy = 0.686). Entropy was notably less than one across all k-class models. Most individuals' latent trajectories were characterised by mid-level (56.40%, ‘mid-decreasing’) or low-level (20.66%, ‘low-decreasing’) paranoia scores that decreased over time. For a smaller percentage (22.92%, ‘high-persistent’), paranoia was higher at baseline and persisted over time. Variability around the growth factors was significant for the mid-decreasing and high-persistent latent classes. For the low-decreasing class, slope factor variance was fixed to zero (to avoid singularity of the information matrix), and intercept variability was nonsignificant. Constraining the slope factor means across classes resulted in worse fit (loglikelihood = −74 342.311, BIC = 148 891.350, entropy = 0.651) than the model with freely estimated slopes. Wald tests of the differences between slopes were significant [mid-decreasing, high-persistent: W (1) = 32.188, p < 0.001; mid-decreasing, low-decreasing: W (1) = 14.139, p < 0.001; low-decreasing, high-persistent: W (1) = 12.704, p < 0.001].

Hallucinations (online Supplementary Tables S21, S22, S23)

The 2-class Model 1A had the lowest BIC of all best-fitting k-class models. Entropy was less than one across all models. The 2-class model suggests a decreasing developmental course of hallucinations across individuals, with significant variability around the growth factors for both classes. Constraining the slope factor means across classes resulted in worse fit (loglikelihood = −54 431.913, BIC = 108 976.593, entropy = 0.763) than the model with freely estimated slopes. For individuals classified in the ‘mid-decreasing’ (45.16%) compared to ‘low-decreasing’ (54.84%) class, hallucinations decreased at a significantly greater rate (W (1) = 905.142, p < 0.001).

Negative symptoms (online Supplementary Tables S31–S33)

The 2-class Model 0 had the lowest BIC of all best-fitting k-class models. Entropy was higher for the 2-class models than for the 3-class models. A homoscedastic 2-class model did not fit the data as well (df = 12, loglikelihood = −55 901.090, BIC = 111 915.528, entropy = 0.557) as the unconstrained model. The 2-class model suggests an increasing developmental course of negative symptoms across individuals. For a marginal majority (55.07%, ‘mid-increasing’), negative symptoms were estimated at a mid-level at baseline that increased over time. For others (45.93%, ‘low-increasing’), negative symptoms were estimated at a lower initial level that increased over time. There was significant growth factor variability, though slope factor variability was nonsignificant for the low-increasing class. Constraining slope factor means across classes resulted in worse fit (loglikelihood = −50 745.690, BIC = 101 642.510, AIC = 101 523.381, entropy = 0.784) than the model with freely estimated slopes. Negative symptoms increased at a greater rate in the mid-increasing class compared to the low-increasing class [W (1) = 18.243, p < 0.001].

Regressions and class-specific means

Whilst we did not identify a ‘persistent’ trajectory for hallucinations, nor for negative symptoms per se, we discuss our regression results in terms of the ‘most elevated’ trajectory class (i.e. high-persistent for paranoia, mid-decreasing for hallucinations, and mid-increasing for negative symptoms) for ease of communication.

GPS (online Supplementary Tables S15, S25, S35)

In the single-predictor models, the GPSs for clinical help seeking (psychiatrist and general practitioner), major depressive disorder, and attention deficit hyperactivity disorder were significantly associated with increased odds of membership in the most elevated trajectory class compared to the reference class across PENS, as predicted. Similarly, an increase in GPS for years of education (and intelligence, for negative symptoms) was associated with decreased odds for being in the most elevated trajectory class for hallucinations and negative symptoms.

Against predictions, for paranoia, an increase in the GPSs for years of education and intelligence was associated with increased odds for being in the most elevated trajectory class. The GPS for autism spectrum disorder was associated with increased odds of membership in the most elevated trajectory class for paranoia and hallucinations, but not for negative symptoms. The GPSs for schizophrenia, obsessive-compulsive disorder, bipolar disorder, and anorexia were not significantly associated with latent trajectory class membership for any PENS.

In the multiple-predictor models, for paranoia, the GPSs for years of education, intelligence, major depressive disorder, autism spectrum disorder, and attention deficit hyperactivity disorder remained significant, predicting increased odds of membership in the high-persistent class compared to the low-decreasing class. For hallucinations, only the GPS for autism spectrum disorder remained significant. For negative symptoms, only the GPS for years of education remained significant.

Family background characteristics (online Supplementary Tables S16, S26, S36)

Family history of schizophrenia and bipolar disorder were both associated with increased odds of being in the most elevated trajectory class for paranoia and hallucinations, but not negative symptoms. As predicted, male sex was associated with membership in the most elevated trajectory class for negative symptoms. As predicted for hallucinations and negative symptoms, lower SES was associated with membership in the most elevated trajectory class. For paranoia, opposite to our hypotheses, higher SES was associated with increased odds of being in the high-persistent trajectory class compared to the low-decreasing class. In the multiple-predictor regressions, all associations remained significant for hallucinations and negative symptoms, and all except the family history of bipolar disorder remained significant for paranoia.

Age 7 characteristics (online Supplementary Tables S17, S27, S37)

For all PENS, more life events, and higher emotional and behavioural problems (SDQ scores) were associated with increased odds of being in the most elevated trajectory class, as was lower educational attainment for hallucinations and negative symptoms, as predicted. However, for paranoia, higher educational attainment was associated with increased odds of being in the most elevated trajectory class. In the multiple-predictor models, SDQ remained significantly associated with PENS class membership. For paranoia, higher educational attainment also remained significant.

Age 22 characteristics

For hallucinations and negative symptoms, class-specific means for educational attainment were lower, and class-specific means for life events, depressive symptoms (MFQ scores) and SDQ scores were higher in the most elevated trajectory class compared to the reference class, as predicted. For paranoia, the same pattern of results was found except that the class-specific mean for educational attainment was higher in the high-persistent class compared to the low-decreasing class.

Table 3 also shows the class-specific means for the GPS, background, and age 7 variables. These results mirror the main findings from the multinominal regression analyses.

Table 3. Characteristics of latent trajectory classes for paranoia, hallucinations, and negative symptoms

Note. N indicates the number of individuals with data contributing to the GMM and not missing on the auxiliary variable. Related and unrelated individuals included, using cluster-robust s.e. For binary variables (1), the mean represents the proportion. Diff. reflects the chi-square value (df 1) of the difference between the means (or proportions), significant at FDR-adjusted q < 0.05 unless indicated by ‘ = ’.

GPS, genome-wide polygenic score; GPSEDU, years of education; GPSIQ, intelligence; GPSPSYCH, ever visited a psychiatrist for nerves, anxiety, tension, or depression; GPSGP, ever visited a general practitioner for nerves, anxiety, tension, or depression; GPSOCD, obsessive-compulsive disorder; GPSMDD, major depressive disorder; GPSASD, autism spectrum disorder; GPSADHD, attention deficit hyperactivity disorder; SES, socioeconomic status; Family SCZ, family history of schizophrenia; Family BIP, family history of bipolar disorder; Ed attainment, educational attainment; SDQ, Strengths and Difficulties Questionnaire total; MFQ, Short Mood and Feeling Questionnaire; Std, standardised.

Discussion

This study investigated trajectories of paranoia, hallucinations, and negative symptoms from mid-adolescence to emerging adulthood in a community sample. We found evidence to suggest that the developmental distribution of these PENS dimensions was best described by multiple latent classes. Across PENS, trajectory classes identified through GMM were largely distinguished by different scores at age 16, but also by different rates of change over time. We found support for the hypothesis that persistence is associated with less favourable scores on both polygenic and phenotypic behavioural and educational attainment measures.

Of the models selected as providing the best representation of the data, a high and persisting latent trajectory class was identified only for paranoia. The percentage of individuals most likely to be assigned to this class for paranoia (~23%) mirrors the ~20% persistence rate estimated through meta-analysis of aggregated psychotic experiences reported across the lifespan, from studies that manually classified individuals (Linscott & van Os, Reference Linscott and van Os2013). Notably, the rate of persistence found in the current study is higher than in previous studies that have estimated trajectories of aggregated psychotic experiences through latent variable modelling, both in adolescence (1–16% for persistent/increasing scores) (Bourque et al., Reference Bourque, Afzali, O'Leary-Barrett and Conrod2017; Lin et al., Reference Lin, Wigman, Nelson, Vollebergh, van Os, Baksheev and Yung2011; Mackie et al., Reference Mackie, Castellanos-Ryan and Conrod2011; Thapar et al., Reference Thapar, Heron, Jones, Owen, Lewis and Zammit2012; Wigman et al., Reference Wigman, Vollebergh, Raaijmakers, Iedema, van Dorsselaer, Ormel and van Os2011c), and adulthood (12%) (Wigman et al., Reference Wigman, van Winkel, Jacobs, Wichers, Derom, Thiery and van Os2011a).

Speculatively, these latent variable modelling estimates of persistence may be attenuated in comparison to our paranoia estimate because they include information on hallucinations as well as paranoia/delusions. That is, whilst it was hypothesised that a persistent class would emerge for hallucinations as well as paranoia, in our study, the 2-class model (that did not include a high-persistent class) was selected; suggesting that the data is best represented by a decreasing developmental course across each trajectory class. Considering the empirically driven constraints on the within-class residual variances in this model, and because our study is the first to estimate latent trajectories of hallucinations as a separate dimension in the community, future research should test whether a high-persistent class for paranoia but not hallucinations is replicated in other community samples of young people using other measures.

We also hypothesised that a persistent trajectory would be identified for negative symptoms. In our study, the 2-class model was selected (because of the better relative model fit and entropy), which suggested an overall pattern of increase from adolescence to emerging adulthood. Our results are the first to estimate latent growth in negative symptoms in the community, though they may be considered in-line with findings from a sample of individuals meeting the criteria for a first episode of psychosis, in which most individuals were classified into subgroups characterised by increasing or stable symptoms (Austin et al., Reference Austin, Mors, Budtz-Jørgensen, Secher, Hjorthøj, Bertelsen and Nordentoft2015). The observed association between male sex and the most elevated negative symptoms trajectory class adds to findings that have reported cross-sectional associations between negative symptoms and male sex (Dominguez et al., Reference Dominguez, Saka, Lieb, Wittchen and van Os2010; Maric et al., Reference Maric, Krabbendam, Vollebergh, de Graaf and van Os2003; Ronald et al., Reference Ronald, Sieradzka, Cardno, Haworth, McGuire and Freeman2014).

The findings of higher GPS for clinical help seeking, major depressive disorder, and attention deficit hyperactivity disorder being associated with membership in the most elevated trajectory class across PENS, suggest that the development of PENS dimensions is at least in part associated with measured genetic variants associated with broad clinical outcomes. The association between polygenic liability for clinical help-seeking (for nerves, anxiety, tension, or depression) and for major depressive disorder, can be considered in-line with both theory and empirical findings suggesting that affective symptoms exacerbate psychotic symptoms in general, and specifically, contribute to the persistence of paranoia (Bird, Waite, Rowsell, Fergusson, & Freeman, Reference Bird, Waite, Rowsell, Fergusson and Freeman2017; Fowler et al., Reference Fowler, Hodgekins, Garety, Freeman, Kuipers, Dunn and Bebbington2012; Freeman & Garety, Reference Freeman and Garety2003; Freeman, Garety, Kuipers, Fowler, & Bebbington, Reference Freeman, Garety, Kuipers, Fowler and Bebbington2002; Reference Freeman, Stahl, McManus, Meltzer, Brugha, Wiles and Bebbington2012).

Our findings of a null association for schizophrenia GPS concur with those from a recent study of aggregated psychotic experiences measured across adolescence and emerging adulthood (Rammos et al., Reference Rammos, Sullivan, Kounali, Jones, Hammerton, Hines and Zammit2021). To the extent that schizophrenia GPS is associated with PENS measured at single time-points (Jones et al., Reference Jones, Stergiakouli, Tansey, Hubbard, Heron, Cannon and Zammit2016; Pain et al., Reference Pain, Dudbridge, Cardno, Freeman, Lu, Lundstrom and Ronald2018), the current results suggest that polygenic liability for schizophrenia may influence the static expression but not the development of PENS. A study that investigated the effect of polygenic liability to schizophrenia GPS on the developmental course of separate dimensions within the negative symptoms construct in schizophrenia found that schizophrenia GPS predicted a more severe course specifically of avolition (Jonas et al., Reference Jonas, Lencz, Li, Malhotra, Perlman, Fochtmann and Kotov2019). This aligns broadly with suggestive cross-sectional GPS findings in the community (Havers, Cardno, Freeman, & Ronald, Reference Havers, Cardno, Freeman and Ronald2022), and with clinical findings and theory suggesting that avolition may be a central feature of the negative symptoms construct (Foussias & Remington, Reference Foussias and Remington2010; Strauss et al., Reference Strauss, Zamani Esfahlani, Sayama, Kirkpatrick, Opler, Saoud and Luthringer2020, Reference Strauss, Bartolomeo and Harvey2021).

Our findings further suggest that family history of psychosis and bipolar disorder, which were associated with the most elevated course of paranoia and hallucinations (though not negative symptoms), is not due to an increased burden of measured polygenic variants for schizophrenia, echoing the findings of Rammos et al. (Reference Rammos, Sullivan, Kounali, Jones, Hammerton, Hines and Zammit2021). A range of evidence suggests that psychotic experiences likely reflect a broad transdiagnostic risk, rather than genetic risk for clinical psychosis, specifically (McGorry & Mei, Reference McGorry and Mei2021; McGrath et al., Reference McGrath, Saha, Al-Hamzawi, Andrade, Benjet, Bromet and Kessler2016).

An unexpected result was that higher years of education and higher intelligence GPS, as well as higher SES and phenotypic educational attainment at ages 7 and 22, were associated with membership in the high-persistent class for paranoia. These findings are hard to interpret in a theoretical context. Whilst broadly similar associations have been found between these measures and paranoia reported at a single time-point in adulthood (Freeman et al., Reference Freeman, McManus, Brugha, Meltzer, Jenkins and Bebbington2011), our findings may further reflect dimension-specific polygenic and behavioural associations specifically for the development of paranoia, which may otherwise be obscured when aggregated measures are used. Replication in other samples using dimensional PENS scales is required to explore this suggestion further and could further test the extent to which the enforced model constraints, and the observed classification error (reflected in the entropy), may have influenced the unexpected results.

Several strengths and limitations of our study should be highlighted. Key strengths are the estimation of trajectories for separate PENS dimensions using data from a community sample, inherently free of treatment confounds, across a period that reflects when psychotic disorders are likely to first emerge (Kessler et al., Reference Kessler, Amminger, Aguilar-Gaxiola, Alonso, Lee and Ustun2007; Maibing et al., Reference Maibing, Pedersen, Benros, Mortensen, Dalsgaard and Nordentoft2015). Further, utilising full information maximum likelihood allowed for the estimation of trajectories using data from all individuals. Nonetheless, only families already responding at 16 were invited to participate at age 17, and the sample at this age was smaller than at ages 16 and 22. A greater number of repeated measures collected over a greater time-period would allow for a more expansive investigation into nonlinear aspects of growth and provide a broader view of the development of PENS. Further, it is highlighted that whilst the multiple latent classes could be representative of underlying subgroups of individuals, they should primarily be considered as statistical approximations. A final limitation is that with newer, larger GWASs, the polygenic score results may change, and this is an important area for future research.

In summary, we modelled latent heterogeneity in the development of paranoia, hallucinations, and negative symptoms in the community. Studying specific PENS dimensions allowed for distinct patterns of growth to be estimated. Our results suggest largely dimension-wide but also dimension-specific polygenic and behavioural associations with the developmental trajectories of PENS. These findings add to a growing body of literature that suggests that a dimension-specific approach may be important for delineating aetiological and developmental pathways for PENS, which in turn may reduce poor outcomes by facilitating more precise intervention and prevention efforts.

Supplementary material

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

Acknowledgements

The authors thank the TEDS participants, and Robert Plomin and Andrew McMillan for the collaboration. LH thanks Eric Thibodeau for providing valuable comments on the statistical modelling used in this study.

Financial support

This work was supported by a UK Medical Research Council grant to AR (G1100559). TEDS is supported by a programme grant to Professor Robert Plomin from the UK Medical Research Council (MR/V012878/1 and previously MR/M021475/1), with additional support from the US National Institutes of Health (AG046938). LH is supported by a UK Economics and Social Research Council studentship. DF is a National Institute for Health Research Senior Investigator. SvS is a recipient of a British Academy Mid-Career Fellowship (MF21\210017) and of a CRISP Jacobs Fellowship.

Conflict of interest

None.

References

Andreasen, N. C. (1982). Negative symptoms in schizophrenia: Definition and reliability. Archives of General Psychiatry, 39(7), 784–788. https://doi.org/10.1001/archpsyc.1982.04290070020005.CrossRefGoogle ScholarPubMed
Angold, A., Costello, E. J., Messer, S. C., & Pickles, A. (1995). Development of a short questionnaire for use in epidemiological studies of depression in children and adolescents. International Journal of Methods in Psychiatric Research, 5(4), 237–249.Google Scholar
Armando, M., Nelson, B., Yung, A. R., Ross, M., Birchwood, M., Girardi, P., & Nastro, P. F. (2010). Psychotic-like experiences and correlation with distress and depressive symptoms in a community sample of adolescents and young adults. Schizophrenia Research, 119(1), 258–265. https://doi.org/10.1016/j.schres.2010.03.001.CrossRefGoogle Scholar
Asparouhov, T., & Muthén, B. (2014). Auxiliary variables in mixture modeling: Three-step approaches using mplus. Structural Equation Modeling: A Multidisciplinary Journal, 21(3), 329–341. https://doi.org/10.1080/10705511.2014.915181.CrossRefGoogle Scholar
Austin, S. F., Mors, O., Budtz-Jørgensen, E., Secher, R. G., Hjorthøj, C. R., Bertelsen, M., … Nordentoft, M. (2015). Long-term trajectories of positive and negative symptoms in first episode psychosis: A 10 year follow-up study in the OPUS cohort. Schizophrenia Research, 168(1), 84–91. https://doi.org/10.1016/j.schres.2015.07.021.CrossRefGoogle Scholar
Bakk, Z., & Vermunt, J. K. (2016). Robustness of stepwise latent class modeling with continuous distal outcomes. Structural Equation Modeling: A Multidisciplinary Journal, 23(1), 20–31. https://doi.org/10.1080/10705511.2014.955104.CrossRefGoogle Scholar
Barragan, M., Laurens, K. R., Navarro, J. B., & Obiols, J. E. (2011). Psychotic-like experiences and depressive symptoms in a community sample of adolescents. European Psychiatry, 26(6), 396–401. https://doi.org/10.1016/j.eurpsy.2010.12.007.CrossRefGoogle Scholar
Bartels-Velthuis, A. A., van de Willige, G., Jenner, J. A., van Os, J., & Wiersma, D. (2011). Course of auditory vocal hallucinations in childhood: 5-year follow-up study. The British Journal of Psychiatry, 199(4), 296–302. https://doi.org/10.1192/bjp.bp.110.086918.CrossRefGoogle ScholarPubMed
Bell, V., Halligan, P. W., & Ellis, H. D. (2006). The Cardiff Anomalous Perceptions Scale (CAPS): A new validated measure of anomalous perceptual experience. Schizophrenia Bulletin, 32(2), 366–377. https://doi.org/10.1093/schbul/sbj014.CrossRefGoogle ScholarPubMed
Benjamini, Y., & Hochberg, Y. (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. Journal of the Royal Statistical Society: Series B (Methodological), 57(1), 289–300. https://doi.org/10.1111/j.2517-6161.1995.tb02031.x.Google Scholar
Bird, J. C., Waite, F., Rowsell, E., Fergusson, E. C., & Freeman, D. (2017). Cognitive, affective, and social factors maintaining paranoia in adolescents with mental health problems: A longitudinal study. Psychiatry Research, 257, 34–39. https://doi.org/10.1016/j.psychres.2017.07.023.CrossRefGoogle ScholarPubMed
Bolck, A., Croon, M., & Hagenaars, J. (2004). Estimating latent structure models with categorical variables: One-step versus three-step estimators. Political Analysis, 12(1), 3–27. https://doi.org/10.1093/pan/mph001.CrossRefGoogle Scholar
Bourque, J., Afzali, M. H., O'Leary-Barrett, M., & Conrod, P. (2017). Cannabis use and psychotic-like experiences trajectories during early adolescence: The coevolution and potential mediators. Journal of Child Psychology and Psychiatry, 58(12), 1360–1369. https://doi.org/10.1111/jcpp.12765.CrossRefGoogle ScholarPubMed
Connell, M., Betts, K., McGrath, J. J., Alati, R., Najman, J., Clavarino, A., … Scott, J. G. (2016). Hallucinations in adolescents and risk for mental disorders and suicidal behaviour in adulthood: Prospective evidence from the MUSP birth cohort study. Schizophrenia Research, 176(2), 546–551. https://doi.org/10.1016/j.schres.2016.06.009.CrossRefGoogle ScholarPubMed
Cosgrave, J., Purple, R. J., Haines, R., Porcheret, K., van Heugten-van der Kloet, D., Johns, L., … Wulff, K. (2021). Do environmental risk factors for the development of psychosis distribute differently across dimensionally assessed psychotic experiences? Translational Psychiatry, 11(1), 1–13. https://doi.org/10.1038/s41398-021-01265-2.CrossRefGoogle ScholarPubMed
Cougnard, A., Marcelis, M., Myin-Germeys, I., Graaf, R. D., Vollebergh, W., Krabbendam, L., … van Os, J. (2007). Does normal developmental expression of psychosis combine with environmental risk to cause persistence of psychosis? A psychosis proneness–persistence model. Psychological Medicine, 37(4), 513–527. https://doi.org/10.1017/S0033291706009731.CrossRefGoogle ScholarPubMed
De Loore, E., Gunther, N., Drukker, M., Feron, F., Sabbe, B., Deboutte, D., … Myin-Germeys, I. (2011). Persistence and outcome of auditory hallucinations in adolescence: A longitudinal general population study of 1800 individuals. Schizophrenia Research, 127(1), 252–256. https://doi.org/10.1016/j.schres.2011.01.015.CrossRefGoogle ScholarPubMed
DeVylder, J. E., Lehmann, M., & Chen, F. (2015). Social and clinical correlates of the persistence of psychotic experiences in the general population. Schizophrenia Research, 169(1), 286–291. https://doi.org/10.1016/j.schres.2015.08.039.CrossRefGoogle ScholarPubMed
Dhossche, D., Ferdinand, R., Ende, J. V. D., Hofstra, M. B., & Verhulst, F. (2002). Diagnostic outcome of self-reported hallucinations in a community sample of adolescents. Psychological Medicine, 32(4), 619–627. https://doi.org/10.1017/S003329170200555X.CrossRefGoogle Scholar
Dominguez, M.-G., Saka, M. C., Lieb, R., Wittchen, H.-U., & van Os, J. (2010). Early expression of negative/disorganized symptoms predicting psychotic experiences and subsequent clinical psychosis: A 10-year study. American Journal of Psychiatry, 167(9), 1075–1082. https://doi.org/10.1176/appi.ajp.2010.09060883.CrossRefGoogle ScholarPubMed
Dominguez, M.-G., Wichers, M., Lieb, R., Wittchen, H.-U., & van Os, J. (2011). Evidence that onset of clinical psychosis is an outcome of progressively more persistent subclinical psychotic experiences: An 8-year cohort study. Schizophrenia Bulletin, 37(1), 84–93. https://doi.org/10.1093/schbul/sbp022.CrossRefGoogle ScholarPubMed
Ericson, M., Tuvblad, C., Raine, A., Young-Wolff, K., & Baker, L. A. (2011). Heritability and longitudinal stability of schizotypal traits during adolescence. Behavior Genetics, 41(4), 499–511. https://doi.org/10.1007/s10519-010-9401-x.CrossRefGoogle ScholarPubMed
Foussias, G., & Remington, G. (2010). Negative symptoms in schizophrenia: Avolition and Occam's razor. Schizophrenia Bulletin, 36(2), 359–369. https://doi.org/10.1093/schbul/sbn094.CrossRefGoogle ScholarPubMed
Fowler, D., Hodgekins, J., Garety, P., Freeman, D., Kuipers, E., Dunn, G., … Bebbington, P. E. (2012). Negative cognition, depressed mood, and paranoia: A longitudinal pathway analysis using structural equation modeling. Schizophrenia Bulletin, 38(5), 1063–1073. https://doi.org/10.1093/schbul/sbr019.CrossRefGoogle ScholarPubMed
Freeman, D., & Garety, P. A. (2003). Connecting neurosis and psychosis: The direct influence of emotion on delusions and hallucinations. Behaviour Research and Therapy, 41(8), 923–947. https://doi.org/10.1016/s0005-7967(02)00104-3.CrossRefGoogle ScholarPubMed
Freeman, D., Garety, P. A., Bebbington, P. E., Smith, B., Rollinson, R., Fowler, D., … Dunn, G. (2005). Psychological investigation of the structure of paranoia in a non-clinical population. The British Journal of Psychiatry, 186(5), 427–435. https://doi.org/10.1192/bjp.186.5.427.CrossRefGoogle Scholar
Freeman, D., Garety, P. A., Kuipers, E., Fowler, D., & Bebbington, P. E. (2002). A cognitive model of persecutory delusions. British Journal of Clinical Psychology, 41(4), 331–347. https://doi.org/10.1348/014466502760387461.CrossRefGoogle ScholarPubMed
Freeman, D., McManus, S., Brugha, T., Meltzer, H., Jenkins, R., & Bebbington, P. (2011). Concomitants of paranoia in the general population. Psychological Medicine, 41(5), 923–936. https://doi.org/10.1017/S0033291710001546.CrossRefGoogle ScholarPubMed
Freeman, D., Stahl, D., McManus, S., Meltzer, H., Brugha, T., Wiles, N., & Bebbington, P. (2012). Insomnia, worry, anxiety and depression as predictors of the occurrence and persistence of paranoid thinking. Social Psychiatry and Psychiatric Epidemiology, 47(8), 1195–1203. https://doi.org/10.1007/s00127-011-0433-1.CrossRefGoogle ScholarPubMed
Goodman, R. (1997). The strengths and difficulties questionnaire: A research note. Journal of Child Psychology and Psychiatry, 38(5), 581–586. https://doi.org/10.1111/j.1469-7610.1997.tb01545.x.CrossRefGoogle ScholarPubMed
Hafeez, D., & Yung, A. R. (2021). Early persistence of psychotic-like experiences in a community sample of adolescents. Early Intervention in Psychiatry, 15(6), 1575–1583. https://doi.org/10.1111/eip.13096.CrossRefGoogle Scholar
Havers, L., Cardno, A., Freeman, D., & Ronald, A. (2022). The latent structure of negative symptoms in the general population in adolescence and emerging adulthood. Schizophrenia Bulletin Open, 3(1), sgac009. https://doi.org/10.1093/schizbullopen/sgac009.CrossRefGoogle ScholarPubMed
Havers, L., Taylor, M. J., & Ronald, A. (2019). Genetic and environmental influences on the stability of psychotic experiences and negative symptoms in adolescence. Journal of Child Psychology and Psychiatry, 60(7), 784–792. https://doi.org/10.1111/jcpp.13045.CrossRefGoogle ScholarPubMed
Healy, C., Brannigan, R., Dooley, N., Coughlan, H., Clarke, M., Kelleher, I., & Cannon, M. (2019). Childhood and adolescent psychotic experiences and risk of mental disorder: A systematic review and meta-analysis. Psychological Medicine, 49(10), 1589–1599. https://doi.org/10.1017/S0033291719000485.CrossRefGoogle ScholarPubMed
Herle, M., Micali, N., Abdulkadir, M., Loos, R., Bryant-Waugh, R., Hübel, C., … De Stavola, B. L. (2020). Identifying typical trajectories in longitudinal data: Modelling strategies and interpretations. European Journal of Epidemiology, 35(3), 205–222. https://doi.org/10.1007/s10654-020-00615-6.CrossRefGoogle ScholarPubMed
Hielscher, E., DeVylder, J., Hasking, P., Connell, M., Martin, G., & Scott, J. G. (2021). Can't get you out of my head: Persistence and remission of psychotic experiences in adolescents and its association with self-injury and suicide attempts. Schizophrenia Research, 229, 63–72. https://doi.org/10.1016/j.schres.2020.11.019.CrossRefGoogle ScholarPubMed
Janssens, M., Boyette, L.-L., Heering, H. D., Bartels-Velthuis, A. A., Lataster, T., Kahn, R. S., … Myin-Germeys, I. (2016). Developmental course of subclinical positive and negative psychotic symptoms and their associations with genetic risk status and impairment. Schizophrenia Research, 174(1), 177–182. https://doi.org/10.1016/j.schres.2016.03.028.CrossRefGoogle ScholarPubMed
Jonas, K. G., Lencz, T., Li, K., Malhotra, A. K., Perlman, G., Fochtmann, L. J., … Kotov, R. (2019). Schizophrenia polygenic risk score and 20-year course of illness in psychotic disorders. Translational Psychiatry, 9(1), 1–8. https://doi.org/10.1038/s41398-019-0612-5.CrossRefGoogle ScholarPubMed
Jones, H. J., Stergiakouli, E., Tansey, K. E., Hubbard, L., Heron, J., Cannon, M., … Zammit, S. (2016). Phenotypic manifestation of genetic risk for schizophrenia during adolescence in the general population. JAMA Psychiatry, 73(3), 221–228. https://doi.org/10.1001/jamapsychiatry.2015.3058.CrossRefGoogle ScholarPubMed
Jung, T., & Wickrama, K. A. (2008). An introduction to latent class growth analysis and growth mixture modeling. Social and Personality Psychology Compass, 2(1), 302–317.CrossRefGoogle Scholar
Kalman, J. L., Bresnahan, M., Schulze, T. G., & Susser, E. (2019). Predictors of persisting psychotic-like experiences in children and adolescents: A scoping review. Schizophrenia Research, 209, 32–39. https://doi.org/10.1016/j.schres.2019.05.012.CrossRefGoogle ScholarPubMed
Kaymaz, N., Drukker, M., Lieb, R., Wittchen, H.-U., Werbeloff, N., Weiser, M., … van Os, J. (2012). Do subthreshold psychotic experiences predict clinical outcomes in unselected non-help-seeking population-based samples? A systematic review and meta-analysis, enriched with new results. Psychological Medicine, 42(11), 2239–2253. https://doi.org/10.1017/S0033291711002911.CrossRefGoogle ScholarPubMed
Kelleher, I., Connor, D., Clarke, M. C., Devlin, N., Harley, M., & Cannon, M. (2012). Prevalence of psychotic symptoms in childhood and adolescence: A systematic review and meta-analysis of population-based studies. Psychological Medicine, 42(9), 1857–1863. https://doi.org/10.1017/S0033291711002960.CrossRefGoogle ScholarPubMed
Kessler, R. C., Amminger, G. P., Aguilar-Gaxiola, S., Alonso, J., Lee, S., & Ustun, T. B. (2007). Age of onset of mental disorders: A review of recent literature. Current Opinion in Psychiatry, 20(4), 359–364. https://doi.org/10.1097/YCO.0b013e32816ebc8c.CrossRefGoogle ScholarPubMed
Kim-Cohen, J., Caspi, A., Moffitt, T. E., Harrington, H., Milne, B. J., & Poulton, R. (2003). Prior juvenile diagnoses in adults with mental disorder: Developmental follow-back of a prospective-longitudinal cohort. Archives of General Psychiatry, 60(7), 709–717. https://doi.org/10.1001/archpsyc.60.7.709.CrossRefGoogle ScholarPubMed
Lin, A., Wigman, J. T. W., Nelson, B., Vollebergh, W. A. M., van Os, J., Baksheev, G., … Yung, A. R. (2011). The relationship between coping and subclinical psychotic experiences in adolescents from the general population – a longitudinal study. Psychological Medicine, 41(12), 2535–2546. https://doi.org/10.1017/S0033291711000560.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
Mackie, C. J., Castellanos-Ryan, N., & Conrod, P. J. (2011). Developmental trajectories of psychotic-like experiences across adolescence: Impact of victimization and substance use. Psychological Medicine, 41(1), 47–58. https://doi.org/10.1017/S0033291710000449.CrossRefGoogle ScholarPubMed
Mackie, C. J., O'Leary-Barrett, M., Al-Khudhairy, N., Castellanos-Ryan, N., Struve, M., Topper, L., & Conrod, P. (2013). Adolescent bullying, cannabis use and emerging psychotic experiences: A longitudinal general population study. Psychological Medicine, 43(5), 1033–1044. https://doi.org/10.1017/S003329171200205X.CrossRefGoogle ScholarPubMed
Maibing, C. F., Pedersen, C. B., Benros, M. E., Mortensen, P. B., Dalsgaard, S., & Nordentoft, M. (2015). Risk of schizophrenia increases after all child and adolescent psychiatric disorders: A nationwide study. Schizophrenia Bulletin, 41(4), 963–970. https://doi.org/10.1093/schbul/sbu119.CrossRefGoogle ScholarPubMed
Maric, N., Krabbendam, L., Vollebergh, W., de Graaf, R., & van Os, J. (2003). Sex differences in symptoms of psychosis in a non-selected, general population sample. Schizophrenia Research, 63(1), 89–95. https://doi.org/10.1016/S0920-9964(02)00380-8.CrossRefGoogle Scholar
McGorry, P. D., & Mei, C. (2021). Clinical staging for youth mental disorders: Progress in reforming diagnosis and clinical care. Annual Review of Developmental Psychology, 3, 15–39.CrossRefGoogle Scholar
McGrath, J. J., Saha, S., Al-Hamzawi, A., Alonso, J., Bromet, E. J., Bruffaerts, R., … Kessler, R. C. (2015). Psychotic experiences in the general population: A cross-national analysis based on 31 261 respondents from 18 countries. JAMA Psychiatry, 72(7), 697–705. https://doi.org/10.1001/jamapsychiatry.2015.0575.CrossRefGoogle ScholarPubMed
McGrath, J. J., Saha, S., Al-Hamzawi, A., Andrade, L., Benjet, C., Bromet, E. J., … Kessler, R. C. (2016). The bidirectional associations between psychotic experiences and DSM-IV mental disorders. American Journal of Psychiatry, 173(10), 997–1006. https://doi.org/10.1176/appi.ajp.2016.15101293.CrossRefGoogle ScholarPubMed
Mehta, P. D., & West, S. G. (2000). Putting the individual back into individual growth curves. Psychological Methods, 5(1), 23–43.CrossRefGoogle ScholarPubMed
Mitchell, C. A. A., Maybery, M. T., Russell-Smith, S. N., Collerton, D., Gignac, G. E., & Waters, F. (2017). The structure and measurement of unusual sensory experiences in different modalities: The Multi-Modality Unusual Sensory Experiences Questionnaire (MUSEQ). Frontiers in Psychology, 8, 1363. https://www.frontiersin.org/article/10.3389/fpsyg.2017.01363.CrossRefGoogle ScholarPubMed
Pain, O., Dudbridge, F., Cardno, A. G., Freeman, D., Lu, Y., Lundstrom, S., … Ronald, A. (2018). Genome-wide analysis of adolescent psychotic-like experiences shows genetic overlap with psychiatric disorders. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics, 177(4), 416–425. https://doi.org/10.1002/ajmg.b.32630.CrossRefGoogle ScholarPubMed
Ram, N., & Grimm, K. J. (2009). Growth mixture modeling: A method for identifying differences in longitudinal change among unobserved groups. International Journal of Behavioral Development, 33(6), 565–576. https://doi.org/10.1177/0165025409343765.CrossRefGoogle ScholarPubMed
Rammos, A., Sullivan, S. A., Kounali, D., Jones, H. J., Hammerton, G., Hines, L. A., … Zammit, S. (2021). Precursors and correlates of transient and persistent longitudinal profiles of psychotic experiences from late childhood through early adulthood. The British Journal of Psychiatry, 220(6), 1–9. https://doi.org/10.1192/bjp.2021.145.Google ScholarPubMed
Rimfeld, K., Malanchini, M., Spargo, T., Spickernell, G., Selzam, S., McMillan, A., … Plomin, R. (2019). Twins early development study: A genetically sensitive investigation into behavioral and cognitive development from infancy to emerging adulthood. Twin Research and Human Genetics, 22(6), 508–513. https://doi.org/10.1017/thg.2019.56.CrossRefGoogle ScholarPubMed
Ronald, A., & Pain, O. (2018). A systematic review of genome-wide research on psychotic experiences and negative symptom traits: New revelations and implications for psychiatry. Human Molecular Genetics, 27(R2), R136–R152. https://doi.org/10.1093/hmg/ddy157.Google ScholarPubMed
Ronald, A., Sieradzka, D., Cardno, A. G., Haworth, C. M. A., McGuire, P., & Freeman, D. (2014). Characterization of psychotic experiences in adolescence using the specific psychotic experiences questionnaire: Findings from a study of 5000 16-year-old twins. Schizophrenia Bulletin, 40(4), 868–877. https://doi.org/10.1093/schbul/sbt106.CrossRefGoogle ScholarPubMed
Roy, M.-A., Maziade, M., Labbé, A., & Mérette, C. (2001). Male gender is associated with deficit schizophrenia: A meta-analysis. Schizophrenia Research, 47(2), 141–147. https://doi.org/10.1016/S0920-9964(99)00231-5.CrossRefGoogle ScholarPubMed
Selzam, S., McAdams, T. A., Coleman, J. R. I., Carnell, S., O'Reilly, P. F., Plomin, R., & Llewellyn, C. H. (2018). Evidence for gene-environment correlation in child feeding: Links between common genetic variation for BMI in children and parental feeding practices. PLoS Genetics, 14(11). https://doi.org/10.1371/journal.pgen.1007757.CrossRefGoogle ScholarPubMed
Selzam, S., Ritchie, S. J., Pingault, J.-B., Reynolds, C. A., O'Reilly, P. F., & Plomin, R. (2019). Comparing within- and between-family polygenic score prediction. American Journal of Human Genetics, 105(2), 351–363.CrossRefGoogle ScholarPubMed
Shakoor, S., McGuire, P., Cardno, A. G., Freeman, D., Plomin, R., & Ronald, A. (2015). A shared genetic propensity underlies experiences of bullying victimization in late childhood and self-rated paranoid thinking in adolescence. Schizophrenia Bulletin, 41(3), 754–763. https://doi.org/10.1093/schbul/sbu142.CrossRefGoogle ScholarPubMed
Sheaves, B., Bebbington, P. E., Goodwin, G. M., Harrison, P. J., Espie, C. A., Foster, R. G., & Freeman, D. (2016). Insomnia and hallucinations in the general population: Findings from the 2000 and 2007 British psychiatric morbidity surveys. Psychiatry Research, 241, 141–146. https://doi.org/10.1016/j.psychres.2016.03.055.CrossRefGoogle ScholarPubMed
Shireman, E. M., Steinley, D., & Brusco, M. J. (2016). Local optima in mixture modeling. Multivariate Behavioral Research, 51(4), 466–481. https://doi.org/10.1080/00273171.2016.1160359.CrossRefGoogle ScholarPubMed
Sijbrandij, J. J., Hoekstra, T., Almansa, J., Peeters, M., Bültmann, U., & Reijneveld, S. A. (2020). Variance constraints strongly influenced model performance in growth mixture modeling: A simulation and empirical study. BMC Medical Research Methodology, 20(1), 1–15.CrossRefGoogle ScholarPubMed
Steenkamp, L. R., Tiemeier, H., Blanken, L. M. E., Hillegers, M. H. J., Kushner, S. A., & Bolhuis, K. (2021). Predicting persistence of hallucinations from childhood to adolescence. The British Journal of Psychiatry, 219(6), 670–677. https://doi.org/10.1192/bjp.2021.115.CrossRefGoogle ScholarPubMed
Stefanis, N. C., Delespaul, P., Henquet, C., Bakoula, C., Stefanis, C. N., & van Os, J. (2004). Early adolescent cannabis exposure and positive and negative dimensions of psychosis. Addiction, 99(10), 1333–1341. https://doi.org/10.1111/j.1360-0443.2004.00806.x.CrossRefGoogle ScholarPubMed
Sterba, S. K. (2014). Fitting nonlinear latent growth curve models with individually varying time points. Structural Equation Modeling: A Multidisciplinary Journal, 21(4), 630–647. https://doi.org/10.1080/10705511.2014.919828.CrossRefGoogle Scholar
Strauss, G. P., Bartolomeo, L. A., & Harvey, P. D. (2021). Avolition as the core negative symptom in schizophrenia: Relevance to pharmacological treatment development. Npj Schizophrenia, 7(1), 1–6. https://doi.org/10.1038/s41537-021-00145-4.CrossRefGoogle ScholarPubMed
Strauss, G. P., Zamani Esfahlani, F., Sayama, H., Kirkpatrick, B., Opler, M. G., Saoud, J. B., … Luthringer, R. (2020). Network analysis indicates that avolition is the most central domain for the successful treatment of negative symptoms: Evidence from the roluperidone randomized clinical trial. Schizophrenia Bulletin, 46(4), 964–970. https://doi.org/10.1093/schbul/sbz141.CrossRefGoogle ScholarPubMed
Thapar, A., Heron, J., Jones, R. B., Owen, M. J., Lewis, G., & Zammit, S. (2012). Trajectories of change in self-reported psychotic-like experiences in childhood and adolescence. Schizophrenia Research, 140(1), 104–109. https://doi.org/10.1016/j.schres.2012.06.024.CrossRefGoogle ScholarPubMed
van de Schoot, R., Sijbrandij, M., Winter, S. D., Depaoli, S., & Vermunt, J. K. (2017). The GRoLTS-checklist: Guidelines for reporting on latent trajectory studies. Structural Equation Modeling: A Multidisciplinary Journal, 24(3), 451–467.CrossRefGoogle Scholar
van Os, J., Linscott, R. J., Myin-Germeys, I., Delespaul, P., & Krabbendam, L. (2009). A systematic review and meta-analysis of the psychosis continuum: Evidence for a psychosis proneness–persistence–impairment model of psychotic disorder. Psychological Medicine, 39(2), 179–195. https://doi.org/10.1017/S0033291708003814.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
van Rossum, I., Dominguez, M.-G., Lieb, R., Wittchen, H.-U., & van Os, J. (2011). Affective dysregulation and reality distortion: A 10-year prospective study of their association and clinical relevance. Schizophrenia Bulletin, 37(3), 561–571. https://doi.org/10.1093/schbul/sbp101.CrossRefGoogle ScholarPubMed
Vermunt, J. K. (2017). Latent class modeling with covariates: Two improved three-step approaches. Political Analysis, 18(4), 450–469. https://doi.org/10.1093/pan/mpq025.CrossRefGoogle Scholar
Wigman, J. T. W., van Winkel, R., Jacobs, N., Wichers, M., Derom, C., Thiery, E., … van Os, J. (2011a). A twin study of genetic and environmental determinants of abnormal persistence of psychotic experiences in young adulthood. American Journal of Medical Genetics Part B: Neuropsychiatric Genetics, 156(5), 546–552. https://doi.org/10.1002/ajmg.b.31193.CrossRefGoogle Scholar
Wigman, J. T. W., van Winkel, R., Raaijmakers, Q. A. W., Ormel, J., Verhulst, F. C., Reijneveld, S. A., … Vollebergh, W. A. M. (2011b). Evidence for a persistent, environment-dependent and deteriorating subtype of subclinical psychotic experiences: A 6-year longitudinal general population study. Psychological Medicine, 41(11), 2317–2329. https://doi.org/10.1017/S0033291711000304.CrossRefGoogle ScholarPubMed
Wigman, J. T. W., Vollebergh, W. A. M., Raaijmakers, Q. A. W., Iedema, J., van Dorsselaer, S., Ormel, J., … van Os, J. (2011c). The structure of the extended psychosis phenotype in early adolescence – A cross-sample replication. Schizophrenia Bulletin, 37(4), 850–860. https://doi.org/10.1093/schbul/sbp154.CrossRefGoogle ScholarPubMed
Yung, A. R., Nelson, B., Baker, K., Buckby, J. A., Baksheev, G., & Cosgrave, E. M. (2009). Psychotic-like experiences in a community sample of adolescents: Implications for the continuum model of psychosis and prediction of schizophrenia. Australian & New Zealand Journal of Psychiatry, 43(2), 118–128. https://doi.org/10.1080/00048670802607188.CrossRefGoogle Scholar
Zavos, H. M. S., Freeman, D., Haworth, C. M. A., McGuire, P., Plomin, R., Cardno, A. G., & Ronald, A. (2014). Consistent etiology of severe, frequent psychotic experiences and milder, less frequent manifestations: A twin study of specific psychotic experiences in adolescence. JAMA Psychiatry, 71(9), 1049–1057. https://doi.org/10.1001/jamapsychiatry.2014.994.CrossRefGoogle ScholarPubMed
Figure 0

Table 1. Growth mixture model fit results for converged models of paranoia, hallucinations, and negative symptoms

Figure 1

Table 2. Parameter estimates of best-fitting model for paranoia, hallucinations, and negative symptoms

Figure 2

Fig. 1. Trajectory plots of best-fitting models for paranoia, hallucinations, and negative symptoms.Note. Lines represent mean trajectories; bands represent 95% confidence intervals. Parameter estimates for the trajectories are reported in Table 2.

Figure 3

Table 3. Characteristics of latent trajectory classes for paranoia, hallucinations, and negative symptoms

Supplementary material: File

Havers et al. supplementary material

Havers et al. supplementary material

Download Havers et al. supplementary material(File)
File 4.2 MB