Network pharmacology and bioinformatics were used to construct a prognostic model and immunoassay of core target genes in the combination of quercetin and kaempferol in the treatment of colorectal cancer

Abstract

Purpose: CRC is a malignant tumor seriously threatening human health. Quercetin and kaempferol are representative components of traditional Chinese medicine (TCM). Previous studies have shown that both quercetin and kaempferol have antitumor pharmacological effects, nevertheless, the underlying mechanism of action remains unclear. To explore the synergy and mechanism of quercetin and kaempferol in colorectal cancer.

Methods: In this study, network pharmacology, and bioinformatics are used to obtain the intersection of drug targets and disease genes. Training gene sets were acquired from the TCGA database, acquired prognostic-related genes by univariate Cox, multivariate Cox, and Lasso-Cox regression models, and validated in the GEO dataset. We also made predictions of the immune function of the samples and used molecular docking to map a model for binding two components to prognostic genes.

Results: Through Lasso-Cox regression analysis, we obtained three models of drug target genes. This model predicts the combined role of quercetin and kaempferol in the treatment and prognosis of CRC. Prognostic genes are correlated with immune checkpoints and immune infiltration and play an adjuvant role in the immunotherapy of CRC.

Conclusion: Core genes are regulated by quercetin and kaempferol to improve the patient's immune system and thus assist in the treatment of CRC.

Introduction

CRC is a swart tumor extremely threatening to human health. In recent years, CRC has not only shown a trend of younger age but also has a lower survival rate for young patients [1]. Due to the insidious early clinical symptoms, about 25% of patients have metastases at initial diagnosis [2]. Metastasis also occurs in 50% to 60% of patients after treatment, so a radical surgical resection of CRC is becoming increasingly difficult [3]. Therefore, the treatment of CRC has become a major challenge for the medical community. CRC treatment includes surgery, chemotherapy, immunotherapy, and targeted therapy. Nevertheless, the prognosis of CRC patients remains suboptimal due to the limitations of surgery, radiation dose limitation, toxic side effects of chemotherapy, and drug resistance [4]. At present, targeted therapy, immunotherapy, and other therapies have made significant progress in CRC treatment, but some adverse reactions will also occur during the treatment process, which will seriously reduce the patient's quality of life [5]. Immune checkpoint inhibitors (ICI) for immunotherapy, and ICI activation can cause some toxic side effects, which is a momentous fling down the gauntlet in ICI clinical application [6]. Currently, the therapeutic goal of CRC is to prolong overall survival and improve patient quality of life.

Some studies have shown that the curative effect of traditional Chinese medicine on liver cancer [7], breast cancer [8], and so on seems to be very promising. Therefore, it is particularly important to clarify the anticancer mechanism of traditional Chinese medicine treatment. In clinical studies, quercetin is used in the treatment of diabetes [9-13], hyperlipidemia [14], and nonalcoholic fatty liver disease (NAFLD) [15-16]. Some preclinical studies have found that quercetin has anticancer effects by inhibiting cell proliferation and metastasis, inducing apoptosis and autophagy [17-20], and can also induce apoptosis in HCT116 and HT29 cells [21-23]. However, there are few relevant clinical studies on quercetin and antitumor. In addition, it has been reported that kaempferol can also reduce the risk of skin and liver cancer [24] and inhibits the growth of CRC [25] by accelerating the necrosis and apoptosis of cells. Quercetin and kaempferol not only play an anti-tumor effect alone but also inhibit the proliferation of colorectal cancer cells when used together [26]. The above studies show that kaempferol and quercetin are potential drugs in CRC treatment.

In the development of TCM, researchers have achieved good preliminary results in revealing the comprehensive effect of multi-channel, multi-target, and multi-component of TCM by referring to the research ideas of network pharmacology [27-29]. However, the pharmacological mechanism of Chinese herbal medicine (CHM) has not been fully elucidated, so the establishment of the CHM database is especially important for network pharmacology research [30] which can be used for mutual verification of network analysis and experiments. Therefore, network pharmacology is being used as a very useful tool to understand drug-target interactions. The purpose of this research was to explore the mechanism of action of quercetin and kaempferol in CRC treatment. The targeting network of quercetin and kaempferol for the treatment of CRC was constructed by the method of network pharmacology, the core genes were screened, and the drug target gene network composed of fifty-seven target genes was obtained. Based on our findings, we can all the better understand the reciprocity between these two drugs and supply theoretical backing for future basic study. Then we analyzed the genes differentially expressed between CRC and para-cancerous tissues to so much better predict the survival time of CRC patients treated with quercetin and kaempferol and to provide a reference for clinical drug use. Univariate Cox together with multivariate Cox regression was used to analyze 49 Differentially expressed genes (DEG) in the drug targeting network and obtain three key prognostic genes. These three prognostic genes were analyzed by the lasso-Cox regression analysis algorithm, and the nomogram predictive map was set up to boost clinical guidance. In addition, the function of genes in this model in tumor occurrence and development and immunologic response was studied by gene cluster enrichment analysis (GSEA).

Finally, to study the immune-related effect of quercetin and kaempferol, we predicted a correlation between the expression of the core genes and immune checkpoint genes and further discussed the effects of quercetin and kaempferol on the immunologic response.

Materials and Methods

Data collection

Based on the PubChem database and TCMSP database, quercetin and kaempferol target genes were obtained. And from the Gene Cards database, Online Mendelian Inheritance in Man (OMIM), Pharmacogenetics and Pharmacogenomics Knowledge Base (PharmGkb), Database of gene-disease associations (DisGeNET), and Treatment Target Database (TTD) were used to obtain CRC disease genes. The training gene set was obtained from TCGA and validation gene sets were obtained from the GEO database. Immunohistochemical images of gene expression in para-cancerous tissues and tumors were obtained from the Human Protein Atlas Database (HPA).

Construction and screening of drug-target interaction network

The targets of quercetin and kaempferol were merged as drug targets, and the acquired genes in five databases were merged as pathogenic gene sets to obtain the two crossover gene sets. Protein-protein interaction (PPI) networks were constructed through the STRING database, and outlier genes were removed. The PPI network was imported into Cytoscape 3.9.1, and the network core genes were screened using the CytoNCA [31] plugin. The screening criteria are betweenness centrality (BC), degree centrality (DC), closeness centrality (CC), eigenvector centrality (EC), local average connectivity (LAC), and network centrality (NC). The filter conditions are all greater than the median value.

Functional and functional enrichment analysis of the drug-target network

To explore the significant functions and cellular pathways of the core genes. The Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analysis were performed using the R package “cluster Profiler” [32]. The pathways were significantly enriched when a P value of < 0.05 and an FDR of < 0.25 were considered.

Core gene analysis

The RNA-seq data in the TCGA database were grouped into tumor groups and para-cancerous groups, heat maps and boxplots of core gene expression levels were plotted, and correlation maps of gene expression were plotted. Undifferentiable genes were removed to obtain the differentially expressed genes (DEGs).

Construction of the prediction model of Lasso-cox regression

Univariate and multivariate Cox regression analyses were used to determine the prognostic value of DEGs (P < 0.05). In this research, we used the R package glmnet [33], using the lasso-cox method to identify CRC prognostic genes. We calculated the optimal cutoff for RiskScore using the R package maxstat [34], based on which we differentiated patients into high and low groups and generated Kaplan-Meier survival curves to evaluate the predictive performance of the associated risk genes.

Validation of the Lasso-Cox regression prediction model

The CRC transcriptome and clinical data of two gene chips in GEO were used to verify the model. We use R-packet pROC [35] for ROC analysis to verify the performance of the model prediction. The risk curve, risk and survival scatter plots, and heat maps of the validation model gene expression are plotted. Forest plots were drawn based on univariate Cox regression analysis. We used the R software package “rms” [36] and the Cox method to establish a nomogram to evaluate the prognostic significance of some features in the sample and predict the survival of patients.

Correlation between prognostic model and tumor immunity

The data on immune cell level were downloaded from TIMER2.0 [37] online immune database, and the bar chart and correlation chart of immune cell content were drawn. To explore the pertinence between the gene expression and immune cells in Lasso regression prognostic model, we used the online analysis tool SangerBox [38] to draw the pertinence map between 3 prognostic genes expression and immune cells, and the pertinence map between 3 prognostic genes expression and immune checkpoint. We selected ESTIMATE [39] in R package IOBR [40] to score the immune infiltrating cells of CRC samples in TCGA. And draw the Kaplan-Meier survival curve of the relationship between the content of stromal cells and immune cells and the survival of patients.

Results

Sample collection

According to inclusion criteria [41]: (a) complete gene expression information (b) complete survival information, the exclusion standard was as follows: (a) incomplete gene expression (b) no survival time and survival status in the clinical information. Transcriptome data and 627 clinical data from 625 CRC organizations and 51 para-cancerous tissues were acquired from the TCGA database, including overall survival time, survival status, age, gender, and other data (Table S1). A sum of 728 transcriptome and clinical information (Table S2) were obtained from GEO's two gene chips, GSE103479 [42] and GSE39582 [43].

Construction of the drug-target interaction network

In this research, we used the study notion of network pharmacology, given the possible loss of genes in individual disease databases, 12549 CRC causal genes were obtained from 5 databases (Figure 1 (A)). 166 quercetin targets and 64 kaempferol targets were obtained from the PubChem database and TCMSP database, and 166 drug targets were removed (Figure 1 (B)). The target genes of quercetin and kaempferol have 155 overlapping genes of CRC (Figure 1 (C)), which are the targets of quercetin and kaempferol for the treatment of CRC. A PPI network containing 154 target proteins (Figure 1 (D)) was obtained through the STRING online database, which was imported into Cytoscape 3.9.1 to obtain the interaction network map between the drug and target genes (Figure 1 (E)). Using plug-in CytoNCA to filter overlapping memes, 57 core memes were obtained (Table S3) (Figure 1 (F)): AKT1, TP53, TNF, IL6, VEGFA, JUN, IL1B, CASP3, PTGS2, HIF1A, MYC, EGFR, EGF, MMP9, CXCL8, CCND1, PTEN, PPARG, CCL2, FOS, IL10, MMP2, ICAM1, NFKBIA, ERBB2, HMOX1, VCAM1, CASP8, BCL2L1, RELA, IFNG, SERPINE1, CDKN2A, STAT1, MAPK8, IL2, NOS3, MAPK1, CDKN1A, IKBKB, PPARA, MMP3, AR, CRP, CXCL10, CAV1, NFE2L2, MPO, PGR, PARP1, IGFBP3, HSPB1, PRKCA, RUNX2, COL1A1, AHR, and RAF1.

Figure 1 (A) A Venn diagram of the five disease-gene databases. (B) Drug target interaction network. (C) Venn chart of the overlapping of drug targets and CRC genes. (D) Protein-protein interplay. (E, F) Screening of the key genes in the network.

Figure 2 (A) GO bubble plot of the network key genes. (B) Bar graph of the network core gene KEGG pathway.

GO and KEGG pathway enrichment analysis

To elucidate the functions and potential actions of drugs versus core genes, we performed GO and KEGG enrichment analyses (Table S4). The GO enrichment results showed that quercetin and kaempferol were enriched in the functions of “regulation of cell death”, “response to oxygen-containing compound”, “response to cytokine”, “apoptotic process”, “response to lipid”, “cytokine-mediated signaling pathway”, “circulatory system development”, “negative regulation of cell death”, “response to abiotic stimulus”, “regulation of cell differentiation”, “membrane microdomain”, “transcription regulator complex”, “chromatin”, “caveola”, “chromosome”, “rna polymerase ii transcription regulator complex”, “death-inducing signaling complex”, “plasma membrane raft”, “plasma membrane protein complex”, “enzyme binding”, “signaling receptor binding”, “transcription factor binding”, “dna binding transcription factor binding”, “cytokine receptor binding”, “rna polymerase ii specific dna binding transcription factor binding”, “identical protein binding”, “protein-containing complex binding”, “kinase binding” and so on. (Figure 2 (A)). KEGG enrichment results show that the main signal pathways of 57 core genes are located at “Pathways in cancer”, “AGE-RAGE signaling pathway in diabetic complications”, “TNF signaling pathway”, “IL-17 signaling pathway”, “Kaposi sarcoma-associated herpesvirus infection”, “Fluid shear stress and atherosclerosis”, “Chagas disease (American trypanosomiasis)”, “Bladder cancer”, “Hepatitis B”, “Pancreatic cancer”, “Human cytomegalovirus infection”, “Hepatitis C”, “Prostate cancer”, “C-type lectin receptor signaling pathway”, “Influenza A”, “Human T-cell leukemia virus 1 infection”, “Toll-like receptor signaling pathway”, “HIF-1 signaling pathway”, “Measles”, “PD-L1 expression and PD-1 checkpoint pathway in cancer”, “Proteoglycans in cancer”, “Yersinia infection”, “Endocrine resistance”, “Relaxin signaling pathway”, “Th17 cell differentiation”, “Apoptosis” and “Colorectal cancer” et al.(Figure 2(B)).

Construction and validation of the Lasso regression prediction model

First of all, to explore the expression levels of the above 57 core genes in normal tissues and tumor samples, we drew a heat map (Figure 3 (A)) and a violin map (Figure 3 (B)). Among them, TNF, RUNX2, PTGS2, MMP2, MAPK8, HSPB1, AKT1, and CRP showed no difference in expression between normal and tumor tissues. Genes with no difference were removed, and 49 DEGs were obtained. Next, to explore the expression correlation of these DEGs, we draw a co-expression association map of genes (Figure 3 (C)). There were 5 groups with a co-expression correlation coefficient greater than 0.7: STAT1 and CXCL10 (correlation coefficient: 0.79), CXCL8 and IL1B (correlation coefficient: 0.78), CCL2 and VCAM1 (correlation coefficient: 0.75), CXCL8 and MMP3 (correlation coefficient: 0.72), IFNG and STAT1 (correlation coefficient: 0.71). Then, we executed univariate Cox regression analysis (Figure 3 (D)) and multivariate Cox prognosis (D) (Figure 3 (E)), and we adopted the R package glmnet, integrating survival time, survival status, and gene expression data for prognostic regression analysis using the lasso-cox method (Figure 3 (F) and Figure 3 (G)). When λ = 0.0383989011310987, 3 prognostic model genes: CDKN2A, MMP3, and SERPINE1 were obtained. The regression coefficients were Coef (CDKN2A) = 0.082101134622246, Coef (MMP3) = -0.0355844623404235, Coef (SERPINE1) = 0.0100737687120257 (Table S5).

The risk score was counted according to the risk score formulary: RiskScore=0.082101134622246 * CDKN2A-0.0355844623404235 * MMP 3 +, 0.0100737687120257 * SERPINE1. The topgallant cutoff value was counted using the R package maxstat as 0.0487909863090454. Based on this, dividing patients into high and low groups (Figure 3 (H)), the risk score can differentiate the survival time and survival status of the high-risk and low-risk populations. The discrepancy in prognosis in both the high and low-risk groups was evaluated by log-rank test, and the prognosis showed a marked difference between the two groups (P < 0.001), and patients in the high-risk group may have shorter survival and worse survival status than those in the low-risk group. We drew ROC curves (Figure 3 (I)) to validate model prediction exactitude and dependability. The AUC value is > 0.6, which proves that the Lasso-Cox prognostic regression curve is credible. Based on this, we drew the risk curve, hazard, survival scatter diagram (Figure 3 (J)), and the expression heat map of prognostic genes (Figure 3 (K)) and analyzed the nexus between diverse risk scores and patient follow up time, survival status and gene expression. With the rise of risk score, the survival ratio of patients reduced significantly (3J), which indicated that the MMP 3 gene was a protective factor, downregulated with the rise of risk score, CDKN2A, and ERPINE1 genes are risk factors and were upregulated with the rise of risk score. In conclusion, we found that the death probability and gene expression level of high-risk patients were higher than those in the low-risk group.

To validate the universality of the Lasso-Cox prediction regression model, we verified the model using the GeneChip from the GEO database (Table S6). The transcriptomic information in the chip was first normalized. Combine gene expression and clinical information to draw survival curves into high-risk and low-risk groups and explore the survival of high/low-risk groups in the validation data set (Figure 4 (A)) (P < 0.001). The ROC curve is between (0.73-0.0.77) AUC, indicating that the pattern is relatively accurate (Figure 4 (B)). The risk curve, hazard, and survival scatter diagram (Figure 4 (C)), and the model gene expression heatmap (Figure 4 (D)) were significantly different in both the high-risk and low-risk groups. These results demonstrate the universality and accuracy of the Lasso-Cox predictive regression model.

We next performed a one-factor outcome analysis (Figure 4 (E)) and found that risk values were highly interrelated with specimen survival status in extensive clinical data. Combining the risk score formulary and survival period to draw the survival prognosis nomogram (Figure 4 (F)) can predict the survival of patients at 1-5 years.

Figure 3 (A) Heatmap of network core gene expression in the TCGA database. (B) Boxplot of differential expression of network core genes. (C) Correlation plot of core differential gene expression. (D) Unigenic Cox prognostic regression analysis of core differential genes. (E) Multigenic Cox prognostic regression analysis of core differential genes. (F, G) Lasso prognostic regression model. (H) Survival curves in the TCGA database. (I) ROC curve to verify the accuracy of risk. (J) Risk curve, risk, and survival scatter plot. (K) Heatmap of model gene expression in the high- and low-risk groups.

Figure 4 (A) Survivorship curves in the GEO database. (B) ROC curve to validate the accuracy of risk. (C) Risk curve, risk, and survival scatter plot. (D) Heatmap of model gene expression in the high- and low-risk groups. (E) Univariate prognostic regression test. (F) Used to forecast survival in patients with the nomogram. (G)The predicted binding model of quercetin and CDKN2A (G), SERPINE1 (H), and MMP3 (I) proteins, and the predicted binding model of kaempferol to MMP3 (J) proteins.

Finally, to explore whether quercetin and kaempferol interact and function with the three key genes, we used Molecular Operating Environment (MOE) 2019 software to draw the binding model of CDKN2A[44] (Figure 4 (G)), SERPINE1[45] (Figure 4 (H)) and MMP3[46] (Figure 4 (I)) to quercetin, and the binding profile of MMP3 (Figure 4 (J)) to kaempferol, demonstrated that quercetin and kaempferol binding to target proteins may produce pharmacological effects.

Figure 5 (A, B, C) GSEA enrichment analysis of CDKN2A, MMP3, and SERPINE1, respectively. (D) Histogram of the multifarious immune cell contents in the TCGA database. (E) Graph of correlations between the levels of multifarious immune cells. (F, G, H) A linear regression graph of various immune cell levels and prognostic gene expression. (I, J, K) Linear regression plots of immune infiltration scores and prognostic gene expression relationships in different samples. (L, M, N) The map of the linear regression between expression levels of multifarious immune checkpoints and prognostic gene expression. (O, P) Survival plots of stromal cells versus immune cells in the samples.

The GSEA cellular pathway enrichment analysis of the model genes

To investigate the association of tumor-related and immune-related pathway genes in the Lasso-Cox prognostic regression model, we performed a GSEA [47] pathway-related enrichment analysis for prognostic genes (Table S7) (The pathways were regarded to be significantly enriched when a P value of < 0.05 and an FDR of < 0.25). The GSEA pathway enrichment results of SERPINE1 (Figure 5 (C)) genes were the most, with 103 significantly related pathways, among which 30 signal pathways were related to tumor development and 18 immune-related signal pathways; MMP3 (Figure 5 (B)) has 33 related pathways, and 7 with both tumor-related and immune-related; CDKN2A (Figure 5 (A)) showed 24 enrichment results, 2 related to tumor, and 1 related to immune. This shows that the three prognostic genes screened above are closely related to tumorigenesis and development and immunotherapy.

Relationship of model genes to tumor immunity

To further explore whether the genes in the model affected the immune response. We downloaded the immune prediction data (Table S8) of CRC samples in TCGA from the TIMER2.0 database and drew the histogram of various immune cell content in some samples (Figure 5 (D)) and the immune cell correlation map (Figure 5 (E)). The relevancy plot shows that T cell CD8+ and Macrophage M0, T cell follicular helper, and Macrophage M0 had a minimal correlation of -0.47; The maximum correlation between B cell memory and T cell gamma delta was 0.73. Next, to explore the relevance between the three prognostic genes and immune cell expression, immune infiltration, and immune checkpoints (Table S9), we drew several scatter plots (Figure 5 (F) -5 (K)). As shown in the image, CDKN2A was significantly associated with T cell CD4, T cell CD8, Neutrophil, Macrophage, and Dendritic cell (P < 0.05); SERPINE1 showed significant associations with B cell, T cell CD4, T cell CD8, Neutrophil, Macrophage and Dendritic cell (P < 0.05); MMP3 were significantly associated with T cell CD8, Neutrophil and Dendritic cell (P < 0.05). The expression of the prognostic genes CDKN2A, SERPINE1, and MMP3 was significantly and positively correlated with the immune infiltration score of the samples (P < 0.001). These results indicate that three prognostic-related genes have a significant association with immune cells.

In addition, we calculated the Pearson correlation of CDKN2A, SERPINE1, and MMP3, and five categories of immune pathway markers, and prognostic gene expression was closely related to immune checkpoints (Figure 5 (L) -5 (N)). Thus, the prognostic genes of the Lasso-Cox prognostic regression model were associated with immune cells and significantly positively associated with the immune checkpoint. In addition, we also performed the survival analysis of the tumor microenvironment of the samples (Table S10) and drew the survival curve of the stromal cells and the immune cells in the samples (Figure 5 (O) -5 (P)).

The results implied that there was no relativity between stromal cell content score and the survival status of the patient, but there was a correlation between immune cell content score and patient survival status, indicating that the more the immune cell content in the tumor tissue, the longer the survival time, that is, the deeper the degree of immune infiltration, the longer the patient's survival time may be.

These results suggest that there is a high correlation between the above three prognostic genes and immune cells, immune infiltration, and immune checkpoints, which can be used as targets for drug development in the immune microenvironment.

Immunohistochemical expression levels in the HPA database

Finally, to validate the expression of the above three prognostic-related genes at the tissue level, we selected the same immunofluorescence antibodies in both normal and tumor samples based on the HPA database (Table S11, Figure S1, S2, S3). And observe Immunohistochemistry images of CDKN2A (Figure 6 (A)-(D)), SERPINE1 (Figure 6 (F)-(I)), MMP3 (Figure 6 (K)-(N)), IOD values in five genes were differentially expressed in normal and cancer tissues (Figure 6 (E), 6 (J), 6 (O)). This suggests that these three genes could serve as therapeutic targets for CRC.

Discussion

At present, tumor treatment methods containing surgical intervention, radiotherapy, chemical therapy, targeted therapy, biological therapy, immunization therapy, and traditional Chinese medicine as complementary and alternative medicine (CAM) therapy also play a significant role [48]. Traditional Chinese medicine has the characteristics of diverse sources, high efficacy, and few toxic and side effects, and is widely welcomed in the world [49-53]. Chinese herbal medicine plants contain a variety of active substances, which can effectively prevent and treat common clinical diseases. Due to the complex mechanism and progression of many diseases, traditional single-target drugs cannot cure some intricate diseases, such as cancer, angiocardiography, and neurological diseases. CHM contains a variety of active compounds, indicating that multi-components of herbal medicine have a synergistic action [54-55]. Therefore, for the therapy of complex diseases for instance cancer, the combination of CHM compounds or combined chemotherapy drugs may be more suitable. It has been shown that a combination of curcumin, resveratrol, and quercetin can work together to treat cancer. For example, resveratrol combined with curcumin significantly bated the proliferation and induced apoptosis of hepatocellular carcinoma cells [56], and quercetin combined with resveratrol significantly increased the reversal of drug resistance of leukemic cells [57]. Curcumin combined with quercetin significantly promoted tumor cell apoptosis in vivo [58].

Figure 6 (A, B, C, D) Immunohistochemical figures of CDKN2A in normal and CRC tissues from the HPA database. (F, G, H, I) Immunohistochemical figures of SERPINE1 in normal and CRC tissues from the HPA database. (K, L, M, N,) Immunohistochemical figures of MMP3 in normal and CRC tissues from the HPA database. (E, J, O) Boxplots of IOD differences in CDKN2A, SERPINE1, and MMP 3 in normal and cancer tissues, respectively.

Compared with monotherapy, Chinese herbal medicine is more suitable for multi-drug combinations. For example, HangAmDan-B1 (HAD-B1), consisting of four herbs—Panax notoginseng, Cordyceps militaris, ginseng, and frankincense—significantly inhibits HCC827-GR cell cycle arrest and accelerates apoptosis compared with monotherapy [59]. Cepharanthine is an alkaloid extracted from Stephania cepharantha Hayata that has anti-inflammatory, and antioxidant effects, and can also induce autophagy, apoptosis, and cell cycle arrest in breast cancer cells [60]. The main active compound of green tea, epigallocatechin-3-gallate (EGCG), has shown preventive effects in different cancer models [61-62]. In addition, Ganoderma lucidum extract has also been used in traditional Chinese medicine to prevent or treat a variety of diseases [63-68].

Studies have claimed that green tea extract (GTE) and Ganodorma lucidum extract (GLE) can inhibit breast cancer cell proliferation, migration, and the combination of the two also has a synergistic effect, and the combination of low concentrations of GLE and GTE significantly inhibits the growth and invasion of MDA-MB-231 metastatic breast cancer [69]. Studies have also shown that the appropriate addition of CHM or CHM compounds to the diet of chemotherapy-resistant patients can enhance immunity, reverse chemotherapy resistance, enhance the toxicity of drugs to cells, reduce side effects, and prolong the life of patients [70]. Astragalus polysaccharides significantly reduced the side effects of vinorelbine and cisplatin in patients with advanced NSCLC and significantly improved their quality of life [71]. In the combination therapy of Chinese herbal medicine, attention needs to be paid to the combination, because not all combinations are safe [72], and there are enormous differences between experimental data and clinical application. Therefore, it is necessary to use some methods to optimize the combination of traditional Chinese medicines, such as the research and analysis method of network pharmacology combined with bioinformatics to study, evaluate and analyze the synergistic effects of the combination of traditional Chinese medicines [73].

In this research, network pharmacology and combined bioinformatics were used to analyze the targets of traditional Chinese medicine extracts quercetin and kaempferol, construct and screen the interaction network between drugs and CRC pathogenic genes, and finally screen three prognostic related genes: CDKN2A, SERPINE1, and MMP3 [74-79].

In recent years, the anticancer activity of quercetin has attracted widespread attention. Pharmacological studies have found that quercetin inhibits tumor growth and induces tumor autophagy and apoptosis by regulating various signaling pathways, for instance, PI3K-Akt signaling pathway, NF-κB signaling pathway, mTOR signaling pathway, MAPK signaling pathway [80], TGF-β1/Smad3/c-MYC pathway [81], EGFR signaling pathway [82] and Jak-STAT signaling pathway [83]. Quercetin can inhibit tumor cell proliferation by inducing apoptosis in glioma [84], melanoma [85], and bladder cancer [86]. It can also inhibit tumor growth by inhibiting angiogenesis [87] or inhibit cell proliferation by inducing autophagy of tumor cells, such as human glioma cells, primary lymphoma cells, ovarian cancer cells, human leukemia cells, etc., such as [88-90]. Cancer cell metastasis is associated with cell migration and invasion and is a major cause of death and poor prognosis in patients [91]. Studies have shown that quercetin can bate the migration and invasion of rectal carcinoma and head and neck squamous cell carcinoma [92-93], lung cancer cells [94], pancreatic cancer cells [95], melanoma cells [96], thereby preventing the metastasis of cancer cells. Not only that, but quercetin can also improve sensitivity and reduce the side effects of radiation and chemotherapy [97].

Quercetin not only exerts an antitumor effect alone, but also exerts an antitumor effect in combination with imatinib, inhibits the activity of cancer cells [98], improves the sensibility of cancer cells to radiotherapy [99], and affects p-glycoprotein (P-gp) to reverse the drug tolerance of cancer cells [100]. The combination of quercetin and gemcitabine also improves outcomes in patients with pancreatic cancer [101]. These results show that quercetin has a good anticancer effect and is used in clinical treatment in combination with other drugs as anticancer sensitizers.

Kaempferol is widely used in medicine for its antioxidant, anti-inflammatory, antibacterial, anticancer, heart, and neuroprotective effects [102]. Kaempferol protects normal cells, prevents normal cells from transforming into tumor cells [103], and inhibits the proliferation of cancer cells by increasing the apoptosis of tumor cells [104]. Not only that, but kaempferol also induces cell death by inducing autophagy [105]. Kaempferol alleviates oxidized LDL ox-LDL-induced apoptosis by participating in various signaling pathways, such as inhibition of AKT signaling-induced cell death in ovarian cancer [106], inhibition of the PI3K/Akt/mTOR pathway, thereby alleviating oxidized LDL ox-LDL-induced apoptosis [107], and mitigating palmitate acid-induced lipid storage, endoplasmic reticulum stress, and pancreatic β cell dysfunction through mTOR-mediated autophagy pathway [108]. It exerts antitumor effects by regulating TNF-α, NF κB pathway [109]. In addition, kaempferol inhibits the conversion of aromatic hydrocarbon receptor AHR [110] and also inhibits cell proliferation in breast and prostate cancer [111]. Kaempferol can reverse tumor resistance by inhibiting P-gp [112-113] and exert adjunctive antitumor effects by increasing the sensitivity of cisplatin and 5-fluorouracil [114-115]. Furthermore, epidemiological studies suggest a correlation between kaempferol intake and the danger of multi-tumors [116-117]. Kaempferol prevents tumor invasion and metastasis. High levels of MMP3 are associated with stronger tumor aggressiveness and poor prognosis. Most flavonoids do not affect MMP3 secretion, but kaempferol significantly inhibits the migration of human breast cancer cells and is dose-dependent [118]. These results reveal that kaempferol plays a crucial role in tumor apoptosis, autophagy, proliferation, drug resistance, invasion, and migration.

Quercetin and kaempferol have many of the same targets, such as AKT1, AHR, BCL2, NR4A1, etc., and AHR can inhibit the development of colorectal cancer [119-121]. In rhabdomyosarcoma, they can inhibit NR4A1-dependent reverse transcription activation and downstream target genes and inhibit the growth of xenograft rhabdomyosarcoma in nude mice [122]. This suggests that in colorectal cancer, the combination of the two may promote the expression of AHR, inhibit the expression of NR4A1, and achieve the synergistic effect of inhibiting CRC. However, due to the interaction between target proteins, after screening by the Cytoscape plugin, most of the target proteins, such as NR4A1, BCL2, CDK1, NOS2, and MMP1, etc., are excluded, but the role of these proteins cannot be ignored, and in the screening process, there are multiple mechanisms of action between protein targets and genes, which cannot be “consistent”.

Although quercetin and kaempferol alone can enhance anti-cancer, and non-specific immunity, and reduce the immunosuppressive effects, side effects, and precancerous lesions, the combination of quercetin and kaempferol can also increase CRC cells [26] of cell cycle arrest, but whether their combined use reduces the therapeutic burden is unclear. We study this through a bioinformatics database to understand the molecular and genetic pathways of these two drugs in vivo. This study preliminarily explored the molecular mechanisms of these two components in the treatment of CRC, laying a foundation for further studies of their biomolecular mechanisms. From the results of the gene function analysis, it can be seen that the combination of quercetin and kaempferol is effective in the treatment of colorectal cancer because they have similar roles in multiple molecular pathways and gene functions, for instance inhibiting the expression of MMP3 [123-128]. However, in CRC, whether kaempferol can inhibit MMP3 activity, thereby affecting the development, progression, and migration of CRC? Although there is no immediate proof at present, we expect to verify it through pharmacological experiments later, and based on the molecular docking prediction map of kaempferol and MMP3, we also doubt that kaempferol might directly interact with MMP3, which could become a novel finding in this study.

Moreover, we applied Lasso-Cox regression analysis to build a prognostic analysis model of Lasso-Cox regression with risk scores for evaluating patients in CRC two-component therapy or planned two-component therapy. The risk score and prognostic survival rate of patients are obtained through this model. We can predict risk scores and prognostic survival rates through the nomogram of this prognostic model to conduct clinical therapeutic strategies. Through bioinformatics analysis, it was found that three prognostic genes of the combination of quercetin and kaempferol in the treatment of CRC were positively correlated with immune infiltration. The expression differences of these three prognostic genes were found by plotting the expression levels of immunohistochemical genes in the cancer group and the normal group in the HPA database. Accordingly, it is reasonable to believe that this predictive model is accurate. This study shows that the same regimen of quercetin with kaempferol may have better results in high-risk patients; that is, immune cell levels may be higher in high-risk patients than in low-risk patients, and the expression of immune checkpoints is more pronounced. Thus, the risk values obtained from this model can guide clinical immunotherapy. In the CRC tumor microenvironment, patients with deeper immune invasion have higher survival rates. In practice, the combination of quercetin and kaempferol may help improve the effectiveness of immunotherapy in patients with CRC.

Conclusions

In conclusion, this study analyzed the molecular mechanism of quercetin and kaempferol in the therapy of CRC and confirmed the synergistic effect of quercetin and kaempferol in the therapy of CRC. Formerly studies showed that these two components have similar effects, there is also overlap between targets, and both can inhibit the role of P-gp to reverse tumor drug resistance and inhibit some signal pathways, such as P13K/AKT signaling pathway, EGFR signaling pathway, mTOR signaling pathway, MAPK signaling pathway, Jak-STAT signaling pathway, and NF-κB signaling pathway et al. Thus, there is reason to believe that quercetin and kaempferol can exert a synergistic effect at the level of molecular pathways.

In this study, we applied modern network pharmacology combined with bioinformatics to analyze the molecular mechanism of quercetin and kaempferol in the treatment of CRC, innovatively established a Lasso-Cox prognostic regression model, obtained three prognosis-related genes, and analyzed the immunotherapy effect of the two drugs. These studies can help to guide clinical application. We hope to next verify the conclusions of this study by using pharmacology, molecular biology, and cell function experiments in later studies to make the conclusions more accurate and dependable. Some deficiencies still exist in this study and are expected to be supplemented in subsequent studies.

Abbreviations

CRC: Colorectal cancer; TCM: Traditional Chinese medicine; TCMSP: Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform; TCGA: The Cancer Genome Atlas; GEO: Gene Expression Omnibus; RNA-seq: RNA sequencing; ICI: Immune checkpoint inhibitors; NAFLD: Nonalcoholic fatty liver disease; CHM : Chinese herbal medicine; GSEA: Gene Set Enrichment Analysis; OMIM: Online Mendelian Inheritance in Man; PharmGkb: Pharmacogenetics and Pharmacogenomics Knowledge Base; DisGeNET: a database of gene-disease associations; TTD: Therapeutic Target Database; HPA: The Human Protein Atlas; PPI: protein-protein interaction; STRING: search tool for retrieval; KEGG: Kyoto Encyclopedia of Genes and Genomes; GO: Gene ontology; ROC: receiver operating characteristic curve; IOBR: Multi-Omics Immuno-Oncology Biological Research; ESTIMATE: Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data; coef: coefficient; TIMER2.0: Tumor Immune Estimation Resource; AUC: Area Under Curve; CAM: Complementary and Alternative Medicine; EGCG: epigallocatechin-3-gallate; GTE: Green tea extract; GLE: Ganodorma lucidum extract; NSCLC: Non-Small Cell Lung Carcinoma; P-gp: P-glycoprotein.

Supplementary figure 2.

Attachment

Supplementary figure 3.

Attachment

Supplementary tables.

Attachment Acknowledgements

We are grateful to the participants for their consent and cooperation in the study. This work was supported by the Basic and Applied Basic Research Foundation of Guangdong Province (No. 2022A1515220217), Medical Scientific Research Foundation of Guangdong Province (No. A2023216; A2022524), Science and Technology Program of Guangzhou (No. 202201010840; 202201010810; 202102080532; 202002030032; 202002020023), Health Commission Program of Guangzhou (20212A010025; 20201A010085), Science and Technology Project of Panyu, Guangzhou (2022-Z04-009; 2022-Z04-090; 2022-Z04-072; 2021-Z04-053; 2020-Z04-026; 2019-Z04-02), Scientific Research project of Guangzhou Panyu Central Hospital (No. 2022Y002; 2021Y004; 2021Y002; ). There was no role of the funding bodies in the design of the study, in the collection, analysis, and interpretation of data, or in writing the manuscript.

Author contributions

Conceptualization, Chenqiong Gu; methodology, Chenqiong Gu; software, Chenqiong Gu and Lindong Tang; formal analysis, Chenqiong Gu and Yinghui Hao; validation, Chenqiong Gu, Jian Shen, and Shanshan Dong; investigation, Chenqiong Gu and Jian Shen; data curation, Chenqiong Gu; writing—original draft preparation, Chenqiong Gu; writing—review and editing, Chenqiong Gu, Jinhua He, and Li Yu; visualization, Chenqiong Gu; supervision, Chenqiong Gu and Jinhua; project administration, Chenqiong Gu and Jinhua He; All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

We acquired the original design data as well as the processed data in the compressed package, including all the files, code, and pictures.

Competing Interests

The authors have declared that no competing interest exists.

References

1. Yelorda KL, Fu SJ, Owens DK. Analysis of Survival Among Adults With Early-Onset Colorectal Cancer. JAMA Netw Open. 2021;4:e2112878

2. Sakin A, Sahin S, Atci MM, Yasar N, Geredeli C, Aribal S. et al. Factors affecting survival in patients with isolated liver-metastatic colorectal cancer treated with local ablative or surgical treatments for liver metastasis. J buon. 2019;24:1801-8

3. Li C, Sun YD, Yu GY, Cui JR, Lou Z, Zhang H. et al. Integrated Omics of Metastatic Colorectal Cancer. Cancer Cell. 2020;38:734-47.e9

4. Akagi T, Inomata M. Essential advances in surgical and adjuvant therapies for colorectal cancer 2018-2019. Ann Gastroenterol Surg. 2020;4:39-46

5. Lacouture M, Sibaud V. Toxic Side Effects of Targeted Therapies and Immunotherapies Affecting the Skin, Oral Mucosa, Hair, and Nails. Am J Clin Dermatol. 2018;19:31-9

6. Andre T, Amonkar M, Norquist JM, Shiu KK, Kim TW, Jensen BV. et al. Health-related quality of life in patients with microsatellite instability-high or mismatch repair deficient metastatic colorectal cancer treated with first-line pembrolizumab versus chemotherapy (KEYNOTE-177): an open-label, randomised, phase 3 trial. Lancet Oncol. 2021;22:665-77

7. Gao L, Wang XD, Niu YY, Duan DD, Yang X, Hao J. et al. Molecular targets of Chinese herbs: a clinical study of hepatoma based on network pharmacology. Sci Rep. 2016;6:24944

8. Liu J, Hao J, Niu Y, Wu X. Network pharmacology-based and clinically relevant prediction of active ingredients and potential targets of Chinese herbs on stage IV lung adenocarcinoma patients. J Cancer Res Clin Oncol. 2021;147:2079-92

9. Zahedi M, Ghiasvand R, Feizi A, Asgari G, Darvish L. Does Quercetin Improve Cardiovascular Risk factors and Inflammatory Biomarkers in Women with Type 2 Diabetes: A Double-blind Randomized Controlled Clinical Trial. Int J Prev Med. 2013;4:777-85

10. Van den Eynde MDG, Geleijnse JM, Scheijen J, Hanssen NMJ, Dower JI, Afman LA. et al. Quercetin, but Not Epicatechin, Decreases Plasma Concentrations of Methylglyoxal in Adults in a Randomized, Double-Blind, Placebo-Controlled, Crossover Trial with Pure Flavonoids. J Nutr. 2018;148:1911-6

11. S. A H, A A, T. O M, T. A A. Quercetin Dampens Postprandial Hyperglycemia in Type 2 Diabetic Patients Challenged with Carbohydrates Load. International Journal of Diabetes Research. 2012;1:32-5

12. Mazloom Z, Abdollahzadeh SM, Dabbaghmanesh MH, Rezaianzadeh A. The Effect of Quercetin Supplementation on Oxidative Stress, Glycemic Control, Lipid Profile and Insulin Resistance in Type 2 Diabetes: A Randomized Clinical Trial. Journal of health sciences and surveillance system. 2014;2:8-14

13. Dhanya R. Quercetin for managing type 2 diabetes and its complications, an insight into multitarget therapy. Biomed Pharmacother. 2022;146:112560

14. Lu TM, Chiu HF, Shen YC, Chung CC, Venkatakrishnan K, Wang CK. Hypocholesterolemic Efficacy of Quercetin Rich Onion Juice in Healthy Mild Hypercholesterolemic Adults: A Pilot Study. Plant Foods Hum Nutr. 2015;70:395-400

15. Prysyazhnyuk VP, Voloshyn OI. Effects of comprehensive treatment with quercetin administration on biochemical blood parameters and pro-and anti-inflammatory cytokines in nonalcoholic fatty liver disease patients. The Pharma Innovation Journal. 2017;6:386-9

16. Salehi B, Machin L, Monzote L, Sharifi-Rad J, Ezzat SM, Salem MA. et al. Therapeutic Potential of Quercetin: New Insights and Perspectives for Human Health. ACS Omega. 2020;5:11849-72

17. Ji Y, Li L, Ma YX, Li WT, Li L, Zhu HZ. et al. Quercetin inhibits growth of hepatocellular carcinoma by apoptosis induction in part via autophagy stimulation in mice. J Nutr Biochem. 2019;69:108-19

18. Li S, Pei Y, Wang W, Liu F, Zheng K, Zhang X. Quercetin suppresses the proliferation and metastasis of metastatic osteosarcoma cells by inhibiting parathyroid hormone receptor 1. Biomed Pharmacother. 2019;114:108839

19. Cheng S, Gao N, Zhang Z, Chen G, Budhraja A, Ke Z. et al. Quercetin induces tumor-selective apoptosis through downregulation of Mcl-1 and activation of Bax. Clin Cancer Res. 2010;16:5679-91

20. Maurya AK, Vinayak M. PI-103 and Quercetin Attenuate PI3K-AKT Signaling Pathway in T- Cell Lymphoma Exposed to Hydrogen Peroxide. PLoS One. 2016;11:e0160686

21. Kim HS, Wannatung T, Lee S, Yang WK, Chung SH, Lim JS. et al. Quercetin enhances hypoxia-mediated apoptosis via direct inhibition of AMPK activity in HCT116 colon cancer. Apoptosis. 2012;17:938-49

22. Kim GT, Lee SH, Kim YM. Quercetin Regulates Sestrin 2-AMPK-mTOR Signaling Pathway and Induces Apoptosis via Increased Intracellular ROS in HCT116 Colon Cancer Cells. J Cancer Prev. 2013;18:264-70

23. Huang S, Zhang Z, Li W, Kong F, Yi P, Huang J. et al. Network Pharmacology-Based Prediction and Verification of the Active Ingredients and Potential Targets of Zuojinwan for Treating Colorectal Cancer. Drug Des Devel Ther. 2020;14:2725-40

24. Imran M, Salehi B, Sharifi-Rad J, Aslam Gondal T, Saeed F, Imran A. et al. Kaempferol: A Key Emphasis to Its Anticancer Potential. Molecules. 2019 24

25. Zhou M, Li J, Luo D, Zhang H, Yu Z, Chen Y. et al. Network Pharmacology and Molecular Docking-Based Investigation: Prunus mume Against Colorectal Cancer via Silencing RelA Expression. Front Pharmacol. 2021;12:761980

26. Carmona SJ, López S, Abia R, Rodríguez-Arcos R, Araujo AJ, Bejarano RG. et al. Combination of quercetin and kaempferol enhances in vitro cytotoxicity on human colon cancer (HCT-116) cells. 2014.

27. Li S. Network Systems Underlying Traditional Chinese Medicine Syndrome and Herb Formula. Current Bioinformatics. 2009;4:188-96

28. Li H, Zhao L, Zhang B, Jiang Y, Wang X, Guo Y. et al. A network pharmacology approach to determine active compounds and action mechanisms of ge-gen-qin-lian decoction for treatment of type 2 diabetes. Evid Based Complement Alternat Med. 2014;2014:495840

29. Ru J, Li P, Wang J, Zhou W, Li B, Huang C. et al. TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J Cheminform. 2014;6:13

30. Zhang GB, Li QY, Chen QL, Su SB. Network pharmacology: a new approach for chinese herbal medicine research. Evid Based Complement Alternat Med. 2013;2013:621423

31. Fu B, Li Y, Shi X, Liu P, Zhang Y, Tian H. NCAPG Promotes Pulmonary Artery Smooth Muscle Cell Proliferation as a Promising Therapeutic Target of Idiopathic Pulmonary Hypertension: Bioinformatics Analysis and Experiment Verification. Int J Mol Sci. 2022 23

32. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012;16:284-7

33. Gui J, Li H. Penalized Cox regression analysis in the high-dimensional and low-sample size settings, with applications to microarray gene expression data. Bioinformatics. 2005 21 13: 3001-8

34. Buergy D, Würschmidt F, Gkika E, Hörner-Rieber J, Knippen S, Gerum S. et al. Stereotactic body radiotherapy of adrenal metastases-A dose-finding study. Int J Cancer. 2022;151:412-21

35. Kohrt HE, Olshen RA, Bermas HR, Goodson WH, Wood DJ, Henry S. et al. New models and online calculator for predicting non-sentinel lymph node status in sentinel lymph node positive breast cancer patients. BMC Cancer. 2008;8:66

36. “1. Liu TT, Li R, Huo C, Li JP, Yao J, Ji XL, et al. Identification of CDK2-Related Immune Forecast Model and ceRNA in Lung Adenocarcinoma, a Pan-Cancer Analysis. Front Cell Dev Biol. 2021;9:682002 ”

37. Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q. et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020;48:W509-w14

38. Shen W, Song Z, Zhong X, Huang M, Shen D, Gao P. et al. Sangerbox: A comprehensive, interaction-friendly clinical bioinformatics analysis platform. iMeta. 2022 1

39. Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, Torres-Garcia W. et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun. 2013;4:2612

40. Zeng D, Ye Z, Shen R, Yu G, Wu J, Xiong Y. et al. IOBR: Multi-Omics Immuno-Oncology Biological Research to Decode Tumor Microenvironment and Signatures. Front Immunol. 2021;12:687975

41. Xiao Y, Zhang G, Wang L, Liang M. Exploration and validation of a combined immune and metabolism gene signature for prognosis prediction of colorectal cancer. Front Endocrinol (Lausanne). 2022;13:1069528

42. Allen WL, Dunne PD, McDade S, Scanlon E, Loughrey M, Coleman H. et al. Transcriptional subtyping and CD8 immunohistochemistry identifies poor prognosis stage II/III colorectal cancer patients who benefit from adjuvant chemotherapy. JCO Precis Oncol. 2018. 2018

43. Marisa L, de Reyniès A, Duval A, Selves J, Gaub MP, Vescovo L. et al. Gene expression classification of colon cancer into molecular subtypes: characterization, validation, and prognostic value. PLoS Med. 2013;10:e1001453

44. Burbidge O, Pastok MW, Hodder SL, Zenkevičiūtė G, Noble MEM, Endicott JA, Itzhaki LS. Nanobodies restore stability to cancer-associated mutants of tumor suppressor protein p16INK4a. bioRxiv. 2021

45. Sillen M, Miyata T, Vaughan DE, Strelkov SV, Declerck PJ. Structural Insight into the Two-Step Mechanism of PAI-1 Inhibition by Small Molecule TM5484. International Journal of Molecular Sciences. 2021;22:1482

46. Durham TB, Klimkowski VJ, Rito CJ, Marimuthu J, Toth JL, Liu C. et al. Identification of Potent and Selective Hydantoin Inhibitors of Aggrecanase-1 and Aggrecanase-2 That Are Efficacious in Both Chemical and Surgical Models of Osteoarthritis. Journal of Medicinal Chemistry. 2014;57:10476-85

47. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA. et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545-50

48. Hu XQ, Sun Y, Lau E, Zhao M, Su SB. Advances in Synergistic Combinations of Chinese Herbal Medicine for the Treatment of Cancer. Curr Cancer Drug Targets. 2016;16:346-56

49. Luan X, Zhang LJ, Li XQ, Rahman K, Zhang WD. Compound-based Chinese medicine formula: From discovery to compatibility mechanism. Journal of Ethnopharmacology. 2020;254:112687 -

50. Wang Q, Jiao L, Wang S, Chen P, Bi L, Zhou D. et al. Maintenance Chemotherapy With Chinese Herb Medicine Formulas vs. With Placebo in Patients With Advanced Non-small Cell Lung Cancer After First-Line Chemotherapy: A Multicenter, Randomized, Double-Blind Trial. Front Pharmacol. 2018;9:1233

51. Tang YC, Zhang Y, Zhou J, Zhi Q, Wu MY, Gong FR. et al. Ginsenoside Rg3 targets cancer stem cells and tumor angiogenesis to inhibit colorectal cancer progression in vivo. Int J Oncol. 2018;52:127-38

52. Zhang Q, Chen X, Luo Y, Ren H, Qiao T. Fuzi Enhances Anti-Tumor Efficacy of Radiotherapy on Lung Cancer. J Cancer. 2017;8:3945-51

53. Liu X, Tian S, Liu M, Jian L, Zhao L. Wogonin inhibits the proliferation and invasion, and induces the apoptosis of HepG2 and Bel7402 HCC cells through NF-κB/Bcl-2, EGFR and EGFR downstream ERK/AKT signaling. Int J Mol Med. 2016;38:1250-6

54. Foungbé S, Kouassi G, Kablan JB, Marcy R. Study of Costus lucanusianus: plant juice, fraction combinations and pharmacologic estimation of natural product total activity. J Ethnopharmacol. 1991;33:221-6

55. Schuster BG. A new integrated program for natural product development and the value of an ethnomedical approach. J Altern Complement Med. 2001;7(Suppl 1):S61-72

56. Zhou SC, liu LH. Effects of Resveratrol Combined with Curcumin on Apoptosis of Liver Cancer Cells. Journal of Medical Molecular Biology. 2018;15:156-160

57. Zhang FH, Yin SM, Niu GM, Yang GL, Ling YW, Li R, Xu JY.Guide of China Medicine. Effects of Quercetin and Resveratrol combined Reversing Drug Resistance of Leukemia Cell Line HL-60-ADM. 2014; 12.

58. Zhang P, Zhang X. Stimulatory effects of curcumin and quercetin on posttranslational modifications of p53 during lung carcinogenesis. Human & experimental toxicology. 2018;37:618-625

59. Song SY, Park JH, Park SJ, Kang IC, Yoo HS. Synergistic Effect of HAD-B1 and Afatinib Against Gefitinib Resistance of Non-Small Cell Lung Cancer. Integr Cancer Ther. 2022;21:15347354221144311

60. Gao S, Li X, Ding X, Qi W, Yang Q. Cepharanthine Induces Autophagy, Apoptosis and Cell Cycle Arrest in Breast Cancer Cells. Cell Physiol Biochem. 2017;41:1633-48

61. Mukhtar H, Ahmad N. Tea polyphenols: prevention of cancer and optimizing health. Am J Clin Nutr. 2000;71:1698S-702S discussion 703S-4S

62. Yang CS, Maliakal P, Meng X. Inhibition of carcinogenesis by tea. Annu Rev Pharmacol Toxicol. 2002;42:25-54

63. Sohretoglu D, Huang S. Ganoderma lucidum Polysaccharides as An Anti-cancer Agent. Anticancer Agents Med Chem. 2018;18:667-74

64. Chang CJ, Lin CS, Lu CC, Martel J, Ko YF, Ojcius DM. et al. Ganoderma lucidum reduces obesity in mice by modulating the composition of the gut microbiota. Nat Commun. 2015;6:7489

65. Seweryn E, Ziała A, Gamian A. Health-Promoting of Polysaccharides Extracted from Ganoderma lucidum. Nutrients. 2021 13

66. Guo C, Guo D, Fang L, Sang T, Wu J, Guo C. et al. Ganoderma lucidum polysaccharide modulates gut microbiota and immune cell function to inhibit inflammation and tumorigenesis in colon. Carbohydr Polym. 2021;267:118231

67. Ren ZL, Wang CD, Wang T, Ding H, Zhou M, Yang N. et al. Ganoderma lucidum extract ameliorates MPTP-induced parkinsonism and protects dopaminergic neurons from oxidative stress via regulating mitochondrial function, autophagy, and apoptosis. Acta pharmacologica Sinica. 2019;40:441-450

68. Chan SW, Tomlinson B, Chan P, Lam CWK. The beneficial effects of Ganoderma lucidum on cardiovascular and metabolic disease risk. Pharm Biol. 2021;59:1161-71

69. Thyagarajan A, Zhu J, Sliva D. Combined effect of green tea and Ganoderma lucidum on invasive behavior of breast cancer cells. Int J Oncol. 2007;30:963-9

70. Zhuang SR, Chiu HF, Chen SL, Tsai JH, Lee MY, Lee HS. et al. Effects of a Chinese medical herbs complex on cellular immunity and toxicity-related conditions of breast cancer patients. Br J Nutr. 2012;107:712-8

71. Guo L, Bai SP, Zhao L, Wang XH. Astragalus polysaccharide injection integrated with vinorelbine and cisplatin for patients with advanced non-small cell lung cancer: effects on quality of life and survival. Med Oncol. 2012;29:1656-62

72. Yap KY, See CS, Chan A. Clinically-relevant chemotherapy interactions with complementary and alternative medicines in patients with cancer. Recent Pat Food Nutr Agric. 2010;2:12-55

73. Li S, Zhang B. Traditional Chinese medicine network pharmacology: theory, methodology and application. Chin J Nat Med. 2013;11:110-20

74. Wang QQ, Zhou YC, Zhou Ge YJ, Qin G, Yin TF, Zhao DY. et al. Comprehensive proteomic signature and identification of CDKN2A as a promising prognostic biomarker and therapeutic target of colorectal cancer. World J Clin Cases. 2022;10:7686-97

75. Shi WK, Li YH, Bai XS, Lin GL. The Cell Cycle-Associated Protein CDKN2A May Promotes Colorectal Cancer Cell Metastasis by Inducing Epithelial-Mesenchymal Transition. Front Oncol. 2022;12:834235

76. Zeng C, Chen Y. HTR1D, TIMP1, SERPINE1, MMP3 and CNR2 affect the survival of patients with colon adenocarcinoma. Oncol Lett. 2019;18:2448-54

77. Guo H, Wang Y, Gou L, Wang X, Tang Y, Wang X. A novel prognostic model based on urea cycle-related gene signature for colorectal cancer. Front Surg. 2022;9:1027655

78. Yuan W, Cai W, Huang X, Peng S. Prognostic value of immune scores in the microenvironment of colorectal cancer. Oncol Lett. 2020;20:256

79. Luo Q, Zhou P, Chang S, Huang Z, Zeng X. Characterization of butyrate-metabolism in colorectal cancer to guide clinical treatment. Sci Rep. 2023;13:5106

80. Lim W, Yang C, Park S, Bazer FW, Song G. Inhibitory Effects of Quercetin on Progression of Human Choriocarcinoma Cells Are Mediated Through PI3K/AKT and MAPK Signal Transduction Cascades. Journal of Cellular Physiology. 2017;232:1428-40

81. An CY, Xie G, Tang XW, Yang XD, Zhang CQ. Effects of quercetin on proliferation and invasion of colon cancer LOVO cells and expression of carcinoembryonic antigen CEA. Chin J Clin Pharmacol Ter. 2013;18:24-29

82. Firdous AB, Sharmila G, Balakrishnan S, RajaSingh P, Suganya S, Srinivasan N. et al. Quercetin, a natural dietary flavonoid, acts as a chemopreventive agent against prostate cancer in an in vivo model by inhibiting the EGFR signaling pathway. Food Funct. 2014;5:2632-45

83. Granato M, Rizzello C, Gilardini Montani MS, Cuomo L, Vitillo M, Santarelli R. et al. Quercetin induces apoptosis and autophagy in primary effusion lymphoma cells by inhibiting PI3K/AKT/mTOR and STAT3 signaling pathways. J Nutr Biochem. 2017;41:124-36

84. Liu Y, Tang ZG, Lin Y, Qu XG, Lv W, Wang GB. et al. Effects of quercetin on proliferation and migration of human glioblastoma U251 cells. Biomed Pharmacother. 2017;92:33-8

85. Zhu JM, Zhang R, Lu XY, Song HJ, Xue H, Liu FM. Effect and mechanism of Quercetin on proliferation of cervical cancer HeLa cells. Chinese Journal of Cancer Prevention and Treatment. 2014;21:1073-1077

86. Chen SN. Study of inhibitory effect and influence on cell cycle of quercetin and its 7-OH derivative bladder carcinoma cell.Lingnan Modern Clinics in Surgery.2016; 16: 204-209

87. Yang F, Jiang X, Song L, Wang H, Mei Z, Xu Z. et al. Quercetin inhibits angiogenesis through thrombospondin-1 upregulation to antagonize human prostate cancer PC-3 cell growth in vitro and in vivo. Oncol Rep. 2016;35:1602-10

88. He Y, Cao X, Guo P, Li X, Shang H, Liu J. et al. Quercetin induces autophagy via FOXO1-dependent pathways and autophagy suppression enhances quercetin-induced apoptosis in PASMCs in hypoxia. Free Radic Biol Med. 2017;103:165-76

89. Liu Y, Gong W, Yang ZY, Zhou XS, Gong C, Zhang TR. et al. Quercetin induces protective autophagy and apoptosis through ER stress via the p-STAT3/Bcl-2 axis in ovarian cancer. Apoptosis. 2017;22:544-57

90. Chang JL, Chow JM, Chang JH, Wen YC, Lin YW, Yang SF. et al. Quercetin simultaneously induces G(0) /G(1) -phase arrest and caspase-mediated crosstalk between apoptosis and autophagy in human leukemia HL-60 cells. Environ Toxicol. 2017;32:1857-68

91. Pezzuto A, Piraino A, Mariotta S. Lung cancer and concurrent or sequential lymphoma: Two case reports with hypersensitivity to bevacizumab and a review of the literature. Oncol Lett. 2015;9:604-8

92. Kee JY, Han YH, Kim DS, Mun JG, Park J, Jeong MY. et al. Inhibitory effect of quercetin on colorectal lung metastasis through inducing apoptosis, and suppression of metastatic ability. Phytomedicine: international journal of phytotherapy and phytopharmacology. 2016;23:1680-1690

93. Chan CY, Lien C.H, Lee MF, Huang CY. Quercetin suppresses cellular migration and invasion in human head and neck squamous cell carcinoma (HNSCC). BioMedicine. 2016;6:15

94. Chang JH, Lai SL, Chen WS, Hung WY, Chow JM, Hsiao M. et al. Quercetin suppresses the metastatic ability of lung cancer through inhibiting Snail-dependent Akt activation and Snail-independent ADAM9 expression pathways. Biochimica et biophysica acta. Molecular cell research. 2017;1864:1746-1758

95. Yu D, Ye T, Xiang Y, Shi Z, Zhang J, Lou B. et al. Quercetin inhibits epithelial-mesenchymal transition, decreases invasiveness and metastasis, and reverses IL-6 induced epithelial-mesenchymal transition, expression of MMP by inhibiting STAT3 signaling in pancreatic cancer cells. OncoTargets and therapy. 2017;10:4719-4729

96. Cao H. et al. Quercetin inhibits HGF/c-Met signaling and HGF-stimulated melanoma cell migration and invasion. Molecular Cancer. 2015

97. Brito AF, Ribeiro M, Abrantes AM, Pires AS, Teixo RJ, Tralhão JG. et al. Quercetin in Cancer Treatment, Alone or in Combination with Conventional Therapeutics?. Curr Med Chem. 2015;22:3025-39

98. Chen DQ. Effect of combined use of quercetin and imatinib on proliferation and apoptosis of SNU-1 cells in vitro[J]. Journal of Gansu University of Chinese Medicine. 2017;34:27-31

99. Wang F. The study of effect of quercetin combined X-ray on U-2OS cell growth and apoptosis. Pharmacology and Clinics of Chinese Materia Medica. 2015;31:30-33

100. Li Cl, Cheng D, Sun ZQ, Liu Y.Study of drug resistance by the Quercetin in human colorectal cancer SW480/OXP cells.Chin J Gastroenterol Hepato.2017; 26:551-554

101. Cao C, Sun L, Mo W, Sun L, Luo J, Yang Z. et al. Quercetin Mediates β-Catenin in Pancreatic Cancer Stem-Like Cells. Pancreas. 2015;44:1334-9

102. Calderón-Montaño JM, Burgos-Morón E, Pérez-Guerrero C, López-Lázaro M. A review on the dietary flavonoid kaempferol. Mini Rev Med Chem. 2011;11:298-344

103. Lee S, Kim YJ, Kwon S, Lee Y, Choi SY, Park J. et al. Inhibitory effects of flavonoids on TNF-alpha-induced IL-8 gene expression in HEK 293 cells. BMB Rep. 2009;42:265-70

104. Seibert H, Maser E, Schweda K, Seibert S, Gülden M. Cytoprotective activity against peroxide-induced oxidative damage and cytotoxicity of flavonoids in C6 rat glioma cells. Food Chem Toxicol. 2011;49:2398-407

105. Kim TW, Lee SY, Kim M, Cheon C, Ko SG. Kaempferol induces autophagic cell death via IRE1-JNK-CHOP pathway and inhibition of G9a in gastric cancer cells. Cell Death Dis. 2018;9:875

106. El-Kott AF, Shati AA, Al-Kahtani MA, Alharbi SA. Kaempferol Induces Cell Death in A2780 Ovarian Cancer Cells and Increases Their Sensitivity to Cisplatin by Activation of Cytotoxic Endoplasmic Reticulum-Mediated Autophagy and Inhibition of Protein Kinase B. Folia Biol (Praha). 2020;66:36-46

107. Nirmala P, Ramanathan M. Effect of kaempferol on lipid peroxidation and antioxidant status in 1,2-dimethyl hydrazine induced colorectal carcinoma in rats. Eur J Pharmacol. 2011;654:75-9

108. Varshney R, Varshney R, Mishra R, Gupta S, Sircar D, Roy P. Kaempferol alleviates palmitic acid-induced lipid stores, endoplasmic reticulum stress and pancreatic β-cell dysfunction through AMPK/mTOR-mediated lipophagy. J Nutr Biochem. 2018;57:212-27

109. Chen SS, Michael A, Butler-Manuel SA. Advances in the treatment of ovarian cancer: a potential role of antiinflammatory phytochemicals. Discov Med. 2012;13:7-17

110. Mukai R, Satsu H, Shimizu M, Ashida H. Inhibition of P-glycoprotein enhances the suppressive effect of kaempferol on transformation of the aryl hydrocarbon receptor. Biosci Biotechnol Biochem. 2009;73:1635-9

111. Wang H, Gao M, Wang J. Kaempferol inhibits cancer cell growth by antagonizing estrogen-related receptor α and γ activities. Cell Biol Int. 2013;37:1190-6

112. To KK, Yu L, Liu S, Fu J, Cho CH. Constitutive AhR activation leads to concomitant ABCG2-mediated multidrug resistance in cisplatin-resistant esophageal carcinoma cells. Mol Carcinog. 2012;51:449-64

113. Nakatsuma A, Fukami T, Suzuki T, Furuishi T, Tomono K, Hidaka S. Effects of kaempferol on the mechanisms of drug resistance in the human glioblastoma cell line T98G. Pharmazie. 2010;65:379-83

114. Luo H, Daddysman MK, Rankin GO, Jiang BH, Chen YC. Kaempferol enhances cisplatin's effect on ovarian cancer cells through promoting apoptosis caused by down regulation of cMyc. Cancer Cell Int. 2010;10:16

115. Ahmed F, Toume K, Ohtsuki T, Rahman M, Sadhu SK, Ishibashi M. Cryptolepine, isolated from Sida acuta, sensitizes human gastric adenocarcinoma cells to TRAIL-induced apoptosis. Phytother Res. 2011;25:147-50

116. Geybels MS, Verhage BA, Arts IC, van Schooten FJ, Goldbohm RA, van den Brandt PA. Dietary flavonoid intake, black tea consumption, and risk of overall and advanced stage prostate cancer. Am J Epidemiol. 2013;177:1388-98

117. Bobe G, Albert PS, Sansbury LB, Lanza E, Schatzkin A, Colburn NH. et al. Interleukin-6 as a potential indicator for prevention of high-risk adenoma recurrence by dietary flavonols in the polyp prevention trial. Cancer Prev Res (Phila). 2010;3:764-75

118. Phromnoi K, Yodkeeree S, Anuchapreeda S, Limtrakul P. Inhibition of MMP-3 activity and invasion of the MDA-MB-231 human invasive breast carcinoma cell line by bioflavonoids. Acta Pharmacol Sin. 2009;30:1169-76

119. Díaz-Díaz CJ, Ronnekleiv-Kelly SM, Nukaya M, Geiger PG, Balbo S, Dator R. et al. The Aryl Hydrocarbon Receptor is a Repressor of Inflammation-associated Colorectal Tumorigenesis in Mouse. Ann Surg. 2016;264:429-36

120. Han H, Davidson LA, Hensel M, Yoon G, Landrock K, Allred C. et al. Loss of Aryl Hydrocarbon Receptor Promotes Colon Tumorigenesis in Apc(S580/+); Kras(G12D/+) Mice. Mol Cancer Res. 2021;19:771-83

121. Han H, Davidson LA, Fan YY, Goldsby JS, Yoon G, Jin UH. et al. Loss of aryl hydrocarbon receptor potentiates FoxM1 signaling to enhance self-renewal of colonic stem and progenitor cells. Embo j. 2020;39:e104319

122. Shrestha R, Mohankumar K, Martin G, Hailemariam A, Lee SO, Jin UH. et al. Flavonoids kaempferol and quercetin are nuclear receptor 4A1 (NR4A1, Nur77) ligands and inhibit rhabdomyosarcoma cell and tumor growth. J Exp Clin Cancer Res. 2021;40:392

123. Qiu L, Luo Y, Chen X. Quercetin attenuates mitochondrial dysfunction and biogenesis via upregulated AMPK/SIRT1 signaling pathway in OA rats. Biomed Pharmacother. 2018;103:1585-91

124. Zhu H, Xiong XG, Lu Y, Wu HC, Zhang ZH, Sun MJ. The mechanism of the anti-inflammatory effect of Oldenlandia diffusa on arthritis model rats: a quantitative proteomic and network pharmacologic study. Ann Transl Med. 2022;10:1098

125. Wang L, Han X, Li H, Lv C, Wang M. The ethyl acetate extract of Wenxia Changfu Formula inhibits the carcinogenesis of lung adenocarcinoma by regulating PI3K-AKT signaling pathway. Sci Rep. 2023;13:4715

126. Wang X, Tan Y, Liu F, Wang J, Liu F, Zhang Q. et al. Pharmacological network analysis of the functions and mechanism of kaempferol from Du Zhong in intervertebral disc degeneration (IDD). J Orthop Translat. 2023;39:135-46

127. Zhu J, Tang H, Zhang Z, Zhang Y, Qiu C, Zhang L. et al. Kaempferol slows intervertebral disc degeneration by modifying LPS-induced osteogenesis/adipogenesis imbalance and inflammation response in BMSCs. Int Immunopharmacol. 2017;43:236-42

128. Huang X, Pan Q, Mao Z, Wang P, Zhang R, Ma X. et al. Kaempferol inhibits interleukin-1β stimulated matrix metalloproteinases by suppressing the MAPK-associated ERK and P38 signaling pathways. Mol Med Rep. 2018;18:2697-704