Genome-wide association study of asthma identifies RAD50-IL13 and HLA-DR/DQ regions

Article Outline

• Abstract
• Methods
  • Study subjects
  • Statistical analysis
• Results
• Discussion
• Acknowledgment
• Fig E1. 
• Fig E2. 
• Fig E3. 
• Fig E4. 
• Table E1. 
• Table E2. 
• Table E3. 
• References
• Copyright

Background

Asthma is a heterogeneous disease that is caused by the interaction of genetic susceptibility with environmental influences. Genome-wide association studies (GWASs) represent a powerful approach to investigate the association of DNA variants with disease susceptibility. To date, few GWASs for asthma have been reported.

Objectives

A GWAS was performed on a population of patients with severe or difficult-to-treat asthma to identify genes that are involved in the pathogenesis of asthma.

Methods

A total of 292,443 single nucleotide polymorphisms (SNPs) were tested for association with asthma in 473 The Epidemiology and Natural History of Asthma: Outcomes and Treatment Regimens (TENOR) cases and 1892 Illumina general population controls. Asthma-related quantitative traits (total serum IgE, FEV₁, forced vital capacity, and FEV₁/forced vital capacity) were also tested in identified candidate regions in 473 TENOR cases and 363 phenotyped controls without a history of asthma to analyze GWAS results further. Imputation was performed in identified candidate regions for analysis with denser SNP coverage.

Results

Multiple SNPs in the RAD50-IL13 region on chromosome 5q31.1 were associated with asthma: rs2244012 in intron 2 of RAD50 (P = 3.04E-07). The HLA-DR/DQ region on chromosome 6p21.3 was also associated with asthma: rs1063355 in the 3′ untranslated region of HLA-DQB1 (P = 9.55E-06). Imputation identified several significant SNPs in the T_H2 locus control region 3′ of RAD50. Imputation also identified a more significant SNP, rs3998159 (P = 1.45E-06), between HLA-DQB1 and HLA-DQA2.

Conclusion

This GWAS confirmed the important role of T_H2 cytokine and antigen presentation genes in asthma at a genome-wide level and the importance of additional investigation of these 2 regions to delineate their structural complexity and biologic function in the development of asthma.

Key words: Asthma, GWAS, RAD50, IL13, HLA-DQB1, TENOR

Abbreviations used: FVC, Forced vital capacity, GC, Genomic control, GWAS, Genome-wide association study, IBS, Identity-by-state, LCR, Locus control region, LD, Linkage disequilibrium, MAF, Minor allele frequency, QC, Quality control, SNP, Single nucleotide polymorphism, TENOR, The Epidemiology and Natural History of Asthma: Outcomes and Treatment Regimens

Asthma is a complex disease that is caused by the interaction of genetic susceptibility with environmental influences. Genome-wide linkage studies, candidate-gene association studies, and genome-wide association studies (GWASs) represent 3 major approaches to investigate the association between genetic variants and disease development.

Genome-wide linkage studies have consistently identified regions linked to asthma or asthma-related traits on chromosome 2q, 5q, 6p, 12q, and 13q.¹ The most highly replicated regions with obvious candidate genes are chromosome 5q31-33 (including IL5, IL13, IL4, CD14, and adrenergic β-2-receptor) and 6p21 (including lymphotoxin-α [or TNFB], TNF, major MHC-II, HLA-DQB1, and HLA-DRB1).² In addition, a recent meta-analysis of genome-wide linkage studies of asthma, bronchial hyperresponsiveness, positive allergen skin prick test, and total IgE identified overlapping regions for multiple phenotypes on chromosomes 5q and 6p as well as 3p and 7p.³ Unfortunately, genome-wide linkage studies can only identify genes with relative strong effects in broad regions that include many genes. Positional cloning studies have identified 6 genes for asthma: a disintegrin and metalloprotease domain 33 on chromosome 20p13,⁴ dipeptidyl-peptidase 10 on 2q14.1,⁵ PHD finger protein 11 on 13q14.11,⁶ neuropeptide S receptor 1 (or GPRA) on 7p14.3,⁷ MHC-I, G (HLA-G) on 6p21.3,⁸ and cytoplasmic FMR1 interacting protein 2 on 5q33.3.⁹

Candidate-gene association studies have identified more than 100 genes for asthma and asthma-related traits.², ¹⁰, ¹¹ Although candidate-gene association studies have identified many genes, only a few have been replicated extensively. Thus, only 14 genes including genes on 5q and 6p (adrenergic β-2-receptor, IL-4 receptor, HLA-DRB1, IL13, CD14, TNF, membrane-spanning 4-domains, subfamily A, member 2 [or FCER1B], IL4, a disintegrin and metalloprotease domain 33, signal transducer and activator of transcription 6, IL4 induced, IL10, HLA-DQB1, glutathione S-transferase π 1, and lymphotoxin-α) have been replicated in more than 20 independent studies.¹⁰ Even for highly replicated genes, replication might be a result of winner's bias and/or loose replication standard (gene as a unit and related phenotypes).

A GWAS is a hypothesis-free approach able to identify novel genes with mild/moderate effects and thus has become the best approach for studying association between genes and common disease phenotypes. To date, only 4 GWASs have been performed for asthma and asthma-related traits.¹² The first GWAS of childhood asthma identified ORM1-like 3 on chromosome 17q12.¹³ The second GWAS of serum YKL-40 levels identified chitinase 3-like 1 on 1q32.¹⁴ The third GWAS was for a related trait, total serum IgE levels, and the most significant single nucleotide polymorphisms (SNPs) are in the Fc fragment of IgE, high-affinity I, receptor for α polypeptide gene (FCER1A) on chromosome 1q23, and the second highest region observed was RAD50 on 5q31.¹⁵ The fourth GWAS of childhood asthma indicated phosphodiesterase 4D, cyclic adenosine monophosphate–specific (phosphodiesterase E3 dunce homolog, Drosophila) (PDE4D) on chromosome 5q12.¹⁶

In this study, we performed a GWAS of asthma in The Epidemiology and Natural History of Asthma: Outcomes and Treatment Regimens (TENOR) population of severe or difficult to treat asthmatics to search for novel genes and to confirm previously identified genes involved in asthma. The purpose of the TENOR study was to investigate the natural history of asthma in a large cohort of well characterized patients with severe or difficult to treat asthma; no treatment intervention was involved, and patients continued to be treated by their asthma specialists.¹⁷, ¹⁸, ¹⁹

Back to Article Outline

Methods 

Study subjects 

The TENOR study was a multicenter observational and longitudinal cohort study of 4756 patients with asthma described as “severe or difficult-to-treat” by their physicians, sponsored by Genentech and Novartis.¹⁷ Subjects were included if they had physician-characterized difficult-to-treat asthma and met additional criteria based on frequency of urgent care visits and/or the use of multiple controller medications. The clinical sites from the original TENOR study were contacted and invited to participate in this study. Sites that agreed were mailed Oragene DNA saliva collection kits (DNA Genotek, Inc, Ontario, Canada), labeled with the TENOR participant identification number. Sites then mailed the kits to participating individuals, who sent their collected samples to the Center for Human Genomics at Wake Forest University School of Medicine. This process was required to maintain anonymity between investigators at Wake Forest University and the study participants. Unfortunately, the TENOR study had ended (end of 2004) before this project started, so it was difficult to recontact participants. A total of 607 samples had sufficient DNA for successful SNP genotyping. Table I shows the demographic data for the TENOR cases and the 2 control populations. The TENOR patients with asthma who were genotyped were similar in characteristics to the larger TENOR cohort.

Table I. Demographics (means ± SDs) of subjects in TENOR and Illumina and phenotyped controls

Illumina controls were used for GWAS. The Wake Forest phenotyped controls were mainly recruited through the NHLBI Collaborative Study on the Genetics of Asthma and the NHLBI Severe Asthma Research Program and were genotyped as a subset of the NHLBI STAMPEED study.

General population controls were obtained by using the Illumina iControlDB client (www.illumina.com) to download genotypes for 3294 white individuals with genotype data available from any of the 3 available HumanHap550 k products (v1, v3, and -2v3). As shown in Table I, only age and sex data are available. Additional control samples for asthma-related quantitative traits were obtained from a separate GWAS for asthma. These 363 phenotyped controls had no personal or family history of asthma and had normal pulmonary function including lack of bronchial hyperresponsiveness or bronchodilator reversibility. Testing also included measures of atopy including total serum IgE levels (Table I). HapMap samples (N = 262) to be used for genetic ancestry check were also downloaded from the iControlDB database (Illumina, Inc) after selecting the HumanHap300_v1 genotyping product.

DNA was isolated by using the protocol described by DNA Genotek, and SNP genotyping was performed by using the Illumina HumanCNV370 BeadChip. The samples were clustered by first applying Illumina's cluster definition, removing samples with call rates less than 0.90, and then reclustering using the samples themselves.

Statistical analysis 

PLINK (version 1.06; http://pngu.mgh.harvard.edu/purcell/plink/)²⁰ was the main software used to perform statistical analysis unless otherwise stated.

Quality control (QC) was applied to cases and controls separately because they were genotyped by using slightly different Illumina products. Genetic ancestry of the TENOR cases was determined using the HapMap 300 k dataset as a reference. Fixed 3 groups clustering and pairwise population concordance of 1.0E-05 based on identity-by-state (IBS) were used to cross-validate ethnic group identity. Subjects were removed if they (1) were not of European white descent, (2) had low genotyping call rates (<95%), (3) were discrepant or ambiguous for genetic sex (heterozygous haploid genotype percentage ≥0.01 or X chromosome homorozygosity F ≥0.9), (4) failed the cryptic relatedness check (PI_HAT > 0.125), or (5) were detected as an outlier (>6 SD for the first or second principal component). After subjects meeting these criteria were deleted, SNPs were deleted if the call rates were low (95%) or were inconsistent with Hardy-Weinberg equilibrium (P < 10E-04). QC was then applied on the subjects and SNPs of merged case-control dataset as done separately. SNPs were also deleted if the minor allele frequency (MAF) was less than 0.05 in cases and controls or the Hardy-Weinberg equilibrium P value was less than .01 in controls only.

Asthma susceptibility was analyzed by comparing the non-Hispanic white TENOR cases to the general population Illumina controls. To reduce population stratification, 4 controls were matched with every 1 case based on pairwise IBS. Principal components were generated by using principal components analysis in EIGENSTRAT (version 3.0; http://genepath.med.harvard.edu/∼reich/Software.htm).²¹ Sex, age, and significant principal components were used as covariates in the logistic additive model. Genomic control (GC) was applied on P values to reduce population stratification further.²² A linear model was analyzed in GWAS-identified candidate regions in 473 TENOR cases and 363 phenotyped controls for asthma-related quantitative traits (total serum IgE, % predicted FEV₁, forced vital capacity [FVC], and FEV₁/FVC).

Haploview (http://www.broad.mit.edu/mpg/haploview/) was used to generate linkage disequilibrium plots.²³ Ninety-five percent CIs on D' were used to define blocks.²⁴ SNAP (version 2.0; http://www.broad.mit.edu/mpg/snap/) was used to generate the association plots.²⁵ Imputation was performed based on HapMap II CEU genotype data²⁶ by using MACH (version 1.0; http://www.sph.umich.edu/csg/abecasis/MaCH/index.html).²⁷ Association of candidate SNPs with nearby gene expression data in lymphocytes was performed based on the GENEVAR dataset (http://www.sanger.ac.uk/humgen/genevar/)²⁸ by using WGAViewer.²⁹

Back to Article Outline

Results 

A total of 607 TENOR cases were genotyped with the HumanCNV370 BeadChip. After removal of nonwhite samples (see this article's Fig E1 in the Online Repository at www.jacionline.org) and removal on the basis of the QC criteria, data from 474 patients with asthma were carried forward to analysis. Of the 3294 Illumina white controls downloaded from iControldb, 3,141 Illumina controls passed QC. After merging 474 TENOR cases with 3141 Illumina controls and evaluating the combined QC metrics, 473 cases and 3106 controls were retained. To reduce population stratification, 4 controls were matched with every 1 case on the basis of pairwise IBS; thus, 473 cases and 1892 Illumina controls were used for GWAS (see Table I for demographics and this article's Fig E2 in the Online Repository at www.jacionline.org). After QC analysis of the 318,075 common SNPs, 292,443 SNPs were retained for the GWAS.

The GWAS of asthma was performed on 292,443 SNPs of 473 TENOR cases and 1892 Illumina controls with sex, age, and significant principal components as covariates in the logistic additive model (Fig 1). GC was applied to P values to reduce population stratification (genomic inflation factor = 1.073 and 1.000 before and after adjustment; see this article's Fig E3 in the Online Repository at www.jacionline.org). In total, 248 SNPs had GC-adjusted P values ≤1.0E-03 (see this article's Table E1 in the Online Repository at www.jacionline.org). Focusing on SNPs with a GC-adjusted P values ≤1.0E-04 and at least 2 neighboring SNPs (±100 kb) with GC-adjusted P values ≤1.0E-03, six regions were identified: RAD50-IL13 on chromosome 5q31.1, HLA-DR/DQ on 6p21.32, low density lipoprotein-related protein 1B on 2q22.1-22.2, sorting nexin 10 on 7p15.2, carbonic anhydrase X on 17q21.33, and potassium inwardly-rectifying channel, subfamily J, member 2 on 17q24.3 (Fig 1; Table E1).

• 
  • View Large Image
  • Download to PowerPoint

• Fig 1. 

  Genome-wide association of 292,443 SNPs in 473 TENOR cases and 1892 Illumina controls. Color scale of the x-axis represents chromosomes. Negative logarithm–transformed GC-adjusted P values are shown on the y-axis.

The RAD50-IL13 region had the strongest evidence for association (Table II; Fig 2, A) with multiple SNPs in this region strongly associated with asthma susceptibility (rs2244012, rs6871536, and rs2897443 in RAD50 ranked highly as 1, 2, and 4, respectively, in this study). rs2244012 in intron 2 of RAD50 had an odds ratio of 1.64 (95% CI, 1.36-1.97; P = 3.04E-07; GC-adjusted P = 7.69E-07). Three SNPs in or near IL13 (rs2243204 [3′ downstream], rs20541 [Arg130Gln], and rs1295686 [intron 3]) were also associated with asthma (P < .001), but in weak LD (0.2 < r² < 0.3) with SNPs in RAD50 (Fig 2, A and B). rs2243300, which is ∼5 kb upstream of IL4, was weakly associated with asthma (P = .0032). Six SNPs downstream of IL5 were not associated with asthma, although they were in weak LD with SNPs in RAD50 (Table II; Fig 2, A and B). Four LD blocks were identified based on 95% CI of D' (Fig 2, B).²⁴ Blocks 1 and 2 were each composed of 3 SNPs downstream of IL5. Block 3 was composed of 4 SNPs in the intron of RAD50. Block 4 was composed of 2 SNPs in IL13.

Table II. Association results of 14 SNPs in T_H2 cytokine locus on chromosome 5

Alleles (M:m), Major allele:minor allele; OR, odds ratio.

Log total IgE, FEV₁/FVC, FEV₁, and FVC are P values of asthma-related quantitative traits.

• 
  • View Large Image
  • Download to PowerPoint

• Fig 2. 

  LD and association plot of 14 SNPs in the T_H2 cytokine locus. A, Association plot: negative logarithm–transformed P values (left) and recombination rate (right). B, LD plot: r² color scheme was used and labeled. 95% CIs on D' were used to set up blocks.

Linear model analysis was performed with the 14 SNPs of the T_H2 cytokine locus in 473 TENOR cases and 363 phenotyped controls for asthma-related quantitative traits (total serum IgE, FEV₁, FVC, and FEV₁/FVC; Table II). Multiple SNPs in each gene—RAD50, IL13, and IL4, but not IL5—showed significant association (P ≤ .05) with asthma-related quantitative traits.

The HLA-DR/DQ region (Table III; Fig 3, A) also showed consistent association with asthma. rs1063355 in the 3′ untranslated region of HLA-DQB1 had an odds ratio of 0.68 (95% CI, 0.58-0.81; P = 9.55E-06; GC-adjusted P = 1.93E-05). Ten of the 46 SNPs in the HLA-DR/DQ region had P ≤.001 (Table III; Fig 3, A). Multiple SNPs in or near butyrophilin-like 2, HLA-DRA, HLA-DRB1, HLA-DQB1, and HLA-DQA2, were strongly associated with asthma (P < 10E-04). LD is complicated in this region when considering all 46 SNPs (data not shown). One LD block composed of 3 SNPs upstream of HLA-DQB1 was formed on the basis of the 95% CI of D' of these 10 SNPs (Fig 3, B).

Table III. Association results of 10 of 46 SNPs (with P ≤ .001) in the HLA-DR/DQ region

Alleles (M:m), Major allele:minor allele; BTNL2, butyrophilin-like 2; OR, odds ratio.

• 
  • View Large Image
  • Download to PowerPoint

• Fig 3. 

  LD and association plot of 10 SNPs in the HLA-DR/DQ region. A, Association plot: negative logarithm–transformed P values (left) and recombination rate (right). B, LD plot: r² color scheme was used and labeled. 95% CIs on D' were used to set up blocks. Only 10 of 46 SNPs (with P ≤ .001) are shown.

Linear model analysis was performed with the 10 SNPs of the HLA-DR/DQ region for asthma-related quantitative traits (Table III). A single SNP, rs1063355, on HLA-DQB1 showed significant association with asthma-related quantitative traits (P = .01, .001, .007, and .05 for total serum IgE, FEV₁, FVC, and FEV₁/FVC, respectively).

Back to Article Outline

Discussion 

The highest associated SNP identified in this study was rs2244012 in intron 2 of RAD50 (P = 3.04E-07). In addition, evidence was observed for association with multiple SNPs in the RAD50-IL13 region for asthma susceptibility and asthma related quantitative traits. The protein encoded by RAD50 is involved in DNA double-strand break repair, and its expression level is constitutively low in most tissues; thus, it has no known function directly related to asthma, although the MER11-RAD50-NBS1 complex has been shown to be involved in somatic hypermutation and gene conversion of immunoglobulin regions.³⁰ On the contrary, other genes (IL4, IL5, and IL13) in the T_H2 cytokines locus are better candidates on the basis of their biologic functions. Three SNPs in IL13 in this study were associated with asthma. IL13 is critical to the pathogenesis of allergen-induced asthma and thus one of the most highly studied and replicated genes in both genome-wide linkage and candidate-gene association studies. rs20541 (Arg130Gln or IL13+4257GA), in the coding region of IL13, was also analyzed in this study and has been shown to be associated with asthma³¹ and total serum IgE levels.³² rs1800925 (IL13-1111CT), in the promoter region of IL13, has been shown to be associated with asthma³³ and total serum IgE levels.³⁴ In a GWAS with total serum IgE levels, 4 SNPs in RAD50 (rs2706347, rs3798135, rs2040704, and rs7737470), have been identified (P < 10E-04).¹⁵ These 4 SNPs in RAD50 were in strong LD with rs1800925 (0.7 < r² < 0.8) and in weak LD with rs20541 (0.2 < r² < 0.3) in IL13.¹⁵ These results are consistent with the results of this study because many of the TENOR patients with asthma were recruited from allergists' offices and the population has increased IgE levels.¹⁸ Because the actual functional SNPs cannot be determined purely by their P values, it is difficult to dissect the association data of RAD50 from IL13 in this study or other genetic studies because of the degree of LD present in this chromosomal region.

In a transgenic mouse study, a T_H2 locus control region (LCR) was identified as the 25-kb fragment at the 3′ end of Rad50.³⁵ An LCR is defined experimentally as regulating the expression of linked genes in a copy number–dependent and tissue-specific manner. The T_H2 LCR is involved in the chromatin configuration to reorganize promoters of IL4, IL5, and IL13 in proximity and coregulation of T_H2 cytokine expression.³⁶ Seven Rad50 DNase I–hypersensitive sites (RHS1-7) were identified, where RHS4-7 formed the core of the LCR.³⁷ LCR-C (RHS7) and LCR-B (RHS6) were possible T_H2 cytokine expression enhancers; LCR-A (RHS6) and LCR-O (RHS5) were likely insulators.³⁸ RHS7 is essential for T_H2 cytokine expression by showing T_H2 specific demethylation after allergen stimulation and intrachromosomal interactions between LCR and the promoters of T_H2 cytokines.³⁹ Furthermore, RHS6, Rad50 promoter (RHS2), and IL5 promoter interacted with IFN-γ (Ifng) on a different chromosome, which suggests an interchromosomal regulation of the expression of T_H1/T_H2 cytokines.⁴⁰ Although all these experiments were done in mouse, the RAD50 sequence is highly conserved in the LCR between human being and mouse. With imputation, multiple significant SNPs were found in the LCR (see this article's Table E2 in the Online Repository at www.jacionline.org): rs3798135 (P = 1.49E-06, in RHS5/LCR-O), rs12653750 (P = 1.49E-06, in RHS6/LCR-A), rs2040704 (P = 1.33E-06, in RHS6/LCR-B), and rs2240032 (P = 6.68E-06, in RHS7/LCR-C). The association of rs2244012 with the expression levels of IL13 in lymphocytes from white adults based on the GENEVAR dataset was not significant (P = 0.176), but may be a result of small sample size.

Because both a previous GWAS for total serum IgE levels and our GWAS of asthma identified RAD50, it appears to be a new candidate gene for asthma. Although it is still possible the signal from RAD50 is purely a result of its LD with the promoter of IL13, RAD50 deserves to be carefully studied when considering T_H2 cytokine locus.

HLA-DR/DQ also showed consistent association with asthma—for example, rs1063355 in the 3′ untranslated region of HLA-DQB1 (P = 9.55E-06), rs2239804 in intron of HLA-DRA (P = 2.80E-05), and rs2516049 5′ upstream of HLA-DRB1 (P = 2.62E-05). HLA-DR/DQ is part of the HLA class II region, which is one of the most gene/variant-dense regions in the human genome and is associated with many diseases.⁴¹ HLA-DQB1 and HLA-DRB1 have been shown to be associated with asthma in multiple independent studies.⁴², ⁴³, ⁴⁴ Genetic variants in the HLA-DR/DQ region have also been shown to be highly associated with HLA-DR/DQ gene expression, indicating that the association of HLA-DR/DQ with disease might be a result of gene expression levels in addition to antigen recognition.⁴⁵, ⁴⁶ The association of rs2516049 with asthma in our study and with the expression levels of HLA-DRB1 (P = 1.25E-04) in lymphocytes from white adults based on GENEVAR dataset indicated that the variant might function through expression level changes (see this article's Fig E4 in the Online Repository at www.jacionline.org).²⁸, ²⁹ Imputation identified a SNP with a more significant P value, rs3998159 (P = 1.45E-06), between HLA-DQB1 and HLA-DQA2 (see this article's Table E3 in the Online Repository at www.jacionline.org). It is difficult to determine the functional genes/SNPs in the HLA-DR/DQ region in our study because of the complicated LD pattern in this region. The long-range LD and haplotype analysis based on the MHC Haplotype Project may solve the issue.⁴⁷

Using a GWAS approach, this study is the first to confirm the association of RAD50-IL13 and HLA-DR/DQ regions with asthma susceptibility, regions that have been identified by multiple candidate-gene association studies and 1 genome-wide association study on total serum IgE levels. Our results weakly replicated the findings of the other GWAS: ORM1-like 3 and gasdermin B (GSDML; rs7216389) with asthma (P = .057); FCER1A (rs2251746) with total serum IgE (P = .040); and chitinase 3-like 1 (rs880633) with FEV₁ (P = .003), FVC (P = .031), and FEV₁/FVC (P = .040). rs1588265 (P = .507) and rs1544791 (P = .678) in PDE4D with asthma were not replicated. GWAS of total serum IgE by Weidinger et al¹⁵ identified several SNPs in RAD50 (P < 10E-04). In our study, the most significant SNP in RAD50 for total serum IgE is rs6871536 (P = 2.61E-03). The geometric mean of total serum IgE in the study by Weidinger et al¹⁵ is 42.41 (95% CI, 39.56-45.47). In our study, the geometric mean of total serum IgE is higher, 48.94 (95% CI, 43.04-55.65). The difference in the total serum IgE distribution and the relatively small sample size in our study may lead to the difference of significant levels between these 2 studies.

The potential for false-negative results could not be avoided in this study because of the relatively small sample size (473 cases), which may also be the reason that although significance levels of 10^-7 (E-07) were observed, no SNP reached the Bonferroni-adjusted multiple test criterion (P = .05/292,443 = 1.71E-07). However, evidence for multiple SNPs was observed in our results in this comprehensively phenotyped relatively homogeneous cohort of patients with difficult-to-treat asthma from the larger TENOR study. Our control datasets (general population and phenotyped controls) both have some limitations. They were both significantly younger (Table I) than TENOR cases, making our results a little conservative because some controls might become asthma cases in the future. Genotyping confirmation and fine-mapping of candidate regions were impossible because the Illumina controls were from a public database, but our approach compensated for this by using imputation. Population stratification was relatively strong between TENOR cases and Illumina 550k controls.

This GWAS confirmed the important role of T_H2 cytokine and antigen presentation genes in asthma at a genome-wide level. Furthermore, these findings will stimulate more comprehensive research (eg, resequencing, long-range LD, epistasis, epigenetics, copy number variant, and function) on these 2 regions because of their functional importance and structural complexity.

Back to Article Outline

We thank Dr Elizabeth J. Ampleford for analytical assistance. We also acknowledge the TENOR/SARP/CSGA/STAMPEED Study Group and the TENOR/SARP/CSGA/STAMPEED participants who contributed DNA samples.

Back to Article Outline

Fig E1. 

• 
  • View Large Image
  • Download to PowerPoint

• Cluster plot of 607 TENOR cases and 262 HapMap samples. PC1 and PC2 are the first and second principle component, respectively.

Back to Article Outline

Fig E2. 

• 
  • View Large Image
  • Download to PowerPoint

• Principle components plot of 473 TENOR cases and 1892 Illumina controls. PC1 and PC2 are the first and second principle component, respectively.

Back to Article Outline

Fig E3. 

• 
  • View Large Image
  • Download to PowerPoint

• GWAS qq-plot of 292,443 SNPs of 473 TENOR cases and 1,892 Illumina controls. Negative logarithm–transformed expected P values are shown on the x-axis. Negative logarithm–transformed observed P values are shown on the y-axis.

Back to Article Outline

Fig E4. 

• 
  • View Large Image
  • Download to PowerPoint

• Association of rs2516049 with the expression of HLA-DRB1 in lymphocytes based on GENEVAR dataset. X-axis represents copy numbers of minor allele of rs2516049. Y-axis represents the expression levels of HLA-DRB1 in lymphocytes.

Back to Article Outline

Table E1. 

Association results of 248 SNPs with GC-adjusted P values ≤1.0E-03

UTR, Untranslated region.

Back to Article Outline

Table E2. 

Association results of 37 imputed SNPs in RAD50

Boldface SNPs were genotyped SNPs; others were imputed.

Back to Article Outline

Table E3. 

Forty-one association results of imputed SNPs in HLA-DR/DQ region (P ≤ 1.0E-04)

BTNL2, Butyrophilin-like 2; UTR, Untranslated region.

Boldface SNPs were genotyped SNPs; others were imputed.

Back to Article Outline

References 

1. Hoffjan S, Ober C. Present status on the genetic studies of asthma. Curr Opin Immunol. 2002;14:709–717
   • View In Article
   • MEDLINE
   • CrossRef
2. Ober C, Hoffjan S. Asthma genetics 2006: the long and winding road to gene discovery. Genes Immun. 2006;7:95–100
   • View In Article
   • MEDLINE
   • CrossRef
3. Denham S, Koppelman GH, Blakey J, Wjst M, Ferreira MA, Hall IP, et al. Meta-analysis of genome-wide linkage studies of asthma and related traits. Respir Res. 2008;9:38
   • View In Article
   • CrossRef
4. Van Eerdewegh P, Little RD, Dupuis J, Del Mastro RG, Falls K, Simon J, et al. Association of the ADAM33 gene with asthma and bronchial hyperresponsiveness. Nature. 2002;418:426–430
   • View In Article
   • MEDLINE
   • CrossRef
5. Allen M, Heinzmann A, Noguchi E, Abecasis G, Broxholme J, Ponting CP, et al. Positional cloning of a novel gene influencing asthma from chromosome 2q14. Nat Genet. 2003;35:258–263
   • View In Article
   • MEDLINE
   • CrossRef
6. Zhang Y, Leaves NI, Anderson GG, Ponting CP, Broxholme J, Holt R, et al. Positional cloning of a quantitative trait locus on chromosome 13q14 that influences immunoglobulin E levels and asthma. Nat Genet. 2003;34:181–186
   • View In Article
   • MEDLINE
   • CrossRef
7. Laitinen T, Polvi A, Rydman P, Vendelin J, Pulkkinen V, Salmikangas P, et al. Characterization of a common susceptibility locus for asthma-related traits. Science. 2004;304:300–304
   • View In Article
   • CrossRef
8. Nicolae D, Cox NJ, Lester LA, Schneider D, Tan Z, Billstrand C, et al. Fine mapping and positional candidate studies identify HLA-G as an asthma susceptibility gene on chromosome 6p21. Am J Hum Genet. 2005;76:349–357
   • View In Article
   • MEDLINE
   • CrossRef
9. Noguchi E, Yokouchi Y, Zhang J, Shibuya K, Shibuya A, Bannai M, et al. Positional identification of an asthma susceptibility gene on human chromosome 5q33. Am J Respir Crit Care Med. 2005;172:183–188
   • View In Article
   • MEDLINE
   • CrossRef
10. Vercelli D. Discovering susceptibility genes for asthma and allergy. Nat Rev Immunol. 2008;8:169–182
    • View In Article
11. Zhang J, Pare PD, Sandford AJ. Recent advances in asthma genetics. Respir Res. 2008;9:4
    • View In Article
    • CrossRef
12. Willis-Owen SA, Cookson WO, Moffatt MF. Genome-wide association studies in the genetics of asthma. Curr Allergy Asthma Rep. 2009;9:3–9
    • View In Article
    • CrossRef
13. Moffatt MF, Kabesch M, Liang L, Dixon AL, Strachan D, Heath S, et al. Genetic variants regulating ORMDL3 expression contribute to the risk of childhood asthma. Nature. 2007;448:470–473
    • View In Article
14. Ober C, Tan Z, Sun Y, Possick JD, Pan L, Nicolae R, et al. Effect of variation in CHI3L1 on serum YKL-40 level, risk of asthma, and lung function. N Engl J Med. 2008;358:1682–1691
    • View In Article
    • CrossRef
15. Weidinger S, Gieger C, Rodriguez E, Baurecht H, Mempel M, Klopp N, et al. Genome-wide scan on total serum IgE levels identifies FCER1A as novel susceptibility locus. PLoS Genet. 2008;4:e1000166
    • View In Article
16. Himes BE, Hunninghake GM, Baurley JW, Rafaels NM, Sleiman P, Strachan DP, et al. Genome-wide association analysis identifies PDE4D as an asthma-susceptibility gene. Am J Hum Genet. 2009;84:581–593
    • View In Article
    • CrossRef
17. Dolan CM, Fraher KE, Bleecker ER, Borish L, Chipps B, Hayden ML, et al. Design and baseline characteristics of the epidemiology and natural history of asthma: Outcomes and Treatment Regimens (TENOR) study: a large cohort of patients with severe or difficult-to-treat asthma. Ann Allergy Asthma Immunol. 2004;92:32–39
    • View In Article
    • Abstract
    • Full-Text PDF (199 KB)
    • CrossRef
18. Borish L, Chipps B, Deniz Y, Gujrathi S, Zheng B, Dolan CM. Total serum IgE levels in a large cohort of patients with severe or difficult-to-treat asthma. Ann Allergy Asthma Immunol. 2005;95:247–253
    • View In Article
    • Abstract
    • Full-Text PDF (348 KB)
    • CrossRef
19. Miller MK, Johnson C, Miller DP, Deniz Y, Bleecker ER, Wenzel SE. Severity assessment in asthma: an evolving concept. J Allergy Clin Immunol. 2005;116:990–995
    • View In Article
    • Abstract
    • Full Text
    • Full-Text PDF (186 KB)
    • CrossRef
20. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–575
    • View In Article
    • CrossRef
21. Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D. Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006;38:904–909
    • View In Article
    • MEDLINE
    • CrossRef
22. Clayton DG, Walker NM, Smyth DJ, Pask R, Cooper JD, Maier LM, et al. Population structure, differential bias and genomic control in a large-scale, case-control association study. Nat Genet. 2005;37:1243–1246
    • View In Article
    • MEDLINE
    • CrossRef
23. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–265
    • View In Article
    • MEDLINE
    • CrossRef
24. Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, et al. The structure of haplotype blocks in the human genome. Science. 2002;296:2225–2229
    • View In Article
    • CrossRef
25. Johnson AD, Handsaker RE, Pulit SL, Nizzari MM, O'Donnell CJ, de Bakker PI. SNAP: a web-based tool for identification and annotation of proxy SNPs using HapMap. Bioinformatics. 2008;24:2938–2939
    • View In Article
    • CrossRef
26. Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, Gibbs RA, et al. A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007;449:851–861
    • View In Article
    • CrossRef
27. Li Y, Abecasis GR. Mach 1.0: rapid haplotype reconstruction and missing genotype inference. Am J Hum Genet. 2006;S79:2290
    • View In Article
28. Stranger BE, Forrest MS, Dunning M, Ingle CE, Beazley C, Thorne N, et al. Relative impact of nucleotide and copy number variation on gene expression phenotypes. Science. 2007;315:848–853
    • View In Article
    • CrossRef
29. Ge D, Zhang K, Need AC, Martin O, Fellay J, Urban TJ, et al. WGAViewer: software for genomic annotation of whole genome association studies. Genome Res. 2008;18:640–643
    • View In Article
    • CrossRef
30. Yabuki M, Fujii MM, Maizels N. The MRE11-RAD50-NBS1 complex accelerates somatic hypermutation and gene conversion of immunoglobulin variable regions. Nat Immunol. 2005;6:730–736
    • View In Article
    • MEDLINE
    • CrossRef
31. Heinzmann A, Mao XQ, Akaiwa M, Kreomer RT, Gao PS, Ohshima K, et al. Genetic variants of IL-13 signalling and human asthma and atopy. Hum Mol Genet. 2000;9:549–559
    • View In Article
    • MEDLINE
    • CrossRef
32. Graves PE, Kabesch M, Halonen M, Holberg CJ, Baldini M, Fritzsch C, et al. A cluster of seven tightly linked polymorphisms in the IL-13 gene is associated with total serum IgE levels in three populations of white children. J Allergy Clin Immunol. 2000;105:506–513
    • View In Article
    • Abstract
    • Full Text
    • Full-Text PDF (83 KB)
    • CrossRef
33. van der Pouw Kraan TC, van Veen A, Boeije LC, van Tuyl SA, de Groot ER, Stapel SO, et al. An IL-13 promoter polymorphism associated with increased risk of allergic asthma. Genes Immun. 1999;1:61–65
    • View In Article
    • MEDLINE
    • CrossRef
34. Liu X, Beaty TH, Deindl P, Huang SK, Lau S, Sommerfeld C, et al. Associations between total serum IgE levels and the 6 potentially functional variants within the genes IL4, IL13, and IL4RA in German children: the German Multicenter Atopy Study. J Allergy Clin Immunol. 2003;112:382–388
    • View In Article
    • Abstract
    • Full Text
    • Full-Text PDF (162 KB)
    • CrossRef
35. Lee GR, Fields PE, Griffin TJ, Flavell RA. Regulation of the Th2 cytokine locus by a locus control region. Immunity. 2003;19:145–153
    • View In Article
    • MEDLINE
    • CrossRef
36. Spilianakis CG, Flavell RA. Long-range intrachromosomal interactions in the T helper type 2 cytokine locus. Nat Immunol. 2004;5:1017–1027
    • View In Article
    • MEDLINE
    • CrossRef
37. Fields PE, Lee GR, Kim ST, Bartsevich VV, Flavell RA. Th2-specific chromatin remodeling and enhancer activity in the Th2 cytokine locus control region. Immunity. 2004;21:865–876
    • View In Article
    • MEDLINE
    • CrossRef
38. Lee DU, Rao A. Molecular analysis of a locus control region in the T helper 2 cytokine gene cluster: a target for STAT6 but not GATA3. Proc Natl Acad Sci U S A. 2004;101:16010–16015
    • View In Article
    • MEDLINE
    • CrossRef
39. Lee GR, Spilianakis CG, Flavell RA. Hypersensitive site 7 of the TH2 locus control region is essential for expressing TH2 cytokine genes and for long-range intrachromosomal interactions. Nat Immunol. 2005;6:42–48
    • View In Article
    • MEDLINE
    • CrossRef
40. Spilianakis CG, Lalioti MD, Town T, Lee GR, Flavell RA. Interchromosomal associations between alternatively expressed loci. Nature. 2005;435:637–645
    • View In Article
    • CrossRef
41. Shiina T, Hosomichi K, Inoko H, Kulski JK. The HLA genomic loci map: expression, interaction, diversity and disease. J Hum Genet. 2009;54:15–39
    • View In Article
    • CrossRef
42. Torio A, Sanchez-Guerrero I, Muro M, Villar LM, Minguela A, Marin L, et al. HLA class II genotypic frequencies in atopic asthma: association of DRB1∗01-DQB1∗0501 genotype with Artemisia vulgaris allergic asthma. Hum Immunol. 2003;64:811–815
    • View In Article
    • MEDLINE
    • CrossRef
43. Lin YC, Lu CC, Su HJ, Shen CY, Lei HY, Guo YL. The association between tumor necrosis factor, HLA-DR alleles, and IgE-mediated asthma in Taiwanese adolescents. Allergy. 2002;57:831–834
    • View In Article
    • MEDLINE
    • CrossRef
44. Cho SH, Kim YK, Oh HB, Jung JW, Son JW, Lee MH, et al. Association of HLA-DRB1(∗)07 and DRB1(∗)04 to citrus red mite (Panonychus citri) and house dust mite sensitive asthma. Clin Exp Allergy. 2000;30:1568–1575
    • View In Article
    • MEDLINE
    • CrossRef
45. Dixon AL, Liang L, Moffatt MF, Chen W, Heath S, Wong KC, et al. A genome-wide association study of global gene expression. Nat Genet. 2007;39:1202–1207
    • View In Article
    • CrossRef
46. Beaty JS, West KA, Nepom GT. Functional effects of a natural polymorphism in the transcriptional regulatory sequence of HLA-DQB1. Mol Cell Biol. 1995;15:4771–4782
    • View In Article
    • MEDLINE
47. Horton R, Gibson R, Coggill P, Miretti M, Allcock RJ, Almeida J, et al. Variation analysis and gene annotation of eight MHC haplotypes: the MHC Haplotype Project. Immunogenetics. 2008;60:1–18
    • View In Article
    • CrossRef