Fast estimation of the difference between two PAM/JTT evolutionary distances in triplets of homologous sequences.
BACKGROUND
The estimation of the difference between two evolutionary distances within a triplet of homologs is a common operation that is used for example to determine which of two sequences is closer to a third one. The most accurate method is currently maximum likelihood over the entire triplet. However, this approach is relatively time consuming.


RESULTS
We show that an alternative estimator, based on pairwise estimates and therefore much faster to compute, has almost the same statistical power as the maximum likelihood estimator. We also provide a numerical approximation for its variance, which could otherwise only be estimated through an expensive re-sampling approach such as bootstrapping. An extensive simulation demonstrates that the approximation delivers precise confidence intervals. To illustrate the possible applications of these results, we show how they improve the detection of asymmetric evolution, and the identification of the closest relative to a given sequence in a group of homologs.


CONCLUSIONS
The results presented in this paper constitute a basis for large-scale protein cross-comparisons of pairwise evolutionary distances.
Dessimoz
Christophe
C
ETH Zurich, Institute of Computational Science, 8092 Zürich, Switzerland. cdessimoz@inf.ethz.ch
Gil
Manuel
M
Schneider
Adrian
A
Gonnet
Gaston H
GH
D000465:Algorithms; D001483:Base Sequence; D002874:Chromosome Mapping / Q000379:methods*; D019143:Evolution, Molecular; D015810:Linkage Disequilibrium / Q000235:genetics*; D008969:Molecular Sequence Data; D011506:Proteins / Q000235:genetics*; D016415:Sequence Alignment / Q000379:methods*; D017422:Sequence Analysis, DNA / Q000379:methods*; D012689:Sequence Homology, Nucleic Acid
D011506:Proteins