RERconverge Update: Runtime Reduction and Analysis Function Overhaul
Discuss this preprint
Start a discussion What are Sciety discussions?Listed in
This article is not in any list yet, why not save it to one of your lists.Abstract
Motivation
Convergent evolution, or the independent acquisition of similar phenotypes in distinct lineages, provides a powerful framework for investigating genomic changes associated with a phenotype. This paper details an update to RERconverge, a powerful R package that tests for associations between gene relative evolutionary rates (RERs) and convergent phenotypes to infer genomic regions associated with traits or selective pressures. We introduce new customizable analysis choices and scalable and efficient algorithms that can process larger genomic datasets, a critical improvement as genomic data become available for more species.
Results
Modifications to core functions in the RERconverge pipeline resulted in an immense speedup (by a factor of up to 28.6). The function that tests for associations between phenotypes and RERs has been expanded to include two new analytical methods for outlier control; we also provide here a summary of the statistical tests users can perform, along with their use cases.
Availability and implementation
The code and walkthrough vignettes for the package are available at https://github.com/nclark-lab/RERconverge .
Contact
Nathan L. Clark nclark@pitt.edu ; Maria Chikina mchikina@pitt.edu
