ORIGINAL RESEARCH article

Front. Psychiatry, 23 January 2023

Sec. Autism

Volume 13 - 2022 | https://doi.org/10.3389/fpsyt.2022.994166

The association between gene polymorphisms in voltage-gated potassium channels Kv2.1 and Kv4.2 and susceptibility to autism spectrum disorder

  • 1. Department of Children’s and Adolescent Health, Public Health College, Harbin Medical University, Harbin, China

  • 2. Department of Preventive Medicine, School of Public Health, Qiqihar Medical University, Qiqihar, China

  • 3. Faculty of Arts and Science, University of Toronto, Toronto, ON, Canada

Abstract

Background:

Autism spectrum disorder (ASD) is a heritable form of neurodevelopmental disorder that arises through synaptic dysfunction. Given the involvement of voltage-gated potassium (Kv) channels in the regulation of synaptic plasticity, we aimed to explore the relationship between the genetic variants in the KCNB1 and KCND2 genes (encoding Kv2.1 and Kv4.2, respectively) and the risk of developing ASD.

Methods:

A total of 243 patients with ASD and 243 healthy controls were included in the present study. Sixty single nucleotide polymorphisms (SNPs) (35 in KCNB1 and 25 in KCND2) were genotyped using the Sequenom Mass Array.

Results:

There were no significant differences in the distribution of allele frequencies and genotype frequencies in KCNB1 between cases and controls. However, the differences were significant in the allelic distribution of KCND2 rs1990429 (pBonferroni < 0.005) and rs7793864 (pBonferroni < 0.005) between the two groups. KCND2 rs7800545 (pFDR = 0.045) in the dominant model and rs1990429 (pFDR < 0.001) and rs7793864 (pFDR < 0.001) in the over-dominant model were associated with ASD risk. The G/A genotype of rs1990429 in the over-dominant model and the G/A–G/G genotype of rs7800545 in the dominant model were correlated with lower severity in the Autism Diagnostic Interview-Revised (ADI–R) restricted repetitive behavior (RRB) domain.

Conclusion:

Our results provide evidence that KCND2 gene polymorphism is strongly associated with ASD susceptibility and the severity of RRB.

1. Introduction

Autism spectrum disorder (ASD) is a neurodevelopmental disease characterized by impaired social interaction and repetitive/stereotyped behaviors. It is reported that the overall prevalence of ASD is 23.0 per 1,000 children aged 8 years in the United States, and the estimated prevalence is 7.0 per 1,000 children aged 6–12 years in China (, ). This alarming situation implies that ASD is being increasingly recognized as a major health burden. Generally, it is assumed that ASD involves strong genetic components, whose heritability is estimated to be 83% (). Our previous study found that 2,174 candidate genes were closely related to ASD, based on whole-genome sequencing. These genes included 14,310 single nucleotide polymorphisms (SNPs), the majority of which play important roles in the regulation of synaptic development and plasticity (). Numerous studies have confirmed that alterations in synaptic plasticity and neuronal excitation/inhibition imbalance contribute to the major etiology of ASD.

Potassium channels, broadly distributed on neuronal cell membranes, are crucial determinants of synaptic plasticity and neuronal excitability. Accumulating evidence shows that genetic mutations in potassium channels inducing alterations of K+ current change the synaptic plasticity and information processing capacity of the brain, thus potentially impairing brain connectivity and function in the early life of individuals with ASD (). Voltage-gated potassium (Kv) channels are important modulators involved in regulating neuronal excitability, mainly including transient outward potassium channels and delayed rectifier potassium channels. Recently, the strong genetic associations between Kv channels and ASD have been partially explained. The expression level of Kv10.2 was significantly reduced in the hippocampus of VPA-induced rats (a model of ASD), and autism-like behaviors, such as stereotypical behaviors and impaired social and exploratory abilities, were effectively ameliorated after upregulation of Kv10.2 expression by lentivirus injection in the hippocampus (). KCNC1 (encoding Kv3.1) knockout mice show deficits in social interaction and hyperactivity (). The Kv7 mutant has been proposed to be potentially pathogenic in autism, owing to altered action potential generation (). The above-mentioned evidence shows that deleterious mutations in Kv channels impact ASD pathogenesis directly or indirectly, which may cause impairments in social interaction and cognitive function.

Kv2.1 (encoded by KCNB1) is the principal delayed rectifier potassium channel, which alters the action potential threshold and firing frequency (). Kv4.2 (encoded by KCND2) is a major predominant transient outward potassium channel, which regulates neuronal signaling by regulating back-propagating action potentials, synaptic integration, and long-term potentiation (). Kv2.1 and Kv4.2 are expressed at high levels throughout the various hippocampal subfields, especially in the CA1 dendritic field (, ). The hippocampus is an important brain region in ASD, which implies that these two Kv channels might be central nodes of dysfunction in ASD. Research has found that 7 of 19 mutants in KCNB1, including S202F, R306C, R312H, W370R, V378A, P385T, and F416L, altered the activity of the voltage-dependence channel, current density, and conductance, which caused epileptic spasms and autism-like developmental phenotypes (). Zhang et al. () showed that strong gating impairment, associated with substitutions of V404L, or V404M in KCND2, increased susceptibility to autism and epileptic seizures. Numerous studies have shown that ASD and epilepsy often co-exist as parallel syndromes. Even in the absence of epileptic seizures, roughly 4–86% of individuals with autism presented abnormal electroencephalogram patterns (). The frequent comorbidity between ASD and epilepsy implies shared underlying neurological abnormalities. It is widely recognized that KCNB1 and KCND2 are the main causes of developmental and epileptic encephalopathy and neurodevelopmental disability (, ). Therefore, we speculated that KCNB1 and KCND2 may be risk genes for autism.

To investigate whether KCNB1 and KCND2 gene polymorphisms are related to ASD risk in the Chinese Han population, we identified the SNPs in the KCNB1 and KCND2 genes between patients with autism and control individuals by allele frequencies, genotype frequencies, and haplotype analyses. Furthermore, we evaluated the relationship between the SNPs and the severity of ASD symptoms. Our results offer supporting evidence for the involvement of Kv channels in the etiology of ASD.

2. Materials and methods

2.1. Participants

The case-control study included 243 pairs of subjects: the ASD patients were recruited from the Child Development and Behavior Research Center of Harbin Medical University, Harbin, China, and the controls were selected from kindergartens, junior schools in Harbin, China. All the subjects were of Chinese Han ethnicity, aged between 2 and 10 years. The diagnosis was confirmed by two professional clinicians according to DSM–V. Most participants with autism were evaluated using Autism Diagnostic Observation Schedule (ADOS) or Autism Diagnostic Interview-Revised (ADI–R), which are currently the gold standard for assessing ASD severity. Participants were excluded if they had other neuropsychiatric disorders or known genetic disorders, such as epilepsy, tuberous sclerosis, intellectual disability, or fragile X syndrome.

The study received ethics approval from the Institutional Review Board of Harbin Medical University for Medical Sciences (Ethics approval number: HMUIRB2012007). The research was performed in accordance with the Declaration of Helsinki principles. All participants or their relatives provided written informed consent.

2.2. Clinical evaluation

The ADOS and ADI–R were assessed to obtain information about autism-specific behaviors and symptoms. Each patient was scored by trained examiners through structured clinical interview with their patients or guardians. Ultimately, of the 243 cases, 166 (68.3%) individuals completed the ADOS assessment and 162 (66.7%) individuals completed the ADI–R assessment. Sample characteristics are provided in Supplementary Table 1. ADOS is a semi-structured observational assessment, which is organized into four modules to assess participants’ social and communication abilities in a standardized context. The calibrated severity score (CSS) is a standardized score of the relative severity of autism-specific behaviors, less influenced by demographic and developmental factors (). This diagnostic algorithm consists of social affect (SA) and restricted repetitive behavior (RRB) domains, which are consistent with DSM–V (). A detailed description of procedures for deriving the ADOS–CSS, SA–CSS, and RRB–CSS can be found in the original study (, ). ADI–R is an anamnestic interview with ASD parents or caregivers, mainly providing information on early development. The interview encompasses three behavioral domains: social, communication, and RRB. Higher scores indicate greater impairment. The intelligence quotient (IQ) scores were derived from the Wechsler Intelligence Scale for Children or the Peabody Picture Vocabulary Test.

2.3. SNP selection and genotyping

Tag SNPs in KCNB1 and KCND2 genes were selected from the Chinese Han in the Beijing panel of HapMap Project Phase II, including 2,000–bp upstream regions at least. Tag SNP selection was performed using the Tagger program incorporated in Haploview v.4.2 software (Broad Institute of Massachusetts Institute of Technology and Harvard, Cambridge, MA, USA) according to the following criteria: minor allele frequency ≥ 0.05 and pairwise tagging (r2 ≥ 0.8). Finally, 57 tag SNPs (32 in KCNB1 and 25 in KCND2 that captured from more than 98 and 216 initially identified SNPs, respectively) were selected for genotyping along with three additional single nucleotide variants in KCNB1 (rs587777850, rs587777849, and rs587777848) that have been linked to autism phenotypes in the previous study (). Mutation information and location for tag SNPs are shown in Supplementary Tables 2, 3.

All the blood samples were collected in venous blood collection tubes containing EDTA. Genomic DNA was extracted from blood cells by the Qiagen QIAamp DNA Mini kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions. Quality assessment and concentration estimation of DNA were done using gel electrophoresis and NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States). SNP genotyping was done using the MassArray platform (Sequenom, San Diego, CA, United States) with primers listed in Supplementary Tables 4, 5. Genotyping quality control for all 60 SNPs was tested using blinded duplicate genotyping for 30 random samples, with reproducibility of 100%.

2.4. Statistical analyses

Hardy-Weinberg equilibrium (HWE) among controls, Linkage disequilibrium (LD) and haplotype construction were determined with Haploview v.4.2 software. Allelic associations were performed with logistic regression using SPSS v.21.0 software (SPSS Inc., Chicago, IL, United States). The odds ratio (OR) and 95% confidence interval (CI) were calculated in different genetic models using SNPstats software1. Genotype associations were estimated according to genetic models of codominant, dominant, recessive, over-dominant, and log-additive, with further adjustment for gender and age. The Akaike information criterion (AIC) and Bayesian information criterion (BIC) were used for model selection, with the lower AIC and BIC indicating the better models. To avoid false-positive results in multiple comparisons, p-values were obtained by either FDR corrected using the Benjamini-Hochberg procedure or Bonferroni correction. Potential associations between genotype and symptom phenotype were evaluated by one-way analysis of variance with adjustment for age, sex, and IQ. All p-values were two tailed, with the significance level set at α = 0.05.

3. Results

3.1. General demographics

There were 243 individuals (39 girls and 204 boys) in the ASD group and 243 (51 girls and 192 boys) in the control group; the mean ages were 5.19 ± 1.96 and 4.96 ± 0.97 years, respectively. The groups did not differ statistically in age (t = 1.643, p = 0.101) and gender (χ2 = 1.964, p = 0.161).

3.2. Links between SNPs and ASD risk

There were a total of 7 SNPs out of the 35 SNPs that were removed in KCNB1: the rs237454, rs490840, and rs802952 call rates were lower than 90%; rs587777850, rs58777749, and rs58777748 showed no dimorphism among participants; rs1961192 deviated from HWE in the control group (p < 0.01; Table 1). In total, 4 out of 25 SNPs were excluded in KCND2: the rs802359 and rs12673992 call rates were lower than 90%; rs2191736 and rs2189977 deviated from HWE in the control group (p < 0.01; Table 2).

TABLE 1

SNPs IDGeneGenomic position (bp)Genic positionReference alleleaCall rate%MAFbHWEb (P)
1rs1051295KCNB149372368intron variant, utr variant 3 primeA95.30.4420.362
2rs6019774KCNB149375641intronT99.40.0600.405
3rs2426154KCNB149385458intronT97.90.4190.089
4rs4810952KCNB149389638intronT98.80.4420.233
5rs1961192KCNB149389770intronG98.80.192<0.001*
6rs9636516KCNB149390663intronA94.00.4550.354
7rs756529KCNB149394471intronA99.00.4151.000
8rs7348799KCNB149405190intronT99.80.0620.449
9rs6067087KCNB149410131intronA97.90.4910.455
10rs6019820KCNB149410790intronG98.60.3990.937
11rs237454KCNB149413854intronA67.30.3790.086
12rs237459KCNB149415093intronT98.20.4400.394
13rs237477KCNB149440911intronT99.00.3240.856
14rs3787318KCNB149441530intronT99.20.0710.661
15rs742759KCNB149444608intronG99.60.1380.994
16rs237478KCNB149449322intronC99.40.4590.198
17rs13044742KCNB149457703intronA96.30.1200.883
18rs572845KCNB149460077intronG99.80.1280.142
19rs610412KCNB149461665intronA99.40.1620.749
20rs6019855KCNB149465446intronT98.80.2690.390
21rs10485612KCNB149468115intronA98.40.0720.672
22rs802950KCNB149469295intronC99.00.1761.000
23rs490840KCNB149469418intronT31.70.3190.437
24rs802952KCNB149469732intronT32.90.1901.000
25rs653070KCNB149470871intronC98.60.3350.551
26rs4809745KCNB149471334intronG98.40.1040.484
27rs552068KCNB149471461intronA99.20.4380.908
28rs6125656KCNB149474242intronG99.00.0670.658
29rs477135KCNB149475917intronT98.20.1790.012
30rs566604KCNB149478053intronA99.40.1640.588
31rs7269864KCNB149480071intronT99.20.1100.283
32rs553213KCNB149481027intronA97.90.3570.480
33rs587777850KCNB149374425intron variant, missenseC99.00.0021.000
34rs587777849KCNB149374439intron variant, missenseG98.80.0001.000
35rs587777848KCNB149374519intron variant, missenseG98.80.0001.000

Characteristics of KCNB1 single nucleotide polymorphisms (SNPs).

MAF, minor allele frequency; HWE, hardy–weinberg equilibrium. *p < 0.01, adetermined by most frequent allele among controls; bamong controls (n = 243). The bold represent the statistically significant correlations (p < 0.05), or meaningful models or values.

TABLE 2

SNPs IDGeneGenomic position (bp)Genic positionReference alleleaCall rateMAFbHWEb (P)
1rs1990429KCND2120293793intron variant, utr variant 3 primeG99.00.1230.033
2rs7800545KCND2120303057intronA99.00.1040.970
3rs2191736KCND2120343577intronA94.90.200<0.001*
4rs17142666KCND2120350103intronG97.50.3310.991
5rs7793864KCND2120367201intronT98.20.0560.922
6rs7810357KCND2120385502intronG98.40.0711.000
7rs7793037KCND2120431403intronA98.40.4850.732
8rs2192373KCND2120452784intronC99.00.0630.758
9rs802359KCND2120509095intronA63.00.4700.010
10rs802372KCND2120529262intronA98.60.4710.200
11rs1527650KCND2120603522intronT98.60.0380.568
12rs10278347KCND2120619780intronC98.20.4710.307
13rs2402539KCND2120628249intronT97.30.4890.317
14rs6979618KCND2120629117intronA98.60.4830.323
15rs7779895KCND2120629812intronC92.80.4170.578
16rs17142875KCND2120632817intronA99.20.4770.070
17rs4727911KCND2120632910intronG98.20.4980.160
18rs7795646KCND2120683673intronA98.40.4540.227
19rs2896298KCND2120687230intronT98.80.4110.034
20rs1072198KCND2120687295intronT98.80.0690.611
21rs17142891KCND2120698554intronG97.70.4180.061
22rs2189977KCND2120704164intronC99.40.367<0.001*
23rs11983106KCND2120719673intronG99.00.0401.000
24rs12673992KCND2120732769intronA39.50.4741.000
25rs727228KCND2120749594utr variant 3 primeT99.20.4160.394

Characteristics of KCND2 single nucleotide polymorphisms (SNPs).

MAF, minor allele frequency; HWE, hardy–weinberg equilibrium. *p < 0.01; adetermined by most frequent allele among controls; bamong controls (n = 243). The bold represent the statistically significant correlations (p < 0.05), or meaningful models or values.

Tables 3, 4 show the allele frequency distributions in cases and controls. Allele frequencies for KCNB1 were not significantly different between the two groups. Five tag SNPs showed statistically significant differences in KCND2 after adjusting for age and sex (p < 0.05). The rs1990429 A allele (OR = 0.53, 95% CI = 0.34–0.82, p = 0.004), rs7800545 G allele (OR = 0.53, 95% CI = 0.33–0.85, p = 0.009), rs7793864 A allele (OR = 0.23, 95% CI = 0.10–0.53, p = 0.001), rs7810357 A allele (OR = 0.50, 95% CI = 0.28–0.89, p = 0.019), and rs6979618 G allele (OR = 0.74, 95% CI = 0.58–0.96, p = 0.024) were associated with a lower risk of ASD. Of particular concern, rs1990429 and rs7793864 remained positive after Bonferroni correction (pBonferroni < 0.005).

TABLE 3

SNPs IDGeneAlleleControl (243)Case (243)P#Adjusted OR (95% CI)#
1rs1051295KCNB1A258 (55.8%)246 (53.0%)0.3691.00
G204 (44.2%)218 (47.0%)1.13 (0.87–1.46)
2rs6019774KCNB1T453 (94.0%)442 (91.3%)0.1351.00
C29 (6.0%)42 (8.7%)1.46 (0.89–2.38)
3rs2426154KCNB1T273 (58.1%)283 (58.7%)0.7461.00
C197 (41.9%)199 (41.3%)0.96 (0.74–1.24)
4rs4810952KCNB1T269 (55.8%)286 (59.8%)0.1771.00
C213 (44.2%)192 (40.2%)0.84 (0.65–1.08)
5rs9636516KCNB1A243 (54.5%)249 (53.2%)0.7891.00
G203 (45.5%)219 (46.8%)1.04 (0.80–1.35)
6rs756529KCNB1A281 (58.5%)298 (61.8%)0.2531.00
G199 (41.5%)184 (38.2%)0.86 (0.66–1.11)
7rs7348799KCNB1T454 (93.8%)445 (91.6%)0.2001.00
C30 (6.2%)41 (8.4%)1.38 (0.84–2.25)
8rs6067087KCNB1A239 (50.9%)244 (50.6%)0.9991.00
G231 (49.1%)238 (49.4%)1.00 (0.78–1.29)
9rs6019820KCNB1G286 (60.1%)304 (63.1%)0.3381.00
A190 (39.9%)178 (36.9%)0.88 (0.68–1.14)
10rs237459KCNB1T269 (56.0%)268 (56.5%)0.8851.00
C211 (44.0%)206 (43.5%)0.98 (0.76–1.27)
11rs237477KCNB1T326 (67.6%)327 (68.1%)0.9821.00
C156 (32.4%)153 (31.9%)1.00 (0.76–1.31)
12rs3787318KCNB1T446 (92.9%)435 (89.9%)0.1161.00
C34 (7.1%)49 (10.1%)1.45 (0.91–2.29)
13rs742759KCNB1G417 (86.2%)404 (83.5%)0.2941.00
A67 (13.8%)80 (16.5%)1.21 (0.85–1.72)
14rs237478KCNB1C261 (54.1%)262 (54.1%)0.9111.00
T221 (45.9%)222 (45.9%)1.02 (0.79–1.31)
15rs13044742KCNB1A419 (88.0%)408 (88.7%)0.6951.00
T57 (12.0%)52 (11.3%)0.92 (0.62–1.38)
16rs572845KCNB1G424 (87.2%)427 (88.2%)0.5711.00
A62 (12.8%)57 (11.8%)0.89 (0.61–1.32)
17rs610412KCNB1A404 (83.8%)401 (82.9%)0.7951.00
C78 (16.2%)83 (17.1%)1.05 (0.74–1.47)
18rs6019855KCNB1T348 (73.1%)344 (71.1%)0.5461.00
G128 (26.9%)140 (28.9%)1.09 (0.82–1.45)
19rs10485612KCNB1A440 (92.8%)434 (90.0%)0.1401.00
G34 (7.2%)48 (10.0%)1.41 (0.89–2.24)
20rs802950KCNB1C394 (82.4%)408 (84.3%)0.4121.00
A84 (17.6%)76 (15.7%)0.87 (0.62–1.22)
21rs653070KCNB1C315 (66.5%)319 (65.9%)0.9761.00
T159 (33.5%)165 (34.1%)1.00 (0.77–1.32)
22rs4809745KCNB1G430 (89.6%)413 (86.8%)0.2161.00
A50 (10.4%)63 (13.2%)1.28 (0.86–1.91)
23rs552068KCNB1A271 (56.2%)259 (53.7%)0.5641.00
C211 (43.8%)223 (46.3%)1.08 (0.84–1.39)
24rs6125656KCNB1G446 (93.3%)443 (91.5%)0.3111.00
A32 (6.7%)41 (8.5%)1.28 (0.79–2.08)
25rs477135KCNB1T391 (82.1%)377 (78.9%)0.2501.00
A85 (17.9%)101 (21.1%)1.21 (0.88–1.67)
26rs566604KCNB1A403 (83.6%)384 (79.3%)0.1111.00
G79 (16.4%)100 (20.7%)1.31 (0.94–1.81)
27rs7269864KCNB1T427 (89.0%)421 (87.0%)0.4041.00
C53 (11.0%)63 (13.0%)1.18 (0.80–1.75)
28rs553213KCNB1A305 (64.3%)292 (61.1%)0.3831.00
G169 (35.7%)186 (38.9%)1.13 (0.86–1.47)

Distribution of allele frequencies of KCNB1 single nucleotide polymorphisms (SNPs) in cases and controls (n, %).

#P value and OR adjusted by age and sex.

TABLE 4

SNPs IDGeneAlleleControl (243)Case (243)P#Adjusted OR (95% CI)#
1rs1990429KCND2G419 (87.7%)450 (93.0%)0.004*1.00
A59 (12.3%)34 (7.0%)0.53 (0.340.82)
2rs7800545KCND2A432 (89.6%)452 (94.2%)0.0091.00
G50 (10.4%)28 (5.8%)0.53 (0.330.85)
3rs17142666KCND2G314 (66.8%)320 (66.9%)0.9201.00
A156 (33.2%)158 (33.1%)0.99 (0.75–1.29)
4rs7793864KCND2T454 (94.2%)475 (98.5%)0.001*1.00
A28 (5.8%)7 (1.5%)0.23 (0.100.53)
5rs7810357KCND2G442 (92.9%)462 (96.2%)0.0191.00
A34 (7.1%)18 (3.8%)0.50 (0.280.89)
6rs7793037KCND2A247 (51.5%)269 (56.5%)0.0961.00
G233 (48.5%)207 (43.5%)0.81 (0.62–1.04)
7rs2192373KCND2C448 (93.7%)459 (94.8%)0.4001.00
T30 (6.3%)25 (5.2%)0.79 (0.46–1.37)
8rs802372KCND2A252 (52.9%)242 (50.2%)0.4011.00
G224 (47.1%)240 (49.8%)1.12 (0.87–1.44)
9rs1527650KCND2T458 (96.2%)466 (96.7%)0.6291.00
G18 (3.8%)16 (3.3%)0.84 (0.42–1.68)
10rs10278347KCND2C252 (52.9%)226 (47.3%)0.0891.00
A224 (47.1%)252 (52.7%)1.25 (0.97–1.61)
11rs2402539KCND2T242 (51.1%)235 (49.8%)0.6831.00
C232 (48.9%)237 (50.2%)1.06 (0.82–1.36)
12rs6979618KCND2A248 (51.7%)282 (59.0%)0.0241.00
G232 (48.3%)196 (41.0%)0.74 (0.580.96)
13rs7779895KCND2C266 (58.3%)250 (56.1%)0.5391.00
A190 (41.7%)196 (43.9%)1.09 (0.83–1.42)
14rs17142875KCND2A253 (52.5%)252 (52.3%)0.9621.00
G229 (47.5%)230 (47.7%)1.01 (0.78–1.30)
15rs4727911KCND2G236 (50.2%)252 (52.1%)0.5471.00
T234 (49.8%)232 (47.9%)0.93 (0.72–1.19)
16rs7795646KCND2A262 (54.8%)252 (52.7%)0.4981.00
G216 (45.2%)226 (47.3%)1.09 (0.85–1.41)
17rs2896298KCND2T286 (59.1%)273 (57.4%)0.6031.00
C198 (40.9%)203 (42.6%)1.07 (0.83–1.39)
18rs1072198KCND2T445 (93.1%)439 (91.1%)0.1881.00
C33 (6.9%)43 (8.9%)1.38 (0.86–2.21)
19rs17142891KCND2G279 (58.4%)269 (57.0%)0.6731.00
A199 (41.6%)203 (43.0%)1.06 (0.82–1.37)
20rs11983106KCND2G459 (96.0%)464 (95.9%)0.9581.00
T19 (4.0%)20 (4.1%)1.02 (0.54–1.94)
21rs727228KCND2T281 (58.5%)272 (56.2%)0.4311.00
A199 (41.5%)212 (43.8%)1.11 (0.86–1.43)

Distribution of allele frequencies of KCND2 single nucleotide polymorphisms (SNPs) in cases and controls (n, %).

#P value and OR adjusted by age and sex; *pBonferroni<0.005. The bold represent the statistically significant correlations (p < 0.05), or meaningful models or values.

Further, we conducted analyses examining whether ASD risk differed according to SNPs in five genetic inheritance models. As shown in Table 5, the genotype frequency of rs477135 in KCNB1 differed significantly between case and control groups (p < 0.05), but it was statistically insignificant after being adjusted by FDR-based correction. In KCND2, there were statistically significant differences in the genotype frequencies of rs1990429, rs7800545, rs7793864, rs7810357, and rs6979618 between the two groups (p < 0.05). Notably, after the p-values were adjusted by FDR-based correction, rs7800545 (OR = 0.49, 95% CI = 0.29–0.83, pFDR = 0.045) in the dominant model and rs1990429 (OR = 0.36, 95% CI = 0.22–0.60, pFDR < 0.001) and rs7793864 (OR = 0.15, 95% CI = 0.06–0.41, pFDR < 0.001) in the over-dominant model were associated with a reduced risk of ASD (for further details see Supplementary Tables 6, 7).

TABLE 5

SNPs IDGeneModelGenotypeControl (243)Case (243)OR (95 CI)PAICBICPFDR
rs477135KCNB1CodominantT/T167 (70.2%)147 (61.5%)1.000.031662.0682.8
A/T57 (23.9%)83 (34.7%)1.63 (1.09–2.44)
A/A14 (5.9%)9 (3.8%)0.70 (0.29–1.68)
DominantT/T167 (70.2%)147 (61.5%)1.000.059663.4680.0
A/T–A/A71 (29.8%)92 (38.5%)1.44 (0.99–2.12)
RecessiveT/T–A/T224 (94.1%)230 (96.2%)1.000.240665.6682.2
A/A14 (5.9%)9 (3.8%)0.60 (0.26–1.43)
Over-dominantT/T–A/A181 (76%)156 (65.3%)1.000.012660.6677.30.336
A/T57 (23.9%)83 (34.7%)1.67 (1.12–2.49)
Log-additive1.20 (0.87–1.64)0.260665.7682.3
rs1990429KCND2CodominantG/G180 (75.3%)212 (87.6%)1.00<0.0001652.6673.4
G/A59 (24.7%)26 (10.7%)0.37 (0.22–0.61)
A/A0 (0%)4 (1.6%)NA (0.00–NA)
DominantG/G180 (75.3%)212 (87.6%)1.00<0.0001659.3676.0
G/A–A/A59 (24.7%)30 (12.4%)0.42 (0.26–0.68)
RecessiveG/G–G/A239 (100%)238 (98.3%)1.000.022666.9683.6
A/A0 (0%)4 (1.6%)NA (0.00–NA)
Over-dominantG/G–A/A180 (75.3%)216 (89.3%)1.00<0.0001655.1671.8< 0.001*
G/A59 (24.7%)26 (10.7%)0.36 (0.22–0.60)
Log-additive0.52 (0.33–0.81)0.0035663.6680.3
rs7800545KCND2CodominantA/A194 (80.5%)214 (89.2%)1.000.024666.9687.8
G/A44 (18.3%)24 (10%)0.49 (0.28–0.83)
G/G3 (1.2%)2 (0.8%)0.59 (0.10–3.58)
DominantA/A194 (80.5%)214 (89.2%)1.000.0064665.0681.70.045*
G/A–G/G47 (19.5%)26 (10.8%)0.49 (0.29–0.83)
RecessiveA/A–G/A238 (98.8%)238 (99.2%)1.000.640672.2688.9
G/G3 (1.2%)2 (0.8%)0.65 (0.11–3.96)
Over-dominantA/A–G/G197 (81.7%)216 (90%)1.000.0076665.3682.0
G/A44 (18.3%)24 (10%)0.49 (0.29–0.84)
Log-additive0.54 (0.33–0.87)0.0093665.6682.3
rs7793864KCND2CodominantT/T213 (88.4%)235 (97.5%)1.00<0.0001654.7675.5
A/T28 (11.6%)5 (2.1%)0.16 (0.06–0.41)
A/A0 (0%)1 (0.4%)NA (0.00–NA)
DominantT/T213 (88.4%)235 (97.5%)1.00<0.0001656.3673.0
A/T–A/A28 (11.6%)6 (2.5%)0.19 (0.08–0.46)
RecessiveT/T–A/T241 (100%)240 (99.6%)1.000.240672.2688.9
A/A0 (0%)1 (0.4%)NA (0.00–NA)
Over-dominantT/T–A/A213 (88.4%)236 (97.9%)1.00<0.0001653.9670.6< 0.001*
A/T28 (11.6%)5 (2.1%)0.15 (0.06–0.41)
Log-additive0.23 (0.10–0.54)<0.0001659.1675.8
rs7810357KCND2CodominantG/G205 (86.1%)222 (92.5%)1.000.040663.0683.9
G/A32 (13.4%)18 (7.5%)0.51 (0.28–0.94)
A/A1 (0.4%)0 (0%)0.00 (0.00–NA)
DominantG/G205 (86.1%)222 (92.5%)1.000.019662.0678.6
G/A–A/A33 (13.9%)18 (7.5%)0.49 (0.27–0.90)
RecessiveG/G–G/A237 (99.6%)240 (100%)1.000.220665.9682.6
A/A1 (0.4%)0 (0%)0.00 (0.00–NA)
Over-dominantG/G–A/A206 (86.5%)222 (92.5%)1.000.028662.6679.3
G/A32 (13.4%)18 (7.5%)0.51 (0.28–0.94)
Log-additive0.49 (0.27–0.88)0.015661.5678.20.079
rs6979618KCND2CodominantA/A68 (28.3%)82 (34.3%)1.000.058665.8686.7
A/G112 (46.7%)118 (49.4%)0.86 (0.57–1.30)
G/G60 (25%)39 (16.3%)0.54 (0.32–0.91)
DominantA/A68 (28.3%)82 (34.3%)1.000.150667.4684.1
A/G–G/G172 (71.7%)157 (65.7%)0.75 (0.51–1.11)
RecessiveA/A–A/G180 (75%)200 (83.7%)1.000.023664.3681.00.097
G/G60 (25%)39 (16.3%)0.59 (0.38–0.93)
Over-dominantA/A–G/G128 (53.3%)121 (50.6%)1.000.620669.3686.0
A/G112 (46.7%)118 (49.4%)1.10 (0.76–1.57)
Log-additive0.75 (0.58–0.97)0.026664.5681.2

KCNB1 and KCND2 single nucleotide polymorphisms (SNPs) in different genetic models associated with autism spectrum disorder (ASD) risk.

OR, odds ratio; CI, confidence interval; AIC, akaike’ information criterion; BIC, bayesian information criterion; pFDR, FDR corrected p value. *p < 0.05. NA, not applicable. The bold represent the statistically significant correlations (p < 0.05), or meaningful models or values.

3.3. Haplotype analysis

The LD analysis revealed that among 28 tag SNPs in the KCNB1 gene, 20 SNPs were in high LD (r2 ≥ 0.8) and were arranged in seven haplotype blocks with frequencies > 5%. Further details for all blocks are shown in Figure 1A. There were no significant differences in the haplotype frequencies between patients and controls (see Table 6).

FIGURE 1

TABLE 6

BlockNo.HaplotypesCase ratio (243)Control ratio (243)POR (95% CI)
Block 11AATc0.5160.5161.00
2GGTc0.3900.4010.4100.89 (0.68–1.17)
3GACc0.0760.0610.2801.32 (0.79–2.21)
Block 21AGd0.4960.5071.00
2GAd0.3840.3970.5300.92 (0.70–1.20)
3GGd0.1180.0940.1901.32 (0.87–2.00)
Block 31TTGCe0.3920.4031.00
2CTGTe0.3120.3200.6301.08 (0.80–1.45)
3TTGTe0.1390.1360.5201.14 (0.77–1.69)
4TCACe0.0880.0720.1501.45 (0.87–2.41)
5TTAC e0.0600.0660.8901.04 (0.61–1.78)
Block 41GAf0.8180.8371.00
2ACf0.1330.1280.7400.94 (0.63–1.38)
Block 51CCg0.6500.6631.00
2CTg0.1800.1620.3501.18 (0.83–1.68)
3ATg0.1630.1730.4500.87 (0.61–1.25)
Block 61GAh0.5370.5611.00
2GCh0.3390.3340.8901.02 (0.77–1.35)
3ACh0.1250.1050.2401.28 (0.85–1.93)
Block 71GTATAi0.6050.6441.00
2GTATGi0.1880.1800.9800.99 (0.70–1.41)
3GAGCGi0.1140.0930.2301.29 (0.85–1.95)
4AAGTGi0.0750.0680.3301.28 (0.78–2.40)

Distribution of KCNB1 haplotypes in cases and controls adjusted by age and sex (frequency more than 5).

cThe order of SNPs in estimated analysis of haplotypes frequency: rs9636516, rs756529, and rs7348799; dThe order of SNPs in estimated analysis of haplotypes frequency: rs6067087 and rs6019820; eThe order of SNPs in estimated analysis of haplotypes frequency: rs237477, rs3787318, rs742759, and r237478; fThe order of SNPs in estimated analysis of haplotypes frequency: rs572845 and rs610412; gThe order of SNPs in estimated analysis of haplotypes frequency: rs802950 and rs653070; hThe order of SNPs in estimated analysis of haplotypes frequency: rs4809745 and rs552068; iThe order of SNPs in estimated analysis of haplotypes frequency: rs6125656, rs477135, rs566604, rs7269864, and rs553213.

There were four tag SNPs in high LD (r2 ≥ 0.8) arranged in two haplotype blocks with frequencies > 5% among 21 studied SNPs in the KCND2 gene. Haplotype blocks were constructed by rs1990429/rs7800545 and rs17142875/rs4727911 (shown in Figure 1B). We found that the rs1990429 A–rs7800545 G haplotype carried a lower risk of ASD in Block 1. However, this was no longer statistically significant after Bonferroni correction. There were no statistically meaningful differences between the two groups in Block 2 (see Table 7).

TABLE 7

BlockNo.HaplotypesCase ratio (243)Control ratio (243)POR (95% CI)
Block 11GAj0.8980.8700.0261.00
2AGj0.0770.0990.58 (0.36–0.93)
Block 21ATk0.4850.4930.7701.00
2GGk0.4750.4771.04 (0.81–1.33)

Distribution of KCND2 haplotypes in cases and controls adjusted by age and sex (frequency more than 5).

jThe order of SNPs in estimated analysis of haplotypes frequency: rs1990429 and rs7800545; kThe order of SNPs in estimated analysis of haplotypes frequency: rs17142875 and rs4727911. The bold represent the statistically significant correlations (p < 0.05), or meaningful models or values.

3.4. Association between genotype and phenotype in ASD

According to the above positive results, we selected three SNPs (rs1990429, rs7800545, and rs7793864) under the optimal genetic model, to evaluate the correlations between KCND2 genotype and ASD symptom severity, including the ADI–R and ADOS domains. The analysis of variance revealed that the scores on the ADOS did not significantly differ between the two groups. The results showed that the G/A genotype of rs1990429 in the over-dominant model and the G/A–G/G genotype of rs7800545 in the dominant model had lower scores for the ADI–R RRB domain after adjusting for age, sex, and IQ (p = 0.001). However, no significant associations were observed for rs7793864 in the over-dominant model (shown in Tables 8, 9).

TABLE 8

SNPsGenetic modelGenotypeSA–CSSFPRRB–CSSFPADOS–CSSFP
rs1990429Over-dominantG/G–A/A7.13 ± 1.660.3510.5546.44 ± 1.883.3630.0696.82 ± 1.590.0540.817
G/A7.39 ± 1.615.56 ± 2.366.78 ± 1.56
rs7800545DominantA/A7.15 ± 1.670.0800.7776.41 ± 1.891.7850.1836.82 ± 1.590.4100.523
G/A–G/G7.00 ± 1.755.72 ± 2.346.56 ± 1.69
rs7793864Over-dominantT/T–A/A7.13 ± 1.691.7240.1916.35 ± 1.950.4340.5116.79 ± 1.610.5360.465
A/T8.00 ± 0.825.50 ± 1.927.25 ± 1.26

The genotype association of KCND2 with ADOS–CSS adjusted by age, sex, and IQ.

ADOS, autism diagnostic observation schedule; SA, social affect; RRB, restricted repetitive behavior; CSS, calibrated severity scores.

TABLE 9

SNPsGenetic modelGenotypeSOCFPVCFPNVCFPRRBFP
rs1990429Over-dominantG/G–A/A22.88 ± 4.500.0420.83817.33 ± 3.640.02200.88310.98 ± 2.810.6110.4365.83 ± 2.7810.950.001
G/A22.33 ± 6.4217.09 ± 4.3011.50 ± 2.903.55 ± 2.33
rs7800545DominantA/A22.96 ± 4.551.1150.29317.38 ± 3.640.0070.93311.05 ± 2.800.7520.3875.81 ± 2.7912.5670.001
G/A–G/G21.67 ± 6.2217.33 ± 4.3611.61 ± 2.793.50 ± 2.26
rs7793864Over-dominantT/T–A/A22.78 ± 4.680.8510.35817.36 ± 3.710.3920.53311.10 ± 2.800.4100.5235.58 ± 2.840.0310.860
A/T23.25 ± 8.2615.00 ± 4.249.75 ± 3.504.50 ± 2.65

The genotype association of KCND2 with ADI–R adjusted by age, sex, and IQ.

ADI–R, autism diagnostic interview–revised; ADI–R SOC, ADI–R social; ADI–R VC, ADI–R communication for verbal; ADI–R NVC, ADI–R communication for non-verbal; ADI–R RRB, ADI–R restricted repetitive behaviors. The bold represent the statistically significant correlations (p < 0.05), or meaningful models or values.

4. Discussion

In this case-control study, we aimed to assess the linkage of polymorphisms of the KCNB1 and KCND2 genes (including tag SNPs and pathogenic SNPs) with the risk of ASD using allele frequencies, genotype frequencies, and haplotype analyses. We also explored the potential relationship between genetic polymorphism and ASD symptom severity. This study found that KCND2 rs1990429, rs7800545, and rs7793864 were associated with ASD risk, of which rs1990429 and rs7800545 were highly correlated with repetitive stereotyped behavior in ASD.

Our study did not find KCNB1 to be associated with ASD risk. Previous research has found that the V378A variant in KCNB1 influenced the expression and localization of Kv2.1 protein, thereby perturbing the voltage-activated current, the ionic selectivity, and the ability of the channel to repolarize (). There was evidence that mutated KCNB1 caused “autism-like” features, including repetitive behaviors and impulsivity, and seizure susceptibility, by affecting the highly conserved structural function of Kv2.1 (). Individuals with KCNB1 pathogenic variants (such as p.Glu43Gly, p.Arg312His, and p.Trp369Arg) were at high risk of ASD and exhibited impaired communication and socialization (). In particular, the missense variants in KCNB1, including rs587777848 (S347R), rs587777849 (T374I), and rs587777850 (G379R), have been shown to be related to language and motor delays (similar to autistic behaviors). These three missense variants were located in the functionally important pore domain of the Kv2.1 protein and produced a loss-of-function effect, resulting in altered ion selectivity and reduced current density at depolarized voltages (). Regrettably, we did not find that rs587777848, rs587777849, and rs587777850 were polymorphic loci or find any significant associations between mutations in KCNB1 and ASD in the Chinese Han population. There are several possible explanations for the apparent differences between our findings and prior studies. ASD is a disorder of great genetic complexity and heterogeneity, making it difficult to delineate the contribution of any single gene to the risk of this disease, where a single candidate gene might be accountable for various disease phenotypes. In addition, differences in the characteristics of participants may be the main reason for the different results. In previous research, patients had comorbid epilepsy, intellectual disability, or developmental delay, while our patients were confined to those with single-incidence ASD among the Chinese Han population. This implies that mutation in KCNB1 may be an etiological factor shared by epilepsy, intellectual disability, and ASD, rather than the specific pathogenic factor in ASD.

KCND2 is located on chromosome 7 (7q31.31). A clinical case report confirmed that genomic copy number loss in this region mediated the core clinical features of ASD, such as stereotypic movements, impairment of social interaction, and poor social skills (, ). This means that KCND2 is located on the most significant susceptibility locus in autism. A study has shown that Kv4.2 protein expression was downregulated in the hippocampal of Fmr1 knockout mice (a common model of ASD), resulting in excess neuronal excitability (). KCND2 knockout mice exhibited delayed synaptic maturation and hippocampal-dependent learning and memory deficits, which suggested a critical role of Kv4.2 in cognition (, ). In addition, Mikhailov et al. () identified three substituted KCND2 variants (N544S, F538S, and R539L) in patients with ASD, which interfered with the expression of Kv4.2 protein. Recent research has also indicated that KCND2 rs10239799 allele C, which was positively selected as an ASD risk allele, was essential in impacting higher-order brain functions such as cognition, behavior, and memory (). Lee et al. () found that the deleterious de novo variant V404M in KCND2 impaired potassium channel inactivation in monozygotic twins with ASD and epilepsy. In vitro studies found that the substitution variant of V404M is located in S6 transmembrane segment, which surrounds the central ion conduction pathway mediating closed state inactivation (). In the Kv4.2 channel with V404M mutant, the inactivation is enhanced directly from preopen closed states, while the pore closure rate is dramatically slowed when the channel opens compared to Kv4.2 WT (, ). However, the pathogenic mutations in KCND2 (N544S, F538S, R539L, and V404M), which have been confirmed to be present in ASD in previous studies, were neither found in the dbSNP database nor detected by MassArray sequencing in this study. Overall, we and others have reported similar findings that KCND2 was crucial for the occurrence of ASD. This is the first time that we found the G/A genotype of rs1990429, G/A–G/G genotype of rs7800545, and A/T genotype of rs7793864 reduced the risk of ASD. To date, there are lack of studies related to biological functionality of these three gene loci. Further research will be required to explore the effects of rs1990429, rs7800545, and rs7793864 in KCND2 on Kv4.2 channel expression and function, and their impact on ASD at the cellular levels and animal models.

To improve the detection capability, haplotype analyses were conducted to estimate the combined influence of genetic variants on ASD risk. Regrettably, although we derived 7 haplotype blocks across 28 SNPs in KCNB1, they did not play a role in ASD. Of the two haplotype blocks detected in KCND2, the AG (rs1990429rs7800545) was a protective factor for ASD, while this statistical difference disappeared after correction for multiple comparisons. In addition, we used ADOS and ADI–R to evaluate the core symptoms of ASD. We found that patients with the G/A genotype of rs1990429 and the G/A–G/G genotype of rs7800545 had lower severity of stereotyped behaviors. Connolly et al. () pointed out that KCND2 was significantly associated with overly serious facial expressions in patients with ASD, which is part of the communication subscale of the Social Responsiveness Scale. However, there is evidence to confirm that reduced Kv4.2 expression altered dendritic spine morphology and density, but did not induce perseverative or repetitive behavior in mice (). In general, there is a lack of consistent association between KCND2 and ASD phenotype.

This study has several limitations that may influence the interpretation of the study results. First, the best model was selected based on the smallest AIC and BIC, which were estimated by maximum likelihood estimation. Given the principle of parsimony, this approach increases the risk of Type I errors. Moreover, this strategy yielded a specific but not very sensitive rule for model selection. Thus, it is necessary to explore the relationship between the risk genotype and the incidence of ASD in multiple models. Second, the results might be influenced by other confounding factors, such as early exposure factors, physical and chemical factors, and nutritional imbalance. Third, our study focused on the Chinese Han population, so further research should include a multi-ethnic population. Finally, studies of the altered function of mutated genes should be validated in animal models.

5. Conclusion

Given its role in synaptic function and plasticity, we examined the role of Kv channels in the ASD risk. Our results support KCND2 as a susceptibility gene for ASD and illustrate that three SNPs (rs1990429, rs7800545, and rs7793864) in the KCND2 gene reduce the risk of developing ASD, among which rs1990429 and rs7800545 alleviate the severity of RRB. These findings will provide important insights into ASD etiopathogenesis and genetic etiology.

Statements

Data availability statement

The datasets presented in this study can be found in dbSNP 156 database with ID 1063426 (https://www.ncbi.nlm.nih.gov/SNP/snp_viewBatch.cgi?sbid=1063426).

Ethics statement

The studies involving human participants were reviewed and approved by the Institutional Review Board of Harbin Medical University for Medical Sciences. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.

Author contributions

MZ provided financial support for the conduct of the research and preparation of the manuscript. MZ and CS designed the study. WX, PG, and FW performed the clinical assessment and experiments. ZL and XY collected the data. MZ and YC analyzed the data. ZL wrote the manuscript. All authors read, revised, and approved the final manuscript.

Funding

This work was supported by Post-doctoral Scientific Research Development Fund of Heilongjiang Province (Grant Number LBH–Z18216).

Acknowledgments

We deeply thank all patients and their parents for their support in this work.

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.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2022.994166/full#supplementary-material

References

Summary

Keywords

autism spectrum disorder, Kv2.1, Kv4.2, single nucleotide polymorphisms, susceptibility

Citation

Liu Z, Yang X, Guo P, Wang F, Xia W, Chen Y, Zou M and Sun C (2023) The association between gene polymorphisms in voltage-gated potassium channels Kv2.1 and Kv4.2 and susceptibility to autism spectrum disorder. Front. Psychiatry 13:994166. doi: 10.3389/fpsyt.2022.994166

Received

14 July 2022

Accepted

30 December 2022

Published

23 January 2023

Volume

13 - 2022

Edited by

Giulia Treccani, Johannes Gutenberg University Mainz, Germany

Reviewed by

Janith Heendeni Pathiranage Don, Translational Genomics Research Institute (TGen), United States; Angela Caruso, National Institute of Health (ISS), Italy

Updates

Copyright

*Correspondence: Mingyang Zou, Caihong Sun,

This article was submitted to Autism, a section of the journal Frontiers in Psychiatry

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics