Front. Pharmacol. Frontiers in Pharmacology Front. Pharmacol. 1663-9812 Frontiers Media S.A. 765638 10.3389/fphar.2021.765638 Pharmacology Original Research Integrated Network Pharmacology and GC-MS–Based Metabolomics to Investigate the Effect of Xiang-Su Volatile Oil Against Menopausal Depression Li et al. Xiang-Su Oil Against Menopausal Depression Li Yao 1 Yang Xinyi 1 Chen Shanshan 1 Wu Lei 2 Zhou Jinyong 3 Jia Keke 1 4 * Ju Wenzheng 1 * Department of Clinical Pharmacology, Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing, China Department of Pharmacy, Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing, China Central Laboratory, Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing, China School of Medicine and Holistic Integrative Medicine, Nanjing University of Chinese Medicine, Nanjing, China

Edited by: Yibin Feng, The University of Hong Kong, Hong Kong SAR, China

Reviewed by: Thangaraj Devadoss, KVSR Siddhartha College of Pharmaceutical Sciences, India

Saeid Abbasi Maleki, Kermanshah University of Medical Sciences, Iran

 *Correspondence: Wenzheng Ju, wzhju333@163.com; Keke Jia, jiakeke668@163.com

This article was submitted to Ethnopharmacology, a section of the journal Frontiers in Pharmacology

 02 12 2021 2021 12 765638 27 08 2021 08 11 2021 Copyright © 2021 Li, Yang, Chen, Wu, Zhou, Jia and Ju. 2021 Li, Yang, Chen, Wu, Zhou, Jia and Ju

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Menopausal depression perplexes a great number of women in later life. Xiangfu-Zisu (Xiang-Su), a traditional Chinese herbal pair composed of rhizomes of Cyperus rotundus L. (Xiangfu) and leaves of Perilla frutescens (L.) Britt. (Zisu), is frequently reported with antidepressant-like effects. The volatile oil from Xiangfu and Zisu has shown good antidepressant action, but its mechanism is still unclear. This study aimed to investigate the pharmacological mechanism of Xiang-Su (XS) volatile oil against menopausal depression through gas chromatography–mass spectrometry (GC-MS)-based network pharmacology and metabolomics. First, ADME screening was performed on actual detected components of XS volatile oil to obtain active constituents, and then duplicates of active constituent–related targets and menopausal depression–related targets were collected. These duplicates were considered as targets for XS volatile oil against menopausal depression, followed by GO and KEGG enrichment analyses. It showed that a total of 64 compounds were identified in XS volatile oil, and 38 active compounds were screened out. 42 overlapping genes between 144 compound-related genes and 780 menopausal depression–related genes were obtained. Results showed that targets of SLC6A4 and SLC6A3, regulation of serotonergic and dopaminergic synapses, were involved in the antidepressant mechanism of XS volatile oil. Next, antidepressant-like effect of XS volatile oil was validated in menopausal rats by ovariectomy (OVX) combined with chronic unpredictable mild stress (CUMS). Behavioral tests, biochemical analysis, and GC-MS–based non-targeted plasma metabolomics were employed to validate the antidepressant effect of XS volatile oil. Experimental evidence demonstrated that XS volatile oil reversed behavioral parameters in the sucrose preference test (SPT), open-field test (OFT), forced swim test (FST), and serum estradiol levels in OVX rats. Furthermore, results of metabolomics indicated that XS volatile oil mainly acts on regulating metabolic pathways of phenylalanine, tyrosine and tryptophan biosynthesis, tyrosine metabolism, and tryptophan metabolism, which were corresponding with the above-predicted results. These data suggest that network pharmacology combined with metabolomics provides deep insight into the antidepressant effect of XS volatile oil, which includes regulating key targets like SLC6A4 and SLC6A3, and pathways of serotonergic and dopaminergic synapses.

 Xiang-Su volatile oil menopause depression network pharmacology metabolomics BK20211394 81903823 2019M661903 2019K163 Natural Science Foundation of Jiangsu Province10.13039/501100004608 National Natural Science Foundation of China10.13039/501100001809 China Postdoctoral Science Foundation10.13039/501100002858 Jiangsu Postdoctoral Research Foundation10.13039/501100010011 Introduction

Menopause is defined as 12 months of amenorrhea following the final menstrual cycle with fluctuations of steroid hormone levels. During the menopause transition, women develop various symptoms such as sleep disturbances, hot flashes, or adverse mood, causing a high risk of depression (Tang et al., 2019, 2020; Zhou et al., 2021). However, current treatment strategies of menopausal depression including selective serotonin reuptake inhibitors (SSRIs) or hormone replacement often exhibit several adverse effects like withdrawal syndromes and sexual dysfunction in SSRIs (Khazaie et al., 2015; Moncrieff, 2019), or increased risk of cardiovascular events and breast cancer due to estrogen replacement (Anagnostis et al., 2019). Therefore, finding safe and effective drugs for the treatment of menopausal depression is of great urgency.

Traditional Chinese medicine (TCM) has a great potential to treat menopausal depression. For example, Jie-Yu Pill showed antidepressant-like effects in mice that experienced ovariectomy (OVX) with chronic unpredictable mild stress (CUMS) (Zhou et al., 2020), which was a common animal model to mimic clinical menopausal depression. Besides, essential oil from medicinal plants also has the potential to relieve depression and secondary depressive symptoms (de Sousa et al., 2017). In China, aromatic botanical drugs, rhizomes of Cyperus rotundus L. (Cyperaceae; Cyperi rhizoma, Xiangfu), leaves of Perilla frutescens (L.) Britt. (Lamiaceae; Perillae folium, Zisu), and formulas that contain these one or two botanical drugs are commonly used in treating depression or menopause-related syndrome. For example, a TCM formula Xiang-su-san (also named Koso-san in Japanese Kampo formula), which contains herbal pair Xiangfu-Zisu (Xiang-Su), had an antidepressant-like effect in mice (Ito et al., 2006). Cyperus rotundus L. is commonly used as a clinical herbal remedy in TCM prescription for treating depression (Zhao et al., 2015). It has also been frequently used to treat clinical gynecology disorders including premenstrual syndrome, primary dysmenorrheal, and polycystic ovary syndrome (Chen et al., 2014a, 2014b; Liao et al., 2018). Another aromatic botanical drug of Perilla frutescens (L.) Britt. showed antidepressant-like activities in CUMS-induced depressive mice (Yi et al., 2013; Ji et al., 2014a). Inspired by both the gynecological and neuroprotective effects of Cyperus rotundus L. and Perilla frutescens (L.) Britt., we speculate that volatile oil from the herbal pair Xiangfu-Zisu has the potential to treat menopausal depression and its following antidepressant-like effect was explored.

Network pharmacology is considered to be an appropriate approach for modern TCM pharmacological research (Zhang et al., 2019). The “compound-proteins/genes-disease” pathways are capable of describing complexities among biological systems, drugs, and diseases from a network perspective, sharing a similar holistic philosophy as TCM, and pointing to a new direction for the prediction of pharmacological mechanisms of TCM. It revealed that Cyperus rotundus L. showed an antidepressant effect by synergistically regulating multiple components, multiple targets, and multiple pathways through network pharmacology analysis (Jia et al., 2019). Another technology, metabolomics, makes a great contribution to understanding the basis of diseases and drug treatment. It shows advantages by integrating information from the final products of interactions among gene expression, protein function, and cellular environment (Rinschen et al., 2019). At present, metabolomics has become an important strategy for determining the antidepressant effect of TCM recipes. A classical TCM formula, Xiaoyaosan, was reported with therapeutic response in depressed patients by the metabolomics approach (Liu et al., 2015). Furthermore, the integrated strategy of metabolomics coupled with network pharmacology is an effective tool in illuminating the antidepressive action of TCM. For example, the antidepressant activity of Huang-Lian Jie-Du decoction was analyzed through network pharmacology combined with the metabolomics approach (Qu et al., 2021). However, one drawback of network pharmacology for TCM was that compounds were from some databases instead of actual identification.

Thus, in this study, we first detected the component composition of Xiang-Su (XS) volatile oil, based on which the absorption, distribution, metabolism, and excretion (ADME) screening was carried out, and network pharmacology was utilized to predict potential bioactive compounds and elucidate the molecular mechanisms for XS volatile oil against menopausal depression. Then the antidepressant-like effect of XS volatile oil was verified in rats induced by OVX with CUMS. Finally, a gas chromatography–mass spectrometry (GC-MS)-based plasma non-targeted metabolomics approach was applied to verify the underlying antidepressant effect of XS volatile oil. The technical strategy of this study is shown in Figure 1. This study will provide experimental evidence and strengthen our understanding of the antidepressant action of XS volatile oil in menopause.

Schematic diagram of the research.

 Materials and Methods Materials and Reagents

The rhizomes of Cyperus rotundus L. (batch No. 180801) and leaves of Perilla frutescens (L.) Britt. (batch No. 190301) collected from Anhui were purchased from Baicaotang of Nanjing University of Chinese Medicine (Nanjing, China), subsequently authenticated by Professor Yu Zhang in Nanjing University of Chinese Medicine (Nanjing, China). Anhydrous sodium sulfate was supplied by Nanjing Chemical Reagent Co., Ltd. (Nanjing, China). Estradiol valerate tablets were purchased from Bayer HealthCare Co., Ltd. (Lot No. 580B, Guangzhou, China). Penicillin was purchased from Harbin Pharmaceutical Group Holding Co., Ltd. (Heilongjiang, China). Tween 80 and isoflurane were purchased from Shanghai Yuanye Bio-Technology Co., Ltd. (Shanghai, China). N, O-bis(trimethylsilyl)trifluoroacetamide (BSTFA), 1,2-13C myristic acid, pyridine, and methoxyamine hydrochloride were purchased from Sigma-Aldrich, Inc. (St. Louis, MO, United States). Methanol (chromatography grade) and n-hexane (chromatography grade) were supplied by Thermo Fisher Scientific Inc. (Waltham, MA, United States). The distilled water was produced by a Milli-Q purification instrument (Milford, MA, United States).

 Preparation and Compound Identification of XS Volatile Oil Preparation of XS Volatile Oil

XS volatile oil consisted of two Chinese botanical drugs, namely, Cyperus rotundus L. and Perilla frutescens (L.) Britt. The volatile oil was extracted by steam distillation as follows. The rhizomes of Cyperus rotundus L. and leaves of Perilla frutescens (L.) Britt. were mixed at the ratio of 1:1 (w/w, with a total weight of 6 kg) and soaked in 8-fold of distilled water for 2 hours at room temperature, and then subjected to hydrodistillation for 3 hours to get the XS volatile oil according to the isolation procedure of volatile oil in Chinese Pharmacopoeia 2020 Edition. Later, the volatile oil was dried with anhydrous sodium sulfate and stored in brown glass at 4°C. The obtained volatile oil was weighted about 35 g. Finally, the yield of this volatile oil was 0.6% (w/w).

 Compound Composition of XS Volatile Oil

The compound composition of XS volatile oil was detected by GC-MS using an Agilent GC 7890B-7000C system (Agilent Technologies Company, United States) fitted with an HP-5 MS capillary column (30.0 m × 250 μm × 0.25 μm, Agilent 19091S-433). GC-MS detection conditions were as follows: carrier gas, 99.999% high-purity helium; flow rate, 1.0 ml/min; sample volume, 1.0 μl (3 mg/ml XS volatile oil of n-hexane solution); injection port and detector temperature, 220°C; split ratio, 20:1; electronic impact, 70 eV; ion source temperature, 230°C; quadrupole rod temperature, 150°C. The oven temperature program was initially set at 50°C for 3 min, ramped at 10°C/min to 140°C, then ramped at 3°C/min to 200°C, and finally ramped at 50°C/min to 230°C and held for 2 min. Compound identification was performed by comparing the spectra in the database of the National Institute of Standards and Technology (NIST).

 Network Pharmacology ADME Screening

Compounds identified from XS volatile oil were converted into the canonical simplified molecular-input line-entry system (SMILES) and screened with human gastrointestinal absorption (HIA), blood–brain barrier (BBB) permeation, and drug-likeness calculated by SwissADME (http://www.swissadme.ch/) (Daina et al., 2017). Parameters of HIA met “high,” BBB met “yes”, and two or more models among five drug-likeness models (Lipinski, Ghose, Veber, Egan, and Muegge) met “yes” were chosen as active compounds with good bioavailability.

 Targets Linked to Identified Compounds or Menopausal Depression

Targets of the identified compounds were predicted by SwissTargetPrediction (http://www.swisstargetprediction.ch/) (Daina et al., 2019). Meanwhile, targets of menopausal depression were obtained and screened with a score ≥5.0 by retrieving the keyword of “menopausal depression” from GeneCards (https://www.genecards.org/), a database integrating all annotated and predicted genes associated with human diseases (Stelzer et al., 2016).

 Network Construction, and GO and KEGG Enrichment Analyses

The overlapping targets between active compounds and menopausal depression were identified, and then the compound–target network was constructed and visualized by Cytoscape 3.6.0. Degree, betweenness centrality, and closeness centrality indicating the topological importance of nodes in the network were analyzed by Network Analyzer Tool in Cytoscape (Otasek et al., 2019). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed using the online functional annotation and enrichment tool DAVID (https://david.ncifcrf.gov/) (Huang et al., 2009), and plotted by the online platform for data analysis and visualization (http://www.bioinformatics.com.cn/). GO terms and KEGG pathways with a p-value < 0.05 were considered statistically significant.

 Animal Study Animals and Treatment

Female SD rats (180–200 g, supplied by Nantong University, Jiangsu, China) were housed and fed with unlimited access to food and water under a 12-hour light/dark cycle except for the following CUMS procedure. After adaptive feed, the rats were randomly divided into six groups as follows: sham group (negative control), OVX group (model control), OVX with estradiol valerate group (positive control, OVX-E2, 0.18 mg/kg), and OVX with different doses of XS volatile oil (OVX-XSL, 10.8 mg/kg; OVX-XSM, 32.4 mg/kg; OVX-XSH, 97.2 mg/kg) groups. The dose of XS volatile oil was converted by a drug-extract ratio of 0.6% (w/w) from the recommended dose of Cyperus rotundus L. and Perilla frutescens (L.) Britt for human in Chinese Pharmacopoeia 2020 Edition, and set in a safe dose range according to references (Jebasingh et al., 2012; Ji et al., 2014a; Kum et al., 2017; Kangwan et al., 2019). XS volatile oil was dissolved in distilled water with 0.5% Tween 80 to form an emulsion with a concentration of 2.0, 6.0, and 18.0 mg/ml, respectively, with gastric perfusion in a volume of 5.4 ml/kg once daily. Meanwhile, distilled water with 0.5% Tween 80 was used for the negative and model control groups. Estradiol valerate tablets used for the positive control group were dissolved in distilled water to a concentration of 0.036 mg/ml of estradiol valerate with gastric perfusion in a volume of 5.0 ml/kg. Bilateral OVX was performed in the rats in model control, positive control, and three XS volatile oil groups, respectively, while pseudo-operation was performed in the sham group. All surgical procedures were performed under isoflurane inhalation anesthesia. After surgery, penicillin was injected for three consecutive days, and all rats were kept in separate cages. Then, all OVX rats were experienced with the CUMS procedure, which was followed by Willner et al. with minor modifications (Willner et al., 1987). Random stressors included food or water deprivation (12 h), cage tilt (12 h), wet bedding (12 h), empty cage (12 h), tail pinch (2 min), physical restraint (2 h), light (24 h), stroboscopic light (12 h), and intermittent light (12 h). The animal study was reviewed and approved by the Ethical Committee of Jiangsu Province Hospital of Chinese Medicine and strictly followed the guidelines for the care and use of laboratory animals (National Research Council Committee for the Update of the Guide for the and Use of Laboratory, 2011). All the experimental procedures are shown in Figure 2.

Schematic representation of the experimental procedure.

 Behavioral Tests

Sucrose preference test (SPT): Before the test, rats were trained for 3 days to consume from two bottles of sucrose solution (1%, w/v) or one bottle of distilled water and one bottle of sucrose. The test was conducted after the rats were water-deprived for 18 hours. Two bottles of water and sucrose were given to rats for 1 hour, then following 1 hour with interchanged placement. Finally, the consumed weight was recorded, and the following formula was used to calculate the sucrose preference rate: sucrose preference rate (%) = sucrose consumption/(water consumption + sucrose consumption) × 100. The SPT was performed according to our previous study with minor modifications (Jing et al., 2019).

Forced swim test (FST): The day before the test, rats were forced to pre-swim for 15 min individually in plastic cylinders filled with water at a temperature of about 25°C to a depth of 30 cm. The cylinder was cleaned carefully after each test. After 24 hours, the immobile time of rats during 5 min of FST was recorded. Immobile time refers to the time during which rats stopped struggling, floated motionlessly, or only tried to keep their heads above water. The test session was recorded by a video camera and analyzed using the SuperFst system (XinRuan Inc., Ltd., Shanghai, China). The FST was performed according to the previous study with minor modifications (Pan et al., 2021).

Open-field test (OFT): The OFT was performed in an open-field apparatus (100 cm × 100 cm). All rats were placed individually in the center of the apparatus. The TopScanHR system (CleverSys Inc, Reston, VA, United States) was used to record the moving distance during 5 min in the apparatus. At each test interval, ethyl alcohol was used to clean the apparatus to avoid excretion interference. The OFT was performed according to the previous study with minor modifications (Pan et al., 2021).

 Biochemical Analysis

Before all rats were sacrificed, blood was first collected from the orbital sinus, and the serum was obtained by centrifugation (3,000 r/min, 15 min). Then the supernatant was separated and stored at −80°C before use. The serum levels of estradiol were measured using an enzyme-linked immunosorbent assay kit (Lot No. ZC-36464, ZCIBIO Technology Co. Ltd., Shanghai, China) according to the manufacturer’s instructions. Briefly, the serum samples were first diluted five times for measurement. Then 50 μl of serum samples as well as standards were added in a 96-well plate precoated with target antibody and followed by adding 100 μl horseradish peroxidase–conjugated antibody. After incubation for 60 min at 37°C, the liquid was discarded. Subsequently, each well was washed with 350 μl washing buffer for five times. 50 μl of substrate A and substrate B were then added into each well, which was next incubated at 37°C in dark for 15 min. Finally, 50 μl of termination solution was added to each well, and the absorbance was measured at 450 nm by a microplate reader. The serum estradiol levels were calculated according to the standard curve.

 Sample Preparation for GC-MS Metabolomics

Blood was also collected in heparin tubes from the abdominal aorta after isoflurane inhalation anesthesia and centrifuged to separate plasma for metabolomic analysis. The optimized method was based on the derivatization method by BSTFA. 50 μl of plasma was accurately taken and added with 200 μl methanol with 12.5 μg/ml 1,2-13C myristic acid as an internal standard. After 3-min vortex and centrifugation at 4°C, 18,000 r/min for 10 min, 150 μl supernatant was separated and concentrated in a low-temperature centrifugal concentrator (50°C, 2 hours) until the solvent was completely evaporated. The residue was dissolved in 45 μl of 10 mg/ml methoxyamine hydrochloride in a pyridine solution. After 5-min vortex and oscillation at 30°C for 1.5 hours (450 r/min), 45 μl of derivatization reagent BSTFA was precisely added. Then the sample was oscillated at 37°C for 0.5 hours (450 r/min) to full derivatization. Finally, the supernatant was collected after centrifugation at 4°C, 18,000 r/min for 10 min. The quality control (QC) sample was prepared by mixing the remaining supernatant after methanol precipitation of all samples, and the following operations were the same as described before. The QCs were injected at regular intervals (every 10 samples) to monitor the stability of the analytical process.

 GC-MS Analysis

The derivatized plasma samples were analyzed on the same equipment as compound identification of XS volatile oil. The detection conditions were as follows: carrier gas; 99.999% high-purity helium; flow rate, 1.0 ml/min; injection volume, 1.0 μl; injector temperature, 240°C; electronic impact, 70 eV; ion source temperature, 300°C; quadrupole rod temperature, 180°C. The column temperature was initially set at 60°C and held for 1 min, increased to 240°C at a rate of 20°C/min, then increased to 260°C at a rate of 5°C/min, and finally increased to 320°C at a rate of 30°C/min and held for 4 min. The solvent delay time was 3.95 min.

 Data Analysis

The GC-MS original data were processed through Agilent MassHunter Workstation software (Qualitative Analysis B.07.00). The database of NIST was used to characterize the metabolites. Then the data were normalized, and names, retention time, peak intensity, and mass of all samples were processed with the Microsoft Excel software (version 2007). Principal component analysis (PCA) and orthogonal partial least square discriminant analysis (OPLS-DA) were performed on the data by SIMCA-P 14.1 software. We also used the online tool Metaboanalyst 5.0 (https://www.metaboanalyst.ca/home.xhtml) to assist in finding metabolite differences and constructing metabolic pathways. All values were presented as the mean ± standard error of the mean (SEM). In the behavioral tests and biochemical analysis, the data were analyzed by one-way ANOVA followed by Fisher’s LSD test through GraphPad 6.0 software. p-value < 0.05 was considered statistically significant.

 Results Compound Identification of XS Volatile Oil and Active Component Screening

The actual chemical composition of XS volatile oil is shown in Table 1, and its chromatographic profile is shown in Supplementary Figure S1. The proportion of each component was obtained by peak area normalization. A total of 64 compounds were identified, accounting for 99.67% area of all the detected chromatographic peaks in XS volatile oil. Among them, perilla aldehyde (C11, 12.08%) from Perilla frutescens (L.) Britt., cyperenone (C54, 27.14%), cyperene (C19, 16.05%), α-cyperone (C60, 9.27%), and dehydrofukinone (C45, 4.95%) from Cyperus rotundus L. were the top 5 major components in this volatile oil, accounting for nearly 70% area of all the peaks. Then, we screened out 38 active compounds according to the parameters of HIA, BBB permeation, and drug-likeness in ADME screening (Table 2). Among these active compounds, perilla aldehyde (C11), dehydrofukinone (C45), cyperenone (C54), and α-cyperone (C60) were also included, which might account for the bioactivity of XS volatile oil.

Chemical composition of XS volatile oil.

No tR/min Compound Formula % Reverse match
C1 8.41 d-Limonene C10H16 1.15 892
C2 9.60 Linalool C10H18O 0.14 878
C3 10.30 L-Pinocarveol C10H16O 0.36 908
C4 11.00 Verbenone C10H14O 0.09 801
C5 11.09 α-Terpineol C10H18O 0.22 904
C6 11.19 (-)-Myrtenol C10H16O 0.76 867
C7 11.40 Berbenone C10H14O 0.17 907
C8 11.64 β-Cyclocitral C10H16O 0.07 846
C9 11.85 p-Cumic aldehyde C10H12O 1.01 908
C10 11.95 1-(Furan-2-yl)-4-methylpentan-1-one C10H14O2 0.47 913
C11 12.44 Perilla aldehyde C10H14O 12.08 920
C12 12.54 8-(1-Methylethylidene)bicycle[5.1.0]octane C11H18 0.24 866
C13 12.72 Perilla alcohol C10H16O 0.35 911
C14 13.15 Cyprotene C14H24 0.47 883
C15 13.62 Epoxycaryophyllene C15H24O 0.13 809
C16 13.86 2,4-Patchouladiene C15H22 1.01 931
C17 13.95 (+)-Cyclosativene C15H24 0.08 893
C18 14.08 Nootkatene C15H22 0.66 858
C19 14.60 Cyperene C15H24 16.05 951
C20 14.93 β-Caryophyllene C15H24 0.20 917
C21 14.97 9,10-dehydro-Isolongifolene C15H22 0.12 725
C22 15.23 Cypera-2,4(15)-diene C15H22 0.56 928
C23 15.50 α-Selinene C15H24 0.26 848
C24 15.75 Rotundene C15H24 2.30 910
C25 16.00 β-Vetispirene C15H22 0.12 850
C26 16.04 β-Gurjunene C15H24 0.27 890
C27 16.32 β-Selinene C15H24 2.23 935
C28 16.43 (+)-Valencene C15H24 0.29 943
C29 16.55 15-Hydroxy-α-muurolene C15H24O 0.12 805
C30 16.62 (+)-Isovalencenol C15H24O 0.12 857
C31 16.90 (-)-Nootkatene C15H22 0.28 887
C32 17.07 γ-Gurjunene C15H24 1.48 822
C33 17.30 Epoxycyperene C15H24O 1.29 888
C34 17.54 α-Calacorene C15H20 1.20 925
C35 17.80 10-Epi-Acora-3, 11-dien-15-al C15H22O 0.12 821
C36 17.99 Spathulenol C15H24O 0.41 780
C37 18.24 Aristol-1(10)-en-9-ol C15H24O 0.27 823
C38 18.52 1, 3-Di(propen-1-yl)adamantane C16H24 0.28 775
C39 18.56 Caryophyllene oxide C15H24O 0.10 902
C40 18.64 (-)-Spathulenol C15H24O 0.41 800
C41 18.81 Eudesma-4(15), 7-dien-1β -ol C15H24O 0.18 802
C42 18.96 4,6-diisopropylidene-8,8-dimethyl-Bicyclo[5.1.0]octan-2-one C16H24O 0.13 752
C43 19.18 cis-α-Copaene-8-ol C15H24O 1.65 837
C44 19.34 4,4,11,11-tetramethyl-7-Tetracyclo[6.2.1.0(3.8)0(3.9)]undecanol C15H24O 0.10 824
C45 19.53 Dehydrofukinone C15H22O 4.95 837
C46 19.65 γ-Gurjunenepoxide-(2) C15H24O 0.28 823
C47 19.89 Caryophylladienol II C15H24O 0.13 845
C48 20.10 13-nor-Eremophil-1(10)-en-11-one C14H22O 0.58 835
C49 20.16 3a,7,7-Trimethyltetrahydro-1H-cyclopropa[c]indene-2,3(1ah,3ah)-dione C13H18O2 0.73 838
C50 20.44 Isoaromadendrene epoxide C15H24O 0.51 835
C51 20.90 Calarene epoxide C15H24O 0.66 818
C52 21.02 Mustakone C15H22O 1.75 893
C53 21.26 Hexahydro-2,5,5-trimethyl-2H-,4a-ethanonaphthalen-8(5H)-one C15H24O 0.38 821
C54 21.58 Cyperenone C15H22O 27.14 923
C55 21.72 α-Costal C15H22O 0.51 875
C56 21.80 (-)-Rotundone C15H22O 0.21 906
C57 21.98 6-Isopropenyl-4,8a-dimethyl-3,5,6,7,8,8a-hexahydro-2(1H)-naphthalenone C15H22O 0.67 834
C58 22.33 Cyclocopacamphan-12-ol C15H24O 0.60 789
C59 22.82 Zizanal C15H22O 0.50 835
C60 23.04 α-Cyperone C15H22O 9.27 923
C61 23.22 Isovelleral C15H20O2 0.21 875
C62 23.40 Aristolone C15H22O 0.14 867
C63 24.50 Nootkatone C15H22O 0.27 913
C64 24.74 Furopelargone A C15H22O2 0.18 811

Compounds of XS volatile oil after ADME screening.

No Compound Formula SMILES
C2 Linalool C10H18O CC(=CCCC(C)(C=C)O)C
C3 L-Pinocarveol C10H16O CC1(C2CC1C(=C)C(C2)O)C
C4 Verbenone C10H14O CC1=CC(=O)C2CC1C2(C)C
C5 α-Terpineol C10H18O CC1=CCC(CC1)C(C)(C)O
C6 (-)-Myrtenol C10H16O CC1(C2CC=C(C1C2)CO)C
C7 Berbenone C10H14O CC1=CC(=O)[C@@H]2CC1C2(C)C
C8 β-Cyclocitral C10H16O CC1=C(C(CCC1)(C)C)C=O
C9 p-Cumic aldehyde C10H12O CC(C)C1=CC=C(C=C1)C=O
C10 1-(Furan-2-yl)-4-methylpentan-1-one C10H14O2 CC(C)CCC(=O)C1=CC=CO1
C11 Perilla aldehyde C10H14O CC(=C)C1CCC(=CC1)C=O
C13 Perilla alcohol C10H16O CC(=C)C1CCC(=CC1)CO
C15 Epoxycaryophyllene C15H24O CC1(CC2C1CCC3(C(O3)CCC2=C)C)C
C30 (+)-Isovalencenol C15H24O CC1CCC=C2C1(CC(=C(C)CO)CC2)C
C33 Epoxycyperene C15H24O CC1CCC2CC34C1(C2(C)C)CCC3(O4)C
C36 Spathulenol C15H24O CC1(C2C1C3C(CCC3(C)O)C(=C)CC2)C
C37 Aristol-1(10)-en-9-ol C15H24O CC1CC(CC2=C(CCC12)C)C=C(C)CO
C39 Caryophyllene oxide C15H24O CC1(CC2C1CCC3(C(O3)CCC2=C)C)C
C40 (-)-Spathulenol C15H24O CC1(C2C1C3C(CCC3(C)O)C(=C)CC2)C
C43 cis-α-Copaene-8-ol C15H24O CC1=CCC2C3C1C2(CC(C3C(C)C)O)C
C45 Dehydrofukinone C15H22O CC1CCCC2=CC(=O)C(=C(C)C)CC12C
C46 γ-Gurjunenepoxide-(2) C15H24O CC1CCC(C=C2C1CCC2C)C3(CO3)C
C47 Caryophylladienol II C15H24O CC1(CC2C1CCC(=C)C(CCC2=C)O)C
C48 13-nor-Eremophil-1(10)-en-11-one C14H22O CC1CCC=C2C1(CC(CC2)C(=O)C)C
C49 3a,7,7-Trimethyltetrahydro-1H-cyclopropa[c]indene-2,3(1ah,3ah)-dione C13H18O2 CC1(CCCC2(C13CC3C(=O)C2=O)C)C
C50 Isoaromadendrene epoxide C15H24O CC1CCC2C1C3C(C3(C)C)CC4C2(O4)C
C51 Calarene epoxide C15H24O CC1CCC2C3(C1(C4C(C4(C)C)CC3)C)O2
C52 Mustakone C15H22O CC1=CC(=O)C2C3C1C2(CCC3C(C)C)C
C54 Cyperenone C15H22O CC1CCC2CC3=C(C(=O)CC13C2(C)C)C
C55 α-Costal C15H22O CC1=CCCC2(C1CC(CC2)C(=C)C=O)C
C56 (-)-Rotundone C15H22O CC1CCC(CC2=C1C(=O)CC2C)C(=C)C
C57 6-Isopropenyl-4,8a-dimethyl-3,5,6,7,8,8a-hexahydro-2(1H)-naphthalenone C15H22O CC1=C2CC(CCC2(CC(=O)C1)C)C(=C)C
C58 Cyclocopacamphan-12-ol C15H24O CC(CO)C1CCC2(C3C1C4C2(C4C3)C)C
C59 Zizanal C15H22O CC1(C2CCC3(C2)C(CCC3C1=C)C=O)C
C60 α-Cyperone C15H22O CC1=C2CC(CCC2(CCC1=O)C)C(=C)C
C61 Isovelleral C15H20O2 CC1(CC2C=C(C3(CC3(C2C1)C)C=O)C=O)C
C62 Aristolone C15H22O CC1CCCC2=CC(=O)C3C(C12C)C3(C)C
C63 Nootkatone C15H22O CC1CC(=O)C=C2C1(CC(CC2)C(=C)C)C
C64 Furopelargone A C15H22O2 CC1CCC(C1C2=C(C=CO2)C(C)C)C(=O)C
Targets and Network Analysis of XS Volatile Oil Against Menopausal Depression

The SMILES of 38 abovementioned active compounds were plotted into SwissTargetPrediction to obtain the potential targets of each compound. Then 144 compound-related genes were collected after removing the duplicates (Supplementary Table S1). Besides, a total of 780 genes related to menopausal depression (score ≥ 5.0) were retrieved from the GeneCards database (Supplementary Table S2). Results showed that there were 42 overlapping genes by matching 144 compound-related genes with 780 disease-related genes (Figure 3A; Supplementary Table S3). A compound–target network of XS volatile oil against menopausal depression was constructed and visualized by importing 38 active compounds and 42 overlapping genes into Cytoscape 3.6.0. This network contained 80 nodes and 216 edges (Figure 3B). The contribution difference of these active compounds and genes to XS volatile oil against menopausal depression was shown in this network. According to the parameters of degree, betweenness centrality, and closeness centrality by topological analysis, cis-α-copaene-8-ol (C43), 1-(furan-2-yl)-4-methylpentan-1-one (C10), isovelleral (C61), calarene epoxide (C51), γ-gurjunenepoxide-(2) (C46), cyperenone (C54), dehydrofukinone (C45), (+)-isovalencenol (C30), α-cyperone (C60), and aristolone (C62) were the top 10 ingredient nodes linked to more targets. In addition, several targets such as cytochrome P450 family 19 subfamily A member 1 (CYP19A1), cytochrome P450 family 17 subfamily A member 1 (CYP17A1), solute carrier family six-member 4 (SLC6A4), and solute carrier family six-member 3 (SLC6A3) were considered as the crucial targets of XS volatile oil against menopausal depression ranked by degree, betweenness centrality, and closeness centrality.

Venn diagram of overlapping genes between compound and menopausal depression–related genes (A) and compound–target network of XS volatile oil against menopausal depression (B).

 GO and KEGG Enrichment Analyses of XS Volatile Oil Against Menopausal Depression

GO and KEGG enrichment analyses were performed on the abovementioned 42 targets for the treatment of XS volatile oil against menopausal depression. As shown in Figure 4A, the top 10 terms in the GO biological process (BP), cellular component (CC), and molecular function (MF) were ranked by p-value (Supplementary Table S4). It showed that these potential targets were mainly located in the neuronal cell body and pre-synapse, regulating molecular functions like steroid binding, enzyme binding, and monoamine transmembrane transporter activity, and participated in biological processes like steroid hormone-mediated signaling pathway and monoamine transport. Furthermore, the results of KEGG analysis indicated that 42 overlapping genes of XS volatile oil against menopausal depression were significantly enriched in 12 signaling pathways (p < 0.05), of which dopaminergic synapse, serotonergic synapse, and ovarian steroidogenesis were involved (Figure 4B; Supplementary Table S5). The crucial genes like CYP19A1 and CYP17A1 were involved in ovarian steroidogenesis signaling pathways, SLC6A4 was related to the serotonergic synapse, and SLC6A3 was related to the dopaminergic synapse. These data suggested that XS volatile oil may exert its antidepressant effect by regulating monoamine transport and pathways of dopaminergic and serotonergic synapses.

GO and KEGG enrichment analyses of 42 targets for XS volatile oil against menopausal depression. (A) Histogram of the top 10 terms in GO enrichment analysis. (B) Bubble chart of 12 signaling pathways in KEGG enrichment analysis. BP, biological process; CC, cellular component; MF, molecular function.

 XS Volatile Oil Prevented Depressive-like Behaviors in OVX Rats

As shown in Figure 5A, OVX resulted in a significant decrease in the ratio of the uterus to body weight (p < 0.001) compared to sham surgery. Behavioral tests of SPT, OFT, and FST were performed to determine the antidepressant effect of XS volatile oil. OVX rats showed depressive-like behaviors, including significantly decreased sucrose consumption in SPT (p < 0.01) (Figure 5B), reduced total traveled distance in OFT (p < 0.001) (Figure 5C), and increased immobile time in FST (p < 0.01) (Figure 5D), compared with the sham group. OVX rats with the administration of XS volatile oil (10.8, 32.4, 97.2 mg/kg) expressed significantly increased sucrose preference (p < 0.01) and traveled distance (10.8 and 32.4 mg/kg, p < 0.001; 97.2 mg/kg, p < 0.01). Additionally, immobile time in the high-dose group of XS volatile oil (97.2 mg/kg) was reduced compared with the OVX rats. These data showed the antidepressant effect of XS volatile oil in OVX rats. The positive drug estradiol (0.18 mg/kg) also presented antidepressant effect by improving sucrose preference in SPT (p < 0.01), increasing total traveled distance in OFT (p < 0.001), and decreasing immobile time in FST (p < 0.05).

Effects of XS volatile oil and estradiol treatment on behavior parameters in OVX rats. (A) Ratio of uterus to body weight. (B) Sucrose preference test (SPT). (C) Open-field test (OFT). (D) Forced swim test (FST). Data were represented as mean ± SEM; n = 6–10 per group. OVX-E2 (estradiol valerate, 0.18 mg/kg), OVX-XSL (10.8 mg/kg), OVX-XSM (32.4 mg/kg), OVX-XSH (97.2 mg/kg). ## p < 0.01, ### p < 0.001 compared with the sham group; * p < 0.05, ** p < 0.01, *** p < 0.001 compared with the OVX group.

 XS Volatile Oil Alleviated Serum Estradiol Deficiency in OVX Rats

The levels of serum estradiol decreased markedly in OVX rats compared with the sham rats (p < 0.001). Administration of XS volatile oil remarkably increased the serum estradiol levels in a dose-dependent manner (32.4 mg/kg, p < 0.01; 97.2 mg/kg, p < 0.001), while estradiol had the same regulatory effect (p < 0.001) in OVX rats (Figure 6).

Effects of XS volatile oil and estradiol treatment on serum estradiol levels in OVX rats. Data were represented as mean ± SEM; n = 8–10 per group. OVX-E2 (estradiol valerate, 0.18 mg/kg), OVX-XSL (10.8 mg/kg), OVX-XSM (32.4 mg/kg), OVX-XSH (97.2 mg/kg). ### p < 0.001 compared with the sham group; ** p < 0.01, *** p < 0.001 compared with the OVX group.

 XS Volatile Oil Improved Metabolic Profiles in OVX Rats

The unsupervised PCA was used to check the quality of data for the metabolomic analysis. As shown in Figure 7, samples from the sham, OVX, OVX-XSH, and QC groups were within the 95% Hotelling’s T-squared ellipse and were separated into clusters (R2X = 0.723, Q2 = 0.375). No outlier was found among these samples. The supervised OPLS-DA was performed to identify the metabolites responsible for the separation between OVX and the other two groups. All groups in the OPLS-DA models have met the 95% Hotelling’s T-squared ellipse and showed clear separation (Figure 8A). The 200 times permutation test was conducted to assess the predictive accuracy and statistical significance. Results from cross-validation suggested the model showed good predictability (Figure 8B). The metabolite peaks with VIP > 1.0 were selected.

Score plot of PCA. n = 8–10 per group, QC (quality control), and OVX-XSH (97.2 mg/kg).

Multivariate statistical analysis of metabolic characters of plasma samples acquired by GC-MS. (A) Score plot of OPLS-DA between different animal groups. (B) Results of cross-validation between different animal groups. n = 8–10 per group.

At the same time, metabolites were retained with p < 0.05 and fold change (FC) > 1.2 or <0.83 through the Metaboanalyst 5.0 database, and the filter rules were referred to other reports (Yuan et al., 2021). When the conditions of FC, p-value, and VIP-value were all satisfied, the metabolites were considered as differentially abundant. Metabolite changes induced by OVX and treatment of XS volatile oil are shown in Table 3 and Figure 9. We found the level of 21 metabolites significantly altered in the OVX group compared with the sham group, among which the levels of 2 metabolites were increased and those of 18 metabolites like l-tyrosine and l-tryptophan were reduced. XS volatile oil obviously recovered the levels of metabolites like l-tyrosine and l-tryptophan in OVX rats. Metabolic pathway analysis of these altered metabolites is shown in Figure 10. It showed that phenylalanine, tyrosine and tryptophan biosynthesis, tyrosine metabolism, and tryptophan metabolism were involved in the antidepressant effect of XS volatile oil in menopausal rats.

Significantly changed metabolites between different groups.

No Metabolites OVX/sham OVX-XSH/OVX
log2(FC) −log10(p) VIP Trend log2(FC) −log10(p) VIP Trend
1 l-Cystine −0.68 3.62 1.25 down 0.74 4.09 1.37 up
2 l-Lactic acid −1.26 5.91 1.34 down 1.29 4.70 1.36 up
3 l-Alanine −2.04 5.72 1.39 down 2.04 4.44 1.38 up
4 Epinephrine 0.46 3.50 1.18 up −0.67 4.32 1.33 down
5 4-Hydroxybutyric acid — — — — 0.83 2.83 1.01 up
6 l-Valine −1.46 4.60 1.37 down 1.49 3.94 1.37 up
7 Urea −0.86 3.94 1.25 down 0.44 1.55 1.08 up
8 l-Proline −1.61 4.39 1.36 down 1.63 3.75 1.34 up
9 l-Serine −0.95 3.80 1.30 down 0.91 2.90 1.28 up
10 l-Threonine −1.16 4.38 1.26 down 0.49 1.55 1.06 up
11 l-Asparagine −0.89 3.25 1.22 down — — — —
12 d-Arabinose 0.86 2.47 1.02 up −1.00 2.99 1.17 down
13 N-Acetylneuraminic acid −2.32 5.83 1.38 down 2.47 4.70 1.37 up
14 l-Glutamine −1.05 3.54 1.27 down 1.04 3.05 1.33 up
15 d-Galactose −0.59 4.38 1.26 down 0.57 2.58 1.14 up
16 l-Lysine −1.58 4.27 1.33 down 1.19 2.52 1.21 up
17 d-Glucose −0.85 2.80 1.19 down 0.75 2.27 1.18 up
18 l-Tyrosine −0.41 1.94 1.02 down 0.53 1.92 1.12 up
19 d-Xylulose −1.24 3.30 1.24 down 1.24 2.99 1.26 up
20 Stearic acid −0.29 3.27 1.14 down — — — —
21 L-Tryptophan −1.14 5.21 1.32 down 0.78 2.25 1.08 up

Discriminating metabolites differentiated in sham, OVX, and OVX-XSH groups. FC > 1.2 or <0.83, p-value < 0.05 and VIP-value > 1.0, n = 8–10 per group, OVX-XSH (97.2 mg/kg).

Metabolic pathways involved in the treatment of XS volatile oil on OVX rats. (1) Phenylalanine, tyrosine, and tryptophan biosynthesis; (2) glycine, serine, and threonine metabolism; (3) aminoacyl-tRNA biosynthesis; (4) tyrosine metabolism; (5) tryptophan metabolism; (6) Pentose and glucuronate interconversions; (7) alanine, aspartate, and glutamate metabolism.

 Discussion

Pathological studies have shown that menopausal syndrome mainly includes an imbalance of neurotransmitters, hormones, cytokines, and immune system (Stefanska et al., 2015), which increased the risk of depression. During menopause, estrogen levels were fluctuated, inducing depression and depressive-like behavior through interactions with neurotrophic factors and the serotonergic system. However, there is a lack of safe and effective drugs to treat menopausal depression.

TCM has been a great treasure in medical practice from Chinese history. The aromatic herbal pair of XS showed potential to relieve menopausal depression in previous studies (Ito et al., 2006; Ji et al., 2014a; Zhao et al., 2015), but its mechanism was poorly understood. With the development of modern science and technology, network pharmacology combined with metabolomics has been a promising strategy to illustrate the pharmacological mechanism of TCM, using which antidepressant actions of Xiaoyaosan and Huang-Lian Jie-Du decoction were illustrated (Liu et al., 2021; Qu et al., 2021). Thus, the present study aimed to study the antidepressant mechanism of XS volatile oil in menopausal rats by integrating network pharmacology with metabolomics.

In this study, a total of 64 constituents were identified from XS volatile oil. Perilla aldehyde (C11) from Perilla frutescens (L.) Britt., and cyperenone (C54), cyperene (C19), α-cyperone (C60), and dehydrofukinone (C45) from Cyperus rotundus L. were accounted for nearly 70% of all the detected chromatographic peaks. 38 active compounds like perilla aldehyde (C11), dehydrofukinone (C45), cyperenone (C54), and α-cyperone (C60) were further obtained by ADME screening, and they were predicted to be acting on 144 targets through the SwissTargetPrediction database. Among these targets, 42 targets were shared between compound-related and menopausal depression–related targets, indicating the possible antidepressant targets of XS volatile oil. Then the compound–target network of XS volatile oil against menopausal depression was constructed. According to the degree, betweenness centrality, and closeness centrality by topological analysis in this network, the top 10 ingredient nodes included dehydrofukinone (C45), cyperenone (C54), and α-cyperone (C60). Moreover, according to the previous studies, the sesquiterpenoid dehydrofukinone had sedative, anesthetic, and anticonvulsant effects through GABAergic or cortisol mechanisms (Garlet et al., 2016; Garlet et al., 2017). It was reported with anxiolytic-like effect in mice due to positive modulation of GABAA receptors and/or inhibition of neuronal calcium influx (Garlet et al., 2019). α-Cyperone was reported to exert antidepressant-like actions in a mice depression model which may be attributed to suppressing NLR family pyrin domain containing 3 (NLRP3) inflammasome (Xia et al., 2020).

Results of network pharmacology also suggested that targets such as CYP19A1, CYP17A1, SLC6A4, and SLC6A3 were the crucial targets of XS volatile oil against menopausal depression ranked by degree, betweenness centrality, and closeness centrality. Go and KEGG enrichment analyses showed that XS volatile oil may exert its antidepressant effect by regulating monoamine transport and pathways of dopaminergic and serotonergic synapses. Moreover, SLC6A4 was involved in the pathway of the serotonergic synapse, and SLC6A3 was involved in the pathway of the dopaminergic synapse. Metabolomic analysis of the hippocampus in a rat model of CUMS-induced depression demonstrated that some altered metabolites were related to the serotonergic synapse, dopaminergic synapse, and glutamatergic synapse, which were involved in the pathology of depression (Gao et al., 2021). SLC6A4 (a serotonin transporter gene) is one of the major determinants of serotonergic neurotransmission. It terminates neurotransmission by transporting serotonin from the synapse into the pre-synaptic nerve terminal. A comprehensive meta-analysis showed that serotonin transporter availability in depressed patients was reduced in key regions of the limbic system compared with healthy controls (Kambeitz and Howes, 2015). SLC6A4 methylation was found to be positively correlated with stress and depression, and possessed the potential for the diagnosis and treatment of major depression (Okada et al., 2014; Park et al., 2019). SLC6A3 (a dopamine transporter gene) affects the function of the dopamine nerve system. In geriatric patients with severe major depressive disorder, a low level of dopamine transporter binding was found in the region of the nucleus accumbens and putamen (Moriya et al., 2020). Alteration of dopamine transporter density in depressed patients with anhedonia showed reduced dopamine concentration in the synaptic cleft (Sarchiapone et al., 2006). CYP19A1-encoding aromatase is responsible for the key step in the synthesis of estradiol, which is widely expressed in the brain (Blakemore and Naftolin, 2016). Depletion of brain estrogen in middle-aged aromatase gene knockout (Ar−/-) mice increased the depressive-like behavior (Ma et al., 2020). All these results suggested that monoamine transport, and serotonergic and dopaminergic synapse pathways may be involved in the antidepressant action of XS volatile oil through network pharmacology analysis.

To verify the antidepressant effect of XS volatile oil, OVX combined with CUMS was established in rats to induce menopausal depression. As previously reported by others, OVX with CUMS resulted in several depressive-like behaviors in animals (Zhang et al., 2020; Zhou et al., 2020). In the present study, OVX induced a reduction of sucrose preference rate in SPT, a decrease of total traveled distance in OFT, and an increase of immobile time in FST, which could be reversed by XS volatile oil. Besides, XS volatile oil also upregulated the OVX-induced reduction of serum estradiol levels in rats. These results demonstrated the antidepressant-like effect of XS volatile oil.

In the following GC-MS–based non-targeted plasma metabolomic analysis, we found the level of 21 metabolites significantly altered in the OVX group compared with the sham group, among which the levels of two metabolites were increased and those of 18 metabolites like l-tyrosine and l-tryptophan were reduced. The results of reduced levels of l-tyrosine and l-tryptophan in this depressive animal model were in accordance with the results of CUMS depressive rats (Han et al., 2019), or the results of OVX-induced menopause (Lee et al., 2016). After XS volatile oil administration, l-tyrosine and l-tryptophan levels were upregulated in OVX rats. Metabolic pathway analysis was performed on these altered metabolites. It revealed that XS volatile oil may exert its antidepressant effect mainly by regulating phenylalanine, tyrosine and tryptophan biosynthesis, tyrosine metabolism, and tryptophan metabolism.

As we all know, amino acids have proven to be associated with the pathophysiological mechanism of depression (Chen et al., 2021). Tyrosine is one of the amino acids involved in dopamine biosynthesis. The dopamine system has been implicated in many different aspects of brain function, including locomotion, affect, and cognition. Many of the depressive symptoms, such as anhedonia and amotivation, have been more consistently associated with dysfunctions in the dopamine system (Grace, 2016). The decrease of dopamine release into the synapse caused altered expression of dopamine receptors within limbic structures, and that was observed in different depression models (Belujon and Grace, 2017). Similarly, tryptophan is a precursor highly essential for the brain regional synthesis of 5-hydroxytryptamine (5-HT, also known as serotonin). Patients with major depression had been confirmed with a deficit in circulation tryptophan, tryptophan uptake, serotonin synthesis and serotonergic neurotransmission (Ruddick et al., 2006; Gibson, 2018), and with decreased levels of brain serotonin and alterations in 5-HT receptors (Maes et al., 2009).

As reported by others, Perilla frutescens (L.) Britt. and Cyperus rotundus L. had the potential to protect against depression. For example, the essential oil of Perilla frutescens (L.) Britt effectively reversed the CUMS-induced reduction of brain-derived neurotrophic factors and improved the alterations of 5-HT concentrations in the hippocampus (Yi et al., 2013; Ji et al., 2014a). Perilla aldehyde, one of the active components in Perilla frutescens (L.) Britt, also exhibited an antidepressant-like effect in the CUMS-induced rat model of depression (Song et al., 2018). In lipopolysaccharide-induced depressive mice, pretreatment with perilla aldehyde reversed the decreased 5-HT levels in the prefrontal cortex, suggesting the antidepressant activity of perilla aldehyde might be related to the alteration of monoaminergic responses (Ji et al., 2014b). Cyperus rotundus L. was also frequently used in TCM prescription for treating depression in the clinic (Zhao et al., 2015). α-Cyperone, one of the active components in Cyperus rotundus L., also showed antidepressant action in a mice depression model (Xia et al., 2020).

The combination of network pharmacology and metabolomics has been an effective strategy to illustrate the antidepressive mechanism of TCM. For example, different efficacy groups of Xiaoyaosan showed synergistic antidepressant effects and contributed to the whole prescription against depression in CUMS rats (Liu et al., 2021). Huang-Lian Jie-Du decoction exhibited antidepressant effects by regulating SLC6A4 and monoamine oxidase A (MAOA) in the tryptophan metabolism of CUMS depressive mice (Qu et al., 2021). In this study, integrated results of metabolomics and network pharmacology analysis indicated that XS volatile oil prevented depression through pathways including phenylalanine, tyrosine and tryptophan biosynthesis, tyrosine metabolism, and tryptophan metabolism, thus affecting serotonergic and dopaminergic synapses.

 Conclusion

In this study, through integrated methods of metabolomics and network pharmacology analysis, we found that XS volatile oil prevented the depressive-like behavior in OVX rats through regulating pathways including phenylalanine, tyrosine and tryptophan biosynthesis, tyrosine metabolism, and tryptophan metabolism to restore serotonergic and dopaminergic synapse. It preliminarily revealed the multi-compounds, multi-targets, and multi-mechanisms of XS volatile oil acting on menopausal depression. However, several limitations must be noted, such as the current results only exhibited metabolite changes in plasma, while changes in brain regions need investigation to further verify the antidepressant effect of XS volatile oil. Moreover, the verification of active components, targets, and pathways for XS volatile oil against menopausal depression could be complemented in the future study, such as detection of components in plasma, mRNA, and protein levels of SLC6A4 and SLC6A3.

 Data Availability Statement

The original data 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 the Ethical Committee of Jiangsu Province Hospital of Chinese Medicine.

 Author Contributions

WJ, KJ, and YL designed the experiments; YL, XY, SC, and KJ performed the experiment; LW and JZ helped in the animal study; YL collected and analyzed the data; YL and KJ wrote and revised the manuscript. All the authors approved the final manuscript.

 Funding

This work was supported by contract/grant sponsor: Leading Talents of scientific research in TCM of Jiangsu Province (No. SLJ0208), the Natural Science Foundation of Jiangsu Province (No. BK20211394), the National Natural Science Foundation of China (No. 81903823), and the China Postdoctoral Science Foundation (No. 2019M661903), Jiangsu Planned Projects for Postdoctoral Research Funds, China (No. 2019K163).

 Conflict of Interest

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

 Publisher’s Note

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

Thanks for the help of Analysis and Measurement Center (School of Pharmacy, Nanjing University of Chinese Medicine, China) and Pharmacological Laboratory (Jiangsu Province Hospital of Chinese Medicine, China) for providing the equipment and technical support.

 Supplementary Material

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

References Anagnostis P. Paschou S. A. Katsiki N. Krikidis D. Lambrinoudaki I. Goulis D. G. (2019). Menopausal Hormone Therapy and Cardiovascular Risk: Where Are We Now? Curr. Vasc. Pharmacol. 17, 564–572. 10.2174/1570161116666180709095348 Belujon P. Grace A. A. (2017). Dopamine System Dysregulation in Major Depressive Disorders. Int. J. Neuropsychopharmacol. 20, 1036–1046. 10.1093/ijnp/pyx056 Blakemore J. Naftolin F. (2016). Aromatase: Contributions to Physiology and Disease in Women and Men. Physiology (Bethesda) 31, 258–269. 10.1152/physiol.00054.2015 Chen H. Xie H. Huang S. Xiao T. Wang Z. Ni X. (2021). Development of Mass Spectrometry-Based Relatively Quantitative Targeted Method for Amino Acids and Neurotransmitters: Applications in the Diagnosis of Major Depression. J. Pharm. Biomed. Anal. 194, 113773. 10.1016/j.jpba.2020.113773 Chen H. Y. Huang B. S. Lin Y. H. Su I. H. Yang S. H. Chen J. L. (2014a). Identifying Chinese Herbal Medicine for Premenstrual Syndrome: Implications from a Nationwide Database. BMC Complement. Altern. Med. 14, 206. 10.1186/1472-6882-14-206 Chen H. Y. Lin Y. H. Su I. H. Chen Y. C. Yang S. H. Chen J. L. (2014b). Investigation on Chinese Herbal Medicine for Primary Dysmenorrhea: Implication from a Nationwide Prescription Database in Taiwan. Complement. Ther. Med. 22, 116–125. 10.1016/j.ctim.2013.11.012 Daina A. Michielin O. Zoete V. (2017). SwissADME: A Free Web Tool to Evaluate Pharmacokinetics, Drug-Likeness and Medicinal Chemistry Friendliness of Small Molecules. Sci. Rep. 7, 42717. 10.1038/srep42717 Daina A. Michielin O. Zoete V. (2019). SwissTargetPrediction: Updated Data and New Features for Efficient Prediction of Protein Targets of Small Molecules. Nucleic Acids Res. 47, W357–W364. 10.1093/nar/gkz382 de Sousa D. P. Silva R. H. N. Silva E. F. D. Gavioli E. C. (2017). Essential Oils and Their Constituents: An Alternative Source for Novel Antidepressants. Molecules 22, 1290. 10.3390/molecules22081290 Gao Y. Mu J. Xu T. Linghu T. Zhao H. Tian J. (2021). Metabolomic Analysis of the Hippocampus in a Rat Model of Chronic Mild Unpredictable Stress-Induced Depression Based on a Pathway Crosstalk and Network Module Approach. J. Pharm. Biomed. Anal. 193, 113755. 10.1016/j.jpba.2020.113755 Garlet Q. I. Pires L. C. Silva D. T. Spall S. Gressler L. T. Bürger M. E. (2016). Effect of (+)-Dehydrofukinone on GABAA Receptors and Stress Response in Fish Model. Braz. J. Med. Biol. Res. 49, e4872. 10.1590/1414-431X20154872 Garlet Q. I. Pires L. D. C. Milanesi L. H. Marafiga J. R. Baldisserotto B. Mello C. F. (2017). (+)-Dehydrofukinone Modulates Membrane Potential and Delays Seizure Onset by GABAa Receptor-Mediated Mechanism in Mice. Toxicol. Appl. Pharmacol. 332, 52–63. 10.1016/j.taap.2017.07.010 Garlet Q. I. Rodrigues P. Barbosa L. B. Londero A. L. Mello C. F. Heinzmann B. M. (2019). Nectandra Grandiflora Essential Oil and its Isolated Sesquiterpenoids Minimize Anxiety-Related Behaviors in Mice through GABAergic Mechanisms. Toxicol. Appl. Pharmacol. 375, 64–80. 10.1016/j.taap.2019.05.003 Gibson E. L. (2018). Tryptophan Supplementation and Serotonin Function: Genetic Variations in Behavioural Effects. Proc. Nutr. Soc. 77, 174–188. 10.1017/S0029665117004451 Grace A. A. (2016). Dysregulation of the Dopamine System in the Pathophysiology of Schizophrenia and Depression. Nat. Rev. Neurosci. 17, 524–532. 10.1038/nrn.2016.57 Han X. M. Qin Y. J. Zhu Y. Zhang X. L. Wang N. X. Rang Y. (2019). Development of an Underivatized LC-MS/MS Method for Quantitation of 14 Neurotransmitters in Rat Hippocampus, Plasma and Urine: Application to CUMS Induced Depression Rats. J. Pharm. Biomed. Anal. 174, 683–695. 10.1016/j.jpba.2019.06.043 Huang D. W. Sherman B. T. Lempicki R. A. (2009). Systematic and Integrative Analysis of Large Gene Lists Using DAVID Bioinformatics Resources. Nat. Protoc. 4, 44–57. 10.1038/nprot.2008.211 Ito N. Nagai T. Yabe T. Nunome S. Hanawa T. Yamada H. (2006). Antidepressant-Like Activity of a Kampo (Japanese Herbal) Medicine, Koso-San (Xiang-Su-San), and its Mode of Action via the Hypothalamic-Pituitary-Adrenal Axis. Phytomedicine 13, 658–667. 10.1016/j.phymed.2006.01.002 Jebasingh D. Jackson D. D. Venkataraman S. Emerald B. S. (2012). Physiochemical and Toxicological Studies of the Medicinal Plant Cyperus Rotundus L (Cyperaceae). Int. J. Appl. Res. Nat. Prod. 5, 1–8. Available at: https://www.researchgate.net/publication/287430015_Physiochemical_and_toxicological_studies_of_the_medicinal_Plant_Cyperus_rotundus_L_Cyperaceae . Ji W. W. Li R. P. Li M. Wang S. Y. Zhang X. Niu X. X. (2014a). Antidepressant-Like Effect of Essential Oil of Perilla Frutescens in a Chronic, Unpredictable, Mild Stress-Induced Depression Model Mice. Chin. J. Nat. Med. 12, 753–759. 10.1016/s1875-5364(14)60115-1 Ji W. W. Wang S. Y. Ma Z. Q. Li R. P. Li S. S. Xue J. S. (2014b). Effects of Perillaldehyde on Alternations in Serum Cytokines and Depressive-like Behavior in Mice after Lipopolysaccharide Administration. Pharmacol. Biochem. Behav. 116, 1–8. 10.1016/j.pbb.2013.10.026 Jia H. Tang C. Liu H. Meng Y. U. Zhang H. Zou Z. (2019). Exploration of Mechanism of Antidepressant of Cyperus Rotundus Based on Network Pharmacology. Drug Eval. Res. 42, 49–55. 10.7501/j.issn.1674-6376.2019.01.007 Jing W. Song S. Sun H. Chen Y. Zhao Q. Zhang Y. (2019). Mahuang-Fuzi-Xixin Decoction Reverses Depression-like Behavior in LPS-Induced Mice by Regulating NLRP3 Inflammasome and Neurogenesis. Neural Plast. 2019, 1571392. 10.1155/2019/1571392 Kambeitz J. P. Howes O. D. (2015). The Serotonin Transporter in Depression: Meta-Analysis of In Vivo and Post Mortem Findings and Implications for Understanding and Treating Depression. J. Affect. Disord. 186, 358–366. 10.1016/j.jad.2015.07.034 Kangwan N. Pintha K. Lekawanvijit S. Suttajit M. (2019). Rosmarinic Acid Enriched Fraction from Perilla Frutescens Leaves Strongly Protects Indomethacin-Induced Gastric Ulcer in Rats. Biomed. Res. Int. 2019, 9514703. 10.1155/2019/9514703 Khazaie H. Rezaie L. Rezaei Payam N. Najafi F. (2015). Antidepressant-induced Sexual Dysfunction during Treatment with Fluoxetine, Sertraline and Trazodone; a Randomized Controlled Trial. Gen. Hosp. Psychiatry 37, 40–45. 10.1016/j.genhosppsych.2014.10.010 Kum C. J. Kim E. Y. Kim J. H. Lee B. Min J. H. Heo J. (2017). Cyperus Rotundus L. Extract Suppresses RANKL-Induced Osteoclastogenesis through NFATc1/c-Fos Downregulation and Prevent Bone Loss in OVX-Induced Osteoporosis Rat. J. Ethnopharmacol. 205, 186–194. 10.1016/j.jep.2017.03.017 Lee M. Y. Kim H. Y. Singh D. Yeo S. H. Baek S. Y. Park Y. K. (2016). Metabolite Profiling Reveals the Effect of Dietary Rubus Coreanus Vinegar on Ovariectomy-Induced Osteoporosis in a Rat Model. Molecules 21, 149. 10.3390/molecules21020149 Liao W. T. Chiang J. H. Li C. J. Lee M. T. Su C. C. Yen H. R. (2018). Investigation on the Use of Traditional Chinese Medicine for Polycystic Ovary Syndrome in a Nationwide Prescription Database in Taiwan. J. Clin. Med. 7, 179. 10.3390/jcm7070179 Liu C. C. Wu Y. F. Feng G. M. Gao X. X. Zhou Y. Z. Hou W. J. (2015). Plasma-Metabolite-Biomarkers for the Therapeutic Response in Depressed Patients by the Traditional Chinese Medicine Formula Xiaoyaosan: A (1)H NMR-Based Metabolomics Approach. J. Affect. Disord. 185, 156–163. 10.1016/j.jad.2015.05.005 Liu X. J. Wang Y. Z. Wei F. X. Lv M. Qu P. Chen S. J. (2021). The Synergistic Anti-depression Effects of Different Efficacy Groups of Xiaoyaosan as Demonstrated by the Integration of Network Pharmacology and Serum Metabolomics. J. Pharm. Biomed. Anal. 197, 113949. 10.1016/j.jpba.2021.113949 Ma L. Xu Y. Zhou J. Li Y. Zhang X. Jiang W. (2020). Brain Estrogen Alters the Effects of the Antidepressant Sertraline in Middle-Aged Female and Male Mice. Mol. Cel. Endocrinol. 516, 110947. 10.1016/j.mce.2020.110947 Maes M. Yirmyia R. Noraberg J. Brene S. Hibbeln J. Perini G. (2009). The Inflammatory & Neurodegenerative (I&ND) Hypothesis of Depression: Leads for Future Research and New Drug Developments in Depression. Metab. Brain Dis. 24, 27–53. 10.1007/s11011-008-9118-1 Moncrieff J. (2019). Persistent Adverse Effects of Antidepressants. Epidemiol. Psychiatr. Sci. 29, e56. 10.1017/S2045796019000520 Moriya H. Tiger M. Tateno A. Sakayori T. Masuoka T. Kim W. (2020). Low Dopamine Transporter Binding in the Nucleus Accumbens in Geriatric Patients with Severe Depression. Psychiatry Clin. Neurosci. 74, 424–430. 10.1111/pcn.13020 National Research Council (US) Committee for the Update of the Guide for the Care and Use of Laboratory Animals (2011). Guide for the Care and Use of Laboratory Animals. Washington, DC: National Academies Press (US). Okada S. Morinobu S. Fuchikami M. Segawa M. Yokomaku K. Kataoka T. (2014). The Potential of SLC6A4 Gene Methylation Analysis for the Diagnosis and Treatment of Major Depression. J. Psychiatr. Res. 53, 47–53. 10.1016/j.jpsychires.2014.02.002 Otasek D. Morris J. H. Bouças J. Pico A. R. Demchak B. (2019). Cytoscape Automation: Empowering Workflow-Based Network Analysis. Genome Biol. 20, 185. 10.1186/s13059-019-1758-4 Pan S. M. Pan Y. Tang Y. L. Zuo N. Zhang Y. X. Jia K. K. (2021). Thioredoxin Interacting Protein Drives Astrocytic Glucose Hypometabolism in Corticosterone‐induced Depressive State. J. Neurochem. 10.1111/jnc.15489 Park C. Rosenblat J. D. Brietzke E. Pan Z. Lee Y. Cao B. (2019). Stress, Epigenetics and Depression: A Systematic Review. Neurosci. Biobehav. Rev. 102, 139–152. 10.1016/j.neubiorev.2019.04.010 Qu S. Y. Li X. Y. Heng X. Qi Y. Y. Ge P. Y. Ni S. J. (2021). Analysis of Antidepressant Activity of Huang-Lian Jie-Du Decoction through Network Pharmacology and Metabolomics. Front. Pharmacol. 12, 619288. 10.3389/fphar.2021.619288 Rinschen M. M. Ivanisevic J. Giera M. Siuzdak G. (2019). Identification of Bioactive Metabolites Using Activity Metabolomics. Nat. Rev. Mol. Cel Biol. 20, 353–367. 10.1038/s41580-019-0108-4 Ruddick J. P. Evans A. K. Nutt D. J. Lightman S. L. Rook G. A. Lowry C. A. (2006). Tryptophan Metabolism in the Central Nervous System: Medical Implications. Expert Rev. Mol. Med. 8, 1–27. 10.1017/S1462399406000068 Sarchiapone M. Carli V. Camardese G. Cuomo C. Di Giuda D. Calcagni M. L. (2006). Dopamine Transporter Binding in Depressed Patients with Anhedonia. Psychiatry Res. 147, 243–248. 10.1016/j.pscychresns.2006.03.001 Song Y. Sun R. Ji Z. Li X. Fu Q. Ma S. (2018). Perilla Aldehyde Attenuates CUMS-Induced Depressive-like Behaviors via Regulating TXNIP/TRX/NLRP3 Pathway in Rats. Life Sci. 206, 117–124. 10.1016/j.lfs.2018.05.038 Stefanska A. Bergmann K. Sypniewska G. (2015). Metabolic Syndrome and Menopause: Pathophysiology, Clinical and Diagnostic Significance. Adv. Clin. Chem. 72, 1–75. 10.1016/bs.acc.2015.07.001 Stelzer G. Rosen N. Plaschkes I. Zimmerman S. Twik M. Fishilevich S. (2016). The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr. Protoc. Bioinformatics 54, 1.30.1–33.30.33. 10.1002/cpbi.5 Tang R. Luo M. Li J. Peng Y. Wang Y. Liu B. (2020). Relationships between Vasomotor Symptoms and Mood in Midlife Urban Chinese Women: Observations in a Prospective Study. J. Clin. Endocrinol. Metab. 105, 3437–3448. 10.1210/clinem/dgaa554 Tang R. Luo M. Li J. Peng Y. Wang Y. Liu B. (2019). Symptoms of Anxiety and Depression Among Chinese Women Transitioning through Menopause: Findings from a Prospective Community-Based Cohort Study. Fertil. Steril. 112, 1160–1171. 10.1016/j.fertnstert.2019.08.005 Willner P. Towell A. Sampson D. Sophokleous S. Muscat R. (1987). Reduction of Sucrose Preference by Chronic Unpredictable Mild Stress, and its Restoration by a Tricyclic Antidepressant. Psychopharmacology (Berl) 93, 358–364. 10.1007/BF00187257 Xia B. Tong Y. Xia C. Chen C. Shan X. (2020). α-Cyperone Confers Antidepressant-like Effects in Mice via Neuroplasticity Enhancement by SIRT3/ROS Mediated NLRP3 Inflammasome Deactivation. Front. Pharmacol. 11, 577062. 10.3389/fphar.2020.577062 Yi L. T. Li J. Geng D. Liu B. B. Fu Y. Tu J. Q. (2013). Essential Oil of Perilla Frutescens-Induced Change in Hippocampal Expression of Brain-Derived Neurotrophic Factor in Chronic Unpredictable Mild Stress in Mice. J. Ethnopharmacol. 147, 245–253. 10.1016/j.jep.2013.03.015 Yuan T. F. Wang S. T. Le J. Li Y. (2021). Steroid Profile Analysis by Liquid Chromatography-Tandem Mass Spectrometry in Second-Trimester Pregnant Women for Trisomy 21 Screening. J. Pharm. Biomed. Anal. 197, 113966. 10.1016/j.jpba.2021.113966 Zhang K. Wang Z. Pan X. Yang J. Wu C. (2020). Antidepressant-Like Effects of Xiaochaihutang in Perimenopausal Mice. J. Ethnopharmacol. 248, 112318. 10.1016/j.jep.2019.112318 Zhang R. Zhu X. Bai H. Ning K. (2019). Network Pharmacology Databases for Traditional Chinese Medicine: Review and Assessment. Front. Pharmacol. 10, 123. 10.3389/fphar.2019.00123 Zhao Y. Q. Teng J. Yang H. J. (2015). Analysis on Medication Regularity of Modern Traditional Chinese Medicines in Treating Melancholia Based on Data Mining Technology. China J. Chin. Mater. Med. 40, 2042–2046. 10.4268/cjcmm20151036 Zhou Q. Wang B. Hua Q. Jin Q. Xie J. Ma J. (2021). Investigation of the Relationship between Hot Flashes, Sweating and Sleep Quality in Perimenopausal and Postmenopausal Women: The Mediating Effect of Anxiety and Depression. BMC Womens Health 21, 293. 10.1186/s12905-021-01433-y Zhou X. D. Yang X. J. Zheng Y. Qin Z. S. Sha W. Chen G. (2020). Jie-Yu Pill, A Proprietary Herbal Medicine, Ameliorates Mood Disorder-like Behavior and Cognitive Impairment in Estrogen-Deprived Mice Exposed to Chronic Unpredictable Mild Stress: Implication for a Potential Therapy of Menopause Syndrome. Front. Psychiatry 11, 579995. 10.3389/fpsyt.2020.579995