Ultraconserved Elements Data Research Articles

Phylogenomic Species Delimitation, Taxonomy, and ‘Bird Guide’ Identification for the Neotropical Ant Genus Rasopone (Hymenoptera: Formicidae)

Abstract Rasopone Schmidt and Shattuck is a poorly known lineage of ants that live in Neotropical forests. Informed by phylogenetic results from thousands of ultraconserved elements (UCEs) and mitochondrial DNA barcodes, we revise the genus, providing a new morphological diagnosis and a species-level treatment. Analysis of UCE data from many Rasopone samples and select outgroups revealed non-monophyly of the genus. Monophyly of Rasopone was restored by transferring several species to the unrelated genus Mayaponera Schmidt and Shattuck. Within Rasopone, species are morphologically very similar, and we provide a ‘bird guide’ approach to identification rather than the traditional dichotomous key. Species are arranged by size in a table, along with geographic range and standard images. Additional diagnostic information is then provided in individual species accounts. We recognize a total of 15 named species, of which the following are described as new species: R. costaricensis, R. cryptergates, R. cubitalis, R. guatemalensis, R. mesoamericana, R. pluviselva, R. politognatha, R. subcubitalis, and R. titanis. An additional 12 morphospecies are described but not formally named due to insufficient material. Rasopone panamensis (Forel, 1899) is removed from synonymy and elevated to species. The following species are removed from Rasopone and made new combinations in Mayaponera: M. arhuaca (Forel, 1901), M. becculata (Mackay and Mackay, 2010), M. cernua (Mackay and Mackay, 2010), M. conicula (Mackay and Mackay, 2010), M. longidentata (Mackay and Mackay, 2010), and M. pergandei (Forel, 1909).

Spider‐specific probe set for ultraconserved elements offers new perspectives on the evolutionary history of spiders (Arachnida, Araneae)

Phylogenomic methods have proven useful for resolving deep nodes and recalcitrant groups in the spider tree of life. Across arachnids, transcriptomic approaches may generate thousands of loci, and target-capture methods, using the previously designed arachnid-specific probe set, can target a maximum of about 1,000 loci. Here, we develop a specialized target-capture probe set for spiders that contains over 2,000 ultraconserved elements (UCEs) and then demonstrate the utility of this probe set through sequencing and phylogenetic analysis. We designed the 'spider-specific' probe set using three spider genomes (Loxosceles, Parasteatoda and Stegodyphus) and ensured that the newly designed probe set includes UCEs from the previously designed Arachnida probe set. The new 'spider-specific' probes were used to sequence UCE loci in 51 specimens. The remaining samples included five spider genomes and taxa that were enriched using Arachnida probe set. The 'spider-specific' probes were also used to gather loci from a total of 84 representative taxa across Araneae. On mapping these 84 taxa to the Arachnida probe set, we captured at most 710 UCE loci, while the spider-specific probe set captured up to 1,547 UCE loci from the same taxon sample. Phylogenetic analyses using maximum likelihood and coalescent methods corroborate most nodes resolved by recent transcriptomic analyses, but not all (e.g. UCE data suggest monophyly of 'symphytognathoids'). Our preferred hypothesis based on topology tests, suggests monophyly of the 'symphytognathoids' (the miniature orb weavers), which in previous studies has only been supported by a combination of morphological and behavioural characters.

Phylogenomic analysis suggests Coreidae and Alydidae (Hemiptera: Heteroptera) are not monophyletic

AbstractNext‐generation sequencing technologies (NGS) allow systematists to amass a wealth of genomic data from non‐model species for phylogenetic resolution at various temporal scales. However, phylogenetic inference for many lineages dominated by non‐model species has not yet benefited from NGS, which can complement Sanger sequencing studies. One such lineage, whose phylogenetic relationships remain uncertain, is the diverse, agriculturally important and charismatic Coreoidea (Hemiptera: Heteroptera). Given the lack of consensus on higher‐level relationships and the importance of a robust phylogeny for evolutionary hypothesis testing, we use a large data set comprised of hundreds of ultraconserved element (UCE) loci to infer the phylogeny of Coreoidea (excluding Stenocephalidae and Hyocephalidae), with emphasis on the families Coreidae and Alydidae. We generated three data sets by including alignments that contained loci sampled for at least 50%, 60%, or 70% of the total taxa, and inferred phylogeny using maximum likelihood and summary coalescent methods. Twenty‐six external morphological features used in relatively comprehensive phylogenetic analyses of coreoids were also re‐evaluated within our molecular phylogenetic framework. We recovered 439–970 loci per species (16%–36% of loci targeted) and combined this with previously generated UCE data for 12 taxa. All data sets, regardless of analytical approach, yielded topologically similar and strongly supported trees, with the exception of outgroup relationships and the position of Hydarinae. We recovered a monophyletic Coreoidea, with Rhopalidae highly supported as the sister group to Alydidae + Coreidae. Neither Alydidae nor Coreidae were monophyletic; the coreid subfamilies Hydarinae and Pseudophloeinae were recovered as more closely related to Alydidae than to other coreid subfamilies. Coreinae were paraphyletic with respect to Meropachyinae. Most morphological traits were homoplastic with several clades defined by few, if any, synapomorphies. Our results demonstrate the utility of phylogenomic approaches in generating robust hypotheses for taxa with long‐standing phylogenetic problems and highlight that novel insights may come from such approaches.

Insight from an ultraconserved element bait set designed for hemipteran phylogenetics integrated with genomic resources

Target enrichment of conserved genomic regions facilitates collecting sequences of many orthologous loci from non-model organisms to address phylogenetic, phylogeographic, population genetic, and molecular evolution questions. Bait sets for sequence capture can simultaneously target thousands of loci, which opens new avenues of research on speciose groups. Current phylogenetic hypotheses on the >103,000 species of Hemiptera have failed to unambiguously resolve major nodes, suggesting that alternative datasets and more thorough taxon sampling may be required to resolve relationships. We use a recently designed ultraconserved element (UCE) bait set for Hemiptera, with a focus on the suborder Heteroptera, or the true bugs, to test previously proposed relationships. We present newly generated UCE data for 36 samples representing three suborders, all seven heteropteran infraorders, 23 families, and 34 genera of Hemiptera and one thysanopteran outgroup. To improve taxon sampling, we also mined additional UCE loci in silico from published hemipteran genomic and transcriptomic data. We obtained 2271 UCE loci for newly sequenced hemipteran taxa, ranging from 265 to 1696 (average 904) per sample. These were similar in number to the data mined from transcriptomes and genomes, but with fewer loci overall. The amount of missing data correlates with greater phylogenetic divergence from taxa used to design the baits. This bait set hybridizes to a wide range of hemipteran taxa and specimens of varying quality, including dried specimens as old as 1973. Our estimated phylogeny yielded topologies consistent with other studies for most nodes and was strongly-supported. We also demonstrate that UCE loci are almost exclusively from the transcribed portion of the genome, thus data can be successfully integrated with existing genomic and transcriptomic resources for more comprehensive phylogenetic sampling, an important feature in the era of phylogenomics. UCE approaches can be used by other researchers for additional studies on hemipteran evolution and other research that requires well resolved phylogenies.

Combining transcriptomes and ultraconserved elements to illuminate the phylogeny of Apidae

Two increasingly popular approaches to reconstruct the Tree of Life involve whole transcriptome sequencing and the target capture of ultraconserved elements (UCEs). Both methods can be used to generate large, multigene datasets for analysis of phylogenetic relationships in non-model organisms. While targeted exon sequencing across divergent lineages is now a standard method, it is still not clear if UCE data can be readily combined with published transcriptomes. In this study, we evaluate the combination of UCEs and transcriptomes in a single analysis using genome-, transcriptome-, and UCE data for 79 bees in the largest and most biologically diverse bee family, Apidae. Using existing tools, we first developed a workflow to assemble phylogenomic data from different sources and produced two large nucleotide matrices of combined data. We then reconstructed the phylogeny of the Apidae using concatenation- and coalescent-based methods, and critically evaluated the resulting phylogenies in the context of previously published genetic, genomic, and morphological data sets. Our estimated phylogenetic trees are robustly supported and largely congruent with previous molecular hypotheses, from deep nodes to shallow species-level phylogenies. Moreover, the combined approach allows us to resolve controversial nodes of the apid Tree of Life, by clarifying the relationships among the genera of orchid bees (Euglossini) and the monophyly of the Centridini. Additionally, we present novel phylogenetic evidence supporting the monophyly of the diverse clade of cleptoparasitic Apidae and the placement of two enigmatic, oil-collecting genera (Ctenoplectra and Tetrapedia). Lastly, we propose a revised classification of the family Apidae that reflects our improved understanding of apid higher-level relationships.

Ultraconserved elements (UCEs) illuminate the population genomics of a recent, high-latitude avian speciation event.

Using a large, consistent set of loci shared by descent (orthologous) to study relationships among taxa would revolutionize among-lineage comparisons of divergence and speciation processes. Ultraconserved elements (UCEs), highly conserved regions of the genome, offer such genomic markers. The utility of UCEs for deep phylogenetics is clearly established and there are mature analytical frameworks available, but fewer studies apply UCEs to recent evolutionary events, creating a need for additional example datasets and analytical approaches. We used UCEs to study population genomics in snow and McKay’s buntings (Plectrophenax nivalis and P. hyperboreus). Prior work suggested divergence of these sister species during the last glacial maximum (∼18–74 Kya). With a sequencing depth of ∼30× from four individuals of each species, we used a series of analysis tools to genotype both alleles, obtaining a complete dataset of 2,635 variable loci (∼3.6 single nucleotide polymorphisms/locus) and 796 invariable loci. We found no fixed allelic differences between the lineages, and few loci had large allele frequency differences. Nevertheless, individuals were 100% diagnosable to species, and the two taxa were different genetically (FST = 0.034; P = 0.03). The demographic model best fitting the data was one of divergence with gene flow. Estimates of demographic parameters differed from published mtDNA research, with UCE data suggesting lower effective population sizes (∼92,500–240,500 individuals), a deeper divergence time (∼241,000 years), and lower gene flow (2.8–5.2 individuals per generation). Our methods provide a framework for future population studies using UCEs, and our results provide additional evidence that UCEs are useful for answering questions at shallow evolutionary depths.

Sequence capture phylogenomics of eyeless Cicurina spiders from Texas caves, with emphasis on US federally-endangered species from Bexar County (Araneae, Hahniidae).

Morphological, mitochondrial, and nuclear phylogenomic data were combined to address phylogenetic and species delimitation questions in cave-limited Cicurina spiders from central Texas. Special effort was focused on specimens and cave locations in the San Antonio region (Bexar County), home to four eyeless species listed as US Federally Endangered. Sequence capture experiments resulted in the recovery of ~200–400 homologous ultra-conserved element (UCE) nuclear loci across taxa, and nearly complete COI mitochondrial DNA sequences from the same set of individuals. Some of these nuclear and mitochondrial sequences were recovered from “standard” museum specimens without special preservation of DNA material, including museum specimens preserved in the 1990s. Multiple phylogenetic analyses of the UCE data agree in the recovery of two major lineages of eyeless Cicurina in Texas. These lineages also differ in mitochondrial clade membership, female genitalic morphology, degree of troglomorphy (as measured by relative leg length), and are mostly allopatric across much of Texas. Rare sympatry was confirmed in Bexar County, where members of the two major clades sometimes co-exist in the same karst feature. Both nuclear phylogenomic and mitochondrial data indicate the existence of undescribed species from the San Antonio region, although further sampling and collection of adult specimens is needed to explicitly test these hypotheses. Our data support the two following species synonymies (Cicurina venii Gertsch, 1992 = Cicurina madla Gertsch, 1992; Cicurina loftini Cokendolpher, 2004 = Cicurina vespera Gertsch, 1992), formally proposed here. Overall, our taxonomy-focused research has many important conservation implications, and again highlights the fundamental importance of robust taxonomy in conservation research.

Phenotypic and Genetic Structure Support Gene Flow Generating Gene Tree Discordances in an Amazonian Floodplain Endemic Species.

Before populations become independent evolutionary lineages, the effects of micro evolutionary processes tend to generate complex scenarios of diversification that may affect phylogenetic reconstruction. Not accounting for gene flow in species tree estimates can directly impact topology, effective population sizes and branch lengths, and the resulting estimation errors are still poorly understood in wild populations. In this study, we used an integrative approach, including sequence capture of ultra-conserved elements (UCEs), mtDNA Sanger sequencing and morphological data to investigate species limits and phylogenetic relationships in face of gene flow in an Amazonian endemic species (Myrmoborus lugubris: Aves). We used commonly implemented species tree and model-based approaches to understand the potential effects of gene flow in phylogenetic reconstructions. The genetic structure observed was congruent with the four recognized subspecies of M. lugubris. Morphological and UCEs data supported the presence of a wide hybrid zone between M. l. femininus from the Madeira river and M. l. lugubris from the Middle and lower Amazon river, which were recovered as sister taxa by species tree methods. When fitting gene flow into simulated demographic models with different topologies, the best-fit model indicated these two taxa as non-sister lineages, a finding that is in agreement with the results of mitochondrial and morphological analyses. Our results demonstrated that failing to account for gene flow when estimating phylogenies at shallow divergence levels can generate topological uncertainty, which can nevertheless be statistically well supported, and that model testing approaches using simulated data can be useful tools to test alternative phylogenetic hypotheses.

Conflicting Evolutionary Histories of the Mitochondrial and Nuclear Genomes in New World Myotis Bats.

The rapid diversification of Myotis bats into more than 100 species is one of the most extensive mammalian radiations available for study. Efforts to understand relationships within Myotis have primarily utilized mitochondrial markers and trees inferred from nuclear markers lacked resolution. Our current understanding of relationships within Myotis is therefore biased towards a set of phylogenetic markers that may not reflect the history of the nuclear genome. To resolve this, we sequenced the full mitochondrial genomes of 37 representative Myotis, primarily from the New World, in conjunction with targeted sequencing of 3648 ultraconserved elements (UCEs). We inferred the phylogeny and explored the effects of concatenation and summary phylogenetic methods, as well as combinations of markers based on informativeness or levels of missing data, on our results. Of the 294 phylogenies generated from the nuclear UCE data, all are significantly different from phylogenies inferred using mitochondrial genomes. Even within the nuclear data, quartet frequencies indicate that around half of all UCE loci conflict with the estimated species tree. Several factors can drive such conflict, including incomplete lineage sorting, introgressive hybridization, or even phylogenetic error. Despite the degree of discordance between nuclear UCE loci and the mitochondrial genome and among UCE loci themselves, the most common nuclear topology is recovered in one quarter of all analyses with strong nodal support. Based on these results, we re-examine the evolutionary history of Myotis to better understand the phenomena driving their unique nuclear, mitochondrial, and biogeographic histories.

Phylogenomic analyses of 539 highly informative loci dates a fully resolved time tree for the major clades of living turtles (Testudines)

Accurate time-calibrated phylogenies are the centerpiece of many macroevolutionary studies, and the relationship between the size and scale of molecular data sets and the density and accuracy of fossil calibrations is a key element of time tree studies. Here, we develop a target capture array specifically for living turtles, compare its efficiency to an ultraconserved element (UCE) dataset, and present a time-calibrated molecular phylogeny based on 539 nuclear loci sequenced from 26 species representing the breadth of living turtle diversity plus outgroups. Our gene array, based on three fully sequenced turtle genomes, is 2.4 times more variable across turtles than a recently published UCE data set for an identical subset of 13 species, confirming that taxon-specific arrays return more informative data per sequencing effort than UCEs. We used our genomic data to estimate the ages of living turtle clades including a mid-late Triassic origin for crown turtles and a mid-Carboniferous split of turtles from their sister group, Archosauria. By specifically excluding several of the earliest potential crown turtle fossils and limiting the age of fossil calibration points to the unambiguous crown lineage Caribemys oxfordiensis from the Late Jurassic (Oxfordian, 163.5–157.3Ma) we corroborate a relatively ancient age for living turtles. We also provide novel age estimates for five of the ten testudine families containing more than a single species, as well as several intrafamilial clades. Most of the diversity of crown turtles appears to date to the Paleogene, well after the Cretaceous-Paleogene mass extinction 66mya.

Ultraconserved elements show utility in phylogenetic inference of Adephaga (Coleoptera) and suggest paraphyly of ‘Hydradephaga’

AbstractThe beetle suborder Adephaga has been the subject of many phylogenetic reconstructions utilizing a variety of data sources and inference methods. However, no strong consensus has yet emerged on the relationships among major adephagan lineages. Ultraconserved elements (UCEs) have proved useful for inferring difficult or unresolved phylogenies at varying timescales in vertebrates, arachnids and Hymenoptera. Recently, a UCE bait set was developed for Coleoptera using polyphagan genomes and a member of the order Strepsiptera as an outgroup. Here, we examine the utility of UCEs for reconstructing the phylogeny of adephagan families, in the first in vitro application a UCE bait set in Coleoptera. Our final dataset included 305 UCE loci for 18 representatives of all adephagan families except Aspidytidae, and two polyphagan outgroups, with a total concatenated length of 83 547 bp. We inferred trees using maximum likelihood analyses of the concatenated UCE alignment and coalescent species tree methods (astral ii, ASTRID, svdquartets). Although the coalescent species tree methods had poor resolution and weak support, concatenated analyses produced well‐resolved, highly supported trees. Hydradephaga was recovered as paraphyletic, with Gyrinidae sister to Geadephaga and all other adephagans. Haliplidae was recovered as sister to Dytiscoidea, with Hygrobiidae and Amphizoidae successive sisters to Dytiscidae. Finally, Noteridae was recovered as monophyletic and sister to Meruidae. Given the success of UCE data for resolving phylogenetic relationships within Adephaga, we suggest the potential for further resolution of relationships within Adephaga using UCEs with improved taxon sampling, and by developing Adephaga‐specific probes.

The impact of GC bias on phylogenetic accuracy using targeted enrichment phylogenomic data

The field of sequence based phylogenetic analyses is currently being transformed by novel hybrid-based targeted enrichment methods, such as the use of ultraconserved elements (UCEs). Rather than analyzing relationships among organisms using a small number of genes, these methods now allow us to evaluate relationships with many hundreds to thousands of individual gene loci. However, the inclusion of thousands of loci does not necessarily overcome the long-standing challenge of incongruence among phylogenetic trees derived from different genes or gene regions. One factor that impacts the level of incongruence in phylogenomic data sets is the level of GC bias. GC rich gene regions are prone to higher recombination rates than AT rich regions, driven by a process referred to as “GC biased gene conversion”. As a result, high GC content can be negatively associated with phylogenetic accuracy, but the extent to which this impacts incongruence among UCEs is currently unstudied. We investigated the impact of GC content on phylogeny reconstruction using in silico captured UCE data for the corbiculate bees (Hymenoptera: Apidae). The phylogeny of this group has been the subject of extensive study, and incongruence among gene trees is thought to be a source of phylogenetic error. We conducted coalescent- and concatenation-based analyses of 810 individual gene loci from all 13 currently available bee genomes, including 8 corbiculate taxa. Both coalescent- and concatenation-based methods converged on a single topology for the corbiculate tribes. In contrast to concatenation, the coalescent-based methods revealed significant topological conflict at nodes involving the orchid bees (Euglossini) and honeybees (Apini). Partitioning the loci by GC content reveals decreasing support for the inferred topology with increasing GC bias. Based on the results of this study, we report the first evidence that GC biased gene conversion may contribute to topological incongruence in studies based on ultraconserved elements.

Phylogenomic Systematics of Ostariophysan Fishes: Ultraconserved Elements Support the Surprising Non-Monophyly of Characiformes.

Ostariophysi is a superorder of bony fishes including more than 10,300 species in 1100 genera and 70 families. This superorder is traditionally divided into five major groups (orders): Gonorynchiformes (milkfishes and sandfishes), Cypriniformes (carps and minnows), Characiformes (tetras and their allies), Siluriformes (catfishes), and Gymnotiformes (electric knifefishes). Unambiguous resolution of the relationships among these lineages remains elusive, with previous molecular and morphological analyses failing to produce a consensus phylogeny. In this study, we use over 350 ultraconserved element (UCEs) loci comprising 5 million base pairs collected across 35 representative ostariophysan species to compile one of the most data-rich phylogenies of fishes to date. We use these data to infer higher level (interordinal) relationships among ostariophysan fishes, focusing on the monophyly of the Characiformes-one of the most contentiously debated groups in fish systematics. As with most previous molecular studies, we recover a non-monophyletic Characiformes with the two monophyletic suborders, Citharinoidei and Characoidei, more closely related to other ostariophysan clades than to each other. We also explore incongruence between results from different UCE data sets, issues of orthology, and the use of morphological characters in combination with our molecular data. [Conserved sequence; ichthyology; massively parallel sequencing; morphology; next-generation sequencing; UCEs.].

Genomic data reveals potential for hybridization, introgression, and incomplete lineage sorting to confound phylogenetic relationships in an adaptive radiation of narrow-mouth frogs.

The microhylid frog genus Kaloula is an adaptive radiation spanning the edge of the Asian mainland and multiple adjacent island archipelagos, with much of the clade's diversity associated with an endemic Philippine radiation. Relationships among clades from the Philippines, however, remain unresolved. With ultraconserved element (UCE) and mitogenomic data, we identified highly supported differences in topology and areas of poor resolution, for each marker set. Using the UCE data, we then identified possible instances of contemporary hybridization, past introgression, and incomplete lineage sorting (ILS) within the Philippine Kaloula. Using a simulation approach, and an estimate of the Philippine Kaloula clade origin (12.7-21.0 mya), we demonstrate that an evolutionary history including inferred instances of hybridization, introgression, and ILS leads to phylogenetic reconstructions that show concordance with results from the observed mitogenome and UCE data. In the process of validating a complex evolutionary scenario in the Philippine Kaloula, we provide the first demonstration of the efficacy of UCE data for phylogenomic studies of anuran amphibians.

Sequence Capture and Phylogenetic Utility of Genomic Ultraconserved Elements Obtained from Pinned Insect Specimens.

Obtaining sequence data from historical museum specimens has been a growing research interest, invigorated by next-generation sequencing methods that allow inputs of highly degraded DNA. We applied a target enrichment and next-generation sequencing protocol to generate ultraconserved elements (UCEs) from 51 large carpenter bee specimens (genus Xylocopa), representing 25 species with specimen ages ranging from 2–121 years. We measured the correlation between specimen age and DNA yield (pre- and post-library preparation DNA concentration) and several UCE sequence capture statistics (raw read count, UCE reads on target, UCE mean contig length and UCE locus count) with linear regression models. We performed piecewise regression to test for specific breakpoints in the relationship of specimen age and DNA yield and sequence capture variables. Additionally, we compared UCE data from newer and older specimens of the same species and reconstructed their phylogeny in order to confirm the validity of our data. We recovered 6–972 UCE loci from samples with pre-library DNA concentrations ranging from 0.06–9.8 ng/μL. All investigated DNA yield and sequence capture variables were significantly but only moderately negatively correlated with specimen age. Specimens of age 20 years or less had significantly higher pre- and post-library concentrations, UCE contig lengths, and locus counts compared to specimens older than 20 years. We found breakpoints in our data indicating a decrease of the initial detrimental effect of specimen age on pre- and post-library DNA concentration and UCE contig length starting around 21–39 years after preservation. Our phylogenetic results confirmed the integrity of our data, giving preliminary insights into relationships within Xylocopa. We consider the effect of additional factors not measured in this study on our age-related sequence capture results, such as DNA fragmentation and preservation method, and discuss the promise of the UCE approach for large-scale projects in insect phylogenomics using museum specimens.

Analysis of a Rapid Evolutionary Radiation Using Ultraconserved Elements: Evidence for a Bias in Some Multispecies Coalescent Methods

Rapid evolutionary radiations are expected to require large amounts of sequence data to resolve. To resolve these types of relationships many systematists believe that it will be necessary to collect data by next-generation sequencing (NGS) and use multispecies coalescent ("species tree") methods. Ultraconserved element (UCE) sequence capture is becoming a popular method to leverage the high throughput of NGS to address problems in vertebrate phylogenetics. Here we examine the performance of UCE data for gallopheasants (true pheasants and allies), a clade that underwent a rapid radiation 10-15 Ma. Relationships among gallopheasant genera have been difficult to establish. We used this rapid radiation to assess the performance of species tree methods, using ∼600 kilobases of DNA sequence data from ∼1500 UCEs. We also integrated information from traditional markers (nuclear intron data from 15 loci and three mitochondrial gene regions). Species tree methods exhibited troubling behavior. Two methods [Maximum Pseudolikelihood for Estimating Species Trees (MP-EST) and Accurate Species TRee ALgorithm (ASTRAL)] appeared to perform optimally when the set of input gene trees was limited to the most variable UCEs, though ASTRAL appeared to be more robust than MP-EST to input trees generated using less variable UCEs. In contrast, the rooted triplet consensus method implemented in Triplec performed better when the largest set of input gene trees was used. We also found that all three species tree methods exhibited a surprising degree of dependence on the program used to estimate input gene trees, suggesting that the details of likelihood calculations (e.g., numerical optimization) are important for loci with limited phylogenetic information. As an alternative to summary species tree methods we explored the performance of SuperMatrix Rooted Triple - Maximum Likelihood (SMRT-ML), a concatenation method that is consistent even when gene trees exhibit topological differences due to the multispecies coalescent. We found that SMRT-ML performed well for UCE data. Our results suggest that UCE data have excellent prospects for the resolution of difficult evolutionary radiations, though specific attention may need to be given to the details of the methods used to estimate species trees.

Comparison of Target-Capture and Restriction-Site Associated DNA Sequencing for Phylogenomics: A Test in Cardinalid Tanagers (Aves, Genus: Piranga).

Restriction-site associated DNA sequencing (RAD-seq) and target capture of specific genomic regions, such as ultraconserved elements (UCEs), are emerging as two of the most popular methods for phylogenomics using reduced-representation genomic data sets. These two methods were designed to target different evolutionary timescales: RAD-seq was designed for population-genomic level questions and UCEs for deeper phylogenetics. The utility of both data sets to infer phylogenies across a variety of taxonomic levels has not been adequately compared within the same taxonomic system. Additionally, the effects of uninformative gene trees on species tree analyses (for target capture data) have not been explored. Here, we utilize RAD-seq and UCE data to infer a phylogeny of the bird genus Piranga The group has a range of divergence dates (0.5-6 myr), contains 11 recognized species, and lacks a resolved phylogeny. We compared two species tree methods for the RAD-seq data and six species tree methods for the UCE data. Additionally, in the UCE data, we analyzed a complete matrix as well as data sets with only highly informative loci. A complete matrix of 189 UCE loci with 10 or more parsimony informative (PI) sites, and an approximately 80% complete matrix of 1128 PI single-nucleotide polymorphisms (SNPs) (from RAD-seq) yield the same fully resolved phylogeny of Piranga We inferred non-monophyletic relationships of Piranga lutea individuals, with all other a priori species identified as monophyletic. Finally, we found that species tree analyses that included predominantly uninformative gene trees provided strong support for different topologies, with consistent phylogenetic results when limiting species tree analyses to highly informative loci or only using less informative loci with concatenation or methods meant for SNPs alone.

A Phylogenomic Perspective on the Radiation of Ray-Finned Fishes Based upon Targeted Sequencing of Ultraconserved Elements (UCEs)

Ray-finned fishes constitute the dominant radiation of vertebrates with over 32,000 species. Although molecular phylogenetics has begun to disentangle major evolutionary relationships within this vast section of the Tree of Life, there is no widely available approach for efficiently collecting phylogenomic data within fishes, leaving much of the enormous potential of massively parallel sequencing technologies for resolving major radiations in ray-finned fishes unrealized. Here, we provide a genomic perspective on longstanding questions regarding the diversification of major groups of ray-finned fishes through targeted enrichment of ultraconserved nuclear DNA elements (UCEs) and their flanking sequence. Our workflow efficiently and economically generates data sets that are orders of magnitude larger than those produced by traditional approaches and is well-suited to working with museum specimens. Analysis of the UCE data set recovers a well-supported phylogeny at both shallow and deep time-scales that supports a monophyletic relationship between Amia and Lepisosteus (Holostei) and reveals elopomorphs and then osteoglossomorphs to be the earliest diverging teleost lineages. Our approach additionally reveals that sequence capture of UCE regions and their flanking sequence offers enormous potential for resolving phylogenetic relationships within ray-finned fishes.