Open Access Highly Accessed Open Badges Research article

Defining the contribution of SNPs identified in asthma GWAS to clinical variables in asthmatic children

Asif S Tulah13, John W Holloway2 and Ian Sayers1*

Author Affiliations

1 Division of Respiratory Medicine, Queen’s Medical Centre, University of Nottingham, Nottingham NG7 2UH, United Kingdom

2 Human Genetics and Medical Genomics, Human Development and Health, Faculty of Medicine, University of Southampton, Southampton, UK

3 Institute of Cellular Medicine, Faculty of Medical Sciences, Newcastle University, Newcastle upon Tyne, UK

For all author emails, please log on.

BMC Medical Genetics 2013, 14:100  doi:10.1186/1471-2350-14-100

The electronic version of this article is the complete one and can be found online at:

Received:14 June 2013
Accepted:24 September 2013
Published:25 September 2013

© 2013 Tulah et al.; licensee BioMed Central Ltd.

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 work is properly cited.



Asthma genome-wide association studies (GWAS) have identified several asthma susceptibility genes with confidence; however the relative contribution of these genetic variants or single nucleotide polymorphisms (SNPs) to clinical endpoints (as opposed to disease diagnosis) remains largely unknown. Thus the aim of this study was to firstly bridge this gap in knowledge and secondly investigate whether these SNPs or those that are in linkage disequilibrium are likely to be functional candidates with respect to regulation of gene expression, using reported data from the ENCODE project.


Eleven of the key SNPs identified in eight loci from recent asthma GWAS were evaluated for association with asthma and clinical outcomes, including percent predicted FEV1, bronchial hyperresponsiveness (BHR) to methacholine, severity defined by British Thoracic Society steps and positive response to skin prick test, using the family based association test additive model in a well characterised UK cohort consisting of 370 families with at least two asthmatic children.


GSDMB SNP rs2305480 (Ser311Pro) was associated with asthma diagnosis (p = 8.9×10-4), BHR (p = 8.2×10-4) and severity (p = 1.5×10-4) with supporting evidence from a second GSDMB SNP rs11078927 (intronic). SNPs evaluated in IL33, IL18R1, IL1RL1, SMAD3, IL2RB, PDE4D, CRB1 and RAD50 did not show association with any phenotype tested when corrected for multiple testing. Analysis using ENCODE data provides further insight into the functional relevance of these SNPs.


Our results provide further support for the role of GSDMB SNPs in determining multiple asthma related phenotypes in childhood asthma including associations with lung function and disease severity.

Asthma; ENCODE; eQTL; GWAS; Clinical endpoints


Asthma is a complex respiratory disease with susceptibility involving the interplay between genetic and environmental factors [1]. The chronic inflammation of the airways is characterised by reversible airflow obstruction and increased hyperresponsiveness of the airways. Genome-wide association studies (GWAS) are currently the preferred method for the identification of genetic determinants in complex multifactorial diseases such as asthma. These studies have the advantage of being able to identify polymorphisms with small effect sizes and localise smaller susceptibility regions due to linkage disequilibrium (LD). The first GWAS for asthma was completed involving 994 childhood-onset asthma patients and 1,243 non-asthma controls and identified single nucleotide polymorphisms (SNPs) within chromosome 17q21 spanning the ORMDL3/GSDMB genes [2]. These SNPs showed highly significant association with childhood asthma and were also replicated in two independent cohorts [2]. A study in an Icelandic cohort showed significant association with asthma and a SNP in IL1RL1[3] and SNPs in PDE4D have also shown association with genome-wide significance with asthma [4]. IL-13 and RAD50 on 5q31.1 and the HLA-DR/DQ locus at 6p21.3 have been implicated by Li et al. [5]. Other polymorphisms spanning IL33, IL18R1, SMAD3, IL2RB and CRB1 have been proposed by other GWAS using large scale discovery and replication approaches [6,7]. However, due to the phenotype used for case definition, which is often a physician diagnosis or self-report of asthma, there remains a need to determine the contribution of these polymorphisms to clinically relevant endpoints in asthma. This will provide a greater insight into the mechanisms altered in individuals carrying the relevant risk alleles.

One problem that arises from GWAS is that LD patterns can result in multiple polymorphisms at a locus being associated with disease phenotype, even if only one of them is causal. A problem therefore exists in the identification of potential causal variants. This can be simply analysed for polymorphisms within coding regions, as they could potentially affect the polypeptide structure and function and those SNPs in regulatory or intronic regions could affect splice site motifs, transcriptional regulation or RNA stability [8,9]. However, many GWAS have implicated polymorphisms in non transcribed regions far away from annotated genes, where their role is likely to be regulatory [10]. The ENCODE project was established to identify all the functional elements such as protein binding sites, transcripts (coding and non coding) and chromatin marks in the human genome. This represents a valuable resource to provide additional meaning to association data provided from GWAS [11,12]. Taken in consideration with LD data from the 1000 Genomes project [13], computational analysis will enable assignment of more meaningful potential function to SNPs identified from GWAS or those SNPs which are in LD with a GWAS associated SNP [10], to help better prioritise SNPs for functional analysis.

The aim of the present study was to complete a SNP for SNP analysis of the key variants identified in the Caucasian population meeting genome-wide significance and/or replication identified 2007–2011 on clinical variables underlying asthma (Table 1). As these SNPs have previously shown genome-wide association, we also sought to investigate the potential functional significance of these eleven SNPs and those in LD (with these SNPs) by mining the ENCODE dataset [14] and also investigating whether they show association with expression quantitative trait loci (eQTL).

Table 1. Key SNPs identified in asthma GWAS included in the current analyses


Asthma cohort characteristics

The asthma family cohort was made up of United Kingdom families from the Southampton and Nottingham areas (n = 370 families were used in this study) (Table 2). From the Southampton area, 341 Caucasian families (n = 1508) with at least two biological siblings with physician diagnosed asthma were recruited. This cohort has been described in detail previously [15]. Baseline FEV1 (forced expiratory volume in one second) was measured as best of three values within 5% performed using Vitalograph dry-wedge bellows spirometer (Vitalograph Ltd, Buckingham, UK). This was determined 14 days after respiratory tract infection or use of bronchodilator or anti-allergic medication. 46 Caucasian families (n = 184) with at least two biological siblings with physician diagnosed asthma were recruited from the Nottingham area [15]. FEV1 was defined as the best of three values. All subjects were classified according to British Thoracic Society Step Guidelines (BTS steps, ranging from step 1 to step 5) based on physician prescribed medication [16]. Ethical approval for the Southampton subjects was obtained from the Southampton and South West Hampshire and the Portsmouth and South East Hampshire Local Research Ethics Committees and for the Nottingham subjects, the Nottingham University Medical School Ethics Committee. Written informed consent for study participation was obtained from participants (from the parent/guardian for the children).

Table 2. Clinical characteristics of UK families

SNP selection and genotyping

Eleven SNPs were chosen that have been identified from asthma GWAS identified 2007–2011 (Table 1). These were in eight loci; ORMDL3/GSDMB, IL33, IL18R/IL1RL1, SMAD3, IL2RB, PDE4D, CRB1 and RAD50. SNPs were genotyped using KASPar technology by KBiosciences (Hertfordshire, UK). For quality control of genotyping data, Chi square was used to test for any deviation between the observed genotype frequencies and the expected, under the Hardy Weinberg Equilibrium.

Association analyses

Association analysis in the family cohorts between the GWAS SNPs and asthma-related traits was conducted using the family based association test (FBAT) software (version 1.5.1) [17] in the additive model. Clinically relevant endpoints including; asthma diagnosis, Forced Expiratory Volume in one second (FEV1% Predicted), bronchial hyperresponsiveness to methacholine (BHR), British Thoracic Society (BTS) defined severity [18] and atopy (positive response to skin prick test to any one of a panel of allergens as outlined [15]). We completed power calculations based on the method of Risch & Merikangas [19]. There was between 0.93 and 0.98 power to detect an association with asthma with a significance level of p = 0.05 and a relative risk of 1.5. There was between 0.51 and 0.69 power to detect an association with p = 0.05 and a relative risk of 1.25 (calculations based on minor allele frequency (min 0.17 to max 0.48 observed). Power curves for these MAF are provided in Additional file 1. To address multiple testing, we used Bonferonni correction, which considered p < 9.1×10-4 significant. This was based on our eleven genotyped SNPs analysed for five phenotypes (0.05/55 total tests = 9.1×10-4 as a threshold p-value).

Additional file 1. Power of the study based on different relative risks for minor allele frequencies of 0.17 and 0.48. Cohort size n = 370 families and p-value of 0.05. Calculation based on the method of Risch & Merikangas.

Format: DOC Size: 175KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Bioinformatic analysis of genotyped SNPs

HaploReg [14] (a web based tool) was used to provide information about the sequence surrounding a SNP. It has been developed to mine chromatin state information from nine human cell lines; GM12878, H1, HMEC, HSMM, HepG2, Huvec, K562, NHEK and NHLF. The tool provides information on SNPs that are in LD with the queried SNP (r2 was set to 1 for these analyses and the population chosen for analyses was of European descent) and also gives information on which of these SNPs lie in enhancer histone marks, DNase sites and protein binding regions. It can also inform on the conservation of the relevant sequence and any changes to regulatory motifs based on the SNP allele changes. Each of the eleven genotyped SNPs were entered into HaploReg and analysis was completed with r2 = 1. A more detailed analysis was completed for any SNP which showed significant association in the population study.

We used the program mRNA by SNP browser (version 1.0.1) to investigate whether any of the GWAS SNPs genotyped as part of this study were associated with regions of known asthma expression quantitative trait loci (eQTL). This program incorporates a database of eQTL from asthma studies and provides a graphical output to browse data from GWAS. It provides association data between 54,675 transcripts and 406,912 SNPs. Statistics are provided for those SNPs with association p < 0.001. We used the program to search for each SNP in turn by inputting either the gene name or the SNP rs number, looking to see whether the association reached genome-wide significance (p < 10-8) [20]. We also used the eQTL browser from the University of Chicago ( webcite) which considers eQTLs from different tissue sources. A second stage of this analysis was to investigate whether these same SNPs showed genome-wide association with eQTL from a genome-wide search in 1,111 human lung samples [21].


Clinical characteristics and genotyping

Clinical characteristics for the asthma families are shown in Table 2. Mean age was 13.3 years and 10.3 years old, mean FEV1 (% Predicted) was 95.4 ± 15.6 and 96.2 ± 14.9 and % atopy (SPT) was 73.9% and 65.1% for sibling 1 and 2 respectively. These asthmatic children had predominantly mild-moderate disease as outlined by BTS guidelines. Genotyping was completed using KASPar (KBioscience) and quality control completed as outlined in previous studies [15]. All SNPs were in Hardy Weinberg equilibrium and allele frequencies in the UK families were similar to those previously reported (Table 1).

Association with asthma, FEV1 and BHR

A GSDMB SNP at the 17q21 locus, rs2305480 (Ser311Pro) was associated with asthma diagnosis, with the T allele conferring risk to develop asthma as previously reported [6] (p = 8.9×10-4, z = -3.23) (Table 3). Another GSDMB SNP rs11078927 (intronic) showed similar findings with asthma but did not reach statistical significance (p = 9.5×10-4, z = -3.23) although this SNP was significantly associated with FEV1 (% Predicted) (p = 8.7×10-4, z = -3.33). Significant association was also observed with rs2305480 (GSDMB, Ser311Pro) and BHR (PC20 methacholine, 4 mg/ml) (p = 8.2×10-4, z = -3.35). No other SNP association survived correction for multiple testing (Table 3).

Table 3. Association of key GWAS SNPs with asthma diagnosis, percent predicted FEV1 and bronchial hyperresponsiveness to methacholine

Association with asthma severity and atopy

Both GSDMB SNPs showed significant association with severity of disease defined by BTS steps (1–5), both with the same direction of association; rs2305480 (p = 1.5×10-4, z = 3.79) and rs11078927 (p = 2.2×10-4, z = -3.70). No other SNP association survived correction for multiple testing (Table 4).

Table 4. Association of key GWAS SNPs with asthma severity defined by BTS steps and atopy

Bioinformatic analysis of SNPs using data from ENCODE

We sought to investigate whether data collected from the ENCODE Consortium [12] could enrich the potential functional significance of our genotyped GWAS associated SNPs and also those SNPs which were in LD (Table 5). The two GSDMB SNPs were found to be in the same LD block (r2 = 1), with 24 other SNPs. rs2244012 (RAD50/intronic) was found to be in an LD block with 50 other SNPs. The other genotyped SNPs were in smaller LD blocks. Most of our genotyped SNPs analysed were found to have at least one SNP (in their LD block) existing in an enhancer histone mark site, DNase site or a region where protein binding consensus sequences exist, highlighting a potential functional effect of these particular SNPs. There was evidence for LD SNPs altering regulatory motifs, with 14 SNPs in the GSDMB LD block and 25 SNPs in the RAD50 LD block the highest numbers (Table 5).

Table 5. Analysis of key GWAS SNPs with data from ENCODE

Analyses of asthma associated GSDMB SNPs using ENCODE data

As the two GSDMB SNPs (rs2305480/Ser311Pro and rs11078927/intronic) showed significant association in our study, we sought to further investigate potential functional implications using the ENCODE dataset. These two GSDMB SNPs themselves were not predicted to be important from a regulatory role, however they were in a LD block with 24 other SNPs. Almost all of the SNPs in this LD block were in non-coding/intronic regions, with the exception of rs10852935, which is a synonymous polymorphism in ZPBP2 and the genotyped rs2305480/Ser311Pro (Table 6). Twenty four of these SNPs were not conserved (according to the GERP conservation score). Three SNPs were located in enhancer histone marks (rs12709365 and rs13380815 both ZPBP2/intronic and identified in cell line GM12878), rs4795399 (GSDMB/intronic) was identified in a HepG2 cell line and also located in a DNase site identified in a human hepatocyte cell line (Huh7.5) and a human leukaemic cell line (CMK). While no SNPs were predicted to occur in a protein binding regions, 14 SNPs were predicted to alter regulatory binding motifs (Table 6).

Table 6. Analysis of the GSDMB LD block (r2 =1) using data from ENCODE

Analysis of asthma eQTL data

We searched for the eleven SNPs genotyped as part of this study in turn using the mRNA by SNP browser resource, looking to see whether the SNP was associated with the effect meeting genome-wide significance (p < 10-8). Information for three of the SNPs (rs230548/GSDMB, rs2886098/CRB1 and rs2244012/RAD50) was not covered by the database and there was no information for SNPs in LD with these. Another two SNPs were also not covered by the database, however there was information for a SNP in high LD. For rs11078927/GSDMB, rs10445308 was in LD (r2 = 1.0) and this SNP showed genome-wide significance with multiple ORMDL3 probes (all C-allele); 223259_at (p = 4.4×10-13, LOD = 11.39, H2 = 14.75) and probe 235136_at (p = 8.0×10-12, LOD = 10.16, H2 = 13.3) with average ORMDL3 (p = 3.9×10-22, LOD = 20.32, H2 = 26.36). For rs3771166/IL18R1, rs54988956 was in LD (r2 = 0.96), but this was not associated with genome-wide significance with any probe. For the other six SNPs, there was no genome-wide significance. Using the University of Chicago eQTL resource, data was available for the two GSDMB SNPs (rs2305480 and rs11078927) which provided some evidence that these were eQTLs and exon-QTLs in lymphoblastoid cell lines, monocytes, liver, brain and T-cells. More specifically, both rs2305480 and rs11078927 were eQTL for IKZF3, ORMDL3 and an exon-QTL for KRT222P, NR1D1 and ORMDL3. Only three of our eleven genotyped SNPs (rs2305480/GSDMB, rs3939286/IL33 and rs744910/SMAD3) were included in the genotyping platform for the genome-wide search for eQTL in 1,111 human lung samples [21]. None of these SNPs were found in the reported data meeting 10% False Discovery Rate (FDR). As a further analysis we also searched this resource for any of the SNPs in high LD (r2 = 1) with our genotyped SNPs. rs1337167/CRB1 met the 10% FDR. This intron SNP was predicted to alter an AP-1 transcription factor binding site according to HaploReg.


GWAS have been successful in identifying susceptibility genes for asthma through the use of large scale population cohorts and replication approaches [2-7]. However, there remains a need to determine the contribution of these polymorphisms to clinically relevant endpoints in asthma as opposed to disease diagnosis. We adopted a two step process where we selected SNPs from key genes in the Caucasian population which have met genome-wide significance and/or with replication, namely eleven SNPs in the following genes; ORMDL3/GSDMB, IL33, IL18R/IL1RL1, SMAD3, IL2RB, PDE4D, CRB1 and RAD50. These SNPs were genotyped in a UK asthma family cohort with at least two affected siblings, to test for association with clinically relevant endpoints in asthma. Our data provide support for the original GSDMB association; however SNPs in the other genes did not show association after correction for multiple testing. Secondly, as these SNPs have previously shown association at a genome-wide level, we also investigated the potential for these SNPs and those in LD to affect chromatin states and alter regulatory motifs/binding sites to try to assign potential function to these associations. This analysis used recently reported data from the ENCODE project. These data support the hypothesis that SNPs associated in GWAS may be tagging SNPs for the actual causative variant, but also provide tentative evidence that these genotyped SNPs may themselves play a functional role e.g. GSDMB rs2305480 (Ser311Pro).

GWAS have the advantage over linkage studies in that they can identify polymorphisms with small effect sizes and localise smaller susceptibility loci as LD only generally spans <500 kb [22]. Since the identification of childhood asthma associated variants on chromosome 17q21 [2], many other GWAS have been completed in different populations implicating variants spanning different genes (see above). One problem with these studies lies in defining the contribution of these associated variants to clinical endpoints in asthma. To this end we have completed a SNP for SNP analysis of key variants identified in the Caucasian population meeting genome-wide significance or where there was replication, focussing on GWAS from 2007–2011 (Table 1). We utilised a well characterised UK family based asthma cohort [15] to test for association with asthma diagnosis and the following clinical endpoints: FEV1 (% predicted), BHR (to methacholine), BTS defined severity and positive SPT to one or more allergens, see [15]. These data show significant association for the two GSDMB SNPs; rs2305480 (Ser311Pro) and rs11078927 (intron) across more than one phenotype, with other phenotypes showing association which does not reach significance after correction for multiple testing. These data both replicate and extend previous findings [2,6] and identify this chromosome 17 locus as containing genetic determinants underlying multiple clinical features of asthma, including FEV1 (% Predicted) and disease severity defined by BTS. These two polymorphisms are in high LD suggesting these associations represent a single locus. The BHR association (rs2305480) is supported by different SNPs in this gene also showing association with this phenotype, albeit in different ethnic populations [23,24]. Interestingly, a recent study has found association between rs2305480 and wheezing phenotypes in asthmatic children, but did not find association with intermediate phenotypes such as atopy or lung function [25].

When interpreting the results of GWAS, the associated variant could be tagging another variant in a larger region of LD [26]. This is particularly important as it is common that only a few variants in a haplotype block are present on genotyping platforms. Methods such as imputation and using genotyping data from the 1000 Genomes project can help to further characterise these patterns [27]. However, once these variants are identified it is still necessary to assign potential function and with many associated variants from GWAS residing in non-transcribed regions, they are likely to play a regulatory role in gene expression. We used HaploReg [14] to screen chromatin data, conservation data and regulatory motif databases which was obtained from the ENCODE project to see if any biological plausible role could be assigned to the SNPs genotyped in this study and those SNPs which were in LD with r2 = 1. Our results show that these SNPs in LD can potentially lead to changes in regulatory element binding motifs affecting the transcription of genes positively or negatively and there is also the suggestion that alterations in enhancer histone marks or DNase sites could be affected. The non-synonymous SNP (Ser 311 Pro) is predicted to be deleterious to protein structure (Polyphen) in Gasdermin B suggesting a potential functional mechanism altering protein function. This may be occurring at the mRNA level (altering the efficiency of splicing) as this polymorphism is predicted to be located in an exonic splicing enchancer site (FuncPred). Both GSDMB SNPs are in high LD with several SNPs spanning ZPBPZ and GSDMB which result in alterations in enhancer and transcription factor binding sites. However the rs2305480 and rs11078927 SNPs themselves are not predicted to have any regulatory role (ENCODE data). Gasdermin-family proteins have been implicated in TGFβ1 signalling and epithelial cell apoptosis [28], with increased airway cell apoptosis reported in severe asthma [29]. The extended chromosome 17 locus has also been implicated previously as containing SNPs associated with asthma spanning three main genes; ZPBP2/GSDMB/ORMDL3 and it is still unclear which genes/polymorphisms underlie these associations. The finding that the GSDMB SNPs including a non-synonymous SNP are particularly relevant to clinical endpoints including lung function, BHR and severity in our childhood asthma cohort is interesting as few recent studies have genotyped this GWAS SNP or investigated clinical outcomes with the same criteria such as BTS used in this study [24,30]. However, a recent study has investigated the effect of rs7216389 (which is in linkage disequilibrium with both GSDMB SNPs in this study) with respect to severe asthma subjects classified with poorly controlled disease whilst taking high doses of inhaled corticosteroids, long-acting bronchodilators and short acting β2 agonists and showed association [31]. Similarly, while not meeting genome-wide association significance, we have previously provided evidence that rs2305480 is associated with severe asthma (p = 5.5×10-5, [32]), although we do acknowledge that the criteria for BTS severity is defined in part by lung function and so this association with severe asthma may be related to a synergistic effect. These data support previous studies that identified the chromosome 17 locus as containing markers of childhood onset disease [2,6,24,30], BHR [24] and elevated total serum IgE [24] in diverse ethnic populations.

Interestingly, no other SNP survived correction for multiple testing for any phenotype analysed, including asthma diagnosis and there was no indication of any other SNP showing a trend towards significance. CRB1 rs2786098 (p = 0.089 in the current study) has shown replicated association in multiple childhood cohorts, but did not reach statistical significance in our study. Perhaps this lack of association may be due to power (MAF 0.22, see earlier). One potential explanation for the lack of association for asthma and related phenotypes for the majority of SNPs is that many of the SNPs/genes were identified in adult asthma populations and the relative contribution to childhood asthma remains unclear at this time. However, in the results not meeting conventional correction for multiple testing, IL33 SNP rs1342326G (25.7 kb from gene transcription start site (TSS)) showed a level of association with atopic asthma defined by positive skin prick test and clinical diagnosis (p = 0.0097, Table 4). These data are in line with the accumulating evidence suggesting a role of IL33 in allergic mechanisms in asthma (Reviewed in [33]) although we cannot exclude that this is a false positive. Similarly, while not surviving correction, suggestive evidence for SNP associations for CRB1/BHR, IL18R/severity and IL33/SPT were apparent (p < 0.05), which may represent true associations.

Studies have shown the potential for GWAS SNPs to be tagging SNPs associated with eQTL [10]. Our analyses did provide some support for this particularly for the two GSDMB SNPs (rs2305480 and rs11078927) when eQTL analyses was completed across multiple cell types and tissues. These eQTL analyses identified regulation of multiple genes in the region, demonstrating the need for further work to define the biology underlying these genetic signals.

It is important to note the limitations of our study, while a SNP for SNP approach focussed to the pivotal SNPs from the GWAS has strengths to directly replicate and extend previous findings we acknowledge the assumption that the linkage disequilibrium pattern between the discovery Caucasian population and our UK population will be the same, which may at least in part explain the inability to detect associations with all variants. Similarly, we did not genotype all GWA significant SNPs identified in each of the regions. We acknowledge the reduced power of the current study using 370 families is a limitation potentially explaining while several SNPs showed nominal significance and did not meet stringent statistical criteria. A cohort of 602 families would be required to provide a power of 0.99 for the lowest MAF observed in the current study (relative risk of 1.5 and significance level of 0.05). We also acknowledge that the design of this study compared to discovery cohorts using case/control design also has reduced power, however we consider the use of an alternative association design i.e. family based to those previously published a strength. Furthermore the lack of assessing gene-environment interactions in GWAS and the current study may also account for not replicating previous observed effects due to differential environmental exposures between study populations. Finally, the potential for winners curse bias is also possible for direct replication of e.g. association with asthma diagnosis.


In summary, this study has provided a greater insight into the relative contribution of asthma GWAS SNPs to clinically relevant endpoints in asthmatic children and further defines the diverse role of GSDMB SNPs particularly in lung function and disease severity.


BHR: Bronchial hyperresponsiveness; BTS: British Thoracic society; CRB1: Crumbs homolog 1 precursor; DENND1B: Denn/madd domain-containing 1B; ENCODE: ENcyclopedia Of DNA elements; eQTL: Expression quantitative trait loci; FBAT: Family based association test; FEV1 (% Predicted): forced expiratory volume in one second (percent predicted); GSDMB: Gasdermin B; GWA: Genome-wide association; GWAS: Genome-wide association study; IgE: Immunoglobulin E; IL18R: Interleukin 18 receptor; IL1RL1: Interleukin 1 receptor like 1; IL2RB IL2: Receptor beta; IL33: Interleukin 33; LD: Linkage disequilibrium; MAF: Minor allele frequency; ORMDL3: Orm1-like protein 3; PDE4D: Phosphodiestase 4D, cAMP specific; Pro: Proline (amino acid); RAD50 S. Cerevisiae: homolog of (RAD50); SPT: Skin prick test; SMAD3: Mothers against decaptentaplegic drosophila homolog 3; SNP: Single nucleotide polymorphism; Ser: Serine (amino acid); TGFβ1: Transforming growth factor beta 1; TSS: Transcription start site; ZPBP2: Zona pellucid binding protein 2 gene.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

IS, JWH and AST designed the study and drafted the manuscript. AST completed the analyses. All authors contributed to the final version of the manuscript. All authors read and approved the final manuscript.


This work was funded by the University of Nottingham.


  1. Yang IA, Savarimuthu S, Kim ST, Holloway JW, Bell SC, Fong KM: Gene-environmental interaction in asthma.

    Curr Opin Allergy Clin Immunol 2007, 7(1):75-82. PubMed Abstract | Publisher Full Text OpenURL

  2. Moffatt MF, Kabesch M, Liang L, Dixon AL, Strachan D, Heath S, Depner M, von Berg A, Bufe A, Rietschel E, et al.: Genetic variants regulating ORMDL3 expression contribute to the risk of childhood asthma.

    Nature 2007, 448(7152):470-473. PubMed Abstract | Publisher Full Text OpenURL

  3. Gudbjartsson DF, Bjornsdottir US, Halapi E, Helgadottir A, Sulem P, Jonsdottir GM, Thorleifsson G, Helgadottir H, Steinthorsdottir V, Stefansson H, et al.: Sequence variants affecting eosinophil numbers associate with asthma and myocardial infarction.

    Nat Genet 2009, 41(3):342-347. PubMed Abstract | Publisher Full Text OpenURL

  4. Himes BE, Hunninghake GM, Baurley JW, Rafaels NM, Sleiman P, Strachan DP, Wilk JB, Willis-Owen SA, Klanderman B, Lasky-Su J, et al.: Genome-wide association analysis identifies PDE4D as an asthma-susceptibility gene.

    Am J Hum Genet 2009, 84(5):581-593. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Li X, Howard TD, Zheng SL, Haselkorn T, Peters SP, Meyers DA, Bleecker ER: Genome-wide association study of asthma identifies RAD50-IL13 and HLA-DR/DQ regions.

    J Allergy Clin Immunol 2010, 125(2):328-335 e311. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Moffatt MF, Gut IG, Demenais F, Strachan DP, Bouzigon E, Heath S, von Mutius E, Farrall M, Lathrop M, Cookson WO: A large-scale, consortium-based genomewide association study of asthma.

    N Engl J Med 2010, 363(13):1211-1221. PubMed Abstract | Publisher Full Text OpenURL

  7. Sleiman PM, Flory J, Imielinski M, Bradfield JP, Annaiah K, Willis-Owen SA, Wang K, Rafaels NM, Michel S, Bonnelykke K, et al.: Variants of DENND1B associated with asthma in children.

    N Engl J Med 2010, 362(1):36-44. PubMed Abstract | Publisher Full Text OpenURL

  8. Birney E, Lieb JD, Furey TS, Crawford GE, Iyer VR: Allele-specific and heritable chromatin signatures in humans.

    Hum Mol Genet 2010, 19(R2):R204-209. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Cheung VG, Spielman RS: Genetics of human gene expression: mapping DNA variants that influence gene expression.

    Nat Rev Genet 2009, 10(9):595-604. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Schaub MA, Boyle AP, Kundaje A, Batzoglou S, Snyder M: Linking disease associations with regulatory information in the human genome.

    Genome Res 2012, 22(9):1748-1759. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. EncodeProjectConsortium: The ENCODE (ENCyclopedia Of DNA Elements) Project.

    Science 2004, 306(5696):636-640. PubMed Abstract | Publisher Full Text OpenURL

  12. Dunham I, Kundaje A, Aldred SF, Collins PJ, Davis CA, Doyle F, Epstein CB, Frietze S, Harrow J, Kaul R, et al.: An integrated encyclopedia of DNA elements in the human genome.

    Nature 2012, 489(7414):57-74. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. 1000GenomesProjectConsortium: A map of human genome variation from population-scale sequencing.

    Nature 2010, 467(7319):1061-1073. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Ward LD, Kellis M: HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants.

    Nucleic Acids Res 2011, 40(Database issue):D930-934. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Barton SJ, Koppelman GH, Vonk JM, Browning CA, Nolte IM, Stewart CE, Bainbridge S, Mutch S, Rose-Zerilli MJ, Postma DS, et al.: PLAUR polymorphisms are associated with asthma, PLAUR levels, and lung function decline.

    J Allergy Clin Immunol 2009, 123(6):1391-1400 e1317. PubMed Abstract | Publisher Full Text OpenURL

  16. BTS: British guidelines on the management of asthma.

    Thorax 2008, 63(Suppl 4):iv1-121. PubMed Abstract | Publisher Full Text OpenURL

  17. Horvath S, Xu X, Laird NM: The family based association test method: strategies for studying general genotype–phenotype associations.

    Eur J Hum Genet 2001, 9(4):301-306. PubMed Abstract | Publisher Full Text OpenURL

  18. BTS: British Guideline on the Managment of Asthma. London: British Thoracic Society; 2009. OpenURL

  19. Risch N, Merikangas K: The future of genetic studies of complex human diseases.

    Science 1996, 273(5281):1516-1517. PubMed Abstract | Publisher Full Text OpenURL

  20. Dixon AL, Liang L, Moffatt MF, Chen W, Heath S, Wong KC, Taylor J, Burnett E, Gut I, Farrall M, et al.: A genome-wide association study of global gene expression.

    Nat Genet 2007, 39(10):1202-1207. PubMed Abstract | Publisher Full Text OpenURL

  21. Hao K, Bosse Y, Nickle DC, Pare PD, Postma DS, Laviolette M, Sandford A, Hackett TL, Daley D, Hogg JC, et al.: Lung eQTLs to help reveal the molecular underpinnings of asthma.

    PLoS Genet 2012, 8(11):e1003029. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Hirschhorn JN, Daly MJ: Genome-wide association studies for common diseases and complex traits.

    Nat Rev Genet 2005, 6(2):95-108. PubMed Abstract | Publisher Full Text OpenURL

  23. Kang MJ, Yu HS, Seo JH, Kim HY, Jung YH, Kim YJ, Kim HJ, Lee SY, Hong SJ: GSDMB/ORMDL3 variants contribute to asthma susceptibility and eosinophil-mediated bronchial hyperresponsiveness.

    Hum Immunol 2012, 73(9):954-959. PubMed Abstract | Publisher Full Text OpenURL

  24. Yu J, Kang MJ, Kim BJ, Kwon JW, Song YH, Choi WA, Shin YJ, Hong SJ: Polymorphisms in GSDMA and GSDMB are associated with asthma susceptibility, atopy and BHR.

    Pediatr Pulmonol 2011, 46(7):701-708. PubMed Abstract | Publisher Full Text OpenURL

  25. Granell R, Henderson AJ, Timpson N, St Pourcain B, Kemp JP, Ring SM, Ho K, Montgomery SB, Dermitzakis ET, Evans DM, et al.: Examination of the relationship between variation at 17q21 and childhood wheeze phenotypes.

    J Allergy Clin Immunol 2013, 131(3):685-694. PubMed Abstract | Publisher Full Text OpenURL

  26. HapMap: A haplotype map of the human genome.

    Nature 2005, 437(7063):1299-1320. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Browning SR: Missing data imputation and haplotype phase inference for genome-wide association studies.

    Hum Genet 2008, 124(5):439-450. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Saeki N, Kim DH, Usui T, Aoyagi K, Tatsuta T, Aoki K, Yanagihara K, Tamura M, Mizushima H, Sakamoto H, et al.: GASDERMIN, suppressed frequently in gastric cancer, is a target of LMO1 in TGF-beta-dependent apoptotic signalling.

    Oncogene 2007, 26(45):6488-6498. PubMed Abstract | Publisher Full Text OpenURL

  29. Cohen L, X E, Tarsi J, Ramkumar T, Horiuchi TK, Cochran R, DeMartino S, Schechtman KB, Hussain I, Holtzman MJ, et al.: Epithelial cell proliferation contributes to airway remodeling in severe asthma.

    Am J Respir Crit Care Med 2007, 176(2):138-145. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Wu H, Romieu I, Sienra-Monge JJ, Li H, del Rio-Navarro BE, London SJ: Genetic variation in ORM1-like 3 (ORMDL3) and gasdermin-like (GSDML) and childhood asthma.

    Allergy 2009, 64(4):629-635. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Binia A, Khorasani N, Bhavsar PK, Adcock I, Brightling CE, Chung KF, Cookson WO, Moffatt MF: Chromosome 17q21 SNP and severe asthma.

    J Hum Genet 2011, 56(1):97-98. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Wan YI, Shrine NR, Soler Artigas M, Wain LV, Blakey JD, Moffatt MF, Bush A, Chung KF, Cookson WO, Strachan DP, et al.: Genome-wide association study to identify genetic determinants of severe asthma.

    Thorax 2012, 67(9):762-768. PubMed Abstract | Publisher Full Text OpenURL

  33. Borish L, Steinke JW: Interleukin-33 in asthma: how big of a role does it play?

    Curr Allergy Asthma Rep 2011, 11(1):7-11. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

Pre-publication history

The pre-publication history for this paper can be accessed here: