Phylogenomic Analysis of Concatenated Ultraconserved Elements Reveals the Recent Evolutionary Radiation of the Fairy Wrasses (Teleostei: Labridae: Cirrhilabrus)

2021 ◽  
Author(s):  
Yi-Kai Tea ◽  
Xin Xu ◽  
Joseph D DiBattista ◽  
Nathan Lo ◽  
Peter F Cowman ◽  
...  
Keyword(s):  
Phylogenetic Signal ◽  
Estimation Error ◽  
Gene Tree ◽  
Molecular Dating ◽  
Morphological Data ◽  
Phylogenomic Analysis ◽  
Gene Trees ◽  
Data Set ◽  
Ultraconserved Elements ◽  
Species Complexes

Abstract The fairy wrasses (genus Cirrhilabrus) are among the most successful of the extant wrasse lineages (Teleostei: Labridae), with their 61 species accounting for nearly 10$\%$ of the family. Although species complexes within the genus have been diagnosed on the basis of coloration patterns and synapomorphies, attempts to resolve evolutionary relationships among these groups using molecular and morphological data have largely been unsuccessful. Here, we use a phylogenomic approach with a data set comprising 991 ultraconserved elements (UCEs) and mitochondrial COI to uncover the evolutionary history and patterns of temporal and spatial diversification of the fairy wrasses. Our analyses of phylogenetic signal suggest that most gene-tree incongruence is caused by estimation error, leading to poor resolution in a summary-coalescent analysis of the data. In contrast, analyses of concatenated sequences are able to resolve the major relationships of Cirrhilabrus. We determine the placements of species that were previously regarded as incertae sedis and find evidence for the nesting of Conniella, an unusual, monotypic genus, within Cirrhilabrus. Our relaxed-clock dating analysis indicates that the major divergences within the genus occurred around the Miocene–Pliocene boundary, followed by extensive cladogenesis of species complexes in the Pliocene–Pleistocene. Biogeographic reconstruction suggests that the fairy wrasses emerged within the Coral Triangle, with episodic fluctuations of sea levels during glacial cycles coinciding with shallow divergence events but providing few opportunities for more widespread dispersal. Our study demonstrates both the resolving power and limitations of UCEs across shallow timescales where there is substantial estimation error in individual gene trees.[Biogeography; concatenation; gene genealogy interrogation; gene trees; molecular dating; summary coalescent; UCEs.]

scholarly journals Optimizing Phylogenomics with Rapidly Evolving Long Exons: Comparison with Anchored Hybrid Enrichment and Ultraconserved Elements

2019 ◽  
Vol 37 (3) ◽  
pp. 904-922 ◽  
Author(s):  
Benjamin R Karin ◽  
Tony Gamble ◽  
Todd R Jackman
Keyword(s):  
Evolutionary Rate ◽  
Estimation Error ◽  
Model Misspecification ◽  
Gene Tree ◽  
Species Tree ◽  
Gene Trees ◽  
Data Set ◽  
Ultraconserved Elements ◽  
Trade Offs ◽  
Anchored Hybrid Enrichment

Abstract Marker selection has emerged as an important component of phylogenomic study design due to rising concerns of the effects of gene tree estimation error, model misspecification, and data-type differences. Researchers must balance various trade-offs associated with locus length and evolutionary rate among other factors. The most commonly used reduced representation data sets for phylogenomics are ultraconserved elements (UCEs) and Anchored Hybrid Enrichment (AHE). Here, we introduce Rapidly Evolving Long Exon Capture (RELEC), a new set of loci that targets single exons that are both rapidly evolving (evolutionary rate faster than RAG1) and relatively long in length (>1,500 bp), while at the same time avoiding paralogy issues across amniotes. We compare the RELEC data set to UCEs and AHE in squamate reptiles by aligning and analyzing orthologous sequences from 17 squamate genomes, composed of 10 snakes and 7 lizards. The RELEC data set (179 loci) outperforms AHE and UCEs by maximizing per-locus genetic variation while maintaining presence and orthology across a range of evolutionary scales. RELEC markers show higher phylogenetic informativeness than UCE and AHE loci, and RELEC gene trees show greater similarity to the species tree than AHE or UCE gene trees. Furthermore, with fewer loci, RELEC remains computationally tractable for full Bayesian coalescent species tree analyses. We contrast RELEC to and discuss important aspects of comparable methods, and demonstrate how RELEC may be the most effective set of loci for resolving difficult nodes and rapid radiations. We provide several resources for capturing or extracting RELEC loci from other amniote groups.

scholarly journals Non-parametric correction of estimated gene trees using TRACTION

2020 ◽  
Vol 15 (1) ◽  
Author(s):  
Sarah Christensen ◽  
Erin K. Molloy ◽  
Pranjal Vachaspati ◽  
Ananya Yammanuru ◽  
Tandy Warnow
Keyword(s):  
Incomplete Lineage Sorting ◽  
Phylogenetic Signal ◽  
Estimation Error ◽  
Single Gene ◽  
Gene Tree ◽  
Species Tree ◽  
Computational Techniques ◽  
Gene Trees ◽  
Species Trees ◽  
Sequencing Data

Abstract Motivation Estimated gene trees are often inaccurate, due to insufficient phylogenetic signal in the single gene alignment, among other causes. Gene tree correction aims to improve the accuracy of an estimated gene tree by using computational techniques along with auxiliary information, such as a reference species tree or sequencing data. However, gene trees and species trees can differ as a result of gene duplication and loss (GDL), incomplete lineage sorting (ILS), and other biological processes. Thus gene tree correction methods need to take estimation error as well as gene tree heterogeneity into account. Many prior gene tree correction methods have been developed for the case where GDL is present. Results Here, we study the problem of gene tree correction where gene tree heterogeneity is instead due to ILS and/or HGT. We introduce TRACTION, a simple polynomial time method that provably finds an optimal solution to the RF-optimal tree refinement and completion (RF-OTRC) Problem, which seeks a refinement and completion of a singly-labeled gene tree with respect to a given singly-labeled species tree so as to minimize the Robinson−Foulds (RF) distance. Our extensive simulation study on 68,000 estimated gene trees shows that TRACTION matches or improves on the accuracy of well-established methods from the GDL literature when HGT and ILS are both present, and ties for best under the ILS-only conditions. Furthermore, TRACTION ties for fastest on these datasets. We also show that a naive generalization of the RF-OTRC problem to multi-labeled trees is possible, but can produce misleading results where gene tree heterogeneity is due to GDL.

scholarly journals Optimizing Phylogenomics with Rapidly Evolving Long Exons: Comparison with Anchored Hybrid Enrichment and Ultraconserved Elements

2019 ◽  
Author(s):  
Benjamin R. Karin ◽  
Tony Gamble ◽  
Todd R. Jackman
Keyword(s):  
Evolutionary Rate ◽  
Estimation Error ◽  
Model Misspecification ◽  
Gene Tree ◽  
Species Tree ◽  
Gene Trees ◽  
Ultraconserved Elements ◽  
Marker Selection ◽  
Trade Offs ◽  
Anchored Hybrid Enrichment

AbstractMarker selection has emerged as an important component of phylogenomic study design due to rising concerns of the effects of gene tree estimation error, model misspecification, and data-type differences. Researchers must balance various trade-offs associated with locus length and evolutionary rate among other factors. The most commonly used reduced representation datasets for phylogenomics are ultraconserved elements (UCEs) and Anchored Hybrid Enrichment (AHE). Here, we introduce Rapidly Evolving Long Exon Capture (RELEC), a new set of loci that targets single exons that are both rapidly evolving (evolutionary rate faster than RAG1) and relatively long in length (greater than 1,500 bp), while at the same time avoiding paralogy issues across amniotes. We compare the RELEC dataset to UCEs and AHE in squamate reptiles by aligning and analyzing orthologous sequences from 17 squamate genomes, composed of ten snakes and seven lizards. The RELEC dataset (179 loci) outperforms AHE and UCEs by maximizing per-locus genetic variation while maintaining presence and orthology across a range of evolutionary scales. RELEC markers show higher phylogenetic informativeness than UCE and AHE loci, and RELEC gene trees show greater similarity to the species tree than AHE or UCE gene trees. Furthermore, with fewer loci, RELEC remains computationally tractable for full Bayesian coalescent species tree analyses. We contrast RELEC to and discuss important aspects of comparable methods, and demonstrate how RELEC may be the most effective set of loci for resolving difficult nodes and rapid radiations. We provide several resources for capturing or extracting RELEC loci from other amniote groups.

scholarly journals A phylogenomic resolution of the sea urchin tree of life

2018 ◽  
Author(s):  
Nicolás Mongiardino Koch ◽  
Simon E. Coppard ◽  
Harilaos A. Lessios ◽  
Derek E. G. Briggs ◽  
Rich Mooi ◽  
...  
Keyword(s):  
De Novo ◽  
Sea Urchins ◽  
Phylogenetic Signal ◽  
Gene Tree ◽  
Tree Of Life ◽  
Morphological Data ◽  
Morphological Evidence ◽  
Phylogenomic Analysis ◽  
Marine Animals ◽  
Sand Dollars

AbstractBackgroundEchinoidea is a clade of marine animals including sea urchins, heart urchins, sand dollars and sea biscuits. Found in benthic habitats across all latitudes, echinoids are key components of marine communities such as coral reefs and kelp forests. A little over 1,000 species inhabit the oceans today, a diversity that traces its roots back at least to the Permian. Although much effort has been devoted to elucidating the echinoid tree of life using a variety of morphological data, molecular attempts have relied on only a handful of genes. Both of these approaches have had limited success at resolving the deepest nodes of the tree, and their disagreement over the positions of a number of clades remains unresolved.ResultsWe performed de novo sequencing and assembly of 17 transcriptomes to complement available genomic resources of sea urchins and produce the first phylogenomic analysis of the clade. Multiple methods of probabilistic inference recovered identical topologies, with virtually all nodes showing maximum support. In contrast, the coalescent-based method ASTRAL-II resolved one node differently, a result apparently driven by gene tree error induced by evolutionary rate heterogeneity. Regardless of the method employed, our phylogenetic structure deviates from the currently accepted classification of echinoids, with neither Acroechinoidea (all euechinoids except echinothurioids), nor Clypeasteroida (sand dollars and sea biscuits) being monophyletic as currently defined. We demonstrate the strength and distribution of phylogenetic signal throughout the genome for novel resolutions of these lineages and rule out systematic biases as possible explanations.ConclusionsOur investigation substantially augments the molecular resources available for sea urchins, providing the first transcriptomes for many of its main lineages. Using this expanded genomic dataset, we resolve the position of several clades in agreement with early molecular analyses but in disagreement with morphological data. Our efforts settle multiple phylogenetic uncertainties, including the position of the enigmatic deep-sea echinothurioids and the identity of the sister clade to sand dollars. We offer a detailed assessment of evolutionary scenarios that could reconcile our findings with morphological evidence, opening up new lines of research into the development and evolutionary history of this ancient clade.

scholarly journals Phylogenetic Conflicts, Combinability, and Deep Phylogenomics in Plants

2019 ◽  
Vol 69 (3) ◽  
pp. 579-592 ◽  
Author(s):  
Stephen A Smith ◽  
Nathanael Walker-Hale ◽  
Joseph F Walker ◽  
Joseph W Brown
Keyword(s):  
Phylogenetic Relationships ◽  
Phylogenetic Signal ◽  
Gene Tree ◽  
Species Tree ◽  
Gene Trees ◽  
Data Set ◽  
Tree Inference ◽  
Tree Methods ◽  
Plant Data ◽  
Inference Methods

Abstract Studies have demonstrated that pervasive gene tree conflict underlies several important phylogenetic relationships where different species tree methods produce conflicting results. Here, we present a means of dissecting the phylogenetic signal for alternative resolutions within a data set in order to resolve recalcitrant relationships and, importantly, identify what the data set is unable to resolve. These procedures extend upon methods for isolating conflict and concordance involving specific candidate relationships and can be used to identify systematic error and disambiguate sources of conflict among species tree inference methods. We demonstrate these on a large phylogenomic plant data set. Our results support the placement of Amborella as sister to the remaining extant angiosperms, Gnetales as sister to pines, and the monophyly of extant gymnosperms. Several other contentious relationships, including the resolution of relationships within the bryophytes and the eudicots, remain uncertain given the low number of supporting gene trees. To address whether concatenation of filtered genes amplified phylogenetic signal for relationships, we implemented a combinatorial heuristic to test combinability of genes. We found that nested conflicts limited the ability of data filtering methods to fully ameliorate conflicting signal amongst gene trees. These analyses confirmed that the underlying conflicting signal does not support broad concatenation of genes. Our approach provides a means of dissecting a specific data set to address deep phylogenetic relationships while also identifying the inferential boundaries of the data set. [Angiosperms; coalescent; gene-tree conflict; genomics; phylogenetics; phylogenomics.]

scholarly journals Phylogenetic signal is associated with the degree of variation in root-to-tip distances

2020 ◽  
Author(s):  
Mezzalina Vankan ◽  
Simon Y.W. Ho ◽  
Carolina Pardo-Diaz ◽  
David A. Duchêne
Keyword(s):  
Sequence Data ◽  
Phylogenetic Signal ◽  
Gene Tree ◽  
Phylogenetic Inference ◽  
Genomic Region ◽  
Published Data ◽  
Data Sets ◽  
Gene Trees ◽  
Species Trees ◽  
Data Set

AbstractThe phylogenetic information contained in sequence data is partly determined by the overall rate of nucleotide substitution in the genomic region in question. However, phylogenetic signal is affected by various other factors, such as heterogeneity in substitution rates across lineages. These factors might be able to predict the phylogenetic accuracy of any given gene in a data set. We examined the association between the accuracy of phylogenetic inference across genes and several characteristics of branch lengths in phylogenomic data. In a large number of published data sets, we found that the accuracy of phylogenetic inference from genes was consistently associated with their mean statistical branch support and variation in their gene tree root-to-tip distances, but not with tree length and stemminess. Therefore, a signal of constant evolutionary rates across lineages appears to be beneficial for phylogenetic inference. Identifying the causes of variation in root-to-tip lengths in gene trees also offers a potential way forward to increase congruence in the signal across genes and improve estimates of species trees from phylogenomic data sets.

scholarly journals Detecting destabilizing species in the phylogenetic backbone of Potentilla (Rosaceae) using low-copy nuclear markers

AoB Plants ◽  
2020 ◽  
Vol 12 (3) ◽  
Author(s):  
Nannie L Persson ◽  
Ingrid Toresen ◽  
Heidi Lie Andersen ◽  
Jenny E E Smedmark ◽  
Torsten Eriksson
Keyword(s):  
Incomplete Lineage Sorting ◽  
Gene Tree ◽  
Data Sets ◽  
Gene Trees ◽  
Lineage Sorting ◽  
Nuclear Markers ◽  
Data Set ◽  
Multispecies Coalescent ◽  
Single Marker ◽  
Named Groups

Abstract The genus Potentilla (Rosaceae) has been subjected to several phylogenetic studies, but resolving its evolutionary history has proven challenging. Previous analyses recovered six, informally named, groups: the Argentea, Ivesioid, Fragarioides, Reptans, Alba and Anserina clades, but the relationships among some of these clades differ between data sets. The Reptans clade, which includes the type species of Potentilla, has been noticed to shift position between plastid and nuclear ribosomal data sets. We studied this incongruence by analysing four low-copy nuclear markers, in addition to chloroplast and nuclear ribosomal data, with a set of Bayesian phylogenetic and Multispecies Coalescent (MSC) analyses. A selective taxon removal strategy demonstrated that the included representatives from the Fragarioides clade, P. dickinsii and P. fragarioides, were the main sources of the instability seen in the trees. The Fragarioides species showed different relationships in each gene tree, and were only supported as a monophyletic group in a single marker when the Reptans clade was excluded from the analysis. The incongruences could not be explained by allopolyploidy, but rather by homoploid hybridization, incomplete lineage sorting or taxon sampling effects. When P. dickinsii and P. fragarioides were removed from the data set, a fully resolved, supported backbone phylogeny of Potentilla was obtained in the MSC analysis. Additionally, indications of autopolyploid origins of the Reptans and Ivesioid clades were discovered in the low-copy gene trees.

scholarly journals Whole Genome Shotgun Phylogenomics Resolves the Pattern and Timing of Swallowtail Butterfly Evolution

2019 ◽  
Vol 69 (1) ◽  
pp. 38-60 ◽  
Author(s):  
Rémi Allio ◽  
Céline Scornavacca ◽  
Benoit Nabholz ◽  
Anne-Laure Clamens ◽  
Felix AH Sperling ◽  
...  
Keyword(s):  
Missing Data ◽  
Shotgun Sequencing ◽  
Whole Genome Shotgun ◽  
Amino Acid Sequences ◽  
Molecular Dating ◽  
Whole Genome ◽  
Gene Trees ◽  
Data Set ◽  
Whole Genomes ◽  
Bayesian Mixture Models

Abstract Evolutionary relationships have remained unresolved in many well-studied groups, even though advances in next-generation sequencing and analysis, using approaches such as transcriptomics, anchored hybrid enrichment, or ultraconserved elements, have brought systematics to the brink of whole genome phylogenomics. Recently, it has become possible to sequence the entire genomes of numerous nonbiological models in parallel at reasonable cost, particularly with shotgun sequencing. Here, we identify orthologous coding sequences from whole-genome shotgun sequences, which we then use to investigate the relevance and power of phylogenomic relationship inference and time-calibrated tree estimation. We study an iconic group of butterflies—swallowtails of the family Papilionidae—that has remained phylogenetically unresolved, with continued debate about the timing of their diversification. Low-coverage whole genomes were obtained using Illumina shotgun sequencing for all genera. Genome assembly coupled to BLAST-based orthology searches allowed extraction of 6621 orthologous protein-coding genes for 45 Papilionidae species and 16 outgroup species (with 32% missing data after cleaning phases). Supermatrix phylogenomic analyses were performed with both maximum-likelihood (IQ-TREE) and Bayesian mixture models (PhyloBayes) for amino acid sequences, which produced a fully resolved phylogeny providing new insights into controversial relationships. Species tree reconstruction from gene trees was performed with ASTRAL and SuperTriplets and recovered the same phylogeny. We estimated gene site concordant factors to complement traditional node-support measures, which strengthens the robustness of inferred phylogenies. Bayesian estimates of divergence times based on a reduced data set (760 orthologs and 12% missing data) indicate a mid-Cretaceous origin of Papilionoidea around 99.2 Ma (95% credibility interval: 68.6–142.7 Ma) and Papilionidae around 71.4 Ma (49.8–103.6 Ma), with subsequent diversification of modern lineages well after the Cretaceous-Paleogene event. These results show that shotgun sequencing of whole genomes, even when highly fragmented, represents a powerful approach to phylogenomics and molecular dating in a group that has previously been refractory to resolution.

scholarly journals DNA Barcoding: Mixed results for big-headed flies (Diptera: Pipunculidae)

Zootaxa ◽  
2007 ◽  
Vol 1423 (1) ◽  
pp. 1-26 ◽  
Author(s):  
JEFFREY H. SKEVINGTON ◽  
CHRISTIAN KEHLMAIER ◽  
GUNILLA STÅHLS
Keyword(s):  
Dna Barcoding ◽  
Sequence Data ◽  
Phylogenetic Signal ◽  
Morphological Characters ◽  
Molecular Data ◽  
Morphological Data ◽  
Base Pairs ◽  
Data Set ◽  
Taxonomic Decision ◽  
Morphological Species

Sequence data from 658 base pairs of mitochondrial cytochrome c oxidase I (cox1) were analysed for 28 described species of Pipunculidae (Diptera) in an effort to test the concept of DNA Barcoding on this family. Two recently revised but distantly related pipunculid lineages with presumed different evolutionary histories were used for the test (Clistoabdominalis Skevington, 2001 and Nephrocerus Zetterstedt, 1838). An effort was made to test the concept using sister taxa and morphologically similar sibling species swarms in these two genera. Morphological species concepts for Clistoabdominalis taxa were either supported by cox1 data or found to be too broad. Most of the discordance could be accounted for after reassessing morphological characters. In these cases, the molecular data were invaluable in assisting taxonomic decision-making. The radiation of Nearctic species of Nephrocerus could not be diagnosed using cox1. The ability of cox1 to recover phylogenetic signal was also tested on Clistoabdominalis. Morphological data for Clistoabdominalis were combined with the molecular data set. The pipunculid phylogeny from molecular data closely resembles the published phylogeny based on morphology. Partitioned Bremer support is used to localize areas of conflict between the datasets.

scholarly journals Phylogenetic conflicts, combinability, and deep phylogenomics in plants

2018 ◽  
Author(s):  
Stephen A. Smith ◽  
Nathanael Walker-Hale ◽  
Joseph F. Walker ◽  
Joseph W. Brown
Keyword(s):  
Phylogenetic Relationships ◽  
Phylogenetic Signal ◽  
Gene Tree ◽  
Species Tree ◽  
Gene Trees ◽  
Data Filtering ◽  
Tree Inference ◽  
Tree Methods ◽  
Inference Methods ◽  
Species Tree Inference

AbstractStudies have demonstrated that pervasive gene tree conflict underlies several important phylogenetic relationships where different species tree methods produce conflicting results. Here, we present a means of dissecting the phylogenetic signal for alternative resolutions within a dataset in order to resolve recalcitrant relationships and, importantly, identify what the dataset is unable to resolve. These procedures extend upon methods for isolating conflict and concordance involving specific candidate relationships and can be used to identify systematic error and disambiguate sources of conflict among species tree inference methods. We demonstrate these on a large phylogenomic plant dataset. Our results support the placement of Amborella as sister to the remaining extant angiosperms, Gnetales as sister to pines, and the monophyly of extant gymnosperms. Several other contentious relationships, including the resolution of relationships within the bryophytes and the eudicots, remain uncertain given the low number of supporting gene trees. To address whether concatenation of filtered genes amplified phylogenetic signal for relationships, we implemented a combinatorial heuristic to test combinability of genes. We found that nested conflicts limited the ability of data filtering methods to fully ameliorate conflicting signal amongst gene trees. These analyses confirmed that the underlying conflicting signal does not support broad concatenation of genes. Our approach provides a means of dissecting a specific dataset to address deep phylogenetic relationships while also identifying the inferential boundaries of the dataset.

Sign in / Sign up

Export Citation Format

Share Document