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 characterization of segmental duplications in genome assemblies

Structural polymorphism and diversity of human segmental duplications

Identifying and Removing Haplotypic Duplication in Primary Genome Assemblies

DNA copy number analysis of fresh and formalin-fixed specimens by shallow whole-genome sequencing with identification and exclusion of problematic regions in the genome assembly

Scanning The Human Genome At Kilobase Resolution

WGDdetector: a Pipeline for Detecting Whole Genome Duplication Events Using the Genome or Transcriptome Annotations

Detecting Dispersed Duplications in High-Throughput Sequencing Data using a Database-Free Approach.

Correcting Genomic Deletion Calls With Complex Boundaries From Next Generation Sequencing Data

Copy Number Variants and Segmental Duplications Show Different Formation Signatures

Robust and Rapid Algorithms Facilitate Large-Scale Whole Genome Sequencing Downstream Analysis in an Integrative Framework

Detecting genomic deletions from high-throughput sequence data with unsupervised learning

Detecting Complex Indels With Wide Length-Spectrum From The Third Generation Sequencing Data

Increased mutation and gene conversion within human segmental duplications

SvABA: genome-wide detection of structural variants and indels by local assembly

LDscaff: LD-based scaffolding of de novo genome assemblies

Reconstructing complex regions of genomes using long-read sequencing technology

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

Klumpy: A Tool to Evaluate the Integrity of Long-Read Genome Assemblies and Illusive Sequence Motifs

DEGAP: Dynamic elongation of a genome assembly path

Segmental Duplication As Potential Biomarkers for Non-Invasive Prenatal Testing of Aneuploidies.

Sprites: Detection Of Deletions From Sequencing Data By Re-Aligning Split Reads