Quantifying the expression of mRNAs in the gonads at the critical stage of molecular sex differentiation stage might help to clarify the regulatory network during early sex differentiation and provide new information on the role of sex-related genes in gonadal function. In this study, transcriptomic analysis of sex-related genes expression profiles in fugu gonads at 60 and 90 days after hatching (dah) was conducted firstly, and a total of 112,504,991 clean reads, encompassing 28.35 Gb of sequences were retrieved. Twenty-three thousand eight hundred ten genes were found to be expressed in juvenile fugu gonads, and we mainly focused on the differentially expressed genes that have the potential to be involved in the gonadal sex differentiation. For 60-dah juveniles, we identified 1014 genes that were upregulated in the ovary and 1570 that were upregulated in the testis. For 90-dah juveniles, we identified 1287 genes that were upregulated in the ovary and 1500 that were upregulated in the testis. The dimorphic expression patterns of 15 genes in gonads at 30 and 40 dah were further investigate using qPCR. Cyp11b and star were expressed at higher levels in XY than in XX, while cyp11a1 and cyp19a1a were expressed at higher levels in XX than in XY at 30 dah. At 40 dah, the levels of gsdf, dmrt1, dmrt3, cyp11c1, star, and hsd3b expression were higher in XY, while the levels of foxl2, cyp19a1a, wnt9b, and foxD4 expression were higher in XX. Sox9, cyp11a1, cyp17a1, cyp17a2, and nr5a2 were expressed at similar levels in XX and XY at 40 dah. This is the first report of gonadal transcriptome of fugu at early sex differentiation stage, and our results provide an archive for further study on molecular mechanism underlying sex differentiation in this species.
Keywords: Differentially expressed genes; Early sex differentiation; Gonadal transcriptome analysis; Takifugu rubripes.