Oxaliplatin induces different cellular and molecular chemoresistance patterns in colorectal cancer cell lines of identical origins
- Piroska Virag†1Email author,
- Eva Fischer-Fodor1,
- Maria Perde-Schrepler1,
- Ioana Brie1, 2,
- Corina Tatomir1,
- Loredana Balacescu1,
- Ioana Berindan-Neagoe1, 2,
- Bogdan Victor1 and
- Ovidiu Balacescu†1
© Virag et al.; licensee BioMed Central Ltd. 2013
Received: 6 March 2013
Accepted: 2 July 2013
Published: 16 July 2013
Cancer cells frequently adopt cellular and molecular alterations and acquire resistance to cytostatic drugs. Chemotherapy with oxaliplatin is among the leading treatments for colorectal cancer with a response rate of 50%, inducing intrastrand cross-links on the DNA. Despite of this drug’s efficiency, resistance develops in nearly all metastatic patients. Chemoresistance being of crucial importance for the drug’s clinical efficiency this study aimed to contribute to the identification and description of some cellular and molecular alterations induced by prolonged oxaliplatin therapy. Resistance to oxaliplatin was induced in Colo320 (Colo320R) and HT-29 (HT-29R) colorectal adenocarcinoma cell lines by exposing the cells to increasing concentrations of the drug. Alterations in morphology, cytotoxicity, DNA cross-links formation and gene expression profiles were assessed in the parental and resistant variants with microscopy, MTT, alkaline comet and pangenomic microarray assays, respectively.
Morphology analysis revealed epithelial-to-mesenchymal transition in the resistant vs parental cells suggesting alterations of the cells’ adhesion complexes, through which they acquire increased invasiveness and adherence. Cytotoxicity measurements demonstrated resistance to oxaliplatin in both cell lines; Colo320 being more sensitive than HT-29 to this drug (P < 0.001). The treatment with oxaliplatin caused major DNA cross-links in both parental cell lines; in Colo320R small amounts of DNA cross-links were still detectable, while in HT-29R not. We identified 441 differentially expressed genes in Colo320R and 613 in HT-29R as compared to their parental counterparts (at least 1.5 -fold up- or down- regulation, p < 0.05). More disrupted functions and pathways were detected in HT-29R cell line than in Colo320R, involving genes responsible for apoptosis inhibition, cellular proliferation and epithelial-to-mesenchymal transition. Several upstream regulators were detected as activated in HT-29R cell line, but not in Colo320R.
Our findings revealed a more resistant phenotype in HT-29R as compared to Colo320R and different cellular and molecular chemoresistance patterns induced by prolonged treatment with oxaliplatin in cell lines with identical origins (colorectal adenocarcinomas).
KeywordsOxaliplatin Colorectal cancer Chemoresistance DNA cross-links Gene expression profile
The last decade has brought major improvements in the treatment of cancers, but in spite of the efficacy of the cytostatic drugs, in time, cells adopt several cellular and molecular alterations, acquiring resistance. In colorectal cancer (CC) the 5 years survival rate remains lower than 10% in patients with metastasis, mainly due to resistance to the cytostatic drugs , regardless of the use of targeted molecular therapies in addition to standard chemotherapeutic regimens. The major treatment for metastatic CC is represented by 5-fluorouracil (5FU) and oxaliplatin (L-OHP). While 5FU inhibits thymidylate synthase during DNA replication , L-OHP acts as a bifunctional alkylating agent, covalently binding DNA and forming platinum-DNA adducts . The intrastrand cross-links formed by L-OHP being the most abundant lesions capable of blocking both replication and transcription of DNA, they are considered to cause the major cytotoxic lesions and being directly involved in the cancer cells death [4, 5].
L-OHP [(1R, 2R)-cyclohexane-1, 2-diamine] (ethanedioato-O, O’) platinum (II), a third generation platinum analogue, is the first compound that have proved to be efficient in the treatment of CC in patients displaying resistance to cisplatin (CDDP) and carboplatin (CBCDA) . Although the response rate to current systemic therapies is 50% resistance develops in almost all patients , limiting the drug’s therapeutic potential. Cells become resistant to platinum-based drugs through reduced cellular uptake, impaired DNA adducts formation, alterations in DNA repair genes such as ERCC1 and XRCC1 and modifications in the levels of copper transporters (ATP7A and ATP7B) [8–10].
Although in the last decade the gene expression profiling of human cancer cells provided valuable insight into the molecular targets of chemoresistance, the mechanisms involved in L-OHP resistance of CC are still poorly understood and the cellular and molecular alterations are not completely recognized.
Our study proposed to identify and describe some of the cellular and molecular alterations that occurred in CC cell lines with induced chemoresistance to L-OHP. In an earlier study, aiming to evaluate the differences in the behavior of the cells selected for L-OHP resistance compared to the sensitive ones, we assessed the cytotoxicity, apoptosis and induction of DNA damages by L-OHP in Colo320 CC cell line. We found lower toxicity, cellular death and fewer DNA damages, in the cells treated previously with L-OHP as compared to the parental ones . In the present study we performed a comparative study on two CC cell lines (Colo320 and HT-29) with identical origins (adenocarcinomas) and their L-OHP resistant counterparts (Colo320R and HT-29R) obtained by prolonged exposure to L-OHP. In addition, we analyzed the cells’ morphological features, DNA cross-links formation and the gene expression profiles. The DNA cross-links induction by L-OHP was determined indirectly with alkaline comet assay (CA), by introducing single strand breaks via ionizing radiation. The reduction of the single strand-breaks produced by the ionizing radiation in the cells treated with L-OHP quantitatively reflected the cross-links induced by the platinum compound. Gene expression profiling and subsequent validations with quantitative real-time PCR (qRT-PCR) were conducted in order to identify the molecular targets and pathways altered during chemoresistance acquirement. Finally, a comparative study was performed between the functions and pathways modulated by L-OHP treatment in the tested cell lines, in order to identify some common patterns of chemoresistance in CCs.
This study demonstrated that prolonged treatment with L-OHP induces several cellular and molecular targets and pathways, which lead to chemoresistance, but these alterations may differ consistently, even if the origins of the tested cell lines were identical. Results converged to the conclusion that CC cells change their morphology and cytotoxicity and react differently to L-OHP therapy by activating genes and upstream regulators that resulted in a more primitive, invasive and migratory, therefore more resistant phenotype. By the identification and description of the cellular and molecular alterations that occurred during the chemoresistance acquiring process the present study aims to contribute to a better understanding of this multifactorial process essential in our attempt to reverse chemoresistance and to identify potential targets for future therapies.
Cytotoxicity of L-OHP in parental and their L-OHP resistant variants (Colo320/Colo320R and HT-29/HT-29R)
7.546 ± 0.5970
20.85 ± 1.069
22.31 ± 2.717
56.80 ± 5.065
Evaluation of the DNA cross-links induced by L-OHP
The exposure of these cells to different concentrations of L-OHP prior to irradiation reduced significantly the length of the comet tails, due to the formation of L-OHP-DNA cross-links. This reduction was significant and dose-dependent in Colo320 (P < 0.0001), for both concentrations of the drug; for HT-29 this reduction was less important (P < 0.001 at 50 μg/ml and P < 0.05 at 100 μg/ml L-OHP) (Figure 2).
Gene expression profiles in the tested cells
Common core set of DE genes modulated by L-OHP in the tested cell lines
Thioredoxin-related transmembrane protein 4
Protein tyrosine phosphatase, receptor type, O
Lectin, galactoside-binding, soluble, 12
Keratin 18, transcript variant 1
ALU8_HUMAN (P39195) Alu subfamily SX sequence contamination warning entry, partial (8%)
PREDICTED: Homo sapiens similar to Keratin, type I cytoskeletal 18 (Cytokeratin-18) (CK-18) (Keratin-18) (K18) (LOC121054), mRNA [XR_019191]
Heat-shock protein beta-1 (HspB1) (Heat shock 27 kDa protein) (HSP 27) (Stress responsive protein 27) (SRP27) (Estrogen-regulated 24 kDa protein) (28 kDa heat shock protein)
PREDICTED: Homo sapiens similar to Keratin, type I cytoskeletal 18 (Cytokeratin-18) (CK-18) (Keratin-18) (K18) (LOC649233), mRNA [XR_018843
Homo sapiens cDNA FLJ11983 fis, clone HEMBB1001337
RR12_SPIMX (P42344) Chloroplast 30S ribosomal protein S12, partial (11%)
Q5U0N8_HUMAN (Q5U0N8) Keratin 18 (Cell proliferation-inducing protein 46), partial (46%) [THC2524582]
KRT18P55 keratin 18 pseudogene 55
PREDICTED: Homo sapiens similar to Keratin, type I cytoskeletal 18 (Cytokeratin-18) (CK-18) (Keratin-18) (K18) (LOC643471), mRNA [XR_018670]
PREDICTED: similar to Keratin, type I cytoskeletal 18 (Cytokeratin-18) (CK-18) (Keratin-18) (K18) (LOC731794), mRNA
Homo sapiens hypothetical protein LOC283666, mRNA (cDNA clone IMAGE:4415549), partial cds
PREDICTED: Homo sapiens similar to Keratin, type I cytoskeletal 18 (Cytokeratin-18) (CK-18) (Keratin-18) (K18) (LOC139060), mRNA [XR_018311]
Homo sapiens insulin-like growth factor 2 (somatomedin A) (IGF2), transcript variant 2, mRNA
ALU1_HUMAN (P39188) Alu subfamily J sequence contamination warning entry, partial (5%) [THC2663167]
Interleukin 1 receptor accessory protein
Growth factor receptor-bound protein 14
RNA binding protein with multiple splicing 2
Rho-related BTB domain containing 1
3-hydroxy-3-methylglutaryl-Coenzyme A synthase 1
Polymerase (RNA) I polypeptide D, 16 kDa
Carbonic anhydrase XII (CA12), transcript variant 1
Aldo-keto reductase family 1, member B10
Interferon induced transmembrane protein 3
Interferon induced transmembrane protein 4 pseudogene
Interferon induced with helicase C domain 1
Complement factor H-related 1
Calcium/calmodulin-dependent protein kinase II inhibitor 1
Sterile alpha motif domain containing 5
N-myc downstream regulated 1
Collagen, type IX, alpha 3
Iroquois homeobox 5
Identification of the biological pathways modulated by L-OHP
Top canonical pathways modulated by L-OHP in Colo320R
Canonical pathways Colo320R
Aryl hydrocarbon receptor signaling
HSPB3, CDKN2A, MYC, FOS, GSTM2, ALDH1A1, NQO1, TGFB2, HSPB1
CDKN2A, MET, FOS, MAP3K6, MAP3K13, ETS2, ELK3
LPS/IL-1 mediated inhibition of RXR function
CHST2, ABCB1, GSTM2, ALDH1A1, SLC27A2, HS3ST1, CPT1C, PLTP, HMGCS1, IL1RAP
Hepatic fibrosis/hepatic stellate cell activation
MET, TIMP1, ACTA2, TGFB2, IGFBP5, IL1RAP, COL3A1
Retinoate biosynthesis I
ALDH1A1, RBP7, AKR1C3
Methylglyoxal degradation III
S-methyl-5′-thioadenosine degradation II
DHRS3, RBP7, LPL
Cell Cycle: G1/S checkpoint regulation
CDKN2A, MYC, HDAC9, TGFB2
The visual cycle
ABCB1, GSTM2, ALDH1A1, FOXO1
Glutamine biosynthesis I
Thyroid cancer signaling
PPARG, MYC, NTF3
Mitochondrial L-carnitine shuttle pathway
Top canonical pathways modulated by L-OHP in HT-29R
Canonical pathways HT-29R
IFIT3, IFIT1, IFITM1, MX1, IRF9, PSMB8, STAT1
Methylglyoxal degradation III
AKR7A3, AKR1C1/AKR1C2, AKR1B10, AKR1B1
Cholesterol biosynthesis I
SQLE, DHCR7, MSMO1, CYP51A1
Cholesterol biosynthesis II (via 24,25-dihydrolanosterol)
SQLE, DHCR7, MSMO1, CYP51A1
Cholesterol biosynthesis III (via Desmosterol)
SQLE, DHCR7, MSMO1, CYP51A1
SCD, TTR, LDLR, SREBF1, AMBP, SERPINA1, PTGS2, TLR3, IL1RAP, CYP51A1, AGT
Mevalonate pathway I
IDI1, HMGCS2, HMGCS1
CKB, CKMT1A/CKMT1B, LCMT2
Activation of IRF by Cytosolic pattern recognition receptors
DHX58, IFIH1, IRF9, STAT1, IFIT2, ISG15
Tryptophan degradation X (Mammalian, via Tryptamine)
ALDH1A3, DDC, MAOA
Maturity onset diabetes of young (MODY) signaling
FOXA2, FABP1, HNF4A
DPEP1, PNPLA3, RARRES3, PTGS2, ALOX5
SCD, NR0B2, CYP2B6, HMGCS2, HNF4A
Colorectal cancer metastasis signaling
VEGFA, MMP7, RND3, TGFB1, FZD3, GNB5, TGFB2, PTGS2, TLR3, STAT1, FNBP1, WNT5A
Aryl hydrocarbon receptor signaling
TGM2, TFF1, NR0B2, ALDH1A3, TGFB1, TGFB2, ALDH18A1,NRIP1
p38 MAPK signaling
DDIT3, TGFB1, TGFB2, MKNK2, RPS6KA2, STAT1, IL1RAP
Choline biosynthesis III
Adenosine nucleotides degradation II
MMP7, CDH3, TGFB1, PPP2R3A, TGFBR3, FZD3, CD44, TGFB2, WNT5A
Identification of upstream regulators induced by L-OHP
Upstream regulators induced by L-OHP treatment in HT-29R
Predicted activation state
p-value of overlap
Target molecules in dataset
DDX60, DDX60L, DEFA1, HIP1, IFI6, IFIT1, IFIT2, IFIT3, IFIT5, IRF9, MT2A, OAS3, OASL, PARP9, PHLDA1, PLSCR1, RARRES3, RUNX2, SAMD9L, SLFN5, STAT1, TGFB1, XAF1
AGT, ALDH1A3,CLIC5, DPP4, HERC6, IFI27, IFI6, IFIT1, IFIT2, IFITM1, IRF9, ISG15, KCTD12, PSMB8, PSMB9, PTGS2, RUNX1, SLFN5, STAT1, TLR3, TYMP, USP18
DHX58, IFI44, IFI44L, IFI6, IFIH1, IFIT1,IFIT2, IFIT3, IFNE, ISG15, MX1, MX2, OASL, PHLDA1, PLK2, PTGS2, STAT1, TLR3, USP18
CXCR4, CYP51A1, DHCR7, HMGCS1, HMGCS2, IDI1, LDLR, PLD1, SCD, SQLE, SREBF1, STARD4,TGFB1
IFIT2, ISG15, LMO4, MX1, OASL, PLD1, PLK2, PTGS2, SH3BP5, SLCO3A1, STAT1, TGFB1,T LR3
IFIH1, IFIT2, MX2, PTGS2, RARRES3, STAT1, TLR3, USP18
IFI44, IFIT1, IFIT2, IFIT3, ISG15, STAT1
BST2, ISG15, MX1, PHLDA1, RND3
CAMK1D, CYP51A1, DHCR7, HMGCS1, HNF4A, IDI1, INSIG1, LDLR, MSMO1, NR0B2, PNPLA3, PTGS2, SCD, SQLE, SREBF1, STARD4,T GFB1, VEGFA
IFI6, IFITM3, IRF9, ISG15, MX1
Validation of microarray results by qRT-PCR
The correlation between microarray and qRT-PCR results
Pearson correlation r (p)
Protein tyrosine phosphatase, receptor type, O
Keratin 18, transcript variant 1
N-myc downstream regulated 1
Apoptosis, caspase activation inhibitor
WNT inhibitory factor 1
Inhibitor of DNA binding 1, dominant negative helix-loop-helix protein
Protein tyrosine phosphatase, receptor type, O
Keratin 18, transcript variant 1
N-myc downstream regulated 1
Transforming growth factor, beta 1
Midkine (neurite growth-promoting factor 2)
Cysteine-rich, angiogenic inducer, 61
Oxaliplatin has antitumoral activity against colorectal, breast, gastric, renal carcinomas and sarcomas  mainly by forming intrastrand cross-links that block DNA replication and transcription. Despite the improvements in the treatment outcome for CC, the tumor cells acquire resistance in time, decreasing the drug’s clinical efficiency.
To address this issue we assessed the morphology, cytotoxicity, DNA cross-links induction and gene expression profiles of two colorectal cell lines with identical origins (adenocarcinoma) with acquired resistance to L-OHP and their parental lines. According to our results the L-OHP resistant cells displayed altered cellular and molecular features as compared to the parental cells. In addition, notable differences were recorded between the functions and pathways modulated by L-OHP in the two tested cell lines.
Some of the morphological alterations we observed here: pseudopodia formation and adoption of fusiforme shape, suggesting an epithelial-to-mesenchymal transition (EMT) and an increased cell-to-cell distance in the HT-29R cells were also identified by Yang et al. in chemoresistant HT-29 cells . Part of the embryonic development, EMT appears to be involved in tumor progression and metastasis [15, 16], a process through which cells switch from the proliferative state to a more primitive, invasive and migratory mode. This conversion was proposed as a potential mechanism through which cells become chemoresistant, being known that the cytostatic drugs are more efficient on the highly proliferative cells . In our study Colo320R cells displayed (a mesenchymal phenotype, but adopted some) different characteristics after prolonged treatment with L-OHP. These cells, usually exposed in suspension, reacted to the prolonged treatment with the cytostatic drug by showing an increased tendency of adherence. Although our findings demonstrated different adaptations of the tested cell lines to the L-OHP treatment, a common feature was obvious: the alteration of the cellular adhesion complexes, suggesting higher invasiveness and attachment capacity.
The IC50 values obtained in the present study revealed that the prolonged treatment of the cells with increasing concentrations of the drug, up to the clinically relevant concentration (2 μmol/l), induced resistance in the treated cells as compared to the parental ones. Our results are comparable to other previous findings on parental and resistant Colo320 cells  and on sensitive Colo320 and HT-29 cells [7, 17].
The CA findings confirmed different behaviors of the tested cell lines to the prolonged treatment to L-OHP. Both parental cell lines were sensitive to 2 Gy of gamma irradiation and displayed consistent DNA damages. The administration of L-OHP prior to irradiation revealed higher cross-links formation in the Colo320 cell line as compared to HT-29. These results are in agreement with the cytotoxicity findings which suggested that Colo320 cell line is more sensitive to L-OHP than HT-29. An intriguing fact was the lack of response to the same dose of ionizing radiation (2 Gy) in the chemoresistant cell groups. Moreover, the higher dose of irradiation (4 Gy) caused DNA lesions only in HT-29R cells, while Colo320R remained unresponsive. These results suggest that acquired resistance to a chemotherapeutic agent could have activated general resistance pathways that impart resistance to multiple agents, including irradiation. Another potential explanation for the lack of the response to irradiation could be the presence of some free radical scavengers that might have contributed to the resistance to irradiation. The redox homeostasis of the cells was previously implicated in chemo-resistance. It is commonly accepted that the sensitivity of tumor cells to L-OHP, can be influenced by gluthatione and gluthatione related enzymes . Thus, gluthatione S-tranferase (GST) catalyses the conjugation of gluthatione to genotoxic compounds, preventing DNA-damage and adduct formation . Gluthatione system limits the cytotoxic activity of oxaliplatin by modifying the production of cellular reactive species (ROS) . If ROS production decreases, also the cytotoxicity of the chemotherapy agent decreases. Additionally, ROS act also as critical mediators in ionizing radiation-induced cell killing . Moreover, some subsets of tumor cells were shown to possess lower levels of ROS and enhanced ROS defense, which contributed to their radioresistance . Therefore, a link between redox homeostasis and cells’ resistance to chemo/radiotherapy is conceivable. Given that this subset of cells, with increased radioresistance was ascertained as cancer stem cells (CSC), we can presume that another cause of the cells’ peculiar response to ionizing irradiation could be the phenotypic characteristics of the cells. Resistant cells in our study have switched their morphology, as our microscopy findings suggest it. It is also possible that this population might have adopted stem-like characters. Literature data sustain that cancer stem cells (CSC), expressing CD133+ marker manifest resistance to irradiation with 2, 5 Gy unlike CD133- non stem-cells . This different behavior was explained by the activation of the DNA damage checkpoint more efficiently in CSC than in tumor cells without stem cells properties, due to the activation of Chk1 and Chk2 checkpoint kinases . We also demonstrated in the present study the phosphorylation of Chk2 and p53, due to activation of ataxia-telangiectasia mutated gene ATM, this being activated by an apoptosis caspase activation inhibitor (AVEN). Activation of AVEN being evident in Colo320R cell line and not in HT-29R may be a possible explanation for the differential behavior of these cells as compared to HT-29R (i.e. HT-29R responded to 4 Gy, unlike Colo320, which remained resistant). The significantly different gene expression profiles of the tested cell lines, sustained also by cytotoxicity suggest that these cells have completely different genomic patterns, which may explain also their different behavior towards gamma irradiation. The fact that in Colo320R cell line we still could detect cross-links confirms once more the higher sensitivity of this cell line to L-OHP as compared to HT-29. In a previous study we observed changes in the comet tails length according to the degree of the treatment with L-OHP, indicating the presence or absence of the DNA cross-links induced by L-OHP on the tested cell lines . In another study we induced standardized DNA strand breaks via ionizing irradiation and compared the cross-links formation capacity of L-OHP vs CDDP in 3 CC cell lines (Colo320, Caco-2 and HT-29) . Our results showed higher cytotoxicity and cross-links formation for L-OHP vs CDDP in the tested cell lines, in spite of lower cellular uptake. In the present study, we determined indirectly the cross-links formation on the irradiation-damaged DNA and demonstrated that L-OHP resistant cells form fewer cross-links than the parental ones, in accordance with the sensitivity of the tested cells to the cytotoxic drug. Similar studies were performed by others treating with CDDP the healthy human leukocytes  and ovarian carcinoma cells  and using L-OHP treatment on lymphocytes and lung carcinoma cell lines . These studies used UV-C, methyl methanesulfonate and ionizing radiations, respectively, as DNA strand-breaking agents and detected DNA cross-links using an in vitro measure of the cells’ chemosensitivity to the tested compounds.
Our microarray data were in agreement with the morphology, cytotxicity and DNA lesions findings showing that the prolonged treatment with L-OHP triggered different patterns in the transcriptional profiles of the two tested cell lines. To our knowledge, there are no similar studies to highlight the differences between the molecular patterns of these two resistant cell lines however there are genomics studies that evaluated the resistance to treatment either in Colo320 or HT-29 . Considering the common origin of these cell lines (adenocarcinomas) and the mechanism of action of L-OHP which blocks DNA replication and transcription through the formation of intra-strand DNA adducts, we would expect at least to some extent, similar molecular and cellular behavior. Surprisingly, our microarray data have revealed only a common core set of 36 genes modulated more than 1.5-fold in both cell lines (p < 0.05) of which just 27 genes exhibited similar profiles (Table 2). These results could be partly explained by the distinct morphology (suspension vs. adherent) and by the intrinsic differences of the two cell lines which emphasize the complexity of the processes that control the resistance acquirement to this cytostatic drug. Our data revealed that L-OHP modulates genes involved in the regulation of some critical mechanisms including DNA replication, cell death and survival, cellular growth and proliferation, cellular movement and cell-to-cell signaling and interaction.
The microarray analysis showed upregulation of keratin 18 (KRT18) and protein tyrosine phosphatase receptor type O (PTPRO), both being involved in apoptosis. The microarray results validated by qRT-PCR confirmed a significant overexpression of these genes in both HT-29R and Colo320R (Table 6). KRT18 was previously identified as being upregulated in colon carcinoma cells . Increased KRT18 expression has been reported to inhibit cytokine-induced death of cervical cancer cells  but there are no evidences about the role of KRT18 in L-OHP-induced resistance in CC. PTPRO is a member of family of receptor-type protein tyrosine phosphatases with multiple tissue-specific functions including inhibition of cell proliferation and promoting of apoptosis. PTPRO was identified as a target gene of Wnt/β-catenin signaling  and a novel regulator of ERBB2 signaling for mammary epithelial transformation . Ramaswamy et al. observed increased expression of PTPRO in breast cancer following the treatment with tamoxifen . In CC there are no studies describing the implication of PTPRO in drug resistance, but this gene was found to be methylated in colon tumors .
The core set of common DE genes also included some members of interferon - inducible transmembrane gene (IFTIM), whose transmembrane proteins are involved in the homotypic cell adhesion functions of interferon (IFN) . We identified significant upregulation of IFITM3, IFITM4P and IFIH1 in HT29R and downregulation of these genes in Colo320R (Table 2, Class C). The overexpression of IFTIM3 is related to an increased proliferation and metastasis of human colon cancer cells. Andreu et al. identified high endogenous levels of IFITM3 in HT29 cells with APC mutated gene . The authors demonstrated that induction of wild-type APC causes a reduction on IFTIM3 genes within 24 hours. In another study, Ghaleb et al. demonstrated that IFITM3 transcription is dependent on activation of Wnt/β-catenin signaling, in intestinal epithelium . This study appears to be in concordance with our results. Analyzing the canonical pathways for both cell lines we noticed an increased activity for Wnt/β-catenin signaling in HT29R but not in Colo320R (Tables 3, 4). These findings support the morphological observations which suggest an epithelial-to-mesenchymal transition in HT-29R cells.
N-myc downstream regulated 1 (NDRG1) gene had a conflicting expression in the two cell lines, being overexpressed in Colo320R and underexpressed in HT-29R (Table 2, Class D). qRT-PCR confirmed upregulation of NDRG1 in Colo320R and downregulation in HT-29R as a result of prolonged treatment with L-OHP (Table 6). The protein encoded by NDRG1 is implicated in p53-mediated caspase activation and apoptosis. Strzelczyk et al. showed correlation between low levels of NDRG1 gene expression and poor prognosis and survival for patients with CC . These results could suggest that lower level of NDRG1 in HT29R than in Colo320R could be related to a more resistant phenotype.
In response to treatment with cytostatic drugs, cells undergo apoptosis according to the drug-induced DNA damage and the cells’ capacity of DNA repair and survival. In Colo320R, the apoptotic process was mediated by genes involved in caspase modulation and cell cycle regulation. Our results showed that apoptosis caspase activation inhibitor (AVEN), Galectin-3 (LGALS3) and nucleolar protein 3 (NOL3) were overexpressed in this cell line. AVEN represents an activator for ataxia-telangiectasia mutated gene (ATM) which has an important role in the repair of DNA breaks . Cell-cycle arrest induced by DNA damage depends on activation of ATM protein kinase, which phosphorylates cell-cycle effectors such as CHEK2 and p53 in order to inhibit cell-cycle progression. LGALS3 and NOL3 are known as downregulators of the enzyme activities of caspase 2, caspase 8 and tumor protein p53. LGALS3 is involved in the resistance of human colon cancers by blocking the death-inducing signaling complex (DISC) formation and recruitment of the apoptosis-initiating protease, procaspase-8 . Conversely, the increased expression of NOL3 reduced the TRAIL-induced apoptosis in SW480 CC cells . We observed an inhibition of cyclin-dependent kinase inhibitor 2A (CDKN2A) and WNT inhibitory factor 1 (WIF1) tumor suppressor genes in Colo320R following the L-OHP treatment. The lack of function of these genes was associated with tumor cell progression [42, 43].
In addition to the inhibition of apoptosis, our results pointed out to the activation of the mechanisms involved in promoting cell survival and tumor progression in Colo320R. We observed overexpression of insulin-like growth factor 2 (IGF2), mitogen-activated protein kinase kinase kinase 6 (MAP3K6), FBJ murine osteosarcoma viral oncogene homolog (FOS), inhibitor of DNA binding 1 genes (ID1), involved especially in signal transduction on MAP kinase cascade. Our data are in agreement with the literature data concerning the role and implication of these molecules in CC [44–46]. We also noticed downregulation of the peroxisome proliferator-activated receptor γ (PPARG) in Colo320R. Inhibition of PPARG promotes the cell proliferation and leads to the expression of c-myc and cyclin D1 genes as well as of the beta-catenin protein in the colon epithelium .
As we underlined above, the genes evidenced with microarray analysis were quite different between the two tested cell lines, therefore we expected apoptosis in HT-29R to be modulated by a different set of genes compared to those identified for Colo320R. The ineffective induction of apoptosis in HT-29R cell line was mediated by genes involved in Bcl-2 modulation. One of these genes, serum/glucocorticoid regulated kinase 1 (SGK1) is involved in phosphorylation and inactivation of the apoptotic transcription factor forkhead box O3 (FKHRL1), that upregulates death receptor components such as tumor necrosis factor receptor superfamily, member 10 b (TNFRSF10B) and proapoptotic Bcl-2 proteins such as Pim . Another member of Bcl-2 family upregulated in HT-29R cell line, XIAP associated factor 1 (XAF1), has an important role in modulation of apoptosis in tumor cells by inhibiting the caspase-3 activity .
We also noticed that the L-OHP resistance in HT-29R was promoted by the overexpression of some important modulators involved in cell proliferation such as: midkine (MDK), cysteine-rich, angiogenic inducer 61 (CYR61), proliferating cell nuclear antigen (PCNA), transforming growth factor beta 1 (TGFB1), spleen tyrosine kinase (SYK) and prostaglandin-endoperoxide synthase 2 (PTGS2). Upregulation of MDK was correlated with tumor progression in oral squamous cell carcinoma , but to our knowledge, MDK was not found to be expressed in CC. CYR61 has multiple roles in tumor growth, adhesion and migration, its role as positive growth-regulator in CC being previously described . PTGS2 and PCNA represent two important molecules for the progression of CC and treatment strategy  and increased levels of PTGS2 were associated with enhanced tumor cell proliferation and tumorigenesis .
Subsequently, we paid attention to upstream regulators in attempt to explain the different phenotypes induced by L-OHP in the two tested CC cell lines. Although statistically significant results were not found in Colo320R, we identified 10 upstream regulators in HT29R, eight out of them with “activated” and two with “inhibited” predictions (Table 5). Among the top predicted regulators in HT29R dataset, the signal transducer and activator of transcription 1 (STAT1) appears to be significantly activated (z-score = 3.592) by L-OHP treatment, its targets being highly enriched in the data (p = 1.47E-08). STAT1 represents an important activator of transcription in CC . The activation of STAT1 in HT-29R was associated with an increased transcriptional activity in a large number of associated interferon-inducible genes (IFITM1, IFIT3, IFIT1, IFI6, IFI27) (Figure 5). We observed that STAT1 is in turn activated by IFNE, NFATC2 and TGM2. The upregulation of STAT1 mediated by interferon epsilon (IFNE) was described on cervical cancer cells , but to our knowledge, there are no studies describing the role of IFNE in CC. Nuclear factor of activated T cells (NFATC2) is a transcription factor with an important role in the transcriptional regulation of the immune response. Its role in CC in promoting carcinoma migration and invasion was previously demonstrated in in vitro and in vivo studies [56, 57]. Our data confirm the implication of NAFTC2 in HT-29R CC cell line as an important upstream regulator (z-score = 3.019, p = 3.84E-05). Transglutaminase 2 (TGM2) is an enzyme involved in cell proliferation, differentiation and apoptosis and could mediate chemoresistance in cancer cells . TGM2 was proposed by Miyoshi et al. as a predictive marker for prognosis and therapeutic target in CC . In our study, TGM2 was the most notably activated upstream regulator identified in HT29R cell lines (z-score = 4.762, p = 8.9E-09) (Table 5). Recent studies have shown that TGM2 promotes drug resistance and invasion by inducing a stem cell-like phenotype in ovarian and breast cancer [60, 61]. Our results showed that overexpression of TGM2 and NAFTC2 induced TGFB1 activation. High level of TGFB1 expression in HT-29R obtained by microarray was confirmed by qRT-PCR (Table 6). TGFB1 is a well-characterized inducer of EMT in ovarian cancer and human squamous cell carcinoma cells, resulting in increased cell migration and invasion [62, 63]. Considering all these findings, upregulation of TGM2, NAFTC2 and TGFB1 in HT29R but not in Colo320R could explain the induction of EMT and acquiring a more resistant phenotype in HT29R.
Another target of TGM2, poly (ADP-ribose) polymerase family member 9 (PARP9) upstream modulator was significantly activated (z-score = 2.433, p = 1.24E-06) in HT-29R and acts as regulator of STAT1. PARP9 was identified as overexpressed in chemoresistant, diffuse large B-cell lymphomas (DLBCLs) , but there is no data concerning the implication of PARP9 in CC.
In our study CC cells adopted several cellular and molecular alterations during the prolonged treatment with L-OHP which led to resistance to this drug. L-OHP resistant cells displayed altered morphologies, higher invasiveness and metastatic capacities, lower cytotoxicities, formed fewer Pt-DNA cross-links and had different gene expression profiles as compared to the sensitive ones. More disrupted functions and pathways were identified in HT-29R than in Colo320R cells, involving genes responsible for apoptosis inhibition, cellular proliferation and epithelial-to-mesenchymal transition. These findings, in agreement with the morphological and cytotoxicity results and the main upstream regulators identified for HT-29R, but not for Colo320R could explain the more resistant phenotype in HT-29R than in Colo320R cell line.
Therefore, we can conclude that prolonged therapy with L-OHP induces different cellular and molecular chemoresistance patterns in CC cells of identical origins (adenocarcinomas). The set of genes modulated by L-OHP and the upstream regulators revealed in our study explain the diverse behavior of the cancer cells to prolonged therapy with L-OHP, moreover could help us to identify some potential means to reverse chemoresistance and consequently to improve the outcome of therapy in CC.
Cell lines and cultures
Colo320 and HT-29 human CC cell lines were obtained from the European Collection of Cell Cultures (ECACC). Colo320 was cultured in RPMI-1640 and HT-29 in McCoy’s 5A Modified Medium, both supplemented with Fetal Calf Serum 10%, L-glutamin and penicillin-streptomycin (Sigma-Aldrich, St. Louis, MO). Experiments were done at 70–80% cell confluence and confirmed in at least three independent experiments.
Development of L-OHP-resistant cell lines
Resistance to L-OHP (Actavis, Bucharest, Romania) was induced by exposing the cells to increasing concentrations of the drug. The initial dose was 0.01 μg/ml and the final concentration (0.87 μg/ml) corresponded to the clinically relevant plasma concentration of L-OHP (2 μmol/l) . The resistant variant of Colo320 (Colo320R) was obtained and described previously . For the HT-29 cell line we used the same procedure, sequentially increasing concentrations of the drug (with 0.05 μg/ml) being added to the cell culture at every second passage. The surviving cells were grown and propagated every 4–5 days. For both cell lines two groups were considered for investigations: parental (Colo320 and HT-29) and cells with induced chemoresistance (Colo320R and HT-29R). All groups were cultivated in specific media, for the chemoresistant cells the culture media being supplemented with L-OHP 0.87 μg/ml at every second passage.
The morphological analysis of the cells was performed with a light microscope (Carl Zeiss MicroImaging GmbH, Gottingen, Germany), with digital photographic capability. The microscopic images of cells were compared with respect to morphological characteristics (shape, polarity, intercellular distances and presence of pseudopodia).
Evaluation of L-OHP cytotoxicity
The MTS/PMS (Promega Corporation, Madison, WI, USA) colorimetric cell proliferation assay was used for the Colo320 and Colo320R suspension-type and MTT test (Sigma-Aldrich, St. Louis, MO, USA) for HT-29 and HT-29R adherent cells, in order to simplify the solubilization procedures for the suspension-type cell lines . Briefly, cells were seeded in triplicate in 96-well flat-bottom plates, at a cell population density of 15x103 and 2x104 for MTS/PMS and MTT assays, respectively. After 24 h, variable concentrations of L-OHP (0.001-300 μg/ml) were added and the cells were incubated for additional 24 h. Absorbances were recorded with an ELISA plate reader, at 490 nm wavelength (Tecan Sunrise, Grödig/Salzburg, Austria). The half maximal inhibitory concentration (IC50) values were calculated as the concentrations corresponding to a 50% reduction of the cellular growth.
Evaluation of cross-links formation
Platinum-DNA cross-links induced by L-OHP on the selected cell lines were determined indirectly with CA. Single strand-breaks were induced in the L-OHP-treated DNA via ionizing radiation as secondary genotoxic agent. We evidenced the L-OHP-DNA cross-links by quantifying the reduction in the single-strand breaks provoked by the ionizing radiations subsequent to L-OHP treatment. Briefly, the selected cell groups were seeded in triplicate in 24-well plates at a cell population density of 4x105. After 24 h incubation, cells were either irradiated with doses of 2 Gy or 4 Gy gamma radiations using a Co60 source (Theratron 1000, Theratronics, Inc., Ottawa, Ontario, Canada), or exposed first to L-OHP (50 or 100 μg/ml) for 2 h and afterward gamma-irradiated. Dose rate of the applied radiation source was 1, 98 Gy/min. The cells were transported from the site of the irradiation to the CA laboratory on ice, in order to avoid DNA repair process. CA was performed under alkaline conditions according to Tice’s protocol . The images of the cell’s nuclei were evaluated for DNA migration (the level of strand breaks) using a fluorescence microscope (Nikon Eclipse E600, Tokyo, Japan), equipped with 510–560 nm excitation and 590 nm barrier filters. 200 cells from each group were scored and assigned to different lesion classes using the Collins’ classification method . For each group a lesion score (LS) was calculated, based on the number of the cells assigned to each comet classes. LS was expressed as arbitrary units (AU) and the detailed description of the calculation of LS was made elsewhere .
Microarray expression profiling
The gene expression profiling of the parental (Colo320 and HT-29) and L-OHP-resistant (Colo320R and HT-29R) cell lines was assessed with 4x44k Whole Human Genome Oligo Microarray G4112F slides (Agilent Technologies, Santa Clara, CA, USA) using two-color design. Three biological replicates of each cell lines, both parental and resistant, were tested in order to identify genes implicated in the L-OHP-resistance acquiring process. Total RNA was isolated using TriReagent (Sigma-Aldrich, St. Louis, MO, USA) and purified with RNeasy Mini Kit (Qiagen, Hilden, Germany). The quality of RNA was assessed with Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). All samples with an RNA integrity number (RIN) greater than 9 and an rRNA 28S/18S ratio higher than 1.8 were considered suitable for further analysis. The microarray probes were synthesized from 200 ng of total RNA and labelled with fluorochromes Cy3 and Cy5 using Low Input Quick Amp Labeling Kit (Agilent Technologies, Santa Clara, CA, USA), according to the manufacturer’s instructions. Technical replicates of each sample were used for hybridisation control (dye-swap design), each sample being labelled with Cy3 and Cy5 in independent experiments. In each array, RNA extracted from L-OHP-resistant cells was compared with RNA extracted from parental cells. Subsequently, all slides were scanned with Agilent Technologies scanner G2505BUS45102867 and quantification of microarray images were done with Feature Extraction software v. 10.5.3 (Agilent Technologies, Santa Clara, CA, USA).
Microarray data analysis was performed in R statistical programming language . Background and foreground intensity ratios were computed taking log2 ratios of intensities for red (R) and green (G) fluorescence channels (M values). No background subtraction was applied due to the weak correlation between background and foreground intensity ratios (ρ < 0.08). Within-array normalization was carried out using Loess regression. Data were further subjected to between-array normalization by quantile method. Median M values (log2(R/G) values) for duplicate spots were computed and used in class comparison analysis to identify changes in gene expression profiles involved in L-OHP-resistance acquiring process of Colo320 and HT-29 cell lines.
Quantitative real-time PCR (qRT-PCR)
The expression levels of the genes selected by microarray were re-evaluated by qRT-PCR using Light Cycler 480 (Roche Applied Science, Penzberg, Germany) with primers (1 μM) and Universal Probe Library (UPL) probes (0.2 μM). In silico design of UPL probes and primers were obtained from Roche Applied Science Software as follow: PTPRO: F-ctatggagacatcactgtggaga, R-tcctgcatctcgtcagca (UPL#6); KRT18: F-tgatgacaccaatatcacacga, R-ctgggcttgtaggcctttta (UPL#63); NDRG1: F-gggtgcagaagggactagg, R-tgctcctggacatcaaactct (UPL#22); ID1:F-gctgctctacgacatgaacg, R-ctcaccttgcggttctgg (UPL#22); WIF1: F-ccagggagacctctgttcaa, R-ttgggttcatggcaggtt (UPL#76); AVEN F-ggtggtccaagaggaagaagt, R-gaaatcatgctgtccaacca (UPL#22); TGFB1 F-gcagcacgtggagctgta, R-cagccggttgctgaggta (UPL#72); MDK: F-ctcttagcggatgcagcac, R- ccgcccttcttcaccttatc (UPL#63); CYR61: F-aagaaacccggatttgtgag, R-gctgcatttcttgcccttt (UPL#66); 18 s rRNA: F-gcaattattccccatgaacg, R-gggacttaatcaacgcacgc (UPL#48). Each reaction was performed in 5 μl of 1:10 (v/v) dilution of the first cDNA strand synthesized with First Strand cDNA Synthesis Kit (Roche Applied Science, Germany) from 1 μg of the total RNA. The cDNA was then amplified with the Light Cycler Taqman Master Kit (Roche Applied Science, Germany) in a final volume of 20 μl. Thermal cycle conditions included 10 minutes at 95°C for enzyme activation followed by 40 cycles of 15 seconds at 95°C, 20 seconds at 55°C and 1 second at 72°C for the amplification step and 30 seconds at 40°C for the cooling step. The fluorescent signals of UPL probes were used to calculate the cycle thresholds. The target genes were normalized to 18S rRNA housekeeping gene and quantified using the comparative threshold cycle (2-ddct) method described by Livak and Schmittgen .
Statistical processing of the cytotoxicity, CA and qRT-PCR results were done using GraphPad Prism software program, version 5.0 (San Diego, CA, USA). Statistical comparison between groups were made by one-way Anova and Bonferroni posttest and by unpaired two-tailed t test for qRT-PCR data (p < 0.05). For the microarray experiment, the correlation between background and foreground intensity ratios (M values) was assessed in R using Spearman’s rank correlation test. Differentially expressed (DE) genes between resistant and parental cell lines were selected with Limma package/R by fitting a linear model to the expression data for each gene and using empirical Bayes methods to moderate the standard errors across genes . A gene was considered differentially expressed if M value was lower than -0.58 or greater than 0.58 (at least 1.5 -fold down- or up-regulation in resistant versus parental cells) and p value adjusted for multiple testing < 0.05 (Benjamini and Hochberg method). Pearson correlation between microarray and qRT-PCR results were performed in GraphPad Prism software program, version 5.0 (San Diego, CA, USA).
Functional profiling was performed using Ingenuity Pathway Analysis (IPA) software (Ingenuity Systems, Redwood City, California) . Accession numbers of DE genes associated with M values were uploaded into the software. Using information stored in the Ingenuity Knowledge Base (IKB), genes were mapped to genetic networks, molecular functions and canonical pathways. The significance of the association between the genes and the molecular functions and the canonical pathways was determined by Fischer’s exact test (p < 0.05). IPA Upstream Regulator Analysis was used to identify key molecules (upstream regulators) which can affect the expression of their target genes and can regulate each other. To predict the activation state of the upstream regulators (“activated” or “inhibited”), a z-score was computed for each of them. The terms “activated” or “inhibited” does not necessary mean that the regulator is literally activated respectively inhibited. An “activated” upstream regulator indicates a molecule expected to be more active in the resistant cell lines than in the parental ones. A p-value less than 0.01 and a z-score greater than 2 (prediction state: “activated”) or smaller than -2 (prediction state: “inhibited”) were considered significant.
PV: conceived, designed and coordinated the study; induced chemoresistance in the tested cell lines; performed the morphology analysis; EFF: performed comet assay; MPS: read and interpreted the comet assay results; drafted the manuscript; IB: participated in the design of the study and performed statistical analysis for the cellular studies; CT: performed cytotoxicity assays; prepared and treated the cells for further studies; BV: irradiated the cells with gamma radiations. LB: performed statistical and bioinformatic analysis of microarray data; IB-N: carried out the RT-PCR study; OB: carried out the microarray study; participated in the study’s design and drafted the manuscript. All authors read and approved the final manuscript.
This work was supported by the National Council of Scientific Research in Higher Education, Romania, Exploratory Research Projects: PN-II-ID-PCE-2011-3-1057 and PN-II-ID-PCE-2011-3-0753.
- Dahan L, Sadok A, Formento JL, Seitz JF, Kovacic H: Modulation of cellular redox state underlies antagonism between oxaliplatin and cetuximab in human colorectal cancer cell lines. Br J Pharmacol. 2009, 158 (2): 610-620.PubMed CentralView ArticlePubMedGoogle Scholar
- Walko CM, Lindley C: Capecitabine: a review. Clin Ther. 2005, 27 (1): 23-44.View ArticlePubMedGoogle Scholar
- Kelland L: The resurgence of platinum-based cancer chemotherapy. Nat Rev Cancer. 2007, 7 (8): 573-584.View ArticlePubMedGoogle Scholar
- Woynarowski JM, Faivre S, Herzig MC, Arnett B, Chapman WG, Trevino AV, Raymond E, Chaney SG, Vaisman A, Varchenko M, Juniewicz PE: Oxaliplatin-induced damage of cellular DNA. Mol Pharmacol. 2000, 58 (5): 920-927.PubMedGoogle Scholar
- Moggs JG, Szymkowski DE, Yamada M, Karran P, Wood RD: Differential human nucleotide excision repair of paired and mispaired cisplatin-DNA adducts. Nucleic Acids Res. 1997, 25 (3): 480-491.PubMed CentralView ArticlePubMedGoogle Scholar
- Misset JL, Bleiberg H, Sutherland W, Bekradda M, Cvitkovic E: Oxaliplatin clinical activity: a review. Crit Rev Oncol Hematol. 2000, 35 (2): 75-93.View ArticlePubMedGoogle Scholar
- Dallas NA, Xia L, Fan F, Gray MJ, Gaur P, 2nd van Buren G, Samuel S, Kim MP, Lim SJ, Ellis LM: Chemoresistant colorectal cancer cells, the cancer stem cell phenotype, and increased sensitivity to insulin-like growth factor-I receptor inhibition. Cancer Res. 2009, 69 (5): 1951-1957.PubMed CentralView ArticlePubMedGoogle Scholar
- Samimi G, Manorek G, Castel R, Breaux JK, Cheng TC, Berry CC, Los G, Howell SB: cDNA microarray-based identification of genes and pathways associated with oxaliplatin resistance. Cancer Chemother Pharmacol. 2005, 55 (1): 1-11.View ArticlePubMedGoogle Scholar
- Desoize B, Madoulet C: Particular aspects of platinum compounds used at present in cancer treatment. Crit Rev Oncol Hematol. 2002, 42 (3): 317-325.View ArticlePubMedGoogle Scholar
- Safaei R: Role of copper transporters in the uptake and efflux of platinum containing drugs. Cancer Lett. 2006, 234: 34-39.View ArticlePubMedGoogle Scholar
- Virag P, Brie I, Fischer-Fodor E, Perde-Schrepler M, Tatomir C, Balacescu O, Irimie A, Postescu ID: Assessment of cytotoxicity, apoptosis and DNA damages in Colo320 colorectal cancer cells selected for oxaliplatin resistance. Cell Biochem Funct. 2011, 29 (5): 351-355.View ArticlePubMedGoogle Scholar
- Data analysis software. http://www.ingenuity.org/,
- Raymond E, Lawrence R, Izbicka E, Faivre S, Von Hoff DD: Activity of oxaliplatin against human tumor colony-forming units. Clin Cancer Res. 1998, 4 (4): 1021-1029.PubMedGoogle Scholar
- Yang AD, Fan F, Camp ER, van Buren G, Liu W, Somcio R, Gray MJ, Cheng H, Hoff PM, Ellis LM: Chronic oxaliplatin resistance induces epithelial-to-mesenchymal transition in colorectal cancer cell lines. Clin Cancer Res. 2006, 12 (14 Pt 1): 4147-4153.View ArticlePubMedGoogle Scholar
- Thiery JP: Epithelial-mesenchymal transitions in tumour progression. Nat Rev Cancer. 2002, 2 (6): 442-454.View ArticlePubMedGoogle Scholar
- Cavallaro U, Christofori G: Cell adhesion and signalling by cadherins and Ig-CAMs in cancer. Nat Rev Cancer. 2004, 4 (2): 118-132. ReviewView ArticlePubMedGoogle Scholar
- Virag P, Perde-Schrepler M, Fischer-Fodor E, Tatomir C, Dorneanu SA, Cernea VI, Irimie A: Superior cytotoxicity and DNA cross-link induction by oxaliplatin versus cisplatin at lower cellular uptake in colorectal cancer cell lines. Anti-Cancer Drugs. 2012, 23 (10): 1032-1038.View ArticlePubMedGoogle Scholar
- El-Akawi Z, Abu-hadid M, Perez R, Glaw J, Zdanowicz J, Creaven PJ, Pendyala L: Altered glutathione metabolism in oxaliplatin resistant ovarian carcinoma cells. Cancer Lett. 1996, 105: 5-14.View ArticlePubMedGoogle Scholar
- Watson MA, Stewart RK, Smith GB, Massey TE, Bell DA: Human glutathione S-transferase P1 polymorphisms: relantionship to lung tissue enzyme activity and population frequency distribution. Carcinogenesis. 1998, 19: 275-280.View ArticlePubMedGoogle Scholar
- Laurent A, Nicco C, Chéreau C, Goulvestre C, Alexandre J, Alves A, Lévy E, Goldwasser F, Panis Y, Soubrane O, Weill B, Batteux F: Controlling tumor growth by modulating endogenous production of reactive oxygen species. Cancer Res. 2005, 65: 948-956.PubMedGoogle Scholar
- Diehn M, Cho RW, Lobo NA, Kalisky T, Dorie MJ, Kulp AN, Qian D, Lam JS, Ailles LE, Wong M, Joshua B, Kaplan MJ, Wapnir I, Dirbas FM, Somlo G, Garberoglio G, Paz B, Shen J, Lau SK, Quake SR, Brown JM, Weissman IL, Clarke MF: Association of reactive oxygen species levels and radioresistance in cancer stem cells. Nature. 2009, 458 (7239): 780-783. 10.1038/nature07733.PubMed CentralView ArticlePubMedGoogle Scholar
- Beier D, Schulz JB, Beier CP: Chemoresistance of glioblastoma cancer stem cells--much more complex than expected. Mol Cancer. 2011, 10: 128-10.1186/1476-4598-10-128.PubMed CentralView ArticlePubMedGoogle Scholar
- McCord AM, Jamal M, Williams ES, Camphausen K, Tofilon PJ: CD 133+ glioblastoma stem-like cells are radiosensitive with a defective DNA damage response compared with established cell lines. Mol Cancer Res. 2009, 15: 5145-5153.Google Scholar
- Bao S, Wu Q, McLendon RE, Hao Y, Shi Q, Hjelmeland AB, Dewhirst MW, Bigner DD, Rich JN: Glioma stem cells promote radioresistance by preferential activation of DNA damage response. Nature. 2006, 444: 756-760.View ArticlePubMedGoogle Scholar
- Hovhannisyan GG, Haroutunyan TS, Arutyunyan RM: Evaluation of cisplatin-DNA crosslinks formation with UV-C application by the alkaline Comet-assay. Exp Oncol. 2004, 26 (3): 240-242.PubMedGoogle Scholar
- Unger FT, Klasen HA, Tchartchian G, de Wilde RL, Witte I: DNA damage induced by cis- and carboplatin as indicator for in vitro sensitivity of ovarian carcinoma cells. BMC Cancer. 2009, 10 (9): 359-View ArticleGoogle Scholar
- Almeida GM, Duarte TL, Steward WP, Jones GDD: Detection of oxaliplatin-induced DNA crosslinks in vitro and in cancer patients using the alkaline comet assay. DNA Repair. 2006, 5 (2): 219-225.View ArticlePubMedGoogle Scholar
- Martinez-Cardús A, Martinez-Balibrea E, Bandrés E, Malumbres R, Ginés A, Manzano JL, Taron M, Garcia-Foncillas J, Abad A: Pharmacogenomic approach for the identification of novel determinants of acquired resistance to oxaliplatin in colorectal cancer. Mol Cancer Ther. 2009, 8 (1): 194-202.View ArticlePubMedGoogle Scholar
- Delouis C, Prochasson P, Laithier M, Brison O: Use of adenoviral E1A protein to analyze K18 promoter deregulation in colon carcinoma cells discloses a role for CtBP1 and BRCA1. BMC Mol Biol. 2005, 6: 8-PubMed CentralView ArticlePubMedGoogle Scholar
- Sullivan BT, Cherry JA, Sakamoto H, Henkes LE, Townson DH, Rueda BR: Cytokeratin 18 expression inhibits cytokine-induced death of cervical cancer cells. Int J Gynecol Cancer. 2010, 20 (9): 1474-1481.PubMedGoogle Scholar
- Kim M, Kim H, Jho EH: Identification of ptpro as a novel target gene of Wnt signaling and its potential role as a receptor for Wnt. FEBS Lett. 2010, 584 (18): 3923-3928.View ArticlePubMedGoogle Scholar
- Yu M, Lin G, Arshadi N, Kalatskaya I, Xue B, Haider S, Nguyen F, Boutros PC, Elson A, Muthuswamy LB, Tonks NK, Muthuswamy SK: Expression profiling during mammary epithelial cell three-dimensional morphogenesis identifies PTPRO as a novel regulator of morphogenesis and ErbB2-mediated transformation. Mol Cell Biol. 2012, 32 (19): 3913-3924.PubMed CentralView ArticlePubMedGoogle Scholar
- Ramaswamy B, Majumder S, Roy S, Ghoshal K, Kutay H, Datta J, Younes M, Shapiro CL, Motiwala T, Jacob ST: Estrogen-mediated suppression of the gene encoding protein tyrosine phosphatase PTPRO in human breast cancer: mechanism and role in tamoxifen sensitivity. Mol Endocrinol. 2009, 23 (2): 176-187.PubMed CentralView ArticlePubMedGoogle Scholar
- Mori Y, Yin J, Sato F, Sterian A, Simms LA, Selaru FM, Schulmann K, Xu Y, Olaru A, Wang S, Deacu E, Abraham JM, Young J, Leggett BA, Meltzer SJ: Identification of genes uniquely involved in frequent microsatellite instability colon carcinogenesis by expression profiling combined with epigenetic scanning. Cancer Res. 2004, 64 (7): 2434-2438.View ArticlePubMedGoogle Scholar
- Lewin AR, Reid LE, McMahon M, Stark GR, Kerr IM: Molecular analysis of a human interferon-inducible gene family. Eur J Biochem. 1991, 199: 417-423.View ArticlePubMedGoogle Scholar
- Andreu P, Colnot S, Godard C, Laurent-Puig P, Lamarque D, Kahn A, Perret C, Romagnolo B: Identification of the IFITM family as a new molecular marker in human colorectal tumors. Cancer Res. 2006, 66 (4): 1949-1955.View ArticlePubMedGoogle Scholar
- Ghaleb AM, Aggarwal G, Bialkowska AB, Nandan MO, Yang VW: Notch inhibits expression of the Krüppel-like factor 4 tumor suppressor in the intestinal epithelium. Mol Cancer Res. 2008, 6: 1920-1927.PubMed CentralView ArticlePubMedGoogle Scholar
- Strzelczyk B, Szulc A, Rzepko R, Kitowska A, Skokowski J, Szutowicz A, Pawelczyk T: Identification of high-risk stage II colorectal tumors by combined analysis of the NDRG1 gene expression and the depth of tumor invasion. Ann Surg Oncol. 2009, 16 (5): 1287-1294.View ArticlePubMedGoogle Scholar
- Guo JY, Yamada A, Kajino T, Wu JQ, Tang W, Freel CD, Feng J, Chau BN, Wang MZ, Margolis SS, Yoo HY, Wang XF, Dunphy WG, Irusta PM, Hardwick JM, Kornbluth S: Aven-dependent activation of ATM following DNA damage. Curr Biol. 2008, 18 (13): 933-942.PubMed CentralView ArticlePubMedGoogle Scholar
- Mazurek N, Byrd JC, Sun Y, Hafley M, Ramirez K, Burks J, Bresalier R: Cell-surface galectin-3 confers resistance to TRAIL by impeding trafficking of death receptors in metastatic colon adenocarcinoma cells. Cell Death Differ. 2011, 10.1038/cdd.2011.123.Google Scholar
- Ao JE, Kuang LH, Zhou Y, Zhao R, Yang CM: Hypoxia-inducible factor 1 regulated ARC expression mediated hypoxia induced inactivation of the intrinsic death pathway in p53 deficient human colon cancer cells. Biochem Biophys Res Commun. 2012, 420 (4): 913-917.View ArticlePubMedGoogle Scholar
- Wang C, Li N, Liu X, Zheng Y, Cao X: A novel endogenous human CaMKII inhibitory protein suppresses tumor growth by inducing cell cycle arrest via p27 stabilization. J Biol Chem. 2008, 283 (17): 11565-11574.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang L, Li MX, Wang L, Li BK, Chen GH, He LR, Xu L, Yuan YF: Prognostic value of Wnt inhibitory factor-1 expression in hepatocellular carcinoma that is independent of gene methylation. Tumour Biol. 2011, 32 (1): 233-240.View ArticlePubMedGoogle Scholar
- Eto N, Miyagishi M, Inagi R, Fujita T, Nangaku M: Mitogen-activated protein 3 kinase 6 mediates angiogenic and tumorigenic effects via vascular endothelial growth factor expression. Am J Pathol. 2009, 174 (4): 1553-1563.PubMed CentralView ArticlePubMedGoogle Scholar
- Cara A, Rabik M, Dolan E: Molecular mechanisms of resistance and toxicity associated with platinating agents. Cancer Treat Rev. 2007, 33 (1): 9-23.View ArticleGoogle Scholar
- Valdimarsdottir G, Goumans MJ, Rosendahl A, Brugman M, Itoh S, Lebrin F, Sideras P, ten Dijke P: Stimulation of Id1 expression by bone morphogenetic protein is sufficient and necessary for bone morphogenetic protein-induced activation of endothelial cells. Circulation. 2002, 106 (17): 2263-2270.View ArticlePubMedGoogle Scholar
- Fujisawa T, Sugiyama M, Tomimoto A, Wada K, Endo H, Takahashi H, Yoneda K, Yoneda M, Inamori M, Saito S, Terauchi Y, Kadowaki T, Tsuchiya N, Nakagama H, Nakajima A: Inhibition of peroxisome proliferator-activated receptor gamma promotes tumorigenesis through activation of the beta-catenin/T cell factor (TCF) pathway in the mouse intestine. J Pharmacol Sci. 2008, 108 (4): 535-544.View ArticlePubMedGoogle Scholar
- Ravi A, Thompson CB: The survival kinases Akt and Pim as potential pharmacological targets. J Clin Invest. 2005, 115 (10): 2618-2624.View ArticleGoogle Scholar
- Yang L, Cao Z, Yan H, Wood WC: Coexistence of high levels of apoptotic signaling and inhibitor of apoptosis proteins in human tumor cells: implication for cancer specific therapy. Cancer Res. 2003, 63 (20): 6815-6824.PubMedGoogle Scholar
- Ota K, Fujimori H, Ueda M, Jono H, Shinriki S, Ota T, Sueyoshi T, Taura M, Taguma A, Kai H, Shinohara M, Ando Y: Midkine expression is correlated with an adverse prognosis and is downregulated by p53 in oral squamous cell carcinoma. Int J Oncol. 2010, 37 (4): 797-804.PubMedGoogle Scholar
- de Angelis PM, Fjell B, Kravik KL, Haug T, Tunheim SH, Reichelt W, Beigi M, Clausen OP, Galteland E, Stokke T: Molecular characterizations of derivatives of HCT116 colorectal cancer cells that are resistant to the chemotherapeutic agent 5-fluorouracil. Int J Oncol. 2004, 24 (5): 1279-1288.PubMedGoogle Scholar
- Zhang YH, Sheng JQ, Geng HG, Wu ZT, Li AQ, Li SR: Effects of cyclooxygenase-2 and proliferating cell nuclear antigen on the onset and development of familial adenomatous polyposis. Ai Zheng. 2009, 28 (11): 1181-1185.PubMedGoogle Scholar
- Tong X, Yin L, Joshi S, Rosenberg DW, Giardina C: Cyclooxygenase-2 regulation in colon cancer cells: modulation of RNA polymerase II elongation by histone deacetylase inhibitors. J Biol Chem. 2005, 280 (16): 15503-15509.View ArticlePubMedGoogle Scholar
- Slattery ML, Lundgreen A, Kadlubar SA, Bondurant KL, Wolff RK: JAK/STAT/SOCS-signaling pathway and colon and rectal cancer. Mol Carcinog. 2011, 10.1002/mc.21841.Google Scholar
- Matsumiya T, Prescott SM, Stafforini DM: IFN-epsilon mediates TNF-alpha-induced STAT1 phosphorylation and induction of retinoic acid-inducible gene-I in human cervical cancer cells. J Immunol. 2007, 179 (7): 4542-4549.View ArticlePubMedGoogle Scholar
- Duque J, Fresno M, Iñiguez MA: Expression and function of the nuclear factor of activated T cells in colon carcinoma cells: involvement in the regulation of cyclooxygenase-2. J Biol Chem. 2005, 280 (10): 8686-8693.View ArticlePubMedGoogle Scholar
- Gerlach K, Daniel C, Lehr HA, Nikolaev A, Gerlach T, Atreya R, Rose-John S, Neurath MF, Weigmann B: Transcription factor NFATc2 controls the emergence of colon cancer associated with IL-6-dependent colitis. Cancer Res. 2012, 72 (17): 4340-4350.View ArticlePubMedGoogle Scholar
- Verma A, Mehta K: Tissue transglutaminase-mediated chemoresistance in cancer cells. Drug Resist Updat. 2007, 10 (4–5): 144-151.View ArticlePubMedGoogle Scholar
- Miyoshi N, Ishii H, Mimori K, Tanaka F, Hitora T, Tei M, Sekimoto M, Doki Y, Mori M: TGM2 is a novel marker for prognosis and therapeutic target in colorectal cancer. Ann Surg Oncol. 2010, 17 (4): 967-972.View ArticlePubMedGoogle Scholar
- Cao L, Shao M, Schilder J, Guise T, Mohammad KS, Matei D: Tissue transglutaminase links TGF-β, epithelial to mesenchymal transition and a stem cell phenotype in ovarian cancer. Oncogene. 2012, 31 (20): 2521-2534.View ArticlePubMedGoogle Scholar
- Kumar A, Xu J, Brady S, Gao H, Yu D, Reuben J, Mehta K: Tissue transglutaminase promotes drug resistance and invasion by inducing mesenchymal transition in mammary epithelial cells. PLoS One. 2010, 5 (10): e13390-10.1371/journal.pone.0013390. Published online 2010 October 12PubMed CentralView ArticlePubMedGoogle Scholar
- Shao M, Cao L, Shen C, Satpathy M, Chelladurai B, Bigsby RM, Nakshatri H, Matei D: Epithelial-to-mesenchymal transition and ovarian tumor progression induced by tissue transglutaminase. Cancer Res. 2009, 69 (24): 9192-9201. 10.1158/0008-5472.CAN-09-1257. Epub,View ArticlePubMedGoogle Scholar
- Saito D, Kyakumoto S, Chosa N, Ibi M, Takahashi N, Okubo N, Sawada S, Ishisaki A, Kamo M: Transforming growth factor-β1 induces epithelial-mesenchymal transition and integrin α3β1-mediated cell migration of HSC-4 human squamous cell carcinoma cells through Slug. J Biochem. 2013, 153 (3): 303-315.View ArticlePubMedGoogle Scholar
- Juszczynski P, Kutok JL, Li C, Mitra J, Aguiar RC, Shipp MA: BAL1 and BBAP are regulated by a gamma interferon-responsive bidirectional promoter and are overexpressed in diffuse large B-cell lymphomas with a prominent inflammatory infiltrate. Mol Cell Biol. 2006, 26 (14): 5348-5359.PubMed CentralView ArticlePubMedGoogle Scholar
- Cory AH, Owen TC, Barltrop JA, Cory JG: Use of an aqueous soluble tetrazolium/formazan assay for cell growth assays in culture. Cancer Commun. 1991, 3 (7): 207-212.PubMedGoogle Scholar
- Tice RR, Agurell E, Anderson D, Burlinson B, Hartmann A, Kobayashi H, Miyamae Y, Rojas E, Ryu JC, Sasaki YF: Single cell gel/comet assay: guidelines for in vitro and in vivo genetic toxicology testing. Environ Mol Mutagen. 2000, 35: 206-221.View ArticlePubMedGoogle Scholar
- Collins AR, Dobson VL, Dušinská M, Kennedy G, Stĕtina R: The comet assay: what it can really tell us?. Mutat Res. 1997, 375: 183-193.View ArticlePubMedGoogle Scholar
- Scharpf RB, Iacobuzio-Donahue CA, Sneddon JB, Parmigiani G: When should one subtract background fluorescence in 2-color microarrays?. Biostatistics. 2007, 8 (4): 695-707.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408.View ArticlePubMedGoogle Scholar
- Smyth GK: Limma: linear models for microarray data, Bioinformatics and computational biology solutions using R and bioconductor. 2005, New York: SpringerGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.