Abstract:AbstractWith the availability of genotyping data of very large samples, there is an increasing need for tools that can efficiently identify genetic relationships among all individuals in the sample. One fundamental measure of genetic relationship of a pair of individuals is identity by descent (IBD), chromosomal segments that are shared among two individuals due to common ancestry. However, the efficient identification of IBD segments among a large number of genotyped individuals is a challenging computational problem. Most existing methods are not feasible for even thousands of individuals because they are based on pairwise comparisons of all individuals and thus scale up quadratically with sample size. Some methods, such as GERMLINE, use fast dictionary lookup of short seed sequence matches to achieve a near-linear time efficiency. However, the number of short seed matches often scales up super-linearly in real population data.In this paper we describe a new approach for IBD detection. We take advantage of an efficient population genotype index, Positional BWT (PBWT), by Richard Durbin. PBWT achieves linear time query of perfectly identical subsequences among all samples. However, the original PBWT is not tolerant to genotyping errors which often interrupt long IBD segments into short fragments. We introduce a randomized strategy by running PBWTs over random projections of the original sequences. To boost the detection power we run PBWT multiple times and merge the identified IBD segments through interval tree algorithms. Given a target IBD segment length, RaPID adjust parameters to optimize detection power and accuracy.Simulation results proved that our tool (RaPID) achieves almost linear scaling up to sample size and is orders of magnitude faster than GERMLINE. At the same time, RaPID maintains a detection power and accuracy comparable to existing mainstream algorithms, GERMLINE and IBDseq. Running multiple times with various target detection lengths over the 1000 Genomes Project data, RaPID can detect population events at different time scales. With our tool, it is feasible to identify IBDs among hundreds of thousands to millions of individuals, a sample size that will become reality in a few years.

Fast simulation of identity-by-descent segments

Identity-by-descent in large samples

Ultra-Fast Identity by Descent Detection in Biobank-Scale Cohorts Using Positional Burrows-Wheeler Transform.

Conflation of short identity-by-descent segments bias their inferred length distribution

Simulating pedigrees ascertained on the basis of observed IBD sharing

Estimating effective population size trajectories from time-series Identity-by-Descent (IBD) segments

Estimating evolutionary and demographic parameters via ARG-derived IBD

Modeling recent positive selection using identity-by-descent segments

Benchmarking and Optimization of Methods for the Detection of Identity-By-Descent in High-Recombining

Rapid detection of identity-by-descent tracts for mega-scale datasets

A renewal theory approach to IBD sharing

Estimation of demography and mutation rates from one million haploid genomes

RaPID-Query for fast identity by descent search and genealogical analysis

Modeling recent positive selection in Americans of European ancestry

Open-source benchmarking of IBD segment detection methods for biobank-scale cohorts

Minimum Segmentation for Pan-genomic Founder Reconstruction in Linear Time

Faster inference of complex demographic models from large allele frequency spectra

A Probabilistic Method for Estimating the Sharing of Identity by Descent for Populations with Migration

Decoding coalescent hidden Markov models in linear time

Accurate genome-wide phasing from IBD data

Sequential importance sampling for multi-resolution Kingman-Tajima coalescent counting