Skip to main content

A graphical assessment of p-values from sliding window haplotype tests of association to identify asthma susceptibility loci on chromosome 11q

Abstract

Background

Past work on asthmatic African American families revealed a strong linkage peak with modest evidence of association on chromosome 11q. Here, we perform tests of association for asthma and a panel of 609 SNPs in African American subjects using a sliding window approach. While efficient in screening a region of dense genotyping, this approach does create some problems: high numbers of tests, assimilating thousands of results, and questions about setting priorities on regions with association signals.

Results

We present a newly developed tool, Graphical Assessment of Sliding P-values or GrASP, which uses color display to indicate the width of the sliding windows, significance of individual tests, density of SNP coverage and location of known genes that simplifies some of these issues, and use it to identify regions of interest in these data.

Conclusion

We demonstrate that GrASP makes it easier to visualize, summarize and prioritize regions of interest from sliding window haplotype analysis, based jointly on the p-value from all the tests from these windows and the building of haplotypes of significance in the region. Using this approach, five regions yielded strong evidence for linkage and association with asthma, including the prior peak linkage region.

Background

To search for causal genes for asthma on chromosome 11q, a region previously highlighted in linkage and association studies on African American families by our group, a dense panel of SNPs was genotyped using an Illumina™ fine mapping panel in these African American subjects. The sliding window approach that has gained much favor in recent years [1–6] to test for association between SNPs and haplotypes in this region was adapted to analyze case-parent trios and case-control data because it is a simple and efficient way to screen a region of dense genotyping. Given an ordered set of markers (1, 2, 3,....n), sliding windows of overlapping haplotypes are tested in sequence, i.e. markers 1-2-3 are treated as a single haplotype, then markers 2-3-4 are treated as a single haplotype, then markers 3-4-5, etc. Haplotypes of varying sizes (i.e. 2-, 3-, 4-SNP haplotypes, etc.) can be assessed. When large numbers of SNPs are used, this approach generates large numbers of non-independent statistical tests, and visualizing and interpreting these multiple tests can become daunting. Currently, it is common practice to use some form of summary measure (such as the minimum p-value of all sliding windows or windows of a particular size that include any one SNP) as a synopsis of each SNP. Here, we present a novel tool called Gr aphical A ssessment of S liding P-values or GrASP that provides a graphical overview of all tests from sliding windows without sub selection.

Asthma, a disease of chronic airway inflammation, shows strong familial aggregation [7–9] but models of inheritance are inconsistent, and linkage studies over the past decade have suggested that causal genes may be located on over two dozen chromosomal regions [10, 11]. Heterogeneity in the genetic control of asthma is obvious from the variable evidence for linkage and/or association across studies and even among families within a single study. Recently, positional cloning studies have identified four novel candidate genes for asthma, including the genes encoding PDH finger protein 11 (PHF11) on chromosome 13q14 [12] dipeptidyl peptidase IV-related protein 3 (DPP10) on chromosome 2q14-q32 [13], a disintegrin and metalloproteinase domain 33 (ADAM33) on chromosome 20p13 [14], and G protein-coupled receptor for asthma susceptibility (GPRA) on 7p [15].

In a genome-wide scan and subsequent fine-mapping analysis, we have previously presented evidence of linkage of asthma to chromosome 11q markers in an African American population, but not Caucasian or Hispanic families [16, 17]. Further, family based association and transmission disequilibrium tests (TDT) showed significant evidence of linkage and linkage disequilibrium for several individual markers in this region. A putative susceptibility locus was estimated to be at map position 70.8 cM with a confidence interval spanning the linkage peak [17]. Several candidate genes in this region were examined for association with asthma: CRTH2, which encodes a receptor for prostaglandin D2 (PGD2) [18]; FCER1B, encoding the high affinity IgE receptor β subunit [19]; and CC10, the Clara cell 10-kD immunomodulator [20]. Family-based and case-control analysis of asthma and two common SNPs (G1544C and G1651A) in the 3' untranslated region of CRTH2 showed significant evidence of linkage in the presence of disequilibrium in this African American and an independent Chinese population [21]. Here, we present further fine-mapping of this region on 11q in the African American families and apply GrASP to these data to illustrate its ability to summarize results from thousands of tests while identify regions of interest based on all single-SNP and haplotype tests of association.

Results

Samples were genotyped at the SNP Center of the Genetic Resources Core Facility (GRCF) using the Illumina™ Fine Mapping Panel #38 spanning 52 Mb on chromosome 11q. The GRCF released genotypes on a total of 705 loci; the overall error rate for the entire data set was 0.55 % (considering duplicated samples – 46 events/8,410 genotypes) and the overall parent-child discordance rate was 0.07 % (4 events/5,612 genotypes). An additional 36 loci had minor allele frequencies < 5%, 8 had 50% Gencall scores < 0.40 [22], 18 markers had reproducibility between the duplicate samples (6 comparison pairs) of < 100%, 7 had Mendelian inconsistencies in families > 2%, and 43 showed significant departure from Hardy Weinberg Equilibrium (p-value < 0.01).

The final data for the family-based analysis included 86 probands, 63 fathers and 75 mothers from the 89 nuclear families resulting in 60 child-father pairs, 73 child-mother pairs or 51 complete trios. The case-control analysis for replication included the 86 probands and 94 control samples. A total of 609 markers on 318 samples were used in the analyses, with an average distance between SNPs of 86119 bp (range: 127 – 3676182 bp). Average linkage disequilibrium (LD) between neighboring SNPs was modest with an average D' = 0.53 and average r2 = 0.11.

Single SNP results

A total of 48 SNPs had p-values < 0.05 in the case-control data and 34 in the family data (6 of which were significant in both data sets). Table 1 presents results for these 6 common SNPs, and for 20 SNPs that had p-values < 0.01 in either the case-control or family data, along with the known genes where they are located.

Table 1 Summary of single SNP results from case-control- and family-based tests.

Selecting regions for haplotype analysis

Considering 609 SNPs and sliding window sizes from 2 – 6 SNPs per window, there are 3030 non-independent tests for each study design (608 2-SNP windows, 607 3-SNP windows, 606 4-SNP windows, 605 5-SNP windows and 604 6-SNP windows) for a total of 6060 tests. Hence, our first approach was to construct systematic rules based on prior information about LD and/or single SNP results to create blocks or segments of contiguous SNPs within which sliding window-based tests could be conducted. The density of SNPs in this fine mapping panel was not optimal for establishing patterns of high LD or constructing large haplotype blocks; although there were 16 blocks of LD identified using the approach of Gabriel et al. [23], most were very small (11 had only 2 SNPs per block, 3 had 3 SNPs, and one each of 5 and 6 SNPs per block). The physical distance across the blocks was also modest (14 were under 25 Kb, one at 136 Kb and only one at 486 Kb). To restrict haplotype tests with sliding windows within these haplotype blocks would have ignored much of the data available here.

Therefore, we began with those single SNPs with marginal levels of significance and expanded outward from these to define segments of contiguous SNPs over which to implement sliding window haplotype analysis. There were 76 SNPs that had a positive signal (p-value < 0.05) in either the case-control or family-based analysis. As a first step, we created a mini-segment including all SNPs within a 1 Mb region centered around each of these 76 SNPs (i.e. 500 Kb on either side). In the next step, these 76 mini-segments were merged wherever two segments overlapped. This resulted in 15 contiguous main segments (covering 498 SNPs), with 15 gaps (covering 111 SNPs) between segments (where no evidence against the null hypothesis was seen). We then performed sliding window haplotype tests within each segment but not across segments, which resulted in 2268 tests for each study design.

Sliding window haplotype results

Figures 1 and 2 are a compact graphical overview of all tests performed on the family data. These figures represent p-values from single SNP TDT tests and overall p-values from haplotypes constructed across sliding windows of sizes 2–6. Each test is represented by a single tile: a box with a black outline which is as long as the number of SNPs in the window, and the p-value of the test is represented by the color of the box. Several important points become immediately apparent and are discussed below.

Figure 1
figure 1

Graphical Assessment of sliding p-values for all tests performed on SNPs 1–303 in the family-based data with prioritized regions of association signal (labelled A-C).

Figure 2
figure 2

Graphical Assessment of sliding p-values for all tests performed on SNPs 304–609 in the family-based data with prioritized regions of association signal (labelled D-E).

Examining p-values for any one SNP (i.e. considering the minimum p-value across all windows that included a SNP), reveal many SNPs had significant results (79 with minimum p-value <0.05, 25 <0.01 and 5 <0.001). Prioritizing regions of interest for further follow up is not a simple task, especially since the minimum p-value should be weighed against the number of tests yielding significant results when setting priorities. GrASP provides a visual presentation of the results that can simplify the identification of regions with both high significance and where multiple tests yielded high significance. Here, we consider regions of highest priority as those that have a building of windows with significance across multiple SNPs.

Relying upon these patterns, five regions depicted in Figures 1 and 2 (and labelled A – E) stand out as "regions of interest". Each of these regions contain SNPs with minimum p-values that are statistically significant, multiple windows in the region showing significance, and have a pattern of building from smaller to larger windows until the strength of apparent association drops off for the largest haplotypes (as would be expected considering the breakdown of the original putative causal haplotype over physical distance).

Figures 3 and 4 present these five regions (labelled A-E) in more detail, including evidence from the case-control analysis, the physical location of SNPs and known genes in this region (both features available in GrASP – either relying on user-provided input or by querying public databases of Entrez). Examining the details of chromosomal regions using this approach can simplify the process of identifying regions of association and prime candidate genes. The five regions with strong patterns in the family data also had signals of association in the case-control data, although the two did not always align perfectly.

Figure 3
figure 3

Detailed overview of the first three regions with strongest association signals (regions A-C in Figure 1) for case-control- and family-based tests along with SNP and known gene locations. Color coding of p-values are as presented in Figure 1.

Figure 4
figure 4

Detailed overview of the final two regions with strongest association signals (regions D-E in Figure 2) for case-control- and family-based tests along with SNP and known gene locations. Color coding of p-values are as presented in Figure 1.

Discussion

In this report we present a graphical tool for haplotype association analysis utilizing data of dense SNPs on chromosome 11q for asthma in an African American population, and identify five regions that may include susceptibility loci for asthma.

High throughput genotyping and dense SNP maps makes the sliding window approach for analyzing haplotypes appealing, and several groups have explored this approach from a statistical [24–29] and applied perspective [1–6, 30]. In this sliding window approach, windows of varying size are examined, generally beginning with the first SNP and moving windows down the framework map one SNP at a time. This strategy is simple and efficient in comprehensively screening a dense region of genotyping for association with the trait of interest. This approach, however, does raise problematic issues: exceedingly high numbers of non-independent tests, the task of assimilating thousands of test results, attempting to create haplotypes across SNPs far from each other, and prioritizing "regions of interest" based on sometimes weak signals from these multiple tests. Here, we provide a novel graphical simplification to some of these issues in our implementation of the Gr aphical A ssessment of S liding P-values (GrASP) analysis tool.

The use of graphical presentations and has proved extremely useful in other analytical strategies for mass amounts of data, including linkage analysis [31], linkage disequilibrium patterns [32], identification of haplotype block structure [23], and most commonly in high-throughput gene expression profiling (e.g. microarray) [33–38]. In many of these diverse applications, color representation can quickly reveal underlying patterns. The use of sliding haplotype windows yields high numbers of tests; in this study sliding windows of 2–6 SNPs over 609 SNPs yielded a maximum of 3030 tests. Limiting regions of sliding windows can reduce the number of tests, but summarizing 2268 tests to identify patterns of association is still difficult. Although using a summary measure, such as the minimum p-value for each SNP, can simplify this issue, it would be ideal to have all test results available for assessment without selecting sub-sets. Here too, as in the case of LD and gene expression patterns, a graphical overview (as implemented in GrASP) is both practical and effective.

The simple assimilation of all the results itself (seen in Figures 1 and 2) overcomes the task of presenting 2268 test results simultaneously, and clearly shows most of the regions have no evidence of statistically significant signals (i.e. green background). However, there are several regions with distinct signals (yellow, orange and red cells), but only a handful of these are compelling (defined by clusters of yellow, orange and red cells).

Setting priorities for regions of association from case-control or case-parent trio data as used here is a key step in identifying putative asthma and atopy genes in this region on chromosome 11q. For any single SNP, up to 21 tests of association are performed (one single SNP test, two 2-SNP windows that include the SNP, three 3-SNP windows, four 4-SNP windows, five 5-SNP windows, and six 6-SNP windows that include the SNP) using this sliding window approach. Relying solely upon minimum p-value for each SNP from these tests clearly ignores the number of tests performed. Graphically examining the strength of the signal for each SNP (the minimum p-value) and the number of tests providing evidence at this SNP as part of increasingly larger haplotypes using GrASP presents both types of information jointly (Figures 1 and 2). Currently one can weigh the evidence in priority regions informally, but the application of formal rules (including aspects such as effect size and haplotype frequencies) could be implemented. Based on this approach, 5 regions can be identified as priorities for subsequent mapping, but our prioritization scheme is not meant to rule out the few additional regions with smaller signals (see Figures 1 and 2). How to narrow the putative region of interest is another major question, and we show that this issue can also be simplified with GrASP. Here, by searching for increasingly significant haplotypes in the areas of interest (as depicted in Figures 3 and 4), the location of the strongest signals given the current genotyping density is easy to determine. This provides the investigator a reasonable starting point for additional genotyping, particularly for SNPs to be typed within and flanking a current block of association.

Finally, as we were unable to use haplotype blocks to guide our haplotype analyses, we substantially reduced the number of sliding window tests by relying on our single SNP results. In these data, all strong signals from the sliding windows were centered at or around a SNP with an individually significant association. However, considering situations with very few single-SNP signals or extremely dense SNP data with high LD, our rules for forming segments (within which sliding window tests are done) may discard too much or too little data. Furthermore, limiting haplotype tests in regions highlighted by single SNPs tests may introduce several issues: missing regions of association that would only have been detected with a haplotype test (i.e. those for which the single-SNP tests did not have sufficient power), and perpetuation of false positives especially when limiting haplotype tests only to SNPs in tight LD with those single SNPs yielding a significant association. Considering issues such as these, it is more appropriate that the sliding window approach should consider the density of SNPs at hand first analyzing smaller windows and then proceeding to larger windows where there are signals from the smaller windows rather than restricting windows to regions based on single-SNP results. Given our SNP density (here relatively sparse for fine mapping), it appears that the windowing approach should be implemented only for windows of size 2–4. The 5- and 6-SNP windows should be tested only where the smaller windows reveal positive associations. This in itself would have reduced the number of tests by an additional 864 (38%) for each analysis type. GrASP now makes this strategy to reduce the number of tests for sliding window approach easy to implement, and thereby efficient.

The five regions of priority contain many known genes as illustrated in Figures 3 and 4. Although most of these genes are not immediately recognized as candidates in the pathophysiology of asthma and atopy, a few are notable. Two genes (ELF5 and EHS) from region A encode members of the ETS transcription factor subfamily [39], show epithelial-specific expression, and their gene products may act as a transcriptional repressors. These two genes have been suggested as candidate genes for asthma in a genome-wide linkage analysis in Tristan da Cunha [40, 41], and variants of the EHS gene were shown to be weakly associated with asthma susceptibility in a Caucasian case-control study [42, 43]. The gene LRP5 in region D encodes the low density lipoprotein receptor-related protein 5 and has been linked to the susceptibility to osteoprosis and IDDM [44]. It is important to note that the peak fine-mapping marker, D11S1337, is situated within this locus [17]. Interestingly, about 400 kb away from LRP5, the TCIRG1 (gene which encodes a vacuolar proton pump H(+)-ATPase) is also included in a region of significant association. Two different gene products through alternative splicing have been found for TCIRG1, one of which, TIRC7, encodes a protein essential in T-cell activation, and is almost exclusively expressed in the immune system [45].

Conclusion

In conclusion, we have demonstrated that the use of GrASP can serve as a step towards simplifying multiple issues in the sliding window approach for haplotype analysis: summarizing test results, setting priorities for defining regions of interest and narrowing these regions of interest. Here, GrASP is used on fine-mapping SNP data, but it may also assist greatly in assessing SNPs within candidate genes of interest and even more so in genome-screen SNP data analysis. We have identified five regions with strong evidence for association with asthma in these African American subjects which contain many known genes including a few noteworthy candidates. Most compelling, is the strong evidence from region four that also includes our peak linkage marker reported previously.

Methods

Study population

Study protocols for the Collaborative Study on the Genetics of Asthma (CSGA) were approved by the Institutional Review Boards of Johns Hopkins University School of Medicine and Howard University College of Medicine. Written informed consent was obtained from each of the study subjects. A total of 89 African American case-parent trios were included[17, 21]. The ascertainment scheme and clinical characteristics of these families have been described in detail elsewhere [16]. Briefly, families were ascertained through an asthmatic sib-pair and asthmatic probands in case-parent trios. In the case-control design, 89 unrelated individuals were selected from among asthmatic probands of 49 multiplex African American families and 40 case-parent trios. In addition, 95 randomly selected healthy individuals were included as a replication study in a case-control study design. We should note that while our sample size is limited, these data are a subset of the original families that provided compelling evidence for linkage and association in this region of 11q, and therefore a good sample for fine-mapping.

Genotyping

Genetic analysis was conducted on either genomic DNAs or DNA amplified directly from blood using multiple displacement amplification (MDA) technology (Molecular Staging Inc.) [46, 47]. Genotyping was performed at the SNP Center of the Genetic Resources Core Facility (GRCF) at the McKusick/Nathans Institute of Genetic Medicine, Johns Hopkins School of Medicine. Genotypes were generated on a BeadLab 1000 system [22] using the Illumina™ Fine Mapping Panel #38 (Illumina OPA ID = FM78) spanning 52 Mb on chromosome 11q (35–87 Mb – using the Build 34 map).

Statistical methods

At the single-SNP level, the family-based association test (FBAT) as described by Rabinowitz and Laird [48] with the empirical variance-covariance estimator was used for the family-based sample and in the case-control design, χ2 tests were used to test for independence between single SNPs and disease status.

For the haplotype analysis, sliding windows of 2 – 6 SNPs each were considered as the approach allows for comprehensive screening of the whole region of genotyping. In the family-based analysis, haplotype tests were performed using methods described by Horvath et al [49] which assumes no recombination between markers, and is analogous to the single marker multi-allelic statistic described by Rabinowitz and Laird [48] allowing for weighting all phased genotypes (haplotypes) possible in the family. An EM algorithm which maximizes the likelihood of the phased haplotype frequencies based on all observed family genotypes and computed under the null hypothesis was used to estimate haplotype frequencies used to obtain these weights. For the case-control data, group-wise haplotype frequencies were estimated using an EM algorithm [50]. χ2 values were calculated for each haplotype versus all others, and empiric significance was assessed using 10,000 permutations of the observed case-control data. The "omnibus test" was evaluated to detect differences in overall haplotype frequency profiles between asthmatics and controls [50]. For both the family- and case-control analyses only the overall test statistic was considered, i.e. one single overall test per sliding window, and not the individual tests of deviation from expected for each haplotype.

Graphical overview of sliding windows

Gr aphical A ssessment of S liding P-values or GrASP, a tool to graphically present and assess p-values from sliding window haplotype tests, was implemented. This program can summarize thousands of p-values from sliding window tests of variable sizes in a simple graphic that uses varying levels of a user-specified color to indicate the width of the sliding windows and the level of statistical significance attained. With this tool it now becomes a simple task to identify regions/blocks of interest from these sliding windows, based jointly on the absolute p-value of the tests and the building of haplotypes of significance in the region. GrASP is executed as an Excel macro, and is written in Excel's built-in version of Visual Basic for Applications. It uses as input, a summary file that contains the test results stored in simple format (SNP order, first SNP of each window and p-values from windows of varying sizes) that is easily constructed from analysis output considering the repetitive nature of the output results. GrASP also illustrates SNP location and the location of known genes in the region as a physical track that is drawn either using user-provided input, or by querying the public databases of Entrez. GrASP is freely available for use at: http://research.nhgri.nih.gov/GrASP/.

References

  1. Hellmig S, Mascheretti S, Renz J, Frenzel H, Jelschen F, Rehbein JK, Folsch U, Hampe J, Schreiber S: Haplotype analysis of the CD11 gene cluster in patients with chronic Helicobacter pylori infection and gastric ulcer disease. Tissue Antigens. 2005, 65: 271-274. 10.1111/j.1399-0039.2005.00362.x.

    Article  PubMed  CAS  Google Scholar 

  2. Gibson F, Froguel P: Genetics of the APM1 locus and its contribution to type 2 diabetes susceptibility in French Caucasians. Diabetes. 2004, 53: 2977-2983.

    Article  PubMed  CAS  Google Scholar 

  3. Skipper L, Wilkes K, Toft M, Baker M, Lincoln S, Hulihan M, Ross OA, Hutton M, Aasly J, Farrer M: Linkage disequilibrium and association of MAPT H1 in Parkinson disease. Am J Hum Genet. 2004, 75: 669-677. 10.1086/424492.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  4. Nakajima T, Wooding S, Sakagami T, Emi M, Tokunaga K, Tamiya G, Ishigami T, Umemura S, Munkhbat B, Jin F, Guan-Jun J, Hayasaka I, Ishida T, Saitou N, Pavelka K, Lalouel JM, Jorde LB, Inoue I: Natural selection and population history in the human angiotensinogen gene (AGT): 736 complete AGT sequences in chromosomes from around the world. Am J Hum Genet. 2004, 74: 898-916. 10.1086/420793.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  5. Innan H, Padhukasahasram B, Nordborg M: The pattern of polymorphism on human chromosome 21. Genome Res. 2003, 13: 1158-1168. 10.1101/gr.466303.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  6. Paschou P, Feng Y, Pakstis AJ, Speed WC, DeMille MM, Kidd JR, Jaghori B, Kurlan R, Pauls DL, Sandor P, Barr CL, Kidd KK: Indications of linkage and association of Gilles de la Tourette syndrome in two independent family samples: 17q25 is a putative susceptibility region. Am J Hum Genet. 2004, 75: 545-560. 10.1086/424389.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  7. Ã…berg N: Familial occurrence of atopic disease: genetic versus environmental factors. Clin Exp Allergy. 1993, 23: 829-34. 10.1111/j.1365-2222.1993.tb00260.x.

    Article  PubMed  Google Scholar 

  8. Cooke RA, Jr AVV: Human Sensitization. J Immunol. 1916, 1: 210-305.

    Google Scholar 

  9. Gerrard JW, Ko CG, Vickers P, Gerrard CD: The familial incidence of allergic disease. Annals of Allergy. 1976, 36: 10-15.

    PubMed  CAS  Google Scholar 

  10. Barnes KC: Atopy and asthma genes--where do we stand?. Allergy. 2000, 55 (9): 803-17. 10.1034/j.1398-9995.2000.00123.x.

    Article  PubMed  CAS  Google Scholar 

  11. Ober C, Moffatt MF: Contributing factors to the pathobiology. The genetics of asthma. Clin Chest Med. 2000, 21: 245-61. 10.1016/S0272-5231(05)70264-1.

    Article  PubMed  CAS  Google Scholar 

  12. Zhang Y, Leaves NI, Anderson GG, Ponting CP, Broxholme J, Holt R, Edser P, Bhattacharyya S, Dunham A, Adcock IM, Pulleyn L, Barnes PJ, Harper JI, Abecasis G, Cardon L, White M, Burton J, Matthews L, Mott R, Ross M, Cox R, Moffatt MF, Cookson WO: Positional cloning of a quantitative trait locus on chromosome 13q14 that influences immunoglobulin E levels and asthma. Nat Genet. 2003, 34: 181-186. 10.1038/ng1153.

    Article  PubMed  CAS  Google Scholar 

  13. Allen M, Heinzmann A, Noguchi E, Abecasis G, Broxholme J, Ponting CP, Bhattacharyya S, Tinsley J, Zhang Y, Holt R, Jones EY, Lench N, Carey A, Jones H, Dickens NJ, Dimon C, Nicholls R, Baker C, Xue L, Townsend E, Kabesch M, Weiland SK, Carr D, Von Mutius E, Adcock IM, Barnes PJ, Lathrop GM, Edwards M, Moffatt MF, Cookson WO: Positional cloning of a novel gene influencing asthma from chromosome 2q14. Nat Genet. 2003, 35: 258-263. 10.1038/ng1256.

    Article  PubMed  CAS  Google Scholar 

  14. Van Eerdewegh P, Little RD, Dupuis J, Del Mastro RG, Falls K, Simon J, Torrey D, Pandit S, McKenny J, Braunschweiger K, Walsh A, Liu Z, Hayward B, Folz C, Manning SP, Bawa A, Saracino L, Thackston M, Benchekroun Y, Capparell N, Wang M, Adair R, Feng Y, Dubois J, FitzGerald MG, Huang H, Gibson R, Allen KM, Pedan A, Danzig MR, Umland SP, Egan RW, Cuss FM, Rorke S, Clough JB, Holloway JW, Holgate ST, Keith TP: Association of the ADAM33 gene with asthma and bronchial hyperresponsiveness. Nature. 2002, 418: 426-430. 10.1038/nature00878.

    Article  PubMed  CAS  Google Scholar 

  15. Laitinen T, Polvi A, Rydman P, Vendelin J, Pulkkinen V, Salmikangas P, Makela S, Rehn M, Pirskanen A, Rautanen A, Zucchelli M, Gullsten H, Leino M, Alenius H, Petays T, Haahtela T, Laitinen A, Laprise C, Hudson TJ, Laitinen LA, Kere J: Characterization of a common susceptibility locus for asthma-related traits. Science. 2004, 304: 300-304. 10.1126/science.1090010.

    Article  PubMed  CAS  Google Scholar 

  16. Xu J, Meyers DA, Ober C, Blumenthal MN, Mellen B, Barnes KC, King RA, Lester LA, Howard TD, Solway J, Langefeld CD, Beaty TH, Rich SS, Bleecker ER, Cox NJ: 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. 2001, 68: 1437-1446. 10.1086/320589.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  17. Huang SK, Mathias RA, Ehrlich E, Plunkett B, Liu X, Cutting GR, Wang XJ, Li XD, Togias A, Barnes KC, Malveaux F, Rich S, Mellen B, Lange E, Beaty TH: Evidence for asthma susceptibility genes on chromosome 11 in an African-American population. Hum Genet. 2003, 113: 71-75. 10.1007/s00439-003-0965-x.

    Article  PubMed  CAS  Google Scholar 

  18. Hirai H, Tanaka K, Yoshie O, Ogawa K, Kenmotsu K, Takamori Y, Ichimasa M, Sugamura K, Nakamura M, Takano S, Nagata K: Prostaglandin D2 selectively induces chemotaxis in T helper type 2 cells, eosinophils, and basophils via seven-transmembrane receptor CRTH2. J Exp Med. 2001, 193: 255-261. 10.1084/jem.193.2.255.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  19. Hizawa N, Yamaguchi E, Jinushi E, Konno S, Kawakami Y, Nishimura M: Increased total serum IgE levels in patients with asthma and promoter polymorphisms at CTLA4 and FCER1B. J Allergy Clin Immunol. 2001, 108: 74-9. 10.1067/mai.2001.116119.

    Article  PubMed  CAS  Google Scholar 

  20. Chen LC, Zhang Z, Myers AC, Huang SK: Cutting edge: altered pulmonary eosinophilic inflammation in mice deficient for Clara cell secretory 10-kDa protein. J Immunol. 2001, 167: 3025-3028.

    Article  PubMed  CAS  Google Scholar 

  21. Huang JL, Gao PS, Mathias RA, Yao TC, Chen LC, Kuo ML, Hsu SC, Plunkett B, Togias A, Barnes KC, Stellato C, Beaty TH, Huang SK: Sequence variants of the gene encoding chemoattractant receptor expressed on Th2 cells (CRTH2) are associated with asthma and differentially influence mRNA stability. Hum Mol Genet. 2004, 13: 2691-2697. 10.1093/hmg/ddh279.

    Article  PubMed  CAS  Google Scholar 

  22. Fan JB, Oliphant A, Shen R, Kermani BG, Garcia F, Gunderson KL, Hansen M, Steemers F, Butler SL, Deloukas P, Galver L, Hunt S, McBride C, Bibikova M, Rubano T, Chen J, Wickham E, Doucet D, Chang W, Campbell D, Zhang B, Kruglyak S, Bentley D, Haas J, Rigault P, Zhou L, Stuelpnagel J, Chee MS: Highly parallel SNP genotyping. Cold Spring Harb Symp Quant Biol. 2003, 68: 69-78. 10.1101/sqb.2003.68.69.

    Article  PubMed  CAS  Google Scholar 

  23. Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, Liu-Cordero SN, Rotimi C, Adeyemo A, Cooper R, Ward R, Lander ES, Daly MJ, Altshuler D: The structure of haplotype blocks in the human genome. Science. 2002, 296: 2225-2229. 10.1126/science.1069424.

    Article  PubMed  CAS  Google Scholar 

  24. Zhao H, Pfeiffer R, Gail MH: Haplotype analysis in population genetics and association studies. Pharmacogenomics. 2003, 4: 171-178. 10.1517/phgs.4.2.171.22636.

    Article  PubMed  Google Scholar 

  25. Beckman L, Fischer C, Obreiter M, Rabes M, Chang-Claude J: Haplotype sharing analysis using Mantel statistics for combined genetic effects. BMC Genet. 2005, 6: S70-10.1186/1471-2156-6-S1-S70.

    Article  Google Scholar 

  26. Durrant C, Morris A: Linkage disequilibrium mapping via cladistic analysis of phase-unknown genotypes and inferred haplotypes in the Genetic Analysis Workshop 14 simulated data. BMC Genet. 2005, 6: S100-10.1186/1471-2156-6-S1-S100.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Meng Z, Zaykin DV, Xu CF, Wagner M, Ehm MG: Selection of genetic markers for association analyses, using linkage disequilibrium and haplotypes. Am J Hum Genet. 2003, 73: 115-130. 10.1086/376561.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  28. Lin S, Chakravarti A, Cutler DJ: Exhaustive allelic transmission disequilibrium tests as a new approach to genome-wide association studies. Nat Genet. 2004, 36: 1181-1188. 10.1038/ng1457.

    Article  PubMed  CAS  Google Scholar 

  29. Cheng R, Ma JZ, Wright FA, Lin S, Gao X, Wang D, Elston RC, Li MD: Nonparametric disequilibrium mapping of functional sites using haplotypes of multiple tightly linked single-nucleotide polymorphism markers. Genetics. 2003, 164: 1175-1187.

    PubMed  CAS  PubMed Central  Google Scholar 

  30. Chiu YF, Liu X, Tsai YY: A comparison in association and linkage genome-wide scans for alcoholism susceptibility genes using single-nucleatide polymorphisms. BMC Genet. 2005, 6: S89-10.1186/1471-2156-6-S1-S89.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Pugh EW, Mandal DM, Wilson AF: A graphical approach for presenting linkage results from a genomic screen. Genet Epidemiol. 1995, 12: 807-812. 10.1002/gepi.1370120646.

    Article  PubMed  CAS  Google Scholar 

  32. Abecasis GR, Cookson WO: GOLD--graphical overview of linkage disequilibrium. Bioinformatics. 2000, 16: 182-3. 10.1093/bioinformatics/16.2.182.

    Article  PubMed  CAS  Google Scholar 

  33. Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  34. Yi Y, Mirosevich J, Shyr Y, Matusik R, George ALJ: Coupled analysis of gene expression and chromosomal location. Genomics. 2005, 85: 401-412. 10.1016/j.ygeno.2004.11.011.

    Article  PubMed  CAS  Google Scholar 

  35. Lingjaerde OC, Baumbusch LO, Liestol K, Glad IK, Borresen-Dale AL: CGH-Explorer: a program for analysis of array-CGH data. Bioinformatics. 2005, 21: 821-822. 10.1093/bioinformatics/bti113.

    Article  PubMed  CAS  Google Scholar 

  36. Awad IA, Rees CA, Hernandez-Boussard T, Ball CA, Sherlock G: Caryoscope: an Open Source Java application for viewing microarray data in a genomic context. BMC Bioinformatics. 2004, 5: 151-10.1186/1471-2105-5-151.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Chi B, DeLeeuw RJ, Coe BP, MacAulay C, Lam WL: SeeGH--a software tool for visualization of whole genome array comparative genomic hybridization data. BMC Bioinformatics. 2004, 5: 13-10.1186/1471-2105-5-13.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Dresen IM, Husing J, Kruse E, Boes T, Jockel KH: Software packages for quantitative microarray-based gene expression analysis. Curr Pharm Biotechnol. 2003, 4: 417-437. 10.2174/1389201033377436.

    Article  PubMed  CAS  Google Scholar 

  39. Sharrocks AD, Brown AL, Ling Y, Yates PR: The ETS-domain transcription factor family. Int J Biochem Cell Biol. 1997, 29: 1371-1387. 10.1016/S1357-2725(97)00086-1.

    Article  PubMed  CAS  Google Scholar 

  40. Slutsky AS, Zamel N: Genetics of asthma: the University of Toronto Program. University of Toronto Genetics of Asthma Research Group. Am J Respir Crit Care Med. 1997, 156: S130-S132.

    Article  PubMed  CAS  Google Scholar 

  41. Sreekumar GP, Binnie CG, Carden EG, Sharma S, Sprankle CS, Sherman PA, Xue ZY, Fling ME, Reisman M, Blumenthal M, Vestbo J, Sundy J, Pericak-Vance MA, Brewster S, Ehm MG, Burns DK, Wagner MJ: Association of asthma and related phenotypes with ASTH-I/J genes on chromosome 11p in Caucasians. Am J Hum Genet. 2000

    Google Scholar 

  42. Baron RM, Palmer LJ, Tantisira K, Gabriel S, Sonna LA, Le L, Hallock A, Libermann TA, Drazen JM, Weiss ST, Silverman ES: DNA sequence variants in epithelium-specific ETS-2 and ETS-3 are not associated with asthma. Am J Respir Crit Care Med. 2002, 166: 927-932. 10.1164/rccm.200201-048OC.

    Article  PubMed  Google Scholar 

  43. Silverman ES, Baron RM, Palmer LJ, Le L, Hallock A, Subramaniam V, Riese RJ, McKenna MD, Gu X, Libermann TA, Tugores A, Haley KJ, Shore S, Drazen JM, Weiss ST: Constitutive and cytokine-induced expression of the ETS transcription factor ESE-3 in the lung. Am J Respir Cell Mol Biol. 2002, 27: 697-704.

    Article  PubMed  CAS  Google Scholar 

  44. http://www.ncbi.nlm.nih.gov/entrez/dispomim.cgi?id=603506. 2006

  45. Utku N, Heinemann T, Tullius SG, Bulwin GC, Beinke S, Blumberg RS, Beato F, Randall J, Kojima R, Busconi L, Robertson ES, Schulein R, Volk HD, Milford EL, Gullans SR: Prevention of acute allograft rejection by antibody targeting of TIRC7, a novel T cell membrane protein. Immunity. 1998, 9: 509-518. 10.1016/S1074-7613(00)80634-2.

    Article  PubMed  CAS  Google Scholar 

  46. Barker DL, Hansen MS, Faruqi AF, Giannola D, Irsula OR, Lasken RS, Latterich M, Makarov V, Oliphant A, Pinter JH, Shen R, Sleptsova I, Ziehler W, Lai E: Two methods of whole-genome amplification enable accurate genotyping across a 2320-SNP linkage panel. Genome Res. 2004, 14: 901-907. 10.1101/gr.1949704.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  47. http://www.molecularstaging.com. 2006

  48. Rabinowitz D, Laird N: A unified approach to adjusting association tests for population admixture with arbitrary pedigree structure and arbitrary missing marker information. Hum Hered. 2000, 50: 211-23. 10.1159/000022918.

    Article  PubMed  CAS  Google Scholar 

  49. Horvath S, Xu X, Lake SL, Silverman EK, Weiss ST, Laird NM: Family-based tests for associating haplotypes with general phenotype data: application to asthma genetics. Genet Epidemiol. 2004, 26: 61-69. 10.1002/gepi.10295.

    Article  PubMed  Google Scholar 

  50. Fallin D, Schork NJ: Accuracy of haplotype frequency estimation for biallelic loci, via the expectation-maximization algorithm for unphased diploid genotype data. Am J Hum Genet. 2000, 67: 947-959. 10.1086/303069.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors would like to thank Dr. Stephen Anyaibe, Maria Stockton, Eva Ehrlich, Beverly Plunkett and Deborah Gause for their work in clinical recruitment and laboratory efforts as part of the Collaborative Study on the Genetics of Asthma, Drs. Alan Scott and David Valle for support and encouragement in the development of GrASP, and Drs. Kimberly Doheny, Ya-Yu Tsai and Ivy McMullen for their help with the SNP data. This research was supported in part by the Intramural Research Program of the National Human Genome Research Institute, National Institutes of Health. This work is supported, in part, by an NIH grant (RO1 AI-052468) and a research grant from Philip Morris USA Inc. KCB was supported in part by the Mary Beryl Patch Turnbull Scholar Program.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Rasika A Mathias.

Additional information

Authors' contributions

RAM conceived of the graphical method, participated in the acquisition of data and family recruitment, performed statistical analysis and wrote the manuscript. PG carried out part of the molecular genetic studies and helped draft the manuscript. JLG did all the programming for GrASP and helped draft the manuscript. AFW and EWP helped design GrASP and helped draft the manuscript. PFH, GMD, FJM and AT participated in the design of the study and recruitment of families for the study and helped draft the manuscript. KCB and THB participated in study design, recruitment of families for the study, guidance on statistical analysis and helped draft the manuscript. SKH was involved in study design, molecular genetic studies and helped draft the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Mathias, R.A., Gao, P., Goldstein, J.L. et al. A graphical assessment of p-values from sliding window haplotype tests of association to identify asthma susceptibility loci on chromosome 11q. BMC Genet 7, 38 (2006). https://doi.org/10.1186/1471-2156-7-38

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2156-7-38

Keywords