Genotypic dispersion and covariance between relatives under linkage disequilibrium

A random mating diploid population under linkage disequilibrium is considered. In the case of two diallelic loci, the problem about condition and joint distributions of genotypes of relatives being in arbitrary genetic relations is solved. Formulae of the partitioning of genotypic variance and covariance between relatives with respect to a polygenic character are inferred (in the case of many characters - of genotypic covariance matrix).     

The trimmed-haplotype test for linkage disequilibrium

Single-marker linkage-disequilibrium (LD) methods cannot fully describe disequilibrium in an entire chromosomal region surrounding a disease allele. With the advent of myriad tightly linked microsatellite markers, we have an opportunity to extend LD analysis from single markers to multiple-marker haplotypes. Haplotype analysis has increased statistical power to disclose the presence of a disease locus in situations where it correctly reflects the historical process involved. For maximum efficiency, evidence of LD ought to come not just from a single haplotype, which may well be rare, but in addition from many similar haplotypes that could have descended from the same ancestral founder but have been trimmed in succeeding generations. We present such an analysis, called the "trimmed-haplotype method." We focus on chromosomal regions that are small enough that disequilibrium in significant portions of them may have been preserved in some pedigrees and yet that contain enough markers to minimize coincidental occurrence of the haplotype in the absence of a disease allele: perhaps regions 1-2 cM in length. In general, we could have no idea what haplotype an ancestral founder carried generations ago, nor do we usually have a precise chromosomal location for the disease-susceptibility locus. Therefore, we must search through all possible haplotypes surrounding multiple locations. Since such repeated testing obliterates the sampling distribution of the test, we employ bootstrap methods to calculate significance levels. Trimmed-haplotype analysis is performed on family data in which genotypes have been assembled into haplotypes. It can be applied either to conventional parent-affected-offspring triads or to multiplex pedigrees. We present a method for summarizing the LD evidence, in any pedigree, that can be employed in trimmed-haplotype analysis as well as in other methods.     

Inferences about linkage disequilibrium.

Existing theory for inferences about linkage disequilibrium is restricted to a measure defined on gametic frequencies. Unless gametic frequencies are directly observable, they are inferred from genotypic frequencies under the assumption of random union of gametes. Primary emphasis in this paper is given to genotypic data, and disequilibrium coefficients are defined for all subsets of two or more of the four genes, two at each of two loci, carried by an individual. Linkage disequilibrium coefficients are defined for genes within and between gametes, and methods of estimating and testing these coefficients are given for gametic data. For genotypic data, when coupling and repulsion double heterozygotes cannot be distinguished. Burrows' composite measure of linkage disequilibrium is discussed. In particular, the estimate for this measure and hypothesis tests based on it are compared to the usual maximum likelihood estimate of gametic linkage disequilibrium, and corresponding likelihood ratio or contingency chi-square tests. General use of the composite measure, whether or not random union of gametes is an appropriate assumption, is recommended. Attention is given to small samples, where the non-normality of gene frequencies will have greatest effect on methods of inference based on normal theory. Even tools such as Fisher's z-transformation for the correlation of gene frequencies are found to perform quite satisfactorily.     

Robust transmission/disequilibrium test for incomplete family genotypes

Several solutions have been proposed to extend the transmission disequilibrium test (TDT) to include cases with missing parental genotype. However, completion of the missing parental genotype may bias the test if the underlying missing data mechanism is informative. Furthermore, all these solutions resolve the problem of missing parental genotype, while offspring with missing genotypes are typically ignored. We propose here an extension to the TDT, called robust TDT (rTDT), able to handle incomplete genotypes on both parents and children and that does not rest on any assumption about the missing data mechanism. rTDT returns minimum and maximum values of TDT that are consistent with all the possible completions of the missing data. We also show that, in some situations, rTDT can achieve both greater power and greater significance than the popular TDT analysis of incomplete data. rTDT is applied to a database of markers of susceptibility to Crohn's disease and it shows that only 2 of the 11 markers originally associated with the phenotype do not depend on assumptions about the missing data mechanism.     

Incorporating interference into linkage analysis for experimental crosses

The phenomenon of interference in genetic recombination is well-known and studied in a wide variety of organisms. Multilocus linkage analysis, which makes use of recombination patterns among all genetic markers simultaneously, is routinely used with data on humans and experimental organisms to build genetic maps. It is also used to try to determine the genes involved in traits of interest, such as common diseases. Most linkage analyses performed today ignore the occurrence of genetical interference. We present an extension to the Lander-Green algorithm for experimental crosses (backcross and intercross) to incorporate crossover interference according to the chi2 model. Simulation results show the impact of using this model on the accuracy of estimated genetic maps.     

The extent of linkage disequilibrium in beef cattle breeds using high-density SNP genotypes

Background

The extent of linkage disequilibrium (LD) between molecular markers impacts genome-wide association studies and implementation of genomic selection. The availability of high-density single nucleotide polymorphism (SNP) genotyping platforms makes it possible to investigate LD at an unprecedented resolution. In this work, we characterised LD decay in breeds of beef cattle of taurine, indicine and composite origins and explored its variation across autosomes and the X chromosome.

Findings

In each breed, LD decayed rapidly and r² was less than 0.2 for marker pairs separated by 50 kb. The LD decay curves clustered into three groups of similar LD decay that distinguished the three main cattle types. At short distances between markers (< 10 kb), taurine breeds showed higher LD (r² = 0.45) than their indicine (r² = 0.25) and composite (r² = 0.32) counterparts. This higher LD in taurine breeds was attributed to a smaller effective population size and a stronger bottleneck during breed formation. Using all SNPs on only the X chromosome, the three cattle types could still be distinguished. However for taurine breeds, the LD decay on the X chromosome was much faster and the background level much lower than for indicine breeds and composite populations. When using only SNPs that were polymorphic in all breeds, the analysis of the X chromosome mimicked that of the autosomes.

Conclusions

The pattern of LD mirrored some aspects of the history of breed populations and showed a sharp decay with increasing physical distance between markers. We conclude that the availability of the HD chip can be used to detect association signals that remained hidden when using lower density genotyping platforms, since LD dropped below 0.2 at distances of 50 kb.     

Linkage detection adaptive to linkage disequilibrium: the disequilibrium maximum-likelihood-binomial test for affected-sibship data

It has been demonstrated in the literature that the transmission/disequilibrium test (TDT) has higher power than the affected-sib-pair (ASP) mean test when linkage disequilibrium (LD) is strong but that the mean test has higher power when LD is weak. Thus, for ASP data, it seems clear that the TDT should be used when LD is strong but that the mean test or other linkage tests should be used when LD is weak or absent. However, in practice, it may be difficult to follow such a guideline, because the extent of LD is often unknown. Even with a highly dense genetic-marker map, in which some markers should be located near the disease-predisposing mutation, strong LD is not inevitable. Besides the genetic distance, LD is also affected by many factors, such as the allelic heterogeneity at the disease locus, the initial LD, the allelic frequencies at both disease locus and marker locus, and the age of the mutation. Therefore, it is of interest to develop methods that are adaptive to the extent of LD. In this report, we propose a disequilibrium maximum-binomial-likelihood (DMLB) test that incorporates LD in the maximum-binomial-likelihood (MLB) test. Examination of the corresponding score statistics shows that this method adaptively combines two sources of information: (a) the identity-by-descent (IBD) sharing score, which is informative for linkage regardless of the existence of LD, and (b) the contrast between allele-specific IBD sharing score, which is informative for linkage only in the presence of LD. For ASP data, the proposed test has higher power than either the TDT or the mean test when the extent of LD ranges from moderate to strong. Only when LD is very weak or absent is the DMLB slightly less powerful than the mean test; in such cases, the TDT has essentially no power to detect linkage. Therefore, the DMLB test is an interesting approach to linkage detection when the extent of LD is unknown.     

Improved power by use of a weighted score test for linkage disequilibrium mapping

Association studies offer an exciting approach to finding underlying genetic variants of complex human diseases. However, identification of genetic variants still includes difficult challenges, and it is important to develop powerful new statistical methods. Currently, association methods may depend on single-locus analysis--that is, analysis of the association of one locus, which is typically a single-nucleotide polymorphism (SNP), at a time--or on multilocus analysis, in which multiple SNPs are used to allow extraction of maximum information about linkage disequilibrium (LD). It has been shown that single-locus analysis may have low power because a single SNP often has limited LD information. Multilocus analysis, which is more informative, can be performed on the basis of either haplotypes or genotypes. It may lose power because of the often large number of degrees of freedom involved. The ideal method must make full use of important information from multiple loci but avoid increasing the degrees of freedom. Therefore, we propose a method to capture information from multiple SNPs but with the use of fewer degrees of freedom. When a set of SNPs in a block are correlated because of LD, we might expect that the genotype variation among the different phenotypic groups would extend across all the SNPs, and this information could be compressed into the low-frequency components of a Fourier transform. Therefore, we develop a test based on weighted Fourier transformation coefficients, with more weight given to the low-frequency components. Our simulation results demonstrate the validity and substantially higher power of the proposed method compared with other common methods. This method provides an additional tool to existing methods for identification of causative genetic variants underlying complex diseases.     

Insights into recombination from patterns of linkage disequilibrium in humans

An ability to predict levels of linkage disequilibrium (LD) between linked markers would facilitate the design of association studies and help to distinguish between evolutionary models. Unfortunately, levels of LD depend crucially on the rate of recombination, a parameter that is difficult to measure. In humans, rates of genetic exchange between markers megabases apart can be estimated from a comparison of genetic and physical maps; these large-scale estimates can then be interpolated to predict LD at smaller ("local") scales. However, if there is extensive small-scale heterogeneity, as has been recently proposed, local rates of recombination could differ substantially from those averaged over much larger distances. We test this hypothesis by estimating local recombination rates indirectly from patterns of LD in 84 genomic regions surveyed by the SeattleSNPs project in a sample of individuals of European descent and of African-Americans. We find that LD-based estimates are significantly positively correlated with map-based estimates. This implies that large-scale, average rates are informative about local rates of recombination. Conversely, although LD-based estimates are based on a number of simplifying assumptions, it appears that they capture considerable information about the underlying recombination rate or at least about the ordering of regions by recombination rate. Using LD-based estimators, we also find evidence for homologous gene conversion in patterns of polymorphism. However, as we demonstrate by simulation, inferences about gene conversion are unreliable, even with extensive data from homogeneous regions of the genome, and are confounded by genotyping error.     

Expected behavior of conditional linkage disequilibrium.

The ubiquitousness of RFLPs in the human genome has greatly helped the mapping of human disease genes, and it has been suggested that population measures of association between disease and marker loci could help with this mapping. For rare diseases, random samples are taken from within disease genotypes in order to obtain reasonable sample sizes, but this sampling strategy requires a modification of the usual measures of association. We present theoretical predictions for the mean and variance of such a modified measure, under the assumption that the disease gene is maintained at a constant low frequency in the population. The coefficient of variation of this modified measure is large enough that caution is needed in using the measure to locate disease genes, and, furthermore, the coefficient of variation cannot be made arbitrarily small by increasing sample size. The modified association measure is calculated for recently published data on cystic fibrosis.     

The transmission/disequilibrium test for linkage on the X chromosome

The transmission/disequilibrium test (TDT), which detects linkage between a marker and disease loci in the presence of linkage disequilibrium, was introduced by Spielman et al. The original TDT requires families in which the genotypes are known for both parents and for at least one affected offspring, and this limits its applicability to diseases with late onset. The sib-TDT, or S-TDT, which utilizes families with affected and unaffected siblings, was introduced as an alternative method, by Spielman and Ewens, and the TDT and S-TDT can be combined in an overall test (i.e., a combined-TDT, or C-TDT). The TDT statistics described so far are for autosomal chromosomes. We have extended these TDT methods to test for linkage between X-linked markers and diseases that affect either males only or both sexes. For diseases of late onset, when parental genotypes are often unavailable, the X-linkage C-TDT may allow for more power than is provided by the X-linkage TDT alone.     

A program to test linkage disequilibrium between loci in autotetraploid species

No software currently implements a test of linkage disequilibrium in autotetraploid species. We propose a program, LD4X that performs a Fisher's exact test between pairs of alleles at two loci. All combinations of alleles from two loci are treated in turn. If two alleles of a pair of loci have a nonrandom distribution, the markers are considered as linked. The program was tested on a set of microsatellite markers in synthetic alfalfa populations.     

Detecting linkage disequilibrium in bacterial populations.

The distribution of the number of pairwise differences calculated from comparisons between n haploid genomes has frequently been used as a starting point for testing the hypothesis of linkage equilibrium. For this purpose the variance of the pairwise differences, VD, is used as a test statistic to evaluate the null hypothesis that all loci are in linkage equilibrium. The problem is to determine the critical value of the distribution of VD. This critical value can be estimated either by Monte Carlo simulation or by assuming that VD is distributed normally and calculating a one-tailed 95% critical value for VD, L, L = EVD + 1.645 sqrt(VarVD), where E(VD) is the expectation of VD, and Var(VD) is the variance of VD. If VD (observed) > L, the null hypothesis of linkage equilibrium is rejected. Using Monte Carlo simulation we show that the formula currently available for Var(VD) is incorrect, especially for genetically highly diverse data. This has implications for hypothesis testing in bacterial populations, which are often genetically highly diverse. For this reason we derive a new, exact formula for Var(VD). The distribution of VD is examined and shown to approach normality as the sample size increases. This makes the new formula a useful tool in the investigation of large data sets, where testing for linkage using Monte Carlo simulation can be very time consuming. Application of the new formula, in conjunction with Monte Carlo simulation, to populations of Bradyrhizobium japonicum, Rhizobium leguminosarum, and Bacillus subtilis reveals linkage disequilibrium where linkage equilibrium has previously been reported.     

A test for linkage and association in general pedigrees: the pedigree disequilibrium test

Family-based tests of linkage disequilibrium typically are based on nuclear-family data including affected individuals and their parents or their unaffected siblings. A limitation of such tests is that they generally are not valid tests of association when data from related nuclear families from larger pedigrees are used. Standard methods require selection of a single nuclear family from any extended pedigrees when testing for linkage disequilibrium. Often data are available for larger pedigrees, and it would be desirable to have a valid test of linkage disequilibrium that can use all potentially informative data. In this study, we present the pedigree disequilibrium test (PDT) for analysis of linkage disequilibrium in general pedigrees. The PDT can use data from related nuclear families from extended pedigrees and is valid even when there is population substructure. Using computer simulations, we demonstrated validity of the test when the asymptotic distribution is used to assess the significance, and examined statistical power. Power simulations demonstrate that, when extended pedigree data are available, substantial gains in power can be attained by use of the PDT rather than existing methods that use only a subset of the data. Furthermore, the PDT remains more powerful even when there is misclassification of unaffected individuals. Our simulations suggest that there may be advantages to using the PDT even if the data consist of independent families without extended family information. Thus, the PDT provides a general test of linkage disequilibrium that can be widely applied to different data structures.     

A sibship test for linkage in the presence of association: the sib transmission/disequilibrium test.

Linkage analysis with genetic markers has been successful in the localization of genes for many monogenic human diseases. In studies of complex diseases, however, tests that rely on linkage disequilibrium (the simultaneous presence of linkage and association) are often more powerful than those that rely on linkage alone. This advantage is illustrated by the transmission/disequilibrium test (TDT). The TDT requires data (marker genotypes) for affected individuals and their parents; for some diseases, however, data from parents may be difficult or impossible to obtain. In this article, we describe a method, called the "sib TDT" (or "S-TDT"), that overcomes this problem by use of marker data from unaffected sibs instead of from parents, thus allowing application of the principle of the TDT to sibships without parental data. In a single collection of families, there might be some that can be analyzed only by the TDT and others that are suitable for analysis by the S-TDT. We show how all the data may be used jointly in one overall TDT-type procedure that tests for linkage in the presence of association. These extensions of the TDT will be valuable for the study of diseases of late onset, such as non-insulin-dependent diabetes, cardiovascular diseases, and other diseases associated with aging.     

The posterior probability of linkage allowing for linkage disequilibrium and a new estimate of disequilibrium between a trait and a marker

The posterior probability of linkage (PPL) statistic has been developed as a method for the rigorous accumulation of evidence for or against linkage allowing for both intra- and inter-sample heterogeneity. To date, the method has assumed linkage equilibrium between alleles at the trait locus and the marker locus. We now generalize the PPL to allow for linkage disequilibrium (LD), by incorporating variable phase probabilities into the underlying linkage likelihood. This enables us to recover the marginal posterior density of the recombination fraction, integrating out nuisance parameters of the trait model, including the locus heterogeneity (admixture) parameter, as well as a vector of LD parameters. The marginal posterior density can then be updated across data subsets or new data as they become available, while allowing parameters of the trait model to vary between data sets. The method applies immediately to general pedigree structures and to markers with multiple alleles. In the case of SNPs, the likelihood is parameterized in terms of the standard single LD parameter D'; and it therefore affords a mechanism for estimation of D' between the marker and the trait, again, without fixing the parameters of the trait model and allowing for updating across data sets. It is even possible to allow for a different associated allele in different populations, while accumulating information regarding the strength of LD. While a computationally efficient implementation for multi-allelic markers is still in progress, we have implemented a version of this new LD-PPL for SNPs and evaluated its performance in nuclear families. Our simulations show that LD-PPLs tend to be larger than PPLs (stronger evidence in favor of linkage/LD) with increased LD level, under a variety of generating models; while in the absence of linkage and LD, LD-PPLs tend to be smaller than PPLs (stronger evidence against linkage). The estimate of D' also behaves well even in relatively small, heterogeneous samples.     

The quantitative trait linkage disequilibrium test: a more powerful alternative to the quantitative transmission disequilibrium test for use in the absence of population stratification

Linkage analysis based on identity-by-descent allele-sharing can be used to identify a chromosomal region harboring a quantitative trait locus (QTL), but lacks the resolution required for gene identification. Consequently, linkage disequilibrium (association) analysis is often employed for fine-mapping. Variance-components based combined linkage and association analysis for quantitative traits in sib pairs, in which association is modeled as a mean effect and linkage is modeled in the covariance structure has been extended to general pedigrees (quantitative transmission disequilibrium test, QTDT). The QTDT approach accommodates data not only from parents and siblings, but also from all available relatives. QTDT is also robust to population stratification. However, when population stratification is absent, it is possible to utilize even more information, namely the additional information contained in the founder genotypes. In this paper, we introduce a simple modification of the allelic transmission scoring method used in the QTDT that results in a more powerful test of linkage disequilibrium, but is only applicable in the absence of population stratification. This test, the quantitative trait linkage disequilibrium (QTLD) test, has been incorporated into a new procedure in the statistical genetics computer package SOLAR. We apply this procedure in a linkage/association analysis of an electrophysiological measurement previously shown to be related to alcoholism. We also demonstrate by simulation the increase in power obtained with the QTLD test, relative to the QTDT, when a true association exists between a marker and a QTL.