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

Fifth mass extinction event triggered the diversification of the largest family of freshwater gastropods (Caenogastropoda: Truncatelloidea: Hydrobiidae)

The fifth mass extinction event (MEE) at the Cretaceous-Palaeogene (K-Pg) boundary 66 million years ago (Ma) led to massive species loss but also triggered the diversification of higher taxa. Five models have been proposed depending on whether this diversification occurred before, during or after the K-Pg boundary and the rate of species accumulation. While the effects of the K-Pg MEE on vertebrate evolution are relatively well understood, the impact on invertebrates, particularly in freshwater ecosystems, remains controversial. One example is the hyperdiverse Hydrobiidae-the most species-rich family of freshwater gastropods. Whereas some studies place its origin in the Jurassic or even Carboniferous, most fossil records postdate the K-Pg event. We therefore used robustly time-calibrated multi-locus phylogenies of >400 species representing >100 hydrobiid genera to unravel its evolutionary history and patterns of diversification. We found that the family started diversifying shortly after the K-Pg boundary (∼60 Ma; 95% highest posterior density 52-69 Ma). Lineage richness gradually increased to the present and phylogenetic diversity until ∼25 Ma. These findings suggest that diversification was not initially driven by ecological opportunity. Combining the two criteria of timing and rate of diversification, a soft-explosive diversification model of aquatic vertebrates best fits the patterns observed. We also show that most higher hydrobiid taxa (i.e. subfamilies) diversified from the Middle Oligocene to Middle Miocene (i.e. 12-28 Ma). Two of the 15 major clades delimited are described here as new subfamilies (i.e. Bullaregiinae n. subfam. and Pontobelgrandiellinae n. subfam.), whose members are restricted to subterranean waters. Our results are an important contribution to understanding how the fifth MEE has shaped evolution and patterns of biodiversity in continental aquatic systems. Given the high extinction risks faced by many hydrobiids today, they also emphasise the need to study the biodiversity of vulnerable ecosystems.

Mitophylogeny of Pangasiid Catfishes and its Taxonomic Implications for Pangasiidae and the Suborder Siluroidei

Pangasiidae (catfish order: Siluriformes) comprises 30 valid catfish species in four genera: Pangasius, Pangasianodon, Helicophagus, and Pseudolais. Their systematics are frequently revised due to the addition of newly described species. Although Pangasiidae is known to be a monophyletic family, the generic and phylogenetic relationships among the taxa are poorly resolved. This study characterized three newly obtained complete mitogenomes of Mekong River catfishes from Vietnam (Pangasius mekongensis, Pangasius krempfi, and Pangasianodon hypophthalmus), as well as the inter-and intrafamilial relationships of the Pangasiidae and catfish families in Siluroidei. The genomic features of their mitogenomes were similar to those of previously reported pangasiids, including all regulatory elements, extended terminal associated sequences (ETAS), and conserved sequence blocks (CSBs) (CSB-1, CSB-2, CSB-3, and CSBs, A to F) in the control region. A comprehensive phylogeny constructed from datasets of multiple 13 PCG sequences from 117 complete mitogenomes of 32 recognized siluriform families established Pangasiidae as monophyletic and a sister group of Austroglanididae. The [Pangasiidae + Austroglanididae] + (Ictaluridae + Cranoglanididae) + Ariidae] clade is a sister to the "Big Africa" major clade of Siluriformes. Furthermore, both phylogenies constructed from the single barcodes (83 partial cox1 and 80 partial cytB, respectively) clearly indicate genus relationships within Pangasiidae. Pangasianodon was monophyletic and a sister to the (Pangasius + Helicophagus + Pseudolais) group. Within the genus Pangasius, P. mekongensis was placed as a sister taxon to P. pangasius. Pangasius sanitwongsei was found to be related to and grouped with Pangasianodon, but in single-gene phylogenies, it was assigned to the Pangasius + Helicophagus + Pseudolais group. The datasets in this study are useful for studying pangasiid systematics, taxonomy and evolution.

Characterization of the male-specific region containing the candidate sex-determining gene in Amur catfish (Silurus asotus) using third-generation- and pool-sequencing data

Amur catfish (Silurus asotus) is an ecologically and economically important fish species in Asia. Here, we assembled the female and male Amur catfish genomes, with genome sizes of 757.15 and 755.44 Mb, respectively, at the chromosome level using nanopore and Hi-C technologies. Consistent with the known diploid chromosome count, both genomes contained 29 chromosome-size scaffolds covering 98.80 and 98.73 % of the complete haplotypic assembly with scaffold N50 of 28.87 and 27.29 Mb, respectively. The female (n = 40) and male (n = 40) pools were re-sequenced. Comparative analysis of sequencing and re-sequencing data from both sexes confirmed the presence of an XX/XY sex determination system in Amur catfish and revealed Chr5 as the sex chromosome containing an approximately 400 kb Y-specific region (MSY). Gene annotation revealed a male-specific duplicate of amhr2, namely amhr2y, in MSY, which is male-specific in different wild populations and expressed only in the testes. Amur catfish shared partially syntenic MSY and amhr2y genes with the southern catfish (S. meridionalis, Chr24), which were located on different chromosomes. High sequence divergence between amhr2y and amhr2 and high sequence similarity with amhr2y were observed in both species. These results indicate the common origin of the sex-determining (SD) gene and transition of amhr2y in the two Silurus species. Accumulation of repetitive elements in the MSY of both species may be the main driver of the transition of amhr2y. Overall, our study provides valuable catfish genomic resources. Moreover, determination of amhr2y as the candidate SD gene in Amur catfish provides another example of amhr2 as the SD gene in fish.

Effective practices for thermal tolerance polygon experiments using mottled catfish Corydoras paleatus

Temperature is an important environmental factor that affects how organisms allocate metabolic resources to physiological processes. Laboratory experiments that determine absolute thermal limits for representative species are important for understanding how fishes are affected by climate change. Critical Thermal Methodology (CTM) and Chronic Lethal Methodology (CLM) experiments were utilized to construct a complete thermal tolerance polygon for the South American fish species, Mottled catfish (Corydoras paleatus). Mottled catfish showed Chronic Lethal Maxima (CLMax) of 34.9 ± 0.52 °C and Chronic Lethal Minima (CLMin) of 3.8 ± 0.08 °C. Fish were chronically acclimated (∼2 weeks) to 6 temperatures ranging from 7.2 ± 0.05 °C →32.2 ± 0.16 °C (7 °C, 12 °C, 17 °C, 22 °C, 27 °C, and 32 °C), and CTM used to estimate upper and lower acute temperature tolerance. Linear regressions of Critical Thermal Maxima (CTMax) and Minima (CTMin) data with each acclimation temperature were used along with CLMax and CLMin to create a complete thermal tolerance polygon. The highest CTMax was 38.4 ± 0.60 °C for fish acclimated to 32.2 ± 0.16 °C, and the lowest CTMin was 3.36 ± 1.84 °C for fish acclimated to 7.2 ± 0.05 °C. Mottled catfish have a polygon measuring 785.7°C2, and the slope of the linear regressions showed the species gained 0.55 °C and 0.32 °C of upper and lower tolerance per degree of acclimation temperature, respectively. We compared slopes of CTMax or CTMin regression lines to each other using a set of comparisons between 3, 4, 5, or 6 acclimation temperatures. Our data demonstrated that 3 acclimation temperatures were as sufficient as 4 → 6 to pair with estimates of chronic upper and lower thermal limits for accurately determining a complete thermal tolerance polygon. Construction of this species' complete thermal tolerance polygon provides a template for other researchers. The following is sufficient to generate a complete thermal tolerance polygon: Three chronic acclimation temperatures that are spread somewhat evenly across a species' thermal range, include an estimation of CLMax and CLMin, and are followed by CTMax and CTMin measurements.

Multi-locus phylogeny of the catfish genus Ictalurus Rafinesque, 1820 (Actinopterygii, Siluriformes) and its systematic and evolutionary implications

BACKGROUND

Ictalurus is one of the most representative groups of North American freshwater fishes. Although this group has a well-studied fossil record and has been the subject of several morphological and molecular phylogenetic studies, incomplete taxonomic sampling and insufficient taxonomic studies have produced a rather complex classification, along with intricate patterns of evolutionary history in the genus that are considered unresolved and remain under debate.

RESULTS

Based on four loci and the most comprehensive taxonomic sampling analyzed to date, including currently recognized species, previously synonymized species, undescribed taxa, and poorly studied populations, this study produced a resolved phylogenetic framework that provided plausible species delimitation and an evolutionary time framework for the genus Ictalurus.

CONCLUSIONS

Our phylogenetic hypothesis revealed that Ictalurus comprises at least 13 evolutionary units, partially corroborating the current classification and identifying populations that emerge as putative undescribed taxa. The divergence times of the species indicate that the diversification of Ictalurus dates to the early Oligocene, confirming its status as one of the oldest genera within the family Ictaluridae.

Exceptional levels of species discovery ameliorate inferences of the biogeography and diversification of an Afrotropical catfish family

Endeavours in species discovery, particularly the characterisation of cryptic species, have been greatly aided by the application of DNA molecular sequence data to phylogenetic reconstruction and inference of evolutionary and biogeographic processes. However, the extent of cryptic and undescribed diversity remains unclear in tropical freshwaters, where biodiversity is declining at alarming rates. To investigate how data on previously undiscovered biodiversity impacts inferences of biogeography and diversification dynamics, we generated a densely sampled species-level family tree of Afrotropical Mochokidae catfishes (220 valid species) that was ca. 70 % complete. This was achieved through extensive continental sampling specifically targeting the genus Chiloglanis a specialist of the relatively unexplored fast-flowing lotic habitat. Applying multiple species-delimitation methods, we report exceptional levels of species discovery for a vertebrate genus, conservatively delimiting a staggering ca. 50 putative new Chiloglanis species, resulting in a near 80 % increase in species richness for the genus. Biogeographic reconstructions of the family identified the Congo Basin as a critical region in the generation of mochokid diversity, and further revealed complex scenarios for the build-up of continental assemblages of the two most species rich mochokid genera, Synodontis and Chiloglanis. While Syndontis showed most divergence events within freshwater ecoregions consistent with largely in situ diversification, Chiloglanis showed much less aggregation of freshwater ecoregions, suggesting dispersal as a key diversification process in this older group. Despite the significant increase in mochokid diversity identified here, diversification rates were best supported by a constant rate model consistent with patterns in many other tropical continental radiations. While our findings highlight fast-flowing lotic freshwaters as potential hotspots for undescribed and cryptic species diversity, a third of all freshwater fishes are currently threatened with extinction, signifying an urgent need to increase exploration of tropical freshwaters to better characterise and conserve its biodiversity.

Phylogenetic relationships of the North American catfishes (Ictaluridae, Siluriformes): Investigating the origins and parallel evolution of the troglobitic species

Insular habitats have played an important role in developing evolutionary theory, including natural selection and island biogeography. Caves are insular habitats that place extreme selective pressures on organisms due to the absence of light and food scarcity. Therefore, cave organisms present an excellent opportunity for studying colonization and speciation in response to the unique abiotic conditions that require extreme adaptations. One vertebrate family, the North American catfishes (Ictaluridae), includes four troglobitic species that inhabit the karst region bordering the western Gulf of Mexico. The phylogenetic relationships of these species have been contentious, and conflicting hypotheses have been proposed to explain their origins. The purpose of our study was to construct a time-calibrated phylogeny of Ictaluridae using first-occurrence fossil data and the largest molecular dataset on the group to date. We test the hypothesis that troglobitic ictalurids have evolved in parallel, thus resulting from repeated cave colonization events. We found that Prietella lundbergi is sister to surface-dwelling Ictalurus and that Prietella phreatophila + Trogloglanis pattersoni are sister to surface-dwelling Ameiurus, suggesting that ictalurids colonized subterranean habitats at least twice in evolutionary history. The sister relationship between Prietella phreatophila and Trogloglanis pattersoni may indicate that these two species diverged from a common ancestor following a subterranean dispersal event between Texas and Coahuila aquifers. We recovered Prietella as a polyphyletic genus and recommend P. lundbergi be removed from this genus. With respect to Ameiurus, we found evidence for a potentially undescribed species sister to A. platycephalus, which warrants further investigation of Atlantic and Gulf slope Ameiurus species. In Ictalurus, we identified shallow divergence between I. dugesii and I. ochoterenai, I. australis and I. mexicanus, and I. furcatus and I. meridionalis, indicating a need to reexamine the validity of each species. Lastly, we propose minor revisions to the intrageneric classification of Noturus including the restriction of subgenus Schilbeodes to N. gyrinus (type species), N. lachneri, N. leptacanthus, and N. nocturnus.

Mitochondrial DNA Corroborates the Genetic Variability of Clarias Catfishes (Siluriformes, Clariidae) from Cameroon

The airbreathing walking catfish (Clariidae: Clarias) comprises 32 species that are endemic to African freshwater systems. The species-level identification of this group is challenging due to their complex taxonomy and polymorphism. Prior to this study, the biological and ecological studies were restricted to a single species, Clarias gariepinus, resulting in a biased view of their genetic diversity in African waters. Here, we generated the 63-mitochondrial Cytochrome c oxidase subunit 1 (COI) gene sequences of Clarias camerunensis and Clarias gariepinus from the Nyong River in Cameroon. Both C. camerunensis and C. gariepinus species maintained adequate intra-species (2.7% and 2.31%) and inter-species (6.9% to 16.8% and 11.4% to 15.1%) genetic distances with other Clarias congeners distributed in African and Asian/Southeast Asian drainages. The mtCOI sequences revealed 13 and 20 unique haplotypes of C. camerunensis and C. gariepinus, respectively. The TCS networks revealed distinct haplotypes of C. camerunensis and shared haplotypes of C. gariepinus in African waters. The multiple species delimitation approaches (ABGD and PTP) revealed a total of 20 and 22 molecular operational taxonomic units (MOTUs), respectively. Among the two Clarias species examined, we found more than one MOTU in C. camerunensis, which is consistent with population structure and tree topology results. The phylogeny generated through Bayesian Inference analysis clearly separated C. camerunensis and C. gariepinus from other Clarias species with high posterior probability supports. The present study elucidates the occurrence of possible cryptic diversity and allopatric speciation of C. camerunensis in African drainages. Further, the present study confirms the reduced genetic diversity of C. gariepinus across its native and introduced range, which might have been induced by unscientific aquaculture practices. The study recommends a similar approach to the same and related species from different river basins to illuminate the true diversity of Clarias species in Africa and other countries.

Mitochondriomics of Clarias Fishes (Siluriformes: Clariidae) with a New Assembly of Clarias camerunensis: Insights into the Genetic Characterization and Diversification

The mitogenome of an endemic catfish Clarias camerunensis was determined from the Cameroon water. This circular mitogenome was 16,511 bp in length and comprised 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a single AT-rich control region. The heavy strand accommodates 28 genes, whereas the light strand is constituted by ND6 and eight transfer RNA (tRNA) genes. The C. camerunensis mitochondrial genome is AT biased (56.89%), as showcased in other Clarias species. The comparative analyses revealed that most of the Clarias species have 6 overlapping and 11 intergenic spacer regions. Most of the PCGs were initiated and terminated with the ATG start codon and TAA stop codon, respectively. The tRNAs of C. camerunensis folded into the distinctive cloverleaf secondary structure, except trnS1. The placement of the conserved domains in the control region was similar in all the Clarias species with highly variable nucleotides in CSB-I. Both maximum likelihood and Bayesian-based matrilineal phylogenies distinctly separated all Clarias species into five clades on the basis of their known distributions (South China, Sundaland, Indochina, India, and Africa). The TimeTree analysis revealed that the two major clades (Indo-Africa and Asia) of Clarias species might have diverged during the Paleogene (≈28.66 MYA). Our findings revealed the separation of Indian species (C. dussumieri) and African species (C. camerunensis and Clarias gariepinus) took place during the Paleogene, as well as the South Chinese species (Clarias fuscus) and Sundaland species (Clarias batrachus) splits from the Indochinese species (Clarias macrocephalus) during the Neogene through independent colonization. This pattern of biotic relationships highlights the influence of topography and geological events in determining the evolutionary history of Clarias species. The enrichment of mitogenomic data and multiple nuclear loci from their native range or type locality will confirm the true diversification of Clarias species in African and Asian countries.

The Evolution of Feeding Mechanics in the Danioninae, or Why Giant Danios Don't Suck Like Zebrafish

By linking anatomical structure to mechanical performance we can improve our understanding of how selection shapes morphology. Here we examined the functional morphology of feeding in fishes of the subfamily Danioninae (order Cypriniformes) to determine aspects of cranial evolution connected with their trophic diversification. The Danioninae comprise three major lineages and each employs a different feeding strategy. We gathered data on skull form and function from species in each clade, then assessed their evolutionary dynamics using phylogenetic-comparative methods. Differences between clades are strongly associated with differences in jaw protrusion. The paedomorphic Danionella clade does not use jaw protrusion at all, members of the Danio clade use jaw protrusion for suction production and prey capture, and members of the sister clade to Danio (e.g., Devario and Microdevario) use jaw protrusion to retain prey after capture. The shape of the premaxillary bone is a major determinant of protrusion ability, and premaxilla morphology in each of these lineages is consistent with their protrusion strategies. Premaxilla shapes have evolved rapidly, which indicates that they have been subjected to strong selection. We compared premaxilla development in giant danio (Devario aequipinnatus) and zebrafish (Danio rerio) and discuss a developmental mechanism that could shift danionine fishes between the feeding strategies employed by these species and their respective clades. We also identified a highly integrated evolutionary module that has been an important factor in the evolution of trophic mechanics within the Danioninae.

Insights into chromosomal evolution and sex determination of Pseudobagrus ussuriensis (Bagridae, Siluriformes) based on a chromosome-level genome

Pseudobagrus ussuriensis is an aquaculture catfish with significant sexual dimorphism. In this study, a chromosome-level genome with a size of 741.97 Mb was assembled for female P. ussuriensis. A total of 26 chromosome-level contigs covering 97.34% of the whole-genome assembly were obtained with an N50 of 28.53 Mb and an L50 of 11. A total of 24,075 protein-coding genes were identified, with 91.54% (22,039) genes being functionally annotated. Based on the genome assembly, four chromosome evolution clusters of catfishes were identified and the formation process of P. ussuriensis chromosomes was predicted. A total of 55 sex-related quantitative trait loci (QTLs) with a phenotypic variance explained value of 100% were located on chromosome 8 (chr08). The QTLs and other previously identified sex-specific markers were located in a sex-determining region of 16.83 Mb (from 6.90 to 23.73 Mb) on chr08, which was predicted as the X chromosome. The sex-determining region comprised 554 genes, with 135 of which being differently expressed between males and females/pseudofemales, and 16 candidate sex-determining genes were screened out. The results of this study provided a useful chromosome-level genome for genetic, genomic and evolutionary studies of P. ussuriensis, and also be useful for further studies on sex-determination mechanism analysis and sex-control breeding of this fish.

New phylogenetic insights into the African catfish families Mochokidae and Austroglanididae

Several hundred catfish species (order: Siluriformes) belonging to 11 families inhabit Africa, of which at least six families are endemic to the continent. Although four of those families are well-known to belong to the 'Big-Africa clade', no previous study has addressed the phylogenetic placement of the endemic African catfish family Austroglanididae in a comprehensive framework with molecular data. Furthermore, interrelationships within the 'Big-Africa clade', including the most diverse family Mochokidae, remain unclear. This study was therefore designed to help reconstruct inter- and intrarelationships of all currently valid mochokid genera, to infer their position within the 'Big Africa clade' and to establish a first molecular phylogenetic hypothesis of the relationships of the enigmatic Austroglanididae within the Siluriformes. We assembled a comprehensive mitogenomic dataset comprising all protein coding genes and representing almost all recognized catfish families (N = 33 of 39) with carefully selected species (N = 239). We recovered the monophyly of the previously identified multifamily clades 'Big Asia' and 'Big Africa' and determined Austroglanididae to be closely related to Pangasiidae, Ictaluroidea and Ariidae. Mochokidae was recovered as the sister group to a clade encompassing Auchenoglanididae, Claroteidae, Malapteruridae and the African Schilbeidae, albeit with low statistical support. The two mochokid subfamilies Mochokinae and Chiloglanidinae as well as the chiloglanid tribe Atopochilini were recovered as reciprocally monophyletic. The genus Acanthocleithron forms the sister group of all remaining Mochokinae, although with low support. The genus Atopodontus is the sister group of all remaining Atopochilini. In contrast to morphological reconstructions, the monophyly of the genus Chiloglanis was strongly supported in our analysis, with Chiloglanis macropterus nested within a Chiloglanis sublineage encompassing only other taxa from the Congo drainage. This is an important result because the phylogenetic relationships of C. macropterus have been controversial in the past, and because we and other researchers assumed that this species would be resolved as sister to most or all other members of Chiloglanis. The apparent paraphyly of Synodontis with respect to Microsynodontis provided an additional surprise, with Synodontis punu turning out to be the sister group of the latter genus.

An ancient truncated duplication of the anti‐Mullerian hormone receptor type 2 gene is a potential conserved master sex determinant in the Pangasiidae catfish family

The evolution of sex determination (SD) in teleosts is amazingly dynamic, as reflected by the variety of different master sex-determining genes identified. Pangasiids are economically important catfishes in South Asian countries, but little is known about their SD system. Here, we generated novel genomic resources for 12 Pangasiids and characterized their SD system. Based on a Pangasianodon hypophthalmus chromosome-scale genome assembly, we identified an anti-Müllerian hormone receptor type Ⅱ gene (amhr2) duplication, which was further characterized as being sex-linked in males and expressed only in testes. These results point to a Y chromosome male-specific duplication (amhr2by) of the autosomal amhr2a. Sequence annotation revealed that the P. hypophthalmus Amhr2by is truncated in its N-terminal domain, lacking the cysteine-rich extracellular part of the receptor that is crucial for ligand binding, suggesting a potential route for its neofunctionalization. Reference-guided assembly of 11 additional Pangasiids, along with sex-linkage studies, revealed that this truncated amhr2by duplication is a male-specific conserved gene in Pangasiids. Reconstructions of the amhr2 phylogeny suggested that amhr2by arose from an ancient duplication/insertion event at the root of the Siluroidei radiation that is dated to ~100 million years ago. Together these results bring multiple lines of evidence supporting that amhr2by is an ancient and conserved master sex-determining gene in Pangasiids, a finding that highlights the recurrent use of the transforming growth factor β pathway, which is often used for the recruitment of teleost master SD genes, and provides another empirical case towards firther understanding of dynamics of SD systems.

Limitations of phylogenomic data can drive inferred speciation rate shifts

Biodiversity analyses of phylogenomic timetrees have produced many high-profile examples of shifts in the rate of speciation across the tree of life. Temporally correlated events in ecology, climate, and biogeography are frequently invoked to explain these rate shifts. In a re-examination of 15 genomic timetrees and 25 major published studies of the pattern of speciation through time, we observed an unexpected correlation between the timing of reported rate shifts and the information content of sequence alignments. Here, we show that the paucity of sequence variation and insufficient species sampling in phylogenomic datasets are the likely drivers of many inferred speciation rate shifts, rather than the proposed biological explanations. Therefore, data limitations can produce predictable but spurious signals of rate shifts even when speciation rates may be similar across taxa and time. Our results suggest that the reliable detection of speciation rate shifts requires the acquisition and assembly of long phylogenomic alignments with near-complete species sampling and accurate estimates of species richness for the clades of study.

A chromosome-level genome assembly of the striped catfish (Pangasianodon hypophthalmus)

Striped catfish (Pangasianodon hypophthalmus), belonging to the Pangasiidae family, has become an economically important fish with wide cultivation in Southeast Asia. Owing to the high-fat trait, it is always considered as an oily fish. In our present study, a high-quality genome assembly of the striped catfish was generated by integration of Illumina short reads, Nanopore long reads and Hi-C data. A 731.7-Mb genome assembly was finally obtained, with a contig N50 of 3.5 Mb, a scaffold N50 of 29.5 Mb, and anchoring of 98.46% of the assembly onto 30 pseudochromosomes. The genome contained 36.9% repeat sequences, and a total 18,895 protein-coding genes were predicted. Interestingly, we identified a tandem triplication of fatty acid binding protein 1 gene (fabp1; thereby named as fabp1-1, fabp1-2 and fabp1-3 respectively), which may be related to the high fat content in striped catfish. Meanwhile, the FABP1-2 and -3 isoforms differed from FABP1-1 by several missense mutations including R126T, which may affect the fatty acid binding properties. In summary, we report a high-quality chromosome-level genome assembly of the striped catfish, which provides a valuable genetic resource for biomedical studies on the high-fat trait, and lays a solid foundation for practical aquaculture and molecular breeding of this international teleost species.

Sequencing an F1 hybrid of Silurus asotus and S. meridionalis enabled the assembly of high-quality parental genomes

Genome complexity such as heterozygosity may heavily influence its de novo assembly. Sequencing somatic cells of the F1 hybrids harboring two sets of genetic materials from both of the paternal and maternal species may avoid alleles discrimination during assembly. However, the feasibility of this strategy needs further assessments. We sequenced and assembled the genome of an F1 hybrid between Silurus asotus and S. meridionalis using the SequelII platform and Hi-C scaffolding technologies. More than 300 Gb raw data were generated, and the final assembly obtained 2344 scaffolds composed of 3017 contigs. The N50 length of scaffolds and contigs was 28.55 Mb and 7.49 Mb, respectively. Based on the mapping results of short reads generated for the paternal and maternal species, each of the 29 chromosomes originating from S. asotus and S. meridionalis was recognized. We recovered nearly 94% and 96% of the total length of S. asotus and S. meridionalis. BUSCO assessments and mapping analyses suggested that both genomes had high completeness and accuracy. Further analyses demonstrated the high collinearity between S. asotus, S. meridionalis, and the related Pelteobagrus fulvidraco. Comparison of the two genomes with that assembled only using the short reads from non-hybrid parental species detected a small portion of sequences that may be incorrectly assigned to the different species. We supposed that at least part of these situations may have resulted from mitotic recombination. The strategy of sequencing the F1 hybrid genome can recover the vast majority of the parental genomes and may improve the assembly of complex genomes.

Characterization of the complete mitochondrial genome of Macrotocinclus affinis (Siluriformes; Loricariidae) and phylogenetic studies of Siluriformes

The vertebrate mitochondrial genome is typically circular molecules made up of 14,000 to 16,000 bp, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (12 s rRNA and 16 s rRNA) and a control region. Compared with nuclear DNA, mitochondrial DNA has a higher mutation rate, so it is one of the most effective and reliable molecular markers in fish phylogeny. Macrotocinclus affinis was the only species in Macrotocinclus (it was classified as Otocinclus in the past) and currently lacks genetic information. Most of the current researches are based on the mitochondrial Cytb gene and RAG1 and RAG2 nuclear genes to study the phylogenetic analysis of Siluriformes. So, the study provides the characteristic features of the Macrotocinclus affinis mitochondrial genome and this is the first time that the phylogenetic relationship of Siluriformes has been reconstructed based on COI. We aimed to sequence the entire mitochondrial genome of Macrotocinclus affinis using conventional PCR techniques and to clarify its phylogenetic status in Siluriformes by using the COI sequence of mitochondria. In this study, we sequenced the whole mitochondrial genome of this species yielding a 16,632 bp circular assembly composed of the typical vertebrate mitochondrial features. It contains 13 protein-coding genes, two rRNA genes, 22 tRNA genes, a putative control region, and one origin of replication on the light-strand. The overall base composition includes A (30.07%), T (24.43%), C (29.43%) and G (16.01%). The genome composition is A + T biased (54.5%), and exhibits AT-skew (0.1036) and GC-skew (-0.2962). Moreover, the 13 PCGs encode 3850 amino acids in total. The result of the phylogenetic tree supports Macrotocinclus affinis has a closest relationship with Otocinclus cf. hoppei far. These results will help to understand the characteristics of the mitochondrial genome of Macrotocinclus affinis and provide molecular basis for the evolutionary relationship of Loricariidae.

Tapeworms (Cestoda) of Ictalurid Catfishes (Siluriformes) in North America: Redescription of Type Species of Two Genera and Proposal of Essexiellinae n. Subfam

Ictalurid catfishes (Siluriformes) in North America harbor proteocephalid tapeworms of the subfamily Corallobothriinae. Type species of 2 of 3 genera of these tapeworms from ictalurids are redescribed, based on museum and newly collected material. Essexiella fimbriata (Essex, 1928) is typified mainly by a wide, umbrella-shaped scolex with a metascolex formed by numerous folds of tissue, anteriorly directed suckers without sphincters, vitellarium bent inwards posteriorly, "flower-shaped" uterus (with anterior, lateral, and posterior diverticula), and a conspicuously pre-equatorial genital atrium. Verified host records of this cestode are only from 3 species of Ictalurus Rafinesque, 1820. Megathylacoides giganteum (Essex, 1928), which seems to be specific to the channel catfish (Ictalurus punctatus), possesses a globular scolex, with a weakly developed metascolex formed by tissue folds posterior to the suckers, anterolaterally directed suckers with large semilunar sphincters, proglottids that are widest at the level of the genital atrium at the anterior third of the proglottid, and uterine diverticula that do not reach the vitelline follicles laterally. A new subfamily, Essexiellinae Scholz and Barčák, is proposed to accommodate species of EssexiellaScholz, de Chambrier, Mariaux and Kuchta, 2011 (type genus), MegathylacoidesJones, Kerley and Sneed, 1956, and CorallotaeniaFreze, 1965 from ictalurid catfishes in the Nearctic Region. These tapeworms possess a metascolex, medullary genital organs, uterus lined with numerous chromophilic cells, pre-equatorial genital atrium, and uterine development of type 2. The new subfamily was monophyletic in all molecular phylogenetic analyses, being most closely related to 3 Neotropical proteocephalids from the redtail catfish, Phractocephalus hemioliopterus (Bloch and Schneider, 1801), but distant from all remaining proteocephalid tapeworms from freshwater fishes in North America.

Chromosomal Evolution in Aspredinidae (Teleostei, Siluriformes): Insights on Intra- and Interspecific Relationships with Related Groups

The family Aspredinidae comprises a clade of complex systematic relationships, both from molecular and morphological approaches. In this study, conventional and molecular cytogenetic studies coupled with nucleotide sequencing were performed in 6 Aspredininae species (Amaralia hypsiura, Bunocephalus cf. aloikae, Bunocephalus amaurus, Bunocephalus aff. coracoideus, Bunocephalus verrucosus, and Platystacus cotylephorus) from different locations of the Amazon hydrographic basin. Our results showed highly divergent diploid numbers (2n) among the species, ranging from 49 to 74, including the occurrence of an XX/X0 sex chromosome system. A neighbor-joining phylogram based on the cytochrome c oxidase I (COI) showed that Bunocephalus coracoideus is not a monophyletic clade, but closely related to B. verrucosus. The karyotypic data associated with COI suggest an ancestral karyotype for Aspredinidae with a reduced 2n, composed of bi-armed chromosomes and a trend toward chromosomal fissions resulting in higher diploid number karyotypes, mainly composed of acrocentric chromosomes. Evolutionary relationships were discussed under a phylogenetic context with related species from different Siluriformes families. The karyotype features and chromosomal diversity of Aspredinidae show an amazing differentiation, making this family a remarkable model for investigating the evolutionary dynamics in siluriforms as well as in fish as a whole.

Comparative chromosomal mapping of microsatellite repeats reveals divergent patterns of accumulation in 12 Siluridae (Teleostei: Siluriformes) species

The freshwater family Siluridae occurs in Eurasia and is especially speciose in South and Southeast Asia, representing an important aquaculture and fishery targets. However, despite the restricted cytogenetic data, a high diploid number variation (from 2n=40 to 92) characterizes this fish group. Considering the large genomic divergence among its species, silurid genomes have experienced an enormous diversification throughout their evolutionary history. Here, we aim to investigate the chromosomal distribution of several microsatellite repeats in 12 Siluridae species and infer about their possible roles in the karyotype evolution that occurred in this group. Our results indicate divergent patterns of microsatellite distribution and accumulation among the analyzed species. Indeed, they are especially present in significant chromosome locations, such as the centromeric and telomeric regions, precisely the ones associated with several kinds of chromosomal rearrangements. Our data provide pieces of evidence that repetitive DNAs played a direct role in fostering the chromosomal differentiation and biodiversity in this fish family.

Draft Genome Assembly of the Freshwater Apex Predator Wels Catfish (Silurus glanis) Using Linked-Read Sequencing

The wels catfish (Silurus glanis) is one of the largest freshwater fish species in the world. This top predator plays a key role in ecosystem stability, and represents an iconic trophy-fish for recreational fishermen. S. glanis is also a highly valued species for its high-quality boneless flesh, and has been cultivated for over 100 years in Eastern and Central Europe. The interest in rearing S. glanis continues to grow; the aquaculture production of this species has almost doubled during the last decade. However, despite its high ecological, cultural and economic importance, the available genomic resources for S. glanis are very limited. To fulfill this gap we report a de novo assembly and annotation of the whole genome sequence of a female S. glanis. The linked-read based technology with 10X Genomics Chromium chemistry and Supernova assembler produced a highly continuous draft genome of S. glanis: ∼0.8Gb assembly (scaffold N₅₀ = 3.2 Mb; longest individual scaffold = 13.9 Mb; BUSCO completeness = 84.2%), which included 313.3 Mb of putative repeated sequences. In total, 21,316 protein-coding genes were predicted, of which 96% were annotated functionally from either sequence homology or protein signature searches. The highly continuous genome assembly will be an invaluable resource for aquaculture genomics, genetics, conservation, and breeding research of S. glanis.

Helminth eggs from early cretaceous faeces

The exceptional fossil site of Las Hoyas (upper Barremian, Cuenca, Spain) yields abundant small to medium vertebrate coprolites, hindering the search for parasites. We studied the contents of 29 coprolites that were previously classified into distinct morphotypes. Several parasitic eggs were retrieved from two of these coprolites, confirming the second record of digenea trematode eggs and nematode (ascaridid) eggs from an Early Cretaceous locality. The cylindrical coprolite containing anisakid eggs was likely produced by a crocodylomorph as the parasite host, whereas the bump-headed lace coprolite indicates the role of a fish as an intermediary or definitive host of the trematodes and ascaridids. These trace and body fossils show that the Las Hoyas 126-129 Ma lacustrine ecosystem documents the early connection between basal Gonorynchiformes fish and digenetic trematodes.

Complete mitochondrial genome and assembled DNA barcoding analysis of Lutjanus fulgens (Valenciennes, 1830) and its comparison with other Lutjanus species

Lutjanus fulgens (Valenciennes, 1830) is a teleost species classified under the family Lutjanidae which is a native of the Eastern Atlantic Ocean. Though highly commercialized due to its abundance and good taste, the production output has declined in recent years. This is an indication of the need for effective management and conservation measures. However, accurate species identification will ensure strategic management and conservation measure. DNA-based species identification has proven its reliability in this regard via precise species identification. Several researchers have confirmed the accuracy of DNAbarcode as a species identification tool as well as species phylogeny analysis based on both the complete mitogenome and COI gene. Currently, nine specimens of L. fulgens were sampled from Ghana and subjected to DNA-based analysis, namely, complete mitochondrial DNAand COI gene (DNA barcoding) analyses. The mitogenomic result revealed that L. fulgens is made up of a 16,500 base pairs (bp) mtDNA which consists of 22 transfer RNAs, 13 protein-coding genes, and two ribosomal RNAs (GenBank Accession Number: MN398650). Furthermore, a sequence polymorphism analysis of the COIgene (MN986442-MN986450) detected two haplotypes. These haplotypes were both collected from the same fish landing site which suggests a possible cryptic linage diversity in the L. fulgens population at Vodza. According to the phylogeny examination, a close taxonomic relationship exists between L. fulgens and Lutjanus buccanella caused by a recent evolution termed as sympatric speciation. This study serves as a novel study for this species, building the foundation for future molecular-based study for this species and as a DNA barcode reference data.

Trinigyrus spp. (Monogenea: Dactylogyridae) from Brazilian catfishes: new species, molecular data and new morphological contributions to the genus

This study describes two new species, Trinigyrus anthus n. sp. and Trinigyrus carvalhoi n. sp., from gills of Hypostomus spp. from the Upper Paraná River basin, Brazil. Trinigyrus peregrinus is redescribed based on examination of its holotype, paratypes and new material of specimens parasitizing Pterygoplichthys ambrosettii, also from the Upper Paraná River basin, Brazil. New morphological features were included in the diagnosis of the genus, such as the presence of a sclerotized border on the anchor base, and a weakly sclerotized fringe on the base of the male copulatory organ (MCO). Trinigyrus anthus n. sp. differs from other congeners by the shape of the MCO, presenting an enlarged base with sclerotized fringes resembling flower petals. Trinigyrus carvalhoi n. sp. and T. peregrinus are similar but can be differentiated from each other mainly by the sclerotization of the vagina (absent in the new species), and the morphology of the MCO (C-shaped versus one counterclockwise circle, respectively). For the first time, gene sequences of Trinigyrus spp. from Brazil were obtained (partial ribosomal 28S and mitochondrial cytochrome c oxidase I (mtCOI)). The genetic divergences among the new species and T. peregrinus varied from 2 to 3% (6‒18 pb) based on sequences of 28S ribosomal DNA (rDNA), and 6-7% (83‒92 pb) using mtCOI. Phylogenetic analyses based on partial 28S rDNA revealed that Trinigyrus, Heteropriapulus and Unilatus formed a monophyletic and well-supported clade of monogeneans from Neotropical freshwater loricariids, suggesting a closer relationship among these dactylogyrids and their hosts.

Origins of Afrotropical freshwater fishes

The Afrotropics house a diverse freshwater ichthyofauna with > 3000 species, almost all of which are endemic. Recent progress in dated phylogenetics and palaeontology of several groups of Afrotropical freshwater fishes (AFFs) has allowed the testing of palaeoecology- and palaeogeography-based hypotheses explaining their early presence in Africa. Seven hypotheses were tested for 37 most-inclusive monophyletic groups of AFFs. Results indicated that ten lineages originated from direct, but asynchronous, marine-to-freshwater shifts. These lineages contribute < 2% to the current AFF species richness. Eleven lineages colonized the Afrotropics from the Orient after the Afro-Arabian plate collided with Eurasia in the early Oligocene. These lineages contribute ~20% to the total diversity. There are seven sister relationships between Afrotropical and Neotropical taxa. For only three of them (4% of the species diversity), the continental drift vicariance hypothesis was not rejected. Distributions of the other four younger trans-Atlantic lineages are better explained by post-drifting long-distance dispersal. In those cases, I discuss the possibility of dispersal through the Northern Hemisphere as an alternative to direct trans-Atlantic dispersal. The origins of ten AFF lineages, including the most species-rich Pseudocrenilabrinae (> 1100 species), are not yet established with confidence.

Genomic Organization of Repetitive DNA Elements and Extensive Karyotype Diversity of Silurid Catfishes (Teleostei: Siluriformes): A Comparative Cytogenetic Approach

The catfish family Siluridae contains 107 described species distributed in Asia, but with some distributed in Europe. In this study, karyotypes and other chromosomal characteristics of 15 species from eight genera were examined using conventional and molecular cytogenetic protocols. Our results showed the diploid number (2n) to be highly divergent among species, ranging from 2n = 40 to 92, with the modal frequency comprising 56 to 64 chromosomes. Accordingly, the ratio of uni- and bi-armed chromosomes is also highly variable, thus suggesting extensive chromosomal rearrangements. Only one chromosome pair bearing major rDNA sites occurs in most species, except for Wallago micropogon, Ompok siluroides, and Kryptoterus giminus with two; and Silurichthys phaiosoma with five such pairs. In contrast, chromosomes bearing 5S rDNA sites range from one to as high as nine pairs among the species. Comparative genomic hybridization (CGH) experiments evidenced large genomic divergence, even between congeneric species. As a whole, we conclude that karyotype features and chromosomal diversity of the silurid catfishes are unusually extensive, but parallel some other catfish lineages and primary freshwater fish groups, thus making silurids an important model for investigating the evolutionary dynamics of fish chromosomes.

A draft genome of the striped catfish, Pangasianodon hypophthalmus, for comparative analysis of genes relevant to development and a resource for aquaculture improvement

Background

The striped catfish, Pangasianodon hypophthalmus, is a freshwater and benthopelagic fish common in the Mekong River delta. Catfish constitute a valuable source of dietary protein. Therefore, they are cultured worldwide, and P. hypophthalmus is a food staple in the Mekong area. However, genetic information about the culture stock, is unavailable for breeding improvement, although genetics of the channel catfish, Ictalurus punctatus, has been reported. To acquire genome sequence data as a useful resource for marker-assisted breeding, we decoded a draft genome of P. hypophthalmus and performed comparative analyses.

Results

Using the Illumina platform, we obtained both nuclear and mitochondrial DNA sequences. Molecular phylogeny using the mitochondrial genome confirmed that P. hypophthalmus is a member of the family Pangasiidae and is nested within a clade including the families Cranoglanididae and Ictaluridae. The nuclear genome was estimated at approximately 700 Mb, assembled into 568 scaffolds with an N50 of 14.29 Mbp, and was estimated to contain ~ 28,600 protein-coding genes, comparable to those of channel catfish and zebrafish. Interestingly, zebrafish produce gadusol, but genes for biosynthesis of this sunscreen compound have been lost from catfish genomes. The differences in gene contents between these two catfishes were found in genes for vitamin D-binding protein and cytosolic phospholipase A₂, which have lost only in channel catfish. The Hox cluster in catfish genomes comprised seven paralogous groups, similar to that of zebrafish, and comparative analysis clarified catfish lineage-specific losses of A5a, B10a, and A11a. Genes for insulin-like growth factor (IGF) signaling were conserved between the two catfish genomes. In addition to identification of MHC class I and sex determination-related gene loci, the hypothetical chromosomes by comparison with the channel catfish demonstrated the usefulness of the striped catfish genome as a marker resource.

Conclusions

We developed genomic resources for the striped catfish. Possible conservation of genes for development and marker candidates were confirmed by comparing the assembled genome to that of a model fish, Danio rerio, and to channel catfish. Since the catfish genomic constituent resembles that of zebrafish, it is likely that zebrafish data for gene functions is applicable to striped catfish as well.

Electronic supplementary material

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

Mitochondrial genome of the yellow catfish Pelteobagrus fulvidraco and insights into Bagridae phylogenetics

The mitochondrial genome (mitogenome) can provide important information for understanding phylogenetic analysis and molecular evolution. Herein, we amplified the complete mitogenome sequence of Pelteobagrus fulvidraco. The mitogenome was 16,526 bp in length and included 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes and a non-coding control region (D-loop). Both the organization and location of genes in the mitogenome were consistent with those from Siluriformes fishes previously published in GenBank. The phylogenetic relationships based on Bayesian inference (BI) and Maximum likelihood (ML) methods showed that P. fulvidraco has close relationships with Pelteobagrus eupogon and Tachysurus intermedius, suggesting that P. fulvidraco belongs to Tachysurus. This study provides evidence that Tachysurus, Pseudobagrus and Leiocassis do not form monophyly, but that these three genera form a monophyletic group. Our results provide reference for further phylogenetic research of the Bagridae species.

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.].

Complete mitochondrial genome of near threatened butter Catfish Ompok bimaculatus (Siluriformes: Siluridae)

The complete mitochondrial genome of Ompok bimaculatus was obtained, using illumina high-throughput NextSeq 500 with 2 × 150 bp sequencing of mitochondrial DNA. The genome of O. bimaculatus was 16,482 bp in length (GenBank Accession No. KY887474) comprised of 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and a control region i.e. D-loop. In present mitogenome, 9 SSR were identified and validated in silico and secondary structures of all the 22 tRNA were predicted. The arrangement of genes was found identical to other siluriformes fish mitogenomes available in NCBI database. Phylogenetic relationship with closely related species were established which provide useful insights into taxonomic status of the species.

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.