Impact of FecB Mutation on Ovarian DNA Methylome in Small-Tail Han Sheep

Single-cell RNA-seq reveals the effects of the FecB mutation on the transcriptome profile in ovine cumulus cells

Genetic variations in the ovine ovulation rate, which are associated with the FecB mutation, provide useful models by which to explore the mechanisms regulating the development of mammalian antral follicles. In order to study the effects of the FecB mutation on cumulus cell differentiation, preovulatory follicles were aspirated and cumulus cells were isolated from three FecB genotypes (homozygous, heterozygous and wild type) of Small Tail Han (STH) sheep superstimulated with FSH. Transcriptome information from tens of thousands of cumulus cells was determined with the 10 × Genomics single-cell RNA-seq technology. Under the superovulation treatment, the observed number of preovulatory follicles in the ovaries of FecB carriers was still significantly higher than that in the wild-type (P < 0.05). The expression patterns of cumulus cells differed between FecB carriers and wild-type ewes. The screened cumulus cells could also be further divided into different cell clusters, and the differentiation states and fates of each group of cumulus cells also remained different, which supports the notion that heterogeneity in gene expression is prevalent in single cells. The oxidative phosphorylation pathway was significantly enriched in differentially expressed genes among the cell differentiation branch nodes of cumulus cells and among the differentially expressed genes of cumulus cells from the three genotypes. Combined with the important role of oxidative phosphorylation in the maturation of COCs, we suggest that the oxidative phosphorylation pathway of cumulus cells plays a crucial role in the differentiation process of cumulus cells and the mutation effect of the FecB gene.

Key mRNAs and lncRNAs of pituitary that affect the reproduction of FecB + + small tail han sheep

BACKGROUND

The pituitary directly regulates the reproductive process through follicle-stimulating hormone (FSH) and luteinizing hormone (LH). Transcriptomic research on the pituitaries of ewes with different FecB (fecundity Booroola) genotypes has shown that some key genes and lncRNAs play an important role in pituitary function and sheep fecundity. Our previous study found that ewes with FecB + + genotypes (without FecB mutation) still had individuals with more than one offspring per birth. It is hoped to analyze this phenomenon from the perspective of the pituitary transcriptome.

RESULTS

The 12 Small Tail Han Sheep were equally divided into polytocous sheep in the follicular phase (PF), polytocous sheep in the luteal phase (PL), monotocous sheep in the follicular phase (MF), and monotocous sheep in the luteal phase (ML). Pituitary tissues were collected after estrus synchronous treatment for transcriptomic analysis. A total of 384 differentially expressed genes (DEGs) (182 in PF vs. MF and 202 in PL vs. ML) and 844 differentially expressed lncRNAs (DELs) (427 in PF vs. MF and 417 in PL vs. ML) were obtained from the polytocous-monotocous comparison groups in the two phases. Functional enrichment analysis showed that the DEGs in the two phases were enriched in signaling pathways known to play an important role in sheep fecundity, such as calcium ion binding and cAMP signaling pathways. A total of 1322 target relationship pairs (551 pairs in PF vs. MF and 771 pairs in PL vs. ML) were obtained for the target genes prediction of DELs, of which 29 DEL-DEG target relationship pairs (nine pairs in PF vs. MF and twenty pairs in PL vs. ML). In addition, the competing endogenous RNA (ceRNA) networks were constructed to explore the regulatory relationships of DEGs, and some important regulatory relationship pairs were obtained.

CONCLUSION

According to the analysis results, we hypothesized that the pituitary first receives steroid hormone signals from the ovary and uterus and that VAV3 (Vav Guanine Nucleotide Exchange Factor 3), GABRG1 (Gamma-Aminobutyric Acid A Receptor, Gamma 1), and FNDC1 (Fibronectin Type III Domain Containing 1) played an important role in this process. Subsequently, the reproductive process was regulated by gonadotropins, and IGFBP1 (Insulin-like Growth Factor Binding Protein 1) was directly involved in this process, ultimately affecting litter size. In addition, TGIF1 (Transforming Growth Factor-Beta-Induced Factor 1) and TMEFF2 (Transmembrane Protein With EGF Like And Two Follistatin Like Domains 2) compensated for the effect of the FecB mutation and function by acting on TGF-β/SMAD signaling pathway, an important pathway for sheep reproduction. These results provided a reference for understanding the mechanism of multiple births in Small Tail Han Sheep without FecB mutation.

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.