Cathepsin D knockdown regulates biological behaviors of granulosa cells and affects litter size traits in goats

Cathepsin D knockdown regulates biological behaviors of granulosa cells and affects litter size traits in goats 组织蛋白酶D沉默调节颗粒细胞的生物学行为并影响山羊的产仔性状 first-author https://orcid.org/0000-0001-6078-938X corresp Xiang CHEN, XChen2@gzu.edu.cn Xiang CHEN, XChen2@gzu.edu.cn XChen2@gzu.edu.cn https://orcid.org/0000-0001-6495-6851 Key Laboratory of Plateau Mountain Animal Genetics, Breeding and Reproduction of Ministry of Education, Key Laboratory of Animal Genetics, Breeding and Reproduction of Guizhou Province, College of Animal Science, Guizhou University, Guiyang550025, China Key Laboratory of Plateau Mountain Animal Genetics, Breeding and Reproduction of Ministry of Education, Key Laboratory of Animal Genetics, Breeding and Reproduction of Guizhou Province, College of Animal Science, Guizhou University, Guiyang550025, China Cathepsin D (CTSD), the major lysosomal aspartic protease that is widely expressed in different tissues, potentially regulates the biological behaviors of various cells. Follicular granulosa cells are responsive to the increase of ovulation number, hence indirectly influencing litter size. However, the mechanism underlying the effect of CTSD on the behaviors of goat granulosa cells has not been fully elucidated. This study used immunohistochemistry to analyze CTSD localization in goat ovarian tissues. Moreover, western blotting was applied to examine the differential expression of CTSD in the ovarian tissues of monotocous and polytocous goats. Subsequently, the effects of CTSD knockdown on cell proliferation, apoptosis, cell cycle, and the expression of candidate genes of the prolific traits, including bone morphogenetic protein receptor IB ( BMPR-IB ), follicle-stimulating hormone ( FSHR ), and inhibin α ( INHA ), were determined in granulosa cells. Results showed that CTSD was expressed in corpus luteum, follicle, and granulosa cells. Notably, CTSD expression in the monotocous group was significantly higher than that in the polytocous group. In addition, CTSD knockdown could improve granulosa cell proliferation, inhibit cell apoptosis, and significantly elevate the expression of proliferating cell nuclear antigen ( PCNA ) and B cell lymphoma 2 ( Bcl-2 ), but it lowered the expression of Bcl-2-associated X ( Bax ) and caspase-3 . Furthermore, CTSD knockdown significantly reduced the ratios of cells in the G0/G1 and G2/M phases but substantially increased the ratio of cells in the S phase. The expression levels of cyclin D2 and cyclin E were elevated followed by the obvious decline of cyclin A1 expression. However, the expression levels of BMPR-IB, FSHR, and INHA clearly increased as a result of CTSD knockdown. Hence, our findings demonstrate that CTSD is an important factor affecting the litter size trait in goats by regulating the granulosa cell proliferation, apoptosis, cell cycle, and the expression of candidate genes of the prolific trait. 目的 2 探究组织蛋白酶D（Cathepsin D，CTSD）对颗粒细胞生物学行为及山羊产羔性状的调控机制。 创新点 2 以中国贵州省优良地方品种黔北麻羊为对象，以CTSD为候选基因，以卵泡颗粒细胞为模型，首次探明了CTSD沉默对山羊颗粒细胞生物学行为的影响及其对山羊产羔性状的调控机理，对进一步挖掘、创新和利用中国山羊品种资源等方面具有重要意义。 方法 2 使用免疫组化法对CTSD在卵巢中进行定位分析；应用蛋白质印迹法（western blotting）技术探究CTSD在单、多羔山羊卵巢中的表达差异；使用细胞计数试剂（cell counting kit-8，CCK-8）技术检测CTSD沉默对细胞增殖的影响；使用流氏细胞仪检测CTSD沉默对颗粒细胞凋亡及周期的影响；使用定量逆转录聚合酶链式反应（qRT-PCR）技术探究CTSD沉默对细胞增殖因子（增殖细胞核抗原（proliferating cell nuclear antigen，PCNA））、细胞凋亡标志基因（B细胞白血病/淋巴瘤2（B cell leukemia/lymphoma 2， Bcl-2 ）、Bcl-2相关X（Bcl-2-associated X， Bax ）、半胱天冬酶3（ caspase-3 ））、细胞周期蛋白（细胞周期蛋白A1（ cyclin A1 ）、细胞周期蛋白D2（ cyclin D2 ）、细胞周期蛋白E（ cyclin E ））及多胎性状候选基因（骨形态发生蛋白受体IB（bone morphogenetic protein receptor, type 1B， BMPR-IB ）、促卵泡素受体（follicle stimulating hormone receptor， FSHR ）、抑制素α（inhibin subunitα， INHA ））表达的影响；使用western blotting检测CTSD敲低对多胎性状候选基因 BMPR-IB 、 FSHR 、 INHA 在翻译水平表达的影响。 结论 2 CTSD沉默可能通过一系列因素的表达水平（包括 PCNA 、 Bcl-2 、 Bax 、 caspase-3 、 cyclin A1 、 cyclin D2 和 cyclin E ）来调控颗粒细胞的增殖、凋亡和细胞周期进展，进而影响颗粒细胞的生物学功能。此外，CTSD可能通过调控多产性状的候选基因，包括 BMPR-IB 、 FSHR 和 INHA 来影响卵泡发育和排卵，从而间接影响山羊的产仔数。 Cathepsin D (CTSD) Litter size trait Granulosa cells Cell apoptosis Cell cycle Cell proliferation 组织蛋白酶D（CTSD） 产羔性状 颗粒细胞 细胞凋亡 细胞周期 细胞增殖 1 Introduction 1 Multiple factors co-regulate the reproductive cap‑acity of goats; however, the reproductive characteristics of goat breeds or lines directly affect their reproductive phenotype and gene selection bias ( Lv et al., 2014 Lv et al., 2014 ; Bowman et al., 2018 Bowman et al., 2018 ). Of note, the "Qianbei Ma" goat is a fine local breed of Guizhou Province in China, with a litter size rate up to 196%. It exhibits early sexual maturity and year-round estrus cycle; thus it is considered as an ideal model for studying the reproductive traits of goats ( Cai et al., 2014 Cai et al., 2014 ; Sui et al., 2015 Sui et al., 2015 ). Litter size is a critical goat reproductive trait, and the ovary plays a vital role in lambing. Moreover, genetic differences of the ovary play a crucial role in variations in follicle development, ovulation, the secretion of female hormones, and litter size ( Soede et al., 2011 Soede et al., 2011 ; Pan et al., 2012 Pan et al., 2012 ; Hsueh et al., 2015 Hsueh et al., 2015 ). Therefore, the appropriate development of the ovary is essential for the reproductive capability of mammals, including goats. The follicle, which is the basic unit driving ovarian function, including growth and development ( McGee and Hsueh, 2000 McGee and Hsueh, 2000 ), continuously matures until ovulation. Its growth, development, and survival are regulated by multiple secretory hormones and growth factors ( Gougeon, 1996 Gougeon, 1996 ; Martins et al., 2014 Martins et al., 2014 ). For instance, granulosa cells have been shown to control follicular growth, development, and atresia by regulating proliferation, differentiation, cell cycle, and apoptosis ( Kranc et al., 2017 Kranc et al., 2017 ). Previously, different studies established the presence of bone morphogenetic protein receptor IB (BMPR-IB), follicle-stimulating hormone (FSHR) and inhibin α (INHA) in granulosa cells and associated them with prolific goat traits and ovulation number ( Bradford, 1972 Bradford, 1972 ; Mahdavi et al., 2014 Mahdavi et al., 2014 ; Cui et al., 2019 Cui et al., 2019 ; Xia et al., 2020 Xia et al., 2020 ). Notably, ovulation is also related to tissue remodeling; the follicular theca fuses with the tunica albuginea, and then thins. Afterwards, the rupture of follicle is triggered by a proteolytic enzyme in follicular fluid, which induces ovulation ( Boots and Jungheim, 2015 Boots and Jungheim, 2015 ). The lysosome contains proteolytic enzymes, such as cathepsin D (CTSD), which is a member of the lysosomal aspartic protease family. Aspartic protease is a proteolytic enzyme in the pepsin family and is a glycoprotein. CTSD was investigated in depth and compared to other aspartic proteases ( Dubey and Luqman, 2017 Dubey and Luqman, 2017 ). It was also reported as an important regulatory factor affecting the different stages of endogenous and exogenous cell apoptosis pathways ( Minarowska et al., 2007 Minarowska et al., 2007 ). CTSD induces cell apoptosis under cytotoxic conditions. Besides, it is involved in mitosis, which prompts the activation of leukocytes and regulates vascular wall permeability ( Minarowska et al., 2007 Minarowska et al., 2007 ). Moreover, the proteolytic activity of this enzyme is regulated by various factors in the lysosome, including pH value, metabolites, hormones, growth factors, and specific inhibitors ( Liu et al., 2016 Liu et al., 2016 ; Francisco et al., 2019 Francisco et al., 2019 ). In addition, CTSD exerts a significant role in the hydrolysis and conversion of various proteins in cell processes, and it is the ubiquitous lysosomal endoprotease of mammalian cells ( Minarowska et al., 2007 Minarowska et al., 2007 ). Simultaneously, CTSD also functions in protein degradation ( Ge et al., 2018 Ge et al., 2018 ). In fish, CTSD mediates the processing of vitellin in the oocyte; therefore, it is essential for living egg production ( Brooks et al., 1997 Brooks et al., 1997 ). Brooks et al. (1997) Brooks et al. (1997) have reported that CTSD messenger RNA (mRNA) is expressed in both ovarian and non-ovarian tissues (including liver, muscle, spleen, and testicles). During ovarian development, the highest expression level of CTSD mRNA occurs at the early stage of vitellogenesis, and then this level gradually declines ( de Stasio et al., 1999 de Stasio et al., 1999 ). Thus, it was suggested that CTSD contributes to the comprehensive growth and development of ovarian tissue, and has an essential function in determining animal reproductive traits ( de Stasio et al., 1999 de Stasio et al., 1999 ). A recent sequencing study identified the differential expression of CTSD in the ovaries of monotocous and polytocous sheep ( Feng et al., 2018 Feng et al., 2018 ). Our previous transcriptome sequencing findings showed that, in comparison with monotocous goats, CTSD expression was significantly down-regulated in the ovarian tissues of polytocous goats (unpublished data). Therefore, CTSD regulates follicular development and biological function in goat granulosa cells. However, no study has yet evaluated the effect of CTSD on the biological behavior of goat granulosa cells or the litter size trait. In this work, we (1) explored the correlation between CTSD expression and litter size in goats and the effects of CTSD on the biological behavior of granulosa cells, as well as the location of CTSD in goat ovarian tissue; (2) analyzed the differential expression of CTSD in ovarian tissues between monotocous and polytocous "Qianbei Ma" goats; and (3) revealed the effects of CTSD knockdown on cell proliferation, apoptosis, cell cycle, and the expression of candidate genes of the prolific trait in granulosa cells. 2 Materials and methods 1 2.1 Source of experimental animals 2 Ten "Qianbei Ma" goats (aged 36 months) with an initial average body weight of (38.6±2.8) kg (mean±standard deviation (SD)) were obtained for the experiment. Ewes with two consecutive litter size records were divided into two groups: a polytocous group ( n =5; litter size=3) and a monotocous group ( n =5; litter size=1). Both groups were housed under conditions of constant temperature (25 ℃) and fixed light/dark cycle (12 h/12 h) at the site of Fuxing Husbandry Co., Ltd. (Guizhou, China) and were provided with feed and water ad libitum. The housing had a high "bed" fence and all goats were clinically healthy. Before the beginning of the experiment, estrus was synchronously inducted in ewes. Then, a sponge suppository was implanted into the ewes' vagina for 13 d, after which the animals were injected 0.5 mg cloprostenol, and the sponge was retrieved. Two days after injection, the ewes were checked for signs of estrus at 10:00 am and 4:00 pm. After confirming the presence of estrus, the ewes were individually slaughtered, and ovarian samples were collected and immediately preserved in liquid nitrogen awaiting RNA and protein extraction. Furthermore, the right ovaries of polytocous "Qianbei Ma" goats were collected and washed with sterile phosphate-buffered saline (PBS) containing antibiotics (penicillin and streptomycin) for further storage and the isolation of follicular granulosa cells ( > 2.5 mm). Subsequently, the right ovarian tissues used for immunohistochemistry were fixed in a 4% (volume fraction) paraformaldehyde (PFA) solution and paraffin-embedded. 2.2 Immunohistochemistry 2 The immunohistochemistry protocol was performed according to that described previously ( Peng et al., 2016 Peng et al., 2016 ). In short, paraffin sections with 5-‍μm thickness were attached to microscope slides, heated at 65 ℃ for 2 h, dewaxed in xylene, and then rehydrated in a graded series of ethanol. Subsequently, these sections were placed in 0.1 mol/L sodium citrate buffer (pH 6.0), boiled in the microwave for 15 min, and then cooled to 37 ℃. After rinsing with PBS, the activity of endogenous peroxidase was blocked by immersing the sections in 3% (mass fraction) hydrogen peroxide (H 2 O 2 ) for 15 min. To block the non-specific staining process, the sections were immersed in a blocking solution (5% (mass fraction) bull serum albumin (BSA) in TBSTw (including Tris-HCl, NaCl, and Tween 20)) at 37 ℃ for 15 min. The sections were then immersed for 12 h at 4 ℃ with rabbit anti-CTSD polyclonal antibody (1:‍2000 (volume ratio, the same below), Bioss, bs-1500R, Beijing, China) in TBSTw. After rinsing with PBS, the sections were immersed with horse radish peroxidase (HRP)‍-labelled secondary antibody (A0208, Beyotime Biotechnology, Shanghai, China) at 37 ℃ for 2 h and then with diaminobenzidine for 15 min. The primary antibody was replaced with non-specific rabbit immunoglobulin G (IgG) to act as a negative control (NC). Digital images of immunohistochemically stained sections were captured through three light source microscopes (80i, Nikon, Japan). 2.3 RNA interference 2 In order to down-regulate CTSD expression, we used four CTSD short hairpin RNAs (shRNAs), including sh1-CTSD, sh2-CTSD, sh3-CTSD, and sh4-CTSD, and negative shRNA control (sh-NC), which were designed and synthesized by GenePharma (Shanghai, China). The information of sh-CTSD and sh-NC fragments is shown in Table 1 Table 1 . Subsequently, we extracted the sh-CTSD plasmids and the corresponding control plasmid (TIANGEN Biotech, Beijing, China) and stored them at -‍20 ℃ until transfection. We verified the suppression efficiency of the CTSD gene at the transcription and translation levels. shRNA: short hairpin RNA; CTSD: cathepsin D; sh-CTSD: CTSD shRNA; sh-NC: negative shRNA control. Fragment Fragment sequences (5'→3') sh1-CTSD Forward: CACCGGATCCCACTGCAGAAGTTCATTCAAGAGATGAACTTCTGCAGTGGGATCCTTTTTTG Reverse: GATCCAAAAAAGGATCCCACTGCAGAAGTTCATCTCTTGAATGAACTTCTGCAGTGGGATCC sh2-CTSD Forward: CACCGCTCAAGAACTACATGGACGCTTCAAGAGAGCGTCCATGTAGTTCTTGAGCTTTTTTG Reverse: GATCCAAAAAAGCTCAAGAACTACATGGACGCTCTCTTGAAGCGTCCATGTAGTTCTTGAGC sh3-CTSD Forward: CACCGGGTCCACCACAAATACAACATTCAAGAGATGTTGTATTTGTGGTGGACCCTTTTTTG Reverse: GATCCAAAAAAGGGTCCACCACAAATACAACATCTCTTGAATGTTGTATTTGTGGTGGACCC sh4-CTSD Forward: CACCGCCTGTCTTCGACAACCTGATTTCAAGAGAATCAGGTTGTCGAAGACAGGCTTTTTTG Reverse: GATCCAAAAAAGCCTGTCTTCGACAACCTGATTCTCTTGAAATCAGGTTGTCGAAGACAGGC sh-NC Forward: CACCGTTCTCCGAACGTGTCACGTTTCAAGAGAACGTGACACGTTCGGAGAATTTTTTG Reverse: GATCCAAAAAATTCTCCGAACGTGTCACGTTCTCTTGAAACGTGACACGTTCGGAGAAC 2.4 Cell culture and plasmid transfection 2 The preparation procedure followed the method proposed by Peng et al. (2016) Peng et al. (2016) . The ovarian tissue was placed on ultra-clean workbench. Then, we pierced the follicle with a 10-gauge needle, which permitted follicular fluid flow into Dulbecco's modified Eagle medium/nutrient mixture F-12 (DMEM/F-12). This was followed by 10-min centrifugation at 1000 r/min. After discarding the supernatant, the pellet was mixed with DMEM/F-12 medium by pipetting and transferred to a cell culture flask for culture (37 ℃, 5% CO 2 ). Subsequently, Lipofectamine TM 2000 (Thermo Fisher Technology, Massachusetts, USA) was used to transiently transfect granulosa cells with 15 μg sh1-CTSD to sh4-CTSD or sh-NC. Cells were harvested 48 h post-transfection ( Zhang et al., 2018 Zhang et al., 2018 ). 2.5 RNA extraction and qRT-PCR analysis 2 TRIzol ® reagent (Invitrogen, Grand Island, USA) was used for total RNA extraction from cultured follicular granulosa cells. Exactly 500 ng of total RNA was obtained for complementary DNA (cDNA) first-strand synthesis by the StarScript II First-strand cDNA Kit (GenStar, Beijing, China). We performed quantitative real-time polymerase chain reaction (qRT-PCR) based on the SYBR fluorescent dye method to quantify CTSD mRNA abundance in follicular granulosa cells (GenStar). A list of the PCR primers used is provided in Table 2 Table 2 . For all reactions, a 20-μL final volume of reaction mixture was prepared, comprising of 10 μL 2× SYBR mixture, 100 nmol/L of each primer, and 1 μL cDNA. The cycling parameters were as follows: one cycle at 95 ℃ for 2 min, followed by 40 cycles at 95 ℃ for 15 s, 60 ℃ for 30 s, and 72 ℃ for 30 s, extended by one cycle at 72 ℃ for 5 min. CTSD : cathepsin D; BMPR-IB : bone morphogenetic protein receptor IB; INHA : inhibin α; FSHR : follicle-stimulating hormone; Bcl-2: B cell lymphoma 2; Bax : Bcl-2-associated X; PCNA : proliferating cell nuclear antigen; PCR: polymerase chain reaction. Gene Primer sequence (5'→3') Product (bp) Amplification efficiency (%) GenBank ID CTSD Forward: TGACTCGCATGTAGGTTGC Reverse: TGGAGACGCTGTAGGAA 185 98.15 XM_018043344.1 BMPR-IB Forward: GCTCTTGGTCCTCATCATTTTATTC Reverse: ATGTAAGTTTCGTCCTGTTCTAACC 99 100.87 XM_013964509.2 INHA Forward: GACCAAGATGTCTCCCAG Reverse: CAGTATGGAACCACAGCT 166 98.98 NM_001285606.1 FSHR Forward: TGTTATGTCCCTCCTTGTGCTC Reverse: CGCTTGGCTATCTTGGTGTCA 129 101.15 NM_001285636.1 Bcl-2 Forward: TCTTTGAGTTCGGAGGGGTC Reverse: TGGTGGAAGAGTGTGCTGTG 188 102.39 XM_018039337.1 Bax Forward: CATCGGAGATGAATTGGACAGTAA Reverse: GGCCTTGAGCACCAGTTTGC 178 100.26 XM_013971446.2 caspase-3 Forward: CCTGGACTGTGGTATTGAGA Reverse: TAACCCGAGTAAGAATGTGC 218 99.12 XM_018041755.1 cyclin A1 Forward: GACAGTGCCAACTACCAACC Reverse: GGTATTTCAAGAACGGGTCA 131 96.77 XM_018056659.1 cyclin D2 Forward: GAGCAGAAGTGCGAAGAGGAGG Reverse: TTGATGGAGTTGTCGGTGTAAATG 191 104.54 XM_005680985.3 cyclin E Forward: GATGTCGGCTGCTTAGAAT Reverse: CACCACTGATACCCTGAAAC 104 103.31 XM_018062248.1 PCNA Forward: GAAGAAAGTGCTGGAGGC Reverse: TCGGAGCGAAGGGTTA 129 99.93 XM_005688167.3 β-actin Forward: AGATGTGGATCAGCAAGCAG Reverse: CCAATCTCATCTCGTTTTCTG 139 97.23 NM_001314342.1 2.6 Western blotting 2 Based on a previously described method ( Ao et al., 2019 Ao et al., 2019 ), follicular granulosa cells grown to 90% confluency and ovarian tissues were lysed with radio immunoprecipitation assay (RIPA) buffer (Beyotime Biotechnology, Shanghai, China). The concentration of extracted total protein was determined by the BCA protein assay kit (Beyotime Biotechnology). Sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and transmembrane electrophoresis were performed with the Mini-PROTEAN Tube Cell Instrument (Bio-Rad, Hercules, CA). The membranes were incubated with primary antibodies, including anti-CTSD (1:2000; Bioss, bs-1615R, Beijing, China), anti-BMPR-IB (1:1500; Bioss, bs-6639R), anti-FSHR (1:1500; Bioss, bs-0895R), anti-INHA (1:1000, Bioss, bs-1032R), and anti-‍β‍-‍actin (1:3000; Affinity Biosciences, AF7018, USA) ( Tang et al., 2019 Tang et al., 2019 ). Next, the membranes were subjected to further incubation with HRP-conjugated goat anti-rabbit secondary antibody (A0208, Beyotime Biotechnology). The immunoreactive bands were visualized using an electrochemiluminescence (ECL) kit (Beyotime Biotechnology). The membrane images were documented using a chemiluminescence imaging system (Bio-Bad, USA), whereas ImageJ (V1.48) was adopted to quantify the densitometric analysis result for the bands. All blots were normalized to β‍-actin to assess the relative expression of target proteins. 2.7 Cell proliferation 2 Using the cell counting kit-8 (CCK-8; Beyotime Biotechnology), the proliferation of granulosa cells was evaluated after CTSD knockdown ( Zhang et al., 2019 Zhang et al., 2019 ). Granulosa cells transfected with sh-NC or sh-CTSD were seeded into 96-well plates and incubated for increasing durations (0, 6, 12, 24, 48, and 72 h) before adding CCK-8 solution (10 µL) into each well for 3 h. Absorbance measurements (at 450 nm) were taken using a microplate reader (Thermo Fisher Technology). Data for the CCK-8 assays were evaluated using a Student's t -test, and were presented as mean±SD. P < 0.05 was the significance threshold. 2.8 Cell apoptosis and cell cycle 2 The cell apoptosis analysis kit (C1052, Beyotime Biotechnology) was used to detect apoptosis while rigorously following the manufacturer's instructions. Briefly, the harvested cells were washed with PBS, and then resuspended in 500 mL annexin-binding buffer. Annexin V-fluoresceine isothiocyanate (FITC) (4 mL) and propidium iodide (pI) (8 mL) were subsequently added, followed by 5 min incubation at room temperature in the dark. Cells were then subjected to flow cytometry and analyzed by FlowJo software (FACS Calibur, Becton Dickinson, CA, USA). Moreover, the goat granulosa cells (1×10 5 cells/mL) were detected using a cell cycle kit (KGA512, KeyGEN Biotech, Jiangsu, China). These cells were transfected with sh-NC or sh-CTSD for 48 h. The treated cells were harvested, washed with PBS, and then fixed in 75% (volume fraction) ethanol at 4 ℃ for 16 h. Following centrifugation to remove residual ethanol, cells were washed once with PBS, and incubated in 2 mL staining buffer (A) and 10 mL of reagent B for 20 min at room temperature in the dark. The stained cells were examined using a flow cytometer (Guava easyCyte 6-2L; Millipore, USA). Data were analyzed using the ModFit LT 5.0 software (Verity Software House, USA). 2.9 Statistical analysis 2 SPSS 19.0 software was employed to analyze the differences in cell proliferation, apoptosis, and cycle data between the sh-CTSD group and the sh-NC group. A Student' s t -test was used to evaluate the differences in CTSD protein expression levels between the monotocous and polytocous ovaries. All data represent the mean±SD. The significance threshold was set as P < 0.05. In this study, three biological repeats and three technical repeats were performed to ensure the accuracy of experimental data. 3 Results 1 3.1 Location of CTSD in goat ovarian tissues 2 Immunohistochemistry analysis was performed to determine the expression pattern and distribution of CTSD in ovarian tissues. CTSD was revealed to exist in oocytes from all levels of follicles ( Figs. 1 Figs. 1 a‍–‍1e). Additionally, the immunoreaction of CTSD protein was found in corpus luteum and follicular granulosa cells. For all the detected antibodies, the reaction of control group ( Fig. 1 Fig. 1 f) confirmed the nonspecific staining process. 3.2 Expression of CTSD in goat ovarian tissues 2 Western blotting was employed to confirm the differential expression of CTSD in ovarian tissues of goats with different litter sizes, i.e., to compare CTSD expression in the ovary of monotocous and polytocous goats. We found greater CTSD protein expression in the ovary of monotocous goats than in the ovary of polytocous goats ( P < 0.01; Fig. 2 Fig. 2 ). 3.3 Suppression of CTSD and expression of candidate genes of the prolific trait in goats 2 To elucidate the influence of CTSD on the goat litter size trait, CTSD expression in granulosa cells was decreased through transfection with specific shRNA. We estimated the changes of candidate genes of the prolific trait in granulosa cells, including BMPR-IB , FSHR , and INHA . Results revealed that, among the four oligonucleotide sequences of sh-CTSD, sh3-CTSD produced the most optimal interference effect on CTSD at the transcriptional level ( Fig. 3 Fig. 3 a). Simultaneously, it could significantly ( P < 0.01) inhibit the expression of CTSD protein at the translational level ( Fig. 3 Fig. 3 c). In addition, the downregulation of CTSD expression in granulosa cells could significantly improve the expression of BMPR-IB, FSHR, and INHA at both transcriptional and translational levels (Figs. 3b 3b and 3d 3d ). 3.4 Inhibition of CTSD, cell proliferation and apoptosis 2 In order to investigate the effects of CTSD on the biological behavior of follicular granulosa cells, and its expression levels of related genes, including proliferating cell nuclear antigen ( PCNA ), B cell lymphoma 2 ( Bcl-2 ), Bcl-2-associated X ( Bax ), and caspase-3 , specific shRNA was adopted to knock down CTSD expression in granulosa cells. The impacts of CTSD expression on granulosa cell proliferation and apoptosis were evaluated via the combination of CCK-8 colorimetry and flow cytometry. Compared to the granulosa cells treated with sh-NC, the proliferation activity of granulosa cells was increased after the knockdown of CTSD expression. Besides, at 12, 24, 48, and 72 h after shRNA transfection, CTSD knockdown significantly improved the proliferation of granulosa cells ( Fig. 4 Fig. 4 a) and increased the expression levels of PCNA genes (Fig. ‍4b). Compared with sh-NC, sh3-CTSD demonstrated a stronger inhibition effect on granulosa cell apoptosis after transfection ( P < 0.01;Fig. ‍4c), wherein the apoptosis rate was decreased by 1.3% ( Fig. 4 Fig. 4 d) and the level of cell anti-apoptotic factor Bcl-2 was significantly higher in granulosa cells after CTSD knockdown. Meanwhile, the expression levels of cell pro-apoptotic factors, Bax and caspase-3 , were obviously reduced ( P < 0.05; Fig. 4 Fig. 4 e). 3.5 Inhibition of CTSD on the cycle of follicular granulosa cells 2 The effects of CTSD knockdown on follicular granulosa cell cycle were explored through annexin-V and PI staining of the granulosa cells treated with sh3-CTSD and subsequent detection by flow cytometry. We found that CTSD knockdown could significantly lower the ratio of granulosa cells in the G0/G1 and G2/M phases ( P < 0.05), whereas the ratio of cells in the S phase was obviously increased ( P < 0.01) (Figs. 4f 4f and 4g 4g ). In addition, to confirm the effect of CTSD on granulosa cell cycle-related factors, the in vitro expression levels of cyclin A1, cyclin D2, and cyclin E were analyzed. The results following CTSD knockdown demonstrated significantly lower expression of the conversion factor cyclin A1 from the S phase to the G2/M phase. In contrast, the expression levels of a crucial cell cycle protein cyclin D2 in the G1 phase and the key conversion factor cyclin E in cells from the G0/G1 phase to the S phase were significantly and extremely high ( P < 0.01; Fig. 4 Fig. 4 h). 4 Discussion 1 Notably, specificity was detected in the gene expression of different tissues and cells. Since the ovary contains a mixture of different cells, the expression of the candidate genes in these cells is not the same ( Zhang et al., 2015 Zhang et al., 2015 ). Morais et al. (2012) Morais et al. (2012) described the expression pattern of CTSD in the ovary, and detected CTSD expression in the follicular cells at each stage of follicular growth and development. Moreover, the immunoreaction of CTSD was found in ovarian tissues of rats ( Zhang et al., 2008 Zhang et al., 2008 ). In the present study, we demonstrated that CTSD exists in the prim‑ordial, primary, and secondary follicles of "Qianbei Ma" goats and triggers positive immunoreaction in the corpus luteum and granulosa cells. This phenomenon verifies that CTSD is widely distributed in each part of goat ovary, whereas the positive expression of CTSD in each level of follicles demonstrates its potential contribution to the growth and development of goat follicle. Furthermore, the immunological competences of CTSD at all levels of follicle were distinct, whereas the immunocompetence in corpus luteum and big antral follicle before ovulation was higher. This result indicates that CTSD might be involved in follicular atresia in goats, as reported in other studies with different animal models ( de Stasio et al., 1999 de Stasio et al., 1999 ; Leyria et al., 2018 Leyria et al., 2018 ; Sales et al., 2019 Sales et al., 2019 ). Therefore, the detected expression pattern of CTSD in the goat ovary might offer novel perspectives for exploring the CTSD function and the regulatory mechanism of reproductive traits in goats. This study demonstrated that CTSD expression in the ovaries of monotocous "Qianbei Ma" goats was significantly higher compared with the ovary of poly‑tocous goats, as detected through western blotting. Similar results were reported by Feng et al. (2018) Feng et al. (2018) . In recent years, the comparison of gene expression pattern between the ovaries of females with high or low reproductive capacity has aided in screening for genes related to litter size ( Ling et al., 2014 Ling et al., 2014 ; Hu et al., 2020 Hu et al., 2020 ). The differential expression of CTSD gene in the ovarian tissues of monotocous and polytocous goats indicated that CTSD might have a high correl‑ation with goat reproductive capacity, especially litter size. Previous research on the litter size trait has indicated that the differences in this trait are determined by hormone and genetic related genes ( Yang et al., 2020 Yang et al., 2020 ). Numerous studies proved that the expression of BMPR-IB , FSHR , and INHA positively correlates with prolificacy in goats, and they may promote ovulation by regulating follicle development ( Bradford, 1972 Bradford, 1972 ; Mahdavi et al., 2014 Mahdavi et al., 2014 ; Cui et al., 2019 Cui et al., 2019 ; Xia et al., 2020 Xia et al., 2020 ). In this study, we established that knocking down the expression of CTSD in granulosa cells significantly elevated the expression of BMPR-IB , FSHR , and INHA at both transcriptional and translational levels. Thus, we inferred that CTSD expression could inhibit the expression of BMPR-IB, FSHR, and INHA in granulosa cells, consequently impeding follicle growth and development, and in turn indirectly reducing the number of ovulations. Follicular growth is characterized by the proliferation of granulosa cells and the enlargement of follicular volume. Meanwhile, follicle development is related to cell apoptosis, cell cycle, and the formation of vascular wall ( Chou and Chen, 2018 Chou and Chen, 2018 ). CTSD is an extensively expressed aspartic protease associated with several biological processes, including cell proliferation and apoptosis ( Tulone et al., 2007 Tulone et al., 2007 ). The PCNA is a widely active cell proliferation marker in all cell types (including follicular granulosa cells), whose expression potentially improves cell proliferation ( Strzalka and Ziemienowicz, 2011 Strzalka and Ziemienowicz, 2011 ). Herein, the CCK-8 method was employed to reveal that CTSD knockdown in granulosa cells could enhance granulosa cell proliferation. Simultaneously, PCNA expression was significantly increased, indicating that CTSD expression could inhibit cell proliferation. Thus, we hypothesized that CTSD-induced cell proliferation might be associated with the expression level of PCNA . Cell cycle regulation is extremely important for cell proliferation and differentiation, with cyclins being the critical regulatory factors for cell cycle progression ( Chen et al., 2016 Chen et al., 2016 ). Through flow cytometry, we indicated that CTSD knockdown in granulosa cells significantly lowered the ratios of cells in the G0/G1 and G2/M phases, while it increased the ratio of cells in the S phase, indicating that CTSD knockdown could block granulosa cells in the S phase. We thus demonstrated that CTSD knockdown might improve granulosa cell proliferation, which was consistent with the CCK-8 detection results. Simultaneously, the expression levels of the important cell cycle regulatory factors cyclin D2 and cyclin E were significantly elevated. On the other hand, the expression of cyclin A1 , which has been shown to play a crucial role in cell cycle progression, was obviously decreased ( Fuchimoto et al., 2001 Fuchimoto et al., 2001 ; Pradeep et al., 2002 Pradeep et al., 2002 ). Therefore, we speculated that CTSD potentially regulates the proliferation and differentiation of granulosa cells by varying the expression levels of cyclin A1 , cyclin D2 , and cyclin E , subsequently altering the progression and other functions of granulosa cells. Cell apoptosis, a programmed cell death mechanism with distinctive characteristics, involves unique factors, such as members of the Bcl-2 protein family (including anti-apoptosis member Bcl-2 and pro-apoptosis member Bax). Notably, caspase-3 of the caspase factor family is the most important effector molecule in cell apoptosis, and is expressed in ovarian tissue and related cells ( Greenfeld et al., 2007 Greenfeld et al., 2007 ; Han et al., 2017 Han et al., 2017 ). According to a previous research report, cell apoptosis is the primary mechanism for ovarian follicular atresia, and thus plays a crucial function in the ovarian follicle ( Liu et al., 2018 Liu et al., 2018 ). This study established that the inhibition of CTSD activity could decrease granulosa cell apoptosis and elevate the expression of anti-apoptosis factor Bcl-2 . In contrast, the expression of pro-apoptosis factors Bax and caspase-3 was reduced. These findings indicated that, by changing the ratio between Bax and Bcl-2 in the Bcl-2 family, CTSD may cause mitochondrial dysfunction and further release the cell terminal apoptotic factor caspase-3, inducing the apoptosis of granulosa cells. Importantly, cell apoptosis induced by lysosomal cathepsin induces follicular atresia in mammals, whereas the increased CTSD activity potentially causes atresia in follicles ( Gwon et al., 2017 Gwon et al., 2017 ). The present study demonstrated that CTSD knockdown lowers granulosa cell apoptosis, and therefore, high CTSD expression is the essential reason for low productivity in goats. 5 Conclusions 1 This study demonstrated that CTSD knockdown potentially regulates granulosa cell prolifer‑ation, apoptosis, and cell cycle progression through a range of factors, including PCNA , Bcl-2 , Bax , caspase-3 , cyclin A1 , cyclin D2 , and cyclin E , which in turn affects the biological functions of granulosa cells. Besides, CTSD might affect follicle development and ovulation by regulating the candidate genes of the prolific trait, including BMPR-IB , FSHR , and INHA , thus indirectly influencing goat litter size. References Ao Z , Wu X , Zhou J , et al. , 2019 . Cloned pig fetuses exhibit fatty acid deficiency from impaired placental transport . Mol Reprod Dev , 86 ( 11 ): 1569 - 1581 . https://doi.org/10.1002/mrd.23242 https://doi.org/10.1002/mrd.23242 Cloned pig fetuses exhibit fatty acid deficiency from impaired placental transport Boots CE , Jungheim ES , 2015 . Inflammation and human ovarian follicular dynamics . Semin Reprod Med , 33 ( 4 ): 270 - 275 . https://doi.org/10.1055/s-0035-1554928 https://doi.org/10.1055/s-0035-1554928 Inflammation and human ovarian follicular dynamics Bowman LL , Kondrateva ES , Timofeyev MA , et al. , 2018 . Temperature gradient affects differentiation of gene expression and SNP allele frequencies in the dominant Lake Baikal zooplankton species . Mol Ecol , 27 ( 11 ): 2544 - 2559 . https://doi.org/10.1111/mec.14704 https://doi.org/10.1111/mec.14704 Temperature gradient affects differentiation of gene expression and SNP allele frequencies in the dominant Lake Baikal zooplankton species Bradford GE , 1972 . Genetic control of litter size in sheep . J Reprod Fertil Suppl , 15 : 23 - 41 . Genetic control of litter size in sheep Brooks S , Tyler CR , Carnevali O , et al. , 1997 . Molecular characterisation of ovarian cathepsin D in the rainbow trout, Oncorhynchus mykiss . Gene , 201 ( 1-2 ): 45 - 54 . https://doi.org/10.1016/s0378-1119(97)00426-5 https://doi.org/10.1016/s0378-1119(97)00426-5 Molecular characterisation of ovarian cathepsin D in the rainbow trout, Oncorhynchusmykiss Cai HF , Chen Z , Luo WX , 2014 . Associations between polymorphisms of the GFI1B gene and growth traits of indigenous Chinese goats . Genet Mol Res , 13 ( 1 ): 872 - 880 . https://doi.org/10.4238/2014.February.13.5 https://doi.org/10.4238/2014.February.13.5 Associations between polymorphisms of the GFI1B gene and growth traits of indigenous Chinese goats Chen FL , Wang N , Yang DQ , et al. , 2016 . Herp depletion arrests the S phase of the cell cycle and increases estradiol synthesis in mouse granulosa cells . J Reprod Dev , 62 ( 2 ): 159 - 166 . https://doi.org/10.1262/jrd.2015-120 https://doi.org/10.1262/jrd.2015-120 Herp depletion arrests the S phase of the cell cycle and increases estradiol synthesis in mouse granulosa cells Chou CH , Chen MJ , 2018 . The effect of steroid hormones on ovarian follicle development . Vitam Horm , 107 : 155 - 175 . https://doi.org/10.1016/bs.vh.2018.01.013 https://doi.org/10.1016/bs.vh.2018.01.013 The effect of steroid hormones on ovarian follicle development Cui ZF , Liu LB , Zhao XL , et al. , 2019 . Analysis of expression and single nucleotide polymorphisms of INHA gene associated with reproductive traits in chickens . BioMed Res Int , 2019 : 8572837 . https://doi.org/10.1155/2019/8572837 https://doi.org/10.1155/2019/8572837 Analysis of expression and single nucleotide polymorphisms of INHA gene associated with reproductive traits in chickens de Stasio R , Borrelli L , Kille P , et al. , 1999 . Isolation, characterization and molecular cloning of cathepsin D from lizard ovary: changes in enzyme activity and mRNA expression throughout ovarian cycle . Mol Reprod Dev , 52 ( 2 ): 126 - 134 . https://doi.‍org/10.1002/(SICI)1098-2795(199902)52:‍2 https://doi.‍org/10.1002/(SICI)1098-2795(199902)52:‍2 < 126: :AID-MRD2 > 3.0.CO;2-O Isolation, characterization and molecular cloning of cathepsin D from lizard ovary: changes in enzyme activity and mRNA expression throughout ovarian cycle Dubey V , Luqman S , 2017 . Cathepsin D as a promising target for the discovery of novel anticancer agents . Curr Cancer Drug Targets , 17 ( 5 ): 404 - 422 . https://doi.org/10.2174/1568009616666161229145115 https://doi.org/10.2174/1568009616666161229145115 Cathepsin D as a promising target for the discovery of novel anticancer agents Feng X , Li FZ , Wang F , et al. , 2018 . Genome-wide differential expression profiling of mRNAs and lncRNAs associated with prolificacy in Hu sheep . Biosci Rep , 38(2):BSR20171350. https://doi.org/10.1042/BSR20171350 https://doi.org/10.1042/BSR20171350 Genome-wide differential expression profiling of mRNAs and lncRNAs associated with prolificacy in Hu sheep Francisco CC , Luis CLJ , Marina EBJ , et al. , 2019 . Effect of temperature and pH on the secondary structure and denaturation process of jumbo squid hepatopancreas cathepsin D . Protein Pept Lett , 26 ( 7 ): 532 - 541 . https://doi.org/10.2174/0929866526666190405124353 https://doi.org/10.2174/0929866526666190405124353 Effect of temperature and pH on the secondary structure and denaturation process of jumbo squid hepatopancreas cathepsin D Fuchimoto DI , Mizukoshi A , Schultz RM , et al. , 2001 . Posttranscriptional regulation of cyclin A1 and cyclin A2 during mouse oocyte meiotic maturation and preimplantation development . Biol Reprod , 65 ( 4 ): 986 - 993 . https://doi.org/10.1095/biolreprod65.4.986 https://doi.org/10.1095/biolreprod65.4.986 Posttranscriptional regulation of cyclin A1 and cyclin A2 during mouse oocyte meiotic maturation and preimplantation development Ge LH , Xu YS , Xia WS , et al. , 2018 . Synergistic action of cathepsin B, L, D and calpain in disassembly and degradation of myofibrillar protein of grass carp . Food Res Int , 109 : 481 - 488 . https://doi.org/10.1016/j.foodres.2018.04.067 https://doi.org/10.1016/j.foodres.2018.04.067 Synergistic action of cathepsin B, L, D and calpain in disassembly and degradation of myofibrillar protein of grass carp Gougeon A , 1996 . Regulation of ovarian follicular development in primates: facts and hypotheses . Endocr Rev , 17 ( 2 ): 121 - 155 . https://doi.org/10.1210/edrv-17-2-121 https://doi.org/10.1210/edrv-17-2-121 Regulation of ovarian follicular development in primates: facts and hypotheses Greenfeld CR , Babus JK , Furth PA , et al. , 2007 . BAX is involved in regulating follicular growth, but is dispensable for follicle atresia in adult mouse ovaries . Reproduction , 133 ( 1 ): 107 - 116 . https://doi.org/10.1530/REP-06-0144 https://doi.org/10.1530/REP-06-0144 BAX is involved in regulating follicular growth, but is dispensable for follicle atresia in adult mouse ovaries Gwon SH , Kim HK , Baek HJ , et al. , 2017 . Cathepsin B & D and the survival of early embryos in red spotted grouper, Ephinephelus akaara . Dev Reprod , 21 ( 4 ): 457 - 466 . https://doi.org/10.12717/DR.2017.21.4.457 https://doi.org/10.12717/DR.2017.21.4.457 Cathepsin B & D and the survival of early embryos in red spotted grouper, Ephinephelusakaara Han P , Xin HY , Peng JY , et al. , 2017 . Identification and expression of X-linked inhibitor of apoptosis protein during follicular development in goat ovary . Theriogenology , 98 : 30 - 35 . https://doi.org/10.1016/j.theriogenology.2017.04.048 https://doi.org/10.1016/j.theriogenology.2017.04.048 Identification and expression of X-linked inhibitor of apoptosis protein during follicular development in goat ovary Hsueh AJW , Kawamura K , Cheng Y , et al. , 2015 . Intraovarian control of early folliculogenesis . Endocr Rev , 36 ( 1 ): 1 - 24 . https://doi.org/10.1210/er.2014-1020 https://doi.org/10.1210/er.2014-1020 Intraovarian control of early folliculogenesis Hu HY , Jia Q , Zhou B , et al. , 2020 . Comparative analysis of the ovarian transcriptome reveals novel insights into fertility differences in Large White sows . Genes Genomics , 42 ( 7 ): 715 - 725 . https://doi.org/10.1007/s13258-020-00926-8 https://doi.org/10.1007/s13258-020-00926-8 Comparative analysis of the ovarian transcriptome reveals novel insights into fertility differences in Large White sows Kranc W , Budna J , Kahan R , et al. , 2017 . Molecular basis of growth, proliferation, and differentiation of mammalian follicular granulosa cells . J Biol Regul Homeost Agents , 31 ( 1 ): 1 - 8 . Molecular basis of growth, proliferation, and differentiation of mammalian follicular granulosa cells Leyria J , Fruttero LL , Ligabue-Braun R , et al. , 2018 . DmCatD, a cathepsin D-like peptidase of the hematophagous insect Dipetalogaster maxima (Hemiptera: Reduviidae): purification, bioinformatic analyses and the significance of its interaction with lipophorin in the internalization by developing oocytes . J Insect Physiol , 105 : 28 - 39 . https://doi.org/10.1016/j.jinsphys.2018.01.002 https://doi.org/10.1016/j.jinsphys.2018.01.002 DmCatD, a cathepsin D-like peptidase of the hematophagous insect Dipetalogastermaxima (Hemiptera: Reduviidae): purification, bioinformatic analyses and the significance of its interaction with lipophorin in the internalization by developing oocytes Ling YH , Xiang H , Li YS , et al. , 2014 . Exploring differentially expressed genes in the ovaries of uniparous and multiparous goats using the RNA-Seq (Quantification) method . Gene , 550 ( 1 ): 148 - 153 . https://doi.org/10.1016/j.gene.2014.08.008 https://doi.org/10.1016/j.gene.2014.08.008 Exploring differentially expressed genes in the ovaries of uniparous and multiparous goats using the RNA-Seq (Quantification) method Liu JH , Yang M , Jing L , et al. , 2018 . Silica nanoparticle exposure inducing granulosa cell apoptosis and follicular atresia in female Balb/c mice . Environ Sci Pollut Res , 25 ( 4 ): 3423 - 3434 . https://doi.org/10.1007/s11356-017-0724-5 https://doi.org/10.1007/s11356-017-0724-5 Silica nanoparticle exposure inducing granulosa cell apoptosis and follicular atresia in female Balb/c mice Liu JL , Yang L , Tian HY , et al. , 2016 . Cathepsin D is involved in the oxygen and glucose deprivation/reperfusion-induced apoptosis of astrocytes . Int J Mol Med , 38 ( 4 ): 1257 - 1263 . https://doi.org/10.3892/ijmm.2016.2709 https://doi.org/10.3892/ijmm.2016.2709 Cathepsin D is involved in the oxygen and glucose deprivation/reperfusion-induced apoptosis of astrocytes Lv FH , Agha S , Kantanen J , et al. , 2014 . Adaptations to climate-mediated selective pressures in sheep . Mol Biol Evol , 31 ( 12 ): 3324 - 3343 . https://doi.org/10.1093/molbev/msu264 https://doi.org/10.1093/molbev/msu264 Adaptations to climate-mediated selective pressures in sheep Mahdavi M , Nanekarani S , Hosseini SD , 2014 . Mutation in BMPR-IB gene is associated with litter size in Iranian Kalehkoohi sheep . Anim Reprod Sci , 147 ( 3-4 ): 93 - 98 . https://doi.org/10.1016/j.anireprosci.2014.04.003 https://doi.org/10.1016/j.anireprosci.2014.04.003 Mutation in BMPR-IB gene is associated with litter size in Iranian Kalehkoohi sheep Martins FS , Saraiva MVA , Magalhães-Padilha DM , et al. , 2014 . Presence of growth hormone receptor (GH-R) mRNA and protein in goat ovarian follicles and improvement of in vitro preantral follicle survival and development with GH . Theriogenology , 82 ( 1 ): 27 - 35 . https://doi.org/10.1016/j.theriogenology.2014.02.019 https://doi.org/10.1016/j.theriogenology.2014.02.019 Presence of growth hormone receptor (GH-R) mRNA and protein in goat ovarian follicles and improvement of invitro preantral follicle survival and development with GH McGee EA , Hsueh AJW , 2000 . Initial and cyclic recruitment of ovarian follicles . Endocr Rev , 21 ( 2 ): 200 - 214 . https://doi.org/10.1210/edrv.21.2.0394 https://doi.org/10.1210/edrv.21.2.0394 Initial and cyclic recruitment of ovarian follicles Minarowska A , Minarowski L , Karwowska A , et al. , 2007 . Regulatory role of cathepsin D in apoptosis . Folia Histochem Cytobiol , 45 ( 3 ): 159 - 163 . Regulatory role of cathepsin D in apoptosis Morais RDVS , Thomé RG , Lemos FS , et al. , 2012 . Autophagy and apoptosis interplay during follicular atresia in fish ovary: a morphological and immunocytochemical study . Cell Tissue Res , 347 ( 2 ): 467 - 478 . https://doi.org/10.1007/s00441-012-1327-6 https://doi.org/10.1007/s00441-012-1327-6 Autophagy and apoptosis interplay during follicular atresia in fish ovary: a morphological and immunocytochemical study Pan ZX , Zhang JB , Lin F , et al. , 2012 . Expression profiles of key candidate genes involved in steroidogenesis during follicular atresia in the pig ovary . Mol Biol Rep , 39 ( 12 ): 10823 - 10832 . https://doi.org/10.1007/s11033-012-1976-2 https://doi.org/10.1007/s11033-012-1976-2 Expression profiles of key candidate genes involved in steroidogenesis during follicular atresia in the pig ovary Peng JY , Gao KX , Xin HY , et al. , 2016 . Molecular cloning, expression analysis, and function of decorin in goat ovarian granulosa cells . Domest Anim Endocrinol , 57 : 108 - 116 . https://doi.org/10.1016/j.domaniend.2016.05.006 https://doi.org/10.1016/j.domaniend.2016.05.006 Molecular cloning, expression analysis, and function of decorin in goat ovarian granulosa cells Pradeep PK , Li XL , Peegel H , et al. , 2002 . Dihydrotestosterone inhibits granulosa cell proliferation by decreasing the cyclin D2 mRNA expression and cell cycle arrest at G1 phase . Endocrinology , 143 ( 8 ): 2930 - 2935 . https://doi.org/10.1210/endo.143.8.8961 https://doi.org/10.1210/endo.143.8.8961 Dihydrotestosterone inhibits granulosa cell proliferation by decreasing the cyclin D2 mRNA expression and cell cycle arrest at G1 phase Sales CF , Melo RMC , Pinheiro APB , et al. , 2019 . Autophagy and Cathepsin D mediated apoptosis contributing to ovarian follicular atresia in the Nile tilapia . Mol Reprod Dev , 86 ( 11 ): 1592 - 1602 . https://doi.org/10.1002/mrd.23245 https://doi.org/10.1002/mrd.23245 Autophagy and Cathepsin D mediated apoptosis contributing to ovarian follicular atresia in the Nile tilapia Soede NM , Langendijk P , Kemp B , 2011 . Reproductive cycles in pigs . Anim Reprod Sci , 124 ( 3-4 ): 251 - 258 . https://doi.org/10.1016/j.anireprosci.2011.02.025 https://doi.org/10.1016/j.anireprosci.2011.02.025 Reproductive cycles in pigs Strzalka W , Ziemienowicz A , 2011 . Proliferating cell nuclear antigen (PCNA): a key factor in DNA replication and cell cycle regulation . Ann Bot , 107 ( 7 ): 1127 - 1140 . https://doi.org/10.1093/aob/mcq243 https://doi.org/10.1093/aob/mcq243 Proliferating cell nuclear antigen (PCNA): a key factor in DNA replication and cell cycle regulation Sui MX , Wang HH , Wang ZW , 2015 . Molecular cloning, polymorphisms, and expression analysis of the RERG gene in indigenous Chinese goats . Genet Mol Res , 14 ( 4 ): 14936 - 14946 . https://doi.org/10.4238/2015.November.24.1 https://doi.org/10.4238/2015.November.24.1 Molecular cloning, polymorphisms, and expression analysis of the RERG gene in indigenous Chinese goats Tang HL , Huang XJ , Wang J , et al. , 2019 . circKIF4A acts as a prognostic factor and mediator to regulate the progression of triple-negative breast cancer . Mol Cancer , 18 : 23 . https://doi.org/10.1186/s12943-019-0946-x https://doi.org/10.1186/s12943-019-0946-x circKIF4A acts as a prognostic factor and mediator to regulate the progression of triple-negative breast cancer Tulone C , Uchiyama Y , Novelli M , et al. , 2007 . Haematopoietic development and immunological function in the absence of cathepsin D . BMC Immunol , 8 : 22 . https://doi.org/10.1186/1471-2172-8-22 https://doi.org/10.1186/1471-2172-8-22 Haematopoietic development and immunological function in the absence of cathepsin D Xia Y , Wang Q , He XD , et al. , 2020 . Cloning and expression analysis of the follicle-stimulating hormone receptor (FSHR) gene in the reproductive axis of female yaks ( Bos grunniens ) . Domest Anim Endocrinol , 70 : 106383 . https://doi.org/10.1016/j.domaniend.2019.07.011 https://doi.org/10.1016/j.domaniend.2019.07.011 Cloning and expression analysis of the follicle-stimulating hormone receptor (FSHR) gene in the reproductive axis of female yaks (Bosgrunniens) Yang ZW , Yang XY , Liu GB , et al. , 2020 . Polymorphisms in BMPR-IB gene and their association with litter size trait in Chinese Hu sheep . Anim Biotechnol, ahead-of-print. https://doi.org/10.1080/10495398.2020.1789158 https://doi.org/10.1080/10495398.2020.1789158 Polymorphisms in BMPR-IB gene and their association with litter size trait in Chinese Hu sheep Zhang H , Wu ZM , Yang YP , et al. , 2019 . Catalpol ameliorates LPS-induced endometritis by inhibiting inflammation and TLR4/NF-κB signaling . J Zhejiang Univ-Sci B (Biomed & Biotechnol) , 20 ( 10 ): 816 - 827 . https://doi.org/10.1631/jzus.B1900071 https://doi.org/10.1631/jzus.B1900071 Catalpol ameliorates LPS-induced endometritis by inhibiting inflammation and TLR4/NF-κB signaling Zhang P , Wang J , Lang H , et al. , 2018 . Knockdown of CREB1 promotes apoptosis and decreases estradiol synthesis in mouse granulosa cells . Biomed Pharmacother , 105 : 1141 - 1146 . https://doi.org/10.1016/j.biopha.2018.06.101 https://doi.org/10.1016/j.biopha.2018.06.101 Knockdown of CREB1 promotes apoptosis and decreases estradiol synthesis in mouse granulosa cells Zhang T , Rawson DM , Tosti L , et al. , 2008 . Cathepsin activities and membrane integrity of zebrafish ( Danio rerio ) oocytes after freezing to -196 ℃ using controlled slow cooling . Cryobiology , 56 ( 2 ): 138 - 143 . https://doi.org/10.1016/j.cryobiol.2008.01.002 https://doi.org/10.1016/j.cryobiol.2008.01.002 Cathepsin activities and membrane integrity of zebrafish (Daniorerio) oocytes after freezing to -196 ℃ using controlled slow cooling Zhang XD , Huang L , Wu T , et al. , 2015 . Transcriptomic analysis of ovaries from pigs with high and low litter size . PLoS ONE , 10 ( 10 ): e0139514 . https://doi.org/10.1371/journal.pone.0139514 https://doi.org/10.1371/journal.pone.0139514 Transcriptomic analysis of ovaries from pigs with high and low litter size Zhinan ZHOU and Xiang CHEN conceived and designed the experimental program. Zhinan ZHOU, Min ZHU, Weiwei WANG, and Zheng AO wrote and edited the manuscript. Zhinan ZHOU, Wen TANG, and Lei HONG performed the experiments. Zhinan ZHOU and Jiafu ZHAO revised the manuscript. All authors have read and approved the final manuscript, and therefore, have full access to all data in this study and take responsibility for the integrity and security of the data. This study was supported by the National Key R&D Program of China (No. 2017YFD0501904), the National Natural Science Foundation of China (No. 31760652), and the Major Projects of Science and Technology in Guizhou Province (QK-major projects 2016-3002), China. 10.1631/jzus.B2100366 1673-1581（2021）11-0893-13 This study was supported by the National Key R&D Program of China (No. 2017YFD0501904), the National Natural Science Foundation of China (No. 31760652), and the Major Projects of Science and Technology in Guizhou Province (QK-major projects 2016-3002), China. 2021-04-18 2021-05-27 2021-11-15 2022-05-09 2021-10-28 Copyright © Zhejiang University and Springer-Verlag GmbH Germany, part of Springer Nature 2021 2021 浙江大学出版社 Zhinan ZHOU, Xiang CHEN, Min ZHU, Weiwei WANG, Zheng AO, Jiafu ZHAO, Wen TANG, and Lei HONG declare that they have no conflict of interest. The animal handling procedures in this study adhered to the Animal Welfare Guidelines of the Animal Protection and Use Committee of Guizhou University, Guiyang, China (approval number: EGZU-2017T010). All animal handling procedures were performed to ensure minimal suffering. Compliance with ethics guidelines 浙江大学学报（英文版）（B辑：生物医学和生物技术） 11 22