Genetic association analysis using sibship data: a multilevel model approach

Family based association study (FBAS) has the advantages of controlling for population stratification and testing for linkage and association simultaneously. We propose a retrospective multilevel model (rMLM) approach to analyze sibship data by using genotypic information as the dependent variable. Simulated data sets were generated using the simulation of linkage and association (SIMLA) program. We compared rMLM to sib transmission/disequilibrium test (S-TDT), sibling disequilibrium test (SDT), conditional logistic regression (CLR) and generalized estimation equations (GEE) on the measures of power, type I error, estimation bias and standard error. The results indicated that rMLM was a valid test of association in the presence of linkage using sibship data. The advantages of rMLM became more evident when the data contained concordant sibships. Compared to GEE, rMLM had less underestimated odds ratio (OR). Our results support the application of rMLM to detect gene-disease associations using sibship data. However, the risk of increasing type I error rate should be cautioned when there is association without linkage between the disease locus and the genotyped marker.     

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.     

Reconstructing parental genotypes when testing for linkage in the presence of association

Various family-based association methods have recently been proposed that allow testing for linkage in the presence of linkage disequilibrium between a marker and a disease even if there is only incomplete parental-genotype information. For some families, it may be possible to reconstruct missing parental genotypes from the genotypes of their offspring. Treating such a reconstructed family as if parental genotypes have been typed, however, can introduce bias. The reconstruction-combined transmission/disequilibrium test (RC-TDT) and its X-chromosomal counterpart, XRC-TDT, employ parental-genotype reconstruction and correct for the biases involved in this reconstruction without relying on population marker allele frequencies. For the two tests, exact P values can be obtained by numerically calculating the convolution of the null distributions corresponding to the families in the sample.     

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.     

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.     

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.     

A Discordant-Sibship Test for Disequilibrium and Linkage: No Need for Parental Data

The sibship disequilibrium test (SDT) is designed to detect both linkage in the presence of association and association in the presence of linkage (linkage disequilibrium). The test does not require parental data but requires discordant sibships with at least one affected and one unaffected sibling. The SDT has many desirable properties: it uses all the siblings in the sibship; it remains valid if there are misclassifications of the affectation status; it does not detect spurious associations due to population stratification; asymptotically it has a chi2 distribution under the null hypothesis; and exact P values can be easily computed for a biallelic marker. We show how to extend the SDT to markers with multiple alleles and how to combine families with parents and data from discordant sibships. We discuss the power of the test by presenting sample-size calculations involving a complex disease model, and we present formulas for the asymptotic relative efficiency (which is approximately the ratio of sample sizes) between SDT and the transmission/disequilibrium test (TDT) for special family structures. For sib pairs, we compare the SDT to a test proposed both by Curtis and, independently, by Spielman and Ewens. We show that, for discordant sib pairs, the SDT has good power for testing linkage disequilibrium relative both to Curtis''s tests and to the TDT using trios comprising an affected sib and its parents. With additional sibs, we show that the SDT can be more powerful than the TDT for testing linkage disequilibrium, especially for disease prevalence >.3.     

Transmission/Disequilibrium Tests Incorporating Unaffected Offspring

We propose a new method for family-based tests of association and linkage called transmission/disequilibrium tests incorporating unaffected offspring (TDTU). This new approach, constructed based on transmission/disequilibrium tests for quantitative traits (QTDT), provides a natural extension of the transmission/disequilibrium test (TDT) to utilize transmission information from heterozygous parents to their unaffected offspring as well as the affected offspring from ascertained nuclear families. TDTU can be used in various study designs and can accommodate all types of independent nuclear families with at least one affected offspring. When the study sample contains only case-parent trios, the TDTU is equivalent to TDT. Informative-transmission disequilibrium test (i-TDT) and generalized disequilibrium test(GDT) are another two methods that can use information of both unaffected offspring and affected offspring. In contract to i-TDT and GDT, the test statistic of TDTU is simpler and more explicit, and can be implemented more easily. Through computer simulations, we demonstrate that power of the TDTU is slightly higher compared to i-TDT and GDT. All the three methods are more powerful than method that uses affected offspring only, suggesting that unaffected siblings also provide information about linkage and association.     

Linkage analysis and association analysis in the presence of linkage using age at onset of COGA alcoholism data

Complex disease mapping usually involves a combination of linkage and association techniques. Linkage analysis can scan the entire genome in a few hundred tests. Association tests may involve an even greater number of tests. However, association tests can localize the susceptibility genes more accurately. Using a recently developed combined linkage and association strategy, we analyzed a subset of the Collaborative Study on the Genetics of Alcoholism (COGA) data for the Genetic Analysis Workshop 14 (GAW14). In this analysis, we first employed linkage analysis based on frailty models that take into account age of onset information to establish which regions along the chromosome are likely to harbor disease susceptibility genes for alcohol dependence. Second, we used an association analysis by exploiting linkage disequilibrium to narrow down the peak regions. We also compare the methods with mean identity-by-descent tests and transmission/disequilibrium tests that do not use age of onset information.     

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.     

Association studies under general disease models

There is great expectation that the levels of association found between genetic markers and disease status will play a role in the location of disease genes. This expectation follows from regarding association as being proportional to linkage disequilibrium and therefore inversely related to recombination value. For disease genes with more than two alleles, the association measure is instead a weighted average of linkage disequilibria, with the weights depending on allele frequencies and genotype susceptibilities at the disease loci. There is no longer a simple relationship, even in expectation, with recombination. We adopt a general framework to examine association mapping methods which helps to clarify the nature of case-control and transmission/disequilibrium-type tests and reveals the relationship between measures of association and coefficients of linkage disequilibrium. In particular, we can show the consequences of additive and nonadditive effects at the trait locus on the behavior of these tests. These concepts have a natural extension to marker haplotypes. The association of two-locus marker haplotypes with disease phenotype depends on a weighted average of three-locus disequilibria (two markers with each disease locus). It is likely that these two-marker analyses will provide additional information in association mapping studies.     

Sibling-based tests of linkage and association for quantitative traits.

The transmission/disequilibrium test (TDT) developed by Spielman et al. can be a powerful family-based test of linkage and, in some cases, a test of association as well as linkage. It has recently been extended in several ways; these include allowance for implementation with quantitative traits, allowance for multiple alleles, and, in the case of dichotomous traits, allowance for testing in the absence of parental data. In this article, these three extensions are combined, and two procedures are developed that offer valid joint tests of linkage and (in the case of certain sibling configurations) association with quantitative traits, with use of data from siblings only, and that can accommodate biallelic or multiallelic loci. The first procedure uses a mixed-effects (i.e., random and fixed effects) analysis of variance in which sibship is the random factor, marker genotype is the fixed factor, and the continuous phenotype is the dependent variable. Covariates can easily be accommodated, and the procedure can be implemented in commonly available statistical software. The second procedure is a permutation-based procedure. Selected power studies are conducted to illustrate the relative power of each test under a variety of circumstances.     

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

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

Impact of non-ignorable missingness on genetic tests of linkage and/or association using case-parent trios

The transmission/disequilibrium test was introduced to test for linkage disequilibrium between a marker and a putative disease locus using case-parent trios. However, parental genotypes may be incomplete in such a study. When parental information is non-randomly missing, due, for example, to death from the disease under study, the impact on type I error and power under dominant and recessive disease models has been reported. In this paper, we examine non-ignorable missingness by assigning missing values to the genotypes of affected parents. We used unrelated case-parent trios in the Genetic Analysis Workshop 14 simulated data for the Danacaa population. Our computer simulations revealed that the type I error of these tests using incomplete trios was not inflated over the nominal level under either recessive or dominant disease models. However, the power of these tests appears to be inflated over the complete information case due to an excess of heterozygous parents in dyads.     

Transmission test for linkage disequilibrium: the insulin gene region and insulin-dependent diabetes mellitus (IDDM).

A population association has consistently been observed between insulin-dependent diabetes mellitus (IDDM) and the "class 1" alleles of the region of tandem-repeat DNA (5'' flanking polymorphism [5''FP]) adjacent to the insulin gene on chromosome 11p. This finding suggests that the insulin gene region contains a gene or genes contributing to IDDM susceptibility. However, several studies that have sought to show linkage with IDDM by testing for cosegregation in affected sib pairs have failed to find evidence for linkage. As means for identifying genes for complex diseases, both the association and the affected-sib-pairs approaches have limitations. It is well known that population association between a disease and a genetic marker can arise as an artifact of population structure, even in the absence of linkage. On the other hand, linkage studies with modest numbers of affected sib pairs may fail to detect linkage, especially if there is linkage heterogeneity. We consider an alternative method to test for linkage with a genetic marker when population association has been found. Using data from families with at least one affected child, we evaluate the transmission of the associated marker allele from a heterozygous parent to an affected offspring. This approach has been used by several investigators, but the statistical properties of the method as a test for linkage have not been investigated. In the present paper we describe the statistical basis for this "transmission test for linkage disequilibrium" (transmission/disequilibrium test [TDT]). We then show the relationship of this test to tests of cosegregation that are based on the proportion of haplotypes or genes identical by descent in affected sibs. The TDT provides strong evidence for linkage between the 5''FP and susceptibility to IDDM. The conclusions from this analysis apply in general to the study of disease associations, where genetic markers are usually closely linked to candidate genes. When a disease is found to be associated with such a marker, the TDT may detect linkage even when haplotype-sharing tests do not.     

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.     

What is the significance of a significant TDT?

This note summarizes the development of the transmission/disequilibrium test (TDT). The initial purpose of the TDT procedure was to test for linkage between a genetic marker and a disease susceptibility locus when an association had been found between the two. An association between disease and marker had sometimes been taken to imply linkage. An association could, however, be due to population stratification even in the absence of linkage. In contrast, the outcome of the TDT is not affected by such stratification. Furthermore, when linkage is not in doubt, the TDT can, in some cases, also provide a test of association between marker and disease. We discuss these various matters in this paper.     

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.     

Linkage of inflammatory bowel disease to human chromosome 6p

Inflammatory bowel disease (IBD) is characterized by a chronic relapsing intestinal inflammation. IBD is subdivided into Crohn disease and ulcerative colitis phenotypes. Given the immunologic dysregulation in IBD, the human-leukocyte-antigen region on chromosome 6p is of significant interest. Previous association and linkage analysis has provided conflicting evidence as to the existence of an IBD-susceptibility locus in this region. Here we report on a two-stage linkage and association analysis of both a basic population of 353 affected sibling pairs (ASPs) and an extension of this population to 428 white ASPs of northern European extraction. Twenty-eight microsatellite markers on chromosome 6 were genotyped. A peak multipoint LOD score of 4.2 was observed, at D6S461, for the IBD phenotype. A transmission/disequilibrium test (TDT) result of P=.006 was detected for D6S426 in the basic population and was confirmed in the extended cohort (P=.004; 97 vs. 56 transmissions). The subphenotypes of Crohn disease, ulcerative colitis, and mixed IBD contributed equally to this linkage, suggesting a general role for the chromosome 6 locus in IBD. Analysis of five single-nucleotide polymorphisms in the TNFA and LTA genes did not reveal evidence for association of these important candidate genes with IBD. In summary, we provide firm linkage evidence for an IBD-susceptibility locus on chromosome 6p and demonstrate that TNFA and LTA are unlikely to be susceptibility loci for IBD.     

Linkage analysis of human leukocyte antigen (HLA) markers in familial psoriasis: strong disequilibrium effects provide evidence for a major determinant in the HLA-B/-C region.

Although psoriasis is strongly associated with certain human leukocyte antigens (HLAs), evidence for linkage to HLA markers has been limited. The objectives of this study were (1) to provide more definitive evidence for linkage of psoriasis to HLA markers in multiplex families; (2) to compare the major HLA risk alleles in these families with those determined by previous case-control studies; and (3) to localize the gene more precisely. By applying the transmission/disequilibrium test (TDT) and parametric linkage analysis, we found evidence for linkage of psoriasis to HLA-C, -B, -DR, and -DQ, with HLA-B and -C yielding the most-significant results. Linkage was detectable by parametric methods only when marker-trait disequilibrium was considered. Case-control association tests and the TDT identified alleles belonging to the EH57.1 ancestral haplotype as the major risk alleles in our sample. Among individuals carrying recombinant ancestral haplotypes involving EH57. 1, the class I markers were retained selectively among affecteds four times more often than among unaffecteds; among the few affected individuals carrying only the class II alleles from the ancestral haplotype, all but one also carried Cw6. These data show that familial and "sporadic" psoriasis share the same risk alleles. They also illustrate that substantial parametric linkage information can be extracted by accounting for linkage disequilibrium. Finally, they strongly suggest that a major susceptibility gene resides near HLA-C.