Statistical Methods For Detecting Genetic Risk Factors of a Disease with Applications to Genome-Wide Association Studies

Doctoral thesis English OPEN
Ali, Fadhaa ; Ministry of Higher Education and Scientific Research- Republic of Iraq (2015)
  • Subject: QA276

This thesis aims to develop various statistical methods for analysing the data derived from genome wide association studies (GWAS).\ud The GWAS often involves genotyping individual human genetic variation, using high-throughput genome-wide single nucleotide polymorphism (SNP) arrays, in thousands of individuals and testing for association between those variants and a given disease under the assumption of common disease/common variant.\ud Although GWAS have identified many potential genetic factors in the genome that affect the risks to complex\ud diseases, there is still much of the genetic heritability that remains unexplained. The power of\ud detecting new genetic risk variants can be improved by considering multiple genetic variants simultaneously with novel statistical methods.\ud Improving the analysis of the GWAS data has received much attention from statisticians and other scientific researchers over the past decade. \ud \ud There are several challenges arising in analysing the GWAS data. First, determining the risk SNPs might be difficult due to non-random correlation between SNPs that can inflate type I and II errors in statistical inference. When a group of SNPs are considered together in the context of haplotypes/genotypes, the distribution of the haplotypes/genotypes is sparse, which makes it difficult to detect risk haplotypes/genotypes in terms of disease penetrance. \ud \ud In this work, we proposed four new methods to identify risk haplotypes/genotypes based on their frequency differences between cases and controls. To evaluate the performances of our methods, we simulated datasets under wide range of scenarios according to both retrospective and prospective designs. \ud \ud In the first method, we first reconstruct haplotypes by using unphased genotypes, followed by clustering and thresholding the inferred haplotypes into risk and non-risk groups with a two-component binomial-mixture model. In the method, the parameters were estimated by using the modified Expectation-Maximization algorithm, where the maximisation step was replaced the posterior sampling of the component parameters. We also elucidated the relationships between risk and non-risk haplotypes under different modes of inheritance and genotypic relative risk.\ud \ud In the second method, we fitted a three-component mixture model to genotype data directly, followed by an odds-ratio thresholding.\ud \ud In the third method, we combined the existing haplotype reconstruction software PHASE and permutation method to infer risk haplotypes. \ud \ud In the fourth method, we proposed a new way to score the genotypes by clustering and combined it with a logistic regression approach to infer risk haplotypes. \ud \ud The simulation studies showed that the first three methods outperformed the multiple testing method of (Zhu, 2010) in terms of average specificity and sensitivity (AVSS) in all scenarios considered. The logistic regression methods also outperformed the standard logistic regression method. \ud \ud We applied our methods to two GWAS datasets on coronary artery disease (CAD) and hypertension (HT), detecting several new risk haplotypes and recovering a number of the existing disease-associated genetic variants in the literature.
  • References (11)
    11 references, page 1 of 2

    [1] Balding, D.J., Bishop, M. and Cannings, C. (2007). Handbook of Statistical Genetics. John Wiley and Sons, Ltd.

    [2] Binder, H., Müller, T., Schwender, H., Golka, K., Steffens, M., Hengstler, J. G., Ickstadt, K. and Schumacher, M. (2012). Cluster-localized sparse logistic regression for SNP data. Statistical Applications in Genetics and Molecular Biology, 11(4).

    [3] Breslow, N. (1996). Statistics in epidemiology: The case-control study. Journal of the American Statistical Association, 91(433), 14-28.

    [4] Breslow, N. and Zhao, L. (1988). Logistic-Regression for Stratified Case Control Studies. Biometrics, 44(3), 891-899.

    [21] Hindorff, L. A., Sethupathy, P., Junkins, H. A., Ramos, E. M., Mehta, J. P., Collins, F. S. and Manolio, T. A. (2009). Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proceedings of the National Academy of Sciences of the United States of America, 106(23), 9362-9367.

    [22] Huang, Y. H., Lee, M. H., Chen, W. J. and Hsiao, C. K (2011). Using an Uncertainty-Coding matrix in bayesian regression models for HaplotypeSpecific risk detection in family association studies. Plos One, 6(7).

    [23] Hudson, R. R. (2002). Generating samples under a Wright-Fisher neutral model. Bioinformatics, 18, 337-8.

    [24] Igo, R. P., Jr., Li, J. and Goddard, K. A. B. (2009). Association Mapping by Generalized Linear Regression With Density-Based Haplotype Clustering. Genetic Epidemiology, 33(1), 16-26.

    [25] Kang, G., Yue, W., Zhang, J., Huebner, M., Zhang, H., Ruan, Y., Lu, T., Ling, Y., Zuo, Y. and Zhang, D., (2008). Two-stage designs to identify the effects of SNP combinations on complex diseases. Journal of Human Genetics, 53(8), 739-746.

    [26] Karlis, D. and Xekalaki, E. (2003). Choosing initial values for the EM algorithm for finite mixtures. Comput. Stat. & Data Ana. , 41, 577-590.

  • Metrics
    No metrics available
Share - Bookmark