Abstract
Keywords
Epilepsy Drugresistance Genetic Polymorphism SCN1A ABCB1 SCN1B MLPA
1. Background
Generally, the clinicians have to use anti-epileptic drugs (AEDs) by trial and error, therefore a broad range of doses are prescribed before determining the final maintenance dose. The SNPs analysis as an informative molecular tool can be used to improve the phenotype-genotype correlation, along with the drug dosage decisions in drug resistant patients with epilepsy.
There are various proteins, including ion channels, drug transporters and receptors which have a critical role in drug resistance (1-5). Some reported SNPs in the SCN1A gene (rs590478, rs8191987, rs3812718, and rs2126152) are associated with epilepsy drug resistance (6, 7).
On the other hand, a well known polymorphism (rs1045642 C > T) in exon 26 of the MDR1 gene that induces the transport of antiepileptic drugs shows clinical associations with drug resistance (6, 8, 9). It is proposed that a multigeneic interaction between the genes involved in both ABCB1 as well as the SCN1A genes could lead to a pharmaco-resistance to the AEDs (6, 7, 10, 11).
We aimed to investigate whether there is a significant difference between the SCN1A and ABCB1 polymorphism’s frequency in epileptic patients and the control groups or not. This would confirm the hypothesis that genetic variation plays a role in determining the inter-individual differences in clinical signs, along with their individual responses to anti-epileptic drugs.
2. Methods
In this study, patients with an Iranian ancestry, were referred by the epilepsy clinic of the Mofid hospital, Shahid Beheshty University, Tehran, Iran, to our laboratory. All the patients were classified according to their histories, clinical examinations, and workups based on the last revised terminology and concepts published by the international league against epilepsy (ILAE) (12). The written informed consent forms were signed by the parents. Extensive clinical data were collected and stored in a database. Para-clinical workups and genetic counseling were done.
34 probands, clinically diagnosed as one type of idiopathic refractory epilepsy (IRE) spectrum, were selected. All patients were tested for their multidrug resistance during their regular treatment of epilepsy. Refractory seizures were defined as “not response to treatment”. Furthermore, patients with IRE did not achieve immunity against seizures after being treated with 3 - 5 commonly used AEDs.
Although some of the patients only showed trivial responses to the broad spectrum of AEDs, such as topiramate, which share multiple mechanisms of antiepileptic actions, the control group that included 81 patients with symptomatic, idiopathic or sporadic form of febrile or afebrile seizures were treated using general AEDs.
DNA was extracted from blood samples using the commercial iNTRON G-spin genomic DNA extraction kit.
2.1. SCN1A SNPs Analysis
We selected two Tag SNPs, including 1067A > G (rs2298771) in exon 16 and intronic SNP, and rs7601520 A > G in intron 15 of SCN1A gene due to its high rate of heterozygosity within the general Iranian population (13). In order to conduct the screening of the SCN1A gene, the SNPs that we designed were intronic primers and amplification refractory mutation system (ARMS) PCR primers by the Primer3 online software. These SNPs were amplified, analyzed by the ARMS-PCR and confirmed by direct sequencing. We used a tetra primer ARMS-PCR for the identification of the two SNPs in the SCN1A gene to evaluate the genotype of the samples (Figure 1).
The SCN1A- Exon 16 SNPs Primer Details
2.2. ABCB1 (MDR1) SNP Analysis
A coding SNP (cSNP), 3435C > T (rs1045642), in the ABCB1 gene was screened by a published PCR-RFLP method (8). To ensure the accuracy of the PCR-RFLP and ARMS-PCR methods, parts of the samples (about 15%) were randomly rechecked by direct sequencing.
The final data association was assessed by statistical software SPSS, where 95% of the confidence interval calculated and P value < 0.05 were considered significant.
3. Results
Genotyping of SCN1A SNPs showed a similar genotype, allelic frequency, and a heterozygosity pattern for the two analyzed SNPs, including rs7601520, rs2298771 compared to the control group and reference heterozygosity adopted from NCBI/SNPs (Tables 1, 2, 3).
Frequency of SCN1A and ABCB1 (MDR1) SNPs in Epileptic and Control Groups
SNPs | Gene | Allele | Control | Patient | Reference* Heterozygosity | HGVS Name |
---|---|---|---|---|---|---|
rs7601520 | SCN1A | AA | 36.81 (0.443) | 12.34 (0.353) | 0.333 | NM_006920.4: c.2914-41G > C |
AG | 25.81 (0.31) | 16.34 (0.47) | NT_005403.16:g.17102499G > A | |||
GG | 20.81 (0.247) | 6.34 (0.177) | ||||
rs2298771 | SCN1A | AA | 37.81 (0.46) | 13.34 (0.38) | 0.328 | NM_006920.4: c.3166G > A |
AG | 24.81 (0.293) | 18.34 (0.54) | NP_008851.3: p. Ala1056Thr | |||
GG | 20.81 (0.247) | 3.34 (0.09) | NT_005403.16:g.17102206C > T | |||
rs1045642 | ABCB1 | CC | 18.81 (0.23) | 5.34 (0.147) | 0.377 | NM_000927.3: c.3435T > C |
CT | 43.81 (0.523) | 16.34 (0.47) | NT_007933.14:g.12372921A > G | |||
TT | 20.81 (0.247) | 13.34 (0.383) |
SCN1A, rs7601520, Population Diversity Obtained from NCBI/SNP
Sample Ascertainment | Genotype Detail | Alleles | ||||||
---|---|---|---|---|---|---|---|---|
Population | Individual Group | Chrom. Sample Cnt. | A/A | A/G | G/G | HWP | A | G |
AFD_EUR_PANEL | European | 48 | 0.208 | 0.625 | 0.167 | 0.251 | 0.521 | 0.479 |
HapMap-CEU | European | 116 | 0.414 | 0.448 | 0.138 | 0.638 | 0.362 | |
European total frequency | 164 | 0.311 | 0.537 | 0.152 | 0.579 | 0.420 | ||
AFD_AFR_PANEL | African American | 44 | 0.773 | 0.182 | 0.045 | 0.294 | 0.864 | 0.136 |
HapMap-YRI | Sub-Saharan African | 118 | 0.763 | 0.220 | 0.017 | 0.873 | 0.127 | |
African American total frequency | 162 | 0.768 | 0.221 | 0.031 | 0.868 | 0.132 | ||
AFD_CHN_PANEL | Asian | 48 | 0.750 | 0.208 | 0.042 | 0.439 | 0.854 | 0.146 |
HapMap-HCB | Asian | 90 | 0.756 | 0.244 | 0.878 | 0.122 | ||
HapMap-JPT | Asian | 88 | 0.750 | 0.227 | 0.023 | 0.864 | 0.136 | |
Asian total frequency | 226 | 0.752 | 0.228 | 0.022 | 0.865 | 0135 | ||
Control group | Iranian | 81 | 0.443 | 0.310 | 0.247 | 0.598 | 0.402 | |
Epileptic patients | Iranian | 34 | 0.353 | 0.470 | 0.177 | 0.588 | 0.412 |
SCN1A, rs2298771, Population Diversity Obtained from NCBI/SNP
Sample Ascertainment | Genotype Detail | Alleles | ||||||
---|---|---|---|---|---|---|---|---|
Population | Individual Group | Chrom. Sample Cnt. | A/A | A/G | G/G | HWP | A | G |
HapMap-CEU | European | 120 | 0.417 | 0.450 | 0.133 | 0.642 | 0.358 | |
AFD_EUR_PANEL | European | 48 | 0.208 | 0.625 | 0.167 | 0.251 | 0.521 | 0.479 |
HapMap-CEU | European | 120 | 0.417 | 0.450 | 0.133 | 1.000 | 0.642 | 0.358 |
European total frequency | 288 | 0.347 | 0.508 | 0.145 | 0.602 | 0.398 | ||
AFD_AFR_PANEL | African American | 44 | 0.773 | 0.182 | 0.045 | 0.294 | 0.864 | 0.136 |
HapMap-YRI | Sub-Saharan African | 120 | 0.767 | 0.217 | 0.017 | 0.875 | 0.125 | |
HapMap-YRI | Sub-Saharan African | 120 | 0.767 | 0.217 | 0.017 | 1.000 | 0.875 | 0.125 |
African American total frequency | 284 | 0.768 | 0.205 | 0.027 | 0.872 | 0.388 | ||
AFD_CHN_PANEL | Asian | 48 | 0.750 | 0.208 | 0.042 | 0.439 | 0.854 | 0.146 |
HapMap-HCB | Asian | 90 | 0.756 | 0.244 | 0.878 | 0.122 | ||
HapMap-JPT | Asian | 90 | 0.733 | 0.222 | 0.044 | 0.844 | 0.156 | |
HapMap-JPT | Asian | 88 | 0.750 | 0.227 | 0.023 | 1.000 | 0.864 | 0.136 |
HapMap-HCB | Asian | 90 | 0.756 | 0.244 | 0.371 | 0.878 | 0.122 | |
Asian total frequency | 406 | 0.749 | 0.229 | 0.022 | 0.229 | 0.136 | ||
Control group | Iranian | 81 | 0.46 | 0.293 | 0.247 | 0.606 | 0.394 | |
Epileptic patients | Iranian | 34 | 0.38 | 0.540 | 0.090 | 0.648 | 0.352 |
The genotype frequencies for rs7601520 were 0.443/0.31/0.247 in the control groups and 0.353/0.47/0.177 in patients, respectively for AA/AG/GG. This SNP shows 0.47 of heterozygosity in the affected probands compared to the 0.333 of the reference heterozygosity.
The second SNP (rs2298771) analysis showed a genotype frequency of 0.46/0.293/0.247 in the control groups and 0.38/0.54/0.09 in patients with a heterozygosity of 0.54 compared to the 0.328 of the reference heterozygosity. The allelic frequency of rs7601520 polymorphism was A; 0.588, G; 0.412 in epileptic patients and A; 0.598, G; 0.402 in the control groups.
The allelic frequencies for rs2298771 were A; 0.606, G; 0.394 in the control groups and A; 0.65, G; 0.35 in the patients.
The ABCB1 polymorphism genotyping showed a frequency of 0.23/0.523/0.247 in the control groups and 0.147/0.47/0.383 in the probands, respectively for CC/CT/TT. Furthermore, the allelic frequency was C; 0.492, T; 0.508 in the control groups compared to C; 0.382, T; 0.618 in the patients (Table 4).
ABCB1, rs1045642, Population Diversity Obtained from NCBI/SNPa
Sample Ascertainment; rs1045642; c.3435T > C | Genotype Detail | Alleles | ||||||
---|---|---|---|---|---|---|---|---|
Population | Individual Group | Chrom. Sample Cnt. | C/C | C/T | T/T | HWP | C | T |
AFD_EUR_PANEL | European | 48 | 0.125 | 0.500 | 0.375 | 0.752 | 0.375 | 0.625 |
Autosome | European | 94 | 0.149 | 0.766 | 0.085 | 0.001 | 0.532 | 0.468 |
HapMap-CEU | European | 116 | 0.155 | 0.603 | 0.241 | 0.457 | 0.543 | |
EGP_CEPH-PANEL | European | 44 | 0.091 | 0.727 | 0.182 | 0.050 | 0.455 | 0.545 |
HapMap-CEU | European | 118 | 0.153 | 0.627 | 0.220 | 0.466 | 0.534 | |
HapMap-CEU | European | 120 | 0.150 | 0.617 | 0.233 | 0.458 | 0.542 | |
European total frequency | 0.137 | 0.640 | 0.223 | 0.457 | 0.553 | |||
AFD_AFR_PANEL | African American | 46 | 0.696 | 0.304 | 0.403 | 0.848 | 0.152 | |
HapMap-YRI | Sub-Saharan African | 118 | 0.797 | 0.186 | 0.017 | 0.890 | 0.110 | |
EGP_YORUB-PANEL | Sub-Saharan African | 24 | 0.583 | 0.417 | 0.371 | 0.792 | 0.208 | |
EGP_AD-PANEL | African American | 28 | 0.857 | 0.143 | 1.000 | 0.929 | 0.071 | |
HapMap-YRI | Sub-Saharan African | 120 | 0.783 | 0.200 | 0.017 | 0.883 | 0.117 | |
HapMap-YRI | Sub-Saharan African | 120 | 0.783 | 0.200 | 0.017 | 0.883 | 0.117 | |
African origins mix | 456 | 0.751 | 0.241 | 0.008 | 0.871 | 0.129 | ||
HapMap-HCB | Asian | 90 | 0.400 | 0.400 | 0.200 | 0.600 | 0.400 | |
HapMap-JPT | Asian | 90 | 0.244 | 0.556 | 0.200 | 0.522 | 0.478 | |
HapMap-HCB | Asian | 90 | 0.400 | 0.400 | 0.200 | 0.600 | 0.400 | |
HapMap-JPT | Asian | 88 | 0.250 | 0.545 | 0.205 | 0.523 | 0.477 | |
EGP_ASIAN-PANEL | Asian | 44 | 0.273 | 0.409 | 0.318 | 0.403 | 0.477 | 0.523 |
HapMap-HCB | Asian | 90 | 0.400 | 0.400 | 0.200 | 0.600 | 0.400 | |
HapMap-JPT | Asian | 90 | 0.244 | 0.556 | 0.200 | 0.522 | 0.478 | |
AFD_CHN_PANEL | Asian | 48 | 0.292 | 0.667 | 0.042 | 0.050 | 0.625 | 0.375 |
Asian origins total frequency | 630 | 0.317 | 0.516 | 0.167 | 0.575 | 0.425 | ||
Control group | Iranian | 81 | 0.230 | 0.523 | 0.247 | 0.492 | 0.508 | |
Epileptic patients | Iranian | 34 | 0.147 | 0.470 | 0.383 | 0.382 | 0.618 |
The rs1045642 heterozygosity in patients was 0.47 compared to the 0.377 of the reference heterozygosity.
4. Discussion
In the present study, the Three Tag SNPs from the SCN1A and ABCB1 genes were analyzed to determine whether each SNP was associated with AED resistance or not. It is believed that the impairments of the voltage-gated sodium channels may be a potential mechanism for the general resistance to the AED therapy. The published studies provide apparent conflicting evidence for the association between the ABCB1 & SCN1A genes and AED resistance (14, 15).
The statistical analysis showed no significant differences in any of the frequencies; neither could the logistic analysis detect any significant SNP–SNP interaction in the epileptic probands compared to the control groups. Comparing the studied SNPs, one by one, with the other previously reported normal populations, a genotype similarity between the Iranian and European samples was found, while distinct differences between our samples and the Asian populations further confirm ethnicity variations (See Tables 1 - 4)
These results do not rule out the contribution of the voltage-dependent sodium channels in the pathogenesis of AED resistance under any circumstances. It could be argued that the number and characteristics of the variables (SNPs tested) were not sufficient to support an intrinsic value in dissecting the gene-to-gene interactions underlying the AED resistance. The tested SNPs in this study were selected based on their theoretical functional value. It is possible that upon the selection of some coding SNPs with a stronger functional effect, different results could have been obtained.
The comparison of the SNPs data (NCBI/SNP) in the ABCB1 (rs1045642) and SCN1A (rs2298771) genes shows that not only the Iranian epileptic patients but also the normal control groups have similarities in their frequencies with the European population.
Interestingly, comparing these data to other populations did not bring out any similarity between the Iranian population and the Asian or African-American origins, pointing towards the fact that it was unexpected. (See Tables 2 - 4)
These results confirmed the previously published haplotype data which focused on the origin of the Iranian population, highlighted the similarities between the Iranian population and European-Caucasian population (16).
The results could be an evidence for the fact that the AED resistance has a complex mechanism, including cellular pathways responsible for drug transportation, receptors, drug targets, and ion channels (17, 18).
On the other hand, current data confirm that we have a unique mixture of gene pool, suggesting similar population based programs be conducted in order to clarify the Iranian haplotypes in fields of molecular medical research, especially pharmacogenetics.
This study did not demonstrate a significant association between polymorphism and drug-resistance epilepsy of the candidates with the SCN1A and ABCB1 genes. However, we think that the variation of SNPs in the SCN1A and the ABCB1 genes could affect the susceptibility of individuals to the drug-resistant forms of epilepsy.
Acknowledgements
References
-
1.
Remy S, Beck H. Molecular and cellular mechanisms of pharmacoresistance in epilepsy. Brain. 2006;129(Pt 1):18-35. [PubMed ID: 16317026]. https://doi.org/10.1093/brain/awh682.
-
2.
Ragsdale DS, Avoli M. Sodium channels as molecular targets for antiepileptic drugs. Brain Res Brain Res Rev. 1998;26(1):16-28. [PubMed ID: 9600622].
-
3.
Abe T, Seo T, Ishitsu T, Nakagawa T, Hori M, Nakagawa K. Association between SCN1A polymorphism and carbamazepine-resistant epilepsy. Br J Clin Pharmacol. 2008;66(2):304-7. [PubMed ID: 18489610]. https://doi.org/10.1111/j.1365-2125.2008.03203.x.
-
4.
Lakhan R, Kumari R, Misra UK, Kalita J, Pradhan S, Mittal B. Differential role of sodium channels SCN1A and SCN2A gene polymorphisms with epilepsy and multiple drug resistance in the north Indian population. Br J Clin Pharmacol. 2009;68(2):214-20. [PubMed ID: 19694741]. https://doi.org/10.1111/j.1365-2125.2009.03437.x.
-
5.
Ohmori I, Kahlig KM, Rhodes TH, Wang DW, George AJ. Nonfunctional SCN1A is common in severe myoclonic epilepsy of infancy. Epilepsia. 2006;47(10):1636-42. [PubMed ID: 17054685]. https://doi.org/10.1111/j.1528-1167.2006.00643.x.
-
6.
Tate SK, Depondt C, Sisodiya SM, Cavalleri GL, Schorge S, Soranzo N, et al. Genetic predictors of the maximum doses patients receive during clinical use of the anti-epileptic drugs carbamazepine and phenytoin. Proc Natl Acad Sci U S A. 2005;102(15):5507-12. [PubMed ID: 15805193]. https://doi.org/10.1073/pnas.0407346102.
-
7.
Loscher W, Klotz U, Zimprich F, Schmidt D. The clinical impact of pharmacogenetics on the treatment of epilepsy. Epilepsia. 2009;50(1):1-23. [PubMed ID: 18627414]. https://doi.org/10.1111/j.1528-1167.2008.01716.x.
-
8.
Vahab SA, Sen S, Ravindran N, Mony S, Mathew A, Vijayan N, et al. Analysis of genotype and haplotype effects of ABCB1 (MDR1) polymorphisms in the risk of medically refractory epilepsy in an Indian population. Drug Metab Pharmacokinet. 2009;24(3):255-60. [PubMed ID: 19571437].
-
9.
Hoffmeyer S, Burk O, von Richter O, Arnold HP, Brockmoller J, Johne A, et al. Functional polymorphisms of the human multidrug-resistance gene: multiple sequence variations and correlation of one allele with P-glycoprotein expression and activity in vivo. Proc Natl Acad Sci U S A. 2000;97(7):3473-8. [PubMed ID: 10716719]. https://doi.org/10.1073/pnas.050585397.
-
10.
Iannetti P, Parisi P, Spalice A, Ruggieri M, Zara F. Addition of verapamil in the treatment of severe myoclonic epilepsy in infancy. Epilepsy Res. 2009;85(1):89-95. [PubMed ID: 19303743]. https://doi.org/10.1016/j.eplepsyres.2009.02.014.
-
11.
Heinzen EL, Yoon W, Tate SK, Sen A, Wood NW, Sisodiya SM, et al. Nova2 interacts with a cis-acting polymorphism to influence the proportions of drug-responsive splice variants of SCN1A. Am J Hum Genet. 2007;80(5):876-83. [PubMed ID: 17436242]. https://doi.org/10.1086/516650.
-
12.
Engel JJ. Report of the ILAE classification core group. Epilepsia. 2006;47(9):1558-68. [PubMed ID: 16981873]. https://doi.org/10.1111/j.1528-1167.2006.00215.x.
-
13.
Ebrahimi A, Houshmand M, Tonekaboni SH, Fallah Mahboob Passand MS, Zainali S, Moghadasi M. Two novel mutations in SCN1A gene in Iranian patients with epilepsy. Arch Med Res. 2010;41(3):207-14. [PubMed ID: 20682179]. https://doi.org/10.1016/j.arcmed.2010.04.007.
-
14.
Depondt C. The potential of pharmacogenetics in the treatment of epilepsy. Eur J Paediatr Neurol. 2006;10(2):57-65. [PubMed ID: 16531088]. https://doi.org/10.1016/j.ejpn.2005.11.009.
-
15.
Siddiqui A, Kerb R, Weale ME, Brinkmann U, Smith A, Goldstein DB, et al. Association of multidrug resistance in epilepsy with a polymorphism in the drug-transporter gene ABCB1. N Engl J Med. 2003;348(15):1442-8. [PubMed ID: 12686700]. https://doi.org/10.1056/NEJMoa021986.
-
16.
Farjadian S, Ota M, Inoko H, Ghaderi A. The genetic relationship among Iranian ethnic groups: an anthropological view based on HLA class II gene polymorphism. Mol Biol Rep. 2009;36(7):1943-50. [PubMed ID: 18979226]. https://doi.org/10.1007/s11033-008-9403-4.
-
17.
Evans WE, McLeod HL. Pharmacogenomics--drug disposition, drug targets, and side effects. N Engl J Med. 2003;348(6):538-49. [PubMed ID: 12571262]. https://doi.org/10.1056/NEJMra020526.
-
18.
Mann MW, Pons G. Various pharmacogenetic aspects of antiepileptic drug therapy: a review. CNS Drugs. 2007;21(2):143-64. [PubMed ID: 17284096].