# The Absence of CYP3A5*3 Is a Protective Factor to Anticonvulsants Hypersensitivity Reactions: A Case-Control Study in Brazilian Subjects

Luciana Kase Tanno, Daniel Shikanai Kerr, Bernardo dos Santos, Leda Leme Talib, Célia Yamaguti, Helcio Rodrigues, Wagner Farid Gattaz, Jorge Kalil

► To cite this version:

Luciana Kase Tanno, Daniel Shikanai Kerr, Bernardo dos Santos, Leda Leme Talib, Célia Yamaguti, et al.. The Absence of CYP3A5*3 Is a Protective Factor to Anticonvulsants Hypersensitivity Reactions: A Case-Control Study in Brazilian Subjects. PLoS ONE, 2015, 10 (8), pp.e0136141. 10.1371/journal.pone.0136141. hal-03796156

HAL Id: hal-03796156

https://hal.science/hal-03796156

Submitted on 4 Oct 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# The Absence of CYP3A5*3 Is a Protective Factor to Anticonvulsants Hypersensitivity Reactions: A Case-Control Study in Brazilian Subjects

Luciana Kase Tanno¹* , Daniel Shikanai Kerr²,⁵• , Bernardo dos Santos⁴ , Leda Leme Talib²,⁵ , Célia Yamaguti³ , Helcio Rodrigues³ , Wagner Farid Gattaz²,⁵ , Jorge Kalil¹,³

1 Clinical Immunology and Allergy Division, University of São Paulo, São Paulo, Brazil, 2 Laboratory of Neuroscience - LIM-27 Institute of Psychiatry, University of São Paulo, São Paulo, Brazil, 3 Laboratory of Immunology - LIM-19 Clinical Immunology and Allergy Division, University of São Paulo, São Paulo, Brazil, 4 Nursing School of the University of São Paulo, São Paulo, Brazil, 5 Center for Interdisciplinary Research on Applied Neurosciences (NAPNA), University of São Paulo, São Paulo, Brazil

• These authors contributed equally to this work.

* luciana.tanno@gmail.com

<figure><img src="image_2.png" /></figure>

CrossMark
click for updates

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Tanno LK, Kerr DS, dos Santos B, Talib LL, Yamaguti C, Rodrigues H, et al. (2015) The Absence of CYP3A5*3 Is a Protective Factor to Anticonvulsants Hypersensitivity Reactions: A Case-Control Study in Brazilian Subjects. PLoS ONE 10(8): e0136141. doi:10.1371/journal.pone.0136141

**Editor:** Marie-Pierre Dubé, Université de Montréal, CANADA

**Received:** April 2, 2015

**Accepted:** July 30, 2015

**Published:** August 20, 2015

**Copyright:** © 2015 Tanno et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** Data are available as a supplementary table. Further inquiries should be addressed to the authors.

**Funding:** This work received funding from Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP - http://fapesp.br/) under the grant number 2011/22748-1. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

## Abstract

Although aromatic anticonvulsants are usually well tolerated, they can cause cutaneous adverse drug reactions in up to 10% of patients. The clinical manifestations of the antiepileptics-induced hypersensitivity reactions (AHR) vary from mild skin rashes to severe cutaneous drug adverse reactions which are related to high mortality and significant morbidity. Genetic polymorphisms in cytochrome P450 genes are associated with altered enzymatic activity and may contribute to the risk of AHR. Here we present a case-control study in which we genotyped SNPs of CYP2C19, 2C9 and 3A5 of 55 individuals with varying severities of AHR, 83 tolerant, and 366 healthy control subjects from São Paulo, Brazil. Clinical characterization was based on standardized scoring systems and drug patch test. All in vivo investigation followed the ENDA (European Network of Drug Allergy) recommendations. Genotype was determined by real time PCR using peripheral blood DNA as a template. Of all 504 subjects, 65% were females, 45% self-identified as Afro-American, 38% as Caucasian and 17% as having non-African mixed ascendancy. Amongst 55 subjects with AHR, 44 had severe cutaneous drug adverse reactions. Of the 46 drug patch tests performed, 29 (63%) were positive. We found a strong association between the absence of CYP3A5*3 and tolerant subjects when compared to AHR (p = 0.0002, OR = 5.28 [CI95% 2.09–14.84]). None of our groups presented positive association with CYP2C19 and 2C9 polymorphisms, however, both SNPs contributed to separation of cases and tolerants in a Classification and Regression Tree. Our findings indicate that drug metabolism genes can contribute in the tolerability of antiepileptics. CYP3A5*3 is the most prevalent CYP3A5 allele associated with reduced enzymatic function. The current study provides evidence that normal CYP3A5 activity might be a protective factor to aromatic antiepileptics-induced hypersensitivity reactions in Brazilian subjects.
## Introduction

Aromatic anticonvulsants (ARA) are important drugs effective in the treatment of epilepsy, trigeminal neuralgia, and bipolar disorder [1,2]. Although well tolerated by the majority, they can cause cutaneous adverse drug reactions (ADRs) in up to 10% of patients [1]. The clinical manifestations of the antiepileptics-induced hypersensitivity reactions (AHR) vary from mild skin rashes to severe cutaneous drug adverse reactions (SCARs), such as Stevens-Johnson syndrome (SJS), toxic epidermal necrolysis (TEN) and drug reaction with eosinophilia and systemic symptoms (DRESS) [3]. These SCARs are related to high mortality and significant morbidity. The most frequently associated ARA are phenobarbital (PB), phenytoin (PHT), and carbamazepine (CBZ) [4].

The mechanisms of aromatic anticonvulsant-induced SCARs are not fully understood, but evidences suggest an immune-mediated etiology. The discovery of drug-specific T-cells in hypersensitive individuals is consistent with this hypothesis [5]. Drug patch tests (PTs) can reproduce delayed T-cell mediated hypersensitivity to drugs and are useful on determining the offending drug [5,6]. Family and twins studies suggest a genetic predisposition to these reactions [7]. Metabolism of ARA to reactive metabolites is thought to be involved in the pathogenesis of SCARs [8]. The detoxification of ARA occurs mostly in the liver where they are converted in more soluble end products which can usually be secreted by the kidneys. Cytochrome P450 enzymes, such as CYP2C19, CYP2C9 and CYP3A5, are the main enzymes in the hydroxylation of ARA and are expressed mainly in the liver as well as extra-hepatic tissues including the skin [9–13] However, in these metabolic pathways, some reactive intermediate metabolites, such as arene-oxides, can be formed [14,15]. Several Single Nucleotide Polymorphisms (SNPs) have major consequences in the expression and activity of these enzymes, however the results are still not comparable in different populations worldwide [16–18]. These informations are particularly relevant for pharmacogenetics strategies, in which the knowledge of CYP genotype may lead to individualized drug dosing and improved therapeutics.

The association between the CYP2C19, CYP2C9 and CYP3A5 allele frequencies and aromatic anticonvulsant-induced SCARs has not been studied in Brazilian subjects yet. The CYP2C19 and CYP2C9 variants have been described as being risk factors for the development of SCARs in Thai and Japanese populations [16,17]. The CYP3A5*3/*3 genotype was associated with increased half-life of CBZ in African-Americans, but no significant association was observed in Caucasians [18,19]. In addition, recent studies have demonstrated the utility of PTs as a method for the identification of both SJS/TEN and DRESS induced by ARA [20,21]. Therefore, the major goal of this study was to analyze the association between aromatic anticonvulsant-induced hypersensitivity reactions and polymorphisms of CYP2C19, CYP2C9 and CYP3A5 in Brazilian subjects and explore the clinical characteristics and results of PTs of patients with SCARs due to PB, PHT and CBZ.

## Materials and Methods

### Ethics statement

The study was reviewed and approved by the ethic local committee. All participating subjects were informed about the contents and the aims of the study and gave their written consent. After written informed consent signed, genomic DNA was extracted from peripheral blood, according to the established protocols.

### Subjects

All the 504 subjects were recruited from the University of São Paulo Faculty of Medicine Clinics Hospital, being 55 cases, 83 tolerants, and 366 healthy control subjects. Case subjects were
patients with a diagnosis of hypersensitivity reactions due to PB, PHT and CBZ within 1–8 week after the drug intake, including SJS, TEN and DRESS. We excluded subjects from the same family or with history of consanguinity. All the subjects were from São Paulo, Brazil, and the ethnicity of each participant was inferred from self-reported data about their four grand-parents ancestry utilizing Brazilian Institute of Geography and Statistics—IBGE—classification (White, Black, Brown, Indian and Yellow). Accordingly, we classified each individual by his/ her ancestry in African-American (at least one grandparent identified as Black), Caucasian or others (non-African descendant mixed population).

## Case subjects

The DRESS group consisted of 32 patients and SJS group of 12, defined by the RegiSCAR [22–25]. We recruited 11 cases of maculopapular exanthema (MPE), defined as rash without systemic symptoms, which required the discontinuation of CBZ within 3 months after the initiation of drug therapy [23].

## Control subjects

Tolerants group consisted of 83 subjects taking ARA for at least 6 months with no clinical or biochemical signs of hypersensitivity. All the subjects were from São Paulo, Brazil.

We included samples from 366 of healthy volunteers. We did not screen either set of population controls for carbamazepine-related adverse drug reactions.

## Drug Patch Test

The PTs were performed 6 weeks to 6 months after complete healing of the adverse drug reaction, and at least 1 month after discontinuation of systemic corticosteroids. PTs were applied on the upper back, using Finn Chambers with culprit ARA. The concentrations, readings and interpretation were performed according to the European Network of Drug Allergy recommendations [26,27]. The test was also performed in 10 control subjects, all resulted negative. The medications or conditions that could affect the interpretation of the test followed the guidelines of the International Contact Dermatitis Group.

## Genotyping

DNA was extracted from peripheral blood by salting-out [28]. SNPs characterizing CYP2C19*2, CYP2C19*3, CYP2C19*17, CYP2C9*2, CYP2C9*3 and CYP3A5*3 alleles (Table 1) were genotyped using TaqMan allele-specific polymerase chain reactions (Life Technologies). Amplification reactions were as follow: TaqMan PCR Mastermix 1x/µL, TaqMan SNP genotyping assay 1x/µL, genomic DNA 1ng/µL, ultrapure water to complete 7µL volume. rs12248560 required 500mM of Betaine (Sigma Aldrich) in the final reaction.

Allele discrimination was evaluated in a Line Gene 9600 (BIOER Technology CO.) comparing amplification curves and fluorescence levels before and after amplification (45 cycles of 15 seconds at 95°C and 1min at 60°C).

The candidate CYP genes were chosen according to the described association to the aromatic anticonvulsant-reactions [14,15,18,19,29] in different populations and the importance in the metabolism of these drugs (Table 1).

According to the alleles presented, each subject was classified with normal (EM), decrease (IM/PM) or increased (UM) function for each enzyme (Table 2).
<table><caption>Table 1. Candidate CYP genes, alleles, polymorphisms identification and enzymatic activity.</caption><thead><tr><th colspan="3">CYP2C19</th></tr><tr><th>Allele</th><th>Ref SNP ID</th><th>Enzymatic activity</th></tr></thead><tbody><tr><td>*2</td><td>rs4244285</td><td>Null</td></tr><tr><td>*3</td><td>rs4986893</td><td>Null</td></tr><tr><td>*17</td><td>rs12248560</td><td>Increased</td></tr><tr><td>*17</td><td>rs11188072</td><td>Increased</td></tr><tr><th colspan="3">CYP2C9</th></tr><tr><td>*2</td><td>rs1057910</td><td>Decreased</td></tr><tr><td>*3</td><td>rs1799853</td><td>Decreased</td></tr><tr><th colspan="3">CYP3A5</th></tr><tr><td>*3</td><td>rs776746</td><td>Decreased</td></tr></tbody></table>

doi:10.1371/journal.pone.0136141.t001

## Statistical Analysis

The data is available as a supporting information table (S1 Table). Descriptive statistical methods (mean, median, and frequency) were applied to analyze the demographic data. Analyses of the difference in frequencies across groups (sex and clinical manifestation) were performed with the chi-squared test or analysis of variance. CYP phenotypes were compared between diagnostic groups (Healthy, Tolerant and Cases) with Pearson's Chi-squared test. Significance was set at p<0.05. The CYP phenotypes were also entered in a Classification and Regression Tree (CART) model to predict tolerant cases. Variable selection was based on Gini impurity index and pruning was conducted based on 10-fold cross-validation to avoid over fitting. All analyses were conducted with the aid of the R software version 3.2.0.

## Results

### Clinical characteristics of aromatic anticonvulsant—induced hypersensitivity reactions patients and control subjects

There were no significant differences in the demographic data between groups as shown in Table 3.

<table><caption>Table 2. Candidate CYP genes, corresponding predicted phenotype and enzymatic activity according to the allele combination (EM = normal metabolizers, IM and PM = decrease function, UM and UMH = increase function).</caption><thead><tr><th>Gene</th><th>Predicted Phenotype</th><th>Enzymatic activity of allele combination</th></tr></thead><tbody><tr><td rowspan="4">CYP2C19</td><td>EM</td><td>Normal+Normal; Increased+Decreased</td></tr><tr><td>IM</td><td>Normal+Decreased</td></tr><tr><td>UM</td><td>Increased+Increased</td></tr><tr><td>UM heterozygote</td><td>Increased+Normal</td></tr><tr><td rowspan="4">CYP2C9</td><td>PM</td><td>Decreased+Decreased</td></tr><tr><td>EM</td><td>Normal+Normal</td></tr><tr><td>IM</td><td>Normal+Decreased</td></tr><tr><td>PM</td><td>Decreased+Decreased</td></tr><tr><td rowspan="3">CYP3A5</td><td>EM</td><td>Normal+Normal</td></tr><tr><td>IM</td><td>Normal+Decreased</td></tr><tr><td>PM</td><td>Decreased+Decreased</td></tr></tbody></table>

doi:10.1371/journal.pone.0136141.t002
<table><caption>Table 3. Demographic data of cases and controls subjects.</caption><thead><tr><th>Variables<br>Categorical Variables</th><th>Cases (N = 55)<br>N (%)</th><th>Tolerants (N = 85)<br>N (%)</th><th>Healthy Controls (N = 366)<br>N (%)</th><th>Total (506)<br>N (%)</th></tr></thead><tbody><tr><td><strong>Gender</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Female</td><td>42 (76)</td><td>47 (55)</td><td>293 (80)</td><td>382 (75)</td></tr><tr><td>Male</td><td>13 (24)</td><td>38 (45)</td><td>73 (20)</td><td>124 (25)</td></tr><tr><td><strong>Age</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Young adults (18–40 years)</td><td>23 (42)</td><td>26 (30)</td><td>147 (40)</td><td>196 (39)</td></tr><tr><td>Adults (40–60 years)</td><td>25 (46)</td><td>44 (52)</td><td>179 (49)</td><td>248 (49)</td></tr><tr><td>Elderly(&gt; 61 years)</td><td>7 (12)</td><td>15 (18)</td><td>40 (11)</td><td>62 (12)</td></tr><tr><td><strong>Ethnicity</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Caucasian</td><td>14 (26)</td><td>20 (24)</td><td>161 (44)</td><td>195 (38)</td></tr><tr><td>Afro-American</td><td>29 (53)</td><td>41 (48)</td><td>154 (42)</td><td>224 (45)</td></tr><tr><td>Others</td><td>12 (21)</td><td>24 (28)</td><td>51 (14)</td><td>87 (17)</td></tr></tbody></table>

doi:10.1371/journal.pone.0136141.t003

The severities of cutaneous manifestation of the reactions varied from MPE to exfoliative dermatitis. As shown in Table 4, 32 patients were validated as DRESS, 12 as SJS and 11 as MPE. CBZ was involved in 82% of the reactions, both mild and SCARs. The late onset of the reaction was important in cases of DRESS, but happened in low proportion in SJS. No MPE happened after 4 weeks of drug intake. SCARs showed to be systemic reactions with internal organ involvement requiring long-term hospitalization, mainly in DRESS cases but also in SJS (Table 4). Most of the cases, both mild and severe, were treated with systemic corticosteroids and anti-histamines. Intravenous immunoglobulin (IgIV) was exclusively indicated to SJS cases. We found no association between ethnicity and clinical status (p = 0.231).

Of all 46 PTs, 29 (63%) were positive and 63% of positive results were in DRESS. The PTs was not performed in 9 cases due to the difficulties on the clinical follow-up or if the patient had contraindications to the procedure. Apparently grade of positivity of PTs did not correlated with the intensity of reaction nor the internal organ involvement, however, positive reactions occurred in the more severe cases. There was no association between positivity of the tests and skin color or reported ethnicity. None of the patients presented adverse reactions during the in vivo procedure.

## CYP3A5, CYP2C9 and CYP2C19 genotyping and aromatic anticonvulsant—induced hypersensitivity reactions

On the basis of our hypothesis that the polymorphisms of CYP2C19, CYP2C9 and CYP3A5 are associated with aromatic anticonvulsant—induced hypersensitivity reactions, we genotyped the CYP2C19*2, CYP2C19*3, CYP2C19*17, CYP2C9*2, CYP2C9*3 and CYP3A5*3 alleles in all subjects included in the study (Tables 1 and 2). The allele frequencies and its distribution among the enzymatic phenotype and the studied groups are shown in Table 5.

We showed a strong association between the absence of CYP3A5*3 and tolerant subjects when compared to cases (p = 0.0002, OR = 5.28 [CI95% 2.09–14.84]). However, we observed homogeneous distribution of variants of CYP2C19 and 2C9 in our case and control groups (Table 5). The classification and regression tree showed that non-decreased function of CYP3A5 (absence of *3 allele; EM), was the most important factor to classify tolerant individuals (Fig 1). On the basis of this result, the absence of CYP3A5*3 has a putative sensitivity of 81.1% as a protective factor to aromatic antiepileptics-hypersensitivity reactions, whereas the presence of IM or PM phenotypes is split in 54.32% cases and 45.67% tolerants.
<table><caption>Table 4. Clinical patterns of aromatic anticonvulsant-induced reactions and demographic characteristics (DRESS = Drug Reaction with eosinophilia and systemic symptoms, SJS = Stevens-Johnson Syndrome, MPE = Maculopapular exanthema, IgIV = Intravenous immunoglobulin).</caption><thead><tr><th></th><th>DRESS (N:32)</th><th>SJS (N:12)</th><th>MPE (N:11)</th><th>TOTAL (N:55)</th></tr></thead><tbody><tr><td>Mean age</td><td>49.7</td><td>48.3</td><td>47.7</td><td>43.6</td></tr><tr><td>Sex (M/F)</td><td>26/6</td><td>7/5</td><td>9/2</td><td>42/13</td></tr><tr><td><strong>Onset of the reaction</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>&lt; 4 weeks (%)</td><td>13 (40)</td><td>11 (92)</td><td>11 (100)</td><td>35 (64)</td></tr><tr><td>&gt; 4 weeks (%)</td><td>19 (60)</td><td>1 (8)</td><td>0 (0)</td><td>20 (36)</td></tr><tr><td><strong>Drug involved</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Phenobarbital</td><td>4 (13)</td><td>2 (17)</td><td>1 (9)</td><td>7 (13)</td></tr><tr><td>Phenytoin</td><td>2 (6)</td><td>1 (8)</td><td>0 (0)</td><td>3 (5)</td></tr><tr><td>Carbamazepine</td><td>26 (81)</td><td>9 (75)</td><td>10 (91)</td><td>45 (82)</td></tr><tr><td><strong>Clinical manifestations</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Cutaneous (%)</td><td>32 (100)</td><td>12 (100)</td><td>11 (100)</td><td>55 (100)</td></tr><tr><td>Gastrointestinal (%)</td><td>32 (100)</td><td>6 (50)</td><td>0 (0)</td><td>38 (69)</td></tr><tr><td>Respiratory (%)</td><td>5 (15)</td><td>4 (34)</td><td>0 (0)</td><td>9 (16)</td></tr><tr><td>Fever (%)</td><td>32 (100)</td><td>11 (91)</td><td>2 (18)</td><td>45 (82)</td></tr><tr><td>Mean time of hospitalization (days)</td><td>20</td><td>12</td><td>0</td><td>24</td></tr><tr><td><strong>Treatment</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Anti-histamines (%)</td><td>32 (100)</td><td>8 (67)</td><td>11 (100)</td><td>51 (93)</td></tr><tr><td>Systemic corticosteroids (%)</td><td>32 (100)</td><td>12 (100)</td><td>8 (72)</td><td>52 (95)</td></tr><tr><td>IgIV (%)</td><td>0 (0)</td><td>5 (41)</td><td>0 (0)</td><td>5 (9)</td></tr><tr><td><strong>Drug patch test (N: 46)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Positive</td><td>18 (72)</td><td>7 (58)</td><td>4 (45)</td><td>29 (63)</td></tr><tr><td>Negative</td><td>7 (28)</td><td>5 (42)</td><td>5 (55)</td><td>17 (37)</td></tr><tr><td>Not Done</td><td>7</td><td>0</td><td>2</td><td>9 (16)</td></tr></tbody></table>

doi:10.1371/journal.pone.0136141.t004

Given the fact that CYP2C19, CYP2C9 and CYP3A5 enzymes are expressed in the skin, we analyzed the association between genotype and the results of PT. Even without statistically significant associations, we observed a predominance of low function phenotype of CYP3A5, normal function of CYP2C9, and increased function of CYP2C19 in positive results.

<table><caption>Table 5. Phenotypic frequencies for cytochrome P450 (CYP) 2C9, 2C19 and 3A5 in healthy controls, cases and tolerants (EM = normal metabolizers, IM and PM = decrease function, UM and UMH = increase function) as predicted by SNP genotype.</caption><thead><tr><th rowspan="2">Gene (SNP)</th><th rowspan="2">Enzymatic Phenotype Group</th><th colspan="3">Phenotypic Frequency (%)</th><th rowspan="2">P Value</th></tr><tr><th>Healthy controls (%)</th><th>Cases (%)</th><th>Tolerants (%)</th></tr></thead><tbody><tr><td rowspan="3">CYP2C9 (rs1057910, rs1799853)</td><td>EM</td><td>251 (71.1)</td><td>43 (81.1)</td><td>54 (65.9)</td><td rowspan="3">0.083</td></tr><tr><td>IM</td><td>97 (27.5)</td><td>8 (15.1)</td><td>25 (30.5)</td></tr><tr><td>PM</td><td>5 (1.4)</td><td>2 (3.8)</td><td>3 (3.7)</td></tr><tr><td rowspan="3">CYP2C19 (rs4244285, rs4986893, rs12248560, rs11188072)</td><td>EM</td><td>163 (48.4)</td><td>24 (47.1)</td><td>41 (51.3)</td><td rowspan="3">0.945</td></tr><tr><td>IM/PM</td><td>79 (23.4)</td><td>12 (23.5)</td><td>20 (25)</td></tr><tr><td>UM/UMH</td><td>95 (28.2)</td><td>15 (29.4)</td><td>19 (23.8)</td></tr><tr><td rowspan="3">CYP3A5 (rs776746)</td><td>EM</td><td>37 (10.6)</td><td>8 (14.8)</td><td>38 (46.3)</td><td rowspan="3">0.0005<br>OR = 4.9</td></tr><tr><td>IM</td><td>122 (34.9)</td><td>23 (42.6)</td><td>19 (23.2)</td></tr><tr><td>PM</td><td>191 (54.6)</td><td>23 (42.6)</td><td>25 (30.5)</td></tr></tbody></table>

doi:10.1371/journal.pone.0136141.t005
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Classification and regression tree: CYP3A5 phenotype discriminates between cases and tolerants (EM = normal metabolizers, IM and PM = decrease function, UM and UMH = increase function).</figcaption></figure>

doi:10.1371/journal.pone.0136141.g001

Discussion

This study is the first and most broad case-control study ever conducted towards the associa-
tion of CYP alleles and aromatic antiepileptics-induced hypersensitivity reactions in Brazilian
subjects. Our data add to the growing evidence of different CYP alleles role in immune-medi-
ated ADRs, such as drug-induced hypersensitivity [5,6]. Genetic variation in genes involved in
their metabolism may result in the increase of arene oxide, which potentially initiates immune
responses [8]. On the other hand, the immune response to CBZ in AHR has been associated
with the drug itself and not its metabolites [6].

Genetic polymorphisms of CYP2C19 are highly present in Asian population and the most
common variant resulting in poor metabolism are CYP2C19*2 and CYP19*3 [16,30]. The
allele frequency of CYP2C9*2 and CYP2C9*3 vary among different populations, but seems to
be more frequent in Caucasians. The CYP3A5 gene is highly polymorphic and the allele
CYP3A5*3 is the most prevalent related to the loss of activity. The prevalence of this allele
appears to be higher in Caucasians when compared to African Americans. The CYP3A5*3/*3
genotype was associated with increased half-life of CBZ in African-Americans, but no signifi-
cant association was observed in Caucasians [16,29,31].

We found that the absence of the CYP3A5*3 is an important protective factor of the full
spectrum of aromatic antiepileptics-induced hypersensitivity reactions in Brazilian subjects.
The CYP3A5*3 has been shown to be related to the decreased function of this enzyme. The
association of drug tolerance with normal activity of CYP3A5 supports previous works indicat-
ing that immune response to CBZ in AHR is due to the drug itself and not its metabolytes
[6,32,33]. Whereas we demonstrated that the absence of this allele is related with increased tolerability of these drugs, the similar proportion of CYP3A5*3 in both cases and controls with IM/PM phenotypes underlined that it cannot be considered an isolated risk factor for the development of antiepileptic hypersensitivity reactions.

However, the presence of the CYP3A5*3 allele is neither necessary nor sufficient for the development of these reactions. Although CYP3A5*3 is the most prevalent CYP3A5 allele with reduced function, it is not the only one [19]. Thus, it remains to be evaluated if cases in our sample that were considered with normal CYP3A5 function (EM) are actually carriers of rarer non-functional alleles. Most studies published in the last years, which have been searching for genetic risk factors for the development of drug-induced hypersensitivity reactions, showed the predictive value of HLA alleles [34,35]. The current study now provides evidence suggesting that the absence of the CYP3A5*3 is a protective factor to aromatic antiepileptics-induced hypersensitivity reactions in Brazilian subjects. The worldwide distribution of CYP3A5*3 provides a remarkable example of population diversity, with allele frequencies ranging from 0.14 among sub-Saharan Africans to 0.95 in European populations [18]. The Brazilian population, in excess of 195 million people, has major ancestral roots in Europe, Africa and America, nevertheless, we didn't find significant association between ethnicity and clinical status. It has been recently reported that the genetic background of the Brazilian population is less heterogeneous than previously believed, with a high percentage (60-77%) of caucasian background [36]. This study is the first to show the association between CYP3C5*3 and aromatic anticonvulsant-induced reactions in subjects from São Paulo, Brazil. CYP3C5*3 may have a potential role in future screening before prescribing aromatic anticonvulsant, especially CBZ. Since we do not know the responsiveness of the healthy population to anticonvulsants we calculated negative (NPV) and positive predictive value (PPV) for the absence of CYP3A5*3 to predict tolerability considering cases and tolerants subjects. In this scenario it shows a NPV of 51% and PPV of 82%. According to our data, testing would prevent 4 cases of adverse reaction in every 5 prescriptions of anticonvulsants (NNT = 4 [2.1-5.5]).

We did not find association between CYP2C19 and CYP2C9 variants and cases in our study. The CYP2C9 and CYP2C19 variants have been described as risk factors for the development of SCARs in Thai and Japanese populations [16,17]. In addition, CYP2C9 has been shown to be associated with PHT metabolism and elimination, but not with PB or CBZ and CYP2C19 may be the major metabolic pathway of PB [15]. Most cases of this study were SCARs due to CBZ in a Brazilian cohort of mixed ethnic background, which might explain the discrepancies with the literature. Patch Tests are the only in vivo diagnostic method available to SCARs up to now, with positive reactions in 60-100% of the cases studied [20,21]. Nevertheless, the reactivity depends on the type of the cutaneous adverse drug reactions. The positivity of PTs in our study varied according to the clinical pattern, being more positive (72%) in cases of DRESS. The CYP2C9, CY2C19 and CYP3A5 enzymes have been described as being expressed in human skin [37]. Here, we hypothesized that they can influence in the positivity of the ARA patch test, however we found no association between functional CYP3A5 genotype in individuals with positive PTs. CYP activation is a possible mechanism for allergic contact dermatitis (ACT), a complex syndrome representing immunological responses to cutaneous and/or systemic exposure to protein-reactive chemicals. Recently, many drug-induced hypersensitivity reactions such as DRESS and SJS have been classified in a similar pathophysiologic scenario as the ACT.

We are aware that the current findings may not be taken as isolated genetic factors to determine lower risk on developing AHR. Considering the ongoing knowledge on the specific HLA alleles as risk factors for SCARs, we would expect that combined analysis of CYP and HLA data could help identify subjects with higher or lower risk of developing AHR. HLA-A*3101
has been associated with AHR in European population while HLA-B*1502 in Chinese. Given the relative rarity of those two alleles in the Brazilian population [38] and the ethnic specificity they present we do not believe that they would substantially modify the present analysis. However, further studies are necessary in order to rule out this possibility.

Some limitations of the present study are similar to the majority of the association studies linking genetic factors and antiepileptics-induced cutaneous adverse reactions, such as the limited number of cases, the population included were exclusively of subjects living in São Paulo and the ancestry was based on self-identification. Due to the self-identified admixed ancestry in our cohort, we believe that the current data may require replication by future multicentric multinational collaborations with accurate phenotypic characterization studies in populations of other ancestry subjects.

We therefore suggest that consideration be given to the absence of CY3A5*3 as a genetic preventive marker of aromatic anticonvulsants hypersensitivity reactions.

## Supporting Information

S1 Table. Data utilized in this study in a CSV format file. Ethnicity (AA—afro american, C —caucasian, O—others), Clinical_Pattern (0—Healthy, 1—Case, 2—Tolerant), Clinical_Pattern_Severity (0—Healthy, 1—Severe Reaction, 2—Tolerant, 3—Mild reaction), Typekip Of Reaction (0—Healthy, 2—Tolerant, 3—Mild reaction, 4—DRESS, 5—Steven Johnson), Contact (P—Positive, N—Negative, ND—Not Done).
(CSV)

## Acknowledgments

This study had a financial support provided by FAPESP (Fundação de Amparo à Pesquisa do Estado de São Paulo), under the grants number 2011/22748-1 and 2013/01352-8.

## Author Contributions

Conceived and designed the experiments: LKT JK DSK. Performed the experiments: LKT CY DSK. Analyzed the data: LKT BS DSK. Contributed reagents/materials/analysis tools: LKT CY HR JK DSK WFG LLT. Wrote the paper: LKT DSK BS LLT CY HR WFG JK. Designed the study, analyzed and interpreted the data, and wrote the manuscript: LKT. Designed the study, analyzed the data and reviewed the paper: JK. Designed the study, proceeded the experimental analysis and wrote the paper: DSK. Was responsible for the statistical analysis and reviewed the paper: BS. Helped in the interpretation of the data and with the revision of the manuscript: LLT CY HR WFG.

## References

1. Marson AG, Al-Kharusi AM, Alwaidh M, Appleton R, Baker GA, Chadwick DW, et al. The SANAD study of effectiveness of valproate, lamotrigine, or topiramate for generalised and unclassifiable epilepsy: an unblinded randomised controlled trial. Lancet. 2007; 369: 1016–1026. doi: 10.1016/S0140-6736(07) 60461-9 PMID: 17382828

2. Wiffen PJ, Derry S, Moore RA, McQuay HJ. Carbamazepine for acute and chronic pain in adults. Cochrane Database Syst Rev. 2011; Issue 1. Available: http://www.mrw.interscience.wiley.com/ cochrane/clsysrev/articles/CD005451/frame.html

3. Roujeau JC, Stern RS. Severe adverse cutaneous reactions to drugs. N Engl J Med. 1994; 331: 1272– 85. doi: 10.1056/NEJM199411103311906 PMID: 7794310

4. Roujeau JC, Kelly JP, Naldi L, Rzany B, Stern RS, Anderson T, et al. Medication use and the risk of Stevens-Johnson syndrome or toxic epidermal necrolysis. N Engl J Med. 1995; 333: 1600–7. doi: 10.1056/ NEJM199512143332404 PMID: 7477195
5. Posadas SJ, Pichler WJ. Delayed drug hypersensitivity reactions—new concepts. Clin Exp Allergy. 2007; 37: 989–999. doi: 10.1111/j.1365-2222.2007.02742.x PMID: 17581192

6. Naisbitt DJ, Britschgi M, Wong G, Farrell J, Depta JPH, Chadwick DW, et al. Hypersensitivity reactions to carbamazepine: characterization of the specificity, phenotype, and cytokine profile of drug-specific T cell clones. Mol Pharmacol. 2003; 63: 732–41. Available: http://www.ncbi.nlm.nih.gov/pubmed/12606784 PMID: 12606784

7. Edwards SG, Hubbard V, Aylett S, Wren D. Concordance of primary generalised epilepsy and carbamazepine hypersensitivity in monozygotic twins. Postgraduate medical journal. 1999. pp. 680–681. doi: 10.1136/pgmj.75.889.680 PMID: 10621883

8. Spielberg S, Gordon G, Blake D, Mellits E, Bross D. Anticonvulsant toxicity in vitro: possible role of arene oxides. J Pharmacol Exp Ther. 1981; 217: 386–389. Available: http://jpet.aspetjournals.org/content/217/2/386.short PMID: 7229980

9. Green VJ, Pirmohamed M, Kitteringham NR, Gaedigk A, Grant DM, Boxer M, et al. Genetic analysis of microsomal epoxide hydrolase in patients with carbamazepine hypersensitivity. Biochem Pharmacol. 1995; 50: 1353–1359. doi: 10.1016/0006-2952(95)02009-8 PMID: 7503783

10. Madden S, Maggs JL, Park BK. Bioactivation of carbamazepine in the rat in vivo: Evidence for the formation of reactive arene oxide(s). Drug Metab Dispos. 1996; 24: 469–479. PMID: 8801063

11. Lertratanangkoon K, Horning MG. Metabolism of carbamazepine. Drug Metab Dispos. 1982; 10: 1–10. PMID: 6124375

12. Gaedigk A, Spielberg SP, Grant DM. Characterization of the microsomal epoxide hydrolase gene in patients with anticonvulsant adverse drug reactions. Pharmacogenetics. 1994; 4: 142–153. PMID: 7920694

13. Baron JM, Wiederholt T, Heise R, Merk HF, Bickers DR. Expression and function of cytochrome p450-dependent enzymes in human skin cells. Curr Med Chem. 2008; 15: 2258–2264. doi: 10.2174/092986708785747535 PMID: 18781947

14. Thorn CF, Leckband SG, Kelsoe J, Leeder JS, Müller DJ, Klein TE, et al. PharmGKB summary: carbamazepine pathway. Pharmacogenet Genomics. 2011; 21: 906–10. PMID: 21738081

15. Thorn CF, Whirl-Carrillo M, Leeder JS, Klein TE, Altman RB. PharmGKB summary: phenytoin pathway. Pharmacogenet Genomics. 2012; 22: 466–70. PMID: 22569204

16. Ozawa S, Soyama A, Saeki M, Fukushima-Uesaka H, Itoda M, Koyano S, et al. Ethnic differences in genetic polymorphisms of CYP2D6, CYP2C19, CYP3As and MDR1/ABCB1. Drug Metab Pharmacokinetic. 2004; 19: 83–95. doi: 10.2133/dmpk.19.83 PMID: 15499174

17. Goldstein JA, Ishizaki T, Chiba K, de Morais SM, Bell D, Krahn PM, et al. Frequencies of the defective CYP2C19 alleles responsible for the mephenytoin poor metabolizer phenotype in various Oriental, Caucasian, Saudi Arabian and American black populations. Pharmacogenetics. 1997; 7: 59–64. PMID: 9110363

18. International HapMap Project [Internet]. 2014. Available: http://hapmap.ncbi.nlm.nih.gov/

19. Hustert E, Haberl M, Burk O, Wolbold R, He YQ, Klein K, et al. The genetic determinants of the CYP3A5 polymorphism. Pharmacogenetics. 2001; 11: 773–779. PMID: 11740341

20. Barbaud A, Collet E, Milpied B, Assier H, Staumont D, Avenel-Audran M, et al. A multicentre study to determine the value and safety of drug patch tests for the three main classes of severe cutaneous adverse drug reactions. Br J Dermatol. 2013; 168: 555–562. doi: 10.1111/bjd.12125 PMID: 23136927

21. Santiago E, Gonçalo M, Vieira R, Coelho S, Figueiredo A. Epicutaneous patch testing in drug hypersensitivity syndrome (DRESS). Contact Dermatitis. 2010; 62: 47–53. doi: 10.1111/j.1600-0536.2009.01659.x PMID: 20136879

22. Mockenhaupt M. The RegiSCAR Project. In: Inclusion Criteria [Internet]. 2014. Available: http://www.regiscar.org/pdf/inclusioncriteriaforscar.pdf

23. Pirmohamed M, Friedmann PS, Molokhia M, Loke YK, Smith C, Phillips E, et al. Phenotype standardization for immune-mediated drug-induced skin injury. Clin Pharmacol Ther. 2011; 89: 896–901. doi: 10.1038/clpt.2011.79 PMID: 21562486

24. Lonjou C, Thomas L, Borot N, Ledger N, de Toma C, LeLouet H, et al. A marker for Stevens-Johnson syndrome...: ethnicity matters. Pharmacogenomics J. 2006; 6: 265–8. doi: 10.1038/sj.tpj.6500356 PMID: 16415921

25. Auquier-Dunant A, Mockenhaupt M, Naldi L, Correia O, Schröder W, Roujeau J-C. Correlations between clinical patterns and causes of erythema multiforme majus, Stevens-Johnson syndrome, and toxic epidermal necrolysis: results of an international prospective study. Arch Dermatol. 2002; 138: 1019–24. Available: http://www.ncbi.nlm.nih.gov/pubmed/12164739 PMID: 12164739

26. Demoly P, Adkinson NF, Brockow K, Castells M, Chiriac AM, Greenberger PA, et al. International Consensus on drug allergy. Allergy. 2014; 69: 420–37. doi: 10.1111/all.12350 PMID: 24697291
27. Brockow K, Garvey LH, Aberer W, Atanaskovic-Markovic M, Barbaud A, Bilo MB, et al. Skin test concentrations for systemically administered drugs—An ENDA/EAACI Drug Allergy Interest Group position paper. Allergy Eur J Allergy Clin Immunol. 2013; 68: 702–712. doi: 10.1111/all.12142

28. Laitinen J, Samarut J, Holtta E. A nontoxic and versatile protein salting-out method for isolation of DNA. Biotechniques. 1994; 17: 316, 318, 320–2. Available: http://www.ncbi.nlm.nih.gov/pubmed/7980935 PMID: 7980935

29. Sim SC, Ingelman-Sundberg M. The Human Cytochrome P450 (CYP) Allele Nomenclature website: a peer-reviewed database of CYP variants and their associated effects. Hum Genomics. 2010; 4: 278–81. Available: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3525213&tool = pmcentrez&rendertype = abstract PMID: 20511141

30. Xie HG, Prasad HC, Kim RB, Stein CM. CYP2C9 allelic variants: Ethnic distribution and functional significance. Advanced Drug Delivery Reviews. 2002. pp. 1257–1270. doi: 10.1016/S0169-409X(02) 00076-5 PMID: 12406644

31. Suarez-Kurtz G, Vargens DD, Santoro AB, Hutz MH, De Moraes ME, Pena SDJ, et al. Global pharmacogenomics: Distribution of CYP3A5 polymorphisms and phenotypes in the Brazilian population. PLoS One. 2014; 9. doi: 10.1371/journal.pone.0083472

32. Wu Y, Sanderson JP, Farrell J, Drummond NS, Hanson A, Bowkett E, et al. Activation of T cells by carbamazepine and carbamazepine metabolites. J Allergy Clin Immunol. 2006; 118: 233–41. doi: 10. 1016/j.jaci.2006.03.005 PMID: 16815161

33. Lee A-Y, Choi J, Chey W-Y. Patch testing with carbamazepine and its main metabolite carbamazepine epoxide in cutaneous adverse drug reactions to carbamazepine. Contact Dermatitis. 2003; 48: 137–9. Available: http://www.ncbi.nlm.nih.gov/pubmed/12755726 PMID: 12755726

34. McCormack M, Alfirevic A, Bourgeois S, Farrell JJ, Kasperaviči tė D, Carrington M, et al. HLA-A*3101 and carbamazepine-induced hypersensitivity reactions in Europeans. N Engl J Med. 2011; 364: 1134–1143. doi: 10.1056/NEJMoa1013297 PMID: 21428769

35. Hung S-I, Chung W-H, Jee S-H, Chen W-C, Chang Y-T, Lee W-R, et al. Genetic susceptibility to carbamazepine-induced cutaneous adverse drug reactions. Pharmacogenet Genomics. 2006; 16: 297–306. PMID: 16538176

36. Pena SDJ, Di Pietro G, Fuchshuber-Moraes M, Genro JP, Hutz MH, Kehdy F de SG, et al. The genomic ancestry of individuals from different geographical regions of Brazil is more uniform than expected. PLoS One. 2011; 6: e17063. doi: 10.1371/journal.pone.0017063 PMID: 21359226

37. Nishimura M, Yaguti H, Yoshitsugu H, Naito S, Satoh T. Tissue distribution of mRNA expression of human cytochrome P450 isoforms assessed by high-sensitivity real-time reverse transcription PCR. Yakugaku Zasshi. 2003; 123: 369–375. PMID: 12772594

38. González-Galarza FF, Takeshita LY, Santos EJ, Kempson F, Maia MH, da Silva AL, et al. Allele frequency net 2015 update: new features for HLA epitopes, KIR and disease and HLA adverse drug reaction associations. Nucleic Acids Res. 2015 Jan; 43 (Database issue):D784–8. doi: 10.1093/nar/ gku1166 PMID: 25414323