Linkage disequilibrium structure and its impact on the localization of a candidate functional mutation

We have used the unblinded MG1/Q1 Genetic Analysis Workshop 12 simulated data as a model system for investigating the use of linkage disequilibrium structure and simple genotype-phenotype associations to identify candidate functional mutations within a gene of interest. Analysis of the pattern of pairwise linkage disequilibrium indicated three groups of single-nucleotide polymorphisms for which the linkage disequilibrium was high between sites within a group, but lower between sites of different groups. Using linear regression to predict levels of the trait Q1 showed that the known functional site, 5782, was usually not the best genetic predictor of Q1, but sites belonging to the same group as 5782 (i.e., group 2) were always included in the prediction model. In 49 out of the 50 replicates, the functional site was not the best predictor of the trait. Finally, more detailed analyses demonstrate that the relationship between the adjusted R2 for the marker in the prediction model and its disequilibrium with 5782 was linear with the intercept at the origin and terminating at the R2 value for the known functional mutation when the disequilibrium is maximal. These data indicate that simple association studies will not identify the functional mutation, but rather will identify candidate functional mutations that are in very tight linkage disequilibrium with the functional mutation.