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).     

Power calculations for the transmission/disequilibrium and affected sib pair tests using elementary probability methods

The transmission/disequilibrium test (TDT) and the affected sib pair test (ASP) both test for the association of a marker allele with some conditions. Here, we present methods for calculating the probability of detecting the association (power) for a study examining a fixed number of families for suitability for the study and for calculating the number of such families to be examined. Both calculations use a genetic model for the association. The model considered posits a bi-allelic marker locus that is linked to a bi-allelic disease locus with a possibly nonzero recombination fraction between the loci. The penetrance of the disease is an increasing function of the number of disease alleles. The TDT tests whether the transmission by a heterozygous parent of a particular allele at a marker locus to an affected offspring occurs with probability greater than 0.5. The ASP tests whether transmission of the same allele to two affected sibs occurs with probability greater than 0.5. In either case, evidence that the probability is greater than 0.5 is evidence for association between the marker and the disease. Study inclusion criteria (IC) can greatly affect the necessary sample size of a TDT or ASP study. IC considered by us include a randomly selected parent at least one parent or both parents required to be heterozygous. It also allows a specified minimum number of affected offspring to be required (TDT only). We use elementary probability calculations rather than complex mathematical manipulations or asymptotic methods (large sample size approximations) to compute power and requisite sample size for a proposed study. The advantages of these methods are simplicity and generality.     

The power of the transmission disequilibrium test (TDT) with both case-parent and control-parent trios

The transmission disequilibrium test (TDT) customarily uses affected children and their parents (often case-parent trios, TDTD). Control-parent trios are necessary to guard against spurious significant results due to segregation distortion but are not generally utilized in the identification of disease susceptibility loci (DSL). Controls are often easy to recruit and the TDT can easily be extended to include control-parent trios into the analyses with unrelated case-parent trios. We present an extension of the TDT (TDTDC) that incorporates unrelated cases and controls and their parents into a single analysis. We develop a simple and accurate analytical method for computing the statistical power of various TDT (e.g. the TDTD, TDTDC, TDTDC and TDTC that employ control-parent trios only) under any genetic model. We investigated the power of these TDT, and particularly compared the relative power of the TDTD and TDTDC. We found that the TDTDC is almost always more powerful than the TDTC and TDTD. The relative power of the TDTDC and TDTD depends largely upon a number of parameters identified in the study. This study provides a basis for efficient use of control-parent trios in DSL identification.     

Allowing for missing parents in genetic studies of case-parent triads

In earlier work, my colleagues and I described a log-linear model for genetic data from triads composed of affected probands and their parents. This model allows detection of and discrimination between effects of an inherited haplotype versus effects of the maternal haplotype, which presumably would be mediated by prenatal factors. Like the transmission disequilibrium test (TDT), the likelihood-ratio test (LRT) based on this model is not sensitive to associations that are due to genetic admixture. When used as a method for testing for linkage disequilibrium, the LRT can be regarded as an alternative to the TDT. When one or both parents are missing, the resulting incomplete triad must be discarded to ensure validity of the TDT, thereby sacrificing information. By contrast, when the problem is set in a likelihood framework, the expectation-maximization algorithm allows the incomplete triads to contribute their information to the LRT without invalidation of the analysis. Simulations demonstrate that much of the lost statistical power can be recaptured by means of this missing-data technique. In fact, power is reasonably good even when no triad is complete-for example, when a study is designed to include only mothers of cases. Information from siblings also can be incorporated to further improve the statistical power when genetic data from parents or probands are missing.     

Gene-environment interactions and the affected-sib-pair designs

In genetic mapping of complex traits, the affected-sib-pair method (ASP) and the transmission disequilibrium test (TDT) are two methods of choice. The major appeal of both ASP and TDT is that they do not require the knowledge of mode of inheritance underlying the trait in question. The relative ease and economy for data collection also is the reason for their popularity. The basic idea of the ASP is to identify genes or chromosomal regions through identifying genetic similarity based on phenotypical similarity. TDT, on the other hand, detects susceptibility genes through detecting unusual transmission patterns in families. Since phenotypic similarity can also be caused by environmental similarity, we investigate how the presence of gene-environment interaction (GEI) affects the power of both methods. For a simple one-locus-one-risk-factor model, our results indicate that, in the presence of GEI, methods developed based on marginal penetrance functions (i.e. ignoring the risk factor) can give spurious results. The triangular restriction on allele-sharing probability may no longer be valid. If the environment effect is strong, using exposure-discordant affected sib pairs may have advantage over other designs. Above all, a genetic model involving both genetic and environmental factors behaves differently from a single-locus genetic model.     

The power of genomic control

Although association analysis is a useful tool for uncovering the genetic underpinnings of complex traits, its utility is diminished by population substructure, which can produce spurious association between phenotype and genotype within population-based samples. Because family-based designs are robust against substructure, they have risen to the fore of association analysis. Yet, if population substructure could be ignored, this robustness can come at the price of power. Unfortunately it is rarely evident when population substructure can be ignored. Devlin and Roeder recently have proposed a method, termed "genomic control" (GC), which has the robustness of family-based designs even though it uses population-based data. GC uses the genome itself to determine appropriate corrections for population-based association tests. Using the GC method, we contrast the power of two study designs, family trios (i.e., father, mother, and affected progeny) versus case-control. For analysis of trios, we use the TDT test. When population substructure is absent, we find GC is always more powerful than TDT; furthermore, contrary to previous results, we show that as a disease becomes more prevalent the discrepancy in power becomes more extreme. When population substructure is present, however, the results are more complex: TDT is more powerful when population substructure is substantial, and GC is more powerful otherwise. We also explore general issues of power and implementation of GC within the case-control setting and find that, economically, GC is at least comparable to and often less expensive than family-based methods. Therefore, GC methods should prove a useful complement to family-based methods for the genetic analysis of complex traits.     

Screening and replication using the same data set: testing strategies for family-based studies in which all probands are affected

For genome-wide association studies in family-based designs, we propose a powerful two-stage testing strategy that can be applied in situations in which parent-offspring trio data are available and all offspring are affected with the trait or disease under study. In the first step of the testing strategy, we construct estimators of genetic effect size in the completely ascertained sample of affected offspring and their parents that are statistically independent of the family-based association/transmission disequilibrium tests (FBATs/TDTs) that are calculated in the second step of the testing strategy. For each marker, the genetic effect is estimated (without requiring an estimate of the SNP allele frequency) and the conditional power of the corresponding FBAT/TDT is computed. Based on the power estimates, a weighted Bonferroni procedure assigns an individually adjusted significance level to each SNP. In the second stage, the SNPs are tested with the FBAT/TDT statistic at the individually adjusted significance levels. Using simulation studies for scenarios with up to 1,000,000 SNPs, varying allele frequencies and genetic effect sizes, the power of the strategy is compared with standard methodology (e.g., FBATs/TDTs with Bonferroni correction). In all considered situations, the proposed testing strategy demonstrates substantial power increases over the standard approach, even when the true genetic model is unknown and must be selected based on the conditional power estimates. The practical relevance of our methodology is illustrated by an application to a genome-wide association study for childhood asthma, in which we detect two markers meeting genome-wide significance that would not have been detected using standard methodology.     

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.     

The score statistic of the LD-lod analysis: detecting linkage adaptive to linkage disequilibrium

We study the properties of a modified lod score method for testing linkage that incorporates linkage disequilibrium (LD-lod). By examination of its score statistic, we show that the LD-lod score method adaptively combines two sources of information: (a) the IBD sharing score which is informative for linkage regardless of the existence of LD and (b) the contrast between allele-specific IBD sharing scores which is informative for linkage only in the presence of LD. We also consider the connection between the LD-lod score method and the transmission-disequilibrium test (TDT) for triad data and the mean test for affected sib pair (ASP) data. We show that, for triad data, the recessive LD-lod test is asymptotically equivalent to the TDT; and for ASP data, it is an adaptive combination of the TDT and the ASP mean test. We demonstrate that the LD-lod score method has relatively good statistical efficiency in comparison with the ASP mean test and the TDT for a broad range of LD and the genetic models considered in this report. Therefore, the LD-lod score method is an interesting approach for detecting linkage when the extent of LD is unknown, such as in a genome-wide screen with a dense set of genetic markers.     

Association mapping in structured populations

The use, in association studies, of the forthcoming dense genomewide collection of single-nucleotide polymorphisms (SNPs) has been heralded as a potential breakthrough in the study of the genetic basis of common complex disorders. A serious problem with association mapping is that population structure can lead to spurious associations between a candidate marker and a phenotype. One common solution has been to abandon case-control studies in favor of family-based tests of association, such as the transmission/disequilibrium test (TDT), but this comes at a considerable cost in the need to collect DNA from close relatives of affected individuals. In this article we describe a novel, statistically valid, method for case-control association studies in structured populations. Our method uses a set of unlinked genetic markers to infer details of population structure, and to estimate the ancestry of sampled individuals, before using this information to test for associations within subpopulations. It provides power comparable with the TDT in many settings and may substantially outperform it if there are conflicting associations in different subpopulations.     

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.     

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.     

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.     

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.     

Combining evidence for association from transmission disequilibrium and case-control studies using single-nucleotide polymorphisms

The aim of the present analysis is to combine evidence for association from the two most commonly used designs in genetic association analysis, the case-control design and the transmission disequilibrium test (TDT) design. The cases here are affected offspring from nuclear families and are used in both the case-control and TDT designs. As a result, inference from these designs is not independent. We applied a simple logistic regression method for combining evidence for association from case-control and TDT designs to single-nucleotide polymorphism data purchased on a region on chromosome 3, replicate 1 of the Aipotu population. Combining the evidence from the case-control and TDT designs yielded a 5-10% reduction in the standard errors of the relative risk estimates. The authors did not know the results before the analyses were conducted.     

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.     

Handling Missing Data in Transmission Disequilibrium Test in Nuclear Families with One Affected Offspring

The Transmission Disequilibrium Test (TDT) compares frequencies of transmission of two alleles from heterozygote parents to an affected offspring. This test requires all genotypes to be known from all members of the nuclear families. However, obtaining all genotypes in a study might not be possible for some families, in which case, a data set results in missing genotypes. There are many techniques of handling missing genotypes in parents but only a few in offspring. The robust TDT (rTDT) is one of the methods that handles missing genotypes for all members of nuclear families [with one affected offspring]. Even though all family members can be imputed, the rTDT is a conservative test with low power. We propose a new method, Mendelian Inheritance TDT (MITDT-ONE), that controls type I error and has high power. The MITDT-ONE uses Mendelian Inheritance properties, and takes population frequencies of the disease allele and marker allele into account in the rTDT method. One of the advantages of using the MITDT-ONE is that the MITDT-ONE can identify additional significant genes that are not found by the rTDT. We demonstrate the performances of both tests along with Sib-TDT (S-TDT) in Monte Carlo simulation studies. Moreover, we apply our method to the type 1 diabetes data from the Warren families in the United Kingdom to identify significant genes that are related to type 1 diabetes.     

A log-linear approach to case-parent-triad data: assessing effects of disease genes that act either directly or through maternal effects and that may be subject to parental imprinting.

We describe a log-linear method for analysis of case-parent-triad data, based on maximum likelihood with stratification on parental mating type. The method leads to estimates of association parameters, such as relative risks, for a single allele, and also to likelihood ratio chi2 tests (LRTs) of linkage disequilibrium. Hardy-Weinberg equilibrium need not be assumed. Our simulations suggest that the LRT has power similar to that of the chi2 "score" test proposed by Schaid and Sommer and that both can outperform the transmission/disequilibrium test (TDT), although the TDT can perform better under an additive model of inheritance. Because a restricted version of the LRT is asymptotically equivalent to the TDT, the proposed test can be regarded as a generalization of the TDT. The method that we describe generalizes easily to accommodate maternal effects on risk and, in fact, produces powerful and orthogonal tests of the contribution of fetal versus maternal genetic factors. We further generalize the model to allow for effects of parental imprinting. Imprinting effects can be fitted by a simple, iterative procedure that relies on the expectation-maximization algorithm and that uses standard statistical software for the maximization steps. Simulations reveal that LRT tests for detection of imprinting have very good operating characteristics. When a single allele is under study, the proposed method can yield powerful tests for detection of linkage disequilibrium and is applicable to a broader array of causal scenarios than is the TDT.     

Use of susceptibility scoring in conjunction with the genotypic transmission disequilibrium test

We explored the utility of selecting a genetically predisposed subgroup to increase the finding of a genetic signal in the Genetic Analysis Workshop 14 Collaborative Study on the Genetics of Alcoholism dataset. A subgroup of affected probands with low environmental risk exposures was defined using a susceptibility score calculated from an environmental risk model. Thirty-nine probands with highly positive scores were selected, along with their parents, for use in a genotypic transmission disequilibrium test (TDT) test. We compared the results of the genotypic TDT in this subgroup to the TDT results using all probands and their parents. For some markers, the susceptibility scoring approach resulted in smaller p-values, while for other markers, evidence for a genetic signal weakened. Further explorations into genetic and environmental population characteristics that benefit from this approach are warranted.     

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.