Quantitative proteomics analysis to assess protein expression levels in the ovaries of pubescent goats

Background Changes in the abundance of ovarian proteins play a key role in the regulation of reproduction. However, to date, no studies have investigated such changes in pubescent goats. Herein we applied isobaric tags for relative and absolute quantitation (iTRAQ) and liquid chromatography–tandem mass spectrometry to analyze the expression levels of ovarian proteins in pre-pubertal (n = 3) and pubertal (n = 3) goats. Results Overall, 7,550 proteins were recognized; 301 (176 up- and 125 downregulated) were identified as differentially abundant proteins (DAPs). Five DAPs were randomly selected for expression level validation by Western blotting; the results of Western blotting and iTRAQ analysis were consistent. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis indicated that DAPs were enriched in olfactory transduction, glutathione metabolism, and calcium signaling pathways. Besides, gene ontology functional enrichment analysis revealed that several DAPs enriched in biological processes were associated with cellular process, biological regulation, metabolic process, and response to stimulus. Protein–protein interaction network showed that proteins interacting with CDK1, HSPA1A, and UCK2 were the most abundant. Conclusions We identified 301 DAPs, which were enriched in olfactory transduction, glutathione metabolism, and calcium signaling pathways, suggesting the involvement of these processes in the onset of puberty. Further studies are warranted to more comprehensively explore the function of the identified DAPs and aforementioned signaling pathways to gain novel, deeper insights into the mechanisms underlying the onset of puberty. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08699-y.

development of animals [4]. Gonadotropin-releasing hormone (GnRH) neurons in the arcuate nucleus of the hypothalamus synthesize and secrete GnRH, and during the onset of puberty, their activity and consequently GnRH secretion are increased, which is characterized by pulsed GnRH release [5]. GnRH stimulates gonadotropic cells in the anterior pituitary gland to secrete folliclestimulating hormone and luteinizing hormone. These hormones act on the ovaries via blood circulation to facilitate the secretion of sex hormones, thereby promoting the rapid development of reproductive organs as well as sexual characteristics and playing a pivotal role in the onset of animal puberty [6,7]. Besides, the ovaries in turn regulate the hypothalamic-pituitary-ovarian axis by positive and negative feedback via the secreted sex hormones, further affecting puberty [8].
The mechanism underlying the onset of puberty is complex and reportedly involves several factors, including neuroendocrinological, genetic, and environmental factors [9]. The specific regulatory mechanism nevertheless remains unclear. Nguyen et al. performed adipose tissue proteomic analyses to study puberty in Brahman heifers; 51 significantly differentially abundant proteins (DAPs) were identified between pre-and post-pubertal heifers. These DAPs were enriched in estrogen signaling and PI3K-Akt signaling pathways, which are known integrators of metabolism and reproduction [10]. Further, Fortes et al. investigated how the uterine tissue and its secretion changes in relation to puberty in Brahman heifers. Using a combination of proteomics and transcriptomics, they identified 258 DAPs in the uterine fluid of pre-and post-pubertal cows [11]. Similarly, Ye et al. applied proteomics to examine the hypothalamus of pre-pubertal and pubertal female goats, which led to the identification of 69 DAPs. These proteins were enriched in the MAPK, RAS, and PI3K-Akt-mTOR signaling pathways, indicating that the identified DAPs and their related signaling pathways are crucial in regulating puberty in goats [12]. Collectively, the results of such studies indicate that proteomics can be effectively used to assess the hypothalamus as well as other tissues to explore the mechanisms underlying puberty.
Although the expression of retinol-binding protein 4 in the ovaries of female mice has been found to remain unchanged until puberty, a significant increase was noted during puberty [13]; moreover, retinol-binding protein 4 evidently plays a chief role in the transportation and storage of cells in follicular fluid [14]. Tahir et al. used whole ovaries of six pre-and six post-pubertal Brahman heifers to perform differential abundance analyses of protein profiles between the two physiological states, identifying 32 steroidogenesis-associated DAPs [15]. An increase in steroidogenesis is observed during puberty [16]. A study identified 36 proteins in ovarian follicular fluid of goats at different developmental stages [17]. Overall, these findings suggest that ovarian proteins play a fundamental role in pubertal development in animals.
Goats play an important socioeconomic role in the lives of people considering that they are a source of, for example, milk and meat. Puberty is a pivotal stage in female animal development; it marks the first occurrence of ovulation and the onset of reproductive capability [18]. A study reported that lncRNA is related to ovarian steroid hormone synthesis, oogenesis, and oocyte maturation during the growth and development of goats [19]; moreover, the levels of hormones in goats have been found to significantly change during puberty [20]. However, to date, no studies have used proteomics to assess changes in the expression level of proteins in the ovaries of pre-pubertal and pubertal goats. Herein we thus aimed to identify proteins and pathways associated with the onset of puberty by measuring protein abundance levels in the ovaries of prepubertal and pubertal goats. Figure 1 shows the experimental design and workflow. We believe our results should enhance our understanding of the key factors and mechanisms that regulate puberty in goats.

Protein Identification
Overall, 7,550 proteins were identified and quantified by isobaric tags for relative and absolute quantitation (iTRAQ) proteomics, of which 2,623 were uncharacterized proteins and 4,927 were proteins with known functions. In total, 301 DAPs were identified: 176 of them were up-and 125 were downregulated (Fig. 2a, b). On comparing data pertaining to pubertal with pre-pubertal goat ovaries, MHC class II antigen showed the highest relative upregulation, while DUF4537 domain-containing protein showed the highest relative downregulation. Tables 1 and 2 show the top 15 up-and downregulated proteins, respectively, in pubertal ovaries compared with pre-pubertal ones.

Cluster analysis of DAPs
Euclidean distance and hierarchical clustering were used to cluster DAPs, and dynamic changes in DAPs were compared (Fig. 4). We found significant differences in protein abundance intensity in the ovaries of pubertal and prepubertal goats, indicative of differences in protein abundance levels between pubertal and pre-pubertal goat ovaries.

Gene Ontology (GO) Functional Enrichment Analysis.
To further analyze the functions of the 301 DAPs, they were classified into the following major categories based on their GO annotations: biological process, cellular component, and molecular function. DAPs enriched in biological processes were mainly associated with cellular process, metabolic process, biological regulation, regulation of biological process, and response to stimulus (Fig. 5). Further, the five most abundant terms in the cellular component category were cell part, cell, organelle, membrane, and extracellular region, and those in the molecular function category were binding, catalytic activity, molecular function regulator, transcription regulator activity, and transporter activity.

Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis
To further explore the signaling pathways involved in regulating puberty onset in goats, DAPs were subjected to KEGG [21] pathway enrichment analysis (Fig. 6). The 301 DAPs were mapped to 260 KEGG pathways, including glutathione metabolism, olfactory transduction, and calcium signaling pathways, indicating their involvement with pubertal onset in goats.

Protein-Protein Interaction (PPI) network analysis
The interaction relationships among DAPs in the enriched signaling pathways (KEGG analysis) were found in the STRING database [22], and a PPI network was constructed (Fig. 7), which revealed some key interactions among DAPs. In the PPI network, nodes represent proteins and edges denote predicted functional associations; 222 nodes and 192 edges were found to be interconnected. CDK1, recombinant heat shock 70 kDa protein 1A (HSPA1A), and uridine cytidine kinase 2 (UCK2) occupied the central position in the PPI network and were observed to interact with other DAPs as a hub.

Discussion
Puberty is a phase of anatomical and physiological development, leading to sexual maturity and reproductive capacity [23]; however, the specific mechanisms underlying puberty remain unclear. Proteins are the material basis of life [24], so changes in the abundance of ovarian proteins are bound to affect pubertal onset. Therefore, we herein performed iTRAQ-LC-MS/MS analysis to compare the expression levels of ovarian proteins in prepubertal (n = 3) and pubertal (n = 3) goats so as to identify core proteins involved in puberty regulation. This led to the identification of 7,550 proteins, of which 301 DAPs met the screening standard.
Among the DAPs identified in this study, MHC class II antigen (fragment) was the most upregulated one in the ovaries of pubertal goats. MHC class II molecules are encoded by the MHC gene locus, which is located on chromosome 6 in humans and chromosome 17 in mice [25]. These molecules play a key regulatory role in several immune responses and are involved in most autoimmune diseases [26]. In addition, MHC genes reportedly influence mating preferences and odor in house mice [27] and even humans [28]. We believe that MHC class II antigen (fragment) affects puberty onset in goats by causing alterations in odor. RCAN1, previously known as DSCR1, was one of the most significantly upregulated proteins in this study. As it is located on human chromosome 21, RCAN1 has been postulated to contribute to mental retardation in Down syndrome. It is evidently involved in transcriptional regulation and/or signal transduction [29]; further, it is expressed in neurons and participates in local protein synthesis during neuronal stimulation [30]. Site-specific phosphorylation of RCAN1 has been reported to regulate bidirectional synaptic plasticity [31]. The downregulation of RCAN1 expression during radial migration of cortical neurons in rats has been observed to impair neural progenitor cell proliferation, resulting in radial migration defects [32]. Besides, both in vitro and in vivo studies have shown that RCAN1 overexpression promotes neuronal loss [33,34]. Nakata et al. found that growth hormone (GH) administration induced RCAN1 mRNA expression in the liver of hypophysectomized rats [35]. It is a known fact that the secretion rate and amount of GH increase during puberty [36], and in vitro experiment results have indicated that GH directly stimulates GnRH secretion from GT1-7 cells [37]. Therefore, it seems that RCAN1 indirectly regulates puberty onset via GH.
In this study, IGF-1 expression levels were also significantly upregulated in the ovaries of goats at puberty. IGF-1, a polypeptide hormone, is an important modulator of GH activity, and it plays a pivotal role in cell differentiation and proliferation [38] as well as in bone growth and development [39]. In addition, in the developing central nervous system, IGF-1 mediates changes in morphology, synaptic efficacy, and cellular organization, enabling cells of the central nervous system to respond to continuous external and internal stimuli [40]. IGF-1 is reportedly the upstream regulator of KiSS-1 [41]; moreover, IGF-1 is an important factor in reproductive and neuroendocrine function regulation, and it also participates in GnRH synthesis and release [42]. Intracerebroventricular administration of IGF-1 has been reported to promote GnRH secretion by activating KiSS-1 expression in the brain of pre-pubertal female rats, advancing the onset of puberty [43]. According to our findings, we believe that this phenomenon could be due to IGF-1 changing the sensitivity of KiSS-1 neurons to stimuli or changing the intensity and duration of neuronal synaptic activity. Further studies are warranted to explore the roles of RCAN1 and IGF-1 in the onset of puberty in goats.
Herein the enrichment of DAPs in olfactory transition, glutathione metabolism, and calcium signaling pathways suggested the involvement of these pathways in the onset of puberty in goats. Olfactory signals evidently play a key role in the onset of mammalian puberty [44]. The mammalian olfactory system can effectively detect pheromones; externally released steroids in mammals have been reported to function as pheromones and provide direct information pertaining to internal hormone state [45]. Olfactory sensory neurons are the first order neurons in the olfactory system and the initial site of odor detection [46]; the signals generated are transmitted to the hypothalamus, which affect reproduction as well as fertility [47]. Several pheromones affect pubertal onset in animals [48]. The olfactory system influences the timing of puberty in mice via GnRH neurons [49]. Furthermore, male scents can reportedly accelerate the onset of puberty in females, promoting uterine growth and the first estrus [50]. Glutathione is a ubiquitous tripeptide composed of glutamic acid, cysteine, and glycine, and it plays a chief role in cell proliferation regulation [51,52]. Glutamatecysteine ligase is the rate-limiting step in glutathione synthesis [53]. Glutamate dehydrogenase catalyzes glutamate synthesis; a study reported the hypothalamic content of glutamate to significantly increase in female rats at first proestrus [54], which may be related to the synthesis of glutathione, considering that the content of glutathione reaches its peak during this time [55]. In addition, subcutaneous administration of glutathione has been found to predate the onset of puberty in immature female rats, and the weight of the pituitary gland and ovary was found to significantly increase [56]. Similarly, intracerebroventricular administration of glutathione was reported to increase plasma levels of follicle-stimulating hormone [57], which is known to play a fundamental role in puberty [58]. From an overall perspective, it appears that glutathione and its metabolism play a chief role in regulating the onset of puberty.
Luteinizing hormone reportedly acts via protein kinase A and C to modulate T-type calcium currents and intracellular calcium transients in mice Leydig cells [59]. Calcium, a critical intracellular messenger, regulates numerous signal transduction pathways and is associated with steroid production [60]. In the excitable cells of the nervous system, an increase in free calcium levels evidently initiates neurotransmitter release [61], which consequently affects the onset of puberty in animals [62,63]. In this study, the expression level of phosphodiesterase (PDE) in the calcium signaling pathway was significantly upregulated. It is notable that the onset of puberty corresponds to the maturation of oocytes [64]. In most mammals, oocytes enter the early stage of meiosis during the fetal stage and their development then stops until ovulation or atresia [65]. The preovulatory gonadotropin surge leads to a decrease in cyclic guanosine monophosphate levels in oocytes and an increase in intracellular PDE activity, which leads to a decrease in cyclic adenosine monophosphate (cAMP) levels and recovery of meiosis [66,67]. PDE activity regulation is a complex, important approach to control su pra mo lec ula r co mp lex sy na ps e vir ion vir ion pa rt sy na ps e pa rt bin din g ca tal yti c ac tiv ity mo lec ula r fun cti on reg ula tor tra ns cri pti on reg ula tor ac tiv ity sig na l tra ns du ce r ac tiv ity str uc tur al mo lec ule ac tiv ity tra ns po rte r ac tiv ity mo lec ula r tra ns du ce r ac tiv ity mo lec ula r ca rrie r ac tiv ity an tio xid an t ac tiv ity the level of intracellular second messenger cAMP [68]. In addition, gonadotropins activate adenylyl cyclase in oocytes and facilitate its interaction with G protein to catalyze the conversion of adenosine triphosphate to cAMP [69,70]. In mammalian oocytes, cAMP maintains meiotic arrest by inhibiting maturation-promoting factor activation and stimulating cAMP-dependent protein kinase A [71]. This prevents oocytes from resuming meiosis prematurely, and they thus obtain normal development ability [72]. Herein we found that PDE and adenylate cyclase type 4 expression levels were significantly upregulated in goat ovaries during puberty. This could be a method to regulate intracellular cAMP levels so that the meiosis of oocytes and eventually the onset of puberty can occur normally. Our PPI network revealed that proteins interacting with CDK1, HSPA1A, and UCK2 were the most abundant. While HSPA1A and UCK2 have been discussed in the context of diseases such as cancer [73,74], their involvement in puberty remains to be reported. CDK1 is a cell cycle-related enzyme that drives cell division and is the main regulator of the cell cycle process [75,76]. In addition, mammalian cells express other cyclin kinases and phosphorylate cellular proteins to drive the cell cycle process [77,78]. Intriguingly, CDK1 can replace other CDKs and is sufficient to drive the mammalian cell cycle [79]. Artificial knockout of CDK1 has been found to lead to early embryo death, and conditional knockout of CDK1 in viable mice causes cell proliferation disorder [80]. Besides, CDK1 is essential and sufficient to drive the resumption of meiosis in mouse oocytes, and the process of meiosis is regulated by CDK1 activity [81,82]. In combination with the results of these studies, our findings imply that CDK1 plays a prominent role in mammalian reproduction.

Conclusions
To summarize, 301 DAPs were identified on analyzing ovarian protein expression levels in pubertal and pre-pubertal goats. We believe that the DAPs identified herein directly or indirectly regulate pubertal onset. Besides, olfactory conduction, glutathione metabolism, and calcium signaling pathways seem to play a significant role in the onset of puberty in goats. Further studies are nevertheless warranted to explore the specific molecular mechanisms underlying the onset of puberty.

Experimental Design.
Quantitative proteomics techniques based on iTRAQ were used to identify and characterize DAPs between prepubertal and pubertal goats. Goat ovaries were collected and suspended in protein lysate; ovarian proteins were then enzymatically digested with trypsin to obtain peptides. The hydrolyzed peptides were subsequently labeled with iTRAQ labels. Peptide separation and high-performance LC were performed, followed by MS and bioinformatics analyses. To verify the accuracy of our proteomics data, the expression levels of several proteins were assessed by WB.

Ovary collection
This study was authorized and endorsed by the Animal Care and Use Committee of Anhui Agricultural University. All experimental procedures involving goats were of puberty onset in female goats was defined as the date of the first estrus detected by male goats and by monitoring changes in the appearance of the vulva [83]. Rams test conditions were performed twice daily at 08:00 and 16:00. The cunnus of pubescent goats was inflamed and mature follicles were observed in the ovaries. The animals were deeply anesthetized by intramuscularly injecting 1% pentobarbital (60 mg/kg, Solarbio, P8410, China). The left ovary was then surgically removed and immediately frozen in liquid nitrogen, and they were stored at − 80 °C until needed.

Protein extraction
The ovaries were transferred into six 5-mL centrifuge tubes. An appropriate amount of 1 × cocktail was then added, followed by the addition of Ethylene Diamine Tetraacetic Acid (Ameresco, Shanghai, China) and sodium dodecyl sulfate (SDS, Solarbio, Beijing, China) lysis buffer 3. Finally, two 5-mm magnetic beads were added to the centrifuge tubes with tweezers to enable thorough grinding. After incubation on ice for 5 min, a tissue lyser (JXFSTPRP, Shanghai, China) was used to break and lyse the ovaries. The samples were then centrifuged at 25,000 × g and 4 °C for 15 min.
To the supernatant thus obtained, dithiothreitol (Bio-Rad, Hercules, CA, USA) was added at a final concentration of 10 mM, followed by incubation in a water bath at 56 °C for 1 h. Subsequently, iodoacetamide (Sigma-Aldrich, St. Louis, MO, USA) was added at a final concentration of 55 mM, and the suspension was incubated in the dark for 45 min. Five times the volume of precooled acetone (Lingfeng, Shanghai, China) was then added, followed by incubation at − 20 °C for 2 h and centrifugation at 25,000 × g and 4 °C for 15 min. The supernatant thus obtained was discarded, and a suitable amount of lysis buffer 3 without SDS was added to the precipitate to promote protein dissolution through a tissue lyser. After centrifugation at 25,000 × g and 4 °C for 15 min, the supernatant was obtained, which represented the protein extract.

Proteolysis and peptide labeling
To a 1.5-mL centrifuge tube, 100 μg protein was added. As per trypsin (μg):substrate protein (μg) ratio of 1:20, an enzyme solution was then added, and the solution was vortexed and then centrifuged at low speed for 1 min, followed by incubation at 37 °C for 4 h. The digested peptide solution was then subjected to desalination and freeze-dried. Peptides were reconstituted in 0.5 M tetraethylammonium bromide and processed with 8-plex iTRAQ reagent (Applied Biosystems), according to manufacturer instructions. Briefly, iTRAQ labeling reagents were thawed to room temperature, and 50 μL isopropyl alcohol was added to each reagent tube. The tube was then vortexed and shaken at low speed. Samples were labeled as follows: sample puberty 1 (113 tags), 2 (114 tags), and 3 (115 tags) and sample prepuberty 1 (118 tags), 2 (119 tags), and 3 (121 tags). After labeling, all samples were incubated at room temperature for 2 h. The labeled peptide mixtures were then pooled and dried by vacuum centrifugation.

Peptide Separation and High-Performance Liquid Chromatography.
The LC-20AB liquid phase system (Shimadzu, Kyoto, Japan) was used for liquid phase separation of samples. We used a Gemini C18 separation column (4.6 × 250 mm, 5 μm). iTRAQ-labeled peptide mixtures were reconstituted in mobile phase A (5% acetonitrile, pH 9.8) and injected; gradient elution was used at a flow rate of 1 mL/ min: 5% mobile phase B (95% acetonitrile, pH 9.8) for 10 min, 5%-35% mobile phase B for 40 min, 35%-95% mobile phase B for 1 min, mobile phase B for 3 min, and 5% mobile phase B for 10 min. Elution was monitored by measuring absorbance at 214 nm, and fractions were collected every 1 min. The eluted peptides were pooled into 20 fractions, followed by freeze drying.
The separation was performed using the Easy-nLC ™ 1200 system (Thermo Fisher Scientific, San Jose, CA). The freeze-dried peptide samples were reconstituted in mobile phase A (2% acetonitrile and 0.1% formic acid), centrifuged at 20,000 × g for 10 min, and the supernatant thus obtained was collected. The samples then entered a tandem self-packed C18 column (internal diameter, 75 μm; column size, 1.9 μm; column length, 25 cm) and were separated at a flow rate of 200 nL/min via the following gradient: 0-3 min, 5% mobile phase B (80% acetonitrile, 0.1% formic acid); 3-45 min, linear increase of mobile phase B from 8 to 44%; 45-50 min, linear increase of mobile phase B from 44 to 60%; 50-53 min, linear increase of mobile phase B from 60 to 100%; and 53-60 min, 80% mobile phase B. The nanoliter liquid phase separation end was directly connected to the mass spectrometer.

Mass spectrometry
The peptides separated in liquid phase were ionized by the nano-ESI source and then transferred to an Orbitrap Exploris 480 (Thermo Fisher Scientific, San Jose, CA) tandem mass spectrometer in data-dependent acquisition mode. The main settings were as follows: ion source voltage, 2.1 kV; scanning range of primary MS, 350-1,600 m/z with resolution of 60,000; and initial m/z of secondary MS, 100 m/z with resolution of 15,000. The screening conditions for the parent ions of secondary fragmentation were as follows: the charge was from 2 + to 7 + , and the peak strength was greater than 50,000 in the first 1 s. The ion fragmentation mode was higher energy collision-induced dissociation; fragment ions were detected in the Orbitrap, and the dynamic exclusion time was set at 30 s. The auto gain control was set to Level 1 1E6 and Level 2 1E5.

DAPs screening
The automated software IQuant [84], which integrates Mascot Percolator [85] algorithm, was used quantitatively analyze peptide samples labeled with isobaric tags. To assess the confidence of peptides, peptide-spectrum matches were pre-filtered at a false discovery rate (FDR) of 1%. Subsequently, based on the "parsimony principle, " identified peptide sequences were assembled into a set of confident proteins. To control the rate of false positives at the protein level, a protein FDR at 1%, which is based on Picked protein FDR strategy [86], was also estimated after protein inference (protein-level FDR ≤ 0.01). When the ratio of pubertal ovaries compared with pre-pubertal ovaries was 1.5 or 0.67, proteins were considered to up-or downregulated.

Western blot
Samples were obtained from the same ovaries as those used for iTRAQ-LC-MS/MS analysis, and protein concentration was determined using the BCA protein detection kit. After adding a protein loading buffer, the sample was boiled for 10 min at 100℃. Subsequently, 20 μg protein was subjected to SDS-polyacrylamide gel electrophoresis (Bio-Rad, Hercules, CA, USA), and protein bands thus obtained were transferred to a polyvinylidene fluoride (Merck Millipore, Billerica, MA, USA) membrane. The protein was sealed at room temperature for 2 h with 5% skimmed milk powder and then incubated with each of the following antibodies for overnight at 4℃: rabbit anti-RCAN1 (1:500, T73889, Abmart), rabbit anti-IGF1 (1:500, PA3635, Abmart), rabbit anti-DHCR24 (1:1000, TD12472, Abmart), rabbit anti-STMN1 (1:500, PA4263, Abmart), and rabbit anti-CDK1 (1:500, T55176, Abmart). The membranes were washed thrice with Tris-buffered saline with Tween (BioSharp, Hefei, China) for 10 min each time, and then a secondary antibody, i.e., horseradish peroxidase-labeled goat anti-rabbit immunoglobulin G (1:3,000, GB23303, Servicebio), was added, followed by incubation at room temperature for 1 h. After washing thrice with Trisbuffered saline with Tween, the membrane was finally analyzed using a chemiluminescence imaging analyzer. The test strips were excised, and their gray value was analyzed by the image analysis software ImageJ. β-Tubulin (1:2,500; Proteintech, Wuhan, China) was used as the internal control; the relative gray value of target proteins was obtained using the ratio of the gray values of target proteins and β-tubulin.

Bioinformatics
Hierarchical Cluster [87] and Java TreeView [88] were used for hierarchical clustering analysis. All genes corresponding to DAPs were subjected to GO annotation [89] and KEGG [90] pathway analyses. Gene symbols of DAPs were submitted to the STRING database [91], and Cytoscape v3.2.1 [92] was used to generate and analyze PPI networks.

Statistical analysis
GraphPad Prism 8 (GraphPad Software, San Diego, CA) was used for histogram plotting. Values represent mean ± SEM, and Student's t-test was used to compare relative expression levels of ovarian proteins between pubertal and pre-pubertal goats. IBM SPSS Statistics 26 (SPSS Inc., Chicago, USA) was used for data analysis. P < 0.05 indicates statistical significance.

Statement
This study was authorized and endorsed by the Animal Care and Use Committee of Anhui Agricultural University and reported in accordance with the ARRIVE guidelines. All experimental procedures involving goats were