Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Genome-Wide Association Study in Immunocompetent Patients with Delayed Hypersensitivity to Sulfonamide Antimicrobials

  • Jennifer M. Reinhart,

    Affiliation Department of Medical Sciences, School of Veterinary Medicine, University of Wisconsin-Madison, Madison, Wisconsin, United States of America

  • Alison Motsinger-Reif,

    Affiliation Bioinformatics Research Center, Department of Statistics, North Carolina State University, Raleigh, North Carolina, United States of America

  • Allison Dickey,

    Affiliation Bioinformatics Research Center, Department of Statistics, North Carolina State University, Raleigh, North Carolina, United States of America

  • Steven Yale,

    Current address: North Florida Regional Medical Center, Gainesville, Florida, United States of America

    Affiliation Marshfield Clinic Research Foundation, Marshfield, Wisconsin, United States of America

  • Lauren A. Trepanier

    lauren.trepanier@wisc.edu

    Affiliation Department of Medical Sciences, School of Veterinary Medicine, University of Wisconsin-Madison, Madison, Wisconsin, United States of America

Abstract

Background

Hypersensitivity (HS) reactions to sulfonamide antibiotics occur uncommonly, but with potentially severe clinical manifestations. A familial predisposition to sulfonamide HS is suspected, but robust predictive genetic risk factors have yet to be identified. Strongly linked genetic polymorphisms have been used clinically as screening tests for other HS reactions prior to administration of high-risk drugs.

Objective

The purpose of this study was to evaluate for genetic risk of sulfonamide HS in the immunocompetent population using genome-wide association.

Methods

Ninety-one patients with symptoms after trimethoprim-sulfamethoxazole (TMP-SMX) attributable to “probable” drug HS based on medical record review and the Naranjo Adverse Drug Reaction Probability Scale, and 184 age- and sex-matched patients who tolerated a therapeutic course of TMP-SMX, were included in a genome-wide association study using both common and rare variant techniques. Additionally, two subgroups of HS patients with a more refined clinical phenotype (fever and rash; or fever, rash and eosinophilia) were evaluated separately.

Results

For the full dataset, no single nucleotide polymorphisms were suggestive of or reached genome-wide significance in the common variant analysis, nor was any genetic locus significant in the rare variant analysis. A single, possible gene locus association (COL12A1) was identified in the rare variant analysis for patients with both fever and rash, but the sample size was very small in this subgroup (n = 16), and this may be a false positive finding. No other significant associations were found for the subgroups.

Conclusions

No convincing genetic risk factors for sulfonamide HS were identified in this population. These negative findings may be due to challenges in accurately confirming the phenotype in exanthematous drug eruptions, or to unidentified gene-environment interactions influencing sulfonamide HS.

Introduction

Potentiated sulfonamides, such as sulfamethoxazole (SMX) in combination with trimethoprim (TMP), are effective antibiotics for the treatment of urinary infections, bronchitis, pneumonia, and methicillin-resistant S. aureus (MRSA) infections. [14] TMP-SMX is also the drug of choice for the prevention of opportunistic protozoal infections, such as toxoplasmosis and Pneumocystis jiroveci pneumonia, in immunocompromised patients, particularly in those with AIDS. [58]

Despite this broad-spectrum of activity, clinical use of TMP-SMX is limited by the development of idiosyncratic, delayed-type hypersensitivity (HS) reactions including fever and cutaneous drug eruptions and, less commonly, multi-organ dysfunction. [9] In fact, TMP-SMX is the leading cause of cutaneous hypersensitivity reactions, and the most common culprit in severe bullous skin eruptions such as Stevens-Johnson syndrome and toxic epidermal necrolysis, which carry up to 30% mortality. [1014] Hypersensitivity reactions occur in approximately 3% of TMP-SMX-treated patients in the general population. [15,16] Both clinical reports and ex vivo drug challenge studies have suggested that this syndrome may be familial; thus a genetic basis for sulfonamide HS has been proposed. [17,18]

Pharmacogenetic risk for sulfonamide HS has only been evaluated in a limited number of studies. Unfortunately, most of these studies were in HIV-positive patients, for whom acquired, rather than genetic, risk is more likely. [1923] Genetic studies have been sparse in immunocompetent patients and have focused predominately on enzymes responsible for SMX biotransformation or on glutathione pathways, which can neutralize reactive drug metabolites. Two reports suggested that “slow” acetylator NAT2 genotypes were associated with sulfonamide HS. [24,25] However, in a larger study in our laboratory, we found no such association despite full resequencing of NAT2 coding region. [26] In that same study we also investigated the sulfamethoxazole detoxification genes, CYB5A and CYB5R3, but polymorphism frequencies were low and were not significantly different between sulfonamide HS and tolerant patients. [26]

CYP2C9*2/*3 variant alleles have been associated with decreased production of the reactive, SMX-hydroxylamine (SMX-HA) metabolite in human liver microsomes; [27] however, allele frequency variation has not been investigated in immunocompetent patients with sulfonamide HS. [20] Further, an increased risk for cutaneous drug eruptions overall was reported for patients with GSTM1 and GSTT1 null genotypes, but this study population was mixed, with sulfonamide HS patients representing only 4/36 affected patients. [28] Finally, an association between an HLA-A30 haplotype and sulfonamide-induced skin eruptions was found in a Turkish population, but HLA genotypes have yet to be studied in other immunocompetent populations. [12]

Beyond these few candidate-gene studies, no other drug biotransformation, redox, or immunoregulatory genes have been evaluated for an association with sulfonamide HS. Genome-wide association studies (GWAS) have successfully identified genetic targets associated with a variety of other adverse drug reactions when a candidate gene approach has proved insufficient. [2932] Therefore, the purpose of this study was to screen for genetic markers of sulfonamide HS in an immunocompetent population, compared to drug tolerant controls, using a GWAS design.

Methods

Patient Identification

Patients with delayed sulfonamide HS and patients tolerant of a course of sulfonamide antibiotics were identified through the Marshfield Clinic Research Foundation’s Personalized Medicine Research Project (PMRP), a cohort of over 20,000 patients who receive their medical care at the Marshfield Clinic, Marshfield, WI. These patients have previously provided informed consent for the use of their samples for biomedical research, linked to de-identified medical records data. [33] Medical records were searched electronically for a history of TMP-SMX administration or for a diagnosis of sulfonamide HS. Marshfield Clinic Research Foundation staff individually reviewed medical records, using a structured abstraction form to identify patients with sulfonamide HS. Each case was adjudicated to ensure consistency and accuracy. The abstraction form included the following eligibility criteria: (1) administration of TMP-SMX for at least 5 days prior to the adverse event; [9] (2) documentation of one or more new clinical signs after starting TMP-SMX, including fever with or without eosinophilia, skin rash, increases in liver enzyme activities, hyperbilirubinemia, blood dyscrasias (anemia, leukopenia or thrombocytopenia), pneumonitis, myocarditis, aseptic meningitis, polyarthritis, acute interstitial nephritis, toxic epidermal necrolysis, or Stevens-Johnson syndrome; [9] (3) lack of other clinical explanation for the adverse event; and (4) resolution of clinical signs with discontinuation of TMP-SMX alone. Patients with only gastrointestinal symptoms such as nausea, vomiting or diarrhea, [9] or with acute anaphylactoid reactions, [34,35] were excluded. Because some forms of immunosuppression, in particular AIDS, lead to a high acquired risk of SMX hypersensitivity, apparently independent of genotype, [20] immunocompromised patients, including those with HIV infection or undergoing immunosuppressive chemotherapy, were not eligible. These criteria together were designed to yield a score of 6 or more, or “probable” adverse reaction, using the Naranjo Adverse Drug Reaction scale. [36]

Control patients (“tolerant;” TOL) within the PMRP that were prescribed TMP-SMX were enrolled sequentially from medical records in random order, to provide a 2:1 match with the HS patients for race, gender, and decade of age at sulfonamide treatment. TOL patients must have been prescribed a course of TMP-SMX at a standard therapeutic daily dosage for at least 10 days, with adequate follow-up in the medical record to indicate that the drug was taken and tolerated without adverse event. Clinical and demographic variables, including body weight, dosage, duration of treatment, and reason for TMP-SMX prescription were also abstracted. Patient data were provided to the investigators in a de-identified format, and therefore the study protocol was reviewed and granted exemption from federal regulations by the UW Health Sciences human subjects Minimal Risk Institutional Review Board. One investigator (SY) from the PRMP did have direct access to patient information as a part of the validation process but the review board was aware of this at the time the exemption was granted.

Genotyping, Quality Control, and Data Pre-analysis

DNA samples from 99 HS and 198 TOL patients from the PMRP were genotyped using one or more of the following platforms: HumanCoreExome, Illumina 660Hg18, Illumina 660Hg19, or Infinium (Illumina; San Diego, CA). Genotype-calls were performed in Genome Studio software (Illumina; San Diego, CA). For the HumanCoreExome, 660Hg18, and 660Hg19 platforms, the data was formatted for use in PLINK (http://pngu.mgh.harvard.edu/purcell/plink/). [37,38] The genotyped Infinium raw data files existed in the A/B format and the single nucleotide polymorphisms (SNPs) were converted to their corresponding alleles to match the PLINK format of the other platforms. SNPs in the Infinium dataset whose genotypes were denoted as either an indel (12,602) or had ambiguous strand designations (1,286) were removed. PLINK formatted files for the HumanCoreExome, 660Hg18, and 660Hg19 platforms are available on dbGaP, along with the raw text files for the Infinium platform.

For each platform, called data were filtered for quality control (QC) in PLINK. No patients were removed from the study due to low genotyping (efficiency < 90%). SNPs were excluded if the genotyping efficiency was < 98% or if the minor allele frequency (MAF) was < 0.01. [39] Trend tests for Hardy-Weinberg equilibrium were performed and variants with extreme deviations were also excluded from analyses. Q-Q plots were used to assess overall quality of the GWAS data, and the quality of any significant SNPs was evaluated manually after association analysis.

Data from the four platforms were combined to form a single dataset. However, few SNPs were shared among all four platforms. Therefore, to maximize the number of SNPs evaluated, only data from the HumanCoreExome and Infinium platforms were included in the final analysis. This resulted in the exclusion of 10 patients (4 HS, 6 TOL) who were only genotyped on the Illumina Hg660 systems. Imputation was performed on the combined dataset using SHAPEIT (https://mathgen.stats.ox.ac.uk/genetics_software/shapeit/shapeit.html#home) and IMPUTE2 (https://mathgen.stats.ox.ac.uk/impute/impute_v2.html) software. The EUR group of individuals (n = 379) in the 1000 Genomes Project reference panel (June 2014 release; 1000genomes.org) was used to phase the haplotypes because the majority of the patients in this study were Caucasian. 420 SNPs in the dataset were either missing in the reference panel or had strand incompatibilities and were removed. The number of genotyped SNPs used in the imputation was 226,516 and the resulting number of imputed SNPs that had an info score > 0.7 and that were used in the downstream analyses was 10,329,316.

Including relatives within the same association study can skew results, so the KING program (http://people.virginia.edu/~wc9c/KING/) was used to assess relatedness between subjects. Patients were removed if a relatedness of the 3rd degree or less was identified, resulting in the exclusion of an additional 10 patients (3 HS, 7 TOL). Further population substructure was assessed using EIGENSOFT software (http://www.hsph.harvard.edu/alkes-price/software/) with a linkage disequilibrium pruned set (62,653) from the genotyped SNPs; principal component analysis identified two clear outliers (1 HS, 1 TOL) from the population, which were removed. The final dataset included 275 patients: 91 HS and 184 TOL.

Common Variant Analysis

Multivariable logistic regression analysis was performed on the genotyped SNPs using PLINK and using SNP dosages for the imputed SNPs. 2-bit dummy-encoding was used to ensure that no genetic model assumptions were made. [39] The first four principle components, determined in EIGENSOFT, were included as covariates to minimize spurious associations due to underlying population substructure. The population was almost entirely Caucasian (99%), so race was not included as a covariate. Similarly, sex and age were not found to be significant covariates (coefficient p-value < 0.05), so were not included in the final model. Before calculating the genomic inflation factor (λ) for the combined genotyped and imputed common variant results, any SNPs with a MAF < 0.01 were removed along with any imputed SNPs that also existed in the genotyped set. This resulted in 226,505 genotyped SNPs and 8,066,491 imputed SNPs. λ was 1.03 and the p-values were adjusted to correct for it. The Q-Q plot is shown in S1 Fig. All SNPs were assessed in the final analysis, but those within candidate genes previously hypothesized to be involved in sulfonamide HS pathogenesis (Table 1) were of particular interest. Based on the number of comparisons including both the genotyped and the imputed SNPs, p ≤ 6.03 x 10−7 was considered suggestive of significance and p ≤ 6.03 x 10−9 was considered significant at the genome-wide level. The Bonferroni correction was used for multiple comparison corrections at an alpha level of 0.05.

thumbnail
Table 1. Candidate genes within the GWAS hypothesized to be involved in the pathogenesis of sulfonamide HS, and the mechanistic rationale for inclusion of each candidate gene.

https://doi.org/10.1371/journal.pone.0156000.t001

Assuming a multiplicative model of inheritance and a risk-allele frequency of 0.5, an a priori power calculation demonstrated that an available sample of 99 HS and 198 TOL patients would yield at least 80% power to detect a genetic risk model with a relative risk of ≥ 3.0, while maintaining an overall family-wise error rate of 10−8 (representing a correction at a genome-wide level).

Rare Variant Analysis

The present study was not powered to detect a significant difference in allele frequencies for SNPs with low MAFs. [38] However, associations between low-MAF calls and various diseases have been previously identified, emphasizing the potential importance of these rare variants. [4042] Therefore, we characterized the influence of rare variants using a combined multivariate and collapsing gene-level approach. [43,44] All SNPs with MAF < 0.03 within an individual gene were collapsed into a single covariate and each gene was included in a multivariate logistical regression model. 17,793 loci were included in this analysis and candidate genes were specifically evaluated. Several rare-variant burden tests (ind, prop, weight) and variance-component tests (skat, skat-o) were used including those that weighted for possible functional significance of individual SNPs. [45, 46] For each locus included in the analysis, p-values from the separate rare-variant tests were combined to yield a single p-value using Lancaster’s method. [47] Gene-level associations with model adjusted q-value ≤ 0.05 were considered significant.

Subgroup Analyses

Because definitive diagnosis of sulfonamide HS is difficult and misclassification can occur, HS patients matching a more refined clinical phenotype were included in two separate subgroup analyses. The first subgroup (FEV-RASH) included patients exhibiting both a fever (body temperature > 98.8°F) and cutaneous drug eruption (rash) (n = 16). The second subgroup (FEV-RASH-EOS) included patients exhibiting fever, rash, and documented peripheral eosinophilia (eosinophils > 0.5 M/l) on complete blood count (n = 8). Each subgroup was compared with the full control group (n = 184) in both a common and a rare variant analysis. The genotyped (226,516) and imputed SNPs (10,329,316) that were used as input for the full set common and rare variant analyses were also used as input for the subgroup analyses. After quality filtering the common variant analysis results, 226,444 genotyped and 8,077,026 imputed SNPs remained for subgroup FEV-RASH. λ was 0.87 and the resulting Q-Q plot can be seen in S2 Fig. 17,766 loci were included in the rare variant analysis. For subgroup FEV-RASH-EOS, 226,415 genotyped and 8,078,605 imputed SNPs remained after the common variant analysis results were quality filtered. λ was 0.68 and the resulting Q-Q plot can be seen in S3 Fig. 17,778 loci were included in the rare variant analysis in this subgroup. Candidate genes were specifically evaluated in both the common and rare variant analyses.

Results

Patient Population

Ninety-one HS and 184 TOL patients were included in the final analysis. Inclusion and exclusion of subjects are represented in Fig 1. Patient demographics are summarized in Table 2. HS patients were predominantly Caucasian (98%) and female (82%), with a median age of sulfonamide treatment of 38 years (range 1–80). Most of the HS patients had exanthematous (95%) or bullous (1%) drug eruptions. Other clinical manifestations included fever (20%), eosinophilia (17%), thrombocytopenia (11%), neutropenia (5%), anemia (4%), or elevated liver enzymes/hyperbilirubinemia (3%). Matched TOL patients were predominantly Caucasian (99%) and female (83%), with a median age of sulfonamide treatment of 41 years (range 8–88).

thumbnail
Table 2. Demographic information for sulfonamide hypersensitive (HS) and drug-tolerant (TOL) patients.

Continuous data are presented as mean ± standard deviation (range). FEV-RASH = hypersensitive subgroup with fever and rash; FEV-RASH-EOS = hypersensitive subgroup with fever, rash, and eosinophilia. Diagnosis represents the underlying rationale for TMP/SMX prescription. (Note: body weight and therefore total daily dose in mg/kg were not available for 58 patients.)

https://doi.org/10.1371/journal.pone.0156000.t002

Common Variant Analysis

The Manhattan plot of the full dataset is shown in Fig 2. No SNP was suggestive of, or reached, genome-wide significance. Similarly, no SNP in any candidate gene reached significance. The SNP with the lowest p-value was rs160978 (p = 4.98 x 10−6) located in an intergenic region of chromosome 5.

thumbnail
Fig 2. Manhattan plot for common variant analysis of full set of sulfonamide hypersensitive (HS) vs. drug tolerant (TOL) patients.

No SNPs were suggestive of (p ≤ 6.03 x 10−7) or reached genome-wide significance (p ≤ 6.03 x 10−9).

https://doi.org/10.1371/journal.pone.0156000.g002

Rare Variant Analysis

Of the 17,793 loci included in this gene-level analysis, there were no statistically significant differences between HS and TOL patients in the full data set. In particular, there were no significant differences for any of the candidate genes. The gene with the lowest q-value was FNBP1 (q = 0.2558), located on chromosome 9, which codes for formin binding protein 1.

Subgroup Analyses

The 16 patients included in subgroup FEV-RASH were predominantly Caucasian (15/16) and female (10/16) with a median age of 42 years (range 1–70). In addition to fever and rash, 7 had blood dyscrasias (neutropenia, thrombocytopenia, and/or anemia) and 1 had bullous skin eruptions. In the common variant analysis, no SNP was suggestive of or reached genome-wide significance (Fig 3). However, in the rare variant analysis, there was a statistically significant difference between the TOL and FEV-RASH groups for the gene COL12A1 (q = 0.0167) on chromosome 6, which codes for type XII α1 collagen.

thumbnail
Fig 3. Manhattan plot for common variant analysis of the HS subgroup with rash and fever (FEV-RASH) vs. all TOL patients.

No SNPs were suggestive of (p ≤ 6.02 x 10−7) or reached genome-wide significance (p ≤ 6.02 x 10−9).

https://doi.org/10.1371/journal.pone.0156000.g003

The 8 patients in subgroup FEV-RASH-EOS were predominantly Caucasian (7/8) with a median age of 48 years (range 24–70). Males (4/8) and females (4/8) were equally distributed. Six of these patients had blood dyscrasias and 1 exhibited bullous skin eruptions. In the common variant analysis, no SNP was suggestive of, or reached, genome-wide significance (Fig 4). In the rare variant analysis, there were no significant differences between TOL and FEV-RASH-EOS patients for any locus, including COL12A1.

thumbnail
Fig 4. Manhattan plot for common variant analysis of HS subgroup with fever, rash and eosinophilia (FEV-RASH-EOS) vs. all TOL patients.

No SNPs were suggestive of (p ≤ 6.02 x 10−7) or reached genome-wide significance (p ≤ 6.02 x 10−9).

https://doi.org/10.1371/journal.pone.0156000.g004

No SNP or locus within the candidate gene set was significant for either subgroup in either the common or rare variant analysis. The top 100,000 SNP associations in the common variant analyses for the full data set and each subgroup are reported in S1 Zip. All associations in the rare variant analyses for the full data set and each subgroup are reported in S2 Zip.

Discussion

In the past 15 years, many pharmacogenetic studies have used novel, genome-wide techniques to identify previously unknown genetic risk factors for drug HS reactions. [2932] Many of the strongest associations have been found for the HLA loci, but other metabolic, transporter, and drug-target genes have also been implicated. [2932,4855] Such associations are particularly important to the rapidly growing field of personalized medicine; in fact, the FDA now recommends prospective genotyping for some of these implicated drugs. [56] Although a genetic basis for sulfonamide HS has long been suspected, little work has been done in immunocompetent patients and none at the genomic level. [9] Therefore, the aim of this study was to identify possible genetic risk factors for sulfonamide HS in the general population using a GWAS approach. Unfortunately, the present study did not demonstrate any genetic associations for sulfonamide HS. These negative findings could result from insufficient study sensitivity, inadequate clinical phenotyping, or represent a true lack of high impact genetic effects for this drug HS syndrome in this population.

The power of a GWAS lies in its ability to simultaneous assess millions of variants across tens of thousands of genes. This provides superior sensitivity over a traditional candidate gene approach in which only one or few genes are investigated. Candidate studies also require an index of suspicion on the part of the investigators that a gene may be mechanistically important to the trait of interest, whereas a genome-wide approach is free from such investigator bias. However, the ability of any study to detect significant differences between groups relies on its sample size and, because a GWAS involves millions of comparisons, large sample sizes are usually required to maintain adequate study sensitivity. The present study was powered to detect a 3-fold increase in relative risk. Pharmacogenetic studies of drug outcomes are typically powered to detect only variants with large genetic effects, because polymorphisms with small effects are unlikely to be adequately predictive to impact clinical decision making. [57] Although the present study was sufficiently powered to detect a relative risk of ≥ 3.0, genotyping more HS and TOL patients would have increased sensitivity and may have allowed us to detect significant risk of lesser magnitude.

Because it is likely that not all of the SNPs are independent, our significance cutoff value of 6.03 x 10−9 may have been overly conservative. If the significance cutoff value were based only on the genotyped SNPs, the p-value would increase from 6.03 x 10−9 to 2.21 x 10−7. Since the smallest p-value found in this study was 4.98 x 10−6, the study conclusions remain the same whether the genotyped SNPs alone, or the genotyped and the imputed SNPs, are used in calculating the significance cutoff value.

The rare variant analysis used in this study accounted for the genetic effects of low MAF SNPs that would have been missed in a traditional GWAS approach. However, this technique does not detect genetic interactions. Development of methods to detect interactions at the SNP, gene, and pathway level is an ongoing area of bioinformatics research. [5860] Study sensitivity may also have been improved had we used a different modality for genotyping. This GWAS was performed using SNP arrays, which included hundreds of thousands of known variants. Combined with imputation, several million genetic markers were used in this study. However, next-generation, whole-genome sequencing exponentially increases the number of possible SNPs included and also assesses non-SNP variants (e.g. indels, inversions, transpositions). Whole-genome sequencing would also allow direct identification of the causative variant, rather than relying on linkage disequilibrium to identify an area of the genome, which then must be re-sequenced.

Our negative findings emphasize the difficulties in studying patients with a diagnosis of drug hypersensitivity as manifested by exanthematous rash, which can have other etiologies that could be misdiagnosed as a drug reaction. We attempted to minimize this by direct medical record review and use of a validated adverse drug reaction scale. Despite these measures, errors and misinterpretation of medical records can still occur. For example, upon secondary review of the cases, it was discovered that one TOL subject failed to complete the minimum 10-day course of TMP/SMX due to GI upset. Given that our analyses did not reveal any significant findings, this oversight is unlikely to affect our results, but it does highlight the difficulties in reviewing medical record information in retrospective studies. Furthermore, drug rechallenge remains the gold standard for confirming adverse drug reactions, [36,61] and this was not performed in most patients in our HS group for ethical reasons. In addition, only banked DNA samples were available for most patients, so we were unable to include other potential measures of causality, such as the lymphocyte transformation test or the in vitro cytotoxicity assay. [16,62,63] These additional biomarkers may have refined our HS population and possibly provided a more clinically uniform phenotype to study.

Because of concerns about phenotype, we also analyzed a subgroup of patients with fever and rash (FEV-RASH) and a subgroup of patients with characteristics of the DRESS (drug reaction with eosinophilia and systemic signs) syndrome (FEV-RASH-EOS). [64] A single gene, COL12A1, was significant for the FEV-RASH group in the rare variant analysis. This gene encodes for type XII collagen, a regulator protein mediating the interaction between type I collagen and the extracellular matrix. [65] Although COL12A1 has been associated with certain familial myopathies, [66] it is not known to be an immunoreactive protein, and there does not seem to be a logical role for a collagen protein in our current understanding of sulfonamide HS pathogenesis. Thus, this modest association was likely due to a Type I error, particularly since the subgroup analyses were underpowered and at risk for sample size bias.

The lack of identifiable genetic associations for sulfonamide HS could be due to the methodical concerns discussed above, but could also represent a true lack of genetic effect for this syndrome. Although a few studies do support a familial inheritance pattern, [17,18] most of the evidence is anecdotal. Recent methods developments have allowed for an estimate of the overall phenotypic variance from GWAS data, but the current study was underpowered to get a reliable estimate. [67] Such estimates are often a challenge in pharmacogenomics studies. [68] With a potentially heritable trait, a genome-wide linkage study among members of an affected family appears an attractive option to minimize enrollment and maximize results. However, it is unusual for all family members to have been treated with the same drug. Furthermore, ethical considerations preclude prospective phenotyping by drug challenge in persons who may have a hereditary predisposition to sulfonamide HS.

Epigenetic and environmental factors may also play a role in the development of sulfonamide HS reactions. For example, patients infected with HIV have an incidence of sulfonamide HS of 20–57% compared to 3% in the general population. [6972] This significant overrepresentation, along with the fact that risk for a reaction increases with declining immune status in AIDS patients, [70,7376] implies that HIV is an acquired risk factor for sulfonamide HS. Several cases of TMP-SMX-induced HS reactions have also been reported during recrudescence of human herpesvirus-6 (HHV6) infection. [7779] Although a cause and effect relationship has not been established, it has been suggested that HHV6 reactivation induces pro-inflammatory cytokines, which might predispose to a cell-mediated immune response, constituting another possible acquired risk factor for sulfonamide HS even in the immunocompetent population. [78,79]

In conclusion, we did not identify any convincing genetic associations for sulfonamide HS, manifested primarily as delayed onset of exanthematous drug eruption with or without fever and eosinophilia, in immunocompetent patients in our Caucasian-American population. These negative findings highlight the need for careful population phenotyping. Clinical use of assays for drug specific T cells or in vitro cytotoxicity may be useful, but new biomarkers for confirming HS reactions are also needed. Large, prospective, multicenter studies would allow for real-time evaluation of phenotype rather than relying on medical record adjudication. Additional studies are currently underway in our laboratory to further investigate the immunopathogenesis of sulfonamide HS.

Supporting Information

S1 Fig. Q-Q plot for common variant analysis of full set of sulfonamide hypersensitive (HS) vs. drug tolerant (TOL) patients (λ = 1.03).

https://doi.org/10.1371/journal.pone.0156000.s001

(TIF)

S2 Fig. Q-Q plot for common variant analysis of HS subgroup with rash and fever (FEV-RASH) vs. all TOL patients (λ = 0.87).

https://doi.org/10.1371/journal.pone.0156000.s002

(TIF)

S3 Fig. Q-Q plot for common variant analysis of HS subgroup with rash, fever, and eosinophilia (FEV-RASH-EOS) vs. all TOL patients (λ = 0.68).

https://doi.org/10.1371/journal.pone.0156000.s003

(TIF)

S1 Zip. Text files including the top 100,000 SNP associations for the common variant analyses for the full data set as well as both subgroups (FEV-RASH and FEV-RASH-EOS).

https://doi.org/10.1371/journal.pone.0156000.s004

(ZIP)

S2 Zip. Text files including all associations for the rare variant analyses for the full data set as well as both subgroups (FEV-RASH and FEV-RASH-EOS).

https://doi.org/10.1371/journal.pone.0156000.s005

(ZIP)

Acknowledgments

The authors would like to thank Terrie Kitchner, Rob Strenn, and the Marshfield Clinic Research Foundation for data abstraction, as well as the University of Wisconsin Biotechnology Center Gene Expression Center for providing the Illumina iScan genotyping services. The authors would also like to thank Skylar Marvel from the North Carolina State University Bioinformatics Research Center for his assistance with data analyses.

Author Contributions

Conceived and designed the experiments: LT AM-R. Performed the experiments: SY AD AM-R. Analyzed the data: JR AD AM-R LT. Contributed reagents/materials/analysis tools: AM-R AD SY. Wrote the paper: JR LT AM-R AD.

References

  1. 1. Avgeri SG, Matthaiou DK, Dimopoulos G, Grammatikos AP, Falagas ME. Therapeutic options for Burkholderia cepacia infections beyond co-trimoxazole: a systematic review of the clinical evidence. Int J Antimicrob Ag. 2009;33(5):394–404.
  2. 2. Drekonja DM, Johnson JR. Urinary tract infections. Prim Care. 2008;35(2):345–367. pmid:18486719
  3. 3. Korbila IP, Manta KG, Siempos II, Dimopoulos G, Falagas ME. Penicillins vs trimethoprim-based regimens for acute bacterial exacerbations of chronic bronchitis: meta-analysis of randomized controlled trials. Can Fam Physician. 2009;55(1):60–67. pmid:19155372
  4. 4. Drekonja DM, Traynor LM, DeCarolis DD, Crossley KB, Johnson JR. Treatment of non-life-threatening methicillin-resistant Staphylococcus aureus infections with alternative antimicrobial agents: a 2-year retrospective review. Diagn Microbiol Infect Dis. 2009;63(2):201–207. pmid:19026510
  5. 5. Gupta D, Zachariah A, Roppelt H, Patel AM, Gruber BL. Prophylactic antibiotic usage for Pneumocystis jirovecii pneumonia in patients with systemic lupus erythematosus on cyclophosphamide: a survey of US rheumatologists and the review of literature. J Clin Rheumatol. 2008;14(5):267–272. pmid:18679133
  6. 6. Rodriguez M, Sifri CD, Fishman JA. Failure of low-dose atovaquone prophylaxis against Pneumocystis jiroveci infection in transplant recipients. Clin Infect Dis. 2004;38(8):e76–78. pmid:15095235
  7. 7. Seger RA. Modern management of chronic granulomatous disease. Br J Haematol. 2008;140(3):255–266. pmid:18217895
  8. 8. Thomas M, Rupali P, Woodhouse A, Ellis-Pegler R. Good outcome with trimethoprim 10 mg/kg/day-sulfamethoxazole 50 mg/kg/day for Pneumocystis jirovecii pneumonia in HIV infected patients. Scand J Infect Dis. 2009;41(11–12):862–868. pmid:19922070
  9. 9. Cribb AE, Lee BL, Trepanier LA, Spielberg SP. Adverse reactions to sulphonamide and sulphonamide-trimethoprim antimicrobials: clinical syndromes and pathogenesis. Adverse Drug React Toxicol Rev. 1996;15(1):9–50. pmid:8920632
  10. 10. Morelli JG, Tay YK, Rogers M, Halbert A, Krafchik B, Weston WL. Fixed drug eruptions in children. J Pediatr. 1999;134(3):365–367. pmid:10064679
  11. 11. Patel RM, Marfatia YS. Clinical study of cutaneous drug eruptions in 200 patients. Indian J Dermatol Venereol Leprol. 2008;74(4):430.
  12. 12. Ozkaya-Bayazit E, Akar U. Fixed drug eruption induced by trimethoprim-sulfamethoxazole: evidence for a link to HLA-A30 B13 Cw6 haplotype. J Am Acad Dermatol. 2001;45(5):712–717. pmid:11606921
  13. 13. Ortiz JE, Horn MS, Peterson HD. Toxic epidermal necrolysis-case report and review of the literature. Ann Plast Surg. 1982;9(3):249–253. pmid:7137822
  14. 14. Saka B, Kombaté K, Mouhari-Toure A, Akakpo S, Tchangaï-Walla K, Pitché P. [Stevens-Johnson syndrome and toxic epidermal necrolysis in a teaching hospital in Lomé, Togo: retrospective study of 89 cases]. Med Trop (Mars). 2010;70(3):255–258.
  15. 15. Wolkenstein P, Charue D, Laurent P, Revuz J, Roujeau JC, Bagot M. Metabolic predisposition to cutaneous adverse drug reactions. Role in toxic epidermal necrolysis caused by sulfonamides and anticonvulsants. Arch Dermatol. 1995;131(5):544–551. pmid:7741540
  16. 16. Neuman MG, Malkiewicz IM, Shear NH. A novel lymphocyte toxicity assay to assess drug hypersensitivity syndromes. Clin Biochem. 2000;33(7):517–524. pmid:11124336
  17. 17. Shear NH, Spielberg SP, Grant DM, Tang BK, Kalow W. Differences in Metabolism of Sulfonamides Predisposing to Idiosyncratic Toxicity. Ann Intern Med. 1986;105(2):179–184. pmid:3729200
  18. 18. Mohanasundaram J, Rajaram S, Mohanasundaram S. Familial tendency in hypersensitivity reactions to co-trimoxazole. J Indian Med Assoc. 1998;96(1):21–22. pmid:9601187
  19. 19. Deloménie C, Grant DM, Mathelier-Fusade P, Jacomet C, Leynadier F, Jacqz-Aigrain E, et al. N-acetylation genotype and risk of severe reactions to sulphonamides in AIDS patients. Br J Clin Pharmacol. 1994;38(6):581–582. pmid:7888297
  20. 20. Pirmohamed M, Alfirevic A, Vilar J, Stalford A, Wilkins EG, Sim E, et al. Association analysis of drug metabolizing enzyme gene polymorphisms in HIV-positive patients with co-trimoxazole hypersensitivity. Pharmacogenetics. 2000;10(8):705–713. pmid:11186133
  21. 21. Wolkenstein P, Loriot M-A, Flahault A, Cadilhac M, Caumes E, Eliaszewicz M, et al. Association analysis of drug metabolizing enzyme gene polymorphisms in AIDS patients with cutaneous reactions to sulfonamides. J Invest Dermatol. 2005;125(5):1080–1082. pmid:16297214
  22. 22. Alfirevic A, Vilar FJ, Alsbou M, Jawaid A, Thomson W, Ollier WE, et al. TNF, LTA, HSPA1L and HLA-DR gene polymorphisms in HIV-positive patients with hypersensitivity to cotrimoxazole. Pharmacogenomics. 2009;10(4):531–540. pmid:19374512
  23. 23. Wang D, Curtis A, Papp AC, Koletar SL, Para MF. Polymorphism in glutamate cysteine ligase catalytic subunit (GCLC) is associated with sulfamethoxazole-induced hypersensitivity in HIV/AIDS patients. BMC Med Genomics. 2012;5(1):32.
  24. 24. Zielińska E, Niewiarowski W, Bodalski J. The arylamine N-acetyltransferase (NAT2) polymorphism and the risk of adverse reactions to co-trimoxazole in children. E J Clin Pharmacol. 1998;54(9–10):779–785.
  25. 25. Wolkenstein P, Carrière V, Charue D, Bastuji-Garin S, Revuz J, Roujeau JC, et al. A slow acetylator genotype is a risk factor for sulphonamide-induced toxic epidermal necrolysis and Stevens-Johnson syndrome. Pharmacogenetics. 1995;5(4):255–258. pmid:8528274
  26. 26. Sacco JC, Abouraya M, Motsinger-Reif A, Yale SH, McCarty CA, Trepanier LA. Evaluation of polymorphisms in the sulfonamide detoxification genes NAT2, CYB5A, and CYB5R3 in patients with sulfonamide hypersensitivity. Pharmacogenet Genomics. 2012;22(10):733–740. pmid:22850190
  27. 27. Gill HJ, Tjia JF, Kitteringham NR, Pirmohamed M, Back DJ, Park BK. The effect of genetic polymorphisms in CYP2C9 on sulphamethoxazole N-hydroxylation. Pharmacogenetics. 1999;9(1):43–53. pmid:10208642
  28. 28. Ateş NA, Tursen U, Tamer L, Kanik A, Derici E. Glutathione S-transferase polymorphisms in patients with drug eruption. Arch Dermatol Res. 2004;295:429–433. pmid:14740231
  29. 29. Daly AK, Donaldson PT, Bhatnagar P, Shen Y, Pe'er I, Floratos A, et al. HLA-B*5701 genotype is a major determinant of drug-induced liver injury due to flucloxacillin. Nat Genet. 2009;41(7):816–819. pmid:19483685
  30. 30. Del-Aguila JL, Beitelshees AL, Cooper-Dehoff RM, Chapman AB, Gums JG, Bailey K, et al. Genome-wide association analyses suggest NELL1 influences adverse metabolic response to HCTZ in African Americans. Pharmacogenomics J. 2014;14(1):35–40. pmid:23400010
  31. 31. Lucena MI, Molokhia M, Shen Y, Urban TJ, Aithal GP, Andrade RJ, et al. Susceptibility to amoxicillin-clavulanate-induced liver injury is influenced by multiple HLA class I and II alleles. Gastroenterology. 2011;141(1):338–347. pmid:21570397
  32. 32. SEARCH Collaborative Group, Link E, Parish S, Armitage J, Bowman L, Heath S, et al. SLCO1B1 variants and statin-induced myopathy—a genomewide study. N Engl J Med. 2008;359(8):789–799. pmid:18650507
  33. 33. McCarty GA. Electronic medical record embedded functions facilitate attending-trainee care oversight and education. AMIA Annu Symp Proc. 2005:1049. pmid:16779336
  34. 34. Harle DG, Baldo BA, Wells JV. Drugs as allergens: detection and combining site specificities of IgE antibodies to sulfamethoxazole. Mol Immunol. 1988;25(12):1347–1354. pmid:3237218
  35. 35. Gruchalla RS, Sullivan TJ. Detection of human IgE to sulfamethoxazole by skin testing with sulfamethoxazoyl-poly-L-tyrosine. J Allergy Clin Immun. 1991;88(5):784–792. pmid:1955637
  36. 36. Naranjo CA, Busto U, Sellers EM, Sandor P, Ruiz I, Roberts EA, et al. A method for estimating the probability of adverse drug reactions. Clin Pharmacol Ther. 1981;30(2):239–245. pmid:7249508
  37. 37. Laurie CC, Doheny KF, Mirel DB, Pugh EW, Bierut LJ, Bhangale T, et al. Quality control and quality assurance in genotypic data for genome-wide association studies. Genet Epidemiol. 2010;34(6):591–602. pmid:20718045
  38. 38. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–575. pmid:17701901
  39. 39. Putnam PP, Zhang G, Wilsey PA. A comparison study of succinct data structures for use in GWAS. BMC Bioinformatics. 2013;14(1):369.
  40. 40. Azzopardi D, Dallosso AR, Eliason K, Hendrickson BC, Jones N, Rawstorne E, et al. Multiple rare nonsynonymous variants in the adenomatous polyposis coli gene predispose to colorectal adenomas. Cancer Res. 2008;68(2):358–363. pmid:18199528
  41. 41. Cohen JC, Kiss RS, Pertsemlidis A, Marcel YL, McPherson R, Hobbs HH. Multiple rare alleles contribute to low plasma levels of HDL cholesterol. Science. 2004;305(5685):869–872. pmid:15297675
  42. 42. Walsh T, McClellan JM, McCarthy SE, Addington AM, Pierce SB, Cooper GM, et al. Rare structural variants disrupt multiple genes in neurodevelopmental pathways in schizophrenia. Science. 2008;320(5875):539–543. pmid:18369103
  43. 43. Li B, Leal SM. Methods for detecting associations with rare variants for common diseases: application to analysis of sequence data. Am J Hum Genet. 2008;83(3):311–321. pmid:18691683
  44. 44. Lee S, Abecasis GR, Boehnke M, Lin X. Rare-variant association analysis: study designs and statistical tests. Am J Hum Genet. 2014;95(1):5–23. pmid:24995866
  45. 45. Wu MC, Lee S, Cai T, Li Y, Boehnke M, Lin X. Rare-variant association testing for sequencing data with the sequence kernel association test. Am J Hum Genet. 2011;89(1):82–93. pmid:21737059
  46. 46. Lee S, Wu MC, Lin X. Optimal tests for rare variant effects in sequencing association studies. Biostatistics. 2012;13(4): 762–775. pmid:22699862
  47. 47. Lancaster HO. The combination of probabilities: an application of orthonormal functions. Aust J Stat. 1961;3(1):20–33.
  48. 48. Pavlos R, Mallal S, Phillips E. HLA and pharmacogenetics of drug hypersensitivity. Pharmacogenomics. 2012;13(11):1285–1306. pmid:22920398
  49. 49. Quiralte J, Sánchez-García F, Torres MJ, Blanco C, Castillo R, Ortega N, et al. Association of HLA-DR11 with the anaphylactoid reaction caused by nonsteroidal anti-inflammatory drugs. J Allergy Clin Immunol. 1999;103(4):685–689. pmid:10200020
  50. 50. Hung SI, Chung WH, Liou LB, Chu CC, Lin M, Huang HP, et al. HLA-B*5801 allele as a genetic marker for severe cutaneous adverse reactions caused by allopurinol. Proc Natl Acad Sci. 2005;102(11):4134–4139. pmid:15743917
  51. 51. Mallal S, Nolan D, Witt C, Masel G, Martin AM, Moore C, et al. Association between presence of HLA-B*5701, HLA-DR7, and HLA-DQ3 and hypersensitivity to HIV-1 reverse-transcriptase inhibitor abacavir. Lancet. 2002;359(9308):727–732. pmid:11888582
  52. 52. Chung WH, Hung SI, Hong HS, Hsih MS, Yang LC, Ho HC, et al. Medical genetics: a marker for Stevens-Johnson syndrome. Nature. 2004;428(6982):486–486. pmid:15057820
  53. 53. Rizzo R, Spaggiari F, Indelli M, Lelli G, Baricordi OR, Rimessi P, et al. Association of CYP1B1 with hypersensitivity induced by taxane therapy in breast cancer patients. Breast Cancer Res Treat. 2010;124(2):593–598. pmid:20632082
  54. 54. Westphal GA, Schnuch A, Schulz TG, Reich K, Aberer W, Brasch J, et al. Homozygous gene deletions of the glutathione S-transferases M1 and T1 are associated with thimerosal sensitization. Int Arch Occup Environ Health. 2000;73(6):384–388. pmid:11007341
  55. 55. Takeuchi F, McGinnis R, Bourgeois S, Barnes C, Eriksson N, Soranzo N, et al. A genome-wide association study confirms VKORC1, CYP2C9, and CYP4F2 as principal genetic determinants of warfarin dose. PLoS Genet. 2009;5(3).
  56. 56. Ziagen [package insert]. Research Triangle Park, NC: GlaxoSmithKline; 2008.
  57. 57. Nelson MR, Bacanu S-A, Mosteller M, Li L, Bowman CE, Roses AD, et al. Genome-wide approaches to identify pharmacogenetic contributions to adverse drug reactions. Pharmacogenomics J. 2009;9(1):23–33. pmid:18301416
  58. 58. Li P, Guo M, Wang C, Liu X, Zou Q. An overview of SNP interactions in genome-wide association studies. Brief Funct Genomics. 2015;14(2):143–155. pmid:25241224
  59. 59. Xu J, Yuan Z, Ji J, Zhang X, Li H, Wu X, et al. A powerful score-based test statistic for detecting gene-gene co-association. BMC Genet. 2016;17:31. pmid:26822525
  60. 60. Mooney MA, Nigg JT, McWeeney SK, Wilmot B. Functional and genomic context in pathway analysis of GWAS data. Trends Genet. 2014;30(9):390–400. pmid:25154796
  61. 61. Bousquet PJ, Gaeta F, Bousquet-Rouanet L, Lefrant JY, Demoly P, Romano A. Provocation tests in diagnosing drug hypersensitivity. Curr Pharm Des. 2008;14:2792–2802. pmid:18991698
  62. 62. Rieder MJ, Uetrecht J, Shear NH, Cannon M, Miller M, Spielberg SP. Diagnosis of sulfonamide hypersensitivity reactions by in vitro rechallenge with hydroxylamine metabolites. Ann Intern Med. 1989;110(4):286–289. pmid:2913915
  63. 63. Reilly TP, Bellevue FH, Woster PM, Svensson CK. Comparison of the in vitro cytotoxicity of hydroxylamine metabolites of sulfamethoxazole and dapsone. Biochem Pharmacol. 1998;55(6):803–810. pmid:9586952
  64. 64. Kardaun SH, Sekula P, Valeyrie-Allanore L, Liss Y, Chu CY, Creamer D, et al. Drug reaction with eosinophilia and systemic symptoms (DRESS): an original multisystem adverse drug reaction. Results from the prospective RegiSCAR study. Br J Dermatol. 2013;169(5):1071–1080. pmid:23855313
  65. 65. Koch M, Bohrmann B, Matthison M. Large and small splice variants of collagen XII: differential expression and ligand binding. J Cell Biol. 1995;130(4):1005–1014. pmid:7642694
  66. 66. Hicks D, Farsani GT, Laval S, Collins J, Sarkozy A, Martoni E, et al. Mutations in the collagen XII gene define a new form of extracellular matrix-related myopathy. Human Molecular Genetics. 2014;23(9):2353–2363. pmid:24334769
  67. 67. Yang J, Less SH, Goddard ME, Visscher PM. GCTA: a tool for Genome-wide complex trait analysis. Am J Hum Genet. 2011;88(1):76–82. pmid:21167468
  68. 68. Motsinger-Reif AA, Jorgenson E, Relling MV, Kroetz DL, Weinshilboum R, Cox NJ, et al. Genome-wide association studies in pharmacogenomics: successes and lessons. Pharmacogenet Genomics. 2013;23(8):383–94. pmid:20639796
  69. 69. Gordin FM, Simon GL, Wofsy CB, Mills J. Adverse reactions to trimethoprim-sulfamethoxazole in patients with the acquired immunodeficiency syndrome. Ann Intern Med. 1984;100(4):495–499. pmid:6230976
  70. 70. Hennessy S, Strom BL, Berlin JA, Brennan PJ. Predicting cutaneous hypersensitivity reactions to cotrimoxazole in HIV-infected individuals receiving primary Pneumocystis carinii pneumonia prophylaxis. J Gen Intern Med. 1995;10(7):380–386. pmid:7472686
  71. 71. Walmsley SL, Khorasheh S, Singer J, Djurdjev O. A randomized trial of N-acetylcysteine for prevention of trimethoprim-sulfamethoxazole hypersensitivity reactions in Pneumocystis carinii pneumonia prophylaxis (CTN 057). J Acquir Immune Defic Syndr Hum Retrovirol. 1998;19(5):498–505. pmid:9859964
  72. 72. Medina I, Mills J, Leoung G, Hopewell PC, Lee B, Modin G, et al. Oral therapy for Pneumocystis carinii pneumonia in the acquired immunodeficiency syndrome. A controlled trial of trimethoprim-sulfamethoxazole versus trimethoprim-dapsone. N Engl J Med. 1990;323(12):776–782. pmid:2392131
  73. 73. Carr A, Swanson C, Penny R, Cooper DA. Clinical and laboratory markers of hypersensitivity to trimethoprim-sulfamethoxazole in patients with Pneumocystis carinii pneumonia and AIDS. J Infect Dis. 1993;167(1):180–185. pmid:8380290
  74. 74. Eliaszewicz M, Flahault A, Roujeau JC, Fillet AM. Prospective evaluation of risk factors of cutaneous drug reactions to sulfonamides in patients with AIDS. J Am Acad Dermatol. 2002;47:40–46. pmid:12077579
  75. 75. Rabaud C, Charreau I, Izard S, Raffi F, Meiffredy V, Leport C, et al. Adverse reactions to cotrimoxazole in HIV-infected patients: predictive factors and subsequent HIV disease progression. Scand J Infect Dis. 2001;33(10):759–764. pmid:11728044
  76. 76. Ryan C, Madalon M, Wortham DW, Graziano FM. Sulfa hypersensitivity in patients with HIV infection: onset, treatment, critical review of the literature. Wisc Med J. 1998;97(5):23–27.
  77. 77. Hubiche T, Milpied B, Cazeau C, Taïeb A, Léauté-Labrèze C. Association of immunologically confirmed delayed drug reaction and human herpesvirus 6 viremia in a pediatric case of drug-induced hypersensitivity syndrome. Dermatology. 2011;222(2):140–141. pmid:21389685
  78. 78. Morimoto T, Sato T, Matsuoka A. Trimethoprim-sulfamethoxazole-induced hypersensitivity syndrome associated with reactivation of human herpesvirus-6. Int Med. 2006;45(2)101–105.
  79. 79. Neuman MG, McKinney KK, Nanau RM, Kong V, Malkiewicz I, Mazulli T, et al. Drug-induced severe adverse reaction enhanced by human herpes virus-6 reactivation. Transl Res. 2013;161(5):430–440. pmid:23333110