An axis of genetic heterogeneity in autism is indexed by age at diagnosis and is associated with varying developmental and mental health profiles ================================================================================================================================================= * Xinhe Zhang * Jakob Grove * Yuanjun Gu * Cornelia K. Buus * Lea K. Nielsen * Sharon A.S. Neufeld * Mahmoud Koko * Daniel S Malawsky * Emma Wade * Ellen Verhoef * Anna Gui * Laura Hegemann * APEX consortium * iPSYCH Autism Consortium * PGC-PTSD Consortium * Daniel H. Geschwind * Naomi R. Wray * Alexandra Havdahl * Angelica Ronald * Beate St. Pourcain * Elise B. Robinson * Thomas Bourgeron * Simon Baron-Cohen * Anders D. Børglum * Hilary C. Martin * Varun Warrier ## Abstract There is growing recognition that earliest signs of autism need not clearly manifest in the first three years of life. To what extent is this variation in developmental trajectories associated with age at autism diagnosis? Does the genetic profile of autism vary with age at autism diagnosis? Using longitudinal data from four birth cohorts, we demonstrate that two different trajectories of socio-emotional behaviours are associated with age at diagnosis. We further demonstrate that the age at autism diagnosis is partly heritable (h2SNP = 0.12, s.e.m = 0.01), and is associated with two moderately correlated (rg = 0.38, s.e.m = 0.07) autism polygenic factors. One of these factors is associated with earlier diagnosis of autism, lower social and communication abilities in early childhood. The second factor is associated with later autism diagnosis, increased socio-emotional difficulties in adolescence, and has moderate to high positive genetic correlations with Attention-Deficit/Hyperactivity Disorder, mental health conditions, and trauma. Overall, our research identifies an axis of heterogeneity in autism, indexed by age at diagnosis, which partly explains heterogeneity in autism and the profiles of co-occurring neurodevelopmental and mental health profiles. Our findings have important implications for how we conceptualise autism and provide one model to explain some of the diversity within autism. ## Main Autism is a term used to describe a group of conditions characterised by difficulties in social-communication, unusually restricted and repetitive interests, and sensory differences1. Ever since its earliest descriptions in the 1940s2,3, autism has been thought of as a condition that typically emerges and is diagnosed in early childhood. However, recent studies demonstrate that more autistic individuals are now receiving an autism diagnosis from mid-childhood onwards than in early childhood4–6. One factor that may explain these findings is a shift in the conceptualisation of autism over time. There is a growing recognition that the signs of autism may not clearly manifest in the first three years of life1,7–9, which has been recognised by the changes to the diagnostic criteria for autism by DSM-5 and ICD-11. Supporting this, several studies have demonstrated that a subset of children who do not initially meet the criteria for an autism diagnosis later receive a diagnosis7,10–14 These findings pose a series of fundamental questions regarding the aetiology of autism. For instance, given the substantial heritability of developmental phenotypes15–17, to what extent does the genetic profile of autism vary with age at receiving an autism diagnosis? How does the developmental variation in the emergence of autism features contribute to age at autism diagnosis, and consequently, the genetic heterogeneity within autism? Is the higher prevalence of mental health diagnosis among autistic individuals diagnosed later in life18,19 partly due to genetic factors? We address these questions using multiple epidemiological and genetic datasets. Using longitudinal data from birth cohorts, we demonstrate that individual differences in socio-behavioural trajectories are associated with age at autism diagnosis. We further demonstrate that age at autism diagnosis is heritable, and this heritability can be partly explained by two correlated polygenic autism factors. The two polygenic factors are differentially associated with developmental and mental health profiles, and partly explain the genetic heterogeneity of autism. We provide a summary of the study and address potential questions regarding the implications of the findings in the **Supplementary FAQs**. ### Socio-behavioural trajectories are linked to varying age at autism diagnosis We first investigated the association between variable developmental trajectories and age at autism diagnosis using four birth cohorts. This included Growing Up in Ireland (GUI, born in 1998), Millennium Cohort Study (MCS, 2000), and Longitudinal Study of Australian Children: Kindergarten cohort (LSAC-K, 1999) and Birth cohort (LSAC-B, 2003) (**Supplementary Table 1, Extended Figure 1, Supplementary Note 1**). All four cohorts collected longitudinal information on socio-behaviour measured using the caregiver-reported Strengths and Difficulties Questionnaire (SDQ)20 and its subscales, and autism diagnosis in data collection sweeps (hereafter “sweeps”) at different ages. The SDQ is widely used, has excellent psychometric properties21–23, and is largely invariant across age, sex, and different populations24–26, meaning that it is measuring the same latent trait across these demographic variables. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/02/2024.07.31.24311279/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2024/08/02/2024.07.31.24311279/F1) Figure 1: Trajectory analyses in three of the four birth cohorts. A - C: Mean SDQ total scores in autistic individuals diagnosed in childhood (blue) and adolescence (orange), and individuals without an autism diagnosis (green) in the MCS (A), LSAC-B (B), and LSAC-K (C) cohorts. Grey regions indicate 95% confidence intervals. D - F: Longitudinal growth mixture models of SDQ total scores among autistic individuals, demonstrating the presence of two groups (green indicating early childhood emergent latent class and purple indicating late childhood emergent latent class) in the MCS (D), LSAC-B (E) and LSAC-K (F) cohorts. G - I: Stacked bar charts providing the proportion of individuals who had been diagnosed as autistic at specific ages, categorised by membership in the latent classes identified from the growth mixture models in MCS (G), LSAC-B (H), and LSAC-K (I) cohorts. Darker colours indicate males and lighter colours indicate females. P-values (inset) are from chi-square tests comparing the distribution of age at autism diagnosis between the two latent classes. Results from GUI have not been plotted here and are available in **Supplementary Figure 6**. Given increasing number of autistic individuals being diagnosed in adolescence4,5, we wondered if there are broad differences in the trajectories of SDQ total and subscale scores among autistic individuals diagnosed before the ages of 9 - 11 (childhood diagnosed group, N = 39 - 118 across cohorts, **Supplementary Table 1**) and after (adolescent diagnosed group, N = 27 -73 across cohorts) (**Methods**). This age cutoff period corresponds to the onset of puberty, the transition from primary to secondary school, and the beginning of an increase of incidence in diagnosis of autism in girls27,28,29. The specific cutoff age was cohort-dependent, as different birth cohorts collected information on autism diagnosis at different ages. We used Latent Growth Curve Models to linearly model the trajectories of SDQ total score and subscales in all four cohorts for both childhood and adolescent diagnosed groups. Across four cohorts, Latent Growth Curve Models identified different trajectories of SDQ total scores between the childhood and adolescent diagnosed groups (Mean scores in **Figure 1A-C, Supplementary Figures 1-6, Supplementary Table 2**). Similar results were obtained for peer relationship problems and prosocial behaviours in all four cohorts (**Supplementary Figures 4 and 5**). Compared with individuals without an autism diagnosis at any time point, the childhood diagnosed autistic group had higher difficulties in early childhood that remained relatively stable or gently declined in adolescence. Compared to the childhood diagnosed autistic group, the adolescent diagnosed autistic group had fewer difficulties during early childhood, but difficulties increased in later childhood and persisted into adolescence. In MCS, we ran a series of sensitivity analyses to check the robustness of the above results. We obtained consistent results: (1) in an expanded sample of MCS which included autistic individuals with co-occurring Attention-Deficit/Hyperactivity Disorder (ADHD) and inconsistent reports of an autism diagnosis (MCS-expanded, **Methods**, **Supplementary Table 2**); (2) after imputing missing SDQ scores (MCS-imputed) (**Supplementary Table 3, Supplementary Note 2**); and (3) when restricting to males. In all four birth cohorts, models stratified by age at diagnosis were generally a better fit to the data than sex-stratifed models (**Supplementary Table 2**), suggesting that the results do not primarily reflect sex differences in age at autism diagnosis (**Supplementary Note 3**). To assess the specificity for autism, we ran Latent Growth Curve Models on SDQ total scores and subscales with children with an ADHD but without an autism diagnosis in the MCS cohort (N = 89). Children with ADHD diagnosed in childhood and adolescence differed only nominally in the slopes of the hyperactivity/inattention (P = 0.026) and prosocial behaviour subscales (P = 0.029) (**Supplementary Table 4 and Supplementary Figures 7 and 8**). In MCS, compared to adolescent diagnosed children with ADHD, adolescent diagnosed autistic children had a steeper increase in peer relationship problems (P = 5.77×10−3) and emotional symptoms (P = 0.012) across development. However, these results must be interpreted cautiously given the low number of children with only ADHD in MCS. Recognising that the age at diagnosis threshold used to categorise autistic individuals into two groups is to some extent arbitrary, we used Growth Mixture Models (GMMs) to identify latent trajectories of SDQ total and subscale scores among autistic individuals in all four cohorts. GMMs do not require *a priori* grouping but identify subgroups based on longitudinal changes in SDQ scores . Across three of the four birth cohorts, GMMs identified two-trajectory models as being optimal for SDQ total scores and a majority of subscale scores (**Supplementary Table 5**). The exception was the GUI cohort, where a one-trajectory model was optimal, likely due to fewer sweeps (three) for SDQ scores. Amongst the other cohorts, one latent trajectory was characterised by difficulties in early childhood, which remained stable or gently declined with age (early childhood emergent latent class). The other latent trajectory was characterised by fewer difficulties in early childhood which increased in late childhood and adolescence (late childhood emergent latent class) (**Figure 1 D-F**). Autistic individuals in the early childhood emergent latent class were more likely to be diagnosed in childhood compared to autistic individuals in the late childhood emergent latent class in MCS (P = 1.43×10−4, chi-square test) and LSAC-B (P = 0.022, chi-square test) (**Figure 1 G-I, Supplementary Table 6**), but this difference was not significant in LSAC-K, possibly because age 11 was the earliest when an autism diagnosis was recorded. Similar results were obtained for all the SDQ subscales in MCS, LSAC-B and LSAC-K except for conduct and peer relationship problems in LSAC-K, where no distinct trajectory groups were identified (**Supplementary Table 5, Supplementary Figures 1-5**). In MCS, these results were largely consistent: (1) in the expanded sample after including individuals with co-occurring ADHD and inconsistent autism diagnoses (**Supplementary Table 4**); (2) after imputation (**Supplementary Table 3, Supplementary Note 4**); and (3) when restricting to only males (**Supplementary Table 5**). In contrast, although we obtained two latent classes based on SDQ and subscale trajectories among individuals with a diagnosis of ADHD but not autism, these latent classes were not significantly associated with age at ADHD diagnosis (**Supplementary Table 4, Supplementary Figures 7 and 8**). In all three cohorts, sex ratio in the late childhood emergent latent class of SDQ total scores compared to the early childhood emergent class was statistically similar (male: female ratio = 3.92 - 2.11; P > 0.05, chi-square tests, **Supplementary Table 6**). After accounting for sex, individuals in the late childhood emergent latent class were more likely than those in the early childhood emergent class to have higher depressive symptoms measured using the Short Mood and Feelings Questionnaire30 (MCS-C, P = 3.84×10−4; LSAC-B, P = 9.82×10−13; LSAC-K, P = 1.86×10−6), have higher rates of diagnosed anxiety (LSAC-B, P = 5.54×10−7; LSAC-K, P = 3.99×10−3) and depression (LSAC-B, P = 3.67×10−4; LSAC-K, P = 4.69×10−3), were more likely to self-harm (LSAC-B, P = 0.018), or have higher rates of suicidal ideation (MCS, P = 2.50×10−3) (**Supplementary Table 7**). Given the significant association between age at autism diagnosis and the GMM latent class membership in LSAC-B and MCS, we wondered if these socio-behavioural trajectories explain any variance in age at autism diagnosis after accounting for socio-demographic factors (e.g., sex, ethnic minority, SES, living area deprivation) and child’s cognitive aptitude, all associated with age at autism diagnosis31,32. Multiple linear regression models indicated that latent class membership from the GMM, socio-demographic factors, and child’s cognitive aptitude together accounted for 17.4% (LSAC-B) - 35.0% (MCS) of the variance in age at autism diagnosis (**Supplementary Table 8**). Latent classes of SDQ total and subscale scores alone explained 9.9% (LSAC-B) - 30.3% (MCS) of total variance. In MCS-expanded, the full model explained 14.8% of the variance, and SDQ total and subscale scores latent class memberships explained 10.0% of the variance. In MCS-imputed, the full model explained 59.8% of the variance and the latent class memberships of SDQ total and subscale scores explained 56.6% of the variance (**Supplementary Table 3**). We further investigated if the effects of socio-demographic factors and cognitive aptitude on age at diagnosis were partly mediated by SDQ latent classes, but did not identify a significant mediated effect. This suggests that the effects of the demographic variables on age at autism diagnosis is largely independent of the effects of SDQ latent classes (**Supplementary Table 8, Supplementary Note 5**). ### Age at autism diagnosis is partly genetic The above analyses demonstrate that variation in socio-behavioural trajectories, measured using the SDQ, is associated with age at autism diagnosis. Previous research has demonstrated that developmental variation in socio-behavioural profiles is partly explained by genetic factors33–38. A corollary of this is that genetic factors may also be associated with age at autism diagnosis. We tested this in a US-based cohort of autistic individuals (SPARK39 : Nmax = 17,105) where we identified a significant heritability estimated from single nucleotide polymorphisms (SNP-based heritability) for age at autism diagnosis (h2SNP = 0.11, s.e.m = 0.01). This heritability did not significantly attenuate after accounting for the child’s developmental phenotypes (age of walking, age at first words, and intellectual disability [ID]), parental socio-economic status that may correlate with greater parental awareness and ability to access diagnostic services, and neighbourhood socioeconomic deprivation which may index availability of healthcare services (**Figure 2A, Supplementary Table 9**). ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/02/2024.07.31.24311279/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2024/08/02/2024.07.31.24311279/F2) Figure 2: Genetic correlates of age at autism diagnosis A. SNP-based heritability of age at autism diagnosis with and without the inclusion of covariates. B. Genetic correlation between age at autism diagnosis genome-wide association studies (GWAS) from SPARK and iPSYCH and psychiatric, neurodevelopmental, and cognitive phenotypes. The Autism GWAS used is the Grove et al., 2019 GWAS. C. Association between PGS for selected neurodevelopmental, cognitive and psychiatric phenotypes and age at autism diagnosis in SPARK. Estimates provided after correcting for ID, sex, age at recruitment into the study, and 10 genetic principal components. Sample sizes provided in inset. D. For those significant in C, association of PGS for schizophrenia, educational attainment, depression, cognitive aptitude and ADHD with age at autism diagnosis in baseline models and after correcting for intellectual disability (ID), developmental (dev.) milestones, socio-economic status (SES) and deprivation. Sample sizes provided in parenthesis. For all plots, points indicate the estimate, whiskers indicate 95% confidence intervals. For plots B-D points with an asterisk (*) indicate significant associations with Benjamini-Yekutieli (BY) adjustment. For plot A, an asterisk (*) indicates significance at P < 0.05 as no multiple testing adjustment is needed for the sensitivity analyses, and whiskers indicate 95% confidence intervals. We used genetic correlation and polygenic score (PGS) analyses to better characterise the genetics of age at autism diagnosis (Nmax = 18,809). Later age at autism diagnosis was significantly positively genetically correlated with ADHD40, and negatively with schizophrenia41 (**Figure 2B and 2C, Supplementary Table 10**). In SPARK, PGS for ADHD and educational attainment remained significantly associated with age at diagnosis after accounting for ID, developmental milestones, socio-economic status and neighbourhood deprivation, and trio status (i.,e, two parents and one child, **Figure 2D, Supplementary Table 11**). As females are significantly more likely to be diagnosed as autistic later than males27, we investigated if there is an interaction effect between sex and the PGS for four phenotypes significantly associated with age at autism diagnosis. We did not identify any significant PGS-by-sex interaction effects (**Supplementary Table 12**). Parents play an important role in recognising autistic features in their children and pursuing an autism diagnosis. Consequently, we wondered if the associations between ADHD and educational attainment (EA) PGS on age at autism diagnosis were due to parental indirect genetic effects (where parental genetics impacts child’s age at autism diagnosis via parental behaviours) or child’s direct genetic effects. In 6,554 trios we observed significant direct effects of ADHD (P = 0.015) and significant indirect effects of EA PGS (P = 6.54×10−4) on increasing age at autism diagnosis respectively (**Supplementary Table 13**). These indirect effects must be interpreted cautiously as they are not immune to confounding, including participation bias. Previous research has demonstrated that autistic individuals are enriched for *de novo* and rare inherited protein truncating or missense variants in genes intolerant to loss of function mutations (constrained genes)42,43. Using trios (N = 6,206), we investigated if rare *de novo* or inherited protein truncating variants or missense variants in highly constrained genes were associated with age at autism diagnosis. We found no significant associations with age at autism diagnosis for either type of variants (**Supplementary Table 11**). This may possibly reflect later autism diagnosis in some carriers of *de novo* mutations due to diagnostic overshadowing by co-occurring ID or global developmental delay. We examined the generalisability of our findings in the Danish iPSYCH cohort (18,965 autistic individuals, mean diagnosis age 10.98 years, SD = 4.64; compared to SPARK’s 4.97 years, SD = 3.28, **Extended Figure 2**). A GWAS of age at autism diagnosis in iPSYCH (N = 19,931) showed similar LDSC-based SNP heritability (h2SNP = 0.10, s.e.m = 0.03) to SPARK (h2SNP = 0.09, s.e.m = 0.03) and moderate genetic correlation between cohorts (rg = 0.51, s.e.m = 0.19, P = 7.56×10−3). The iPSYCH GWAS had positive genetic correlation with depression44 and negative correlation with educational attainment45 and cognitive aptitude46. Differences in genetic correlations may be due to varying age distributions and potential participation bias in SPARK. ![Extended Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/02/2024.07.31.24311279/F3.medium.gif) [Extended Figure 1:](http://medrxiv.org/content/early/2024/08/02/2024.07.31.24311279/F3) Extended Figure 1: Schematic diagram of the cohorts included in the study and the ages when data was collected for SDQ scores (dots) and autism diagnosis (in boxes). MCS = Millennium Cohort Study; GUI = Growing up in Ireland (cohort ’98); LSAC-B = Longitudinal Study of Australian Children (Birth cohort); LSAC-K = Longitudinal Study of Australian Children (Kindergarten cohort). Sample sizes and the year of initial SDQ data collection for each cohort are shown on the ordinate axis. The age cutoff used in the Latent Growth Curve Models for each cohort is indicated by a red line. ![Extended Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/02/2024.07.31.24311279/F4.medium.gif) [Extended Figure 2:](http://medrxiv.org/content/early/2024/08/02/2024.07.31.24311279/F4) Extended Figure 2: Distribution of age at autism diagnosis in SPARK and iPSYCH. Frequency histograms of age at autism diagnosis in iPSYCH and SPARK. Sample sizes have been provided in the inset.**Extended** Figure 3**: SNP heritability for age at diagnosis stratified autism GWAS.** SNP heritability (points) by age at autism diagnosis for varying levels of autism prevalence. Shaded regions, 95% confidence intervals. Each vertical line indicates the best guess autism prevalence for each age at diagnosis stratified autism GWAS. SNP heritability and associated standard error (in parentheses) of autism at the best guess prevalence estimate provided in the top left corner of each facet. ### Characterising the genetic relationship between age at autism diagnosis and autism The previous findings collectively demonstrate that age at autism diagnosis is heritable but with complex genetic correlations that vary by cohort. Subsequently, we sought to characterise the genetic relationship between age at autism diagnosis and autism. Across phenotypes such as schizophrenia47 and depression48, the age at diagnosis/onset is largely negatively genetically correlated with the phenotype itself49. This indicates that earlier diagnosis/onset is associated with greater polygenic propensity for the condition. However, with autism, we observed variable genetic correlations between age at autism diagnosis and different GWAS of autism, including a nominally significant positive genetic correlation with the females-only iPSYCH50 autism GWAS and the SPARK age at autism diagnosis GWAS (**Figure 3A, Supplementary Table 14**). In addition, both age of diagnosis GWAS had moderate negative genetic correlations with both the Psychiatric Genomics Consortium 2017 (PGC-2017) autism GWAS51 and a GWAS of autism in SPARK52, consistent with the average age at autism diagnosis in PGC-2017 and SPARK being lower than that of iPSYCH. The pattern of genetic correlations between age at autism diagnosis and the various autism GWAS does not align well with differences in sex-ratio among the GWAS (**Supplementary Figure 9**), but does align reasonably well with the median age at diagnosis for the autism GWAS. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/02/2024.07.31.24311279/F5.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2024/08/02/2024.07.31.24311279/F5) Figure 3: Genetic correlates of age at diagnosis stratified autism GWAS. A. Genetic correlation between age at autism diagnosis in SPARK and different autism GWAS. Sample sizes are provided in **Supplementary Table 14**. B. Association between age at autism diagnosis PGS for iPSYCHbefore11 and iPSYCHafter10 in the SPARK cohort. Estimates provided after correcting for ID, developmental (dev.) milestones, socio-economic status (SES) and deprivation. C. Over-transmission of PGS for iPSYCHbefore11 and iPSYCHafter10 from parents to autistic children in the SPARK cohort. Estimates provided for unstratified and sex-stratified analyses. Children’s PGS have been standardised to parental mean PGS, with the line at zero indicating no over-transmission. D. Association between autism diagnosed in childhood and adolescence and PGS for iPSYCHbefoere 11and iPSYCHafter10 GWAS in the MCS cohort. For all plots, points indicate the estimate, whiskers indicate 95% confidence intervals. For graphs A, C, and D, points with an asterisk (*) indicate significant associations with Benjamini-Yekutieli adjustment. N indicates sample size. For graph B, points with an asterisk (*) indicate significant association after Bonferroni correction within each sensitivity analysis. N indicates sample size. For graph D, sample sizes are provided as Nautistic/Nnonautistic. Given these varying genetic correlations, we wondered if the polygenic signal for age at autism diagnosis reflects a mixture of different age-dependent polygenic traits. To test this, we conducted GWAS of autism within the iPSYCH dataset, stratifying participants into two groups: those diagnosed before age 11 (iPSYCHbefore11, Nautistic = 9,500) and those diagnosed at age 11 or later (iPSYCHafter10, Nautistic = 9,231). This roughly coincided with the age window of 9 - 11 where we find an increase in SDQ difficulties in the adolescent diagnosed group. We identified moderate positive genetic correlation (rg = 0.70, s.e.m = 0.06) between the two GWAS, which was significantly less than 1 (P = 3.68×10−7). To provide further resolution based on age at diagnosis, we also generated two additional, smaller GWAS of autism stratified by age at diagnosis in iPSYCH: autism diagnosed before age nine (iPSYCHbefore9, Nautistic = 5,451) and autism diagnosed after age 11 (iPSYCHafter11, Nautistic = 8,260). We used the same group of unrelated individuals without an autism diagnosis as controls for all four GWAS (Ncontrol = 36,667). Both genetic correlation and PGS association identified positive shared genetics between iPSYCHafter10 autism and age at autism diagnosis (**Figure 3A and B, Supplementary Tables 14 and 15**), confirming the validity of the age of diagnosis GWAS. Further sensitivity analyses confirmed that the PGS association could not be explained by: (1) Diagnostic overshadowing due to ADHD; (2) Trio status; or (3) Changes in diagnostic criteria between DSM-IV and DSM-5. Furthermore, in SPARK, both sets of PGS were over-transmitted from parents to their autistic children with a larger over-transmission of iPSYCHbefore11 compared to iPSYCHafter10 PGS (P = 1.28×10−3) (**Supplementary Table 15, Figure 3C**). This is consistent with 90% of autistic individuals in SPARK being diagnosed before age 10. We observed consistent results after stratifying by sex and ID (**Supplementary Table 15**). Finally, in MCS, iPSYCHbefore11 but not iPSYCHafter10 was associated with autism diagnosed before age 11 (**Supplementary Table 16, Figure 3D**). Taken together, this suggests that although both sets of PGS are associated with autism, their effects on autism vary by age at diagnosis. Across a range of prevalence estimates, iPSYCHbefore11 had moderately higher SNP-based heritability (h2 = 0.18, s.e.m = 0.02) compared to the iPSYCHafter10 GWAS (h2 = 0.13, s.e.m = 0.01) (**Supplementary Table 17, Extended Figure 3**). The heritability of iPSYCHbefore11 was statistically similar to SPARK (h2 =0.19, s.e.m = 0.03) and PGC-2017 (h2 = 0.20, s.e.m = 0.02), suggesting that stratifying by age at diagnosis identifies similar SNP-based heritability between iPSYCH, SPARK, and PGC-2017. ![Extended Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/02/2024.07.31.24311279/F6.medium.gif) [Extended Figure 3:](http://medrxiv.org/content/early/2024/08/02/2024.07.31.24311279/F6) Extended Figure 3: SNP heritability for age at diagnosis stratified autism GWAS. SNP heritability (points) by age at autism diagnosis for varying levels of autism prevalence. Shaded regions, 95% confidence intervals. Each vertical line indicates the best guess autism prevalence for each age at diagnosis stratified autism GWAS. SNP heritability and associated standard error (in parentheses) of autism at the best guess prevalence estimate provided in the top left corner of each facet. We further investigated whether genetics from age-stratified GWAS supported trajectory modelling findings. In the MCS cohort (N = 6,142 - 5,135), multivariate linear mixed effect models with age-by-PGS interaction showed that iPSYCHafter10 PGS, but not iPSYCHbefore11, was significantly associated with increasing emotional symptoms (BY adjusted P = 7.11×10−4), peer relationship problems (BY adjusted P = 1.75×10−8), SDQ total scores (BY adjusted P = 2.87×10−3), and decreasing prosocial behaviours (FDR adjusted P = 0.030) with age (**Supplementary Table 18, Extended Figure 4**). Consistent results were obtained when including weights to account for sampling bias (**Supplementary Table 18**). ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/02/2024.07.31.24311279/F7.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2024/08/02/2024.07.31.24311279/F7) Figure 4: Genetic correlation heatmaps between different GWAS of autism. Genetic correlation heatmaps of (A) all GWAS of autism, and (B) GWAS of autism after excluding GWAS not stratified by age at diagnosis in SPARK and iPSYCH. GWAS have been ordered based on hierarchical clustering of the genetic correlations. Asterisks (*) indicate significant genetic correlations after Benjamini-Yekutieli adjustment. ![Extended Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/02/2024.07.31.24311279/F8.medium.gif) [Extended Figure 4:](http://medrxiv.org/content/early/2024/08/02/2024.07.31.24311279/F8) Extended Figure 4: Cross-sectional association between PGS for age at diagnosis stratified autism GWAS and socio-behavioural traits measured at different ages. A. Association between PGS for iPSYCHafter10 and iPSYCHbefore11 and scores on the SDQ total and subscales in (A) the MCS cohort at six ages (3 - 17) and (B) ALSPAC at five ages (7 - 17). For all plots, points indicate the estimate, whiskers indicate 95% confidence intervals. In age-by-PGS interaction analyses in ALSPAC (N = 7,172 - 4,977), iPSYCHbefore11 PGS was associated only with decreasing hyperactivity/inattention (BY adjusted P = 4.51×10−3), while iPSYCHafter10 PGS showed nominal associations with increasing peer relationship problems and SDQ total scores with age (**Supplementary Table 18**). In both cohorts, iPSYCHafter10 PGS showed larger increases in effect on SDQ total and peer relationship problems scores from childhood to adolescence compared to iPSYCHbefore11 PGS. Entropy balancing did not alter these findings, suggesting that the differences in results are unlikely to be due to ascertainment differences between cohorts, but may instead reflect secular trends in mental health trajectories53 or variations in the developmental periods analysed. ### Two correlated autism polygenic factors are associated with differing age at autism diagnosis The above age at diagnosis stratified analyses suggested that the genetic signal underlying age at autism is likely a mixture of two or more genetic signals, with varying effects on socio-behavioural trajectories. However, any age-based cutoff for diagnosis is inherently arbitrary. Recognising this, we sought to understand the latent genetic structures in autism using different autism GWAS and their relationship with age at autism diagnosis by modelling the genetic covariances among the different autism GWAS. To enable this and provide greater resolution based on age at diagnosis, we generated three additional age at diagnosis stratified GWAS of autism in SPARK using (unscreened) non-autistic parents and siblings as controls (Ncontrol = 24,965). The three GWAS were: SPARK, diagnosed before age 6 (SPARKbefore6; Nautistic = 14,578); (2) SPARK, diagnosed before age 11 (SPARKbefore11, Nautistic = 18,719); and (3) SPARK, diagnosed after age 10 (SPARKafter10, Nautistic = 3,358). Next, we generated genetic correlations among all the GWAS of autism we had access to. We observed genetic correlations ranging from 0.04 (s.e.m = 0.14) to 0.98 (s.e.m = 0.01) (**Figure 4A, Supplementary Table 19**). This was not entirely explained by cohort differences or sex differences (**Supplementary Note 3**). Hierarchical clustering of the genetic correlations identified two broad, overlapping clusters (**Figure 4A**), one comprising GWAS of autism with predominantly childhood diagnosed individuals, and another comprising GWAS of autism with a large fraction of individuals diagnosed in adolescence or later. This pattern became clearer when excluding GWAS not stratified by age at diagnosis in SPARK and iPSYCH (**Figure 4B**). This is indicative of genetic heterogeneity indexed by age at autism diagnosis. To formally test whether the varying genetic correlation patterns among the different GWAS of autism emerge from different age at diagnosis correlated latent genetic traits, we modelled the genetic covariance using genomicSEM54. GenomicSEM uses structural equation models to identify latent factors. We investigated if two genetic latent traits underlie this heterogeneity, and compared it against five alternative models including a common-factor autism model (**Supplementary Table 20**). Using six minimally overlapping GWAS for autism with wide variation in age at autism diagnosis, we identified a correlated two-factor model that was the most parsimonious and fit the data best (Akaike information criterion [AIC]: 30.09, confirmatory fit index [CFI]: 1, standardised root mean residual [SRMR]: 0.039, **Figure 5A**). ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/02/2024.07.31.24311279/F9.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2024/08/02/2024.07.31.24311279/F9) Figure 5: Two genetic latent factors in autism. A. Path diagram illustrating the two-correlated-genetic-factor models for autism, using six minimally overlapping autism GWAS datasets. F1 = Factor 1, F2 = Factor 2. One-headed arrows depict the regression relationship pointing from the independent variables to the dependent variables. The numbers are the regression coefficients of the factor loadings, with the standard errors provided in parentheses. Covariance between variables are represented as two-headed arrows linking the variables. The numbers on the two-headed arrows can be interpreted as genetic correlation estimates with the standard errors provided in parentheses. Residual variances are represented using a two-headed arrow connecting the residual variable (u) to itself. Standard errors are provided in parentheses. B. Genetic correlation between the two autism factors and a range of mental health, neurodevelopmental, and cognitive traits. Points indicate the estimate, whiskers indicate 95% confidence intervals, and points with an asterisk (*) indicate significant associations with Benjamini-Yekutieli adjustment. Two asterisks (**) indicate phenotypes where the difference in genetic correlation between earlier and later diagnosed autism is statistically significant at P < 0.05. Factor 1 (Earlier diagnosed autism factor) was defined by the GWAS with predominantly early childhood diagnosed individuals (PGC-2017, SPARKbefore6). Factor 2 (Later diagnosed autism factor) was defined primarily by GWAS with adolescent or adult diagnosed individuals (iPSYCHafter10, FinnGen, and SPARKafter10). The cross loading of iPSYCHbefore9 suggests that Factor 2 may impact behaviours in mid/late childhood as well, leading to a diagnosis before age nine. The two factors had a moderate genetic correlation (rg = 0.38, s.e.m = 0.07). Factor 1 was negatively genetically correlated with both the age at autism GWAS, whilst Factor 2 positively genetically correlated with age at autism diagnosis in SPARK (**Figure 5B**), confirming that age at autism diagnosis is linked to genetic heterogeneity in autism. Sensitivity analyses using partly different GWAS identified a two-correlated-factor model as the best fitting model, with similar moderate genetic correlations between the two factors (rg = 0.39 (0.08) - 0.52 (0.09), **Supplementary Table 20, Supplementary Figure 10).** ### Earlier and later diagnosed autism genetic factors are associated with different mental health profiles GenomicSEM analyses revealed at least two correlated autism genetic factors. We wondered if these factors are differentially genetically correlated with cognitive, psychiatric and neurodevelopmental traits. Given the higher prevalence of mental health diagnoses in later-diagnosed autistic individuals18,19, we hypothesised this might partly stem from differences in shared genetics between these autism factors and other mental health and cognitive phenotypes. We tested this hypothesis using genetic correlation analyses. The earlier diagnosed autism factor (Factor 1) was positively correlated with educational attainment and cognitive aptitude but had modest genetic correlations with measures of trauma and ADHD (**Supplementary Table 21, Figure 5B**). The later diagnosed autism factor (Factor 2) had lower genetic correlation with educational attainment but had significant and higher positive genetically correlations with a range of other mental health conditions, including internalising disorders/problems, trauma and related sequelae (Depression, PTSD, childhood maltreatment, and suicide attempts) and ADHD. The iPSYCH female-stratified autism GWAS (iPSYCHfemales) had significantly smaller (P = 0.028) genetic correlation with Factor 1 than the iPSYCH male-stratified autism GWAS (iPSYCHmales), consistent with epidemiological observations that autistic females are, on average, diagnosed later than males27. Sensitivity analyses using age of diagnosis stratified GWAS from iPSYCH and SPARK yielded largely consistent genetic correlation results (**Supplementary Note 6**). We wondered if the genetic signal for later diagnosed autism can be explained by diagnostic misclassification. However, decomposition of the iPSYCH autism genetic signal using genomicSEM indicated that later diagnosed autism cannot be entirely explained by diagnostic misclassification, although genetic effects of ADHD accounted for some of the genetic variance in later diagnosed autism (**Supplementary Note 3**). Accounting for ADHD’s genetic effects revealed attenuated but significant moderate genetic correlations between iPSYCHafter10 and mental health conditions, suggesting shared genetics with ADHD do not fully explain the elevated correlation between later diagnosed autism and mental health phenotypes (**Supplementary Note 6**). We further wondered if there is genetic overlap between the earlier diagnosed autism factor and developmental milestones, as delays in developmental milestones are the earliest indicator that a child may be autistic55. The earlier but not later diagnosed autism factor was positively and significantly genetically correlated with greater difficulties in social behaviour at age three (**Supplementary Table 22, Figure 6A**). Supporting these findings, PGS for iPSYCHbefore11 but not iPSYCHafter10 GWAS was associated with greater difficulties in social communication (gestures) at 15 months (**Figure 6B, Supplementary Table 23**). Nevertheless, the genetic correlations and regression coefficients did not statistically differ between the two GWAS. Furthermore, PGS for neither iPSYCHbefore11 nor iPSYCHafter10 were associated with later attainment of developmental milestones among autistic individuals from SPARK. However, these results may reflect collider bias as autistic carriers of rare genetic variants within the cohort have lower autism PGS and substantially delayed developmental milestones50 (**Supplementary Table 24, Figure 6C**). ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/02/2024.07.31.24311279/F10.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2024/08/02/2024.07.31.24311279/F10) Figure 6: Association between age at diagnosis stratified autism and developmental milestones. A. Genetic correlation between earlier and later diagnosed autism genetic factors and a range of developmental phenotypes. mos = months. B. Association between PGS for iPSYCHbefore11 and iPSYCHafter10 and social communication skills at 15 months in ALSPAC. N is the sample size. C. Association between PGS for iPSYCHbefore11 and iPSYCHafter10 and developmental milestones among autistic individuals in SPARK. N is the sample size. For all plots, points indicate the estimate, whiskers indicate 95% confidence intervals, and points with an asterisk (*) indicates significant associations after Benjamini-Yekutieli adjustment. For all plots, positive values indicate greater difficulties/delays. ## Discussion Understanding why some autistic individuals receive a diagnosis earlier than others has been an important scientific and societal question. Here we show, using multiple methods and datasets, that some of the variability in age at autism diagnosis is linked to differences in socio-behavioural trajectories and associated genetic profiles among autistic individuals. The genetic variation associated with age at autism diagnosis is also associated with the genetic heterogeneity within autism. This relationship can partly explain the often contradictory patterns of genetic correlations between autism and various cognitive, neurodevelopmental, and psychiatric phenotypes across different autism GWAS. This axis of genetic heterogeneity within autism, indexed by age at autism diagnosis, is not fully explained by several other factors that may influence age at autism diagnosis including sex, co-occurring ID and developmental delays, changes to the diagnostic criteria, cohort differences, diagnostic misclassification, or parental factors influencing diagnosis (**Supplementary Note 3**). Our findings are consistent with the wider literature that demonstrates that genetic influences on social-communication vary across development in the general population37,38. However, extending this line of enquiry, we show that genetic influences on clinically diagnosed autism too vary based on age at diagnosis. Modelling the genetic correlations among various GWAS of autism identified two correlated autism polygenic factors that explained the data better than the alternative models considered, including a single autism factor model (**Figure 5**). Notably, the genetic correlation between the factors is 0.38, which is similar to the genetic correlation between depression and schizophrenia56. It is likely that other dimensions contribute to heterogeneity in autism, including potentially further genetic differences based on age at diagnosis. For example, a significant proportion of the variation in the FinnGen autism GWAS was not explained by either of the two factors (**Figure 5A**) and the correlation between some autism GWAS (e.g., iPSYCHfemales and SPARKbefore6) is even lower than the genetic correlation observed between earlier and later diagnosed autism factors. Importantly, the low genetic correlation between the factors strongly suggests that later diagnosed autism is not merely a broader autism subtype or that they emerge from the same underlying latent genetic distribution. This is further supported by differences in the patterns of genetic correlation between the two autism genetic factors and other phenotypes, and trajectory modelling of SDQ total and subscale scores. Of the two autism polygenic factors, we found some evidence to suggest that the earlier diagnosed autism genetic factor was associated with social communication difficulties in early childhood (**Figure 6**). However, neither of the two autism polygenic factors are predominantly driven by a subset of children with co-occurring ID or neurodevelopmental delays. The heritability of age at autism diagnosis does not attenuate when accounting for ID or developmental delays, and there is no association between either factor and neurodevelopmental conditions or delays in developmental milestones (for example, significant delays in walking or first words). Our findings imply that profound neurodevelopmental disabilities among autistic individuals may be aetiologically distinct from the two autism polygenic factors. Elevated polygenic propensity for the later diagnosed autism genetic factor may lead to less clear difficulties in early childhood and thus may not be picked up by caregivers as reasons to pursue diagnosis or support. These are consistent with findings that children who do not initially meet the criteria for an autism diagnosis may later meet them7,10, and with parental reports of on average “milder” autism features among later diagnosed autistic individuals57. In contrast, both the later diagnosed genetic autism factor (**Figure 5**) and the late emergent latent class of SDQ total scores (**Figure 1**) are associated with greater mental health problems, particularly internalising difficulties, self-harm, and correlates of trauma. The later diagnosed autism genetic factor is also associated with a larger increase in socio-behavioural difficulties in adolescence. Yet again, this is consistent with the epidemiological findings of greater mental health difficulties among later diagnosed autistic individuals18,19. However, our findings demonstrate that part of the epidemiological findings may be explained by genetic heterogeneity in autism indexed by age at autism diagnosis. How exactly these genetic differences lead to greater mental health problems remains to be resolved. The variation in genetic correlation between ADHD and autism stratified by age of diagnosis is particularly noteworthy. Older GWAS of autism (including the PGC-2017) were not significantly genetically correlated with ADHD51,58,59 whereas more recent GWAS for autism have moderate genetic correlations with ADHD60. Genetic correlation analyses (**Figure 5**) indicate that the genetic correlation with ADHD increases with a later diagnosis of autism. We confirmed this using within-family analyses: autistic individuals diagnosed before age five do not over-inherit PGS for ADHD (**Supplementary Table 13c**). However, even within ADHD, there is genetic heterogeneity based on age at diagnosis61. ADHD diagnosed in childhood shows a larger genetic correlation with autism and lower genetic correlation with internalising disorders62 compared to late-diagnosed ADHD. These results suggest a complex genetic interplay between autism and ADHD that is dependent on age at diagnosis. Several sensitivity analyses indicate that our findings are not primarily capturing sex differences (**Supplementary Note 3**). However, given that autistic females are, on average, diagnosed later than males4,5, research that investigates sex and gender differences in both autism and co-occurring conditions63 needs to account for genetic confounding associated with age at autism diagnosis. Findings that may be thought to reflect sex differences may additionally reflect differences in age at diagnosis. For example, the higher prevalence of mental health problems in autistic females compared to males64 attenuates when restricting to autistic individuals diagnosed before age five18. Although we explored the impact of several additional clinical and demographic factors on age at autism diagnosis, these account for less than 50% of the variance in age at diagnosis (**Supplementary Table 8** and **Supplementary Note 5**). Notably, there is substantial variation across the datasets explored, highlighting that age at diagnosis of autism is immensely complex, varying across geography and time. Local cultural factors, access to healthcare, gender bias, stigma, and camouflaging, all of which are difficult to measure, likely have an impact on who receives a diagnosis and when. Of interest is camouflaging, which has been hypothesised as one reason for later diagnosis65 particularly among autistic females. Our results cannot be fully explained by camouflaging. Although children as young as seven years of age may camouflage66, it is unlikely that infants can camouflage behaviours or developmental milestones in infancy or early toddlerhood. Furthermore, even among autistic individuals, there is variation in camouflaging67,68. Although we are unable to explicitly test the impact of camouflaging, our results are consistent with the correlates of camouflaging. For instance, it is known that autistic females are more likely to camouflage69, and higher camouflaging is associated with greater mental health difficulties70, and later autism diagnosis65,71. In conclusion, using genetic data and longitudinal analyses of birth cohorts, we identify an axis of heterogeneity in autism which is indexed by age at autism diagnosis. This axis of heterogeneity partly explains the varying genetic correlations among the different GWAS of autism and between autism and various mental health conditions. These findings provide further support to the hypothesis that the umbrella term “autism” describes multiple phenomena with differing aetiologies, developmental trajectories, and correlations with mental health conditions. These findings have implications for how we conceptualise neurodevelopment more broadly, and for understanding diagnosis, sex and gender differences, and co-occurring health profiles in autism. ## Methods ### A note on terminology We use the term autistic to refer to people who have an autism diagnosis72. We use non-autistic to refer to people who do not have an autism diagnosis. We use sex to refer to sex assigned at birth, and use the terms males and females to refer to sex. ### Analyses of birth cohorts Study design and participants We used four population-based birth cohorts that vary in ages covered and when data were collected (**Extended Figure 1**). Briefly, the four cohorts included are the UK-based Millennium Cohort Study (MCS)73, the Ireland-based Growing Up in Ireland (GUI) Child cohort (aka Cohort 98’)74, and the Australia-based Longitudinal Study of Australian Children - Birth (LSAC-B) and Kindergarten (LSAC-K) cohorts75. All children included in the cohorts were born in the 21st century. Further details about the cohorts are provided in the **Supplementary Note 1.** As indicated in **Supplementary Table 1**, these cohorts were adopted for their longitudinal nature, the nationally representative cohort members, and the availability of data on behavioural profiles and neurodevelopmental diagnosis allowing for cross-country comparisons and generalisation76. ### Measures #### Autism and ADHD diagnosis and age at diagnosis In all cohorts, across multiple sweeps, the main caregiver was asked if the participant had a diagnosis of autism (**Extended Figure 1**). For age at diagnosis, we used the age at the sweep when participants first reported being diagnosed as autistic in every cohort, to maximise sample sizes and ensure consistency across cohorts for effective comparisons. For instance, if a participant first reported an autism diagnosis at the age 11 sweep, we considered age at diagnosis to be 11 years. Although the specific age at diagnosis was provided for LSAC-B and LSAC-K, we opted not to use this, as we identified errors in some reports where months and years of diagnosis were swapped or not reported. In MCS, we had reports of both autism and ADHD diagnoses, allowing us to conduct several sensitivity analyses. For our primary analyses, we included a narrowly defined sample of children with consistently reported autism diagnoses by primary and proxy caregivers (when both were available) and no other reported neurodevelopmental diagnosis (particularly ADHD). To assess the generalizability of our results and increase the sample size, we then expanded the sample to include all children with any reported diagnosis of autism. This expanded sample, which we refer to as “MCS-expanded”, included cases regardless of whether the diagnoses were consistent across sweeps or caregivers, and included those with co-occurring ADHD. Additionally, we imputed the independent variables and covariates for autistic individuals with missing information, as detailed below. We refer to this sample as “MCS-imputed”. Finally, to assess the specificity of the trajectories for autism, we conducted analyses among children who had a consistent ADHD diagnosis but no diagnosis of autism. We refer to this sample as “MCS-ADHD”. #### Strengths and Difficulties Questionnaire (SDQ) We used the SDQ to capture social, emotional, and behavioural profiles of participants, with repeated measures from 3 to 18 years across cohorts (**Extended Figure 1**). SDQ comprises 25 statements that respondents are asked to rate on a 3-point Likert scale (“not true”, “somewhat true”, and “certainly true”) based on the child’s symptoms or behaviours over the past six months. There are five subscales, each containing five items, which assess emotional symptoms, conduct problems, hyperactivity-inattention, peer relationship problems, and prosocial behaviours respectively20. The first four subscales assess difficulties, and the total score ranges from 0 to 40, with higher scores indicating more significant difficulties. The fifth subscale represents strengths and has a total score ranging from 0 to 10, with higher scores indicating more prosocial behaviours. The SDQ demonstrates good test-retest reliability and criterion validity across countries21–23. Each subscale of SDQ has been found to exhibit correlations with diagnosis of autism and ADHD77. Its five-factor structure (each subscale as a factor) has shown consistency and invariance across age, sex, and ethnic background21,25. The SDQ captures several core features of mental health and neurodevelopmental conditions, including autism and ADHD78. Only children with complete data of SDQ across all sweeps were included in the analyses, except for imputation analyses. #### Socio-demographic measures Socio-demographic measures were included as covariates to account for their impact on age at diagnosis in each cohort (**Supplementary Table 25**). Specific measures and available information vary across cohorts, but we generally included sex, ethnic background, maternal age at delivery, child’s cognitive aptitude, household socio-economic status (SES), and deprivation level of the living area, to account for some factors that may impact the age when someone receives an autism diagnosis32,79. Only subsets of children in the complete-SDQ samples, with complete data for these socio-demographic factors, were included in the respective analyses, resulting in a further reduction in sample sizes. In MCS, although various census classifications for ethnic groups were available, we opted to use a binary indicator to identify non-white ethnic minorities. This approach was chosen to maintain consistency with other cohorts. Ethnicity data were not collected in either LSAC cohort. Instead, visible ethnic minority status was determined primarily by parental country of birth and the language(s) spoken at home.80 Maternal age at delivery was collected only in MCS. In other cohorts, we used maternal age (in years) at first sweep of data collection to reflect the variation in maternal age at delivery. For cognitive aptitudes and other socio-demographic factors, including SES and deprivation, we adopted summary scores using principal component analysis (PCA) to capture information measured by diverse scales. However, as including more social factors lead to smaller sample sizes, we prioritised factors based on the availability of data among already limited autistic samples. Information on socio-demographic factors in each cohort, including variables and scales included in PCA and resulting sample sizes, can be found in **Supplementary Table 25**. Note, in the birth cohorts, no autistic child was identified as having intellectual disability (ID), defined as scoring two or more standard deviations below the mean value of the first principal component score (‘g’ factor) derived from multiple cognitive aptitude measures in corresponding cohort. ### Statistical analyses Following the participant selection process, we used two methods to model the longitudinal trajectories of SDQ total and subscale scores. In the first analyses, we *a priori* defined two groups of autistic individuals - childhood diagnosed (diagnosed before ages 9 - 11, depending on the cohort), and adolescent diagnosed (diagnosed after the ages of 9 - 11, depending on the cohort). This age period was chosen as there is epidemiological evidence showing increased autism incidence among females during this window27 and because trajectory analyses have identified increasing autism-related traits in a subset of individuals after this period 81. We were also limited in choosing alternate cutoffs due to the absence of information on both SDQ and autism diagnosis at earlier and later time points in some of the cohorts, and the relatively low sample sizes of the resulting groups. Anyone who had no report of an autism diagnosis were included in the general population group. For MCS in particular, children with neither autism nor ADHD diagnosis were included in the general population group. We used linear Latent Growth Curve Models (LGCM) to identify the latent trajectories of SDQ total and subscale scores in the three groups (childhood diagnosed, adolescent diagnosed, and the general population) for all cohorts. Each linear model included a latent intercept to represent the initial level of the outcome variable, and a linear latent slope to represent the mean rate of change over time. As sensitivity analyses, quadratic models were also fitted in MCS, MCS-expanded, and LSAC-B, in which quadratic time scores were assigned across sweeps respectively to capture this nonlinear change over time. However, quadratic models for most subscales among the three cohorts either did not converge or demonstrated Heywood cases, i.e., negative variance estimates, in slope terms This likely stems from insufficient statistical power due to the small sample sizes or model misspecification, which hinders meaningful theoretical interpretation. Therefore, we decided to use linear models. Given the well-known sex differences in diagnosis age27, we also applied the same models stratified by sex, i.e., estimating latent intercept and slope for each sex, within the autistic samples. All LGCM were fitted under the structural equation modelling framework using the *lavaan* package in R82. In parallel, we conducted Growth Mixture Models (GMM) to identify if there were latent groups of autistic individuals based on their trajectories of SDQ total and subscale scores. GMM assumes that the sample being studied consists of multiple mixed effects models, each capturing a subgroup trajectory with shared intercept and slope83. We fitted models with one to four groups for each subscale and SDQ total scores in each cohort, using the *lcmm* package in R84. The optimal number of latent trajectories were then determined by comparing fit indices, including Bayesian Information Criterion (BIC) values, classification quality measure (entropy), and substantive interpretation. Models with lower BIC values and higher entropy are favoured85. Also, models identifying subgroups with less than 5% of the sample size were not considered for poor statistical reliability and limited practical significance86. Multiple regressions were subsequently conducted to investigate the association between individuals’ age at diagnosis (the outcome variable) and their SDQ total and subscale latent class memberships identified in optimal GMM, as well as other socio-demographic covariates. We did not detect any multicollinearity among the variables using variance inflation factors. Estimates of coefficients of predictors and corresponding p-values were interpreted to determine which factors contribute to differing age at autism diagnosis. Additionally, considering the limited sample sizes and the number of explanatory variables included, the relative importance of each predictor was assessed using dominance analysis87. We employed the *misty*88 package in R for this analysis, using a correlation matrix extracted from the fitted model via the *lavInspect* function from the *lavaan* package82. This approach leverages the correlation matrix to consider not only individual predictors but also the correlations among them, providing a more comprehensive assessment of their relative importance89. To examine potential causal pathways, mediation analyses were conducted, allowing socio-demographic factors to indirectly influence the age at diagnosis through their effects on latent class memberships identified in the optimal GMM. Using structural equation modelling in the *lavaan* package82, both direct and indirect effects were assessed, with their significance calculated using bootstrapping analysis. Further details are provided in **Supplementary Note 5**. To investigate the specificity of our findings to autism, we conducted LGCM, GMM, and sequential mediation analyses in individuals with ADHD but without a co-occurring autism diagnosis in the MCS cohort (N = 89, **Supplementary Table 4**). ADHD diagnoses were available in the same sweeps as autism diagnoses, reported at age 5,7,11, and 14. Carers were asked the following question: ‘Has a doctor or other health professional ever told you that had Attention Deficit Hyperactivity Disorder (ADHD)?’. In autistic individuals, using the GMM based latent classes of the SDQ total scores, we used multiple regression to investigate the association with mental health phenotypes in MCS, LSAC-K, and LSAC-B. We included sex as a covariate. #### Imputation To assess the impact of missingness, we applied *softImpute*, to impute missing data for all children with an autism diagnosis reported by any carer in any sweep in the MCS cohort (N = 623, **Supplementary Table 3**). Given the longitudinal nature of data collections for SDQ subscale scores and some cognitive aptitude measures, SoftImpute was chosen for its computational efficiency in handling large-scale matrices through low-rank approximation, effectively preserving underlying structure of input data. To enhance imputation quality and reduce bias, we included related auxiliary variables in the imputation process, along with SDQ subscale scores in all available sweeps (**Supplementary Table 3**). Further information is provided in **Supplementary Note 2**. ### SPARK cohort: Genotyping, quality control and imputation We used data from the SPARK cohort39 iWES2 v1 dataset (released in Feb 2022) which included data from 70,487 autistic individuals and their families. All participants were genotyped using the Illumina Global Screening Array (GSA_24v2-0_A2). To avoid false positives due to fine-scale population stratification, we restricted the analyses to individuals of genetically-inferred European ancestries (N = 51,869 autistic and non-autistic participants), which was provided by the SPARK consortium. From this, we restricted to individuals with genotyping rate > 98%, individuals without sex mismatches and excess heterozygosity (3 standard deviations from the mean heterozygosity), and where trio data was available, trios with fewer than 5% Mendelian errors, resulting in 47,170 autistic and non-autistic individuals. We included genetic variants with minor allele frequency > 1%, genotyping rate > 95%, and that were in Hardy Weinberg Equilibrium (HWE p-value > 1×10−6), resulting in 518,189 SNPs. We used this quality controlled genotype data for imputation, calculating genetic principal components, and inferring relatedness among individuals. We inferred genetic relatedness using KING90. For genetic principal component analysis, we pruned SNPs for linkage disequilibrium (LD) (maximum r2 = 0.1) and removed the human leukocyte antigen (HLA) region. Using PC-AiR91, we first calculated principal components (PCs) in genetically unrelated individuals and then projected the PCs onto related individuals. We imputed genotypes using the TOPMED imputation panel92 on the Michigan imputation server (v1.7.3)93 using Minimac493 and after phasing using Eagle v2.594. Post imputation, variants were converted from GRCh38/hg38 to GRCh37/hg19 using liftOver. We restricted downstream analyses only to variants with minor allele frequency > 0.1% and with an imputation R2 > 0.6. ### SPARK cohort: Association analyses #### PGS association analyses Polygenic scores (PGS) were calculated using PRScs95 which uses a Bayesian shrinkage prior. PGS were calculated for autism (iPSYCH only dataset, N = 19,870 autistic individuals and 39,078 non-autistic individuals)50,96, ADHD40, bipolar disorder97, major depressive disorder44, schizophrenia41, educational attainment98, and cognitive aptitude46, autism diagnosed before age 11 (iPSYCHbefore11), and autism diagnosed after age 10 (iPSYCHafter10). The latter two GWAS were generated using the iPSYCH201596 cohort, details of which are provided below. For simplicity we refer to this cohort as iPSYCH throughout. We ran association analyses between PGS and age at autism diagnosis (converted to years in all analyses) in the quality controlled dataset. We excluded individuals older than 22 to focus on those who had an autism diagnosis using either the DSM-IV99 or DSM-5, retaining a maximum of 18,809 autistic individuals for PGS analyses. This criteria also allowed us to focus on individuals who received their diagnosis in childhood or adolescence, as older adults may have missed an earlier diagnosis of autism due to secular changes in social attitudes towards autism. For psychiatric conditions, we ran separate linear regressions with age at autism diagnosis and the aforementioned PGS. The baseline model included ID (caregiver reported), sex, and the first 10 genetic principal components as covariates. For schizophrenia, ADHD, depression, cognitive aptitude, educational attainment, iPSYCHbefore11 and iPSYCHafter10, we ran sensitivity analyses by sequentially including age at walking and age at first words (developmental milestones), parental occupation, highest parental education, and household income (together, socio-economic status or SES), and national area deprivation percentile (deprivation) as covariates. We also included trio status in the baseline model to account for potential participation bias. Additionally, for iPSYCHbefore11 and iPSYCHafter10 we included any diagnosis of an attentional or behavioural disorder as a covariate in the baseline model to account for diagnostic overshadowing. For the PGS with iPSYCHbefore11 and iPSYCHafter10, we ran sensitivity analyses after stratifying by sex. We tested if the effects varied by sex by including a PGS by sex interaction term in the baseline model. We tested for direct and indirect genetic effects of ADHD and educational attainment (EA) PGS in two ways. First, we generated pseudocontrols for complete trios in Plink 1.9100,101, and calculated PGS separately for autistic individuals and pseudocontrols. We regressed the effects of the ADHD and EA PGS for autistic individuals and pseudocontrols on age at autism diagnosis. We included sex, ID (caregiver reported), and the first 10 genetic principal components as covariates. The direct effects were calculated by subtracting the effects of the untransmitted PGS (indirect effect) from the transmitted PGS (total effects). Standard errors and p-values were calculated by bootstrapping 10,000 times as done previously102. Second, in complete trios, we used polygenic transmission disequilibrium tests (pTDT)103 to calculate the deviation of PGS from the parental mean PGS for ADHD and EA. We checked for over-transmission of PGS stratified by sextiles based on age at autism diagnosis. We used pTDT to investigate if there is an over-transmission of PGS for iPSYCHbefore11 and iPSYCHafter10 among autistic individuals in the SPARK cohort. In the SPARK cohort, we obtained data for age at achieving nine developmental milestones (in months) for autistic individuals. For all milestones, we excluded individuals who were greater than five median absolute deviations from the median. We ran multiple linear regression with PGS for iPSYCHbefore11 and iPSYCHafter10 in which we included sex, age at recruitment into the study, and the first 10 genetic principal components as covariates. #### Rare high-impact de novo variants and inherited variantes We identified rare (minor allele frequency [MAF] < 0.1%) *de novo* and inherited variants in complete trios from SPARK as previously described104. We identified high impact protein truncating variants by restricting variants in loss-of-function observed/expected upper bound fraction (LOEUF)105 highly constrained decile (LOEUF < 0.37) that were annotated as either frameshift, stop gained, or start lost; and had a loss-of-function transcript effect estimator (LOFTEE) “high confidence annotation”. To identify high-impact *de novo* missense variants, we restricted to variants in LOEUF highly constrained genes (LOEUF < 0.37), and had an MPC (missense badness, PolyPhen-2, and constraint) score106 > 2. All variants were rare, with an allele frequency < 0.1% in SPARK and gnomAD. We ran regression analyses separately for high-impact *de novo* and inherited protein truncating variants and missense variants, and additionally by combining both protein truncating and missense variants. We included sex and age at recruitment into the study as covariates for analyses with *de novo* variants, and additionally the first 10 genetic principal components for analyses with inherited variants. ### GWAS for age at autism diagnosis and age at diagnosis stratified autism We generated a GWAS of age at autism diagnosis (in years) on the quality controlled dataset from SPARK, restricting it to individuals who were under 22 years of age (N = 18,809), and SNPs with a MAF > 1%. GWAS was generated using FastGWA107 with age at recruitment into the study, sex, ID, and the first 10 genetic principal components as covariates. In iPSYCH, we generated an additional GWAS of age at autism diagnosis (in years) in a quality controlled dataset of unrelated individuals with sex and ID included as covariates using FastGWA107, restricting to SNPs with an MAF > 1%. To keep it consistent with SPARK, we excluded individuals who were diagnosed after age 22, resulting in 18,965 individuals. Briefly, pre-imputation quality control of the iPSYCH data was performed using the Ricopili pipeline108, prephased using Eagle v.2.3.5, and imputated using Minimac3109, using the downloadable version of the Haplotype Reference Consortium (HRC)110 (accession no. EGAD00001002729). Further details of quality control and imputation are provided in Als et al., 202344. We additionally generated three age at autism diagnosis stratified GWAS in SPARK using (unscreened) non-autistic parents and siblings as controls (Ncontrol = 24,965). The three GWAS were: (1) SPARK, diagnosed before age 6 (SPARKbefore6; Nautistic = 14,578); (2) SPARK, diagnosed before age 11 (SPARKbefore11, Nautistic = 18,719); and (3) SPARK, diagnosed after age 10 (SPARKafter10, Nautistic = 3,358). For these analyses, we did not restrict it to individuals under 22 to increase sample size. GWAS were generated using quality controlled SNPs with a MAF > 1% using FastGWA-GLMM111. We included age at recruitment into the study (to account for parental controls potentially lacking autism diagnoses due to historical diagnostic changes), sex and the first 10 genetic principal components as covariates. Fast-GWA GLMM can account for relatedness and fine-scale population stratification even in family-based samples like SPARK. Although inclusion of unscreened related individuals as controls can decrease heritability and statistical power to identify loci112, we used the GWAS to primarily conduct genetic correlation and related analyses. To ensure the robustness of these models we: (1) confirmed that the attenuation ratio for all GWAS was not significantly greater than 1; (2) generated an additional GWAS of SPARK without stratifying by age at autism diagnosis using the same methods and confirmed a high genetic correlation (rg = 0.92, s.e.m = 0.17) with a previous SPARK GWAS52 which used a case-pseudocontrol approach; and (3) in the genomicSEM analyses, ran sensitivity analyses using a trio-based SPARK GWAS in lieu of the age at diagnosis stratified GWAS from SPARK and confirmed our findings. We generated four age at diagnosis stratified GWAS of autism in iPSYCH cohort96. The primary GWAS used in the analyses were GWAS of autism diagnosed before age 11 (iPSYCHbefore11: 9,500 autistic and 36,667 non-autistic individuals) individuals and autism diagnosed after age 10 (iPSYCHafter10: 9,231 autistic and 36,667 non-autistic individuals). We chose this age cutoff to divide the iPSYCH cohort into two subgroups with similar sample sizes and because age coincided with the window in which we observe an increase in SDQ scores in the birth cohorts, and which is associated with an increase in diagnosis of females in epidemiological samples27. We conducted two additional GWAS with smaller sample sizes: GWAS of autism diagnosed before age nine (iPSYCHbefore9: 5,451 autistic and 36,667 non-autistic individuals) and after age 11 (iPSYCHafter11: 8,260 autistic and 36,667 non-autistic individuals). These were used in sensitivity analyses. For the last two GWAS, we also conducted GWAS after excluding individuals born after 1994, to ensure that all autistic individuals received a diagnosis using either DSM-IV or DSM-5 criteria. However, we observed high genetic correlation between the GWAS when using the full sample and when restricting the sample to those born after 1994, suggesting that changes in the diagnostic criteria do not substantially impact the genetic analyses. To increase sample size and statistical power, we conducted all downstream analyses without excluding autistic individuals born before 1994. All individuals included in these GWAS from iPSYCH were born between May 1980 and December 2008 to mothers who were living in Denmark. GWAS was conducted on individuals of European ancestry, with the first 10 genetic principal components included as covariates using logistic regression as provided in PLINK. ### Heritability, genetic correlation, and genomicSEM Heritability analyses for age at autism diagnosis were conducted using a single-component genome-wide complex trait analysis with genomic-relatedness-based restricted maximum likelihood approach (GCTA-GREML)113,114 in unrelated autistic individuals using the quality controlled genetic data in SPARK. We estimated SNP-based heritability first after including sex, age, and the first ten genetic principal components as covariates. We ran sensitivity analyses after sequentially including ID, developmental milestones, SES and area deprivation as covariates. SNP-based heritability for the age at diagnosis stratified autism GWAS from iPSYCH was calculated using linkage disequilibrium score regression coefficient (LDSC)58,115 using linkage disequilibrium scores from the north-west European population. We converted observed scale heritability estimates to liability scale estimates using a range of autism lifetime prevalence estimates, including a “best guess” autism lifetime prevalence estimate for each of the age-stratified autism GWAS. We conducted genetic correlation analyses using LDSC, using linkage disequilibrium scores from the north-west European populations. For genomicSEM54 analyses, we first conducted genetic correlation analyses among fourteen different autism GWAS using LDSC. This included a multi-ancestry case-pseudocontrol GWAS in SPARK52 (6,222 case-pseudocontrol pairs); GWAS from FinnGen (Data Release - r10)116 (646 cases and 301,879 controls), the PGC-2017 autism GWAS51 (7,387 cases and 8,567 controls), seven GWAS from iPSYCH, and age at diagnosis stratified GWAS from SPARK. The iPSYCH GWAS included an unstratified (19,870 autistic individuals [15,025 males and 4,845 females] and 39,078 controls) and sex-stratified GWAS50, and four age at diagnosis stratified GWAS as mentioned earlier. Subsequently, we restricted to six GWAS with minimal sample overlap, without high genetic correlation (rg > 0.95), and with wide variation in age at diagnosis to conduct genomicSEM analyses using autosomes. Using the patterns of genetic correlations observed we tested an age at diagnosis related correlated two-factor model. We additionally tested: (1) a single-factor model; (2) a correlated two-factor “geography” model where three US-based autism GWAS loaded onto one factor, and three Europe-based autism GWAS loaded onto a second factor; (3) a bifactor model based on age at diagnosis; (4) a bifactor model based on the geography of the cohorts; and (5) a hierarchical factor model based on age at diagnosis. The two-factor model was chosen as it had lower RMSEA and higher CFI and was more parsimonious than the bifactor model. We ran sensitivity analyses using different GWAS of autism as input and confirmed that the two-correlated-factor model was the best fitting model of the models tested. ### Analyses in ALSPAC and MCS #### Genetic quality control We obtained quality controlled and imputed genotype data from ALSPAC117–119. Further details about the cohort are provided in the **Supplementary Note 1.** Briefly, ALSPAC children were genotyped using the Illumina HumanHap550 quad chip genotyping platforms by 23andme. Individuals were excluded due to sex mismatches, excess heterozygosity, missingness > 3%, and insufficient sample replication (Identical-By-Descent [IBD] < 0.8). After multidimensional scaling, and comparison with Hapmap II (release 22), only individuals of genetically inferred European ancestries were retained. SNPs with low frequency (MAF < 1%), poor genotyping (call rate < 95%) and deviations from Hardy-Weinberg equilibrium (P < 5×10−7) were removed. 9,115 subjects and 500,527 SNPs passed quality control. Genotypes were phased using ShapeIT, and imputation was done using the Haplotype Reference Consortium panel using the Michigan imputation server. After imputation, we further removed low frequency SNPs (MAF < 1%). Further details of the quality control and imputation of ALSPAC are provided here: [https://proposals.epi.bristol.ac.uk/alspac\_omics\_data\_catalogue.html#org89bb79b](https://proposals.epi.bristol.ac.uk/alspac_omics_data_catalogue.html#org89bb79b). Genome-wide genotype data was generated by Sample Logistics and Genotyping Facilities at Wellcome Sanger Institute and LabCorp (Laboratory Corporation of America) using support from 23andMe. We also obtained quality controlled and imputed data from MCS. Briefly, MCS samples were genotyped using the Illumina Global Screening Array120. Individuals were excluded due to sex mismatches, excess heterozygosity, and missingness > 2%. We identified European samples using the GenoPred pipeline121 ([https://github.com/opain/GenoPred](https://github.com/opain/GenoPred)). SNPs with low frequency (MAF < 1%), poor genotyping (call rate < 97%) and deviations from Hardy-Weinberg equilibrium (p < 1×10−6) were removed. Imputation was conducted using Minimac493 using the TOPMED reference panel92 in the Michigan imputation server93. Post imputation, SNPs with an imputation R2 INFO score < 0.8, with > 3% missing, and with a MAF < 1% were excluded. Further details are available here: [https://cls-genetics.github.io/docs/MCS.html](https://cls-genetics.github.io/docs/MCS.html) PGS for both ALSPAC and MCS were calculated in individuals of genetically inferred European ancestries. Genetic principal components were calculated for both cohorts using PC-AiR as described earlier. We calculated PGS for iPSYCHbefore11 and iPSYCHafter10 and used these in all analyses in the MCS to keep it consistent with analyses in SPARK where we could only use the iPSYCH GWAS to avoid overlap between the training and testing sample. #### Association with SDQ We obtained scores on the SDQ total and subscales for six ages in the MCS and five ages in ALSPAC. We ran cross-sectional analysis at each age using multiple linear regression with PGS for iPSYCHbefore11 and iPSYCHafter10, with sex, age, and the first 10 genetic principal components as covariates. Additionally, we ran multiple linear mixed effects regression using *lme4* package in R122, fitting a PGS by age interaction term to investigate if the effects of PGS on SDQ change over time. To investigate if the differences in association between MCS and ALSPAC were due to differences in ascertainment between the two cohorts, we matched ALSPAC to MCS using entropy balancing123 and re-ran the PGS association analyses. Entropy balancing is a reweighting technique that ensures the covariate distributions are identical between groups. This method uses optimisation algorithms to assign weights to individuals such that the weighted average of the covariates in ALSPAC (the larger genotyped cohort) matches that of MCS (the smaller genotyped cohort), minimising confounding biases and increasing comparability. We used the child’s biological sex, maternal age at delivery, and maternal highest educational qualification at first data collection in each cohort as matching factors. We considered using propensity score matching with a 1:1 ratio to obtain a well-balanced subsample of ALSPAC, yet this approach would have resulted in substantial data loss and potential risk of residual confounding due to limited covariate selection. Therefore, we opted for entropy balancing to retain a larger sample size in the ALSPAC cohort. Entropy balancing was conducted using the *ebal* package in R124. #### Association with developmental milestones and autism diagnosis In ALSPAC, we obtained understanding of simple phrases (e.g., “do you want that”, or “come here”) and gesture scores from the Macarthur-Bates Communicative Development Inventories125 at 15 months of age. We conducted multiple linear regression using PGS for iPSYCHbefore11 and iPSYCHafter10, with sex, age, and the first 10 genetic principal components as covariates. Autism diagnosis in the MCS was obtained using parent/caregiver reports of autism/asperger syndrome diagnosis by a doctor at ages 5, 7, 11, and 14. We identified individuals with an autism diagnosis at age 7 or earlier, age 11 or earlier, or between ages 11 and 14. We conducted Firth’s bias-reduced multiple logistic regression (*logistf* package in R) using PGS for iPSYCHbefore11 and iPSYCHafter10, with sex, age and the first 10 genetic principal components covariates. ## Supporting information Supplementary Notes and FAQs [[supplements/311279_file02.pdf]](pending:yes) Supplementary Tables [[supplements/311279_file03.xlsx]](pending:yes) ## Code availability * Lavaan (LGCM): [https://lavaan.ugent.be/tutorial/growth.html](https://lavaan.ugent.be/tutorial/growth.html) * lcmm(GMM): [https://github.com/CecileProust-Lima/lcmm](https://github.com/CecileProust-Lima/lcmm) * Softimpute: [https://cran.r-project.org/web/packages/softImpute/softImpute.pdf](https://cran.r-project.org/web/packages/softImpute/softImpute.pdf) * SPARK quality control, imputation and GWAS: [https://github.com/vwarrier/SPARK\_iWES2\_imputation/](https://github.com/vwarrier/SPARK_iWES2_imputation/) * Bespoke genetic analyses code: [https://github.com/vwarrier/autism\_agediagnosis/](https://github.com/vwarrier/autism_agediagnosis/) * PRScs: [https://github.com/getian107/PRScs](https://github.com/getian107/PRScs) * fastGWA and GCTA: [https://yanglab.westlake.edu.cn/software/gcta/#Overview](https://yanglab.westlake.edu.cn/software/gcta/#Overview) * GenomicSEM:[https://github.com/GenomicSEM/GenomicSEM](https://github.com/GenomicSEM/GenomicSEM) * LDSC: [https://github.com/bulik/ldsc](https://github.com/bulik/ldsc) * KING: [https://www.kingrelatedness.com/manual.shtml](https://www.kingrelatedness.com/manual.shtml) * Plink: [https://www.cog-genomics.org/plink/2.0/](https://www.cog-genomics.org/plink/2.0/) * PC-AiR: [https://github.com/UW-GAC/GENESIS](https://github.com/UW-GAC/GENESIS) * Lme4: [https://github.com/lme4/lme4/](https://github.com/lme4/lme4/) * Logistf: [https://cran.r-project.org/web/packages/logistf/index.html](https://cran.r-project.org/web/packages/logistf/index.html) ## Data availability * SPARK autism GWAS: [https://bitbucket.org/steinlabunc/spark\_asd\_sumstats/src](https://bitbucket.org/steinlabunc/spark_asd_sumstats/src) * Finngen autism GWAS: [https://www.finngen.fi/en/access\_results](https://www.finngen.fi/en/access_results) * iPSYCH autism GWAS (unstratified, sex-stratified and age at diagnosis stratified) can be obtained from Anders Borglum and Jakob Grove. * Psychiatric GWAS summary stats: [https://pgc.unc.edu/](https://pgc.unc.edu/) * GWAS educational attainment: [https://thessgac.com/](https://thessgac.com/) * GWAS cognitive aptitude: [https://cncr.nl/research/summary\_statistics/](https://cncr.nl/research/summary_statistics/) * For ALSPAC, the study website contains details of all the data that is available through a fully searchable data dictionary and variable search tool“: [http://www.bristol.ac.uk/alspac/researchers/our-data/](http://www.bristol.ac.uk/alspac/researchers/our-data/) * For MCS, data can be obtain after application through the UK Data Service: [https://beta.ukdataservice.ac.uk/datacatalogue/series/series?id=2000031](https://beta.ukdataservice.ac.uk/datacatalogue/series/series?id=2000031) * Summary statistics for the SPARK based GWAS will be made upon publication. ## Ethics declarations ADB received speakers’ fee from Lundbeck. The authors declare no competing interests. ## Author contributions XZ and VW conducted most of the analyses, with the remainder being conducted by JG, YG, CKB, and LKN. MK, EV, AG, LH, AH, AR, BSt.P, and ADB provided summary statistics for various analyses. SASN, DSM, and EMW carried out data preparation and quality control, with assistance from KES, VKC, PD, SL, TM, MK, SA and DB. VW supervised the analyses and directed the study with inputs from HCM and JG. VW and XZ wrote the initial draft with input from HCM. DSM, EW, DHG, NRW, EBR, TB, and SBC provided intellectual input. All authors read and commented on the final manuscript. ## APEX Consortium Deep Adhya, Carrie Allison, Bonnie Ayeung, Rosie Bamford, Simon Baron-Cohen, Richard Bethlehem, Tal Biron-Shental, Graham Burton, Wendy Cowell, Jonathan Davies, Joanna Davis, Dori Floris, Alice Franklin, Lidia Gabis, Daniel Geschwind, Ramin Ali Marandi Ghoddousi, David M. Greenberg, Yuanjun Gu, Alexandra Havdahl, Alexander Heazell, Rosemary Holt, Matthew Hurles, Yumnah Khan, Meng-Chuan Lai, Madeline Lancaster, Michael Lombardo, Hilary Martin, Jose Gonzalez Martinez, Jonathan Mill, Mahmoud Musa, Kathy Niakan, Adam Pavlinek, Lucia Dutan Polit, Marcin Radecki, David Rowitch, Jenifer Sakai, Laura Sichlinger, Deepak Srivastava, Alexandros Tsompanidis, Florina Uzefovsky, Varun Warrier, Elizabeth Weir, Xinhe Zhang. ## iPSYCH Autism working group Anders Borglum, Jonas Bybjerg-Grauholm, Jakob Grove, David M. Hougaard, Ole Mors, Preben Bo Mortensen, Merete Nordentoft and Thomas Werge. ## PGC-PTSD consortium Caroline M. Nievergelt, Adam X. Maihofer, Elizabeth G. Atkinson, Chia-Yen Chen, Karmel W. Choi, Jonathan R. I. Coleman, Nikolaos P. Daskalakis, Laramie E. Duncan, Renato Polimanti, Cindy Aaronson, Ananda B. Amstadter, Soren B. Andersen, Ole A. Andreassen, Paul A. Arbisi, Allison E. Ashley-Koch, S. Bryn Austin, Esmina Avdibegoviç, Dragan Babić, Silviu-Alin Bacanu, Dewleen G. Baker, Anthony Batzler, Jean C. Beckham, Sintia Belangero, Corina Benjet, Carisa Bergner, Linda M. Bierer, Joanna M. Biernacka, Laura J. Bierut, Jonathan I. Bisson, Marco P. Boks, Elizabeth A. Bolger, Amber Brandolino, Gerome Breen, Rodrigo Affonseca Bressan, Richard A. Bryant, Angela C. Bustamante, Jonas Bybjerg-Grauholm, Marie Bækvad-Hansen, Anders D. Børglum, Sigrid Børte, Leah Cahn, Joseph R. Calabrese, Jose Miguel Caldas-de-Almeida, Chris Chatzinakos, Sheraz Cheema, Sean A. P. Clouston, Lucía Colodro-Conde, Brandon J. Coombes, Carlos S. Cruz-Fuentes, Anders M. Dale, Shareefa Dalvie, Lea K. Davis, Jürgen Deckert, Douglas L. Delahanty, Michelle F. Dennis, Frank Desarnaud, Christopher P. DiPietro, Seth G. Disner, Anna R. Docherty, Katharina Domschke, Grete Dyb, Alma Džubur Kulenović, Howard J. Edenberg, Alexandra Evans, Chiara Fabbri, Negar Fani, Lindsay A. Farrer, Adriana Feder, Norah C. Feeny, Janine D. Flory, David Forbes, Carol E. Franz, Sandro Galea, Melanie E. Garrett, Bizu Gelaye, Joel Gelernter, Elbert Geuze, Charles F. Gillespie, Slavina B. Goleva, Scott D. Gordon, Aferdita Goçi, Lana Ruvolo Grasser, Camila Guindalini, Magali Haas, Saskia Hagenaars, Michael A. Hauser, Andrew C. Heath, Sian M. J. Hemmings, Victor Hesselbrock, Ian B. Hickie, Kelleigh Hogan, David Michael Hougaard, Hailiang Huang, Laura M. Huckins, Kristian Hveem, Miro Jakovljević, Arash Javanbakht, Gregory D. Jenkins, Jessica Johnson, Ian Jones, Tanja Jovanovic, Karen-Inge Karstoft, Milissa L. Kaufman, James L. Kennedy, Ronald C. Kessler, Alaptagin Khan, Nathan A. Kimbrel, Anthony P. King, Nastassja Koen, Roman Kotov, Henry R. Kranzler, Kristi Krebs, William S. Kremen, Pei-Fen Kuan, Bruce R. Lawford, Lauren A. M. Lebois, Kelli Lehto, Daniel F. Levey, Catrin Lewis, Israel Liberzon, Sarah D. Linnstaedt, Mark W. Logue, Adriana Lori, Yi Lu, Benjamin J. Luft, Michelle K. Lupton, Jurjen J. Luykx, Iouri Makotkine, Jessica L. Maples-Keller, Shelby Marchese, Charles Marmar, Nicholas G. Martin, Gabriela A. Martínez-Levy, Kerrie McAloney, Alexander McFarlane, Katie A. McLaughlin, Samuel A. McLean, Sarah E. Medland, Divya Mehta, Jacquelyn Meyers, Vasiliki Michopoulos, Elizabeth A. Mikita, Lili Milani, William Milberg, Mark W. Miller, Rajendra A. Morey, Charles Phillip Morris, Ole Mors, Preben Bo Mortensen, Mary S. Mufford, Elliot C. Nelson, Merete Nordentoft, Sonya B. Norman, Nicole R. Nugent, Meaghan O’Donnell, Holly K. Orcutt, Pedro M. Pan, Matthew S. Panizzon, Gita A. Pathak, Edward S. Peters, Alan L. Peterson, Matthew Peverill, Robert H. Pietrzak, Melissa A. Polusny, Bernice Porjesz, Abigail Powers, Xue-Jun Qin, Andrew Ratanatharathorn, Victoria B. Risbrough, Andrea L. Roberts, Alex O. Rothbaum, Barbara O. Rothbaum, Peter Roy-Byrne, Kenneth J. Ruggiero, Ariane Rung, Heiko Runz, Bart P. F. Rutten, Stacey Saenz de Viteri, Giovanni Abrahão Salum, Laura Sampson, Sixto E. Sanchez, Marcos Santoro, Carina Seah, Soraya Seedat, Julia S. Seng, Andrey Shabalin, Christina M. Sheerin, Derrick Silove, Alicia K. Smith, Jordan W. Smoller, Scott R. Sponheim, Dan J. Stein, Synne Stensland, Jennifer S. Stevens, Jennifer A. Sumner, Martin H. Teicher, Wesley K. Thompson, Arun K. Tiwari, Edward Trapido, Monica Uddin, Robert J. Ursano, Unnur Valdimarsdóttir, Miranda Van Hooff, Eric Vermetten, Christiaan H. Vinkers, Joanne Voisey, Yunpeng Wang, Zhewu Wang, Monika Waszczuk, Heike Weber, Frank R. Wendt, Thomas Werge, Michelle A. Williams, Douglas E. Williamson, Bendik S. Winsvold, Sherry Winternitz, Christiane Wolf, Erika J. Wolf, Yan Xia, Ying Xiong, Rachel Yehuda, Keith A. Young, Ross McD Young, Clement C. Zai, Gwyneth C. Zai, Mark Zervas, Hongyu Zhao, Lori A. Zoellner, John-Anker Zwart, Terri deRoon-Cassini, Sanne J. H. van Rooij, Leigh L. van den Heuvel, AURORA Study, Estonian Biobank Research Team, FinnGen Investigators, HUNT All-In Psychiatry, Murray B. Stein, Kerry J. Ressler and Karestan C. Koenen ## Acknowledgements This research was supported by funding from the Simons Foundation for Autism Research Initiative, the Wellcome Trust (214322\Z\18\Z), Horizon-Europe R2D2-MH (grant agreement number 101057385), and UKRI (10063472). For the purpose of open access, we have applied a CC BY public copyright licence to any author-accepted manuscript version arising from this submission. S.B.-C. also received funding from the Autism Centre of Excellence at Cambridge, the Templeton World Charitable Fund, the MRC and the National Institute for Health Research Cambridge Biomedical Research Centre. The research was supported by the National Institute for Health Research Applied Research Collaboration East of England. Any views expressed are those of the author(s) and not necessarily those of the funder. Some of the results leading to this publication have received funding from the Innovative Medicines Initiative 2 Joint Undertaking under grant agreement no. 777394 for the project AIMS-2-TRIALS. This joint undertaking receives support from the European Union’s Horizon 2020 research and innovation program and the EFPIA and Autism Speaks, Autistica and the SFARI. The iPSYCH team was supported by grants from the Lundbeck Foundation (R102-A9118, R155-2014-1724 and R248-2017-2003), the NIMH (1R01MH124851-01 to A.D.B.), and EU’s Horizon Europe program (R2D2-MH; grant agreement no. 101057385 to A.D.B.). The Danish National Biobank resource was supported by the Novo Nordisk Foundation. High-performance computer capacity for handling and statistical analysis of iPSYCH data on the GenomeDK HPC facility was provided by the Center for Genomics and Personalized Medicine and the Centre for Integrative Sequencing, iSEQ, Aarhus University, Denmark (grant to A.D.B.). The UK Medical Research Council and Wellcome (Grant ref: 217065/Z/19/Z) and the University of Bristol provide core support for ALSPAC. This publication is the work of the authors and the authors will serve as guarantors for the contents of this paper. A comprehensive list of grants funding is available on the ALSPAC website ([http://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgements.pdf](http://www.bristol.ac.uk/alspac/external/documents/grant-acknowledgements.pdf)). R2D2-MH has been funded by Horizon Europe [grant agreement no. 101057385], by UK Research and Innovation (UKRI) under the UK government’s Horizon Europe funding guarantee [grant no.10039383] and by the Swiss State Secretariat for Education, Research and Innovation (SERI) under contract number 22.00277 We are grateful to the Centre for Longitudinal Studies (CLS), UCL Social Research Institute, for the use of these data and to the UK Data Service for making them available. However, neither CLS nor the UK Data Service bear any responsibility for the analysis or interpretation of these data. This paper uses unit record data from Growing Up in Australia, the Longitudinal Study of Australian Children. The study is conducted in partnership between the Department of Social Services (DSS), the Australian Institute of Family Studies (AIFS) and the Australian Bureau of Statistics (ABS). The findings and views reported in this paper are those of the author and should not be attributed to DSS, AIFS or the ABS. Growing Up in Ireland (GUI) has been funded by the Government of Ireland through the Department of Children, Equality, Disability, Integration and Youth (DCEDIY) and the Central Statistics Office (CSO). Results in this report are based on analysis of data from Research Microdata Files provided by the Central Statistics Office (CSO). Neither the CSO nor the DCEDIY take any responsibility for the views expressed or the outputs generated from these analyses. We are extremely grateful to all the families who took part in this study, the midwives for their help in recruiting them, and the whole ALSPAC team, which includes interviewers, computer and laboratory technicians, clerical workers, research scientists, volunteers, managers, receptionists and nurses. We thank Alex Kwong, Tamsin Ford, Will Mandy, and Andrew Grotzinger for helpful discussions. * Received July 31, 2024. * Revision received July 31, 2024. * Accepted August 2, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.Lord, C. et al. Autism spectrum disorder. Nat Rev Dis Primers 6, 5 (2020). 2. 2.Kanner, L. Autistic disturbances of affective contact. Nervous Child: Journal of Psychopathology, Psychotherapy, Mental Hygiene, and Guidance of the Child 2 217–50 (1943). 3. 3.Asperger, H. ‘Autistic psychopathy’ in childhood. in Autism and Asperger syndrome (ed. Frith, U.) 37–92 (Cambridge University Press, Cambridge, 1944). 4. 4.Schendel, D. E. & Thorsteinsson, E. Cumulative Incidence of Autism Into Adulthood for Birth Cohorts in Denmark, 1980-2012. JAMA 320, 1811–1813 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.2018.11328&link_type=DOI) 5. 5.Russell, G. et al. Time trends in autism diagnosis over 20 years: a UK population-based cohort study. J. Child Psychol. Psychiatry 63, 674–682 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/jcpp.13505&link_type=DOI) 6. 6.Jensen, C. M., Steinhausen, H.-C. & Lauritsen, M. B. Time trends over 16 years in incidence-rates of autism spectrum disorders across the lifespan based on nationwide Danish register data. J. Autism Dev. Disord. 44, 1808–1818 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10803-014-2053-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24554161&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 7. 7.Ozonoff, S. et al. Diagnosis of Autism Spectrum Disorder After Age 5 in Children Evaluated Longitudinally Since Infancy. J. Am. Acad. Child Adolesc. Psychiatry 57, 849–857.e2 (2018). 8. 8.Mandy, W. et al. Mental health and social difficulties of late-diagnosed autistic children, across childhood and adolescence. J. Child Psychol. Psychiatry 63, 1405–1414 (2022). 9. 9.May, T., Brignell, A. & Williams, K. Parent-reported autism diagnostic stability and trajectories in the Longitudinal Study of Australian Children. Autism Res. 14, 773–786 (2021). 10. 10.Allison, C. et al. Quantitative Checklist for Autism in Toddlers (Q-CHAT). A population screening study with follow-up: the case for multiple time-point screening for autism. BMJ Paediatr Open 5, e000700 (2021). 11. 11.Davidovitch, M., Levit-Binnun, N., Golan, D. & Manning-Courtney, P. Late diagnosis of autism spectrum disorder after initial negative assessment by a multidisciplinary team. J. Dev. Behav. Pediatr. 36, 227–234 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/DBP.0000000000000133&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25651066&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 12. 12.Avlund, S. H. et al. Factors Associated with a Delayed Autism Spectrum Disorder Diagnosis in Children Previously Assessed on Suspicion of Autism. J. Autism Dev. Disord. 51, 3843–3856 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10803-020-04849-x&link_type=DOI) 13. 13.Bazelmans, T. et al. Mid-childhood autism sibling recurrence in infants with a family history of autism. Autism Res. (2024) doi:10.1002/aur.3182. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/aur.3182&link_type=DOI) 14. 14.Landa, R. J., Reetzke, R., Holingue, C. B., Herman, D. & Hess, C. R. Diagnostic Stability and Phenotypic Differences Among School-Age Children Diagnosed With ASD Before Age 2. Front. Psychiatry 13, 805686 (2022). 15. 15.Austerberry, C., Mateen, M., Fearon, P. & Ronald, A. Heritability of Psychological Traits and Developmental Milestones in Infancy: A Systematic Review and Meta-analysis. JAMA Netw Open 5, e2227887 (2022). 16. 16.Hegemann, L. et al. Genetic and phenotypic heterogeneity in early neurodevelopmental traits in the Norwegian Mother, Father and Child Cohort Study. Mol. Autism 15, 25 (2024). 17. 17. St Pourcain, B., et al. Variability in the common genetic architecture of social-communication spectrum phenotypes during childhood and adolescence. Mol. Autism 5, 18 (2014). 18. 18.Rødgaard, E.-M., Jensen, K., Miskowiak, K. W. & Mottron, L. Autism comorbidities show elevated female-to-male odds ratios and are associated with the age of first autism diagnosis. Acta Psychiatr. Scand. 144, 475–486 (2021). 19. 19.Jadav, N. & Bal, V. H. Associations between co-occurring conditions and age of autism diagnosis: Implications for mental health training and adult autism research. Autism Res. 15, 2112–2125 (2022). 20. 20.Goodman, R. Strengths and Difficulties Questionnaire. Child and Adolescent Psychiatry and Mental HealthEuropean Journal of Psychological AssessmentPsychological AssessmentSchool Psychology QuarterlyClinical Psychologist doi:10.1037/t00540-000. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1037/t00540-000&link_type=DOI) 21. 21.Stone, L. L. et al. The Strengths and Difficulties Questionnaire: psychometric properties of the parent and teacher version in children aged 4-7. BMC Psychol 3, 4 (2015). 22. 22.Kovacs, S. & Sharp, C. Criterion validity of the Strengths and Difficulties Questionnaire (SDQ) with inpatient adolescents. Psychiatry Res. 219, 651–657 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.psychres.2014.06.019&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25048754&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 23. 23.Borg, A.-M., Kaukonen, P., Salmelin, R., Joukamaa, M. & Tamminen, T. Reliability of the strengths and difficulties questionnaire among Finnish 4-9-year-old children. Nord. J. Psychiatry 66, 403–413 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3109/08039488.2012.660706&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22397524&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000310078400006&link_type=ISI) 24. 24.Speyer, L. G., Auyeung, B. & Murray, A. L. Longitudinal Invariance of the Strengths and Difficulties Questionnaire Across Ages 4 to 16 in the ALSPAC Sample. Assessment 30, 1884–1894 (2023). 25. 25.Murray, A. L., Speyer, L. G., Hall, H. A., Valdebenito, S. & Hughes, C. A Longitudinal and Gender Invariance Analysis of the Strengths and Difficulties Questionnaire Across Ages 3, 5, 7, 11, 14, and 17 in a Large U.K.-Representative Sample. Assessment 29, 1248–1261 (2022). 26. 26.Woerner, W. et al. The Strengths and Difficulties Questionnaire overseas: Evaluations and applications of the SDQ beyond Europe. Eur. Child Adolesc. Psychiatry 13, ii47–ii54 (2004). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00787-004-2008-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15243786&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 27. 27.Dalsgaard, S. et al. Incidence Rates and Cumulative Incidences of the Full Spectrum of Diagnosed Mental Disorders in Childhood and Adolescence. JAMA Psychiatry 77, 155–164 (2020). 28. 28.Dahl, R. E., Allen, N. B., Wilbrecht, L. & Suleiman, A. B. Importance of investing in adolescence from a developmental science perspective. Nature 554, 441–450 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature25770&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29469094&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 29. 29.Foulkes, L. & Blakemore, S.-J. Studying individual differences in human adolescent brain development. Nat. Neurosci. 21, 315–323 (2018). 30. 30.Angold, A., Costello, E., Messer, S. & Pickles, A. Development of a short questionnaire for use in epidemiological studies of depression in children and adolescents: Factor composition and structure across development. Int. J. Methods Psychiatr. Res. 5, 251–262 (1995). [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995TQ48000003&link_type=ISI) 31. 31.Hrdlicka, M. et al. Predictors of age at diagnosis in autism spectrum disorders: the use of multiple regression analyses and a classification tree on a clinical sample. Eur. Child Adolesc. Psychiatry 33, 1171–1177 (2024). 32. 32.Brett, D., Warnell, F., McConachie, H. & Parr, J. R. Factors Affecting Age at ASD Diagnosis in UK: No Evidence that Diagnosis Age has Decreased Between 2004 and 2014. Journal of Autism and Developmental Disorders vol. 46 1974–1984 Preprint at doi:10.1007/s10803-016-2716-6 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10803-016-2716-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27032954&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 33. 33.St Pourcain, B., et al. Heritability and genome-wide analyses of problematic peer relationships during childhood and adolescence. Hum. Genet. 134, 539–551 (2015). 34. 34.Knafo, A. & Plomin, R. Prosocial behavior from early to middle childhood: genetic and environmental influences on stability and change. Dev. Psychol. 42, 771–786 (2006). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1037/0012-1649.42.5.771&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16953685&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000240530500002&link_type=ISI) 35. 35.Stergiakouli, E. et al. Shared genetic influences between dimensional ASD and ADHD symptoms during child and adolescent development. Mol. Autism 8, 18 (2017). 36. 36.Jami, E. S. et al. Genome-wide Association Meta-analysis of Childhood and Adolescent Internalizing Symptoms. J. Am. Acad. Child Adolesc. Psychiatry 61, 934–945 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaac.2021.11.035&link_type=DOI) 37. 37.St Pourcain, B., et al. ASD and schizophrenia show distinct developmental profiles in common genetic overlap with population-based social communication difficulties. Mol. Psychiatry (2017) doi:10.1038/mp.2016.198. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/mp.2016.198&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28044064&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 38. 38.St Pourcain, B., et al. Developmental Changes Within the Genetic Architecture of Social Communication Behavior: A Multivariate Study of Genetic Variance in Unrelated Individuals. Biol. Psychiatry 83, 598–606 (2018). 39. 39.SPARK Consortium. Electronic address: pfeliciano{at}simonsfoundation.org & SPARK Consortium. SPARK: A US Cohort of 50,000 Families to Accelerate Autism Research. Neuron 97, 488–493 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2018.01.015&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29420931&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 40. 40.Demontis, D. et al. Genome-wide analyses of ADHD identify 27 risk loci, refine the genetic architecture and implicate several cognitive domains. Nat. Genet. 55, 198–208 (2023). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/S41588-022-01285-8&link_type=DOI) 41. 41.Trubetskoy, V. et al. Mapping genomic loci implicates genes and synaptic biology in schizophrenia. Nature 604, 502–508 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/S41586-022-04434-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 42. 42.Fu, J. M. et al. Rare coding variation provides insight into the genetic architecture and phenotypic context of autism. Nat. Genet. 54, 1320–1331 (2022). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35982160&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 43. 43.Satterstrom, F. K. et al. Large-Scale Exome Sequencing Study Implicates Both Developmental and Functional Changes in the Neurobiology of Autism. Cell 180, 568–584.e23 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2019.12.036&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 44. 44.Als, T. D. et al. Depression pathophysiology, risk prediction of recurrence and comorbid psychiatric disorders using genome-wide analyses. Nat. Med. 29, 1832–1844 (2023). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-023-02352-1&link_type=DOI) 45. 45.Lee, J. J. et al. Gene discovery and polygenic prediction from a genome-wide association study of educational attainment in 1.1 million individuals. Nat. Genet. 50, 1112–1121 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0147-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30038396&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 46. 46.Savage, J. E. et al. Genome-wide association meta-analysis in 269,867 individuals identifies new genetic and functional links to intelligence. Nat. Genet. 1 (2018). 47. 47.Sada-Fuente, E. et al. Common genetic variants contribute to heritability of age at onset of schizophrenia. Transl. Psychiatry 13, 201 (2023). 48. 48.Harder, A. et al. Genetics of age-at-onset in major depression. Transl. Psychiatry 12, 124 (2022). 49. 49.49. Feng, Y.-C. A., et al. Findings and insights from the genetic investigation of age of first reported occurrence for complex disorders in the UK Biobank and FinnGen. bioRxiv (2020) doi:10.1101/2020.11.20.20234302. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1101/2020.11.20.20234302&link_type=DOI) 50. 50.Warrier, V. et al. Genetic correlates of phenotypic heterogeneity in autism. Nat. Genet. (2022) doi:10.1038/s41588-022-01072-5. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-022-01072-5&link_type=DOI) 51. 51.The Autism Spectrum Disorders Working Group of The Psychiatric Genomics Consortium. Meta-analysis of GWAS of over 16,000 individuals with autism spectrum disorder highlights a novel locus at 10q24.32 and a significant overlap with schizophrenia. Mol. Autism 8, 21 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13229-017-0137-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28540026&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 52. 52.Matoba, N. et al. Common genetic risk variants identified in the SPARK cohort support DDHD2 as a candidate risk gene for autism. Transl. Psychiatry 10, 265 (2020). 53. 53.Patalay, P. & Gage, S. H. Changes in millennial adolescent mental health and health-related behaviours over 10 years: a population cohort comparison study. Int. J. Epidemiol. 48, 1650–1664 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyz006&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 54. 54.Grotzinger, A. D. et al. Genomic structural equation modelling provides insights into the multivariate genetic architecture of complex traits. Nature Human Behaviour 1 (2019). 55. 55.Dawson, G., Rieder, A. D. & Johnson, M. H. Prediction of autism in infants: progress and challenges. Lancet Neurol. 22, 244–254 (2023). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s1474-4422(22)00407-0&link_type=DOI) 56. 56.Wray, N. R. et al. Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nat. Genet. 50, 668–681 (2018). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0090-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29700475&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 57. 57.Sheldrick, R. C., Maye, M. P. & Carter, A. S. Age at First Identification of Autism Spectrum Disorder: An Analysis of Two US Surveys. J. Am. Acad. Child Adolesc. Psychiatry 56, 313–320 (2017). 58. 58.Bulik-Sullivan, B. K. et al. An atlas of genetic correlations across human diseases and traits. Nat. Genet. 47, 1236–1241 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3406&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26414676&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 59. 59.Cross-Disorder Group of the Psychiatric Genomics Consortium et al. Genetic relationship between five psychiatric disorders estimated from genome-wide SNPs. Nat. Genet. 45, 984–994 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2711&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23933821&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 60. 60.Grove, J. et al. Identification of common genetic risk variants for autism spectrum disorder. Nat. Genet. 51, 431–444 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0344-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30804558&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 61. 61.Rajagopal, V. M. et al. Differences in the genetic architecture of common and rare variants in childhood, persistent and late-diagnosed attention-deficit hyperactivity disorder. Nat. Genet. 54, 1117–1124 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-022-01143-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=35927488&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 62. 62.Breunig, S. et al. Examining Differences in the Genetic and Functional Architecture of Attention-Deficit/Hyperactivity Disorder Diagnosed in Childhood and Adulthood. Biol Psychiatry Glob Open Sci 4, 100307 (2024). 63. 63.Martini, M. I. et al. Sex Differences in Mental Health Problems and Psychiatric Hospitalization in Autistic Young Adults. JAMA Psychiatry 79, 1188–1198 (2022). 64. 64.Dworzynski, K., Ronald, A., Bolton, P. & Happé, F. How different are girls and boys above and below the diagnostic threshold for autism spectrum disorders? J. Am. Acad. Child Adolesc. Psychiatry 51, 788–797 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jaac.2012.05.018&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22840550&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 65. 65.Milner, V. et al. Does camouflaging predict age at autism diagnosis? A comparison of autistic men and women. Autism Res. 17, 626–636 (2024). 66. 66.Dean, M., Harwood, R. & Kasari, C. The art of camouflage: Gender differences in the social behaviors of girls and boys with autism spectrum disorder. Autism 21, 678–689 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1362361316671845&link_type=DOI) 67. 67.Hull, L. et al. Development and Validation of the Camouflaging Autistic Traits Questionnaire (CAT-Q). J. Autism Dev. Disord. 49, 819–833 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10803-018-3792-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 68. 68.Cook, J., Crane, L., Hull, L., Bourne, L. & Mandy, W. Self-reported camouflaging behaviours used by autistic adults during everyday social interactions. Autism 26, 406–421 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/13623613211026754&link_type=DOI) 69. 69.Milner, V., Mandy, W., Happé, F. & Colvert, E. Sex differences in predictors and outcomes of camouflaging: Comparing diagnosed autistic, high autistic trait and low autistic trait young adults. Autism 27, 402–414 (2023). 70. 70.Ross, A., Grove, R. & McAloon, J. The relationship between camouflaging and mental health in autistic children and adolescents. Autism Res. 16, 190–199 (2023). 71. 71.Perry, E., Mandy, W., Hull, L. & Cage, E. Understanding Camouflaging as a Response to Autism-Related Stigma: A Social Identity Theory Approach. J. Autism Dev. Disord. 52, 800–810 (2022). 72. 72.Roman-Urrestarazu, A., Dumas, G. & Warrier, V. Naming Autism in the Right Context. JAMA Pediatr. 176, 633–634 (2022). 73. 73.Connelly, R. & Platt, L. Cohort profile: UK Millennium Cohort Study (MCS). Int. J. Epidemiol. 43, 1719–1725 (2014). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyu001&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24550246&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 74. 74.Layte, R. & McCrory, C. Growing Up in Ireland: Maternal Health Behaviours and Child Growth in Infancy. Research Series (2015). 75. 75.Clifford, S. A., Davies, S., Wake, M. & Child Health CheckPoint Team. Child Health CheckPoint: cohort summary and methodology of a physical health and biospecimen module for the Longitudinal Study of Australian Children. BMJ Open 9, 3–22 (2019). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYm1qb3BlbiI7czo1OiJyZXNpZCI7czoxMToiOS9TdXBwbF8zLzMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOC8wMi8yMDI0LjA3LjMxLjI0MzExMjc5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 76. 76.Anderson, E. R. Analyzing change in short-term longitudinal research using cohort-sequential designs. J. Consult. Clin. Psychol. 61, 929–940 (1993). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1037/0022-006X.61.6.929&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8113494&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1993MN43800004&link_type=ISI) 77. 77.Russell, G., Rodgers, L. R. & Ford, T. The strengths and difficulties questionnaire as a predictor of parent-reported diagnosis of autism spectrum disorder and attention deficit hyperactivity disorder. PLoS One 8, e80247 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0080247&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24312466&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 78. 78.Grasso, M., Lazzaro, G., Demaria, F., Menghini, D. & Vicari, S. The Strengths and Difficulties Questionnaire as a Valuable Screening Tool for Identifying Core Symptoms and Behavioural and Emotional Problems in Children with Neuropsychiatric Disorders. Int. J. Environ. Res. Public Health 19, (2022). 79. 79.Roman-Urrestarazu, A. et al. Association of Race/Ethnicity and Social Disadvantage With Autism Prevalence in 7 Million School Children in England. JAMA Pediatr. 175, e210054 (2021). 80. 80.Terhaag, S., Fitzsimons, E., Daraganova, G. & Patalay, P. Sex, ethnic and socioeconomic inequalities and trajectories in child and adolescent mental health in Australia and the UK: findings from national prospective longitudinal studies. J. Child Psychol. Psychiatry 62, 1255–1267 (2021). 81. 81.Pender, R., Fearon, P., St Pourcain, B., Heron, J. & Mandy, W. Developmental trajectories of autistic social traits in the general population. Psychol. Med. 1–9 (2021). 82. 82.Rosseel, Y. Lavaan: An R package for structural equation modeling and more. Version 0.5--12 (BETA). J. Stat. Softw. 48, 1–36 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v048.i04&link_type=DOI) 83. 83.Hoeksma, J. B. & Kelderman, H. On growth curves and mixture models. Infant Child Dev. 15, 627–634 (2006). 84. 84.Proust-Lima, C., Philipps, V. & Liquet, B. Estimation of Extended Mixed Models Using Latent Classes and Latent Processes: The R Package lcmm. J. Stat. Softw. 78, 1–56 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v078.i02&link_type=DOI) 85. 85.Ram, N. & Grimm, K. J. Growth Mixture Modeling: A Method for Identifying Differences in Longitudinal Change Among Unobserved Groups. Int. J. Behav. Dev. 33, 565–576 (2009). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/0165025409343765&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23885133&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000271142100009&link_type=ISI) 86. 86.Grimm, K. J., Mazza, G. L. & Davoudzadeh, P. Model Selection in Finite Mixture Models: A k-Fold Cross-Validation Approach. Struct. Equ. Modeling 24, 246–256 (2017). 87. 87.Tonidandel, S. & LeBreton, J. M. Relative Importance Analysis: A Useful Supplement to Regression Analysis. J. Bus. Psychol. 26, 1–9 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10869-010-9204-3&link_type=DOI) 88. 88.Leaf, J. B. et al. Increasing social skills and pro-social behavior for three children diagnosed with autism through the use of a teaching package. Res. Autism Spectr. Disord. 3, 275–289 (2009). 89. 89.Azen, R. & Budescu, D. V. Comparing Predictors in Multivariate Regression Models: An Extension of Dominance Analysis. J. Educ. Behav. Stat. 31, 157–180 (2006). 90. 90.Manichaikul, A. et al. Robust relationship inference in genome-wide association studies. Bioinformatics 26, 2867–2873 (2010). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btq559&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20926424&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000283919800010&link_type=ISI) 91. 91.Conomos, M. P., Miller, M. B. & Thornton, T. A. Robust inference of population structure for ancestry prediction and correction of stratification in the presence of relatedness. Genet. Epidemiol. 39, 276–293 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/gepi.21896&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25810074&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 92. 92.Taliun, D. et al. Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program. Nature 590, 290–299 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-021-03205-y&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33568819&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 93. 93.Das, S. et al. Next-generation genotype imputation service and methods. Nat. Genet. 48, 1284–1287 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3656&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27571263&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 94. 94.Loh, P.-R., Palamara, P. F. & Price, A. L. Fast and accurate long-range phasing in a UK Biobank cohort. Nat. Genet. 48, 811–816 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3571.&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27270109&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 95. 95.Ge, T., Chen, C.-Y., Ni, Y., Feng, Y.-C. A. & Smoller, J. W. Polygenic prediction via Bayesian regression and continuous shrinkage priors. Nat. Commun. 10, 1776 (2019). 96. 96.Bybjerg-Grauholm, J. et al. The iPSYCH2015 Case-Cohort sample: updated directions for unravelling genetic and environmental architectures of severe mental disorders. Preprint at doi:10.1101/2020.11.30.20237768. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4xMS4zMC4yMDIzNzc2OHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDgvMDIvMjAyNC4wNy4zMS4yNDMxMTI3OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 97. 97.Mullins, N. et al. Genome-wide association study of more than 40,000 bipolar disorder cases provides new insights into the underlying biology. Nat. Genet. 53, 817–829 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-021-00857-4&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34002096&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 98. 98.Okbay, A. et al. Polygenic prediction of educational attainment within and between families from genome-wide association analyses in 3 million individuals. Nat. Genet. 54, 437–449 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588--01016-z&link_type=DOI) 99. 99.Dsm-Iv. Diagnostic and Statistical Manual of Mental Disorders, Fourth *Edition*. (American Psychiatric Publishing, Inc., Washington, DC, 1994). 100.100.Purcell, S. et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575 (2007). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/519795&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17701901&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 101.101.Chang, C. C. et al. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 4, 7 (2015). 102.102.Demange, P. A. et al. Estimating effects of parents’ cognitive and non-cognitive skills on offspring education using polygenic scores. Nat. Commun. 13, 4801 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-022-32003-x&link_type=DOI) 103.103.Weiner, D. J. et al. Polygenic transmission disequilibrium confirms that common and rare variation act additively to create risk for autism spectrum disorders. Nat. Genet. 49, 978–985 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3863&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28504703&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 104.104.Koko, M. E. et al. Contribution of autosomal rare and de novo variants to sex differences in autism. medRxiv 2024–2004 (2024). 105.105.Karczewski, K. J. et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 581, 434–443 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2308-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32461654&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 106.106.Samocha, K. E., et al. Regional missense constraint improves variant deleteriousness prediction. bioRxiv 148353 (2017) doi:10.1101/148353. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czo4OiIxNDgzNTN2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA4LzAyLzIwMjQuMDcuMzEuMjQzMTEyNzkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 107.107.Jiang, L. et al. A resource-efficient tool for mixed model association analysis of large-scale data. Nat. Genet. 51, 1749–1755 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588–019–0530–8&link_type=DOI) 108.108.Lam, M. et al. RICOPILI: Rapid Imputation for COnsortias PIpeLIne. Bioinformatics 36, 930–933 (2020). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 109.109.Howie, B., Fuchsberger, C., Stephens, M., Marchini, J. & Abecasis, G. R. Fast and accurate genotype imputation in genome-wide association studies through pre-phasing. Nat. Genet. 44, 955–959 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2354&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22820512&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 110.110.Iglesias, A. I. et al. Haplotype reference consortium panel: Practical implications of imputations with large reference panels. Hum. Mutat. 38, 1025–1032 (2017). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/humu.23247&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 111.111.Jiang, L., Zheng, Z., Fang, H. & Yang, J. A generalized linear mixed model association tool for biobank-scale data. Nat. Genet. 53, 1616–1621 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-021-00954-4&link_type=DOI) 112.112.Peyrot, W. J., Boomsma, D. I., Penninx, B. W. J. H. & Wray, N. R. Disease and polygenic architecture: Avoid trio design and appropriately account for unscreened control subjects for common disease. Am. J. Hum. Genet. 98, 382–391 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2015.12.017&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26849113&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 113.113.Yang, J., Lee, S. H., Goddard, M. E. & Visscher, P. M. GCTA: a tool for genome-wide complex trait analysis. Am. J. Hum. Genet. 88, 76–82 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2010.11.011&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21167468&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 114.114.Yang, J. et al. Common SNPs explain a large proportion of the heritability for human height. Nat. Genet. 42, 565–569 (2010). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.608&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20562875&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000279242400007&link_type=ISI) 115.115.Bulik-Sullivan, B. K. et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat. Genet. 47, 291–295 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3211&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25642630&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 116.116.Kurki, M. I. et al. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature 613, 508–518 (2023). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-02205473-8&link_type=DOI) 117.117.Golding, J., Pembrey, M., Jones, R. & ALSPAC Study Team. ALSPAC--the Avon Longitudinal Study of Parents and Children. I. Study methodology. Paediatr. Perinat. Epidemiol. 15, 74–87 (2001). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1365-3016.2001.00325.x&link_type=DOI) 118.118.Boyd, A. et al. Cohort Profile: The ‘Children of the 90s’—the index offspring of the Avon Longitudinal Study of Parents and Children. Int. J. Epidemiol. 42, 111–127 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dys064&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22507743&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000316699300012&link_type=ISI) 119.119.Fraser, A. et al. Cohort Profile: the Avon Longitudinal Study of Parents and Children: ALSPAC mothers cohort. Int. J. Epidemiol. 42, 97–110 (2013). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dys066&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22507742&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000316699300011&link_type=ISI) 120.120.Fitzsimons, E. et al. Collection of genetic data at scale for a nationally representative population: the UK Millennium Cohort Study. Longit. Life Course Stud. 13, 169–187 (2021). 121.121.Pain, O. et al. Evaluation of polygenic prediction methodology within a reference-standardized framework. PLoS Genet. 17, e1009021 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1009021&link_type=DOI) 122.122.Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting Linear Mixed-Effects Models Using lme4. J. Stat. Softw. 67, 1–48 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.18637/jss.v067.i0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24227677&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F02%2F2024.07.31.24311279.atom) 123.123.Hainmueller, J. Entropy Balancing for Causal Effects: A Multivariate Reweighting Method to Produce Balanced Samples in Observational Studies. Polit. Anal. 20, 25–46 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/pan/mpr025&link_type=DOI) 124.124.ebal package - RDocumentation. [https://rdocumentation.org/packages/ebal/versions/0.1-8](https://rdocumentation.org/packages/ebal/versions/0.1-8). 125.125.Fenson, L. et al. MacArthur-Bates Communicative Development Inventories, Second Edition. doi:10.1037/t11538-000. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1037/t11538-000&link_type=DOI)