Jun Rao1,2,3†
Xing Wang4†
Tanxiu Chen5†
Mingzi Mo6Conglong Xu7
So-Yi Chang2,3
Ssu-Wei Hsu2,3Xiaoqun Han6
Ching-Hsien Chen2,3*
Zhi Zheng4*- 1Jiangxi Cancer Hospital & Institute, The Second Affiliated Hospital of Nanchang Medical College, Jiangxi Clinical Research Center for Cancer, Nanchang, Jiangxi, China
- 2Division of Nephrology, Department of Internal Medicine, University of California, Davis, Davis, CA, United States
- 3Comprehensive Cancer Center, University of California, Davis, Sacramento, CA, United States
- 4Jiangxi Provincial People’s Hospital, The First Affiliated Hospital of Nanchang Medical College, Nanchang, Jiangxi, China
- 5Institute of Neurology and Department of Neurology, Jiangxi Academy of Clinical Medical Sciences, The First Affiliated Hospital, Jiangxi Medical College, Nanchang University, Nanchang, Jiangxi, China
- 6Yichun University, Yichun, China
- 7Jiangxi Jingde Traditional Chinese Medicine Co., Ltd., Jingdezhen, Jiangxi, China
Background: Hedysari Radix (HR), commonly known as Hong-Qi in Chinese, is a traditional Chinese herbal medicine recognized for possessing anti-inflammatory and anti-tumor properties. While the polysaccharides in HR have been extensively studied, other HR metabolites and their potential anti-tumor properties remain largely unknown.
Methods: We employed a multi-omics strategy integrating metabolomics, network analysis, proteomics, phosphoproteomics, and molecular docking to identify HR metabolites with anti-colorectal cancer (CRC) property and investigate underlying mechanisms.
Results: Using mass spectrometry-based metabolomics, we identified 1,292 metabolites across eight processed HR products. Key metabolites including medicarpin, formononetin, naringenin, and quercetin were validated via the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP). Notably, formononetin-derived metabolites were significantly enriched during HR processing. The metabolite-metabolite correlation analysis revealed key compounds such as flavonoids and formononetin. Subsequent network analysis combined with label-free data-independent acquisition (DIA) proteomics and phosphoproteomics in colon cancer cells identified 194 potential targets, 291 differentially expressed proteins, and 1,535 phosphorylated proteins that were regulated by formononetin. Cell-surface enzymes carbonic anhydrase IX (CA9) and membrane metalloendopeptidase (MME) were consistently identified in different analyses as key targets, and molecular docking results confirmed their strong binding to formononetin. Bioinformatics analyses further revealed significant enrichment of cancer-associated pathways, including PI3K-Akt, Hippo, HIF-1 signaling, and cholesterol metabolism upon formononetin treatment.
Conclusion: The findings provide novel insights into the HR metabolome and reveal the multi-targeting roles of formononetin in CRC development, laying the foundation for developing new CRC therapeutic strategies.
GRAPHICAL ABSTRACT |
1 Introduction
Colorectal cancer (CRC) is the third most common malignancy and the second leading cause of cancer-related deaths worldwide (Morgan et al., 2023). Global CRC incidence is projected to reach 3.2 million by 2040, posing a significant public health challenge. Current CRC treatments include surgery, radiotherapy, chemotherapy, targeted therapy, and immunotherapy. Patients with stage I-III CRC typically undergo radical surgery, and patients with metastatic stage IV CRC receive systemic therapies including conventional chemotherapeutic agents (e.g., 5-Fluorouracil, Irinotecan, Oxaliplatin), targeted therapeutic drugs (e.g., Aflibercept, Bevacizumab, Cetuximab, Panitumumab), and immune checkpoint inhibitors (e.g., Pembrolizumab, Nivolumab) (El Bali et al., 2021; Groelly et al., 2023; Miyamoto et al., 2023). However, these treatments often face problems such as immune suppression, chemotherapy resistance, toxic side effects, and high recurrence rates, making it necessary to develop more effective and less toxic alternatives.
Traditional Chinese medicine (TCM) has been used for thousands of years in China and other Asian countries to treat various diseases (Huang et al., 2024). As a complementary therapy for cancers including CRC, lung cancer, breast cancer, and ovarian cancer, TCM has been demonstrated to be able to enhance chemotherapy efficacy while reducing side effects and toxicity (Wan et al., 2022; Wei et al., 2023; Wang et al., 2024; Yu et al., 2024). TCM has shown benefits in regulating immune responses and protecting colonic mucosal barrier in stage II and III CRC patients, making it a potential long-term treatment strategy for late stage CRC (Feng et al., 2021; Wang et al., 2021; Wan et al., 2022). Hedysari Radix (Hedysarum polybotrys Hand.-Mazz. [Fabaceae], or HR), a well-known TCM, is widely used due to its diverse antioxidant, immunomodulatory, anti-inflammatory, and anti-tumor effects (Mo et al., 2022). While HR’s anti-tumor effects in lung, gastric, liver, and ovarian cancers are mostly credited to its polysaccharides, other bioactive compounds in HR like terpenoids and flavonoids also demonstrated significant anti-tumor properties and attracted broad research interests.
Advances in “omics” technologies, such as metabolomics and proteomics, have significantly contributed to deeper understanding of therapeutic mechanisms in many diseases including cancer (Zhu et al., 2022). Among them, metabolomics has been instrumental in identifying small molecule bioactive ingredients in TCM to elucidate their effects (Wang et al., 2017). In an integrated metabolomics and network pharmacology analysis (Qu et al., 2022), it was revealed that dandelion could suppress tumor progression in triple-negative breast cancer (TNBC) by modulating cell cycle signaling and metabolic pathways. Studies using proteomics and phosphoproteomics have focused on identifying molecular targets of TCMs to facilitate the development of novel cancer therapeutic agents (Wang et al., 2019; Song et al., 2022; Xia et al., 2022). One TCM ingredient, triptolide, has been demonstrated to be able to suppress CRC proliferation by regulating key proteins and pathways including PI3K/AKT and Hedgehog signaling (Song et al., 2022).
In this study, we integrated metabolomics, network analysis, proteomics, and phosphoproteomics to investigate the anti-tumor mechanisms of HR. We investigated eight heat-processed HR products and SW620 colon cancer cells to identify potentially key active metabolites and their molecular targets. Additionally, we analyzed metabolic changes during the heat-processing (called paozhi in Chinese) of HR to elucidate the therapeutic effects generated from the processing against CRC.
2 Materials and methods
2.1 Sample collection
Eight heat-processed HR products (labeled JZ, JPT, XFT, TQT, JJ, SFT, MDYY, and JD) were purchased from different manufacturers in China. Unprocessed JD HR was also obtained and designated SJD HR. Each HR product was prepared in triplicate, air-dried, flash-frozen in liquid nitrogen, and stored at −80 °C until analysis.
2.2 Metabolomics analysis
Targeted metabolomics analysis focusing on secondary metabolites was conducted following the standards outlined in the ConPhyMP statement (Supplementary Tables S1, S2) (Heinrich et al., 2022; Wang et al., 2022; Sun et al., 2023). Metabolites were extracted with 1,200 μL pre-cooled (−20 °C) 70% aqueous methanol containing 2-chlorophenylalanine (1 ppm; purity 98%) as an internal standard. After vortexing and centrifugation, the supernatant was filtered through a 0.22 μm microporous membrane (ANPEL, Shanghai, China) and transferred to an injection vial. Metabolomics analysis was performed using an ultra-high-performance liquid chromatography system (Nexera X2 UPLC, Shimadzu, Tokyo, Japan) coupled with an electrospray ionization-triple quadrupole-linear ion trap mass spectrometer (4500 QTRAP mass spectrometer, Applied Biosystems, Waltham, MA, USA). Chromatographic separation was achieved using an Agilent SB-C18 column (1.8 µm, 2.1 mm × 100 mm) with mobile solution A (0.1% formic acid in water) and solution B (0.1% formic acid in acetonitrile). Gradient elution was performed at a flow rate of 0.35 mL/min with an injection volume of 2 μL. The gradient program was: 95% A/5% B at 0 min, linearly ramped to 5% A/95% B at 9.0 min and held until 10.0 min, then returned to 95% A/5% B at 11.1 min and maintained until 14.0 min. The elution was directly introduced into the mass spectrometer with key operational parameters as follows: turbo spray ion source, source temperature 550 °C, ion spray voltage 5500 V (positive) and 4500 V (negative), ion source gas I 50 psi, gas II 60 psi, curtain gas 25 psi (Wang et al., 2022; Sun et al., 2023). The declustering potential and collision energy were individually optimized for each transition. Metabolite identification and quantification were conducted by searching the Metware Database (MWDB, MetWare Biological Co., Ltd., Wuhan, China). Chromatographic peaks were integrated and corrected using MultiQuant software (AB SCIEX, Framingham, MA, USA). SIMCA-P 18.0 software was used for principal component analysis (PCA), while significantly changed metabolites were determined using the orthogonal partial least squares discriminant analysis (OPLS-DA), followed by an independent t-test (SPSS 22.0 software) (Wan et al., 2022). VIP (Variable Importance in Projection) values from the OPLS-DA model reflected the contribution of each variable, with VIP values >1 indicating significance. The natural variation of the content of each secondary metabolite in HR was presented in a Z-score plot generated using R software (version 4.0.2) with three SD as the cutoff. Metabolite-metabolite correlation analysis was also performed using R software for Pearson’s product-moment correlation (Pearson’s r), p-values were calculated using the cor.test function.
2.3 Network analysis for target selection
The Simplified Molecular Input Line Entry System (SMILES) structure of formononetin was retrieved from PubChem (https://pubchem.ncbi.nlm.nih.gov/) to predict its potential molecular targets based on its structural properties and bioactivity profiles. Target prediction was performed in SwissTargetPrediction (http://www.swisstargetprediction.ch/), TargetNet (http://targetnet.scbdd.com/), and PharmMapper (http://lilab-ecust.cn/pharmmapper/submitfile.html), which integrate multiple computational approaches to achieve comprehensive coverage of plausible targets. Duplicate entries were removed to retain unique targets. Targets related to CRC were compiled by querying DisGeNET (https://www.disgenet.org/), GeneCards (https://www.genecards.org/), and OMIM (https://www.omim.org/) databases using the keyword “Colorectal cancer” to select genes with documented relevance to CRC pathogenesis. Middle targets between formononetin and CRC were identified using Venny 2.1.0 (https://bioinfogp.cnb.csic.es/tools/venny/) to generate a Venn diagram including potential therapeutic targets of formononetin in CRC. The overlapped targets were then selected and subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses via the Bioinformatics online platform (https://www.bioinformatics.com.cn/) to elucidate the biological processes and signaling pathways involved.
2.4 Cell culture and cell counting kit-8 (CCK-8) assay
Human SW620 colon cancer cells (ATCC, Manassas, VA) were cultured in RPMI-1640 medium supplemented with 10% fetal bovine serum at 37 °C in a 5% CO2 incubator. Cells were seeded in 96-well plates (2,000 cells/well) and treated with formononetin (catalog #F408902, Aladdin, Shanghai, China) at final concentrations of 0, 5, 10, 20, 40, and 80 μM for 48 h. Cell viability was assessed using CCK-8 reagent (catalog #DB884, Dojindo, Kumamoto, Japan). After 1 hour of incubation with CCK-8, absorbance at 450 nm was measured using a microplate spectrophotometer.
2.5 Proteomic analysis
Proteomic profiling was performed using a label-free DIA approach to identify protein changes in formononetin-treated SW620 cells as previously reported (Wang et al., 2020; Wan et al., 2022). Proteins were extracted, reduced, and digested using the Filter-Aided Sample Preparation (FASP) method and analyzed on an EASY-nLC 1,200 system (Thermo Fisher Scientific, San Jose, CA, USA) coupled to a Q Exactive HF mass spectrometer (Thermo Fisher Scientific, San Jose, CA, USA). Chromatographic separation was performed on a reversed-phase C18 column with a binary gradient mobile solutions comprising 0.1% formic acid in water (A) and 0.1% formic acid in 80% acetonitrile (B). The gradient program was: 6%–10% B (0–5 min), 10%–30% B (5–47 min), 30%–45% B (47–55 min), and 95% B (55–60 min). MS/MS analysis in DIA mode was set at a resolution of 30,000, 32 isolation windows (first mass: 200 m/z), automatic maximum injection time, an AGC target of 3e6, and normalized collision energy of 28% (Wang et al., 2020). DIA raw data was processed using Spectronaut 18.0 (Biognosys) for protein identification. Differentially expressed proteins were defined by a fold change >1.2 and p value <0.05. GO and KEGG enrichment analyses of these proteins were conducted using OmicsBean (http://www.omicsbean.cn). The raw MS data have been deposited to the ProteomeXchange Consortium (https://proteomecentral.proteomexchange.org) via the iProX partner repository with the dataset identifier PXD059226 (Chen et al., 2022).
2.6 Phosphoproteomic analysis
Phosphoproteomic analysis was performed on formononetin-treated SW620 cells using DIA-based methods. Protein extraction, digestion, and phosphopeptide enrichment were carried out using the High-Select™ Fe-NTA Kit (Thermo Fisher Scientific, San Jose, CA, USA) following the manufacturer’s instructions (Song et al., 2022). Briefly, lyophilized peptides were reconstituted in 200 μL of immobilized metal affinity chromatography (IMAC) binding/wash buffer. IMAC spin columns were activated by centrifugation (1,000 × g, 30 s) and equilibrated with two washes of the same buffer. After loading the samples, the resin was gently resuspended and incubated for 30 min with manual mixing at 10-min intervals. Following incubation, the columns were washed three times with binding/wash buffer and once with HPLC-grade water. Phosphopeptides were eluted twice with 100 μL elution buffer each time (1,000 × g, 30 s each), and the two elution were combined and dried under vacuum. The eluted phosphopeptides were analyzed on an EASY-nLC 1,200 system coupled to a Q Exactive HF mass spectrometer under the same parameters as in the proteomic analysis. Phosphorylation site abundances were normalized by dividing the abundance of each site by its corresponding protein abundance (Fan et al., 2020). Phosphorylation sites exhibiting significant alterations (fold change >1.2 and p < 0.05) were subjected to functional annotation through GO and KEGG pathway enrichment analyses.
2.7 Molecular docking
Molecular docking simulations were conducted using MOE software (version 2019.0102). The molecular structures of formononetin and four reference drugs (5-Fluorouracil, Capecitabine, Oxaliplatin, and Irinotecan) were retrieved from the PubChem database. The structures of target proteins CA9, MME, and YAP1 were obtained from the Protein Data Bank (https://www.rcsb.org/). Protein structures were processed in MOE, and semi-flexible docking simulations were performed to evaluate interactions between formononetin and the selected targets.
3 Results
3.1 Metabolite profiling of HR
A targeted UPLC-ESI-QTRAP-MS/MS metabolomic analysis focusing on secondary metabolites was conducted on eight heat-processed HR products from different manufacturers. We identified 1,292 metabolites (Supplementary Table S3), classified into eight major categories: flavonoids, phenolic acids, terpenoids, alkaloids, lignans and coumarins, and others. Flavonoids constituted the largest category (29.41%), followed by phenolic acids (18.89%), others (12.54%), alkaloids (12.31%), terpenoids (12.31%), and lignans and coumarins (10.91%) (Figure 1A). There were 26 metabolites that accounted for more than 1% each of the total relative content of the HR metabolome, such as 3-hydroxy-1-methylpyrrolidin-2-one, hypaphorine, formononetin, and ferulic acid. Notably, formononetin and its five derivatives including formononetin-7-O-glucoside and formononetin-7-O-(6″-malonyl)glucoside constituted approximately 4.22% of the total HR metabolome, while the most abundant single metabolite, 3-hydroxy-1-methylpyrrolidin-2-one, was 3.06% of the HR metabolome.
Figure 1. (A) Classification of the 1,292 detected metabolites in HR according to UPLC-ESI-QTRAP-MS/MS analysis. (B) PCA score plot generated from the metabolomic data across the detected samples. (C) Z-score plot of the 1,292 detected metabolites in the HR samples. Metabolites labeled in red or green represent those that are increased or decreased, respectively, in SJD HR.
We further searched the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) to screen for metabolites with potential drug activities (Li et al., 2023). The selection criteria were oral bioavailability (OB) ≥30%, drug-likeness (DL) ≥0.18, molecular weight (MW) ≤500, logP ≤5, hydrogen bond donors (nOHNH) ≤5, and hydrogen bond acceptors (nOH) ≤10. The criteria were consistent with standard TCMSP thresholds. OB ≥ 30% ensures sufficient systemic exposure after oral administration, DL cutoff of 0.18 corresponds to the mean value of known drug molecules in DrugBank. Using these criteria, we identified 11 metabolites with potential drug activity (Table 1), and six of them—medicarpin, calycosin, formononetin, naringenin, quercetin, and vestitol—were detected in our metabolomic profiling of HR.
Table 1. Predicted metabolites with potential activity in HR satisfying OB ≥ 30%, DL ≥ 0.18, MW ≤ 500, miLogP ≤ 5, nOHNH ≤ 5, and nOH ≤ 10.
3.2 Metabolic changes during heat processing of HR
To examine metabolic alterations induced by heat-processing (called paozhi in Chinese) of HR, we performed UPLC-ESI-QTRAP-MS/MS analysis on SJD, the unprocessed form of JD HR. The PCA plot demonstrated clear separation between the three replicates of SJD HR and the other heat-processed HR samples including heat processed JD HR (Figure 1B). Processed HR samples clustered more tightly than SJD HR, indicating increased metabolic consistency post-processing. Z-score analysis identified 11 metabolites in SJD HR as outliers relative to processed HR products, including rhamnetin-3-O-rutinoside-5-O-rhamnoside, chrysoeriol-5,7-di-O-glucoside, physcion-8-O-(6-acetyl)-glucoside, and trans-5-O-(p-coumaroyl)shikimate (Figure 1C). Six metabolites had increased levels post-processing, and three of them were flavonoids, including Rhamnetin-3-O-rutinoside-5-O-rhamnoside and chrysoeriol-5,7-di-O-glucoside whose absolute Z-scores were bigger than 5.
The metabolic differences between the unprocessed and processed JD HR were also investigated using supervised OPLS-DA, and 108 differentially expressed metabolites (VIP >1, p ≤ 0.05) were identified to contribute to the separation of the two groups in OPLS-DA plot (Figure 2). Among these metabolites, 29 were significantly upregulated post-processing, approximately 34% of which were flavonoids, including acetyl wistin, bonannione A, and hesperidin (Supplementary Table S4). Although the total formononetin content in JD HR did not change significantly after heat processing, two formononetin derivatives, 8-methoxy acetyl ononin and formononetin acetyl glucoside, were significantly upregulated. These results indicated that heat processing significantly changed the HR metabolome and particularly increased the contents of flavonoid derivatives.
Figure 2. (A) Score plot from the OPLS-DA model generated from the metabolomic data of JD HR and its original material, SJD HR. (B) S-plot from the OPLS-DA model generated from the metabolomic data of JD HR and its original material, SJD HR. Metabolites with a VIP value greater than one in the OPLS-DA model are labeled in red, while metabolites with a VIP value less than one are labeled in green. (C) Metabolites-metabolites correlation/significance in HR. X and Y-axes were categorized into metabolites. Both r and p values of the correlations were shown in distinct colors. (D) The regulatory metabolic network based on formononetin-related significant correlations (r2 ≥ 0.81; p ≤ 0.05). Metabolites were represented as nodes, displayed different colors in different categories. The relations among metbolites were represented as edges. The positive correlations were displayed in red while the negative ones were displayed in blue. 6-hydroxy-2-(3-methoxybenzylidene)-1-benzofuran-3(2H)-one*: 6H1B-one; 2-(3,4-Dimethoxybenzylidene)-6-Hydroxy-1-Benzofuran-3(2H)-One*: 2–6H1B-one; 2-acetyl-7-methoxy-3H-benzo [f]chromen-3-one: 2A7M-one; 3,4-Dihydro-4-(4-hydroxy-3-methoxyphenyl)-3-(hydroxymethyl)-6,7-dimethoxy-(3R,4S)-2-naphthalenecarboxaldehyde: 3,4-D-N; 4-Aminobenzoic acid: 4-A acid; 3-Aminosalicylic acid: 3-A acid; 1,3-dihydroxy-6-methoxy-7-methylanthraquinone: 1,3-D-M; 2-Methyl-1,3,6-trihydroxy-9,10-anthraquinone: 2-M-anthraquinone; 2-hydroxy-3-hydroxymethyl-anthraquinone: 2-H-anthraquinone.
3.3 Metabolite-metabolite correlation analysis
To reveal the metabolite regulatory network in HR, we performed correlation analysis among the 1,292 detected HR metabolites. The resulted heatmap included 833,986 correlations ranging from −0.9939 to 1 (Figure 2C). When setting the threshold of r2 ≥ 0.49 (r ≥ 0.7 or r ≤ −0.7) and p ≤ 0.05, 88,955 significant correlations remained including 67,457 positive and 21,498 negative correlations. Flavonoids dominated these significant metabolite-metabolite correlations, with 41,174 significant correlations among 380 flavonoids (32,739 positive and 8,435 negative). Next to flavonoids, there were 30,908 significant correlations among 244 phenolic acids, 20,583 significant correlations among 141 lignans and coumarins, 18,625 significant correlations among 159 terpenoids, and 18,520 significant correlations among 159 alkaloids. When screening highly significant correlations with formononetin by r2 ≥ 0.81, 26 correlations were identified (Figure 2D). Four of these 26 correlations were negative, including three alkaloids (hypaphorine, valerine, salicylamide) and one terpenoid (6′-O-Glucosylaucubin). Formononetin had high positive correlation with ten flavonoids, including prunetin, afrormosin, biochanin A, and genkwanin. The other 12 high correlations with formononetin were five lignans and coumarins, four phenolic acids, and three quinones.
3.4 Identification of formononetin and CRC-associated targets via network analysis
Since it has been established via Zinc database analysis (https://zinc.docking.org/) that formononetin is not a pan-assay interference compound (Bolz et al., 2021), we performed network analysis to search formononetin’s putative molecular targets against CRC. We retrieved 220 predicted targets from SwissTargetPrediction, TargetNet, and PharmMapper databases. A total of 14,050 CRC-related genes were obtained from GeneCards, OMIM, and DisGeNET. The Venn diagram revealed 194 overlapping genes between predicted formononetin targets and CRC-related genes (Figures 3A,B). GO enrichment analysis of the 194 genes identified 517 enriched terms: 358 biological processes (BP), 43 cellular components (CC), and 116 molecular functions (MF). The top 10 enriched terms in each category were shown in Supplementary Figure S1, highlighting key processes such as peptidyl-tyrosine phosphorylation, cytosol localization, receptor complex formation, cytoplasmic signaling, steroid binding, zinc ion binding, and protein kinase activity. KEGG pathway enrichment analysis of the 194 genes identified 114 significantly enriched pathways. The top 30 pathways (Figure 3C) included several cancer-related signaling pathways such as pathways in cancer, MAPK signaling, Ras signaling, and PI3K-Akt signaling. These results suggested that formononetin may exert anti-CRC effects by modulating multiple oncogenic signaling pathways.
Figure 3. (A) Screening of overlapping targets between formononetin and CRC. (B) Details of the 194 overlapping targets in the network analysis. (C) KEGG pathway enrichment analysis of the 194 targets.
3.5 Proteomic changes in colon cancer cells treated with formononetin
We evaluated the inhibitory effect of formononetin on SW620 colon cancer cells by measuring cell viability using the CCK-8 assay. Compared to the untreated control cells, formononetin treatment significantly reduced SW620 cell viability at tested concentrations (5, 10, 20, 40, and 80 μM), although the effect was not strictly dose-dependent (Figure 4A). The concentration of 10 μM was selected for further proteomics study since it gave pronounced inhibitory effect without inflicting overwhelming cytotoxicity. DIA proteomics identified 7,381 proteins, of which 291 were significantly altered upon formononetin treatment, and the fold changes ranged from 0.0019 to 8.2382 (Figure 4B; Supplementary Table S5). Specifically, 105 proteins were significantly upregulated, while 186 proteins were significantly downregulated including CCDC18, CA9, IGF1R, and MME.
Figure 4. (A) Viability of SW620 cells significantly inhibited by different concentrations of formononetin solution at 48 h **p < 0.01. (B) Volcano plot generated from the proteomic data of formononetin-treated SW620 cells and controls. (C) Pathway enrichment analysis of 291 differentially expressed proteins in formononetin-treated SW620 cells.
GO enrichment analysis of these 291 differentially expressed proteins identified 3,457 enriched BP terms, 372 CC terms, and 404 MF terms. The top ten enriched GO terms in each category were displayed in Supplementary Figure S2, including single-organism cellular process, multicellular organismal process, cellular anatomical entity, cytoplasmic part, protein binding, and carboxylic acid binding. KEGG pathway enrichment analysis revealed 17 significantly enriched pathways, including metabolic pathways, Hippo signaling, HIF-1 signaling, nitrogen metabolism, one-carbon pool by folate, and cholesterol metabolism (Figure 4C). KEGG pathway enrichment analyses of upregulated and downregulated proteins were also conducted and shown in Supplementary Tables S6, S7. These findings indicated that formononetin caused broad proteomic changes involving multiple metabolic and oncogenic pathways in SW620 cells.
3.6 Phosphoproteomic alterations in colon cancer cells treated with formononetin
Further investigation into formononetin’s molecular mechanisms of cancer cell inhibition was conducted using DIA-based phosphoproteomic analysis. We quantified 19,038 phosphorylation sites, among which 2,587 sites in 1,535 phosphorylated proteins were significantly altered (Figure 5A; Supplementary Table S8). These phosphorylation sites comprised 79.1% phosphoserine, 16.4% phosphothreonine, and 4.5% phosphotyrosine modifications. Among them, 1,244 phosphorylation sites were significantly upregulated (involving proteins such as VAPA, CA9, YY1, FOXK2, SQSTM1), while 1,343 sites were downregulated (in proteins such as AKT2, YWHAB, MME, YAP1). GO enrichment analysis of the phosphoproteins identified 7,832 significantly enriched terms (p < 0.05): 6,104 BP, 898 MF, and 830 CC. The most enriched terms were cellular component organization or biogenesis, protein binding, poly(A) RNA binding, intracellular organelle part, and nuclear compartmentalization (Supplementary Figure S3).
Figure 5. (A) Volcano plot generated from the phosphoproteomic data of formononetin-treated SW620 cells and controls. (B) Pathway enrichment analysis of 1,535 significantly altered phosphorylated proteins in formononetin-treated SW620 cells. (C) Comparison of significantly enriched pathways common among network, proteomics, and phosphoproteomics analyses. (D) Comparison of common proteins (targets) generated from network, proteomics, and phosphoproteomics analyses.
KEGG pathway enrichment analysis highlighted 41 significantly enriched pathways (p < 0.05), such as cell cycle, autophagy, Hippo signaling, mTOR signaling, AMPK signaling, glycolysis/gluconeogenesis, lysine degradation, and insulin signaling (Figure 5B). While there was no commonly enriched pathways across network analysis, proteomics, and phosphoproteomics, 19 pathways were enriched by at least two analyses, including Hippo signaling and lysine degradation by both proteomics and phosphoproteomics, and proteoglycans in cancer and glycolysis/gluconeogenesis by both network analysis and phosphoproteomics (Figure 5C). These findings highlighted formononetin’s multifaceted regulatory effects on key oncogenic and metabolic pathways.
3.7 Target verification via molecular docking
To further validate key molecular targets of formononetin, we compared proteins identified by network analysis, proteomics, and phosphoproteomics (Figure 5D). Two targets, CA9 and MME, were consistently identified across all three analyses, suggesting their central roles in formononetin-mediated anti-cancer effects. Additionally, 50 proteins were common between proteomic and phosphoproteomic datasets (Supplementary Table S9), while 11 targets, including FYN and DNMT1, were shared between network analysis and phosphoproteomics. We confirmed the binding affinity of formononetin to CA9, MME, and YAP1 (a key regulator of Hippo signaling) by performing molecular docking. Control drugs 5-Fluorouracil, Capecitabine, Oxaliplatin, and Irinotecan were included for comparison. The binding energy scores for formononetin-CA9, formononetin-MME, and formononetin-YAP1 in molecular docking analysis were −4.8162, −6.1931, and −4.3183 kcal/mol, respectively. These scores were higher than those for Capecitabine and Irinotecan, but lower than those for 5-Fluorouracil and Oxaliplatin (Supplementary Table S10). The complete docking analysis results demonstrated that formononetin had strong binding affinity to CA9, MME, and YAP1, indicating they were potential formononetin targets (Supplementary Figure S4).
4 Discussion
Hedysari Radix (HR) is well-known as a TCM for its diverse pharmacological properties, including antioxidant, anti-tumor, anti-inflammatory, lipid-lowering, and antihypertensive activities. Similar to another traditional Chinese medicine Astragalus membranaceus (Fisch.) Bunge [Fabaceae], HR has been traditionally used for its tonic effects on “Qi” (Bai et al., 2020; Tsai et al., 2022). Comprehensive understandings of specialized secondary metabolites in TCMs are crucial in understanding TCM’s pharmacological efficacy, and the secondary metabolite profile in HR has not been completely investigated. In this study, we applied targeted UPLC-ESI-QTRAP-MS/MS analysis to profile secondary metabolites in eight heat-processed HR products. A total of 1,292 secondary metabolites including flavonoids, phenolic acids, alkaloids, and terpenoids were identified. Compared to previous studies (Bai et al., 2020; Tsai et al., 2022), our findings expanded the inventory of secondary metabolome of HR, demonstrating the power of metabolomic technologies in elucidating TCM’s molecular composition.
Among the identified metabolites, 380 belonged to the flavonoid group, which were well-documented for their anti-cancer properties in clinical and preclinical studies (Khan et al., 2021). Several flavonoids with known anti-tumor activities, including formononetin, prunetin, kaempferide, medicarpin, and afrormosin, were detected in HR (Wu et al., 2017; Li et al., 2020; Khan et al., 2021; Köksal Karayildirim et al., 2021; Kim et al., 2022; Chen et al., 2023), and each of them were more than 0.5% of the total metabolite content. Besides flavonoids, abundant phenolic acid metabolites such as ferulic acid, isoferulic acid, and acetyl resveratrol were also identified and had reported anti-cancer effects (Tino et al., 2016; Gao et al., 2018; Long et al., 2020). Further validation using TCMSP identified 11 predicted metabolites with potential activity in HR, and six of them including formononetin, naringenin, medicarpin, quercetin, were consistently detected in different HR products.
Heat-processing (pronounced Paozhi in Chinese) is an ancient traditional Chinese medicine processing technique to enhance therapeutic efficacy and safety of TCMs. We investigated metabolic variations between unprocessed HR (SJD HR) and other heat-processed HR products including JD HR. The PCA plot of metabolomic data revealed that heat-processing improved HR homogeneity in different HR materials. Moreover, the levels of several pharmacologically important metabolites, such as coumarin, vicenin-3, and rhamnetin-3-O-rutinoside-5-O-rhamnoside, were significantly elevated after heat-processing (Liu et al., 2023; Taherkhani et al., 2023; Todorov et al., 2023). Approximately 8% of the detected metabolites were significantly altered (p ≤ 0.05) between SJD HR and JD HR. Notably, 34% of the significantly increased metabolites post-processing were flavonoids, including two formononetin derivatives. So heat processing seemed to enhance HR’s pharmacological potency by increasing levels of flavonoids such as formononetin derivatives. These findings support the long-standing hypothesis that heat-processing could enhance TCM’s efficacy from a molecular basis.
Metabolite-metabolite correlation analysis has been widely used to elucidate key regulatory pathways (Toubiana et al., 2012; Ji et al., 2017), so we also conducted a comprehensive metabolite-metabolite correlation analysis of HR secondary metabolites. A large number of significant correlations were revealed, and over 75% of them were positive correlations. Flavonoids dominated the significant metabolite-metabolite correlations, followed by phenolic acids, lignans and coumarins, terpenoids, and alkaloids. Most of flavonoid-related correlations were positive, while the majority of negative correlations were associated with terpenoids or alkaloids. One of the significant correlations involving flavonoids were between formononetin and biochanin A. These two metabolites have been identified as the most abundant flavonoids in M. truncatula roots, and they share structural similarities and both have hepatoprotective and anti-ulcer activities (Sharma and Kabra, 2025). Formononetin is the precursor of afrormosin, and significant correlation between them were also identified (Al-Ani and Dewick, 1980).
Formononetin is notable for its antitumorigenic properties, which have been demonstrated by both in vitro and in vivo studies. The antitumorigenic properties were reported to function through regulation of multiple oncogenic signaling pathways, such as PI3K/AKT, MAPK, and apoptotic regulators Bax, Bcl-2, and caspase-3 (Tay et al., 2019; Aliya et al., 2023). Recent studies have shown that formononetin could ameliorate DSS-induced colitis by inhibiting MAPK/PPAR-γ/NF-κB/ROS signaling and suppress colitis-associated colon carcinogenesis through regulating lipid metabolism and inhibiting the mTORC2/Akt axis (Cao et al., 2025; Liu et al., 2025). The molecular mechanisms underlying formononetin’s anti-CRC activity, however, remain largely unexplored. In the current study, we conducted an in-depth mechanistic investigation of formononetin’s effects on colon cancer cells by integrating network analysis, proteomics, and phosphoproteomics. Our findings revealed that a broad range of proteins were affected by formononetin treatment to inhibit SW620 colorectal cancer cell’s proliferation. Specifically, 50 proteins and their corresponding phosphorylation sites were identified as key participants in formononetin-mediated anti-tumor activity. Among these proteins, MME, CDH1, APC, and CA9 have not been previously reported to be involved in anti-cancer regulations. MME and CA9 are particularly important since their regulation by formononetin were consistently observed across network analysis, proteomics, and phosphoproteomics.
CA9 is typically expressed at low levels in normal tissues (primarily gastrointestinal tract) but highly upregulated in multiple cancers including CRC (Ronca and Supuran, 2024). Its overexpression correlates with poor prognosis and chemoradiotherapy resistance in CRC patients, making it a potentially important therapeutic target (Kivela et al., 2001; Guedj et al., 2011). The proteomic data revealed that the level of CA9 protein in SW620 cells was downregulated while its phosphorylation at Ser448 was significantly increased following formononetin treatment. Given that dephosphorylation at Ser448 is required for full enzymatic activity, these findings suggested that formononetin may exert anti-cancer effects by disrupting CA9 enzyme function (Ditte et al., 2011). Similarly, MME has also been implicated in tumorigenesis of various cancers including CRC (Haase et al., 2017). Both MME protein levels and its phosphorylation were significantly altered following formononetin treatment in SW620 cells. Molecular docking analysis confirmed that formononetin could strongly interact with CA9 and MME, indicating their critical roles in formononetin’s suppression of colon cancer cells.
Several studies have demonstrated that formononetin could inhibit colon carcinoma cell growth in vitro and in vivo (Huang et al., 2015; Wang et al., 2018). Our previous study showed that the Yiqi Sanjie (YQSJ) formula, which contained formononetin, was effective in treating stage III CRC patients and AOM/DSS-induced CRC model mice (Wan et al., 2022; Zheng et al., 2022). Several cancer-related pathways regulated by YQSJ were also significantly enriched in the current study, such as cholesterol metabolism and Hippo signaling. The Hippo signaling pathway was significantly enriched in both proteomics and phosphoproteomics analyses. This pathway is highly conserved and pivotal in tumorigenesis, organ size control, and cell proliferation, making it an attractive cancer therapy target (Jin et al., 2023). While the protein levels of YAP1, the central effector of Hippo signaling pathway, did not change following formononetin treatment, its Thr110 phosphorylation was significantly reduced. Since YAP1 phosphorylation is critical for its nuclear localization and transcription, the decreased Thr110 phosphorylation suggested that formononetin might inhibit oncogenic signaling through affecting YAP1’s phosphorylation (Piccolo et al., 2014; Lin et al., 2023). Besides Hippo signaling, 14 other pathways were also significantly enriched in both network and phosphoproteomics analyses, including cancer-associated pathways such as insulin signaling and cell cycle regulation (Aliya et al., 2023). Cholesterol metabolism and HIF-1 signaling pathways in CRC cancer cells were also affected by formononetin treatment, and HIF-1 signaling is closely linked to CA9 function and hypoxic adaptation in tumors. These results highlighted formononetin’s effects in disrupting multiple oncogenic pathways and its potential as a valid cancer therapeutic agent.
To our knowledge, this is the first comprehensive study integrating proteomics and phosphoproteomics to systematically elucidate HR and formononetin’s anti-CRC mechanisms, and potential formononetin targets CA9 and MME were identified. There were several limitations of this study. First, network pharmacology and molecular docking analyses were computational and required further experimental confirmation. These evidences of formononetin’s target engagement were preliminary and not direct. Future investigations should be conducted to validate these predictions according to established pharmacological standards such as the “four pillars of best practice” (Heinrich et al., 2022). Secondly, mechanistic analyses were limited to a single CRC cell line (SW620), potentially restricting generalizability of the results. The response of SW620 cells to formononetin was not strictly dose-dependent, which may stem from cell line-specific sensitivity or biphasic (hormetic) effect, highlighting the need for additional experiments using more cell lines (e.g., SW1116, HCT116) across extended time points (e.g., 24 and 72 h) and wider treatment concentration range (e.g., 1–5 μmol/L) (Jodynis-Liebert and Kujawska, 2020). Finally, the physiological relevance of identified formononetin targets and related pathways need to be validated in vivo such as in animal models, which is essential to substantiate formononetin’s therapeutic potential in CRC.
5 Conclusion
In summary, this study provided a comprehensive metabolomic characterization of HR, revealed metabolic changes in HR after heat processing. The anti-tumor mechanisms of formononetin in HR were explored, and its potential regulation of multiple targets and signaling pathways were revealed. In particular, regulation of HIF-1-CA9 and Hippo-YAP1 signaling pathways were highlighted as the key mechanisms of formononetin’s anti-CRC activity. These findings offer new insights into the secondary metabolites in HR and the anti-tumor properties of formononetin, provide a foundation for developing novel therapeutic strategies for CRC.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
JR: Writing – review and editing, Writing – original draft, Funding acquisition, Conceptualization. XW: Methodology, Investigation, Validation, Writing – review and editing. TC: Writing – review and editing, Formal Analysis, Investigation. MM: Writing – review and editing, Methodology, Investigation, Validation. CX: Conceptualization, Writing – review and editing, Resources. S-YC: Formal Analysis, Writing – review and editing, Validation, Investigation. S-WH: Methodology, Writing – review and editing, Formal Analysis. XH: Conceptualization, Writing – review and editing, Supervision. C-HC: Supervision, Conceptualization, Writing – review and editing, Writing – original draft. ZZ: Conceptualization, Methodology, Writing – review and editing, Funding acquisition, Resources.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This work was supported by the “Five-Level Progressive” Talent Cultivation Project of Jiangxi Cancer Hospital & Institute (No. WCDJ2024JQ02), the Jiangxi Provincial Natural Science Foundation (Project Nos. 20224BAB206114, 20232BAB216074, and 20242BAB20254), the Jiangxi Province Talent Leader Plan of Major Disciplines (20232BCJ22057), grants from the Department of Science and Technology of Jiangxi Province (No. 20212BAG70040), and the Jiangxi Provincial Health and Health Commission Science and Technology Plan Project (202510022).
Conflict of interest
Author CX was employed by Jiangxi Jingde Traditional Chinese Medicine Co., Ltd.
The remaining author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2025.1649691/full#supplementary-material
References
Al-Ani, H. A. M., and Dewick, P. M. (1980). Isoflavone biosynthesis in onobrychis viciifolia: Formononetin and texasin as precursors of afrormosin. Phytochemistry 19, 2337–2339. doi:10.1016/S0031-9422(00)91022-7
Aliya, S., Alhammadi, M., Park, U., Tiwari, J. N., Lee, J. H., Han, Y. K., et al. (2023). The potential role of formononetin in cancer treatment: an updated review. Biomed. Pharmacother. 168, 115811. doi:10.1016/j.biopha.2023.115811
Bai, H. Y., Zheng, W. H., Han, S., Bao, F., Sun, L. L., Zhang, K. X., et al. (2020). Metabolomic determination of specialized metabolites using liquid chromatography-tandem mass spectrometry in the traditional Chinese medicines astragali radix and Hedysari Radix. Nat. Prod. Commun. 15, 1934578X19901192. doi:10.1177/1934578X19901192
Bolz, S. N., Adasme, M. F., and Schroeder, M. (2021). Toward an understanding of pan-assay interference compounds and promiscuity: a structural perspective on binding modes. J. Chem. Inf. Model. 61, 2248–2262. doi:10.1021/acs.jcim.0c01227
Cao, S., Lv, B., Tai, Y., Zuo, H. X., Xing, Y., Surh, Y. J., et al. (2025). Formononetin ameliorates DSS-Induced colitis by inhibiting the MAPK/PPAR-γ/NF-κB/ROS signaling pathways. Toxicol. Appl. Pharmacol. 496, 117239. doi:10.1016/j.taap.2025.117239
Chen, T., Ma, J., Liu, Y., Chen, Z., Xiao, N., Lu, Y., et al. (2022). iProX in 2021: connecting proteomics data sharing with big data. Nucleic Acids Res. 50, D1522–D1527. doi:10.1093/nar/gkab1081
Chen, Y., Yin, L., Hao, M., Xu, W., Gao, J., Sun, Y., et al. (2023). Medicarpin induces G1 arrest and mitochondria-mediated intrinsic apoptotic pathway in bladder cancer cells. Acta Pharm. 73, 211–225. doi:10.2478/acph-2023-0016
Ditte, P., Dequiedt, F., Svastova, E., Hulikova, A., Ohradanova-Repic, A., Zatovicova, M., et al. (2011). Phosphorylation of carbonic anhydrase IX controls its ability to mediate extracellular acidification in hypoxic tumors. Cancer Res. 71, 7558–7567. doi:10.1158/0008-5472.CAN-11-2520
El Bali, M., Bakkach, J., and Bennani Mechita, M. (2021). Colorectal cancer: from genetic landscape to targeted therapy. J. Oncol. 2021, 9918116. doi:10.1155/2021/9918116
Fan, Y., Cheng, Y., Li, Y., Chen, B., Wang, Z., Wei, T., et al. (2020). Phosphoproteomic analysis of neonatal regenerative myocardium revealed important roles of checkpoint kinase 1 via activating mammalian target of rapamycin C1/ribosomal protein S6 kinase b-1 pathway. Circulation 141, 1554–1569. doi:10.1161/CIRCULATIONAHA.119.040747
Feng, W., Liu, J., Huang, L., Tan, Y., and Peng, C. (2021). Gut microbiota as a target to limit toxic effects of traditional Chinese medicine: implications for therapy. Biomed. Pharmacother. 133, 111047. doi:10.1016/j.biopha.2020.111047
Gao, J., Yu, H., Guo, W., Kong, Y., Gu, L., Li, Q., et al. (2018). The anticancer effects of ferulic acid is associated with induction of cell cycle arrest and autophagy in cervical cancer cells. Cancer Cell Int. 18, 102. doi:10.1186/s12935-018-0595-y
Groelly, F. J., Fawkes, M., Dagg, R. A., Blackford, A. N., and Tarsounas, M. (2023). Targeting DNA damage response pathways in cancer. Nat. Rev. Cancer 23, 78–94. doi:10.1038/s41568-022-00535-5
Guedj, N., Bretagnol, F., Rautou, P. E., Deschamps, L., Cazals-Hatem, D., Bedossa, P., et al. (2011). Predictors of tumor response after preoperative chemoradiotherapy for rectal adenocarcinomas. Hum. Pathol. 42, 1702–1709. doi:10.1016/j.humpath.2011.01.015
Haase, G., Gavert, N., Brabletz, T., and Ben-Ze'ev, A. (2017). A point mutation in the extracellular domain of L1 blocks its capacity to confer metastasis in colon cancer cells via CD10. Oncogene 36, 1597–1606. doi:10.1038/onc.2016.329
Heinrich, M., Jalil, B., Abdel-Tawab, M., Echeverria, J., Kulić, Ž., McGaw, L. J., et al. (2022). Best practice in the chemical characterisation of extracts used in pharmacological and toxicological research-the ConPhyMP-guidelines. Front. Pharmacol. 13, 953205. doi:10.3389/fphar.2022.953205
Huang, J., Xie, M., Gao, P., Ye, Y., Liu, Y., Zhao, Y., et al. (2015). Antiproliferative effects of formononetin on human colorectal cancer via suppressing cell growth in vitro and in vivo. Process Biochem. 50, 912–917. doi:10.1016/j.procbio.2015.03.001
Huang, N., Huang, W., Wu, J., Long, S., Luo, Y., and Huang, J. (2024). Possible opportunities and challenges for traditional Chinese medicine research in 2035. Front. Pharmacol. 15, 1426300. doi:10.3389/fphar.2024.1426300
Ji, Y., Rao, J., Rong, X., Lou, S., Zheng, Z., and Lu, Y. (2017). Metabolic characterization of human aqueous humor in relation to high myopia. Exp. Eye Res. 159, 147–155. doi:10.1016/j.exer.2017.03.004
Jin, F., Ni, X., Yu, S., Jiang, X., Shi, X., Zhou, J., et al. (2023). The Ethyl acetate extract from Celastrus orbiculatus suppresses non-small-cell lung cancer by activating Hippo signaling and inhibiting YAP nuclear translocation. Phytomedicine 114, 154761. doi:10.1016/j.phymed.2023.154761
Jodynis-Liebert, J., and Kujawska, M. (2020). Biphasic dose-response induced by phytochemicals: experimental evidence. J. Clin. Med. 9, 718. doi:10.3390/jcm9030718
Khan, H., Belwal, T., Efferth, T., Farooqi, A. A., Sanches-Silva, A., Vacca, R. A., et al. (2021). Targeting epigenetics in cancer: therapeutic potential of flavonoids. Crit. Rev. Food Sci. Nutr. 61, 1616–1639. doi:10.1080/10408398.2020.1763910
Kim, H., Han, M., Shin, S. A., An, J., Ahn, M. J., Lee, J. H., et al. (2022). Afrormosin exerts an anticancer effect via MAPK and AKT signaling pathways in B16F10 cells. Appl. Biol. Chem. 65, 71. doi:10.1186/s13765-022-00743-5
Kivela, A. J., Saarnio, J., Karttunen, T. J., Kivelä, J., Parkkila, A. K., Pastorekova, S., et al. (2001). Differential expression of cytoplasmic carbonic anhydrases, CA I and II, and membrane-associated isozymes, CA IX and XII, in normal mucosa of large intestine and in colorectal tumors. Dig. Dis. Sci. 46, 2179–2186. doi:10.1023/a:1011910931210
Köksal Karayildirim, Ç., Nalbantsoy, A., and Karabay Yavaşoğlu, N. Ü. (2021). Prunetin inhibits nitric oxide activity and induces apoptosis in urinary bladder cancer cells via CASP 3 and TNF-α genes. Mol. Biol. Rep. 48, 7251–7259. doi:10.1007/s11033-021-06719-w
Li, H. L., Li, S. M., Luo, Y. H., Xu, W. T., Zhang, Y., Zhang, T., et al. (2020). Kaempferide induces G0/G1 phase arrest and apoptosis via ROS-mediated signaling pathways in A549 human lung cancer cells. Nat. Prod. Commun. 15 (7), 1934578X20935226. doi:10.1016/j.lwt.203.114451
Li, T., Li, W., Guo, X., Tan, T., Xiang, C., and Ouyang, Z. (2023). Unraveling the potential mechanisms of the anti-osteoporotic effects of the Achyranthes bidentata-Dipsacus asper herb pair: a network pharmacology and experimental study. Front. Pharmacol. 14, 1242194. doi:10.3389/fphar.2023.1242194
Lin, Q., Cao, J., Yu, J., Zhu, Y., Shen, Y., Wang, S., et al. (2023). YAP-mediated trophoblast dysfunction: the common pathway underlying pregnancy complications. Cell Commun. Signal 21, 353. doi:10.1186/s12964-023-01371-2
Liu, Y., Li, Y., Liu, Z., Wang, L., Bi, Z., Sun, C., et al. (2023). Integrated transcriptomic and metabolomic analysis revealed altitude-related regulatory mechanisms on flavonoid accumulation in potato tubers. Food Res. Int. 170, 112997. doi:10.1016/j.foodres.2023.112997
Liu, M., Yang, C., Peng, X., Zheng, S., He, H., Wang, W., et al. (2025). Formononetin suppresses colitis-associated colon cancer by targeting lipid synthesis and mTORC2/Akt signaling. Phytomedicine 142, 156665. doi:10.1016/j.phymed.2025.156665
Long, Z., Feng, G., Zhao, N., Wu, L., and Zhu, H. (2020). Isoferulic acid inhibits human leukemia cell growth through induction of G2/M-phase arrest and inhibition of Akt/mTOR signaling. Mol. Med. Rep. 21, 1035–1042. doi:10.3892/mmr.2020.10926
Miyamoto, Y., Ogawa, K., Ohuchi, M., Tokunaga, R., and Baba, H. (2023). Emerging evidence of immunotherapy for colorectal cancer. Ann. Gastroenterol. Surg. 7, 216–224. doi:10.1002/ags3.12633
Mo, X., Guo, D., Jiang, Y., Chen, P., and Huang, L. (2022). Isolation, structures and bioactivities of the polysaccharides from Radix Hedysari: a review. Int. J. Biol. Macromol. 199, 212–222. doi:10.1016/j.ijbiomac.2021.12.095
Morgan, E., Arnold, M., Gini, A., Lorenzoni, V., Cabasag, C. J., Laversanne, M., et al. (2023). Global burden of colorectal cancer in 2020 and 2040: incidence and mortality estimates from GLOBOCAN. Gut 72, 338–344. doi:10.1136/gutjnl-2022-327736
Piccolo, S., Dupont, S., and Cordenonsi, M. (2014). The biology of YAP/TAZ: hippo signaling and beyond. Physiol. Rev. 94, 1287–1312. doi:10.1152/physrev.00005.2014
Qu, J., Ke, F., Liu, Z., Yang, X., Li, X., Xu, H., et al. (2022). Uncovering the mechanisms of dandelion against triple-negative breast cancer using a combined network pharmacology, molecular pharmacology and metabolomics approach. Phytomedicine 99, 153986. doi:10.1016/j.phymed.2022.153986
Ronca, R., and Supuran, C. T. (2024). Carbonic anhydrase IX: an atypical target for innovative therapies in cancer. Biochim. Biophys. Acta Rev. Cancer 1879, 189120. doi:10.1016/j.bbcan.2024.189120
Sharma, N., and Kabra, A. (2025). Formononetin: pharmacological properties and therapeutic potential. Naunyn Schmiedeb. Arch. Pharmacol. 398, 15053–15073. doi:10.1007/s00210-025-04247-z
Song, X., He, H., Zhang, Y., Fan, J., and Wang, L. (2022). Mechanisms of action of triptolide against colorectal cancer: insights from proteomic and phosphoproteomic analyses. Aging (Albany NY) 14, 3084–3104. doi:10.18632/aging.203992
Sun, Y., Zhou, M., Luo, L., Pan, H., Zhang, Q., and Yu, C. (2023). Metabolic profiles, bioactive compounds and antioxidant activity of rosehips from Xinjiang, China. LWT 174, 114451. doi:10.1016/j.lwt.2023.114451
Taherkhani, A., Khodadadi, P., Samie, L., Azadian, Z., and Bayat, Z. (2023). Flavonoids as strong inhibitors of MAPK3: a computational drug discovery approach. Int. J. Anal. Chem. 2023, 8899240. doi:10.1155/2023/8899240
Tay, K. C., Tan, L. T., Chan, C. K., Hong, S. L., Chan, K. G., Yap, W. H., et al. (2019). Formononetin: a review of its anticancer potentials and mechanisms. Front. Pharmacol. 10, 820. doi:10.3389/fphar.2019.00820
Tino, A. B., Chitcholtan, K., Sykes, P. H., and Garrill, A. (2016). Resveratrol and acetyl-resveratrol modulate activity of VEGF and IL-8 in ovarian cancer cell aggregates via attenuation of the NF-κB protein. J. Ovarian Res. 99, 84. doi:10.1186/s13048-016-0293-0
Todorov, L., Saso, L., and Kostova, I. (2023). Antioxidant activity of coumarins and their metal complexes. Pharmaceuticals 16, 651. doi:10.3390/ph16050651
Toubiana, D., Semel, Y., Tohge, T., Beleggia, R., Cattivelli, L., Rosental, L., et al. (2012). Metabolic profiling of a mapping population exposes new insights in the regulation of seed metabolism and seed, fruit, and plant relations. PLoS Genet. 8, e1002612. doi:10.1371/journal.pgen.1002612
Tsai, Y. C., Lin, M. K., Peng, W. H., Tseng, C. K., Lee, M. S., Yang, B. C., et al. (2022). Comparison of the immunomodulatory effect of TCM formulas containing either Astragali Radix or with this replaced by Hedysari Radix. Nat. Prod. Commun. Nat. Prod. Commun. 17, 1934578X221142797. doi:10.1177/1934578X221142797
Wan, X., Tou, F., Zeng, J., Chen, X., Li, S., Chen, L., et al. (2022). Integrative analysis and identification of key elements and pathways regulated by Traditional Chinese Medicine (Yiqi Sanjie formula) in colorectal cancer. Front. Pharmacol. 13, 1090599. doi:10.3389/fphar.2022.1090599
Wang, M., Chen, L., Liu, D., Chen, H., Tang, D. D., and 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
Wang, A. L., Li, Y., Zhao, Q., and Fan, L. Q. (2018). Formononetin inhibits colon carcinoma cell growth and invasion by microRNA-149-mediated EphB3 downregulation and inhibition of PI3K/AKT and STAT3 signaling pathways. Mol. Med. Rep. 17, 7721–7729. doi:10.3892/mmr.2018.8857
Wang, Z., Kim, U., Jiao, Y., Li, C., Guo, Y., Ma, X., et al. (2019). Quantitative proteomics combined with affinity MS revealed the molecular mechanism of ginsenoside antitumor effects. J. Proteome Res. 18, 2100–2108. doi:10.1021/acs.jproteome.8b00972
Wang, G., Wang, Y., Zhang, L., Cai, Q., Lin, Y., Lin, L., et al. (2020). Proteomics analysis reveals the effect of Aeromonas hydrophila sirtuin CobB on biological functions. J. Proteomics 225, 103848. doi:10.1016/j.jprot.2020.103848
Wang, Y., Zhang, X., Li, J., Zhang, Y., Guo, Y., Chang, Q., et al. (2021). Sini decoction ameliorates colorectal cancer and modulates the composition of gut microbiota in mice. Front. Pharmacol. 12, 609992. doi:10.3389/fphar.2021.609992
Wang, J. R., Song, X. H., Li, L. Y., Gao, S. J., Shang, F. H., Zhang, X. M., et al. (2022). Metabolomic analysis reveals dynamic changes in secondary metabolites of Sophora japonica L. during flower maturation. Front. Plant Sci. 13, 916410. doi:10.3389/fpls.2022.916410
Wang, Y., Xie, L., Liu, F., Ding, D., Wei, W., and Han, F. (2024). Research progress on traditional Chinese medicine-induced apoptosis signaling pathways in ovarian cancer cells. J. Ethnopharmacol. 319, 117299. doi:10.1016/j.jep.2023.117299
Wei, Z., Chen, J., Zuo, F., Guo, J., Sun, X., Liu, D., et al. (2023). Traditional Chinese Medicine has great potential as candidate drugs for lung cancer: a review. J. Ethnopharmacol. 300, 115748. doi:10.1016/j.jep.2022.115748
Wu, Y., Zhang, X., Li, Z., Yan, H., Qin, J., and Li, T. (2017). Formononetin inhibits human bladder cancer cell proliferation and invasiveness via regulation of miR-21 and PTEN. Food Funct. 8, 1061–1066. doi:10.1039/c6fo01535b
Xia, Z., Li, M., Hu, M., Lin, Y., Atteh, L. L., Fu, W., et al. (2022). Phosphoproteomics reveals that cinobufotalin promotes intrahepatic cholangiocarcinoma cell apoptosis by activating the ATM/CHK2/p53 signaling pathway. Front. Oncol. 12, 982961. doi:10.3389/fonc.2022.982961
Yu, Q., Xu, C., Song, J., Jin, Y., and Gao, X. (2024). Mechanisms of traditional Chinese medicine/natural medicine in HR-positive Breast Cancer: a comprehensive Literature review. J. Ethnopharmacol. 319, 117322. doi:10.1016/j.jep.2023.117322
Zheng, Z., Wei, Q., Wan, X., Zhong, X., Liu, L., Zeng, J., et al. (2022). Correlation analysis between trace elements and colorectal cancer metabolism by integrated serum proteome and metabolome. Front. Immunol. 13, 921317. doi:10.3389/fimmu.2022.921317
Keywords: anti-tumor, CA9, colorectal cancer, formononetin, hedysari radix, MME
Citation: Rao J, Wang X, Chen T, Mo M, Xu C, Chang S-Y, Hsu S-W, Han X, Chen C-H and Zheng Z (2026) Multi-omics study of the anti-colorectal cancer mechanisms of formononetin in Hedysari Radix. Front. Pharmacol. 16:1649691. doi: 10.3389/fphar.2025.1649691
Received: 18 June 2025; Accepted: 10 December 2025;
Published: 06 January 2026.
Edited by:
Stalin Antony, University of Electronic Science and Technology of China, ChinaReviewed by:
Anjana Dwivedi, Birla Institute of Technology, Mesra, IndiaHuan Liu, Changchun University of Chinese Medicine, China
Copyright © 2026 Rao, Wang, Chen, Mo, Xu, Chang, Hsu, Han, Chen and Zheng. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhi Zheng, emhlbmd4aWFfMjAwN0AxNjMuY29t; Ching-Hsien Chen, amNoY2hlbkB1Y2RhdmlzLmVkdQ==
†These authors have contributed equally to this work