FCERI and Histamine Metabolism Gene Variability in Selective Responders to NSAIDS

The high-affinity IgE receptor (Fcε RI) is a heterotetramer of three subunits: Fcε RIα, Fcε RIβ, and Fcε RIγ (αβγ2) encoded by three genes designated as FCER1A, FCER1B (MS4A2), and FCER1G, respectively. Recent evidence points to FCERI gene variability as a relevant factor in the risk of developing allergic diseases. Because Fcε RI plays a key role in the events downstream of the triggering factors in immunological response, we hypothesized that FCERI gene variants might be related with the risk of, or with the clinical response to, selective (IgE mediated) non-steroidal anti-inflammatory (NSAID) hypersensitivity. From a cohort of 314 patients suffering from selective hypersensitivity to metamizole, ibuprofen, diclofenac, paracetamol, acetylsalicylic acid (ASA), propifenazone, naproxen, ketoprofen, dexketoprofen, etofenamate, aceclofenac, etoricoxib, dexibuprofen, indomethacin, oxyphenylbutazone, or piroxicam, and 585 unrelated healthy controls that tolerated these NSAIDs, we analyzed the putative effects of the FCERI SNPs FCER1A rs2494262, rs2427837, and rs2251746; FCER1B rs1441586, rs569108, and rs512555; FCER1G rs11587213, rs2070901, and rs11421. Furthermore, in order to identify additional genetic markers which might be associated with the risk of developing selective NSAID hypersensitivity, or which may modify the putative association of FCERI gene variations with risk, we analyzed polymorphisms known to affect histamine synthesis or metabolism, such as rs17740607, rs2073440, rs1801105, rs2052129, rs10156191, rs1049742, and rs1049793 in the HDC, HNMT, and DAO genes. No major genetic associations with risk or with clinical presentation, and no gene-gene interactions, or gene-phenotype interactions (including age, gender, IgE concentration, antecedents of atopy, culprit drug, or clinical presentation) were identified in patients. However, logistic regression analyses indicated that the presence of antecedents of atopy and the DAO SNP rs2052129 (GG) were strongly related (P < 0.001 and P = 0.005, respectively) with selective hypersensitivity to ibuprofen. With regard to patients with selective hypersensitivity to ASA, men were more prone to develop such a reaction than women (P = 0.011), and the detrimental DAO SNP rs10156191 in homozygosity increased the risk of developing such hypersensitivity (P = 0.039).

Introduction

Type B drug-induced hypersensitivity reactions (DHR) occur only in susceptible individuals with a frequency of 5–10% of all adverse drug reactions (Khan and Solensky, 2010). In general, these reactions are severe and occasionally may be life-threatening. In recent years evidence has accumulated to support the notion that DHR may be caused by various mechanisms and is the result of a complex multifactorial and multigenic process (Pirmohamed, 2006). Regarding culprit drugs, non-steroidal anti-inflammatory drugs (NSAIDs) are among the most frequent causes of DHR together with antibiotics (Cornejo-Garcia et al., 2009; Doña et al., 2011, 2012; Kowalski et al., 2011). Hypersensitivity reactions to a single NSAID (selective reactions) are the result of an immunological mechanism, which is either IgE-mediated in acute reactions or T cell-mediated in delayed reactions. Clinical manifestations include cutaneous reactions, respiratory reactions and anaphylaxis, which may appear with different onset times after drug intake. Histamine release from mast cells after IgE receptor activation plays a relevant role in allergic inflammation and in the development of clinical symptoms (Mita et al., 2001; Kowalski et al., 2011).

A priori, genetic variations affecting both components (expression and/or function) of this signaling pathway, including high affinity IgE receptor and histamine metabolizing enzymes, could modify the risk of developing hypersensitivity to NSAIDs, and/or the presentation of clinical manifestations. In fact, recent studies have demonstrated an association between several allergic diseases including drug hypersensitivity and polymorphisms in these genes (García-Martín et al., 2006, 2007a; Kennedy et al., 2008; Gervasini et al., 2010; Maintz et al., 2011).

The high-affinity IgE receptor (Fcε RI) is a heterotetramer of three subunits: Fcε RIα, Fcε RIβ, and Fcε RIγ (αβγ₂) encoded by three genes designated as FCER1A, FCER1B (MS4A2), and FCER1G, respectively. The genes coding for all FcεRI subunits are polymorphic and some of these polymorphisms have been associated with several atopic diseases (MacGlashan et al., 1998, 1999; Saini et al., 1999; Weidinger et al., 2008; Zhang et al., 2010; Li et al., 2014). The FCER1A gene is located on chromosome 1q23 (1:159259504-159278014; 1000 Genomes, 2015), and it encodes a protein with two extracellular IgE-like domains with ligand-binding sites (Potaczek and Kabesch, 2012). The FCER1B gene is located on chromosome 11q12-13 (11:59855734–59863444; 1000 Genomes, 2015), and it encodes a 244-amino acid protein with a non-canonic intracellular ITAM domain. Functionally, it is a signal-augmenting subunit. The FCER1G gene is located on chromosome 1q23 (1:161185024–161190489; 1000 Genomes, 2015) and encodes an 86-amino acid protein. Functionally, it is a signal-transducing subunit and plays an essential role in the induction of mast cell degranulation and survival (Manikandan et al., 2012).

With regard to histamine homeostasis, histamine biosynthesis is catalyzed by the enzyme L-histidine decarboxylase (HDC, E.C. 4.1.1.22). Two enzymes participate in the degradation of histamine: Histamine N-methyltransferase (HNMT, E.C. 2.1.1.8) and diamine oxidase (DAO; E.C. 1.4.3.6; García-Martín et al., 2009). There is high individual variability in histamine metabolism which is, at least in part, genetically determined, although other factors, such as gender (García-Martín et al., 2007b) play a major role in the variability of histamine metabolism. An association between polymorphisms of histamine metabolizing enzymes and the clinical presentation of allergic diseases has been demonstrated (García-Martín et al., 2006, 2007a; Kennedy et al., 2008; Gervasini et al., 2010; Maintz et al., 2011).

The HDC gene is located in chromosome 15q21-q22 and spans around 24 kilobases. The HNMT gene is located in chromosome 2q22.1. The human DAO gene spans ~10 kb and is located in chromosome 7q34-q36. A promoter SNP and three common non-synonymous SNPs have been identified in Caucasian individuals, and the functional effects of these SNPs in enzyme activity have been studied in detail (Ayuso et al., 2007; García-Martín et al., 2007b).

In order to investigate the role of genetic factors in the risk of developing selective NSAID hypersensitivity, both isolated and combined, given that the clinical outcome depends on the interplay of IgE response and the consequent release of mediators, in this study we analyzed functional polymorphisms in high affinity IgE receptors and genes related to histamine metabolism in a large group of well-phenotyped patients suffering from single-NSAID-induced hypersensitivity reactions. Our aim was to elucidate the putative association of these gene polymorphisms with several clinical phenotype parameters, including gender, age, antecedents of atopy, culprit drug, IgE levels and clinical presentation, and to analyze putative gene-gene interactions. The final goal was to identify biomarkers which alone or combined with pharmacogenomics or metabolic biomarkers (Agúndez et al., 2009a,b, 2015; Cornejo-García et al., 2012; Martinez et al., 2014; García-Martín et al., 2015a; Ariza et al., 2016) might be useful in prevention, diagnosis and/or proper management of these patients.

Patients and Methods

Study Population

We studied a cohort of 899 individuals, consisting of 314 unrelated Spanish patients with single-NSAID-induced hypersensitivity reactions and 585 unrelated healthy controls. Written consent for participation was obtained for all participants. All the patients who were invited to participate in the study agreed to do so. Of the patients, 145 were recruited from the Allergy service, Hospital Infanta Leonor (Madrid, Spain), 100 were recruited from the Allergy service, Regional Hospital (Málaga, Spain) and 69 were recruited from the Allergy Department, Infanta Cristina Hospital (Badajoz, Spain). Diagnosis was carried out as previously described (Doña et al., 2011). After confirming good tolerance to a full therapeutic dose of ASA, incremental doses of the culprit drug were given until the therapeutic dose for the analgesic/NSAIDs implicated in the study was achieved. In the case of selective responses to ASA, confirmation of the selective response was made by challenge with indomethacin as reported elsewhere (Doña et al., 2011).

Participants' characteristics are summarized in Table 1. The NSAIDs responsible for the reaction are summarized in Table 2. Clinical presentation distributed according to gender and culprit drug are shown in Table 3.

TABLE 1

Table 1. Characteristics of the study group.

TABLE 2

Table 2. Culprit drug for selective NSAID-induced hypersensitivity.

TABLE 3

Table 3. Clinical presentation of selective NSAID-induced hypersensitivity.

The healthy controls were recruited from staff and medical students of the Hospitals and the Universities participating in the study, and were ethnically matched with patients (all were unrelated Spanish individuals). All control subjects tolerated the NSAIDs most frequently involved in selective hypersensitivity, as shown in Table 2. Specifically, all control individuals previously received metamizole, ibuprofen, diclofenac, paracetamol, and ASA, and experienced no hypersensitivity or other ADRs. Individuals with the above-mentioned characteristics were asked to participate, and 97% of them agreed. A medical history was obtained and an examination was performed for each participant to exclude pre-existing disorders. Individuals with familial (up to second-degree relatives) or personal antecedents of allergic, atopic, or autoimmune diseases were excluded from the control group to avoid confounders. The protocol for this study was in accordance with the Declaration of Helsinki and its subsequent revisions and was approved by the respective Ethics Committees of the participating Hospitals.

Genotype Analysis

Genomic DNA was obtained from peripheral leukocytes and purified in accordance with standard procedures. The SNPs analyzed were selected according to allele frequencies (over 0.01) in the study population, and either functional or clinical relevance, in line with published evidence (Preuss et al., 1998; García-Martín et al., 2006, 2007a, 2009; Ayuso et al., 2007; Maintz et al., 2011; Amo et al., 2016), and the public 1000 genomes database release of 17 Nov. 2015. Genotyping was performed by TaqMan assays (Life Technologies, Alcobendas, Madrid, Spain). Details of the TaqMan probes and the allele frequencies in Caucasian individuals are summarized in Table 4. We studied nine FCER1 SNPs, as well as SNPs corresponding to genes involved in histamine synthesis (HDC) or metabolism (HNMT and DAO). All these histamine-related SNPs have demonstrated functional and/or clinical implications (García-Martín et al., 2006, 2007a,b, 2009, 2015b; Ayuso et al., 2007; Gervasini et al., 2010; Agúndez et al., 2012).

TABLE 4

Table 4. SNPs analyzed in this study.

Detection was carried out by means of real-time PCR (qPCR) in an Eppendorf realplex thermocycler using fluorescent probes. The amplification conditions were as follows: After a denaturation time of 10 min at 96°C, 45 cycles of 92°C 15 s 60°C 90 s were carried out and fluorescence was measured at the end of each cycle and at endpoint. All samples were determined in triplicate and genotypes were assigned both by means of gene identification software (RealPlex 2.0, Eppendorf), and by analysis of the reference cycle number for each fluorescence curve, calculated using the CalQPlex algorithm (Eppendorf).

Statistical Analyses

Statistical power for each SNP was evaluated with a genetic model to analyze the frequency for the minor allele with an odds ratio (OR) value = 1.5 (α = 0.05) based on the allele frequencies observed in the control group. Table 4 shows the statistical power for each SNP analyzed. For most SNPs the statistical power was very high. In some cases, because of the low minor allele frequency observed, the statistical power was not sufficient to detect an OR = 1.5 but was sufficient to detect an OR = 2.0 with a bilateral power of more than 80% (Table 4), with the single exception of the HDC SNP rs2073440 T/G, whose power was sufficient to detect an OR = 2.1.

SNPStats software (Solé et al., 2006) was used to calculate allele and genotype frequencies, to analyze the Hardy-Weinberg equilibrium, and to determine linkage disequilibrium statistics and haplotype frequency estimation. Multiple comparison adjustment was done by using the False Discover Rate (FDR) correction (http://www.sdmproject.com/utilities/?show=FDR). Analyses of association with a response variable (culprit drugs and reactions) based on logistic regression were performed using SPSS 21.0 for Windows. For these analyses, we also determined the odds ratio and the corresponding 95% confidence intervals. The Hosmer-Lemeshow goodness of fit test for logistic regression was used. For the IgE response variable a multiple linear regression model was calculated. The results were considered as statistically significant when the p-value was less than 0.05.

Results

The percentage of women was slightly higher in the cases, as compared to the control individuals (Table 1), although no statistically significant gender differences between cases and controls were present (Chi-square P = 0.191). Conversely, age was lower in the controls than in the cases. However, age is not a key factor in this study, as all the control individuals were tolerant to NSAIDs and the odds are extremely low that they may eventually develop selective NSAID-induced hypersensitivity.

The most common culprit drugs for selective NSAID-induced hypersensitivity were metamizole, ibuprofen, diclofenac, paracetamol, and ASA (Table 2). Although the frequency for metamizole-induced hypersensitivity was higher in women than in men, OR = 1.42 and 95% confidence interval (CI) 0.84–2.41; P = 0.171, the frequency difference was not statistically significant, and neither were the differences for gender-related frequencies for the rest of the NSAIDs included in Table 2. These frequencies correspond to those previously described by our group (Doña et al., 2011; Blanca-López et al., 2016a,b).

The most frequent clinical presentation was urticaria + angioedema, followed by anaphylaxis, exanthema and mixed pattern. No gender-related differences in clinical presentation were observed (Table 3). Clinical presentation, however, was strongly related with the culprit drug: The clinical presentation urticaria + angioedema was particularly frequent when the culprit drug was ibuprofen (OR = 3.50, 95% CI = 1.99–6.19); P <0.001. In contrast, when the causative drug was metamizole or diclofenac, the most common clinical presentation was anaphylaxis (OR = 2.65, 95% CI = 1.58–4.44; P < 0.001) and (OR = 2.89, 95% CI = 1.34–6.28, P = 0.003), respectively. When the culprit drug was paracetamol or ASA the most frequent presentation was urticaria + angioedema, although for these two drugs the association of the drug with clinical presentation was not statistically significant. These phenotypic features correspond to those previously reported for selective NSAID hypersensitivity patients (Cornejo-Garcia et al., 2009; Doña et al., 2011; García-Martín et al., 2015a).

FCER1 genotyping results are summarized in Table 5. We checked the codominant, dominant, recessive, overdominant and additive models and the best fit was obtained with the recessive model. FCER1 SNPs did not show statistically significant differences when patients and controls were compared in any of the genetic models analyzed. The genotyping results related to histamine synthesis and metabolism genes are summarized in Table 6. Once again, the best fit for histamine metabolism genes was obtained with the recessive model. Statistically significant differences on comparing cases and control subjects were identified for the DAO SNP rs10156191, which caused decreased enzyme activity (Ayuso et al., 2007). The observed difference was related to the frequency of homozygous individuals for the minor allele, which obtained a marginal significance that was not observed when the allele frequency (instead of the genotypes) was analyzed. When correction for multiple comparisons was carried out by using FDR, the P value for the recessive model for the DAO rs10156191 SNP was not significant (corrected P = 0.294), whereas the P-value for carriers of the minor allele for the HDC SNP rs2073440 remained significant (corrected P-value = 0.021).

TABLE 5

Table 5. FCER1 SNPs analyzed in this study.

TABLE 6

Table 6. SNPs related with histamine synthesis and degradation analyzed in this study.

Table 7 shows the statistically significant interaction of the genotypes studied and gender. The association of the FCER1A SNP rs2427837 with the risk of developing NSAID-induced hypersensitivity showed a positive interaction with gender, the association being stronger in women (Table 7). A statistically significant, genotype-gender interaction was observed for two other FCER1 genotypes (rs2251746 in women and rs11587213 in men). When FDR correction for multiple comparisons (both genders and all SNPs) was made, the corrected P-values remained significant: P = 0.045 for the genetic associations with the SNPs rs2427837 and rs2251746 and P = 0.035 for the SNP rs11587213.

TABLE 7

Table 7. Statistically significant gender-related risk associations.

Logistic regression analyses were carried out with separate models for each clinical presentation, as shown in Table 3, by comparing between cases all genotypes, gender, age, IgE concentration, and antecedents of atopy. In addition, we analyzed putative associations with response, stratifying patients into two groups: Single NSAID-induced urticarial/angioedema or anaphylaxis (SNIUAA) and single NSAID-induced delayed reactions (SNIDR). No significant associations were identified. In addition, logistic regression analyses were carried out with separate models for each culprit drug (only drugs with 19 or more cases were included, as shown in Table 2), by comparing between cases all genotypes, gender, age, IgE concentration, and antecedents of atopy.

For ibuprofen, the Hosmer–Lemeshow (HL) goodness of fit test was equal to 0.80 and revealed that age [P = 0.032; OR (95% CI) = 0.91 (0.83–0.99)], the presence of atopy antecedents [P < 0.001; OR = 19.61 (4.13–90.90)] and the DAO SNP rs2052129 (GG); [P = 0.005; OR = 13.25 (2.14–81.84)] were related to the risk of developing hypersensitivity. With regard to metamizole, although linear regression analysis suggested association with the absence of atopy (P = 0.008; OR = 5.80, 95% CI = 1.58–21.28), the HL goodness of fit test was equal to 0.001, and therefore the significance of these findings is limited. For ASA, we identified significant associations with the SNPs rs10156191 (TT) [P = 0.035, OR = 44.59; 95% CI = 1.22–1630.56] and with gender (for men, P = 0.011, OR = 25.64; 95% CI = 2.13–333.33) with an HL goodness of fit test equal to 0.957. No significant associations were observed with paracetamol or diclofenac. IgE levels did not show any association with phenomic or genomic markers. We did not identify any additional significant associations, although we cannot rule out association with other culprit drugs because the subgroup sizes were not sufficiently large to reach statistical significance.

Discussion

NSAID-induced hypersensitivity type B adverse reactions are mediated by immunological and non-immunological mechanisms. Two major clinical phenotypes have been described: Selective NSAID hypersensitivity, which is drug-specific and an IgE-mediated mechanism, and cross-intolerance in which chemically non-related NSAIDs induce the reaction (Kowalski et al., 2013). Because selective hypersensitivity is an IgE-mediated mechanism, we analyzed genetic variations at the high-affinity IgE receptor, which has been shown to be related with allergic disorders (MacGlashan et al., 1998, 1999; Saini et al., 1999; Weidinger et al., 2008; Zhang et al., 2010; Li et al., 2014; Amo et al., 2016). It is to be noted that, despite the large body of published evidence supporting association of FCER1 SNPs with allergic diseases, this is the first study to analyze the putative role of FCER1 SNPs in selective NSAID hypersensitivity. In addition, we previously identified genetic factors related to cross intolerance (Agúndez et al., 2012) and, of these, one non-synonymous DAO gene variation, designated as rs10156191, was overrepresented among cross-intolerant patients, thus providing the basis for a detailed study on the role of genetic variations in histamine metabolism in patients with selective hypersensitivity to NSAIDs.

FCER1 genotypes have been linked to ASA-intolerant asthma. In a study carried out on 126 Korean patients with ASA-intolerant asthma, Palikhe and co-workers analyzed six FCER1 SNPs, five of which were also analyzed in our study (Palikhe et al., 2008a,b). They identified a weak association of the FCER1G rs11587213 SNP with ASA-intolerant asthma, the patients showing increased frequency for the AA genotype. The same study reported an association of the two FCER1A SNPs rs2427827 and rs2251746 with specific IgE levels and an association of the FCER1G SNP 11587213 with both total IgE and specific IgE levels.

Our findings do not support an association of these FCER1A genotypes with the risk of developing selective hypersensitivity to ASA, or any other NSAID included in this study, and neither do they support an association with total IgE levels. Potential discrepancies between our study and that of Palikhe and co-workers may arise from the differences in clinical presentation, differences in the number of patients, and the different ethnic origin of patients and controls. In fact, in the study by Palikhe et al. (2008a,b) the allele frequencies observed in Korean individuals differ considerably from those reported in this study and in the 1000 genomes website for individuals of Caucasian descent (see Table 4).

Additional clinical associations for FCER1A SNPs are the putative association of the SNP rs2298804 with the risk of developing systemic lupus erythematosus in a study carried out in China (Yang et al., 2013), and the association of the SNP rs2298805 with the risk of developing chronic spontaneous urticaria in Chinese individuals (Guo et al., 2015). These SNPs are ethnic-specific as they have only been identified in Oriental individuals, but they do not occur in Caucasian individuals according to the 1000 genomes website. Additional clinical associations for FCER1B SNPs, all related to the SNP rs569108, include increased risk of developing asthma in Chinese individuals (Ramphul et al., 2014; Hua et al., 2016), and atopic allergy in individuals from Philippines (de Guia et al., 2015). A meta-analysis of 24 studies also supports association of the FCER1B rs569108 SNP with asthma, although the risk seems to be restricted to East-Asian individuals (Yang et al., 2014).

Additional clinical associations for FCER1G include a weak effect on food sensitization, which is associated with the interactive effect of the FCER1G rs2070901 SNP with other SNPs in the IL4, FCER1B, and CYP24A1 genes and cord blood 25(OH) D (Liu et al., 2011).

Strengths of this study include a high number of patients with selective NSAID-hypersensitivity (n = 314). Moreover, the clinical phenotypes of these patients, including the proportion of each gender, ages, culprit drugs, and patients' clinical presentations correspond to those described previously among Spaniards, thus indicating that the patient group is representative. The number of patients and controls is sufficiently high to obtain a good statistical power, which is required to obtain conclusive evidence. Limitations in this study include a low number of patients for some subgroups according to the culprit drug (Table 2), the younger age of control individuals as compared to patients, and the low frequency for some of the SNPs analyzed which were, nevertheless, included in the study because of their functional or clinical impact (Table 4).

The results of this study do not support a major association of FCER1 genotypes in the risk of developing selective NSAID-induced hypersensitivity. Similar genotype distributions and allele frequencies were observed among patients and controls, and the genotypes and frequencies correspond to those previously reported among Spanish patients (Agúndez et al., 2012; Amo et al., 2016). Similarly, no major differences in histamine-metabolizing genes were observed, with the exception of a marginally significant lower frequency of homozygous variant genotypes corresponding to the DAO gene variation rs10156191 among patients, which was not significant on comparing allele frequencies (Table 6). The rest of the histamine-metabolism SNPs did not show significant differences when patients and controls were compared, and the genotypes and allele frequencies corresponded with those previously reported for Spanish individuals (García-Martín et al., 2006, 2008, 2015b; Agúndez et al., 2012). Although some gender-related risk associations were identified (summarized in Table 7), in all cases these were due to differences in allele frequencies, but no statistically significant differences for genotypes, in any of the genetic models analyzed, were identified.

In addition to FCER1A genes, we analyzed histamine-metabolism genes because these genes, alone or interacting with FCER1A genes, may be involved in, and hence modulate, the events that occur downstream of reaction triggering. Again, to our knowledge no previous studies addressing the role of histamine-metabolism in these genes in selective NSAID hypersensitivity have been carried out. We previously reported that the non-synonymous variant on the diamine oxidase gene, rs10156191, which causes decreased metabolic capacity, was significantly associated with cross-intolerance to NSAIDs (OR, 1.7; 95% CI, 1.3–2.1; Pc = 0.0003; (Agúndez et al., 2012)). Conversely, in the present study we did not find any association of this SNP with overall selective hypersensitivity. Regarding HNMT, it has been reported that the 939A>G polymorphism, which lowers HNMT enzymatic activity by decreasing HNMT mRNA stability, is associated with aspirin intolerant chronic urticaria (Kim et al., 2009). This aside, no other studies analyzing the possible role of gene-related histamine metabolism in NSAID-hypersensitivity have been published. The lack of association of selective NSAID-hypersensitivity with functional histamine-metabolizing SNPs observed in our study is somewhat unexpected because the culprit NSAIDs implicated in both selective and cross hypersensitivity reactions are the same. However, it should be taken into consideration that the mechanisms involved in cross-intolerance are completely different from those involved in selective hypersensitivity (Agúndez et al., 2012) and this may explain the differences in the linkage of histamine-related genes with the clinical entities.

Phenotype-genotype interaction, however, may be relevant to selective NSAID hypersensitivity: We analyzed by linear regression the putative interaction of all genotypes as well as phenotypic factors. The most relevant phenotypic factor was previous history of atopy, which was strongly related to hypersensitivity to ibuprofen. Age was related to hypersensitivity to ibuprofen. Gender was related to the risk of developing hypersensitivity to ASA, with men showing increased risk, which is unexpected because drug allergy is more frequent in women (Doña et al., 2011). This drug-specific phenotypic feature deserves further investigation. This study also revealed the association of the DAO rs2052129 GG genotype with hypersensitivity to ibuprofen and a weak association of the DAO SNP rs10156191 in homozygosity for the detrimental allele (TT) to an increased risk of developing selective hypersensitivity to ASA, but not to other NSAIDs. The fact that we observed a statistically significant association after multiple regression with only 19 patients with selective hypersensitivity to ASA suggests that the association is strong. Although we should be cautious with regard to this association because of the low number of individuals carrying the TT genotype, further studies are warranted. Additional studies focusing on genes known to be related to organ-specific NSAID-induced hypersensitivity, such as hepatotoxicity (Lucena et al., 2008, 2010; Andrade et al., 2009; Agúndez et al., 2011), should also be conducted in subgroups of selective responders to NSAIDs stratified according to the culprit drug.

Author Contributions

Conceived and designed the experiments: EG Performed the experiments: GA and GE. Analyzed the data: JA and EG. Wrote the paper: EG and JA. Patient assessment: MB, JG, CC, JC, MT, CrM, NB, and GC. Acquisition of data: GA, JG, CC, JC, MT, CaM, NB, GC, GE, MB, JA, and EG. Statistical analysis: AR, JA, and EG. All authors participated in the critical review of the manuscript.

Funding

This study was financed by grants PI12/00241, PI12/00324, PI15/00303, RETICS RD12/0013/0002, and RETICS RD16/0006/0004 from Fondo de Investigación Sanitaria, Instituto de Salud Carlos III, Spain, and GR15026 from Junta de Extremadura, Spain. Financed in part with FEDER funds from the European Union. The authors are grateful to Professor James McCue for assistance in language editing.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

1000 Genomes (2015). A Deep Catalog of Human Genetic Variation. Available online at: http://browser.1000genomes.org/index.html (Accessed July 4, 2016).

Agúndez, J. A., Ayuso, P., Cornejo-García, J. A., Blanca, M., Torres, M. J., Doña, I., et al. (2012). The diamine oxidase gene is associated with hypersensitivity response to non-steroidal anti-inflammatory drugs. PLoS ONE 7:e47571. doi: 10.1371/journal.pone.0047571

PubMed Abstract | CrossRef Full Text | Google Scholar

Agúndez, J. A., García-Martín, E., and Martínez, C. (2009b). Genetically based impairment in CYP2C8- and CYP2C9-dependent NSAID metabolism as a risk factor for gastrointestinal bleeding: is a combination of pharmacogenomics and metabolomics required to improve personalized medicine? Expert Opin. Drug Metab. Toxicol. 5, 607–620. doi: 10.1517/17425250902970998

PubMed Abstract | CrossRef Full Text | Google Scholar

Agúndez, J. A., Lucena, M. I., Martínez, C., Andrade, R. J., Blanca, M., Ayuso, P., et al. (2011). Assessment of nonsteroidal anti-inflammatory drug-induced hepatotoxicity. Expert Opin. Drug Metab. Toxicol. 7, 817–828. doi: 10.1517/17425255.2011.574613

PubMed Abstract | CrossRef Full Text | Google Scholar

Agúndez, J. A., Martínez, C., Pérez-Sala, D., Carballo, M., Torres, M. J., and García-Martín, E. (2009a). Pharmacogenomics in aspirin intolerance. Curr. Drug Metab. 10, 998–1008. doi: 10.2174/138920009790711814

PubMed Abstract | CrossRef Full Text | Google Scholar

Agúndez, J. A., Mayorga, C., and García-Martin, E. (2015). Drug metabolism and hypersensitivity reactions to drugs. Curr. Opin. Allergy Clin. Immunol. 15, 277–284. doi: 10.1097/ACI.0000000000000174

PubMed Abstract | CrossRef Full Text | Google Scholar

Amo, G., García-Menaya, J., Campo, P., Cordobés, C., Plaza Serón, M. C., Ayuso, P., et al. (2016). A nonsynonymous FCER1B SNP is associated with risk of developing allergic rhinitis and with IgE levels. Sci. Rep. 6:19724. doi: 10.1038/srep19724

PubMed Abstract | CrossRef Full Text | Google Scholar

Andrade, R. J., Agúndez, J. A., Lucena, M. I., Martínez, C., Cueto, R., and García-Martín, E. (2009). Pharmacogenomics in drug induced liver injury. Curr. Drug Metab. 10, 956–970. doi: 10.2174/138920009790711805

PubMed Abstract | CrossRef Full Text | Google Scholar

Ariza, A., García-Martín, E., Salas, M., Montańez, M. I., Mayorga, C., Blanca-Lopez, N., et al. (2016). Pyrazolones metabolites are relevant for identifying selective anaphylaxis to metamizole. Sci. Rep. 6:23845. doi: 10.1038/srep23845

PubMed Abstract | CrossRef Full Text | Google Scholar

Ayuso, P., García-Martín, E., Martínez, C., and Agúndez, J. A. (2007). Genetic variability of human diamine oxidase: occurrence of three nonsynonymous polymorphisms and study of their effect on serum enzyme activity. Pharmacogenet. Genomics 17, 687–693. doi: 10.1097/FPC.0b013e328012b8e4

PubMed Abstract | CrossRef Full Text | Google Scholar

Blanca-López, N., Pérez-Alzate, D., Andreu, I., Doña, I., Agúndez, J. A., García-Martín, E., et al. (2016b). Immediate hypersensitivity reactions to ibuprofen and other arylpropionic acid derivatives. Allergy 71, 1048–1056. doi: 10.1111/all.12855

PubMed Abstract | CrossRef Full Text | Google Scholar

Blanca-López, N., Pérez-Sánchez, N., Agúndez, J. A., García-Martin, E., Torres, M. J., Cornejo-Garcia, J. A., et al. (2016a). Allergic reactions to metamizole: immediate and delayed responses. Int. Arch. Allergy Immunol. 169, 223–230. doi: 10.1159/000444798

PubMed Abstract | CrossRef Full Text | Google Scholar

Cornejo-Garcia, J. A., Blanca-López, N., Doña, I., Andreu, I., Agúndez, J. A., Carballo, M., et al. (2009). Hypersensitivity reactions to non-steroidal anti-inflammatory drugs. Curr. Drug Metab. 10, 971–980. doi: 10.2174/138920009790711841

PubMed Abstract | CrossRef Full Text | Google Scholar

Cornejo-García, J. A., Jagemann, L. R., Blanca-López, N., Dońa, I., Flores, C., Guéant-Rodriguez, R. M., et al. (2012). Genetic variants of the arachidonic acid pathway in non-steroidal anti-inflammatory drug-induced acute urticaria. Clin. Exp. Allergy 42, 1772–1781. doi: 10.1111/j.1365-2222.2012.04078.x

PubMed Abstract | CrossRef Full Text | Google Scholar

de Guia, R. M., Echavez, M. D., Gaw, E. L., Gomez, M. R., Lopez, K. A., Mendoza, R. C., et al. (2015). Multifactor-dimensionality reduction reveals interaction of important gene variants involved in allergy. Int. J. Immunogenet. 42, 182–189. doi: 10.1111/iji.12200

PubMed Abstract | CrossRef Full Text | Google Scholar

Doña, I., Blanca-López, N., Cornejo-García, J. A., Torres, M. J., Laguna, J. J., Fernández, J., et al. (2011). Characteristics of subjects experiencing hypersensitivity to non-steroidal anti-inflammatory drugs: patterns of response. Clin. Exp. Allergy 41, 86–95. doi: 10.1111/j.1365-2222.2010.03651.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Doña, I., Blanca-López, N., Torres, M. J., García-Campos, J., Garcia-Nuñez, I., Gomez, F., et al. (2012). Drug hypersensitivity reactions: response patterns, drug involved, and temporal variations in a large series of patients. J. Invest. Allergol. Clin. Immunol. 22, 363–371.

Google Scholar

García-Martín, E., Ayuso, P., Martínez, C., and Agúndez, J. A. (2007b). Improved analytical sensitivity reveals the occurrence of gender-related variability in diamine oxidase enzyme activity in healthy individuals. Clin. Biochem. 40, 1339–1341. doi: 10.1016/j.clinbiochem.2007.07.019

PubMed Abstract | CrossRef Full Text | Google Scholar

García-Martín, E., Ayuso, P., Martínez, C., Blanca, M., and Agúndez, J. A. (2009). Histamine pharmacogenomics. Pharmacogenomics 10, 867–883. doi: 10.2217/pgs.09.26

PubMed Abstract | CrossRef Full Text | Google Scholar

García-Martín, E., Esguevillas, G., Blanca-López, N., García-Menaya, J., Blanca, M., Amo, G., et al. (2015a). Genetic determinants of metamizole metabolism modify the risk of developing anaphylaxis. Pharmacogenet. Genomics 25, 462–464. doi: 10.1097/fpc.0000000000000157

PubMed Abstract | CrossRef Full Text | Google Scholar

García-Martín, E., García-Menaya, J., Sánchez, B., Martínez, C., Rosendo, R., and Agundez, J. A. (2007a). Polymorphisms of histamine-metabolizing enzymes and clinical manifestations of asthma and allergic rhinitis. Clin. Exp. Allergy 37, 1175–1182. doi: 10.1111/j.1365-2222.2007.02769.x

PubMed Abstract | CrossRef Full Text | Google Scholar

García-Martín, E., Martínez, C., Serrador, M., Alonso-Navarro, H., Ayuso, P., Navacerrada, F., et al. (2015b). Diamine oxidase rs10156191 and rs2052129 variants are associated with the risk for migraine. Headache 55, 276–286. doi: 10.1111/head.12493

PubMed Abstract | CrossRef Full Text | Google Scholar

García-Martín, E., Mart´inez, C., Serrador, M., Alonso-Navarro, H., Navacerrada, F., Agúndez, J. A., et al. (2008). Histamine-N-methyl transferase polymorphism and risk for migraine. Headache 48, 1343–1348. doi: 10.1111/j.1526-4610.2007.01056.x

PubMed Abstract | CrossRef Full Text | Google Scholar

García-Martín, E., Mendoza, J. L., Martinez, C., Taxonera, C., Urcelay, E., Ladero, J. M., et al. (2006). Severity of ulcerative colitis is associated with a polymorphism at diamine oxidase gene but not at histamine N-methyltransferase gene. World J. Gastroenterol. 12, 615–620. doi: 10.3748/wjg.v12.i4.615

PubMed Abstract | CrossRef Full Text | Google Scholar

Gervasini, G., Agúndez, J. A., García-Menaya, J., Martínez, C., Cordobés, C., Ayuso, P., et al. (2010). Variability of the L-Histidine decarboxylase gene in allergic rhinitis. Allergy 65, 1576–1584. doi: 10.1111/j.1398-9995.2010.02425.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, A., Zhu, W., Zhang, C., Wen, S., Chen, X., Chen, M., et al. (2015). Association of FCER1A genetic polymorphisms with risk for chronic spontaneous urticaria and efficacy of nonsedating H1-antihistamines in Chinese patients. Arch. Dermatol. Res. 307, 183–190. doi: 10.1007/s00403-014-1525-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Hua, L., Zuo, X. B., Bao, Y. X., Liu, Q. H., Li, J. Y., Lv, J., et al. (2016). Four-locus gene interaction between IL13, IL4, FCER1B, and ADRB2 for asthma in Chinese Han children. Pediatr. Pulmonol. 51, 364–371. doi: 10.1002/ppul.23322

PubMed Abstract | CrossRef Full Text | Google Scholar

Kennedy, M. J., Loehle, J. A., Griffin, A. R., Doll, M. A., Kearns, G. L., Sullivan, J. E., et al. (2008). Association of the histamine N-methyltransferase C314T (Thr105Ile) polymorphism with atopic dermatitis in Caucasian children. Pharmacotherapy 28, 1495–1501. doi: 10.1592/phco.28.12.1495

PubMed Abstract | CrossRef Full Text | Google Scholar

Khan, D. A., and Solensky, R. (2010). Drug allergy. J. Allergy Clin. Immunol. 125, S126–S137. doi: 10.1016/j.jaci.2009.10.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, S. H., Kang, Y. M., Kim, S. H., Cho, B. Y., Ye, Y. M., Hur, G. Y., et al. (2009). Histamine N-methyltransferase 939A>G polymorphism affects mRNA stability in patients with acetylsalicylic acid-intolerant chronic urticaria. Allergy 64, 213–221. doi: 10.1111/j.1398-9995.2008.01795.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kowalski, M. L., Asero, R., Bavbek, S., Blanca, M., Blanca-Lopez, N., Bochenek, G., et al. (2013). Classification and practical approach to the diagnosis and management of hypersensitivity to nonsteroidal anti-inflammatory drugs. Allergy 68, 1219–1232. doi: 10.1111/all.12260

PubMed Abstract | CrossRef Full Text | Google Scholar

Kowalski, M. L., Makowska, J. S., Blanca, M., Bavbek, S., Bochenek, G., Bousquet, J., et al. (2011). Hypersensitivity to nonsteroidal anti-inflammatory drugs (NSAIDs) - classification, diagnosis and management: review of the EAACI/ENDA(#) and GA2LEN/HANNA*. Allergy 66, 818–829. doi: 10.1111/j.1398-9995.2011.02557.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Gibson, A. W., and Kimberly, R. P. (2014). Human FcR polymorphism and disease. Curr. Top. Microbiol. Immunol. 382, 275–302. doi: 10.1007/978-3-319-07911-0_13

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X., Wang, G., Hong, X., Wang, D., Tsai, H. J., Zhang, S., et al. (2011). Gene-vitamin D interactions on food sensitization: a prospective birth cohort study. Allergy 66, 1442–1448. doi: 10.1111/j.1398-9995.2011.02681.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lucena, M. I., Andrade, R. J., Martínez, C., Ulzurrun, E., García-Martín, E., Borraz, Y., et al. (2008). Spanish group for the study of drug-induced liver disease, glutathione S-transferase m1 and t1 null genotypes increase susceptibility to idiosyncratic drug-induced liver injury. Hepatology 48, 588–596. doi: 10.1002/hep.22370

PubMed Abstract | CrossRef Full Text | Google Scholar

Lucena, M. I., García-Martín, E., Andrade, R., J., Martínez, C., Stephens, C., Agundez, J. A., et al. (2010). Mitochondrial superoxide dismutase and glutathione peroxidase in idiosyncratic drug-induced liver injury. Hepatology 52, 303–312. doi: 10.1002/hep.23668

PubMed Abstract | CrossRef Full Text | Google Scholar

MacGlashan, D. Jr., Lichtenstein, L. M., McKenzie-White, J., Chichester, K., Henry, A. J., Sutton, B. J., et al. (1999). Upregulation of FcepsilonRI on human basophils by IgE antibody is mediated by interaction of IgE with FcepsilonRI. J. Allergy Clin. Immunol. 104, 492–498. doi: 10.1016/S0091-6749(99)70399-4

PubMed Abstract | CrossRef Full Text | Google Scholar

MacGlashan, D. Jr., McKenzie-White, J., Chichester, K., Bochner, B. S., Davis, F. M., Schroeder, J. T., et al. (1998). In vitro regulation of FcepsilonRIalpha expression on human basophils by IgE antibody. Blood 91, 1633–1643.

PubMed Abstract | Google Scholar

Maintz, L., Yu, C. F., Rodriguez, E., Baurecht, H., Bieber, T., Illig, T., et al. (2011). Association of single nucleotide polymorphisms in the diamine oxidase gene with diamine oxidase serum activities. Allergy 66, 893–902. doi: 10.1111/j.1398-9995.2011.02548.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Manikandan, J., Kothandaraman, N., Hande, M. P., and Pushparaj, P. N. (2012). Deciphering the structure and function of FcepsilonRI/mast cell axis in the regulation of allergy and anaphylaxis: a functional genomics paradigm. Cell. Mol. Life Sci. 69, 1917–1929. doi: 10.1007/s00018-011-0886-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Martinez, C., Andreu, I., Amo, G., Miranda, M. A., Esguevillas, G., Torres, M. J., et al. (2014). Gender and functional CYP2C and NAT2 polymorphisms determine the metabolic profile of metamizole. Biochem. Pharmacol. 92, 457–466. doi: 10.1016/j.bcp.2014.09.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Mita, H., Endoh, S., Kudoh, M., Kawagishi, Y., Kobayashi, M., Taniguchi, M., et al. (2001). Possible involvement of mast-cell activation in aspirin provocation of aspirin-induced asthma. Allergy 56, 1061–1067. doi: 10.1111/j.1398-9995.2001.00913.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Palikhe, N. S., Kim, S. H., Cho, B. Y., Ye, Y. M., Hur, G. Y., and Park, H. S. (2008b). Association of three sets of high-affinity IgE receptor (FcepsilonR1) polymorphisms with aspirin-intolerant asthma. Respir. Med. 102, 1132–1139. doi: 10.1016/j.rmed.2008.03.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Palikhe, N. S., Kim, S. H., and Park, H. S. (2008a). What do we know about the genetics of aspirin intolerance? J. Clin. Pharm. Ther. 33, 465–472. doi: 10.1111/j.1365-2710.2008.00961.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pirmohamed, M. (2006). Genetic factors in the predisposition to drug-induced hypersensitivity reactions. AAPS J. 8, E20–E26. doi: 10.1208/aapsj080103

PubMed Abstract | CrossRef Full Text | Google Scholar

Potaczek, D. P., and Kabesch, M. (2012). Current concepts of IgE regulation and impact of genetic determinants. Clin. Exp. Allergy 42, 852–871. doi: 10.1111/j.1365-2222.2011.03953.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Preuss, C. V., Wood, T. C., Szumlanski, C. L., Raftogianis, R. B., Otterness, D. M., Girard, B., et al. (1998). Human histamine N-methyltransferase pharmacogenetics: common genetic polymorphisms that alter activity. Mol. Pharmacol. 53, 708–717.

PubMed Abstract | Google Scholar

Ramphul, K., Lv, J., Hua, L., Liu, Q. H., Fang, D. Z., Ji, R. X., et al. (2014). Single nucleotide polymorphisms predisposing to asthma in children of Mauritian Indian and Chinese Han ethnicity. Braz. J. Med. Biol. Res. 47, 394–397. doi: 10.1590/1414-431x20143751

PubMed Abstract | CrossRef Full Text | Google Scholar

Saini, S. S., MacGlashan, D. W. Jr., Sterbinsky, S. A., Togias, A., Adelman, D. C., Lichtenstein, L. M., et al. (1999). Down-regulation of human basophil IgE and FC epsilon RI alpha surface densities and mediator release by anti-IgE-infusions is reversible in vitro and in vivo. J. Immunol. 162, 5624–5630.

PubMed Abstract | Google Scholar

Solé, X., Guin´o, E., Valls, J., Iniesta, R., and Moreno, V. (2006). SNPStats: a web tool for the analysis of association studies. Bioinformatics 22, 1928–1929. doi: 10.1093/bioinformatics/btl268

PubMed Abstract | CrossRef Full Text | Google Scholar

Weidinger, S., Gieger, C., Rodriguez, E., Baurecht, H., Mempel, M., Klopp, N., et al. (2008). Genome-wide scan on total serum IgE levels identifies FCER1A as novel susceptibility locus. PLoS Genet. 4:e1000166. doi: 10.1371/journal.pgen.1000166

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, H. J., Zheng, L., Zhang, X. F., Yang, M., and Huang, X. (2014). Association of the MS4A2 gene promoter C-109T or the 7th exon E237G polymorphisms with asthma risk: a meta-analysis. Clin. Biochem. 47, 605–611. doi: 10.1016/j.clinbiochem.2014.01.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Lu, M. M., Lu, Y. W., Feng, C. C., Leng, R. X., Pan, H. F., et al. (2013). Sex-specific differences in the relationship between the single-nucleotide polymorphism rs2298804 of FCER1A and the susceptibility to systemic lupus erythematosus in a Chinese Han population. Clin. Exp. Dermatol. 38, 410–416. doi: 10.1111/ced.12035

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Duan, S., Lin, X., Zhang, W., Meng, N., Zhao, L., et al. (2010). SNPs in the FCER1A gene region show no association with allergic rhinitis in a Han Chinese population. PLoS ONE 5:e15792. doi: 10.1371/journal.pone.0015792

PubMed Abstract | CrossRef Full Text | Google Scholar