Abstract:Abstract Motivation Sparse regularized regression methods are now widely used in genome-wide association studies (GWAS) to address the multiple testing burden that limits discovery of potentially important predictors. Linear mixed models (LMMs) have become an attractive alternative to principal components (PCs) adjustment to account for population structure and relatedness in high-dimensional penalized models. However, their use in binary trait GWAS rely on the invalid assumption that the residual variance does not depend on the estimated regression coefficients. Moreover, LMMs use a single spectral decomposition of the covariance matrix of the responses, which is no longer possible in generalized linear mixed models (GLMMs). Results We introduce a new method called pglmm, a penalized GLMM that allows to simultaneously select genetic markers and estimate their effects, accounting for between-individual correlations and binary nature of the trait. We develop a computationally efficient algorithm based on penalized quasi-likelihood estimation that allows to scale regularized mixed models on high-dimensional binary trait GWAS. We show through simulations that when the dimensionality of the relatedness matrix is high, penalized LMM and logistic regression with PC adjustment fail to select important predictors, and have inferior prediction accuracy compared to pglmm. Further, we demonstrate through the analysis of two polygenic binary traits in a subset of 6731 related individuals from the UK Biobank data with 320K SNPs that our method can achieve higher predictive performance, while also selecting fewer predictors than a sparse regularized logistic lasso with PC adjustment. Availability and implementation Our Julia package PenalizedGLMM.jl is publicly available on github: https://github.com/julstpierre/PenalizedGLMM. Supplementary information Supplementary data are available at Bioinformatics online.

HighDimMixedModels.jl: Robust High Dimensional Mixed Models across Omics Data

Expectile regression for analyzing heteroscedasticity in high dimension

An Improved Algorithm for High-Dimensional Continuous Threshold Expectile Model with Variance Heterogeneity

Application of a genomic model for high-dimensional chemometric analysis

glmmPen: High Dimensional Penalized Generalized Linear Mixed Models

Sparse high-dimensional linear mixed modeling with a partitioned empirical Bayes ECM algorithm

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

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

High-dimensional variable selection accounting for heterogeneity in regression coefficients across multiple data sources

Non-Concave Penalization in Linear Mixed-Effects Models and Regularized Selection of Fixed Effects

SCAD-Ridge Penalized Likelihood Estimators for Ultra-high Dimensional Models

Variable Selection for High-dimensional Generalized Linear Models using an Iterated Conditional Modes/Medians Algorithm

Spike-and-Slab LASSO Generalized Additive Models and Scalable Algorithms for High-Dimensional Data Analysis

An enriched approach to combining high-dimensional genomic and low-dimensional phenotypic data

hdm: High-Dimensional Metrics

The use of the EM algorithm for regularization problems in high-dimensional linear mixed-effects models

Fast selection of nonlinear mixed effect models using penalized likelihood

High-dimensional genome-wide association study and misspecified mixed model analysis

Ultra high-dimensional semiparametric longitudinal data analysis

Inference for High-Dimensional Linear Mixed-Effects Models: A Quasi-Likelihood Approach

Partial least square based approaches for high-dimensional linear mixed models