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).
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).
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) |
Frequency of SCN1A and ABCB1 (MDR1) SNPs in Epileptic and Control Groups
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, 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 |
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 |
SCN1A, rs2298771, Population Diversity Obtained from NCBI/SNP
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).
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 |
ABCB1, rs1045642, Population Diversity Obtained from NCBI/SNPa
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.