- 1Central Laboratory, The First Affiliated Hospital of Hebei North University, Zhangjiakou, China
- 2Department of General Surgery, The First Affiliated Hospital of Hebei North University, Zhangjiakou, China
- 3Department of Obstetrics and Gynecology, Shenzhen University General Hospital, Shenzhen, China
- 4Guangdong Key Laboratory for Biomedical Measurements and Ultrasound Imaging, School of Biomedical Engineering, Shenzhen University Health Science Center, Shenzhen, China
- 5Shenzhen Key Laboratory, Shenzhen University General Hospital, Shenzhen, China
- 6Department of Internal Medicine of Traditional Chinese Medicine, The First Clinical Medical School, Guangzhou University of Chinese Medicine, Guangzhou, China
Patients with colonic adenocarcinoma (COAD) are at relatively high risk of SARS-CoV-2 infection. However, there is a lack of medical strategies to treat COVID-19/COAD comorbidity. Puerarin, a natural product, is a known antiviral, antitumor, and immunomodulatory effect. Therefore, we hypothesised that puerarin could be used to treat COVID-19/COAD patients. Based on network pharmacology and bioinformatics analysis, the potential targets and pharmacological mechanisms of puerarin in COVID-19/COAD were identified. By intersecting therapeutic target genes for puerarin, COVID-19-related genes and COAD-related genes, 42 target genes of puerarin that could potentially treat COVID-19/COAD comorbidity were obtained. By using the 42 potential target genes to construct the protein-protein interaction (PPI) network, we obtained five core target genes, namely RELA, BCL2, JUN, FOS, and MAPK1. The results of bioinformatics analysis revealed that puerarin could be able to treat COVID-19/COAD comorbidity through apoptosis, antiviral, antioxidant, NF-κB signaling pathway, MAPK signaling pathway, IL-17 signaling pathway, TNF signaling pathway, and HIF-1 signaling pathway etc. This study found that puerarin has the potential to treat COVID-19/COAD patients and that the therapeutic target genes obtained in the study may provide clues for the treatment of COVID19/COAD comorbidity.
Introduction
Since the emergence of the SARS-CoV-2 virus, COVID-19 has spread rapidly worldwide. The latest data from the World Health Organization show that the cumulative number of confirmed cases and deaths worldwide has been more than 500 million and six million (updated on 05 May 2022). With a global daily average of about 500,000 patients with COVID-19, the situation remains critical for epidemic prevention. Similar to SARS-CoV, SARS-CoV-2 was a highly pathogenic coronavirus, which were mainly transmitted through aerosols or droplets, leading to the development of respiratory diseases and even death (Xu et al., 2020). Due to the small size of droplets and aerosols carrying the SARS-CoV-2 virus, it is difficult to achieve the expected preventive effect (Jayaweera et al., 2020). The common clinical symptoms in most SARS-CoV-2 infected patients were cough and fever, with partial patients also presenting with sore throat, diarrhea, and loss of taste (Giacomelli et al., 2020; Li et al., 2020). However, a few of patients with COVID-19 probably developed acute respiratory failure or acute respiratory distress syndrome, which might even lead to death (Chen et al., 2020). The emergence of the SARS-CoV-2 virus has caused enormous economic and health losses to people around the world. Through the continuous efforts of scientists from around the world, several vaccines and drugs have been used for COVID-19 (Liu et al., 2020; Lazarus et al., 2021). However, there was no significant change in the number of additional COVID-19 patients per day.
Both malignant tumor and COVID-19 are major public health problems worldwide. Colorectal cancer, a common gastrointestinal malignancy, had the third-highest incidence and mortality rate of all cancers in the United States, for both men and women (Siegel et al., 2022). With the rapid development of modern medicine, significant breakthroughs have been made in basic research, early diagnosis, chemotherapy, and targeted therapy for colorectal cancer, yet the survival rate of colorectal cancer patients has not been significantly improved (Durinikova et al., 2021). Colon adenocarcinoma (COAD) is a common malignant tumor among all types of colon cancer, with a high mortality rate and high recurrence rate (Liu, 2021). Due to the lack of specific early diagnostic symptoms and markers, most patients were already at an advanced stage at the time of diagnosis, and the high metastatic rate characteristic of advanced COAD was a major factor in the high mortality rate of COAD patients (Brown and Solomon, 2018). Therefore, exploring the mechanism of COAD pathogenesis, progression, and metastasis and finding effective biomarkers and therapeutic targets can greatly help to improve the survival rate of COAD patients. The latest study found that the vast majority of patients with malignancies are chronically immunosuppressed and therefore more susceptible to SARS-CoV-2 infection than healthy individuals, often resulting in a poorer prognosis (Kamboj and Sepkowitz, 2009; Liang et al., 2020). If patients with COAD are infected with the SARS-CoV-2 virus, it may lead to a worse prognosis. Therefore, it is necessary to find drugs that can be used to treat COVID-19/COAD patients.
Traditional Chinese medicine (TCM), derived from more than 2,000 years of clinical experience, has been a crucial role in the treatment of a variety of epidemics and has made significant contributions to the health of the Chinese people. With the outbreak of COVID-19, combinations of TCM treatment protocols have been widely used in China, achieving better clinical outcomes. In China, more than four-fifths of COVID-19 patients have received TCM treatment (Ren et al., 2020). Pueraria Lobata, one of the commonly used TCMs, has be reportedly used for anti-inflammatory, immunomodulatory, anti-cancer and so on (Wang et al., 2020). Puerarin, as the main active ingredient of Pueraria Lobata, has similar pharmacological effects to Pueraria Lobata. As a potent inhibitor of influenza virus neuraminidase, puerarin could reduce influenza virus titers, decrease respiratory inflammatory responses and diminish mortality (Wang et al., 2021). Puerarin nanosuspension showed good anticancer activity and low toxicity and has the potential to become a therapeutic medicine for patients with COAD (Wang et al., 2013). On this basis, we hypothesized that puerarin might be used for COVID-19/COAD patients. Therefore, network pharmacology and other bioinformatic methods were employed to explore the therapeutic targets and pharmacological mechanisms of puerarin in COVID-19/COAD and to elucidate the potential therapeutic value of puerarin.
Materials and Methods
Flowchart of the Study
Figure 1 depicts the entire flow of this study. Firstly, targets associated with COVID-19 were identified though the Gene Expression Omnibus (GEO) dataset and multiple disease-target related databases; then targets associated with COAD were identified by using the TCGA-COAD dataset, the GTEx-colon data dataset, and multiple disease-target related databases; next, multiple drug-target related databases were used to find the therapeutic targets of puerarin. The COVID-19-related targets, COAD-related targets, and therapeutic targets of Puerarin were intersected. Finally, bioinformatics analysis was used to explore the pharmacological mechanisms of puerarin in COVID-19/COAD comorbidity.
FIGURE 1. Flowchart. This study followed the analytical flow of the flowchart to explore the target genes and molecular mechanisms of puerarin in the treatment of COAD/COVID-19.
Identification of Targets Associated With COVID-19 or COAD
Transcriptome sequencing data from the COVDID-19-related datasets (GSE157103, GSE152075, GSE171110) were downloaded though the GEO database (https://www.ncbi.nlm.nih.gov/geo/). The differentially expressed genes (DEGs) associated with COVDID-19 were filtered according to the filtering criteria (p < 0.05 and absolute value of log fold change >1) by using the “limma” package in R. In addition, we also searched seven authoritative databases containing CTD (http://ctdbase.org/), OMIM (https://omim.org/), DisGeNET (https://www.disgenet.org/), GeneCards (https://www.genecards.org), NCBI (https://www.ncbi.nlm.nih.gov/), TTD (http://db.idrblab.net/), pharmaGKB (https://www.pharmgkb.org/) to further explore genes associated with COVID-19. Gene names obtained from the different databases were harmonized. We considered the genes that appeared at least twice as COVID-19-related genes.
Transcriptome sequencing data from the TCGA-COAD and GTEx-CO datasets were obtained via the UCSC Xena (http://xena.ucsc.edu/) database. DEGs associated with COAD (p < 0.05 and absolute value of log fold change >1) were screened by the R language “limma” package. In addition, CTD, OMIM, DisGeNET, GeneCards, NCBI, TTD, and pharmaGKB databases were used to obtain COAD-associated genes. We considered the genes that appeared at least twice as COAD-related genes.
Identification of Puerarin Therapeutic Targets
To fully explore the therapeutic targets of puerarin, we searched seven authoritative databases related to drug targets: 1) PharmMapper (http://www.lilab-ecust.cn/pharmmapper/); 2) CTD (http://ctdbase.org/); 3) SwissTargetPrediction (http://www.swisstargetprediction.ch/); 4) Drug Gene Interaction Database (DGIdb, https://www.dgidb.org/); 5) TCMSP (https://tcmspw.com/tcmsp.php); 6) STITCH (http://stitch.embl.de/); 7) ECTM (http://www.tcmip.cn/ETCM/). To ensure the accuracy of subsequent analyses, targets that appeared at least twice in the database were considered to be therapeutic target genes for puerarin.
Construction of PPI Network and Analysis of Core Genes
We referred to the intersection of therapeutic target genes for puerarin, COVID-19-related genes and COAD-related genes as target genes of puerarin in COVID-19/COAD comorbidity. The protein-protein interaction (PPI) network was constructed by using the GeneMANIA database (http://genemania.org/). The GeneMANIA database is one of the recognized bioinformatics databases that helps to explain the molecular mechanisms involved in the target proteins. Cytoscape software was used to visualize the degree of target genes and to identify the core genes for puerarin treatment of COVID-19/COAD.
Bioinformatics Analysis
The uniport database was used to convert gene names to gene IDs for subsequent bioinformatic analysis. We called the “enrichGO” and “enrichKEGG” functions in R to perform Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis on the target genes for puerarin treatment of COVID-19/COAD, and the “ggplot2” package in R to visualise the analysis results in R (only displayed the 20 results with the lowest q values). Finally, we used the “circlize” package in R to construct the circos plots showing the results of the top five GO and KEGG enrichment analyses with the corresponding genes.
Results
Identification Results of COVID-19-Related and COAD-Related Targets
The study contained 3 GEO datasets associated with COVID-19: GSE152075, GSE157103, and GSE171110. In the GSE157103 dataset, we included sequencing data from 484 nasopharyngeal swab samples (containing 430 COVID-19 positive patients and 54 COVID-19 negative healthy individuals). In the GSE157103 dataset, we included sequencing data from 126 plasma samples (containing 100 COVID-19 positive patients and 26 COVID-19 negative healthy individuals). In the GSE171110 dataset, we included sequencing data from 54 whole blood samples (containing 44 COVID-19 positive patients and 10 COVID-19 negative healthy individuals). The DEGs in the GEO dataset were searched though the “limma” package in R. 2334, 1316, 3302 DEGs were obtained from the GSE152075, GSE157103, GSE171110 datasets respectively (Figures Figure2A–C). In addition, CTD, OMIM, DisGeNET, GeneCards, NCBI, TTD, and pharmaGKB databases were used to find 9,879, 2, 1843, 4,600, 346, 68, and 0 genes associated with COVID-19, respectively (Figure Figure3A, Supplementary Table S1). In the current study, we considered target genes that appeared in at least two databases to be COVID-19-related genes, resulting in 5968 COVID-19-related genes.
FIGURE 2. DEGs from COAD patients or COVID-19 patients (p < 0.05 and absolute value of log fold change >1). The red and blue dots represented genes with increased or decreased expression in COVID-19 patients or COAD patients, respectively, while the black dots represented genes with no significant difference in expression between COVID-19 or COAD patients. (A) Volcano plot of DEGs in the GSE152075 dataset. (B) Volcano plot of DEGs in the GSE157103 dataset. (C) Volcano plot of DEGs in the GSE171110 dataset. (D) Volcano plot of DEGs in the TCGA-COAD dataset and the GTEx-CO dataset.
FIGURE 3. Potential target genes of puerarin in the treatment of COVID-19/COAD comorbidity. (A) Genes associated with COVID-19 from 3 GEO datasets and six databases. (B) Genes associated with COAD from the TCGA-COAD dataset, GTEx-CO dataset, and seven databases. (C) Therapeutic target genes of Puerarin from seven databases. (D) Venn diagram of therapeutic target genes for puerarin, COVID-19-related genes and COAD-related genes.
In this study, the TCGA and GTEx datasets were used to identify DEGs in COAD patients. We combined the data from TCGA-COAD and GTEx-colon and included sequencing data from a total of 349 normal colon samples and 471 COAD samples. The data were analyzed by using the “limma” package in R, resulting in 5,818 DEGs (Figure Figure2D). Through CTD, OMIM, DisGeNET, GeneCards, NCBI, TTD, and pharmaGKB databases, we obtained 29,159, 54, 452, 8,675, 2,260, 3, 56 genes associated with COAD (Figure Figure3B, Supplementary Table S2). We referred to target genes that appeared in at least two databases as COAD-related genes, resulting in 11,755 COAD-related genes.
Identification Results of Puerarin Target Genes
To obtain therapeutic target genes for puerarin, seven databases containing PharmMapper, CTD, SwissTargetPrediction, DGIdb, TCMSP, STITCH, and ECTM were searched, resulting in 268, 124, 102, 2, 53, 26, and six therapeutic target genes respectively (Figure 3C, Supplementary Table S3). The therapeutic target genes were screened by using the presence of at least two databases as a screening condition, resulting in 57 therapeutic targets for puerarin.
Therapeutic Target Genes of Puerarin in COVID-19/COAD
We intersected the therapeutic target genes of puerarin, COAD-related genes, and COVID-19-related genes, and presented them by using Venn diagram (Figure 3D, Supplementary Table S4). We defined these 42 intersecting genes as gene set 1 and used them for subsequent bioinformatics analysis.
PPI Network
We used 42 genes from gene set 1 for PPI network construction. The GeneMANIA database, based on large amounts of genomics and proteomics data, is commonly used for the prediction of gene function and the construction of interactions between genes. In this study, we uploaded gene set 1 to the GeneMANIA database. The results of the analysis were shown in Figure 4A (Supplementary Table S5). The circular nodes represented genes, and the different colors of the node represented gene enriched in different functional analysis results. The different colored lines between the nodes represented the different interactions. The mechanism of puerarin in COVID-19/COAD comorbidity might be related to apoptotic, chemical stress, and oxidative stress. By using Cytoscape software to visualize the degree of each gene in the PPI network (Figure 4B), we found that RELA, BCL2, JUN, FOS, and MAPK1 had the highest degree, suggesting that the five genes might be core targets of puerarin in COVID-19/COAD comorbidity.
FIGURE 4. PPI network of potential therapeutic target genes of puerarin in the treatment of COVID-19/COAD. (A) The PPI network was derived from the GeneMANIA database. The different relationships between nodes were indicated by the different colored connecting lines. Nodes were enriched in different functions indicated by the colors of the nodes. (B) PPI networks constructed from nodal degree values. The degree of the node was proportional to the depth of the node color.
Results of Bioinformatics Analysis
To investigate the pharmacological mechanism of puerarin in the treatment of COVID-19/COAD, we performed bioinformatic analysis of 42 genes in gene set 1. The “cluster profile” package in R was used to perform enrichment analysis of the BP, CC, and MF aspects of GO (Figures 5A–C, Supplementary Table S6). The results of the top five BP, CC, MF, and KEGG enrichment analyses and the corresponding genes were presented using circus plots (Figures 6A–D). The results of the BP enrichment analysis mainly included regulation of apoptotic signaling pathway, response to the metal ion, response to lipopolysaccharide, response to molecule of bacterial origin, regulation of neuron death, productive system development, peptidyl-serine modification, and so on. The results of the CC enrichment analysis found that gene set 1 was mainly involved in caveola, plasma membrane raft, myelin sheath, mitochondrial outer membrane, pseudopodium, organelle outer membrane, protein-DNA complex, and so on. The results of the MF enrichment analysis found that gene set 1 mainly involved MAP kinase activity, phosphatase binding, NF-κB binding, ubiquitin-protein ligase binding, coenzyme binding, histone deacetylase binding, tyrosine kinase activity, etc.
FIGURE 5. GO and KEGG enrichment analysis of potential therapeutic target genes of puerarin for COVID-19/COAD (only the 20 pathways with the smallest p-values were shown). (A), (C) and (E) The results of BP, CC, and MF term enrichment analysis, respectively. (B), (D) and (F) The results of correlation analysis of BP, CC, MF term enrichment analysis, respectively. (G) The results of KEGG enrichment analysis. (H)The results of correlation analysis of KEGG enrichment analysis.
FIGURE 6. Top five enrichment analysis results and corresponding genes. (A), (B) and (C)The top five results of BP, CC, MF enrichment analysis and corresponding genes. (D) The top five results of KEGG enrichment analysis and the corresponding genes.
In addition, the results of the KEGG pathway enrichment analysis mainly included Apoptosis, Colorectal cancer, IL-17 signaling pathway, TNF signaling pathway, Human immunodeficiency virus 1 infection, HIF-1 signaling pathway (Figure 5D, Supplementary Table S7).
Discussion
The SARS-CoV-2 virus has been spreading globally for more than 2 years and has not yet shown signs of improvement, yet we still lack effective drugs against COVID-19. However, in China, protocols that combine Chinese medicine with Western medicine to treat patients with COVID-19 have achieved excellent results. Clinical studies have found that TCM treatment can reduce the rate of severe disease, decrease the use of antibiotics and alleviate clinical symptoms including fever and cough (Huang et al., 2021; Zhang et al., 2021). Cumulatively, one in twenty people worldwide has been infected with the SARS-CoV-2 virus, and patients with COAD are more susceptible to SARS-CoV-2 virus infection than healthy individuals due to their chronically depressed immune system. Based on the antiviral and anticancer properties of puerarin, we combined bioinformatics and network pharmacology approaches to investigate the target genes and mechanisms of puerarin in the treatment of COVID-19/COAD (Wei et al., 2014).
We analyzed data from open databases and defined the 42 potential target genes of puerarin in the treatment of COVID-19/COAD as gene set 1. Through the GeneMANIA database to construct a PPI network and perform functional analysis, we found that gene set 1 mainly has the functions of chemical stress, apoptosis, and oxidative stress. Chemical stress in the body might trigger adverse drug reactions, and adverse reactions during the treatment of COAD patients often led to the inability of patients to tolerate subsequent treatment (Mutter et al., 2015). Apoptosis isthe normal physiological process that occurs within the body in a state of relative equilibrium. Disruption of this balance may lead to disease. A high number of apoptotic cells, for example, could lead to degenerative disorders, whereas a low number of apoptotic cells could contribute to tumor growth and progression (Wong, 2011). But apoptosis is a double-edged sword, and induction of apoptosis in tumor cells is one of the therapeutic mechanisms for tumor patients. Abnormal differentiation and apoptosis of T cells, one of the hallmarks of severe COVID-19, led to immune dysfunction in patients (Leonardi and Proenca, 2020). Persistent oxidative stress might lead to chronic inflammation, which in turn led to the development of diseases including cancer, lung disease, and diabetes (Reuter et al., 2010). Studies have confirmed that oxidative stress is positively related to the incidence of COAD (Barrett et al., 2017). SARS-CoV-2 virus induces the abnormal oxidative stress response in humans, which in turn can lead to the development of the respiratory syndrome and even death in patients. Therefore, drugs that improve the oxidative stress status of patients may improve prognosis and reduce mortality. We screened the genes based on their degree values and identified RELA, BCL2, JUN, FOS, and MAPK1 as possible core target genes for puerarin in treatment of COVID-19/COAD. The NF-κB signaling pathway is involved in several biological processes such as cellular transformation, and immune response in the body. Abnormal NF-κB directly contributed to tumorigenesis, growth, and metastasis, while NF-κB activity was controlled by phosphorylation of upstream genes containing RELA (Lu and Yarbrough, 2015). Inhibition of NF-κB pathway activity reduced the likelihood of severe disease in COVID-19 patients (Hariharan et al., 2021). As an important factor regulating apoptosis, BCL2 has been associated with the development and progression of cancer and autoimmune diseases, and several drugs targeting BCL2 have been developed (Siddiqui et al., 2015). Clinical studies have found that serum concentrations of BCL2 in COVID-19 patients correlate with patient prognosis, with high BCL2 concentrations predicting a better prognosis (Lorente et al., 2021). Activator protein 1, consisting of Jun-Jun dimer or Jun-Fos dimer, regulated the transcription of multiple genes required for biological processes and the development and progression of many tumors, including COAD (Ashida et al., 2005). As one of the complexes of activator protein 1, c-Jun was involved in the crosstalk, integration, and amplification of intercellular signaling pathways, which were involved in cellular activities including survival, apoptosis, proliferation, and tumorigenesis (Meng and Xia, 2011). MAPK signaling pathway controlled cell survival, senescence, and drug resistance by integrating extracellular signals (Braicu et al., 2019). MAPK1 was found to play an important role in the invasion and proliferation of COAD cell line SW480 (Yang et al., 2018). The severe inflammatory response that occurred in severe COVID-19 patients might be associated with activation of the p38 MAKP signaling pathway (Grimes and Grimes, 2020).
The results of bioinformatics analysis revealed that the mechanism of puerarin in COVID-19/COAD might be associated with NF-κB, MAPK, IL-17, TNF, and HIF-1 signaling pathway. Elevated pro-inflammatory cytokines, including IL-17 and TNF-α, were the main causes of disease deterioration in COVID-19 patients (Shibabaw, 2020). IL-17 was involved in the activation of multiple signaling pathways resulting to the emergence of multiple cytokines (e.g. TNF-α, IL-1β) and chemokines, which in turn led to cytokine storm (Ryzhakov et al., 2011; Veldhoen, 2017; Badawi, 2020). Inflammatory cytokines IL-17 and TNF-α upregulated PD-L1 expression in COAD cell line HCT116 through activation of NF-κB signaling pathways, leading to COAD progression (Wang et al., 2017). HIF-1 induced an adaptive response in tumor cells under localized hypoxic conditions in cancer patients and was associated with poor patient prognosis (Lu et al., 2016). Activation of HIF signaling pathway was related to tumor cell invasion and tumor cell drug resistance (Brocato et al., 2014). HIF-1α signaling pathway was involved in SARS-CoV-2 virus infection and subsequent inflammatory response, often leading to poor prognosis in COVID-19 patients (Tian et al., 2021).
Conclusion
In this study, a series of bioinformatics analysis were employed to explore that the potential treatment mechanisms of puerarin in COVID-19/COAD comorbidity. These mechanisms may be related to apoptosis, antiviral, antioxidant, IL-17, TNF, and HIF-1 signaling pathway. The present study provides the preliminary basis and directions for the study of puerarin in COVID-19/COAD. However, as the study is based on transcriptome sequencing data and multiple databases, experimental validation is still lacking. In the follow-up studies, we will use new coronavirus-infected immune cell lines to co-culture with COAD cell lines. The therapeutic effect of puerarin will be observed by simulating the in vivo condition of COVID-19/COAD patients.
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
JX, HW, and WL developed the concept of the project. The data was collected and evaluated by WL, XL, and YY with the help of QM and XW. All authors reviewed and discussed the results and contributed to the paper preparation. JX, HW, WL, and XL wrote the manuscript. All authors have read and approved the final manuscript.
Funding
This study was supportive by Hebei Funding Project for High-level Talents (A202101062); Training Project for Excellent Clinical Medical Talents of Hebei Provincial Department of Finance (361009); Hebei Provincial Health Commission Medical Science Research Key Project Program (20210954); Shenzhen Key Laboratory Foundation (ZDSYS20200811143757022); Project of Department of education of Guangdong province (No. 2021KQNCX074)
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.921517/full#supplementary-material
Supplementary Table S1 | Therapeutic target genes of Puerarin.
Supplementary Table S2 | COVID-19-related genes.
Supplementary Table S3 | COAD-related genes.
Supplementary Table S4 | Potential target genes of puerarin in the treatment of COVID-19/COAD.
Supplementary Table S5 | The results of PPI network.
Supplementary Table S6 | The results of GO enrichment analysis.
Supplementary Table S7 | The results of KEGG enrichment analysis.
Abbreviations
COAD, Colon Adenocarcinoma; PPI, Protein-Protein Interaction; TCM, Traditional Chinese Medicine; GEO, Gene Expression Omnibus; DEGs, Differentially Expressed Genes; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes.
References
Ashida, R., Tominaga, K., Sasaki, E., Watanabe, T., Fujiwara, Y., Oshitani, N., et al. (2005). AP-1 and Colorectal Cancer. Inflammopharmacology 13 (1-3), 113–125. doi:10.1163/156856005774423935
Badawi, A. (2020). Hypercytokinemia and Pathogen-Host Interaction in COVID-19. J. Inflamm. Res. 13, 255–261. doi:10.2147/jir.S259096
Barrett, C. W., Short, S. P., and Williams, C. S. (2017). Selenoproteins and Oxidative Stress-Induced Inflammatory Tumorigenesis in the Gut. Cell Mol. Life Sci. 74 (4), 607–616. doi:10.1007/s00018-016-2339-2
Braicu, C., Buse, M., Busuioc, C., Drula, R., Gulei, D., Raduly, L., et al. (2019). A Comprehensive Review on MAPK: A Promising Therapeutic Target in Cancer. Cancers (Basel) 11 (10), 1618. doi:10.3390/cancers11101618
Brocato, J., Chervona, Y., and Costa, M. (2014). Molecular Responses to Hypoxia-Inducible Factor 1α and beyond. Mol. Pharmacol. 85 (5), 651–657. doi:10.1124/mol.113.089623
Brown, K. G. M., and Solomon, M. J. (2018). Progress and Future Direction in the Management of Advanced Colorectal Cancer. Br. J. Surg. 105 (6), 615–617. doi:10.1002/bjs.10759
Chen, N., Zhou, M., Dong, X., Qu, J., Gong, F., Han, Y., et al. (2020). Epidemiological and Clinical Characteristics of 99 Cases of 2019 Novel Coronavirus Pneumonia in Wuhan, China: a Descriptive Study. Lancet 395 (10223), 507–513. doi:10.1016/s0140-6736(20)30211-7
Durinikova, E., Buzo, K., and Arena, S. (2021). Preclinical Models as Patients' Avatars for Precision Medicine in Colorectal Cancer: Past and Future Challenges. J. Exp. Clin. Cancer Res. 40 (1), 185. doi:10.1186/s13046-021-01981-z
Giacomelli, A., Pezzati, L., Conti, F., Bernacchia, D., Siano, M., Oreni, L., et al. (2020). Self-reported Olfactory and Taste Disorders in Patients With Severe Acute Respiratory Coronavirus 2 Infection: A Cross-Sectional Study. Clin. Infect. Dis. 71 (15), 889–890. doi:10.1093/cid/ciaa330
Grimes, J. M., and Grimes, K. V. (2020). p38 MAPK Inhibition: A Promising Therapeutic Approach for COVID-19. J. Mol. Cell Cardiol. 144, 63–65. doi:10.1016/j.yjmcc.2020.05.007
Hariharan, A., Hakeem, A. R., Radhakrishnan, S., Reddy, M. S., and Rela, M. (2021). The Role and Therapeutic Potential of NF-Kappa-B Pathway in Severe COVID-19 Patients. Inflammopharmacology 29 (1), 91–100. doi:10.1007/s10787-020-00773-9
Huang, K., Zhang, P., Zhang, Z., Youn, J. Y., Wang, C., Zhang, H., et al. (2021). Traditional Chinese Medicine (TCM) in the Treatment of COVID-19 and Other Viral Infections: Efficacies and Mechanisms. Pharmacol. Ther. 225, 107843. doi:10.1016/j.pharmthera.2021.107843
Jayaweera, M., Perera, H., Gunawardana, B., and Manatunge, J. (2020). Transmission of COVID-19 Virus by Droplets and Aerosols: A Critical Review on the Unresolved Dichotomy. Environ. Res. 188, 109819. doi:10.1016/j.envres.2020.109819
Kamboj, M., and Sepkowitz, K. A. (2009). Nosocomial Infections in Patients with Cancer. Lancet Oncol. 10 (6), 589–597. doi:10.1016/s1470-2045(09)70069-5
Lazarus, J. V., Ratzan, S. C., Palayew, A., Gostin, L. O., Larson, H. J., Rabin, K., et al. (2021). A Global Survey of Potential Acceptance of a COVID-19 Vaccine. Nat. Med. 27 (2), 225–228. doi:10.1038/s41591-020-1124-9
Leonardi, A. J., and Proenca, R. B. (2020). Akt-Fas to Quell Aberrant T Cell Differentiation and Apoptosis in Covid-19. Front. Immunol. 11, 600405. doi:10.3389/fimmu.2020.600405
Li, Q., Guan, X., Wu, P., Wang, X., Zhou, L., Tong, Y., et al. (2020). Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia. N. Engl. J. Med. 382 (13), 1199–1207. doi:10.1056/NEJMoa2001316
Liang, W., Guan, W., Chen, R., Wang, W., Li, J., Xu, K., et al. (2020). Cancer Patients in SARS-CoV-2 Infection: a Nationwide Analysis in China. Lancet Oncol. 21 (3), 335–337. doi:10.1016/s1470-2045(20)30096-6
Liu, W. (2021). Long Non-coding RNA VPS9D1-AS1 Promotes Growth of Colon Adenocarcinoma by Sponging miR-1301-3p and CLDN1. Hum. Cell 34 (6), 1775–1787. doi:10.1007/s13577-021-00604-1
Liu, X., Liu, C., Liu, G., Luo, W., and Xia, N. (2020). COVID-19: Progress in Diagnostics, Therapy and Vaccination. Theranostics 10 (17), 7821–7835. doi:10.7150/thno.47987
Lorente, L., Martín, M. M., González-Rivero, A. F., Pérez-Cejas, A., Argueso, M., Perez, A., et al. (2021). Blood Concentrations of Proapoptotic sFas and Antiapoptotic Bcl2 and COVID-19 Patient Mortality. Expert Rev. Mol. Diagn 21 (8), 837–844. doi:10.1080/14737159.2021.1941880
Lu, X., and Yarbrough, W. G. (2015). Negative Regulation of RelA Phosphorylation: Emerging Players and Their Roles in Cancer. Cytokine Growth Factor Rev. 26 (1), 7–13. doi:10.1016/j.cytogfr.2014.09.003
Lu, Y., Wang, B., Shi, Q., Wang, X., Wang, D., and Zhu, L. (2016). Brusatol Inhibits HIF-1 Signaling Pathway and Suppresses Glucose Uptake under Hypoxic Conditions in HCT116 Cells. Sci. Rep. 6, 39123. doi:10.1038/srep39123
Meng, Q., and Xia, Y. (2011). c-Jun, at the Crossroad of the Signaling Network. Protein Cell 2 (11), 889–898. doi:10.1007/s13238-011-1113-3
Mutter, F. E., Park, B. K., and Copple, I. M. (2015). Value of Monitoring Nrf2 Activity for the Detection of Chemical and Oxidative Stress. Biochem. Soc. Trans. 43 (4), 657–662. doi:10.1042/bst20150044
Ren, J.-l., Zhang, A.-H., and Wang, X.-J. (2020). Traditional Chinese Medicine for COVID-19 Treatment. Pharmacol. Res. 155, 104743. doi:10.1016/j.phrs.2020.104743
Reuter, S., Gupta, S. C., Chaturvedi, M. M., and Aggarwal, B. B. (2010). Oxidative Stress, Inflammation, and Cancer: How Are They Linked? Free Radic. Biol. Med. 49 (11), 1603–1616. doi:10.1016/j.freeradbiomed.2010.09.006
Ryzhakov, G., Lai, C. C., Blazek, K., To, K. W., Hussell, T., and Udalova, I. (2011). IL-17 Boosts Proinflammatory Outcome of Antiviral Response in Human Cells. J. Immunol. 187 (10), 5357–5362. doi:10.4049/jimmunol.1100917
Shibabaw, T. (2020). Inflammatory Cytokine: IL-17A Signaling Pathway in Patients Present with COVID-19 and Current Treatment Strategy. J. Inflamm. Res. 13, 673–680. doi:10.2147/jir.S278335
Siddiqui, W. A., Ahad, A., and Ahsan, H. (2015). The Mystery of BCL2 Family: Bcl-2 Proteins and Apoptosis: an Update. Arch. Toxicol. 89 (3), 289–317. doi:10.1007/s00204-014-1448-7
Siegel, R. L., Miller, K. D., Fuchs, H. E., and Jemal, A. (2022). Cancer Statistics, 2022. CA Cancer J. Clin. 72 (1), 7–33. doi:10.3322/caac.21708
Tian, M., Liu, W., Li, X., Zhao, P., Shereen, M. A., Zhu, C., et al. (2021). HIF-1α Promotes SARS-CoV-2 Infection and Aggravates Inflammatory Responses to COVID-19. Signal Transduct. Target Ther. 6 (1), 308. doi:10.1038/s41392-021-00726-w
Veldhoen, M. (2017). Interleukin 17 Is a Chief Orchestrator of Immunity. Nat. Immunol. 18 (6), 612–621. doi:10.1038/ni.3742
Wang, H. X., Zeng, M. S., Ye, Y., Liu, J. Y., and Xu, P. P. (2021). Antiviral Activity of Puerarin as Potent Inhibitor of Influenza Virus Neuraminidase. Phytother. Res. 35 (1), 324–336. doi:10.1002/ptr.6803
Wang, S., Zhang, S., Wang, S., Gao, P., and Dai, L. (2020). A Comprehensive Review on Pueraria: Insights on its Chemistry and Medicinal Value. Biomed. Pharmacother. 131, 110734. doi:10.1016/j.biopha.2020.110734
Wang, X., Yang, L., Huang, F., Zhang, Q., Liu, S., Ma, L., et al. (2017). Inflammatory Cytokines IL-17 and TNF-α Up-Regulate PD-L1 Expression in Human Prostate and Colon Cancer Cells. Immunol. Lett. 184, 7–14. doi:10.1016/j.imlet.2017.02.006
Wang, Y., Ma, Y., Zheng, Y., Song, J., Yang, X., Bi, C., et al. (2013). In Vitro and In Vivo Anticancer Activity of a Novel Puerarin Nanosuspension against Colon Cancer, with High Efficacy and Low Toxicity. Int. J. Pharm. 441 (1-2), 728–735. doi:10.1016/j.ijpharm.2012.10.021
Wei, S. Y., Chen, Y., and Xu, X. Y. (2014). Progress on the Pharmacological Research of Puerarin: a Review. Chin. J. Nat. Med. 12 (6), 407–414. doi:10.1016/s1875-5364(14)60064-9
Wong, R. S. (2011). Apoptosis in Cancer: from Pathogenesis to Treatment. J. Exp. Clin. Cancer Res. 30 (1), 87. doi:10.1186/1756-9966-30-87
Xu, J., Zhao, S., Teng, T., Abdalla, A. E., Zhu, W., Xie, L., et al. (2020). Systematic Comparison of Two Animal-To-Human Transmitted Human Coronaviruses: SARS-CoV-2 and SARS-CoV. Viruses 12 (2), 244. doi:10.3390/v12020244
Yang, Y., Li, X. J., Li, P., and Guo, X. T. (2018). MicroRNA-145 Regulates the Proliferation, Migration and Invasion of Human Primary Colon Adenocarcinoma Cells by Targeting MAPK1. Int. J. Mol. Med. 42 (6), 3171–3180. doi:10.3892/ijmm.2018.3904
Keywords: puerarin, COVID-19, COAD, comorbidity, network pharmacology, bioinformatics
Citation: Liang W, Li X, Yao Y, Meng Q, Wu X, Wang H and Xue J (2022) Puerarin: A Potential Therapeutic for Colon Adenocarcinoma (COAD) Patients Suffering From SARS-CoV-2 Infection. Front. Pharmacol. 13:921517. doi: 10.3389/fphar.2022.921517
Received: 16 April 2022; Accepted: 06 May 2022;
Published: 23 May 2022.
Edited by:
Yongye Huang, Northeastern University, ChinaReviewed by:
Bowang Chen, Chinese Academy of Medical Sciences and Peking Union Medical College, ChinaYuting Chen, Genesis Healthcare, Japan
Copyright © 2022 Liang, Li, Yao, Meng, Wu, Wang and Xue. 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: Weizheng Liang, am1iYjEyMDNAMTI2LmNvbQ==; Hao Wang, aGFvd2FuZzA4MDZAZ21haWwuY29t; Jun Xue, eHVlanVuQGhlYmVpbnUuZWR1LmNu
†These authors have contributed equally to this work