Adverse drug reactions are a major public health issue as they represent an important cause of morbidity and mortality 1 . Skin lesions are frequent expressions of adverse drug reactions with, for some drugs, up to 10% of cutaneous reactions observed 2 3 . SJS and TEN are severe cutaneous adverse reactions characterized by the development of acute exanthema which progresses towards limited (in SJS) or widespread (in TEN) blistering and erosion of the skin and mucous membranes 4 . SJS and TEN are thus considered to be two different forms of the same disease with TEN representing the most severe form 5 . The incidence of SJS/TEN is estimated to be of approximately 1-2 patients per million individuals per year 6 . It is thus a very rare disease but with a high morbidity and mortality (reaching up to 45% in TEN) that requires intensive treatment.

SJS/TEN is not associated with a single drug or a single group of drugs but several different drugs have been involved. A limited number of drugs however are more often associated with the disease: antibacterial sulfonamides (especially sulfamethoxazole), allopurinol, which is the most frequent drug involved 7 , carbamazepine, lamotrigine, phenobarbital, phenytoine, non-steroidal anti-inflammatory drugs (NSAIDs) of the oxicam type and neviparine 4 8 . Only a small number of individuals exposed to these "highly suspected" drugs develop the disease and a genetic susceptibility has been suggested 9 10 11 12 .

An association with HLA was reported more than 20 years ago 12 13 . More recently, studies in the Han Chinese population have involved the HLA-B locus with very strong drug-specific associations: the HLA-B*1502 allele was found in all carbamazepine-induced SJS/TEN patients 14 and the HLA-B*5801 in all allopurinol-induced SJS/TEN patients 15 . An investigation of HLA-B associations in European samples did not detect the association between the HLA-B*1502 allele and carbamazepine-induced SJS/TEN but did report a strong but not complete association between HLA-B*5801 and the allopurinol-induced disease 16 17 . Furthermore, it has also been suggested that HLA genetic predisposition may not be the same for SJS and TEN and that it might thus be important to take into account the disease severity in association tests 18 .

Apart from HLA, several other candidate genes have been tested for association with SJS/TEN, mostly genes involved in the immune response, in the inflammation process or in drug metabolism 18 19 20 21 . However, no association with these genes has been consistently found and it has been suggested that rather than focusing on candidate genes, a genome-wide association study might provide more insights into the genetic susceptibility of adverse drug reactions 21 .

In the context of the RegiSCAR project (European Registry of Severe Cutaneous Adverse Reactions to Drugs and Collection of Biological Samples) DNA of 563 cases of SJS/TEN was collected. This is the largest available sample of SJS/TEN patients in the world with accurate medical information regarding the severity of the disease and the history of drug intake. These patients were genotyped at the Centre National de Genotypage (CNG) using Illumina 317 K chips and a Genome Wide Association Study (GWAS) was conducted against controls selected from the CNG European Reference Control Panel 22 .

The clinical details of the 424 patients remaining after quality-control are summarized in Table 1: 61.3% are female, 47.2% have a SJS phenotype, 34.4% a SJS/TEN overlap phenotype and 18.4% a TEN phenotype. A majority of the patients was sampled from France (110 patients) and Germany (277 patients). History of drug intake was carefully monitored and a subgroup of 57 patients (13.4%) with potentially allopurinol-induced reactions was identified. There is no difference in phenotype distribution between males and females but there are some differences depending on the drug and the country where samples were collected. There is a lack of TEN among allopurinol-induced patients (test of homogeneity of the three disease forms in allopurinol induced versus non-allopurinol induced partients: chi-square = 7.89, p-value = 0.019) and a lack of TEN among patients from Germany as compared to patients from France (chi-square = 23.68, p-value = 7.2 10-6). The two variables however are not independent since allopurinol is more used in Germany than in some other countries and in particular France.

For association testing, the 424 patients were compared to the 1,881 controls from the reference panel originating from France (1,228 individuals) and Germany (653 individuals) who were also the closest controls to the cases on the PCA plot (Figure 1). We then tested each SNP for association controlling on the first two PCs of the PCA using Eigenstrat 24 . Results are presented in Figure 3. One genome-wide significant signal was detected in the MHC region on chromosome 6 where 6 SNPs showed False Discovery Rate adjusted q-values below 5% (see Figure 4 and Table 2). The most significant SNP, rs9469003, is located at position 31,515,807, ~85 kb upstream of the HLA-B locus and has an OR of 1.73 (95%CI = [1.44;2.08]). The five other SNPs are all located in a more telomeric region, 250 kb apart from rs9469003, in between positions 31,114,834 and 31,250,224 (see Figure 5). There is very limited Linkage Disequilibrium (LD) between these five SNPs and the top one and the six SNPs defined 13 haplotypes with a frequency of 1% or above. The haplotype CACGAC formed by the risk allele at each locus showed the strongest association with the disease with an OR of 2.84 (95%CI = [2.03; 3.98]) that is significantly higher than the one observed for rs9469003 only. To dissect further this haplotypic association, the haplotype conditioning test implemented in Plink v1.06 was used and we found that none of the 6 SNPs by itself has an independent effect nor explains the haplotype association.

Using a multinomial logistic regression approach, we found that the association with haplotype CACGAC was much stronger in the subgroup of patients with an allopurinol induced disease (ORallopurinol = 7.77, 95%CI = [4.66; 12.98]) than in other patients (ORotherdrugs = 1.92, 95%CI = [1.40;2.64]). The equality of odds-ratio (ORallopurinol = ORotherdrugs) was strongly rejected (p = 6.56 × 10^-7). Results were the same in France and Germany a shown in Additional File 1 Tables S1 and S2.

To determine how much of the observed signal could be explained by the known association between allopurinol-induced SJS/TEN and HLA-B*5801, we took advantage of the availability of HLA-B two-digit resolution genotypes for 74 of the 424 patients 16 . Among the 74 HLA-B genotyped patients, 11 were carriers of an HLA-B58 allele and had an allopurinol-induced disease. The frequency of the CACGAC haplotype is increased in HLA-B58 carriers as compared to non HLA-B58 carriers (32.44% versus 10.44%, p-value of the test adjusted on the top 2 PC = 0.0052) but the linkage disequilibrium is not complete. Note however that given the small sample size of HLA-B typed individuals, the power was quite low to detect this difference (power of 60% at a type-one error rate level of 5%).

Discovering genes involved in severe cutaneous adverse reactions and especially SJS/TEN is a major challenge for pharmacogenetics as these reactions, when not fatal, are a sword of Damocles for those who already had the disease and fear to take any drug. Despite important efforts, only genes located in the HLA region have been identified so far. One possible explanation for this lack of success is the limited sample sizes of patients available. Indeed SJS/TEN is fortunately a very rare disease and in most studies, sample sizes rarely exceed one or two hundreds of patients, making it difficult to investigate more than a few candidate genes. Through a collaborative effort, the RegiSCAR group was able to collect detailed medical information and DNA of more than half a thousand of patients from Europe who were genotyped on Illumina 317 K chips. By comparing their genotypes to the ones of 1,881 controls genetically matched for the country of origin, we were able to study the SJS/TEN association at a genome-wide level on European samples. Apart from six SNPs located in the HLA region, no other locus was found to be associated with the disease at a high enough significance level to ensure it is not false positive due to the multiple tests performed. Sample sizes of half a thousand are not enough to ensure a good power to detect common variants with small effects (OR in the range between 1.1 and 1.3) similar to those identified in several common diseases such as diabetes or different cancers 29 where larger samples are easier to collect. However, given the 424 patients remaining after Quality Control, the study was powerful enough to detect common variants with modest effects (the power exceeded 80% to detect variants with an allele frequency above 15% conferring ORs above 1.7 under a multiplicative model). The fact that only the HLA region is detected, suggests that there might not be any other common variant that confer a substantial increase in disease risk and could thus be of interest as a predictive factor. These power computations might however be too optimistic as they are based on the assumption that the same genetic variants might be involved in cases of SJS/TEN associated to different drugs. If we consider that there might exist some levels of genetic heterogeneity depending on the drugs involved as it is the case for HLA associations, the power can be dramatically reduced. For example, considering only the 57 allopurinol-exposed cases, the power to detect a variant of frequency 15% conferring an OR of 2 is then only 5% at a nominal type-one error rate value of 10^-6. On the other hand, the power to detect a similar effect than the one observed at rs9469003 (i.e., risk allele of frequency of 15% with an associated OR of ~4 among allopurinol-induced SJS/TEN) is 99%, showing that the study was not underpowered to detect strong effects even if they were restricted to a small subset of patients.

It is true however that we might have missed some variants with important effects, especially those with minor allele frequencies below 5% that are not well covered by SNP-chips and it could thus be of interest to investigate the association with rare variants. If the "common disease-common variant hypothesis" was believed to explain the susceptibility to multifactorial diseases 30 31 , results from GWAS have shown that common variants only explain a minor part of most common disease heritabilities. It has been suggested that rare variants could explain at least part of this "missing heritability" and indeed, rare variants have been found in several diseases 32 33 34 35 36 . Interestingly these rare variants are often functional variants with a direct impact on the protein functionality and they usually confer a much stronger increase in disease risk than common variants. They are also more likely to be affected by some moderate levels of negative selection 37 38 . For diseases such as adverse drug reactions, it is not unlikely that the genes and variants involved could be under selection. Indeed if most of the drugs have been introduced too recently to be directly responsible for selective pressures, they are often derived from nutrients that have been consumed by humans for a long time. Thus, genes involved in drug metabolism can potentially be involved in natural selection and this was confirmed by previous studies 39 40 41 . However, the implication of drug metabolism genes in SJS/TEN remains to be established.

Concerning the association with HLA detected in this sample, we confirmed that it is drug-specific with the strongest association found for the group of patients where allopurinol is suspected to be the cause of the disease. However the disease association was still detectable after exclusion of all patients exposed to allopurinol (whether considered as causing the reaction or not) (OR = 2.13, 95% CI = [1.41; 3.23]). This suggests that different HLA alleles are probably involved depending on the drug, and the 6 SNPs identified could be those that were in linkage disequilibrium (LD) with several of these HLA alleles. To further explore this hypothesis, we tried to impute HLA-B genotypes from the SNP data 42 but we were not successful because of the presence of multiple rare alleles at this locus such as HLA-B58 itself that is only present in one of the European families genotyped by de Bakker et al. 43 precluding the possibility to find appropriate tags for this allele using this sample. To overcome the problem, we tried to determine tag SNPs using our own HLA-B genotyped sample that was enriched in HLA-B58 but we did not find any SNP or group of SNPs that was correlated enough with the HLA-B58 allele (the maximum r² value was below 20%). This poor HLA-B tagging ability of the SNPs could certainly explain why the association found here in relation with allopurinol is much weaker than the one reported in studies on European samples using resolved HLA-B alleles 16 44 . However, even in these two latter studies, the association with HLA-B*5801 is much weaker than the one found in Han Chinese 15 or Thai populations 45 where all patients with allopurinol-induced SJS/TEN are carriers of this allele. It is thus possible that the HLA genetic determinants involved in SJS/TEN are not the HLA-B alleles themselves but some loci in LD with them. Further studies in this genomic region would certainly be necessary to better understand the mechanisms involved.

On a methodological point of view, this study illustrates how Reference Control Panels could be used to test for association when only cases are genotyped. Only patients were genotyped here and controls were selected from the CNG European Control Panel 22 based on their country of origin that was either France and Germany, since those were the two main countries where the cases come from. We also tried to use the whole panel of controls and correct for population stratification using Principal Component Analysis (PCA) and an adjustment on up to 20 PCs instead of the 2 used here. We found that including controls that were more genetically distant to the cases does not really improve the signal in the HLA region and instead could lead to false positive results in other genomic regions. This is somewhat in contradiction with the results recently reported by Zhuang et al. 46 and will probably need to be further investigated. We also used some alternative approaches where controls were genetically matched to cases using either a distance based on the identity by state as described by Guan et al. 47 or a distance derived from the top PCs of the PCA performed to identify axes of variations similar to the one proposed by Luca et al. 48 . Tests accounting for the matching units were then performed and we found that only the HLA region pointed out with significance levels in the same order of magnitude as the ones obtained here. These matching strategies might however be of interest as they could allow the analysis of more heterogeneous samples by including for example those non-European RegiSCAR cases that were excluded here and matching them to individuals from the CEPH-Human Genome Diversity Panel 49 .

Our study confirms the involvement of genetic variants located in the HLA region in the susceptibility to SJS/TEN in European samples, especially in association with allopurinol. No difference is seen depending on the disease severity and no other locus reaches genome-wide association in this sample that is also the largest one collected so far.

The authors declare that they have no competing interests.

MM, MS, JCR, AH, ML designed the study, DZ, ML supervised the genotyping, LN, YL, JCR, MM, PS were responsible for clinical data collection, RK, EG performed the statistical analysis, EG wrote the manuscript and all authors contributed to and approved the final draft.