title: Evidence of Statistical Inconsistency of Phylogenetic Methods in the Presence of Multiple Sequence Alignment Uncertainty.
abstract: Evolutionary studies usually use a two-step process to investigate sequence data. Step one estimates a multiple sequence alignment (MSA) and step two applies phylogenetic methods to ask evolutionary questions of that MSA. Modern phylogenetic methods infer evolutionary parameters using maximum likelihood or Bayesian inference, mediated by a probabilistic substitution model that describes sequence change over a tree. The statistical properties of these methods mean that more data directly translates to an increased confidence in downstream results, providing the substitution model is adequate and the MSA is correct. Many studies have investigated the robustness of phylogenetic methods in the presence of substitution model misspecification, but few have examined the statistical properties of those methods when the MSA is unknown. This simulation study examines the statistical properties of the complete two-step process when inferring sequence divergence and the phylogenetic tree topology. Both nucleotide and amino acid analyses are negatively affected by the alignment step, both through inaccurate guide tree estimates and through overfitting to that guide tree. For many alignment tools these effects become more pronounced when additional sequences are added to the analysis. Nucleotide sequences are particularly susceptible, with MSA errors leading to statistical support for long-branch attraction artifacts, which are usually associated with gross substitution model misspecification. Amino acid MSAs are more robust, but do tend to arbitrarily resolve multifurcations in favor of the guide tree. No inference strategies produce consistently accurate estimates of divergence between sequences, although amino acid MSAs are again more accurate than their nucleotide counterparts. We conclude with some practical suggestions about how to limit the effect of MSA uncertainty on evolutionary inference.
authors.affiliation: Faculty of Life Sciences, University of Manchester, United Kingdom Evolutionary Biology, Evolutionary Biology Centre, Uppsala University, Sweden.
authors.forename: A S
authors.initials: AS
authors.lastname: Md Mukarram Hossain
authors.affiliation: Faculty of Life Sciences, University of Manchester, United Kingdom.
authors.forename: Benjamin P
authors.initials: BP
authors.lastname: Blackburne
authors.affiliation: Faculty of Life Sciences, University of Manchester, United Kingdom.
authors.forename: Abhijeet
authors.initials: A
authors.lastname: Shah
authors.affiliation: Evolutionary Biology, Evolutionary Biology Centre, Uppsala University, Sweden simon.whelan@ebc.uu.se.
authors.forename: Simon
authors.initials: S
authors.lastname: Whelan
mesh_terms: D016477:Artifacts; D015233:Models, Statistical; D010802:Phylogeny; D016415:Sequence Alignment / Q000379:methods*; D035501:Uncertainty
keywords: maximum likelihood; multiple sequence alignment; phylogenetic inference; sequence divergence; tree inference