Phylogeny of the Elopomorpha (Teleostei): Evidence from six nuclear and mitochondrial markers

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.

The First Genome-Wide Survey of Shortbelly Eel (Dysomma anguillare Barnard, 1923) to Provide Genomic Characteristics, Microsatellite Markers and Complete Mitogenome Information

Dysomma anguillare is a demersal eel widespread distributing in tropical waters of the Indo-West Pacific and Atlantic. As an important component of the coastal fishery and marine ecosystem, the lack of genomic information for this species severely restricts the progress of relevant researches. In this study, the abecedarian genome-wide characteristics and phylogenetic relationships analyses were carried out based on next-generation sequencing (NGS) platform. The revised genome size was approximately 1 310 Mb, with the largest scaffold length reaching 23 878 bp through K-mer (K = 17) method. The heterozygosity, repetitive rate and average GC content were about 0.94%, 51.93% and 42.23%, respectively. A total of 1 160 104 microsatellite motifs were identified from the de novo assembled genome of D. anguillare, in which dinucleotide repeats accounted for the largest proportion (592 234, 51.05%), the highest occurrence frequency (14.58%) as well as the largest relative abundance (379.27/Mb). The high-polymorphic and moderate-polymorphic loci composed around 73% of the total single sequence repeats (SSRs), showing a latent capacity for subsequent population genetic structure and genetic diversity appraisal researches. Another byproduct of whole-genome sequencing, the double-stranded and circular mitogenome (16 690 bp) was assembled to investigate the evolutionary relationships of D. anguillare. The phylogenic tree constructed with maximum likelihood (ML) method showed that D. anguillare was closely related to Synaphobranchidae species, and the molecular systematic results further supported classical taxonomy status of D. anguillare.

Genomic microsatellite characteristics analysis of Dysommaanguillare (Anguilliformes, Dysommidae), based on high-throughput sequencing technology

Microsatellite loci were screened from the genomic data of Dysommaanguillare and their composition and distribution were analysed by bioinformatics for the first time. The results showed that 4,060,742 scaffolds with a total length of 1,562 Mb were obtained by high-throughput sequencing and 1,160,104 microsatellite loci were obtained by MISA screening, which were distributed on 770,294 scaffolds. The occurrence frequency and relative abundance were 28.57% and 743/Mb, respectively. Amongst the six complete microsatellite types, dinucleotide repeats accounted for the largest proportion (592,234, 51.05%), the highest occurrence frequency (14.58%) and the largest relative abundance (379.27/Mb). A total of 1488 microsatellite repeats were detected in the genome of D.anguillare, amongst which the hexanucleotide repeat motifs were the most abundant (608), followed by pentanucleotide repeat motifs (574), tetranucleotide repeat motifs (232), trinucleotide repeat motifs (59), dinucleotide repeat motifs (11) and mononucleotide repeat motifs (4). The abundance of microsatellites of the same repeat type decreased with the increase of copy numbers. Amongst the six types of nucleotide repeats, the preponderance of repeated motifs are A (191,390, 43.77%), CA (150,240, 25.37%), AAT (13,168, 14.05%), CACG (2,649, 8.14%), TAATG (119, 19.16%) and CCCTAA (190, 19.16%, 7.65%), respectively. The data of the number, distribution and abundance of different types of microsatellites in the genome of D.anguillare were obtained in this study, which would lay a foundation for the development of high-quality microsatellite markers of D.anguillare in the future.

Facial and opercular muscles in the Anguilliformes (Elopomorpha: Teleostei): Comparative anatomy and phylogenetic implications for the basal position of Protanguilla

The teleost order Anguilliformes, true eels, comprises more than 1000 described species in 20 families, commonly known as eels, congers, morays, and gulper eels. Comprehensive studies of Anguilliformes are limited, resulting in a lack of consensus for morphology-based phylogenetic hypotheses. A detailed morphological analysis of the cephalic and opercular myology offers a promising new source of characters to help elucidate the intrarelationships of Anguilliformes. Our study is the most extensive myological analysis for the group and includes 97 terminal taxa, with representatives from each of the 20 families of Anguilliformes plus outgroup clades. Results demonstrate that muscle characters inform phylogenetic relationships within Anguilliformes, and we propose two new synapomorphies for all extant members, including Protanguilla palau, the "living fossil"-adductor mandibulae originating on the parietal (vs. restricted to suspensorium) and segmentum mandibularis absent (vs. present). Exceptions for the first condition characterize highly modified saccopharyngoids, and for the second one, Notacanthidae. More importantly, we suggest three new synapomorphies for the remaining extant anguilliforms (except in highly modified saccopharyngoids)-adductor mandibulae originates on the frontals (vs. frontals naked), adductor mandibulae stegalis is separated from the rictalis (vs. ricto-stegalis fused into a single piece), and the levator operculi inserts on the lateral surface of the opercle (vs. medial surface of the opercle). Our phylogenetic optimization strongly corroborates the hypothesis that Protanguilla is the sister group of all other extant eels. A further goal of this paper is to clearly document the substantive conflicts between the available molecular data and the extensive and diverse morphological evidence.

Comparative Mitogenomic Analysis of Two Snake Eels Reveals Irregular Gene Rearrangement and Phylogenetic Implications of Ophichthidae

The family Ophichthidae has the largest number and the most various species (about 359 valid species) in the order Anguilliformes worldwide. Both morphological and molecular characteristics have been used to assess their taxonomic status. However, due to the ambiguous morphological features, molecular data such as mitochondrial DNA sequences have been implemented for the correct identification and classification of these fishes. In this study, the gene arrangement and structure characteristics of two Ophichthidae mitochondrial genomes were investigated for the first time. The total mitogenome lengths of O. evermanni and O. erabo were 17,759 bp and 17,856 bp, respectively. Comparing with the ancestral mitochondrial gene order, the irregular gene rearrangement happened between ND6 and tRNA-Pro (P) genes with another similar control region emerging between tRNA-Thr (T) and ND6 genes, which could be explained by the tandem duplication and random loss (TDRL) model appropriately. ML phylogenetic tree demonstrated that the family Ophichthidae was monophyletic origin, but genus Ophichthus might be polyphyletic because of the confused cluster relationships among different species.

Two New Moray Eels of Genera Diaphenchelys and Gymnothorax (Anguilliformes: Muraenidae) from Taiwan and the Philippines

Two new moray eels of the genera Diaphenchelys and Gymnothorax from Taiwan and the Philippines are described. Diaphenchelys laimospila sp. nov. is described based on two specimens that represent the third species and a new geographic record of the genus. It can be distinguished from the other two congeners by the number of cephalic sensory pores, vertebral formula, morphometric measurements, and the coloration pattern. Gymnothorax pseudokidako sp. nov. is a muraenid with a dark brown body covered by pale snowflake-like blotches. It differs from the most similar species Gymnothorax kidako (Temminck and Schlegel) by having a relatively short tail (50.5-53.0% vs. 52.9-56.4% of TL), more dentary teeth (17-26 vs. 16-20), fewer total vertebrae (134-139 vs. 137-143), and the absence of white margin on anal fin (vs. prominent white margin). These two new species were also confirmed by molecular analyses, the mitochondrial COI gene (593 bp) for D. laimospila, and the nuclear EGR3 gene (767 bp) for G. pseudokidako.

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.

Resolving the Early Divergence Pattern of Teleost Fish Using Genome-Scale Data

Regarding the phylogenetic relationship of the three primary groups of teleost fishes, Osteoglossomorpha (bonytongues and others), Elopomorpha (eels and relatives), Clupeocephala (the remaining teleost fish), early morphological studies hypothesized the first divergence of Osteoglossomorpha, whereas the recent prevailing view is the first divergence of Elopomorpha. Molecular studies supported all the possible relationships of the three primary groups. This study analyzed genome-scale data from four previous studies: 1) 412 genes from 12 species, 2) 772 genes from 15 species, 3) 1,062 genes from 30 species, and 4) 491 UCE loci from 27 species. The effects of the species, loci, and models used on the constructed tree topologies were investigated. In the analyses of the data sets (1)–(3), although the first divergence of Clupeocephala that left the other two groups in a sister relationship was supported by concatenated sequences and gene trees of all the species and genes, the first divergence of Elopomorpha among the three groups was supported using species and/or genes with low divergence of sequence and amino-acid frequencies. This result corresponded to that of the UCE data set (4), whose sequence divergence was low, which supported the first divergence of Elopomorpha with high statistical significance. The increase in accuracy of the phylogenetic construction by using species and genes with low sequence divergence was predicted by a phylogenetic informativeness approach and confirmed by computer simulation. These results supported that Elopomorpha was the first basal group of teleost fish to have diverged, consistent with the prevailing view of recent morphological studies.

Novel gene rearrangement in the mitochondrial genome of Muraenesox cinereus and the phylogenetic relationship of Anguilliformes

The structure and gene sequence of the fish mitochondrial genome are generally considered to be conservative. However, two types of gene arrangements are found in the mitochondrial genome of Anguilliformes. In this paper, we report a complete mitogenome of Muraenesox cinereus (Anguilliformes: Muraenesocidae) with rearrangement phenomenon. The total length of the M. cinereus mitogenome was 17,673 bp, and it contained 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNA genes, and two identical control regions (CRs). The mitochondrial genome of M. cinereus was obviously rearranged compared with the mitochondria of typical vertebrates. The genes ND6 and the conjoint trnE were translocated to the location between trnT and trnP, and one of the duplicated CR was translocated to the upstream of the ND6. The tandem duplication and random loss is most suitable for explaining this mitochondrial gene rearrangement. The Anguilliformes phylogenetic tree constructed based on the whole mitochondrial genome well supports Congridae non-monophyly. These results provide a basis for the future Anguilliformes mitochondrial gene arrangement characteristics and further phylogenetic research.

Molecular cytogenetics insights in two pelagic big-game fishes in the Atlantic, the tarpon, Megalops atlanticus (Elopiformes: Megalopidae), and the sailfish, Istiophorus platypterus (Istiophoriformes: Istiophoridae)

Abstract Some pelagic and usually large sized fishes are preferential targets for sport and commercial fishing. Despite their economic importance, cytogenetic data on their evolutionary processes and management are very deficient, especially due to logistical difficulties. Here, information for two of such charismatic species, the tarpon, Megalops atlanticus (Elopiformes: Megalopidae), and the sailfish, Istiophorus platypterus (Istiophoriformes: Istiophoridae), both with a wide Atlantic distribution, were provided. Cytogenetic data were obtained using conventional methods (Giemsa staining, Ag-NORs technique, and C-banding), base-specific fluorochrome staining and fluorescence in situ hybridization (FISH) with rDNA probes. Megalops atlanticus has 2n = 50 chromosomes, all acrocentric ones (NF = 50), while Istiophorus platypterus has 2n = 48 chromosomes, 2m + 2st + 44a (NF = 52). Megalops atlanticus populations from the South Atlantic and Caribbean share identical karyotypic patterns, likely associated with gene flow between them. In turn, I. platypterus presents karyotype similarities with phylogenetically close groups, such as Carangidae. The chromosomal characteristics of these species highlight their independent evolutionary paths. Additionally, the current data contribute to knowledge of new aspects of pelagic fish fauna and will support further comparative studies with congeneric species, clarifying evolutionary karyotype trends of these fish groups.

Diversity, phylogeny, and historical biogeography of large-eye seabreams (Teleostei: Lethrinidae)

The large-eye seabreams or Monotaxinae is one of two subfamilies in the Lethrinidae, a family of perch-like coral reef fishes. Despite its widespread occurrence and its commercial interest in the tropical Indo-West Pacific (IWP), this subfamily has traditionally been considered a taxonomically difficult group. Based on 268 samples collected from all 15 known large-eye seabream species throughout their distribution ranges, we investigated the taxonomic diversity and phylogenetic relationships in the subfamily. From the results of multiple analyses on four gene markers, we confirmed the monophyly of all four genera in the subfamily (Gnathodentex, Gymnocranius, Monotaxis and Wattsia). We confirmed the occurrence of two species in the genus Monotaxis. We reported 15 delimited species within the most speciose genus Gymnocranius, four of which are potentially new species. The time-calibrated phylogenetic reconstruction enabled us to clarify the evolutionary history of the large-eye seabreams and to infer past patterns of species distribution. The most recent common ancestor to the Monotaxinae likely occurred in the central IWP ca. 32 million years ago. A burst of species diversification likely took place during the Mid- to Late Miocene, coinciding with tectonic change in the central IWP region. This gave rise to most extant lineages in Gymnocranius. The observed geographic distribution patterns in the subfamily most likely point to the central IWP as the area of origin and diversification. This was followed by multiple events of centrifugal range expansion towards either the Indian Ocean or the western Pacific Ocean, or both. Our results thus provide new support for S. Ekman's center-of-origin hypothesis.

Gonadotropin-inhibitory hormone in teleosts: New insights from a basal representative, the eel

Since its discovery in birds, gonadotropin-inhibitory hormone (GnIH) has triggered investigation in the other groups of vertebrates. In the present study, we have identified a single gnih gene in the European eel (Anguilla anguilla), a representative species of a basal group of teleosts (Elopomorphs). We have also retrieved a single gnih gene in Osteoglossomorphs, as well as in more recently emerged teleosts, Clupeocephala. Phylogeny and synteny analyses allowed us to infer that one of the two gnih paralogs emerged from the teleost-specific whole genome duplication (TWGD or 3R), would have been lost shortly after the 3R, before the emergence of the basal groups of teleosts. This led to the presence of a single gnih in extant teleosts as in other vertebrates. Two gnih paralogs were still found in some teleost species, such as in salmonids, but resulting from the additional whole genome duplication that specifically occurred in this lineage (4R). Eel gnih was mostly expressed in the diencephalon part of the brain, as analyzed by quantitative real-time PCR. Cloning of eel gnih cDNA confirmed that the sequence of the GnIH precursor encoded three putative mature GnIH peptides (aaGnIH-1, aaGnIH-2 and aaGnIH-3), which were synthesized and tested for their direct effects on eel pituitary cells in vitro. Eel GnIH peptides inhibited the expression of gonadotropin subunits (lhβ, fshβ, and common a-subunit) as well as of GnRH receptor (gnrh-r2), with no effect on tshβ and gh expression. The inhibitory effect of GnIH peptides on gonadotropic function in a basal teleost is in agreement with an ancestral inhibitory role of GnIH in the neuroendocrine control of reproduction in vertebrates.

Identification and stable expression of vitellogenin receptor through vitellogenesis in the European eel

In teleosts, vitellogenin (Vtg) is a phospholipoglycoprotein synthesized by the liver, released into the blood circulation and incorporated into the oocytes via endocytosis mediated by the Vtg receptor (VTGR) to form the yolk granules. The VTGR is crucial for oocyte growth in egg-laying animals but is also present in non-oviparous vertebrates, such as human. The VTGR belongs to the low-density lipoprotein receptor superfamily (LDLR) and is also named very-low-density lipoprotein receptor (VLDLR). In this study, we identified and phylogenetically positioned the VTGR of a basal teleost, the European eel, Anguilla anguilla. We developed quantitative real-time PCR (qRT-PCR) and investigated the tissue distribution of vtgr transcripts. We compared by qRT-PCR the ovarian expression levels of vtgr in juvenile yellow eels and pre-pubertal silver eels. We also analyzed the regulation of ovarian vtgr expression throughout vitellogenesis in experimentally matured eels. The Vtg plasma level was measured by homologous ELISA experimental maturation. Our in silico search and phylogenetical analysis revealed a single vtgr in the European eel, orthologous to other vertebrate vtgr. The qRT-PCR studies revealed that vtgr is mainly expressed in the ovary and also detected in various other tissues such as brain, pituitary, gill, fat, heart, and testis, suggesting some extra-ovarian functions of VTGR. We showed that vtgr is expressed in ovaries of juvenile yellow eels with no higher expression in pre-pubertal silver eels nor in experimentally matured eels. This suggests that vtgr transcription already occurs during early pre-vitellogenesis of immature eels and is not further activated in vitellogenic oocytes. European eel Vtg plasma level increased throughout experimental maturation in agreement with previous studies. Taken together, these results suggest that vtgr transcript levels may not be a limiting step for the uptake of Vtg by the oocyte in the European eel.

Complete mitochondrial genome of Ophichthus brevicaudatus reveals novel gene order and phylogenetic relationships of Anguilliformes

Generally, a teleostean group possesses only one type or a set of similar mitochondrial gene arrangement. However, two types of gene arrangement have been identified in the mitochondrial genomes (mitogenomes) of Anguilliformes. Here, a newly sequenced mitogenome of Ophichthus brevicaudatus (Anguilliformes; Ophichthidae) was presented. The total length of the O. brevicaudatus mitogenome was 17,773 bp, and it contained 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNA (tRNA) genes, and two identical control regions (CRs). The gene order differed from that of the typical vertebrate mitogenomes. The genes ND6 and the conjoint trnE were translocated to the location between trnT and trnP, and one of the duplicated CR was translocated to the upstream of the ND6. The duplication-random loss model was adopted to explain the gene rearrangement events in this mitogenome. The most comprehensive phylogenetic trees of Anguilliformes based on complete mitogenome was constructed. The non-monophyly of Congridae was well supported, whereas the non-monophyly of Derichthyidae and Chlopsidae was not supported. These results provide insight into gene arrangement features of anguilliform mitogenomes and lay the foundation for further phylogenetic studies on Anguilliformes.

Phylogenetic interrelationships of the eel families Derichthyidae and Colocongridae (Elopomorpha: Anguilliformes) based on the pectoral skeleton

A cladistic analysis of the eel families Derichthyidae and Colocongridae is herein proposed for the first time on the basis of morphological data. We discovered dozens of new phylogenetic characters derived from a detailed analysis of the pectoral skeleton, an anatomical system neglected by most previous studies. Our maximum parsimony analysis indicates that Colocongridae sensu lato is paraphyletic, with its two constituent genera Coloconger and Congriscus appearing as successive sister groups of derichthyids. Monophyly of the family Derichthyidae, which has been questioned by some studies, is herein strongly supported by 10 unambiguous synapomorphies. We also stress the importance of the appendicular skeleton as a useful source of phylogenetic information for the resolution of systematic problems within Anguilliformes.

Molecular ontogeny of larval immunity in European eel at increasing temperatures

Temperature is a major factor that modulates the development and reactivity of the immune system. Only limited knowledge exists regarding the immune system of the catadromous European eel, Anguilla anguilla, especially during the oceanic early life history stages. Thus, a new molecular toolbox was developed, involving tissue specific characterisation of 3 housekeeping genes, 9 genes from the innate and 3 genes from the adaptive immune system of this species. The spatial pattern of immune genes reflected their function, e.g. complement component c3 was mainly produced in liver and il10 in the head kidney. Subsequently, the ontogeny of the immune system was studied in larvae reared from hatch to first-feeding at four temperatures, spanning their thermal tolerance range (16, 18, 20, and 22 °C). Expression of some genes (c3 and igm) declined post hatch, whilst expression of most other genes (mhc2, tlr2, il1β, irf3, irf7) increased with larval age. At the optimal temperature, 18 °C, this pattern of immune-gene expression revealed an immunocompromised phase between hatch (0 dph) and teeth-development (8 dph). The expression of two of the studied genes (mhc2, lysc) was temperature dependent, leading to increased mRNA levels at 22 °C. Additionally, at the lower end of the thermal spectrum (16 °C) immune competency appeared reduced, whilst close to the upper thermal limit (22 °C) larvae showed signs of thermal stress. Thus, protection against pathogens is probably impaired at temperatures close to the critical thermal maximum (CT_max), impacting survival and productivity in hatcheries and natural recruitment.

Rhodopsin gene evolution in early teleost fishes

Rhodopsin mediates an essential step in image capture and is tightly associated with visual adaptations of aquatic organisms, especially species that live in dim light environments (e.g., the deep sea). The rh1 gene encoding rhodopsin was formerly considered a single-copy gene in genomes of vertebrates, but increasing exceptional cases have been found in teleost fish species. The main objective of this study was to determine to what extent the visual adaptation of teleosts might have been shaped by the duplication and loss of rh1 genes. For that purpose, homologous rh1/rh1-like sequences in genomes of ray-finned fishes from a wide taxonomic range were explored using a PCR-based method, data mining of public genetic/genomic databases, and subsequent phylogenomic analyses of the retrieved sequences. We show that a second copy of the fish-specific intron-less rh1 is present in the genomes of most anguillids (Elopomorpha), Hiodon alosoides (Osteoglossomorpha), and several clupeocephalan lineages. The phylogenetic analysis and comparisons of alternative scenarios for putative events of gene duplication and loss suggested that fish rh1 was likely duplicated twice during the early evolutionary history of teleosts, with one event coinciding with the hypothesized fish-specific genome duplication and the other in the common ancestor of the Clupeocephala. After these gene duplication events, duplicated genes were maintained in several teleost lineages, whereas some were secondarily lost in specific lineages. Alternative evolutionary schemes of rh1 and comparison with previous studies of gene evolution are also reviewed.

Studies of Myxidium giardi Cépède, 1906 infections in Icelandic eels identifies a genetically diverse clade of myxosporeans that represents the Paramyxidium n. g. (Myxosporea: Myxidiidae)

BACKGROUND

The myxosporean Myxidium giardi Cépède, 1906 was described infecting the kidney of the European eel, Anguilla anguilla (L.), having spindle-shaped myxospores and terminal sub-spherical polar capsules. Since then, numerous anguillid eels globally have been documented to have similar Myxidium infections. Many of these have been identified using the morphological features of myxospores or by the location of infection in the host, and some have been subsequently synonymised with M. giardi. Therefore, it is not clear whether M. giardi is a widely distributed parasite, infecting numerous species of eels, in multiple organs, or whether some infections represent other, morphologically similar but different species of myxosporeans. The aim of the present study was to assess the status of M. giardi infections in Icelandic eels, and related fish hosts in Malaysia and to use spore morphology and molecular techniques to evaluate the diversity of myxosporeans present.

RESULTS

The morphologies of the myxospores from Icelandic eels were very similar but the overall dimensions were significantly different from the various tissue locations. Myxospores from the kidney of the Malaysian tarpon, Megalops cyprinoides (Broussonet), were noticeably smaller. However, the SSU rDNA sequences from the different tissues locations in eels, were all very distinct, with percentage similarities ranging from 92.93% to as low as 89.8%, with the sequence from Malaysia being even more dissimilar. Molecular phylogenies consistently placed these sequences together in a clade that we refer to as the Paramyxidium clade that is strongly associated with the Myxidium clade (sensu stricto). We erect the genus Paramyxidium n. g. (Myxidiidae) to accommodate these histozoic taxa, and transfer Myxidium giardi as Paramyxidium giardi Cépède, 1906 n. comb. as the type-species.

CONCLUSIONS

There is not a single species of Myxidium (M. giardi) causing systemic infections in eels in Iceland. There are three species, confirmed with a robust phylogeny, one of which represents Paramyxidium giardi n. comb. Additional species probably exist that infect different tissues in the eel and the site of infection in the host fish is an important diagnostic feature for this group (Paramyxidium n. g. clade). Myxospore morphology is generally conserved in the Paramyxidium clade, although actual spore dimensions can vary between some species. Paramyxidium spp. are currently only known to infect fishes from the Elopomorpha.

A review of the systematic biology of fossil and living bony-tongue fishes, Osteoglossomorpha (Actinopterygii: Teleostei)

ABSTRACT The bony-tongue fishes, Osteoglossomorpha, have been the focus of a great deal of morphological, systematic, and evolutionary study, due in part to their basal position among extant teleostean fishes. This group includes the mooneyes (Hiodontidae), knifefishes (Notopteridae), the abu (Gymnarchidae), elephantfishes (Mormyridae), arawanas and pirarucu (Osteoglossidae), and the African butterfly fish (Pantodontidae). This morphologically heterogeneous group also has a long and diverse fossil record, including taxa from all continents and both freshwater and marine deposits. The phylogenetic relationships among most extant osteoglossomorph families are widely agreed upon. However, there is still much to discover about the systematic biology of these fishes, particularly with regard to the phylogenetic affinities of several fossil taxa, within Mormyridae, and the position of Pantodon. In this paper we review the state of knowledge for osteoglossomorph fishes. We first provide an overview of the diversity of Osteoglossomorpha, and then discuss studies of the phylogeny of Osteoglossomorpha from both morphological and molecular perspectives, as well as biogeographic analyses of the group. Finally, we offer our perspectives on future needs for research on the systematic biology of Osteoglossomorpha.

Reconsidering pectoral girdle and fin morphology in Anguillidae (Elopomorpha: Anguilliformes)

The morphology of the pectoral girdle and fin of Anguillidae is reconsidered via the inclusion of skeletal components that have previously been unassessed. For example, the pectoral girdle and fin in Anguilla were erroneously reported to lack a scapular bone, a cartilaginous scapulocoracoid plate and a cartilaginous propterygium. The pectoral morphology of Anguilla is also compared with the closely related genus Nemichthys, including additional data on the anatomy of this eel family.

Phylogenetic analysis shows the general diversification pattern of deep-sea notacanthiforms (Teleostei: Elopomorpha)

The Notacanthiformes is an ancient group of deep-sea ray-finned fishes comprising 27 species in two families; Halosauridae and Notacanthidae. Although many studies have tried to reconstruct the phylogenetic relationships among the major clades of Elopomorpha, little is known about the evolutionary history of notacanthiforms. Molecular and morphological data were used to test previous hypotheses regarding the phylogenetic relationships among notacanthiform taxa, and to unravel the origin and evolution of this group. The molecular analyses of notacanthids showed similar results to those previously obtained employing osteological data, which proposed the existence of the Lipogenyinae (Lipogenys) and Notacanthinae (Notacanthus + Polyacanthonotus) subfamilies. Nevertheless, when the external morphology data is considered Lipogenys is more related to Notacanthus than Polyacanthonotus. The analyses could not fully resolve the inner relationships of the halosaurids. The time-calibrated tree of the order Notacanthiformes shows a long process of diversification spanning from the upper Cretaceous, to 50 million years after the K-Pg extinction, with the gradual emergence of all the modern families and genera of the group. This is the first specific phylogeny of the order Notacanthiformes, combining different analyses and data in order to obtain a wider perspective of the evolution and diversification of this group of fishes.

Rhodopsin gene copies in Japanese eel originated in a teleost-specific genome duplication

BACKGROUND

Gene duplication is considered important to increasing the genetic diversity in animals. In fish, visual pigment genes are often independently duplicated, and the evolutionary significance of such duplications has long been of interest. Eels have two rhodopsin genes (rho), one of which (freshwater type, fw-rho) functions in freshwater and the other (deep-sea type, ds-rho) in marine environments. Hence, switching of rho expression in retinal cells is tightly linked with eels' unique life cycle, in which they migrate from rivers or lakes to the sea. These rho genes are apparently paralogous, but the timing of their duplication is unclear due to the deep-branching phylogeny. The aim of the present study is to elucidate the evolutionary origin of the two rho copies in eels using comparative genomics methods.

RESULTS

In the present study, we sequenced the genome of Japanese eel Anguilla japonica and reconstructed two regions containing rho by de novo assembly. We found a single corresponding region in a non-teleostean primitive ray-finned fish (spotted gar) and two regions in a primitive teleost (Asian arowana). The order of ds-rho and the neighboring genes was highly conserved among the three species. With respect to fw-rho, which was lost in Asian arowana, the neighboring genes were also syntenic between Japanese eel and Asian arowana. In particular, the pattern of gene losses in ds-rho and fw-rho regions was the same as that in Asian arowana, and no discrepancy was found in any of the teleost genomes examined. Phylogenetic analysis supports mutual monophyly of these two teleostean synteny groups, which correspond to the ds-rho and fw-rho regions.

CONCLUSIONS

Syntenic and phylogenetic analyses suggest that the duplication of rhodopsin gene in Japanese eel predated the divergence of eel (Elopomorpha) and arowana (Osteoglossomorpha). Thus, based on the principle of parsimony, it is most likely that the rhodopsin paralogs were generated through a whole genome duplication in the ancestor of teleosts, and have remained till the present in eels with distinct functional roles. Our result indicates, for the first time, that teleost-specific genome duplication may have contributed to a gene innovation involved in eel-specific migratory life cycle.

Morphology and Species Composition of Southern Adriatic Sea Leptocephali Evaluated Using DNA Barcoding

Leptocephali are the characteristic larvae of the superorder Elopomorpha that are difficult to identify at the species level. In this study, we used DNA barcoding (i.e. short genetic sequences of DNA used as unique species tags) coupled with classical taxonomic methods to identify leptocephali in the southern Adriatic Sea. This information will provide an assessment of the biodiversity of the eel larvae in this region. A total of 2,785 leptocephali were collected, and using external morphology were assigned to seven morphotypes: Ariosoma balearicum, Conger conger, Gnathophis mystax, Facciolella sp., Nettastoma melanurum, Dalophis imberbis and Chlopsis bicolor. Collectively, these seven morphotypes are considered to be a good proxy for the Anguilliformes community (the main order of the Elopomorpha) in the southern Adriatic Sea (to date, seven families and sixteen species have been recorded in this region). Interestingly, the higher number of G. mystax larvae collected suggests an increased abundance of this genus. To validate the morphological identifications, we sequenced 61 leptocephali (at a 655 bp fragment from the cytochrome oxidase subunit 1 mitochondrial region) and developed barcode vouchers for the seven morphotypes. Using genetic information from reference databases, we validated three of these morphotypes. Where reference sequences were unavailable, we generated barcodes for both adult and juvenile forms to provide additional genetic information. Using this integrated approach allowed us to characterize a new species of Facciolella in the Adriatic Sea for the first time. Moreover, we also revealed a lack of differentiation, at the species level, between G. mistax and G. bathytopos, a western Atlantic Ocean species. Our morphological and barcode data have been published in the Barcoding of the Adriatic Leptocephali database. This work represents the first contribution to a wider project that aims to create a barcode database to support the assessment of leptocephali diversity in the Mediterranean Sea.

A giant anguilliform leptocephalus Thalassenchelys foliaceus Castle & Raju is a junior synonym of Congriscus maldivensis (Norman 1939)

A single specimen of giant leptocephalus Thalassenchelys foliaceus Castle & Raju 1975 was caught in subtropical waters of the western North Pacific Ocean. Mitochondrial coI gene sequence divergence between T. foliaceus and Congriscus maldivensis (Norman 1939) was 0·64 ± 0·27% (mean ± s.e.), and the myomere and vertebral counts of these species were similar, indicating T. foliaceus is a junior synonym of C. maldivensis.

Three nuclear and two membrane estrogen receptors in basal teleosts, Anguilla sp.: Identification, evolutionary history and differential expression regulation

Estrogens interact with classical intracellular nuclear receptors (ESR), and with G-coupled membrane receptors (GPER). In the eel, we identified three nuclear (ESR1, ESR2a, ESR2b) and two membrane (GPERa, GPERb) estrogen receptors. Duplicated ESR2 and GPER were also retrieved in most extant teleosts. Phylogeny and synteny analyses suggest that they result from teleost whole genome duplication (3R). In contrast to conserved 3R-duplicated ESR2 and GPER, one of 3R-duplicated ESR1 has been lost shortly after teleost emergence. Quantitative PCRs revealed that the five receptors are all widely expressed in the eel, but with differential patterns of tissue expression and regulation. ESR1 only is consistently up-regulated in vivo in female eel BPG-liver axis during induced sexual maturation, and also up-regulated in vitro by estradiol in eel hepatocyte primary cultures. This first comparative study of the five teleost estradiol receptors provides bases for future investigations on differential roles that may have contributed to the conservation of multiple estrogen receptors.

Horizontal functional gene transfer from bacteria to fishes

Invertebrates can acquire functional genes via horizontal gene transfer (HGT) from bacteria but fishes are not known to do so. We provide the first reliable evidence of one HGT event from marine bacteria to fishes. The HGT appears to have occurred after emergence of the teleosts. The transferred gene is expressed and regulated developmentally. Its successful integration and expression may change the genetic and metabolic repertoire of fishes. In addition, this gene contains conserved domains and similar tertiary structures in fishes and their putative donor bacteria. Thus, it may function similarly in both groups. Evolutionary analyses indicate that it evolved under purifying selection, further indicating its conserved function. We document the first likely case of HGT of functional gene from prokaryote to fishes. This discovery certifies that HGT can influence vertebrate evolution.

Histozoic myxosporeans infecting the stomach wall of elopiform fishes represent a novel lineage, the Gastromyxidae

BACKGROUND

Traditional studies on myxosporeans have used myxospore morphology as the main criterion for identification and taxonomic classification, and it remains important as the fundamental diagnostic feature used to confirm myxosporean infections in fish and other vertebrate taxa. However, its use as the primary feature in systematics has led to numerous genera becoming polyphyletic in subsequent molecular phylogenetic analyses. It is now known that other features, such as the site and type of infection, can offer a higher degree of congruence with molecular data, albeit with its own inconsistencies, than basic myxospore morphology can reliably provide.

METHODS

Histozoic gastrointestinal myxosporeans from two elopiform fish from Malaysia, the Pacific tarpon Megalops cyprinoides and the ten pounder Elops machnata were identified and described using morphological, histological and molecular methodologies.

RESULTS

The myxospore morphology of both species corresponds to the generally accepted Myxidium morphotype, but both had a single nucleus in the sporoplasm and lacked valvular striations. In phylogenetic analyses they were robustly grouped in a discrete clade basal to myxosporeans, with similar shaped myxospores, described from gill monogeneans, which are located at the base of the multivalvulid clade. New genera Gastromyxum and Monomyxum are erected to accommodate these myxosporean taxa from fish and gill monogeneans respectively. Each are placed in a new family, the Gastromyxidae with Gastromyxum as the type genus and Monomyxidae with Monomyxum as the type genus.

CONCLUSIONS

To improve modern systematics of the myxosporeans it is clear that a combination of biological, ecological, morphological and molecular data should be used in descriptive studies, and the naming and redistribution of taxa and genera is going to be necessary to achieve this. Here we demonstrate why some Myxidium-shaped myxospores should not be included in the family Myxidiidae, and create two new families to accommodate them based on their site of infection, host biology / ecology, DNA sequence data and morphological observations. Subsequent descriptive works need to follow a similar course if we are going to create a prevailing and workable systematic structure for the Myxosporea.

Demonstration of the Coexistence of Duplicated LH Receptors in Teleosts, and Their Origin in Ancestral Actinopterygians

Pituitary gonadotropins, FSH and LH, control gonad activity in vertebrates, via binding to their respective receptors, FSHR and LHR, members of GPCR superfamily. Until recently, it was accepted that gnathostomes possess a single FSHR and a single LHR, encoded by fshr and lhcgr genes. We reinvestigated this question, focusing on vertebrate species of key-phylogenetical positions. Genome analyses supported the presence of a single fshr and a single lhcgr in chondrichthyans, and in sarcopterygians including mammals, birds, amphibians and coelacanth. In contrast, we identified a single fshr but two lhgcr in basal teleosts, the eels. We further showed the coexistence of duplicated lhgcr in other actinopterygians, including a non-teleost, the gar, and other teleosts, e.g. Mexican tetra, platyfish, or tilapia. Phylogeny and synteny analyses supported the existence in actinopterygians of two lhgcr paralogs (lhgcr1/ lhgcr2), which do not result from the teleost-specific whole-genome duplication (3R), but likely from a local gene duplication that occurred early in the actinopterygian lineage. Due to gene losses, there was no impact of 3R on the number of gonadotropin receptors in extant teleosts. Additional gene losses during teleost radiation, led to a single lhgcr (lhgcr1 or lhgcr2) in some species, e.g. medaka and zebrafish. Sequence comparison highlighted divergences in the extracellular and intracellular domains of the duplicated lhgcr, suggesting differential properties such as ligand binding and activation mechanisms. Comparison of tissue distribution in the European eel, revealed that fshr and both lhgcr transcripts are expressed in the ovary and testis, but are differentially expressed in non-gonadal tissues such as brain or eye. Differences in structure-activity relationships and tissue expression may have contributed as selective drives in the conservation of the duplicated lhgcr. This study revises the evolutionary scenario and nomenclature of gonadotropin receptors, and opens new research avenues on the roles of duplicated LHR in actinopterygians.

Duplicated leptin receptors in two species of eel bring new insights into the evolution of the leptin system in vertebrates

Since its discovery in mammals as a key-hormone in reproduction and metabolism, leptin has been identified in an increasing number of tetrapods and teleosts. Tetrapods possess only one leptin gene, while most teleosts possess two leptin genes, as a result of the teleost third whole genome duplication event (3R). Leptin acts through a specific receptor (LEPR). In the European and Japanese eels, we identified two leptin genes, and for the first time in vertebrates, two LEPR genes. Synteny analyses indicated that eel LEPRa and LEPRb result from teleost 3R. LEPRb seems to have been lost in the teleost lineage shortly after the elopomorph divergence. Quantitative PCRs revealed a wide distribution of leptins and LEPRs in the European eel, including tissues involved in metabolism and reproduction. Noticeably, leptin1 was expressed in fat tissue, while leptin2 in the liver, reflecting subfunctionalization. Four-month fasting had no impact on the expression of leptins and LEPRs in control European eels. This might be related to the remarkable adaptation of silver eel metabolism to long-term fasting throughout the reproductive oceanic migration. In contrast, sexual maturation induced differential increases in the expression of leptins and LEPRs in the BPG-liver axis. Leptin2 was strikingly upregulated in the liver, the central organ of the reproductive metabolic challenge in teleosts. LEPRs were differentially regulated during sexual maturation, which may have contributed to the conservation of the duplicated LEPRs in this species. This suggests an ancient and positive role of the leptin system in the vertebrate reproductive function. This study brings new insights on the evolutionary history of the leptin system in vertebrates. Among extant vertebrates, the eel represents a unique case of duplicated leptins and leptin receptors as a result of 3R.

Phylogenetic analysis of molecular and morphological data highlights uncertainty in the relationships of fossil and living species of Elopomorpha (Actinopterygii: Teleostei)

Elopomorpha is one of the three main clades of living teleost fishes and includes a range of disparate lineages including eels, tarpons, bonefishes, and halosaurs. Elopomorphs were among the first groups of fishes investigated using Hennigian phylogenetic methods and continue to be the object of intense phylogenetic scrutiny due to their economic significance, diversity, and crucial evolutionary status as the sister group of all other teleosts. While portions of the phylogenetic backbone for Elopomorpha are consistent between studies, the relationships among Albula, Pterothrissus, Notacanthiformes, and Anguilliformes remain contentious and difficult to evaluate. This lack of phylogenetic resolution is problematic as fossil lineages are often described and placed taxonomically based on an assumed sister group relationship between Albula and Pterothrissus. In addition, phylogenetic studies using morphological data that sample elopomorph fossil lineages often do not include notacanthiform or anguilliform lineages, potentially introducing a bias toward interpreting fossils as members of the common stem of Pterothrissus and Albula. Here we provide a phylogenetic analysis of DNA sequences sampled from multiple nuclear genes that include representative taxa from Albula, Pterothrissus, Notacanthiformes and Anguilliformes. We integrate our molecular dataset with a morphological character matrix that spans both living and fossil elopomorph lineages. Our results reveal substantial uncertainty in the placement of Pterothrissus as well as all sampled fossil lineages, questioning the stability of the taxonomy of fossil Elopomorpha. However, despite topological uncertainty, our integration of fossil lineages into a Bayesian time calibrated framework provides divergence time estimates for the clade that are consistent with previously published age estimates based on the elopomorph fossil record and molecular estimates resulting from traditional node-dating methods.

A multi-gene dataset reveals a tropical New World origin and Early Miocene diversification of croakers (Perciformes: Sciaenidae)

Widely distributed groups of living animals, such as the predominantly marine fish family Sciaenidae, have always attracted the attention of biogeographers to document the origins and patterns of diversification in time and space. In this study, the historical biogeography of the global Sciaenidae is reconstructed within a molecular phylogenetic framework to investigate their origin and to test the hypotheses explaining the present-day biogeographic patterns. Our data matrix comprises six mitochondrial and nuclear genes in 93 globally sampled sciaenid species from 52 genera. Within the inferred phylogenetic tree of the Sciaenidae, we identify 15 main and well-supported lineages; some of which have not been recognized previously. Reconstruction of habitat preferences shows repeated habitat transitions between marine and euryhaline environments. This implies that sciaenids can easily adapt to some variations in salinity, possibly as the consequence of their nearshore habitats and migratory life history. Conversely, complete marine/euryhaline to freshwater transitions occurred only three times, in South America, North America and South Asia. Ancestral range reconstruction analysis concomitant with fossil evidence indicates that sciaenids first originated and diversified in the tropical America during the Oligocene to Early Miocene before undergoing two range expansions, to Eastern Atlantic and to the Indo-West Pacific where a maximum species richness is observed. The uncommon biogeographic pattern identified is discussed in relation to current knowledge on origin of gradients of marine biodiversity toward the center of origin hypothesis in the Indo-West Pacific.

Deep sequencing of the olfactory epithelium reveals specific chemosensory receptors are expressed at sexual maturity in the European eel Anguilla anguilla

Vertebrate genomes encode a diversity of G protein-coupled receptor (GPCR) that belong to large gene families and are used by olfactory systems to detect chemical cues found in the environment. It is not clear however, if individual receptors from these large gene families have evolved roles that are specific to certain life stages. Here, we used deep sequencing to identify differentially expressed receptor transcripts in the olfactory epithelia (OE) of freshwater, seawater and sexually mature male eels (Anguilla anguilla). This species is particularly intriguing because of its complex life cycle, extreme long-distance migrations and early-branching position within the teleost phylogeny. In the A. anguillaOE, we identified full-length transcripts for 13, 112, 6 and 38 trace amine-associated receptors, odorant receptors (OR) and type I and type II vomeronasal receptors (V1R and V2R). Most of these receptors were expressed at similar levels at different life stages and a subset of OR and V2R-like transcripts was more abundant in sexually mature males suggesting that ORs and V2R-like genes are important for reproduction. We also identified a set of GPCR signal transduction genes that were differentially expressed indicating that eels make use of different GPCR signal transduction genes at different life stages. The finding that a diversity of chemosensory receptors is expressed in the olfactory epithelium and that a subset is differentially expressed suggests that most receptors belonging to large chemosensory gene families have functions that are important at multiple life stages, while a subset has evolved specific functions at different life stages.

FASconCAT-G: extensive functions for multiple sequence alignment preparations concerning phylogenetic studies

BACKGROUND

Phylogenetic and population genetic studies often deal with multiple sequence alignments that require manipulation or processing steps such as sequence concatenation, sequence renaming, sequence translation or consensus sequence generation. In recent years phylogenetic data sets have expanded from single genes to genome wide markers comprising hundreds to thousands of loci. Processing of these large phylogenomic data sets is impracticable without using automated process pipelines. Currently no stand-alone or pipeline compatible program exists that offers a broad range of manipulation and processing steps for multiple sequence alignments in a single process run.

RESULTS

Here we present FASconCAT-G, a system independent editor, which offers various processing options for multiple sequence alignments. The software provides a wide range of possibilities to edit and concatenate multiple nucleotide, amino acid, and structure sequence alignment files for phylogenetic and population genetic purposes. The main options include sequence renaming, file format conversion, sequence translation between nucleotide and amino acid states, consensus generation of specific sequence blocks, sequence concatenation, model selection of amino acid replacement with ProtTest, two types of RY coding as well as site exclusions and extraction of parsimony informative sites. Convieniently, most options can be invoked in combination and performed during a single process run. Additionally, FASconCAT-G prints useful information regarding alignment characteristics and editing processes such as base compositions of single in- and outfiles, sequence areas in a concatenated supermatrix, as well as paired stem and loop regions in secondary structure sequence strings.

CONCLUSIONS

FASconCAT-G is a command-line driven Perl program that delivers computationally fast and user-friendly processing of multiple sequence alignments for phylogenetic and population genetic applications and is well suited for incorporation into analysis pipelines.

Multiple thyrotropin β-subunit and thyrotropin receptor-related genes arose during vertebrate evolution

Thyroid-stimulating hormone (TSH) is composed of a specific β subunit and an α subunit that is shared with the two pituitary gonadotropins. The three β subunits derive from a common ancestral gene through two genome duplications (1R and 2R) that took place before the radiation of vertebrates. Analysis of genomic data from phylogenetically relevant species allowed us to identify an additional Tshβ subunit-related gene that was generated through 2R. This gene, named Tshβ2, present in cartilaginous fish, little skate and elephant shark, and in early lobe-finned fish, coelacanth and lungfish, was lost in ray-finned fish and tetrapods. The absence of a second type of TSH receptor (Tshr) gene in these species suggests that both TSHs act through the same receptor. A novel Tshβ sister gene, named Tshβ3, was generated through the third genomic duplication (3R) that occurred early in the teleost lineage. Tshβ3 is present in most teleost groups but was lostin tedraodontiforms. The 3R also generated a second Tshr, named Tshrb. Interestingly, the new Tshrb was translocated from its original chromosomic position after the emergence of eels and was then maintained in its new position. Tshrb was lost in tetraodontiforms and in ostariophysians including zebrafish although the latter species have two TSHs, suggesting that TSHRb may be dispensable. The tissue distribution of duplicated Tshβs and Tshrs was studied in the European eel. The endocrine thyrotropic function in the eel would be essentially mediated by the classical Tshβ and Tshra, which are mainly expressed in the pituitary and thyroid, respectively. Tshβ3 and Tshrb showed a similar distribution pattern in the brain, pituitary, ovary and adipose tissue, suggesting a possible paracrine/autocrine mode of action in these non-thyroidal tissues. Further studies will be needed to determine the binding specificity of the two receptors and how these two TSH systems are interrelated.