# Social and non-social autism symptoms and trait domains are genetically dissociable

Varun Warrier, Roberto Toro, Hyejung Won, Claire S. Leblond, Freddy Cliquet, Richard Delorme, Ward de Witte, Janita Bralten, Bhismadev Chakrabarti, Anders D. Børglum, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Varun Warrier, Roberto Toro, Hyejung Won, Claire S. Leblond, Freddy Cliquet, et al., Social and non-social autism symptoms and trait domains are genetically dissociable. Communications Biology, 2019, 2 (1), pp.328. 10.1038/s42003-019-0558-4. pasteur-03261098

HAL Id: pasteur-03261098

https://pasteur.hal.science/pasteur-03261098

Submitted on 15 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

COMMUNICATIONS
BIOLOGY

ARTICLE

https://doi.org/10.1038/s42003-019-0558-4

OPEN

# Social and non-social autism symptoms and trait domains are genetically dissociable

Varun Warrier<fnref n="1" />, Roberto Toro ID ², Hyejung Won³, Claire S. Leblond², Freddy Cliquet ID ², Richard Delorme²,⁴,
Ward De Witte⁵, Janita Bralten⁵,⁶, Bhismadev Chakrabarti ID <fnref n="1" />,⁷, Anders D. Børglum ID ⁸,⁹,¹⁰,
Jakob Grove ID ⁸,⁹,¹⁰,¹¹, Geert Poelmans⁵, the 23andMe Research Team#, David A. Hinds ID ¹²,¹³,
Thomas Bourgeron²,¹³ & Simon Baron-Cohen<fnref n="1" />,¹³

The core diagnostic criteria for autism comprise two symptom domains – social and communication difficulties, and unusually repetitive and restricted behaviour, interests and activities. There is some evidence to suggest that these two domains are dissociable, though this hypothesis has not yet been tested using molecular genetics. We test this using a genome-wide association study (N = 51,564) of a non-social trait related to autism, systemising, defined as the drive to analyse and build systems. We demonstrate that systemising is heritable and genetically correlated with autism. In contrast, we do not identify significant genetic correlations between social autistic traits and systemising. Supporting this, polygenic scores for systemising are significantly and positively associated with restricted and repetitive behaviour but not with social difficulties in autistic individuals. These findings strongly suggest that the two core domains of autism are genetically dissociable, and point at how to fractionate the genetics of autism.

<footnote n="1">¹Autism Research Centre, Department of Psychiatry, University of Cambridge, Cambridgeshire, UK. ²Human Genetics and Cognitive Functions, Institut Pasteur, UMR3571 CNRS, Université de Paris, Paris, France. ³Department of Genetics and Neuroscience Center, University of North Carolina, Chapel Hill, NC 27599, USA. ⁴Child and Adolescent Psychiatry Department, Robert Debré Hospital, Paris, France. ⁵Department of Human Genetics, Radboud University Medical Center, Nijmegen, The Netherlands. ⁶Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, The Netherlands. ⁷Centre for Autism, School of Psychology and Clinical Language Sciences, University of Reading, Reading, UK. ⁸The Lundbeck Foundation Initiative for Integrative Psychiatric Research, iPSYCH, Aarhus, Denmark. ⁹Centre for Integrative Sequencing, iSEQ, Aarhus University, Aarhus, Denmark. ¹⁰Department of Biomedicine - Human Genetics, Aarhus University, Aarhus, Denmark. ¹¹Bioinformatics Research Centre, Aarhus University, Aarhus, Denmark. ¹²23andMe Inc., Mountain View, CA 94043, USA. ¹³These authors contributed equally: David A. Hinds, Thomas Bourgeron, Simon Baron-Cohen. #Members of the 23andMe Research Team are listed at the end of the paper. Correspondence and requests for materials should be addressed to V.W. (email: vw260@medschl.cam.ac.uk) or to S.B.-C. (email: sb205@cam.ac.uk)</footnote>
The core diagnostic criteria of autism comprises two symptom domains: difficulties in social interactions and communication (the social domain), and unusually repetitive and restricted behaviour and stereotyped interests (the non-social domain)¹. Multiple lines of evidence suggest that these two domains are dissociable²,³. First, factor and principal component analysis of autism and autistic traits have predominantly identified two factors—a social and a non-social factor⁴⁻⁹. Second, investigations of autistic traits in large cohorts have demonstrated a positive phenotypic correlation between different social traits and different non-social traits separately, but only a limited correlation between social and non-social traits⁹⁻¹². Third, twin genetic correlations between social and non-social symptom domains in autism are low, although both social and non-social trait domains are highly heritable in neurotypical¹³,¹⁴ or autistic twins¹⁵. Fourth, difficulties in social and non-social domains can occur independently of each other¹⁶,¹⁷, which has been used to subgroup individuals on the spectrum based on the two domains¹⁸. This suggests that the genetic and phenotypic architecture of autism consists of at least two broadly dissociable domains. This has implications for genetic, biological, and clinical studies of autism, since most studies have investigated autism as if it is a unitary condition³. The idea that social and non-social symptom domains are dissociable is unsurprising given their very different nature, and very different underlying neurology and cognitive processes: one related to interpreting animate motion and mental states (theory of mind) and the other related to recognising inanimate objects, events or patterns (systemising)³. Nevertheless, a diagnosis of autism is only given when the social and non-social symptom domains cluster together.

However, to date, there has been limited molecular genetic evidence in support of this dissociability hypothesis, partly due to the limited large-scale research on the genetics of social and non-social domains. Most genetic research into the social and non-social domains has been primarily through linkage and genome-wide association studies (GWAS) in relatively small samples of autistic individuals and the general population (*N* < 5K)¹⁹⁻²⁵. This has precluded a detailed molecular genetic investigation of the social and non-social domains associated with autism. Given currently available sample sizes with phenotypic information, investigating the genetics of the social and non-social domains in autistic individuals is difficult. However, several studies have demonstrated that the underlying liability for autism is normally distributed in the general population²⁶⁻²⁹. Factor analyses have failed to identify discontinuities between clinical autism and autistic traits in the general population³⁰. Autistic traits are heritable³¹⁻³³, are elevated in family members of autistic individuals compared to the general population³⁴,³⁵, and are transmitted intergenerationally³⁶,³⁷. Factor analysis of autistic traits measures have also identified two different factors in both the general population and autistic individuals—one linked to the social domain, and another linked to the non-social domain, mirroring the factor structure of clinical autism domains⁶,⁹,³⁰,³⁸. Studies have further demonstrated moderate to high shared genetics between the extremes of the liability distribution and the rest of the distribution¹⁴,³⁹⁻⁴¹. One twin study investigated the bivariate genetic correlation between research and clinical autism diagnosis and autistic traits and identified high genetic correlations (0.7 < *r*<sub>g</sub> < 0.89)⁴². Validating this, studies have identified modest shared genetics between autism and autistic traits⁴³⁻⁴⁵. Taken together, there is considerable evidence to suggest that autism represents the extreme end of the autistic traits continuum.

While a few studies have investigated the genetics of traits contributing to the social domains such as social and

communications difficulties¹⁹,⁴⁴,⁴⁵, empathy⁴⁶, and emotion recognition⁴⁷, there have been limited studies investigating the genetics of the non-social domain²⁵,⁴⁸. Neither of these studies have replicably identified significant variants associated with the non-social domain, primarily because of the relatively modest sample sizes of the GWAS. An alternate approach is to investigate the genetics of non-social traits related to autism in the typical population, maximising the sample size. To better understand the genetics of a non-social trait related to autism, we investigate the genetics of systemising measured using a 75-item well validated, self-report measure called the Systemising Quotient-Revised (SQ-R) (see ‘Methods’ section). Systemising involves identifying input-operation-output (or if-and-then) relationships in order to analyse and build systems, and to understand the laws that govern specific systems⁴⁹. The hyper-systemising theory of autism proposes that autistic individuals, on average, have superior attention to detail, and a stronger drive to systemise compared to individuals in the general population⁴⁹. This has been validated in several studies⁵⁰,⁵¹ including a recent study in more 650,000 individuals, including 36,000 autistic individuals¹². Several lines of evidence suggest that autistic individuals have at least intact if not superior systemising. The idea was noted in the earliest papers describing autism by both Asperger⁵² and Kanner⁵³, although these early papers do not use the term ‘systemising’ but instead comment on strong interests in pattern recognition, the need for order and predictability, excellent memory for facts, and a strong focus on objects and understanding how things work. Further, autistic adults, on average, score higher on the SQ-R compared to individuals in the general population¹⁰,⁵¹, a profile also observed in autistic children⁵⁴. Several items in the SQ-R specifically measure circumscribed interests and insistence on sameness, two of the items mentioned in the DSM-5, and several of these items map onto items on the Autism Spectrum Quotient (AQ), a well validated measure of autistic traits²⁷ (see Supplementary Note). Because systems follow rules, they repeat. A fascination with systems may thus manifest as unusually repetitive behaviour. And because systems depend on precise variables, a fascination with systems may also manifest as unusually narrow interests in autism.

The present study has two aims: first, to investigate the genetic architecture of a non-social trait linked to autism (systemising); and second, to investigate whether social and non-social traits related to autism, measured in the general population, are genetically dissociable.

## Results

**GWAS results.** We first conducted a GWAS of systemising (*N* = 51,564) measured using the SQ-R. Following this, and using data from GWAS of social traits genetically correlated with autism (GWAS of self-reported empathy (*N* = 46,861)⁴⁶, and GWAS of social relationship satisfaction⁵⁵ measured using friendship (*N*<sub>effective</sub> = 164,112) and family relationship (*N*<sub>effective</sub> = 158,116) satisfaction scales), we investigated whether the social and non-social traits related to autism are genetically dissociable in the general population. A flow chart of the study design is provided in Fig. 1.

Systemising was measured in the 23andMe sample (*N* = 51,564) using scores from the SQ-R¹⁰. Scores on the SQ-R were normally distributed, with a mean of 71 ± 21 out of 150. As hypothesised based on previous research¹⁰,¹²,⁵¹, males (76.5 ± 20), on average, scored higher than females (65.4 ± 20.6) (*P* < 0.001, Cohen's *d* = 0.54, Supplementary Fig. 1). Given the significant sex differences in scores, we conducted a non-stratified and sex-stratified GWAS for the SQ-R. Genome-wide association analyses identified three significant loci (Fig. 2, Supplementary Data 1 and
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Schematic diagram of the study. We conducted a GWAS of the SQ-R (N = 51,564) and quantified SNP heritability ($h_{SNP}^2$), quantified genetic correlations with multiple phenotypes, and conducted polygenic score analyses. In addition, we conducted sex-stratified GWAS of the SQ-R, and investigated $h_{SNP}^2$ within sex and genetic correlation between males and females. Finally, we investigated the clustering of all phenotypes that are genetically correlated with autism, and whether the social and the non-social phenotypes associated with autism are genetically correlated</figcaption></figure>

Supplementary Fig. 2). Two of these were significant in the non-stratified GWAS: rs4146336 on chromosome 3 ($P = 2.58 \times 10^{-8}$) and rs1559586 on chromosome 18 ($P = 4.78 \times 10^{-8}$). The third significant locus was in the males-only GWAS (rs8005092 on chromosome 14, $P = 3.74 \times 10^{-8}$). rs8005092 and rs1559586 lie in regions of high genetic recombination. Linkage-disequilibrium score regression (LDSR) intercept suggested that there was minimal inflation due to population stratification (Fig. 2). Fine-mapping of the three regions identified 14 credible SNPs (see ‘Methods’ section). None of the SNPs overlapped with fetal brain eQTL. However, two of these SNPs mapped onto two genes—*LSAMP* and *PTMAP8*, both of chromosome 3—using chromatin interaction data in the fetal brain. Of these, *LSAMP* is a neuronal adhesion molecule in the limbic system of the developing brain. In addition, gene-based analysis identified four significant genes *SDCCAG8*, *ZSWIM6*, *ZNF574*, and *FUT8* (Supplementary Data 2). Of these, mutations in *ZSWIM6* cause a neurodevelopmental disorder with, in some cases, co-morbid autism and unusually repetitive movements and behaviour⁵⁶. As supporting analyses, we investigated the direction of effect for all independent SNPs with $P < 1 \times 10^{-6}$ in the non-stratified SQ-R GWAS in GWAS of autism⁵⁷, educational attainment⁵⁸, and cognitive aptitude⁵⁹. Five out of six SNPs tested had concordant effect direction in the GWAS for educational attainment and GWAS for cognitive aptitude ($P = 0.21$, two-sided binomial sign test for each comparison). Similarly, four out of five SNPs tested had concordant effect direction in the GWAS for autism (Supplementary Data 3a) ($P = 0.37$, two-sided binomial sign test). For these three phenotypes, we additionally assessed effect direction concordance using binomial sign test at less stringent *P*-value thresholds in the SQ-R GWAS, after LD-based clumping ($P < 1$, 0.5, 0.1 and $1 \times 10^{-4}$). Binomial sign test was statistically significant at three of the four *P*-value thresholds ($P = 1$, 0.5 and

0.1) for all three phenotypes but not statistically significant at $P = 1E-4$, presumably due to the low statistical power (Supplementary Data 3b). In addition, we tested effect direction concordance ($P < 1 \times 10^{-6}$) in a GWAS (N = 1981) of ‘insistence on sameness’, a phenotype similar to systemising (see ‘Methods’ section). Four out of five SNPs had a concordant effect direction including the two SNPs with $P < 5 \times 10^{-8}$ in the non-stratified SQ-R GWAS (P = 0.37, two-sided binomial sign test).

**Genetic correlation between the SQ-R and other phenotypes.** Additive SNP-based heritability ($h_{SNP}^2$) calculated using LDSR was $0.12 \pm 0.012$ for the SQ-R ($P = 1.2 \times 10^{-20}$). Despite small but significant sex differences in the SQ-R scores, there was no significant difference in $h_{SNP}^2$ between males and females (P = 0.34) (Supplementary Fig. 3 and Supplementary Data 4), which was strengthened by the high genetic correlation between males and females ($1 \pm 0.17$; P = $3.91 \times 10^{-10}$), suggesting a similar polygenic architecture between sexes. The per-SNP effect for the most significant SNPs was small, suggesting a highly polygenic architecture ($R^2 = 0.001–0.0002$%, after correcting for winner’s curse, Supplementary Data 5).

Partitioned heritability for functional categories identified significant enrichment for evolutionary conserved regions, transcription start sites, fetal DNase hyper-sensitivity sites, and H3 lysine 27 acetylation (H3K27ac), suggesting a prominent role for regulatory and conserved genomic regions in systemising (Supplementary Data 6). Partitioning heritability based on tissue-specific active chromatin marks identified a significant enrichment for brain specific chromatin signatures. Notably, this enrichment was significant in both adult and fetal brain specific active chromatin marks (Supplementary Data 7 and Supplementary Fig. 4). Enrichment for genes expressed in the brain was high
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 2** Manhattan and QQ-plots for the three GWAS. Manhattan plot for the three SQ-R GWAS: non-stratified (a), males-only (c), females-only (e). Significant SNPs are highlighted in red. QQ-plots for the three SQ-R GWAS: non-stratified (b), males-only (d), females-only (f). SQ-R non-stratified (N = 51,564): λGC = 1.10, LSDR intercept = 0.99, SQ-R males-only (N = 26,063): λGC = 1.06, LSDR intercept = 0.99, SQ-R females-only (N = 25,501): λGC = 1.05, LSDR intercept = 1.01</figcaption></figure>

but failed to reach statistical significance after correcting for the multiple tests conducted (Supplementary Fig. 5 and Supplementary Data 8).

We identified a significant positive genetic correlation between the SQ-R and autism as well as measures of intelligence (cognitive aptitude and educational attainment) (Supplementary Data 9 and Fig. 3a). Of all the psychiatric conditions tested (see ‘Methods’ section), SQ-R was only significantly genetically correlated with autism ($r_g = 0.26 \pm 0.06$; $P = 3.35 \times 10^{-5}$), demonstrating the relative specificity of the correlation of the SQ-R to autism. Notably, the absolute magnitude of the genetic correlation between autism and the SQ-R is similar to the genetic correlation between autism and self-reported empathy (measured using the Empathy Quotient (EQ$^{60}$): $r_g = -0.27 \pm 0.07$) and scores on the Social and Communication Disorders Checklist (SCDC$^{61}$): $r_g = 0.27 \pm 0.13$). Controlling for the genetic effects of educational attainment on the SQ-R GWAS using genome-wide inferred statistics (GWIS) (see ‘Methods’ section) attenuated the genetic correlation with autism only modestly, suggesting that the SQ-R scores are genetically correlated with autism independently of the genetic effects of education (Fig. 3b and Supplementary Data 10). We validated this using genomic structural equation modelling

(GSEM) (see ‘Methods’ section) using both educational attainment and cognitive aptitude (Fig. 3c). Further, the SQ-R was not genetically correlated with any of the social measures related to autism—friendship and family relationship satisfaction, scores on a self-report measure of empathy (the EQ), and the scores on the Social and Communication Disorders Checklist (SCDC), which is a measure of social and communication difficulties (see Supplementary Note for how these traits map onto social domains in autism). Estimates of genetic correlations between SQ-R scores and the various social traits are also small, suggesting that there is limited shared genetics between social autism traits and the SQ-R.

**Genetic correlations between social/non-social traits and psychiatric conditions.** To understand the genetic relationship between the SQ-R and autism in a broader context, we evaluated the genetic correlations between multiple phenotypes with evidence of significant genetic correlation with autism (15 phenotypes in total, see ‘Methods’ section for a list of phenotypes included). Clustering highlighted three broad clusters: a social cluster, a psychiatric cluster, and an intelligence cluster (Fig. 4a and Supplementary Tables 11 and 12). The SQ-R clusters closely with measures of intelligence, but while educational attainment
<figure><img src="image_11.png" /><figcaption>**Fig. 3** Genetic correlation between the SQ-R and other phenotypes, and GWIS and GSEM estimates between SQ, educational attainment and cognitive aptitude. **a** Genetic correlations between the SQ-R and multiple other phenotypes provided. The bars represent 95% confidence intervals. Sample sizes and PMID are provided in Supplementary Data 9. The following genetic correlations were significant after Bonferroni correction: autism ($r_g$ = 0.26 ± 0.06; $P$ = 3.35 × 10⁻⁵, $N$ = 46,350), years of schooling ($r_g$ = 0.13 ± 0.03; $P$ = 4.73 × 10⁻⁵, $N$ = 293,723), college completion ($r_g$ = 0.18 ± 0.05; $P$ = 1.30 × 10⁻³, $N$ = 95,427), and cognitive aptitude ($r_g$ = 0.19 ± 0.04; $P$ = 2.35 × 10⁻⁵, $N$ = 78,308). **b** Results of the GWIS analysis. Red lines represent genetic correlation with the SQ-R, blue lines represent genetic correlations with the SQ-R independent of the genetic effects of educational attainment. The bars represent 95% confidence intervals. **c** Path diagrams providing the results of the standardised SEM models to investigate whether the SQ-R is genetically correlated with autism independent of the genetic effects of cognitive aptitude ($CA_g$) and educational attainment ($EA_g$). GWIS genome-wide inferred statistics, GSEM genomic structural equation modelling</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Fig. 4** Genetic correlation of autism and related traits, and genetic correlations between social and non-social traits and multiple psychiatric conditions. **a** Correlogram of genetic correlations between all phenotypes that are genetically correlated with autism. Please note the upper and lower triangle are identical. Asterisk (provided only in the lower triangle) represents significant correlations after Bonferroni correction. Genetic correlations have been clustered using hierarchical clustering. Colour provides the magnitude of genetic correlation. **b** Genetic correlation between empathy, friendship satisfaction, and systemising with nine psychiatric conditions. Only autism was significantly genetically correlated with all three phenotypes. Full results are present in Supplementary Data 11</figcaption></figure>
and cognitive aptitude are significantly genetically correlated with multiple social traits and psychiatric conditions, the SQ-R is only genetically correlated with autism.

Given that the two major domains of autism as identified by the DSM-5 are persistent difficulties in social interaction and communication and unusually restrictive, stereotyped, and repetitive interests¹, we hypothesised that the combination of significant negative genetic correlation with social traits (friendship satisfaction and empathy) and significant positive genetic correlation with SQ-R would be uniquely associated with autism (see 'Methods' section). Indeed, across the nine psychiatric conditions for which we had summary GWAS statistics, this combination was uniquely observed for autism (Fig. 4b, Supplementary Data 13).

**Validation in additional cohorts.** Given that our current analysis focussed on the general population, we sought to investigate whether polygenic scores from the SQ-R were associated with social and non-social autism domains in 2221 autistic individuals from the Simons simplex collection (see 'Methods' section). We hypothesised that SQ-R may be significantly associated with the non-social domain in autism, but not associated with the social domain in autism. Polygenic scores for SQ-R were significantly associated with scores on the Repetitive Behaviour Scale-Revised (RBS-R) (beta = 0.052 ± 0.02, *P* = 0.013), but not on the social and communication subscale of ADOS-G (beta = -0.00099 ± 0.018, *P* = 0.95) after adjusting for multiple test (Bonferroni alpha = 0.025). We validated this in 426 additional individuals of which 401 had a diagnosis of autism with RBS-R scores from the EU-AIMS LEAP, AGRE, and Paris cohorts. Here, we identified a concordant effect direction for polygenic scores of the SQ-R (beta = 0.02 ± 0.05, *P* = 0.65), although the results were not significant potentially due to the small sample size. Inverse-variance meta-analysis of the discovery and the validation cohorts marginally improved the significance of the association (beta = 0.047 ± 0.018, *P* = 0.010), and the results remained statistically significant (Bonferroni alpha = 0.025). In a separate sample of 475 autistic individuals from the AGRE cohort, polygenic scores for the SQ-R were not associated with the social and communication subscale of ADOS-G (beta = -0.046 ± 0.04, *P* = 0.24). Meta-analysis of the two cohorts did not produce a statistically significant result (beta = -0.008 ± 0.016, *P* = 0.60) (see Power calculations in the Supplementary Note). We note that the lack of association between the polygenic scores for the SQ-R and the ADOS-G social and communication subscale is not indicative of absence of shared genetics, but rather indicative of lower shared genetics between the SQ-R and the ADOS-G social and communication subscale than that between the RBS-R and the SQ-R.

Finally, to further validate the results in autistic individuals, we conducted bivariate genetic correlations on scores on the RBS-R and the ADOS-G social and communication subscale in 2989 individuals from the SSC, AGRE, EU-AIMS LEAP and Paris cohorts (2964 autistic individuals). Both the RBS-R ($h_{SNP}^2$ = 0.11 ± 0.11, *P* = 0.15) and the ADOS-G social and communication subscale ($h_{SNP}^2$ = 0.26 ± 0.10, *P* = 0.004) had modest $h_{SNP}^2$, though only the latter was statistically significant. We identified a small genetic correlation ($r_g$ = 0.15 ± 0.46, *P* = 0.74), which was not statistically different from 0. Given the small sample size, the genetic correlation is unlikely to be statistically significant. However, the effect was small and statistically less than 1 (*P* = 0.034, one-tailed *t*-test).

## Discussion

Here we present, to our knowledge, the largest GWAS of a non-social trait related to autism in the general population—systemising, measured using the SQ-R. We demonstrate that systemising

is heritable and genetically correlated with autism. Associated loci are enriched in genomic regions containing brain chromatin signatures and we identify three genome-wide significant loci, but these must be replicated in an independent cohort. Despite the modest sample size, our GWAS is well-powered to investigate genetic correlations between various phenotypes including social traits related to autism, as the Z-score of the $h_{SNP}^2$ is above the recommended threshold of four⁶². We identify high sign concordance of the top SNPs in genetically correlated traits, enrichment for active chromatin marks in fetal and adult brain, and significant polygenic score association with the RBS-R. Polygenic score analysis suggests that the shared genetics between systemising and the non-social domain of autism is considerably higher than the shared genetics between systemising and the social domain of autism. In addition, using a smaller sample of autistic individuals, we provide preliminary evidence that the social and non-social domains in autistic individuals have low shared genetics. Our results highlight the need to collect deeper clinical and cognitive information in autistic individuals to better understand the phenotypic heterogeneity in autism.

Most studies model autism and autistic traits as a single phenotype. This has likely arisen because of the difficulties in recruiting and phenotyping sufficient numbers of autistic people. Our study suggests that both in the general population and in autistic individuals, social and non-social autistic traits and symptom domains are genetically dissociable. This may to some extent explain why, compared to GWAS of other psychiatric conditions of roughly similar sample sizes⁵⁷,⁶³⁻⁶⁵, the largest GWAS of autism to date has identified fewer loci. One possible explanation is statistical signal-attenuation because of the underlying heterogeneity. However, this does not necessarily suggest that systemising, or the other individual trait domains are less complex. For instance, we observe similar $h_{SNP}^2$ for SQ-R, self-reported empathy⁴⁶, and the largest and most recent GWAS of autism⁵⁷.

It is important to investigate whether these domains are dissociable in a larger cohort of autistic individuals and identify potential convergence of the two domains in gene expression networks in the developing brain. Our results confirm the need to rethink our understanding of autism as existing along a single dimension³,⁶⁶. We hypothesise that the dissociation of the two domains will extend to other research modalities in studies of autism and autistic traits. It is important to note that, while our results demonstrate two broadly dissociable autistic trait domains in the general population and in autistic individuals, more research is needed to identify other potentially dissociable domains and to investigate whether this dissociability is driven by different designs of phenotypic instruments (e.g. self-report vs informant report). For example, our research does not make a distinction between communication and social interaction abilities, or between sensory difficulties and repetitive behaviours, and future molecular genetic studies may identify varying levels of overlap between these domains. The same principle applies to other research modalities (neuroimaging, cognitive studies, hormonal assays, etc.) investigating the biology of autism and autistic traits. These different symptom domains of autism may contribute to different co-morbidities. Our results identify shared genetics between the social traits related to autism and psychiatric conditions such as schizophrenia and depression, but limited shared genetics between the SQ-R and these conditions.

## Methods

**Participants.** The current study included participants from 23andMe (primary GWAS - SQ-R), from ALSPAC (GWAS of scores on the Social and Communication Disorders Checklist (SCDC)) and autistic individuals from the Simons
Simplex Collection (SSC), the Autism Genetic Resource Exchange (AGRE), and the EU-AIMS LEAP and PARIS cohorts.

**23andMe participants.** Research participants in the GWAS of the SQ-R were from 23andMe and are described in detail elsewhere<sup>67,68</sup>. All participants provided informed consent and answered surveys online according to a human subjects' research protocol, which was reviewed and approved by Ethical & Independent Review Services, an external AAHRPP-accredited private institutional review board (http://www.eandireview.com). All participants completed the online version of the SQ-R on the 23andMe participant portal. Only participants who were primarily of European ancestry (97% European Ancestry) were selected for the analysis using existing methods<sup>69</sup>. Unrelated individuals were selected using a segmental identity-by-descent algorithm<sup>70</sup>. A total of 51,564 participants completed the SQ-R (males = 26,063, and females = 25,501).

**ALSPAC participants.** ALSPAC is a longitudinal cohort which recruited pregnant mothers in the Avon region of the UK. The ALSPAC cohort comprises 14,541 initial pregnancies from women in Avon resulting in a total of 13,988 children who were alive at 1 year of age. Children were enrolled in additional phases, described in greater detail elsewhere<sup>71</sup>. This study received ethical approval from the ALSPAC Law-and-Ethics Committee, and the Cambridge Human Biology Research Ethics Committee. Written informed consent was obtained from parent or a responsible legal guardian for the child to participate. Assent was obtained from the child participants where possible. We conducted a GWAS of scores on the SCDC in 5,421 individuals from ALSPAC.

**Other cohorts.** We included data from four cohorts to conduct polygenic score and bivariate genetic correlation analyses. The SSC (n = 2221 unrelated autistic individuals) consists of simplex autistic families, and are described elsewhere<sup>72</sup>. The AGRE cohort (n = 482 unrelated autistic individuals) consists of multiplex autism families, details of which are provided elsewhere<sup>73</sup>. In addition, we included 401 individuals (including 25 neurotypical individuals) from the EU-AIMS LEAP<sup>74</sup> and Paris<sup>75</sup> cohorts. Across all cohorts, we included only unrelated individuals, who were predominantly of European Ancestry as defined by genetic principal components (5 SD deviations above or below the mean of PC1 and PC2 from the HapMap CEU population).

Additionally, we also included data from 1981 unrelated individuals (1000 males, 981 females) from the Nijmegen Biomedical Study (NBS) to provide support for the independent SNPs with $P < 1 \times 10^{-6}$ in the non-stratified GWAS. Participants were asked the question: 'It upsets me if my daily routine is disturbed', which is related to a non-social domain of autism, and is similar to an item in the Autism Spectrum Quotient. Further information including genotyping and quality control is provided elsewhere<sup>43</sup>. Genetic association for the top SNPs were conducted using age, sex, and the first five genetic principal components as covariates using linear regression.

**Phenotypes.** The primary phenotype for this study is the SQ-R, which was used to conduct a GWAS in participants from 23andMe. The SQ-R is self-report measure of systemising drive, or interest in rule-based patterns<sup>10</sup>. The SQ-R taps a variety of domains of systemising, such as interest in mechanical (e.g. car engines), abstract (e.g. mathematics), natural (e.g. the weather), motor (e.g. knitting), and collectible (e.g. stamp collecting) systems. There are 75 items on the SQ-R, with a maximum score of 150 and a minimum score of 0. Scores on the test are normally distributed<sup>10</sup>. The SQ-R has good cross-cultural stability and good psychometric properties with Cronbach's alpha ranging from 0.79 to 0.94 in different studies<sup>76</sup>. Test-retest reliability available in a Dutch sample indicated high reliability of 0.79 (Pearson correlation)<sup>76</sup>. This was supported by another study in 4058 individuals which identified high internal cohesion<sup>77</sup>. Exploratory followed by confirmatory factor analysis using Rasch modelling suggests that the SQ-R is unidimensional<sup>77</sup>. A sex difference has been observed in multiple studies with males, on average, scoring significantly higher than females<sup>10,51</sup>. Criterion validity shows that the SQ-R has a modest but significant correlation with the mental rotation test ($r = 0.25, P = 0.013$), as well as its subscales<sup>78</sup>. Autistic individuals, on average, score higher on the SQ-R in multiple different studies<sup>10,51,79</sup>. Further, the SQ-R also predicts autistic traits, with a combination of the SQ-R and the Empathy Quotient predicting as much as 75% of the variance on the autism spectrum quotient, a measure of autistic traits<sup>10</sup>. The SQ-R has been validated using a short form in a very large population of 600,000 controls and 36,000 autistic individuals<sup>12</sup>.

In addition, we used the following secondary phenotypes: SCDC in ALSPAC, ADOS-G social and communication scores and the RBS-R in the other cohorts. We also used a single question which is a measure of 'insistence on sameness' in the NBS cohort.

The SCDC is a questionnaire that measures difficulties in verbal and nonverbal communication, and social interaction including reciprocal social interaction<sup>61</sup>. The questionnaire consists of 12 questions, with scores ranging from 0 to 24, and with higher scores reflecting difficulties in social interaction and communication. The SCDC has good internal consistency (0.93) and good test-retest reliability (0.81)<sup>61</sup>. The SCDC has reasonable specificity and sensitivity in distinguishing autistic from control individuals<sup>80</sup>. Previous research has demonstrated that the

SCDC is genetically correlated with autism<sup>44,45,57</sup>. We conducted a GWAS of SCDC to investigate whether it is genetically correlated with SQ-R in this study. We used mother-reported SCDC scores on children aged 8. Although SCDC has been measured at different ages in the ALSPAC cohort, we chose SCDC scores measured at age 8 as this has the largest sample size and has high $h_{SNP}^2$<sup>19</sup> ($h^2 = 0.24 \pm 0.07$).

We chose two measures of social and non-social traits in autistic individuals. For the social trait, we used the social and communication domain scores from the ADOS-G, a widely used instrument for diagnosing and assessing autism in four cohorts (SSC, AGRE, EU-AIMS LEAP, and Paris). Participants completed one of the following ADOS-G modules<sup>81</sup>: 1 (used for children with little or no phrase speech), 2 (for children with non-fluent speech), 3 (verbally fluent children), and 4 (verbally fluent adolescents and adults). For this study, we used the raw totals of the scores from the social domain and the communication domain, combined. Scores for all four modules range from 0 to 24. The ADOS-G has high overall internal consistency, and high test-retest reliability for the social and communication subscales<sup>81</sup>. The choice for combining the social and communication domain scores were informed by factor analysis which suggested that the two domains contribute to one underlying factor<sup>82</sup>.

In contrast to the social and communication domain, the restricted and repetitive behaviour domain of the ADOS-G has poor test-retest reliability ($r < 0.6$) and a smaller range of scores (0–8) as it captures fewer repetitive and restrictive behaviour<sup>81</sup>. Hence, for this study, we used scores on the RBS-R<sup>83</sup>. The RBS-R is a measure developed to specifically measure restricted and repetitive behaviours in autistic individuals and captures stereotyped, self-injurious, sameness, compulsive, ritualistic, and restricted behaviour<sup>84</sup>, and has high inter-rater reliability and internal consistency<sup>84</sup>. The RBS-R comprises 43 questions with scores ranging from 0 to 3 for each item based on a Likert scale.

'Insistence on sameness' in the NBS cohort was measured using a single item: 'It upsets me if my daily routine is disturbed'. This is related to a non-social domain of autism, and is again similar to an item in the Autism Spectrum Quotient. Participants were asked to indicate on a 4-point Likert scale 'definitely agree', 'slightly agree', 'slightly disagree', 'definitely disagree'.

**Genotyping, imputation, and quality control and genetic association in the 23andMe cohort.** Details of genotyping, imputation and quality control in the 23andMe cohort are provided elsewhere<sup>47</sup>. Briefly, unrelated participants were included if they had a call rate of >98.5%, and were of primarily European ancestry (97% European ancestry). A total of 1,030,430 SNPs (including InDels) were genotyped or imputed. SNPs were excluded if: they failed the Hardy-Weinberg equilibrium test at $P < 10^{-20}$; had a genotype rate of <90%; they failed the parent-offspring transmission test using trio data in the larger 23andMe research participant database; or if allele frequencies were significantly different from the European 1000 Genomes reference data ($\chi^2$ test, $P < 10^{-20}$). Phasing was conducted using Beagle (version 3.3.1)<sup>85</sup> in batches of 8000–9000 individuals. This was followed by imputation against all-ethnicity 1000 Genomes haplotypes (excluding monomorphic and singleton sites) using Minimac<sup>86</sup>. Genetic association analyses were restricted to SNPs with a minor-allele frequency > 1%. After quality control, 9,955,952 SNPs (imputed and genotyped) were included in the GWAS.

Our primary analysis was an additive model of genetic effects and was conducted using a linear regression with age, sex, and the first five ancestry principal components included as covariates. In addition, given the modest sex difference, we also conducted sex-stratified analyses. SNPs were considered significant at a genome-wide threshold of $P < 5 \times 10^{-8}$. Leading SNPs were identified after LD-pruning using Plink ($r^2 > 0.8$). Winner's curse correction was conducted using an FDR-based shrinking<sup>87</sup>.

We calculated variance explained by first standardising the regression estimates and then squaring the estimates. This is equivalent to:

$$R^2 = B_j^2 \frac{2(MAF_j)(1 - MAF_j)}{\sigma_j^2},$$

where $R^2$ is the proportion of variance explained for SNP<sub>j</sub>, $B_j^2$ is the non-standardised regression coefficient, MAF is the minor-allele frequency for SNP<sub>j</sub>, and $\sigma_j^2$ is the variance of SQ. Further details of this formula are provided in the Supplementary Note.

**Genotyping, imputation, and quality control and genetic association in the ALSPAC.** The SCDC<sup>61</sup> scores were calculated from children of the 90 s (ALSPAC cohort)<sup>71</sup>, in children aged 8. In total, SCDC scores were available on $N = 7,825$ children. From this, we removed individuals for whom complete SCDC scores were not available. After excluding related individuals and individuals with no genetic data, data was available on a total of $N = 5,421$ unrelated individuals.

Participants were genotyped using the Illumina® HumanHap550 quad chip by Sample Logistics and Genotyping Facilities at the Wellcome Sanger Institute and LabCorp (Laboratory Corporation of America) using support from 23andMe. Individuals were excluded based on gender mismatches, high missingness (>3%), and disproportionate heterozygosity. We restricted subsequent analyses to individuals of European descent (CEU), which were identified by multi-
dimensional scaling analysis and compared with Hapmap II (release 22).

Individuals were also removed if cryptic relatedness, assessed using identity by descent, was >0.1. Genotyped SNPs were filtered out if they had >5% missingness, violated Hardy-Weinberg equilibrium (P < 1 × 10⁻⁶), and had a minor-allele frequency < 1%, resulting in a total of 526,688 genotyped SNPs. Haplotypes were estimated using data from mothers and children using ShapeIT (v2.r644)⁸⁸. Imputation was performed using Impute2 V2.2.2⁸⁹ against the 1000 genomes reference panel (Phase 1, Version 3). Imputed SNPs were excluded from all further analyses if they had a minor-allele frequency < 1% and info < 0.8. After quality control, there were 8,282,911 genotyped and imputed SNPs that were included in subsequent analyses.

Dosage data from BGEN files were converted using hard-calls, with calls with uncertainty > 0.1 treated as missing data. Post-imputation, we excluded SNPs that deviated from Hardy-Weinberg equilibrium (P < 1 × 10⁻⁶), with minor-allele frequency < 0.01 and missing call rates > 2%. We further excluded individuals with genotype missing rates > 5%. The SCDC score was not normally distributed so we log-transformed the scores and ran regression analyses using the first two ancestry principal components and sex as the covariates using Plink 2.0 (ref. ⁹⁰).

The log-transformed SCDC scores (henceforth, SCDC scores) had a modest but significant h²SNP as quantified using LSDR (h²SNP = 0.12 ± 0.05). LSDR intercept (0.99) suggested that there was no inflation in GWAS estimates due to population stratification. The λGC was 1.013. We replicated the previously identified genetic correlation with autism⁵⁷ (constrained intercept) using our SCDC GWAS (r_g = 0.45 ± 0.18, P = 0.01). In addition, we also identified a negative genetic correlation between educational attainment⁵⁸ and SCDC (r_g = -0.30 ± 0.11, P = 0.007).

**Genomic inflation factor, heritability, and functional enrichment for the SQ-R GWAS.** LSDR⁹¹,⁹² was used to calculate for inflation in test statistics due to unaccounted population stratification. Heritability was calculated using LSDR using the north-west European LD scores. Difference in heritability between males and females was quantified using:

$$Z = \frac{h_{\text{males}}^2 - h_{\text{females}}^2}{\sqrt{SE_{\text{males}}^2 + SE_{\text{females}}^2}},$$

where Z is the Z-score for the difference in heritability for a trait, (h²males - h²females) is the difference h²SNP estimate in males and females, and SE is the standard errors for heritability. Two-tailed P-values were calculated and reported as significant if P < 0.05.

For the primary GWAS (non-stratified analyses), we conducted functional annotation using FUMA⁹³. We restricted our analyses to the non-stratified analyses due to the high genetic correlation between the sexes and the low statistical power of the sex-stratified GWAS. We conducted gene-based association analyses using MAGMA⁹⁴ within FUMA and report significant genes after using a stringent Bonferroni corrected P < 0.05. In addition, we conducted enrichment for tissue specific expression and pathway analyses within FUMA. For the significant SNPs, we investigated enrichment for eQTLs using brain tissues in the BRAINEAC and GTEx⁹⁵ database within FUMA. We further conducted partitioned heritability for tissue-specific active chromatin marks and baseline functional categories using extended methods in LSDR⁹⁶.

**Hi-C-based annotations of fine mapped loci.** We fine mapped three genome-wide significant loci (index SNPs: rs4146336 and rs1559586 or SQ; rs8005092 for SQ-R males) to obtain credible SNPs. First, we selected SNPs with P < 0.01 that are located in the LD region (r² > 0.6) with an index SNP. LD structure within a locus was constructed by calculating correlations between SNPs within a locus (1KG v20130502). CAVIAR⁹⁷ was then applied to the summary association statistics and LD structure for each index SNP to generate potentially causal (credible) SNPs with a posterior probability of 0.95. In total, we identified 14 credible SNPs from the three GWS loci.

For each locus, candidate genes were identified by mapping credible SNPs based on physical interactions in foetal brain as previously described⁹⁸. One locus (index SNP rs4146336) was mapped to two genes, LSAMP and PTMAP8, indicating that two credible SNPs (rs13066948 and rs11713893) located in this locus physically interact with these genes.

**Genetic correlation.** For all phenotypes, we performed genetic correlation without constraining the intercept using LSDR. We identified significant genetic correlations using a Bonferroni adjusted P-value < 0.05. For the primary genetic correlation analysis with SQ-R, we included psychiatric conditions⁵⁷,⁶³,⁹⁹⁻¹⁰², personality traits¹⁰³⁻¹⁰⁵, measures of intelligence⁵⁸,⁵⁹,¹⁰⁶,¹⁰⁷, and social traits related to autism⁴⁶,⁵⁵ including scores on the SCDC, as previous research has investigated the phenotypic correlation between these domains and systemising¹⁰⁷,¹⁰⁸,¹¹⁰.

To understand the correlation between systemising and various phenotypes that have been genetically correlated with autism, we used GWAS data from 15 phenotypes including autism. 10 of these phenotypes (cognitive aptitude⁵⁹, educational attainment⁵⁸, tiredness¹¹³, neuroticism¹⁰³, subjective wellbeing¹⁰³, schizophrenia¹¹⁴, major depression¹⁰², depressive symptoms¹⁰³, ADHD⁶³,

and chronotype¹¹⁵), have been previously reported to be significantly genetically correlated with autism out of 234 phenotypes tested using LDHub⁶² (P < 2.1 × 10⁻⁴). We excluded college degree from this list, as previous work has identified near perfect genetic correlation between educational attainment and college degree⁵⁸. In addition, we included data from friendship satisfaction⁵⁵, family relationship satisfaction⁵⁵, systemising, and self-reported empathy⁴⁶, all of which are also significantly genetically correlated with autism with P < 2.1 × 10⁻⁴. These four additional phenotypes were not included in the previous paper which investigated genetic correlations with autism. Details of sample sizes with PMIDs/ DOIs are provided in Supplementary Data 12. Cross trait genetic correlations were computed for all 15 phenotypes, and results were corrected for multiple testing using Bonferroni correction. A correlogram was created after using hierarchical clustering to cluster the phenotypes.

To investigate whether the combination of negative genetic correlation social traits and positive genetic correlation for non-social traits is specific to autism, we conducted a genetic correlation between all psychiatric conditions for which we had access to summary GWAS statistics (ADHD⁶³, Anxiety¹¹⁶, Autism⁵⁷, Anorexia¹⁰¹, Bipolar Disorder⁹⁹, Major Depressive Disorder¹⁰², OCD¹¹⁷,¹¹⁸, PTSD¹¹⁹, and Schizophrenia¹¹⁴) and SQ-R, self-reported empathy measured using the EQ⁴⁶ and friendship satisfaction⁵⁵. We chose friendship satisfaction and self-reported empathy as representative of social traits as these are the most relevant to the social domain of autism for which we had access to GWAS summary statistics. The EQ is a short, 40-item self-report measure of empathy, which has been widely used and has good psychometric properties⁶⁰,¹²⁰. In addition, differences in aspects of empathy compared to the neurotypical population have been widely reported in autism⁵⁰,⁵¹,¹²¹, and empathy is one of the items in measures such as ADOS-G. Friendship satisfaction was selected as difficulties in making friends is listed as a criteria for an autism diagnosis in the DSM-5¹.

**GWIS and GSEM.** To investigate whether the SQ-R is genetically correlated with autism independent of the genetic effects of educational attainment, we constructed a unique SQ-R phenotype after conditioning on the genetic effects of educational attainment using GWIS¹²². GWIS takes into account the genetic covariance between the two phenotypes to calculate the unique component of the phenotypes as a function of the genetic covariance and the h²SNP. Before performing GWIS, we standardised the beta coefficients for the SQ-R GWAS by using the following formula:

$$\hat{B}_{\text{std}} = \hat{B} \sqrt{\frac{2(\text{MAF})(1 - \text{MAF})}{\sigma_y^2}},$$

where Bstd is the standardised regression coefficients, B is the regression coefficient obtained from the non-standardised GWAS, MAF is the minor-allele frequency, σ²y is the variance of the SQ-R. This equation is explained in detail in the Supplementary Note. We conducted GWIS using only educational attainment as we were unclear whether the GWAS of cognitive aptitude⁵⁹ was conducted on a standardised phenotype. Further, there is a high genetic correlation between cognitive aptitude and educational attainment. In addition to GWIS, to validate the findings, we conducted GSEM¹²³, a complementary but independent method. GSEM uses the genetic correlations and covariances calculated using LSDR after accounting for sample overlap.

**Polygenic scores in the SSC, AGRE, EU-AIMS LEAP, and Paris cohorts.** We generated polygenic scores for SQ-R (mean weighted score of all the alleles that contribute to higher systemising) in 2221 probands from the SSC (Discovery dataset). We downloaded genotype data from the SSC from SFARI base (https://www.sfari.org/resource/sfari-base/). Individuals were genotyped on three different platforms: Illumina Omni2.5, Illumina 1Mv3, or Illumina 1Mv1. Informed consent or assent was obtained from all participants. In addition, the research team obtained ethical approval from the Cambridge Human Biology Research Ethics Committee to access and analyse the de-identified data from the SSC. We conducted stringent quality control and imputation separately for each platform. The full pipeline is available here: https://github.com/autism-research-centre/SSC_liftover_imputation. Briefly, individuals were excluded if they had: a genotyping rate < 95%, excessive or low heterozygosity (less or more than 3 SD from the mean), mismatched reported and genetic sex, and families with Mendelian errors > 5%. We further removed SNPs that significantly deviated from Hardy-Weinberg equilibrium (P < 1 × 10⁻⁶), had Mendelian errors in >10% of the families, and SNPs that were not genotyped in >10% of the families. We then conducted multidimensional scaling using the HapMap3 phase 3 population using the unrelated individuals CEU and TSI populations as representatives of the European population. This was conducted only in the parents to retain unrelated individuals for multi-dimensional scaling. Genetic principal components were calculated using only SNPs with minor-allele frequency > 5%, and pruning the SNPs in Plink using an r² of 0.2. We excluded families from further downstream analyses if either one the parents were greater or less than 5 standard deviations from the means of the first two genetic principal components calculated using only the unrelated individuals in HapMap3 CEU and TSI populations. Quality control was done using Plink v 1.9 and R. Phasing and imputation were conducted using the Michigan
Imputation Server (https://imputationserver.sph.umich.edu/start.html) using the 1000 genomes Phase 3 v5 as the reference panel.

Polygenic scores were generated using PRSice2 (https://choishingwan.github.io/PRSice/) for the SQ-R using the non-stratified GWAS data. We calculated the mean polygenic score for each of the 2221 probands in the SSC, after clumping SNPs using an R² threshold of 0.1. Prior to generating polygenic scores, we confirmed that the probands were not related to each other using identity by descent PI-HAT > 0.15 as a relatedness cut-off. We used a P-value threshold of 1 as previous research on educational attainment, subjective wellbeing and social relationship satisfaction, all suggest that the maximum variance explained is at a threshold of 1 (refs. ⁵⁸,¹⁰³). This is expected for highly polygenic traits where many SNPs incrementally contribute to the variance explained¹²⁴. Polygenic scoring was done using standardised scores on two different phenotypes as the dependent variable (RBS-R and the social and communication domain of the ADOS-G). We included sex, platform, the first 15 genetic principal components and standardised full-scale IQ as covariates. In addition, for the analysis of ADOS-G, we included the ADOS-G module as a covariate. Linear regression was conducted in R. A total of 135,233 SNPs were included in the polygenic score analyses after clumping and thresholding.

To validate the polygenic scores, we conducted additional polygenic score analysis using data combined from the AGRE, EU-AIMS LEAP and Paris cohorts. We followed similar quality control and imputation procedures to the SSC cohort. Given that this dataset was a mix of related and unrelated individuals, we chose unrelated individuals using a genomic relationship matrix (GRM) as provided in GCTA (grm-cutoff 0.05)¹²⁵. To calculate GRMs, we included only SNPs with minor-allele frequency > 1%. Scripts are provided here: https://github.com/vwarrier/PARIS_LEAP_analysis. Polygenic scores were calculated using PRSice2 as described for the SSC data. Given the differences in dataset, polygenic scores were calculated separately for the AGRE dataset, and the EU-AIMS LEAP and Paris datasets combined. For each regression, we included sex and the first ten genetic principal components (standardised). The dependent variables were standardised scores on the RBS-R (N = 426) and the ADOS-G social and communication subscale (N = 475). IQ information was unavailable for most individuals, and hence we did not include IQ as a covariate. We combined the results of the EU-AIMS LEAP and Paris cohorts, and the AGRE dataset using inverse-variance weighted fixed-effect meta-analysis using the formula below:

$$ 
\begin{aligned}
w_i &= 1/\text{SE}_{i}^2 \\
\text{SE}_{\text{meta}} &= \sqrt{1/\text{SE}_{i} w_i} , \\
\text{Beta}_{\text{meta}} &= \Sigma_{i} \beta_{i} w_{i} / \Sigma_{i} w_{i}
\end{aligned}
 $$

where βᵢ is the standardised regression coefficient of the polygenic scores, SEᵢ is the associated standard error, and wᵢ is the weight.

**Bivariate GREML.** We conducted bivariate genetic correlation using GCTA GREML to test the genetic correlation between the ADOS social and communication domains and the RBS-R scores. We created a GRM after including autistic individuals from the SSC, AGRE, EU-AIMS LEAP, and Paris cohorts. We excluded SNPs and individuals using the same quality control pipeline as applied to the SSC dataset outlined in the section above. We further restricted our analysis only to SNPs with a minor-allele frequency > 1%. We excluded related individuals (–grm-cutoff 0.05) resulting in a total of 2989 individuals. Of this, 2652 individuals had scores for the ADOS social and communication domain and 2550 individuals had scores on the RBS-R. We included sex and the first ten genetic principal components as covariates.

**Reporting summary.** Further information on research design is available in the Nature Research Reporting Summary linked to this article.

**Data availability**

The SQ-R GWAS results are available from 23andMe. The full set of summary statistics can be made available to qualified investigators who enter into an agreement with 23andMe that protects participant confidentiality. Interested investigators should email dataset-request@23andme.com for more information. Top SNPs (n = 10,000) can be visualised here: https://ghfc.pasteur.fr.Data. for ALSPAC can be requested here: http://www.bristol.ac.uk/alspac/researchers/access/. Data from the Simons Simplex Collection can be requested here: https://www.sfari.org/resource/sfari-base/. Summary GWAS statistics were downloaded from the PGC consortium: http://www.med.unc.edu/pgc/results-and-downloads. Data for chronotype was downloaded from http://www.t2diabetesgenes.org/data/. Data for self-reported tiredness was downloaded from http://www.ccace.ac.uk/node/335. Additional source data are available in Supplementary Data 1–13.

**Code availability**

Genomic-SEM: https://github.com/MichelNivard/GenomicSEM, GWIS: https://sites.google.com/site/mgnivard/gwis, Plink: https://www.cog-genomics.org/plink2/, PRSice2: https://choishingwan.github.io/PRSice/, CAVIAR: http://genetics.cs.ucla.edu/caviar/, Michigan Imputation Server: https://imputationserver.sph.umich.edu/index.html, custom code for quality control of the SSC and the other cohorts can be downloaded from https://github.com/autism-research-centre/SSC_liftover_imputation (https://doi.org/10.5281/zenodo.3342561) and from https://github.com/vwarrier/ PARIS_LEAP_analysis (https://doi.org/10.5281/zenodo.3342569).

Received: 20 December 2018 Accepted: 19 July 2019
Published online: 03 September 2019

**References**

1. American Psychiatric Association. *Diagnostic and Statistical Manual of Mental Disorders* 5th edn (American Psychiatric Publishing, Arlington, VA, 2013).

2. Happé, F. & Ronald, A. The 'Fractionable Autism Triad': a review of evidence from behavioural, genetic, cognitive and neural research. *Neuropsychol. Rev.* **18**, 287–304 (2008).

3. Happé, F., Ronald, A. & Plomin, R. Time to give up on a single explanation for autism. *Nat. Neurosci.* **9**, 1218–1220 (2006).

4. Shuster, J., Perry, A., Bebko, J. & Toplak, M. E. Review of factor analytic studies examining symptoms of autism spectrum disorders. *J. Autism Dev. Disord.* **44**, 90–110 (2014).

5. Mandy, W. P. L. & Skuse, D. H. Research review: what is the association between the social-communication element of autism and repetitive interests, behaviours and activities? *J. Child Psychol. Psychiatry* **49**, 795–808 (2008).

6. Palmer, C. J., Paton, B., Enticott, P. G. & Hohwy, J. 'Subtypes' in the presentation of autistic traits in the general adult population. *J. Autism Dev. Disord.* **45**, 1291–1301 (2015).

7. Guthrie, W., Swineford, L. B., Wetherby, A. M. & Lord, C. Comparison of DSM-IV and DSM-5 factor structure models for toddlers with autism spectrum disorder. *J. Am. Acad. Child Adolesc. Psychiatry* **52**, 797–805.e2 (2013).

8. Frazier, T. W., Youngstrom, E. A., Kubu, C. S., Sinclair, L. & Rezai, A. Exploratory and confirmatory factor analysis of the autism diagnostic interview-revised. *J. Autism Dev. Disord.* **38**, 474–480 (2008).

9. Grove, R., Baillie, A., Allison, C., Baron-Cohen, S. & Hoekstra, R. A. Empathizing, systemizing, and autistic traits: Latent structure in individuals with autism, their parents, and general population controls. *J. Abnorm. Psychol.* **122**, 600–609 (2013).

10. Wheelwright, S. J. et al. Predicting autism spectrum quotient (AQ) from the systemizing quotient-revised (SQ-R) and empathy quotient (EQ). *Brain Res.* **1079**, 47–56 (2006).

11. Svedholm-Häkkinen, A. M., Halme, S. & Lindeman, M. Empathizing and systemizing are differentially related to dimensions of autistic traits in the general population. *Int. J. Clin. Heal. Psychol.* **18**, 35–42 (2018).

12. Greenberg, D. M., Warrier, V., Allison, C. & Baron-Cohen, S. Testing the empathizing-systemizing theory of sex differences and the extreme male brain theory of autism in half a million people. *Proc. Natl. Acad. Sci. USA* **115**, 12152–12157 (2018).

13. Ronald, A., Happe, F. & Plomin, R. The genetic relationship between individual differences in social and nonsocial behaviours characteristic of autism. *Dev. Sci.* **8**, 444–458 (2005).

14. Ronald, A. et al. Genetic heterogeneity between the three components of the autism spectrum: a twin study. *J. Am. Acad. Child Adolesc. Psychiatry* **45**, 691–699 (2006).

15. Dworzynski, K., Happé, F., Bolton, P. & Ronald, A. Relationship between symptom domains in autism spectrum disorders: a population based twin study. *J. Autism Dev. Disord.* **39**, 1197–1210 (2009).

16. Norbury, C. F. Practitioner review: social (pragmatic) communication disorder conceptualization, evidence and clinical implications. *J. Child Psychol. Psychiatry* **55**, 204–216 (2014).

17. Uljarević, M., Evans, D. W., Alvares, G. A. & Whitehouse, A. J. O. Short report: relationship between restricted and repetitive behaviours in children with autism spectrum disorder and their parents. *Mol. Autism* **7**, 29 (2016).

18. Georgiades, S. et al. Investigating phenotypic heterogeneity in children with autism spectrum disorder: a factor mixture modeling approach. *J. Child Psychol. Psychiatry* **54**, 206–215 (2013).

19. 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).

20. Cantor, R. M. et al. ASD restricted and repetitive behaviours associated at 17q21.33: genes prioritized by expression in fetal brains. *Mol. Psychiatry* **23**, 993–1000 (2017).

21. Alarcón, M., Cantor, R. M., Liu, J., Gilliam, T. C. & Geschwind, D. H. Evidence for a language quantitative trait locus on chromosome 7q in multiplex autism families. *Am. J. Hum. Genet.* **70**, 60–71 (2002).

22. Lowe, J. K., Werling, D. M., Constantino, J. N., Cantor, R. M. & Geschwind, D. H. Social responsiveness, an autism endophenotype: genomewide significant
linkage to two regions on chromosome 8. *Am. J. Psychiatry* **172**, 266–275 (2015).

23. Shao, Y. et al. Fine mapping of autistic disorder to chromosome 15q11-q13 by use of phenotypic subtypes. *Am. J. Hum. Genet.* **72**, 539–548 (2003).

24. Cannon, D. S. et al. Genome-wide linkage analyses of two repetitive behavior phenotypes in Utah pedigrees with autism spectrum disorders. *Mol. Autism* **1**, 3 (2010).

25. Tao, Y. et al. Evidence for contribution of common genetic variants within chromosome 8p21.2-8p21.1 to restricted and repetitive behaviors in autism spectrum disorders. *BMC Genome* **17**, 163 (2016).

26. Ruzich, E. et al. Measuring autistic traits in the general population: a systematic review of the autism-spectrum quotient (AQ) in a nonclinical population sample of 6,900 typical adult males and females. *Mol. Autism* **6**, 2 (2015).

27. Baron-Cohen, S., Wheelwright, S. J., Skinner, R., Martin, J. & Clubley, E. The autism-spectrum quotient (AQ): evidence from Asperger syndrome/high-functioning autism, males and females, scientists and mathematicians. *J. Autism Dev. Disord.* **31**, 5–17 (2001).

28. Posserud, M.-B., Lundervold, A. J. & Gillberg, C. Autistic features in a total population of 7-9-year-old children assessed by the ASSQ (autism spectrum screening questionnaire). *J. Child Psychol. Psychiatry* **47**, 167–175 (2006).

29. Constantino, J. N. & Todd, R. D. Autistic traits in the general population. *Arch. Gen. Psychiatry* **60**, 524 (2003).

30. Constantino, J. N. et al. The factor structure of autistic traits. *J. Child Psychol. Psychiatry* **45**, 719–726 (2004).

31. Hoekstra, R. A. et al. Heritability of autistic traits in the general population. *Arch. Pediatr. Adolesc. Med.* **161**, 372–377 (2007).

32. Ronald, A. & Hoekstra, R. A. Autism spectrum disorders and autistic traits: a decade of new twin studies. *Am. J. Med. Genet. Part B Neuropsychiatr. Genet.* **156**, 255–274 (2011).

33. de Zeeuw, E. L., van Beijsterveldt, C. E. M., Hoekstra, R. A., Bartels, M. & Boomsma, D. I. The etiology of autistic traits in preschoolers: a population-based twin study. *J. Child Psychol. Psychiatry* **58**, 893–901 (2017).

34. Wheelwright, S. J., Auyeung, B., Allison, C. & Baron-Cohen, S. Defining the broader, medium and narrow autism phenotype among parents using the autism spectrum quotient (AQ). *Mol. Autism* **1**, 10 (2010).

35. Sucksmith, E., Roth, I. & Hoekstra, R. A. Autistic traits below the clinical threshold: re-examining the broader autism phenotype in the 21st century. *Neuropsychol. Rev.* **21**, 360–389 (2011).

36. Sasson, N. J., Lam, K. S., Parlier, M., Daniels, J. L. & Piven, J. Autism and the broad autism phenotype: familial patterns and intergenerational transmission. *J. Neurodev. Disord.* **5**, 11 (2013).

37. Constantino, J. N. & Todd, R. D. Intergenerational transmission of subthreshold autistic traits in the general population. *Biol. Psychiatry* **57**, 655–660 (2005).

38. Murray, A. L., McKenzie, K., Kuenssberg, R. & Booth, T. Do the autism spectrum quotient (AQ) and autism spectrum quotient short form (AQ-S) primarily reflect general ASD traits or specific ASD traits? A bi-factor analysis. *Assessment* **24**, 444–457 (2017).

39. Lundström, S. et al. Autism spectrum disorders and autistic-like traits. *Arch. Gen. Psychiatry* **69**, 46 (2012).

40. Robinson, E. B. et al. Evidence that autistic traits show the same etiology in the general population and at the quantitative extremes (5, 2.5, and 1%). *Arch. Gen. Psychiatry* **68**, 1113–1121 (2011).

41. Ronald, A., Happé, F., Price, T. S., Baron-Cohen, S. & Plomin, R. Phenotypic and genetic overlap between autistic traits at the extremes of the general population. *J. Am. Acad. Child Adolesc. Psychiatry* **45**, 1206–1214 (2006).

42. Colvert, E. et al. Heritability of autism spectrum disorder in a UK population-based twin sample. *JAMA Psychiatry* **72**, 415–423 (2015).

43. Bralten, J. et al. Autism spectrum disorders and autistic traits share genetics and biology. *Mol. Psychiatry* **23**, 1205–1212 (2018).

44. St Pourcain, B. et al. ASD and schizophrenia show distinct developmental profiles in common genetic overlap with population-based social communication difficulties. *Mol. Psychiatry* **23**, 263–270 (2017).

45. Robinson, E. B. et al. Genetic risk for autism spectrum disorders and neuropsychiatric variation in the general population. *Nat. Genet.* **48**, 552–555 (2016).

46. Warrier, V. et al. Genome-wide analyses of self-reported empathy: correlations with autism, schizophrenia, and anorexia nervosa. *Transl. Psychiatry* **8**, 35 (2018).

47. Warrier, V. et al. Genome-wide meta-analysis of cognitive empathy: heritability, and correlates with sex, neuropsychiatric conditions and cognition. *Mol. Psychiatry* **23**, 1402–1409 (2017).

48. Cantor, R. M. et al. ASD restricted and repetitive behaviors associated at 17q21.33: genes prioritized by expression in fetal brains. *Mol. Psychiatry* **23**, 993–1000 (2018).

49. Baron-Cohen, S. The hyper-systemizing, assortative mating theory of autism. *Prog. Neuropsychopharmacol. Biol. Psychiatry* **30**, 865–872 (2006).

50. Auyeung, B. et al. The Children’s empathy quotient and systemizing quotient: sex differences in typical development and in autism spectrum conditions. *J. Autism Dev. Disord.* **39**, 1509–1521 (2009).

51. Baron-Cohen, S. et al. Attenuation of typical sex differences in 800 adults with autism vs. 3,900 controls. *PLoS ONE* **9**, e102251 (2014).

52. Asperger, H. in *Autism and Asperger Syndrome* (ed. Frith, U.) 37–92 (Cambridge University Press, Cambridge, 1944).

53. Kanner, L. Autistic disturbances of affective contact. *Nerv. Child J. Psychopathol. Psychother. Ment. Hyg. Guid. Child* **2**, 217–250 (1943).

54. Auyeung, B. et al. Foetal testosterone and the child systemizing quotient. *Eur. J. Endocrinol.* **155**, S123–S130 (2006).

55. Warrier, V., the 23andMe Research team, Bourgeron, T. & Baron-Cohen, S. Genome-wide association study of social relationship satisfaction: significant loci and correlations with psychiatric conditions. *bioRxiv* 196071, https://doi.org/10.1101/196071 (2017).

56. Palmer, E. E. et al. A recurrent de novo nonsense variant in ZSWIM6 results in severe intellectual disability without frontonasal or limb malformations. *Am. J. Hum. Genet.* **101**, 995–1005 (2017).

57. Grove, J. et al. Identification of common genetic risk variants for autism spectrum disorder. *Nat. Genet.* **51**, 431–444 (2019).

58. Okbay, A. et al. Genome-wide association study identifies 74 loci associated with educational attainment. *Nature* **533**, 539–542 (2016).

59. Sniekers, S. et al. Genome-wide association meta-analysis of 78,308 individuals identifies new loci and genes influencing human intelligence. *Nat. Genet.* https://doi.org/10.1038/ng.3869 (2017).

60. Baron-Cohen, S. & Wheelwright, S. J. The empathy quotient: an investigation of adults with Asperger syndrome or high functioning autism, and normal sex differences. *J. Autism Dev. Disord.* **34**, 163–175 (2004).

61. Skuse, D. H., Mandy, W. P. L. & Scourfield, J. Measuring autistic traits: heritability, reliability and validity of the social and communication disorders checklist. *Br. J. Psychiatry* **187**, 568–572 (2005).

62. Zheng, J. et al. LD Hub: a centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis. *Bioinformatics* **33**, 272–279 (2016).

63. Demontis, D. et al. Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder. *Nat. Genet.* **51**, 63–75 (2019).

64. Ripke, S. et al. Genome-wide association analysis identifies 13 new risk loci for schizophrenia. *Nat. Genet.* **45**, 1150–1159 (2013).

65. Stahl, E. et al. Genome wide association study identifies 30 loci associated with bipolar disorder. *bioRxiv* 173062, https://doi.org/10.1101/173062 (2018).

66. Baron-Cohen, S. The extreme male brain theory of autism. *Trends Cogn. Sci.* **6**, 248–254 (2002).

67. Tung, J. Y. et al. Efficient replication of over 180 genetic associations with self-reported medical data. *PLoS ONE* **6**, e24373 (2011).

68. Do, C. B. et al. Web-based genome-wide association study identifies two novel loci and a substantial genetic component for Parkinson's disease. *PLoS Genet.* **7**, e1002141 (2011).

69. Eriksson, N. et al. Novel associations for hypothyroidism include known autoimmune risk loci. *PLoS ONE* **7**, e34442 (2012).

70. Henn, B. M. et al. Cryptic distant relatives are common in both isolated and cosmopolitan genetic samples. *PLoS ONE* **7**, e34267 (2012).

71. 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).

72. Fischbach, G. D. & Lord, C. The Simons simplex collection: a resource for identification of autism genetic risk factors. *Neuron* **68**, 192–195 (2010).

73. Geschwind, D. H. et al. The autism genetic resource exchange: a resource for the study of autism and related neuropsychiatric conditions. *Am. J. Hum. Genet.* **69**, 463–466 (2001).

74. Charman, T. et al. The EU-AIMS Longitudinal European Autism Project (LEAP): clinical characterisation. *Mol. Autism* **8**, 27 (2017).

75. Durand, C. M. et al. Mutations in the gene encoding the synaptic scaffolding protein SHANK3 are associated with autism spectrum disorders. *Nat. Genet.* **39**, 25–27 (2007).

76. Groen, Y., Fuermaier, A. B. M., Den Heijer, A. E., Tucha, O. & Althaus, M. The empathy and systemizing quotient: The psychometric properties of the Dutch version and a review of the cross-cultural stability. *J. Autism Dev. Disord.* **45**, 2848–2864 (2015).

77. Allison, C., Baron-Cohen, S., Stone, M. H. & Muncer, S. J. Rasch modeling and confirmatory factor analysis of the systemizing quotient-revised (SQ-R) scale. *Span. J. Psychol.* **18**, E16 (2015).

78. Ling, J., Burton, T. C., Salt, J. L. & Muncer, S. J. Psychometric analysis of the systemizing quotient (SQ) scale. *Br. J. Psychol.* **100**, 539–552 (2009).

79. Wakabayashi, A. et al. Empathizing and systemizing in adults with and without autism spectrum conditions: cross-cultural stability. *J. Autism Dev. Disord.* **37**, 1823–1832 (2007).
80. Bölte, S., Westerwald, E., Holtmann, M., Freitag, C. M. & Poustka, F. Autistic traits and autism spectrum disorders: the clinical validity of two measures presuming a continuum of social communication skills. *J. Autism Dev. Disord.* **41**, 66–72 (2011).

81. Lord, C. et al. The autism diagnostic observation schedule—generic: a standard measure of social and communication deficits associated with the spectrum of autism. *J. Autism Dev. Disord.* **30**, 205–223 (2000).

82. Gotham, K., Risi, S., Pickles, A. & Lord, C. The autism diagnostic observation schedule: revised algorithms for improved diagnostic validity. *J. Autism Dev. Disord.* **37**, 613–627 (2007).

83. Bodfish, J. W., Symons, F. J., Parker, D. E. & Lewis, M. H. Varieties of repetitive behavior in autism: comparisons to mental retardation. *J. Autism Dev. Disord.* **30**, 237–243 (2000).

84. Lam, K. S. L. & Aman, M. G. The repetitive behavior scale-revised: independent validation in individuals with autism spectrum disorders. *J. Autism Dev. Disord.* **37**, 855–866 (2007).

85. Browning, S. R. & Browning, B. L. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies By use of localized haplotype clustering. *Am. J. Hum. Genet.* **81**, 1084–1097 (2007).

86. Fuchsberger, C., Abecasis, G. R. & Hinds, D. A. minimac2: faster genotype imputation. *Bioinformatics* **31**, 782–784 (2015).

87. Bigdeli, T. B. et al. A simple yet accurate correction for winner's curse can predict signals discovered in much larger genome scans. *Bioinformatics* **32**, 2598–2603 (2016).

88. Delaneau, O., Marchini, J. & Zagury, J.-F. A linear complexity phasing method for thousands of genomes. *Nat. Methods* **9**, 179–181 (2011).

89. Howie, B. N., Donnelly, P. & Marchini, J. A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. *PLoS Genet.* **5**, e1000529 (2009).

90. 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).

91. Bulik-Sullivan, B. K. et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. *Nat. Genet.* **47**, 291–295 (2015).

92. Bulik-Sullivan, B. K. et al. An atlas of genetic correlations across human diseases and traits. *Nat. Genet.* **47**, 1236–1241 (2015).

93. Watanabe, K., Taskesen, E., van Bochoven, A. & Posthuma, D. Functional mapping and annotation of genetic associations with FUMA. *Nat. Commun.* **8**, 1826 (2017).

94. de Leeuw, C. A., Mooij, J. M., Heskes, T. & Posthuma, D. MAGMA: generalized gene-set analysis of GWAS data. *PLoS Comput. Biol.* **11**, 1–19 (2015).

95. Ardlie, K. G. et al. The Genotype-Tissue Expression (GTEx) pilot analysis: multitissue gene regulation in humans. *Science* **348**, 648–660 (2015).

96. Finucane, H. K. et al. Partitioning heritability by functional annotation using genome-wide association summary statistics. *Nat. Genet.* **47**, 1228–1235 (2015).

97. Hormozdiari, F., Kostem, E., Kang, E. Y., Pasaniuc, B. & Eskin, E. Identifying causal variants at loci with multiple signals of association. *Genetics* **198**, 497–508 (2014).

98. Won, H. et al. Chromosome conformation elucidates regulatory relationships in developing human brain. *Nature* **538**, 523–527 (2016).

99. Sklar, P. et al. Large-scale genome-wide association analysis of bipolar disorder identifies a new susceptibility locus near ODZ4. *Nat. Genet.* **43**, 977–983 (2011).

100. Pardiñas, A. F. et al. Common schizophrenia alleles are enriched in mutation-intolerant genes and in regions under strong background selection. *Nat. Genet.* **50**, 381–389 (2018).

101. Duncan, L. E. et al. Significant locus and metabolic genetic correlations revealed in genome-wide association study of anorexia nervosa. *Am. J. Psychiatry* **174**, 850–858 (2017).

102. 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).

103. Okbay, A. et al. Genetic variants associated with subjective well-being, depressive symptoms, and neuroticism identified through genome-wide analyses. *Nat. Genet.* **48**, 624–633 (2016).

104. de Moor, M. H. M. et al. Meta-analysis of genome-wide association studies for personality. *Mol. Psychiatry* **17**, 337–349 (2012).

105. van den Berg, S. M. et al. Meta-analysis of genome-wide association studies for extraversion: findings from the genetics of personality consortium. *Behav. Genet.* **46**, 170–182 (2016).

106. Rietveld, C. A. et al. Common genetic variants associated with cognitive performance identified using the proxy-phenotype method. *Proc. Natl. Acad. Sci. USA* **111**, 13790–13794 (2014).

107. Rietveld, C. A. et al. GWAS of 126,559 individuals identifies genetic variants associated with educational attainment. *Science* **340**, 1467–1471 (2013).

108. Hambrook, D., Tchanturia, K., Schmidt, U., Russell, T. & Treasure, J. Empathy, systemizing, and autistic traits in anorexia nervosa: a pilot study. *Br. J. Clin. Psychol.* **47**, 335–339 (2008).

109. Russell-Smith, S. N., Bayliss, D. M., Maybery, M. T. & Tomkinson, R. L. Are the autism and positive schizotypy spectra diametrically opposed in empathizing and systemizing? *J. Autism Dev. Disord.* **43**, 695–706 (2013).

110. Courty, A. et al. Levels of autistic traits in anorexia nervosa: a comparative psychometric study. *BMC Psychiatry* **13**, 222 (2013).

111. Nettle, D. Empathizing and systemizing: what are they, and what do they contribute to our understanding of psychological sex differences? *Br. J. Psychiatry* **98**, 237–255 (2007).

112. Wakabayashi, A. & Kawashima, H. Is empathizing in the E-S theory similar to agreeableness? The relationship between the EQ and SQ and major personality domains. *Pers. Individ. Differ.* **76**, 88–93 (2015).

113. Deary, V. et al. Genetic contributions to self-reported tiredness. *Mol. Psychiatry* **23**, 609–620 (2018).

114. Ripke, S. et al. Biological insights from 108 schizophrenia-associated genetic loci. *Nature* **511**, 421–427 (2014).

115. Jones, S. E. et al. Genome-wide association analyses in 128,266 individuals identifies new morningness and sleep duration loci. *PLoS Genet.* **12**, e1006125 (2016).

116. Otowa, T. et al. Meta-analysis of genome-wide association studies of anxiety disorders. *Mol. Psychiatry*. https://doi.org/10.1038/mp.2015.197 (2016).

117. Mattheisen, M. et al. Genome-wide association study in obsessive-compulsive disorder: results from the OCGAS. *Mol. Psychiatry* **20**, 337–344 (2015).

118. International Obsessive Compulsive Disorder Foundation Genetics Collaborative (IOCDF-GC) and OCD Collaborative Genetics Association Studies (OCGAS), P. D. et al. Revealing the complex genetic architecture of obsessive-compulsive disorder using meta-analysis. *Mol. Psychiatry* **23**, 1181–1188 (2018).

119. Duncan, L. E. et al. Largest GWAS of PTSD (N = 20,070) yields genetic overlap with schizophrenia and sex differences in heritability. *Mol. Psychiatry* **23**, 666–673 (2017).

120. Allison, C., Baron-Cohen, S., Wheelwright, S. J., Stone, M. H. & Muncer, S. J. Psychometric analysis of the empathy quotient (EQ). *Pers. Individ. Dif.* **51**, 829–835 (2011).

121. Kok, F. M., Groen, Y., Becke, M., Fuermaier, A. B. M. & Tucha, O. Self-reported empathy in adult women with autism spectrum disorders—a systematic mini review. *PLoS ONE* **11**, e0151568 (2016).

122. Nieuwboer, H. A., Pool, R., Dolan, C. V., Boomsma, D. I. & Nivard, M. G. GWIS: Genome-Wide Inferred Statistics for functions of multiple phenotypes. *Am. J. Hum. Genet.* **99**, 917–927 (2016).

123. Grotzinger, A. D. et al. Genomic structural equation modelling provides insights into the multivariate genetic architecture of complex traits. *Nat. Hum. Behav.* **1**, https://doi.org/10.1038/s41562-019-0566-x (2019).

124. Boyle, E. A., Li, Y. I. & Pritchard, J. K. An expanded view of complex traits: from polygenic to omnigenic. *Cell* **169**, 1177–1186 (2017).

125. 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).

Acknowledgements

We are grateful to Michel Nivard and Beate St. Pourcain for their help with the analytical methods. We thank the research participants and employees of 23andMe for making this work possible. We are 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. V.W. was funded by St. John's College, Cambridge, and the Cambridge Commonwealth Trust. This study was funded by grants to SBC from the Medical Research Council, the Wellcome Trust, the Autism Research Trust, the Templeton World Charity Foundation, and to T.B. from the Institut Pasteur, the CNRS, the INSERM, The Fondamental Foundation, the APHP, the BioPsy Labex and the University Paris Diderot. The research was conducted in association with the National Institute for Health Research (NIHR) Collaboration for Leadership in Applied Health Research and Care East of England at Cambridgeshire and Peterborough NHS Foundation Trust. We also received support from the NIHR Cambridge Biomedical Research Centre. We acknowledge with gratitude the generous support of Drs Dennis and Mireille Gillings in strengthening the collaboration between S.B.-C. and T.B., and between Cambridge University and the Institut Pasteur. The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health. Data obtained from 23andMe was supported by the National Human Genome Research Institute of the National Institutes of Health (grant number R44HG006981). The UK Medical Research Council and Wellcome (grant ref: 102215/2/13/2) and the University of Bristol provide core support for ALSPAC. GWAS data were generated by Sample Logistics and Genotyping Facilities at Wellcome Sanger Institute and LabCorp (Laboratory Corporation of America) using support from 23andMe. This publication is the work of the authors who will serve as guarantors for the content of this paper. The iPSYCH (The Lundbeck Foundation Initiative for Integrative Psychiatric Research) team acknowledges funding from The Lundbeck Foundation (grant nos. R102-A9118 and R155-2014-1724), the Stanley Medical Research Institute, the European Research Council
(project no.: 294838), the Novo Nordisk Foundation for supporting the Danish National Biobank resource, and grants from Aarhus and Copenhagen Universities and University Hospitals, including support to the iSEQ Center, the GenomeDK HPC facility, and the CIRRAU Center. The project leading to this application has received funding from the Innovative Medicines Initiative 2 Joint Undertaking (JU) under grant agreement no. 777394. The JU receives support from the European Union's Horizon 2020 research and innovation programme and EFPIA and AUTISM SPEAKS, Autistica, SFARI. We thank the iPSCH-Broad Autism Group and the EU-AIMS LEAP group for sharing data. A full list of the authors and affiliations in the iPSYCH-Broad autism group and the EU-AIMS LEAP group is provided in the Supplementary Information.

**Author contributions**

V.W., T.B., S.B.-C., and D.A.H. conceived and designed the analysis. C.S.L., F.C., R.D., W.D.W., J.B., A.D.B., J.G., G.P., 23andMe Research Team, and D.A.H. collected or contributed the data or analysis tools. V.W., R.T., H.W., F.C., and W.D.W. performed the analysis. V.W., R.T., B.C., D.A.H., T.B., and S.B.C. wrote the paper. D.A.H., T.B., and S. B.-C. supervised the analysis.

**Additional information**

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s42003-019-0558-4.

**Competing interests:** D.A.H. and the 23andMe Research Team are employees of 23andMe, Inc. The remaining authors declare no competing interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_13.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019

**The 23andMe Research Team**

Michelle Agee¹², Babak Alipanahi¹², Adam Auton¹², Robert K. Bell¹², Katarzyna Bryc¹², Sarah L. Elson¹², Pierre Fontanillas¹², Nicholas A. Furlotte¹², Karen E. Huber¹², Aaron Kleinman¹², Nadia K. Litterman¹², Jennifer C. McCreight¹², Matthew H. McIntyre¹², Joanna L. Mountain¹², Carrie A.M. Northover¹², Steven J. Pitts¹², J. Fah Sathirapongsasuti¹², Olga V. Sazonova¹², Janie F. Shelton¹², Suyash Shringarpure¹², Chao Tian¹², Joyce Y. Tung¹², Vladimir Vacic¹² & Catherine H. Wilson¹²