Alignment Modulates Ancestral Sequence Reconstruction Accuracy

Abstract

Accurate reconstruction of ancestral states is a critical evolutionary analysis when studying ancient proteins and comparing biochemical properties between parental or extinct species and their extant relatives. It relies on multiple sequence alignment (MSA) which may introduce biases, and it remains unknown how MSA methodological approaches impact ancestral sequence reconstruction (ASR). Here, we investigate how MSA methodology modulates ASR using a simulation study of various evolutionary scenarios. We evaluate the accuracy of ancestral protein sequence reconstruction for simulated data and compare reconstruction outcomes using different alignment methods. Our results reveal biases introduced not only by aligner algorithms and assumptions, but also tree topology and the rate of insertions and deletions. Under many conditions we find no substantial differences between the MSAs. However, increasing the difficulty for the aligners can significantly impact ASR. The MAFFT consistency aligners and PRANK variants exhibit the best performance, whereas FSA displays limited performance. We also discover a bias towards reconstructed sequences longer than the true ancestors, deriving from a preference for inferring insertions, in almost all MSA methodological approaches. In addition, we find measures of MSA quality generally correlate highly with reconstruction accuracy. Thus, we show MSA methodological differences can affect the quality of reconstructions and propose MSA methods should be selected with care to accurately determine ancestral states with confidence.

Evolution of Predicted Acid Resistance Mechanisms in the Extremely Acidophilic Leptospirillum Genus

“…Gene duplication and gene diversification events were identified using a combination of phylogenetic inference based on alignments of families (cluster analysis) [123] and calculations of D n /D s [94][95][96]. Multiple examples of gene duplication events were detected.…”

Non-conservation of folding rates in the thioredoxin family reveals degradation of ancestral unassisted-folding

“…There could be some doubt, therefore, that the S74G replacement actually occurred in the line of descent that led to E. coli thioredoxin. This is of particular interest given discussions from the literature (Eick et al, 2017;Groussin et al, 2015;Vialle et al, 2018;Williams et al, 2006) that ancestral sequence reconstruction can potentially be biased from uncertainties in the process. In this case, however, the identification of serine as the ancestral state at position 74 is quite robust.…”

Benchmarking tree and ancestral sequence inference for B cell receptor sequences

“…In addition, indels are not treated as the informative characters they are in mainstream phylogenetics software; rather, they are treated as missing data. Benchmarking phylogenetic tools would also require benchmarking the alignment step, which has an effect on ancestral sequence reconstruction accuracy (73). Nevertheless, this will be another important focus for future tool development and ancestral sequence reconstruction benchmarking within the field of BCR phylogenetic reconstruction.…”