A maximum pseudo-likelihood approach for estimating species trees under the coalescent model
Open accessAbstract
BACKGROUND: Several phylogenetic approaches have been developed to estimate species trees from collections of gene trees. However, maximum likelihood approaches for estimating species trees under the coalescent model are limited. Although the likelihood of a species tree under the multispecies coalescent model has already been derived by Rannala and Yang, it can be shown that the maximum likelihood estimate (MLE) of the species tree (topology, branch lengths, and population sizes) from gene trees under this formula does not exist. In this paper, we develop a pseudo-likelihood function of the species tree to obtain maximum pseudo-likelihood estimates (MPE) of species trees, with branch lengths of the species tree in coalescent units. RESULTS: We show that the MPE of the species tree is statistically consistent as the number M of genes goes to infinity. In addition, the probability that the MPE of the species tree matches the true species tree converges to 1 at rate O(M -1). The simulation results confirm that the maximum pseudo-likelihood approach is statistically consistent even when the species tree is in the anomaly zone. We applied our method, Maximum Pseudo-likelihood for Estimating Species Trees (MP-EST) to a mammal dataset. The four major clades found in the MP-EST tree are consistent with those in the Bayesian concatenation tree. The bootstrap supports for the species tree estimated by the MP-EST method are more reasonable than the posterior probability supports given by the Bayesian concatenation method in reflecting the level of uncertainty in gene trees and controversies over the relationship of four major groups of placental mammals. CONCLUSIONS: MP-EST can consistently estimate the topology and branch lengths (in coalescent units) of the species tree. Although the pseudo-likelihood is derived from coalescent theory, and assumes no gene flow or horizontal gene transfer (HGT), the MP-EST method is robust to a small amount of HGT in the dataset. In addition, increasing the number of genes does not increase the computational time substantially. The MP-EST method is fast for analyzing datasets that involve a large number of genes but a moderate number of species.
The claims and statistics of this paper haven’t been checked yet.
Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.
Integrity signals
Registry
Cites
Statistics
Claims the paper makes
The claims and statistics of this paper haven’t been checked yet.
Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.
Claims other papers make about this one
Each is checked against this paper’s own text: does it say what the citing paper says it does?
Comparative genomics analyses (species-tree inference, comparative annotation, selection detection) require genome alignments.
“Comparative genomics analyses, including species-tree inference7,8, comparative annotation9,10, and selection detection11,12, require genome alignments.”
From Progressive Cactus is a multiple-genome aligner for the thousand-genome era (2020) · cites this paper as [7]
The Methods section (page not available) describes estimating gene trees from multilocus sequences and resampling sites in aligned sequences, supporting the use of alignments for species-tree inference. But it says MP-EST constructs a species tree from rooted gene trees; it does not establish that genome alignments are required. This citation supports the species-tree portion of the sentence, not the separately cited annotation and selection examples. (LLM, conf=0.85)
Experimental — accuracy not yet measured
Self-check
Percentages, sums and reported differences re-computed from the paper’s own numbers.
The self-check needs the paper’s full text.
Cites
Evolutionary trees from DNA sequences: A maximum likelihood approach.
Felsenstein, J. (1981) · Journal of Molecular Evolution · doi:10.1007/bf01734359 · Cited by 14,571 · see citing works
Genome-scale approaches to resolving incongruence in molecular phylogenies.
Rokas, A., Williams, B. L., King, N., & Carroll, S. B. (2003) · Nature · doi:10.1038/nature02053 · Cited by 1,512 · see citing works
Nonparametric estimates of standard error: The jackknife, the bootstrap and other methods
EFRON, B. (1981) · Biometrika · doi:10.1093/biomet/68.3.589 · Cited by 1,399 · see citing works
Flagged first, then by citations
Research works printed but not linked: 62
- Jennings WB Edwards SV Speciational history of Australian grass finches (Poephila) inferred from thirty gene trees Evolution 2005 59 2033 2047 16261740
- Pollard DA Iyer VN Moses AM Eisen MB Widespread discordance of gene trees with species tree in Drosophila: Evidence for incomplete lineage sorting Plos Genet 2006 2 1634 1647 PMC1626107 17132051
- Brumfield RT Jernigan RW McDonald DB Braun MJ Evolutionary implications of divergent clines in an avian (Manacus: Aves) hybrid zone Evolution 2001 55 2070 2087 11761066
- Huelsenbeck JP Bull JJ Cunningham CW Combining data in phylogenetic analysis Trends Ecol Evol 1996 11 152 158 10.1016/0169-5347(96)10006-9 21237790 · DOI
- de Queiroz A Gatesy J The supermatrix approach to systematics Trends Ecol Evol 2007 22 34 41 10.1016/j.tree.2006.10.002 17046100 · DOI
- Mossel E Roch S Incomplete lineage sorting: consistent phylogeny estimation from multiple Loci IEEE/ACM Transactions on Computational Biology and Bioinformatics 2007 7 166 171 10.1109/TCBB.2008.66 20150678 · DOI
- Meng C Kubatko LS Detecting hybrid speciation in the presence of incomplete lineage sorting using gene tree incongruence: a model Theor Popul Biol 2009 75 35 45 10.1016/j.tpb.2008.10.004 19038278 · DOI
- Maddison WP Knowles LL Inferring phylogeny despite incomplete lineage sorting Syst Biol 2006 55 21 30 10.1080/10635150500354928 16507521 · DOI
- Liu L BEST: Bayesian estimation of species trees under the coalescent model Bioinformatics 2008 24 2542 2543 10.1093/bioinformatics/btn484 18799483 · DOI
- Liu L Pearl DK Species trees from gene trees: Reconstructing Bayesian posterior distributions of a species phylogeny using estimated gene tree distributions Syst Biol 2007 56 504 514 10.1080/10635150701429982 17562474 · DOI
- Liu L Yu L Pearl DK Edwards SV Estimating species phylogenies using coalescence times among sequences Syst Biol 2009 58 468 477 10.1093/sysbio/syp031 20525601 · DOI
- Ewing GB Ebersberger I Schmidt HA von Haeseler A Rooted triple consensus and anomalous gene trees BMC Evol Biol 2008 8 118 10.1186/1471-2148-8-118 18439266 PMC2409437 · DOI
- Bryant D Berry V A structured family of clustering and tree construction methods Adv Appl Math 2001 27 705 732 10.1006/aama.2001.0758 · DOI
- Day WHE McMorris FR Wilkinson M Explosions and hot spots in supertree methods J Theor Biol 2008 253 345 348 10.1016/j.jtbi.2008.03.024 18472112 · DOI
- Bininda-Emonds ORP Bryant HN Properties of matrix representation with parsimony analyses Syst Biol 1998 47 497 508 12066690
- Cotton JA Wilkinson M Majority-rule supertrees Syst Biol 2007 56 445 452 10.1080/10635150701416682 17558966 · DOI
- Page RDM Extracting species trees from complex gene trees: Reconciled trees and vertebrate phylogeny Mol Phylogenet Evol 2000 14 89 106 10.1006/mpev.1999.0676 10631044 · DOI
- Bonizzoni P Della Vedova G Dondi R Reconciling gene trees to a species tree 2003 Berlin: Springer
- Page RDM Charleston MA From gene to organismal phylogeny: Reconciled trees and the gene tree species tree problem Mol Phylogenet Evol 1997 7 231 240 10.1006/mpev.1996.0390 9126565 · DOI
- Slowinski JB Knight A Rooney AP Inferring species trees from gene trees: A phylogenetic analysis of the elapidae (Serpentes) based on the amino acid sequences of venom proteins Mol Phylogenet Evol 1997 8 349 362 10.1006/mpev.1997.0434 9417893 · DOI
- Steel M Rodrigo A Maximum likelihood supertrees Syst Biol 2008 57 243 250 10.1080/10635150802033014 18398769 · DOI
- Ane C Larget B Baum DA Smith SD Rokas A Bayesian estimation of concordance among gene trees Mol Biol Evol 2007 24 412 426 10.1093/molbev/msl170 17095535 · DOI
- Rannala B Yang ZH Bayes estimation of species divergence times and ancestral population sizes using DNA sequences from multiple loci Genetics 2003 164 1645 1656 12930768 PMC1462670
- Nielsen R Maximum likelihood estimation of population divergence times and population phylogenies under the infinite sites model Theor Popul Biol 1998 53 143 151 10.1006/tpbi.1997.1348 9615473 · DOI
- Degnan JH Salter LA Gene tree distributions under the coalescent process Evolution 2005 59 24 37 15792224
- Takahata N Gene genealogy in three related populations: consistency probability between gene and population trees Genetics 1989 122 957 966 2759432 PMC1203770
- Degnan JH Rosenberg NA Gene tree discordance, phylogenetic inference, and the multispecies coalescent Trends Ecol Evol 2009 24 332 340 10.1016/j.tree.2009.01.009 19307040 · DOI
- Efromovich S Kubatko LS Coalescent time distributions in trees of arbitrary size Stat Appl Genet Mol Biol 2008 7 Article 2 18241192
- Carstens BC Knowles LL Estimating species phylogeny from gene-tree probabilities despite incomplete lineage sorting: an example from Melanoplus grasshoppers Syst Biol 2007 56 400 411 10.1080/10635150701405560 17520504 · DOI
- Liu L Yu L Pearl DK Maximum Tree: a consistent estimator of the species tree J Math Biol 2010 60 95 106 10.1007/s00285-009-0260-0 19283383 · DOI
- Kubatko LS Carstens BC Knowles LL STEM: species tree estimation using maximum likelihood for gene trees under coalescence Bioinformatics 2009 25 971 973 10.1093/bioinformatics/btp079 19211573 · DOI
- Edwards SV Is a new and general theory of molecular systematics emerging? Evolution 2009 63 1 19 10.1111/j.1558-5646.2008.00549.x 19146594 · DOI · Paper page
- Liu L Pearl DK Brumfield RT Edwards SV Estimating species trees using multiple-allele DNA sequence data Evolution 2008 62 2080 2091 10.1111/j.1558-5646.2008.00414.x 18462214 · DOI
- Degnan JH Rosenberg NA Discordance of species trees with their most likely gene trees Plos Genetics 2006 2 762 768 10.1371/journal.pgen.0020068 PMC1464820 16733550 · DOI
- Steel M The Complexity of Reconstructing Trees from Qualitative Characters and Subtrees Journal of Classification 1992 9 91 116 10.1007/BF02618470 · DOI
- Degnan JH DeGiorgio M Bryant D Rosenberg NA Properties of consensus methods for inferring species trees from gene trees Syst Biol 2009 58 35 54 10.1093/sysbio/syp008 20525567 PMC2909780 · DOI
- Pamilo P Nei M Relationships between Gene Trees and Species Trees Mol Biol Evol 1988 5 568 583 3193878
- Wang J A pseudo-likelihood method for estimating effective population size from temporally spaced samples Genet Res 2001 78 243 257 10.1017/S0016672301005286 11865714 · DOI
- Besag J Statistical analysis of non-lattice data The Statistician 1975 24 179 195 10.2307/2987782 · DOI · Paper page
- Feller W An Introduction to Probability Theory and Its Applications New York: Wiley 1968 1
- Seo T-K Calculating Bootstrap Probabilities of Phylogeny Using Multilocus Sequence Data Mol Biol Evol 2008 25 960 971 10.1093/molbev/msn043 18281270 · DOI
- Soltis PS Soltis DE Applying the bootstrap in phylogeny reconstruction Stat Sci 2003 18 256 267 10.1214/ss/1063994980 · DOI
- Liu L Yu L Phybase: an R package for species tree analysis Bioinformatics 2010 20156990
- Casella G Berger RL Statistical Inference 2002 Duxbury Press
- Robinson DR Foulds LR Comparison of phylogenetic trees Mathematical Biosciences 1981 53 131 147 10.1016/0025-5564(81)90043-2 · DOI
- Guindon S Gascuel O A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood Syst Biol 2003 52 696 704 10.1080/10635150390235520 14530136 · DOI
- Driskell AC Ane C Burleigh JG McMahon MM O'Meara BC Sanderson MJ Prospects for building the tree of life from large sequence databases Science 2004 306 1172 1174 10.1126/science.1102036 15539599 · DOI
- Huelsenbeck JP Ronquist F MRBAYES: Bayesian inference of phylogenetic trees Bioinformatics 2001 17 754 755 10.1093/bioinformatics/17.8.754 11524383 · DOI · Paper page
- Springer MS Burk-Herrick A Meredith R Eizirik E Teeling E O'Brien SJ Murphy WJ The adequacy of morphology for reconstructing the early history of placental mammals Syst Biol 2007 56 673 684 10.1080/10635150701491149 17661234 · DOI
- Sanderson MJ Shaffer HB Troubleshooting phylogenies Annu Rev Ecol Syst 2002 49 72 10.1146/annurev.ecolsys.33.010802.150509 · DOI
- Cummings MP Handley SA Myers DS Reed DL Rokas A Winka K Comparing bootstrap and posterior probability values in the four-taxon case Syst Biol 2003 52 477 487 10.1080/10635150390218213 12857639 · DOI
- Douady CJ Delsuc F Boucher Y Doolittle WF Douzery EJ Comparison of Bayesian and maximum likelihood bootstrap measures of phylogenetic reliability Mol Biol Evol 2003 20 248 254 10.1093/molbev/msg042 12598692 · DOI
- Hallstrom BM Kullberg M Nilsson MA Janke A Phylogenomic data analyses provide evidence that Xenarthra and Afrotheria are sister groups Mol Biol Evol 2007 24 2059 2068 10.1093/molbev/msm136 17630282 · DOI
- Shoshani J McKenna MC Higher taxonomic relationships among extant mammals based on morphology, with selected comparisons of results from molecular data Mol Phylogenet Evol 1998 9 572 584 10.1006/mpev.1998.0520 9668007 · DOI
- Kriegs JO Churakov G Kiefmann M Jordan U Brosius J Schmitz J Retroposed elements as archives for the evolutionary history of placental mammals PLoS Biol 2006 4 e91 10.1371/journal.pbio.0040091 16515367 PMC1395351 · DOI
- Kjer KM Honeycutt RL Site specific rates of mitochondrial genomes and the phylogeny of eutheria BMC Evol Biol 2007 7 8 10.1186/1471-2148-7-8 17254354 PMC1796853 · DOI
- Nikolaev S Montoya-Burgos JI Margulies EH Rougemont J Nyffeler B Antonarakis SE Early history of mammals is elucidated with the ENCODE multiple species sequencing data PLoS Genet 2007 3 e2 10.1371/journal.pgen.0030002 17206863 PMC1761045 · DOI
- Churakov G Kriegs JO Baertsch R Zemann A Brosius J Schmitz J Mosaic retroposon insertion patterns in placental mammals Genome Res 2009 19 868 875 10.1101/gr.090647.108 19261842 PMC2675975 · DOI
- Hallstrom BM Janke A Resolution among major placental mammal interordinal relationships with genome data imply that speciation influenced their earliest radiations BMC Evol Biol 2008 8 162 10.1186/1471-2148-8-162 18505555 PMC2435553 · DOI
- Nishihara H Maruyama S Okada N Retroposon analysis and recent geological data suggest near-simultaneous divergence of the three superorders of mammals Proc Natl Acad Sci USA 2009 106 5235 5240 10.1073/pnas.0809297106 19286970 PMC2655268 · DOI
- Eckert AJ Carstens BC Does gene flow destroy phylogenetic signal? The performance of three methods for estimating species phylogenies in the presence of gene flow Mol Phylogenet Evol 2008 49 832 842 10.1016/j.ympev.2008.09.008 18845264 · DOI
- Liu L Yu L Kubatko LS Pearl DK Edwards SV Coalescent methods for estimating multilocus phylogenetic trees Mol Phylogenet Evol 2009 53 320 328 10.1016/j.ympev.2009.05.033 19501178 · DOI
Cited by
Sorted by citations
Whole-genome analyses resolve early branches in the tree of life of modern birds
Jarvis, E. D., Mirarab, S., Aberer, A. J., Li, B., Houde, P., Li, C., Ho, S. Y. W., Faircloth, B. C., Nabholz, B., Howard, J. T., et al. (2014) · Science · doi:10.1126/science.1253451 · Cited by 1,901 · see citing works
Self-citationA comprehensive phylogeny of birds (Aves) using targeted next-generation DNA sequencing
Prum, R. O., Berv, J. S., Dornburg, A., Field, D. J., Townsend, J. P., Lemmon, E. M., & Lemmon, A. R. (2015) · Nature · doi:10.1038/nature15697 · Cited by 1,668 · see citing works
Progressive Cactus is a multiple-genome aligner for the thousand-genome era
Armstrong, J., Hickey, G., Diekhans, M., Fiddes, I. T., Novak, A. M., Deran, A., Fang, Q., Xie, D., Feng, S., Stiller, J., et al. (2020) · Nature · doi:10.1038/s41586-020-2871-y · Cited by 586 · see citing works
Large-scale ruminant genome sequencing provides insights into their evolution and distinct traits
Chen, L., Qiu, Q., Jiang, Y., Wang, K., Lin, Z., Li, Z., Bibi, F., Yang, Y., Wang, J., Nie, W., et al. (2019) · Science · doi:10.1126/science.aav6202 · Cited by 441 · see citing works
Phylogenomics reveals three sources of adaptive variation during a rapid radiation
Pease, J. B., Haak, D. C., Hahn, M. W., & Moyle, L. C. (2016) · PLOS Biology · doi:10.1371/journal.pbio.1002379 · Cited by 432 · see citing works
Statistical binning enables an accurate coalescent-based estimation of the avian tree
Mirarab, S., Bayzid, M. S., Boussau, B., & Warnow, T. (2014) · Science · doi:10.1126/science.1250463 · Cited by 257 · see citing works
- Authors
- Liu, L., Yu, L., & Edwards, S. V.
- Affiliations
- Delaware State University (US)
- Georgia Southern University (US)
- Harvard University (US)
- Corresponding author
- Liu, Liang
- Venue
- BMC Evolutionary Biology
- Published
- 2010
- Volume
- 10
- Issue
- 1
- Pages
- 302
- ISSN
- 1471-2148 (ISSN-L)
- PMC
- PMC2976751
- OpenAlex
- W2125024684
- Cited by
- 607
- References
- 3
- Citations by year
- 2026: 17 · 2025: 26 · 2024: 21 · 2023: 29 · 2022: 37 · 2021: 48 · 2020: 41 · 2019: 62 · 2018: 60 · 2017: 53 · 2016: 61 · 2015: 55 · 2014: 58 · 2013: 24 · 2012: 11
- FWCI
- 4.73
- Type
- journal article
- Language
- en
- Open access
- gold
- doi.org
- Licence
- cc-by
- Locations
- BMC Evolutionary Biology (publishedVersion, cc-by)
- PubMed Central (publishedVersion, cc-by)
- PubMed (publishedVersion)
- Digital Access to Scholarship at Harvard (DASH) (Harvard University) (submittedVersion, cc-by)
- DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion)
- Funders
- National Science Foundation (0743616)
- Division of Environmental Biology
- Topics
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Phylogenetic Studies
- Life Sciences › Biochemistry, Genetics and Molecular Biology › Genetics › Genetic diversity and population structure
- Physical Sciences › Earth and Planetary Sciences › Paleontology › Evolution and Paleontology Studies
- Keywords
- species tree estimation, multispecies coalescent, gene trees, statistical consistency, placental mammals, horizontal gene transfer, bootstrap support, gene flow
- MeSH terms
- Models, Theoretical, Phylogeny, Likelihood Functions, Gene Transfer, Horizontal
About these checks
Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.
Statistics and arithmetic are recomputed, not judged.
Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.
What counts as a claim ›How we measure confidence ›How we check this ›
Are you an author, or do you see a problem with this page? Request a review