Abstract:Background: Developing binary classification rules based on SNP observations has been a major challenge for many modern bioinformatics applications, e.g., predicting risk of future disease events in complex conditions such as cancer. Small-sample, high-dimensional nature of SNP data, weak effect of each SNP on the outcome, and highly non-linear SNP interactions are several key factors complicating the analysis. Additionally, SNPs take a finite number of values which may be best understood as ordinal or categorical variables, but are treated as continuous ones by many algorithms. Methods: We use the theory of high dimensional model representation (HDMR) to build appropriate low dimensional glass-box models, allowing us to account for the effects of feature interactions. We compute the second order HDMR expansion of the log-likelihood ratio to account for the effects of single SNPs and their pairwise interactions. We propose a regression based approach, called linear approximation for block second order HDMR expansion of categorical observations (LABS-HDMR-CO), to approximate the HDMR coefficients. We show how HDMR can be used to detect pairwise SNP interactions, and propose the fixed pattern test (FPT) to identify statistically significant pairwise interactions. Results: We apply LABS-HDMR-CO and FPT to synthetically generated HAPGEN2 data as well as to two GWAS cancer datasets. In these examples LABS-HDMR-CO enjoys superior accuracy compared with several algorithms used for SNP classification, while also taking pairwise interactions into account. FPT declares very few significant interactions in the small sample GWAS datasets when bounding false discovery rate (FDR) by 5%, due to the large number of tests performed. On the other hand, LABS-HDMR-CO utilizes a large number of SNP pairs to improve its prediction accuracy. In the larger HAPGEN2 dataset FTP declares a larger portion of SNP pairs used by LABS-HDMR-CO as significant. Conclusion: LABS-HDMR-CO and FPT are interesting methods to design prediction rules and detect pairwise feature interactions for SNP data. Reliably detecting pairwise SNP interactions and taking advantage of potential interactions to improve prediction accuracy are two different objectives addressed by these methods. While the large number of potential SNP interactions may result in low power of detection, potentially interacting SNP pairs, of which many might be false alarms, can still be used to improve prediction accuracy.

Efficient penalized generalized linear mixed models for variable selection and genetic risk prediction in high-dimensional data

A Non-Parametric Method for Building Predictive Genetic Tests on High-Dimensional Data

Combining Sparse Group Lasso and Linear Mixed Model Improves Power to Detect Genetic Variants Underlying Quantitative Traits

glmmPen: High Dimensional Penalized Generalized Linear Mixed Models

A comprehensive analysis comparing linear and generalized linear models in detecting adaptive SNPs

A generalized linear mixed model association tool for biobank-scale data

Analyzing Association Mapping in Pedigree-Based GWAS Using a Penalized Multitrait Mixed Model.

A Penalized Multi-trait Mixed Model for Association Mapping in Pedigree-based GWAS

High dimensional model representation of log likelihood ratio: binary classification with SNP data

An Efficient Unified Model for Genome-Wide Association Studies and Genomic Selection

A Bayesian Framework for Generalized Linear Mixed Modeling Identifies New Candidate Loci for Late-Onset Alzheimer's Disease

Efficient and powerful familywise error control in genome-wide association studies using generalized linear models

Improving GWAS discovery and genomic prediction accuracy in biobank data

Hierarchical Generalized Linear Mixed Model for Genome-wide Association Analysis

Efficient computation of high-dimensional penalized generalized linear mixed models by latent factor modeling of the random effects

Improved polygenic prediction by Bayesian multiple regression on summary statistics

Ordered multinomial regression for genetic association analysis of ordinal phenotypes at Biobank scale

Effective Genetic Risk Prediction Using Mixed Models

Hierarchical selection of genetic and gene by environment interaction effects in high-dimensional mixed models

High-throughput and Efficient Multilocus Genome-Wide Association Study on Longitudinal Outcomes.

Winner's Curse Correction and Variable Thresholding Improve Performance of Polygenic Risk Modeling Based on Genome-Wide Association Study Summary-Level Data