Skip to main content

ORIGINAL RESEARCH article

Front. Neurol., 09 October 2014
Sec. Epilepsy

MDR-1 and MRP2 Gene Polymorphisms in Mexican Epileptic Pediatric Patients with Complex Partial Seizures

\r\n      David Escalante-Santiago,David Escalante-Santiago1,2 \r\n      Iris Anglica Feria-RomeroIris Angélica Feria-Romero2 \r\n      Rosa María Ribas-AparicioRosa María Ribas-Aparicio1 \r\n      Dario Rayo-MaresDario Rayo-Mares3 \r\n      Pietro FagiolinoPietro Fagiolino4 \r\n      Marta VzquezMarta Vázquez4 \r\n      Consuelo Escamilla-NúezConsuelo Escamilla-Núñez5 \r\n      Israel Grijalva-OteroIsrael Grijalva-Otero2 \r\n      Miguel Angel Lpez-GarcíaMiguel Angel López-García2Sandra Orozco-Surez*Sandra Orozco-Suárez2*
  • 1Programa de Biomedicina y Biotecnología Molecular, Escuela Nacional de Ciencias Biológicas del Instituto Politécnico Nacional, Mexico City, Mexico
  • 2Unidad de Investigación Médica en Enfermedades Neurológicas, Hospital de Especialidades, Centro Médico Nacional Siglo XXI, Instituto Mexicano del Seguro Social, Mexico City, Mexico
  • 3Neurología, Hospital de Pediatría, Centro Médico Nacional Siglo XXI, Instituto Mexicano del Seguro Social, Mexico City, Mexico
  • 4Departamento de Ciencias Farmacéuticas de la Facultad de Química, Universidad de la República, Montevideo, Uruguay
  • 5Instituto Nacional de Salude Publica, Cuernavaca, Mexico

Although the Pgp efflux transport protein is overexpressed in resected tissue of patients with epilepsy, the presence of polymorphisms in MDR1/ABCB1 and MRP2/ABCC2 in patients with antiepileptic-drugs resistant epilepsy (ADR) is controversial. The aim of this study was to perform an exploratory study to identify nucleotide changes and search new and reported mutations in patients with ADR and patients with good response (CTR) to antiepileptic drugs (AEDs) in a rigorously selected population. We analyzed 22 samples In Material and Methods, from drug-resistant patients with epilepsy and 7 samples from patients with good response to AEDs. Genomic DNA was obtained from leukocytes. Eleven exons in both genes were genotyped. The concentration of drugs in saliva and plasma was determined. The concentration of valproic acid in saliva was lower in ADR than in CRT. In ABCB1, five reported SNPs and five unreported nucleotide changes were identified; rs2229109 (GA) and rs2032582 (AT and AG) were found only in the ADR. Of six SNPs associated with the ABCC2 that were found in the study population, rs3740066 (TT) and 66744T > A (TG) were found only in the ADR. The strongest risk factor in the ABCB1 gene was identified as the TA genotype of rs2032582, whereas for the ABCC2 gene the strongest risk factor was the T allele of rs3740066. The screening of SNPs in ACBC1 and ABCC2 indicates that the Mexican patients with epilepsy in this study display frequently reported ABCC1 polymorphisms; however, in the study subjects with a higher risk factor for drug resistance, new nucleotide changes were found in the ABCC2 gene. Thus, the population of Mexican patients with AED-resistant epilepsy (ADR) used in this study exhibits genetic variability with respect to those reported in other study populations; however, it is necessary to explore this polymorphism in a larger population of patients with ADR.

Introduction

Epilepsy is a chronic non-communicable disorder of the brain that affects approximately 70 million people worldwide (1, 2). The administration of anti-epileptic drugs (AEDs) is the treatment of choice. In 30% of patients with epilepsy, seizures persist despite polytherapy with more than one AED, adequate monitoring of serum drug levels and the watchful care of a neurologist (3). Hence, the timely identification of patients who do not respond favorably to AEDs is crucial because these patients are affected both by the treatment and by continuing seizures. Although the pathophysiological basis of drug resistance in epilepsy is not clear, some mechanisms for this resistance have been proposed. Following seizures, the overexpression of drug-transporter proteins in the blood brain barrier (BBB) could be mediated by the nuclear receptor PXR (4); alternatively, polymorphisms in these proteins could abrogate their function. The main function of carrier proteins is to allow the passage of compounds through biological barriers such as the BBB. Using ATP, the ABC transporter family of proteins actively transports a wide variety of compounds, including toxins and xenobiotics, across cell membranes. ABC transporter proteins have two functional domains. The first domain is anchored to the cell membrane by its alpha-helical structure; the second, which is a nucleotide-binding domain (NBD), binds ATP, which provides the energy necessary for the conformational change that result in transport of the compound across the membrane. All ABC transporters have at least two transmembrane domains and two NBDs that are highly conserved between species (5).

MDR-1 (P-glycoprotein or P-gp) is an ATP-dependent efflux pump protein that controls the flow of toxins and drugs such as felbamate (FBM), gabapentin (GBP), lamotrigine (LTG), phenytoin (PTH), and topiramate (TPM) into the brain (69). The gene encoding MDR-1 is ABCB1, which is located on chromosome 7q21.12 and consists of 209,461 base pairs (bp) and 29 exons (NCBI reference sequence NG_011513.1). The most extensively studied single-nucleotide polymorphisms (SNPs) in the ABCB1 gene are rs1128503 (C1236T), rs2032582 (G2677T/A), and rs1045642 (C3435T); the latter has received increased attention as a critical variable in resistance to AEDs (1016) MRP2 is another efflux transporter protein present in the BBB. Its gene, ABCC2, is located on chromosome 10q24 and contains 69,000 bp representing 33 exons (NG_011798.1). The MRP2 transporter protein recognizes CBZ, LTG, and FBM; its most relevant polymorphisms are rs2273697 (G1249A) and rs3740066 (C3972T). The latter polymorphism reduces the transport of carbamazepine (1720). The aim of this study was to perform a non-inferential exploratory study to identify reported nucleotide changes and explore new ones in a rigorously selected Mexican population of patients with AED-resistant epilepsy (ADR) and patients with good response to AEDs (CTR).

Materials and Methods

Patients and Sample Collection

This project was approved by the National Research Committee of the Instituto Mexicano del Seguro Social (Mexican Social Security Institute, IMSS). An observational study with 22 patients with ADR (cases) and 7 patients with good response to AEDs (controls) was performed, while avoiding inbreeding between the patient’s biological parents.

The inclusion criteria for ADR were patients with demonstrable epileptic focus through EEG and radiological focal structural lesion to discard and classified as resistant to pharmacological treatment, treated with two or more drugs (Table 1) at appropriate doses, serum levels within therapeutic range for at least 6 months of continuous treatment, and under the supervision of a neurologist pediatrician, with a frequency of three seizures per month, who attend Neurology Service, Hospital of Pediatrics, National Medical Center, Century XXI, who is of Mexican origin like their parents and grandparents, from 1 to 16 years of age, both genders, whereas for asymptomatic control patients was seizure-free for at least 6 months before the study.

TABLE 1
www.frontiersin.org

Table 1. Demographics and clinical data of patients with AEDs-resistant epilepsy and patients with good response to AEDs.

Both groups had complex partial seizures, showed similar clinical variables, were previously selected under medical supervision, and gave written informed consent. Exclusion criteria were all types of seizures that are not complex partial, little adherence to treatment.

A total of 15 ml of blood was collected from each patient via antecubital vein using an EDTA tube, and two fractions of 0.5–2.0 ml of saliva were collected upon stimulation with citric acid. Samples that were used to monitor drug concentrations were taken before the first dose in the morning.

Amplification and Sequencing

Genomic DNA was extracted from leukocytes in patient blood samples using a commercial kit (Genomic DNA Purification kit, Thermo Scientific®, USA), according to the supplier’s recommendations. The polymorphisms (SNPs) from exons 1 (rs9282564), 10 (rs2229109), 11 (rs1128503), 13 (rs28381902), 20 (rs2032582), and 25 (rs1045642) were located in GenBank accession CCDS5608.1, corresponding to the ABCB1 gene, and exons 10 (rs2273697), 25 (rs8187692, rs17222723), 28 (rs3740066), 29 (rs7899457), and 32 (rs8187710) were located in GenBank accession CCDS7484.1, corresponding to the ABCC2 gene. These polymorphisms were selected based on their frequency in the Hispanic or Mexican population, as noted in NCBI GenBank. Each amplification reaction was performed with 5 μl (5 ng/μl) of genomic DNA in a 50 μl total reaction volume containing 5 μl of 10× reaction buffer, MgCl2 (concentration depending on the primer set), 1 μl 10 mMdNTPs (Thermo Scientific®, USA), 1 μl of each 10 mM flanking primer (Invitrogen®, USA), and 1 μl of Taq polymerase (Thermo Scientific®, USA). PCR was performed in a DNA engine system®thermocycler (Bio-Rad®, USA) with a cycle program at 94°C for 3 min, 38 cycles of 94°C for 30 s, Th for 35 s, 72°C for 35 s, and one extension cycle of 10 min at 72°C. The primers, MgCI2 concentration, and Th for each amplification reaction are listed in Table 2.

TABLE 2
www.frontiersin.org

Table 2. Primers used on ABCB1 and ABCC2 fragments amplification.

The amplification products were purified with the GeneJET Gel Extraction kit (Thermo Scientific®, USA) and sequenced in a 5 μl total reaction volume using 100 ng of amplicon, Fw or Rv primer, and Big Dye 3.1 sequencing Ready Reaction Kit (Applied Biosystems, USA) according to the manufacturer’s recommendations. Sequencing was performed in a DNA engine system®thermocycler (Bio-Rad®, USA) with a cycle program at 94°C for 1 min, 30 cycles of 94°C for 10 s, 50°C for 5 s, 60°C for 4 min, and finally the fragments were analyzed in Automated Sequencing ABI Prism 310 Gene Analyzer (Applied Biosystems, USA), previous purification with DyeEx 2.0 Spin Kit (Qiagen, Germany).

Drug Monitoring

Plasma samples were used for VPA and CBZ detection by immunofluorescence with a polarized light system (FPIA, Dade Behring, USA), according to the supplier’s recommendations.

Saliva samples were used to detection concentration of CBZ and PHT (determined by Axsym®, Abbot Laboratories®, USA) according to the supplier’s recommendations, VPA, and LEV [determined by high performance liquid chromatography (HPLC)]. For VPA, 30 μl of an internal standard solution (98% octanoic acid diluted 1/500 with methanol), 200 μl of phosphoric acid, and 5 ml of cyclohexane were added to 1 ml of saliva and vortexed for 2 min. After centrifugation, a back-extraction was performed by adding 100 μl of 1% TEA to the organic phase supernatant. The mixture was vortexed and centrifuged. Finally, 20 μl of the aqueous phase was directly injected into a Dionex Ultimate 3000 Series. A Phenomenex Luna 5 μm CN column (15 cm × 4.6 mm) was used as the stationary phase. The mobile phase was a mixture of 40 mM phosphate potassium, monobasic, pH 3.4: Acetonitrile (90:10). The column compartment was kept at 40°C, and the wavelength of detection was 210 nm. For LEV, 50 μl of internal standard solution (diethylbarbituric acid, 12 mg/l) was added to 500 μl of saliva. For the extraction, 1.5 ml of ethyl acetate was added, and the mixture was vortexed for 1 min. After centrifugation, the supernatant was separated and dried under nitrogen stream at 37–40°C. Dry residue was dissolved with 100 μl of mobile phase, and 20 μl was injected into a Shimadzu LC-6A.A Phenomenex® Luna C18 column (5 μm, 100 Å, 150 mm × 4.6 mm), which was used as a reversed stationary phase. The mobile phase was a mixture of 50 mM phosphate buffer and acetonitrile (72:8).

Statistical Analysis

Clinical data from patients were analyzed using non-parametric and unpaired tests (p < 0.05); allelic and genotype relative frequencies were analyzed using the Fisher exact test in Graph Pad Prism software version 4.00 (GraphpadPrisma Software Inc.® San Diego, CA, USA). Hardy–Weinberg equilibrium (HWE) values were derived from chi-square tests, and the association between odds ratios (ORs) and confidence intervals was determined using software available at http://ihg.gsf.de/cgi-bin/hw/hwa1.pl.

Results

Clinical Data

Because previous studies suggest that external factors such as gender, age of onset, number of seizures, and other factors influence the response to AED treatment, the selection of a homogeneous population for the study was an important aspect of this work (21, 22). There was no significant difference in male: female frequency between the two groups used in the study. Neither the average age of the patients (10.57 ± 1.56 for patients with good response to AEDs and 10.23 ± 0.78 for patients with ADR; p = 0.42) or the average age of onset of seizures (4.57 ± 1.43 years in patients with good response to AEDs and 4.38 ± 0.81 patients with ADR; p = 0.34) differed significantly. Furthermore, the average seizure-free time in patients with good response to AEDs was 1.9 years, while the average time with seizures in patients with ADR was 5.8 years with a frequency of 63.4 seizures per month. The AEDs VPA, CBZ, LTG, PHT, and LEV were used with relative frequencies of 0.71, 0.14, 0.29, 0.0, and 0.14, respectively, in patients with good response to AEDs and with relative frequencies of 0.59, 0.09, 0.0, 0.18, and 0.23, respectively, in patients with ADR. The most commonly used AED was VPA, but VPA and CBZ were prescribed to both groups (Table 1). The daily dose of VPA (814.3 ± 212.1 mg/day in patients with good response to AEDs and 598.0 ± 129.7 mg/day in patients with ADR; p = 0.206) was similar in the two groups. The two groups showed no significant difference in mean plasma VPA concentrations (35.17 ± 18.04 mg/l patients with good response to AEDs and 31.62 ± 8.03 mg/l in patients with ADR; p = 0.42); however, the average saliva drug concentration was lower in patients with ADR (1.68 ± 0.59 mg/l patients with good response to AEDs and 0.57 ± 0.12 mg/l in patients with ADR; p = 0.01).

Polymorphisms in the ABCB1 and ABCC2 Genes in Hispanic and Mexican Populations

In this study, chromatograms from amplicons were analyzed to identify expected polymorphisms and to determine the relative frequencies of alleles and genotypes, as described in Table 3. For ABCB1, only the polymorphisms rs2229109, rs1128503, rs28381902, rs2032582, and rs1045642, but not rs9282564, were detected, whereas rs3740066 but not rs2273697, rs8187692, rs17222723, rs7899457, or rs8187710 were found in ABCC2 in our patients, unlike the Hispanic/Mexican populations initially selected by the NCBI GenBank. In addition to these known polymorphisms, we identified rs2214102, rs2235047, and rs2235048 of ABCB1 by sequencing the entire fragment including the exon and parts of the introns.

TABLE 3
www.frontiersin.org

Table 3. Genotype and allelic relative frequencies of SNPs, MAF value, amino acid changes, and OR found in ABCB1 and ABCC2 genes.

Another comparison was made between the relative frequencies of the minor allele (MAF) in our population and in the population represented by the NCBI GenBank data. The most significant similarities were found in rs1128503 and rs3740066 in ADR and in rs1045642 and rs2235048 in CTR (Table 3).

An important aspect of genetic alteration is whether a nucleotide change results in a change in an amino acid residue and whether this, in turn, structurally modifies the protein. In this respect, the SNPs rs2229109 and rs2032582 are of interest due to the possible modifications that they may produce in protein MDR1 (Table 3).

The relative frequencies of the studied alleles were also compared with their frequencies in the Hispanic and Mexican reference populations reported in the NCBI GenBank (Table 4). In our population, SNP rs2229109 in exon 10 of ABCB1 occurred with a relative frequency of 0.09 and was only present in heterozygous form in the ADR group (Table 3). Although a reference for the Mexican population was not available, we observed that our results showed a twofold higher minor allele frequency with respect to sample ss48292345 from the Hispanic population (0.043) in a similar number of individuals (Table 4). The SNPs rs1128503 and rs2032582 in exons 11 and 20, respectively, of the ABCB1 gene showed similar relative frequencies in the CTR (0.71) and ADR (0.82) patient groups; however, there was a difference in the frequencies of homozygotes and heterozygotes in the CTR group (Table 3). Furthermore, there were no Mexican individuals with AG heterozygous mutations in SNP rs2032582, and the AT heterozygous mutation was absent from both populations (Table 4).

TABLE 4
www.frontiersin.org

Table 4. Genotype and allelic population frequencies of SNPs found in ABCB1 and ABCC2 genes.

In studies of patients with ADR, SNP rs1045642 in exon 25 has been shown to be relevant in association studies with drug resistant; in our population, the relative frequencies of this SNP were similar in the CTR and ARE groups at both the allelic and the genotypic levels. Furthermore, the SNPs rs2235047 and rs2235048 were identified in an intron between exons 25 and 26. The heterozygous form of rs2235047 showed a higher frequency ADR than in CTR (0.400 vs. 0.286), as seen in Table 3; however, the most outstanding difference in this SNP was in comparison to reference populations, relative to which it was more than three times as frequent (Table 4); this is more significant in our population than rs1045642. The relative frequency of SNP rs2235048 was the same as that of rs1045642; when SNP rs1045642 showed the TT, TC, or CC genotype, we found the genotype CC, CT, or TT, respectively, in the position of SNP rs2235048. This association may be relevant to our population, which so far has not shown a clear similarity to reference populations.

SNP rs3740066 in exon 28 of the ABCC2 gene occurred with a relative frequency of 0.13 as a heterozygous polymorphism in the controlled group and with a relative frequency of 0.41 as homozygous and heterozygous polymorphisms in the drug-resistant group; the observed heterozygous ratio (ADR: CTR) was 2:1. In addition, the unexpected SNP rs2214102 located in the 5′UTR region of the ABCB1 gene was present in the homozygous form in both groups; only 5% of the drug-resistant patients had the heterozygous form. However, this finding is not relevant because it is similar to a previous report in a Mexican population.

Some other differences observed in this study are worth noting. At SNP rs2032582, the allelic and genotypic distribution in the controlled group was more similar to reported populations in NCBI GenBank than was that of the ADR; in particular, the heterozygous A > T mutation was not found in other populations similar to ours. Furthermore, although allele A of SNP rs2229109 was present at a lower frequency, our population showed an allele frequency more than 2.5-fold greater than the reported MAF (0.045 vs. 0.017).

New Polymorphisms of the ABCC2 Gene in Mexican Epileptic Pediatric Patients

By sequencing amplicons from exons 28 and 29 of the ABCB1 and ABCC2 genes, we found other single-nucleotide changes (SNCs) that have not been reported in the NCBI GenBank database (Table 2, RefSeqGene NG_011798.1). The SNC located in position 66744T > G adjacent to rs3740066 in exon 28 is present as a homozygous wild-type or heterozygous mutation, allowing the formation of four different triplets: ATC, ATT, AGT, and AGC. The ATC and ATT triplets, which show relative frequencies of 0.929 and 0.682 for triplet ATC and 0.071 and 0.205 for triplet ATT in CTR and ADR, respectively, are not associated with an amino acid change (Ile1324Ile). The AGT and AGC triplets produce an amino acid change (Ile1324Ser, Figure 1A); these occur at relative frequencies of 0.068 and 0.045, respectively, but only in drug-resistant patients.

FIGURE 1
www.frontiersin.org

Figure 1. Chromatograms of single-nucleotide changes found in the ABCC2 gene from Mexican pediatric patients with epilepsy. (A) The 66744 T > G change in located in exon 28. (B) The 67967 C > A change in located in exon 29. (C) The 68049 T > A change in located in exon 29. (D) The 68072 C > A change in located in exon 29. (E) The 68088 G > C change in located in intron 29. The position indicates the nucleotide change in the ABCC2 gene (RefSeqGene NG_011798.1).

For the SNC 67967C > A in exon 29, no amino acid change (Ser1342Ser) was observed due to nucleotide substitution (Figure 1B); however, SNC 68072C > A resulted in the amino acid change Pro1381Thr (Figure 1D). This allele was present as a heterozygous mutation in both groups (1.0 in the CTR group, 0.86 in the ADR group). The SNC 68049T > A in exon 29 resulted in a Leu1373His change (Figure 1C); this allele was heterozygous (0.43 in the CTR group, 0.32 in the ADR group), but the wild-type allele predominated in both groups (>50%). These nucleotide changes are located in the region of the gene coding for the fourth domain of the protein, corresponding to its cytoplasmic region. In contrast, the SNC 68088G > C (Figure 1E), which is located between exons 29 and 30, occurred in our study population as a homozygous mutation with a relative frequency of 0.29 in the CTR group and a relative frequency of 0.05 in ADR group.

Discussion

The pediatric epilepsy patients admitted to the Neurology Service, Hospital of Pediatrics, National Medical Center, Century XXI, were referred from family clinics and general hospitals for inadequate control of epileptic seizures. In most cases, the patients received dosage adjustments or changes in drug regimen, drug monitoring, and electroencephalographic studies and were followed closely for adequate control of their seizures. However, some patients continued to interictal epileptiform discharges on EGG and/or the occurrence of seizures, and some gradually deteriorated. The results obtained with a small group of patients who did not respond to at least three different drug regimens but were followed for 4.7 ± 3.8 (CI = 1–14) years are presented in this work.

Fagiolino et al. (23) suggested that the quantitation of drugs in saliva can be used to identify drug-resistant patients. The hypothesis presented by these authors is based on the drug concentration circulating in free plasma, which irrigates the glandular capillary tissue, resulting in equilibrium with the aqueous spaces. Considering our results with respect to this hypothesis, we found a lower concentration of free VPA in drug-resistant patients, suggesting a possible deficiency in drug transport in these patients.

The overexpression of P-gp and MRP2 in resected tissue from patients who have undergone surgery for epilepsy suggests that this phenomenon may be conditional upon the presence of polymorphisms in the genes ABCB1 and ABCC2 and not simply to exposure to xenobiotics that recognize these transporter proteins. However, the association between polymorphisms of the ABCB1 and ABCC2 genes and anti-epileptic drug-resistant epilepsy is controversial; SNPs rs1128503 (C1236T), rs2032582 (G2677T/A), and especially rs1045642 (C3435T) in the ABCB1 gene have not shown a clear relationship to anti-epileptic drug-resistant epilepsy from different populations or ethnic groups (13, 2427). This may be because, as has been reported, SNP rs1045642 is more relevant as a haplotype with rs1128503 and rs2032582. Zimprichet et al. (28), studying European patients with temporal lobe epilepsy, demonstrated a relationship between the CGC haplotype at SNPs rs1128503, rs2032582, and rs1045642, respectively, and resistance to several AEDs. In another report, Seo et al. (10) demonstrated that Japanese with anti-epileptic drug-resistant epilepsy showed a tendency to possess the haplotype TTT at SNPs rs1128503, rs2032582, and rs1045642 when CBZ was employed. To enable comparison with those reports, haplotype frequency was calculated in our population. Both haplotypes (CGC and TTT) were found in our population. Within each group, the two haplotypes occurred with similar frequencies (CGC = 0.619, TTT = 0.523 in ADR and CGC = 0.714, TTT = 0.714 in CTR), but both were present at higher frequencies in the CTR group.

For patients who exhibit a mutated allele, OR values can be used to represent relative risk factors among the ADR and CTR patients. The polymorphisms or nucleotide changes, in descending order of risk factor magnitude, are rs2214102 > rs2032582 [T ↔ A] > rs2229109 > rs1128503 > rs2032582 [T ↔ G] > rs2235047 in ABCB1, and rs3740066 > 66744 T > G in ABCC2; relative to heterozygote genotype rs2032582 [T ↔ A] > rs2214102 > rs2229109 > rs2032582 [T ↔ G] > rs11 28503 > rs2235047 in ABCB1, and 66744 T > G > rs3740066 in ABCC2; finally, with regard to the homozygote genotype, rs1128503 > rs2032582 [T ↔ G] in ABCB1 and rs3740066 in ABCC2 (Table 3).

Although the number of patients is lower in our work, the results are similar to studies with larger numbers of patients (25, 29, 30), where the lack of association of ABCB1 polymorphisms with anti-epileptic drug resistance (indicated by us as a risk factor) it has been pointed out, including two Indian studies (31, 32) and a meta-analysis (26). In contrast, Taur et al. (33) and Siddiqui et al. (34) reported association with anti-epileptic drug resistance. On the other hand, we find that the ABCC2 is a greater risk factor for increased anti-epileptic drug resistance than ABCB1, although in previous studies, we only share the SNP rs3740066, not representing a major risk factor as the high activity promoter variant c-24C > T; these data suggest that ABCC2 polymorphism may affect the response to AEDs.

This study employed a different strategy, rigorously selecting pediatric patients with ADR and patients with good seizure control to decrease the potential influence of confounding variables associated with the disease, to develop a study utilizing a non-inferential exploratory search based on the polymorphisms already reported and on unreported new characteristics of the Mexican population. Because our population has been little studied, this will allow us to conduct association studies with new SNPs and AEDs used to treat different forms of epilepsy. To avoid sequencing entire genes, we selected exons containing polymorphisms that have been previously reported in Mexican and Hispanic populations, including the P-gp and MRP2 proteins, which are associated with drug resistance. The most important outcome of this study was the identification of five new changes in the ABCC2; these are 66744T > G, 67967C > A, 68072C > A, 68049C > A, and G68088C. The nucleotide change 66744T > G was found only in ADR, and the risk factor associated with it was one of the highest found in the study. Moreover, the genetic variability in SNP rs2032582 in ADR was different and unexpected based on previous reports in the Mexican population in that the A allele was found. This result requires replication in a larger pool of Mexican patients.

Conflict of Interest Statement

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.

References

1. Ngugi AK, Bottomley C, Kleinschmidt I, Sander JW, Newton CR. Estimation of the burden of active and life-time epilepsy: a meta-analytic approach. Epilepsia (2010) 51:883–90. doi: 10.1111/j.1528-1167.2009.02481.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

2. World Health Organization (WHO). Epilepsy: Fact Sheet No 999. (2012). Available from: http://www.who.int/mediacentre/factsheets/fs999/en/index.html

3. Kwan P, Baum L, Wong V, Ng PW, Lui CH, Sin NC, et al. Association between ABCB1 C3435T polymorphism and drug-resistant epilepsy in Han Chinese. Epilepsy Behav (2007) 11:112–7. doi:10.1016/j.yebeh.2007.04.013

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

4. Stepien KM, Tomaszewski M, Tomaszewska J, Czuczwar SJ. The multidrug transporter P-glycoprotein in pharmacoresistance to antiepileptic drugs. Pharmacol Rep (2012) 64:1011–9. doi:10.1016/S1734-1140(12)70900-3

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

5. Aller SG, Yu J, Ward A, Weng Y, Chittaboina S, Zhuo R, et al. Structure of P-glycoprotein reveals a molecular basis for poly-specific drug binding. Science (2009) 323:1718–22. doi:10.1126/science.1168750

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

6. Rivers F, O’Brien TJ, Callaghan R. Exploring the possible interaction between anti-epilepsy drugs and multidrug efflux pumps; in vitro observations. Eur J Pharmacol (2008) 598:1–8. doi:10.1016/j.ejphar.2008.09.014

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

7. Löscher W, Klotz U, Zimprich F, Schmidt D. The clinical impact of pharmacogenetics on the treatment of epilepsy. Epilepsia (2009) 50:1–23. doi:10.1111/j.1528-1167.2008.01716.x

CrossRef Full Text

8. Moerman L, Wyffels L, Slaets D, Raedt R, Boon P, De Vos F. Antiepileptic drugs modulate P-glycoproteins in the brain: a mice study with 11C-desmethylloperamide. Epilepsy Res (2011) 94:18–25. doi:10.1016/j.eplepsyres.2010.12.013

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

9. Aronica E, Sisodiya SM, Gorter JA. Cerebral expression of drug transporters in epilepsy. Adv Drug Deliv Rev (2012) 64:919–29. doi:10.1016/j.addr.2011.11.008

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

10. Seo T, Ishitsu T, Ueda N, Nakada N, Yurube K, Ueda K, et al. ABCB1 polymorphisms influence the response to antiepileptic drugs in Japanese epilepsy patients. Pharmacogenomics (2006) 7:551–61. doi:10.2217/14622416.7.4.551

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

11. Kwan P, Arzimanoglou A, Berg AT, Brodie MJ, Allen Hauser W, Mathern G, et al. Definition of drug resistant epilepsy: consensus proposal by the ad hoc task force of the ILAE commission on therapeutic strategies. Epilepsia (2010) 51:1069–77. doi:10.1111/j.1528-1167.2009.02397.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

12. Nurmohamed L, Garcia-Bournissen F, Buono RJ, Shannon MW, Finkelstein Y. Predisposition to epilepsy does the ABCB1 gene play a role? Epilepsia (2010) 51:1882–5. doi:10.1111/j.1528-1167.2010.02588.x

CrossRef Full Text

13. Haerian BS, Lim KS, Mohamed EH, Tan HJ, Tan CT, Raymond AA, et al. Lack of association of ABCB1 and PXR polymorphisms with response to treatment in epilepsy. Seizure (2011) 20:387–94. doi:10.1016/j.seizure.2011.01.008

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

14. Haerian BS, Lim KS, Tan HJ, Mohamed EH, Tan CT, Raymond AA, et al. Association between ABCB1 polymorphism and response to sodium valproate treatment in Malaysian epilepsy patients. Epileptic Disord (2011) 13:65–75. doi:10.1684/epd.2011.0419

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

15. Kumari R, Lakhan R, Garg RK, Kalita J, Misra UK, Mittal B. Pharmacogenomic association study on the role of drug metabolizing, drug transporters and drug target gene polymorphisms in drug-resistant epilepsy in a north Indian population. Ind J Hum Genet (2011) 1:S32–40. doi:10.4103/0971-6866.80357

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

16. Sayyah M, Kamgarpour F, Maleki M, Karimipoor M, Gharagozli K, Shamshiri AR. Association analysis of intractable epilepsy with C3435T and G2677T/A ABCB1 gene polymorphisms in Iranian patients. Epileptic Disord (2011) 13:155–65. doi:10.1684/epd.2011.0443

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

17. Sterjev Z, Trencevska GK, Cvetkovska E, Petrov I, Kuzmanovski I, Ribarska JT, et al. The association of C3435T single-nucleotide polymorphism, Pgp-glycoprotein gene expression levels and carbamazepine maintenance dose in patients with epilepsy. Neuropsychiatr Dis Treat (2012) 8:191–6. doi:10.2147/NDT.S28285

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

18. Kim WJ, Lee JH, Yi J, Cho YJ, Heo K, Lee SH, et al. A nonsynonymous variation in MRP2/ABCC2 is associated with neurological adverse drug reactions of carbamazepine in patients with epilepsy. Pharmacogenet Genomics (2010) 20:249–56. doi:10.1097/FPC.0b013e328338073a

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

19. Ufer M, von Stülpnagel C, Muhle H, Haenisch S, Remmler C, Majed A, et al. Impact of ABCC2 genotype on antiepileptic drug response in Caucasian patients with childhood epilepsy. Pharmacogenet Genomics (2011) 21:624–30. doi:10.1097/FPC.0b013e3283498131

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

20. Qu J, Zhou BT, Yin JY, Xu XJ, Zhao YC, Lei GH, et al. ABCC2 polymorphisms and haplotype are associated with drug resistance in Chinese epileptic patients. CNS NeurosciTher (2012) 18:647–51. doi:10.1111/j.1755-5949.2012.00336.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

21. Ramos-Lizana J, Aguirre-Rodríguez J, Aguilera-López P, Cassinello-García E. Recurrence risk after a first remote symptomatic unprovoked seizure in childhood: a prospective study. Dev Med Child Neurol (2009) 51:68–73. doi:10.1111/j.1469-8749.2008.03124.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

22. Martínez-Juárez IE, López-Zapata R, Gómez-Arias B, Bravo-Armenta E, Romero-Ocampo L, Estévez-Cruz Z, et al. Epilepsia farmacorresistente: uso de la nueva definición y factores de riesgo relacionados. Estudio en población mexicana de un centro de tercer nivel. Rev Neurol (2012) 54:159–66.

Pubmed Abstract | Pubmed Full Text

23. Fagiolino P, Vázquez M, Maldonado C. Aspectos farmacocinéticos del tratamiento antiepiléptico y su monitoreo mediante el uso de saliva. In: Beas Zárate C, Ureña Guerrero M, et al., editors. En: Tópicos de Actualización en Neurobiología: Excitotoxicidad y Cognición en Enfermedades Neurogenerativas: Aspectos Básicos, Clínicos y Sociales. Guadalajara: Universidad de Guadalajara (2010). p. 381–400.

24. Chen P, Yan Q, Xu H, Lu A, Zhao P. The effects of ABCC2 G1249A polymorphism on the risk of resistance to antiepileptic drugs: a meta-analysis of the literature. Genet Test Mol Biomarkers (2014) 18:106–11. doi:10.1089/gtmb.2013.0362

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

25. Grover S, Bala K, Sharma S, Gourie-Devi M, Baghel R, Kaur H, et al. Absence of a general association between ABCB1 genetic variants and response to antiepileptic drugs in epilepsy patients. Biochimie (2010) 92:1207–12. doi:10.1016/j.biochi.2010.04.008

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

26. Haerian BS, Roslan H, Raymond AA, Tan CT, Lim KS, Zulkifli SZ, et al. ABCB1 C3435T polymorphism and the risk of resistance to antiepileptic drugs in epilepsy: a systematic review and meta-analysis. Seizure (2010) 19:339–46. doi:10.1016/j.seizure.2010.05.004

CrossRef Full Text

27. Dong L, Luo R, Tong Y, Cai X, Mao M, Yu D. Lack of association between ABCB1 gene polymorphisms and pharmacoresistant epilepsy: an analysis in a western Chinese pediatric population. Brain Res (2011) 1391:114–24. doi:10.1016/j.brainres.2011.03.028

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

28. Zimprich F, Sunder-Plassmann R, Stogmann E, Gleiss A, Dal-Bianco A, Zimprich A, et al. Association of an ABCB1 gene haplotype with pharmacoresistance in temporal lobe epilepsy. Neurology (2004) 63:1087–9. doi:10.1212/01.WNL.0000141021.42763.F6

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

29. Tan NC, Heron SE, Scheffer IE, Pelekanos JT, McMahon JM, Vears DF, et al. Failure to confirm association of a polymorphism in ABCB1 with multi-drug resistant epilepsy. Neurology (2004) 63:1090–2. doi:10.1212/01.wnl.0000137051.33486.c7

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

30. Szoeke C, Sills GJ, Kwan P, Petrovski S, Newton M, Hitiris N, et al. Multidrug-resistance genotype (ABCB1) and seizure recurrence in newly treated epilepsy: data from International Pharmacogenetic Cohorts. Epilepsia (2009) 50:1689–96. doi:10.1111/j.1528-1167.2009.02059.x

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

31. Lakhan R, Misra UK, Kalita J, Pradhan S, Gogtay NJ, Singh MK, et al. No association of ABCB1 polymorphisms with drug-refractory epilepsy in a north Indian population. Epilepsy Behav (2009) 14:78–82. doi:10.1016/j.yebeh.2008.08.019

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

32. 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:255–60. doi:10.2133/dmpk.24.255

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

33. Taur SR, Kulkarni NB, Gandhe PP, Thelma BK, Ravat SH, Gogtay N, et al. Association of polymorphisms of CYP2C9, CYP2C19, and ABCB1, and activity of P-glycoprotein with response to anti-epileptic drugs. J Postgrad Med (2014) 60(3):265–9. doi:10.4103/0022-3859.138739

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

34. Siddiqui A, Kerb R, Weale ME, Brinkmann U, Smith A, Goldstein DB. Association of multidrug resistance in epilepsy with a polymorphism in the drug-transporter gene ABCB1. N Engl J Med (2003) 348:1442–8. doi:10.1056/NEJMoa021986

Pubmed Abstract | Pubmed Full Text | CrossRef Full Text

Keywords: drug-resistant, epilepsy, ABCB1, ABCC2, Mexican patients, single-nucleotide polymorphisms, anti-epileptic drugs

Citation: Escalante-Santiago D, Feria-Romero IA, Ribas-Aparicio RM, Rayo-Mares D, Fagiolino P, Vázquez M, Escamilla-Núñez C, Grijalva-Otero I, López-García MA and Orozco-Suárez S (2014) MDR-1 and MRP2 gene polymorphisms in Mexican epileptic pediatric patients with complex partial seizures. Front. Neurol. 5:184. doi: 10.3389/fneur.2014.00184

Received: 24 July 2014; Accepted: 09 September 2014;
Published online: 09 October 2014.

Edited by:

Fernando Cendes, University of Campinas, Brazil

Reviewed by:

Marilisa M Guerreiro, University of Campinas, Brazil
Rodrigo Secolin, University of Campinas, Brazil

Copyright: © 2014 Escalante-Santiago, Feria-Romero, Ribas-Aparicio, Rayo-Mares, Fagiolino, Vázquez, Escamilla-Núñez, Grijalva-Otero, López-García and Orozco-Suárez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Sandra Orozco-Suárez, Unidad de Investigación Médica en Enfermedades Neurológicas, Hospital de Especialidades, Centro Médico Nacional Siglo XXI, IMSS. Av Cuauhtémoc No. 330-4, Col. Doctores, 06720-25, Mexico city, Mexico e-mail: sorozco5@hotmail.com

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.