Abstract
Circular RNAs (circRNAs) play an important role in regulating the mammalian reproductive system, especially testicular development and spermatogenesis. However, their functions in testicular development and spermatogenesis in the Qianbei Ma goat, the Guizhou endemic breed are still unclear. In this study, tissue sectioning and circRNAs transcriptome analysis were conducted to compare the changes of morphology and circular RNAs gene expression profile at four different developmental stages (0Y, 0-month-old; 6Y, 6-month-old; 12Y, 12-month-old; 18Y, 18-month-old). The results showed that the circumferences and area of the seminiferous tubule gradually increased with age, and the lumen of the seminiferous tubule in the testis differentiated significantly. 12,784 circRNAs were detected from testicular tissues at four different developmental stages by RNA sequencing, and 8,140 DEcircRNAs (differentially expressed circRNAs) were found in 0Y vs. 6Y, 6Y vs. 12Y, 12Y vs. 18Y and 0Y vs. 18Y, 0Y vs. 12Y, 6Y vs. 18Y Functional enrichment analysis of the source genes showed that they were mainly enriched in testicular development and spermatogenesis. In addition, the miRNAs and mRNAs associated with DECircRNAs in 6 control groups were predicted by bioinformatics, and 81 highly expressed DECircRNAs and their associated miRNAs and mRNAs were selected to construct the ceRNA network. Through functional enrichment analysis of the target genes of circRNAs in the network, some candidate circRNAs related to testicular development and spermatogenesis were obtained. Such as circRNA_07172, circRNA_04859, circRNA_07832, circRNA_00032 and circRNA_07510. These results will help to reveal the mechanism of circRNAs in testicular development and spermatogenesis, and also provide some guidance for goat reproduction.