title: Performance of a blockwise approach in variable selection using linkage disequilibrium information.
abstract: BACKGROUND
  Genome-wide association studies (GWAS) aim at finding genetic markers that are significantly associated with a phenotype of interest. Single nucleotide polymorphism (SNP) data from the entire genome are collected for many thousands of SNP markers, leading to high-dimensional regression problems where the number of predictors greatly exceeds the number of observations. Moreover, these predictors are statistically dependent, in particular due to linkage disequilibrium (LD). We propose a three-step approach that explicitly takes advantage of the grouping structure induced by LD in order to identify common variants which may have been missed by single marker analyses (SMA). In the first step, we perform a hierarchical clustering of SNPs with an adjacency constraint using LD as a similarity measure. In the second step, we apply a model selection approach to the obtained hierarchy in order to define LD blocks. Finally, we perform Group Lasso regression on the inferred LD blocks. We investigate the efficiency of this approach compared to state-of-the art regression methods: haplotype association tests, SMA, and Lasso and Elastic-Net regressions.
  
  
  RESULTS
  Our results on simulated data show that the proposed method performs better than state-of-the-art approaches as soon as the number of causal SNPs within an LD block exceeds 2. Our results on semi-simulated data and a previously published HIV data set illustrate the relevance of the proposed method and its robustness to a real LD structure. The method is implemented in the R package BALD (Blockwise Approach using Linkage Disequilibrium), available from http://www.math-evry.cnrs.fr/publications/logiciels .
  
  
  CONCLUSIONS
  Our results show that the proposed method is efficient not only at the level of LD blocks by inferring well the underlying block structure but also at the level of individual SNPs. Thus, this study demonstrates the importance of tailored integration of biological knowledge in high-dimensional genomic studies such as GWAS.
authors.affiliation: Laboratoire de Mathématiques et Modélisation d'Evry (LaMME), Université d'Evry-Val-d'Essonne/UMR CNRS 8071/ENSIIE/USC INRA, Evry, France. alia.dehman@genopole.cnrs.fr.
authors.forename: Alia
authors.initials: A
authors.lastname: Dehman
authors.affiliation: Laboratoire de Mathématiques et Modélisation d'Evry (LaMME), Université d'Evry-Val-d'Essonne/UMR CNRS 8071/ENSIIE/USC INRA, Evry, France. christophe.ambroise@genopole.cnrs.fr.
authors.forename: Christophe
authors.initials: C
authors.lastname: Ambroise
authors.affiliation: Laboratoire de Mathématiques et Modélisation d'Evry (LaMME), Université d'Evry-Val-d'Essonne/UMR CNRS 8071/ENSIIE/USC INRA, Evry, France. pierre.neuvial@genopole.cnrs.fr.
authors.forename: Pierre
authors.initials: P
authors.lastname: Neuvial
mesh_terms: D000465:Algorithms; D005819:Genetic Markers / Q000235:genetics; D055106:Genome-Wide Association Study / Q000379:methods*; D006239:Haplotypes / Q000235:genetics*; D006801:Humans; D015810:Linkage Disequilibrium; D008962:Models, Theoretical; D020641:Polymorphism, Single Nucleotide / Q000235:genetics*
chemical_list: D005819:Genetic Markers