The Power of Alignment-Free Histogram-based Functions: a Comprehensive Genome Scale Experimental Analysis – Version 1

Abstract

Motivation: Alignment-free (AF, for short) distance/similarity functions are a key tool for sequence analysis. Experimental studies on real datasets abound and, to some extent, there are also studies regarding their control of false positive rate (Type I error). However, assessment of their power, i.e., their ability to identify true similarity, has been limited to some members or variants of the $D_2$ family and the asymptotic theoretic results have been complemented by experimental studies on short sequences, not adequate for current genome-scale applications. Such a State of the Art is methodologically problematic, since information regarding a key feature such as power is either missing or limited.

Results: By concentrating on histogram-based AF functions, we perform the first coherent and uniform evaluation of the power of those functions, involving also Type I error for completeness. The experiments carried out are extensive, as we use two Alternative models of important genomic features (CIS Regulatory Modules and Horizontal Gene Transfer), sequence lengths from a few thousand to millions and different values of $k$. As a result, and using power, we provide a characterization of those AF functions that is novel and informative. Indeed, we identify weak and strong points of each function considered, which may be used as a guide to choose one for analysis tasks. In synthesis, and quite remarkably, of the fifteen functions that we have considered, only four stand out. Finally, in order to encourage the use of our methodology for validation of future AF functions, the Big Data platform supporting it is public.

Avatar
Raffaele Giancarlo
Full professor

Professor of Algorithms and PI of the project: Research Unit Univ. Palermo