CATE: A fast and scalable CUDA implementation to conduct highly parallelized evolutionary tests on large scale genomic data.

Abstract

Statistical tests for molecular evolution provide quantifiable insights into the selection pressures that govern a genomes evolution. Increasing sample sizes used for analysis leads to higher statistical power. However, this requires more computational nodes or longer computational time. CATE (CUDA Accelerated Testing of Evolution) is a computational solution to this problem comprised of two main innovations. The first is a file organization system coupled with a novel search algorithm and the second is a large-scale parallelization of algorithms using both GPU and CPU. CATE is capable of conducting evolutionary tests such as Tajimas D, Fu and Lis, and Fay and Wus test statistics, McDonald-Kreitman Neutrality Index, Fixation Index, and Extended Haplotype Homozygosity. CATE is magnitudes faster than standard tools with benchmarks estimating it being on average over 180 times faster. For instance, CATE processes all 54,849 human genes for all 22 autosomal chromosomes across the five super populations present in the 1000 Genomes Project in less than thirty minutes while counterpart software took 3.62 days. This proven framework has the potential to be adapted for GPU-accelerated large-scale parallel analyses of many evolutionary and genomic analyses.