Introduction
Eating disorders (EDs), obsessive-compulsive disorder (OCD), and schizophrenia (SCZ) are serious psychiatric illnesses (American Psychiatric Association, 2013; Keshaviah et al., Reference Keshaviah, Edkins, Hastings, Krishna, Franko, Herzog, Thomas, Murray and Eddy2014; Pearlson, Reference Pearlson2000; Treasure, Duarte, & Schmidt, Reference Treasure, Duarte and Schmidt2020; World Health Organization, 2008) that impact an individual’s overall well-being. Epidemiological, clinical, and genetic evidence support an overlap between EDs, OCD, and SCZ (Cederlöf et al., Reference Cederlöf, Thornton, Baker, Lichtenstein, Larsson, Rück, Bulik and Mataix-Cols2015; Cross-Disorder Group of the Psychiatric Genomics Consortium, 2019; Foulon, Reference Foulon2003; Godart, Flament, Perdereau, & Jeammet, Reference Godart, Flament, Perdereau and Jeammet2002; Halmi et al., Reference Halmi, Tozzi, Thornton, Crow, Fichter, Kaplan, Keel, Klump, Lilenfeld, Mitchell, Plotnicov, Pollice, Rotondo, Strober, Woodside, Berrettini, Kaye and Bulik2005; Hudson & Jonas, Reference Hudson and Jonas1984, Hudson, Hiripi, Pope, & Kessler, Reference Hudson, Hiripi, Pope and Kessler2007; Khalil, Hachem, & Richa, Reference Khalil, Hachem and Richa2011; Morylowska-Topolska et al., Reference Morylowska-Topolska, Ziemiński, Molas, Gajewski, Flis, Stelmach and Karakuła-Juchnowicz2017; Swinbourne & Touyz, Reference Swinbourne and Touyz2007; Yilmaz et al., Reference Yilmaz, Halvorsen, Bryois, Yu, Thornton, Zerwas, Micali, Moessner, Burton, Zai, Erdman, Kas, Arnold, Davis, Knowles, Breen, Scharf, Nestadt, Mathews and Crowley2020).
Epidemiological observations show that EDs and OCD commonly co-occur in the clinic (Godart et al., Reference Godart, Flament, Perdereau and Jeammet2002; Halmi et al., Reference Halmi, Tozzi, Thornton, Crow, Fichter, Kaplan, Keel, Klump, Lilenfeld, Mitchell, Plotnicov, Pollice, Rotondo, Strober, Woodside, Berrettini, Kaye and Bulik2005; Hudson et al., Reference Hudson, Hiripi, Pope and Kessler2007; Swinbourne & Touyz, Reference Swinbourne and Touyz2007). It is estimated that between 20–60% of patients with EDs have a lifetime history of OCD (Godart et al., Reference Godart, Flament, Curt, Perdereau, Lang, Venisse, Halfon, Bizouard, Loas, Corcos, Jeammet and Fermanian2003; Halmi et al., Reference Halmi, Tozzi, Thornton, Crow, Fichter, Kaplan, Keel, Klump, Lilenfeld, Mitchell, Plotnicov, Pollice, Rotondo, Strober, Woodside, Berrettini, Kaye and Bulik2005; Kaye, Bulik, Thornton, Barbarich, & Masters, Reference Kaye, Bulik, Thornton, Barbarich and Masters2004) and that between 3 and 13% of patients with OCD have a lifetime history of an ED (du Toit, van Kradenburg, Niehaus, & Stein, Reference du Toit, van Kradenburg, Niehaus and Stein2001; Fireman, Koran, Leventhal, & Jacobson, Reference Fireman, Koran, Leventhal and Jacobson2001; LaSalle et al., Reference LaSalle, Cromer, Nelson, Kazuba, Justement and Murphy2004; D. A. Pinto, Mancebo, Eisen, Pagano, & Rasmussen, Reference Pinto, Mancebo, Eisen, Pagano and Rasmussen2006; Sallet et al., Reference Sallet, de Alvarenga, Ferrão, de Mathis, Torres, Marques, Hounie, Fossaluza, Rosario, Fontenelle, Petribu and Fleitlich-Bilyk2010). Studies have also shown that in both clinical and population samples, SCZ psychotic symptoms and EDs are comorbid (Foulon, Reference Foulon2003; Hudson & Jonas, Reference Hudson and Jonas1984; Khalil et al., Reference Khalil, Hachem and Richa2011; Morylowska-Topolska et al., Reference Morylowska-Topolska, Ziemiński, Molas, Gajewski, Flis, Stelmach and Karakuła-Juchnowicz2017; Solmi, Mascarell, Zammit, Kirkbride, & Lewis, Reference Solmi, Mascarell, Zammit, Kirkbride and Lewis2019a).
EDs, particularly Anorexia Nervosa (AN), OCD, and SCZ, also share phenotypic traits such as perfectionism (Bardone-Cone et al., Reference Bardone-Cone, Wonderlich, Frost, Bulik, Mitchell, Uppala and Simonich2007; Bulik et al., Reference Bulik, Tozzi, Anderson, Mazzeo, Aggen and Sullivan2003), neuroticism (Cassin & von Ranson, Reference Cassin and von Ranson2005; Lilenfeld, Reference Lilenfeld2011; Samuels et al., Reference Samuels, Nestadt, Bienvenu, Riddle, Liang, Hoehn-Saric, Grados and Cullen2000), obsessive thoughts and repetitive behaviors (Braun, Sunday, & Halmi, Reference Braun, Sunday and Halmi1994; George Hsu, Kaye, & Weltzin, Reference George Hsu, Kaye and Weltzin1993; Halmi et al., Reference Halmi, Eckert, Marchi, Sampugnaro, Apple and Cohen1991; Serpell, Livingstone, Neiderman, & Lask, Reference Serpell, Livingstone, Neiderman and Lask2002). In addition to these shared characteristics, patients with both AN and SCZ often present with cognitive disturbances (Morylowska-Topolska et al., Reference Morylowska-Topolska, Ziemiński, Molas, Gajewski, Flis, Stelmach and Karakuła-Juchnowicz2017). While this disturbance is generalized to the entire cognitive processes in SCZ, AN patients show impaired reasoning in regard to body image, body weight, and eating behavior (Morylowska-Topolska et al., Reference Morylowska-Topolska, Ziemiński, Molas, Gajewski, Flis, Stelmach and Karakuła-Juchnowicz2017; Powers, Simpson, & McCormick, Reference Powers, Simpson and McCormick2005). It is thought that the observed maladaptive behaviors in ED, OCD, and SCZ reflect altered brain structure or function. Adult AN patients have been shown to have a reduced total intracranial volume (TIV) (Van den Eynde et al., Reference Van den Eynde, Suda, Broadbent, Guillaume, Van den Eynde, Steiger, Israel, Berlim, Giampietro, Simmons, Treasure, Campbell and Schmidt2012), grey matter volume (GMV), and an increased cerebrospinal fluid (CSF) volume compared to controls (Titova, Hjorth, Schiöth, & Brooks, Reference Titova, Hjorth, Schiöth and Brooks2013). Similarly, reduced GMV and cortical thickness have been identified as hallmarks of SCZ (Howes, Cummings, Chapman, & Shatalina, Reference Howes, Cummings, Chapman and Shatalina2023), while impaired functioning of the frontal cortex is commonly observed in patients with OCD (van den Heuvel et al., Reference van den Heuvel, Veltman, Groenewegen, Cath, van Balkom, van Hartskamp, Barkhof and van Dyck2005) Figure 1 highlights shared and disorder-specific cognitive and illness traits across AN, OCD, and SCZ. It is currently unclear whether the cognitive characteristics observed in AN - that are similar to those seen in OCD and SCZ – are state (i.e. follow the clinical diagnosis: the effect of the illness) or trait (i.e. precede clinical diagnosis) manifestations.
Anorexia nervosa, obsessive compulsive disorder, and schizophrenia characteristics (white circles) and transdiagnostic traits (grey shaded circles). Within each circle, disorder symptoms are identified in bold, cognitive characteristics in italics, and neural characteristics in regular text.
Evidence from genetic studies in the last 5–10 years suggests that these three psychiatric conditions (AN, OCD, SCZ) have a common biological underpinning, suggesting that an overlapping set of genetic markers confers risk for the three disorders. Previous investigations identified the AN-OCD genetic correlation as the 2nd highest among eight analyzed psychiatric conditions (r = 0.50) and reported a moderate genetic correlation between OCD and SCZ (r = 0.35) (Cross-Disorder Group of the Psychiatric Genomics Consortium, 2019). A recent study (Lu et al., Reference Lu, Ploner, Birgegård, Bulik and Bergen2024) further provided evidence for a genetic overlap between AN and SCZ, identifying 10 genomic loci associated with both disorders. This suggests that genetic risk for one of these psychiatric disorders can increase risk for another.
Polygenic risk scores (PRSs) index an individual’s genetic propensity for a specific disease (Lewis & Vassos, Reference Lewis and Vassos2020). Results from our group indicate that both PRS for AN and PRS for OCD can inform about AN cognition and behaviours in adolescents in a population-based cohort (Yilmaz et al., Reference Yilmaz, Schaumberg, Halvorsen, Goodman, Brosof, Crowley, Mathews, Mattheisen, Breen, Bulik, Micali and Zerwas2023). This suggests that genetic risk specific to a given psychiatric disorder can inform about common phenotypic traits.
The literature on AN and psychiatric comorbidities mainly focuses on clinical manifestations of these disorders, disregarding shared symptoms and behaviours across disorders that precede clinical onset (Yilmaz et al., Reference Yilmaz, Schaumberg, Halvorsen, Goodman, Brosof, Crowley, Mathews, Mattheisen, Breen, Bulik, Micali and Zerwas2023). Moreover, the majority of studies include more severe cases (e.g. clinical studies, case–control studies), which prevents generalization to the general population. However, given that 85% of AN cases are said to begin before the age of 20, with almost all cases starting before the age of 25 (Herpertz-Dahlmann, Reference Herpertz-Dahlmann2009), it is important to understand how psychiatric traits are expressed in a population at risk (i.e. 8–21 years old) without a clinical diagnosis.
The present study aims to investigate whether polygenic risk for AN, OCD, and SCZ associates with disorder-specific and shared clinical, neurocognitive, and cortical phenotypes in a non-clinical population sample. Given shared genetics between AN, OCD, and SCZ, we hypothesized that PRSs for AN, OCD, and SCZ are associated with clinical, neurocognitive, and cortical characteristics seen in these disorders. Moreover, we hypothesized that disease-specific polygenic risk to be more informative about disease-specific cognitive and clinical phenotypes compared to non-disease-specific PRS. Lastly, in order to identify cortical markers of AN phenotypes, we performed an additional exploratory analysis focusing on the relationship between AN clinical phenotypes and cortical morphometry. Based on previous clinical reports of cortical changes in AN patients (Frank et al., Reference Frank2015), we expected AN clinical phenotypes to associate with cortical markers. Given the sex differences in psychiatric disorders (Yang et al., Reference Yang, Fang, Arnberg, Kuja-Halkola, D’Onofrio, Larsson, Brikell, Chang, Andreassen, Lichtenstein, Valdimarsdóttir and Lu2024), we performed exploratory sex-stratified analyses for all reported investigations. Considering the cross-sectional design of the study, we hereunder use ‘state’ to refer to one-time measurements and ‘trait’ to refer to phenotypes that are genetically determined.
Methods and materials
Participants
The Philadelphia Neurodevelopmental Cohort (PNC) is a population-based study of youth living in the greater Philadelphia (United States) area (Calkins et al., Reference Calkins, Merikangas, Moore, Burstein, Behr, Satterthwaite, Ruparel, Wolf, Roalf, Mentch, Qiu, Chiavacci, Connolly, Sleiman, Gur, Hakonarson and Gur2015) that aims to describe how genetics impacts brain development and cognitive functioning in adolescents (Satterthwaite et al., Reference Satterthwaite, Elliott, Ruparel, Loughead, Prabhakaran, Calkins, Hopson, Jackson, Keefe, Riley, Mensh, Sleiman, Verma, Davatzikos, Hakonarson, Gur and Gur2014).
Between 2006 and 2012, 50,293 participants were recruited by the Children’s Hospital of Philadelphia Center for Applied Genomics through a network of non-psychiatric pediatric clinics. To be included in the study, participants had to be between 8 and 21 years old, provide written informed consent to be re-contacted for future studies, be proficient in English, and not have any significant developmental delay or physical conditions that would interfere with study completion (Calkins et al., Reference Calkins, Merikangas, Moore, Burstein, Behr, Satterthwaite, Ruparel, Wolf, Roalf, Mentch, Qiu, Chiavacci, Connolly, Sleiman, Gur, Hakonarson and Gur2015). Overall, 19,161 individuals met the inclusion criteria; 9,428 (49.2%) of these completed clinical and neurocognitive assessments. Data cleaning was performed only on genotyped individuals (N = 8,139); 650 individuals were removed due to relatedness, and 2,760 individuals were removed due to non-European ancestry (2,050 individuals of African descent and about 700 admixed individuals). A total of 4,729 participants with genotypic data and at least one clinical or neurocognitive measure remained eligible for analysis.
Measures
Clinical measures
All participants underwent a clinical assessment using a computerized structured interview (Kiddie Schedule for Affective Disorders and Schizophrenia (K-SADS) Family Study Interview (Kaufman et al., Reference Kaufman, Birmaher, Brent, Rao, Flynn, Moreci, Williamson and Ryan1997) that was administered to an informant (caregiver or legal guardian), for those between 8 and 10 years old, to an informant and probands, for those between 11 and 17 years old, and to probands only for those aged 18–21 years old. The interview assessed demographics, medical history, and psychopathology. For our study, we focused on the items assessing AN (2 items), OCD (based on previous studies on this cohort, dimension scores were generated (Bralten et al., Reference Bralten, Widomska, Witte, Yu, Mathews, Scharf, Buitelaar, Crosbie, Schachar, Arnold, Lemire, Burton, Franke and Poelmans2020), and SCZ (5 items assessing auditory, visual, olfactory, tactile hallucinations and delusional beliefs) (Calkins et al., Reference Calkins, Merikangas, Moore, Burstein, Behr, Satterthwaite, Ruparel, Wolf, Roalf, Mentch, Qiu, Chiavacci, Connolly, Sleiman, Gur, Hakonarson and Gur2015). Questions about symptom frequency, symptom/episode duration, onset/offset of symptoms/episodes were excluded from the analyses. See Supplementary Table 1 for a list of the K-SADS items included.
Cognitive measures
Cognitive performance was assessed using the computerized neurocognitive battery (CNB) (Gur et al., Reference Gur, Richard, Hughett, Calkins, Macy, Bilker, Brensinger and Gur2010). The CNB assesses accuracy and speed in abstraction and mental flexibility, attention, working memory, episodic memory (e.g. word, face, spatial), language reasoning, spatial processing, sensorimotor, motor speed, and emotion identification (see Gur et al., Reference Gur, Richard, Hughett, Calkins, Macy, Bilker, Brensinger and Gur2010 for a complete review). The tests cover four cognitive domains (executive control, episodic memory, complex cognition, and social cognition) (Roalf et al., Reference Roalf, Ruparel, Verma, Elliott, Gur and Gur2013). A standardized reading test (Wide Range Achievement Test, WRAT4 (Wilkinson & Robertson, Reference Wilkinson and Robertson2006) was also administered prior to the CNB to test participants’ Intelligence Quotient (IQ) and their ability to complete the full battery.
The neurocognitive scores were computed following previous reports (Moore, Reise, Gur, Hakonarson, & Gur, Reference Moore, Reise, Gur, Hakonarson and Gur2015). Raw accuracy and speed scores were converted to z-scores based on the sample’s mean and standard deviation. Speed values were computed by multiplying by −1 the z-scores for median response time. Accuracy scores and speed values were then summed to create an efficiency score. For example, an individual with an accuracy score of 2.50 (i.e. very accurate) and a speed score of −2.50 (i.e. very slow) had an overall efficiency score of 0. Domain-specific scores were calculated as the average of different efficiency scores. The executive functioning domain score was calculated as the average of three efficiency scores: abstraction, attention, and working memory. The episodic memory domain score was computed as the average of verbal, face, and spatial memory efficiency scores. The complex cognition score reflected the average of language and analogical reasoning, nonverbal reasoning, as well as spatial processing efficiency scores. The social cognition domain score was computed as the average of emotion processing, emotion differentiation, and age differentiation efficiency scores (Supplementary Table 2).
Cortical measures
A subsample of participants (N = 1,598) underwent a Magnetic Resonance Imaging (MRI) scan on a 3 T Siemens TIMTrio scanner at the Hospital of the University of Pennsylvania (see Satterthwaite et al., Reference Satterthwaite, Elliott, Ruparel, Loughead, Prabhakaran, Calkins, Hopson, Jackson, Keefe, Riley, Mensh, Sleiman, Verma, Davatzikos, Hakonarson, Gur and Gur2014 for all the details regarding the neuroimaging protocol). A sub-analysis was carried out only on participants with complete data (i.e. genetic, clinical, neurocognitive, and neuroimaging data) (n = 627). Of the 627 subjects with complete data, 626 had a Magnetization Prepared Rapid Gradient Echo Imaging (MPRAGE) T1-weighted structural MRI image available, acquired with Repetition Time (TR) = 1810 ms, Echo Time (TE) = 3.5 ms, slice thickness = 1 mm, flip angle = 9°, field-of-view (FoV) = 180 mm (Satterthwaite et al., Reference Satterthwaite, Elliott, Ruparel, Loughead, Prabhakaran, Calkins, Hopson, Jackson, Keefe, Riley, Mensh, Sleiman, Verma, Davatzikos, Hakonarson, Gur and Gur2014). The MPRAGE image was used to extract GMV, white matter volume (WMV), CSF, TIV, cortical thickness, and gyrification.
Genotype data
The Center for Applied Genomics performed the genotyping of the entire PNC cohort (Satterthwaite et al., Reference Satterthwaite, Elliott, Ruparel, Loughead, Prabhakaran, Calkins, Hopson, Jackson, Keefe, Riley, Mensh, Sleiman, Verma, Davatzikos, Hakonarson, Gur and Gur2014) (for a complete overview of the genotyping pipeline, see Glessner et al., Reference Glessner, Wang, Cai, Korvatska, Kim, Wood, Zhang, Estes, Brune, Bradfield, Imielinski, Frackelton, Reichert, Crawford, Munson, Sleiman, Chiavacci, Annaiah, Thomas and Hakonarson2009; D. Pinto et al., Reference Pinto, Pagnamenta, Klei, Anney, Merico, Regan, Conroy, Magalhaes, Correia, Abrahams, Almeida, Bacchelli, Bader, Bailey, Baird, Battaglia, Berney, Bolshakova, Bölte and Piven2010). The genotyping was performed on one of the three Illumina platforms (HumanHap550, HumanHap610, or OmniExpress v2 (Robinson et al., Reference Robinson, Kirby, Ruparel, Yang, McGrath, Anttila, Neale, Merikangas, Lehner, Sleiman, Daly, Gur, Gur and Hakonarson2015). All data were retrieved through the National Institute of Mental Health’s Database of Genotypes and Phenotypes (dbGaP) (dbGaP, 2014) (PNC).
Missing data
Given that less than 20% of the data per variable was missing at random (range: 0.7–4.5%), multivariate imputation by chained equations (MICE) (mice package in R (Buuren & Groothuis-Oudshoorn, Reference Buuren and Groothuis-Oudshoorn2011) was used to estimate missing values for the clinical and neurocognitive variables. A prediction matrix was created to set which variables inform the imputation of which variable (Supplementary Table 8). Thus, clinical variables only informed the imputation of other clinical variables, while neurocognitive variables only informed the imputation of other neurocognitive variables. Biological variables such as sex, height, weight, and PRSs were not imputed. No data imputation was run on the neuroimaging dataset. Dichotomous clinical variables were imputed using the Bayesian logistic regression method (Rubin, Reference Rubin1987), while numeric cognitive variables were imputed using the predictive mean matching method. Five imputations were first run on the entire sample as the default setting in the mice package. The number of imputations was stepwise increased until the largest fraction of missing information (fmi) was 100 times smaller than the number of imputations (Madley-Dowd, Hughes, Tilling, & Heron, Reference Madley-Dowd, Hughes, Tilling and Heron2019). Seven imputations were needed to satisfy this criterion. As a second stage, imputations were run separately on female- and male-only datasets. Eight imputations were required to meet the fmi condition. Summary scores (e.g. obsession, compulsion, executive component, episodic memory, complex cognition, and social cognition scores) were created after the imputation of the raw variables.
The quality of the imputation was assessed via diagnostic plots (e.g. plots of the mean and standard deviation of the imputed values against the iteration number).
Statistical analyses
Polygenic risk scores
PRSs were computed using the polygenic risk score continuous shrinkage (PRS-CS) software (Ge, Chen, Ni, Feng, & Smoller, Reference Ge, Chen, Ni, Feng and Smoller2019). The discovery sample consisted of previously released GWAS summary statistics (Freeze 2 AN GWAS (Watson et al., Reference Watson, Yilmaz, Thornton, Hübel, Coleman, Gaspar, Bryois, Hinney, Leppä, Mattheisen, Medland, Ripke, Yao, Giusti-Rodríguez, Hanscombe, Purves, Adan, Alfredsson, Ando and Bulik2019), Freeze 1 OCD GWAS (Arnold et al., Reference Arnold, Askland, Barlassina, Bellodi, Bienvenu, Black, Bloch, Brentani, Burton, Camarena, Cappi, Cath, Cavallini, Conti, Cook, Coric, Cullen, Cusi and Davis2018), and Freeze 3 SCZ GWAS (Schizophrenia Working Group of the Psychiatric Genomics Consortium, 2022) and the target sample consisted of genotyped PNC participants of European ancestry.
Neuroimaging
All MPRAGE structural scans were preprocessed using the standard pipeline included in the Computational Anatomy Toolbox v. 12 (CAT12) (https://neuro-jena.github.io/cat/) and the Statistical Parametric Mapping (SPM12) toolbox (https://www.fil.ion.ucl.ac.uk/spm/software) in MATLAB (The Mathworks). Given the age range of our subjects (i.e. 8–21), a customized tissue probability map created with the CerebroMatic toolbox (Wilke, Altaye, Holland, & Consortium, Reference Wilke, Altaye, Holland and Consortium2017) was used during preprocessing. The toolbox contains a dataset of more than 1000 brain scans of subjects between the ages of 13 months to 75 years old. The customization was done based on participants’ age and sex, as well as the strength of the MRI magnetic field. GMV, WMV, CSF, and TIV were extracted from CAT12 for each participant, while cortical thickness and gyrification values were extracted from each region-of-interest (ROIs) after cortical parcellation based on the Destrieux Atlas (Destrieux, Fischl, Dale, & Halgren, Reference Destrieux, Fischl, Dale and Halgren2010). Cortical thickness and gyrification were averaged across all ROIs to get a mean value per participant.
PRS analyses
Regression models with clinical or neurocognitive variables as dependent variables and PRS as independent variables were used to test for association. Logistic regressions were employed for the binary dependent variables (e.g. body image distortion, binge eating (BE), hoarding, and perfectionism) based on genetic risk (i.e. PRS). For the dependent variables derived by summing across different item variables (e.g., obsession scores), negative binomial regressions were employed to account for over-dispersed count data. For continuous dependent variables (i.e., cognitive traits), linear regression models were run. The regression models were run on each imputed dataset. Each dataset had a different value for the imputed data (non-missing values remain unchanged) because of the random component. As a result, the regression parameter estimates slightly differed between the datasets. The results were then combined to account for the variation in parameter estimates.
Neuroimaging data analyses
Exploratory analyses were run on a subset of the sample with available neuroimaging data (N = 626, 309 females), Mean age = 14.47(Standard Deviation (SD) = 3.64). Linear regression models with averaged cortical parameters as dependent variables and PRS as independent variables tested the association between PRS and cortical features. Moreover, unadjusted and adjusted (for age and parental socio-economic status (SES), calculated as the mean of parents’ education) linear regression models with AN clinical traits as dichotomous independent variables and cortical parameters as continuous dependent variables were used to model the association between AN traits and cortical characteristics.
Regression models, including PRSs, were corrected for the first four principal components. The analyses were first run on the entire sample and then on the sex-stratified datasets. Each type of analysis (e.g., mental health, neurocognitive, cortical) was corrected for multiple comparisons using the Benjamini–Hochberg procedure (False-Discovery Rate (FDR) (Benjamini & Hochberg, Reference Benjamini and Hochberg1995). Given the exploratory nature of our research, we reported both unadjusted and FDR-adjusted analyses. All analyses were run on R version 4.1.1 (R Core Team, 2021). Brain and genetics computations were performed at the University of Geneva using the Baobab HPC service.
Results
Socio-demographic characteristics
The mean age of our sample (2343 females and 2373 males) was 13.81 years (Standard Deviation = 3.67; Range = 14). Table 1 describes participant characteristics.
Demographics and clinical information
Mental health phenotypes
A higher AN PRS was significantly associated with body image distortion (pFDR = 0.02); a one standard deviation increase in the AN PRS corresponded to 1.21 (95% Confidence intervals [CI] = 1.08, 1.35) increased odds of endorsing body image distortion. Additionally, some results passed the threshold of nominal significance (puncorrected < = 0.05): a higher AN PRS associated with compulsive behavior (β = 0.10; 95% CI = 0.01, 0.19; puncorrected = 0.03), a higher SCZ PRS associated with BE (OR = 1.12; 95% CI = 1.01, 1.23; puncorrected = 0.03) as well as perfectionism (OR = 1.10; 95% CI = 1.02, 1.19; puncorrected = 0.02). Moreover, the OCD PRS was associated with psychotic traits (β = 0.07; 95% CI = 0.00, 0.15; puncorrected = 0.05; Table 2).
Associations between PRSs and mental health traits (N = 4,729): Results from logistic/linear regressions
AN: Anorexia Nervosa; BE: Binge Eating; OCD: Obsessive Compulsive Disorder; SCZ: Schizophrenia; FDR: False-Discovery Rate.
* p < = 0.05.
Associations between AN PRS and body image distortion, and SCZ PRS and perfectionism were observed in the female-only dataset (OR = 1.24; 95% CI = 1.08, 1.41; pFDR = 0.04; OR = 1.15; 95% CI = 1.04, 1.28; puncorrected = 0.01, respectively). The association between AN PRS and perfectionism showed nominal significance amongst females (OR = 1.12; 95% CI = 1.00, 1.25; puncorrected = 0.04) (Supplementary Table 3), while a nominally significant association between SCZ PRS and BE was seen in males (OR = 1.21; 95%CI =1.04, 1.40; puncorrected = 0.01) (Supplementary Table 3).
Cognitive traits
No associations passed the threshold of significance (nominal or corrected for multiple comparisons) regardless of biological sex (Table 3).
Associations between PRSs and cognitive measures (N = 4729): Results from linear regression
AN: Anorexia Nervosa; OCD: Obsessive Compulsive Disorder; PRS: Polygenic Risk Score; SCZ: Schizophrenia.
Cortical features
There was a nominal significant association between the AN PRS and CSF volume (β = 4.43; 95% CI =1.40, 7.47; puncorrected < 0.01) (Table 4). Trends in tested associations were observed in the females only dataset: OCD PRS and TIV (β = −17.08; 95%CI = −30.25, −3.90; puncorrected = 0.01) and OCD PRS and GMV, respectively (β = −8.93; 95%CI = −16.72, −1.15; puncorrected = 0.03); and in the males only dataset: AN PRS and TIV (β = 17.55; 95%CI =2.37, 32.73; puncorrected = 0.02), AN PRS and CSF (β = 5.99; 95% CI =1.65, 10.33; puncorrected = 0.01) (Supplementary Table 5).
Associations between PRSs and cortical phenotypes (N = 626): Results from linear regression
AN, anorexia nervosa; CSF, cerebrospinal fluid; GMV, grey matter volume; OCD, obsessive compulsive disorder; PRS, polygenic risk score; SCZ, schizophrenia; TIV, total intracranial volume; WMV, white matter volume.
* p < = 0.05.
In the exploratory analyses focusing on the associations between ED traits (i.e. body image distortion and BE) and brain measures, there was a significant association between body image distortion and reduced GMV (β = −37.56; 95% CI = −60.19, −14.94; pFDR = 0.01) and reduced mean cortical thickness (β = −0.05; 95%CI = −0.08, −0.01; pFDR = 0.04). Body image disturbance was nominally associated with reduced TIV (β = −47.75; 95%CI = −88.73, −6.77; puncorrected = 0.02). BE also showed a nominally significant association with reduced mean gyrification (β = −0.28; 95%CI = −0.50, −0.06; puncorrected = 0.01) (Table 5). Only the association between body image distortion and reduced GMV remained significant after adjustment for SES and age (Table 5). Nominal associations between body image disturbance and GMV and BE and cortical gyrification were observed in the female-only dataset (β = −23.81; 95%CI = −47.32, −0.29; puncorrected = 0.05; β = −0.37; 95%CI = −0.67, −0.07, puncorrected = 0.02, respectively). The results remained significant after correction for SES and age (Supplementary Table 6).
Associations between ED traits and cortical phenotypes (N = 626): Results from unadjusted and adjusteda linear regression analyses
a Adjusted for age and socioeconomic status.
BE: Binge Eating; CSF: Cerebrospinal Fluid; GMV: Grey Matter Volume; TIV: Total Intracranial Volume; WMV: White Matter Volume.
* p < = 0.05.
Amongst males, nominal associations were observed between body image distortion and TIV (β = 71.88; 95%CI = 3.67, 140.09; puncorrected = 0.04) and WMV, respectively (β = 43.88; 95%CI = 9.05, 78.72; puncorrected = 0.01). Only the association between body image distortion and WMV persisted after adjusting for SES and age (Supplementary Table 7).
Discussion
In line with our hypotheses, in this large population-based study, we observed that PRSs for three psychiatric disorders are associated with clinical and cortical characteristics seen in these disorders. The AN PRS was significantly associated with an AN trait in a non-clinical youth population: those with a high AN PRS were at a higher risk of expressing body image distortion. Interestingly, in our sex-stratified analyses, we observed this association amongst females, suggesting that females and males might differ in their genetic propensity for body image distortion. It is also possible that these associations are largely due to sex imbalances (females: 14,898 cases and 27,545 controls vs males: 447 cases and 20,347 controls) in the original GWAS from which we derived our PRS (Watson et al., Reference Watson, Diemer, Zerwas, Gustavson, Knudsen, Torgersen, Reichborn-Kjennerud and Bulik2019). Overall, the findings agree with previous work that showed that AN PRS can significantly inform about AN traits in the general population (Curtis et al., Reference Curtis, Colodro-Conde, Medland, Gordon, Martin, Wade and Cohen-Woods2024; Yilmaz et al., Reference Yilmaz, Schaumberg, Halvorsen, Goodman, Brosof, Crowley, Mathews, Mattheisen, Breen, Bulik, Micali and Zerwas2023).
The results also confirmed our second hypothesis that disease-specific PRS can better inform about disease-specific phenotypes than related PRS. For example, AN PRS informed better than any other psychiatric PRS about AN-specific symptomatology (i.e. body image distortion). Nevertheless, we also observed trends in associations that cross diagnostic boundaries. For example, SCZ PRS associated with BE and perfectionism, and AN PRS associated with compulsive behavior. These findings are in line with previous findings (Lundgren, Rempfer, Brown, Goetz, & Hamera, Reference Lundgren, Rempfer, Brown, Goetz and Hamera2010; Ramacciotti et al., Reference Ramacciotti, Paoli, Catena, Ciapparelli, Dell’Osso, Schulte and Garfinkel2004; Solmi, Mascarell, Zammit, Kirkbride, & Lewis, Reference Solmi, Mascarell, Zammit, Kirkbride and Lewis2019b) and add to the literature that supports an association between EDs and SCZ (Bulik-Sullivan et al., Reference Bulik-Sullivan, Finucane, Anttila, Gusev, Day, Loh, Consortium, Consortium, Duncan, Perry, Patterson, Robinson, Daly, Price and Neale2015; Duncan et al., Reference Duncan, Yilmaz, Walters, Goldstein, Anttila, Bulik-Sullivan, Ripke, Thornton, Hinney, Daly, Sullivan, Zeggini, Breen and Bulik2017; Zhang et al., Reference Zhang, Larsen, Kuja-Halkola, Thornton, Yao, Larsson, Lichtenstein, Petersen, Bulik and Bergen2021). We observed some sex-specific trends in PRS associations, with the AN and SCZ PRSs showing marginal evidence of an association with perfectionism amongst females and the SCZ PRS nominally associating with BE in males. These cross-disorder associations might result from spurious associations when testing multiple outcomes or might point to higher-order psychopathological factors like factor p. While nominally significant, these associations require further validation for a conclusive statement on the transdiagnostic effects of genetic risk.
Our study failed to replicate previous associations between psychiatric PRS and neurocognitive traits (Hatzimanolis et al., Reference Hatzimanolis, Bhatnagar, Moes, Wang, Roussos, Bitsios, Stefanis, Pulver, Arking, Smyrnis, Stefanis and Avramopoulos2015; Mallet, Strat, Dubertret, & Gorwood, Reference Mallet, Strat, Dubertret and Gorwood2020). One possible explanation for the lack of results is that the cognitive alterations seen in clinical populations, on which the PRSs are computed, are extreme and might not be present at the population level, or they might be state markers of clinical diagnoses. Another possible explanation is that the chosen neurocognitive traits might not be good trait markers for the disorders. One study looking at the most suitable cognitive markers for SCZ identified a measure of sustained attention and a measure of general cognitive function as the best candidates, based on the shared genetic effect with SCZ (McCarthy et al., Reference McCarthy, Badcock, Clark, Knowles, Cadby, Melton, Morgan, Blangero, Moses, Glahn and Jablensky2017).
The lack of association between AN PRS and global cortical volume is in line with recent reports showing a nominal association between AN PRS and reduced caudate volume but no associations with global cortical measures (e.g. cortical thickness, surface area) (Westwater et al., Reference Westwater, Mallard, Warrier, Bethlehem, Scheinost, Grillon, Fletcher, Seidlitz and Ernst2023). While we did not observe any significant associations between PRSs and cortical features, we did observe a significant association between body image distortion and GMV in the whole group and in the female-only group. The lack of significant associations between the included PRSs and cortical characteristics, in the presence of an association between mental health traits and cortical characteristics (specifically for AN), suggests that brain alterations might be state markers and not trait markers. This again supports previous reports of restoration of cortical thinning to a normal level in AN patients following weight gain, further suggesting that brain alterations follow the clinical diagnosis and are not neural traits of the disorder (Bernardoni et al., Reference Bernardoni, King, Geisler, Stein, Jaite, Nätsch, Tam, Boehm, Seidel, Roessner and Ehrlich2016; King et al., Reference King, Geisler, Ritschel, Boehm, Seidel, Roschinski, Soltwedel, Zwipp, Pfuhl, Marxen, Roessner and Ehrlich2015; Nickel et al., Reference Nickel, Joos, Tebartz van Elst, Matthis, Holovics, Endres, Zeeck, Hartmann, Tüscher and Maier2018).
Our analysis failed to replicate previous reports of associations between PRS SCZ and impaired cognitive functioning (Germine et al., Reference Germine, Robinson, Smoller, Calkins, Moore, Hakonarson, Daly, Lee, Holmes, Buckner, Gur and Gur2016; Hagenaars et al., Reference Hagenaars, Harris, Davies, Hill, Liewald, Ritchie, Marioni, Fawns-Ritchie, Cullen, Malik, Worrall, Sudlow, Wardlaw, Gallacher, Pell, McIntosh, Smith, Gale and Deary2016; Liebers et al., Reference Liebers, Pirooznia, Seiffudin, Musliner, Zandi and Goes2016) and between OCD PRS and brain region functioning (Heinzel et al., Reference Heinzel, Kaufmann, Grützmann, Klawohn, Riesel, Bey, Heilmann-Heimbach, Weinhold, Ramirez, Wagner and Kathmann2021). The lack of associations between OCD and SCZ PRS and cognitive or cortical traits can be seen in the light of our methodology, for example, the choice of cognitive measures and focus on structural brain morphometrics. To the best of our knowledge, no other investigation, looked at the association between genetic risk for OCD and SCZ, indexed as PRSs, and global cognitive and cortical measures in a sample of young adults.
Strengths and limitations
The present study provides a comprehensive view of the relationship between three psychiatric disorders PRSs and several phenotypes and neural characteristics (i.e. clinical, neurocognitive, and neurocortical) in a non-clinical population during an important time window for psychiatric disorder development.
The current study has many methodological strengths that are worth mentioning. First, the approaches employed in the data analysis aimed to maximize data availability (e.g. imputation of missing data). Second, our analyses were data-driven, without stringent a priori hypotheses. This allowed the exploration of all the available data, without bias. This approach allowed the expansion of the current knowledge outside the previously tested hypotheses. Additionally, a tailored brain template was used, adequate for our sample (i.e., age and gender), for our neuroimaging analysis. This allowed to capture of cortical changes linked to development that might have otherwise been missed.
Several limitations need to be acknowledged. First, despite relying on the K-SADS, a validated instrument assessing psychiatric disorders, the assessment for ED symptoms was fairly short (i.e., 2 clinical items). Moreover, we were unable to incorporate a measure of body mass index (BMI) into the computation of a variable that can better reflect the AN symptomatology, as we were limited by the availability of anthropometric data (i.e., 20% missing for height and 19% missing for weight). Second, we were limited in our power in the males-only dataset by the unbalance of AN symptoms between females and males (i.e., low prevalence of AN symptoms in males) (Supplementary Table 1). Third, the AN PRS, a measure derived mainly from females with an AN diagnosis (Watson, Diemer, et al., Reference Watson, Diemer, Zerwas, Gustavson, Knudsen, Torgersen, Reichborn-Kjennerud and Bulik2019), was tested against male-reported AN symptoms, which might explain our inability to detect any significant associations between AN PRS and AN traits in the male sample. Additionally, while we did try to account for the age heterogeneity in the neuroimaging analyses, we acknowledge that the sample age range spans important developmental time points (pre-teens, teens, and young youth) that are worth investigating in greater detail. Additionally, the cross-sectional design of the study limits our ability to make inferences about state vs trait in the true sense of the word. While we do not have any information on the temporality of behaviours (i.e., when a certain behavior occurs), we assumed that since genetic predisposition remains constant throughout a lifetime that the absence of an association between PRS and any phenotype could suggest that the specific characteristic (cognitive, cortical, or clinical) is a consequence of the illness itself. Finally, we limited our analysis to individuals of European descent in the absence of a genome-wide association study (GWAS) of our traits for individuals of African descent and PRS methods that deal well with admixed individuals. However, we showed in our sensitivity analysis (Supplementary Table 9) that PRS SCZ performed similarly in the two ancestral groups (European and African), suggesting that not including individuals of African descent into our main analysis did not significantly alter the results and their interpretation.
It is likely that our ability to link mental health traits to genetic propensity will be improved by a better understanding of the genetics of AN, OCD, and SCZ that will come with larger, more ancestrally diverse GWAS. To our knowledge, a new GWAS of AN and binge eating disorder, that should advance our understanding of the underlying biology of these disorders, is underway. Future studies should aim to replicate the present findings in larger longitudinal cohorts, as the results of the current work might have been diluted by the heterogeneity introduced by the large age range and by the sample size. Given the observed differences in genetic risk between males and females, future research should consider stratifying the analysis by biological sex. This would address previously expressed concerns that both sex and developmental timing influence the developmental pathways to psychiatric risk (Yilmaz et al., Reference Yilmaz, Schaumberg, Halvorsen, Goodman, Brosof, Crowley, Mathews, Mattheisen, Breen, Bulik, Micali and Zerwas2023) and that genetic risk for AN might have different effects on cortical measures in men and women (Leehr et al., Reference Leehr, Opel, Werner, Redlich, Repple, Grotegerd, Dohm, Goltermann, Böhnlein, Amare, Sindermann, Förster, Meinert, Enneking, Richter, Hahn, Baune and Dannlowski2019).
Lastly, future neuroimaging studies investigating the relation between PRS and cortical parameters should use all the available information. Given that the analyzed data originated from a non-clinical population and no a priori hypotheses regarding specific brain regions were formulated, our neuroimaging analyses only looked at global and averaged cortical measures such as TIV, cortical thickness, and GMV. However, it is well established that AN patients show brain changes in specific regions involved in eating-related functions such as the insula and the orbitofrontal cortex (Frank, Reference Frank2015). Further studies should also focus on ROIs to better understand what role they play in these disorders. Both approaches (e.g., brain region specific and whole brain) are needed and complementary as they can advance our understanding of the effect of psychopathology at the macro and micro scale.
In conclusion, our study shows that genetic risk for AN, quantified through a PRS, can inform about AN traits, while genetic risk for AN, OCD, and SCZ did not associate with cognitive or cortical markers. Furthermore, our results suggest that brain alterations (i.e., reduced GMV and reduced cortical thickness) associated with AN traits are likely to be the effect of illness (i.e., state markers) and not an inherent trait (i.e., trait markers). This study sheds light on the shared genetic risk between AN, OCD, and SCZ and potential transdiagnostic traits.
Supplementary material
The supplementary material for this article can be found at http://doi.org/10.1017/S0033291725101839.
Funding statement
This work was partly supported by a Laureate Grant Award from the Novo Nordisk Foundation (Grant No. NNF22OC0071010, to N.M.).
Open access
