Abstract:There is an increasing demand for software that can handle an arbitrary number of linked markers in forensic genetics; primarily with application to inference of relationships and direct matching but also in applications such as ancestry inference and mixture interpretation. With the emergence of sequencing technologies, denser sets of SNP markers are generated and analyzed. Additionally, sequence data of low quality and quantity DNA generate uncertainty about the underlying true genotype. We provide an efficient implementation of a general model for pedigree likelihood computations with genetic marker data using a three-layered approach. The top and first layer is the population model where allele frequencies and population substructure are accounted for. The second layer is the inheritance model which efficiently handles linked markers using an IBD model. The third and bottom layer is the observational level where we model the likelihood of the true genotype given underlying reads as well as parameters for errors.We exemplify the utility of our implementation as well as provide validation according to guidelines established by the ISFG using a combination of two published SNP panels. We demonstrate that computations are feasible for panels encompassing 10,000 markers and we argue that, due to the properties of the underlying algorithm, extending the number of markers will result in a linear increase in computation time. In addition we study the impact of parameters used in our model and suggest some guidelines pertaining to their values. The results demonstrate that a probabilistic model for low coverage sequence read data is needed instead of relying on an a threshold based genotype and applying our general model for inference of relationships on a real case can be superior, i.e. higher information content, to other methods relying on either fixed genotypes with low quality sequence data or simple pair wise relationship tests.In summary, the implementation, FamLink2 (freely available at https://famlink.se ), can jointly handle genetic linkage, genotype uncertainty and population substructure for an arbitrary pedigree with data for any number of individuals. Whereas the current study will focus on calculations disregarding mutations, FamLink2 has the ability to model mutations for certain built-in pedigrees.

SNPest: a probabilistic graphical model for estimating genotypes

A general approach to single-nucleotide polymorphism discovery

An empirical Bayes method for genotyping and SNP detection using multi-sample next-generation sequencing data.

Distangsd: Fast and Accurate Inference of Genetic Distances for Next-Generation Sequencing Data

SNP Calling Using Genotype Model Selection on High-Throughput Sequencing Data

Genotyping Polyploids from Messy Sequencing Data

A Fast and Flexible Statistical Model for Large-Scale Population Genotype Data: Applications to Inferring Missing Genotypes and Haplotypic Phase

SNP Genotype Imputation in Forensics—A Performance Study

Large-scale Genotyping of Complex DNA

Simpute: an Efficient Solution for Dense Genotypic Data

A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data

FamLink2 - A comprehensive tool for likelihood computations in pedigrees analyses involving linked DNA markers accounting for genotype uncertainties

Genomic relationships computed from either next-generation sequence or array SNP data

Simpute: A Simple Genotype Imputation Method

A Probabilistic Model to Predict Clinical Phenotypic Traits from Genome Sequencing

GGTyper: genotyping complex structural variants using short-read sequencing data

MaCH: using sequence and genotype data to estimate haplotypes and unobserved genotypes

FastImpute: A Baseline for Open-source, Reference-Free Genotype Imputation Methods -- A Case Study in PRS313

A Model of Higher Accuracy for the Individual Haplotyping Problem Based on Weighted Snp Fragments and Genotype with Errors

Single-sample SNP Detection by Empirical Bayesian Method Using Next Generation Sequencing Data