Comprehensive analysis of long noncoding RNA and mRNA expression patterns in sheep testicular maturation

Long noncoding RNAs (LncRNAs) have been identified as important regulators of testis development; however, their expression patterns and roles in sheep are not yet clear. Thus, we used stranded specific RNA-seq to profile the testis transcriptome (lncRNAs and mRNAs) in premature and mature sheep. Hormone levels and the testis index were examined, and histological analyses were performed at five stages of testis development, 5-day-old (D5), 3-month-old (3M), 6-month-old (6M), 9-month-old (9M), and 2-year-old (2Y), the results of which indicate a significant difference in hormone levels and testis morphometries between the 3M and 9M stages (P < 0.05). Based on the comparison between 3M and 9M samples, we found 1,118 differentially expressed (DE) lncRNAs and 7,253 DE mRNAs in the testes, and qRT-PCR analysis showed that the results correlated well with the transcriptome data. Furthermore, we constructed lncRNA-protein-coding gene interaction networks. Forty-seven DE lncRNA-targeted genes enriched for male reproduction were obtained by cis- and trans-acting; 51 DE lncRNAs and 45 cis-targets, 2 DE lncRNAs and 2 trans-targets were involved in this network. Of these, 5 lncRNAs and their targets, PRKCD, NANOS3, SERPINA5, and CYP19A1, were enriched for spermatogenesis and male gonad development signaling pathways. We further examined the expression levels of 5 candidate lncRNAs and their target genes during testis development. Lastly, the interaction of lncRNA TCONS__00863147 and its target gene PRKCD was validated in vitro in sheep Leydig cells. This study provides a valuable resource for further study of lncRNA function in sheep testis development and spermatogenesis.