- Methodology article
- Open Access
Combining identity by descent and association in genetic case-control studies
- Qingrun Zhang^{1},
- Shuang Wang^{2} and
- Jurg Ott^{1, 3}Email author
https://doi.org/10.1186/1471-2156-9-42
© Zhang et al; licensee BioMed Central Ltd. 2008
- Received: 25 December 2007
- Accepted: 05 July 2008
- Published: 05 July 2008
Abstract
Background
In human case-control association studies, one of the chi-square tests typically carried out is based on a 2 × 3 table of genotypes (homogeneity of three genotype frequencies in case and control individuals). We formulate the two degrees of freedom associated with a given genotype distribution in terms of two biologically relevant parameters, (1) the probability F that an individual's two alleles are identical by descent (IBD) and (2) the frequency p of one of the alleles.
Results
Imposing the restriction, F ≥ 0, makes some of the genotype frequencies invalid thereby reducing noise. We propose a new statistical association test, the FP test, by focusing on allele frequency differences between case and control individuals while allowing for suitable IBD probabilities. Power calculations show that (1) the practice of generally carrying out two association tests (allele and genotype test) has an increased type I error and (2) our test is more powerful than conventional genotype and allele tests under recessive trait inheritance, and at least as powerful as these conventional tests under dominant inheritance.
Conclusion
For dominant and recessive modes of inheritance, any apparent power gain by an allele test when carried out in conjunction with a genotype test tends to be purchased entirely by an increased rate of false positive results due to omission of a multiple testing correction. As an alternative to these two standard association tests, our FP test represents a convenient and more powerful alternative.
Keywords
- Control Individual
- Association Test
- Recessive Trait
- Genotype Test
- Allele Test
Background
In their well-known paper on homozygosity mapping published 20 years ago, Lander and Botstein [1] showed that fewer than a dozen unrelated inbred children should suffice to map a recessive trait given a dense map of genetic markers. They recommended that one should search for extended regions of homozygosity shared by a set of inbred individuals. The power of this approach was demonstrated by the mapping of a gene for a rare recessive trait in a genome-wide investigation of only three distantly related patients [2]: On chromosome 18, five of the six disease chromosomes shared a region of 19 cM in length.
Even individuals seemingly collected at random from the population tend to exhibit extended regions of allele sharing [3], which may be interpreted as the consequence of the mating of (distantly) related individuals. Such homozygosity represents autozygosity, that is, the sharing of two alleles that are copies of a single ancestral allele [3].
In recent years, researchers have shown renewed interest in extended segments of homozygosity and have generally done so by focusing on segments of specific lengths [4–7]. Our purpose here is to focus on individual SNPs rather than on genomic segments of arbitrary lengths, and to allow for biologically reasonable values of IBD while testing for allele frequency differences. The reason for this is not so much biological plausibility but rather statistical power: By disregarding parameter values that are unlikely to be of importance, we reduce the effect of statistical noise and thereby gain accuracy and power. It will be seen below that our approach is somewhat analogous to the "possible triangle" method [8] in affected sib-pair linkage analysis and is expected to lead to similar increases in power.
Results
Statistical model for SNP genotype frequencies
Genotype parametrization
Genotype | Frequency |
---|---|
AA | q_{1} = Fp + (1 - F)p^{2} |
AB | q_{2} = 2p(1 - p)(1 - F) |
BB | q_{3} = F(1 - p) + (1 - F)(1 - p)^{2} |
Sum | 1 |
The MLEs of q_{1}, q_{2}, and q_{3} are simply the proportions of individuals with given genotypes. Because of the invariance property of MLEs, the functions p and F of q_{1} and q_{2} (equation 1) are also MLEs. The inbreeding coefficient F may be viewed as an indicator of how far the genotype frequencies deviate from HWE. If q_{2} is smaller than expected under HWE then F > 0.
The inbreeding coefficient F in human populations is known to be small and positive. In North America, it is generally much smaller than 0.01 [9] but in isolated populations may reach value of around 0.10 [10]. However, in samples of individuals affected with a heritable trait, the inbreeding coefficient may be even higher. For example, for the SNP most strongly associated with age-related macular degeneration (AMD) [11], rs380390, application of equations (1) leads to an estimated inbreeding coefficient of 0.13 in case individuals. Such increases occur because of enrichment of a disease genotype due to ascertainment of case individuals, and/or they may be a consequence of the fact that individuals affected with a heritable trait may be distantly related. Conversely, the estimated (unrestricted) F value for rs380390 was -0.07 in control individuals.
So far, the expressions for the genotype frequencies in Table 1 simply amount to a specific formulation of the two df's associated with the three frequencies. Several other such parametrizations have previously been proposed [12]. The current transformation (1) of genotype frequencies does not by itself yield any new insights. However, in the next section we impose restrictions on the range of parameter values, which will result in a new test.
Statistical test
We want to test the null hypothesis H_{0} of no association versus the alternative hypothesis H_{1} of association. Under H_{0}, allele frequencies and F values are the same in case and control individuals while under H_{1}, allele frequencies may be different between cases and controls and so may be F values. Thus, this parametrization allows testing for allele frequency differences but does so by working with genotype frequencies. To make our test as powerful as possible, we restrict estimates of inbreeding coefficients to non-negative values as these are expected to occur preferentially under H_{1}. As outlined in detail in the Methods section, we formulate this test as a likelihood ratio (LR) test, the FP test. The log likelihood for case individuals is given by $\mathrm{log}\phantom{\rule{0.1em}{0ex}}[{L}_{a}({p}_{a},{F}_{a})]={\displaystyle {\sum}_{i=1}^{3}{n}_{i}\mathrm{log}\phantom{\rule{0.1em}{0ex}}({q}_{i})}$ where n_{1}, n_{2}, and n_{3} are the respective numbers of case individuals with genotypes AA, AB, and BB, the q_{i} are functions of F and p (Table 1), p_{ a }and F_{ a }are the parameter values in case individuals. For control individuals, the log likelihood log [L_{ b }(p_{ b }, F_{ b })] is obtained in an analogous manner. The test statistic is T = 2{log [L_{ a }(p_{ a }, F_{ a })] + log [L_{ b }(p_{ b }, F_{ b })] - log [L_{ c }(p_{ c }, F_{ c })]}, where the subscript c refers to the combined data (case and control individual), that is, common parameter values. For unrestricted parameter values, the test statistic T has an asymptotic chi-square distribution under H_{0}. However, because of the conditions imposed (F ≥ 0), T does not follow a chi-square distribution. Therefore, as outlined in the Methods section, we compute associated significance levels numerically by computer-based permutation testing. Another reason for applying permutation tests is as follows. This single test may replace the conventional two association tests (allele and genotype test) in genome-wide association studies, where 100,000s of markers are used. Test results for these markers are not independent, which is optimally taken into account in permutation tests but would be difficult to capture analytically.
Power calculations
To evaluate the performance of our new test with existing tests, we carry out power calculations under a recessive and a dominant model of disease inheritance, where we assume a functional SNP fully associated with the disease variant. Model parameters (penetrances and disease allele frequencies) are calibrated to predict a trait prevalence of 5% for each model. The proportion of affected individuals in the population whose disease is due to the given gene is fixed at 10%. The "strength" of a model is measured by the penetrance ratio, γ, where γ = 1 corresponds to the null hypothesis.
We compare three tests, our FP test, the chi-square genotype test based on a 2 × 3 table of SNP genotypes versus case and control individuals, and the chi-square allele test based on a 2 × 2 table of SNP alleles. Power calculations are carried out for a type I error (rate of false positive results) of 0.05 and assumed numbers of observations of 100 case and 100 control individuals. In practice, most researchers carry out both, the allele and the genotype test, and emphasize whichever result has a smaller p-value without correcting for the effects of multiple testing inherent in this procedure. Thus, we capture the statistical properties of this practice by formulating a test statistic, MaxGA, which is the smaller of the two p-values associated with the genotype and allele tests. Under our model assumptions, if a result is declared significant whenever either the allele test or genotype test is significant, this practice has a type I error of 0.076 for recessive traits and 0.059 for dominant traits. In our power calculations, of course, a type I error of 0.05 is imposed for the MaxGA test statistic, that is, a critical limit is chosen such that the "power" of the MaxGA statistic is equal to 0.05 under H_{0}.
Under both dominant and recessive models, the MaxGA statistic has essentially the same power as the genotype test. Thus, any superiority of the allele test, if carried out at the same time as the genotype test, is wiped out by the multiple testing correction. As the usual practice is to do both, an allele and a genotype test, the potentially higher power of the allele test is fallacious as it is entirely purchased by an increased rate of false positive results.
We also considered Risch's genotype relative risk model [14], in which the penetrances for the three disease genotypes are given by f, γf, and γ^{2}f so that this model may be viewed as being intermediate between dominant and recessive inheritance. For our assumptions on disease inheritance, this model always predicts unrestricted parameter values of F = 0 in case and slightly negative F values in control individuals, where the latter would be estimated to be zero under our restrictions. Thus, the FP test is expected to essentially default to the genotype test under Risch's penetrance structure but this is not further pursued here.
Application to published data
Test results for observed data
Data | SNP | p _{genotype} | p _{FP} | p _{allele} | F _{case} | F _{control} |
---|---|---|---|---|---|---|
AMD | rs380390 | 0.0380 | 0.0090 | 0.0056 | 0.215 | -0.073 |
rs10272438 | 1.0000 | 0.9068 | 0.0194 | 0.733 | 0.611 | |
AMD HK | rs10490924 | 0.0002 | 0.0002 | 0.0002 | 0.243 | -0.062 |
rs10504152 | 0.0058 | 0.1286 | 0.2222 | 0.132 | -0.271 | |
rs584244 | 0.1824 | 0.0996 | 0.1010 | -0.011 | -0.149 | |
PD | rs9952724 | 0.0004 | 0.0002 | 1.0000 | 0.788 | 0.022 |
rs850084 | 0.0022 | 0.0002 | 0.9932 | 0.828 | 0.243 | |
rs10963676 | 0.0058 | 0.0004 | 0.0072 | 0.817 | 0.086 | |
rs4746675 | 0.0062 | 0.0004 | 1.0000 | 0.839 | 0.048 | |
rs557074 | 0.0068 | 0.0012 | 1.0000 | 0.736 | 0.029 | |
rs1504212 | 0.0088 | 0.0014 | 1.0000 | 0.494 | -0.023 | |
rs12364577 | 0.0174 | 0.0020 | 1.0000 | 0.519 | 0.014 | |
rs1468375 | 0.0240 | 0.0042 | 0.0002 | 0.452 | -0.040 |
For each of the three studies in Table 3, all SNPs are listed that achieved an experiment-wise significance level of 0.05 or less in either one of the three association tests. Of the 13 resulting SNPs, 11 show a smaller p-value for the FP test than the genotype test and one SNP shows the same p-value. These results clearly demonstrate the usefulness of our new association test. As expected, observed unrestricted F values are larger in case than control individuals and in the latter are often negative.
Discussion
As mentioned in the introduction, researchers often look for genomic regions of increased homozygosity or autozygosity by sliding a window of fixed length across the genome. Our test offers an elegant alternative to such windows of fixed and arbitrary lengths. We propose to work with scan statistics as previously developed [21]. This method also employs a window of a fixed length (fixed number of SNPs) and determines the maximum of the sum of test statistics for all such windows in the genome, which is the scan statistic of the given length. What sets this approach apart from ad hoc approaches is that it applies different window sizes from 1 up to a specified maximum length and estimates optimal window length by maximum likelihood. An updated version of our scanstat program is available that incorporates the FP test statistic [22]. This implementation allows users to determine the most significant stretch of continguous markers with high values of the FP statistic, which we interpret as a genomic region of high IBD.
It is interesting to note unrestricted values of F predicted by our disease models. For example, for γ = 5, the recessive model predicts 0.316 in cases and -0.025 in controls. This explains why our FP test has higher power than conventional tests for recessive traits: Inbreeding coefficients tend to be strongly positive in cases and only slightly negative in controls (this implies strong deviations from HWE in cases). Our parameter restrictions disallow negative F values, which reduces "noise" in the determination of significance levels. On the other hand, the dominant model predicts unrestricted F values of -0.041 in cases and 0.002 in controls. These values are only slightly different from 0 and will become non-negative in the FP test, that is, the FP essentially defaults to the genotype test with only a slight advantage over it.
The null distribution of the FP test statistic (under F ≥ 0) is not known and would be difficult to obtain, particularly for large numbers of markers whose test results are non-independent.
Conclusion
At least for the recessive and dominant models considered here, our FP test is more powerful than allele and genotype tests. Thus, it represents an attractive alternative to these conventional tests. A potential disadvantage of the FP test might be that it requires permutation testing for an appropriate determination of p-values. However, permutation testing is one of the best approaches to correct for multiple testing in genome-wide association studies and is often carried out anyway, so the FP test does not represent an additional burden.
Methods
The resulting restricted parameter space corresponds to the area marked "F > 0" in Figure 3, that is, the area between the convex solid line and the diagonal. The restricted parameter space is only 2/3 of the surface of the unrestricted parameter space.
For power calculations, we assume disease models with two alleles and three genotypes, dd, Dd, and DD. The respective penetrances are f_{1}, f_{2}, and f_{3}, where we set f_{2} = f_{1} for recessive models and f_{2} = f_{3} for dominant models. The "strength" of a model is measured by the penetrance ratio, γ = f_{3}/f_{1}, which is very approximately equal to the odds ratio. Thus, we have three genetic parameters, p, f_{1}, and γ, which predict trait prevalence as
K = γf_{1}p^{2} + f_{1}(1 - p^{2}) for recessive traits and as
K = q^{2}f_{1} + (1-q^{2})γf_{1}, q = 1 - p, for dominant traits.
Also, the proportion of genetic cases among all affected individuals is
Q = γp^{2}/(γp^{2} + 1 - p^{2}) for recessive traits and
Q = (1 -q^{2})γ/[(1 - q^{2})γ + q^{2}] for dominant traits.
Fixing K = 0.05 and Q = 0.10 leaves one free parameter, which we vary to generate power curves. For a fixed set of genetic parameter values, Bayesian calculations yield conditional genotype frequencies in case and control individuals, from which random samples (replicates) are drawn.
For each of dominant and recessive models, with a value of the penetrance ratio, γ = 1, critical limits for test statistics are chosen so as to make the type I error for each of them equal to 0.05. That is, the critical limits are chosen such that the proportion of randomization samples exceeding this limit is equal to 0.05 (in other words, we are using the 95^{th} percentile of the computer-generated null distribution of the test statistic as the critical limit). Then power is determined for penetrance ratios ranging from 1 through 8. All power calculations were carried out based on 5,000 replicates.
Declarations
Acknowledgements
This work was supported by China NSFC grants, project numbers 30730057 (JO) and 30700442 (QRZ), and by grant MH44292 (JO) from the U.S. National Institute of Mental Health. This study used data from the SNP Database at the NINDS Human Genetics Resource Center DNA and Cell Line Repository [23], as well as clinical data. The original genotyping was performed in the laboratories of Drs. Singleton and Hardy, (NIA, LNG), Bethesda, MD USA.
Authors’ Affiliations
References
- Lander ES, Botstein D: Homozygosity mapping: a way to map human recessive traits with the DNA of inbred children. Science. 1987, 236 (4808): 1567-1570. 10.1126/science.2884728.View ArticlePubMedGoogle Scholar
- Houwen RH, Baharloo S, Blankenship K, Raeymaekers P, Juyn J, Sandkuijl LA, Freimer NB: Genome screening by searching for shared segments: mapping a gene for benign recurrent intrahepatic cholestasis. Nat Genet. 1994, 8 (4): 380-386. 10.1038/ng1294-380.View ArticlePubMedGoogle Scholar
- Broman KW, Weber JL: Long homozygous chromosomal segments in reference families from the centre d'Etude du polymorphisme humain. Am J Hum Genet. 1999, 65 (6): 1493-1500. 10.1086/302661.PubMed CentralView ArticlePubMedGoogle Scholar
- Levran O, Attwooll C, Henry RT, Milton KL, Neveling K, Rio P, Batish SD, Kalb R, Velleuer E, Barral S, Ott J, Petrini J, Schindler D, Hanenberg H, Auerbach AD: The BRCA1-interacting helicase BRIP1 is deficient in Fanconi anemia. Nat Genet. 2005, 37 (9): 931-3. Epub 2005 Aug 21.. 10.1038/ng1624.View ArticlePubMedGoogle Scholar
- Gibson J, Morton NE, Collins A: Extended tracts of homozygosity in outbred human populations. Hum Mol Genet. 2006, 15 (5): 789-95. Epub 2006 Jan 25.. 10.1093/hmg/ddi493.View ArticlePubMedGoogle Scholar
- Woods CG, Cox J, Springell K, Hampshire DJ, Mohamed MD, McKibbin M, Stern R, Raymond FL, Sandford R, Malik Sharif S, Karbani G, Ahmed M, Bond J, Clayton D, Inglehearn CF: Quantification of homozygosity in consanguineous individuals with autosomal recessive disease. Am J Hum Genet. 2006, 78 (5): 889-96. Epub 2006 Mar 21.. 10.1086/503875.PubMed CentralView ArticlePubMedGoogle Scholar
- Li LH, Ho SF, Chen CH, Wei CY, Wong WC, Li LY, Hung SI, Chung WH, Pan WH, Lee MT, Tsai FJ, Chang CF, Wu JY, Chen YT: Long contiguous stretches of homozygosity in the human genome. Hum Mutat. 2006, 27 (11): 1115-1121. 10.1002/humu.20399.View ArticlePubMedGoogle Scholar
- Holmans P: Asymptotic properties of affected-sib-pair linkage analysis. Am J Hum Genet. 1993, 52 (2): 362-374.PubMed CentralPubMedGoogle Scholar
- Vogel F, Motulsky AG: Human genetics: problems and approaches. 1997, Berlin; New York , Springer, xxxvi, 851-3rd completely rev.View ArticleGoogle Scholar
- Agarwala R, Schaffer AA, Tomlin JF: Towards a complete North American Anabaptist Genealogy II: analysis of inbreeding. Hum Biol. 2001, 73 (4): 533-545. 10.1353/hub.2001.0045.View ArticlePubMedGoogle Scholar
- Klein RJ, Zeiss C, Chew EY, Tsai JY, Sackler RS, Haynes C, Henning AK, SanGiovanni JP, Mane SM, Mayne ST, Bracken MB, Ferris FL, Ott J, Barnstable C, Hoh J: Complement factor H polymorphism in age-related macular degeneration. Science. 2005, 308 (5720): 385-9. Epub 2005 Mar 10.. 10.1126/science.1109557.PubMed CentralView ArticlePubMedGoogle Scholar
- Zou GY, Donner A: The merits of testing Hardy-Weinberg equilibrium in the analysis of unmatched case-control data: a cautionary note. Ann Hum Genet. 2006, 70 (Pt 6): 923-933.PubMedGoogle Scholar
- Online manual for FP test. [http://www.genemapping.cn/FPtest.html]
- Risch N, Merikangas K: The future of genetic studies of complex human diseases. Science. 1996, 273 (5281): 1516-1517. 10.1126/science.273.5281.1516.View ArticlePubMedGoogle Scholar
- Dewan A, Liu M, Hartman S, Zhang SS, Liu DT, Zhao C, Tam PO, Chan WM, Lam DS, Snyder M, Barnstable C, Pang CP, Hoh J: HTRA1 promoter polymorphism in wet age-related macular degeneration. Science. 2006, 314 (5801): 989-92. Epub 2006 Oct 19. 10.1126/science.1133807.View ArticlePubMedGoogle Scholar
- Fung HC, Scholz S, Matarin M, Simon-Sanchez J, Hernandez D, Britton A, Gibbs JR, Langefeld C, Stiegert ML, Schymick J, Okun MS, Mandel RJ, Fernandez HH, Foote KD, Rodriguez RL, Peckham E, De Vrieze FW, Gwinn-Hardy K, Hardy JA, Singleton A: Genome-wide genotyping in Parkinson's disease and neurologically normal controls: first stage analysis and public release of data. Lancet Neurol. 2006, 5 (11): 911-916. 10.1016/S1474-4422(06)70578-6.View ArticlePubMedGoogle Scholar
- Simon-Sanchez J, Scholz S, Del Mar Matarin M, Fung HC, Hernandez D, Gibbs JR, Britton A, Hardy J, Singleton A: Genomewide SNP assay reveals mutations underlying Parkinson disease. Hum Mutat. 2007, 9: 9-Google Scholar
- Nijenhuis A, Wilf HS: Combinatorial algorithms for computers and calculators. 1978, New York , Academic Press, xv, 302 p.-2dGoogle Scholar
- Online manual for sumstat program. [http://www.genemapping.cn/sumstat.html]
- Press WH, Teukolsky SA, Vetterling WT, Flannery BP: Numerical recipes 3rd edition: The art of scientific computing. 2007, Cambridge, UK; New York , Cambridge University Press, xxi, 1235-3rdGoogle Scholar
- Hoh J, Ott J: Scan statistics to scan markers for susceptibility genes. Proc Natl Acad Sci U S A. 2000, 97 (17): 9615-9617. 10.1073/pnas.170179197.PubMed CentralView ArticlePubMedGoogle Scholar
- Online manual for scanstat program. [http://www.genemapping.cn/scanstat.html]
- Coriell NINDS Collection. [http://ccr.coriell.org/ninds]
Copyright
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.