Evolutionary history of Otophysi (Teleostei), a major clade of the modern freshwater fishes: Pangaean origin and Mesozoic radiation

A phylogenomic approach to reconstruct interrelationships of main clupeocephalan lineages with a critical discussion of morphological apomorphies

BACKGROUND

Previous molecular studies on the phylogeny and classification of clupeocephalan fishes revealed numerous new taxonomic entities. For re-analysing these taxa, we perform target gene capturing and subsequent next generation sequencing of putative ortholog exons of major clupeocephalan lineages. Sequence information for the RNA bait design was derived from publicly available genomes of bony fishes. Newly acquired sequence data comprising > 800 exon sequences was subsequently used for phylogenetic reconstructions.

RESULTS

Our results support monophyletic Otomorpha comprising Alepocephaliformes. Within Ostariophysi, Gonorynchiformes are sister to a clade comprising Cypriniformes, Characiformes, Siluriformes and Gymnotiformes, where the interrelationships of Characiformes, Siluriformes and Gymnotiformes remain enigmatic. Euteleosts comprise four major clades: Lepidogalaxiiformes, Protacanthopterygii, Stomiatii, and Galaxiiformes plus Neoteleostei. The monotypic Lepidogalaxiiformes form the sister-group to all remaining euteleosts. Protacanthopterygii, comprising Argentini-, Esoci- and Salmoniformes, is sister to Stomiatii (Osmeriformes and Stomiatiformes) and Galaxiiformes plus Neoteleostei.

CONCLUSIONS

Several proposed monophyla defined by morphological apomorphies within the Clupeocephalan phylogeny are confirmed by the phylogenetic estimates presented herein. However, other morphologically described groups cannot be reconciled with molecular phylogenies. Thus, numerous morphological apomoprhies of supposed monophyla are called into question. The interpretation of suggested morphological synapomorphies of otomorph fishes is strongly affected by the inclusion of deep-sea inhabiting, and to that effect morphologically adapted Alepocephaliformes. Our revision of these potential synapomorphies, in the context that Alepocephaliformes are otomorph fishes, reveals that only a single character of the total nine characters proposed as synapomorphic for the group is clearly valid for all otomorphs. Three further characters remain possible apomorphies since their status remains unclear in the deep-sea adapted Alepocephaliformes showing developmental lag and lacking a swim bladder. Further, our analysis places Galaxiiformes as sister group to neoteleosts, which contradicts some previous molecular phylogenetic studies. This needs further investigation from a morphological perspective, as suggested synapomophies for several euteleostean lineages are challenged or still lacking. For the verification of results presented herein, a denser phylogenomic-level taxon sampling should be applied.

A history of the biogeography of Amazonian fishes

ABSTRACT The history of knowledge about Amazonian biogeography is as rich as its fish community, and a fascinating theme of study in itself. Several current paradigms and controversies about Amazonian fish biogeography are rooted in principles dating from the second half of the 18th to the first half of the 19th centuries. The present work establishes a relationship between current biogeographical ideas and their old predecessors, on the basis of a chronologically-oriented historical continuity analysis. The advent of evolutionary theory has not contributed significantly to a transformation of the knowledge on the biogeography of Amazonian fishes. On the other hand, the two main schools of biogeographical thought (dispersalist and vicariant) had major implications on how Amazonian fish distribution is interpreted. The process was gradual and many hypotheses have combined elements from each of the two schools. Chronologically, practically the entire tradition of Amazonian biogeography takes place within the evolutionary paradigm, although its founder Louis Agassiz was himself an anti-evolutionist. The birth of Amazonian biogeography is Agassiz´s travel in Amazon. That document makes it clear that the author did not consider dispersal as a valid explanation for the biogeographical patterns he found. Later, Carl Eigenmann helps to spread the dispersalist tradition as a model for biogeographical explanations in fish distributions, a phase which lasted until the late 20th century. A major shift occurs with the contributions of Marylin Weitzman, Stanley Weitzman and Richard Vari, who associated the temporal framework of phylogenetic hypotheses with distribution patterns, underscoring the predictive power of vicariant biogeography. The present-day paradigm begins with the work of John Lundberg and attempts to incorporate geomorphological and phylogenetic information into integrative biogeographical hypotheses. Some emblematic problems have delayed proposition of general hypotheses on the phylogenetic biogeography of South American fishes, such as the poor state of knowledge of their species-level systematics; an incomplete distributional record for most species and sparse or non-existent data on the phylogenetic history of most supraspecific taxa. Such drawbacks are now being corrected at a fast pace. Recent advances on geographical distribution and an increasing number of phylogenetic hypotheses will allow unprecedented large-scale biogeographic analyses, including those based on event models and Bayesian inference. Thus, the biogeography of South American fishes, especially Amazonian ones, should soon experiment a new age of progress. The success of that new phase will depend on its ability to recognize and segregate multiple overlapping temporal layers of hydrological changes, and to develop analytical tools that can deal with temporal mixing.

Otomorphs (= otocephalans or ostarioclupeomorphs) revisited

ABSTRACT A morphological revision is presented here on the cohort Otomorpha, a clade currently interpreted as the most primitive among the large supercohort Clupeocephala. Otomorpha is a morphologically heterogeneous group represented by clupei forms , alepocephaliforms, and ostariophysans (gonorynchiforms, cypriniforms, characiforms, siluriforms, and gymnoti forms) that inhabit various marine and freshwater environments worldwide. Otomorphs have a long (ca. 145 Ma) and diverse fossil record. They are the largest fish teleostean clade worldwide, as well as the largest of the Neotropical Region. While molecular studies strongly confirm the monophyly of Otomorpha, most potential morphological synapomorphies of the group become homoplastic largely due to the peculiar morphological character states (either losses or transformations) present in alepocephaliforms. The fusion of haemal arches with their respective vertebral centra anterior to preural centrum 2 stands as an unambiguous synapomorphy of the clade. The ankylosis or fusion of the extrascapular and parietal bones, and silvery areas associated with the gas bladder are also interpreted as synapomorphies, although they are homoplastic characters mainly due to secondary losses or further transformations of the morphological features in the alepocephaliforms.

Air- breathing in fish: Air- breathing organs and control of respiration: Nerves and neurotransmitters in the air-breathing organs and the skin

In fishes, exploitation of aerial gas exchange has evolved independently many times, involving a variety of air-breathing organs. Indeed, air-breathing occurs in at least 49 known families of fish (Graham, 1997). Many amphibious vertebrates, at some stage of their development are actually trimodal breathers that use various combinations of respiratory surfaces to breath both water (skin and/or gill) and air (skin and/or lung). The present review examines the evolutionary implications of air-breathing organs in fishes and the morphology of the peripheral receptors and the neurotransmitter content of the cells involved in the control of air-breathing. Control of breathing, whether gill ventilation or air-breathing, is influenced by feedback from peripheral and/or central nervous system receptors that respond to changes in PO2, PCO2 and/or pH. Although the specific chemoreceptors mediating the respiratory reflexes have not been conclusively identified, studies in water-breathing teleosts have implicated the neuroepithelial cells (NECs) existing in gill tissues as the O2 sensitive chemoreceptors that initiate the cardiorespiratory reflexes in aquatic vertebrates. Some of the air-breathing fishes, such as Protopterus, Polypterus and Amia have been shown to have NECs in the gills and/or lungs, although the role of these receptors and their innervation in the control of breathing is not known. NECs have been also reported in the specialized respiratory epithelia of accessory respiratory organs (ARO's) of some catfish species and in the gill and skin of the mudskipper Periophthalmodon schlosseri. Unlike teleosts matching an O2-oriented ventilation to ambient O2 levels, lungfishes have central and peripheral H+/CO2 receptors that control the acid-base status of the blood.

Whole Genome Sequencing of the Pirarucu (Arapaima gigas) Supports Independent Emergence of Major Teleost Clades

The Pirarucu (Arapaima gigas) is one of the world's largest freshwater fishes and member of the superorder Osteoglossomorpha (bonytongues), one of the oldest lineages of ray-finned fishes. This species is an obligate air-breather found in the basin of the Amazon River with an attractive potential for aquaculture. Its phylogenetic position among bony fishes makes the Pirarucu a relevant subject for evolutionary studies of early teleost diversification. Here, we present, for the first time, a draft genome version of the A. gigas genome, providing useful information for further functional and evolutionary studies. The A. gigas genome was assembled with 103-Gb raw reads sequenced in an Illumina platform. The final draft genome assembly was ∼661 Mb, with a contig N50 equal to 51.23 kb and scaffold N50 of 668 kb. Repeat sequences accounted for 21.69% of the whole genome, and a total of 24,655 protein-coding genes were predicted from the genome assembly, with an average of nine exons per gene. Phylogenomic analysis based on 24 fish species supported the postulation that Osteoglossomorpha and Elopomorpha (eels, tarpons, and bonefishes) are sister groups, both forming a sister lineage with respect to Clupeocephala (remaining teleosts). Divergence time estimations suggested that Osteoglossomorpha and Elopomorpha lineages emerged independently in a period of ∼30 Myr in the Jurassic. The draft genome of A. gigas provides a valuable genetic resource for further investigations of evolutionary studies and may also offer a valuable data for economic applications.

Phylogeny of the hillstream loach genus Mesonoemacheilus reveals widespread diversification through ancient drainage connections in the Western Ghats Biodiversity Hotspot

Rivers draining the Western Ghats (WG) mountain ranges in peninsular India harbor an exceptionally diverse, unique and evolutionarily distinct assemblage of lower vertebrates with high levels of endemism, attributed to their evolution and potentially long history of isolation during the Late Cretaceous or Early Tertiary. A molecular phylogeny of hillstream loaches of the genus Mesonoemacheilus endemic to the WG revealed the presence of four clades which we designate as 'species groups'. A consensus of various species delimitation methods indicates the likelihood of 'at least' seven more undescribed species within Mesonoemacheilus. Molecular clock analysis dates the basal clade around 38 mya in the Paleogene, and subsequent diversification in the Neogene and Quaternary periods resulting in the current genetic diversity. Biogeographic analysis suggests that vicariance events which separated the rivers on either side of the two geological barriers/gaps, the Palghat and Shencottah, in the Neogene, as well as range contractions and cladogenetic events contributed to the current patterns of diversity and distribution of this genus. Our results also provide preliminary indications on the interconnections and faunal exchange between historical river drainages in the WG region.

Back to the roots: Reducing evolutionary rate heterogeneity among sequences gives support for the early morphological hypothesis of the root of Siluriformes (Teleostei: Ostariophysi)

Catfishes (Teleostei: Siluriformes) are a highly diverse order within Ostariophysi that is distributed worldwide. At the base of this clade emerge three lineages with well-defined monophylies: Diplomystidae, Loricarioidei, and Siluroidei. Morphological phylogeny studies place the Diplomystidae as the earliest branching of these three lineages, but studies based on molecular phylogenetics consistently find the fast-evolving Loricarioidei instead. The high lineage evolutionary rate heterogeneity in this order and the fact that the lineage placed closest to the root in the molecular phylogenies is fast evolving, including many long branches, raises the possibility that the discrepancy between morphological and molecular phylogenies may be the result of a long branch attraction inference artifact. We test this hypothesis by using a 10-gene dataset to evaluate the arrangement of the three main siluriform lineages, and apply the LS³ and LS⁴ taxon sequence subsampling methods to reduce evolutionary rate heterogeneity among lineages. The initial and complete dataset supports the basal branching of Loricarioidei as in all previous molecular phylogenies, but once lineage rate heterogeneity is reduced with LS³ or LS⁴ through the removal of sequences disrupting homogeneity, the phylogeny shows Diplomystidae as the earliest branching group, with high supports, as proposed by morphological phylogeny. The result obtained with LS³, however, introduces the misplacement of one of the species with the highest amount of missing data, Scoloplax sp. Because the sequence sub-selection criterion of LS⁴ has been optimized to reduce data removal, the phylogeny resulting from the LS⁴-processed data is in agreement with the known intra-lineage relationships in addition to supporting the morphologically-based rooting hypothesis. Our results are the first instance in which a consensus between molecular and morphological phylogeny is reached concerning the root of this order.

Trunk dental tissue evolved independently from underlying dermal bony plates but is associated with surface bones in living odontode-bearing catfish

Although oral dental tissue is a vertebrate attribute, trunk dental tissue evolved in several extinct vertebrate lineages but is rare among living species. The question of which processes trigger dental-tissue formation in the trunk remains open, and would shed light on odontogenesis evolution. Extra-oral dental structures (odontodes) in the trunk are associated with underlying dermal bony plates, leading us to ask whether the formation of trunk bony plates is necessary for trunk odontodes to emerge. To address this question, we focus on Loricarioidei: an extant, highly diverse group of catfish whose species all have odontodes. We examined the location and cover of odontodes and trunk dermal bony plates for all six loricarioid families and 17 non-loricarioid catfish families for comparison. We inferred the phylogeny of Loricarioidei using a new 10-gene dataset, eight time-calibration points, and noise-reduction techniques. Based on this phylogeny, we reconstructed the ancestral states of odontode and bony plate cover, and find that trunk odontodes emerged before dermal bony plates in Loricarioidei. Yet we discovered that when bony plates are absent, other surface bones are always associated with odontodes, suggesting a link between osteogenic and odontogenic developmental pathways, and indicating a remarkable trunk odontogenic potential in Loricarioidei.

Evolutionary emergence of the rac3b/rfng/sgca regulatory cluster refined mechanisms for hindbrain boundaries formation

Developmental programs often rely on parallel morphogenetic mechanisms that guarantee precise tissue architecture. While redundancy constitutes an obvious selective advantage, little is known on how novel morphogenetic mechanisms emerge during evolution. In zebrafish, rhombomeric boundaries behave as an elastic barrier, preventing cell intermingling between adjacent compartments. Here, we identify the fundamental role of the small-GTPase Rac3b in actomyosin cable assembly at hindbrain boundaries. We show that the novel rac3b/rfng/sgca regulatory cluster, which is specifically expressed at the boundaries, emerged in the Ostariophysi superorder by chromosomal rearrangement that generated new cis-regulatory interactions. By combining 4C-seq, ATAC-seq, transgenesis, and CRISPR-induced deletions, we characterized this regulatory domain, identifying hindbrain boundary-specific cis-regulatory elements. Our results suggest that the capacity of boundaries to act as an elastic mesh for segregating rhombomeric cells evolved by cooption of critical genes to a novel regulatory block, refining the mechanisms for hindbrain segmentation.

Characterization and gene expression analysis of pacu (Piaractus mesopotamicus) inducible nitric oxide synthase (iNOS) following Aeromonas dhakensis infection

Nitric oxide (NO) is an important effector molecule which is involved in a myriad of biological processes, including immune responses against pathogens such as parasites, virus and bacteria. During the inflammatory processes in vertebrates, NO is produced by the inducible nitric oxide synthase (iNOS) enzyme in practically all nucleated cells to suppress or kill intracellular pathogens. The aim of the present study was to characterize the full coding region of the iNOS gene of pacu (Piaractus mesopotamicus), an economically and ecologically important South American fish species, and to analyze mRNA expression levels following intraperitoneal infection with the pathogenic bacterium Aeromonas dhakensis by means of quantitative real time PCR (qPCR). The results showed that the pacu iNOS transcript is 3237 bp in length, encoding a putative protein composed of 1078 amino acid residues. The amino acid sequence showed similarities ranging from 69.03% to 94.34% with other teleost fish and 57.70% with the human iNOS, with all characteristic domains and cofactor binding sites of the enzyme detected. Phylogenetic analysis showed that the iNOS from the red-bellied piranha, another South American characiform, was the closest related sequence to the pacu iNOS. iNOS transcripts were constitutively detected in the liver, spleen and head kidney, and there was a significant upregulation in the liver and spleen at 12, 24 and 48 h after infection with A. dhakensis. No significant variations were observed in the head kidney during the periods analyzed. These results show that iNOS expression was induced by A. dhakensis infection and suggest that this enzyme may be involved in the response to this bacterium in pacu.

Phylogenomic Systematics of Ostariophysan Fishes: Ultraconserved Elements Support the Surprising Non-Monophyly of Characiformes

Ostariophysi is a superorder of bony fishes including more than 10,300 species in 1100 genera and 70 families. This superorder is traditionally divided into five major groups (orders): Gonorynchiformes (milkfishes and sandfishes), Cypriniformes (carps and minnows), Characiformes (tetras and their allies), Siluriformes (catfishes), and Gymnotiformes (electric knifefishes). Unambiguous resolution of the relationships among these lineages remains elusive, with previous molecular and morphological analyses failing to produce a consensus phylogeny. In this study, we use over 350 ultraconserved element (UCEs) loci comprising 5 million base pairs collected across 35 representative ostariophysan species to compile one of the most data-rich phylogenies of fishes to date. We use these data to infer higher level (interordinal) relationships among ostariophysan fishes, focusing on the monophyly of the Characiformes-one of the most contentiously debated groups in fish systematics. As with most previous molecular studies, we recover a non-monophyletic Characiformes with the two monophyletic suborders, Citharinoidei and Characoidei, more closely related to other ostariophysan clades than to each other. We also explore incongruence between results from different UCE data sets, issues of orthology, and the use of morphological characters in combination with our molecular data. [Conserved sequence; ichthyology; massively parallel sequencing; morphology; next-generation sequencing; UCEs.].

Phylogenomic Perspective on the Relationships and Evolutionary History of the Major Otocephalan Lineages

The phylogeny of otocephalan fishes is the subject of broad controversy based on morphological and molecular evidence. The primary unresolved issue pertaining to this lineage relates to the origin of Characiphysi, especially the paraphyly of Characiformes. The considerable uncertainty associated with this lineage has precluded a greater understanding of the origin and evolution of the clade. Herein, a phylogenomic approach was applied to resolve this debate. By analyzing 10 sets of transcriptomic data generated in this study and 12 sets of high-throughput data available in public databases, we obtained 1,110 single-copy orthologous genes (935,265 sites for analysis) from 22 actinopterygians, including 14 otocephalan fishes from six orders: Clupeiformes, Gonorynchiformes, Cypriniformes, Siluriformes, Characiformes, and Gymnotiformes. Based on a selection of 125 nuclear genes screened from single-gene maximum likelihood (ML) analyses and sequence bias testing, well-established relationships among Otocephala were reconstructed. We suggested that Gymnotiformes are more closely related to Characiformes than to Siluriformes and Characiformes are possibly paraphyletic. We also estimated that Otocephala originated in the Early-Late Jurassic, which postdates most previous estimations, and hypothesized scenarios of the early historical biogeographies of major otocephalan lineages.

A local duplication of the Melanocortin receptor 1 locus in Astyanax

In this study, we report evidence of a novel duplication of Melanocortin receptor 1 (Mc1r) in the cavefish genome. This locus was discovered following the observation of excessive allelic diversity in a ∼820 bp fragment of Mc1r amplified via degenerate PCR from a natural population of Astyanax aeneus fish from Guerrero, Mexico. The cavefish genome reveals the presence of two closely related Mc1r open reading frames separated by a 1.46 kb intergenic region. One open reading frame corresponds to the previously reported Mc1r receptor, and the other open reading frame (duplicate copy) is 975 bp in length, encoding a receptor of 325 amino acids. Sequence similarity analyses position both copies in the syntenic region of the single Mc1r locus in 16 representative craniate genomes spanning bony fish (including Astyanax) to mammals, suggesting we discovered tandem duplicates of this important gene. The two Mc1r copies share ∼89% sequence similarity and, within Astyanax, are more similar to one another compared to other melanocortin family members. Future studies will inform the precise functional significance of the duplicated Mc1r locus and if this novel copy number variant may have adaptive significance for the Astyanax lineage.

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.

Reducing the information gap on Loricarioidei (Siluriformes) mitochondrial genomics

Background

The genetic diversity of Neotropical fish fauna is underrepresented in public databases. This distortion is evident for the order Siluriformes, in which the suborders Siluroidei and Loricarioidei share equivalent proportion of species, although far less is known about the genetics of the latter clade, endemic to the Neotropical Region. Recently, this information gap was evident in a study about the structural diversity of fish mitochondrial genomes, and hampered a precise chronological resolution of Siluriformes. It has also prevented molecular ecology investigations about these catfishes, their interactions with the environment, responses to anthropogenic changes and potential uses.

Results

Using high-throughput sequencing, we provide the nearly complete mitochondrial genomes for 26 Loricariidae and one Callichthyidae species. Structural features were highly conserved. A notable exception was identified in the monophyletic clade comprising species of the Hemiancistrus, Hypostomini and Peckoltia-clades, a ~60 nucleotide-long deletion encompassing the seven nucleotides at the 3′ end of the Conserved Sequence Block (CSB) D of the control region. The expression of mitochondrial genes followed the usual punctuation pattern. Heteroplasmic sites were identified in most species. The retrieved phylogeny strongly corroborates the currently accepted tree, although bringing to debate the relationship between Schizolecis guntheri and Pareiorhaphis garbei, and highlighting the low genetic variability within the Peckoltia-clade, an eco-morphologically diverse and taxonomically problematic group.

Conclusions

Herein we have launched the use of high-throughput mitochondrial genomics in the studies of the Loricarioidei species. The new genomic resources reduce the information gap on the molecular diversity of Neotropical fish fauna, impacting the capacity to investigate a variety of aspects of the molecular ecology and evolution of these fishes. Additionally, the species showing the partial CSB-D are candidate models to study the replication and transcription of vertebrate mitochondrial genome.

Electronic supplementary material

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

Testing for the Occurrence of Selective Episodes During the Divergence of Otophysan Fishes: Insights from Mitogenomics

How natural selection shapes biodiversity constitutes a topic of renewed interest during the last few decades. The division Otophysi comprises approximately two-thirds of freshwater fish diversity and probably underwent an extensive adaptive radiation derived from a single invasion of the supercontinent Pangaea, giving place to the evolution of the main five Otophysan lineages during a short period of time. Little is known about the factors involved in the processes that lead to lineage diversification among this group of fishes and identifying directional selection acting over protein-coding genes could offer clues about the processes acting on species diversification. The main objective of this study was to explore the otophysan mitochondrial genome evolution, in order to account for the possible signatures of selective events in this lineage, and to explore for the functional connotations of these molecular substitutions. Mainly, three different approaches were used: the "ω-based" BS-REL and MEME methods, implemented in the DATAMONKEY web server, and analysis of selection on amino acid properties, implemented in the software TreeSAAP. We found evidence of selective episodes along several branches of the evolutionary history of othophysan fishes. Analyses carried out using the BS-REL algorithm suggest episodic diversifying selection at basal branches of the otophysan lineage, which was also supported by analyses implemented in MEME and TreeSAAP. These results suggest that throughout the Siluriformes radiation, an important number of adaptive changes occurred in their mitochondrial genome. The metabolic consequences and ecological correlates of these molecular substitutions should be addressed in future studies.

The liver transcriptome of suckermouth armoured catfish (Pterygoplichthys anisitsi, Loricariidae): Identification of expansions in defensome gene families

Pterygoplichthys is a genus of related suckermouth armoured catfishes native to South America, which have invaded tropical and subtropical regions worldwide. Physiological features, including an augmented resistance to organic xenobiotics, may have aided their settlement in foreign habitats. The liver transcriptome of Pterygoplichthys anisitsi was sequenced and used to characterize the diversity of mRNAs potentially involved in the responses to natural and anthropogenic chemicals. In total, 66,642 transcripts were assembled. Among the identified defensome genes, cytochromes P450 (CYP) were the most abundant, followed by sulfotransferases (SULT), nuclear receptors (NR) and ATP binding cassette transporters (ABC). A novel expansion in the CYP2Y subfamily was identified, as well as an independent expansion of the CYP2AAs. Two expansions were also observed among SULT1. Thirty-two transcripts were classified into twelve subfamilies of NR, while 21 encoded ABC transporters. The diversity of defensome transcripts sequenced herein could contribute to this species' resistance to organic xenobiotics.

Genome-wide interrogation advances resolution of recalcitrant groups in the tree of life

Much progress has been achieved in disentangling evolutionary relationships among species in the tree of life, but some taxonomic groups remain difficult to resolve despite increasing availability of genome-scale data sets. Here we present a practical approach to studying ancient divergences in the face of high levels of conflict, based on explicit gene genealogy interrogation (GGI). We show its efficacy in resolving the controversial relationships within the largest freshwater fish radiation (Otophysi) based on newly generated DNA sequences for 1,051 loci from 225 species. Initial results using a suite of standard methodologies revealed conflicting phylogenetic signal, which supports ten alternative evolutionary histories among early otophysan lineages. By contrast, GGI revealed that the vast majority of gene genealogies supports a single tree topology grounded on morphology that was not obtained by previous molecular studies. We also reanalysed published data sets for exemplary groups with recalcitrant resolution to assess the power of this approach. GGI supports the notion that ctenophores are the earliest-branching animal lineage, and adds insight into relationships within clades of yeasts, birds and mammals. GGI opens up a promising avenue to account for incompatible signals in large data sets and to discern between estimation error and actual biological conflict explaining gene tree discordance.

Comparative histology in the liver and spleen of three species of freshwater teleost

ABSTRACT We describe and compare the histology of liver and spleen ofGeophagus brasiliensis (Perciformes), Hypostomus francisci (Siluriformes) and Hoplias aff. malabaricus (Characiformes), tropical freshwater fishes. InG. brasiliensisandH. aff. malabaricusthe hepatocytes were arranged in tubular form whereas in H. franciscithey cord-like. In all species, hepatocytes presented glycogen, but in G. brasiliensis and H. aff. malabaricus they showed strong stained for hemossiderin in the cytoplasm. InG. brasiliensis and H. aff. malabaricus, melanomacrophage centres (MMCs) were associated to hepatic structures and only in G. brasiliensis was observed intrahepatic exocrine pancreas. The spleen, in all species, was characterized by red and white pulp without boundary between the two regions, but only in H. francisci was recorded nodular organization in splenic parenchyma. The G. brasiliensisandH. aff. malabaricuspresented in the white pulp MMCs linked mainly to ellipsoids. Besides, we observed large MMCs in the spleen in relation to liver of G. brasiliensis and H. aff. malabaricus. In liver, highest values of reticular fibers and collagen were observed inG. brasiliensis. In spleen, highest values of reticular fibers and collagen were recorded inH. aff. malabaricusandH. francisci, respectively. Histological differences confirm the hypothesis that the phylogenetic distance is reflected in liver and spleen.

The complete mitochondrial genome of the Asian stinging catfish, Heteropneustes fossilis (Siluriformes, Heteropneustidae) and its comparison with other related fish species

The complete mitogenome of Heteropneustes fossilis is described using Ion Torrent (PGM sequencer), which showed it was 16,489 bp in size comprising 13 mRNAs, 22 tRNAs, 2 rRNA genes, and 858 bp as D-Loop control region, along with gene order and organization, being similar to most of the other related Siluriformes fish mitogenome of NCBI databases. The 20 RNAs were packed into a typical cloverleaf structure. The mitogenome in the present study has 99% similarity to the complete mitogenome sequence of H. fossilis mitogenome reported earlier and also would be helpful in understanding the population genetics, phylogenetics, and evolution of catfishes.

A Time-Calibrated Mitogenome Phylogeny of Catfish (Teleostei: Siluriformes)

A very significant part of the world's freshwater ichthyofauna is represented by ancient, exceptionally diverse and cosmopolitan ray-finned teleosts of the order Siluriformes. Over the years, catfish have been established as an exemplary model for probing historical biogeography at various scales. Yet, several tantalizing gaps still exist in their phylogenetic history, timeline and mode of diversification. Here, we re-examine the phylogeny of catfish by assembling and analyzing almost all publicly available mitogenome data. We constructed an ingroup matrix of 62 full-length mitogenome sequences from 20 catfish families together with four cypriniform outgroups, spanning 15,557 positions in total. Partitioned maximum likelihood analyses and Bayesian relaxed clock dating using fossil age constraints provide some useful and novel insights into the evolutionary history of this group. Loricarioidei are recovered as the first siluriform group to diversify, rendering Neotropics the cradle of the order. The next deepest clade is the South American Diplomystoidei placed as a sister group to all the remaining Siluroidei. The two multifamilial clades of "Big Asia" and "Big Africa" are also recovered, albeit nodal support for the latter is poor. Within "Big Asia", Bagridae are clearly polyphyletic. Other interfamilial relationships, including Clariidae + Heteropneustidae, Doradidae + Auchenipteridae and Ictaluridae + Cranoglanididae are robustly resolved. Our chronogram shows that siluriforms have a Pangaean origin, at least as far back as the Early Cretaceous. The inferred timeline of the basal splits corroborates the "Out-of-South America" hypothesis and accords well with the fossil record. The divergence of Siluroidei most likely postdated the final separation of Africa and South America. An appealing case of phylogenetic affinity elaborated by biogeographic dispersal is exemplified by the Early Paleogene split between the Southeast Asian Cranoglanididae and Ictaluridae, with the latter radiating into North America's freshwater realm by Eocene. The end of Cretaceous probably concludes the major bout of diversification at the family level while with the dawn of the Cenozoic a prolific radiation is evident at the generic level.

A second visual rhodopsin gene, rh1-2, is expressed in zebrafish photoreceptors and found in other ray-finned fishes

Rhodopsin (rh1) is the visual pigment expressed in rod photoreceptors of vertebrates that is responsible for initiating the critical first step of dim-light vision. Rhodopsin is usually a single copy gene; however, we previously discovered a novel rhodopsin-like gene expressed in the zebrafish retina, rh1-2, which we identified as a functional photosensitive pigment that binds 11-cis retinal and activates in response to light. Here, we localized expression of rh1-2 in the zebrafish retina to a subset of peripheral photoreceptor cells, which indicates a partially overlapping expression pattern with rh1 We also expressed, purified and characterized Rh1-2, including investigation of the stability of the biologically active intermediate. Using fluorescence spectroscopy, we found the half-life of the rate of retinal release of Rh1-2 following photoactivation to be more similar to that of the visual pigment rhodopsin than to the non-visual pigment exo-rhodopsin (exorh), which releases retinal around 5 times faster. Phylogenetic and molecular evolutionary analyses show that rh1-2 has ancient origins within teleost fishes, is under similar selective pressure to rh1, and likely experienced a burst of positive selection following its duplication and divergence from rh1 These findings indicate that rh1-2 is another functional visual rhodopsin gene, which contradicts the prevailing notion that visual rhodopsin is primarily found as a single copy gene within ray-finned fishes. The reasons for retention of this duplicate gene, as well as possible functional consequences for the visual system, are discussed.

The complete mitochondrial genome of the threatened Neotropical catfish Lophiosilurus alexandri (Silurifomes: Pseudopimelodidae) and phylogenomic analysis indicate monophyly of Pimelodoidea

Lophiosilurus alexandri is an endemic catfish from the São Francisco River Basin (Brazil) popularly known as pacamã, which has economic potential for aquaculture farming. The mitochondrial genome was sequenced for the threatened Neotropical catfish L. alexandri. Assembly into scaffolds using MIRA and MITObim software produced the whole, circularized mitochondrial genome, which comprises 16,445 bp and presents the typical gene arrangement of Teleostei mitochondria. A phylogenomic analysis was performed after the concatenation of all proteins obtained from whole mitogenomes of 20 Siluriformes and two outgroups. The results confirmed the monophyly of nine families of catfishes and also clustered L. alexandri as a sister group to the family Pimelodidae, thus confirming the monophyly of the superfamily Pimelodoidea. This is the first mitochondrial phylogenomics study for Pimelodoidea and the first mitogenome described for the Pseudopimelodidae family, representing an important resource for phylogeography, evolutionary biology, and conservation genetics studies in Neotropical fishes.

Evolution of opercle bone shape along a macrohabitat gradient: species identification using mtDNA and geometric morphometric analyses in neotropical sea catfishes (Ariidae)

Transitions between the marine and freshwater macrohabitat have occurred repeatedly in the evolution of teleost fishes. For example, ariid catfishes have moved from freshwater to marine environments, and vice versa. Opercles, a skeletal feature that has been shown to change during such transitions, were subjected to 2D geometric morphometric analyses in order to investigate evolutionary shape changes during habitat transition in ariid catfishes and to test the influence of habitat on shape changes. A mtDNA marker, which proved useful in previous studies, was used to verify species identities. It greatly improved the assignment of specimens to a species, which are difficult to assign by morphology alone. The application of a mtDNA marker confirmed the occurrence of Notarius biffi in Central America, South of El Salvador. Molecular identification together with principal component analysis (PCA) and further morphological inspection of neurocrania indicated the existence of a cryptic species within Bagre pinnimaculatus. Principal component (PC) scores of individual specimens clustered in morphospace by genus rather than by habitat. Strong phylogenetic structure was detected using a permutation test of PC scores of species means on a phylogenetic tree. Calculation of Pagel's λ suggested that opercle shape evolved according to a Brownian model of evolution. Yet canonical variate analysis (CVA) conducted on the habitat groups showed significant differences in opercle shapes among freshwater and marine species. Overall, opercle shape in tropical American Ariidae appears to be phylogenetically constrained. This verifies the application of opercle shape as a taxonomic tool for species identification in fossil ariid catfishes. At the same time, adaptation to freshwater habitats shows characteristic opercle shape trajectories in ariid catfishes, which might be used to detect habitat preferences in fossils.

Mitogenomic perspectives on the origin of Tibetan loaches and their adaptation to high altitude

Tibetan loaches are the largest group of Tibetan fishes and are well adapted to the Tibetan Plateau. To investigate the origin of Tibetan loaches and their adaptations to the Tibetan Plateau, we determined 32 complete mitochondrial genomes that included 29 Tibetan loach species, two Barbatula species and Schistura longus. By combining these newly determined sequences with other previously published mitochondrial genomes, we assembled a large mitogenomic data set (11,433 bp) of 96 species in the superfamily Cobitoidea, to investigate the phylogenetic status of the genus Triplophysa. The resulting phylogeny strongly supported that the genus Triplophysa forms a monophyletic group within Nemacheilidae. Our molecular dating time suggests that the lineage leading to the Tibetan loaches and other loaches diverged approximately 23.5 Ma, which falls within the period of recent major uplifts of the Tibetan Plateau in the Early Miocene. Selection analyses revealed that the mitochondrial protein-coding genes of Tibetan loaches have larger ratios of nonsynonymous to synonymous substitutions than do those of non-Tibetan loaches, indicating that Tibetan loaches accumulated more nonsynonymous mutations than non-Tibetan loaches and exhibited rapid evolution. Two positively selected sites were identified in the ATP8 and ND1 genes.

Mitochondrial genomes of the South American electric knifefishes (Order Gymnotiformes)

Three complete mitochondrial genomes of South American electric fishes (Gymnotiformes), derived from high-throughput RNA sequencing (RNA-Seq), are reported herein. We report the complete mitochondrial genome of the bluntnose knifefish Brachyhypopomus n.sp. VERD, determined from newly sequenced data. We also provide the complete mitochondrial genomes for Sternopygus arenatus and the electric eel Electrophorus electricus, assembled from previously published transcriptome data. The mitochondrial genomes of Brachyhypopomus n.sp. VERD, Sternopygus arenatus and Electrophorus electricus have 13 protein-coding genes, 1 D-loop, 2 ribosomal RNAs and 22 transfer RNAs, and are 16,547, 16,667 and 16,906 bp in length, respectively. Phylogenetic analysis of the eight available mitochondrial genomes of gymnotiform fishes shows Apteronotus to be the sister lineage of other gymnotiformes, contradicting the “Sinusoidea” hypothesis that Apteronotidae and Sternopygidae are sister taxa.

Molecular evolution of cryptochromes in fishes

Circadian rhythmicity is an endogenous biological cycle of about 24h, which exists in cyanobacteria and fungi, plants and animals. Circadian rhythms improve the adaptability of organisms in both constant and changing environments. The cryptochrome (CRY) is a key element of the circadian system in various animal groups including fishes. We studied evolution of cryptochromes in the phylogenetically and ecologically diverse fish taxa. The phylogenetic tree of fish Cry features two major clades: Cry1 and Cry2. Teleosts possess extra copies of Cry1 due to the genome duplication, which resulted in 3 main paralogous subfamilies (1A, 1B and 1C). Cry1 experienced further diversification through additional duplications in some taxa. 1A of Cry1 is more conserved than the other paralogs (dN=0.010 ± 0.003, π=0.119 ± 0.058). The analysis of selection indicated that, while the Cry homologs in fish evolved under the different levels of selection pressure, strong purifying selection (average ω=0.017) dominated in their evolution.

Structural and functional evolution of 2',3'-cyclic nucleotide 3'-phosphodiesterase

2',3'-cyclic nucleotide 3'-phosphodiesterase (CNPase) is an abundant membrane-associated enzyme within the vertebrate myelin sheath. While the physiological function of CNPase still remains to be characterized in detail, it is known - in addition to its in vitro enzymatic activity - to interact with other proteins, small molecules, and membrane surfaces. From an evolutionary point of view, it can be deduced that CNPase is not restricted to myelin-forming cells or vertebrate tissues. Its evolution has involved gene fusion, addition of other small segments with distinct functions, such as membrane attachment, and possibly loss of function at the polynucleotide kinase-like domain. Currently, it is unclear whether the enzymatic function of the conserved phosphodiesterase domain in vertebrate myelin has a physiological role, or if CNPase could actually function - like many other classical myelin proteins - in a more structural role. This article is part of a Special Issue entitled SI: Myelin Evolution.

A Comparative Transcriptome Analysis between Wild and Albino Yellow Catfish (Pelteobagrus fulvidraco)

Body colours are important and striking features for individual survival and reproductive success, in particular in vertebrates where mating behaviour and mate preference may be strongly influenced by non-normal phenotypes. Pigmentation disorders may be generated by disruption of one or many independent genes as well as by environmental factors. The first discovery of albino yellow catfish (Pelteobagrus fulvidraco Richardson) with golden skin colour from fish farms in China provides us valuable material to study the molecular mechanism underlying the abnormalities of pigmentation. In this study, transcriptome sequencing of fin tissues corresponding to the distinct body colours, wild type and mutant albino yellow catfish, were performed using Illumina sequencing technology. Based on next-generation sequencing technology and de novo assembly, we generated a transcriptome of P. fulvidraco. A number of genes differentially expressed between the wild types and albinos were identified, suggesting their contribution to the different phenotypes and fitness. However, non-synonymous mutations result from single nucleotide substitutions residing in coding regions may not contribute to such differences. Based on the high-throughput expression data generated for the two different types of P. fulvidraco, we found that alterations of expression pattern may be more common than non-synonymous mutations. The transcriptome of P. fulvidraco will be an invaluable resource for subsequent comparative genomics and evolutionary analyses of this economically important fish.

'Fish' (Actinopterygii and Elasmobranchii) diversification patterns through deep time

Actinopterygii (ray-finned fishes) and Elasmobranchii (sharks, skates and rays) represent more than half of today's vertebrate taxic diversity (approximately 33000 species) and form the largest component of vertebrate diversity in extant aquatic ecosystems. Yet, patterns of 'fish' evolutionary history remain insufficiently understood and previous studies generally treated each group independently mainly because of their contrasting fossil record composition and corresponding sampling strategies. Because direct reading of palaeodiversity curves is affected by several biases affecting the fossil record, analytical approaches are needed to correct for these biases. In this review, we propose a comprehensive analysis based on comparison of large data sets related to competing phylogenies (including all Recent and fossil taxa) and the fossil record for both groups during the Mesozoic-Cainozoic interval. This approach provides information on the 'fish' fossil record quality and on the corrected 'fish' deep-time phylogenetic palaeodiversity signals, with special emphasis on diversification events. Because taxonomic information is preserved after analytical treatment, identified palaeodiversity events are considered both quantitatively and qualitatively and put within corresponding palaeoenvironmental and biological settings. Results indicate a better fossil record quality for elasmobranchs due to their microfossil-like fossil distribution and their very low diversity in freshwater systems, whereas freshwater actinopterygians are diverse in this realm with lower preservation potential. Several important diversification events are identified at familial and generic levels for elasmobranchs, and marine and freshwater actinopterygians, namely in the Early-Middle Jurassic (elasmobranchs), Late Jurassic (actinopterygians), Early Cretaceous (elasmobranchs, freshwater actinopterygians), Cenomanian (all groups) and the Paleocene-Eocene interval (all groups), the latter two representing the two most exceptional radiations among vertebrates. For each of these events along with the Cretaceous-Paleogene extinction, we provide an in-depth review of the taxa involved and factors that may have influenced the diversity patterns observed. Among these, palaeotemperatures, sea-levels, ocean circulation and productivity as well as continent fragmentation and environment heterogeneity (reef environments) are parameters that largely impacted on 'fish' evolutionary history, along with other biotic constraints.

The complete mitogenome of the Sheatfish Pterocryptis cochinchinensis (Siluriformes: Siluridae) and phylogenetic implications

To better understand the phylogenetic status of the Sheatfish, Pterocryptis cochinchinensis, we determined its complete mitogenome sequence using polymerase chain reaction and the direct sequencing method. The complete mitogenome sequence was 16,501bp in length and contained 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and 1 control region (D-loop), the gene was identical to typical vertebrate. This was the first report of the complete mitogenome sequence in genus Pterocryptis. Phylogenetic based on mitogenome sequences and cytb gene suggested that the Sheatfish formed the most basal branch having sister relationship with the clade containing all other analyzed genus Silurus fishes.

Complete mitochondrial genome sequence of Heteropneustes fossilis obtained by paired end next generation sequencing

In the present study, the complete mitochondrial genome sequence of Heteropneustes fossilis is reported using massive parallel sequence technology. The complete mitogenome of H. fossilis is obtained by de novo assembly of paired end Illumina sequences using CLC Genomics Workbench version 7.0.4, which is 16,489 bp in length. It comprised of 13 protein- coding genes, 22 tRNAs, 2 rRNA genes and a putative control region along with the gene order and organization, being similar to most of vertebrates. The mitogenome in the present study has 99% similarity to the complete mitogneome sequence of H. Fossilis, as reported earlier. Phylogenetic analysis of Siluriformes depicted that Heteropneustids were closer to Clariids. The mitogenome sequence of H. fossilis contributes better understanding of population genetics, phylogenetics and evolution of Indian catfish species.

Complete mitochondrial genome of the Neotropical catfish Pseudoplatystoma magdaleniatum (Siluriformes, Pimelodidae)

The Neotropical freshwater fish Pseudoplatystoma magdaleniatum is a trans-Andean species that belongs to the family of long-whiskered catfishes (family Pimelodidae). In this study, the complete mitochondrial genome of P. magdaleniatum was sequenced using the MiSeq Illumina platform. The complete circular mitogenome is 16,568 bp in length, exhibiting an average GC content of 44.19% and codes for 13 proteins, 2 ribosomal RNA genes and 22 transfer RNA genes. Additionally, it exhibits perfect synteny and similar length with the mitogenome of Pimelodus pictus.

Development of the splanchnocranium in Prochilodus argenteus (Teleostei: Characiformes) with a discussion of the basal developmental patterns in the Otophysi

Development of the mandibular, hyoid and gill arches, which constitute the splanchnocranium, are described for Prochilodus argenteus, order Characiformes, one of the basal lineages of the Otophysi. Development was examined from just hatched larvae through juveniles using whole specimens cleared and counterstained for cartilage and bone as well as histological preparations. Observations are compared with the developmental trends reported for Cypriniformes, the basalmost clade of the Otophysi. Shortened developmental sequences for Prochilodus compared to the cypriniform Catostomus were discovered in the ontogeny of the ceratohyals, ceratobranchials 1-5, epibranchials 1-4 and the symplectic portion of the hyosymplectic. Prochilodus also differs from Catostomus in having the basihyal plus the anterior copula appearing at different stages of ontogeny rather than simultaneously. Contrary to previous assumptions, developmental information indicates that hypobranchial 4 as well as likely basibranchial 5 are present in Prochilodus. Various developmental patterns in Prochilodus considered basal for the Otophysi, the predominant component of the Ostariophysi, are likely conserved from patterns prevalent in basal groups in the Actinopterygii.

Mitochondrial genome of the Trans-Andean shovelnose catfish Sorubim cuspicaudus (Siluriformes, Pimelodidae)

The Trans-Andean shovelnose catfish Sorubim cuspicaudus is the largest species within the genus Sorubim. In this work, the pyrosequencing technology was used to obtain the complete mitochondrial genome sequence of S. cuspicaudus. The 16,544 bp molecule contains 13 protein-coding genes, 22 tRNAs, 2 ribosomal RNAs and exhibit perfect synteny with other South-American catfishes.

First complete mitochondrial genome for any anostomid fish: Leporinus piavussu, a recently described piracema species

The first complete mitochondrial genome for the Anostomidae fish family has been announced. Piracema fishes are the potamodromous migratory species from South America, which undergo lengthy and dramatic yearly reproductive upstream runs. The piracema species Leporinus piavussu has recently been described after a long misidentification history. Its mitogenome, assembled using NGS data and Sanger sequencing, consists in a 16,682 bp circular molecule (GenBank accession number KM886569). The exact sequence and position of 37 mitochondrial genes and the control region is established. A possible case of heteroplasmy was found with NGS and corroborated by Sanger sequencing. These results will positively contribute to the debate about this group's taxonomy, evolution and conservation.

Complete mitogenome of the Neotropical fish Brycon henni, Eigenmann 1913 (Characiformes, Bryconidae)

Brycon henni (Characiformes, Bryconidae) is a trans-Andean omnivorous fish considered as promising species in aquaculture of temperate zones. In this work we report for first time the complete mitochondrial genome of B. henni. This mitogenome spans 16,885 bp and encodes 13 proteins, 22 tRNAs and 2 ribosomal RNAs. Additionally, the synteny is conserved with others species of the order Characiformes.

Mitochondrial genome of the Neotropical detritivore fish Curimata mivartii Steindachner 1878 (Characiformes, Curimatidae)

Curimata mivartii is a detritivorous species Neotropical, which is of great importance to local fisheries and riverine ecosystems. In this work, the mitochondrial genome of C. mivartii was completely sequenced using a combination of 454 FLX(+) platform and Sanger/capillary sequencing. The mitogenome is 16,705 bp in length, encodes 13 proteins, 22 tRNAs, 2 ribosomal RNAs and exhibit perfect synteny with others Characiformes.

Drivers of parasite sharing among Neotropical freshwater fishes

Because host-parasite interactions are so ubiquitous, it is of primary interest for ecologists to understand the factors that generate, maintain and constrain these associations. Phylogenetic comparative studies have found abundant evidence for host-switching to relatively unrelated hosts, sometimes related to diversification events, in a variety of host-parasite systems. For Monogenoidea (Platyhelminthes) parasites, it has been suggested that the co-speciation model alone cannot explain host occurrences, hence host-switching and/or non-vicariant modes of speciation should be associated with the origins and diversification of several monogenoid taxa. The factors that shape broad patterns of parasite sharing were investigated using path analysis as a way to generate hypotheses about the origins of host-parasite interactions between monogenoid gill parasites and Neotropical freshwater fishes. Parasite sharing was assessed from an interaction matrix, and explanatory variables included phylogenetic relationships, environmental preferences, biological traits and geographic distribution for each host species. Although geographic distribution of hosts and host ecology are important factors to understand host-parasite interactions, especially within host lineages that share a relatively recent evolutionary history, phylogeny had the strongest overall direct effect on parasite sharing. Phylogenetic contiguity of host communities may allow a 'stepping-stone' mode of host-switching, which increases parasite sharing. Our results reinforce the importance of including evolutionary history in the study of ecological associations, including emerging infectious diseases risk assessment.