Abstract
GABAA receptor subunit genes GABRB3, GABRA5, and GABRG3 located on chromosome 15q11-q13 have been implicated in the etiology of autistic spectrum disorders (ASD). This study intended to investigate the possible role of single-nucleotide polymorphisms (SNPs) present in GABRB3 (rs2081648 and rs1426217), GABRA5 (rs35586628), and GABRG3 (rs208129) genes in ASD susceptibility and symptom-based and developmental phenotypes of ASD in Chinese Han children and adolescents. 99 ASD patients and 231 age- and gender- frequency-matched typical developing (TD) controls were tested by TaqMan® genotyping assay. Symptom-based phenotypes were evaluated by Childhood Autism Rating Scale (CARS) and Autism Behavior Checklist (ABC), and developmental phenotypes were assessed by Early Childhood Development Questionnaire (ECDQ) in ASD patients. Three haplotypes and global χ 2 test of all SNPs demonstrated significant associations between ASD and TD groups. Besides, GABRB3 rs2081648, GABRA5 rs35586628, and GABRG3 rs208129 polymorphisms were associated with symptom-based deficits in social interaction, sensorimotor and somatosensory coordination, visual response, imitation, activity level, language expression and adaptability. Developmental abnormalities in late emergences of social interaction and fine motor were detected in GABRB3 rs2081648 polymorphism. Overall results indicated that gene synergy may participate in ASD pathogenesis, and GABAA receptor gene polymorphisms can predict symptom-based and developmental deficits in ASD individuals.
Similar content being viewed by others
Introduction
Autistic spectrum disorders (ASD) are a constellation of heterogeneous neurodevelopmental disorders characterized by early-onset deficits in social communication and interaction, and by restricted and repetitive behaviors, interests, or activities1. It is reported that the prevalence of ASD is on rise from 1/2000 in 1970s to 1/200 in 2000 s2 and a recent survey indicated that approximately 1 in 68 children has ASD, with a male: female ratio of 4.5:13. ASD are typically considered as a group of life-long disorders, with heavy care and financial burdens on families and society. Nevertheless, there are still no doubtless pharmacological treatments to alleviate the core deficits in individuals with ASD. The obtainment of better targeted, ASD-specific therapies will only be possible according to a better understanding of ASD pathogenesis.
Genetics has a key role in the etiology of ASD, in conjunction with developmentally early environmental factors. However, these effects have not yet been confirmed. Emerging evidence suggests that increased excitatory and reduced inhibitory neurotransmission may form a final common pathway in ASD4. Brain hyperexcitability with abnormalities in GABA transmission5, 6 and lower brain gamma-aminobutyric acid (GABA) levels7,8,9, as well as reduced gamma-aminobutyric acid type A (GABAA) receptors in the superior frontal cortex of ASD patients10 have been assumed to underlie the symptoms of ASD.
GABA is the predominant inhibitory neurotransmitter in the adult brain, mainly acting via an intricate series of ionotropic GABAA receptors (ligand-gated chloride channels) on the postsynaptic neuron. The GABAA receptors are composed of 19 different subunits (α1–6, β1–3, γ1–3, δ, ε, θ, π, ρ1–3) arranged around a central pore and mediate the majority of fast synaptic inhibition in the brain. Activation of GABAA receptors is associated with impaired long-term potentiation (LTP) and impaired learning in vivo 11. Although the GABAA receptors contain the sites for several therapeutic drugs and agents, such as benzodiazepines, steroids, and anesthetics12, evidence supporting the benefit of GABAergic drugs in ASD is inconclusive and limited. GABAA receptor subtypes may offer the promise of a new CNS pharmacology beyond classical benzodiazepines on the basis of the regulation of cognitive behavior by α5 GABAA receptors13.
Converging genetic evidence specifically implicate the involvement of a cluster of GABAA receptor subunit genes (GABRB3, GABRA5 and GABRG3) located on the 15q11-q13 (β3, α5 and γ3 subunits) region in the pathogenesis of ASD14, 15. These three GABAA receptor subunit genes are physically positioned in the region on chromosome 15q which is most commonly reported loci of chromosomal abnormalities documented in patients with ASD16, 17, including deletions and duplications of chromosome 15q11-q13. For instance, maternal deletion of 15q11-q13 is responsible for Angelman syndrome which is characterized by impaired language and speech development, movement disorder, and mental retardation, while paternal deletion of 15q11.2-q12 is related to Prader-Willi syndrome which is characterized by hypotonia, short stature, and obesity. Both Angelman and Prader-Willi syndromes are liable to have ASD18.
Recently, both genomic linkage screens19,20,21,22,23 and linkage disequilibrium (LD)24,25,26,27,28,29 analyses have implicated GABRB3 as an excellent candidate gene for ASD. Nevertheless, the association between GABRB3 and ASD has not been universally identified30,31,32,33,34. Symptom-based phenotypes of ASD have offered some evidence for association with the GABRB3 region. Affected individuals with ASD having high insistence-on-sameness scores exhibited a higher linkage signal to the GABRB3 region among families21. GABRB3 deficient mice showed occasional epilepsy, hyper-responsive to human contact, diminished nurturing behaviors, hyperactive run in tight circles, poor motor skills, electroencephalographic abnormalities, impaired social and exploratory behaviors, hypoplasia of cerebellar vermis and learning and memory deficits35,36,37,38,39,40. However, the association between the GABRB3 region and savant skills is still contradictory with the positive41 and negative42 findings in ASD patients.
In addition to evidence implicating the GABRB3 region, there exists relative fewer supports for association between ASD and the GABRA5 region27, 43, as well as between ASD and the GABRG3 region44. Most findings were negative in the GABRA5 and GABRG3 regions with ASD15, 28, 32, 34, 45, 46. However, Kim et al.46 found a nominally significant association between one single-nucleotide polymorphism (SNP) of GABRA5 and a symptom-based phenotype defined as ‘relative failure to initiate or sustain conversational interchange’ in the form of inflexible language behavior47. Moreover, GABRA5 deficient mice exhibited deficits in short-term memory when the task became increasingly more difficult. Reduced expression and function of GABRA5 may cause neurodevelopmental changes that contribute to ASD-like behaviors48.
There exists most family-based association analysis studies of GABAA receptor genes on chromosome 15q11-q13 for reports of the positive association between maternal interstitial duplication and ASD, as well as significant associations and linkage studies in chromosomally normal ASD families41, but case-control studies are rare and little is known about the degree to which genetic polymorphisms underlie symptom-based and developmental variability in ASD patients. To replenish ongoing efforts to describe allele associations at GABRB3, GABRA5 and GABRG3 in ASD family trios, we sought to utilize a case-control association study to observe relevant clues of ASD pathogenesis, as well as analyzing the associations between the GABRB3 (rs2081648 and rs1426217), GABRA5 (rs35586628), and GABRG3 (rs208129) gene polymorphisms and symptom-based and developmental deficits of ASD patients in Chinese Han children and adolescents.
Results
Association analysis of GABRB3, GABRA5, and GABRG3 SNPs
No significant deviations from the Hardy-Weinberg equilibrium were found in both ASD patients and typical developing (TD) controls for four SNPs, including two GABRB3 SNPs (rs2081648 and rs1426217), one GABRA5 SNP (rs35586628), as well as one GABRG3 SNP (rs208129). The genotypic and allelic frequencies of the four SNPs between the ASD patients and TD controls demonstrated no statistical differences (p > 0.05) (Table 1).
Linkage disequilibrium (LD) analysis showed two GABRB3 SNPs (rs2081648 and rs1426217) were in strong disequilibrium in ASD patients (D′ = 0.795) and in positive disequilibrium in TD controls (D′ = 0.532). Four SNPs of GABRB3 (rs2081648 and rs1426217), GABRA5 (rs35586628), and GABRG3 (rs208129) formed eleven effective haplotypes in both ASD and TD groups, and three haplotypes were associated with ASD (Table 2). The global p was 0.007 (χ 2 = 24.33, df = 10) among four SNPs between ASD and TD groups.
GABRB3, GABRA5, and GABRG3 SNPs and ASD symptom-based phenotypes
For the CARS total scores deviated from a normal distribution (p < 0.05), we conducted a Kruskal-Wallis H test to analyze the relationships between CARS total scores and the genotypes of four SNPs, and found no evidence for any associations. Nevertheless, an ordinal polytomous logistic regression analysis showed some significant associations between ratings in items of CARS and genotypes of GABRB3 SNP rs2081648, GABRA5 SNP rs35586628, and GABRG3 SNP rs208129 in ASD patients (while adjusting for sex, age, and IQ). ORs (95% CI) for the ‘4 = severely abnormal’ among different genotypes were shown in Table 3. ASD patients with a TT genotype of GABRB3 SNP rs2081648 trended to suffer from severe abnormality of ‘visual response’ compared to those with a TC genotype. Moreover, ASD participants with CC and TT genotypes of GABRA5 SNP rs35586628 were likely to perform worse in ‘verbal communication’ than those with a CT genotype, and ASD individuals with TA and TT genotypes of GABRG3 SNP rs208129 showed worse performance in ‘imitative behavior’ and ‘activity level’ than those with an AA genotype.
Scores of each ABC category and ABC total scores conformed to a normal distribution and homogeneity of variance (p > 0.05), but no significant associations between genotypes of four SNPs and all ABC scores were observed after the one-way ANCOVA with sex, age, and IQ as covariates (p > 0.05). However, a binary logistic regression analysis found some significant associations between ratings in items of ABC and genotypes of GABRB3 SNP rs2081648, GABRA5 SNP rs35586628, and GABRG3 SNP rs208129 in ASD patients (while adjusting for sex, age, and IQ). ORs (95% CI) for ‘1 = yes’ of existing relevant symptoms among different genotypes were shown in Table 4.
GABRB3, GABRA5, and GABRG3 SNPs and developmental phenotypes of ASD
An ordinal polytomous logistic regression analysis detected significant associations between ratings in two items of ECDQ and genotypes of GABRB3 rs2081648 in ASD patients (while adjusting for sex, age, and IQ) (n = 9 TT, n = 31 TC, and n = 27 CC). The OR for ‘3 = late’ of ‘smiling to his/her mother’ in those with a TT genotype compared to those with a CC genotype was 4.39 (95% CI = 1.43–13.52, p = 0.010), those with a TC genotype compared to those with a CC genotype was 0.21 (95% CI = 0.09–0.48, p = 0.000), and those with a TT genotype compared to those with a TC genotype was 20.86 (95% CI = 3.39–128.50, p = 0.001). We found that ASD patients with a TT genotype trended to show social interaction in older ages. In addition, the OR for ‘3 = late’ of ‘grasping things by himself/herself’ in those carrying a TT genotype compared with those carrying a CC genotype was 2.26 (95% CI = 0.86–5.91, p = 0.098), those carrying a TC genotype compared with those carrying a CC genotype was 0.28 (95% CI = 0.13–0.63, p = 0.002), and those carrying a TT genotype compared with those carrying a TC genotype was 8.02 (95% CI = 1.59–40.56, p = 0.012). As a result, ASD participants with a TC genotype developed earlier in fine motor.
Discussion
To our knowledge, this is the first age- and gender- frequency-matched case-control study designed to investigate the association between the three GABAA receptor genes, GABRB3, GABRA5, and GABRG3, and ASD in Chinese Han children and adolescents, as well as to detect relevant symptom-based and developmental deficits associated with the three GABAA receptor gene polymorphisms in Chinese Han children and adolescents with ASD.
Previous studies have investigated the association of the three GABAA receptor genes with ASD, but the contradictory findings demand subsequent replications in different populations. In this study, we found no positive evidence of the associations between any single SNP of the three GABAA receptor genes and ASD. Similar to our result, McCauley et al.27 demonstrated the negative evidence for the associations between three SNPs (rs2081648, rs1426217, and rs208129) and ASD, as well as Kim et al.46 reported no significant supports for four SNPs (rs2081648, rs1426217, rs35586628, and rs208129) of ASD mainly in Caucasian population with family trios. However, Kim et al.29 found an allele at SNP rs2081648 exhibited prior transmission in Korean trios. These differences may originate from the multiple demographic characteristics of the ASD populations (age, gender, ethnicity, sample size, diagnostic criteria, etc.) and different investigative methods (family-based association, case-control association, etc). The statistical power of our case-control association study is approximately 0.65 (α = 0.05) under the assumptions of 0.28% prevalence of ASD in Tianjin of China49 and genotypic relative risk of 1.834. Thus small effect might not be detected because of the sample size.
Moreover, we detected the TD controls had lower D′ than ASD patients, which predicted that TD individuals might have more chance to restructure randomly in accordance with genetic equilibrium than ASD individuals. Further studies can apply family-based association test (FBAT) in Chinese Han family trios of ASD and TD children and adolescents with the three GABAA receptor genes to discover more significant evidence on this difference. In addition, our study revealed that individuals with a CATT haplotype and a TGCA haplotype were more likely to have ASD, but individuals with a TACA haplotype had more potential to be TD controls. Our findings supported that the etiology of ASD may include gene-gene interactions rather than the effect of one single gene50.
For the limited evidence on the association between GABAA receptor genes and characteristic phenotypes of ASD, we further observed some significant associations between three SNPs of the three GABAA receptor genes and phenotypic features of ASD without GABRB3 SNP rs1426217. In the symptom-based phenotypes, as a supplement for the detections of associations between GABRB3 region and high degree of insistence-on-sameness21, as well as savant skills41 in previous studies, the ASD patients in our study with a TT genotype of GABRB3 SNP rs2081648 trended to perform worse in ‘visual response’ of CARS, and those with a CC genotype of GABRB3 SNP rs2081648 were more likely to show a deficit in social interaction according to the ‘relating’ part of ABC (38 = (R) Has not developed any friendships). Moreover, ASD patients with a CC genotype of GABRA5 SNP rs35586628 exhibited a deficit in language expression or cognitive functioning basing on the ‘verbal communication’ of CARS and the ‘language’ part of ABC (37 = (L) Cannot point to more than five named objects), and the others with a CT genotype showed sensorimotor and somatosensory abnormalities (13 = (R) Does not reach out when reached for; 29 = (S) Sometimes shows no ‘startle response’ to loud noise), and rigid language patterns (48 = (L) Repeats sounds or words over and over) according to the ‘relating’, ‘sensory’ and ‘language’ parts of ABC. Our finding provided more supports for the relationships between genotypes of GABRA5 and core features of ASD than Kim et al.46, who only reported a significant association between one GABRA5 SNP rs2075716 and ‘relative failure to initiate or sustain conversational interchange’. Furthermore, the ASD patients with a TA genotype of GABRG3 SNP rs208129 showed the highest potential to bear the most severe deficits in ‘imitative behavior’ and ‘activity level’ of CARS, as well as to display insistence-on-sameness or worse adaptability (14 = (s′) Strong reactions to changes in routine/environment), somatosensory disability (26 = (S) Sometimes painful stimuli such as bruises, cuts, and injections evoke no reaction) and difficulties in social communication (56 = (L) Uses at least 15 but less than 30 spontaneous phrases daily to communicate) in the ‘social and self-help’, ‘sensory’ and ‘language’ parts of ABC. Above all, the heterozygotes of GABRA5 SNP rs35586628 and GABRG3 SNP rs208129 showed the primary potential to increase the probability for ASD patients to undertake relevant core deficits, learning disabilities, and especially the decreased sensory sensitivity in voice and pain conformed abnormal sensory processing as one of the core features of individuals with ASD1. In the developmental phenotypes, ASD patients with a TT genotype of GABRB3 SNP rs2081648 were likely to show negative effects on initial social interaction and fine motor with the late emergences of ‘smiling to his/her mother’ and ‘grasping things by himself/herself’. As far as we know, there has been no previous studies intending to distinguish different phenotypes of early childhood development based on SNPs in children and adolescents with ASD. However, during the fetus period, the GABAA receptor subunit cluster (containing GABRB3, GABRA5, and GABRG3) on chromosome 15q11-q13 region acts as a developmental role in GABAergic signaling for building neuronal connectivity, and plays a pivotal role in the maintenance of inhibitory tone in the adult brain41. Therefore, any abnormalities in the three genes are possible to cause developmental deficits in the whole life. Subsequent researches should pay more attention to the associations between polymorphisms of GABAA receptor genes and early childhood development in both ASD patients and TD controls to obtain better and precise clues of early diagnosis for ASD patients, so as to intervene at an earlier time.
In conclusion, our findings are supportive of the fact that different subunits of GABAA receptor subtypes may work together in the pathogenesis of ASD, as well as play differential roles in the symptom-based and developmental deficits in Chinese Han children and adolescents with ASD. To date, very few studies have explored the role of these genetic variants of GABRB3, GABRA5, and GABRG3 genes in both ASD susceptibility and phenotypes, replication studies with a larger sample size are required to confirm our findings in different races.
Materials and Methods
Participants
The study protocol was approved by the Tianjin Medical University Institutional Review Board (Ethics Committee of Tianjin Medical University) and written informed consent was obtained from school principals, parents and/or caregivers to allow the collection and genetic analysis of DNA samples after a complete and extensive description of the study, in accordance with the Declaration of Helsinki. All participants enrolled into this study were Chinese Han children and adolescents from Tianjin, China. We recruited 99 ASD patients (79 males and 20 females) from five local special education schools and 231 age- and gender- frequency-matched TD controls (185 males and 46 females) from four local mainstream schools with a male to female ratio of 4:1 between 2 and 18 years of age. The diagnosis of ASD was made by qualified and experienced psychiatrists basing on the criteria of Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition (DSM-IV)51 and was further identified using the Childhood Autism Rating Scales (CARS)52. A CARS total score of ≥ 30 is the cutoff for distinguishing children ‘at risk’ for autism from pervasive developmental disorder-not otherwise specified (PDD-NOS)53, 54, and a CARS total score of ≥ 25.5 is indicative of ASD54. Meanwhile, the Autism Behavior Checklist (ABC)55 was also administered for the assessment of their symptom-based phenotypes and Early Childhood Development Questionnaire (ECDQ) was used to evaluate their developmental phenotypes. In addition, intelligent quotient (IQ) was determined using the Chinese Binet Scale (Binet)56 by clinical pediatric psychologists. ASD subjects were included with a full-scale IQ of ≥ 30 and TD controls were included with an IQ in the average range. Participant characteristics are presented in Table 5. In the present study, none of the participants in any group were taking drugs. Exclusion criteria for both groups included medical history of epileptic seizures, head trauma, and known psychiatric, neurological or genetic disorders, and TD controls were with a score of < 7 on the Clancy Autism Behavior Scale (CABS)57.
Phenotypic Assessment
Childhood Autism Rating Scale (CARS)
The Childhood Autism Rating Scale (CARS) is a 15-item behavior-based clinical evaluation that includes 14 domains plus one category of general impressions rated on a four-point scale (1 = appropriate for age; 2 = mildly abnormal; 3 = moderately abnormal; 4 = severely abnormal) according to interaction and observation. The 15 items are as follows: ‘relating to people’, ‘imitative behavior’, ‘emotional response’, ‘body use’, ‘object use’, ‘adaptation to change’, ‘visual response’, ‘listening response’, ‘perceptive response’, ‘fear or anxiety’, ‘verbal communication’, ‘non-verbal communication’, ‘activity level’, ‘level and consistency of intellective relations’, and ‘general impressions’55. The CARS has been proved to have a high degree of internal consistency, inter-rater and test-retest reliability, good discriminant validity, and high criterion-related validity58. It is widely used by psychiatrists to identify children with autism and as a further measurement of the severity of this disease. Total scores can range from a low of 15 to a high of 60; scores of less than 30 indicate that the individual is in the non-autistic range, scores between 30 and 36.5 indicate mild/moderate autism, and scores from 37 to 60 indicate severe autism53. In this study, all ASD patients got this assessment.
Autism Behavior Checklist (ABC)
The Autism Behavior Checklist (ABC) is a behavior checklist that consists of 57 items in 5 categories: ‘sensory (S) (9 items)’, ‘relating (R) (12 items)’, ‘body and object use (B) (12 items)’, ‘language (L) (13 items)’, and ‘social and self-help (s′) (11 items)’. Each item corresponds to a single score referring to a single symptomatological area55. For each item, we set the ratings as a two-point scale (0 = no; 1 = yes). The scale utilizes an observer’s rating of a series of typical autistic behaviors in a certain subject and provides advice for educational intervention. In this study, all ASD patients got this assessment.
Early Childhood Development Questionnaire (ECDQ)
The initial ages of early childhood development were assessed by parents or caregivers according to the normal range of age in different developmental parts of the ECDQ. It contains 9 items, including ‘smiling to his/her mother’, ‘grasping things by himself/herself’, ‘sitting by himself/herself’, ‘walking by himself/herself’, ‘calling daddy or mummy’, ‘speaking phrases’, ‘controlling defecate and urinate’, ‘stopping wetting the bed’, ‘wearing by himself/herself’. The developmental stage was evaluated on the ratings of a three-point scale (1: early, 2: normal, 3: late) for each participant. 67 parents or caregivers of ASD patients completed this assessment.
Sample Collection
Participants were previously informed to fast in the morning before sampling to avoid the influence of food and/or drink. DNA samples were obtained by an experienced technician by taking a sterile swab and rubbing it against the inside of participants’ cheeks for 1 min and placing it into a labeled sterile 2-ml Falcon tube. This was done for both cheeks.
DNA Extraction
DNA was extracted from the swabs using the Swab Gen DNA Kit (CW0530, CWBIO, Beijing, China) according to the manufacturer’s suggestion. The DNA solution for the current experiment was then stored at 4 °C, and the redundant stock was stored at −20 °C. The concentration of each DNA sample was determined using a NanoDrop® ND-2000 spectrophotometer (Thermo Scientific).
SNP Selection and Genotyping
Through the SNP database (http://www.ncbi.nlm.nih.gov/snp/), we selected four SNPs of the GABAA receptor genes with relatively high allele frequency and heterozygosity in Asian populations, including two GABRB3 SNPs (rs2081648, Assay ID: C___2911917_10; rs1426217, Assay ID: C___2901088_10), one GABRA5 SNP (rs35586628, Assay ID: C____252720_10), as well as one GABRG3 SNP (rs208129, Assay ID: C___2665692_10).
All SNPs were genotyped by TaqMan® genotyping assay. The TaqMan probes were ordered from the Assays on Demand system of the Applied Biosystems (Applied Biosystems, Foster City, CA, USA). Genotyping was performed in 96-well plates in 5-μl system containing 2.5 μl of TaqMan® Genotyping Master Mix (Applied Biosystems, Foster City, CA, USA), 0.125 μl of 40 × TaqMan probe (Applied Biosystems, Foster City, CA, USA), 1.375 μl ddH2O and 1 μl (5–20ng) of genomic DNA using Roche 480 PCR System (Roche Applied Science, Penzberg, Germany) in accordance with the manufacturer’s instructions. The first row of each 96-deep-well plate contained the negative controls, which consisted of 2.5 μl of TaqMan® Genotyping Master Mix (Applied Biosystems, Foster City, CA, USA), 0.125 μl of 40 × TaqMan probe (Applied Biosystems, Foster City, CA, USA) and 2.375 μl ddH2O. PCR parameters were as follows: enzyme activation: 95 °C for 10 min, 45 cycles of amplification: 95 °C for 15 s and 60 °C for 1 min and cooling at 40 °C for 10 s. Results from the amplified PCR products were viewed using the Roche LightCycler®480 II Real-Time PCR System (Roche Applied Science, Penzberg, Germany). Analysis of the SNP genotypes was performed using a Roche LightCycler®480 Sequence Detection System (Roche Diagnostics GmbH, Penzberg, Germany).
Statistical Analysis
In both ASD patients and TD controls, Haploview version 4.259, 60 was conducted to assess the Hardy-Weinberg equilibrium for genotypic frequency distributions by χ 2 test and the linkage disequilibrium (LD) analysis among SNPs by D′. A value of D′ > 0.5 was defined as a positive linkage and D′ > 0.75 were considered as a strong linkage. Power and Sample Size Calculation (http://biostat.mc.vanderbilt.edu/wiki/Main/PowerSampleSize) was employed to perform a power analysis. Haplotype counts and frequencies in both ASD and TD groups were evaluated using SHEsis (http://analysis.bio-x.cn/SHEsisMain.htm)61, a frequency of < 0.03 in both groups was ignored in the association analysis. The other analyses were conducted by SPSS version 17.0 (IBM Corporation New York, USA) and SAS version 9.3.2 (SAS Institute Inc, Cary, NC). Comparisons of genotypic and allelic frequencies between ASD patients and TD controls were performed with the χ 2 test. The Kolmogorov-Smirnov test was performed to check the distribution of our data. The Levene’s test was utilized to determine the homogeneity of variance. Kruskal-Wallis H tests were utilized for comparisons of CARS total scores among genotypes of four SNPs. A one-way analysis of covariance (ANCOVA) with sex, age, and IQ as covariates was applied for comparisons of ABC scores among genotypes of four SNPs. While adjusting for sex, age, and IQ, an ordinal polytomous logistic regression analysis was performed to discover the associations between ratings in items of CARS as well as ECDQ and genotypes of four SNPs, and a binary logistic regression analysis was utilized to find the associations between ratings in items of ABC and genotypes of four SNPs in ASD patients. Odds ratios (ORs) were expressed with 95% confidence intervals (95% CI). Post-hoc Bonferroni correction for multiple comparisons was performed. A p value of < 0.05 was described as significant.
Data Availability
The datasets generated and analyzed during the current study are not publicly available due to the public availability would violate the privacy of participants, especially for the special individuals with ASD. Data are from the Department of Maternal, Child and Adolescent Health, School of Public Health, Tianjin Medical University and are available from the corresponding author (email: zhangxin@tmu.edu.cn) on a reasonable request.
References
American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders: DSM-V. Arlington, TX: American Psychiatric Association (2013).
Brian, J., Doyle-Thomas, K., Baribeau, D. & Anagnostou, E. Novel mechanisms and treatment approaches in autism spectrum disorder. Discov. Med 22, 47–54 (2016).
CDC. Prevalence and Characteristics of Autism Spectrum Disorder Among Children Aged 8 Years–Autism and Developmental Disabilities Monitoring Network, 11 Sites, United States, 2012. MMWR Surveill Summ. 65, 1–23 (2016).
Fatemi, S. H. The hyperglutamatergic hypothesis of autism. Prog Neuropsychopharmacol Biol. Psychiatry 32, 911, author reply 912–913 (2008).
Pizzarelli, R. & Cherubini, E. Alterations of GABAergic signaling in autism spectrum disorders. Neural. Plast. 2011, 297153 (2011).
Yizhar, O. et al. Neocortical excitation/inhibition balance in information processing and social dysfunction. Nature 477, 171–178 (2011).
Cochran, D. M. et al. Relationship among Glutamine, gamma-Aminobutyric Acid, and Social Cognition in Autism Spectrum Disorders. J. Child. Adolesc. Psychopharmacol. 25, 314–322 (2015).
Omura, Y. et al. Early Detection of Autism (ASD) by a Non-invasive Quick Measurement of Markedly Reduced Acetylcholine & DHEA and Increased beta-Amyloid (1–42), Asbestos (Chrysotile), Titanium Dioxide, Al, Hg & often Coexisting Virus Infections (CMV, HPV 16 and 18), Bacterial Infections etc. in the Brain and Corresponding Safe Individualized Effective Treatment. Acupunct. Electrother. Res. 40, 157–187 (2015).
Schur, R. R. et al. Brain GABA levels across psychiatric disorders: A systematic literature review and meta-analysis of (1) H-MRS studies. Hum. Brain Mapp. 37, 3337–3352 (2016).
Fatemi, S. H. et al. Downregulation of GABAA receptor protein subunits alpha6, beta2, delta, epsilon, gamma2, theta, and rho2 in superior frontal cortex of subjects with autism. J. Autism Dev. Disord. 44, 1833–1845 (2014).
Whissell, P. D. et al. Acutely increasing deltaGABA(A) receptor activity impairs memory and inhibits synaptic plasticity in the hippocampus. Front Neural Circuits 7, 146 (2013).
Reynolds, D. S. et al. Sedation and anesthesia mediated by distinct GABA(A) receptor isoforms. J. Neurosci. 23, 8608–8617 (2003).
Rudolph, U. & Mohler, H. GABAA receptor subtypes: Therapeutic potential in Down syndrome, affective disorders, schizophrenia, and autism. Annu. Rev. Pharmacol. Toxicol. 54, 483–507 (2014).
Veenstra-Vanderweele, J., Cook, E. Jr. & Lombroso, P. J. Genetics of childhood disorders: XLVI. Autism, part 5: genetics of autism. J. Am. Acad. Child. Adolesc. Psychiatry 42, 116–118 (2003).
Ashley-Koch, A. E. et al. An analysis paradigm for investigating multi-locus effects in complex disease: examination of three GABA receptor subunit genes on 15q11-q13 as risk factors for autistic disorder. Ann. Hum. Genet. 70, 281–292 (2006).
Wolpert, C. M. et al. Three probands with autistic disorder and isodicentric chromosome 15. Am. J. Med. Genet. 96, 365–372 (2000).
Schroer, R. J. et al. Autism and maternally derived aberrations of chromosome 15q. Am. J. Med. Genet. 76, 327–336 (1998).
Zafeiriou, D. I., Ververi, A., Dafoulis, V., Kalyva, E. & Vargiami, E. Autism spectrum disorders: the quest for genetic syndromes. Am. J. Med. Genet. B Neuropsychiatr. Genet 162B, 327–366 (2013).
Liu, J. et al. A genomewide screen for autism susceptibility loci. Am. J. Hum. Genet. 69, 327–340 (2001).
Shao, Y. et al. Genomic screen and follow-up analysis for autistic disorder. Am. J. Med. Genet. 114, 99–105 (2002).
Shao, Y. et al. Fine mapping of autistic disorder to chromosome 15q11-q13 by use of phenotypic subtypes. Am. J. Hum. Genet. 72, 539–548 (2003).
Philippe, A. et al. Genome-wide scan for autism susceptibility genes. Paris Autism Research International Sibpair Study. Hum. Mol. Genet. 8, 805–812 (1999).
Tsang, K. M. et al. A genome-wide survey of transgenerational genetic effects in autism. PLoS One 8, e76978 (2013).
Cook, E. H. Jr. et al. Linkage-disequilibrium mapping of autistic disorder, with 15q11-13 markers. Am. J. Hum. Genet. 62, 1077–1083 (1998).
Martin, E. R. et al. Analysis of linkage disequilibrium in gamma-aminobutyric acid receptor subunit genes in autistic disorder. Am. J. Med. Genet. 96, 43–48 (2000).
Buxbaum, J. D. et al. Association between a GABRB3 polymorphism and autism. Mol Psychiatry 7, 311–316 (2002).
McCauley, J. L. et al. A linkage disequilibrium map of the 1-Mb 15q12 GABA(A) receptor subunit cluster and association to autism. Am. J. Med. Genet. B Neuropsychiatr. Genet 131B, 51–59 (2004).
Curran, S. et al. An association analysis of candidate genes on chromosome 15 q11-13 and autism spectrum disorder. Mol. Psychiatry 11, 709–713 (2006).
Kim, S. A., Kim, J. H., Park, M., Cho, I. H. & Yoo, H. J. Association of GABRB3 polymorphisms with autism spectrum disorders in Korean trios. Neuropsychobiology 54, 160–165 (2006).
Maestrini, E. et al. Serotonin transporter (5-HTT) and gamma-aminobutyric acid receptor subunit beta3 (GABRB3) gene polymorphisms are not associated with autism in the IMGSA families. The International Molecular Genetic Study of Autism Consortium. Am. J. Med. Genet. 88, 492–496 (1999).
Salmon, B. et al. Absence of linkage and linkage disequilibrium to chromosome 15q11-q13 markers in 139 multiplex families with autism. Am. J. Med. Genet. 88, 551–556 (1999).
Ma, D. Q. et al. Identification of significant association and gene-gene interaction of GABA receptor subunit genes in autism. Am. J. Hum. Genet. 77, 377–388 (2005).
Nurmi, E. L. et al. Linkage disequilibrium at the Angelman syndrome gene UBE3A in autism families. Genomics 77, 105–113 (2001).
Tochigi, M. et al. No evidence for significant association between GABA receptor genes in chromosome 15q11-q13 and autism in a Japanese population. J. Hum. Genet. 52, 985–989 (2007).
Homanics, G. E. et al. Mice devoid of gamma-aminobutyrate type A receptor beta3 subunit have epilepsy, cleft palate, and hypersensitive behavior. Proc. Natl. Acad. Sci. USA 94, 4143–4148 (1997).
DeLorey, T. M. et al. Mice lacking the beta3 subunit of the GABAA receptor have the epilepsy phenotype and many of the behavioral characteristics of Angelman syndrome. J. Neurosci. 18, 8505–8514 (1998).
Krasowski, M. D., Rick, C. E., Harrison, N. L., Firestone, L. L. & Homanics, G. E. A deficit of functional GABA(A) receptors in neurons of beta 3 subunit knockout mice. Neurosci. Lett. 240, 81–84 (1998).
Huntsman, M. M., Porcello, D. M., Homanics, G. E., DeLorey, T. M. & Huguenard, J. R. Reciprocal inhibitory connections and network synchrony in the mammalian thalamus. Science 283, 541–543 (1999).
DeLorey, T. M., Sahbaie, P., Hashemi, E., Homanics, G. E. & Clark, J. D. Gabrb3 gene deficient mice exhibit impaired social and exploratory behaviors, deficits in non-selective attention and hypoplasia of cerebellar vermal lobules: a potential model of autism spectrum disorder. Behav. Brain Res. 187, 207–220 (2008).
DeLorey, T. M. GABRB3 gene deficient mice: a potential model of autism spectrum disorder. Int. Rev.Neurobiol. 71, 359–382 (2005).
Nurmi, E. L. et al. Exploratory subsetting of autism families based on savant skills improves evidence of genetic linkage to 15q11-q13. J. Am. Acad. Child. Adolesc. Psychiatry 42, 856–863 (2003).
Ma, D. Q. et al. Ordered-subset analysis of savant skills in autism for 15q11-q13. Am. J. Med. Genet. B Neuropsychiatr. Genet 135b, 38–41 (2005).
Fatemi, S. H. et al. mRNA and protein levels for GABAAalpha4, alpha5, beta1 and GABABR1 receptors are altered in brains from subjects with autism. J. Autism Dev. Disord. 40, 743–750 (2010).
Menold, M. M. et al. Association analysis of chromosome 15 gabaa receptor subunit genes in autistic disorder. J. Neurogenet. 15, 245–259 (2001).
Hogart, A., Nagarajan, R. P., Patzel, K. A., Yasui, D. H. & Lasalle, J. M. 15q11-13 GABAA receptor genes are normally biallelically expressed in brain yet are subject to epigenetic dysregulation in autism-spectrum disorders. Hum. Mol. Genet. 16, 691–703 (2007).
Kim, S. J. et al. Transmission disequilibrium testing of the chromosome 15q11-q13 region in autism. Am. J.Med. Genet. B Neuropsychiatr. Genet 147b, 1116–1125 (2008).
Georgiades, S. et al. Structure of the autism symptom phenotype: A proposed multidimensional model. J Am Acad Child Adolesc Psychiatry 46, 188–196 (2007).
Zurek, A. A. et al. alpha5GABAA receptor deficiency causes autism-like behaviors. Ann Clin Transl Neurol 3, 392–398 (2016).
Huang, J. P. et al. Prevalence and early signs of autism spectrum disorder (ASD) among 18–36 month-old children of Tianjin in China. Biomed. Environ.Sci. 27, 453–461 (2014).
Schaefer, G. B. Clinical Genetic Aspects of ASD Spectrum Disorders. Int. J. Mol. Sci. 17 (2016).
American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders: DSM-IV. Washington, DC: American Psychiatric Association (1996).
Schopler, E., Reichler, R. J. & Renner, B. R. The Childhood Autism Rating Scale for Diagnostic Screening and Classification of Autism. Irvington, NY: Irvington (1986).
Schopler, E., Reichler, R. & Rochen Renner, B. The childhood autism rating scale. Los Angeles: Western Psychological Services (1988).
Chlebowski, C., Green, J. A., Barton, M. L. & Fein, D. Using the childhood autism rating scale to diagnose autism spectrum disorders. J. Autism Dev. Disord. 40, 787–799 (2010).
Rellini, E., Tortolani, D., Trillo, S., Carbone, S. & Montecchi, F. Childhood Autism Rating Scale (CARS) and Autism Behavior Checklist (ABC) correspondence and conflicts with DSM-IV criteria in diagnosis of autism. J. Autism Dev. Disord. 34, 703–708 (2004).
Wu, T. M. The guide to Chinese Binet Test, 3rd ed (in Chinese). Beijing: Beijing University Press (1982).
Clancy, H., Dugdale, A. & Rendle-Short, J. The diagnosis of infantile autism. Dev. Med. Child. Neurol 11, 432–442 (1969).
Nordin, V., Gillberg, C. & Nyden, A. The Swedish version of the Childhood Autism Rating Scale in a clinical setting. J. Autism Dev. Disord. 28, 69–75 (1998).
Barrett, J. C., Fry, B., Maller, J. & Daly, M. J. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21, 263–265 (2005).
Wigginton, J. E., Cutler, D. J. & Abecasis, G. R. A note on exact tests of Hardy-Weinberg equilibrium. Am. J. Hum. Genet. 76, 887–893 (2005).
Shi, Y. Y. & He, L. SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci. Cell Res. 15, 97–98 (2005).
Acknowledgements
This research was supported by grants from the National Natural Science Foundation of China (No. 81072313). The authors gratefully thank all the subjects, their schools and families for their support and participation.
Author information
Authors and Affiliations
Contributions
X.Z. conceived of the study, participated in its design and wrote the protocol; Y.S., W.D. and Y.H. collected the samples; X.G. and X.D. performed the experiments; L.G. participated in interpretation of the data; S.Y. collected the samples, performed the experiments, undertook the statistical analysis and wrote the first draft of the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Yang, S., Guo, X., Dong, X. et al. GABAA receptor subunit gene polymorphisms predict symptom-based and developmental deficits in Chinese Han children and adolescents with autistic spectrum disorders. Sci Rep 7, 3290 (2017). https://doi.org/10.1038/s41598-017-03666-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-017-03666-0
This article is cited by
-
Gamma-Aminobutyric Acid Type A Receptor Variants are Associated with Autism Spectrum Disorders
Journal of Molecular Neuroscience (2023)
-
SNP-PCR genotyping links alterations in the GABAA receptor (GABRG3: rs208129) and RELN (rs73670) genes to autism spectrum disorder among peadiatric Iraqi Arabs
Molecular Biology Reports (2022)
-
GABA Receptor SNPs and Elevated Plasma GABA Levels Affect the Severity of the Indian ASD Probands
Journal of Molecular Neuroscience (2022)
-
An Association Study of Gamma-Aminobutyric Acid Type A Receptor Variants and Susceptibility to Autism Spectrum Disorders
Journal of Autism and Developmental Disorders (2021)
-
Identification of primary copy number variations reveal enrichment of Calcium, and MAPK pathways sensitizing secondary sites for autism
Egyptian Journal of Medical Human Genetics (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
