Copyright © 2007 The American Society of Human Genetics. All rights reserved.
The American Journal of Human Genetics, Volume 81, Issue 2, 234-242, 1 August 2007
doi:10.1086/519221
Article
George Ayodoa, Alkes L. Pricea, Alon Keinana, Arthur Ajwangb, Michael F. Otienob, Alloys S.S. Oragob, Nick Pattersona and David Reicha,
, 
a From the Department of Genetics, Harvard Medical School, and Broad Institute of Harvard and MIT, Boston
b Department of Pre-Clinical Sciences, Kenyatta University, Nairobi
Address for correspondence and reprints: Dr. David Reich, Harvard Medical School, Department of Genetics, 77 Avenue Louis Pasteur, New Research Building, Boston, MA 02115Abstract
Statistical power to detect disease variants can be increased by weighting candidates by their evidence of natural selection. To demonstrate that this theoretical idea works in practice, we performed an association study of 10 putative resistance variants in 471 severe malaria cases and 474 controls from the Luo in Kenya. We replicated associations at HBB (P=.0008) and CD36 (P=.03) but also showed that the same variants are unusually differentiated in frequency between the Luo and Yoruba (who historically have been exposed to malaria) and the Masai and Kikuyu (who have not been exposed). This empirically demonstrates that combining association analysis with evidence of natural selection can increase power to detect risk variants by orders of magnitude—up to P=.000018 for HBB and P=.00043 for CD36.
| Haplotype-Based Association Analysis via Variance-Components Score Test The American Journal of Human Genetics, Volume 81, Issue 5, 1 November 2007, Pages 927-938 Jung-Ying Tzeng and Daowen Zhang Abstract Haplotypes provide a more informative format of polymorphisms for genetic association analysis than do individual single-nucleotide polymorphisms. However, the practical efficacy of haplotype-based association analysis is challenged by a trade-off between the benefits of modeling abundant variation and the cost of the extra degrees of freedom. To reduce the degrees of freedom, several strategies have been considered in the literature. They include (1) clustering evolutionarily close haplotypes, (2) modeling the level of haplotype sharing, and (3) smoothing haplotype effects by introducing a correlation structure for haplotype effects and studying the variance components (VC) for association. Although the first two strategies enjoy a fair extent of power gain, empirical evidence showed that VC methods may exhibit only similar or less power than the standard haplotype regression method, even in cases of many haplotypes. In this study, we report possible reasons that cause the underpowered phenomenon and show how the power of the VC strategy can be improved. We construct a score test based on the restricted maximum likelihood or the marginal likelihood function of the VC and identify its nontypical limiting distribution. Through simulation, we demonstrate the validity of the test and investigate the power performance of the VC approach and that of the standard haplotype regression approach. With suitable choices for the correlation structure, the proposed method can be directly applied to unphased genotypic data. Our method is applicable to a wide-ranging class of models and is computationally efficient and easy to implement. The broad coverage and the fast and easy implementation of this method make the VC strategy an effective tool for haplotype analysis, even in modern genomewide association studies. Abstract | | |
| Intragenic Cis and Trans Modification of Genetic Susceptibility in DYT1 Torsion Dystonia The American Journal of Human Genetics, Volume 80, Issue 6, 1 June 2007, Pages 1188-1193 Neil J. Risch, Susan B. Bressman, Geetha Senthil and Laurie J. Ozelius Abstract A GAG deletion in the DYT1 gene is a major cause of early-onset dystonia, but clinical disease expression occurs in only 30% of mutation carriers. To gain insight into genetic factors that may influence penetrance, we evaluated three DYT1 single-nucleotide polymorphisms, including D216H, a coding-sequence variation that moderates the effects of the DYT1 GAG deletion in cellular models. We tested DYT1 GAG-deletion carriers with (n=119) and without (n=113) clinical signs of dystonia and control individuals (n=197) and found the frequency of the 216H allele to be increased in GAG-deletion carriers without dystonia and to be decreased in carriers with dystonia, compared with the control individuals. Analysis of haplotypes demonstrated a highly protective effect of the H allele in trans with the GAG deletion; there was also suggestive evidence that the D216 allele in cis is required for the disease to be penetrant. Our findings establish, for the first time, a clinically relevant gene modifier of DYT1. Abstract | | |
| Case-Control Association Testing with Related Individuals: A More Powerful Quasi-Likelihood Score Test The American Journal of Human Genetics, Volume 81, Issue 2, 1 August 2007, Pages 321-337 Timothy Thornton and Mary Sara McPeek Abstract We consider the problem of genomewide association testing of a binary trait when some sampled individuals are related, with known relationships. This commonly arises when families sampled for a linkage study are included in an association study. Furthermore, power to detect association with complex traits can be increased when affected individuals with affected relatives are sampled, because they are more likely to carry disease alleles than are randomly sampled affected individuals. With related individuals, correlations among relatives must be taken into account, to ensure validity of the test, and consideration of these correlations can also improve power. We provide new insight into the use of pedigree-based weights to improve power, and we propose a novel test, the MQLS test, which, as we demonstrate, represents an overall, and in many cases, substantial, improvement in power over previous tests, while retaining a computational simplicity that makes it useful in genomewide association studies in arbitrary pedigrees. Other features of the MQLS are as follows: (1) it is applicable to completely general combinations of family and case-control designs, (2) it can incorporate both unaffected controls and controls of unknown phenotype into the same analysis, and (3) it can incorporate phenotype data about relatives with missing genotype data. The methods are applied to data from the Genetic Analysis Workshop 14 Collaborative Study of the Genetics of Alcoholism, where the MQLS detects genomewide significant association (after Bonferroni correction) with an alcoholism-related phenotype for four different single-nucleotide polymorphisms: tsc1177811 (P=5.9×10−7), tsc1750530 (P=4.0×10−7), tsc0046696 (P=4.7×10−7), and tsc0057290 (P=5.2×10−7) on chromosomes 1, 16, 18, and 18, respectively. Three of these four significant associations were not detected in previous studies analyzing these data. Abstract | | |