ORIGINAL RESEARCH article

Front. Endocrinol., 19 October 2022

Sec. Experimental Endocrinology

Volume 13 - 2022 | https://doi.org/10.3389/fendo.2022.994642

Yin Huo Tang, a traditional Chinese herbal formula, relives ovariectomy and empty bottle stimulation-induced menopause-like symptoms in mice

  • 1. Department of Traditional Chinese Medicine, Peking University Third Hospital, Beijing, China

  • 2. Department of Neuroscience, Institute for Chinese Medicine, Heilongjiang University of Chinese Medicine, Harbin, China

  • 3. Department of Integrated Traditional Chinese and Western Medicine, College of Medicine, Yangzhou University, Yangzhou, China

  • 4. Department of Pharmacy, Wuxi Higher Health Vocational Technology School, Wuxi, China

  • 5. National Institute on Drug Dependence and Beijing Key Laboratory of Drug Dependence, Peking University, Beijing, China

  • 6. Research Institute of Medicine and Pharmacy, Qiqihar Medical University, Qiqihar, China

Abstract

Background:

Yin Huo Tang (YHT), a traditional Chinese herbal formula, is effectively used for the clinical treatment of menopause-like symptoms in China. This study aimed to investigate its efficacy on menopause-like symptoms in mice using behavioral tests and histopathological assessment, and to determine its possible mechanism of action based on network pharmacology.

Methods:

Liquid chromatography-mass spectrometry (LC-MS) technology was used to identify the potential active ingredients of YHT. In mice, menopause-like symptoms were induced by combination of bilateral ovariectomy and empty bottle stimulation. The mice were then treated with the YHT aqueous extract for three weeks. Behavior, sleep state, body weight, organ index, and histomorphology were analyzed separately. Additionally, network pharmacology and molecular docking were used to predict the mechanisms underlying the action of YHT. Finally, serum estradiol was quantified to preliminarily verify the results of network pharmacology.

Results:

YHT not only improved the behavior of mice (attack and explore behavior reduced; modify behavior increased) but also ameliorated the sleep state (sleep time increased and incubation time reduced). YHT reduced body weight, increased uterine weight, and improved the histomorphology of some organs. Network pharmacology and molecular docking analyses revealed that the estrogen signaling pathway might play a key role in attenuating menopause-like symptoms. Furthermore, YHT treatment reversed the reduction in serum estradiol levels.

Conclusions:

YHT alleviates menopause-like symptoms in a mouse model, providing a rationale for using it as a potential therapeutic strategy.

Introduction

Menopause is the moment women stop menstruating, defined as 12 consecutive months following the last menstrual period. The average age at which women experience menopause is 45–55 years (1). Generally, menopause occurs with physical symptoms (hot flashes, sweating, and poor circulation), emotional symptoms (depressive mood, lack of energy, concentration difficulties, sleep disorders, and emotional instability), and organic disorders (local urogenital atrophy, skin atrophy, and development of osteoporosis) with ample variations (2). The pathogenesis of the menopausal syndrome is mainly associated with a decrease in estrogen levels (3). Hormone replacement therapy (HRT) is an effective treatment method for the menopausal syndrome, which relieves menopausal discomfort and improves the quality of life; however, the use of HRT is restricted by its potential adverse effects (4, 5). Thus, finding an effective treatment method with fewer side effects is important for women with menopausal syndrome.

In China, traditional Chinese medicine (TCM) has become an important therapy for women with menopausal syndrome. Yin Huo Tang (YHT), a classic Chinese medicine formula, was originally recorded in Bian Zheng Lu, which was an ancient Chinese medical book written by Chen Shiduo during the Qing Dynasty. This formula consists of five commonly used herbs: 1) Rehmannia glutinosa (Gaetn.) Libosch. ex Fisch. et Mey. (Shu Dihuang); 2) Morinda officinalis How. (Ba Jitian); 3) Ophiopogon japonicus (L.f.) Ker-Gawl. (Mai Dong); 4) Poria cocos(Schw.)Wolf (Fu Ling); 5) Schisandra chinensis (Wu Weizi). YHT has been clinically used to treat patients with Yin deficiency syndrome such as female menopause syndrome (68). In an animal study, YHT demonstrated protective effects against postmenopausal depression and the SIRT1/PGC-1α pathway may be involved in this action (9). The therapeutic potential of YHT requires further confirmation, whereas its possible mechanism of action requires further exploration in animal models.

In this study, we investigated the pharmacological effects of YHT in a mouse model of menopausal syndrome. We assessed the behavior, sleep state, body weight, and organ index, and analyzed the histomorphology and serum estradiol levels in mice. Moreover, we employed network pharmacology and molecular docking to reveal the possible mechanisms underlying the protective effects of YHT. The results of this study provide new insights into this field and provide a scientific basis for the application of YHT in menopausal syndrome.

Methods

Animals

In total, 155 female ICR mice were employed at the age of five weeks and the study was performed in a randomized and blinded manner. Animals were housed at an ambient temperature of 22°C with a relative humidity of 60% and a 12 h light-dark cycle. They were provided water ad libitum, except during the testing period. Food was provided ad libitum throughout the experiments. This study was conducted in strict accordance with the local institutional guidelines for animal experimentation. The experimental protocol was approved by the Animal Care and Ethics Committee of the Heilongjiang University of Chinese Medicine. During experiments, all efforts were made to humanize the experiments for animals.

YHT herbal materials

Dried YHT, consisting of five herbal components (Table 1) was purchased from Beijing Tongrentang Pharmaceutical (Harbin, China). We slightly adjusted the doses of the herbs used in this study, as needed. The crude individual herbal components of YHT were cut into small pieces and then mixed at a ratio of 3:2:2:1:1. They were then suspended in distilled water for two hours to soften the raw materials and facilitate the extraction of water-soluble ingredients in subsequent steps. The herbs were boiled for one hour and passed through filter paper. The process lasted three times. The recommended daily intake of YHT is 45 g/day. Converting the clinically equivalent dose for a 70 kg adult to that for a 20 g mouse by body surface area ratio, the reference dose for mice was 5.85 g/kg/day (low dose in the study) (10). The filtrate was condensed to the required concentrations and stored at 4°C for later use.

Table 1

Chinese nameBotanical nameCommon nameGenusFamilyWeight (g)Part used
Shu DihuangRehmannia glutinosa (Gaert.) Libosch. ex Fisch. et Mey.dried rehmannia rootRehmanniaScrophulariaceae15Root
Ba JitianMorinda officinalis How.morinda officinalisMorindaRubiaceae10Root
Mai DongOphiopogon japonicus (L.f.) Ker-Gawl.radix ophiopogonisOphiopogonLiliaceae10Root
Fu LingPoria cocos(Schw.) WolftuckahoePoriaPolyporaceae5Sclerotia
Wu WeiziSchisandra chinensisfruit of Chinese magnoliavineSchisandraMagnoliaceae5Fruit

Contents of Yin Huo Tang.

LC-MS analysis of YHT

Liquid chromatography-mass spectrometry (LC-MS) analysis was performed using an ACQUITY UHPLC system (Waters Corporation, Milford, USA) coupled with an AB SCIEX Triple TOF 5600 System (AB SCIEX, Framingham, USA). The concentration of YHT analyzed was 0.15 g/mL (standard concentration: 45 g YHT/300 mL water). An ACQUITY UPLC HSS T3 column (2.1×100 mm, 1.8 μm) was employed in both the positive and negative modes. The binary gradient elution system consisted of (A) water (containing 0.1% formic acid, v/v) and (B) acetonitrile (containing 0.1% formic acid, v/v), and separation was achieved using the following gradient: 0 min, 5% B; 2 min, 5% B; 4 min, 25% B; 8 min, 50% B; 10 min, 80% B; 14 min, 100% B; 15 min, 100% B; 15.1 min, 5%; 16 min, 5%B. The flow rate and the injection volume were 0.35 mL/min and 2 μL, respectively. The column temperature was maintained at 45°C, and the samples were kept at 4°C. Data acquisition was performed in the full scan mode (m/z range from 100 to 1000) combined with the IDA mode. The MS parameters were as follows: ESI ion source temperature: 550°C (+) and 550°C (−); ion spray voltage: 5500 V (+) and 4500 V (−); curtain gas: 35 PSI; nebulizer gas: 55 PSI; auxiliary gas: 55 PSI; declustering potential: 80V (+) and 80V (−); collision energy: 10 eV (+) and −10 eV (−); interface heater temperature: 550°C (+) and 550°C (−). The LC-MS data were processed using Progenesis QI V2.3 (Nonlinear, Dynamics, Newcastle, UK) for baseline filtering, peak identification, integration, retention time correction, peak alignment, normalization, and analysis.

Drug administration and experimental protocol

The mice were acclimatized for one week before the experiment. Subsequently, they were randomly assigned to the following eight groups: control group, sham group (sham operation), sham + empty bottle stimulation group, model group (combination of ovariectomy and empty bottle stimulation), positive control group (Progynova, 0.13 mg/kg, Bayer HealthCare), YHT low-dose group (YHT-L, 5.85 g/kg), YHT middle-dose group (YHT-M, 11.7 g/kg), and YHT high-dose group (YHT-H, 23.4 g/kg) (10). All drugs were administrated daily via the intragastric route. The mice in the control group did not undergo any intervention throughout the experiments.

After the habituation period, mice in the sham + empty bottle stimulation, model, Progynova, YHT-L, YHT-M, and YHT-H groups were trained to drink water at a specific time for 14 days. Furthermore, mice in these groups (except the sham + empty bottle stimulation group) underwent bilateral ovariectomy on the 11th day. From the 15th to the 35th day, the mice in these groups underwent empty bottle stimulation. On the 36th day, all the animals were sacrificed for further analysis.

Bilateral ovariectomy

This model is described in detail in our previously published study (11). Briefly, mice were anesthetized, and the operation area was shaved and disinfected using 75% ethanol. A small incision (1-1.2 cm) was then made along the ventral midline between the midpoint of the leading edge of the pubic bone and the navel. The uterus and the side of the uterine horn were found, and the ovary was identified along the uterine horn. Subsequently, the ovarian mesenteric blood vessels were ligated with silk threads, and the ovaries were removed with a pair of small forceps and scissors. Finally, the incision was closed with sutures.

Empty bottle stimulation

Empty bottle stimulation was then used to create a model of chronic emotional stress (11). For two weeks, mice were trained to drink water from 08:00 h to 08:10 h and 20:00 h to 20:10 h by allowing them access to water bottles only during these time periods. To induce emotional stress, mice were given empty water bottles irregularly during one of the two watering periods for three weeks after training. Mice in the control and sham groups were provided water ad libitum throughout the experiment.

Body weight and organ index evaluation

All mice were weighed once per week throughout the experiment. After three weeks of treatment, all mice were sacrificed, and the uteri, adrenal glands, spleens, thymi, and pituitary glands were carefully dissected, washed with cold sterile physiological saline, and weighed. Their weights relative to their final body weights were calculated as organ indices. After weighing, the uteri and vaginas were fixed with 10% neutral buffered formalin solution for subsequent histological assessment.

Histopathological assessment

Hematoxylin and eosin (H&E) staining was performed for the histological assessment. The uteri and vaginas were fixed with 10% neutral buffered formalin solution, embedded in paraffin wax, and stained with H&E in 5-µm sections using a standard protocol. All images were captured using a microscope (ECLIPSE 50i, Nikon, Japan). The white balance was adjusted for all images.

Pentobarbital-induced sleep test

The pentobarbital-induced sleep test was performed with a slight modification to a previous report (12). Briefly, mice were administered with pentobarbital sodium (50 mg/kg, i.p.) to induce sleep. The latency of the loss of the righting reflex after pentobarbital administration (“onset” time) and the total sleeping time (the time between the loss and the recovery of the righting reflex, “duration”) were determined for each mouse. The mouse was considered awake if it could right itself (return to the upright position). The time between the loss and recovery of the righting reflex was defined as sleep time. The entire process of the test was recorded using a digital video recorder for analysis.

Behavioral analysis

Behavioral assessments were performed as described in our previous study (11). The mice were tested in a cage with only an empty bottle. The behavior of each mouse was recorded using a digital video recorder for 10 min. The observed behaviors were divided into three patterns: attacking behavior (biting or attacking the empty bottle and cage), exploring behavior (moving around or frequently visiting the place of the bottle), and grooming behavior (combing fur or touching the face by themselves). Each behavioral pattern was counted 10 times during the 10-minutes testing (once per minute). Once an animal exhibited a particular pattern of the three patterns, the score of that pattern was recorded as 1; otherwise, it was recorded as 0. The score for each behavioral pattern ranged from 0 to 10. The analysis was performed by an independent investigator who was blinded to group assignments.

Enzyme linked immunosorbent assay

Serum estradiol levels were measured using a commercially available enzyme-linked immunosorbent assay (ELISA) kit. The ELISA kit for mouse estradiol was purchased from Nanjing Jiancheng Bioengineering Institute (catalog no. E-20380, Nanjing, China), and used following the manufacturer’s instructions.

Network pharmacology analysis

The target proteins of each ingredient in YHT were obtained from the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) (13). Menopause-related targets were collected from the GeneCards database (14). The targets of YHT and menopause were introduced into Venny 2.1.0 to identify overlapping targets. Afterward, the ingredients-targets-disease network and protein-protein interaction (PPI) network were constructed using the Cytoscape V3.9.0 software. Gene Ontology analysis (including biological processes, cellular components, and molecular functions) was performed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) and Bioinformatics online tool (www.bioinformatics.com.cn). The functions of common genes were analyzed using the Metascape database (15).

Molecular docking

We first obtained the three-dimensional structures of ingredients and target proteins from the TCMSP and Protein Data Bank (PDB) databases, respectively (16). Water molecules and ligands of the target proteins were removed using PyMOL software. Next, the target proteins were hydrogenated and converted to PDBQT format. AutoDock software was used to obtain molecular docking results. Finally, the PyMOL software was used to visualize the results.

Statistical analysis

Data were presented as means ± standard errors of the mean. The mean number of data in individual groups was compared using the one-way/two-way analysis of variance (ANOVA) followed by Tukey’s test post-hoc analysis. Statistical work was performed using SPSS version 17.0 and a p-value < 0.05 was considered statistically significant.

Results

Compounds characterization of YHT

As shown in Table 2, 93 compounds were characterized using LC-MS, including galactitol, betaine, glucosamine, and L-asparagine. The base peaks in the positive and negative ion modes from the YHT analysis are presented in Figure 1. Among the 93 compounds, 56 were detected in positive ion mode and 37 in negative ion mode. For further study, the TCMSP database was used to predict the corresponding potential targets of the 93 compounds.

Table 2

No.RT (min)Ion modeIdentificationScoreError (ppm)KEGG IDFormula
10.71negGalactitol47.3-3.60C01697C6H14O6
20.71posBetaine40.3-0.06C00719C5H11NO2
30.71posGlucosamine50.7-0.55C00329C6H13NO5
40.71posL-Asparagine53.5-0.11C00152C4H8N2O3
50.73posSucrose53.1-0.58C00089C12H22O11
60.75posCreatine38.80.31C00300C4H9N3O2
70.75posL-Proline57.41.61C00148C5H9NO2
80.78negMaleic acid37.4-7.12C01384C4H4O4
90.84posAdenine46.3-0.10C00147C5H5N5
100.87posCitric acid43.3-0.38C00158C6H8O7
110.87posGuanine54.5-0.23C00242C5H5N5O
121.05negUridine54.20.14C00299C9H12N2O6
131.06posL-Tyrosine48.70.38C00082C9H11NO3
141.16posL-Lysine38.7-0.21C00047C6H14N2O2
151.16negQuinic acid51.51.30C00296C7H12O6
161.24negGuanosine52.51.33C00387C10H13N5O5
171.27negAucubin51.71.54C09771C15H22O9
181.33posAdenosine57.20.32C00212C10H13N5O4
191.35negAscorbic acid47.6-3.83C01041C6H8O6
201.40posPhenylethylamine37.53.37C05332C8H11N
211.50negFumaric acid39.70.07C00122C4H4O4
221.55posHippuric acid44.50.59C01586C9H9NO3
231.61posHordenine44.50.29C06199C10H15NO
241.63negGallic acid56.5-3.81C01424C7H6O5
252.02negKojic acid46.2-2.71C14516C6H6O4
262.11posHydroquinone44.82.84C00530C6H6O2
272.11posL-Phenylalanine54.50.33C00079C9H11NO2
282.22negShikimic acid48.4-4.29C00493C7H10O5
292.24posBiotin43.31.20C00120C10H16N2O3S
302.25negCatalpol52.12.18C09773C15H22O10
312.31negBenzoic acid37.9-5.07C00539C7H6O2
322.68pos4-Methylcatechol38.92.43C06730C7H8O2
332.93negChlorogenic acid541.90C00852C16H18O9
342.95neg(-)-Epigallocatechin51.61.63C12136C15H14O7
353.05posPimelic acid39.10.02C02656C7H12O4
363.42posGentisic acid39.20.48C00628C7H6O4
373.48posGenistein36.51.16C06563C15H10O5
383.58negL-Tryptophan53.8-1.58C00078C11H12N2O2
393.62posEphedrine38.90.72C01575C10H15NO
403.62posDaidzin37.60.14C10216C21H20O9
413.72negAconitic acid46.2-3.57C00417C6H6O6
423.76posEpicatechin46.3-0.13C09727C15H14O6
433.83posPiperine36.71.33C03882C17H19NO3
443.87negProtocatechuic acid48.1-5.19C00230C7H6O4
453.89negtrans-Cinnamic acid38.7-2.99C10438C9H8O2
463.89posp-Anisic acid45.60.40C02519C8H8O3
473.95posCinnamic acid42.70.48C10438C9H8O2
483.99posCitraconic acid38.7-0.18C02226C5H6O4
494.02posIsoferulic acid54.10.25C10470C10H10O4
504.03negTrehalose38.8-0.13C01083C12H22O11
514.08posNaringin38.71.28C09789C27H32O14
524.11negGeniposidic acid55.41.34C11673C16H22O10
534.12pos4-Hydroxybenzoic acid39.31.39C00156C7H6O3
544.15negIsomangiferin45.12.84C16979C19H18O11
554.18negScopoletin50.5-2.84C01752C10H8O4
564.22posApigenin41.6-0.43C01477C15H10O5
574.26negCaffeic acid56.6-3.83C01481C9H8O4
584.33posPhenylacetaldehyde53.32.11C00601C8H8O
594.36negHarpagoside41.12.05C09783C24H30O11
604.37posConiferin39.1-1.04C00761C16H22O8
614.41posSebacic acid39.2-0.16C08277C10H18O4
624.48negSyringic acid56.2-2.52C10833C9H10O5
634.49pos3-O-Feruloylquinic acid50.4-0.20C02572C17H20O9
644.74negSuberic acid53.3-3.71C08278C8H14O4
654.82negMethyl jasmonate38.91.63C11512C13H20O3
664.82posVanillin54.2-0.44C00755C8H8O3
674.84posTyramine390.73C00483C8H11NO
684.94posQuercitrin57-0.31C01750C21H20O11
695.03negIsorhamnetin541.80C10084C16H12O7
705.04posNobiletin38.7-1.32C10112C21H22O8
715.06posPapaverine36.4-0.62C06533C20H21NO4
725.12negBiochanin A39.11.63C00814C16H12O5
735.18negAzelaic acid55.4-1.24C08261C9H16O4
745.40negCinnamaldehyde38.7-4.53C00903C9H8O
755.42posOnonin57.3-1.21C10509C22H22O9
765.42posHesperetin37.9-0.68C01709C16H14O6
775.42negSalicylic acid55-7.04C00805C7H6O3
785.82negFormononetin57.30.87C00858C16H12O4
795.88posKaempferol38.7-1.27C05903C15H10O6
805.91negQuercetin51.51.01C00389C15H10O7
815.94negVanillic acid52.1-4.42C06672C8H8O4
826.10posCatechol38.81.82C15571C6H6O2
836.14posProtopine38.8-1.20C05189C20H19NO5
846.41posGlutaric acid38.80.04C00489C5H8O4
856.54posSclareol37-0.78C09183C20H36O2
867.01posSinapic acid39.6-0.16C00482C11H12O5
877.56negScoparone46.51.18C09311C11H10O4
888.63posCurdione42.7-0.85C17493C15H24O2
899.77posGamma-Linolenic acid390.16C06426C18H30O2
909.85posAlpha-Linolenic acid39.3-0.24C06427C18H30O2
9110.17posEstradiol50.2-0.13C00951C18H24O2
9211.20posDibutyl phthalate44.3-0.17C14214C16H22O4
9311.93posOleic acid39.5-0.07C00712C18H34O2

Identified chemical compounds within YHT by LC-MS.

RT, retention time; pos, positive; neg, negative.

Figure 1

Effects of YHT on body weight and organ index in menopausal mice

The experimental workflow used in this study is shown in Figure 2A. The estrous cycle stages were identified based on the cell types present, as shown in Figure 2B. Mice in the control group presented different stages of the estrous cycle. Mice in the ovariectomy groups only presented the diestrus stage, which is characterized by a large number of leukocytes (shown in Figures 2B–D). Body weight and organ indices were measured to identify the effects of YHT on the mice. The results showed that the body weight at week five was significantly higher in the model group than in the control group (Figure 2C). Compared with the model group, the body weights at week four and week five were significantly lower in the progynova, YHT-M, and YHT-H groups (Figure 2C). No significant difference in body weight was observed between the model and YHT-L groups throughout the experiment. Additionally, the organ index of the uterus was lower in the model group than that in the control group (Figure 2D). The uterine index in the progynova, YHT-M, and YHT-H groups were higher than that in the model group (Figure 2D). The adrenal gland index was significantly lower in the model group than that in the control group (Figure 2E). Neither progynova nor YHT reversed the decrease in adrenal gland index in menopausal mice. Our results also showed that the organ indices of the spleen, thymus, and pituitary gland did not show any significant difference in any of the groups (Figures 2F–H).

Figure 2

Histological effects of YHT on uterus and vagina in menopausal mice

To identify the effects of YHT on organ structure, the histomorphology of the uterus and vagina was assessed using H&E staining (Figures 3A, B). In menopausal mice, the lining of the uterus (endometrium) became thinner and smoother, and the endometrial glands decreased significantly compared to the control mice (Figures 3C, D). The histomorphology of the uterus in the Progynova, YHT-M, and YHT-H groups was similar to that of the control group. However, in the YHT-L group, the histomorphology of the uterus did not change much compared to the model group (Figures 3C, D). Compared to the control group, the vaginal epithelium was atrophied in the model group (Figure 3E). In the Progynova and YHT-M groups, vaginal epithelium atrophy was reversed. However, no significant difference in the histomorphology of the vagina was observed between the YHT-L and YHT-H groups and the model group (Figure 3E).

Figure 3

Effects of YHT on anxiety-like behavior and sleeping behavior in menopausal mice

Mice in the model group exhibited significant attacking behaviors compared with mice in the control group (Figure 4A). However, compared with the model group, attacking behaviors in all other groups were significantly decreased. In the model group, the mice demonstrated more exploratory behaviors than that in the control group (Figure 4B). However, the other groups showed less exploratory behaviors compared to the model group. In addition, grooming behaviors in the model group were reduced compared with those in the control group (Figure 4C). However, grooming behaviors in the other groups were greater than those in the model group.

Figure 4

To examine the effects of YHT on sleep behavior, all the mice were injected with pentobarbital (50 mg/kg). The results demonstrated that the onset time in the model group was longer than that in the control group (Figure 4D). In the YHT-L and YHT-H groups, the onset time decreased significantly compared to that in the model group. However, no significant difference was found in the onset time between the Progynova and YHT-M groups and the model group (Figure 4D). Moreover, compared to the control group, the model group exhibited a significant decrease in duration (Figure 4E). The duration increased greatly in the YHT-L and YHT-M groups compared with that in the model group. However, the duration in the Progynova and YHT-H groups was not significantly different from that in the model group (Figure 4E).

Ingredients-targets-disease network construction and PPI analysis

Using the TCMSP database, we identified 433 potential targets of YHT. Next, we predicted a total of 820 menopause-related targets from the GeneCards database (Relevance score > 1.00). After overlap, 146 intersectional genes were identified as potential key targets of YHT in the treatment of menopause (Figure 5A). The ingredients-targets-disease network was then constructed by using Cytoscape V3.9.0 software (Figure 5B). This network includes 274 nodes and 1129 edges. YHT active ingredients are represented by orange nodes, menopause targets by blue nodes, and YHT itself by green nodes. The number of links connected to other nodes indicates the importance of a node in the network (17). Topological analysis of degree revealed that the top 10 active ingredients were: quercetin, genistein, kaempferol, apigenin, estradiol, isorhamnetin, formononetin, nobiletin, epicatechin, and maleic acid.

Figure 5

The PPI network included 146 nodes and 3101 edges, and the average number of neighbors was 42.5 (Figure 5C). Topological analysis revealed that the top 10 core targets in the network were: ALB, AKT1, IL6, INS, TP53, VEGFA, JUN, IL1B, ESR1, and CASP3. Network topology analysis of the PPI network is shown in Figure 5D. Core targets were obtained by filtering twice with scores for degree centrality (DC), betweenness centrality (BC), and closeness centrality (CC).

GO and pathways enrichment analysis

We obtained 1006 GO items from the DAVID, including 782 biological process (BP), 79 cellular component (CC), and 145 molecular function (MF) terms. The top 20 most enriched items are shown in a bar chart (Figures 6A–C). The BPs were mainly related to the response to hypoxia, positive regulation of cell proliferation, aging, negative regulation of the apoptotic process, and response to estradiol. CCs were mainly related to the extracellular space, extracellular region, caveola, macromolecular complex, and platelet alpha granule lumen. The MFs were mainly related to enzyme binding, protein binding, cytokine activity, and heme binding. Pathway enrichment analysis revealed that the most significantly enriched pathways were cancer, response to hormones, and AGE-RAGE signaling pathways (Figure 6D).

Figure 6

Molecular docking analysis

Molecular docking was used to evaluate the binding abilities of the core compounds and protein targets of YHT. The lower the binding affinity, the better the binding activity between the compound and protein (18). As shown in Figure 7, AKT1, ALB, EGFR, ESR1, and PTGS2 generally exhibited good docking performance with the core active compounds of YHT (apigenin, estradiol, formononetin, genistein, and quercetin). In particular, Estradiol showed a high binding affinity of −11.1 kcal/mol and −10.2 kcal/mol in ALB and EGFR, respectively. Quercetin showed a high binding affinity of −10.4 kcal/mol in AKT1. Interestingly, estradiol showed a high binding affinity of −9.1 kcal/mol in ESR1, which further confirmed the potential key role of estradiol-related pathways.

Figure 7

Experimental validation and potential pathways prediction

To validate the results of network pharmacology and molecular docking, we measured serum levels of estradiol. The estradiol concentration in the model group was lower than that in the control group (Figure 8A). However, estradiol concentrations in the Progynova, YHT-M, and YHT-H groups were significantly higher than those in the model group. No significant differences were observed between the YHT-L and model groups (Figure 8A). Finally, the KEGG PATHWAY database was used to predict the potential pathways involved in the effects of YHT. According to the above analysis, the estradiol signaling pathway is considered to be a critical pathway underlying the protective effects of YHT (Figure 8B). Multiple molecules, such as ESR1, MAPK1, and AKT1, may participate in the therapeutic mechanism of YHT.

Figure 8

Discussion

Several women experience a range of menopause-like symptoms related to changes in sex hormone levels during menopause (19). The most common symptoms include depressive disorders, sleep disturbances, sexual dysfunction, muscle pain, osteoporosis, and characteristic hot flashes (20). HRT is believed to be the most effective intervention for menopausal symptoms (21). However, recent evidence has shown that long-term HRT may have negative consequences on women’s health (22, 23). In this study, bilateral ovariectomy and empty bottle stimulation were combined to induce menopause-like symptoms in the mice. Our study provides evidence that YHT administration in animals partly relieved menopause-like symptoms and provides a rationale for using the YHT as a potential therapeutic strategy to ameliorate menopausal syndrome.

YHT has already been used clinically to treat menopausal symptoms by many TCM doctors, with good effects (6, 7). This ancient herbal formula was created by Chen Shiduo, an outstanding physician during the Qing dynasty in China. According to TCM theory, YHT was applied to treat the syndrome of hyperactivity due to Yin deficiency, and symptoms of this syndrome were similar to female menopausal symptoms. Therefore, it is theoretically feasible to treat the menopausal symptoms of YHT. In addition, progynova was selected as the positive control agent because it showed a potent effect in many menopausal animal studies (24, 25).

The ovariectomy-induced menopausal model has been widely used in many studies. Empty bottle stimulation was initially used to induce an emotional stress model (26). Animals that underwent empty bottle stimulation showed a state of irritability, which was similar to the clinical situation of patients with an anxiety disorder (26, 27). In our previous study, we combined ovariectomy and empty bottle stimulation to induce menopause-like symptoms in mice (28). This model expressed typical menopause-like symptoms such as reduced estrogen level and sleep quality, increased body weight, and aggressive behaviors.

Being overweight and obese causes not only the physical inconvenience but also secondary diseases such as diabetes, coronary heart disease, hypertension, and stroke (29). The risk of overweight or obesity induced by estrogen deficiency increases significantly in post-menopausal females (30). In our study, the body weight of menopausal mice increased significantly compared to that of control mice. However, YHT showed an obvious improvement in weight gain. Given that estrogen plays an important role in body weight regulation, this mechanism may be related to an increase in estrogen levels (31). Menopause results in loss of ovarian function and an atrophic state of genital organs, such as the uterus and vagina (32). The reduction in the uterine index due to ovariectomy was reversed by YHT. In addition, the adrenal gland, spleen, thymus, and pituitary gland indices did not change in this study. The morphological results were consistent with those of the uterine index. YHT exhibited good effects in inhibiting endometrial atrophy and vaginal epithelial atrophy induced by ovariectomy and empty bottle stimulation.

Sleep disturbances are common during menopause (33). Menopausal women frequently complain of having difficulty in sleep initiation or sleep maintenance (34). However, the mechanisms underlying sleep disturbances that arise during menopause remain unclear. In this study, YHT significantly shortened the process of falling asleep and increased the duration of continuous sleep. However, based on our results, progynova did not improve sleep disturbance induced by ovariectomy and empty bottle stimulation. This suggests that purely supplemented estradiol did not present sufficient efficacy to relieve sleep disturbances during menopause.

With the decline in ovarian steroid levels, the peri- and postmenopausal periods represent a window of vulnerability for the emergence of anxiety symptoms and disorders (35). Many women experience profound anxiety attacks that start occurring while experiencing menopausal symptoms. The connection between anxiety attacks and menopause remains clear. Mice showed significant anxiety-like behaviors after ovariectomy and empty-bottle stimulation, and anxiety was significantly attenuated by YHT. Current data indicate that YHT may be as effective as estradiol in relieving anxiety during menopause.

Estrogens are female sex hormones involved in many physiological processes, especially reproductive development and function. Estradiol, the main form of circulating estrogen, is secreted by ovaries (36). The decline of estradiol at menopause often leads to functional disorders that affect the quality of life (37). In our study, YHT reversed the decline in estradiol induced by menopausal syndrome. This may partly explain why YHT relieved the atrophy of sexual organs and anxiety.

We obtained 93 active ingredients of YHT using LC-MS analysis. Among them, quercetin, genistein, kaempferol, apigenin, and estradiol were considered the core ingredients. Quercetin has been reported to increase the antioxidant capacity of ovaries in menopausal rats and in vitro (38). Genistein exhibits estrogenic activity and improves uterine morphology in mice (39). Another study demonstrated that kaempferol might exert a protective effect against post-menopausal bone loss (40). Recently, a study revealed that apigenin alleviates xerostomia via estrogen receptor α-mediated upregulation of aquaporin 5 activation (41). PPI analysis showed that ALB, AKT1, IL6, and INS were the key targets of YHT in treating menopause-like symptoms in mice. ALB, AKT1, IL6, and INS indicate albumin, RAC-alpha serine/threonine-protein kinase, interleukin-6, and insulin, respectively. Early menopause is associated with higher risk of urinary albumin-creatinine ratio elevation in postmenopausal women (42). The expression of Akt1 dysregulated in menopausal rats (43). In addition, the increased deposition of visceral fat promotes the release of interleukin-6 in the postmenopause period (44). Increased inflammation of beta cells impairs insulin secretion in a post-menopausal diabetic rat model (45).

BP enrichment and pathway enrichment analysis showed that the estradiol signaling pathway might play a key role in the protective effects of YHT. The docking results also showed that estradiol and ESR1 had a good docking performance. The predicted results were verified by ELISA of the serum estradiol levels. We found that YHT reversed the decrease in serum estradiol levels in menopausal mice. Our results are similar to those of a previous study (46). In that study, YHT promoted the release of endogenous estrogen in the brain and improved depression-like behavior induced by chronic unpredictable mild stress in ovariectomized mice.

In this study, we explored the optimal dosage of YHT for attenuation of menopause-like symptoms in mice. Interestingly, we did not observe a distinct dose-dependent effect of YHT on physiological and psychological symptoms in mice with menopausal syndrome among the three experimental dosages used. However, we found that mice in the YHT middle-dose group showed more improvement than that in the low-dose group in some indicators (such as body weight, uterus index, endometrial thickness, gland number, the thickness of vaginal epithelium, and estradiol concentration). The effect of the middle dosage of YHT was found to be the best among the three doses tested, which was probably due to the toxicity of YHT at high dosages and suggests the importance of exploring the optimal dosage in the clinic setting (47).

Our study has some limitations that need to be discussed. First, this study investigated only the pharmacological properties of YHT in menopausal animals. Further research needed to investigate how YHT can improve menopausal syndrome. Second, network pharmacology technology only provides clues based on existing research results, which indicates that the innovativeness of this study is relatively limited. Omics technologies such as transcriptomics and proteomics are needed to explore the targets of YHT for treating menopause syndrome.

Conclusion

In conclusion, YHT improved the bilateral ovariectomy and empty bottle stimulation-induced menopause-like symptoms in mice. Our findings suggest that YHT significantly ameliorates both physiological and psychological symptoms in a novel mouse model through a multitarget and multi-pathway approach. Animal experiments have validated that the estradiol-mediated pathway may play a key role in this process. Considering its significant therapeutic potency, we concluded that YHT has the potential to improve menopause-like symptoms in middle-aged and older women. Further studies on the potential mechanisms and compositional analysis of YHT are needed in the future.

Funding

This work was supported by National Natural Science Foundation of China (No. 81873108, 82003975), Natural Science Foundation of Heilongjiang Province (No. H2015-068), Project of Qiqihar Academy of Medical Sciences (No. QMSI2020L-05), and Excellent Creative Talents Support Program of Heilongjiang University of Chinese Medicine (No. 2018RCQ08).

Publisher’s note

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

Statements

Data availability statement

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

Ethics statement

The animal study was reviewed and approved by Committee on the Ethics of Heilongjiang University of Chinese Medicine.

Author contributions

YY: data collection, data analysis, manuscript writing. BZ: data collection, data analysis. YL: data collection. H-DX: writing - review and editing. X-ML: writing - review and editing. S-MH: supervision. RW: conceptualization, supervision. DL: project administration. All authors contributed to the article and approved the submitted version.

Conflict of interest

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

Supplementary material

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

References

  • 1

    DavisSRBaberRJ. Treating menopause - MHT and beyond. Nat Rev Endocrinol (2022) 18:490502. doi: 10.1038/s41574-022-00685-4

  • 2

    PatelBSDW. Menopause review: Emerging treatments for menopausal symptoms. Best Pract Res Clin Obstet Gynaecol (2022) 81:134–44. doi: 10.1016/j.bpobgyn.2021.10.010

  • 3

    SantoroNRoecaCPetersBANeal-PerryG. The menopause transition: Signs, symptoms, and management options. J Clin Endocrinol Metab (2021) 106:115. doi: 10.1210/clinem/dgaa764

  • 4

    FloresVAPalLMansonJE. Hormone therapy in menopause: Concepts, controversies, and approach to treatment. Endocr Rev (2021) 42:720–52. doi: 10.1210/endrev/bnab011

  • 5

    RozenbergSDi PietrantonioVVandrommeJGillesC. Menopausal hormone therapy and breast cancer risk. Best Pract Res Clin Endocrinol Metab (2021) 35:101577. doi: 10.1016/j.beem.2021.101577

  • 6

    ChenPWangMWuJ. Li′s YinHuo decoction in treating 58 cases of insomnia during the menopause period. West J Tradit Chin Med (2014) 27:57–8. doi: 10.3969/j.issn.1004-6852.2014.08.022

  • 7

    ChenPWangMHeF. LiShi YinHuo tang combined with classical prescription in treating 108 female patients with menopause insomnia. West J Tradit Chin Med (2016) 29:66–7. doi: 10.3969/j.issn.1004-6852.2016.12.022

  • 8

    WangMZhouTHuK. Application of yinhuo decoction in coping with side effects of radiotherapy and chemotherapy of malignant tumor. Acta Chin Med (2019) 34:1156–8. doi: 10.16368/j.issn.1674-8999.2019.06.275

  • 9

    XuHXingSLeiXYiJLiuSDuYet al. A famous Chinese medicine formula: Yinhuo decoction antagonizes the damage of corticosterone to PC12 cells and improves depression by regulating the SIRT1/PGC-1alpha pathway. BioMed Res Int (2022) 2022:3714857. doi: 10.1155/2022/3714857

  • 10

    LiuDLinSLiYZhouTHuKLiQ. Network pharmacology and experimental verification to explore the potential mechanism of yin-Huo-Tang for lung adenocarcinoma recurrence. Drug Des Devel Ther (2022) 16:375–95. doi: 10.2147/DDDT.S343149

  • 11

    YeYLiuCLiuXHuangS. Ovariectomy changes the response to antidepressant drugs in tail suspension test in mice. Gynecol Endocrinol (2016) 32:986–90. doi: 10.1080/09513590.2016.1197197

  • 12

    EgashiraNGotoYIbaHKawanakaRTakahashiRTaniguchiCet al. Kamishoyosan potentiates pentobarbital-induced sleep in socially isolated, ovariectomized mice. J Ethnopharmacol (2021) 281:114585. doi: 10.1016/j.jep.2021.114585

  • 13

    RuJLiPWangJZhouWLiBHuangCet al. TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J Cheminform (2014) 6:13. doi: 10.1186/1758-2946-6-13

  • 14

    SafranMDalahIAlexanderJRosenNIny SteinTShmoishMet al. GeneCards version 3: The human gene integrator. Database (Oxford) (2010) 2010):baq020. doi: 10.1093/database/baq020

  • 15

    ZhouYZhouBPacheLChangMKhodabakhshiAHTanaseichukOet al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun (2019) 10:1523. doi: 10.1038/s41467-019-09234-6

  • 16

    Ww PDBC. Protein data bank: the single global archive for 3D macromolecular structure data. Nucleic Acids Res (2019) 47:D520–D8. doi: 10.1093/nar/gky949

  • 17

    WeiFSongYGongAPanCZhuangYZhangXet al. Investigating the molecular mechanism of xijiao dihuang decoction for the treatment of SLE based on network pharmacology and molecular docking analysis. BioMed Res Int (2022) 2022:5882346. doi: 10.1155/2022/5882346

  • 18

    WangGZhouBWangZMengYLiuYYaoXet al. Pharmacological mechanisms underlying the anti-asthmatic effects of modified guomin decoction determined by network pharmacology and molecular docking. Front Mol Biosci (2021) 8:644561. doi: 10.3389/fmolb.2021.644561

  • 19

    RoccaWAMielkeMMGazzuola RoccaLStewartEA. Premature or early bilateral oophorectomy: A 2021 update. Climacteric (2021) 24:466–73. doi: 10.1080/13697137.2021.1893686

  • 20

    PerichTUssherJ. Stress predicts depression symptoms for women living with bipolar disorder during the menopause transition. Menopause (2021) 29:231–5. doi: 10.1097/GME.0000000000001894

  • 21

    StutePEversheimHOrtius-LechnerDMayMFeigC. Care reality of menopausal women in Germany: healthcare research using quantitative (SHI claims data) and qualitative (survey) data collection. Arch Gynecol Obstet (2022) 306:513–21. doi: 10.1007/s00404-022-06457-9

  • 22

    BotteriEXuZStoerNCMishraGD. Menopausal hormone therapy and melanoma risk in the Australian longitudinal study on women's health. Maturitas (2022) 160:13. doi: 10.1016/j.maturitas.2022.01.004

  • 23

    HonigbergMCTrinderMNatarajanP. Lipoprotein(a), menopausal hormone therapy, and risk of coronary heart disease in postmenopausal individuals. JAMA Cardiol (2022) 7:565–8. doi: 10.1001/jamacardio.2022.0716

  • 24

    HajializadehZKhaksariMNajafipourHSanjariMMahaniFDRaji-AmirhasaniA. Substitution of calorie restriction for protective effects of estrogen on cardiometabolic risk factors and oxidative stress in obese postmenopausal rat model. Life Sci (2022) 294:120367. doi: 10.1016/j.lfs.2022.120367

  • 25

    WeiMMahadyGBLiuDZhengZSLuY. Astragalin, a flavonoid from morus alba (Mulberry) increases endogenous estrogen and progesterone by inhibiting ovarian granulosa cell apoptosis in an aged rat model of menopause. Molecules (2016) 21:675. doi: 10.3390/molecules21050675

  • 26

    ShaoFLinWWangWWashingtonWCJr.ZhengL. The effect of emotional stress on the primary humoral immunity of rats. J Psychopharmacol (2003) 17:179–83. doi: 10.1177/0269881103017002005

  • 27

    YeXFuangfooTLiT. The effect of empty bottle stimulation on the bahaviours of rats and evaluation of open-filed test system. B Health Sci Technol (2015) 13:75–9.

  • 28

    YeYLiuCZWangRZhangYWZhangBCuiYet al. A new animal model for menopausal transition: combination of ovariectomy and empty bottle stimulation. Gynecol Endocrinol (2018) 34:840–4. doi: 10.1080/09513590.2018.1451835

  • 29

    KimJYChoSMYooYLeeTKimJK. Association between stroke and abdominal obesity in the middle-aged and elderly Korean population: KNHANES data from 2011-2019. Int J Environ Res Public Health (2022) 19:6140. doi: 10.3390/ijerph19106140

  • 30

    KimMHKimEJChoiYYHongJYangWM. Lycium chinense improves post-menopausal obesity via regulation of PPAR-gamma and estrogen receptor-alpha/beta expressions. Am J Chin Med (2017) 45:269–82. doi: 10.1142/S0192415X17500173

  • 31

    Ludgero-CorreiaAJr.AguilaMBMandarim-De-LacerdaCAFariaTS. Effects of high-fat diet on plasma lipids, adiposity, and inflammatory markers in ovariectomized C57BL/6 mice. Nutrition (2012) 28:316–23. doi: 10.1016/j.nut.2011.07.014

  • 32

    ZhouSZhaoLYiTWeiYZhaoX. Menopause-induced uterine epithelium atrophy results from arachidonic acid/prostaglandin E2 axis inhibition-mediated autophagic cell death. Sci Rep (2016) 6:31408. doi: 10.1038/srep31408

  • 33

    DarehzereshkiSDehghaniFEnjezabB. Mindfulness-based stress reduction group training improves of sleep quality in postmenopausal women. BMC Psychiatry (2022) 22:254. doi: 10.1186/s12888-022-03869-4

  • 34

    KravitzHMZhaoXBrombergerJTGoldEBHallMHMatthewsKAet al. Sleep disturbance during the menopausal transition in a multi-ethnic community sample of women. Sleep (2008) 31:979–90.

  • 35

    SiegelAMMathewsSB. Diagnosis and treatment of anxiety in the aging woman. Curr Psychiatry Rep (2015) 17:93. doi: 10.1007/s11920-015-0636-3

  • 36

    WallEHHewittSCCaseLKLinCYKorachKSTeuscherC. The role of genetics in estrogen responses: a critical piece of an intricate puzzle. FASEB J (2014) 28:5042–54. doi: 10.1096/fj.14-260307

  • 37

    ValeraMCGourdyPTremollieresFArnalJF. From the women's health initiative to the combination of estrogen and selective estrogen receptor modulators to avoid progestin addition. Maturitas (2015) 82:274–7. doi: 10.1016/j.maturitas.2015.07.012

  • 38

    WangJQianXGaoQLvCXuJJinHet al. Quercetin increases the antioxidant capacity of the ovary in menopausal rats and in ovarian granulosa cell culture in vitro. J Ovarian Res (2018) 11:51. doi: 10.1186/s13048-018-0421-0

  • 39

    BARRichardsonKAYangSPatelSFlawsJANowakRA. Effects of chronic dietary exposure to phytoestrogen genistein on uterine morphology in mice. J Agric Food Chem (2021) 69:1693–704. doi: 10.1021/acs.jafc.0c07456

  • 40

    WattelAKamelSMentaverriRLorgetFProuilletCPetitJPet al. Potent inhibitory effect of naturally occurring flavonoids quercetin and kaempferol on in vitro osteoclastic bone resorption. Biochem Pharmacol (2003) 65:3542. doi: 10.1016/s0006-2952(02)01445-4

  • 41

    WeiWCaoTPathakJLLiuXMaoTWatanabeNet al. Apigenin, a single active component of herbal extract, alleviates xerostomia via ERalpha-mediated upregulation of AQP5 activation. Front Pharmacol (2022) 13:818116. doi: 10.3389/fphar.2022.818116

  • 42

    YaoMLiuHLiBLiuYMuY. The relationship between earlier onset of natural menopause and elevated urinary albumin-creatinine ratio in postmenopausal Chinese women. Diabetes Metab Syndr Obes (2021) 14:847–56. doi: 10.2147/DMSO.S292041

  • 43

    SyedAARezaMISinghPHusainADadgeSGayenJR. Polyphenolic-rich cissus quadrangularis extract ameliorates insulin resistance by activating AdipoR1 in peri-/post-menopausal rats. Exp Gerontol (2022) 159:111681. doi: 10.1016/j.exger.2021.111681

  • 44

    SinatoraRVChagasEFBMatteraFOPMellemLJSantosAPereiraLPet al. Relationship of inflammatory markers and metabolic syndrome in postmenopausal women. Metabolites (2022) 12:73. doi: 10.3390/metabo12010073

  • 45

    HabibiPAhmadiaslNNourazarianAYousefiH. Swimming exercise improves SIRT1, NF-kappaB, and IL-1beta protein levels and pancreatic tissue injury in ovariectomized diabetic rats. Horm Mol Biol Clin Investig (2022) 43:345–52. doi: 10.1515/hmbci-2021-0069

  • 46

    XuHQiuQDuYLiuSHuangS. Exploration on effect of yinhuotang on depression-like behavior in mice with bilateral ovariectomy induced by chronic stress based on e 2 /ERβ signaling pathway. Chin J Exp Tradit Med Form (2021) 27:10–8. doi: 10.13422/j.cnki.syfjx.20211903

  • 47

    YeYLiQPanCSYanLSunKWangXYet al. QiShenYiQi inhibits tissue plasminogen activator-induced brain edema and hemorrhage after ischemic stroke in mice. Front Pharmacol (2021) 12:759027. doi: 10.3389/fphar.2021.759027

Summary

Keywords

menopause, Yin Huo Tang, ovariectomy, empty bottle stimulation, estradiol

Citation

Ye Y, Zhang B, Li Y, Xu H-D, Liu X-M, Huang S-M, Wang R and Li D (2022) Yin Huo Tang, a traditional Chinese herbal formula, relives ovariectomy and empty bottle stimulation-induced menopause-like symptoms in mice. Front. Endocrinol. 13:994642. doi: 10.3389/fendo.2022.994642

Received

15 July 2022

Accepted

06 October 2022

Published

19 October 2022

Volume

13 - 2022

Edited by

Wensheng Li, Sun Yat-sen University, China

Reviewed by

Farhath Sultana, Icahn School of Medicine at Mount Sinai, United States; María Belén Hapon, CONICET Instituto de Medicina y Biologia Experimental de Cuyo (IMBECU), Argentina

Updates

Copyright

*Correspondence: Rui Wang, ; Dong Li,

†These authors have contributed equally to this work and share first authorship

This article was submitted to Experimental Endocrinology, a section of the journal Frontiers in Endocrinology

Disclaimer

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

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics