Rapidly Registering Identity-by-Descent Across Ancestral Recombination Graphs.
The genomes of remotely related individuals occasionally contain long segments that are identical by descent (IBD). Sharing of IBD segments has many applications in population and medical genetics, and it is thus desirable to study their properties in simulations. However, no current method provides a direct, efficient means to extract IBD segments from simulated genealogies. Here, we introduce computationally efficient approaches to extract ground-truth IBD segments from a sequence of genealogies, or equivalently, an ancestral recombination graph. Specifically, we use a two-step scheme, where we first identify putative shared segments by comparing the common ancestors of all pairs of individuals at some distance apart. This reduces the search space considerably, and we then proceed by determining the true IBD status of the candidate segments. Under some assumptions and when allowing a limited resolution of segment lengths, our run-time complexity is reduced from O(n(3) log n) for the naïve algorithm to O(n log n), where n is the number of individuals in the sample.
1 Department of Computer Science, Columbia University , New York, New York.
Shuo
S
Yang
3 Braun School of Public Health, Faculty of Medicine, Hebrew University, Jerusalem, Israel .
Shai
S
Carmi
1 Department of Computer Science, Columbia University , New York, New York.|2 Department of Systems Biology, Columbia University , New York, New York.
Itsik
I
Pe'er
D000465:Algorithms; D019295:Computational Biology / Q000379:methods*; D005828:Genetics, Population; D016678:Genome
ancestral recombination graphs; identity-by-descent; population genetics; simulation