Suk-Ling Ma1†
Lu Hua Chen2,3†
Chi-Chiu Lee4
Kelly Y. C. Lai1Se-Fong Hung1Chun-Pan Tang4Ting-Pong Ho5Caroline Shea6Flora Mo6
Timothy S. H. Mak2
Pak-Chung Sham2,5*‡
Patrick W. L. Leung3*‡- 1Department of Psychiatry, The Chinese University of Hong Kong, Hong Kong, China
- 2Centre for PanorOmic Sciences - Genomics and Bioinformatics Cores, The University of Hong Kong, Hong Kong, China
- 3Department of Psychology, The Chinese University of Hong Kong, Hong Kong, China
- 4Kwai Chung Hospital, Hospital Authority, Hong Kong, China
- 5Department of Psychiatry, The University of Hong Kong, Hong Kong, China
- 6Alice Ho Miu Ling Nethersole Hospital, Hospital Authority, Hong Kong, China
Background: Recent findings indicated a high comorbidity between attention-deficit/hyperactivity disorder (ADHD) and autism spectrum disorder (ASD), as well as shared genetic influences on them. The latter might contribute at least partly to the former clinical scenario. This study aimed at investigating whether SHANK genes were potential pleiotropic genes to the two said disorders, underlying their genetic overlap.
Methods: This study recruited 298 boys with ADHD (including 256 family trios of 1 ADHD boy and his 2 biological parents), 134 boys with ASD, 109 boys with both ADHD and ASD, and 232 typically developing boys as community controls. They were aged between 6 and 11 years old.
Results: There was no significant difference in allele frequency of a number of single nucleotide polymorphisms (SNPs) in SHANK2/SHANK3 between the three clinical groups (ADHD, ASD, and ADHD + ASD) and between the two control groups (community controls and pseudo-controls), respectively. The three clinical groups and the two control groups were thus, respectively, combined. A comparison between the two aggregated samples identified significant evidence of disease association for three SHANK2 SNPs with both ADHD and ASD, even after multiple testing correction: rs11236616 (OR = 0.762, permuted p = 0.0376), rs7106631 (OR = 0.720, permuted p = 0.0034), and rs9888288 (OR = 0.770, permuted p = 0.0407). Comparisons among individual groups pointed to a similar trend of findings.
Conclusion: SHANK2 could be considered a potential pleiotropic gene underlying the genetic overlap between ADHD and ASD. This might contribute partly to their high comorbidity in the afflicted children.
Introduction
Attention-deficit/hyperactivity disorder (ADHD) and autism spectrum disorder (ASD) are two separate common neurodevelopmental disorders in children. However, their diagnoses were previously hierarchical such that a diagnosis of ASD precluded a diagnosis of ADHD. The abandonment of this hierarchical diagnostic arrangement in the current psychiatric classification, Diagnostic and Statistical Manual of Mental Disorders-Fifth Edition (DSM-5) (American Psychiatric and Association, 2013) means that the same individual can receive both diagnoses. Subsequent studies have reported high comorbidity rates between the two disorders. For example, 20–50% of ADHD children met diagnostic criteria for ASD, while 30–80% of ASD children also met diagnostic criteria for ADHD (Simonoff et al., 2008; Sinzig et al., 2009). This high comorbidity cannot be explained by diagnostic or phenotypic overlap. The two disorders were found conceptually distinctive and empirically separable by a two-factor common pathway model, while the elevated ASD traits in ADHD children could not be explained by ADHD or other associated behavioral symptoms (Grzadzinski et al., 2011; Ronald et al., 2014). Instead, the high comorbidity of ADHD and ASD is likely to be suggestive of potential sharing of underlying etiologies, including genetic predisposition.
Family-based and twin studies have consistently found that both ADHD and ASD are highly heritable, each with genetic influences accounting for approximately 70–80% of the phenotypic variance and a significant portion of them being shared between the two disorders. The latter was implicated by findings of co-aggregation of ADHD and ASD in the same individuals or in the same families (Reiersen et al., 2008; Mulligan et al., 2009; Ghirardi et al., 2018; Miller et al., 2019). Recently, molecular genetic approaches, such as genome-wide association studies (GWAS), had begun to identify single nucleotide polymorphisms (SNPs) or copy number variants (CNVs) statistically associated with both ASD and ADHD (Rommelse et al., 2010; Martin et al., 2014; Zhao and Nyholt, 2017; Cross-Disorder Group of the Psychiatric Genomics Consortium, 2019). Those genes associated with these implicated SNPs/CNVs are possible pleiotropic genes which contribute to the genetic overlap between ASD and ADHD. Yet, most of the existing studies on ADHD and ASD were conducted at North American and European sites with a relatively small number of ethnic Chinese participants by proportion (Li et al., 2014; Demontis et al., 2019). This present study will be conducted with ethnic Chinese children in Hong Kong.
One strategy for identifying pleiotropic genes for ADHD and ASD with a modest sample size is to take a candidate gene approach, e.g., starting with genes found to be associated with one disorder and testing their association with the other. SHANK genes are good candidates for this strategy. They had established association with ASD, reported in both meta-analysis (Leblond et al., 2014) and other independent studies with ASD cohorts and knock-out mice (Peça et al., 2011; Leblond et al., 2014; Connolly et al., 2017). In one animal study, treatment of SHANK2–/– knock-out mice through restoring NMDA receptor function improved the mice’s ASD traits of impaired social interaction (Won et al., 2012). SHANK gene family consists of three members (SHANK1, SHANK2, and SHANK3), which encodes scaffolding proteins required for the proper formation and functioning of neuronal synapses (Monteiro and Feng, 2017). This protein family interacts with a number of synaptic proteins in the postsynaptic density (PSD) of glutamate receptors (NMDA receptors and AMPA receptors). Specifically, SH3 and multiple ankyrin repeat domains protein 2 encoded by SHANK2, and SH3 and multiple ankyrin repeat domains protein 3 encoded by SHANK3 are the two proteins first expressed at the developing PSD. These proteins have been consistently identified to be important for cognitive and emotional well-being, including social behavior (Guilmatre et al., 2014). Thus, there are good theoretical grounds for SHANK genes to be associated with ASD.
Though there are fewer reported association between SHANK genes and ADHD, our proposition is not entirely without basis. First, paternal balanced translocation of SHANK3 was reported to result in both ASD and ADHD (Durand et al., 2007; Moessner et al., 2007). Second, SHANK2–/– knock-out mice displayed both hyperactive and autistic-like behavior (Schmeisser et al., 2012). Based on these lines of existing evidence, the present study aims at testing SHANK2 and SHANK3 for pleiotropic effects by examining their common variants (SNPs) for association with both ASD and ADHD.
Materials and Methods
Participants
The gender ratios of children with ADHD and children with ASD in the clinic population of Hong Kong were around 10 (male) to 1 (female), as reported in a local study (Yuen, 2017) and in an European study (Nøvik et al., 2006). The marked male preponderance was probably due to detection and referral biases (i.e., fewer girls being detected in the communities and/or referred to clinics), which likely exaggerated the already imbalanced community gender ratios in the eventual clinic population (Rutherford et al., 2016; Young et al., 2020). Thus, only boys were included in this study. Otherwise, having a small number of girls in the sample would make the interpretation of the combined boy/girl results difficult, while separate analysis for gender would prove untenable unless a very huge sample was planned for. Participants consisted of a clinical sample and a community sample. The clinical sample included three groups of boys aged 6–11 years with a primary clinical diagnosis of either ADHD (combined type), ASD or both according to DSM-5: (1) boys with a primary clinical diagnosis of ADHD (combined type) without ASD (ADHD group), (2) boys with a primary clinical diagnosis of ASD without ADHD (ASD group), and (3) boys with dual clinical diagnoses of ADHD and ASD (ADHD + ASD group). Exclusion criteria included the clinical diagnosis of psychosis, bipolar disorder, intellectual disability, Tourette’s or persistent motor or vocal tic disorder, or major neurological or medical illnesses. Operationally, with respect to intellectual disability, children with IQ < 80 were excluded, given potential measurement errors. A community control group of typically developing boys was also recruited. Exclusion criteria, in addition to those for the clinical sample, were parental reports of signs of ADHD, ASD, or other specific learning problems. Research assistants were trained to ask for symptoms of these disorders/problems to parents without using standardized rating scales, as open reports by local parents were generally considered as credible and reliable. This was because local Chinese schooling system was highly demanding both in terms of academic attainment and discipline, the latter including social cooperation among children in group activities, so that children with ASD, ADHD and other specific learning problems were not difficult to be spotted and would be referred for formal assessment in public services. Children whose parents indicated concerns on the above-listed problems or had received queries/complaints/referral for assessment from schools or others on these problems were excluded from the present study.
Procedure
This study was approved by the Joint Chinese University of Hong Kong-New Territories East Cluster Clinical Research Ethics Committee and complied with the latest version of the Declaration of Helsinki. The clinical sample was recruited from three child psychiatric clinics in Hong Kong. The community control sample was recruited via local mainstream schools or advertisements in social media. The interested parents were first asked two sets of questions to screen whether their children, including the sibling(s) of the child to be participating, had ADHD, ASD, or other specific learning problems. Children, whose parents gave a positive reply to one of the two sets of questions, were declined participation in this study. Parents of the clinical and community control samples provided written informed consent for their children to participate in the study.
Additionally, for the ADHD group, designated funding was available to examine the DNA of the two biological parents for a family-based study of ADHD without the involvement of ASD. Written informed consent was thus obtained from these parents for their own participation. However, this family-based study of ADHD is beyond the scope of this manuscript to be reported here.
Genomic DNA Extraction
Venous blood was collected from the ADHD group and their biological fathers and mothers. Genomic DNA was isolated using QIAamp DNA Blood Mini Kit (Qiagen, Hilden, Germany) in accordance to a standard protocol. Buccal swabs were collected from the ASD and ADHD + ASD children using Catch-All Buccal Swab (Epicentre Biotechnologies, Madison, United States). Saliva was collected from the community control group using Oragene DNA Collection Kit (DNA Genoteck, Ontario, Canada). Genomic DNA from collected buccal swabs or saliva was isolated by prepIT DNA extraction kit (DNA Genoteck, Ontario, Canada) using a standard procedure.
SNPs Selection, Genotyping, and Quality Control
Tag SNP strategy was used to select the SNPs in SHANK2 and SHANK3 for genotyping. A cluster-based algorithm which clustered SNPs by r2 was used (Carlson et al., 2004). The algorithm determined the pairwise r2 between all SNPs and clustered correlated SNPs with r2 over 0.8 into bins. Each tagSNP selected represented all the SNPs in the same bin. Fifteen tagSNPs for SHANK2 and eight tagSNPs for SHANK3 were, respectively, selected by a R2 algorithm among SNPs with minor allele frequency (MAF) of at least 10% in Asians (data from HapMap project Phase I- I- I) (Tang et al., 2006). Genotyping was performed using MassARRAY iPLEX platform (Sequenom, San Diego, United States). SNPs with MAF less than 5%, missing rates higher than 10%, Hardy-Weinberg Equilibrium (HWE) p-values smaller than 0.1% or above 5% Mendel error rates were excluded. Altogether, three SNPs were excluded with one SNP (rs7103613) in SHANK2 deviating from HWE distribution (p < 0.1%), one SNP (rs13055562) in SHANK3 having a high missing rate (>10%), and one SNP (rs4980613) in SHANK2 having a >5% Mendel error rate.
Statistical Analysis
From the ADHD family trios’ data, we constructed matched pseudo-controls by an in-house R program, which was based on transmitted and non-transmitted parental alleles using the Haplotype-Relative-Risk (HRR) principle (Falk and Rubinstein, 1987). According to this principle, for each SNP locus, the pair of parental non-transmitted alleles was regarded as the genotype of the pseudo-controls. The pseudo-controls could be considered as real random controls and were added to the community typically developing controls to increase the sample size of the controls and thus the statistical power of our significance tests.
The differences in allele frequency between the clinical and control groups were assessed by the χ2 test using the PLINK software (Purcell et al., 2007). Haploview v4.2 was employed for visualizing the linkage disequilibrium (LD) pattern. Multiple testing correction was performed by the permutation test with 10,000 times. In this study, permuted <0.05 was regarded as statistically significant. Gene-based association analysis was performed on KGG software (a systemic biological knowledge-based mining system for genome-wide genetic studies), using its GATES function (Li et al., 2010).
Results
Participants Recruited
The sample sizes of the three clinical groups were, respectively, 298 boys for the ADHD group (including 256 family trios of 1 ADHD boy and his biological father and mother), 134 boys for the ASD group and 109 boys for the ADHD + ASD group. The community control group had 232 typically developing boys. The mean ages of the ADHD group, ASD group, ADHD + ASD group, and community control group were 9.3 years (SD = 1.5), 8.0 years (SD = 1.8), 8.5 years (SD = 1.4) and 9.2 years (SD = 1.7), respectively. There was no significant age difference between the ADHD and control groups, as well as between the ASD and ADHD + ASD groups. However, the former two groups were significantly older than the latter two by about 1 year (p = 0.001). Two hundred and fifty-six pseudo-controls were constructed based on 256 family trios’ genotyping data.
SNP-Level Association Analysis
Thirteen SNPs in SHANK2 and seven SNPs in SHANK3 were qualified for data analysis after quality control. Their allele frequency is listed in Table 1. A series of planned hypothesis-driven comparisons were executed. First, we tested whether there were significant differences in the allele frequency between the three clinical groups of ADHD, ASD, and ADHD + ASD, given our hypothesis on SHANK2/3 as pleiotropic genes, and between the two control groups, i.e., community controls and pseudo-controls, given the latter considered as real random controls. The results showed no significant difference in the allele frequency of the SNPs between the three clinical groups and between the two control groups (Table 1). Given this lack of significant difference among the clinical and control groups, their data were, respectively, merged to form one single aggregated clinical sample to be compared against one single aggregated control sample. This aggregation of groups had the advantage of enlarging the comparison group size and thus increasing statistical power. Assuming an OR of 1.5 and minor allele frequency of 0.4, the statistical power of the comparison between the aggregated clinical sample and the aggregated controls was 90% at a critical p-value of 0.001. Before aggregation, the estimated statistical power was, respectively, 72, 36, and 28% for ADHD, ASD, and ADHD + ASD groups vs. the aggregated controls (Genetic Power Calculator)1. Significant association with both ADHD and ASD was found for six SHANK2 SNPs, namely, rs7113016 (p = 0.0273, OR = 0.819), rs1073294 (p = 0.0075, OR = 1.275), rs11236616 (p = 0.0034, OR = 0.762), rs7106631 (p = 0.0002, OR = 0.720), rs10899158 (p = 0.0483, OR = 1.194), rs9888288 (p = 0.0037, OR = 0.770) (Table 1). Three of them survived the multiple testing correction and remained statistically significant (Table 2), i.e., rs11236616 (permuted p = 0.0376), rs7106631 (permuted p = 0.0034), and rs9888288 (permuted p = 0.0407). Significant protective effects were found from the A allele of rs11236616 as well as the T allele of both rs7106631 and rs9888288, all of which decreased disease risk of both ADHD and ASD by around 20–30%. These percentages were calculated based on ORs of those three SNPs. The ORs for rs11236616, rs7106631 and rs9888288 were, respectively, 0.762, 0.720, and 0.770 (Table 2). Assuming approximate equivalence between OR and relative risk for low-frequency diseases, an OR of 0.72 represented a reduction of 0.28 in terms of risk. Thus, ORs between 0.72 and 0.77 represented a risk reduction of around 20–30%. Further genotypic association analysis exhibited even stronger protective effects of rs11236616 (AA vs. AG + GG, OR = 0.549, permuted p = 0.0376) and rs7106631 (TT vs. TG + GG, OR = 0.515, permuted p = 0.0013) in a recessive model (Table 2). No significant result was detected for SHANK3 SNPs by allelic association analysis (Table 1).
Table 1. Allelic association analysis of SNPs in SHANK genes.
Table 2. Significant SNPs in SHANK2 (ADHD/ASD/ADHD + ASD vs. control/pseudo-control).
We then proceeded to finer comparison between groups (Table 1). Each clinical group was individually compared to the combined control sample. We expected fewer instances of significant differences, given the reduced sample size and statistical power (see above power computation). Yet, the overall trend of the results was largely identical to that of the above. Regarding the ADHD group, four of the six above-named SNPs were found to be significantly linked to ADHD with a fifth one (rs7113016) at marginal significance (p = 0.0685). The ASD and ADHD + ASD groups had even smaller sample sizes (134 and 109, respectively, as compared to 298 of the ADHD group) and thus reported fewer instances of statistical significance. Although the ASD group reported only one SNP (rs11236616) among the six SNPs identified above to be significantly associated to ASD, there were three more SNPS among them (rs7113016, rs7106631 and rs9888288) reaching marginal significance (p < 0.08). The ADHD + ASD group reported two SNPs of significant association to ADHD and ASD among the above six named SNPs (rs7106631 and rs10899158).
Finally, we created two enlarged clinical samples by diagnosis, i.e., adding ADHD and ADHD + ASD groups together to form one enlarged ADHD sample, whilst adding ASD and ADHD + ASD groups to form one enlarged ASD sample. When these two enlarged clinical samples were, respectively, compared to the aggregated control sample, the results came close to the findings involving all three clinical groups combined. For the enlarged ADHD sample, only one of the six above named SNPs (rs7113016) marginally missed the statistical significance (p = 0.0714) associated with ADHD. For the enlarged ASD sample (with a smaller sample size of 243 vs. 407 of the ADHD sample), among the above-named six SNPs, two were significantly associated with ASD (rs7106631 and rs10899158) with three more at marginal significance (rs7113016, rs1073294, and rs11236616) (p < 0.08) (Table 1).
Gene-Level Association Analysis
Gene-level association analysis integrates the p-values of a number of SNPs within a gene and provides an overall p-value for the entire gene. As shown in Table 3, at the gene-level, SHANK2, compared to the aggregated control sample, was associated with the ADHD group (p = 0.0082), the enlarged ADHD sample (i.e., involving both ADHD and ADHD + ASD groups) (p = 0.0034), and finally the aggregated clinical sample (p = 0.0031 and permuted p = 0.0061). Though not entirely unexpected, the smaller sample size and lower statistical power were unable to pick up the association between SHANK2 and ASD, represented by the ASD group or the enlarged ASD sample (involving both ASD and ADHD + ASD groups). As expected, no evidence for disease association was identified for SHANK3.
Table 3. Gene-level association analysis.
Discussion
In terms of the allele frequency of SHANK2 SNPs, our study finds no significant difference between the three clinical groups, i.e., ADHD, ASD and ADHD + ASD groups, and between the two control groups, respectively. Instead, the three clinical groups combined differ significantly (multiple testing corrected) from the two control groups combined with respect to three SHANK2 SNPs (rs7106631, rs11236616, and rs9888288). Significant protective effects are found from the A allele of rs11236616 as well as the T allele of both rs7106631 and rs9888288, all of which decrease disease risk of both ADHD and ASD by around 20–30%. Subsequent analysis with each individual clinical group or with the two enlarged clinical samples by diagnosis shows a largely similar trend of findings, despite that the actual results are impacted by the reduced sample size and statistical power of the comparison groups. To the best of the authors’ knowledge, there is no previous study reporting association of SHANK2 SNPs with both ASD and ADHD. This new finding suggests SHANK2 to be a potential pleiotropic gene to the two disorders. However, at the gene-level analysis, the comparisons between various ASD groups and the aggregated controls are non-significant (see Table 3). These results may be related to the fact that the p-value at the gene-level is an aggregate of p-values of the SNPs within a gene. The non-significant SNPs within a gene may mask the significance of a few, particularly when the statistical power is quite low.
In this study, the three strongest-linked SNPs to ADHD and ASD (rs7106631, rs11236616, and rs9888288) were all clustered in the same intronic region of the SHANK2 gene (see Figure 1). Although no direct biochemical function could be predicted for three of them, we identified 17 other SNPs which were in LD (r2 > 0.800) with them using the LD information from the Chinese subpopulation (CHB and CHS) in the 1000G project (The 1000 Genomes Project Consortium, 2015). Three of these 17 SNPs (rs10899152, rs7945377, and rs7111763) were in strong LD with rs7106631 (r2 = 0.953, 0.857, and 0.841, respectively), the latter being the one SNP identified in this study to be strongest in its association with both ADHD and ASD. Interestingly, the three SNPs associated with this SNP (rs7106631) had been classified as functionally important by RegulomeDB database, affecting transcription factor binding sites or regulating gene expression (Boyle et al., 2012). Thus, this identified SNP of our study, rs7106631, can potentially be involved in the transcriptional regulation or gene expression of SHANK2 in both ADHD and ASD. Further analysis by expression quantitative trait loci (eQTL), using FUMA2, suggested that SHANK2 rs1893121, which was in strong LD with the currently identified rs7106631 (r2 = 0.949), was significantly expressed in brain hippocampus (p = 0.046) and basal ganglia (p = 0.037). Future research should explore and define the functional involvement of our identified SNPs, particularly rs7106631, in these two disorders.
Figure 1. –log (P) graph and pair-wise LD structure of selected SNPs in SHANK2. Above graph shows association analysis results of 13 SNPs; each blue dot represents the –log (p-value) of each SHANK2 SNP. Below graph represents the pair-wise LD structure of SNPs in SHANK2; the number in each cell is the R2 between two SNPs with darker color representing stronger R2.
The considerable interrelated nature of psychiatric disorders at the genetic level is speculated to reflect in phenotypic or endophenotypic overlap between disorders, particularly the latter. For example, it would be difficult to find phenotypic overlap between ADHD and schizophrenia; yet, they shared neurocognitive deficits in the endophenotypes of executive functions (The Brainstorm Consortium, 2018). The same possibility may exist for ADHD and ASD. As noted above, their phenotypes were found to be conceptually distinctive and empirically separable (Grzadzinski et al., 2011; Ronald et al., 2014). However, the two disorders were found to share deficits in various components of executive functions, such as working memory, response inhibition, processing speed, fluency and concept formation (Craig et al., 2016; Karalunas et al., 2018). These shared neurocognitive dysfunctions are highly heritable. The involvement of SHANK genes in NMDA and AMPA glutamate receptors, which are important in defining synaptic plasticity and cognitive functioning (Guilmatre et al., 2014; Monteiro and Feng, 2017), suggests that SHANK genes may play a crucial role in memory and executive dysfunctions found in a wide range of neuropsychiatric disorders, including ADHD and ASD. Thus, one future line of investigation is to examine whether SHANK genes represent at least parts of the genetic underpinnings responsible for the development of those shared executive function deficits in ADHD and ASD.
On the contrary to SHANK2, we do not detect any disease associations with SHANK3, whose LD pattern is shown in Figure 2. Previous studies on SHANK3 with ASD have mainly been conducted with its rare variants (Durand et al., 2007; Moessner et al., 2007; Waga et al., 2011; Boyle et al., 2012; Boccuto et al., 2013; Cochoy et al., 2015; Nemirovsky et al., 2015). Regarding investigation with common variants as in this study, Connolly et al. (2017) explored two ASD family GWAS datasets and found significant maternal genetic effect for SHANK3 rs5770820 only in the Simons Simplex Consortium dataset, but not in the Autism Genome Project dataset. Furthermore, disease association with SHANK3 SNPs was also not found or replicated in one Chinese study, which also investigated rs5770820 (Qin et al., 2009). The discrepant findings may be attributable to different study methodology and focus. For example, different from traditional association analysis, the maternal genetic effect of SHANK3 on ASD detected in Connolly et al.’s study (2017) specifically investigated intrauterine phenotypic development.
Figure 2. Pair-wise LD structure of selected SNPs in SHANK3. The graph represents the pair-wise LD structure of SNPs in SHANK3; the number in each cell is the R2 between two SNPs with darker color representing stronger R2.
In contrast to ASD, SHANK2 has not been reported to be associated to ADHD in GWAS. In fact, ADHD GWAS have not been very successful in detecting a sizable number of significant SNPs in ADHD until recently by Demontis et al. (2019). This ADHD GWAS in its publication, including supplementary results, reported 12 genome-wide significant loci, but none of them was related to SHANK2. However, its sample sizes of 20,183 ADHD participants and 35,191 controls are not considered to be particularly large by comparison to those involved in other common psychiatric disorders, e.g., >450,000 participants in depression (Cross-Disorder Group of the Psychiatric Genomics Consortium, 2019). Thus, the above-cited ADHD sample size may not give sufficient power to identify a larger pool of loci associated with ADHD, which perhaps could include SHANK2. Furthermore, as mentioned above, the majority of the participants of these GWAS were from European-ancestry populations with relatively few ethnic Chinese samples by comparison (Li et al., 2014; Demontis et al., 2019).
This study identifies SHANK2 as the genetic overlap between ADHD and ASD. Yet, there must also be other factors leading to their separate development. Reiersen et al. (2008) found that maternal smoking during pregnancy (SDP) interacted with the 7-repeat allele of dopamine receptor D4 (DRD4) gene on the development of ADHD, but not on the development of ASD traits in ADHD children. Interestingly, SDP had been found to be a marker of genetic predisposition than an independent environmental etiological factor of ADHD (Thapar et al., 2009; Skoglund et al., 2014). Indeed, one recent review suggested a genetic overlap between ADHD and SDP (van Amsterdam et al., 2018). Future studies could investigate in more details the shared genes between ADHD and SDP as those contributing to the separate development of ADHD and ASD, in addition to SHANK2 which contributes to their joint development.
This study has a number of limitations. First, this study only includes boys as participants. It is thus unsure whether our current findings are generalizable to girls. Second, this study relies on clinical diagnosis than on standardized diagnostic interviews. However, while the latter can provide standardization in assessment, clinical diagnosis remains the key criterion against which most, if not all, standardized diagnostic tools are initially validated. Third, in terms of functionality, there are no known functions for those SHANK2 SNPs which this study identifies. However, we do find our identified SNPs in strong LD with other SNPs involved in the transcriptional regulation or gene expression of SHANK2. Fourth, we have only genotyped common variants in the present study; subsequent exome sequencing of this Chinese children sample may uncover potential rare variants which drive the presently observed associations with the SHANK2 SNPs. Finally, our new finding that SHANK2 could be a potential pleiotropic gene for ADHD and ASD should be viewed with caution, if not suspicion. Particularly, our modest overall sample size may give rise to potential chance results, despite surviving multiple testing correction. Furthermore, the gene-level analysis is not entirely supportive with respect to ASD. However, even recent large ASD GWAS studies also failed to identify SHANK genes (The Autism Spectrum Disorders Working Group of The Psychiatric Genomics Consortium, 2017; Grove et al., 2019), despite the fairly well-replicated association of SHANK genes with ASD in previous association and animal studies, cited above (Peça et al., 2011; Won et al., 2012; Leblond et al., 2014; Connolly et al., 2017). Perhaps, these recent ASD GWAS studies are still underpowered, as we are. It is obviously too early to conclude on our new finding and yet dismissing it altogether may not be constructive either. Our study is very much in line with the recent efforts to detect pleiotropic genes and their mechanisms among multiple psychiatric disorders (Cross-Disorder Group of the Psychiatric Genomics Consortium, 2019). These efforts help researchers and clinicians to understand the common phenomenon of high comorbidity among psychiatric disorders. Future study should proceed to confirming our current findings with a much larger independent sample, as well as examining the functional aberrations of our identified SNPs/gene, using both in vivo and in vitro models.
In clinical practice, screening ADHD in children with ASD and vice versa should be routine. Because of their genetic overlap, they are more likely to be comorbid with each other. The previous hierarchical, mutually exclusive diagnostic practice had seriously misguided clinicians to miss the diagnosis of ADHD in ASD children. Such omission is most unfortunate, since ADHD is a treatable disorder by drug and behavioral intervention. If untreated, it will incur significant psychosocial impairment to a developing child. Our finding also has a significant impact on genetic counseling. In the past, the awareness of the elevated comorbidity between ASD and ADHD was not high. Based on the new knowledge from our study, families seeking genetic counseling will benefit from a better understanding of the genetic interconnectedness of these two disorders in order to make informed decisions.
Data Availability Statement
The dataset of this study which supports the findings and conclusions of the article will be available from the corresponding authors upon reasonable request.
Ethics Statement
The studies involving human participants were reviewed and approved by the Joint Chinese University of Hong Kong-New Territories East Cluster Clinical Research Ethics Committee. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
S-LM, C-CL, KL, S-FH, C-PT, T-PH, P-CS, and PL: conceptualization and study design. C-PT, CS, FM, and PL: fieldwork. LC, TM, and P-CS: formal analysis. LC: writing—original draft. S-LM, LC, C-CL, KL S-FH, C-PT, T-PH, CS, FM, TM, P-CS, and PL: writing—review and editing. All authors contributed to the article and approved the submitted version.
Funding
The part of the work on ADHD was supported by a General Research Fund (GRF) of the Hong Kong Research Grants Council (RGC) (RGC 449511), and the part of the work on ASD was supported by a Direct Grant from the Faculty of Social Science (4052042) as well as by a Seed Grant from the Institute of Innovative Medicine (8303104), the Chinese University of Hong Kong (CUHK). Funding sources play no role in study design, data collection, analysis, interpretation, manuscript preparation, or submission.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We are very grateful to the child and adolescent psychiatric services of three public hospitals in Hong Kong, namely, Kwai Chung Hospital, Alice Ho Mui Ling Nethersole Hospital, and United Christian Hospital, which assist in the recruitment of ADHD and ASD children. Special thanks must also be accorded to the hard work of our four former research assistants, Sharon Cheung, Fungyee Ching, Fanny Lau, and Fanny Fang. Finally, we are very thankful to the children and parents who so willingly participate in this project.
Footnotes
References
American Psychiatric and Association (2013). Diagnostic and Statistical Manual of Mental Disorders : DSM-5, 5th Edn. Washington, D.C.: American Psychiatric Association.
Boccuto, L., Lauri, M., Sarasua, S. M., Skinner, C. D., Buccella, D., Dwivedi, A., et al. (2013). Prevalence of SHANK3 variants in patients with different subtypes of autism spectrum disorders. Eur. J. Hum. Genet. 21, 310–316. doi: 10.1038/ejhg.2012.175
Boyle, A. P., Hong, E. L., Hariharan, M., Cheng, Y., Schaub, M. A., Kasowski, M., et al. (2012). Annotation of functional variation in personal genomes using RegulomeDB. Genome Res. 22, 1790–1797. doi: 10.1101/gr.137323.112
Carlson, C. S., Eberle, M. A., Rieder, M. J., Yi, Q., Kruglyak, L., and Nickerson, D. A. (2004). Selecting a maximally informative set of single-nucleotide polymorphisms for association analyses using linkage disequilibrium. Am. J. Hum. Genet. 74, 106–120. doi: 10.1086/381000
Cochoy, D. M., Kolevzon, A., Kajiwara, Y., Schoen, M., Pascual-Lucas, M., Lurie, S., et al. (2015). Phenotypic and functional analysis of SHANK3 stop mutations identified in individuals with ASD and/or ID. Mol. Autism 6:23. doi: 10.1186/s13229-015-0020-5
Connolly, S., Anney, R., Gallagher, L., and Heron, E. A. (2017). A genome-wide investigation into parent-of-origin effects in autism spectrum disorder identifies previously associated genes including SHANK3. Eur. J. Hum. Genet. 25, 234–239. doi: 10.1038/ejhg.2016.153
Craig, F., Margari, F., Legrottaglie, A. R., Palumbi, R., de Giambattista, C., and Margari, L. (2016). A review of executive function deficits in autism spectrum disorder and attention-deficit/hyperactivity disorder. Neuropsychiatric. Dis. Treat. 12, 1191–1202. doi: 10.2147/NDT.S104620
Cross-Disorder Group of the Psychiatric Genomics Consortium. (2019). Genomic relationships, novel loci, and pleiotropic mechanisms across eight psychiatric disorders. Cell 179, 1469–1482. doi: 10.1016/j.cell.2019.11.020
Demontis, D., Walters, R. K., Martin, J., Mattheisen, M., Als, T. D., Agerbo, E., et al. (2019). Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder. Nat. Genet. 51, 63–75. doi: 10.1038/s41588-018-0269-7
Durand, C. M., Betancur, C., Boeckers, T. M., Bockmann, J., Chaste, P., Fauchereau, F., et al. (2007). Mutations in the gene encoding the synaptic scaffolding protein SHANK3 are associated with autism spectrum disorders. Nat. Genet. 39, 25–27. doi: 10.1038/ng1933
Falk, C., and Rubinstein, P. (1987). Haplotype relative risks: an easy reliable way to construct a proper control sample for risk calculations. Ann. Hum. Genet. 51, 227–233. doi: 10.1111/j.1469-1809.1987.tb00875.x
Ghirardi, L., Brikell, I., Kuja-Halkola, R., Freitag, C. M., Franke, B., Asherson, P., et al. (2018). The familial co-aggregation of ASD and ADHD: a register-based cohort study. Mol. Psychiatry 23, 257–262. doi: 10.1038/mp.2017.17
Grove, J., Ripke, S., Als, T. D., Mattheisen, M., Walters, R. K., Won, H., et al. (2019). Identification of common genetic risk variants for autism spectrum disorder. Nat. Genet. 51, 431–444. doi: 10.1038/s41588-019-0344-8
Grzadzinski, R., Di Martino, A., Brady, E., Mairena, M. A., O’Neale, M., Petkova, E., et al. (2011). Examining autistic traits in children with ADHD: does the autism spectrum extend to ADHD? J. Autism Dev. Disord. 41, 1178–1191. doi: 10.1007/s10803-010-1135-3
Guilmatre, A., Huguet, G., Delorme, R., and Bourgeron, T. (2014). The emerging role of SHANK genes in neuropsychiatric disorders. Dev. Neurobiol. 74, 113–122. doi: 10.1002/dneu.22128
Karalunas, S. L., Hawkey, E., Gustafsson, H., Miller, M., Langhorst, M., Cordova, M., et al. (2018). Overlapping and distinct cognitive impairments in attention-deficit/hyperactivity and autism spectrum disorder without intellectual disability. J. Abnorm. Child Psychol. 46, 1705–1716. doi: 10.1007/s10802-017-0394-2
Leblond, C. S., Nava, C., Polge, A., Gauthier, J., Huguet, G., Lumbroso, S., et al. (2014). Meta-analysis of SHANK mutations in autism spectrum disorders: a gradient of severity in cognitive impairments. PLoS Genet. 10:e1004580. doi: 10.1371/journal.pgen.1004580
Li, M. X., Sham, P. C., Cherny, S. S., and Song, Y. Q. (2010). A knowledge-based weighting framework to boost the power of genome-wide association studies. PLoS One 5:e14480. doi: 10.1371/journal.pone.0014480
Li, Z., Chang, S. H., Zhang, L. Y., Gao, L., and Wang, J. (2014). Molecular genetic studies of ADHD and its candidate genes: a review. Psychiatry Res. 219, 10–24. doi: 10.1016/j.psychres.2014.05.005
Martin, J., Cooper, M., Hamshere, M. L., Pocklington, A., Scherer, S. W., Kent, L., et al. (2014). Biological overlap of attention-deficit/hyperactivity disorder and autism spectrum disorder: evidence from copy number variants. J. Am. Acad. Child Adolesc. Psychiatry 53, 761–770. doi: 10.1016/j.jaac.2014.03.004
Miller, M., Musser, E. D., Young, G. S., Olson, B., Steiner, R. D., and Nigg, J. T. (2019). Sibling recurrence risk and cross-aggregation of attention-deficit/hyperactivity disorder and autism spectrum disorder. JAMA Pediatr. 173, 147–152. doi: 10.1001/jamapediatrics.2018.4076
Moessner, R., Marshall, C. R., Sutcliffe, J. S., Skaug, J., Pinto, D., Vincent, J., et al. (2007). Contribution of SHANK3 mutations to autism spectrum disorder. Am. J. Hum. Genet. 81, 1289–1297. doi: 10.1086/522590
Monteiro, P., and Feng, G. (2017). SHANK proteins: roles at the synapse and in autism spectrum disorder. Nat. Rev. Neurosci. 18, 147–157. doi: 10.1038/nrn.2016.183
Mulligan, A., Anney, R. J. L., O’Regan, M., Chen, W., Butler, L., Fitzgerald, M., et al. (2009). Autism symptoms in attention-deficit/hyperactivity disorder: a familial trait which correlates with conduct, oppositional defiant, language and motor disorders. J. Autism Dev. Disord. 39, 197–209. doi: 10.1007/s10803-008-0621-3
Nemirovsky, S. I., Córdoba, M., Zaiat, J. J., Completa, S. P., Vega, P. A., González-Morón, D., et al. (2015). Whole genome sequencing reveals a de novo SHANK3 mutation in familial autism spectrum disorder. PLoS One 10:e0116358. doi: 10.1371/journal.pone.0116358
Nøvik, T. S., Hervas, A., Ralston, S. J., Dalsgaard, S., Rodrigues, P. R., Lorenzo, M. J., et al. (2006). Influence of gender on attention-deficit/hyperactivity disorder in Europe - ADORE. Eur. Child Adolesc. Psychiatry 15, (Suppl. 1), 1–24. doi: 10.1007/s00787-006-1003-z
Peça, J., Feliciano, C., Ting, J. T., Wang, W., Wells, M. F., Venkatraman, T. N., et al. (2011). Shank3 mutant mice display autistic-like behaviours and striatal dysfunction. Nature 472, 437–442. doi: 10.1038/nature09965
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A., Bender, D., et al. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575. doi: 10.1086/519795
Qin, J., Jia, M., Wang, L., Lu, T., Ruan, Y., Liu, J., et al. (2009). Association study of SHANK3 gene polymorphisms with autism in Chinese Han population. BMC Med. Genet. 10:61. doi: 10.1186/1471-2350-10-61
Reiersen, A. M., Constantino, J. N., Grimmer, M., Martin, N. G., and Todd, R. D. (2008). Evidence for shared genetic influences on self-reported ADHD and autistic symptoms in young adult Australian twins. Twin Res. Hum. Genet. 11, 579–585. doi: 10.1375/twin.11.6.579
Rommelse, N. N. J., Franke, B., Geurts, H. M., Hartman, C. A., and Buitelaar, J. K. (2010). Shared heritability of attention-deficit/hyperactivity disorder and autism spectrum disorder. Eur. Child Adolesc. Psychiatry 19, 281–295. doi: 10.1007/s00787-010-0092-x
Ronald, A., Larsson, H., Anckarsäter, H., and Lichtenstein, P. (2014). Symptoms of autism and ADHD: a swedish twin study examining their overlap. J. Abnorm. Psych. 132, 440–451. doi: 10.1037/a0036088
Rutherford, M., McKenzie, K., Johnson, T., Catchpole, C., O’Hare, A., McClure, I., et al. (2016). Gender ratio in a clinical population sample, age of diagnosis and duration of assessment in children and adults with autism spectrum disorder. Autism 20, 628–634.
Schmeisser, M. J., Ey, E., Wegener, S., Bockmann, J., Stempel, A. V., Kuebler, A., et al. (2012). Autistic-like behaviours and hyperactivity in mice lacking ProSAP1/Shank2. Nature 486, 256–260. doi: 10.1038/nature11015
Simonoff, E., Pickles, A., Charman, T., Chandler, S., Loucas, T., and Baird, G. (2008). Psychiatric disorders in children with autism spectrum disorders: prevalence, comorbidity, and associated factors in a population-derived sample. J. Am. Acad. Child Adolesc. Psychiatry 47, 921–929. doi: 10.1097/CHI.0b013e318179964f
Sinzig, J., Walter, D., and Doepfner, M. (2009). Attention deficit/hyperactivity disorder in children and adolescents With autism spectrum disorder: symptom or syndrome? J. Atten. Disord. 13, 117–126. doi: 10.1177/1087054708326261
Skoglund, C., Chen, Q., D’Onofrio, B. M., Lichtenstein, P., and Larsson, H. (2014). Familial confounding of the association between maternal smoking during pregnancy and ADHD in offspring. J. Child Psychol. Psychiatry 55, 61–68. doi: 10.1111/jcpp.12124
Tang, N. L. S., Pharoah, P. D. P., Ma, S. L., and Easton, D. F. (2006). Evaluation of an algorithm of tagging SNPs selection by linkage disequilibrium. Clin. Biochem. 39, 240–443. doi: 10.1016/j.clinbiochem.2005.11.014
Thapar, A., Rice, F., Hay, D., Boivin, J., Langley, K., Van Den Bree, M., et al. (2009). Prenatal smoking might not cause attention-deficit/hyperactivity disorder: evidence from a novel design. Biol. Psychiatry 66, 722–727. doi: 10.1016/j.biopsych.2009.05.032
The 1000 Genomes Project Consortium (2015). A global reference for human genetic variation. Nature 526, 68–74. doi: 10.1038/nature15393
The Autism Spectrum Disorders Working Group of The Psychiatric Genomics Consortium (2017). Meta-analysis of GWAS of over 16,000 individuals with autism spectrum disorder highlights a novel locus at 10q24.32 and a significant overlap with schizophrenia. Mol. Autism 8:21. doi: 10.1186/s13229-017-0137-9
The Brainstorm Consortium (2018). Analysis of shared heritability in common disorders of the brain. Science 360, 1313. doi: 10.1126/science.aap8757
van Amsterdam, J., van der Velde, B., Schulte, M., and van den Brink, W. (2018). Causal factors of increased smoking in ADHD: a systematic review. Subst. Use Misuse 53, 432–445. doi: 10.1080/10826084.2017.1334066
Waga, C., Okamoto, N., Ondo, Y., Fukumura-Kato, R., Goto, Y. I., Kohsaka, S., et al. (2011). Novel variants of the SHANK3 gene in Japanese autistic patients with severe delayed speech development. Psychiatr. Genet. 21, 208–211. doi: 10.1097/YPG.0b013e328341e069
Won, H., Lee, H. R., Gee, H. Y., Mah, W., Kim, J. I., Lee, J., et al. (2012). Autistic-like social behaviour in Shank2-mutant mice improved by restoring NMDA receptor function. Nature 486, 261–265. doi: 10.1038/nature11208
Young, S., Adamo, N., Asgeirsdottir, B. B., Branney, P., Beckett, M., Colley, W., et al. (2020). Females with ADHD: an expert consensus statement taking a lifespan approach providing guidance for the identification and treatment of attention-deficit/hyperactivity disorder in girls and women. BMC Psychiatr. 20:404. doi: 10.1186/s12888-020-02707-9
Yuen, E. (2017). A Public Child and Adolescent Psychiatric Clinic in Hong Kong: Characteristics of the Child and Adolescent Patients Seen. Unpublished MSSc thesis, Hong Kong: Department of Psychology, The Chinese University of Hong Kong.
Keywords: SHANK genes, ADHD, ASD, genetic overlap, pleiotropic gene
Citation: Ma S-L, Chen LH, Lee C-C, Lai KYC, Hung S-F, Tang C-P, Ho T-P, Shea C, Mo F, Mak TSH, Sham P-C and Leung PWL (2021) Genetic Overlap Between Attention Deficit/Hyperactivity Disorder and Autism Spectrum Disorder in SHANK2 Gene. Front. Neurosci. 15:649588. doi: 10.3389/fnins.2021.649588
Received: 05 January 2021; Accepted: 06 April 2021;
Published: 27 April 2021.
Edited by:
Adriana Lori, Emory University, United StatesReviewed by:
Lu Liu, Peking University Sixth Hospital, ChinaJubao Duan, University of Chicago, United States
Copyright © 2021 Ma, Chen, Lee, Lai, Hung, Tang, Ho, Shea, Mo, Mak, Sham and Leung. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Pak-Chung Sham, pcsham@hku.hk; Patrick W. L. Leung, pleung@cuhk.edu.hk
†These authors have contributed equally to this work and share first authorship
‡These authors have contributed equally to this work and share co-corresponding authorship