Abstract
BACKGROUND: The Testis is an important reproductive organ in male mammals and the site for spermatogenesis, androgen synthesis, and secretion. Non-coding RNAs (ncRNAs) play an important regulatory role in various biological processes. However, the regulatory role of ncRNAs in the development of yak testes and spermatogenesis remains largely unclear. RESULT: In this study, we compared the expression profiles of circular RNAs (circRNAs), microRNAs (miRNAs), and messenger RNAs (mRNAs) in yak testicular tissue samples collected at 6 months (Y6M), 18 months (Y18M), and 4 years (Y4Y). Using RNA sequencing (RNA-Seq), we observed a significant difference in the expression patterns of ncRNAs in the samples collected at different testicular development stages. Twenty-two differentially expressed (DE) circRNAs, 69 DE miRNAs, and 64 DE mRNAs were detected in Y6M, Y18M, and Y4Y testicular samples, respectively. The results of gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that the source genes of DE circRNAs, predicted target genes of DE miRNAs, and DE mRNAs were specifically associated with signaling pathways and GO terms that were related to sperm synthesis, sperm vitality, and testicular development, such as cell cycle, Wnt signaling pathway, MAPK signaling pathway, GnRH signaling pathway, and spermatogenesis. The analysis of the circRNA-miRNA-mRNA network revealed that some DE ncRNAs, including miR-574, miR-449a, CDC42, and CYP11A1, among others, may be involved in testicular spermatogenesis. Concurrently, various circRNA-miRNA interaction pairs were observed. CONCLUSION: Our findings provide a database of circRNAs, miRNAs, and mRNAs expression profiles in testicular tissue of yaks at different developmental stages and a detailed understanding of the regulatory network of ncRNAs in yak testicular development and provide data that can help elucidate the molecular mechanisms underlying yak testicular development.