Despite extensive research efforts, the patho-aetiology of psychosis spectrum disorders (PSD) remains poorly understood. A significant challenge stems from the broad range of clinical presentations within PSD. Notably, sex differences appear to contribute to this heterogeneity, as evidenced by dissimilarities between men and women in terms of the onset, progression, prognosis and postulated aetiology of PSD. Reference Pence, Pries, Ferrara, Rutten, van Os and Guloksuz1,Reference Ferrara, Curtarello, Gentili, Domenicano, Vecchioni and Zese2
Sex/gender differences in mental health trajectories during puberty
Many of the differences become apparent during sensitive developmental phases such as puberty. During this period, sex-specific (biological) and gender-specific (social and/or cultural) processes occur that have a central role in the emergence of distinct symptom profiles. Boys seem to exhibit greater externalising issues, especially at the ages of 6 to 11 years. Reference Rescorla, Achenbach, Ivanova, Dumenci, Almqvist and Bilenberg3 On the other hand, pre-pubertal girls and boys exhibit similarities in internalising problems, with disparities becoming apparent during and after puberty. Particularly from the age of 12 years, girls manifest more internalising symptoms such as depression compared with boys. Reference Rescorla, Achenbach, Ivanova, Dumenci, Almqvist and Bilenberg3,Reference Morken, Viddal, Von Soest and Wichstrøm4 These differences may be at least partly be attributed to increased vulnerability in girls during this period. Reference Morken, Viddal, Von Soest and Wichstrøm4 Similarly, PSD tends to surface during adolescence and young adulthood, predominantly in men, Reference Ochoa, Usall, Cobo, Labad and Kulkarni5 and clinical presentation is often more severe in men. Reference Ferrara, Curtarello, Gentili, Domenicano, Vecchioni and Zese2 Therefore, exploration of the underlying neurobiology during maturation may lead to opportunities for intervention.
The role of puberty, early menarche and oestrogen in the development of psychosis spectrum disorders
Extensive research has established a link between early puberty and elevated externalising and internalising psychopathology in both males and females. Reference Ullsperger and Nikolas6 By contrast, there has been limited research to explore the relationship with PSD, with some studies indicating that early menarche (age at first period) may be associated with fewer negative symptoms, improved functioning and later onset in female patients with schizophrenia. Reference Hochman and Lewine7–Reference Ruiz, Blanco, Santander and Miranda11 A recent study using UK Biobank data identified 619 phenotypes associated with genetic risk for early menarche. Notably, it identified a positive association with variables related to depression and a negative association with distress stemming from psychotic experiences. Reference Magnus, Guyatt, Lawn, Wyss, Trajanoska and Küpers12 Early menarche may mark vulnerability to psychopathologies such as depression and anxiety, while potentially offering protection against PSD, particularly in the context of oestrogen and neural maturation. Reference Damme, Hernandez and Mittal13,Reference Damme, Ristanovic, Vargas and Mittal14 The ‘oestrogen hypothesis’ suggests that oestrogens, especially oestradiol, protect females against PSD. Levels of oestradiol, one of three oestrogens, increase during female maturation. Reference Frederiksen, Johannsen, Andersen, Albrethsen, Landersoe and Petersen15 This may modify mechanisms that are potentially linked to PSD, including mitochondrial function, dopamine activity and stress-related systems. Reference Gogos, Sbisa, Sun, Gibbons, Udawela and Dean16 Moreover, menstrual phases with higher oestradiol levels are linked to reduced symptoms of PSD, Reference Handy, Greenfield, Yonkers and Payne17 whereas menopause (with decreased oestradiol levels) is linked to increased incidence of PSD in females. Reference Culbert, Thakkar and Klump18 Preliminary research suggests that oestradiol treatment may improve symptoms in male and female PSD patients. Reference Gogos, Sbisa, Sun, Gibbons, Udawela and Dean16 However, limited results are available on the effects of oestradiol and puberty timing on psychotic experiences in children and adolescents.
The role of genomic and exposomic factors in the development of psychosis spectrum disorders
PSD has a complex aetiology, shaped by the genome and exposome (the environmental exposures an individual encounters over their lifetime). Reference Pries, Dal Ferro, van Os, Delespaul, Kenis and Lin19 Multiple genes constitute genomic liability, and various environmental factors, including childhood adversity and cannabis use, contribute to exposomic liability. Reference Pries, Dal Ferro, van Os, Delespaul, Kenis and Lin19,Reference Pries, Moore, Visoki, Sotelo, Barzilay and Guloksuz20 We have previously shown in a case–control study that the interplay between genomic and exposomic liabilities significantly amplifies the risk of developing schizophrenia, Reference Pries, Dal Ferro, van Os, Delespaul, Kenis and Lin19 as well as demonstrating that it increased the risk of experiencing significantly distressing psychotic experiences among child and adolescent participants in the Adolescent Brain Cognitive Development (ABCD) Study. Reference Di Vincenzo, Prachason, Sampogna, Arias-Magnasco, Lin and Pries21 Another ABCD cohort study found that early puberty, combined with heightened environmental risk represented by a latent factor, was linked to increased rule-breaking behaviours in both males and females, as well as greater depressive symptoms in females. Reference Vijayakumar, Youssef, Bereznicki, Dehestani, Silk and Whittle22 These findings demonstrate the importance of considering genomic and exposomic factors to understand psychosis development during puberty.
Study objectives and approach
In the present study, we investigated the impact of puberty on psychosis expression in female adolescents of the ABCD Study. To comprehensively understand puberty factors and triangulate our findings, we examined puberty status, timing and hormone levels, focusing on variables specific to females. Furthermore, as an explorative secondary aim, we investigated whether these processes might vary in the context of exposomic and genomic risk of schizophrenia.
Method
Participants
Data were extracted from the 5.1 release (DOI: 10.15154/z563-zd24) of the ABCD Study, encompassing assessments from 1 September 2016 to 15 January 2022. This ongoing longitudinal study in the USA consists of 11 868 children who were initially 9 to 10 years old at baseline. The recruitment process spanned 21 sites, employing a multi-stage probability sampling method to establish a school sample characterised by demographic diversity. Reference Garavan, Bartsch, Conway, Decastro, Goldstein and Heeringa23 The data collection was approved by the centralised institutional review board (IRB) at the University of California San Diego and local research site institutional review boards, and data access was facilitated by the National Institute of Mental Health National Data Archive. Written informed consent was obtained from participating parents and/or caregivers, along with assent from the youth.
The sample for this analysis was restricted to female participants who had completed the evaluation of psychotic experiences at least once across the three annual assessment points (i.e. baseline, first follow-up and second follow-up assessments). Individuals using hormonal contraception were excluded from analyses involving oestradiol (n = 24). In addition, as polygenic risk score performance varies across different populations owing to an imbalance favouring European ancestry, Reference Martin, Kanai, Kamatani, Okada, Neale and Daly24 only participants of European descent with high-quality genotyping data were included in the genetic analyses, consistent with prior research. Reference Di Vincenzo, Prachason, Sampogna, Arias-Magnasco, Lin and Pries21,Reference Karcher, Paul, Johnson, Hatoum, Baranger and Agrawal25
Measurements
Distressing psychotic experiences
The Prodromal Questionnaire–Brief Child Version, which has been previously validated for a school-age population, Reference Loewy, Pearson, Vinogradov, Bearden and Cannon26 was used to assess psychotic experiences at baseline and during follow-up assessments. Psychotic experiences, such as unusual thought content and perceptual abnormalities within the past month, were evaluated using a 21-item scale. Participants indicated the level of distress they experienced associated with the psychotic experiences on a five-point Likert scale (1 = Not very bothered; 2 = Slightly bothered; 3 = Moderately bothered; 4 = Very much bothered; 5 = Extremely bothered). In the current study, a binary variable indicating significantly distressing psychotic experiences was created to assess a clinically relevant manifestation of psychotic experiences, in line with previous reports in this data-set. Reference Di Vincenzo, Prachason, Sampogna, Arias-Magnasco, Lin and Pries21 Psychotic experiences were considered to be present if participants endorsed at least one psychotic experience accompanied by significant psychological distress, with a scoring threshold of ≥3 of the five points.
Perceived pubertal status and age at menarche
We used caregiver-rated measures for perceived puberty status and age at menarche, as caregiver ratings are considered more accurate than adolescent self-reports. At each assessment, puberty stages and, if applicable, age at menarche were recorded.
To evaluate perceived puberty status, we employed the Pubertal Development Scale, Reference Carskadon and Acebo27 with ratings provided by both the primary caregiver and the adolescent participant. To assess puberty status, five questions were used to inquire about growth spurt, body hair, skin, breast development and menarche. The physical markers were rated on a four-point Likert scale (1 = has not begun yet, 2 = barely begun, 3 = definitely begun, 4 = seems complete). Menarche was assessed first with a yes/no response, followed by a question to determine the age. Following previous research, we used a pubertal category score. Reference Carskadon and Acebo27 The scores for body hair growth and breast development were summed, and menarche was used to categorise as follows: pre-pubertal, score = 2 and no menarche; early pubertal, score = 3 and no menarche; mid-pubertal, score ≥ 3 and no menarche; late pubertal, score ≤7 and menarche; post-pubertal, score = 8 and menarche. Reference Carskadon and Acebo27 Given the young age of the cohort, there were very few individuals in the post-pubertal stage. Therefore, we combined the late and post-pubertal categories for the current analyses.
Age at menarche was evaluated separately to assess the timing of puberty. Implausible responses (e.g. 111) were treated as missing. When multiple reports were available across the three assessment points, we used the earliest recorded assessment. In cases of incomplete caregiver reports on age at menarche, adolescent self-reports were used as a supplement (n = 126). For analyses involving age at menarche, we included only individuals who had experienced their first period before the second follow-up assessment, and we standardised available data for these analyses.
Hormone measures
Salivary samples were obtained to measure oestradiol concentrations at baseline and each subsequent follow-up assessment. Comprehensive information on the collection process and sample preparation is available in an earlier report. Reference Uban, Horton, Jacobus, Heyser, Thompson and Tapert28 In short, saliva samples were collected by research assistants using the passive drool method. Participants were instructed to refrain from eating, chewing gum or drinking for 30 min before collection and not to have major meals 60 min before collection. After collection, saliva samples were directly cooled on ice and then sent to Salimetrics. Reference Uban, Horton, Jacobus, Heyser, Thompson and Tapert28 To avoid multiple freeze–thaw cycles, samples were assayed in duplicate within 1 day.
In the pre-processing steps, we adhered to the procedures outlined in the previous report. Reference Herting, Uban, Gonzalez, Baker, Kan and Thompson29 Only participants who provided a valid saliva sample were considered for hormonal analyses. We excluded participants if there was a discrepancy between their recorded salivary sex and sex reported at birth. In addition, any measures that were affected by hormone quality issues or fell outside the sensitivity limits of the assay were excluded (see ref. 29 for details). When two replicates were available, the mean value was used. Otherwise, the single value was used. Hormone data were standardised for each time point separately.
Following established procedures within the ABCD Study Reference Chaku and Barry30 and to address inconsistencies in time assessments, we recoded measures as missing under the following conditions: when the start time of data collection occurred before the wake-up time, when the end time of collection occurred before the start time, or when the freeze time preceded the end time of collection. If caregivers or participants reported that the participants were using hormonal contraceptives, they were excluded from the hormonal analyses. To ensure alignment with the cohort’s protocol and recommendations, Reference Chaku and Barry30 we restricted sample collection start times to between 07:00 and 19:00.
Exposome score for schizophrenia (ES-SCZ)
The ES-SCZ was computed based on a previous report. Reference Pries, Lage-Castellanos, Delespaul, Kenis, Luykx and Lin31 This score has previously been used in the ABCD Study. Reference Di Vincenzo, Prachason, Sampogna, Arias-Magnasco, Lin and Pries21 Nine environmental factors (emotional neglect, physical neglect, emotional abuse, physical abuse, sexual abuse, cannabis use, winter birth, hearing impairment and bullying) at baseline and two follow-up assessments were extracted from the ABCD Study data-set. More information on these exposures can be found in the Supplementary Material available at https://doi.org/10.1192/bjp.2025.36 and in a previous report. Reference Di Vincenzo, Prachason, Sampogna, Arias-Magnasco, Lin and Pries21 Binary variables (0 = absent; 1 = present) indicating lifetime exposure for each environmental factor were generated as detailed in the Supplementary Material. An aggregated weighted score was calculated by summing the nine exposures multiplied by their weighted risks (log odds) for schizophrenia. Reference Pries, Lage-Castellanos, Delespaul, Kenis, Luykx and Lin31
Polygenic risk score for schizophrenia (PRS-SCZ)
PRS-SCZ was constructed for participants who passed genetic and sample quality control. For comprehensive quality control steps and principal component analysis related to ancestry, please refer to the Supplementary Material. We generated PRS-SCZ using data from the most recent schizophrenia genome-wide association study (GWAS; European subsample) based on 53 386 cases and 77,258 controls, Reference Trubetskoy, Pardinas, Qi, Panagiotaropoulou, Awasthi and Bigdeli32 applying a Bayesian framework method that utilised continuous shrinkage on single-nucleotide polymorphism effect sizes. This method is robust to varying genetic architectures, provides substantial computational advantages and enables multivariate modelling of local linkage disequilibrium patterns. Reference Ge, Chen, Ni, Feng and Smoller33 We used the 1000 Genomes Project European Sample (https://github.com/getian107/PRScs) as a disequilibrium reference panel. To compute posterior effect sizes, we used the default settings (Supplementary Material). After calculation of posterior effect sizes, PRS-SCZ was calculated using the ‘--score’ function and the SUM modifier in PLINK1.9. Reference Purcell, Neale, Todd-Brown, Thomas, Ferreira and Bender34 After quality control, 742 011 variants were used in the PRS-SCZ calculation. PRS-SCZ data from a subsample of 2775 European females were used and PRS-SCZ was standardised.
Statistical analysis
All analyses in the current study were performed using Stata (release 18). 35 We applied separate longitudinal mixed logistic regression models with oestradiol concentration (baseline: n = 5163) or puberty status (baseline: n = 5456) as the independent variable and psychotic experiences as the dependent variable, using three assessment points: baseline and two annual follow-up assessments. To account for multiple testing in the longitudinal analyses, Bonferroni correction was applied (0.05/2 = 0.025). Age at menarche was analysed cross-sectionally in individuals who had reached menarche (n = 2947) at the second follow-up, with psychotic experiences at the second follow-up serving as the outcome variable. Missing values for the covariates can be found in Supplementary Table 2.
All models included random intercepts for site and family membership. The longitudinal analyses also accounted for individual assessments to manage multiple assessments per participant and further incorporated a random slope for assessment points. To investigate the effects of puberty factors within the context of environmental and genetic vulnerability to schizophrenia, we tested for multiplicative interactions with ES-SCZ and PRS-SCZ, respectively.
In alignment with previous literature, the models were adjusted for two sets of covariates: (1) age; and (2) age, family income, parental education, body mass index and ethnicity (except for genetic analyses, which were limited to a European subsample; for details, see the Supplementary Material). Hormonal analyses included additional methodological adjustments (wake-up time on collection day, start time of collection, duration of collection, and time from collection to freezing) and physiological adjustments (caffeine intake and engagement in vigorous physical exercise within the past 12 h; yes/no). Genetic analyses were adjusted for ten principal components.
Results
In this study, there were 5673 participants with data for baseline assessment, and 5348 and 5203 for 1-year and 2-year follow-up assessments, respectively. Table 1 presents sample characteristics at different time points, Supplementary Table 1 details the covariates and Supplementary Table 2 provides information on missing values.
Sample characteristics at different time points
ES-SCZ, exposome score for schizophrenia; NA, not applicable.
a Restricted to individuals with age ≥ age at menarche at the second follow-up assessment.
The cross-sectional analyses revealed that earlier age at menarche was significantly associated with a greater likelihood of reporting psychotic experiences at the second follow-up assessment (model 1: odds ratio 0.68, 95% CI: 0.59 to 0.78, P < 0.001; model 2: odds ratio 0.77, 95% CI: 0.67 to 0.88, P < 0.001). The longitudinal analyses indicated that increased standardised oestradiol concentration was associated with increased likelihood of psychotic experiences (model 1: odds ratio 1.08, 95% CI: 1.01 to 1.16, P = 0.020). The association was trend-significant (P = 0.096) in model 2 (Table 2).
Main effects of puberty factors and vulnerability to schizophrenia on psychotic experiences
Bold text indicates statistically significant results. ES-SCZ, exposome score for schizophrenia; PRS, polygenic risk score for schizophrenia; NA, not applicable.
a Additionally adjusted for salivary collection covariates.
b This analysis was conducted at the second follow-up assessment to align with the analysis of age at menarche.
c Genetic analyses were conducted within the European subsample and additionally adjusted for ten principal components.
According to the longitudinal analyses of puberty status, being in later puberty stages (compared with the pre-pubertal reference group) was associated with increased likelihood of experiencing psychotic experiences. Specifically, individuals in the mid-pubertal (model 1: odds ratio 1.41, 95% CI: 1.18 to 1.69, P < 0.001) and the late or post-pubertal stages (model 1: odds ratio 2.23, 95% CI: 1.74 to 2.86, P < 0.001) had more psychotic experiences. The difference between the pre-pubertal and early pubertal stages was not statistically significant (model 1: odds ratio 1.17, 95% CI: 0.96 to 1.41, P = 0.114). Follow-up contrast analyses further indicated differences between the early pubertal and mid-pubertal stages (model 1: χ2(1) = 5.31, P = 0.021), the early pubertal and late or post-pubertal stages (model 1: χ²(1) = 31.31, P < 0.001) and the mid-pubertal and late or post-pubertal stages (model 1: χ²(1) = 26.92, P < 0.001).
Results for model 2 were similar (Fig. 1 and Table 2); however, the difference between the pre-pubertal (reference group) and mid-pubertal stages was no longer significant. Follow-up contrast analyses further indicated differences between the early pubertal and late or post-pubertal stages (model 2: χ²(1) = 6.73, P = 0.010), as well as between the mid-pubertal and late or post-pubertal stages (model 2: χ²(1) = 9.32, P = 0.002). The difference between the early pubertal and mid-pubertal stages was no longer significant (model 2: χ²(1) = 0.18, P = 0.670).
Effects of perceived puberty status on psychotic experiences. Model 1: adjusted for age; model 2: adjusted for age, body mass index, family income, parental education and ethnicity.
ES-SCZ was significantly associated with psychotic experiences in the cross-sectional analyses (model 1: odds ratio 1.61, 95% CI: 1.39 to 1.85, P < 0.001) and longitudinal analyses (model 1: odds ratio 1.86, 95% CI: 1.69 to 2.05, P < 0.001). The results converged with those of model 2 (Table 2). Furthermore, PRS-SCZ was significantly associated with psychotic experiences in the cross-sectional analyses (model 1: odds ratio 1.32, 95% CI: 1.05 to 1.66, P = 0.017) and the longitudinal analyses (model 1: odds ratio 1.36, 95% CI: 1.18 to 1.58; P < 0.001). The results for model 2 were similar; however, the cross-sectional analyses using model 2 showed a trend-significant association (P = 0.060, Table 2). There were no statistically significant interactions of ES-SCZ or PRS-SCZ with puberty factors (Table 3).
Interaction effects of puberty factors and vulnerability to schizophrenia on psychotic experiences
ES-SCZ, exposome score for schizophrenia; PRS-SCZ: polygenic risk score for schizophrenia; NA, not applicable.
a This analysis was conducted at the second follow-up assessment to align with the analysis of age at menarche.
b Additionally adjusted for salivary collection covariates.
c Genetic analyses were conducted within the European subsample and additionally adjusted for ten principal components.
Discussion
To the best of our knowledge, this study represents the first comprehensive exploration of the impact of puberty factors (timing, perceived status and oestradiol concentration) on psychotic experiences in a female adolescent general population, particularly within the context of exposomic and genomic susceptibility to schizophrenia. Earlier age at menarche and higher oestradiol concentration were associated with a greater likelihood of psychotic experiences. Furthermore, being in a later perceived puberty stage was associated with psychotic experiences. Both increased genetic and environmental vulnerability to schizophrenia were associated with the likelihood of psychotic experiences. However, we found no interactions between puberty factors and exposomic or genomic risk for schizophrenia with respect to the effects of these factors on psychosis.
The finding of the present study that earlier age at menarche was associated with an increased likelihood of experiencing psychotic experiences during adolescence aligns with studies linking early puberty to internalising and externalising symptoms. Reference Ullsperger and Nikolas6 However, there is also some evidence that early menarche may be protective against severe outcomes in patients with PSD, although the relationship remains inconclusive in the literature. Reference Hochman and Lewine7-Reference Magnus, Guyatt, Lawn, Wyss, Trajanoska and Küpers12 Notably, most previous studies have involved small sample sizes and primarily focused on patients and older individuals. The relationship between puberty timing and psychosis may differ in younger individuals from the general population compared with adults with severe clinical conditions. Although early puberty might represent a transdiagnostic risk factor for various mental health issues during adolescence, it could potentially still be protective in patients with PSD, especially when it is linked to oestradiol concentration. It is also plausible that as time progresses, immediate sociobiological effects diminish, making way for more enduring neurobiological modifications. Furthermore, some studies suggest that both early and late puberty may be associated with increased risk of psychosis in the adult general population. Reference Kaiser and Gruzelier36 The current sample comprised relatively young individuals, with a maximum age of 14 at the second follow-up assessment. Furthermore, only a subsample had undergone menarche by the second follow-up or had information on age at menarche at this assessment point. Consequently, a comprehensive follow-up that includes individuals with later age at menarche is necessary.
The current findings indicate that higher oestradiol concentrations are associated with a greater likelihood of experiencing psychotic experiences. This may indicate that higher oestradiol concentrations are not universally protective against psychosis across the lifespan. This divergence from the oestrogen hypothesis, in combination with the finding on age at menarche, could be attributed to the complex nature of puberty, in which factors have both immediate- and long-term effects, influencing neurobiological mechanisms that may be activated later in life. Notably, a recent ABCD study Reference Damme, Hernandez and Mittal13 demonstrated an association between menarche status (present versus not present) at baseline and first follow-up and the incidence of psychotic experiences. Individuals who were post-menarche at both time points exhibited the highest likelihood of psychotic experiences, followed by those transitioning from pre- to post-menarche and, finally, individuals who were pre-menarche at both time points. That study also reported that heightened hippocampal connectivity, which has previously been linked to increased oestradiol exposure and neurological protection in high-risk individuals, Reference Damme, Ristanovic, Vargas and Mittal14 was associated with fewer psychotic experiences only in those who remained pre-menarche at both time points. These findings, combined with ours, raise important questions about when high oestradiol concentration indicates vulnerability or a protective effect and how the timing and status of puberty affect the individual.
The current study also found that being in a later stage of perceived puberty, as opposed to the pre-pubertal phase, was associated with more reports of psychotic experiences, indicating a specific vulnerability during the later stages. These findings align with research showing that puberty is a critical vulnerability period during which many mental health problems emerge. Notably, although being in a later stage of puberty increased the risk of psychotic experiences in this sample, studies have suggested that older age is generally associated with a reduced risk of such experiences. Younger individuals report psychotic experiences more frequently, with a decline in prevalence observed from childhood through adolescence and into adulthood. Reference Yates, Lang, Peters, Wigman, McNicholas and Cannon37,Reference Kelleher, Connor, Clarke, Devlin, Harley and Cannon38 Although experiencing psychotic experiences is linked to the development of PSD, these experiences may not independently indicate a risk of later PSD in many individuals. Therefore, future studies should aim to understand how multiple factors – such as the timing of puberty and the occurrence of psychotic experiences during different puberty stages – affect the development of more severe clinical outcomes later in life. Notably, the age at which first psychotic experiences occur may also interact with puberty factors to affect severe outcomes, adding further complexity to this relationship.
The effects of sex (biological) and gender (sociocultural) on mental health during puberty are intricately intertwined. In addition to biological changes, researchers have theorised that the effects of puberty on internalising and externalising problems in girls could be attributed to different psychological mechanisms. Reference Pfeifer and Allen39 For instance, developmental readiness suggests that individuals are not prepared for changes during puberty. Maturational deviance posits that early or late developers diverge from the typical developmental timeline, which can lead to adverse social reactions. For instance, research suggests that teachers may have lower academic expectations of girls displaying early secondary sex characteristics. Reference Carter, Mustafaa and Leath40 In the current study, we investigated both visible physical developmental changes, which may provoke sociocultural reactions, and hormonal levels, which are more non-visual biological changes. This is crucial, as a recent ABCD study found that although hormone levels and perceived physical markers of puberty showed moderate correlations, they represented distinct aspects of development. Reference Herting, Uban, Gonzalez, Baker, Kan and Thompson29 Another crucial aspect is brain development during puberty, which can be influenced by hormonal fluctuations, although the mechanisms underlying this influence are not fully understood. Future research should explore how puberty factors relate to PSD symptoms through brain adaptations.
This study confirms that exposomic and genomic risk factors for schizophrenia influence psychotic experiences in female general-population adolescents. These findings support previous research on the critical roles of environmental and genetic risk factors across the psychosis spectrum. Reference Di Vincenzo, Prachason, Sampogna, Arias-Magnasco, Lin and Pries21 However, no interactions were found between genetic or environmental vulnerability to schizophrenia and puberty factors. These findings should be interpreted cautiously owing to the limited sample size for these analyses and the young age of participants. Previous studies Reference Vijayakumar, Youssef, Bereznicki, Dehestani, Silk and Whittle22 have indicated an interaction between environmental risk and puberty with respect to effects on mental health problems, underscoring the multifactorial nature of such problems. Future follow-up analyses will be needed as more clinically relevant phenotypes emerge. Although current GWAS evidence indicates that sex-specific PRS-SCZ is unlikely to differ significantly from general PRS-SCZ, Reference Trubetskoy, Pardinas, Qi, Panagiotaropoulou, Awasthi and Bigdeli32 future advances in this field may offer valuable insights. In addition, an important area for research would be to explore how sex-specific vulnerabilities influence help-seeking behaviour of those at risk of psychosis.
This study represents a comprehensive evaluation of the role of puberty factors in psychotic experiences in female adolescents. Nonetheless, there are several limitations to consider. First, the young age of the participants may have limited the findings on interaction effects, as many individuals had not yet completed all stages of puberty or provided information on age at menarche. It is possible that these effects will become more evident as participants age and more severe clinical phenotypes emerge. Likewise, the analyses involving age at menarche were restricted to individuals who had experienced menarche before the second follow-up assessment; therefore, they could have overlooked different effects that might emerge among those with a later age at menarche. Furthermore, individuals with later menarche may not have reached advanced puberty in this cohort; this could have affected associations between puberty stages and psychotic experiences. However, significant differences among pre-menarche stages (pre- versus mid-pubertal, early versus mid-pubertal) suggest that pubertal changes may influence psychotic experiences risk beyond menarche status. Second, the small sample size, particularly concerning individuals with information on age at menarche, follow-up assessments and genetic data, might have led to low statistical power for interaction effects. Third, following prior research, Reference Di Vincenzo, Prachason, Sampogna, Arias-Magnasco, Lin and Pries21 we limited genetic analyses to the European subpopulation to optimise PRS-SCZ performance with ancestry-matched GWAS data. This improves precision for the European group but limits generalisability to those of non-European ancestries. Future studies should adopt approaches that include diverse ancestry groups to address these limitations. Fourth, the ABCD Study solely focused on salivary measures of the free form of oestradiol in females, constraining our exploration of the broader oestrogenic landscape. Research on the connection between oestrogen and PSD has predominantly focused on oestradiol during reproductive years, and certain protective effects have been attributed to this form of oestrogen. A previous ABCD study linked cognitive functioning to hormone profiles, including testosterone, dehydroepiandrosterone and oestradiol. Given the importance of cognition in PSD, this is a valuable approach that could also be applied to psychosis in future studies. Fifth, various factors, such as other hormones (e.g. cortisol), daily hormonal fluctuations and menstrual cyclicity may influence hormone levels. We aligned our methodology with previous covariate approaches Reference Uban, Horton, Jacobus, Heyser, Thompson and Tapert28,Reference Chaku and Barry30 and used available variables. However, future studies could benefit from exploring additional factors that may affect oestradiol concentrations and should aim to gain deeper insights into the relationship between puberty factors and psychosis, particularly in the context of more severe outcomes.
Finally, the findings of the current study regarding the important influence of puberty on the risk of psychosis expression during adolescence emphasise the need for comprehensive follow-up research that examines both the immediate- and long-term impacts of puberty on mental health outcomes across diverse populations. It is crucial for such research to explore male-specific aspects of puberty to gain a more complete understanding of sex-specific effects.
Supplementary material
Supplementary material is available online at https://doi.org/10.1192/bjp.2025.36
Acknowledgements
Data used in the preparation of this article were obtained from the ABCD Study (https://abcdstudy.org), held in the National Institute of Mental Health Data Archive. This is a multisite, longitudinal study designed to recruit more than 11 500 children aged 9–10 years and follow them over 10 years into early adulthood. The ABCD Study is supported by the National Institutes of Health and additional federal partners under grant numbers U01DA041022, U01DA041025, U01DA041028, U01DA041048, U01DA041089, U01DA041093, U01DA041106, U01DA041117, U01DA041120, U01DA041134, U01DA041148, U01DA041156, U01DA041174, U24DA041123 and U24DA041147. A full list of supporters is available at https://abcdstudy.org/nih-collaborators. A listing of participating sites and a complete listing of the study investigators can be found at https://abcdstudy.org/principalinvestigators.html. ABCD consortium investigators designed and implemented the study and/or provided data but did not necessarily participate in the analysis of the data or the writing of this report. This manuscript reflects the views of the authors and may not reflect the opinions or views of the National Institutes of Health or ABCD consortium investigators. The ABCD data repository grows and changes over time. The ABCD data used in this report were from https://dx.doi.org/10.15154/1523041 (data release 4.0) and https://dx.doi.org/10.15154/z563-zd24 (data release 5.1).
Author contributions
L.-K.P. conceptualised the study and its design, performed data analyses, interpreted the results and wrote the first draft of the manuscript. T.P. contributed to the analyses by precleaning the data, reviewed the content and provided final approval for publication. A.A.-M. prepared the genetic data for analyses, reviewed the content and provided final approval for publication. B.D.L. prepared the genetic data for analyses, reviewed the content and provided final approval for publication. B.P.F.R. contributed to the acquisition of the data, reviewed the content and provided final approval for publication. S.G. contributed to the acquisition of funding and data, development of the study design, and drafting and reviewing the manuscript, as well as providing final approval.
Funding
S.G. is supported by the Ophelia research project, ZonMw grant 636340001. B.P.F.R. was funded by a Vidi award (91718336) from The Netherlands Scientific Organisation. S.G., B.R., L.-K.P., B.D.L. and A.A.-M. are supported by the YOUTH-GEMs project, funded by the European Union’s Horizon Europe programme under grant agreement number 101057182.
Declaration of interest
S.G. serves as a member of the British Journal of Psychiatry editorial board. However, he was not involved in the review or decision-making process for this manuscript.
Open access
eLetters
No eLetters have been published for this article.