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

The durability of natural infection and vaccine-induced immunity against future infection by SARS-CoV-2.

The durability of vaccine-mediated immunity to SARS-CoV-2, the durations to breakthrough infection, and the optimal timings of booster vaccination are crucial knowledge for pandemic response. Here, we applied comparative evolutionary analyses to estimate the durability of immunity and the likelihood of breakthrough infections over time following vaccination by BNT162b2 (Pfizer-BioNTech), mRNA-1273 (Moderna), ChAdOx1 (Oxford-AstraZeneca), and Ad26.COV2.S (Johnson & Johnson/Janssen). We evaluated anti-Spike (S) immunoglobulin G (IgG) antibody levels elicited by each vaccine relative to natural infection. We estimated typical trajectories of waning and corresponding infection probabilities, providing the distribution of times to breakthrough infection for each vaccine under endemic conditions. Peak antibody levels elicited by messenger RNA (mRNA) vaccines mRNA-1273 and BNT1262b2 exceeded that of natural infection and are expected to typically yield more durable protection against breakthrough infections (median 29.6 mo; 5 to 95% quantiles 10.9 mo to 7.9 y) than natural infection (median 21.5 mo; 5 to 95% quantiles 3.5 mo to 7.1 y). Relative to mRNA-1273 and BNT1262b2, viral vector vaccines ChAdOx1 and Ad26.COV2.S exhibit similar peak anti-S IgG antibody responses to that from natural infection and are projected to yield lower, shorter-term protection against breakthrough infection (median 22.4 mo and 5 to 95% quantiles 4.3 mo to 7.2 y; and median 20.5 mo and 5 to 95% quantiles 2.6 mo to 7.0 y; respectively). These results leverage the tools from evolutionary biology to provide a quantitative basis for otherwise unknown parameters that are fundamental to public health policy decision-making.

Ancient fish lineages illuminate toll-like receptor diversification in early vertebrate evolution.

Since its initial discovery over 50 years ago, understanding the evolution of the vertebrate RAG- mediated adaptive immune response has been a major area of research focus for comparative geneticists. However, how the evolutionary novelty of an adaptive immune response impacted the diversity of receptors associated with the innate immune response has received considerably less attention until recently. Here, we investigate the diversification of vertebrate toll-like receptors (TLRs), one of the most ancient and well conserved innate immune receptor families found across the Tree of Life, integrating genomic data that represent all major vertebrate lineages with new transcriptomic data from Polypteriformes, the earliest diverging ray-finned fish lineage. Our analyses reveal TLR sequences that reflect the 6 major TLR subfamilies, TLR1, TLR3, TLR4, TLR5, TLR7, and TLR11, and also currently unnamed, yet phylogenetically distinct TLR clades. We additionally recover evidence for a pulse of gene gain coincident with the rise of the RAG-mediated adaptive immune response in jawed vertebrates, followed by a period of rapid gene loss during the Cretaceous. These gene losses are primarily concentrated in marine teleost fish and synchronous with the mid Cretaceous anoxic event, a period of rapid extinction for marine species. Finally, we reveal a mismatch between phylogenetic placement and gene nomenclature for up to 50% of TLRs found in clades such as ray-finned fishes, cyclostomes, amphibians, and elasmobranchs. Collectively, these results provide an unparalleled perspective of TLR diversity and offer a ready framework for testing gene annotations in non-model species.

A highly diverse set of novel immunoglobulin-like transcript (NILT) genes in zebrafish indicates a wide range of functions with complex relationships to mammalian receptors.

Multiple novel immunoglobulin-like transcripts (NILTs) have been identified from salmon, trout, and carp. NILTs typically encode activating or inhibitory transmembrane receptors with extracellular immunoglobulin (Ig) domains. Although predicted to provide immune recognition in ray-finned fish, we currently lack a definitive framework of NILT diversity, thereby limiting our predictions for their evolutionary origin and function. In order to better understand the diversity of NILTs and their possible roles in immune function, we identified five NILT loci in the Atlantic salmon (Salmo salar) genome, defined 86 NILT Ig domains within a 3-Mbp region of zebrafish (Danio rerio) chromosome 1, and described 41 NILT Ig domains as part of an alternative haplotype for this same genomic region. We then identified transcripts encoded by 43 different NILT genes which reflect an unprecedented diversity of Ig domain sequences and combinations for a family of non-recombining receptors within a single species. Zebrafish NILTs include a sole putative activating receptor but extensive inhibitory and secreted forms as well as membrane-bound forms with no known signaling motifs. These results reveal a higher level of genetic complexity, interindividual variation, and sequence diversity for NILTs than previously described, suggesting that this gene family likely plays multiple roles in host immunity.

Infection by SARS-CoV-2 with alternate frequencies of mRNA vaccine boosting.

One of the most consequential unknowns of the COVID-19 pandemic is the frequency at which vaccine boosting provides sufficient protection from infection. We quantified the statistical likelihood of breakthrough infections over time following different boosting schedules with messenger RNA (mRNA)-1273 (Moderna) and BNT162b2 (Pfizer-BioNTech). We integrated anti-Spike IgG antibody optical densities with profiles of the waning of antibodies and corresponding probabilities of infection associated with coronavirus endemic transmission. Projecting antibody levels over time given boosting every 6 months, 1, 1.5, 2, or 3 years yielded respective probabilities of fending off infection over a 6-year span of >93%, 75%, 55%, 40%, and 24% (mRNA-1273) and >89%, 69%, 49%, 36%, and 23% (BNT162b2). Delaying the administration of updated boosters has bleak repercussions. It increases the probability of individual infection by SARS-CoV-2, and correspondingly, ongoing disease spread, prevalence, morbidity, hospitalization, and mortality. Instituting regular, population-wide booster vaccination updated to predominant variants has the potential to substantially forestall-and with global, widespread uptake, eliminate-COVID-19.

Placing human gene families into their evolutionary context.

Following the draft sequence of the first human genome over 20 years ago, we have achieved unprecedented insights into the rules governing its evolution, often with direct translational relevance to specific diseases. However, staggering sequence complexity has also challenged the development of a more comprehensive understanding of human genome biology. In this context, interspecific genomic studies between humans and other animals have played a critical role in our efforts to decode human gene families. In this review, we focus on how the rapid surge of genome sequencing of both model and non-model organisms now provides a broader comparative framework poised to empower novel discoveries. We begin with a general overview of how comparative approaches are essential for understanding gene family evolution in the human genome, followed by a discussion of analyses of gene expression. We show how homology can provide insights into the genes and gene families associated with immune response, cancer biology, vision, chemosensation, and metabolism, by revealing similarity in processes among distant species. We then explain methodological tools that provide critical advances and show the limitations of common approaches. We conclude with a discussion of how these investigations position us to gain fundamental insights into the evolution of gene families among living organisms in general. We hope that our review catalyzes additional excitement and research on the emerging field of comparative genomics, while aiding the placement of the human genome into its existentially evolutionary context.

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.

On the relationship between extant innate immune receptors and the evolutionary origins of jawed vertebrate adaptive immunity.

For over half a century, deciphering the origins of the genomic loci that form the jawed vertebrate adaptive immune response has been a major topic in comparative immunogenetics. Vertebrate adaptive immunity relies on an extensive and highly diverse repertoire of tandem arrays of variable (V), diversity (D), and joining (J) gene segments that recombine to produce different immunoglobulin (Ig) and T cell receptor (TCR) genes. The current consensus is that a recombination-activating gene (RAG)-like transposon invaded an exon of an ancient innate immune VJ-bearing receptor, giving rise to the extant diversity of Ig and TCR loci across jawed vertebrates. However, a model for the evolutionary relationships between extant non-recombining innate immune receptors and the V(D)J receptors of the jawed vertebrate adaptive immune system has only recently begun to come into focus. In this review, we provide an overview of non-recombining VJ genes, including CD8ß, CD79b, natural cytotoxicity receptor 3 (NCR3/NKp30), putative remnants of an antigen receptor precursor (PRARPs), and the multigene family of signal-regulatory proteins (SIRPs), that play a wide range of roles in immune function. We then focus in detail on the VJ-containing novel immune-type receptors (NITRs) from ray-finned fishes, as recent work has indicated that these genes are at least 50 million years older than originally thought. We conclude by providing a conceptual model of the evolutionary origins and phylogenetic distribution of known VJ-containing innate immune receptors, highlighting opportunities for future comparative research that are empowered by this emerging evolutionary perspective.

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

Transcriptome annotation reveals minimal immunogenetic diversity among Wyoming toads, Anaxyrus baxteri.

Briefly considered extinct in the wild, the future of the Wyoming toad (Anaxyrus baxteri) continues to rely on captive breeding to supplement the wild population. Given its small natural geographic range and history of rapid population decline at least partly due to fungal disease, investigation of the diversity of key receptor families involved in the host immune response represents an important conservation need. Population decline may have reduced immunogenetic diversity sufficiently to increase the vulnerability of the species to infectious diseases. Here we use comparative transcriptomics to examine the diversity of toll-like receptors and major histocompatibility complex (MHC) sequences across three individual Wyoming toads. We find reduced diversity at MHC genes compared to bufonid species with a similar history of bottleneck events. Our data provide a foundation for future studies that seek to evaluate the genetic diversity of Wyoming toads, identify biomarkers for infectious disease outcomes, and guide breeding strategies to increase genomic variability and wild release successes.

Molecular Biology and Evolution of Cancer: From Discovery to Action.

Cancer progression is an evolutionary process. During this process, evolving cancer cell populations encounter restrictive ecological niches within the body, such as the primary tumor, circulatory system, and diverse metastatic sites. Efforts to prevent or delay cancer evolution-and progression-require a deep understanding of the underlying molecular evolutionary processes. Herein we discuss a suite of concepts and tools from evolutionary and ecological theory that can inform cancer biology in new and meaningful ways. We also highlight current challenges to applying these concepts, and propose ways in which incorporating these concepts could identify new therapeutic modes and vulnerabilities in cancer.

Holosteans contextualize the role of the teleost genome duplication in promoting the rise of evolutionary novelties in the ray-finned fish innate immune system.

Over 99% of ray-finned fishes (Actinopterygii) are teleosts, a clade that comprises half of all living vertebrate species that have diversified across virtually all fresh and saltwater ecosystems. This ecological breadth raises the question of how the immunogenetic diversity required to persist under heterogeneous pathogen pressures evolved. The teleost genome duplication (TGD) has been hypothesized as the evolutionary event that provided the substrate for rapid genomic evolution and innovation. However, studies of putative teleost-specific innate immune receptors have been largely limited to comparisons either among teleosts or between teleosts and distantly related vertebrate clades such as tetrapods. Here we describe and characterize the receptor diversity of two clustered innate immune gene families in the teleost sister lineage: Holostei (bowfin and gars). Using genomic and transcriptomic data, we provide a detailed investigation of the phylogenetic history and conserved synteny of gene clusters encoding diverse immunoglobulin domain-containing proteins (DICPs) and novel immune-type receptors (NITRs). These data demonstrate an ancient linkage of DICPs to the major histocompatibility complex (MHC) and reveal an evolutionary origin of NITR variable-joining (VJ) exons that predate the TGD by at least 50 million years. Further characterizing the receptor diversity of Holostean DICPs and NITRs illuminates a sequence diversity that rivals the diversity of these innate immune receptor families in many teleosts. Taken together, our findings provide important historical context for the evolution of these gene families that challenge prevailing expectations concerning the consequences of the TGD during actinopterygiian evolution.

A comprehensive phylogeny of birds (Aves) using targeted next-generation DNA sequencing.

Although reconstruction of the phylogeny of living birds has progressed tremendously in the last decade, the evolutionary history of Neoaves--a clade that encompasses nearly all living bird species--remains the greatest unresolved challenge in dinosaur systematics. Here we investigate avian phylogeny with an unprecedented scale of data: >390,000 bases of genomic sequence data from each of 198 species of living birds, representing all major avian lineages, and two crocodilian outgroups. Sequence data were collected using anchored hybrid enrichment, yielding 259 nuclear loci with an average length of 1,523 bases for a total data set of over 7.8 × 10(7) bases. Bayesian and maximum likelihood analyses yielded highly supported and nearly identical phylogenetic trees for all major avian lineages. Five major clades form successive sister groups to the rest of Neoaves: (1) a clade including nightjars, other caprimulgiforms, swifts, and hummingbirds; (2) a clade uniting cuckoos, bustards, and turacos with pigeons, mesites, and sandgrouse; (3) cranes and their relatives; (4) a comprehensive waterbird clade, including all diving, wading, and shorebirds; and (5) a comprehensive landbird clade with the enigmatic hoatzin (Opisthocomus hoazin) as the sister group to the rest. Neither of the two main, recently proposed Neoavian clades--Columbea and Passerea--were supported as monophyletic. The results of our divergence time analyses are congruent with the palaeontological record, supporting a major radiation of crown birds in the wake of the Cretaceous-Palaeogene (K-Pg) mass extinction.

Transcriptome Ortholog Alignment Sequence Tools (TOAST) for phylogenomic dataset assembly.

BACKGROUND: Advances in next-generation sequencing technologies have reduced the cost of whole transcriptome analyses, allowing characterization of non-model species at unprecedented levels. The rapid pace of transcriptomic sequencing has driven the public accumulation of a wealth of data for phylogenomic analyses, however lack of tools aimed towards phylogeneticists to efficiently identify orthologous sequences currently hinders effective harnessing of this resource. RESULTS: We introduce TOAST, an open source R software package that can utilize the ortholog searches based on the software Benchmarking Universal Single-Copy Orthologs (BUSCO) to assemble multiple sequence alignments of orthologous loci from transcriptomes for any group of organisms. By streamlining search, query, and alignment, TOAST automates the generation of locus and concatenated alignments, and also presents a series of outputs from which users can not only explore missing data patterns across their alignments, but also reassemble alignments based on user-defined acceptable missing data levels for a given research question. CONCLUSIONS: TOAST provides a comprehensive set of tools for assembly of sequence alignments of orthologs for comparative transcriptomic and phylogenomic studies. This software empowers easy assembly of public and novel sequences for any target database of candidate orthologs, and fills a critically needed niche for tools that enable quantification and testing of the impact of missing data. As open-source software, TOAST is fully customizable for integration into existing or novel custom informatic pipelines for phylogenomic inference. Software, a detailed manual, and example data files are available through github carolinafishes.github.io.

Optimal Rates for Phylogenetic Inference and Experimental Design in the Era of Genome-Scale Data Sets.

With the rise of genome-scale data sets, there has been a call for increased data scrutiny and careful selection of loci that are appropriate to use in an attempt to resolve a phylogenetic problem. Such loci should maximize phylogenetic information content while minimizing the risk of homoplasy. Theory posits the existence of characters that evolve at an optimum rate, and efforts to determine optimal rates of inference have been a cornerstone of phylogenetic experimental design for over two decades. However, both theoretical and empirical investigations of optimal rates have varied dramatically in their conclusions: spanning no relationship to a tight relationship between the rate of change and phylogenetic utility. Herein, we synthesize these apparently contradictory views, demonstrating both empirical and theoretical conditions under which each is correct. We find that optimal rates of characters-not genes-are generally robust to most experimental design decisions. Moreover, consideration of site rate heterogeneity within a given locus is critical to accurate predictions of utility. Factors such as taxon sampling or the targeted number of characters providing support for a topology are additionally critical to the predictions of phylogenetic utility based on the rate of character change. Further, optimality of rates and predictions of phylogenetic utility are not equivalent, demonstrating the need for further development of comprehensive theory of phylogenetic experimental design. [Divergence time; GC bias; homoplasy; incongruence; information content; internode length; optimal rates; phylogenetic informativeness; phylogenetic theory; phylogenetic utility; phylogenomics; signal and noise; subtending branch length; state space; taxon and character sampling.].

Assessing phylogenetic information to reveal uncertainty in historical data: An example using Goodeinae (Teleostei: Cyprinodontiformes: Goodeidae).

A major emerging challenge to resolution of a stable phylogenetic Tree of Life has been incongruent inference among studies. Given the increasing ubiquity of incongruent studies, analyzing the predicted phylogenetic utility and quantitative evidence regarding contributions toward resolution of commonly-used markers in historical studies over the last decade represents an important, yet neglected, component of phylogenetics. Here we examine the phylogenetic utility of two sets of commonly-used legacy markers for understanding the evolutionary relationships among goodeines, a group of viviparous freshwater fishes endemic to central Mexico. Our analyses reveal that the validity of existing inferences is compromised by both lack of information and substantially biased patterns of nucleotide substitution. Our analyses demonstrate that many of the evolutionary relationships of goodeines remain uncertain - despite over a century of work. Our results provide an updated baseline of critically needed areas of investigation for the group and underscore the importance of quantifying phylogenetic information content as a fundamental step towards eroding false confidence in results based on weak and biased evidence.

Implications of lemuriform extinctions for the Malagasy flora.

Madagascar's lemurs display a diverse array of feeding strategies with complex relationships to seed dispersal mechanisms in Malagasy plants. Although these relationships have been explored previously on a case-by-case basis, we present here the first comprehensive analysis of lemuriform feeding, to our knowledge, and its hypothesized effects on seed dispersal and the long-term survival of Malagasy plant lineages. We used a molecular phylogenetic framework to examine the mode and tempo of diet evolution, and to quantify the associated morphological space occupied by Madagascar's lemurs, both extinct and extant. Using statistical models and morphometric analyses, we demonstrate that the extinction of large-bodied lemurs resulted in a significant reduction in functional morphological space associated with seed dispersal ability. These reductions carry potentially far-reaching consequences for Malagasy ecosystems, and we highlight large-seeded Malagasy plants that appear to be without extant animal dispersers. We also identify living lemurs that are endangered yet occupy unique and essential dispersal niches defined by our morphometric analyses.

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.

Eyes Wide Shut: the impact of dim-light vision on neural investment in marine teleosts.

Understanding how organismal design evolves in response to environmental challenges is a central goal of evolutionary biology. In particular, assessing the extent to which environmental requirements drive general design features among distantly related groups is a major research question. The visual system is a critical sensory apparatus that evolves in response to changing light regimes. In vertebrates, the optic tectum is the primary visual processing centre of the brain and yet it is unclear how or whether this structure evolves while lineages adapt to changes in photic environment. On one hand, dim-light adaptation is associated with larger eyes and enhanced light-gathering power that could require larger information processing capacity. On the other hand, dim-light vision may evolve to maximize light sensitivity at the cost of acuity and colour sensitivity, which could require less processing power. Here, we use X-ray microtomography and phylogenetic comparative methods to examine the relationships between diel activity pattern, optic morphology, trophic guild and investment in the optic tectum across the largest radiation of vertebrates-teleost fishes. We find that despite driving the evolution of larger eyes, enhancement of the capacity for dim-light vision generally is accompanied by a decrease in investment in the optic tectum. These findings underscore the importance of considering diel activity patterns in comparative studies and demonstrate how vision plays a role in brain evolution, illuminating common design principles of the vertebrate visual system.

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