Skip to main content

Suppression of cancer stemness p21-regulating mRNA and microRNA signatures in recurrent ovarian cancer patient samples

Abstract

Background

Malignant ovarian disease is characterised by high rates of mortality due to high rates of recurrent chemoresistant disease. Anecdotal evidence indicates this may be due to chemoresistant properties of cancer stem cells (CSCs). However, our understanding of the role of CSCs in recurrent ovarian disease remains sparse. In this study we used gene microarrays and meta-analysis of our previously published microRNA (miRNA) data to assess the involvement of cancer stemness signatures in recurrent ovarian disease.

Methods

Microarray analysis was used to characterise early regulation events in an embryonal carcinoma (EC) model of cancer stemness. This was then compared to our previously published microarray data from a study of primary versus recurrent ovarian disease. In parallel, meta-analysis was used to identify cancer stemness miRNA signatures in tumor patient samples.

Results

Microarray analysis demonstrated a 90% difference between gene expression events involved in early regulation of differentiation in murine EC (mEC) and embryonic stem (mES) cells. This contrasts the known parallels between mEC and mES cells in the undifferentiated and well-differentiated states. Genelist comparisons identified a cancer stemness signature set of genes in primary versus recurrent data, a subset of which are known p53-p21 regulators. This signature is present in primary and recurrent or in primary alone but essentially never in recurrent tumors specifically. Meta-analysis of miRNA expression showed a much stronger cancer stemness signature within tumor samples. This miRNA signature again related to p53-p21 regulation and was expressed prominently in recurrent tumors. Our data indicate that the regulation of p53-p21 in ovarian cancer involves, at least partially, a cancer stemness component.

Conclusion

We present a p53-p21 cancer stemness signature model for ovarian cancer. We propose that this may, at least partially, differentially regulate the p53-p21 mechanism in ovarian disease. Targeting CSCs within ovarian cancer represents a potential therapeutic avenue.

Background

Ovarian cancer is the leading gynecological malignancy, affecting more than 200,000 women per annum worldwide [1, 2]. This is largely due to high rates of chemoresistant recurrence associated with the disease. Primary ovarian cancer develops silently, with most patients symptom-free, only presenting at an advanced stage. Treatment of primary disease generally consists of surgical removal of the malignancy in combination with platinum-based treatments. In recent years, chemotherapeutic agent carboplatin has proved successful in eliminating primary malignancy while reducing side effects for the patient [Reviewed in [3]]. Mechanistically, platinum-based drugs bind nucleotides within the DNA backbone, causing cross-linking. In response, cells activate DNA-repair mechanisms that ultimately result in apoptosis. Today, the majority of primary ovarian malignancies are successfully treated, where up to 80% of women will recover [2]. The remaining 20% may be explained by late presentation of the disease by asymptomatic women. Alarmingly, up to 80% of these survivors will develop chemoresistant terminal recurrent disease within two years, which is accepted as the main factor in fatality rates. We have previously used comparative microarray analysis to demonstrate that primary and recurrent disease have substantially different gene and microRNA (miRNA) expression profiles [4, 5], which we continue in this study.

Current treatment of recurrent disease, which is similar to treatment of primary disease, has proved ineffective. Thus, recurrent disease must be fully characterised and novel therapeutic approaches developed. One such approach involves targeting cancer cells with stemness properties. These cancer stem cells (CSCs) have been described in ovarian cancer [Reviewed in [6]] and have several properties with relevance to recurrent ovarian cancer. CSCs are sufficient to regenerate malignancy in vivo via extensive self-renewal and differentiation. Tumor regeneration from CSCs is remarkably efficient, where a single CSC is often sufficient to re-establish disease [7, 8]. CSCs proliferate well in the hypoxic conditions found in the tumor microenvironment [9, 10]. As they differentiate, CSCs quickly develop neo-vasculature to fuel further tumorigenesis. Perhaps the most alarming aspect of CSCs is their uninhibited proliferation in the presence of chemotherapeutic agents. It is broadly accepted that CSCs play a role in most, if not all, primary malignancies. Theoretically, the persistence of a single CSC post-intervention could be sufficient to explain chemoresistant recurrence. However, the role of CSCs in recurrent ovarian disease is poorly understood. Ultimately we must develop methods of targeting specific CSC populations as part of a combined anti-cancer strategy.

Many studies have demonstrated the presence of CSCs in ovarian malignancy [6]. However, establishing ovarian CSC models in culture has proved challenging. In this study we employed an embryonal carcinoma (EC) model of cancer stemness. Originally derived from malignant teratomas that can develop in the ovary, EC cells are the original and best characterised CSC model [11–14]. We have previously shown high relevance between EC cells and ovarian serous carcinoma patient samples at the miRNA level [15]. Pluripotent EC cells can differentiate into cells representing all three germ layers and are considered the malignant equivalent of embryonic stem (ES) cells [11–14]. Nullipotent EC cells can avoid differentiation in vivo to generate poorly-differentiated, highly-malignant tumors [11–14]. Comparison of ES cells with pluripotent and nullipotent EC cells can establish mechanisms required for functional malignant differentiation. The cells are so similar that EC cells are used as an easily cultured model of ES biology, reflecting the difficulty of targeting CSCs without damaging non-malignant stem cell populations [16–18].

In this study we first used gene microarrays to assess upstream regulation of differentiation in murine EC (mEC) and mES cells. Our analysis describes aberrant regulation of differentiation in EC cells. Subsequently, we compared mEC genelists to our previously published primary versus recurrent tumor sample data [5]. We described the presence of a cancer stemness p53-p21 regulatory mechanism in ovarian tumor samples. This mechanism is employed by primary disease and suppressed in recurrent disease. Subsequently, we conducted a meta-analysis of our previously published human EC (hEC) and tumor sample miRNA data [15, 6]. We report that cancer stemness signature miRNAs are more relevant to ovarian cancer than cancer stemness signature genes. We detail substantial recruitment of stemness signature miRNAs by recurrent disease. Thus recurrent tumors suppress and activate stemness signature genes and miRNAs respectively. Our analysis indicates that cancer stemness mechanisms are specifically and differentially regulated in primary and recurrent ovarian malignancy, with obvious implications for treatment.

Methods

Cell Culture

Murine ES (ES-E14TG2a) and EC cells (pluripotent 'SCC-PSA1' and nullipotent 'Nulli-SCC1') were purchased from ATCC, cultured on murine irradiated fibroblasts in DMEM supplemented with 10% foetal bovine serum, 4 mM L-glutamine (Invitrogen) and 100 U/ml of penicillin/streptomycin (Invitrogen Corporation, Carlsbad, CA, USA) and spontaneously differentiated via removal of feeder layer. Human EC cells were retinoic acid-differentiated as previously described [15].

Tumor Samples

Tumor sample data was previously published [5, 6]. Briefly, two cohorts of primary and recurrent samples were assessed. Cohort 1 contained 5 primary and recurrent serous papillary adenocarcinomas (Grade 3). Cohort 2 contained 3 paired ovarian cancers from the same patient but with different histologies: papillary serous, mixed mullerian and clear cell carcinomas.

Microarray Analysis

RNA was isolated using the RNeasy kit (Qiagen, West Sussex, UK) as per manufacturer's protocol. Digoxigenin-UTP labelled cRNA was synthesized via the Chemiluminescent RT-IVT Labelling Kit v2.0 (Life Technologies, Foster City, CA, USA) and hybridized to Mouse Genome Survey arrays (Life Technologies) as per manufacturers' instructions. Data was filtered to a signal/noise ration threshold > 3 in at least one sample using R and further analysed using Spotfire® (Life Technologies). Genelists were generated using cut-offs of 0.05 (p-Value) and ± 2.0 (fold change). Functional relationships were analysed using DAVID [19, 20]. Pathways associations of predicted targets of miRNAs highlighted were generated using DIANA miRPath [21] using cut-offs of ≥ 2 genes per pathway and p-value ≤ 0.05.

qPCR Analysis

  1. 2

    μg total RNA was used to synthesis cDNA using the High Capacity cDNA Archive Kit (Life Technologies) as per manufacturer's instructions. Microarrays were validated using 36 pre-designed TaqMan assays (Life Technologies). Gene expression values were generated using the 2^-ddCt method [22]. microRNA was isolated using the mirVANA kit (Ambion) and miRNA TaqMan qPCR (Life Technologies) [23] analysis carried out as previously described [5, 6, 15]. Data plotted represents the mean value across a minimum of n = 3. Error bars represent standard error of the mean.

Results

Microarray analysis of early mEC and mES differentiation

It is well established that ES and EC cells express similar gene profiles in the undifferentiated and well-differentiated (one week or later) states [16, 11, 17, 18]. In contrast, our understanding of the earlier, upstream regulation of differentiation is sparse. We hypothesized that comparison of early differentiation of mES and mEC cells would identify cancer-specific differences in upstream regulation of stem cell differentiation. Addressing this we used microarray analysis to assay early (three day) differentiation of mES and mEC cells.

Microarray data was validated through qPCR analysis, showing excellent correlation (Figure 1). An overview of the number of differentially expressed genes in pluripotent (SCC-PSA1) and nullipotent (Nulli-SCC) mEC and mES cells is shown in Table 1. At cut-offs of 0.05 (p-Value) and ± 2.0 (fold change) SCC-PSA1 cells alter the expression of 724 genes: 202 upregulated and 522 downregulated at fold change levels between +18 and -18 (Additional File 1). Top ten SCC-PSA1 genes are characterised by receptor activity and growth and differentiation/development roles (Table 2). Noteworthy events include upregulation of apoptosis (Bcl)-related gene Bid3 and downregulation of Cav2 tumor suppressor [24] and metastasis-linked Nupr1 [25]. Functional relationship analysis identified upregulation of developmental pathways and downregulation of transcription regulation processes and Toll-Like Receptor (TLR), Interleukin-2 and cancer pathways (Additional File 1).

Figure 1
figure 1

Validation of mEC and mES microarray data. Microarray data was validated through qPCR (TaqMan) analysis of a group of 36 genes. Data presented represents gene expression change in differentiated cells compared to undifferentiated and shows good correlation for mES (A) and PSA-SCC1 (B) datasets.

Table 1 An overview of the numbers and percentage overlap of differentially expressed genes (D/U) during early differentiation of mES and mEC cells.
Table 2 Top ten genes differentially expressed (D/U) during early differentiation of mES and mEC cells.

Nulli-SCC cells responded to differentiation stimuli through the upregulation of 185 and downregulation of 152 genes at levels from -6.3 to 14.0 fold (Additional File 2). Top ten genes included signal transducers and regulators of development/differentiation and malignancy (Table 2). Notable genes include hypoxia and tumor growth regulator Loxl2 [26] and tumor suppressor Serpini2 [27]. Interestingly Ssa2 is downregulated, a gene that is commonly expressed on the surface of apoptotic cells. Functional analysis identified upregulation of signal transduction regulators and downregulation of growth regulators (Additional File 2).

Upstream differentiation of mES cells is characterized by substantial levels of upregulation: 554 upregulated and 832 downregulated genes at levels of 232 to -68 fold (Additional File 3). Top ten genes are populated with receptors and developmental regulators (Table 2). Tll1 is linked to cardiac development, the first organised system formed during embryogenesis. Notably, a key RNAi gene, Eif2c4, is upregulated during differentiation, perhaps reflective of involvement of the RISC complex [28]. Upregulated mES genes regulate development, signalling and gene expression while downregulated genes regulate morphogenesis, particularly growth factor binding. Stemness-linked pathways such as Wnt-catenin and Hedgehog signalling were upregulated while signalling pathways including TLR and TGF-ß were downregulated (Additional File 3).

Aberrant upstream regulation of differentiation in mEC cells

A comparison of mES and mEC early differentiation genelists is summarised in Table 1 and detailed in additional files 1, 2 and 3. In contrast to documented undifferentiated and well-differentiated comparisons, 90% of the mES genelist differed to the mEC genelist at this earlier time point (Table 1). Similarly, almost 70% of the SCC-PSA1 genelist differed from the mES genelist (Table 1). Functional relationship analysis indicates that quite different mechanisms are activated during early differentiation of mEC and mES cells. This included mES-specific upregulation of p53 signaling pathway genes (Additional File 3). There is very little overlap between Nulli-SCC and the other cell types (Additional Files 1, 2 and 3). Only four genes are upregulated by SCC-PSA1 and downregulated by Nulli-SCC cells, while only two are downregulated by both cell types. The downregulation of symporters, signal-transducing membrane proteins, which are upregulated by pluripotent cells, may indicate a potential counteraction of differentiation. Upstream regulation of differentiation represents a substantial difference between these cell types, supporting our hypothesis. While similar genes maintain the self-renewal state in each cell, different mechanisms are employed to regulate the early events in differentiation.

A SCC-PSA1 p53 mechanism is expressed in primary and maintained in recurrent tumors

We have previously published microarray analysis of primary versus recurrent tumor samples [5, 6]. The study contained two cohorts. Cohort 1 represents a group of matched primary and recurrent tumors while cohort 2 represents primary and recurrent tumor samples from the same patients. In this study, raw microarray data from the primary versus recurrent study was reanalysed in an identical fashion to mES and mEC data described above (Additional File 4). Primary versus recurrent disease and mEC genelists were then compared. Genes altered similarly in mEC and mES data were not considered to be cancer-specific and were removed from this analysis. Comparison of mEC and tumor data identified 16 SCC-PSA1 genes expressed in tumor samples (Figure 2, Table 3). These genes group into those that are A) expressed in both primary and recurrent tumors and B) those expressed in primary but not recurrent tumors (Table 3). Many of these genes have links to stemness and malignancy. Tmprss2 is a transmembrane signalling protein that is upregulated in prostate cancer [29]. Cthrc1 is a Smad2/3 (TGF-β signalling) inhibiting Wnt signalling modulator that is differentially expressed in invasive breast cancer and several solid tumors [30]. Nkx3-1 is a metastatic marker transcription factor expressed in prostate cancer [31]. Pdgfc is a cisplatin-associated growth factor [32]. Col4a5 is linked to several cancers while Plaur is a regulator of tissue reorganisation [33]. Ndufs6 is an oxidative phosphorylation enzyme linked to cervical cancer [34]. Pdzk1 is linked to oestrogen-sensitivity in breast and ovarian cancers [35]. Sdsl is a cancer-specific metabolic enzyme [36]. Only one gene, Gata2, an endodermal differentiation marker [37] was upregulated by mEC cells and expressed higher in primary tumors than in recurrent.

Figure 2
figure 2

Expression of SCC-PSA1 stemness signature mRNAs in recurrent tumors. Sixteen genes were identified as differentially expressed in SCC-PSA1 mEC cells and in tumor data. Data is presented as the percentage change in gene expression in differentiated compared to undifferentiated SCC-PSA1 cells (blue) and in primary tumor samples compared to recurrent (red). Gene Expression values are detailed in table 3 (Microarray data p-Value ≤ 0.05).

Table 3 Percentage gene expression of mEC-specific genes expressed in primary versus recurrent tumor samples (Group A expressed similarly in primary and recurrent samples).

When scrutinised, we noted that several of the genes highlighted above have been defined as p53 regulators in various models, as now described. Dusp26 is a p53-inhibiting phosphatase that negatively regulates proliferation of epithelial cells [38]. Stemness gene Sox4 is a p16 and p53 regulator in cancer cells [39] while Hsf2 is a regulator of p53 stability [40]. Hoxb2 has been linked to p205 regulation of p53 and is a well known regulator of EC differentiation [41]. Collectively, our analysis indicates that both primary and recurrent ovarian tumors express this 'p53-regulating stemness signature'.

A NULLI-SCC p21 mechanism is suppressed by recurrent tumors

Despite the reduced genelist size, 14 Nulli-SCC genes were expressed in A) both primary and recurrent tumors or B) primary tumors only (Figure 3, Table 3). These genes related to apoptosis/cellular proliferation, signaling and regulation. Dgcr8 (Pasha) is a key miRNA biosynthesis gene [42], while Tirap is a regulator of TLR signaling [43]. TNF-family related Ltbr and hypoxia-linked Egln3 are apoptosis regulators [44, 45]. Gpr6 is a development regulator expressed in umbilical cord cells [46]. Ndufab1 is a TGF-β signaling related NADPH enzyme [47]. Slc15a1 is involved in drug absorption in the small intestine and has been linked to several cancers and metastasis [48]. Coupled with this is the recurrent suppression of apoptosis regulators Bnip3 and Stau2 [49, 50]. Notably, two p21 regulators are expressed higher in primary tumors compared to recurrent: Cask mediates the expression of p21 to control cell proliferation [51] while Pak6 is a p21 interacting kinase that is a required for chemoresistance in prostate cells [52]. Collectively, an EC cancer stemness signature expressed in tumor samples is linked to maintained p53 regulation and suppression of p53's main target, p21, in recurrent disease.

Figure 3
figure 3

Expression of Nulli-SCC stemness signature mRNAs in recurrent tumors. Fourteen genes were identified as differentially expressed in Nulli-SCC mEC cells and in tumor data. Data is presented as the percentage change in gene expression in differentiation-stimulated compared to undifferentiated Nulli-SCC cells (blue) and in primary tumor samples compared to recurrent (red). Gene Expression values are detailed in table 3 (Microarray data p-Value ≤ 0.05).

Recruitment of cancer stemness signature miRNAs during recurrence

Having identified gene level overlaps, we next conducted overlap meta-analysis of our previously published miRNA data for primary and recurrent patient samples and human EC (hEC) early (three day) differentiation [5, 6, 15]. The earlier study identified cancer stemness signature miRNAs: those miRNAs involved in the differentiation of hEC cells. Specifically, our previous tumor study highlighted 60 miRNAs (52 up and 8 downregulated) in recurrent disease [6]. Of these, 55 miRNAs (92%) are expressed in hEC cells (Additional File 5). 21 recurrent disease-specific miRNAs are linked to differentiation of pluripotent NTera2 hEC cells (Figure 4). We have previously shown that nullipotent 2102Ep hEC cells express a large number of miRNAs at substantially higher levels than NTera2 cells [15]. Here we report that 26 (43%) recurrent disease-specific miRNAs are expressed at higher levels in 2102Ep cells than in NTera2 (Figure 5). Thus, development of recurrent tumors involves recruitment of cancer stemness signature miRNAs. Specific examples include miR-9, which is the most downregulated miRNA in recurrent tumors and is > 1000% higher expressed in undifferentiated 2102Ep cells compared to NTera2, and miR-206, which is in the top ten miRNAs upregulated by recurrent tumors and downregulated during NTera2 differentiation. Molecular pathway relationships between predicted gene targets of the miRNAs highlighted were identified using DIANAmirPATH (Additional file 5). While little pathway overlap was observed in gene array data, miRNA data showed strong pathway associations. Pathway analysis highlighted alteration of several cancer pathways (miRs-10b, -100, -106, -107, -128 and let-7g) as well as Wnt and TGF-β stemness signaling pathways (mirs-10b, -100, -106, -128a and 137). Finally, we assessed the expression of p53-p21 regulating miRNAs in these datasets. Two miRNAs, miRs-106a and b, are validated targets of p21 [52] that are upregulated in recurrent disease and expressed in hEC cells. Notably, miR-106b expression in 2102Ep cells is double that of NTera2 cells. In contrast, miR-155, the only validated p53-regulating miRNA, is unaltered in recurrent tumors. We note that the p53 signaling pathway was highlighted for let-7g and miRs-106b and -107 in pathway analysis (Additional file 5). In overview, we find that miRNAs linked to 2102Ep malignancy are highly relevant to primary and recurrent tumors.

Figure 4
figure 4

Expression of NTera2 stemness signature miRNAs in recurrent tumors. Twenty one miRNAs were identified as differentially expressed in differentiated NTera2 hEC cells and in tumor data. Data is presented as the percentage change in miRNA expression in differentiated NTera2 hEC cells (blue) and in recurrent tumors compared to primary (red). Values represent the mean of at least n = 3 and error bars the standard error of the mean.

Figure 5
figure 5

Expression of 2102Ep stemness signature miRNAs in recurrent tumors. Twenty six miRNAs were identified as differentially expressed in undifferentiated 2102Ep compared to undifferentiated NTera2 hEC cells and in tumor data. Data is presented as log10 (fold change). miRNAs presented showed altered expression in undifferentiated 2102Ep cells compared to undifferentiated NTera2 (blue) and in recurrent tumors compared to primary (red). Values represent the mean of at least n = 3 and error bars the standard error of the mean.

Discussion

Although CSCs are obvious suspects in the development of recurrent ovarian malignancy, a relationship has yet to be established or described in detail. Anecdotal evidence includes altered regulation of Notch3 in chemoresistant ovarian disease and the clear parallel between epithelial-mesenchymal transition (EMT) and CSC differentiation mechanisms [53, 54]. In this study we conducted microarray and meta-analysis of mRNA and miRNA expression in primary and recurrent tumor samples and an EC model of cancer stemness. Our analysis reiterates that development of primary and recurrent ovarian disease involves quite different mechanisms: thousands of genes are differentially expressed. At the gene level, recurrent tumors appear to repress a cancer stemness signature related to p53-p21 regulation. In parallel, recurrent tumors recruit a population of miRNAs with close links to the development of highly malignant, poorly-differentiated tumors from nullipotent hEC cells.

Different genetic profiles are employed by primary and recurrent ovarian tumors [5, 6]. In this study we demonstrate that malignant stem cell differentiation genes are expressed in either primary tumors or both primary and recurrent tumors but essentially never in recurrent tumors-specifically. Some CSC mechanisms are similarly employed in primary and recurrent tumorigenesis. In addition, an obvious implication of our study is that CSCs that survive chemotherapy to repopulate recurrent disease can do so using different mechanisms than those employed in primary disease. Functional relationship analysis indicated that these stemness signature genes have a particular relevance to cellular proliferation and apoptosis. Several of the genes highlighted are known 'p53-p21 signaling regulators'. Mechanistically this relates to regulation of p53-p21 processes, where p53 regulation is enhanced and p21 regulation no longer required in recurrent tumors. This is supported by increased expression of p21 repressing miRNAs in recurrent tumors and strong predicted targeting of p53 signaling genes by tumor-specific miRNAs. Altered p53-p21 regulation is the primary mechanism through which cancers avoid apoptosis and stimulate cellular proliferation. Predictably, we did not find loss of p53 or p21 in recurrent disease (data not shown). It appears that p53-p21 regulation is required at both stages of ovarian malignancy. In Figure 6 we present a schematic to illustrate the p53-p21 regulators highlighted in out study. We propose that these genes and miRNAs regulate p53-p21 signaling, at least partially, in primary and recurrent disease (Figure 6). Indeed, this is likely to be a component of a larger mechanism. This p53-p21 regulating component appears to play a role in primary tumors that is not used during recurrence. We refer to this as a p53-p21 regulating mechanism within the cancer stemness signature (genes altered during differentiation of EC cells but not by ES cells). As a key tumorigenesis component, differential regulation of stemness-linked p53-p21 mechanisms in primary and recurrent disease is an important outcome of this study and will be the subject of ongoing analysis.

Figure 6
figure 6

Suppression of cancer stemness p53-p21 regulation in recurrent tumors. Primary disease is characterised by the expression of p53-regulating stemness signature genes Sox4 and Sdsl, which is continued and enhanced with Dusp26 expression in the recurrency. In contrast, p21 regulating stemness signature genes Pak6 and Cask are expressed in primary disease and suppressed in recurrency. This process is paralleled by recruitment of stemness signature miRNAs by recurrent disease.

It is well established that EC and ES cells are highly similar in the undifferentiated and well-differentiated states [16–18]. This illustrates the significant challenges to the concept of targeting CSCs in a manner that does not harm the non-malignant stem cell pool. In this study we have identified upstream regulation of differentiation as a substantial difference between EC and ES cells, supporting our hypothesis. While downregulated mEC and mES genes displayed similarity, upregulated SCC-PSA1 genes were almost 90% specific to malignancy. This supports a model where normal and malignant stem cells employ similar mechanisms to maintain the self-renewal state. The different phenotypes developing from differentiation, therefore, are related to activation of specific malignant or non-malignant genes. Both cell types alter genes related to similar processes: receptor-mediated signalling of development/differentiation. Thus the differentiation of malignant and non-malignant cells is driven by a divergent group of genes. It is noteworthy that the primary-recurrent genetic switch contained an equally strong Nulli-SCC cell signature, despite the much reduced genelist. Nulli-SCC cells avoid differentiation through maintained levels of gene and miRNA expression to generate highly malignant tumors [11]. While a small number of molecular events take place in these cells response to differentiation, these appear to have a particular relevance to the difference between primary and recurrent disease. Stemness genes are never expressed by recurrent disease only, suggesting a less stem-like profile. These genes have a particular relevance to cellular proliferation and apoptosis, including p53-p21 regulation. Of particular note is the downregulation in Nulli-SCC cells of TLR signaling adapter Tirap, a gene that is constantly expressed in primary and recurrent disease. TLR signaling has received increased attention in both cancer and stemness studies in recent years [55]. In summary, recurrent disease appears to have more correlation with nullipotent cells rather than EC cells with good differential potential.

Recurrent tumor development involves the suppression of twice as many genes as are specifically activated (Cohort 1). This indicates that recurrent malignancy does not require a substantial number of mechanisms employed by primary tumors. Specifically, angiogenesis and development genes are turned off by recurrent disease as malignancy genes are turned on. The upregulation of polycystic ovary-associated gene Fabp4 and ovarian cancer gene Prkcbp1 may be of particular importance. There was little overlap between genes altered in cohort 1 and cohort 2, which altered genes more associated with malignancy and less with differentiation. Functional relationship analysis revealed that recurrent disease no longer requires homeostasis or stimulus response processes while upregulating catalytic activity and protein binding process. In general, recurrent disease behaves more as a developing cancer rather than the chemical stress responses required by primary disease.

Conclusion

CSCs targeting is a potential avenue through which treatment of recurrent, chemoresistant ovarian cancer may be improved. This is complicated by the similarities between cancer and non-cancer stem cells and our poor understanding of recurrent ovarian disease. We have identified the early events of stem cell differentiation as a key area of difference between cancer and non-cancer stem cells. Furthermore, we have highlighted the association of a p53-p21 related cancer stemness signature within ovarian disease. Our data suggests that a stem cell involved in development of recurrent disease employs different mechanisms of tumorigenesis. Our study suggests that it may be possible to target early differentiation events in CSCs without damaging non-cancer stem cells, which would have broad implications for treatments. Our data indicates that such therapies should be independently tailored for primary and recurrent ovarian disease. CSC targeting during treatment of primary disease is likely to have a negative impact on recurrent tumorigenesis. CSC targeting in recurrent disease should be developed with consideration to independent mechanisms. Development of strategies to achieve this will continue in our group.

Abbreviations

CSC:

cancer stem cell

EC:

embryonal carcinoma

h:

human

mEC:

murine embryonal carcinoma (cells)

mES:

murine embryonic stem (cells)

miRNA:

microRNA

P/R:

Primary compared to recurrent

qPCR:

quantitative polymerase chain reaction

U/D:

Undifferentiated compared to differentiated

References

  1. Jemal A, Siegel R, Ward E, Hao Y, Xu J, Murray T, Thun MJ: Cancer statistics, 2008. CA 2008, 58: 71–96.

    PubMed  Google Scholar 

  2. Cannistra SA: Cancer of the ovary. N Engl j med 2004,351(24):2519–29. 10.1056/NEJMra041842

    Article  CAS  PubMed  Google Scholar 

  3. Muggia F: Platinum compounds 30 years after the introduction of cisplatin: implications for the treatment of ovarian cancer. Gyne oncol 2009, 112: 275–81. 10.1016/j.ygyno.2008.09.034

    Article  CAS  Google Scholar 

  4. Laios A, O'Toole SA, Flavin R, Martin C, Ring M, Gleeson N, D'Arcy T, McGuinness EP, Sheils O, Sheppard BL, O' Leary JJ: An integrative model for recurrence in ovarian cancer. Molec cancer 2008, 7: 8. 10.1186/1476-4598-7-8

    Article  Google Scholar 

  5. Laios A, O'Toole S, Flavin R, Martin C, Kelly L, Ring M, Finn SP, Barrett C, Loda M, Gleeson N, D'Arcy T, McGuinness E, Sheils O, Sheppard B, O' Leary J: Potential role of miR-9 and miR-223 in recurrent ovarian cancer. Molec cancer 2008, 7: 35. 10.1186/1476-4598-7-35

    Article  Google Scholar 

  6. Berry NB, Bapat SA: Ovarian cancer plasticity and epigenomics in the acquisition of a stem-like phenotype. J Ov Res 2008, 1: 8. 10.1186/1757-2215-1-8

    Article  Google Scholar 

  7. Kleinsmith LJ, Pierce GB: Multipotentiality of single embryonal carcinoma cells. Cancer Res 1964, 24: 1544–1552.

    CAS  PubMed  Google Scholar 

  8. Al-Hajj M, Wicha MS, Benito-Hernandez A, Morrison SJ, Clarke MF: Prospective identification of tumorigenic breast cancer cells. PNAS, USA 2003, 100: 3983–3988. 10.1073/pnas.0530291100

    Article  CAS  Google Scholar 

  9. Silvan U, Diez-Torre A, Arluzea A, Andrade R, Silio M, Arechaga J: Hypoxia and pluripotency in embryonic and embryonal carcinoma stem cell biology. Differentiation 2009, 78: 159–168. 10.1016/j.diff.2009.06.002

    Article  CAS  PubMed  Google Scholar 

  10. Keith B, Simon SC: Hypoxia-inducible factors, stem cells and cancer. Cell 2007, 129: 465–472. 10.1016/j.cell.2007.04.019

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  11. Andrews PW, Martin MM, Bahrami AR, Damjaov I, Gokhale P, Draper JS: Embryonic stem (ES) cells and embryonal carcinoma (EC) cells: opposite sides of the same coin. Biochemical Society Transactions 2005,33(Pt 6):1526–30.

    CAS  PubMed  Google Scholar 

  12. Andrews PW, Fenderson B, Hakomori S: Human embryonal carcinoma cells and their differentiation in culture. Int j andrology 1987,10(1):95–104. 10.1111/j.1365-2605.1987.tb00170.x

    Article  CAS  Google Scholar 

  13. Andrews PW, Damjanov I, Simon D, Banting GS, Carlin C, Dracopoli NC, Fogh J: Pluripotent embryonal carcinoma clones derived from the human teratocarcinoma cell line Tera-2. Differentiation in vivo and in vitro. Lab Invest 1984,50(2):147–162.

    CAS  PubMed  Google Scholar 

  14. Andrews PW: From teratocarcinomas to embryonic stem cells. Phil Trans R Soc Lond, B 2002, 357: 405–417. 10.1098/rstb.2002.1058

    Article  Google Scholar 

  15. Gallagher MF, Flavin RJ, Elbaruni SA, McInerney JK, Smyth PC, Salley YM, Vencken SF, O'Toole SA, Laios A, Lee MY, Denning K, Li J, Aherne ST, Lao KQ, Martin CM, Sheils OM, O'Leary JJ: Regulation of microRNA biosynthesis and expression in 2102Ep embryonal carcinoma stem cells is mirrored in ovarian serous adenocarcinoma patients. J ov res 2009, 2: 19. 10.1186/1757-2215-2-19

    Article  Google Scholar 

  16. Josephson R, Ording CJ, Liu Y, Shin S, Lakshmipathy U, Toumadje A, Love B, Chesnut JD, Andrews PW, Rao MS, Auerbach JM: Qualification of embryonal carcinoma 2102Ep as a reference line for human embryonic stem cell research. Stem cells 2007, 25: 427–446.

    Article  Google Scholar 

  17. Schwartz CM, Spivak CE, Baker SC, McDaniel TK, Loring JF, Nguyen C, Chrest FJ, Wersto R, Arenas E, Zeng X, Freed WJ, Rao MS: NTera2: a model system to study dopaminergic differentiation of human embryonic stem cells. Stem cells dev 2005, 14: 517–534. 10.1089/scd.2005.14.517

    Article  CAS  PubMed  Google Scholar 

  18. Skotheim RI, Lind GE, Monni O, Nesland JM, Abeler VM, Fosså SD, Duale N, Brunborg G, Kallioniemi O, Andrews PW, Lothe RA: Differentiation of human embryonal carcinomas in vitro and in vivo reveals expression profiles relevant to normal development. Cancer res 2005, 65: 5588–5560. 10.1158/0008-5472.CAN-05-0153

    Article  CAS  PubMed  Google Scholar 

  19. Huang DW, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID Bioinformatics Resources. Nat Protocols 2009,4(1):44–57.

    Article  CAS  Google Scholar 

  20. Dennis G Jr, 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

    Article  PubMed  Google Scholar 

  21. Papadopoulos GL, Alexiou P, Maragkakis M, Reczko M, Hatzigeorgiou AG: DIANA-mirPath: Integrating human and mouse microRNAs in pathways. Bioinfor 2009., 10: (1093/bioinformatics/btp29)

    Google Scholar 

  22. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2 -ddCt method. Methods 2001, 25: 402–8. 10.1006/meth.2001.1262

    Article  CAS  PubMed  Google Scholar 

  23. Lao K, Xu NL, Yeung V, Chen C, Livak KJ, Straus NA: Multi-plexing RT-PCR for the detection of multiple miRNA species in small samples. Bioc Bioph Res Comm 2006, 343: 85–9. 10.1016/j.bbrc.2006.02.106

    Article  CAS  Google Scholar 

  24. Han SY, Druck T, Huebner K: Candidate tumor suppressor genes at FRAG7G are coamplified with MET and do not suppress malignancy in gastric cancer. Genomics 2003,18(2):105–7.

    Article  Google Scholar 

  25. Cano CE, Hamidi T, Sandi MJ, Iovanna JL: Nupr1: the swiss-knife of cancer. J cell physiol 2011,266(6):1439–43.

    Article  Google Scholar 

  26. Bocca C, Bozzo F, Cannito S, Parola M, Migilietta A: Celecoxib inactivates epithelial-mesenchymal transition stimulated by hypoxia and/or epidermal growth factor in colon cancer cells. Molec carcinog 2011. (ePub)

    Google Scholar 

  27. Regg H: The role of serpins in the surveillance of the secretory pathway. Cell molec life sci 2007,64(21):2763–70. 10.1007/s00018-007-7157-0

    Article  Google Scholar 

  28. Yoda M, Kawamata T, Paroo Z, Ye X, Iwasaki S, Liu Q, Tomari Y: ATP-dependent human RISC assembly pathways. Nat struc molec biol 2010,17(1):17–23. 10.1038/nsmb.1733

    Article  CAS  Google Scholar 

  29. Schmidt C: Unraveling the prostate genome. J nat cancer inst 2011,103(9):711–3. 10.1093/jnci/djr167

    Article  PubMed  Google Scholar 

  30. LeClair RJ, Durmus T, Wang Q, Pyagay P, Terzic A, Lindner V: Cthrc1 is a novel inhibitor of transforming growth factor-beta signaling and neointimal lesion formation. Circ res 2007,100(6):826–33. 10.1161/01.RES.0000260806.99307.72

    Article  CAS  PubMed  Google Scholar 

  31. Gurel B, Ali TZ, Montgomery EA, Begum S, Hicks J, Goggins M, Eberhart CG, Clark DP, Bieberich CJ, Epstein JI, De Marzo AM: NKX3.1 is a marker of prostatic origin in metastatic tumors. Am j surg path 2010,34(8):1097–105. 10.1097/PAS.0b013e3181e6cbf3

    Article  PubMed Central  PubMed  Google Scholar 

  32. Yamano Y, Uzawa K, Saito K, Nakashima D, Kasamatsu A, Koike H, Kouzu Y, Shinozuka K, Nakatani K, Negoro K, Fujita S, Tanzawa H: Identification of cisplatin-resistance related genes in head and neck squamous cell carcinoma. Internat j cancer 2010,126(2):437–49. 10.1002/ijc.24704

    Article  CAS  Google Scholar 

  33. Furlan F, Galbiati C, Jorgensen NR, Jensen JE, Mrak E, Rubinacci A, Tolatta F, Verde P, Blasi F: Urokinase plasminogen activator receptor affects bone homeostasis by regulating osteoblast and osteoclast function. J bone min res 2007,22(9):1387–96. 10.1359/jbmr.070516

    Article  CAS  Google Scholar 

  34. Scotto L, Narayan G, Nandula SV, Subramaniyam S, Kaufmann AM, Wright JD, Pothuri B, Mansukhani M, Schneider A, Arias-Pulido H, Murty VV: Integrative genomics analysis of chromosome 5p gain in cervical cancer reveals target over-expressed genes, including Drosha. Molec Cancer 2008, 7: 58. 10.1186/1476-4598-7-58

    Article  Google Scholar 

  35. Walker G, MacLeod K, Williams AR, Cameron DA, Smyth JF, Langdon SP: Estrogen-regulated gene expression predicts response to endocrine therapy in patients with ovarian cancer. Gynecol oncol 2007,106(3):461–8. 10.1016/j.ygyno.2007.05.009

    Article  CAS  PubMed  Google Scholar 

  36. Safavy A, Raisch KP, Matusiak D, Bhatnagar S, Helson L: Single-drug multigland conjugates: synthesis and preliminary cytotoxicity evaluation of a paclitaxel-dipeptide "scorpion" molecule. Bioconj chem 2006,17(3):565–70. 10.1021/bc050224c

    Article  CAS  Google Scholar 

  37. Peiffer I, Belhomme D, Barbet R, Haydont V, Zhou YP, Fortunel NO, Li M, Hatzfeld A, Fabini JN, Hatzfeld JA: Simultaneous differentiation of endothelial and trophoblastic cells derived from human embryonic stem cells. Stem cells dev 2007,16(3):393–402. 10.1089/scd.2006.0013

    Article  CAS  PubMed  Google Scholar 

  38. Shang X, Vasudevan SA, Yu Y, Ge N, Ludwig AD, Wesson CL, Wang K, Burlingame SM, Zhao YJ, Rao PH, Lu X, Russell HV, Okcu MF, Hicks MJ, Shohet JM, Donehower LA, Nuchtern JG, Yang J: Dual-specificity phosphatase 26 is a novel p53 phosphatase and inhibits p53 tumor suppressor functions in human neuroblastoma. Oncogene 2010,29(35):4938–46. 10.1038/onc.2010.244

    Article  CAS  PubMed  Google Scholar 

  39. Hur W, Rhim H, Jung CK, Kim JD, Bae SH, Jang LW, Yang JM, Oh ST, Kim DG, Wang HJ, Lee SB, Yoon SK: Sox4 overexpression regulates the p53-mediated apoptosis in hepatocellular carcinoma clinical implications and functional analysis in vitro. Carcinog 2010,31(7):1298–307. 10.1093/carcin/bgq072

    Article  CAS  Google Scholar 

  40. Lecomte S, Desmots F, Le Masson F, Le Goff P, Mechel D, Christians ES, Le Drean Y: Roles of heat shock factor 1 and 2 in response to proteasome inhibition: consequences on p53 stability. Oncogene 2010,29(29):4216–24. 10.1038/onc.2010.171

    Article  CAS  PubMed  Google Scholar 

  41. Asefa B, Dermott JM, Kaldis P, Stefanisko K, Garfinkel JR: p205, a potential tumor suppressor, inhibits cell proliferation via multiple pathways of cell cycle. FEBS letters 2006,580(5):1205–14. 10.1016/j.febslet.2006.01.032

    Article  CAS  PubMed  Google Scholar 

  42. Yeom KH, Lee Y, Han J, Suh MR, Kim VN: Characterisation of DGCR8/Pasha, the essential cofactor for Drosha in primary miRNA processing. Nucl acid res 2006,34(16):4622–9. 10.1093/nar/gkl458

    Article  CAS  Google Scholar 

  43. Mellett M, Atzei P, Jackson R, O'Neill LA, Moynagh PN: Mal mediates TLR-induced activation of CREB and expression of IL-10. J Immunol 2011,186(8):4925–35. 10.4049/jimmunol.1002739

    Article  CAS  PubMed  Google Scholar 

  44. Chen MC, Hwang MJ, Chou YC, Chen WH, Cheng G, Nakano H, Luh TY, Mai SC, Hsieh SL: The role of apoptosis signal-regulating kinase 1 in lymphotoxin-beta receptor-mediated cell death. J biol chem 2003,278(18):16073–81. 10.1074/jbc.M208661200

    Article  CAS  PubMed  Google Scholar 

  45. Tennant DA, Gottlieb E: HIF prolyl hydroxylase-3 mediates alpha-ketoglutarate-induced apoptosis and tumor suppression. J molec med (Berlin) 2010,88(8):839–49. 10.1007/s00109-010-0627-0

    Article  CAS  Google Scholar 

  46. Uhlenbrock K, Huber J, Ardati A, Busch AE, Kostenis E: Fluid sheer stress differentially regulates gpr3, gpr6 and grp12 expression in human umbilical vein endothelial cells. Cell physiol biochem 2003,13(2):75–84. 10.1159/000070251

    Article  CAS  PubMed  Google Scholar 

  47. Zhang X, Azhar G, Helms S, Zhong Y, Wei JY: Identification of a subunit of NADH-dehydrogenase as a p49/STRAP-binding protein. BMC Cell biology 2003, 9: 8.

    Article  Google Scholar 

  48. Chen M, Singh A, Xiao F, Dringenberg U, Wang J, Engelhardt R, Teruva S, Rubio-Aliaga I, Nassl AM, Kottra G, Daniel H, Seidler U: Gene ablation for PEPT1 in mice abolishes the effects of dipeptides on small intestinal fluid absorption, short-circuit current and intracellular pH. Am j physiol: Gastro liver physiol 2010,229(1):G265–74.

    Google Scholar 

  49. Jurasz P, Yurkova N, Kirshenbaum L, Stewart DJ: VEGF masks BNIP3-mediated apoptosis of hypoxic endothelial cells. Angiogenesis 2011,14(92):199–207.

    Article  CAS  PubMed  Google Scholar 

  50. Ramasamy S, Wang H, Quach HN, Sampath K: Zebrafish Straufen1 and Straufen2 are required for the survival and migration of primordial germ cells. Dev biol 2006,292(2):393–406. 10.1016/j.ydbio.2006.01.014

    Article  CAS  PubMed  Google Scholar 

  51. Sun R, Su Y, Zhao X, Qi J, Luo X, Yang Z, Yao Y, Luo X, Xia Z: Human calcium/calmodulin-dependent serine protein kinase regulates the expression of p21 via the E2A transcription factor. Biochemical j 2009,419(2):457–66. 10.1042/BJ20080515

    Article  CAS  Google Scholar 

  52. Jaffer ZM, Chernoff J: p21-activated kinases: three more join the Pak. Intl j biochem cell biol 2002,34(7):713–7. 10.1016/S1357-2725(01)00158-3

    Article  CAS  Google Scholar 

  53. Ivanovska I, Ball AS, Diaz RL, Magnus JF, Kibukawa M, Schelter JM, Kobayashi SV, Lim L, Burchard J, Jackson AL, Linsley PS, Cleary MA: MicroRNAs in the miR-106b family regulate p21/CDKN1A and promote cell cycle progression. Molec cell biol 2008, 28: 2167–2174. 10.1128/MCB.01977-07

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  54. Ahmed N, Abubaker K, Findlay J, Quinn M: Epithelial mesenchymal transition and cancer stem cell-like phenotypes facilitate chemoresistance in recurrent ovarian cancer. Curr cancer drug targets 2010, 10: 268–278. 10.2174/156800910791190175

    Article  CAS  PubMed  Google Scholar 

  55. Park JT, Chen X, Trope CG, Davidons B, Shih IM, Wang TL: Notch3 Overexpression Is Related to the Recurrence of Ovarian Cancer and Confers Resistance to Carboplatin. Am J Path 2010. Epub Jul 29

    Google Scholar 

  56. Hennessy EJ, Parker AE, O'Neill LA: Targeting Toll-like receptors: emerging therapeutics. Nat rev drug disc 2010, 9: 293–307. 10.1038/nrd3203

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The authors wish to acknowledge the support of Cancer Research Ireland and The Emer Casey Foundation.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Michael F Gallagher.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

MFG designed and carried out the comparisons and was the main author of the manuscript. CCBBH, MFG and PCS carried out microarray experiments and analysis. BF conducted statistical analysis. AL, SAO'T, RJF, SAE, CMM, OMS were involved in the original experiments and contributed to data analysis and interpretation. CDS contributed to manuscript design. JJO'L oversaw the project. All authors have read and approved the final manuscript.

Electronic supplementary material

13048_2011_82_MOESM1_ESM.XLS

Additional file 1: Microarray and functional relationship analysis of early differentiation response of SCC-PSA1 cells. This file contains genelists of genes upregulated and downregulated by SCC-PSA1 cells in response to differentiation stimulus. This file also contains results of analysis to identify functional relationships between these genes. (XLS 114 KB)

13048_2011_82_MOESM2_ESM.XLS

Additional file 2: Microarray and functional relationship analysis of early differentiation response of NULLI-SCC cells. This file contains genelists of genes upregulated and downregulated by NULLI-SCC cells in response to differentiation stimulus. This file also contains results of analysis to identify functional relationships between these genes. (XLS 56 KB)

13048_2011_82_MOESM3_ESM.XLS

Additional file 3: Microarray and functional relationship analysis of early differentiation response of mES cells. This file contains genelists of genes upregulated and downregulated by mES cells in response to differentiation stimulus. This file also contains results of analysis to identify functional relationships between these genes. (XLS 195 KB)

13048_2011_82_MOESM4_ESM.XLS

Additional file 4: Microarray and functional relationship analysis of primary versus recurrent tumor samples. This file contains genelists of genes upregulated and downregulated in primary versus recurrent tumor samples. This file also contains results of analysis to identify functional relationships between these genes. (XLS 234 KB)

13048_2011_82_MOESM5_ESM.XLS

Additional file 5: Comparison of microRNA expression in hEC and primary versus recurrent tumor samples. This file details the relative expression patterns and levels of microRNAs in hEC cells and primary versus recurrent tumor samples and the pathway associations of genes targeted by these microRNAs. (XLS 40 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Gallagher, M.F., Heffron, C.C., Laios, A. et al. Suppression of cancer stemness p21-regulating mRNA and microRNA signatures in recurrent ovarian cancer patient samples. J Ovarian Res 5, 2 (2012). https://doi.org/10.1186/1757-2215-5-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1757-2215-5-2

Keywords