HOTAIR and its surrogate DNA methylation signature indicate carboplatin resistance in ovarian cancer
Genome Medicine volume 7, Article number: 108 (2015)
Understanding carboplatin resistance in ovarian cancer is critical for the improvement of patients’ lives. Multipotent mesenchymal stem cells or an aggravated epithelial to mesenchymal transition phenotype of a cancer are integrally involved in pathways conferring chemo-resistance. Long non-coding RNA HOTAIR (HOX transcript antisense intergenic RNA) is involved in mesenchymal stem cell fate and cancer biology.
We analyzed HOTAIR expression and associated surrogate DNA methylation (DNAme) in 134 primary ovarian cancer cases (63 received carboplatin, 55 received cisplatin and 16 no chemotherapy). We validated our findings by HOTAIR expression and DNAme analysis in a multicentre setting of five additional sets, encompassing 946 ovarian cancers. Chemo-sensitivity has been assessed in cell culture experiments.
HOTAIR expression was significantly associated with poor survival in carboplatin-treated patients with adjusted hazard ratios for death of 3.64 (95 % confidence interval [CI] 1.78–7.42; P < 0.001) in the discovery and 1.63 (95 % CI 1.04–2.56; P = 0.032) in the validation set. This effect was not seen in patients who did not receive carboplatin (0.97 [95 % CI 0.52–1.80; P = 0.932]). HOTAIR expression or its surrogate DNAme signature predicted poor outcome in all additional sets of carboplatin-treated ovarian cancer patients while HOTAIR expressors responded preferentially to cisplatin (multivariate interaction P = 0.008).
Non-coding RNA HOTAIR or its more stable DNAme surrogate may indicate the presence of a subset of cells which confer resistance to carboplatin and can serve as (1) a marker to personalise treatment and (2) a novel target to overcome carboplatin resistance.
Late stage presentation aside, carboplatin resistance in ovarian cancer is the key obstacle to improving survival in this disease . The observation that re-treatment with platinum-based drugs 6–12 months after primary response proved to be successful in a certain percentage of patients  is consistent with the idea that platinum sensitivity can be modulated by both cancer cell-autonomous and non-autonomous factors. For both these factors stromal/mesenchymal differentiation is crucially important. Epithelial–mesenchymal transition (EMT) in ovarian cancer cells is associated with platinum resistance [3–5]. On the other hand the cancer cell-autonomous tumour stroma and mesenchymal stem cells (MSCs) — mainly recruited from the bone marrow  — might play an important role in ovarian cancer biology [7, 8]. Recently, bone marrow-derived MSCs and embryonic fibroblasts, but not more extensively differentiated stromal cells, have been shown to induce platinum resistance in ovarian cancer .
Long non-coding RNAs are known to epigenetically remodel chromatin states and influence gene transcription in normal and cancer tissue towards stromal/mesenchymal differentiation [10–13]. Aberrant expression of non-coding RNAs has been observed in numerous diseases, including cancer , yet their precise contribution to disease aetiology and biology is far from clear. HOX antisense transcript intergenic RNA (HOTAIR), transcribed from the HOXC locus, represses transcription by recruiting polycomb repressive complex 2 (PRC2) to specific polycomb group target (PCGT) genes, in particular to those normally targeted by PRC2 in embryonic fibroblasts .
In stem cells, PCGTs are repressed through PRC2 occupancy and PCGTs important for specialised cell identities become de-repressed upon differentiation [15, 16]. We and others have shown that the promoters of these stem cell PCGTs become methylated and silenced in cancer [17–20]. It was recently reported that the expression of HOTAIR is increased in various cancer entities and that high levels of expression correlate with cancer invasiveness, metastases and poor prognosis [10, 21]. A recent systematic review of 19 papers (including a total of 2255 patients) demonstrates consistently that HOTAIR expression is a poor prognostic marker across a large set of cancers . It is unclear, however, whether HOTAIR is associated with an aberrant DNA methylation profile in cancer and whether this robust DNA-based imprint mediates resistance to specific drugs.
Here we tested the hypothesis that HOTAIR RNA expression, or a HOTAIR-associated DNA methylation (DNAme) signature, as surrogates for mesenchymal differentiation, serve as markers for carboplatin resistance in primary ovarian cancer.
Ovarian cancer data sets
We analyzed six different data sets, details of which are provided in Additional file 1.
The first data set consisted of primary ovarian cancer samples (n = 134, 24–87 years, median 62.7 years at diagnosis; Additional file 2) treated at the Innsbruck Medical University, denoted “INNSBRUCK”. Clinicopathological features are shown in Table 1. The study was approved by the ethical committee of the Medical University Innsbruck (reference number UN4044). For the majority of patients exemption from obtaining informed consent was received as the majority of ovarian cancer patients were dead at the time the application was evaluated. The median survival time was 3.8 years. HOTAIR expression was measured in all 134 samples. DNA methylation data are available as Additional file 3.
The second data set consisted of primary ovarian cancer samples (n = 175, 21–83 years, median 60.0 years at diagnosis; Additional file 4) treated at the University Medical Center in Groningen , denoted “GRONINGEN”. The median survival time in this set, which consisted only of stage III/IV patients, was 2.1 years. For these 175 samples, we measured HOTAIR expression. Of the 175 samples, 157 received carboplatin only, whilst 18 received cisplatin instead. For 114 of these samples there were matched mRNA array expression profiles available (Operon Human v.3 ~ 35 K 70-mer two-color oligonucleotide arrays, Gene Expression Omnibus accession [GEO:GSE13876]). Patients gave informed consent for collection and storage of tissue samples in a tissue bank for future research. All relevant patient data were retrieved and transferred into an anonymous, password-protected database. The patients’ identity was protected by study-specific, unique patient codes and their true identity was only known to two dedicated data managers. According to Dutch regulations, these precautions meant no further institutional review board approval was needed.
The third data set consisted of primary ovarian cancer samples (n = 342, serous ovarian cancers, median 58 years at diagnosis; 316 received carboplatin-based chemotherapy with the rest received cisplatin or were untreated), analyzed within The Cancer Genome Atlas (TCGA) program, and for which Illumina Infinium 27 k DNAme data were publicly available  (Additional file 5), denoted “TCGA”. The median survival time was 2.6 years.
The fourth data set consisted of primary ovarian cancer samples from three European Cancer centres (Leuven, Oslo, Rotterdam; n = 206, median 58 years at diagnosis; Additional file 6), denoted “EUROPE”. The median survival time was 3.4 years. For this data set, Illumina Infinium 450 k DNAme data were available for 121 carboplatin- and 85 cisplatin-treated patients. Data are deposited in the GEO, accession [GEO:GSE72021]. The study from Rotterdam has been approved by the local medical ethics committee (MEC-2008-183), performed in accordance with the Code of Conduct of the Federation of Medical Scientific Societies in the Netherlands. The Regional Committee for Medical Research Ethics in Norway approved the study (for patients diagnosed before 2007, exemption from obtaining informed consent was received as the majority of ovarian cancer patients were dead at the time the application was evaluated; patients diagnosed after 2007 signed general consent allowing for use of the tumours for research purposes). Written informed consent for the use of tumour tissue and prospective clinical data collection was obtained from all patients and approved by the Leuven ethics committee.
The fifth data set consisted of primary ovarian cancer samples from Bergen (n = 49) with 40 receiving carboplatin and 9 untreated (no chemotherapy), denoted “BERGEN”. Patients were included in the study after written informed consent, approved by the Regional Research Ethics Committee in Medicine. For samples from this cohort, we measured HOTAIR expression.
The sixth data set consisted of primary ovarian cancer samples from Rochester-Mayo (n = 174), denoted “ROCHESTER-MAYO”. All 174 patients received carboplatin and for these samples we measured DNAme using Illumina 450 k beadarrays. All patients gave informed consent and the Mayo Clinic Institutional Review Board approved the study. The data are available from Dr Ellen Goode at the Department of Health Sciences Research, Mayo Clinic, Rochester, USA, upon request.
Our research conformed to the Helsinki Declaration.
Total RNA was extracted by the acid guanidium thiocyanate-phenol-chloroform method . Reverse transcription of RNA was performed as previously described . Primers and probes for HOTAIR were designed using Primer Express (Applied Biosystems, Foster City, CA, USA). Samples in which HOTAIR was not amplified by real-time PCR after 45 cycles were classified as negative (HOTAIR-ve; Additional file 1).
DNA methylation analysis
DNA was isolated from tissue samples using the Qiagen DNeasy Blood and Tissue Kit (Qiagen Ltd, UK, 69506) and 600 ng was bisulphite converted using the Zymo Methylation Kit (Zymo Research Inc, USA, D5004/8). Genome-wide methylation analysis was performed using the Illumina Infinium Methylation 27 K or 450 K beadchip (Illumina Inc., USA, WG-311-1201 and WG-314-1003). Analysis and quality control were performed as previously described [19, 27, 28].
HOTAIR overexpression in ovarian cancer cell lines
The SKOV3IP cells were stably transduced with HOTAIR and LacZ constructs, kindly provided by Dr Chang (Stanford)  and single clones of HOTAIR/LacZ overexpressing cells were used for experiments (Additional file 1). Cells were treated with cisplatin (0.5–18 μM) or carboplatin (10–160 μM) for 3 days and analysed by the cell survival MTT assay (Sigma).
To test for differences in categorical variables, we used the Chi square test. Impact of HOTAIR expression on ovarian cancer survival was ascertained using log rank test and Kaplan Meier curves. To demonstrate the effect of HOTAIR expression on DNAme patterns, we first performed univariate analysis to rank CpGs according to their association with HOTAIR expression (Additional file 1). We used 10-fold internal cross-validations to identify an optimal HOTAIR-associated DNAme signature consisting of 67 CpGs (Additional file 7) at an estimated false discovery rate (FDR) of approximately 0.17. To build a single-sample classifier from this signature, a cutoff was optimized using receiver operating characteristic (ROC) analysis. This same cutoff was then used to assign samples from independent cohorts into two groups exhibiting high and low correlations with the HOTAIR DNAme signature. For the PCGT enrichment analysis we relaxed the threshold of the DNAme signature to include the top 500 ranked CpGs (FDR < 0.3), and divided the 500 into the 233 which were hypermethylated and the 267 which were hypomethylated in high HOTAIR expressors.
HOTAIR expression in primary ovarian cancer is not associated with clinicopathological features
We analyzed HOTAIR expression in 134 primary ovarian cancer samples (INNSBRUCK) and found that 72 were positive and 62 were negative for HOTAIR. Consistent with our findings that only a subset of cells may express HOTAIR in ovarian cancer tissue, the cycle threshold (ct) values in the positive samples were very low (mean ct 37.4) compared with the TBP reference gene (mean ct 27.4). On assessing the clinicopathological characteristics of the sample set, survival status was the only characteristic that was significantly associated with HOTAIR expression, i.e., those patients whose tumours were HOTAIR + ve had a poor outcome compared with HOTAIR-ve tumours (P = 0.023, Table 1).
Association between HOTAIR expression and poor outcome is restricted to carboplatin-treated patients
In order to test the hypothesis that HOTAIR is linked to carboplatin resistance we analyzed survival in patients who received carboplatin both alone or as part of a combination therapy (‘carboplatin’ group) compared with those who received cisplatin or no chemotherapy (‘no carboplatin group’) in the INNSBRUCK set. We note that there was no difference in any clinicopathological feature or survival between the two groups (Additional file 2). HOTAIR expression was significantly associated with both risk of relapse (hazard ratio (HR) 4.46 [P < 0.001] and 3.38 [P = 0.003] in uni- and multivariate analysis, respectively) and of death (HR 4.02 [P < 0.001] and 3.64 [P < 0.001] in uni- and multivariate analysis, respectively) in the carboplatin group (Fig. 1a; Additional file 8), whereas HOTAIR expression was not associated with survival in the ‘no carboplatin’ group (Fig. 1b; Additional file 8).
Only 34 % of patients who received carboplatin and whose tumours expressed HOTAIR survived the first 3 years post-diagnosis, whereas 85 % of those patients who received carboplatin, but whose tumours were HOTAIR-ve, survived during the same time period (Fig. 1a). Nearly half of the carboplatin-treated patients also received paclitaxel, and in both groups, i.e., single agent carboplatin and combined carboplatin/paclitaxel, HOTAIR was significantly associated with poor outcome (log rank P value = 0.006 and 0.003, respectively), again indicating the interaction of HOTAIR with carboplatin but not with cisplatin or paclitaxel.
To validate these findings we analyzed 175 ovarian cancer samples from Groningen  and 49 samples from Bergen and again confirmed that HOTAIR expression is a poor prognostic factor specifically in carboplatin-treated patients (Fig. 1c, d; Additional files 9 and 10).
HOTAIR expression is associated with a DNAme profile enriched for PCGTs and associated with multipotent MSCs
As HOTAIR is known to modulate chromatin, in particular at PCGTs, we asked if DNAme differed between HOTAIR-expressing and non-expressing ovarian cancer samples. Since HOTAIR expression was associated with carboplatin resistance, we restricted the analysis to the 63 carboplatin-treated patients in the INNSNRUCK set. Of these 63 samples, 35 expressed HOTAIR whilst 28 did not. We identified a 67-CpG DNAme signature representing a statistically significant association with HOTAIR expression. Boxplots of beta methylation values of the 67 CpGs confirmed the relatively large differences in methylation between HOTAIR expressors and non-expressors (Additional file 11). We observed that many of the top CpGs mapped to PCGTs (Fig. 2a) and found that PCGTs (defined as PRC2 targets in both human embryonic stem cells and human embryonic fibroblasts, but not PRC2 targets in breast cancer cells) were highly enriched among CpGs hypermethylated in HOTAIR expressors (Additional file 12). HOTAIR expression was also associated with lower expression of PCGTs, in particular human embryonic fibroblast PCGTs, in the GRONINGEN set , supporting the role of HOTAIR as an epigenetic regulator of MSCs in ovarian cancer (Additional file 13). By correlating the 67-CpG DNAme signature to the methylation profile of any given sample, a correlation score was obtained which can be viewed as a DNA-based surrogate for HOTAIR expression. We first evaluated this correlation score in early passage (multipotent) MSCs, late passage (more differentiated) MSCs, reprogrammed MSCs, embryonic stem cells  and ovarian cancer cell lines, demonstrating that the HOTAIR DNAme signature is likely to be a surrogate marker for either multipotent MSCs within the ovarian cancer tissue, or for ovarian cancer cells with an increased tendency to undergo EMT (Fig. 2b; Additional file 14).
The HOTAIR-DNAme signature predicts survival outcome
Given the relationship between HOTAIR expression and carboplatin resistance, we next checked that the HOTAIR DNAme signature would be similarly associated with outcome in the 63 carboplatin-treated subgroup in the INNSBRUCK set. As expected, the score obtained by correlating the sample-specific DNAme profile to the HOTAIR DNAme signature was predictive of carboplatin resistance in univariate as well as in multivariate Cox-regression analyses adjusted for age, stage and size of residual tumour (Additional file 15). In order to build a single-sample classifier we also optimized a cutoff on the correlation score of the ovarian cancer tissue samples using ROC analysis to ensure approximately 80 % sensitivity and 80 % specificity between the DNAme signature and HOTAIR expression (Additional file 16). Dividing the samples into two groups based on this cutoff further confirmed a significant difference in survival rates of these two subgroups of patients, i.e., samples with a high DNAme signature score had a HR of 3.07 (P = 0.001) for death relative to those with a low surrogate score (Fig. 2c).
Validation of HOTAIR DNAme signature in three independent large cohorts
To validate our HOTAIR DNAme signature, we first tested it in an independent large data set of serous ovarian cancers from TCGA set, consisting of 316 patients who received carboplatin-based therapy and 26 who received cisplatin. For each of the carboplatin-treated patients, we computed a carboplatin resistance score by correlating the DNAme profile of the tumour to the previously determined 67-CpG DNAme HOTAIR signature. This score predicted outcome in both univariate as well as multivariate Cox-regression analyses (Additional file 15), and was a much stronger predictor of outcome than scores constructed using random signatures (Additional file 17). Using the previously determined cutoff to assign samples into high and low DNAme signature score groups further demonstrated the robustness of the predictive classifier in the carboplatin-treated subgroup (Fig. 2d). Of note, the prediction obtained using the HOTAIR DNAme signature outperforms the classification obtained using either mRNA or microRNA expression predictors, as reported in TCGA study . In contrast to the carboplatin-treated group, the HOTAIR DNAme signature was not predictive of outcome in the 26 patients in TCGA set who did not receive carboplatin (Fig. 2e).
We further tested the HOTAIR DNAme signature in an independent set of 121 carboplatin-treated patients (EUROPE set) with DNAme profiles generated using a different assay (Illumina Infinium Human Methylation 450 k). In order to more rigorously assess the specificity of the signature, it was also tested in a further 85 cisplatin-treated patients from the same EUROPE cohort which were also profiled with the same 450 k technology. Once again, the DNAme signature-based surrogate scores for HOTAIR expression correctly predicted carboplatin resistance (Fig. 2f; Additional file 18), with no association observed in the cisplatin-treated subgroup (Fig. 2g; Additional file 18).
Further strengthening the robustness of the HOTAIR DNAme signature, we found it to be predictive of carboplatin resistance in an additional independent set of 174 carboplatin-treated ovarian cancer patients (ROCHESTER-MAYO set), which had also been profiled with Illumina 450 k DNAme bead arrays (Additional file 19).
HOTAIR expressors respond preferentially to cisplatin-based chemotherapy
In all data sets analysed, we observed a consistent trend for HOTAIR expression, or DNAme-based surrogate HOTAIR expression, to be associated preferentially with cisplatin response, although statistical significance was not observed in individual data sets. Thus, to investigate this further we used a meta-analysis approach and asked if the type of chemotherapy received was associated with a different response in HOTAIR expressors compared with non-expressors (Table 2). In all data sets the risk of death was lower in HOTAIR expressors who received cisplatin compared with those who received carboplatin-based therapies, whereas for non-HOTAIR expressors the opposite pattern was observed (Table 2; Additional files 20 and 21). Using a combined probability Fisher test in a meta-analysis over all data sets, we found a highly significant interaction between chemotherapy type received and HOTAIR expression in dictating response to treatment (P < 0.001; Table 2), which was retained in multivariate analysis (P = 0.008; Table 2).
Effect of HOTAIR expression and platinum sensitivity in ovarian cancer cell line
In order to test whether HOTAIR expression modulates response to carboplatin we used SKOV3IP ovarian cancer cells, which are sensitive to platinum-based chemotherapy and do not express HOTAIR. We observed that overexpression of HOTAIR in this particular cell line reduces only sensitivity to carboplatin but not cisplatin (Fig. 3). The half maximal inhibitory concentration (IC50) shifted from 30 to 60 μM for carboplatin and from 3 to 3.5 μM for cisplatin in LacZ and HOTAIR-expressing SKOV3IP cells, respectively. HOTAIR expression on two other cell lines had either no effect (A2780) or increased cisplatin sensitivity (OVCAR8) (Additional file 22).
Here we have shown that women with ovarian tumours expressing HOTAIR RNA, or an equivalent HOTAIR-associated DNAme signature, experience a poorer survival outcome post carboplatin-based chemotherapy compared with HOTAIR-ve tumours.
These data offer solid evidence for two novel concepts. First, we demonstrate that not only a non-coding RNA, but importantly also a presumed downstream effect, represented by a specific DNAme signature, reproducibly predicts cancer outcome. This concept is appealing, because (i) a DNA-based biomarker is more stable and (ii) would provide a spatially and temporally more comprehensive surrogate for cancer biology compared with a snap-shot RNA assessment. Second, our data very much challenge the dogma that cisplatin and carboplatin have the same effect on ovarian cancer and that parameters that allow for discrimination between patients that benefit form carboplatin and those that benefit form cisplatin would not exist.
A recent study from Roodhart et al.  highlights the complexity of the chemo-resistance response which integrally involves cancer cell non-autonomous factors. The authors used a mouse model and demonstrated that cells with a multilineage differentiation potential, such as embryonic fibroblasts, become activated during treatment with platinum analogs and secrete factors systemically that protect tumour cells against platinum chemotherapeutics. Roodhart et al.  demonstrated this for both carboplatin and cisplatin whereas our data — based on 1080 human ovarian cancer samples — demonstrate that HOTAIR and its surrogate DNAme signature are associated only with carboplatin resistance and not cisplatin resistance.
Although in the past cisplatin and carboplatin have been considered to be nearly identical drugs and prospective randomized trials showed similar overall survival for cisplatin- and carboplatin-based regimens in ovarian cancer [30, 31], both drugs have substantially different side-effect profiles with higher rates of nausea, vomiting and renal toxicity for cisplatin and thrombocytopenia for carboplatin . Moreover, cisplatin and carboplatin are known to differ in their cell membrane transport characteristics . This is noteworthy because cell membrane transport proteins have been shown to be critical determinants of platinum drug sensitivity/resistance, possibly as a result of secreted factors from MSCs which can affect transport characteristics. Lending further credence to our data, in cervical cancer (a disease which is known to expresses high levels of HOTAIR ) carboplatin has been repeatedly reported to be a less effective platinum analog than cisplatin [35–37].
Although we provide unprecedented strong evidence for a non-coding RNA and its DNAme surrogate signature to be a predictive and prognostic marker in ovarian cancer, there are several limitations to our study. First, although in situ hybridisation for HOTAIR has been successfully established for cell lines, we were not able (despite substantial efforts; data not shown) to determine the specific subset of HOTAIR-expressing cells within the bulk tumour sample. Hence, at this stage we can only state that HOTAIR expression serves as an excellent surrogate for “mesenchyme-ness” of a cancer, not knowing whether this reflects the presence of MSCs in the tumour stroma or whether it is a reflection of the number of cancer cells that have undergone EMT, or a combination of both. Second, we don’t provide a mechanistic model as to why and how HOTAIR modulates carboplatin and cisplatin response differently. Whereas recent evidence shows that HOTAIR promotes proliferation by modulating cell cycle and apoptosis , no evidence exists to show that this has an impact on platinum resistance. Our findings are consistent with the view that HOTAIR modulates the epigenome at the level of the DNA methylome in both cancer cells and tumour stroma. Although still speculative, it is likely that the differential response of carboplatin and cisplatin is due to underlying differences in how the MSC biology of tumour stroma and the EMT characteristics of cancer cells affect the two drugs. Support for this view comes from recent data (unpublished): we have analysed MSCs and modulated expression of HOTAIR. Overexpression and knockdown of HOTAIR inhibited or stimulated, respectively, in vitro differentiation of MSCs. Modification of HOTAIR expression evoked consistent effects on gene expression, particularly in polycomb group target genes and genes involved in cancer. Furthermore, overexpression and knockdown of HOTAIR resulted in DNAme changes that are enriched in HOTAIR binding sites.
Despite these limitations, our findings have a number of immediate clinical implications. We provide a solid rationale for prospective randomized clinical trials — ideally in a neo-adjuvant setting — to assess whether the HOTAIR DNAme signature is an appropriate tool to stratify women with ovarian cancer (and possibly also other cancers) into groups which benefit preferentially from cisplatin or from carboplatin treatment. Strategies to reduce HOTAIR activity (i.e., by intra-peritoneal HOTAIR small interfering RNA) may lead to a novel strategy to (re)sensitize cancers to chemotherapy.
Our data demonstrate that HOTAIR and its surrogate DNAme signature play a crucially important role in ovarian cancer biology and provide novel leads to revisit the clinically important field of platinum resistance in this disease.
false discovery rate
Gene Expression Omnibus
- HOTAIR :
HOX transcript antisense intergenic RNA
mesenchymal stem cell
polycomb group target
polycomb repressor complex 2
receiver operating characteristic
The Cancer Genome Atlas.
Markman M. Antineoplastic agents in the management of ovarian cancer: current status and emerging therapeutic strategies. Trends Pharmacol Sci. 2008;29:515–9.
Colombo N, Gore M. Treatment of recurrent ovarian cancer relapsing 6–12 months post platinum-based chemotherapy. Crit Rev Oncol Hematol. 2007;64:129–38.
Marchini S, Fruscio R, Clivio L, Beltrame L, Porcu L, Fuso NI, et al. Resistance to platinum-based chemotherapy is associated with epithelial to mesenchymal transition in epithelial ovarian cancer. Eur J Cancer. 2013;49:520–30.
Helleman J, Smid M, Jansen MP, van der Burg ME, Berns EM. Pathway analysis of gene lists associated with platinum-based chemotherapy resistance in ovarian cancer: the big picture. Gynecol Oncol. 2010;117:170–6.
Gupta N, Xu Z, El-Sehemy A, Steed H, Fu Y. Notch3 induces epithelial-mesenchymal transition and attenuates carboplatin-induced apoptosis in ovarian cancer cells. Gynecol Oncol. 2013;130:200–6.
Kidd S, Spaeth E, Watson K, Burks J, Lu H, Klopp A, et al. Origins of the tumor microenvironment: quantitative assessment of adipose-derived and bone marrow-derived stroma. PLoS One. 2012;7:e30563.
McLean K, Gong Y, Choi Y, Deng N, Yang K, Bai S, et al. Human ovarian carcinoma-associated mesenchymal stem cells regulate cancer stem cells and tumorigenesis via altered BMP production. J Clin Invest. 2011;121:3206–19.
Lis R, Touboul C, Raynaud CM, Malek JA, Suhre K, Mirshahi M, et al. Mesenchymal cell interaction with ovarian cancer cells triggers pro-metastatic properties. PLoS One. 2012;7:e38340.
Roodhart JM, Daenen LG, Stigter EC, Prins HJ, Gerrits J, Houthuijzen JM, et al. Mesenchymal stem cells induce resistance to chemotherapy through the release of platinum-induced fatty acids. Cancer Cell. 2011;20:370–83.
Gupta RA, Shah N, Wang KC, Kim J, Horlings HM, Wong DJ, et al. Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature. 2010;464:1071–6.
Padua AC, Fonseca AS, Muys BR, de Barros E, Lima Bueno R, Burger MC, et al. Brief report: The lincRNA Hotair is required for epithelial-to-mesenchymal transition and stemness maintenance of cancer cell lines. Stem Cells. 2013;31:2827–32.
Wang KC, Chang HY. Molecular mechanisms of long noncoding RNAs. Mol Cell. 2011;43:904–14.
Xu ZY, Yu QM, Du YA, Yang LT, Dong RZ, Huang L, et al. Knockdown of long non-coding RNA HOTAIR suppresses tumor invasion and reverses epithelial-mesenchymal transition in gastric cancer. Int J Biol Sci. 2013;9:587–97.
Tsai MC, Spitale RC, Chang HY. Long intergenic noncoding RNAs: new links in cancer progression. Cancer Res. 2011;71:3–7.
Bernstein BE, Mikkelsen TS, Xie X, Kamal M, Huebert DJ, Cuff J, et al. A bivalent chromatin structure marks key developmental genes in embryonic stem cells. Cell. 2006;125:315–26.
Lee TI, Jenner RG, Boyer LA, Guenther MG, Levine SS, Kumar RM, et al. Control of developmental regulators by Polycomb in human embryonic stem cells. Cell. 2006;125:301–13.
Ohm JE, McGarvey KM, Yu X, Cheng L, Schuebel KE, Cope L, et al. A stem cell-like chromatin pattern may predispose tumor suppressor genes to DNA hypermethylation and heritable silencing. Nat Genet. 2007;39:237–42.
Schlesinger Y, Straussman R, Keshet I, Farkash S, Hecht M, Zimmerman J, et al. Polycomb-mediated methylation on Lys27 of histone H3 pre-marks genes for de novo methylation in cancer. Nat Genet. 2007;39:232–6.
Teschendorff AE, Menon U, Gentry-Maharaj A, Ramus SJ, Weisenberger DJ, Shen H, et al. Age-dependent DNA methylation of genes that are suppressed in stem cells is a hallmark of cancer. Genome Res. 2010;20:440–6.
Widschwendter M, Fiegl H, Egle D, Mueller-Holzner E, Spizzo G, Marth C, et al. Epigenetic stem cell signature in cancer. Nat Genet. 2007;39:157–8.
Qiu JJ, Lin YY, Ye LC, Ding JX, Feng WW, Jin HY, et al. Overexpression of long non-coding RNA HOTAIR predicts poor patient prognosis and promotes tumor metastasis in epithelial ovarian cancer. Gynecol Oncol. 2014;134:121–8.
Deng Q, Sun H, He B, Pan Y, Gao T, Chen J, et al. Prognostic value of long non-coding RNA HOTAIR in various cancers. PLoS One. 2014;9:e110059.
Crijns AP, Fehrmann RS, de Jong S, Gerbens F, Meersma GJ, Klip HG, et al. Survival-related profile, pathways, and transcription factors in ovarian cancer. PLoS Med. 2009;6:e24.
Bell D, Berchuck A, Birrer M, Chien J, Cramer DW, Dao F, et al. Integrated genomic analyses of ovarian carcinoma. Nature. 2011;474:609–15.
Widschwendter M, Berger J, Hermann M, Muller HM, Amberger A, Zeschnigk M, et al. Methylation and silencing of the retinoic acid receptor-beta2 gene in breast cancer. J Natl Cancer Inst. 2000;92:826–32.
Muller HM, Fiegl H, Goebel G, Hubalek MM, Widschwendter A, Muller-Holzner E, et al. MeCP2 and MBD2 expression in human neoplastic and non-neoplastic breast tissue and its association with oestrogen receptor status. Br J Cancer. 2003;89:1934–9.
Teschendorff AE, Menon U, Gentry-Maharaj A, Ramus SJ, Gayther SA, Apostolidou S, et al. An epigenetic signature in peripheral blood predicts active ovarian cancer. PLoS One. 2009;4:e8274.
Teschendorff AE, Zhuang J, Widschwendter M. Independent surrogate variable analysis to deconvolve confounding factors in large-scale microarray profiling studies. Bioinformatics. 2011;27:1496–505.
Koch CM, Reck K, Shao K, Lin Q, Joussen S, Ziegler P, et al. Pluripotent stem cells escape from senescence-associated DNA methylation changes. Genome Res. 2013;23:248–59.
Gore ME. ICON2 trial. UK consensus statement on standard practice for chemotherapy in ovarian cancer. Lancet. 1999;353:587–8.
ICON2 Consortium. ICON2: randomised trial of single-agent carboplatin against three-drug combination of CAP (cyclophosphamide, doxorubicin, and cisplatin) in women with ovarian cancer. ICON Collaborators. International Collaborative Ovarian Neoplasm Study. Lancet. 1998;352:1571–6.
Ardizzoni A, Boni L, Tiseo M, Fossella FV, Schiller JH, Paesmans M, et al. Cisplatin- versus carboplatin-based chemotherapy in first-line treatment of advanced non-small-cell lung cancer: an individual patient data meta-analysis. J Natl Cancer Inst. 2007;99:847–57.
Burger H, Zoumaro-Djayoon A, Boersma AW, Helleman J, Berns EM, Mathijssen RH, et al. Differential transport of platinum compounds by the human organic cation transporter hOCT2 (hSLC22A2). Br J Pharmacol. 2010;159:898–908.
Jing L, Yuan W, Ruofan D, Jinjin Y, Haifeng Q. HOTAIR enhanced aggressive biological behaviors and induced radio-resistance via inhibiting p21 in cervical cancer. Tumour Biol. 2014;36:3611–9.
Arseneau J, Blessing JA, Stehman FB, McGehee R. A phase II study of carboplatin in advanced squamous cell carcinoma of the cervix (a Gynecologic Oncology Group Study). Invest New Drugs. 1986;4:187–91.
Weiss GR, Green S, Hannigan EV, Boutselis JG, Surwit EA, Wallace DL, et al. A phase II trial of carboplatin for recurrent or metastatic squamous carcinoma of the uterine cervix: a Southwest Oncology Group study. Gynecol Oncol. 1990;39:332–6.
McGuire III WP, Arseneau J, Blessing JA, DiSaia PJ, Hatch KD, Given Jr FT, et al. A randomized comparative trial of carboplatin and iproplatin in advanced squamous carcinoma of the uterine cervix: a Gynecologic Oncology Group study. J Clin Oncol. 1989;7:1462–8.
Qiu JJ, Wang Y, Ding JX, Jin HY, Yang G, Hua KQ. The long non-coding RNA HOTAIR promotes the proliferation of serous ovarian Cancer cells through the regulation of cell cycle arrest and apoptosis. Exp Cell Res. 2015;333:238–48.
The research leading to these results has received funding from the European Union’s Seventh Framework Programme (FP7/2007-2013) under grant agreement number 305428 (Project EpiFemCare) and was strongly supported by funds from The Eve Appeal (http://www.eveappeal.org.uk/) and a grant from the UCLH/UCL Comprehensive Biomedical Research Center and was undertaken at UCLH/UCL, which received a proportion of its funding from the Department of Health NIHR Biomedical Research Centers funding scheme.
A.E.T. was supported by a Heller Research Fellowship and by the Shanghai Institute for Biological Sciences/Chinese Academy of Sciences. R.G.J. was supported by a MRC Career Development Fellowship. H.B.S was supported by Helse Vest, Research Council of Norway and The Norwegian Cancer Society (Harald Andersens legacy). B.D. and C.G.T. were supported by the Inger and John Fredriksen Foundation for Ovarian Cancer Research. L.D. was supported by grant number RO1CA 133117 from the United States National Institutes of Health. H.M.H. is supported by a translational research fellowship from Dutch Cancer Society (KWF 2013-5869) This project was supported by the European Network Translational Research in Gynaecological Oncology (ENTRIGO) of the European Society of Gynaecological Oncology (ESGO).
No funding bodies had any role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Arturo Orjalo was an employee of Biosearch Technologies, Novato, CA, USA and has moved to Genentech, South San Francisco, CA, USA after submission of this manuscript. The other authors declare that they have no competing interests associated with the content of this paper. None of the authors have a private company associated with the content of this manuscript.
MW was responsible for the concept of the study and design of the experiments, obtained funding and drafted the manuscript. AET did the statistical analyses. HF, ELG, BLF, EMJJB, EW, HBS, GBAW, AGJvdZ, BD, CGT, SL and IV were responsible for data collection and provided materials and data for experiment conduction. AJ, S-HL, MK, WW, HF, KC, IE, LD, AO, HMH, LN, AK, WY, RGJ, IV, HC and IJJ contributed to molecular studies or in vitro experiments or the interpretation of these data. All authors contributed to data analysis/interpretation and manuscript preparation. All authors have approved the final version of this manuscript.
Additional description of methods and sample sets. (PDF 221 kb)
Clinicopathological characteristics of ovarian cancer patients (INNSBRUCK set) for which HOTAIR expression and DNAme have been done. (PDF 189 kb)
Normalised beta-value DNAme data matrix that was used to identify the 67 HOTAIR -associated CpGs. This data matrix has dimensions 5000 CpGs (the 5000 most variable positions) by 134 samples; samples are annotated according to whether they are HOTAIR + ve/-ve and type of chemotherapy received. (XLSX 9331 kb)
Clinicopathological features of patients (GRONINGEN set) stratified according to HOTAIR expression. (PDF 161 kb)
Clinicopathological characteristics of ovarian cancer patients from TCGA set. (PDF 178 kb)
Clinicopathological characteristics of ovarian cancer patients from the EUROPE set. (PDF 177 kb)
Ten-fold internal cross-validations to identify an optimal DNAme signature. Upper panel shows the total misclassification error (y-axis) as a function of the shrinkage threshold (x-axis) used. Lower panel shows the misclassification error for each phenotype (1 = low HOTAIR expression, 2 = high HOTAIR expression) as a function of the same shrinkage threshold. The optimal minimal classifier was found at a threshold of approximately 1.47, corresponding to a 67-CpG signature at an estimated false discovery rate (FDR) of approximately 0.17 (not shown). The FDR was estimated using a permutation scheme as implemented in the pamr R-package, and the relatively low FDR (only about 17 % of the 67 CpGs are expected to be false positives) demonstrates the presence of a genuine DNAme signal associated with HOTAIR expression. (PDF 13 kb)
Hazard ratios for relapse and death for patients who received carboplatin treatment (a) and patients who did not receive carboplatin treatment (b) in the original INNSBRUCK set. (PDF 183 kb)
Hazard ratios (HR) for death of patients who received carboplatin-based treatment in the GRONINGEN set and had RNA available to analyze HOTAIR. (PDF 138 kb)
Kaplan-Meier survival estimates in patients from the BERGEN set who received no chemotherapy (untreated group, n = 9) (a) or carboplatin-based chemotherapy (n = 40) (b) and stratified according to HOTAIR expression. Patients (n = 49) treated in Bergen (Norway) and whose cancers had >25 % stroma component were analyzed; 8, 2, 34 and 5 had stage 1, 2, 3 and 4 disease, respectively; 32, 6, 9 and 2 had a serous, mucinous, endometrioid and clear cell cancer, respectively. In the untreated group, significantly more patients had stage 1 disease (44 % versus 10 % in the carboplatin group) and no residual disease after primary surgery (75 % versus 38 % in the carboplatin group). The top tertile expressing samples were deemed as high (positive) HOTAIR expressors and compared with low/absent (negative) HOTAIR expressors. There is no significant difference between high and low HOTAIR expressors with regards to grade, stage or residual disease. Comparing HOTAIR-positive with HOTAIR-negative patients, the hazard ratio is 2.55 (95 % confidence interval 1.14–5.68), P value 0.018. (PDF 244 kb)
Boxplots of beta methylation values of the 67 CpGs (INNBSRUCK set) which demonstrate the largest difference between HOTAIR -negative (indicated as “ 1 ” and green boxes ) and HOTAIR -positive (indicated as “ 2 ” and red boxes ) ovarian cancer samples. (PDF 144 kb)
Enrichment odds ratios (OR) with 95 % confidence intervals for PCGTs according to different definitions among CpGs undergoing significant hyper- and hypomethylation with HOTAIR (HA) expression (top 500 CpGs) in ovarian cancer samples (INNSBRUCK set). PCGTs are defined as genes associated with SUZ12, EZH2 and H3K27me3 (TriplePCGT) or any one of these factors (SinglePCGT) in human embryonic stem cells (Lee et al. ) or as the ~850 genes that gain H3K27me3 upon overexpression of HOTAIR in MDAMB321 breast cancer cells as described by Gupta et al.  (PRC2-MDAMB231) or PRC2 targets in human embryonic fibroblast (PRC2-hEF) (Bracken et al. 2006). (PDF 5 kb)
Expression analyses for the GRONINGEN set. Of the 175 samples from the GRONINGEN set, 114 samples had matched array expression profiles (Operon Human v3 ~ 35 K 70-mer two-color oligonucleotide arrays, GSE13876, Crijns et al. [23 ]) and HOTAIR expression (36 HOTAIR -ve, 78 HOTAIR + ve). We correlated the expression of genes on the array to HOTAIR expression and selected genes with a t-statistic P value < 0.05 (1147 genes). From these we then selected those genes which are PRC2 targets in MDAMB231 breast cancer cell lines (PRC2-MDAMB231; Gupta et al. ) or in human embryonic fibroblasts (PRC2-hEF; Bracken et al. 2006). There were 42 and 97 such PRC2-MDAMB231 and PRC2-hEF genes, respectively. The average expression (AvExp) over these enriched PRC2-MDAMB231 and PRC2-hEF genes was then computed for each sample separately and these values are compared between HOTAIR-ve and HOTAIR + ve samples. Wilcoxon rank sum test P value is given. (PDF 5 kb)
HOTAIR DNAme score in ovarian cancer cell lines and MSCs depending on number of passages after stable HOTAIR transfection (ovarian cancer cells) or after starting in vitro culture (MSCs). The HOTAIR DNAme score is the Pearson correlation coefficient between the 67-CpG HOTAIR DNAme signature and the corresponding DNAme profile of each of the cell lines. Whereas neither HOTAIR expression nor passage number had an impact on the correlation coefficient in ovarian cancer cell lines, early passage (multipotent) MSCs (irrespective of whether they had been irradiated or not) showed a much stronger association with the 67-CpG HOTAIR signature than higher passage (senescent) MSCs (from the same individuals). (PDF 128 kb)
Performance of the DNAme signature in the primary INNSBRUCK carboplatin-treated ovarian cancer set and TCGA set. (PDF 215 kb)
Predicting HOTAIR expression with the 67-CpG DNA methylation signature. A ROC analysis was used to optimize the cutoff of the correlation score between the 67-CpG DNAme signature and the corresponding DNAme of each of the 63 ovarian cancer samples (INNSBRCUK set, carboplatin-treated subgroup). A correlation score of −0.16 ensures approximately 80 % sensitivity and 80 % specificity between the 67-CpG DNAme signature and HOTAIR RNA expression. (PDF 4 kb)
Comparison of the predictive score of our HOTAIR DNA methylation signature (red) to those obtained from 1000 randomized signatures (green) (randomized signatures have the same number of CpGs as the original signature). To further substantiate our results, we compared the predictive score of our HOTAIR DNAme signature to a randomized signature (same number of CpGs as in the original signature) obtained by randomly permuting the methylation profiles of CpGs in TCGA data. By comparing the observed Cox score to the ones obtained by a large number of such randomizations, we evaluated the predictive significance of the HOTAIR DNAme signature against the background probability. The randomization procedure showed that in only 2 of the 10,000 runs (P < 0.001) the Cox score was more significant than the observed (unpermutated) score. Thus, the selection of predictive CpGs in our primary ovarian cancer set identified CpGs more likely to be predictive of carboplatin response in the independent TCGA data. (PDF 3 kb)
Performance of the DNAme signature in the carboplatin (a) and cisplatin (b) treated EUROPE set. (PDF 212 kb)
Performance of the DNAme signature in the carboplatin-treated set from the ROCHESTER-MAYO set. (PDF 201 kb)
Kaplan-Meier survival estimates for patients whose tumours expressed HOTAIR (a, c) and whose tumours did not express HOTAIR (b, d) in the INNSBRUCK (a, b) and GRONINGEN sets (c, d). Survival analysis was performed according to the type of chemotherapy patients received: carboplatin monotherapy (Carbo-Mono), carboplatin-paclitaxel (Carbo-Pacli), carboplatin-cyclophosphamide (Carbo-Cyclo) or cisplatin-based chemotherapy (Cisplat). Hazard ratios (HR) and P values were calculated comparing cisplatin- with carboplatin-based chemotherapy regimens. (PDF 281 kb)
Kaplan-Meier survival estimates in patients in TCGA set whose tumours had a high (a) and a low (b) HOTAIR DNA methylation score. Survival analysis was performed according to the type of chemotherapy patients received: carboplatin based (Carboplatin), carboplatin-cisplatin (Carbo-Cisp) or cisplatin chemotherapy (Cisplatin). Hazard ratios (HR) and P values were calculated comparing cisplatin- with carboplatin-based (non-cisplatin-containing) chemotherapy regimens. (PDF 8 kb)
Chemosensitivity of A2780 and OVCAR8 ovarian cancer cells which are stably transfected with LacZ (control) or HOTAIR. Cells were treated with cisplatin (0.5–25 μM) or carboplatin (10–160 μM) for 3 days and analysed by the cell survival MTT assay (Sigma). (PDF 397 kb)
About this article
Cite this article
Teschendorff, A.E., Lee, SH., Jones, A. et al. HOTAIR and its surrogate DNA methylation signature indicate carboplatin resistance in ovarian cancer. Genome Med 7, 108 (2015). https://doi.org/10.1186/s13073-015-0233-4
- Ovarian Cancer
- Ovarian Cancer Tissue
- Primary Ovarian Cancer
- HOTAIR Expression