Genome-wide association study reveals the genetic basis of growth trait in yellow catfish with sexual size dimorphism

Chr23-miR-200s and Dmrt1 Control Sexually Dimorphic Trade-Off Between Reproduction and Growth in Zebrafish

In animals, a trade-off exists between reproduction and growth, which are the most fundamental traits. Males and females exhibit profound differences in reproduction and growth in fish species. However, the precise molecular mechanism governing this phenomenon is still not clear. Here, we uncovered that chr23-miR-200s and dmrt1 knockout specifically caused an impairment in reproduction and an increase in body growth in female and male zebrafish, respectively. Chr23-miR-200s and Dmrt1 directly regulate the stat5b gene by targeting its 3'UTR and promoter. The loss of stat5b completely abolished the elevated growth performance in chr23-miR-200s-KO or dmrt1-/- zebrafish. Moreover, the dmrt1 transgenic zebrafish had significantly lower body length and body weight than the control males, accompanied by a significant reduction in stat5b expression in the liver of transgenic fish. In summary, our study proposes a regulatory model elucidating the roles of chr23-miR-200s and Dmrt1 in controlling the sexually dimorphic trade-off between reproduction and growth.

Genome-wide association study for growth traits with 1066 individuals in largemouth bass (Micropterus salmoides)

The largemouth bass is a native species of North America that was first introduced to mainland China in the 1980s. In recent years, it has been extensively farmed in China due to its high meat quality and broad adaptability. In this study, we collected growth trait data from 1,066 largemouth bass individuals across two populations. We generated an average of approximately 7× sequencing coverage for these fish using Illumina sequencers. From the samples, we identified 2,695,687 SNPs and retained 1,809,116 SNPs for further analysis after filtering. To estimate the number of genome-wide effective SNPs, we performed LD pruning with PLINK software and identified 77,935 SNPs. Our GWAS revealed 15 SNPs associated with six growth traits. We identified a total of 24 genes related to growth, with three genes-igf1, myf5, and myf6-directly associated with skeletal muscle development and growth, located near the leading SNP on chromosome 23. Other candidate genes are involved in the development of tissues and organs or other physiological processes. These findings provide a valuable set of SNPs and genes that could be useful for genetic breeding programs aimed at enhancing growth in largemouth bass.

Multi-omics analysis identifies sex-specific hepatic protein-metabolite networks in yellow catfish (Pelteobagrus fulvidraco) exposed to chronic hypoxia

Hypoxia disrupts the endocrine system of teleosts. The liver plays important roles in the endocrine system, energy storage, and metabolic processes. The aim of this study was to investigate the sex-specific hepatic response of yellow catfish under chronic hypoxia at the multi-omics level. Common hepatic responses in both sexes included the HIF-1 signaling pathway, glycolysis/gluconeogenesis, and steroid biosynthesis. Hypoxia dysregulated primary bile acid biosynthesis, lipid metabolism, and vitellogenin levels in female fish. Endoplasmic reticulum function in females also tended to be disrupted by hypoxia, as evidenced by significantly enriched pathways, including ribosome, protein processing in the endoplasmic reticulum, and RNA degradation. Other pathways, including the TCA cycle, oxidative phosphorylation, and Parkinson's and Huntington's disease, were highly enriched by hypoxia in male fish, suggesting that mitochondrial function was dysregulated. In both sexes of yellow catfish, the cell cycle was arrested and apoptosis was inhibited under chronic hypoxia. Multi-omics suggested that SLC2A5, CD209, LGMN, and NEDD8 served as sex-specific markers in these fish under chronic hypoxia. Our results provide insights into hepatic adaptation to chronic hypoxia and facilitate our understanding of sex-specific responses in fish.

Identification and Functional Analysis of Cynoglossus semilaevis Z-Linked E3 Ubiquitin Ligase rnf34

The high proportion of males in C. semilaevis hinders their industrial development. The genetic ZW individual can become a pseudomale by sex reversal. And the pseudomale can produce Z-sperm (with epigenetic information to cause sex reversal) while W-sperm is absent, which leads to an even higher male proportion in offspring. Recently, with the development of transcriptomic technologies, research on spermatogenesis in C. semilaevis has been focused on the ubiquitination pathway. In this study, we analyzed the function of the ubiquitin ligase rnf34 gene on the Z chromosome. A qPCR experiment showed that its expression level in the gonad was the highest among different tissues. In the ovary, the expression gradually increased with development from 40 days post-hatching (dph) to 1.5 years post-hatching (yph). In the testis, rnf34 showed increased expression from 40 dph to 6 months post-hatching (mpf) and stabilized up until 1.5 ypf. In situ hybridization showed that the mRNA of rnf34 was mainly distributed in the germ cells of the testis and the ovary. In vivo siRNA-mediated knockdown of the rnf34 gene in male fish affected the expression of a series of genes related to sex differentiation and spermatogenesis. These results provide genetic data on the molecular mechanisms of gonadal development and spermatogenesis in C. semilaevis.

Genome-Wide Association Analysis of Heat Tolerance in F2 Progeny from the Hybridization between Two Congeneric Oyster Species

As the world's largest farmed marine animal, oysters have enormous economic and ecological value. However, mass summer mortality caused by high temperature poses a significant threat to the oyster industry. To investigate the molecular mechanisms underlying heat adaptation and improve the heat tolerance ability in the oyster, we conducted genome-wide association analysis (GWAS) analysis on the F2 generation derived from the hybridization of relatively heat-tolerant Crassostrea angulata ♀ and heat-sensitive Crassostrea gigas ♂, which are the dominant cultured species in southern and northern China, respectively. Acute heat stress experiment (semi-lethal temperature 42 °C) demonstrated that the F2 population showed differentiation in heat tolerance, leading to extremely differentiated individuals (approximately 20% of individuals die within the first four days with 10% survival after 14 days). Genome resequencing and GWAS of the two divergent groups had identified 18 significant SNPs associated with heat tolerance, with 26 candidate genes located near these SNPs. Eleven candidate genes that may associate with the thermal resistance were identified, which were classified into five categories: temperature sensor (Trpm2), transcriptional factor (Gata3), protein ubiquitination (Ube2h, Usp50, Uchl3), heat shock subfamily (Dnajc17, Dnaja1), and transporters (Slc16a9, Slc16a14, Slc16a9, Slc16a2). The expressional differentiation of the above genes between C. gigas and C. angulata under sublethal temperature (37 °C) further supports their crucial role in coping with high temperature. Our results will contribute to understanding the molecular mechanisms underlying heat tolerance, and provide genetic markers for heat-resistance breeding in the oyster industry.

A chromosome-level genome assembly of the darkbarbel catfish Pelteobagrus vachelli

The darkbarbel catfish (Pelteobagrus vachelli), an economically important aquaculture species in China, is extensively employed in hybrid yellow catfish production due to its superior growth rate. However, information on its genome has been limited, constraining further genetic studies and breeding programs. Leveraging the power of PacBio long-read sequencing and Hi-C technologies, we present a high-quality, chromosome-level genome assembly for the darkbarbel catfish. The resulting assembly spans 692.10 Mb, with an impressive 99.9% distribution over 26 chromosomes. The contig N50 and scaffold N50 are 13.30 Mb and 27.55 Mb, respectively. The genome is predicted to contain 22,109 protein-coding genes, with 96.1% having functional annotations. Repeat elements account for approximately 35.79% of the genomic landscape. The completeness of darkbarbel catfish genome assembly is highlighted by a BUSCO score of 99.07%. This high-quality genome assembly provides a critical resource for future hybrid catfish breeding, comparative genomics, and evolutionary studies in catfish and other related species.

Comparative transcriptome profiles of four sexually size dimorphic fish

Sexual size dimorphism is widespread in fish species. Although sex growth differences in multiple species have been studied successively, the commonalities of regulatory mechanisms across sexually dimorphic species are unknown. In this study, we performed RNA-seq analysis of four representative fish (loach, half-smooth tongue sole, yellow catfish, and Nile tilapia) with significant growth differences between females and males. Clean reads were identified from four fish species, ranging from 45,718,052 to 57,733,120. Following comparison transcriptome analysis, there were 1,132 and 1,108, 1,290 and 1,102, 4,732 and 4,266, 748 and 192 differentially expressed genes (DEGs) in the brain and muscle of loach, half-smooth tongue sole, yellow catfish, and Nile tilapia, respectively. Furthermore, the expression levels were validated by quantitative real-time PCR (qRT-PCR). Comparative transcriptome profiles of four fish described here will provide fundamental information for further studies on the commonalities of sexually size dimorphic fish in regulating growth differences between females and males.