Recent advances in human quantitative-trait-locus mapping: comparison of methods for discordant sibling pairs

Extreme discordant sibling pairs (EDSPs) are theoretically powerful for the mapping of quantitative-trait loci (QTLs) in humans. EDSPs have not been used much in practice, however, because of the need to screen very large populations to find enough pairs that are extreme and discordant. Given appropriate statistical methods, another alternative is to use moderately discordant sibling pairs (MDSPs)--pairs that are discordant but not at the far extremes of the distribution. Such pairs can be powerful yet far easier to collect than extreme discordant pairs. Recent work on statistical methods for QTL mapping in humans has included a number of methods that, though not developed specifically for discordant pairs, may well be powerful for MDSPs and possibly even EDSPs. In the present article, we survey the new statistics and discuss their applicability to discordant pairs. We then use simulation to study the type I error and the power of various statistics for EDSPs and for MDSPs. We conclude that the best statistic(s) for discordant pairs (moderate or extreme) is (are) to be found among the new statistics. We suggest that the new statistics are appropriate for many other designs as well-and that, in fact, they open the way for the exploration of entirely novel designs.     

Efficiency robust tests for mapping quantitative trait loci using extremely discordant sib pairs

In 1972, Haseman and Elston proposed a pioneering regression method for mapping quantitative trait loci using randomly selected sib pairs. Recently, the statistical power of their method was shown to be increased when extremely discordant sib pairs are ascertained. While the precise genetic model may not be known, prior information that constrains IBD probabilities is often available. We investigate properties of tests that are robust against model uncertainty and show that the power gain from further constraining IBD probabilities is marginal. The additional linkage information contained in the trait values can be incorporated by combining the Haseman-Elston regression method and a robust allele sharing test.     

A powerful and robust new linkage statistic for discordant sibling pairs

Previously, Szatkiewicz and colleagues evaluated the performance of a wide variety of statistics for quantitative-trait-locus linkage, using discordant sibling pairs. They found that the most powerful statistics, in general, were a score statistic and a "composite statistic." However, whereas these two statistics have equal power under ideal conditions, each has limitations that reduce its power in certain circumstances. The score statistic depends on estimates of trait parameters and can lose a lot of power if those estimates are incorrect. The composite statistic is not sensitive to trait-parameter estimates but does depend on arbitrary weights that must be chosen on the basis of the ascertainment scheme. In this report, we elucidate the algebraic relationship between the score and composite statistics and then use that relationship to suggest a new statistic that combines the best properties of both. We call our new statistic the "robust discordant pair" (RDP) statistic. We report simulation studies to show that the RDP statistic does, indeed, have all of the strengths and none of the weaknesses of the score and composite statistics.     

Nonparametric simulation-based statistics for detecting linkage in general pedigrees.

We present here four nonparametric statistics for linkage analysis that test whether pairs of affected relatives share marker alleles more often than expected. These statistics are based on simulating the null distribution of a given statistic conditional on the unaffecteds' marker genotypes. Each statistic uses a different measure of marker sharing: the SimAPM statistic uses the simulation-based affected-pedigree-member measure based on identity-by-state (IBS) sharing. The SimKIN (kinship) measure is 1.0 for identity-by-descent (IBD) sharing, 0.0 for no IBD status sharing, and the kinship coefficient when the IBD status is ambiguous. The simulation-based IBD (SimIBD) statistic uses a recursive algorithm to determine the probability of two affecteds sharing a specific allele IBD. The SimISO statistic is identical to SimIBD, except that it also measures marker similarity between unaffected pairs. We evaluated our statistics on data simulated under different two-locus disease models, comparing our results to those obtained with several other nonparametric statistics. Use of IBD information produces dramatic increases in power over the SimAPM method, which uses only IBS information. The power of our best statistic in most cases meets or exceeds the power of the other nonparametric statistics. Furthermore, our statistics perform comparisons between all affected relative pairs within general pedigrees and are not restricted to sib pairs or nuclear families.     

Powerful allele sharing statistics for nonparametric linkage analysis

Nonparametric linkage analysis is widely used to map susceptibility genes for complex diseases. This paper introduces six nonparametric statistics for measuring marker allele sharing among the affected members of a pedigree. We compare the power of these new statistics and three previous statistics to detect linkage with Mendelian diseases having recessive, additive, and dominant modes of inheritance. The nine statistics represent all possible combinations of three different IBD scoring functions and three different schemes for sampling genes among affecteds. Our results strongly suggest that the statistic T(rec)(blocks) is best for recessive traits, while the two statistics T(kin)(pairs) and T(all)(kin) vie for best for an additive trait. The best statistic for a dominant trait is less clear. The statistics T(kin)(pairs) and T(all)(kin) are equally promising for small sibships, but in extended pedigrees the statistics T(dom)(blocks) and T(dom)(pairs) appear best. For a complex trait, we advocate computing several of these statistics.     

Powerful regression-based quantitative-trait linkage analysis of general pedigrees

We present a new method of quantitative-trait linkage analysis that combines the simplicity and robustness of regression-based methods and the generality and greater power of variance-components models. The new method is based on a regression of estimated identity-by-descent (IBD) sharing between relative pairs on the squared sums and squared differences of trait values of the relative pairs. The method is applicable to pedigrees of arbitrary structure and to pedigrees selected on the basis of trait value, provided that population parameters of the trait distribution can be correctly specified. Ambiguous IBD sharing (due to incomplete marker information) can be accommodated in the method by appropriate specification of the variance-covariance matrix of IBD sharing between relative pairs. We have implemented this regression-based method and have performed simulation studies to assess, under a range of conditions, estimation accuracy, type I error rate, and power. For normally distributed traits and in large samples, the method is found to give the correct type I error rate and an unbiased estimate of the proportion of trait variance accounted for by the additive effects of the locus-although, in cases where asymptotic theory is doubtful, significance levels should be checked by simulations. In large sibships, the new method is slightly more powerful than variance-components models. The proposed method provides a practical and powerful tool for the linkage analysis of quantitative traits.     

Nonpaternity in Linkage Studies of Extremely Discordant Sib Pairs

An approach commonly used to increase statistical power in linkage studies is the study of extremely discordant sibling pairs. This design is powerful under both additive and dominant-gene models and across a wide range of allele frequencies. A practical problem with the design is that extremely discordant pairs that are ostensibly full sibs may be half sibs. Although estimates vary, the population rates of such nonpaternity may be as high as 5%-10%. The proportion in discordant pairs may be much higher. The present article explores this potential inflation as a function of the resemblance of sib pairs and the criteria for discordance used for selection.     

On the power for linkage detection using a test based on scan statistics

We analyze some aspects of scan statistics, which have been proposed to help for the detection of weak signals in genetic linkage analysis. We derive approximate expressions for the power of a test based on moving averages of the identity by descent allele sharing proportions for pairs of relatives at several contiguous markers. We confirm these approximate formulae by simulation. The results show that when there is a single trait-locus on a chromosome, the test based on the scan statistic is slightly less powerful than that based on the customary allele sharing statistic. On the other hand, if two genes having a moderate effect on a trait lie close to each other on the same chromosome, scan statistics improve power to detect linkage.     

A linkage strategy for detection of human quantitative-trait loci. II. Optimization of study designs based on extreme sib pairs and generalized relative risk ratios.

We are concerned here with practical issues in the application of extreme sib-pair (ESP) methods to quantitative traits. Two important factors-namely, the way extreme trait values are defined and the proportions in which different types of ESPs are pooled, in the analysis-are shown to determine the power and the cost effectiveness of a study design. We found that, in general, combining reasonable numbers of both extremely discordant and extremely concordant sib pairs that were available in the sample is more powerful and more cost effective than pursuing only a single type of ESP. We also found that dividing trait values with a less extreme threshold at one end or at both ends of the trait distribution leads to more cost-effective designs. The notion of generalized relative risk ratios (the lambda methods, as described in the first part of this series of two articles) is used to calculate the power and sample size for various choices of polychotomization of trait values and for the combination of different types of ESPs. A balance then can be struck among these choices, to attain an optimum design.     

Recent advances in human quantitative-trait-locus mapping: comparison of methods for selected sibling pairs

During the past few years, there has been a great deal of new work on methods for mapping quantitative-trait loci by use of sibling pairs and sibships. There are several new methods based on linear regression, as well as several more that are based on score statistics. In theory, most of the new methods should be relatively robust to violations of distributional assumptions and to selected sampling, but, in practice, there has been little evaluation of how the methods perform on selected samples. We survey most of the new regression-based statistics and score statistics and propose a few minor variations on the score statistics. We use simulation to evaluate the type I error and the power of all of the statistics, considering (a) population samples of sibling pairs and (b) sibling pairs ascertained on the basis of at least one sibling with a trait value in the top 10% of the distribution. Most of the statistics have correct type I error for selected samples. The statistics proposed by Xu et al. and by Sham and Purcell are generally the most powerful, along with one of our score statistic variants. Even among the methods that are most powerful for "nice" data, some are more robust than others to non-Gaussian trait models and/or misspecified trait parameters.     

Selection and subsequent analysis of sib pair data for QTL detection.

Haseman and Elston (1972) developed a robust regression method for the detection of linkage between a marker and a quantitative trait locus (QTL) using sib pair data. The principle underlying this method is that the difference in phenotypes between pairs of sibs becomes larger as they share a decreasing number of alleles at a particular QTL identical by descent (IBD) from their parents. In this case, phenotypically very different sibs will also on average share a proportion of alleles IBD at any marker linked to the QTL that is lower than the expected value of 0.5. Thus, the deviation of the proportion of marker alleles IBD from the expected value in pairs of sibs selected to be phenotypically different (i.e. discordant) can provide a test for the presence of a QTL. A simple regression method for QTL detection in sib pairs selected for high phenotypic differences is presented here. The power of the analytical method was found to be greater than the power obtained using the standard analysis when samples of sib pairs with high phenotypic differences were used. However, the use of discordant sib pairs was found to be less powerful for QTL detection than alternative selective genotyping schemes based on the phenotypic values of the sibs except with intense selection, when its advantage was only marginal. The most effective selection scheme overall was the use of sib pairs from entire families selected on the basis of high within-family variance for the trait in question. There is little effect of selection on QTL position estimates, which are in good agreement with the simulated values. However, QTL variance estimates are biased to a greater or lesser degree, depending on the selection method.     

Combined linkage and association mapping of quantitative trait loci by multiple markers

Using multiple diallelic markers, variance component models are proposed for high-resolution combined linkage and association mapping of quantitative trait loci (QTL) based on nuclear families. The objective is to build a model that may fully use marker information for fine association mapping of QTL in the presence of prior linkage. The measures of linkage disequilibrium and the genetic effects are incorporated in the mean coefficients and are decomposed into orthogonal additive and dominance effects. The linkage information is modeled in variance-covariance matrices. Hence, the proposed methods model both association and linkage in a unified model. On the basis of marker information, a multipoint interval mapping method is provided to estimate the proportion of allele sharing identical by descent (IBD) and the probability of sharing two alleles IBD at a putative QTL for a sib-pair. To test the association between the trait locus and the markers, both likelihood-ratio tests and F-tests can be constructed on the basis of the proposed models. In addition, analytical formulas of noncentrality parameter approximations of the F-test statistics are provided. Type I error rates of the proposed test statistics are calculated to show their robustness. After comparing with the association between-family and association within-family (AbAw) approach by Abecasis and Fulker et al., it is found that the method proposed in this article is more powerful and advantageous based on simulation study and power calculation. By power and sample size comparison, it is shown that models that use more markers may have higher power than models that use fewer markers. The multiple-marker analysis can be more advantageous and has higher power in fine mapping QTL. As an application, the Genetic Analysis Workshop 12 German asthma data are analyzed using the proposed methods.     

The power to detect linkage disequilibrium with quantitative traits in selected samples

Results from power studies for linkage detection have led to many ongoing and planned collections of phenotypically extreme nuclear families. Given the great expense of collecting these families and the imminent availability of a dense diallelic marker map, the families are likely to be used in allelic-association as well as linkage studies. However, optimal selection strategies for linkage may not be equally powerful for association. We examine the power to detect linkage disequilibrium for quantitative traits after phenotypic selection. The results encompass six selection strategies that are in widespread use, including single selection (two designs), affected sib pairs, concordant and discordant pairs, and the extreme-concordant and -discordant design. Selection of sibships on the basis of one extreme proband with high or low trait scores provides as much power as discordant sib pairs but requires the screening and phenotyping of substantially fewer initial families from which to select. Analysis of the role of allele frequencies within each selection design indicates that common trait alleles generally offer the most power, but similarities between the marker- and trait-allele frequencies are much more important than the trait-locus frequency alone. Some of the most widespread selection designs, such as single selection, yield power gains only when both the marker and quantitative trait loci (QTL) are relatively rare in the population. In contrast, discordant pairs and the extreme-proband design provide power for the broadest range of QTL-marker-allele frequency differences. Overall, proband selection from either tail provides the best balance of power, robustness, and simplicity of ascertainment for family-based association analysis.     

Theoretical and empirical power of regression and maximum-likelihood methods to map quantitative trait loci in general pedigrees

Both theoretical calculations and simulation studies have been used to compare and contrast the statistical power of methods for mapping quantitative trait loci (QTLs) in simple and complex pedigrees. A widely used approach in such studies is to derive or simulate the expected mean test statistic under the alternative hypothesis of a segregating QTL and to equate a larger mean test statistic with larger power. In the present study, we show that, even when the test statistic under the null hypothesis of no linkage follows a known asymptotic distribution (the standard being chi(2)), it cannot be assumed that the distribution under the alternative hypothesis is noncentral chi(2). Hence, mean test statistics cannot be used to indicate power differences, and a comparison between methods that are based on simulated average test statistics may lead to the wrong conclusion. We illustrate this important finding, through simulations and analytical derivations, for a recently proposed new regression method for the analysis of general pedigrees to map quantitative trait loci. We show that this regression method is not necessarily more powerful nor computationally more efficient than a maximum-likelihood variance-component approach. We advocate the use of empirical power to compare trait-mapping methods.     

Strategies for mapping heterogeneous recessive traits by allele-sharing methods.

We investigate strategies for detecting linkage of recessive and partially recessive traits, using sibling pairs and inbred individuals. We assume that a genomewide search is being conducted and that locus heterogeneity of the trait is likely. For sibling pairs, we evaluate the efficiency of different statistics under the assumption that one does not know the true degree of recessiveness of the trait. We recommend a sibling-pair statistic that is a linear compromise between two previously suggested statistics. We also compare the power of sibling pairs to that of more distant relatives, such as cousins. For inbred individuals, we evaluate the power of offspring of different types of matings and compare them to sibling pairs. Over a broad range of trait etiologies, sibling pairs are more powerful than inbred individuals, but for traits caused by very rare alleles, particularly in the case of heterogeneity, inbred individuals can be much more powerful. The models we develop can also be used to examine specific situations other than those we look at. We present this analysis in the idealized context of a dense set of highly polymorphic markers. In general, incorporation of real-world complexities makes inbred individuals, particularly offspring of distant relatives, look slightly less useful than our results imply.     

Identity-by-descent estimation and mapping of qualitative traits in large, complex pedigrees

Computing identity-by-descent sharing between individuals connected through a large, complex pedigree is a computationally demanding task that often cannot be done using exact methods. What I present here is a rapid computational method for estimating, in large complex pedigrees, the probability that pairs of alleles are IBD given the single-point genotype data at that marker for all individuals. The method can be used on pedigrees of essentially arbitrary size and complexity without the need to divide the individuals into separate subpedigrees. I apply the method to do qualitative trait linkage mapping using the nonparametric sharing statistic S(pairs). The validity of the method is demonstrated via simulation studies on a 13-generation 3028-person pedigree with 700 genotyped individuals. An analysis of an asthma data set of individuals in this pedigree finds four loci with P-values <10(-3) that were not detected in prior analyses. The mapping method is fast and can complete analyses of approximately 150 affected individuals within this pedigree for thousands of markers in a matter of hours.     

Mapping quantitative trait loci with extreme discordant sib pairs: sampling considerations.

Elsewhere we have proposed the use of extreme discordant sib pairs (EDSPs) for mapping quantitative trait loci in humans. Here we present sample sizes necessary to achieve a given level of power with this study design, as well as the number of sibs that need to be screened to obtain the required sample. Further, we present simple formulas for adjusting sample sizes to account for variable significance levels and power, as well as the density and informativeness of linkage markers in a multipoint sib-pair analysis. We conclude that with EDSPs, the most powerful study design, the smallest genetic effect detectable with a realistic sample size is approximately 10% of the variance of the trait.     

Nonparametric trend statistic incorporating dispersion differences in sib pair linkage for quantitative traits

The Haseman-Elston (HE) regression method and its extensions are widely used in genetic studies for detecting linkage to quantitative trait loci (QTL) using sib pairs. The principle underlying the simple HE regression method is that the similarity in phenotypes between two siblings increases as they share an increasing number of alleles identical by descent (IBD) from their parents at a particular marker locus. In such a procedure, similarity was identified with the locations, that is, means of groups of sib pairs sharing 0, 1, and 2 alleles IBD. A more powerful, rank-based nonparametric test to detect increasing similarity in sib pairs is presented by combining univariate trend statistics not only of locations, but also of dispersions of the squared phenotypic differences of two siblings for three groups. This trend test does not rely on distributional assumptions, and is applicable to the skewed or leptokurtic phenotypic distributions, in addition to normal or near normal phenotypic distributions. The performances of nonparametric trend statistics, including nonparametric regression slope, are compared with the HE regression methods as genetic linkage strategies.     

On the use of double haploids for detecting QTL in outbred populations

The power to detect quantitative trait loci (QTL) using the double haploid (DH), full-sib (FS) and hierarchical (HI) designs implemented in outbred fish populations was assessed for interval mapping using deterministic methods. The predictions were tested using simulation. The DH design was most efficient for the range of designs and parameters considered and was most beneficial when the FS design was not very powerful. The difference between the design was largest for a low amount of residual genetic variation. Accounting for an increase of the environmental variance due to the genetic constitution of the double haploid progeny changed the magnitude of the power, but the ranking of the designs remained the same. As large full sib family sizes can be obtained in fish, the practical value of HI designs as a strategy for increasing the power of QTL mapping experiments is limited when compared with the FS design. Overall, the results suggested that the DH design could be a very useful tool for QTL mapping in fish, and of particular importance when the effect of the QTL is low and the residual genetic variation from other chromosomes can be controlled by using multiple markers.     

Haplotype sharing analysis using mantel statistics

OBJECTIVE: The potential value of haplotypes has attracted widespread interest in the mapping of complex traits. Haplotype sharing methods take the linkage disequilibrium information between multiple markers into account, and may have good power to detect predisposing genes. We present a new approach based on Mantel statistics for spacetime clustering, which is developed in order to improve the power of haplotype sharing analysis for gene mapping in complex disease. METHODS: The new statistic correlates genetic similarity and phenotypic similarity across pairs of haplotypes for case-only and case-control studies. The genetic similarity is measured as the shared length between haplotypes around a putative disease locus. The phenotypic similarity is measured as the mean-corrected cross-product based on the respective phenotypes. We analyzed two tests for statistical significance with respect to type I error: (1) assuming asymptotic normality, and (2) using a Monte Carlo permutation procedure. The results were compared to the chi(2) test for association based on 3-marker haplotypes. RESULTS: The results of the type I error rates for the Mantel statistics using the permutational procedure yielded pointwise valid tests. The approach based on the assumption of asymptotic normality was seriously liberal. CONCLUSION: Power comparisons showed that the Mantel statistics were better than or equal to the chi(2) test for all simulated disease models.