Abstract:Abstract Motivation Estimating the individual inbreeding coefficient and pairwise kinship is an important problem in human genetics (e.g. in disease mapping) and in animal and plant genetics (e.g. inbreeding design). Existing methods, such as sample correlation-based genetic relationship matrix, KING, and UKin, are either biased, or not able to estimate inbreeding coefficients, or produce a large proportion of negative estimates that are difficult to interpret. This limitation of existing methods is partly due to failure to explicitly model inbreeding. Since all humans are inbred to various degrees by virtue of shared ancestries, it is prudent to account for inbreeding when inferring kinship between individuals. Results We present “Kindred,” an approach that estimates inbreeding and kinship by modeling latent identity-by-descent states that accounts for all possible allele sharing—including inbreeding—between two individuals. Kindred used non-negative least squares method to fit the model, which not only increases computation efficiency compared to the maximum likelihood method, but also guarantees non-negativity of the kinship estimates. Through simulation, we demonstrate the high accuracy and non-negativity of kinship estimates by Kindred. By selecting a subset of SNPs that are similar in allele frequencies across different continental populations, Kindred can accurately estimate kinship between admixed samples. In addition, we demonstrate that the realized kinship matrix estimated by Kindred is effective in reducing genomic control values via linear mixed model in genome-wide association studies. Finally, we demonstrate that Kindred produces sensible heritability estimates on an Australian height dataset. Availability and implementation Kindred is implemented in C with multi-threading. It takes vcf file or stream as input and works seamlessly with bcftools. Kindred is freely available at https://github.com/haplotype/kindred.

biMM: efficient estimation of genetic variances and covariances for cohorts with high-dimensional phenotype measurements

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

lme4breeding: enabling genetic evaluation in the era of genomic data

Multivariate Mixed Models in Ecology and Evolutionary biology: Inferences and implementation in R

Asymptotically Exact Fit for Linear Mixed Model in Genetic Association Studies

Comparison of REML methods for the study of phenome-wide genetic variation

Principal component analysis revisited: fast multitrait genetic evaluations with smooth convergence

Optimal Estimation Of Genetic Relatedness In High-Dimensional Linear Models

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

LARGE-SCALE MULTIVARIATE SPARSE REGRESSION WITH APPLICATIONS TO UK BIOBANK

SCAMPI: A scalable statistical framework for genome-wide interaction testing harnessing cross-trait correlations

Estimating and accounting for unobserved covariates in high dimensional correlated data

Multidimensional Stochastic Process Model and its Applications to Analysis of Longitudinal Data with Genetic Information

A high-performance computing toolset for relatedness and principal component analysis of SNP data

Coupled mixed model for joint genetic analysis of complex disorders with two independently collected data sets

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

Estimation of Jacquard's genetic identity coefficients with bi-allelic variants by constrained least-squares

Subsampling Technique to Estimate Variance Component for UK-Biobank Traits

Estimation of Jacquard’s genetic identity coefficients with bi-allelic variants by constrained least-squares

Optimal Estimation of Co-heritability in High-dimensional Linear Models

Estimation of inbreeding and kinship coefficients via latent identity-by-descent states