An extension of the transmission disequilibrium test incorporating imprinting

The recombination rates in meioses of females and males are often different. Some genes that affect development and behavior in mammals are known to be imprinted, and >1% of all mammalian genes are believed to be imprinted. When the gene is imprinted and the recombination fractions are sex specific, the conventional transmission disequilibrium test (TDT) is shown to be still valid for testing for linkage. The power function of the TDT is derived, and the effect of the degree of imprinting on the power of the TDT is investigated. It is learned that imprinting has little effect on the power when the female and male recombination rates are equal. On the basis of case-parents trios, the transmissions from the heterozygous fathers/mothers to their affected children are separated as paternal and maternal, and two TDT-like statistics, TDT(p) and TDT(m), are consequently constructed. It is found that the TDT(p) possesses a higher power than the TDT for maternal imprinting genes, and the TDT(m) is more powerful than the TDT for paternal imprinting genes. On the basis of the parent-of-origin effects test statistic (POET), a novel statistic, TDT incorporating imprinting (TDTI) is proposed to test for linkage in the presence of linkage disequilibrium, which is shown to be more powerful than the TDT when parent-of-origin effects are significant but slightly less powerful than the TDT when parent-of-origin effects are negligible. The validity of the TDT and TDTI is assessed by simulation. The power approximation formulas for the TDT and TDTI are derived and the simulation results show that they are accurate. The simulation study on power comparison shows that the TDTI outperforms the TDT for imprinted genes. The improvement can be substantial in the case of complete paternal/maternal imprinting.     

The future is now - will the real disease gene please stand up?

The transmission/disequilibrium test (TDT) [Spielman et al.: Am J Hum Genet 1993;52:506-516] has been postulated as the future of gene mapping for complex diseases, provided one is able to genotype a dense enough map of markers across the genome. Risch and Merikangas [Science 1996;273:1516-1517] suggested a million-marker screen in affected sibpair (ASP) families, demonstrating that the TDT is a more powerful test of linkage than traditional linkage tests based on allele-sharing when there is also association between marker and disease alleles. While the future of genotyping has arrived, successes in family-based association studies have been modest. This is often attributed to excessive false positives in candidate gene studies. This problem is only exacerbated by the increasing numbers of whole genome association (WGA) screens. When applied in ASPs, the TDT statistic, which assumes transmissions to siblings are independent, is not expected to have a constant variance in the presence of variable linkage. This results in generally more extreme statistics, hence will further aggravate the problem of having a large number of positive results to sort through. So an important question is how many positive TDT results will show up on a chromosome containing a disease gene due only to linkage, and will they obfuscate the true disease gene location. To answer this question we combined theory and computer simulations. These studies show that in ASPs the normal version of the TDT statistic has a mean of 0 and a variance of 1 in unlinked regions, but has a variance larger than 1 in linked regions. In contrast, the pedigree disequilibrium test (PDT) statistic adjusts for correlation between siblings due to linkage and maintains a constant variance of 1 at unassociated markers irrespective of linkage. The TDT statistic is generally larger than the PDT statistic across linked regions. This is true for unassociated as well as associated markers. To compare the two tests we ranked both statistics at the disease locus, or an associated marker, among statistics at all other markers. The TDT did better job than PDT placing the score of the associated marker near the top. Though, strictly speaking, the TDT in ASPs should be interpreted as a test of linkage and not a test of association, there is a good chance that if a marker stands out, the marker is associated as well as linked. In conclusion, our results suggest that TDT is an effective screening tool for WGA studies, especially in multiplex families.     

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.     

How do homozygous parents affect TDT as a test for association?

The traditional transmission disequilibrium test (TDT) (Spielman et al., 1993) is a powerful test for association only in the presence of linkage. Since allele transmissions from homozygous parents do not carry any information on linkage, the TDT statistic uses data only on heterozygous parents. However, homozygous parents carry information on association between alleles at a marker locus and a disease locus. In this article, we explore whether inclusion of homozygous parents increases the power to detect association. The resultant test statistic follows a chi(2) distribution with 2 degrees of freedom. Monte-Carlo simulations are included to compare the performance of this test with the traditional TDT under different disease models.     

A test of transmission/disequilibrium for quantitative traits in pedigree data, by multiple regression.

The transmission/disequilibrium (TD) test (TDT), proposed, by Spielman et al., for binary traits is a powerful method for detection of linkage between a marker locus and a disease locus, in the presence of allelic association. As a test for linkage disequilibrium, the TDT makes the assumption that any allelic association present is due to linkage. Allison proposed a series of TD-type tests for quantitative traits and calculated their power, assuming that the marker locus is the disease locus. All these tests assume that the observations are independent, and therefore they are applicable, as a test for linkage, only for nuclear-family data. In this report, we propose a regression-based TD-type test for linkage between a marker locus and a quantitative trait locus, using information on the parent-to-offspring transmission status of the associated allele at the marker locus. This method does not require independence of observations, thus allowing for analysis of pedigree data as well, and allows adjustment for covariates. We investigate the statistical power and validity of the test by simulating markers at various recombination fractions from the disease locus.     

Transmission disequilibrium test with discordant sib pairs when parents are available

The transmission disequilibrium test (TDT) has been employed to map disease susceptibility loci (DSL), while being immune to the problem of population admixture. The customary TDT test (TDT(D)) was developed for affected child(ren) and their parents and was most often applied to case-parent trios. Recently, the TDT has been extended to the situations when (1) parents are not available but affected and nonaffected sibs from each family are available, (2) unrelated control-parent trios are available for combined analyses with case-parent trios (TDT(DC)), and (3) large pedigrees. For many diseases, affected children in the case-parent trios enlisted into the TDT(D) have unaffected sibs who can be recruited. We present an extension of the TDT by effectively incorporating one unaffected sib of each of the affected children in the case-parent trios into a single analysis (TDT(DS), where DS denotes discordant sib pairs). We have developed a general analytical method for computing the statistical power of the TDT(DS) under any genetic model, the accuracy of which is validated by computer simulations. We compare the power of the TDT(D), TDT(DC), and TDT(DS) under a range of parameter space and genetic models. We find that the TDT(DS) is generally more powerful than the TDT(DC) and TDT(D), particularly when the disease is prevalent (>30%) in the population. The relative power of the TDT(D) and the TDT(DS) largely depends upon the allele frequencies and genetic effects at the DSL, whereas the recombination rate, the degree of linkage disequilibrium, and the marker allele frequencies have little effect. Importantly, the TDT(DS) not only may be more powerful, it also has the advantage of being able to test for segregation distortion that may yield false linkage/association in the TDT(D).     

X-APL: an improved family-based test of association in the presence of linkage for the X chromosome

Family-based association methods have been developed primarily for autosomal markers. The X-linked sibling transmission/disequilibrium test (XS-TDT) and the reconstruction-combined TDT for X-chromosome markers (XRC-TDT) are the first association-based methods for testing markers on the X chromosome in family data sets. These are valid tests of association in family triads or discordant sib pairs but are not theoretically valid in multiplex families when linkage is present. Recently, XPDT and XMCPDT, modified versions of the pedigree disequilibrium test (PDT), were proposed. Like the PDT, XPDT compares genotype transmissions from parents to affected offspring or genotypes of discordant siblings; however, the XPDT can have low power if there are many missing parental genotypes. XMCPDT uses a Monte Carlo sampling approach to infer missing parental genotypes on the basis of true or estimated population allele frequencies. Although the XMCPDT was shown to be more powerful than the XPDT, variability in the statistic due to the use of an estimate of allele frequency is not properly accounted for. Here, we present a novel family-based test of association, X-APL, a modification of the test for association in the presence of linkage (APL) test. Like the APL, X-APL can use singleton or multiplex families and properly infers missing parental genotypes in linkage regions by considering identity-by-descent parameters for affected siblings. Sampling variability of parameter estimates is accounted for through a bootstrap procedure. X-APL can test individual marker loci or X-chromosome haplotypes. To allow for different penetrances in males and females, separate sex-specific tests are provided. Using simulated data, we demonstrated validity and showed that the X-APL is more powerful than alternative tests. To show its utility and to discuss interpretation in real-data analysis, we also applied the X-APL to candidate-gene data in a sample of families with Parkinson disease.     

Linkage analysis of alcohol dependence using both affected and discordant sib pairs

The basic idea of affected-sib-pair (ASP) linkage analysis is to test whether the inheritance pattern of a marker deviates from Mendelian expectation in a sample of ASPs. The test depends on an assumed Mendelian control distribution of the number of marker alleles shared identical by descent (IBD), i.e., 1/4, 1/2, and 1/4 for 2, 1, and 0 allele(s) IBD, respectively. However, Mendelian transmission may not always hold, for example because of inbreeding or meiotic drive at the marker or a nearby locus. A more robust and valid approach is to incorporate discordant-sib-pairs (DSPs) as controls to avoid possible false-positive results. To be robust to deviation from Mendelian transmission, here we analyzed Collaborative Study on the Genetics of Alcoholism data by modifying the ASP LOD score method to contrast the estimated distribution of the number of allele(s) shared IBD by ASPs with that by DSPs, instead of with the expected distribution under the Mendelian assumption. This strategy assesses the difference in IBD sharing between ASPs and the IBD sharing between DSPs. Further, it works better than the conventional LOD score ASP linkage method in these data in the sense of avoiding false-positive linkage evidence.     

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.     

The affected-/discordant-sib-pair design can guarantee validity of multipoint model-free linkage analysis of incomplete pedigrees when there is marker-marker disequilibrium

Genomewide linkage studies are tending toward the use of single-nucleotide polymorphisms (SNPs) as the markers of choice. However, linkage disequilibrium (LD) between tightly linked SNPs violates the fundamental assumption of linkage equilibrium (LE) between markers that underlies most multipoint calculation algorithms currently available, and this leads to inflated affected-relative-pair allele-sharing statistics when founders' multilocus genotypes are unknown. In this study, we investigate the impact that the degree of LD, marker allele frequency, and association type have on estimating the probabilities of sharing alleles identical by descent in multipoint calculations and hence on type I error rates of different sib-pair linkage approaches that assume LE. We show that marker-marker LD does not inflate type I error rates of affected sib pair (ASP) statistics in the whole parameter space, and that, in any case, discordant sib pairs (DSPs) can be used to control for marker-marker LD in ASPs. We advocate the ASP/DSP design with appropriate sib-pair statistics that test the difference in allele sharing between ASPs and DSPs.     

Marker selection for the transmission/disequilibrium test, in recently admixed populations.

Recent admixture between genetically differentiated populations can result in high levels of association between alleles at loci that are <=10 cM apart. The transmission/disequilibrium test (TDT) proposed by Spielman et al. (1993) can be a powerful test of linkage between disease and marker loci in the presence of association and therefore could be a useful test of linkage in admixed populations. The degree of association between alleles at two loci depends on the differences in allele frequencies, at the two loci, in the founding populations; therefore, the choice of marker is important. For a multiallelic marker, one strategy that may improve the power of the TDT is to group marker alleles within a locus, on the basis of information about the founding populations and the admixed population, thereby collapsing the marker into one with fewer alleles. We have examined the consequences of collapsing a microsatellite into a two-allele marker, when two founding populations are assumed for the admixed population, and have found that if there is random mating in the admixed population, then typically there is a collapsing for which the power of the TDT is greater than that for the original microsatellite marker. A method is presented for finding the optimal collapsing that has minimal dependence on the disease and that uses estimates either of marker allele frequencies in the two founding populations or of marker allele frequencies in the current, admixed population and in one of the founding populations. Furthermore, this optimal collapsing is not always the collapsing with the largest difference in allele frequencies in the founding populations. To demonstrate this strategy, we considered a recent data set, published previously, that provides frequency estimates for 30 microsatellites in 13 populations.     

Unbiased application of the transmission/disequilibrium test to multilocus haplotypes

When the transmission/disequilibrium test (TDT) is applied to multilocus haplotypes, a bias may be introduced in some families for which both parents have the same heterozygous genotype at some locus. The bias occurs because haplotypes can only be deduced from certain offspring, with the result that the transmissions of the two parental haplotypes are not independent. We obtain an unbiased TDT for individual haplotypes by calculating the correct variance for the transmission count within a family, using information from multiple siblings if they are available. An existing correction for dependence between siblings in the presence of linkage is retained. To obtain an unbiased multihaplotype TDT, we must either count transmissions from one randomly chosen parent or count all transmissions and estimate the significance level empirically. Alternatively, we may use missing-data techniques to estimate uncertain haplotypes, but these methods are not robust to population stratification. An illustration using data from the insulin-gene region in type 1 diabetes shows that the validity and power of the TDT may vary by an order of magnitude, depending on the method of analysis.     

The transmission/disequilibrium test: history, subdivision, and admixture.

Disease association with a genetic marker is often taken as a preliminary indication of linkage with disease susceptibility. However, population subdivision and admixture may lead to disease association even in the absence of linkage. In a previous paper, we described a test for linkage (and linkage disequilibrium) between a genetic marker and disease susceptibility; linkage is detected by this test only if association is also present. This transmission/disequilibrium test (TDT) is carried out with data on transmission of marker alleles from parents heterozygous for the marker to affected offspring. The TDT is a valid test for linkage and association, even when the association is caused by population subdivision and admixture. In the previous paper, we did not explicitly consider the effect of recent history on population structure. Here we extend the previous results by examining in detail the effects of subdivision and admixture, viewed as processes in population history. We describe two models for these processes. For both models, we analyze the properties of (a) the TDT as a test for linkage (and association) between marker and disease and (b) the conventional contingency statistic used with family data to test for population association. We show that the contingency test statistic does not have a chi 2 distribution if subdivision or admixture is present. In contrast, the TDT remains a valid chi 2 statistic for the linkage hypothesis, regardless of population history.     

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 transmission/disequilibrium test approach to screen for quantitative trait loci in two selected lines of large white pigs

Pedigree and marker data from a multiple-generation pig selection experiment have been analysed to screen for loci affecting quantitative traits (QTL). Pigs from a base population were selected either for low backfat thickness at fixed live weight (L-line) or high live weight at fixed age (F-line). Selection was based on single-trait own performance and DNA was available on selected individuals only. Genotypes for three marker loci with known positions on chromosome 4 were available. The transmission/disequilibrium test (TDT) was originally described in human genetics to test for linkage between a genetic marker and a disease-susceptibility locus, in the presence of association. Here, we adapt the TDT to test for linkage between a marker and QTL favoured by selection, and for linkage disequilibrium between them in the base population. The a priori unknown distribution of the test statistic under the null hypothesis, no linkage, was obtained via Monte Carlo simulation. Significant TDT statistics were found for markers AFABP and SW818 in the F-line, indicating the presence of a closely linked QTL affecting growth performance. In the L-line, none of the markers studied showed significance. This study emphasizes the potential of the TDT as a quick and simple approach to screen for QTL in situations where marker genotypes are available on selected individuals. The results suggest that previously identified QTL in crosses of genetically diverse breeds may also segregate in commercial selection lines.     

Assessing whether an allele can account in part for a linkage signal: the Genotype-IBD Sharing Test (GIST)

To fine map genes, investigators often test for disease-marker association in chromosomal regions with evidence for linkage. Given a marker allele tentatively associated with disease, one would ask if this allele, or one in linkage disequilibrium (LD) with it, could account in part for the observed linkage signal. This question can be addressed by determining if families selected on the basis of the presence of the tentatively associated allele show stronger evidence of linkage as measured by increased allele sharing identical by descent (IBD) by affected family members. However, common selection strategies can be biased for or against linkage in the marker region, even given no disease-marker association. We define unbiased selection schemes and extend the definition to allow weighted selection on the basis of all genotyped family members. For affected-sibship data, we describe three genotype-based weight variables, corresponding to dominant, recessive, and additive models. We then introduce a test for association of a family weight variable with excess IBD sharing. This test allows us to determine if the linkage signal in a region can be attributed in part to the presence of a marker allele, either because of direct involvement in disease etiology or because of LD with a predisposing genetic variant. For samples of 500 affected sib pairs, the tests are powerful in detection of genotype-IBD sharing association, even for disease models with sib relative risk as low as lambda S=1.1, or when evidence for linkage is absent because of sampling variation. This makes our method a new tool for detecting linkage as well as association, especially in regions harboring a candidate gene. We have implemented these methods in the software package GIST (Genotype-IBD Sharing Test).     

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.     

Inferring linkage disequilibrium between a polymorphic marker locus and a trait locus in natural populations

Three approaches are proposed in this study for detecting or estimating linkage disequilibrium between a polymorphic marker locus and a locus affecting quantitative genetic variation using the sample from random mating populations. It is shown that the disequilibrium over a wide range of circumstances may be detected with a power of 80% by using phenotypic records and marker genotypes of a few hundred individuals. Comparison of ANOVA and regression methods in this article to the transmission disequilibrium test (TDT) shows that, given the genetic variance explained by the trait locus, the power of TDT depends on the trait allele frequency, whereas the power of ANOVA and regression analyses is relatively independent from the allelic frequency. The TDT method is more powerful when the trait allele frequency is low, but much less powerful when it is high. The likelihood analysis provides reliable estimation of the model parameters when the QTL variance is at least 10% of the phenotypic variance and the sample size of a few hundred is used. Potential use of these estimates in mapping the trait locus is also discussed.     

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.     

Application and interpretation of transmission/disequilibrium tests: transmission of HLA-DQ haplotypes to unaffected siblings in 526 families with type 1 diabetes

It is widely believed that, if a genetic marker shows a transmission distortion in patients by the transmission/disequilibrium test (TDT), then a transmission distortion in healthy siblings would be seen in the opposite direction. This is also the case in a complex disease. Furthermore, it has been suggested that replacing the McNemar statistics of the TDT with a test of heterogeneity between transmissions to affected and unaffected children could increase the power to detect disease association. To test these two hypotheses empirically, we analyzed the transmission of HLA-DQA1-DQB1 haplotypes in 526 Norwegian families with type 1 diabetic children and healthy siblings, since some DQA1-DQB1 haplotypes represent major genetic risk factors for type 1 diabetes. Despite the strong positive and negative disease associations with particular DQ haplotypes, we observed no significant deviation from 50% for transmission to healthy siblings. This could be explained by the low penetrance of susceptibility alleles, together with the fact that IDDM loci also harbor strongly protective alleles that can override the risk contributed by other loci. Our results suggest that, in genetically complex diseases, detectable distortion in transmission to healthy siblings should not be expected. Furthermore, the original TDT seems more powerful than a heterogeneity test.