Identification of Circular RNAs of Testis and Caput Epididymis and Prediction of Their Potential Functional Roles in Donkeys

Circular RNAs (circRNAs) are a class of noncoding RNAs with a covalently closed loop. Studies have demonstrated that circRNA can function as microRNA (miRNA) sponges or competing endogenous RNAs. Although circRNA has been explored in some species and tissues, the genetic basis of testis development and spermatogenesis in donkeys remain unknown. We performed RNA-seq to detect circRNA expression profiles of adult donkey testes. Length distribution and other characteristics were shown a total of 1971 circRNAs were differentially expressed and 12,648 and 6261 circRNAs were detected from the testis and caput epididymis, respectively. Among these circRNAs, 1472 circRNAs were downregulated and 499 circRNAs were upregulated in the testis. Moreover, KEGG pathway analyses and Gene Ontology were performed for host genes of circRNAs. A total of 39 upregulated circRNA host genes were annotated in spermatogenesis terms, including PIWIL2, CATSPERD, CATSPERB, SPATA6, and SYCP1. Other host genes were annotated in the focal adhesion, Rap1 signaling pathway. Downregulated expressed circRNA host genes participated in the TGF-β signaling pathway, GnRH signaling pathway, estrogen signaling pathway, and calcium signaling pathway. Our discoveries provide a solid foundation for identifying and characterizing critical circRNAs involved in testis development or spermatogenesis.