Objective Sturgeons are considered living fossils, and have a very high conservation and economic value. Studies on the molecular mechanism of sturgeon gonadal development and sex differentiation would not only aid in understanding vertebrate sex determination but also benefit sturgeon aquaculture. Piwi-interacting RNAs (piRNAs) have been shown to function in germline or gonadal development. In this study, we performed small RNA deep sequencing and microarray hybridization to identify potential sturgeon piRNAs. Methods Male and female sturgeon gonads were collected and used for small RNA sequencing on an Illumina HiSeq platform with the validation of piRNA expression by microarray chip. The program Bowtie and k-mer scheme were performed to filter small RNA reads and discover potential sturgeon piRNAs. A known piRNA database, the coding sequence (CDS), 5′ and 3′ untranslated region (UTR) database of the A. Schrenckii transcriptome, Gene Ontology (GO) database and KEGG pathway database were searched subsequently to analyze the potential bio-function of sturgeon piRNAs. Results A total of 875,679 putative sturgeon piRNAs were obtained, including 93 homologous to known piRNAs and hundreds showing sex-specific and sex-biased expression. Further analysis showed that they are predominant in both the ovaries and testes and those with a sex-specific expression pattern are nearly equally distribution between sexes. This may imply a relevant role in sturgeon gonadal development. KEGG pathway and GO annotation analyses indicated that they may be related to sturgeon reproductive processes. Conclusion Our study provides the first insights into the gonadal piRNAs in a sturgeon species and should serve as a useful resource for further elucidation of the gene regulation involved in the sex differentiation of vertebrates. These results should also facilitate the technological development of early sex identification in sturgeon aquaculture
