Find research datasets worth reusing

Data from: Phylogenomic analyses of Crassiclitellata support major Northern and Southern Hemisphere clades and a Pangaean origin for earthworms

Background: Earthworms (Crassiclitellata) are a diverse group of annelids of substantial ecological and economic importance. Earthworms are primarily terrestrial infaunal animals, and as such are probably rather poor natural dispersers. Therefore, the near global distribution of earthworms reflects an old and likely complex evolutionary history. Despite a long-standing interest in Crassiclitellata, relationships among and within major clades remain unresolved. Methods: In this study, we evaluate crassiclitellate phylogenetic relationships using 38 new transcriptomes in combination with publicly available transcriptome data. Our data include representatives of nearly all extant earthworm families and a representative of Moniligastridae, another terrestrial annelid group thought to be closely related to Crassiclitellata. We use a series of differentially filtered data matrices and analyses to examine the effects of data partitioning, missing data, compositional and branch-length heterogeneity, and outgroup inclusion. Results and discussion: We recover a consistent, strongly supported ingroup topology irrespective of differences in methodology. The topology supports two major earthworm clades, each of which consists of a Northern Hemisphere subclade and a Southern Hemisphere subclade. Divergence time analysis results are concordant with the hypothesis that these north-south splits are the result of the breakup of the supercontinent Pangaea. Conclusions: These results support several recently proposed revisions to the classical understanding of earthworm phylogeny, reveal two major clades that seem to reflect Pangaean distributions, and raise new questions about earthworm evolutionary relationships.

Data from: Anchored phylogenomics improves the resolution of evolutionary relationships in the rapid radiation of Protea L.

PREMISE OF THE STUDY: Estimating phylogenetic relationships in relatively recent evolutionary radiations is challenging, especially if short branches associated with recent divergence result in multiple gene tree histories. We combine anchored enrichment next-generation sequencing with species tree analyses to produce a robust estimate of phylogenetic relationships in the genus Protea (Proteaceae), an iconic radiation in South Africa. METHODS: We sampled multiple individuals within 59 out of 112 species of Protea and 6 outgroup species for a total of 163 individuals, and obtained sequences for 498 low-copy, orthologous nuclear loci using anchored phylogenomics. We compare several approaches for building species trees, and explore gene tree–species tree discrepancies to determine whether poor phylogenetic resolution reflects a lack of informative sites, incomplete lineage sorting, or hybridization. KEY RESULTS: Phylogenetic estimates from species tree approaches are similar to one another and recover previously well-supported clades within Protea, in addition to providing well-supported phylogenetic hypotheses for previously poorly resolved intrageneric relationships. Individual gene trees are markedly different from one another and from species trees. Nonetheless, analyses indicate that differences among gene trees occur primarily concerning clades supported by short branches. CONCLUSIONS: Species tree methods using hundreds of nuclear loci provided strong support for many previously unresolved relationships in the radiation of the genus Protea. In cases where support for particular relationships remains low, these appear to arise from few informative sites and lack of information rather than strongly supported disagreement among gene trees.

Data from: Phylogenomics of Aplacophora (Mollusca, Aculifera) and a solenogaster without a foot

Recent molecular phylogenetic investigations strongly supported placement of the shell-less, worm-shaped aplacophoran molluscs (Solenogastres and Caudofoveata) and chitons (Polyplacophora) in a clade called Aculifera, which is the sister taxon of all other molluscs. Thus, understanding the evolutionary history of aculiferan molluscs is important for understanding early molluscan evolution. In particular, fundamental questions about evolutionary relationships within Aplacophora have long been unanswered. Here, we supplemented the paucity of available data with transcriptomes from 25 aculiferans and conducted phylogenomic analyses on datasets with up to 525 genes and 75,914 amino acid positions. Our results indicate that aplacophoran taxonomy requires revision as several traditionally recognized groups are non-monophyletic. Most notably, Cavibelonia, the solenogaster taxon defined by hollow sclerites, is polyphyletic, suggesting parallel evolution of hollow sclerites in multiple lineages. Moreover, we describe Apodomenia enigmatica sp. nov., a bizarre new species that appears to be a morphological intermediate between Solenogastres and Caudofoveata. This animal is not a missing link however; molecular and morphological studies show that it is a derived solenogaster that lacks a foot, mantle cavity, and radula. Taken together, these results shed light on the evolutionary history of Aplacophora and reveal a surprising degree of morphological plasticity within the group.

Data from: Squamate Conserved Loci (SqCL): a unified set of conserved loci for phylogenomics and population genetics of squamate reptiles

The identification of conserved loci across genomes, along with advances in target capture methods and high-throughput sequencing, has helped spur a phylogenomics revolution by enabling researchers to gather large numbers of homologous loci across clades of interest with minimal upfront investment in locus design. Target capture for vertebrate animals is currently dominated by two approaches – anchored hybrid enrichment (AHE) and ultraconserved elements (UCE) – and both approaches have proven useful for addressing questions in phylogenomics, phylogeography, and population genomics. However, these two sets of loci have minimal overlap with each other; moreover, they do not include many traditional loci that that have been used for phylogenetics. Here, we combine across UCE, AHE, and traditional phylogenetic gene locus sets to generate the Squamate Conserved Loci (SqCL) set, a single integrated probe set that can generate high-quality and highly complete data across all three loci types. We use these probes to generate data for 44 phylogenetically-disparate taxa that collectively span approximately 33% of terrestrial vertebrate diversity. Our results generated an average of 4.29 Mb across 4709 loci per individual, of which an average of 2.99 Mb was sequenced to high enough coverage (≥10×) to use for population genetic analyses. We validate the utility of these loci for both phylogenomic and population genomic questions, provide a comparison among these locus sets of their relative usefulness, and suggest areas for future improvement.

opencc-zeroDec 2017View details →

Data from: Phylogenomics clarifies repeated evolutionary origins of inbreeding and fungus farming in bark beetles (Curculionidae, Scolytinae)

Bark and ambrosia beetles (Curculionidae, Scolytinae) display a conspicuous diversity of unusual genetic and ecological attributes and behaviors. Reconstructing the evolution of Scolytinae, particularly the large and ecologically significant tribe Cryphalini (pygmy borers), has long been problematic. These challenges have not adequately been addressed using morphological characters, and previous research has used only DNA sequence data from small numbers of genes. Through a combination of anchored hybrid enrichment, low-coverage draft genomes, and transcriptomes, we addressed these challenges by amassing a large molecular phylogenetic dataset for bark and ambrosia beetles. The resulting DNA sequence data from 251 protein coding genes (114,276 bp of nucleotide sequence data) support inference of the first robust phylogeny of Scolytinae, with a special focus on the species rich tribe Cryphalini and its close relatives. Key strategies, including inbreeding mating systems and fungus farming, evolved repeatedly across Scolytinae. We confirm 12 of 16 hypothesized origins of fungus farming, 6 of 8 origins of inbreeding polygyny and at least 11 independent origins of a super-generalist host range. These three innovations are statistically correlated, but their appearance within lineages was not necessarily simultaneous. Additionally, the evolution of extreme host plant generalism often preceded, rather than succeeded, fungus farming. Of the high-diversity tribes of Scolytinae, only Xyleborini is monophyletic, Corthylini is paraphyletic and Cryphalini is highly polyphyletic. Cryphalini sensu stricto is part of a clade containing the genera Hypothenemus, Cryphalus and Trypophloeus, and the tribe Xyloterini. Stegomerus and Cryptocarenus (Cryphalini) are part of a clade otherwise containing all Corthylini. Several other genera, including Ernoporus and Scolytogenes (Cryphalini), make up a distantly related clade. Several of the genera of Cryphalini are also intermixed. For example, Cryphalus and Hypocryphalus are intermingled, as well as Ernoporicus, Ptilopodius and Scolytogenes. Our data are consistent with widespread polyphyly and paraphyly across Scolytinae and within Cryphalini, and provides new insights into the evolution of inbreeding mating systems and fungus farming in the species rich and ecologically significant weevil subfamily Scolytinae.

opencc-zeroDec 2014View details →

Data from: The challenges of resolving a rapid, recent radiation: empirical and simulated phylogenomics of Philippine shrews

Phylogenetic relationships in recent, rapid radiations can be difficult to resolve due to incomplete lineage sorting and reliance on genetic markers that evolve slowly relative to the rate of speciation. By incorporating hundreds to thousands of unlinked loci, phylogenomic analyses have the potential to mitigate these difficulties. Here, we attempt to resolve phylogenetic relationships among eight shrew species (genus Crocidura) from the Philippines, a phylogenetic problem that has proven intractable with small (< 10 loci) data sets. We sequenced hundreds of ultraconserved elements and whole mitochondrial genomes in these species and estimated phylogenies using concatenation, summary coalescent, and hierarchical coalescent methods. The concatenated approach recovered a maximally supported and fully resolved tree. In contrast, the coalescent-based approaches produced similar topologies, but each had several poorly supported nodes. Using simulations, we demonstrate that the concatenated tree could be positively misleading. Our simulations also show that the tree shape we tend to infer, which involves a series of short internal branches, is difficult to resolve, even if substitution models are known and multiple individuals per species are sampled. As such, the low support we obtained for backbone relationships in our coalescent-based inferences reflects a real and appropriate lack of certainty. Our results illuminate the challenges of estimating a bifurcating tree in a rapid and recent radiation, providing a rare empirical example of a nearly simultaneous series of speciation events in a terrestrial animal lineage as it spreads across an oceanic archipelago.

opencc-zeroDec 2017View details →

Data from: Phylogenomics, life history and morphological evolution of ophiocomid brittlestars

Brittlestars in the family Ophiocomidae are large and colourful inhabitants of tropical shallow water habitats across the globe. Here we use targeted capture and next-generation sequencing to generate robust phylogenomic trees for 39 of the 43 species in order to test the monophyly of existing genera. The large genus Ophiocoma, as currently constituted, is paraphyletic on our trees and required revision. Four genera are recognised herein: an expanded Ophiomastix (now including Ophiocoma wendtii, O. occidentalis, O. endeani, O. macroplaca, and Ophiarthrum spp), Ophiocomella (now including the non-fissiparous Ophiocoma pumila, aethiops and valenciae) and Breviturma (now including Ophiocoma pica, O. pusilla, O. paucigranulata and O. longispina) and a restricted Ophiocoma. The resulting junior homonym Ophiomastix elegans is renamed O. brocki. The genus Ophiomastix exhibits relatively high rates of morphological disparity compared to other lineages. Ophiomastix flaccida and O. (formerly Ophiarthrum) pictum have divergent mitochondrial genomes, characterised by gene-order rearrangements, strand recoding, enriched GT base composition, and a corresponding divergence of nuclear mitochondrial protein genes. The new phylogeny indicates that larval and developmental transitions occurred rarely. Larval culture trials show that species with abbreviated lecithotrophic larval development occur only within Ophiomastix, although the possible monophyly of these species is obscured by the rapid early radiation within this genus. Asexual reproduction by fission is limited to one species-complex within Ophiocomella, also characterised by elevated levels of allelic heterozygosity, and which has achieved a relatively rapid global distribution. The crown ages of the new genera considerably predate the closure of the Tethyan seaway and all four are distributed in both the Atlantic and Indo-Pacific Oceans. Two species pairs appear to reflect the closure of the Panama Seaway, although their fossil-calibrated node ages (12–14 ± 6 my), derived from both concatenated sequence and multispecies coalescent analyses, considerably predate the terminal closure event. Ophiocoma erinaceus has crossed the East Pacific barrier and is recorded from Clipperton Island, SW of Mexico.

Data from: Restructuring higher taxonomy using broad-scale phylogenomics: the living Ophiuroidea

The power and throughput of next-generation sequencing is instigating a major transformation in our understanding of evolution and classification of life on our planet. The new trees of life are robust and comprehensive. Here we provide a landmark phylogeny of the living ophiuroids and use it as the basis for a major revision of the higher classification of this class of marine invertebrates. We used an exon-capture system to generate a 1484 exon (273 kbp) data-matrix from DNA extracted from ethanol-preserved museum samples. We successfully obtained an average of 90% of our target sequence from 576 species spread across the known taxonomic diversity. The topology of the major lineages was robust to taxon sampling, exon-sampling, models and methods. However, estimates of node age were much less precise, varying by about a quarter of mean age. We used a combination of phylogenetic distinctiveness and temporal-banding to guide our revision of the family-level classification. Empirically, we determined that limiting family crown age to 110 ± 10 Ma (mid Cretaceous) selected phylogenetically distinct nodes while minimising disruption to the existing taxonomy. The resulting scheme of 32 families and six orders considerably expands the number of higher taxa. The families are generally longitudinally widespread across the world's oceans, although 17 are largely confined to temperate and equatorial latitudes and six to relatively shallow water (less than 1000 m depth).

Data from: Phylogenomic analysis of a rapid radiation of misfit fishes (Syngnathiformes) using ultraconserved elements

Phylogenetics is undergoing a revolution as large-scale molecular datasets reveal unexpected but repeatable rearrangements of clades that were previously thought to be disparate lineages. One of the most unusual clades of fishes that has been found using large-scale molecular datasets is an expanded Syngnathiformes including traditional long-snouted syngnathiform lineages (Aulostomidae, Centriscidae, Fistulariidae, Solenostomidae, Syngnathidae), as well as a diverse set of largely benthic-associated fishes (Callionymoidei, Dactylopteridae, Mullidae, Pegasidae) that were previously dispersed across three orders. The monophyly of this surprising clade of fishes has been upheld by recent studies utilizing both nuclear and mitogenomic data, but the relationships among major lineages within Syngnathiformes remain ambiguous; previous analyses have inconsistent topologies and are plagued by low support at deep divergences between the major lineages. In this study, we use a dataset of ultraconserved elements (UCEs) to conduct the first phylogenomic study of Syngnathiformes. UCEs have been effective markers for resolving deep phylogenetic relationships in fishes and, combined with increased taxon sampling, we expected UCEs to resolve problematic syngnathiform relationships. Overall, UCEs were effective at resolving relationships within Syngnathiformes at a range of evolutionary timescales. We find consistent support for the monophyly of traditional long-snouted syngnathiform lineages (Aulostomidae, Centriscidae, Fistulariidae, Solenostomidae, Syngnathidae), which better agrees with morphological hypotheses than previously published topologies from molecular data. This result was supported by all Bayesian and maximum likelihood analyses, was robust to differences in matrix completeness and potential sources of bias, and was highly supported in coalescent-based analyses in ASTRAL when matrices were filtered to contain the most phylogenetically informative loci. While Bayesian and maximum likelihood analyses found support for a benthic-associated clade (Callionymidae, Dactylopteridae, Mullidae, and Pegasidae) as sister to the long-snouted clade, this result was not replicated in the ASTRAL analyses. The base of our phylogeny is characterized by short internodes separating major syngnathiform lineages and is consistent with the hypothesis of an ancient rapid radiation at the base of Syngnathiformes. Syngnathiformes therefore present an exciting opportunity to study patterns of morphological variation and functional innovation arising from rapid but ancient radiation.

opencc-zeroDec 2013View details →

Data from: Extensive introgression in a malaria vector species complex revealed by phylogenomics

Introgressive hybridization is now recognized as a widespread phenomenon, but its role in evolution remains contested. Here we use newly available reference genome assemblies to investigate phylogenetic relationships and introgression in a medically important group of Afrotropical mosquito sibling species. We have identified the correct species branching order to resolve a contentious phylogeny, and show that lineages leading to the principal vectors of human malaria were among the first to split. Pervasive autosomal introgression between these malaria vectors means that only a small fraction of the genome, mainly on the X chromosome, has not crossed species boundaries. Our results suggest that traits enhancing vectorial capacity may be gained through interspecific gene flow, including between non-sister species.

Data from: Whole genome shotgun phylogenomics resolves the pattern and timing of swallowtail butterfly evolution

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 non-biological 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 6,621 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 dataset (760 orthologs and 12% missing data) indicate a mid-Cretaceous origin of Papilionoidea around 99.2 million years ago (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.

Data from: Conflicting phylogenomic signals reveal a pattern of reticulate evolution in a recent high-Andean diversification (Asteraceae: Astereae: Diplostephium)

High-throughput sequencing is helping biologists to overcome the difficulties of inferring the phylogenies of recently diverged taxa. The present study analyzes the phylogenetic signal of genomic regions with different inheritance patterns using genome skimming and ddRAD-seq in a species-rich Andean genus (Diplostephium) and its allies. We analyzed the complete nuclear ribosomal cistron, the complete chloroplast genome, a partial mitochondrial genome, and a nuclear-ddRAD matrix separately with phylogenetic methods. We applied several approaches to understand the causes of incongruence among datasets, including simulations and the detection of introgression using the D-statistic (ABBA-BABA test). We found significant incongruence among the nuclear, chloroplast, and mitochondrial phylogenies. The strong signal of hybridization found by simulations and the D-statistic among genera and inside the main clades of Diplostephium indicate reticulate evolution as a main cause of phylogenetic incongruence. Our results add evidence for a major role of reticulate evolution in events of rapid diversification. Hybridization and introgression confound chloroplast and mitochondrial phylogenies in relation to the species tree as a result of the uniparental inheritance of these genomic regions. Practical implications regarding the prevalence of hybridization are discussed in relation to the phylogenetic method.

opencc-zeroDec 2014View details →

Data from: Phylogenomics of horned lizards (genus: Phrynosoma) using targeted sequence capture data

New genome sequencing techniques are enabling phylogenetic studies to scale-up from using a handful of loci to hundreds or thousands of loci from throughout the genome. In this study, we use targeted sequence capture (TSC) data from 540 ultraconserved elements and 44 protein-coding genes to estimate the phylogenetic relationships among all 17 species of horned lizards in the genus Phrynosoma. Previous molecular phylogenetic analyses of Phrynosoma based on a few nuclear genes, restriction site associated DNA (RAD) sequencing, or mitochondrial DNA (mtDNA) have produced conflicting relationships. Some of these conflicts are likely the result of rapid speciation at the start of Phrynosoma diversification, whereas other examples of gene tree discordance appear to be caused by active and residual traces of hybridization. Concatenation and coalescent-based species tree phylogenetic analyses of these new TSC data support the same topology, and a divergence dating analysis suggests that the Phrynosoma crown group is up to 30 million years old. The new phylogenomic tree supports the recognition of four main clades within Phrynosoma, including Anota (P. mcallii, P. solare, and the P. coronatum complex), Doliosaurus (P. modestum, P. goodei, and P. platyrhinos), Tapaja (P. ditmarsi, P. douglasii, P. hernandesi, and P. orbiculare), and Brevicauda (P. braconnieri, P. sherbrookei, and P. taurus). The phylogeny provides strong support for the relationships among all species of Phrynosoma and provides a robust new framework for conducting comparative analyses.

Data from: First plastid phylogenomic study reveals potential cyto-nuclear discordance in the evolutionary history of Ficus L. (Moraceae)

Standard Sanger chloroplast markers provide limited information to resolve species level relationships within plants, in particular within large genera. Figs (Ficus L., Moraceae) compose one of the 50 largest genera of angiosperms with ∼750 species occurring in the tropics and subtropics worldwide. Figs, in addition to being a keystone food resource in rainforests, are well-known for the mutualistic interactions with their pollinating wasps. It is regarded as a model system for understanding co-evolution dating back more than 75 million years. However, despite significant taxon sampling, combinations of low copy nuclear, nuclear ribosomal and chloroplast regions have not been able to confidently resolve relationships among major groups of figs. Using a high throughput sequencing approach we attempted to resolve the major lineages of Ficus based on plastome data. In this study, we show that the use of a de novo assembled plastome from within the genus provides less ambiguity and higher coverage across the 59 Ficus and 6 outgroup plastome assemblies compared to using the nearest available reference plastome outside the genus resulting in improved resolution and higher support of the phylogenetic relationships within Ficus inferred from plastome data. Chloroplast genome data confidently resolved relationships among major groups of figs and largely support current understanding based on nuclear sequence data including passively pollinated Neotropical section Pharmacosycea as sister lineage to all other Ficus. However, conflicts between the new plastome topology and previous nuclear studies are observed for both individual species as well as relationships among some sections at deeper levels. Conflicts could be caused by lack of resolution in the nuclear data or may indicate potential cyto-nuclear discordance as previously observed in an African lineage of Ficus.

opencc-zeroJul 2019View details →

Data from: Chloroplast phylogenomic data support Paleocene - Eocene amphi-Pacific early radiation for the Asian Palmate core Araliaceae

Traditional phylogenies based on analysis of multiple genes have failed to obtain a well-resolved evolutionary history for the backbone of the Asian Palmate group of Araliaceae, the largest clade of the family. In this study we applied the genome skimming approach of next-generation sequencing to address whether the lack of resolution at the base of the Asian Palmate tree is due to molecular sampling error or the footprint of an ancient radiation. Twenty-nine complete plastid genomes of Araliaceae (17 newly sequenced) were analyzed (RAxML, Beast, Lagrange, BioGeoBears) to provide the first phylogenomic reconstruction of the group (95% of genera included). As a result, the early divergences of the Asian Palmate group have been clarified but the backbone of its core is not totally resolved, with short internal branches pointing to an ancient radiation scenario. East Asia is inferred as the most likely ancestral area for the Asian Palmate group (early Upper Cretaceous and the Paleocene) from which an early colonization of the Neotropics is inferred during the late Paleocene – early Eocene. The radiation of the core took place during the Paleocene – early Eocene, most likely in the context of the Boreotropical hypothesis. Recurrent episodes of southward migration (to the tropics) coupled with northern latitude local extinctions (promoting lineages geographic isolation) and northward expansion (promoting lineages contact that erased the trace of preceding geographic isolation) are hypothesized to have linked to the alternation of the cold and warm periods of the Eocene.

Data from: A phylogenomic rodent tree reveals the repeated evolution of masseter architectures

Understanding the number of times a trait has evolved is a necessary foundation for comprehending its potential relationships with selective regimes, developmental constraints, and evolutionary diversification. Rodents make up over 40% extant of mammalian species and their ecological and evolutionary success has been partially attributed to the increase in biting efficiency that resulted from a forward shift of one or two portions of the masseter muscle from the zygomatic arch onto the rostrum. This forward shift has occurred in three discrete ways, but the number of times it has occurred has never been explicitly quantified. We estimated an ultrametric phylogeny, the first to include all rodent families, using thousands of ultraconserved elements. We examined support for evolutionary relationships among the five rodent suborders and then incorporated relevant fossils, fit models of character evolution, and used stochastic character mapping to determine that a portion of the masseter muscle has moved forward onto the rostrum at least seven times (with one reversal) during the ca. 70 MY history of rodents. Combined, the repeated evolution of this key innovation, its increasing prevalence through time, and the species diversity of clades with this character underscores the adaptive value of improved biting efficiency and the relative ease with which some advantageous traits arise.

opencc-zeroDec 2014View details →

Data from: Phylogenomic analyses resolve an ancient trichotomy at the base of Ischyropsalidoidea (Arachnida, Opiliones) despite high levels of gene tree conflict and unequal minority resolution frequencies

Phylogenetic resolution of ancient rapid radiations has remained problematic despite major advances in statistical approaches and DNA sequencing technologies. Here we report on a combined phylogenetic approach utilizing transcriptome data in conjunction with Sanger sequence data to investigate a tandem of ancient divergences in the harvestmen superfamily Ischyropsalidoidea (Arachnida, Opiliones, Dyspnoi). We rely on Sanger sequences to resolve nodes within and between closely related genera, and use RNA-seq data from a subset of taxa to resolve a short and ancient internal branch. We use several analytical approaches to explore this succession of ancient diversification events, including concatenated and coalescent-based analyses and maximum likelihood gene trees for each locus. We evaluate the robustness of phylogenetic inferences using a randomized locus sub-sampling approach, and find congruence across these methods despite considerable incongruence across gene trees. Incongruent gene trees are not recovered in frequencies expected from a simple multispecies coalescent model, and we reject incomplete lineage sorting as the sole contributor to gene tree conflict. Using these approaches we attain robust support for higher-level phylogenetic relationships within Ischyropsalidoidea.

Data from: Assessing combinability of phylogenomic data using Bayes Factors

With the rapid reduction in sequencing costs of high-throughput genomic data, it has become commonplace to use hundreds of genes to infer phylogeny of any study system. While sampling a large number of genes has given us a tremendous opportunity to uncover previously unknown relationships and improve phylogenetic resolution, it also presents us with new challenges when the phylogenetic signal is confused by differences in the evolutionary histories of sampled genes. Given the incorporation of accurate marginal likelihood estimation methods into popular Bayesian software programs, it is natural to consider using the Bayes Factor (BF) to compare different partition models in which genes within any given partition subset share both tree topology and edge lengths. We explore using marginal likelihood to assess data subset combinability when data subsets have varying levels of phylogenetic discordance due to deep coalescence events among genes (simulated within a species tree), and compare the results with our recently-described phylogenetic informational dissonance index (D) estimated for each data set. BF effectively detects phylogenetic incongruence, and provides a way to assess the statistical significance of D values. We use BFs to assess data combinability using an empirical data set comprising 56 plastid genes from the green algal order Volvocales. We also discuss the potential need for calibrating BFs and demonstrate that BFs used in this study are correctly calibrated.

Data from: Conserved non-exonic elements: a novel class of marker for phylogenomics

Noncoding markers have a particular appeal as tools for phylogenomic analysis because, at least in vertebrates, they appear less subject to strong variation in GC content among lineages. Thus far, ultraconserved elements (UCEs) and introns have been the most widely used noncoding markers. Here we analyze and study the evolutionary properties of a new type of noncoding marker, conserved non-exonic elements (CNEEs), which consists of noncoding elements that are estimated to evolve slower than the neutral rate across a set of species. Although they often include UCEs, CNEEs are distinct from UCEs because they are not ultraconserved, and, most importantly, the core region alone is analyzed, rather than both the core and its flanking regions. Using a data set of 16 birds plus an alligator outgroup, and ∼3600 - ∼3800 loci per marker type, we found that although CNEEs were less variable than bioinformatically-derived UCEs or introns and in some cases exhibited a slower approach to branch resolution as determined by phylogenomic subsampling, the quality of CNEE alignments was superior to those of the other markers, with fewer gaps and missing species. Phylogenetic resolution using coalescent approaches was comparable among the three marker types, with most nodes being fully and congruently resolved. Comparison of phylogenetic results across the three marker types indicated that one branch, the sister group to the passerine+falcon clade, was resolved differently and with moderate (> 70%) bootstrap support between CNEEs and UCEs or introns. Overall, CNEEs appear to be promising as phylogenomic markers, yielding phylogenetic resolution as high as for UCEs and introns but with fewer gaps, less ambiguity in alignments and with patterns of nucleotide substitution more consistent with the assumptions of commonly used methods of phylogenetic analysis.