Expression dynamics, relationships, and transcriptional regulations of diverse transcripts in mouse spermatogenic cells

Among all tissues of the metazoa, the transcritpome of testis displays the highest diversity and specificity. However, its composition and dynamics during spermatogenesis have not been fully understood. Here, we have identified 20,639 message RNAs (mRNAs), 7,168 long non-coding RNAs (lncRNAs) and 15,101 circular RNAs (circRNAs) in mouse spermatogenic cells, and found many of them were specifically expressed in testes. lncRNAs are significantly more testis-specific than mRNAs. At all stages, mRNAs are generally more abundant than lncRNAs, and linear transcripts are more abundant than circRNAs. We showed that the productions of circRNAs and piRNAs were highly regulated instead of random processes. Based on the results of a small-scale functional screening experiment using cultured mouse spermatogonial stem cells, many evolutionarily conserved lncRNAs are likely to play roles in spermatogenesis. Typical classes of transcription factor binding sites are enriched in the promoters of testis-specific m/lncRNA genes. Target genes of CREM and RFX2, 2 key TFs for spermatogenesis, were further validated by using ChIP-chip assays and RNA-seq on RFX2-knockout spermatogenic cells. Our results contribute to the current understanding of the transcriptomic complexity of spermatogenic cells and provide a valuable resource from which many candidate genes may be selected for further functional studies.