Genome-wide Association Study Identifies Genomic Loci of Sex Determination, Gonadal Weight and Gonadosomatic Index Traits in Takifugu Bimaculatus
Mei Wang,Leibin Li,Huanlin Lin,Zhixiong Zhou,Bo Liu,Jianxing Zhong,Fei Pu,Yue Shi,Tao Zhou,Peng Xu
DOI: https://doi.org/10.1016/j.aquaculture.2021.737389
IF: 5.135
2021-01-01
Aquaculture
Abstract:Takifugu bimaculatus is a cultured marine teleost fish, distributed ranging from the southern Yellow Sea to the South China Sea. Understanding the mechanism of sex determination is essential for sex control breeding biotechnologies of T. bimaculatus. However, there lacks a rapid and accurate sex specific marker for the T. bimaculatus. In this study, whole genome resequencing was conducted on 71 individuals (including 49 females and 22 males), 331,190 high-quality polymorphic loci (34,435 INDEL, 296,755 single-nucleotide polymorphisms (SNPs)) were obtained through genome-wide association analysis (GWAS) by GEMMA (Genome-wide Efficient Mixed Model Association algorithm). 1403 loci on chromosome 7 were identified as sex significantly associated loci (P < 1.51 x 10(-7)), and a 12 Mb sex-associated region on Chr7 was narrowed down by F-st, scanning. It was worth noting that the SNP, SNP1100879 on Chr7 at position 8,764,381 was identified to be sex determination in T. bimaculatus. 13 pairs of primers for sex identification Marker-F/R were designed and it was validated by PCR. However, only one pair of primers (SNP-7F/R) was sex-linked, in which 95% of female genotypes are homozygous CC and 88% of male genotypes are heterozygous CT, suggesting a XX/XY sex determination system in T. bimaculatus. In addition, 12 sex-related candidate genes were identified around 306 SNPs, including amhr2, dzip1, patz1, brdt, hmga2, gli1, inhbb, mmp19, wnt2b, taf4 and so on. Because of the significant difference of gonadal weight (GW) and gonadosomatic index (GSI) between females and males of T. bimaculatus, GWAS on GW and GSI of different genders were also conducted independently. Besides, we identified 11 SNPs, located on Chr7, Chr18, Chr21, and Chr2, showing genome-wide significant association with GW and GSI and five genes including cyp27b1, hexb, inhba, cgba and shbg, involved in oogenesis, ovarian follicle development, and regulation of estradiol secretion. Then, we performed the gene ontology (GO) and KEGG signaling pathway (KEGG) and Protein-protein interaction (PPI) network analyses of the candidate genes. amhr2, gli1, ihbb and dzip1 genes were identified as hub genes. In the present study, a sex-linked marker of T. bimaculatus was identified for the first time using whole-genome resequencing technology, and its limited versatility provided important clues for deciphering the sex determination mechanism in T. bimaculatus.