Evolution of rhodopsin in flatfishes (Pleuronectiformes) is associated with depth and migratory behavior

Visual signals are involved in many fitness-related tasks and are therefore essential for survival in many species. Aquatic organisms are ideal systems to study visual evolution, as the high diversity of spectral properties in aquatic environments generates great potential for adaptation to different light conditions. Flatfishes are an economically important group, with over 800 described species distributed globally, including halibut, flounder, sole, and turbot. The diversity of flatfish species and wide array of environments they occupy provides an excellent opportunity to understand how this variation translates to molecular adaptation of vision genes. Using models of molecular evolution, we investigated how the light environments inhabited by different flatfish lineages have shaped evolution in the rhodopsin gene, which is responsible for mediating dim-light visual transduction. We found strong evidence for positive selection in rhodopsin, and this was correlated with both migratory behavior and several fundamental aspects of habitat, including depth and freshwater/marine evolutionary transitions. We also identified several mutations that likely affect the wavelength of peak absorbance of rhodopsin, and outline how these shifts in absorbance correlate with the response to the light spectrum present in different habitats. This is the first study of rhodopsin evolution in flatfishes that considers their extensive diversity, and our results highlight how ecologically-driven molecular adaptation has occurred across this group in response to transitions to novel light environments.

The gustatory stalk of the Remo flounder exemplifies how complex evolutionary novelties may arise

The appearance of evolutionary novelties is a central issue in biology. Since Darwin's theory, difficulties in explaining how novel intricate body parts arose have often been used by creationists and other deniers to challenge evolution. Here, we describe the gustatory stalk of the Remo flounder (Oncopterus darwinii), an anatomically and functionally complex organ presumably used as a chemoreceptor probe to detect prey buried in the substrate. We demonstrate that the gustatory stalk is derived from the first dorsal-fin ray, which acquired remarkable modifications in its external morphology, integument, skeleton, muscles, and nerves. Such structural innovations are echoed in both functional and ecological specializations. We reveal that the gustatory stalk arose through the gradual accumulation of changes that evolved at different levels of the phylogenetic tree of ray-finned fishes. At least five preconditions arose in nodes preceding Oncopterus darwinii. This finding constitutes an interesting example of how evolution can deeply remodel body parts to perform entirely new functions. In this case, a trivial support structure primitively used for swimming became a sophisticated sensory tool to uncover hidden prey.

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.

Insights into the Mitochondrial Genetic Makeup and Miocene Colonization of Primitive Flatfishes (Pleuronectiformes: Psettodidae) in the East Atlantic and Indo-West Pacific Ocean

The mitogenomic evolution of the Psettodes flatfishes is still poorly known from their range distribution in eastern Atlantic and Indo-West Pacific Oceans. The study delves into the matrilineal evolutionary pathway of these primitive flatfishes, with a specific focus on the complete mitogenome of the Psettodes belcheri species, as determined through next-generation sequencing. The mitogenome in question spans a length of 16,747 base pairs and comprises a total of 37 genes, including 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region. Notably, the mitogenome of P. belcheri exhibits a bias towards AT base pairs, with a composition of 54.15%, mirroring a similar bias observed in its close relative, Psettodes erumei, which showcases percentages of 53.07% and 53.61%. Most of the protein-coding genes commence with an ATG initiation codon, except for Cytochrome c oxidase I (COI), which initiates with a GTG codon. Additionally, four protein-coding genes commence with a TAA termination codon, while seven others exhibit incomplete termination codons. Furthermore, two protein-coding genes, namely NAD1 and NAD6, terminate with AGG and TAG stop codons, respectively. In the mitogenome of P. belcheri, the majority of transfer RNAs demonstrate the classical cloverleaf secondary structures, except for tRNA-serine, which lacks a DHU stem. Comparative analysis of conserved blocks within the control regions of two Psettodidae species unveiled that the CSB-II block extended to a length of 51 base pairs, surpassing the other blocks and encompassing highly variable sites. A comprehensive phylogenetic analysis using mitochondrial genomes (13 concatenated PCGs) categorized various Pleuronectiformes species, highlighting the basal position of the Psettodidae family and showed monophyletic clustering of Psettodes species. The approximate divergence time (35-10 MYA) between P. belcheri and P. erumei was estimated, providing insights into their separation and colonization during the early Miocene. The TimeTree analysis also estimated the divergence of two suborders, Psettodoidei and Pleuronectoidei, during the late Paleocene to early Eocene (56.87 MYA). The distribution patterns of Psettodes flatfishes were influenced by ocean currents and environmental conditions, contributing to their ecological speciation. In the face of climate change and anthropogenic activities, the conservation implications of Psettodes flatfishes are emphasized, underscoring the need for regulated harvesting and adaptive management strategies to ensure their survival in changing marine ecosystems. Overall, this study contributes to understanding the evolutionary history, genetic diversity, and conservation needs of Psettodes flatfishes globally. However, the multifaceted exploration of mitogenome and larger-scale genomic data of Psettodes flatfish will provide invaluable insights into their genetic characterization, evolutionary history, environmental adaptation, and conservation in the eastern Atlantic and Indo-West Pacific Oceans.

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.

Peptide mass fingerprinting of preserved collagen in archaeological fish bones for the identification of flatfish in European waters

Bones of Pleuronectiformes (flatfish) are often not identified to species due to the lack of diagnostic features on bones that allow adequate distinction between taxa. This hinders in-depth understanding of archaeological fish assemblages and particularly flatfish fisheries throughout history. This is especially true for the North Sea region, where several commercially significant species have been exploited for centuries, yet their archaeological remains continue to be understudied. In this research, eight peptide biomarkers for 18 different species of Pleuronectiformes from European waters are described using MALDI-TOF MS and liquid chromatography tandem mass spectrometry data obtained from modern reference specimens. Bone samples (n = 202) from three archaeological sites in the UK and France dating to the medieval period (ca seventh-sixteenth century CE) were analysed using zooarchaeology by mass spectrometry (ZooMS). Of the 201 that produced good quality spectra, 196 were identified as flatfish species, revealing a switch in targeted species through time and indicating that ZooMS offers a more reliable and informative approach for species identification than osteological methods alone. We recommend this approach for future studies of archaeological flatfish remains as the precise species uncovered from a site can tell much about the origin of the fish, where people fished and whether they traded between regions.

Comparative mitogenomics and phylogenetics of the family Carangidae with special emphasis on the mitogenome of the Indian Scad Decapterus russelli

Carangids are abundant and commercially important marine fish that contribute to a significant portion of the fisheries in many parts of the world. In the present study, we characterized the complete mitogenome of the Indian scad, Decapterus russelli and performed a comprehensive comparative mitogenomic analysis of the family Carangidae. The comparative mitogenomics provided valuable insights into the structure, variability, and features of the coding and non-coding regions that evolved across species over millions of years. The structural features of tRNAs revealed changes in the frequency of mismatched and wobble base pairs, which is reflected in the base composition of H and L strands. The highly conserved sequence motif of the mTERF binding site in carangids over the ~ 400 MYA of their divergence demonstrated the functional importance of these sites. The control region of carangids was characterized by the presence of discontinuous repeat units with a high rate of sequence divergence in the form of base substitutions, insertions, and deletions. The maintenance of secondary structures in the control region independent of the rapid evolution of primary structure suggested the effect of selective constraints on their maintenance. Maximum likelihood (ML) and Bayesian inference (BI) phylogeny revealed a similar topology consistent with previous taxonomic studies. The extant carangids diverged through the evolutionary events experienced during the Cretaceous, Paleogene, and Neogene periods.

Phylogenetics of archerfishes (Toxotidae) and evolution of the toxotid shooting apparatus

Archerfishes (Toxotidae) are variously found in the fresh- and brackish-water environments of Asia Pacific and are well known for their ability to shoot water at terrestrial prey. These shots of water are intended to strike their prey and cause it to fall into the water for capture and consumption. While this behavior is well known, there are competing hypotheses (blowpipe vs. pressure tank hypothesis) of how archerfishes shoot and which oral structures are involved. Current understanding of archerfish shooting structures is largely based on two species, Toxotes chatareus and T. jaculatrix. We do not know if all archerfishes possess the same oral structures to shoot water, if anatomical variation is present within these oral structures, or how these features have evolved. Additionally, there is little information on the evolution of the Toxotidae as a whole, with all previous systematic works focusing on the interrelationships of the family. We first investigate the limits of archerfish species using new and previously published genetic data. Our analyses highlight that the current taxonomy of archerfishes does not conform to the relationships we recover. Toxotes mekongensis and T. siamensis are placed in synonymy of T. chatareus, Toxotes carpentariensis is recognized as a species and removed from synonymy of T. chatareus, and the genus Protoxotes is recognized for T. lorentzi based on the results of our analyses. We then take an integrative approach, using a combined analysis of discrete hard- and soft-tissue morphological characters with genetic data, to construct a phylogeny of the Toxotidae. Using the resulting phylogenetic hypothesis, we then characterize the evolutionary history and anatomical variation within the archerfishes. We discuss the variation in the oral structures and the evolution of the mechanism with respect to the interrelationships of archerfishes, and find that the oral structures of archerfishes support the blowpipe hypothesis but soft-tissue oral structures may also play a role in shooting. Finally, by comparing the morphology of archerfishes to their sister group, we find that the Leptobramidae has relevant shooting features in the oral cavity, suggesting that some components of the archerfish shooting mechanism are examples of co-opted or exapted traits.

Molecular phylogeny of the threadfin fishes (Polynemidae) using ultraconserved elements

Threadfins (Teleostei: Polynemidae) are a group of fishes named for their elongated and threadlike pectoral-fin rays. These fishes are commonly found in the world's tropical and subtropical waters, and are an economically important group for people living in these regions, with more than 100,000 t harvested in recent years. However, we do not have a detailed understanding of polynemid evolutionary history such that these fishes can be monitored, managed and conserved as an important tropical food source. Recent studies hypothesize at least one genus of threadfins is polyphyletic, and no studies have focused on generating a hypothesis of relationship for the Polynemidae using DNA sequences. In this study, we analyse a genomic dataset of ultraconserved-element and mitochondrial loci to construct a phylogeny of the Polynemidae. We recover the threadfins as a clade sister to flatfishes, with the most taxonomically rich genus, Polydactylus, being resolved as polyphyletic. When comparing our dataset to data from previous studies, we find that a few recent broad-scale phylogenies of fishes have incorporated mislabelled, misidentified or chimeric terminals into their analyses, impacting the relationships of threadfins they recover. We highlight these problematic sequences, providing revised identifications based on the data sequenced in this study. We then discuss the intrarelationships of threadfins, highlighting morphological or ecological characters that support the clades we recover.

Investigations into the ancestry of the Grape-eye Seabass (Hemilutjanus macrophthalmos) reveal novel limits and relationships for the Acropomatiformes (Teleostei: Percomorpha)

Abstract For 175 years, an unremarkable bass, the Grape-eye Seabass (Hemilutjanus macrophthalmos), has been known from coastal waters in the Eastern Pacific. To date, its phylogenetic placement and classification have been ignored. A preliminary osteological examination of Hemilutjanus hinted that it may have affinities with the Acropomatiformes. To test this hypothesis, we conducted a phylogenetic analysis using UCE and Sanger sequence data to study the placement of Hemilutjanus and the limits and relationships of the Acropomatiformes. We show that Hemilutjanus is a malakichthyid, and our results corroborate earlier studies that have resolved a polyphyletic Polyprionidae; accordingly, we describe Stereolepididae, new family, for Stereolepis. With these revisions, the Acropomatiformes is now composed of the: Acropomatidae; Banjosidae; Bathyclupeidae; Champsodontidae; Creediidae; Dinolestidae; Epigonidae; Glaucosomatidae; Hemerocoetidae; Howellidae; Lateolabracidae; Malakichthyidae; Ostracoberycidae; Pempheridae; Pentacerotidae; Polyprionidae; Scombropidae; Stereolepididae, new family; Symphysanodontidae; Synagropidae; and Schuettea. Finally, using our new hypothesis, we demonstrate that acropomatiforms repeatedly evolved bioluminescence and transitioned between shallow waters and the deep sea.

Interrogating Phylogenetic Discordance Resolves Deep Splits in the Rapid Radiation of Old World Fruit Bats (Chiroptera: Pteropodidae)

The family Pteropodidae (Old World fruit bats) comprises $>$200 species distributed across the Old World tropics and subtropics. Most pteropodids feed on fruit, suggesting an early origin of frugivory, although several lineages have shifted to nectar-based diets. Pteropodids are of exceptional conservation concern with $>$50% of species considered threatened, yet the systematics of this group has long been debated, with uncertainty surrounding early splits attributed to an ancient rapid diversification. Resolving the relationships among the main pteropodid lineages is essential if we are to fully understand their evolutionary distinctiveness, and the extent to which these bats have transitioned to nectar-feeding. Here we generated orthologous sequences for $>$1400 nuclear protein-coding genes (2.8 million base pairs) across 114 species from 43 genera of Old World fruit bats (57% and 96% of extant species- and genus-level diversity, respectively), and combined phylogenomic inference with filtering by information content to resolve systematic relationships among the major lineages. Concatenation and coalescent-based methods recovered three distinct backbone topologies that were not able to be reconciled by filtering via phylogenetic information content. Concordance analysis and gene genealogy interrogation show that one topology is consistently the best supported, and that observed phylogenetic conflicts arise from both gene tree error and deep incomplete lineage sorting. In addition to resolving long-standing inconsistencies in the reported relationships among major lineages, we show that Old World fruit bats have likely undergone at least seven independent dietary transitions from frugivory to nectarivory. Finally, we use this phylogeny to identify and describe one new genus. [Chiroptera; coalescence; concordance; incomplete lineage sorting; nectar feeder; species tree; target enrichment.].

Exon-capture data and locus screening provide new insights into the phylogeny of flatfishes (Pleuronectoidei)

There is an extensive collection of literature on the taxonomy and phylogenetics of flatfishes (Pleuronectiformes) that extends over two centuries, but consensus on many of their evolutionary relationships remains elusive. Phylogenetic uncertainty stems from highly divergent results derived from morphological and genetic characters, and between various molecular datasets. Deciphering relationships is complicated by rapid diversification early in the Pleuronectiformes tree and an abundance of studies that incompletely and inconsistently sample taxa and genetic markers. We present phylogenies based on a genome-wide dataset (4,434 nuclear markers via exon-capture) and wide taxon sampling (86 species spanning 12 of 16 families) of the largest flatfish suborder (Pleuronectoidei). Nine different subsets of the data and two tree construction approaches (eighteen phylogenies in total) are remarkably consistent with other recent molecular phylogenies, and show strong support for the monophyly of all families included except Pleuronectidae. Analyses resolved a novel phylogenetic hypothesis for the family Rhombosoleidae as being within the Pleuronectoidea rather than the Soleoidea, and failed to support the subfamily Hippoglossinae as a monophyletic group. Our results were corroborated with evidence from previous phylogenetic studies to outline regions of persistent phylogenetic uncertainty and identify groups in need of further phylogenetic inference.

Phylogenetic analyses of ray-finned fishes (Actinopterygii) using collagen type I protein sequences

Ray-finned fishes (Actinopterygii) are the largest and most diverse group of vertebrates, comprising over half of all living vertebrate species. Phylogenetic relationships between ray-finned fishes have historically pivoted on the study of morphology, which has notoriously failed to resolve higher order relationships, such as within the percomorphs. More recently, comprehensive genomic analyses have provided further resolution of actinopterygian phylogeny, including higher order relationships. Such analyses are rightfully regarded as the 'gold standard' for phylogenetics. However, DNA retrieval requires modern or well-preserved tissue and is less likely to be preserved in archaeological or fossil specimens. By contrast, some proteins, such as collagen, are phylogenetically informative and can survive into deep time. Here, we test the utility of collagen type I amino acid sequences for phylogenetic estimation of ray-finned fishes. We estimate topology using Bayesian approaches and compare the congruence of our estimated trees with published genomic phylogenies. Furthermore, we apply a Bayesian molecular clock approach and compare estimated divergence dates with previously published genomic clock analyses. Our collagen-derived trees exhibit 77% of node positions as congruent with recent genomic-derived trees, with the majority of discrepancies occurring in higher order node positions, almost exclusively within the Percomorpha. Our molecular clock trees present divergence times that are fairly comparable with genomic-based phylogenetic analyses. We estimate the mean node age of Actinopteri at ∼293 million years (Ma), the base of Teleostei at ∼211 Ma and the radiation of percomorphs beginning at ∼141 Ma (∼350 Ma, ∼250-283 Ma and ∼120-133 Ma in genomic trees, respectively). Finally, we show that the average rate of collagen (I) sequence evolution is 0.9 amino acid substitutions for every million years of divergence, with the α3 (I) sequence evolving the fastest, followed by the α2 (I) chain. This is the quickest rate known for any vertebrate group. We demonstrate that phylogenetic analyses using collagen type I amino acid sequences generate tangible signals for actinopterygians that are highly congruent with recent genomic-level studies. However, there is limited congruence within percomorphs, perhaps due to clade-specific functional constraints acting upon collagen sequences. Our results provide important insights for future phylogenetic analyses incorporating extinct actinopterygian species via collagen (I) sequencing.

Large-scale sequencing of flatfish genomes provides insights into the polyphyletic origin of their specialized body plan

The evolutionary and genetic origins of the specialized body plan of flatfish are largely unclear. We analyzed the genomes of 11 flatfish species representing 9 of the 14 Pleuronectiforme families and conclude that Pleuronectoidei and Psettodoidei do not form a monophyletic group, suggesting independent origins from different percoid ancestors. Genomic and transcriptomic data indicate that genes related to WNT and retinoic acid pathways, hampered musculature and reduced lipids might have functioned in the evolution of the specialized body plan of Pleuronectoidei. Evolution of Psettodoidei involved similar but not identical genes. Our work provides valuable resources and insights for understanding the genetic origins of the unusual body plan of flatfishes.

New insights into the phylogeny and evolution of lady beetles (Coleoptera: Coccinellidae) by extensive sampling of genes and species

Ladybirds (family Coccinellidae) are one of the most diverse groups of beetles and globally comprise over 6000 species. Despite their scientific and economic significance, the taxonomy of Coccinellidae remains unstable, and we still know little about their evolutionary history. By using a small number of genes, previous phylogenetic analyses have not reliably resolved the relationships among major ladybird lineages. In this study, we sequenced 94 nuclear protein-coding genes for 214 species of Coccinellidae and 14 outgroups, covering 90 genera and 35 tribes. We found that nucleotide compositional heterogeneity is present among ladybird tribes so that phylogenetic inference at the amino acid level is more reliable than at the DNA level. Based on the maximum likelihood analyses of the amino acid dataset, we recognize three subfamilies in Coccinellidae: Microweiseinae, Monocoryninae stat. nov., and Coccinellinae. The subfamily relationships are strongly supported as (Microweiseinae, (Monocoryninae stat. nov., Coccinellinae)). The tribes of ladybirds are mostly monophyletic, except Ortaliini, Sticholotidini, Scymnini, and Coccidulini. The phylogenetic relationships among tribes of Coccinellinae are still not well resolved, with many nodes weakly supported. Our divergence time analysis suggests that the crown group of extant lady beetles arose in the Early Cretaceous ~ 143 million years ago (Mya) and experienced a rapid diversification during the Late Cretaceous (120-70 Mya). We hypothesize that the boom of angiosperms in the Late Cretaceous promoted the diversification of herbivorous Sternorrhyncha insects, especially aphids, which in turn drove the rapid radiation of predatory lady beetles. In summary, our work provides a comprehensive time-calibrated phylogeny of Coccinellidae that provides a sound framework for revising their classification and understanding the origin of their biodiversity.

Revisiting the phylogeny of the genus Lolliguncula Steenstrup 1881 improves understanding of their biogeography and proves the validity of Lolliguncula argus Brakoniecki & Roper, 1985

The biogeography of American loliginid squids has been improved in recent years, but certain key taxa have been missing. Given that the most accurate phylogenies and estimates of divergence times of common ancestors depend heavily on good taxonomic coverage we have reanalyzed the genus Lolliguncula in light of new samples that increase the geographic and taxonomic coverage. New sequences were produced using standard methods to update an existing dataset for COI, 16S and Rhodopsin markers. Data was analyzed using various species delimitation methods, rigorous phylogenetic analyses and estimates of divergence times between clades. Within Lolliguncula we recover five monophyletic lineages that relate to the known species L. argus, L. diomedeae, L. panamensis, L. brevis North Atlantic and L. brevis South Atlantic. Except when using low divergence thresholds in ABGD, species delimitation methods only identify four of these lineages as distinct species, grouping L. argus and L. diomedeae as a single species. However, considering the reciprocal monophyly, recent divergence time estimate and morphological diagnoses we refrain from synonymizing L. argus within L. diomedeae, considering them very recently diverged species. The biogeography of the American loliginids is discussed, wherein basal cladogenesis in both Lolliguncula and Doryteuthis occur between the Atlantic and Pacific about 45 mya, with subsequent speciation around 20 mya associated with seafloor changes during the formation of the Caribbean. The recent speciation between L. argus and L. diomedeae is associated to oceanic environmental changes associated with glaciation, deep sea cooling and tropical upwelling.

Flatfishes colonised freshwater environments by acquisition of various DHA biosynthetic pathways

The colonisation of freshwater environments by marine fishes has historically been considered a result of adaptation to low osmolality. However, most marine fishes cannot synthesise the physiologically indispensable fatty acid, docosahexaenoic acid (DHA), due to incomplete DHA biosynthetic pathways, which must be adapted to survive in freshwater environments where DHA is poor relative to marine environments. By analysing DHA biosynthetic pathways of one marine and three freshwater-dependent species from the flatfish family Achiridae, we revealed that functions of fatty acid metabolising enzymes have uniquely and independently evolved by multi-functionalisation or neofunctionalisation in each freshwater species, such that every functional combination of the enzymes has converged to generate complete and functional DHA biosynthetic pathways. Our results demonstrate the elaborate patchwork of fatty acid metabolism and the importance of acquiring DHA biosynthetic function in order for fish to cross the nutritional barrier at the mouth of rivers and colonise freshwater environments.

Matsushita et al. demonstrate the evolution of DHA biosynthetic mechanisms in four species of flatfish as some of them colonised freshwater environments. Their analyses show independent changes to the biosynthetic pathways as a way to overcome the lack of exogenous DHA that would typically be available from prey in the marine environment.

imPhy: Imputing Phylogenetic Trees with Missing Information Using Mathematical Programming

Advances in modern genomics have allowed researchers to apply phylogenetic analyses on a genome-wide scale. While large volumes of genomic data can be generated cheaply and quickly, data missingness is a non-trivial and somewhat expected problem. Since the available information is often incomplete for a given set of genetic loci and individual organisms, a large proportion of trees that depict the evolutionary history of a single genetic locus, called gene trees, fail to contain all individuals. Data incompleteness causes difficulties in data collection, information extraction, and gene tree inference. Furthermore, identifying outlying gene trees, which can represent horizontal gene transfers, gene duplications, or hybridizations, is difficult when data is missing from the gene trees. The typical approach is to remove all individuals with missing data from the gene trees, and focus the analysis on individuals whose information is fully available - a huge loss of information. In this work, we propose and design an optimization-based imputation approach to infer the missing distances between leaves in a set of gene trees via a mixed integer non-linear programming model. We also present a new research pipeline, imPhy, that can (i) simulate a set of gene trees with leaves randomly missing in each tree, (ii) impute the missing pairwise distances in each gene tree, (iii) reconstruct the gene trees using the Neighbor Joining (NJ) and Unweighted Pair Group Method with Arithmetic Mean (UPGMA) methods, and (iv) analyze and report the efficiency of the reconstruction. To impute the missing leaves, we employ our newly proposed non-linear programming framework, and demonstrate its capability in reconstructing gene trees with incomplete information in both simulated and empirical datasets. In the empirical datasets apicomplexa and lungfish, our imputation has very small normalized mean square errors, even in the extreme case where 50 percent of the individuals in each gene tree are missing. Data, software, and user manuals can be found at https://github.com/yasuiniko/imPhy.

Comprehensive Transcriptome Analysis Reveals Insights into Phylogeny and Positively Selected Genes of Sillago Species

Sillago species lives in the demersal environments and face multiple stressors, such as localized oxygen depletion, sulfide accumulation, and high turbidity. In this study, we performed transcriptome analyses of seven Sillago species to provide insights into the phylogeny and positively selected genes of this species. After de novo assembly, 82,024, 58,102, 63,807, 85,990, 102,185, 69,748, and 102,903 unigenes were generated from S. japonica, S. aeolus, S. sp.1, S. sihama, S. sp.2, S. parvisquamis, and S. sinica, respectively. Furthermore, 140 shared orthologous exon markers were identified and then applied to reconstruct the phylogenetic relationships of the seven Sillago species. The reconstructed phylogenetic structure was significantly congruent with the prevailing morphological and molecular biological view of Sillago species relationships. In addition, a total of 44 genes were identified to be positively selected, and these genes were potential participants in the stress response, material (carbohydrate, amino acid and lipid) and energy metabolism, growth and differentiation, embryogenesis, visual sense, and other biological processes. We suspected that these genes possibly allowed Sillago species to increase their ecological adaptation to multiple environmental stressors.

A Novel Test for Absolute Fit of Evolutionary Models Provides a Means to Correctly Identify the Substitution Model and the Model Tree

A novel test is described that visualizes the absolute model-data fit of the substitution and tree components of an evolutionary model. The test utilizes statistics based on counts of character state matches and mismatches in alignments of observed and simulated sequences. This comparison is used to assess model-data fit. In simulations conducted to evaluate the performance of the test, the test estimator was able to identify both the correct tree topology and substitution model under conditions where the Goldman-Cox test-which tests the fit of a substitution model to sequence data and is also based on comparing simulated replicates with observed data-showed high error rates. The novel test was found to identify the correct tree topology within a wide range of DNA substitution model misspecifications, indicating the high discriminatory power of the test. Use of this test provides a practical approach for assessing absolute model-data fit when testing phylogenetic hypotheses.

The Prevalence and Impact of Model Violations in Phylogenetic Analysis

In phylogenetic inference, we commonly use models of substitution which assume that sequence evolution is stationary, reversible, and homogeneous (SRH). Although the use of such models is often criticized, the extent of SRH violations and their effects on phylogenetic inference of tree topologies and edge lengths are not well understood. Here, we introduce and apply the maximal matched-pairs tests of homogeneity to assess the scale and impact of SRH model violations on 3,572 partitions from 35 published phylogenetic data sets. We show that roughly one-quarter of all the partitions we analyzed (23.5%) reject the SRH assumptions, and that for 25% of data sets, tree topologies inferred from all partitions differ significantly from topologies inferred using the subset of partitions that do not reject the SRH assumptions. This proportion increases when comparing trees inferred using the subset of partitions that rejects the SRH assumptions, to those inferred from partitions that do not reject the SRH assumptions. These results suggest that the extent and effects of model violation in phylogenetics may be substantial. They highlight the importance of testing for model violations and possibly excluding partitions that violate models prior to tree reconstruction. Our results also suggest that further effort in developing models that do not require SRH assumptions could lead to large improvements in the accuracy of phylogenomic inference. The scripts necessary to perform the analysis are available in https://github.com/roblanf/SRHtests, and the new tests we describe are available as a new option in IQ-TREE (http://www.iqtree.org).

Digging for the spiny rat and hutia phylogeny using a gene capture approach, with the description of a new mammal subfamily

Next generation sequencing (NGS) and genomic database mining allow biologists to gather and select large molecular datasets well suited to address phylogenomics and molecular evolution questions. Here we applied this approach to a mammal family, the Echimyidae, for which generic relationships have been difficult to recover and often referred to as a star phylogeny. These South-American spiny rats represent a family of caviomorph rodents exhibiting a striking diversity of species and life history traits. Using a NGS exon capture protocol, we isolated and sequenced ca. 500 nuclear DNA exons for 35 species belonging to all major echimyid and capromyid clades. Exons were carefully selected to encompass as much diversity as possible in terms of rate of evolution, heterogeneity in the distribution of site-variation and nucleotide composition. Supermatrix inferences and coalescence-based approaches were subsequently applied to infer this family's phylogeny. The inferred topologies were the same for both approaches, and support was maximal for each node, entirely resolving the ambiguous relationships of previous analyses. Fast-evolving nuclear exons tended to yield more reliable phylogenies, as slower-evolving sequences were not informative enough to disentangle the short branches of the Echimyidae radiation. Based on this resolved phylogeny and on molecular and morphological evidence, we confirm the rank of the Caribbean hutias - formerly placed in the Capromyidae family - as Capromyinae, a clade nested within Echimyidae. We also name and define Carterodontinae, a new subfamily of Echimyidae, comprising the extant monotypic genus Carterodon from Brazil, which is the closest living relative of West Indies Capromyinae.

Post-Cretaceous bursts of evolution along the benthic-pelagic axis in marine fishes

Ecological opportunity arising in the aftermath of mass extinction events is thought to be a powerful driver of evolutionary radiations. Here, we assessed how the wake of the Cretaceous-Palaeogene (K-Pg) mass extinction shaped diversification dynamics in a clade of mostly marine fishes (Carangaria), which comprises a disparate array of benthic and pelagic dwellers including some of the most astonishing fish forms (e.g. flatfishes, billfishes, remoras, archerfishes). Analyses of lineage diversification show time-heterogeneous rates of lineage diversification in carangarians, with highest rates reached during the Palaeocene. Likewise, a remarkable proportion of Carangaria's morphological variation originated early in the history of the group and in tandem with a marked incidence of habitat shifts. Taken together, these results suggest that all major lineages and body plans in Carangaria originated in an early burst shortly after the K-Pg mass extinction, which ultimately allowed the occupation of newly released niches along the benthic-pelagic habitat axis.

CURatio: Genome-wide phylogenomic analysis method using ratios of total branch lengths

Evolutionary hypotheses provide important underpinnings of biological and medical sciences, and comprehensive, genome-wide understanding of evolutionary relationships among organisms are needed to test and refine such hypotheses. Theory and empirical evidence clearly indicate that phylogenies (trees) of different genes (loci) should not display precisely matching topologies. The main reason for such phylogenetic incongruence is reticulated evolutionary history of most species due to meiotic sexual recombination in eukaryotes, or horizontal transfers of genetic material in prokaryotes. Nevertheless, many genes should display topologically related phylogenies, and should group into one or more (for genetic hybrids) clusters in poly-dimensional "tree space". Unusual evolutionary histories or effects of selection may result in "outlier" genes with phylogenies that fall outside the main distribution(s) of trees in tree space. We present a new phylogenomic method, CURatio, which uses ratios of total branch lengths in gene trees to help identify phylogenetic outliers in a given set of ortholog groups from multiple genomes. An advantage of CURatio over other methods is that genes absent from and/or duplicated in some genomes can be included in the analysis. We conducted a simulation study under the coalescent model, and showed that, given sufficient species depth and topological difference, these ratios are significantly higher for the "outlier" gene phylogenies. Also, we applied CURatio to a set of annotated genomes of the fungal family, Clavicipitaceae, and identified alkaloid biosynthesis genes as outliers, probably due to a history of duplication and loss. The source code is available at https://github.com/QiwenKang/CURatio, and the empirical data set on Clavicipitaceae and simulated data set are available at Mendeley https://data.mendeley.com/datasets/mrxts7wjrr/1.

Systematics of Clupeiformes and testing for ecological limits on species richness in a trans-marine/freshwater clade

ABSTRACT Clupeiformes (herring, sardines, shad, anchovies and allies) are a globally distributed clade with nearly 400 marine, freshwater, and diadromous species. Although best known as filter feeding fishes that form large schools, this group occupies a diverse array of trophic guilds and habitats. Theory suggests that species richness in clades is modulated by ecological limits, which results in diversity-dependent clade growth, a pattern that most clades exhibit. As a trans-marine/freshwater clade that has undergone repeated transitions between marine and freshwaters, Clupeiformes are an excellent system for investigating the interplay between ecological diversity and macroevolutionary dynamics. In this study we review the systematics of Clupeiformes and explore discordance in phylogenetic relationships and divergence times between mitochondrial and nuclear loci. We then use comparative methods to test whether ecological limits regulate diversity in Clupeiformes. We find discordance in phylogenetic relationships at various taxonomic scales, but also considerable agreement between genomes. Our results suggest that trans-marine/freshwater clades are able to circumvent ecological limits on clade growth at regional, but not on local scales. Our study demonstrates that phylogenies are a critical link between ecology and macroevolutionary dynamics, and suggests habitat transitions can play a key role in shaping diversity patterns, particularly in the neotropics.

How the Central American Seaway and an Ancient Northern Passage Affected Flatfish Diversification

While the natural history of flatfish has been debated for decades, the mode of diversification of this biologically and economically important group has never been elucidated. To address this question, we assembled the largest molecular data set to date, covering > 300 species (out of ca. 800 extant), from 13 of the 14 known families over nine genes, and employed relaxed molecular clocks to uncover their patterns of diversification. As the fossil record of flatfish is contentious, we used sister species distributed on both sides of the American continent to calibrate clock models based on the closure of the Central American Seaway (CAS), and on their current species range. We show that flatfish diversified in two bouts, as species that are today distributed around the equator diverged during the closure of CAS, whereas those with a northern range diverged after this, hereby suggesting the existence of a postCAS closure dispersal for these northern species, most likely along a trans-Arctic northern route, a hypothesis fully compatible with paleogeographic reconstructions.

Reanalyzing the Palaeoptera problem - The origin of insect flight remains obscure

The phylogenetic relationships of the winged insect lineages - mayflies (Ephemeroptera), damselflies and dragonflies (Odonata), and all other winged insects (Neoptera) - are still controversial with three hypotheses supported by different datasets: Palaeoptera, Metapterygota and Chiastomyaria. Here, we reanalyze available phylogenomic data with a focus on detecting confounding and alternative signal. In this context, we provide a framework to quantitatively evaluate and assess incongruent molecular phylogenetic signal inherent in phylogenomic datasets. Despite overall support for the Palaeoptera hypothesis, we also found considerable signal for Chiastomyaria, which is not easily detectable by standardized tree inference approaches. Analyses of the accumulation of signal across gene partitions showed that signal accumulates gradually. However, even in case signal only slightly supported one over the other hypothesis, topologies inferred from large datasets switch from statistically strongly supported Palaeoptera to strongly supported Chiastomyaria. From a morphological point of view, Palaeoptera currently appears to be the best-supported hypothesis; however, recent analyses were restricted to head characters. Phylogenetic approaches covering all organ systems including analyses of potential functional or developmental convergence are still pending so that the Palaeoptera problem has to be considered an open question in insect systematics.

Improving phylogenetic inference of core Chlorophyta using chloroplast sequences with strong phylogenetic signals and heterogeneous models

Phylogenetic relationships within the green algal phylum Chlorophyta have proven difficult to resolve. The core Chlorophyta include Chlorophyceae, Ulvophyceae, Trebouxiophyceae, Pedinophyceae and Chlorodendrophyceae, but the relationships among these classes remain unresolved and the monophyly of Ulvophyceae and Trebouxiophyceae are highly controversial. We analyzed a dataset of 101 green algal species and 73 protein-coding genes sampled from complete and partial chloroplast genomes, including six newly sequenced ulvophyte genomes (Blidingia minima NIES-1837, Ulothrix zonata, Halochlorococcum sp. NIES-1838, Scotinosphaera sp. NIES-154, Caulerpa brownii and Cephaleuros sp. HZ-2017). We applied the Tree Certainty (TC) score to quantify the level of incongruence between phylogenetic trees in chloroplast genomic datasets, and show that the conflicting phylogenetic trees of core Chlorophyta stem from the most GC-heterogeneous sites. With removing the most GC-heterogeneous sites, our chloroplast phylogenomic analyses using heterogeneous models consistently support monophyly of the Chlorophyceae and of the Trebouxiophyceae, but the Ulvophyceae was resolved as polyphyletic. Our analytical framework provides an efficient approach to reconstruct the optimal phylogenetic relationships by minimizing conflicting signals.

Flatfish monophyly refereed by the relationship of Psettodes in Carangimorphariae

Background

The monophyly of flatfishes has not been supported in many molecular phylogenetic studies. The monophyly of Pleuronectoidei, which comprises all but one family of flatfishes, is broadly supported. However, the Psettodoidei, comprising the single family Psettodidae, is often found to be most closely related to other carangimorphs based on substantial sequencing efforts and diversely analytical methods. In this study, we examined why this particular result is often obtained.

Results

The mitogenomes of five flatfishes were determined. Select mitogenomes of representative carangimorph species were further employed for phylogenetic and molecular clock analyses. Our phylogenetic results do not fully support Psettodes as a sister group to pleuronectoids or other carangimorphs. And results also supported the evidence of long-branch attraction between Psettodes and the adjacent clades. Two chronograms, derived from Bayesian relaxed-clock methods, suggest that over a short period in the early Paleocene, a series of important evolutionary events occurred in carangimorphs.

Conclusion

Based on insights provided by the molecular clock, we propose the following evolutionary explanation for the difficulty in determining the phylogenetic position of Psettodes: The initial diversification of Psettodes was very close in time to the initial diversification of carangimorphs, and the primary diversification time of pleuronectoids, the other suborder of flatfishes, occurred later than that of some percomorph taxa. Additionally, the clade of Psettodes is long and naked branch, which supports the uncertainty of its phylogenetic placement. Finally, we confirmed the monophyly of flatfishes, which was accepted by most ichthyologists.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4788-5) contains supplementary material, which is available to authorized users.

Comprehensive phylogeny of ray-finned fishes (Actinopterygii) based on transcriptomic and genomic data

Our understanding of phylogenetic relationships among bony fishes has been transformed by analysis of a small number of genes, but uncertainty remains around critical nodes. Genome-scale inferences so far have sampled a limited number of taxa and genes. Here we leveraged 144 genomes and 159 transcriptomes to investigate fish evolution with an unparalleled scale of data: >0.5 Mb from 1,105 orthologous exon sequences from 303 species, representing 66 out of 72 ray-finned fish orders. We apply phylogenetic tests designed to trace the effect of whole-genome duplication events on gene trees and find paralogy-free loci using a bioinformatics approach. Genome-wide data support the structure of the fish phylogeny, and hypothesis-testing procedures appropriate for phylogenomic datasets using explicit gene genealogy interrogation settle some long-standing uncertainties, such as the branching order at the base of the teleosts and among early euteleosts, and the sister lineage to the acanthomorph and percomorph radiations. Comprehensive fossil calibrations date the origin of all major fish lineages before the end of the Cretaceous.

Revised classification of the righteye flounders (Teleostei: Pleuronectidae) based on multilocus phylogeny with complete taxon sampling

Members of the family Pleuronectidae are common representatives of the marine benthic fauna inhabiting northern regions of the Atlantic and Pacific oceans. The most recent comprehensive classification of the family, based entirely on morphological synapomorphies, recognized five subfamilies, 23 genera, and 61 extant species. However, several subsequent molecular studies have shown that many synapomorphic characters discovered in the morphological study might represent homoplasies, thereby questioning the reliance on these characters with the warning that they may provide misleading information for testing other morphology-based evolutionary hypotheses. In the present study, we propose a comprehensive taxonomic reassessment of the family Pleuronectidae based on the molecular phylogeny reconstructed from four nuclear and three mitochondrial loci and represented by complete taxon sampling of all but one valid species currently assigned to this family. To check for robustness of the phylogenetic hypothesis, we analyzed the effect of base compositional heterogeneity on phylogenetic signal for each locus and compared six different gene partitioning schemes. The final dataset, comprising 14 partitions and 154 individuals, was used to reconstruct phylogenetic trees in RAxML, MrBayes and BEAST2. Alternative topologies for several questionable nodes were compared using Bayes factors. The topology with the highest marginal likelihood was selected as the final phylogenetic tree for inferring pleuronectid relationships and character evolution. Based on our results, we recognize the Pleuronectidae comprising five subfamilies, 24 genera and 59 species. Our new phylogeny comprises five major monophyletic groups within the family, which we define as the subfamilies within the family: Atheresthinae, Pleuronichthyinae, Microstominae, Hippoglossinae and Pleuronectinae. Taxonomic composition of most of these subfamilies is different from that proposed in previous classifications. We also re-assess hypotheses proposed in earlier studies regarding intra-relationships of species of each lineage. Results of the current study contribute to better understanding of the evolutionary relationships of pleuronectid flatfishes based on molecular evidence, and they also provide the framework towards future comprehensive morphological revision of constituent lineages within the family Pleuronectidae.

The impact of phenotypic and molecular data on the inference of Colletotrichum diversity associated with Musa

Developing a comprehensive and reliable taxonomy for the Colletotrichum gloeosporioides species complex will require adopting data standards on the basis of an understanding of how methodological choices impact morphological evaluations and phylogenetic inference. We explored the impact of methodological choices in a morphological and molecular evaluation of Colletotrichum species associated with banana in Brazil. The choice of alignment filtering algorithm has a significant impact on topological inference and the retention of phylogenetically informative sites. Similarly, the choice of phylogenetic marker affects the delimitation of species boundaries, particularly if low phylogenetic signal is confounded with strong discordance, and inference of the species tree from multiple-gene trees. According to both phylogenetic informativeness profiling and Bayesian concordance analyses, the most informative loci are DNA lyase (APN2), intergenic spacer (IGS) between DNA lyase and the mating-type locus MAT1-2-1 (APN2/MAT-IGS), calmodulin (CAL), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), glutamine synthetase (GS), β-tubulin (TUB2), and a new marker, the intergenic spacer between GAPDH and an hypothetical protein (GAP2-IGS). Cornmeal agar minimizes the variance in conidial dimensions compared with potato dextrose agar and synthetic nutrient-poor agar, such that species are more readily distinguishable based on phenotypic differences. We apply these insights to investigate the diversity of Colletotrichum species associated with banana anthracnose in Brazil and report C. musae, C. tropicale, C. theobromicola, and C. siamense in association with banana anthracnose. One lineage did not cluster with any previously described species and is described here as C. chrysophilum.

Resolving Relationships among the Megadiverse Butterflies and Moths with a Novel Pipeline for Anchored Phylogenomics

The advent of next-generation sequencing technology has allowed for thecollection of large portions of the genome for phylogenetic analysis. Hybrid enrichment and transcriptomics are two techniques that leverage next-generation sequencing and have shown much promise. However, methods for processing hybrid enrichment data are still limited. We developed a pipeline for anchored hybrid enrichment (AHE) read assembly, orthology determination, contamination screening, and data processing for sequences flanking the target "probe" region. We apply this approach to study the phylogeny of butterflies and moths (Lepidoptera), a megadiverse group of more than 157,000 described species with poorly understood deep-level phylogenetic relationships. We introduce a new, 855 locus AHE kit for Lepidoptera phylogenetics and compare resulting trees to those from transcriptomes. The enrichment kit was designed from existing genomes, transcriptomes, and expressed sequence tags and was used to capture sequence data from 54 species from 23 lepidopteran families. Phylogenies estimated from AHE data were largely congruent with trees generated from transcriptomes, with strong support for relationships at all but the deepest taxonomic levels. We combine AHE and transcriptomic data to generate a new Lepidoptera phylogeny, representing 76 exemplar species in 42 families. The tree provides robust support for many relationships, including those among the seven butterfly families. The addition of AHE data to an existing transcriptomic dataset lowers node support along the Lepidoptera backbone, but firmly places taxa with AHE data on the phylogeny. Combining taxa sequenced for AHE with existing transcriptomes and genomes resulted in a tree with strong support for (Calliduloidea $+$ Gelechioidea $+$ Thyridoidea) $+$ (Papilionoidea $+$ Pyraloidea $+$ Macroheterocera). To examine the efficacy of AHE at a shallow taxonomic level, phylogenetic analyses were also conducted on a sister group representing a more recent divergence, the Saturniidae and Sphingidae. These analyses utilized sequences from the probe region and data flanking it, nearly doubled the size of the dataset; resulting trees supported new phylogenetics relationships, especially within the Saturniidae and Sphingidae (e.g., Hemarina derived in the latter). We hope that our data processing pipeline, hybrid enrichment gene set, and approach of combining AHE data with transcriptomes will be useful for the broader systematics community.

Incomplete Lineage Sorting in Mammalian Phylogenomics

The impact of incomplete lineage sorting (ILS) on phylogenetic conflicts among genes, and the related issue of whether to account for ILS in species tree reconstruction, are matters of intense controversy. Here, focusing on full-genome data in placental mammals, we empirically test two assumptions underlying current usage of tree-building methods that account for ILS. We show that in this data set (i) distinct exons from a common gene do not share a common genealogy, and (ii) ILS is only a minor determinant of the existing phylogenetic conflict. These results shed new light on the relevance and conditions of applicability of ILS-aware methods in phylogenomic analyses of protein coding sequences.

Multi-gene phylogeny of jacks and pompanos (Carangidae), including placement of monotypic vadigo Campogramma glaycos

In this study, the phylogenetic trees of jacks and pompanos (Carangidae), an ecologically and morphologically diverse, globally distributed fish family, are inferred from a complete, concatenated data set of two mitochondrial (cytochrome c oxidase I, cytochrome b) loci and one nuclear (myosin heavy chain 6) locus. Maximum likelihood and Bayesian inferences are largely congruent and show a clear separation of Carangidae into the four subfamilies: Scomberoidinae, Trachinotinae, Naucratinae and Caranginae. The inclusion of the carangid sister lineages Coryphaenidae (dolphinfishes) and Rachycentridae (cobia), however, render Carangidae paraphyletic. The phylogenetic trees also show with high statistical support that the monotypic vadigo Campogramma glaycos is the sister to all other species within the Naucratinae.

Maximizing Power in Phylogenetics and Phylogenomics: A Perspective Illuminated by Fungal Big Data

Since its original inception over 150 years ago by Darwin, we have made tremendous progress toward the reconstruction of the Tree of Life. In particular, the transition from analyzing datasets comprised of small numbers of loci to those comprised of hundreds of loci, if not entire genomes, has aided in resolving some of the most vexing of evolutionary problems while giving us a new perspective on biodiversity. Correspondingly, phylogenetic trees have taken a central role in fields that span ecology, conservation, and medicine. However, the rise of big data has also presented phylogenomicists with a new set of challenges to experimental design, quantitative analyses, and computation. The sequencing of a number of very first genomes presented significant challenges to phylogenetic inference, leading fungal phylogenomicists to begin addressing pitfalls and postulating solutions to the issues that arise from genome-scale analyses relevant to any lineage across the Tree of Life. Here we highlight insights from fungal phylogenomics for topics including systematics and species delimitation, ecological and phenotypic diversification, and biogeography while providing an overview of progress made on the reconstruction of the fungal Tree of Life. Finally, we provide a review of considerations to phylogenomic experimental design for robust tree inference. We hope that this special issue of Advances in Genetics not only excites the continued progress of fungal evolutionary biology but also motivates the interdisciplinary development of new theory and methods designed to maximize the power of genomic scale data in phylogenetic analyses.

Phylogenetic classification of bony fishes

BACKGROUND

Fish classifications, as those of most other taxonomic groups, are being transformed drastically as new molecular phylogenies provide support for natural groups that were unanticipated by previous studies. A brief review of the main criteria used by ichthyologists to define their classifications during the last 50 years, however, reveals slow progress towards using an explicit phylogenetic framework. Instead, the trend has been to rely, in varying degrees, on deep-rooted anatomical concepts and authority, often mixing taxa with explicit phylogenetic support with arbitrary groupings. Two leading sources in ichthyology frequently used for fish classifications (JS Nelson's volumes of Fishes of the World and W. Eschmeyer's Catalog of Fishes) fail to adopt a global phylogenetic framework despite much recent progress made towards the resolution of the fish Tree of Life. The first explicit phylogenetic classification of bony fishes was published in 2013, based on a comprehensive molecular phylogeny ( www.deepfin.org ). We here update the first version of that classification by incorporating the most recent phylogenetic results.

RESULTS

The updated classification presented here is based on phylogenies inferred using molecular and genomic data for nearly 2000 fishes. A total of 72 orders (and 79 suborders) are recognized in this version, compared with 66 orders in version 1. The phylogeny resolves placement of 410 families, or ~80% of the total of 514 families of bony fishes currently recognized. The ordinal status of 30 percomorph families included in this study, however, remains uncertain (incertae sedis in the series Carangaria, Ovalentaria, or Eupercaria). Comments to support taxonomic decisions and comparisons with conflicting taxonomic groups proposed by others are presented. We also highlight cases were morphological support exist for the groups being classified.

CONCLUSIONS

This version of the phylogenetic classification of bony fishes is substantially improved, providing resolution for more taxa than previous versions, based on more densely sampled phylogenetic trees. The classification presented in this study represents, unlike any other, the most up-to-date hypothesis of the Tree of Life of fishes.

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

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