Pedigree linkage disequilibrium mapping of quantitative trait loci

In this paper, we propose to use pedigrees of any size and any types of relatives in joint high-resolution linkage disequilibrium (LD) and linkage mapping of quantitative trait loci (QTL) by variance component models. Two or multiple markers can be simultaneously used in modeling association with the trait locus, instead of using one marker a time in the analysis. The proposed method can provide a unified result by using two or multiple markers in the modeling. This may avoid the complications of different results obtained from the separate analysis of marker by marker. The models simultaneously incorporate both linkage and LD information. The measures of LD are modeled by mean coefficients, and linkage information is modeled by variance-covariance matrix. Using analytical formulas to calculate the regression coefficients, the genetic effects are shown to be decomposed into additive and dominance components. The noncentrality parameter approximations of test statistics of LD are provided to make power calculations. Power and type I error rates are explored to investigate the merit of the proposed method by both the analytical formulas and simulations. Comparing with the association between-family and association within-family ('AbAw') approach of Fulker and Abecasis et al, it is evident that the method proposed in this article is more powerful. The method is applied to investigate the relation between polymorphisms in the angiotensin 1-converting enzyme (ACE) genes and circulating ACE levels, with a better result than that of the 'AbAw' approach. Moreover, two markers I/D and 4656(CT)3/2 can fully interpret association with the trait locus at a 0.01 significance level, which provides a unique result for the ACE data.