Weighting affected sib pairs by marker informativity

For the analysis of affected sib pairs (ASPs), a variety of test statistics is applied in genomewide scans with microsatellite markers. Even in multipoint analyses, these statistics might not fully exploit the power of a given sample, because they do not account for incomplete informativity of an ASP. For meta-analyses of linkage and association studies, it has been shown recently that weighting by informativity increases statistical power. With this idea in mind, the first aim of this article was to introduce a new class of tests for ASPs that are based on the mean test. To take into account how much informativity an ASP contributes, we weighted families inversely proportional to their marker informativity. The weighting scheme is obtained by use of the de Finetti representation of the distribution of identity-by-descent values. We derive the limiting distribution of the weighted mean test and demonstrate the validity of the proposed test. We show that it can be much more powerful than the classical mean test in the case of low marker informativity. In the second part of the article, we propose a Monte Carlo simulation approach for evaluating significance among ASPs. We demonstrate the validity of the simulation approach for both the classical and the weighted mean test. Finally, we illustrate the use of the weighted mean test by reanalyzing two published data sets. In both applications, the maximum LOD score of the weighted mean test is 0.6 higher than that of the classical mean test.     

Determining informativity of marker typing for genetic counseling in a pedigree

Summary For the situation of a Mendelian disease linked to a genetic marker, a new method is described that allows evaluating for genetic counseling the information potentially available from the linked marker before the marker data are actually obtained, that is, prior to drawing blood for marker typing. For a consultand in a family pedigree, the method determines the risk distribution (small families) or an approximation to it (larger families) and calculates the probability that the risk will deviate beyond certain limits from its a priori value, which exists without marker data, for example, that the risk will be smaller than 0.10 or larger than 0.90. The method was applied here to a pedigree of 15 individuals for which analytical calculations would be difficult to carry out.     

A Bayesian approach for applying Haseman-Elston methods

The main goal of this paper is to couple the Haseman-Elston method with a simple yet effective Bayesian factor-screening approach. This approach selects markers by considering a set of multigenic models that include epistasis effects. The markers are ranked based on their marginal posterior probability. A significant improvement over our previously proposed Bayesian variable selection methodology is a simple Metropolis-Hasting algorithm that requires minimum tuning on the prior settings. The algorithm, however, is also flexible enough for us to easily incorporate our hypotheses and avoid computational pitfalls. We apply our approach to the microsatellite data of Collaborative Studies on Genetics of Alcoholism using the coded values for the ALDX1 variable as our response.     

Information perspectives of the Haseman-Elston method

The Haseman-Elston regression approach can be viewed as a special case in a larger framework in which independent sibling pairs are ascertained, perhaps based on their phenotypes. Using likelihood arguments in an idealized setting, expressions are obtained for the inherent linkage information contained in the phenotype distributions. The results provide insight and analytic results that apply to general phenotype distributions and can be used to judge the effects of various data transformations and selective sampling of extreme phenotypes.     

Matrix representation of the Haseman-Elston method.

The Hasemann-Elston method of linkage detection is based on the probabilities of a sib pair having 0, 1, or 2 alleles identical by descent (IBD) at a marker and a trait locus. These probabilities form a 3x3 matrix. Here, the characteristic values and characteristic vectors of this matrix were used to clarify the structure of the equations and to simplify calculations. As examples, the regression coefficients were derived for three genetic systems: a trait and a marker, two epistatic traits and two markers, and one trait locus and two markers. The last model was studied under the assumption of no crossover interference, the expression for allele IBD sharing at a trait locus was derived as a function of allele IBD sharing at two marker loci, and the regression is shown to be non-linear.     

Extension of the Haseman-Elston regression model to longitudinal data

We propose an extension to longitudinal data of the Haseman and Elston regression method for linkage analysis. The proposed model is a mixed model having several random effects. As response variable, we investigate the sibship sample mean corrected cross-product (smHE) and the BLUP-mean corrected cross product (pmHE), comparing them with the original squared difference (oHE), the overall mean corrected cross-product (rHE), and the weighted average of the squared difference and the squared mean-corrected sum (wHE). The proposed model allows for the correlation structure of longitudinal data. Also, the model can test for gene x time interaction to discover genetic variation over time. The model was applied in an analysis of the Genetic Analysis Workshop 13 (GAW13) simulated dataset for a quantitative trait simulating systolic blood pressure. Independence models did not preserve the test sizes, while the mixed models with both family and sibpair random effects tended to preserve size well.     

Weighting improves the "new Haseman-Elston" method

Elston et al. [Genet Epidemiol, in press] apply the results of Wright [Am J Hum Genet 1997;60:740-742] and Drigalenko [Am J Hum Genet 1998;63:1242-1245] to extend the traditional Haseman-Elston regression scheme [Haseman and Elston, Behav Genet 1972;2:3-19] to include not only linkage information contained in the sib pair's squared difference, but also information in their mean-corrected squared sum. The new algorithm detects linkage to a quantitative trait locus by modelling sib pair trait covariance as a function of identity-by-descent status. We demonstrate why this new estimator is suboptimal and can in some cases be inferior to the original Haseman-Elston method. We also describe a simple approach to estimation which improves on this new Haseman-Elston method by incorporating variance-based weights into the test statistic while staying within the linear modelling framework. In support of our theoretical claim, we conduct both a sib pair simulation and an application to GAW 10 sib pair data showing that our new estimator is superior to both the old and new Haseman-Elston schemes currently implemented in the analysis package S.A.G.E. 4.0.     

A modified revisited Haseman-Elston method to further improve power

The original and revisited Haseman-Elston methods are simple robust methods to detect linkage, but neither is uniformly optimal in terms of power. In this report, we propose a simple modification of the revisited Haseman-Elston method that retains the simplicity and robustness properties, but increases its power. We demonstrate theoretically that the modification can be more powerful than the optimally weighted Haseman-Elston method when the sibship mean can be correctly specified. We then examine the properties of this modification by simulation when the sibship mean is replaced by its best linear unbiased predictor. The simulation results indicate that this modification maintains good control over type I error, even in the case of larger sibships, and that the empirical power of this modification is similar to that of the optimally weighted Haseman-Elston method in most cases.     

Testing the robustness of the new Haseman-Elston quantitative-trait loci-mapping procedure

Variance components (VC) techniques have emerged as among the more powerful methods for detection of quantitative-trait loci (QTL) in linkage analysis. Allison et al. found that, with particularly marked leptokurtosis in the phenotypic distribution and moderate-to-high residual sibling correlation, maximum likelihood (ML) VC methods may produce a severe excess of type I errors. The new Haseman-Elston (NHE) method is a least-squares-based VC method for mapping of QTL in sib pairs (Elston et al.). Using simulation, we investigate the robustness of the NHE to marked nonnormality, by means of the same distributions and worst-case conditions identified by Allison et al. for the ML approach (i.e., 100 pairs; high residual sibling correlation). Results showed that, when marked nonnormality is present, the NHE can be used without severe type I error-rate inflation, even at very small alpha levels.     

Equivalence between Haseman-Elston and variance-components linkage analyses for sib pairs

The Haseman-Elston regression method offers a simpler alternative to variance-components (VC) models, for the linkage analysis of quantitative traits. However, even the "revisited" method, which uses the cross-product--rather than the squared difference--in sib trait values, is, in general, less powerful than VC models. In this report, we clarify the relative efficiencies of existing Haseman-Elston methods and show how a new Haseman-Elston method can be constructed to have power equivalent to that of VC models. This method uses as the dependent variable a linear combination of squared sums and squared differences, in which the weights are determined by the overall trait correlation between sibs in a population. We show how this method can be used for both the selection of maximally informative sib pairs for genotyping and the subsequent analysis of such selected samples.     

A unified Haseman-Elston method for testing linkage with quantitative traits

The Haseman and Elston (H-E) method uses a simple linear regression to model the squared trait difference of sib pairs with the shared allele identical by descent (IBD) at marker locus for linkage testing. Under this setting, the squared mean-corrected trait sum is also linearly related to the IBD sharing. However, the resulting slope estimate for either model is not efficient. In this report, we propose a simple linkage test that optimally uses information from the estimates of both models. We also demonstrate that the new test is more powerful than both the traditional one and the recently revisited H-E methods.     

X-linked extension of the revised Haseman-Elston algorithm for linkage analysis in sib pairs

Haseman and Elston (H-E) proposed a regression-based robust test of linkage between a marker and an autosomal quantitative trait locus, using the squared sib pair trait difference as a dependent variable and the proportion of alleles shared identical by descent by the sib pair as an independent variable. Several authors have proposed improvement of the original H-E's seminal work by using an optimal linear combination of squared sum and squared difference as the dependent variable. In this paper, we extend Haseman and Elston's sib pair method to an X-linked locus. We give a general formulation of the complete regression model and details of the regression coefficients in terms of variance components. Simulation results are presented to describe the power of this technique for a theoretical best case scenario.     

Analysis of heregulin symmetry by weighted evolutionary tracing

Heregulins are members of the protein family of EGF-like growth and differentiation factors. The primary cell-surface targets of heregulins are heterodimers of the EGF-receptor homolog HER2 with either HER3 or HER4. We used a weighted evolutionary trace analysis to identify structural features that distinguish the EGF-like domain (hrg) of heregulins from other members of the EGF family. In this analysis, each amino acid sequence is weighted according to its uniqueness and the variability in each position is assigned by an amino acid substitution matrix. Conserved residues in heregulin that are variable in other EGF-like domains are considered possible specificity-conferring residues. This analysis identifies two clusters of residues at the foot of the boot-shaped hrg domain. The residues in one cluster are recruited from the N-terminus; those in the other are from the ohm-loop region and show a weak sequence similarity to the N-terminal residues at the opposite side of the boot. The remaining residues with high conservation scores distribute themselves into these two distinct surfaces on hrg. This pseudo-twofold symmetry and the presence of two distinct interfaces may reflect the preference of hrg for heterodimeric versus homodimeric HER complexes.     

Novel RFLPs and microsatellite repeats increase informativity at four loci mapping to Xq22-q25

Three microsatellites and five restriction fragment length polymorphisms (RFLPs) have been detected at loci DXS275, DXS97, DXS11 and DXS100. All of the reported polymorphisms are in linkage equilibrium with existing polymorphisms at these loci. However, two RFLPs at DXS275 and two at DXS97 appear to be in linkage disequilibrium with each other. Increased informativity at these loci have enabled exclusion of Xq22-q25 as a candidate region for X-linked spina bifida and anencephaly, and should aid in the mapping of other genes in the region.     

Comparative phylogenetic informativity of single ribosomal cluster regions in freshwater horsehair worms (Gordiacea,Nematomorpha)

The informativity of some regions of the ribosomal cluster (rDNA) in freshwater horsehair worms (Nematomorpha) has been analyzed to trace the nucleotide differences for discrimination of species and genera of these parasites using available sequences from the NCBI GenBank and our data.