Are flatfishes (Pleuronectiformes) monophyletic?

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.

Phylogenetic relationships and the origin of New World soles (Teleostei: Pleuronectiformes: Achiridae): The role of estuarine habitats

Even though the monophyletic status of Achiridae has been supported by morphological and molecular data, the interrelationships within the representatives of this family are poorly resolved. In the present study, we carried out the most complete molecular phylogenetic analysis of this group, encompassing all genera and employing both nuclear (Rhodopsin, Recombination activator [Rag 1], Mixed - lineage Leukemia [MLL] and Early Growth Response Protein 3 [EGR3]) and mitochondrial (Cytochrome C Oxidase Subunit I [COI], Cytochrome B [CytB], ATPase 6.8, 16S and 12S RNAr) genes. All topologies based on Maximum Likelihood, Bayesian inferences and Bayesian Inference of the Multispecies Coalescent confirmed the monophyletism of Achiridae, in spite of some incongruences in relation to Achirus mucuri, A. lineatus, Apionichthys finis and Trinectes microphthalmus. In fact, Achirus and Trinectes proved to be non-monophyletic genera while Hypoclinemus mentalis was closely related to A. achirus, suggesting this species should be reevaluated. We provided evidence that Achiridae has first arisen in estuaries (about 23.5 million years ago) and some lineages have evolved independently to either marine or freshwater habitats. Furthermore, we propose a diversification scenario of New World soles involving at least two events of marine incursions during Miocene and Pliocene - Pleistocene associated with natural geographic barriers (Victoria-Trindade chain), the width and exposure of continental shelf and headwater capture along the Amazon basin. Finally, the evolutionary dependence of Achirid soles on estuaries, characterized as highly dynamic environments, has probably driven the recent divergence of many species of Achiridae.

Integration of Maps Enables a Cytogenomics Analysis of the Complete Karyotype in Solea senegalensis

The Pleuronectiformes order, which includes several commercially-important species, has undergone extensive chromosome evolution. One of these species is Solea senegalensis, a flatfish with 2n = 42 chromosomes. In this study, a cytogenomics approach and integration with previous maps was applied to characterize the karyotype of the species. Synteny analysis of S. senegalensis was carried out using two flatfish as a reference: Cynoglossus semilaevis and Scophthalmus maximus. Most S. senegalensis chromosomes (or chromosome arms for metacentrics and submetacentrics) showed a one-to-one macrosyntenic pattern with the other two species. In addition, we studied how repetitive sequences could have played a role in the evolution of S. senegalensis bi-armed (3, and 5–9) and acrocentric (11, 12 and 16) chromosomes, which showed the highest rearrangements compared with the reference species. A higher abundance of TEs (Transposable Elements) and other repeated elements was observed adjacent to telomeric regions on chromosomes 3, 7, 9 and 16. However, on chromosome 11, a greater abundance of DNA transposons was detected in interstitial BACs. This chromosome is syntenic with several chromosomes of the other two flatfish species, suggesting rearrangements during its evolution. A similar situation was also found on chromosome 16 (for microsatellites and low complexity sequences), but not for TEs (retroelements and DNA transposons). These differences in the distribution and abundance of repetitive elements in chromosomes that have undergone remodeling processes during the course of evolution also suggest a possible role for simple repeat sequences in rearranged regions.

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.

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.

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.

Cytogenomics Unveil Possible Transposable Elements Driving Rearrangements in Chromosomes 2 and 4 of Solea senegalensis

Cytogenomics, the integration of cytogenetic and genomic data, has been used here to reconstruct the evolution of chromosomes 2 and 4 of Solea senegalensis. S. senegalensis is a flat fish with a karyotype comprising 2n = 42 chromosomes: 6 metacentric + 4 submetacentric + 8 subtelocentric + 24 telocentric. The Fluorescence in situ Hybridization with Bacterial Artificial Chromosomes (FISH-BAC) technique was applied to locate BACs in these chromosomes (11 and 10 BACs in chromosomes 2 and 4, respectively) and to generate integrated maps. Synteny analysis, taking eight reference fish species (Cynoglossus semilaevis, Scophthalmus maximus, Sparus aurata, Gasterosteus aculeatus, Xiphophorus maculatus, Oryzias latipes, Danio rerio, and Lepisosteus oculatus) for comparison, showed that the BACs of these two chromosomes of S. senegalensis were mainly distributed in two principal chromosomes in the reference species. Transposable Elements (TE) analysis showed significant differences between the two chromosomes, in terms of number of loci per Mb and coverage, and the class of TE (I or II) present. Analysis of TE divergence in chromosomes 2 and 4 compared to their syntenic regions in four reference fish species (C. semilaevis, S. maximus, O. latipes, and D. rerio) revealed differences in their age of activity compared with those species but less notable differences between the two chromosomes. Differences were also observed in peaks of divergence and coverage of TE families for all reference species even in those close to S. senegalensis, like S. maximus and C. semilaevis. Considered together, chromosomes 2 and 4 have evolved by Robertsonian fusions, pericentric inversions, and other chromosomal rearrangements mediated by TEs.

Novel Gene Rearrangement and the Complete Mitochondrial Genome of Cynoglossus monopus: Insights into the Envolution of the Family Cynoglossidae (Pleuronectiformes)

Cynoglossus monopus, a small benthic fish, belongs to the Cynoglossidae, Pleuronectiformes. It was rarely studied due to its low abundance and cryptical lifestyle. In order to understand the mitochondrial genome and the phylogeny in Cynoglossidae, the complete mitogenome of C. monopus has been sequenced and analyzed for the first time. The total length is 16,425 bp, typically containing 37 genes with novel gene rearrangements. The tRNA-Gln gene is inverted from the light to the heavy strand and translocated from the downstream of tRNA-Ile gene to its upstream. The control region (CR) translocated downstream to the 3'-end of ND1 gene adjoining to inverted to tRNA-Gln and left a 24 bp trace fragment in the original position. The phylogenetic trees were reconstructed by Bayesian inference (BI) and maximum likelihood (ML) methods based on the mitogenomic data of 32 tonguefish species and two outgroups. The results support the idea that Cynoglossidae is a monophyletic group and indicate that C. monopus has the closest phylogenetic relationship with C. puncticeps. By combining fossil records and mitogenome data, the time-calibrated evolutionary tree of families Cynoglossidae and Soleidae was firstly presented, and it was indicated that Cynoglossidae and Soleidae were differentiated from each other during Paleogene, and the evolutionary process of family Cynoglossidae covered the Quaternary, Neogene and Paleogene periods.

The distribution of the recessus orbitalis across flatfishes (order: Pleuronectiformes)

The recessus orbitalis is an accessory organ of flatfishes functioning in the protrusion of the eyes. This character, along with cranial asymmetry and a forward insertion of the dorsal fin, have been considered synapomorphies for the Pleuronectiformes. New dissections and examination of images taken in the wild show that the recessus orbitalis is present in all representatives of Pleuronectoidei examined but is absent in the single species of Psettoidei dissected. Psettoidei, the most primitive pleuronectiform lineage, contains three recognized species; thus, the absence of the recessus orbitalis in this whole lineage is unclear without further dissections. Ancestral character estimation at the family level for the recessus orbitalis indicates that the recessus orbitalis was likely absent in the common ancestor of Pleuronectiformes but was most likely present in the common ancestor of the Pleuronectoidei. Given that so few species of flatfishes have been assessed for the recessus orbitalis to date, additional characterization of the distribution of the recessus orbitalis across flatfishes will further inform what states this character may have and if it is a synapomorphy of Pleuronectiformes or simply a derived character state of Pleuronectoidei.

Shared ecological traits influence shape of the skeleton in flatfishes (Pleuronectiformes)

In the age of phylogenetic comparative methods, evolutionary biologists have been able to explore evolutionary trends in form in unique and extraordinarily diverse groups of animals. Pleuronectiformes, commonly known as flatfishes, is a diverse and specialized order of fishes that have remarkable asymmetry induced by ocular migration and a benthic life style. Although flatfishes are unique from other fishes, species within the group are morphologically diverse. The origin of ocular migration has been a primary focus of research; however, little is known about overall shape diversification among the flatfishes. In this study, we use integrative methods to examine how body shape evolved within the flatfishes. Shape was quantified from X-rays using geometric morphometrics for 389 individuals across 145 species. The most recent and robust phylogeny was overlaid onto the morphospace and phylogenetic signal was calculated to ascertain convergence in the morphospace. In addition, phylogenetic linear models were employed to determine if ecological traits were correlated with shape and if size had an effect on overall body shape. Results revealed that the majority of variation evolved recently, within the past 15–10-million-years in the middle Miocene, and is highly variable within the flatfishes. These changes are best summarized by body depth, jaw length and medial fin length. Dorsal and anal fin length are correlated, which may be due to the unique mode of locomotion used by flatfishes. A phylogenetic linear model and phylomorphospace analysis suggested that several ecological traits are correlated with shape, which indicates an ecological role in the diversification of flatfishes.

Evidence for a Robertsonian fusion in Solea senegalensis (Kaup, 1858) revealed by zoo-FISH and comparative genome analysis

Background

Solea senegalensis (Kaup, 1858) is a commercially important flatfish species, belonging to the Pleuronectiformes order. The taxonomy of this group has long been controversial, and the karyotype of the order presents a high degree of variability in diploid number, derived from chromosomal rearrangements such as Robertsonian fusions. Previously it has been proposed that the large metacentric chromosome of S. senegalensis arises from this kind of chromosome rearrangement and that this is a proto-sex chromosome.

Results

In this work, the Robertsonian origin of the large metacentric chromosome of S. senegalensis has been tested by the Zoo-FISH technique applied to two species of the Soleidae family (Dicologlossa cuneata and Dagetichthys lusitanica), and by comparative genome analysis with Cynoglossus semilaevis. From the karyotypic analysis we were able to determine a chromosome complement comprising 2n = 50 (FN = 54) in D. cuneata and 2n = 42 (FN = 50) in D. lusitanica. The large metacentric painting probe gave consistent signals in four acrocentric chromosomes of the two Soleidae species; and the genome analysis proved a common origin with four chromosome pairs of C. semilaevis. As a result of the genomic analysis, up to 61 genes were annotated within the thirteen Bacterial Artificial Chromosome clones analysed.

Conclusions

These results confirm that the large metacentric chromosome of S. senegalensis originated from a Robertsonian fusion and provide new data about the chromosome evolution of S. senegalensis in particular, and of Pleuronectiformes in general.

Electronic supplementary material

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

Highly dense linkage maps from 31 full-sibling families of turbot (Scophthalmus maximus) provide insights into recombination patterns and chromosome rearrangements throughout a newly refined genome assembly

Highly dense linkage maps enable positioning thousands of landmarks useful for anchoring the whole genome and for analysing genome properties. Turbot is the most important cultured flatfish worldwide and breeding programs in the fifth generation of selection are targeted to improve growth rate, obtain disease resistant broodstock and understand sex determination to control sex ratio. Using a Restriction-site Associated DNA approach, we genotyped 18,214 single nucleotide polymorphism in 1,268 turbot individuals from 31 full-sibling families. Individual linkage maps were combined to obtain a male, female and species consensus maps. The turbot consensus map contained 11,845 markers distributed across 22 linkage groups representing a total normalised length of 3,753.9 cM. The turbot genome was anchored to this map, and scaffolds representing 96% of the assembly were ordered and oriented to obtain the expected 22 megascaffolds according to its karyotype. Recombination rate was lower in males, especially around centromeres, and pairwise comparison of 44 individual maps suggested chromosome polymorphism at specific genomic regions. Genome comparison across flatfish provided new evidence on karyotype reorganisations occurring across the evolution of this fish group.

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.

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.

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.

Identifying Optimal Models of Evolution

Most phylogenetic methods are model-based and depend on models of evolution designed to approximate the evolutionary processes. Several methods have been developed to identify suitable models of evolution for phylogenetic analysis of alignments of nucleotide or amino acid sequences and some of these methods are now firmly embedded in the phylogenetic protocol. However, in a disturbingly large number of cases, it appears that these models were used without acknowledgement of their inherent shortcomings. In this chapter, we discuss the problem of model selection and show how some of the inherent shortcomings may be identified and overcome.

Phylogenomic analysis of carangimorph fishes reveals flatfish asymmetry arose in a blink of the evolutionary eye

Background

Flatfish cranial asymmetry represents one of the most remarkable morphological innovations among vertebrates, and has fueled vigorous debate on the manner and rate at which strikingly divergent phenotypes evolve. A surprising result of many recent molecular phylogenetic studies is the lack of support for flatfish monophyly, where increasingly larger DNA datasets of up to 23 loci have either yielded a weakly supported flatfish clade or indicated the group is polyphyletic. Lack of resolution for flatfish relationships has been attributed to analytical limitations for dealing with processes such as nucleotide non-stationarity and incomplete lineage sorting (ILS). We tackle this phylogenetic problem using a sequence dataset comprising more than 1,000 ultraconserved DNA element (UCE) loci covering 45 carangimorphs, the broader clade containing flatfishes and several other specialized lineages such as remoras, billfishes, and archerfishes.

Results

We present a phylogeny based on UCE loci that unequivocally supports flatfish monophyly and a single origin of asymmetry. We document similar levels of discordance among UCE loci as in previous, smaller molecular datasets. However, relationships among flatfishes and carangimorphs recovered from multilocus concatenated and species tree analyses of our data are robust to the analytical framework applied and size of data matrix used. By integrating the UCE data with a rich fossil record, we find that the most distinctive carangimorph bodyplans arose rapidly during the Paleogene (66.0–23.03 Ma). Flatfish asymmetry, for example, likely evolved over an interval of no more than 2.97 million years.

Conclusions

The longstanding uncertainty in phylogenetic hypotheses for flatfishes and their carangimorph relatives highlights the limitations of smaller molecular datasets when applied to successive, rapid divergences. Here, we recovered significant support for flatfish monophyly and relationships among carangimorphs through analysis of over 1,000 UCE loci. The resulting time-calibrated phylogeny points to phenotypic divergence early within carangimorph history that broadly matches with the predictions of adaptive models of lineage diversification.

Electronic supplementary material

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

The complete mitochondrial genome of Pseudorhombus dupliocellatus (Pleuronectiformes: Paralichthyidae)

The Pseudorhombus dupliocellatus belongs to family Paralichthyidae of Pleuronectiformes. In this study, the complete mitochondrial genome of P. dupliocellatus is determined and described. The mitogenome is 16 621 bp in length and consists of 13 protein-coding genes, 22 tRNAs, 2 rRNAs, a control region, and a L-strand replication origin. The arrangement of the mitogenome is identical to that of the typical teleost. The overall base composition is 26.9%, 25.3%, 31.0%, and 16.8% for A, T, C, and G, respectively, with a slight bias on A+T content (52.2%). The phylogenetic tree of 13 species all in Pleuronectiformes demonstrated that P. dupliocellatus, as well as the other Paralichthyidae fishes containing Paralichthys olivaceus and Pseudorhombus cinnamoneus, clustered in a clade and had a closer relationship with Pleuronectidae species than Bothidae ones. This study is expected to contributing to the systematic evolution of Paralichthyidae and further Pleuronectiformes.

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.

Concerted Evolution of Duplicate Control Regions in the Mitochondria of Species of the Flatfish Family Bothidae (Teleostei: Pleuronectiformes)

Mitogenomes of flatfishes (Pleuronectiformes) exhibit the greatest diversity of gene rear-rangements in teleostean fishes. Duplicate control regions (CRs) have been found in the mito-genomes of two flatfishes, Samariscus latus (Samaridae) and Laeops lanceolata (Bothidae), which is rare in teleosts. It has been reported that duplicate CRs have evolved in a concerted fashion in fishes and other animals, however, whether concerted evo-lution exists in flatfishes remains unknown. In this study, based on five newly sequenced and six previously reported mitogenomes of lefteye flounders in the Bothidae, we explored whether duplicate CRs and concerted evolution exist in these species. Results based on the present study and previous reports show that four out of eleven bothid species examined have duplicate CRs of their mitogenomes. The core regions of the duplicate CRs of mitogenomes in the same species have identical, or nearly identical, sequences when compared to each other. This pattern fits the typical characteristics of concerted evolution. Additionally, phylogenetic and ancestral state reconstruction analysis also provided evidence to support the hypothesis that duplicate CRs evolved concertedly. The core region of concerted evolution is situated at the conserved domains of the CR of the mitogenome from the termination associated sequences (TASs) to the conserved sequence blocks (CSBs). Commonly, this region is con-sidered to regulate mitochondrial replication and transcription. Thus, we hypothesize that the cause of concerted evolution of the duplicate CRs in the mtDNAs of these four bothids may be related to some function of the conserved sequences of the CRs during mitochondrial rep-lication and transcription. We hope our results will provide fresh insight into the molecular mechanisms related to replication and evolution of mitogenomes.

Molecular phylogenetics of sculpins of the subfamily Oligocottinae (Cottidae)

The sculpin subfamily Oligocottinae includes 18-20 species of nearshore benthic fishes with a diverse array of reproductive strategies. As a first step toward understanding the evolution of that diversity, we conducted a phylogenetic study based on DNA sequences from eight genomic regions from 31 sculpin species aimed at testing monophyly and relationships of the Oligocottinae. Representatives from the perciform families Agonidae, Cottidae, Hemitripteridae, Hexagrammidae, Psychrolutidae, and Rhamphocottidae served as outgroups. The sequence data were analyzed in maximum likelihood and Bayesian phylogenetic inference frameworks. Results of these analyses show that a systematic revision of the group is warranted. The genus Clinocottus is a polyphyletic assemblage of three distinct lineages, which should be indicated by resurrection of the subgenera Blennicottus, Clinocottus, and Oxycottus; Leiocottus hirundo is more closely related to Clinocottus analis than C. analis is related to any other member of Clinocottus; the composition of the tribe Oligocottini should be revised to include only the genera Oligocottus, Clinocottus, and Orthonopias; and the genus Sigmistes should be removed from the subfamily Oligocottinae.

Mitochondrial genomic investigation of flatfish monophyly

We present the first study to use whole mitochondrial genome sequences to examine phylogenetic affinities of the flatfishes (Pleuronectiformes). Flatfishes have attracted attention in evolutionary biology since the early history of the field because understanding the evolutionary history and patterns of diversification of the group will shed light on the evolution of novel body plans. Because recent molecular studies based primarily on DNA sequences from nuclear loci have yielded conflicting results, it is important to examine phylogenetic signal in different genomes and genome regions. We aligned and analyzed mitochondrial genome sequences from thirty-nine pleuronectiforms including nine that are newly reported here, and sixty-six non-pleuronectiforms (twenty additional clade L taxa [Carangimorpha or Carangimorpharia] and forty-six secondary outgroup taxa). The analyses yield strong support for clade L and weak support for the monophyly of Pleuronectiformes. The suborder Pleuronectoidei receives moderate support, and as with other molecular studies the putatively basal lineage of Pleuronectiformes, the Psettodoidei is frequently not most closely related to other pleuronectiforms. Within the Pleuronectoidei, the basal lineages in the group are poorly resolved, however several flatfish subclades receive consistent support. The affinities of Lepidoblepharon and Citharoides among pleuronectoids are particularly uncertain with these data.

Fine mapping and evolution of the major sex determining region in turbot (Scophthalmus maximus)

Fish sex determination (SD) systems are varied, suggesting evolutionary changes including either multiple evolution origins of genetic SD from nongenetic systems (such as environmental SD) and/or turnover events replacing one genetic system by another. When genetic SD is found, cytological differentiation between the two members of the sex chromosome pair is often minor or undetectable. The turbot (Scophthalmus maximus), a valuable commercial flatfish, has a ZZ/ZW system and a major SD region on linkage group 5 (LG5), but there are also other minor genetic and environmental influences. We here report refined mapping of the turbot SD region, supported by comparative mapping with model fish species, to identify the turbot master SD gene. Six genes were located to the SD region, two of them associated with gonad development (sox2 and dnajc19). All showed a high association with sex within families (P = 0), but not at the population level, so they are probably partially sex-linked genes, but not SD gene itself. Analysis of crossovers in LG5 using two families confirmed a ZZ/ZW system in turbot and suggested a revised map position for the master gene. Genetic diversity and differentiation for 25 LG5 genetic markers showed no differences between males and females sampled from a wild population, suggesting a recent origin of the SD region in turbot. We also analyzed associations with markers of the most relevant sex-related linkage groups in brill (S. rhombus), a closely related species to turbot; the data suggest that an ancient XX/XY system in brill changed to a ZZ/ZW mechanism in turbot.

Consolidation of the genetic and cytogenetic maps of turbot (Scophthalmus maximus) using FISH with BAC clones

Bacterial artificial chromosomes (BAC) have been widely used for fluorescence in situ hybridization (FISH) mapping of chromosome landmarks in different organisms, including a few in teleosts. In this study, we used BAC-FISH to consolidate the previous genetic and cytogenetic maps of the turbot (Scophthalmus maximus), a commercially important pleuronectiform. The maps consisted of 24 linkage groups (LGs) but only 22 chromosomes. All turbot LGs were assigned to specific chromosomes using BAC probes obtained from a turbot 5× genomic BAC library. It consisted of 46,080 clones with inserts of at least 100 kb and <5 % empty vectors. These BAC probes contained gene-derived or anonymous markers, most of them linked to quantitative trait loci (QTL) related to productive traits. BAC clones were mapped by FISH to unique marker-specific chromosomal positions, which showed a notable concordance with previous genetic mapping data. The two metacentric pairs were cytogenetically assigned to LG2 and LG16, and the nucleolar organizer region (NOR)-bearing pair was assigned to LG15. Double-color FISH assays enabled the consolidation of the turbot genetic map into 22 linkage groups by merging LG8 with LG18 and LG21 with LG24. In this work, a first-generation probe panel of BAC clones anchored to the turbot linkage and cytogenetical map was developed. It is a useful tool for chromosome traceability in turbot, but also relevant in the context of pleuronectiform karyotypes, which often show small hardly identifiable chromosomes. This panel will also be valuable for further integrative genomics of turbot within Pleuronectiformes and teleosts, especially for fine QTL mapping for aquaculture traits, comparative genomics, and whole-genome assembly.