A Novel SNP in the Promoter Region of IGF1 Associated With Yunshang Black Goat Kidding Number via Promoting Transcription Activity by SP1

IGF1, a member of the insulin-like growth factor (IGF) superfamily, is also known as the growth-promoting factor (somatomedin C). IGF1 is involved in vertebrate growth and development, immunity, cell metabolism, reproduction, and breeding. However, there are relatively few studies on the relationship between IGF1 and goat reproduction. In this study, a new transcription factor SP1 bound to the IGF1 g. 64943050T>C promoted granulosa cell (GC) proliferation. A mutation g.64943050T>C located in the promoter region of IGF1 was identified. Association analysis revealed that the kidding number in the first and second litters and the average number of first three litters of the CC genotype (2.206 ± 0.044, 2.254 ± 0.056, and 2.251 ± 0.031) were significantly higher than those in the TC genotype (1.832 ± 0.049, 1.982 ± 0.06, and 1.921 ± 0.034) and TT genotype (1.860 ± 0.090, 1.968 ± 0.117, and 1.924 ± 0.062) (p < 0.05). The kidding number in the third litter of the CC genotype (2.355 ± 0.057) was significantly higher than that in the TT genotype (2.000 ± 0.107) (p < 0.05). Then, the function of this mutation was validated by the dual-luciferase reporter assay and EMSA. The results showed that the luciferase activity of IGF1-mutant-C was significantly higher than that of IGF1-Wild-T (p < 0.05). The EMSA also showed that the binding ability of IGF1-mutant-C was higher than that of IGF1-Wild-T (p < 0.05). Subsequently, the transcription factor SP1 was predicted to bind to the mutation of IGF1 (g.64943050T>C). Overexpression of SP1 promotes the expression of IGF1 in the primary granulosa cells (GCs). The results of the CCK-8 assay and the expression of GC proliferation factors (CDK4, cyclin D1, and cyclin D2) demonstrated that SP1 promoted GC proliferation by regulating IGF1 expression. Our results suggested that the IGF1 g.64943050T>C was significantly associated with the kidding number of Yunshang black goats, and SP1 as a transcription factor of IGF1 binding to the mutation T>C regulated the expression of IGF1. Furthermore, SP1 promoted goat GC proliferation by regulating the expression of IGF1, which provides a new insight for the goat fertility trait.

Introduction

The carcasses of goats are very lean compared to those of sheep at the same age (Stanišić et al., 2012). Nutritionally, goat meat is an important source of high-quality protein, healthy fats, low calorie intramuscular fats, saturated fats, and sodium. In addition, goat meat has high levels of iron, potassium, and essential amino acids (Mcmillin and Brock, 2005; Givens et al., 2006; Lee et al., 2008; Horcada et al., 2012). However, goats are more difficult to rear than sheep, and their growth rate is generally slower and smaller than that of sheep, making goat meat scarce and significantly more expensive than sheep (Brady et al., 2020). Increasing the kidding number of goats is the easiest and most direct way to increase goat meat production. Compared to traditional breeding methods, molecular genetic breeding has several advantages, including high accuracy, operability, cost saving, and the ability to quickly and efficiently improve goat prolificacy traits.

Insulin-like growth factor 1 (IGF1), known as growth-promoting factor (somatomedin C), has insulin-like functions and significant diversity in biological effects because of its molecular structure being extremely similar to insulin (Ma et al., 2015). IGF1 is conserved across species, and IGF1 can be selectively spliced into multiple transcripts encoding different precursor IGF1 proteins. In goats, IGF1 is located on chromosome 5 and contains six exons (Górecki et al., 2007; Xiao et al., 2009). IGF1 is an important component of the insulin-like growth factor (IGF) superfamily and is commonly expressed in mammalian gonadal axis tissues such as the hypothalamus, anterior pituitary, ovary, oviduct, and uterus (Bach and Bondy, 1992; Mikawa et al., 1995; Adam et al., 2000; Hunter et al., 2004; Filby and Tyler, 2007). IGF1 is mainly secreted by the liver, and GH (growth hormone) secreted by the pituitary acts on the GHR (growth hormone receptor) of the liver to regulate the release of IGF1. IGF1 acts as a systemic growth factor and also acts on the hypothalamus as a negative feedback regulator of growth hormone secretion. Meanwhile, both growth hormone and insulin-like growth factor are affected by binding proteins (GHBP and IGFBP, respectively) (Hunter et al., 2004). IGF1 acts on IGF1 receptors and is generally involved in vertebrate growth and development, organismal immunity, cellular metabolism, and reproduction (Cannata et al., 2010; Ni et al., 2013; Hu and Liu, 2014; Kindler et al., 2016). The receptors of IGF1 are found in almost every organ of the mammalian reproductive system, regulating reproductive system processes from the center to the organs. IGF1, as a hormonal regulator, acts directly on GnRH neurons, stimulating their transcriptional translation. It also acts synergistically with GnRH to promote follicular growth and steroid hormone synthesis, thereby controlling the initiation of the first estrus and reproductive activity in adult mammals (Lucy, 2000; Miller and Gore, 2001). Previous studies have also shown that polymorphisms in IGF1 are significantly associated with lambing numbers in sheep and goats (He et al., 2012; Thomas et al., 2016; Naicy et al., 2017; Tao et al., 2021). This suggests that IGF1 plays an important role in regulating the reproductive performance of females.

In this study, we detected an SNP in the promoter region of the IGF1, determined the relationship between mutations in the promoter region of the IGF1 and kidding number of goats, predicted the transcription factor SP1 bound to this mutation locus, and explored the potential molecular mechanisms between the expression of transcription factors SP1 and IGF1 and reproductive performance in goats. These results are important for elucidating the molecular characteristics of IGF1, understanding its function in regulating kidding number of goats, and comprehensively exploring the molecular mechanism of the IGF family in regulating the fertility of goats.

Materials and Methods

Animals and Sample Collection

Blood was collected from 400 Yunshang black goats (2–5 years old). All the goats were raised in Honghe Hani and Yi Autonomous Prefectures in Yunnan Province, China in the same environment. Genomic DNA was extracted using the phenol chloroform method. DNA sample concentration was detected by Nano Drop 2000 (Thermo Fisher, USA), and DNA quality was detected by 1.2% agarose gel electrophoresis. Ten high- and low- yielding Yunshang black goats were selected, and after euthanasia, ovarian tissues were extracted, put into liquid nitrogen immediately, and then stored at −80°C. RNA was extracted from the tissue samples using the Animal Tissue RNA Extraction Kit (Tiangen, China), and RNA concentration was detected using Nano Drop 2000, and RNA quality was detected by 1.2% agarose gel electrophoresis.

Reverse Transcription, Real-Time Quantitative PCR, and Genotyping

The cDNA was synthesized using the TaKaRa Reverse Transcription Kit (TaKaRa, Japan). Reverse transcription system: 1.0 μl reverse transcription mixed primer, 1 μl Oligo dT primer, 1 μl random primer, 4 μl fluorescence quantification buffer, 1 μl total RNA, 12 μl double-distilled water, and operated on ice throughout. Reaction conditions: 37°C for 15 min; 85°C for 5 s. The cDNA product after the completion of reverse transcription was diluted, and PCR was performed with the housekeeping gene RPL19, and the quality was qualified and stored at −20°C for the detection of gene mRNA expression.

The total volume of the PCR system was 20 μl: 10 μl of Ex Taq primer II, 0.8 μl of the upstream primer, 0.8 μl of the downstream primer, 2 μl of the cDNA template, and 6.4 ul of double-distilled water. RPL19 was used as the control gene, and the primer sequences are shown in Table 1.

TABLE 1

TABLE 1. Primers used in this study.

The IGF1 was typed using the KASP method using a 384-well plate with a total volume of 10.14 µl. The reaction system included 2.5 µl DNA template, 2.5 µl 2×Master Mix, 0.07 µl primer mix, 2.57 µl working solution, and 2.5 µl RNase-free water blown and mixed and then performed instantaneously at 3,000 rpm. After PCR, the wells were removed, and the data were scanned.

Bioinformatics Analysis

The binding sites of transcription factors at the location of the mutant site and the transcription factors bound were predicted using Jaspar (http://jaspar.genereg.net/).

Plasmid Construction and Extraction

We constructed plasmids for the IGF1 wild type (IGF1-Wild-T) and mutant type (IGF1-Mutant-C), and we inserted the +/-100bp base before and after the mutation site into the luciferase reporter vector PGL3-basic plasmid to test whether the fragment possesses initiation activity. The overexpression plasmid of transcription factor SP1 was constructed using the pIRES2-EGFP plasmid, and we inserted the CDS region of transcription factor SP1 into the pIRES2-EGFP plasmid. For SP1 interference, we used Sigma’s pLKO.1-puro, which is a shRNA vector, and we designed three shRNAs and experimentally verified the three shRNA sequences shown in Supplementary Table S1. The plasmids were extracted using the QIAGEN Plasmid Midi Kit (QIAGEN, Germany).

GC Isolation and Culture

To isolate primary granulosa cells, we extracted fresh goat ovarian follicles in a sterile Hank’s balanced salt solution. The ovarian cells were rinsed repeatedly with PBS, the cells were filtered through a 200-mesh cell sieve and collected by centrifugation, and then the primary cells were cultured in 6-cm dishes for 48 h, during which the cell culture medium was changed once to ensure cell survival and the cells were transferred to 10-cm dishes for subsequent experiments. The HEK293T cell lines were obtained from previous cultures in the laboratory, and the cell culture conditions were 37°C, 5% CO₂, and 95% air saturation humidity. Both goat pellet cells and HEK293T cell lines were cultured using DMEM (Thermo Fisher, USA) with 10% fetal bovine serum and 100 U/ml penicillin/streptomycin. HEK293T cell lines were used to assay luciferase activity.

Transfection and Dual-Luciferase Reporter Assay

Transfection was performed using Lip 2000 (Thermo Fisher, USA), and a dual-luciferase activity assay was performed on 24-well plates, and HEK293T cells were transfected when their density grew to about 70%. Opti-DMEM (Thermo Fisher, USA) was used for transfection, and the cells were changed 6 h after transfection and were collected and processed after 48 h of incubation to determine luciferase activity. Goat primary pellet cells were used for the gene expression assay. The cells were transfected in 6-well plates until the pellet cell density reached about 70%, and the cells were changed after 6 h. The cells were collected after 48 h of incubation for RNA extraction.

Ethynyldeoxyuridine Incorporation Assay

Primary GCs were cultured in 96-well plates. Cell proliferation was detected after transfection with the CCK-8 Cell Proliferation and Cytotoxicity Assay Kit (Solarbio, China). After transfection treatment of the cells, the solution was changed after 6 h, and 10 µl of the CCK-8 reagent was added to each well. After incubation in the cell incubator for 2 h, the optical density values were detected at 450 nm with a microplate reader (Thermo Fisher, USA) at 0, 6, 12, 24, and 48 h. EdU was used to determine the proliferation status of granulocytes using the EdU-488 Cell Proliferation Assay Kit (Beyotime, China) according to the manufacturer’s guidelines. Using a 6-well plate, the final concentration of EdU was added to each well at 10 µM, and incubation was continued for 6 h in a cell culture incubator. The, the cells were washed three times with PBS (Thermo Fisher, USA), then fixed with a fixative for 15 min at room temperature, the fixative was removed, and the cells were rinsed three times using washing solution. 1 ml of 0.3% Triton X-100 PBS was added and incubated for 15 min at room temperature to increase cell membrane permeability of the cell membrane. In addition, 500 µl of the prepared Click reaction solution was added to each well, and the cells were incubated at room temperature for 30 min in the dark to observe and quantify the number of cells stained with EdU. Three fields were randomly selected for statistical analysis.

Electrophoretic Mobility Shift Assay

The nuclear extracts were prepared from goat ovaries using the Nuclear Protein Extraction Kit (Beyotime, China) according to the manufacturer’s instructions. The EMSA test uses the EMSA kit (Thermo Fisher, USA), and biotin-labeled probes were synthesized by the company (Saville, China), and the biotin-labeled probe sequences are shown in Supplementary Table S2. All experimental steps were performed according to the manufacturer’s instructions.

Statistical Analysis

PopGenev 1.31 software (http://www.ualberta.ca/∼fyeh/fyeh) was used. Allele frequencies, heterozygosity, and polymorphic information content were calculated. The method of association analysis between different genotypes and the number of lambing were analyzed using a general linear model. All results are shown as mean ± SEM. Statistical analyses were performed using IBM SPSS Statistics 25 statistical software (SPSS Inc., USA) using independent sample t-test.

Results

Association Between IGF1 g.64943050T>C Mutation and Litter Size and Population Genetic Analysis

In our study, the IGF1 g.64943050T>C mutation locus was present in three genotypes (TT, TC, and CC) in Yunshang black goat, with a detection rate of >90% in 400 samples, in which the dominant genotype was CC and the dominant gene was C. Population genetic statistics showed that the mutation was moderately polymorphic in Yunshang black goat (0.25<PIC < 0.5). The chi-square fitness test showed that the distribution of this SNP was in Hardy–Weinberg equilibrium (p > 0.05) (Table 2).

TABLE 2

TABLE 2. Population polymorphism analysis of locus in goat.

The association analysis between g.64943050T>C and the kidding number of Yunshang black goats was performed, and the results are shown in Table 3. The kidding number of goats with the CC genotype (2.206 ± 0.044, 2.254 ± 0.056) individuals was significantly higher than that of the TC genotype (1.832 ± 0.049, 1.982 ± 0.061) and TT genotype (1.860 ± 0.090, 1.968 ± 0.117) (p < 0.05). In the third litter, the number of kids born of the CC genotype (2.355 ± 0.057) was significantly higher than that of the TT genotype (2.000 ± 0.107) (p < 0.05). The mean value of kidding number in all three litters was also significantly higher in the CC genotype (2.251 ± 0.031) than in the TC (1.921 ± 0.034) and TT (1.924 ± 0.062) genotypes. These results suggest that the IGF1 g.64943050T>C mutation was involved in the regulation of the kidding number in Yunshang black goats.

TABLE 3

TABLE 3. Least squares means and standard errors of litter size in black goats with different genotypes.

In addition, we quantified the expression of IGF1 in the ovaries of Yunshang black goats from pure individuals with different genotypes, and the results showed that the IGF1 expression was significantly higher in the individuals with the C mutation than in those with wild type T (Figure 1) (p < 0.05).

FIGURE 1

FIGURE 1. Detection of the expression of different genotypes of IGF1 in goat ovaries. Bars represent the mean ± SEM of at least three repeats. *p < 0.05; **p < 0.01.

Effect of g.64943050T>C Mutation on the Promoter Activity of IGF1

To investigate the effect of g.64943050T>C mutation on the promoter activity of IGF1, the IGF1-Wild-T and IGF1-Mutant-C were constructed and transfected into HEK293T cell lines. The results showed that the luciferase activity of the IGF1-Mutant-C group was significantly higher than that in the IGF1-Wild-T group (Figure 2). The EMSA analysis showed that the bio-probe IGF1-Mutant-C significantly bound more proteins and formed more DNA–protein complexes compared to the bio-probe IGF1-Wild-T (Figure 3).

FIGURE 2

FIGURE 2. Effect of mutations on the promoter activity of IGF1. IGF1-Mutant-C represents the mutant type, IGF1-Wild-T represents the wild type, and PGL3-basic represents control. They were transfected into 293T cells. The results are expressed as the mean ± SEM (n = 3) in arbitrary units based on firefly luciferase activity normalized against Renilla luciferase activity. A t-test was conducted using SPSS 25.0 to detect the differences. Bars represent the mean ± SEM of at least three repeats. *p < 0.05; **p < 0.01.

FIGURE 3

FIGURE 3. Binding of different biological probes to ovarian nucleoprotein. (A) Binding of the probe to the nucleoprotein. These different probes were incubated with ovary nuclear extracts, and DNA–protein complexes were visualized by autoradiography. DNA–protein complex: DNA and transcription factor SP1 binding complex. (B) Sample gray value.

Mutation of g.64943050T>C Adds a New Transcription Factor Binding Site

Jaspar online software was used to predict the transcription factor binding site of g.64943050T>C mutation of the IGF1 promoter region. The result showed that transcription factor SP1 was one of the transcription factors that obtained the highest score (Supplementary Figure S2).

Transcription Factor SP1 Affects the Expression Abundance of IGF1

To investigate the role of transcription factor SP1 in the transcription process of IGF1, the overexpression and interference vectors for transcription factor SP1 were constructed by pIRES2-SP1 and pLK0.1 puro-SP1, respectively. One of the best interference efficiencies of the interference vector was chosen, and the transfection efficiency is shown in Figure 4A. The RT-qPCR results showed that the overexpression of SP1 in goat GCs transfected with pIRES2-SP1 was significantly higher than that in the control group (p < 0.01) (Figure 4C). Correspondingly, the expression of IGF1 was also significantly higher in goat granulosa cells transfected with pIRES2-SP1. Similarly, the expression levels of SP1 and IGF1 were significantly lower in goat GCs transfected with pLK0.1 puro-SP1 than those in the control group (p < 0.05) (Figure 4D). These results suggest that the transcription factor SP1 is involved in regulating the expression of IGF1 in goat GCs.

FIGURE 4

FIGURE 4. Transcription factor SP1 affects the expression of IGF1. (A) Transfection efficiency of the vector in granular cells. (B) Three SP1 interference efficiency graphs were represented using RT-qPCR to detect the expression level of SP1. (C,D) SP1 overexpression and interference vectors were transfected into goat granulosa cells, and the expression levels of SP1 and IGF1 were detected by RT-qPCR. Bars represent the mean ± SEM of at least three repeats. *p < 0.05, **p < 0.01.

Transcription Factor SP1 Affects the Proliferation of GCs

To further investigate the effect of transcription factor SP1 on goat reproduction, we transfected pIRES2-SP1 and pLK0.1 puro-SP1 into goat ovarian GCs to detect the effect of SP1 on GC proliferation. The results of both CCK-8 and EdU showed that the proliferation rate of GCs transfected with pIRES2-SP1 was higher in GCs than in the control group (Figures 5A,E). Similarly, the proliferation rate of GCs transfected with pLK0.1 puro-SP1 was lower than that of the control group (Figures 5B,E). To further verify the effect of transcription factor SP1 on the proliferation of goat GCs, we detected the expression of the cell proliferation factors using RT-qPCR. The results showed that the expression of CDK4, Cyclin-D1, and Cyclin-D2 was significantly higher in GCs transfected with pIRES2-SP1 than that in the control group (Figure 5C), and similarly, the cell proliferation factors were all significantly lower in GCs transfected with pLK0.1 puro-SP1 than that in the control group (Figure 5D). All these results indicated that the transcription factor SP1 had an effect on goat GC proliferation.

FIGURE 5

FIGURE 5. Transcription factor SP1 regulates the proliferation of granulosa cells. (A) After SP1 was overexpressed, the cell counting kit-8 (CCK-8) was used to determine the proliferation curve of granulosa cells. (B) After interfering with SP1, Cell Counting Kit-8 (CCK-8) was used to determine the proliferation curve of granular cells. (C) Expression abundance of cell proliferation–related genes (CDK4, cyclinD1, and cyclinD2) was detected after overexpression of SP1 by RT-qPCR. (D) Expression abundance of cell proliferation–related genes (CDK4, cyclinD1, and cyclinD2) was detected after interference with SP1 by RT-qPCR. (E) Granulosa cells stained positive for 5-ethynyl-2-deoxyuridine (EdU) were detected with an EdU kit after SP1 overexpression and inhibition, EdU (red), Hoechst (blue). (F,G) Fold change of granulosa cell proliferation rate after SP1 overexpression and inhibition. Bars represent the mean ± SEM of at least three repeats. *p < 0.05; **p < 0.01.

Discussion

IGF1, as a member of the IGF superfamily, is an important growth factor that plays an important role in a variety of physiological processes, including reproduction, fetal development, and growth in mammals (Shen et al., 2005; Velazquez et al., 2008). Studies have shown that IGF1 is closely associated with gestation time (Echternkamp et al., 2004), twin ovulation (Velazquez et al., 2005), preimplantation embryonic development (Patton et al., 2007), and conception rate (Yilmaz et al., 2005), while it has also been reported that polymorphisms in IGF1 are associated with reproductive traits in multiple species, such as sheep (Li et al., 2006), cattle (Andrade et al., 2008; Kim et al., 2009), chicken (Bian et al., 2008), and humans (Vella et al., 2008). In addition, IGF1 plays an important role in the regulation of many hormones that are essential for reproduction. Mammalian reproduction is mainly controlled by gonadotropins (LH and FSH), and IGF1 promotes gonadotropin function by enhancing the function of gonadotropin receptors (Scaramuzzi et al., 2006), while IGF1 also stimulates ovarian function and steroidogenesis through the secretion of gonadotropin-releasing hormone (Luginbuehl et al., 2013). IGF1 is able to act directly through its receptors in the ovary. Activation of IGF1R stimulates ovarian estradiol and luteinizing hormone production (Demeestere et al., 2004).

In this study, a mutation located in the promoter region of IGF1 was identified in goats. Previous studies have shown that polymorphisms in the 5’ regulatory region and the CDS region of the IGF1 were associated with lambing number in four sheep breeds, namely, the Small-Tail Han sheep, the Hu sheep, Texel, and Dorset ewes, which were significantly associated with lambing number (He et al., 2012). In this study, it was found that the IGF1g.64943050T>C locus was significantly associated with the litter size produced in the first three litters of goats. In the first two litters, the litter size produced by CC individuals was 0.374 and 0.224 higher than that of TC and TT individuals (Table 3). In the third litter, the litter size of the CC type individuals was 0.335 higher than that of the TC type individuals (Table 3). Similarly, the average number of kids born in the first three litters was significantly higher in the CC genotype than that in the TC and TT genotypes. This suggests that IGF1 could be a candidate molecular marker for high-fertility goat breeding. Furthermore, the expression of IGF1 in the ovaries of homozygote individuals with different genotypes was detected, and the results showed that IGF1 was differentially expressed between high- and low-yielding goat groups, which indicated that IGF1 was an important candidate gene related to the goat prolificacy trait (Figure 1).

To further investigate the effect of mutations of the IGF1 5′ regulatory region on reproductive performance of goats, we predicted the transcription factors bound near the mutation site and found that the mutation of T>C created a new binding site for transcription factor SP1. Gene expression is a complex process, and transcription initiation is crucial for gene expression (de Klerk and t Hoen, 2015; Haberle and Stark, 2018). The transcription factors play an important role in eukaryotic transcription initiation. Transcription factor SP1 (specificity protein 1), a commonly expressed transcription factor, is widely expressed in mammals, and the family members include SP2, SP3, and SP4 (Suske, 1999). Transcription factor SP1 was initially found in HeLa cell lines and is bound to multiple GGGCGG sequences (GC box) in the 21 bp repeat element of SV40 and activated in vitro transcription of the SV40 early promoter (Dynan and Tjian, 1983a; Dynan and Tjian, 1983b; Gidoni et al., 1984; Gidoni et al., 1985). Several studies have demonstrated that the transcription factor SP1 is important for cell growth, differentiation, apoptosis, and oncogenesis (Kadonaga et al., 1987; Vizcaíno et al., 2015). We speculated that the IGF1 T>C mutation leading to promoter activity change may be associated with the newly generated binding site of transcription factor SP1. Therefore, overexpression and interference treatment of transcription factor SP1 were constructed into vectors and then transfected into goat GCs. The results showed that the expression of IGF1 was significantly increased in the goat GCs with the SP1 overexpressed group (Figure 4C), while the expression level of IGF1 was significantly reduced in the goat GCs with the siR-SP1 group (Figure 4D). This indicated that the expression of transcription factor SP1 significantly influenced the expression of IGF1 in goat GCs.

The results of EdU and CCK-8 transfected with overexpressed SP1 and siR-SP1 suggested that the transcription factor SP1 may affect the proliferation of goat GCs by influencing the expression of IGF1 and, thus, the proliferation of GCs. It has also been previously shown that IGF1 increased the rate of reproductive hormone production and had an effect on cell growth with specific reproductive structures (McMurtry, 1998). In mammals and birds, IGF1 has been shown to promote mitosis in ovarian GCs, luteal GCs, and follicular membrane cells (McMurtry, 1998). In antral follicles, IGFs stimulate granulosa cell proliferation and steroidogenesis in most mammals (Silva et al., 2009). These studies further corroborate our results. CDK4, Cyclin-D1, and Cyclin-D2 are the common cell proliferation factors (Jiang et al., 2021). In this study, the results showed that the expression of CDK4, Cyclin-D1, and Cyclin-D2 was significantly higher in GCs transfected with pIRES2-SP1 than that in the control group (Figure 5C), and the opposite trend was observed in GCs transfected with pLK0.1 puro-SP1. Overall, the transcription factor SP1 influenced the proliferation of goat GCs by regulating IGF1.

In summary, we identified a T>C mutation in the promoter region of IGF1, which resulted in an increase in IGF1 promoter activity. The mutation also created a new binding site for transcription factor SP1, and the increased activity of the IGF1 promoter may be associated with SP1. Transcription factor SP1 affected goat GC proliferation by influencing the expression of IGF1, and IGF1 was associated with ovulation and the kidding number of goats (Figure 6). This study provides a new insight for goat breeding.

FIGURE 6

FIGURE 6. Flow chart of the IGF1 g.64943050T>C mutation regulating goat reproduction through the transcription factor SP1.

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

The animal study was reviewed and approved by All experimental procedures mentioned in the present study were approved by the Science Research Department (in charge of animal welfare issues) of the Institute of Animal Sciences, Chinese Academy of Agricultural Sciences (IAS-CAAS) (Beijing, China). Ethical approval was provided by the Animal Ethics Committee of IASCAAS (No. IAS2021-25).

Author Contributions

Conceptualization and methodology: KL, YL, and LT; Experimental verification and data analysis: KL and YL; Investigation: YJ, QH, and RL; writing—original draft preparation: KL and YL; Writing—review and editing: KL, YL, and MC; Project administration: QH, XH, and MC; Funding acquisition, XH and MC.

Funding

This work was financially supported by the National Natural Science Foundation of China Youth Program (32102509), the Agricultural Science and Technology Innovation Program of China (CAAS-ZDRW202106 and ASTIP-IAS13), the China Agriculture Research System of MOF and MARA (CARS-38), the Major Science and Technology Project of Yunnan Province (202102AE090039), and the Basic Research Foundation Key Project of Yunnan Province (202001AS070002).

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.

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/fcell.2022.873095/full#supplementary-material

Supplementary Figure S1 | KASP typing results of IGF1. Blue represents the TT genotype, green represents the TC genotype, and red represents the CC genotype.

Supplementary Figure S2 | Transcription factor prediction results.

Supplementary Table S1 | Transcription factor SP1 interferes with the shRNA sequence.

Supplementary Table S2 | EMSA probe sequence.

References

Adam, C., Gadd, T., Findlay, P., and Wathes, D. (2000). IGF-I Stimulation of Luteinizing Hormone Secretion, IGF-Binding Proteins (IGFBPs) and Expression of mRNAs for IGFs, IGF Receptors and IGFBPs in the Ovine Pituitary Gland. J. Endocrinol. 166, 247–254. doi:10.1677/joe.0.1660247

PubMed Abstract | CrossRef Full Text | Google Scholar

Andrade, P. C., Grossi, D. A., Paz, C. C. P., Alencar, M. M., Regitano, L. C. A., and Munari, D. P. (2008). Association of Aninsulin-like Growth Factor 1gene Microsatellite with Phenotypic Variation and Estimated Breeding Values of Growth Traits in Canchim Cattle. Anim. Genet. 39, 480–485. doi:10.1111/j.1365-2052.2008.01755.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bach, M. A., and Bondy, C. A. (1992). Anatomy of the Pituitary Insulin-like Growth Factor System. Endocrinology 131, 2588–2594. doi:10.1210/endo.131.6.1280202

PubMed Abstract | CrossRef Full Text | Google Scholar

Bian, L. H., Wang, S. Z., Wang, Q. G., Zhang, S., Wang, Y. X., and Li, H. (2008). Variation at the Insulin-like Growth Factor 1 Gene and its Association with Body Weight Traits in the Chicken. J. Anim. Breed. Genet. 125, 265–270. doi:10.1111/j.1439-0388.2008.00739.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Brady, M., Refshauge, G., Robertson, S., Atkinson, T., Allworth, B., and Hernandez‐Jover, M. (2020). An Exploratory Study to Investigate Animal Health and Reproductive Wastage Among Australian Meat Goat Producers. Aust. Vet. J. 98, 602–609. doi:10.1111/avj.13033

PubMed Abstract | CrossRef Full Text | Google Scholar

Cannata, D., Lann, D., Wu, Y., Elis, S., Sun, H., Yakar, S., et al. (2010). Elevated Circulating IGF-I Promotes Mammary Gland Development and Proliferation. Endocrinology 151, 5751–5761. doi:10.1210/en.2010-0792

PubMed Abstract | CrossRef Full Text | Google Scholar

de Klerk, E., and ‘t Hoen, P. A. C. (2015). Alternative mRNA Transcription, Processing, and Translation: Insights from RNA Sequencing. Trends Genet. 31, 128–139. doi:10.1016/j.tig.2015.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Demeestere, I., Gervy, C., Centner, J., Devreker, F., Englert, Y., and Delbaere, A. (2004). Effect of Insulin-like Growth Factor-I during Preantral Follicular Culture on Steroidogenesis, In Vitro Oocyte Maturation, and Embryo Development in Mice1. Biol. Reprod. 70, 1664–1669. doi:10.1095/biolreprod.103.023317

PubMed Abstract | CrossRef Full Text | Google Scholar

Dynan, W. S., and Tjian, R. (1983a). Isolation of Transcription Factors that Discriminate between Different Promoters Recognized by RNA Polymerase II. Cell 32, 669–680. doi:10.1016/0092-8674(83)90053-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Dynan, W. S., and Tjian, R. (1983b). The Promoter-specific Transcription Factor Sp1 Binds to Upstream Sequences in the SV40 Early Promoter. Cell 35, 79–87. doi:10.1016/0092-8674(83)90210-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Echternkamp, S. E., Roberts, A. J., Lunstra, D. D., Wise, T., and Spicer, L. J. (2004). Ovarian Follicular Development in Cattle Selected for Twin Ovulations and Births12. J. Anim. Sci. 82, 459–471. doi:10.2527/2004.822459x

PubMed Abstract | CrossRef Full Text | Google Scholar

Filby, A. L., and Tyler, C. R. (2007). Cloning and Characterization of cDNAs for Hormones And/or Receptors of Growth Hormone, Insulin-like Growth Factor-I, Thyroid Hormone, and Corticosteroid and the Gender-, Tissue-, and Developmental-specific Expression of Their mRNA Transcripts in Fathead Minnow (Pimephales promelas). Gen. Comp. Endocrinol. 150, 151–163. doi:10.1016/j.ygcen.2006.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Gidoni, D., Dynan, W. S., and Tjian, R. (1984). Multiple Specific Contacts between a Mammalian Transcription Factor and its Cognate Promoters. Nature 312, 409–413. doi:10.1038/312409a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Gidoni, D., Kadonaga, J. T., Barrera-Saldaña, H., Takahashi, K., Chambon, P., and Tjian, R. (1985). Bidirectional SV40 Transcription Mediated by Tandem Sp1 Binding Interactions. Science 230, 511–517. doi:10.1126/science.2996137

PubMed Abstract | CrossRef Full Text | Google Scholar

Givens, D. I., Kliem, K. E., and Gibbs, R. A. (2006). The Role of Meat as a Source of N−3 Polyunsaturated Fatty Acids in the Human Diet. Meat Sci. 74, 209–218. doi:10.1016/j.meatsci.2006.04.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Górecki, D. C., Beręsewicz, M., and Zabłocka, B. (2007). Neuroprotective Effects of Short Peptides Derived from the Insulin-like Growth Factor 1. Neurochem. Int. 51, 451–458. doi:10.1016/j.neuint.2007.04.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Haberle, V., and Stark, A. (2018). Eukaryotic Core Promoters and the Functional Basis of Transcription Initiation. Nat. Rev. Mol. Cel Biol 19, 621–637. doi:10.1038/s41580-018-0028-8

PubMed Abstract | CrossRef Full Text | Google Scholar

He, J. N., Zhang, B. Y., Chu, M. X., Wang, P. Q., Feng, T., Cao, G. L., et al. (2012). Polymorphism of Insulin-like Growth Factor 1 Gene and its Association with Litter Size in Small Tail Han Sheep. Mol. Biol. Rep. 39, 9801–9807. doi:10.1007/s11033-012-1846-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Horcada, A., Ripoll, G., Alcalde, M. J., Sañudo, C., Teixeira, A., and Panea, B. (2012). Fatty Acid Profile of Three Adipose Depots in Seven Spanish Breeds of Suckling Kids. Meat Sci. 92, 89–96. doi:10.1016/j.meatsci.2012.04.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, F., and Liu, F. (2014). Targeting Tissue-specific Metabolic Signaling Pathways in Aging: the Promise and Limitations. Protein Cell 5, 21–35. doi:10.1007/s13238-013-0002-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Hunter, M. G., Robinson, R. S., Mann, G. E., and Webb, R. (2004). Endocrine and Paracrine Control of Follicular Development and Ovulation Rate in Farm Species. Anim. Reprod. Sci. 82-83, 461–477. doi:10.1016/j.anireprosci.2004.05.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, M., Xu, S., Bai, M., and Zhang, A. (2021). The Emerging Role of MEIS1 in Cell Proliferation and Differentiation. Am. J. Physiology-Cell Physiol. 320, C264–c269. doi:10.1152/ajpcell.00422.2020

CrossRef Full Text | Google Scholar

Kadonaga, J. T., Carner, K. R., Masiarz, F. R., and Tjian, R. (1987). Isolation of cDNA Encoding Transcription Factor Sp1 and Functional Analysis of the DNA Binding Domain. Cell 51, 1079–1090. doi:10.1016/0092-8674(87)90594-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, E.-S., Shi, X., Cobanoglu, O., Weigel, K., Berger, P. J., and Kirkpatrick, B. W. (2009). Refined Mapping of Twinning-Rate Quantitative Trait Loci on Bovine Chromosome 5 and Analysis of Insulin-like Growth Factor-1 as a Positional Candidate Gene1. J. Anim. Sci. 87, 835–843. doi:10.2527/jas.2008-1252

PubMed Abstract | CrossRef Full Text | Google Scholar

Kindler, J. M., Pollock, N. K., Laing, E. M., Jenkins, N. T., Oshri, A., Isales, C., et al. (2016). Insulin Resistance Negatively Influences the Muscle-dependent IGF-1-Bone Mass Relationship in Premenarcheal Girls. J. Clin. Endocrinol. Metab. 101, 199–205. doi:10.1210/jc.2015-3451

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, J. H., Kouakou, B., and Kannan, G. (2008). Chemical Composition and Quality Characteristics of Chevon from Goats Fed Three Different post-weaning Diets. Small Ruminant Res. 75, 177–184. doi:10.1016/j.smallrumres.2007.10.003

CrossRef Full Text | Google Scholar

Li, M. H., Adamowicz, T., Switonski, M., Ammosov, I., Ivanova, Z., Kiselyova, T., et al. (2006). Analysis of Population Differentiation in North Eurasian Cattle (Bos taurus) Using Single Nucleotide Polymorphisms in Three Genes Associated with Production Traits. Anim. Genet. 37, 390–392. doi:10.1111/j.1365-2052.2006.01479.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lucy, M. C. (2000). Regulation of Ovarian Follicular Growth by Somatotropin and Insulin-like Growth Factors in Cattle. J. Dairy Sci. 83, 1635–1647. doi:10.3168/jds.S0022-0302(00)75032-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Luginbuehl, V., Zoidis, E., Meinel, L., von Rechenberg, B., Gander, B., and Merkle, H. P. (2013). Impact of IGF-I Release Kinetics on Bone Healing: a Preliminary Study in Sheep. Eur. J. Pharmaceutics Biopharmaceutics 85, 99–106. doi:10.1016/j.ejpb.2013.03.004

CrossRef Full Text | Google Scholar

Ma, D. K. G., Stolte, C., Krycer, J. R., James, D. E., and O’Donoghue, S. I. (2015). SnapShot: Insulin/IGF1 Signaling. Cell 161, 948. doi:10.1016/j.cell.2015.04.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Mcmillin, K. W., and Brock, A. P. (2005). Production Practices and Processing for Value-Added Goat Meat. J. Anim. Sci. 83, E57–E68.

Google Scholar

McMurtry, J. P. (1998). Nutritional and Developmental Roles of Insulin-like Growth Factors in Poultry. J. Nutr. 128, 302s–305s. doi:10.1093/jn/128.2.302S

PubMed Abstract | CrossRef Full Text | Google Scholar

Mikawa, S., Yoshikawa, G.-i., Aoki, H., Yamano, Y., Sakai, H., and Komano, T. (1995). Dynamic Aspects in the Expression of the Goat Insulin-like Growth Factor-I (IGF-I) Gene: Diversity in Transcription and post-transcription. Biosci. Biotechnol. Biochem. 59, 87–92. doi:10.1271/bbb.59.87

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, B. H., and Gore, A. C. (2001). Alterations in Hypothalamic Insulin-like Growth Factor-I and its Associations with Gonadotropin Releasing Hormone Neurones during Reproductive Development and Ageing. J. Neuroendocrinol 13, 728–736. doi:10.1046/j.1365-2826.2001.00686.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Naicy, T., Venkatachalapathy, T., Aravindakshan, T., Raghavan, K. C., Mini, M., and Shyama, K. (2017). cDNA Cloning, Structural Analysis, SNP Detection and Tissue Expression Profile of the IGF1 Gene in Malabari and Attappady Black Goats of India. J. Genet. 96, 307–312. doi:10.1007/s12041-017-0762-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Ni, F., Sun, R., Fu, B., Wang, F., Guo, C., Tian, Z., et al. (2013). IGF-1 Promotes the Development and Cytotoxic Activity of Human NK Cells. Nat. Commun. 4, 1479. doi:10.1038/ncomms2484

PubMed Abstract | CrossRef Full Text | Google Scholar

Patton, J., Kenny, D. A., McNamara, S., Mee, J. F., O’Mara, F. P., Diskin, M. G., et al. (2007). Relationships Among Milk Production, Energy Balance, Plasma Analytes, and Reproduction in Holstein-Friesian Cows. J. Dairy Sci. 90, 649–658. doi:10.3168/jds.S0022-0302(07)71547-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Scaramuzzi, R. J., Campbell, B. K., Downing, J. A., Kendall, N. R., Khalid, M., Muñoz-Gutiérrez, M., et al. (2006). A Review of the Effects of Supplementary Nutrition in the Ewe on the Concentrations of Reproductive and Metabolic Hormones and the Mechanisms that Regulate Folliculogenesis and Ovulation Rate. Reprod. Nutr. Dev. 46, 339–354. doi:10.1051/rnd:2006016

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, W., Wisniowski, P., Denne, S. C., Boyle, D. W., and Liechty, E. A. (2005). Anabolic Effects of Insulin and IGF-I in the Ovine Fetus Are Reduced by Prolonged Maternal Fasting. Am. J. Physiology-Endocrinology Metab. 288, E907–E913. doi:10.1152/ajpendo.00551.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Silva, J. R. V., Figueiredo, J. R., and van den Hurk, R. (2009). Involvement of Growth Hormone (GH) and Insulin-like Growth Factor (IGF) System in Ovarian Folliculogenesis. Theriogenology 71, 1193–1208. doi:10.1016/j.theriogenology.2008.12.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Stanišić, N., Žujović, M., Tomić, Z., Maksimović, N., Bijelić, Z., Ivanović, S., et al. (2012). The Effects of Crossing Balkan and Saanen Goat Breeds on Carcass Traits and Certain Quality Parameters of Kid Meat. Ann. Anim. Sci. 12, 53–62. doi:10.2478/v10220-012-0004-8

CrossRef Full Text | Google Scholar

Suske, G. (1999). The Sp-Family of Transcription Factors. Gene 238, 291–300. doi:10.1016/s0378-1119(99)00357-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Tao, L., He, X., Jiang, Y., Liu, Y., Ouyang, Y., Shen, Y., et al. (2021). Genome-Wide Analyses Reveal Genetic Convergence of Prolificacy between Goats and Sheep. Genes 12, 480. doi:10.3390/genes12040480

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomas, N., Venkatachalapathy, T., Aravindakshan, T., and Raghavan, K. C. (2016). Molecular Cloning, SNP Detection and Association Analysis of 5′ Flanking Region of the Goat IGF1 Gene with Prolificacy. Anim. Reprod. Sci. 167, 8–15. doi:10.1016/j.anireprosci.2016.01.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Velazquez, M. A., Newman, M., Christie, M. F., Cripps, P. J., Crowe, M. A., Smith, R. F., et al. (2005). The Usefulness of a Single Measurement of Insulin-like Growth Factor-1 as a Predictor of Embryo Yield and Pregnancy Rates in a Bovine MOET Program. Theriogenology 64, 1977–1994. doi:10.1016/j.theriogenology.2005.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Velazquez, M. A., Spicer, L. J., and Wathes, D. C. (2008). The Role of Endocrine Insulin-like Growth Factor-I (IGF-I) in Female Bovine Reproduction. Domest. Anim. Endocrinol. 35, 325–342. doi:10.1016/j.domaniend.2008.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Vella, A., Bouatia-Naji, N., Heude, B., Cooper, J. D., Lowe, C. E., Petry, C., et al. (2008). Association Analysis of the IGF1 Gene with Childhood Growth, IGF-1 Concentrations and Type 1 Diabetes. Diabetologia 51, 811–815. doi:10.1007/s00125-008-0970-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Vizcaíno, C., Mansilla, S., and Portugal, J. (2015). Sp1 Transcription Factor: A Long-Standing Target in Cancer Chemotherapy. Pharmacol. Ther. 152, 111–124. doi:10.1016/j.pharmthera.2015.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, S., Li, S., Zhang, J., Zhang, S., Dai, L., Bao, Y., et al. (2009). Cloning and Characterization of Class 1 and Class 2 Insulin-like Growth Factor-I mRNA in Songliao Black Pig. Mol. Biol. Rep. 36, 415–421. doi:10.1007/s11033-007-9195-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Yilmaz, A., Davis, M. E., Hines, H., and Chung, H. (2005). Detection of Two Nucleotide Substitutions and Putative Promoters in the 5' Flanking Region of the Ovine IGF-I Gene. J. Appl. Genet. 46, 307–309.

PubMed Abstract | Google Scholar