Yu Wang, Hei-Ying Jin, Ming-Zhi Fang, Xiao-Feng Wang, Hao Chen, Shu-Liang Huang, De-Song Kong, Min Li,Xiu Zhang, Yu Sun, Shui-Ming Wang
Abstract
Key words: Epigallocatechin gallate; Dimethylhydrazine; Colorectal cancer; Aberrant crypt foci; Mitogen-activated protein kinase; The Cancer Genome Atlas
Colorectal cancer (CRC) has now become one of the most common malignancies of the digestive tract in the world[1]. In 2018, there were 1096601 new cases of colon cancer and 704376 new cases of rectal cancer in 185 major countries all over the world.In the same year, about 551269 and 310394 people died of colon and rectal cancer,respectively[2]. Even so, the morbidity and mortality of CRC have been declining slowly and steadily, which could be because of emergence of advanced detection methods, such as colonoscopy, flexible sigmoidoscopy, fecal occult blood test, fecal DNA examination, and barium enema X-ray examination. These have become part of routine cancer screening to detect the occurrence of early intestinal cancer and even polyps[3-5]. Treatment is as important as diagnosis. Although we have many effective treatment methods for CRC, including surgery, chemotherapy, radiotherapy, targeted therapy, and the recently popular immunotherapy techniques, these treatment methods are associated with many adverse effects that might necessitate treatment interruption or cause disease progression[6].
Green tea has gained popularity in the daily health care[7,8]. It contains more than ten effective polyphenols, including catechin, gallic acid, caffeic acid, epicatechin,gallocatechin, catechin gallate, epigallocatechin, epicatechin gallate, and epigallocatechin gallate (EGCG)[8]. When combined with traditional radiation and chemotherapeutic agents, synergistic effects with reduced toxicity have been reported[9]. EGCG, the most abundant polyphenol (50%) in green tea, has the highest activity[10]. Recent studies have shown that EGCG can inhibit tumor cell proliferation and metastasis and promote tumor cell apoptosis through various signaling pathwaysin vivoandin vitro[11-13]. EGCG probably exerted an anticancer effect,i.e., it halted the proliferation of PC-9 and A431 cells by arresting their cell cycle in the G2/M and G0/G1 phases, respectively[14,15]. EGCG also inhibits the proliferation of LoVo cells and induces their apoptosis by shifting the cell population from the S and M/G2 phases to the G1 phase[16]. It inhibits the growth of CRC cells by inhibiting the activation of receptor tyrosine kinases. In addition, EGCG inhibits the DNA methyltransferase of tumor cells through various signaling pathways, inhibit tumor cell proliferation, and promote apoptosis[17-19].
In this study, we aimed to investigate the inhibitory effects of EGCG on dimethylhydrazine (DMH)-induced CRC in rats. Through network pharmacology analysis, we also aimed to identify EGCG-associated targets and potential pathways in CRC.
Eighty six-week-old male Wistar rats weighing 180-200 g were purchased from Shanghai Laboratory Animal Center, CAS (SLACCAS) (production license number:SCXK 2007-0011). These rats were reared under SPF conditions (Jiangsu Pharmaceutical Biotechnology Co., Ltd.), fed cobalt-60 radiation-sterilized granulated feed, and adapted to laboratory conditions for two weeks. DMH [Aladdin Reagent(Shanghai) Co., Ltd.] was first prepared into a solution at a concentration of 4 mg/mL.EDTA (27 mg) was added to this solution, and the pH was adjusted to 6.5 with 0.1 mol NaOH. EGCG (Hangzhou Hetian Biotechnology Co., Ltd.) was dissolved in distilled water to a concentration of 50 mg/mL. For histological analysis, 4% (w/v)methylene blue solution was diluted in five volumes of distilled water (0.8% w/v). All experiments were conducted in accordance with the requirements of the Chinese Affidavit of the Approval of Animal Ethics and Welfare.
Eighty SPF Wistar rats were randomly divided into five groups as follows: Group 1,Vehicle (vehicle; saline, s.c., twice weekly) group (n= 10); Group 2, DMH (40 mg/kg,s.c., twice weekly) control group (n= 10); Group 3, DMH (40 mg/kg, s.c., twice weekly) + low-dose EGCG (50 mg/kg, p.o., once daily) group (n= 20); Group 4, DMH(40 mg/kg, s.c., twice weekly) + medium-dose EGCG (100 mg/kg, p.o., once daily)group (n= 20); Group 5, DMH (40 mg/kg, s.c., twice weekly) + high-dose EGCG (200 mg/kg, p.o., once daily) group (n= 20). DMH injection and EGCG administration both continued for eight weeks. Biweekly body weight measurements were performed to adjust the dosage of DMH. Five rats in each group were killed at week 12 and week 20.
Measurement of body weight and other physiological parameters: Each rat was weighed every two weeks to observe the change in body weight. Ascites formation was determined. The changes in anal and intestinal erosion, redness, swelling, and stool in rats were also recorded and observed. Tumor volume was measured using a Vernier caliper according to the following formula:
Gross tumor volume (mm3)[20]: V = (W2× L)/2 (V is tumor volume, W is tumor width, L is tumor length).
Tumor formation rate and tumor inhibition rate were also calculated using the following formulas:
Tumor formation rate (%) = (number of tumorigenic rats)/(total number of rats) ×100%.
Tumor inhibition rate (%)[21]= (C-T)/C × 100% (C is average weight of control group tumor, and T is average weight of experimental group tumor).Counting of aberrant crypt foci: The intestinal tissues of rats were isolated and washed with physiological saline. The whole intestine was opened up longitudinally,spread out flat, fixed in 10% neutral buffered formalin solution for 24 h, and stained with 0.2% methylene blue for 3-5 min. The intestinal tract was divided into small intestinal and colorectal parts. The number and size of aberrant crypt foci (ACF) were observed under an optical microscope at 40-200 × magnification.
Histopathology: The intestinal tissues were fixed in 95% ethanol for 24 h, dehydrated with a gradient ethanol series, and embedded in paraffin in blocks. The blocks were sliced into 3 μm thickness using a microtome, dewaxed with xylene, hydrated with a gradient ethanol series, and stained with hematoxylin and eosin. The histopathological changes in small intestinal and colorectal parts were observed by a double blinded histopathologist under an optical microscope.
The data of Canonical SMILES on EGCG was obtained from the PubChem database(https://pubchem.ncbi.nlm.nih.gov/). We finally found totally 87 targets of EGCG through Swiss Target Prediction Database (http://www.swisstargetprediction.ch/).
About 7787 targets of CRC were acquired from the Online Mendelian Inheritance in Man Database (http://omim.org/) and GeneCard (https://www.genecards.org/).
Using the collective targets of EGCG and CRC, a protein-protein interaction (PPI)regulatory network was established through Cytoscape (http://www.cytoscape.org/)and String (https://string-db.org/).
The expression profile of CRC was retrieved from The Cancer Genome Atlas (TCGA)portal (http://tcgadata.nci.nih.gov), and the dysregulated genes were analyzed using the edger R package.
To clarify the function and biological pathways of differentially expressed methylation loci from our data, Kyoto Encyclopedia of Genes and Genomes (KEGG,http://www.kegg.jp/) and the Gene ontology (GO, http://www.geneontology.org)pathway enrichment analyses were performed using the scripts in R.P< 0.05 in GO terms and KEGG was considered significantly enriched by genes.
All data were analyzed by theχ2test or Fisher's test using SPSS 20.0 (IBM, Armonk,NY, United States), and are represented as the mean ± SD. Multiple groups were compared by one-way ANOVA andP< 0.05 was considered to indicate a significant difference.
The average weight of the rats in all five groups was at the same baseline. Biweekly body weight measurements showed that rats in the DMH group (Group 2)significantly gained weight compared with rats in the DMH + high-dose EGCG group(Group 5) especially at week 6 (P< 0.01) and week 12 (P< 0.05). Ascites formation in rats in the DMH group (Group 2) at week 6 was significantly (P< 0.01) higher than that in the DMH + high-dose EGCG group (Group 5) (Table 1; Figure 1).
One rat in the DMH group (Group 2) died from an ulcer at the subcutaneous injection site; two rats died in the DMH + low-dose EGCG group (Group 3) from unexplained weight loss; and one rat in the DMH + medium-dose EGCG group died from a fight with other rats (Group 4). All dead rats pathologically manifested intestinal flatulence and small intestine dilatation, accompanied by ascites formation,ileal hemorrhage, necrosis, or colon darkening.
The stools of rats in the vehicle group (Group 1) were normal. Those of rats in the other groups were thin or loose. Four rats in the DMH group (Group 2) presented moderate anal erosion with ascites or bloody ascites, and abundant mesenteric nodules, intestinal wall stiffness, and fold disorder or enlargement in the intestinal mucosa were observed (Figure 2).
Tumor volumes in Groups 4 and 5 were significantly (P< 0.05) smaller than thosein the DMH group (Group 2), as shown in Table 2. The tumor volumes in the DMH +medium-dose EGCG group (Group 4) and DMH + high-dose EGCG group (Group 5)were significantly (P< 0.05) smaller than that in the DMH group (Group 2).
Table 1 Average weight (g) of each group at different time points (mean ± SD)
Observation under a light microscope showed that the intestinal mucosa of the vehicle group (Group 1) was intact, with a few orderly crypts, and no malignant tumor was observed. However, varying degrees of abundant ACF appeared in Groups 2-5, presenting as small abscesses with mucosal vascular congestion and inflammatory cell infiltration (Figure 3). The ACF count for the DMH + high-dose EGCG group (Group 5) was significantly (P< 0.01) lower than that in the DMH group(Group 2) (Table 2).
Histopathological analysis showed that 811 intestinal tumors of different sizes were found in the intestines of all rats, among which 178 were cancerous tumors (21.9%)and 633 were non-cancerous tumors (78.1%) determined through pathological examination (Figure 4).
At week 12, the total tumor formation rates in the DMH group (Group 2), DMH +low-dose EGCG group (Group 3), DMH + medium-dose EGCG group (Group 4), and DMH + high-dose EGCG group (Group 5) were 50.13%, 44.27%, 40.18%, and 32.75%,respectively. The total tumor formation rates in Group 4 and Group 5 were significantly (P< 0.01) lower than that in Group 2. It was lower in Group 5 than in Group 3 (P< 0.05) as well. The number of cancerous tumors was lower in Group 5 than in Group 2 (P< 0.05). Tumor inhibition rates in Group 3, Group 4, and Group 5 were 10.26%, 12.19%, and 13.43%, respectively. There was no significant difference in the tumor inhibition rate between each of these three groups (P> 0.05).
At week 20, the tumor formation rates in Groups 2-5 were 77.84%, 63.13%, 51.96%,and 40.12%, respectively. The total tumor formation rates in Group 3-5 were significantly (P< 0.01) reduced compared with that in Group 2. The incidence of cancerous and non-cancerous tumors in Group 3-5 was lower than that in Group 2 (P< 0.01). Group 5 had a statistically significant lower tumor formation rate than Group 3 (P< 0.05). Tumor inhibition rates in Groups 3-5 were 19.33%, 23.58%, and 34.53%,respectively. The tumor inhibition rate in Group 5 was significantly higher than that in Group 3 (P< 0.05) (Table 3, Figure 5 and Figure 6).
The 2-D and 3-D structures are shown in Figure 7A and 7B, respectively. Based on the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform databases, 78 target genes of EGCG were constructed into a protein-protein regulatory network composed of 88 connected nodes (Figure 7C).
Online Mendelian Inheritance in Man and GeneCard databases were used to identify 7787 targets of CRC for constructing a protein-protein regulation network (Figure 8A). Thus, 78 collective targets of EGCG and CRC were obtained (Figure 8B) and used to construct the protein-protein regulation network (Figure 9A). Then, the 78 target genes were incorporated into a PPI network using String (Figure 9B); the top 30 genes are shown in Figure 9C.
To narrow down the targets of EGCG in CRC, we then analyzed the dysregulated CRC genes in TCGA (Figure 10A) and found that 28 dysregulated genes in CRC were involved in the collective targets of EGCG and CRC (Figure 10B).
Figure 1 Change curves of body weights of rats in different groups. DMH: Dimethylhydrazine; EGCG: Epigallocatechin gallate.
KEGG and GO analyses: KEGG analysis showed that the dysregulated genes were enriched in hsa05210 (CRC), hsa04115 (p53 signaling pathway), and hsa04151 (PI3KAkt signaling pathway) (Figure 10C). GO analysis showed that the dysregulated genes were enriched in GO:0043124 (negative regulation of I-kappaB kinase/NFkappaB signaling pathway), GO:2001244 (positive regulation of intrinsic apoptotic signaling pathway), and GO:0043409 [negative regulation of mitogen-activated protein kinase (MAPK) cascade] (Figure 10D).
The need for alternative treatments for CRC is mounting. Because green tea, and especially, EGCG in green tea have shown promise in cancer treatment, we sought to evaluate the anticancer mechanisms of EGCG in experimentally induced CRC in rats.EGCG was effective in reducing tumor induction, tumor volume, ascites formation,and ACF counts at all three doses, with the high dose (200 mg/kg) being the most effective.
In our study, we observed that except for the vehicle group rats, the stools of rats in other groups showed different degrees of abnormalities. An abnormal stool sample was characterized by thin consistency, shapelessness, looseness, and having one or more concave holes. Such abnormal stools were noted mostly in the DMH group, and the abnormality of stools was least obvious in the EGCG high-dose group.
In the DMH group, four rats had moderate erosion in the anus. Anatomical analysis revealed varying degrees of ascites or bloody ascites, a large number of tumors of varying sizes in the mesentery, and disordered and thickened intestinal folds in all groups except the control group. The average body weight for all five groups of rats was the same at baseline. Since week 6, the body weight was significantly reduced in the EGCG high-dose group (Group 5) compared with that in the DMH group (P< 0.01). At week 12, the body weight in Group 2 was higher than that in the other groups. At 8 wk, Groups 3 and 4 presented a sharp rise in body weight, and at week 20, body weights in both these groups were higher than that in Group 2. At week 12, the ascites production and tumor formation in the DMH group were greater than those in the other groups, and at this time, the rats in the DMH group were in the compensatory stage, and there was no significant loss of weight caused by cachexia. Anatomical analysis revealed that the DMH group had more ascites than the other groups at week 12. After that, rats in the DMH group gradually developed cachexia, had poor nutritional status, and showed fatigue. Weight loss due to cachexia was more prominent than overall weight gained due to increase in tumor mass; therefore, weight loss in the DMH group was significant. Groups 3 and 4 showed lower ascites and tumor formation and better health status than Group 2,indicating the anticancer effect of EGCG. In Group 5, the tumor suppression effect of EGCG was significant, and there was a reduction in tumor and ascites formation,resulting in lower body weights in this group.
Figure 2 Formation of tumor in the intestine of different groups. A: Dimethylhydrazine (DMH) group; B: DMH +low-dose epigallocatechin gallate (EGCG) group; C: DMH + medium-dose EGCG group; D: DMH + high-dose EGCG group; E: The whole intestine with tumors; F: Intestinal inflation and obstruction with tumors.
Morphological observation of ACF revealed that the intestinal epithelial cells in the DMH group were larger, more in number, and presented with crypts that changed from round to elliptical shape, increasing the crypt’s peripheral space. The size of the nucleus of the mucosal epithelial cells around ACF increased, the proportion of nucleoplasm increased, the number of goblet cells decreased, and the population of cells in the mitotic phase increased. These changes are the cytological features of precancerous lesions. The process of intestinal carcinogenesis involves conversion of normal glandular epithelium to ACF, adenoma, and finally to tumor. All three stages are characterized by abnormal cell proliferation, inflammatory cell infiltration, and gene mutations[22,23]. It is currently recognized that the pathogenesis of CRC follows the pathway of ACF → adenoma → CRC[24]. Therefore, observing the cell morphology of ACF is crucial for judging whether it is a precancerous lesion. This process is consistent with the occurrence of human CRC.
To determine if the anticancer effects of EGCG are dose-dependent, we administered EGCG at three doses: 50 mg/kg, 100 mg/kg, and 200 mg/kg. From the data, we can see that with the extension of EGCG treatment duration and the increase in dose, the effect of EGCG in inhibiting tumor formation increased. Surprisingly,only the minimum dose of EGCG of 50 mg/kg had a significant inhibitory effect on tumor formation.
We used a network pharmacology approach to identify 78 collective targets of EGCG and CRC. After constructing the PPI network, we found that MYC, epidermal growth factor receptor (EGFR), vascular endothelial growth factor A, CASP3, IL6, and others had a close association with CRC. Overexpression of EIF5A2 promotes invasiveness of CRC by upregulating MTA1 through C-myc to induce epithelialmesenchymal transition[25]. Besides, EGFR was reported to mediate the re-activation of MAPK signaling, which increases the sensitivity ofBRAFmutant CRC to vemurafenib[26]. Vascular endothelial growth factor A, CASP3, and IL6 are closely related to CRC[27-29]. Masahito’s experiments investigated the effects of EGCG on VEGF/vascular EGFR (VEGFR) axis activity and the expression of hypoxia-induciblefactor (HIF)-1α, suggesting that the factor promotes angiogenesis by increasing VEGF levels in CRC. EGCG treatment decreased the mRNA expression ofHIF-1α, VEGF,EGF, insulin-like growth factor (IGF)-1,IGF-2, andHeregulinwithin 3 h in SW837 CRC cells, whereas it decreased the expression of VEGFR-2, p-VEGFR-2, p-ERK, p-IGF-1 receptor, and p-Akt proteins within 6 h[30]. This suggests that EGCG can inhibit HIF-1α and several major growth factors related to the VEGF/VEGFR axis to arrest the growth of CRC cells. Therefore, it is necessary to find the target genes of EGCG and CRC intersection through network pharmacology. In the future, we will conduct cellular molecular research on relevant target genes and target pathways through cell experiments, so as to further explore the molecular mechanisms of EGCG in the treatment of CRC.
Table 2 Aberrant crypt foci in different groups
In conclusion, we have successfully induced CRC in rats using DMH and found that EGCG at low and high doses can inhibit tumor formation. Moreover, the effect of EGCG positively correlates with dose and treatment duration. This reveals the clinical treatment potential of EGCG for CRC. Through bioinformatics analysis, we found 78 target genes that were common to EGCG and CRC. Through KEGG and GO analyses,we found that EGCG could regulate the signaling pathways of CRC, p53, and PI3KAkt and I-kappaB kinase/NF-kappaB signal pathways, apoptosis signal pathways,and MAPK cascades, respectively. Our findings reveal the potential mechanisms of action and pathways of EGCG in CRC. Furtherin vitrostudies using cell lines would be beneficial to verify the target genes and pathways of EGCG in CRC. One limitation of our study is that our CRC induction was short because of the fact that the tumor formation rate was low. Future studies can be conducted using azoxymethane, which is more potent and stable than DMH and has ability to induce significant tumor formation within a short time.
Table 3 Number and diameter of tumor in each group at 12th and 20th wk
Figure 3 lntestinal crypt abscess formation, vascular congestion of the mucosa and submucosa, and infiltration of inflammatory cells. A: Hematoxylin-eosin staining, × 200; B: Hematoxylin-eosin staining, × 40.
Figure 4 Pathological changes in different groups. A: Dimethylhydrazine (DMH) + low-dose epigallocatechin gallate (EGCG) group showed moderate to severe dysplasia and local carcinogenesis (Hematoxylin-eosin staining, × 100); B: DMH + low-dose EGCG group showed moderate to severe dysplasia and local carcinogenesis (Hematoxylin-eosin staining, × 400); C: DMH + medium-dose EGCG group showed moderate to severe dysplasia (Hematoxylin-eosin staining, × 100);D: DMH + medium-dose EGCG group showed moderate to severe dysplasia (Hematoxylin-eosin staining, ×400); E: DMH + high-dose EGCG group showed chronic inflammation (Hematoxylin-eosin staining, × 100); F: DMH + high-dose EGCG group showed chronic inflammation (Hematoxylin-eosin staining, × 400); G: DMH group showed signet-ring cell carcinoma (Hematoxylin-eosin staining, × 100); H: DMH group showed signet-ring cell carcinoma (Hematoxylin-eosin staining, × 400).
Figure 5 Tumor formation in the whole intestine, small intestine, and colorectum A: Week 12; B: Week 20. DMH: Dimethylhydrazine; EGCG: Epigallocatechin gallate.
Figure 6 Tumor formation in total, non-cancerous, and cancerous tumors. A: Week 12; B: Week 20. aP < 0.05, bP < 0.01. DMH: Dimethylhydrazine; EGCG:Epigallocatechin gallate.
Figure 7 Protein-protein regulation network established using targets of epigallocatechin gallate. A: 2-D structure of epigallocatechin gallate (EGCG); B: 3-D structure of EGCG; C: 87 candidate targets of EGCG were identified to establish the protein-protein regulation network. The yellow nodes represent the protein targets of EGCG, and the red node stands for the EGCG (For interpretation of the references to color in this figure legend, please refer to the web version of this article).EGCG: Epigallocatechin gallate.
Figure 8 ldentified collective targets of epigallocatechin gallate and colorectal cancer. A: Protein-protein regulation network established using 287 targets of colorectal cancer (CRC). The pink nodes represent the protein targets, and the yellow node stands for the CRC; B: Identified 78 collective targets of epigallocatechin gallate and CRC. CRC: Colorectal cancer; EGCG: Epigallocatechin gallate.
Figure 9 Protein-protein interaction network of 78 collective targets of epigallocatechin gallate and colorectal cancer. A: The 78 collective targets of epigallocatechin gallate and colorectal cancer were used to establish a protein-protein regulation network. The pink nodes represent the protein targets, red node stands for colorectal cancer, and the green nod represents epigallocatechin gallate; B: The protein-protein interaction network of the 78 collective targets; C: The top 30 targets with more influence were showed. CRC: Colorectal cancer; EGCG: Epigallocatechin gallate.
Figure 10 Potential pathways that epigallocatechin gallate may affect. A: Heatmap of dysregulated genes in colorectal cancer expression profile from The Cancer Genome Atlas; B: Twenty-eight collective genes that were dysregulated among the 78 collective targets of epigallocatechin gallate and colorectal cancer; C:Pathways that the 28 genes were enriched by Kyoto Encyclopedia of Genes and Genomes analysis; D: Biological process that 28 genes were enriched by Gene ontology analysis. CRC: Colorectal cancer; EGCG: Epigallocatechin gallate.
“Tea polyphenols” is a general term for more than 30 polyphenols in tea, of which epigallocatechin gallate (EGCG) is the main compound. In recent years, a large number of studies have demonstrated that EGCG is a highly effective, biologically active, non-toxic,naturally extracted compound. Its antitumor effect and mechanism of action are constantly being clarified, and they have become a hot topic of research today. However, most of the existing studies only use cell experiments to explain its inhibitory effect.
The aim of our research was to explore whether EGCG can inhibit the occurrence of dimethylhydrazine (DMH)-induced rat intestinal cancerin vivoin rats. We combined bioinformatic data analysis technology to explain the common targets and significant differential genes of EGCG and colorectal cancer (CRC) for the first time and analyze the mechanism of EGCG in CRC at the molecular gene level.
To, using animal experiments and network pharmacological analysis, discover the inhibitory mechanism of EGCG on CRC cell proliferation and clarify its pharmacological targets.
DMH (40 mg/kg) was used to induce CRC in rats twice weekly for 8 wk. Then, different doses of EGCG were administered (50, 100, or 200 mg/kg for 8 wk). Rats were euthanized at week 12 and week 20 to observe aberrant crypt foci and tumor formation. We used bioinformatic analysis to predict the pathways altered by EGCG in CRC.
We discovered that different doses of EGCG had inhibitory effects on aberrant crypt foci and tumors in rats and, as the dose increased, EGCG was more effective in the suppression of tumor formation. Using bioinformatics, we identified 78 target genes that were common to EGCG and CRC and 28 of the 78 genes were identified to have significant differential expression. Using KEGG and GO analyses, we found that EGCG regulated CRC, p53, PI3K-Akt, I-kappaB kinase/NF-kappaB, apoptosis, and mitogen-activated protein kinase signaling pathways.
EGCG can inhibit DMH-induced CRC tumor formation in rats. Moreover, the effect of EGCG positively correlated with dose and treatment duration. This reveals the clinical treatment potential of EGCG for CRC. Our findings reveal the potential mechanisms of action and pathways of EGCG in CRC.
In the future, azoxymethane, which is more potent and stable than DMH, could be used as it has the ability to induce tumor formation in a short time. Additionally,in vitrostudies using cell lines would be beneficial to verify the target genes and pathways of EGCG in CRC.
World Journal of Gastroenterology2020年17期