Sex determination mechanisms in teleost fish broadly differ from mammals and birds, with sex chromosomes that are far less differentiated and recombination often occurring along the length of the X and Y chromosomes, posing major challenges for the identification of specific sex determination genes. Here, we take an innovative approach of comparative genome analysis of the genomic sequences of the X chromosome and newly sequenced Y chromosome in the channel catfish. Using a YY channel catfish as the sequencing template, we generated, assembled, and annotated the Y genome sequence of channel catfish. The genome sequence assembly had a contig N50 size of 2.7 Mb and a scaffold N50 size of 26.7 Mb. Genetic linkage and GWAS analyses placed the sex determination locus within a genetic distance less than 0.5 cM and physical distance of 8.9 Mb. However, comparison of the channel catfish X and Y chromosome sequences showed no sex-specific genes. Instead, comparative RNA-Seq analysis between females and males revealed exclusive sex-specific expression of an isoform of the breast cancer anti-resistance 1 (BCAR1) gene in the male during early sex differentiation. Experimental knockout of BCAR1 gene converted genetic males (XY) to phenotypic females, suggesting BCAR1 as a putative sex determination gene. We present the first Y chromosome sequence among teleost fish, and one of the few whole Y chromosome sequences among vertebrate species. Comparative analyses suggest that sex-specific isoform expression through alternative splicing may underlie sex determination processes in the channel catfish, and we identify BCAR1 as a potential sex determination gene.

Spotted scat (Scatophagus argus) can tolerate a wide range of salinity fluctuations. It is a good model for studying environmental salinity adaptation. Lipid metabolism plays an important role in salinity adaptation in fish. To elucidate the mechanism of lipid metabolism in the osmoregulation, the liver transcriptome was analyzed after 22 d culture with a salinity of 5 ppt (Low-salinity group: LS), 25 ppt (Control group: Ctrl), and 35 ppt (High-salinity group: HS) water by using RNA sequencing (RNA-seq) in spotted scat. RNA-seq analysis showed that 1276 and 2768 differentially expressed genes (DEGs) were identified in the LS vs. Ctrl and HS vs. Ctrl, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the pathways of steroid hormone biosynthesis, steroid biosynthesis, glycerophospholipid metabolism, glycerolipid metabolism, and lipid metabolism were significantly enriched in the LS vs. Ctrl. The genes of steroid biosynthesis (sqle, dhcr7, and cyp51a1), steroid hormone biosynthesis (ugt2a1, ugt2a2, ugt2b20, and ugt2b31), and glycerophospholipid metabolism (cept1, pla2g4a, and ptdss2) were significantly down-regulated in the LS vs. Ctrl. The pathways related to lipid metabolisms, such as fatty acid metabolism, fatty acid biosynthesis, peroxisome proliferator-activated receptor (PPAR) signaling pathway, adipocytokine signaling pathway, fatty acid degradation, and unsaturated fatty acid biosynthesis, were significantly enriched in the HS vs. Ctrl. The genes of unsaturated fatty acid biosynthesis (scd1, hacd3, fads2, pecr, and elovl1) and adipocytokine signaling pathway (g6pc1, socs1, socs3, adipor2, pck1, and pparα) were significantly up-regulated in the HS vs. Ctrl. These results suggest that the difference in liver lipid metabolism is important to adapt to low- and high-salinity stress in spotted scat, which clarifies the molecular regulatory mechanisms of salinity adaptation in euryhaline fish.

Silver sillago (Sillago sihama) is a commercially important marine fish species in East Asia. In this study, we compared the transcriptome response to hypoxia stress in the gill tissue of S. sihama. The fish were divided into four groups, such as 1 h of hypoxia (hypoxia1h, DO = 1.5 ± 0.1 mg/L), 4 h of hypoxia (hypoxia4h, DO = 1.5 ± 0.1 mg/L), 4 h of reoxygen (reoxygen4h, DO = 8.0 ± 0.2 mg/L) after 4 h of hypoxia (DO = 1.5 mg/L), and normoxia or control (DO = 8.0 ± 0.2 mg/L) groups. Compared to the normoxia group, a total of 3550 genes were identified as differentially expressed genes (DEGs) (log₂foldchange > 1 and padj < 0.05), including 1103, 1451 and 996 genes in hypoxia1h, hypoxia4h and reoxygen4h groups, respectively. Only 247 DEGs were differentially co-expressed in all treatment groups. According to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, DEGs were significantly enriched in steroid biosynthesis, biosynthesis of amino acids, glutathione metabolism and metabolism of xenobiotics by cytochrome P450, ferroptosis and drug metabolism—cytochrome P450 pathways. Of these, the cytochrome P450 (CYP) and glutathione S-transferase (GST) gene families were widely expressed. Our study represents the insights into the underlying molecular mechanisms of hypoxia stress.

Growth hormone (GH) has been considered as a candidate gene for growth traits in fish. In this study, polymorphisms of the GH gene were evaluated for associations with growth traits in 282 Siniperca chuatsi individuals. Using directly sequencing, four single nucleotide polymorphisms (SNPs) were identified in GH gene, with two mutations in intron 4 (g.4940A>C, g.4948A>T), one mutation in exon 5 (g.5045T>C) and one in intron 5 (g.5234T>G). Notably, three of them were significantly associated with growth performance, particularly for g.4940A>C which was highly correlated with all the four growth traits. In conclusion, our results demonstrated that these SNPs in GH gene could influence growth performance of S.chuatsi and could be used for marker-assisted selection (MAS) in this species.

Chronic heat stress can have detrimental effects on the survival of fish. This study aimed to investigate the impact of prolonged high temperatures on the growth, antioxidant capacity, apoptosis, and transcriptome analysis of Hong Kong catfish (Clarias fuscus). By analyzing the morphological statistics of C. fuscus subjected to chronic high-temperature stress for 30, 60, and 90 days, it was observed that the growth of C. fuscus was inhibited compared to the control group. The experimental group showed a significant decrease in body weight and body length compared to the control group after 60 and 90 days of high-temperature stress (p < 0.05, p < 0.01). A biochemical analysis revealed significant alterations in the activities of three antioxidant enzymes superoxide dismutase activity (SOD); catalase activity (CAT); glutathione peroxidase activity (GPx), the malondialdehyde content (MDA), and the concentrations of serum alkaline phosphatase (ALP); Aspartate aminotransferase (AST); and alanine transaminase (ALT) in the liver. TUNEL staining indicated stronger apoptotic signals in the high-temperature-stress group compared to the control group, suggesting that chronic high-temperature-induced oxidative stress, leading to liver tissue injury and apoptosis. Transcriptome analysis identified a total of 1330 DEGs, with 835 genes being upregulated and 495 genes being downregulated compared to the control group. These genes may be associated with oxidative stress, apoptosis, and immune response. The findings elucidate the growth changes in C. fuscus under chronic high temperature and provide insights into the underlying response mechanisms to a high-temperature environment.