Next Article in Journal
Similar Features, Different Behaviors: A Comparative In Vitro Study of the Adipogenic Potential of Stem Cells from Human Follicle, Dental Pulp, and Periodontal Ligament
Next Article in Special Issue
Pharmacogenetics to Avoid Adverse Reactions in Cardiology: Ready for Implementation?
Previous Article in Journal
Identification of ROBO2 as a Potential Locus Associated with Inhaled Corticosteroid Response in Childhood Asthma
Previous Article in Special Issue
A Pharmacogenetic Study of CYP2C19 in Acute Coronary Syndrome Patients of Colombian Origin Reveals New Polymorphisms Potentially Related to Clopidogrel Therapy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Association of HLA-B*51:01, HLA-B*55:01, CYP2C9*3, and Phenytoin-Induced Cutaneous Adverse Drug Reactions in the South Indian Tamil Population

1
Department of Clinical Pharmacy, Faculty of Pharmaceutical Sciences, Prince of Songkla University, Songkhla, Hat Yai 90110, Thailand
2
Department of Immunology, School of Biological Sciences, Madurai Kamaraj University, Madurai 625021, Tamil Nadu, India
3
Division of Pharmacogenomics and Personalized Medicine, Department of Pathology, Faculty of Medicine, Ramathibodi Hospital, Mahidol University, Bangkok 10400, Thailand
4
Laboratory for Pharmacogenomics, SomdechPhraDebaratana Medical Center (SDMC), Ramathibodi Hospital, Bangkok 10400, Thailand
5
Department of Neurology, Meenakshi Mission Hospital and Research Center, Madurai 625107, Tamil Nadu, India
6
Department of Pharmaceutical Chemistry, Faculty of Pharmaceutical Sciences, Prince of Songkla University, Songkhla, Hat Yai 90110, Thailand
*
Author to whom correspondence should be addressed.
J. Pers. Med. 2021, 11(8), 737; https://doi.org/10.3390/jpm11080737
Submission received: 1 July 2021 / Revised: 25 July 2021 / Accepted: 26 July 2021 / Published: 28 July 2021
(This article belongs to the Special Issue Pharmacogenetics to Avoid Adverse Drug Reactions)

Abstract

:
Phenytoin (PHT) is one of the most commonly reported aromatic anti-epileptic drugs (AEDs) to cause cutaneous adverse reactions (CADRs), particularly severe cutaneous adverse reactions (SCARs). Although human leukocyte antigen (HLA)-B*15:02 is associated with PHT-induced Steven Johnson syndrome/toxic epidermal necrosis (SJS/TEN) in East Asians, the association is much weaker than it is reported for carbamazepine (CBZ). In this study, we investigated the association of pharmacogenetic variants of the HLA B gene and CYP2C9*3 with PHT-CADRs in South Indian epileptic patients. This prospective case-controlled study included 25 PHT-induced CADRs, 30 phenytoin-tolerant patients, and 463 (HLA-B) and 82 (CYP2C9*3) normal-controls from previous studies included for the case and normal-control comparison. Six SCARs cases and 19 mild-moderate reactions were observed among the 25 cases. Pooled data analysis was performed for the HLA B*51:01 and PHT-CADRs associations. The Fisher exact test and multivariate binary logistic regression analysis were used to identify the susceptible alleles associated with PHT-CADRs. Multivariate analysis showed that CYP2C9*3 was significantly associated with overall PHT-CADRs (OR = 12.00, 95% CI 2.759–84.87, p = 003). In subgroup analysis, CYP2C9*3 and HLA B*55:01 were found to be associated with PHT-SCARs (OR = 12.45, 95% CI 1.138–136.2, p = 0.003) and PHT-maculopapular exanthema (MPE) (OR = 4.041, 95% CI 1.125–15.67, p = 0.035), respectively. Pooled data analysis has confirmed the association between HLA B*51:01/PHT-SCARs (OR = 6.273, 95% CI 2.24–16.69, p = <0.001) and HLA B*51:01/PHT-overall CADRs (OR = 2.323, 95% CI 1.22–5.899, p = 0.037). In this study, neither the case nor the control groups had any patients with HLA B*15:02. The risk variables for PHT-SCARs, PHT-overall CADRs, and PHT-MPE were found to be HLA B*51:01, CYP2C9*3, and HLA B*55:01, respectively. These alleles were identified as the risk factors for the first time in the South Indian Tamil population for PHT-CADRs. Further investigation is warranted to establish the clinical relevance of these alleles in this population with larger sample size.

1. Introduction

Phenytoin (PHT) is still the most effective treatment for generalized tonic-clonic seizures (GTCS) despite newer anti-epileptic drugs’ (AED) availability [1,2]. However, cutaneous adverse drug reactions (CADRs) may limit its use; the estimated relative risk of PHT-severe cutaneous adverse reactions (SCARs) was reported to be 13% [3]. The prevalence of CADRs ranges from 2–5% in India [4]. Further, PHT and carbamazepine (CBZ) are the worst offenders of CADRs, with incidence rates of 13 and 18%, respectively [5]. The contributing factors to CADRs can be both genetic and non-genetic. In 2004, the association between HLA B*15:02 and CBZ-Steven Johnson syndrome/toxic epidermal necrosis (SJS/TEN) was reported among the Han Chinese populations [6]. Later this association has been confirmed with PHT-SJS/TEN in the Thai and Chinese Asian population [7,8]. As a result, the FDA issued a warning for HLA B*15:02/PHT-SJS/TEN cross-reactivity. However, the strength of this association is weaker than CBZ-SJS/TEN and not demonstrated well enough in many populations [9,10].
Polymorphisms in genes that encode drug-metabolizing enzymes, in addition to HLA-B alleles, play a key role in the initiation of CADRs by slowing drug metabolism. PHT is metabolized primarily by CYP 450 in phase I reactions and mostly by UDP-glucuronosyl transferase (UGT) in phase II reactions. The enzymes CYP2C9 and CYP2C19 are responsible for 90% and 10% of the metabolism of PHT, respectively [11].In CYP2C9, 2* and 3* variants are responsible for reduced PHT clearance [12,13].
Due to distinct waves of immigration, dissimilar genetic patterns in the Indian population have been extensively documented: North Indians are genetically closer to Caucasians, whilst Central Indians are closer to Asians than the European group [14,15]. The South Indians, or Dravidians, who are distributed on the southern side of India (peninsular India), are the original inhabitants of the Indian sub-continent; thus, the distribution of polymorphic alleles is not homogenous [16,17]. For example, the frequency of HLA B*15:02 is lower (0%) among West Coast Parsi and higher (6%) among Pawra in the Khandesh region [18,19]. Similarly, the frequency of CYP2C9*3 varies by location in India: a higher frequency (7–9%) is seen in the Dravidian population, whereas the prevalence is relatively low (0–2%) in the North Indian population [20].
A few studies in India have reported the association of HLA alleles and AEDs induced CADRs [21,22,23,24], and one study from South India reported a high serum PHT concentration in healthy volunteers carrying CYP2C9*3 [19]. There has been little research into the genetic risk factors for PHT-CADRs in the South Indian Tamils, who inhabit primarily in Tamil Nadu, India, and also scantily in Malaysia, Sri Lanka, Singapore, and Mauritius. We aimed to look at the relationship between HLA-B alleles, CYP2C9*3, and PHT-CADRs in the South Indian Tamil population for the first time in this study.

2. Methodology

2.1. Study Design and Settings

This study was conducted as a prospective case-controlled study over a period of 13-months in the Neurology-Outpatient Department (OPD), Neurology ward, and Neurology-Intensive Care Unit (ICU) of Meenakshi Mission Hospital and Research Center (MMHRC) in Madurai, Tamil Nadu, India. This study was approved by the Institutional Ethics Committee Board of MMHRC (MMHRC/IEC/07/2018). DNA analysis of both cases and tolerant controls was performed at the Immunology Department, School of Biological Sciences, Madurai Kamaraj University, Madurai, Tamil Nadu, India (Figure 1).

2.2. Participants

Patients who had reported CADRs within 12 weeks of using PHT were included in the study. Patients on PHT for more than 3 months and had no signs or symptoms of CADRs were considered tolerant controls. The patients who signed the consent form and agreed to give 3 mL of whole blood were included. Patients experiencing other PHT-ADRs and with established skin problems, such as psoriasis and contact dermatitis, were eliminated. Normal-control data was retrieved only from South Indian populations, but no such restrictions were kept for pooled data analysis.
The dermatologist at MMHRC diagnosed all of the PHT-CADRs as well as additional diagnostic criteria, such as the temporal relationship with phenytoin, clinical morphology of the skin, and mucosal and systemic involvement. SJS/TEN, acute generalized exanthematous pustulosis (AGEP), exfoliative dermatitis (ED), drug rash with eosinophilia and systemic symptoms (DRESS) reactions were considered as SCARs. The remaining reactions with no/less mucosal involvement were considered to be mild-moderate reactions in this study [25]. Maculopapular exanthema was defined as a rash characterized by fine pink macules/papules/lesions on the skin with no mucosal or systemic involvement (MPE) [26]. Acneiform drug eruption was defined as a monomorphic eruption without comedones (AFDE) [27].
Fixed drug eruption (FDE) was defined as a single round and oral, sharply demarcated, red-lined lesion with a diameter of 1–10 cm [28]. PHT-lichenoid drug eruption was defined as skin lesions characterized by scaling and hypertrophic pigmentation, generally in combination with oral eruption (LDE) [29]. Patients with rapidly developing dark purpuric macules, atypical target lesions, blisters accompanied by mucosal and skin detachment were diagnosed with SJS, according to Roujeau’s criteria [30]. The RegiSCAR criteria were used to diagnose DRESS and DHS, which included an acute skin rash with at least one internal organ lymphadenopathy; hematologic abnormalities, such as eosinophilia and atypical lymphocytosis; and fever [25]. Patients with erythematic inflammatory skin disease, scaling on the cutaneous surface of the skin, thickened skin, itching, swollen lymph nodes, fever, and fluid loss were considered as PHT-ED [31].
Once the patient was identified as a case or control, the consent form was handed to them along with a patient information document created in the local language of Tamil. The normal-control data for the HLA-B association study was obtained from Leenam Dedhia et al. (2015), who investigated HLA diversity and its significance in South Indians [32]. The normal control data for CYP2C9*3 and PHT-CADRs association testing is from the published literature (Nahar R et al. 2013) [33].

2.3. Causality Assessment

The causality of PHT-CADRs was assessed using Naranjo’s scale, with patients scoring 9 (definite) and 5–8 (probable) included in this study [34]. Patients reported with DRESS were included if the RegiSCAR score was definite >5 and probable (4–5) [35]. The ALDEN score was used to assess the drug’s causality with SJS: patients scoring very probable >6 or probable 4–5 were included [36].

2.4. Genotyping

2.4.1. DNA Extraction/HLA-B Genotyping

DNA was extracted from 3mL of peripheral blood using the salting-out method [37,38]. A UV spectrophotometer was used to measure the concentration and purity of DNA by measuring its optical density (OD) at 260 nm. The polymerase chain reactions-sequence-specific primer technique was used to genotype HLA-B. (PCR-SSP) (Applied Biosystems Verti-Thermal cycler, Thermo Fisher Scientific, Waltham, MA, USA) [39].

2.4.2. CYP2C9*3 Sequencing

The reference DNA sequence of the target variant CYP2C9*3 was retrieved (rs1057910), and the genomic DNA was amplified using the selected forward primer from chromosome passion at 94981018 to 94981037 bp (GTGCATCTGTAACCATCCTC) and the reverse complementary primer from chromosome passion at 94981455 to 94981476 bp (GAGTTATGCACTTCTCTCACCC). The PCR DNA was purified according to the manufacturer’s protocol by MinElute PCR Purification Kit (Cat. No. 28006, Qiagen, Valencia, CA, USA) and sequenced using a 3500 automatic DNA segmentation analyser (3730 DNA analyzer, Applied Biosystems, Thermo Fisher Scientific, Waltham, MA, USA). Sequence scanner software was used to obtain sequential sequencing.

2.5. Statistical Analysis

To compare the demographics and clinical features of the case and control groups, a Student’s t-test was performed. The Fisher exact or Pearson chi-square tests were used if the demographics are categorical variables. The OR was calculated to see if there was an association between specific pharmacogenetic risk factors and PHT-CADRs. The Fisher exact test was used to get the p-value. Woolf’s logit method was used for any cells in a contingency table that had zero. The bivariate analysis was carried out on alleles with a prevalence of more than 5% in the case group. To examine the association of various risk alleles with PHT-induced CADRs, PGx variants that exhibit significant association in bivariate analysis were included in a multivariate binary logistic regression analysis.
We used pooled data analysis for the HLA B*51:01 to boost the study’s power. The HLA B*51:01 case and tolerant control data were obtained from one of the North Indian studies that reported the relationship between HLA alleles and AED-induced CADRs [20].
All the statistical analysis was performed using GraphPad Prism 8. After Bonferroni correction, a p-value 0.008 (<0.05/6—two-tailed) was considered significant.

3. Results

3.1. Patient Demographics

This study included 30 PHT-tolerant and 25 PHT-induced CADRs epileptic cases in this study. Six (23.07%) of the 25 cases were SCARs, including two cases of SJS, three cases of DRESS, and one case of ED. The remaining 19 cases (76.92%) were mild-moderate reactions, comprised of 15 MPE, 2 AFDE, 1 LDE, and 1 FDE. The case group included 14 males and 11 females with a mean age of 40.60 ± 18.15 years, while the PHT-tolerant group included 18 males and 12 females with a mean age of 36.21 ± 14.71 years. Epilepsy, seizure, cerebrovascular accidents, and CNS infections were the common indications for PHT in both the case and tolerant groups (Table 1) (Supplementary Table S1).
The CYP2C9 normal-control data from the previous study included 82 healthy people (40 males and 42 females) from Dravidian or South Indian populations, such as Tamils (25), Andhra Pradesh (32), and Kerala [35]. This study included the data of 463 HLA-B normal-controls (Tamils) that was retrieved from past literature [34]. A total of 52 PHT-CADRs patients (31 MPE, 7 SJS/TEN, 8 DRESS, and 5 FDE) and 100 PHT-tolerant control patients with ages ranging from 6 to 72 were included in the pooled analysis. There were 22/30 and 41/59 females and males in their case and tolerant groups, respectively [20].

3.2. Clinical Features of PHT-CADRs

The itching was more common among mild-moderate reactions (13/20) than SCARs (1/6). Maculopapular rash/exanthematous rash or lesions/skinredness/burning sensation or warmth while touching were the most common cutaneous clinical manifestations in mild-moderate reactions (6/6 in SCARs and 13/20). Papules, pustules, blisters, and erythema were other serious features of skin reactions that were more common with SCARs (4/6) than mild to moderate reactions (5/20). The most commonly impacted mucosal sites were the mouth, eyes, genitals, and anogenital mucosa. SCARs were the only ones that showed systemic involvement. The most commonly affected systems were the liver and hematological systems. All three DRESS patients showed lymphadenopathy, abnormal lymphocytes, and eosinophilia. In one ED patient, neutrophilic leukocytosis was observed (Table 1).
The onset latency period ranged from 7 to 42 days, with a mean of 21.7 days. In Naranjo’s causation assessment, all cases received a likely score of ≥5. According to the RegiSCAR and ALDEN criteria, all DRESS and SJS patients had a definite (>5) and very probable (>6) connection.

3.3. Frequency of HLA-B Alleles in PHT-Cases and Tolerant Controls

HLA-B genotyping data from PHT-CADRs cases showed higher frequencies (>5%) of HLA B*40:01 (40%), HLA B*55:01 (20%), HLA B*51:01 (18%), and HLA B*07:02 (10%) alleles and lower frequencies (<5) of HLA B*57:01 (4%), HLA B*52:01(2%), HLA B*15:01 (2%), HLA B*07:01, and HLA B*35:01(2%).In tolerant-controls, the following alleles were observed more frequently (>5%): HLA B*40:01 (28.33%), HLA B*55:01 (11.66%), HLA B*51:01 (11.66%), and HLA B*07:02 (6.66%), HLA B*52:01(6.66%), HLA B*15:01 (11.66%), and HLA B*35:01 (13.33%). On the other hand, HLA B*07:01 and HLA B*54:01 were reported at lower frequencies (Figure 2 and Supplementary Table S1).

3.4. Bivariate Analysis of HLA B Alleles and PHT-CADRs

3.4.1. SCARs

All six SCARs patients carried HLA B*51:01 allele against 7 out of 30 PHT-tolerant controls. In sub-group analysis, the association between PHT-DRESS was well demonstrated with both tolerant (OR = 21.93, 95% CI 1.013–474.9, p = 0.022) and normal controls (OR = 42.54, 95% CI 2.114–855.8, 0.003). The association between PHT-SJS and HLA B*51:01 was marginal with the PHT-tolerant control group (p = 0.072) and stronger with normal control group (p = 0.021), respectively. HLA B*40:01 was reported in all DRESS (3) and ED (1) cases against 14 out of 30 tolerant patients (28.33%), three of which were found to be homozygous. The association between PHT/DRESS and HLA B*40:01 was insignificant as the comparison made with PHT-tolerant group (p = 0.227) but found to be significant as the cases compared with normal control group (p = 0.0001). In this study, no other HLA B alleles were identified to be related to PHT-SCARs. The HLA B*51:01 genetic marker to PHT-SCAR has a significant positive/negative predictive value (46/100).

3.4.2. Mild-Moderate Reactions

The HLA B*40:01 allele was detected in 14 patients (73.68%) in the mild-moderate reaction group, and two PHT-MPE patients were homozygous for it. When compared to normal control, the association between HLA B*40:01/PHT-mild-moderate reactions, particularly MPE (OR = 10.45, 95% CI 3.299–29.34, p = 0.0001) and AFDE (OR = 25.65, 95% CI 1.172–552.8, p = 0.028), was found to be significant. The HLA B*55:01 allele was the second most common among the mild-moderate reaction group. Only PHT-MPE patients had a greater frequency of this allele (33.33%), and homozygosity was found in one patient. The association between HLA B*55:01 and PHT/MPE was stronger in both PHT-tolerant (OR = 4.929, 95% CI 1.322–17.33, p = 0.022*) and normal control group (OR = 204.0, 95% CI 27.72–2229, p < 0.0001*) with high positive/negative predictive value (56/79).

3.4.3. Overall PHT-Induced CADRs

HLA B*51:01, HLA B*40:01, and HLA B*55:01 alleles were shown to have a significant association with overall PHT-CADRs when compared to normal controls, (HLA B*51:01; OR = 3.43, 95% CI 1.422–8.652, p = 0.017; HLA B*40:01, OR = 13.44, 95% CI 4.849–33.14, p = <0.0001; HLA B*55:01;OR = 76.50, 95% CI 10.40–841.9, p < 0.0001), but none of them were found to have a significant association with PHT-overall CADRs when compared to PHT-tolerant controls. The PHT-tolerant group showed higher frequencies of HLA B*15:01 and HLA B*35:01 than the case and normal control groups. For these alleles, the negative association with PHT-CADRs was found to be significant. (HLA B*15:01; OR = 0.1369, 95% CI of 0.0117–0.9204, p = 0.0593 and HLA B*35:01; OR = 0.1146, 95% CI 0.0099–0.7273, p = 0.0307) (Table 2).

3.4.4. CYP2C9*3

In the current study, the CYP2C9*3 (AC) genotype was found in 12 out of 25 PHT-CADRs cases, but only two patients in the PHT-tolerant group had this allele, and the homozygous (CC) genotype was not found in either case or the tolerant control group (Table 3). The association testing was performed between case vs. PHT-tolerant and cases vs normal-control [23]. In both comparisons, the analysis revealed a stronger association between CYP2C9*3 and overall PHT-CADRs (case vs. PHT-tolerant; OR = 12.92: 95% CI 2.777–61.46, p = 0.0006, case vs. normal control; OR = 5.385, 95% CI 1.917–13.67 p = 0.0017). SCARs (OR = 26.00, 95% CI2.855–1720, p = 0.0043) and mild-moderate reactions (OR = 9.455; 95% CI1.628–47.25, p = 0.0086) both showed a positive susceptibility association with CYP2C9*3 in sub-group analysis. All three DRESS patients had this mutant allele in SCARs, and the relationship was found to be significant (OR = 74.20, 95% CI 2.922–1884, p = 0.002). This allele was observed in 5 MPE, 1 AFDE, 1 LDE, and 1 FDE patients with mild-moderate reactions. The MPE was found to have a substantial association withCYP2C9*3 (OR = 6.500, 95% CI 1.008–35.01, p = 0.039) (Table 2).

3.5. Pooled Data Analysis

Prior to pooling, the association between HLA B*51:01 and PHT-SCARs was shown to be stronger with broad CI (OR = 40.73, 95% CI 2.045–811.3, p = 0.0009) but not with overall PHT-CADRs (OR = 0.848, 95% CI 0.5871–0.377, p = 0.377).The number of CADRs increased to 75 after pooling, with 21 patients carrying HLA B*51:01, while the number of tolerant controls climbed to 130, with 15 HLA B*51:01 carriers.HLA B*51:01 and PHT-overall CADRs were found to have a significant relationship (OR = 2.323, 95% CI1.122–5.899, p = 0.037). Even after pooling the data, the connection between HLA B* 51:01 and mild-moderate reactions is weaker (p = 0.07) (Table 3).

3.6. Multivariate Binary Logistic Regression Analysis

The multivariate binary logistic regression analysis showed a stronger association between CYP2C9*3 and PHT-induced all kinds of CADRs, and it was the only predictor that was related to both severe and mild-moderate reactions (OR = 4.041, 95% CI 1.125–15.6, p = 0.035). In subgroup analysis, PHT-MPE was shown to be highly associated with HLA B*55:01 (OR = 12.00, 95%CI 2.759–84.82, p = 0.003), and CYP2C9*3 was the sole predictor variable with a significant association with PHT-SCARs (OR = 12.00, 95%CI 2.759–84.82, p = 0.003) (Table 4).

4. Discussion

In Asian populations, HLA B*51:01 has been related to PHT-induced CADRs, particularly with PHT-SCARs [40,41].A recent multicenter East Asian study that investigated the genetic predictors of PHT-hypersensitivity reported that concurrent testing of HLA B* 13:01/HLA B* 15:02/HLA B*51:01 and CYP2C9*3 would help in identifying individuals at risk of developing PHT-CADRs [42]. In India, one of the North Indian studies reported HLA B*57:01 (OR = 11.00, 95% CI: 1.41–85.81, p = 0.05) and HLA B*51:01 (OR = 6.90, 95% CI: 1.38–34.29, p = 0.007) as risk factors for CBZ-SJS and PHT-DRESS, respectively [20]. In the current study, the HLA B*51:01 allele was found to be strongly associated with PHT-SCARs, especially DRESS, and the pooled data analysis corroborated this relationship with both SCARs and overall CADRs. However, this allele was presented in all of our SCAR patients (6/6), whereas only 4 of 15 patients in pooled data had this allele, which could be attributable to distribution variations. In Mumbai (central India), for example, it is more frequent among Patels (19.60%) and Iyers (17.60%), whereas it is less common among Marathas (4.84%) [43] and North Indian Hindus (3.5%) (Lucknow) [44]. Its distribution among South Indians, particularly Tamils, ranges between 8–12.5%.
A study from South India reported that HLA B*07 was the most common allele (6–13%) in the HLA-B gene, and its association with cervical cancer along with HLA DQ8 was found to be significant. The next common alleles reported in South Indians were the split antigens of the HLA B5 serotype, HLA B* 51 (8–12.5%) and HLA B* 52 (5–10%). Their association with different vasculitides was reported (HLA B*51 and Behchet’sdiseases, HLA B* 52, and Takayasu’s arteritis) [45].
In this study, the association between HLA B*40:01 and PHT-mild moderate reactions in particular, MPE, was found to be stronger when compared to normal healthy controls, confirming previous findings (Sukesm et al., 2020) [46], which confirmed HLA B*40:01 as a risk factor for PHT-induced MPE (OR 3.647; 95% CI, 1.193–11.147; p = 0.023).This allele could be a drug-specific HLA genetic marker for PHT-MPE. However, a study with a larger cohort is needed to confirm this finding.
HLA B*55:01 was not shown to be susceptible to PHT-induced CADRs in any previous association studies that looked into the relationship between HLA and AEDs. The current study is the first to confirm this association in the South Indian-Tamil population. However, a Han-Chinese study also showed a correlation between HLA B*55:01 and LTG-induced MPE (OR = 24.78, 95% CI 1.50–408.76, p = 0.08) [47]. This genotype has been linked to penicillin hypersensitivity and nevirapine-induced SCARs in addition to AEDs [48,49]. In this study, HLA B*15:01 and HLA B*35:01 were found in a higher percentage of PHT-tolerant people than in cases and normal controls. This finding is consistent with other Indian studies that indicated a higher prevalence of the HLA B*15:01 allele in the control group than in the case group [23,24], whereas HLA B*35:01 was associated with LTG-induced MPE [50]. In the Tamil population, HLA B*15:01 and HLA B*35:01 may be protective alleles for PHT-CADRs.
The mutant allele (AC) was present in 48 and 7.14 percent of cases and tolerant-control groups, respectively, with no homozygosity (CC), which is similar to a study that found no frequency of CC genotype in the South Indian Dravidian community [19]. The current investigation found a substantial link between the CYP2C9*3 heterozygous condition and PHT-induced CADRs. In this investigation, patients with CYP2C9*3 alleles were almost 13 times more vulnerable to PHT-CADRs than the tolerant group. This finding is comparable to that of a Thai study, which found that patients with CYP2C9*3 have a 14.5 times higher incidence of PHT-SCARs. A subgroup analysis of this study within the SCARs group reveals a higher correlation between CYP2C9*3 and PHT-induced DRESS, which was also corroborated in another Thai investigation [46,51,52].
There may be some limitations to this research. We did not rule out patients who were using CYP2C9 inhibitors, which could have contributed to the rise in PHT levels. Despite this, only three patients were prescribed VPA (CYP2C9 inhibitor), and no other known CYP2C9 inhibitors were prescribed in this group. Although a few studies have found a relationship between CYP2C9*2 and PHT-CADRs, we did not investigate this allele in this study because it is extremely rare (1–2%) in the South Indian population. Another limitation of our research is the small sample size. The rare outcome of interest is the reason for it. In addition to genetic defects, clinical and non-clinical factors may play a role in the initiation of PHT-CADRs, and these should be examined alongside genetic variants.

5. Conclusions

CYP2C9*3 and HLA B*51:01 were found to be associated with PHT-SCARs and PHT-DRESS. On the other hand, PHT-mild/moderate cutaneous reactions are linked with HLA B*55:01 and HLA B*40:01 in this study. This is the first study in South India, specifically among Tamils, to show a correlation between HLA B*51:01, HLA B*55:01, and CYP2C9*3 alleles and PHT-CADRs. These alleles can be employed as genetic markers to identify individuals who are susceptible to PHT-CADRs and to ensure that PHT is as safe as possible for Tamil epileptic patients. Furthermore, our findings highlight the necessity of including the HLA B*5101 and CYP2C9*3 alleles into a pre-emptive genetic testing panel for Asians with PHT-CADRs.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/jpm11080737/s1, Figure S1: Frequently observed HLA B alleles among PHT-CADRs and PHT Tolerant; Table S1 Clinical characteristics and HLA-B and CYP2C9*3 genotyping of the patients with PHT-CADRs.

Author Contributions

S.J. was responsible for concept, study design, data collection, statistical analysis, and drafting and revising the manuscript. S.P. was responsible for concept and study design, supervised the study, revised the manuscript content, and approved the final manuscript. K.B. contributed to phenotype data of the study participants, laboratory consultation, and critically revised the manuscript content. C.S. contributed to study design, and supervised the study. T.C.V.A. contributed samples and case diagnosis. B.C. contributed to study design and revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was done as part of a PhD study that was funded by the Graduate School, Prince of Songkla University, Hat Yai, Thailand-90110.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board (or Ethics Committee) of MEENAKSHI MISSION HOSPITAL AND RESEARCH CENTER (MMHRC/IEC/07/2018).

Informed Consent Statement

Informed consent was obtained from all subjects.

Acknowledgments

We would like to thank the staff members of the Neurology Department, Madurai Meenakshi Mission Hospital and Research Center, and staff members of the Department of Immunology, School of Biological Sciences, Madurai Kamaraj University Madurai, TN, India, for their assistance in conducting the study.

Conflicts of Interest

All authors declare no conflicts of interest.

Ethics Approval:

The study protocol was approved by the Institutional Ethical Committee of the Meenakshi Mission Hospital and Research Center. (MMHRC/IEC/07/2018).

Abbreviations

ALDENAlgorithm of drug causality for epidermal necrolysis
AEDsAnti-epileptic drugs
CBZCarbamazepine
PHTPhenytoin
DRESSDrug reaction with eosinophilia and systemic symptoms
HLAHuman leukocyte antigen
MPEMaculopapular exanthema
NPVNegative predictive values
SCARsSevere cutaneous adverse drug reactions
EDExfoliative dermatitis
AFDEAcneiform drug eruption
FDEFixed drug eruption
LDELichenoid drug eruption
CADRsCutaneous adverse drug reactions
VPAValproic acid
TENToxic epidermal necrolysis

References

  1. Twardowschy, C.A.; Germiniani, F.; Werneck, L.C.; Silvado, C.; De Paola, L. Pearls & Oysters: Soft-tissue necrosis as a result of intravenous leakage of phenytoin. Neurology 2009, 73, e94–e95. [Google Scholar] [PubMed]
  2. Glick, T.H.; Workman, T.P.; Gaufberg, S.V. Preventing phenytoin intoxication: Safer use of a familiar anticonvulsant. J. Fam. Pr. 2004, 53, 197–202. [Google Scholar]
  3. Błaszczyk, B.; Szpringer, M.; Czuczwar, S.J.; Lasoń, W. Single centre 20 years survey of antiepileptic drug-induced hypersensitivity reactions. Pharmacol. Rep. 2013, 65, 399–409. [Google Scholar] [CrossRef]
  4. Chatterjee, S.; Ghosh, A.; Barbhuiya, J.; Dey, S. Adverse cutaneous drug reactions: A one year survey at a dermatology outpatient clinic of a tertiary care hospital. Indian J. Pharmacol. 2006, 38, 429–431. [Google Scholar] [CrossRef]
  5. Patel, T.; Barvaliya, M.; Sharma, D.; Tripathi, C. A systematic review of the drug-induced Stevens-Johnson syndrome and toxic epidermal necrolysis in Indian population. Indian J. Dermatol. Venereol. Leprol. 2013, 79, 389–398. [Google Scholar] [CrossRef]
  6. Chung, W.H.; Hung, S.I.; Hong, H.S. Medical genetics: A marker for Stevens-Johnson syndrome. Nature 2004, 428, 486. [Google Scholar] [CrossRef] [PubMed]
  7. Karnes, J.H.; Rettie, A.E.; Somogyi, A.A.; Huddart, R.; Fohner, A.E.; Formea, C.M.; Lee, M.T.M.; Llerena, A.; Whirl-Carrillo, M.; Klein, T.E.; et al. Clinical Pharmacogenetics Implementation Consortium (CPIC) guideline for CYP2C9 and HLA-B genotypes and phenytoin dosing: 2020 Update. Clin. Pharmacol. Ther. 2020, 109, 302–309. [Google Scholar] [CrossRef]
  8. Sukasem, C.; Sririttha, S.; Chaichan, C.; Nakkrut, T.; Satapornpong, P.; Jaruthamsophon, K.; Jantararoungtong, T.; Koomdee, N.; Medhasi, S.; Oo-Puthinan, S.; et al. Spectrum of cutaneous adverse reactions to aromatic antiepileptic drugs and human leukocyte antigen genotypes in Thai patients and meta-analysis. Pharm. J. 2021, 1–9. [Google Scholar] [CrossRef]
  9. Locharernkul, C.; Loplumlert, J.; Limotai, C.; Korkij, W.; Desudchit, T.; Tongkobpetch, S.; Kangwanshiratada, O.; Hirankarn, N.; Suphapeetiporn, K.; Shotelersuk, V. Carbamazepine and phenytoin induced Stevens-Johnson syndrome is associated with HLA-B*1502 allele in Thai population. Epilepsia 2008, 49, 2087–2091. [Google Scholar] [CrossRef]
  10. Chang, C.-C.; Ng, C.-C.; Too, C.-L.; Choon, S.-E.; Lee, C.-K.; Chung, W.-H.; Hussein, S.H.; Lim, K.-S.; Murad, S. Association of HLA-B*15:13 and HLA-B*15:02 with phenytoin-induced severe cutaneous adverse reactions in a Malay population. Pharm. J. 2016, 17, 170–173. [Google Scholar] [CrossRef]
  11. Thorn, C.F.; Whirl-Carrillo, M.; Leeder, J.S.; Klein, T.E.; Altman, R.B. PharmGKB summary: Phenytoin pathway. Pharm. Genom. 2012, 22, 466–470. [Google Scholar] [CrossRef] [PubMed]
  12. Chung, W.-H.; Chang, W.-C.; Lee, Y.-S.; Wu, Y.-Y.; Yang, C.-H.; Ho, H.-C.; Chen, M.-J.; Lin, J.-Y.; Hui, R.C.-Y.; Ho, J.-C.; et al. Genetic variants associated with phenytoin-related severe cutaneous adverse reactions. J. Am. Med. Assoc. 2014, 312, 525–534. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Rosemary, J.; Surendiran, A.; Rajan, S.; Shashindran, C.H.; Adithan, C. Influence of the CYP2C9 AND CYP2C19 polymorphisms on phenytoin hydroxylation in healthy individuals from south India. IndianJ. Med. Res. 2006, 123, 665–670. [Google Scholar]
  14. Agrawal, S.; Srivastava, S.K.; Borkar, M.; Chaudhuri, T.K. Genetic affinities of north and northeastern populations of India: Inference from HLA-based study. Tissue Antigens 2008, 72, 120–130. [Google Scholar] [CrossRef]
  15. Joshi, C.G.; Patel, J.S.; Patel, M.M.; Koringa, P.G.; Shah, T.M.; Patel, A.K.; Tripathi, A.K.; Mathew, A.; Rajapurkar, M.M. Human leukocyte antigen alleles, genotypes and haplotypes frequencies in renal transplant donors and recipients from West Central India. Indian J. Hum. Genet. 2013, 19, 219–232. [Google Scholar] [CrossRef] [Green Version]
  16. Coon, C.S. Human populations. The New Encyclopedia Britannica, 15th ed.; Encyclopaedia Britannica: Chicago, IL, USA, 1983; Volume 14, pp. 839–848. [Google Scholar]
  17. Mittal, K.K.; Naik, S.; Sansonetti, N.; Cowherd, R.; Kumar, R.; Wong, D.M. The HLA antigens in Indian Hindus. Tissue Antigens 1982, 20, 223–226. [Google Scholar] [CrossRef] [PubMed]
  18. Shankarkumar, U.; Hla, A.B.C. Allele frequencies in a Parsi population from Western India. Hum. Immunol. 2004, 65, 992–993. [Google Scholar] [CrossRef]
  19. Shankarkumar, U.; Hla, A.B.C. Allele frequencies in a Pawra population from Khandesh, India. Hum. Immunol. 2004, 65, 958. [Google Scholar] [CrossRef]
  20. Nizamuddin, S.; Dubey, S.; Singh, S.; Sharma, S.; Mishra, A.; Harish, K.; Joshi, J.; Thangaraj, K. Genetic polymorphism of Cytochrome-P450-2C9 (CYP2C9) in Indian populations. bioRxiv 2017. [Google Scholar] [CrossRef] [Green Version]
  21. Ihtisham, K.; Ramanujam, B.; Srivastava, S.; Mehra, N.K.; Kaur, G.; Khanna, N.; Jain, S.; Kumar, S.; Kaul, B.; Samudrala, R.; et al. Association of cutaneous adverse drug reactions due to antiepileptic drugs with HLA alleles in a North Indian population. Seizure 2019, 66, 99–103. [Google Scholar] [CrossRef] [Green Version]
  22. Devi, K. The association of HLA B*15:02 allele and Stevens-Johnson syndrome/toxic epidermal necrolysis induced by aromatic anticonvulsant drugs in a South Indian population. Int. J. Dermatol. 2017, 57, 70–73. [Google Scholar] [CrossRef]
  23. Ramanujam, B.; Ihtisham, K.; Kaur, G.; Srivastava, S.; Mehra, N.K.; Khanna, N.; Singh, M.; Tripathi, M. Spectrum of cutaneous adverse reactions to Levetiracetam and human Leukocyte antigen typing in North-Indian patients. J. Epilepsy Res. 2016, 6, 87–92. [Google Scholar] [CrossRef] [PubMed]
  24. Aggarwal, R.; Sharma, M.; Modi, M.; Garg, V.K.; Salaria, M. HLA-B*15 02 is associated with carbamazepine induced Stevens-Johnson syndrome in North Indian population. Hum. Immunol. 2014, 75, 1120–1122. [Google Scholar] [CrossRef]
  25. Patel, T.K.; Thakkar, S.H.; Sharma, D.C. Cutaneous adverse drug reactions in Indian population: A systematic review. Indian Dermatol. Online J. 2014, 5, S76–S86. [Google Scholar] [CrossRef] [PubMed]
  26. Shi, Y.W.; Min, F.L.; Zhou, D.; Qin, B.; Wang, J.; Hu, F.Y.; Cheung, Y.K.; Zhou, J.H.; Hu, X.S.; Zhou, J.Q.; et al. HLA-A*24:02 as a common risk factor for antiepileptic drug-induced cutaneous adverse reactions. Neurology 2017, 88, 2183–2191. [Google Scholar] [CrossRef] [PubMed]
  27. Du-Thanh, A.; Kluger, N.; Bensalleh, H.; Guillot, B. Drug-induced acneiform eruption. Am. J. Clin. Dermatology 2011, 12, 233. [Google Scholar] [CrossRef]
  28. Rita, V.; Nilofar, D.; Nidhi, J.; Rochit, S. A study of adverse cutaneous drug reactions (ACDR) owing to antiepileptics at a rural-based tertiary-care center, Gujarat. Natl. J. Physiol. Pharm.Pharmacol. 2016, 6, 140. [Google Scholar]
  29. Roji, M.; Sebastian, M.; Lucca, J.M.; Pss, R. Phenytoin induced oral Lichenoid eruption and melasma: A case report. Indian J. Pharm. Pr. 2018, 11, 55–57. [Google Scholar] [CrossRef] [Green Version]
  30. Roujeau, J.C. The spectrum of Stevens-Johnson syndrome and toxic epidermal necrolysis: A clinical classification. J.Invest. Dermatol. 1994, 102, 28S–30S. [Google Scholar] [CrossRef]
  31. Voora, L.; Shastry, C.; Bhandari, R.; Sukeerthi, D.; Rawal, K.B.; Chand, S. Phenytoin-induced Erythroderma. J. Young Pharm. 2019, 11, 320–321. [Google Scholar] [CrossRef] [Green Version]
  32. Dedhia, L.; Gadekar, S.; Mehta, P.; Parekh, S. HLA haplotype diversity in the South Indian population and its relevance. Indian J. Transplant. 2015, 9, 138–143. [Google Scholar] [CrossRef]
  33. Nahar, R.; Deb, R.; Saxena, R.; Puri, R.D.; Verma, I.C. Variability in CYP2C9 allele frequency: A pilot study of its predicted impact on warfarin response among healthy South and North Indians. Pharmacol. Rep. 2013, 65, 187–194. [Google Scholar] [CrossRef]
  34. Naranjo, C.; Busto, U.; Sellers, E.; Sandor, P.; Ruiz, I.; Roberts, E.; Janecek, E.; Domecq, C.; Greenblatt, D.J. A method for estimating the probability of adverse drug reactions. Clin. Pharmacol. Ther. 1981, 30, 239–245. [Google Scholar] [CrossRef] [PubMed]
  35. Kardaun, S.H.; Sidoroff, A.; Valeyrie-Allanore, L.; Halevy, S.; Davidovici, B.B.; Mockenhaupt, M.; Roujeau, J.C. Variability in the clinical pattern of cutaneous side-effects of drugs with systemic symptoms: Does a DRESS syndrome really exist? Br.J. Dermatol. 2007, 156, 609–611. [Google Scholar] [CrossRef] [PubMed]
  36. Sassolas, B.; Haddad, C.; Mockenhaupt, M.; Dunant, A.; Liss, Y.; Bork, K.; Haustein, U.F.; Vieluf, D.; Roujeau, J.C.; Le Louet, H. ALDEN, an algorithm for assessment of drug causality in Stevens—Johnson syndrome and toxic epidermal necrolysis: Comparison with case-control analysis. Clin. Pharmacol. Ther. 2010, 88, 60–68. [Google Scholar] [CrossRef]
  37. Welsh, K.; Bunce, M. Molecular typing for the MHC with PCR-SSP. Rev. Immunogenet. 1999, 1, 157–176. [Google Scholar]
  38. Bunce, M.; O’neill, C.M.; Barnardo, M.C.N.M.; Krausa, P.; Browning, M.J.; Morris, P.J.; Welsh, K.I. Phototyping: Comprehensive DNA Typing for HLA A, B, C, DRB1, DRB3, DRB4, and DRB5 & DQB1 by PCR with 144 primers mixes utilizing sequence specific primers (PCR-SSP). Tissue Antigens 1995, 46, 355–367. [Google Scholar]
  39. Tonks, S.; Marsh, S.G.E.; Bunce, M.; Bodmer, J.G. Molecular typing for HLA class I using ARMS-PCR: Further developments following the 12th international histocompatibility workshop. Tissue Antigens 1999, 53, 175–183. [Google Scholar] [CrossRef]
  40. Kaniwa, N.; Sugiyama, E.; Saito, Y.; Kurose, K.; Maekawa, K.; Hasegawa, R.; Furuya, H.; Ikeda, H.; Takahashi, Y.; Muramatsu, M.; et al. Specific HLA types are associated with antiepileptic drug-induced Stevens-Johnson syndrome and toxic epidermal necrolysis in Japanese subjects. Pharmacogenomics 2013, 14, 1821–1831. [Google Scholar] [CrossRef]
  41. Manuyakorn, W.; Likkasittipan, P.; Wattanapokayakit, S.; Suvichapanich, S.; Inunchot, W.; Wichukchinda, N.; Khongkhatithuml, C.; Thampratankul, L.; Kamchaisatian, W.; Benjaponpitak, S.; et al. Association of HLA genotypes with phenytoin induced severe cutaneous adverse drug reactions in Thai children. Epilepsy Res. 2020, 162, 106321. [Google Scholar] [CrossRef]
  42. Su, S.-C.; Chen, C.-B.; Chang, W.-C.; Wang, C.-W.; Fan, W.-L.; Lu, L.-Y.; Nakamura, R.; Saito, Y.; Ueta, M.; Kinoshita, S.; et al. HLA alleles and CYP2C9*3 as predictors of Phenytoin hypersensitivity in East Asians. Clin. Pharmacol. Ther. 2018, 105, 476–485. [Google Scholar] [CrossRef]
  43. Shankarkumar, U.; Ghosh, K.; Colah, R.B.; Gorakshakar, A.C.; Gupte, S.C.; Mohanty, D. HLA antigen distribution in selected caste groups from Mumbai, Maharastra, India. J. Hum. Ecol. 2002, 13, 209–215. [Google Scholar] [CrossRef]
  44. Subramanian, V.S.; Selvaraj, P.; Narayanan, P.R.; Prabhakar, R.; Damodaran, C. Distribution of HLA (class I and class II) antigens in the native Dravidian Hindus of Tamil Nadu, south India. Gene Geogr. Comput. Bull. Hum. Gene Freq. 1995, 9, 15–24. [Google Scholar]
  45. Narayan, S.; Maiers, M.; Halagan, M.; Sathishkannan, A.; Naganathan, C.; Madbouly, A.; Periathiruvadi, S. Human leukocyte antigen (HLA)-A, -B, -C, -DRB1 and -DQB1 haplotype frequencies from 2491 cord blood units from Tamil speaking population from Tamil Nadu, India. Mol. Biol. Rep. 2018, 45, 2821–2829. [Google Scholar] [CrossRef] [PubMed]
  46. Sukasem, C.; Sririttha, S.; Tempark, T.; Klaewsongkram, J.; Rerkpattanapipat, T.; Puangpetch, A.; Boongird, A.; Chulavatnatol, S. Genetic and clinical risk factors associated with phenytoin-induced cutaneous adverse drug reactions in Thai population. Pharmacoepidemiol. Drug Saf. 2020, 29, 565–574. [Google Scholar] [CrossRef]
  47. Shi, Y.-W.; Min, F.-L.; Liu, X.-R.; Zan, L.-X.; Gao, M.-M.; Yu, M.-J.; Liao, W.-P. Hla-B alleles and lamotrigine-induced cutaneous adverse drug reactions in the Han Chinese population. Basic Clin. Pharmacol. Toxicol. 2011, 109, 42–46. [Google Scholar] [CrossRef] [PubMed]
  48. Krebs, K.; Bovijn, J.; Lepamets, M.; Censin, J.C.; Jürgenson, T.; Särg, D.; Abner, E.; Laisk, T.; Luo, L.; Skotee, L.; et al. Genome-wide study identifies association between HLA-B*55:01 and penicillin allergy. Am. J. Hum. Genet. 2020, 107, 612. [Google Scholar] [CrossRef] [PubMed]
  49. Pavlos, R.; McKinnon, E.J.; Ostrov, D.A.; Peters, B.; Buus, S.; Koelle, D.; Chopra, A.; Schutte, R.; Rive, C.; Redwood, A.; et al. Shared peptide binding of HLA Class I and II alleles associate with cutaneous nevirapine hypersensitivity and identify novel risk alleles. Sci. Rep. 2017, 7, 8653. [Google Scholar] [CrossRef]
  50. Peña-Balderas, A.M.; López-Revilla, R. Pharmacogenetics of adverse cutaneous reactions to lamotrigine. Rev. Mex. Neuroci. 2019, 20, 200–206. [Google Scholar] [CrossRef] [Green Version]
  51. Tassaneeyakul, W.; Prabmeechai, N.; Sukasem, C.; Kongpan, T.; Konyoung, P.; Chumworathayi, P.; Tiamkao, S.; Khunarkornsiri, U.; Kulkantrakorn, K.; Saksit, N.; et al. Associations between HLA class I and cytochrome P450 2C9 genetic polymorphisms and phenytoin-related severe cutaneous adverse reactions in a Thai population. Pharmacogenet. Genomics. 2016, 26, 225–234. [Google Scholar] [CrossRef]
  52. Yampayon, K.; Sukasem, C.; Limwongse, C.; Chinvarun, Y.; Tempark, T.; Rerkpattanapipat, T.; Kijsanayotin, P. Influence of genetic and non-genetic factors on phenytoin-induced severe cutaneous adverse drug reactions. Eur. J. Clin. Pharmacol. 2017, 73, 855–865. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Flowchart of study design.
Figure 1. Flowchart of study design.
Jpm 11 00737 g001
Figure 2. Frequently observed HLA B alleles among PHT-CADRs and PHT-tolerant controls.
Figure 2. Frequently observed HLA B alleles among PHT-CADRs and PHT-tolerant controls.
Jpm 11 00737 g002
Table 1. Patient demographics and clinical manifestations of PHT-CADRs and PHT-tolerant group.
Table 1. Patient demographics and clinical manifestations of PHT-CADRs and PHT-tolerant group.
ParametersPHT-Induced CADRsN (%)Total
(25)
Tolerant
(30)
p Value
Severe (6)MPE and other Mild-Moderate Reactions c (19)
Age<2002 (33.33)1 (5.26)3 (12)5 (16.66)0.71
21–4001 (16.66)8 (42.10)9 (36)11 (36.66)0.80
41–602 (33.33)6 (31.57)8 (32)12 (40)0.26
>602 (33.33)3 (15.78)5 (20)02 (6.66)0.12
Mean age ± SD 44.5 ± 21.7739.75 ± 1.41442.15 ± 3.3536.21 ± 14.710.36 a
No of comorbidities (Mean) 2 ± 0.001 ± 1.051.54 ± 0.82037 ± 0.750.01
GenderMale04 (66.66)10 (52.63)14 (56)18 (60)0.59 b
Female02 (33.33)0911 (44)12 (40)
Onset latency (days) 9–267–427–42--
IndicationsSeizure3 (50)9 (47.36)12 (48)11 (36.66)-
Epilepsy2 (33.33)8 (42.10)10 (40)16 (53.33)-
CVA1 (16.66)1 (5.26)02 (8)03 (10)-
Others01 (5.26)01 (4)--
Social HistoryAlcohol (Yes)1 (16.66)6 (31.57)7 (28)8 (26.66)>0.99
Smoking (Yes)03 (15.78)3 (12)1 (3.33)0.31
Causality scoresNaranjo’s75–75–7--
RegiSCARS>5->5--
ALDEN7–8 7–8--
CM/SkinItching1 (16.66)13 (68.42)14 (56)--
Pap. Rash/Redness6 (100)19 (100)25 (25)--
Papules/Pustules4 (66.66)04 (21.05)08 (32)--
Blistering/Peeling2 (33.33)002 (8)--
Erythema P/P3 (50)003 (12)--
CM/MucosaEye6 (100)2 (10.52)08 (32)--
Oral Mucosa6 (100)006 (24)--
Genital4 (66.66)004 (16)--
Anogenital2 (33.33)002 (8)--
F.Edema4 (66.66)2 (10.52)06 (24)--
CM/Systemic
abnormalities
Fever6 (100)3 (10.52)09 (36)--
Liver
AST IU/L (Mean ± SD)SJS/DRESS/ED96.50 ± 9.192/183.7 ± 96.66/185----
ALT IU/L(Mean ± SD)SJS/DRESS/ED72.00 ± 12.73/241.3 ± 139.9/131----
Alk.PhosphataseSJS/DRESS/ED117.5 ± 17.68/331.0 ± 128.7/401----
Haematological abn.
WBC (Mean ± SD)SJS/DRESS/ED5200 ± 424.3/11200 ± 2458/8000/µL----
Lymphocytes (M %)SJS/DRESS/ED17%/56%/-----
Eosinophil (M %)SJS/DRESS/ED09/12/-----
LymphadnopathySJS/DRESS/EDAbsent/3/Absent----
AEDs CombinationPHT+VPA2 (33.33)2 (10.00)4 (15.38)6 (18.75)>0.999
PHT+LEV2 (33.33)5 (25.00)7 (26.92)4 (12.50)0.1935
PHT+VPA+LEV-1 (5.00)1 (03.84)-0.4576
Others2 (33.33)2 (10.00)4 (15.38)0-
SJS, Steven Johnson syndrome; DRESS, drug reactions eosinophilia systemic syndrome; ED, exfoliative dermatitis; MPR, maculo-papular rash; c FDE, fixed drug eruption; LDE, lichenoid drug eruption; AFDE, acneiform drug eruption were in others category, a p value was calculated by independent Student’s t test, b p-value calculated by Fisher’s exact test.
Table 2. Bivariate analyses of association between HLA B alleles and CYP2C9*3 with PHT-CADRs.
Table 2. Bivariate analyses of association between HLA B alleles and CYP2C9*3 with PHT-CADRs.
HLA B AllelesPhenotype
(no of Cases)
Number of CasesCase vs. TolerantCase vs. Healthy PopulationPPV/NPVSensitivity/Specificity
Case
n = 25
Tolerant
n = 30HLA B/28CYP2C9*3
General Population R1 & R2
n = 463/82
OR (95% CI)p Value aOR (95% CI)p Value a
N (%)N (%)N (%)
PHT-SCARs
B*51/51: 01SJS (2)2 (100)7 (23.33)19 (13.86)15.67 (0.6743–64.0)0.0730.38 (1.405–657.1)0.0222/100100/76
DRESS (3)3 (100) 21.93 (1.013–474.9)0.0242.54 (2.114–855.8)0.00330/100100/76
ED (1)1 (100) 9.40 (0.3452–56.0)0.2518.23 (0.716–463.8)0.1412.5/100100/76
B* 40/40:01DRESS (3)3 (100)14 (46.66)22 (16.05)7.96 (0.3787–67.5)0.2235.93 (1.794–719.9)0.00517/100100/53
ED (1)1 (100) 3.41 (0.1288–0.49)0.4815.40 (0.607–390.2)0.166.6/100100/53
B* 5101SCARs (6)6 (100)7 (23.33)19 (13.86)40.73 (2.045–811.3)0.0009296.3 (16.11–5450)<0.00146/100100/76
B* 4001SCARs (6)4 (66.66)14 (46.66)22 (16.05)6.571 (1.157–37.69)0.05740.09 (8.70–213.6)0.000136/9266/76
PHT-Mild-Moderate
B* 40/40:01MPE 15)10 (40)14 (28.33)22 (16.05)2.286( 0.6742–7.887)0.34210.45 (3.299–29.34)<0.00141/7666/53
AFDE (2)2 (50) 5.690 (0.2519–128.5)0.48325.65 (1.172–552.8)0.0212.5/100100/53
LDE (1)1 (50) 3.414 (0.1288–90.49)0.48315.40 (0.607–390.2)0.166.6/100100/53
FDE (1)1 (50) 3.414 (0.1288–90.49)0.48315.40 (0.607–390.2)0.166.6/100100/53
B* 55/55:01MPE (15)9 (33.3)07 (11.66)1 (0.729)4.929 (1.322–17.33)0.022204.0 (27.72–2229)<0.00156/7960/76
B* 40/40:01Mild-Mode CADRs1914 (73.68)14 (46.66)22 (4.75)3.200 (0. 885–10.34)0.08056.13 (18.42–147.5)<0.00150/7073/53
B* 55/55:01Mild-Mode CADRs199 (47.3)07 (23.33)1 (0.21)2.957 (0.922–10.83)0.119415.8 (51.92–454)<0.00156/6947/76
PHT-Over all CADR
B*51/51:01 CADRs9 (36)7 (23.33)19 (13.8)1.848 (0.5871–5.905)0.3773.493 (1.422–8.652)0.0156/5836/76
B*40/40:01 CADRs18 (72)14 (46.66)22 (16.05)2.939 (0.9870–8.842)0.09813.44 (4.849–33.14)<0.00156/6972/53
B*55/55:01 CADRs10 (40)07 (23.33)1 (0.72)1.848 (0.5871–5.905)0.37776.50 (10.40–841.9)<0.00158/6040/76
B*57/57:01 CADRs02 (8)04 (13.33)6 (4.37)0.5652 (0.1013–2.643)0.6771.899 (0.369–8.216)0.3533/538/86
B*52/52:01 CADRs02 (8)04 (13.33)11 (8.02)0.5652 (0.1013–2.643)0.6770.996 (0.209–3.973)>0.9933/538/86
B*07/07:02 CADRs05 (20)04 (13.33)26 (18.9)1.625 (0.3967–5.800)0.7161.067 (0.4067–.909)>0.9955/5620/86
B*15/15:01 CADRs01 (4)07 (23.33)11 (8.02)0.1369 (0.0117–0.9204)0.0590.477(0.042–3.005)0.6912.5/484/76
B*35/35:01 CADRs01 (4)08 (26.66)22 (16.05)0.1146 (0.009–0.7273)0.0300.2178 (0.020–.370)0.2011/474/73
CYP2C9*3 carriers
CYP2C9*3Severe CADRs (6)4 (66.6)2 (7.14)12 (14.63)26.00 (2.855–1726)0.000611.69 (2.386–63.63)0.00966/9266/92
CYP2C9*3DRESS(3)3 (100)2 (7.14)12 (14.63)74.20 (2.922–1884)0.00239.48 (1.920–811.9)0.00460/100100/92
CYP2C9*3M-Moder.
(19)
8 (42.2)2 (7.14)12 (14.63)9.455 (1.628–47.25)0.0082.917 (0.9178–0.969)0.0280/7042/92
CYP2C9*3MPE(15)5 (33.33)2 (7.14)12 (14.63)6.500 (1.008–35.01)0.0395.385 (1.917–13.67)0.1371/7233/92
CYP2C9*3Overall (25)12 (48)2 (7.14)12 (14.63)12.92 (2.777–61.46)0.0044.242 (1.378–12.89)0.00185/6648/92
SJS, Steven Johnson syndrome; DRESS, drug reactions eosinophilia systemic syndrome; ED, exfoliative dermatitis; MPE, maculo-papular eruption; FDE, fixed drug eruption; LDE, lichenoid drug eruption; AFD, acneiform drug eruption; a Fisher’s exact test was used for p calculation, R1 & R2 Percentages of HLA B alleles in Tamil population was taken from Leenam Dedhia et al. 2015 [22], and CYP2C9*3 distribution data was taken from Nahar R et al. 2013 [23] After Bonferroni correction, p-value 0.008 (<0.05/6—two tailed) were considered significant.
Table 3. Pooled data analysis for HLA B*51:01 and PHT-CADRs.
Table 3. Pooled data analysis for HLA B*51:01 and PHT-CADRs.
AllelesPhenotypePresent StudyLiteraturesTotal after PoolingBefore PoolingAfter Pooling PPV/NPV
CaseControlCaseControlCaseControlOR (95%CI)p ValueOR (95%CI)p Value
HLA B*51:01R1SCARs6/67/304/158/10011/2115/13040.73 (2.045–811.3)0.00096.273 (2.24–16.69)0.000846/100 a
47/8 8b
MM3/197/309/358/10012/5415/1300.6161 (0.1563–2.603)0.712.190 (0.939–4.97)0.0730/58 a
44/73 b
Overall CADRs9/257/3013/508/10022/7515/1300.848 (0.5871–0.377)0.372.323 (1.122–5.899)0.0336/76 a
59/68 b
R1 Percentage of HLA B*51:01 in PHT-CADRs and tolerant control data was taken from Ihtisham et al. (2019), a Positive and negative predictive value before pooling, b Positive and negative predictive value after pooling; SCARs, severe cutaneous adverse drug reactions; OR, odds ratio; CI, confidence interval; MM reac, mild-moderate reactions that includes: MPE, maculo-papular Exanthema; FDE, fixed drug eruption; LDE, lichenoid drug eruption; AFDE, acneiform drug eruption.
Table 4. Multivariate binary logistic regression analysis of PHT-CADRs and tolerant-control.
Table 4. Multivariate binary logistic regression analysis of PHT-CADRs and tolerant-control.
Predictor VariablesβSEOR95% CIp Value
PHT- All types of CADRs/tolerant Control
CYP2C9*32.480.83122.759–84.870.003
HLA B*51:010.480.671.610.429–6.1040.47
PHT-SCARs/tolerant Control
CYP2C9*32.521.2212.41.138–136.20.003
HLA B*51:012.041.457.700.447–133.00.16
HLA B*40:01−0.331.340.710.52–9.8950.80
MPE/tolerant Control
HLA B*55:011.390.664.041.125–15.670.03
HLAB*40:010.390.781.470.318–6.8730.61
CYP2C9*31.571.054.890.617–37.790.13
β regression coefficient, SE, standard error; OR, odds ratio; CI, confidence interval;MPE, maculo-papularexanthema.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

John, S.; Balakrishnan, K.; Sukasem, C.; Anand, T.C.V.; Canyuk, B.; Pattharachayakul, S. Association of HLA-B*51:01, HLA-B*55:01, CYP2C9*3, and Phenytoin-Induced Cutaneous Adverse Drug Reactions in the South Indian Tamil Population. J. Pers. Med. 2021, 11, 737. https://doi.org/10.3390/jpm11080737

AMA Style

John S, Balakrishnan K, Sukasem C, Anand TCV, Canyuk B, Pattharachayakul S. Association of HLA-B*51:01, HLA-B*55:01, CYP2C9*3, and Phenytoin-Induced Cutaneous Adverse Drug Reactions in the South Indian Tamil Population. Journal of Personalized Medicine. 2021; 11(8):737. https://doi.org/10.3390/jpm11080737

Chicago/Turabian Style

John, Shobana, Karuppiah Balakrishnan, Chonlaphat Sukasem, Tharmarajan Chinnathambi Vijay Anand, Bhutorn Canyuk, and Sutthiporn Pattharachayakul. 2021. "Association of HLA-B*51:01, HLA-B*55:01, CYP2C9*3, and Phenytoin-Induced Cutaneous Adverse Drug Reactions in the South Indian Tamil Population" Journal of Personalized Medicine 11, no. 8: 737. https://doi.org/10.3390/jpm11080737

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop