Search tips
Search criteria 


Logo of plosonePLoS OneView this ArticleSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)
PLoS One. 2013; 8(2): e56179.
Published online 2013 February 14. doi:  10.1371/journal.pone.0056179
PMCID: PMC3572953

Integration of Mouse and Human Genome-Wide Association Data Identifies KCNIP4 as an Asthma Gene

Huiping Zhang, Editor


Asthma is a common chronic respiratory disease characterized by airway hyperresponsiveness (AHR). The genetics of asthma have been widely studied in mouse and human, and homologous genomic regions have been associated with mouse AHR and human asthma-related phenotypes. Our goal was to identify asthma-related genes by integrating AHR associations in mouse with human genome-wide association study (GWAS) data. We used Efficient Mixed Model Association (EMMA) analysis to conduct a GWAS of baseline AHR measures from males and females of 31 mouse strains. Genes near or containing SNPs with EMMA p-values <0.001 were selected for further study in human GWAS. The results of the previously reported EVE consortium asthma GWAS meta-analysis consisting of 12,958 diverse North American subjects from 9 study centers were used to select a subset of homologous genes with evidence of association with asthma in humans. Following validation attempts in three human asthma GWAS (i.e., Sepracor/LOCCS/LODO/Illumina, GABRIEL, DAG) and two human AHR GWAS (i.e., SHARP, DAG), the Kv channel interacting protein 4 (KCNIP4) gene was identified as nominally associated with both asthma and AHR at a gene- and SNP-level. In EVE, the smallest KCNIP4 association was at rs6833065 (P-value 2.9e-04), while the strongest associations for Sepracor/LOCCS/LODO/Illumina, GABRIEL, DAG were 1.5e-03, 1.0e-03, 3.1e-03 at rs7664617, rs4697177, rs4696975, respectively. At a SNP level, the strongest association across all asthma GWAS was at rs4697177 (P-value 1.1e-04). The smallest P-values for association with AHR were 2.3e-03 at rs11947661 in SHARP and 2.1e-03 at rs402802 in DAG. Functional studies are required to validate the potential involvement of KCNIP4 in modulating asthma susceptibility and/or AHR. Our results suggest that a useful approach to identify genes associated with human asthma is to leverage mouse AHR association data.


Asthma is a common chronic respiratory disease with a rise in prevalence over the past decades, affecting over 25 million Americans and 300 million people world-wide [1], [2]. Asthma is characterized by airway hyperresponsiveness (AHR), a trait distinguished by increased airway smooth muscle contractility in response to certain exposures. The mouse is commonly used to model the genetics of human diseases because mice are physiologically similar to humans, the mouse genome has been sequenced, and many tools exist that allow for direct testing of genetic alterations in mice [3], [4]. Although attempts to recreate all of the features of human asthma in mice have not been successful [5], many of the allergic asthmatic responses in mouse are similar to the responses observed in humans [6], [7]. The most common asthma phenotype studied in the mouse is AHR, and several regions that have been associated with mouse AHR are homologous with genomic regions linked with asthma-related phenotypes in human cohorts. For example, two linkage studies [8], [9] in mouse identified AHR quantitative trait loci (QTL) on chromosome 7 that are homologous with human AHR QTL on chromosome 19q observed in Hutterites and Chinese individuals [10], [11]. Mouse AHR QTL observed on chromosome 17 in three studies [9], [12], [13] overlap with human chromosome 6p QTLs identified in four studies of asthma and allergic phenotypes [14], [15], [16], [17]. Thus, it is likely that mouse and human share genetic variants that predispose both species to asthma-related phenotypes, and mouse AHR is an appropriate phenotype to identify some of these shared genetic variants.

Genome-wide association analysis of mouse inbred strains can provide a cost effective complement to traditional QTL methods [18], [19]. Such methods leverage the fact that laboratory mouse strains, whose lineage is well known, are nearly isogenic due to their inbred origin [20] and therefore, share large haplotype blocks that are identical by descent. Assuming that phenotypic differences observed between mouse strains are due to different variants inherited within ancestral haplotypes, regions associated with a trait can be narrowed by comparison of haplotypes across many strains. The ability to narrow such regions efficiently is made possible by genetic resources and bioinformatics tools, including dense SNP sets for mouse strains [21], [22] and association algorithms [18], [19]. Studies using mouse strain surveys have been successful in identifying genes relevant to human disease, such as HDL levels in blood [18]. A previous in silico QTL study of baseline AHR in 36 female mouse strains identified eight possible regions that could underlie differences in AHR among strains [23]. The results of this study have not been extensively validated or evaluated in humans.

The genetic basis of asthma has been widely studied in humans [24]. Most recently, multi-center, multi-cohort genome-wide association studies (GWAS) of asthma have been carried out in Europeans and diverse North American populations [25], [26]. These and other GWAS have begun to identify loci (e.g. IKZF3-ZPBP2-GSDMB-ORMDL3 locus, HLA-DQ, IL1RL1, IL18RL1, IL33, TSLP, SLC22A5, SMAD3, and RORA) that are consistently associated with asthma at statistical thresholds that leave little doubt the results are truly significant. However, it is likely that some of the findings that do not meet genome-wide significance levels in these GWAS represent true associations that are biologically relevant for asthma. The question of how to distinguish true associations among the false positive ones is challenging. In this work, we address this challenge by using mouse baseline AHR association measures and searching among previously published nominally significant EVE Consortium GWAS meta-analysis results [26] to identify genes that may be associated with asthma. After attempting to confirm our results in three additional asthma GWAS and two AHR GWAS, we found that the Kv channel interacting protein 4 (KCNIP4) is likely to be related to asthma and AHR.


Figure 1 is an overview of our study design. Measures of baseline AHR in 31 mouse strains found that there was up to a 4.2-fold change in AHR slope across strains, with the smallest slope (i.e. 0.40 SEM 0.10) corresponding to the C57BL/6J strain and the largest slope (i.e. 1.67 SEM 0.20) corresponding to the KK/H1J strain [Figure 2]. The complete AHR phenotype results are freely available as project “Berndt1” in the Mouse Phenome Database of the Jackson Laboratory. The EMMA results for AHR slope did not contain any P-values that were significant after correction for multiple comparisons (i.e. <0.05/281,300 = 1.8e-07) [Figure 3]. Nonetheless, some genomic regions had P-values that were nominally significant and may indicate true associations with AHR slopes. The set of 227 mouse genes within 50 kb of SNPs with EMMA P-values <0.001 was selected for further study.

Figure 1
Study overview.
Figure 2
Baseline AHR slope phenotype distribution among 31 mouse strains.
Figure 3
Manhattan plot of mouse AHR slope EMMA results.

Of the 227 top mouse genes, 145 had human orthologs, according to the human-mouse HomoloGene map utilized. Results from the EVE GWAS based on all samples and specific to race/ethnic groups were obtained for all SNPs within 50 kb of the set of 145 genes. A subset of 35 genes had at least one SNP with a P-value <0.001 in the EVE European American, African and Caribbean American, Latino American and/or all subjects’ GWAS. The results in the vicinity of these 35 genes were explored more carefully to ensure that 1) the AHR slope phenotype distribution approximately matched the genotype distribution across mouse strains, 2) the mouse association results were based on SNPs with minor allele frequency >2/31, and 3) the association results were supported by a region of LD in both mouse and human. Two regions that included a total of five genes met these criteria: 1) KCNIP4, a 1220.19 kb-long gene that was covered (including 50 kb on either end) by 1,766 EVE combined meta-analysis SNPs, and 2) PDZD2/GOLPH3/MTMR12/ZFR, a 732 kb-long region that was covered (including 50 kb on either end) by 697 EVE European American analysis SNPs [Figure 4]. The smallest mouse EMMA p-values for these two regions were 9.1e-04 for Kcnip4 and 8.4e-04 for Pdzd2/Golph3/Mtmr12/Zfr. The genotype distribution for mouse SNPs with these P-values in the corresponding regions are shown in Figure 5. The EVE P-values supporting the KCNIP4 association were from the combined meta-analysis, and the lowest P-value (2.9e-04) was for SNP rs6833065. The EVE P-values supporting the PDZD2/GOLPH3/MTMR12/ZFR region of association were from the European American results, and the lowest P-value (5.5e-04) was for SNP rs17526969.

Figure 4
Top regions of association in mouse and human.
Figure 5
Mouse genotypes for SNPs with P-value <0.001 in the top two association regions.

We next attempted to replicate the asthma association findings in two human cohorts: Sepracor/LOCCS/LODO/Illumina and GABRIEL. We found that there was gene level replication for KCNIP4, and, considering the genomic interval spanning four genes, some regional replication for the PDZD2/GOLPH3/MTMR12/ZFR association [Figure 6]. Each asthma GWAS had a different top hit in each region, and only KCNIP4 showed evidence of SNP-level replication at a nominally significant threshold [Table 1]. Based on the replication results in GABRIEL and Sepracor/LOCCS/LODO/Illumina, we attempted to replicate the KCNIP4 findings in DAG. The lowest EVE P-value of 2.9e-04 at rs6833065 did not replicate in Sepracor/LOCCS/LODO/Illumina or DAG, and data for this SNP was not available in GABRIEL. The lowest GABRIEL association at rs469177 (P-value 1.0e-03) was nominally significant in EVE and DAG, but not Sepracor/LOCCS/LODO/Illumina. The lowest Sepracor/LOCCS/LODO/Illumina association at rs7664617 (P-value 1.5e-03) was nominally significant in EVE, but data for this SNP was not available in GABRIEL or DAG. The lowest combined P-value across the four asthma GWAS was 1.1e-04 for rs4697177, while the three SNPs with nominally significant P-values in 3/4 GWAS were rs4697177, rs6448072, and rs6856781. The allelic frequencies and odds ratios for the top KCNIP4 SNPs are shown in Table 2 and Table 3, respectively.

Figure 6
Top regions of association in two independent human asthma GWAS.
Table 1
Top KCNIP4 SNP Associations.
Table 2
Allele Frequencies of Top KCNIP4 SNP Associations.
Table 3
Effect Sizes of Top KCNIP4 SNP Associations.

Because the mouse phenotype measured, baseline AHR, is not a precise surrogate of human asthma, we also attempted to measure the association of top EMMA SNPs with another related phenotype for which GWAS data are available: human AHR among asthmatics. Of the top two regions, KCNIP4 and PDZD2/GOLPH3/MTMR12/ZFR, only KCNIP4 had gene-level replication in the SHARP AHR GWAS at a nominally significant level in the same regions as the asthma GWAS studies [Figure 7]. At a SNP level, there was no overlap between the strongest associations with asthma vs. AHR in either the SHARP or DAG AHR GWAS [Table 1].

Figure 7
Top regions of association in a human AHR GWAS.


Many risk-modifying loci have been identified via human GWAS studies for a wide range of complex diseases, including asthma [27]. Over time, the number of loci identified has grown, due in part to the greater statistical power of studies conducted. The primary way in which increased statistical power has been achieved is to conduct large-scale meta-analyses based on cohorts from many smaller studies. While the number of associations that pass genome-wide significance thresholds has increased, it is likely that some true associations that represent biologically important processes underlying diseases are present among the nominally significant ones. Approaches to identify true associations that are not among results meeting genome-wide significance levels in GWAS include attempts to replicate nominally significant findings in independent populations [26], [28], to conduct functional studies of top-ranked associations [29], [30], and to use gene-based [31] and pathway-based [32] methods. In this work, we have attempted to use a mouse genome-wide association analysis based on a strain survey of baseline AHR to identify asthma associations in humans.

The top two regions identified based on having mouse EMMA P-values <0.001 and EVE meta-analysis P-values <0.001 were KCNIP4 and PDZD2/GOLPH3/MTMR12/ZFR. While there was some evidence of gene- or regional-level replication for both sites among the GABRIEL and Sepracor/LOCCS/LODO/Illumina asthma GWAS studies, only the KCNIP4 gene had evidence of replication at individual SNPs. Additionally, only KCNIP4 SNPs were nominally associated with AHR among asthmatics in the SHARP AHR GWAS, indicating that this gene may be related to both baseline AHR (i.e. AHR in the absence of asthma) and asthma. Thus, we selected KCNIP4 as the final candidate for replication in DAG. We found that there was additional gene-level and SNP-level replication for KCNIP4 with asthma in DAG, and gene-level replication for KCNIP4 with AHR in DAG.

In addition to the nominal associations we observed, other studies have observed KCNIP4 SNP associations with asthma-related phenotypes. First, a small GWAS of toluene diisocyanate-induced asthma in Koreans (84 cases and 263 controls) found that KCNIP4 SNP rs4697192 had a P-value of 6.11e-05 for association under a recessive model [33]. Second, a family-based GWAS of change in lung function in response to glucocorticoid therapy among 118 asthmatics found that KCNIP4 SNP rs4282162 had a P-value of 0.028 under a recessive model of inheritance, but attempts to replicate this result in three independent populations were unsuccessful [30]. Neither of these SNPs was among those considered in Table 1. In a follow-up study of the first EVE meta-analysis GWAS, all SNPs with meta-analysis P-values <0.001 in the combined sample or <1e-04 in the race/ethnic-specific analyses were selected for replication [28]. Five KCNIP4 SNPs, including rs6833065 that is in Table 1, were among the SNPs genotyped in nine different replication populations of diverse origin (4 European American, 4 African American, 5 Latino American), comprised of 7202 cases, 6426 controls, and 507 case-parent trios. None of the KCNIP4 SNPs replicated the primary EVE findings. Taken together, the KCNIP4 association results from the current study and three previous ones suggest that if the nominally significant associations found represent a true relationship between this gene and asthma, then there is not a clear single functional variant underlying this relationship. Rather, there are weak effects observed at multiple sites that implicate the gene. It is possible that more thorough genotyping of KCNIP4 in the populations where single or few SNPs did not replicate primary findings would uncover regions of association that further support gene- or SNP-level replication. Future studies may also be able to investigate whether KCNIP4 variants are associated with AHR among non-asthmatics, as the mouse data found that KCNIP4 variants were associated with baseline AHR.

Based on the mouse and human association results, the Kv channel interacting protein 4 (KCNIP4) emerged as the gene most likely to be related to asthma and AHR. Current information for this gene was gathered with the AceView tool [34]. KCNIP4 is found on human chromosome 4, at 4p15.32, covering 1220.19 kb, from 21950417 to 20730227 (NCBI 37, August 2010). The protein encoded by KCNIP4 was first identified in a study searching for proteins that interacted with presenilins, which are linked to early-onset Alzheimer’s disease [35]. KCNIP4 was found to interact with the voltage-gated potassium channel subunit Kv4, as well as the presenilin PSEN2, and to be mostly expressed in brain. Based on its structure, KCNIP4 belongs to the recoverin branch of calcium binding proteins, characterized by having 4 EF-hand motifs, the first of which does not bind calcium. In addition to being expressed in brain, KCNIP4 is expressed in other tissues, including human skeletal muscle. Previous studies of peptides that bind to the EF-hands of calmodulin and other proteins that modulate intracellular calcium levels have found that these calcium sensors may be therapeutic targets for asthma [36]. For example, intratracheal pretreatment of guinea pigs with a calcium-like peptide (i.e. CALP1) that binds the EF-hands of calmodulin prior to ovalbumin challenge, prevented the development of AHR and attenuated the increased radical production by alveolar inflammatory cells of ovalbumin-challenged guinea pigs [37]. Thus, one possible mechanism by which KCNIP4 may influence asthma and AHR is via modulation of intracellular calcium.

Mouse is a species that is commonly used in the study of asthma genetics. Most genome-scale association studies of asthma in mouse have been linkage studies that identified QTL for AHR [8], [9], [12], [13]. More recently, in silico QTL approaches have been employed. In a study that characterized baseline AHR in female mice from 36 strains and then used haplotype association mapping (HAM) to identify chromosomal regions associated with AHR, eight peaks on six chromosomes (i.e., 3, 5, 8, 12, 13, 14) with nominal association p-values were suggested as interesting for additional study [23]. The peak on chromosome 13 at 34.7 Mb was found to coincide with previous mouse AHR QTL results and was explored further to identify 29 candidate genes. Our mouse results do not support seven of the eight peaks found in the previous paper, but they do support the previous association on chromosome 13: our EMMA results had a peak of 25 SNPs with P-value 5.2e-04 in the range chr13[ratio]36611438–36707895. Human associations for this region most strongly support the phenylalanyl-tRNA synthetase 2, mitochondrial (FARS2) gene as involved in asthma, with an EVE Latino American p-value of 6.4e-04 at rs9502304. However, the associations at this SNP were not significant among European Americans or African Americans, leading to a combined EVE meta-analysis p-value of 0.07 and exclusion as a top gene candidate in the current study. For the current study, some of the phenotype data reported by Leme et al [23] was utilized, but we extended this data by including AHR measures from male mice. Overall, we included data for 31 strains that had available AHR slope measures in mice of both genders and available genotype data. In addition to utilizing more phenotype data, we used a newer association method (i.e. EMMA) and more complete strain genotypic data than did the previous study. The primary advantage of EMMA vs. HAM is that EMMA takes the known genetic similarity between strains into account to control for population structure effects. The genome-wide genotype data for mouse strains used in the current study were for 281,300 SNPs, rather than 70,000 used in the previous report. Despite these improvements, our mouse association data still suffers from limitations. First, we are limited to a relatively small number of strains over which the range of AHR does not vary greatly. While the variability of AHR within strains seems to be significantly lower than that across strains, it may be the case that additional measures within strains would result in more precise phenotypes that have less within-strain variability and that increasing the number of strains phenotyped would provide a greater range of AHR variability across strains. Second, the phenotype studied, baseline AHR, may not be the best phenotype to characterize asthma susceptibility. For example, measures of AHR after mice are sensitized to ovalbumin may result in greater phenotype heterogeneity across strains and/or less heterogeneity within strains. Such limitations likely led the current mouse study to have insufficient statistical power to detect associations passing multiple comparisons correction. In addition to increasing the number of strains characterized and including additional phenotypes, future mouse genetic studies of asthma may be greatly improved with resources from the Collaborative Cross [38].

Translating mouse association results to the human genome is challenging at a genome-wide scale because of incomplete knowledge of the human/mouse genome map. We chose a gene-based analysis because the orthologous regions available in HomoloGene are the most accurately mapped ones, and also because such genes are most likely to be functionally relevant to phenotypes that are expressed in both mouse and human. It is possible that mouse genes excluded by the HomoloGene mapping are determinants of baseline AHR in mice. Thus, future studies with improved mouse/human genome maps may uncover additional regions of interest.

In summary, integration of genome-wide association results for a mouse strain survey of baseline AHR with human asthma and AHR GWAS results suggests that KCNIP4 is a gene related to both asthma and AHR. Functional studies are required to validate the potential involvement of this gene in asthma and to understand how specific SNPs may modulate asthma susceptibility and/or AHR.


Mouse Airway Responsiveness Measures

A survey of AHR was conducted in males of 33 mouse strains, to augment a survey of AHR that was previously conducted in females of 36 mouse strains [23]. AHR was measured by gathering a baseline resistance reading for 6–10 mice per strain and gender, followed by a reading after administration of saline, and sequentially increasing concentrations of methacholine (1, 3, 10, 30 mg/ml) administered via nebulizer through a tracheostomy [23]. AHR was quantified as the slope of the resistance vs. log-transformed methacholine concentration.

Mouse in silico Association Mapping

Of the 33 mouse strains with AHR phenotype data in males and the 36 mouse strains with AHR phenotype data in females, a subset of 31 strains overlapped and had genome-wide genotype data available. This subset of 31 mouse strains with AHR measures in both males and females and with available genotype data was used to conduct association analyses. The association of SNPs with AHR was measured using Efficient Mixed Model Association (EMMA) software [19], which uses a linear mixed model with a variance component using a kinship matrix that is based on the genetic similarity between strains to control for population structure effects. The genotype data for the mouse strains was part of the “Subspecific Origin and Haplotype Diversity in the Laboratory Mouse” project of the Jackson Laboratory available at [39]. Autosomal chromosome association results obtained with male and female mice were utilized. Monomorphic and singleton SNPs were excluded, and the final dataset contained association results for 281,300 SNPs.

Primary Human Asthma GWAS

The primary group of subjects consisted of participants of the EVE consortium, which comprises nine research teams in the USA with genome-wide association data for diverse asthma cohorts [26]. Briefly, 4,867 European Americans (1,486 cases, 1,539 controls, and 620 trios), 4,644 Latinos (606 cases, 792 controls, and 1,082 trios), and 3,447 African Americans (1,154 cases, 1,054 controls, and 413 trios) recruited through clinics and health systems in the U.S., Puerto Rico, Mexico, and Barbados were used in a meta-analysis of GWAS of asthma. The previously generated meta-analysis results for >2 million SNPs, based on genotyped or HapMap Phase 2 imputed data, for all subjects, as well as those limited to European American, African American, and Latino American groups were utilized. The test statistic to assess association of SNPs with asthma was based on a linear combination of normally distributed test statistics weighted by the square root of individual study sample sizes. Significance of this test statistic was assessed using standard normal approximations [26], [28].

Replication Human Asthma GWAS

1) Sepracor/LOCCS/LODO/Illumina

This cohort consisted of 531 non-Hispanic white cases with mild to severe asthma from three adult asthma populations: (1) a medication trial conducted by Sepracor, Inc., US [40], [41]; (2) the Leukotriene Modifier or Corticosteroid Salmeterol (LOCCS) study [42]; and (3) the Effectiveness of Low Dose Theophylline as an Add-on Treatment in Asthma (LODO) trial [43]. Cases were matched with 660 population controls obtained from Illumina’s IconDB resource ( using the genetic matching (GEM) algorithm [44]. Genome-wide genotyping of cases was performed on the Illumina 610 quad platform. Genome-wide genotyping of controls was performed on the Illumina HumanHap 550K v3 platform. Genotypes of SNPs of interest that were not captured by the genotype data were inferred using imputation with the Markov Chain Haplotyping software (MaCH) [45] based on HapMap Phase 2 Release 22 data [46]. The ratio of the empirically observed dosage variance to the expected (binomial) dosage variance for imputed SNPs utilized was greater than 0.3, indicating good quality of imputation. Dosage data was used to compute association statistics with PLINK [47]. The corresponding genetic inflation factor was 1.04, demonstrating minimal population stratification.


This European consortium-based GWAS of asthma consisted of 10,365 persons with physician-diagnosed asthma and 16,110 unaffected persons, all of whom were matched for ancestry [25]. Publicly available results for the random-effects pooled analysis that tested for association using the entire study population with asthma were used [Available at].

3) Dutch asthma GWAS (DAG)

This cohort was comprised of 920 DAG subjects with doctor-diagnosed asthma and documented AHR and 985 controls [48], [49]. Genotyping was performed using the Hapmap 317K platform or Illumina 370 Duo Chip. Tests of association were performed via logistic regression using PLINK.



A GWAS of AHR, quantified as the natural log of the dosage of methacholine causing a 20% drop in FEV1, was performed with 994 non-Hispanic white asthmatic subjects from the Childhood Asthma Management Program (CAMP) [50], and subsets of clinical trials within the Childhood Asthma Research and Education (CARE) network [51], and the Asthma Clinical Research Network (ACRN) [52] participating in the NHLBI SNP Health Association Resource Asthma Resource project (SHARP) [53]. Genome-wide SNP genotyping for SHARP subjects was performed by Affymetrix, Inc. (Santa Clara, CA) using the Affymetrix Genome-Wide Human SNP Array 6.0. Imputation of all SNPs available in HapMap Phase 2 Release 22 CEU data using MaCH [45] was performed. The primary GWAS was based on imputed SNPs with minor allele frequency (MAF) >0.05 and a ratio of empirically observed dosage variance to the expected (binomial) dosage variance greater than 0.3, indicating good quality of imputation. The association of SNPs with LnPC20 was measured with a linear regression model using dosage data as implemented in PLINK.

2) Dutch asthma GWAS (DAG)

This cohort was comprised of 650 DAG subjects with doctor-diagnosed asthma and documented AHR [48], [49]. All subjects had smoking history and steroid use data available at the time of the AHR test. Participants with a history of AHR but in remission during the test were excluded. Remission was defined as not on steroids and without 20% or greater fall in FEV1 during the AHR test. The AHR test was conducted using histamine or methacholine as a stimulus. AHR was quantified as the difference between FEV1 at baseline and at the dose step at which a 20% or greater FEV1 drop was achieved, divided by the dose of stimulant used (slope). Because two protocols were used, one with a 30-second tidal breathing method and a second with a 2-minute tidal breathing phase, the AHR slopes measured with the 30-second tidal breathing method were divided by 4, in order to compensate for the 4 times greater duration of administration of stimulus. Slope values were log-transformed so that they would follow a normal distribution. Genotyping was performed using the Hapmap 317K platform or Illumina 370 Duo Chip. Tests of AHR association were performed via linear regression, with smoking and inhaled/oral steroid use as covariates using PLINK.

Each study was approved by its respective Institutional Review Board, which ensured that all procedures followed were in accordance with the ethical standards of the responsible committee on human experimentation. Informed consent was obtained for all study participants.

Integration of Mouse and Human GWAS Results

Mouse gene coordinates were mapped to the human genome using the mouse comparative homology map data corresponding to Mouse Build 37.2, Human Build 37.2 and HomoloGene file from April 15, 2011 available at [54]. After selection of orthologous genes, the list of genes was restricted based on having an EVE meta-analysis P-value <0.001 in one of the race/ethnic-specific analyses or in the combined GWAS. Plots of human association results near specific genes were created using LocusZoom with the hg18/HapMap Phase II CEU GenomeBuild/LD Population [55]. Combined P-values with replication populations were computed using Fisher’s combined probability method [56] where hypothesis tests in replication populations had one-sided alternatives, based on the direction of the association in EVE, so that SNPs with association tests in opposite directions would not produce inappropriately small P-values. Analyses were carried out using the R computing environment [57].


We thank all CAMP subjects for their ongoing participation in this study. We acknowledge the CAMP investigators and research team, supported by the National Heart, Lung and Blood Institute (NHLBI), for collection of CAMP Genetic Ancillary Study data. All work on data collected from the CAMP Genetic Ancillary Study was conducted at the Channing Laboratory of the Brigham and Women’s Hospital under appropriate CAMP policies and human subject’s protections. We acknowledge the following GALA 1 investigators: Michael A. LeNoir, MD (Bay Area Pediatrics, Oakland, CA, USA), Harold J. Farber, MD (Section of Pulmonology, Baylor College of Medicine, Houston, TX, USA), Rajesh Kumar, MD, MPH (Division of Allergy and Immunology, Children's Memorial Hospital, Chicago, IL, USA), Pedro C. Avila, MD (Department of Medicine, Northwestern University, Chicago, IL, USA), Kelley Meade, MD (Children's Hospital of Oakland, Oakland, CA, USA), Denise Serebrisky, MD (Jacobi Medical Center and Department of Pediatrics, Albert Einstein College of Medicine, Bronx, NY, USA), Shannon Thyne, MD (Department of Medicine, University of California, San Francisco, San Francisco, CA, USA), William Rodriguez-Cintron, MD (Veterans Affairs Medical Center, San Juan, Puerto Rico), Jose R. Rodriguez-Santana, MD (Centro de Neumologia Pediatrica, San Juan, Puerto Rico), Luisa N. Borrell, DDS, PhD (Department of Health Sciences, Lehman College, City University of New York, New York, NY, USA).

Funding Statement

The CAMP Genetics Ancillary Study is supported by U01 HL075419, U01 HL65899, P01 HL083069, R01 HL086601, and T32 HL07427 from the NHLBI, National Institutes of Health (NIH). Additional support was provided by NIH U10 HL064287, U10 HL064288, U10 HL064295, U10 HL064305, U10 HL064307, U01 HL064313, RC2 HL101487, R01 HL087699 and U01 HL65899, an NIH Pharmacogenomics Research Network (PGRN) – RIKEN Center for Genomic Medicine (CGM) Global Alliance and by funding from the BioBank Japan project that was supported by the Ministry of Education, Culture, Sports, Sciences and Technology of the Japanese government. We thank the American Lung Association Asthma Clinical Research Centers (ALA-ACRC) for use of LOCCS and LODO study samples. B.E.H. was supported by NIH K99 HL105663. L.K.W. is funded by the American Asthma Foundation, the Fund for Henry Ford Hospital, and by NIH R01AI079139 and R01AI061774 from the National Institute of Allergy and Infectious Diseases (NIAID). K.C.B. was supported in part by the Mary Beryl Patch Turnbull Scholar Program. S.J.L. is supported by the Division of Intramural Research, National Institute of Environmental Health Sciences. D.A.G. was supported by the Netherlands Asthma Foundation grant AF (AF 95.09, AF 98.48, AF and AF and a grant from the University Medical Center Groningen. GHK was supported by a Ter Meulen Fund grant from the Royal Netherlands Academy of Arts and Sciences. D.S.P. was supported by a Royal Netherlands Academy of Arts and Sciences Honorary Professorship. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.


1. Akinbami L, Moorman J, Bailey C, Zahran H, King M, et al.. (2012) Trends in asthma prevalence, health care use, and mortality in the United States, 2001–2010. Hyattsville, MD: National Center for Health Statistics.
2. Global Strategy for the Diagnosis, Management and Prevention of COPD, Global Initiative for Chronic Obstructive Lung Disease (GOLD) 2007. Available: [PubMed]
3. Abiola O, Angel JM, Avner P, Bachmanov AA, Belknap JK, et al. (2003) The nature and identification of quantitative trait loci: a community's view. Nat Rev Genet 4: 911–916 [PMC free article] [PubMed]
4. Peters LL, Robledo RF, Bult CJ, Churchill GA, Paigen BJ, et al. (2007) The mouse as a model for human biology: a resource guide for complex trait analysis. Nat Rev Genet 8: 58–69 [PubMed]
5. Zosky GR, Larcombe AN, White OJ, Burchell JT, Janosi TZ, et al. (2008) Ovalbumin-sensitized mice are good models for airway hyperresponsiveness but not acute physiological responses to allergen inhalation. Clin Exp Allergy 38: 829–838 [PubMed]
6. Wills-Karp M, Ewart SL (1997) The genetics of allergen-induced airway hyperresponsiveness in mice. Am J Respir Crit Care Med 156: S89–96 [PubMed]
7. Henderson WR Jr, Tang LO, Chu SJ, Tsao SM, Chiang GK, et al. (2002) A role for cysteinyl leukotrienes in airway remodeling in a mouse asthma model. Am J Respir Crit Care Med 165: 108–116 [PubMed]
8. Ewart SL, Kuperman D, Schadt E, Tankersley C, Grupe A, et al. (2000) Quantitative trait loci controlling allergen-induced airway hyperresponsiveness in inbred mice. Am J Respir Cell Mol Biol 23: 537–545 [PubMed]
9. De Sanctis GT, Singer JB, Jiao A, Yandava CN, Lee YH, et al. (1999) Quantitative trait locus mapping of airway responsiveness to chromosomes 6 and 7 in inbred mice. Am J Physiol 277: L1118–1123 [PubMed]
10. Ober C, Tsalenko A, Parry R, Cox NJ (2000) A second-generation genomewide screen for asthma-susceptibility alleles in a founder population. Am J Hum Genet 67: 1154–1162 [PubMed]
11. Xu X, Fang Z, Wang B, Chen C, Guang W, et al. (2001) A genomewide search for quantitative-trait loci underlying asthma. Am J Hum Genet 69: 1271–1277 [PubMed]
12. De Sanctis GT, Merchant M, Beier DR, Dredge RD, Grobholz JK, et al. (1995) Quantitative locus analysis of airway hyperresponsiveness in A/J and C57BL/6J mice. Nat Genet 11: 150–154 [PubMed]
13. Zhang Y, Lefort J, Kearsey V, Lapa e SilvaJR, Cookson WO, et al. (1999) A genome-wide screen for asthma-associated quantitative trait loci in a mouse model of allergic asthma. Hum Mol Genet 8: 601–605 [PubMed]
14. Yokouchi Y, Nukaga Y, Shibasaki M, Noguchi E, Kimura K, et al. (2000) Significant evidence for linkage of mite-sensitive childhood asthma to chromosome 5q31–q33 near the interleukin 12 B locus by a genome-wide search in Japanese families. Genomics 66: 152–160 [PubMed]
15. Xu J, Meyers DA, Ober C, Blumenthal MN, Mellen B, et al. (2001) Genomewide screen and identification of gene-gene interactions for asthma-susceptibility loci in three U.S. populations: collaborative study on the genetics of asthma. Am J Hum Genet 68: 1437–1446 [PubMed]
16. Daniels SE, Bhattacharrya S, James A, Leaves NI, Young A, et al. (1996) A genome-wide search for quantitative trait loci underlying asthma. Nature 383: 247–250 [PubMed]
17. Wjst M, Fischer G, Immervoll T, Jung M, Saar K, et al. (1999) A genome-wide search for linkage to asthma. German Asthma Genetics Group. Genomics 58: 1–8 [PubMed]
18. Pletcher MT, McClurg P, Batalov S, Su AI, Barnes SW, et al. (2004) Use of a dense single nucleotide polymorphism map for in silico mapping in the mouse. PLoS Biol 2: e393. [PMC free article] [PubMed]
19. Kang HM, Zaitlen NA, Wade CM, Kirby A, Heckerman D, et al. (2008) Efficient control of population structure in model organism association mapping. Genetics 178: 1709–1723 [PubMed]
20. Yang H, Bell TA, Churchill GA, Pardo-Manuel de Villena F (2007) On the subspecific origin of the laboratory mouse. Nat Genet 39: 1100–1107 [PubMed]
21. Frazer KA, Eskin E, Kang HM, Bogue MA, Hinds DA, et al. (2007) A sequence-based variation map of 8.27 million SNPs in inbred mouse strains. Nature 448: 1050–1053 [PubMed]
22. Yang H, Ding Y, Hutchins LN, Szatkiewicz J, Bell TA, et al. (2009) A customized and versatile high-density genotyping array for the mouse. Nat Methods 6: 663–666 [PMC free article] [PubMed]
23. Leme AS, Berndt A, Williams LK, Tsaih SW, Szatkiewicz JP, et al. (2010) A survey of airway responsiveness in 36 inbred mouse strains facilitates gene mapping studies and identification of quantitative trait loci. Mol Genet Genomics 283: 317–326 [PMC free article] [PubMed]
24. Ober C, Hoffjan S (2006) Asthma genetics 2006: the long and winding road to gene discovery. Genes Immun 7: 95–100 [PubMed]
25. Moffatt MF, Gut IG, Demenais F, Strachan DP, Bouzigon E, et al. (2010) A large-scale, consortium-based genomewide association study of asthma. N Engl J Med 363: 1211–1221 [PubMed]
26. Torgerson DG, Ampleford EJ, Chiu GY, Gauderman WJ, Gignoux CR, et al. (2011) Meta-analysis of genome-wide association studies of asthma in ethnically diverse North American populations. Nat Genet 43: 887–892 [PMC free article] [PubMed]
27. Manolio TA (2010) Genomewide association studies and assessment of the risk of disease. N Engl J Med 363: 166–176 [PubMed]
28. Myers RA, Himes BE, Gignoux CR, Yang JY, Gauderman WJ, et al.. (2012) Further Replication Studies of the EVE Consortium Meta-Analysis Identifies Two Asthma Risk Loci in European Americans. J Allergy Clin Immunol In press. [PMC free article] [PubMed]
29. Himes BE, Jiang X, Hu R, Wu AC, Lasky-Su JA, et al. (2012) Genome-Wide Association Analysis in Asthma Subjects Identifies SPATS2L as a Novel Bronchodilator Response Gene. PLoS Genet 8: e1002824. [PMC free article] [PubMed]
30. Tantisira KG, Lasky-Su J, Harada M, Murphy A, Litonjua AA, et al. (2011) Genomewide association between GLCCI1 and response to glucocorticoid therapy in asthma. N Engl J Med 365: 1173–1183 [PMC free article] [PubMed]
31. Hong MG, Reynolds CA, Feldman AL, Kallin M, Lambert JC, et al. (2011) Genome-wide and gene-based association implicates FRMD6 in Alzheimer disease. Hum Mutat 33: 521–529 [PMC free article] [PubMed]
32. Eleftherohorinou H, Hoggart CJ, Wright VJ, Levin M, Coin LJ (2011) Pathway-driven gene stability selection of two rheumatoid arthritis GWAS identifies and validates new susceptibility genes in receptor mediated signalling pathways. Hum Mol Genet 20: 3494–3506 [PubMed]
33. Kim SH, Cho BY, Park CS, Shin ES, Cho EY, et al. (2009) Alpha-T-catenin (CTNNA3) gene was identified as a risk variant for toluene diisocyanate-induced asthma by genome-wide association analysis. Clin Exp Allergy 39: 203–212 [PubMed]
34. Thierry-Mieg D, Thierry-Mieg J (2006) AceView: a comprehensive cDNA-supported gene and transcripts annotation. Genome Biol 7 Suppl 1: S12 11–14. [PMC free article] [PubMed]
35. Morohashi Y, Hatano N, Ohya S, Takikawa R, Watabiki T, et al. (2002) Molecular cloning and characterization of CALP/KChIP4, a novel EF-hand protein interacting with presenilin 2 and voltage-gated potassium channel subunit Kv4. J Biol Chem 277: 14965–14975 [PubMed]
36. Ten Broeke R, Blalock JE, Nijkamp FP, Folkerts G (2004) Calcium sensors as new therapeutic targets for asthma and chronic obstructive pulmonary disease. Clin Exp Allergy 34: 170–176 [PubMed]
37. Ten Broeke R, Brandhorst MC, Leusink-Muis T, Villain M, De Clerck F, et al. (2003) Ca2+ sensors modulate asthmatic symptoms in an allergic model for asthma. Eur J Pharmacol 476: 151–157 [PubMed]
38. Churchill GA, Airey DC, Allayee H, Angel JM, Attie AD, et al. (2004) The Collaborative Cross, a community resource for the genetic analysis of complex traits. Nat Genet 36: 1133–1137 [PubMed]
39. Yang H, Wang JR, Didion JP, Buus RJ, Bell TA, et al. (2011) Subspecific origin and haplotype diversity in the laboratory mouse. Nat Genet 43: 648–655 [PMC free article] [PubMed]
40. Silverman ES, Palmer LJ, Subramaniam V, Hallock A, Mathew S, et al. (2004) Transforming growth factor-beta1 promoter polymorphism C-509T is associated with asthma. Am J Respir Crit Care Med 169: 214–219 [PubMed]
41. Silverman ES, Baron RM, Palmer LJ, Le L, Hallock A, et al. (2002) Constitutive and cytokine-induced expression of the ETS transcription factor ESE-3 in the lung. Am J Respir Cell Mol Biol 27: 697–704 [PubMed]
42. Peters SP, Anthonisen N, Castro M, Holbrook JT, Irvin CG, et al. (2007) Randomized comparison of strategies for reducing treatment in mild persistent asthma. N Engl J Med 356: 2027–2039 [PubMed]
43. American Lung Association Asthma Clinical Research Centers (2007) Clinical trial of low-dose theophylline and montelukast in patients with poorly controlled asthma. Am J Respir Crit Care Med 175: 235–242 [PubMed]
44. Luca D, Ringquist S, Klei L, Lee AB, Gieger C, et al. (2008) On the use of general control samples for genome-wide association studies: genetic matching highlights causal variants. Am J Hum Genet 82: 453–463 [PubMed]
45. Willer CJ, Sanna S, Jackson AU, Scuteri A, Bonnycastle LL, et al. (2008) Newly identified loci that influence lipid concentrations and risk of coronary artery disease. Nat Genet 40: 161–169 [PubMed]
46. International HapMap Consortium (2005) A haplotype map of the human genome. Nature 437: 1299–1320 [PMC free article] [PubMed]
47. 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 [PubMed]
48. Koppelman GH, Meyers DA, Howard TD, Zheng SL, Hawkins GA, et al. (2009) Identification of PCDH1 as a novel susceptibility gene for bronchial hyperresponsiveness. Am J Respir Crit Care Med 180: 929–935 [PMC free article] [PubMed]
49. Vonk JM, Postma DS, Maarsingh H, Bruinenberg M, Koppelman GH, et al. (2010) Arginase 1 and arginase 2 variations associate with asthma, asthma severity and beta2 agonist and steroid response. Pharmacogenet Genomics 20: 179–186 [PubMed]
50. Childhood Asthma Management Program Research Group (1999) The Childhood Asthma Management Program (CAMP): design, rationale, and methods. Control Clin Trials 20: 91–120 [PubMed]
51. Guilbert TW, Morgan WJ, Krawiec M, Lemanske RF Jr, Sorkness C, et al. (2004) The Prevention of Early Asthma in Kids study: design, rationale and methods for the Childhood Asthma Research and Education network. Control Clin Trials 25: 286–310 [PubMed]
52. Szefler SJ, Chinchilli VM, Israel E, Denlinger LC, Lemanske RF Jr, et al. (2012) Key observations from the NHLBI Asthma Clinical Research Network. Thorax 67: 450–455 [PMC free article] [PubMed]
53. Lasky-Su J, Himes BE, Raby BA, Klanderman BJ, Sylvia JS, et al. (2012) HLA-DQ strikes again: Genome-wide association study further confirms HLA-DQ in the diagnosis of asthma among adults. Clin Exp Allergy 42: 1724–1733 [PubMed]
54. Wheeler DL, Barrett T, Benson DA, Bryant SH, Canese K, et al. (2005) Database resources of the National Center for Biotechnology Information. Nucleic Acids Res 33: D39–45 [PMC free article] [PubMed]
55. Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, et al. (2010) LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics 26: 2336–2337 [PMC free article] [PubMed]
56. Fisher RA (1948) Combining independent tests of significance. Am Statistician 2: 30
57. R Development Core Team (2008) R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing.

Articles from PLoS ONE are provided here courtesy of Public Library of Science