1. Introduction
Epilepsy is a neurological disorder that is described by intermittent fits and seizures [Reference Nazish, Ullah and Ullah1]. Epilepsy was characterized theoretically in 2005 as a disorder of the brain described by repeatedly occurring epileptic seizures [Reference Fisher, Acevedo and Arzimanoglou2].Several antiepileptic drugs have been developed to treat epileptic seizures, but not all the patients respond to the drugs. Approximately 20% of the epileptic patients do not respond to drugs and are considered to be drug resistant and continue to have debilitating seizures [Reference Margari, Legrottaglie and Vincenti3].
Drug-resistant epilepsy is defined as consistency of epileptic seizures despite recommended use of antiepileptic drugs [Reference Kwan, Arzimanoglou and Berg4]. Regardless of the fact that many advances in this field have been made to develop new antiepileptic drugs with improved efficacy, there is no significant decrease in the proportion of patients with drug-resistant epilepsy [Reference Chen, Brodie, Kwan and Kwan5]. The mechanism for the development of drug resistance is not fully understood [Reference Sánchez, Herranz and Leno6]. Multiple hypotheses were proposed to understand the underlying mechanism of drug resistance. The target and transporter hypothesis are two well-understood hypotheses for explaining drug resistance [Reference Kasperavičiūtė7,Reference Rogawski and Johnson8].
Voltage-gated ion channels have an important role in triggering epilepsy, and also, they act as targets for a number of antiepileptic drugs [Reference Kumari, Lakhan, Garg, Kalita and Mittal9]. Voltage-gated sodium channels play a crucial role in generating action potential as well as in membrane excitability. Voltage-gated sodium channels consist of α and ß subunits. Any noticeable change or anomaly in subunits of sodium channels alter their activation; i.e., they are activated at a slower rate. Because of slower activation, the membrane remains depolarized for a longer period of time and can be the cause of epilepsy and generation of epileptic seizures [Reference Alekov, Rahman, Mitrovic, Lerche and Lerche10,Reference Vilin11].
Neuronal voltage-gated sodium channels act as targets for a number of antiepileptic drugs (AEDs), i.e., carbamazepine, phenytoin, and valproate. Aside from their job in nerve conduction and the process of epileptogenesis, these voltage-gated sodium channels are likewise perceived as the significant focuses with respect to AED viability [Reference Rogawski and Löscher12,Reference Xie, Dale, John, Cater, Peakman and Clare13]. The sodium channel genes SCN1A and SCN2A encode the α subunit of voltage-gated sodium channels, and thus, mutations in these genes can be potential source of drug resistance in epilepsy [Reference Lakhan, Kumari, Misra, Kalita, Pradhan, Mittal and Mittal14,Reference Meisler15].
Single-nucleotide polymorphisms (SNPs) are the most common variations in the human genome which cause changes in efficacy, suitability, and duration of drug action, and they can be considered as possible factors for drug-resistant epilepsy [Reference Löscher, Klotz, Schmidt and Schmidt16]. SCN1A and SCN2A are reported to be associated with effectiveness of drug therapy whether it is mono or multiple antiepileptic drug therapy, but the results of different investigations are contradicted. For instance, SCN1A (rs2298771) is found to be significantly associated with response to antiepileptic drugs in epileptic patients [Reference Zhou, Zhou and Yin17], while some studies reported that the SCN1Ars2298771is not linked with responsiveness to drugs, and similar is the case with SCN2A [Reference Haerian, Baum, Kwan, Tan, Mohamed and Mohamed18].
SCN1Ars3812718 G > A SNP is found to be linked with regular and maximum use of antiepileptic drugs carbamazepine and phenytoin [Reference Abe, Seo, Ishitsu, Nakagawa, Nakagawa and Nakagawa19–Reference Tate, Singh and Hung21], but there can be population variations as this association was not found to be significant in the Austrian population [Reference Zimprich, Stogmann and Bonelli22].
Keeping in view the unsure role of SCN1A SNPs, we, for the first time in the Pakistani population, analyzed the association pattern of these SNPs in drug-resistant epilepsy.
2. Materials and Methods
2.1. Subjects/Participants
The study was conducted on 201 individuals, consisting of 101 patients and 100 controls. The patients were further accounted in two groups: 42 were drug resistant, and 59 were drug responders. All the drug-resistant patients were clinically diagnosed by neurophysiologists. The main diagnostic criterion was continuation of seizures after following at least two recommended and appropriately chosen therapeutic regimens. The average age of drug responders was 5.5 ± 2.2 years, while the average age of drug resistant was 6.5 ± 2.4 years. Controls were healthy individuals having no physical illness or any history of neurological disorder. The mean age of controls was 6.1 ± 2.3 years. All the subjects of this study were of the same ethnic region to exclude potential biases.
This study was ethically approved by the Bioethics Committee of the University of Punjab, Lahore, Pakistan.
2.2. Genotyping
Blood samples of all participants were collected in EDTA-coated tubes (Catalogue No.17091619). Genomic DNA was isolated using an organic method [Reference Sambrook and Maniatis23]. Exon 6(730 G > T, 735G > T, and 736 A > T), exon 26(rs794726754, rs1057518252, rs121918809, rs12191792, and rs121917932), rs10167228, and rs2298771 of the SCN1A gene and rs17183814 of SCN2A were selected for mutational analysis. The targeted SNPs were amplified by polymerase chain reaction using reported primers (Table 1). The primers were optimized at 57°C, 59.5°C, 61°C, 58°C, and 58.5°C, respectively. The targeted amplicons were checked by performing gel electrophoresis after amplification by PCR. The cycling conditions of PCR were as follows: denaturation at 95°C for 5 mins followed by denaturation at 95°C for 30 sec, annealing for 45 sec at the optimized temperatures for selected SNP, and extension at 72°C for 45 sec. Sanger’s sequencing was performed from a commercial source (Base Asia Singapore) to detect any change in the DNA sequence of the patient and control group.
2.3. Genotypic and Statistical Analysis
The sequences were analyzed by performing BLAST and were also visualized through BioEdit software (7.2.0) for detection of any mutation. Association analysis between controls and drug-resistant epilepsy patients for rs2298771 was conducted by applying the chi-square test through SHEsis online (http://analysis.bio-x.cn/myAnalysis.php). The distribution of genotypic and allelic frequencies was checked by Hardy–Weinberg Equilibrium (HWE). Association analysis between drug-responder epileptic patients and drug-resistant epileptic patients was also performed to check the distribution of wild and mutant allele.
2.4. In Silico Analysis
Search Tool for the Retrieval of Interacting Genes (STRING) was used to study protein-protein interaction for the selected genes. PolyPhen-2 software was used to check the pathogenicity of the studied polymorphism. PolyPhen-2 predicts the pathogenicity of a polymorphism based on pathogenicity score. If the pathogenicity score is near 1, the polymorphism is predicted to be damaging, and if the pathogenicity score is near to 0, the polymorphism is predicted to be benign (Figure 1).
3. Results
The clinical data of the subjects are presented in Table 2. In this study, there was no association between SCN2Ars17183814 and DRE. Exon 6, 26, and rs10167228 polymorphisms of SCN1A were also not associated with drug-resistant epilepsy. SCN1Ars2298771 polymorphism was observed in patients of DRE. No significant association was shown between rs2298771 polymorphism and DRE. The distribution of frequencies was also checked by HWE, and it was not violated. The distribution of genotypic and allelic frequencies of cases and controls is given in Table 3. The frequency of GG was higher in controls than in cases, while that of GA was higher in cases. No significant difference was shown in the distribution of genotypic and allelic frequencies between drug-resistant and drug-responsive patients (Table 4). The p value for genotypic distribution is 0.14 which showed there is no significant difference between two groups. The results are presented in Table 4. The chromatogram showing mutation is shown in Figure 2.
GTC : generalized tonic clonic.
Data are presented as number (N) and frequency ∗p < 0.01, ∗∗p < 0.001, and ∗∗∗p < 0.0001. OR: odds ratio; CI : confidence interval.
Data are presented as number (N) and frequency. OR: odds ratio; CI : confidence interval.
3.1. In Silico Analysis
Figure 3 shows the protein network highlighting functional association of SCN1A with other proteins. A node is showing a protein, and the edge is showing predicted functional association. Different line colors show the evidence types for the association. Red lines show the fusion evidence, yellow lines show miming evidence, blue line show database evidence, and black lines show coexpression evidence.
Figure 1 shows the pathogenicity score of rs2298771. According to pathogenicity score by PolyPhen-2 which is 0, r2298771 is likely to be benign.
4. Discussion
About twenty percent of epilepsy patients worldwide show resistance to antiepileptic drugs which leads to various psychological problems, i.e., depression and anxiety. Continuous exposure to ineffective drugs can affect the quality of life of drug-resistant patients [Reference Wirrell28]. Most of the antiepileptic drugs work by affecting sodium channels and blocking them, so the genes encoding these channels are the most obvious candidates for studying genetic polymorphisms affecting drug response [Reference Depondt29–Reference Szoeke, Newton and Wood31].
Sodium channel genes are the most important targets for antiepileptic drugs affecting their potency. In this study, the association of sodium channel genes SCN1A and SCN2A polymorphisms with drug-resistant epilepsy was studied. The results showed no significant association of SCN1Ars2298771 with DRE, and also, no association was observed between SCN2A and drug resistance.
A significant association of SCN1Ars2298771 with DRE was observed in the north Indian population where the frequency of the AG genotype was significantly higher in patients suffering from drug-resistant epilepsy than in healthy control individuals, but a weak association was observed between SCN2Ars17183814 and drug-resistant epilepsy [Reference Lakhan, Kumari, Misra, Kalita, Pradhan, Mittal and Mittal14].
In many studies, intronic polymorphisms of SCN1A have also been described to be significantly associated with response to AEDs. Contradictory to our results, a study conducted on 120 epilepsy pediatric patients showed a significant association between intronic SNPs of the SCN1A gene and drug resistance, but no association was observed between exonic SNPs of SCN1A and drug-resistant epilepsy. Three intronic SNPs of SCN1A, rs6730344, rs6732655, and rs10167228, were shown to be potential risk factors for developing drug resistance [Reference Margari, Legrottaglie and Vincenti3].
In accordance with results of our study, a study conducted on epileptic patients to investigate the effect of valproic acid (VPA), one of the widely used drug for treatment of epilepsy, revealed that SCN1Ars2298771 and SCN2Ars17183814 do not significantly affect response to VPA. Hence, these polymorphisms were not detected to be significantly linked with DRE in the Chinese population [Reference Shi, Zhu and Li32].
The lack of association of SCN2Ars17183814 with response to antiepileptic drugs in Chinese epileptic patients was also confirmed in [Reference Li, Zhang and Wu33]. A study conducted on the Scottish population concluded a weak association of SCN2Ars17183814 with DRE [Reference Sills, Mohanraj and Butler34].
This study showed no mutations in exon 6 of the SCN1A gene. The previous studies have reported mutations in exon 6 of SCN1A which lead to truncation of channel protein [Reference Morimoto, Mazaki and Nishimura24]. No significant mutation in exon 26 of SCN1A was observed in this study which is consistent with the fact that there are no mutations reported in exon 26 of the SCN1A gene causing drug-resistant epilepsy. However, mutational analysis showed that frameshift mutation in exon 26 of SCN2A can lead to premature stop codon which can affect the protein structure leading to myoclonic epilepsy [Reference Claes, Del-Favero, Ceulemans, Lagae, Van Broeckhoven and De Jonghe25].
5. Conclusions
Voltage-gated ion channels have an important role in triggering epilepsy, and also, they act as targets for a number of antiepileptic drugs. The current study has shown that rs2298771 of the SCN1A gene may be involved in the epilepsy susceptibility but not in the occurrence of drug-resistant epilepsy in the Pakistani population. However, more studies from different ethnicities in Pakistan and worldwide are required to evaluate the role of sodium-channel genes in developing drug-resistant epilepsy so that potential targets can be explored for the antiepileptic drugs.
Data Availability
The data used to support the findings of this study are available upon request from the corresponding author.
Conflicts of Interest
The authors declare no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors are thankful to all the participants of this study. This work was supported by the Higher Education Commission, Pakistan (No. 8369/Punjab/NRPU/R&D/HEC/2017) under the National Research Program for Universities.