Gonadal transcriptomic analysis of yellow catfish (Pelteobagrus fulvidraco): identification of sex-related genes and genetic markers

Yellow catfish (Pelteobagrus fulvidraco) has been recognized as a vital freshwater aquaculture species in East and Southeast Asia. In addition to its commercial interest, it is also attracted much attention because of its value in studying sex-determination mechanisms. A comprehensive gonadal transcriptome analysis is believed to provide a resource for genome annotation, candidate gene identification, and molecular marker development. Herein, we performed a de novo assembly of yellow catfish gonad transcriptome by high-throughput Illumina sequencing. A total of 82,123 contigs were obtained, ranging from 351 to 21,268 bp, and N50 of 2,329 bp. Unigenes of 21,869 in total were identified. Of these, 229 and 1,188 genes were found to be specifically expressed in XY gonad tissue for 1 yr and 2 yr old yellow catfish, respectively; correspondingly, 51 and 40 genes were identified in XX gonad tissue at those two stages. Gene ontology and KEGG analysis were conducted and classified all contigs into different categories. A large number of unigenes involved in sex determination were identified, as well as microsatellites and SNP variants. The expression patterns of sex-related genes were then validated by quantitative real-time PCR (qRT-PCR) suggesting the high reliability of RNA-Seq results. In this study, the transcriptome of yellow catfish gonad was first sequenced, assembled, and characterized; it provides a valuable genomic resource for better understanding of yellow catfish sex determination as well as development of molecular markers, thereby assisting in the production of monosex yellow catfish for aquaculture.