The phylogenetic position of toadfishes (order Batrachoidiformes) in the higher ray-finned fish as inferred from partitioned Bayesian analysis of 102 whole mitochondrial genome sequences

Mitogenomic Characterization and Phylogenetic Insights of the Ornamental Sail-Fin Molly (Poecilia velifera) in Non-Native Indonesian Waters

The ornamental fish Poecilia velifera (Sail-fin molly, Poeciliidae) has spread widely to various non-native ecosystems around the world, far from its native habitat in the Yucatan Peninsula, Mexico. Despite the availability of some partial mitochondrial and nuclear genetic information, the complete mitogenomic structure and its variation remain unknown for this species, which is essential for a comprehensive genetic characterization and detailed phylogenetic investigation. This study applied next-generation sequencing to generate the de novo mitogenome of morphologically identified P. velifera from a non-native brackish water ecosystem in Banten Province, Indonesia. The resulting mitogenome was 16,627 bp in length and encompassed 13 protein-coding genes (PCGs), 22 transfer RNAs, two ribosomal RNAs, and a non-coding control region (CR). The result enhances our understanding of the genetic makeup of P. velifera compared to its congeners. Furthermore, the identified nucleotide variations within the conserved blocks of the CR region could provide insights into the functional role of this non-coding region. Bayesian phylogenetic inference using concatenated PCGs distinguished P. velifera from its congeners and showed monophyletic clustering of Poecilia in the family Poeciliidae, consistent with earlier evolutionary hypotheses. This first mitogenome of P. velifera paves the way for using multiple mitochondrial markers in species identification and understanding population structure in the near future. In addition, looking into the genetic evidence of this ornamental species in a non-native ecosystem, the study emphasizes the importance of strict quarantine regulations to protect Indonesia's native fish species.

Elucidating the Mitogenomic Blueprint of Pomadasys perotaei from the Eastern Atlantic: Characterization and Matrilineal Phylogenetic Insights into Haemulid Grunts (Teleostei: Lutjaniformes)

The parrot grunt fish, Pomadasys perotaei, has a limited distribution in the Eastern Atlantic Ocean and is an important species in marine capture fisheries across several West African countries. Despite its ecological and economic significance, the mitogenomic information for this species is lacking. This study utilized next-generation sequencing to generate the de novo mitogenome of P. perotaei from Eastern Atlantic. The resulting mitogenome is 16,691 base pairs and includes 13 protein-coding genes (PCGs), 22 transfer RNAs, two ribosomal RNAs, and an AT-rich control region (CR). Most of the PCGs exhibit nonsynonymous (Ka) and synonymous (Ks) substitution rates of less than '1', indicating strong negative selection across haemulid fishes. The control region of Pomadasys species contains four conserved domains, as seen in other teleost's, with polymorphic nucleotides that can be used to study population structures through the amplification of short mitochondrial gene fragments. Additionally, Bayesian phylogenetic analysis based on PCGs revealed a non-monophyletic clustering pattern of Pomadasys within the haemulid matrilineal tree. Overall, the structural characterization and phylogenetic analysis enhance our understanding of the genetic composition and evolutionary history of Pomadasys species from the Indo-West Pacific and Eastern Atlantic Oceans.

Exploring the phylogeny and depth evolution of cusk eels and their relatives (Ophidiiformes: Ophidioidei)

With 289 known species in 51 genera, the ophidiiform family Ophidiidae together with their relatives from the Carapidae (36 species in eight genera) of the same suborder Ophidioidei dominate the deep sea, but some occur also in shallow water habitats. Despite their high species diversity in the deep sea and wide bathymetric distributions, their phylogenetic relationships and evolution remain unexplored due in part to sampling difficulties. Thanks to the biodiversity exploratory program entitled "Tropical Deep-Sea Benthos" and joint efforts between Taiwan and French teams for sampling from different localities across the Indo-West Pacific over the last two decades, we are able to compile comprehensive datasets for investigations. In this study, 59 samples representing 36 of 59 known ophidioid genera are selected and used to construct a multi-gene dataset to infer the phylogenetic relationships of ophidioid fishes and their relatives. Our results reveal that the Ophidiidae forms a paraphyletic group with respect to the Carapidae. The four main clades of Ophidioidei resolved are the (1) clade comprising species from the subfamily Brotulinae; (2) clade that includes species in the genera Acanthonus and Xyelacyba; (3) clade grouping Hypopleuron caninum with species from the family Carapidae; and (4) clade containing the species in the subfamily Brotulotaenilinae, Neobythitinae (in part), and Ophidiinae. Accordingly, we suggest the following new revisions based on our results and proposed morphological diagnoses. The subfamily Brotulinae should be elevated to the family level. The genera Xyelacyba and probably Tauredophidium (unsampled in this study) should be included in the newly established family Acanthonidae with Acanthonus. The families Carapidae and Ophidiidae are re-defined. Our time-calibrated phylogenetic and ancestral depth reconstructions enable us to clarify the evolutionary history of ophidiiform fishes and infer past patterns of species distributions at different depths. While Ophidiiformes is inferred to have originated in shallow waters around 96.25 million years ago (Mya), the common ancestor to the Ophidioidei is inferred to have invaded the deep sea around 90.22 Mya, the dates coinciding with the global anoxic event of the OAE2. The observed bathymetric distribution patterns in Ophidioidei most likely point to the mesopelagic zone as the center of origin and diversification. This was followed by multiple events of depth transitions or range expansions towards either shallower waters or greater depth zones, which were likely triggered by past climate changes during the Paleogene-Neogene.

Long-read de novo genome assembly of Gulf toadfish (Opsanus beta)

BACKGROUND

The family Batrachoididae are a group of ecologically important teleost fishes with unique life histories, behavior, and physiology that has made them popular model organisms. Batrachoididae remain understudied in the realm of genomics, with only four reference genome assemblies available for the family, with three being highly fragmented and not up to current assembly standards. Among these is the Gulf toadfish, Opsanus beta, a model organism for serotonin physiology which has recently been bred in captivity.

RESULTS

Here we present a new, de novo genome and transcriptome assemblies for the Gulf toadfish using PacBio long read technology. The genome size of the final assembly is 2.1 gigabases, which is among the largest teleost genomes. This new assembly improves significantly upon the currently available reference for Opsanus beta with a final scaffold count of 62, of which 23 are chromosome scale, an N50 of 98,402,768, and a BUSCO completeness score of 97.3%. Annotation with ab initio and transcriptome-based methods generated 41,076 gene models. The genome is highly repetitive, with ~ 70% of the genome composed of simple repeats and transposable elements. Satellite DNA analysis identified potential telomeric and centromeric regions.

CONCLUSIONS

This improved assembly represents a valuable resource for future research using this important model organism and to teleost genomics more broadly.

Reassessment of the phylogenetic position of the spiny-scale pricklefish Hispidoberyx ambagiosus (Teleostei: Hispidoberycidae) based on comparative morphology

The spiny-scale pricklefish Hispidoberyx ambagiosus Kotlyar, 1981, the sole member of the family Hispidoberycidae, is known from only a few specimens collected from tropical waters of the Indo-West Pacific region. Because its phylogenetic position has not been well investigated, the present study describes the osteology, myology and other morphological features of H. ambagiosus, and reassesses the phylogenetic position of the species, and its relationships with related taxa. Many significant characters, including Tominaga's organ, were newly discovered in H. ambagiosus and related taxa. Following a phylogenetic analysis of characters in 80 transformation series, three most parsimonious trees were obtained, with H. ambagiosus inferred as forming a monophyletic group together with Barbourisiidae, Cetomimidae, Gibberichthyidae, Rondeletiidae and Stephanoberycidae. Within this clade, H. ambagiosus was inferred as a sister taxon of a clade including the latter three of the aforementioned families. It is considered that H. ambagiosus retains many primitive features, having fewer derived characters than other species in the clade including the six families.

Novel insights into adaptive evolution based on the unusual AT-skew in Acheilognathus gracilis mitogenome and phylogenetic relationships of bitterling

Acheilognathus gracilis, a bitterling species, distribute in lower reaches of Yangtze River. They are identified as the top-priority bitterling species for conservation as having high evolutionary distinctiveness and are at risk of extinction. In present study, we first sequenced the complete mitogenome of A. gracilis and analyzed its phylogenetic position using 13 PCGs. The A. gracilis mitogenome is 16,774 bp in length, including 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, a control region and the origin of the light strand replication. The overall base composition of A. gracilis in descending order is T 27.9 %, A 27.7 %, C 26.1 % and G 18.3 %, shows a unusual AT-skew with slightly negative. Further investigation revealed A. gracilis uses excess T over A in NADH dehydrogenase 5 (nd5), whereas the most of other bitterlings are biased toward to use A not T, implying there is likely to be unique strategy of adaptive evolution in A. gracilis. We also compared 13 PCGs of 30 bitterling mitogenomes and the results exhibit highly conservative. Phylogenetic trees constructed by 13 PCGs strongly support the monophyly of Acheilognathus and the paraphyly of Rhodeus and Tanakia. Current results will provide valuable information for follow-up research on conservation of species facing with serious population decline and can provide novel insights into the phylogenetic analysis and evolutionary biology research.

Mitogenomic Characterization and Phylogenetic Placement of African Hind, Cephalopholis taeniops: Shedding Light on the Evolution of Groupers (Serranidae: Epinephelinae)

The global exploration of evolutionary trends in groupers, based on mitogenomes, is currently underway. This research extensively investigates the structure of and variations in Cephalopholis species mitogenomes, along with their phylogenetic relationships, focusing specifically on Cephalopholis taeniops from the Eastern Atlantic Ocean. The generated mitogenome spans 16,572 base pairs and exhibits a gene order analogous to that of the ancestral teleost's, featuring 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an AT-rich control region. The mitogenome of C. taeniops displays an AT bias (54.99%), aligning with related species. The majority of PCGs in the mitogenome initiate with the start codon ATG, with the exceptions being COI (GTG) and atp6 (TTG). The relative synonymous codon usage analysis revealed the maximum abundance of leucine, proline, serine, and threonine. The nonsynonymous/synonymous ratios were <1, which indicates a strong negative selection among all PCGs of the Cephalopholis species. In C. taeniops, the prevalent transfer RNAs display conventional cloverleaf secondary structures, except for tRNA-serine (GCT), which lacks a dihydrouracil (DHU) stem. A comparative examination of conserved domains and sequence blocks across various Cephalopholis species indicates noteworthy variations in length and nucleotide diversity. Maximum likelihood, neighbor-joining, and Bayesian phylogenetic analyses, employing the concatenated PCGs and a combination of PCGs + rRNAs, distinctly separate all Cephalopholis species, including C. taeniops. Overall, these findings deepen our understanding of evolutionary relationships among serranid groupers, emphasizing the significance of structural considerations in mitogenomic analyses.

The complete mitochondrial genome of the Mexican-endemic cavefish Ophisternon infernale (Synbranchiformes, Synbranchidae): insights on patterns of selection and implications for synbranchiform phylogenetics

Ophisternoninfernale is one of the 200+ troglobitic fish species worldwide, and one of the two cave-dwelling fishes endemic to the karstic aquifer of the Yucatán Peninsula, Mexico. Because of its elusive nature and the relative inaccessibility of its habitat, there is virtually no genetic information on this enigmatic fish. Herein we report the complete mitochondrial genome of O.infernale, which overall exhibits a configuration comparable to that of other synbranchiforms as well as of more distantly related teleosts. The K_A/K_S ratio indicates that most mtDNA PCGs in synbranchiforms have evolved under strong purifying selection, preventing major structural and functional protein changes. The few instances of PCGs under positive selection might be related to adaptation to decreased oxygen availability. Phylogenetic analysis of mtDNA comparative data from synbranchiforms and closely related taxa (including the indostomid Indostomusparadoxus) corroborate the notion that indostomids are more closely related to synbranchiforms than to gasterosteoids, but without rendering the former paraphyletic. Our phylogenetic results also suggest that New World species of Ophisternon might be more closely related to Synbranchus than to the remaining Ophisternon species. This novel phylogenetic hypothesis, however, should be further tested in the context of a comprehensive systematic study of the group.

The Emerging Phylogenetic Perspective on the Evolution of Actinopterygian Fishes

The emergence of a new phylogeny of ray-finned fishes at the turn of the twenty-first century marked a paradigm shift in understanding the evolutionary history of half of living vertebrates. We review how the new ray-finned fish phylogeny radically departs from classical expectations based on morphology. We focus on evolutionary relationships that span the backbone of ray-finned fish phylogeny, from the earliest divergences among teleosts and nonteleosts to the resolution of major lineages of Percomorpha. Throughout, we feature advances gained by the new phylogeny toward a broader understanding of ray-finned fish evolutionary history and the implications for topics that span from the genetics of human health to reconsidering the concept of living fossils. Additionally, we discuss conceptual challenges that involve reconciling taxonomic classification with phylogenetic relationships and propose an alternate higher-level classification for Percomorpha. Our review highlights remaining areas of phylogenetic uncertainty and opportunities for comparative investigations empowered by this new phylogenetic perspective on ray-finned fishes.

Complete vertebrate mitogenomes reveal widespread repeats and gene duplications

BACKGROUND

Modern sequencing technologies should make the assembly of the relatively small mitochondrial genomes an easy undertaking. However, few tools exist that address mitochondrial assembly directly.

RESULTS

As part of the Vertebrate Genomes Project (VGP) we develop mitoVGP, a fully automated pipeline for similarity-based identification of mitochondrial reads and de novo assembly of mitochondrial genomes that incorporates both long (> 10 kbp, PacBio or Nanopore) and short (100-300 bp, Illumina) reads. Our pipeline leads to successful complete mitogenome assemblies of 100 vertebrate species of the VGP. We observe that tissue type and library size selection have considerable impact on mitogenome sequencing and assembly. Comparing our assemblies to purportedly complete reference mitogenomes based on short-read sequencing, we identify errors, missing sequences, and incomplete genes in those references, particularly in repetitive regions. Our assemblies also identify novel gene region duplications. The presence of repeats and duplications in over half of the species herein assembled indicates that their occurrence is a principle of mitochondrial structure rather than an exception, shedding new light on mitochondrial genome evolution and organization.

CONCLUSIONS

Our results indicate that even in the "simple" case of vertebrate mitogenomes the completeness of many currently available reference sequences can be further improved, and caution should be exercised before claiming the complete assembly of a mitogenome, particularly from short reads alone.

Development of the caudal-fin skeleton reveals multiple convergent fusions within Atherinomorpha

Background

The caudal fin of teleosts is a highly diverse morphological structure and a valuable source of information for comparative analyses. Within the Atherinomorpha a high variation of conditions of the caudal-fin skeleton can be found. These range from complex but basal configurations to simple yet derived configurations. When comparing atherinomorph taxa, it is often difficult to decide on the homology of skeletal elements if only considering adult specimens. However, observing the development of caudal-fin skeletons allows one to evaluate complex structures, reveal homologies and developmental patterns, and even reconstruct the grundplan of the examined taxa.

Results

We studied the development of the caudal-fin skeleton in different atheriniform, beloniform and cyprinodontiform species using cleared and stained specimens. Subsequently we compared the development to find similarities and differences in terms of 1) which structures are formed and 2) which structures fuse during ontogeny. For many structures, i.e., the parhypural, the epural(s), the haemal and neural spines of the preural centra and the uroneural, there were either no or only minor differences visible between the three taxa. However, the development of the hypurals revealed a high variation of fusions within different taxa that partly occurred independently in atheriniforms, beloniforms and cyprinodontiforms. Moreover, comparing the development of the ural centra exposed two ways of formation of the compound centrum: 1) in atheriniforms and the beloniforms Oryzias and Hyporhamphus limbatus two ural centra develop and fuse during ontogeny while 2) in cyprinodontiforms and Exocoetidae (Beloniformes) only a single ural centrum is formed during ontogeny.

Conclusions

We were able to reconstruct the grundplan of the developmental pattern of the caudal-fin skeleton of the Atheriniformes, Beloniformes and Cyprinodontiformes as well as their last common ancestors. We found two developmental modes of the compound centrum within the Atherinomorpha, i.e., the fusion of two developing ural centra in atheriniforms and beloniforms and the development of only one ural centrum in cyprinodontiforms. Further differences and similarities for the examined taxa are discussed, resulting in the hypothesis that the caudal-fin development of a last common ancestor to all atherinomorphs is very much similar to that of extant atheriniforms.

Confronting Sources of Systematic Error to Resolve Historically Contentious Relationships: A Case Study Using Gadiform Fishes (Teleostei, Paracanthopterygii, Gadiformes)

Reliable estimation of phylogeny is central to avoid inaccuracy in downstream macroevolutionary inferences. However, limitations exist in the implementation of concatenated and summary coalescent approaches, and Bayesian and full coalescent inference methods may not yet be feasible for computation of phylogeny using complicated models and large data sets. Here, we explored methodological (e.g., optimality criteria, character sampling, model selection) and biological (e.g., heterotachy, branch length heterogeneity) sources of systematic error that can result in biased or incorrect parameter estimates when reconstructing phylogeny by using the gadiform fishes as a model clade. Gadiformes include some of the most economically important fishes in the world (e.g., Cods, Hakes, and Rattails). Despite many attempts, a robust higher-level phylogenetic framework was lacking due to limited character and taxonomic sampling, particularly from several species-poor families that have been recalcitrant to phylogenetic placement. We compiled the first phylogenomic data set, including 14,208 loci (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$>$\end{document}2.8 M bp) from 58 species representing all recognized gadiform families, to infer a time-calibrated phylogeny for the group. Data were generated with a gene-capture approach targeting coding DNA sequences from single-copy protein-coding genes. Species-tree and concatenated maximum-likelihood (ML) analyses resolved all family-level relationships within Gadiformes. While there were a few differences between topologies produced by the DNA and the amino acid data sets, most of the historically unresolved relationships among gadiform lineages were consistently well resolved with high support in our analyses regardless of the methodological and biological approaches used. However, at deeper levels, we observed inconsistency in branch support estimates between bootstrap and gene and site coefficient factors (gCF, sCF). Despite numerous short internodes, all relationships received unequivocal bootstrap support while gCF and sCF had very little support, reflecting hidden conflict across loci. Most of the gene-tree and species-tree discordance in our study is a result of short divergence times, and consequent lack of informative characters at deep levels, rather than incomplete lineage sorting. We use this phylogeny to establish a new higher-level classification of Gadiformes as a way of clarifying the evolutionary diversification of the order. We recognize 17 families in five suborders: Bregmacerotoidei, Gadoidei, Ranicipitoidei, Merluccioidei, and Macrouroidei (including two subclades). A time-calibrated analysis using 15 fossil taxa suggests that Gadiformes evolved \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\sim $\end{document}79.5 Ma in the late Cretaceous, but that most extant lineages diverged after the Cretaceous–Paleogene (K-Pg) mass extinction (66 Ma). Our results reiterate the importance of examining phylogenomic analyses for evidence of systematic error that can emerge as a result of unsuitable modeling of biological factors and/or methodological issues, even when data sets are large and yield high support for phylogenetic relationships. [Branch length heterogeneity; Codfishes; commercial fish species; Cretaceous-Paleogene (K-Pg); heterotachy; systematic error; target enrichment.]

The caudal skeleton of Batrachoidiformes (Teleostei: Percomorphacea): a study of morphological diversity, intraspecific variation, and phylogenetic inferences

The caudal-fin skeleton is a primary data source for systematics of fishes, with characters from this complex being proposed as synapomorphies at many taxonomic levels. Batrachoidiformes is recognized as monophyletic, although intraordinal relationships are unclear. Likewise, interrelationships of Batrachoidiformes to other percomorphs are not well established. The caudal skeleton of Batrachoidiformes has not been thoroughly studied and is poorly represented in recent phylogenetic analyses. In this study, we examined the caudal-fin skeleton of 55 of the 82 species and 22 of the 23 genera of Batrachoidiformes, emphasizing the detection of intraspecific variation to recognize morphological characters with phylogenetic significance. Intraspecific variation is high, especially in the shape of epurals and the parhypural flange. A dorsal prezygapophysis on the first ural centrum and the acute articular edge of the parhypural flange are interpreted as putative synapomorphies of Porichthyinae. The anterior epural supporting two procurrent fin rays is found only in some Halophryninae, but is absent in Allenbatrachus, Batrachomeus, Batrichthys and Halophryne. Among Batrachoidiformes, a hypurapophysis-like process on the first ural centrum is found in Thalassophryninae and Barchatus, Batrichthys, Bifax, Chatrabus, Colletteichthys, Halobatrachus, Perulibatrachus and Riekertia. Caudal-fin ray counts are phylogenetically informative at several taxonomic levels. Distal caudal cartilages are described for Batrachoidiformes for the first time.

Phylogenetic relationships within the primitive acanthomorph fish genus Polymixia, with changes to species composition and geographic distributions

The genus Polymixia is the only survivor of a Late Cretaceous marine fish radiation and is often said to be the most primitive living acanthomorph (i.e., Polymixia possesses the greatest number of primitive character states for Acanthomorpha). Recent studies, including this one, place Polymixia as the sister to all other Paracanthopterygii. Despite its importance, most species of Polymixia are extremely difficult to discriminate on the basis of morphology. As a result, the number of valid species is uncertain. Moreover, there has never been a phylogenetic analysis of the genus. Thus, a molecular phylogenetic study was needed to clarify species boundaries and to resolve relationships within the genus. Tissue or DNA samples backed by museum vouchers were obtained for most species, with additional samples from new geographic areas representing specimens with distinctively different meristics and uncertain identifications. Seven loci (five nuclear and two mitochondrial) were sequenced, from which Bayesian and maximum-likelihood trees were generated. Results reveal nine species-level clades, of which five represent previously known species (Polymixia berndti, P. japonica, P. longispina, P. lowei, and P. nobilis). Surprisingly, results also reveal four previously unknown species-level clades, one close to P. lowei, one close to P. nobilis, and two new species clades related to P. japonica. The species clades are distinguished by their phylogenetic histories, sequence differences, geographic distributions, and morphologies. The clade containing P. berndti is recovered as the sister to all other species of Polymixia. Its genetic variability suggests that it might contain two or more species and it is referred to here as a “species complex”. Polymixia nobilis, the type species, was previously thought to be restricted to the Atlantic, but is now shown to be widespread in the Pacific and possibly in the Indian Ocean. Specimens from waters off Australia identified as P. busakhini actually belong to P. nobilis. In contrast, P. japonica is confirmed only in the area near Japan and the East China Sea; other more distant records are misidentifications. Wide (antipodal) geographic distributions are seen in several clades, including P. nobilis, the P. berndti species complex, and the P. japonica species group. The new phylogeny helps explain the evolution of some morphological characters previously used to distinguish groups of species, particularly dorsal-fin soft-ray count, shape of rows of scale ctenii, and number of pyloric caeca.

Phylogenetic relationships among tribes of the green lacewing subfamily Chrysopinae recovered based on mitochondrial phylogenomics

Chrysopidae (green lacewings) is the second largest family in Neuroptera, and it includes medium-size lacewings largely recognized by the presence of golden-colored eyes, bright green bodies and delicate wings with dense venation patterns. The subfamily Chrysopinae includes 97% of the species diversity in the family and it is currently divided into four tribes: Ankylopterygini, Belonopterygini, Chrysopini and Leucochrysini. Here we sequenced and annotated the nearly complete mitochondrial genomes of four species of each these tribes: Abachrysa eureka, Italochrysa insignis, Leucochrysa pretiosa, Parankyloteryx sp. We then reconstructed the phylogenetic relationships with estimated divergence times among tribes of Chrysopinae based on the mt genomic data. Our results suggest that Chrysopinae sans Nothancyla verreauxi evolved as two reciprocally monophyletic lineages formed by stem members of the tribes Leucochrysini plus Belonopterygini on one hand, and the stem members of Ankylopterygini plus Chrysopini on the other. Our estimations of divergence times place the diversification of stem Chrysopinae into the extant tribes during the Middle Jurassic to Late Cretaceous. The relatively young ages previously estimated for the green lacewing divergences were probably underestimated due to false inferences of homology between non-sister taxa that are later correctly identified as homoplasy after more taxa are added.

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.

Fin modules: an evolutionary perspective on appendage disparity in basal vertebrates

BACKGROUND

Fishes are extremely speciose and also highly disparate in their fin configurations, more specifically in the number of fins present as well as their structure, shape, and size. How they achieved this remarkable disparity is difficult to explain in the absence of any comprehensive overview of the evolutionary history of fish appendages. Fin modularity could provide an explanation for both the observed disparity in fin configurations and the sequential appearance of new fins. Modularity is considered as an important prerequisite for the evolvability of living systems, enabling individual modules to be optimized without interfering with others. Similarities in developmental patterns between some of the fins already suggest that they form developmental modules during ontogeny. At a macroevolutionary scale, these developmental modules could act as evolutionary units of change and contribute to the disparity in fin configurations. This study addresses fin disparity in a phylogenetic perspective, while focusing on the presence/absence and number of each of the median and paired fins.

RESULTS

Patterns of fin morphological disparity were assessed by mapping fin characters on a new phylogenetic supertree of fish orders. Among agnathans, disparity in fin configurations results from the sequential appearance of novel fins forming various combinations. Both median and paired fins would have appeared first as elongated ribbon-like structures, which were the precursors for more constricted appendages. Among chondrichthyans, disparity in fin configurations relates mostly to median fin losses. Among actinopterygians, fin disparity involves fin losses, the addition of novel fins (e.g., the adipose fin), and coordinated duplications of the dorsal and anal fins. Furthermore, some pairs of fins, notably the dorsal/anal and pectoral/pelvic fins, show non-independence in their character distribution, supporting expectations based on developmental and morphological evidence that these fin pairs form evolutionary modules.

CONCLUSIONS

Our results suggest that the pectoral/pelvic fins and the dorsal/anal fins form two distinct evolutionary modules, and that the latter is nested within a more inclusive median fins module. Because the modularity hypotheses that we are testing are also supported by developmental and variational data, this constitutes a striking example linking developmental, variational, and evolutionary modules.

New insights on the sister lineage of percomorph fishes with an anchored hybrid enrichment dataset

Percomorph fishes represent over 17,100 species, including several model organisms and species of economic importance. Despite continuous advances in the resolution of the percomorph Tree of Life, resolution of the sister lineage to Percomorpha remains inconsistent but restricted to a small number of candidate lineages. Here we use an anchored hybrid enrichment (AHE) dataset of 132 loci with over 99,000 base pairs to identify the sister lineage of percomorph fishes. Initial analyses of this dataset failed to recover a strongly supported sister clade to Percomorpha, however, scrutiny of the AHE dataset revealed a bias towards high GC content at fast-evolving codon partitions (GC bias). By combining several existing approaches aimed at mitigating the impacts of convergence in GC bias, including RY coding and analyses of amino acids, we consistently recovered a strongly supported clade comprised of Holocentridae (squirrelfishes), Berycidae (Alfonsinos), Melamphaidae (bigscale fishes), Cetomimidae (flabby whalefishes), and Rondeletiidae (redmouth whalefishes) as the sister lineage to Percomorpha. Additionally, implementing phylogenetic informativeness (PI) based metrics as a filtration method yielded this same topology, suggesting PI based approaches will preferentially filter these fast-evolving regions and act in a manner consistent with other phylogenetic approaches aimed at mitigating GC bias. Our results provide a new perspective on a key issue for studies investigating the evolutionary history of more than one quarter of all living species of vertebrates.

The Wide Distribution and Change of Target Specificity of R2 Non-LTR Retrotransposons in Animals

Transposons, or transposable elements, are the major components of genomes in most eukaryotes. Some groups of transposons have developed target specificity that limits the integration sites to a specific nonessential sequence or a genomic region to avoid gene disruption caused by insertion into an essential gene. R2 is one of the most intensively investigated groups of sequence-specific non-LTR retrotransposons and is inserted at a specific site inside of 28S ribosomal RNA (rRNA) genes. R2 is known to be distributed among at least six animal phyla even though its occurrence is reported to be patchy. Here, in order to obtain a more detailed picture of the distribution of R2, we surveyed R2 using both in silico screening and degenerate PCR, particularly focusing on actinopterygian fish. We found two families of the R2C lineage from vertebrates, although it has previously only been found in platyhelminthes. We also revealed the apparent movement of insertion sites of a lineage of actinopterygian R2, which was likely concurrent with the acquisition of a 28S rRNA-derived sequence in their 3' UTR. Outside of actinopterygian fish, we revealed the maintenance of a single R2 lineage in birds; the co-existence of four lineages of R2 in the leafcutter bee Megachile rotundata; the first examples of R2 in Ctenophora, Mollusca, and Hemichordata; and two families of R2 showing no target specificity. These findings indicate that R2 is relatively stable and universal, while differences in the distribution and maintenance of R2 lineages probably reflect characteristics of some combination of both R2 lineages and host organisms.

Structure and variation of the mitochondrial genome of fishes

Background

The mitochondrial (mt) genome has been used as an effective tool for phylogenetic and population genetic analyses in vertebrates. However, the structure and variability of the vertebrate mt genome are not well understood. A potential strategy for improving our understanding is to conduct a comprehensive comparative study of large mt genome data. The aim of this study was to characterize the structure and variability of the fish mt genome through comparative analysis of large datasets.

Results

An analysis of the secondary structure of proteins for 250 fish species (248 ray-finned and 2 cartilaginous fishes) illustrated that cytochrome c oxidase subunits (COI, COII, and COIII) and a cytochrome bc1 complex subunit (Cyt b) had substantial amino acid conservation. Among the four proteins, COI was the most conserved, as more than half of all amino acid sites were invariable among the 250 species. Our models identified 43 and 58 stems within 12S rRNA and 16S rRNA, respectively, with larger numbers than proposed previously for vertebrates. The models also identified 149 and 319 invariable sites in 12S rRNA and 16S rRNA, respectively, in all fishes. In particular, the present result verified that a region corresponding to the peptidyl transferase center in prokaryotic 23S rRNA, which is homologous to mt 16S rRNA, is also conserved in fish mt 16S rRNA. Concerning the gene order, we found 35 variations (in 32 families) that deviated from the common gene order in vertebrates. These gene rearrangements were mostly observed in the area spanning the ND5 gene to the control region as well as two tRNA gene cluster regions (IQM and WANCY regions). Although many of such gene rearrangements were unique to a specific taxon, some were shared polyphyletically between distantly related species.

Conclusions

Through a large-scale comparative analysis of 250 fish species mt genomes, we elucidated various structural aspects of the fish mt genome and the encoded genes. The present results will be important for understanding functions of the mt genome and developing programs for nucleotide sequence analysis. This study demonstrated the significance of extensive comparisons for understanding the structure of the mt genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3054-y) contains supplementary material, which is available to authorized users.

Characterization of the Complete Mitochondrial Genome Sequence of the Globose Head Whiptail Cetonurus globiceps (Gadiformes: Macrouridae) and Its Phylogenetic Analysis

The particular environmental characteristics of deep water such as its immense scale and high pressure systems, presents technological problems that have prevented research to broaden our knowledge of deep-sea fish. Here, we described the mitogenome sequence of a deep-sea fish, Cetonurus globiceps. The genome is 17,137 bp in length, with a standard set of 22 transfer RNA genes (tRNAs), two ribosomal RNA genes, 13 protein-coding genes, and two typical non-coding control regions. Additionally, a 70 bp tRNA(Thr)-tRNA(Pro) intergenic spacer is present. The C. globiceps mitogenome exhibited strand-specific asymmetry in nucleotide composition. The AT-skew and GC-skew values in the whole genome of C. globiceps were 0 and -0.2877, respectively, revealing that the H-strand had equal amounts of A and T and that the overall nucleotide composition was C skewed. All of the tRNA genes could be folded into cloverleaf secondary structures, while the secondary structure of tRNA(Ser(AGY)) lacked a discernible dihydrouridine stem. By comparing this genome sequence with the recognition sites in teleost species, several conserved sequence blocks were identified in the control region. However, the GTGGG-box, the typical characteristic of conserved sequence block E (CSB-E), was absent. Notably, tandem repeats were identified in the 3' portion of the control region. No similar repetitive motifs are present in most of other gadiform species. Phylogenetic analysis based on 12 protein coding genes provided strong support that C. globiceps was the most derived in the clade. Some relationships however, are in contrast with those presented in previous studies. This study enriches our knowledge of mitogenomes of the genus Cetonurus and provides valuable information on the evolution of Macrouridae mtDNA and deep-sea fish.

Phylogenetic placement of enigmatic percomorph families (Teleostei: Percomorphaceae)

Percomorphs are a large and diverse group of spiny-finned fishes that have come to be known as the "bush at the top" due to their persistent lack of phylogenetic resolution. Recently, the broader Euteleost Tree of Life project (EToL) inferred a well-supported phylogenetic hypothesis that groups the diversity of percomorphs into nine well-supported series (supraordinal groups): Ophidiaria, Batrachoidaria, Gobiaria, Syngnatharia, Pelagiaria, Anabantaria, Carangaria, Ovalentaria, and Eupercaria. The EToL also provided, for the first time, a monophyletic definition of Perciformes - the largest order of vertebrates. Despite significant progress made in accommodating the diversity of percomorph taxa into major clades, some 62 families (most previously placed in "Perciformes", as traditionally defined) were not examined by the EToL. Here, we provide evidence for the phylogenetic affinities of 10 of those 62 families, seven of which have largely remained enigmatic. This expanded taxonomic sampling also provides further support for the nine EToL supraordinal series. We examined sequences from 21 genes previously used by the EToL and added two fast-evolving mitochondrial markers in an attempt to increase resolution within the rapid percomorph radiations. We restricted the taxonomic sampling to 1229 percomorph species, including expanded sampling from recent studies. Results of maximum likelihood analysis revealed that bathyclupeids (Bathyclupeidae), galjoen fishes (Dichistiidae), kelpfishes (Chironemidae), marblefishes (Aplodactylidae), trumpeters (Latridae), barbeled grunters (Hapalogenyidae), slopefishes (Symphysanodontidae), and picarel porgies (formerly Centracanthidae) are members of the series Eupercaria ("new bush at the top"). The picarel porgies and porgies (Sparidae) are now placed in the same family (Sparidae). Our analyses suggest a close affinity between the orders Spariformes (including Lethrinidae, Nemipteridae and Sparidae) and Lobotiformes (including the tripletails or Lobotidae, the barbeled grunters, and tigerperches or Datnioididae), albeit support for this group is low. None of the newly examined families belong in the order Perciformes, as recently defined. Finally, we confirm results from other recent studies that place the Australasian salmons (Arripidae) within Pelagiaria, and the false trevallies (Lactariidae) close to flatfishes, jacks, and trevallies, within Carangaria.

Are 100 enough? Inferring acanthomorph teleost phylogeny using Anchored Hybrid Enrichment

BACKGROUND

The past decade has witnessed remarkable progress towards resolution of the Tree of Life. However, despite the increased use of genomic scale datasets, some phylogenetic relationships remain difficult to resolve. Here we employ anchored phylogenomics to capture 107 nuclear loci in 29 species of acanthomorph teleost fishes, with 25 of these species sampled from the recently delimited clade Ovalentaria. Previous studies employing multilocus nuclear exon datasets have not been able to resolve the nodes at the base of the Ovalentaria tree with confidence. Here we test whether a phylogenomic approach will provide better support for these nodes, and if not, why this may be.

RESULTS

After using a novel method to account for paralogous loci, we estimated phylogenies with maximum likelihood and species tree methods using DNA sequence alignments of over 80,000 base pairs. Several key relationships within Ovalentaria are well resolved, including 1) the sister taxon relationship between Cichlidae and Pholidichthys, 2) a clade containing blennies, grammas, clingfishes, and jawfishes, and 3) monophyly of Atherinomorpha (topminnows, flyingfishes, and silversides). However, many nodes in the phylogeny associated with the early diversification of Ovalentaria are poorly resolved in several analyses. Through the use of rarefaction curves we show that limited phylogenetic resolution among the earliest nodes in the Ovalentaria phylogeny does not appear to be due to a deficiency of data, as average global node support ceases to increase when only 1/3rd of the sampled loci are used in analyses. Instead this lack of resolution may be driven by model misspecification as a Bayesian mixed model analysis of the amino acid dataset provided good support for parts of the base of the Ovalentaria tree.

CONCLUSIONS

Although it does not appear that the limited phylogenetic resolution among the earliest nodes in the Ovalentaria phylogeny is due to a deficiency of data, it may be that both stochastic and systematic error resulting from model misspecification play a role in the poor resolution at the base of the Ovalentaria tree as a Bayesian approach was able to resolve some of the deeper nodes, where the other methods failed.

Deep-sea bigscales, pricklefishes, gibberfishes and whalefishes (Teleostei: Stephanoberycoidei) off Brazil: new records, range extensions for the south-western Atlantic Ocean and remarks on the taxonomy of Poromitra

The Stephanoberycoidei includes 23 genera and c. 94 species of deep-sea teleosts commonly known as bigscales, pricklefishes, gibberfishes and whalefishes. Stephanoberycoidei is one of the least known groups of deep-sea fishes, in spite of their apparent relative abundance in meso and bathypelagic depths. Nine species of the Stephanoberycoidei are reported here for the first time in Brazilian waters, and most of them represent new range extensions for the south-western Atlantic Ocean. Those species are Melamphaes polylepis, Melamphaes typhlops, Poromitra sp. and Scopeloberyx robustus (Melamphaidae), Acanthochaenus luetkenii and Stephanoberyx monae (Stephanoberycidae), Rondeletia bicolor and Rondeletia loricata (Rondeletiidae) and Gyrinomimus sp. (Cetomimidae). Occurrences of the pricklefish Scopelogadus mizolepis (Melamphaidae), the gibberfish Gibberichthys pumilus (Gibberichthyidae) and the velvet whalefish Barbourisia rufa (Barbourisiidae) are confirmed in the Brazilian exclusive economic zone, but previously published records of Poromitra capito and Melamphaes simus (Melamphaidae) in the region most likely represent misidentifications. Validities of the recently described Poromitra kukuevi and Poromitra indooceanica are discussed in light of new specimens of the genus collected in the south-western Atlantic Ocean. An identification key for the 13 species of Stephanoberycoidei reported off Brazil is also provided.

The infrabranchial musculature and its bearing on the phylogeny of percomorph fishes (Osteichthyes: Teleostei)

The muscles serving the ventral portion of the gill arches ( = infrabranchial musculature) are poorly known in bony fishes. A comparative analysis of the infrabranchial muscles in the major percomorph lineages reveals a large amount of phylogenetically-relevant information. Characters derived from this anatomical system are identified and discussed in light of current hypotheses of phylogenetic relationships among percomorphs. New evidence supports a sister-group relationship between the Batrachoidiformes and Lophiiformes and between the Callionymoidei and Gobiesocoidei. Investigated data also corroborate the existence of two monophyletic groups, one including the Pristolepididae, Badidae, and Nandidae, and a second clade consisting of all non-amarsipid stromateiforms. New synapomorphies are proposed for the Atherinomorphae, Blenniiformes, Lophiiformes, Scombroidei (including Sphyraenidae), and Gobiiformes. Within the latter order, the Rhyacichthyidae and Odontobutidae are supported as the successive sister families of all remaining gobiiforms. The present analysis further confirms the validity of infrabranchial musculature characters previously proposed to support the grouping of the Mugiliformes with the Atherinomorphae and the monophyly of the Labriformes with the possible inclusion of the Pholidichthyiformes. Interestingly, most hypotheses of relationships supported by the infrabranchial musculature have been advanced by preceding anatomists on the basis of distinct data sources, but were never recovered in recent molecular phylogenies. These conflicts clearly indicate the current unsatisfactory resolution of the higher-level phylogeny of percomorphs.

Mitogenomic circumscription of a novel percomorph fish clade mainly comprising "Syngnathoidei" (Teleostei)

Percomorpha, comprising about 60% of modern teleost fishes, has been described as the "(unresolved) bush at the top" of the tree, with its intrarelationships still being ambiguous owing to huge diversity (>15,000 species). Recent molecular phylogenetic studies based on extensive taxon and character sampling, however, have revealed a number of unexpected clades of Percomorpha, and one of which is composed of Syngnathoidei (seahorses, pipefishes, and their relatives) plus several groups distributed across three different orders. To circumscribe the clade more definitely, we sampled several candidate taxa with reference to the previous studies and newly determined whole mitochondrial genome (mitogenome) sequences for 16 percomorph species across syngnathoids, dactylopterids, and their putatively closely-related fishes (Mullidae, Callionymoidei, Malacanthidae). Unambiguously aligned sequences (13,872 bp) from those 16 species plus 78 percomorphs and two outgroups (total 96 species) were subjected to partitioned Bayesian and maximum likelihood analyses. The resulting trees revealed a highly supported clade comprising seven families in Syngnathoidei (Gasterosteiformes), Dactylopteridae (Scorpaeniformes), Mullidae in Percoidei and two families in Callionymoidei (Perciformes). We herein proposed to call this clade "Syngnathiformes" following the latest nuclear DNA studies with some revisions on the included families.

Phylogeny and biogeography of a shallow water fish clade (Teleostei: Blenniiformes)

BACKGROUND

The Blenniiformes comprises six families, 151 genera and nearly 900 species of small teleost fishes closely associated with coastal benthic habitats. They provide an unparalleled opportunity for studying marine biogeography because they include the globally distributed families Tripterygiidae (triplefin blennies) and Blenniidae (combtooth blennies), the temperate Clinidae (kelp blennies), and three largely Neotropical families (Labrisomidae, Chaenopsidae, and Dactyloscopidae). However, interpretation of these distributional patterns has been hindered by largely unresolved inter-familial relationships and the lack of evidence of monophyly of the Labrisomidae.

RESULTS

We explored the phylogenetic relationships of the Blenniiformes based on one mitochondrial (COI) and four nuclear (TMO-4C4, RAG1, Rhodopsin, and Histone H3) loci for 150 blenniiform species, and representative outgroups (Gobiesocidae, Opistognathidae and Grammatidae). According to the consensus of Bayesian Inference, Maximum Likelihood, and Maximum Parsimony analyses, the monophyly of the Blenniiformes and the Tripterygiidae, Blenniidae, Clinidae, and Dactyloscopidae is supported. The Tripterygiidae is the sister group of all other blennies, and the Blenniidae is the sister group of the remaining blennies. The monophyly of the Labrisomidae is supported with the exclusion of the Cryptotremini and inclusion of Stathmonotus, and we elevate two subgenera of Labrisomus to establish a monophyletic classification within the family. The monophyly of the Chaenopsidae is supported with the exclusion of Stathmonotus (placed in the Stathmonotini) and Neoclinus and Mccoskerichthys (placed in the Neoclinini). The origin of the Blenniiformes was estimated in the present-day IndoPacific region, corresponding to the Tethys Sea approximately 60.3 mya. A largely Neotropical lineage including the Labrisomidae, Chaenopsidae and Dactyloscopidae (node IV) evolved around 37.6 mya when the Neotropics were increasingly separated from the IndoPacific, but well before the closure of the Tethys Sea.

CONCLUSIONS

Relationships recovered in this study are similar to those of earlier analyses within the Clinidae and Chaenopsidae, and partially similar within the Blenniidae, but tripterygiid relationships remain poorly resolved. We present the first comprehensive phylogenetic hypothesis for a monophyletic Labrisomidae with five tribes (Labrisomini, Mnierpini, Paraclinini, Stathmonotini and Starksiini). Global distributions of blenny genera included in our analysis support the evolution of a largely Neotropical clade whose closest relatives (clinids and cryptotremines) are temperate in distribution.

A comprehensive description and evolutionary analysis of 22 grouper (perciformes, epinephelidae) mitochondrial genomes with emphasis on two novel genome organizations

Groupers of the family Epinephelidae are a diverse and economically valuable group of reef fishes. To investigate the evolution of their mitochondrial genomes we characterized and compared these genomes among 22 species, 17 newly sequenced. Among these fishes we identified three distinct genome organizations, two of them never previously reported in vertebrates. In 19 of these species, mitochondrial genomes followed the typical vertebrate canonical organization with 13 protein-coding genes, 22 tRNAs, two rRNAs, and a non-coding control region. Differing from this, members of genus Variola have an extra tRNA-Ile between tRNA-Val and 16S rRNA. Evidence suggests that this evolved from tRNA-Val via a duplication event due to slipped strand mispairing during replication. Additionally, Cephalopholisargus has an extra tRNA-Asp in the midst of the control region, likely resulting from long-range duplication of the canonical tRNA-Asp through illicit priming of mitochondrial replication by tRNAs. Along with their gene contents, we characterized the regulatory elements of these mitochondrial genomes' control regions, including putative termination-associated sequences and conserved sequence blocks. Looking at the mitochondrial genomic constituents, rRNA and tRNA are the most conserved, followed by protein-coding genes, and non-coding regions are the most divergent. Divergence rates vary among the protein-coding genes, and the three cytochrome oxidase subunits (COI, II, III) are the most conserved, while NADH dehydrogenase subunit 6 (ND6) and the ATP synthase subunit 8 (ATP8) are the most divergent. We then tested the phylogenetic utility of this new mt genome data using 12 protein-coding genes of 48 species from the suborder Percoidei. From this, we provide further support for the elevation of the subfamily Epinephelinae to family Epinephelidae, the resurrection of the genus Hyporthodus, and the combination of the monotypic genera Anyperodon and Cromileptes to genus Epinephelus, and Aethaloperca to genus Cephalopholis.

Phylogeny and tempo of diversification in the superradiation of spiny-rayed fishes

Spiny-rayed fishes, or acanthomorphs, comprise nearly one-third of all living vertebrates. Despite their dominant role in aquatic ecosystems, the evolutionary history and tempo of acanthomorph diversification is poorly understood. We investigate the pattern of lineage diversification in acanthomorphs by using a well-resolved time-calibrated phylogeny inferred from a nuclear gene supermatrix that includes 520 acanthomorph species and 37 fossil age constraints. This phylogeny provides resolution for what has been classically referred to as the "bush at the top" of the teleost tree, and indicates acanthomorphs originated in the Early Cretaceous. Paleontological evidence suggests acanthomorphs exhibit a pulse of morphological diversification following the end Cretaceous mass extinction; however, the role of this event on the accumulation of living acanthomorph diversity remains unclear. Lineage diversification rates through time exhibit no shifts associated with the end Cretaceous mass extinction, but there is a global decrease in lineage diversification rates 50 Ma that occurs during a period when morphological disparity among fossil acanthomorphs increases sharply. Analysis of clade-specific shifts in diversification rates reveal that the hyperdiversity of living acanthomorphs is highlighted by several rapidly radiating lineages including tunas, gobies, blennies, snailfishes, and Afro-American cichlids. These lineages with high diversification rates are not associated with a single habitat type, such as coral reefs, indicating there is no single explanation for the success of acanthomorphs, as exceptional bouts of diversification have occurred across a wide array of marine and freshwater habitats.

The phylogenetic significance of colour patterns in marine teleost larvae

Ichthyologists, natural-history artists, and tropical-fish aquarists have described, illustrated, or photographed colour patterns in adult marine fishes for centuries, but colour patterns in marine fish larvae have largely been neglected. Yet the pelagic larval stages of many marine fishes exhibit subtle to striking, ephemeral patterns of chromatophores that warrant investigation into their potential taxonomic and phylogenetic significance. Colour patterns in larvae of over 200 species of marine teleosts, primarily from the western Caribbean, were examined from digital colour photographs, and their potential utility in elucidating evolutionary relationships at various taxonomic levels was assessed. Larvae of relatively few basal marine teleosts exhibit erythrophores, xanthophores, or iridophores (i.e. nonmelanistic chromatophores), but one or more of those types of chromatophores are visible in larvae of many basal marine neoteleosts and nearly all marine percomorphs. Whether or not the presence of nonmelanistic chromatophores in pelagic marine larvae diagnoses any major teleost taxonomic group cannot be determined based on the preliminary survey conducted, but there is a trend toward increased colour from elopomorphs to percomorphs. Within percomorphs, patterns of nonmelanistic chromatophores may help resolve or contribute evidence to existing hypotheses of relationships at multiple levels of classification. Mugilid and some beloniform larvae share a unique ontogenetic transformation of colour pattern that lends support to the hypothesis of a close relationship between them. Larvae of some tetraodontiforms and lophiiforms are strikingly similar in having the trunk enclosed in an inflated sac covered with xanthophores, a character that may help resolve the relationships of these enigmatic taxa. Colour patterns in percomorph larvae also appear to diagnose certain groups at the interfamilial, familial, intergeneric, and generic levels. Slight differences in generic colour patterns, including whether the pattern comprises xanthophores or erythrophores, often distinguish species. The homology, ontogeny, and possible functional significance of colour patterns in larvae are discussed. Considerably more investigation of larval colour patterns in marine teleosts is needed to assess fully their value in phylogenetic reconstruction.

Molecular phylogeny and new classification of the genera Eulophias and Zoarchias (PISCES, Zoarcoidei)

Morphological and osteological studies of the Zoarcoidei group have previously been undertaken, but the group (especially the genera Eulophias and Zoarchias) still remains enigmatic. Therefore, we conducted molecular phylogenetic studies on the two genera Eulophias and Zoarchias using two mitochondrial (16S rRNA and COI) and two nuclear genes (RAG2 and RNF213). Our phylogenetic analysis supported the monophyly of the suborder level of the Zoarcoidei, but rejected the previous morphology- and osteology-based classification hypotheses regarding the two genera. Conflict between mtDNA and nDNA phylogenies within the genus Eulophias implies that the genus shows a complicated relationship such as hybridization in the process of the evolutionary history. The genetic distances between the Eulophias (or Zoarchias) and other Zoarcoidei spp. were the greatest, showing different family-level affiliations. In addition, the mtDNA topology showed the two genera were clearly separated from each other as well as from the families Stichaeidae and Zoarcidae. Considering the new molecular phylogeny, we suggest a new classification for the two genera: (1) Eulophias belongs to a new family named as the Eulophiidae; (2) Zoarchias belongs to the family Neozoarcidae (sensu Radchenko et al., 2012b) rather than to Stichaeidae and Zoarcidae.

Mitogenomic sequences and evidence from unique gene rearrangements corroborate evolutionary relationships of myctophiformes (Neoteleostei)

Background

A skewed assemblage of two epi-, meso- and bathypelagic fish families makes up the order Myctophiformes – the blackchins Neoscopelidae and the lanternfishes Myctophidae. The six rare neoscopelids show few morphological specializations whereas the divergent myctophids have evolved into about 250 species, of which many show massive abundances and wide distributions. In fact, Myctophidae is by far the most abundant fish family in the world, with plausible estimates of more than half of the oceans combined fish biomass. Myctophids possess a unique communication system of species-specific photophore patterns and traditional intrafamilial classification has been established to reflect arrangements of photophores. Myctophids present the most diverse array of larval body forms found in fishes although this attribute has both corroborated and confounded phylogenetic hypotheses based on adult morphology. No molecular phylogeny is available for Myctophiformes, despite their importance within all ocean trophic cycles, open-ocean speciation and as an important part of neoteleost divergence. This study attempts to resolve major myctophiform phylogenies from both mitogenomic sequences and corroborating evidence in the form of unique mitochondrial gene order rearrangements.

Results

Mitogenomic evidence from DNA sequences and unique gene orders are highly congruent concerning phylogenetic resolution on several myctophiform classification levels, corroborating evidence from osteology, larval ontogeny and photophore patterns, although the lack of larval morphological characters within the subfamily Lampanyctinae stands out. Neoscopelidae is resolved as the sister family to myctophids with Solivomer arenidens positioned as a sister taxon to the remaining neoscopelids. The enigmatic Notolychnus valdiviae is placed as a sister taxon to all other myctophids and exhibits an unusual second copy of the tRNA-Met gene – a gene order rearrangement reminiscent of that found in the tribe Diaphini although our analyses show it to be independently derived. Most tribes are resolved in accordance with adult morphology although Gonichthyini is found within a subclade of the tribe Myctophini consisting of ctenoid scaled species. Mitogenomic sequence data from this study recognize 10 reciprocally monophyletic lineages within Myctophidae, with five of these clades delimited from additional rearranged gene orders or intergenic non-coding sequences.

Conclusions

Mitogenomic results from DNA sequences and unique gene orders corroborate morphology in phylogeny reconstruction and provide a likely scenario for the phylogenetic history of Myctophiformes. The extent of gene order rearrangements found within the mitochondrial genomes of myctophids is unique for phylogenetic purposes.

The complete mitochondrial genome of the three-spot seahorse, Hippocampus trimaculatus (Teleostei, Syngnathidae)

The complete mitochondrial genome of the three-spot seahorse was sequenced using a polymerase chain reaction-based method. The total length of mitochondrial DNA is 16,535 bp and includes 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region. The mitochondrial gene order of the three-spot seahorse also conforms to the distinctive vertebrate mitochondrial gene order. The base composition of the genome is A (32.7%), T (29.3%), C (23.4%), and G (14.6%) with an A + T-rich hallmark as that of other vertebrate mitochondrial genomes.

Characterization of comparative genome-derived simple sequence repeats for acanthopterygian fishes

Simple sequence repeats (SSRs) have become one of the most popular molecular markers for population genetic studies. The application of SSR markers has often been limited to source species because SSR loci are too labile to be maintained in even closely related species. However, a few extremely conserved SSR loci have been reported. Here, we tested for the presence of conserved SSR loci in acanthopterygian fishes, which include over 14 000 species, by comparing the genome sequences of four acanthopterygian fishes. We also examined the comparative genome-derived SSRs (CG-SSRs) for their transferability across acanthopterygian fishes and their applicability to population genetic analysis. Forty-six SSR loci with conserved flanking regions were detected and examined for their transferability among seven nonacanthopterygian and 27 acanthopterygian fishes. The PCR amplification success rate in nonacanthopterygian fishes was low, ranging from 2.2% to 21.7%, except for Lophius litulon (Lophiiformes; 80.4%). Conversely, the rate in most acanthopterygian fishes exceeded 70.0%. Sequencing of these 46 loci revealed the presence of SSRs suitable for scoring while fragment analysis of 20 loci revealed polymorphisms in most of the acanthopterygian fishes. Population genetic analysis of Cottus pollux (Scorpaeniformes) and Sphaeramia orbicularis (Perciformes) using CG-SSRs showed that these populations did not deviate from linkage equilibrium or Hardy-Weinberg equilibrium. Furthermore, almost no loci showed evidence of null alleles, suggesting that CG-SSRs have strong resolving power for population genetic analysis. Our findings will facilitate the use of these markers in species in which markers remain to be identified.

Visceral anatomy of ocean sunfish (Mola mola (L., 1758), Molidae, Tetraodontiformes) and angler (Lophius piscatorius (L., 1758), Lophiidae, Lophiiformes) investigated by non-invasive imaging techniques

The purpose of this work is to examine the gross visceral anatomy of ocean sunfish and angler using non-invasive imaging techniques: computed tomography imaging (CT) and magnetic resonance imaging (MRI). Similarities and differences in the internal organisation of these two species are verified. Both species lack a swimbladder and present a significant asymmetry in the hepatic lobes, an elongated bile duct terminating close to the stomach, a compact thyroid embedded in a blood lacuna, and very reduced brain and spinal cord. These observations are important in regard to the close relationships between Tetraodontiformes and Lophiiformes, established by several molecular works, but not yet confirmed by morpho-anatomical data. However the occurrence of these features has to be examined in other taxa before phylogenetic hypotheses are proposed.

Integrating multi-origin expression data improves the resolution of deep phylogeny of ray-finned fish (Actinopterygii)

The actinopterygians comprise nearly one-half of all extant vertebrate species and are very important for human well-being. However, the phylogenetic relationships among certain groups within the actinopterygians are still uncertain, and debates about these relationships have continued for a long time. Along with the progress achieved in sequencing technologies, phylogenetic analyses based on multi-gene sequences, termed phylogenomic approaches, are becoming increasingly common and often result in well-resolved and highly supported phylogenetic hypotheses. Based on the transcriptome sequences generated in this study and the extensive expression data currently available from public databases, we obtained alignments of 274 orthologue groups for 26 scientifically and commercially important actinopterygians, representing 17 out of 44 orders within the class Actinopterygii. Using these alignments and probabilistic methods, we recovered relationships between basal actinopterygians and teleosts, among teleosts within protacanthopterygians and related lineages, and also within acanthomorphs. These relationships were recovered with high confidence.

Egg surface structure of the freshwater toadfish Thalassophryne amazonica (Teleostei: Batrachoididae) with information on its distribution and natural habitat

The egg surface structure of Thalassophryne amazonica, a freshwater toadfish from the Amazon basin is described. Eggs of this species show a remarkable, highly unusual system of parallel ridges and intermittent grooves that originate at the equator of the egg and lead to the micropylar pit, at which they end in a spiralling pattern. A similar egg surface structure has so far been described only from a group of Asian anabantoid percomorphs, obviously not closely related to Thalassophryne. This egg surface pattern may enhance fertilization success by guiding sperm to the micropyle. We review museum records for T. amazonica, present an updated map of its occurrence in the Amazon basin, and provide information on its habitat.

Molecular phylogenetics of squirrelfishes and soldierfishes (Teleostei: Beryciformes: Holocentridae): reconciling more than 100 years of taxonomic confusion

Squirrelfishes and soldierfishes (Holocentridae) are among the most conspicuous species in the nocturnal reef fish community. However, there is no clear consensus regarding their evolutionary relationships, which is reflected in a complicated taxonomic history. We collected DNA sequence data from multiple single copy nuclear genes and one mitochondrial gene sampled from over fifty percent of the recognized holocentrid species and infer the first species-level phylogeny of the Holocentridae. Our results strongly support the monophyly of the clades Myripristinae (soldierfishes) and Holocentrinae (squirrelfishes). The molecular phylogenies differ with regard to previous hypotheses of relationships within the Myriprisitinae, resolving a clade of cryptic reef associated and deep water non-reef dwelling lineages (Corniger+Plectrypops+Ostichthys) that is the sister lineage to a monophyletic Myripristis. Within Holocentrinae, Neoniphon and Sargocentron are strongly supported as paraphyletic, while Holocentrus is nested within Sargocentron. Using Bayesian ancestral state reconstruction methods, we demonstrate the taxonomically diagnostic characters for Neoniphon and Sargocentron likely represent character states with a complex evolutionary history that is not reflective of shared common ancestry. We propose a new classification for Holocentrinae, recognizing four lineages that are treated as genera: Sargocentron Fowler, 1904, Holocentrus Scopoli, 1777, Flameo Jordan and Evermann, 1898, and Neoniphon Castelnau, 1875.

The evolution of pharyngognathy: a phylogenetic and functional appraisal of the pharyngeal jaw key innovation in labroid fishes and beyond

The perciform group Labroidei includes approximately 2600 species and comprises some of the most diverse and successful lineages of teleost fishes. Composed of four major clades, Cichlidae, Labridae (wrasses, parrotfishes, and weed whitings), Pomacentridae (damselfishes), and Embiotocidae (surfperches); labroids have been an icon for studies of biodiversity, adaptive radiation, and sexual selection. The success and diversification of labroids have been largely attributed to the presence of a major innovation in the pharyngeal jaw apparatus, pharyngognathy, which is hypothesized to increase feeding capacity and versatility. We present results of large-scale phylogenetic analyses and a survey of pharyngeal jaw functional morphology that allow us to examine the evolution of pharyngognathy in a historical context. Phylogenetic analyses were based on a sample of 188 acanthomorph (spiny-rayed fish) species, primarily percomorphs (perch-like fishes), and DNA sequence data collected from 10 nuclear loci that have been previously used to resolve higher level ray-finned fish relationships. Phylogenies inferred from this dataset using maximum likelihood, Bayesian, and species tree analyses indicate polyphyly of the traditional Labroidei and clearly separate Labridae from the remainder of the traditional labroid lineages (Cichlidae, Embiotocidae, and Pomacentridae). These three "chromide" families grouped within a newly discovered clade of 40 families and more than 4800 species (>27% of percomorphs and >16% of all ray-finned fishes), which we name Ovalentaria for its characteristic demersal, adhesive eggs with chorionic filaments. This fantastically diverse clade includes some of the most species-rich lineages of marine and freshwater fishes, including all representatives of the Cichlidae, Embiotocidae, Pomacentridae, Ambassidae, Gobiesocidae, Grammatidae, Mugilidae, Opistognathidae, Pholidichthyidae, Plesiopidae (including Notograptus), Polycentridae, Pseudochromidae, Atherinomorpha, and Blennioidei. Beyond the discovery of Ovalentaria, this study provides a surprising, but well-supported, hypothesis for a convict-blenny (Pholidichthys) sister group to the charismatic cichlids and new insights into the evolution of pharyngognathy. Bayesian stochastic mapping ancestral state reconstructions indicate that pharyngognathy has evolved at least six times in percomorphs, including four separate origins in members of the former Labroidei, one origin in the Centrogenyidae, and one origin within Beloniformes. Our analyses indicate that all pharyngognathous fishes have a mechanically efficient biting mechanism enabled by the muscular sling and a single lower jaw element. However, a major distinction exists between Labridae, which lacks the widespread, generalized percomorph pharyngeal biting mechanism, and all other pharyngognathous clades, which possess this generalized biting mechanism in addition to pharyngognathy. Our results reveal a remarkable history of pharyngognathy: far from a single origin, it appears to have evolved at least six times, and its status as a major evolutionary innovation is reinforced by it being a synapomorphy for several independent major radiations, including some of the most species rich and ecologically diverse percomorph clades of coral reef and tropical freshwater fishes, Labridae and Cichlidae. [Acanthomorpha; Beloniformes; Centrogenyidae; key innovation; Labroidei; Ovalentaria; pharyngeal jaws; Perciformes.].

A comparative study of nemertean complete mitochondrial genomes, including two new ones for Nectonemertes cf. mirabilis and Zygeupolia rubens, may elucidate the fundamental pattern for the phylum Nemertea

BACKGROUND

The mitochondrial genome is important for studying genome evolution as well as reconstructing the phylogeny of organisms. Complete mitochondrial genome sequences have been reported for more than 2200 metazoans, mainly vertebrates and arthropods. To date, from a total of about 1275 described nemertean species, only three complete and two partial mitochondrial DNA sequences from nemerteans have been published. Here, we report the entire mitochondrial genomes for two more nemertean species: Nectonemertes cf. mirabilis and Zygeupolia rubens.

RESULTS

The sizes of the entire mitochondrial genomes are 15365 bp for N. cf. mirabilis and 15513 bp for Z. rubens. Each circular genome contains 37 genes and an AT-rich non-coding region, and overall nucleotide composition is AT-rich. In both species, there is significant strand asymmetry in the distribution of nucleotides, with the coding strand being richer in T than A and in G than C. The AT-rich non-coding regions of the two genomes have some repeat sequences and stem-loop structures, both of which may be associated with the initiation of replication or transcription. The 22 tRNAs show variable substitution patterns in nemerteans, with higher sequence conservation in genes located on the H strand. Gene arrangement of N. cf. mirabilis is identical to that of Paranemertes cf. peregrina, both of which are Hoplonemertea, while that of Z. rubens is the same as in Lineus viridis, both of which are Heteronemertea. Comparison of the gene arrangements and phylogenomic analysis based on concatenated nucleotide sequences of the 12 mitochondrial protein-coding genes revealed that species with closer relationships share more identical gene blocks.

CONCLUSION

The two new mitochondrial genomes share many features, including gene contents, with other known nemertean mitochondrial genomes. The tRNA families display a composite substitution pathway. Gene order comparison to the proposed ground pattern of Bilateria and some lophotrochozoans suggests that the nemertean ancestral mitochondrial gene order most closely resembles the heteronemertean type. Phylogenetic analysis proposes a sister-group relationship between Hetero- and Hoplonemertea, which supports one of two recent alternative hypotheses of nemertean phylogeny.

Pervasive indels and their evolutionary dynamics after the fish-specific genome duplication

Insertions and deletions (indels) in protein-coding genes are important sources of genetic variation. Their role in creating new proteins may be especially important after gene duplication. However, little is known about how indels affect the divergence of duplicate genes. We here study thousands of duplicate genes in five fish (teleost) species with completely sequenced genomes. The ancestor of these species has been subject to a fish-specific genome duplication (FSGD) event that occurred approximately 350 Ma. We find that duplicate genes contain at least 25% more indels than single-copy genes. These indels accumulated preferentially in the first 40 my after the FSGD. A lack of widespread asymmetric indel accumulation indicates that both members of a duplicate gene pair typically experience relaxed selection. Strikingly, we observe a 30-80% excess of deletions over insertions that is consistent for indels of various lengths and across the five genomes. We also find that indels preferentially accumulate inside loop regions of protein secondary structure and in regions where amino acids are exposed to solvent. We show that duplicate genes with high indel density also show high DNA sequence divergence. Indel density, but not amino acid divergence, can explain a large proportion of the tertiary structure divergence between proteins encoded by duplicate genes. Our observations are consistent across all five fish species. Taken together, they suggest a general pattern of duplicate gene evolution in which indels are important driving forces of evolutionary change.