Abstract
Long non-coding RNA (LncRNA) have been identified as important regulators in the hypothalamic-pituitary-ovarian axis associated with sheep prolificacy. However, their expression pattern and potential roles in the pituitary are yet unclear. To explore the potential mRNAs and lncRNAs that regulate the expression of the genes involved in sheep prolificacy, we used stranded specific RNA-seq to profile the pituitary transcriptome (lncRNA and mRNA) in high prolificacy (genotype FecB BB, litter size = 3; H) and low prolificacy sheep (genotype FecB B+; litter size = 1; L). Our results showed that 57 differentially expressed (DE) lncRNAs and 298 DE mRNAs were found in the pituitary between the two groups. The qRT-PCR results correlated well with the RNA-seq results. Moreover, functional annotation analysis showed that the target genes of the DE lncRNAs were significantly enriched in pituitary function, hormone-related pathways as well as response to stimulus and some other terms related to reproduction. Furthermore, a co-expression network of lncRNAs and target genes was constructed and reproduction related genes such as SMAD2, NMB and EFNB3 were included. Lastly, the interaction of candidate lncRNA MSTRG.259847.2 and its target gene SMAD2 were validated in vitro of sheep pituitary cells. These differential mRNA and lncRNA expression profiles provide a valuable resource for understanding the molecular mechanisms underlying Hu sheep prolificacy.