| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Regular Articles |
Department of Genetics and Genomics, Boston University School of Medicine, Boston, Massachusetts
| ABSTRACT |
|---|
|
|
|---|
B-regulated antiapoptotic factors. Experiments using HeLa cells in our lab and prior animal model studies performed elsewhere confirm that inhibitors of these respective pathways super-additively enhance CPTs cytotoxicity, suggesting their potential as targets for adjuvant therapy with CPT. | INTRODUCTION |
|---|
|
|
|---|
Camptothecin (CPT) is a plant alkaloid derived from the tree, Camptotheca acuminata. Its derivatives have shown promise in clinical trials for use against small cell lung carcinoma (3) and other malignancies, including cervical cancer (4) . The compound exhibits high specificity for topoisomerase I (topo I), which it inhibits by preventing the religation step of this enzymes DNA relaxing activity. In the presence of CPT, topo I forms a covalent intermediate with the 3'-phosphate end of the broken DNA strand (5) . This single-strand nick lesion is believed to convert into an irreversible double-strand break after it encounters the replication fork machinery in S phase. This unusual form of DNA damage induces a signaling cascade that is propagated by ATM-family kinases in human cells (6) , which phosphorylate numerous targets involved in transcription, DNA repair, cell cycle arrest, and apoptosis.
A comprehensive survey of the entire transcriptome after CPT treatment via the Affymetrix U133A Chip has not been reported, although one study examined a small number of genes after CPT treatment (7)
. Moreover, it remains unclear whether such information could be of therapeutic value. In the following study, we present data generated using the Affymetrix U133A chip (22,283 probe sets) from total RNA collected from HeLa cells treated with CPT at a dose (10 µM) and time (8 h) sufficient to induce rapid apoptosis. We observed significant 3-fold changes in expression of 683 genes, including novel transcripts and several putative p53-regulated genes. Consistent with the induction of certain p53-responsive genes, we observed increased protein levels of p53, phosphorylation on serine 15, and S-phase arrest after CPT treatment. Moreover, we detected striking evidence of nuclear factor (NF)-
B transactivation and epidermal growth factor receptor (EGFR) survival signaling (via HB-EGF). We reasoned that changes in the gene expression pattern might represent regulation of pathways important for CPT survival. The data were used to generate the hypothesis that inhibition of the NF-
B and EGFR pathways would increase the sensitivity of HeLa cells to CPT. Administration of an inhibitor of this pathway, AG1478, yielded a superadditive increase in apoptosis when administered in combinations with CPT. Inhibition of NF-
B signaling has been shown to dramatically increase tumor cell sensitivity to CPT-11 (8
, 9)
. These results suggest both that CPT may be effective in mobilizing a p53 response in cervical carcinomas and that data from microarray experiments can form the basis for rational hypotheses in the treatment of cancer.
| MATERIALS AND METHODS |
|---|
|
|
|---|
40,000 cells/cm2 (
75% confluence) using an RNAeasy mini-kit (Qiagen, Valencia, CA).
Chemicals.
All chemicals and reagents, unless otherwise specified, were purchased from Sigma-Aldrich (St. Louis, MO).
Cell Culture.
HeLa cells were purchased from the American Type Culture Collection (Manassas, VA) and maintained at 37°C in 5% CO2 in DMEM (Mediatech, Herndon, VA) containing 10% fetal bovine serum supplemented with penicillin/streptomycin (Hyclone, Logan, UT). Experiments listed here were conducted on cells that had been passaged three to seven times following arrival from American Type Culture Collection.
Apoptosis Response and p53 Activation.
HeLa cells were plated at a density of 106 cells/6 cm-diameter tissue culture dish and incubated 16 h. Culture medium was then removed and replaced with media containing varying doses of CPT and/or inhibitors. At times of data collection, both floating (apoptotic and/or necrotic cells) and adherent cell material from each dish were collected on ice using cell lifters and centrifugation (5 min x 200 g). Next, cell pellets were resuspended in ice-cold PBS, centrifuged again (5 min x 200 g), and lysed in 200 µl of lysis buffer [1% Igepal, 0.5% deoxycholate, 2 mM DTT, 150 mM NaCl, and 20 mM HEPES (pH 7.4); supplemented with phosphatase and protease inhibitors]. Protein concentrations in lysates were determined by the Bio-Rad assay, and 100 µg/lane samples were immediately prepared for loading onto 12% acrylamide SDS-PAGE by boiling (5 min) in 1% SDS sample buffer. Completed gels were then electrotransferred to nitrocellulose membranes, which were stained with Ponceau S and cut in half across protein bands corresponding to
Mr 35,000. Membrane halves corresponding to higher molecular weight proteins were Western blotted with primary antibodies directed against whole p53 protein (Ab-12; Oncogene Research Products, San Diego, CA) and serine-15-phosphorylated p53 (Cell Signaling Technology, Beverly, MA); membranes corresponding to the lower molecular weight proteins were Western blotted with an antibody directed against the Mr 23,000 caspase-specific cleavage product of poly(ADP-ribose) polymerase (R&D Systems, Minneapolis, MN). These Western blot procedures were conducted according to the directions of the manufacturer, and immunoreactive bands were visualized by chemiluminescent exposure to instant film such that subsaturation exposures (
1 min) were obtained. These film images were then digitally scanned and converted to .tiff images by Adobe PhotoShop with the "Autolevels" (Shft+Ctrl+L) image function selected.
Scanned Western blot images were then analyzed with the Scion Image software package freely available from ScionCorp.1 The immunoreactive bands were highlighted, and pixel counts of each band were obtained via the Integrated Density analysis option. Integrated densities were linearly normalized such that the gel band corresponding to the positive control in the experiment (10 µM CPT at 8 h.) was assigned an immunoreactivity value of 100.
WST-1 Viability Assay.
HeLa cells were plated at a density of 10,000 cells/well on a 96-well tissue culture plate. Twenty-four h later, media (DMEM +10% FCS) containing varying concentrations of CPT +/- AG1478 (500 nM) were added. Twelve h later, 10 µl/well of the viability stain, WST-1 (Roche Biochemicals, Indianapolis, IN), were added for 30 min. Colorimetric readings were obtained by plate reader at 450 nm. The average of data points from wells containing WST-1-stained media, alone (with no cells), was subtracted from the raw data from wells containing cells. Percent decline in viability/well was calculated as
![]() |
RNA Labeling and Hybridization.
Using a poly-dT primer incorporating a T7 promoter, double-stranded cDNA was synthesized from 10 µg of total RNA using a Superscript cDNA Synthesis kit (Invitrogen, Carlsbad, CA). Double-stranded cDNA was purified by phenol/chloroform extraction. The aqueous phase was isolated using Phase-Lock Gel Heavy (Brinkmann Instruments, Westbury, NY) and the cDNA ethanol precipitated. Subsequently, biotin-labeled cRNA was generated from the double-stranded cDNA template though in vitro transcription with T7 polymerase using a BioArray High Yield RNA Transcript Labeling kit (Enzo Diagnostics, Farmingdale, NY). The biotinylated cRNA was purified using RNeasy affinity columns (Qiagen). Biotinylated cRNA (20 µg) was fragmented in 40 mM Tris-acetate (pH 8.1), 100 mM KOAc, and 30 mM MgOAc for 35 min at 94°C to an average size of 35200 bases. Fragmented, biotinylated cRNA (10 µg), along with hybridization controls (Affymetrix, Santa Clara, CA), were hybridized to Affymetrix Human Genome U133A GeneChip arrays containing probes for 22,283 transcripts. The arrays were hybridized for 16 h at 45°C and 60 rpm. After hybridization, arrays were washed and stained according to the standard Antibody Amplification for Eukaryotic Targets protocol (Affymetrix). The stained GeneChip arrays were scanned at 488 nm using a G2500AGeneArray Scanner (Agilent, Palo Alto, CA) and Microarray Suite 5.0 software (Affymetrix).
Data Quantification and Normalization.
After data acquisition, the scanned images were quantified using Microarray Suite 5.0 software (Affymetrix) yielding a signal intensity for each probe on the GeneChip. The signal intensities from the 22 probes for each gene were then used to determine an overall expression level according to algorithms implemented in Microarray Suite 5.0. The arrays were then linearly scaled to an average expression level of 500 units on each chip. Scaling factors varied from 2.29 to 4.30. Upon publication, complete microarray data for each sample will be made publicly available through National Center for Biotechnology Informations Gene Expression Omnibus.2
| RESULTS |
|---|
|
|
|---|
Data Analysis.
For each gene, fold-change (FC) and statistical significance of differential expression (p) were calculated. FC was calculated using the average signal from the two groups. Statistical significance was calculated using a t-Test implemented in Cyber-T.3
Cyber-T uses a Bayesian estimate of the variance within treatment groups based on a prior distribution obtained from an estimate of the variance from genes at a similar expression level. This provides increased power to identify differentially expressed genes because the variability within groups is estimated based on a larger number of data points.
To correct significance estimates for multiple hypothesis testing we used a modified version the false discovery rate (FDR) method of Benjamini and Hochberg (10)
, which calculates a new test statistic that estimates the portion of type I errors within a group of genes meeting a significance cutoff. This modified FDR is the quotient of the number of unchanged genes expected at a given significance cutoff over the number of genes detected at that significance cutoff.
![]() |
Two assumptionsthat changed genes would not have p > 0.90 and that unchanged genes would have an even distribution of Ps between 0 and 1were used to estimate the number of unchanged genes. A total of 1,412 genes had p > 0.90; thus, it was estimated that 14,120 genes were not differentially expressed between the treatment groups. Because we have assumed that these unchanged genes will have an even distribution of p values, p * 14,120 genes that were not truly differentially expressed can be expected to meet a significance cutoff p by chance alone. Using this model of type I errors, the set of FDRs associated with each gene was calculated as:
![]() |
The FDR of a gene was then defined as the minimal FDR at which time the gene appears significant.
To define a set of differentially expressed genes, we imposed cutoffs based on statistical significance and FC. Probe sets with FDR < 0.05 were sufficiently reproducible to indicate differential expression (of any magnitude). These probe sets can be expected to have a type I (false positive) error rate of 5%. A total of 2991 probe sets were identified at FDR < 0.05, indicating that we expect
2840 of these genes to in fact be differentially expressed. Genes with FC > 3 were considered sufficiently changed to indicate biological meaning. A total of 1009 genes, the expression of which was 3-fold changed as a result of CPT treatment, was detected. The criterion for biological significance, a change of
3-fold, is clearly somewhat arbitrary and results in the exclusion of many genes meeting the significance threshold from additional analysis. Combining the FC a significance cutoffs, we generated a list of 188 genes, the expression of which was reproducibly (FDR < 0.05) up-regulated
3-fold, and 495 genes, the expression of which was reproducibly down-regulated
3-fold.
Each gene was annotated based on a complete download of the NetAffx database (Affymetrix; database build date March 1, 2003). The relationship among Affymetrix Gene ID numbers and other gene descriptors such as gene symbols, Unigene clusters, GenBank designations, and so on, is summarized in Table 2A
for clarity. We will discuss each unique Affymetrix probe set as if it were a separate gene although it is possible that two probe sets could detect different splice variant RNAs transcribed from a single Unigene cluster.
|
|
|
B-Regulated Genes.
Given our data suggesting that CPT treatment may induce genes involved in inflammation, we examined whether NF-
B-mediated transcription might be responsible because this factor comprises a central node in cellular responses to both inflammation and DNA damage. We analyzed a list of 56 putative NF-
B-responsive genes (11)
and determined whether CPT-induced changes in expression were observed among the corresponding probe sets. This list of NF-
B-inducible genes was disproportionately changed (3.8-fold overrepresentation; Tables 2B
and 3A
) and up-regulated in comparison with the total list of probe sets (yellow dots in Fig. 1B
; Table 3A
). Moreover, a Students t-Test comparison (p = 4.1 x 10-5) between the group of all probe set FCs (n = 22,283, mean FC = -0.0), and our list of NF-
B annotated probe set FCs (n = 56, mean FC = 1.7) strongly suggested an overall positive directionality in NF-
B-regulated gene levels. Thus, within 8 h, CPT treatment engages an NF-
B-regulated transcriptional network. The importance of the NF-
B pathway in mediating CPT-resistance has been noted (8
, 9
, 11)
; see "Discussion").
|
p53-Regulated Genes.
The duration (8 h) and dosage (10 µM) of CPT treatment were sufficient to induce poly(ADP-ribose) polymerase cleavage, an early marker of apoptosis, among the HeLa cell samples (Fig. 2A)
. Because apoptosis induced by DNA damaging agents is largely p53 dependent (13
, 14)
, we asked whether a p53 response was mobilized in HeLa cells where p53 function is generally considered to be abolished by association with human papillomavirus type E6 (15)
. We detected dramatically increased phosphorylation of p53 on serine 15 after CPT treatment (Fig. 2B)
, as well as higher protein level (data not shown), indicating that topo I blockage in HeLa cells results in greater p53 transactivation and increased protein stability. To determine whether CPT treatment influences expression of putative p53-responsive genes, we examined a list of 31 genes that are regulated by p53 and queried their expression status among the Affymetrix probe sets. Gene expression changes in this list were, overall, positive compared with the total list of probe sets (Table 2B
, Fig. 1C
). Moreover, several of these genes were among the highest-fold induced probe sets such as DTR, DKK1, ATF3, and GADD45 (Table 3B)
. Therefore, it is likely that CPT treatment induces a partial or cell-specific p53 response in HeLa.
|
B (8)
, Bcl-2 family (19, 20, 21)
, and histone deacetylation (22)
. Specifically, we report increased expression of the gene encoding HB-EGF (DTR), which is an EGFR ligand and a survival factor for HeLa cells. The highly specific EGF receptor inhibitor, AG1478 (23
, 24)
, is known to inhibit ATP binding to EGFR and is sufficient to abolish the antiapoptotic effect of HB-EGF at a concentration of 500 nM in HeLa (data not shown). At this concentration, AG1478 exerts an apparent superadditive effect on apoptosis together with CPT (Fig. 3A)
|
| DISCUSSION |
|---|
|
|
|---|
Although p53-mediated increase in Bax (32
, 33)
or PUMA (BBC3; Ref. 34
) level is sometimes regarded as a major mechanism for DNA damage-induced cell death, we did not find significant induction of these genes (Table 3B)
. However, a robust increase in the level of PMAIP1 (Noxa) was detected, suggesting that this BH3-only protein may be involved in apoptosis initiation via its disruption of antiapoptotic Bcl-2 family members and subsequent stimulation of caspase-9 (35)
. Other potentially pro-apoptotic changes in the transcriptome include an increase in level of GADD45 message, which may participate in c-Jun NH2-terminal kinase activation (36)
and an increase in DKK1, which disrupts survival signals mediated by the WNT signaling pathway (37)
. An approximate 2-fold increase in the levels of tumor necrosis factor (TNF)-
and one of its receptors, TNFR2, suggests that caspase-8 and Bid-mediated apoptosis could play an additional role in the induction of CPT-mediated apoptosis. TNF-
or similar ligands are likely to be active after CPT because significant up-regulation of TNFAIP6 is observed.
In addition to apoptosis, the TNFR family of receptors also propagates proinflammatory and prosurvival signals mediated by NF-
B transactivation. Our finding of increase (4.5-fold) in the transcript level of TNFRSF10D (TRAILRA) implies that its corresponding protein may be transmitting signals via NF-
B without the caspase-recruiting TRADD domain that is characteristic of other TNFR family members (38)
. Consistent with an increase in signaling via NF-
B and other inflammatory factors, we observe a 3-fold increase in expression of NF-
B-2/RelB, an almost 6-fold increase in CD14, a nearly 4-fold increase in NFATC1, and a 3.5-fold increase in AOC1 (VAP-1). In addition, we detect increases in levels of NF-
B-inducible genes such as BDKRB1, CXCL1, IL6, IL8, BIRC3 (cIAP2), and PTX3. Although NF-
B is known to be involved in responses to genotoxic stress (39)
, it is likely that this response attenuates the parallel apoptotic response, possibly allowing the cell time to repair the DNA or undergo a checkpoint-mediated arrest.
At least seven NF-
B-inducible genes have been implicated in a survival (or antiapoptotic) pathway (8)
. Of these, only the BIRC3 (c-IAP2) gene is elevated (2.7-fold) after CPT treatment. Expression of BIRC3 has been shown to be sufficient to protect against etoposide-induced apoptosis in HT1080I cells (40)
. Several other genes implicated in NF-
B and/or apoptotic pathways are also up-regulated by CPT exposure. These include NF-
B family member NF-
B-2 (p52/p100) and the NF-
B inhibitor NF-
B-IA (I
B
). In addition, the IRF1 and TNFAIP3 genes are also up-regulated and reported to regulate apoptosis. Thus, our data demonstrate that up-regulation of a large number of genes controlled by NF-
B occurs in response to CPT exposure. Blockage of the NF-
B response using a dominant negative form of the NF-
B inhibitor IkB
delivered with an adenovirus induces rapid apoptotic tumor regression in mice treated with the CPT derivative CPT-11 (9)
. This finding serves to validate the paradigm of examining CPT-responsive genes as potential targets for adjuvant therapy.
In addition to p53 and NF-
B, other transcriptional responses to double-strand DNA breaks are likely to occur via MAFF, ATF3, LHX6, and GBX2 because their message levels are dramatically increased after 8 h of CPT treatment in HeLa cells. Interestingly, ATF3 is believed to repress transcription at the DNA binding elements of both cAMP-responsive element binding protein (41)
and p53 (42)
, suggestive of a possible autoregulatory circuit involving p53. Additional complex regulatory networks may be engaged because potent induction of DTR (HB-EGF), dramatic decline in the level of RasGAP, and increases in the levels of immediate early genes Fos1L and Jun (AP-1) suggest that signaling through mitogen-activated protein kinase family pathways may be elevated. Prolonged mitogen-activated protein kinase activation, in contrast to the transient stimulations that produce cell cycle progression, could generate the opposite effect with cell cycle arrest and/or differentiation, as observed in PC12 (43
, 44)
and other epithelial cells (45)
. In agreement with this idea, we observe highly significant induction of RIS1 (Table 1C)
, which correlates with Ras-mediated senescence but not DNA damage in fibroblasts (46)
.
Recently, it was demonstrated that a chromatin remodeling complex could replace H2A with an H2A variant called H2A.Z in vitro and in vivo (47) , suggesting that the variant histones might be incorporated into damaged DNA. In addition, a homologue of Saccharomyces cerevisiae Sir2, SIRT4, is also highly elevated in expression after CPT treatment. Sir2 is a NAD-dependent histone deacetylase that silences transcription in telomeric regions (48) . We have previously shown that posttranslational modifications of the nucleosome are required for repair of DNA damage caused by CPT (12) .
If changes in the transcriptome after CPT treatment correlate with protein levels of the corresponding genes, the discovery of specific, CPT-inducible drug targets by microarray technology may enable adjunctive treatments that are beneficial to patients who are administered topo I inhibitors. Our finding of increased DTR (HB-EGF) level suggests that its gene product may activate its cognate receptor, ErbB1/EGFR, and generate downstream antiapoptotic signals via diverging pathways that could be Akt/PKB-dependent (49
, 50)
and/or independent (50
, 51)
. This possible mechanism is consistent with a previous study from the Aaronson group, which demonstrated that DNA damage (mitomycin C) directed at certain cell lines causes p53-dependent induction of HB-EGF mRNA and protein level; this autocrine factor, in turn, seems to account for a significant fraction of Akt/PKB and mitogen-activated protein kinase activity induced by artificial p53 expression (50)
. Direct blockage of contrary EGFR-mediated survival signals might therefore increase the predominantly apoptotic effect of CPT in HeLa. To test this, we examined whether DTR (HB-EGF) is sufficient to protect HeLa cells against CPT-mediated apoptosis and determined that a highly specific inhibitor of the EGFR (AG1478) potentiates it. Currently, strategies to block growth and survival signals emanating from receptors of the ErbB family (Herceptin, C225, and Iressa) are under clinical investigation (16
, 52)
. Similarly, the robust CPT-mediated increase in SAT transcript level led us to speculate that alteration of polyamine metabolism would be effective at increasing CPT-mediated apoptosis. The effect of polyamine depletion on apoptosis may differ dramatically depending on the cell type and nature of insult (53, 54, 55)
, and thus, additional experiments to elucidate whether pharmacological inhibition of synthesisor catabolismof spermine/spermidine via agents such as difluoromethylornithine may be justified. In addition, our finding of dramatic induction of certain variant histones (Fig. 1)
prompts speculation that small molecule inhibitors of histone modification may prove effective in augmenting a response to topo I directed chemotherapy, perhaps by modulation of the proposed apoptotic histone code (56)
.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Requests for reprints: Michael F. Christman. Phone: (617) 414-1636; Fax: (617) 414-1646; E-mail: mfc{at}bu.edu; Internet address: http://www.gg.bu.edu.
1 Internet address: http://www.scioncorp.com. ![]()
2 Internet address: http://www.ncbi.nlm.nih.gov/geo/. ![]()
3 Internet address: http://www6.unito.it/genex/cybert/. ![]()
4 Internet address: http://david.niaid.nih.gov/david/ease.htm. ![]()
Received 7/ 8/03. Revised 1/14/04. Accepted 1/15/04.
| REFERENCES |
|---|
|
|
|---|
B as a therapeutic target in cancer. Trends Mol Med, 8: 385-9, 2002.[CrossRef][Medline]
B. Nat Med, 5: 412-7, 1999.[CrossRef][Medline]
B/Rel family of transcription factors in oncogenic transformation and apoptosis. Mutat Res, 437: 231-43, 1999.[CrossRef][Medline]
-Difluoromethylornithine as treatment for metastatic breast cancer patients. Clin Cancer Res, 5: 3438-44, 1999.
B and protects against TRAIL-mediated apoptosis, yet retains an incomplete death domain. Immunity, 7: 813-20, 1997.[CrossRef][Medline]
B: control of oncogenesis and cancer therapy resistance. Biochim Biophys Acta, 1470: M55-62, 2000.[Medline]
B antiapoptosis: induction of TRAF1 and TRAF2 and c-IAP1 and c-IAP2 to suppress caspase-8 activation. Science (Wash. DC), 281: 1680-3, 1998.This article has been cited by other articles:
![]() |
Y. Gong, J. Yue, X. Wu, X. Wang, J. Wen, L. Lu, X. Peng, B. Qiang, and J. Yuan NSPc1 is a cell growth regulator that acts as a transcriptional repressor of p21Waf1/Cip1 via the RARE element Nucleic Acids Res., December 4, 2006; 34(21): 6158 - 6169. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. A. Byers, J. P. Price, J. J. Cooper, Q. Li, and D. H. Price HEXIM2, a HEXIM1-related Protein, Regulates Positive Transcription Elongation Factor b through Association with 7SK J. Biol. Chem., April 22, 2005; 280(16): 16360 - 16367. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |