Expression profiles of oviductal mRNAs and lncRNAs in the follicular phase and luteal phase of sheep (Ovis aries) with 2 fecundity gene (FecB) genotypes

FecB (also known as BMPR1B) is a crucial gene in sheep reproduction, which has a mutation (A746G) that was found to increase the ovulation rate and litter size. The FecB mutation is associated with reproductive endocrinology, such mutation can control external estrous characteristics and affect follicle-stimulating hormone during the estrous cycle. Previous researches showed that the FecB mutation can regulate the transcriptomic profiles in the reproductive-related tissues including hypothalamus, pituitary, and ovary during the estrous cycle of small-tailed Han (STH) sheep. However, little research has been reported on the correlation between FecB mutation and the estrous cycle in STH sheep oviduct. To investigate the coding and noncoding transcriptomic profiles involved in the estrous cycle and FecB in the sheep oviduct, RNA sequencing was performed to analyze the transcriptomic profiles of mRNAs and long noncoding RNAs (lncRNAs) in the oviduct during the estrous cycle of STH sheep with mutant (FecBBB) and wild-type (FecB++) genotypes. In total, 21,863 lncRNAs and 43,674 mRNAs were screened, the results showed that mRNAs had significantly higher expression levels than the lncRNAs, and the expression levels of these screened transcripts were lower in the follicular phase than they were in the luteal phase. Among them, the oviductal glycoprotein gene (OVGP1) had the highest expression level. In the comparison between the follicular and luteal phases, 57 differentially expressed (DE) lncRNAs and 637 DE mRNAs were detected, including FSTL5 mRNA and LNC_016628 lncRNA. In the comparison between the FecBBB and FecB++ genotypes, 26 DE lncRNAs and 421 DE mRNAs were detected, including EEF1D mRNA and LNC_006270 lncRNA. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes functional enrichment analysis indicated that the DE mRNAs were enriched mainly in terms related to reproduction such as the tight junction, SAGA complex, ATP-binding cassette, nestin, and Hippo signaling pathway. The interaction network between DE lncRNAs and DE mRNAs indicated that LNC_018420 may be the key regulator in sheep oviduct. Together, our results can provide novel insights into the oviductal transcriptomic function against a FecB mutation background in sheep reproduction.