RNA-Seq analysis of ovary and testis reveals potential differentially expressed transcripts associated with gonadal unsynchronization development in Onychostoma macrolepis

The Onychostoma macrolepis (O. macrolepis) is a rare and endangered wild species. Their endangered extinction might be due to their low fertility. To further illustrate the molecular mechanism of gonad development of the male and female O. macrolepis, the present study carried out de novo testicular and ovarian transcriptome sequencing. By comparing ovary and testis, 30,869 differentially expressed unigenes (9870 in female, 20999 in male) were identified. In addition, KEGG and GO analysis suggested that the Hedgehog signaling pathway have important roles in testis maintenance and spermatogenesis, whereas the Hippo signaling pathway and Wnt signaling pathway are likely to participate in ovary maintenance. RT-qPCR analysis results were consistent with transcriptome sequencing that all of gender differentiation-related genes (FOXL2, GDF9, WNT4, CYP19A1, SOX9 and GATA4), temperature-enriched genes (NOVA1, CTGF and NR4A1), clock-related genes (PER2, PER3, CRY1, CRY2, BMAL1 and CIPC) were significantly differential expression in testis compared with ovaries. Taken together, these results revealed a potential molecular mechanism that low fertility of the O. macrolepis might strong correlate with the gonadal dyssynchrony development of the male and female, which might provide theoretical basis and technical support for artificial reproduction and germplasm resource protection of the O. macrolepis.