2020-06-13T00:27:34Zhttps:/www.ncbi.nlm.nih.gov/pmc/oai/oai.cgi
oai:pubmedcentral.nih.gov:4217472004-06-13carcinogpmc-open
J Carcinog Journal of Carcinogenesis 1477-3163 BioMed Central London PMC421747 PMC421747 421747 15140256 1477-3163-3-8 15140256 10.1186/1477-3163-3-8 Research Time- and dose-dependent effects of curcumin on gene expression in human colon cancer cells van Erk Marjan J 1 3 marjan.vanerk@wur.nl Teuling Eva 1 e.teuling@erasmusmc.nl Staal Yvonne CM 1 Y.Staal@GRAT.unimaas.nl Huybers Sylvie 1 sylviehuybers@hotmail.com van Bladeren Peter J 3 4 peter.van-bladeren@rdls.nestle.com Aarts Jac MMJG 1 3 jac.aarts@wur.nl van Ommen Ben 2 vanommen@voeding.tno.nl Wageningen University, Division of Toxicology, Tuinlaan 5, 6703 HE Wageningen, the Netherlands TNO Nutrition and Food Research, PO Box 360, 3700 AJ Zeist, the Netherlands Wageningen University/TNO Centre for Food Toxicology, PO Box 8000, 6700 EA Wageningen, the Netherlands Nestlé Research Centre, PO Box 44, CH-1000 Lausanne 26, Switzerland 2004 12 5 2004 3 8 8 28 1 2004 12 5 2004 Copyright © 2004 van Erk et al; licensee BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL. 2004 van Erk et al; licensee BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL. Background

Curcumin is a spice and a coloring food compound with a promising role in colon cancer prevention. Curcumin protects against development of colon tumors in rats treated with a colon carcinogen, in colon cancer cells curcumin can inhibit cell proliferation and induce apoptosis, it is an anti-oxidant and it can act as an anti-inflammatory agent. The aim of this study was to elucidate mechanisms and effect of curcumin in colon cancer cells using gene expression profiling.

Methods

Gene expression changes in response to curcumin exposure were studied in two human colon cancer cell lines, using cDNA microarrays with four thousand human genes. HT29 cells were exposed to two different concentrations of curcumin and gene expression changes were followed in time (3, 6, 12, 24 and 48 hours). Gene expression changes after short-term exposure (3 or 6 hours) to curcumin were also studied in a second cell type, Caco-2 cells.

Results

Gene expression changes (>1.5-fold) were found at all time points. HT29 cells were more sensitive to curcumin than Caco-2 cells. Early response genes were involved in cell cycle, signal transduction, DNA repair, gene transcription, cell adhesion and xenobiotic metabolism. In HT29 cells curcumin modulated a number of cell cycle genes of which several have a role in transition through the G2/M phase. This corresponded to a cell cycle arrest in the G2/M phase as was observed by flow cytometry. Functional groups with a similar expression profile included genes involved in phase-II metabolism that were induced by curcumin after 12 and 24 hours. Expression of some cytochrome P450 genes was downregulated by curcumin in HT29 and Caco-2 cells. In addition, curcumin affected expression of metallothionein genes, tubulin genes, p53 and other genes involved in colon carcinogenesis.

Conclusions

This study has extended knowledge on pathways or processes already reported to be affected by curcumin (cell cycle arrest, phase-II genes). Moreover, potential new leads to genes and pathways that could play a role in colon cancer prevention by curcumin were identified.

Background

Curcumin (diferuloylmethane) is a spice and a coloring agent derived from the root of the plant Curcuma longa to which colon cancer-preventive properties have been attributed. It is present in curry and mustard, and it is used extensively in Asian countries, also in traditional medicine. The low incidence of colon cancer in Asian countries could be related to low meat intake, but also to the regular use of curcumin in the diet [1].

In rats and mice curcumin has a profound effect on colon carcinogenesis. In rats treated with colon carcinogen azoxymethane (AOM), for example, consuming a diets with 2000 ppm curcumin resulted in a significant reduction of the number of aberrant crypt foci (ACF) after 9 weeks [2] and after 52 weeks incidence and multiplicity of colon adenocarcinomas were significantly reduced [3]. Similarly, incidence and multiplicity of adenomas was decreased in AOM-treated rats fed a diet with 8 or 16 g/kg curcumin for 45 weeks [4]. Also in a mouse model with a mutation in the APC gene curcumin reduced the number of colon tumors [5] or the multiplicity of colon adenomas [6]. Moreover, curcumin was found to protect against development of colon cancer during both the initiation and the promotion stage in AOM-treated rats and AOM-treated mice [7,8]. In addition to the effect on colon cancer, curcumin showed anticancer effects in intestinal cancer, stomach cancer and hepatocellular carcinoma [9-11].

In vitro, curcumin caused a dose-dependent decrease in cell proliferation in colon cancer cells and the cells accumulated in the G2/M phase [12-14]. Apoptosis was increased in colon cancer cells in response to curcumin [13,14]. Also in other cell lines, e.g. breast cancer cells, prostate cancer cells and leukaemia cells, curcumin inhibited cell proliferation and induced apoptosis [15-19]. In addition to inhibition of cell proliferation and increased apoptosis, many mechanisms have been proposed to explain the anti-carcinogenic effect of curcumin, including its anti-inflammatory and antioxidant activity, induction of phase-II detoxification enzymes, inhibition of cyclooxygenase 2 (COX-2), effect on AP-1 and NFκB transcription factors, inhibition of matrix metalloproteinase (MMP), effect on protein kinases and more [20-22].

In rats absorption of curcumin from the intestine was reported to be about 60% [23]. Curcumin and metabolites formed in intestine and liver are mostly excreted in the faeces [24,25]. As the colon is exposed to both curcumin and its metabolites, it is a likely target for the anticarcinogenic activity of these compounds. Moreover, the fact that humans were able to consume up to 8 grams of curcumin per day without toxic effects [26] makes curcumin a very interesting chemopreventive agent.

New techniques, like multiple gene expression analysis using microarrays, allow for a more comprehensive study of the effects and mechanisms of food components. By using cDNA microarrays the expression of thousands of genes can be studied in one experiment. The power of these techniques lies not only in the fact that many genes can be studied in one experiment, but also in the possibility to identify leads to pathways and mechanisms, as opposed to the extrapolation of assessments of a few genes. In addition, multiple gene expressions can be applied as 'fingerprint' biomarkers.

The aim of our study was to use new and genome-wide information on the gene expression profile induced by curcumin in colon cancer cells to elucidate mechanisms involved in the cancer-preventive action of curcumin. Therefore, HT29 colon cancer cells were exposed to two concentrations of curcumin and gene expression changes were measured at five exposure time points ranging from 3 hours to 48 hours. In addition, the response to short-term exposure to curcumin (3 and 6 hours) was also studied in Caco-2 cells. Furthermore, in HT29 cells changes in expression of cell cycle genes in response to curcumin were related to changes in cell cycle distribution. The time- and concentration-dependent changes in gene expression in HT29 and Caco-2 cells are reported and unique findings and observed similarities are discussed in relation to data from other (microarray) studies to gain more insight into mechanisms of cancer prevention by curcumin.

Methods Cell culture

HT29 cells (ATCC, Rockville, USA) were grown in McCoy's 5 A medium with L-glutamine (Gibco BRL Life Technologies) with 10% FCS. Caco-2 cells (ATCC, Rockville, USA) were grown in DMEM with 25 mM HEPES, without sodium pyruvate, with 4500 mg/l glucose, with pyridoxine (Gibco BRL Life Technologies), to which was added 10% FCS, 1% non-essential amino acids (Gibco BRL Life Technologies) and 2% penicillin/streptomycin (Gibco BRL Life Technologies). Cells were maintained at 37°C in a humidified atmosphere containing 5% CO2. Curcumin (Sigma, St. Louis, MO, USA) was dissolved in DMSO. At 70% confluence HT29 cells and Caco-2 cells were exposed to 30 μmol/L (referred to as low concentration) or 100 μmol/L curcumin (referred to as high concentration) for 3 or 6 hours. In addition, HT29 cells were exposed to 25 μmol/L (referred to as low concentration) or 100 μmol/L curcumin (referred to as high concentration) for 12, 24 or 48 hours. The final DMSO concentration in the medium was 0.1%.

Cell cycle analysis using flow cytometry

HT29 cells were plated at a density of 3 × 105 cells in 25 cm2 culture flasks. After 24 hours, cells were exposed to curcumin. After the exposure period cells were trypsinized and collected in the original medium (to include floating cells in the analysis). Cells were pelleted by centrifugation at 500 × g at 4°C and washed with PBS. Cells were resuspended in PBS, ice-cold ethanol (75%) was added to the cells while vortexing and cells were incubated on ice for 1 hour to fixate the cells. Cells were washed with PBS and finally PBS with propidium iodide (Sigma, St. Louis, MO, USA; 50 μg/ml) and Rnase A (Qiagen, Hilden, Germany; 0.1 mg/ml) was added. After incubation for 30 min in the dark, cells were analysed using an Epics XL-MCL flow cytometer (Beckman Coulter). P-values for difference between cell cycle distribution in treated cells and in untreated cells were calculated using a Student's t-test. Also, after trypsinization cells were counted using a Bürker-Turk counting chamber.

RNA isolation

After exposure total RNA was isolated from the cells using Trizol (Life Technologies S.A., Merelbeke, Belgium) according to the manufacturer's protocol. RNA clean-up and Dnase digestion was performed using the RNeasy mini kit (Qiagen, Hilden, Germany). RNA was checked for purity and stability by gel electrophoresis and UV spectrometry. Absorption at 260 and 280 nm was measured and RNA quantity and A260/A280 ratio were calculated. Only RNA samples with A260/A280 ratio > 1.6 were used in further experiments. Similar to other studies [27,28], RNA from two or three exposure experiments was pooled before labelling and hybridisation to reduce possible bias from single exposure.

Quantitative real-time polymerase chain reaction

2 μg of total RNA was reverse transcribed into cDNA using 250 ng random hexamer primers and 200 units M-MLV reverse transcriptase (Invitrogen Life Technologies, Breda, the Netherlands) in a final volume of 20 μl. The same batch of cDNA was used for all real-time PCR experiments. Real-time PCR was performed using an iCycler PCR machine (Biorad, Veenendaal, the Netherlands) and the QuantiTect SYBR Green PCR kit (Qiagen, Hilden, Germany). For all types of amplicon reactions were performed in a final volume of 25 μl containing 400 nM of each primer, 5 μl of diluted cDNA preparation (equivalent to 50 or 5 ng of starting RNA) and 1× QuantiTect SYBR Green Master Mix (Qiagen, Hilden, Germany). GAPDH and beta actin expression was measured to correct for differences in efficiency during reverse transcription. The primer pairs used to amplify beta-actin, GAPDH, AKR1C1 and EGR1 were as follows: beta-actin upstream: 5'-CAC CCC GTG CTG CTG AC-3', downstream: 5'-CCA GAG GCG TAC AGG GAT AG-3'; GAPDH upstream: 5'-TGC ACC ACC AAC TGC TTA GC-3', downstream: 5'-GGC ATG GAC TGT GGT CAT GAG-3'; AKR1C1 upstream: 5'-CGT GGG AGG CCG TGG AGA AG-3', downstream: 5'-GCT GCC TGC GGT TGA AGT TGG-3'; EGR1 upstream: 5'-GTC CCC GCT GCA GAT CTC T-3', downstream: 5'-CTC CAG CTT AGG GTA GTT GTC CAT-3'.

Before real-time PCR analyses the specificity of the used primer pairs and the correct size of the obtained PCR products was checked by gel electrophoresis. At the end of all real-time PCR runs, a melt curve peak analysis was performed to ensure amplification of only the correct product.

Amplification efficiency for each gene was calculated using a dilution series from a mixture of all cDNAs, as recommended in Pfaffl et al. [29]. A threshold was chosen so that for all samples in the dilution series amplification was in the exponential phase and the correlation coefficient (as found for the linear relation between threshold cycle value and logarithm of starting quantity) was maximal. This threshold level was also applied to determine the threshold cycle value for all individual samples. Each sample was measured in duplo at two different starting concentrations of cDNA (equivalent to 50 or 5 ng of starting RNA).

cDNA microarray preparation

A set of 4069 sequence-verified human cDNA clones from the I.M.A.G.E. consortium was purchased (Research Genetics, U.S.A.) as PCR products. The cDNA was amplified by PCR with forward (5'-CTGCAAGGCGATTAAGTTGGGTAAC-3') and reverse (5'-GTGAGCGGATAACAATTTCACACAGGAAACAGC-3') primers. The primers contained a 5'-C6-aminolinker (Isogen Bioscience, Maarsen, The Netherlands) to facilitate cross-linking to the aldehyde-coated glass microscope slides. PCR products were checked by electrophoresis on a 1% agarose gel. PCR products were purified by isopropanol precipitation and washing in 70% ethanol, and were dissolved in 3 × SSC. The clones were spotted on CSS-100 silylated aldehyde glass slides (TeleChem, Sunnyvale, CA, USA) in a controlled atmosphere. To reduce free aldehyde residues, slides were blocked with borohydride after spotting and drying. Slides were stored in the dark and dust-free until further use.

Labelling and hybridisation

Total RNA (25 μg) was labelled using CyScribe first-strand cDNA labelling kit (Amersham Biosciences, Freiburg, Germany). During reverse transcription of the RNA, Cy3- or Cy5-labelled dUTP was built into the cDNA. After incubating the reaction for 1.5 h at 42°C, RNA was hydrolysed by adding NaOH. After neutralisation, free nucleotides were removed from the solution using AutoSeq G50 columns (Amersham Biosciences, Freiburg, Germany).

Each sample was labelled twice, once with Cy3 and once with Cy5. cDNA from cells exposed to curcumin (treated sample) was hybridised to the microarray in competition with cDNA from cells exposed to DMSO only (untreated sample), while one of the samples was labelled with Cy3 and the other with Cy5. Hybridisations were repeated with a dye swap.

Before hybridisation, Cy3- and Cy5-labeled cDNAs were mixed and human cot-1 DNA (3 μg, Life Technologies S.A., Merelbeke, Belgium), yeast tRNA (100 μg, Life Technologies S.A., Merelbeke, Belgium) and poly(dA-dT) (20 μg, Amersham Biosciences, Freiburg, Germany) were added to avoid non-specific binding. The hybridisation mix was dissolved in 30 μl Easyhyb hybridisation buffer (Roche Diagnostics, Mannheim, Germany), denatured for 1.5 min at 100°C and incubated for 30 minutes at 42°C.

Before adding the hybridisation mix to the slides, slides were prehybridised in prehybridisation buffer (5 × SSC, 0.1 % SDS and 10 mg/ml bovine serum albumin) for 45 minutes at 42°C, washed in milliQ water, washed with isopropanol and dried.

After pipetting the hybridisation mix on the slides, the slides were covered with a plastic coverslip and hybridised overnight in a slide incubation chamber (Corning, Life Sciences, Schiphol, the Netherlands) submerged in a 42°C waterbath. After hybridisation, slides were washed by submersion and agitation in 0.5 × SSC with 0.2% SDS and in 0.5 × SSC. Then, slides were firmly shaken in 0.2 × SSC and put on a rotation plateau for 10 min. This step was repeated once and slides were dried quickly by centrifugation at 700 rpm.

Slides were scanned with a ScanArray Express confocal laser scanner (Perkin Elmer Life Sciences, USA) and Imagene 4.0 (Biodiscovery Inc., Los Angeles, USA) was used to extract data from the images, with automatic flagging of weak or negative signals and spots with non-homogeneous signal.

Data analysis

Data were imported into SAS Enterprise guide V2 (SAS Institute Inc., Cary, USA). Spots with a signal/background ratio less than 1.5 or spots that were flagged by the Imagene software were not included in the data analysis. For each spot, local background intensity was subtracted from mean signal intensity and the expression ratio was calculated by dividing background-corrected signal intensity of the treated sample by the background-corrected signal intensity of the untreated sample. Expression ratios were then log transformed (base 2), normalised per slide using an intensity-dependent method (Lowess) [30] and scaled. Data were transferred to Microsoft Excel 97 (Microsoft Corporation, USA). Expression ratios of duplo (dye swap) arrays were combined and an average expression ratio was calculated, provided that an expression ratio was present for both arrays. In addition, genes that showed specific dye bias at all timepoint were excluded from further analysis and when looking for early or late response genes only genes that did not show dye bias were included.

Since it is impossible to report results for all genes on the array for all time points, we chose to look at genes with a relatively large change in expression over time. Among these genes we identified early response genes (which were more than 1.5 fold induced or repressed after 3 or 6 hours of exposure) and genes that respond at later time points. In addition, we identified genes with similar responses to curcumin within a pathway and we looked at pathways known to be important in colon carcinogenesis, like cell cycle control and apoptosis. Clustering methods like K-means clustering were used to identify genes with a similar expression profile over the different time points.

Principal Component Analysis

To analyze and visualize the results of the microarray experiments Principal Component Analysis (PCA) was used, which is a well-known pattern recognition method in the field of multivariate data analysis. Data analysis was performed using the Matlab software (The MathWorks, Inc., 1984–2001) version 6.1.0.450 (R12.1). For PCA analysis the procedure from the PLS Toolbox was used (Version 2.0.1b 1999, Eigenvector Research, Inc., 1995–1999). Average expression ratios were used in the PCA analysis. Only genes with less than three missing values among the different time points were included. Before PCA analysis, data sets were mean-centered.

Results

Curcumin exposure resulted in gene expression changes in HT29 cells at all time points. However, after exposure to curcumin for 48 hours only 20–30 genes were up- or downregulated more than 1.5 fold, while after for example 3 or 12 hours about 130 genes were up- or downregulated more than 1.5 fold. The total set of gene expression data was submitted to principal component analysis (PCA), identifying the two major components within the total variation between samples, which are then visualized in a two-dimensional plot in which the expression profiles at each of the time points are projected. This PCA plot (Figure 1) nicely visualizes the time-dependent changes in the gene expression patterns. Interestingly, the time-dependent shift of the overall gene expression patterns is similar for the low and the high curcumin concentration. The points representing the gene expression patterns after exposure to curcumin for 48 hours are plotted quite distantly from the gene expression patterns at the other time points, also indicating that the expression patterns of the cells exposed to curcumin for 48 hours are different from the gene expression patterns at the other time points.

Principal Component Analysis (PCA). The gene expression pattern at each time point is visualized as a dot in this two-dimensional graph. The axes show the scores of the gene expression profiles of HT29 cells at the specific time points after exposure to curcumin in principal component (PC) 1 and 2 (which explain the largest part of the variance). The variance explained by the PCs is indicated between parentheses. Open circles and dotted arrows indicate exposure to the low curcumin concentration, filled circles and solid arrows indicate exposure to the high curcumin concentration.

Early response genes were defined as genes that were differentially expressed after exposure to curcumin for 3 or 6 hours. A selection of the early response genes in HT29 cells is listed in table 1 (low concentration) and table 2 (high concentration). Early response genes were involved in several processes, including cell cycle control, signal transduction, DNA repair, transcription regulation, cell adhesion and xenobiotic metabolism. Both curcumin concentrations caused an increase in expression of genes involved in DNA repair, e.g. MLH1, MSH3 and ERCC2 (Tables 1 and 2). Upregulated signal transduction genes included STAT3 and STAT5b (Table 1) and some genes of the MAPK signal transduction pathway (MAP3K10, MAP4K2; table 1 and 2). Some other MAPK signal transduction genes were downregulated by the low curcumin concentration (Table 1). Expression of a group of genes involved in cell adhesion and protein binding was induced by short-term exposure to curcumin, including annexin (Table 1) and integrin genes (Table 2). Several genes involved in xenobiotic metabolism were downregulated after short-term exposure to the high concentration of curcumin, namely GSTT2, GSTM4, CYP1B1 (Table 2). Expression of GCLC, involved in glutathione synthesis, was upregulated after 6 hours (Table 2).

Early response genes in HT29 cells after exposure to low concentration of curcumin (30 μmol/L).

Fold change (treated/control)

Accession number Gene Name Gene Symbol 3 h 6 h Involved in
R10662 mutL homolog 1, colon cancer, nonpolyposis type 2 (E. coli) MLH1 1.90 1.04 DNA repair
AA421716 mutS homolog 3 (E. coli) MSH3 1.11 1.53 DNA repair
R54492 excision repair cross-complementing rodent repair deficiency, complementation group 2 ERCC2 1.63 1.20 DNA repair
AA399410 signal transducer and activator of transcription 3 (acute-phase response factor) STAT3 1.56 1.19 Signal transduction/transcription regulation
AA280647 signal transducer and activator of transcription 5B STAT5B 1.67 -1.04 Signal transduction/transcription regulation
AA434420 protein tyrosine phosphatase, non-receptor type 9 PTPN9 1.74 1.64 Signal transduction
H01340 mitogen-activated protein kinase kinase kinase 10 MAP3K10 1.83 1.42 Signal transduction
AA425826 mitogen-activated protein kinase kinase 2 MAP2K2 -1.72 -1.03 Signal transduction
T94169 mitogen-activated protein kinase 8 MAPK8 -1.55 1.16 Signal transduction
R82176 MAD, mothers against decapentaplegic homolog 7 MADH7 -1.77 NA1 Signal transduction
N71159 metastasis associated 1 MTA1 1.68 NA Transcription regulation
AA448256 metal-regulatory transcription factor 1 MTF1 -1.34 -1.66 Transcription regulation
AA465236 forkhead box O3A FOXO3A -1.43 -2.01 Transcription regulation
AA775415 SMT3 suppressor of mif two 3 homolog 2 (yeast) SMT3H2 -1.77 -1.57 Nuclear transport
T54121 cyclin E1 CCNE1 1.63 NA Cell cycle
N74285 CDC5 cell division cycle 5-like CDC5L 1.51 1.05 Cell cycle
AA789328 cyclin-dependent kinase (CDC2-like) 10 CDK10 1.51 NA Cell cycle
AA877595 cyclin-dependent kinase inhibitor 2A (melanoma, p16, inhibits CDK4) CDKN2A 1.49 1.86 Cell cycle
AA488324 BUB1 budding uninhibited by benzimidazoles 1 homolog beta (yeast) BUB1B -1.95 NA Cell cycle
R19158 serine/threonine kinase 6 STK6 -1.25 -1.80 Cell cycle
AA071486 serine/threonine kinase 12 STK12 -1.50 -1.83 Cell cycle
AA082943 cyclin G1 CCNG1 -1.54 1.24 Cell cycle
W51794 matrix metalloproteinase 3 (stromelysin 1, progelatinase) MMP3 -1.12 1.52 Breakdown of extracellular matrix
N33214 matrix metalloproteinase 14 (membrane-inserted) MMP14 -1.53 NA Breakdown of extracellular matrix
AA406571 carcinoembryonic antigen-related cell adhesion molecule 1 (biliary glycoprotein) CEACAM1 1.76 1.18 Cell adhesion
AA464982 annexin A11 ANXA11 1.75 1.13 Protein binding
AA419015 annexin A4 ANXA4 NA 1.79 Protein binding
AA856874 furin (paired basic amino acid cleaving enzyme) FURIN 1.87 1.48 Proprotein convertase
R26186 protein phosphatase 1, catalytic subunit, beta isoform PPP1CB 1.64 1.65 Phosphatase activity
N28497 protein phosphatase 2 (formerly 2A), regulatory subunit A (PR 65), beta isoform PPP2R1B 1.91 NA Phosphatase activity
AA877845 LIM domain kinase 2 LIMK2 1.88 1.56 Protein-protein interactions
AA634028 major histocompatibility complex, class II, DP alpha 1 HLA-DPA1 -2.06 NA Immune-related
T63324 major histocompatibility complex, class II, DQ alpha 1 HLA-DQA1 1.60 1.17 Immune-related
H95960 secreted protein, acidic, cysteine-rich (osteonectin) SPARC -1.13 -1.79 Collagen-binding Calcium-binding
AA456291 Developmentally regulated GTP binding protein 2 DRG2 1.67 1.47 GTP-binding protein
AA430178 RAN binding protein 2-like 1 RANBP2L1 -1.85 -1.11 GTP-binding protein

1 NA: not available

Early response genes in HT29 cells after exposure to high concentration of curcumin (100 μmol/L).

Fold change (treated/control)

Accession number Gene Name Gene Symbol 3 h 6 h Involved in
AA421716 mutS homolog 3 (E. coli) MSH3 1.52 1.47 DNA repair
R91570 signal transducer and activator of transcription 4 STAT4 1.17 -1.57 Signal transduction
H01340 mitogen-activated protein kinase kinase kinase 10 MAP3K10 1.81 1.33 Signal transduction
R50953 mitogen-activated protein kinase kinase kinase kinase 2 MAP4K2 1.65 NA1 Signal transduction
R38343 protein tyrosine phosphatase, receptor type, G PTPRG 1.56 1.05 Signal transduction
AA486533 early growth response 1 EGR1 8.66 1.29 Transcription regulation
AA600217 activating transcription factor 4 (tax-responsive enhancer element B67) ATF4 1.94 NA Transcription regulation
R83270 TGFB-induced factor (TALE family homeobox) TGIF -1.97 -1.06 Transcription regulation
AA733038 polymerase (RNA) I polypeptide C, 30 kDa POLR1C -1.66 -1.09 Transcription regulation
AA857163 amphiregulin (schwannoma-derived growth factor) AREG 1.98 1.23 Growth factor (cell cycle)
R45059 vascular endothelial growth factor VEGF 1.78 NA Growth factor (cell cycle)
R54846 fibroblast growth factor receptor 1 FGFR1 1.78 NA Cell growth (signal transduction)
AA877595 cyclin-dependent kinase inhibitor 2A (melanoma, p16, inhibits CDK4) CDKN2A 1.72 NA Cell cycle
AA450265 proliferating cell nuclear antigen PCNA -1.59 -1.09 Cell cycle
N74956 insulin-like growth factor binding protein 7 IGFBP7 -1.55 -1.39 Cell proliferation
AA180742 tubulin, alpha 1 (testis specific) TUBA1 -1.53 -1.09 Cytoskeleton
AA480995 methylene tetrahydrofolate dehydrogenase (NAD+ dependent), methenyltetrahydrofolate cyclohydrolase MTHFD2 1.79 1.56 Single carbon unit transfer
H23187 carbonic anhydrase II CA2 -1.09 -1.69 Single carbon unit transfer
AA037229 integrin, beta 3 (platelet glycoprotein IIIa, antigen CD61) ITGB3 -1.21 -1.50 Cell adhesion
AA485668 Integrin, beta 4 ITGB4 1.82 1.21 Cell adhesion
W56754 Integrin, beta 8 ITGB8 1.51 -1.08 Cell adhesion
R48796 integrin, alpha L (antigen CD11A (p180), lymphocyte function-associated antigen ITGAL 1.50 1.19 Cell adhesion
AA159577 mucin 5, subtype B, tracheobronchial MUC5B 1.84 1.27 Extracellular matrix formation
AA490208 glutathione S-transferase theta 2 GSTT2 -1.56 -1.48 Xenobiotic metabolism
AA486570 glutathione S-transferase M4 GSTM4 -1.13 -1.54 Xenobiotic metabolism
AA448157 cytochrome P450, family 1, subfamily B, polypeptide 1 CYP1B1 -1.74 -1.33 Xenobiotic metabolism
H56069 glutamate-cysteine ligase, catalytic subunit GCLC -1.16 1.52 Glutathione synthesis
AA759046 dual specificity phosphatase 2 DUSP2 2.31 1.24 Inactivation MAPK family kinases
AA187349 ferredoxin 1 FDX1 -1.57 NA Electron transfer, metabolism
AA450205 translocation protein 1 TLOC1 -1.54 -1.54 Protein transport
H13691 major histocompatibility complex, class II, DM beta HLA-DMB 1.51 -1.01 Immune-related
T70057 immunoglobulin J polypeptide IGJ -1.91 1.02 Immune-related
AA496359 immediate early protein ETR101 1.54 1.10 Unknown
AA425687 DEAD/H (Asp-Glu-Ala-Asp/His) box polypeptide 1 DDX1 -1.56 NA Unknown

1 NA: not available

Also, genes involved in transcription regulation were induced or repressed by short-term exposure to curcumin. Among these were transcription factors such as activating transcription factor 4 (ATF4) and early growth response 1 (EGR1) (Table 2). One of the target genes regulated by ATF4 is asparagine synthetase (ASNS) [31]. Figure 2 shows the expression profile of ATF4 and ASNS in response to curcumin. Induction of ATF4 expression at the early time points is followed by induction of expression of ASNS at the same and later time points. EGR1, a transcription factor involved in cell growth regulation and tumor suppression [32], was the most upregulated early response gene. In contrast to this strong upregulation after exposure to the high concentration of curcumin, this gene was not induced by the low curcumin concentration (Figure 3A). Expression of EGR1 was also measured by real-time RT-PCR. The induction factor for EGR1 after 3 hours of exposure to the high curcumin concentration was even higher when measured with real-time PCR than with the cDNA microarray (23.0 ± 2.73 when normalised to beta-actin and 30.5 ± 2.96 when normalised to GAPDH versus 8.7 as measured on cDNA microarray). Similarly, downregulation of EGR1 after 12 hours (2-fold) and upregulation after 24 hours of exposure (1.4-fold) were confirmed by real-time RT-PCR analysis of the effects on EGR1 expression (data not shown).

Expression profile of activating transcription factor 4 (ATF4, black bars) and asparagine synthetase (ASNS, grey bars) in response to exposure to curcumin. A: low concentration, B: high concentration.

Expression profile of genes in response to curcumin; grey bars indicated the low concentration and black bars indicate the high concentration. EGR1: early growth response 1; PLK: polo-like kinase; TP53: tumor protein p53; CA2: carbonic anhydrase 2; MTHFD2: methylene tetrahydrofolate dehydrogenase/cyclohydrolase; AKT1: protein kinase B/Akt; PLAU: urokinase-type plasminogen activator; PLAUR: urokinase-type plasminogen activator receptor; TM4SF1: transmembrane 4 superfamily member 1; TM4SF4: transmembrane 4 superfamily member 4.

In addition to the functional groups mentioned above several genes involved in the cell cycle or cell growth were among the early response genes. For example several growth factors (AREG, VEGF, FGFR1) were upregulated three hours after exposure to the high concentration of curcumin (Table 2). At the same time point, expression of cyclin-dependent kinase inhibitor p16INK4 (CDKN2A) was upregulated (Table 2). Expression of cell growth related genes PCNA and IGFBP7 was downregulated by the high concentration of curcumin (Table 2). Also after exposure to the low curcumin concentration some cell cycle genes were downregulated, like BUB1B (a mitotic checkpoint gene), STK6, STK12 and cyclin G1 (CCNG1). Upregulated cell cycle genes included CDK10, CDC5L and cyclin E1 (CCNE1), but also cell cycle inhibitor p16INK4 (CDKN2A) (Table 1). Also at later time points differential expression of genes involved in cell cycle or cell growth was found. One of the most strongly downregulated genes in HT29 cells was polo-like kinase (PLK), after exposure to the high curcumin concentration for 24 hours (fig. 3B). PLK is a cell cycle gene involved in spindle assembly. It is expressed at a higher level in colorectal cancer than in normal colon tissue [33]. Downregulation of PLK has been shown to inhibit cell growth in cancer cells [34]. Several histone genes (H3F3A, HIST1H4C) were downregulated by curcumin, especially at the 12 h and 24 h time points (data not shown). Histone deacetylase (HDAC1) was downregulated 1.7-fold by the high concentration of curcumin after 24 hours. Other genes involved in cell cycle control that were differentially expressed after 12 or 24 hour exposure to high concentration curcumin were retinoblastoma 1 (RB1), MAD2L1, BUB1, cyclin G1 (downregulated 1.6 to 2-fold). In response to exposure to the low concentration of curcumin for 12 hours expression of cyclin H was downregulated (1.5-fold) and expression of cyclin A2 was upregulated (1.7-fold). In contrast to cell cycle-related genes, only a few genes involved in apoptosis were differentially expressed in response to curcumin. The most striking effect was observed with programmed cell death 2 (PDCD2), which was downregulated 1.8-fold by the high concentration of curcumin after exposure for 12 hours.

When looking in detail into the gene expression profiles at the different time-points, several genes known to be involved in colon carcinogenesis were found that responded to curcumin exposure. Expression changes of these genes are shown in figure 3. Protein expression of urokinase-type plasminogen activator (PLAU) was reported to be higher in colon tumours than in normal colon tissue [35]. In our study, curcumin downregulated expression of PLAU in HT29 cells (Figure 3). In addition, expression of urokinase-type plasminogen activator receptor (PLAUR) was upregulated by curcumin (Figure 3). Recently, it was shown that members of the transmembrane 4 superfamily could play a role in colon cancer [36]. Two members of this family, TM4SF1 and TM4SF4, showed differential responses to exposure to curcumin (Figure 3 and 3). Expression of carbonic anhydrase 2 (CA2) was found to be downregulated in colon tumors [37,38]. In this study, expression of CA2 was downregulated in response to exposure to curcumin (Figure 3). Expression of protein kinase B (AKT1) was upregulated by the low curcumin concentration after 3 hours and by the high curcumin concentration after 24 hours (Figure 3) Protein kinase B/Akt is an important part of signal transduction pathways regulating for example apoptosis. Curcumin downregulated expression of tumor protein p53 (TP53) (Figure 3), a transcription factor that has a role in cell cycle control [39]. Expression of methylene tetrahydrofolate dehydrogenase/cyclohydrolase (MTHFD2) was upregulated at early time points and downregulated at the later time points. This gene has a role in the maintenance of the single carbon metabolic pool carried by folate. Folate could be involved in colorectal carcinogenesis through DNA methylation and polymorphisms in genes involved in folate metabolism [40].

Clustering methods were used to identify groups of genes with similar expression profile across the different time points. Several functional groups of genes with a similar expression profile were found. Figure 4 shows the response of the cluster of genes involved in phase-II metabolism (biotransformation of electrophilic mutagens and other toxic compounds) to exposure to curcumin at the different time points. Curcumin induced expression of these genes, with a maximum induction after exposure for 12 or 24 hours. After 48 hours no differential expression is seen. The most upregulated gene is aldo-keto reductase family 1 member C1 (AKR1C1). The expression changes of AKR1C1 in response to curcumin were confirmed with real-time RT-PCR (Figure 5). This gene, also known as dihydrodiol dehydrogenase, was recently shown to be underexpressed in 50% of gastric cancers [41]. Previously it was also shown that gene expression of AKR1C1 could be upregulated in HT29 cells by phenolic antioxidants [42]. Isothiocyanates induced protein expression of AKR1C1, NAD(P)H:quinone oxidoreductase 1 (NQO1) and heavy subunit of glutamate-cysteine ligase (GCLC) in a colon cancer cell line [43]. These genes were also upregulated by curcumin in this study. In addition, NQO2 and epoxide hydrolase 1 (EPHX1) were also upregulated (Figure 4). Several other redox-sensitive genes were upregulated after short-term exposure to the high concentration of curcumin, like AREG, ATF4, EGR1, FGFR1 (Table 2) [44]. However, not all phase-II genes were upregulated by curcumin. In this study several GSTs were downregulated, e.g. GSTT2 and GSTM4 (table 2), GSTZ1 (after exposure for 24 hours). In addition, GCLM (the light regulatory subunit of gamma-glutamylcysteine synthetase) was downregulated by curcumin with maximal downregulation by the low concentration of curcumin after 12 hours (1.8-fold) and by the high concentration of curcumin after 24 hours (1.9-fold). Gamma-glutamyl hydrolase (GGH) was downregulated 1.5-fold by the high curcumin concentration after 24 hours. Several phase-I metabolism genes were downregulated by curcumin, for example expression of CYP1B1 was inhibited by the high curcumin concentration after 3 hours (Table 2) and expression of some other cytochrome P450 genes was downregulated after exposure for 12 hours to the high curcumin concentration (data not shown). Expression of the aryl hydrocarbon receptor (AHR) was also slightly downregulated at this point (1.4-fold).

Response of genes involved in phase II metabolism to exposure to curcumin; A: low concentration, B: high concentration. AKR1C1: aldo-keto reductase family 1 member C1; NQO1: NAD(P)H dehydrogenase, quinone 1; NQO2: NAD(P)H dehydrogenase, quinone 2; TXNRD1: thioredoxin reductase 1; EPHX1: epoxide hydrolase 1.

AKR1C1 gene expression changes in response to curcumin [A: low concentration, B: high concentration], measured with cDNA microarray (black bars), and measured by real time RT-PCR normalized to beta-actin expression (white bars) or to GAPDH expression (grey bars). Expression is plotted as induction factors (level in treated cells compared to untreated cells) (mean ± standard deviation).

A second group of genes that showed a similar response to curcumin exposure in this study were metallothionein genes. All were downregulated by curcumin, specifically at the 12 and 24 hour time points (Figure 6). In addition, the metallothionein downregulation was preceded by a downregulation of MTF1, the transcription factor regulating metallothionein expression, by the low curcumin concentration (Table 1). In human colon cancer cells metallothionein was found to be a possible proliferation marker, as it is higher expressed in proliferating cells and its expression is cell cycle regulated [45]. Expression of metallothionein in colon tumors is associated with poor prognosis [46].

Response of metallothionein (MT) genes 1B-1L to exposure to curcumin; A: low concentration, B: high concentration.

Another group of genes that responded to curcumin exposure in a similar fashion were tubulin genes. These genes were downregulated after 3 hours, but were upregulated after 48 hours of exposure to the high curcumin concentration (Figure 7). Alpha tubulin has been reported to be differentially expressed in normal colon compared to colon tumors [47] and alpha tubulin was downregulated during differentiation of HT29-D4 cells [48].

Response of tubulin genes in HT29 cells to exposure to curcumin (high concentration); TUBA1: tubulin, alpha 1; TUBA3: tubulin, alpha 3; TUBB2: tubulin, beta 2.

Several proteasome genes were upregulated by the high concentration of curcumin after 12 or 24 hours: PSMA1, PSMA7, PSMB2 after 12 hours and PSMB6, PSMC4, PSMD2 after 24 hours (data not shown).

In addition to HT29 cells, expression changes at early time points of exposure to curcumin were also studied in Caco-2 cells. Fewer genes were differentially expressed in these cells than in HT29 cells, between 20 and 50 genes were up- or downregulated more than 1.5 fold after short-term exposure to curcumin. The most striking response in Caco-2 cells was a strong downregulation of CYP1A1 expression after exposure to both concentrations of curcumin for 3 and 6 hours (Table 3). CYP1A1 expression did not change in HT29 cells in response to curcumin, however CYP1B1 expression was downregulated by the high concentration of curcumin (Table 2). Some genes showed a similar response in Caco-2 cells as in HT29 cells, for example EGR1 and AKR1C1. However, upregulation in Caco-2 cells (Table 3) was not as strong as in HT29 cells (Figure 4 and 3). Expression of TLOC1 and POLR1C was downregulated in Caco-2 cells (Table 3) as well as in HT29 cells after short-term exposure (Table 1 and 2). Also upregulation of ETR101 and an integrin gene in Caco-2 cells is consistent with the results in HT29 cells. In Caco-2 cells several transcription factors were upregulated by curcumin, like NR2F1, NR4A1 and MTF1, whereas transcription factor FOSB and transcription-related genes HNRPA0 and POLR1C were downregulated (Table 3).

Early response genes in Caco-2 cells after exposure to curcumin.

Fold change (treated/control)
30 μmol/L 100 μmol/L

Accession number Gene Name Gene Symbol 3 h 6 h 3 h 6 h Involved in
AA418907 cytochrome P450, family 1, subfamily A, polypeptide 1 CYP1A1 -2.93 -3.00 -2.89 -6.39 Xenobiotic metabolism
R93124 aldo-keto reductase family 1, member C1 (dihydrodiol dehydrogenase 1) AKR1C1 1.05 1.60 -1.15 1.18 Xenobiotic metabolism
AA486533 early growth response 1 EGR1 1.11 1.03 2.49 2.30 Gene transcription
AA452909 nuclear receptor subfamily 2, group F, member 1 NR2F1 -1.40 NA1 -1.20 1.83 Gene transcription
N94487 nuclear receptor subfamily 4, group A, member 1 NR4A1 -1.19 1.09 1.45 1.08 Gene transcription
AA448256 metal-regulatory transcription factor 1 MTF1 0.28 0.34 0.09 0.74 Gene transcription
T61948 FBJ murine osteosarcoma viral oncogene homolog B FOSB -1.37 1.21 -1.22 -1.62 Gene transcription
AA599176 heterogeneous nuclear ribonucleoprotein A0 HNRPA0 -1.05 -1.31 -1.12 -1.58 Gene transcription
AA733038 polymerase (RNA) I polypeptide C, 30 kDa POLR1C -1.22 -1.28 -1.34 -1.73 Gene transcription
AA464970 phospholipase C, beta 2 PLCB2 NA 1.93 1.08 1.05 Signal transduction
AA450205 translocation protein 1 TLOC1 1.15 -1.36 -1.20 -1.62 Protein transport
AA419177 solute carrier family 7, member 5 SLC7A5 -1.27 -1.54 -1.56 -1.99 Transport
AA425451 integrin, alpha E ITGAE NA 1.43 1.54 NA Cell adhesion
AA888148 tubulin, beta 2 TUBB2 -1.25 -1.25 -1.57 -1.97 Cytoskeleton
AA496359 immediate early protein ETR101 -1.09 -1.02 1.77 1.52 unknown

1 NA: not available

When HT29 cells were exposed to the low curcumin concentration for 6 hours, the number of cells was decreased compared to cells not exposed to curcumin (Figure 8). However, after 24 hours a difference with cells that were not exposed to curcumin was no longer seen.

Number of HT29 cells after exposure to curcumin at different time points. Black bars: untreated cells, grey bars: low concentration, white bars: high concentration. Values are mean ± standard deviation.

Exposure of the cells to the high concentration of curcumin caused a decrease in cell number at all time points (Figure 8). After exposure to curcumin for 3 or 6 hours a significant decrease in the percentage of cells in the G1 phase of the cell cycle and a significant increase in the percentage of cells in the G2/M phase of the cell cycle was observed (Figure 9). This effect was dose-dependent. At the 24 hour time point the percentage of cells in the G1 phase was decreased significantly and the percentages of cells in both the S and the G2/M phase were increased significantly (Figure 9). When absolute cell numbers were recalculated from the percentages, the number of cells in G1 and S phase decreased significantly after exposure for 3 hours (only for the high curcumin concentration) or 6 hours. The absolute number of cells in G2/M phase remained constant. After exposure to the low concentration of curcumin for 24 hours, the absolute number of cells in G1 phase decreased significantly, but the absolute number of cells in S and G2/M phase increased significantly (data not shown).

Cell cycle distribution of HT29 cells exposed to curcumin for 3, 6 or 24 hours, measured with propidium iodide assay. Black bars: untreated cells, grey bars: low concentration, white bars: high concentration. Values are mean ± standard deviation (n = 3 or 4). *: significantly different from untreated cells, P < 0.05

Discussion

In this study the effect of curcumin on gene expression in HT29 colon cancer cells is determined at different time points. In addition, gene expression changes in response to short term exposure to curcumin were also studied in Caco-2 cells. The concentrations of curcumin (in the micromolar range) used in this study reflect the in vivo situation. Based on almost complete faecal excretion of curcumin and its metabolites [24,25] and assuming a volume of one liter in the stomach and a 10-fold dilution from stomach to colon, the concentration of curcumin in the intestinal lumen can be as high as 270 μmol/L after consuming one gram of curcumin, with a meal and/or as supplements. A similar concentration range for curcumin in colon was reported by Wortelboer et al. [49]. The highest concentration used in this study was 100 μmol/L. Exposure of HT29 cells to this curcumin concentration resulted in a decrease in cell number and floating cells were seen in the culture flasks after 24 or 48 hours. Therefore, it is possible that gene expression changes after exposure to 100 μmol/L curcumin for longer time periods (24, 48 hours) were related to toxic effects of the compound.

The study set-up with a time series of exposure to curcumin allowed us to identify early response genes and gene expression changes over time. Early response genes were identified; these were involved in DNA repair, signal transduction, transcription regulation, cell adhesion, xenobiotic metabolism. Also, genes involved in cell cycle control were up- or downregulated by short-term exposure to curcumin. Additionally, we have shown that exposure of HT29 cells to curcumin for 3 or 6 hours results in a decrease in percentage of cells in G1 phase and an increase in percentage of cells in G2/M phase. Other studies also reported a cell cycle arrest in the G2/M phase in colon cancer cells after curcumin exposure (similar concentrations, in micromolar range), but only exposure periods of 12 hours or longer were studied [12-14]. In addition, after 24 hours an increase in the percentage of cells in the S phase was observed. Hanif et al. [12] also found an increase in the percentage of cells in the S-phase in HT29 cells after exposure to curcumin for 24 hours. Interestingly, the initial decrease in cell number after short-term exposure to the low concentration of curcumin was no longer seen after 24 hours. At this time point, exposure to the low curcumin concentration had resulted in a significant increase in the absolute number of cells in the S and G2/M phase when compared to cells that were not exposed to curcumin. This indicates that the cells may have overcome the cell growth-inhibiting effect of curcumin at this time point. Consistent with this cell-physiological observation, changes in expression of cell cycle genes were not found after exposure to the low curcumin concentration for 24 hours.

The downregulation of expression of PLK by curcumin as found in this study could be involved in the G2/M arrest. Recently it was shown that downregulation of PLK in colon cancer cell line SW480 resulted in an increase in the percentage of cells in the G2/M phase [50]. Also in HeLa cells, depletion of PLK resulted in G2/M phase arrest and apoptosis [51]. Also other genes involved in transition through the G2/M phase of the cell cycle were downregulated by curcumin, like BUB1B and MAD2L1. Downregulation of HDAC1 by curcumin could possibly also be related to the S/G2/M phase arrest, since exposure of leukaemia cells to histone deacetylase inhibitors resulted in a G2/M phase arrest [52]. In cancer cells inhibition of histone deacetylase activity can cause growth arrest and apoptosis, and thus inhibit carcinogenesis [53]. Recently, it was reported that HDAC1 is associated with PCNA [54]. PCNA, involved in DNA replication, was also downregulated by curcumin.

It is known that dietary curcumin can induce activity of antioxidant or phase-II enzymes in livers of rats treated with 1,2-dimethylhydrazine (DMH) or in liver and kidney of ddY mice [55,56]. In this study, it was found that curcumin also has an effect on phase-II enzyme gene expression level in colonic cells. Induction of phase-II genes can be a mechanism to protect against development of cancer [57]. Expression of most of these genes is induced through the antioxidant response element (also known as EpRE) [58,59]. In human bronchial epithelial cells curcumin increased EpRE-binding activity and changed the composition of the EpRE-binding complex [60]. One of the transcription factors present in the EpRE-binding complex is Nrf2. Curcumin exposure resulted in an increase in nuclear Nrf2 content in bronchial epithelial cells [60] and in increased expression of Nrf2 in renal epithelial cells [61]. Unfortunately, Nrf2 was not represented on the cDNA microarray used in this study. Possibly, transcription factor ATF4 can form a dimer with Nrf2 [62], suggesting a role in EpRE-mediated effects. Indeed, curcumin was found to upregulate expression of ATF4 at the early time points, consistent with its role as an initiator of the curcumin effect through gene transcription modulation. Curcumin downregulated expression of several GST genes. This corresponds to an earlier observation that GSTP1 expression was downregulated by curcumin in leukaemia cells [63]. Overall, the effect of curcumin on the phase-II and redox-sensitive genes could indicate an increased level of protection of the cells against oxidative stress, consistent with its function as an anti-oxidant. In addition to its effect on phase-II biotransformation genes, it is suggested that curcumin can inhibit activation of carcinogens by cytochrome P450 enzymes [64]. In Caco-2 cells curcumin caused a strong downregulation of CYP1A1 gene expression. In HT29 cells curcumin downregulated expression of the CYP1B1 and aryl hydrocarbon receptor genes.

Changes in expression of cell cycle-related genes together with induced cell cycle arrest as well as induction of phase-II genes can be mechanisms of colon cancer prevention. Differential expression of tubulin genes and of proteasome genes in response to curcumin could also be interesting in view of the anticarcinogenic effect of curcumin. Microtubules formed by tubulin are important for spindle formation during cell division [65]. Protein degradation by proteasomes plays a role in cell cycle control and apoptosis and proteasome inhibitors are tested for use in anticancer therapy [66]. Furthermore, downregulation of p53 expression could also play a role. Rodrigues et al. reported that in HT29 cells mutated p53 is highly expressed [67]. Similar to HT29 cells, Caco-2 cells also contain a mutated p53 gene [68]. Overexpression of p53 is found in colon tumors and is associated with low chances on disease-free survival [69]. Recently, it was described that EGR1 is required for p53 tumor suppression [70]. EGR1 was the gene most induced by curcumin at the early time points. In HT29 cells treated with a mitosis-inhibiting peptide expression of EGR1 increased already after 20 minutes [71]. Therefore, also in the case of curcumin early upregulation of EGR1 could be related to growth inhibition. Interestingly, endothelial cells contrasted with HT29 cells in that induced EGR1 expression in these cells was found to be suppressed by curcumin [72]. Another interesting observation in our study is the differential expression of p16(INK4) (CDKN2A), RB1 and p53 in response to curcumin. Functions of these genes, working together in a signalling network regulating cell cycle, are often impaired in cancer cells [73]. Other interesting pathways that were influenced by curcumin are MAPK signal transduction and DNA repair.

In the study by Mariadason et al. [74] gene expression changes in SW620 cells in response to curcumin (25 μM) were measured at time points between 30 minutes and 48 hours. Of the 6253 genes analysed by Mariadason et al. [74], 1350 genes are present in our study. In this subset, genes that were more than 1.5-fold up- or downregulated were identified and compared with our data. Several genes showed a similar response to curcumin in SW620 and HT29 cells (i.e. up- or downregulation in both cell types), for example activating transcription factor 4 (ATF4; up), vascular endothelial growth factor (VEGF; up), proteasome subunit PSMD8 (up), glutathione reductase (GSR; up), adenosine kinase (ADK; up), RAN binding protein 2-like 1 (RANBP2L1; down), prostaglandin D synthase (PTGDS; up), signal sequence receptor alpha (SSR1; down), ribonuclease 4 (RNASE4; down). Although both studies used colon cancer cells, there were also differences in response to curcumin between SW620 and HT29 cells. However, the fact that cell lines differ in their response to a food component is not surprising, since we found considerable differences between twelve human colon cancer cell lines when their RNA expression profiles under standard optimal culture conditions were compared (unpublished results). Moreover, the fact that SW620 colon cancer cells were derived from a metastasis of colon cancer, whereas HT29 cells were derived from a primary colon tumor, implies that substantial physiological differences exist between these cell lines, which could contribute to the differences in gene expression profile in response to curcumin. In our study expression of fewer genes was changed in Caco-2 cells than in HT29 cells after short-term exposure to curcumin. Consistent with this observation, it was reported earlier that HT29 cells are more sensitive to curcumin than Caco-2 cells [75].

Curcumin is an anti-inflammatory agent and can act as a natural non-steroidal anti-inflammatory drug (NSAID). Some of the genes differentially expressed in response to curcumin in our study were also differentially expressed in human colon cancer cells lines in response to other non-steroidal compounds with anti-inflammatory action (NSAIDs) such as aspirin or sulindac [76]. For example, expression of ASNS, ATF4 and MTHFD2 was upregulated by sulindac and expression of BUB1B and PDCD2 was downregulated by aspirin. Similar changes in gene expression were found in our study in response to curcumin. In a microarray study with rat colon carcinoma cells, it was found that a large proportion of the genes differentially expressed in response to aspirin were also differentially expressed in the same direction in response to butyrate [27]. In invasive human colonocytes butyrate inhibited urokinase plasminogen activator (uPA) activity, and downregulated PCNA and TP53 levels after exposure for 12–18 hours [77]. In our study a similar response was found, as curcumin reduced expression of PCNA, TP53 and PLAU (uPA). This indicates that there may be some overlap in response and in mechanism of action between different NSAIDs like curcumin, aspirin and sulindac, but also between NSAIDs and butyrate. However, both aspirin and sulindac increased expression of several metallothionein genes, which were downregulated by curcumin in our study [76].

Two studies looked at the effect of curcumin and demethoxycurcumin on gene expression in human umbilical vein endothelial cells using microarrays [78,79]. In these cells curcumin or demethoxycurcumin upregulated expression of cyclin-dependent kinase inhibitor 2D (p19INK4D) and downregulated expression of PCNA, HDAC2, MAP2K1 and PLAU. This is consistent with our results, where curcumin exposure resulted in an increased expression of cyclin-dependent kinase inhibitor 2A (CDKN2A) and a decreased expression of PCNA, HDAC1, MAP2K2 and PLAU. It is interesting to see that human cell lines derived from different origin (colon and umbilical vein) show similar responses to curcumin.

Conclusions

In conclusion, this study describes changes in gene expression profiles in colon cancer cells in response to exposure to curcumin and relates these gene expression changes to functional and physiological processes. To our knowledge, this is the first study that uses a genomics approach to investigate in detail the mechanisms of effects of curcumin in colon cancer cells. Studying both early and later time points allowed us to identify genes that changed in expression with time of exposure. Some known effects of curcumin were confirmed (G2/M cell cycle arrest, induction of phase-II genes) and the existing knowledge was extended with extra information (e.g. time points of the observed changes and genes involved or linked to these physiological effects). Also, potential new leads to mechanisms explaining the biological activity of curcumin were identified, for example the effect on tubulin genes and differential expression of p16(INK4)/TP53/RB1. Studying expression changes of thousands of genes has provided increased insight into the mechanism of action of curcumin in colon cancer cells, helping us to understand how this compound can protect against development of colon cancer.

Authors' contributions

ME, ET, YS and SH carried out the cell culture, RNA isolation, microarray studies and flow cytometry studies. ME carried out the real-time PCR, the microarray data analysis and drafted the manuscript. PB, JA and BO participated in design of the study and of the manuscript. All authors read and approved the manuscript.

Acknowledgements

The authors thank Elles Klein Koerkamp for assistance in the flow cytometry analyses, Dr. Sabina Bijlsma for assistance in principal component analysis and Dr. Frank Schuren for setting-up the microarray facility. This study was financially supported by EU-grant QLRT-1999-00706.

Mohandas KM Desai DC Epidemiology of digestive tract cancers in India. V. Large and small bowel Indian J Gastroenterol 1999 18 118 121 10407566 Rao CV Simi B Reddy BS Inhibition by dietary curcumin of azoxymethane-induced ornithine decarboxylase, tyrosine protein kinase, arachidonic acid metabolism and aberrant crypt foci formation in the rat colon Carcinogenesis 1993 14 2219 2225 8242846 Rao CV Rivenson A Simi B Reddy BS Chemoprevention of colon carcinogenesis by dietary curcumin, a naturally occurring plant phenolic compound Cancer Res 1995 55 259 266 7812955 Pereira MA Grubbs CJ Barnes LH Li H Olson GR Eto I Juliana M Whitaker LM Kelloff GJ Steele VE Lubet RA Effects of the phytochemicals, curcumin and quercetin, upon azoxymethane-induced colon cancer and 7,12-dimethylbenz[a]anthracene- induced mammary cancer in rats Carcinogenesis 1996 17 1305 1311 8681447 Mahmoud NN Carothers AM Grunberger D Bilinski RT Churchill MR Martucci C Newmark HL Bertagnolli MM Plant phenolics decrease intestinal tumors in an animal model of familial adenomatous polyposis Carcinogenesis 2000 21 921 927 10783313 10.1093/carcin/21.5.921 Perkins S Verschoyle RD Hill K Parveen I Threadgill MD Sharma RA Williams ML Steward WP Gescher AJ Chemopreventive Efficacy and Pharmacokinetics of Curcumin in the Min/+ Mouse, a Model of Familial Adenomatous Polyposis Cancer Epidemiol Biomarkers Prev 2002 11 535 540 12050094 Kawamori T Lubet R Steele VE Kelloff GJ Kaskey RB Rao CV Reddy BS Chemopreventive effect of curcumin, a naturally occurring anti- inflammatory agent, during the promotion/progression stages of colon cancer Cancer Res 1999 59 597 601 9973206 Huang MT Lou YR Ma W Newmark HL Reuhl KR Conney AH Inhibitory effects of dietary curcumin on forestomach, duodenal, and colon carcinogenesis in mice Cancer Res 1994 54 5841 5847 7954412 Chuang SE Cheng AL Lin JK Kuo ML Inhibition by curcumin of diethylnitrosamine-induced hepatic hyperplasia, inflammation, cellular gene products and cell-cycle- related proteins in rats Food Chem Toxicol 2000 38 991 995 11038236 10.1016/S0278-6915(00)00101-0 Collett GP Robson CN Mathers JC Campbell FC Curcumin modifies Apc(min) apoptosis resistance and inhibits 2-amino 1- methyl-6-phenylimidazo[4,5-b]pyridine (PhIP) induced tumour formation in Apc(min) mice Carcinogenesis 2001 22 821 825 11323404 10.1093/carcin/22.5.821 Singh SV Hu X Srivastava SK Singh M Xia H Orchard JL Zaren HA Mechanism of inhibition of benzo[a]pyrene-induced forestomach cancer in mice by dietary curcumin Carcinogenesis 1998 19 1357 1360 9744529 10.1093/carcin/19.8.1357 Hanif R Qiao L Shiff SJ Rigas B Curcumin, a natural plant phenolic food additive, inhibits cell proliferation and induces cell cycle changes in colon adenocarcinoma cell lines by a prostaglandin-independent pathway J Lab Clin Med 1997 130 576 584 9422331 Chen H Zhang ZS Zhang YL Zhou DY Curcumin inhibits cell proliferation by interfering with the cell cycle and inducing apoptosis in colon carcinoma cells Anticancer Res 1999 19 3675 3680 10625938 Moragoda L Jaszewski R Majumdar AP Curcumin induced modulation of cell cycle and apoptosis in gastric and colon cancer cells Anticancer Res 2001 21 873 878 11396178 Choudhuri T Pal S Agwarwal ML Das T Sa G Curcumin induces apoptosis in human breast cancer cells through p53- dependent Bax induction FEBS Lett 2002 512 334 340 11852106 10.1016/S0014-5793(02)02292-5 Deeb D Xu YX Jiang H Gao X Janakiraman N Chapman RA Gautam SC Curcumin (Diferuloyl-Methane) Enhances Tumor Necrosis Factor-related Apoptosis-inducing Ligand-induced Apoptosis in LNCaP Prostate Cancer Cells Mol Cancer Ther 2003 2 95 103 12533677 Mehta K Pantazis P McQueen T Aggarwal BB Antiproliferative effect of curcumin (diferuloylmethane) against human breast tumor cell lines Anticancer Drugs 1997 8 470 481 9215611 Pan MH Chang WL Lin-Shiau SY Ho CT Lin JK Induction of apoptosis by garcinol and curcumin through cytochrome c release and activation of caspases in human leukemia HL-60 cells J Agric Food Chem 2001 49 1464 1474 11312881 10.1021/jf001129v Mukhopadhyay A Banerjee S Stafford LJ Xia C Liu M Aggarwal BB Curcumin-induced suppression of cell proliferation correlates with down- regulation of cyclin D1 expression and CDK4-mediated retinoblastoma protein phosphorylation Oncogene 2002 21 8852 8861 12483537 10.1038/sj.onc.1206048 Leu TH Maa MC The molecular mechanisms for the antitumorigenic effect of curcumin Curr Med Chem Anti-Canc Agents 2002 2 357 370 12678737 Lin JK Pan MH Lin-Shiau SY Recent studies on the biofunctions and biotransformations of curcumin Biofactors 2000 13 153 158 11237176 Aggarwal BB Kumar A Bharti AC Anticancer potential of curcumin: preclinical and clinical studies Anticancer Res 2003 23 363 398 12680238 Ravindranath V Chandrasekhara N Absorption and tissue distribution of curcumin in rats Toxicology 1980 16 259 265 7423534 10.1016/0300-483X(80)90122-5 Ravindranath V Chandrasekhara N Metabolism of curcumin--studies with [3H]curcumin Toxicology 1981 22 337 344 7342372 10.1016/0300-483X(81)90027-5 Ireson CR Jones DJ Orr S Coughtrie MW Boocock DJ Williams ML Farmer PB Steward WP Gescher AJ Metabolism of the cancer chemopreventive agent curcumin in human and rat intestine Cancer Epidemiol Biomarkers Prev 2002 11 105 111 11815407 Cheng AL Hsu CH Lin JK Hsu MM Ho YF Shen TS Ko JY Lin JT Lin BR Ming-Shiang W Yu HS Jee SH Chen GS Chen TM Chen CA Lai MK Pu YS Pan MH Wang YJ Tsai CC Hsieh CY Phase I clinical trial of curcumin, a chemopreventive agent, in patients with high-risk or pre-malignant lesions Anticancer Res 2001 21 2895 2900 11712783 Germann A Dihlmann S Hergenhahn M Von Knebel Doeberitz M Koesters R Expression profiling of CC531 colon carcinoma cells reveals similar regulation of beta-catenin target genes by both butyrate and aspirin Int J Cancer 2003 106 187 197 12800193 10.1002/ijc.11215 Kindermann B Doring F Pfaffl M Daniel H Identification of Genes Responsive to Intracellular Zinc Depletion in the Human Colon Adenocarcinoma Cell Line HT-29 J Nutr 2004 134 57 62 14704293 Pfaffl MW Horgan GW Dempfle L Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR Nucleic Acids Res 2002 30 e36. 11972351 10.1093/nar/30.9.e36 Yang YH Dudoit S Luu P Lin DM Peng V Ngai J Speed TP Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation Nucleic Acids Res 2002 30 e15. 11842121 Siu F Bain PJ LeBlanc-Chaffin R Chen H Kilberg MS ATF4 is a mediator of the nutrient-sensing response pathway that activates the human asparagine synthetase gene J Biol Chem 2002 277 24120 24127 11960987 10.1074/jbc.M201959200 Liu C Calogero A Ragona G Adamson E Mercola D EGR-1, the reluctant suppression factor: EGR-1 is known to function in the regulation of growth, differentiation, and also has significant tumor suppressor activity and a mechanism involving the induction of TGF-beta1 is postulated to account for this suppressor activity Crit Rev Oncog 1996 7 101 125 9109500 Takahashi T Sano B Nagata T Kato H Sugiyama Y Kunieda K Kimura M Okano Y Saji S Polo-like kinase 1 (PLK1) is overexpressed in primary colorectal cancers Cancer Sci 2003 94 148 152 12708489 Spankuch-Schmitt B Wolf G Solbach C Loibl S Knecht R Stegmuller M von Minckwitz G Kaufmann M Strebhardt K Downregulation of human polo-like kinase activity by antisense oligonucleotides induces growth inhibition in cancer cells Oncogene 2002 21 3162 3171 12082631 10.1038/sj.onc.1205412 Papadopoulou S Scorilas A Yotis J Arnogianaki N Plataniotis G Agnanti N Talieri M Significance of urokinase-type plasminogen activator and plasminogen activator inhibitor-1 (PAI-1) expression in human colorectal carcinomas Tumour Biol 2002 23 170 178 12218297 10.1159/000064033 Hashida H Takabayashi A Tokuhara T Hattori N Taki T Hasegawa H Satoh S Kobayashi N Yamaoka Y Miyake M Clinical significance of transmembrane 4 superfamily in colon cancer Br J Cancer 2003 89 158 167 12838318 10.1038/sj.bjc.6601015 Kivela AJ Saarnio J Karttunen TJ Kivela J Parkkila AK Pastorekova S Pastorek J Waheed A Sly WS Parkkila TS Rajaniemi H Differential expression of cytoplasmic carbonic anhydrases, CA I and II, and membrane-associated isozymes, CA IX and XII, in normal mucosa of large intestine and in colorectal tumors Dig Dis Sci 2001 46 2179 2186 11680594 10.1023/A:1011910931210 Kitahara O Furukawa Y Tanaka T Kihara C Ono K Yanagawa R Nita ME Takagi T Nakamura Y Tsunoda T Alterations of gene expression during colorectal carcinogenesis revealed by cDNA microarrays after laser-capture microdissection of tumor tissues and normal epithelia Cancer Res 2001 61 3544 3549 11325815 Levine AJ p53, the cellular gatekeeper for growth and division Cell 1997 88 323 331 9039259 10.1016/S0092-8674(00)81871-1 Pufulete M Emery PW Sanders TA Folate, DNA methylation and colo-rectal cancer Proc Nutr Soc 2003 62 437 445 14506892 Lee S Baek M Yang H Bang YJ Kim WH Ha JH Kim DK Jeoung DI Identification of genes differentially expressed between gastric cancers and normal gastric mucosa with cDNA microarrays Cancer Lett 2002 184 197 206 12127692 10.1016/S0304-3835(02)00197-0 Burczynski ME Lin HK Penning TM Isoform-specific induction of a human aldo-keto reductase by polycyclic aromatic hydrocarbons (PAHs), electrophiles, and oxidative stress: implications for the alternative pathway of PAH activation catalyzed by human dihydrodiol dehydrogenase Cancer Res 1999 59 607 614 9973208 Bonnesen C Eggleston IM Hayes JD Dietary indoles and isothiocyanates that are generated from cruciferous vegetables can both stimulate apoptosis and confer protection against DNA damage in human colon cell lines Cancer Res 2001 61 6120 6130 11507062 Allen RG Tresini M Oxidative stress and gene regulation Free Radic Biol Med 2000 28 463 499 10699758 10.1016/S0891-5849(99)00242-7 Nagel WW Vallee BL Cell cycle regulation of metallothionein in human colonic cancer cells Proc Natl Acad Sci U S A 1995 92 579 583 7831333 Hishikawa Y Kohno H Ueda S Kimoto T Dhar DK Kubota H Tachibana M Koji T Nagasue N Expression of metallothionein in colorectal cancers and synchronous liver metastases Oncology 2001 61 162 167 11528256 10.1159/000055368 van Belzen N Dinjens WN Eussen BH Bosman FT Expression of differentiation-related genes in colorectal cancer: possible implications for prognosis Histol Histopathol 1998 13 1233 1242 9810514 van Belzen N Diesveld MP van der Made AC Nozawa Y Dinjens WN Vlietstra R Trapman J Bosman FT Identification of mRNAs that show modulated expression during colon carcinoma cell differentiation Eur J Biochem 1995 234 843 848 8575443 Wortelboer HM Usta M van der Velde AE Boersma MG Spenkelink B van Zanden JJ Rietjens IM van Bladeren PJ Cnubben NH Interplay between MRP inhibition and metabolism of MRP inhibitors: the case of curcumin Chem Res Toxicol 2003 16 1642 1651 14680379 10.1021/tx034101x Spankuch-Schmitt B Bereiter-Hahn J Kaufmann M Strebhardt K Effect of RNA silencing of polo-like kinase-1 (PLK1) on apoptosis and spindle formation in human cancer cells J Natl Cancer Inst 2002 94 1863 1877 12488480 10.1093/jnci/94.24.1863 Liu X Erikson RL Polo-like kinase (Plk)1 depletion induces apoptosis in cancer cells Proc Natl Acad Sci U S A 2003 100 5789 94 12732729 10.1073/pnas.1031523100 Peart MJ Tainton KM Ruefli AA Dear AE Sedelies KA O'Reilly LA Waterhouse NJ Trapani JA Johnstone RW Novel mechanisms of apoptosis induced by histone deacetylase inhibitors Cancer Res 2003 63 4460 4471 12907619 Kim DH Kim M Kwon HJ Histone deacetylase in carcinogenesis and its inhibitors as anti-cancer agents J Biochem Mol Biol 2003 36 110 119 12542981 Milutinovic S Zhuang Q Szyf M Proliferating cell nuclear antigen associates with histone deacetylase activity, integrating DNA replication and chromatin modification J Biol Chem 2002 277 20974 20978 11929879 10.1074/jbc.M202504200 Devasena T Rajasekaran KN Menon VP Bis-1,7-(2-hydroxyphenyl)-hepta-1,6-diene-3,5-Dione (a curcumin analog) ameliorates dmh-induced hepatic oxidative stress during colon carcinogenesis Pharmacol Res 2002 46 39 45 12208119 10.1016/S1043-6618(02)00043-9 Iqbal M Sharma SD Okazaki Y Fujisawa M Okada S Dietary supplementation of curcumin enhances antioxidant and phase II metabolizing enzymes in ddY male mice: possible role in protection against chemical carcinogenesis and toxicity Pharmacol Toxicol 2003 92 33 38 12710595 Talalay P Chemoprotection against cancer by induction of phase 2 enzymes Biofactors 2000 12 5 11 11216505 Hayes JD McMahon M Molecular basis for the contribution of the antioxidant responsive element to cancer chemoprevention Cancer Lett 2001 174 103 113 11689285 10.1016/S0304-3835(01)00695-4 Surh YJ Cancer chemoprevention with dietary phytochemicals Nat Rev Cancer 2003 3 768 780 14570043 10.1038/nrc1189 Dickinson DA Iles KE Zhang H Blank V Forman HJ Curcumin alters EpRE and AP-1 binding complexes and elevates glutamate-cysteine ligase gene expression Faseb J 2003 17 473 475 12514113 Balogun E Hoque M Gong P Killeen E Green CJ Foresti R Alam J Motterlini R Curcumin activates the haem oxygenase-1 gene via regulation of Nrf2 and the antioxidant-responsive element Biochem J 2003 371 887 895 12570874 10.1042/BJ20021619 He CH Gong P Hu B Stewart D Choi ME Choi AM Alam J Identification of activating transcription factor 4 (ATF4) as an Nrf2-interacting protein. Implication for heme oxygenase-1 gene regulation J Biol Chem 2001 276 20858 20865 11274184 10.1074/jbc.M101198200 Duvoix A Morceau F Delhalle S Schmitz M Schnekenburger M Galteau MM Dicato M Diederich M Induction of apoptosis by curcumin: mediation by glutathione S-transferase P1-1 inhibition Biochem Pharmacol 2003 66 1475 1483 14555224 10.1016/S0006-2952(03)00501-X Thapliyal R Maru GB Inhibition of cytochrome P450 isozymes by curcumins in vitro and in vivo Food Chem Toxicol 2001 39 541 547 11346483 10.1016/S0278-6915(00)00165-4 Nakaseko Y Yanagida M Cell biology. Cytoskeleton in the cell cycle Nature 2001 412 291 292 11460149 10.1038/35085684 Adams J The proteasome: structure, function, and role in the cell Cancer Treat Rev 2003 29 Suppl 1 3 9 12738238 10.1016/S0305-7372(03)00081-1 Rodrigues NR Rowan A Smith ME Kerr IB Bodmer WF Gannon JV Lane DP p53 mutations in colorectal cancer Proc Natl Acad Sci U S A 1990 87 7555 7559 1699228 Gartel AL Feliciano C Tyner AL A new method for determining the status of p53 in tumor cell lines of different origin Oncol Res 2003 13 405 408 12725531 Diez M Medrano M Muguerza JM Ramos P Hernandez P Villeta R Martin A Noguerales F Ruiz A Granell J Influence of tumor localization on the prognostic value of P53 protein in colorectal adenocarcinomas Anticancer Res 2000 20 3907 3912 11268475 Krones-Herzig A Adamson E Mercola D Early growth response 1 protein, an upstream gatekeeper of the p53 tumor suppressor, controls replicative senescence Proc Natl Acad Sci U S A 2003 100 3233 3238 12629205 10.1073/pnas.2628034100 Reichelt WH Yndestad A Wright MS Elgjo K Haug T Reichelt KL The colon mitosis-inhibitor pyroglutamyl-histidyl-glycine alters expression of immediate-early cancer-related genes in HT-29 cells Anticancer Res 2003 23 1229 1234 12820376 Pendurthi UR Rao LV Suppression of transcription factor Egr-1 by curcumin Thromb Res 2000 97 179 189 10674404 10.1016/S0049-3848(99)00148-6 Lowe SW Sherr CJ Tumor suppression by Ink4a-Arf: progress and puzzles Curr Opin Genet Dev 2003 13 77 83 12573439 10.1016/S0959-437X(02)00013-8 Mariadason JM Corner GA Augenlicht LH Genetic reprogramming in pathways of colonic cell maturation induced by short chain fatty acids: comparison with trichostatin A, sulindac, and curcumin and implications for chemoprevention of colon cancer Cancer Res 2000 60 4561 4572 10969808 Jiang MC Yang-Yen HF Yen JJ Lin JK Curcumin induces apoptosis in immortalized NIH 3T3 and malignant cancer cell lines Nutr Cancer 1996 26 111 120 8844727 Iizaka M Furukawa Y Tsunoda T Akashi H Ogawa M Nakamura Y Expression profile analysis of colon cancer cells in response to sulindac or aspirin Biochem Biophys Res Commun 2002 292 498 512 11906190 10.1006/bbrc.2002.6648 Emenaker NJ Calaf GM Cox D Basson MD Qureshi N Short-chain fatty acids inhibit invasive human colon cancer by modulating uPA, TIMP-1, TIMP-2, mutant p53, Bcl-2, Bax, p21 and PCNA protein expression in an in vitro cell culture model J Nutr 2001 131 3041S 6S 11694645 Kim JH Shim JS Lee SK Kim KW Rha SY Chung HC Kwon HJ Microarray-based analysis of anti-angiogenic activity of demethoxycurcumin on human umbilical vein endothelial cells: crucial involvement of the down-regulation of matrix metalloproteinase Jpn J Cancer Res 2002 93 1378 1385 12495478 Park MJ Kim EH Park IC Lee HC Woo SH Lee JY Hong YJ Rhee CH Choi SH Shim BS Lee SH Hong SI Curcumin inhibits cell cycle progression of immortalized human umbilical vein endothelial (ECV304) cells by up-regulating cyclin- dependent kinase inhibitor, p21WAF1/CIP1, p27KIP1 and p53 Int J Oncol 2002 21 379 383 12118335