Skip to main content

Identification of intestinal flora-related key genes and therapeutic drugs in colorectal cancer



Colorectal cancer (CRC) is a multifactorial tumor and a leading cause of cancer-specific deaths worldwide. Recent research has shown that the alteration of intestinal flora contributes to the development of CRC. However, the molecular mechanism by which intestinal flora influences the pathogenesis of CRC remains unclear. This study aims to explore the key genes underlying the effect of intestinal flora on CRC and therapeutic drugs for CRC.


Intestinal flora-related genes were determined using text mining. Based on The Cancer Genome Atlas database, differentially expressed genes (DEGs) between CRC and normal samples were identified with the limma package of the R software. Then, the intersection of the two gene sets was selected for enrichment analyses using the tool Database for Annotation, Visualization and Integrated Discovery. Protein interaction network analysis was performed for identifying the key genes using STRING and Cytoscape. The correlation of the key genes with overall survival of CRC patients was analyzed. Finally, the key genes were queried against the Drug-Gene Interaction database to find drug candidates for treating CRC.


518 genes associated with intestinal flora were determined by text mining. Based on The Cancer Genome Atlas database, we identified 48 DEGs associated with intestinal flora, including 25 up-regulated and 23 down-regulated DEGs in CRC. The enrichment analyses indicated that the selected genes were mainly involved in cell–cell signaling, immune response, cytokine-cytokine receptor interaction, and JAK-STAT signaling pathway. The protein–protein interaction network was constructed with 13 nodes and 35 edges. Moreover, 8 genes in the significant cluster were considered as the key genes and chemokine (C-X-C motif) ligand 8 (CXCL8) correlated positively with the overall survival of CRC patients. Finally, a total of 24 drugs were predicted as possible drugs for CRC treatment using the Drug-Gene Interaction database.


These findings of this study may provide new insights into CRC pathogenesis and treatments. The prediction of drug-gene interaction is of great practical significance for exploring new drugs or novel targets for existing drugs.

Peer Review reports


Colorectal cancer (CRC) has the third-highest incidence rate and the second-highest mortality rate among all types of cancers worldwide, according to the 2018 global cancer statistics [1]. The World Health Organization estimates that the prevalence of CRC will rise to 3093 million by 2040 [2]. Although great effort has been made to improve clinical treatment over the past few decades, CRC still poses an enormous threat to human health. For example, the median overall survival rarely exceeds 6.8 months for patients with metastatic CRC [3, 4], and the 5-year survival rate remains only 10% [5]. Therefore, targeting novel pathways is indispensable for further improvement in the prognosis of patients with CRC.

CRC is generally considered to be a multifactorial disease involving diet, inflammatory process, genetic alteration, and environmental factor [6, 7]. The human intestinal tract contains 10 to 100 trillion florae, which is 10 times more than the number of total human body cells. There is increasing evidence that patients with CRC harbor a distinct microbiota. Compared with healthy individuals, the probiotics, including Bifidobacterium and Lactobacillus acidophilus, are decreased while the pathogenic bacteria, including Bacteroides/Prevotella and Enterococcus faecalis, are increased in patients with CRC [8,9,10,11]. Previous studies have confirmed the role of intestinal flora as an environmental factor for the carcinogenesis of CRC [12, 13]. Intestinal flora dysfunction can induce abnormal immune reactions, resulting in a special immune microenvironment in colorectal tissue and producing carcinogenic metabolites that induce DNA damage and gene mutations in host cells [14]. For example, alterations in the microbiota can drive the upregulation of interleukin-17c in intestinal epithelial cells during intestinal tumorigenesis. Microbiota-driven IL-17c promotes cell survival and tumorigenesis by inducing the expression of Bcl-2 and Bcl-xl in intestinal epithelial cells [15]. The enzymes produced by Clostridium bacteria induce increased secondary bile acids that lead to ROS production and DNA damage. These changes can lead to non-integer multiples, KRAS mutations, and micronuclei [16]. However, the current knowledge on the mechanism by which intestinal flora influences the pathogenesis of CRC remains scarce.

Due to the existence of tens of thousands of biomedical journals, the field of biomedicine is flooded with new articles, leading to information overload. Therefore, there is increasing interest in text mining, an efficient approach that enables the identification of biologically relevant entities, such as genes and diseases, complex biological relations, and comprehensive networks. Researches in multifarious fields have been performed using text mining, ranging from pattern-matching methods to molecular events extraction [17,18,19]. Yet to date, the approaches of text mining have not addressed how intestinal flora has an impact on the development of CRC.

In this study, we extracted the genes associated with intestinal flora using text mining. Then, we generated common genes by combining these extracted genes and differentially expressed genes (DEGs) between CRC and normal samples. With further analyses of the functional enrichment and protein–protein interaction (PPI), we identified 8 potential target genes. Finally, candidate drugs for CRC treatment were derived using the Drug-Gene Interaction database (DGIdb).


Gene collection

The web server GenCLiP 3 ( was used to perform text mining. The inclusion criteria for screening were: (i) literature published from 1980 to 2020; (ii) human genes that co-occur with “gut microbiota” or “intestinal flora” in sentence (iii) genes in MEDLINE. Subsequently, a list of genes associated with intestinal flora was extracted. To identify the DEGs between CRC and normal samples from The Cancer Genome Atlas (TCGA) database, the matrix data of gene expression levels were analyzed with the limma package of the R software. Log2 (fold change) > 2 or < − 2 and false discovery rate < 0.05 were set as the threshold for filtering DEGs. VennDiagram package was used to identify the intersection of the two aforementioned gene sets.

Functional enrichment analyses

Using the online tool Database for Annotation, Visualization and Integrated Discovery (DAVID), enrichment analyses of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were performed for the intersection genes. P value < 0.05 was considered as statistically significant.

PPI network analysis

The PPI network of the selected common genes was constructed using the online STRING ( and then was visualized using Cytoscape v3.7.2. A confidence score of 0.4 or higher was set as the cut-off criterion. The Molecular Complex Detection (MCODE) plug-in of Cytoscape was used to classify the significant node clusters with the default. The genes in the most significant cluster were selected as key genes for overall survival and drug-gene interaction analyses.

Correlation of the key genes with overall survival

Gene Expression Profiling Interactive Analysis (GEPIA) is a web-based tool that is freely available to users [20]. Based on TCGA data, GEPIA contains multiple modules, including analyses of the relationship between gene expression and overall survival. We used the GEPIA to analyze the correlation of the key genes with overall survival of CRC patients.

Drug-gene interaction

The screened key genes were served as the potential targets for searching drugs. DGIdb database (v3.0,, a freely available database, was used to identify known and potential drug-gene interactions. Finally, the FDA approved antineoplastic drugs were considered as the potential drugs for CRC treatment.


Identification of intestinal flora-related DEGs in CRC

According to the inclusion criteria, 518 genes were identified within 1110 published papers. To explore whether intestinal flora influences the pathogenesis of CRC by regulating the gene expression, we then identified the DEGs between CRC and normal samples from the TCGA database. Compared with the normal samples, there were 1665 and 2235 DEGs in colon adenocarcinoma (COAD) and rectal adenocarcinoma (READ), respectively. Moreover, there were 25 up-regulated genes (Fig. 1a) and 23 down-regulated genes (Fig. 1b) in common between the DEGs and intestinal flora-related genes. The 48 common genes listed in Table 1 were selected for further analysis.

Fig. 1

Intersection between intestinal flora-related genes and DEGs in CRC. a Venn diagrams showing the overlapping up-regulated genes among intestinal flora-related genes, DEGs in COAD, and DEGs in READ. b Venn diagrams showing the overlapping down-regulated genes among intestinal flora-related genes, DEGs in COAD, and DEGs in READ. CRC, colorectal cancer; DEGs, differentially expressed genes; COAD, colon adenocarcinoma; READ, rectal adenocarcinoma

Table 1 Genes in common between the DEGs and intestinal flora-related genes

Functional enrichment analyses

Using the DAVID tool, we identified 55 GO terms and 5 KEGG terms in which the 48 intestinal flora-related DEGs enriched significantly (P < 0.05). The results showed the top 5 significant enrichment terms for biological processes, cellular component and molecular function, and 5 KEGG pathway terms (Fig. 2). As listed in Table 2, in the biological processes annotation, selected DEGs are mainly involved in cell–cell signaling, immune response, inflammatory response, digestion, and G-protein coupled receptor signaling pathway. In the cellular component annotation, DEGs are mainly involved in extracellular space, extracellular region, plasma membrane, extracellular exosome, and apical plasma membrane. As for molecular function, DEGs are mainly involved in hormone activity, growth factor activity, cytokine activity, receptor binding, and chemokine activity.

Fig. 2

Functional enrichment analyses of the 48 intestinal flora-related DEGs. The horizontal axis denotes different KEGG pathways and GO terms, including molecular function (MF), cellular component (CC), and biological process (BP). The vertical axis denotes the number of enriched genes and significant differences. DEGs, differentially expressed genes; GO, gene ontology; KEGG, Kyoto Gene and Genomic Encyclopedia

Table 2 The top five GO enrichment terms of the intestinal flora-related DEGs

Of the 5 enriched pathways, the two most significant pathways were cytokine-cytokine receptor interaction (P = 2.22E−05) and JAK-STAT signaling pathway (P = 0.006). Additional highly enriched relevant pathways were rheumatoid arthritis, fat digestion and absorption, and TNF signaling pathway (Table 3).

Table 3 The KEGG pathway enrichment terms of the intestinal flora-related DEGs

PPI network construction and survival analysis

From the KEGG pathway analysis, 13 genes were selected. With these 13 genes, the PPI network analysis was performed using STRING and then was visualized using Cytoscape v3.7.2 (Fig. 3a). There were 13 nodes and 35 edges in the PPI network (PPI enrichment P < 1.0E−16). Moreover, the most significant cluster network (score = 7.143) was created using MCODE plug-in, consisting of 25 edges and 8 nodes/genes (CSF2, CXCL1, MMP3, CXCL10, IL23A, CXCL8, IL22, and IL6R) (Fig. 3b).

Fig. 3

PPI network and survival analysis. a Based on the STRING online database, 13 intestinal flora-related DEGs were filtered into PPI network. b The most significant module from the PPI network. c Kaplan–Meier curves of overall survival for patients grouped by expression level of CXCL8. PPI, protein–protein interaction; DEGs, differentially expressed genes; CXCL8, chemokine (C-X-C motif) ligand 8

Further, we used the GEPIA to analyze the correlation of the 8 key genes with overall survival of CRC patients. The result showed that only chemokine (C-X-C motif) ligand 8 (CXCL8) was closely related to the overall survival of patients with CRC (Fig. 3c).

Drug‑gene interaction

The 8 potential key genes in the PPI network were selected for drug-gene interaction analysis. A list of 24 drugs meeting the requirements for CRC treatment was compiled, comprising only the antineoplastic drugs that had been approved by the FDA (Table 4). As listed in Table 4, the potential targets of these drugs include CSF2, CXCL8, IL6R, and CXCL10, and 79.2% (19/24) of the drugs target CSF2 and CXCL8. Several drugs have been applied either alone or in combination in clinical trials or treatments for CRC, such as cetuximab, oxaliplatin, bevacizumab, temozolomide, and interferon alfa-2b.

Table 4 Candidate drugs targeting key genes


The intestinal flora plays significant roles in the formation and development of CRC via producing carcinogenic toxins and metabolites, causing intestinal dysbiosis, and altering immune response [14]. However, current knowledge of the specific mechanism by which intestinal flora influences the pathogenesis of CRC remains limited. In this study, we utilized text mining to extract the intestinal flora-related genes. With further analyses of the functional enrichments and PPI network, we identified 8 key genes associated with intestinal flora and the development of CRC. Finally, candidate drugs targeting 4 key genes (CSF2, CXCL8, IL6R, and CXCL10) were derived using the DGIdb database.

To explore whether intestinal flora contributes to the development of CRC via regulating gene expression, we identified the DEGs between CRC and normal samples. The interaction of the DEGs and intestinal flora-related genes extracted by text mining was selected for further functional enrichments and PPI analyses. Functional enrichment of the 48 selected DEGs in the GO biological process and KEGG pathways highlight their role in immune response, inflammatory response, and intestinal function. Previous studies have confirmed that alteration of intestinal flora can promote CRC development by inducing inflammation, immune suppression, and attacking the gut barrier system [21,22,23]. Moreover, intestinal flora has been proven to promote tumorigenesis and chemoresistance of CRC via regulating gene expression [24,25,26]. Therefore, the mechanism of intestinal flora contributing to the development of CRC is complex and multifaceted.

Using the MCODE plug-in, we created the most significant cluster network, including 8 genes. 87.5% (7/8) of these genes enriched in the cytokine-cytokine receptor interaction pathway, indicating that interactions between the intestinal flora and immune system contribute to the development of CRC.

Drug resistance has been a Gordian knot in the treatment of cancer. Therefore, we identified a list of 24 drugs with the potential therapeutic efficacy against CRC. Among the 8 key genes, the potential gene targets of these drugs are CSF2, CXCL8, IL6R, and CXCL10, and most of the drugs were CSF2 and CXCL8 inhibitors. As a cytokine, CSF2 can stimulate the recruitment and maturation of dendritic cells to induce protective immunity and then exert anti-tumor effects [27]. Yet in the tumor microenvironment, CSF2 is often up-regulated and suppress the immune response, resulting in a poor prognosis for patients [27, 28]. Another cytokine, CXCL8, has been proven to be associated with chemoresistance of CRC [29]. The increased expression of CXCL8 induced by anti-cancer drugs, such as doxorubicin and cisplatin, can upregulate the expression of ATP-binding cassette transporters, resulting in poor chemotherapeutic response [29]. Moreover, a high level of CXCL8 predicted poor overall survival in patients with CRC [30], which is contrary to our findings in this study. This discrepancy probably attributes to the expression characteristics of CXCL8 in different cells within CRC tissues. The tumor microenvironment is a complex system that contains infiltrating immune cells, epithelial cells, and fibroblasts, as well as tumor cells. In addition, CXCL8 is a pro-inflammatory cytokine produced by tumor cells, neutrophils, and endothelial cells [31]. Oladipo et al. found that CRC patients with CXCL8 positivity in the tumor-infiltrating cells had a significantly improved prognosis compared with patients with negativity [32]. Therefore, we speculated that CXCL8 expression in infiltrating cells was dominant in CRC tissues that were obtained from the TCGA database. Among the listed drugs, cetuximab was considered as a prospective drug for CRC therapy thanks to its ability to decrease CXCL8 expression other than inhibiting the epidermal growth factor receptor [33, 34].

IL6R is the receptor of IL-6, which forms a dimer with glycoprotein-130. IL-6 binds the IL6R to initiate the IL-6 signaling that transduces intracellular signals via activation of the JAK-STAT3 pathway [35]. Research has suggested that the IL-6 signaling pathway plays an important role in the development and chemoresistance of various cancers, including CRC [36,37,38]. Accordingly, IL6R has been proposed as a promising target for CRC treatment. IL6R antagonist antibody, tocilizumab, could significantly reduce viability and enhance the apoptosis of CRC cells by blocking the IL-6/STAT3 pathway [39].

CXCL10 is a member of interferon-inducible proteins, which is increasingly being considered as a pro-tumorigenic factor in various cancers, including CRC [40]. Besides, elevated serum CXCL10 was associated with liver metastasis and poor prognosis in CRC [41]. Thus, CXCL10 may be a potential therapeutic target for CRC.


In conclusion, we presented a novel method to explore the molecular mechanism underlying the effect of intestinal flora on CRC. Importantly, we identified 8 potential key genes and 24 candidate drugs. Ten of the 24 drugs have not been tested in CRC, which not only provides a theoretical basis for new trials but also provides new insights into targeting drug discovery. However, our study has some limitations. First, a limitation of the present study is its retrospective nature. All the findings were generated based on the published literature and public database. More prospective studies should be required to verify our findings. Second, it was not clear whether the decrease in probiotics or the increase in pathogenic bacteria caused the DEGs in CRC. Therefore, further mechanistic study of intestinal flora-related genes is encouraged.

Availability of data and materials

A list of genes associated with intestinal flora was extracted using the publicly available GenCLiP3 database ( The microarray studies of CRC analyzed during the current study were available in The Cancer Genome Atlas (TCGA) (dataset ID: TCGA-COAD,; dataset ID: TCGA-READ, The online tool DAVID ( was used to perform functional enrichment analyses. The PPI network of the selected common genes was constructed using the online STRING ( Gene Expression Profiling Interactive Analysis (GEPIA) tool ( to analyze the correlation of the key genes with overall survival of CRC patients. DGIdb database ( was used to identify known and potential drug-gene interactions.



Colorectal cancer


Differentially expressed genes


The Cancer Genome Atlas


Chemokine (C-X-C motif) ligand 8


Protein–protein interaction


Drug-gene interaction database


Database for annotation, visualization and integrated discovery


Enrichment analyses of gene ontology


Kyoto encyclopedia of genes and genomes


Molecular complex detection


Gene expression profiling interactive analysis


Colon adenocarcinoma


Rectal adenocarcinoma


  1. 1.

    Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.

    Article  Google Scholar 

  2. 2.

    Marcuello M, Vymetalkova V, Neves RPL, Duran-Sanchon S, Vedeld HM, Tham E, et al. Circulating biomarkers for early detection and clinical management of colorectal cancer. Mol Aspects Med. 2019;69:107–22.

    CAS  Article  Google Scholar 

  3. 3.

    Chen EX, Jonker DJ, Loree JM, Kennecke HF, Berry SR, Couture F et al. Effect of Combined immune checkpoint inhibition vs. best supportive care alone in patients with advanced colorectal cancer: The Canadian Cancer Trials Group CO.26 Study. JAMA Oncol 2020, 6(6):831–838.

  4. 4.

    Kim RD, Azad NS, Morse MA, Poplin E, Mahipal A, Tan B Jr, et al. Phase II study of ensituximab, a novel chimeric monoclonal antibody, in adults with unresectable. Metastatic Colorectal Cancer Clin Cancer Res. 2020;26(14):3557–64.

    CAS  PubMed  Google Scholar 

  5. 5.

    Lafitte M, Sirvent A, Roche S. Collagen kinase receptors as potential therapeutic targets in metastatic colon cancer. Front Oncol. 2020;10:125.

    Article  Google Scholar 

  6. 6.

    Ibanez-Sanz G, Diez-Villanueva A, Alonso MH, Rodriguez-Moranta F, Perez-Gomez B, Bustamante M, et al. Risk model for colorectal cancer in spanish population using environmental and genetic factors: results from the MCC-Spain study. Sci Rep. 2017;7:43263.

    CAS  Article  Google Scholar 

  7. 7.

    Manzat-Saplacan RM, Mircea PA, Balacescu L, Chira RI, Berindan-Neagoe I, Balacescu O. Can we change our microbiome to prevent colorectal cancer development? Acta Oncol. 2015;54(8):1085–95.

    Article  Google Scholar 

  8. 8.

    Kosumi K, Hamada T, Koh H, Borowsky J, Bullman S, Twombly TS, et al. The amount of bifidobacterium genus in colorectal carcinoma tissue in relation to tumor characteristics and clinical outcome. Am J Pathol. 2018;188(12):2839–52.

    CAS  Article  Google Scholar 

  9. 9.

    Li SC, Lin HP, Chang JS, Shih CK. Lactobacillus acidophilus-Fermented Germinated Brown Rice Suppresses Preneoplastic Lesions of the Colon in Rats. Nutrients 2019, 11(11).

  10. 10.

    Sobhani I, Tap J, Roudot-Thoraval F, Roperch JP, Letulle S, Langella P, et al. Microbial dysbiosis in colorectal cancer (CRC) patients. PLoS ONE. 2011;6(1):e16393.

    CAS  Article  Google Scholar 

  11. 11.

    Geravand M, Fallah P, Yaghoobi MH, Soleimanifar F, Farid M, Zinatizadeh N, et al. Investigation of enterococcus faecalis population in patients with polyp and colorectal cancer in comparison of healthy individuals. Arq Gastroenterol. 2019;56(2):141–5.

    Article  Google Scholar 

  12. 12.

    Lin C, Cai X, Zhang J, Wang W, Sheng Q, Hua H, et al. Role of gut microbiota in the development and treatment of colorectal cancer. Digestion. 2019;100(1):72–8.

    Article  Google Scholar 

  13. 13.

    Song M, Chan AT, Sun J. Influence of the gut microbiome, diet, and environment on risk of colorectal cancer. Gastroenterology. 2020;158(2):322–40.

    CAS  Article  Google Scholar 

  14. 14.

    Si H, Yang Q, Hu H, Ding C, Wang H, Lin X. Colorectal cancer occurrence and treatment based on changes in intestinal flora. Semin Cancer Biol 2020.

  15. 15.

    Song X, Gao H, Lin Y, Yao Y, Zhu S, Wang J, et al. Alterations in the microbiota drive interleukin-17C production from intestinal epithelial cells to promote tumorigenesis. Immunity. 2014;40(1):140–52.

    CAS  Article  Google Scholar 

  16. 16.

    Ajouz H, Mukherji D, Shamseddine A. Secondary bile acids: an underrecognized cause of colon cancer. World J Surg Oncol. 2014;12:164.

    Article  Google Scholar 

  17. 17.

    Liu H, Hunter L, Keselj V, Verspoor K. Approximate subgraph matching-based literature mining for biomedical events and relations. PLoS ONE. 2013;8(4):e60954.

    CAS  Article  Google Scholar 

  18. 18.

    Verspoor KM, Heo GE, Kang KY, Song M. Establishing a baseline for literature mining human genetic variants and their relationships to disease cohorts. BMC Med Inform Decis Mak. 2016;16(Suppl 1):68.

    Article  Google Scholar 

  19. 19.

    Kveler K, Starosvetsky E, Ziv-Kenet A, Kalugny Y, Gorelik Y, Shalev-Malul G, et al. Immune-centric network of cytokines and cells in disease context identified by computational mining of PubMed. Nat Biotechnol. 2018;36(7):651–9.

    CAS  Article  Google Scholar 

  20. 20.

    Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–102.

    CAS  Article  Google Scholar 

  21. 21.

    Tilg H, Adolph TE, Gerner RR, Moschen AR. The intestinal microbiota in colorectal cancer. Cancer Cell. 2018;33(6):954–64.

    CAS  Article  Google Scholar 

  22. 22.

    Sina C, Kemper C, Derer S. The intestinal complement system in inflammatory bowel disease: Shaping intestinal barrier function. Semin Immunol. 2018;37:66–73.

    CAS  Article  Google Scholar 

  23. 23.

    Wu J, Li Q, Fu X. Fusobacterium nucleatum contributes to the carcinogenesis of colorectal cancer by inducing inflammation and suppressing host immunity. Transl Oncol. 2019;12(6):846–51.

    Article  Google Scholar 

  24. 24.

    Zhang Y, Wang XL, Zhou M, Kang C, Lang HD, Chen MT, et al. Crosstalk between gut microbiota and Sirtuin-3 in colonic inflammation and tumorigenesis. Exp Mol Med. 2018;50(4):21.

    Article  Google Scholar 

  25. 25.

    Cremonesi E, Governa V, Garzon JFG, Mele V, Amicarella F, Muraro MG, et al. Gut microbiota modulate T cell trafficking into human colorectal cancer. Gut. 2018;67(11):1984–94.

    CAS  Article  Google Scholar 

  26. 26.

    Zhang S, Yang Y, Weng W, Guo B, Cai G, Ma Y, et al. Fusobacterium nucleatum promotes chemoresistance to 5-fluorouracil by upregulation of BIRC3 expression in colorectal cancer. J Exp Clin Cancer Res. 2019;38(1):14.

    Article  Google Scholar 

  27. 27.

    Xu Z, Zhang Y, Xu M, Zheng X, Lin M, Pan J, et al. Demethylation and overexpression of CSF2 are involved in immune response, chemotherapy resistance, and poor prognosis in colorectal cancer. Onco Targets Ther. 2019;12:11255–69.

    CAS  Article  Google Scholar 

  28. 28.

    Lee YY, Wu WJ, Huang CN, Li CC, Li WM, Yeh BW, et al. CSF2 Overexpression is associated with STAT5 phosphorylation and poor prognosis in patients with urothelial carcinoma. J Cancer. 2016;7(6):711–21.

    CAS  Article  Google Scholar 

  29. 29.

    Du J, He Y, Li P, Wu W, Chen Y, Ruan H. IL-8 regulates the doxorubicin resistance of colorectal cancer cells via modulation of multidrug resistance 1 (MDR1). Cancer Chemother Pharmacol. 2018;81(6):1111–9.

    CAS  Article  Google Scholar 

  30. 30.

    Cheng XS, Li YF, Tan J, Sun B, Xiao YC, Fang XB, et al. CCL20 and CXCL8 synergize to promote progression and poor survival outcome in patients with colorectal cancer by collaborative induction of the epithelial-mesenchymal transition. Cancer Lett. 2014;348(1–2):77–87.

    CAS  Article  Google Scholar 

  31. 31.

    Alfaro C, Sanmamed MF, Rodriguez-Ruiz ME, Teijeira A, Onate C, Gonzalez A, et al. Interleukin-8 in cancer pathogenesis, treatment and follow-up. Cancer Treat Rev. 2017;60:24–31.

    CAS  Article  Google Scholar 

  32. 32.

    Oladipo O, Conlon S, O’Grady A, Purcell C, Wilson C, Maxwell PJ, et al. The expression and prognostic impact of CXC-chemokines in stage II and III colorectal cancer epithelial and stromal tissue. Br J Cancer. 2011;104(3):480–7.

    CAS  Article  Google Scholar 

  33. 33.

    Janjigian YY, Smit EF, Groen HJ, Horn L, Gettinger S, Camidge DR, et al. Dual inhibition of EGFR with afatinib and cetuximab in kinase inhibitor-resistant EGFR-mutant lung cancer with and without T790M mutations. Cancer Discov. 2014;4(9):1036–45.

    CAS  Article  Google Scholar 

  34. 34.

    Paul T, Schumann C, Rudiger S, Boeck S, Heinemann V, Kachele V, et al. Cytokine regulation by epidermal growth factor receptor inhibitors and epidermal growth factor receptor inhibitor associated skin toxicity in cancer patients. Eur J Cancer. 2014;50(11):1855–63.

    CAS  Article  Google Scholar 

  35. 35.

    Harrison SC, Smith AJ, Jones GT, Swerdlow DI, Rampuri R, Bown MJ, et al. Interleukin-6 receptor pathways in abdominal aortic aneurysm. Eur Heart J. 2013;34(48):3707–16.

    CAS  Article  Google Scholar 

  36. 36.

    Tseng-Rogenski SS, Hamaya Y, Choi DY, Carethers JM. Interleukin 6 alters localization of hMSH3, leading to DNA mismatch repair defects in colorectal cancer cells. Gastroenterology. 2015;148(3):579–89.

    CAS  Article  Google Scholar 

  37. 37.

    Rokavec M, Oner MG, Li H, Jackstadt R, Jiang L, Lodygin D, et al. IL-6R/STAT3/miR-34a feedback loop promotes EMT-mediated colorectal cancer invasion and metastasis. J Clin Invest. 2014;124(4):1853–67.

    CAS  Article  Google Scholar 

  38. 38.

    Moon SU, Kang MH, Sung JH, Kim JW, Lee JO, Kim YJ, et al. Effect of Smad3/4 on chemotherapeutic drug sensitivity in colorectal cancer cells. Oncol Rep. 2015;33(1):185–92.

    CAS  Article  Google Scholar 

  39. 39.

    Wang J, Zhou J, Jiang C, Zheng J, Namba H, Chi P, et al. LNRRIL6, a novel long noncoding RNA, protects colorectal cancer cells by activating the IL-6-STAT3 pathway. Mol Oncol. 2019;13(11):2344–60.

    CAS  Article  Google Scholar 

  40. 40.

    Liu M, Guo S, Stiles JK. The emerging role of CXCL10 in cancer (Review). Oncol Lett. 2011;2(4):583–9.

    CAS  Article  Google Scholar 

  41. 41.

    Toiyama Y, Fujikawa H, Kawamura M, Matsushita K, Saigusa S, Tanaka K, et al. Evaluation of CXCL10 as a novel serum marker for predicting liver metastasis and prognosis in colorectal cancer. Int J Oncol. 2012;40(2):560–6.

    CAS  PubMed  Google Scholar 

Download references


Not applicable.


Not applicable.

Author information




YL and TZ conceived and designed the study and prepared the manuscript. JZ, HZ, and FL were responsible for the literature search, data visualization, and analysis. ZS reviewed and revised the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Yezhou Li or Tiancheng Zhao.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhang, J., Zhang, H., Li, F. et al. Identification of intestinal flora-related key genes and therapeutic drugs in colorectal cancer. BMC Med Genomics 13, 172 (2020).

Download citation


  • Colorectal cancer
  • Intestinal flora
  • Text mining
  • Key genes
  • Drugs