- Open Access
Power of the 2-locus TDT for testing the interaction of two susceptibility genes
© Kotti et al; licensee BioMed Central Ltd. 2007
- Published: 18 December 2007
We recently proposed a new strategy: 2-locus TDT for detecting two susceptibility genes through their interaction in trio families. We apply our method to two candidate genes, A and C, on the Genetic Analysis Workshop 15 (GAW15) simulated rheumatoid arthritis data and study the power to identify an interactive effect of these genes.
This study was performed with full knowledge of the answers.
- Rheumatoid Arthritis
- Likelihood Ratio Test
- Restricted Model
- Genetic Analysis Workshop
- Rheumatoid Arthritis Susceptibility
Most multifactorial diseases result from complex interactions among multiple genes with modest effects and various environmental factors. Current strategies usually aim to detect the individual effect of each factor.
We have recently proposed a new strategy: the 2-locus TDT (transmission-disequilibrium test), which tests for interaction between two candidate genes . Like the one-locus TDT , it applies to trio families including the patient and his two parents. Our method is used in the context of studying good candidate pathways. Indeed, we assume prior information that the disease under study is probably caused by particular genes recognized as key elements of the metabolic pathway. We apply the 2-locus TDT to Problem 3 simulated Genetic Analysis Workshop 15 (GAW15) rheumatoid arthritis (RA) data and compute its power for detecting the interaction of two susceptibility genes.
The GAW15 Problem 3 consisted of 100 replicates of 1500 nuclear families with two affected sibs. From each replicate, we extracted 1500 trio families: one affected (the index case) and his two parents.
We focus the analyses on two candidate genes in the RA susceptibility: Locus A on chromosome 16 and Locus C on chromosome 6. Locus A has an indirect effect on RA hazard. Its implication on the susceptibility to the disease depends on the HLA-DRB1 locus located on chromosome 6. Locus C is in strong linkage disequilibrium (LD) with HLA-DRB1. Both genes are involved directly or indirectly in the RA outcome through the HLA-DRB1 gene and thus do not act independently. We apply the 2-locus TDT to the SNPs indicated as functional in the two candidate genes in each replicate of 1500 trios.
We denote A1, A2 and C1, C2 the SNP alleles of genes A and C, respectively. For each patient, there are nine possible genotypes that we denote AiAjCkCl, where Ai, Aj and Ck, Cl are the alleles inherited from the two parents at the loci A and C, respectively. Subscripts i, j, k and l can be either 1 or 2. A2 and C2 are considered to be the high risk alleles and A2A2C2C2 are the reference genotype.
Penetrance estimation under the restricted model
The 2-locus TDT is achieved in two steps. In the first step, the parameters (μ, ν, ϕ, and π) of the restricted model and ε ijkl of the general model are estimated by their maximum likelihood estimator (MLE). In the second one, the restricted model is compared to the general one by a maximum likelihood ratio test.
The estimation of these parameters requires specifying the allele frequencies at each locus. This was achieved by the affected family-based controls (AFBAC) method, using the parental alleles untransmitted to the patients from which the family was ascertained . Here, the families were not ascertained through an affected patient but through a patient with an affected sib. Our parameter estimation must also take into account this specific mode of ascertainment and consequently was achieved by using the marker association segregation chi-square (MASC) program .
This ratio follows, under the null hypothesis, a chi-square distribution with four degrees of freedom (df) corresponding to the difference in the number of parameters between the two models.
For each replicate of 1500 trios, the LRT was computed and the significance of the test was assessed in order to evaluate the power of our test in the present situation.
We illustrate the variation of the LRT and p-values on two contrasting replicates: 9 and 27. In the Replicate 27, there is a strong evidence for interaction: LRT = 14.35 and p = 7 × 10-4. In contrast, we did not find interaction in Replicate 9 because the LRT = 2.49 and p = 0.29.
Replicate 9: Matrix of relative penetrances for the general model
Replicate 27: Matrix of relative penetrances for the general model
In this paper, we studied the performance of the 2-locus TDT for detecting the interaction of the genes A and C in the susceptibility to RA simulated data and show that interaction may be detected with a power of 66% at a 5% level.
The 2-locus TDT searches for a departure from multiplicativity in the two-locus penetrance matrix. This departure cannot be observed without an effect of the two genes. Therefore, the 2-locus TDT may also be used for detecting the effect of two candidate genes through their interaction.
More generally, the 2-locus TDT may be particularly interesting in the situation of weak effect of one or both genes and strong interaction. In such a situation, traditional methods will fail because they will stop in the absence of noticeable main effect. In fact, they ignore the possibility that the interactive effect of multilocus functional genetic variants may play a larger role than the single-locus effect in the determinism of one trait .
Our test applies to candidate genes already recognized as key elements of a metabolic pathway thought to be involved in the disease pathogenesis. There are difficulties in good candidate choice and the use of tagSNPs as markers for most genetic variations in the human genome has become a cornerstone of genetic research. There are, however, reasons to think that the expectations they have raised are somewhat inflated and may lead to unwarranted neglect of other strategies. In particular, we believe in particular that identification of biologically relevant pathways followed by candidate gene strategies may be quite effective at identifying new variants and their interaction in disease susceptibility .
This work was supported by a grant from the Fondation pour la Recherche Médicale.
This article has been published as part of BMC Proceedings Volume 1 Supplement 1, 2007: Genetic Analysis Workshop 15: Gene Expression Analysis and Approaches to Detecting Multiple Functional Loci. The full contents of the supplement are available online at http://www.biomedcentral.com/1753-6561/1?issue=S1.
- Kotti S, Bickeboller H, Clerget-Darpoux F: Strategy for detecting susceptibility genes with weak or no marginal effect. Hum Hered. 2007, 63: 85-92. 10.1159/000099180.View ArticlePubMedGoogle Scholar
- Spielman RS, McGinnis RE, Ewens WJ: Transmission test for linkage disequilibrium: the insulin gene region and insulin-dependent diabetes mellitus (IDDM). Am J Hum Genet. 1993, 52: 506-516.PubMed CentralPubMedGoogle Scholar
- Thomson G: Mapping disease genes: family-based association studies. Am J Hum Genet. 1995, 57: 487-498.PubMed CentralPubMedGoogle Scholar
- Clerget-Darpoux F, Babron MC, Prum B, Lathrop GM, Deschamps I, Hors J: A new method to test genetic models in HLA associated diseases: the MASC method. Ann Hum Genet. 1988, 52: 247-258. 10.1111/j.1469-1809.1988.tb01102.x.View ArticlePubMedGoogle Scholar
- Bourgain C, Genin E, Cox N, Clerget-Darpoux F: Are genome-wide association studies all that we need to dissect the genetic component of complex human diseases?. Eur J Hum Genet. 2007, 15: 260-263. 10.1038/sj.ejhg.5201753.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.