Modular Characteristics and Mechanism of Action of Herbs for Endometriosis Treatment in Chinese Medicine: A Data Mining and Network Pharmacology–Based Identification

Endometriosis is a common benign disease in women of reproductive age. It has been defined as a disorder characterized by inflammation, compromised immunity, hormone dependence, and neuroangiogenesis. Unfortunately, the mechanisms of endometriosis have not yet been fully elucidated, and available treatment methods are currently limited. The discovery of new therapeutic drugs and improvements in existing treatment schemes remain the focus of research initiatives. Chinese medicine can improve the symptoms associated with endometriosis. Many Chinese herbal medicines could exert antiendometriosis effects via comprehensive interactions with multiple targets. However, these interactions have not been defined. This study used association rule mining and systems pharmacology to discover a method by which potential antiendometriosis herbs can be investigated. We analyzed various combinations and mechanisms of action of medicinal herbs to establish molecular networks showing interactions with multiple targets. The results showed that endometriosis treatment in Chinese medicine is mainly based on methods of supplementation with blood-activating herbs and strengthening qi. Furthermore, we used network pharmacology to analyze the main herbs that facilitate the decoding of multiscale mechanisms of the herbal compounds. We found that Chinese medicine could affect the development of endometriosis by regulating inflammation, immunity, angiogenesis, and other clusters of processes identified by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. The antiendometriosis effect of Chinese medicine occurs mainly through nervous system–associated pathways, such as the serotonergic synapse, the neurotrophin signaling pathway, and dopaminergic synapse, among others, to reduce pain. Chinese medicine could also regulate VEGF signaling, toll-like reporter signaling, NF-κB signaling, MAPK signaling, PI3K-Akt signaling, and the HIF-1 signaling pathway, among others. Synergies often exist in herb pairs and herbal prescriptions. In conclusion, we identified some important targets, target pairs, and regulatory networks, using bioinformatics and data mining. The combination of data mining and network pharmacology may offer an efficient method for drug discovery and development from herbal medicines.

Introduction

Endometriosis is a common benign condition in women of reproductive age, which is characterized by inflammation, compromised immunity, hormone dependence, and neuroangiogenesis, etc. (Rogers et al., 2009). The main symptoms of endometriosis include dysmenorrhea, endometriosis-associated pain, and infertility (Brown and Farquhar, 2014; Bedaiwy et al., 2017). In addition to surgical treatment, conventional drug therapies includes nonsteroidal anti-inflammatory drugs (NSAIDs) and other drugs, such as androgens, aromatase inhibitors, selective progesterone receptor modulators, oral contraceptives, and gonadotropin releasing hormone (GnRH) agonists (Mihalyi et al., 2006; Dunselman et al., 2014; Legendre et al., 2018). Although these drugs play a certain therapeutic role, the main effect is only to relieve symptoms, delay recurrence, and bring certain side effects. Chinese medicine has attracted extensive attention for its ability to treat complex diseases due to its moderate treatment effect and lower side effect. Given the complex pathogenesis of the condition and limited treatment effect, patients with endometriosis often explore Chinese medicines to treat the primary lesions or control the symptoms (Flower et al., 2012). Chinese medicine has significant advantages in treating gynecological disorders such as endometriosis associated dysmenorrhea, chronic pelvic pain, abnormal uterine bleeding, and infertility (Zhou et al., 2019). However, due to the differences in clinical experience and the complexity of Chinese medicine, prescriptions differ greatly from each other. In addition, current Traditional Chinese Medicine (TCM)–related research is mainly focused on experience summarization (Su et al., 2014). These lack the in-depth study of the rules and mechanisms of medication in prescriptions. We urgently need to summarize the rules for the administration of Chinese medicine and analyze the mechanism of action of the herbs, in order to provide the basis for the optimization of drug use and the screening of new compounds.

Data mining approaches and network pharmacology have been utilized to discover the underlying pathogenesis of the condition (Ren et al., 2019). Data mining approaches can analysis the frequency of herbal medications, formula ontology, and changes in formula patterns derived from a knowledge graph (Leem et al., 2018; You et al., 2019). Network pharmacology is a holistic approach to understanding the function and behavior of a biological system, at the systems level in the context of biological networks. Network pharmacology, when combined with traditional pharmacology, can predict the target profiles and pharmacological actions of herbal compounds to reveal drug-gene-disease comodule associations, screen synergistic multicompounds, and elucidate herbal ingredients with their related properties and compound-target and target-disease relationships (Zeng et al., 2017).

In this study, data mining was used to analyze the application rules of TCM in the treatment of endometriosis. We also proposed a computational systems pharmacology method and statistical analysis to determine the associated molecular mechanisms. We analyzed an endometriosis biomarker network and identified biologically relevant functional modules of Chinese medicine targets in endometriosis. These findings would provide a better understanding of the mechanism of multitarget regulation of endometriosis in Chinese medicine and shed light on the screening of core compounds. This in turn has facilitated the discovery of novel effective drugs.

Materials and Methods

Workflow of Data Mining and Network Pharmacology Approach

The research framework of this study is presented in Figure 1A. It is summarized as follows: (1) Base datasets of herbal medicines with antiendometriosis effects were constructed. These datasets of Chinese herbal medicines were screened from the literature. (2) The main ingredients and targets were obtained. The targets of endometriosis treatment were derived from the Genecard dataset, National Center for Biotechnology Information (NCBI), Genbank, and the Online Mendelian Inheritance in Man (OMIM) databases, disease-related drug targets from the Therapeutic Target Database (TTD), DrugBank databases, the search tool for interactions of chemicals (STITCH) dataset (http://stitch.embl.de/) (Szklarczyk et al., 2016), etc. The Gene Ontology (GO) database, Kyoto Encyclopedia of Genes and Genomes (KEGG), and the Database for Annotation, Visualization, and Integrated Discovery (DAVID) were all explored. (3) Molecular networks and pathways were analyzed, mutual target proteins were mined and the regulatory network was explored. (4) The network of targets for the condition was constructed. Cytoscape, an open source software project for the integration of biomolecular interaction networks with high-throughput expression data, was applied and other molecular states were combined in a unified conceptual framework (Shannon et al., 2003).

FIGURE 1

Figure 1 Workflow of the systems pharmacology approach. (A) Method (1) Endometriosis-associated genes were integrated from three different sources: the Genecard database, PubMed database, Genbank, and the Online Mendelian Inheritance in Man (OMIM) database. (2) Construction of protein-protein interaction (PPI) topological modules. Gene Ontology (GO) and pathway enrichment functional analyses were conducted for the endometriosis and Chinese medicine targets modules. (3) Disease-related conventional medication targets from the Therapeutic Target Database (TTD) and DrugBank databases, guidelines, and published literature. (4) KEGG and GO analyses of the targets of major compounds. (5) Analysis of the targets of endometriosis and pathways that Chinese herbs could regulate. (B) Workflow of literatures search for data mining.

Search Methods

We conducted an electronic search of the Chinese and English databases from their inception to 31/12/2018. We searched the Chinese literatures in China National Knowledge Infrastructure (CNKI) database and WanFang database. The English databases included Embase and the PubMed database. All types and formulations of herbal medicines were considered, including extracts, decoctions, pills, and natural compounds. The outcomes of clinical studies were not considered.

The following search terms were used either individually or in combination: “Chinese medicine,” “traditional medicine,” “herbal medicine,” “oriental medicine,” “herb,” “plant,” “prescription,” “decoction,” and “endometriosis,” among others. All prescriptions for endometriosis treatment were selected and the drug records were complete. Formulations of pills, powders, and decoctions were all evaluated. The name of each herbal formula was extracted, as well as its composition of medicinal herbs, its origin (name of article or ancient literature), author, publication year, and internal/external application. Exclusion criteria included the following: missing prescription; the application of acupuncture, massage, or other external treatments. The method of literature searching is presented in Figure 1B.

TCM and Compound Targets

The names of Chinese herb included the Chinese pinyin and full plant name with taxonomic validation (species name/scientific name). The validated information of herbs/plants names with taxonomic validation were collected from Kew Royal Botanic Garden (https://mpns.science.kew.org/mpns-portal/), the Plant List (http://www.theplantlist.org), and the Plants of the World Online(http://plantsoftheworldonline.org/). Chemical ingredients were compiled from the TCM Systems Pharmacology (TCMSP) Database (http://tcmspw.com/tcmsp.php), TCM Integrative Database (TCMID) (http://www.megabionet.org/tcmid/), universal natural product database (UNPD) (Gu et al., 2013), TCM-MeSH(http://mesh.tcm.microbioinformatics.org/), STITCH database (http://stitch.embl.de), and DrugBank (https://www.drugbank.ca). Finally, the structures of relevant natural compounds in the medicinal plants were retrieved from the NCBI Pubchem dataset (https://www.ncbi.nlm.nih.gov/pccompound/).

Active Ingredients Screening

The absorption, distribution, metabolism, excretion, and toxicity (ADMET) modeling as a tool for rational drug design has significant effect in new drug discovery (Wang et al., 2015). The druggability of each candidate was analyzed according to its oral bioavailability (OB) and drug-likeness (DL) indices, as recommended by the TCMSP database. High OB values is often an important consideration for the development of bioactive molecules as therapeutic agents. DL evaluation is used in drug design to evaluate whether a compound is chemically suitable for use as a drug, and how drug-like a molecule is with respect to parameters that affect its pharmacodynamic and pharmacokinetic profiles, which ultimately impact its absorption, distribution, metabolism, and excretion properties. The molecules with OB ≥ 25% and DL ≥ 0.18 could be considered to exhibit relatively better pharmacological properties. We further selected the major compounds based on literatures to identify the potential therapeutic effect. Although some compounds, such as essential oils, have lower DL values, they are selected because the effect had been experimental verification.

Endometriosis-Related Targets and Endometriosis-Associated Drugs Targets

We also collected targets of endometriosis. To obtain the most reliable results, the following steps were taken: (A) endometriosis-associated genes were integrated from three different sources: the Genecard database, Malacards database (https://www.malacards.org/), literature on PubMed, Genbank, and the OMIM database (http://www.omim.org/). (B) Disease-related drug targets were compiled from the TTD (http://bidd.nus.edu.sg/group/cjttd/), DrugBank databases (http://www.drugbank.ca/, version: 3.0), PharmGKB (https://www.pharmgkb.org), The Drug-Gene Interaction database (DGIdb) (http://dgidb.org/search_interactions), as well as the European Society of Human Reproduction and Embryology (ESHRE) guidelines (Dunselman et al., 2014). We used those drug-target interactions to determine the treatments for which the targets were human genes/proteins. (C) Associated targets and structures of human proteins were determined. The targets and proteins were researched in the UniProt (https://www.uniprot.org).

Data Analysis

Data Mining Analysis

We used “A Framework for Automated Knowledge Graph Construction Towards Traditional Chinese Medicine” developed by the Second Affiliated Hospital of Guangzhou University of Chinese Medicine to analyze the Chinese herbs for endometriosis in the literature (Weng et al., 2017). In the early stage of data set, we used more than 1,100 ancient literatures in the “Chinese Medical Code” dataset to extract knowledge of Chinese medicine. And we further constructed a large-scale knowledge graph of TCM. It contains 13 first-class knowledge classifications, 116 second-class knowledge classifications, more than 59,000 core knowledge elements, more than 0.5 million expand core knowledge elements, and more than 215,000 knowledge lists of herbal medicine and prescriptions. Then we used the recurrent neural network (RNN) model (Cho et al., 2014) to establish the generalized Chinese medicine knowledge meta-semantic representation. And we used fine-tune's transfer learning method, a way of machine learning, to learn the core knowledge elements of gynecological disease in dataset, and further obtained Chinese medicine knowledge meta-semantic representation vector. Finally, the hierarchical clustering results network and semantic representation could be visualized. Based on the above research, we selected relevant prescriptions in endometriosis treatment for analysis.

A. Descriptive analysis and frequent itemset mining: the frequency and pattern was used to calculate the type and frequency of each Chinese medicinal compound (Chen et al., 2008). The TCM paired drugs discovered by frequent itemset analysis (Feng et al., 2006).

B. Knowledge graph and RNNs: this was conducted using a hierarchical clustering algorithm, according to the nature, flavor, meridian tropism, and main efficacy of medicinal compounds. We established the knowledge graph of RNNs. RNNs are temporal-based neural networks that capture long sequences of inputs using the internal memory which has been used in medical research (Zhang Y. et al., 2019).

C. Association rule mining (ARM): The compatibility rules of the couplet medicinal/group medicinal were extracted according to association rule analysis. The strength of the association rule was measured in terms of its support, confidence, and lift. ARM was used to demonstrate the prescription (Chen et al., 2008).

ARM and RNNs analysis could further reflect the compatibility rules. Based on the above, combined with the results of previous literature reviews, further screening of TCMs, drug pairs, and compound pharmacological analyses could be performed.

Construction of Gene Enrichment Analysis

To better understand the processes associated with endometriosis and Chinese medicine treatment, we performed GO terms and KEGG pathway enrichment analyses. The GO enrichment analysis provides three structured networks of defined terms to describe gene attributes. Enriched GO terms are classified according to biological process (BP), molecular function (MF), and cellular component (CC). The KEGG (http://www.genome.jp/kegg/) is a database for large-scale systematic analysis of molecular interaction networks of genes or proteins (Kanehisa and Goto, 2000). The DAVID bioinformatics resources consist of an integrated biological knowledge base and analytic tools, aimed at systematically extracting biological meaning from large gene or protein lists. We used the web-based search engine, DAVID, to determine over-represented GO terms and KEGG pathways with thresholds of an enrichment score > 2, count > 5, and P < 0.05 and analyzed endometriosis-related pathways and GO terms. Venn diagram and bubble graphs of were performed using the OmicShare tools, a free online platform for data analysis (http://www.omicshare.com/tools). Significant pathway terms of KEGG were mapped into a bubble graph. The big and higher bubbles represent those highly significantly enriched pathway terms.

Protein-Protein Interaction Network Analysis

Protein-protein interaction (PPI) networks included information on the biological processes and molecular functions of cells (Vidal et al., 2011). We used the online search tool for recurring instances of neighboring genes (STRING, Version 9.1) (http://www.string-db.org) to predict the interactions. The Cytoscape software 3.7.2 (http://cytoscape.org/) was used to visualize networks. To identify crucial relationships in the PPI network, potentially overlapping modules that were densely connected were subsequently identified using Molecular Complex Detection (MCODE) plugin in the Cytoscape program. The MCODE plugin was used to re‐analyze the clusters among the network according to the k‐core = 2. A value of P < 0.01 was considered the significant threshold. And Cytoscape plugin cytohubba (Chin et al., 2014) was used to identify central elements of biological networks.

Network Construction of KEGG Pathway

To better analysis the holistic mechanism of Chinese herbs in endometriosis treatment, the subnetworks pathway was compiled by following the procedures: All targets of Chinese herbs in endometriosis treatment were submitted to an online tool KEGG Search Pathway (https://www.genome.jp/kegg/tool/map_pathway1.html). Based on the mechanism of endometriosis, multiple pathways were integrated and overlapped according to cross-talk targets in these maps. Based on the cross-talk of the pathways, we further constructed a targets-pathways network of Chinese medicine treatment.

The Connectiviy Map Analysis

The Connectivity Map (CMap) is an online pharmacogenomic database using gene-expression signatures from cultured cells treated individually with various chemicals to connect small molecules, genes, and disease(Chin et al., 2014). The expression profiles GSE25628 were downloaded from the Gene Expression Omnibus (GEO) database of patients both with and without endometriosis. The GSE25628 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE25628) gene expression profiles comprised 22 samples (GPL571 Affymetrix Human Genome U133A 2.0 Array), with eight ectopic samples and eight eutopic endometrium samples from several affected women in the proliferative phase, and normal health donors as the control. The limma package in the R software was used to compare the differentially expressed genes (DEGs) between eutopic endometrium samples group and ectopic endometrium samples group, and the DEGs between eutopic endometrium samples group and endometrium in normal health women group. In addition, these two sets of DEG profiles were used to predict compounds by the Connectivity Map O2 (CMap) (https://portals.broadinstitute.org/cmap/, Update September 12, 2017) (Lamb et al., 2006; Lamb, 2007). The representing similarity was calculated, ranging from −1 to 1. A positive connectivity score indicated that the drug was able to induce the input signature in human cell lines. The negative connectivity score was indicated potential therapeutic value. After rank ordering of all instances, the connectivity score of various instances were filtered according to the P-value < (0.05) to determine the natural compounds for further consideration.

Network Construction

We established a network analysis of the relevant compounds targets. The components of the target networks for the herbs were constructed using the Cytoscape software 3.7.2. The nodes in each network were evaluated based on three indices: degree, node betweenness, and node closeness. Degree indicated the number of edges between a single node and other nodes in a network. Node closeness represented the inverse of the sum of the distance from one node to other nodes. The importance of a node in a network was indicated by the values of these indices, with higher values indicating greater importance.

Results

Chinese Medicines Commonly Used for Endometriosis

In the past, we have studied the prescriptions of Chinese medicine for endometriosis through literature monographs. A total of 551 literatures were screened: 513 literatures in Chinese databases (including CNKI and WanFang databases) and 38 literatures from the PubMed and Embase database. Among those screened, 315 literatures were clinical studies, and 234 literatures were experimental studies (25 literatures were in vitro, 210 literatures were in vivo, and one literature was both in vitro and in vivo). The major TCM literatures were clinical observational studies. Fifteen literatures were clinical randomized controlled studies. Some literatures contain two or three prescriptions according to different syndromes in TCM. A total of 615 medicinal prescriptions were selected from the literature screened (Figure 1B). The most published prescriptions included those for Guizhi Fuling Wan (GFW), Jiawei Sanleng Wan, Neiyi Fang, and Shaofu Zuyu decoction.

The total 231 kinds of Chinese herbs were selected in literatures. The top 20 herbs were determined using frequency statistics and are presented in Table 1. The validated information of major herbs, including location, used part, famliy, genus are showed in Supplementary Table 1. We analyzed the classification, medicinal properties, and medicinal taste of all common drugs for endometriosis. Analysis from the prescriptions for endometriosis indicated a higher proportion and frequency of deficiency-tonifying herbs, blood-activating herbs, dampness-draining diuretic herbs, heat-clearing herbs, and digestant herbs (Figure 2A). Curcuma phaeocaulis Valeton (Ezhu), Sparganium stoloniferum (Buch.-Ham. ex Graebn.) Buch.-Ham. ex Juz.[Typhaceae] (Sanleng), and Corydalis yanhusuo (Y.H.Chou & Chun C.Hsu) W.T.Wang ex Z.Y.Su & C.Y.Wu(Yanhusuo) are both effective in promoting blood circulation to alleviate pain and activate qi to resolve stagnation in Chinese theory. The treatment of endometriosis with basic theory of Chinese medicine has always been based on promoting blood circulation to minimize blood stasis and alleviate pain.

TABLE 1

Table 1 Frequency, module, and meridian tropism of the top 20 medicinal herbs.

FIGURE 2

Figure 2 Data mining of Chinese herbs from literature. (A) Ratio of herbs of different classifications. Traditional Chinese medicine includes tonics, heat-clearing medicines, and blood-activating and stasis-removing medicines. (B–D) Ratios of the taste and action of Chinese herbs. Among the Chinese herbs evaluated, most were used as warm drugs, mainly targeting the liver and spleen. (E, F) Knowledge Graphs of major Chinese herbs. The size of a node represents the frequency of Traditional Chinese medicine in our gynecological dataset of Chinese medicine. The node represents the knowledge element, Chinese medicine. Cosine distance of the knowledge meta-semantic representation vector, a feature of clustering, was evaluated. Based on the association rules, frequent itemsets and recurrent neural networks (RNNs), the knowledge graphs were constructed, which were divided into four clusters and nine clusters, respectively. The knowledge graphs may provide the basis for the discovery of relevant prescriptions. Cluster analysis could also provide a reference for prescription guidelines.

The characteristics of four properties, five tastes, and channel tropism were evaluated. The analysis of herbal properties is presented in Figures 2B–D. Chinese herbs that were considered warm and slight were given priority. The treatment of endometriosis is mainly based on invigorating the spleen, tonifying the liver to strengthen healthy qi, and supplementation with blood-activating herbs to eliminate pathogenic factors.

Frequent itemset mining is a crucial data mining task with numerous applications in discovery, including recommendation and classification, among others (Soltani and Akbarzadeh-T, 2014). By considering both the frequent itemsets and association rules, we inferred the relevant drug pairs and combinations, such as S. stoloniferum (Sanleng)-C. phaeocaulis(Ezhu); C. phaeocaulis-Paeonia lactiflora Pall(Chishao); and C. yanhusuo (Yanhusuo)-Angelica sinensis (Oliv.)Diels (Danggui), among others. The main roles of these herbs were to promote blood circulation and relieve pain. Table 2 shows the herbs commonly used in combination, as determined by frequent itemset mining analysis. The frequent item sets were showed in Supplementary Table 2.

TABLE 2

Table 2 The major Chinese herbs pairs/prescriptions in endometriosis treatment.

To uncover latent knowledge, unsupervised machine learning and RNNs were applied to generate the knowledge graph, by calculating the distance between knowledge element vectors, according to preset categories. We set various parameters, and obtained multiple knowledge graphs. We obtained a ratio of nine clusters/four clusters of knowledge maps according to the clinical medication rules of Chinese medicine. Different clusters were formed by setting parameters and applying the combination rules of TCM. These clusters may reflect the law of compatibility and combination of clinical prescriptions in TCM. Figures 2E, F demonstrate the knowledge graphs for endometriosis herbs. The classifications reflect the possible combinations of TCM prescription medications. Based on the knowledge graphs and association rule mining, synergistic herbal combinations could be derived; however, further analysis of the associated mechanism is required.

GO terms and KEGG Pathway Analysis of Endometriosis-Associated Genes

We collected a total of 1,289 endometriosis-related targets from Genecard, Genbank, as well as the OMIM records of genomic databases of human diseases. The Venn diagram in Figure 3A reflects common genes from different databases. The genes significantly associated with endometriosis were tested for functional enrichment, including the relevant pathways and GO terms. The results of the GO analysis indicated that “inflammatory response,” “innate immune response,” and “chemokine production” were the most significant terms related to endometriosis in the BP category. Furthermore, the significant MF terms related to endometriosis were “cytokine activity,” “growth factor activity,” and “receptor binding.” The results of GO enrichement showed in Figure 3B.

FIGURE 3

Figure 3 Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analysis of endometriosis-associated genes. (A) Venn diagram showing analysis of endometriosis genes from different datasets. MMP2, PTEN, LHB, BIRC5, EMX2, ENDO1, KRAS, NR5A1, AHRR, CCL11, MMP14, CCL5, ESR2, ARID1A, and MMP9 were the shared genes from the Genecard, Genbank, and OMIM datasets. (B) GO enrichment of endometriosis-associated genes. (C) KEGG pathway of endometriosis genes. Bubble chart for significantly enriched pathway terms from KEGG analysis. The size of each bubble is reflective of the relative ratio of the number of targets hitting each pathway over the total number of targets. P indicates the statistical significance of the P-values; greater numbers indicate greater significance. (D) Hub genes of endometriosis-associated genes, determined by cytoHubba in Cytoscape 3.7.2. (E) Cluster analysis of endometriosis-associated genes, determined by Molecular Complex Detection (MCODE) in Cytoscape 3.7.2.

The majority of the pathways were related to cytokine-cytokine receptor interactions, the PI3K-Akt signaling pathway, focal adhesion, the TNF signaling pathway, and HIF-1 signaling pathway, among others (Figure 3C). The KEGG pathways of endometriosis associated genes are showed in Supplementary Table 3. These signal pathways can affect the cell proliferation, migration and invasion of endometriosis.

We used the cytoHubba in Cytoscape to find the hub genes. The hub genes showed in Figure 3D. And we further used MCODE plugin densely connected regions to analyze the major modules in biological processes. The results showed that the genes associated with endometriosis could be divided into 11 clusters (Figure 3E). The PPI network of the major endometriosis genes is showed in Supplementary Figure 1. The genes in cluster 1, which is dominated by chemokines, are related to the chemokine signaling pathway and cytokine-cytokine receptor interactions. The genes in cluster 2 participate in receptor ligand activity, growth factor activity, HIF-1 signaling, PI3K-Akt signaling, and the Ras signaling pathway. The genes in cluster 3 are mainly interleukins, suggesting pathways related to inflammation. In addition, cluster 4 is related to kinase proteins. Cluster 7 is related to hormone regulation. Cluster 9 is closely related to tumor transcription regulation. Cluster 6 and cluster 10 are related to cell adhesion, intercellular adhesion molecules, vascular cell adhesion molecules, and matrix metalloproteinases. The pathogenesis of endometriosis is related to chemokine regulation, adhesion, invasion, angiogenesis, inflammation, immunity, and hormone regulation. These findings provide the genetic background of endometriosis for further network pharmacology research.

Overall Ingredients and Targets of Common Chinese Herbal Treatments

Chinese herbal medicines each contain dozens or even hundreds of various ingredients. Thus, the compilation of an ingredient database is crucial. Some Chinese medicines are animal-derived compounds that were not included in the network pharmacological analysis. Typha angustifolia L. [Typhaceae] (Puhang) and Faeces Togopteri (Wulingzhi) represented one herbal pair (HP). The source of F. Togopteri is Trogopterus xanthipes feces; thus, we did not analysis this HP. The single herbs, HPs, and prescriptions were evaluated, to determine the characteristics of Chinese medicine more comprehensively. In order to further corroborate our predictions, a literature review was conducted to determine whether the compounds were already experimentally validated for any associated therapeutic effects. We analyzed gene enrichment for the major Chinese herbs, as hub herb associated targets. Compounds and targets from of major Chinese herbs for endometriosis treatment are listed in Supplementary Table 4. The Chinese herbs associated targets are listed in Supplementary Table 5.

Salvia miltiorrhiza Bunge[Lamiaceae] (Dangshen)

More than 200 compounds from Salvia miltiorrhiza Bunge were found, including lipophilic diterpenoids, water-soluble phenolic acids, and other constituents. The active ingredients of S. miltiorrhiza include tanshinone I, tanshinone IIA, salvianolic acid, and dihydrotanshinone, among others (MEIm et al., 2019). A total of 196 targets were possibly related to S. miltiorrhiza, and 63 targets could be associated with endometriosis genes. Based on the KEGG pathway analysis, the main pathways were those associated with cancer, calcium signaling, VEGF signaling, T cell receptor signaling, progesterone-mediated oocyte maturation, apoptosis, and p53 signaling pathway, among others that were related to endometriosis (P < 0.05). Compounds such as ursolic acid, rosmarinic acid, ferulic acid, caffeic acid, tanshinone IIA, protocatechuic acid, and tetramethylpyrazine could evidently regulate pain-associated targets to reduce pain, based on the network pharmacological analysis.

Angelica Sinensis (Oliv.)/Angelica Sinensis Var. Wilsonii (H.Wolff) Z.H.Pan & M.F.Watson (Danggui)

Angelica sinensis var. wilsonii (H.Wolff) Z.H.Pan & M.F.Watson is predominantly known as a treatment for intractable gynecological disorders. Ultrahigh performance liquid chromatography-tandem mass spectrometry (UHPLC-MS/MS) showed that A. sinensis contains eight components, including ferulic acid, senkyunolide A, butylphthalide, ligustilide, butylidenephalide, senkyunolide I, senkyunolide H, and levistolide A (Gui and Zheng, 2019). The volatile oil of A. sinensis, has evident antiinflammatory activities (Zhong et al., 2016). A total of 86 targets were possibly related to A. Sinensis, and 27 targets were associated with endometriosis genes. The KEGG enrichment showed associated signaling pathways in neuroactive ligand-receptor interactions, calcium signaling, TNF signaling, cGMP-PKG signaling, and the estrogen signaling pathway.

Corydalis yanhusuo (Y.H.Chou & Chun C.Hsu) (Yanhusuo)

Corydalis yanhusuo (Y.H.Chou & Chun C.Hsu) could effectively attenuate acute inflammatory and neuropathic pain. The main alkaloid contents and composition of C. yanhusuo includes protopine, α-allocryptopine, tetrahydrocolumbamine, coptisine, palmatine, berberine, dehydrocorydaline D, L-tetrahydropalmatine, tetrahydroberberine, corydaline, and tetrahydrocoptisine, as determined by the high performance liquid chromatography-diode array detector (HPLC-DAD) method. A total of 165 targets were possibly related to C. yanhusuo, and 47 targets were associated with endometriosis genes. The KEGG enrichment analysis revealed pathways associated with cancer, calcium signaling, neuroactive ligand-receptor interactions, VEGF signaling, apoptosis, T cell receptor signaling, and B cell receptor signaling pathway, among others.

Ligusticum chuanxiong S.H.Qiu, Y.Q.Zeng, K.Y.Pan, Y.C.Tang & J.M.Xu/Conioselinum anthriscoides ‘Chuanxiong’ (Chuanxiong)

Conioselinum anthriscoides ‘Chuanxiong’ is used to regulate menstruation and relieve pain in multiple diseases. A total of 174 ingredients of L. chuanxiong and 22 compounds demonstrated favorable bioavailability (Chen Z. et al., 2018). A total of 95 targets may be related to L. chuanxiong, and 26 targets were associated with endometriosis genes. Based on the KEGG enrichment analysis, we obtained the following associated pathways: neuroactive ligand-receptor interaction, calcium signaling, pathways in cancer, and VEGF signaling pathway. We inferred that fumarine, isocorypalmine, and fagarine I could regulate the 5-hydroxytryptamine receptor to reduce pain in endometriosis.

Astragalus mongholicus Bunge (Huangqi)

Astragalus mongholicus Bunge is one of the most popular traditional medicinal herbs with several pharmacological activities, including hematopoietic, antiinflammatory, and immunological. The main components of A. mongholicus were hederagenin, kumatakenin, isorhamnetin, 3,9-di-O-methylnissolin, calycosin, 7-O-methylisomucronulatol, formononetin, quercetin, and betulinic acid (Fu et al., 2014). A total of 263 targets were possibly related to A. membranaceus, and 68 targets were associated with endometriosis genes. The KEGG enrichment analysis showed the following pathways: cancer, neuroactive ligand-receptor interactions, calcium signaling, endometrial cancer, p53 signaling, T cell receptor signaling, metabolism of xenobiotics by cytochrome P450, toll-like receptor signaling, and VEGF signaling pathway, among others.

Cyperus rotundus L. (Xiangfu)

Cyperus rotundus L., a widely distributed perennial sedge has a relatively higher concentration of active ingredients in the form of essential oils, phenolic acids, ascorbic acids, and flavonoids in the tuber and rhizomes. C. rotundus is widely used in many disorders such as inflammation, diabetes, diarrhea, tumors, among others (Pirzada et al., 2015). A total of 246 targets were possibly related to C. rotundus, and 74 targets were associated with endometriosis genes. The KEGG enrichment analysis showed the following pathways: PI3K-Akt signaling pathway, neuroactive ligand-receptor interaction, TNF signaling pathway, MAPK signaling pathway, and HIF-1 signaling pathway, among others.

Commiphora myrrha (T.Nees) Engl. (Moyao)

Commiphora myrrha (T.Nees) Engl., the Commiphora species, known as “myrrh,” are characterized by resinous exudates from the bark of plants. They are used in the treatment of trauma, arthritis, and fractures, and exert antiproliferative, antioxidant, antiinflammatory, and antibacterial effects (Shen et al., 2012). The compounds of C. myrrha include terpenoids, steroids, flavonoids, sugars, and lignans, among others. Furanosesquiterpenes, such as furanoelemanes, furanoeudesmanes, and furanogermacranes, are compounds with analgesic effects. A total of 213 targets were possibly related to C. myrrha, and 104 targets were associated with endometriosis genes. The KEGG enrichment analysis showed the following pathways: TNF signaling, HIF-1 signaling, toll-like receptor signaling, and PI3K-Akt signaling, among others.

Carthamus tinctorius L. (Honghua)

Carthamus tinctorius L. can invigorate blood circulation and has been recently shown to have antioxidant, analgesic, antiinflammatory, and antidiabetic properties. Carthamidin, carthamone, carvacrol, and isocarthamidin luteolin are the main constituents of C. tinctorius L. (Asgarpanah and Kazemivash, 2013). Furthermore, hydroxyethylcarthamin, hydroxysafflor yellow A (HSYA), safflor yellow B (SYB), safflomin A, safflomin B, safflomin C, isosafflomin C, safflor yellow A (SYA), and precarthamin, among others, have also been reported as constituents (Yue et al., 2013). A total of 74 targets were possibly related to C. tinctorius, and 37 targets were associated with endometriosis genes. The KEGG enrichment analysis showed the following significant pathways: cancer, FoxO signaling, endometrial cancer, and drug metabolism - cytochrome P450, among others.

HPs for Endometriosis Treatment

Curcuma phaeocaulis (Ezhu) and Sparganium Stoloniferum (Buch.-Ham. ex Graebn.) Buch.-Ham. ex Juz. (Sanleng)

Curcuma phaeocaulis Valeton is the most widely used species in HPs for endometriosis in Chinese medicine, and its essential oils are widely applied in the treatment of tumors in China. In C. Rhizoma, curcumol, bisdemethoxycurcumin, (4S,5S)-germacrone-1,4-diepoxide, aromadendrene, hederagenin, epoxycaryophyllene, and calarene showed relatively higher levels of OB and DL (Zhou Y. et al., 2016). S. stoloniferum is used for its hematopoietic functions, antiinflammatory activity, and immunological properties. Notable constituents include trans-gondoic acid, hederagenin, beta-sitosterol, formononetin, stigmasterol, and epibetulinic acid. A total of 130 targets were possibly related to C. phaeocaulis and S.stoloniferum. The GO enrichment analysis showed that C. phaeocaulis and S.stoloniferum could regulate apoptosis, cell death, and cell proliferation. The KEGG analysis showed significant pathways in cancer, neuroactive ligand-receptor interaction, calcium signaling, the cGMP-PKG signaling pathway, PI3K-Akt signaling pathway, apoptosis, and serotonergic synapse, among others. A total of 38 targets could be associated with endometriosis genes. Topological analysis showed the major targets included: TP53, SRC, TNF, VEGFA, PIK3CA, IL8, and EGFR, among others, which are also the core genes that cause endometriosis.

Core Prescriptions-GFW

GFW was the core prescription used in data mining for endometriosis treatment, as well as the classical prescription in Chinese medicine. A total of 565 herbal compounds could be identified from databases. There were 230 targets in Cinnamomum cassia (L.) J.Presl [Lauraceae](Guizhi); 78 targets in Paeonia suffruticosa Andrews [Paeoniaceae](Mudanpi); 135 targets in Paeonia lactiflora Pall.(Chishao); 67 targets in Prunus persica (L.) Batsch(Taoren); and 55 targets in Poria cocos(Schw.)Wolf(Fuling). In Cinnamomum cassia, taxifolin; beta-sitosterol; sitosterol; catechin; ent-epicatechin; and peroxyergosterol all showed effective OB and DL levels. Protocatechuic acid, coumarin, cinnamyl alcohol, 2-methoxy cinnamic acid, cinnamic acid, and cinnamaldehyde are considered the core components of C. cassia, and the core herbs in such prescriptions. The target proteins in C. cassia were focused on neurological disease, inflammatory disease, cancer, cellular growth and proliferation, cell signaling, and molecular transport. P. lactiflora is used for pain and blood stasis, and has hematopoietic functions, antiinflammatory activity, and immunological properties. And a total of 94 of its targets were associated with endometriosis treatment. Furthermore, 70 targets in P. suffruticosa, 16 targets in P. cocos and 10 targets in P. persica could be associated with endometriosis genes.

In the present study, we obtained 521 targets (Supplementary Table 5), which may be regulated by the above mentioned Chinese herbs derived from databases. These targets were involved in kinase pathways, angiogenesis, inflammation, immunity, and other modules. The main pathways included those of neuroactive ligand-receptor interactions, toll-like receptor signaling, metabolism of xenobiotics by cytochrome P450, VEGF signaling, apoptosis, drug metabolism, endometrial cancer, and the calcium signaling pathway, among others.

The venn diagram of the common targets between Chinese herbs and conventional treatment drugs is showed in Figure 4A. The results of GO enrichment are presented in Figure 4B. The KEGG pathway annotation and KEGG pathways enrichment are presented in Figures 4C, D. The major KEGG pathways of major Chinese herbs, herbs pairs and prescription are showed in Supplementary Table 6. Furthermore, following cytoHubba analysis, we found that VEGFA, MAPK3, JUN, AKT1, TP53, IL6, ALB, INS, MAPK1, MMP9, PTGS2, among others could be core genes in the regulation of the screened herbs. Using MCODE analysis, we determined the modules of various herbs (Figure 4E). The results of the modules were related to the pathogenesis of endometriosis.

FIGURE 4

Figure 4 Multiregulation effect of all targets in selected Chinese medicines. These traditional Chinese medicines are mainly therapeutic medicines and official prescription medicines, which could regulate multiple pathways involving nerves, inflammation, and immunity. These actions reflect the multitarget regulatory mechanisms of Traditional Chinese medicine. (A) Venn diagram showing the common targets between Chinese herbs and conventional treatment drugs. (B) Gene Ontology (GO) enrichment analysis of Chinese herbs and associated targets. (C) Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in annotation. (D) KEGG pathway enrichment of the targets of Chinese herbs. (E) Module analysis of Chinese herbs, determined by Molecular Complex Detection (MCODE) of the Cytoscape software. These targets could be divided into 13 clusters. The MCODE_2 (pathways in cancer and response to reactive oxygen species); MCODE_3 (cytochrome P450 - arranged by substrate type); MCODE_6 (cellular responses to stress), MCODE_7 and MCODE_12 (apoptosis); MCODE_9 (adenylate cyclase-activating G protein-coupled receptor signaling pathway, cAMP-mediated signaling pathway); and MCODE_11 (VEGFR2 mediated cell proliferation, Ras signaling pathway) were all related to endometriosis pathogenesis.

Overlap of Chinese Herbs and Endometriosis Disease Modules

We further analyzed the targets of endometriosis-associated genes that coincided with those of Chinese herbs. A total of 170 targets of Chinese herbal medicines coincided with those of genes associated with endometriosis. The targets of endometriosis-associated genes that coincided with those of Chinese herbs are showed in Supplementary Table 7. Chinese herbal treatment for endometriosis could regulate the biological processes of “regulation of apoptosis,” “regulation of cell death,” and “response to hormone stimulus,” among others. Figure 5A shows the GO enrichment of herb-associated targets. The GO enrichment of endometriosis-associated genes that coincided with Chinese herbs are showed in Supplementary Table 8.

FIGURE 5

Figure 5 Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses of the targets of Chinese medicine in endometriosis treatment. (A) GO analysis of the targets of endometriosis-associated genes that coincided with those of Chinese herbs. (B) KEGG enrichment of endometriosis-associated genes that coincided with those of Chinese herbs. (C) Pathways and targets network of Chinese medicine in endometriosis treatment. (D) Module analysis of the targets of Chinese medicine in endometriosis treatment, determined by Molecular Complex Detection (MCODE) of the Cytoscape software.

We also identified associated pathways, such as those of cancer, PI3K-Akt signaling, MAPK signaling, FoxO signaling, focal adhesion, HIF-1 signaling, Ras signaling, TNF signaling, estrogen signaling, toll-like receptor signaling, and VEGF signaling pathway, among others. To determine the functions of Chinese herbs, significant pathway terms in KEGG analysis were mapped onto a bubble graph (Figure 5B). Larger and higher bubbles in the Figure represent the highly, significantly enriched pathway terms. The network of major pathways and targets is presented in Figure 5C. Table 3 shows the major KEGG pathways of Chinese herbs in endometriosis treatment.

TABLE 3

Table 3 The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of endometriosis-associated genes that coincided with Chinese herbs.

The MCODE plugin was used to analyze the modules. Cluster 1 was related to endocrine resistance, prolactin signaling pathway, and pathways in cancer. Cluster 2 was related to signaling by receptor tyrosine kinases and EGFR tyrosine kinase inhibitor resistance. Cluster 3 was related to responses to steroid hormones. Cluster 4 was related to the Pathway Interaction Database ceramide pathway. Cluster 5 was related to steroid hormone biosynthesis, cytochrome P450-arranged by substrate type, and estrogen metabolic process. Cluster 6 was related to the adenylate cyclase-activating G protein-coupled receptor signaling pathway. Cluster 7 was related to cellular detoxification in GO enrichment (Figure 5D). Module analysis also facilitated the discovery of potential sub-modules of TCM in the regulation of endometriosis.

Target-Pathway Network

We constructed a KEGG map of Chinese medicine associated targets to reflect the regulation of Chinese medicine on endometriosis genes. TCM affects the development of endometriosis through multiple interactions such as inflammation, immunity, angiogenesis, and kinase pathways. Multiple pathways were integrated and overlapped, based on cross-talk targets.The common targets showed the largest overlap with the cancer signaling pathway. The multiregulation map of the KEGG pathway is shown in Figure 6. By mapping the targets to related pathways, we found that endometriosis treatment is mostly related to four function modules of the mechanism pathways. The cross-talk pathways network of Chinese herbs in endometriosis treatment is showed in Supplementary Figure 2. Therefore, we gained a deeper understanding of the mechanisms of these pathways.

FIGURE 6

Figure 6 Pathway of multiregulatory effects in Chinese medicine. Target network showing multiregulation in Chinese medicine. The reconstructed Kyoto Encyclopedia of Genes and Genomes (KEGG) map reflects the multipath, multitarget regulatory pathway of Chinese medicine. Target names are shown in rectangles. The pathway network with a blue background is a pathway module that is directly related to endometriosis in Chinese medicine. The green background indicates the network of Chinese medicine that regulates endometriosis-associated pain and estrogen regulatory pathways. The red border indicates the direct target for Chinese medicine treatment in endometriosis. The black border indicates the indirect target. The arrows in the figure indicate the upstream and downstream genes in pathway actions of the targets. Interactions between the pathways were noted. Each pathway target synergistically regulates endometriosis pain, inflammation and immunity, adhesion, invasion, and angiogenesis.

Pain relieving function is the major module in Chinese treatment, which is related to the nervous system in KEGG categories, such as neuroactive ligand-receptor interaction, serotonergic synapse, GABAergic synapse, cholinergic synapse, neurotrophin signaling pathway, dopaminergic synapse, calcium signaling pathway, glutamatergic synapse, and cAMP signaling pathway, among others.

The major targets, acetylcholinesterase (AChE), adrenoceptor beta, HTR2A prostaglandin endoperoxide synthase (PTGS), brain-derived neurotrophic factor (BDNF), serotonin receptors (HTR2A, HTR3A, HTR1A, and HTR2C) are all related to the regulation of neural receptors. The regulatory mechanism of these targets can directly affect various neurotransmitter receptors and the synthesis of neurotransmitters to modulate pain. Chinese medicine could also indirectly affect pain through downstream immune inflammatory factors (IL6, IL10, and TNF), and immuno-inflammatory pathways, such as cytokine-cytokine receptor interactions, the TNF signaling, VEGF signaling, HIF-1 signaling, toll-like receptor signaling, and PI3K-Akt signaling pathway. Estrogen could also regulated of the central serotoninergic system in pain.

The second module was associated with cell growth and death, invasion, adhesion, and angiogenesis-related signaling pathways in endometriosis. Chinese medicine could regulate invasion and adhesion-related signaling pathways, induce apoptosis, and inhibit cell proliferation, through processes such as apoptosis, P53 signaling, PI3K-Akt signaling, Wnt signaling, ECM-receptor interaction, and focal adhesion, among others. The angiogenesis and the tumor-related functions also showed considerable overlap, which included HIF signaling and the VEGF signaling pathway.

The third module was associated with inflammation and the immune system. Related processes included cytokine-cytokine receptor interaction, toll-like receptor signaling, T cell receptor signaling, B cell receptor signaling, NF-κB signaling, and IL-17 signaling pathway, which are also important in regulating inflammation and immune responses.

Moreover, kinase signaling pathways in endometriosis could affect the proliferation and differentiation of endometriosis cells, which could be potential targets for non-hormonal therapeutics. The associated module could be related to the canonical IKKβ/NFκB pathway, MAPK pathways, the PI3K/AKT/mTOR pathway, and the AMPK signaling pathway (McKinnon et al., 2016). These results suggest that multiple targets could affect various pathways to regulate the pathological processes of endometriosis.

Target Comparison of Chinese Herbs/Natural Compounds and Conventional Drugs

In order to clarify the similarities and differences between the related targets of Chinese medicine and conventional drugs for endometriosis, we considered the targets of conventional therapeutic drugs for comparison. At present, the treatment options for endometriosis include NSAIDs and hormonal drugs (progestogens, dienogest, and GnRH agonists). A total of 85 protein targets of these drugs were identified from DrugBank and the TTD database of conventional drugs. The drug targets of endometriosis conventional treatment are listed in Supplementary Table 9. Conventional drugs for the treatment of endometriosis are mainly aimed at the nervous system and hormonal pathway. The core pathways of these targets included steroid hormone biosynthesis, retinol metabolism, drug metabolism-cytochrome P450, and pentose and glucuronate interconversions, among others.

A total of 17 Chinese medicine associated targets (AR, BCHE, BCL2, CYP19A1, CYP1A2, CYP2C8, CYP3A4, ESR1, NR1I2, NR3C2, PGR, PLAT, PPARA, PPARG, PTGS1, PTGS2, and THBD) coincided with the targets of conventional drugs for endometriosis treatment. The enrichment of overlapping genes was mainly focused on the regulation of hormones and pain relief. The mechanism of action of Chinese herbs was similar to that of current conventional treatments. Furthermore, Chinese medicine and its compounds provided more possibilities for multitarget therapy.

Regulatory Effects of TCM on Related Targets of Endometriosis Pain

Endometriosis-related pain has been the main focus of TCM research. Thus, we further analyzed the related pathways of pain, to provide a basis for the discovery of effective pain-relieving compounds. The analgesic effect of TCM is related to the regulation of neurotransmitters and related pain factors. The serotonergic synapse was the main network of neurotransmitter regulation, and its KEGG map is presented in Figure 7A. Figure 7B shows the network of targets network of pain associated pathways.

FIGURE 7

Figure 7 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and sub-network analysis of traditional Chinese medicine for the endometriosis pain module. The pain module in endometriosis showed multiple neurotransmitters, nervous system pathways, and direct involvement in pain regulation. (A) KEGG mapper of the serotonergic synapse. (B) The network of targets network of pain associated pathways. Pathway of multineurotransmitter and inflammation, and immunity regulation network in endometriosis-associated pain. Purple nodes represent pathways, and blue nodes represent targets. (C) The protein-protein interaction (PPI) network of pain associated genes in Chinese medicine treatment. (D) Network of major compounds and targets in endometriosis-associated pain treatment. The middle green node represents pain treatment-related compounds, and the surrounding nodes are composed of different colors, which represent the target points enriched, according to the pathways. And the yellow node represents pathways.

We also selected the targets of Chinese medicines for the treatment of endometriosis from the pain database. This database is a comprehensive network of contextualized PPIs specifically associated with pain, that has been created through study of the pain interactome (Jamieson et al., 2014). The results showed that a total of 61 genes were involved in pain prediction, which were also related to endometriosis. These included PTGES, PTGS2, PTGS1, BDNF, TNF, IL6, ESR1, IL10, MMP9, and MAPK, among others (Figure 7C). These compounds, which have similar targets to those of NSAIDs and opioid analgesics, may be considered active compounds in Chinese medicine for the treatment of various pain-related diseases in the future. Figure 7D shows the network of the “nervous system pathways-targets-compounds”.

We found that Chinese medicines with analgesic effects include polysaccharides, saponins, alkaloids, flavonoids, terpenoids and others. Compounds such as ursolic acid, rosmarinic acid, ferulic acid, caffeic acid, tanshinone IIA, and oleanolic acid, in S. miltiorrhiza could regulate pain targets. Interestingly, we also found that many of the compounds that regulate pain originate from volatile oils/essential oils molecules (de Cassia et al., 2017). Terpene compounds are the main components of volatile oils. Moreover, the combination of multiple compounds may further play a synergistic effect in Chinese medicine. These volatile oil compounds have 109 targets. The volatile oil compounds associated targets were showed in Table 4. We found that the DL value of volatile oils was low. Thus, the pharmaceutical activity of these compounds still needs further experimental verification.

TABLE 4

Table 4 Volatile oils/essential oils from core herbs to associated pain.

Related Small Molecule Drugs Screening by CMap Analysis

We used expression profile data of CMap analysis for drug discovery. By querying CMap, we screened compounds in likely drug targets showing a similar gene expression profile with the desired and chemopreventive conditions. The DEGs were downloaded from the GSE25628 expression file. This dataset contained three sets of sample data, and we chose the comparison of two sets of samples for analysis. The endometriosis eutopic groups were compared with the non-endometriosis eutopic group (healthy women group), to include 52 upregulated and 169 downregulated DEGs (P < 0.05 and |logFC| > 2). A total of 67 upregulated and four downregulated DEGs with ectopic endometrium and eutopic endometrium were detected among the patients with endometriosis.

As the results, among the drugs or natural compounds identified in Chinese medicine and plants, genistein, atractyloside, naringenin, canadine, ursolic acid, and lycorine showed a higher negative correlation with, and greater potential to effectively treat endometriosis. The natural compounds with highly significant correlations from CMap analysis results are listed in Table 5. Based on the results of CMap analysis, we identified potential natural compounds that could be beneficial in the treatment of endometriosis, and provided a basis for further drug discovery.

TABLE 5

Table 5 The results of Connectivity Map (CMap) analysis.

Discovery of Potential Compounds

Based on network pharmacology prediction, we identified some compounds that may have therapeutic effects, including polyphenolic compounds, sesquiterpenes, terpenoids, flavonoids, alkaloids, polysaccharides, and steroid glycosides. These compounds have better biological activity and are common in the literature regarding the treatment of endometriosis and diseases with similar target pathways. The Venn diagram of comomn targets of major Chinese herbs in endometriosis treatment was established by FunRich, an open access standalone functional enrichment and interaction network analysis tool (Figure 8A) (Pathan et al., 2015). The Venn analysis suggests that the targets of these Chinese herbs are similar, which means that these herbs could play synergistic role. Furthermore, compounds such as ursolic acid, rosmarinic acid, ferulic acid, tanshinone IIA, and oleic acid could regulate pain. Quercetin, salviolone, acetic acid, formononetin, luteolin, hederagenin, tanshinone II A, palmitic acid, cryptotanshinone, rutin, and curcumol could regulate the inflammatory or immunomodulatory response. These compounds could regulate core genes related to endometriosis and with more experimental verification, and could be potentially used in the treatment of endometriosis. The potential effective compounds are listed in Table 6, and the compounds and targets network is presented in Figure 8B. The potential compounds and targets are listed in Supplementary Table 10.

FIGURE 8

Figure 8 Discovery of potential compounds. (A) Venn diagram of major Chinese herbs associated targets. The venn diagram reflects the similarity of these Chinese medicines in regulating targets. (B) Network analysis to decipher the synergistic mechanisms of the compound-target network. We selected natural compounds that may have therapeutic effects for endometriosis and combined that information with molecular docking results. The compounds were derived from data mining of common Chinese herbs. Relationship between formula targets and drug targets. The network was constructed with compounds sharing the same targets between compounds. The data represents the ratio of formula targets and drug targets in biological processes. Blue nodes represent targets and orange nodes represent compounds. The deeper the color of a node, the more frequently it appears. Deeper colors may be the main targets and compounds. The bigger font size, the greater its degree based on undirected network analyisis. The major compounds were showed in Table 6.

TABLE 6

Table 6 The potential effective compounds and targets in endoemtriosis treatment.

Discussion

In this study, data mining and network pharmacology were used to study the rules and potential mechanisms of Chinese medicine in endometriosis. Our research shows that the compatibility of TCM is favorable, and the herbs currently used are mainly blood stasis drugs; the mechanism of action entails multiple pathways. Chinese medicine could provide a new source of screening for potential therapeutic drugs. Such Chinese herbs could also be used to treat other diseases, especially those associated with pain.

Guidelines of Chinese herbs in Endometriosis

Endometriosis is the most common cause of secondary dysmenorrhea. Blood stasis syndrome (BSS) is an important associated pathology in Chinese medicine. Based on the data mining results, most of the common prescriptions for the treatment of endometriosis originate from classic prescriptions in ancient Chinese treatment of dysmenorrhea. These discovered herbs combinations and prescriptions could provide guidance for clinical treatment. The use of data mining methods to screen the core prescriptions of TCM also provides a new method for discovery of TCM.

Machine learning and deep learning methods provide a new method to extract informative medicinal and pharmaceutical research (Liang et al., 2019). Data mining could provide a way to explore the empirical knowledge to promote the development of Chinese medicine from individualized empirical knowledge to large-scale evidence-based medicine (Zhou et al., 2010). The “Framework for Automated Knowledge Graph Construction Towards Traditional Chinese Medicine” in our work could also aid the discovery of pharmaceutic rule of Chinese medicines. These tools are crucial to establish the effectiveness of Chinese medicine/herbs in endometriosis clinical treatment.

We selected some commonly used herbs and HPs to provide a source for network pharmacology. To provide better clinical evidence, the quality of modern clinical research in Chinese medicine is in urgent need of improvement. The mechanisms of animal-derived Chinese medicine, as well as some herbal Chinese medicines, need further analysis and exploration.

The Major Mechanism of TCM

In this study, we discovered potential compounds for the treatment of endometriosis. Studies on phytochemicals have been mostly focused on polyphenolic compounds, sesquiterpenes, terpenoids, flavonoids, alkaloids, polysaccharides, and steroid glycosides. The molecular biological approach has further helped to elucidate the mechanisms of HPs and their pharmacological actions. TCM has advantages in the regulation of endometriosis pain, can directly regulate neurotransmitters, and indirectly regulate pain by modulating inflammation and immunity.

Chinese herbs may interfere with the development of endometriosis, during the processes of invasion, adhesion, angiogenesis, immune inflammation, and oxidative stress, among other processes, such as inhibition of the epithelial-mesenchymal transition. The results of the current in vivo and in vitro experiments are consistent with the analysis of network pharmacology. This suggests that TCM could reduce ectopic lesions and inhibit the development of endometriosis by intervening inflammation, immunity, and inducing apoptosis.

The potential effective compounds based on network pharmacology have also been validated in experiments. The effects of compounds in endometriosis treatment typically induce apoptosis to reduce lesions through various pathways. For example, nerolidol could reduce the average volume of lesions in rats with endometriosis (Melekoglu et al., 2018). Beta-caryophyllene could suppress the growth of endometrial implants in a rat model of endometriosis, without affecting fertility by inducing apoptosis in the luminal epithelium of cysts and in the endothelial cells of blood vessels (Abbas et al., 2013). Tanshinone IIA could also reduce the expression of 14-3-3ζ in ectopic endometrial stromal cells but has no effect on apoptosis (Wan et al., 2015).

To analyze the role of various modules, we discuss only the main modules of Chinese medicine treatment in endometriosis.

The Pain Relief Function

Increased levels of proinflammatory cytokines have been observed from endometriotic lesions, based on examination of the peritoneal fluid of women with endometriosis. These cytokines are also responsible for excessive sensory innervation and the development of chronic pelvic pain (Kobayashi et al., 2014; Zhang J. et al., 2015). The results showed that Chinese herbs could reduce pain associated with the regulation of neurotransmitters (5-hydroxytryptamine, dopamine, and nervous system associated pathways (Wang et al., 2018).

Compounds such as fumarine isocorypalmine fagarine I, in C. anthriscoides ‘Chuanxiong’ and C. yanhusuo could reduce pain via HTR3A, HTR2C, HTR2A. Furthermore, PTGS, also known as cyclooxygenase, is the key enzyme in prostaglandin biosynthesis. It acts both as a dioxygenase and a peroxidase, and may be a core target of TCM for pain. Prostaglandin E2 (PGE2) production by PTGS2/cyclooxygenase 2 (COX-2) is known to be a critical inflammatory factor in endometriosis-associated pain (Sacco et al., 2012).

Various volatile oils, especially terpenes and derivative compounds in Chinese herbs could relieve endometriosis associated pain by regulating prostaglandins and opioidreceptors, in a similar manner to steroidal antiinflammatory drugs (Sacco et al., 2012) or prescribed opioid drugs, cannabis, and cannabinoids (Guimaraes et al., 2014). Chinese medicine compounds could also alleviate pain by regulating inflammatory, immunity, and other related neuronal receptor pathways. For example, C. myrrha could exert analgesic effects, depending on the presence of a biologically active sesquiterpene with a furandadiene skeleton (Shen et al., 2012). Corydalmine, L-tetrahydropalmatine, and demethylcorydalmatine in Curcuma aromatica Salisb. (Curcumae Radix,Yujin) and Curcuma phaeocaulis Valeton(Curcumae Rhizoma,Ezhu) alleviate mechanical hyperalgesia in models of chronic inflammatory and neuropathic pain in mice (Zhou H. et al., 2016).

In addition, BDNF may be a key factor in chronic pelvic pain of endometriosis, and is an effector of estrogen in endometrial cells, thus affecting cell growth (Dong et al., 2017; Ding et al., 2018). Furthermore, BDNF is the potential target of oleic acid, based on the results of network pharmacology.

Chinese medicine could also relieve pain by regulating the inflammatory immune processes and neurotransmitters. Ferulic acid could significantly attenuate vincristine-induced behavioral alterations, as well as electrophysiological and histopathological changes (Vashistha et al., 2017). C. yanhusuo effectively attenuates neuropathic pain, and could exhibit prominent dopamine receptor antagonistic properties(Wang et al., 2016). The herbs that can eliminate blood stasis could also regulate pain through multiple cross-talk pathways. Based on the results of data mining, we found that Chinese medicines often include several blood stasis-associated herbs that can also relieve pain. Further screening of pain-related HPs and effective combinations of compounds could identify new drugs for the treatment of endometriosis and other diseases associated with pain.

Regulatory Hormone Function

Estrogen-mediated associated signaling pathways have important implications for the pathogenesis of endometriosis (Tang et al., 2019). Estrogenreceptor (ER) β could inhibit TNF-α-induced apoptosis and interact with components of the cytoplasmic inflammasome, and finally enhances epithelial-mesenchymal transition (Han et al., 2015). Moreover, we found that estrogen pathway and multiple pathways have cross-talk. Estrogen pathway is related to apoptosis and neural and inflammatory regulatory pathways. Estrogen also play an important role in the regulation of the central serotoninergic system in pain (Paredes et al., 2019). Aromatase is considered a key therapeutic target for regulating the biosynthesis of local estrogen in endometriosis.

At present, the main targets of drug treatment include aromatase, steroid hormone biosynthesis, the estrogen signaling pathway, and drug metabolism-cytochrome P450 pathway. A. sinensis could act in conjunction with estrogen signaling pathways. Genistein could regulate estrogen receptor-α and estrogen receptor-β, and reduce the TNF-α, IL-6, VEGF, and HIF-1α levels in ectopic lesions in the endometriosis murine model (Sutrisno et al., 2018). Chinese medicine may play a synergistic role in the regulation of estrogen and serotonin. The regulatory effect of Chinese herbs/compounds on the estrogen receptor still needs further experimental verification, as its agonistic and inhibitory effects cannot be determined through network pharmacology.

The Antiinvasive, Antiadhesion, and Angiogenesis Functions

Invasion, adhesion, and angiogenesis are the core mechanisms of endometriosis (Gagne, 2003), which was also the core module of our Chinese medicine intervention. Various molecular compounds (genistein, oleanolic acid, luteolin, kaempferol, and ursolic acid) have been involved in antiinvasive activity [via decreased expression of intercellular adhesion molecule-1 (sICAM-1), vascular cell adhesion molecule-1 (VCAM-1), and matrix metalloproteinases (MMPs)] associated with endometriosis (Kuessel et al., 2017). For example, formononetin, an isoflavone from A. membranaceus, showed potential effects on the induction of apoptosis, and suppression of migration and invasion (Zhang J. et al., 2018).

The VEGF pathway is a core pathway in Chinese medicine to regulate angiogenesis. Formononetin also could also promote the proliferation and migration of human umbilical vein endothelial cells (HUVECs) by upregulating VEGF and activating extracellular signal-regulated kinase (ERK) (Liang et al., 2018). Tanshinone IIA could induce angiogenesis by inducing the VEGF/VEGFR2 pathway and CD146 (melanoma adhesion molecule) and regulating angiogenic function in HUVECs (Zhang et al., 2016). Curcumin could reduce ectopic endometrial microvascular density (MVD) and VEGF protein expression in an endometriosis rat model (Zhang et al., 2011). Traditional Chinese herbal medicine for the promotion of blood circulation and removal of blood stasis may become a crucial target for antiangiogenesis in endometriosis.

Antiinflammatory Effect and Regulation of Immune Function

The results of network pharmacology and related molecular biology experiments suggest that the toll-like reporter signal pathway and NF-κB signal pathway could be core pathways of TCM in inflammatory immune regulation. Differential expression of antiinflammatory cytokines (IL-10, TGF-β, TLR4 and NF-κB) occurs in women with endometriosis, and this can promote survival, growth, invasion, differentiation, angiogenesis, and immune escape of the endometriotic lesions (Zhou et al., 2019). Rutin, quercetin, quercitrine, catechin, kaempferol, epicatechin, caffeine, and theobromine are also active compounds that are present in several other types of Chinese herbs. These compounds could form complexes with NF-κB and IκB in the cytoplasm (Yuniwati et al., 2018). Curcumin, epigallocatechin gallate, luteolin, quercetin, resveratrol, caffeic acid phenethyl ester, xanthohumol, genistein, and berberine could inhibit TLR4 activation (Chen C. et al., 2018).

Regulating the Kinase Signalling Pathways Function

Imbalance of protein kinases can lead to uncontrolled cell proliferation by stimulating the tumor formation process, leading to kinase-dependent tumor growth. Chinese medicine could play an active role in regulating the kinase signaling pathways. These phenomena are also related to the occurrence and development of endometriosis. The kinase signaling pathway represents a viable target for the treatment of endometriosis (McKinnon et al., 2016). Major signaling pathways, such as the PI3K-Akt signaling pathway, generally promote survival through the inhibition of proapoptotic factors and activation of antiapoptotic factors.

The PI3K-Akt-mTOR pathway is one of the major signaling pathways in endometriosis. Furthermore, mTOR is a key kinase downstream of PI3K-Akt, which regulates tumor cell proliferation, growth, survival and angiogenesis. The expression of the PI3K-Akt pathway (PI3K, PTEN, Akt, and p-Akt) is related to the severity of endometriosis (Madanes et al., 2019). Natural pentacyclic triterpenoids and their semisynthetic derivatives are involved in cell cycle arrest, apoptosis, and autophagy, triggered by the effects of triterpenoids on TGF-β and HER cell surface receptors. They also induce the PI3KAkt-mTOR and IKK/NF-kB signaling pathways, and the STAT3 and MAPK pathways (Markov et al., 2017). Luteolin, could inhibit the development of endometriosis by regulating the expression of PI3K/Akt signaling, the MAPK signaling pathway, and CCNE1 protein in human VK2/E6E7 and End1/E6E7 cell lines and in the endometriosis mouse model (Yi et al., 2018). Naringenin could induce apoptosis in human endometriosis cells by regulating the MAPK and Akt pathways (Park et al., 2017). Research on the mechanism of TCM treatment in endometriosis based on the kinase pathway can also be a main direction for future research.

The Anticancer Function

Based on the results of gene enrichment analysis of endometriosis, we can also determine that the related genes are associated with many cancer pathways. Endometriosis shows tumor-like behavior. Neoplasia is also a rare but significant complication of endometriosis. It usually occurs in the ovaries, and can often arise in younger patients (Matias-Guiu and Stewart, 2018). Curcumol could inhibit MMP-9 via the c-Jun N-terminal kinase (JNK) 1/2 and NF-κB signaling pathways in breast cancer cells (Ning et al., 2016). Ellagic acid is a naturally occurring polyphenolic compound with strong antioxidant and anticancer properties, which could also inhibit tumor cell migration, extracellular matrix invasion, and angiogenesis (Ceci et al., 2018).

β-elemene, a compound extracted from C. phaeocaulis, C. cassia, and C. myrrha, has proven broad-spectrum antitumor activity and is an effective treatment for several types of tumors (Bi et al., 2018). β-caryophyllene and β-caryophyllene oxide-based natural compounds possess anticancer and analgesic properties (Fidyt et al., 2016). Betulinic acid and betulin are lupane-type pentacyclic triterpenoids with multiple bioactivities, particularly antitumor effects (Zhang D. et al., 2015).

Phytosterols, such as stigmasterol, β-sitosterol, and spinasterol, which are closely related to endometriosis-related targets, could exert anticancer effects by reducing cell cycle progression, inducing apoptosis, and inhibiting tumor metastasis (Shahzad et al., 2017). Hydroxysafflor yellow A could inhibit hepatocellular carcinoma through the inhibition of MMP-2, MMP-9, and p38MAPK signaling pathways in HepG2 cells (Zhang J. et al., 2019). These tumor-related targets are closely related to endometriosis and are the core targets of current research in endometriosis treatment.

Based on the results of gene enrichment, we could more comprehensively consider the regulatory effects of Chinese medicine on neurotransmitter-associated pathways, focal adhesion, PI3K-Akt-mTOR signaling, toll-like receptor signaling, VEGF signaling, and MAPK signaling pathway, among others, which may be used to treat the uterine endometriosis-associated core pathway. The mechanism of Chinese medicine in endometriosis is presented in Figure 9.

FIGURE 9

Figure 9 Potential mechanism of action of Chinese herbs/compounds in endometriosis. Based on the results of our network pharmacological analysis, polyphenolic compounds, sesquiterpenes, terpenoids, flavonoids, alkaloids, polysaccharides, and steroid glycosides are commonly used for the treatment of endometriosis. Traditional Chinese medicine for treating endometriosis includes regulating estrogen receptors, regulating nerve-related receptors, regulating inflammatory and immunity, inhibiting invasion, adhesion, and angiogenesis, affecting the kinase signalling pathways. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that the treatment of endometriosis with Chinese medicine is regulated by multiple cross-talk pathways. Such treatments could act directly on neurotransmitter-related pathways and the estrogen pathway, which is similar to the action of conventional treatment. Chinese medicine also could regulate adhesion, invasion, and angiogenesis by regulating PI3K-Akt signaling, toll-like receptor signaling, NF-κB signaling, and the MAPK signaling pathway, among others.

The Combining Rule and Synergy Effects of HPs and Prescriptions

The formation of Chinese formulas composed of multiple herbs is not random. Data mining provides the basis for our discovery of drug combinations and core prescriptions. The effects of a combination of TCM compounds differ from those of a single natural compound and could play a synergistic role (Yang et al., 2014). HPs play a key role in the investigation of general herbal compatibility and their indispensable effects (Wang et al., 2012). Different compounds/herbs may regulate either the same or different target in various pathways, as well as regulate the enzymes and improve oral drug bioavailability. Systematic identification of the druggable interactions could provide a way to determine the effective mechanisms and new potential drugs for endometriosis treatment (Jiang et al., 2017).

C. phaeocaulis and S. stoloniferum(Ezhu and Sanleng), as a basic HP was first recorded as Sanleng Wan in “Jingyan Liangfang” during the Qing dynasty. The pairs possess extensive activities, including anticancer activity, antiinflammatory activity, and antithrombotic activity (Xu et al., 2015). This HP has shown therapeutic effects on hysteromyoma via the regulation of multiple metabolic pathways. This was demonstrated by UPLC-Q-TOF-MS analysis of a hysteromyoma rat model comprising 16 potential biomarkers from serum and 18 potential biomarkers from urine (Li et al., 2019).

The five herbs in GFW prescriptions appear to play different roles in complementary ways. In addition, GFW could regulate toll-like receptor signaling, TNF signaling, NF-kappa B signaling, neuroactive ligand-receptor interaction, and affect the occurrence and development of endometriosis based on network pharmacology. The pharmacological network functions of C. ramulus target proteins were focused on neurological, inflammatory disease, cancer, cellular growth and proliferation, cell signaling, and molecular transport. Cinnamon extract, active compounds in C. ramulus, could exert antiinflammatory effects by influencing the TLR2 and TLR4 signaling pathways (Schink et al., 2018). Major compounds of P. lactiflora Pall. have been classified as flavonoids, hydrolysable tannins (polyphenols), and monoterpene glycosides. Paeoniflorin could inhibit the plantar incision-induced microglia TLR4/matrix metalloproteinase (MMP-9)/MMP2/IL-1β signaling pathway and suppress postoperative pain in mice (Fan et al., 2018). Paeoniflorigenone could induce apoptosis and has antiproliferative effects in the cancer cell lines (Huang et al., 2017).Paeonol has antipyretic, analgesic, antitumor, antioxidation, antiinflammatory, and immunoregulatory properties (Zhang L. et al., 2019). Paeonol could regulate TLR4/MyD88/NF-κB signaling pathway in lipopolysaccharides-induced acute lung injury C57BL/6 mice (Wang et al., 2019).

These results suggest that the herbal combinations with GFW could have synergistic effects, and may induce apoptosis, antiinflammatory and regulate immunity through the Toll-like receptor signaling, NF-κB signaling pathway.

Other commonly used Chinese medicine prescriptions, such as the ShaoFu Zuyu decoction could reduce the size of ectopic lesions in rats with endometriosis, inhibit cell proliferation, promote apoptosis, and reduce microvessel density and HIF-1α expression (Huang et al., 2016). Jiawei Foshou San, composed of ligustrazine, ferulic acid, and tetrahydropalmatine, could inhibit the progression of endometriosis by regulating epithelial-mesenchymal transformation, based on network pharmacology and experimental verification (Tang et al., 2014; Chen Y. et al., 2018). Network pharmacology could facilitate the analysis of the roles of the synergism effect of compouds or Chinese medicine (Yuan et al., 2017).

Most Chinese herbs/natural compounds may work in a synergistic manner. More research should be focussing on combinatorial effects from HPs/compounds and not just single compounds in drug discovery for complex diseases (Thomford et al., 2018). TCM provides valuable resources for research on complex diseases such as endometriosis. We should adhere to the concept of HPs and prescriptions of TCM. While we still need to use data mining, computer-aided drug discovery and other methods to find more effective herb/medical plant/compound combinations.

This study presented the systemic pharmacological mechanisms of antiendometriosis herbal medicines. A well-designed herbal formula composed of several herbs may synergistically enhance the treatment of specific symptoms. These symptoms may be manifestations of the functional genome. Therefore, using data mining and network pharmacological methods to construct a suitable herbal formula and compound combination could have positive significance.

However, there are still several problems need to be solved. The methodological quality of clinical trials of Chinese herbs in endometriosis needs to be improved. More rigorous research is required to accurately assess the potential roles in endometriosis treatment. Some Chinese medicines are derived from plants of different family and genus. In addition, there are several kinds of names,such as plant names, latin names, and English names, etc., which makes it difficult to confirm medical plants. Accurate scientific nomenclature for plants is essential in the study of ethnopharmacology (Rivera et al., 2014). TCM databases, such as TCMSP, TCMID, TCM-MeSH, are still not comprehensively included in some medical plants/drug targets. Thus, our study used multiple databases to supplement the available information. Network pharmacology analysis needs more evidence to evaluate the potential pharmacological effects of compounds/herbs. We need to further establish comprehensive databases of Chinese medicine and natural compounds to improve compound-related information (Zhang R. et al., 2019).

We should also focus on the pharmacology/bioactivity of bioactive preparations (Heinrich et al., 2020). However, it is confusing that several compounds with low OB, DL in the database have reported many pharmacological values. Such as Salvianolic acid A (OB=2.96%) has several pharmacological actions in antithrombosis, antifibrosis (Xu et al., 2018). In the field of target analyses of network pharmacology, it is difficult to distinguish whether a compound has an effect of inhibition or agitation. The amounts of each bioactive compound in herbs and TCM formulas are very low and Complex Interactions are not fully understood (Eng et al., 2019).

Thus, regarding the compounds and gene targets in the network pharmacology database of Chinese medicine, and the large amount of available information, we still need to conduct more comprehensive evaluations of their therapeutic action based on the literature and in-depth in vivo and in vitro studies, as well as toxic or side effects. We also need to identify ways to improve the current methods. Chinese medicine research requires modernization, such as computational methods (computer-aided drug design) (Chen X. et al., 2018), metabolomics (Wang et al., 2017), mass spectrometry (Zhang A. et al., 2019), and new multiinformational-based profiling approaches integrating taxonomic or bioactivity data (Wolfender et al., 2019), which would present additional techniques by which new compounds can be discovered.

Conclusions

Endometriosis is a common and difficult disease in gynecology. The mechanism of TCM in the treatment of endometriosis is realized through functional modules, such as inhibiting inflammation, enhancing the immune response, regulating angiogenesis-related pathways, inhibiting the epithelial-mesenchymal transition, and inducing apoptosis. The use of data mining combined with bioinformatics techniques could help us understand the associated targets and pathways networks, candidate genes in endometriosis, and the mechanisms of action of Chinese medicine. Combined techniques may also offer an efficient method of drug discovery and development from herbal medicines.

However, there are still several potential limitations in the study and need to improve. First,the mechanism of endometriosis needs to be further clarified based on molecular biology and multiomics technology. Moreover,the antagonistic or agonistic effects of compounds–targets pairs should be further clarified. It is essential to improve the methods in data mining and network pharmacology. Further in vitro or in vivo experiments should be performed to validate the predicted herbs/compounds and targets.

Data Availability Statement

All datasets generated for this study are included in the article/Supplementary Material.

Author Contributions

All authors were responsible for the study concept and design: WZ and JWu have contributed equally for this work. WZ drafted the paper. JWu help drafted the paper. WZ, JWu, and JWa participated in the literature search. WZ, JWu, and HW participated in data mining. WZ, JG, and TW participated in network pharmacology analysis. LC and XL supervised the study. All authors approved the final paper.

Funding

This study is supported by the National Natural Science Foundation of China (81574008), the National Natural Science Foundation of China (61871141), Special Research Project for Traditional Chinese Medicine of Guangdong Hospital of TCM (NO. YN2016ML05), and Situ Yi Renowned Medical Heritage Studio.

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.

Acknowledgments

The authors thank the whole team for their assistance.

Supplementary Material

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

Supplementary Table 1 | The validated information of major herbs, including location, used part, famliy, genus and medical plant reference.

Supplementary Table 2 | The frequent item sets of major Chinese herbs.

Supplementary Table 3 | KEGG pathways of endometriosis associated genes.

Supplementary Table 4 | The major compounds and targets from of major Chinese herbs for endometriosis treatment.

Supplementary Table 5 | The Chinese herbs associated targets names.

Supplementary Table 6 | The major KEGG pathways of major Chinese herbs, herbs pairs and prescription.

Supplementary Table 7 | The targets of endometriosis-associated genes that coincided with those of Chinese herbs.

Supplementary Table 8 | The GO enrichment of endometriosis-associated genes that coincided with Chinese herbs.

Supplementary Table 9 | Drug targets of endometriosis conventional treatment.

Supplementary Table 10 | The potential compounds and targets.

Supplementary Figure 1 | The PPI network of the major endometriosis genes.

Supplementary Figure 2 | The cross-talk pathways network of Chinese herbs in endometriosis treatment.

Abbreviations

ADMET, absorption, distribution, metabolism, excretion and toxicity; BP, biological process; BSS, blood stasis syndrome; CC, cellular component; DEGs, differentially expressed genes; DL, drug-likeness; EMT, epithelial-to-mesenchymal transition; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; MCODE, molecular complex detection; MF, molecular function; OB, oral bioavailability; OMIM, Online Mendelian Inheritance in Man; PPI, protein-protein interactions; RCSB PDB, Research Collaboratory for Structural Bioinformatics Protein Data Bank; TCMSP, Traditional Chinese Medicine Systems Pharmacology; TTD, Therapeutic Target Database; UHPLC-QTOFMS, ultrahigh-performance liquid chromatography tandem quadruple time of flight mass spectrometry.

References

Abbas, M. A., Taha, M. O., Zihlif, M. A., Disi, A. M. (2013). beta-Caryophyllene causes regression of endometrial implants in a rat model of endometriosis without affecting fertility. Eur. J. Pharmacol. 702, 12–19. doi: 10.1016/j.ejphar.2013.01.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Asgarpanah, J., Kazemivash, N. (2013). Phytochemistry, pharmacology and medicinal properties of Carthamus tinctorius L. Chin. J. Integr. Med. 19, 153–159. doi: 10.1007/s11655-013-1354-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Bedaiwy, M. A., Alfaraj, S., Yong, P., Casper, R. (2017). New developments in the medical treatment of endometriosis. Fertil. Steril. 107, 555–565. doi: 10.1016/j.fertnstert.2016.12.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Bi, Y. H., Zhang, L. H., Chen, S. J., Ling, Q. Z. (2018). Antitumor Mechanisms of Curcumae Rhizoma based on network pharmacology. Evid. Based Complement Alternat. Med. 2018, 4509892. doi: 10.1155/2018/4509892

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, J., Farquhar, C. (2014). Endometriosis: an overview of cochrane reviews. Cochrane Database Syst. Rev. 3, CD009590. doi: 10.1002/14651858.CD009590.pub2

CrossRef Full Text | Google Scholar

Ceci, C., Lacal, P. M., Tentori, L., De Martino, M. G., Miano, R., Graziani, G. (2018). Experimental evidence of the antitumor, antimetastatic and antiangiogenic activity of ellagic acid. Nutrients 10, E1756. doi: 10.3390/nu10111756

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, F. P., Kung, Y. Y., Chen, Y. C., Jong, M. S., Chen, T. J., Chen, F. J., et al. (2008). Frequency and pattern of Chinese herbal medicine prescriptions for chronic hepatitis in Taiwan. J. Ethnopharmacol. 117, 84–91. doi: 10.1016/j.jep.2008.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C. Y., Kao, C. L., Liu, C. M. (2018). The cancer prevention, anti-inflammatory and anti-oxidation of bioactive phytochemicals targeting the TLR4 signaling pathway. Int. J. Mol. Sci. 19, E2729. doi: 10.3390/ijms19092729

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Zheng, C., Wang, C., Guo, Z., Gao, S., Ning, Z., et al. (2018). Systems-mapping of herbal effects on complex diseases using the network-perturbation signatures. Front. Pharmacol. 9, 1174. doi: 10.3389/fphar.2018.01174

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Wei, J., Zhang, Y., Sun, W., Li, Z., Wang, Q., et al. (2018). Anti-endometriosis mechanism of Jiawei Foshou San based on network pharmacology. Front. Pharmacol. 9, 811. doi: 10.3389/fphar.2018.00811

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Z., Zhang, C., Gao, F., Fu, Q., Fu, C., He, Y., et al. (2018). A systematic review on the rhizome of Ligusticum chuanxiong Hort. (Chuanxiong). Food Chem. Toxicol. 119, 309–325. doi: 10.1016/j.fct.2018.02.050

PubMed Abstract | CrossRef Full Text | Google Scholar

Chin, C. H., Chen, S. H., Wu, H. H., Ho, C. W., Ko, M. T., Lin, C. Y. (2014). cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst. Biol. 8 Suppl 4, S11. doi: 10.1186/1752-0509-8-S4-S11

PubMed Abstract | CrossRef Full Text | Google Scholar

Cho, K., Van Merrienboer, B., Gülcehre, C., Bahdanau, D., Bougares, F., Schwenk, H., et al. (2014). Learning phrase representations using RNN encoder-decoder for statistical machine translation. Comput. Sci ArXiv, abs/1406.1078. doi: 10.3115/v1/D14-1179

CrossRef Full Text | Google Scholar

de Cassia, D. S. E. S., Lima, T. C., Da, N. F., de Brito, A., de Sousa, D. P. (2017). Analgesic-like activity of essential oil constituents: an update. Int. J. Mol. Sci. 18. E2392 doi: 10.3390/ijms18122392

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, S., Zhu, T., Tian, Y., Xu, P., Chen, Z., Huang, X., et al. (2018). Role of brain-derived neurotrophic factor in endometriosis pain. Reprod. Sci. 25, 1045–1057. doi: 10.1177/1933719117732161

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, F., Zhang, Q., Kong, W., Chen, J., Ma, J., Wang, L., et al. (2017). Regulation of endometrial cell proliferation by estrogen-induced BDNF signaling pathway. Gynecol. Endocrinol. 33, 485–489. doi: 10.1080/09513590.2017.1295439

PubMed Abstract | CrossRef Full Text | Google Scholar

Dunselman, G. A., Vermeulen, N., Becker, C., Calhaz-Jorge, C., D'Hooghe, T., De Bie, B., et al. (2014). ESHRE guideline: management of women with endometriosis. Hum. Reprod. 29, 400–412. doi: 10.1093/humrep/det457

PubMed Abstract | CrossRef Full Text | Google Scholar

Eng, Y. S., Lee, C. H., Lee, W. C., Huang, C. C., Chang, J. S. (2019). Unraveling the molecular mechanism of traditional chinese medicine: formulas against acute airway viral infections as examples. Molecules 24, E3505. doi: 10.3390/molecules24193505

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, Y. X., Hu, L., Zhu, S. H., Han, Y., Liu, W. T., Yang, Y. J., et al. (2018). Paeoniflorin attenuates postoperative pain by suppressing Matrix Metalloproteinase-9/2 in mice. Eur. J. Pain 22, 272–281. doi: 10.1002/ejp.1116

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, Y., Wu, Z., Zhou, X., Zhou, Z., Fan, W. (2006). Knowledge discovery in traditional Chinese medicine: state of the art and perspectives. Artif. Intell. Med. 38, 219–236. doi: 10.1016/j.artmed.2006.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Fidyt, K., Fiedorowicz, A., Strzadala, L., Szumny, A. (2016). β-caryophyllene and β-caryophyllene oxide-natural compounds of anticancer and analgesic properties. Cancer Med. 5, 3007–3017. doi: 10.1002/cam4.816

PubMed Abstract | CrossRef Full Text | Google Scholar

Flower, A., Liu, J. P., Lewith, G., Little, P., Li, Q. (2012). Chinese herbal medicine for endometriosis. Cochrane Database Syst. Rev. 5, CD006568. doi: 10.1002/14651858.CD006568.pub3

CrossRef Full Text | Google Scholar

Fu, J., Wang, Z., Huang, L., Zheng, S., Wang, D., Chen, S., et al. (2014). Review of the botanical characteristics, phytochemistry, and pharmacology of Astragalus membranaceus (Huangqi). Phytother. Res. 28, 1275–1283. doi: 10.1002/ptr.5188

PubMed Abstract | CrossRef Full Text | Google Scholar

Gagne, D. (2003). Levels of vascular endothelial growth factor (VEGF) in serum of patients with endometriosis. Hum. Reprod. 18, 1674–1680. doi: 10.1093/humrep/deg326

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu, J., Gui, Y., Chen, L., Yuan, G., Lu, H. Z., Xu, X. (2013). Use of natural products as chemical library for drug discovery and network pharmacology. PloS One 8, e62839. doi: 10.1371/journal.pone.0062839

PubMed Abstract | CrossRef Full Text | Google Scholar

Gui, Q., Zheng, J. (2019). Simultaneous determination of eight components in Angelica sinensis based on UHPLC-ESI-MS/MS method for quality evaluation. Biomed. Chromatogr. 33, e4326. doi: 10.1002/bmc.4326

PubMed Abstract | CrossRef Full Text | Google Scholar

Guimaraes, A. G., Serafini, M. R., Quintans-Junior, L. J. (2014). Terpenes and derivatives as a new perspective for pain treatment: a patent review. Expert Opin. Ther. Pat. 24, 243–265. doi: 10.1517/13543776.2014.870154

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, S. J., Jung, S. Y., Wu, S. P., Hawkins, S. M., Park, M. J., Kyo, S., et al. (2015). Estrogen receptor beta modulates apoptosis complexes and the inflammasome to drive the pathogenesis of endometriosis. Cell 163, 960–974. doi: 10.1016/j.cell.2015.10.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Heinrich, M., Appendino, G., Efferth, T., Furst, R., Izzo, A. A., Kayser, O., et al. (2020). Best practice in research - Overcoming common challenges in phytopharmacological research. J. Ethnopharmacol. 246, 112230. doi: 10.1016/j.jep.2019.112230

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, X., Su, S., Duan, J. A., Sha, X., Zhu, K. Y., Guo, J., et al. (2016). Effects and mechanisms of Shaofu-Zhuyu decoction and its major bioactive component for cold - stagnation and blood - stasis primary dysmenorrhea rats. J. Ethnopharmacol. 186, 234–243. doi: 10.1016/j.jep.2016.03.067

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, Y., Ohno, O., Suenaga, K., Miyamoto, K. (2017). Apoptosis-inducing activity and antiproliferative effect of Paeoniflorigenone from moutan cortex. Biosci. Biotechnol. Biochem. 81, 1106–1113. doi: 10.1080/09168451.2017.1300517

PubMed Abstract | CrossRef Full Text | Google Scholar

Jamieson, D. G., Moss, A., Kennedy, M., Jones, S., Nenadic, G., Robertson, D. L., et al. (2014). The pain interactome: connecting pain-specific protein interactions. Pain 155, 2243–2252. doi: 10.1016/j.pain.2014.06.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, L., Tang, C., Rao, J., Xue, Q., Wu, H., Wu, D., et al. (2017). Systematic identification of the druggable interactions between human protein kinases and naturally occurring compounds in endometriosis. Comput. Biol. Chem. 71, 136–143. doi: 10.1016/j.compbiolchem.2017.10.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa, M., Goto, S. (2000). KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28, 27–30. doi: 10.1093/nar/28.1.27

PubMed Abstract | CrossRef Full Text | Google Scholar

Kobayashi, H., Yamada, Y., Morioka, S., Niiro, E., Shigemitsu, A., Ito, F. (2014). Mechanism of pain generation for endometriosis-associated pelvic pain. Arch. Gynecol. Obstet. 289, 13–21. doi: 10.1007/s00404-013-3049-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuessel, L., Wenzl, R., Proestling, K., Balendran, S., Pateisky, P., Yotova, S., et al. (2017). Soluble VCAM-1/soluble ICAM-1 ratio is a promising biomarker for diagnosing endometriosis. Hum. Reprod. 32, 770–779. doi: 10.1093/humrep/dex028

PubMed Abstract | CrossRef Full Text | Google Scholar

Lamb, J., Crawford, E. D., Peck, D., Modell, J. W., Blat, I. C., Wrobel, M. J., et al. (2006). The connectivity map: using gene-expression signatures to connect small molecules, genes, and disease. Science 313, 1929–1935. doi: 10.1126/science.1132939

PubMed Abstract | CrossRef Full Text | Google Scholar

Lamb, J. (2007). The connectivity map: a new tool for biomedical research. Nat. Rev. Cancer 7, 54–60. doi: 10.1038/nrc2044

PubMed Abstract | CrossRef Full Text | Google Scholar

Leem, J., Jung, W., Kim, Y., Kim, B., Kim, K. (2018). Exploring the combination and modular characteristics of herbs for alopecia treatment in traditional Chinese medicine: an association rule mining and network analysis study. BMC Complement Altern. Med. 18, 204. doi: 10.1186/s12906-018-2269-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Legendre, G., Delbos, L., Hudon, E., Chabbert-Buffet, N., Geoffron, S., Sauvan, M., et al. (2018). New medical treatments for painful endometriosis: CNGOF-HAS endometriosis guidelines. Gynecol Obstet. Fertil. Senol. 46, 256–263. doi: 10.1016/j.gofs.2018.02.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, W., Hong, B., Li, Q., Li, Z., Bi, K. (2019). An integrated serum and urinary metabonomic research of Rhizoma Curcumae-Rhizoma Sparganii drug pair in hysteromyoma rats based on UPLC-Q-TOF-MS analysis. J. Ethnopharmacol. 231, 374–385. doi: 10.1016/j.jep.2018.11.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, C., Zhou, A., Sui, C., Huang, Z. (2018). The effect of formononetin on the proliferation and migration of human umbilical vein endothelial cells and its mechanism. Biomed. Pharmacother. 111, 86–90. doi: 10.1016/j.biopha.2018.12.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, Z., Liu, J., Ou, A., Zhang, H., Li, Z., Huang, J. X. (2019). Deep generative learning for automated EHR diagnosis of traditional Chinese medicine. Comput. Methods Programs BioMed. 174, 17–23. doi: 10.1016/j.cmpb.2018.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Madanes, D., Bilotas, M. A., Baston, J. I., Singla, J. J., Meresman, G. F., Baranao, R. I., et al. (2019). PI3K/AKT pathway is altered in the endometriosis patient's endometrium and presents differences according to severity stage. Gynecol. Endocrinol., 1–5. doi: 10.1080/09513590.2019.1680627

CrossRef Full Text | Google Scholar

Markov, A. V., Zenkova, M. A., Logashenko, E. B. (2017). Modulation of Tumour-related signaling pathways by natural Pentacyclic Triterpenoids and their semisynthetic derivatives. Curr. Med. Chem. 24, 1277–1320. doi: 10.2174/0929867324666170112115313

PubMed Abstract | CrossRef Full Text | Google Scholar

Matias-Guiu, X., Stewart, C. (2018). Endometriosis-associated ovarian neoplasia. Pathology 50, 190–204. doi: 10.1016/j.pathol.2017.10.006

PubMed Abstract | CrossRef Full Text | Google Scholar

McKinnon, B. D., Kocbek, V., Nirgianakis, K., Bersinger, N. A., Mueller, M. D. (2016). Kinase signalling pathways in endometriosis: potential targets for non-hormonal therapeutics. Hum. Reprod. Update 22, 382–403. doi: 10.1093/humupd/dmv060

PubMed Abstract | CrossRef Full Text | Google Scholar

MEIm, X. D., Cao, Y. F., Che, Y. Y., Li, J., Shang, Z. P., Zhao, W. J., et al. (2019). Danshen: a phytochemical and pharmacological overview. Chin. J. Nat. Med. 17, 59–80. doi: 10.1016/S1875-5364(19)30010-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Melekoglu, R., Ciftci, O., Eraslan, S., Cetin, A., Basak, N. (2018). The beneficial effects of nerolidol and hesperidin on surgically induced endometriosis in a rat model(). Gynecol. Endocrinol. 34, 1–6. doi: 10.1080/09513590.2018.1469611

PubMed Abstract | CrossRef Full Text | Google Scholar

Mihalyi, A., Simsa, P., Mutinda, K. C., Meuleman, C., Mwenda, J. M., D'Hooghe, T. M. (2006). Emerging drugs in endometriosis. Expert Opin. Emerg. Drugs 11, 503–524. doi: 10.1517/14728214.11.3.503

PubMed Abstract | CrossRef Full Text | Google Scholar

Ning, L., Ma, H., Jiang, Z., Chen, L., Li, L., Chen, Q., et al. (2016). Curcumol suppresses breast cancer cell metastasis by inhibiting MMP-9 Via JNK1/2 and Akt-dependent NF-kappaB signaling pathways. Integr. Cancer Ther. 15, 216–225. doi: 10.1177/1534735416642865

PubMed Abstract | CrossRef Full Text | Google Scholar

Paredes, S., Cantillo, S., Candido, K. D., Knezevic, N. N. (2019). An Association of serotonin with pain disorders and its modulation by estrogens. Int. J. Mol. Sci. 20, E5729. doi: 10.3390/ijms20225729

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, S., Lim, W., Bazer, F. W., Song, G. (2017). Naringenin induces mitochondria-mediated apoptosis and endoplasmic reticulum stress by regulating MAPK and AKT signal transduction pathways in endometriosis cells. Mol. Hum. Reprod. 23, 842–854. doi: 10.1093/molehr/gax057

PubMed Abstract | CrossRef Full Text | Google Scholar

Pathan, M., Keerthikumar, S., Ang, C. S., Gangoda, L., Quek, C. Y., Williamson, N. A., et al. (2015). FunRich: an open access standalone functional enrichment and interaction network analysis tool. Proteomics 15, 2597–2601. doi: 10.1002/pmic.201400515

PubMed Abstract | CrossRef Full Text | Google Scholar

Pirzada, A. M., Ali, H. H., Naeem, M., Latif, M., Bukhari, A. H., Tanveer, A. (2015). Cyperus rotundus L.: traditional uses, phytochemistry, and pharmacological activities. J. Ethnopharmacol. 174, 540–560. doi: 10.1016/j.jep.2015.08.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, L., Zheng, X., Liu, J., Li, W., Fu, W., Tang, Q., et al. (2019). Network pharmacology study of traditional Chinese medicines for stroke treatment and effective constituents screening. J. Ethnopharmacol. 242, 112044. doi: 10.1016/j.jep.2019.112044

PubMed Abstract | CrossRef Full Text | Google Scholar

Rivera, D., Allkin, R., Obon, C., Alcaraz, F., Verpoorte, R., Heinrich, M. (2014). What is in a name? The need for accurate scientific nomenclature for plants. J. Ethnopharmacol. 152, 393–402. doi: 10.1016/j.jep.2013.12.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Rogers, P. A., Donoghue, J. F., Walter, L. M., Girling, J. E. (2009). Endometrial angiogenesis, vascular maturation, and lymphangiogenesis. Reprod. Sci. 16, 147–151. doi: 10.1177/1933719108325509

PubMed Abstract | CrossRef Full Text | Google Scholar

Sacco, K., Portelli, M., Pollacco, J., Schembri-Wismayer, P., Calleja-Agius, J. (2012). The role of prostaglandin E2 in endometriosis. Gynecol. Endocrinol. 28, 134–138. doi: 10.3109/09513590.2011.588753

PubMed Abstract | CrossRef Full Text | Google Scholar

Schink, A., Naumoska, K., Kitanovski, Z., Kampf, C. J., Frohlich-Nowoisky, J., Thines, E., et al. (2018). Anti-inflammatory effects of cinnamon extract and identification of active compounds influencing the TLR2 and TLR4 signaling pathways. Food Funct. 9, 5950–5964. doi: 10.1039/c8fo01286e

PubMed Abstract | CrossRef Full Text | Google Scholar

Shahzad, N., Khan, W., Md, S., Ali, A., Saluja, S. S., Sharma, S., et al. (2017). Phytosterols as a natural anticancer agent: Current status and future perspective. Biomed. Pharmacother. 88, 786–794. doi: 10.1016/j.biopha.2017.01.068

PubMed Abstract | CrossRef Full Text | Google Scholar

Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi: 10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, T., Li, G. H., Wang, X. N., Lou, H. X. (2012). The genus Commiphora: a review of its traditional uses, phytochemistry and pharmacology. J. Ethnopharmacol. 142, 319–330. doi: 10.1016/j.jep.2012.05.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Soltani, A., Akbarzadeh-T, M. R. (2014). Confabulation-inspired association rule mining for rare and frequent itemsets. IEEE Trans. Neural Netw. Learn Syst. 25, 2053–2064. doi: 10.1109/TNNLS.2014.2303137

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, S. Y., Muo, C. H., Sung, F. C., Morisky, D. E. (2014). Reduction of surgery rate in endometriosis patients who take Chinese medicine: a population-based retrospective cohort study. Complement. Ther. Med. 22, 632–639. doi: 10.1016/j.ctim.2014.06.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Sutrisno, S., Aprina, H., Simanungkalit, H. M., Andriyani, A., Barlianto, W., Sujuti, H., et al. (2018). Genistein modulates the estrogen receptor and suppresses angiogenesis and inflammation in the murine model of peritoneal endometriosis. J. Tradit. Complement Med. 8, 278–281. doi: 10.1016/j.jtcme.2017.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Szklarczyk, D., Santos, A., von Mering, C., Jensen, L. J., Bork, P., Kuhn, M. (2016). STITCH 5: augmenting protein-chemical interaction networks with tissue and affinity data. Nucleic Acids Res. 44, D380–D384. doi: 10.1093/nar/gkv1277

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, Q., Shang, F., Wang, X., Yang, Y., Chen, G., Chen, Y., et al. (2014). Combination use of ferulic acid, ligustrazine and tetrahydropalmatine inhibits the growth of ectopic endometrial tissue: a multi-target therapy for endometriosis rats. J. Ethnopharmacol. 151, 1218–1225. doi: 10.1016/j.jep.2013.12.047

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, Z. R., Zhang, R., Lian, Z. X., Deng, S. L., Yu, K. (2019). Estrogen-receptor expression and function in female reproductive disease. Cells 8, E1123. doi: 10.3390/cells8101123

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomford, N. E., Senthebane, D. A., Rowe, A., Munro, D., Seele, P., Maroyi, A., et al. (2018). Natural products for drug discovery in the 21st Century: innovations for novel drug discovery. Int. J. Mol. Sci. 19, E1578. doi: 10.3390/ijms19061578

PubMed Abstract | CrossRef Full Text | Google Scholar

Vashistha, B., Sharma, A., Jain, V. (2017). Ameliorative potential of ferulic acid in vincristine-induced painful neuropathy in rats: an evidence of behavioral and biochemical examination. Nutr. Neurosci. 20, 60–70. doi: 10.1179/1476830514Y.0000000165

PubMed Abstract | CrossRef Full Text | Google Scholar

Vidal, M., Cusick, M. E., Barabasi, A. L. (2011). Interactome networks and human disease. Cell 144, 986–998. doi: 10.1016/j.cell.2011.02.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Wan, L., Zou, Y., Wan, L., Wang, L., Huang, M., Wu, J., et al. (2015). Tanshinone IIA inhibits the proliferation, migration and invasion of ectopic endometrial stromal cells of adenomyosis via 14-3-3ζ downregulation. Arch. Gynecol. Obstet. 292, 1301–1309. doi: 10.1007/s00404-015-3766-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Hu, Y., Tan, W., Wu, X., Chen, R., Cao, J., et al. (2012). Compatibility art of traditional Chinese medicine: from the perspective of herb pairs. J. Ethnopharmacol. 143, 412–423. doi: 10.1016/j.jep.2012.07.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Xing, J., Xu, Y., Zhou, N., Peng, J., Xiong, Z., et al. (2015). In silico ADME/T modelling for rational drug design. Q. Rev. Biophys. 48, 488–515. doi: 10.1017/S0033583515000190

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Zhang, Y., Wang, Z., Gong, N., Kweon, T. D., Vo, B., et al. (2016). The Antinociceptive properties of the Corydalis yanhusuo extract. PloS One 11, e0162875. doi: 10.1371/journal.pone.0162875

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M., Chen, L., Liu, D., Chen, H., Tang, D. D., Zhao, Y. Y. (2017). Metabolomics highlights pharmacological bioactivity and biochemical mechanism of traditional Chinese medicine. Chem. Biol. Interact. 273, 133–141. doi: 10.1016/j.cbi.2017.06.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, R., Han, L., Gao, Q., Chen, D., Wang, Y., Zhang, X., et al. (2018). Progress on active analgesic components and mechanisms of commonly used traditional Chinese medicines: a comprehensive review. J. Pharm. Pharm. Sci. 21, 437–480. doi: 10.18433/jpps30212

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, F., Zhu, M., Jiang, N., Zhang, M., Feng, L., Jia, X. (2019). Paeonol ameliorates lipopolysaccharides-induced acute lung injury by regulating TLR4/MyD88/NF-kappaB signaling pathway. Pharmazie 74, 101–106. doi: 10.1691/ph.2019.8860

PubMed Abstract | CrossRef Full Text | Google Scholar

Weng, H., Liu, Z., Yan, S. (2017). A framework for automated knowledge graph construction towards traditional Chinese medicine. Int. Conf. Health Inf. Sci. 10594, 170–181. doi: 10.1007/978-3-319-69182-4_18

CrossRef Full Text | Google Scholar

Wolfender, J. L., Litaudon, M., Touboul, D., Queiroz, E. F. (2019). Innovative omics-based approaches for prioritisation and targeted isolation of natural products - new strategies for drug discovery. Nat. Prod. Rep. 36, 855–868. doi: 10.1039/c9np00004f

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, G. L., Geng, D., Xie, M., Teng, K. Y., Tian, Y. X., Liu, Z. Z., et al. (2015). Chemical composition, antioxidative and anticancer activities of the essential oil: Curcumae Rhizoma-Sparganii Rhizoma, a traditional herb pair. Molecules 20, 15781–15796. doi: 10.3390/molecules200915781

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, J., Wei, K., Zhang, G., Lei, L., Yang, D., Wang, W., et al. (2018). Ethnopharmacology, phytochemistry, and pharmacology of Chinese Salvia species: a review. J. Ethnopharmacol. 225, 18–30. doi: 10.1016/j.jep.2018.06.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Zhang, Z., Li, S., Ye, X., Li, X., He, K. (2014). Synergy effects of herb extracts: pharmacokinetics and pharmacodynamic basis. Fitoterapia 92, 133–147. doi: 10.1016/j.fitote.2013.10.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Yi, F., Li, L., Xu, L. J., Meng, H., Dong, Y. M., Liu, H. B., et al. (2018). In silico approach in reveal traditional medicine plants pharmacological material basis. Chin. Med. 13, 33. doi: 10.1186/s13020-018-0190-0

PubMed Abstract | CrossRef Full Text | Google Scholar

You, X., Xu, Y., Huang, J., Zhi, Y., Wu, F., Ding, F., et al. (2019). A data mining-based analysis of medication rules in treating bone marrow suppression by Kidney-Tonifying method. Evid. Based Complement Alternat. Med. 2019, 1907848. doi: 10.1155/2019/1907848

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, H., Ma, Q., Cui, H., Liu, G., Zhao, X., Li, W., et al. (2017). how can synergism of traditional medicines benefit from network pharmacology? Molecules 22, E1135. doi: 10.3390/molecules22071135

PubMed Abstract | CrossRef Full Text | Google Scholar

Yue, S., Tang, Y., Li, S., Duan, J. A. (2013). Chemical and biological properties of quinochalcone C-glycosides from the florets of Carthamus tinctorius. Molecules 18, 15220–15254. doi: 10.3390/molecules181215220

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuniwati, C., Ramli, N., Purwita, E., Yusnaini, Y., Nurdahliana, N., Miko, A., et al. (2018). Molecular docking for active compounds of Scurrula Atropurpurea as anti-inflammatory candidate in endometriosis. Acta Inform Med. 26, 254–257. doi: 10.5455/aim.2018.26.254-257

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, L., Yang, K., Liu, H., Zhang, G. (2017). A network pharmacology approach to investigate the pharmacological effects of Guizhi Fuling Wan on uterine fibroids. Exp. Ther. Med. 14, 4697–4710. doi: 10.3892/etm.2017.5170

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Cao, H., Hu, Y. Y., Wang, H., Zhang, C. J. (2011). Inhibitory effect of curcumin on angiogenesis in ectopic endometrium of rats with experimental endometriosis. Int. J. Mol. Med. 27, 87–94. doi: 10.3892/ijmm.2010.552

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, D. M., Xu, H. G., Wang, L., Li, Y. J., Sun, P. H., Wu, X. M., et al. (2015). Betulinic acid and its derivatives as potential antitumor agents. Med. Res. Rev. 35, 1127–1155. doi: 10.1002/med.21353

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Echeverry, S., Lim, T. K., Lee, S. H., Shi, X. Q., Huang, H. (2015). Can modulating inflammatory response be a good strategy to treat neuropathic pain? Curr. Pharm. Des. 21, 831–839. doi: 10.2174/1381612820666141027115508

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Q., Liu, X., Guo, S. (2016). Progressive development of endometriosis and its hindrance by anti-platelet treatment in mice with induced endometriosis. Reprod. Biomed. Online 34, 124–136. doi: 10.1016/j.rbmo.2016.11.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, A., Sun, H., Wang, X. (2018). Mass spectrometry-driven drug discovery for development of herbal medicine. Mass Spectrom. Rev. 37, 307–320. doi: 10.1002/mas.21529

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Liu, L., Wang, J., Ren, B., Zhang, L., Li, W. (2018). Formononetin, an isoflavone from Astragalus membranaceus inhibits proliferation and metastasis of ovarian cancer cells. J. Ethnopharmacol. 221, 91–99. doi: 10.1016/j.jep.2018.04.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Li, J., Song, H., Xiong, Y., Liu, D., Bai, X. (2019). Hydroxysafflor yellow A suppresses angiogenesis of hepatocellular carcinoma through inhibition of p38 MAPK phosphorylation. Biomed. Pharmacother. 109, 806–814. doi: 10.1016/j.biopha.2018.09.086

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L., Li, D. C., Liu, L. F. (2019). Paeonol: pharmacological effects and mechanisms of action. Int. Immunopharmacol. 72, 413–421. doi: 10.1016/j.intimp.2019.04.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, R., Zhu, X., Bai, H., Ning, K. (2019). Network pharmacology databases for traditional Chinese medicine: review and assessment. Front. Pharmacol. 10, 123. doi: 10.3389/fphar.2019.00123

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Lin, H., Yang, Z., Wang, J., Sun, Y., Xu, B., et al. (2019). Neural network-based approaches for biomedical relation classification: a review. J. Biomed. Inform. 99, 103294. doi: 10.1016/j.jbi.2019.103294

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhong, L. J., Hua, Y. L., Ji, P., Yao, W. L., Zhang, W. Q., Li, J., et al. (2016). Evaluation of the anti-inflammatory effects of volatile oils from processed products of Angelica sinensis radix by GC-MS-based metabolomics. J. Ethnopharmacol. 191, 195–205. doi: 10.1016/j.jep.2016.06.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, X., Chen, S., Liu, B., Zhang, R., Wang, Y., Li, P., et al. (2010). Development of traditional Chinese medicine clinical data warehouse for medical knowledge discovery and decision support. Artif. Intell. Med. 48, 139–152. doi: 10.1016/j.artmed.2009.07.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, H. H., Wu, D. L., Gao, L. Y., Fang, Y., Ge, W. H. (2016). L-Tetrahydropalmatine alleviates mechanical hyperalgesia in models of chronic inflammatory and neuropathic pain in mice. Neuroreport 27, 476–480. doi: 10.1097/WNR.0000000000000560

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y., Xie, M., Song, Y., Wang, W., Zhao, H., Tian, Y., et al. (2016). Two traditional Chinese medicines Curcumae Radix and Curcumae Rhizoma: an ethnopharmacology, phytochemistry, and pharmacology review. Evid. Based Complement Alternat. Med. 2016, 4973128. doi: 10.1155/2016/4973128

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, W. J., Yang, H. L., Shao, J., Mei, J., Chang, K. K., Zhu, R., et al. (2019). Anti-inflammatory cytokines in endometriosis. Cell. Mol. Life Sci. 76, 2111–2132. doi: 10.1007/s00018-019-03056-x

PubMed Abstract | CrossRef Full Text | Google Scholar