Regulatory subunits of PKA define an axis of cellular proliferation/differentiation in ovarian cancer cells
BMC Medical Genomics volume 1, Article number: 43 (2008)
The regulatory subunit of cAMP-dependent protein kinase (PKA) exists in two isoforms, RI and RII, which distinguish the PKA isozymes, type I (PKA-I) and type II (PKA-II). Evidence obtained from a variety of different experimental approaches has shown that the relative levels of type I and type II PKA in cells can play a major role in determining the balance between cell growth and differentiation. In order to characterize the effect of PKA type I and type II regulatory subunits on gene transcription at a global level, the PKA regulatory subunit genes for RIα and RIIβ were stably transfected into cells of the ovarian cancer cell line (OVCAR8).
RIα transfected cells exhibit hyper-proliferative growth and RIIβ transfected cells revert to a relatively quiescent state. Profiling by microarray revealed equally profound changes in gene expression between RIα, RIIβ, and parental OVCAR cells. Genes specifically up-regulated in RIα cells were highly enriched for pathways involved in cell growth while genes up-regulated in RIIβ cells were enriched for pathways involved in differentiation. A large group of genes (~3600) was regulated along an axis of proliferation/differentiation between RIα, parental, and RIIβ cells. RIα/wt and RIIβ/wt gene regulation was shown by two separate and distinct gene set analytical methods to be strongly cross-correlated with a generic model of cellular differentiation.
Overexpression of PKA regulatory subunits in an ovarian cancer cell line dramatically influences the cell phenotype. The proliferation phenotype is strongly correlated with recently identified clinical biomarkers predictive of poor prognosis in ovarian cancer suggesting a possible pivotal role for PKA regulation in disease progression.
The critical role of cAMP acting as a second messenger and exerting control over the regulation of cell growth and differentiation in a wide variety of cell types has been well established [1–3]. Experimental evidence has shown that the selective modulation of two isoforms of cAMP-dependent protein kinase (PKA-I and PKA-II) act as positive and negative intracellular regulators , respectively, of cell growth. PKA-I is only transiently overexpressed in normal cells in response to the physiologic stimuli of cell proliferation while, in contrast, it is constitutively overexpressed in cancer cells and this over-expression is associated with poor prognosis in many different human cancers. The disruption of the normal balance between PKA isozymes is highly associated with tumorigenesis and tumor growth.
Global gene expression profiling by microarray has demonstrated that antisense suppression of RIα in PC3M prostate and LS-174T colon carcinoma cells, exogenously treated with RIα antisense oligonucleotides, simultaneously up-regulates RIIβ and down-regulates a wide range of genes involved in cell proliferation and transformation . Conversely, the vector-mediated overexpression of RIIβ exhibits induction of differentiation genes along with the suppression of cell proliferation and can lead to a reversion of tumor phenotype. Thus, switching of PKA isozymes can cause tumor cells to undergo a phenotypic reversion of malignancy.
To investigate the molecular mechanisms of this phenomenon, ovarian carcinoma cells (OVCAR-8) were transfected with human genes encoding PKA RIα or RIIβ subunits. Ovarian cancer cells were chosen as a useful model system since cAMP-signaling has already been shown to be of vital importance for the normal functioning of the ovary . RIIβ is hormonally up-regulated during follicular development with expression of RIIβ reaching a peak level at a highly differentiated state of follicle development in response to a luteinizing hormone (LH) surge and elevated intracellular concentrations of cAMP. The importance of the sequence of these events led us to expect a high sensitivity of ovarian cells to any modulations of the cAMP- dependent pathway, despite possible deregulation of this system in an established ovarian cancer cell line. Moreover, previous studies had showed the possibility of ovarian cancer cell growth inhibition with the use of RIα antisense oligonucleotides .
OVCAR-8, human ovarian cancer cells, were obtained from DCT-Tumor Repository (NCI – Frederick Cancer Research Center). Tissue culture reagents were purchased from (Invitrogen, Inc., Carlsbad, CA). Monoclonal antibody for RIα, RIIα, and RIIβ were purchased from BD Biosciences Pharmingen (San Diego, CA). Polyclonal antibodies for RAB25 were kindly provided by Dr. K. Cheng from MD Anderson Cancer Center, University of Taxes. All other antibodies were purchased from Santa Cruz Biotechnology, Inc. (Santa Cruz, CA). DOTAP was purchased from Roche Applied Science (Indianapolis, IN). PKA inhibitor H89 and Protease Inhibitor Cocktail I were obtained from EMD Biosciences (Darmstadt, Germany).
Cell culture and treatment
OVCAR-8 cells were grown in RPMI Medium 1640 supplemented with 10% heat-inactivated fetal bovine serum, MEM non-essential amino acids, and antibiotic-antimycotic, in a humidified incubator (95% air and 5% CO2) at 37°C. The stable transfectants, containing retroviral vectors OT1521 or OT1529 with the internal inducible mouse metallothionine-1 (MT-1) promoter, and genes encoding PKA subunits RIα and RIIβ were obtained as previously described . For maximal induction of PKA genes without cytotoxicity, cells were treated with 60 uM ZnSO4 for 6 days prior to the start of the experiment. To exclude other effects of Zn++, parental cells were treated with ZnSO4 in the same manner as transfected cells.
Cell morphology studies
Whole-cell morphology was determined as previously described . Briefly, cells were washed with PBS, fixed with 70% methanol and stained with Giemsa (BioRad, Hercules, CA) according to the manufacturer's instructions. After staining, the cells were viewed under the microscope.
Western blot analysis of proteins in parental and transfected OVCAR-8 cells was performed as described earlier . Briefly, cells were lysed by homogenization in 20 mM Tris-HCl, pH 7.5, 100 mM NaCl, 5 mM MgCl2, 1% NP40, 0.5% sodium deoxycholate, Protease Inhibitor Cocktail I with subsequent centrifugation at 10 000 rpm, 10 min, 4°C. Equal amounts of protein lysate were subjected to SDS-PAGE, transferred to nitrocellulose membranes, and probed with antibodies as indicated in figure legends. Complexes were visualized with the appropriate horseradish peroxidase-conjugated secondary antibody and developed by enhanced chemiluminescence procedure (Santa Cruz Biotechnology, Santa Cruz, CA).
RNA purification, northern blot analysis, microarray probe labeling
Total RNA was extracted using the Trizol Reagent method (Invitrogen, Carlsbad, California 92008, cat. no. 15596-026). Northern blot analysis was performed as described earlier . Additional purification was performed on RNAeasy columns (Qiagen, Valencia, CA 913555, cat. no. 74104). The quality of total RNA samples was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA).
RNA samples were labeled according to the manufacturers recommended protocols. In short, 0.5 μg of total RNA from each sample was labeled by using the Illumina RNA Amplification Kit (Ambion, Austin, TX 78744-1832, cat. no. I1755). Single stranded RNA (cRNA) was generated and labeled by incorporating biotin-16-UTP (Roche Diagnostics GmbH, Mannheim, Germany, cat. no. 11388908910). 0.85 ugs of biotin-labeled cRNA was hybridized (16 hours) to Illumina's Sentrix HumanRef-8 Expression BeadChips (Illumina, San Diego, CA 92121-1975, cat.no. 11201828). The hybridized biotinylated cRNA was detected with streptavidin-Cy3 and quantitated using Illumina's BeadStation 500GX Genetic Analysis Systems scanner.
For Affymetrix, 5 μg of total RNA from each sample was converted into double-stranded cDNA by using the SuperScript Choice system (Invitrogen) with an oligo-dT primer containing the T7 RNA polymerase promoter. The double-stranded cDNA was used in an in vitro transcription reaction using the BioArray RNA transcript labeling kit (Enzo Life Sciences, Inc., Farmingdale, NY 11735). 20 ugs of biotin-labeled cRNA was fragmented and hybridized (16 hours) to Affymetrix GeneChip® Human Genome U133 Plus 2.0 Arrays. The hybridized biotinylated cRNA was detected with phycoerythrin-streptavidin and quantitated by scanning (Affymetrix GeneChip Scanner).
Quantitative RT-PCR (QRT-PCR) analysis
Reverse transcription was performed using total RNA isolated from tissue or cells and processed with Applied Biosystems (Foster City, CA) High-Capacity cDNA Archive kit first-strand synthesis system for RT-PCR according to the manufacturer's protocol. QRT-PCR was performed using the TaqMan assay system from Applied Biosystems. All PCR amplifications were carried out in duplicate on an ABI Prism® 7300 Sequence Detection System, using a fluorogenic 5' nuclease assay (TaqMan® probes). Probes and primers were designed and synthesized by Applied Biosystems. Relative gene expressions are calculated by using the 2-ΔΔCt method, in which Ct indicates cycle threshold, the fractional cycle number where the fluorescent signal reaches detection threshold . The normalized ΔCt value of each sample is calculated using a total of 3 endogenous control genes (gapdh, actb, and pgk1). Fold change values are presented as average fold change = 2-(averageΔΔCt) for genes in treated relative to control samples. Error bars represent the SEM for multiple biological replicates.
Array data analysis
Preliminary analysis of Illumina data was performed using Illumina BeadStudio software which returns the trimmed mean average intensity for each single gene probe type (non-normalized). Any gene consistently below a background threshold level of D = .98 for all samples was eliminated from further analysis. This background filter resulted in the removal of approximately 55% of all the genes on the Illumina array. Z-transformation for normalization was performed on the remaining 11,048 genes for each Illumina sample/array on a stand-alone basis  [see Additional File 1].
Affymetrix data in the form of CEL files was normalized by the RMA method [12, 13] using the web-based Bioconductor software [12, 14]. Probes close to or at background were eliminated if any probe stayed consistently in the lower (average) quartile for all assayed samples.
Significant changes in gene expression were calculated for either RIA or RIIB versus control by Z test [15, 16]. In addition, the mean difference of all calculated changes in gene expression was expressed in units of standard deviation from the average change of all genes for that comparison and referred to as a Z ratio. Z ratios are a direct measure of the likelihood that an observed change is an outlier in an otherwise normal distribution and their use allows for the detection of large changes with high variance which might otherwise be missed by conventional significance testing. Significant genelists were calculated for all categories (RIA or RIIB, up or down, separately for both Affymetrix and Illumina data) by selecting genes which satisfied significance threshold criteria of Z test p- values less than or equal to 0.001 (10-3), a false discovery rate less than or equal to 0.1 , and a fold change ± 2 or greater, or a Z ratio value greater than ± 3.0 [see Additional File 2]. Proliferation/differentiation genes were computed by the simple criteria that their average gene expression was either RIα > Control > RIIβ (proliferation) or RIα < Control < RIIβ (differentiation) [see Additional File 3].
In order to combine data from both Affymetrix and the Illumina gene expression arrays, we first mapped the gene annotations for each platform to the Human Genome Organization (HUGO) human gene symbols and names. Only genes with a positive identification to a HUGO gene were retained and this list was trimmed further to remove genes consistently at or below background as measured for either platform (as described above). Finally, duplicate gene probes were averaged leaving a core group of 7,170 genes common to both platforms, identified and used in subsequent analyses.
The myotube formation dataset was derived from a total of 48 samples corresponding to a time course of nine days of cells in culture following serum withdrawal (John Hall, personal communication). A minimum of 5 replicates at each time point were tested. In addition, cells were tested with and without gamma interferon as an internal control. A total of 8,446 genes were identified as being present and which had a clear functional annotation and were retained for further analysis. Significant changes in gene expression were calculated by the criteria outlined above by either comparing the average gene expression of the entire time course (day 1–day 9) or the interferon treated cells to the average gene expression at baseline.
Hierarchical clustering was performed using the Cluster and TreeView software programs, developed at Stanford University . The clustering algorithm was set to complete linkage clustering using the uncentered Pearson correlation.
GSMA, GSEA, and functional annotation
Gene Set Matrix Analysis (GSMA)  was performed using the median differences for differentially expressed genes tested against genesets derived from a variety of sources (for example, the Pathway genelists were originally obtained from the Gene Set Enrichment Analysis (GSEA) website maintained by the Broad Institute @ http://www.broad.mit.edu/gsea.
PAGE calculations  were automatically derived using a custom script for the GSMA algorithm implemented in JMP (v6.0) the statistical analysis software from the SAS Institute (Cary, NC), according to the formula:
where Sm is the median of Z ratio values of genes for a given gene set and the size of the given gene set is m. The median of total Z ratio values (μ) and standard deviation of total Z ratio values (σ) of a given microarray data set were calculated for all genes between two experimental groups.
A variation of this procedure is the substitution of different gene expression change measures such as fold changes, logratios, or p-values, so, for example, both z ratios and the log of fold changes were independently tested in GSMA for both Affymetrix and Illumina measurements [see Additional File 4]. These measurements were shown to give highly concordant results at the gene set level.
Z ratios were used as the gene expression change metric to test the OVCAR significantly regulated genelists by GSMA against the myoblast-myotube dataset. The top 250 genes (by z ratio) of each significant genelist for either RIα or RIIβ regulated genes were used for both GSMA and GSEA analysis.
Overexpression of PKA subunits
To study the role of PKA isozymes in the regulation of cell proliferation in vitro, OVCAR-8 human ovarian cancer cells were transduced with retroviral vectors containing human PKA subunit genes, including the coding sequences for RIα, RIIα, RIIβ, and the catalytic subunit Cα. As shown in Figure 1A, a marked increase was seen for the mRNA levels of PKA subunits in transfectants, while in cells receiving the control vector, the subunit mRNA levels remained unchanged compared with the parental cells. Similarly, Western blot analysis (as shown in Figure 1B) confirmed the specific overexpression of PKA subunit proteins as compared with the parental cells.
We next examined the growth properties of the PKA subunit gene transfectants in vitro and determined that in monolayer culture, as expected, RIIβ transfectants showed a dramatic growth inhibition while, in contrast, RIα transfectants grew at a much faster rate than parental control cells (Figure 1C–F). These results were not surprising in light of the fact that it had been previously demonstrated that blockade of PKA-I protein with RIα antisense oligonucleotides causes arrest of tumor cell growth, induces apoptosis, inhibits tyrosine kinase signaling, and blocks changes in cell morphology [7, 23–27] while, conversely, preferential expression of PKA-II is found in normal non-proliferating tissues as well as in growth-arrested cells [28, 29].
RIα and RIIβ PKA subunit transfected cells were chosen for further characterization based upon their distinct and diametrically opposed phenotypes relative to each other and to their parental (OVCAR) cell line. The cellular morphology of RIα, parental, and RIIβ cells appeared to define an axis which in the parental to RIα direction leads to uncontrolled cell growth and proliferation and in the parental to RIIβ direction leads to apparent cellular differentiation and quiescence. We set out to determine whether this axis could be defined at the molecular level, on the basis of changes in global gene expression.
Microarray analysis of subunit gene expression
Three biological replicates for each of the two subunit cell lines (RIα and RIIβ) and the parental OVCAR cells were used for microarray measurement analysis using several different platforms including Illumina, Agilent, and Affymetrix . The overall gene expression results were shown to be highly correlated between platforms as illustrated, for example, by the similar results in calculations of selected regulated pathways obtained using either Affymetrix or Illumina (Figure 2C, Figure 3). For the purposes of this report, only the microarray measurements generated using Illumina BeadArrays will be used for all data analysis and results unless otherwise specifically indicated. A total of 11,048 genes were called present out of 24,350 genes tested on the arrays and were used as the basis group for further analysis. All present genes are displayed in the heat map in Figure 2B constructed by processing the data using unsupervised hierarchical clustering . In this view, only the gene expression values have been re-sorted and the sample order has been held constant in order to visibly demonstrate both the reproducibility of replicates (r2 > 0.99)  within samples as well as the dramatic changes in gene expression by sample type. A subset of the results from a preliminary pathway analysis of all the data by Gene Set Matrix Analysis (GSMA)  focusing on the global regulation of either RIα or RIIβ versus parental is displayed in Figure 2C. Multiple pathways are differentially regulated between RIα and RIIβ cell types including genes related to apoptosis, energy metabolism (mitochondrial pathway), mRNA processing and splicing, and cell cycle. Interestingly, regulated genes of the RIα and RIIβ cell types versus parental were shown to be concordant with results from a series of experiments in which fifty fibroblast cultures derived from ten anatomic sites were cultured in 10% fetal bovine serum and a stereotyped gene expression program identifying genes that were reproducibly induced or repressed in fibroblasts in response to serum was identified . These genes (Fig. 2C, CHANG_SERUM_RESPONSE_UP or DN) were enriched in RIα and RIIβ cells for up and down-regulation matching serum stimulated (RIα) or serum depressed (RIIβ) gene expression. Furthermore, Chang et al. experimentally derived a subset of serum activated genes which eliminated the contributions of genes directly related to cell proliferation (i.e. directly related to cell cycle progression). Testing of RIα and RIIβ regulated genes against this core group of genes (SERUM_FIBROBLAST_CORE_UP or DN) gave similar results as before, demonstrating that the changes in gene expression induced by these two PKA subunits went well beyond just simple alterations of cell cycle. The genes induced in the fibroblast serum-response program are expressed in tumors, by the tumor cells themselves, and by tumor-associated fibroblasts. They are evident at an early clinical stage and predict increased risk of metastasis and death in breast, lung, and gastric carcinomas . Further concordance of RIα and RIIβ regulated genes was demonstrated with other empirically derived gene expression signatures related to cancer and cancer progression, including, for example, the down-regulation of genes in RIIβ cells related to DNA replication (E2F1_DNA_UP) , as well as to genes down-regulated following successful treatment of sarcoma cell lines with Ecteinascidin 743 (ET743_SARCOMA_DN) , and to p53-dependent down-regulation of gene expression in the ovarian cancer cell line 2774 following the ectopic expression of p21(CDKN1A) (P21_P53_ANY or MIDDLE_DN) . Conversely, the up-regulation of genes in RIα cells was shown to be concordant with genes up-regulated in mouse skin by the phorbol ester carcinogen, TPA (TPA_SKIN_UP) , with genes down regulated in human fibroblasts due to old-age and progeria (OLDAGE_DN) , with genes associated with the overexpression of human telomerase in human mammary epithelial cells (HMECs) leading to a diminished requirement for exogenous growth factors (SMITH_HTERT_UP) , and to 30 genes whose up-regulation most clearly defined the progression from MM1 to MM4 subgroups in patients with multiple myeloma, correlating to clinical parameters of poor prognosis for this cohort (ZHAN_MULTIPLE_MYELOMA_SUBCLASSES_DIFF) . Figure 3 illustrates a breakout of one of these correlating gene expression signatures on a gene by gene basis (VERNELL_PRB_CLSTR1) which catalogues a group of genes which are up-regulated by E2F and down-regulated by pRB and p16 in human osteosarcoma cells. GSMA gene set analysis indicated that as a group, these genes were down regulated in RIIβ cells and, as shown in Fig. 3, include the specific down-regulation of genes involved in DNA replication and repair (PCNA, FANCA), and control of the cell cycle, and cell division (CDCA5).
Identification and validation of key regulatory genes
Several candidate genes identified in previous studies of ovarian cancer progression were identified in the PKA OVCAR model system by microarray and validated by RT-PCR (Figure 4). These genes include RAB25, a member of the RAS oncogene family, recently implicated in the regulation of cell proliferation and apoptosis in ovarian cancer cells  and with reports that tumor cells overexpressing the RAB25 protein were more aggressive and associated with a poorer clinical outcome . A high level of RAB25 protein in patients with either breast or ovarian cancer was associated with an almost 50% reduction in five year survival rates. RAB25 mRNA levels were elevated between 10- (microarray) to 100-fold (RT-PCR), in RIα transfected OVCAR cells relative to either RIIβ transfected or parental OVCAR cells.
The loss of a tumor suppressor gene, tumor suppressor candidate 3 (TUSC3), was seen in RIα cells compared to either RIIβ transfected or to parental OVCAR cells. Loss of heterozygosity on chromosomal band 8p22 and decreased gene expression of TUSC3 has been associated with an increase of metastatic potential in prostate, colorectal, and ovarian cancer [41–43]. In particular, TUSC3 showed significantly lower expression in grade 3 primary ovarian carcinoma tumors compared with tumors of lower grade or compared with normal controls . Taken together, the dramatic up-regulation of RAB25 and down-regulation of TUSC3 in RIα cells appears to faithfully mimic the activity profiles of proven markers of ovarian cancer clinical progression.
Several other interesting genes displayed direct counter-regulation between RIα and RIIβ cells, including the human baculoviral IAP repeat-containing 3 (BIRC3) gene mRNA. The BIRC3 gene codes for a protein whose function includes the antagonism of the activation of apoptosis-promoting ICE-like proteases. Elevated levels of BIRC3 have been shown by others to promote tumor cell survival [44, 45], and the down-regulation of BIRC3 in RIIβ cells may indicate the release of a block to apoptosis which may account for, in part, for the decreased cell number in this cell type. Counter-regulation in the opposite direction was displayed by the peroxisome proliferator-activated receptor gamma (PPARG) gene, which was highly up-regulated in RIIβ cells and mildly down-regulated in RIα cells. PPARG is a gene which has been shown to play a pivotal role in the processes of cellular differentiation, adipogenesis, and several reports connect PPARG status with neoplastic processes suggesting that PPARG may act as a tumor suppressor for some tissues and in some cellular contexts [46, 47].
PKA regulated genes were tested directly against an unrelated in vitro model system of cellular proliferation/differentiation for confirmation that these data were representative of a generic cellular response. A gene expression profile for a time course of human myotube formation (Figure 5A) was used for comparison and will be described in detail elsewhere. Examination of the myotube gene expression profile shows that the major programmatic changes start almost immediately following serum withdrawal and remain relatively constant throughout the time course.
Genelists were derived from the OVCAR data using the top 250 genes (by size of fold change) significantly regulated for either RIα or RIIβ cells versus parental cells as previously described  and these genelists were tested against all calculated changes in gene expression between day 2–9 and day 0–1 from the myotube dataset for concordance by Gene Set Matrix Analysis (GSMA) analysis . The results as shown in Figure 5B, demonstrated a highly significant correlation for RIα up-regulated genes ((p < 2E-4) and RIIβ down-regulated genes (p < 1E-11) and the gene expression differences which distinguish myoblasts from myotube cells. Genes up-regulated in rapidly dividing myocytes correlate well with genes up-regulated by RIα or down-regulated by RIIβ in OVCAR cells. These apparently unrelated, and certainly biologically distinct, cell types are connected by a common theme at the level of gene expression.
As a cross check to the GSMA calculations, we decided to test the same dataset/genelist combination using another increasingly popular gene set analysis method, Gene Set Enrichment Analysis (GSEA) [48–50]. In this instance, genelists derived from significantly regulated genes in the myotube/myoblast dataset were used as positive controls and OVCAR PKA subunit significant genelists derived from three different commercial platforms  were tested simultaneously against the entire myoblast/myotube dataset. GSEA, without supervision, grouped all the submitted genelists in the predicted direction (Table 1), i.e. RIIβ up-regulated and RIα down-regulated genelists were associated with myotube formation (MYO_UP), while RIα up-regulated and RIIβ down-regulated genelists were associated with myoblasts (MYO_DOWN). In GSEA, the difference data is rank ordered from the largest positive change to the largest negative change (the grey profile at the bottom of each GSEA graph – Figure 6). Genes from the submitted genelists are mapped to this distribution (blue bar code) and enrichment for any particular gene list at either end of the distribution is statistically calculated. The Illumina derived OVCAR gene lists, in particular, performed consistently well, with three of the four lists surviving a stringent permutation estimate of significance (NOM p-val) followed by a multiple comparison correction (FDR q-val < 0.1). The functional annotation of the core group of genes (highlighted in Fig. 6 by red circles) corresponding to the overlap between RIα and myoblast up-regulated genes included genes involved in energy metabolism, mRNA processing, and ribosome formation. Functional annotation of the core group of genes corresponding to the overlap between RIβ and myotube up-regulated genes included genes involved in the formation of extracellular matrix, focal adhesion and cell-cell communication. In the ovarian cancer cell model we have used, a shift in the relative amounts of the RIα and the RIIβ regulatory subunits of PKA is sufficient to trigger a programmatic shift that both morphologically and at the level of gene expression looks very much like a generic program of cellular proliferation and differentiation.
We chose to characterize RIα and RIIβ transfected ovarian cancer cells where the relative levels of Type I and Type II PKA isozymes have been deliberately shifted with obvious and profound consequences (Fig. 1C–E). It appeared likely that these effects would also be paralleled at the level of gene expression and they were indeed reproducible and profound (Fig. 2B). In order to directly investigate the apparent antagonistic effects of RIα and RIIβ in ovarian cancer cells at the gene expression level, a large group of genes was identified solely on the basis that their gene expression progressed through the RIα-parental-RIIβ axis and these genes were tested for biological relevancy by looking for patterns of functional enrichment of genes involved in proliferation (RIα > parental > RIIβ), or differentiation (RIα < parental < RIIβ). We performed gene set analysis on these data and showed that the association of RIα and RIIβ with proliferation/differentiation genes was consistent with the resulting patterns of GSMA enrichment scores, a subset of which is featured in Fig. 2C. One particularly striking example of this enrichment was the precise overlap of RIα and RIIβ with up- and down-regulation in the Chang serum response data . The Chang study investigated aspects of cancer invasion and metastasis that appear to overlap with the gene expression pattern of wound healing (modeled by response to serum exposure) since genes induced in the fibroblast serum-response program are often found to be coordinately up-regulated in many human tumors as well . Our data tends to support their hypothesis of a connection since RIα-induced gene expression not only correlates well with a positive response to serum but also simultaneously significantly up-regulates genes, like RAB25 (Fig. 4), independently associated with advanced metastasis in ovarian cell cancer.
Over-expression of RAB25 in both breast and ovarian cancer cells reportedly decreased apoptosis and increased proliferation of these cells in culture and increased their aggressiveness in vivo . RAB25 was among a group of genes whose overexpression distinguished ovarian/primary peritoneal serous carcinoma (OC/PPC) from diffuse peritoneal malignant mesothelioma (DMPM), two highly aggressive tumor types which are closely related, both morphologically and histogenetically . All evidence to date suggests that increased RAB25 gene expression (both by increases in DNA copy number as well as accompanying increases in mRNA expression levels ) is specifically associated with ovarian cancer progression and our results demonstrate that elevated levels of the RIα PKA regulatory subunit can lead directly to elevated levels of RAB25 in ovarian cancer cells. Similarly, the TUSC gene encodes a protein with oligosaccahryl transferase activity which was originally identified as a homozygous deletion in metastatic prostate cancer . TUSC3 has since been shown to be hypermethylated in acute lymphoblastoid leukemia (ALL) cells , as well as in cervical intraepitelial neoplasia , and to be specifically associated with loss of heterozygosity (LOH) on chromosomal band 8p22, a common event in several epithelial tumors including ovarian carcinoma . Significantly lower expression of TUSC3 was correlated with an increase in clinical grade severity in a study of 58 primary ovarian carcinoma tissues . The overexpression of RIα PKA regulatory subunit leads directly to significantly decreased levels of TUSC3 in the OVCAR PKA subunit model system.
We have accumulated evidence from cell growth and morphology, from pathway involvement in global gene expression patterns, and from key regulatory genes which act as biomarkers for clinical ovarian cancer progression, that the transfection of PKA regulatory subunits in an ovarian cancer cell line directly controls cell fate along an axis which we have characterized as proliferation/differentiation (RIα-parental-RIIβ). We tested this characterization directly by correlating significant changes in gene expression derived from the OVCAR model system against the entire dataset of gene expression changes generated from an in vitro model system of human myotube formation using both a parametric (GSMA) and a non-parametric (GSEA) approach. Both analytical techniques returned results, without supervision, which indicated a significant overlap between the two model systems.
While it is always hazardous to extrapolate without hesitation from immortalized cell lines in tissue culture into the behavior of cells in an intact environment under various physiological conditions and stresses, cell lines, however, are still widely used as in vitro models in cancer research because they are relatively easy to handle, and are a renewable resource that can be grown in almost infinite quantities . In general, they exhibit a relatively high degree of homogeneity and preserve in vitro, the genetic aberrations unique to their parent histology from which they were derived despite the fact that they are prone to a measurable level of genotypic and phenotypic drift during their continual culture . Despite these caveats we have demonstrated that an established ovarian epithelial cancel line, OVCAR8, can be manipulated by the selective overexpression of PKA regulatory subunits into states, which at the level of gene expression, are clearly associated with either progression or differentiation. We show that these states are both general (canonical pathway regulation) as well as specifically correlated with patterns of gene expression related to cancer progression as reported from other publicly available microarray studies. We provide evidence on a gene-by-gene basis that biomarkers specifically associated with poor prognosis in ovarian cancer (elevated levels of RAB25, decreased levels of TUSC3) are dramatically elevated by RIα over expression in a cell line which may well be acting, in this case, as a stage-specific surrogate to a full blown ovarian cancer cell malignancy. Finally, we provide evidence that this stage-specific transition can not only be triggered solely by the differential expression of PKA regulatory subunits but also correlates well with an, otherwise, unrelated model system of proliferation/differentiation generated by the mapping of gene expression changes between myoblast and myotube cells.
The overexpression of PKA regulatory subunits in an ovarian cancer cell model has been used to show that cells enter either a hyperproliferative (RIα) or a quasi-quiescent (RIIβ) phenotypic state which, at the level of gene expression, mimics changes in gene expression associated with good or poor prognosis in ovarian cancer clinical outcomes. These patterns of gene expression are also clearly related to a more generalized pattern of gene expression in other model systems of cellular proliferation and differentiation. Taken altogether, these observations support the conclusion, already suggested by a great deal of previous work, that whether causative or in part as a reaction to other cellular events, the state and condition of PKA remains a fundamental determinant of cell fate.
Gene Set Enrichment Analysis
Gene Set Matrix Analysis
Parametric Analysis of Gene Enrichment.
Cho-Chung YS: Hypothesis. Cyclic AMP and its receptor protein in tumor growth regulation in vivo. J Cyclic Nucleotide Res. 1980, 6 (3): 163-177.
Pastan IH, Johnson GS, Anderson WB: Role of cyclic nucleotides in growth control. Annu Rev Biochem. 1975, 44: 491-522. 10.1146/annurev.bi.44.070175.002423.
Puck TT: Genetic regulation of growth control: role of cyclic AMP and cell cytoskeleton. Somat Cell Mol Genet. 1987, 13 (4): 451-457. 10.1007/BF01534947.
Cho-Chung YS: Role of cyclic AMP receptor proteins in growth, differentiation, and suppression of malignancy: new approaches to therapy. Cancer Res. 1990, 50 (22): 7093-7100.
Cho YS, Kim MK, Cheadle C, Neary C, Becker KG, Cho-Chung YS: Antisense DNAs as multisite genomic modulators identified by DNA microarray. Proc Natl Acad Sci USA. 2001, 98 (17): 9819-9823. 10.1073/pnas.171314398.
Richards JS: Hormonal control of gene expression in the ovary. Endocr Rev. 1994, 15 (6): 725-751. 10.1210/er.15.6.725.
Alper O, Hacker NF, Cho-Chung YS: Protein kinase A-Ialpha subunit-directed antisense inhibition of ovarian cancer cell growth: crosstalk with tyrosine kinase signaling pathway. Oncogene. 1999, 18 (35): 4999-5004. 10.1038/sj.onc.1202830.
Nesterova M, Yokozaki H, McDuffie E, Cho-Chung YS: Overexpression of RII beta regulatory subunit of protein kinase A in human colon carcinoma cell induces growth arrest and phenotypic changes that are abolished by site-directed mutation of RII beta. Eur J Biochem. 1996, 235 (3): 486-494. 10.1111/j.1432-1033.1996.00486.x.
Neary CL, Nesterova M, Cho YS, Cheadle C, Becker KG, Cho-Chung YS: Protein kinase A isozyme switching: eliciting differential cAMP signaling and tumor reversion. Oncogene. 2004, 23 (54): 8847-8856. 10.1038/sj.onc.1208165.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.
Cheadle C, Vawter MP, Freed WJ, Becker KG: Analysis of microarray data using Z score transformation. J Mol Diagn. 2003, 5 (2): 73-81.
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 (2): 185-193. 10.1093/bioinformatics/19.2.185.
Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31 (4): e15-10.1093/nar/gng015.
Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4 (2): 249-264. 10.1093/biostatistics/4.2.249.
Cheadle C, Cho-Chung YS, Becker KG, Vawter MP: Application of z-score transformation to Affymetrix data. Appl Bioinformatics. 2003, 2 (4): 209-217.
Nadon R, Woody E, Shi P, Rghei N, Hubschle H, Susko E, Ramm P: Statistical inference in array genomics. Microarrays for the Neurosciences. Edited by: Geschwind D, Gregg J. 2002, Cambridge: MIT Press
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. JRSS-B. 1995, 57: 289-300.
Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95 (25): 14863-14868. 10.1073/pnas.95.25.14863.
Cheadle C, Watkins T, Fan J, Williams MA, Georas S, Hall J, Rosen A, Barnes KC: Gene Set Matrix Analysis, an automated method for rapid hypothesis testing of gene expression data. Bioinformatics and Biology Insights. 2007, I: 1-14.
Kim SY, Volsky DJ: PAGE: parametric analysis of gene set enrichment. BMC Bioinformatics. 2005, 6: 144-10.1186/1471-2105-6-144.
Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4 (5): P3-10.1186/gb-2003-4-5-p3.
Hosack DA, Dennis G, Sherman BT, Lane HC, Lempicki RA: Identifying biological themes within lists of genes with EASE. Genome Biol. 2003, 4 (10): R70-10.1186/gb-2003-4-10-r70.
Srivastava RK, Srivastava AR, Seth P, Agrawal S, Cho-Chung YS: Growth arrest and induction of apoptosis in breast cancer cells by antisense depletion of protein kinase A-RI alpha subunit: p53-independent mechanism of action. Mol Cell Biochem. 1999, 195 (1–2): 25-36. 10.1023/A:1006990231186.
Wang H, Cai Q, Zeng X, Yu D, Agrawal S, Zhang R: Antitumor activity and pharmacokinetics of a mixed-backbone antisense oligonucleotide targeted to the RIalpha subunit of protein kinase A after oral administration. Proc Natl Acad Sci USA. 1999, 96 (24): 13989-13994. 10.1073/pnas.96.24.13989.
Nesterova M, Cho-Chung YS: Oligonucleotide sequence-specific inhibition of gene expression, tumor growth inhibition, and modulation of cAMP signaling by an RNA-DNA hybrid antisense targeted to protein kinase A RIalpha subunit. Antisense Nucleic Acid Drug Dev. 2000, 10 (6): 423-433.
Cho-Chung YS, Nesterova M, Becker KG, Srivastava R, Park YG, Lee YN, Cho YS, Kim MK, Neary C, Cheadle C: Dissecting the circuitry of protein kinase A and cAMP signaling in cancer genesis: antisense, microarray, gene overexpression, and transcription factor decoy. Ann N Y Acad Sci. 2002, 968: 22-36.
Cho YS, Kim MK, Tan L, Srivastava R, Agrawal S, Cho-Chung YS: Protein kinase A RIalpha antisense inhibition of PC3M prostate cancer cell growth: Bcl-2 hyperphosphorylation, Bax up-regulation, and Bad-hypophosphorylation. Clin Cancer Res. 2002, 8 (2): 607-614.
Cho-Chung YS, Pepe S, Clair T, Budillon A, Nesterova M: cAMP-dependent protein kinase: role in normal and malignant growth. Crit Rev Oncol Hematol. 1995, 21 (1–3): 33-61. 10.1016/1040-8428(94)00166-9.
Tortora G, Ciardiello F: Protein kinase A as target for novel integrated strategies of cancer therapy. Ann N Y Acad Sci. 2002, 968: 139-147.
Cheadle C, Becker KG, Cho-Chung YS, Nesterova M, Watkins T, Wood W, Prabhu V, Barnes KC: A rapid method for microarray cross platform comparisons using gene expression signatures. Mol Cell Probes. 2007, 21 (1): 35-46. 10.1016/j.mcp.2006.07.004.
Chang HY, Sneddon JB, Alizadeh AA, Sood R, West RB, Montgomery K, Chi JT, Rijn van de M, Botstein D, Brown PO: Gene expression signature of fibroblast serum response predicts human cancer progression: similarities between tumors and wounds. PLoS Biol. 2004, 2 (2): E7-10.1371/journal.pbio.0020007.
Kalma Y, Marash L, Lamed Y, Ginsberg D: Expression analysis using DNA microarrays demonstrates that E2F-1 up-regulates expression of DNA replication genes including replication protein A2. Oncogene. 2001, 20 (11): 1379-1387. 10.1038/sj.onc.1204230.
Martinez N, Sanchez-Beato M, Carnero A, Moneo V, Tercero JC, Fernandez I, Navarrete M, Jimeno J, Piris MA: Transcriptional signature of Ecteinascidin 743 (Yondelis, Trabectedin) in human sarcoma cells explanted from chemo-naive patients. Mol Cancer Ther. 2005, 4 (5): 814-823. 10.1158/1535-7163.MCT-04-0316.
Wu Q, Kirschmeier P, Hockenberry T, Yang TY, Brassard DL, Wang L, McClanahan T, Black S, Rizzi G, Musco ML, Mirza A, Liu S: Transcriptional regulation during p21WAF1/CIP1-induced apoptosis in human ovarian cancer cells. J Biol Chem. 2002, 277 (39): 36329-36337. 10.1074/jbc.M204962200.
Schlingemann J, Hess J, Wrobel G, Breitenbach U, Gebhardt C, Steinlein P, Kramer H, Furstenberger G, Hahn M, Angel P, Lichter P: Profile of gene expression induced by the tumour promotor TPA in murine epithelial cells. Int J Cancer. 2003, 104 (6): 699-708. 10.1002/ijc.11008.
Ly DH, Lockhart DJ, Lerner RA, Schultz PG: Mitotic misregulation and human aging. Science. 2000, 287 (5462): 2486-2492. 10.1126/science.287.5462.2486.
Smith LL, Coller HA, Roberts JM: Telomerase modulates expression of growth-controlling genes and enhances cell proliferation. Nat Cell Biol. 2003, 5 (5): 474-479. 10.1038/ncb985.
Zhan F, Hardin J, Kordsmeier B, Bumm K, Zheng M, Tian E, Sanderson R, Yang Y, Wilson C, Zangari M, Anaissie E, Morris C, Muwalla F, van Rhee F, Fassas A, Crowley J, Tricot G, Barlogie B, Shaughnessy J: Global gene expression profiling of multiple myeloma, monoclonal gammopathy of undetermined significance, and normal bone marrow plasma cells. Blood. 2002, 99 (5): 1745-1757. 10.1182/blood.V99.5.1745.
Fan Y, Xin XY, Chen BL, Ma X: Knockdown of RAB25 expression by RNAi inhibits growth of human epithelial ovarian cancer cells in vitro and in vivo. Pathology. 2006, 38 (6): 561-567. 10.1080/00313020601024037.
Cheng KW, Lahad JP, Kuo WL, Lapuk A, Yamada K, Auersperg N, Liu J, Smith-McCune K, Lu KH, Fishman D, Gray JW, Mills GB: The RAB25 small GTPase determines aggressiveness of ovarian and breast cancers. Nat Med. 2004, 10 (11): 1251-1256. 10.1038/nm1125.
Pils D, Horak P, Gleiss A, Sax C, Fabjani G, Moebus VJ, Zielinski C, Reinthaller A, Zeillinger R, Krainer M: Five genes from chromosomal band 8p22 are significantly down-regulated in ovarian carcinoma: N33 and EFA6R have a potential impact on overall survival. Cancer. 2005, 104 (11): 2417-2429. 10.1002/cncr.21538.
Bova GS, MacGrogan D, Levy A, Pin SS, Bookstein R, Isaacs WB: Physical mapping of chromosome 8p22 markers and their homozygous deletion in a metastatic prostate cancer. Genomics. 1996, 35 (1): 46-54. 10.1006/geno.1996.0321.
Andersen CL, Wiuf C, Kruhoffer M, Korsgaard M, Laurberg S, Orntoft TF: Frequent occurrence of uniparental disomy in colorectal cancer. Carcinogenesis. 2007, 28 (1): 38-48. 10.1093/carcin/bgl086.
Gordon GJ, Appasani K, Parcells JP, Mukhopadhyay NK, Jaklitsch MT, Richards WG, Sugarbaker DJ, Bueno R: Inhibitor of apoptosis protein-1 promotes tumor cell survival in mesothelioma. Carcinogenesis. 2002, 23 (6): 1017-1024. 10.1093/carcin/23.6.1017.
Dai Z, Zhu WG, Morrison CD, Brena RM, Smiraglia DJ, Raval A, Wu YZ, Rush LJ, Ross P, Molina JR, Otterson GA, Plass C: A comprehensive search for DNA amplification in lung cancer identifies inhibitors of apoptosis cIAP1 and cIAP2 as candidate oncogenes. Hum Mol Genet. 2003, 12 (7): 791-801. 10.1093/hmg/ddg083.
Grommes C, Landreth GE, Heneka MT: Antineoplastic effects of peroxisome proliferator-activated receptor gamma agonists. Lancet Oncol. 2004, 5 (7): 419-429. 10.1016/S1470-2045(04)01509-8.
Koeffler HP: Peroxisome proliferator-activated receptor gamma and cancers. Clin Cancer Res. 2003, 9 (1): 1-9.
Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstrale M, Laurila E, Houstis N, Daly MJ, Patterson N, Mesirov JP, Golub TR, Tamayo P, Spiegelman B, Lander ES, Hirschhorn JN, Altshuler D, Groop LC: PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003, 34 (3): 267-273. 10.1038/ng1180.
Sweet-Cordero A, Mukherjee S, Subramanian A, You H, Roix JJ, Ladd-Acosta C, Mesirov J, Golub TR, Jacks T: An oncogenic KRAS2 expression signature identified by cross-species gene-expression analysis. Nat Genet. 2005, 37 (1): 48-55.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.
Cheadle C, Becker KG, Cho-Chung YS, Nesterova M, Watkins T, Wood W, Prabhu V, Barnes KC: A rapid method for microarray cross platform comparisons using gene expression signatures. Mol Cell Probes. 2006
Cheng KW, Lahad JP, Gray JW, Mills GB: Emerging role of RAB GTPases in cancer and human disease. Cancer Res. 2005, 65 (7): 2516-2519. 10.1158/0008-5472.CAN-05-0573.
Davidson B, Zhang Z, Kleinberg L, Li M, Florenes VA, Wang TL, Shih Ie M: Gene expression signatures differentiate ovarian/peritoneal serous carcinoma from diffuse malignant peritoneal mesothelioma. Clin Cancer Res. 2006, 12 (20 Pt 1): 5944-5950. 10.1158/1078-0432.CCR-06-1059.
Cheng KW, Lu Y, Mills GB: Assay of Rab25 function in ovarian and breast cancers. Methods Enzymol. 2005, 403: 202-215. 10.1016/S0076-6879(05)03017-X.
Scholz C, Nimmrich I, Burger M, Becker E, Dorken B, Ludwig WD, Maier S: Distinction of acute lymphoblastic leukemia from acute myeloid leukemia through microarray-based DNA methylation analysis. Ann Hematol. 2005, 84 (4): 236-244. 10.1007/s00277-004-0969-1.
Kekeeva TV, Zhevlova AI, Podisov Iu I, Solov'eva Iu V, Zaletataev DV, Nemtsova MV: [Aberrant methylation of tumor suppressor genes and allelic imbalance in cervical intraepitelial neoplasia]. Mol Biol (Mosk). 2006, 40 (2): 224-230.
Burdall SE, Hanby AM, Lansdown MR, Speirs V: Breast cancer cell lines: friend or foe?. Breast Cancer Res. 2003, 5 (2): 89-95. 10.1186/bcr577.
Greshock J, Nathanson K, Martin AM, Zhang L, Coukos G, Weber BL, Zaks TZ: Cancer cell lines as genetic models of their parent histology: analyses based on array comparative genomic hybridization. Cancer Res. 2007, 67 (8): 3594-3600. 10.1158/0008-5472.CAN-06-3674.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/1/43/prepub
This research was supported (in part) by the Intramural Research Program of the NIH, National Institute on Aging (KGB). KCB was supported in part by the Mary Beryl Patch Turnbull Scholar Program.
The authors declare that they have no competing interests.
MN performed the cell culture, transfections, Western and Northern blots. CC and TW performed the microarray assays. CC performed the QRT-PCR, carried out the statistical analysis, and drafted the manuscript. JCH and AR contributed the myotube data. KCB and YSC-C participated in the design of the study. All authors read and approved the final manuscript.
Chris Cheadle, Maria Nesterova contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Cheadle, C., Nesterova, M., Watkins, T. et al. Regulatory subunits of PKA define an axis of cellular proliferation/differentiation in ovarian cancer cells. BMC Med Genomics 1, 43 (2008). https://doi.org/10.1186/1755-8794-1-43