==== Front Sci Rep Sci Rep Scientific Reports 2045-2322 Nature Publishing Group UK London 79013 10.1038/s41598-020-79013-7 Article Samul-tang ameliorates oocyte damage due to cyclophosphamide-induced chronic ovarian dysfunction in mice Kim Jihyun You Sooseong ethink33@kiom.re.kr grid.418980.c0000 0000 8749 5149Clinical Medicine Division, Korea Institute of Oriental Medicine, 1672 Yuseongdae-ro, Yuseong-gu, Daejeon, 34054 Republic of Korea 14 12 2020 14 12 2020 2020 10 2192522 6 2020 3 12 2020 © The Author(s) 2020Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.Samul-tang (SM), a traditional herbal medicine, has been used to treat menstrual irregularities and infertility in women. However, the cellular and molecular mechanisms underlying the effects of SM remain elusive. We investigated the potential protective effect of SM against chronic ovarian dysfunction and used bioinformatics analysis to identify its underlying mechanism in a mouse model of cyclophosphamide (CP)-induced diminished ovarian reserve. Female C57BL/6 mice were intraperitoneally injected with CP three times a week, followed by oral administration of distilled water (CP group) or SM (CP + SM group) for 4 weeks. Four weeks later, the effect of SM was assessed by ovarian tissue histological analysis, steroid hormone measurement, oocyte quality, and mRNA and microRNA microarray analysis in the ovaries. Although SM administration did not prevent CP-induced follicle loss in mice, the quality of oocytes was better in CP + SM mice than in CP mice. Gene expression analysis revealed that the expression of fertilisation- and ovarian follicle development-related genes was altered by CP treatment but normalized after SM administration. Further bioinformatics analysis showed possible interactions between differentially expressed mRNAs and microRNAs. Therefore, we demonstrated the protective effects of SM on ovarian function and oocyte maturation against CP-induced damage via multiple epigenetic mechanisms. Subject terms Developmental biologyGeneticsMolecular medicinehttp://dx.doi.org/10.13039/501100003718Korea Institute of Oriental MedicineKSN2013240You Sooseong issue-copyright-statement© The Author(s) 2020 ==== Body Introduction Extensive research is currently being conducted to develop alternative protective options against various devastating impacts of chemotherapy, surgery, and endocrine-disrupting chemicals on reproductive health1. In particular, herbal formulas have been proven to prevent and treat many diseases by promoting recovery and ameliorating impaired functions2. Samul-tang (SM), consisting of equal proportions of four herbs, Paeonia lactiflora, Ligusticum striatum, Rehmannia glutinosa, and Angelica gigas, is a well-known mixed herbal medicine. In in vivo and in vitro human and rodent models, SM has been reported to exert biological protective effects against haematological disorders3, oxidative stress4, apoptotic cell death5, and inflammation6. Furthermore, SM is commonly prescribed to women with gynaecological disorders, such as irregular menstruation and postmenopausal syndrome, which often occur in individuals with diminished ovarian reserve (DOR)7. DOR is defined as a decrease in the quantity and quality of mature oocytes, leading to infertility8. Women with DOR often have low blood flow to the ovaries, which has a negative consequence for competence9. The beneficial effect of SM on blood circulation has been reported10, even though the molecular mechanisms of action remain elusive. The aetiology of DOR includes previous ovarian surgery, exposure to toxic chemicals, chemotherapy, and radiation11. Repeated monthly chemotherapy administration results in folliculogenesis failure and severely impairs ovarian function12. In particular, cyclophosphamide (CP) exposure is associated with a long-term risk of premature ovarian failure and early menopause13,14. Further, CP exposure induces oxidative stress in spermatogenetic cells, granulosa cells and oocytes, thereby inducing organ toxicity15. Herbal medicine was found to attenuate CP-induced organ toxicity via anti-oxidative and anti-apoptosis pathways16. With the growing numbers of young female cancer survivors and the aim of fertility preservation, the maintenance of reproductive potential after cancer therapy has become a paramount concern. Although herbal medicines are traditionally known to improve fertility-related parameters, scientific evidence on their use in adjuvant therapy is limited. Therefore, in this study, we investigated the potential protective effect of SM and its underlying molecular mechanism in mice with a CP-induced chronic impairment of ovarian function. Results SM reduced CP-induced impairment of hormonal regulation The mice were monitored weekly throughout the study (Fig. 1A). The body weights of the CP-treated mice rapidly decreased during the week of CP administration and recovered only slightly after the cessation of CP exposure (Fig. 1B). SM administration resulted in a progressive, although not statistically significant, increase in body weight.Figure 1 Body weight changes and serum hormone levels in mice at 4 weeks after Samul-tang (SM) administration. (A) Mice were intraperitoneally injected with saline (n = 4, Con group) or cyclophosphamide (CP) on days 1, 3 and 5 in the first week; mice injected with CP were then orally administered distilled water (n = 5, CP group) or 5 g/kg SM (n = 5, CP + SM group) five times a week for 4 weeks. Four weeks after SM administration, the mice were weighed, and hormonal assessment was performed. (B) Body weight changes. (C) Serum levels of follicle-stimulating hormone (FSH), luteinising hormone (LH), and anti-Müllerian hormone (AMH), and the LH/FSH ratio. Statistical analysis was performed using the Student’s t-test. Arrows and arrowheads indicate the days of CP injection and SM administration, respectively. NS not significant. We next measured the serum follicle-stimulating hormone (FSH), luteinising hormone (LH), and anti-Müllerian hormone (AMH) levels in the Con, CP, and CP + SM mice. The levels of FSH and LH did not differ between CP and CP + SM mice. The CP mice had low LH levels, resulting in a decreased LH/FSH ratio (Fig. 1C). However, SM administration restored the LH/FSH ratio to a value comparable with that in the control (Con) mice. The serum AMH levels in the CP + SM mice were not markedly different from those in the CP and Con mice. SM did not prevent CP-induced follicle loss To investigate whether SM exerts a protective effect on ovarian follicle growth under conditions of CP exposure, we performed a histological analysis of the ovaries excised from the mice at 4 weeks after SM administration (Fig. 2A). The cross-sectional equivalent diameters of the ovaries were similar among the three mouse groups (Fig. 2B). Histological changes induced by CP exposure were observed in the ovarian tissues (Fig. 2C). Compared to the Con mice, the CP and CP + SM mice showed lower total numbers of follicles at different stages (Fig. 2D). Furthermore, the latter two groups showed significantly higher proportions of growing (secondary and preovulatory) follicles (P < 0.05; Fig. 2E). The follicular growth was not fount to differ between the CP and CP + SM mice. Thus, SM administration did not prevent the reduction of the primordial follicle pool and the CP-induced burnout phenomenon17.Figure 2 Histological analysis of ovarian follicles in mice at 4 weeks after Samul-tang (SM) administration. (A) Mice were intraperitoneally injected with saline (n = 4, Con group) or cyclophosphamide (CP) on days 1, 3, and 5 in the first week; mice injected with CP were then orally administered distilled water (n = 4, CP group) or 5 g/kg SM (n = 5, CP + SM group) five times a week for 4 weeks. Four weeks after SM administration, ovaries from all three groups were assessed histologically. (B) Cross-sectional equivalent diameter of the ovaries. (C) Representative histological images of ovaries from mice in the three experimental groups. (D) Numbers of ovarian follicles in different stages and the total number of ovarian follicles; and (E) proportions of primordial, primary, and growing (secondary and preovulatory) follicles in the Con, CP, and CP + SM mice. Statistical analysis was performed using the Student’s t-test. Arrows and arrowheads indicate the days of CP injection and SM administration, respectively. NS not significant. SM ameliorated CP-induced impairment of oocyte quality To investigate whether SM exerts a protective effect on oogenesis under conditions of CP exposure, the mice were hormonally superovulated (Fig. 3A). Oocytes were collected from the oviducts at 18 h post-hCG administration, and their quantity and quality were assessed. Oocytes from CP mice showed ZP defect as loosely compacted and enlarged (Fig. 3B). The numbers of total oocytes and mature metaphase II (MII) oocytes retrieved from the CP mice were significantly lower than those retrieved from the Con mice (P < 0.05; Fig. 3C). The number of MII oocytes, with normal chromosomes and well-organised spindle alignments, was significantly higher in the CP + SM mice than in the CP mice displaying chromosomal abnormalities and spindle misalignments (Fig. 3C).Figure 3 Quality and quantity of mouse oocytes retrieved at 4 weeks after Samul-tang (SM) administration. (A) Mice were injected with saline (n = 5, Con group) or cyclophosphamide (CP) on days 1, 3, and 5 in the first week; mice injected with CP were then orally administered distilled water (n = 5, CP group) or 5 g/kg SM (n = 5, CP + SM group) five times a week for 4 weeks. Four weeks after SM administration, the mice were superovulated via intraperitoneal injections of 5 IU of pregnant mare serum gonadotropin (PMSG) and 5 IU of human chorionic gonadotropin (hCG) at 48 h after PMSG injection. (B) Oocytes retrieved at 18 h after hCG injection from the Con, CP, and CP + SM mice. (C) Numbers of oocytes, MII oocytes, and MII oocytes with normal chromosomal and spindle alignment, retrieved from the three mouse groups. (D) Areas of oocytes, the ooplasm, and perivitelline space (PVS) in the three mouse groups. Statistical analysis was performed using the Student’s t-test. Scale bars indicate 100 µm. Arrows and arrowheads indicate the days of CP injection and SM administration, respectively. NS not significant. Moreover, CP treatment led to an increase in oocyte size and an enlargement of the perivitelline space (PVS). These effects were reversed upon SM administration so that the CP + SM mice had PVS and ooplasm areas similar to those in the Con mice (Fig. 3D). These results indicate that SM can protect the quality of superovulated oocytes even at 4 weeks after the completion of SM administration. SM reversed CP-induced changes in mRNA expression in mouse ovaries We next performed a microarray analysis to compare the mRNA expression patterns in the ovulated ovaries from the Con, CP, and CP + SM mice (Fig. 4A). Hierarchical clustering analysis showed marked differences among the three mouse groups (Fig. 4B), with 33,793 genes displaying significantly different expression levels in at least one comparison set.Figure 4 Hierarchical clustering and analysis of differentially expressed mRNAs. (A) Microarray was performed to compare the gene expression in ovulated ovaries of Con (n = 5), CP (n = 5), and CP + SM mice (n = 5). (B) Hierarchical clustering among the mRNA expression profiles. (C) Venn diagram showing 89 differentially expressed mRNAs in the three groups, with a fold-change > 1.5 and P < 0.05. Con saline-treated control mice; CP cyclophosphamide-treated mice; CP + SM CP-treated mice orally administered Samul-tang. Arrows and arrowheads indicate the day of CP injection and SM administration, respectively. Of the 33,793 genes identified, 89 showed significantly altered expression, with fold changes > 1.5 (P < 0.05). The numbers of differentially expressed genes (DEGs) between the three mouse group pairs are presented in a Venn diagram (Fig. 4C). Comparison of the Con and CP mice data revealed 13 DEGs, including the Ereg and Rhox8 genes, which are involved in ovarian folliculogenesis and oocyte ovulation (Fig. 4C, Supplementary Fig. S1A). Compared to the Con mice, both CP and CP + SM mice showed significant increases in the expression of three genes, Hspa1b, Hspa1a, and Mterf1b (Fig. 4C, Supplementary Fig. S1B). Gene expression was validated by qPCR using fluorescent probe-based Taqman assays (Supplementary Fig. S1). We found that expression patterns of Rhox8 and Mterf1b were different compared to microarray data. These genes are involved in the regulation of mitotic spindle assembly and protein refolding. These results, indicating CP-induced differential gene expression, might explain the presence of abnormal oocytes, with disarranged spindle alignment, resulting in fertilisation failure. Compared to the CP mice, the CP + SM mice showed 44 significant DEGs, of which 31 (70.5%) were upregulated and 13 (29.5%) were downregulated (Supplementary Tables S1 and S2). Interestingly, nine genes overlapped between the CP vs Con and CP + SM vs CP pairs (Fig. 4C). Of these, the zona pellucida glycoprotein 2 and 3 genes (Zp2 and Zp3), which play an important role in fertilisation, were downregulated in the CP mice, but their expression levels were restored by SM administration (Fig. 5A). Thirty DEGs which were differentially expressed only in CP + SM vs CP mice included NLR family pyrin domain containing 5 (Nlrp5), bone morphogenetic protein 15 (Bmp15), 2′-5′ oligoadenylate synthetase 1D (Oas1d), and oocyte-specific homeobox 1 and 5 (Obox1 and Obox5) (Fig. 4C). These genes are known to be involved in fertilisation, folliculogenesis, and oocyte maturation (Fig. 5A–C). The expression of these genes was validated by qPCR with fluorescent probe-based Taqman assays (Fig. 6A–C). We found that the expression pattern of Obox5 was different compared to microarray data. CP treatment significantly altered the expression of these five genes in the ovaries, whereas their expression levels in the ovaries of the CP + SM mice were comparable to those in the ovaries of the Con mice. These findings suggest that SM administration ameliorated the CP-induced impairment of ovarian function via genetic regulation. The five DEGs that overlapped between the CP + SM vs Con and CP + SM vs CP pairs were associated with immunoglobulin formation and immune responses (Fig. 4C).Figure 5 Differentially expressed mRNAs involved in ovarian function. Microarray fluorescence intensities indicate the expression levels of (A) Zp2, Zp3, and Nlrp5 genes, involved in fertilisation, (B) Bmp15 and Oas1d genes, involved in ovarian follicle development, and (C) Obox1 and Obox5 genes, involved in oocyte maturation, in ovaries of the Con, CP, and CP + SM mice. Statistical analysis was performed using the Student’s t-test. Con saline-treated control mice; CP cyclophosphamide-treated mice; CP + SM CP-treated mice orally administered Samul-tang. Figure 6 Validation of the expression of differentially expressed mRNAs involved in ovarian function. qPCR was performed to validate the expression of (A) Zp2, Zp3, and Nlrp5 genes involved in fertilisation, (B) Bmp15 and Oas1d genes, involved in ovarian follicle development, and(C) Obox1 and Obox5 genes, involved in oocyte maturation, in the ovaries of the Con, CP, and CP + SM mice. Statistical analysis was performed using the Student’s t-test. Con: saline-treated control mice; CP cyclophosphamide-treated mice; CP + SM CP-treated mice orally administered Samul-tang. Gene ontology (GO) analysis of the CP and CP + SM mouse datasets revealed that the DEGs were enriched in biological processes such as ovarian follicle development and fertilisation (Fig. 7A), in molecular functions such as structural constituents of the egg coat and acrosin binding (Fig. 7B), and in cellular components such as the egg coat and secretory and circulating IgA and IgM immunoglobulin complexes (Fig. 7C).Figure 7 Gene Ontology (GO) and KEGG analysis of differentially expressed mRNAs. Targets were subjected to pathway analysis using KEGG and subsequently classified based on their enrichment in (A) biological process, (B) molecular function, and (C) cellular component GO terms. Asterisks (*) indicate the most important biological processes altered by SM administration compared to those in CP mice. Con saline-treated control mice; CP cyclophosphamide-treated mice; CP + SM CP-treated mice orally administered Samul-tang. SM reversed CP-induced changes in the expression of microRNAs in mouse ovaries A small number of microRNAs were differentially expressed in the three mouse groups. A comprehensive microarray tool showed the presence of 998.6 ± 37.4, 1,014.8 ± 26.8, and 1,052.2 ± 17.3 microRNAs in the Con, CP, and CP + SM mice, respectively (Supplementary Fig. S2A). Volcano plots showing differential microRNA expression among the three mouse groups are presented in Supplementary Fig. S2B. Hierarchical clustering analysis and a volume plot revealed that the expression of six microRNAs was significantly different (fold-change > 1.5 and P < 0.05) between the CP + SM and CP mice (Fig. 8A,B). Interestingly, CP treatment decreased the expression of miR-365-1-5p and miR-574-5p and increased that of miR-3472 and miR-667-3p, whereas SM administration restored the expression levels of these microRNAs to those in the Con mice (Fig. 8C). Compared to the Con and CP mice, the CP + SM mice showed significant differential expression of two microRNAs, miR-200b-3p and miR-665-3p (Fig. 8C). Collectively, these findings indicate that SM administration may ameliorate the CP-induced impairment of ovarian function via epigenetic regulation.Figure 8 Hierarchical clustering and analysis of differentially expressed microRNAs. Microarray analysis was performed to compare the microRNA expression in ovulated ovaries of the CP and CP + SM mice. (A) Hierarchical clustering among the microRNA expression profiles. (B) Volume plot of log2 fold changes between the CP and CP + SM mice, showing six significantly differentially expressed microRNAs with fold changes > 1.5 and P < 0.05. (C) Microarray fluorescence intensities indicate the expression levels of miR-365-1-5p, miR-3472, miR-200b-3p, miR-665-3p, miR-667-3p, and miR-574-5p in the Con, CP, and CP + SM mice. CP cyclophosphamide-treated mice; CP + SM CP-treated mice orally administered Samul-tang. Integrative analysis of differentially expressed mRNAs and microRNAs Further enrichment analysis using the TargetScan webtool identified possible direct interactions between differentially expressed mRNAs and microRNAs (Fig. 9). We found that the Bmp15 mRNA harboured binding sites for miR-200b-3p and miR-665-3p and that the Oas1 family transcript harboured a binding site for miR-667-3p, of which expressions were changed by SM administration after CP exposure.Figure 9 Binding sites between differentially expressed mRNAs and microRNAs. Binding sites in (A) Bmp15 for miR-200b-3p and miR-665-3p and (B) Oas1 for miR-667-3p, as predicted by TargetScan analysis. Discussion Traditional medicine-based treatment strategies for infertility vary between patients. SM has been traditionally used as an adjuvant herb in Korea to treat menstrual irregularities and infertility in women18. However, it remains unclear whether herbal treatments may improve clinical outcomes in patients with infertility. Therefore, in this study, we investigated the effects of SM as a fertoprotective adjuvant in a CP-induced DOR mouse model. SM administration restored the LH/FSH ratio, which is known to be associated with follicular depletion19,20. These findings suggested that SM administration initiated the hormonal feedback regulatory system. Interestingly, the oocyte quality in mice significantly improved after SM administration without the protection of impaired physiological conditions and the loss of ovarian follicles. CP treatment is known to induce oxidative stress in granulosa cells and oocytes during folliculogenesis and ovulation. Furthermore, CP treatment induces ovulation of premature oocytes, thereby resulting in poor embryo quality21. SM has reportedly shown radical-scavenging activities against the hydroxyl radical and nitric oxide22, which are known to deteriorate the MII oocyte spindle and cause poor oocyte quality23. Our results demonstrated that SM exerted protective effects on oocytes with respect to chromosomal normality and spindle alignment, supposedly via its antioxidant activity. SM administration inhibited the CP-induced alterations in the spatiotemporal dynamics of oocyte maturation and aided in complete oocyte maturation. To elucidate the underlying mechanism of the protective effects of SM, we performed microarray analysis to investigate the mRNA and microRNA expression in ovaries. SM administration reversed CP-induced alterations in the expression of genes related to the fertilisation potential, such as those involved in sperm-egg recognition, negative regulation of fertilisation, and regulation of spindle assembly. The DEGs identified between the CP and CP + SM groups were found to be involved in ovarian follicle development, binding of sperm to ZP, and microtubule nucleation. Of these, CP induced significant decreases in the expression of Zp2, Zp3, and Nlrp5, which was restored to normal levels after SM administration. Reduced mRNA levels of the Zp2 and Zp3 genes result in a low fertilisation rate and IVF failure in humans and mice24,25 Moreover, Nlrp5 expression has been reported to decrease with maternal age, and Nlrp5-deficient human oocytes show mitochondrial depletion, resulting in the arrest of embryonic maturation26,27. However, the microarray and qPCR data of Rhox8, Mterf1b, and Obox5 were in disagreement. These variations may be attributed to the filtering condition of microarray data for measures of quality including fold-change > 1.5 and P < 0.05 as previously described28. Further investigation is necessary for identifying the genetic regulation of SM. Importantly, the expression of oocyte-specific genes was restored in CP-treated mice after SM administration. BMP15, OAS1D, OBOX1, and OBOX5 are well-known oocyte-specific growth factors involved in ovarian folliculogenesis and oocyte maturation29,30, and downregulation of their encoding genes has been suggested to compromise fertility in women undergoing cancer treatment. Bmp15 is known to regulate ovulation by preventing the overactivation of primordial follicles, and the expression of this gene is used as a marker of oocyte quality31. The reversal of the CP-induced decrease in the Bmp15 expression by SM administration improved the quality of ovulated oocytes. Although little is known about the function of Obox family genes, their transcription is induced during zygotic gene activation, suggesting functional compensation between Obox family members32,33. Interestingly, among the differentially expressed microRNAs identified in our study, three microRNAs (miR-200b-3p, miR-665-3p, and miR-667-3p) have binding sites on either Bmp15 or Oas1 mRNAs. These results indicate that changes in the expression of microRNAs and their targets may be involved in the regulation of different fertility-related processes, from ovulation to fertilisation. The expression of miR-365-1-5p, miR-3472, miR-574-5p, and miR-667-3p was restored in CP-treated mice at 4 weeks after SM administration. Little is known about the roles of miR-365-1-5p, miR-3472, and miR-667-3p. Meanwhile, one of the known functions of miR-574-5p is tumour suppression via the inhibition of MMP3 expression, which is important for tissue remodelling and extracellular matrix maintenance during the follicular and luteal phases34. Thus, SM-mediated restoration of the miR-574-5p expression may play a role in the protective effects of SM on folliculogenesis and ovulation in CP-treated mice. These changes in microRNA expression may be involved in the complex epigenetic regulation of ovarian function. Our study did not perform the functional analysis of differentially expressed mRNAs and microRNAs after SM administration in chronically impaired ovarian function. However, this is the first report on the effect of SM on CP-induced ovarian dysfunction using integrative microarray analysis. It provides scientific evidence on its use in adjuvant therapy. The comprehensive microarray analysis performed in this study showed that the adverse effect of CP on ovarian function was sustained for 8 weeks after CP exposure, which is sufficient time to replace CP-damaged primordial and growing follicles surrounding oocytes. This time is considered equivalent to more than 6 months in humans35,36. Women subjected to anticancer therapy are advised to avoid conception for 6 months to up to 2 years to prevent adverse effects on the foetus37. For application in DOR patients after CP treatment, it should have a protective effect against chronic ovarian dysfunction. Collectively, the differential mRNA and microRNA expression analysis and subsequent integrative analysis performed in this study highlight the protective effect of SM on ovarian function and oocyte maturation under chronically impaired conditions. This protective effect remained for at least 4 weeks after the termination of SM administration and led to an increase in the number of mature oocytes. Therefore, SM may be used to increase the oocyte yield in IVF patients with DOR after CP treatment, potentially improving the chances of IVF success. Methods Mice All the experiments and analyses were conducted in accordance with the relevant guidelines and regulations. Experimental animal protocols were approved by the Institutional Animal Care and Use Committee at the Korea Institute of Oriental Medicine, Daejeon, Korea (approval number 19–019). Female 8-week-old C57BL/6 mice (Nara Biotech, Pyeongtaek, Korea) were housed under specific pathogen-free conditions. The mice were administered an intraperitoneal injection of saline (control or Con group) or an equal volume of 75 mg/kg CP (Sigma–Aldrich, St. Louis, MO, USA) at days 1, 3, and 5 for one week. To investigate the effect of SM, powdered SM (Hanpoong, Iksan, Korea) was dissolved in distilled water, and the CP-treated mice were orally administered either distilled water (CP group) or 5 g/kg SM (CP + SM group) five times a week for 4 weeks using feeding needles. Mice from the three groups were separated within the same cage to synchronise their hormonal cycles38. All mice were sacrificed 4 weeks after the completion of SM administration. This post-treatment interval was selected because it provided sufficient time for the development of preantral follicles from the newly recruited primordial follicles35. Ovaries were removed and immediately placed into 4% paraformaldehyde (Biosesang, Seongnam, Korea) or liquid nitrogen for histological or microarray analysis, respectively. Enzyme-linked immunosorbent assay (ELISA) for hormonal assessment At 4 weeks after SM administration, blood samples were collected from the mice, and sera were separated and stored at − 80 °C until analysis. The serum concentrations of FSH, LH, and AMH were measured using hormone-specific ELISA kits from Cusabio Biotech Co. (Wuhan, China), Enzo Life Sciences (Farmingdale, NY, USA), and MyBioSource (San Diego, CA, USA), respectively, according to standard protocols and the manufacturers’ instructions. For FSH, both intra- and inter-assay coefficients of variation (CVs) were < 15%, with a sensitivity of 2.5 mIU/mL. For LH, the intra- and inter-assay CVs were 7% and 15%, respectively, and the functional sensitivity was 5.2 mIU/mL. For AMH, the inter-assay CV was < 10%, with a sensitivity of 0.19 ng/mL. Histological assessment of ovarian follicles At 4 weeks after SM administration, both ovaries were serially sectioned to obtain 5-μm-thick tissue sections, which were stained with haematoxylin and eosin. The cross-sectional equivalent diameters of ovaries were measured using NIS-Elements BR version 4.60.00 (Nikon Instruments, Tokyo, Japan) according to the manufacturer’s guidelines. Primordial, primary, secondary, and preovulatory follicles, with visible oocytes, were counted in every fifth stained section to avoid counting the same follicle twice. The follicular stages were classified as previously described39: primordial follicles, with a single flat layer of granulosa cells surrounding the oocyte; primary follicles, with a single cuboidal granulosa cell layer; secondary follicles, with at least two granulosa cell layers and a theca layer; and preovulatory follicles, with a complete antrum and theca layer. Assessment of oocyte quantity and quality At 4 weeks after SM administration, the mice were superovulated via an intraperitoneal injection of 5 IU of pregnant mare serum gonadotropin (Prospec, Rehovot, Israel), followed by an injection of 5 IU of human chorionic gonadotropin (hCG; Prospec) 48 h later. Oocytes were collected at 18 h post-hCG injection into the preincubated human tubal fluid medium (Irvine Scientific, CA, USA). Oocytes were fixed with 4% paraformaldehyde (Biosesang) and permeabilised with 0.5% Triton X-100 (Sigma–Aldrich) for 10 min, followed by blocking with PBS containing 3% bovine serum albumin (GenDEPOT, TX, USA) and incubation with a rabbit anti-α-tubulin antibody (1:200; Cell Signaling Technology, MA, USA). Oocytes were mounted on slides using the VECTASHIELD antifade mounting medium with 4′,6-diamidino-2-phenylindole (Vector Laboratories, Peterborough, UK) to visualise the chromosomes and were observed under a fluorescent microscope (BX51; Olympus, Tokyo, Japan). Oocytes with well-organised bipolar spindles and tightly aligned chromosomes at the metaphase were scored as normal. The oocyte quality was evaluated by measuring morphometric parameters, including the areas of the whole oocyte, ooplasm, and PVS, using a computed image analysis system (Nikon Instruments). Microarray analysis of mRNA and microRNA expression Four weeks after SM administration, the ovaries were collected from the mice post-ovulation, and total RNA was extracted using the RNeasy mini kit (Qiagen, Venlo, Netherlands) according to the manufacturer’s instructions. The purity and integrity of the extracted RNA were evaluated using a NanoDrop ND-1000 UV–Vis spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). All the samples showed a high purity [optical density (OD)260/OD280 > 1.80] and integrity (RNA integrity number > 7.0). Hybridisation with GeneChip Mouse Gene 2.0 ST arrays (Affymetrix) was controlled using the GeneChip Command Console software. MicroRNA expression was analysed using the Affymetrix GeneChip microRNA array 4.0, based on MirBase release 20. The Affymetrix Expression Console 1.4 software was used for basic data extraction (CEL files) and quality control metrics. A fold-change value > 1.5 and a P-value < 0.05 were used as thresholds to identify DEGs. Integrative analysis of mRNA and microRNA expression profiles Gene enrichment and functional annotation analysis for a significant probe list were performed using GO (http://geneontology.org) and the Kyoto Encyclopedia of Genes and Genomes (KEGG; http://kegg.jp) to identify the potential functions of DEGs in the biological process, molecular function, and cellular component categories. Data analyses and visualisation of DEGs were performed using R 3.5.1 (www.r-project.org). An integrative analysis was performed using the TargetScan webtool 7.2 (www.targetscan.org, accessed March 2018) to find binding sites for critical microRNAs and DEGs. Validation of selected differentially expressed gene expression in the ovaries To confirm mRNA microarray results, validation was performed on significant genes of interest using real-time quantitative PCR (qPCR). The synthesis of complementary DNAs from extracted total RNA was performed using iScript cDNA Synthesis kit (Bio-Rad Laboratories, Hercules, CA) according to the manufacturer’s instructions. The qPCR was performed in a final reaction volume of 20 µL using QuantStudio 6 Flex Real-time PCR System with fluorescent probe-based Taqman assays according to the manufacturer’s instructions (Thermo Fisher Scientific). The cycle threshold was normalised and compared using Gapdh as the internal standards. Statistical analysis The statistical significance of differences between two groups was determined by the Student’s t-test using GraphPad Prism version 8.4.0 (GraphPad Software, Inc., La Jolla, CA, USA). Differences with P < 0.05 were considered statistically significant. Supplementary Information Supplementary Information. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Supplementary Information The online version contains supplementary material available at 10.1038/s41598-020-79013-7. Acknowledgements This work was supported by a project grant from the Korea Institute of Oriental Medicine (KSN2013240). Author contributions J.K.: conception, experiments, writing, and revision of the manuscript. S.Y.: experiments, revision of the manuscript, and correspondence. Data availability The datasets generated and/or analysed during the current study are available from the corresponding author on reasonable request. Competing interests The authors declare no competing interests. ==== Refs References 1. Spears N Ovarian damage from chemotherapy and current approaches to its protection Hum. Reprod. Update. 2019 25 673 693 10.1093/humupd/dmz027 31600388 2. Sanders K Moran Z Shi Z Paul R Greenlee H Natural products for cancer prevention: clinical update 2016 Semin. Oncol. Nurs. 2016 32 215 240 10.1016/j.soncn.2016.06.001 27539278 3. Choi SC Yae HK Effects of Samultang on blood components, coagulation and platelet aggregation in rats J. Cancer Prev. 2012 17 110 115 4. Lee I-C Kim B-H Kim M-K Biological activity and single oral dose toxicity of Dansam-samultang Lab. Anim. Res. 2010 26 385 391 10.5625/lar.2010.26.4.385 5. Lee GW Kim MS Water extract of samultang reduces apoptotic cell death by H2O2-induced oxidative injury in SK-N-MC cells Korean J. Physiol. Pharmacol. 2009 13 139 145 10.4196/kjpp.2009.13.3.139 19885029 6. Choi ES Samul-Tang regulates cell cycle and migration of vascular smooth muscle cells against TNF-α stimulation Evid. Based Complement. Alternat. Med. 2018 2018 102497 10.1155/2018/1024974 7. Huang G-C Tsai Y-Z Lee C-J Chang H-Y Wang C-C Elucidation of the effects of Sisi-Wu Tang on menstrual disorder patterns through activation of aromatase and antioxidation Evid. Based Complement. Alternat. Med. 2019 2019 4761651 10.1155/2019/4761651 30949219 8. Rasool S Shah D Fertility with early reduction of ovarian reserve: the last straw that breaks the Camel's back Fertil. Res. Pract. 2017 3 15 10.1186/s40738-017-0041-1 29046817 9. Yuan H Efficiency of DHEA on diminished ovarian reserve patients undergoing IVF-ET based on ultrasonography Int. J. Clin. Exp. Med. 2016 9 14337 14344 10. Sun J To unveil the molecular mechanisms of Qi and blood through systems biology-based investigation into Si-Jun-Zi-Tang and Si-Wu-Tang formulae Sci. Rep. 2016 2016 34328 10.1038/srep34328 11. Georgescu ES Goldberg JM Plessis SSD Agarwal A Present and future fertility preservation strategies for female cancer patients Obstet. Gynecol. Surv. 2008 63 725 732 10.1097/OGX.0b013e318186aaea 18928577 12. Morgan S Anderson RA Gourley C Wallace WH Spears N How do chemotherapeutic agents damage the ovary? Hum. Reprod. Update. 2012 18 525 535 10.1093/humupd/dms022 22647504 13. Mahajan N Fertility preservation in female cancer patients: an overview J. Hum. Reprod. 2015 8 3 13 10.4103/0974-1208.153119 14. Marder W Fisseha S Ganser MA Somers EC Ovarian damage during chemotherapy in autoimmune diseases: broad health implications beyond fertility Clin. Med. Insights. Reprod. Health. 2012 2012 9 18 23970822 15. Cengiz M Cyclophosphamide induced oxidative stress, lipid per oxidation, apoptosis and histopathological changes in rats: Protective role of boron J. Trace. Elem. Med. Biol. 2020 62 126574 10.1016/j.jtemb.2020.126574 32516632 16. Cengiz M Kutlu HM Peker Cengiz B Ayhancı A Escin attenuates oxidative damage, apoptosis and lipid peroxidation in a model of cyclophosphamide-induced liver damage Drug. Chem. Toxicol. 2020 10.1080/01480545.2020.1810262 33280440 17. Kalich-Philosoph L Cyclophosphamide triggers follicle activation and "burnout"; AS101 prevents follicle loss and preserves fertility Sci. Transl. Med. 2013 5 185ra162 10.1126/scitranslmed.3005402 18. Fang Z Evaluating the pharmacological mechanism of Chinese medicine Si-Wu-Tang through multi-level data integration PLoS ONE 2013 8 e72334 10.1371/journal.pone.0072334 24223693 19. Lee JE Yoon SH Kim HO Min EG Correlation between the serum luteinizing hormone to folliclestimulating hormone ratio and the anti-Mullerian hormone levels in normo-ovulatory women J. Korean Med. Sci. 2015 30 296 300 10.3346/jkms.2015.30.3.296 25729253 20. Bedoschi G Navarro PA Oktay K Chemotherapy-induced damage to ovary: mechanisms and clinical impact Future Oncol. 2016 12 2333 2344 10.2217/fon-2016-0176 27402553 21. Jeelani R Cyclophosphamide and acrolein induced oxidative stress leading to deterioration of metaphase II mouse oocyte quality Free. Radic. Biol. Med. 2017 110 11 18 10.1016/j.freeradbiomed.2017.05.006 28499912 22. Nam MN Lee AY Sin SM Goo Y-M Cho EJ Protective effect of Samultang and its four herbal plants against reactive oxygen species in vitro and cellular system Korean J. Agri. Sci. 2019 46 427 437 23. Shaeib F Banerjee J Maitra D Diamond MP Abu-Soud HM Impact of hydrogen peroxide-driven Fenton reaction on mouse oocyte quality Free. Radic. Biol. Med. 2012 58 154 159 10.1016/j.freeradbiomed.2012.12.007 23261938 24. Rankin TL Defective zonae pellucidae in Zp2-null mice disrupt folliculogenesis, fertility and development Development 2001 128 1119 1126 11245577 25. Liu C Targeted disruption of the mZP3 gene results in production of eggs lacking a zona pellucida and infertility in female mice Proc. Natl. Acad. Sci. U. S. A. 1996 93 5431 5436 10.1073/pnas.93.11.5431 8643592 26. Sharov AA Effects of aging and calorie restriction on the global gene expression profiles of mouse testis and ovary BMC Biol. 2008 6 24 10.1186/1741-7007-6-24 18522719 27. Fernandes R NLRP5 mediates mitochondrial function in mouse oocytes and embryos Biol. Reprod. 2012 86 138 149 10.1095/biolreprod.111.093583 22357545 28. Morey JS Fau RJ Dolah V Microarray validation: factors influencing correlation between oligonucleotide microarrays and real-time PCR Biol. Proced. Online 2006 8 175 193 10.1251/bpo126 17242735 29. Urzúa U The ovarian transcriptome of reproductively aged multiparous mice: candidate genes for ovarian cancer protection Biomolecules 2020 10 113 10.3390/biom10010113 30. Racki WJ Richter JD CPEB controls oocyte growth and follicle development in the mouse Development 2006 133 4527 4537 10.1242/dev.02651 17050619 31. Sanfins A Rodrigues P Albertini DF GDF-9 and BMP-15 direct the follicle symphony J. Assist. Reprod. Genet 2018 35 1741 1750 10.1007/s10815-018-1268-4 30039232 32. Stein P Medvedev S Schultz RM OBOX proteins are recruited during oocyte maturation and are essential for early development in mouse Biol. Reprod. 2012 87 210 10.1093/biolreprod/87.s1.210 33. Cheng W-C Hsieh-Li HM Yeh Y-J Li H Mice lacking the Obox6 homeobox gene undergo normal early embryonic development and are fertile Dev. Dyn. 2007 236 2636 2642 10.1002/dvdy.21261 17676645 34. Zheng J Zhou Y Li X-J Hu J-M MiR-574-3p exerts as a tumor suppressor in ovarian cancer through inhibiting MMP3 expression Eur. Rev. Med. Pharmacol. Sci. 2019 23 6839 6848 31486483 35. McGee EA Hsueh AJ Initial and cyclic recruitment of ovarian follicles Endocr. Rev. 2000 21 200 214 10.1210/edrv.21.2.0394 10782364 36. Smith MF Ricke WA Bakke LJ Dow MPD Smith GW Ovarian tissue remodeling: role of matrix metalloproteinases and their inhibitors Mol. Cell. Endocrinol. 2002 191 45 56 10.1016/S0303-7207(02)00054-0 12044918 37. Lawrenz B Pregnancy after breast cancer: case report and review of the literature Arch. Gynecol. Obstet. 2011 283 837 843 10.1007/s00404-010-1829-y 21221981 38. McClintock MK Social control of the ovarian cycle and the function of estrous synchrony Am. Zool. 1981 21 243 256 10.1093/icb/21.1.243 39. Pedersen T Peters H Proposal for a classification of oocytes and follicles in the mouse ovary J. Reprod. Fertil. 1968 17 555 557 10.1530/jrf.0.0170555 5715685