Abstract:Abstract Background Advancements in statistical methods and sequencing technology have led to numerous novel discoveries in human genetics in the past two decades. Among phenotypes of interest, most attention has been given to studying genetic associations with continuous or binary traits. Efficient statistical methods have been proposed and are available for both types of traits under different study designs. However, for multinomial categorical traits in related samples, there is a lack of efficient statistical methods and software. Results We propose an efficient score test to analyze a multinomial trait in family samples, in the context of genome-wide association/sequencing studies. An alternative Wald statistic is also proposed. We also extend the methodology to be applicable to ordinal traits. We performed extensive simulation studies to evaluate the type-I error of the score test, Wald test compared to the multinomial logistic regression for unrelated samples, under different allele frequency and study designs. We also evaluate the power of these methods. Results show that both the score and Wald tests have a well-controlled type-I error rate, but the multinomial logistic regression has an inflated type-I error rate when applied to family samples. We illustrated the application of the score test with an application to the Framingham Heart Study to uncover genetic variants associated with diabesity, a multi-category phenotype. Conclusion Both proposed tests have correct type-I error rate and similar power. However, because the Wald statistics rely on computer-intensive estimation, it is less efficient than the score test in terms of applications to large-scale genetic association studies. We provide computer implementation for both multinomial and ordinal traits.

Blindly using Wald's test can miss rare disease-causal variants in case-control association studies.

Rare variant association tests for ancestry-matched case-control data based on conditional logistic regression

A stepwise likelihood ratio test procedure for rare variant selection in case–control studies

A Shrinkage Method for Testing the Hardy–Weinberg Equilibrium in Case‐Control Studies

A Comparison of the Likelihood Ratio Test and the Variance‐Stabilising Transformation‐Based Tests for Detecting Association of Rare Variants

A nonparametric test for association with multiple loci in the retrospective case-control study.

Pooled Association Tests for Rare Genetic Variants: A Review and Some New Results

Block-based Association Tests for Rare Variants Using Kullback–Leibler Divergence

Nonparametric Tests of Associations with Disease Based on U‐Statistics

Efficient Utilization of Rare Variants for Detection of Disease-Related Genomic Regions

Kullback–Leibler Distance Methods for Detecting Disease Association with Rare Variants from Sequencing Data

Comparison of dimension reduction based logistic regression models for case control genome wide association study： principal components analysis vs. partial least squares

Single-Locus Genetic Association Analysis By Ordinal Tests

Testing Genetic Association with Rare Variants in Admixed Populations.

Enhancing the Power to Detect Low-Frequency Variants in Genome-Wide Screens

Association analysis of rare variants with quantitative trait based on minimum P-value

Gains in Power for Exhaustive Analyses of Haplotypes Using Variable-Sized Sliding Window Strategy: a Comparison of Association-Mapping Strategies.

Adjustment of Familial Relatedness in Association Test for Rare Variants

Meta-analysis of Gene-Level Tests for Rare Variant Association.

Genetic association tests in family samples for multi-category phenotypes

[Family-Based Association Tests for Rare Variants].