- Research article
- Open Access
Gene expression analysis of cell death induction by Taurolidine in different malignant cell lines
BMC Cancer volume 10, Article number: 595 (2010)
The anti-infective agent Taurolidine (TRD) has been shown to have cell death inducing properties, but the mechanism of its action is largely unknown. The aim of this study was to identify potential common target genes modulated at the transcriptional level following TRD treatment in tumour cell lines originating from different cancer types.
Five different malignant cell lines (HT29, Chang Liver, HT1080, AsPC-1 and BxPC-3) were incubated with TRD (100 μM, 250 μM and 1000 μM). Proliferation after 8 h and cell viability after 24 h were analyzed by BrdU assay and FACS analysis, respectively. Gene expression analyses were carried out using the Agilent -microarray platform to indentify genes which displayed conjoint regulation following the addition of TRD in all cell lines. Candidate genes were subjected to Ingenuity Pathways Analysis and selected genes were validated by qRT-PCR and Western Blot.
TRD 250 μM caused a significant inhibition of proliferation as well as apoptotic cell death in all cell lines. Among cell death associated genes with the strongest regulation in gene expression, we identified pro-apoptotic transcription factors (EGR1, ATF3) as well as genes involved in the ER stress response (PPP1R15A), in ubiquitination (TRAF6) and mitochondrial apoptotic pathways (PMAIP1).
This is the first conjoint analysis of potential target genes of TRD which was performed simultaneously in different malignant cell lines. The results indicate that TRD might be involved in different signal transduction pathways leading to apoptosis.
Taurolidine (TRD) - a derivate of the aminosulfoacid Taurin - has been clinically used for many years in peritonitis and catheter related blood stream infections due to its anti-microbial and anti-inflammatory properties [1–3]. Recently it has been shown, that TRD also exerts anti-proliferative and anti-neoplastic activity in vitro as well as in vivo [4, 5]. TRD has been reported to inhibit proliferation and to induce programmed cell death in a variety of cell lines derived from malignant tumours e.g. glioblastoma [6, 7], melanoma [8, 9], mesothelioma [10, 11], colon carcinoma [12, 13], squamous cell oesophageal carcinoma  and sarcoma [15, 16]. Recently, favourable pharmacokinetic and safety data for TRD have been reported following systemic application in healthy volunteers  as well as in patients with locally advanced gastric carcinoma and glioblastoma [18–20].
However, cell death inducing mechanisms of TRD remain to be fully elucidated. Both the mitochondrial dependent pathway [7, 10, 21–23] as well as the death receptor associated pathways have been reported for TRD [24, 25, 16, 14]. Since the majority of information about TRD effects is provided from studies with one single cell line, there is a lack of a comprehensive view across several cell lines of different malignancies. So far, only two publications have addressed the changes in gene expression following TRD exposure to malignant cells using cDNA microarray techniques [14, 16]. The aim of this study was therefore, to analyse gene expression by microarray analyses simultaneously in different malignant cell lines - to identify potential TRD target genes which displayed conjoint regulation in all cell lines.
Cell lines and culture conditions
Five different human neoplastic cancer cell lines were used for this experiment: HT29 colon carcinoma (CLS Cell Lines Service, Eppelheim, Germany), Chang Liver (HeLa contaminant, CLS Cell Lines Service, Eppelheim, Germany), HT1080 fibrosarcoma (ATCC - LGC Standards GmbH, Wesel, Germany), AsPC-1 pancreas carcinoma (CLS Cell Lines Service, Eppelheim, Germany) and BxPC-3 pancreas carcinoma (ATCC - LGC Standards GmbH, Wesel, Germany). Chang Liver cells were maintained with Dulbecco's Modified Eagle Medium (DMEM) - Hams's F12, whereas HT1080 cells were cultured in modified Eagle's medium (MEM). The remaining cell lines (HT29, AsPC-1, BxPC-3) were maintained in RPMI 1640 (Biowest, Nuaille, France). All cultures were supplemented with 10% fetal bovine serum, supplemented with penicillin (100 U/ml), streptomycin (100 μg/ml) and 2 mM L-Glutamine (Biowest, Nuaille, France). AsPC-1 and HT1080 cells were further supplemented with 1 mM Sodium Pyruvate. Cells were grown as subconfluent monolayer and cultured in 25 cm2 flasks at 37°C and 5% CO2 in a humidified atmosphere.
TRD (Taurolin®) ultrapure powder (kindly provided by Geistlich Pharma AG, Wolhusen, Switzerland) was dissolved in a Povidon 5% solution (K16 Povidon, generously provided by Geistlich Pharma AG, Wolhusen, Switzerland) and sterile filtered to achieve the respective TRD concentrations. A Povidon 5% solution in equal volume served as a control for TRD treatment.
BrdU proliferation assay
Cells were seeded to a density of 3 × 106 cells/well in 6-well plates (growth area 9.6 cm2/well) and incubated for 18-24 hours under the above mentioned culture conditions to obtain a subconfluent monolayer. Subsequently, cells were washed and incubated for another 2 hours before reagents were added to the culture medium. To examine the dose-response of TRD regarding its anti-proliferative activity, cells were incubated with increasing concentrations of TRD (100 μM, 250 μM, and 1000 μM) and Povidon 5% as control for 8 h and submitted to BrdU proliferation assay (5-bromo-2'-deoxyuridine)-ELISA (Roche Applied Science, Mannheim, Germany) according to the manufacturer's instructions. Our group has recently shown, that the TRD concentration range applied in this experiment (100 μM, 250 μM and 1000 μM) is suitable for examining anti-neoplastic activity in these five cell lines . Based on the incorporation of the thymidine analogue BrdU during DNA synthesis, the amount of newly synthesised DNA and thus of proliferation cells is detected using a microplate absorbance reader Sunrise™ (Tecan trading AG, Switzerland) after applying anti-BrdU conjugated with peroxidase and enhancing a specific substrate reaction. BrdU assays were performed with 8 replicates of three independent experiments with consecutive passages. The incubation time of 8 h has been shwon to be appropriate for the BrdU proliferation assay in previous experiments (data not shown).
Flow Cytometry Analysis
All five cell lines were incubated for 24 h with TRD 250 μM and Povidon 5% as control and submitted to FACS analysis. Our group has recently shown, that TRD 250 μM is sufficient to induce cell death in the respective cell lines . However, for the current study, FACS analysis for this particular concentration (250 μM) was performed in 9-11 independent experiments with 4-5 consecutive passages. In brief, floating cells were collected together with the supernatant and adherent cells were harvested by trypsinisation. Cells were sedimented by centrifugation, resuspended and fixed in 195 μl binding buffer (Bender MedSystems, Vienna, Austria). Cell density in the cell suspension was adjusted to 2 × 103 cells/μl. Subsequently, 5 μl Annexin V-FITC (BD Biosciences, Heidelberg, Germany) was added to the cell suspension followed by gently vortexing and incubation for 10 min at room temperature in the dark. Thereafter, the cell suspension was centrifuged followed by resuspension in 190 μl binding buffer before 10 μl Propidiumiodide (Bender MedSystems, Vienna, Austria) was added. Cells were analyzed immediately using a FACS (fluorescence activated cell sorting) flow cytometer (FACS Calibur BD Biosciences, Heidelberg, Germany) for Annexin V-FITC and PI binding. For each measurement, 20.000 cells were counted. Dot plots and histograms were analyzed by CellQuest Pro software (BD Biosciences, Heidelberg, Germany). Annexin V positive cells were considered apoptotic; Annexin V and PI positive cells were identified as necrotic. Annexin V and PI negative cells were termed viable. All experiments were repeated 9-11 times with 4-5 consecutive passages. The incubation time of 24 h has been shwon to be appropriate for the FACS analysis in previous experiments lines .
cDNA microarray analysis
Total RNA was purified from the cells after incubation with the different substances (Povidon 5%; TRD 100 μM, 250 μM and 1000 μM) for 6 h using the RNeasy KIT from Qiagen (Hilden, Germany), as specified by the manufacturer. RNA integrity was assessed using standard denaturing agarose gel electrophoresis. For microarray analyses, we used the Agilent Array platform employing the manufacturer's standard protocols for sample preparation and microarray hybridization. Briefly, total RNA (500 ng) from each sample was amplified and transcribed into fluorescent cRNA following the manufacturer's Quick Amp Labeling protocol (Version 5.7, Agilent Technologies). Since we aimed at an intensity-based analysis of the two-colour microarrays as suggested by 't Hoen et al. , RNA samples from individual series of analyses comprising of RNAs derived from one specific cell line under the various TRD treatments and control conditions (Povidon 5%) were either Cy-5 (HT29, Chang Liver and HT1080) or Cy-3 (BxPC3, AsPC1) labelled. Labelled samples were hybridized towards the Whole Genome Oligo Array (4 × 44 Kk, product no. G4412T; Agilent Technologies) and following the washing steps the arrays were scanned using the Agilent Scanner G2505B. Agilent's Feature Extraction software 10.7.3.1 (Feature Extraction Protocol GE2_ 105-Dec08) was used to analyze acquired array images. Subsequent data processing was performed using the GeneSpringGX11.01 software package (Agilent Technologies). Following quantile normalization of the raw data, data sets were further reduced by filtering the raw signal intensity value to ≥50 in at least one out of the 20 samples analyzed. A one-way ANOVA model followed by Tukey's HSD (Honestly Significant Difference) test was used to test the hypothesis that there was no difference in expression between the TRD treatment group and the Povidon control group and multiplicity correction (Benjamini-Hochberg) was included to control the false discovery rate (FDR) at 0.05%. In a pair-wise comparison for differentially expressed genes identified between POV and TRD 250 μM treated cell lines by ANOVA a subset of genes was identified that displayed a conjoint regulation in all five cell lines after TRD 250 μM treatment with a mean increase or decrease ≥ 2 fold compared to control treatment (Povidon 5%). This subset of genes was subjected to Ingenuity Pathways Analysis (Ingenuity® Systems, Redwood City, CA, USA). Genes associated with biological functions in Ingenuity's Knowledge Base were considered for the analysis. Right-tailed Fisher's exact test was used to calculate a p-value determining the probability that each biological function assigned to that data set is due to chance alone.
Real-time PCR for microarray data validation
Microarray data validation was performed for the selected candidate genes ATF3, GADD34, TRAF6, PMAIP1 as well as EGR1 for TRD 250 μM and Povidon 5% as control. Total RNA (2 μg) was reverse transcribed using the High Capacity cDNA Archive Kit (Applied Biosystems). RNA isolation was performed from cells harvested after 6 h following treatment. cDNA was synthesized using 2 μg of total RNA, oligo(dT)18 primers and 200 U of MMLV reverse transcriptase (Promega, Mannheim, Germany) following the manufacturer's protocol. qRT-PCR was performed using a SYBR Green I reaction mixture containing 75 mM Tris-HCl (pH 8.8), 20 mM ammonium sulfate, 0.01% (v/v) Tween 20, 2 mM magnesium chloride (all Sigma-Aldrich, Munich, Germany), 1 μl of a 600-fold dilution of SYBR Green I (BioWhittaker, Rockland, ME, USA), 2.5 U Taq polymerase (NEB, Frankfurt a.M., Germany), 0.2 mM dNTP (Promega, Mannheim, Germany) and 0.2 μM of forward and reverse primer (QIAgen, Hilden, Germany) in a final reaction volume of 20 μl. Reactions were run on a DNA Engine Opticon®2 cycler (MJ Research, Waltham, MA, USA). The cycling conditions consisted of 3 min initial denaturation at 94°C and 40 cycles of 94°C for 30 sec, 60°C for 30 sec, 72°C for 30 sec and 80°C for 3 sec. Fluorescence was measured at the last step of each cycle. Melting curves were obtained after each PCR run and showed single PCR products. cDNAs were run in duplicates, including non-RT (without reverse transcriptase) and no-template controls. PCR efficiencies were determined using serial dilutions of a cDNA derived from a cell line. Expression levels for genes of interest and for housekeeping genes were measured in independent PCR runs. Expression ratios were calculated as described by Pfaffl  including the geometric mean expression of the housekeeping genes GAPD and PPIA to normalize the expression data for the gene of interest. Primers used are available through the authors upon request.
To validate the findings of changed gene expression on protein level, Western Blots were performed for selected target genes whenever appropriate antibodies were available. Total protein was purified from the cells after incubation for 8 h with Taurolidine 250 μM and Povidon 5% as control treatment. For this purpose floating cells were collected together with the supernatant and sedimented by centrifugation. After removal of the supernatant, adherent cells were added after harvesting with a cell scrapper and Cell Culture Lysis Reagent (Promega Corporation, Mannheim, Germany). Probes were incubated with a total amount of 50 μl Cell Culture Lysis Reagent (Promega Corporation, Mannheim, Germany) each for 1 h on ice. The cell remnants were sedimented by centrifugation leading to supernatants containing the purified protein. Protein concentration in each sample was determined using a Micro BCA protein assay kit (Pierce, Rockford, IL, USA). Aliquots that contained equal amount of total proteins from each sample were separated in SDS-polyacrylamide gels and electrophoretically transferred onto PVDF membranes (Bio-Rad Laboratories, Munich, Germany), which were subsequently probed for 1 h - 3 h at room temperature with ATF3 (C-19) or GADD34/PPP1R15A (S-20) antibodies (Santa Cruz Biotechnology, Heidelberg, Germany). Blots were further incubated for 1 h with secondary horseradish peroxidase (HRP)-conjugated anti-rabbit IgG (GE Healthcare, Little Chalfont, UK), and then developed using an enhanced chemiluminescence detection system (ECL Plus Western Blotting Detection System, GE Healthcare, Little Chalfont, UK) according to the manufacturers' instructions. After finishing the western blot for the respective target genes, each membrane was stripped with stripping buffer (100 mM NaOH, 2% SDS, 0.5% DTT, 1 h for 55C) and incubated with GAPDH antibody (FL-335) (Santa Cruz Biotechnology, Heidelberg, Germany) as loading control.
Results of FACS-analysis (percentage of viable, apoptotic and necrotic cells) as well as results of BrdU assay (percentage of proliferating cells) are expressed as means ± SEM. Comparison between experimental groups (control, TRD 100 μM, TRD 250 μM and TRD 1000 μM) was performed using one-way ANOVA followed by Tukey's post-hoc test. P-values ≤ 0.05 were considered as statistically significant and indicated in the figures as follows: *** p ≤ 0.001, ** p ≤ 0.01, * p ≤ 0.05.
TRD uniformly inhibits proliferation of all cell lines in a dose dependent manner
As indicated in figure 1, incubation with TRD in increasing concentrations (100 μM, 250 μM and 1000 μM) for 8 hours resulted in an dose dependent reduction of proliferating cells - as measured by BrdU assay. In all cell lines except HT29, even the lowest TRD concentration (100 μM) was capable of inhibiting proliferation leading to values of proliferating cells ranging between 53.4% (± 11.0%) (Chang Liver) and 93.8% (± 8.1%) (BxPC-3) which was significantly lower compared to untreated controls (100%) (Povidon 5%) (figure 1). However, the following TRD concentration of 250 μM significantly inhibited proliferation in all five cell lines. The pronounced inhibition mediated by TRD 250 μM resulted in values of proliferating cells between 12.7% (± 0.56%) for Chang liver and 56.1% (±11.7) for HT29 (figure 1a-b). The maximal dose of 1000 μM led to further inhibition of proliferation in all cell lines. As a result, dose response for cell proliferation was proportional in all cell lines (figure 1).
TRD induces apoptotic cell death in all cell lines at a concentration of 250 μM
As summarized in figure 2, incubation of the five cell lines with TRD 250 μM for 24 hours resulted in a significant reduction of viable cells compared to control treatment with Povidon 5% as evaluated by FACS analysis for Annexin V-FITC and PI. The TRD concentration of 250 μM was choosen in this experiment, since we could previously show that TRD 250 μM was the lowest among several concentrations (100 μM, 250 μM and 1000 μM) that significantly inhibited cell viability in all cell lines . Furthermore, TRD 250 μM was the lowest concentration that significantly inhibited proliferation in all cell lines as measured by BrdU assay (figure 1). As shown in figure 2a, cell viability following incubation with TRD 250 μM was varying between 69.5% (± 3.2%) for HT29 and 35.7% (± 1.7%) for Chang Liver cells. The significant reduction of cell viability by TRD 250 M was paralleled by a significant increase of apoptotic cells in all cell lines (figure 2b). The apoptotic effect was ranging between 20.1% (± 1.5%) for BxPC-3 and 31.0% (± 2.4%) for AsPC-1) (figure 2b). Additionally, there was also a significant increase in necrosis in Chang Liver, HT1080 and BxPC-3 cells (figure 2c). Representative dot plots obtained by FACS analysis for incubation with TRD 250 μM and control treatment (Povidon 5%) in the five respective cell lines are presented in figure 3.
TRD 250 μM leads to conjoint regulation of different genes in all five malignant cell lines
In order to answer the question whether TRD treatment may modulate signalling pathways common in several different tumour types rather than being tumour cell type specific we subjected the above described different cell lines (HT29 (colon), Chang Liver (liver), HT1080 (fibrosarcoma), AsPC-1, BxPC-3 (both pancreas)) to global gene expression analyses using Agilent microarray technology. A comparison of the expression profiles in the TRD treatment (250 μM) vs. control treatment (Povidon 5%) revealed 592 probe sets that were significantly up- or down regulated in all five cell lines with a mean change ≥ 2 fold compared to the control treatment group (additional file 1, table S1). A total of 450 probe sets were significantly down-regulated in all cell lines with a mean decrease of 3.1 fold (95% CI 2.97 - 3.18) from control treatment to TRD 250 μM treatment. The remaining 142 genes were up-regulated with a mean increase of 7.3 fold (95% CI 6.97 - 7.72) from control treatment to TRD 250 μM treatment. From these genes we identified the top high-level functions by Ingenuity global function analysis (p < 0.05). 55 genes out of 592 probe sets were associated with at least one of the three functions "Cell Death", "Cell Growth and Proliferation", and "Cell Cycle" (figure 4). The function "Cell Death" was most frequent and associated with 36 genes, followed by "Cell Cycle (29 genes) and "Cell Growth and Proliferation" (29 genes) (figure 4) (additional file 1, table S1). The distribution and overlapping of target genes within the three functions is illustrated in figure 4. Potential TRD regulated genes common to all five cell lines that increased or decreased ≥ 4 fold and that were identified by Ingenuity Pathways Analysis are further summarized in table 1. There was a striking up-regulation of several transcription factors (TFs) like EGR1, ATF3, FOSB, FOS, SNAI1, NRA2 and HES1. Other groups of genes comprised cell cycle regulators (GADD45A, GADD45B, E2F3 and CDKN2B) as well as genes involved in the ubiquitin pathway (UBE2B, TRAF6), the endoplasmic reticulum response (PPP1R15A), the mitochondrial apoptotic pathway (PMAIP1) or the death receptor pathway (FADD) (table 1). The complete microarray data for all probe sets with the respective normalised values (all cell lines, Povidon 5%, TRD 100 μM, TRD 250 μM, TRD 1000 μM) are provided in additional file 2, table S2.
Confirmation of altered gene expression following TRD 250 μM incubation by qRT-PCR and western blot
Messenger RNA expression of five selected TRD regulated genes was examined by quantitative real-time RT-PCR in all cell lines (figure 5). The gene with the strongest regulation during microarray analysis was EGR1 (Early growth response 1) with a mean fold change (FC) of 45.9 (table 1) which belongs to the EGR zinc-finger family of TFs  and is known to have anti-neoplastic activity . This intensive induction could be confirmed by qRT-PCR, which revealed an up-regulation ranging from 8-fold to 145-fold (figure 5a). The second strongest regulation observed in the microarray analysis was referred to ATF3 (Activating transcription factor 3) (mean FC of 39.9) (table 1). ATF 3 is a member of the CREB family of transcription factors and is activated by various stimuli . This result could also be confirmed by qRT-PCR, which resulted in a 4-fold to 48-fold increase (figure 5b). Two genes of special interest were also validated by qRT-PCR: PPP1R15A (Protein phosphatase 1, regulatory (inhibitory) subunit 15A; syn. GADD34), representing an important protein phosphatase involved in cell death pathways [32, 33] as well as PMAIP1 (Phorbol-12-myristate-13-acetate-induced protein 1; syn. NOXA), a pro-apoptotic mitochondrial protein of the Bcl-2 family, which is known to be involved in the intrinsic mitochondrial apoptotic pathway . Although the expressional changes in microarray experiments were lower than for the two other validated genes (mean FC of 9.0 and 10.2, respectively), comparable results could be obtained by qRT-PCR for PPP1R15A (3-fold to 18 fold increase) and PMAIP1 (2 fold to 10 fold increase) (figure 5c+d). TRAF6 (TNF receptor associated factor 6), a member of the E3 ubiquitin ligase family which can specificly target different proteins for Lys63 ubiquitination [35, 36] was significantly down-regulated both in microarray analysis (mean FC of 4.8) and in qRT-PCR (5 to 9 fold decrease) (figure 5e). To further confirm altered gene expression of potential TRD target genes, western blots were performed when appropriate antibodies were available. As indicated in figure 6a, all cell lines displayed a pronounced upregulation of ATF3 protein, which was mirrored by the gene expression data as previously mentioned. PPP1R15A was also expressional regulated on protein level (figure 6b). However, in HT1080 cells we could not observe a convincing increase in PPP1R15A protein (figure 6b), which was also less pronounced on mRNA level (figure 5c).
Taurolidine (TRD) represents an interesting anti-neoplastic agent with a potential perspective in oncologic pharmacotherapy. Although it has already been intravenously applied to patients with advanced gastric cancer  and glioblastoma during pilot studies with promising results [19, 18], the mechanisms of action remains to be elucidated in detail [4, 5]. Microarray technology offers a powerful tool for investigating the cellular responses to anti-neoplastic substances like TRD because it monitors thousands of genes simultaneously. Although two recent publications from our group have already focussed on microarray derived transcriptional profiling of TRD treatment in two single cell lines [16, 14], the current project provides the first microarray analysis which was performed simultaneously in five different and heterogeneous cell lines of four different malignancies to identify potential common TRD target genes.
Microarray experiments were focussed on a TRD concentration of 250 μM. We have recently shown in a study analysing TRD effects by FACS analysis in five cell lines identical with the current study (HT29, Chang Liver, HT1080, AsPC-1 and BxPC-3) that TRD inhibits cell viability over a broad range of concentrations (100 μM, 250 μM and 1000 μM) . In the latter study, TRD 250 μM was the lowest concentration that significantly inhibited cell viability in all cell lines . This fact is now supported by the current study, since TRD 250 μM significantly induced apoptotic cell death in our experiments - as measured by FACS analysis (figure 2+3). Furthermore, TRD 250 μM was the lowest concentration that significantly inhibited proliferation in all cell lines as measured by BrdU assay (figure 1). The effectiveness of TRD 250 μM as powerful anti-neoplastic concentration in cell culture experiments is consistent with several studies that evaluated the anti-proliferative and cell death inducing activity of TRD over a broad range of concentrations (50 μM - 2000 μM) in different cell lines using proliferation assays [37, 38, 11, 39, 10, 16] or quantitative FACS analysis [16, 14, 13, 11]. As indicated in figure 1 and 2, our results showed quite distinct and cell line specific differences in susceptibility towards TRD induced cell death. One can only speculate about the differences in the anti-neoplastic activity of TRD. However, none of the cell lines tested in our study was resistent to TRD induced cell death.
Using cRNA Microarray assay combined with Ingenuity pathways analysis we identified several TRD regulated genes that might represent potential target genes for TRD in malignant cells. Among the 21 genes with ≥ 4 fold expressional change (table 1), we focussed in our validation experiments on 5 genes, that represent - in our eyes - genes of special interest: the transcription factors (TFs) EGR1 and ATF3, the ER-stress response gene PPP1R15A, the pro-apoptotic mitochondrial protein PMAIP1 and the E3 ubiquitin ligase TRAF6. EGR1 (early growth response 1) representing the gene exhibiting the strongest induction upon TRD treatment was confirmed by qRT-PCR. EGR1 belongs to the EGR zinc-finger family of TFs  which can provide anti-neoplastic activity in numerous tumour cell lines by inducing tumour suppressors like PTEN, TGFβ1, p53, fibronectin or ATF3 [40, 30, 41]. Furthermore, several anti-neoplastic or chemopreventive agents have recently been shown to target EGR1 expression and EGR1 driven pathways in colorectal cancer [42–44]. Furthermore, our microarray analyses revealed a strong up-regulation of the pro-apoptotic transcription factor ATF3 (table 1), whereas other EGR1 target genes with tumour suppressor function like PTEN, TGFβ1, p53 and fibronectin did not show any significant regulation (additional file 2, table S2). The common induction of of EGR1 and ATF3 expression in all lines tested might indicate that activation of the EGR1/ATF3 signalling axis is an important mode of regulation how TRD mediates its growth suppressive action and/or apoptosis induction. The importance of EGR1/ATF3 axis in death signalling is supported by recent studies that described an EGR1 driven activation of the pro-apoptotic transcription factor ATF3 in colorectal cancer cells following treatment with anti-neoplastic substances [44, 45]. ATF3 itself belongs to the CREB family of transcription factors and is regarded as a stress-inducible gene due to its activation by variety of stimuli e.g. radiation , chemotherapeutic agents [46, 47] oxidative stress [48, 49] or ER stress . Although the physiologic functions of ATF3 are not well understood, there is some evidence that ATF3 functions as a pro-apoptotic gene in cancer cells. ATF3 has been shown to enhance the apoptotic effect of curcumin in squamous cell carcinoma cells  and to mediate nitric oxide-induced apoptosis in pancreatic beta cells . Furthermore, ATF3 has been shown to promote UV induced apoptosis and cell cycle arrest in fibroblasts by binding the promoter and repressing the transcription of cyclin D1 -a key regulator of the G1-S checkpoint . We can only speculate about the TRD driven stimulation of ATF3, but oxidative stress provides a possible pathway. Previous studies have presented first evidence for involvement of TRD mediated reactive oxygen species (ROS) production in malignant cells [7, 25, 26] and cell death induced by TRD has been shown to be reversible by application of radical scavengers like N-Acetylcysteine [7, 54, 25, 10, 26]. Another possible cell death associated pathway resulting in activation of ATF3 by TRD is the ER stress response pathway [55, 56], since TRD has been reported to inhibit protein translation in tumour cells at an early stage . Together with the above mentioned production of ROS, this might result in disturbance of ER homeostasis leading to ER stress response. Indicative of an involvement of the ER stress response pathway, we identified two genes that were up-regulated by TRD and known to be activated during ER stress, e.g. PPP1R15A (syn. GADD34) and the above mentioned ATF3. PPP1R15A represents an interesting potential TRD target gene, since PPP1R15A confers both apoptotic signalling [32, 33] and inhibition of the translational silencing during ER stress via binding the catalytic subunit of protein phosphatase 1 leading to dephosphorylation of eIF2α [57, 58]. Another strongly TRD regulated gene is PMAIP1 (syn. NOXA), a pro-apoptotic member of the BH3-only protein family, which is known to be involved in the intrinsic mitochondrial apoptotic pathway . PMAIP1 exerts its pro-apoptotic function mainly by neutralizing anti-apoptotic members of the Bcl2 family through direct interaction, which has been proven for Bcl-2 related protein A1 (BCL2A1) and myeloid cell leukemia sequence 1 (MCL1) [59, 60]. Interestingly, besides Bcl-2 like protein 11 (BCL2L11), which was 3-fold down-regulated, none of the remaining Bcl-2 family members (including MCL1 and BCL2A1) revealed a significant expressional change after TRD treatment in our experiments (additional file 2, table S2). Among the genes that were significantly down-regulated following TRD treatment, TRAF6 was considered as an especially important candidate gene. TRAF6 represents a member of the TRAF family of proteins that is involved in both the TNF receptor and the interleukin-1 receptor (IL-1R)/Toll-like receptor (TLR) superfamily signal transduction and is referred to innate and adaptive immunity, bone metabolism and developmental pathways [61, 62]. Recently it has been shown that TRAF6 possesses an ubiquitin ligase activity domain (E3 ligase) that confers a site specific (Lys63) ubiquitination of different target proteins [35, 36]. Interestingly, the protein kinase AKT - activated and transmitting survival signals in many tumours - has also been reported to be an ubiquitination target of TRAF6. Ubiquitination of AKT by TRAF6 is essential for its membrane recruitment and phosphorylation following growth factor activation . Furthermore, inhibition of AKT phosphorylation has been described following TRD application in mesothelioma cells . Although we did not find significant changes in the expression of AKT (additional file 2, table S2) and did not analyse AKT phosphorylation, the reduced expression of TRAF6 in our experiments provides a potential hint to the inhibition of the oncogenic AKT pathway by TRD. However, the involvement of the above mentioned pathways in TRD mediated cell death has to be addressed in further studies.
This study provides the first identification of potential common target genes modulated by TRD treatment of cell lines originating from four different cancer types. Among the TRD regulated genes, we identified pro-apoptotic transcription factors, cell cycle regulators and proteins involved in ubiquitination, endoplasmic reticulum response and mitochondrial apoptotic pathways. Our results indicate that TRD is involved in different signal transduction pathways leading to multifaceted cell death mechanisms. Further studies are necessary to address the potential target genes in functional assays.
Activating transcription factor 3
Cyclin dependent kinase inhibitor 2B
E2F transcription factor 3
Early growth response 1
v-fos FBJ murine osteosarcoma viral oncogene homolog
FBJ murine osteosarcoma viral oncogene homolog B
Growth arrest and DNA-damage inducible 34
Growth arrest and DNA-damage inducible, alpha
Growth arrest and DNA-damage inducible, beta
Hairy and enhancer of split 1
Nuclear receptor subfamily 4, group A, member 2
Phorbol-12-myristate-13-acetate-induced protein 1 (Syn. NOXA)
Protein phosphatase 1, regulatory (inhibitory) subunit 15A (Syn. GADD34)
Snail homolog 1
TNF receptor-associated factor 6
Ubiquitin-conjugating enzyme E2B.
Baker DM, Jones JA, Nguyen-Van-Tam JS, Lloyd JH, Morris DL, Bourke JB, Steele RJ, Hardcastle JD: Taurolidine peritoneal lavage as prophylaxis against infection after elective colorectal surgery. Br J Surg. 1994, 81: 1054-1056. 10.1002/bjs.1800810743.
Simon A, Ammann RA, Wiszniewsky G, Bode U, Fleischhack G, Besuden MM: Taurolidine-citrate lock solution (TauroLock) significantly reduces CVAD-associated grampositive infections in pediatric cancer patients. BMC Infect Dis. 2008, 8: 102-10.1186/1471-2334-8-102.
Koldehoff M, Zakrzewski JL: Taurolidine is effective in the treatment of central venous catheter-related bloodstream infections in cancer patients. Int J Antimicrob Agents. 2004, 24: 491-495. 10.1016/j.ijantimicag.2004.06.006.
Jacobi CA, Menenakos C, Braumann C: Taurolidine--a new drug with anti-tumor and anti-angiogenic effects. Anticancer Drugs. 2005, 16: 917-921. 10.1097/01.cad.0000176502.40810.b0.
Neary PM, Hallihan P, Wang JH, Pfirrmann RW, Bouchier-Hayes DJ, Redmond HP: The Evolving Role of Taurolidine in Cancer Therapy. Ann Surg Oncol. 2010, 17 (4): 1135-43. 10.1245/s10434-009-0867-9.
Stendel R, Stoltenburg-Didinger G, Al Keikh CL, Wattrodt M, Brock M: The effect of taurolidine on brain tumor cells. Anticancer Res. 2002, 22: 809-814.
Stendel R, Biefer HR, Dekany GM, Kubota H, Munz C, Wang S, Mohler H, Yonekawa Y, Frei K: The antibacterial substance taurolidine exhibits anti-neoplastic action based on a mixed type of programmed cell death. Autophagy. 2009, 5: 194-210. 10.4161/auto.5.2.7404.
Braumann C, Jacobi CA, Rogalla S, Menenakos C, Fuehrer K, Trefzer U, Hofmann M: The tumor suppressive reagent taurolidine inhibits growth of malignant melanoma--a mouse model. J Surg Res. 2007, 143: 372-378. 10.1016/j.jss.2007.01.041.
Sun BS, Wang JH, Liu LL, Gong SL, Redmond HP: Taurolidine induces apoptosis of murine melanoma cells in vitro and in vivo by modulation of the Bcl-2 family proteins. J Surg Oncol. 2007, 96: 241-248. 10.1002/jso.20827.
Opitz I, Sigrist B, Hillinger S, Lardinois D, Stahel R, Weder W, Hopkins-Donaldson S: Taurolidine and povidone-iodine induce different types of cell death in malignant pleural mesothelioma. Lung Cancer. 2007, 56: 327-336. 10.1016/j.lungcan.2007.01.024.
Aceto N, Bertino P, Barbone D, Tassi G, Manzo L, Porta C, Mutti L, Gaudino G: Taurolidine and oxidative stress: a rationale for local treatment of mesothelioma. Eur Respir J. 2009, 34: 1399-1407. 10.1183/09031936.00102308.
Braumann C, Schoenbeck M, Menenakos C, Kilian M, Jacobi CA: Effects of increasing doses of a bolus injection and an intravenous long-term therapy of taurolidine on subcutaneous (metastatic) tumor growth in rats. Clin Exp Metastasis. 2005, 22: 77-83. 10.1007/s10585-005-4695-4.
Chromik AM, Daigeler A, Hilgert C, Bulut D, Geisler A, Liu V, Otte JM, Uhl W, Mittelkotter U: Synergistic effects in apoptosis induction by taurolidine and TRAIL in HCT-15 colon carcinoma cells. J Invest Surg. 2007, 20: 339-348. 10.1080/08941930701772157.
Daigeler A, Chromik AM, Geisler A, Bulut D, Hilgert C, Krieg A, Klein-Hitpass L, Lehnhardt M, Uhl W, Mittelkötter U: Synergistic apoptotic effects of taurolidine and TRAIL on squamous carcinoma cells of the esophagus. Int J Oncol. 2008, 32: 1205-1220. 10.3892/ijo_32_6_1205.
Walters DK, Muff R, Langsam B, Gruber P, Born W, Fuchs B: Taurolidine: a novel anti-neoplastic agent induces apoptosis of osteosarcoma cell lines. Invest New Drugs. 2007, 25: 305-312. 10.1007/s10637-007-9052-9.
Daigeler A, Brenzel C, Bulut D, Geisler A, Hilgert C, Lehnhardt M, Steinau HU, Flier A, Steinstraesser L, Klein-Hitpass L, et al: TRAIL and Taurolidine induce apoptosis and decrease proliferation in human fibrosarcoma. J Exp Clin Cancer Res. 2008, 27: 82-10.1186/1756-9966-27-82.
Gong L, Greenberg HE, Perhach JL, Waldman SA, Kraft WK: The pharmacokinetics of taurolidine metabolites in healthy volunteers. J Clin Pharmacol. 2007, 47: 697-703. 10.1177/0091270007299929.
Braumann C, Winkler G, Rogalla P, Menenakos C, Jacobi CA: Prevention of disease progression in a patient with a gastric cancer-re-recurrence. Outcome after intravenous treatment with the novel antineoplastic agent taurolidine. Report of a case. World J Surg Oncol. 2006, 4: 34-10.1186/1477-7819-4-34.
Stendel R, Picht T, Schilling A, Heidenreich J, Loddenkemper C, Janisch W, Brock M: Treatment of glioblastoma with intravenous taurolidine. First clinical experience. Anticancer Res. 2004, 24: 1143-1147.
Stendel R, Scheurer L, Schlatterer K, Stalder U, Pfirrmann RW, Fiss I, Mohler H, Bigler L: Pharmacokinetics of taurolidine following repeated intravenous infusions measured by HPLC-ESI-MS/MS of the derivatives taurultame and taurinamide in glioblastoma patients. Clin Pharmacokinet. 2007, 46: 513-524. 10.2165/00003088-200746060-00005.
Darnowski JW, Goulette FA, Cousens LP, Chatterjee D, Calabresi P: Mechanistic and antineoplastic evaluation of taurolidine in the DU145 model of human prostate cancer. Cancer Chemother Pharmacol. 2004, 54: 249-258. 10.1007/s00280-004-0806-1.
Han Z, Ribbizi I, Pantazis P, Wyche J, Darnowski J, Calabresi P: The antibacterial drug taurolidine induces apoptosis by a mitochondrial cytochrome c-dependent mechanism. Anticancer Res. 2002, 22: 1959-1964.
Rodak R, Kubota H, Ishihara H, Eugster HP, Konu D, Mohler H, Yonekawa Y, Frei K: Induction of reactive oxygen intermediates-dependent programmed cell death in human malignant ex vivo glioma cells and inhibition of the vascular endothelial growth factor production by taurolidine. J Neurosurg. 2005, 102: 1055-1068. 10.3171/jns.2005.102.6.1055.
Stendel R, Scheurer L, Stoltenburg-Didinger G, Brock M, Mohler H: Enhancement of Fas-ligand-mediated programmed cell death by taurolidine. Anticancer Res. 2003, 23: 2309-2314.
Rodak R, Kubota H, Ishihara H, Eugster H-P, Könü D, Möhler H, Yonekawa Y, Frei K: Induction of reactive oxygen intermediates-dependent programmed cell death in human malignant ex vivo glioma cells and inhibition of the vascular endothelial growth factor production by taurolidine. J Neurosurg. 2005, 102: 1055-1068. 10.3171/jns.2005.102.6.1055.
Chromik AM, Daigeler A, Bulut D, Flier A, May C, Harati K, Roschinsky J, Sulberg D, Ritter PR, Mittelkotter U, et al: Comparative analysis of cell death induction by Taurolidine in different malignant human cancer cell lines. J Exp Clin Cancer Res. 2010, 29: 21-10.1186/1756-9966-29-21.
t Hoen PA, Turk R, Boer JM, Sterrenburg E, de Menezes RX, van Ommen GJ, den Dunnen JT: Intensity-based analysis of two-colour microarrays enables efficient and flexible hybridization designs. Nucleic Acids Res. 2004, 32: e41-10.1093/nar/gnh038.
Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.
Gashler A, Sukhatme VP: Early growth response protein 1 (Egr-1): prototype of a zinc-finger family of transcription factors. Prog Nucleic Acid Res Mol Biol. 1995, 50: 191-224. full_text.
Baron V, Adamson ED, Calogero A, Ragona G, Mercola D: The transcription factor Egr1 is a direct regulator of multiple tumor suppressors including TGFbeta1, PTEN, p53, and fibronectin. Cancer Gene Ther. 2006, 13: 115-124. 10.1038/sj.cgt.7700896.
Fan F, Jin S, Amundson SA, Tong T, Fan W, Zhao H, Zhu X, Mazzacurati L, Li X, Petrik KL, et al: ATF3 induction following DNA damage is regulated by distinct signaling pathways and over-expression of ATF3 protein suppresses cells growth. Oncogene. 2002, 21: 7488-7496. 10.1038/sj.onc.1205896.
Hollander MC, Sheikh MS, Yu K, Zhan Q, Iglesias M, Woodworth C, Fornace AJ: Activation of Gadd34 by diverse apoptotic signals and suppression of its growth inhibitory effects by apoptotic inhibitors. Int J Cancer. 2001, 96: 22-31. 10.1002/1097-0215(20010220)96:1<22::AID-IJC3>3.0.CO;2-K.
Adler HT, Chinery R, Wu DY, Kussick SJ, Payne JM, Fornace AJ, Tkachuk DC: Leukemic HRX fusion proteins inhibit GADD34-induced apoptosis and associate with the GADD34 and hSNF5/INI1 proteins. Mol Cell Biol. 1999, 19: 7050-7060.
Ploner C, Kofler R, Villunger A: Noxa: at the tip of the balance between life and death. Oncogene. 2008, 27 (Suppl 1): S84-92. 10.1038/onc.2009.46.
Lamothe B, Besse A, Campos AD, Webster WK, Wu H, Darnay BG: Site-specific Lys-63-linked tumor necrosis factor receptor-associated factor 6 auto-ubiquitination is a critical determinant of I kappa B kinase activation. J Biol Chem. 2007, 282: 4102-4112. 10.1074/jbc.M609503200.
Sorrentino A, Thakur N, Grimsby S, Marcusson A, von Bulow V, Schuster N, Zhang S, Heldin CH, Landstrom M: The type I TGF-beta receptor engages TRAF6 to activate TAK1 in a receptor kinase-independent manner. Nat Cell Biol. 2008, 10: 1199-1207. 10.1038/ncb1780.
Braumann C, Henke W, Jacobi CA, Dubiel W: The tumor-suppressive reagent taurolidine is an inhibitor of protein biosynthesis. Int J Cancer. 2004, 112: 225-230. 10.1002/ijc.20393.
McCourt M, Wang JH, Sookhai S, Redmond HP: Taurolidine inhibits tumor cell growth in vitro and in vivo. Ann Surg Oncol. 2000, 7: 685-691. 10.1007/s10434-000-0685-6.
Nici L, Monfils B, Calabresi P: The effects of taurolidine, a novel antineoplastic agent, on human malignant mesothelioma. Clin Cancer Res. 2004, 10: 7655-7661. 10.1158/1078-0432.CCR-0196-03.
Nair P, Muthukkumar S, Sells SF, Han SS, Sukhatme VP, Rangnekar VM: Early growth response-1-dependent apoptosis is mediated by p53. J Biol Chem. 1997, 272: 20131-20138. 10.1074/jbc.272.32.20131.
Thiel G, Cibelli G: Regulation of life and death by the zinc finger transcription factor Egr-1. J Cell Physiol. 2002, 193: 287-292. 10.1002/jcp.10178.
Baek SJ, Kim JS, Moore SM, Lee SH, Martinez J, Eling TE: Cyclooxygenase inhibitors induce the expression of the tumor suppressor gene EGR-1, which results in the up-regulation of NAG-1, an antitumorigenic protein. Mol Pharmacol. 2005, 67: 356-364. 10.1124/mol.104.005108.
Baek SJ, Wilson LC, Hsi LC, Eling TE: Troglitazone, a peroxisome proliferator-activated receptor gamma (PPAR gamma) ligand, selectively induces the early growth response-1 gene independently of PPAR gamma. A novel mechanism for its anti-tumorigenic activity. J Biol Chem. 2003, 278: 5845-5853. 10.1074/jbc.M208394200.
Cho KN, Sukhthankar M, Lee SH, Yoon JH, Baek SJ: Green tea catechin (-)-epicatechin gallate induces tumour suppressor protein ATF3 via EGR-1 activation. Eur J Cancer. 2007, 43: 2404-2412. 10.1016/j.ejca.2007.07.020.
Yamaguchi K, Lee SH, Kim JS, Wimalasena J, Kitajima S, Baek SJ: Activating transcription factor 3 and early growth response 1 are the novel targets of LY294002 in a phosphatidylinositol 3-kinase-independent pathway. Cancer Res. 2006, 66: 2376-2384. 10.1158/0008-5472.CAN-05-1987.
Mallory JC, Crudden G, Oliva A, Saunders C, Stromberg A, Craven RJ: A novel group of genes regulates susceptibility to antineoplastic drugs in highly tumorigenic breast cancer cells. Mol Pharmacol. 2005, 68: 1747-1756.
De Angelis PM, Svendsrud DH, Kravik KL, Stokke T: Cellular response to 5-fluorouracil (5-FU) in 5-FU-resistant colon cancer cell lines during treatment and recovery. Mol Cancer. 2006, 5: 20-10.1186/1476-4598-5-20.
Allen-Jennings AE, Hartman MG, Kociba GJ, Hai T: The roles of ATF3 in glucose homeostasis. A transgenic mouse model with liver dysfunction and defects in endocrine pancreas. J Biol Chem. 2001, 276: 29507-29514. 10.1074/jbc.M100986200.
Hai T, Hartman MG: The molecular biology and nomenclature of the activating transcription factor/cAMP responsive element binding family of transcription factors: activating transcription factor proteins and homeostasis. Gene. 2001, 273: 1-11. 10.1016/S0378-1119(01)00551-0.
Jiang HY, Wek SA, McGrath BC, Lu D, Hai T, Harding HP, Wang X, Ron D, Cavener DR, Wek RC: Activating transcription factor 3 is integral to the eukaryotic initiation factor 2 kinase stress response. Mol Cell Biol. 2004, 24: 1365-1377. 10.1128/MCB.24.3.1365-1377.2004.
Yan C, Jamaluddin MS, Aggarwal B, Myers J, Boyd DD: Gene expression profiling identifies activating transcription factor 3 as a novel contributor to the proapoptotic effect of curcumin. Mol Cancer Ther. 2005, 4: 233-241.
Hartman MG, Lu D, Kim M-L, Kociba GJ, Shukri T, Buteau J, Wang X, Frankel WL, Guttridge D, Prentki M, et al: Role for activating transcription factor 3 in stress-induced beta-cell apoptosis. Mol Cell Biol. 2004, 24: 5721-5732. 10.1128/MCB.24.13.5721-5732.2004.
Lu D, Wolfgang CD, Hai T: Activating transcription factor 3, a stress-inducible gene, suppresses Ras-stimulated tumorigenesis. J Biol Chem. 2006, 281: 10473-10481. 10.1074/jbc.M509278200.
Aceto N, Bertino P, Barbone D, Tassi G, Manzo L, Porta C, Mutti L, Gaudino G: Taurolidine and oxidative stress: a rationale for local treatment of mesothelioma. Eur Respir J. 2009, 34: 1399-1407. 10.1183/09031936.00102308.
Healy SJ, Gorman AM, Mousavi-Shafaei P, Gupta S, Samali A: Targeting the endoplasmic reticulum-stress response as an anticancer strategy. Eur J Pharmacol. 2009, 625: 234-246. 10.1016/j.ejphar.2009.06.064.
Kim I, Xu W, Reed JC: Cell death and endoplasmic reticulum stress: disease relevance and therapeutic opportunities. Nat Rev Drug Discov. 2008, 7: 1013-1030. 10.1038/nrd2755.
Novoa I, Zeng H, Harding HP, Ron D: Feedback inhibition of the unfolded protein response by GADD34-mediated dephosphorylation of eIF2alpha. J Cell Biol. 2001, 153: 1011-1022. 10.1083/jcb.153.5.1011.
Novoa I, Zhang Y, Zeng H, Jungreis R, Harding HP, Ron D: Stress-induced gene expression requires programmed recovery from translational repression. EMBO J. 2003, 22: 1180-1187. 10.1093/emboj/cdg112.
Oda E, Ohki R, Murasawa H, Nemoto J, Shibue T, Yamashita T, Tokino T, Taniguchi T, Tanaka N: Noxa, a BH3-only member of the Bcl-2 family and candidate mediator of p53-induced apoptosis. Science. 2000, 288: 1053-1058. 10.1126/science.288.5468.1053.
Shibue T, Taniguchi T: BH3-only proteins: integrated control point of apoptosis. Int J Cancer. 2006, 119: 2036-2043. 10.1002/ijc.21751.
Inoue J, Gohda J, Akiyama T: Characteristics and biological functions of TRAF6. Adv Exp Med Biol. 2007, 597: 72-79. full_text.
Chung JY, Lu M, Yin Q, Lin SC, Wu H: Molecular basis for the unique specificity of TRAF6. Adv Exp Med Biol. 2007, 597: 122-130. full_text.
Yang W-L, Wang J, Chan C-H, Lee S-W, Campos AD, Lamothe B, Hur L, Grabiner BC, Lin X, Darnay BG, Lin H-K: The E3 ligase TRAF6 regulates Akt ubiquitination and activation. Science. 2009, 325: 1134-1138. 10.1126/science.1175065.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/10/595/prepub
The authors thank Prof. Dr. W. E. Schmidt (Department of Internal Medicine I, St. Josef Hospital, Ruhr-University of Bochum) as well as Prof. Dr. A. Mügge (Department of Internal Medicine II, St. Josef Hospital, Ruhr-University of Bochum) for generously supporting cell culture experiments and FACS analysis. Furthermore, they thank Kirsten Mros and Rainer Lebert (Gastrointestinal Research Laboratory, St. Josef Hospital, Ruhr-University of Bochum) for technical assistance. This study was supported by FoRUM AZ F472-2005 and FoRUM AZ F544-2006 from the Ruhr-University Bochum, Germany.
AMC received financial support by Geistlich Pharma (Suisse) for laboratory experiments. All other authors declare that they have no competing interests.
AMC and SAH conceived of the study and its design, coordinated the experiments, carried out the statistical analysis and drafted the manuscript. AD and AF supervised the cell culture experiments and microarray experiments. DB was responsible for adjusting the FACS analysis and helped to draft the manuscript. CM, KH and JR carried out the cell culture experiments. DS helped with the statistical analysis and revised manuscript. DW, UM and WU participated in the design and coordination of the study and revised the manuscript. All authors have read and approved the final manuscript.
Ansgar M Chromik, Stephan A Hahn contributed equally to this work.
Electronic supplementary material
Authors’ original submitted files for images
About this article
Cite this article
Chromik, A.M., Hahn, S.A., Daigeler, A. et al. Gene expression analysis of cell death induction by Taurolidine in different malignant cell lines. BMC Cancer 10, 595 (2010). https://doi.org/10.1186/1471-2407-10-595
- Ingenuity Pathway Analysis
- Chang Liver
- Malignant Cell Line
- Chang Liver Cell