Skip to main content
Advertisement
  • Loading metrics

Polygenic risk for autism spectrum disorder associates with anger recognition in a neurodevelopment-focused phenome-wide scan of unaffected youths from a population-based cohort

Abstract

The polygenic nature and the contribution of common genetic variation to autism spectrum disorder (ASD) allude to a high degree of pleiotropy between ASD and other psychiatric and behavioral traits. In a pleiotropic system, a single genetic variant contributes small effects to several phenotypes or disorders. While analyzed broadly, there is a paucity of research studies investigating the shared genetic information between specific neurodevelopmental domains and ASD. We performed a phenome-wide association study of ASD polygenetic risk score (PRS) against 491 neurodevelopmental subdomains ascertained in 4,309 probands from the Philadelphia Neurodevelopmental Cohort (PNC) who lack an ASD diagnosis. Our main analysis calculated ASD PRS in 4,309 PNC probands using the per-SNP effects reported in a recent genome-wide association study of ASD in a case-control design. In a high-resolution manner, our main analysis regressed ASD PRS against 491 neurodevelopmental phenotypes with age, sex, and ten principal components of ancestry as covariates. Follow-up analyses included in the regression model PRS derived from brain-related traits genetically correlated with ASD. Our main finding demonstrated that 11-17-year old probands with the highest ASD genetic risk were able to identify angry faces (R2 = 1.06%, p = 1.38 × 10−7, pBonferroni-corrected = 1.9 × 10−3). This ability replicated in older probands (>18 years; R2 = 0.55%, p = 0.036) and persisted after covarying with other psychiatric disorders, brain imaging traits, and educational attainment (R2 = 0.2%, p = 0.019). We also detected several suggestive associations between ASD PRS and emotionality and connectedness with others. These data (i) indicate how genetic liability to ASD may influence neurodevelopment in the general population, (ii) reinforce epidemiological findings of heightened ability of ASD cases to predict certain social psychological events based on increased systemizing skills, and (iii) recapitulate theories of imbalance between empathizing and systemizing in ASD etiology.

Author summary

Large-scale genetic studies have identified many regions of the genome associated with autism spectrum disorder that are considered common in the general population. We investigated how the additive effects of these genetic variations associate with neurodevelopment in youths who lack an ASD diagnosis to better understand how genetic risk for ASD may contribute to other aspects of mental health. We uncovered a relationship between greater genetic risk for ASD and more accurate recognition of angry emotions in others, which persists after considering genetic associations with other psychiatric disorders, educational attainment, and brain region volume. This finding is consistent with existing theories of the relationship between ASD genetic liability and a person’s ability to build generalizable and impulse driven models for responding to social phenomena.

Introduction

Autism spectrum disorder (ASD) describes a group of pervasive neurodevelopmental disorders characterized by impaired social and communication skills. ASD typically manifests as a heterogeneous combination of repetitive and restrictive behaviors along with intellectual capabilities ranging from above average intelligence quotient to intellectual disability [1]. ASD affects between 1–1.5% of the general population and is diagnosed more frequently in males than females [2]. ASD often causes serious social impairment for affected individuals, although some function quite normally.

The polygenic nature of ASD and the relatively common frequencies of most ASD risk alleles pose several population-level questions. For example, why genetic liability to ASD correlates with better cognitive function [3]. We might suspect that the effects of ASD risk alleles in unaffected individuals might relate to social cognition, psychopathology, executive control, memory, and/or sensorimotor capabilities. Here we investigate which of these features are associated with the genetic risk for ASD.

Grove, et al. [1] reported an estimated SNP-based observed-scale ASD heritability of ~12% in Europeans. There were robust genetic and phenotype correlations between ASD, cognitive ability, educational attainment, and several behavioral traits [1]. The highly polygenic nature and relatively high contribution of common genetic variation to ASD allude to a high degree of pleiotropy between ASD and other psychiatric and behavioral traits. In a pleiotropic system, a single genetic variant contributes small effects to several phenotypes or disorders. While analyzed broadly, there is a paucity of research studies investigating the shared genetic information between specific neurodevelopmental domains and ASD. Under the hypothesis that polygenic risk for ASD is associated with features of neurotypical development, we attempt to identify these domains in the present study. One way to identify such relationships is a phenome-wide association study (PheWAS). PheWAS are a way to statistically evaluate the pleiotropic nature of a risk variant in relation to other similar or disparately related phenotypes [48]. The PheWAS approach to characterizing pleiotropy also has been applied to investigate how the polygenic architecture of a phenotype associates with other phenotypes [58]. To our knowledge, phenome-wide scans of ASD polygenic risk has not been performed previously.

Here we conducted a phenome-wide scan, testing the relationship between ASD polygenic risk and hundreds of neurodevelopmental phenotypes in young probands from the Philadelphia Neurodevelopmental Cohort (PNC) without an ASD diagnosis [912]. Based on high-resolution polygenic risk scoring (PRS), we observed a significant positive association between ASD PRS and the ability to recognize anger. Conversely, we report a suggestive negative relationship between ASD PRS and the ability to correctly distinguish the age of others. These results together with other suggestive associations highlight several features of human neurodevelopment in the young, such as emotional intelligence, as key pathophysiological targets for ASD etiology.

Results

An overview of the methodological approach used in this study is provided in Fig 1.

thumbnail
Fig 1. Analysis flowchart with relevant research questions in bold text.

https://doi.org/10.1371/journal.pgen.1009036.g001

PNC participant inclusion and description

The Philadelphia Neurodevelopmental Cohort (PNC) consists of youths aged 8–21 stratified by PNC study-organizer-defined age groups (Table 1). PNC was designed to study the genomics of complex pediatric disorders but is not enriched for any specific disorder; the cohort is considered generally healthy. All participants underwent clinical assessment, including a neuropsychiatric structured interview and review of electronic medical records. They were also administered a neuroscience-based computerized neurocognitive battery and a subsample underwent neuroimaging.

thumbnail
Table 1. PNC sample description.

Samples sizes for each age group (defined by Philadelphia Neurodevelopmental Cohort) after quality control and selection of unrelated European samples.

https://doi.org/10.1371/journal.pgen.1009036.t001

This study aimed to discover attributes of neurodevelopment associated with genome-wide additive risk for ASD. To avoid associations confounded by the presence of ASD cases in the target data, PNC probands were excluded from this study if they or their informant answered “yes” to the question “Autism or Pervasive Developmental Disorder—Do/did you have this problem?” (MED291), resulting in removal of fewer than 5% of each proband group [13]. Note that the PNC ASD probands are insufficiently powered to detect an association between ASD PRS and PEITANG (N = 184).

All study participants between the ages of 8 and 18 provided assent and parental approval to participate in the PNC. Probands over 18 years of age provided consent to PNC participation.

Neurodevelopmental trait prediction

ASD PRS was regressed against 491 neurodevelopmental phenotypes in youths lacking an ASD diagnosis. These participants were analyzed in three age groups defined in the data structure of the PNC (young, middle, and adult; Table 1) and consistent with age-of-diagnosis statistics for ASD in the United States (young proband, N = 1,035, age 8–10; middle probands, N = 2,499, age 11–17; and adult probands, N = 775 age ≥18) (Fig 2, Table 1, and S1 Table). Note that age groups defined within the PNC data structure reasonably correspond to different neurodevelopmental timeframes such as early childhood or middle-childhood to adolescence. To maintain consistent nomenclature with the PNC reporting, we refer to age-stratified cohorts by the PNC-assigned proband grouping.

thumbnail
Fig 2. Autism PRS association with PNC phenotypes.

Overview of best-fit models for autism spectrum disorder (ASD) associating with neurodevelopmental phenotypes in adult (AP, Nphenotypes = 324), middle (MP, Nphenotypes = 490), and young (YP, Nphenotypes = 311) probands of the Philadelphia Neurodevelopmental Cohort. (A) Maximum phenotypic variance explained (R2) by the best model fit for each trait given genetic liability to ASD. The dashed horizontal line represents the nominal significance threshold. (B-D) The relationship between binned ASD polygenic risk scores for AP, MP, and YP participants and the most significant phenotype predicted by ASD genetic liability: (B) SIP033 Structural Interview for Prodromal Symptoms: “Has anyone pointed out to you that you are less emotional or connected to people than you used to be?”; (C) PEITANG: Number of correct responses to anger trials during completion of The Penn Emotional Identification Test (PEIT) for recognizing angry emotions; (D) PADT_SAME_PC: Percent of correct responses to test trials with no age difference during completion of the Penn Age Differentiation Test for detecting which face in a face pair appears older. Note the lowest quartile in figures B-D represents the referent.

https://doi.org/10.1371/journal.pgen.1009036.g002

ASD PRS was significantly associated with the emotional intelligence phenotype PEITANG (Penn Emotional Identification Test (PEIT) recognition of angry emotions; z-score = 5.28, R2 = 1.06%, p = 1.38 × 10−7, P-value threshold = 0.4487, NSNPs = 51,452; Fig 2) after applying Bonferroni multiple-testing correction accounting for the number of phenotypes and PRS tested (p = 0.05/13,787 independent tests = 3.63 × 10−6) in subjects aged 11–17 years. The PEIT measures a proband’s ability to recognize and identify a face as either happy, sad, angry, fearful, or neutral in a multiple-choice testing format normalizing the age, sex, and race of faces shown to each proband [14]. The PNC reports PEIT accuracy and speed for each emotion individually and all emotions combined as separate phenotype items. The association of ASD PRS with PEITANG in the middle proband group nominally replicated with adult probands of mean age 19 ± 1.2 years (z-score = 2.10, R2 = 0.55%, p = 0.036). Combining the three age-stratified samples, a stronger association of ASD PRS with PEITANG phenotype was observed (z-score = 5.69, R2 = 0.70%, p = 1.37 × 10−8). When binned by quartiles, there was a positive relationship between increased ASD PRS and the ability to recognize angry faces. The highest quartile of ASD PRS distribution had a 53% increase in correct recognition of angry emotions relative to lowest quartile (βq4vs.q1 = 0.427, pq4vs.q1 = 1.42 × 10−6; Fig 2).

The PEITANG trait was at least nominally significantly correlated with all other emotion recognition accuracies (0.007 ≤ r2 ≤ 0.378, 2.20 × 10−16 ≤ p ≤ 1.86 × 10−5) and reaction times in happy, sad, angry, and fearful emotion trials (0.001 ≤ r2 ≤ 0.007, 2.87 × 10−8 ≤ p ≤ 0.038; S2 Table). To evaluate the relationship between ASD PRS and other emotion recognition tasks independent of anger recognition accuracy, PEITANG was included as a covariate in these models. In models covaried for PEITANG, ASD PRS was significantly associated with the ability to recognize happiness (z-score = -2.41, R2 = 0.221%, p = 0.016, pdiff = 6.41 × 10−4) and neutral emotions (z-score = -2.05, R2 = 0.159%, p = 0.041, pdiff = 0.433; Fig 3). ASD PRS was associated with reaction time for anger trials (z-score = -2.08, R2 = 0.138%, p = 0.038), fear trials (z-score = -2.06, R2 = 0.161%, p = 0.018), and happy trials (z-score = 2.21, R2 = 0.356%, p = 0.027; Fig 3); however, covarying with PEITANG did not significantly change the relationship between ASD PRS and reaction times for any trials (p>0.05). A related PNC trait, the Penn Emotion Differentiation Test (PEDT), detects a participant’s ability to decide which of two faces is more severe. All PEDT trials were weakly associated with ASD PRS (S1 Fig).

thumbnail
Fig 3. Predicting emotion recognition with ASD PRS.

Association between facial emotion capabilities of the middle proband group (ages 11 to 17) of the Philadelphia Neurodevelopmental Cohort using polygenic risk (PRS) for autism spectrum disorder. The x-axis shows correctness and response time for each facial emotion using two neurocognitive instruments: The Penn Emotion Identification Test (shaded colors: original results; tinted colors: results covaried for the effects of PEITANG). Significance is indicated by * for p < 0.05 and ** for FDR Q < 0.05.

https://doi.org/10.1371/journal.pgen.1009036.g003

Traits associated with ASD PRS using a suggestive threshold are detailed in the S1 Text, S2 and S3 Figs, and S3 Table.

Influence of related brain traits

We next tested whether the association between ASD PRS and PEITANG was attributed to the shared genetic liability of ASD with other brain-related phenotypes. To do this we considered the association between PEITANG and PRS calculated from brain imaging phenotypes [15], other psychiatric disorders (i.e., attention deficit hyperactivity disorder (ADHD) [16], anorexia nervosa [17], bipolar disorder [18], major depressive disorder [19], Tourette syndrome [20], obsessive compulsive disorder [21], and schizophrenia [22]), and educational attainment [23] (Table 2 and S4 Table). Among these traits, we selected fourteen brain imaging phenotypes and four psychiatric disorders that were significantly genetically correlated with ASD after multiple testing correction [24]. To remove genetically redundant traits, covariates for PRS were selected from genetically correlated trait pairs based on which had the highest SNP-based heritability z-score (S4 Fig). This approach identified six brain imaging phenotypes (left cerebellar white matter volume, mean diffusion tensor mode (MO) in sagittal stratum, mean L1 (i.e., strength of diffusion along the L1 principal axis of diffusion magnetic resonance imaging (dMRI)) in anterior limb of internal capsule (right), mean L3 (i.e., strength of diffusion along the L1 principal axis of dMRI) in uncinate fasciculus (left), right hemisphere central sulcus area, and left hemisphere cuneus gyrus thickness, two psychiatric disorders (ADHD and schizophrenia), and educational attainment. For complete description of brain image acquisition and brain mapping, refer to Miller, et al. 2016 [25] and Elliot, et al. 2018 [15]. ASD PRS was significantly associated with the ability to recognize angry emotions in the middle proband group when considering one covariate at a time and when including all covariates in the model (Fig 4). The effects of ASD PRS on the relationship between each covariate and PEITANG also was considered. At matched GWS thresholds, the PRS for schizophrenia, educational attainment, and all six brain imaging phenotypes maintained significant associations with PEITANG after covarying for ASD PRS (Table 2). The relationship between the ADHD PRS and PEITANG was not independent of the effects of ASD PRS.

thumbnail
Fig 4. Robustness of the ASD and PEITANG relationship.

Association between the ability to recognize angry faces in the middle proband group (ages 11–17) using the polygenic risk for autism spectrum disorder (ASD) covaried for age, sex, ten principal components, and the PRS for four brain imaging phenotypes, attention deficit hyperactivity disorder, schizophrenia, and educational attainment. (A) PRS p-value across a range of genome-wide significance (GWS) thresholds; the maximum PRS before and after covarying for all brain and psychiatry traits are labeled. (B) The positive correlation between quartiles of ASD polygenic risk and the number of correct responses to anger recognition trials in the Penn Emotional Intelligence Test; the lowest PRS quartile represents the referent.

https://doi.org/10.1371/journal.pgen.1009036.g004

thumbnail
Table 2. Relationship between PEITANG and PRS of additional phenotypes.

Best model fit PRS of psychiatric disorder, educational attainment, and brain image-derived phenotype polygenic risk scores associated with the PEITANG phenotype in the Philadelphia Neurodevelopmental Cohort middle proband group (aged 11–17 years old).

https://doi.org/10.1371/journal.pgen.1009036.t002

Informant perception

Self-reported behavioral attributes (i.e., the person’s perception of their own trait or behavior) may deviate from others’ perception of that same attribute [26]. In the middle proband age group, we evaluated associations between ASD PRS, informant perception of traits, and the differences between informant-reported and self-reported traits. After multiple testing correction based on the number of phenotypes tested (N = 507; FDR Q < 0.05), 24 informant-reported phenotypes were significantly associated with ASD PRS. The most significantly associated informant-reported phenotype was PHB014 (informant-reported PHB014: z-score = -3.86, R2 = 1.20%, p = 1.18 × 10−4) and this relationship agreed with the proband-reported effect direction (PEITANG-covaried proband-reported PHB014: z-score = -3.01, R2 = 0.156%, p = 0.001, pdiff = 3.25 × 10−4; S5 Table). The PHB014 trait contains answers to the question “Thinking about all of the time that you were afraid of (insert worst fear), whether or not you actually faced it, how long did this fear last (Months)?

Discussion

ASD PRS was at least nominally associated with several aspects of emotion recognition in healthy subjects aged 11 to 17 years old (S5 Fig). Thus, our results tie together the genetic risk for ASD pathology and related traits ascertained in the normal population. We showed that ASD PRS was positively associated with the ability to recognize angry faces (i.e., PEITANG). Emotion recognition deficits are typically considered hallmark attributes of ASD cases [27] although this is often contested in the literature on the basis of different diagnostic instrument accuracy, ASD case severity heterogeneity, and eye movement tracking frequency and accuracy [28]. After covarying all other emotion recognition measure with PEITANG, ASD PRS remained associated with poor recognition of happy and neutral faces in others.

We hypothesize that this paradoxical relationship stems from positive associations between ASD and the ability to predict social phenomena more generally. Recently characterized by Gollwitzer and Bargh [29], people perceptive skills (e.g., “People prefer to interact with people who are different than them, rather than similar to them?” True-False) and person perceptive skills (e.g., “This person appears happy when interacting with others.” True-False) share some epidemiological correlates but are considered distinct social psychological skillsets. Gollwitzer, et al. demonstrated that the inability of ASD cases to perceive individual intention, emotion, and thoughts is uniquely distinct (and indeed oppositely correlated) from the ability of ASD cases to predict social phenomena more generally [30]. Importantly, this observation was observed across data from 104 countries suggesting that culturally-driven social norms do not confound this skill. The authors propose the skill of systemizing as a likely origin for their observations. Systemizing describes the ability to predict generalizable rules and regularities of a system and is positively associated with cognitive ability and cognitive curiosity. Our detected increased ability of individuals with high genetic risk for ASD, but no ASD diagnosis, to recognize negative emotions in others is reflective of their superior ability to build generalizable “if a, then b” rules (e.g., “if the face is red, the brow is furrowed, and the jaw is clenched, then this person is angry”). After covarying for schizophrenia, ADHD, educational attainment, and six brain imaging phenotypes, ASD PRS remained positively associated with the PEITANG trait. This observation suggests that although PRS for other traits (e.g., schizophrenia) independently predicted recognition of angry faces with greater magnitude than ASD PRS, the relationship between ASD PRS and the ability to recognize anger in others is independent of the effects of highly genetically correlated traits with greater sample sizes. PRS for schizophrenia, educational attainment, cerebellar regions, the uncinate fasciculus, and cuneus were positively associated with recognizing angry faces. This is unsurprising given the relationship between education and psychiatric conditions and perhaps reinforces human and animal models of the role of oxytocin, stress, and anger on the brain [31]. This persistent association between genetic risk for several brain phenotypes, including strong association with neuropsychiatric conditions, suggests that anger identification may be a common feature across neuropsychiatry more generally than just ASD. However, most of the effects observed appear to be independent. This suggests a complex scenario where the association of genetic risk across psychiatric disorders with anger recognition is due to independent mechanisms rather than to a common shared pathway. Additionally, there are notable power differences among GWAS used to generate the PRS. Indeed, the stronger effect observed for SCZ PRS is likely due to it being the most informative GWAS dataset used to generate the PRS.

Three previous studies have investigated the genetic relationships between emotional intelligence, cognitive empathy, and ASD [3234]. These previous analyses relied on ASD GWAS data from substantially smaller sample sizes than those used herein (N = 46,351 vs. N = 10,610 in [3234]). Nevertheless, our data replicate previous nominally significant findings between ASD PRS and negative emotion recognition tasks (e.g., sadness and fear). To our knowledge, this is the first observation of a relationship between ASD PRS (based on a total sample size over 4 times larger than previous investigations) and angry emotion recognition using genetic data in 11-17-year-old youths.

It is intriguing that the strong significant relationship between ASD PRS and PEITANG in 11-17-year-old probands replicated in adult probands but not young probands. According to our power calculation (see Materials and methods), both adult proband and young proband groups have sample sizes large enough to replicate the PEITANG association at nominal significance level (p < 0.05), but not considering a multiple testing correction significance level (p < 3.63 × 10−6; required N = 1,070). One possible explanation for lack of replication across development may be related to how sex balance across proband groups distributes the relationship between puberty, hormones, and social cognition. Females are better at emotion recognition tasks than males [35,36] and there is evidence that late-pubescent and post-pubescent individuals are generally better able to recognize emotions than those in early puberty or prepubescents (e.g., the PNC young proband group) [37]. Though adjusted for both age and sex, our analysis includes a larger number of females in both the middle (52% female) and adult (55% female) proband groups relative to the young proband group (44% female). The greater proportion of females in adult and middle proband age groups and the earlier onset of puberty in females likely increases the number of late-pubescent and even post-pubescent individuals in the middle and adult proband groups. We hypothesize that in combination, later onset puberty in males relative to females and the greater number of males in the YP group contributed to masking the relationship between ASD PRS and PEITANG. An additional explanation may relate to documented developmental age-specific correlative and polygenic predictive effects of ASD and other psychiatric disorders on social communication deficits [38,39]. PEITANG is not a social communication deficit and the association between ASD PRS and communication deficits appears to have a negative linear relationship with age. Therefore, our findings reinforce that age is an important moderator of the effect of ASD PRS on behavior in unaffected individuals.

This study has four main limitations. First, the PNC evaluates broad aspects of human neurodevelopment and with respect to social cognition only investigates facial stimuli at the exclusion of gesture, vocalization, and complex social situations. Future work may investigate the relationship between ASD PRS and higher resolution measures of social cognition, including those related to people perception. Second, ASD cases were removed prior to PRS but the remaining proband groups likely contain disproportionate ratios of disorders genetically correlated with ASD which may bias the observed predictions in a direction that favors ASD PRS strongly predicting a PNC trait (i.e., false positives). Given the strength of the relationship between ASD PRS and PEITANG, this potential bias likely lies within the phenotypes predicted by ASD PRS with nominal significance (Fig 1). Even though these biases may influence the specific traits predicted by ASD PRS in each proband group, the strongest relationship (ASD PRS and PEITANG) was replicable across age groups. Third, we were not able to find an adequately powered sample independent from the PNC cohort to externally replicated the PEITANG association (see Materials and methods). Though replicated in the PNC adult proband group, our results have not been replicated externally and should be interpreted as such until verified in a sufficiently large cohort of individuals tested with the PEIT. Future work should aim to confirm this relationship in the middle proband age group as well as across developmental age groups. Fourth, unfortunately we did not identify a publicly available dataset using the PEIT that also included adult participants to investigate how genetic liability to ASD is associated with anger recognition and other neurodevelopmental traits across the lifespan.

In this study, we demonstrated that ASD PRS is associated with the ability to recognize angry faces in ASD-unaffected individuals; however, ASD PRS was not associated with the ability to differentiate the intensity of emotions. These results were observed in a cohort without enrichment for specific neurodevelopmental phenotypes and self- and informant-reported ASD cases were removed from our analyses [13]. These data support further investigation of how ASD risk alleles relate to psychopathology and neurodevelopment in healthy and ASD patients. Considering epidemiological data regarding emotion recognition, future work may investigate how ASD PRS predicts correct recognition and reaction time for emotion recognition and differentiation tasks in ASD cases and the relationship between ASD genetic risk and features of generalized social cognition skills (e.g., social loafing, misattribution, deindividuation, self-serving bias). Additionally, our findings appear to indicate that ASD PRS association is independent from the genetic risk of other psychiatric and brain-related traits. Further studies based on multivariate statistics will be needed to expand this result, characterizing the molecular mechanisms shared across the psychiatric spectrum and those that are disease-specific.

Materials and methods

Genetic data

Genome-wide association study (GWAS) summary statistics for ASD (18,382 cases and 27,969 controls of European descent from Grove, et al. 2019) and additional psychiatric disorders (S4 Table) were obtained via the Psychiatric Genomics Consortium (PGC; available at https://www.med.unc.edu/pgc/results-and-downloads). GWAS summary statistics for 3,144 brain image-derived phenotypes were obtained from the Oxford Brain Imaging Genetics project (BIG; available at http://big.stats.ox.ac.uk/) [15,25].

Neurodevelopmental phenotype and genotype data for 9,267 youths aged 8–21 were obtained from the PNC (Neurodevelopmental Genomics: Trajectories of Complex Phenotypes, dbGaP phs000607.v3.p2). Comprehensive details of the clinical and cognitive assessments have been reported previously [912]. Briefly, the Children’s Hospital of Philadelphia (CHOP) and University of Pennsylvania deeply phenotyped youth participants visiting CHOP or an affiliated clinic for routine visit and volunteered to participate in genomic studies of complex pediatric disorders. The cohort is considered generally healthy and is not enriched for any specific disorder, behavior, or trait.

Clinical testing for each participant included (1) GOASSES (a modified version of the Kiddie-Schedule for Affective Disorders and Schizophrenia) to identify timeline of life events, demographic and medical history, Global Assessment of Functioning, and general interviewer observations, (2) a psychopathology symptom and criterion-related assessment of mood disorders, anxiety disorders, behavioral disorders, psychosis spectrum, eating disorders, suicidal thinking and behavior, and treatment history, and (3) an abbreviated form of the Family Interview for Genetics Studies to assess major domains of psychopathology in the proband’s first-degree relatives.

A neurocognitive battery was performed for each PNC proband. The battery yields measures of accuracy and speed for domains of executive control functions (abstraction, attention, and working memory), episodic memory (verbal, facial, and spatial processing), social cognition (emotion identification, emotion intensity differentiation, and age differentiation), and sensorimotor and motor speed. These neurodevelopmental domains are evaluated using the Penn Conditioning Exclusion Test, Penn Continuous Performance Test, Letter N-back Test, Penn Word Memory Test, Penn Face Memory Test, Visual Object Learning Test, Penn Verbal Reasoning, Penn Line Orientation Test, Penn Emotion Identification Test, and Penn Emotion Differentiation Test.

For a complete list of neurodevelopmental domains assessed and the specific features of each domain tested in the neurocognitive battery, see https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000607.v3.p2.

The PNC assigned each participant to a proband age group based on age at recruitment and sampling: young, middle, or adult (Table 1). ASD in the United States often is diagnosed after age 9 and rarely in youth or young adulthood, even if the individual is symptomatic. The PNC-defined age groups coincide with this trend so we maintained the age-stratified architecture of PNC probands and investigated each age group independently (Table 1). Phenotypes were included in our analysis if they had ≥ 500 participants and effective sample size ≥ 50 (i.e., the number of cases for each phenotype necessary to detect an effect weighted by the proportion of the binary categories with respect to that variable).

GWAS quality control and trait inclusion

Pre-imputation quality control was performed in plink-v1.9 following PGC analysis pipelines specifically designed to handle large-scale datasets consisting of multiple genotyping platforms (see https://sites.google.com/a/broadinstitute.org/ricopili/preimputation-qc). Individuals of European descent were confirmed with genetic information via principle component analysis using SNPs shared across all bead-chips (38,739 SNPs) and the 1000 Genomes Project reference panel for populations with European ancestry (N = 503). For sample pairs with cryptic relatedness (PI-HAT > 0.2), the sample with more informative phenotypes was retained. Imputation was performed for 4,309 unrelated individuals of European ancestry using SHAPEIT for pre-phasing [40], IMPUTE2 for imputation [41], and the human 1000 Genomes Project Phase 3 as a reference panel. Samples were stratified by proband group assignment in the PNC for PRS analyses (Table 1).

Polygenic risk scoring and phenome-wide association

PRS based on ASD GWAS summary statistics were regressed against phenotypes in three PNC proband groups. PRS was performed in PRSice v2 [42] with the default settings (i.e., SNPs were clumped in 250-kb windows on either side of target SNPs using clump-r2 and clump-p thresholds of 0.1 and 1, respectively). High resolution PRS were computed using GWAS variants with GWAS association p-values ranging from p  =  0.0001 to p  = 1 in 0.000005 increments. PRS were weighted by their effect size in the independent GWAS (e.g., ASD). All PRS were covaried for age, sex, and the first ten PCs. The most significant PRS also was covaried using the PRS of UKB BIG brain imaging phenotypes and other psychiatric disorders using the same PRS procedure described above (Table 2 and S4 Table). To select additional phenotypes to include as covariates in our PRS calculation, SNP-based observed-scale heritability and genetic correlation for PGC psychiatric disorders and BIG brain imaging phenotypes were determined using the Linkage Disequilibrium Score Regression (LDSC) method (24). Covariates from significantly genetically correlated trait pairs were selected for PRS based on which phenotype from each trait pair had the higher SNP-based heritability z-score.

To our knowledge there is no sample overlap between PNC participants and those samples included in the GWAS of ASD or other phenotypes used herein.

To evaluate how traits vary with increased ASD polygenic risk, the individual genetic liabilities to ASD were binned into four quartiles and regressed against PNC traits using a generalized linear model including the same covariates as the respective PRS prediction. The phenotype values for each quartile were then compared to the reference quartile one-by-one with quartile status of ASD risk as a predictor of target phenotype in a regression. We applied a multiple testing correction accounting for the number of independent tests performed in each age group, as determined by spectral decomposition of non-parametric correlations (Spearman’s rho) between PRS tested and PNC phenotypes [43]. The nominal significance threshold of 0.05 was adjusted for the sum of independent tests from each age group (13,787 independent tests). Additionally, we considered a suggestive significance threshold based on False Discovery Rate (FDR Q = 5%) accounting for the number of phenotypes tested (N = 491).

Statistical power

The avengeme R package (Additive Variance Explained and Number of Genetics Effects Method of Estimation) was used to estimate the sample size needed to reach an adequate statistical power to observe significant associations surviving a nominal significance threshold (p<0.05) and the multiple testing correction significance threshold (p<3.63×10−6). The power calculation was conducted considering the PEITANG (Penn Emotion Identification Test correct responses to anger trials) phenotype as the reference, because it showed an association with autism spectrum disorder (ASD) polygenic risk score (PRS) that survived multiple testing correction. Based on a weighted effects model (i.e., the number of risk alleles carried are weighted by per-SNP effect estimates), and number of SNPs contributing to the model, a sample of 281 participants is required to achieve 80% power to detect a nominally significant association (p<0.05) between ASD PRS and PEITANG. With the same model specifics, a sample size of 1,070 is required to achieve 80% power to detect a significant association between ASD PRS and PEITANG that survives multiple testing correction (p<3.63x10-6).

Phenotype correlation

Spearman’s correlation between phenotype measurements of the PNC were calculated in R studio using the rcorr function of the Hmisc library. Correlation p-values were adjusted for the number of tests performed using FDR Q < 0.05.

Supporting information

S1 Text. Supplemental results–Neurodevelopmental trait association.

https://doi.org/10.1371/journal.pgen.1009036.s001

(DOCX)

S1 Fig. Prediction of facial emotion differentiation capabilities of the middle proband group (ages 11 to 17) of the Philadelphia Neurodevelopmental Cohort using polygenic risk (PRS) for autism spectrum disorder.

The x-axis shows correctness and response time for each facial emotion using two neurocognitive instruments: The Penn Emotion Differentiation Test Correctness and Reaction Time. Significance is indicated by * for p < 0.05.

https://doi.org/10.1371/journal.pgen.1009036.s002

(DOCX)

S2 Fig. Fifty-nine traits phenotypically correlated with PEITANG and at least nominally predicted by polygenic risk (PRS) for autism spectrum disorder (ASD) in the Philadelphia Neurodevelopmental Cohort middle proband group (age 11–17).

Neurodevelopmental traits are labeled if the correlation with PEITANG and PRS survive multiple testing correction; these are PEIT_CR: Penn Emotion Identification Test total correct responses for all trials; GAD018: Generalized Anxiety Disorder: Did you feel any of the following physical symptoms when you worried the most: irritability (feeling easily annoyed)?; PHB012: Specific Phobia: Thinking about all of the time that you were afraid of (insert worst fear), whether or not you actually faced it, how long did this fear last? (Days); PHB014: Specific Phobia: Thinking about all of the time that you were afraid of (insert worst fear), whether or not you actually faced it, how long did this fear last? (Months).

https://doi.org/10.1371/journal.pgen.1009036.s003

(DOCX)

S3 Fig. Z-score converted difference in polygenic risk score (PRS) coefficient for 30 phenotypes suggestively predicted by polygenic risk for autism spectrum disorder in the middle proband group of the Philadelphia Neurodevelopmental Cohort.

Two differences were evaluated: (1; circles) the coefficient difference between original best-fit PRS and the same PRS covaried for the effects of PEITANG (i.e., the same genome-wide significance threshold in both instances) and (2; triangles) the coefficient difference between original best-fit PRS and the best-fit PRS after covarying for the effects of PEITANG (i.e., the GWS may be different in each instance). Note that in several instances the covaried best-fit PRS is the same as the matched PRS from the original test. Filled-in shapes indicate statistically significant differences between original and covaried PRS coefficients (false discovery rate Q < 0.05).

https://doi.org/10.1371/journal.pgen.1009036.s004

(DOCX)

S4 Fig. Selection of brain image-derived phenotypes as covariates in the polygenic risk score analyses between autism spectrum disorder (ASD) and neuropsychiatric traits in the young.

(A) Genetic correlation between ASD 3,199 brain imaging phenotypes from the Brain Imaging Genetics project and (B) genetic correlation between 14 brain imaging phenotypes nominally genetically correlated with ASD.

https://doi.org/10.1371/journal.pgen.1009036.s005

(DOCX)

S5 Fig. Summary of the relationship between autism spectrum disorder (ASD) polygenic risk score (PRS) and anger recognition (PEITANG; Penn Emotion Identification Test correct responses to anger trials) across Philadelphia Neurodevelopmental Cohort (PNC) proband age groups.

https://doi.org/10.1371/journal.pgen.1009036.s006

(DOCX)

S1 Table. Phenotypes meeting suggestive threshold for significant prediction by polygenic risk for autism spectrum disorder in the young (YP), middle (MP), and adult (AP) proband groups.

https://doi.org/10.1371/journal.pgen.1009036.s007

(DOCX)

S2 Table. Spearman’s correlation (r2) and significance (p-value) between PEITANG and other measures recorded by the Penn Emotion Identification Test (PEIT).

https://doi.org/10.1371/journal.pgen.1009036.s008

(DOCX)

S3 Table. Traits significantly correlated (r2 and correlation p) with PEITANG in the middle proband group of the Philadelphia Neurodevelopmental Cohort and at least nominally associated with polygenic risk (PRS) for autism spectrum disorder (PRS R2, z-score, and PRS p).

https://doi.org/10.1371/journal.pgen.1009036.s009

(DOCX)

S4 Table. Samples size and SNP-based observed-scale heritability (h2) for genome-wide association study summary statistics used in this study.

https://doi.org/10.1371/journal.pgen.1009036.s010

(DOCX)

S5 Table. Phenotypes associated with polygenic risk for autism spectrum disorder based on perceived behavior of middle age group (11 to 17 years old) informants from the Philadelphia Neurodevelopmental Cohort.

https://doi.org/10.1371/journal.pgen.1009036.s011

(DOCX)

Acknowledgments

The authors thank the Philadelphia Neurodevelopmental Cohort (PNC) data collection and curation team.

References

  1. 1. Grove J, Ripke S, Als TD, Mattheisen M, Walters RK, Won H, et al. Identification of common genetic risk variants for autism spectrum disorder. Nat Genet. 2019;51(3):431–44. pmid:30804558
  2. 2. Kim YS, Leventhal BL, Koh YJ, Fombonne E, Laska E, Lim EC, et al. Prevalence of autism spectrum disorders in a total population sample. Am J Psychiatry. 2011;168(9):904–12. pmid:21558103
  3. 3. Polimanti R, Gelernter J. Widespread signatures of positive selection in common risk alleles associated to autism spectrum disorder. PLoS Genet. 2017;13(2):e1006618. pmid:28187187
  4. 4. Easey K, Haan E, Schellas L, Sallis H, Wootton R, Munafò M, et al. P11 The association of alcohol PRS on mental health phenotypes: a PheWAS in the avon longitudinal study of parents and children (ALSPAC). Journal of Epidemiology and Community Health. 2019;73(Suppl 1):A76–A7.
  5. 5. Fritsche LG, Beesley LJ, VandeHaar P, Peng RB, Salvatore M, Zawistowski M, et al. Exploring various polygenic risk scores for skin cancer in the phenomes of the Michigan genomics initiative and the UK Biobank with a visual catalog: PRSWeb. PLoS Genet. 2019;15(6):e1008202. pmid:31194742
  6. 6. Leppert B, Millard LA, Riglin L, Smith GD, Thapar A, Tilling K, et al. A cross-disorder MR-pheWAS of 5 major psychiatric disorders in UK Biobank. bioRxiv. 2019:634774.
  7. 7. Richardson TG, Harrison S, Hemani G, Davey Smith G. An atlas of polygenic risk score associations to highlight putative causal relationships across the human phenome. Elife. 2019;8.
  8. 8. Zheutlin AB, Dennis J, Karlsson Linner R, Moscati A, Restrepo N, Straub P, et al. Penetrance and Pleiotropy of Polygenic Risk Scores for Schizophrenia in 106,160 Patients Across Four Health Care Systems. Am J Psychiatry. 2019;176(10):846–55. pmid:31416338
  9. 9. Calkins ME, Merikangas KR, Moore TM, Burstein M, Behr MA, Satterthwaite TD, et al. The Philadelphia Neurodevelopmental Cohort: constructing a deep phenotyping collaborative. J Child Psychol Psychiatry. 2015;56(12):1356–69. pmid:25858255
  10. 10. Calkins ME, Moore TM, Merikangas KR, Burstein M, Satterthwaite TD, Bilker WB, et al. The psychosis spectrum in a young U.S. community sample: findings from the Philadelphia Neurodevelopmental Cohort. World Psychiatry. 2014;13(3):296–305. pmid:25273303
  11. 11. Glessner JT, Reilly MP, Kim CE, Takahashi N, Albano A, Hou C, et al. Strong synaptic transmission impact by copy number variations in schizophrenia. Proc Natl Acad Sci U S A. 2010;107(23):10584–9. pmid:20489179
  12. 12. Gur RC, Calkins ME, Satterthwaite TD, Ruparel K, Bilker WB, Moore TM, et al. Neurocognitive growth charting in psychosis spectrum youths. JAMA Psychiatry. 2014;71(4):366–74. pmid:24499990
  13. 13. Merikangas AK, Calkins ME, Bilker WB, Moore TM, Gur RC, Gur RE. Parental Age and Offspring Psychopathology in the Philadelphia Neurodevelopmental Cohort. J Am Acad Child Adolesc Psychiatry. 2017;56(5):391–400. pmid:28433088
  14. 14. Kohler CG, Turner T, Stolar NM, Bilker WB, Brensinger CM, Gur RE, et al. Differences in facial expressions of four universal emotions. Psychiatry Res. 2004;128(3):235–44. pmid:15541780
  15. 15. Elliott LT, Sharp K, Alfaro-Almagro F, Shi S, Miller KL, Douaud G, et al. Genome-wide association studies of brain imaging phenotypes in UK Biobank. Nature. 2018;562(7726):210–6. pmid:30305740
  16. 16. Demontis D, Walters RK, Martin J, Mattheisen M, Als TD, Agerbo E, et al. Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder. Nat Genet. 2019;51(1):63–75. pmid:30478444
  17. 17. Duncan L, Yilmaz Z, Gaspar H, Walters R, Goldstein J, Anttila V, et al. Significant Locus and Metabolic Genetic Correlations Revealed in Genome-Wide Association Study of Anorexia Nervosa. Am J Psychiatry. 2017;174(9):850–8. pmid:28494655
  18. 18. Bipolar Disorder and Schizophrenia Working Group of the Psychiatric Genomics Consortium (2018) Genomic Dissection of Bipolar Disorder and Schizophrenia, Including 28 Subphenotypes. Cell.173(7):1705–15 e16. pmid:29906448
  19. 19. Wray NR, Ripke S, Mattheisen M, Trzaskowski M, Byrne EM, Abdellaoui A, et al. Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nat Genet. 2018;50(5):668–81. pmid:29700475
  20. 20. Yu D, Sul JH, Tsetsos F, Nawaz MS, Huang AY, Zelaya I, et al. Interrogating the Genetic Determinants of Tourette's Syndrome and Other Tic Disorders Through Genome-Wide Association Studies. Am J Psychiatry. 2019;176(3):217–27. pmid:30818990
  21. 21. International Obsessive Compulsive Disorder Foundation Genetics Collaborative (IOCDF-GC) and OCD Collaborative Genetics Association Studies (OCGAS. Revealing the complex genetic architecture of obsessive-compulsive disorder using meta-analysis. Mol Psychiatry. 2018;23(5):1181–8. pmid:28761083
  22. 22. Schizophrenia Working Group of the Psychiatric Genomics Consortium (2014) Biological insights from 108 schizophrenia-associated genetic loci. Nature.511(7510):421–7. pmid:25056061
  23. 23. Lee JJ, Wedow R, Okbay A, Kong E, Maghzian O, Zacher M, et al. Gene discovery and polygenic prediction from a genome-wide association study of educational attainment in 1.1 million individuals. Nat Genet. 2018;50(8):1112–21. pmid:30038396
  24. 24. Bulik-Sullivan B, Finucane HK, Anttila V, Gusev A, Day FR, Loh PR, et al. An atlas of genetic correlations across human diseases and traits. Nat Genet. 2015;47(11):1236–41. pmid:26414676
  25. 25. Miller KL, Alfaro-Almagro F, Bangerter NK, Thomas DL, Yacoub E, Xu J, et al. Multimodal population brain imaging in the UK Biobank prospective epidemiological study. Nat Neurosci. 2016;19(11):1523–36. pmid:27643430
  26. 26. Baker LA, Jacobson KC, Raine A, Lozano DI, Bezdjian S. Genetic and environmental bases of childhood antisocial behavior: a multi-informant twin study. J Abnorm Psychol. 2007;116(2):219–35. pmid:17516756
  27. 27. Lai MC, Lombardo MV, Baron-Cohen S. Autism. Lancet. 2014;383(9920):896–910. pmid:24074734
  28. 28. Kou J, Le J, Fu M, Lan C, Chen Z, Li Q, et al. Comparison of three different eye-tracking tasks for distinguishing autistic from typically developing children and autistic symptom severity. bioRxiv. 2019:547505.
  29. 29. Gollwitzer A, Bargh JA. Social psychological skill and its correlates. Social Psychology. 2018;49(2):88–102.
  30. 30. Gollwitzer A, Marte C, McPartland JC, Bargh JA. Autism spectrum traits predict higher social psychological skill. Proc Natl Acad Sci U S A. 2019.
  31. 31. Smith KW, Balkwill LL, Vartanian O, Goel V. Syllogisms delivered in an angry voice lead to improved performance and engagement of a different neural system compared to neutral voice. Front Hum Neurosci. 2015;9:273. pmid:26029089
  32. 32. Coleman JRI, Lester KJ, Keers R, Munafo MR, Breen G, Eley TC. Genome-wide association study of facial emotion recognition in children and association with polygenic risk for mental health disorders. Am J Med Genet B Neuropsychiatr Genet. 2017;174(7):701–11. pmid:28608620
  33. 33. Warrier V, Baron-Cohen S. Genetic contribution to 'theory of mind' in adolescence. Sci Rep. 2018;8(1):3465. pmid:29472613
  34. 34. Warrier V, Grasby KL, Uzefovsky F, Toro R, Smith P, Chakrabarti B, et al. Genome-wide meta-analysis of cognitive empathy: heritability, and correlates with sex, neuropsychiatric conditions and cognition. Mol Psychiatry. 2018;23(6):1402–9. pmid:28584286
  35. 35. Lausen A, Broering C, Penke L, Schacht A. Hormonal and modality specific effects on males’ emotion recognition ability. bioRxiv. 2019:791376.
  36. 36. Thompson AE, Voyer D. Sex differences in the ability to recognise non-verbal displays of emotion: a meta-analysis. Cogn Emot. 2014;28(7):1164–95. pmid:24400860
  37. 37. Mancini G, Biolcati R, Agnoli S, Andrei F, Trombini E. Recognition of Facial Emotional Expressions Among Italian Pre-adolescents, and Their Affective Reactions. Front Psychol. 2018;9:1303. pmid:30123150
  38. 38. St Pourcain B, Robinson EB, Anttila V, Sullivan BB, Maller J, Golding J, et al. ASD and schizophrenia show distinct developmental profiles in common genetic overlap with population-based social communication difficulties. Mol Psychiatry. 2018;23(2):263–70. pmid:28044064
  39. 39. Thapar A, Riglin L. The importance of a developmental perspective in Psychiatry: what do recent genetic-epidemiological findings show? Mol Psychiatry. 2020.
  40. 40. Delaneau O, Marchini J, Zagury JF. A linear complexity phasing method for thousands of genomes. Nat Methods. 2011;9(2):179–81. pmid:22138821
  41. 41. Howie B, Marchini J, Stephens M. Genotype imputation with thousands of genomes. G3 (Bethesda). 2011;1(6):457–70.
  42. 42. Euesden J, Lewis CM, O'Reilly PF. PRSice: Polygenic Risk Score software. Bioinformatics. 2015;31(9):1466–8. pmid:25550326
  43. 43. Nyholt DR. A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other. Am J Hum Genet. 2004;74(4):765–9. pmid:14997420