Statistical framework to determine indel-length distribution

Elya Wygoda, Gil Loewenthal, Asher Moshe, Michael Alburquerque, Itay Mayrose, Tal Pupko*

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

Abstract

Motivation: Insertions and deletions (indels) of short DNA segments, along with substitutions, are the most frequent molecular evolutionary events. Indels were shown to affect numerous macro-evolutionary processes. Because indels may span multiple positions, their impact is a product of both their rate and their length distribution. An accurate inference of indel-length distribution is important for multiple evolutionary and bioinformatics applications, most notably for alignment software. Previous studies counted the number of continuous gap characters in alignments to determine the best-fitting length distribution. However, gap-counting methods are not statistically rigorous, as gap blocks are not synonymous with indels. Furthermore, such methods rely on alignments that regularly contain errors and are biased due to the assumption of alignment methods that indels lengths follow a geometric distribution. Results: We aimed to determine which indel-length distribution best characterizes alignments using statistical rigorous methodologies. To this end, we reduced the alignment bias using a machine-learning algorithm and applied an Approximate Bayesian Computation methodology for model selection. Moreover, we developed a novel method to test if current indel models provide an adequate representation of the evolutionary process. We found that the best-fitting model varies among alignments, with a Zipf length distribution fitting the vast majority of them.

Original languageEnglish
Article numberbtae043
JournalBioinformatics
Volume40
Issue number2
DOIs
StatePublished - 1 Feb 2024

Funding

FundersFunder number
Edmond J. Safra Center for Bioinformatics
Israel Science Foundation2818/21
Tel Aviv University

    Fingerprint

    Dive into the research topics of 'Statistical framework to determine indel-length distribution'. Together they form a unique fingerprint.

    Cite this