1
|
Li W, Koshkarov A, Tahiri N. Comparison of phylogenetic trees defined on different but mutually overlapping sets of taxa: A review. Ecol Evol 2024; 14:e70054. [PMID: 39119174 PMCID: PMC11307105 DOI: 10.1002/ece3.70054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2023] [Revised: 07/03/2024] [Accepted: 07/10/2024] [Indexed: 08/10/2024] Open
Abstract
Phylogenetic trees represent the evolutionary relationships and ancestry of various species or groups of organisms. Comparing these trees by measuring the distance between them is essential for applications such as tree clustering and the Tree of Life project. Many distance metrics for phylogenetic trees focus on trees defined on the same set of taxa. However, some problems require calculating distances between trees with different but overlapping sets of taxa. This study reviews state-of-the-art distance measures for such trees, covering six major approaches, including the constraint-based Robinson-Foulds (RF) distance RF(-), the completion-based RF(+), the generalized RF (GRF), the dissimilarity measure, the vectorial tree distance, and the geodesic distance in the extended Billera-Holmes-Vogtmann tree space. Among these, three RF-based methods, RF(-), RF(+), and GRF, were examined in detail on generated clusters of phylogenetic trees defined on different but mutually overlapping sets of taxa. Additionally, we reviewed nine related techniques, including leaf imputation methods, the tree edit distance, and visual comparison. A comparison of the related distance measures, highlighting their principal advantages and shortcomings, is provided. This review offers valuable insights into their applicability and performance, guiding the appropriate use of these metrics based on tree type (rooted or unrooted) and information type (topological or branch lengths).
Collapse
Affiliation(s)
- Wanlin Li
- Department of Computer ScienceUniversity of SherbrookeSherbrookeQuebecCanada
| | - Aleksandr Koshkarov
- Department of Computer ScienceUniversity of SherbrookeSherbrookeQuebecCanada
| | - Nadia Tahiri
- Department of Computer ScienceUniversity of SherbrookeSherbrookeQuebecCanada
| |
Collapse
|
2
|
Frost LA, Bedoya AM, Lagomarsino LP. Artifactual Orthologs and the Need for Diligent Data Exploration in Complex Phylogenomic Datasets: A Museomic Case Study from the Andean Flora. Syst Biol 2024; 73:308-322. [PMID: 38170162 DOI: 10.1093/sysbio/syad076] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2022] [Revised: 11/20/2023] [Accepted: 01/02/2024] [Indexed: 01/05/2024] Open
Abstract
The Andes mountains of western South America are a globally important biodiversity hotspot, yet there is a paucity of resolved phylogenies for plant clades from this region. Filling an important gap in our understanding of the World's richest flora, we present the first phylogeny of Freziera (Pentaphylacaceae), an Andean-centered, cloud forest radiation. Our dataset was obtained via hybrid-enriched target sequence capture of Angiosperms353 universal loci for 50 of the ca. 75 spp., obtained almost entirely from herbarium specimens. We identify high phylogenomic complexity in Freziera, including the presence of data artifacts. Via by-eye observation of gene trees, detailed examination of warnings from recently improved assembly pipelines, and gene tree filtering, we identified that artifactual orthologs (i.e., the presence of only one copy of a multicopy gene due to differential assembly) were an important source of gene tree heterogeneity that had a negative impact on phylogenetic inference and support. These artifactual orthologs may be common in plant phylogenomic datasets, where multiple instances of genome duplication are common. After accounting for artifactual orthologs as source of gene tree error, we identified a significant, but nonspecific signal of introgression using Patterson's D and f4 statistics. Despite phylogenomic complexity, we were able to resolve Freziera into 9 well-supported subclades whose evolution has been shaped by multiple evolutionary processes, including incomplete lineage sorting, historical gene flow, and gene duplication. Our results highlight the complexities of plant phylogenomics, which are heightened in Andean radiations, and show the impact of filtering data processing artifacts and standard filtering approaches on phylogenetic inference.
Collapse
Affiliation(s)
- Laura A Frost
- Shirley C. Tucker Herbarium, Department of Biological Sciences, Louisiana State University, Life Science Annex Building A257, Baton Rouge, LA 70803, USA
- Biology Department, University of South Alabama, 5871 USA N Dr, Mobile, AL 36688, USA
| | - Ana M Bedoya
- Shirley C. Tucker Herbarium, Department of Biological Sciences, Louisiana State University, Life Science Annex Building A257, Baton Rouge, LA 70803, USA
| | - Laura P Lagomarsino
- Shirley C. Tucker Herbarium, Department of Biological Sciences, Louisiana State University, Life Science Annex Building A257, Baton Rouge, LA 70803, USA
| |
Collapse
|
3
|
Balaban M, Jiang Y, Zhu Q, McDonald D, Knight R, Mirarab S. Generation of accurate, expandable phylogenomic trees with uDance. Nat Biotechnol 2024; 42:768-777. [PMID: 37500914 PMCID: PMC10818028 DOI: 10.1038/s41587-023-01868-8] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2022] [Accepted: 06/20/2023] [Indexed: 07/29/2023]
Abstract
Phylogenetic trees provide a framework for organizing evolutionary histories across the tree of life and aid downstream comparative analyses such as metagenomic identification. Methods that rely on single-marker genes such as 16S rRNA have produced trees of limited accuracy with hundreds of thousands of organisms, whereas methods that use genome-wide data are not scalable to large numbers of genomes. We introduce updating trees using divide-and-conquer (uDance), a method that enables updatable genome-wide inference using a divide-and-conquer strategy that refines different parts of the tree independently and can build off of existing trees, with high accuracy and scalability. With uDance, we infer a species tree of roughly 200,000 genomes using 387 marker genes, totaling 42.5 billion amino acid residues.
Collapse
Affiliation(s)
- Metin Balaban
- Bioinformatics and Systems Biology Graduate Program, University of California San Diego, La Jolla, CA, USA
| | - Yueyu Jiang
- Department of Electrical and Computer Engineering, University of California San Diego, La Jolla, CA, USA
| | - Qiyun Zhu
- Biodesign Center for Fundamental and Applied Microbiomics, Arizona State University, Tempe, AZ, USA
- School of Life Sciences, Arizona State University, Tempe, AZ, USA
| | - Daniel McDonald
- Department of Pediatrics, University of California San Diego, La Jolla, CA, USA
| | - Rob Knight
- Department of Pediatrics, University of California San Diego, La Jolla, CA, USA
- Department of Computer Science and Engineering, Jacobs School of Engineering, University of California San Diego, La Jolla, CA, USA
- Department of Bioengineering, University of California San Diego, La Jolla, CA, USA
- Center for Microbiome Innovation, Jacobs School of Engineering, University of California San Diego, La Jolla, CA, USA
| | - Siavash Mirarab
- Department of Electrical and Computer Engineering, University of California San Diego, La Jolla, CA, USA.
- Department of Computer Science and Engineering, Jacobs School of Engineering, University of California San Diego, La Jolla, CA, USA.
- Center for Microbiome Innovation, Jacobs School of Engineering, University of California San Diego, La Jolla, CA, USA.
| |
Collapse
|
4
|
Han Y, Molloy EK. Quartets enable statistically consistent estimation of cell lineage trees under an unbiased error and missingness model. Algorithms Mol Biol 2023; 18:19. [PMID: 38041123 PMCID: PMC10691101 DOI: 10.1186/s13015-023-00248-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2023] [Accepted: 11/19/2023] [Indexed: 12/03/2023] Open
Abstract
Cancer progression and treatment can be informed by reconstructing its evolutionary history from tumor cells. Although many methods exist to estimate evolutionary trees (called phylogenies) from molecular sequences, traditional approaches assume the input data are error-free and the output tree is fully resolved. These assumptions are challenged in tumor phylogenetics because single-cell sequencing produces sparse, error-ridden data and because tumors evolve clonally. Here, we study the theoretical utility of methods based on quartets (four-leaf, unrooted phylogenetic trees) in light of these barriers. We consider a popular tumor phylogenetics model, in which mutations arise on a (highly unresolved) tree and then (unbiased) errors and missing values are introduced. Quartets are then implied by mutations present in two cells and absent from two cells. Our main result is that the most probable quartet identifies the unrooted model tree on four cells. This motivates seeking a tree such that the number of quartets shared between it and the input mutations is maximized. We prove an optimal solution to this problem is a consistent estimator of the unrooted cell lineage tree; this guarantee includes the case where the model tree is highly unresolved, with error defined as the number of false negative branches. Lastly, we outline how quartet-based methods might be employed when there are copy number aberrations and other challenges specific to tumor phylogenetics.
Collapse
Affiliation(s)
- Yunheng Han
- Department of Computer Science, University of Maryland, College Park, MD, USA
| | - Erin K Molloy
- Department of Computer Science, University of Maryland, College Park, MD, USA.
- University of Maryland Institute for Advanced Computer Studies, College Park, MD, USA.
| |
Collapse
|
5
|
Wang YW, McKeon MC, Elmore H, Hess J, Golan J, Gage H, Mao W, Harrow L, Gonçalves SC, Hull CM, Pringle A. Invasive Californian death caps develop mushrooms unisexually and bisexually. Nat Commun 2023; 14:6560. [PMID: 37875491 PMCID: PMC10598064 DOI: 10.1038/s41467-023-42317-z] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2023] [Accepted: 10/05/2023] [Indexed: 10/26/2023] Open
Abstract
Canonical sexual reproduction among basidiomycete fungi involves the fusion of two haploid individuals of different mating types, resulting in a heterokaryotic mycelial body made up of genetically different nuclei. Using population genomics data and experiments, we discover mushrooms of the invasive and deadly Amanita phalloides can also be homokaryotic; evidence of sexual reproduction by single, unmated individuals. In California, genotypes of homokaryotic mushrooms are also found in heterokaryotic mushrooms, implying nuclei of homokaryotic mycelia are also involved in outcrossing. We find death cap mating is controlled by a single mating type locus, but the development of homokaryotic mushrooms appears to bypass mating type gene control. Ultimately, sporulation is enabled by nuclei able to reproduce alone as well as with others, and nuclei competent for both unisexuality and bisexuality have persisted in invaded habitats for at least 17 but potentially as long as 30 years. The diverse reproductive strategies of invasive death caps are likely facilitating its rapid spread, suggesting a profound similarity between plant, animal and fungal invasions.
Collapse
Affiliation(s)
- Yen-Wen Wang
- Department of Botany, University of Wisconsin-Madison, Madison, WI, USA.
| | - Megan C McKeon
- Department of Genetics, University of Wisconsin-Madison, Madison, WI, USA
- Department of Biomolecular Chemistry, University of Wisconsin-Madison, Madison, WI, USA
| | | | | | - Jacob Golan
- Department of Botany, University of Wisconsin-Madison, Madison, WI, USA
| | - Hunter Gage
- Department of Biomolecular Chemistry, University of Wisconsin-Madison, Madison, WI, USA
| | - William Mao
- Department of Botany, University of Wisconsin-Madison, Madison, WI, USA
| | - Lynn Harrow
- Department of Botany, University of Wisconsin-Madison, Madison, WI, USA
- Department of Chemistry, University of Wisconsin-Madison, Madison, WI, USA
| | - Susana C Gonçalves
- Centre for Functional Ecology, Department of Life Sciences, University of Coimbra, Coimbra, Portugal
| | - Christina M Hull
- Department of Biomolecular Chemistry, University of Wisconsin-Madison, Madison, WI, USA
- Department of Medical Microbiology & Immunology, University of Wisconsin-Madison, Madison, WI, USA
| | - Anne Pringle
- Department of Botany, University of Wisconsin-Madison, Madison, WI, USA
- Department of Bacteriology, University of Wisconsin-Madison, Madison, WI, USA
| |
Collapse
|
6
|
Cerca J, Cotoras DD, Santander CG, Bieker VC, Hutchins L, Morin-Lagos J, Prada CF, Kennedy S, Krehenwinkel H, Rominger AJ, Meier J, Dimitrov D, Struck TH, Gillespie RG. Multiple paths toward repeated phenotypic evolution in the spiny-leg adaptive radiation (Tetragnatha; Hawai'i). Mol Ecol 2023; 32:4971-4985. [PMID: 37515430 DOI: 10.1111/mec.17082] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2022] [Revised: 07/05/2023] [Accepted: 07/07/2023] [Indexed: 07/30/2023]
Abstract
The repeated evolution of phenotypes provides clear evidence for the role of natural selection in driving evolutionary change. However, the evolutionary origin of repeated phenotypes can be difficult to disentangle as it can arise from a combination of factors such as gene flow, shared ancestral polymorphisms or mutation. Here, we investigate the presence of these evolutionary processes in the Hawaiian spiny-leg Tetragnatha adaptive radiation, which includes four microhabitat-specialists or ecomorphs, with different body pigmentation and size (Green, Large Brown, Maroon, and Small Brown). We investigated the evolutionary history of this radiation using 76 newly generated low-coverage, whole-genome resequenced samples, along with phylogenetic and population genomic tools. Considering the Green ecomorph as the ancestral state, our results suggest that the Green ecomorph likely re-evolved once, the Large Brown and Maroon ecomorphs evolved twice and the Small Brown evolved three times. We found that the evolution of the Maroon and Small Brown ecomorphs likely involved ancestral hybridization events, while the Green and Large Brown ecomorphs likely evolved through novel mutations, despite a high rate of incomplete lineage sorting in the dataset. Our findings demonstrate that the repeated evolution of ecomorphs in the Hawaiian spiny-leg Tetragnatha is influenced by multiple evolutionary processes.
Collapse
Affiliation(s)
- José Cerca
- Berkeley Evolab, Department of Environmental Science, Policy, and Management, UC Berkeley, Berkeley, California, USA
- Frontiers in Evolutionary Zoology, Natural History Museum, University of Oslo, Oslo, Norway
- Department of Natural History, NTNU University Museum, Norwegian University of Science and Technology, Trondheim, Norway
- Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, Oslo, Norway
| | - Darko D Cotoras
- Department of Terrestrial Zoology, Senckenberg Research Institute and Natural History Museum, Frankfurt am Main, Germany
- Department of Entomology, California Academy of Sciences, San Francisco, California, USA
| | - Cindy G Santander
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Vanessa C Bieker
- Department of Natural History, NTNU University Museum, Norwegian University of Science and Technology, Trondheim, Norway
| | - Leke Hutchins
- Berkeley Evolab, Department of Environmental Science, Policy, and Management, UC Berkeley, Berkeley, California, USA
| | - Jaime Morin-Lagos
- Department of Natural History, NTNU University Museum, Norwegian University of Science and Technology, Trondheim, Norway
| | - Carlos F Prada
- Grupo de Investigación de Biología y Ecología de Artrópodos, Facultad de Ciencias, Universidad del Tolima, Tolima, Colombia
| | - Susan Kennedy
- Department of Biogeography, Trier University, Trier, Germany
| | | | - Andrew J Rominger
- School of Biology and Ecology, University of Maine, Orono, Maine, USA
| | - Joana Meier
- Department of Zoology, University of Cambridge, Cambridge, UK
- Tree of Life Programme, Sanger Institute, Hinxton, UK
| | - Dimitar Dimitrov
- Department of Natural History, University Museum of Bergen, University of Bergen, Bergen, Norway
| | - Torsten H Struck
- Frontiers in Evolutionary Zoology, Natural History Museum, University of Oslo, Oslo, Norway
| | - Rosemary G Gillespie
- Berkeley Evolab, Department of Environmental Science, Policy, and Management, UC Berkeley, Berkeley, California, USA
| |
Collapse
|
7
|
Nuñez LP, Gray LN, Weisrock DW, Burbrink FT. The Phylogenomic and Biogeographic History of the Gartersnakes, Watersnakes, and Allies (Natricidae: Thamnophiini). Mol Phylogenet Evol 2023:107844. [PMID: 37301486 DOI: 10.1016/j.ympev.2023.107844] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2022] [Revised: 06/01/2023] [Accepted: 06/03/2023] [Indexed: 06/12/2023]
Abstract
North American Thamnophiini (gartersnakes, watersnakes, brownsnakes, and swampsnakes) are an ecologically and phenotypically diverse temperate clade of snakes representing 61 species across 10 genera. In this study, we estimate phylogenetic trees using ∼3,700 ultraconserved elements (UCEs) for 76 specimens representing 75% of all Thamnophiini species. We infer phylogenies using multispecies coalescent methods and time calibrate them using the fossil record. We also conducted ancestral area estimation to identify how major biogeographic boundaries in North America affect broadscale diversification in the group. While most nodes exhibited strong statistical support, analysis of concordant data across gene trees reveals substantial heterogeneity. Ancestral area estimation demonstrated that the genus Thamnophis was the only taxon in this subfamily to cross the Western Continental Divide, even as other taxa dispersed southward toward the tropics. Additionally, levels of gene tree discordance are overall higher in transition zones between bioregions, including the Rocky Mountains. Therefore, the Western Continental Divide may be a significant transition zone structuring the diversification of Thamnophiini during the Neogene and Pleistocene. Here we show that despite high levels of discordance across gene trees, we were able to infer a highly resolved and well-supported phylogeny for Thamnophiini, which allows us to understand broadscale patterns of diversity and biogeography.
Collapse
Affiliation(s)
- Leroy P Nuñez
- Department of Herpetology, American Museum of Natural History, New York, NY, USA; Richard Gilder Graduate School, American Museum of Natural History, New York, NY, USA.
| | - Levi N Gray
- Fort Collins Science Center, United States Geological Survey, Guam, USA
| | - David W Weisrock
- Department of Biology, University of Kentucky, Lexington, KY, USA
| | - Frank T Burbrink
- Department of Herpetology, American Museum of Natural History, New York, NY, USA
| |
Collapse
|
8
|
Wang YW, McKeon MC, Elmore H, Hess J, Golan J, Gage H, Mao W, Harrow L, Gon ßalves SC, Hull CM, Pringle A. Invasive Californian death caps develop mushrooms unisexually and bisexually. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.01.30.525609. [PMID: 36778337 PMCID: PMC9915504 DOI: 10.1101/2023.01.30.525609] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Canonical sexual reproduction among basidiomycete fungi involves the fusion of two haploid individuals of different sexes, resulting in a heterokaryotic mycelial body made up of genetically different nuclei 1 . Using population genomics data, we discovered mushrooms of the deadly invasive Amanita phalloides are also homokaryotic, evidence of sexual reproduction by single individuals. In California, genotypes of homokaryotic mushrooms are also found in heterokaryotic mushrooms, implying nuclei of homokaryotic mycelia also promote outcrossing. We discovered death cap mating is controlled by a single mating-type locus ( A. phalloides is bipolar), but the development of homokaryotic mushrooms appears to bypass mating-type gene control. Ultimately, sporulation is enabled by nuclei able to reproduce alone as well as with others, and nuclei competent for both unisexuality and bisexuality have persisted in invaded habitats for at least 17 but potentially as long as 30 years. The diverse reproductive strategies of invasive death caps are likely facilitating its rapid spread, revealing a profound similarity between plant, animal and fungal invasions 2,3 .
Collapse
|
9
|
Zaharias P, Warnow T. Recent progress on methods for estimating and updating large phylogenies. Philos Trans R Soc Lond B Biol Sci 2022; 377:20210244. [PMID: 35989607 PMCID: PMC9393559 DOI: 10.1098/rstb.2021.0244] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2021] [Accepted: 01/07/2022] [Indexed: 12/20/2022] Open
Abstract
With the increased availability of sequence data and even of fully sequenced and assembled genomes, phylogeny estimation of very large trees (even of hundreds of thousands of sequences) is now a goal for some biologists. Yet, the construction of these phylogenies is a complex pipeline presenting analytical and computational challenges, especially when the number of sequences is very large. In the past few years, new methods have been developed that aim to enable highly accurate phylogeny estimations on these large datasets, including divide-and-conquer techniques for multiple sequence alignment and/or tree estimation, methods that can estimate species trees from multi-locus datasets while addressing heterogeneity due to biological processes (e.g. incomplete lineage sorting and gene duplication and loss), and methods to add sequences into large gene trees or species trees. Here we present some of these recent advances and discuss opportunities for future improvements. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.
Collapse
Affiliation(s)
- Paul Zaharias
- Department of Computer Science, University of Illinois Urbana-Champaign, Urbana, IL 61801, USA
| | - Tandy Warnow
- Department of Computer Science, University of Illinois Urbana-Champaign, Urbana, IL 61801, USA
| |
Collapse
|
10
|
Gatesy J, Springer MS. Phylogenomic Coalescent Analyses of Avian Retroelements Infer Zero-Length Branches at the Base of Neoaves, Emergent Support for Controversial Clades, and Ancient Introgressive Hybridization in Afroaves. Genes (Basel) 2022; 13:1167. [PMID: 35885951 PMCID: PMC9324441 DOI: 10.3390/genes13071167] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2022] [Revised: 06/20/2022] [Accepted: 06/21/2022] [Indexed: 01/25/2023] Open
Abstract
Retroelement insertions (RIs) are low-homoplasy characters that are ideal data for addressing deep evolutionary radiations, where gene tree reconstruction errors can severely hinder phylogenetic inference with DNA and protein sequence data. Phylogenomic studies of Neoaves, a large clade of birds (>9000 species) that first diversified near the Cretaceous−Paleogene boundary, have yielded an array of robustly supported, contradictory relationships among deep lineages. Here, we reanalyzed a large RI matrix for birds using recently proposed quartet-based coalescent methods that enable inference of large species trees including branch lengths in coalescent units, clade-support, statistical tests for gene flow, and combined analysis with DNA-sequence-based gene trees. Genome-scale coalescent analyses revealed extremely short branches at the base of Neoaves, meager branch support, and limited congruence with previous work at the most challenging nodes. Despite widespread topological conflicts with DNA-sequence-based trees, combined analyses of RIs with thousands of gene trees show emergent support for multiple higher-level clades (Columbea, Passerea, Columbimorphae, Otidimorphae, Phaethoquornithes). RIs express asymmetrical support for deep relationships within the subclade Afroaves that hints at ancient gene flow involving the owl lineage (Strigiformes). Because DNA-sequence data are challenged by gene tree-reconstruction error, analysis of RIs represents one approach for improving gene tree-based methods when divergences are deep, internodes are short, terminal branches are long, and introgressive hybridization further confounds species−tree inference.
Collapse
Affiliation(s)
- John Gatesy
- Division of Vertebrate Zoology, American Museum of Natural History, New York, NY 10024, USA
| | - Mark S. Springer
- Department of Evolution, Ecology, and Organismal Biology, University of California, Riverside, CA 92521, USA;
| |
Collapse
|
11
|
Mai U, Mirarab S. Completing gene trees without species trees in sub-quadratic time. Bioinformatics 2022; 38:1532-1541. [PMID: 34978565 DOI: 10.1093/bioinformatics/btab875] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2021] [Revised: 11/27/2021] [Accepted: 12/30/2021] [Indexed: 02/03/2023] Open
Abstract
MOTIVATION As genome-wide reconstruction of phylogenetic trees becomes more widespread, limitations of available data are being appreciated more than ever before. One issue is that phylogenomic datasets are riddled with missing data, and gene trees, in particular, almost always lack representatives from some species otherwise available in the dataset. Since many downstream applications of gene trees require or can benefit from access to complete gene trees, it will be beneficial to algorithmically complete gene trees. Also, gene trees are often unrooted, and rooting them is useful for downstream applications. While completing and rooting a gene tree with respect to a given species tree has been studied, those problems are not studied in depth when we lack such a reference species tree. RESULTS We study completion of gene trees without a need for a reference species tree. We formulate an optimization problem to complete the gene trees while minimizing their quartet distance to the given set of gene trees. We extend a seminal algorithm by Brodal et al. to solve this problem in quasi-linear time. In simulated studies and on a large empirical data, we show that completion of gene trees using other gene trees is relatively accurate and, unlike the case where a species tree is available, is unbiased. AVAILABILITY AND IMPLEMENTATION Our method, tripVote, is available at https://github.com/uym2/tripVote. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Uyen Mai
- Department of Computer Science and Engineering, University of California San Diego, San Diego, CA 92093, USA
| | - Siavash Mirarab
- Department of Electrical and Computer Engineering, University of California San Diego, San Diego, CA 92093, USA
| |
Collapse
|
12
|
Liu B, Warnow T. Scalable Species Tree Inference with External Constraints. J Comput Biol 2022; 29:664-678. [PMID: 35196115 DOI: 10.1089/cmb.2021.0543] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Species tree inference is a basic step in biological discovery, but discordance between gene trees creates analytical challenges and large data sets create computational challenges. Although there is generally some information available about the species trees that could be used to speed up the estimation, only one species tree estimation method that addresses gene tree discordance-ASTRAL-J, a recent development in the ASTRAL family of methods-is able to use this information. Here we describe two new methods, NJst-J and FASTRAL-J, that can estimate the species tree, given a partial knowledge of the species tree in the form of a nonbinary unrooted constraint tree. We show that both NJst-J and FASTRAL-J are much faster than ASTRAL-J and we prove that all three methods are statistically consistent under the multispecies coalescent model subject to this constraint. Our extensive simulation study shows that both FASTRAL-J and NJst-J provide advantages over ASTRAL-J: both are faster (and NJst-J is particularly fast), and FASTRAL-J is generally at least as accurate as ASTRAL-J. An analysis of the Avian Phylogenomics Project data set with 48 species and 14,446 genes presents additional evidence of the value of FASTRAL-J over ASTRAL-J (and both over ASTRAL), with dramatic reductions in running time (20 hours for default ASTRAL, and minutes or seconds for ASTRAL-J and FASTRAL-J, respectively).
Collapse
Affiliation(s)
- Baqiao Liu
- Department of Computer Science, University of Illinois Urbana-Champaign, Urbana, Illinois, USA
| | - Tandy Warnow
- Department of Computer Science, University of Illinois Urbana-Champaign, Urbana, Illinois, USA
| |
Collapse
|
13
|
Tilic E, Stiller J, Campos E, Pleijel F, Rouse GW. Phylogenomics resolves ambiguous relationships within Aciculata (Errantia, Annelida). Mol Phylogenet Evol 2021; 166:107339. [PMID: 34751138 DOI: 10.1016/j.ympev.2021.107339] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2021] [Revised: 10/23/2021] [Accepted: 10/27/2021] [Indexed: 10/20/2022]
Abstract
Aciculata (Eunicida + Phyllodocida) is among the largest clades of annelids, comprising almost half of the known diversity of all marine annelids. Despite the group's large size and biological importance, most phylogenomic studies on Annelida to date have had a limited sampling of this clade. The phylogenetic placement of many clades within Phyllodocida in particular has remained poorly understood. To resolve the relationships within Aciculata we conducted a large-scale phylogenomic analysis based on 24 transcriptomes (13 new), chosen to represent many family-ranked taxa that have never been included in a broad phylogenomic study. Our sampling also includes several enigmatic taxa with challenging phylogenetic placement, such as Histriobdella, Struwela, Lacydonia, Pilargis and the holopelagic worms Lopadorrhynchus, Travisiopsis and Tomopteris. Our robust phylogeny allows us to name and place some of these problematic clades and has significant implications on the systematics of the group. Within Eunicida we reinstate the names Eunicoidea and Oenonoidea. Within Phyllodocida we delineate Phyllodociformia, Glyceriformia, Nereidiformia, Nephtyiformia and Aphroditiformia. Phyllodociformia now includes: Lacydonia, Typhloscolecidae, Lopadorrhynchidae and Phyllodocidae. Nephtyiformia includes Nephtyidae and Pilargidae. We also broaden the delineation of Glyceriformia to include Sphaerodoridae, Tomopteridae and Glyceroidea (Glyceridae + Goniadidae). Furthermore, our study demonstrates and explores how conflicting, yet highly supported topologies can result from confounding signals in gene trees.
Collapse
Affiliation(s)
- Ekin Tilic
- Scripps Institution of Oceanography, UC San Diego, La Jolla, CA, USA; Institute of Evolutionary Biology and Animal Ecology, University of Bonn, Germany; Marine Biological Section, Department of Biology, University of Copenhagen, DK-2100 Copenhagen, Denmark.
| | - Josefin Stiller
- Centre for Biodiversity Genomics, Section for Ecology and Evolution, Department of Biology, University of Copenhagen, DK-2100 Copenhagen, Denmark
| | - Ernesto Campos
- Facultad de Ciencias, Universidad Autónoma de Baja California. Ensenada, Baja California, México
| | - Fredrik Pleijel
- Department of Marine Sciences, University of Gothenburg, Tjärnö, Sweden
| | - Greg W Rouse
- Scripps Institution of Oceanography, UC San Diego, La Jolla, CA, USA.
| |
Collapse
|
14
|
Tilic E, Sayyari E, Stiller J, Mirarab S, Rouse GW. More is needed—Thousands of loci are required to elucidate the relationships of the ‘flowers of the sea’ (Sabellida, Annelida). Mol Phylogenet Evol 2020; 151:106892. [DOI: 10.1016/j.ympev.2020.106892] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2020] [Revised: 06/09/2020] [Accepted: 06/10/2020] [Indexed: 11/17/2022]
|