Genetic liability for schizophrenia and childhood psychopathology in the general population

Background Genetic liability for schizophrenia is associated with psychopathology in early life. It is not clear if these associations are time-dependent during childhood, nor if they are specific across different forms of psychopathology. Methods Using genotype and questionnaire data on children (N = 15,105) from the Norwegian Mother, Father, and Child Cohort Study (MoBa), we tested associations between schizophrenia polygenic risk scores and measures of childhood emotional and behavioural problems for developmental stability and domain specificity. We then sought to identify symptom profiles, across development and domains, associated with elevated schizophrenia polygenic liability. Results. We found evidence for developmental stability in associations between schizophrenia polygenic risk scores and emotional and behavioural problems, with the latter being mediated via the rate of change in symptoms between 18 months and 5 years specifically (B = 0.032; 95% CI 0.007-0.057). At age 8, associations with emotional and behavioural psychopathology were found to be better explained by a model of symptom-specific polygenic risk score effects, rather than effects mediated via a general p factor or by domain-specific factors. Overall, individuals with higher schizophrenia polygenic risk scores were more likely (OR= 1.310 [95% CIs: 1.122-1.528]) to have increasing behavioural and emotional symptoms in early childhood, followed by relatively elevated symptoms of conduct disorder, oppositional defiant disorder, hyperactivity and inattention in middle childhood. Interpretation Schizophrenia-associated alleles are linked to specific patterns of early-life psychopathology. The associations are small, but findings of this nature can help us better understand the developmental emergence of schizophrenia.


Research in context
Evidence before this study Genetic liability to schizophrenia, conferred incrementally by many genetic variants of small effect, is associated with various forms of psychopathology -both symptoms and diagnosesin the general population during adulthood. To get an overview of the relevant evidence for how genetic liability to schizophrenia manifests in psychopathology in childhood, we performed a systematic search of the published literature using the Ovid MEDLINE and PsychINFO databases, for English-language peer-reviewed journal articles published prior to 28 January 2020. We found 7 studies of core relevance (i.e., containing assessment of symptoms or diagnoses of psychopathology in pre-adolescent children), with a further 9 studies examining outcomes potentially related to psychopathology (such as brain structure, cognitive performance, and social functioning). Of the 7 core studies, 3 used clinical samples to demonstrate that polygenic risk scores for schizophrenia are higher among children with ADHD diagnoses than controls, and among cases of the rare child-onset form of schizophrenia than their healthy siblings. The remainder of studies all examined symptoms of psychopathology among children in the general population, finding modest but robust associations of schizophrenia polygenic risk scores with emotional and behavioural problems measured from 3 years of age, as well as with symptoms of depression, ADHD, anxiety, oppositional defiant disorder and conduct disorder in middle childhood.

Added value of this study
In this study, we present a set of analyses designed to improve our understanding of the nature of associations between schizophrenia risk alleles and childhood psychopathology. Specifically, we employ an approach that aims not just to quantify, but also to explore how the effects of schizophrenia risk manifest across childhood, and across different domains of psychopathology. We find evidence that effects of schizophrenia polygenic risk scores on symptoms of emotional and behavioural problems in early childhood are stable, influencing the overall level and rates of change in symptoms, rather than age-specific (i.e., transient or developmental). We also find evidence of specificity in the effects of schizophrenia polygenic risk scores on different domains of psychopathology in 8-year-old children. Overall, we find that higher schizophrenia polygenic risk scores are associated with a developmental symptom profile comprising elevated and increasing symptoms of behavioural problems and increasing levels of emotional problems in early childhood, as well as particularly elevated symptoms of conduct disorder, inattention, hyperactivity, and oppositional defiant disorder in middle childhood.
Implications of all the available evidence Findings of our study align with a growing body of evidence that the effects of schizophrenia risk alleles on psychopathology begin early in life, and influence the likelihood of children experiencing difficulties across development. While previous work has largely found similar effects of schizophrenia polygenic risk scores across different domains of childhood psychopathology, indicating that such effects may be mediated by a hypothetical latent 'general psychopathology' or 'p' factor, our results suggest that domain-and even symptomlevel specificity may emerge by middle childhood. We may be able to improve our understanding of processes underpinning the emergence of schizophrenia later in life by paying attention to nuances in the ways that genetic risk for schizophrenia manifests across childhood and into adolescence.

Methods
Using genotype and questionnaire data on children (N = 15,105) from the Norwegian Mother, Father, and Child Cohort Study (MoBa), we tested associations between schizophrenia polygenic risk scores and measures of childhood emotional and behavioural problems for developmental stability and domain specificity. We then sought to identify symptom profiles -across development and domains -associated with elevated schizophrenia polygenic liability.

Outcomes
We found evidence for developmental stability in associations between schizophrenia polygenic risk scores and emotional and behavioural problems, with the latter being mediated via the rate of change in symptoms between 18 months and 5 years specifically (β slope = 0.032; 95% CI 0.007 -0.057). At age 8, associations with emotional and behavioural psychopathology were found to be better explained by a model of symptom-specific polygenic risk score effects, rather than effects mediated via a general "p" factor or by domain-specific factors. Overall, individuals with higher schizophrenia polygenic risk scores were more likely (OR=

Introduction
Genetic risk for schizophrenia is highly polygenic (1,2), representing an additive combination of many common genetic variants with small effects. Large-scale genome-wide association studies (GWAS) (3) discover these effects, with currently more than 176 genetic loci identified as conferring risk for schizophrenia (4). Effect sizes and standard errors from GWAS can be combined across many variants in a polygenic risk score (PRS) (5), to examine how common genetic liability for schizophrenia is associated with traits and behaviours in the general population (6).
Schizophrenia is most commonly diagnosed in late adolescence or early adulthood (7). PRS-based approaches can be used to study early-life manifestations of genetic risk for schizophrenia, with the potential to inform as to how and why the disorder emerges later in life. Previous studies have shown that genetic liability for schizophrenia is modestly associated with a range of childhood outcomes, including infant neuromotor development (8), early neurocognitive and behavioural development (9), sleep problems (10) and social cognition (11) -as well as with measures of psychopathology (including symptoms of anxiety, depression, attention deficit hyperactivity disorder, and conduct problems) across childhood (12)(13)(14)(15). These associations appear to persist into adolescence (16,17) as well as potentially diversifying further (e.g., into disordered eating (18) and cannabis use (19)).
Although it seems clear that early life manifestations of schizophrenia genetic liability across the psychopathological "phenome" are relatively diverse, it is not well understood how this diversity arises. For example, it could be that schizophrenia-associated genes have transient, developmentally-varying effects; or that early effects in specific domains trigger developmental cascades that encompass a wider range of behaviours. It could simply be the case that these genes have highly generalized effects on behaviour; or that the environmental and maturational restrictions of childhood make their expression more diffuse than is observed in adulthood. To narrow down the various possible explanations, it is necessary to investigate associations in two additional dimensions: developmental time, and phenotypic space. That is, to ask to what degree effects are developmentally stable (influencing behaviour consistently across development) vs. age-specific (transient or emerging at a specific point), and to explore how broadly vs. selectively they influence different behaviours and symptoms.
In the current study, we use data from the population-based pregnancy cohort: the Norwegian Mother, Father, and Child Cohort Study (MoBa) (20) to investigate the extent to which manifestations of genetic liability for schizophrenia in childhood emotional and behavioural problems are: i) developmentally stable vs. unstable (age-specific); and ii) domain-general vs. domainspecific. Additionally, we investigate genomic prediction, using schizophrenia PRS, of individuals' latent symptom profiles across development and a broad range of symptom domains.
The analyses we present are exploratory, with the aim of identifying the most parsimonious of a set of models describing the covariation between schizophrenia PRS and questionnaire measures of childhood emotional and behavioural problems. However, we also had two broad hypotheses. First, we expected developmentally-stable associations between schizophrenia-associated genetic variants and symptoms of psychopathology across early and middle-childhood, based on both work with schizophrenia PRS (12,13,15) and findings of the stable genetic influence on childhood psychopathology (21). Second, we expected associations later in childhood to be primarily general, mediated via a latent, general "p" (for "psychopathology") factor (22). The "p" factor model offers a parsimonious explanation for comorbidity and shared genetic influence between symptom domains (23,24) and schizophrenia PRS have been shown to explain variance in the latent "p" component of the model in childhood (25) and adolescence (26).

Study sample
The Norwegian Mother, Father and Child Cohort Study (MoBa) is a population-based pregnancy cohort study conducted by the Norwegian Institute of Public Health (27). Participants were recruited from all hospitals and obstetric units in Norway during 1999-2008. The women consented to participation in 41% of the pregnancies. The cohort now includes 114,500 children, 95,200 mothers and 75,200 fathers. The current study is based on version 12 of the quality-assured data files released in January 2019.
The establishment and data collection in MoBa were previously based on a license from the Norwegian Data Protection Agency and approval from The Regional Committee for Medical Research Ethics (REK) and is now based on regulations under the Norwegian Health Registry Act. The current study was approved by REK (2016/1702).
Blood samples were obtained from children (umbilical cord) at birth. Genotyping of the entire MoBa cohort is ongoing. We used genotype data from children in 17,000 randomly selected trios, and after quality control the analytic sample was a genotyped sub-set (N = 15,105) of children (20). Details about the processing of the genetic data are outlined in the Appendix.

Measures
We used measurements of children's symptoms of emotional and behavioural psychopathology from maternal questionnaires collected when children were aged 18 months, 3 years, 5 years, and 8 years. Fifteen items from the Child Behaviour Checklist (CBCL) (28), which could be sub-divided into emotional problems (5 items) and behavioural problems (10 items), were included for longitudinal analyses across early childhood (18 months to 5 years). Three instruments measuring six domains of emotional and behavioural psychopathology in middle childhood (at the 8-year data collection) were included for the specificity analyses. The 13-item Short Mood and Feelings Questionnaire (sMFQ) (29) measured symptoms of depression, a 5-item short form of the Screen for Child Anxiety Related Disorders (SCARED) (30) measured symptoms of anxiety, and a 34-item version of the Rating Scale for Disruptive Behaviour Disorders (RS-DBD (31)) measured symptoms of conduct problems (CD), oppositional defiant disorder (ODD), hyperactivity, and inattention. Details about measure selection and psychometrics are in the appendix.

Polygenic risk scores
We calculated PRS for schizophrenia using PRSice2 (32), based on European samples from the most recent (2) Psychiatric Genomics Consortium schizophrenia genome-wide association study (GWAS). PRS can be calculated using effect estimates for all variants in common between the discovery (i.e., GWAS) sample and target sample, for variants whose effects in the GWAS had a pvalue below a specified threshold. Typically, PRS are created at a range of p-value thresholds (0-1), reflecting the expectation that the polygenic signal will gradually increase as more variants (either with weaker effects or a lower frequency in the population) are included, up until the point additional variants contribute only statistical noise to the score. We opted to use the thresholds from the polygenic scoring analysis presented in the original GWAS (2). These were: p < 0.001, p < 0.01, p < 0.05, p < 0.1, p < 0.2, p < 0.5, p < 1. Full details of the parameters and procedure used in generating the PRS are in the Appendix.

Analyses
Developmental modelling . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted May 5, 2020.
Latent growth models were used to explore variability in the development of emotional and behavioural problems between the ages of 18 months and 5 years in the context of genetic liability for schizophrenia. In these models, variance in observed variables is explained by two latent growth factors: a latent intercept, which is specified to load equally at all waves (here at 1.5, 3, and 5 years), and a latent slope, loadings for which are fixed proportional to their temporal distance from the first wave of measurement (here the loadings are 0, 1.5, and 3). For each schizophrenia PRS (based on the different thresholds), we formally compared the fit of models specifying an effect on growth parameters (intercept and slope) vs age-specific residuals in emotional and behavioural problems. This model fitting indicated to what extent the association between schizophrenia PRS and emotional/behavioural problems is developmentally stable (i.e., via latent growth factors), agespecific (i.e., via residuals), or null. Informal comparisons were also made to ascertain whether stable effects could be primarily ascribed to either the latent intercept or slope factor (for full details of model fitting and comparison strategy, see appendix).
"p" factor modelling Next, we specified models describing covariation among items from the 8-year psychopathology variables with 6 domain-specific factors and one overall "p" factor. We compared their fit to the data when schizophrenia PRS were allowed to influence, respectively: symptom-specific residuals, domain-specific factors, the domain-general "p" factor, as well specified to have no effect. All versions of a given model were nested, so formal tests of fit (χ 2 difference tests) were used to ascertain the best-fitting model. For both the developmental and "p" factor modelling, we examined the average R 2 value among the outcomes and present results for the PRS threshold at which R 2 was maximised.

Latent profile analysis
Finally, we incorporated both the developmental modelling and "p" factor modelling approaches into a latent profile analysis. In latent profile analysis, a categorical latent variable is used to assign individuals in a sample into one of a pre-specified number of profiles based on their pattern of scores across multiple observed continuous variables. For our analysis, profile membership was informed by individuals' estimated scores on the continuous latent growth factors from the emotional and behavioural problems latent growth models, and their observed scores on compiled scale versions of the 6 psychopathology domains measured at 8 years. We then included schizophrenia PRS as predictors of profile assignment. The full model is displayed in Figure 1. We specified models with 2, 3, 4, 5, and 6 profiles respectively, and assessed them using a combination of standard criteria (Vuong-Lo-Mendell-Rubin likelihood ratio test (33,34), entropy, fit indices).
All modelling for the first two parts of the analyses was carried out in R version 3.4.4 using the lavaan package (35)

Results
Descriptive statistics for all study variables are available in Table 1, alongside an exploration of the extent of selection bias among our genotyped sub-sample of MoBa and of selective attrition with respect to the longitudinal analyses (both of which were observed at low levels) is included in the Appendix.
Schizophrenia genetic liability and developmental stability in childhood psychopathology . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted May 5, 2020. Linear growth models provided an acceptable fit to the data from the measures of emotional and behavioural problems (for details, see appendix). We found evidence of associations between schizophrenia genetic liability and emotional and behavioural problems, with the null model rejected at most p-value thresholds in each domain. Moreover, we observed evidence in favour of developmentally stable rather than age-specific associations, with models incorporating PRS effects on the latent growth factors preferred in both domains (see appendix for genetic model-fitting details). For emotional problems, we could not formally distinguish whether these effects were primarily mediated via the intercept (i.e., predicting the overall level of symptoms) or slope (i.e., rate of change in symptoms over time) in the best-fitting model, and point estimates were similar for both (β intercept = 0.019, 95% CI: -0.005-0.043; β slope = 0.019, 95% CI: -0.010-0.048; PRS threshold: p< 0.1). In contrast, for behavioural problems, these effects were mediated via the slope factor alone at the most predictive threshold (β slope = 0.032; 95% CI 0.007-0.057; threshold: p<1), meaning that schizophrenia PRS predicted rate of change in symptoms of behavioural problems across early childhood. However, in both cases the associations were small, explaining ~0.1% variance in the growth factors. Parameter estimates and 95% CIs for best-fitting models are displayed in the Appendix.

Schizophrenia genetic liability and phenotypic specificity in childhood psychopathology
A structural model of domain-specific factors and one overarching "p" factor provided a good fit to item-level data from the 8-year psychopathology measures (CFI=0.98, TLI=0.97, RMSEA=0.04). Details of the model fitting process are presented in the appendix. We found that schizophrenia PRS were significant predictors of variability in symptoms of psychopathology at 8 years, and that this prediction was maximised at the PRS threshold including all variants (p<1). At this threshold, neither the model incorporating a PRS effect on the general "p" factor nor the model incorporating domain-specific effects fit the data sufficiently well to be accepted. Instead, the preferred model allowed the PRS to influence item-specific residuals -i.e., explaining symptom-specific variation. Figure 2 shows standardized beta coefficients for these symptom-level relationships.
The pattern of results in Figure 2 show that there is substantial heterogeneity, in terms of schizophrenia PRS associations, between symptoms -even within domains. This is particularly evident for symptoms within the depression and inattention domains. As with the developmental modelling, effect sizes were very small, with PRS explaining a maximum of 0.4% variance in symptoms. Parameter estimates and 95% CIs for the alternate versions of the model, which were preferred for PRS comprising variants significant only at less stringent p-value thresholds, are presented in the appendix.

Schizophrenia genetic liability predicting developmental symptom profiles
Scores from all scales (CBCL measures of emotional and behavioural problems at 18 months, 3 years, and 5 years and six 8-year psychopathology domains) were incorporated into a single model for the latent profile analysis. Fit statistics for these models are presented in the appendix. Entropy (reflecting the overall certainty with which individuals could be assigned to symptom profiles) was good (~0.8) for each version of the model, but the Vuong-Lo-Mendell-Rubin test statistic for the 5profile model indicated that it offered no significant improvement of fit on the simpler 4-profile model. Figure 3 shows the 4 symptom profiles. Individuals assigned to profile 1 (5.7% of sample) had increasing symptoms of emotional problems between 18 months and 5 years and elevated anxiety symptoms at age 8. Profile 2 (7.9%) was characterised by moderate, stable behavioural problems symptomatology in early childhood, and moderate symptoms of ODD, hyperactivity, and inattention at 8 years. Profile 3 was considered the normative symptom profile, as a large majority (84.7%) could be assigned to it, and it was characterised by decreasing problems developmentally and low levels of symptoms at 8 years. Profile 4 was the least populous (1.8% of the sample) and most symptomatic profile, characterised by increasing/elevated symptoms of behavioural problems and . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted May 5, 2020. increasing symptoms of emotional problems preceding relatively elevated symptoms of conduct disorder, ODD, hyperactivity and inattention at 8 years.
The probabilities of classification into profiles 1, 2, and 4 relative to the normative profile are shown in Figure 4 as a function of PRS (at the p <0.05 threshold), alongside odds ratios of PRS on probability of classification into a given profile rather than the normative profile. Profile 4 was better defined than profiles 1 or 2 (the density of high probabilities of assignment into this profile in the bottom-left panel of Figure 4 is greater than in the two panels above, and there is greater separation between individuals classified in this profile and those with normative profile). The average PRS was higher in profile 4 than in all other profiles, and indeed individuals with higher PRS were more likely to be assigned to profile 4 than the normative profile (OR= 1.310 [95% CIs: 1.122-1.528], bottom-right of figure). There was no strong evidence of associations between schizophrenia PRS and probability of assignment to either of the other profiles (relative to the normative profile).

Discussion
The results indicate that associations between schizophrenia risk-associated alleles and symptoms of early childhood (18 months to 5 years) psychopathology tend to be developmentally stable, rather than transient or emergent. In middle childhood (8 years) we found that, contrary to expectations, associations between schizophrenia genetic risk and symptoms of psychopathology may be better explained by symptom-specific effects than by prediction of a latent "p" factor indexing general psychopathology. Finally, we showed that by combining information about symptoms of psychopathology across time and domains, associations with schizophrenia genetic liability can be summarised in the form of a characteristic symptom profile.
Our results are consistent with findings of widespread associations between schizophrenia genetic liability and symptoms of childhood psychopathology. There are a number of possible explanations for such associations. They may arise because of pleiotropy; that is, direct effects of schizophreniarelated genetic variants on other behaviours and symptoms. They may represent early manifestations of schizophrenia liability, where environmental or maturational restrictions mean that the classic "schizophrenia phenotype" (and its sub-clinical analogues) cannot be fully expressed until adolescence or later. A further alternative is that (some) childhood manifestations of genetic liability for schizophrenia are actually part of a causal process in the development of the disorder. This could involve intrinsic, developmental processes (e.g., specific thought patterns becoming "grooved" over time) or extrinsic causal processes, where genetic risk is mediated via the environment. An example of such a process would be the putatively causal link between cannabis use and schizophrenia (38,39). Cannabis is an environmental exposure, but its use is associated with genetic variants linked to impulsivity (40) and risk-taking (41). Assuming a demonstrable causal relationship between cannabis use and the development of schizophrenia, the environment (cannabis use) would mediate the effects of these genetic variants on schizophrenia; and impulsivity and risk-taking earlier in development would be on the causal pathway.
We applied structural models to associations between schizophrenia PRS and symptoms of childhood psychopathology to provide context that might render some of these competing explanations more or less likely. Using developmental models, we found that associations between schizophrenia PRS and emotional and behavioural symptomatology were best explained via effects on developmentally stable growth processes. This is in line with evidence on the stability of genetic influences, in general, across childhood (21). We had limited power to investigate the nature of these stable associations, but found evidence that, for behavioural problems, schizophreniaassociated alleles are primarily predictive of individual differences in their propensity to change over time in early childhood. Such a pattern may be more consistent with an aetiological model implicating environmentally mediated genetic effects and active or evocative gene-environment correlations, than direct pleiotropy. Future work incorporating measured environmental mediators may be able to establish whether these mechanisms are involved in the childhood manifestation of schizophrenia genetic liability.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted May 5, 2020. The "p" factor structure of childhood psychopathology has been largely supported by previous genetically-informative work indicating that genetic prediction of symptoms is primarily mediated via a domain-general factor (e.g., (23)). Schizophrenia genetic risk has been shown to associate with psychopathology via this route (14,16). In apparent contrast, we found that models with symptomspecific effects were preferred, despite the cost of estimating many more parameters. There is likely some overfitting involved in these models, which is why we do not interpret the individual symptomlevel effects. Nonetheless, the heterogeneity of these effects -even within domains -is clear and strengthens the conclusion that genetic liability for schizophrenia may influence childhood psychopathology in a more specific and less generalized manner than has been shown previously.
It is likely that two aspects of our approach explain these differences. First, we used item-level data, providing more scope for heterogeneity -from both signal and noise -to emerge, whereas scalelevel data has often been used in the past. Second, we tested schizophrenia PRS at different thresholds -including thresholds at which all or almost all SNPs are included in the score. Previous work has often selected a single threshold score, typically p< 0.05 or lower. We found that lower threshold scores were likely to have smaller associations that, partly for reasons of statistical power, led to more parsimonious models being preferred. However, the amount of variance explained in symptoms continued to increase for scores created with less stringent thresholds suggesting that a single-score approach to using PRS potentially leads to meaningful signal from SNPs being missed. Using scores at higher thresholds both increases power and is consistent with the theory underpinning the polygenic model of genetic effects (5).
The small effect sizes for schizophrenia PRS on measures of childhood psychopathology in our sample are consistent with the literature (9,12,15). Our latent profile analytic approach demonstrates one way in which combining information across time and symptom domains can draw out more meaningful effect sizes (i.e., we observed effects equivalent to a 30% increase in the probability of displaying a particular symptom profile per standard deviation increase in schizophrenia PRS). However, we acknowledge the inherent bounds on the utility of effects explaining so little variance in childhood and adolescent outcomes. Using PRS for adult psychiatric disorders to identify children and adolescents at substantially increased, personal risk for behavioural and emotional difficulties and basing targeted prevention strategies on this information is currently not a realistic prospect. If the true effects captured by PRS are of the magnitude that the evidence so far suggests, PRS will have more value as tools to help clarify the mechanisms by which symptoms emerge, are maintained, differentiate, exacerbate or improve developmentally, than as individual-level predictors of disorder risk.
Our study is subject to several limitations. Small effect sizes for schizophrenia PRS on childhood traits mean that, even with a sample of several thousand children, we have limited power for differentiating similar parameterisations (e.g., on latent growth processes) and therefore opted against stratifying by sex. Secondly, the genotyped sample we use is subject to some selection effects, and the MoBa sample as a whole is (like all cohort studies) affected by selective attrition as the propensity to drop out of studies over time is linked to poorer health, which is also reflected genetically (42). Furthermore, the measures used to index symptoms of emotional and behavioural psychopathology are relatively brief and available only as maternal reports. Corroboration by additional reporters and using clinical interviews or diagnoses from linked health care registries would help to mitigate this limitation in future work.
Overall, our findings suggest that schizophrenia-related genetic variants are associated with symptoms of psychopathology stably from as early as 18 months of age, and with considerable phenotypic specificity by middle childhood. Further work -such as expanding models longitudinally and to incorporate measured environments and testing within-family polygenic prediction -can help to refine and triangulate developmental explanations for the emergence of schizophrenia later in life.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted May 5, 2020. is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted May 5, 2020.  is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted May 5, 2020. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted May 5, 2020. .

Figure 1. Latent profile analysis to ascertain developmental and domain-specific profiles associated of psychopathology symptoms associated with schizophrenia genetic effects
Note-boxes represent observed variables and circles model-estimated latent variables; I = intercept factor, which loads equally on observed variables at all waves; C = categorical latent variable, sub-dividing the sample into a specified number of classes (we tested models with 2, 3, 4, 5, and 6)  . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted May 5, 2020. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted May 5, 2020.  is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted May 5, 2020.

Figure 4. Density of relative probabilities of classification into the different profiles by schizophrenia PRS and odds ratios for prediction of symptom profile classification by PRS
Note-contour plots show probabilities of classification into a specific profile (1, 2, or 4) relative to probabilities of classification into the normative profile (3) for individuals ultimately assigned to each of these profiles (indicated by colour coding), as a function of PRS; diamond marker indicates within-profile means (both classification probability and PRS) and the shaded region around them shows the 95% CIs in either dimension (i.e., vertically for probability, horizontally for PRS); profile 3 is used as reference category for logistic regression; PRS threshold used is p<0.05, selected as it maximizes the OR for profile 4 vs. profile 3; but the pattern of results is consistent at other thresholds (see Appendix) . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted May 5, 2020. . https://doi.org/10.1101/2020.05.01.20086215 doi: medRxiv preprint