Molecular Cancer Therapeutics Funding
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

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Olivero, M.
Right arrow Articles by Di Renzo, M. F.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Olivero, M.
Right arrow Articles by Di Renzo, M. F.
Mol Cancer Ther. 2006;5:1126-1135
© 2006 American Association for Cancer Research

Research Articles: Targets

Genes regulated by hepatocyte growth factor as targets to sensitize ovarian cancer cells to cisplatin

Martina Olivero1, Tina Ruggiero1, Silvia Saviozzi3, Andrea Rasola2, Nadia Coltella1, Stefania Crispi4, Ferdinando Di Cunto5, Raffaele Calogero3 and Maria Flavia Di Renzo1

1 Laboratory of Cancer Genetics and 2 Division of Molecular Oncology, Institute for Cancer Research and Treatment, University of Torino School of Medicine, Candiolo (Turin), Italy; 3 Department of Clinical and Biological Sciences, University of Torino, Orbassano (Turin), Italy; 4 Gene Expression Core-Human Molecular Genetics Laboratory, Institute of Genetics and Biophysics "A.B.T.," Consiglio Nazionale delle Ricerche, Naples, Italy; and 5 Department of Genetics, Biology, and Biochemistry, University of Torino, Turin, Italy

Requests for reprints: Maria Flavia Di Renzo, Laboratory of Cancer Genetics, Institute for Cancer Research and Treatment, University of Torino School of Medicine, SP 142, KM 3.95, 10060, Candiolo (Torino), Italy. Phone: 39-11-993-3343; Fax: 39-11-993-3524. E-mail: mariaflavia.direnzo{at}ircc.it


    Abstract
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Advanced ovarian cancers are initially responsive to chemotherapy with platinum drugs but develop drug resistance in most cases. We showed recently that hepatocyte growth factor (HGF) enhances death of human ovarian cancer cell lines treated with cisplatin (CDDP) and that this effect is mediated by the p38 mitogen-activated protein kinase. In this work, we integrated genome-wide expression profiling, in silico data survey, and functional assays to identify transcripts regulated in SK-OV-3 ovarian cancer cells made more responsive to CDDP by HGF. Using oligonucleotide microarrays, we found that HGF pretreatment changes the transcriptional response to CDDP. Quantitative reverse transcription-PCR not only validated all the 15 most differentially expressed genes but also confirmed that they were primarily modulated by the combined treatment with HGF and CDDP and reversed by suppressing p38 mitogen-activated protein kinase activity. Among the differentially expressed genes, we focused functional analysis on two regulatory subunits of the protein phosphatase 2A, which were down-modulated by HGF plus CDDP. Decrease of each subunit by RNA interference made ovarian cancer cells more responsive to CDDP, mimicking the effect of HGF. In conclusion, we show that HGF and CDDP modulate transcription in ovarian cancer cells and that this transcriptional response is involved in apoptosis regulation. We also provide the proof-of-concept that the identified genes might be targeted to either increase the efficacy of chemotherapeutics or revert chemotherapy resistance. [Mol Cancer Ther 2006;5(5):1126–35]


    Introduction
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Drug resistance remains a major clinical challenge for cancer treatment. Patients suffering from advanced ovarian carcinoma are in most cases initially responsive to the combined cytoreductive and chemotherapy treatment. However, they later experience disease relapse due to eventual tumor recurrence and emergence of drug-resistant tumor cells. There are several mechanisms by which tumor cells resist to cytotoxic agents (1, 2). One is related to resistance to apoptosis, a type of cell death that is triggered in response to common chemotherapy and radiotherapy regimens. Evasion from apoptosis is also one of the fundamental hallmarks of cancer (3). A corollary of this is that if the apoptotic pathway inactivated in tumor development is the same as, or overlaps, that leading to cell death by drug, most cancers would be expected to be resistant to drug-induced apoptosis from their onset. Therefore, therapeutic interventions that can lower the threshold for apoptosis of tumor cells could become useful approaches to treat cancer when used either as a single agent or in combination with other therapeutic modalities (4).

In ovarian cancer, apoptosis is impaired by the antiapoptotic action of several oncogenes, which include growth factor receptors and other kinases that transduce signals from the membrane to the nucleus (57). We have shown recently (8) that the hepatocyte growth factor (HGF) sensitizes several ovarian carcinoma cell lines to cisplatin (CDDP) and Taxol, commonly used as first-line chemotherapy in advanced ovarian cancer. HGF increases cancer cell death at very low drug doses. This was a novel and surprising finding. HGF elicits a distinctive biological program known as "invasive growth" by orchestrating cell survival, proliferation, and motility through activation of its receptor, encoded by the MET oncogene, and of several downstream signaling pathways, such as the extracellular signal-regulated kinase 1/2 mitogen-activated protein kinases (MAPK) and phosphatidylinositol 3-kinase/AKT and the AKT substrate mammalian target of rapamycin (9, 10). Nevertheless, we have shown that activation of these survival pathways does not hinder the ability of HGF to enhance drug-dependent apoptosis in ovarian cancer cells. In these cells, HGF simultaneously activates the p38 MAPK, which is further activated by drugs and leads to cell death (11).

As we have shown that in ovarian cancer cell lines HGF-dependent sensitization to drugs required long-term exposure (8), we hypothesized that HGF effect could be coupled to transcriptional regulation of apoptosis-related genes. To explore this possibility, we have studied the transcriptional targets of HGF and CDDP in an ovarian cancer cell line. The identified molecules were sought as regulators of ovarian cancer cell apoptotic death and can be prospectively viewed as targets for therapy. Here, we show the identification of transcripts modulated in cells committed to CDDP-induced apoptosis by HGF. Among them, we found subunits of protein phosphatase 2A (PP2A) that we targeted to enhance ovarian cancer cell response to the drug.


    Materials and Methods
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Flow Cytometry Analysis of Apoptosis Induction
Experiments were done on SK-OV-3 ovarian carcinoma cells, which showed resistance to platinum drugs (12). Exponentially growing cells were cultured for 48 hours in the presence of pure recombinant HGF (R&D Systems, Minneapolis, MN) at the concentration of 50 ng/mL or control medium. Apoptosis was then induced by adding fresh medium, with or without HGF, supplemented with 10 µmol/L CDDP. Flow cytometry recordings of several independent apoptotic changes were done by a single-tube analysis as described (8, 13).

Microarray Sample Preparation
Total RNA was extracted and purified from SK-OV-3 cell lines at each time point using the Concert Cytoplasmic RNA Purification Reagent (Invitrogen, Carlsbad, CA) as suggested by the manufacturer. RNAs were then quantified and inspected by Bioanalyzer (Agilent Technologies, Waldbrom, Germany) analysis. cRNAs were generated and hybridized on 16 HGU133a Affymetrix (Santa Clara, CA) DNA chips according to the Affymetrix protocol. The chips were scanned with a specific scanner (Affymetrix) to generate digitized image data files.

Microarray Data Analysis
Two of the time-course experiments, done independently (TR5 and TR7), were analyzed. From digitized image data files of raw data, background-normalized image data (CEL files) were generated. Microarray quality control and statistical validation was done using Bioconductor (14). The presence of hybridization/construction artifacts was evaluated with the fitPLM function (Bioconductor package affyPLM). The probes (PM) intensity distribution was evaluated using hist function (Bioconductor package affy). Only two arrays, CDDP 6 hours and CDDP/HGF 6 hours in TR7 experiment were discarded due to low quality. TR5 and TR7 experiments were then analyzed separately. Probe set intensities were obtained by robust multiarray average (15), and normalization was done according to quantiles method (16). The HGU133a 22283 probe sets in TR5 and TR7 were filtered separately to have an interquantile range for each probe set greater than 0.4 (17). This filtering yielded 3,893 probe sets in common between the two experiments. TR5 and TR7, represented as ratio with respect to probe sets mean across the experiments, were combined in a unique exprSet object. The identification of differentially expressed genes associated to the HGF effect was addressed using linear modeling approach (18). For assessing differential expression, an empirical Bayes method (18) was used to moderate the SE of the estimated log-fold changes together with the BH correction (19) of the false discovery rate. Four hundred six probe sets (334 Entrez Genes) were identified as differentially expressed as consequence of HGF effect (adjusted P ≤ 0.05). For each of the two experiments, fold change variation versus the time point 0 of the CDDP treatment was calculated followed by the calculation of fold change variation depending on HGF pretreatment (e.g., HGF/CDDP 6 hours – CDDP 6 hours).

Principal component analysis and hierarchical clustering (ST, Euclidean distance, average clustering, and 5,000 jackknife resampling steps) were done using TMEV 3.1 software.6

Coexpression Analysis
The coexpression analysis was done essentially as described (20). Briefly, we downloaded all the human cDNA microarray experiments deposited in the Stanford Microarray Database (4,129 experiments on May 2005) and identified the probes corresponding to the 15 top ranked differentially expressed genes. For each probe, we then calculated the Pearson correlation coefficient with each probe in the database normalized for the number of common experiments. The normalized Pearson coefficient was then used to rank all the probes in the database, and the top 1% of the ranked database was selected as the list of coexpressed probes. We chose a 1% cutoff as we showed previously that on average most of the biologically significant coexpression is concentrated in this interval (20). We thus obtained a total of 89 coexpression lists that were merged to evaluate the recurrence of probes representing the 334 differentially expressed genes. The observed frequency was then compared with that expected by chance, and an overrepresentation P was calculated ({chi}2 test).

Quantitative Reverse Transcription-PCR
cDNA was synthesized from 1 µg RNA using MMLV-RT(H–) enzyme (Promega, Madison, WI). Quantitative reverse transcription-PCR (RT-PCR) was done on an ABI PRISM 7900HT Sequence Detection System (PE Biosystems, Foster City, CA) in 384-well plates assembled by Biorobot 8000 (Qiagen, Hilden, Germany) using a final volume of 20 µL. All quantitative RT-PCR (qPCR) mixtures contained 20 ng retrotranscribed RNA, 1x SYBR Green PCR Master Mix (2x, Applera, Foster City, CA), and 150 µmol/L of each target-specific primer. The primer sequences are reported in Supplementary Table S1B.7 The expression of each target genes was evaluated using a relative quantification approach (–{Delta}{Delta}CT method; ref. 21) with cyclophilin A as internal reference. The consistency of cyclophilin A (NM_021130) was assessed in a study aimed at defining the most stable expressed genes of seven putative reference genes. This study consisted of qPCR of the seven genes on the same samples (data not shown) followed by analysis with the genNorm software (22).

Short Hairpin RNA Stable Expression
SK-OV-3 and amphotropic phoenix retroviral producer cells were both from American Type Culture Collection (Manassas, VA). For production of retroviral stocks, phoenix cells were seeded at 4 x 106 per 10-cm dish and transfected using CaPO4 precipitation of 15 µg pRetroSuper-short hairpin RNA (shRNA) plasmid DNA. The following sequences were subcloned: PPP2R1B (NM_002716) AGGACCCGAAGTGAATTGT and PPP2R5E (NM_006246) CCAGTTATTAGGGGGTTAA. Forty-eight hours after infection, cells were selected using 1 µg/mL puromycin; typically 95% to 100% of the cells were infected and drug resistant.


    Results
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
HGF Sensitizes Ovarian Cancer Cells to CDDP-Dependent Apoptosis via the p38 MAPK
In response to CDDP, the SK-OV-3 ovarian cancer cells enter the apoptotic program as assessed using a multiparametric cytofluorimetric assay (8, 13), which simultaneously measures hallmarks of early, intermediate, and late phases of apoptosis: decrease of mitochondrial membrane potential, cell shrinkage, increase in granularity, exposure of phosphatidylserine at cell surface, and plasma membrane permeabilization to propidium iodide. We showed previously that SK-OV-3 cells express the receptor for the HGF encoded by the MET oncogene, which is activated by HGF (8, 23). As shown in Fig. 1A , a 48-hour long pretreatment with 50 ng/mL HGF caused a marked sensitization to the subsequent treatment with 10 µmol/L CDDP, with a net increase of the measured apoptotic features after 48 to 72 hours. Notably, this sensitization required pretreatment with HGF, as its coincubation with the drug was ineffective (Fig. 1A). These data confirmed what we observed previously (8) and established the time course of the increased death of SK-OV-3 cells caused by the combined treatment with HGF and CDDP.


Figure 1
View larger version (55K):
[in this window]
[in a new window]
 
Figure 1. HGF synergizes with CDDP in activating apoptosis through p38 MAPK in ovarian cancer cells. Apoptosis of ovarian carcinoma cells (SK-OV-3) was measured with a multiparametric fluorescence-activated cell sorting analysis. Bottom, cells displaying mitochondrial depolarization (reduced tetramethylrhodamine methyl ester staining); right quadrant, cells exposing phosphatidylserine on their surface (increased Annexin V-FITC staining). Numbers, percentages of healthy cells (H, delimited in the quadrant) and dead, propidium iodide–permeable cells (D, not shown in the diagrams). Cells undergoing apoptosis and thus outside the quadrant should be added to the healthy cells to give 100%. A, time course (hours) of the cell apoptotic response to CDDP after 48 h of treatment with recombinant purified HGF (50 ng/mL). Top right, response of cells simultaneously exposed to CDDP and HGF and examined after 72 h (+HGF). Representative experiment of the seven (TR1–TR7) that measured the apoptosis achievement by cells from which RNA was extracted. B, apoptotic response to CDDP treatment (72 h) with or without HGF pretreatment (50 ng/mL for 48 h) of cells stably carrying and expressing the DN-p38 MAPK transgene. Representative experiment of the five that measured the apoptosis achievement by cells from which RNA was extracted. In all experiments, SK-OV-3 cells were exposed to 10 µmol/L CDDP.

 
We have also observed previously that either HGF or CDDP, alone or in combination, activate the p38 MAPK (11), which is activated by cell stress and is involved in cancer cell apoptosis induced by several chemotherapeutics, including platinum drugs and taxanes (24, 25).

We stably expressed a dominant-negative form of p38 MAPK (DN-p38) using a lentiviral vector carrying the double T180A/Y182F Flag-tagged p38 MAPK mutant (11). Lentiviral vectors drive the stable transgene random integration in cell genomic DNA, thus allowing repetitive studies of the bulk unselected cell population. We transduced SK-OV-3 and assessed using Western blot and fluorescence-activated cell sorting analysis with anti-Flag antibodies (data not shown) that ~90% of cells expressed the DN-p38. We showed previously that the DN-p38 MAPK abolished the wild-type p38 MAPK activation without interfering with other HGF-dependent signaling pathways (11). Remarkably, we found that the inhibition of the p38 MAPK pathway almost abolished apoptosis elicited by CDDP and enhanced by HGF in SK-OV-3 cells (Fig. 1B). These data indicate that p38 MAPK activation is required and has a central role in the induction and modulation of apoptosis driven by CDDP and enhanced by HGF in ovarian cancer cells. Thus, SK-OV-3 cells stably expressing the DN-p38 have been used as a model to establish a functional link between cell transcriptional response to HGF and CDDP and induction of apoptosis (see below).

HGF Modulates Transcription in Ovarian Cancer Cells Committed to CDDP-Dependent Apoptosis
As pretreatment is required to obtain maximal sensitization of ovarian carcinoma cells to CDDP, we assumed that HGF mainly enhances the apoptotic response by modulating gene transcription. Thus, we studied the expression profiles of ovarian carcinoma cells treated with HGF and CDDP alone or in combination. To perform genome-wide profiling, we first used microarray experiments and then validated the differentially expressed genes using real-time qPCR. As shown in Fig. 1, in both cells treated with CDDP alone or with CDDP and HGF, even the earliest apoptotic changes, such as the decrease of mitochondrial membrane potential, become manifest not earlier than after 48 hours of CDDP exposure (Fig. 1A). Therefore, time-course experiments were done, pretreating SK-OV-3 cells with HGF, exposing HGF-treated cells to CDDP, and preparing RNA after 6, 12, and 24 hours. Control cells were either exposed to control treatment or exposed to HGF alone or CDDP alone for 6, 12, and 24 hours. We analyzed the expression profiles of cells from three independent experiments after confirming apoptosis achievement by a sample aliquot of cells kept in culture for 48 and 72 hours; a representative control experiment is shown in Fig. 1A. Each of the three time-course experiments, called TR5, TR6, and TR7, included eight experimental points as follows: cells treated with control medium for 48 hours (0-hour CDDP), cells treated with HGF for 48 hours (HGF/0-hour CDDP), cells pretreated with control medium and thereafter with CDDP for 6 hours (6-hour CDDP), 12 hours (12-hour CDDP), and 24 hours (24-hour CDDP), and cells pretreated with HGF and thereafter with CDDP for 6 hours (HGF/6-hour CDDP), 12 hours (HGF/12-hour CDDP), and 24 hours (HGF/24-hour CDDP). mRNAs from cells at each of the above points were analyzed. Two of the three time-course experiments (TR5 and TR7) were analyzed using HGU133a Affymetrix oligonucleotide microarrays. The third experiment (TR6) was added for the further validation of results with qPCR (see below).

Principal component analysis (26) was used to investigate the overall effects of HGF and CDDP in TR5 and TR7 experiments. This technique represents the global gene expression changes between samples as coordinates in a multidimensional space in which the axes are called the "principal components." In TR5 and TR7 experiments, the first and second components account for the greatest part of the variability within the data set (73.5%); thus, the graphical representation of the first two components (Fig. 2 ) is sufficient to show the differences between samples. Figure 2 shows that HGF pretreatment above all strongly influenced the overall gene expression variation, as experiments simply clustered based on this factor (above the dashed line in Fig. 2). Conversely, the effect of CDDP was time dependent, although the two experiments were only partially synchronized.


Figure 2
View larger version (10K):
[in this window]
[in a new window]
 
Figure 2. HGF above all regulates transcription in ovarian cancer cells sensitized to CDDP-dependent apoptosis. Microarray experiments are analyzed using principal component analysis that is a data reduction method to identify patterns in a data matrix and reflects the variance within the experimental data. Principal component analysis was done on overall expression data of cells either not treated or treated with HGF for 48 h before exposure to CDDP for the indicated hours (0, 6, 12, and 24 h). Cells not treated with HGF are indicated by closed symbols, whereas cells treated with HGF are indicated by open symbols. All eight time point conditions of the TR5 experiment (diamonds) were examined, whereas only six of these conditions of the TR7 experiment (circles) were evaluated, as the 6-h treatment with CDDP was not taken into consideration (see Materials and Methods). The first two principal components, PC1 and PC2, account for 58% and 15.5% of the total variance, respectively. HGF effect is visible in both experiments, allowing the separation of the plus HGF (+HGF) with respect to the minus (-HGF) condition, positioned above and below the dashed line, respectively. The effect of increasing length of CDDP treatment is visible as projection on PC1. Trends are similar within the two experiments, but they are only partially synchronized.

 
Therefore, we focused the subsequent statistical analysis on the identification of probe sets differentially expressed in cells undergoing CDDP-dependent apoptosis as a consequence of the HGF pretreatment. Using linear modeling analysis (18), we identified 406 probe sets associated to the HGF pretreatment. Unsupervised hierarchical clustering of these 406 probe sets (Fig. 3A ), corresponding to 334 genes, again showed HGF effect on CDDP-regulated transcripts in both TR5 and TR7 experiments. Subclustering related to the timing of CDDP effect can be also observed (Fig. 3A), more evidently in one cluster (blue bar in Fig. 3A).


Figure 3
View larger version (45K):
[in this window]
[in a new window]
 
Figure 3. HGF regulates the differential expression of 334 genes in ovarian cancer cells sensitized to CDDP-dependent apoptosis. Microarray data clustering shows the HGF- and time-dependent transcriptional response of ovarian cancer cells undergoing apoptosis in response to CDDP. SK-OV-3 cells were either not treated or treated with HGF for 48 h before exposure to CDDP for the indicated hours (0, 6, 12, and 24 h). Cells not treated with HGF are indicated as CDDP, whereas cells pretreated with HGF are indicated as HGF/CDDP. All eight experimental points of the TR5 experiment were examined, whereas only six of these conditions of the TR7 experiment were evaluated, as the 6-h treatment with CDDP was not taken into consideration (see Materials and Methods). A, gene-oriented unsupervised hierarchical clustering of the 406 probe sets, corresponding to 334 genes, associated to the HGF effect. Before clustering, gene expression data were mean centered, replacing each value by the value subtracted of the mean of the row that value belongs. B, two-dimensional unsupervised hierarchical clustering of the top ranked 18 probe sets, corresponding to 15 genes, further characterized in qPCR experiments. This nomenclature is a step forward from the National Center for Biotechnology Information LocusLink, with both a major increase in taxonomic scope and an improved access through the many tools associated with National Center for Biotechnology Information Entrez (49). E.G., Entrez Gene.

 
For qPCR validation, we selected the 18 probe sets, of 406, which were characterized by the largest differential expression (>2-fold change in both TR5 and TR7 experiments; Fig. 3B). The 15 genes corresponding to these probe sets are listed in Fig. 3B and described in detail in Supplementary Table S1A. Differential expression of all these 15 transcripts was confirmed by qPCR not only in the TR5 and TR7 experiments but also in the TR6 experiment similarly carried out (see above) but not analyzed with microarrays. The qPCR validation of each of the 15 genes is shown in Fig. 4 and Supplementary Figs. S1 and S2.7 It is noteworthy that all the 15 top ranked differentially expressed transcripts showed a quantitative regulation, which paralleled that found in microarrays. However, by comparing the fold change in expression measured using microarray analysis (Fig. 3B) and that measured using qPCR (calculated by the equation {Delta}{Delta}CT; Fig. 4; Supplementary Figs. S1 and S2),7 the differential expression came out to be two to four times higher using qPCR probably due to the higher sensitivity of the latter technique. qPCR also confirmed the consistency of transcript regulation in all the three time-course experiments (Fig. 4; Supplementary Figs. S1 and S2).7


Figure 4
View larger version (17K):
[in this window]
[in a new window]
 
Figure 4. p38 MAPK activity mediates gene expression regulation by HGF in ovarian cancer cells undergoing CDDP-dependent apoptosis. qPCR validation of 6 of the 15 most differentially expressed genes identified with microarray analysis (the other 9 genes are shown in Supplementary Figs. S1 and S2; available online at http://mct.aacrjournals.org/). Expression variations of the gene were measured in wild-type (wt) SK-OV-3 cells pretreated for 48 h either with pure recombinant HGF ({blacksquare}) or with control medium ({square}) and then cultured with 10 µmol/L CDDP for the indicated times (0, 6, 12, and 24 h). Transcripts were also measured in SK-OV-3 cells stably expressing the DN-p38 MAPK transgene (DN-p38) and similarly treated. {triangleup}, control medium; {blacktriangleup}, +HGF. Points, mean expression variations measured at each time point using mRNAs prepared from three experiments carried out independently; bars, SD. The mean fold change in expression of the target gene in treated cells at each time point versus untreated cells at time 0 was calculated using the formula: –{Delta}CT = –[(CT,targetCT,cyclophilin A)time x – (CT,targetCT,cyclophilin A)time 0].

 
Furthermore, qPCR allowed a better classification of the 15 top ranked differentially expressed genes. Some were mainly regulated by the combined treatment with HGF and CDDP (e.g., RAMP and RRM2; Fig. 4). Others were already modulated by HGF alone (e.g., MAT2A and EIF5; Fig. 4); in the latter case, in fact, consistent differences in expression were already measured after 48 hours HGF treatment (i.e., before the addition of CDDP; 0-hour point in Fig. 4). Using qPCR, we confirmed that 12 of 15 top ranked genes were down-modulated by CDDP plus HGF (Figs. 3 and 4; Supplementary Figs. S1 and S2).7 Two of the 15 (CGB and GMDS; Fig. 4; Supplementary Fig. S1)7 were down-modulated by CDDP but reversed by HGF. Only 1 of 15 genes (KRT14; Fig. 4) was consistently and strongly up-regulated by both HGF alone and HGF plus CDDP. Supplementary Table S27 lists all the genes regulated >2-fold by HGF alone after 48 hours of treatment. This list included the above-mentioned MAT2A and EIF5 and the other genes (KRT14, ASNS, CGB, FBX09, SLC35E2, LCMT2, and PHDLA1), which we validated as regulated by HGF in qPCR (Fig. 4; Supplementary Figs. S1 and S2).7

Transcriptional Regulation by HGF and CDDP Relies on p38 MAPK Activation
As p38 MAPK regulates several transcription factors (2730), we asked if HGF and CDDP transcription modulation requires p38 MAPK activation, which mediates their apoptotic effect. The bulk unselected SK-OV-3 cell population, stably expressing the DN-p38, was treated with either CDDP alone or CDDP and HGF, and the 15 top ranked transcripts were quantified using qPCR at each time point after treatment. The actual apoptotic response in the time-course experiments was confirmed on a sample aliquot of cells kept in culture for 72 hours (Fig. 1B). As shown in Fig. 4 and Supplementary Figs. S1 and S2,7 all the 15 transcripts were differently regulated in cells where p38 MAPK was inactive. In most cases, inactivation of the p38 MAPK reverted the transcript regulation by HGF and CDDP. These data show that the HGF-dependent transcriptional regulation of ovarian cancer cells committed to CDDP-induced apoptosis mostly relies on p38 MAPK activation.

Identification of Genes Coregulated with the Transcriptional Targets of HGF and CDDP
The finding that transcription of 15 genes was consistently modulated by HGF and CDDP and reverted by DN-p38 raised the interesting possibility that they could be part of a wider transcriptional response functionally related to p38-dependent apoptosis. Furthermore, we noticed that the degree of differential expression of the 15 genes assessed using microarray analysis was remarkably lower than that measured with qPCR. Altogether, data suggested that a wider transcriptional response could include other genes possibly underestimated by the microarray analysis. If this was the case, the genes belonging to the group of the 15 top ranked differentially expressed genes, as well as many of the other genes modulated to a lesser extent by the HGF pretreatment, would be expected to undergo significant coregulation under unrelated experimental conditions.

To test this hypothesis, and to identify other genes possibly belonging to the same transcriptional module, we decided to perform a coexpression analysis using a collection of 4,129 cDNA microarray experiments deposited in the Stanford Microarray Database.

In particular, for each of the above 15 genes, we obtained coexpression ranked lists containing the 1% probes of the database showing the highest expression profile similarity (see Materials and Methods; see ref. 20). Consistent with our hypothesis, in these lists, the recurrence of the 15 top ranked genes was much higher than expected by chance (cumulative P = 2.7 x 10–6, {chi}2 test), thus ruling out that the expression regulation observed in our settings was merely due to cell death.

We then calculated an overrepresentation P for each gene in the obtained coexpression lists and used this number as an independent filter to select those genes that, among the 334 differentially expressed genes, were above all recurrent (threshold P < 0.01). The selected 94 genes, listed in Supplementary Table S3,7 are the most likely to be important in the functional response induced by the HGF and CDDP combined treatment.

Functional Evaluation of Genes Differentially Expressed in Ovarian Cancer Cells Committed to Apoptosis
We reasoned that genes regulated by HGF and CDDP in cells undergoing apoptosis could be directly targeted to modulate ovarian cancer cell response to CDDP. Among the 94 differentially expressed genes that are also recurrently coexpressed with the top ranked 15 (Supplementary Table S3),7 we focused on two regulatory subunits (PPP2R1B and PPP2R5E) of the PP2A (described in Supplementary Table S1A).7 Indeed, their unexpected role as survival phosphatases had been also suggested by MacKeigan et al. (31), who carried out a genome-wide screening of kinases and phosphatases using a small interfering RNA library. qPCR analysis confirmed that these PP2A subunits were down-modulated by HGF and CDDP and that down-modulation was reversed in cells expressing DN-p38 (Supplementary Fig. S2).7

Therefore, we investigated if down-modulation of either the PPP2R1B or the PPP2R5E subunit of the PP2A was alone able to commit ovarian cancer cells to CDDP-mediated cell death. shRNAs specific for either PPP2R1B or PPP2R5E were stably expressed in SK-OV-3 cells by retroviral vectors. After selection and propagation in culture, the bulk population showed a mean 60% and 70% reduction of the PPP2R5E and PPP2R1B subunit transcript, respectively, as assessed with qPCR (data not shown). Figure 5 shows that lessening this subunit expression did not induce cell death per se but sensitized ovarian cancer cells to CDDP-mediated apoptosis. Data show that cells expressing the shRNAs specific to these PP2A subunits are stably ready to respond to CDDP. Remarkably, phosphatase subunit knocking down increased the number of cells in each apoptotic phase, either early (cells displaying depolarized mitochondria), intermediate (Annexin V–positive cells), or late (propidium iodide–positive cells), although the intensity of Annexin V staining was lowered with respect to control conditions. The latter data indicate that apoptosis signaling was boosted from its very early steps, but the propagation of plasma membrane lipid scrambling could not be enhanced by the absence of phosphatase activity.


Figure 5
View larger version (34K):
[in this window]
[in a new window]
 
Figure 5. Silencing of PP2A subunit genes sensitizes ovarian cancer cells to CDDP. Apoptosis of ovarian carcinoma cells (SK-OV-3) was measured with a multiparametric fluorescence-activated cell sorting analysis as described in Fig. 1 legend. Bottom, cells displaying mitochondrial depolarization (reduced tetramethylrhodamine methyl ester staining); right quadrant, cells exposing phosphatidylserine on their surface (increased Annexin V-FITC staining). Numbers, percentages of healthy cells (H, delimited in the quadrant) and dead, propidium iodide–permeable cells (D, not shown in the diagrams). Wild-type cells or cells stably expressing interfering shRNA specific for either the PPP2R1B or the PPP2R5E subunits of the PP2A or a control unrelated shRNA were treated for 72 h with 10 µmol/L CDDP. Representative experiments of the four.

 

    Discussion
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 
In this work, we combined different approaches to identify genes whose expression mediate the ability of HGF to sensitize ovarian cancer cells to drug-induced apoptosis. The aim was to find genes whose modulation was per se able to predispose cells to better respond to drugs. We here provide proof-of-concepts that we have been able to identify genes that might be targeted to increase the response of ovarian cancer cells to CDDP.

We combined in vitro and in silico analyses to identify the transcripts modulated in cells committed by HGF to respond to CDDP and relying on the p38 MAPK activation as we have shown previously that p38 MAPK mediates HGF and CDDP combined effects (11). The association between each of the most differentially expressed genes and cell death regulation is confirmed by the description reported in Supplementary Table S1A.7

Genes modulated only in response to the combined treatment with HGF and CDDP, but not by the single treatment, are particularly interesting as they could be better related to HGF-mediated sensitization to CDDP. In fact, this is a distinctive response of ovarian cancer cells as shown in several ovarian cancer cell lines (8, 11). In this group, those above all down-modulated were PANK3, RAMP, RRM2, and TDG. Most of them have been already associated to cell proliferation and apoptosis and even to cancer cell resistance to chemotherapy as described in detail in Supplementary Table S1A (3235).7 Among the most differentially expressed genes, we identified a second group regulated by HGF alone and further modulated, in most cases, after cell treatment with CDDP. This group includes MAT2A, ASNS, EIF5, FBX09, LCMT2, SLC35E2, and PHDLA1, which have been associated to cell proliferation and survival, as described in Supplementary Table S1A (36, 37).7

Data suggest that each of the 15 top ranked genes is pertinent to apoptosis mediated by HGF and CDDP. In fact, the regulation of their expression was reverted by the DN-p38 MAPK, which abolished the HGF-enhanced apoptotic response to CDDP (11). This establishes a functional link between genes regulated by HGF and ovarian cancer cell sensitization to CDDP and allowed the identification of the regulated genes as suitable functional targets. As expected, because each of the 15 top ranked genes responded differently to HGF and CDDP, the reversion induced by the DN-p38 MAPK resulted in different outcome. This further supports the idea that the 15 top ranked genes belong to different transcriptional modules.

The qPCR validation not only fully confirmed the reliability of microarray experiments but also showed that microarray analysis led to underestimate the degree of expression regulation. Thus, we assumed that some of the 334 genes identified by this analysis as modestly differentially expressed might be concurrently and consistently regulated in cells undergoing death. More importantly, we could have missed genes either linked to more physiologic cellular functions than response to drug or involved in important biological processes or more targetable, such as those associated to enzymatic activities. The coexpression test, set up to use the huge number of microarray experiments available online, helped us to refine the selection of HGF-regulated genes. Among the genes that passed this filter, we found two subunits of the PP2A, which were down-modulated in cells treated with HGF and CDDP in a p38 MAPK-dependent way. Repression of these genes sensitized cells to CDDP, mimicking HGF effect and confirming their role as survival phosphatases, in agreement with the results reported in a wide, small interfering RNA–based, screen of human kinases and phosphatases (31).

PP2A is a family of abundantly expressed serine/threonine phosphatases implicated in a multitude of cellular functions (38). The PP2A core structure consists of a catalytic C and a scaffolding A subunits. The core A/C dimer can recruit a third regulatory B subunit. Each PP2A subunit has at least two isoforms. The differential association of all these subunits gives rise to an extensive subset of oligomeric A/B/C holoenzymes, allowing the fine-tuning of PP2A activity and substrate specificity (39). Several lines of evidence concur in classifying the PP2A complex as a tumor suppressor protein. Conversely, as other authors (31), we found that the PPP2R1B and PPP2R5E subunits are "survival" proteins. This is in line with the discovery that many other phosphatases protect cells from death (31) and with other reports that highlight the complexity of the PP2A-regulated activity (40). The PPP2R1B gene down-modulated by HGF and CDDP encodes the A/PR65 PP2A scaffolding subunit. It has been shown recently (41) that whereas suppression of the endogenous PPP2R1A (the {alpha} isoform of the same subunit) by 50% leads to cell transformation, further suppression resulted in cell cycle arrest and apoptosis. These findings are also consistent with the demonstration that in rat cells (42) and Drosophila S2 cells (43) a minimal level of PP2A scaffolding A subunit is required for cell survival. Furthermore, we also found the down-modulated PPP2R5E gene, which encodes the B56{varepsilon} subunit. As each specific B subunit might change PP2A enzymatic activity or even direct each holoenzyme PP2A complex to distinct intracellular location (44), it is not surprising that the reduction of the B56{varepsilon} subunit alone might sensitize cells to CDDP-induced apoptosis. Interestingly, it has been shown that another noncatalytic PP2A subunit is required to repress apoptosis in murine cells (45).

Here, we describe the transcriptional response of ovarian cancer cells made more susceptible to CDDP-mediated death by pretreatment with HGF. The genes regulated by HGF and the pathways in which they are involved might be envisaged as prospective therapeutic targets. We provide a proof-of-concept, as we made an ovarian cancer cell line more susceptible to CDDP by suppressing the PP2A subunits, which we found down-modulated by HGF and CDDP in that cell line. Therefore, this work proposes a valuable protocol to identify potential targets that should be validated in several additional settings and possibly tailored to ovarian cancer subsets. We found that most of the genes that we considered most likely to be important in the apoptotic response induced by the HGF and CDDP combined treatment are not reported in the signatures that have been proposed to classify either chemoresistant or chemosensitive ovarian cancers (12, 35, 4648). This was not surprising, as the steady-state expression level of the identified genes is not expected to mark ovarian cancer cells sensitivity to CDDP. Nevertheless, HGF-regulated genes are of potential value to modify response to CDDP, the most active chemotherapeutic in ovarian cancer treatment. The potential use of low-dose chemotherapy in chemotherapy-sensitive cancers, like ovarian carcinomas, is important, as lower dosage is easily attainable and less likely to cause toxicity in patients. In addition, in ovarian cancer, chemotherapy resistance is the major cause of secondary failure of conventional chemotherapy. Therefore, the discovery of novel pathways and molecules involved in the apoptotic response to drugs is valuable, provided that they are targetable either to kill cells or to revert chemotherapy resistance.


    Acknowledgments
 
We thank Vincenzo De Sio, Raffaella Albano, and Solange Tienga for technical help, Elena Grassi for the enthusiastic help in data mining, and Dr. Enzo Medico for providing the PPP2R1B and PPP2R5E specific shRNA viruses.


    Footnotes
 
Grant support: Italian Ministry of University and Research, Regione Piemonte, and Italian Association for Cancer Research (M.F. Di Renzo), Ministero Salute and Italian Association for Cancer Research (R. Calogero), and Regione Piemonte Ricerca Scientifica Applicata fellowship (A. Rasola).

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.

Note: The present address for T. Ruggiero is Gene Expression Regulation Laboratory, Experimental Oncology, Istituto Nazionale Ricerca sul Cancro, Genova, Italy. The present address for A. Rasolas is Department of Biomedical Sciences, University of Padova, Padova, Italy.

6 http://www.tigr.org/software/. Back

7 Supplementary materials for this article are available at Molecular Cancer Therapeutics Online (http://mct.aacrjournals.org/). Back

Received 1/ 9/06; revised 3/10/06; accepted 3/24/06.


    References
 Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 

  1. Agarwal R, Kaye SB. Ovarian cancer: strategies for overcoming resistance to chemotherapy. Nat Rev Cancer 2003;3:502–16.[CrossRef][Medline]
  2. Brown JM, Attardi LD. The role of apoptosis in cancer development and treatment response. Nat Rev Cancer 2005;5:231–7.[CrossRef][Medline]
  3. Hanahan D, Weinberg RA. The hallmarks of cancer. Cell 2000;100:57–70.[CrossRef][Medline]
  4. Green DR, Kroemer G. Pharmacological manipulation of cell death: clinical applications in sight? J Clin Invest 2005;115:2610–7.[CrossRef][Medline]
  5. Shayesteh L, Lu Y, Kuo WL, et al. PIK3CA is implicated as an oncogene in ovarian cancer. Nat Genet 1999;21:99–102.[CrossRef][Medline]
  6. Yuan ZQ, Sun M, Feldman RI, et al. Frequent activation of AKT2 and induction of apoptosis by inhibition of phosphoinositide-3-OH kinase/Akt pathway in human ovarian cancer. Oncogene 2000;19:2324–30.[CrossRef][Medline]
  7. Cuello M, Ettenberg SA, Clark AS, et al. Down-regulation of the erbB-2 receptor by trastuzumab (herceptin) enhances tumor necrosis factor-related apoptosis-inducing ligand-mediated apoptosis in breast and ovarian cancer cell lines that overexpress erbB-2. Cancer Res 2001;61:4892–900.[Abstract/Free Full Text]
  8. Rasola A, Anguissola S, Ferrero N, et al. Hepatocyte growth factor sensitizes human ovarian carcinoma cell lines to paclitaxel and cisplatin. Cancer Res 2004;64:1744–50.[Abstract/Free Full Text]
  9. Birchmeier C, Birchmeier W, Gherardi E, Vande Woude GF. Met, metastasis, motility and more. Nat Rev Mol Cell Biol 2003;4:915–25.[CrossRef][Medline]
  10. Trusolino L, Comoglio PM. Scatter-factor and semaphorin receptors: cell signalling for invasive growth. Nat Rev Cancer 2002;2:289–300.[CrossRef][Medline]
  11. Coltella N, Rasola A, Nano E, et al. p38 MAPK turns hepatocyte growth factor to a death signal that commits ovarian cancer cells to chemotherapy-induced apoptosis. Int J Cancer 2006;118:2981–90.[CrossRef][Medline]
  12. Roberts D, Schick, J, Conway S, et al. Identification of genes associated with platinum drug sensitivity and resistance in human ovarian cancer cells. Br J Cancer 2005;92:1149–58.[CrossRef][Medline]
  13. Rasola A, Gramaglia D, Boccaccio C, Comoglio PM. Apoptosis enhancement by the HIV-1 Nef protein. J Immunol 2001;166:81–8.[Abstract/Free Full Text]
  14. Gentleman RC, Carey VJ, Bates DM, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol 2004;5:R80–1.[CrossRef][Medline]
  15. Irizarry RA, Hobbs B, Collin F, et al. Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics 2003;4:249–64.[Abstract]
  16. Bolstad BM, Irizarry RA, Astrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 2003;19:185–93.[Abstract/Free Full Text]
  17. von Heydebreck A, Huber W, Gentleman R. Differential expression with the Bioconductor Project. Bioconductor Project Working Papers 2004; Working Paper 7. Available from: http://www.bepress.com/bioconductor/paper7.
  18. Wettenhall JM, Smyth GK. limmaGUI: a graphical user interface for linear modeling of microarray data. Bioinformatics 2004;20:3705–6.[Abstract/Free Full Text]
  19. Reiner A, Yekutieli D, Benjamini Y. Identifying differentially expressed genes using false discovery rate controlling procedures. Bioinformatics 2003;19:368–75.[Abstract/Free Full Text]
  20. Pellegrino M, Provero P, Silengo L, Di Cunto F. CLOE: identification of putative functional relationships among genes by comparison of expression profiles between two species. BMC Bioinformatics 2004;5:179.[CrossRef][Medline]
  21. 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:402–8.[CrossRef][Medline]
  22. Vandesompele J, De Preter K, Pattyn F, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol 2002;3:RESEARCH0034.[Medline]
  23. Maggiora P, Lorenzato A, Fracchioli S, et al. The RON and MET oncogenes are co-expressed in human ovarian carcinomas and cooperate in activating invasiveness. Exp Cell Res 2003;288:382–9.[CrossRef][Medline]
  24. Losa JH, Parada Cobo C, Viniegra JG, Sanchez-Arevalo Lobo VJ, Ramon y Cajal S, Sanchez-Prieto R. Role of the p38 MAPK pathway in cisplatin-based therapy. Oncogene 2003;22:3998–4006.[CrossRef][Medline]
  25. Olson JM, Hallahan AR. p38 MAP kinase: a convergence point in cancer therapy. Trends Mol Med 2004;10:125–9.[CrossRef][Medline]
  26. Raychaudhuri S, Stuart JM, Altman RB. Principal components analysis to summarize microarray experiments: application to sporulation time series. In: Altman RB, et al, editors. Pacific Symposium on Biocomputing 2000. World Sci. Pub. Co.; 2000, p. 655–6.
  27. Wang XZ, Ron D. Stress-induced phosphorylation and activation of the transcription factor CHOP (GADD153) by p38 MAP kinase. Science 1996;272:1347–9.[Abstract]
  28. Tan Y, Rouse J, Zhang A, Cariati S, Cohen P, Comb MJ. FGF and stress regulate CREB and ATF-1 via a pathway involving p38 MAP kinase and MAPKAP kinase-2. EMBO J 1996;15:4629–42.[Medline]
  29. Morton S, Davis RJ, Cohen P. Signalling pathways involved in multisite phosphorylation of the transcription factor ATF-2. FEBS Lett 2004;572:177–83.[CrossRef][Medline]
  30. Poizat C, Puri PL, Bai Y, Kedes L. Phosphorylation-dependent degradation of p300 by doxorubicin-activated p38 mitogen-activated protein kinase in cardiac cells. Mol Cell Biol 2005;25:2673–87.[Abstract/Free Full Text]
  31. MacKeigan JP, Murphy LO, Blenis J. Sensitized RNAi screen of human kinases and phosphatases identifies new regulators of apoptosis and chemoresistance. Nat Cell Biol 2005;7:591–600.[CrossRef][Medline]
  32. Zhou B, Mo X, Liu X, Qiu W, Yen Y. Human ribonucleotide reductase M2 subunit gene amplification and transcriptional regulation in a homogeneous staining chromosome region responsible for the mechanism of drug resistance. Cytogenet Cell Genet 2001;95:34–42.[CrossRef][Medline]
  33. Zhou BS, Tsai P, Ker R, et al. Overexpression of transfected human ribonucleotide reductase M2 subunit in human cancer cells enhances their invasive potential. Clin Exp Metastasis 1998;16:43–9.[CrossRef][Medline]
  34. Duxbury MS, Ito H, Zinner MJ, Ashley SW, Whang EE. RNA interference targeting the M2 subunit of ribonucleotide reductase enhances pancreatic adenocarcinoma chemosensitivity to gemcitabine. Oncogene 2004;23:1539–48.[CrossRef][Medline]
  35. Macleod K, Mullen P, Sewell J, et al. Altered ErbB receptor signaling and gene expression in cisplatin-resistant ovarian cancer. Cancer Res 2005;65:6789–800.[Abstract/Free Full Text]
  36. Martinez-Chantar ML, Corrales FJ, Martinez-Cruz LA, et al. Spontaneous oxidative stress and liver tumors in mice lacking methionine adenosyltransferase 1A. FASEB J 2002;16:1292–4.[Abstract/Free Full Text]
  37. Paneda C, Gorospe I, Herrera B, Nakamura T, Fabregat I, Varela-Nieto I. Liver cell proliferation requires methionine adenosyltransferase 2A mRNA up-regulation. Hepatology 2002;35:1381–91.[CrossRef][Medline]
  38. Janssens V, Goris J, Van Hoof C. PP2A: the expected tumor suppressor. Curr Opin Genet Dev 2005;15:34–41.[CrossRef][Medline]
  39. Cohen P. The structure and regulation of protein phosphatases. Annu Rev Biochem 1989;58:453–508.[CrossRef][Medline]
  40. Van Hoof C, Goris J. Phosphatases in apoptosis: to be or not to be, PP2A is in the heart of the question. Biochim Biophys Acta 2003;1640:97–104.[Medline]
  41. Chen W, Arroyo JD, Timmons JC, Possemato R, Hahn WC. Cancer-associated PP2A A{alpha} subunits induce functional haploinsufficiency and tumorigenicity. Cancer Res 2005;65:8183–92.[Abstract/Free Full Text]
  42. Strack S, Cribbs JT, Gomez L. Critical role for protein phosphatase 2A heterotrimers in mammalian cell survival. J Biol Chem 2004;279:47732–9.[Abstract/Free Full Text]
  43. Li X, Scuderi A, Letsou A, Virshup DM. B56-associated protein phosphatase 2A is required for survival and protects from apoptosis in Drosophila melanogaster. Mol Cell Biol 2002;22:3674–84.[Abstract/Free Full Text]
  44. Dagda RK, Zaucha JA, Wadzinski BE, Strack S. A developmentally regulated, neuron-specific splice variant of the variable subunit Bß targets protein phosphatase 2A to mitochondria and modulates apoptosis. J Biol Chem 2003;278:24976–85.[Abstract/Free Full Text]
  45. Kong M, Fox CJ, Mu J, et al. The PP2A-associated protein {alpha}4 is an essential inhibitor of apoptosis. Science 2004;306:695–8.[Abstract/Free Full Text]
  46. Lu KH, Patterson AP, Wang L, et al. Selection of potential markers for epithelial ovarian cancer with gene expression arrays and recursive descent partition analysis. Clin Cancer Res 2004;10:3291–300.[Abstract/Free Full Text]
  47. Jazaeri AA, Awtrey CS, Chandramouli GV, et al. Gene expression profiles associated with response to chemotherapy in epithelial ovarian cancers. Clin Cancer Res 2005;11:6300–10.[Abstract/Free Full Text]
  48. Spentzos D, Levine DA, Kolia S, et al. Unique gene expression profile based on pathologic response in epithelial ovarian cancer. J Clin Oncol 2005;23:7911–8.[Abstract/Free Full Text]
  49. Maglott D, Ostell J, Pruitt KD, Tatusova T. Entrez Gene: gene-centered information at NCBI. Nucleic Acids Res 2005;33:D54–8.[Abstract/Free Full Text]




This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Olivero, M.
Right arrow Articles by Di Renzo, M. F.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Olivero, M.
Right arrow Articles by Di Renzo, M. F.


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