Abstract:<p>To understand how organisms adapt to their environment, a gene‐environmental association (GEA) analysis is commonly conducted. GEA methods based on mixed models, such as linear latent factor mixed models (LFMM) and LFMM2, have grown in popularity for their robust performance in terms of power and computational speed. However, it is unclear how the assumption of a Gaussian distribution for the response variables influences model performance. In this paper, we develop a generalized linear model (GLM) that allows for non‐Gaussian distribution in the genotypic response variables, and treatment of multiallelic nucleotide polymorphisms. Moreover, this multinomial logistic regression model (MLR) is combined with an admixture‐based model or principal components analysis to correct for population structure (MLR‐ADM and MLR‐PC). Using simulations, we evaluate the type I error, false discovery rates (FDR), and power to detect selected SNPs, to guide model choice and best practices. With genomic control, MLR‐PC and LFMM2 have similar type I error, FDRs, and power when analyzing biallelic SNPs, while dramatically outperforming models not accounting for population structure. Differences in performance occur under continuous population structure where MLR‐PC outperforms LFMM/LFMM2, especially when a larger number of clusters or triallelic SNPs are analyzed. The Human Genome Diversity Project (HGDP) dataset shows that both MLR‐PC and LFMM2 control the inflation of P‐values. Analysis of the 1000 Genome Project Phase III dataset illustrates that MLR‐PC and LFMM2 produce consistent results for most significant SNPs, while MLR‐PC discovered additional SNPs corresponding to certain genes, suggesting MLR‐PC may be a useful alternative to GEA inference.</p>

Efficient Genomic Control for Mixed Model Associations in Large-scale Population

Optimizing genomic control in mixed model associations with binary diseases

An Improved Linear Mixed Model for Multivariate Genome-Wide Association Studies

A SUPER powerful method for genome wide association study.

Mixed Linear Model Approaches of Association Mapping for Complex Traits Based on Omics Variants

Association Test Between Haplotypes and Longitudinal Traits in Complex Pedigrees.

Regression-based Approach for Testing the Association Between Multi-Region Haplotype Configuration and Complex Trait

Efficient and Accurate Multiple-Phenotype Regression Method for High Dimensional Data Considering Population Structure

A Sparse Graph-Structured Lasso Mixed Model for Genetic Association with Confounding Correction

Efficient Multivariate Analysis Algorithms for Longitudinal Genome-wide Association Studies

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

Mixed Linear Model Approaches for Analyzing Genetic Models of Complex Quantitative Traits

Asymptotically Exact Fit for Linear Mixed Model in Genetic Association Studies

A rapid epistatic mixed-model association analysis by linear retransformations of genomic estimated values.

Rapid epistatic mixed-model association studies by controlling multiple polygenic effects

FastBiCmrMLM: a fast and powerful compressed variance component mixed logistic model for big genomic case-control genome-wide association study

Detect and Adjust for Population Stratification in Population-Based Association Study Using Genomic Control Markers: an Application of Affymetrix Genechip® Human Mapping 10K Array

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

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

Hierarchical Generalized Linear Mixed Model for Genome-wide Association Analysis

Tradeoffs of Linear Mixed Models in Genome-wide Association Studies