Effect of the Booroola fecundity (FecB) gene on the reproductive performance of ewes under assisted reproduction

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.

Genetic variation in the BLM gene and its expression in the ovaries is closely related to kidding number in goats

Bloom (BLM) helicase plays an important role in DNA replication and the maintenance of genome integrity. BLM protein deficiency, which plays a vital role in the sperm-egg union and germ-cell development during reproduction, can lead to severe DNA damage in goats. However, the effect of BLM protein deficiency on goat litter size has not been reported. Herein, we studied the association between the genetic variation in the BLM gene and the number of kids per litter in Guizhou white goats. We explored differences in the expression of the BLM protein in the follicles of single and multi-kid nanny goats. We also analyzed the effects of dysregulated BLM gene expression on the proliferation and apoptosis of ovarian granulosa cells and the expression of genes related to follicle development in goats. Five single nucleotide polymorphism (SNP) loci, including the non-synonymous mutations g.38179 A > G, g.40626 G > C and g.89621 T > G; the intron synonymous mutation g.56961 G > A and the exon synonymous mutation g.65796 C > T were found in the BLM gene. All SNPs loci were in Hardy-Weinberg equilibrium, and correlation analysis showed that the g.65796 C > T and g.89621 T > G loci polymorphism was strongly associated with litter size in the first three litters (P < 0.05). The diplogenotype Hap 2/2 (AAGGAACCTT) showed no significant difference in litter size between different births, indicating that the diploid genotype is stable in different litter sizes. Bioinformatics analysis showed that three non-synonymous mutation loci (p.T488A, p.A662S, and p.S1373A) could affect BLM protein stability, and mutations in p.T488A and p.S1373A led to changes in amino acid polarity and associated interactions. qPCR results showed that the expression level of the BLM gene in the uterus and ovaries of TT genotype nanny goats was significantly higher than that of GG genotype nanny goats. Indirect immunofluorescence assay (IF) showed that the BLM protein was significantly overexpressed in both the primordial and growing follicles of nanny goats with multiple kids (P < 0.01). Disrupting BLM gene expression in the ovarian granulosa cells down-regulated the expression of the Cyp19A1 gene. It also significantly inhibited the proliferation of follicles and induces early apoptosis of the granulosa cells. These findings confirm that polymorphism in the BLM gene is closely related to the littering traits of Guizhou white goats, and it affects the reproductive performance of nanny goats by regulating the development of the oocytes and granulosa cells. This work provides new evidence on the regulatory effect of the BLM gene on the litter size of nanny goats.

Proteomic Analysis Identifies Distinct Protein Patterns for High Ovulation in FecB Mutant Small Tail Han Sheep Granulosa Cells

The Booroola fecundity (FecB) mutation in the bone morphogenetic protein receptor type 1B (BMPR1B) gene increases ovulation in sheep. However, its effect on follicular maturation is not fully understood. Therefore, we collected granulosa cells (GCs) at a critical stage of follicle maturation from nine wild-type (WW), nine heterozygous FecB mutant (WB), and nine homozygous FecB mutant (BB) Small Tail Han sheep. The GCs of three ewes were selected at random from each genotype and consolidated into a single group, yielding a total of nine groups (three groups per genotype) for proteomic analysis. The tandem mass tag technique was utilized to ascertain the specific proteins linked to multiple ovulation in the various FecB genotypes. Using a general linear model, we identified 199 proteins significantly affected by the FecB mutation with the LIMMA package (p < 0.05). The differential abundance of proteins was enriched in pathways related to cholesterol metabolism, carbohydrate metabolism, amino acid biosynthesis, and glutathione metabolism. These pathways are involved in important processes for GC-regulated 'conservation' of oocyte maturation. Further, the sparse partial least-squares discriminant analysis and the Fuzzy-C-mean clustering method were combined to estimate weights and cluster differential abundance proteins according to ovulation to screen important ovulation-related proteins. Among them, ZP2 and ZP3 were found to be enriched in the cellular component catalog term "egg coat", as well as some apolipoproteins, such as APOA1, APOA2, and APOA4, enriched in several Gene Ontology terms related to cholesterol metabolism and lipoprotein transport. A higher abundance of these essential proteins for oocyte maturation was observed in BB and WB genotypes compared with WW ewes. These proteins had a high weight in the model for discriminating sheep with different FecB genotypes. These findings provide new insight that the FecB mutant in GCs improves nutrient metabolism, leading to better oocyte maturation by altering the abundance of important proteins (ZP2, ZP3, and APOA1) in favor of increased ovulation or better oocyte quality.

miRNA-mRNA analysis of sheep adrenal glands reveals the network regulating reproduction

Background

The adrenal gland participates in the process of sheep reproduction. MicroRNAs (miRNAs), endogenous small noncoding RNAs, regulate gene expression at the posttranscriptional level. However, the miRNA-mRNA network profile of adrenal glands relating to reproduction in sheep is still not well-studied. As sheep with FecB^BB genotype show higher lambing number compare with the sheep with FecB⁺⁺ genotype. This research aims to compare gene expression by small RNA-seq in adrenal tissues at follicular (F) and luteal (L) phases in FecB^BB (MM) and FecB⁺⁺ (ww) sheep. After analysis of gene expression, significant differentially expressed microRNAs (DEMs) and corresponding target genes were identified.

Results

A total of 180 miRNAs were found in this study, of which 19 DEMs were expressed in the four comparison groups (MM_F_A vs. MM_L_A, MM_F_A vs. ww_F_A, MM_L_A vs. ww_L_A, ww_F_A vs. ww_L_A). Subsequently, 354 target genes of 19 DEMs were predicted by integrated analysis. Cluster analysis was performed by K_means_cluster, and the expression patterns of these DEMs were separated into four subclusters. Functional analysis of target genes was performed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). The results indicated that the target genes were involved mainly in the Notch signaling pathway, signal transduction, cell communication, innate immune response and amino acid metabolism. Specifically, the Notch signaling pathway, biosynthetic process and metabolic process of pyrimidine nucleotide and amino acid metabolism appear to play key regulatory roles in the sheep fertility trait. Furthermore, miRNA-mRNA interaction networks were constructed by differentially expressed genes combined with our previous study of transcriptome data. The results showed that several key genes, including TDRD3, ANAPC7, CCNL2, BRD2 and MUT, were related to the transformation from the follicular phase to the luteal phase. PLAC8L1, NFAT5, DDX24 and MBD1 were related to the high fecundity of small tail Han sheep.

Conclusions

In this study, the miRNA transcriptome profile was identified, and miRNA-mRNA interaction networks were constructed in adrenal gland tissue of small tail Han sheep, the interaction between miR-370-3p and its targets were considered to play a major role in the reproduction regulation process. The results enriched the number of known miRNAs in adrenal glands and provided novel ideas and further information to demonstrate posttranscriptional regulation mechanisms at follicular and luteal phases in different genotypes of small tail Han sheep.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01060-y.

Are Copy Number Variations within the FecB Gene Significantly Associated with Morphometric Traits in Goats?

The Booroola fecundity (FecB) gene is a major fertility-related gene first identified in Booroola sheep. Numerous studies have investigated whether the FecB gene is a major fecundity gene in goats or whether there are other genes that play a critical role in goat fertility. Nevertheless, little attention has been paid to the role of the FecB gene in the body morphometric traits of goats, despite the positive relationship discerned between litter size and growth. We identified five copy number variations (CNVs) within the FecB gene in 641 goats, including 318 Shaanbei white cashmere (SBWC) goats, 203 Guizhou Heima (GZHM) goats, and 120 Nubian goats, which exhibited different distributions among these populations. Our results revealed that these five CNVs were significantly associated with goat morphometric traits (p < 0.05). The normal type of CNV3 was the dominant type and displayed superior phenotypes in both litter size and morphometric traits, making it an effective marker for goat breeding. Consequently, LD blocks in the region of 10 Mb upstream and downstream from FecB and potential transcription factors (TFs) that could bind with the CNVs were analyzed via bioinformatics. Although no significant LD block was detected, our results illustrated that these CNVs could bind to growth-related TFs and indirectly affect the growth development of the goats. We identified potential markers to promote litter size and growth, and we offer a theoretical foundation for further breeding work.

TIMP1 may affect goat prolificacy by regulating biological function of granulosa cells

Tissue inhibitor of metalloproteinase 1 (TIMP1) is associated with animal reproductive processes, such as follicular growth, ovulation, luteinization, and embryo development in mammals. The purposes of this study were to explore the expression and localization of TIMP1 in the ovarian tissues and determine the effect of TIMP1 on the function of granulosa cells and the association of TIMP1 with lambing-related genes of the goats. Immunohistochemical analysis showed that TIMP1 protein was strongly expressed by granulosa cells. Enzyme-linked immunosorbent assay (ELISA) results showed that TIMP1 overexpression promoted the secretion of estradiol of granulosa cells after 12, 24, and 48 h of transfection. Moreover, in vitro experiments indicated that TIMP1 had the ability to promote the cell proliferation and elevate the transcriptional levels of four genes associated with goat prolificacy, including BMPR-1B, BMP15, GDF9, and FSHB, in granulosa cells. In conclusion, TIMP1 could be an important molecule in regulating reproductive performance of the goats by affecting estrogen secretion and cell proliferation, as well as the expression of lambing-related genes of granulosa cells in the goats.

Comparative Transcriptomic Analysis of Hu Sheep Pituitary Gland Prolificacy at the Follicular and Luteal Phases

The pituitary gland directly regulates the reproduction of domestic animals. Research has increasingly focused on the potential regulatory mechanism of non-coding RNA in pituitary development. Little is known about the differential expression pattern of lncRNAs in Hu sheep, a famous sheep breed with high fecundity, and its role in the pituitary gland between the follicular phase and luteal phase. Herein, to identify the transcriptomic differences of the sheep pituitary gland during the estrus cycle, RNA sequencing (RNA-Seq) was performed. The results showed that 3529 lncRNAs and 16,651 mRNAs were identified in the pituitary gland. Among of them, 144 differentially expressed (DE) lncRNA transcripts and 557 DE mRNA transcripts were screened in the follicular and luteal phases. Moreover, GO and KEGG analyses demonstrated that 39 downregulated and 22 upregulated genes interacted with pituitary functions and reproduction. Lastly, the interaction of the candidate lncRNA XR_001039544.4 and its targeted gene LHB were validated in sheep pituitary cells in vitro. LncRNA XR_001039544.4 and LHB showed high expression levels in the luteal phase in Hu sheep. LncRNA XR_001039544.4 is mainly located in the cytoplasm, as determined by FISH analysis, indicating that XR_001039544.4 might act as competing endogenous RNAs for miRNAs to regulate LHB. LncRNA XR_001039544.4 knockdown significantly inhibited LH secretion and cell proliferation. LncRNA XR_001039544.4 may regulate the secretion of LH in the luteal-phase pituitary gland via affecting cell proliferation. Taken together, these findings provided genome-wide lncRNA- and mRNA-expression profiles for the sheep pituitary gland between the follicular and luteal phases, thereby contributing to the elucidation of the molecular mechanisms of pituitary function.

Effects of FecB Mutation on Estrus, Ovulation, and Endocrine Characteristics in Small Tail Han Sheep

The Booroola fecundity gene (FecB) has a mutation that was found to increase the ovulation rate and litter size in Booroola Merino sheep. This mutation is also associated with the fecundity of small-tail han (STH) sheep, an important maternal breed used to produce hybrid offspring for mutton production in China. Previous research showed that the FecB gene affects reproduction in STH sheep, based on litter size records. However, the effects of this gene on estrus, ovulation, and endocrine characteristics in these sheep remain unclear. Here, we analyzed the traits mentioned earlier and compared them among the three FecB genotypes of STH ewes using estrus synchronization. Overall, 53 pluriparous ewes were selected from among 890 STH ewes and subjected to FecB genotyping for experiments to characterize estrous and ovulation rates. FecB heterozygous (+B) ewes presented an earlier onset of estrus (42.9 ± 2.2 h) and a shorter estrous cycle (17.2 ± 0.2 days) (P ≤ 0.05). The ovulation rates increased with the increasing copy number of the B allele (P ≤ 0.01). Ovulation time showed no significant differences among the three FecB genotypes. The serum concentrations of follicle-stimulating hormone (FSH), luteinizing hormone, estrogen (E2), and progesterone (P4) were measured in 19 of the ewes. Serum concentrations of E2 and FSH dramatically varied around the time of behavioral estrus. In FecB mutant homozygous (BB) ewes, E2 concentration had two peaks, which were higher (P ≤ 0.05) than those of ++ genotypes. FSH concentration of BB ewes was higher (P ≤ 0.05) than that of the ++ ewes just after estrus. The expression of the estrogen receptor 1 (ESR1) gene in the +B genotype was higher than in the other genotypes. Based on the data for the reproductive performance of STH ewes with the three FecB genotypes, our study suggests that the development of follicles in ewes with the B allele is dependent on the response to FSH regulated by E2 in the early stage. +B ewes, exhibiting moderate ovulation and litter size and a shorter estrous cycle, can be highly recommended in sheep crossbreeding systems for commercial mutton production. Moreover, this study provides useful information to conserve better and use the genetic resources of STH sheep in China.

Detection of mRNA Expression and Copy Number Variations Within the Goat Fec B Gene Associated With Litter Size

The Booroola fecundity (Fec ^B ) gene, as the first major fecundity gene identified in Booroola sheep, has attracted careful attention. So far, previous research have uncovered the Fec^B mutation (Q249R) as the main mutation by virtue of which sheep exhibits multiple lambing phenomena. This mutation is now being intensively studied and widely used. However, such effect of the Fec^B mutation has not been applied to goats, and similar types of the Fec ^B gene in goats still need to be studied. Thus, the current study attempted to verify potential mutations in the goat Fec ^B gene as well as investigate their functions related to fecundity. First, Fec ^B expression was investigated in six different goat tissues, and we found that Fec ^B expression was highest in the mammary gland, followed by the ovary. Next, the influence of the Fec ^B gene was analyzed from a new perspective, where five potential copy number variations (CNVs) (CNV1-5) within the Fec ^B gene were identified for the first time, and then their effects on litter size were measured. Our results point out that CNV3 (P = 3.44E-4) and CNV5 (P = 0.034) could significantly influence the litter size of goats. Identically, the combination genotype of CNV3 and CNV5 which consisted of their dominant genotypes was also significantly associated with goat litter size (P = 7.80E-5). Hence, CNV3 and CNV5 could serve as potential DNA molecular markers applied to DNA editing and DNA microarray. Additionally, the abovementioned study has laid a theoretical foundation for the detection of potential fertility-related quantitative trait loci within the goat Fec ^B gene.

Survey of the relationship between polymorphisms within the BMPR1B gene and sheep reproductive traits

The BMPRIB gene is one of the main genes that can be used as a molecular genetic marker for the early selection of highly productive ewes. It is well-documented that the p.Q249R (g.746A > G) is the first mutation in the kinase domain of the BMPR1B gene that is highly related to increased ovulation rate and litter size. It is likely that the presence of the p.Q249R mutation in the sheep population is one of the factors contributing to the outstanding productivity of the sheep. Moreover, in recent years, researchers have been explored other polymorphisms in the BMPR1B gene with respect to reproductive traits in sheep. Therefore, we carried out the current study to evaluate the association between polymorphisms in this gene and sheep litter size from all appropriate studies. As a result, among 41 polymorphisms in the ovine BMPRIB gene, eight variants, including p.Q249R (g.746A > G), g.29362047T > C, g.29427689G > A, BMPR1B-2 (ss:1960972599), g.29382337G > A, g.29382340G > A, rs1092293287 (10 bp insertion/deletion) and g.29380965A > G were found to be associated with litter size in sheep. This systematic analysis presents the most current data evidence for BMPRIB polymorphisms, highlighting the need for further large-scale studies to determine more important variants.

Combined approaches identify known and novel genes associated with sheep litter size and non-seasonal breeding

Improvement of ewe reproduction is considerable by appropriately increasing litter size and sustaining non-seasonal breeding. However, their genetic makeups have not been entirely elucidated. Genome-wide analyses of 821 individuals were performed by combining three genomic approaches (genome-wide association study, XP-nSL, and runs of homozygosity). Consequently, 35 candidate genes including three domestication genes (TSHR, GTF2A1, and KITLG) were identified. Other than the FecB mutation at BMPR1B, we described a significant association of a missense mutation rs406686139 at seasonal lambing-associated TSHR gene with litter size. Some promising novel genes may be relevant for sheep reproduction by multitude biological processes, such as FETUB functioning in fertilization, HNRNPA1 in oogenesis, DCUN1D1 in spermatogenesis, and HRG in fertility outcome. The present study suggests that improvement of ewe reproduction is attributed to selective breeding, and casts light on the genetic basis and improvement of sheep reproduction.

The ovarian follicle of ruminants: the path from conceptus to adult

This review resulted from an international workshop and presents a consensus view of critical advances over the past decade in our understanding of follicle function in ruminants. The major concepts covered include: (1) the value of major genes; (2) the dynamics of fetal ovarian development and its sensitivity to nutritional and environmental influences; (3) the concept of an ovarian follicle reserve, aligned with the rise of anti-Müllerian hormone as a controller of ovarian processes; (4) renewed recognition of the diverse and important roles of theca cells; (5) the importance of follicular fluid as a microenvironment that determines oocyte quality; (6) the 'adipokinome' as a key concept linking metabolic inputs with follicle development; and (7) the contribution of follicle development to the success of conception. These concepts are important because, in sheep and cattle, ovulation rate is tightly regulated and, as the primary determinant of litter size, it is a major component of reproductive efficiency and therefore productivity. Nowadays, reproductive efficiency is also a target for improving the 'methane efficiency' of livestock enterprises, increasing the need to understand the processes of ovarian development and folliculogenesis, while avoiding detrimental trade-offs as greater performance is sought.

Transcriptome Analysis Reveals Differentially Expressed Genes and Long Non-coding RNAs Associated With Fecundity in Sheep Hypothalamus With Different FecB Genotypes

Small-tailed Han sheep, with different FecB genotypes, manifest distinct ovulation rates and fecundities, which are due to differences in reproductive hormones secreted by the hypothalamic-pituitary-ovarian axis. Nevertheless, the function of the hypothalamus against a FecB mutant background on increasing ovulation rate is rarely reported. Therefore, we determined the expression profiles of hypothalamus tissue collected from six wild-type (WW) and six FecB mutant homozygous (BB) ewes at the follicular and luteal phases by whole-transcriptome sequencing. We identified 53 differentially expressed mRNAs (DEGs) and 40 differentially expressed long non-coding RNAs (DELs) between the two estrus states. Functional annotation analysis revealed that one of the DEGs, PRL, was particularly enriched in the hypothalamic function, hormone-related, and reproductive pathways. The lncRNA-target gene interaction networks and KEGG analysis in combination suggest that the lncRNAs LINC-676 and WNT3-AS cis-acting on DRD2 and WNT9B in different phases may induce gonadotropin-releasing hormone (GnRH) secretion. Furthermore, there were differences of regulatory elements and WNT gene family members involved in the follicular-luteal transition in the reproductive process between wild-type (WNT7A) and FecB mutant sheep (WNT9B). We combined the DEG and DEL data sets screened from different estrus states and genotypes. The overlap of these two sets was identified to select the mRNAs and lncRNAs that have major effects on ovulation. Among the overlapping molecules, seven DEGs and four DELs were involved in the follicular-luteal transition regulated by FecB mutation. Functional annotation analysis showed that two DEGs (FKBP5 and KITLG) were enriched in melanogenesis, oxytocin, and GnRH secretion. LINC-219386 and IGF2-AS were highly expressed in the BB ewes compared with WW ewes, modulating their target genes (DMXL2 and IGF2) to produce more GnRH during follicular development, which explains why mutated ewes produced more mature follicles. These results from expression profiling of the hypothalamus with the FecB mutation at different estrus states provide new insights into how the hypothalamus regulates ovulation under the effect of the FecB mutation.

A novel (ATC)n microsatellite locus is associated with litter size in an indigenous Chinese pig

Simple sequence repeats (SSRs) are an important part of the genome and have become powerful auxiliary DNA markers in animal breeding using marker-assisted selection (MAS). Based on previous sequencing data of Qinghai Bamei pigs, a total of three novel candidate SSR loci were analysed in this study. Time-of-flight mass spectrometry (TOF-MS) was used for SSR genotyping, and association analyses between SSRs and the litter size of Qinghai Bamei sows was also performed. The results of genotyping showed that the (ATC)_n -P1, (AC)_n -P2 and (AC)_n -P3 loci had 2, 3 and 18 genotypes, respectively; 2, 3 and 8 alleles were also identified at these loci. Except for the (AC)_n -P2 locus, the polymorphism information content (PIC) values of other loci were greater than 0.25. Association analyses indicated that only the (ATC)_n -P1 locus was significantly associated with the litter size of Qinghai Bamei sows (p = .047). Compared to 189-/189- genotype, individuals with the 189-/195- genotype had the senior litter size, which was 9.04 ± 0.21. Our results enrich the data on SSRs in Qinghai Bamei pigs and indicate that (ATC)_n -P1 is a candidate locus for MAS in the pig industry.

Deletion mutation within the goat PPP3CA gene identified by GWAS significantly affects litter size

The protein phosphatase 3 catalytic subunit α (PPP3CA) gene is a high reproduction traits candidate gene for goats as revealed by a genome-wide association study. The aim of this work was to explore the genetic variations of the goat PPP3CA as well as to evaluate the genetic effects on litter size. Three novel insertions/deletions (indels) within the goat PPP3CA were found and their minor allelic frequencies (MAF) were 0.105, 0.066, and 0.042, respectively. The results showed that only the 20bp indel polymorphism was significantly associated with litter size in Shaanbei white cashmere goats (P<0.05) and individuals with deletion/deletion (DD) genotypes demonstrated the junior phenotypes when compared with those with other genotypes. These findings suggested that the 20bp indel is a potential DNA marker for selecting superior individuals in marker-assisted selection for breeding concerning fecundity in goats.