Association analysis of single nucleotide polymorphisms in the LIFR gene with lambing number in sheep

Lambing number trait in sheep is a complex trait controlled by multiple genes, and low lambing number has a severe economic impact on the sheep industry. Previous studies have shown that Sparc/Osteonectin, cwcv, and kazal-like domain proteoglycan 1 (SPOCK1), A disintegrin and metalloproteinase with thrombospondin-like repeats 1 (ADAMTS1), heparin-binding epidermal growth factor-like growth factor (HBEGF), and leukemia inhibitory factor receptor (LIFR) are involved in mammalian reproduction. However, the effects of these genes on lambing number in sheep are still unclear. In this study, single nucleotide polymorphism (SNP) loci at the above four genes were genotyped in five sheep breeds (two single-born sheep breeds and three multiple-born sheep breeds with a total of 768 sheep) using the Sequenom MassARRAY^® SNP assay, and their associations with the lambing number in small-tailed Han sheep were analyzed. The results showed that a total of six SNP loci (c.^*1633A>G, c.^*1388T>C, c.^*1095G>A, c.1847A>G, c.^*2403A>G, and c.^*127T>C) were identified in the four genes, SPOCK1, ADAMTS1, HBEGF, and LIFR. All the above SNPs had three genotypes in five sheep breeds. The population genetic analysis of SNPs of the four genes and the association analysis with lambing number showed that the polymorphism information content of the five sheep breeds of small-tailed Han sheep, Hu sheep, Cele black sheep, Sunite sheep, and Bamei mutton sheep was between 0 and 0.37, and a few sheep breeds were in Hardy-Weinberg equilibrium (p>0.05). The c.^*127 T>C locus of the LIFR gene may be affected by natural or artificial selection in these sheep breeds. Additionally, the association analysis between the c.^*127 T>C locus of LIFR gene and the lambing number of small-tailed Han sheep showed that the c.^*127 T>C locus of the LIFR gene was significantly associated with the lambing number of the second, third, and average parity of smalltailed Han sheep (p < 0.05). The lambing number of ewes with the LIFR CC genotype was significantly lower than that with the TT and CT genotypes (p < 0.05). The LIFR protein interaction network was also predicted and found to interact with the reported ciliary neurotrophic factor (CNTF), cardiotrophinlike cytokine factor 1 (CLCF1), interleukin 6 cytokine family signal transducer (IL6ST), and ciliary neurotrophic factor receptor (CNTFR) proteins. In conclusion, the c.^*127 T>C locus of LIFR gene can be used as a candidate genetic molecular marker to increase lambing numbers in polytocous sheep.

1 Introduction

Sheep are important economic animals as they provide valuable agricultural products to humans. Increasing the reproductive performance of sheep is one of the most effective ways to enhance the sheep industry and increase its economic benefits (1). Traditional breeding methods face challenges in rapidly increasing lambing numbers due to the low heritability. Compared to traditional breeding methods, molecular marker-assisted selection breeding methods significantly shorten the breeding cycle and improve the success rate of breeding (2). Therefore, molecular marker-assisted selection breeding is significant in improving the lambing number of sheep. In particular, a few studies have found that general transcription factor IIA subunit 1 (GTF2A1) (3), GATA binding protein 4 (GATA4) (4), growth differentiation factor 9 (GDF9) (5), mutations at loci such as prolactin receptor (PRLR), insulin like growth factor 1 (IGF1), and lymphoid enhancer-binding factor (LEF) (6) have been associated with reproductive traits in sheep, leading to differences in lambing numbers. Assisted selection using these molecular markers can improve the reproductive efficiency of sheep breeds, leading to higher productivity and improved economic returns.

Sparc/Osteonectin, cwcv, and kazal-like domain proteoglycan 1 (SPOCK1) gene inhibits the apoptosis of tumor cells and is located in the nucleus and mitochondria of tumor cells. However, the physiological function of this gene has not been determined (7). It has been reported that SPOCK1 can regulate cell proliferation, migration, and apoptosis, and promote the development of several cancers (8, 9). Based on genome-wide association study (GWAS) correlation analyses, SPOCK1 was identified as a new candidate gene in beef cattle (10) and human (11) menstrual cycle, suggesting that SPOCK1 gene could be involved in animal reproductive life (12). It has been found that SPOCK1 inhibits matrix metalloproteinase-2 (MMP-2) activation, which is a key factor in promoting endometrial menstrual rupture and the initiation of menstrual bleeding (13, 14). SPOCK1 gene associated with female reproductive tract traits in polytocous sheep has the potential to regulate fertility and reproduction in the species (15). A disintegrin and metalloproteinase with thrombospondin-like repeats 1 (ADAMTS1) is located in granulosa cells (16) and plays an important role in the ovary through the progesterone receptor (PR)-dependent pathway (17). It is upregulated in preovulatory follicles of humans, mice, sheep, and cattle (18–20). Heparin-binding epidermal growth factor-like growth factor (HBEGF) inhibits apoptosis through the phosphatidylinositol-3 kinase/serine/threonine kinase (PI3K/Akt) signaling pathway (21), which is involved in early placental development and protects trophoblast cells from hypoxia-induced damage during pregnancy (22). Studies have also shown that HBEGF can be well-positioned to transmit paracrine near-secretory signals to blastocysts, which is related to the activation of uterine blastocysts during pregnancy (23, 24). HBEGF is a key regulator of embryo implantation in rodents, humans, and non-primates (22, 25). These results suggest that SPOCK1, ADAMST1, and HBEGF are closely related to the regulation of ovarian function, embryo implantation, and hormone synthesis, and have an effect on lambing numbers in sheep.

Leukemia inhibitory factor receptor (LIFR) belongs to the family of type I cytokine receptors and regulates cell proliferation and differentiation by activating the Janus Kinase/Signal Transducer and Activator of Transcription (JAK/STAT) signaling pathway. Mice deficient in LIFR show failure of implantation and metamorphosis, leading to embryonic death, suggesting that both the embryo and the uterus are potential targets for LIFR signaling (26, 27). Studies on the association between leukemia inhibitory factor (LIF) and lambing number in pigs have shown that LIFR plays an important role in embryo implantation (28). LIF is a key factor in regulating embryo implantation in many mammals (29). LIFR is expressed in the endometrium, preimplantation uterus, and placenta during embryo implantation, as well as in later stages of development (30). Although there have been limited studies on LIFR in sheep, it plays a crucial role in cell proliferation and embryonic development, affects reproduction, and can influence reproductive traits in sheep.

The above four candidate genes (SPOCK1, ADAMTS1, HBEGF, and LIFR) play a potential role in reproduction and so could be involved in sheep reproduction, resulting in differences in lambing number (15–18, 22, 29). SNPs in the above four genes were also found to be associated with reproductive traits in sheep in the pre-sequencing data. The obtained SNP annotation results and Fixation Index (Fst) were integrated and screened according to the following conditions: SNPs in introns with Fst > 0.15; non-synonymous mutation with Fst >0.05. In the present study, Sequenom MassARRAY^® SNP technology was used to detect the distribution of these four genes in five sheep breeds, categorized as polytocous and monotocous, and population genetics analysis was performed. Subsequently, bioinformatics analysis tools were used to predict the tertiary structure of candidate genes and their interacting proteins. The results provide new ideas and valuable genetic markers for sheep breeding.

2 Materials and methods

2.1 Animal preparation and sample collection

Ewes of five sheep breeds were included in this study, three of which were polytocous breeds, namely, small-tailed Han sheep (n = 384), Cele black sheep (n = 96), and Hu sheep (n = 96), and two were monotocous breeds, namely, Sunite sheep (n = 96), and Bamei mutton sheep (n = 96). The details of the five sheep breeds are displayed in Table 1. Blood samples were collected from all the sheep via the jugular vein, anticoagulated with EDTA, and stored at −20°C. The blood samples were collected at the same time from all the sheep that were uniformly bred, and details of the lambing numbers of small-tailed Han Sheep up to 3 years of age were recorded (refer to Supplementary Table S1). The other four breeds in this experiment were either all monotocous or polytocous breeds. However, the small-tailed Han sheep population is unique in that it encompasses individuals that produce both single and multiple lambs. Given this diversity, we have used the four breeds with established reproductive traits as reference points, and the small-tailed Han sheep with specific lambing records as our validation cohort. This strategy enables us to pinpoint potential genetic variations associated with lambing number, even when direct phenotypic data for all breeds are not available.

Table 1

Table 1. Information on ewes used in this study.

2.2 DNA extraction

DNA extraction from the sheep blood was performed using a DNA extraction kit (Tiangen Biochemical Technology Co., Ltd., Beijing, China). The steps are described in the instructions of the “Blood/Cell/Tissue Genomic DNA Extraction Kit” (catalog number: DP304). The concentration of the extracted DNA samples was determined using a Nanodrop 2000 (Thermo Fisher, USA), and the quality of the extracted DNA samples was evaluated by 1.5% agarose gel electrophoresis.

2.3 Genotyping

According to the sequence of sheep SPOCK1, ADAMTS1, HBEGF, and LIFR in GenBank ARS-UI_Ramb_v2.0 (accession number: XM_042250908.1, XM_004002797.5, XM_004003101.6, and XM_060400477.1), the single base extension primers of four loci (c.^*1633 A>G, c.^*1388T>C, and c.^*1095G>A in SPOCK1, c.1847 A>G in ADAMTS1, c.^*2403 A>G in HBEGF, and c.^*127 T>C in LIFR) were designed by MassARRAY Assay Design v.3.1. The gene primers used in this study are shown in Table 2. These primers were synthesized by Beijing Compass Biotechnology Co., Ltd. (Beijing, China). Finally, Sequenom MassARRAY^® SNP technology was used to genotype the four loci of the polytocous groups (small-tailed Han sheep, Cele black sheep, and Hu sheep) and the monotocous groups (Sunite sheep and Bamei mutton sheep). For detailed information about the system and program, we referred to the previous research of our laboratory (31).

Table 2

Table 2. Primer sequence information of the Sequenom MassARRAY^® SNP genotyping.

2.4 Polymerase chain reaction amplification and sanger sequencing

According to the results of genotyping and lambing number analysis of small-tailed Han sheep, the LIFR gene on chromosome 16 (Oar_v3.1) of sheep was selected for Polymerase Chain Reaction (PCR). Three DNA samples were randomly selected from each genotype for Sanger sequencing. Based on the annotation of LIFR SNP rs162302066 (Oar_v3.1) in the Ensembl gene annotation database, a 505-bp sequence was selected as the center to design primers for PCR amplification. The primer sequence is shown in Table 3. PCR amplification products were detected using agarose gel electrophoresis before Sanger sequencing. The primer amplification and Sanger sequencing were performed by Beijing Sangon Biotech Co., Ltd. The PCR reaction system included 8.91 μl of double-distilled water (ddH₂O), 1 μl each of upstream and downstream primers (10 μM), 1.59 μl of template DNA ( ≤ 1 μg), and 12.5 μl of 2 × Master Mix. The amplification cycle included 30 cycles of pre-denaturation at 94°C for 3 min, followed by denaturation at 94°C for 30 s, annealing at 55°C for 30 s, extension at 72°C for 1 min, and final extension at 72°C for 5 min.

Table 3

Table 3. Sequence information of PCR amplification primer for the LIFR gene.

2.5 Statistical analysis

The formulas for calculating allele and genotype frequencies, polymorphism information content (PIC), heterozygosity (He), and effective number of alleles (Ne) using genotyping data are as follows (32):

$\begin{array}{l}PIC=1-{\displaystyle \sum _{i=1}^n{p}_i^2}-{\displaystyle \sum _{i=1}^{n-1}{\displaystyle \sum _{j=i+1}^{n}2p_i^2{p}_j^2}}\\ HE=1-{\displaystyle \sum _{i=1}^n{p}_i^2}\\ NE=1/{\displaystyle \sum _{i=1}^n{p}_i^2}\end{array}$

where n is the number of alleles, p_i is the allele frequency of the ith allele, and p_j is the allele frequency of the jth allele.

The chi-squared test was used to detect whether the genotype distribution of each locus was in Hardy–Weinberg equilibrium. SPSS 26 (analysis of variance) was used to analyze the correlation between the four genotypes and lambing number. The following adjusted linear model was used:

$\begin{array}{l}y=\text{μ }+P+G1+G2+G3+G4+e\end{array}$

The least square mean of multiple comparisons of lambing number among different genotypes of small-tailed Han sheep was used. In the above equation, y is the phenotypic value (lambing number), μ is the population mean, P is the fixed parity effect, G1 is the fixed effect of candidate SNP in SPOCK1, G2 is the fixed effect of candidate SNP in ADAMTS1, G3 is the fixed effect of candidate SNP in HBEGF, G4 is the fixed effect of candidate SNP in LIFR, and e is the random error effect of each observation. At the same time, we also examined the interaction effect between each loci, in order to check whether the interaction between candidate loci needs to be supplemented in the model. However, no significant interaction was found between the candidate loci. Analysis was performed using R software (avo, Version 4.0.3).

2.6 Protein interaction network predicted by STRING database

For SNP loci significantly associated with lambing number, we predicted the protein interaction networks LIFR from the STRING database v.12.0 (https://cn.string-db.org) (33), which is a useful tool for predicting protein interactions.

3 Results

3.1 Genotyping and population genetic analysis of SNPs in candidate genes SPOCK1, ADAMTS1, HBEGF, and LIFR

The six SNPs (c.^*1633 T>C, c.^*1388T>C, and c.^*1095G>A in SPOCK1, c.^*1847 A>G in ADAMTS1, c.^*2403 T>C in HBEGF, and c.^*127 T>C in LIFR) were genotyped in five sheep breeds. The results showed that there were three genotypes in the SNP loci of the four genes, which could be visually distinguished by the MassARRAY^® SNP system (Figure 1).

Figure 1

Figure 1. Genotyping results of candidate SNPs in SPOCK1, ADAMTS1, HBEGF, and LIFR genes using MassARRAY^® SNP system: (A) c.*1633T>C of SPOCK1, (B) c.*1388T>C of SPOCK1, (C) c.*1095G>A of SPOCK1, (D) c.1847A>G of ADAMTS1, (E) c.*2403T>C of HBEGF, and (F) c.*127T>C of LIFR.

The results of population genetic analysis of five sheep breeds (refer to Table 4) showed that the c.^*1633 T>C and c.^*1388T>C loci of SPOCK1 gene showed low polymorphism in small-tailed Han sheep and Bamei mutton sheep (PIC <0.25), and moderate polymorphism in Cele black sheep, Hu sheep, and Sunite sheep (0.25 ≤ PIC <0.5). The c.^*1095G>A locus of the SPOCK1 gene showed low polymorphism in Bamei mutton sheep (PIC <0.25), and moderate polymorphism in small-tailed Han sheep, Cele black sheep, Hu sheep, and Sunite sheep (0.25 ≤ PIC <0.5). The c.1847A>G locus of the ADAMTS1 gene showed low polymorphism in five sheep breeds (PIC <0.25). The c.^*2403 T>C locus of the HBEGF gene showed low polymorphism (PIC <0.25) in small-tailed Han sheep, Hu sheep, Sunite sheep, and Bamei mutton sheep, and moderate polymorphism (0.25 ≤ PIC <0.5) in Cele black sheep. The c.^*127 T>C locus of the LIFR gene showed moderate polymorphism in five sheep breeds (0.25 ≤ PIC <0.5). According to the Hardy–Weinberg equilibrium, the c.^*1633 T>C, c.^*1388T>C, and c.^*1095G>A loci of SPOCK1, the c.1847 A>G locus of ADAMTS1, and the c.^*2403 T>C locus of HBEGF were in Hardy–Weinberg equilibrium in five sheep breeds (p > 0.05). The c.^*127 T>C locus of the LIFR gene was in Hardy–Weinberg equilibrium in small-tailed Han sheep, Cele black sheep, Hu sheep, and Bamei mutton sheep (p > 0.05), and deviated from Hardy–Weinberg equilibrium in Sunite sheep (p < 0.05).

Table 4

Table 4. Population polymorphism analysis of four candidate loci in five sheep breeds.

3.2 Genotyping and allele frequencies of SNPs in the four genes of polytocous and monotocous sheep

According to the lambing number, the five breeds were divided into two categories: polytocous sheep (small-tailed Han sheep, Cele black sheep, and Hu sheep) and monotocous sheep (Sunite sheep and Bamei mutton sheep). Genotyping of six SNPs in the five sheep breeds is shown in Table 5. There was no significant difference in genotype frequency and allele frequency of c.^*1633 T>C locus and c.^*1388T>C locus of the SPOCK1 gene, c.1847 A>G locus of the ADAMTS1 gene, and c.^*2403 T>C locus of the HBEGF gene between polytocous and monotocous sheep breeds (p > 0.05). The genotype and allele frequencies of c.^*1095G>A locus of the SPOCK1 gene and c.^*127 T>C locus of the LIFR gene were significantly different between polytocous and monotocous sheep breeds (p < 0.05).

Table 5

Table 5. Genotype and allele frequencies of SNP loci in populations of sheep with different reproductive capacities.

3.3 Association analysis of SNPs in SPOCK1, ADAMTS1, HBEGF, and LIFR with lambing number in small-tailed Han sheep

The association analysis of SNPs in SPOCK1, ADAMTS1, HBEGF, and LIFR genes with lambing number in Small-tailed Han sheep is shown in Table 6. There was a significant association between the c.^*127 T>C locus of LIFR and the second, third, and average lambing number of small-tailed Han sheep (p < 0.05). The lambing number of ewes with the CC genotype was significantly lower than that with the TT and CT genotypes (p < 0.05). The results showed that the mutation site of LIFR c.^*127 T>C could be involved in the regulation of lambing number in small-tailed Han sheep to a certain extent, while the other loci were not significantly associated with the lambing number of small-tailed Han sheep (p > 0.05).

Table 6

Table 6. Least squares mean and standard errors of the lambing number in small-tailed Han sheep with different genotypes.

3.4 SNP site sequencing of the LIFR candidate gene

Association analysis between candidate genes and the lambing number of small-tailed Han sheep showed that there was a significant difference between the LIFR candidate gene and lambing number in this breed. Genotyping of LIFR SNP c.^*127T>C revealed the presence of three genotypes of this locus in the small-tailed Han sheep population, namely, TT (209, 54.4%), CT (147, 38.3%), and CC (28, 7.3%), with a detection rate of >99% in 384 samples (Figure 1). The results of PCR amplification and Sanger sequencing further illustrated the existence of LIFR SNP c.^*127T>C and the reliability of the typing results (Figure 2).

Figure 2

Figure 2. Electrophoresis results of PCR products and sequencing results for the LIFR SNP c.*127T>C locus.

3.5 Structure and interaction network analysis of the LIFR protein

To further explore the biological function of the LIFR gene in reproductive traits in sheep, a network of closely related proteins was constructed using the STRING website. Among the 10 proteins that interacted most strongly with LIFR (Figure 3), the key proteins that were closely related to lambing number, such as ciliary neurotrophic factor (CNTF), cardiotrophin-like cytokine factor 1 (CLCF1), interleukin 6 cytokine family signal transducer (IL6ST), and ciliary neurotrophic factor receptor (CNTFR), were identified based on interacting protein information (refer to Table 7) and were found to be closely related to reproduction, such as development and embryo implantation (34–36). These findings further suggest that LIFR is an important factor influencing lamb numbers in sheep.

Figure 3

Figure 3. Protein networks closely related to LIFR in sheep.

Table 7

Table 7. Information on the 10 proteins that interact most strongly with LIFR.

4 Discussion

We previously screened four candidate genes (SPOCK1, ADAMTS1, HBEGF, and LIFR) associated with lambing number in sheep by GWAS analysis and selective signal scanning (15, 22, 37). Their association with lambing number in sheep needs to be further clarified. In this study, the four candidate genes mentioned above were validated by amplification using Sequenom MassARRAY^® SNP technology.

In this study, specifically, mutation loci including c.^*1633T>C, c.^*1388T>C, and c.^*1095G>A for SPOCK1, c.1847A>G for ADAMTS1, and c.^*2403T>C for HBEGF exhibited low polymorphism and moderate polymorphism. This suggests that the level of genetic variation of these loci varies in different breeds, which may have an impact on the genetic diversity and adaptability of the breed. On the other hand, the c.^*127T>C mutation of LIFR showed moderate polymorphism, suggesting that the level of genetic variation at this locus is relatively high. This polymorphism may contribute to the phenotypic diversity observed in sheep populations. Meanwhile, analyses of Hardy–Weinberg equilibrium revealed interesting patterns. The c.^*127T>C mutation in LIFR was found to significantly deviate from Hardy–Weinberg equilibrium (p < 0.05) in a population of Sunite sheep, suggesting that certain evolutionary forces, such as selection, mutation, migration, or non-random mating, may act on these loci. In contrast, c.^*1633T>C, c.^*1388T>C, and c.^*1095G>A of SPOCK1, c.1847A>G of ADAMTS1, and c.^*2403T>C of HBEGF were all found to be in Hardy–Weinberg equilibrium. Additionally, the c.^*127T>C mutation of LIFR in the other four sheep breeds was also in Hardy–Weinberg equilibrium. This suggests that these loci are relatively stable under current population conditions and are not subject to strong evolutionary pressures. In addition, we found that the genotype frequencies and allele frequencies of the c.^*1095G>A locus of the SPOCK1 gene and the c.^*127 T>C locus of the LIFR gene differed significantly (p < 0.05) between polytocous and monotocous sheep. This implies that these loci are significantly associated with lambing traits. In contrast, the genotype frequencies and allele frequencies of the c.^*1633 T>C locus and c.^*1388 T>C locus of the SPOCK1 gene, the c.1847 A>G locus of the ADAMTS1 gene, and the c.^*2403 T>C locus of the HBEGF gene were not significantly different between multiparous and singleton sheep (p > 0.05), which may indicate that these loci are not significantly associated with lambing traits in sheep. In addition to genetic analyses of the identified SNPs, a review of the existing literature on LIFR function provided valuable insights.

In the present study, we also carried out an association analysis of the SNP loci of SPOCK1, ADAMTS1, HBEGF, and LIFR genes with the lambing number of small-tailed Han sheep and found that the c.^*127 T>C locus of the LIFR gene was also significantly correlated with the second and third lambing numbers and the average lambing number of the small-tailed Han sheep. The number of lambs in ewes with TT and CT genotypes was significantly higher than that in ewes with the CC genotype. This is in agreement with previous findings (28), which suggests that LIFR may be a potential genetic marker influencing the number of lambs in sheep reproduction. SNPs in the LIFR gene were successfully identified by Sanger sequencing. The presence of LIFR SNPs in the small-tailed Han sheep population was established, indicating the presence of the LIFR SNP c.^*127 T>C and the reliability of the typing results. In addition, as a transmembrane receptor of leukemia inhibitory factor (LIF), LIFR is widely considered to be a key factor in regulating mammalian embryo implantation (29). Moreover, LIF can inhibit the differentiation of mouse embryonic stem cells (ES), suggesting that this cytokine regulates uterine function and embryonic development. Histologically, the lack of well-organized glial cells and labyrinthine layers in the placenta suggests that embryonic lethality is due to changes in trophoblast cell differentiation and that LIFR signaling plays an important role in promoting embryonic trophoblast cell differentiation (38, 39). Moreover, LIFR transmits LIF signals through downstream channels (40). Disruption of LIFR results in decreased embryo survival. LIFR may be a ligand for the activation of Signal Transducer and Activator of Transcription 3 (STAT3), which has been shown to play a role in mammary degenerative apoptosis (41). LIF/STAT3 signaling is consistent with the induction of apoptosis in mammary epithelial cells (42), suggesting that increased expression of LIFR in granulosa cells leads to increased STAT3 phosphorylation. Normal embryo implantation and development, as well as follicular atresia, have an important impact on mammalian fertility and are also key factors affecting the number of lambs produced by the small-tailed Han sheep. This further suggests that LIFR is involved in the reproductive process of sheep and can influence the number of lambs produced by sheep by affecting embryo development, which is consistent with the result of this study that LIFR can affect reproductive traits in sheep.

The majority of physiological processes are influenced by multiple proteins. Therefore, to better understand the role of LIFR genes in reproductive traits in sheep, a protein–protein interaction network was established to deepen our understanding of complex traits such as reproduction. Based on the protein–protein interaction network provided by the STRING database, we predicted 10 proteins that interact closely with LIFR. The proteins Janus kinase (JAK1), Janus kinase 2 (JAK2), and Janus kinase 3 (JAK3) are important members of the JAK family, which are involved in the signaling of a variety of cytokines and growth factors, and are essential for immune response, hematopoiesis, and metabolic regulation. Phosphoinositide-3-kinase regulatory subunit 3 (PIK3R3) is a key regulator of the PI3K signaling pathway, affecting cell growth and survival. However, four proteins, CNTF, CLCF1, IL6ST, and CNTFR, were related to reproduction, such as development and embryo implantation (35, 43, 44). Furthermore, all belong to the interleukin 6 (IL6) family (45, 46). Its stimulation by interferon tau (IFNT) and regulation by Progesterone 4 (P4) in a complex physiological phase and cell-specific manner is associated with the acceptability of intrauterine implantation in many species, including humans, rodents, pigs, cattle, mink and sheep (45, 47). CNTF plays an important role in the production of luteinizing hormone (LH) and prolactin (PRL) before ovulation in rats. The absence of CNTF will affect the production of LH, which in turn affects the ovulation rate (48). IL6 mediates early embryo implantation and placenta development through ovarian steroid hormones, estrogen, and progesterone (49). It indicates that CNTF, CLCF1, IL6ST, and CNTFR are involved in the regulation of reproductive processes. In conclusion, LIFR genes are associated with the regulation of lambing numbers in sheep. Thus, these networks further confirm that the c.^*127 T>C locus of the LIFR gene can be used as a potential molecular marker for the lambing number trait in small-tailed Han sheep. However, the interactions between these genes and how they synergistically regulate embryo implantation, follicular development, and ovulation remain to be further explored. It was also found that LIFR can influence lambing number in sheep through its function in embryonic development. the significant correlation between the c.^*127 T>C locus of the LIFR gene and lambing number suggests that this SNP could serve as a valuable molecular marker in sheep breeding programs, thereby improving the economics of sheep production.

5 Conclusion

In this study, the association between SNPs of SPOCK1, ADAMTS1, HBEGF, and LIFR and lambing number in small-tailed Han sheep was analyzed. The results showed that the genotype frequency and allele frequency of the LIFR gene c.^*127 T>C locus were significantly different between polytocous and monotocous sheep breeds. There was a significant association between the c.^*127T>C locus of the LIFR gene and lambing number in small-tailed Han sheep. It is speculated that there may be a correlation between this locus and the lambing number trait in sheep. Therefore, this gene locus can be used as a potential genetic marker for sheep breeding.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

All animals were authorized by the Department of Scientific Research (responsible for animal welfare), Institute of Animal Science, Chinese Academy of Agricultural Sciences (IAS-CAAS; Beijing, China). In addition, the IAS-CAAS Animal Ethics Committee approved the ethics of animal survival (No. IAS 2020-63) on 27 April 2020. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

YW: Data curation, Formal analysis, Methodology, Supervision, Validation, Writing – review & editing. YX: Investigation, Methodology, Validation, Visualization, Writing – original draft. RZ: Software, Writing – original draft. KL: Software, Writing – original draft. YL: Conceptualization, Funding acquisition, Project administration, Writing – review & editing. MC: Conceptualization, Funding acquisition, Project administration, Resources, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was funded by National Natural Science Foundation of China (32302719), Central Public-interest Scientific Institution Basal Research Fund (Y2024YJ08), Science Projects in Henan Province (242102111009), Young Talent Lifting Project in Henan Province (2023HYTP026), Postdoctoral Fellowship Program of CPSF (GZC20233053), Science and Technology Achievement Demonstration Project of Xinjiang Uygur Autonomous Region (2022NC116), Agricultural Science and Technology Innovation Program of China (ASTIP-IAS13), and China Agriculture Research System of MOF and MARA (CARS-38).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2025.1629162/full#supplementary-material

References

1. Zhang Z, Sui Z, Zhang J, Li Q, Zhang Y, Wang C, et al. Identification of signatures of selection for lambing number and pubertal initiation in two sheep populations. Animals. (2022) 12:2520. doi: 10.3390/ani12192520

PubMed Abstract | Crossref Full Text | Google Scholar

2. Raina VS, Kour A, Chakravarty AK, Vohra V. Marker-assisted selection vis-à-vis bull fertility: coming full circle-a review. Mol Biol Rep. (2020) 47:9123–33. doi: 10.1007/s11033-020-05919-0

PubMed Abstract | Crossref Full Text | Google Scholar

3. Yue L, Li W, Pan Y, Lan X, Zhang Q, Pan C. Polymorphism of GTF2A1 gene is associated with lambing number in sheep. Theriogenology. (2023) 208:194–200. doi: 10.1016/j.theriogenology.2023.06.004

PubMed Abstract | Crossref Full Text | Google Scholar

4. Wang Z, Wang X, Lan X, Zhu H, Qu L, Pan C. Polymorphism within the GATA binding protein 4 gene is significantly associated with goat lambing number. Anim Biotechnol. (2023) 34:4291–300. doi: 10.1080/10495398.2022.2147533

PubMed Abstract | Crossref Full Text | Google Scholar

5. An XP, Hou JX, Zhao HB, Li G, Bai L, Peng JY, et al. Polymorphism identification in goat GNRH1 and GDF9 genes and their association analysis with lambing number. Anim Genet. (2013) 44:234–8. doi: 10.1111/j.1365-2052.2012.02394.x

PubMed Abstract | Crossref Full Text | Google Scholar

6. El-Shorbagy HM, Abdel-Aal ES, Mohamed SA, El-Ghor AA. Association of PRLR, IGF1, and LEP genes polymorphism with milk production and lambing number in Egyptian Zaraibi goat. Trop Anim Health Prod. (2022) 54:321. doi: 10.1007/s11250-022-03316-2

PubMed Abstract | Crossref Full Text | Google Scholar

7. Váncza L, Horváth A, Seungyeon L, Rókusz A, Dezso K, Reszegi A, et al. SPOCK1 Overexpression suggests poor prognosis of ovarian cancer. Cancers. (2023) 15:2037. doi: 10.3390/cancers15072037

PubMed Abstract | Crossref Full Text | Google Scholar

8. Du Z, Lin Z, Wang Z, Liu D, Tian D, Xia L. SPOCK1 overexpression induced by platelet-derived growth factor-BB promotes hepatic stellate cell activation and liver fibrosis through the integrin α5β1/PI3K/Akt signaling pathway. Lab Invest. (2020) 100:1042–56. doi: 10.1038/s41374-020-0425-4

PubMed Abstract | Crossref Full Text | Google Scholar

9. Song X, Han P, Liu J, Wang Y, Li D, He J, et al. Up-regulation of SPOCK1 induces epithelial-mesenchymal transition and promotes migration and invasion in esophageal squamous cell carcinoma. J Mol Histol. (2015) 46:347–56. doi: 10.1007/s10735-015-9627-2

PubMed Abstract | Crossref Full Text | Google Scholar

10. Fortes MR, Reverter A, Zhang Y, Collis E, Nagaraj SH, Jonsson NN, et al. Association weight matrix for the genetic dissection of puberty in beef cattle. Proc Natl Acad Sci U S A. (2010) 107:13642–7. doi: 10.1073/pnas.1002044107

PubMed Abstract | Crossref Full Text | Google Scholar

11. Liu YZ, Guo YF, Wang L, Tan LJ, Liu XG, Pei YF, et al. Genome-wide association analyses identify SPOCK as a key novel gene underlying age at menarche. PLoS Genet. (2009) 5:e1000420. doi: 10.1371/journal.pgen.1000420

PubMed Abstract | Crossref Full Text | Google Scholar

12. Dvornyk V, Waqar ul H. Genetics of age at menarche: a systematic review. Hum Reprod Update. (2012) 18:198–210. doi: 10.1093/humupd/dmr050

PubMed Abstract | Crossref Full Text | Google Scholar

13. Lin YW, Wen YC, Hsiao CH, Lai FR, Yang SF, Yang YC, et al. Proteoglycan SPOCK1 as a poor prognostic marker promotes malignant progression of clear cell renal cell carcinoma via triggering the snail/slug-MMP-2 axis-mediated epithelial-to-mesenchymal transition. Cells. (2023) 12:352. doi: 10.3390/cells12030352

PubMed Abstract | Crossref Full Text | Google Scholar

14. Irwin JC, Kirk D, Gwatkin RB, Navre M, Cannon P, Giudice LC. Human endometrial matrix metalloproteinase-2, a putative menstrual proteinase hormonal regulation in cultured stromal cells and messenger RNA expression during the menstrual cycle. J Clin Invest. (1996) 97:438–47. doi: 10.1172/JCI118433

PubMed Abstract | Crossref Full Text | Google Scholar

15. Dolebo AT, Khayatzadeh N, Melesse A, Wragg D, Rekik M, Haile A, et al. Genome-wide scans identify known and novel regions associated with prolificacy and reproduction traits in a sub-Saharan African indigenous sheep (Ovis aries). Mamm Genome. (2019) 30:339–52. doi: 10.1007/s00335-019-09820-5

PubMed Abstract | Crossref Full Text | Google Scholar

16. Doyle KM, Russell DL, Sriraman V, Richards JS. Coordinate transcription of the ADAMTS-1 gene by luteinizing hormone and progesterone receptor. Mol Endocrinol. (2004) 18:2463–78. doi: 10.1210/me.2003-0380

PubMed Abstract | Crossref Full Text | Google Scholar

17. Robker RL, Russell DL, Yoshioka S, Sharma SC, Lydon JP, O'Malley BW, et al. Ovulation: a multi-gene, multi-step process. Steroids. (2000) 65:559–70. doi: 10.1016/S0039-128X(00)00114-8

Crossref Full Text | Google Scholar

18. Hernández-Delgado P, Felix-Portillo M, Martínez-Quintana JA. ADAMTS proteases: importance in animal reproduction. Genes. (2023) 14:1181. doi: 10.3390/genes14061181

PubMed Abstract | Crossref Full Text | Google Scholar

19. Freimann S, Ben-Ami I, Dantes A, Armon L, Ben Ya'cov-Klein A, Ron-El R, et al. Differential expression of genes coding for EGF-like factors and ADAMTS1 following gonadotropin stimulation in normal and transformed human granulosa cells. Biochem Biophys Res Commun. (2005) 333:935–43. doi: 10.1016/j.bbrc.2005.04.177

PubMed Abstract | Crossref Full Text | Google Scholar

20. Robker RL, Russell DL, Espey LL, Lydon JP, O'Malley BW, Richards JS. Progesterone-regulated genes in the ovulation process: ADAMTS-1 and cathepsin L proteases. Proc Natl Acad Sci U S A. (2000) 97:4689–94. doi: 10.1073/pnas.080073497

PubMed Abstract | Crossref Full Text | Google Scholar

21. Fang L, Li G, Liu G, Lee SW, Aaronson SA. p53 induction of heparin-binding EGF-like growth factor counteracts p53 growth suppression through activation of MAPK and PI3K/Akt signaling cascades. EMBO J. (2001) 20:1931–9. doi: 10.1093/emboj/20.8.1931

PubMed Abstract | Crossref Full Text | Google Scholar

22. Jessmon P, Leach RE, Armant DR. Diverse functions of HBEGF during pregnancy. Mol Reprod Dev. (2009) 76:1116–27. doi: 10.1002/mrd.21066

PubMed Abstract | Crossref Full Text | Google Scholar

23. Armant DR, Wang J, Liu Z. Intracellular signaling in the developing blastocyst as a consequence of the maternal-embryonic dialogue. Semin Reprod Med. (2000) 18:273–87. doi: 10.1055/s-2000-12565

PubMed Abstract | Crossref Full Text | Google Scholar

24. Lim H, Das SK, Dey SK. erbB genes in the mouse uterus: cell-specific signaling by epidermal growth factor (EGF) family of growth factors during implantation. Dev Biol. (1998) 204:97–110. doi: 10.1006/dbio.1998.9072

PubMed Abstract | Crossref Full Text | Google Scholar

25. Xie H, Wang H, Tranguch S, Iwamoto R, Mekada E, Demayo FJ, et al. Maternal heparin-binding-EGF deficiency limits pregnancy success in mice. Proc Natl Acad Sci U S A. (2007) 104:18315–20. doi: 10.1073/pnas.0707909104

PubMed Abstract | Crossref Full Text | Google Scholar

26. Stewart CL, Kaspar P, Brunet LJ, Bhatt H, Gadi I, Köntgen F, et al. Blastocyst implantation depends on maternal expression of leukaemia inhibitory factor. Nature. (1992) 359:76–9. doi: 10.1038/359076a0

PubMed Abstract | Crossref Full Text | Google Scholar

27. Song H, Lim H, Das SK, Paria BC, Dey SK. Dysregulation of EGF family of growth factors and COX-2 in the uterus during the preattachment and attachment reactions of the blastocyst with the luminal epithelium correlates with implantation failure in LIF-deficient mice. Mol Endocrinol. (2000) 14:1147–61. doi: 10.1210/me.14.8.1147

Crossref Full Text | Google Scholar

28. Yuan Z, Zhang J, Li W, Wang W, Li F, Yue X. Association of polymorphisms in candidate genes with the lambing number in two sheep breeds. Animals. (2019) 9:958. doi: 10.3390/ani9110958

PubMed Abstract | Crossref Full Text | Google Scholar

29. Ma B, Cui H, Wang X, Feng W, Zhang J, Chen N, et al. IFNT-induced IRF1 enhances bovine endometrial receptivity by transactivating LIFR. J Reprod Immunol. (2024) 163:104212. doi: 10.1016/j.jri.2024.104212

PubMed Abstract | Crossref Full Text | Google Scholar

30. Modric T, Kowalski AA, Green ML, Simmen RC, Simmen FA. Pregnancy-dependent expression of leukaemia inhibitory factor (LIF), LIF receptor-beta and interleukin-6 (IL-6) messenger ribonucleic acids in the porcine female reproductive tract. Placenta. (2000) 21:345–53. doi: 10.1053/plac.1999.0493

PubMed Abstract | Crossref Full Text | Google Scholar

31. Zhou M, Pan Z, Cao X, Guo X, He X, Sun Q, et al. Single nucleotide polymorphisms in the HIRA gene affect lambing number in small tail han sheep. Animals. (2018) 8:71. doi: 10.3390/ani8050071

PubMed Abstract | Crossref Full Text | Google Scholar

32. Li J, Gong Y, Wang X, He X, He X, Chu M, et al. Screening of litter-size-associated SNPs in NOX4, PDE11A and GHR genes of sheep. Animals. (2024) 14:767. doi: 10.3390/ani14050767

PubMed Abstract | Crossref Full Text | Google Scholar

33. Szklarczyk D, Kirsch R, Koutrouli M, Nastou K, Mehryary F, Hachilif R, et al. The STRING database in 2023: protein-protein association networks and functional enrichment analyses for any sequenced genome of interest. Nucleic Acids Res. (2023) 51:D638–46. doi: 10.1093/nar/gkac1000

PubMed Abstract | Crossref Full Text | Google Scholar

34. Crisponi L, Buers I, Rutsch F. CRLF1 and CLCF1 in development, health and disease. Int J Mol Sci. (2022) 23:992. doi: 10.3390/ijms23020992

PubMed Abstract | Crossref Full Text | Google Scholar

35. Wright AF, Carothers AD. CNTF in the embryo. Nat Genet. (1994) 7:460. doi: 10.1038/ng0894-460a

PubMed Abstract | Crossref Full Text | Google Scholar

36. Schwerd T, Krause F, Twigg SRF, Aschenbrenner D, Chen YH, Borgmeyer U, et al. A variant in IL6ST with a selective IL-11 signaling defect in human and mouse. Bone Res. (2020) 8:24. doi: 10.1038/s41413-020-0098-z

PubMed Abstract | Crossref Full Text | Google Scholar

37. Pan Z, Li S, Liu Q, Wang Z, Zhou Z, Di R, et al. Whole-genome sequences of 89 Chinese sheep suggest role of RXFP2 in the development of unique horn phenotype as response to semi-feralization. Gigascience. (2018) 7:giy019. doi: 10.1093/gigascience/giy019

PubMed Abstract | Crossref Full Text | Google Scholar

38. Ware CB, Horowitz MC, Renshaw BR, Hunt JS, Liggitt D, Koblar SA, et al. Targeted disruption of the low-affinity leukemia inhibitory factor receptor gene causes placental, skeletal, neural and metabolic defects and results in perinatal death. Development. (1995) 121:1283–99. doi: 10.1242/dev.121.5.1283

PubMed Abstract | Crossref Full Text | Google Scholar

39. Takahashi Y, Carpino N, Cross JC, Torres M, Parganas E, Ihle JN. SOCS3: an essential regulator of LIF receptor signaling in trophoblast giant cell differentiation. EMBO J. (2003) 22:372–84. doi: 10.1093/emboj/cdg057

PubMed Abstract | Crossref Full Text | Google Scholar

40. Ebrahimi B, Viswanadhapalli S, Pratap UP, Rahul G, Yang X, Pitta Venkata P, et al. Pharmacological inhibition of the LIF/LIFR autocrine loop reveals vulnerability of ovarian cancer cells to ferroptosis. NPJ Precis Oncol. (2024) 8:118. doi: 10.1038/s41698-024-00612-y

PubMed Abstract | Crossref Full Text | Google Scholar

41. Kritikou EA, Sharkey A, Abell K, Came PJ, Anderson E, Clarkson RW, et al. A dual, non-redundant, role for LIF as a regulator of development and STAT3-mediated cell death in mammary gland. Development. (2003) 130:3459–68. doi: 10.1242/dev.00578

PubMed Abstract | Crossref Full Text | Google Scholar

42. Abell K, Bilancio A, Clarkson RW, Tiffen PG, Altaparmakov AI, Burdon TG, et al. Stat3-induced apoptosis requires a molecular switch in PI(3)K subunit composition. Nat Cell Biol. (2005) 7:392–8. doi: 10.1038/ncb1242

PubMed Abstract | Crossref Full Text | Google Scholar

43. Namiki T, Terakawa J, Karakama H, Noguchi M, Murakami H, Hasegawa Y, et al. Uterine epithelial Gp130 orchestrates hormone response and epithelial remodeling for successful embryo attachment in mice. Sci Rep. (2023) 13:854. doi: 10.1038/s41598-023-27859-y

PubMed Abstract | Crossref Full Text | Google Scholar

44. Awonuga AO, Camp OG, Abu-Soud HM, Rappolee DA, Puscheck EE, Diamond MP. Determinants of embryo implantation: roles of the endometrium and embryo in implantation success. Reprod Sci. (2023) 30:2339–48. doi: 10.1007/s43032-023-01224-w

PubMed Abstract | Crossref Full Text | Google Scholar

45. Auernhammer CJ, Melmed S. Leukemia-inhibitory factor-neuroimmune modulator of endocrine function. Endocr Rev. (2000) 21:313–45. doi: 10.1210/edrv.21.3.0400

PubMed Abstract | Crossref Full Text | Google Scholar

46. Trouillas M, Saucourt C, Guillotin B, Gauthereau X, Taupin JL, Moreau JF, et al. The LIF cytokine: towards adulthood. Eur Cytokine Netw. (2009) 20:51–62. doi: 10.1684/ecn.2009.0148

PubMed Abstract | Crossref Full Text | Google Scholar

47. Kimber SJ. Leukaemia inhibitory factor in implantation and uterine biology. Reproduction. (2005) 130:131–45. doi: 10.1530/rep.1.00304

PubMed Abstract | Crossref Full Text | Google Scholar

48. Watanobe H, Habu S. Ciliary neurotrophic factor, a gp130 cytokine, regulates preovulatory surges of luteinizing hormone and prolactin in the rat. Neuroendocrinology. (2001) 74:281–7. doi: 10.1159/000054695

PubMed Abstract | Crossref Full Text | Google Scholar

49. Yoo I, Han J, Kim M, Jang H, Sa S, Choi SH, et al. Expression and regulation of interleukin 6 and its receptor at the maternal-conceptus interface during pregnancy in pigs. Theriogenology. (2017) 96:85–91. doi: 10.1016/j.theriogenology.2017.04.007

PubMed Abstract | Crossref Full Text | Google Scholar