Dataset for genome sequencing and de novo assembly of the Vietnamese bighead catfish (Clarias macrocephalus Günther, 1864)

The complete mitochondrial genome of the blackskin catfish (Clarias meladerma: Clariidae) from Rokan River, Riau, Indonesia

Clarias meladerma Bleeker, 1846, a native catfish species in Indonesia belonging to the family Clariidae. The present study the complete mitochondrial genome sequence of C. meladerma from the Rokan River was sequenced by using next-generation sequencing, and its phylogenetic relationship was analyzed. The mitochondrial genome comprises 13 protein-coding genes (PCGs), 22 tRNA genes, and two rRNA genes, with a total length of 16,808 bp. The mitogenome of C. meladerma exhibits a base composition of 32.49% adenine, 25.75% thymine, 14.51% guanine, and 27.25% cytosine. Phylogenetic analysis indicated that C. meladerma has the same clade with C. macrocephalus, C. batrachus, and C. fucus. In essence, the findings of this study lay down a genetic foundation for future investigations into C. meladerma.

The Utilization of Reference-Guided Assembly and In Silico Libraries Improves the Draft Genome of Clarias batrachus and Culter alburnus

Long-read sequencing technologies can generate highly contiguous genome assemblies compared to short-read methods. However, their higher cost often poses a significant barrier. To address this, we explore the utilization of mapping-based genome assembly and reference-guided assembly as cost-effective alternative approaches. We assess the efficacy of these approaches in improving the contiguity of Clarias batrachus and Culter alburnus draft genomes. Our findings demonstrate that employing an iterative mapping strategy leads to a reduction in assembly errors. Specifically, after three iterations, the Mismatches per 100 kbp value for the C. batrachus genome decreased from 2447.20 to 2432.67, reaching a minimum of 2422.67 after two iterations. Additionally, the N50 value for the C. batrachus genome increased from 362,143 to 1,315,126 bp, with a maximum of 1,315,403 bp after two iterations. Furthermore, we achieved Mismatches per 100 kbp values of 3.70 for the reference-guided assembly of C. batrachus and 0.34 for C. alburnus. Correspondingly, the N50 value for the C. batrachus and C. alburnus genomes increased from 362,143 bp and 3,686,385 bp to 2,026,888 bp and 43,735,735 bp, respectively. Finally, we successfully utilized the improved C. batrachus and C. alburnus genomes to compare genome studies using the combined approach of Ragout and Ragtag. Through a comprehensive comparative analysis of mapping-based and reference-guided genome assembly methods, we shed light on the specific contributions of reference-guided assembly in reducing assembly errors and improving assembly continuity and integrity. These advancements establish reference-guided assembly and the utilization of in silico libraries as a promising and suitable approach for comparative genomics studies.

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.

Chromosome-level assembly and annotation of the blue catfish Ictalurus furcatus, an aquaculture species for hybrid catfish reproduction, epigenetics, and heterosis studies

Background

The blue catfish is of great value in aquaculture and recreational fisheries. The F₁ hybrids of female channel catfish (Ictalurus punctatus) × male blue catfish (Ictalurusfurcatus) have been the primary driver of US catfish production in recent years because of superior growth, survival, and carcass yield. The channel–blue hybrid also provides an excellent model to investigate molecular mechanisms of environment-dependent heterosis. However, transcriptome and methylome studies suffered from low alignment rates to the channel catfish genome due to divergence, and the genome resources for blue catfish are not publicly available.

Results

The blue catfish genome assembly is 841.86 Mbp in length with excellent continuity (8.6 Mbp contig N50, 28.2 Mbp scaffold N50) and completeness (98.6% Eukaryota and 97.0% Actinopterygii BUSCO). A total of 30,971 protein-coding genes were predicted, of which 21,781 were supported by RNA sequencing evidence. Phylogenomic analyses revealed that it diverged from channel catfish approximately 9 million years ago with 15.7 million fixed nucleotide differences. The within-species single-nucleotide polymorphism (SNP) density is 0.32% between the most aquaculturally important blue catfish strains (D&B and Rio Grande). Gene family analysis discovered significant expansion of immune-related families in the blue catfish lineage, which may contribute to disease resistance in blue catfish.

Conclusions

We reported the first high-quality, chromosome-level assembly of the blue catfish genome, which provides the necessary genomic tool kit for transcriptome and methylome analysis, SNP discovery and marker-assisted selection, gene editing and genome engineering, and reproductive enhancement of the blue catfish and hybrid catfish.

Endogenic upregulations of HIF/VEGF signaling pathway genes promote air breathing organ angiogenesis in bimodal respiration fish

Air-breathing has evolved independently serval times with a variety of air-breathing organs (ABOs) in fish. The physiology of the air-breathing in bimodal respiration fish has been well understood, while studies on molecular mechanisms of the character are very limited. In the present study, we first determined the gill indexes of 110 fish species including 25 and 85 kinds of bimodal respiration fishes and non-air-breathing fishes, respectively. Then combined with histological observations of gills and ABOs/non-ABOs in three bimodal respiration fishes and two non-air breathing fishes, we found that the bimodal respiration fish was always of a degeneration gill and a well-vascularized ABO. Meanwhile, a comparative transcriptome analysis of posterior intestines, namely a well vascularized ABO in Misgurnus anguillicaudatus and a non-ABO in Leptobotia elongata, was performed to expound molecular variations of the air-breathing character. A total of 5,003 orthologous genes were identified. Among them, 1,189 orthologous genes were differentially expressed, which were enriched in 14 KEGG pathways. More specially, the expressions of hemoglobin genes and various HIF/VEGF signaling pathway genes were obviously upregulated in the ABO of M. anguillicaudatus. Moreover, we found that HIF-1α, VEGFAa, and MAP2K1 were co-expressed dramatically higher in ABOs of bimodal respiration fishes than those of non-ABOs of non-air-breathing fishes. These results indicated that the HIF/VEGF pathway played an important role in ABO angiogenesis/formation to promote fish to do aerial respiration. This study will contribute to our understanding of molecular mechanisms of air-breathing in fish.

Comparative Genomic and Transcriptomic Analyses Revealed Twenty-Six Candidate Genes Involved in the Air-Breathing Development and Function of the Bighead Catfish Clarias macrocephalus

The bighead catfish (Clarias macrocephalus) and channel catfish (Ictalurus punctatus) are freshwater species in the Siluriformes order. C. macrocephalus has both gills and modified gill structures serving as an air-breathing organ (ABO), while I. punctatus does not possess such an organ, and cannot breathe in air, providing an excellent model for studying the molecular basis of ABO development in teleost fish. To investigate the critical time window for the development of air-breathing function, seven development stages were selected based on hypoxia challenge results, and RNA-seq was performed upon C. macrocephalus to compare with the non-air-breathing I. punctatus. Five-hundred million reads were generated and 25,239 expressed genes were annotated in C. macrocephalus. Among those, 8675 genes were differentially expressed across developmental stages. Comparative genomic analysis identified 1458 C. macrocephalus specific genes, which were absent in I. punctatus. Gene network and protein-protein interaction analyses identified 26 key hub genes involved in the air-breathing function. Three top candidate genes, mb, ngb, hbae, are mainly associated with oxygen carrying, oxygen binding, and heme binding activities. Our study provides a rich data set for exploring the genomic basis of air-breathing function in C. macrocephalus and offers insights into the adaption to hypoxic environments.