Considering decoupled phenotypic diversification between ontogenetic phases in macroevolution: An example using Triggerfishes (Balistidae).

Across the Tree of Life, most studies of phenotypic disparity and diversification have been restricted to adult organisms. However, many lineages have distinct ontogenetic phases that differ from their adult forms in morphology and ecology. Focusing disproportionately on the evolution of adult forms unnecessarily hinders our understanding of the pressures shaping evolution over time. Non-adult disparity patterns are particularly important to consider for coastal ray-finned fishes, which often have juvenile phases with distinct phenotypes. These juvenile forms are often associated with sheltered nursery environments, with phenotypic shifts between adults and juvenile stages that are readily apparent in locomotor morphology. Whether this ontogenetic variation in locomotor morphology reflects a decoupling of diversification dynamics between life stages remains unknown. Here we investigate the evolutionary dynamics of locomotor morphology between adult and juvenile triggerfishes. We integrate a time-calibrated phylogenetic framework with geometric morphometric approaches and measurement data of fin aspect ratio and incidence, and reveal a mismatch between morphospace occupancy, the evolution of morphological disparity, and the tempo of trait evolution between life stages. Collectively, our results illuminate how the heterogeneity of morpho-functional adaptations can decouple the mode and tempo of morphological diversification between ontogenetic stages.

Phylogenomics of the Ancient and Species-Depauperate Gars Tracks 150 Million Years of Continental Fragmentation in the Northern Hemisphere.

Reconstructing deep-time biogeographic histories is limited by the comparatively recent diversification of most extant lineages. Ray-finned fishes, which include nearly half of all living vertebrates, are no exception. Although most lineages of ray-finned fishes radiated around the Cretaceous-Paleogene boundary, a handful of ancient, species-poor clades still persist. These lineages can illuminate very old biogeographic trends, but their low species richness can also limit the reconstruction of these patterns. The seven extant species of gars distributed in freshwater habitats in North America and Cuba are an old clade with a fossil record spanning over 150 million years of Earth history. Using a genomic data set of DNA sequences of 1105 exons for the seven living species and an updated morphological matrix of all extant and extinct taxa, we infer the phylogenetic relationships of gars and test how divergence times and biogeographic reconstructions are influenced by sequential and joint estimation and the effect on these inferences when using different taxon sets based on fossil completeness. Our analyses consistently show that the two extant gar genera Atractosteus and Lepisosteus diverged approximately 105 million years ago and many of the inferred divergences in the gar time-calibrated phylogeny closely track major Mesozoic tectonic events, including the separation of the Americas, the expansion of the early Atlantic, and the Cretaceous reorganization of North American river systems. The crown clades Atractosteus and Lepisosteus originated in the Cenozoic of eastern North America, implying that this region has served as both the origin of extant gar diversity and the refugium of this iconic ancient lineage. These results exemplify how combining phylogenomics with the fossil record provides congruence around the evolutionary history of ancient clades like gars and can reveal long-lost biogeographic patterns. [Lepisosteidae; fossilized birth death; biogeography; phylogenomics; paleontology.].

An inverse latitudinal gradient in speciation rate for marine fishes.

Far more species of organisms are found in the tropics than in temperate and polar regions, but the evolutionary and ecological causes of this pattern remain controversial1,2. Tropical marine fish communities are much more diverse than cold-water fish communities found at higher latitudes3,4, and several explanations for this latitudinal diversity gradient propose that warm reef environments serve as evolutionary 'hotspots' for species formation5-8. Here we test the relationship between latitude, species richness and speciation rate across marine fishes. We assembled a time-calibrated phylogeny of all ray-finned fishes (31,526 tips, of which 11,638 had genetic data) and used this framework to describe the spatial dynamics of speciation in the marine realm. We show that the fastest rates of speciation occur in species-poor regions outside the tropics, and that high-latitude fish lineages form new species at much faster rates than their tropical counterparts. High rates of speciation occur in geographical regions that are characterized by low surface temperatures and high endemism. Our results reject a broad class of mechanisms under which the tropics serve as an evolutionary cradle for marine fish diversity and raise new questions about why the coldest oceans on Earth are present-day hotspots of species formation.

Widespread sympatry in a species-rich clade of marine fishes (Carangoidei).

A universal paradigm describing patterns of speciation across the tree of life has been debated for decades. In marine organisms, inferring patterns of speciation using contemporary and historical patterns of biogeography is challenging due to the deficiency of species-level phylogenies and information on species' distributions, as well as conflicting relationships between species' dispersal, range size and co-occurrence. Most research on global patterns of marine fish speciation and biogeography has focused on coral reef or pelagic species. Carangoidei is an ecologically important clade of marine fishes that use coral reef and pelagic environments. We used sequence capture of 1314 ultraconserved elements (UCEs) from 154 taxa to generate a time-calibrated phylogeny of Carangoidei and its parent clade, Carangiformes. Age-range correlation analyses of the geographical distributions and divergence times of sister species pairs reveal widespread sympatry, with 73% of sister species pairs exhibiting sympatric geographical distributions, regardless of node age. Most species pairs coexist across large portions of their ranges. We also observe greater disparity in body length and maximum depth between sympatric relative to allopatric sister species. These and other ecological or behavioural attributes probably facilitate sympatry among the most closely related carangoids.

Phylogeography, hybridization, and species discovery in the Etheostoma nigrum complex (Percidae: Etheostoma: Boleosoma).

The history of riverine fish diversification is largely a product of geographic isolation. Physical barriers that reduce or eliminate gene flow between populations facilitate divergence via genetic drift and natural selection, eventually leading to speciation. For freshwater organisms, diversification is often the product of drainage basin rearrangements. In young clades where the history of isolation is the most recent, evolutionary relationships can resemble a tangled web. One especially recalcitrant group of freshwater fishes is the Johnny Darter (Etheostoma nigrum) species complex, where traditional taxonomy and molecular phylogenetics indicate a history of gene flow and conflicting inferences of species diversity. Here we assemble a genomic dataset using double digest restriction site associated DNA (ddRAD) sequencing and use phylogenomic and population genetic approaches to investigate the evolutionary history of the complex of species that includes E. nigrum, E. olmstedi, E. perlongum, and E. susanae. We reveal and validate several evolutionary lineages that we delimit as species, highlighting the need for additional work to formally describe the diversity of the Etheostoma nigrum complex. Our analyses also identify gene flow among recently diverged lineages, including one instance involving E. susanae, a localized and endangered species. Phylogeographic structure within the Etheostoma nigrum species complex coincides with major geologic events, such as parallel divergence in river basins during Pliocene inundation of the Atlantic coastal plain and multiple northward post-glacial colonization routes tracking river basin rearrangements. Our study serves as a nuanced example of how low dispersal rates coupled with geographic isolation among disconnected river systems in eastern North America has produced one of the world's freshwater biodiversity hotspots.

Diversification of the sleepers (Gobiiformes: Gobioidei: Eleotridae) and evolution of the root gobioid families.

Eleotridae (sleepers) and five smaller families are the earliest diverging lineages within Gobioidei. Most inhabit freshwaters in and around the Indo-Pacific, but Eleotridae also includes species that have invaded the Neotropics as well as several inland radiations in the freshwaters of Australia, New Zealand, and New Guinea. Previous efforts to infer phylogeny of these families have been based on sets of mitochondrial or nuclear loci and have yielded uncertain resolution of clades within Eleotridae. We expand the taxon sampling of previous studies and use genomic data from nuclear ultraconserved elements (UCEs) to infer phylogeny, then calibrate the hypothesis with recently discovered fossils. Our hypothesis clarifies ambiguously resolved relationships, provides a timescale for divergences, and indicates the core crown Eleotridae diverged over a short period 24.3-26.3 Ma in the late Oligocene. Within Eleotridae, we evaluate diversification dynamics with BAMM and find evidence for an overall slowdown in diversification over the past 35 Ma, but with a sharp increase 3.5 Ma in the genus Mogurnda, a clade of brightly colored species found in the freshwaters of Australia and New Guinea.

Introgression and Species Delimitation in the Longear Sunfish Lepomis megalotis (Teleostei: Percomorpha: Centrarchidae).

Introgression and hybridization are major impediments to genomic-based species delimitation because many implementations of the multispecies coalescent framework assume no gene flow among species. The sunfish genus Lepomis, one of the world's most popular groups of freshwater sport fish, has a complicated taxonomic history. The results of ddRAD phylogenomic analyses do not provide support for the current taxonomy that recognizes two species, Lepomis megalotis and Lepomis peltastes, in the L. megalotis complex. Instead, evidence from phylogenomics and phenotype warrants recognizing six relatively ancient evolutionary lineages in the complex. The introgressed and hybridizing populations in the L. megalotis complex are localized and appear to be the result of secondary contact or rare hybridization events between nonsister species. Segregating admixed populations from our multispecies coalescent analyses identifies six species with moderate to high genealogical divergence, whereas including admixed populations drives all but one lineage below the species threshold of genealogical divergence. Segregation of admixed individuals also helps reveal phenotypic distinctiveness among the six species in morphological traits used by ichthyologists to discover and delimit species over the last two centuries. Our protocols allow for the identification and accommodation of hybridization and introgression in species delimitation. Genomic-based species delimitation validated with multiple lines of evidence provides a path towards the discovery of new biodiversity and resolving long-standing taxonomic problems.[ddRAD; genealogical divergence index; hybridization; integrative species delimitation; phylogeny; secondary contact; systematics; taxonomy.].

Erratum: Early members of 'living fossil' lineage imply later origin of modern ray-finned fishes.

This corrects the article DOI: 10.1038/nature23654.

Early members of 'living fossil' lineage imply later origin of modern ray-finned fishes.

Modern ray-finned fishes (Actinopterygii) comprise half of extant vertebrate species and are widely thought to have originated before or near the end of the Middle Devonian epoch (around 385 million years ago). Polypterids (bichirs and ropefish) represent the earliest-diverging lineage of living actinopterygians, with almost all Palaeozoic taxa interpreted as more closely related to other extant actinopterygians than to polypterids. By contrast, the earliest material assigned to the polypterid lineage is mid-Cretaceous in age (around 100 million years old), implying a quarter-of-a-billion-year palaeontological gap. Here we show that scanilepiforms, a widely distributed radiation from the Triassic period (around 252-201 million years ago), are stem polypterids. Importantly, these fossils break the long polypterid branch and expose many supposedly primitive features of extant polypterids as reversals. This shifts numerous Palaeozoic ray-fins to the actinopterygian stem, reducing the minimum age for the crown lineage by roughly 45 million years. Recalibration of molecular clocks to exclude phylogenetically reassigned Palaeozoic taxa results in estimates that the actinopterygian crown lineage is about 20-40 million years younger than was indicated by previous molecular analyses. These new dates are broadly consistent with our revised palaeontological timescale and coincident with an interval of conspicuous morphological and taxonomic diversification among ray-fins centred on the Devonian-Carboniferous boundary. A shifting timescale, combined with ambiguity in the relationships of late Palaeozoic actinopterygians, highlights this part of the fossil record as a major frontier in understanding the evolutionary assembly of modern vertebrate diversity.

Periodic Environmental Disturbance Drives Repeated Ecomorphological Diversification in an Adaptive Radiation of Antarctic Fishes.

AbstractThe ecological theory of adaptive radiation has profoundly shaped our conceptualization of the rules that govern diversification. However, while many radiations follow classic early-burst patterns of diversification as they fill ecological space, the longer-term fates of these radiations depend on many factors, such as climatic stability. In systems with periodic disturbances, species-rich clades can contain nested adaptive radiations of subclades with their own distinct diversification histories, and how adaptive radiation theory applies in these cases is less clear. Here, we investigated patterns of ecological and phenotypic diversification within two iterative adaptive radiations of cryonotothenioid fishes in Antarctica's Southern Ocean: crocodile icefishes and notoperches. For both clades, we observe evidence of repeated diversification into disparate regions of trait space between closely related taxa and into overlapping regions of trait space between distantly related taxa. We additionally find little evidence that patterns of ecological divergence are correlated with evolution of morphological disparity, suggesting that these axes of divergence may not be tightly linked. Finally, we reveal evidence of repeated convergence in sympatry that suggests niche complementarity. These findings reflect the dynamic history of Antarctic marine habitats and may guide hypotheses of diversification dynamics in environments characterized by periodic disturbance.

Phylogenomic Species Delimitation Dramatically Reduces Species Diversity in an Antarctic Adaptive Radiation.

Application of genetic data to species delimitation often builds confidence in delimitations previously hypothesized using morphological, ecological, and geographic data and frequently yields recognition of previously undescribed cryptic diversity. However, a recent critique of genomic data-based species delimitation approaches is that they have the potential to conflate population structure with species diversity, resulting in taxonomic oversplitting. The need for an integrative approach to species delimitation, in which molecular, morphological, ecological, and geographic lines of evidence are evaluated together, is becoming increasingly apparent. Here, we integrate phylogenetic, population genetic, and coalescent analyses of genome-wide sequence data with investigation of variation in multiple morphological traits to delimit species within the Antarctic barbeled plunderfishes (Artedidraconidae: Pogonophryne). Pogonophryne currently comprises 29 valid species, most of which are distinguished solely by variation in the ornamentation of the mental barbel that projects from the lower jaw, a structure previously shown to vary widely within a single species. However, our genomic and phenotypic analyses result in a dramatic reduction in the number of distinct species recognized within the clade, providing evidence to support the recognition of no more than six species. We propose to synonymize 24 of the currently recognized species with five species of Pogonophryne. We find genomic and phenotypic evidence for a new species of Pogonophryne from specimens collected in the Ross Sea. Our findings represent a rare example in which the application of molecular data provides evidence of taxonomic oversplitting on the basis of morphology, clearly demonstrating the utility of an integrative species delimitation framework.[ddRADseq; multispecies coalescent; Notothenioidei; SNPs; Southern Ocean.].

Accelerated Diversification Explains the Exceptional Species Richness of Tropical Characoid Fishes.

The Neotropics harbor the most species-rich freshwater fish fauna on the planet, but the timing of that exceptional diversification remains unclear. Did the Neotropics accumulate species steadily throughout their long history, or attain their remarkable diversity recently? Biologists have long debated the relative support for these museum and cradle hypotheses, but few phylogenies of megadiverse tropical clades have included sufficient taxa to distinguish between them. We used 1288 ultraconserved element loci spanning 293 species, 211 genera, and 21 families of characoid fishes to reconstruct a new, fossil-calibrated phylogeny and infer the most likely diversification scenario for a clade that includes a third of Neotropical fish diversity. This phylogeny implies paraphyly of the traditional delimitation of Characiformes because it resolves the largely Neotropical Characoidei as the sister lineage of Siluriformes (catfishes), rather than the African Citharinodei. Time-calibrated phylogenies indicate an ancient origin of major characoid lineages and reveal a much more recent emergence of most characoid species. Diversification rate analyses infer increased speciation and decreased extinction rates during the Oligocene at around 30 Ma during a period of mega-wetland formation in the proto-Orinoco-Amazonas. Three species-rich and ecomorphologically diverse lineages (Anostomidae, Serrasalmidae, and Characidae) that originated more than 60 Ma in the Paleocene experienced particularly notable bursts of Oligocene diversification and now account collectively for 68% of the approximately 2150 species of Characoidei. In addition to paleogeographic changes, we discuss potential accelerants of diversification in these three lineages. While the Neotropics accumulated a museum of ecomorphologically diverse characoid lineages long ago, this geologically dynamic region also cradled a much more recent birth of remarkable species-level diversity. [Biodiversity; Characiformes; macroevolution; Neotropics; phylogenomics; ultraconserved elements.].

Hidden species diversity in an iconic living fossil vertebrate.

Ancient, species-poor lineages persistently occur across the Tree of life. These lineages are likely to contain unrecognized species diversity masked by the low rates of morphological evolution that characterize living fossils. Halecomorphi is a lineage of ray-finned fishes that diverged from its closest relatives before 200 Ma and is represented by only one living species in eastern North America, the bowfin, Amia calva Linnaeus. Here, we use double digest restriction-site-associated DNA sequencing and morphology to illuminate recent speciation in bowfins. Our results support the delimitation of a second living species of Amia, with the timing of diversification dating to the Plio-Pleistocene. This delimitation expands the species diversity of an ancient lineage that is integral to studies of vertebrate genomics and development, yet is facing growing conservation threats driven by the caviar fishery.

Genomic and phenotypic divergence informs translocation strategies for an endangered freshwater fish.

Translocation, the movement of organisms for conservation purposes, can result in unintended introgression if genetic material flows between populations in new ways. The Bluemask Darter Etheostoma akatulo is a federally endangered species of freshwater fish inhabiting the Caney Fork River system and three of its tributaries (Collins River, Rocky River, and Cane Creek) in Tennessee. The current conservation strategy for Bluemask Darters involves translocating the progeny of broodstock from the Collins River (in the west) to the Calfkiller River (in the east) where the species had been extirpated. In this study, we use ddRAD sequence data from across the extant range to assess this translocation strategy in light of population structure, phylogeny, and demography. We also include museum specimen data to assess morphological variation among extant and extirpated populations. Our analyses reveal substantial genetic and phenotypic disparities between a western population in the Collins River and an eastern population encompassing the Rocky River, Cane Creek, and upper Caney Fork, the two of which shared common ancestry more than 100,000 years ago. Furthermore, morphological analyses classify 12 of 13 Calfkiller River specimens with phenotypes consistent with the eastern population. These results suggest that current translocations perturb the evolutionary boundaries between two delimited populations. Instead, we suggest that repopulating the Calfkiller River using juveniles from the Rocky River could balance conflicting signatures of demography, diversity, and divergence. Beyond conservation, the microgeographic structure of Bluemask Darter populations adds another puzzle to the phylogeography of the hyperdiverse freshwater fishes in eastern North America.

Phylogeny and time scale of diversification in the fossil-rich sunfishes and black basses (Teleostei: Percomorpha: Centrarchidae).

Species of the North American freshwater fish lineage Centrarchidae are apex predators in their habitats and are among the world's most popular sport fishes. Centrarchids boast a rich fossil record that extends from the latest Eocene to the Pleistocene. To investigate the phylogeny and timing of diversification of Centrarchidae, we deploy a dataset of DNA sequences of 16 nuclear genes sampled from nearly all of the recognized and undescribed species. We also utilize previously published morphological datasets to assess the phylogenetic placement of one of the oldest known centrarchid fossils, †Plioplarchus whitei. A Bayesian multispecies coalescent species tree analysis provides insight on relationships that evaded resolution in earlier studies, such as the relationships of Acantharchus pomotis, the resolution of a clade consisting of species previously synonymized under the Spotted Bass, Micropterus punctulatus, and a clade of recently described species previously considered populations of the Redeye Bass, Micropterus coosae. This new molecular phylogeny and the inclusion of †P. whitei and other centrarchid fossils in the tip-dated fossilized birth-death analysis results in a new hypothesis of the timing of diversification in Centrarchidae that contextualizes the ages of centrarchid fossils to the timing of speciation among the extant species. In addition to providing new temporal perspectives on the diversification of freshwater fishes in North America, this study may close of the chapter of centrarchid phylogeny inferred using Sanger-sequenced genes, as the use of genomic-scale datasets becomes mainstream in the phylogenetics of fishes.

Phylogenomic Signatures of Ancient Introgression in a Rogue Lineage of Darters (Teleostei: Percidae).

Evolutionary history is typically portrayed as a branching phylogenetic tree, yet not all evolution proceeds in a purely bifurcating manner. Introgressive hybridization is one process that results in reticulate evolution. Most known examples of genome-wide introgression occur among closely related species with relatively recent common ancestry; however, we present evidence for ancient hybridization and genome-wide introgression between major stem lineages of darters, a species-rich clade of North American freshwater fishes. Previous attempts to resolve the relationships of darters have been confounded by the uncertain phylogenetic resolution of the lineage Allohistium. In this study, we investigate the phylogenomics of darters, specifically the relationships of Allohistium, through analyses of approximately 30,000 RADseq loci sampled from 112 species. Our phylogenetic inferences are based on traditional approaches in combination with strategies that accommodate reticulate evolution. These analyses result in a novel phylogenetic hypothesis for darters that includes ancient introgression between Allohistium and other two major darter lineages, minimally occurring 20 million years ago. Darters offer a compelling case for the necessity of incorporating phylogenetic networks in reconstructing the evolutionary history of diversification in species-rich lineages. We anticipate that the growing wealth of genomic data for clades of non-model organisms will reveal more examples of ancient hybridization, eventually requiring a re-evaluation of how evolutionary history is visualized and utilized in macroevolutonary investigations.

A phylogenomic framework for pelagiarian fishes (Acanthomorpha: Percomorpha) highlights mosaic radiation in the open ocean.

The fish clade Pelagiaria, which includes tunas as its most famous members, evolved remarkable morphological and ecological variety in a setting not generally considered conducive to diversification: the open ocean. Relationships within Pelagiaria have proven elusive due to short internodes subtending major lineages suggestive of rapid early divergences. Using a novel sequence dataset of over 1000 ultraconserved DNA elements (UCEs) for 94 of the 286 species of Pelagiaria (more than 70% of genera), we provide a time-calibrated phylogeny for this widely distributed clade. Some inferred relationships have clear precedents (e.g. the monophyly of 'core' Stromateoidei, and a clade comprising 'Gempylidae' and Trichiuridae), but others are unexpected despite strong support (e.g. Chiasmodontidae + Tetragonurus). Relaxed molecular clock analysis using node-based fossil calibrations estimates a latest Cretaceous origin for Pelagiaria, with crown-group families restricted to the Cenozoic. Estimated mean speciation rates decline from the origin of the group in the latest Cretaceous, although credible intervals for root and tip rates are broad and overlap in most cases, and there is higher-than-expected partitioning of body shape diversity (measured as fineness ratio) between clades concentrated during the Palaeocene-Eocene. By contrast, more direct measures of ecology show either no substantial deviation from a null model of diversification (diet) or patterns consistent with evolutionary constraint or high rates of recent change (depth habitat). Collectively, these results indicate a mosaic model of diversification. Pelagiarians show high morphological disparity and modest species richness compared to better-studied fish radiations in contrasting environments. However, this pattern is also apparent in other clades in open-ocean or deep-sea habitats, and suggests that comparative study of such groups might provide a more inclusive model of the evolution of diversity in fishes.

Phylogenetic analysis of Antarctic notothenioids illuminates the utility of RADseq for resolving Cenozoic adaptive radiations.

Notothenioids are a clade of ∼120 species of marine fishes distributed in extreme southern hemisphere temperate near-shore habitats and in the Southern Ocean surrounding Antarctica. Over the past 25 years, molecular and morphological approaches have redefined hypotheses of relationships among notothenioid lineages as well as their relationships among major lineages of percomorph teleosts. These phylogenies provide a basis for investigation of mechanisms of evolutionary diversification within the clade and have enhanced our understanding of the notothenioid adaptive radiation. Despite extensive efforts, there remain several questions concerning the phylogeny of notothenioids. In this study, we deploy DNA sequences of ∼100,000 loci obtained using RADseq to investigate the phylogenetic relationships of notothenioids and to assess the utility of RADseq loci for lineages that exhibit divergence times ranging from the Paleogene to the Quaternary. The notothenioid phylogenies inferred from the RADseq loci provide unparalleled resolution and node support for several long-standing problems including, (1) relationships among species of Trematomus, (2) resolution of Indonotothenia cyanobrancha as the sister lineage of Trematomus, (3) the deep paraphyly of Nototheniidae, (4) the paraphyly of Lepidonotothen s.l., (5) paraphyly of Artedidraco, and 6) the monophyly of the Bathydraconidae. Assessment of site rates demonstrates that RADseq loci are similar to mtDNA protein coding genes and exhibit peak phylogenetic informativeness at the time interval during which the major Antarctic notothenioid lineages originated and diversified. In addition to providing a well-resolved phylogenetic hypothesis for notothenioids, our analyses quantify the predicted utility of RADseq loci for Cenozoic phylogenetic inferences.

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.].