BAIAP2 exhibits association to childhood ADHD especially predominantly inattentive subtype in Chinese Han subjects

Abstract

Background

Attention-deficit/hyperactivity disorder (ADHD) is a common chronic neurodevelopmental disorder with a high heritability. Much evidence of hemisphere asymmetry has been found for ADHD probands from behavioral level, electrophysiological level and brain morphology. One previous research has reported possible association between BAIAP2, which is asymmetrically expressed in the two cerebral hemispheres, with ADHD in European population. The present study aimed to investigate the association between BAIAP2 and ADHD in Chinese Han subjects.

Methods

A total of 1,397 ADHD trios comprised of one ADHD proband and their parents were included for family-based association tests. Independent 569 ADHD cases and 957 normal controls were included for case-control studies. Diagnosis was performed according to the DSM-IV criteria. Nine single nucleotide polymorphisms (SNPs) of BAIAP2 were chosen and performed genotyping for both family-based and case-control association studies.

Results

Transmission disequilibrium tests (TDTs) for family-based association studies showed significant association between the CA haplotype comprised by rs3934492 and rs9901648 with predominantly inattentive type (ADHD-I). For case-control study, chi-square tests provided evidence for the contribution of SNP rs4969239, rs3934492 and rs4969385 to ADHD and its two clinical subtypes, ADHD-I and ADHD-C. However, only the associations for ADHD and ADHD-I retained significant after corrections for multiplicity or logistic regression analyses adjusting the potential confounding effect of gender and age.

Conclusions

These above results indicated the possible involvement of BAIAP2 in the etiology of ADHD, especially ADHD-I.

Background

Attention-deficit/hyperactivity disorder (ADHD) is a common neurodevelopmental disorder. The high heritability of approximately 0.76 [1] has suggested the important role of genetic factors in its etiology.

The typical development of human cerebral asymmetry has profound effect on the normal lateralization of cognitive and motor functions, such as language and handedness. Several studies have suggested the possible link of the disturbance of cerebral asymmetry with the pathogenesis of neurodevelopmental disorders [2], such as schizophrenia [3], autism [4] and ADHD [5].

The non-right-handedness, especially mixed-handedness, might be a risk indicator for the symptom severity of inattention in ADHD [6]. One of our previous behavioral studies indicated that ADHD children showed an atypical pattern of right hemisphere in conflict control task compared to controls [7]. Another study demonstrated that the direction of spatially asymmetrical interference effects in ADHD was opposite to controls, indicating disruption within right hemisphere attentional networks [8]. The structural imaging findings in ADHD has demonstrated volumetric reductions in total and right cerebral, right caudate, cerebellar and corpus callosum [9]. ADHD adults had thinner cortex in the cortical networks, especially in the right hemisphere involving inferior parietal lobule, dorsolateral prefrontal and anterior cingulate cortices [10]. In addition, some evidence has showed that the critical feature of ADHD is the delayed maturation not only of prefrontal cortical thickness but also of cortical surface area, especially the right hemispheric lobes [11, 12]. The degree of rightward volumetric asymmetry in caudate nucleus might predict the severity of inattentive symptoms of ADHD [13]. An fMRI study of ADHD children showed under-activation of the right caudate nucleus and inferior parieta cortex [14].

The gene BAIAP2, which is located on 17q25 and encodes brain-specific angiogenesis inhibitor 1-associated protein 2 (BAIAP2), has been suggested to be involved in cerebral asymmetry [15]. Recently, a genetic study of adulthood ADHD in two independent European populations suggested an association of BAIAP2 with ADHD, supplying genetic evidence of abnormal left-right brain asymmetries with this disorder [16]. In addition, a genome-wide association study (GWAS) has indicated nominal association between BAIAP2 and its isoform BAIAP2L1 with ADHD [17]. Moreover, methamphetamine, one of psychostimulants being considered as first-line pharmacological treatments for ADHD patients, has enhanced the expression of BAIAP2 in rat cerebral cortices [18]. BAIAP2 has also been found to confer risk for autism spectrum disorders (ASD), which shared some genetic risk factors with ADHD [19].

For the potentially important role of cerebral asymmetry in the pathophysiology of ADHD, BAIAP2 has been suggested to be one of novel candidate genes for ADHD, but need more work for replication [20]. Our present study is to investigate the relationship between BAIAP2 and ADHD in Chinese Han subjects.

Methods

Subjects

All ADHD cases were recruited from child psychiatric clinics of Peking University Institute of Mental Health. Diagnosis was performed according to the DSM-IV criteria by experienced child psychiatrists, using the Clinical Diagnostic Interview Scale (CDIS) [21], which was translated into Chinese by our groups before [22, 23]. The CDIS assesses the three DSM-IV subtypes of ADHD, including ADHD inattentive type (ADHD-I), ADHD hyperactive-impulsive type (ADHD-HI) and ADHD combined type (ADHD-C). It was also used to evaluate comorbidities, including oppositional defiant disorder (ODD), tic disorder (TD), learning disorder (LD), etc. The diagnosis of LD is just based on a brief parents’ report of the general academic achievement, but not including detailed assessment for reading, writing or arithmetic abilities. So dyslexia, dysgraphia or dyscalculia was not defined in our current study. All cases were of Chinese Han descent, with age between 6 and 16 years, and full-scale estimated IQ > 70. Any major neurological disorders, a diagnosis of schizophrenia, pervasive development disorder, epilepsy, mental retardation or other brain disorders were excluded. Finally, a total of 1,966 ADHD probands were included. Among these cases, 1,397 ADHD probands, along with their parents, constituted trios for family-based association analyses. The other independent 569 ADHD cases were included for case-control study.

The control sample consisted of 957 subjects from local elementary schools, healthy blood donors from the blood center of the First Hospital of Peking University, and healthy volunteers at our institute. ADHD, other major psychiatric disorders, family history of psychosis, severe physical diseases and substance abuse were excluded (more details have been described in our previously published article [24]). Demographic and clinical characteristics of both ADHD cases and control sample are shown in Table 1.

Table 1 Demographic and clinical characteristics of sample

All subjects were treated according to the Declaration of Helsinki and this work was approved by the Ethics Committee of Peking University Health Science Center. Written informed consent was obtained from each subjects or parents of children.

DNA isolation

Peripheral blood was collected for all subjects included in this study. Then genomic DNA was extracted following the standard protocols using E.Z.N.A.™ Blood DNA Kits (Omega Bio-tek Inc., Doraville, GA).

SNP selection and genotyping

For selection of single nucleotide polymorphisms (SNPs), Haploview version 4.2 was used to pick up tag SNPs based on the CHB database from Hapmap (http://www.hapmap.org). Threshold limit was set of r² > 0.80 for all SNPs with minor allele frequency (MAF) > 0.25. Thirteen tag SNPs were chosen with these criteria. However, only eight tag SNPs were included for genotyping and further association analyses. The other five SNPs were not selected for the location at the same haplotype block or strong to moderate linkage disequilibrium (LD) with those selected SNPs (Table 2). We also included an additional SNP rs4969385 located on intron 6, which was found to be potentially associated with adult ADHD in Spanish population [16].

Table 2 Selection of SNPs

Genotyping for all SNPs were carried out on an ABI 7900-HT instrument (Applied Biosystems, Foster City, USA), using Taqman allelic genotyping assays [25] and following the standard protocol as described by our previous report [26]. The SDS version 2.3 software was used for genotype identification. For quality control, firstly, two to four non-template controls (NTC) were set on each 384-well plate with no genotype called. Secondly, 3% samples were selected randomly and genotyped for the same SNP by different experimenters, indicating the concordance rate of 100%. Call rates for all SNPs were ranging from 98.8% to 99.5%.

Statistical analyses

Calculation of MAF and analysis of Hardy-Weinberg equilibrium (HWE) were performed using Haploview 4.2 software for cases and controls separately, showing no departure from HWE for any SNP (P-value from 0.178 to 0.921). The MAFs for all SNPs were from 0.162 to 0.499. Haploview software was also used to estimate the linkage disequilibrium (LD) of all SNPs and generate haplotype blocks which will be used for multi-marker haplotype-based association tests [27]. As shown by Figure 1, there were three haplotype blocks generated including Block 1 comprised of SNP 1–4, Block 2 comprised of SNP 5–6 and Block 3 comprised of SNP 8–9. However, the SNP7 rs4076037 was not included into any block and subsequently not included for multi-marker analyses. For Block 1, only three of seven generated haplotypes (ACAA, GAGG and ACGG) were used for haplotype analyses, which captured 85.3% of the genetic variance in the investigated four SNPs, while the other four haplotypes with frequencies < 0.1 were not included.

Figure 1

Linkage disequilibrium (LD) plot of BAIAP2 estimated by Haploview software. It displayed LD value in D’ (i.e., 99 indicated D’ of 0.99. Squares with no number indicates a D’ of 1). Generated haplotypes were listed on the right for each of three LD blocks with frequencies.

For family-based association study, TDT tests were conducted using Haploview to investigate whether there was biased-transmission of alleles or haplotypes in ADHD trios. For case-control study, chi-square tests were conducted using Haploview for both single-marker and multi-marker haplotype-based analyses, to compare the frequencies of alleles and haplotypes between ADHD cases and controls. In addition, comparison of the genotype frequencies was also performed for case-control statistics using SPSS software. Once there was nominally significant association appeared either for alleles or genotypes under additive model (nominal P < 0.05), further analyses for genotypes based on dominant and recessive models were conducted [16]. As indicated in Table 1, the comparisons of distribution of gender and age between two groups showed statistical significance. Then, we further developed logistic regression analyses to control the potential confounding effect of gender and age.

For multiple testing corrections of single-marker analyses, Bonferroni corrections were conducted. We firstly used SNPSpD software [28] to consider the inter-correlation of SNPs for multiple testing of SNPs in LD with each other. For our current study, the effective number of independent marker loci yielded by SNPSpD was 7. Then, for TDT tests, taking into account three clinical subtypes and the effective number of SNPs (n = 7), the experiment-wide significance threshold required to keep type I error rate at 5% was set at P < 0.0024. For case-control study, considering the effective number of SNPs (n = 7), three clinical subtype, and the comparison of genotype and allele frequencies, the adjusted significance was set at P < 0.0012. For multi-marker association tests, haplotype analysis was only conducted when nominal association existed in the single-marker analyses and significance was corrected by 5,000 permutations using Haploview software.

The minimum statistically genetic power for TDT tests and case-control studies was calculated using the Genetic Power Calculator software (http://pngu.mgh.harvard.edu/~purcell/gpc/), taking the lowest MAF of 0.162 and assuming the prevalence of 0.05, significance level of 0.05 and odd ratio (OR) of 1.5. Under above settings, our sample showed the minimum statistically genetic power of 80.9% for TDT statistics with about 500 efficient trios, and 74.4% for 569 cases vs. 957 controls under additive model.

Results

Family-based association tests

Single-marker analyses

In the general ADHD trios, we did not find any biased transmission of any allele for any SNP (all P > 0.05, Table 3). Further analyses were conducted for three subtypes separately. For ADHD-I trios, TDT tests showed that the C allele of SNP rs3934492 was over-transmitted with nominal P-value of 0.030 (Table 3). No biased-transmission was found either in ADHD-HI or in ADHD-C trios (all P values > 0.05, Table 3).

Table 3 TDT tests for ADHD trios and its three clinical subtypes ^a

Multi-marker haplotype-based association tests

Haplotype analyses were only conducted for ADHD-I trios. The CA haplotype of Block 2 coded by rs3934492 and rs9901648 was over-transmitted (χ²=11.18, nominal P=8.0e-4, empirical P=0.005) and the GA haplotype was potentially under-transmitted (χ²=4.66, nominal P = 0.031, empirical P=0.207) (Table 4).

Table 4 TDT tests for haplotypes of Block2 with ADHD-I trios (n = 706)

Case-control studies

Single-marker analyses

For allelic analyses, nominal association was displayed with ADHD for SNP rs4969385 (χ²=5.59, P = 0.018, OR = 1.27 [1.04-1.55]) (Table 5). What’s more, the potential confounding factors of subtypes were considered by analyses in ADHD-I and ADHD-C separately, while ADHD-HI was not analyzed due to its small sample size. Two SNPs showed association with ADHD-I: rs4969239 (χ²=6.58, P = 0.010, OR = 1.29 [1.06-1.58]) and rs3934492 (χ²=4.28, P = 0.039, OR = 1.22 [1.01-1.48]). One SNP rs4969385 showed association with ADHD-C (χ²=4.14, P = 0.042, OR = 1.31 [1.01-1.71]) (Table 5).

Table 5 Analyses of alleles and genotypes for 569 ADHD (295 ADHD-I, 252 ADHD-C) and controls (n = 957)^a

Genotypic analyses also support association between above SNPs with ADHD, ADHD-I and ADHD-C (Table 5). However, only the association of rs4969385 with ADHD (P = 0.001, OR = 2.96 [1.51-5.81]) and ADHD-I (P = 0.001, OR = 3.48 [1.64-7.39]) retained significant after corrections for multiplicity, but not for ADHD-C (P = 0.040, OR = 2.64 [1.13-6.17]) (Table 5). Logistic regression analyses also showed retained association of BAIAP2 with ADHD (P < 0.05) and ADHD-I (P < 0.05) after adjusting the potential effect of gender and age, but not for ADHD-C (P = 0.086) (Table 6).

Table 6 Logistic regression analyses of genotype distribution between ADHD and controls ^a

Multi-marker haplotype-based association tests

Based on the association from single-marker analyses, we performed haplotype analyses for ADHD, ADHD-I and ADHD-C. As indicated in Table 7, the CT haplotype of Block3 showed nominal association with ADHD in general (χ²=5.44, nominal P = 0.019, empirical P = 0.097). For ADHD-I, the GA haplotype of Block2 showed lower frequency (χ²=3.98, nominal P = 0.046, empirical P = 0.217) in cases than controls, while the CA haplotype showed higher frequency (χ²=5.20, nominal P = 0.023, empirical P = 0.110). For ADHD-C, the CT heplotype of Block3 showed hither frequency in cases (χ²=3.96, nominal P = 0.047, empirical P = 0.215). However, none of above association remained significant after permutations.

Table 7 Haplotype analyses for 569 ADHD versus 957 controls

Discussion

Our results showed that BAIAP2 was associated with childhood ADHD of Chinese Han descent, especially for the predominantly inattentive type (ADHD-I). Analyses for case-control studies indicated different alleles and genotypes distribution between ADHD and controls, and the association was only retained for ADHD and ADHD-I after multiple corrections or adjusting the confounding effect of gender and age, but not for ADHD-C. From TDT tests for trios, the significant association was only indicated for ADHD-I from haplotype analyses.

In some extent, our current findings are consistant with previous reports. Ribasés et al. [16] have analyzed six functional candidate genes showing differential expression between hemispheres in ADHD and normal controls, but only found an association between BAIAP2 with ADHD. However, in their study, the association of BAIAP2 with ADHD was only observed in adults, which has also been replicated in an independent population, but not children. The possible explanation for this discrepancy may be that the genetic load for ADHD children may be lower than ADHD adults, while the sample size in the study by Ribasés et al. [16] may be not enough to detect the relatively weak association in ADHD children. Further investigation in ADHD adults of Chinese Han descent especially by follow-up studies may promote our understanding of the explicit effect of BAIAP2 on ADHD.

In our present study, the SNP rs4969385 was the only associated one with ADHD and ADHD-I after corrections for multiple testing. This SNP has also been reported in adult samples by Ribasés et al. [16]. This replication suggested further exploration of this SNP and its related functional variants in the predisposition to ADHD. However, we note that the direction of the effect in our current study was opposite to the one observed by Ribasés et al. [16]. In our study, the minor ‘T’ allele showed risk for ADHD, while the major ‘C’ allele did in Ribasés’s report. Previous studies has also reported similar phenomenon between different ethnicities [29, 30]. In addition, another SNP rs8079626, showing none association in our current study, has been reported to be of nominal association with adult ADHD in German sample by Ribasés et al. [16]. In their study, the ‘A’ allele showed higher frequencies in German adult ADHD sample than controls. In our study, for the ‘G’ allele, the transmission was more than non-transmission from TDT tests and its frequency was higher in ADHD probands versus controls, although not achieving statistically significant difference. When checking the LD between rs8079626 with the best two SNPs of our current analyses, we found its low LDs with rs3934492 (D’ = 0.231, r² = 0.04) and rs4969385 (D’ = 0.072, r² = 0.002). So we can conclude that the SNP rs8079626 was not associated with ADHD in Chinese Han subjects. This discrepancy may be greatly due to the ethnic difference, that its allele distribution in our current analyses of Chinese Han subjects (62.8% of G, 37.2% of A) was adverse to that of German samples (31.0% of G, 69.0% of A) from the study by Ribasés et al. [16]. This adverse distribution is according with the Hapmap database (62.2% of G for CHB, 28.3% of G for CEU). Taken together, it is very important to consider potential pathogenic genetic variants in different ethnic populations for ADHD etiological studies.

Another interesting finding is the specific association of BAIAP2 with ADHD-I subtype. Three subtypes of ADHD were defined in DSM-IV including predominantly inattentive, predominantly hyperactive/impulsive and the combined subtypes [31]. A growing body of literature has addressed that inattentive and combined subtypes may be distinct disorders [32, 33]. Attention-deficit disorder (ADD; inattentive, without hyperactivity) was different from ADHD on genetic basis, comorbidity, related brain region, etc. (for a review, see [33]). These findings have some effect on the revision of ADHD subtypes in DSM-V [34]. Our previous studies on molecular genetics have revealed the potential subtype-specific genes associated with ADHD in the Chinese Han population. COMT, 5-HT1B, MAOA, CHRNA4, SYP and DDC exhibited associations primarily with ADHD-I, whereas HTR2C, 5-HT1D, ADRA2C, DRD3 and NET1 were associated mainly with ADHD-C [24, 35–38]. The family-based study of NET1 and ADHD strongly suggested that the haplotype blocks within different regions of NET1 show divergent association based on sex and subtype [39]. Taken the current finding together with previous reports, considering subtypes in ADHD genetic studies is very important to reduce the heterogeneity and may help us to explore the hidden existed genes associated with ADHD.

Our findings need to be considered in light of some limitations. Firstly, twenty-seven genes have been identified involving in the left-right asymmetric cortical development in humans [15]. We only included BAIAP2 in this study and whether other genes played the role for the pathogenesis of ADHD remained unclear. Secondly, we have set MAF > 0.25 for tag SNPs selection and included an additional reported SNP rs4969385 with MAF of 0.162. We could not preclude the rare genetic contribution of uncommon SNPs (MAF <0.162) to the occurrence of ADHD. Thirdly, we have not screened those parents for ADHD diagnosis that we could not explain the genetic mechanisms of BAIAP2 in the etiology of ADHD clearly. Fourthly, ADHD was associated with high rates of comorbidity including dyslexia (reading disorder). In the existing literature, evidence has support the common risk neurobiological phenotype of atypical cerebral asymmetry (ACA) and shared ACA genes for ADHD and dyslexia [40]. The diagnosis instrument as described above in the present study lacked the ability to differ ADHD from dyslexia. Lastly, we did not collect data on the handedness, which is also strongly correlated with cerebral asymmetry. Then, we could not preclude the confounding effect of the diagnosis of dyslexia and handedness on our observed associations.

Conclusions

The present family-based and case-control association studies in our Chinese Han populations provide further evidence for the role of BAIAP2 in the predisposition to ADHD, especially ADHD-I. Further studies should investigate the involvement of hemispheric asymmetry genes modulating the difference of left-right hemisphere in ADHD and related cognition.