Skip to main content
Advertisement
Browse Subject Areas
?

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

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Unifying Candidate Gene and GWAS Approaches in Asthma

  • Sven Michel,

    Affiliation Center for Pediatrics, Clinic for Pediatric Pneumology, Allergology and Neonatology, Hannover Medical School, Hannover, Germany

  • Liming Liang,

    Affiliation Center for Statistical Genetics, Department of Biostatistics, University of Michigan, Ann Arbor, Michigan, United States of America

  • Martin Depner,

    Affiliation University Children's Hospital, Ludwig Maximilian's University, Munich, Germany

  • Norman Klopp,

    Affiliation Institute of Epidemiology, Helmholtz Centre Munich, Neuherberg, Germany

  • Andreas Ruether,

    Affiliation Institute of Clinical Molecular Biology, University Hospital Schleswig-Holstein, Campus Kiel, Kiel, Germany

  • Ashish Kumar,

    Affiliation Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, United Kingdom

  • Michaela Schedel,

    Affiliation Center for Pediatrics, Clinic for Pediatric Pneumology, Allergology and Neonatology, Hannover Medical School, Hannover, Germany

  • Christian Vogelberg,

    Affiliation University Children's Hospital, Technical University Dresden, Dresden, Germany

  • Erika von Mutius,

    Affiliation University Children's Hospital, Ludwig Maximilian's University, Munich, Germany

  • Andrea von Berg,

    Affiliation Research Institute for the Prevention of Allergic Diseases, Children's Department, Marien-Hospital, Wesel, Germany

  • Albrecht Bufe,

    Affiliation Department of Experimental Pneumology, Ruhr-University, Bochum, Germany

  • Ernst Rietschel,

    Affiliation University Children's Hospital, University of Cologne, Cologne, Germany

  • Andrea Heinzmann,

    Affiliation University Children's Hospital, Albert Ludwigs University, Freiburg, Germany

  • Otto Laub,

    Affiliation Kinder- und Jugendarztpraxis Laub, Rosenheim, Germany

  • Burkhard Simma,

    Affiliation Children's Department, University Teaching Hospital, Landeskrankenhaus Feldkirch, Feldkirch, Austria

  • Thomas Frischer,

    Affiliation University Children's Hospital Vienna, Vienna, Austria

  • Jon Genuneit,

    Affiliation Institute of Epidemiology, Ulm University, Ulm, Germany

  • Ivo G. Gut,

    Affiliations Centre National de Génotypage, Institut Génomique, Commissariat à l'Énergie Atomique, Evry, France, Centro Nacional de Analisis Genomico, Barcelona, Spain

  • Stefan Schreiber,

    Affiliation Institute of Clinical Molecular Biology, University Hospital Schleswig-Holstein, Campus Kiel, Kiel, Germany

  • Mark Lathrop,

    Affiliation Centre National de Génotypage, Institut Génomique, Commissariat à l'Énergie Atomique, Evry, France

  • Thomas Illig,

    Affiliation Institute of Epidemiology, Helmholtz Centre Munich, Neuherberg, Germany

  •  [ ... ],
  • Michael Kabesch

    kabesch.michael@mh-hannover.de

    Affiliation Center for Pediatrics, Clinic for Pediatric Pneumology, Allergology and Neonatology, Hannover Medical School, Hannover, Germany

  • [ view all ]
  • [ view less ]

Abstract

The first genome wide association study (GWAS) for childhood asthma identified a novel major susceptibility locus on chromosome 17q21 harboring the ORMDL3 gene, but the role of previous asthma candidate genes was not specifically analyzed in this GWAS. We systematically identified 89 SNPs in 14 candidate genes previously associated with asthma in >3 independent study populations. We re-genotyped 39 SNPs in these genes not covered by GWAS performed in 703 asthmatics and 658 reference children. Genotyping data were compared to imputation data derived from Illumina HumanHap300 chip genotyping. Results were combined to analyze 566 SNPs covering all 14 candidate gene loci. Genotyped polymorphisms in ADAM33, GSTP1 and VDR showed effects with p-values <0.0035 (corrected for multiple testing). Combining genotyping and imputation, polymorphisms in DPP10, EDN1, IL12B, IL13, IL4, IL4R and TNF showed associations at a significance level between p = 0.05 and p = 0.0035. These data indicate that (a) GWAS coverage is insufficient for many asthma candidate genes, (b) imputation based on these data is reliable but incomplete, and (c) SNPs in three previously identified asthma candidate genes replicate in our GWAS population with significance after correction for multiple testing in 14 genes.

Introduction

Genome wide association studies (GWAS) have recently identified genetic susceptibility for many complex diseases by genotyping hundreds of thousands of single nucleotide polymorphisms (SNPs) across the genome in very large sets of patients and controls [1]. While GWAS discovered novel disease loci, many previous candidate genes identified by hypothesis driven approaches were not replicated by GWAS at a genome wide significance level. Thus, the question arose in many complex disorders if associations with previous candidate genes represent false positive results. Alternatively it was suggested that GWAS SNP chips insufficiently covered these genes of interest. Imputation was recently introduced to predict allelic status of SNPs not covered by direct genotyping in GWAS, extrapolating information from neighbouring SNPs for which genotyping data and linkage disequilibrium (LD) from HapMap was available. The value of these approaches is not yet determined. We used our GWAS data set on asthma [2] to study these open questions specifically for candidate genes in asthma.

In childhood asthma, the first GWAS [2] to which we contributed a discovery population of app. 1.400 German children, captured a major susceptibility signal from chromosome 17q21 and a few less suggestive hits in other regions of the genome, none of which was overlapping with the numerous asthma candidate genes previously identified by hypothesis driven approaches or positional cloning [3].

Therefore, we developed bioinformatics tools to systematically identify asthma candidate genes from the literature and genotyped these genes for previously suggested asthma related SNPs. We compared these results with the initial GWAS genotyping and imputation data derived from our original GWAS genotyping [2]. We used the combination of these approaches to cover all candidate genes in depth with tagging SNPs and performed a comprehensive analysis of previous candidate gene loci.

Results

We implemented a script to scan an assemblage of abstracts for specific gene names or collections of gene names as depicted in figure 1a. The script searches for single words or specific groups of words in a given text. Gene names or a group of several names for the same gene were searched in a collection of abstracts. For every gene name or cluster of gene names respectively a separate file was produced that consists of all abstracts a gene of the cluster is mentioned in. Gene names and symbols were obtained from the Human Genome Organization - Gene Nomenclature Committee [4]. We scanned all 28,472 approved symbols in the HUGO database as of October 22nd 2009. If available we clustered these gene names with their previous symbols and their attached aliases so that overall 98743 specific gene names were searched against a selection of abstracts published in the pubmed database (http://www.pubmedcentral.nih.gov/) as of October 22nd 2009. For the current study we searched for the terms “asthma” and “association” and “SNP” or “polymorphism” or “genetic”, retrieving 1,341 abstracts. With the script 825 cluster of gene names, hence genes, were extracted. Details on the use of the bioinformatics search tool for genetic association data from the literature are available upon request from the authors. All abstracts were then examined to contain significant associations (nominal p-value of <0.05) with asthma. For that purpose, different forms of asthma like childhood, adult or atopic asthma were not discriminated and reports were not weighted for study size. Associations with asthma severity or intermediate phenotypes of asthma were not considered and the analysis was restricted to studies in English language. All genes for which genetic associations with asthma were reported at the gene level in a discovery dataset and at least three independent populations (not necessarily published independently) were studied further. Direction of association and associated SNP in a gene could vary (loose replication) as indicated in table 1. Genes discovered by further GWAS studies [5], [6] were excluded in this study as their discovery relied on SNPs present at currently used SNP chips.

thumbnail
Figure 1.

a) Schematic depiction of literature query process. * genes not further investigated were (a) ORMDL3, originally discovered in this GWA; (b) discovered by other GWAS: CHI3L1, IL1RL1, MYB, WDR36 (c) NOS1 where CA-repeats and not SNPs were reported; (d) gene deletions for GSTM1, FLG. b) Number of populations that show an association between SNPs in the 14 selected candidate genes and asthma diagnosis. (based on literature query from October 22nd 2009).

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

thumbnail
Table 1. Comparison of risk allele, frequency of risk allele in cases and controls and effect size in our study population and previously published studies for candidate gene polymorphisms with a p- value of <0.05 in our population (for complete results on all tested SNPSs see table S2).

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

The most prevalent gene symbol in the retrieved abstracts for genetic association with asthma was IL4 (91 citations) followed by IL13 (75 citations) and TNF (73 citations). Significant associations of SNPs in IL4 with asthma were reported in 15 different populations. The following 14 genes distributed across the genome were associated with asthma diagnosis in more than three independent study populations (figure 1b): ADAM33, CCL5, CD14, DPP10, EDN1, GPR154, GSTP1, IL12B, IL13, IL4, IL4R, PTGDR, TNF, and VDR.

In these genes, asthma association had been reported for a total of 89 SNPs (table 2 and table S1); only 19 of these had been directly genotyped by the Illumina Sentrix HumanHap300 BeadChip used in our GWAS. Using linkage disequilibrium (LD) based on HapMap [7], a further 14 of these 89 SNPs were indirectly covered by tagging SNPs from chip-genotyping. Thus, only a total of approximately 37% of previously suggested SNPs in asthma candidate genes were captured directly or indirectly by GWAS genotyping. None of the SNPs in ADAM33, CCL5, CD14, or IL4 (table 2) were covered by the original GWAS genotyping.

thumbnail
Table 2. Number of previously associated asthma candidate gene SNPs and coverage by 1st generation GWAS genotyping, imputation and re-genotyping.

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

We aimed at re-genotyping all additional candidate gene SNPs and performed imputation based on GWAS genotyping. From 56 SNPs not covered be GWAS genotyping, genotyping assays could be designed for 45 of which 39 passed quality control. For 58 of the 89 candidate-gene derived SNP imputation was successful. Thus, for 11 of the 89 SNPs only imputation but no genotyping data and for 55 candidate SNPs (with previous association) imputation and genotyping data were available. In cases where both imputation and genotyping results were available for a specific SNP, genotyping results were always used for further analyses. Using GWAS genotyping data, re-genotyping and imputation in combination, data for 83 out of 89 candidate gene SNPs from the literature could be analyzed successfully.

SNPs in the genes ADAM33 (rs528557, p = 0.000019) and VDR (rs1540339, p = 0.0014), previously associated with asthma, were also associated with asthma in this dataset at a significance level of 0.0035 corrected for multiple testing. Previously associated SNPs in EDN1, IL13, IL4R and TNF replicated at a nominal significance level between 0.05≥p>0.0035. Effect size and effect directions were compared to previously published results in table 1.

For 241 SNPs in the 14 candidate genes both imputation and direct genotyping data were available irrespective of previously reported associations (table S2). Effect size and p-value of effects correlated tightly (figure 2), and R2 is given in table S2. Direction of effects changed 12 times, mostly when effects were close to 1. However, in one case (rs2280788, CCL5) imputation implied an odds ratio (OR) of 0.75 when re-genotyping showed an effect size of 1.22. For IL4, imputation suggested rs2070874 (C-33T) to be significantly associated with asthma (p = 0.005) but the association was not significant when re-genotyped by matrix-assisted laser desorption/ionization time-of flight (MALDI-TOF) (p = 0.054). An opposite outcome was observed for IL4R rs2057768, when an effect was significant by genotyping but not by imputation (p = 0.009 versus p = 0.2). This suggests that imputation is a reasonably reliable tool to estimate SNP effects of asthma candidate genes based on 1st generation GWAS genotyping data, when effect sizes are beyond 1.25 in a population of app. 1,000 cases. For detailed comparison between re-genotyped and inputated SNPs see figure 3. However, imputation based on GWAS BeadChip genotyping data could only increase coverage of the previously reported SNPs in asthma candidate genes to app. 60%; missing numerous SNPs of specific interest.

thumbnail
Figure 2. Correlation between effect sizes observed by concomitant genotyping and imputation of candidate gene SNPs for which both measures were available.

(red indicates change of effect direction between genotyping and imputation). For the comparison based on HapMapII (HM2) data 241 SNPs and for HapMapIII (HM3) 216 SNPs were available.

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

thumbnail
Figure 3. Comparison of effect sizes and p-values derived from regenotyped SNPs with data based on imputation.

Log transferred p-values are shown on the y-axis. Each individual SNP effect is depicted as a circle where the diameter of the circle reflects the respective odds ratios. The dark grey circles represent the regenotyping the light grey circles stand for the imputation based on HapMap2 and the transparent show the imputation data based on HapMap3. Significance level are marked as solid (p = 0.05) and dashed grey lines (p = 0.0035).

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

Thus, we next combined data from 1st generation GWAS genotyping by BeadChip, re-genotyping by MALDI-TOF and TaqMan® MGB, and imputation based on this combined set of genotypes to comprehensively saturate all 14 candidate gene loci under study with 566 tagging SNPs as well as previously associated SNPs (figure 4 and table S2). This way, further associations below the 0.0035 threshold were detected for SNPs in GSTP1 and VDR. Further SNPs in DPP10, IL12B and IL4 were now found to be associated with asthma in our population, but only at the nominal significance level of 0.05 (table S2).

thumbnail
Figure 4. Association of 566 tagging SNPs in 14 candidate genes with asthma.

Genes are positioned along the x-axis in alphabetical order. Each dot represents SNPs in their relative position within the gene. Log transferred p-values are shown on the y-axis. Each individual SNP effect is depicted as a circle where the diameter of the circle reflects the respective odds ratios. Significance level are marked as solid (p = 0.05) and dashed green lines (p = 0.0035).

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

Discussion

We systematically investigated the role of SNPs from previously suggested asthma candidate genes with more than three positive association reports in the literature. In a large German population for which 1st generation GWAS was also available, imputation was provided based on Illumina Sentrix HumanHap300 BeadChip data and compared to re-genotyping of candidate gene SNPs by MALDI-TOF and TaqMan® techniques. Thus, an in depth analysis of these previously suggested candidate genes was performed. The study was powered to detect asthma odds ratios of 1.5 at a significance level of alpha = 0.05 with a power of 95% for SNPs with a minor allele frequency of >0.1 and for odds ratios of 1.3 with a power of 80%. SNPs in three of 14 investigated candidate genes showed associations at a significance level corrected for multiple testing (p<0.0035), and further 7 genes were associated at a significance level of p<0.05. Effect directions were consistent for associated SNPs in our population compared to previous publications in populations of similar ancestry (table 1).

We expected to observe less than five positive associations at a nominal significance level of p<0.05 under the null hypothesis of no association (83×0.05 = 4.15). We observed 9 SNPs associated with asthma at this level, twice more than expected by chance. We corrected for multiple testing by multiplying the nominal significance level of alpha = 0.05 by the number of genes studied (N = 14). This seems appropriate for the purpose of replication testing. When extending the analyses to all 566 SNPs in the respective gene locus, this level of correction may be too permissive. Also, one has to keep in mind that novel associations found with the broader approach of saturating tagging SNPs lack the level of replication which the SNPs drawn from the literature have already achieved. Therefore, the associations found with rs614080 and rs6591256 in GSTP1, rs17764121, rs1396936, rs966859, rs17823623, rs12997234 and rs13034936 in DPP10, rs3213102 in IL12B and rs2243274 in IL4 need to be treated with caution until a sufficient level of replication has been achieved for these SNPs in other populations.

Imputation was compared to re-genotyping in our dataset to assess the usefulness of imputation to answer such a specific question as to what extent certain candidate gene SNPs play a role in asthma in a GWAS data set. Clearly, accuracy of imputation depends on coverage of the specific locus under investigation. With 1st generation GWAS SNP chips, gray and black areas with low or no SNP coverage exist in the genome (table S3). Interestingly, this is also the case for well known genes associated with many diseases: For CD14 and IL4, no SNP was genotyped by the Illumina Sentrix HumanHap300 BeadChip. This has improved for later versions of GWAS genotyping chips (table S3). But still, the coverage of putatively functional SNPs previously associated with disease is insufficient for many genes even in next generation GWAS chips and even when tagging SNP approaches based on LD were taken into account (table S3).

LD for our imputation was derived from HapMap2 version 2006 as well as HapMap3 version 2009. With further improvements of HapMap and the finalisation of the 1000 Genomes project, where 1000 human genomes are currently being sequenced (www.1000genomes.org), imputation information also on existing SNP sets will improve and re-analysis of existing GWAS data may be useful. Despite it's known limitations, imputation proved valuable and helpful in this project. Reliability of results was satisfactory when compared to real genotyping results in those cases where information from both sources were acquired. However, re-genotyping of major susceptibility signals found by imputation may still be useful as this represents the gold standard for ascertainment of allelic status of SNPs.

Effects of candidate gene SNPs on asthma risk are modest, not exceeding a change in risk by more than app. 50% (odds ratios between 0.58 and 1.50) which is within the expected range for genetic effects in complex diseases and not significantly different from susceptibility signals reported by GWAS in atopic diseases [8][10]: The chromosome 17q21 susceptibility locus shows an asthma effect size of OR 1.11–1.52 [2], and a further GWAS on asthma report ORs of 0.88 respectively [9]. Odds ratios for asthma identified through intermediate phenotypes by GWAS were 0.54 [5] and 1.07–1.16 [6]. A recently reported susceptibility locus for atopic dermatitis on 11q13.5 shows an overall effect of OR 1.22 [8]. These data indicate that a number of previously suggested asthma candidate gene SNPs influence asthma susceptibility at a modest level, comparable to those genes identified through GWAS in allergy genetics. Effects can be detected in large study populations but effect size may be minor and significance of effects and effect direction may vary in smaller studies or be influenced by environmental factors in specific populations (table 1). Highly significant p-values of associations in initial GWAS reports may imply but not guarantee replication of associations in further cohorts. At the same time, biological candidates such as those identified through hypothesis driven candidate gene approaches may not yield p-values significant after correction for multiple testing in GWAS but may replicate consistently over many populations tested. In some cases, as shown in table 1, replication of gene effects on asthma is on a gene level rather than on the individual SNP level, indicating either that the putatively causal SNP may not be typed or that gene by environment effects may modify effects; definitions of asthma phenotypes may vary between studies and different genetic models may have been applied in different studies [11]. Publication bias may be present for candidate genes and could have influenced the selection of genes exceeding the threshold of 3 positive replications. For all mentioned genes and susceptibility signals, systematic replications are thus required before final conclusions can be drawn. Taken all considerations into account, the strongest evidence for true association between a previous candidate SNP and asthma in this population exists for rs528557 in ADAM33, while the signal from chromosome 17q21 is the most replicated finding from GWAS studies so far.

Our analyses suggest that GWAS is a powerful tool to discover novel susceptibility loci for complex diseases and it has recently been applied very successfully. However, GWAS have limitations that may be overcome with improvements in SNP coverage and imputation. At this point GWAS are not designed to finally exclude potential candidate genes from a role in a specific disease. This is particularly the case for complex diseases where a single study can not cover all potential phenotypes, all possible gene-environment interactions and all different genetic models. Replication in as many populations as possible is still important, as well as the determination of biological function for candidate genes and SNPs in a specific disease. This is the challenge for the post-GWAS era.

Methods

Population description

An overview of the population selection for this study is given in figure 5 and described briefly here. From 2001 to 2007, the Multicenter Asthma Genetic in Childhood (MAGIC) study recruited 1255 children in seven German and Austrian clinical asthma centers located in Bochum, Cologne, Feldkirch, Freiburg, Munich/Rosenheim, Vienna and Wesel. A pediatric pulmonologist diagnosed asthma and allergies according to clinical guidelines and objective measures such as lung function tests, clinical examination, history, and allergy testing in 865 children (mean age 11 years, range 3–17 years, 64% male). Children without asthma but other conditions (such as inflammatory bowel disease) were also recruited in another leg of the study in the same clinical centers (n = 390) but have not been included in any analysis so far. All 655 children from the MAGIC study that were included in this and previous analyses [2] were diagnosed with asthma, were of German origin and had sufficient high quality DNA available at the time the genome wide analysis was performed in 2006.

thumbnail
Figure 5. Study Design and population selection strategy.

(ISAAC  =  International Study for Asthma and Allergies in Childhood; MAGICS  =  Multicentre Asthma Genetics in Childhood Study; GWAS  =  Genome wide Association Study; QC  =  quality control).

https://doi.org/10.1371/journal.pone.0013894.g005

As a reference, a random sample consisting of 800 German children was drawn from the cross sectional International Study of Asthma and Allergies in Childhood phase II (ISAAC II). 767 passed DNA quality control for GWAS and of these, 73 were asthmatics and 694 were non-asthmatics. The ISAAC II population had been described in great detail before. In brief, the study assessed the prevalence of asthma and allergies in schoolchildren aged 9–11 years in Munich and Dresden in 1995 and 1996 [12]. Asthma status was based on the parent's report of a physician's diagnosis of asthma at least once, or of spastic or asthmatic bronchitis more than once in self-administered questionnaires.

Written informed consent was obtained from all parents of children included in these studies and study methods were approved by the respective ethics committees. All children were of German origin to control for population admixture [2].

Genotyping

Technical details on GWAS genotyping of our study population using the Illumina Sentrix HumanHap300 BeadChip have been published elsewhere [2]. For additional in depth re-genotyping of candidate gene SNPs, matrix-assisted laser desorption/ionization time-of flight (MALDI-TOF) mass spectrometry (Sequenom Inc., San Diego, USA) was used as previously described [13]. SNPs failing MALDI-TOF genotyping were genotyped by a TaqMan® MGB biallelic discrimination system as previously described [14]. Deviations from Hardy–Weinberg Equilibrium (HWE) in the control group were assessed for quality control of genotyping procedures. Detailed information including specific amplification and extension primers for all genotyping procedures are provided by the authors upon request.

Statistical analyses and quality control

Statistical analyses were performed with the PLINK [15] software package (version 1.04, http://pngu.mgh.harvard.edu/purcell/plink/) and confirmed with the SAS statistical software package (version 9.2; SAS, Inc., Cary, NC). All associations of asthma with SNPs, their accompanying odds ratio and confidence intervals were evaluated by the additive logistic regression model. We restricted the analysis to one genetic model (which is the common model for GWAS) to reduce the problem of multiple testing. Power calculations were performed with the PGA software package [16].

We imputed genotypes for all polymorphic HapMap SNPs (www.hapmap.org) by using a hidden Markov model programmed in MACH [17]. The method combines genotypes from the study samples with the HapMap CEU sample (July 2006 phased haplotype release, called HapMapII and Feb 2009 phased haplotype release, called HapMapIII) and identifies the stretches of haplotype shared between the study samples and the HapMap sample. For each individual, the genotype at the untyped SNP can be summarized by taking (1) the most likely genotype according to the posterior probability of the three possible genotypes and (2) allele dosage, the expected number of copies of the reference allele (a fractional value between 0 and 2). We used the imputed most likely genotypes for association analyses.

Exclusion criteria for further analyses were a low genotyping rate (<95%) with the Illumina Sentrix HumanHap300 BeadChip, incomplete phenotype information, unusual high heterozygosity, missing age information or age above 18 years, genetic identification of cryptic cousin pairs or non-Caucasian origin. Thus, for further analyses 703 asthmatics, of which 414 were males and 658 controls (372 males), remained. Table 3 provides additional clinical information. Background stratification was measured by a genomic control parameter of 1.01. All markers which failed the HWE test in controls (p< = 0.0001), had low minor-allele frequencies (<5%) or low genotyping call rates (<95%) were also excluded. The remaining 300,745 SNPs were used for further analyses. For re-genotyping of candidate gene SNPs in children with sufficient DNA available, genotype call rates of at least 95% were achieved by MALDI-TOF and TaqMan® MGB assays. The software Haploview [18] was used to calculate linkage disequilibrium (LD) patterns and to select tagging SNPs (MAF 0.001, r2>0.8). The chi-square test was applied to analyse deviations from Hardy-Weinberg equilibrium (HWE), with expected frequencies derived from allele frequencies. The alpha level was set to 0.0035 to correct for multiple testing of 14 genes (0.05/14 = 0.0035).

thumbnail
Table 3. Characteristics of the population used for the actual study consisting of two parts, the Multicentre Asthma Genetics in Childhood Study (MAGICS) and the International Study for Asthma and Allergies in Childhood (ISAAC).

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

Supporting Information

Table S1.

Gene symbol, references, total number of SNPs and individual rs numbers of asthma associated SNPs

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

(0.28 MB DOC)

Table S2.

Candidate gene associations with asthma and comparison between genotyping and imputation results. Odds ratios, 95% confidence intervals are given for SNP effects. R2 is given for SNPs were both imputation and genotyping data is available. Bold font indicates previously reported SNPs, standard font is tagging SNPs, cursive font indicates SNPs that are neither tagging nor previously reported SNPs; * p<0.05; **p<0.0035

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

(0.10 MB PDF)

Table S3.

Comparison of candidate gene SNP coverage by GWAS genotyping systems. Number of genotyped SNPs per candidate gene locus from 5 kb upstream of 1st exon to 5 kb after last exon are given and in brackets () genotyped or covered SNPs with previously reported asthma association.

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

(0.04 MB DOC)

Acknowledgments

We thank Sonja Zeilinger and Ilona Dahmen for excellent technical assistance and Christian Schröter for his valuable contributions to data acquisition.

Author Contributions

Conceived and designed the experiments: MK. Performed the experiments: SM NK AR TI. Analyzed the data: SM MD MK. Contributed reagents/materials/analysis tools: LL MD AK MS CV EvM AvB AB ER AH OL BS TF JG IGG SS ML TI MK. Wrote the paper: SM MK.

References

  1. 1. The Wellcome Trust Case Control Consortium (2007) Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature 447: 661–678.
  2. 2. 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.
  3. 3. Vercelli D (2008) Discovering susceptibility genes for asthma and allergy. Nat Rev Immunol 8: 169–182.
  4. 4. Eyre TA, Ducluzeau F, Sneddon TP, Povey S, Bruford EA, et al. (2006) The HUGO Gene Nomenclature Database, 2006 updates. Nucleic Acids Res 34: D319–321.
  5. 5. Ober C, Tan Z, Sun Y, Possick JD, Pan L, et al. (2008) Effect of variation in CHI3L1 on serum YKL-40 level, risk of asthma, and lung function. N Engl J Med 358: 1682–1691.
  6. 6. 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.
  7. 7. (2003) The International HapMap Project. Nature 426: 789–796.
  8. 8. Esparza-Gordillo J, Weidinger S, Folster-Holst R, Bauerfeind A, Ruschendorf F, et al. (2009) A common variant on chromosome 11q13 is associated with atopic dermatitis. Nat Genet 41: 596–601.
  9. 9. 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.
  10. 10. Weidinger S, Gieger C, Rodriguez E, Baurecht H, Mempel M, et al. (2008) Genome-wide scan on total serum IgE levels identifies FCER1A as novel susceptibility locus. PLoS Genet 4: e1000166.
  11. 11. Kabesch M (2009) The art of replication. Thorax 64: 370–371.
  12. 12. Weiland SK, von Mutius E, Hirsch T, Duhme H, Fritzsch C, Werner B, et al. (1999) Prevalence of respiratory and atopic disorders among children in the East and West of Germany five years after unification. Eur Respir J 14: 862–870.
  13. 13. Kormann MS, Carr D, Klopp N, Illig T, Leupold W, et al. (2005) G-Protein-coupled receptor polymorphisms are associated with asthma in a large German population. Am J Respir Crit Care Med 171: 1358–1362.
  14. 14. Haas SL, Ruether A, Singer MV, Schreiber S, Bocker U (2007) Functional P2X7 receptor polymorphisms (His155Tyr, Arg307Gln, Glu496Ala) in patients with Crohn's disease. Scand J Immunol 65: 166–170.
  15. 15. 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. Am J Hum Genet 81: 559–575.
  16. 16. Menashe I, Rosenberg PS, Chen BE (2008) PGA: power calculator for case-control genetic association analyses. BMC Genet 9: 36.
  17. 17. Scott LJ, Mohlke KL, Bonnycastle LL, Willer CJ, Li Y, et al. (2007) A genome-wide association study of type 2 diabetes in Finns detects multiple susceptibility variants. Science 316: 1341–1345.
  18. 18. Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21: 263–265.
  19. 19. Thongngarm T, Jameekornrak A, Limwongse C, Sangasapaviliya A, Jirapongsananuruk O, et al. (2008) Association between ADAM33 polymorphisms and asthma in a Thai population. Asian Pac J Allergy Immunol 26: 205–211.
  20. 20. Van Eerdewegh P, Little RD, Dupuis J, Del Mastro RG, Falls K, et al. (2002) Association of the ADAM33 gene with asthma and bronchial hyperresponsiveness. Nature 418: 426–430.
  21. 21. Howard TD, Postma DS, Jongepier H, Moore WC, Koppelman GH, et al. (2003) Association of a disintegrin and metalloprotease 33 (ADAM33) gene with asthma in ethnically diverse populations. J Allergy Clin Immunol 112: 717–722.
  22. 22. Hirota T, Hasegawa K, Obara K, Matsuda A, Akahoshi M, et al. (2006) Association between ADAM33 polymorphisms and adult asthma in the Japanese population. Clin Exp Allergy 36: 884–891.
  23. 23. Blakey JD, Sayers I, Ring SM, Strachan DP, Hall IP (2009) Positionally cloned asthma susceptibility gene polymorphisms and disease risk in the British 1958 Birth Cohort. Thorax 64: 381–387.
  24. 24. Zhu G, Carlsen K, Carlsen KH, Lenney W, Silverman M, et al. (2008) Polymorphisms in the endothelin-1 (EDN1) are associated with asthma in two populations. Genes Immun 9: 23–29.
  25. 25. Black S, Teixeira AS, Loh AX, Vinall L, Holloway JW, et al. (2009) Contribution of functional variation in the IL13 gene to allergy, hay fever and asthma in the NSHD longitudinal 1946 birth cohort. Allergy 64: 1172–1178.
  26. 26. Heinzmann A, Mao XQ, Akaiwa M, Kreomer RT, Gao PS, et al. (2000) Genetic variants of IL-13 signalling and human asthma and atopy. Hum Mol Genet 9: 549–559.
  27. 27. Hosseini-Farahabadi S, Tavakkol-Afshari J, Rafatpanah H, Farid Hosseini R, Khaje Daluei M (2007) Association between the polymorphisms of IL-4 gene promoter (−590C>T), IL-13 coding region (R130Q) and IL-16 gene promoter (−295T>C) and allergic asthma. Iran J Allergy Asthma Immunol 6: 9–14.
  28. 28. Hytonen AM, Lowhagen O, Arvidsson M, Balder B, Bjork AL, et al. (2004) Haplotypes of the interleukin-4 receptor alpha chain gene associate with susceptibility to and severity of atopic asthma. Clin Exp Allergy 34: 1570–1575.
  29. 29. Zhang H, Zhang Q, Wang L, Chen H, Li Y, et al. (2007) Association of IL4R gene polymorphisms with asthma in Chinese populations. Hum Mutat 28: 1046.
  30. 30. Beghe B, Barton S, Rorke S, Peng Q, Sayers I, et al. (2003) Polymorphisms in the interleukin-4 and interleukin-4 receptor alpha chain genes confer susceptibility to asthma and atopy in a Caucasian population. Clin Exp Allergy 33: 1111–1117.
  31. 31. Zhang W, Zhang X, Qiu D, Sandford A, Tan WC (2007) IL-4 receptor genetic polymorphisms and asthma in Asian populations. Respir Med 101: 186–190.
  32. 32. Ober C, Leavitt SA, Tsalenko A, Howard TD, Hoki DM, et al. (2000) Variation in the interleukin 4-receptor alpha gene confers susceptibility to asthma and atopy in ethnically diverse populations. Am J Hum Genet 66: 517–526.
  33. 33. Li YF, Gauderman WJ, Avol E, Dubeau L, Gilliland FD (2006) Associations of tumor necrosis factor G-308A with childhood asthma and wheezing. Am J Respir Crit Care Med 173: 970–976.
  34. 34. Castro-Giner F, Kogevinas M, Machler M, de Cid R, Van Steen K, et al. (2008) TNFA -308G>A in two international population-based cohorts and risk of asthma. Eur Respir J 32: 350–361.
  35. 35. Wu H, Romieu I, Sienra-Monge JJ, del Rio-Navarro BE, Anderson DM, et al. (2007) Parental smoking modifies the relation between genetic variation in tumor necrosis factor-alpha (TNF) and childhood asthma. Environ Health Perspect 115: 616–622.
  36. 36. Munthe-Kaas MC, Carlsen KL, Carlsen KH, Egeland T, Haland G, et al. (2007) HLA Dr-Dq haplotypes and the TNFA-308 polymorphism: associations with asthma and allergy. Allergy 62: 991–998.
  37. 37. Kim HB, Kang MJ, Lee SY, Jin HS, Kim JH, et al. (2008) Combined effect of tumour necrosis factor-alpha and interleukin-13 polymorphisms on bronchial hyperresponsiveness in Korean children with asthma. Clin Exp Allergy 38: 774–780.
  38. 38. Kumar V, Khosla R, Gupta V, Sarin BC, Sehajpal PK (2008) Differential association of tumour necrosis factor-alpha single nucleotide polymorphism (-308) with tuberculosis and bronchial asthma. Natl Med J India 21: 120–122.
  39. 39. Witte JS, Palmer LJ, O'Connor RD, Hopkins PJ, Hall JM (2002) Relation between tumour necrosis factor polymorphism TNFalpha-308 and risk of asthma. Eur J Hum Genet 10: 82–85.
  40. 40. Kumar A, Gupta V, Changotra H, Sarin BC, Sehajpal PK (2008) Tumor necrosis factor—alpha and transforming growth factor—beta1 polymorphisms in bronchial asthma. Indian J Med Sci 62: 323–330.
  41. 41. Wang TN, Chen WY, Wang TH, Chen CJ, Huang LY, et al. (2004) Gene-gene synergistic effect on atopic asthma: tumour necrosis factor-alpha-308 and lymphotoxin-alpha-NcoI in Taiwan's children. Clin Exp Allergy 34: 184–188.
  42. 42. Jimenez-Morales S, Velazquez-Cruz R, Ramirez-Bello J, Bonilla-Gonzalez E, Romero-Hidalgo S, et al. (2009) Tumor necrosis factor-alpha is a common genetic risk factor for asthma, juvenile rheumatoid arthritis, and systemic lupus erythematosus in a Mexican pediatric population. Hum Immunol 70: 251–256.
  43. 43. Winchester EC, Millwood IY, Rand L, Penny MA, Kessling AM (2000) Association of the TNF-alpha-308 (G—>A) polymorphism with self-reported history of childhood asthma. Hum Genet 107: 591–596.
  44. 44. Chagani T, Pare PD, Zhu S, Weir TD, Bai TR, et al. (1999) Prevalence of tumor necrosis factor-alpha and angiotensin converting enzyme polymorphisms in mild/moderate and fatal/near-fatal asthma. Am J Respir Crit Care Med 160: 278–282.
  45. 45. Moffatt MF, Cookson WO (1997) Tumour necrosis factor haplotypes and asthma. Hum Mol Genet 6: 551–554.
  46. 46. Shin HD, Park BL, Kim LH, Jung JH, Wang HJ, et al. (2004) Association of tumor necrosis factor polymorphisms with asthma and serum total IgE. Hum Mol Genet 13: 397–403.
  47. 47. Albuquerque RV, Hayden CM, Palmer LJ, Laing IA, Rye PJ, et al. (1998) Association of polymorphisms within the tumour necrosis factor (TNF) genes and childhood asthma. Clin Exp Allergy 28: 578–584.
  48. 48. Poon AH, Laprise C, Lemire M, Montpetit A, Sinnett D, et al. (2004) Association of vitamin D receptor genetic variants with susceptibility to asthma and atopy. Am J Respir Crit Care Med 170: 967–973.
  49. 49. Raby BA, Lazarus R, Silverman EK, Lake S, Lange C, et al. (2004) Association of vitamin D receptor gene polymorphisms with childhood and adult asthma. Am J Respir Crit Care Med 170: 1057–1065.