Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2012; 7(4): e35334.
Published online Apr 24, 2012. doi:  10.1371/journal.pone.0035334
PMCID: PMC3335873
Association Analysis of Single Nucleotide Polymorphisms at Five Loci: Comparison between Atopic Dermatitis and Asthma in the Chinese Han Population
Hua-Yang Tang,#1,2,3 Xian-Fa Tang,#1,2,3 Xian-Bo Zuo,2,3 Jin-Ping Gao,1,2,3 Yu-Jun Sheng,1,2,3 Yang Li,1,2,3 Fu-Sheng Zhou,2,3 Xian-Yong Yin,1,2,3 Feng-Li Xiao,1,2,3 Wei-Dong Du,2,3 Sen Yang,1,2,3 Liang-Dan Sun,1,2,3* and Xue-Jun Zhang1,2,3*
1Institute of Dermatology and Department of Dermatology, No.1 Hospital, Anhui Medical University, Hefei, Anhui, China
2Key Laboratory of Dermatology, Anhui Medical University, Ministry of Education, China, Hefei, Anhui, China
3State Key Laboratory Incubation Base of Dermatology, Anhui Medical University, Hefei, Anhui, China
Maurice A. M. van Steensel, Editor
Maastricht University Medical Center, The Netherlands
#Contributed equally.
* E-mail: ahmusld/at/ (LDS); ayzxj/at/ (XJZ)
Conceived and designed the experiments: HYT XFT LDS XJZ. Performed the experiments: FSZ YL JPG YJS. Analyzed the data: XBZ XYY WDD. Contributed reagents/materials/analysis tools: FLX SY. Wrote the paper: HYT.
Received June 22, 2011; Accepted March 14, 2012.
Atopic diseases, such as atopic dermatitis (AD) and asthma, are closely related to clinical phenotypes with hypersensitivity, and often share some similar genetic and pathogenic bases. Our recent GWAS identified three susceptibility gene/loci FLG (rs11204971 and rs3126085), 5q22.1 (rs10067777, rs7701890, rs13360927 and rs13361382) and 20q13.33 (rs6010620) to AD. The effect of these AD associated polymorphisms in asthma is so far unknown. To investigate whether AD relevant genetic variants is identical to asthma and reveal the differences in genetic factors between AD and asthma in Chinese Han population, seven AD associated single nucleotide polymorphisms (SNPs) as well as 3 other SNPs (rs7936562 and rs7124842 at 11q13.5 and rs4982958 at 14q11.2) from our previous AD GWAS were genotyped in 463 asthma patients and 985 controls using Sequenom MassArray system. We found rs4982958 at 14q11.2 was significantly associated with asthma (P = 3.04×10−4, OR = 0.73). We also detected one significant risk haplotype GGGA from the 4 SNPs (rs10067777, rs7701890, rs13360927 and rs13361382) at 5q22.1 in AD cases (Pcorrection = 3.60×10−10, OR = 1.26), and the haplotype was suggestive of risk in asthma cases in this study (P = 0.014, Pcorrection = 0.084, OR = 1.38). These SNPs (rs11204971, rs3126085, rs7936562, rs712484 and rs6010620) at AD susceptibility genes/loci FLG, 11q13.5 and 20q13.33 were not associated with asthma in this study. Our results further comfirmed that 14q11.2 was an important candidate locus for asthma and demonstrated that 5q22.1 might be shared by AD and asthma in Chinese Han population.
Atopic diseases, such as atopic dermatitis (AD), asthma and allergic rhinitis (AR), are closely related to clinical phenotypes that are often observed with hypersensitivity. The mechanisms leading to atopic diseases development remain yet unknown. It is supposed that the cases would attribute to the interacted effect of multiple genetic variations and environmental factors on their pathophysiology. Asthma is documented being associated with atopic dermatitis [1]. Approximately half of AD patients, particularly in the severe AD, will develop into asthma [2].
Previous linkage study, candidate gene and GWAS have identified several susceptibility loci/genes predisposing to AD and asthma. Some common genetic factors such as FLG, ADAM33, and GSDML/ORMDL3 had been observed in patients with AD and asthma [3]. Other genetic investigations illuminated several important candidate genes in Chromosome 14q11.2, which could contribute to the genetic predisposition to specific allergic reactions characterized by asthma and high level of serum IgE [4], [5]. In one previous GWAS of AD in the European populations, 11q13.5 (C11orf30) was highlighted as a susceptibility locus of AD [6], which was also pointed to contain a concurrent involvement of the common genetic factors for AD and asthma [7]. In one recent AD GWAS, we identified two novel susceptibility loci 5q22.1 (rs10067777, rs7701890, rs13360927 and rs13361382), 20q13.33 (rs6010620) and confirmed susceptibility gene FLG (rs11204971 and rs3126085) of AD [8]. The genetic effect of these AD associated factors in asthma is unknown.
In this study, we further analyze the 7 SNPs described above at FLG, 5q22.1, 20q13.33, as well as 3 other SNPs (rs7936562 and rs7124842 at 11q13.5, rs4982958 at 14q11.2) to identify genetic variants for asthma and reveal the common and different genetic factors for AD and asthma in the Chinese Han population.
A total of 463 patients with asthma and 985 controls were enrolled in this study (Table 1). All samples were unrelated and Chinese Han in origin. Patients with asthma were recruited from No. 1 Hospital of Anhui Medical University in China according to the American Thoracic Society (ATS) guidelines [9]: First, the presence of at least two of the following symptoms was examined: recurrent cough, wheezing or dyspnea. Second, we examined for increased airway responsiveness to methacholine or the presence of reversible airflow limitation; the latter refers to 15% variability in the forced expiratory volume in 1 second (FEV1), or in the peak expiratory flow rate with or without an inhaled short-acting β2-agonist. Third, there was the absence of any other pulmonary diseases. Clinical information was collected from the affected individuals through a full clinical checkup by physician specialists. Additional demographic information was collected from both cases and controls through a structured questionnaire. All controls were clinically assessed to be without asthma, other atopic diseases, family history of atopic diseases (including first-, second- and third-degree relatives) or systemic disorders. The study was approved by the Institutional Ethics Committee of Anhui Medical University and was conducted according to Declaration of Helsinki principles. Collection of blood samples and clinical information from cases and controls was undertaken after written informed consent was obtained from all participants. DNA was extracted from peripheral blood lymphocytes using QIAamp DNA Blood kit (Qiagen, Valencia, CA, USA) according to the manufacturer's instructions. The extracted genomic DNAs were analyzed by agarose gel electrophoresis, quantified by spectrophotometer, and stored at −80°C until used.
Table 1
Table 1
Summary information of asthma patients and health controls.
SNP selection and genotyping
We selected 10 SNPs based on our previous AD GWAS and genotyped them in 463 asthma patients. Seven of the 10 SNPs were associated with AD (rs11204971 and rs3126085 in FLG, rs10067777, rs7701890, rs13360927 and rs13361382 at 5q22.1, and rs6010620 at 20q13.33). Additional 2 top SNPs within 11q13.5 (rs7936562 and rs7124842) and one top SNP within 14q11.2 (rs4982958) in discovery samples of our AD GWAS were also selected. The annotated nominal gene was taken from the Illumina SNP database. SNPs were genotyped using the Sequenom MassArray system at State Key Laboratory Incubation Base of Dermatology, Ministry of National Science and Technology, Hefei, Anhui, China. Approximately 15 ng of genomic DNA was used to genotype each sample. Locus-specific PCR and detection primers were designed using the MassARRAY Assay Design 3.0 software (Sequenom). The DNA samples were amplified by multiplex PCR reactions, and the PCR products were then used for locus-specific single-base extension reactions. The resulting products were desalted and transferred to a 384-element SpectroCHIP array. Allele detection was performed using MALDI-TOF MS. The mass spectrograms were analyzed by the MassARRAY Typer software (Sequenom).
Statistical analysis
The SNPs were analyzed for association with the disease by means of comparison of the minor allele frequency in the cases and the controls using PLINK 1.07 software [10]. The level of associated significance was assigned at P values of less than 0.005 after Bonferroni Multiple Testing correction for asthma case-control analysis, and less than 5.0×10−8 for AD GWAS. In addition to the allelic test of association, the genetic models (dominant model, recessive model and additive model) were calculated for the associated SNPs. For the SNPs at 5q22.1 (TMEM232/SLC25A46), haplotype analyses were performed to calculate haplotype frequencies and significance of associations for asthma case-control study by using Haploview v4.2 [11]. The extent of genetic heterogeneity was assessed by using Cochran's Q test and the I2 index [12] between AD and asthma groups. All the SNPs passed the quality control in terms of with a call rate >95%, Hardy–Weinberg equilibrium (HWE) (P>0.01) in the controls. The genetic statistical power for all genotyped SNPs was estimated using CaTS-Power Calculator software [13].
SNP rs4982958 at 14q11.2 was significantly associated with asthma in this study (P = 3.04×10−4, OR = 0.73). Risk allele C of the SNP rs4982958 represented the major allele and allele T was protective one in Chinese Han population for asthma (Table 2). The genotype TT appeared significantly less frequent in asthma cases (5.5%) than in controls (12.9%), higher frequency of the genotype CC was observed in cases (46.1%) than in controls (39.8%). The difference in distribution of rs4982958 genotype between the asthma cases and controls was statistically significant (P = 5.96×10−5) (Table 3). Those SNPs locating within AD susceptibility gene/loci FLG, 5q22.1, 11q13.5 and 20q13.33 did not reach the threshold significant association for asthma (P>0.005, after Bonferroni Multiple Testing correction) (Table 2). For the SNPs in FLG (rs11204971 and rs3126085) and 14q11.2 (rs4982958), we found significant evidence for genetic heterogeneity between AD and asthma groups in this study (Phet≤0.0031, I2≥88.6) (Table 2).
Table 2
Table 2
Summary of association results of 10 SNPs in five loci/genes and heterogeneity test between AD and asthma groups.
Table 3
Table 3
Distribution of genotypes and genetic model analysis for rs4982958 in asthma patients and controls.
Six haplotypes in combination with the 4 SNPs (rs10067777, rs7701890, rs13360927 and rs13361382) were established at 5q22.1 in asthma samples, and 3 haplotypes in AD samples. The integrated haplotype GGGA showed a significant association with AD (Pcorrection = 3.60×10−10, OR = 1.26) and weakly significant in asthma cases and controls (P = 0.014, OR = 1.38) (Table 4), which had no genetic heterogeneity between AD and asthma (Phet = 0.590, I2 = 0). After Bonferroni Multiple Testing correction, haplotype GGGA was no more significant for asthma (Pcorrection = 0.084) (Table 4).
Table 4
Table 4
Haplotype analysis of 4 SNPs within 5q22.1 in AD/asthma patients and controls.
AD and asthma are complex genetic diseases, which would result from integrated effect of multiple genetic and environmental factors on pathophysiology process of the entities. Some common genetic factors being shared in the patients with AD and asthma have been identified [3]. Here we investigated the association of 10 SNPs within FLG, 5q22.1, 11q13.5, 14q11.2 and 20q13.33 for asthma in the Chinese Han population, to evaluate the relationship between AD and asthma in these loci/genes.
At the locus 14q11.2, the five genes: KIAA0323, SDR39U1, CMA1, CTSG and GZMH were located within a single 220 kb LD block surrounding the tag biomarker rs4982958 which was associated with asthma in this study (Figure 1a). Previous studies selected mast cell chymase gene (CMA1) at 14q11.2 as a candidate gene for analysis of AD, asthma and other allergic phenotypes [14]. Mast cell tryptase and chymase are key players in mediating chronic inflammation, as observed in asthmatics of the lungs. Two polymorphisms in CMA1, including promoter polymorphism rs1800875 (−1903G/A) [15][18] and (TG)n(GA)m repeat polymorphism [19][21], were identified being associated with the allergic phenotypes in different populations. However, other studies failed to replicate role of −1903G/A (rs1800875) in AD or asthma [18], [22][24]. Inconsistence in these studies probably was attributable to different populations, variation on phenotype definitions, ethnicity and environmental exposure.
Figure 1
Figure 1
Regional LD plots for associated SNPs.
We observed that the association of rs4982958 with asthma under additive model (P = 1.75×10−5) was of more significance than that under dominant model (P = 2.50×10−2) and recessive model (P = 3.05×10−5) by using a logistic regression analysis (Table 3). Clearly, risk allele C of rs4982958 was the major one and allele T was of protection effect. Compared with C/C, the genotype of homozygote T/T (ORHom = 0.37) showed smaller OR than the heterozygote T/C (ORHet = 0.89) (Table 3). This result may provide more information for futher study on this locus. Rs4982958 located at 10.4 kb upstream of CMA1. Linkage disequilibrium (LD) test showed that the two SNPs (rs4982958 and rs1800875) were moderately correlated with each other (D' = 1.0, r2 = 0.5) in CHB and JPT from 1000 Genomes Project data. Imputation of promoter polymorphism rs1800875 of CMA1 was performed in our AD GWAS samples and then tested for association (P = 1.11×10−2, OR = 0.85). In AD GWAS, rs4982958 at 14q11.2 was replicated in expanded AD samples, we found a weak evidence of association in the Chinese Han population (P = 1.22×10−2, OR = 0.94), but did not reach genome-wide significance (P<5.0×10−8) [8]. Statistical tests showed that the odds ratio of rs4982958 in asthma samples (OR = 0.73) was smaller than in AD (OR = 0.94) and the existence of significant genetic heterogeneity between AD and asthma (Phet = 0.0031, I2 = 88.60) in Chinese Han population.
Our AD GWAS had identified 5q22.1 as a susceptibility locus for AD [8]. In this study, further haplotype analysis of the 4 SNPs (rs10067777, rs7701890, rs13360927 and rs13361382) within 5q22.1 was performed in AD and asthma groups, respectively. Only two haplotypes (AAAG and GGGA) had allele frequencies >5% (Table 4). The risk haplotype GGGA had a consistent association evidence for AD (Pcorrection = 3.60×10−10, OR = 1.26) as single SNP did (Table 4). It further supported the presence of underlying causal variants within this locus for AD. None of the 4 SNPs within 5q22.1 was associated with asthma in such a small samples size and low statistical power (1%~2%) in this study. The risk allele of these SNPs at this locus was the same for AD and asthma. In asthma group, haplotype GGGA was a suggestive risk factor in this study (P = 0.014), although it was barely below significance (P = 0.084) after multiple testing correction. We observed that the odds ratio of haplotype GGGA (OR = 1.38) was much higher than single SNP in asthma cases-controls samples (1.04≤OR≤1.10). Further more, heterogeneity test at 5q22.1 showed that no genetic heterogeneity was found for each SNPs (0.2313≤Phet≤0.5265, 0≤I2≤30.2) (Table 2) and risk haplotype GGGA (Phet = 0.590, I2 = 0) (Table 4) between AD and asthma groups, which means 5q22.1 might have similar effect on AD and asthma. The four SNPs (rs10067777, rs7701890, rs13360927 and rs13361382) located within a 600 kb LD block harboring two genes TMEM232 (FLJ43080) and SLC25A46. Little is known about the function of the two genes. At 5q22.1, TSLP and WDR36 were proposed to be candidate genes for asthma in Caucasian population [25][27]. Nevertheless, we found haplotype GGGA fall into different LD blocks with TSLP and WDR36, which were separated by a recombination hot spot (Figure 1b). Fine-mapping or direct sequencing would be powerful to further search for any potential and associated susceptibility genes at 5q22.1 dedicated to both AD and asthma.
The 5 SNPs in FLG, 11q13.5 and 20q13.33 did not reached the statistical significance in this study, which might either result from low statistical power (1%~15%) (Table 2) or be attributed to 10.6% of these asthma cases with previous history of AD (Table 1). Comprehensive sequencing of entire FLG in asthma patients and larger scale of asthma samples or increased related tagSNPs for 11q13.5 and 20q13.33 would reliably enable to explore the relationship between AD and asthma in these loci/gene.
In conclusion, we compared the identical SNPs from the asthma cohort and our previous AD GWAS. We comfirmed 14q11.2 was an important candidate locus for asthma, and also demonstrated that 5q22.1 might be shared by AD and asthma in Chinese Han population. However, the evidences for the associations will be required to validate in diverse populations and then gain a better understanding on AD and asthma pathogenesis.
We are most grateful to all the members who have so willingly participated in this study, which made this study possible.
Competing Interests: The authors have declared that no competing interests exist.
Funding: This study was funded by General Program of National Natural Science Foundation of China (81071284, 81172838, 31171224, and 31000528), Science and Technological Fund of Anhui Province for Outstanding Youth (1108085J10) and the National Basic Research Program of China (973 Program, 2012CB722404). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
1. Cookson WO, Moffatt MF. The genetics of atopic dermatitis. Curr Opin Allergy Clin Immunol. 2002;2:383–387. [PubMed]
2. Spergel JM. From atopic dermatitis to asthma: the atopic march. Ann Allergy Asthma Immunol. 2010;105:99–106; quiz 107–109, 117. [PubMed]
3. Holloway JW, Yang IA, Holgate ST. Genetics of allergic disease. J Allergy Clin Immunol. 2010;125:S81–94. [PubMed]
4. Moffatt MF, Hill MR, Cornelis F, Schou C, Faux JA, et al. Genetic linkage of T-cell receptor alpha/delta complex to specific IgE responses. Lancet. 1994;343:1597–1600. [PubMed]
5. CSGA. A genome-wide search for asthma susceptibility loci in ethnically diverse populations. The Collaborative Study on the Genetics of Asthma (CSGA). Nat Genet. 1997;15:389–392. [PubMed]
6. Esparza-Gordillo J, Weidinger S, Folster-Holst R, Bauerfeind A, Ruschendorf F, et al. A common variant on chromosome 11q13 is associated with atopic dermatitis. Nat Genet. 2009;41:596–601. [PubMed]
7. Marenholz I, Bauerfeind A, Esparza-Gordillo J, Kerscher T, Granell R, et al. The eczema risk variant on chromosome 11q13 (rs7927894) in the population-based ALSPAC cohort: a novel susceptibility factor for asthma and hay fever. Hum Mol Genet. 2011;20:2443–2449. [PubMed]
8. Sun LD, Xiao FL, Li Y, Zhou WM, Tang HY, et al. Genome-wide association study identifies two new susceptibility loci for atopic dermatitis in the Chinese Han population. Nat Genet. 2011;43:690–694. [PubMed]
9. American Thoracic Society. Standards for the diagnosis and care of patients with chronic obstructive pulmonary disease (COPD) and asthma. Am Rev Respir Dis. 1987;136:225–244. [PubMed]
10. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–575. [PubMed]
11. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–265. [PubMed]
12. Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. British Medical Journal. 2003;327:557–560. [PMC free article] [PubMed]
13. Skol AD, Scott LJ, Abecasis GR, Boehnke M. Joint analysis is more efficient than replication-based analysis for two-stage genome-wide association studies. Nat Genet. 2006;38:209–13. [PubMed]
14. Hart PH. Regulation of the inflammatory response in asthma by mast cell products. Immunol Cell Biol. 2001;79:149–153. [PubMed]
15. Mao XQ, Shirakawa T, Yoshikawa T, Yoshikawa K, Kawai M, et al. Association between genetic variants of mast-cell chymase and eczema. Lancet. 1996;348:581–583. [PubMed]
16. Mao XQ, Shirakawa T, Enomoto T, Shimazu S, Dake Y, et al. Association between variants of mast cell chymase gene and serum IgE levels in eczema. Hum Hered. 1998;48:38–41. [PubMed]
17. Tanaka K, Sugiura H, Uehara M, Sato H, Hashimoto-Tamaoki T, et al. Association between mast cell chymase genotype and atopic eczema: comparison between patients with atopic eczema alone and those with atopic eczema and atopic respiratory disease. Clin Exp Allergy. 1999;29:800–803. [PubMed]
18. Iwanaga T, McEuen A, Walls AF, Clough JB, Keith TP, et al. Polymorphism of the mast cell chymase gene (CMA1) promoter region: lack of association with asthma but association with serum total immunoglobulin E levels in adult atopic dermatitis. Clin Exp Allergy. 2004;34:1037–1042. [PubMed]
19. Sharma S, Rajan UM, Kumar A, Soni A, Ghosh B. A novel (TG)n(GA)m repeat polymorphism 254 bp downstream of the mast cell chymase (CMA1) gene is associated with atopic asthma and total serum IgE levels. J Hum Genet. 2005;50:276–282. [PubMed]
20. Hossny EM, Amr NH, Elsayed SB, Nasr RA, Ibraheim EM. Association of polymorphisms in the mast cell chymase gene promoter region (−1903 g/A) and (TG)n(GA)m repeat downstream of the gene with bronchial asthma in children. J Investig Allergol Clin Immunol. 2008;18:376–381. [PubMed]
21. Hersberger M, Thun GA, Imboden M, Brandstatter A, Waechter V, et al. Association of STR polymorphisms in CMA1 and IL-4 with asthma and atopy: The SAPALDIA Cohort. Hum Immunol 2010 [PubMed]
22. Kawashima T, Noguchi E, Arinami T, Kobayashi K, Otsuka F, et al. No evidence for an association between a variant of the mast cell chymase gene and atopic dermatitis based on case-control and haplotype-relative-risk analyses. Hum Hered. 1998;48:271–274. [PubMed]
23. Forrest S, Dunn K, Elliott K, Fitzpatrick E, Fullerton J, et al. Identifying genes predisposing to atopic eczema. J Allergy Clin Immunol. 1999;104:1066–1070. [PubMed]
24. Pascale E, Tarani L, Meglio P, Businco L, Battiloro E, et al. Absence of association between a variant of the mast cell chymase gene and atopic dermatitis in an Italian population. Hum Hered. 2001;51:177–179. [PubMed]
25. He JQ, Hallstrand TS, Knight D, Chan-Yeung M, Sandford A, et al. A thymic stromal lymphopoietin gene variant is associated with asthma and airway hyperresponsiveness. J Allergy Clin Immunol. 2009;124:222–229. [PubMed]
26. Hunninghake GM, Soto-Quiros ME, Avila L, Kim HP, Lasky-Su J, et al. TSLP polymorphisms are associated with asthma in a sex-specific fashion. Allergy. 2010;65:1566–1575. [PMC free article] [PubMed]
27. Gudbjartsson DF, Bjornsdottir US, Halapi E, Helgadottir A, Sulem P, et al. Sequence variants affecting eosinophil numbers associate with asthma and myocardial infarction. Nat Genet. 2009;41:342–347. [PubMed]
Articles from PLoS ONE are provided here courtesy of
Public Library of Science