==== Front Allergy Asthma Immunol ResAllergy Asthma Immunol ResAAIRAllergy, Asthma & Immunology Research2092-73552092-7363The Korean Academy of Asthma, Allergy and Clinical Immunology; The Korean Academy of Pediatric Allergy and Respiratory Disease 10.4168/aair.2018.10.5.555Original ArticleThe Fas Signaling Pathway Is a Common Genetic Risk Factor for Severe Cutaneous Drug Adverse Reactions Across Diverse Drugs https://orcid.org/0000-0002-6970-3228Park Heung-Woo 12Kim Sang-Heon 3Chang Yoon-Seok 4Kim Sang-Hoon 5Jee Young-Koo 6Lee Ai-Young 7Jang In-Jin 8https://orcid.org/0000-0003-2614-0303Park Hae-Sim 9Min Kyung-Up 121 Department of Internal Medicine, Seoul National University College of Medicine, Seoul, Korea.2 Institute of Allergy and Clinical Immunology, Seoul National University Medical Research Center, Seoul, Korea.3 Department of Internal Medicine, Hanyang University College of Medicine, Seoul, Korea.4 Department of Internal Medicine, Seoul National University Bundang Hospital, Seongnam, Korea.5 Department of Internal Medicine, Eulji University School of Medicine, Daejeon, Korea.6 Department of Internal Medicine, Dankook University College of Medicine, Cheonan, Korea.7 Department of Dermatology, Dongguk University Ilsan Hospital, Goyang, Korea.8 Department of Clinical Pharmacology and Therapeutics, Seoul National University College of Medicine, Seoul, Korea.9 Department of Allergy and Clinical Immunology, Ajou University School of Medicine, Suwon, Korea.Correspondence to Heung-Woo Park, MD, PhD. Department of Internal Medicine, Seoul National University College of Medicine, 103 Daehak-ro, Jongno-gu, Seoul 03080, Korea. Tel: +82-2-2072-0699, Fax: +82-2-742-3291, guinea71@snu.ac.kr9 2018 09 7 2018 10 5 555 561 06 2 2018 25 5 2018 01 6 2018 Copyright © 2018 The Korean Academy of Asthma, Allergy and Clinical Immunology • The Korean Academy of Pediatric Allergy and Respiratory Disease2018The Korean Academy of Asthma, Allergy and Clinical Immunology • The Korean Academy of Pediatric Allergy and Respiratory DiseaseThis is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (https://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.Purpose Human leukocyte antigen (HLA) has been recognized as the most important genetic risk factor for severe cutaneous adverse drug reactions (SCARs) caused by certain drugs. However, cumulated observations suggest the presence of genetic risk factors for SCARs other than drug-specific HLA. We aimed to identify a common genetic risk factor of SCARs across multiple drugs. Methods We performed 2 independent genome-wide association studies (GWASs). A total of 68 and 38 subjects with a diagnosis of SCAR were enrolled in each GWAS. Their allele frequencies were compared to those of healthy subjects in Korea. Results No single nucleotide polymorphism (SNP) with genome-wide significance was found in either GWAS. We next selected and annotated the 200 top-ranked SNPs from each GWAS. These 2 sets of annotated genes were then entered into the web interface of ConsensusPathDB for a pathway-level analysis. The Fas signaling pathway was significantly over-represented in each gene set from the 2 GWASs. Conclusions Our observations suggest that the Fas signaling pathway may be a common genetic risk factor for SCARs across multiple drugs. DrugStevens-Johnson syndrometoxic epidermal necrolysisgenome-wide association studyFas signaling pathway ==== Body INTRODUCTION Severe cutaneous adverse drug reactions (SCARs) characterized by an acute inflammatory reaction of the skin and mucous membranes are very rare, but sometimes fatal.1 SCARs include Stevens-Johnson syndrome (SJS), toxic epidermal necrolysis (TEN), and drug rash with eosinophilia and systemic symptoms (DRESS).1 Thus far, human leukocyte antigen (HLA) has been recognized as the most important genetic risk factor for SCARs caused by certain drugs.2 However, we reported that only 20% of carriers of the HLA-B*58:01 allele experienced an allopurinol-induced SCAR.3 Similarly, the risk of carbamazepine-induced hypersensitivity was 26% in carriers of HLA-A*31:01 in a European population.4 These observations suggest the presence of genetic risk factors for SCARs other than drug-specific HLA. The purpose of this study was to identify common genetic risk factors of SCARs in addition to drug-specific HLA. Cases of SCAR induced by various drugs were collected and a genome-wide association study (GWAS) was performed. As previously mentioned, drug-specific HLA is believed to be the strongest risk factor for SCAR. Therefore, it is possible that genetic variants with modest effects on SCAR are not captured in a GWAS. To overcome this problem, we used the gene set-based approach, which increases the likelihood of identifying the biological pathway involved by utilizing enrichment tools.5 MATERIALS AND METHODS We performed 2 independent GWASs in 2008 (the first) and 2012 (the second). All subjects with SCAR were enrolled by the Adverse Drug Reaction Research Group, Korea. Each study was approved by the Institutional Review Board of the Seoul National University Hospital (H-0506-150-011 and C-1111-100-387) and informed consent was obtained from all study participants. Study subjects and phenotype definitions SJS, SJS/TEN overlap and TEN were diagnosed according to the range of detached surface area (< 10%, 10%–30% and > 30%, respectively).67 DRESS was defined when a subject met at least five of the following criteria: development of a maculopapular rash more than 3 weeks after starting a limited number of drugs, prolonged clinical symptoms for 2 weeks after discontinuation of the causative drug, fever (> 38°C), liver abnormalities or other organ involvement, leukocyte abnormalities such as leukocytosis (>11 × 109 cells/L), atypical lymphocytosis (> 5%), or eosinophilia (>1.5 × 109 cells/L), lymphadenopathy, and human herpes virus 6 reactivation.8 All clinical evaluations were performed by allergy specialists with more than 10 years of experience in drug allergy consultation. A drug was classified as a cause of SCAR when it satisfied the grade of ‘certain’ or ‘probable’ on the World Health Organization-Uppsala Monitoring Center causality assessment system.9 Characteristics of subjects with SCAR are shown in Table 1. Control subjects were selected from a publicly available database to provide a 4:1 match with the case subjects. All control subjects were healthy Koreans and were genotyped on the same chip in each GWAS. Table 1 Characteristics of SCAR subjects Characteristic First GWAS (n = 68) Second GWAS (n = 38) Male 35 (51.5) 35 (51.5) Age (yr) 49.8 ± 16.8 52 ± 19.2 Phenotype SJS/Overlap*/TEN/DRESS 14 (20.6)/12 (17.6)/4 (5.9)/38 (55.9) 7 (18.4)/6 (15.8)/4 (10.5)/21 (55.3) Causative drug Allopurinol 13 (19.1) 5 (13.2) Carbamazepine 18 (26.5) 3 (7.9) Anti-epileptics except carbamazepine 7 (10.3)† 6 (15.8)¶ Acetaminophen 3 (4.4) 3 (7.9) NSAID 6 (8.8)‡ 12 (31.5)** Beta-lactam antibiotics 7 (10.3)§ 9 (23.7)†† Others 14 (20.6)∥ - Data are shown as number (%) or mean ± standard deviation. SCAR, severe cutaneous adverse drug reaction; GWAS, genome-wide association study; SJS, Stevens-Johnson syndrome; TEN, toxic epidermal necrolysis; DRESS, drug rash with eosinophilia and systemic symptom; NSAID, non-steroidal anti-inflammatory drug. *SJS/TEN overlap syndrome; †Including phenytoin (3), valproic acid (2), and lamotrigine (2); ‡Including ibuprofen (4) and mefenamic acid (2); §Including ceftezole (2), cefaclor (2), piperacillin (1), cefazolin (1), and cefpiramide (1); ∥Including anti-tuberculosis medication (6), oriental herbal medication (3), vancomycin (2), sulfasalazine (1), pantoprazole (1), and propylthiouracil (1); ¶Including phenytoin (5) and lamotrigine (1); **Including ibuprofen (6), loxoprofen (3), and mefenamic acid (3); ††Including amoxicillin (2), cefaclor (2), cefbuperazone (1), cefixime (1), cefoxitin (1), cephradine (1), and meropenem (1). Genotyping DNA of subjects included in the first GWAS was genotyped on the Affymetrix 500K Chip (Affymetrix, Santa Clara, CA, USA). DNA of subjects included in the second GWAS was genotyped on the Affymetrix Axiom® Genome-Wide East Asian Ancestry Chip. All single nucleotide polymorphisms (SNPs) that were included in the GWAS had a completion rate of >95%, a minor allele frequency of >0.05, and a Hardy-Weinberg equilibrium P value >0.0001. Statistical analysis Associations between SNPs and the occurrence of SCAR were measured according to a linear regression model, as implemented in PLINK,10 using an additive genetic model. The regression models were adjusted for age and sex. SNPs with P values < 10−8 in the GWAS were considered genome-wide significant. For over-representation analysis, we selected and annotated the 200 top-ranked SNPs from both GWASs using the open-source software package NCBI2R in the statistical and graphical environment of R software (http://www.r-project.org). Only genes that harbored the SNP of interest were selected; 62 genes were identified from the first GWAS and 69 were identified from the second GWAS. Gene lists are provided in Supplementary Tables S1 and S2. The list of genes was then entered into the web interface of ConsensusPathDB (http://cpdb.molgen.mpg.de) for pathway-level interpretation. ConsensusPathDB is a meta-database that integrates different types of functional interactions from heterogeneous interaction data resources.11 To control for multiple tests, 0.05 was set as the limit of Q values, which were calculated using the false discovery rate method. RESULTS The top-ranked SNP in the first GWAS was rs2327661 (P = 6.3 × 10−7) and that in the second GWAS was rs8180036 (P = 4.48 × 10−6). However, none of the SNP reached genome-wide significance. A list of the 200 top-ranked SNPs is provided in Supplementary Table S3 and the Manhattan plots each GWAS are presented in Supplementary Fig. S1. Over-representation analysis of the annotated gene sets by ConsensusPathDB identified 17 significant pathways in the first GWAS and 3 significant pathways in the second GWAS, see Tables 2 and 3. Interestingly, the Fas signaling pathway (the fourth ranked pathway in the first GWAS and the 1st ranked pathway in the second GWAS) was commonly found. Two genes (CASP10 and MAP3K7) in the gene set from the first GWAS and 3 genes (PTPN13, PRKDC and FADD) in the gene set from the second GWAS were significantly over-represented in the Fas signaling pathway, see Supplementary Table S4. Table 2 Pathway-based gene sets identified from the first GWAS P value Q value Pathway Source Overlapped genes* 0.00069918 0.02407312 Propanoate metabolism SMPDB ACACA; ABAT 0.00091873 0.02407312 IL-4 signaling pathway PID, BioCarta HLA-DRA; JAK1 0.00109038 0.02407312 Leishmaniasis - Homo sapiens (human) KEGG HLA-DRA; MAP3K7; JAK1 0.00144439 0.02407312 Fas signaling pathway (cd95) PID, BioCarta CASP10; MAP3K7 0.00244457 0.02575628 A tetrasaccharide linker sequence is required for GAG synthesis Reactome VCAN; GPC5 0.00263548 0.02575628 Alk in cardiac myocytes PID, BioCarta BMP7; MAP3K7 0.00283319 0.02575628 Thrombin signaling and protease-activated receptors PID, BioCarta PLCB1; MAP3K7 0.00346684 0.02636741 IL-17 signaling pathway Wikipathways MAP3K7; JAK1 0.00346684 0.02636741 Mucin type O-glycan biosynthesis - Homo sapiens (human) KEGG GALNTL6; WBSCR17 0.00369144 0.02636741 Propanoate metabolism - Homo sapiens (human) KEGG ACACA; ABAT 0.00446835 0.02978898 Toxoplasmosis - Homo sapiens (human) KEGG HLA-DRA; MAP3K7; JAK1 0.00517687 0.03235546 FAS pathway and stress induction of HSP regulation Wikipathways CASP10; MAP3K7 0.00600657 0.03533279 BMP receptor signaling PID BMP7; MAP3K7 0.00751491 0.0417495 Proton pump inhibitor pathway, Pharmacodynamics PharmGKB PLCB1; AKAP2 0.00816134 0.04295442 Chondroitin sulfate/dermatan sulfate metabolism Reactome VCAN; GPC5 0.00917629 0.04588144 Downstream TCR signaling Reactome HLA-DRA; MAP3K7 0.00966203 0.04600966 Hippo signaling pathway - Homo sapiens (human) KEGG BMP7; GLI2; CTNNA3 Input options: 13 databases, minimum overlap with input list = 2 and P value cutoff = 0.01. Bold styled value means overlapped pathway. GWAS, genome-wide association study. *Overlapped genes between the input and database sets; bold denotes a common pathway identified from the first and second GWASs. Table 3 Pathway-based gene sets identified from the second GWAS P value Q value Pathway Source Overlapped genes* 0.00010214 0.00423877 Fas signaling pathway (cd95) PID, BioCarta PTPN13; FADD; PRKDC 0.00085642 0.02369427 G protein signaling pathways Wikipathways PRKCH; PDE4B; PRKAR2B; PDE8B 0.00214336 0.03557973 Repression of pain sensation by the transcriptional regulator, DREAM PID, BioCarta PRKAR2B; CREM Input options: 13 databases, minimum overlap with input list = 2 and P value cutoff = 0.01. Bold styled value means overlapped pathway. GWAS, genome-wide association study *Overlapped genes between the input and database sets; bold denotes a common pathway identified from the first and second GWASs. DISCUSSION Previous GWASs have identified only SNPs located in the HLA region12 or SNPs in linkage disequilibrium with HLA alleles1314 as genetic risk factors for SCAR. In our study, causative drugs were diverse, and thus we could not adjust for the potential effects of the drug-specific HLA. Considering that the drug-specific HLA is believed to be the strongest genetic risk factor of SCAR, failure to adjust for its effects may partly explain why we did not detect any SNPs with genome-wide significance. It may be true that no single genetic variant confers an increased risk of SCAR greater than the drug-specific HLA. A gene set-based approach is widely used in genetic studies and has shown promising results. Likewise, we found that the gene sets annotated with the 200 top-ranked SNPs from t2 independent GWASs of SCAR induced by diverse drugs were commonly over-represented in the Fas signaling pathway. To the best of our knowledge, this is the first study to show that the Fas signaling pathway is a common genetic risk factor for SCAR across diverse drugs. Over-representation analysis statistically evaluates the fraction of genes in a particular pathway found among a set of genes (input).5 Over-representation analysis can be performed in ConsensusPathDB using gene sets and functional modules derived from incorporating pathway definitions provided by the source databases. For each predefined module, a P value is calculated based on the hypergeometric distribution, which reflects the significance of the observed overlap between the input gene list and the module's members compared to random expectations.15 If the input list is obtained from a case-control study, over-representation analysis may indicate pathways and functional sub-networks that are dysregulated in the disease state. The Fas receptor (CD95) mediates apoptosis via the Fas-ligand (FasL), which is expressed on the surface of other cells. The Fas-FasL interaction plays an important role in the immune system and a lack of this interaction leads to autoimmunity.16 Fas-induced keratinocyte apoptosis plays a key role in the pathogenesis of SCAR.1718 Therefore, our observations suggest that the genetic susceptibility associated with the Fas signaling pathway may be a common (drug-non-specific) risk factor for the development of SCAR in conjunction with the drug-specific HLA, although detailed genetic mechanisms should be investigated in a future study. Although genome-wide significance was not reached, some SNPs mapped to NKAIN2 (Sodium/Potassium Transporting ATPase Interacting 2 gene) and ANGPT2 (Angiopoietin 2 gene) showed significant associations with SCAR both in 2 GWASs. NKAIN2 gene encodes a transmembrane protein that interacts with the beta subunit of a sodium/potassium-transporting ATPase. A previous report showed that a chromosomal translocation involving this gene was a cause of lymphoma.19 Angiopoietin 2 encoded by ANGPT2 is an antagonist of angiopoietin 1 and thus disrupts its vascular remodeling ability which may induce endothelial cell apoptosis.20 However, so far, there has been no report to show the association between SCAR or drug adverse reaction and NKAIN2 or ANGPT2. This finding may emphasize an importance of the gene-set based analysis, although further mechanistic studies need to be performed. There are a few limitations to the present study. First, our genotyping platforms differed between the 2 GWASs and thus SNPs genotyped overlapped poorly, which increased the likelihood that SNPs that were significantly associated with SCAR were not genotyped across the studies. However, it is known that the probabilistic nature of imputed SNPs presents challenges when testing for association of those SNPs.21 Therefore, we used a gene-based approach instead of imputation of our SNP dataset. Secondly, a SNP was annotated to a gene when a gene harbored that SNP within its chromosomal region. A SNP located outside of a gene can be functionally relevant by affecting the expression of genes located adjacently or remotely.22 In addition, we selected only the 200 top-ranked SNPs for analysis. Therefore, it is possible that some functionally important genes were missed and thus not included in the over-representation analysis, which limited the identification of more significant pathways. Secondly, approximately half of the enrolled subjects (subjects with a diagnosis of DRESS) showed skin eruption rather than blister or detachment. The role of Fas-FasL in the pathogenesis of drug eruption is not clear. However, a previous report showed that overlapping expression of Fas and FasL is accompanied by apoptosis in fixed-drug eruption lesions.23 In addition, we found no difference between DRESS cases and all other cases in our over-representation analysis (data not shown). Thus, it is possible that the Fas signaling pathway plays its own role in the pathogenesis of eruption seen in subjects with DRESS. Another possible explanation is that prompt and intensive treatment may have stopped the progression from eruption to blister or detachment in the subjects enrolled in this study. Thirdly, a small number of patients with SCAR was also an important issue before generalizing our results. However, given that SCAR is a very rare event, our findings may provide investigators to prioritize variants for follow-up analysis. A pathway analysis based on the gene set is known to be particularly useful for pilot studies with small sample sizes.24 In conclusion, an over-representation analysis using gene sets from 2 independent GWASs of SCAR induced by diverse drugs identified the Fas signaling pathway as a significant and common pathway. Our observations suggest that the Fas signaling pathway may be a common genetic risk factor for SCARs across diverse drugs. Disclosure: There are no financial or other issues that might lead to conflict of interest. SUPPLEMENTARY MATERIALS Supplementary Table S1 A total of 62 SNPs and their annotated genes in the first GWAS used for an over-representation analysis Supplementary Table S2 A total of 69 SNPs and their annotated genes in the second GWAS used for an over-representation analysis Supplementary Table S3 The top-ranked 200 SNPs identified from the first and second GWAS Supplementary Table S4 Genes significantly over-represented in the Fas signaling pathway Supplementary Fig. S1 The Manhattan plots of each GWAS. (A) The first GWAS. (B) The second GWAS. ==== Refs 1 Chung WH Wang CW Dao RL Severe cutaneous adverse drug reactions J Dermatol 2016 43 758 766 27154258 2 Pavlos R Mallal S Phillips E HLA and pharmacogenetics of drug hypersensitivity Pharmacogenomics 2012 13 1285 1306 22920398 3 Jung JW Song WJ Kim YS Joo KW Lee KW Kim SH HLA-B58 can help the clinical decision on starting allopurinol in patients with chronic renal insufficiency Nephrol Dial Transplant 2011 26 3567 3572 21393610 4 McCormack M Alfirevic A Bourgeois S Farrell JJ Kasperavičiūtė D Carrington M HLA-A*3101 and carbamazepine-induced hypersensitivity reactions in Europeans N Engl J Med 2011 364 1134 1143 21428769 5 Khatri P Sirota M Butte AJ Ten years of pathway analysis: current approaches and outstanding challenges PLoS Comput Biol 2012 8 e1002375 22383865 6 Creamer D Walsh SA Dziewulski P Exton LS Lee HY Dart JK U.K. guidelines for the management of Stevens-Johnson syndrome/toxic epidermal necrolysis in adults 2016 Br J Dermatol 2016 174 1194 1227 27317286 7 French LE Toxic epidermal necrolysis and Stevens Johnson syndrome: our current understanding Allergol Int 2006 55 9 16 17075281 8 Shiohara T Iijima M Ikezawa Z Hashimoto K The diagnosis of a DRESS syndrome has been sufficiently established on the basis of typical clinical features and viral reactivations Br J Dermatol 2007 156 1083 1084 17381452 9 World Health Organization, Uppsala Monitoring Centre The use of the WHO-UMC system for standardized case causality assessment Uppsala The Uppsala Monitoring Centre 2005 10 Purcell S Neale B Todd-Brown K Thomas L Ferreira MA Bender D PLINK: a tool set for whole-genome association and population-based linkage analyses Am J Hum Genet 2007 81 559 575 17701901 11 Kamburov A Wierling C Lehrach H Herwig R ConsensusPathDB--a database for integrating human functional interaction networks Nucleic Acids Res 2009 37 D623 D628 18940869 12 Génin E Schumacher M Roujeau JC Naldi L Liss Y Kazma R Genome-wide association study of Stevens-Johnson Syndrome and Toxic Epidermal Necrolysis in Europe Orphanet J Rare Dis 2011 6 52 21801394 13 Ozeki T Mushiroda T Yowang A Takahashi A Kubo M Shirakata Y Genome-wide association study identifies HLA-A*3101 allele as a genetic risk factor for carbamazepine-induced cutaneous adverse drug reactions in Japanese population Hum Mol Genet 2011 20 1034 1041 21149285 14 Tohkin M Kaniwa N Saito Y Sugiyama E Kurose K Nishikawa J A whole-genome association study of major determinants for allopurinol-related Stevens-Johnson syndrome and toxic epidermal necrolysis in Japanese patients Pharmacogenomics J 2013 13 60 69 21912425 15 Kamburov A Pentchev K Galicka H Wierling C Lehrach H Herwig R ConsensusPathDB: toward a more complete picture of cell biology Nucleic Acids Res 2011 39 D712 D717 21071422 16 Brunner T Wasem C Torgler R Cima I Jakob S Corazza N Fas (CD95/Apo-1) ligand regulation in T cell homeostasis, cell-mediated cytotoxicity and immune pathology Semin Immunol 2003 15 167 176 14563115 17 French LE Trent JT Kerdel FA Use of intravenous immunoglobulin in toxic epidermal necrolysis and Stevens-Johnson syndrome: our current understanding Int Immunopharmacol 2006 6 543 549 16504917 18 Su SC Chung WH Cytotoxic proteins and therapeutic targets in severe cutaneous adverse reactions Toxins (Basel) 2014 6 194 210 24394640 19 Zhao SC Zhou BW Luo F Mao X Lu YJ The structure and function of NKAIN2-a candidate tumor suppressor Int J Clin Exp Med 2015 8 17072 17079 26770299 20 Augustin HG Koh GY Thurston G Alitalo K Control of vascular morphogenesis and homeostasis through the angiopoietin-Tie system Nat Rev Mol Cell Biol 2009 10 165 177 19234476 21 Marchini J Howie B Genotype imputation for genome-wide association studies Nat Rev Genet 2010 11 499 511 20517342 22 Westra HJ Franke L From genome to function by studying eQTLs Biochim Biophys Acta 2014 1842 1896 1902 24798236 23 Choi HJ Ku JK Kim MY Kang H Cho SH Kim HO Possible role of Fas/Fas ligand-mediated apoptosis in the pathogenesis of fixed drug eruption Br J Dermatol 2006 154 419 425 16445769 24 Kao PY Leung KH Chan LW Yip SP Yap MK Pathway analysis of complex diseases for GWAS, extending to consider rare variants, multi-omics and interactions Biochim Biophys Acta 2017 1861 335 353 27888147