Abstract
Normal testicular development plays a crucial role in male reproduction and is the precondition for spermatogenesis. PIWI-interacting RNAs (piRNAs) are novel noncoding RNAs expressed in animal germ cells that form complexes with PIWI family proteins and are involved in germ cell development, differentiation, and spermatogenesis. However, changes in piRNA expression profiles during early testicular development in yak have not been investigated. In this study, we used small RNA sequencing to evaluate the differences and potential functions of piRNA expression profiles in 6-, 18-, and 30-month-old yak testis tissues. Differential expression analysis found 109, 293, and 336 differentially expressed piRNAs in M30 vs. M18, M18 vs. M6, and M30 vs. M6, respectively, and found 30 common differentially expressed piRNAs in the three groups of M6, M18, and M30. In addition, the functional enrichment analysis of differentially expressed piRNAs target genes indicated that they were related to testicular development and spermatogenesis. Finally, we detected the expression of the PIWI protein family in the yak testis at different developmental stages and found that PIWIL1, PIWIL2, PIWIL3, and PIWIL4 were highly expressed in 18- and 30-month-old yak testis and almost not expressed in 6-month-old yak testis. In conclusion, this study summarizes the changes of piRNA expression patterns during the early development of yak testis and provides new clues for the regulatory role of piRNA in yak testis.