The SNP rs3128965 of HLA-DPB1 as a Genetic Marker of the AERD Phenotype

Abstract

Background

Two common clinical syndromes of acetylsalicylic acid (aspirin) hypersensitivity, aspirin-exacerbated respiratory disease (AERD) and aspirin-exacerbated cutaneous disease (AECD), were subjected to a genome-wide association study to identify strong genetic markers for aspirin hypersensitivity in a Korean population.

Methods

A comparison of SNP genotype frequencies on an Affymetrix Genome-Wide Human SNP array of 179 AERD patients and 1989 healthy normal control subjects (NC) revealed SNPs on chromosome 6 that were associated with AERD, but not AECD. To validate the association, we enrolled a second cohort comprising AERD (n = 264), NC (n = 238) and disease-control (aspirin tolerant asthma; ATA, n = 387) groups.

Results

The minor genotype frequency (AG or AA) of a particular SNP, rs3128965, in the HLA-DPB1 region was higher in the AERD group compared to the ATA or NC group (P = 0.001, P = 0.002, in a co-dominant analysis model, respectively). Comparison of rs3128965 alleles with the clinical features of asthmatics revealed that patients harboring the A allele had increased bronchial hyperresponsiveness to inhaled aspirin and methacholine, and higher 15-HETE levels, than those without the A allele (P = 0.039, 0.037, and 0.004, respectively).

Conclusions

This implies the potential of rs3128965 as a genetic marker for diagnosis and prediction of the AERD phenotype.

About the Authors:

Seung-Hyun Kim

* E-mail: [email protected] (SHK); [email protected] (HSP)

Affiliation: Department of Allergy & Clinical Immunology, Ajou University School of Medicine, Suwon, Korea,

Bo-Young Cho

Affiliation: Department of Allergy & Clinical Immunology, Ajou University School of Medicine, Suwon, Korea,

Hyunna Choi

Affiliation: Department of Allergy & Clinical Immunology, Ajou University School of Medicine, Suwon, Korea,

Eun-Soon Shin

Affiliation: DNA link, Inc., Seoul, Korea

Young-Min Ye

Affiliation: Department of Allergy & Clinical Immunology, Ajou University School of Medicine, Suwon, Korea,

Jong-Eun Lee

Affiliation: DNA link, Inc., Seoul, Korea

Hae-Sim Park

* E-mail: [email protected] (SHK); [email protected] (HSP)

Affiliation: Department of Allergy & Clinical Immunology, Ajou University School of Medicine, Suwon, Korea,

Introduction

Hypersensitivity reactions to aspirin manifest as broad-ranging clinical features of allergic reactions including aspirin-exacerbated respiratory disease (AERD), aspirin-exacerbated cutaneous disease (AECD), and occasionally, anaphylaxis [1], [2]. A genetic association study reported several genetic risk factors for aspirin hypersensitivity, with most being involved in arachidonate metabolism or the inflammatory pathway [3], [4] To date, the best genetic marker for AERD is HLA-DPB1*0301, which is also associated with a higher prevalence of chronic rhinosinusitis and a higher leukotriene receptor antagonist dose to control asthmatic symptoms in AERD patients [5]–[7].

To identify additional genetic determinants for increased risk of developing aspirin hypersensitivity in adult patients with underlying asthma or urticaria, we performed a genome-wide association study (GWAS) using an Affymetrix Genome-Wide Human SNP array.

To verify any significant associations, a second replication study was carried out using a larger study cohort in a Korean population.

Here we report a newly identified single-nucleotide polymorphism (SNP), rs3128965, as a genetic marker of the AERD phenotype.

Methods

Study subjects

Patients with aspirin hypersensitivity (n = 390; 179 AERD, 211 AECD) were recruited from the Ajou University Medical Center, Suwon, Korea. AERD patients were diagnosed by a positive lysine-aspirin bronchoprovocation test (L-ASA BPT) from asthmatic patients who showed typical clinical features of asthma [8]. Lys-ASA BPT was performed as described previously with increasing dosages of Lys-ASA (75 mg, 150 mg, 300 mg) [8]. The AERD patients had upper and lower respiratory symptoms without any history of drug-induced cutaneous symptoms, while the ATA patients were defined as asthmatic patients showing a negative response to L-ASA BPT. The AECD patients were defined as having > 6 weeks of daily itchy wheals, a history of chronic urticaria/angioedema aggravated by ingestion of more than two nonsteroidal anti-inflammatory drugs, and showing positive responses in the oral aspirin challenge test. Blood samples and clinical information from patients and controls were collected with written informed consent and with the approval of the Ethics Review Board of the Ajou Medical Center (AJIRB-GEN-GEN-11-304).

Clinical characterization of AERD and AECD patients including rhinosinusitis, pulmonary functional test, and aspirin provocation test were determined as described previously [9]. Serum total IgE and eosinophil cationic protein (ECP) levels were measured using the ImmunoCAP system (Phadia, AB, Uppsala, Sweden). The presence of rhinosinusitis was determined based on paranasal sinus X-rays, CT scans, and/or rhinoscopy. The presence of moderate to chronic rhinosinusitis was defined as a Lund-Mackay score higher than 3 [10]. The 15-hydroxyeicosatetraenoic acid (15-HETE) level was measured in the cell supernatants of peripheral blood leukocytes treated with arachidonic acid (40 µM) at 37°C for 5 min, as described previously [11], [12], and quantified using an EIA kit (Cayman Chemicals, Ann Arbor, MI, USA).

Genome-wide association and the replication studies

Using an Affymetrix genome-wide human SNP array, a GWAS was performed to identify the genetic risk factors for aspirin hypersensitivity in adult patients with underlying asthma or urticaria (AERD, n = 179; AECD, n = 211). Anonymous genotype data from 1989 healthy normal control subjects (NC) were supplied by the Korean National Institute of Health (http://www.cdc.go.kr). The clinical characteristics of the study subjects used in the genome-wide association study are listed in S1 Table.

For the initial analysis of the unrelated case-control samples, we examined the identity-by-state (IBS) counts of the whole-genome genotype data, filtered out samples with <97% genotype call rate, and excluded samples that shared >10% of their genome. There was no evidence of differential genotyping between patients and controls, and no evidence of a population substructure in either group (data not shown). We used the Bayesian Robust Linear Modeling with Mahalanobis distance classifier (BRLMM) genotype calling algorithm to assess the SNPs. SNP markers with low genotyping call rates (<97%), low minor-allele frequencies (<1%) and deviation from Hardy-Weinberg equilibrium (P<0.0001) were filtered out. The association between each SNP and the case-control was assessed using a Cochran–Armitage trend test using PLINK (v1.07) [13]. The quantile-quantile (Q–Q) plot was plotted with the expected distribution of association test statistics under null distribution across the Cochran– Armitage trend observed P-values with R statistical environment version 2.9.0. A Manhattan plot was generated using PLINK. Q-Q plots are shown in S1 Fig. (AERD vs. NC) and S2 Fig. (AECD vs. NC).

Targeted SNP genotyping was also performed using a larger study cohort composing AERD (n = 264), ATA (n = 387), and NC groups (n = 238). The clinical characteristics of the study subjects used in the second-stage validation study are listed in S2 Table. SNP genotyping was performed by the primer extension method using a SNAPshot ddNTP primer extension kit (Applied Biosystems, Foster City, CA, USA). The sequences of the amplifying and extension primers were as follows: forward, 5′-AGCCTCATTTCCTCAAAAA-3′; reverse, 5′-ACAAGCAGGAAAGAGAATGA-3′; and extension, 5′-AGGAAACTTGTGAGAAACCTATGCA-3′. The genotype distributions of the AERD patients and the control subjects in the second-stage replication analyses were assessed using multiple logistic regression models (dominant, recessive and co-dominant), with age and gender as covariates. All statistical analyses were performed using SPSS, version 12.0 (SPSS, Chicago, IL, USA).

Results

After filtering for study subjects with an average call rate of 97%, and markers with high genotyping call rates (>97%), high minor allele frequencies (>1%) and no deviation from Hardy-Weinberg equilibrium (P>0.0001), we included genotype data from 275,862 SNPs from 2379 study subjects (179 AERD, 211 AECD, 1989 NC). The mean age of AERD patients was 43.72 years (±13.09), that of AECD patients was 36.09 years (±11.14), and that of healthy normal controls was 47.67 years (±7.13) (S1 Table). There was a significant difference between the mean age of the case and NC groups (<0.001), but no significant difference in gender proportions between the two groups.

A genome-wide Manhattan plot was drawn using the chromosomal positions of individual SNPs (x-axis) and the negative logarithm of P values calculated using the Cochran–Armitage trend test (y-axis; Fig. 1). The top 20 SNPs associated with AERD are shown in Table 1. A genome-wide Manhattan plot of GWAS in a comparison of AECD with NC is shown in S3 Fig. The top 20 SNPs associated with AERD are shown in S3 Table.

[Figure omitted. See PDF.]

Figure 1. A genome-wide Manhattan plot of GWAS in a comparison of AERD with NC.

Distribution of negative logarithm P values are plotted against chromosomes.

https://doi.org/10.1371/journal.pone.0111220.g001

[Figure omitted. See PDF.]

Table 1. Genome-wide association result for the top 20 single-nucleotide polymorphisms (SNPs) in a comparison of AERD with NC.

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

Two SNPs located on chromosome 6, rs3117230 and rs3128965, showed a strong association with the AERD phenotype and were found to be in perfect linkage disequilibrium (D′ = 0.984, r2 = 0.961). However, none of SNPs listed in Table 1 reached significance at a genome-wide threshold after Bonferroni multiple correction (P<1.8×10−7). The SNP rs3128965, located in the HLA-DPB1 region of chromosome 6, was targeted for a second-stage replication study with a larger cohort due to its specific association with the AERD phenotype, but not the AECD phenotype (Fig. 2). The genome-wide significance of rs3128965 was replicated in the second study cohort (Table 2). The genotype frequency of the minor allele (AG or AA) was significantly higher in the AERD group compared to the disease control (ATA) and NC groups (P = 0.001, P = 0.002, in a co-dominant analysis model, respectively). The significance remained after multiple correction testing.

[Figure omitted. See PDF.]

Figure 2. SNP rs3128965 of HLA-DPB1 as a target SNP for susceptibility to AERD.

Negative logarithm P values are shown around the HLA-DPB1 region. Closed circles indicate the P value of each SNP in a comparison of AERD and NC. Open circles indicate the P value of each SNP in a comparison of AECD and NC.

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

[Figure omitted. See PDF.]

Table 2. Genetic association of the SNP, rs3128965, with AERD phenotype.

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

The clinical features associated with this SNP are shown in Table 3. Adult asthmatic patients with risk allele A showed enhanced bronchial hyperresponsiveness to lysine-aspirin inhaled (% fall of FEV after L-ASA BPT) and methacholine (PC20), and higher 15-HETE levels than those without the A allele (P = 0.039, P = 0.037 and P = 0.004, respectively; Table 3).

[Figure omitted. See PDF.]

Table 3. Clinical features according to the SNP, rs3128965, in asthmatics (264 AERD and 387 ATA).

https://doi.org/10.1371/journal.pone.0111220.t003

Discussion

Several genetic studies on AERD have been conducted [3], [4] and the HLA-DPB1*0301 was reported to have a strong association with the AERD phenotype [5]. The association of DPB1*0301 with AERD was also identified in a Korean population, together with a significant association with the requirement for a leukotriene receptor antagonist [6], [7]. The HLA haplotype HLA-DRB1*1301–DQB1*0609 was identified as a genetic risk factor for AECD in Koreans [14]. Various immunological and genetic mechanisms may contribute to the development of the clinical features of aspirin hypersensitivity.

In this study, rs3128965 was identified as a genetic marker for AERD in an initial GWAS, and its association was validated significance in a subsequent replication study. Recently, two GWAS suggested two novel SNPs, rs7572857 of CEP68 [15] and rs1042151 of HLA-DPB1 [16] as genetic markers for AERD susceptibility in a Korean population. CEP68 rs7572857 was not associated with the AERD phenotype in the present study cohort, while HLA-DPB1 rs1042151 was in perfect linkage disequilibrium with rs3128965, the target SNP in the present study. The replicated association of HLA-DPB1 gene polymorphisms with AERD phenotype suggests that the HLA-DPB1 gene has an important role in AERD susceptibility. Moreover, this study identified the specific influence of HLA-DPB1 rs3128965 on the susceptibility to AERD, but not AECD, suggesting a different genetic mechanism involved in the pathogenesis of AECD.

In addition, rs3128965 showed a significant association with the percentage decrease in FEV1 after L-ASA BPT, which is the gold standard for diagnosis of AERD [17]. This implies that rs3128965 is a reliable genetic marker for diagnosis of aspirin hypersensitivity in adult asthmatic patients.

The significant association of rs3128965 with increased 15-HETE levels, a potential biomarker for AERD patients for in vitro diagnosis [12], [18], and enhanced airway hyperresponsiveness, methacholine PC20, demonstrated that this SNP could be useful in predicting the progression of asthmatic symptoms in AERD patients after exposure to aspirin.

AERD patients are more likely to present with severe or difficult-to-treat asthma because of poor lung function, life threatening exacerbation, and higher requirements for therapeutics such as steroid and leukotriene receptor antagonists (LTRAs) [19], [20]. A previous report suggested use of the HLA-DPB1*0301 marker for predicting LTRA requirement in AERD patients [7]. In the present study, rs3128965 was associated with LTRA requirement (S4 Table), suggesting that this SNP could be a potential marker for predicting favorable response to LTRA treatment in AERD patients.

This is the first GWAS to demonstrate differential genetic susceptibility to two different clinical phenotypes, AERD and AECD. However, this study does have a number of limitations. First, the sample size was relatively small for a GWAS. To exclude false positive associations, follow-up targeted SNP genotyping was performed. However, other candidate SNPs detected on the initial GWAS also require validation by replication studies. Second, this study did not identify any novel genes suggested to confer genetic susceptibility to aspirin, and the limited sample size could not guarantee target SNPs with genome-wide significance. Therefore future replication studies are required.

Taken together, our results validate the use of rs3128965 as a genetic marker to predict the AERD phenotype, both in vivo (greater bronchoconstriction to aspirin exposure) and in vitro (greater production of 15-HETE from peripheral blood cells). These findings also confirm the usefulness of SNP genotyping for diagnosis and prediction of the AERD phenotype.

Supporting Information

[Figure omitted. See PDF.]

S1 Fig.

Q-Q plot of GWAS data from the patients with AERD vs. NC.

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

(TIF)

S2 Fig.

Q-Q plot of GWAS data from the patients with AECD vs. NC.

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

(TIF)

S3 Fig.

The Manhattan plot of GWAS shows P values for AECD vs. NC using a Cochran–Armitage trend test.

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

(TIF)

S1 Table.

Clinical demographics of the study subjects enrolled in the first genome-wide association study.

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

(DOCX)

S2 Table.

Clinical demographics of the study subjects enrolled in the second validation study.

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

(DOCX)

S3 Table.

Genome-wide association result for the top 20 single-nucleotide polymorphisms (SNPs) in a comparison of AECD and NC.

https://doi.org/10.1371/journal.pone.0111220.s006

(DOCX)

S4 Table.

Association of the SNP rs3128965 with leukotriene receptor antagonist (LTRA) requirement in 75 AERD patients.

https://doi.org/10.1371/journal.pone.0111220.s007

(DOCX)

Acknowledgments

We thank the Korean National Institute of Health (http://www.cdc.go.kr) for generous provisions of GWAS data from normal, healthy volunteers.

Author Contributions

Conceived and designed the experiments: SHK HSP. Performed the experiments: SHK BYC HC ESS. Analyzed the data: SHK BYC HC ESS. Contributed reagents/materials/analysis tools: HSP YMY SHK JEL. Wrote the paper: SHK HSP.

Citation: Kim S-H, Cho B-Y, Choi H, Shin E-S, Ye Y-M, Lee J-E, et al. (2014) The SNP rs3128965 of HLA-DPB1 as a Genetic Marker of the AERD Phenotype. PLoS ONE 9(12): e111220. https://doi.org/10.1371/journal.pone.0111220

References

1. Ayuso P, Blanca-Lopez N, Dona I, Torres MJ, Gueant-Rodriguez RM, et al. (2013) Advanced phenotyping in hypersensitivity drug reactions to NSAIDs. Clin Exp Allergy 43:1097–1109.

2. Sanchez-Borges M, Caballero-Fonseca F, Capriles-Hulett A (2013) Aspirin-exacerbated cutaneous disease. Immunol Allergy Clin North Am 33:251–262.

3. Kim SH, Sanak M, Park HS (2013) Genetics of hypersensitivity to aspirin and nonsteroidal anti-inflammatory drugs. Immunol Allergy Clin North Am 33:177–194.

4. Park SM, Park JS, Park HS, Park CS (2013) Unraveling the genetic basis of aspirin hypersensitivity in asthma beyond arachidonate pathways. Allergy Asthma Immunol Res 5:258–276.

5. Dekker JW, Nizankowska E, Schmitz-Schumann M, Pile K, Bochenek G, et al. (1997) Aspirin-induced asthma and HLA-DRB1 and HLA-DPB1 genotypes. Clin Exp Allergy 27:574–577.

6. Choi JH, Lee KW, Oh HB, Lee KJ, Suh YJ, et al. (2004) HLA association in aspirin-intolerant asthma: DPB1*0301 as a strong marker in a Korean population. J Allergy Clin Immunol 113:562–564.

7. Park HS, Kim SH, Sampson AP, Lee KW, Park CS (2004) The HLA-DPB1*0301 marker might predict the requirement for leukotriene receptor antagonist in patients with aspirin-intolerant asthma. J Allergy Clin Immunol 114:688–689.

8. Park HS (1995) Early and late onset asthmatic responses following lysine-aspirin inhalation in aspirin-sensitive asthmatic patients. Clin Exp Allergy 25:38–40.

9. Kim SH, Yang EM, Lee HN, Choi GS, Ye YM, et al. (2010) Association of the CCR3 gene polymorphism with aspirin exacerbated respiratory disease. Respir Med 104:626–632.

10. Hopkins C, Browne JP, Slack R, Lund V, Brown P (2007) The Lund-Mackay staging system for chronic rhinosinusitis: how is it used and what does it predict? Otolaryngol Head Neck Surg 137:555–561.

11. Liu C, Xu D, Liu L, Schain F, Brunnstrom A, et al. (2009) 15-Lipoxygenase-1 induces expression and release of chemokines in cultured human lung epithelial cells. American journal of physiology Lung cellular and molecular physiology 297:L196–203.

12. Kowalski ML, Ptasinska A, Jedrzejczak M, Bienkiewicz B, Cieslak M, et al. (2005) Aspirin-triggered 15-HETE generation in peripheral blood leukocytes is a specific and sensitive Aspirin-Sensitive Patients Identification Test (ASPITest). Allergy 60:1139–1145.

13. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. American journal of human genetics 81:559–575.

14. Kim SH, Choi JH, Lee KW, Kim SH, Shin ES, et al. (2005) The human leucocyte antigen-DRB1*1302-DQB1*0609-DPB1*0201 haplotype may be a strong genetic marker for aspirin-induced urticaria. Clin Exp Allergy 35:339–344.

15. Kim JH, Park BL, Cheong HS, Bae JS, Park JS, et al. (2010) Genome-wide and follow-up studies identify CEP68 gene variants associated with risk of aspirin-intolerant asthma. PLoS One 5:e13818.

16. Park BL, Kim TH, Kim JH, Bae JS, Pasaje CF, et al. (2013) Genome-wide association study of aspirin-exacerbated respiratory disease in a Korean population. Human genetics 132:313–321.

17. Nizankowska-Mogilnicka E, Bochenek G, Mastalerz L, Swierczynska M, Picado C, et al. (2007) EAACI/GA2LEN guideline: aspirin provocation tests for diagnosis of aspirin hypersensitivity. Allergy 62:1111–1118.

18. Korosec P, Tisler U, Bajrovic N, Silar M, Mrhar A, et al. (2011) Acetylsalicylic acid-triggered 15-HETE generation by peripheral leukocytes for identifying ASA sensitivity. Respir Med 105 Suppl 1S81–83.

19. Mascia K, Haselkorn T, Deniz YM, Miller DP, Bleecker ER, et al. (2005) Aspirin sensitivity and severity of asthma: evidence for irreversible airway obstruction in patients with severe or difficult-to-treat asthma. J Allergy Clin Immunol 116:970–975.

20. Wenzel S (2012) Severe asthma: from characteristics to phenotypes to endotypes. Clin Exp Allergy 42:650–658.