Theoretical and Practical Considerations when using Retroelement Insertions to Estimate Species Trees in the Anomaly Zone
- PMID: 34677617
- DOI: 10.1093/sysbio/syab086
Theoretical and Practical Considerations when using Retroelement Insertions to Estimate Species Trees in the Anomaly Zone
Abstract
A potential shortcoming of concatenation methods for species tree estimation is their failure to account for incomplete lineage sorting. Coalescent methods address this problem but make various assumptions that, if violated, can result in worse performance than concatenation. Given the challenges of analyzing DNA sequences with both concatenation and coalescent methods, retroelement insertions (RIs) have emerged as powerful phylogenomic markers for species tree estimation. Here, we show that two recently proposed quartet-based methods, SDPquartets and ASTRAL_BP, are statistically consistent estimators of the unrooted species tree topology under the coalescent when RIs follow a neutral infinite-sites model of mutation and the expected number of new RIs per generation is constant across the species tree. The accuracy of these (and other) methods for inferring species trees from RIs has yet to be assessed on simulated data sets, where the true species tree topology is known. Therefore, we evaluated eight methods given RIs simulated from four model species trees, all of which have short branches and at least three of which are in the anomaly zone. In our simulation study, ASTRAL_BP and SDPquartets always recovered the correct species tree topology when given a sufficiently large number of RIs, as predicted. A distance-based method (ASTRID_BP) and Dollo parsimony also performed well in recovering the species tree topology. In contrast, unordered, polymorphism, and Camin-Sokal parsimony (as well as an approach based on MDC) typically fail to recover the correct species tree topology in anomaly zone situations with more than four ingroup taxa. Of the methods studied, only ASTRAL_BP automatically estimates internal branch lengths (in coalescent units) and support values (i.e., local posterior probabilities). We examined the accuracy of branch length estimation, finding that estimated lengths were accurate for short branches but upwardly biased otherwise. This led us to derive the maximum likelihood (branch length) estimate for when RIs are given as input instead of binary gene trees; this corrected formula produced accurate estimates of branch lengths in our simulation study provided that a sufficiently large number of RIs were given as input. Lastly, we evaluated the impact of data quantity on species tree estimation by repeating the above experiments with input sizes varying from 100 to 100,000 parsimony-informative RIs. We found that, when given just 1000 parsimony-informative RIs as input, ASTRAL_BP successfully reconstructed major clades (i.e., clades separated by branches $>0.3$ coalescent units) with high support and identified rapid radiations (i.e., shorter connected branches), although not their precise branching order. The local posterior probability was effective for controlling false positive branches in these scenarios. [Coalescence; incomplete lineage sorting; Laurasiatheria; Palaeognathae; parsimony; polymorphism parsimony; retroelement insertions; species trees; transposon.].
© The Author(s) 2021. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.
Similar articles
-
ILS-Aware Analysis of Low-Homoplasy Retroelement Insertions: Inference of Species Trees and Introgression Using Quartets.J Hered. 2020 Apr 2;111(2):147-168. doi: 10.1093/jhered/esz076. J Hered. 2020. PMID: 31837265
-
Whole-Genome Analyses Resolve the Phylogeny of Flightless Birds (Palaeognathae) in the Presence of an Empirical Anomaly Zone.Syst Biol. 2019 Nov 1;68(6):937-955. doi: 10.1093/sysbio/syz019. Syst Biol. 2019. PMID: 31135914 Free PMC article.
-
Collapsing dubiously resolved gene-tree branches in phylogenomic coalescent analyses.Mol Phylogenet Evol. 2021 May;158:107092. doi: 10.1016/j.ympev.2021.107092. Epub 2021 Feb 2. Mol Phylogenet Evol. 2021. PMID: 33545272
-
Challenges in Species Tree Estimation Under the Multispecies Coalescent Model.Genetics. 2016 Dec;204(4):1353-1368. doi: 10.1534/genetics.116.190173. Genetics. 2016. PMID: 27927902 Free PMC article. Review.
-
Delimiting Coalescence Genes (C-Genes) in Phylogenomic Data Sets.Genes (Basel). 2018 Feb 26;9(3):123. doi: 10.3390/genes9030123. Genes (Basel). 2018. PMID: 29495400 Free PMC article. Review.
Cited by
-
Quartets enable statistically consistent estimation of cell lineage trees under an unbiased error and missingness model.Algorithms Mol Biol. 2023 Dec 1;18(1):19. doi: 10.1186/s13015-023-00248-w. Algorithms Mol Biol. 2023. PMID: 38041123 Free PMC article.
-
New Insights Into the Relationships Within Subtribe Scorzonerinae (Cichorieae, Asteraceae) Using Hybrid Capture Phylogenomics (Hyb-Seq).Front Plant Sci. 2022 Jul 1;13:851716. doi: 10.3389/fpls.2022.851716. eCollection 2022. Front Plant Sci. 2022. PMID: 35873957 Free PMC article.
-
Homoplasy of Retrotransposon Insertions in Toothed Whales.Genes (Basel). 2023 Sep 21;14(9):1830. doi: 10.3390/genes14091830. Genes (Basel). 2023. PMID: 37761970 Free PMC article.
-
Dollo-CDP: a polynomial-time algorithm for the clade-constrained large Dollo parsimony problem.Algorithms Mol Biol. 2024 Jan 8;19(1):2. doi: 10.1186/s13015-023-00249-9. Algorithms Mol Biol. 2024. PMID: 38191515 Free PMC article.
-
SINE-Based Phylogenomics Reveal Extensive Introgression and Incomplete Lineage Sorting in Myotis.Genes (Basel). 2022 Feb 23;13(3):399. doi: 10.3390/genes13030399. Genes (Basel). 2022. PMID: 35327953 Free PMC article.
Publication types
MeSH terms
Substances
Associated data
LinkOut - more resources
Full Text Sources
Other Literature Sources