Transmission/disequilibrium test based on haplotype sharing for tightly linked markers

Studies using haplotypes of multiple tightly linked markers are more informative than those using a single marker. However, studies based on multimarker haplotypes have some difficulties. First, if we consider each haplotype as an allele and use the conventional single-marker transmission/disequilibrium test (TDT), then the rapid increase in the degrees of freedom with an increasing number of markers means that the statistical power of the conventional tests will be low. Second, the parental haplotypes cannot always be unambiguously reconstructed. In the present article, we propose a haplotype-sharing TDT (HS-TDT) for linkage or association between a disease-susceptibility locus and a chromosome region in which several tightly linked markers have been typed. This method is applicable to both quantitative traits and qualitative traits. It is applicable to any size of nuclear family, with or without ambiguous phase information, and it is applicable to any number of alleles at each of the markers. The degrees of freedom (in a broad sense) of the test increase linearly as the number of markers considered increases but do not increase as the number of alleles at the markers increases. Our simulation results show that the HS-TDT has the correct type I error rate in structured populations and that, in most cases, the power of HS-TDT is higher than the power of the existing single-marker TDTs and haplotype-based TDTs.     

Sample reproducibility of genetic association using different multimarker TDTs in genome-wide association studies: characterization and a new approach

Multimarker Transmission/Disequilibrium Tests (TDTs) are very robust association tests to population admixture and structure which may be used to identify susceptibility loci in genome-wide association studies. Multimarker TDTs using several markers may increase power by capturing high-degree associations. However, there is also a risk of spurious associations and power reduction due to the increase in degrees of freedom. In this study we show that associations found by tests built on simple null hypotheses are highly reproducible in a second independent data set regardless the number of markers. As a test exhibiting this feature to its maximum, we introduce the multimarker 2-Groups TDT (mTDT(2G)), a test which under the hypothesis of no linkage, asymptotically follows a χ2 distribution with 1 degree of freedom regardless the number of markers. The statistic requires the division of parental haplotypes into two groups: disease susceptibility and disease protective haplotype groups. We assessed the test behavior by performing an extensive simulation study as well as a real-data study using several data sets of two complex diseases. We show that mTDT(2G) test is highly efficient and it achieves the highest power among all the tests used, even when the null hypothesis is tested in a second independent data set. Therefore, mTDT(2G) turns out to be a very promising multimarker TDT to perform genome-wide searches for disease susceptibility loci that may be used as a preprocessing step in the construction of more accurate genetic models to predict individual susceptibility to complex diseases.     

一种有效的复杂疾病基因定位的检测法

连锁不平衡（LD）应用于某些复杂疾病基因的定位,近年来发展了许多LD定位方法,除TDT外,大多数LD定位方法须先假定无人群混和,人群混合可增大在疾病基因定位时犯Ⅰ类错误的机率,产生无效结果。此方法利用LD来检测标记位点和疾病敏感位点（DSL）的连锁（有连锁不平衡）相关（有连锁）。分析时采用不相关样本,已知其父母基因型和至少父母之一为杂合子,再将随机样本依基因型不同分类,然后对来自不同类的数据应用有力的统计方法进行单独和联合分析。此LD定位法不仅适用于患病和正常个体,而且有效消除据父母基因分类的样本定位时人群混合的影响,分析结果和模拟结果也表明此方法解决了在检测标记位点和疾病敏感位点之间的连锁和相关时人群混和的问题,但与TDT比,此法在检测的位点为DSL时丙能有效和充分地利用矫正数据,检测位点不是DSL时,此法和TDT法可相互补充更有效地检测连锁的DSL。     

Detecting genome wide haplotype sharing using SNP or microsatellite haplotype data

Genome wide association studies using high throughput technology are already being conducted despite the significant hurdles that need to be overcome (Nat Rev Genet 6:95–108, 2005; Nat Rev Genet 6:109–118, 2005). Methods for detecting haplotype association signals in genome wide haplotype datasets are as yet very limited. Much methodological research has already been devoted to linkage disequilibrium (LD) fine mapping where the focus is the identification of the disease locus rather than the detection of a disease signal. Applications of these approaches to genome wide scanning are limited by the strong model assumptions of the sharing process, which lead to computational complexity. We describe a new algorithm for the initial identification of disease susceptibility loci in genome wide haplotype association studies. Excess sharing of ancestral haplotypes, which indicates the presence of a disease locus, is detected with a simple, easy to interpret, χ ² based statistic. The method allows genome wide scanning for qualitative traits within reasonable computational timeframes and can serve as a first pass analysis prior to the usage of likelihood based methods, providing candidate regions and inferred susceptibility haplotypes. Our method makes no assumptions regarding the population history or the pattern of background LD. Statistical significance is evaluated with permutation tests. The method is illustrated on simulated and real data where it is applied to simple (cystic fibrosis) and complex disease (multiple sclerosis) examples. The statistic has low type I error and greater power to map disease loci over conventional single marker tests for low to moderate levels of LD.     

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.     

Linkage-disequilibrium mapping of disease genes by reconstruction of ancestral haplotypes in founder populations.

Linkage disequilibrium (LD) mapping may be a powerful means for genome screening to identify susceptibility loci for common diseases. A new statistical approach for detection of LD around a disease gene is presented here. This method compares the distribution of haplotypes in affected individuals versus that expected for individuals descended from a common ancestor who carried a mutation of the disease gene. Simulations demonstrate that this method, which we term "ancestral haplotype reconstruction" (AHR), should be powerful for genome screening of phenotypes characterized by a high degree of etiologic heterogeneity, even with currently available marker maps. AHR is best suited to application in isolated populations where affected individuals are relatively recently descended (< approximately 25 generations) from a common disease mutation-bearing founder.     

Nonparametric disequilibrium mapping of functional sites using haplotypes of multiple tightly linked single-nucleotide polymorphism markers

As the speed and efficiency of genotyping single-nucleotide polymorphisms (SNPs) increase, using the SNP map, it becomes possible to evaluate the extent to which a common haplotype contributes to the risk of disease. In this study we propose a new procedure for mapping functional sites or regions of a candidate gene of interest using multiple linked SNPs. Based on a case-parent trio family design, we use expectation-maximization (EM) algorithm-derived haplotype frequency estimates of multiple tightly linked SNPs from both unambiguous and ambiguous families to construct a contingency statistic S for linkage disequilibrium (LD) analysis. In the procedure, a moving-window scan for functional SNP sites or regions can cover an unlimited number of loci except for the limitation of computer storage. Within a window, all possible widths of haplotypes are utilized to find the maximum statistic S* for each site (or locus). Furthermore, this method can be applied to regional or genome-wide scanning for determining linkage disequilibrium using SNPs. The sensitivity of the proposed procedure was examined on the simulated data set from the Genetic Analysis Workshop (GAW) 12. Compared with the conventional and generalized TDT methods, our procedure is more flexible and powerful.     

Transmission-disequilibrium tests for quantitative traits.

The transmission-disequilibrium test (TDT) of Spielman et al. is a family-based linkage-disequilibrium test that offers a powerful way to test for linkage between alleles and phenotypes that is either causal (i.e., the marker locus is the disease/trait allele) or due to linkage disequilibrium. The TDT is equivalent to a randomized experiment and, therefore, is resistant to confounding. When the marker is extremely close to the disease locus or is the disease locus itself, tests such as the TDT can be far more powerful than conventional linkage tests. To date, the TDT and most other family-based association tests have been applied only to dichotomous traits. This paper develops five TDT-type tests for use with quantitative traits. These tests accommodate either unselected sampling or sampling based on selection of phenotypically extreme offspring. Power calculations are provided and show that, when a candidate gene is available (1) these TDT-type tests are at least an order of magnitude more efficient than two common sib-pair tests of linkage; (2) extreme sampling results in substantial increases in power; and (3) if the most extreme 20% of the phenotypic distribution is selectively sampled, across a wide variety of plausible genetic models, quantitative-trait loci explaining as little as 5% of the phenotypic variation can be detected at the .0001 alpha level with <300 observations.     

Transmission/disequilibrium test meets measured haplotype analysis: family-based association analysis guided by evolution of haplotypes

Family data teamed with the transmission/disequilibrium test (TDT), which simultaneously evaluates linkage and association, is a powerful means of detecting disease-liability alleles. To increase the information provided by the test, various researchers have proposed TDT-based methods for haplotype transmission. Haplotypes indeed produce more-definitive transmissions than do the alleles comprising them, and this tends to increase power. However, the larger number of haplotypes, relative to alleles at individual loci, tends to decrease power, because of the additional degrees of freedom required for the test. An optimal strategy would focus the test on particular haplotypes or groups of haplotypes. In this report we develop such an approach by combining the theory of TDT with that of measured haplotype analysis (MHA). MHA uses the evolutionary relationships among haplotypes to produce a limited set of hypothesis tests and to increase the interpretability of these tests. The theory of our approach, called the "evolutionary tree" (ET)-TDT, is developed for two cases: when haplotype transmission is certain and when it is not. Simulations show the ET-TDT can be more powerful than other proposed methods under reasonable conditions. More importantly, our results show that, when multiple polymorphisms are found within the gene, the ET-TDT can be useful for determining which polymorphisms affect liability.     

Mapping multiple sclerosis susceptibility to the HLA-DR locus in African Americans

An underlying complex genetic susceptibility exists in multiple sclerosis (MS), and an association with the HLA-DRB1*1501-DQB1*0602 haplotype has been repeatedly demonstrated in high-risk (northern European) populations. It is unknown whether the effect is explained by the HLA-DRB1 or the HLA-DQB1 gene within the susceptibility haplotype, which are in strong linkage disequilibrium (LD). African populations are characterized by greater haplotypic diversity and distinct patterns of LD compared with northern Europeans. To better localize the HLA gene responsible for MS susceptibility, case-control and family-based association studies were performed for DRB1 and DQB1 loci in a large and well-characterized African American data set. A selective association with HLA-DRB1*15 was revealed, indicating a primary role for the DRB1 locus in MS independent of DQB1*0602. This finding is unlikely to be solely explained by admixture, since a substantial proportion of the susceptibility chromosomes from African American patients with MS displayed haplotypes consistent with an African origin.     

Identification of a psoriasis susceptibility candidate gene by linkage disequilibrium mapping with a localized single nucleotide polymorphism map

Psoriasis is a chronic inflammatory disease of the skin with both genetic and environmental risk factors. Here we describe the creation of a single-nucleotide polymorphism (SNP) map spanning 900-1200 kb of chromosome 3q21, which had been previously recognized as containing a psoriasis susceptibility locus, PSORS5. We genotyped 644 individuals, from 195 Swedish psoriatic families, for 19 polymorphisms. Linkage disequilibrium (LD) between marker and disease was assessed using the transmission/disequilibrium test (TDT). In the TDT analysis, alleles of three of these SNPs showed significant association with disease (P<0.05). A 160-kb interval encompassing these three SNPs was sequenced, and a coding sequence consisting of 13 exons was identified. The predicted protein shares 30-40% homology with the family of cation/chloride cotransporters. A five-marker haplotype spanning the 3' half of this gene is associated with psoriasis to a P value of 3.8<10(-5). We have called this gene SLC12A8, coding for a member of the solute carrier family 12 proteins. It belongs to a class of genes that were previously unrecognized as playing a role in psoriasis pathogenesis.     

Haplotypes at ATM identify coding-sequence variation and indicate a region of extensive linkage disequilibrium

Genetic variation in the human population may lead to functional variants of genes that contribute to risk for common chronic diseases such as cancer. In an effort to detect such possible predisposing variants, we constructed haplotypes for a candidate gene and tested their efficacy in association studies. We developed haplotypes consisting of 14 biallelic neutral-sequence variants that span 142 kb of the ATM locus. ATM is the gene responsible for the autosomal recessive disease ataxia-telangiectasia (AT). These ATM noncoding single-nucleotide polymorphisms (SNPs) were genotyped in nine CEPH families (89 individuals) and in 260 DNA samples from four different ethnic origins. Analysis of these data with an expectation-maximization algorithm revealed 22 haplotypes at this locus, with three major haplotypes having frequencies > or = .10. Tests for recombination and linkage disequilibrium (LD) show reduced recombination and extensive LD at the ATM locus, in all four ethnic groups studied. The most striking example was found in the study population of European ancestry, in which no evidence for recombination could be discerned. The potential of ATM haplotypes for detection of genetic variants through association studies was tested by analysis of 84 individuals carrying one of three ATM coding SNPs. Each coding SNP was detected by association with an ATM haplotype. We demonstrate that association studies with haplotypes for candidate genes have significant potential for the detection of genetic backgrounds that contribute to disease.     

Genetic dissection of the human leukocyte antigen region by use of haplotypes of Tasmanians with multiple sclerosis

Association of multiple sclerosis (MS) with the human leukocyte antigen (HLA) class II haplotype DRB1*1501-DQB1*0602 is the most consistently replicated finding of genetic studies of the disease. However, the high level of linkage disequilibrium (LD) in the HLA region has hindered the identification of other loci that single-marker tests for association are unlikely to resolve. In order to address this issue, we generated haplotypes spanning 14.754 Mb (5 cM) across the entire HLA region. The haplotypes, which were inferred by genotyping relatives of 152 patients with MS and 105 unaffected control subjects of Tasmanian ancestry, define a genomic segment from D6S276 to D6S291, including 13 microsatellite markers integrated with allele-typing data for DRB1 and DQB1. Association to the DRB1*1501-DQB1*0602 haplotype was replicated. In addition, we found that the class I/extended class I region, defined by a genomic segment of approximately 400 kb between MOGCA and D6S265, harbors genes that independently increase risk of, or provide protection from, MS. Log-linear modeling analysis of constituent haplotypes that represent genomic regions containing class I (MOGCA-D6S265), class III (TNFa-TNFd-D6S273), and class II (DRB1-DQB1) genes indicated that having class I and class II susceptibility variants on the same haplotype provides an additive effect on risk. Moreover, we found no evidence for a disease locus in the class III region defined by a 150-kb genomic segment containing the TNF locus and 14 other genes. A global overview of LD performed using GOLD identified two discrete blocks of LD in the HLA region that correspond well with previous findings. We propose that the analysis of haplotypes, by use of the types of approaches outlined in the present article, should make it possible to more accurately define the contribution of the HLA to MS.     

Localization of psoriasis-susceptibility locus PSORS1 to a 60-kb interval telomeric to HLA-C

Recent genome scans have established the presence of a major psoriasis-susceptibility locus in the human leukocyte antigen (HLA) complex on chromosome 6p21.3. To narrow the interval for candidate gene testing, we performed a linkage-disequilibrium analysis of 339 families, with the use of 62 physically mapped microsatellite markers spanning the major histocompatibility complex (MHC). As detected by use of the transmission/disequilibrium test (TDT), individual markers yielded significant linkage disequilibrium across most of the MHC. However, the strongest evidence for marker-trait disequilibrium was found in an approximately 300-kb region extending from the MICA gene to the corneodesmosin gene. Maximum-likelihood haplotypes were constructed across the entire MHC in the original sample and across a 1.2-Mb region of the central MHC in an expanded sample containing 139 additional families. Short (two- to five-marker) haplotypes were subjected to the TDT using a "moving-window" strategy that reduced the variability of TDT P values relative to the single-locus results. Furthermore, the expanded sample yielded a sharp peak of evidence for linkage disequilibrium that spanned approximately 170 kb and that was centered 100 kb telomeric to HLA-C. The 1.2-Mb interval was further dissected by means of recombinant ancestral haplotype analysis. This analysis identified risk haplotype 1 (RH1), which is a 60-kb fragment of ancestral haplotype 57.1, on all identifiable HLA risk haplotypes. One of these haplotypes exhibits significant linkage disequilibrium with psoriasis but does not carry Cw6, which is the HLA allele most strongly associated with the disease. These results demonstrate that RH1 is highly likely to carry the disease allele at PSORS1, and they exclude HLA-C and corneodesmosin with a high degree of confidence.     

Detecting Local Haplotype Sharing and Haplotype Association

A novel haplotype association method is presented, and its power is demonstrated. Relying on a statistical model for linkage disequilibrium (LD), the method first infers ancestral haplotypes and their loadings at each marker for each individual. The loadings are then used to quantify local haplotype sharing between individuals at each marker. A statistical model was developed to link the local haplotype sharing and phenotypes to test for association. We devised a novel method to fit the LD model, reducing the complexity from putatively quadratic to linear (in the number of ancestral haplotypes). Therefore, the LD model can be fitted to all study samples simultaneously, and, consequently, our method is applicable to big data sets. Compared to existing haplotype association methods, our method integrated out phase uncertainty, avoided arbitrariness in specifying haplotypes, and had the same number of tests as the single-SNP analysis. We applied our method to data from the Wellcome Trust Case Control Consortium and discovered eight novel associations between seven gene regions and five disease phenotypes. Among these, GRIK4, which encodes a protein that belongs to the glutamate-gated ionic channel family, is strongly associated with both coronary artery disease and rheumatoid arthritis. A software package implementing methods described in this article is freely available at http://www.haplotype.org.     

Analysis of the MTHFR gene linkage disequilibrium structure and association of polymorphic gene variants with coronary atherosclerosis

Analysis of the genome-specific linkage disequilibrium patterns in certain populations is a highly promising approach to the identification of functional variants that underlie susceptibility to complex diseases. In the present study, the linkage disequilibrium patterns of the methylenetetrahydrofolate reductase gene (MTHFR) were examined in a group of patients with coronary atherosclerosis (coronary artery disease, CAD) and in a control sample from the Russian population. It was demonstrated that in the samples from one population, which were differentiated by the presence or absence of CAD, the MTHFR linkage disequilibrium patterns had similar features. Association of the MTHFR rs7533315 and rs2066462 polymorphisms with CAD was demonstrated. In addition, the evolution of the haplotypes and their role in the formation of CAD in the Russian population was reconstructed. The data on the association between genetic variability in the MTHFR locus and pathogenetically important indices of lipid metabolism were obtained. The high informativeness of the haplotype approach in case-control tests for associations with CAD was demonstrated.     

Tests for linkage and association in nuclear families.

The transmission/disequilibrium test (TDT) originally was introduced to test for linkage between a genetic marker and a disease-susceptibility locus, in the presence of association. Recently, the TDT has been used to test for association in the presence of linkage. The motivation for this is that linkage analysis typically identifies large candidate regions, and further refinement is necessary before a search for the disease gene is begun, on the molecular level. Evidence of association and linkage may indicate which markers in the region are closest to a disease locus. As a test of linkage, transmissions from heterozygous parents to all of their affected children can be included in the TDT; however, the TDT is a valid chi2 test of association only if transmissions to unrelated affected children are used in the analysis. If the sample contains independent nuclear families with multiple affected children, then one procedure that has been used to test for association is to select randomly a single affected child from each sibship and to apply the TDT to those data. As an alternative, we propose two statistics that use data from all of the affected children. The statistics give valid chi2 tests of the null hypothesis of no association or no linkage and generally are more powerful than the TDT with a single, randomly chosen, affected child from each family.     

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.