Volume 6 Supplement 1
Genetic Analysis Workshop 14: Microsatellite and singlenucleotide polymorphism
Haplotypesharing analysis using Mantel statistics for combined genetic effects
 Lars Beckmann^{1}Email author,
 Christine Fischer^{2},
 Markus Obreiter^{1},
 Michael Rabes^{1} and
 Jenny ChangClaude^{1}
DOI: 10.1186/147121566S1S70
© Beckmann et al; licensee BioMed Central Ltd 2005
Published: 30 December 2005
Abstract
We applied a new approach based on Mantel statistics to analyze the Genetic Analysis Workshop 14 simulated data with prior knowledge of the answers. The method was developed in order to improve the power of a haplotype sharing analysis for gene mapping in complex disease. The new statistic correlates genetic similarity and phenotypic similarity across pairs of haplotypes from casecontrol studies. The genetic similarity is measured as the shared length between haplotype pairs around a genetic marker. The phenotypic similarity is measured as the mean corrected crossproduct based on the respective phenotypes. Cases with phenotype P1 and unrelated controls were drawn from the population of Danacaa. Power to detect main effects was compared to the X^{2}test for association based on 3marker haplotypes and a global permutation test for haplotype association to test for main effects. Power to detect gene × gene interaction was compared to unconditional logistic regression. The results suggest that the Mantel statistics might be more powerful than alternative tests.
Background
Recently we proposed a flexible approach to gene mapping of complex diseases, whereby we combine Mantel statistics for spacetime clustering with genetic information obtained from haplotypes [1]. It has been shown that haplotype sharing methods are well suited for mapping such genes [2–5]. Mantel statistics were introduced in 1967 to correlate temporal and spatial distributions of cancer, notably childhood leukemia, in a generalized regression approach [6]. The Mantel statistic M is the sum of the cross product of the spatial similarity X_{ ij }multiplied by the temporal similarity Y_{ ij }across all pairs of cases i and j:
The idea behind this approach is that in the presence of spacetime clustering the values of spatial similarity X_{ ij }correspond to the values of temporal similarity Y_{ ij }for correlated cases i and j.
Methods
Mantel statistics using haplotypes
Here we apply the general approach of Mantel's statistics for spacetime clustering (Equation 1) to correlate genetic and phenotypic similarity, and to test for gene × gene interaction. The first statistic has the form:
where x denotes a genetic marker, and i and j are haplotypes. L_{ ij }(x) denotes the genetic similarity between the haplotypes i and j at x, and is defined as the number of intervals surrounding x that are flanked by markers with the same alleles, i.e., that are identical by state (IBS). The phenotypic similarity for two haplotype copies i and j derived from individuals s_{ i }and s_{ j }is defined as the mean corrected product Y_{ sisj }= (y_{ si } μ)(y_{ si } μ), where y_{ si }and y_{ sj }are the phenotypes of s_{ i }and s_{ j }, and μ denotes the expectation of the phenotype. Here, we chose μ as the sample mean, i.e., μ = 0.5. Concordant pairs of affected and concordant pairs of unaffected individuals have the weights Y_{ sisj }= 0.25, while discordant pairs have the weights Y_{ sisj }= 0.25. Alternative measures of phenotypic similarity were discussed in the framework of sibpair analysis, e.g., the HasemanElston method [7] and the weighted pairwise correlation statistics [8], as well as in familybased association analysis [9]. The summation is over all pairwise comparisons of haplotypes for i ≠ j, where the haplotypes are derived from casecontrol studies.
The second statistic is constructed to test for the combined effect of two loci:
The information of the first locus x is incorporated as the shared length L_{ ij }(x). At the second locus only genotype information is used. The variable z_{ si }is coded in a dominant way, i.e., z_{ si }is 1, if the individual s_{ i }carries at least one mutant allele, and 0 otherwise. The measure of genotypic similarity Z_{ sisj }is then 1, if z_{ si }= z_{ sj }, and 0 otherwise.
The summands of the Mantel statistic are highly correlated, and any statistical procedure to test for significance has to take into account the interrelationship of the data. Here, we use a Monte Carlo permutation approach to test for significance, as proposed by Mantel [6]. For M_{0}(x) the phenotype y_{ si }is permuted over the individuals. The definition of Z is such that M_{1}(x) is the sum over all comparisons of haplotypes from individuals who have the same genotype coding z at the second locus. To derive the null hypothesis of no statistical interaction, the phenotype y_{ si }and the genotype coding z_{ si }at the second locus for individual s_{ i }are permuted jointly over the individuals, and thus the comparisons of haplotypes derived from discordant individuals are incorporated under the null hypothesis.
Statistical tests for comparison
Main effects
We used two alternative tests for power comparison.
1. We applied the X^{2}test for association to 3marker haplotypes. The region of interest was covered by overlapping sliding windows. The haplotypes consisted of 3 consecutive genetic markers. The test was based on a 2xk X^{2}table, with k 1 degrees of freedom, where k denotes the number of haplotypes that occurred in either the case or the control sample. A pvalue was assigned to the marker in the center of the window. Note that no tests were performed for the marginal markers.
2. The haplotype assignment software PHASE [10, 11] performs a global permutation test for significant differences in haplotype frequencies in case and control groups. PHASE tests the null hypothesis that the case and control haplotypes are a random sample from a single set of haplotype frequencies, versus the alternative that cases are more similar to other cases than to controls. Here, this test was based on 100 permutations due to computational burden.
Gene × gene interaction
We compared the test statistic M_{1}(x) using haplotypes to unconditional logistic regression based on the genotypes at 2 genetic markers [12]. The respective genotypes were coded for both the recessive and the dominant model.
Datasets and genetic data
The casecontrol study samples for two different samples sizes were drawn from the population Danacaa to limit the analysis to individuals defined by phenotype P1.
In this dataset, two major genes, D1 and D2, interacted in an epistatic model. Mode of inheritance is dominant for both D1 and D2.
Study samples used in the analysis
Sample  Number of cases/controls  Replicates used  SNPs/microsatellites  Number of markers  

Cases  Controls  
A  200/200  1, 2  3, 4  C01R0045 – 0055  11 
B  200/200  1, 2  3, 4  C01R0045 – 0055, D01S0021 – 0024  15 
C  400/400  1–4, 18  5–8  C01R0045 – 0046, C01R0050 – 0053, C01R0055  7 
D  400/400  1–4, 18  5–8  B01T0555 – 0559, C01R0052, B01T0561 – 0565  11 
E  200/200  1, 2  3, 4  C03R0276 – 0281  6 
F  400/400  1–4, 18  5–8  C03R0276 – 0281  6 
Software
Haplotype pairs assigned to the unrelated individuals were estimated by the use of the PHASE program [10, 11]. PHASE lists the most likely pairs of haplotypes for each individual, together with their posterior probability. The most likely (best) estimate of haplotype pairs was chosen for our analysis. SAS 8.02 (SAS Institute Inc., Cary, NC, USA) was used to test for normality and for logistic regression. All other calculations were performed with software developed within our group. Software for the proposed Mantel statistics is available upon request.
Results
Main effects
Results of the Mantel statistic (x) and the haplotypebased X^{2}_{hap} – test for main effects
Sample  M_{0}(x)  X ^{2} _{hap}  Permutation test global  

Nearest marker (pvalue)  Marker with lowest pvalue (pvalue)  Nearest marker  Marker with lowest pvalue (pvalue)  
Chromosome 1  
A  C01R0052 (0.778)  C01R0054 (0.273)  C01R0052 (0.648)  C01R0047 (0.159)  0.41 
B  C01R0052 (0.042)  C01R0052 (0.042)  (ND^{a})  (ND^{a})  0.09 
C  C01R0052 (0.579)  C01R0046 (0.054)  C01R0052 (0.471)  C01R0045 (0.31)  0.81 
D  C01R0052 (0.068)  C01R0045 (0.014)  C01R0052 (0.095)  C01R0052 (0.095)  0.03 
Chromosome 3  
E  C03R0281 (0.134)  C03R0281 (0.134)  C03R0280^{b} (0.668)  C03R0279 (0.629)  0.68 
F  C03R0281 (0.043)  C03R0280 (0.002)  C03R0280^{b} (0.11)  C03R0279 (0.09)  0.15 
M_{0}(x) did not yield significant results for the markers flanking D2 with small sample size. The most significant SNP in the large sample was C03R0280 (p = 0.002). The X^{2}_{hap}test for association, however, did not produce significant results with either the small or the large samples. The permutation test yielded one globally significant pvalue of 0.03 in the large sample D.
Gene × gene interaction
Results of the Mantel statistic M_{1}(x) to test for gene × gene interaction
Sample, nearest marker  pValue  Marker with lowest pvalue  pValue 

Chromosome 1, C01R0052  0.062  
A  0.122  C01R0048  0.007 
B  0.064  C01R0048  0.001 
C  0.009  C01R0053  0.009 
D  0.031  C01R0045  
Chromosome 3, C03R0281  
E  0.02  C03R0281  0.02 
F  0.003  C03R0281  0.003 
Conclusion
We successfully employed a new approach to map disease predisposing genes in casecontrol studies based on Mantel statistics that correlate genetic and phenotypic similarity. Two types of gene effects involved in complex diseases were considered: main effects and joint effects.
1. The Mantel statistic M_{0}(x) identified the major gene D2 on chromosome 3 given adequate sample size, whereas the alternative methods failed. Major gene D1 on chromosome 1 was simulated without linkage disequilibrium (LD). LD is necessary for haplotype association methods, therefore M_{0}(x)as expecteddid not map D1 correctly.
We acknowledge that the comparison against the X^{2} association test for 3 marker haplotypes is somewhat unfair, but we know of no other standard association test examining longer haplotypes that is not confronted with problems of huge degrees of freedom and sparse data. Additionally, other more sophisticated haplotypebased methods cannot yet be regarded as standard.
2. The Mantel statistic M_{1}(x) accounted for the joint effects of 2 putative disease loci. Taking the combined effects into account, the results were significant for the major genes D1 and D2 and showed lower pvalues than the results obtained when considering main effects only.
These results show that main effects might not be detectable if gene × gene interaction is present and not considered in the analysis. Our proposed method M_{1}(x) revealed significant statistical interaction between the genes analyzed in contrast to the results obtained in the logistic regression model.
The proposed Mantel statistics employ haplotypes from casecontrol data and might not be robust to population stratification. In our analysis, we used samples drawn from the Danacaa population and affection status defined by phenotype P1 to reduce heterogeneity in the data. Population stratification is therefore not a major concern in this analysis. We did not adjust the pvalues for multiple comparisons in this candidate analysis.
Multiple testing is a serious problem especially if all possible gene × gene interactions increase the multiplicity. We solved the problem in the mean time by implementing a stepdown algorithm to take into account multiple testing [13, 14].
Comprehensive power comparisons are currently being carried out to reveal under which conditions our approach is more powerful than alternative methods.
Abbreviations
 GAW14:

Genetic Analysis Workshop 14
 IBS:

Identical by state
 SNP:

Singlenucleotide polymorphism
Declarations
Acknowledgements
This work was supported by a Deutsche Forschungsgemeinschaft grant (CH117/31) (LB, MR, MO). We thank Kati Smit for technical assistance.
Authors’ Affiliations
References
 Beckmann L, Thomas D, Fischer C, ChangClaude J: Haplotype sharing analysis using Mantel statistics. Hum Hered. 2005, 59: 6778. 10.1159/000085221.View ArticlePubMedGoogle Scholar
 Beckmann L, Fischer C, Deck KG, Nolte IM, te Meerman G, ChangClaude J: Exploring haplotype sharing methods in general and isolated populations to detect gene(s) of a complex genetic trait. Genet Epidemiol. 2001, 21 (Suppl 1): S554S559.PubMedGoogle Scholar
 Fischer C, Beckmann L, Majoram P, te Meerman G, ChangClaude J: Haplotype sharing analysis with SNPs in candidate genes: the Genetic Analysis Workshop 12 example. Genet Epidemiol. 2003, 24: 6873. 10.1002/gepi.10207.View ArticlePubMedGoogle Scholar
 Qian D, Thomas DC: Genome scan of complex traits by haplotype sharing correlation. Genet Epidemiol. 2001, 21 (Suppl 1): S582S587.PubMedGoogle Scholar
 Boon M, Nolte IM, Bruinenberg M, Spijker GT, Terpstra P, Raelson J, De Keyser J, Zwanikken CP, Hulsbeek M, Hofstra RM, Buys CH, te Meerman GJ: Mapping of a susceptibility gene for multiple sclerosis to the 51 kb interval between G511525 and D6S1666 using a new method of haplotype sharing analysis. Neurogenetics. 2001, 3: 221230.PubMedGoogle Scholar
 Mantel N: The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27: 209220.PubMedGoogle Scholar
 Wang K, Huang J: A scorestatistic approach for the mapping of quantitativetrait loci with sibships of arbitrary size. Am J Hum Genet. 2002, 70: 412424. 10.1086/338659.PubMed CentralView ArticlePubMedGoogle Scholar
 Commenges D, BeurtonAimar M: Multipoint linkage analysis using the weightedpairwise correlation statistic. Genet Epidemiol. 1999, 17 (Suppl 1): S515S519.View ArticlePubMedGoogle Scholar
 Lunetta KL, Faraone SV, Biederman J, Laird NM: Familybased tests of association and linkage that use unaffected sibs, covariates, and interactions. Am J Hum Genet. 2000, 66: 605614. 10.1086/302782.PubMed CentralView ArticlePubMedGoogle Scholar
 Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001, 68: 978989. 10.1086/319501.PubMed CentralView ArticlePubMedGoogle Scholar
 Stephens M, Donnelly P: A comparison of bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet. 2003, 73: 11621169. 10.1086/379378.PubMed CentralView ArticlePubMedGoogle Scholar
 Gauderman WJ: Sample size requirements for association studies of genegene interaction. Am J Epidemiol. 2002, 155: 478484. 10.1093/aje/155.5.478.View ArticlePubMedGoogle Scholar
 Ge YC, Dudoit S, Speed TP: Resamplingbased multiple testing for microarray data analysis. Test. 2003, 12: 177. 10.1007/BF02595811.View ArticleGoogle Scholar
 Obreiter M, Fischer C, ChangClaude J, Beckman L: SDMinP: a program to control the family wise error rate using stepdown minP adjusted Pvalues. Bioinformatics. 2005, 21: 31833184. 10.1093/bioinformatics/bti480.View ArticlePubMedGoogle Scholar
Copyright
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.