Where are the strongest associations between autistic traits and traits of ADHD? evidence from a community-based twin study.

Autism spectrum conditions (ASC) and attention-deficit/hyperactivity disorder (ADHD) regularly co-occur. Twin studies increasingly indicate that these conditions may have overlapping genetic causes. Less is known about the degree to which specific autistic traits relate to specific behaviours characteristic of ADHD. We hence tested, using the classical twin design, whether specific dimensional autistic traits, including social difficulties, communication atypicalities and repetitive behaviours, would display differential degrees of aetiological overlap with specific traits of ADHD, including hyperactivity/impulsivity and inattention. Parents of approximately 4,000 pairs of 12-year-old twins completed the Childhood Autism Spectrum Test and Conners' Parent Rating Scale. These measures were divided into subscales corresponding to different types of autistic and ADHD behaviours. Twin model fitting suggested that the degree of genetic overlap was particularly strong between communication difficulties and traits of ADHD (genetic correlations = .47-.51), while repetitive behaviours and social difficulties showed moderate (genetic correlations = .12-.33) and modest (.05-.11) genetic overlap respectively. Environmental overlap was low across all subscales (correlations = .01-.23). These patterns were also apparent at the extremes of the general population, with communication difficulties showing the highest genetic overlap with traits of ADHD. These findings indicate that molecular genetic studies seeking to uncover the shared genetic basis of ASC and ADHD would benefit from taking a symptom-specific approach. Furthermore, they could also help to explain why studies of the communication abilities of individuals with ASC and ADHD have produced overlapping findings.

Introduction

Based on their diagnostic criteria, autism spectrum conditions (ASC) and attention-deficit/hyperactivity disorder (ADHD) have little in common. ASC are characterised by atypical social interaction, communication difficulties, and behavioural inflexibility, while ADHD is defined by excessive hyperactivity, impulsivity, and/or inattentiveness [1]. Yet it is common for these conditions to co-occur; ADHD is thought to present in 30–60 % of individuals with ASC [2, 3], while autistic traits appear elevated in individuals with ADHD [4]. It is thus necessary for research to explore why these conditions so often co-occur with one another.

ASC and ADHD are highly heritable [5, 6]. One possible explanation for their co-occurrence could be that their genetic and environmental causes overlap with one another. This assertion has been supported by a multitude of recent twin studies. These studies have reported a moderate degree of genetic overlap between ASC and ADHD, both as dimensional characteristic traits and diagnosed conditions, in middle childhood [7], early adolescence [8-10], early adulthood [11], and adulthood [10]. These studies also report a modest [7–9, 11] to moderate [10] degree of environmental overlap across these traits.

There is converging evidence from family and twin studies that the core autistic traits, social and communication difficulties and restricted repetitive behaviours, have largely different causes [12-16]; that is, certain genes or environments may be associated with social difficulties, for example, but may not cause communication difficulties or rigid behavioural patterns. There is also evidence of some genetic specificity across traits of ADHD [17-19]. As such, it is possible that the degree of aetiological overlap between autistic traits and ADHD behaviours differs across particular types of these traits. Two recent twin studies tested this possibility. A Dutch study reported that attention switching difficulties, rather than social autistic traits, showed particularly strong genetic overlap with traits of ADHD [20]. A study of a Swedish cohort found that repetitive routines and behaviours showed elevated genetic overlap with traits of ADHD [21]. Together, these studies indicate that traits of ADHD are particularly associated with non-social autistic traits.

There have yet to be any twin studies of these associations in extreme-scoring groups. Individuals displaying extreme scores on trait measures typically display scores that are more comparable to clinical groups, thus indicating whether or not we would expect trait findings to extend to clinical groups. For this reason, we conducted a twin study aimed at testing the degree of genetic overlap between specific autistic traits, including social difficulties, communication atypicalities, and repetitive routines and behaviours, and specific traits of ADHD, including hyperactivity/impulsivity and inattention, at the extremes of the general population. We also further tested genetic and environmental overlap between these traits across the general population. We hypothesised that the degree of aetiological overlap between autistic traits and traits of ADHD would differ according to different types of autistic and ADHD trait. We further expected repetitive routines and behaviours to show stronger aetiological overlap with traits of ADHD than social and communicative autistic traits.

Method

Participants

Parents of 5,879 twins participating in the Twins Early Development Study (TEDS) completed and returned questionnaires when twins were aged 12. TEDS is a community-based, population representative study of twins born in England and Wales between 1994 and 1996 [22]. Exclusions were conducted for genetic syndromes, including Fragile X and cystic fibrosis, chromosomal abnormalities, including Down Syndrome and cerebral palsy, extreme perinatal or prenatal difficulties, and missing first contact or zygosity data. Participants with confirmed ASC (N = 80) were included. The final sample comprised 5,689 twin pairs, including 927 monozygotic (MZ) male pairs, 1,124 MZ female pairs, 865 dizygotic (DZ) male pairs, 989 DZ female pairs, and 1,784 DZ opposite-sex pairs. Zygosity was ascertained through DNA testing and parent ratings of twin resemblance [23]. All participants provided written informed consent prior to participation. TEDS has full ethical approval from the King’s College London Research Ethics Committee.

Measures

Autistic traits

Parents completed the Childhood Autism Spectrum Test (CAST [24]), a 30-item questionnaire that enquires about various autistic traits, such as lack of social interest, pronounced interests, and difficulties with conversation. Each question was answered ‘yes’ or ‘no’, meaning the maximum possible score was 30. Scores over 15 predict ASC diagnoses with 100 % sensitivity and 97 % specificity [25]. The measure was divided into three scales based on DSM-IV criteria, which were employed in previous studies [13, 14]: social difficulties (‘Social’), communication atypicalities (‘Communication’), and repetitive, restricted behaviours and interests (RRBI). The number of items, possible range of scores, and Cronbach’s α values for each subscale are shown in Table 1.

Table 1

Descriptive statistics, phenotypic correlations, and twin correlations

Measure	Number of items	Cronbach’s α	Skewness	x¯ Full sample(95 % CI)	x¯ 5 % extremes(95 % CI)	Probands (N)
CAST social	11	0.54	1.42 (−0.01)	1.61 (1.58/1.64)	5.91 (5.80/6.01)	546
CAST communication	12	0.64	1.36 (0.02)	1.90 (1.87/1.93)	7.01 (6.91/7.12)	590
CAST RRBI	7	0.49	.96 (−0.05)	1.41 (1.39/1.44)	4.47 (4.42/4.57)	747
Conners ADHD hyperactivity/impulsivity	9	0.84	1.70 (−0.09)	4.30 (4.22/4.38)	16.85 (16.57/17.14)	656
Conners ADHD inattention	9	0.90	1.38 (−0.32)	5.58 (5.48/5.67)	19.89 (19.63/20.14)	632

Skewness values are first given for untransformed scales; the values in parentheses are for transformed, standardised residuals used in the analyses

CAST Childhood Autism Spectrum Test, CAST social CAST social difficulties subscale, CAST communication CAST communication atypicalities subscale; CAST RRBI CAST repetitive, restricted behaviours and interests subscale

Traits of ADHD

Parents completed the ADHD subscale of the Conners’ Parent Rating Scale (‘Conners ADHD’ [26] ), a measure formed of 18 statements. Each statement described a behaviour characteristic of ADHD; parents rated, on a three-point scale, how true each statement was of their child. The highest possible score was 54. Conners ADHD is a valid assessment of ADHD traits; individuals with ADHD display significantly higher scores on the measure than controls [26]. Consistent with previous recommendations [26], the measure was divided into two subscales, which closely corresponded to DSM-IV criteria for ADHD: Hyperactivity/Impulsivity and Inattention. Possible ranges of scores, number of items, and Cronbach’s α for both subscales are provided in Table 1.

Data analysis

Twin analyses of the full sample

Data were prepared for analysis by first log transforming any skewed subscales (see Table 1), and then regressing the effects of sex and age out of all measures, as is standard behavioural genetic procedure [27]. Analyses were then performed on standardised residual scores using Mx [28].

First, phenotypic associations between the CAST and Conners ADHD subscales were established using phenotypic correlations (rph), correlation coefficients estimated from structural equation models. Second, twin correlations, which are the foundation of twin analysis, were estimated. Cross-trait cross-twin correlations provide one with an initial indication of the extent to which covariance between two traits is influenced by genetic and environmental factors. These involved correlating one twin’s score on a CAST subscale with their co-twin’s score on a Conners ADHD subscale. Cross-trait cross-twin correlations were estimated separately for MZ and DZ twins. These estimates cannot exceed rph between two traits. As MZ twins are assumed to share all of their segregating DNA code with one another, while DZ twins share, on average, 50 %, additive genetic (A) influences on the covariance between two traits are implicated if the MZ cross-trait cross-twin correlation exceeds the DZ estimate. As MZ twins are assumed to be genetically identical, any within pair MZ differences are assumed to be caused by nonshared environment (E), environmental factors that differ across twins in a pair and create differences between them (measurement error is included in this term). Hence, if the MZ cross-trait cross-twin correlation between two traits is less than rph between them, E is implicated on their covariance. Shared environmental (C) influences, on the other hand, are common to both twins in a pair and heighten their similarity; these influences are indicated if the DZ cross-trait cross-twin correlation exceeds half the MZ statistic. Alternatively, non-additive genetic (D) influences are implicated if the DZ cross-trait cross-twin correlation is less than half the MZ estimate. These correlations were estimated from a constrained saturated model (see below).

Finally, multivariate structural equation twin models were fitted to data to formally estimate A, C, D, and E. Cholesky decompositions, presented here as correlated factors solutions [29], were tested. For each trait, A, C (or D, if this is suggested by the twin correlations), and E were estimated. The additive genetic correlation (rA) then estimates A overlap between traits. This estimate falls between 0 and 1; estimates of 1 suggest total genetic overlap across two traits, while estimates of 0 would suggest that the genetic influences on two traits are totally independent of one another. Shared environmental (rC), nonshared environmental (rE), and non-additive genetic (rD) correlations can also be computed, and operate in the same manner. The model then computes bivariate heritability, which estimates the extent to which A can explain rph between two traits, from the univariate A estimates and rA. Equivalent statistics are calculated for C, E, and D.

ACE and ADE models were tested. Parameters in these models can either be equated across sexes, or allowed to differ in magnitude across sexes (quantitative sex differences). Such differences are indicated when male and female twin correlations differ. In order to test for quantitative sex differences, a saturated model that allowed twin correlations to differ by sex was fitted. These estimates were then equated across sexes.

The fit of each model was compared against that of a saturated model of the observed data. For each model, the −2LL fit statistic was computed. Model fit was assessed using the likelihood-ratio test (LRT), which capitalises on the fact that differences in −2LL between two models are χ2 distributed, with degrees of freedom equal to the difference in number of estimated parameters. Significant χ2 results suggest that a given model is a significantly poorer fit than the saturated model. Model fit was further assessed using Akaike’s Information Criteria (AIC) and Bayesian Information Criteria (BIC). Lower, preferably negative, AIC values indicate good model fit, while models with more negative BIC estimates are to be favoured. Selection of the best-fitting model was based on BIC, which performs better than either the LRT or AIC in larger samples [30]; when comparing two models, the model with the more negative BIC estimate is regarded as better-fitting, with differences of 10 or more argued to show that a given model is a good fit [31].

Analysis of the extremes

DeFries–Fulker analysis [32, 33] was used to analyse extreme scores on the CAST and Conners ADHD subscales. DeFries–Fulker analysis is a regression analysis of means. The univariate procedure is designed to estimate group heritability (), which refers to genetic influences on extreme scores [30]. Bivariate DeFries–Fulker analysis is concerned with testing for bivariate heritability (h2.xy); in these analyses, this refers to the extent to which the genes that cause extreme scores on one trait influence continuous scores on another trait [32]. In performing univariate and bivariate DeFries–Fulker analyses, it is possible to estimate genetic correlations between extreme scores on two measures [33].

For the univariate DeFries–Fulker analyses, extreme scorers, termed probands, were selected from the sample used in the twin analyses on the basis of scoring within the highest 5 % of the z-score distribution on a measure of interest. Scores were then transformed so that the mean proband score on the measure of interest was 1, while the control mean was 0. In transforming scores in this manner, one can gain an indication of by examining the mean scores of co-twins of probands. If the DZ co-twin mean is closer to 0 than the MZ co-twin mean, then greater resemblance between MZ twins than DZ twins is suggested. Much like how one interprets twin correlations, this can be taken as evidence of on a trait. Univariate DeFries-Fulker analysis more formally estimates through a regression equation for predicting co-twin scores on a given measure from proband scores and zygosity. The regression coefficient of co-twin scores on zygosity is an estimate of . The estimate of should not exceed the mean scores of MZ co-twins, unless non-additive genetic effects are influential. In this instance, estimates were fixed to equal the mean MZ co-twin score.

In bivariate DeFries–Fulker analysis, probands were selected on the basis of scoring within the highest 5 % of the z-score distribution on a given measure; the measure used to select probands is the selection variable. The outcome variable is the second measure of interest. Thus, the bivariate procedure is directional in that a CAST subscale could be used as the selection variable and a Conners ADHD subscale used as the outcome variable, and vice versa.

Prior to conducting bivariate DeFries–Fulker analysis, the phenotypic associations between the selection and outcome variables were tested using phenotypic group correlations. These were calculated by dividing proband’s mean z-score on the selection variable by their mean z-score on the outcome variable, and can be interpreted in the same manner as correlation coefficients. Subsequently, data were prepared for bivariate analysis by dividing scores on both measures by the zygosity-specific proband mean for the selection variable [34]. Transforming scores in this manner meant that the mean score on the selection variable for probands was 1, while the mean score for controls on the outcome measure was 0. The co-twin means provide an indication of the extent of h2.xy between two measures; if the DZ co-twin mean on the outcome variable is closer to 0 than MZ co-twin means, h2.xy between the selection and outcome variables is indicated. The bivariate DeFries–Fulker regression equation, which predicts co-twin scores on the outcome measure from proband scores on the selection measure and zygosity, was then used to formally estimate h2.xy. The regression coefficient of co-twin scores on zygosity is an estimate of h2.xy. Again, this estimate only exceeds the mean scores on the outcome variable for MZ co-twins if non-additive genetic influences are in operation; as in the univariate analysis, h2.xy was fixed to equal the MZ co-twin mean for the outcome variable in these instances. These analyses were conducted using the CAST subscales as selection variables, and then using the Conners ADHD subscales as selection variables.

The results of the univariate and bivariate DeFries–Fulker analyses were then used to calculate genetic correlations between extreme CAST subscale scores and extreme Conners ADHD subscale scores, using the equation: [35]. β and β represent h2.xy estimates, while β and β are estimates for both measures.

Results

Full sample

Descriptive statistics for all measure subscales are shown in Table 1. Table 2 shows rph estimates. The highest estimate was between Communication and Hyperactivity/Impulsivity (.43), followed by Communication and Inattention (.40). RRBI displayed modest rph with Hyperactivity/Impulsivity (.26) and Inattention (.05), while associations were weakest between Social and hyperactivity/impulsivity (.06) and Inattention (.07).

Table 2

Phenotypic correlations

Conners ADHDCAST subscalesubscale	Conners ADHD hyperactivity/impulsivity	Conners ADHD inattention
CAST social	0.06*	0.07*
CAST communication	0.43*	0.40*
CAST RRBI	0.26*	0.05*

CAST Childhood Autism Spectrum Test, CAST social CAST social difficulties subscale, CAST communication CAST communication atypicalities subscale; CAST RRBI CAST repetitive, restricted behaviours and interests subscale, Conners ADHD Conners’ Parent Rating Scale, Conners ADHD hyperactivity/impulsivity Conners hyperactivity/impulsivity subscale, Conners Inattention Conners inattention subscale

* p < 0 .05

Twin correlations are shown in Table 3. These suggest that each CAST and Conners ADHD subscale was heritable. Cross-trait cross-twin correlations suggested that the covariance between the CAST and Conners ADHD subscales was largely influenced by genetic factors; with the exception of the cross-trait cross-twin correlations between RRBI and Inattention, MZ estimates exceeded DZ estimates. Many of the DZ estimates were also less than half the MZ estimates, indicating that D may have played a role in the covariance between some of the subscales. No MZ cross-trait cross-twin correction exceeded rph, indicating E.

Table 3

Twin correlations

Cross-twin correlations
Measure	Monozygotic twins		Dizygotic twins
	Estimate	95 % CI	Estimate	95 % CI
CAST Social	0.72	0.70–0.74	0.34	0.30–0.38
CAST Communication	0.76	0.74–0.77	0.39	0.35–0.43
CAST RRBI	0.71	0.69–0.73	0.30	0.26–0.34
Conners Hyp/Imp	0.88	0.87–0.89	0.50	0.46–0.53
Conners Inattention	0.78	0.76–0.79	0.37	0.33–0.41

CAST Childhood Autism Spectrum Tes, CAST social CAST social difficulties subscale, CAST communication/CAST comm. CAST communication atypicalities subscale, CAST RRBI CAST repetitive, restricted behaviours and interests subscale, Conners hyp/imp Conners’ Parent Rating Scale (Conners ADHD) hyperactivity/impulsivity subscale, Conners inattention/Conners inatt. Conners ADHD inattention subscale

Fit statistics for the twin models tested are given in Table 4. A saturated model that equated twin correlations across sexes fitted better than one that allowed them to differ by sex, as indicated by a considerably lower BIC estimate for the constrained saturated model (−157,379.80) than for the saturated model with sex differences (−156,485.70). ACE and ADE models were thus tested with equivalent estimates in both sexes. Models included all three CAST and both Conners ADHD subscales. Relative to the saturated model, an ADE correlated factors solution had most negative BIC value (−157,928.84). In dropping D from the model, BIC (−158,004.70) reduced further, hence an AE model was chosen as the best fitting model. Parameter estimates from the AE model are all shown in Table 5. All three CAST subscales were heritable, with A explaining 69–72 % of the variance in each. The same was true for both Conners ADHD subscales (A = 75–86 %). The aetiological correlations between the measure subscales suggested that the strongest degree of genetic overlap was between Communication and both Hyperactivity/Impulsivity (rA = .50) and Inattention (rA = .48). Genetic correlations between the CAST communication difficulties subscale and traits of ADHD were significantly stronger than all other genetic correlations; there was moderate A overlap between RRBI and Hyperactivity/Impulsivity (rA = .31). However, rA was low between all other subscales (rA = .08–.14). Nonshared environmental overlap between the three CAST subscales and Conners ADHD subscales was low. The highest rE estimates were between Communication and Hyperactivity/Impulsivity (rE = 22) and Inattention (rE = .19). This was again significantly higher than the estimates seen between other subscales. RRBI displayed greater E overlap with hyperactivity/impulsivity (rE = .11) than Inattention (rE = −.03), while rE between Social and Conners ADHD subscales was low and largely non-significant.

Table 4

Twin model fit statistics

Model	−2LL	df	Parameters	BIC	Δχ 2	Δdf	p	AIC
Saturated	66,869.83	27,116	130	−157,379.84	–	–	–	–
ACE	67,001.24	27,196	50	−157,910.04	131.41	80	<0.001	−28.59
ADE	66,982.44	27,196	50	−157,928.84	112.61	80	<0.05	−47.39

All submodels are nested within the ADE model, which was a better fit relative to the saturated model than the ACE model

−2LL fit statistic, BIC Bayesian Information Criteria, AIC Akaike’s Information Criteria, Δχ change in −2LL between two models, Chi square distributed

aIndicates best fitting model; p values are corrected for multiple comparisons using the Bonferoni correction

Table 5

Estimates derived from the best fitting correlated factors solution

	A		E
	Estimate	95 % CI	Estimate	95 % CI
Variance Components Estimates
Social	0.72	0.68–0.75	0.28	0.27–0.30
Comm.	0.71	0.69–0.73	0.29	0.27–0.30
RRBI	0.69	0.66–0.70	0.31	0.30–0.33
H/I	0.86	0.84–0.87	0.14	0.12–0.15
Inatt.	0.75	0.74–0.76	0.25	0.22–0.26

All italicised estimates were not significant, as indicated by confidence intervals overlapping with zero

CAST Childhood Autism Spectrum Test, Social social difficulties subscale of the CAST, Communication/Comm. communication atypicalities subscale of the CAST, RRBI repetitive, restricted behaviours and interests subscale of the CAST, Conners ADHD Conners’ Parent Rating Scale, Hyp/imp hyperactivity/impulsivity subscale of conners adhd,Inatt inattention subscale of Conners ADHD, A additive genetic influences,E nonshared environmental influences, r additive genetic correlation, r nonshared environmental correlation

The majority of rph between the CAST and Conners ADHD subscales could be explained by A (see Table 5). Bivariate heritability was estimated as at least 80 % got each rph estimate, with estimates of bivariate heritability peaking at 100 %. E, by contrast, explained no greater than 20 % of each rph estimate. These results did not significantly change when individuals with ASC were removed from the analyses (contact the first author for further details).

Extremes analysis

Results of the DeFries–Fulker analysis are provided in full in the supplementary materials. Descriptive statistics for the CAST and Conners ADHD subscales in the 5 % extreme groups, along with the number of probands in each group, are given in Table 1.

Table S1 shows the phenotypic group correlations, which broadly mirrored the pattern of rph estimates in the full sample. Communication displayed stronger associations with both Conners ADHD subscales (53–60) than RRBI (24–40) and Social (24–34). In the univariate analyses, DZ co-twin means were closer to 0 than MZ co-twin means, indicating for each subscale (see Table S2). Regression analyses confirmed this; estimates, which are shown in Table S3, ranged from .69–.73 for the CAST subscales, and from .70–.85 for the Conners ADHD subscales. In the bivariate analyses, all DZ co-twin means on the outcome measure regressed toward 0 to a greater extent than MZ co-twin means, suggesting h2.xy (see Table S2). Further, bivariate heritability estimates, given in Table S3, indicated that phenotypic group correlations between the CAST and Conners ADHD subscales were almost entirely explained by genetic factors.

From the and h2.xy estimates derived from the DeFries-Fulker analyses, genetic correlations between extreme scores on the CAST and Conners ADHD subscales were calculated. Generally, these were higher than rA estimates derived in the full sample. The highest estimate was between communication and inattention (.65); communication also displayed considerable genetic overlap with hyperactivity/impulsivity (.60). RRBI displayed moderate genetic overlap with both Conners ADHD subscales (.35–.36). Genetic overlap between Social and the two Conners ADHD subscales was weak (.21) to moderate (.45).

Discussion

We tested the degree of genetic and environmental overlap between specific autistic traits and specific traits of ADHD in the general population, including, for the first time, at the extremes of the general population. The strongest genetic associations, in both the full sample and at the extremes, were between communication difficulties characteristic of ASC and traits of ADHD, with little differentiation in the magnitude of this association across different traits of ADHD. There was also a modest degree of phenotypic and genetic overlap between repetitive routines and interests, and traits of ADHD, and least overlap was observed for autistic social difficulties. As expected, the degree of aetiological overlap between autistic traits and traits of ADHD varied across different types of autistic trait.

In contrast with our expectations, and with the two existing twin studies on this subject [20, 21], we did not find that non-social autistic traits showed the strongest association between autistic and ADHD behaviours. It is notable, however, that we did still find, in some instances, moderate genetic overlap between repetitive routines and interests and traits of ADHD, at the extremes and with hyperactivity in the full sample. Hence, our findings do support the notion that traits of ADHD do, at least, show stronger associations with non-social than social autistic traits. The key question is why communication difficulties, specifically, showed the strongest overlap with traits of ADHD.

Item overlap between the CAST Communication subscale and the Conners ADHD subscales represents one possible explanation for why traits of ADHD were so strongly linked with autistic communication difficulties relative to other autistic traits. The CAST Communication subscale is concerned with pragmatic aspects of communication, enquiring specifically about the manner in which the child interacts with other people. The Conners ADHD also asks about these aspects of communication, such as whether the child often struggles to maintain conversations due to their ADHD behaviours. Yet there is considerable evidence that these measures are valid assessments of traits of ASC [25] and ADHD [26], indicating that there can be confidence that these results represent genuine genetic overlap between communication difficulties and traits of ADHD.

From these findings one might predict that individuals with ADHD will be more likely to display communication difficulties similar to those seen in ASC. Indeed, research on the language abilities of individuals with ASC and ADHD are producing somewhat overlapping findings. While structural language skills show considerable variability in both these conditions, pragmatic language difficulties appear to feature prominently in both [36, 37]. This pattern might be expected if it is the case that the genes that are associated with communication difficulties characteristic of ASC also link with traits of ADHD.

On a theoretical level, these findings lend additional support to the fractionable autism triad hypothesis [38, 39]. This hypothesis posits that the core features of ASC are each associated with largely different causal influences, and is supported by twin studies [12-17]. These findings extend this hypothesis, by suggesting that the manner in which core autistic traits relate to traits characteristic of ADHD differs across autistic trait domains. Intriguingly, these findings are similar to those reported in a twin study, also conducted on the TEDS sample, that suggested heightened aetiological overlap between communication difficulties and internalising traits [40], and also reflect those of a recent longitudinal twin study of the phenotypic associations between ASC and ADHD [8]. Thus, this study adds to increasing evidence to suggest that the fractionable autism triad hypothesis extends to how autistic traits relate to traits of other conditions.

Findings such as these have the potential to inform molecular genetic research endeavours targeted at identifying genetic variants that are shared across ASC and other conditions. This is clearly a distal goal; most molecular genetic studies have focused on ASC and ADHD as single constructs. For instance, one study reported that genetic variants, specifically calcium-signalling genes, were associated with both ASC and ADHD [41]. Could these variants be associated, in particular, with communicative aspects of ASC? Research taking a symptom-specific approach to ASC is scarce [e.g. 42]. Encouragingly, a recent paper, using data from the ALSPAC sample, tested the association between polygenic scores for ADHD and various neurodevelopmental difficulties, including autistic traits. Of note, the authors reported that while polygenic scores for ADHD did not significantly associate with social cognition difficulties, they were particularly linked with pragmatic language abilities [43], which is a similar pattern to that seen in the present study. Thus, our findings, along with these independent results, suggest that this is a potentially fruitful future direction for molecular genetic research on ASC and ADHD overlap.

These findings need to be regarded in light of the limitations of this study. Owing to the large sample size, requisite for twin model fitting to be adequately performed and for sufficient probands to be identified for extremes analysis, we were unable to conduct in-depth assessments on participants. However, this need not be viewed as a limitation. There is evidence that the heritability of traits of ASC and ADHD is constant across the general population and at the extremes [5, 6], and we did find similar, albeit slightly stronger, associations at the extremes of our sample. Furthermore, we were able to avoid referral biases that can be present in clinically ascertained samples, which can inflate co-occurrence rates between conditions. One might question whether these findings from twins generalise to singletons. However, autistic traits and traits of ADHD do not appear to be elevated in twins relative to singletons [44, 45]. Future studies should also focus on assessments of ASC and ADHD from multiple raters, rather than parents alone. The internal consistency of the three CAST subscales was relatively low, as noted previously [13, 14]. Particular items failing to contribute to the internal consistency sufficiently did not account for this, since removing any items did not result in discernable improvements in internal consistency. The modest internal consistency is, perhaps, due to the relatively small number of items that comprise each of the subscales, or heterogeneity within each subscale. An important future direction, then, is to assess autistic traits in a manner that maximises reliability, although good construct validity for the CAST subscales has been reported previously [13, 14].

This study demonstrated that the degree to which autistic traits and traits of ADHD share genetic and environmental influences with one another varies by autistic trait domain, showing for the first time that this pattern is also apparent amongst individuals in the general population displaying particularly extreme degrees of traits of ASC and ADHD. Autistic communication difficulties, in particular, were associated with similar genetic influences to traits of ADHD, more so than autistic social or non-social behaviours. These findings from a community sample suggest the importance of taking a trait-specific approach to understanding why ASC co-occurs with other conditions, rather than focusing on the ASC phenotype as a whole. These findings may lead to the prediction that communication difficulties characteristic of ASC will present with increased frequency in individuals with ADHD. Indeed, similar communication difficulties, such as pragmatic difficulties, have been reported in both ASC and ADHD [36, 37]; our findings offer a potential explanation with regard to these findings by suggesting that similar causes are associated with both traits of ADHD and communication difficulties characteristic of ASC.

Supplementary material 1 (DOCX 55 kb)