Genome-Wide Association Studies (GWAS) and Transcriptome Analysis Reveal Male Heterogametic Sex-Determining Regions and Candidate Genes in Northern Snakeheads (Channa argus)

The Northern snakehead (Channa argus) is a significant economic aquaculture species in China. Exhibiting sexual dimorphism in the growth rate between females and males, mono-sex breeding holds substantial value for aquaculture. This study employed GWAS and transcriptome analysis were applied to identify sex determination genomic regions and develop sex-specific markers. A total of 270 single-nucleotide polymorphisms (SNPs) and 31 insertion-deletions (InDels) were identified as being sexually dimorphic through GWAS and fixation index (Fst) scanning. Based on GWAS results, two sex-specific InDel markers were developed, effectively distinguishing genetic sex for XX females, XY males, and YY super-males via (polymerase chain reaction) PCR amplification. A major genomic segment of approximately 115 kb on chromosome 3 (Chr 03) was identified as the sex-determination region. A comparative transcriptome analysis of gonads for three sexes identified 158 overlapping differentially expressed genes (DEGs). Additionally, three sex-related candidate genes were identified near the sex determination region, including id2, sox11, and rnf144a. Further studies are required to elucidate the functions of these genes. Overall, two sex-specific InDel markers support a male heterogametic XX/XY sex-determination system in Northern snakeheads and three candidate genes offer new insights into sex determination and the evolution of sex chromosomes in teleost fish.