 Proceedings
 Open Access
 Published:
HasemanElston weighted by marker informativity
BMC Genetics volume 6, Article number: S50 (2005)
Abstract
In the HasemanElston approach the squared phenotypic difference is regressed on the proportion of alleles shared identical by descent (IBD) to map a quantitative trait to a genetic marker. In applications the IBD distribution is estimated and usually cannot be determined uniquely owing to incomplete marker information. At Genetic Analysis Workshop (GAW) 13, Jacobs et al. [BMC Genet 2003, 4(Suppl 1):S82] proposed to improve the power of the HasemanElston algorithm by weighting for information available from marker genotypes. The authors did not show, however, the validity of the employed asymptotic distribution. In this paper, we use the simulated data provided for GAW 14 and show that weighting HasemanElston by marker information results in increased type I error rates. Specifically, we demonstrate that the number of significant findings throughout the chromosome is significantly increased with weighting schemes. Furthermore, we show that the classical HasemanElston method keeps its nominal significance level when applied to the same data. We therefore recommend to use HasemanElston with marker informativity weights only in conjunction with empirical pvalues. Whether this approach in fact yields an increase in power needs to be investigated further.
Background
The Haseman and Elston (HE) method [1] is the best known approach to map quantitative traits by linkage analyses. It has been extended in many different ways to increase statistical power. For example, Amos et al. [2] introduced a generalized least squares approach and weighted the squared phenotypic differences inversely proportional to their respective variances. Sham and Purcell [3] also proposed to weight phenotypes: they combined the squared trait difference and the squared mean centered sum linearly, the weights given to these two components depending on the overall trait correlation.
A different approach has been followed by Jacobs et al. [4] at the Genetic Analysis Workshop (GAW) 13. They proposed to weight families for HE linkage analyses according to marker informativity (as measured by the difference between the allele sharing at the marker and the allele sharing for a noninformative marker) and applied their method to genome scan data for blood pressure. The anticipated gain in power is substantial. For example, the pvalue dropped from approximately 10^{2} to 10^{4} on chromosome 5 and from 10^{3} to 10^{5} on chromosome 12 for the mean systolic blood pressure. However, Jacobs et al. [4] did not show the validity of their approach, i.e., the asymptotic normality of the proposed test statistic.
Here, we pick up their work. We use the same distance metric as Jacobs et al. as well as a different simplexbased weighting scheme and show that the HE regression with weights according to this measure of marker informativity suffers from inflated type I error rates. We illustrate the effect of the weighted HE regression using data from chromosome 4 of the Aipotu population as available for GAW 14. For the analyses, we required and utilized the correct answers for the simulated data.
Simplex weighting scheme
In the simplex weighting scheme, families are weighted according to their degree of informativity for linkage. The degree of informativity is determined by the Euclidian distance d between the current available identitybydescent (IBD) marker information and the IBD information for uninformative markers.
The calculation of the simplex weights can be illustrated by an equilateral triangle of height 1 (Fig. 1). Viviani's theorem states that, in an equilateral triangle, for any point f = (f_{0}, f_{1}, f_{2}), the sum of the perpendiculars f_{i} from f to the sides of the triangle equals the triangle's height, i.e., f_{0} + f_{1} + f_{2} = 1. Therefore, any point in Figure 1 represents exactly one possible IBD distribution. To compute the Euclidian distance between IBD points f = (f_{0}, f_{1}, f_{2}) and g = (g_{0}, g_{1}, g_{2}), the IBD distributions f and g have to be mapped to the Cartesian coordinates (f_{x}, f_{y}) and (g_{x}, g_{y}), respectively. The required mapping function can easily be deduced from Figure 1 and is given by
If a genetic marker is noninformative, the IBD distribution equals (1/4,1/2,1/4). We therefore define simplex weights w as the Euclidian distance d between a marker with IBD distribution f = (f_{0}, f_{1}, f_{2}) and a noninformative marker by
Jacobs et al. [4] proposed slightly different weights and determined marker informativity by a distance metric D defined as
Other measures of informativity might be preferable but are beyond the scope of this paper.
Classical HE and weighted HE regression
For simplicity, we consider a sample of n independent sib pairs. Then, the classical HE regression ignores a possible dominance effect and fits the linear model
Δ^{2}= X_{ γ }+ ε, where
and the identity matrix Ω. The squared phenotypic difference of the i^{th} sib pair is represented by , and π_{ i }denotes the proportion of alleles IBD (π_{ i }= f_{i 2}+ f_{i 1}/ 2) of the i^{th} sib pair. In practice, a ttest statistic is employed to test the null hypothesis of no linkage. Under H_{0} the test statistic
is asymptotically distributed as t_{n1}, where n  1 denotes the degrees of freedom. The parameter estimators for β and are given by
As soon as weights are introduced, Ω no longer represents the identity matrix. Specifically, we aim to use simplex weights as described above, so that Ω becomes a diagonal matrix with elements w_{ i }(f). Since the IBD distribution is estimated from the available marker data, replaces the true matrix Ω. Test statistics and estimates are obtained from equations (3) and (4) where Ω is replaced by its estimate .
Methods
We illustrate the effect of the weighted HE regression by using data from chromosome 4 of Aipotu. The ethnic group as well as the chromosome was randomly chosen. No positive finding in the microsatellite scan can be expected for chromosome 4.
We employed a development branch of SIBPAL and GENIBD from S.A.G.E. [5] to compute three different asymptotic HE regression models for each of the 100 replicates: the classical HE method as well as HE weighted by weights according to Equations (1) and (2), respectively. GENIBD was utilized to estimate multipoint IBD distributions; estimates were obtained at the marker positions.
We wished to investigate the validity of the asymptotic distribution of the weighted HE methods. This can be achieved by testing whether the number of significant findings across a chromosome is significantly increased. For sake of simplicity, we used the 5% test level for all further investigations. For each genetic marker position, each of the 100 replicates and each weighting scheme, p < 0.05 was tested. Generalized estimating equations with an autoregressive (AR(1)) working correlation structure, using the replicate as class level indicator, was used for each weighting scheme to investigate whether the proportion of significant findings across the chromosome exceeded 5%. This model adequately adjusts for intermarker correlations on a chromosome. The estimated proportion of significant findings across the whole chromosome is reported with its corresponding asymptotic 95% confidence interval (CI).
Furthermore, we want to show that the number of significant findings is greater using the weighted HE methods compared with the classical HE approach. To this end, a Wilcoxon signed rank test was employed. Specifically, we counted the number of microsatellite positions where p < 0.05 for both weighted and the classical HE regression across the whole chromosome. The Wilcoxon signed rank tested was computed across the 100 independent chromosomes. If the weighted HE methods was too liberal, the number of significant findings would be significantly increased compared with the classical HE method.
Results
Table 1 shows that the number of significant results at a marker position across the 100 replicates is at least as high for the weighted HE regressions as for classical HE method. This finding was substantiated by the generalized estimating equation model, which shows that the proportion of significant findings exceeds the nominal test level of 5% for both weighted HE approaches ( = 0.0685; 95% CI: 0.0546–0.0823 for the simplex weighting scheme; = 0.0634; 95% CI: 0.0562–0.0706 for the weighting scheme proposed by Jacobs et al. [4]) but not for the classical HE method ( = 0.0552; 95% CI: 0.0431–0.0673). Furthermore, the number of markers with positive linkage is greater for both weighted HE regressions compared with the classical HE method (p = 3.3 × 10^{7} for the simplex weight (Eq. (1)), p = 7.1 × 10^{6} for the weights proposed by Jacobs et al. [4]). These findings are invariant to increased sample size as achieved by pairwise pooling of replicates (detailed results not shown).
Conclusion
Weighting HE regression models by informativity is an appealing approach. However, some care is needed when applying this approach to real data. If phenotypes are weighted appropriately within families, this may result in substantial gain in statistical power [2, 3]. Instead of weighting phenotypes, an increase in power might also be obtained by weighting according to marker informativity. This approach has been successfully utilized in the context of metaanalyses in which studies have been weighted according to their informativity [6, 7]. The method of Jacobs et al. [4] combines both approaches. Jacobs et al. weighted individual sibpair families in the HE regression according to their marker informativity as indicated above. However, the asymptotic normality of the proposed test statistic was not shown.
In this paper, using simulated data available for GAW14, we demonstrated that the HE method with familywise weights according to marker informativity suffers from inflated type I error levels when we measure informativity as the distance between current IBD marker information and the IBD information for uninformative markers. If pvalues are computed asymptotically from weighted regression models only, the ttest statistic maybe distorted.
We therefore recommend the use of weighting functions in conjunction only with empirically computed pvalues until a theoretical solution to the detected problem is found. Furthermore, it needs to be investigated whether empirical p values in fact yield an increase in statistical power for weighted HE compared to the classical HE regression.
Abbreviations
 CI:

Confidence interval
 GAW:

Genetic Analysis Workshop
 HE:

Haseman and Elston
 IBD:

Identical by descent
References
 1.
Haseman JK, Elston RC: The investigation of linkage between a quantitative trait and a marker locus. Behav Genet. 1972, 2: 319. 10.1007/BF01066731.
 2.
Amos CJ, Elston RC, Wilson AF, BaileyWilson JE: A more powerful robust sibpair test of linkage for quantitative traits. Genetic Epidemiol. 1989, 6: 435449. 10.1002/gepi.1370060306.
 3.
Sham PC, Purcell S: Equivalence between HasemanElston and variancecomponents linkage analyses for sib pairs. Am J Hum Genet. 2001, 68: 15271532. 10.1086/320593.
 4.
Jacobs KB, GrayMcGuire C, Cartier KC, Elston RC: Genomewide linkage scan for genes affecting longitudinal trends in systolic blood pressure. BMC Genet. 2003, 4 (Suppl 1): S8210.1186/147121564S1S82.
 5.
S.A.G.E.: Statistical Analysis for Genetic Epidemiology. 2004, 4.6: [http://darwin.cwru.edu/sage/]
 6.
Loesgen S, Dempfle A, Gölla A, Bickeböller H: Weighting schemes in pooled linkage analysis. Genet Epidemiol. 2001, 21 (Suppl 1): S142S147.
 7.
Dempfle A, Loesgen S: Metaanalysis of linkage studies for complex diseases: an overview of methods and a simulation study. Ann Hum Genet. 2004, 68 (Pt 1): 6983. 10.1046/j.15298817.2003.00061.x.
Acknowledgements
This work was supported by Deutsche Forschungsgemeinschaft grant ZI 591/121.
Author information
Additional information
Authors' contributions
AZ had the original idea for the study. DF and AK did the statistical analyses. RCE was an intellectual consultant on the study. AZ and DF wrote the report. All authors read and approved the final manuscript.
Rights and permissions
About this article
Cite this article
Franke, D., Kleensang, A., Elston, R.C. et al. HasemanElston weighted by marker informativity. BMC Genet 6, S50 (2005) doi:10.1186/147121566S1S50
Published
DOI
Keywords
 Weighting Scheme
 Marker Informativity
 Marker Information
 Genetic Analysis Workshop
 Inflated Type