Whole-Genome Sequencing of the Giant Devil Catfish, Bagarius yarrelli

Chromosome-level genome assembly of the largefin longbarbel catfish (Hemibagrus macropterus)

The largefin longbarbel catfish, Hemibagrus macropterus, is an economically important fish species in southwestern China, with males growing faster than females. This study presents a high-quality chromosome-level genome assembly of the largefin longbarbel catfish, generated by integrating Illumina short reads, PacBio HiFi long reads, and Hi-C data. The assembled genome size was 858.5 Mb, with a contig and scaffold N50 of 5.8 Mb and 28.4 Mb, respectively. A total of 656 contigs were successfully anchored to 30 pseudochromosomes with a BUSCO score of 97.7%, consistent with the number of chromosomes analyzed by karyotype. The genome contained 29.5% repeat sequences, and a predicted total of 26,613 protein-coding genes, of which 25,769 (96.8%) were functionally annotated in different databases. Evolutionary analysis showed that H. macropterus was most closely related to H. wyckioides, with a divergence time of approximately 16.3 million years. Chromosomal syntenic relationships among H. macropterus, H. wyckioides, and Pelteobagrus fulvidraco revealed a one-to-one relationship for most chromosomes, except for break, fission, and inversion of some chromosomes. The first high-quality reference genome will not only provide a valuable genetic resource for the study of sex determination mechanisms and genetic breeding of largefin longbarbel catfish, but also contribute to comparative analyses of genome and chromosome evolution within Siluriformes.

Enlarged fins of Tibetan catfish provide new evidence of adaptation to high plateau

The uplift of the Tibetan Plateau significantly altered the geomorphology and climate of the Euroasia by creating large mountains and rivers. Fishes are more likely to be affected relative to other organisms, as they are largely restricted to river systems. Faced with the rapidly flowing water in the Tibetan Plateau, a group of catfish has evolved greatly enlarged pectoral fins with more numbers of fin-rays to form an adhesive apparatus. However, the genetic basis of these adaptations in Tibetan catfishes remains elusive. In this study, we performed comparative genomic analyses based on the chromosome-level genome of Glyptosternum maculatum in family Sisoridae and detected some proteins with conspicuously high evolutionary rates in particular in genes involved in skeleton development, energy metabolism, and hypoxia response. We found that the hoxd12a gene evolved faster and a loss-of-function assay of hoxd12a supports a potential role for this gene in shaping the enlarged fins of these Tibetan catfishes. Other genes with amino acid replacements and signatures of positive selection included proteins involved in low temperature (TRMU) and hypoxia (VHL) responses. Functional assays reveal that the G. maculatumTRMU allele generates more mitochondrial ATP than the ancestral allele found in low-altitude fishes. Functional assays of VHL alleles suggest that the G. maculatum allele has lower transactivation activity than the low-altitude forms. These findings provide a window into the genomic underpinnings of physiological adaptations that permit G. maculatum to survive in the harsh environment of the Tibetan Himalayas that mirror those that are convergently found in other vertebrates such as humans.

Whole genome assembly of the armored loricariid catfish Ancistrus triradiatus highlights herbivory signatures

The catfish Ancistrus triradiatus belongs to the species-rich family Loricariidae. Loricariids display remarkable traits such as herbivory, a benthic lifestyle, the absence of scales but the presence of dermal bony plates. They are exported as ornamental fish worldwide, with escaped fishes becoming a threat locally. Although genetic and phylogenetic studies are continuously increasing and developmental genetic investigations are underway, no genome assembly has been formally proposed for Loricariidae yet. We report a high-quality genome assembly of Ancistrus triradiatus using long and short reads, and a newly assembled transcriptome. The genome assembly is composed of 9530 scaffolds, including 85.6% of ray-finned fish BUSCOs, and 26,885 predicted protein-coding genes. The genomic GC content is higher than in other catfishes, reflecting the higher metabolism associated with herbivory. The examination of the SCPP gene family indicates that the genes presumably triggering scale loss when absent, are present in the scaleless A. triradiatus, questioning their explanatory role. The analysis of the opsin gene repertoire revealed that gene losses associated to the nocturnal lifestyle of catfishes were not entirely found in A. triradiatus, as the UV-sensitive opsin 5 is present. Finally, most gene family expansions were related to immunity except the gamma crystallin gene family which controls pupil shape and sub-aquatic vision. Thus, the genome of A. triradiatus reveals that fish herbivory may be related to the photic zone habitat, conditions metabolism, photoreception and visual functions. This genome is the first for the catfish suborder Loricarioidei and will serve as backbone for future genetic, developmental and conservation studies.

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.

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 complete mitogenome datasets and their characterization of the European catfish (Silurus glanis)

We present new complete mitogenome sequences of Silurus glanis (S. glanis) from 4 samples such as male and female individuals from two countries (Hungary, Czech Republic). The complete mitochondria were determined from genome sequencing by using Illumina MiSeq platform resulting in long, 300 bp. paired-end reads. De novo assembly was performed resulting in one nod (scaffold) covering the total mitochondria in each sample. The mitochondrial genomes were circular, double-stranded molecules of 16,524 bp in length and consisted of 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes, and 1 control region. These sequences were deposited in the NCBI GeneBank under the accession numbers (MW796040, MW796041, MW796042, MW796043) and compared with the only available S. glanis mitochondrial genome (NC_014261.1) sequenced by unidentified technology and showed 99% similarity. We found in seq1 82, in seq2 82, seq3 83, seq4 82 nucleotide alterations involving 10 protein-coding genes and meaning 29 amino acid substitutions as well.

Genome-Wide Transcriptional Changes of Rhodosporidium kratochvilovae at Low Temperature

Rhodosporidium kratochvilovae strain YM25235 is a cold-adapted oleaginous yeast strain that can grow at 15°C. It is capable of producing polyunsaturated fatty acids. Here, we used the Nanopore Platform to first assemble the R. kratochvilovae strain YM25235 genome into a 23.71 Mb size containing 46 scaffolds and 8,472 predicted genes. To explore the molecular mechanism behind the low temperature response of R. kratochvilovae strain YM25235, we analyzed the RNA transcriptomic data from low temperature (15°C) and normal temperature (30°C) groups using the next-generation deep sequencing technology (RNA-seq). We identified 1,300 differentially expressed genes (DEGs) by comparing the cultures grown at low temperature (15°C) and normal temperature (30°C) transcriptome libraries, including 553 significantly upregulated and 747 significantly downregulated DEGs. Gene ontology and pathway enrichment analysis revealed that DEGs were primarily related to metabolic processes, cellular processes, cellular organelles, and catalytic activity, whereas the overrepresented pathways included the MAPK signaling pathway, metabolic pathways, and amino sugar and nucleotide sugar metabolism. We validated the RNA-seq results by detecting the expression of 15 DEGs using qPCR. This study provides valuable information on the low temperature response of R. kratochvilovae strain YM25235 for further research and broadens our understanding for the response of R. kratochvilovae strain YM25235 to low temperature.

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.

Chromosome-level genome assembly of the Chinese longsnout catfish Leiocassis longirostris

The Chinese longsnout catfish ( Leiocassis longirostris Günther) is one of the most economically important freshwater fish in China. As wild populations have declined sharply in recent years, it is also a valuable model for research on sexual dimorphism, comparative biology, and conservation. However, the current lack of high-quality chromosome-level genome information for the species hinders the advancement of comparative genomic analysis and evolutionary studies. Therefore, we constructed the first high-quality chromosome-level reference genome for L. longirostris. The total genome was 703.19 Mb, with 389 contigs and contig N50 length of 4.29 Mb. Using high-throughput chromosome conformation capture (Hi-C) data, the genome sequences (685.53 Mb) were scaffolded into 26 chromosomes ranging from 17.36 to 43.97 Mb, resulting in a chromosomal anchoring rate for the genome of 97.44%. In total, 23 708 protein-coding genes were identified in the genome. Phylogenetic analysis indicated that L. longirostris and its closest related species P. fulvidraco diverged approximately 26.6 million years ago. This high-quality reference genome of L. longirostris should pave the way for future genomic comparisons and evolutionary research.

The American Paddlefish Genome Provides Novel Insights into Chromosomal Evolution and Bone Mineralization in Early Vertebrates

Sturgeons and paddlefishes (Acipenseriformes) occupy the basal position of ray-finned fishes, although they have cartilaginous skeletons as in Chondrichthyes. This evolutionary status and their morphological specializations make them a research focus, but their complex genomes (polyploidy and the presence of microchromosomes) bring obstacles and challenges to molecular studies. Here, we generated the first high-quality genome assembly of the American paddlefish (Polyodon spathula) at a chromosome level. Comparative genomic analyses revealed a recent species-specific whole-genome duplication event, and extensive chromosomal changes, including head-to-head fusions of pairs of intact, large ancestral chromosomes within the paddlefish. We also provide an overview of the paddlefish SCPP (secretory calcium-binding phosphoprotein) repertoire that is responsible for tissue mineralization, demonstrating that the earliest flourishing of SCPP members occurred at least before the split between Acipenseriformes and teleosts. In summary, this genome assembly provides a genetic resource for understanding chromosomal evolution in polyploid nonteleost fishes and bone mineralization in early vertebrates.

The use of molecular markers in the verification of fish and seafood authenticity and the detection of adulteration

The verification of authenticity and detection of food mislabeling are elements that have been of high importance for centuries. During the last few decades there has been an increasing consumer demand for the verification of food identity and the implementation of stricter controls around these matters. Fish and seafood are among the most easily adulterated foodstuffs mainly due to the significant alterations of the species' morphological characteristics that occur during the different types of processing, which render the visual identification of the animals impossible. Even simple processes, such as filleting remove very important morphological elements and suffice to prevent the visual identification of species in marketed products. Novel techniques have therefore been developed that allow identification of species, the differentiation between species and also the differentiation of individuals that belong to the same species but grow in different populations and regions. Molecular markers have been used during the last few decades to fulfill this purpose and several improvements have been implemented rendering their use applicable to a commercial scale. The reliability, accuracy, reproducibility, and time-and cost-effectiveness of these techniques allowed them to be established as routine methods in the industry and research institutes. This review article aims at presenting the most important molecular markers used for the authentication of fish and seafood. The most important techniques are described, and the results of numerous studies are outlined and discussed, allowing interested parties to easily access and compare information about several techniques and fish/seafood species.

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.

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

Freshwater catfish of the genus Clarias, known as the airbreathing catfish, are widespread and important for food security through small scale inland fisheries and aquaculture. Limited genomic data are available for this important group of fishes. The bighead catfish (Clarias macrocephalus) is a commercial aquaculture species in southeast Asia used for aquaculture and threatened in its natural environment through habitat destruction, over-exploitation and competition from other introduced species of Clarias. Despite its commercial importance and threats to natural populations, public databases do not include any genomic data for C. macrocephalus. We present the first genomic data for the bighead catfish from Illumina sequencing. A total of 128 Gb of sequence data in paired-end 150 bp reads were assembled de novo, generating a final assembly of 883 Mbp contained in 27,833 scaffolds (N₅₀ length: 80.8 kbp) with BUSCO completeness assessments of 96.3% and 87.6% based on metazoan and Actinopterygii ortholog datasets, respectively. Annotation of the genome predicted 21,124 gene sequences, which were assigned putative functions based on homology to existing protein sequences in public databases. Raw fastq reads and the final version of the genome assembly have been deposited in the NCBI (BioProject: PRJNA604477, WGS: JAAGKR000000000, SRA: SRR11188453). The complete C. macrocephalus mitochondrial genome was also recovered from the same sequence read dataset and is available on NCBI (accession: MT109097), representing the first mitogenome for this species. Lastly, we find an expansion of the mb and ora1 genes thought to be associated with adaptations to air-breathing and a semi-terrestrial life style in this genus of catfish.

Whole genome sequencing and comparative genomic analysis of oleaginous red yeast Sporobolomyces pararoseus NGR identifies candidate genes for biotechnological potential and ballistospores-shooting

BACKGROUND

Sporobolomyces pararoseus is regarded as an oleaginous red yeast, which synthesizes numerous valuable compounds with wide industrial usages. This species hold biotechnological interests in biodiesel, food and cosmetics industries. Moreover, the ballistospores-shooting promotes the colonizing of S. pararoseus in most terrestrial and marine ecosystems. However, very little is known about the basic genomic features of S. pararoseus. To assess the biotechnological potential and ballistospores-shooting mechanism of S. pararoseus on genome-scale, the whole genome sequencing was performed by next-generation sequencing technology.

RESULTS

Here, we used Illumina Hiseq platform to firstly assemble S. pararoseus genome into 20.9 Mb containing 54 scaffolds and 5963 predicted genes with a N50 length of 2,038,020 bp and GC content of 47.59%. Genome completeness (BUSCO alignment: 95.4%) and RNA-seq analysis (expressed genes: 98.68%) indicated the high-quality features of the current genome. Through the annotation information of the genome, we screened many key genes involved in carotenoids, lipids, carbohydrate metabolism and signal transduction pathways. A phylogenetic assessment suggested that the evolutionary trajectory of the order Sporidiobolales species was evolved from genus Sporobolomyces to Rhodotorula through the mediator Rhodosporidiobolus. Compared to the lacking ballistospores Rhodotorula toruloides and Saccharomyces cerevisiae, we found genes enriched for spore germination and sugar metabolism. These genes might be responsible for the ballistospores-shooting in S. pararoseus NGR.

CONCLUSION

These results greatly advance our understanding of S. pararoseus NGR in biotechnological potential and ballistospores-shooting, which help further research of genetic manipulation, metabolic engineering as well as its evolutionary direction.