Revisiting Waiting Times in DNA evolution

TL;DR

This paper investigates the waiting times for specific DNA motifs to appear under evolution, comparing existing models and introducing a new approach using clump analysis and automata to improve accuracy, especially for overlapping words.

Contribution

It introduces a novel clump analysis and automaton-based method to better estimate waiting times for overlapping DNA motifs, extending previous models.

Findings

Automata approach confirms previous models for non-overlapping words.

Provides up to 44% correction for highly overlapping words.

Proves quasi-linear behavior of motif appearance probability over large ranges.

Abstract

Transcription factors are short stretches of DNA (or $k$-mers) mainly located in promoters sequences that enhance or repress gene expression. With respect to an initial distribution of letters on the DNA alphabet, Behrens and Vingron consider a random sequence of length $n$ that does not contain a given $k$-mer or word of size $k$. Under an evolution model of the DNA, they compute the probability $\mathfrak{p}_n$ that this $k$-mer appears after a unit time of 20 years. They prove that the waiting time for the first apparition of the $k$-mer is well approximated by $T_n=1/\mathfrak{p}_n$. Their work relies on the simplifying assumption that the $k$-mer is not self-overlapping. They observe in particular that the waiting time is mostly driven by the initial distribution of letters. Behrens et al. use an approach by automata that relaxes the assumption related to words overlaps. Their numerical evaluations confirms the validity of Behrens and Vingron approach for non self-overlapping words, but provides up to 44% corrections for highly self-overlapping words such as $\mathtt{AAAAA}$. We devised an approach of the problem by clump analysis and generating functions; this approach leads to prove a quasi-linear behaviour of $\mathfrak{p}_n$ for a large range of values of $n$, an important result for DNA evolution. We present here this clump analysis, first by language decomposition, and next by an automaton construction; finally, we describe an equivalent approach by construction of Markov automata.