Improved gene tree inference from removing alignment errors both from focal genes and when training substitution models

Abstract

Multiple Sequence Alignment (MSA) is a key step in phylogenetic analysis, and is prone to error. Unfortunately, algorithms that remove likely alignment errors from MSAs sometimes also remove informative residues, making phylogenetic tree inference worse. Here we present a novel MSA cleaning algorithm based on consensus between MSAs using a range of Hidden Markov Models and guide trees, named CLOAK (CLeaning On the basis of Alignment C(K)onsensus). CLOAK is a gentle filter, with a low false positive rate for removal from MSAs according to the BALiBASE benchmarks, while still removing a significant fraction of likely alignment errors. Gentle vs. stringent MSA filtering methods are appropriate for different tasks. We assess methods based on their ability to bring the gene trees of single copy orthologs closer to the accepted species tree. Amino acid substitution models trained on filtered MSAs improve gene tree inference, with stricter filtering methods providing the biggest model improvements. In contrast, it is gentler filtering of single gene MSAs that provides additional improvements to gene tree inference, with CLOAK performing best.