Figures
Abstract
Aspirin-intolerant asthma (AIA) is a rare condition that is characterized by the development of bronchoconstriction in asthmatic patients after ingestion of non-steroidal anti-inflammatory drugs including aspirin. However, the underlying mechanisms of AIA occurrence are still not fully understood. To identify the genetic variations associated with aspirin intolerance in asthmatics, the first stage of genome-wide association study with 109,365 single nucleotide polymorphisms (SNPs) was undertaken in a Korean AIA (n = 80) cohort and aspirin-tolerant asthma (ATA, n = 100) subjects as controls. For the second stage of follow-up study, 150 common SNPs from 11 candidate genes were genotyped in 163 AIA patients including intermediate AIA (AIA-I) subjects and 429 ATA controls. Among 11 candidate genes, multivariate logistic analyses showed that SNPs of CEP68 gene showed the most significant association with aspirin intolerance (P values of co-dominant for CEP68, 6.0×10−5 to 4.0×10−5). All seven SNPs of the CEP68 gene showed linkage disequilibrium (LD), and the haplotype of CEP68_ht4 (T-G-A-A-A-C-G) showed a highly significant association with aspirin intolerance (OR = 2.63; 95% CI = 1.64–4.21; P = 6.0×10−5). Moreover, the nonsynonymous CEP68 rs7572857G>A variant that replaces glycine with serine showed a higher decline of forced expiratory volume in 1s (FEV1) by aspirin provocation than other variants (P = 3.0×10−5). Our findings imply that CEP68 could be a susceptible gene for aspirin intolerance in asthmatics, suggesting that the nonsynonymous Gly74Ser could affect the polarity of the protein structure.
Citation: Kim J-H, Park B-L, Cheong HS, Bae JS, Park JS, Jang AS, et al. (2010) Genome-Wide and Follow-Up Studies Identify CEP68 Gene Variants Associated with Risk of Aspirin-Intolerant Asthma. PLoS ONE 5(11): e13818. https://doi.org/10.1371/journal.pone.0013818
Editor: Matthias Wjst, GSF National Research Center for Environment and Health, Germany
Received: June 29, 2010; Accepted: October 16, 2010; Published: November 3, 2010
Copyright: © 2010 Kim et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: Hyoung Doo Shin received grant number M1-0302-00-0073 from the Korea Science and Engineering Foundation (KOSEF) funded by the Korean government (MEST) (No. 2009-0080157), and a Priority Research Centers Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (2009-0093822). Choon-Sik Park received a grant from the Korea Health 21 R&D Project (A010249). SNP Genetics, Inc. provided the Beadstation500 scanner instrument and BeadStudio 3.0 software used in the research. They also were involved in the study design, data collection and analysis, decision to publish and preparation of the manuscript.
Competing interests: SNP Genetics, Inc. provided the Beadstation500 scanner instrument and BeadStudio 3.0 software used in the research. However, these competing interests do not alter the authors' adherence to all the PLoS ONE policies. No other relevant interests are declared.
Introduction
Aspirin, also known as acetylsalicylic acid, has been used as an analgesic to relieve pain and fever, as well as an anti-inflammatory medication. Up to 20% of asthmatics are sensitive to aspirin and other non-steroidal anti-inflammatory drugs (NSAIDs). Aspirin-intolerant asthma (AIA), as a unique clinical syndrome with acute bronchospasm after ingestion of aspirin or other NSAIDs, was first described in 1922 [1], [2]. AIA is characterized by the triad of aspirin hypersensitivity, bronchial asthma, and chronic rhinosinusitis with nasal polyposis [3], [4]. AIA progression from the upper to the lower respiratory tract is accompanied by persistent asthmatic symptoms with intense eosinophilic infiltration into the upper and lower airways. AIA generally begins at about 30 years of age and occurs more frequently in women [2], [5].
Recent studies have shown that aspirin hypersensitivity is likely to be related with over-production of pro-inflammatory cysteinyl leukotrienes (CysLTs) such as LTC4, LTD4, and LTE4 [2]. CysLTs are generated from arachidonic acid via the leukotriene synthetic pathway or the 5-lipoxygenase pathway. Aspirin blocks the cyclooxygenase pathway that converts arachidonic acid into prostaglandins, and eventually diverts arachidonic acid metabolites to the 5-lipoxygenase pathway [6], [7]. Therefore, CysLT receptors and arachidonic acid cascade-related enzymes/proteins have been considered as main targets for aspirin intolerance in asthmatics. CysLT receptors are selectively antagonized by several leukotriene modifiers, such as montelukast, pranlukast, and zafirlukast [8], [9]; however, clinical studies have demonstrated that the responses to these drugs is incomplete [10], [11], suggesting the presence of an alternative pathway leading to AIA. Previous genetic studies have revealed significant associations between AIA and the polymorphisms of cysteinyl leukotriene receptor 1 [12], cysteinyl leukotriene receptor 2 [13], thromboxane A2 receptor [14], and prostaglandin E2 receptor subtype 2 [15]. The genetic polymorphisms of leukotriene C4 synthase (LTC4S) [16], [17], [18] and arachidonate 5-lipoxygenase (ALOX5) [19], [20] have been reported to be positively or negatively associated with AIA depending on populations and polymorphic loci. In addition to the genes on the arachidonate pathway, genetic variants on TNF [21], HLA [22], [23], [24], TBX21 [25], ACE [26], IL-10/TGF [27], and NLRP3 [28] genes in the pathways of immune response and inflammation are also associated with aspirin hypersensitivity in asthmatics. These data suggest that genetic variants of genes on other pathways may be more relevant to development of aspirin hypersensitivity in asthmatics than was previously thought.
Genome-wide association studies (GWAS) have recently emerged as a technology in predicting genetic variations across the genome that are associated with human diseases and clinical response to drug treatment [29], [30]. Recently, GWAS for asthma and related phenotypes have reported several susceptible genes including ORMDL3, PDE4D and IL1RL1 [31], [32], [33]. In the case of AIA, most of the recent genetic association analyses have been investigated through candidate gene approaches. Therefore, genome-wide association analysis and follow-up study of aspirin intolerance in asthmatics were performed in order to identify novel and comprehensive etiology of AIA.
Results
Characteristics of the Study Subjects
The clinical characteristics and comparisons of AIA, intermediate AIA (AIA-I), and ATA groups are summarized in Table 1. First of all, an aspirin-induced decline in FEV1 of −15% to 68% was observed in all study subjects. In both of genome-wide and follow-up studies, the % declines of FEV1 by aspirin provocation in AIA and AIA-I patients were significantly increased compared to those of ATA controls (P<0.0001), indicating that this study could reflect an association between genetic polymorphisms and aspirin hypersensitivity in asthmatics. In subjects of 2nd stage, the values of predicted FEV1 % and PC20 methacholine were significantly lower in AIA patients than those of ATA controls, whereas the total IgE level was higher in cases than in controls (P<0.01). A significant decrease of body mass index (BMI) in aspirin-intolerant asthmatics was also observed. In addition, the mean age of first medical examination was significantly lower in AIA than in the ATA group.
Genome-Wide Association Analyses
A total of 109,365 SNPs genotype assays were tested on the DNA samples of 80 AIA and 100 ATA subjects using the Illumina's Human-1 Genotyping BeadChip. The overall call rate of individual samples was over 98.0% after passing the genotype quality threshold, and a total of 4,515 SNPs (4.1% of 109,365 SNPs) failed to provide accurate genotype results. In further quality control analysis of the remaining 104,850 SNPs, monomorphic, X-chromosomal and Hardy Weinberg Equilibrium-departed (P<0.001) SNPs were additionally omitted; the final remaining 96,984 SNPs (88.7%) with an the average call rate of 99.9% were selected for analysis.
From results of the allelic association tests for each SNP in the co-dominant model, the 11 SNPs in the gene region that had the most significant association signals were screened for the second stage of follow-up study (Table 2).
Follow-up Study of Second Stage
For the follow-up study, the 11 genetic variants “within gene regions” with the lowest P value (P<0.0001) for risk of aspirin intolerance were selected. Then, 150 common SNPs in 11 candidate genes with minor allele frequency of over 0.05 based on Asian population were selected from the International HapMap Project (http://hapmap.ncbi.nlm.nih.gov/) and genotyped in 163 AIA subjects including AIA-I group, and 429 ATA subjects. In order to include the maximum number of patients with aspirin hypersensitivity, 80 AIA and 100 ATA subjects from the first stage were included in the second stage due to the rareness of the AIA condition. The case/control associations between genotype and aspirin intolerance in asthmatics were analyzed using multivariate logistic analyses adjusted for initial diagnosed age, sex, smoking status, atopy, and BMI. Polymorphisms of the CEP68 gene showed the most significant association with AIA compared to ATA controls in co-dominant model (Table 3, P = 6.0×10−5 to 4.0×10−5). Furthermore, a nonsynonymous SNP (rs7572857G>A) of CEP68 gene that replaces glycine with serine revealed to have the lowest P value among seven SNPs of the gene (OR = 2.63; 95% CI = 1.64–4.21; P = 6.0×10−5 for co-dominant model). This rs7572857G>A was also statistically significant in the dominant model (P = 0.0005), but not in recessive model (Table S1). In further association analysis excluding the GWAS subjects in the second set, the statistically significant association still remained although reduced significances were observed (Table S2).
This study also observed modest associations of other genes with AIA. In particular, the intronic rs4867084G>A of the PDZK3 gene (P = 0.006 for co-dominant model) and rs11060333C>T of the LRRC43 gene (P = 0.009 for co-dominant model) showed significant associations with AIA (Table S3). In addition, several intergenic SNPs showed significant associations with AIA in our GWAS results (P<0.0001; Table 4). Two SNPs, rs139719 and rs7744030, were identified as intergenic SNPs during the first stage of GWAS (build 126 of dbSNP database version), but they have been recently listed as intronic SNPs (build 130 of dbSNP). On the other hand, as for other nearby potential genes in the region where CEP68 is located, the RAB1A gene, despite the lack of reports showing direct relations with aspirin and/or respiratory disease, was also found to be associated with AIA by composing a strong LD with CEP68 (P<0.01; Figure 1).
(A) Associations of other potential genes nearby CEP68. A significant association of SNPs in RAB1A is found from the result of this GWAS (P<0.01). (B) LD Plot between CEP68 and RAB1A. LD coefficient (D') among SNPs of SLC1A4, RAB1A and CEP68 in Korean population.
Pair-wise comparisons among SNPs of CEP68 gene showed tight LDs (Figure 2, Table S4). After haplotypes were inferred using PHASE software, results from logistic analyses for associations between haplotypes and aspirin intolerance showed that haplotype CEP68_ht4 (T-G-A-A-A-C-G) of AIA was more frequent than that of ATA as a control group, showing that the CEP68_ht4 with the nonsynonymous rs7572857 “A” allele at the fourth position was significantly associated with AIA in both the co-dominant (OR = 2.63; 95% CI = 1.64–4.21; P = 6.0×10−5, Table 3) and dominant models (OR = 2.49; 95% CI = 1.49–4.15; P = 0.0005; Table S1) when compared to that of ATA as a control group.
(A) Physical map of CEP68 and its targeted SNPs. (B) Haplotypes of CEP68. Associations of haplotypes with frequency >0.05 are shown in Table 3. (C) LD coefficient (D') among CEP68 SNPs in a Korean population.
Effect of rs7572857 on the Decline of FEV1 by Aspirin Provocation
In regression analyses, the genotypes of several polymorphisms in the CEP68 gene were also significantly associated with the decline of FEV1 by aspirin provocation. Four SNPs, rs2302647, rs2252867, rs2723087, and rs6741255, showed a similar distribution pattern with about two-fold increase in the mean decline of FEV1 for the homozygous rare genotype compared to the homozygous common genotype (Table 5). In the case of rs7572857G>A, however, the nonsynonymous nucleotide substitution showed a different distribution with about a four-fold increase in the mean decline of FEV1 for the homozygote of the rare allele (Figure 3), indicating that rs7572857 could have a stronger influence on the higher decline of FEV1 by aspirin provocation than other variants. Despite our findings on the relationship between SNP genotypes and decline in FEV1 among AIA patients and then among ATA controls, results still showed that rs7572857 has an effect on AIA (P = 0.003; Table S5), but not on ATA.
Distribution presents the mean ± SE values of % decline of FEV1 by aspirin provocation in each genotype of the nonsynonymous SNP (rs7572857G>A) in the CEP68 gene. The homozygote of the rare allele (AA) showed a significant increase in decline of FEV1.
Discussion
To our knowledge, this is the first genome-wide association study for AIA. Our results identified CEP68 as a positive risk factor for the development of aspirin intolerance in asthmatics. Two stages of genome-wide analysis and a follow-up study were used to reduce the number of false-positive SNPs and retain statistical power. In the present study, logistic and regression analyses yielded significant associations between variations of the CEP68 gene and AIA, with the lowest P value at the nonsynonymous SNP rs7572857. In addition, with a tight LD among the SNPs of the CEP68 gene, the haplotype CEP68_ht4 (T-G-A-A-A-C-G) in AIA subjects was more frequent when compared to that of ATA controls and was significantly associated with AIA. In particular, the fourth positioned “A” that was derived from rs7572857G>A was found only in haplotype CEP68_ht4. On the other hand, although we also do not rule out a possible protective effect of rs10496123G>A due to its decreased risk for AIA (OR = 0.81; Table 3), it is considered that haplotype CEP68_ht4, which includes most minor alleles of the significantly associated SNPs and the major G allele of rs10496123G>A, may mainly affect the association with AIA.
Although the inclusion of intermediate AIA in AIA case group also showed significant associations with the SNPs and haplotypes of CEP68 gene compared to ATA controls, a more correct association for genetic polymorphisms and aspirin intolerance is suggested to the comparison between AIA subjects without intermediate AIA group and ATA controls (Table 3). On the other hand, when observing the effects of polymorphisms of CEP68 on the decline of FEV1 by aspirin provocation, the homozygous AA of rs7572857G>A variant showed a different distribution and a more significant increase of FEV1 decline in AIA patients than other variants. In addition, five AIA patients were found to have the homozygote of this rare allele (AA), with 34.80% FEV1 decline; whereas no ATA controls were observed to have the AA genotype (Table S5). This observation suggests that the change from glycine to serine could have an effect on the higher decline of FEV1 by aspirin provocation.
More recently, there have been debates concerning whether GWAS can successfully detect the variants that are associated with diseases. However, GWAS not only has made it possible to predict risk factors that are associated with diseases, but also has discovered additional variants that are associated with many diseases [34], [35]. Although dense genome-wide genotyping chips containing over 600 K and as many as 1 M SNPs have been recently developed, this study used the earlier chip with about 100 K SNPs due to the starting date of our research. Significantly, ours is the first GWAS to identify genetic factors that might influence the occurrence of aspirin intolerance in asthmatics. On the other hand, while GWAS have found many common variants that have modest effects on human common diseases, albeit some may have a key role, whereas the systematic identification of rare variants that are too rare to be detected in the GWAS and confer a substantial risk of the diseases is also required [36]. In terms of our study's limitations, statistical power indicated that our GWAS contained an insufficient sample size. Statistical powers of 72% in the first GWAS stage and 52% after selecting 150 SNPs in the first stage were estimated by QpowR, an interactive power calculator for two stage association studies [37]. In addition, after Bonferroni correction for multiple comparisons, only two SNPs, rs1053744 in SBF1 and rs828616 in DCBLD2, were significantly associated with AIA. However, this study applied a modest threshold (P<0.0001) to select top 11 candidate genes because we focused on the SNPs within a gene region. Moreover, the follow-up study for 11 promising candidate genes increased the power to calculate the association between the CEP68 gene and aspirin intolerance in asthmatics.
Since CEP68 gene (also referred to as KIAA0582) in the chromosome 2p14 has been recently discovered, its functions have not been fully understood, except for its role in centrosome cohesion and epidermal growth factor (EGF) signaling [38], [39]. In addition, given that smoke is a risk factor for asthma, CEP68 has been reported to be differentially expressed following exposure to environmental smoke [40]. However, the functions of the CEP68 protein are not yet fully understood. Thus, our findings provide a new insight into the relationship between CEP68 and aspirin intolerance in asthma. In particular, it is suggested that the Gly74Ser substitution in exon2 of the CEP68 gene could affect the polarity of the protein structure (Figure 4A), or affect the function of the CEP68 protein itself [41], [42]. Although in silico annotation showed that the nonsynonymous rs7572857G>A (Gly74Ser) appeared to be tolerable to diseases, this site was not highly conserved among mammals (Figure S1). Furthermore, in an additional functional site prediction, it was found that amino acids including 74Gly and its nearby positions could be targets for several kinases (Figure 4B). Since SH domains are especially unique sequences for several regulators, the change from glycine to serine might play an important role in correlations with other modulators. Therefore, phenotypic studies investigating the effects of minor A allele of rs7572857 on protein activity of CEP68 in relation to aspirin sensitivity are needed in the future.
The helical wheel projection of Gly74Ser is presented by the Expert Protein Analysis System (ExPASy). The change from 74Gly (A) to 74Ser (B) showed an increase in polarity due to its interactions with nearby polar amino acids. The solid surface energy and the degree of vapor/H2O distribution of glycine and serine were compared. (C) Motifs including 74Gly and nearby 74Gly for the functional site were found using a prediction program (http://elm.eu.org/index.html).
Linkage disequilibrium, as a nonrandom association, has been considered to be profoundly associated with recombination hot spots, quantitative characters, and complex diseases [43]. In addition, experimental evidence has shown an association between haplotype blocks and hot spots of recombination [44]. When compared with other populations from HapMap, the asthmatic subjects in the Korean population showed a stronger LD from the selected CEP68 SNPs than other populations, with a pattern that was closest to that of the Chinese (Figure S2). Given the significant association of haplotype CEP68_ht4 with AIA, the LD of CEP68 might facilitate mapping out genetic factors that function in aspirin metabolism and/or aspirin hypersensitivity in asthma. Therefore, further replications in AIA patients from other populations are required to establish whether or not these structural variations affect AIA.
Other previous reports have suggested associations between genetic polymorphisms and AIA. In one study, 370 SNPs of 63 candidate genes that are involved in the arachidonic acid metabolic cascade were studied for associations with AIA. The results demonstrated that SNPs in the promoter region of the prostaglandin E2 receptor subtype 2 gene were significantly associated with AIA [15]. Several SNPs in the promoter of LTC4S [19], [45] and ALOX5 [20] that encode key enzymes for CysLT synthesis also showed significant associations with AIA [13], [46]. In addition, three promoter polymorphisms in the cysteinyl leukotriene receptor 1 gene were associated with AIA in males [46]. Moreover, polymorphisms of genes that are involved in the pathways of immune response and inflammation are also associated with aspirin hypersensitivity in asthma [21], [25]. This study also confirmed that several other genes, which were associated previously with AIA, showed significant association signals with AIA pathogenesis, especially in the TNF, TGF, HLA-DPB1, ALOX5 and IL-10 genes (Table S6). Despite the many studies that have shown positive associations, a comprehensive association between genetic risks and aspirin hypersensitivity has not yet been discovered. Furthermore, constant revisions to find new risk alleles are needed [35], [47]. Hence, our findings on the first genome-wide and follow-up studies of AIA could be a breakthrough in the discovery of susceptible genes using a different approach.
In conclusion, this study confers a positive association between CEP68 and aspirin-intolerant asthma, suggesting that these findings would be useful for further genetic investigations of NSAIDs metabolism as well as other inflammatory diseases. However, this study does have several limitations. First, this study included only the top 11 genes for the follow-up study due to research budget constraints. Second, our second stage replication analysis was insufficient due to the inclusion of samples from the first stage. In addition, AIA-I patients who were also referred to as positive responders to aspirin hypersensitivity [48], [49] were included in the second round. Although the significances were reduced in the further association analysis without the GWAS subjects due to the high proportion of AIA-I in the second round, the statistically significant association of the SNPs including nonsynonymous rs7572857G>A with AIA still remained (Table S2). Therefore, further replication studies in a larger number of subjects are also needed. Finally, functional studies of the nonsynonymous SNP rs7572857 in the future could provide important insights into the genetic etiology of aspirin intolerance in asthmatics.
Materials and Methods
Study Subjects
The subjects were recruited from Soonchunhyang University, Chung Ang University, Chung Nam University, Chungbuk University, and Seoul National University in Korea in a span of about 6 years (between 2003 and 2008). All subjects were Korean. All patients were diagnosed by a physician and met the definition of asthma set forth in the Global Initiative for Asthma guidelines [50]. All patients had a history of dyspnea and wheezing during the previous 12 months, plus one of the following: (1) >15% increase in FEV1 or >12% increase plus 200 mL following inhalation of a short-acting bronchodilator, (2) <10 mg/mL PC20 methacholine, and (3) >20% increase in FEV1 following 2 weeks of treatment with inhaled steroids and long-acting bronchodilators. Twenty-four common inhalant allergens were used for a skin prick test [26]. Atopy was defined as having a wheal reaction over 3 mm in diameter. Total IgE was measured by the CAP system (Pharmacia Diagnostics, Uppsala, Sweden). The asthmatic patients had experienced no exacerbation of asthma and respiratory tract infection in the 6 weeks preceding the oral aspirin challenge (OAC). The OAC was performed with increasing doses of aspirin using methods slightly modified from those described previously [26], [49]. Changes in FEV1 were followed for 5 hours after the last aspirin challenge dose. Aspirin-induced bronchospasms, as reflected by rate (%) of FEV1 decline, were calculated as the pre-challenge FEV1 minus the post-challenge FEV1 divided by the pre-challenge FEV1. OAC reactions were categorized into three groups as follows: Subjects with 20% or greater decreases in FEV1 or 15% to 19% decreases in FEV1 with naso-ocular or cutaneous reactions were designated as the AIA group, those with 15% to 19% decreases in FEV1 or naso-ocular or cutaneous reactions only were the AIA-I group, and those with less than 15% decreases in FEV1 without naso-ocular or cutaneous reactions were the ATA group. The distribution of FEV1 decline rate of subjects is shown in Figure S3. All subjects provided informed consent, and the protocols were approved by the Institutional Review Board of each hospital.
The First Stage of Genome-Wide SNP Genotyping
About 750 ng of genomic DNA from 80 AIA cases and 100 ATA controls was used to genotype each sample on the Illumina's Human-1 Genotyping BeadChip (Illumina, San Diego, USA). Genotyping was performed according to Illumina's assay manual. Briefly, each sample was processed by whole-genome amplification, fragmentation, precipitation, and resuspension in an appropriate hybridization buffer. Denatured samples were hybridized on the prepared Human-1 Genotyping BeadChip for 16 h at 48°C. Then, the processed beadchip for the single-base extension reaction was stained and imaged on an Illumina Bead Array Reader. To convert fluorescent intensities into SNP genotypes, normalized bead intensity data obtained for each sample were loaded into the Beadstudio 3.0® software (Illumina). SNP clusters for genotype calling of all SNPs were examined using Beadstudio 3.0® software. The overall call rate for all SNPs was 98.0%.
The Second Stage of Follow-up Study
Eleven genes that include the GWAS-driven top SNPs, which are positioned not at intergenic region but at the 5′-untranslated region (UTR), exons, introns, and the 3′-UTR region, were selected for the follow-up study. For the next study, a total of 150 common SNPs within 11 candidate genes, including nearby upstream region (1.5 kb) of each gene, (8 SNPs for the SBF1 gene; 12 for DCBLD2; 15 for WDR21A; 34 for FILIP1; 2 for PDZK3; 15 for LRRC43; 18 for CIITA; 5 for DAF; 7 for ENPP5; 7 for CEP68; 27 for C6, respectively) were genotyped using a total of 163 AIA cases including AIA-I subjects and 429 ATA controls. The SNPs were scanned using BeadExpress® (Illumina, San Diego, USA).
Statistics
For genome-wide analysis, associations of genotype distributions between AIA cases and ATA controls were calculated by logistic analyses. In addition to the genetic homogeneity of the study subjects as described above, a parameter (λ) of genomic control was calculated in this study by dividing median χ2 statistics by 0.456, with an estimate of 1.056 [51]. Power calculations of the first GWAS stage and after selecting candidate SNPs in the first stage were estimated using QpowR program (https://www.msu.edu/~steibelj/JP_files/QpowR.pdf). In the association analysis of the polymorphisms in CEP68, we examined Lewontin's D' (|D'|) and the LD coefficient r2 between all pairs of biallelic loci [52]. Haploview v4.1 software downloaded from the Broad Institute (http://www.broadinstitute.org/mpg/haploview) was used to determine LD of CEP68 [53]. Haplotypes were first estimated using the PHASE software [54], and then computed by logistic analyses using the Statistical Analysis System (SAS) program. Subjects harboring missing genotypes were omitted in the analysis of individual SNPs and haplotypes. Comparisons of genotype distributions between AIA and ATA were carried out with logistic analyses adjusted for initial diagnosed age, sex, smoking status, atopy, and body mass index as multivariates using SAS. The common (C) alleles were used as the referent genotype to the heterozygote and homozygote of the rare (R) allele. The effective numbers of independent marker loci in each gene were calculated to correct for multiple testing using the software SNPSpD (http://genepi.qimr.edu.au/general/daleN/SNPSpD/), which is based on the spectral decomposition (SpD) of matrices of pair-wise LD between SNPs. Significant associations were represented by P value <0.05.
Supporting Information
Table S1.
Comparison of co-dominant, dominant and recessive models for SNPs in the CEP68 gene between AIA and ATA.
https://doi.org/10.1371/journal.pone.0013818.s001
(0.07 MB DOC)
Table S2.
Association analysis without GWAS samples in the second replication round analysis.
https://doi.org/10.1371/journal.pone.0013818.s002
(0.05 MB DOC)
Table S3.
Logistic analyses of 150 SNPs in Top 11 genes between AIA and ATA.
https://doi.org/10.1371/journal.pone.0013818.s003
(0.34 MB DOC)
Table S4.
LD coefficients (D' and r2) among CEP68 polymorphisms.
https://doi.org/10.1371/journal.pone.0013818.s004
(0.03 MB DOC)
Table S5.
Relation of the SNP genotypes in CEP68 to the decline on FEV1 in AIA and in ATA.
https://doi.org/10.1371/journal.pone.0013818.s005
(0.05 MB DOC)
Table S6.
List of other genes that were associated previously with AIA.
https://doi.org/10.1371/journal.pone.0013818.s006
(0.05 MB DOC)
Figure S1.
In silico annotation of nonsynonymous rs7572857G>A (Gly74Ser). Impact on protein function and conservation across species are predicted by the SNPs3D program (http://www.snps3d.org/). (A) The higher scores of entropy and Position Specific Scoring Matrix (PSSM) indicate more tolerable to diseases. (B) The 74th amino acid of CEP68 is not highly conserved among mammals.
https://doi.org/10.1371/journal.pone.0013818.s007
(0.17 MB DOC)
Figure S2.
Comparison of LD from other populations for the selected CEP68 SNPs. The LDs are constructed by seven CEP68 SNPs that are equivalent with this study from HapMap (http://hapmap.ncbi.nlm.nih.gov/index.html.en). CEU, Caucasian; CHB, Chinese, JPT, Japanese; YRI, African.
https://doi.org/10.1371/journal.pone.0013818.s008
(0.12 MB DOC)
Figure S3.
Distribution of FEV1 decline rate. Distributions are calculated from the decline rate of FEV1 by aspirin provocation and its number of subjects with 5% intervals. The distribution of AIA patients with FEV1 decline rate less than 15% is due to the responders to naso-ocular or cutaneous reactions.
https://doi.org/10.1371/journal.pone.0013818.s009
(0.03 MB DOC)
Acknowledgments
The DNA samples were generously provided by the Soonchunhyang University, Bucheon Hospital Biobank, a member of the National Biobank of Korea, supported by the Ministry of Health, Welfare and Family Affairs, Republic of Korea.
Author Contributions
Conceived and designed the experiments: CSP HDS. Performed the experiments: HSC JSB. Analyzed the data: JHK BLP HSC JSB HDS. Contributed reagents/materials/analysis tools: JSP ASJ STU JSC YHK MKK ISC SHC BWC CSP. Wrote the paper: JHK BLP CSP HDS.
References
- 1. Widal F, Abrami P, Lermoyez J (1993) Anaphylaxie et idiosyncrasie. 1992 [Anaphylaxis and idiosyncrasy. 1992]. Allergy Proc 14: 373–376.
- 2. Babu KS, Salvi SS (2000) Aspirin and asthma. Chest 118: 1470–1476.
- 3. Szczeklik A, Stevenson DD (2003) Aspirin-induced asthma: advances in pathogenesis, diagnosis, and management. J Allergy Clin Immunol 111: 913–921.
- 4. Kowalski ML (2007) Aspirin-sensitive rhinosinusitis and asthma. Clin Allergy Immunol 19: 147–175.
- 5. Palikhe NS, Kim SH, Park HS (2008) What do we know about the genetics of aspirin intolerance? J Clin Pharm Ther 33: 465–472.
- 6. Tantisira KG, Drazen JM (2009) Genetics and pharmacogenetics of the leukotriene pathway. J Allergy Clin Immunol 124: 422–427.
- 7. Duroudier NP, Tulah AS, Sayers I (2009) Leukotriene pathway genetics and pharmacogenetics in allergy. Allergy 64: 823–839.
- 8. Drazen JM, Israel E, O'Byrne PM (1999) Treatment of asthma with drugs modifying the leukotriene pathway. N Engl J Med 340: 197–206.
- 9. Leff AR (2001) Regulation of leukotrienes in the management of asthma: biology and clinical therapy. Annu Rev Med 52: 1–14.
- 10. Park JS, Jang AS, Park SW, Lee YM, Uh ST, et al. (2010) Protection of leukotriene receptor antagonist against aspirin-induced bronchospasm in asthmatics. Allergy Asthma Immunol Res 2: 48–54.
- 11. Dahlen SE, Malmstrom K, Nizankowska E, Dahlen B, Kuna P, et al. (2002) Improvement of aspirin-intolerant asthma by montelukast, a leukotriene antagonist: a randomized, double-blind, placebo-controlled trial. Am J Respir Crit Care Med 165: 9–14.
- 12. Kim SH, Ye YM, Hur GY, Lee SK, Sampson AP, et al. (2007) CysLTR1 promoter polymorphism and requirement for leukotriene receptor antagonist in aspirin-intolerant asthma patients. Pharmacogenomics 8: 1143–1150.
- 13. Park JS, Chang HS, Park CS, Lee JH, Lee YM, et al. (2005) Association analysis of cysteinyl-leukotriene receptor 2 (CYSLTR2) polymorphisms with aspirin intolerance in asthmatics. Pharmacogenet Genomics 15: 483–492.
- 14. Kim SH, Choi JH, Park HS, Holloway JW, Lee SK, et al. (2005) Association of thromboxane A2 receptor gene polymorphism with the phenotype of acetyl salicylic acid-intolerant asthma. Clin Exp Allergy 35: 585–590.
- 15. Jinnai N, Sakagami T, Sekigawa T, Kakihara M, Nakajima T, et al. (2004) Polymorphisms in the prostaglandin E2 receptor subtype 2 gene confer susceptibility to aspirin-intolerant asthma: a candidate gene approach. Hum Mol Genet 13: 3203–3217.
- 16. Sanak M, Simon HU, Szczeklik A (1997) Leukotriene C4 synthase promoter polymorphism and risk of aspirin-induced asthma. Lancet 350: 1599–1600.
- 17. Van Sambeek R, Stevenson DD, Baldasaro M, Lam BK, Zhao J, et al. (2000) 5′ flanking region polymorphism of the gene encoding leukotriene C4 synthase does not correlate with the aspirin-intolerant asthma phenotype in the United States. J Allergy Clin Immunol 106: 72–76.
- 18. Sanak M, Szczeklik A (2001) Leukotriene C4 synthase polymorphism and aspirin-induced asthma. J Allergy Clin Immunol 107: 561–562.
- 19. Choi JH, Park HS, Oh HB, Lee JH, Suh YJ, et al. (2004) Leukotriene-related gene polymorphisms in ASA-intolerant asthma: an association with a haplotype of 5-lipoxygenase. Hum Genet 114: 337–344.
- 20. Kim SH, Bae JS, Suh CH, Nahm DH, Holloway JW, et al. (2005) Polymorphism of tandem repeat in promoter of 5-lipoxygenase in ASA-intolerant asthma: a positive association with airway hyperresponsiveness. Allergy 60: 760–765.
- 21. Kim SH, Ye YM, Lee SK, Choi JH, Holloway JW, et al. (2006) Association of TNF-alpha genetic polymorphism with HLA DPB1*0301. Clin Exp Allergy 36: 1247–1253.
- 22. 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.
- 23. 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.
- 24. 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.
- 25. Akahoshi M, Obara K, Hirota T, Matsuda A, Hasegawa K, et al. (2005) Functional promoter polymorphism in the TBX21 gene associated with aspirin-induced asthma. Hum Genet 117: 16–26.
- 26. Kim TH, Chang HS, Park SM, Nam BY, Park JS, et al. (2008) Association of angiotensin I-converting enzyme gene polymorphisms with aspirin intolerance in asthmatics. Clin Exp Allergy 38: 1727–1737.
- 27. Kim SH, Yang EM, Lee HN, Cho BY, Ye YM, et al. (2009) Combined effect of IL-10 and TGF-beta1 promoter polymorphisms as a risk factor for aspirin-intolerant asthma and rhinosinusitis. Allergy 64: 1221–1225.
- 28. Hitomi Y, Ebisawa M, Tomikawa M, Imai T, Komata T, et al. (2009) Associations of functional NLRP3 polymorphisms with susceptibility to food-induced anaphylaxis and aspirin-induced asthma. J Allergy Clin Immunol 124: 779–785 e776.
- 29. Hirschhorn JN, Daly MJ (2005) Genome-wide association studies for common diseases and complex traits. Nat Rev Genet 6: 95–108.
- 30. Roses AD (2008) Pharmacogenetics in drug discovery and development: a translational perspective. Nat Rev Drug Discov 7: 807–817.
- 31. Moffatt MF, Kabesch M, Liang L, Dixon AL, Strachan D, et al. (2007) Genetic variants regulating ORMDL3 expression contribute to the risk of childhood asthma. Nature 448: 470–473.
- 32. Himes BE, Hunninghake GM, Baurley JW, Rafaels NM, Sleiman P, et al. (2009) Genome-wide association analysis identifies PDE4D as an asthma-susceptibility gene. Am J Hum Genet 84: 581–593.
- 33. Gudbjartsson DF, Bjornsdottir US, Halapi E, Helgadottir A, Sulem P, et al. (2009) Sequence variants affecting eosinophil numbers associate with asthma and myocardial infarction. Nat Genet 41: 342–347.
- 34. Goldstein DB (2009) Common genetic variation and human traits. N Engl J Med 360: 1696–1698.
- 35. Kraft P, Hunter DJ (2009) Genetic risk prediction–are we there yet? N Engl J Med 360: 1701–1703.
- 36. Cirulli ET, Goldstein DB (2010) Uncovering the roles of rare variants in common disease through whole-genome sequencing. Nat Rev Genet 11: 415–425.
- 37. Zabaneh D, Balding DJ (2010) A genome-wide association study of the metabolic syndrome in Indian asian men. PLoS One(In Press).
- 38. Graser S, Stierhof YD, Nigg EA (2007) Cep68 and Cep215 (Cdk5rap2) are required for centrosome cohesion. J Cell Sci 120: 4321–4331.
- 39. Chen Y, Low TY, Choong LY, Ray RS, Tan YL, et al. (2007) Phosphoproteomics identified Endofin, DCBLD2, and KIAA0582 as novel tyrosine phosphorylation targets of EGF signaling and Iressa in human cancer cells. Proteomics 7: 2384–2397.
- 40. Rouse RL, Boudreaux MJ, Penn AL (2007) In utero environmental tobacco smoke exposure alters gene expression in lungs of adult BALB/c mice. Environ Health Perspect 115: 1757–1766.
- 41. Weinshilboum R, Wang L (2004) Pharmacogenomics: bench to bedside. Nat Rev Drug Discov 3: 739–748.
- 42. Maheshwari R, Dhathathreyan A (2004) Investigation of surface properties of amino acids: polarity scale for amino acids as a means to predict surface exposed residues in films of proteins. J Colloid Interface Sci 277: 79–83.
- 43. Slatkin M (2008) Linkage disequilibrium–understanding the evolutionary past and mapping the medical future. Nat Rev Genet 9: 477–485.
- 44. Jeffreys AJ, Kauppi L, Neumann R (2001) Intensely punctate meiotic recombination in the class II region of the major histocompatibility complex. Nat Genet 29: 217–222.
- 45. Kawagishi Y, Mita H, Taniguchi M, Maruyama M, Oosaki R, et al. (2002) Leukotriene C4 synthase promoter polymorphism in Japanese patients with aspirin-induced asthma. J Allergy Clin Immunol 109: 936–942.
- 46. Kim SH, Oh JM, Kim YS, Palmer LJ, Suh CH, et al. (2006) Cysteinyl leukotriene receptor 1 promoter polymorphism is associated with aspirin-intolerant asthma in males. Clin Exp Allergy 36: 433–439.
- 47. Altshuler D, Daly MJ, Lander ES (2008) Genetic mapping in human disease. Science 322: 881–888.
- 48.
Stevenson DD (1997) Aspirin-sensitive Asthma. In: Barnes PJ, Grunstein MM, Leff AR, Woolcock AJ, editors. Asthma. New York: Lippincott-Raven. pp. 2025–2040.
- 49. Cormican LJ, Farooque S, Altmann DR, Lee TH (2005) Improvements in an oral aspirin challenge protocol for the diagnosis of aspirin hypersensitivity. Clin Exp Allergy 35: 717–722.
- 50.
Global Initiative for Asthma (GINA) Global strategy for asthma management and prevention. NHLBI/WHO workshop reportBethesda, MD: National Institutes of Health, National Heart, Lung, and Blood Institute, 1995, NIH publication, 1995. pp. 95–3659.
- 51. Devlin B, Roeder K, Wasserman L (2001) Genomic control, a new approach to genetic-based association studies. Theor Popul Biol 60: 155–166.
- 52. Hedrick PW (1987) Gametic disequilibrium measures: proceed with caution. Genetics 117: 331–341.
- 53. Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21: 263–265.
- 54. Stephens M, Smith NJ, Donnelly P (2001) A new statistical method for haplotype reconstruction from population data. Am J Hum Genet 68: 978–989.