The genome of walking catfish Clarias magur (Hamilton, 1822) unveils the genetic basis that may have facilitated the development of environmental and terrestrial adaptation systems in air-breathing catfishes

Fish antifreeze protein origin in sculpins by frameshifting within a duplicated housekeeping gene

Antifreeze proteins (AFPs) are found in a variety of marine cold-water fishes where they prevent freezing by binding to nascent ice crystals. Their diversity (types I, II, III and antifreeze glycoproteins), as well as their scattered taxonomic distribution hint at their complex evolutionary history. In particular, type I AFPs appear to have arisen in response to the Late Cenozoic Ice Age that began ~ 34 million years ago via convergence in four different groups of fish that diverged from lineages lacking this AFP. The progenitor of the alanine-rich α-helical type I AFPs of sculpins has now been identified as lunapark, an integral membrane protein of the endoplasmic reticulum. Following gene duplication and loss of all but three of the 15 exons, the final exon, which encoded a glutamate- and glutamine-rich segment, was converted to an alanine-rich sequence by a combination of frameshifting and mutation. Subsequent gene duplications produced numerous isoforms falling into four distinct groups. The origin of the flounder type I AFP is quite different. Here, a small segment from the original antiviral protein gene was amplified and the rest of the coding sequence was lost, while the gene structure was largely retained. The independent origins of type I AFPs with up to 83% sequence identity in flounder and sculpin demonstrate strong convergent selection at the level of protein sequence for alanine-rich single alpha helices that bind to ice. Recent acquisition of these AFPs has allowed sculpins to occupy icy seawater niches with reduced competition and predation from other teleost species.

DNA Methylation Profiling in Genetically Selected Clarias magur (Hamilton, 1822) Provides Insights into the Epigenetic Regulation of Growth and Development

DNA methylation is an epigenetic alteration that impacts gene expression without changing the DNA sequence affecting an organism's phenotype. This study utilized a reduced representation bisulfite sequencing (RRBS) approach to investigate the patterns of DNA methylation in genetically selected Clarias magur stocks. RRBS generated 249.22 million reads, with an average of 490,120 methylation sites detected in various parts of genes, including exons, introns, and intergenic regions. A total of 896 differentially methylated regions (DMRs) were identified; 356 and 540 were detected as hyper-methylated and hypo-methylated regions, respectively. The DMRs and their association with overlapping genes were explored using whole genome data of magur, which revealed 205 genes in exonic, 210 in intronic, and 480 in intergenic regions. The analysis identified the maximum number of genes enriched in biological processes such as RNA biosynthetic process, response to growth factors, nervous system development, neurogenesis, and anatomical structure morphogenesis. Differentially methylated genes (DMGs) such as myrip, mylk3, mafb, egr3, ndnf, meis2a, foxn3, bmp1a, plxna3, fgf6, sipa1l1, mcu, cnot8, trim55b, and myof were associated with growth and development. The selected DMGs were analyzed using real-time PCR, which showed altered mRNA expression levels. This work offers insights into the epigenetic mechanisms governing growth performance regulation in magur stocks. This work provides a valuable resource of epigenetic data that could be integrated into breeding programs to select high-performing individuals.

A chromosome-level genome of the striated frogfish (Antennarius striatus)

The striated frogfish (Antennarius striatus), a member of the sub-order Antennarioidei within the order Lophiiformes, possesses remarkable adaptations. These include the ability to modulate body coloration for camouflage, utilize bioluminescent esca for predation, and employ elbow-like pectoral fins for terrestrial locomotion, making it a valuable model for studying bioluminescence, adaptive camouflage, fin-to-limb transition, and walking-like behaviors. To better study and contribute to the conservation of the striated frogfish, we obtained the micro-CT image of the pectoral fin bones and generated a high-quality, chromosome-level genome assembly using multiple sequencing technologies. The assembly spans 548.56 Mb with a contig N50 of 21.05 Mb, and 99.35% of the genome is anchored on 24 chromosomes, making it the most complete genome available within Lophiiformes. The genome annotation revealed 28.43% repetitive sequences and 23,945 protein-coding genes. This chromosome-level genome provides valuable genetic resources for frogfish conservation and offers insights into the genetic mechanisms underlying its unique phenotypic evolution. Furthermore, it establishes a foundation for future research on limb development and adaptive camouflage in this species.

The concept of "state fishes of India": Current status, knowledge gaps, and strategic plans for conservation and sustainable utilization

The concept of "state fishes of India" highlights the importance and significance of the prioritized fish species distributed across various states within the country. This review article systematically documents the current status of state fishes from various perspectives, identifies the prevailing knowledge gaps, and also highlights the issues and strategic plans essential for the conservation and sustainable utilization of these valuable genetic resources. A total of 8357 publications were checked for the consolidated information on state fish species, and the appropriate items were selected under eight categories: biology, physiology and nutrition, aquaculture, habitat and environmental parameters, genetics and biotechnology, harvest and postharvest, fish health management, and others. The synthesized information was used to present the current status of research and development on state fish species. The knowledge gaps that are to be addressed are also depicted under the perspectives of fisheries management and conservation aquaculture. Based on the findings, strategic plans for the targeted conservation programmes are proposed and discussed under various in situ and ex situ conservation measures. Further, the departmental processes involved in the declaration, the importance of stakeholder involvement, namely, local communities and policymakers, in fostering effective conservation measures, and planning for utilization of these valuable fish genetic resources are also indicated.

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.

A chromosome-level genome assembly of Hong Kong catfish (Clarias fuscus) uncovers a sex-determining region

BACKGROUND

Hong Kong catfish (Clarias fuscus) is an ecologically and economically important species that is widely distributed in freshwater regions of southern China. Hong Kong catfish has significant sexual growth dimorphism. The genome assembly of the Hong Kong catfish would facilitate study of the sex determination and evolution mechanism of the species.

RESULTS

The first high-quality chromosome-level genome of the Hong Kong catfish was constructed. The total genome was 933.4 Mb, with 416 contigs and a contig N50 length of 8.52 Mb. Using high-throughput chromosome conformation capture (Hi-C) data, the genome assembly was divided into 28 chromosomes with a scaffold N50 length of 36.68 Mb. A total of 23,345 protein-coding genes were predicted in the genome, and 94.28% of the genes were functionally annotated in public databases. Phylogenetic analysis indicated that C. fuscus and Clarias magur diverged approximately 63.7 million years ago. The comparative genome results showed that a total of 60 unique, 353 expanded and 851 contracted gene families were identified in Hong Kong catfish. A sex-linked quantitative trait locus identified in a previous study was located in a sex-determining region of 30.26 Mb (0.02 to 30.28 Mb) on chromosome 13 (Chr13), the predicted Y chromosome. This QTL region contained 785 genes, of which 18 were identified as sex-related genes.

CONCLUSIONS

This study is the first to report the chromosome-level genome assembly of Hong Kong catfish. The study provides an excellent genetic resource that will facilitate future studies of sex determination mechanisms and evolution in fish.

Genome size estimation and its associations with body length, chromosome number and evolution in teleost fishes

Precise estimation of genome size (GS) is vital for various genomic studies, such as deciding genome sequencing depth, genome assembly, biodiversity documentation, evolution, genetic disorders studies, duplication events etc. Animal Genome Size Database provides GS of over 2050 fish species, which ranges from 0.35 pg in pufferfish (Tetraodon nigroviridis) to 132.83 pg in marbled lungfish (Protopterus aethiopicus). The GS of majority of the fishes inhabiting waters of Indian subcontinent are still missing. In present study, we estimated GS of 51 freshwater teleost (31 commercially important, 7 vulnerable and 13 ornamental species) that ranged from 0.58 pg in banded gourami (Trichogaster fasciata) to 1.92 pg in scribbled goby (Awaous grammepomus). Substantial variation in GS was observed within the same fish orders (0.64-1.45 pg in cypriniformes, 0.70-1.41 pg in siluriformes and 0.58-1.92 pg in perciformes). We examined the relationship between the GS, chromosome number and body length across all the fishes. Body length was found to be associated with GS, whereas no relationship was noticed between the GS and the chromosome number. The analysis using ancestral information revealed haploid chromosome number 25, 27 and 24 for the most recent common ancestor of cypriniformes, siluriformes and perciformes, respectively. The study led to generation of new records on GS of 43 fish species and revalidated records for 8 species. The finding is valuable resource for further research in the areas of fish genomics, molecular ecology and evolutionary conservation genetics.

BAC-FISH Based Physical Map of Endangered Catfish Clarias magur for Chromosome Cataloguing and Gene Isolation through Positional Cloning

Construction of a physical chromosome map of a species is important for positional cloning, targeted marker development, fine mapping of genes, selection of candidate genes for molecular breeding, as well as understanding the genome organization. The genomic libraries in the form of bacterial artificial chromosome (BAC) clones are also a very useful resource for physical mapping and identification and isolation of full-length genes and the related cis acting elements. Some BAC-FISH based studies reported in the past were gene based physical chromosome maps of Clarias magur (magur) to understand the genome organization of the species and to establish the relationships with other species in respect to genes' organization and evolution in the past. In the present study, we generated end sequences of the BAC clones and analyzed those end sequences within the scaffolds of the draft genome of magur to identify and map the genes bioinformatically for each clone. A total of 36 clones mostly possessing genes were identified and used in probe construction and their subsequent hybridization on the metaphase chromosomes of magur. This study successfully mapped all 36 specific clones on 16 chromosome pairs, out of 25 pairs of magur chromosomes. These clones are now recognized as chromosome-specific makers, which are an aid in individual chromosome identification and fine assembly of the genome sequence, and will ultimately help in developing anchored genes' map on the chromosomes of C. magur for understanding their organization, inheritance of important fishery traits and evolution of magur with respect to channel catfish, zebrafish and other species.

Variation in immuno-reproductive milieu of testis in Clarias magur from pre-spawning to spawning phase: An indication towards non-canonical role of immune elements in testes

Immune mechanisms are major players in ensuring the normal functioning of testicular functions. However, apart from their role in active defence against pathogens, prior studies have also suggested a possibility for reproduction-related (non-immune) functions of certain immune elements. This study employs a comparative transcriptomics approach followed by network analysis for tracking the variation in the immuno-reproductive milieu of Clarias magur testis in spawning versus pre-spawning phase. The results show a significant modulation of both reproduction and immune-relevant genes in spawning versus pre-spawning phase. The functional enrichment of the upregulated reproduction-relevant gene network also shows immune-related biological processes which indicates a probability of involvement of these candidates in spermatogenesis-related events for switching from pre-spawning to spawning phase. The upregulated immune network is highly dense with 40 hubs, 10 cluster sub-networks and 142 functionally enriched pathways in comparison to its downregulated counterpart with only 5 hubs, 1 cluster and 1 enriched pathway. These findings indicate that the synchronisation in modulation of both reproductive and immune-related factors is critical for progression of testicular events guiding the switch from pre-spawning to spawning phase. The reproductive phase-dependent variation in plasma sex steroid levels and the selected genes for quantitative PCR also corroborated this hypothesis. The study also serves as a preliminary screening step for probable immune candidates that may be involved in reproductive functions of testis in addition to defence.

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.

Aquaculture omics: An update on the current status of research and data analysis

Aquaculture is the fast-growing agricultural sector and has the ability to meet the growing demand for protein nutritional security for future population. In future aquaculture is going to be the major source of fish proteins as capture fisheries reached at its maximum. However, several challenges need to overcome such as lack of genetically improved strains/varieties, lack of species-specific feed/functional feed, round the year availability of quality fish seed, pollution of ecosystems and increased frequencies of disease occurrence etc. In recent years, the continuous development of high throughput sequencing technology has revolutionized the biological sciences and provided necessary tools. Application of 'omics' in aquaculture research have been successfully used to resolve several productive and reproductive issues and thus ensure its sustainability and profitability. To date, high quality draft genomes of over fifty fish species have been generated and successfully used to develop large number of single nucleotide polymorphism markers (SNPs), marker panels and other genomic resources etc in several aquaculture species. Similarly, transcriptome profiling and miRNAs analysis have been used in aquaculture research to identify key transcripts and expression analysis of candidate genes/miRNAs involved in reproduction, immunity, growth, development, stress toxicology and disease. Metagenome analysis emerged as a promising scientific tool to analyze the complex genomes contained within microbial communities. Metagenomics has been successfully used in the aquaculture sector to identify novel and potential pathogens, antibiotic resistance genes, microbial roles in microcosms, microbial communities forming biofloc, probiotics etc. In the current review, we discussed application of high-throughput technologies (NGS) in the aquaculture sector.

Thermal stress induces expression of Nuclear protein and Parkin genes in endangered catfish, Clarias magur

Fluctuation in water temperature can create thermal stress, which may impact many aspects of fish life, such as survival, growth, reproduction, disease occurrence etc. The endangered catfish, Clarias magur, has been reported to survive at higher thermal stress, even though the exact mechanism is unknown. The genes coding for Nuclear protein 1 (Nupr1) and Parkin E3 ubiquitin protein ligase (Park2) have been reported to protect cells from stress-induced damage and death. In this study, we characterized both the genes and assessed their quantitative expression in C. magur. Structural features of both the genes were found similar to a related catfish, Ictalurus punctatus, and model fish zebrafish. The genes were fairly conserved in fishes as observed through phylogenetic analysis. The real time expression of the two stress-associated genes were also assessed in brain, kidney, liver and muscle tissues of C. magur exposed to warm (34 °C) and cold (15 °C) water. RT-PCR analysis revealed up-regulation in the relative expression levels of Nupr1 and Park2 genes at both temperatures with maximum positive fold change during stress to cold water, even though the posteriori Dunnett's test after ANOVA revealed that there were significant differences between the control and challenged groups. The study indicated that Nupr1 gene plays role in muscle tissue at both high and low thermal stress, but at high thermal stress in liver, while Park2 plays role in muscle, brain and kidney at low temperature and in liver at high temperature stress in C. magur. The study has generated first-hand information under warm- and cold water, which pave the way to understand the expression response of these genes to thermal vacillations and to establish evolutionary significance in catfishes and other species.

Molecular characterization and structural dynamics of Aquaporin1 from walking catfish in lipid bilayers

Aquaporin's (AQPs) are the major superfamily of small integral membrane proteins that facilitates transportation of water, urea, ammonia, glycerol and ions across biological cell membranes. Despite of recent advancements made in understanding the biology of Aquaporin's, only few isoforms of aquaporin 1 (AQP1) some of the teleost fish species have been characterized at molecular scale. In this study, we made an attempt to elucidate the molecular mechanism of water transportation in AQP1 from walking catfish (Clarias batrachus), a model species capable of breathing in air and inhabits in challenging environments. Using state-of-the-art computational modelling and all-atoms molecular dynamics simulation, we explored the structural dynamics of full-length aquaporin 1 from walking catfish (CbAQP1) in lipid mimetic bilayers. Unlike AQP1 of human and bovine, structural ensembles of CbAQP1 from MD revealed discrete positioning of pore lining residues at the intracellular end. Snapshots from MD simulation displayed differential dynamics of aromatic/arginine (ar/R) filter and extracellular loop C bridging transmembrane (TM) helix H3 and H4. Distinct conformation of large extracellular loops, loop bridging TM2 domain and HB helix along with positioning of selectivity filter lining residues controls the permeability of water across the bilayer. Moreover, the identified unique and conserved lipid binding sites with 100% lipid occupancy signifies lipid mediated structural dynamics of CbAQP1. All-together, this is the first ever report on structural-dynamics of aquaporin 1 in walking catfish which will be useful to understand the molecular basis of transportation of water and other small molecules under varying degree of hyperosmotic environment.