Starless bias and parameter-estimation bias in the likelihood-based phylogenetic method

Abstract

I analyzed various site pattern combinations in a 4-OTU case to identify sources of starless bias and parameter-estimation bias in likelihood-based phylogenetic methods, and reported three significant contributions. First, the likelihood method is counterintuitive in that it may not generate a star tree with sequences that are equidistant from each other. This behaviour, dubbed starless bias, happens in a 4-OTU tree when there is an excess (i.e., more than expected from a star tree and a substitution model) of conflicting phylogenetic signals supporting the three resolved topologies equally. Special site pattern combinations leading to rejection of a star tree, when sequences are equidistant from each other, were identified. Second, fitting gamma distribution to model rate heterogeneity over sites is strongly confounded with tree topology, especially in conjunction with the starless bias. I present examples to show dramatic differences in the estimated shape parameter Α between a star tree and a resolved tree. There may be no rate heterogeneity over sites (with the estimated Α > 10000) when a star tree is imposed, but Α < 1 (suggesting strong rate heterogeneity over sites) when an (incorrect) resolved tree is imposed. Thus, the dependence of “rate heterogeneity” on tree topology implies that “rate heterogeneity” is not a sequence-specific feature, cautioning against interpreting a small Α to mean that some sites are under strong purifying selection and others not. Thirdly, because there is no existing (and working) likelihood method for evaluating a star tree with continuous gamma-distributed rate, I have implemented the method for JC69 in a self-contained R script for a four-OTU tree (star or resolved), in addition to another R script assuming a constant rate over sites. These R scripts should be useful for teaching and exploring likelihood methods in phylogenetics.

Keywords

maximum likelihood - molecular phylogenetics - rate heterogeneity - starless - star-tree paradox

Publication History

Received: 17 September 2018

Accepted: 03 April 2019

Article published online:
10 May 2021

• References

• 1 Saitou N, Nei M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol Biol Evol 1987; 4: 406-425
• 2 Desper R, Gascuel O. Theoretical foundation of the balanced minimum evolution method of phylogenetic inference and its relationship to weighted least-squares tree fitting. Mol Biol Evol 2004; 21: 587-598
• 3 Xia X. Phylogenetic bias in the likelihood method caused by missing data coupled with Among-Site rate variation: An analytical approach. Basu M, Pan Y, Wang J. Bioinformatics Research and Applications. Springer; 2014: 012-023
• 4 Jukes TH, Cantor CR. Evolution of protein molecules. Munro HN. Mammalian Protein Metabolism. New York: Academic Press; 1969: 021-123
• 5 Hasegawa M, Kishino H. Heterogeneity of tempo and mode of mitochondrial DNA evolution among mammalian orders. Jpn J Genet 1989; 64: 243-258
• 6 Kishino H, Hasegawa M. Evaluation of the maximum likelihood estimate of the evolutionary tree topologies from DNA sequence data, and the branching order in Hominoidea. J Mol Evol 1989; 29: 170-179
• 7 Hasegawa M, Kishino H, Yano T. Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol 1985; 22: 160-174
• 8 Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol 1993; 10: 512-526
• 9 Lanave C, Preparata G, Saccone C. et al. A new method for calculating evolutionary substitution rates. J Mol Evol 1984; 20: 086-093
• 10 Tavarï S. Some probabilistic and statistical problems in the analysis of DNA sequences. Miura RM. Lectures on Mathematics in the Life Sciences. Providence, RI: Amer Math Soc; 1986: 057-086
• 11 Xia X. Deriving transition probabilities and evolutionary distances from substitution rate matrix by probability reasoning. J Genet Genome Res 2017; 3: 031
• 12 Xia X. Nucleotide substitution models and evolutionary distances. Bioinf Cell 2018; 269-314
• 13 Yang Z. PAML 4: Phylogenetic analysis by maximum likelihood. Mol Biol Evol 2007; 24: 1586-1591
• 14 Xia X. Maximum likelihood in molecular phylogenetics. Bioinf Cell 2018; 381-395
• 15 Guindon S, Dufayard JF, Lefort V. et al. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst Biol 2010; 59: 307-321
• 16 Xia X. DAMBE7: New and improved tools for data analysis in molecular biology and evolution. Mol Biol Evol 2018; 35: 1550-1552