Suppression of hypoxia-induced stemness and chemoresistance in ovarian tumors

Suppression of hypoxia-induced stemness and chemoresistance in ovarian tumors | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Suppression of hypoxia-induced stemness and chemoresistance in ovarian tumors Sandra Munoz-Galvan, Eva M. Verdugo-Sivianes, José M. Santos-Pereira, and 2 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-3730407/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 26 Feb, 2024 Read the published version in Journal of Experimental & Clinical Cancer Research → Version 1 posted 5 You are reading this latest preprint version Abstract Background: Hypoxia in solid tumors is an important source of chemoresistance that can determine poor patient prognosis. Such chemoresistance relies on the presence of cancer stem cells (CSCs), and hypoxia promotes their generation through transcriptional activation by HIF transcription factors. Methods. We used OC cell lines, xenograft models, OC patient samples, transcriptional databases, iPSCs and ATAC-seq . Results Here, we show that hypoxia induces CSC formation and chemoresistance in ovarian cancer through transcriptional activation of the PLD2 gene. Mechanistically, HIF-1a activates PLD2 transcription through hypoxia response elements, and both hypoxia and PLD2 overexpression lead to increased accessibility around stemness genes, detected by ATAC-seq, at sites bound by AP-1 transcription factors. This in turn provokes a rewiring of stemness genes, including the overexpression of SOX2 , SOX9 or NOTCH1 . PLD2 overexpression also leads to decreased patient survival, enhanced tumor growth and CSC formation, and increased iPSCs reprograming, confirming its role in dedifferentiation to a stem-like phenotype. Importantly, hypoxia-induced stemness is dependent on PLD2 expression, demonstrating that PLD2 is a major determinant of de-differentiation of ovarian cancer cells to stem-like cells in hypoxic conditions. Finally, we demonstrate that high PLD2 expression increases chemoresistance to cisplatin and carboplatin treatments, both in vitro and in vivo , while its pharmacological inhibition restores sensitivity. Conclusions. Altogether, our work highlights the importance of the HIF-1a-PLD2 axis for CSC generation and chemoresistance in OC and proposes an alternative treatment for patients with high PLD2 expression. ovarian cancer hypoxia phospholipase D therapy resistance stemness Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 BACKGROUND Hypoxia is a common feature of the tumor microenvironment in solid tumors (Harris, 2002 ). Hypoxia is generated by insufficient oxygen diffusion towards parts of the tumor that are not irrigated, although the highly irregular tumor microvasculature may also generate hypoxic regions. Cancer cells under hypoxic conditions undergo a series of transcriptional changes that are induced by hypoxia-inducible factors (HIF), HIF-1, -2 and − 3, with HIF-1 being the best known. HIFs are heterodimeric transcription factors containing bHLH-PAS domains and consist of an α subunit regulated by oxygen levels and a constitutively expressed β subunit (also called ARNT). The HIFa-ARNT dimers bind DNA at specific sequences known as hypoxia-response elements (HREs) to promote target gene expression (Schito and Semenza, 2016 ). In normoxia, HIF-α subunits are hydroxylated at specific proline and asparagine residues by prolyl hydroxylase 2 (PHD2) and are recognized and targeted for degradation by the von Hippel‒Lindau (VHL) tumor suppressor. However, hypoxia leads to the inhibition of such hydroxylation and subsequent HIF-α accumulation, heterodimerization with ARNT and transcriptional activation of target genes (Kaelin and Ratcliffe, 2008 ). The multiple effects of hypoxia on the biology of tumors include preventing apoptosis and promoting proliferation and autophagy, inducing metabolic alterations, and promoting angiogenesis, the epithelial-to-mesenchymal (EMT) transition, invasion and metastasis (Jing et al., 2019 ; Wilson and Hay, 2011 ). Indeed, HIF-1α is overexpressed in multiple tumor types and is associated with a poor prognosis in patients and therapy resistance (Chen et al., 2018a ; Klemba et al., 2020 ; Lin and Koong, 2018 ; Mayer et al., 2008 ; Simiantonaki et al., 2008 ; Sun et al., 2007 ; Ueda et al., 2017 ; Unruh et al., 2003 ). Ovarian cancer (OC) is the most lethal gynecological malignancy (Siegel et al., 2020 ) mainly due to its nonspecific clinical manifestations, which lead to a late diagnosis and high chemoresistance (Torre et al., 2018 ). Similar to other solid tumors, hypoxia is a key modulator of the tumor microenvironment in OC, affecting not only the primary tumor but also the ascitic fluid, which is the main route of spread of these malignancies and shows low oxygen levels (Klemba et al., 2020 ; Munoz-Galvan and Carnero, 2020 ). Therefore, ovarian tumors are highly hypoxia-dependent and this dependency influence the response to treatment. A major cause of chemoresistance is the persistence of a cancer cell subpopulation known as cancer stem cells (CSCs) or tumor-initiating cells. CSCs are able to recapitulate a new tumor since they possess self-renewal and pluripotency properties similar to those of normal stem cells (Batlle and Clevers, 2017 ; Beck and Blanpain, 2013 ). CSCs are resistant to common antitumor therapies, which may indeed cause their enrichment, leading to chemoresistance and relapse (Carnero et al., 2016 ; Maugeri-Sacca et al., 2011 ). Several studies indicate that HIF-1 is required for maintaining CSCs and that its activation in hypoxia leads to the increased expression of stem marker genes in multiple cancer types (Mathieu et al., 2011 ; Peng and Liu, 2015 ). Therefore, hypoxia may promote the generation of CSCs leading to chemoresistance through HIF factors, but our knowledge of HIF targets that may be responsible for CSC activity is limited. We previously showed that the PLD2 gene, encoding phospholipase D2, is overexpressed in colorectal cancer and induces stemness in cancer cells through communication with the tumor microenvironment (Munoz-Galvan et al., 2019b ). Here we find that PLD2 is also overexpressed in ovarian cancer patients, being associated with poor patient survival. We explore a possible connection of PLD2 with hypoxia in OC and demonstrate that PLD2 expression in OC cells is stimulated by hypoxia and that HIF-1α promotes PLD2 transcription through HREs at its promoter and an intronic enhancer. This leads to increased chromatin accessibility around stemness genes, which are in turn overexpressed provoking an enhancement of tumor growth and formation of CSC-like tumorspheres. This is corroborated by induced pluripotent stem cell (iPSC) reprogramming experiments that confirm the role of hypoxia and PLD2 in dedifferentiation of OC cells to stem-like cells. Importantly, the increase in CSC-like features induced by hypoxia relies on PLD2 expression, indicating that the hypoxia-PLD2 axis is a major contributor to tumor stemness. We confirm these findings in transcriptional databases of OC patients, where high PLD2 expression is associated with a transcriptional rewiring of genes involved in the hypoxia response and in the maintenance of stem cells. Finally, we show that PLD2 overexpression causes resistance to platinum-based compounds and propose a new therapy based on pharmacological inhibition of phospholipases D to suppress such chemoresistance. RESULTS PLD2 is overexpressed in OC patients and in ovarian cancer cells under hypoxic conditions. First, we wondered whether PLD2 was overexpressed in OC patients. For this, we analysed PLD2 expression in 7 OC and one ovarian tissue databases using the R2 platform and found that PLD2 expression was significantly higher in tumoral than in non-tumoral samples (Fig. 1 A). This result was confirmed by a comparison between patients and control individuals in 3 databases containing their own controls (GSE18520, GSE4595 and GSE3866) (Fig. 1 B). Next, we wondered whether PLD2 expression was associated with patient survival and analysed 4 OC databases with available overall survival (OS) data (GSE13876, GSE19161, GSE23554 and GSE31245). We split patients in low-risk and high-risk groups based on their OS, but the difference in OS between both groups was statistically significant only in one database, GSE19161 (Fig. 1 C; Supplementary Figure S1 ). However, the expression levels of PLD2 were significantly higher in the high-risk group than in the low-risk group in the 4 databases (Fig. 1 C; Supplementary Figure S1 ), indicating that PLD2 is commonly overexpressed in OC patients and may be associated with decreased patient survival. Similar to other solid tumors, OC shows hypoxic areas, and its main dissemination site, the ascitic fluid, is also characterized by hypoxia (Klemba et al., 2020 ). Since PLD2 expression is increased under hypoxic conditions in colon cancer cells (Liu et al., 2020 ) and is involved in tumor stemness through communication with the microenvironment in colorectal cancer (Munoz-Galvan et al., 2019b ), we wondered whether hypoxia could increase the expression of PLD2 in OC. For this, we selected SKOV3, OVCAR8 and ES-2 OC cell lines for our studies and analysed the expression levels of PLD2 under hypoxic conditions. We found that PLD2 expression was similar among the three cell lines and that oxygen levels of 3% led to a 2-fold increase in the expression of PLD2 in all of them, while the well-known hypoxia target genes LDHA and VEGFA showed a similar increase, validating the results (Fig. 1 D). The hypoxia-induced increase in PLD2 expression was also observed at the protein level by performing immunofluorescence in the three OC cell lines (Fig. 1 E). In addition, we further demonstrated these findings by using the HIF-hydroxylase inhibitor DMOG, which increases HIF levels, generating a hypoxia-like phenotype under normoxic conditions (Supplementary Figure S2A-B). Thus, we treated OC cells with DMOG and found that this treatment resulted in a similar increase in hypoxia marker gene expression as that induced by hypoxia, also leading to the observed increase in PLD2 expression at both the mRNA and protein levels (Fig. 1 D-E). Therefore, we can conclude that PLD2 expression is promoted by hypoxia. HIF-1α activates PLD2 transcription through HREs at promoter and hypoxia-specific enhancer regions. We aimed to understand how hypoxia promotes PLD2 expression. HIF-1α is considered the master transcriptional regulator of the cellular response to hypoxia. It forms a heterodimer with ARNT that binds HREs to control the expression of hypoxia-response genes (Schito and Semenza, 2016 ). To address the possibility that HIF-1α could regulate PLD2 expression at the transcriptional level, we first searched for possible cis-regulatory elements (CREs) near the PLD2 gene that may contain HREs. According to public 3D chromatin conformation experiments (micro-C) in human embryonic stem cells, the PLD2 gene is located within a topologically associating domain (TAD) of 190 kb with a high interaction frequency at the 3D level and is relatively isolated from the neighbouring regions (Fig. 1 F). To identify CREs that may regulate PLD2 expression, we focused on a smaller region of 50 kb surrounding the PLD2 gene with a higher interaction frequency with the PLD2 promoter, which we called the PLD2 regulatory region. We scanned the CREs annotated by ENCODE within the PLD2 regulatory region for the presence of the DNA binding motif of HIF1A with a high score (> 90% relative score). We found 15 out of 35 CREs fulfilling this condition, some of which corresponded to gene promoters and others to enhancers, including the PLD2 promoter and an enhancer in PLD2 intron 12 (Fig. 1 F). These CREs with high-score HIF1A motifs represent putative HREs. To assess the regulatory activity of these two CREs in the PLD2 gene (promoter and putative enhancer) containing the HIF1A motif, we cloned both genomic regions in promoter and enhancer reporter vectors controlling the expression of the luciferase gene. We transfected OC cells with these vectors and measured the luciferase activity under normoxia and hypoxia. We found that the PLD2 promoter was able to activate reporter expression in normoxia and that hypoxia led to a significant increase in luciferase activity (8-fold) (Fig. 1 G), suggesting that the PLD2 promoter responds to hypoxic conditions by increasing PLD2 transcription. However, the enhancer contained within the PLD2 intron was unable to activate reporter expression in normoxia but led to a 22-fold increase in its expression in hypoxia (Fig. 1 G), indicating that this enhancer acts as an HRE in OC cells. Next, we wanted to validate functionally that the hypoxia-induced upregulation of PLD2 expression was indeed mediated by HIF-1α. Therefore, we depleted HIF1A using a small interfering RNA (siRNA) in the three OC cell lines and found that it suppressed the increase in the PLD2 protein levels induced by hypoxia (Fig. 1 H). Additionally, we generated hif1a mutant OC cell lines by transfecting cells with a mutant hif1a allele (Hu et al., 2007 ) that is unable to be hydroxylated and, therefore, is constitutively active even under normoxic conditions. Interestingly, we observed that the PLD2 levels in hif1a mutant cells were as high as those induced by hypoxia on a HIF1α wild-type background in both normoxia and hypoxia, confirming that HIF-1α activates PLD2 expression (Fig. 1 I-J). Altogether, these data indicate that hypoxia induces PLD2 expression in OC cells through transcriptional activation by HIF-1α at HREs in the PLD2 gene. Hypoxia alters the chromatin landscape of ovarian cancer cells in a PLD2-dependent manner. We wondered whether PLD2 expression mediated by HIF-1α could have an impact on hypoxia-induced gene regulation. For this, we first analyzed the effect of hypoxia and PLD2 expression in the epigenomic landscape of OC cells through ATAC-seq experiments in SKOV3 cells under normoxia and hypoxia and altered PLD2 expression under normoxia ( PLD2 overexpression) and hypoxia ( PLD2 depletion). We computationally called open chromatin regions (ATAC peaks) in normoxia and hypoxia and compared both conditions by a differential accessibility analysis; we detected 140 and 102 peaks with increased or decreased accessibility in hypoxia, respectively (Fig. 2 A). The heatmaps and aggregate profiles of these differentially accessible regions (DARs) showed that the peaks with increased accessibility in hypoxia were also more open upon PLD2 overexpression in normoxia, although to a lower extent, and vice versa, with the peaks showing decreased accessibility in hypoxia (Fig. 2 B), suggesting that PLD2 overexpression in normoxia has a similar effect on chromatin accessibility as hypoxia. Indeed, the DARs of EV- versus PLD2 -overexpressing cells under normoxia showed similar changes in accessibility under hypoxia, reinforcing the previous idea (Supplementary Figure S3A-B). Both DARs in normoxia versus hypoxia and control versus PLD2 overexpression were associated with genes enriched in Gene Ontology terms related to the response to hypoxia or well-known functions of PLD2 , respectively (Supplementary Figure S3C-D). Moreover, changes in accessibility induced by hypoxia were suppressed by PLD2 depletion (Fig. 2 B), suggesting that the effect of hypoxia in chromatin accessibility is mediated by PLD2. Similarly, changes in accessibility upon PLD2 overexpression in normoxia were absent in PLD2 -depleted cells in hypoxia (Supplementary Figure S3B), reinforcing the previous idea. Altogether, these data indicate that both hypoxia and PLD2 overexpression induce similar alterations in the chromatin landscape of OC cells and that the effect of hypoxia is mediated by PLD2. Next, we sought to investigate the possible mechanisms driving the alterations in the chromatin accessibility landscape induced by hypoxia and PLD2 overexpression. For this, we first performed motif enrichment analyses of DARs. We found that DARs in hypoxia were enriched in the motifs of the AP-1, ETS and C2H2 zinc finger transcription factor families in the increased accessibility sites and the C2H2 zinc finger and fork head families in the decreased accessibility sites (Supplementary Figure S4A). Similar enrichments were found in DARs upon PLD2 overexpression (Supplementary Figure S4B), reinforcing the idea of a similar effect under both conditions. To further elucidate the possible TFs involved in hypoxia and PLD2-mediated epigenomic changes, we estimated differential TF binding among the conditions based on footprints in ATAC-seq data. Using this approach, we confirmed the increased chromatin binding of AP-1 family transcription factors, such as FOS and JUN, consistent with their implication in the response to hypoxia (Kunz and Ibrahim, 2003 ). (Fig. 2 C-D). We also found other TF families showing increased TF binding in hypoxia, such as homeobox, paired box or fork head TFs, and the C2H2 zinc finger TF ZBTB32, while TF families with decreased binding in hypoxia included members of the bHLH, NF-Y and ETS families, among others (Supplementary Figure S4C). When we compared these TFs with those showing increased binding upon PLD2 overexpression in normoxia, we found that most of these TFs overlapped (Supplementary Figure S4C). In particular, 25 TF motifs of the AP-1 family showed increased binding under both conditions. Interestingly, these AP-1 and 20 more TF motifs with increased binding in hypoxia showed decreased binding upon PLD2 depletion (Supplementary Figure S4C-E), indicating that they are dependent on PLD2 expression. Altogether, these results suggest a function of the AP-1 family of TFs mediating the alterations in the chromatin landscape mediated by hypoxia and PLD2. Since we previously connected PLD2 overexpression with tumor stemness in colorectal cancer, we wondered whether alterations in the chromatin accessibility landscape of OC cells induced by hypoxia and PLD2 could result in increased expression of stemness genes. Therefore, we first selected ATAC peaks falling within the putative regulatory landscapes of the genes associated with stem cell maintenance and proliferation. Clustering of these 3,572 peaks revealed 4 groups with different accessibility levels and behaviours (Supplementary Figure S5A). Among them, Cluster 3 corresponded with peaks with increased accessibility in both hypoxic and PLD2 -overexpressing cells but decreased accessibility in PLD2 -depleted cells. Among the stemness genes associated with this cluster, we found SOX9 , PROM1 , WNT7A or JAG1 (Fig. 2 E; Supplementary Figure S5B). These results indicate that hypoxia promotes chromatin accessibility around stemness genes in a PLD2-dependent manner in OC cells, and suggests that hypoxia and PLD2 could be connected with tumor stemness in OC. High PLD2 expression in OC patients leads to the transcriptomic rewiring of stemness and hypoxia genes To determine whether PLD2 expression in OC patients is related to hypoxia and stemness, we analysed the expression of genes related to these functions in the three OC databases with available expression data from control individuals. First, we selected the genes annotated to the Gene Ontology (GO) term “Response to Hypoxia” and whose expression was significantly correlated with that of PLD2 in OC patients (p 0.2 or <-0.2). Then, we performed hierarchical clustering of patient and control individuals based on the expression levels of these genes (Fig. 3 A; Supplementary Figure S6). In the GSE18520 database, the clustering clearly separated the control individuals (‘non-tumoral’, NT) and a reduced group of patients that we termed ‘Tumoral Cluster 1’ (T1) from most patients who clustered in what we termed ‘Tumoral Cluster 2’ (T2) (Fig. 3 A). The patients at T1 showed a transcriptional profile of hypoxia-related genes more similar to NT and clearly different from T2. In contrast, the clustering in the GSE4095 and GSE38666 databases clearly separated clusters of NT and tumoral (T) individuals, who showed distinct transcriptional profiles of hypoxia-related genes (Supplementary Figure S6). These results suggest that there is a switch in the expression of hypoxia-related genes in OC tumors compared with healthy ovaries. Next, we repeated the hierarchical clustering with the genes annotated to the GO term “Stem cell maintenance” and whose expression was significantly correlated with that of PLD2 in OC patients (p 0.2 or <-0.2). Surprisingly, this clustering based on stem-related genes separated the patients and control individuals into the same clusters as the hypoxia-related genes (Fig. 3 A; Supplementary Figure S6), suggesting a connection between both groups of genes that supports the model of CSC generation induced by hypoxia. Then, we plotted the expression levels of PLD2 in the 3 clusters obtained from the GSE18520 database and found that PLD2 exhibited significantly increased expression in Cluster T2 compared with that in Cluster NT, while Cluster T1 showed similar levels to NT (Fig. 3 B). This observation suggests that a connection exists among PLD2 expression, the hypoxia response and stemness since PLD2 expression is misregulated only in patients showing transcriptional profiles highly different from healthy controls. In addition, we checked the expression of stem-associated and hypoxia-related genes in these clusters. As shown in Fig. 3 B, the expression of PAX8 , a well-known OC marker, was significantly increased in both patient Clusters T1 and T2, similar to other stemness genes, such as SOX9 , SOX17 , EPCAM , PROM1 , CD24 , NOTCH1 and WNT7A . Other genes in this group showed a significant increase in expression only in Cluster T2, coinciding with higher PLD2 levels, including PAX2 , POU5F1 ( OCT4 ), SOX5 , SOX11 , CD34 and TP63 , while the others were unaffected or even significantly reduced, such as KLF4 , NANOG or SOX2 , although the latter showed a nonsignificant increase in Cluster T2 (Fig. 3 B). This finding suggests that there is an OC stemness signature in patients that may be stronger with higher PLD2 expression. In addition, some hypoxia-related genes showed a significant increase in expression in both Clusters T1 and T2, including VEGFA , SLC2A1 ( GLUT1 ), HK2 and NOX4 , or only in T2, including SLC2A4 , NOS1 and MMP14 , and a nonsignificant increase in SLC2A14 ( GLUT14 ) (Fig. 3 B). Similar results were obtained in the NT and T clusters in the GSE4095 and GSE38666 databases (Supplementary Figure S6). Altogether, these results suggest that there is transcriptional rewiring of the expression of hypoxia- and stem-related genes in OC patients with PLD2 overexpression. PLD2 promotes tumorigenesis and CSC-like features in ovarian cancer cells. The ability of hypoxia to induce a CSC-like phenotype in OC cells has been previously observed in several cancer types (Liang et al., 2012 ; Mathieu et al., 2011 ; Peng and Liu, 2015 ). We first aimed to validate these results in our OC cell lines SKOV3, OVCAR8 and ES-2 and found that hypoxia led to significant increases in the number of tumorspheres, which were generated by growing the cells under low-attachment conditions, and in the percentage of holoclones, both of which were used as a proxy for CSCs (Supplementary Figure S7A-B). Next, we analysed the expression of stem cell markers in OC cell lines grown under hypoxic conditions and detected an increase in the mRNA levels of NANOG , CD44 , SOX2 and EPCAM and the percentage of cells containing the surface CSC marker CD133 (Supplementary Figure S7C-D). These results confirm that hypoxia induces a CSC-like phenotype in OC cells. Next, we wondered whether increased PLD2 expression led to an increase in the CSC population in OC cells, as suggested by the chromatin accessibility and gene expression data. Therefore, we first established OC cell lines expressing ectopic PLD2 cDNA or depleted of PLD2 using a short hairpin RNA (shRNA). The expression of PLD2 under these conditions was assessed at the mRNA and protein levels (Fig. 4 A-B and Supplementary Figure S8). We observed that the overexpression of PLD2 in OC cells led to a significant increase in the number of clones generated by the three cell lines, while a significant decrease was detected in the OVCAR8 cells upon PLD2 depletion (Fig. 4 C), suggesting that PLD2 promotes tumor growth. To address this question, we analysed the growth of these cell lines and found that the enhanced PLD2 expression led to a significant increase in proliferation, while the PLD2 depletion generated the opposite effect with statistical significance in all cell lines (Fig. 4 D). This effect was further confirmed in vivo by generating xenograft models of OC cells overexpressing or depleted of PLD2 , showing an increase or decrease in the tumor volume, respectively, 50 days after transplantation (Fig. 4 E). Altogether, these results indicate that PLD2 expression promotes tumorigenesis. Next, we wondered whether PLD2 expression was related to the formation of ovarian CSCs. First, we analysed the formation of different types of colonies, including holoclones, meroclones and paraclones, which are considered stem cells, transit-amplifying cells and differentiated cells, respectively (Barrandon and Green, 1987 ). We found a significant increase in the percentage of holoclones and a significant decrease in the percentage of paraclones in the three cell lines overexpressing PLD2 (Fig. 4 F). A significant decrease in holoclone formation was also observed in the three OC cell lines upon PLD2 depletion. Furthermore, we measured the formation of tumorspheres under low attachment conditions in OC cells overexpressing or depleted of PLD2 . We found that PLD2 overexpression led to a significant increase in the number of tumorspheres, while PLD2 depletion generated the opposite effect (Fig. 4 G), although we did not observe changes in the size of such tumorspheres. These data indicate that PLD2 expression promotes the formation of ovarian CSCs. Finally, we analysed the expression of pluripotency and CSC marker genes in our OC cell lines overexpressing or depleted of PLD2 . We found that the expression of SOX2 , CD44 and EPCAM was significantly increased in the ES-2 and SKOV3 cells overexpressing PLD2 , while NANOG was only increased in ES-2 cells (Fig. 4 H). In OVCAR8 cells, the expression of these genes was increased in the same trend, although in a nonsignificant manner. Then, we measured the expression levels of these pluripotency genes in the tumorspheres extracts. First, we found that PLD2 was highly expressed in the tumorspheres compared with that in the total extracts transfected with only the empty vector (Fig. 4 H), confirming that CSCs indeed have higher expression levels of PLD2 . The expression of stemness genes was also increased in the tumorspheres compared with that in the total extracts, as expected, while they were further upregulated in most cases upon PLD2 overexpression, or downregulated upon PLD2 depletion (Fig. 4 H). Altogether, these data indicate that PLD2 , whose expression is induced by hypoxia, is an important oncogene in OC and that its overexpression leads to increased tumor stemness. The hypoxia-induced stemness of ovarian cancer cells partially depends on PLD2 expression. Thus far, we have showed that both hypoxia and PLD2 expression led to an increase in CSCs in OC cells (Fig. 4 , Supplementary Figure S7) and that PLD2 expression was increased under hypoxic conditions in a HIF-1α-dependent manner (Fig. 1 ). Thus, we wondered whether both phenomena were connected and whether the hypoxia-induced increase in CSC-like features was dependent on PLD2 expression. To address this, we first analysed the expression levels of stemness genes by RT‒qPCR using custom TaqMan Array plates containing probes against a selection of these genes in OC cells. We observed that either hypoxia or PLD2 overexpression in normoxia led to the increased expression of many stemness genes, while PLD2 depletion largely suppressed this increase (Fig. 5A). Indeed, hierarchical clustering of the four analysed conditions showed that the samples corresponding to EV hypoxic cells and PLD2 -overexpressing cells clustered together, while EV normoxic cells and PLD2 -depleted cells clustered separately (Fig. 5B). We confirmed these result by RT‒qPCR of individual representative genes, including SOX2 , NANOG , CD44 and EPCAM , showing that either hypoxia or PLD2 overexpression in normoxia lead to increased expression of stemness genes, while combination of both conditions further increased their expression (Fig. 5C). PLD2 depletion partially suppressed the hypoxia-induced enhancement of expression, with a lower non-statistically significant effect in normoxia, and rescue experiments confirmed the specificity of PLD2 depletion (Fig. 5C). Next, we analyzed the formation of tumorspheres in normoxia and hypoxia with altered PLD2 expression. We observed that either PLD2 overexpression or hypoxia led to a similar increase in the formation of tumorspheres, with only a slightly higher nonsignificant increase when both conditions were combined (Fig. 5D). However, PLD2 depletion caused a partial suppression of the increase in tumorsphere formation under hypoxic conditions in SKOV3 and OVCAR8 cells, which was rescued by overexpressing back PLD2 in shRNA-transfected cells (Fig. 5D and Supplementary Figure S9A), suggesting that PLD2 is partially responsible for the hypoxia-induced stemness. Then, we measured the formation of holoclones, meroclones and paraclones under hypoxic conditions upon PLD2 overexpression or depletion and found that the increase in the percentage of holoclones induced by hypoxia was further enhanced by the PLD2 overexpression, while it was suppressed by PLD2 depletion and rescued back by expressing PLD2 after its depletion (Supplementary Figure S9B). Then, we measured the protein levels of the pluripotency factors Sox2, Sox17, Sox9 and Notch1 (found to correlate with PLD2 in OC patients, Fig. 3 B) by immunofluorescence in tumorspheres to determine whether PLD2 could influence their expression in CSCs. We validated PLD2 protein levels in tumorspheres (Supplementary Figure S9C) and found that either hypoxia or PLD2 overexpression led to an increase in the levels of Sox2, Sox9 and Notch1, while only hypoxia led to an increase in Sox17 protein levels (Fig. 5E). In addition, PLD2 depletion led to a partial suppression of the hypoxia-induced expression of Sox2, Sox9 and Notch1 that was rescued by expressing back PLD2 in these cells (Fig. 5E), suggesting that PLD2 plays a role in the generation of CSCs in hypoxia through these genes. These observations were confirmed at the mRNA level by RT‒qPCR (Fig. 5F). Altogether, these data indicate that PLD2 plays a major role in the induction of the CSC phenotype in hypoxia, promoting the expression of specific stem-related genes, such as SOX2 , SOX9 or NOTCH1 . Finally, we extended the gene expression analyses to EMT genes using TaqMan Arrays to assess whether PLD2 expression may have a role in tumor invasion and metastasis. We found that either hypoxia or PLD2 overexpression in normoxia led to an increase in the expression of many of these genes, but PLD2 depletion was unable to suppress such increase (Supplementary Figure S9A). Indeed, hierarchical clustering of the samples did not exhibit the pattern observed in stemness genes (Supplementary Figure S9B), and results were further validated by RT‒qPCR of particular EMT genes (Supplementary Figure S9C). Consistently, invasiveness assays using Boyden’s chamber showed that both PLD2 overexpression and hypoxia were able to increase invasion, but PLD2 depletion did not have any effect (Supplementary Figure S9D). These results indicate that while the increased expression of stemness genes induced by hypoxia relies on PLD2 overexpression, this is not the case for EMT genes and suggests that PLD2 is a specific mediator of the increase in CSCs induced by hypoxia in OC cells. Hypoxia-mediated reprogramming to induced pluripotent stem cells is dependent on PLD2. We aimed to obtain additional evidence of the contribution of PLD2 to dedifferentiation or reprogramming events mediated by hypoxia that may generate ovarian CSCs from normal OC cells. Therefore, we performed reprogramming experiments of mouse embryonic fibroblasts (MEFs) to induced pluripotent stem cells (iPSCs) in normoxia and hypoxia and upon alteration of PLD2 expression levels (overexpression or depletion). We used a previously published protocol (Yoshida et al., 2009 ) in which MEFs were infected using a HEK293T cell-derived virus that provides OSKM genes and Nanog reporter retroviruses and then cocultured on SNL feeder cells that produce LIF. Then, the samples were incubated with or without hypoxia for 7 days, and the efficiency of iPSC generation was measured for additional 5 days. Cell reprogramming and the acquisition of pluripotency were assessed by colony morphology, alkaline phosphatase and nanog promoter-driven GFP expression analyses (Fig. 6 A) to assess the effect of PLD2 and hypoxia on the efficiency of the reprogramming process and the acquisition of stem cell-like properties. Using this protocol, we found that, as expected, hypoxia led to a significant increase in the generation of iPSCs (Fig. 6 B). Furthermore, we found that PLD2 overexpression in normoxia provoked a similar increase in iPSC generation, consistent with its effect on the generation of CSCs, and that the combination of both hypoxia and PLD2 overexpression further increased iPSC formation (Fig. 6 B). This confirms that high PLD2 expression leads to dedifferentiation processes. Finally, PLD2 depletion in hypoxia suppressed the increased iPSC production induced by hypoxia, and this was recovered by expressing back PLD2 in PLD2 -depleted cells (Fig. 6 B), suggesting that PLD2 is an important mediator in the activation of pluripotency by hypoxic conditions. Overexpression of PLD2 leads to chemotherapy resistance in ovarian tumors. Since we showed that PLD2 overexpression leads to an increase in CSC-like cells in OC and CSCs were previously proposed to be responsible for chemotherapy resistance and tumor relapse, we wondered whether PLD2 overexpression could cause resistance to conventional therapy in ovarian tumors. Therefore, we first analysed the expression levels of PLD2 in our own cohort of OC patients. The immunohistochemistry analyses showed that PLD2 protein levels were higher in tumors than in healthy tissue (Fig. 7 A), and RT‒qPCR revealed that PLD2 mRNA was significantly more abundant in OC patients than in control non-tumoral samples (Fig. 7 B), thus confirming the results observed in the transcriptomic databases (Fig. 1 A-B). Then, we separated our patient samples into those who were sensitive or resistant to platinum-based chemotherapy (without or with tumor relapse within the next 6 months after chemotherapy, respectively) and analysed PLD2 expression levels. Importantly, we found that the resistant patients showed significantly higher expression of PLD2 than the sensitive patients (Fig. 7 C), suggesting that PLD2 overexpression may contribute to resistance to platinum-based therapy. Resistant patients in our cohort showed reduced OS and PFS (Fig. 7 D-E), consistent with the reduced survival of patients with high PLD2 expression (Fig. 7 C). Next, we analysed the effect of PLD2 expression on resistance to platinum compounds in OC cells in vitro . First, cells overexpressing or depleted of PLD2 were treated with increasing concentrations of cisplatin and carboplatin, and the IC 50 values were calculated in each case. We found that PLD2 overexpression led to a significant increase in the IC 50 values, while PLD2 depletion led to only a weak nonsignificant reduction (Fig. 7 F). This finding suggests that higher PLD2 expression causes resistance to platinum-based compounds. We repeated these experiments under hypoxic conditions and found that hypoxia also led to increased IC 50 values (Fig. 7 F). However, PLD2 depletion reduced the hypoxia-induced increase in IC 50 values, suggesting that enhanced resistance to cisplatin and carboplatin in OC cells under hypoxic conditions relies on PLD2 expression. Finally, we performed in vivo analyses to validate our findings by establishing xenograft models from SKOV3 and OVCAR8 cells overexpressing or depleted of PLD2 and analysing tumor growth upon treatment with cisplatin. The control tumors from cells transfected with the empty vector were sensitive to the cisplatin treatment, significantly reducing tumor growth in xenografts from both SKOV3 and OVCAR8 cells (Fig. 7 H). However, tumors overexpressing PLD2 showed higher tumor growth that was not reduced upon cisplatin treatment, indicating both a higher aggressiveness of PLD2 -overexpressing tumors and the resistance of these tumors to cisplatin. However, PLD2 depletion resulted in significantly reduced tumor growth that was further reduced upon cisplatin treatment (Fig. 7 H). Importantly, cisplatin treatment in control tumors led to increased survival, while PLD2 -overexpressing tumors did not exhibit improved survival, and PLD2 -depleted tumors exhibited increased survival independent of cisplatin treatment (Fig. 7 I). These results indicate that the overexpression of PLD2 causes resistance to platinum-based chemotherapy in OC tumors. Combination treatment with cisplatin and a PLD inhibitor suppresses chemotherapy resistance in ovarian cancer. Finally, we wondered whether the increased therapy resistance to platinum-based compounds induced by PLD2 overexpression and hypoxia could be suppressed by the pharmacological inhibition of PLD2. For this, we used the PLD inhibitor (PLDi) 5-Fluoro-2-indolyl des-chlorohalopemide (FIPI), which inhibits the catalytic activity of phospholipases D (Ganesan et al., 2015 ). First, we tested this possibility in vitro by calculating the IC 50 in OC cells treated with cisplatin, PLDi and their combination in normoxia and hypoxia with altered levels of PLD2 . We found that the higher IC 50 to cisplatin in OC cells overexpressing PLD2 or in hypoxia was suppressed by the PLDi (Fig. 7 G). Next, we validated these results in vivo by establishing xenografts of OC cells expressing EV or PLD2 and treating mice with cisplatin, PLDi or their combination. We found that the increased tumor growth provoked by PLD2 overexpression was reduced upon treatment with PLDi (Fig. 7 H). Moreover, although treatment with cisplatin did not reduce the higher tumor growth induced by PLD2 overexpression, its combination with PLDi led to a significant reduction in tumor growth that was stronger than that following treatment with PLDi alone (Fig. 7 H). This finding was confirmed in xenografts from two OC cell lines (OVCAR8 and SKOV3) and led to an increase in survival (Fig. 7 I). Altogether, these results indicate that chemotherapy resistance to cisplatin caused by PLD2 overexpression can be overcome by the pharmacological inhibition of PLD2, suggesting that combined treatment with cisplatin and PLDi is a promising alternative treatment for patients with high PLD2 expression levels. DISCUSSION We show here that the HIF-1α-PLD2 axis is a major player in the chemoresistance of OC by promoting the generation of CSCs under hypoxic conditions. On the one hand, PLD2 expression is regulated by HIF-1α through HREs in its promoter and a hypoxia-specific enhancer; on the other hand, PLD2 expression is required for the full transcriptional and epigenomic rewiring promoted by hypoxia as evidenced using gene expression and chromatin accessibility analyses. In particular, the expression of stem-related genes, the opening of enhancers in the vicinity of these genes, the generation of CSCs induced by hypoxia, and the reprogramming of normal cells to iPSCs rely on normal PLD2 expression levels. These findings indicate that PLD2 is a major player in the response to hypoxia in cancer cells that leads to increased stem cell properties resulting in higher chemoresistance. Hypoxia is a feature of the tumor microenvironment in regions with low oxygen supply that is known to increase the stemness features of cancer cells in several types of cancer (Carnero and Lleonart, 2016 ; Li et al., 2009b; Mathieu et al., 2011 ; Peng and Liu, 2015 ; Schwab et al., 2012 ; Wang et al., 2011 ), including OC in which hypoxia has been shown to increase the stem-like properties of cancer cells (Liang et al., 2012 ). HIF-ARNT can regulate the expression of many genes that promote the hypoxic response (Kaelin and Ratcliffe, 2008 ; Schito and Semenza, 2016 ). Therefore, HIF factors may exert their effect of increasing CSCs via multiple mechanisms. In OC, HIF factors contribute to the upregulation of pluripotency factor genes, such as SOX2 or OCT3/4 (Liang et al., 2012 ; Seo et al., 2016 ), proliferation pathways, such as Notch or Wnt (Chau et al., 2013 ; Seo et al., 2016 ), or epigenetic modulation by affecting chromatin modifiers, such as SIRT1 (Qin et al., 2017 ). Here we show that the OC cell lines ES-2, SKOV3 and OVCAR8 under hypoxia show upregulated expression of PLD2 , which encodes phospholipase D2, in these cells in a Hif-1α-dependent manner (Fig. 1 ), as recently reported in colon cancer (Liu et al., 2020 ). Using DNA binding motif analyses and reporter assays, we also demonstrate that PLD2 expression is regulated at the transcriptional level by HREs located within the PLD2 promoter and a hypoxia-specific enhancer that activates PLD2 transcription in hypoxia, thus providing a mechanistic explanation of Hif-1α-mediated PLD2 overexpression. Therefore, it is likely that the HIf-1α-PLD2 axis works in other solid tumors under hypoxic conditions. The expression of PLD2 is elevated in several cancer types (Carnero et al., 1994 ; Frankel et al., 1999 ; Henkels et al., 2013 ; Munoz-Galvan et al., 2019b ; Saito et al., 2007 ; Song et al., 1991 ; Zheng et al., 2006 ), and here, using public patient databases and our cohort of patients, we show that this is also the case in OC in which it may be related to decreased OS (Fig. 1 ). Furthermore, consistent with the results in OC cell lines, clustering analyses of OC patient gene expression revealed that PLD2 expression is correlated with the rewiring of transcriptomic programs of the response to hypoxia and stem cell maintenance (Fig. 3 ). Indeed, we found two different clusters of patients, one of which showed gene expression patterns more different from healthy controls coinciding with higher PLD2 expression. Although the expression of the OC marker PAX8 and other stemness- and hypoxia-related genes was enhanced in both clusters, other markers showed increased expression only in the patient cluster with high PLD2 expression, suggesting that PLD2 may promote stemness through specific genes or pathways as we show in SOX2 , SOX9 and NOTCH1 , but not SOX17 . These data indicate that a correlation exists between PLD2 expression and highly altered transcriptomic programs of the response to hypoxia and stemness. In addition, we provide evidence that PLD2 is important for the alteration in the epigenomic landscape provoked by hypoxia since both hypoxia and PLD2 overexpression in normoxia lead to similar alterations in chromatin accessibility that are counteracted by PLD2 depletion, including the opening of enhancers in the proximity of genes related to the stem fate, which we show were upregulated (Fig. 2 ). These changes connect hypoxia with stemness gene expression and likely occur through the activation of CREs bound by the AP-1 family of TFs. In a previous study analysing OC tumors, solid metastasis and effusions, higher PLD2 expression was described in effusions rather than solid tumors and metastasis (Harel-Dassa et al., 2017 ). This finding is consistent with our findings since OC effusions, most of which are peritoneal, are characterized by low oxygen levels and a high content in CSCs (Munoz-Galvan and Carnero, 2020 ). Therefore, the hypoxia-PLD2 axis seems to play a major role in the generation of ovarian CSCs. Importantly, PLD2 depletion partially suppresses the effect induced by hypoxia (Fig. 5), suggesting that PLD2 is an important mediator of hypoxia-induced stemness. This finding was also corroborated by reprogramming experiments of MEFs to iPSCs in which PLD2 was required for the increased reprogramming induced by hypoxia. Indeed, PLD2 might also potentiate the response to hypoxia by a positive feedback loop. This hypothesis can be supported since the combination of hypoxia and PLD2 overexpression additively enhanced the stem properties of OC cells. The activation of HIF-1α expression or activity by PLD2 has been reported in endothelial, glioma and renal cancer cells (Ghim et al., 2014 ; Han et al., 2014 ; Toschi et al., 2008 ), but there is evidence of the opposite effect in HEK293 cells (Park et al., 2015 ), suggesting that this feedback loop may work in specific contexts or cell types. Nevertheless, whether the effect of PLD2 on stemness is an autocrine or a paracrine effect remains to be elucidated. In this regard, we previously showed that exosomes from PLD2 -overexpressing colorectal cancer cells induced senescence in stromal fibroblasts (Munoz-Galvan et al., 2019b ), which, in turn, induced WNT pathway activation and increased stemness in tumor cells. In OC, hypoxia-induced exosomes have been involved in increased tumorigenic properties and chemoresistance by several mechanisms. These include exosome-containing oncogenic proteins, such as STAT3 and FAS (Dorayappan et al., 2018 ), microRNAs that altered tumor-associated macrophages (Chen et al., 2018b ; Zhu et al., 2019 ), as well as plasma gelsolin, which induces the conversion of chemosensitive OC cells to chemoresistant cells (Asare-Werehene et al., 2020 ). Finally, a role of PLD1 and PLD2 inducing exosome secretion in OC cells has also been recently proposed (Onallah et al., 2022 ), suggesting that PLD2 may also influence the tumor microenvironment in OC. Several mechanisms have been described to promote chemotherapy resistance under hypoxia in OC, including the upregulation of the ABCG2 transporter gene, which increases drug efflux (He et al., 2019 ; Wang et al., 2016 ), c- KIT overexpression (Chau et al., 2013 ) and high cysteine levels (Nunes et al., 2018 ). CSCs are responsible for chemotherapy resistance (Batlle and Clevers, 2017 ; Beck and Blanpain, 2013 ), and ovarian CSCs were identified sixteen years ago and reported to be chemoresistant (Bapat et al., 2005 ; Hu et al., 2010 ). In agreement with this idea, we previously found several markers linking ovarian CSCs and chemoresistance (Munoz-Galvan et al., 2019a ; Munoz-Galvan et al., 2020 ). Here, using OC cells and xenograft models, we show that the overexpression of PLD2 leads to resistance to platinum-derived compounds, including cisplatin and carboplatin (Fig. 7 ). Moreover, PLD2 expression is higher in patients resistant to platinum-based chemotherapy than in sensitive patients, confirming our results in cell lines and mouse models. How PLD2 provokes such resistance is an intriguing issue, although it is likely that its enzymatic product PA, an important molecule acting as a second messenger in multiple cellular functions (Jang et al., 2012 ), might play some relevant role in avoiding chemotherapy-induced cell damage. Indeed, we previously showed that PA administration has similar effects as PLD2 overexpression (Munoz-Galvan et al., 2019b ). This is particularly relevant in ovarian tumors, which are typically treated with platinum-based compounds and show high rates of chemoresistance, with frequent metastasis in hypoxic environments, such as abdominal ascites. CONCLUSIONS Our findings suggest a model in which hypoxia leads to the transcriptional overexpression of PLD2 in OC, which, in turn, generates PA and induces the generation of chemoresistant ovarian CSCs. Therefore, we propose an alternative treatment based on a combination of cisplatin and the pharmacological inhibition of PLD2. Our in vitro and in vivo results demonstrate that this combined treatment may be useful for patients with high PLD2 expression, who are resistant to conventional therapy with cisplatin alone. METHODS Cell culture. Cells were cultured according to the manufacturer's recommended procedures in McCoy (ES-2 line) or RPMI (SKOV3 and OVCAR8 lines) and incubated at 37°C in 5% CO 2 in a humidified atmosphere. Gene transfer. The gene transfer was performed as previously described (Munoz-Galvan et al., 2020 ). The PLD2 overexpression plasmid was described in (Munoz-Galvan et al., 2019b ). The shRNA against PLD2 was provided by Origene. Proliferation assay. The proliferation assay was performed as previously described (Munoz-Galvan et al., 2020 ). Cytotoxic assay. ES-2, SKOV3 or OVCAR8 cells were seeded and then treated with platinum drugs and/or the PLD inhibitor 5-Fluoro-2-indolyl des-chlorohalopemide (FIPI) at 300nM concentration 24 hours later. After 96 hours, cells were stained with 0.5% crystal violet. Then, the crystal violet was solubilized in 20% acetic acid and quantified at 595 nm absorbance to measure the cell viability. Maintenance of mouse colonies. All experiments involving animals received expressed approval from the IBIS/HUVR Ethical Committee for the Care and Health of Animals. The mice were maintained in the IBIS animal facility according to the facility guidelines, which are based on the Real Decreto 53/2013 and were sacrificed by CO 2 inhalation either using a planned procedure or as a human endpoint when the animals showed significant signs of illness. Colony formation assay and clonal heterogeneity analysis This analysis was performed as previously described in (Munoz-Galvan et al., 2020 ). Briefly, in total, 10 3 cells were seeded onto 10 cm plates, and every condition was evaluated in triplicate. The medium was replaced every 3 days for 12 days, and the colonies were fixed, stained and counted. The values are expressed as the number of observed colonies among the 10 3 seeded cells. To analyze the clonal heterogeneity, 10 2 random colonies were classified in triplicate as having the following phenotypes: holoclone, meroclone and paraclone (Li et al., 2008 ), which are considered stem cells, transit-amplifying cells and differentiated cells, respectively (Barrandon and Green, 1987 ). Sphere-forming assay . In total, 2x 10 3 cells were resuspended in 1 ml of complete MammoCultTM Basal Medium (Stemcell Tech) and seeded in ultralow attachment plates. The cultures were imaged, the tumorspheres were counted, and their diameters were quantified using CellSenseDimension software on Days 2, 3 and 4. Western blot analyses and immunofluorescence were performed according to standard procedures. Information of antibodies and dilutions is shown in Supplementary Table S1 . RT–qPCR. The total RNA was isolated using an RNeasy kit (Qiagen), and cDNA was generated from 1 µg of RNA with MultiScribe Reverse Transcriptase (Applied Biosystems). qPCR was performed using a TaqMan Assay (Applied Biosystems) with probes. The relative mRNA expression was calculated as 2 −∆Ct relative to the ACTB gene. Information of probes is shown in Supplementary Table S1 . Fluorescence-activated cell sorting For FACS staining, live cells were incubated with antibodies for 30 minutes at dilutions specified in the manufacturer's protocols. iPSCs protocol Briefly, mouse embryonic fibroblasts (MEFs) were infected with 4 different retroviral vectors encoding the Yamanaka factors (pMXs-Oct3/4, pMXs-Sox2, pMXs-Klf4, and pMXs-cMyc), an additional lentiviral vector that expresses GFP in cells where the Nanog promoter/enhancer is active (mNanog-pGreenZeo), and the corresponding plasmid overexpressing PLD2, carrying shPLD2 or carrying Ev. Seventy-two hours after the infection, MEFs were seeded on top of an SNL feeder layer in the presence of ES media, and the media was renewed every 24 h. Cell reprogramming and the acquisition of pluripotency were assessed by colony morphology, GFP expression and alkaline phosphatase activity assays to assess the effect of the gene of interest/condition of interest* on the efficiency of the reprogramming process and the acquisition of stem cell-like properties. ATAC-seq. ATAC-seq assays were performed using standard protocols (Buenrostro et al., 2013 ; Fernandez-Minan et al., 2016 ), with minor modifications. Briefly, 70,000 ovarian cancer cells overexpressing PLD2 carrying shPLD2 or carrying Ev growing under normoxic or hypoxic conditions were collected by centrifugation for 5 min at 500 g 4°C. The supernatant was removed, and the cells were washed with PBS. Then, the cells were lysed in 50 µl of lysis buffer (10 mM Tris-HCl pH 7.4, 10 mM NaCl, 3 mM MgCl 2 , 0.1% NP-40, 1x Roche Complete protease inhibitors cocktail) by pipetting up and down. The whole cell lysate was used for TAGmentation, which was centrifuged for 10 min at 500 g 4°C, resuspended in 50 µl of the Transposition Reaction containing 2.5 µl of Tn5 enzyme and TAGmentation Buffer (10 mM Tris-HCl pH 8.0, 5 mM MgCl2, 10% w/v dimethylformamide), and incubated for 30 min at 37°C. Immediately after TAGmentation, DNA was purified using a Minelute PCR Purification Kit (Qiagen) and eluted in 20 µl. Libraries were generated by PCR amplification using NEBNext High-Fidelity 2X PCR Master Mix (NEB). The resulting libraries were multiplexed and sequenced in a HiSeq 4000 paired-end lane, producing 100 M 49-bp paired-end reads per sample. Quantification and statistical analysis. All statistical analyses were performed using GraphPad Prism 4. The distribution of the quantitative variables among different study groups was assessed using parametric (Student’s t test ) or nonparametric (Kruskal–Wallis or Mann–Whitney) tests as appropriate. The experiments were performed a minimum of three times and were performed in independent triplicates each time. The survival data from the patient databases were analyzed by a log-rank Mantel‒Cox statistical test. Analyses of cancer patient databases. We performed meta-analyses using the R2 Genomics analysis and visualization platform ( http://hgserver1.amc.nl ) to analyze the PLD2 expression levels in tumor and nontumor ovarian samples from the databases. The statistical significance of the tumor versus normal samples was assessed (P < 0.05). Patient survival was analyzed using an R2 Genomics analysis and visualization platform ( http://hgserver1.amc.nl ), which was developed by the Department of Oncogenomics of the Academic Medical Center (Amsterdam, Netherlands). Kaplan‒Meier plots showing patient survival were generated using the databases with available survival data with the scan method, which searches for the optimum survival cut-off based on statistical analyses (log-rank test), thereby identifying the most significant expression cut-off. ATAC-seq data analyses. ATAC-seq reads were aligned to the GRCh38 (hg38) human genome assembly using Bowtie2 2.3.5 (Langmead and Salzberg, 2012 ), and pairs separated by more than 2 kb were removed. For ATAC-seq, the Tn5 cutting site was determined as position − 4 (minus strand) or + 5 (plus strand) from each read start, and this position was extended 5 bp in both directions. Reads below 150 bp were considered nucleosome free. The conversion of the SAM alignment files to BAM was performed using SAMtools 1.9 (Li et al., 2009a). The conversion of BAM to BED files and peak analyses, such as overlaps or merges, were carried out using the Bedtools 2.29.2 suite (Quinlan and Hall, 2010 )⁠. The conversion of BED to BigWig files was performed using the genomecov tool from Bedtools and the wigToBigWig utility from UCSC (Haeussler et al., 2019 ). For ATAC-seq, peaks were called using the MACS2 2.1.1.20160309 algorithm (Zhang et al., 2008 ) with an FDR < 0.05 for each replicate and merged into a single pool of peaks that was used to calculate the differentially accessible sites with the DESeq2 1.18.1 package in R 3.4.3 (Love et al., 2014 ); a value < 0.05 was set as the cut-off for statistical significance of the differential accessibility. For visualization purposes, reads were extended 100 bp for ATAC-seq. For the data comparison, all ATAC-seq experiments used were normalized using reads falling into peaks to counteract differences in background levels between experiments and replicates (Santos-Pereira et al., 2019 ). Heatmaps, average profiles and k-means clustering of the ATAC-seq data were generated using the computeMatrix and plotHeatmap tools from the Deeptools 3.5 toolkit (Ramirez et al., 2016 ). The TF motif enrichment was calculated using XSTREME (Grant and Bailey, 2021 ) with the standard parameters. For the gene assignment to ATAC peaks, we used the GREAT 3.0.0 tool (Hiller et al., 2013 )⁠, with the basal plus extension association rule and the standard parameters (5 kb upstream, 1 kb downstream, and 1 Mb maximum extension). For the footprinting analyses, we used TOBIAS 0.12.9 (Bentsen et al., 2020 ). First, we performed bias correction using ATACorrect and calculated the footprint scores with ScoreBigwig, both with the standard parameters. Then, we used BINDetect to determine the differential TF binding of all vertebrate motifs in the JASPAR database (Fornes et al., 2020 ). We considered motifs with a linear fold-change ≥ 15% between conditions differentially bound. Aggregated ATAC-seq signals in the footprints were visualized using PlotAggregate. ATAC-seq data generated in this study is available through the Gene Expression Omnibus (GEO) accession number GSE210599. In vivo xenograft studies. Tumor growth was assayed following the subcutaneous injection of 4x10 6 SKOV3 or OVCAR8 cells that were transfected with a plasmid carrying PLD2 or shRNA against PLD2 in cohorts of five nude mice each that were analyzed weekly. The tumors were measured using callipers. All mice were sacrificed once the growth experiment was completed. In vivo xenograft treatment. Tumors were harvested when they reached 1500 mm 3 , cut into 2×2×2 mm pieces and reimplanted. Mice were randomly allocated to the drug-treated and control-treated (solvent only) groups, and once the tumor reached 20 mm 3 , the mice received the appropriate treatment for 4 weeks (2 doses/week). The mice were monitored daily for signs of distress and weighed twice a week. The tumor size was measured, and the size was estimated according to the following equation: tumor volume = [length x width 2 ]/2. The experiments were terminated when the tumor reached 350 mm 3 or when the clinical endpoint was reached. The drugs cisplatin and carboplatin were obtained from Pharmacy HUVR and were freshly prepared and administered by intraperitoneal injection. We used higher doses in mice, assuming a 70 kg average weight for humans (125 mg/dose in humans) (Munoz-Galvan et al., 2020 ). We administered two doses per week as follows: 3.5 mg/kg of cisplatin (equivalent to 7 mg/kg, averaging 25 g body weights of each mouse) with or without 3 mg/kg of FIPI. We did not observe signs of toxicity. In vivo xenografts from tumorspheres. This assay involved the subcutaneous injection of 1×10 3 cells grown as tumorspheres into the hind legs of 4-week-old female athymic nude mice. The animals were treated as previously described, examined twice a week, incubated for 4 more weeks, and killed, and the tumors were extracted. The tumors were measured using callipers. Patient cohort. The entire procedure was approved by the local ethical committee of the HUVR (CEEA O309-N-15). A cohort of paraffin-embedded tissue samples from 25 patients with ovarian cancer was obtained from the biobank of Hospital Universitario Virgen del Rocío-Instituto de Biomedicina de Sevilla (Sevilla, Spain) for the RNA expression studies and the evaluation of the correlation of the clinicopathological features (see Supplementary Table S2). The samples were obtained from biopsies of patients who were subjected to platinum treatment and who were evaluated for their response according to the RECIST criteria; normal tissue, platinum-resistant tumor samples and platinum-sensitive tumor samples were obtained. The tumor samples were sent to the pathology laboratory for diagnosis and were prepared for storage with formalin fixation and paraffin embedding. The samples were stained with haematoxylin/eosin, and RNA was extracted from the tumor tissue. Declarations Ethics approval and consent to participate: All methods were performed in accordance with the relevant guidelines and regulations of the Institute for Biomedical Research of Seville (IBIS) and University Hospital Virgen del Rocio (HUVR). All animal experiments and the entire procedures of the patient cohort were performed according to the experimental protocol approved by HUVR Animals Ethics (CEI 0309-N-15). Consent for publication: Written consent for publication was obtained from all patients involved in our study Availability of data and materials: The datasets used and/or analysed during the current study are available from the corresponding author upon reasonable request. ATAC-seq data generated in this study is available through the Gene Expression Omnibus (GEO) accession number GSE210599. Competing interests: The authors declare that they have no other competing financial interests. Funding: This research was funded by Grants RTI2018-097455-B-I00 and PID2021-122629OB-I00 funded by MCIN/AEI/10.13039/501100011033 and by “ERDF A way of making Europe”, by the “European Union”. Additional grants from CIBER de Cáncer (CB16/12/00275), from Consejeria de Salud (PI-0397-2017) and Project P18-RT-2501 from 2018 competitive research projects call within the scope of PAIDI 2020—80% co-financed by the European Regional Development Fund (ERDF) from the Regional Ministry of Economic Transformation, Industry, Knowledge and Universities. Junta de Andalucía. Special thanks to the AECC (Spanish Association of Cancer Research) Founding Ref. GC16173720CARR for supporting this work. SMG was funded by a grant from the Fundación AECC. EMV-S and JMS-P are funded by postdoctoral fellowships from Junta de Andalucía (DOC_01655 and DOC_00512, respectively). Acknowledgements : The authors thank the donors and the HUVR-IBiS Biobank (Andalusian Public Health System Biobank and ISCIII-Red de Biobancos PT17/0015/0041) for the human specimens that were used in this study. Author contributions : SMG and AC conceived and designed this study. SMG and EMVS performed the experiments; JMSP analysed the NGS data; PEG collected the clinical data; and SMG and AC analysed and interpreted the data and drafted and edited the manuscript. All authors revised the manuscript References Asare-Werehene, M., K. Nakka, A. Reunov, C.T. Chiu, W.T. Lee, M.R. Abedini, P.W. Wang, D.B. Shieh, F.J. Dilworth, E. Carmona, T. Le, A.M. Mes-Masson, D. Burger, and B.K. Tsang. 2020. The exosome-mediated autocrine and paracrine actions of plasma gelsolin in ovarian cancer chemoresistance. Oncogene 39:1600-1616. Bapat, S.A., A.M. Mali, C.B. Koppikar, and N.K. Kurrey. 2005. Stem and progenitor-like cells contribute to the aggressive behavior of human epithelial ovarian cancer. Cancer Res 65:3025-3029. Barrandon, Y., and H. Green. 1987. Three clonal types of keratinocyte with different capacities for multiplication. Proc Natl Acad Sci U S A 84:2302-2306. Batlle, E., and H. Clevers. 2017. Cancer stem cells revisited. Nat Med 23:1124-1134. Beck, B., and C. Blanpain. 2013. Unravelling cancer stem cell potential. Nat Rev Cancer 13:727-738. Bentsen, M., P. Goymann, H. Schultheis, K. Klee, A. Petrova, R. Wiegandt, A. Fust, J. Preussner, C. Kuenne, T. Braun, J. Kim, and M. Looso. 2020. ATAC-seq footprinting unravels kinetics of transcription factor binding during zygotic genome activation. Nat Commun 11:4267. Buenrostro, J.D., P.G. Giresi, L.C. Zaba, H.Y. Chang, and W.J. Greenleaf. 2013. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat Methods 10:1213-1218. Carnero, A., A. Cuadrado, L. del Peso, and J.C. Lacal. 1994. Activation of type D phospholipase by serum stimulation and ras-induced transformation in NIH3T3 cells. Oncogene 9:1387-1395. Carnero, A., Y. Garcia-Mayea, C. Mir, J. Lorente, I.T. Rubio, and L.L. ME. 2016. The cancer stem-cell signaling network and resistance to therapy. Cancer Treat Rev 49:25-36. Carnero, A., and M. Lleonart. 2016. The hypoxic microenvironment: A determinant of cancer stem cell evolution. Bioessays 38 Suppl 1:S65-74. Chau, W.K., C.K. Ip, A.S. Mak, H.C. Lai, and A.S. Wong. 2013. c-Kit mediates chemoresistance and tumor-initiating capacity of ovarian cancer cells through activation of Wnt/beta-catenin-ATP-binding cassette G2 signaling. Oncogene 32:2767-2781. Chen, A., J. Sceneay, N. Godde, T. Kinwel, S. Ham, E.W. Thompson, P.O. Humbert, and A. Moller. 2018a. Intermittent hypoxia induces a metastatic phenotype in breast cancer. Oncogene 37:4214-4225. Chen, X., J. Zhou, X. Li, X. Wang, Y. Lin, and X. Wang. 2018b. Exosomes derived from hypoxic epithelial ovarian cancer cells deliver microRNAs to macrophages and elicit a tumor-promoted phenotype. Cancer Lett 435:80-91. Dorayappan, K.D.P., R. Wanner, J.J. Wallbillich, U. Saini, R. Zingarelli, A.A. Suarez, D.E. Cohn, and K. Selvendiran. 2018. Hypoxia-induced exosomes contribute to a more aggressive and chemoresistant ovarian cancer phenotype: a novel mechanism linking STAT3/Rab proteins. Oncogene 37:3806-3821. Fernandez-Minan, A., J. Bessa, J.J. Tena, and J.L. Gomez-Skarmeta. 2016. Assay for transposase-accessible chromatin and circularized chromosome conformation capture, two methods to explore the regulatory landscapes of genes in zebrafish. Methods Cell Biol 135:413-430. Fornes, O., J.A. Castro-Mondragon, A. Khan, R. van der Lee, X. Zhang, P.A. Richmond, B.P. Modi, S. Correard, M. Gheorghe, D. Baranasic, W. Santana-Garcia, G. Tan, J. Cheneby, B. Ballester, F. Parcy, A. Sandelin, B. Lenhard, W.W. Wasserman, and A. Mathelier. 2020. JASPAR 2020: update of the open-access database of transcription factor binding profiles. Nucleic Acids Res 48:D87-D92. Frankel, P., M. Ramos, J. Flom, S. Bychenok, T. Joseph, E. Kerkhoff, U.R. Rapp, L.A. Feig, and D.A. Foster. 1999. Ral and Rho-dependent activation of phospholipase D in v-Raf-transformed cells. Biochem Biophys Res Commun 255:502-507. Ganesan, R., M. Mahankali, G. Alter, and J. Gomez-Cambronero. 2015. Two sites of action for PLD2 inhibitors: The enzyme catalytic center and an allosteric, phosphoinositide biding pocket. Biochim Biophys Acta 1851:261-272. Ghim, J., J.S. Moon, C.S. Lee, J. Lee, P. Song, A. Lee, J.H. Jang, D. Kim, J.H. Yoon, Y.J. Koh, C. Chelakkot, B.J. Kang, J.M. Kim, K.L. Kim, Y.R. Yang, Y. Kim, S.H. Kim, D. Hwang, P.G. Suh, G.Y. Koh, Y.Y. Kong, and S.H. Ryu. 2014. Endothelial deletion of phospholipase D2 reduces hypoxic response and pathological angiogenesis. Arterioscler Thromb Vasc Biol 34:1697-1703. Grant, C.E., and T.L. Bailey. 2021. XSTREME: Comprehensive motif analysis of biological sequence datasets. bioRxiv Haeussler, M., A.S. Zweig, C. Tyner, M.L. Speir, K.R. Rosenbloom, B.J. Raney, C.M. Lee, B.T. Lee, A.S. Hinrichs, J.N. Gonzalez, D. Gibson, M. Diekhans, H. Clawson, J. Casper, G.P. Barber, D. Haussler, R.M. Kuhn, and W.J. Kent. 2019. The UCSC Genome Browser database: 2019 update. Nucleic Acids Res 47:D853-D858. Han, S., J. Huh, W. Kim, S. Jeong, S. Min do, and Y. Jung. 2014. Phospholipase D activates HIF-1-VEGF pathway via phosphatidic acid. Exp Mol Med 46:e126. Harel-Dassa, K., S. Yedgar, C.G. Trope, B. Davidson, and R. Reich. 2017. Phospholipase D messenger RNA expression and clinical role in high-grade serous carcinoma. Hum Pathol 62:115-121. Harris, A.L. 2002. Hypoxia--a key regulatory factor in tumour growth. Nat Rev Cancer 2:38-47. He, M., H. Wu, Q. Jiang, Y. Liu, L. Han, Y. Yan, B. Wei, F. Liu, X. Deng, H. Chen, L. Zhao, M. Wang, X. Wu, W. Yao, H. Zhao, J. Chen, and M. Wei. 2019. Hypoxia-inducible factor-2alpha directly promotes BCRP expression and mediates the resistance of ovarian cancer stem cells to adriamycin. Mol Oncol 13:403-421. Henkels, K.M., G.P. Boivin, E.S. Dudley, S.J. Berberich, and J. Gomez-Cambronero. 2013. Phospholipase D (PLD) drives cell invasion, tumor growth and metastasis in a human breast cancer xenograph model. Oncogene 32:5551-5562. Hiller, M., S. Agarwal, J.H. Notwell, R. Parikh, H. Guturu, A.M. Wenger, and G. Bejerano. 2013. Computational methods to detect conserved non-genic elements in phylogenetically isolated genomes: application to zebrafish. Nucleic Acids Res 41:e151. Hu, C.J., A. Sataur, L. Wang, H. Chen, and M.C. Simon. 2007. The N-terminal transactivation domain confers target gene specificity of hypoxia-inducible factors HIF-1alpha and HIF-2alpha. Mol Biol Cell 18:4528-4542. Hu, L., C. McArthur, and R.B. Jaffe. 2010. Ovarian cancer stem-like side-population cells are tumourigenic and chemoresistant. Br J Cancer 102:1276-1283. Jang, J.H., C.S. Lee, D. Hwang, and S.H. Ryu. 2012. Understanding of the roles of phospholipase D and phosphatidic acid through their binding partners. Prog Lipid Res 51:71-81. Jing, X., F. Yang, C. Shao, K. Wei, M. Xie, H. Shen, and Y. Shu. 2019. Role of hypoxia in cancer therapy by regulating the tumor microenvironment. Mol Cancer 18:157. Kaelin, W.G., Jr., and P.J. Ratcliffe. 2008. Oxygen sensing by metazoans: the central role of the HIF hydroxylase pathway. Mol Cell 30:393-402. Klemba, A., L. Bodnar, H. Was, K.K. Brodaczewska, G. Wcislo, C.A. Szczylik, and C. Kieda. 2020. Hypoxia-Mediated Decrease of Ovarian Cancer Cells Reaction to Treatment: Significance for Chemo- and Immunotherapies. Int J Mol Sci 21: Kunz, M., and S.M. Ibrahim. 2003. Molecular responses to hypoxia in tumor cells. Mol Cancer 2:23. Langmead, B., and S.L. Salzberg. 2012. Fast gapped-read alignment with Bowtie 2. Nat Methods 9:357-359. Li, H., X. Chen, T. Calhoun-Davis, K. Claypool, and D.G. Tang. 2008. PC3 human prostate carcinoma cell holoclones contain self-renewing tumor-initiating cells. Cancer Res 68:1820-1825. Li, H., B. Handsaker, A. Wysoker, T. Fennell, J. Ruan, N. Homer, G. Marth, G. Abecasis, R. Durbin, and S. Genome Project Data Processing. 2009a. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25:2078-2079. Li, Z., S. Bao, Q. Wu, H. Wang, C. Eyler, S. Sathornsumetee, Q. Shi, Y. Cao, J. Lathia, R.E. McLendon, A.B. Hjelmeland, and J.N. Rich. 2009b. Hypoxia-inducible factors regulate tumorigenic capacity of glioma stem cells. Cancer Cell 15:501-513. Liang, D., Y. Ma, J. Liu, C.G. Trope, R. Holm, J.M. Nesland, and Z. Suo. 2012. The hypoxic microenvironment upgrades stem-like properties of ovarian cancer cells. BMC Cancer 12:201. Lin, S.H., and A.C. Koong. 2018. Breathing New Life Into Hypoxia-Targeted Therapies for Non-Small Cell Lung Cancer. J Natl Cancer Inst 110: Liu, M., K. Du, B. Jiang, and X. Wu. 2020. High Expression of PhospholipaseD2 Induced by Hypoxia Promotes Proliferation of Colon Cancer Cells through Activating NF- kappa Bp65 Signaling Pathway. Pathol Oncol Res 26:281-290. Love, M.I., W. Huber, and S. Anders. 2014. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 15:550. Mathieu, J., Z. Zhang, W. Zhou, A.J. Wang, J.M. Heddleston, C.M. Pinna, A. Hubaud, B. Stadler, M. Choi, M. Bar, M. Tewari, A. Liu, R. Vessella, R. Rostomily, D. Born, M. Horwitz, C. Ware, C.A. Blau, M.A. Cleary, J.N. Rich, and H. Ruohola-Baker. 2011. HIF induces human embryonic stem cell markers in cancer cells. Cancer Res 71:4640-4652. Maugeri-Sacca, M., P. Vigneri, and R. De Maria. 2011. Cancer stem cells and chemosensitivity. Clin Cancer Res 17:4942-4947. Mayer, A., M. Hockel, A. Wree, C. Leo, L.C. Horn, and P. Vaupel. 2008. Lack of hypoxic response in uterine leiomyomas despite severe tissue hypoxia. Cancer Res 68:4719-4726. Munoz-Galvan, S., and A. Carnero. 2020. Targeting Cancer Stem Cells to Overcome Therapy Resistance in Ovarian Cancer. Cells 9: Munoz-Galvan, S., B. Felipe-Abrio, M. Garcia-Carrasco, J. Dominguez-Pinol, E. Suarez-Martinez, E.M. Verdugo-Sivianes, A. Espinosa-Sanchez, L.E. Navas, D. Otero-Albiol, J.J. Marin, M.P. Jimenez-Garcia, J.M. Garcia-Heredia, A.G. Quiroga, P. Estevez-Garcia, and A. Carnero. 2019a. New markers for human ovarian cancer that link platinum resistance to the cancer stem cell phenotype and define new therapeutic combinations and diagnostic tools. J Exp Clin Cancer Res 38:234. Munoz-Galvan, S., B. Felipe-Abrio, E.M. Verdugo-Sivianes, M. Perez, M.P. Jimenez-Garcia, E. Suarez-Martinez, P. Estevez-Garcia, and A. Carnero. 2020. Downregulation of MYPT1 increases tumor resistance in ovarian cancer by targeting the Hippo pathway and increasing the stemness. Mol Cancer 19:7. Munoz-Galvan, S., A. Lucena-Cacace, M. Perez, D. Otero-Albiol, J. Gomez-Cambronero, and A. Carnero. 2019b. Tumor cell-secreted PLD increases tumor stemness by senescence-mediated communication with microenvironment. Oncogene 38:1309-1323. Nunes, S.C., C. Ramos, F. Lopes-Coelho, C.O. Sequeira, F. Silva, S. Gouveia-Fernandes, A. Rodrigues, A. Guimaraes, M. Silveira, S. Abreu, V.E. Santo, C. Brito, A. Felix, S.A. Pereira, and J. Serpa. 2018. Cysteine allows ovarian cancer cells to adapt to hypoxia and to escape from carboplatin cytotoxicity. Sci Rep 8:9513. Onallah, H., S.T. Mannully, B. Davidson, and R. Reich. 2022. Exosome Secretion and Epithelial-Mesenchymal Transition in Ovarian Cancer Are Regulated by Phospholipase D. Int J Mol Sci 23: Park, M.H., S.S. Bae, K.Y. Choi, and S. Min do. 2015. Phospholipase D2 promotes degradation of hypoxia-inducible factor-1alpha independent of lipase activity. Exp Mol Med 47:e196. Peng, G., and Y. Liu. 2015. Hypoxia-inducible factors in cancer stem cells and inflammation. Trends Pharmacol Sci 36:374-383. Qin, J., Y. Liu, Y. Lu, M. Liu, M. Li, J. Li, and L. Wu. 2017. Hypoxia-inducible factor 1 alpha promotes cancer stem cells-like properties in human ovarian cancer cells by upregulating SIRT1 expression. Sci Rep 7:10592. Quinlan, A.R., and I.M. Hall. 2010. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26:841-842. Ramirez, F., D.P. Ryan, B. Gruning, V. Bhardwaj, F. Kilpert, A.S. Richter, S. Heyne, F. Dundar, and T. Manke. 2016. deepTools2: a next generation web server for deep-sequencing data analysis. Nucleic Acids Res 44:W160-165. Saito, M., M. Iwadate, M. Higashimoto, K. Ono, Y. Takebayashi, and S. Takenoshita. 2007. Expression of phospholipase D2 in human colorectal carcinoma. Oncol Rep 18:1329-1334. Santos-Pereira, J.M., L. Gallardo-Fuentes, A. Neto, R.D. Acemel, and J.J. Tena. 2019. Pioneer and repressive functions of p63 during zebrafish embryonic ectoderm specification. Nat Commun 10:3049. Schito, L., and G.L. Semenza. 2016. Hypoxia-Inducible Factors: Master Regulators of Cancer Progression. Trends Cancer 2:758-770. Schwab, L.P., D.L. Peacock, D. Majumdar, J.F. Ingels, L.C. Jensen, K.D. Smith, R.C. Cushing, and T.N. Seagroves. 2012. Hypoxia-inducible factor 1alpha promotes primary tumor growth and tumor-initiating cell activity in breast cancer. Breast Cancer Res 14:R6. Seo, E.J., D.K. Kim, I.H. Jang, E.J. Choi, S.H. Shin, S.I. Lee, S.M. Kwon, K.H. Kim, D.S. Suh, and J.H. Kim. 2016. Hypoxia-NOTCH1-SOX2 signaling is important for maintaining cancer stem cells in ovarian cancer. Oncotarget 7:55624-55638. Siegel, R.L., K.D. Miller, and A. Jemal. 2020. Cancer statistics, 2020. CA Cancer J Clin 70:7-30. Simiantonaki, N., M. Taxeidis, C. Jayasinghe, U. Kurzik-Dumke, and C.J. Kirkpatrick. 2008. Hypoxia-inducible factor 1 alpha expression increases during colorectal carcinogenesis and tumor progression. BMC Cancer 8:320. Song, J.G., L.M. Pfeffer, and D.A. Foster. 1991. v-Src increases diacylglycerol levels via a type D phospholipase-mediated hydrolysis of phosphatidylcholine. Mol Cell Biol 11:4903-4908. Sun, H.C., Z.J. Qiu, J. Liu, J. Sun, T. Jiang, K.J. Huang, M. Yao, and C. Huang. 2007. Expression of hypoxia-inducible factor-1 alpha and associated proteins in pancreatic ductal adenocarcinoma and their impact on prognosis. Int J Oncol 30:1359-1367. Torre, L.A., B. Trabert, C.E. DeSantis, K.D. Miller, G. Samimi, C.D. Runowicz, M.M. Gaudet, A. Jemal, and R.L. Siegel. 2018. Ovarian cancer statistics, 2018. CA Cancer J Clin 68:284-296. Toschi, A., J. Edelstein, P. Rockwell, M. Ohh, and D.A. Foster. 2008. HIF alpha expression in VHL-deficient renal cancer cells is dependent on phospholipase D. Oncogene 27:2746-2753. Ueda, S., T. Saeki, A. Osaki, T. Yamane, and I. Kuji. 2017. Bevacizumab Induces Acute Hypoxia and Cancer Progression in Patients with Refractory Breast Cancer: Multimodal Functional Imaging and Multiplex Cytokine Analysis. Clin Cancer Res 23:5769-5778. Unruh, A., A. Ressel, H.G. Mohamed, R.S. Johnson, R. Nadrowitz, E. Richter, D.M. Katschinski, and R.H. Wenger. 2003. The hypoxia-inducible factor-1 alpha is a negative factor for tumor therapy. Oncogene 22:3213-3220. Wang, W.J., H. Sui, C. Qi, Q. Li, J. Zhang, S.F. Wu, M.Z. Mei, Y.Y. Lu, Y.T. Wan, H. Chang, and P.T. Guo. 2016. Ursolic acid inhibits proliferation and reverses drug resistance of ovarian cancer stem cells by downregulating ABCG2 through suppressing the expression of hypoxia-inducible factor-1alpha in vitro. Oncol Rep 36:428-440. Wang, Y., Y. Liu, S.N. Malek, P. Zheng, and Y. Liu. 2011. Targeting HIF1alpha eliminates cancer stem cells in hematological malignancies. Cell Stem Cell 8:399-411. Wilson, W.R., and M.P. Hay. 2011. Targeting hypoxia in cancer therapy. Nat Rev Cancer 11:393-410. Yoshida, Y., K. Takahashi, K. Okita, T. Ichisaka, and S. Yamanaka. 2009. Hypoxia enhances the generation of induced pluripotent stem cells. Cell Stem Cell 5:237-241. Zhang, Y., T. Liu, C.A. Meyer, J. Eeckhoute, D.S. Johnson, B.E. Bernstein, C. Nusbaum, R.M. Myers, M. Brown, W. Li, and X.S. Liu. 2008. Model-based analysis of ChIP-Seq (MACS). Genome Biol 9:R137. Zheng, Y., V. Rodrik, A. Toschi, M. Shi, L. Hui, Y. Shen, and D.A. Foster. 2006. Phospholipase D couples survival and migration signals in stress response of human cancer cells. J Biol Chem 281:15862-15868. Zhu, X., H. Shen, X. Yin, M. Yang, H. Wei, Q. Chen, F. Feng, Y. Liu, W. Xu, and Y. Li. 2019. Macrophages derived exosomes deliver miR-223 to epithelial ovarian cancer cells to elicit a chemoresistant phenotype. J Exp Clin Cancer Res 38:81. Supplementary Files MunozGalvansupplementary.docx Cite Share Download PDF Status: Published Journal Publication published 26 Feb, 2024 Read the published version in Journal of Experimental & Clinical Cancer Research → Version 1 posted Editorial decision: Major revision 21 Dec, 2023 Reviewers agreed at journal 12 Dec, 2023 Reviewers invited by journal 12 Dec, 2023 Editor assigned by journal 11 Dec, 2023 First submitted to journal 11 Dec, 2023 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-3730407","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":259151158,"identity":"1b9c113d-548f-4008-ae2d-fc72710ca8f2","order_by":0,"name":"Sandra Munoz-Galvan","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAAvUlEQVRIiWNgGAWjYBADHgNmBsYHJGthNiDNGqByNgmiVMpPO2P24MMfBhlzdvZn1Tx/GOTkGwgZfjvH3HAGDwOPZTOP2W3eNgZjgwOEtEjnmEnzSPznMTjMw3abt4EhcQNBh80GaTEAev8w+7NioMPq5xNyGMNtkJYEkBYGM2YeNoYEBoIOu51WJjnjANgvxpJz2yQMNxDSIj87eZsEMMTszfmPP/zw5o+NPMEQQwfERc0oGAWjYBSMAgIAAChIMibuYuqHAAAAAElFTkSuQmCC","orcid":"","institution":"Science Research Council Inistitute of Biomedical Science in Seville: Instituto de Biomedicina de Sevilla","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"Sandra","middleName":"","lastName":"Munoz-Galvan","suffix":""},{"id":259151159,"identity":"c915c699-2911-478f-8514-fb562b2e000a","order_by":1,"name":"Eva M. Verdugo-Sivianes","email":"","orcid":"","institution":"Instituto de Biomedicina de Sevilla","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Eva","middleName":"M.","lastName":"Verdugo-Sivianes","suffix":""},{"id":259151160,"identity":"ffd87239-4467-41ac-9a07-3ae99783d5a7","order_by":2,"name":"José M. Santos-Pereira","email":"","orcid":"","institution":"CABD: Centro Andaluz de Biologia del Desarrollo","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"José","middleName":"M.","lastName":"Santos-Pereira","suffix":""},{"id":259151161,"identity":"2cfa1ff9-5af2-49a9-8869-b03a6fec654e","order_by":3,"name":"Purificación Estevez-García","email":"","orcid":"","institution":"Hospital Universitario Virgen del Rocio","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Purificación","middleName":"","lastName":"Estevez-García","suffix":""},{"id":259151162,"identity":"2cb8b4d5-0188-4845-97bd-720ac22fb03d","order_by":4,"name":"Amancio Carnero","email":"","orcid":"https://orcid.org/0000-0003-4357-3979","institution":"Consejo Superior de Investigaciones Cient�ficas Instituto de Biomedicina de Sevilla: Instituto de Biomedicina de Sevilla","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Amancio","middleName":"","lastName":"Carnero","suffix":""}],"badges":[],"createdAt":"2023-12-09 13:56:14","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-3730407/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-3730407/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1186/s13046-024-02988-y","type":"published","date":"2024-02-26T15:01:38+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":48290644,"identity":"7fdfb004-a004-42a0-91cb-21cf574db373","added_by":"auto","created_at":"2023-12-15 17:28:17","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":507038,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eHypoxia leads to the overexpression of \u003c/strong\u003e\u003cem\u003e\u003cstrong\u003ePLD2\u003c/strong\u003e\u003c/em\u003e\u003cstrong\u003e through HIF-1a in OC. (A)\u003c/strong\u003e \u003cem\u003ePLD2\u003c/em\u003e expression in the OC patient databases GSE18520, GSE40595, GSE38666, GSE12172, GSE9891, GSE63885 and GSE2109. \u003cstrong\u003e(B)\u003c/strong\u003e \u003cem\u003ePLD2\u003c/em\u003e expression in the ovarian cancer patient databases GSE18520, GSE40595, and GSE38666\u003cstrong\u003e. \u003c/strong\u003eIn \u003cstrong\u003e(A)\u003c/strong\u003e and \u003cstrong\u003e(B), \u003c/strong\u003ethe box plots show the \u003cem\u003ePLD2\u003c/em\u003e expression levels in ovarian tumor tissue (blue) or nontumor tissue (red) from patients.\u003cstrong\u003e (C)\u003c/strong\u003e Kaplan–Meier plot showing the overall survival of patients with high (red) or low (green) PLD2 expression levels in the OC patient database GSE19161. The data were analysed with a log-rank test, and the associated P values are shown in the graph. For A to C, expression levels are shown as log\u003csub\u003e2\u003c/sub\u003e transformed values from the R2 database. \u003cstrong\u003e(D)\u003c/strong\u003e \u003cem\u003ePLD2\u003c/em\u003e, \u003cem\u003eLDHA\u003c/em\u003e and\u003cem\u003e VEGFA \u003c/em\u003eexpression levels in SKOV3, OVCAR8 and ES-2 ovarian cancer cell lines under normoxia (black) and hypoxia (white) or in the presence of the HIF hydroxylase inhibitor DMOG (grey), measured by RT-qPCR. The mRNA expression was calculated as 2\u003csup\u003e-∆Ct\u003c/sup\u003e relative to the \u003cem\u003eACTB\u003c/em\u003e gene. \u003cstrong\u003e(E)\u003c/strong\u003e Representative images and quantification of PLD2 protein levels by immunofluorescence in SKOV3, OVCAR8 and ES-2 cells carrying empty vector (Ev) under normoxia (black) and hypoxia (white) or in the presence of the HIF-hydroxylase inhibitor DMOG (grey). Scale bars: 10 μm. \u003cstrong\u003e(F)\u003c/strong\u003e Top, heatmaps showing the microC signal in hESCs in a 0.5-Mb region of chromosome 17, microC hESC boundaries (green), TAD containing the \u003cem\u003ePLD2\u003c/em\u003e gene (red) and genes (blue). Bottom, zoom within the TAD containing the \u003cem\u003ePLD2\u003c/em\u003e gene showing ENCODE cis-regulatory elements (CREs, yellow to red) and those containing HIF1 motifs with relative scores higher than 0.9 (black). \u003cstrong\u003e(G) \u003c/strong\u003eLuciferase activity assay of the \u003cem\u003ePLD2\u003c/em\u003e promoter and putative enhancer in HEK293 cells under hypoxic or normoxic conditions. \u003cstrong\u003e(H)\u003c/strong\u003e Western blot showing PLD2, HIF-1a, and alpha-tubulin protein levels in SKOV3 or OVCAR8 ovarian cancer cells in the presence or absence of a small interfering RNA (siRNA) of HIF-1a. \u003cstrong\u003e(I)\u003c/strong\u003e Western blot showing PLD2, HIF-1a, and alpha-tubulin protein levels in SKOV3 and OVCAR8 cells expressing a \u003cem\u003ehif1a\u003c/em\u003e mutant (HIF mut) or Ev. \u003cstrong\u003e(J)\u003c/strong\u003e Representative images and quantification of PLD2 protein levels by immunofluorescence in SKOV3, OVCAR8 and ES-2 cells expressing a \u003cem\u003ehif1a\u003c/em\u003e mutant (Hif Mut) under normoxia (black) and hypoxia (white). Scale bars: 10 μm. A minimum of three biological were performed per each experiment. The data were compared using Student’s t tests. For D, E and G, asterisks indicate statistical significance with respect to normoxia. For J, asterisks indicate statistical significance with respect to normoxia in panel B. *p \u0026lt; 0.05; **p \u0026lt; 0.01; ***p \u0026lt; 0.001.\u003c/p\u003e","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/rs-3730407/v1/fc87a4ec4fc17102124a9460.png"},{"id":48290641,"identity":"1ad51b14-acc7-4d2d-9287-91250695bed8","added_by":"auto","created_at":"2023-12-15 17:28:17","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":523082,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eHypoxia and PLD2 modify the chromatin landscape of OC cells. (A) \u003c/strong\u003eDifferential analyses of accessibility between OC cells carrying Ev under normoxic or hypoxic conditions from ATAC-seq data (n = 2 biological replicates per condition). The log2 normalized read counts of peaks versus the log2-fold-change of accessibility are plotted. Peaks showing a statistically significant change (P value \u0026lt; 0.05) are highlighted in blue (hypoxia increased peaks) or red (normoxia increased peaks). \u003cstrong\u003e(B)\u003c/strong\u003e Heatmaps plotting normalized ATAC-seq signals at differentially accessible regions (DARs) from (f) in OC cells carrying Ev or expressing \u003cem\u003ePLD2\u003c/em\u003e under normoxia and carrying Ev or expressing \u003cem\u003eshPLD2\u003c/em\u003e under hypoxia. \u003cstrong\u003e(C)\u003c/strong\u003e Differential transcription factor (TF) binding analysis in OC cells carrying Ev under normoxia vs. hypoxia conditions using TOBIAS. Volcano plot showing the differential binding score versus the -log\u003csub\u003e10\u003c/sub\u003e p value. TF motifs with increased (blue) or decreased (red) binding in hypoxia are highlighted. \u003cstrong\u003e(D)\u003c/strong\u003e Clustering of TFs with increased (top) or decreased (bottom) binding in hypoxia. \u003cstrong\u003e(E)\u003c/strong\u003e Tracks with ATAC-seq signals in OC cells carrying Ev in hypoxia or normoxia, expressing \u003cem\u003ePLD2\u003c/em\u003e in normoxia or expressing \u003cem\u003eshPLD2\u003c/em\u003e in hypoxia at the \u003cem\u003eSOX9,\u003c/em\u003e \u003cem\u003ePROM1, and WNT7A loci.\u003c/em\u003e\u003c/p\u003e","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-3730407/v1/a09f2cefc130ef5fc8f52ff8.png"},{"id":48290640,"identity":"7116f8f8-bf06-4612-b905-85ddc109a7e5","added_by":"auto","created_at":"2023-12-15 17:28:17","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":638069,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eHigh PLD2 expression in OC patients leads to a transcriptomic rewiring of hypoxia and stemness genes. (A)\u003c/strong\u003eHeatmaps showing the expression z scores of stem-associated genes or hypoxia response genes whose expression was correlated with \u003cem\u003ePLD2\u003c/em\u003e in the GSE18520 OC patient database. \u003cstrong\u003e(B)\u003c/strong\u003e Expression levels of \u003cem\u003ePLD2\u003c/em\u003e and selected stem-associated or hypoxia response genes correlated with \u003cem\u003ePLD2\u003c/em\u003e in the OC patient database GSE18520. Box plots showing gene expression in patients in tumoral Cluster 1 (T1; light blue), tumoral Cluster 2 (T2; dark blue), or nontumor tissue (NT; red). Box plots representing the centreline, median; box limits, 25th and 75th percentiles. For A, B, C and E, data were compared using Student’s t tests. Asterisks indicate statistical significance with respect to NT tissue. *p \u0026lt; 0.05; **p \u0026lt; 0.01; ***p \u0026lt; 0.001.\u003c/p\u003e","description":"","filename":"3.png","url":"https://assets-eu.researchsquare.com/files/rs-3730407/v1/d3d867750de56f7f51e832b8.png"},{"id":48291416,"identity":"ca6b619c-1238-48d4-a277-cc298b53e7fc","added_by":"auto","created_at":"2023-12-15 17:36:17","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":498855,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cem\u003e\u003cstrong\u003ePLD2\u003c/strong\u003e\u003c/em\u003e\u003cstrong\u003e promotes tumorigenesis and CSC-like features in ovarian cancer. (A)\u003c/strong\u003e Analysis of the expression of \u003cem\u003ePLD2\u003c/em\u003e by RT‒qPCR in SKOV3, OVCAR8 and ES-2 cells carrying an empty vector (Ev) expressing \u003cem\u003ePLD2\u003c/em\u003e or \u003cem\u003eshPLD2\u003c/em\u003e. The mRNA expression was calculated as 2\u003csup\u003e-∆Ct\u003c/sup\u003e relative to the \u003cem\u003eACTB\u003c/em\u003e gene. \u003cstrong\u003e(B)\u003c/strong\u003e Western blot showing the protein levels of PLD2 in cells carrying Ev and expressing \u003cem\u003ePLD2\u003c/em\u003e or \u003cem\u003eshPLD2\u003c/em\u003e. \u003cstrong\u003e(C)\u003c/strong\u003e Quantification of the number of colonies formed by SKOV3, OVCAR8 and ES-2 cells carrying Ev and expressing \u003cem\u003ePLD2\u003c/em\u003e or \u003cem\u003eshPLD2\u003c/em\u003e. \u003cstrong\u003e(D)\u003c/strong\u003e Growth curve of SKOV3, OVCAR8 and ES-2 cells carrying Ev (light green) or expressing \u003cem\u003ePLD2\u003c/em\u003e (dark green) or \u003cem\u003eshPLD2\u003c/em\u003e (blue) represented as the accumulation of the doubling times. \u003cstrong\u003e(E)\u003c/strong\u003e Tumor growth in xenografts from SKOV3 and OVCAR8 cells carrying Ev, overexpressing \u003cem\u003ePLD2\u003c/em\u003e or \u003cem\u003eshPLD2\u003c/em\u003e were injected into female athymic nude mice. Cohorts of 5 mice each were used. \u003cstrong\u003e(F)\u003c/strong\u003e Percentage of paraclones, meroclones and holoclones formed by SKOV3, OVCAR8 and ES-2 cells carrying Ev expressing \u003cem\u003ePLD2\u003c/em\u003e or \u003cem\u003eshPLD2\u003c/em\u003e. At least 200 individual clones were analysed. The averages and SDs of three independent experiments are shown. \u003cstrong\u003e(G)\u003c/strong\u003e Top, Representative images of tumorspheres formed by SKOV3, OVCAR8 and ES-2 cells carrying Ev expressing \u003cem\u003ePLD2\u003c/em\u003e or \u003cem\u003eshPLD2\u003c/em\u003e. Bottom, quantification of the number and size of tumorspheres. Scale bars: 250 μm. \u003cstrong\u003e(H)\u003c/strong\u003e Analysis of the expression by RT‒qPCR of stemness-associated genes and \u003cem\u003ePLD2\u003c/em\u003e in total cell extracts and tumorspheres from SKOV3, OVCAR8 and ES-2 cells carrying Ev and expressing \u003cem\u003ePLD2\u003c/em\u003e or \u003cem\u003eshPLD2\u003c/em\u003e. The mRNA expression was calculated as 2\u003csup\u003e-∆Ct\u003c/sup\u003e relative to the \u003cem\u003eACTB\u003c/em\u003e gene. The average and SD of at least three independent experiments are shown. The data were analysed using Student’s t test. Asterisks indicate statistical\u003c/p\u003e","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/rs-3730407/v1/bd00eaa18e9aa8f223f8f8ad.png"},{"id":48290645,"identity":"cd5508c6-79bc-4ea1-b4ad-c5ab4aef4d7c","added_by":"auto","created_at":"2023-12-15 17:28:17","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":573965,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eThe hypoxia-induced stemness of OC cells depends on PLD2 expression. (A)\u003c/strong\u003e Heatmaps showing the z scores of the expression of stemness-associated genes obtained from TaqMan arrays. Genes are sorted according to decreasing z scores in Ev-carrying cells under normoxia. \u003cstrong\u003e(B)\u003c/strong\u003e Heatmaps showing the z-scores of stemness genes expression levels in SKOV3 cells carrying EV or plasmid overexpressing PLD2 under normoxia and carrying EV or plasmid expressing shPLD2 under hypoxia. Hierarchical clustering of the samples of (A) is shown. \u003cstrong\u003e(C)\u003c/strong\u003e Expression levels of \u003cem\u003eSOX2, NANOG, CD44 \u003c/em\u003eand\u003cem\u003e EPCAM\u003c/em\u003e stemness-associated genes in SKOV3, OVCAR8 and ES-2 cells carrying Ev and expressing \u003cem\u003ePLD2\u003c/em\u003e, \u003cem\u003eshPLD2\u003c/em\u003e or both\u003cem\u003e \u003c/em\u003eunder normoxic or hypoxic conditions. The mRNA expression was calculated as 2\u003csup\u003e-∆Ct\u003c/sup\u003e relative to the \u003cem\u003eACTB\u003c/em\u003e gene.\u003cstrong\u003e (D) \u003c/strong\u003eQuantification of the number and size of tumorspheres formed by SKOV3, OVCAR8 and ES-2 cells carrying Ev and expressing \u003cem\u003ePLD2\u003c/em\u003e, \u003cem\u003eshPLD2\u003c/em\u003e or both. \u003cstrong\u003e(E)\u003c/strong\u003e Top, determination of the Sox9, Notch1, Sox2, and Sox17 protein levels by immunofluorescence in tumorspheres formed by OC cell lines carrying Ev and expressing \u003cem\u003ePLD2\u003c/em\u003e, \u003cem\u003eshPLD2\u003c/em\u003e or both. Bottom, quantification of the percentage of cells with Sox9 and Notch1 expression in tumorspheres formed by ovarian cancer cell lines carrying Ev and expressing \u003cem\u003ePLD2\u003c/em\u003e or \u003cem\u003eshPLD2\u003c/em\u003e. Scale bars: 100 μm. \u003cstrong\u003e(F)\u003c/strong\u003e Expression levels of \u003cem\u003eSOX9, NOTCH1, SOX2 \u003c/em\u003eand\u003cem\u003e SOX17 \u003c/em\u003estemness-associated genes in tumorspheres formed by ovarian cancer cell lines carrying Ev and expressing \u003cem\u003ePLD2\u003c/em\u003e, \u003cem\u003eshPLD2\u003c/em\u003e or both. The mRNA expression was calculated as 2\u003csup\u003e-∆Ct\u003c/sup\u003e relative to the \u003cem\u003eACTB\u003c/em\u003e gene. A minimum of three biological were performed per each experiment. The data were analysed using Student’s t test. Asterisks indicate statistical significance with respect to Ev carrying cells, unless indicated by horizontal lines. *, P \u0026lt; 0.05; **, P \u0026lt; 0.01; ***, P \u0026lt; 0.001.\u003c/p\u003e","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-3730407/v1/0a5a0b5ee72d05614054a8fd.png"},{"id":48290639,"identity":"e4c0d799-0efa-4d1b-9957-bd95871716db","added_by":"auto","created_at":"2023-12-15 17:28:17","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":335010,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eHypoxia-mediated reprogramming to pluripotent stem cells is dependent on PLD2. (A)\u003c/strong\u003e Representative images of induced pluripotent stem cells (iPSCs) by reprogramming mouse embryonic fibroblasts (MEFs) with OSKM genes (Yamanaka factors Oct3/4, Sox2, Klf4 and cMyc) and Nanog reporter retroviruses. Top, bright field microscopy of iPSCs. Medium, immunofluorescence image showing GFP in cells in which \u003cem\u003eNanog\u003c/em\u003e is active. Bottom, Alkaline phosphatase activity. \u003cstrong\u003e(B) \u003c/strong\u003eLeft, quantification of iPSCs generated from MEFs infected with OSKM genes and an additional lentiviral vector that expresses GFP in cells in which the Nanog promoter/enhancer is active and the corresponding plasmid overexpressing \u003cem\u003ePLD2\u003c/em\u003e, \u003cem\u003eshPLD2\u003c/em\u003e, both or carrying Ev under normoxia or hypoxia. Right, quantification of cells with alkaline phosphatase activity at the end of the iPSC generation experiment. A minimum of three biological were performed per each experiment. The data were analysed using Student’s t test. Asterisks indicate statistical significance with respect to Ev carrying cells, unless indicated by horizontal lines. *, P \u0026lt; 0.05; **, P \u0026lt; 0.01; ***, P \u0026lt; 0.001.\u003c/p\u003e","description":"","filename":"6.png","url":"https://assets-eu.researchsquare.com/files/rs-3730407/v1/05ab197864e8ae75e439849c.png"},{"id":48290643,"identity":"dab856e4-8453-4364-bf61-74bccbf48582","added_by":"auto","created_at":"2023-12-15 17:28:17","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":477206,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cem\u003e\u003cstrong\u003ePLD2\u003c/strong\u003e\u003c/em\u003e\u003cstrong\u003e overexpression is associated with resistance to treatment in ovarian cancer. (A) \u003c/strong\u003eRepresentative images of PLD2 immunostaining in ovarian cancer and non-tumoral samples. Scale bars, 50 µm\u003cem\u003e. \u003c/em\u003e\u003cstrong\u003e(B)\u003c/strong\u003e \u003cem\u003ePLD2\u003c/em\u003e expression in ovarian patients from the HUVR-IBIS database. Box plots showing the \u003cem\u003ePLD2\u003c/em\u003e expression levels in ovarian tumor tissue (black) or nontumor tissue (red) from patients. Box plots representing the centreline, median; box limits, 25th and 75th percentiles; whiskers, minimum and maximum values. \u003cstrong\u003e(C)\u003c/strong\u003e Analysis of \u003cem\u003ePLD2\u003c/em\u003e expression by RT‒qPCR in a cohort of OC patients who were sensitive (S; pink) or resistant (R; green) to platinum treatment (HUVR-IBIS). For B and C, the mRNA expression was calculated as 2\u003csup\u003e-∆Ct\u003c/sup\u003e relative to the \u003cem\u003eACTB\u003c/em\u003e gene. \u003cstrong\u003e(D-E)\u003c/strong\u003e Kaplan‒Meier plots showing overall or progression-free survival in patients who were sensitive (pink) or resistant (green) to platinum treatment in the HUVR-IBIS cohort. \u003cstrong\u003e(F)\u003c/strong\u003e Determination of the IC\u003csub\u003e50 \u003c/sub\u003eof cis-platinum and carboplatin in ES-2, SKOV3 and OVCAR8 cells overexpressing \u003cem\u003ePLD2\u003c/em\u003e or \u003cem\u003eshPLD2\u003c/em\u003e or carrying EV in normoxia and hypoxia, in combination or without the PLD inhibitor (PLDi) FIPI in ES-2, SKOV3 and OVCAR8 cells overexpressing \u003cem\u003eshPLD2\u003c/em\u003e or \u003cem\u003ePLD2\u003c/em\u003e or carrying EV in normoxia or hypoxia.\u003cstrong\u003e (G-H) \u003c/strong\u003eDetermination of the tumor volume (G) and survival (H) after treatment with saline, cisplatin, PLDi or both in xenografts of SKOV3 cells expressing \u003cem\u003eshPLD2, PLD2\u003c/em\u003e or EV. A minimum of three biological were performed per each experiment. The data were compared using Student’s t tests. For F to H, asterisks indicate statistical significance with respect to Ev carrying cells or xenografts, unless indicated by horizontal lines. *p \u0026lt; 0.05; **p \u0026lt; 0.01; ***p \u0026lt; 0.001.\u003c/p\u003e","description":"","filename":"7.png","url":"https://assets-eu.researchsquare.com/files/rs-3730407/v1/3eb8b016fe5edaab4044a353.png"},{"id":51958543,"identity":"6bb7736c-1691-430f-9571-dc19298f3cb4","added_by":"auto","created_at":"2024-03-04 15:17:19","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":3988930,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-3730407/v1/094d160e-8fd4-46d0-ac5e-d0113516801d.pdf"},{"id":48290647,"identity":"ebb910f7-44fc-4a1e-8f81-19cbc12ab981","added_by":"auto","created_at":"2023-12-15 17:28:18","extension":"docx","order_by":4,"title":"","display":"","copyAsset":false,"role":"supplement","size":15883554,"visible":true,"origin":"","legend":"","description":"","filename":"MunozGalvansupplementary.docx","url":"https://assets-eu.researchsquare.com/files/rs-3730407/v1/a5aa588eb828703d746cea60.docx"}],"financialInterests":"","formattedTitle":"Suppression of hypoxia-induced stemness and chemoresistance in ovarian tumors","fulltext":[{"header":"BACKGROUND","content":"\u003cp\u003eHypoxia is a common feature of the tumor microenvironment in solid tumors (Harris, \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2002\u003c/span\u003e). Hypoxia is generated by insufficient oxygen diffusion towards parts of the tumor that are not irrigated, although the highly irregular tumor microvasculature may also generate hypoxic regions. Cancer cells under hypoxic conditions undergo a series of transcriptional changes that are induced by hypoxia-inducible factors (HIF), HIF-1, -2 and \u0026minus;\u0026thinsp;3, with HIF-1 being the best known. HIFs are heterodimeric transcription factors containing bHLH-PAS domains and consist of an α subunit regulated by oxygen levels and a constitutively expressed β subunit (also called ARNT). The HIFa-ARNT dimers bind DNA at specific sequences known as hypoxia-response elements (HREs) to promote target gene expression (Schito and Semenza, \u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). In normoxia, HIF-α subunits are hydroxylated at specific proline and asparagine residues by prolyl hydroxylase 2 (PHD2) and are recognized and targeted for degradation by the von Hippel‒Lindau (VHL) tumor suppressor. However, hypoxia leads to the inhibition of such hydroxylation and subsequent HIF-α accumulation, heterodimerization with ARNT and transcriptional activation of target genes (Kaelin and Ratcliffe, \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e2008\u003c/span\u003e). The multiple effects of hypoxia on the biology of tumors include preventing apoptosis and promoting proliferation and autophagy, inducing metabolic alterations, and promoting angiogenesis, the epithelial-to-mesenchymal (EMT) transition, invasion and metastasis (Jing et al., \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2019\u003c/span\u003e; Wilson and Hay, \u003cspan citationid=\"CR72\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). Indeed, HIF-1α is overexpressed in multiple tumor types and is associated with a poor prognosis in patients and therapy resistance (Chen et al., \u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e2018a\u003c/span\u003e; Klemba et al., \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Lin and Koong, \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Mayer et al., \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2008\u003c/span\u003e; Simiantonaki et al., \u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e2008\u003c/span\u003e; Sun et al., \u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Ueda et al., \u003cspan citationid=\"CR68\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Unruh et al., \u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e2003\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eOvarian cancer (OC) is the most lethal gynecological malignancy (Siegel et al., \u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e2020\u003c/span\u003e) mainly due to its nonspecific clinical manifestations, which lead to a late diagnosis and high chemoresistance (Torre et al., \u003cspan citationid=\"CR66\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). Similar to other solid tumors, hypoxia is a key modulator of the tumor microenvironment in OC, affecting not only the primary tumor but also the ascitic fluid, which is the main route of spread of these malignancies and shows low oxygen levels (Klemba et al., \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Munoz-Galvan and Carnero, \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Therefore, ovarian tumors are highly hypoxia-dependent and this dependency influence the response to treatment. A major cause of chemoresistance is the persistence of a cancer cell subpopulation known as cancer stem cells (CSCs) or tumor-initiating cells. CSCs are able to recapitulate a new tumor since they possess self-renewal and pluripotency properties similar to those of normal stem cells (Batlle and Clevers, \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Beck and Blanpain, \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). CSCs are resistant to common antitumor therapies, which may indeed cause their enrichment, leading to chemoresistance and relapse (Carnero et al., \u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Maugeri-Sacca et al., \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). Several studies indicate that HIF-1 is required for maintaining CSCs and that its activation in hypoxia leads to the increased expression of stem marker genes in multiple cancer types (Mathieu et al., \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2011\u003c/span\u003e; Peng and Liu, \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). Therefore, hypoxia may promote the generation of CSCs leading to chemoresistance through HIF factors, but our knowledge of HIF targets that may be responsible for CSC activity is limited.\u003c/p\u003e \u003cp\u003eWe previously showed that the \u003cem\u003ePLD2\u003c/em\u003e gene, encoding phospholipase D2, is overexpressed in colorectal cancer and induces stemness in cancer cells through communication with the tumor microenvironment (Munoz-Galvan et al., \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2019b\u003c/span\u003e). Here we find that \u003cem\u003ePLD2\u003c/em\u003e is also overexpressed in ovarian cancer patients, being associated with poor patient survival. We explore a possible connection of \u003cem\u003ePLD2\u003c/em\u003e with hypoxia in OC and demonstrate that \u003cem\u003ePLD2\u003c/em\u003e expression in OC cells is stimulated by hypoxia and that HIF-1α promotes \u003cem\u003ePLD2\u003c/em\u003e transcription through HREs at its promoter and an intronic enhancer. This leads to increased chromatin accessibility around stemness genes, which are in turn overexpressed provoking an enhancement of tumor growth and formation of CSC-like tumorspheres. This is corroborated by induced pluripotent stem cell (iPSC) reprogramming experiments that confirm the role of hypoxia and PLD2 in dedifferentiation of OC cells to stem-like cells. Importantly, the increase in CSC-like features induced by hypoxia relies on \u003cem\u003ePLD2\u003c/em\u003e expression, indicating that the hypoxia-PLD2 axis is a major contributor to tumor stemness. We confirm these findings in transcriptional databases of OC patients, where high \u003cem\u003ePLD2\u003c/em\u003e expression is associated with a transcriptional rewiring of genes involved in the hypoxia response and in the maintenance of stem cells. Finally, we show that \u003cem\u003ePLD2\u003c/em\u003e overexpression causes resistance to platinum-based compounds and propose a new therapy based on pharmacological inhibition of phospholipases D to suppress such chemoresistance.\u003c/p\u003e"},{"header":"RESULTS","content":"\u003cp\u003e\u003cstrong\u003ePLD2 is overexpressed in OC patients and in ovarian cancer cells under hypoxic conditions.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eFirst, we wondered whether \u003cem\u003ePLD2\u003c/em\u003e was overexpressed in OC patients. For this, we analysed \u003cem\u003ePLD2\u003c/em\u003e expression in 7 OC and one ovarian tissue databases using the R2 platform and found that \u003cem\u003ePLD2\u003c/em\u003e expression was significantly higher in tumoral than in non-tumoral samples (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eA). This result was confirmed by a comparison between patients and control individuals in 3 databases containing their own controls (GSE18520, GSE4595 and GSE3866) (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eB). Next, we wondered whether \u003cem\u003ePLD2\u003c/em\u003e expression was associated with patient survival and analysed 4 OC databases with available overall survival (OS) data (GSE13876, GSE19161, GSE23554 and GSE31245). We split patients in low-risk and high-risk groups based on their OS, but the difference in OS between both groups was statistically significant only in one database, GSE19161 (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eC; Supplementary Figure \u003cspan class=\"InternalRef\"\u003eS1\u003c/span\u003e). However, the expression levels of \u003cem\u003ePLD2\u003c/em\u003e were significantly higher in the high-risk group than in the low-risk group in the 4 databases (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eC; Supplementary Figure \u003cspan class=\"InternalRef\"\u003eS1\u003c/span\u003e), indicating that \u003cem\u003ePLD2\u003c/em\u003e is commonly overexpressed in OC patients and may be associated with decreased patient survival.\u003c/p\u003e\n\u003cp\u003eSimilar to other solid tumors, OC shows hypoxic areas, and its main dissemination site, the ascitic fluid, is also characterized by hypoxia (Klemba et al., \u003cspan class=\"CitationRef\"\u003e2020\u003c/span\u003e). Since \u003cem\u003ePLD2\u003c/em\u003e expression is increased under hypoxic conditions in colon cancer cells (Liu et al., \u003cspan class=\"CitationRef\"\u003e2020\u003c/span\u003e) and is involved in tumor stemness through communication with the microenvironment in colorectal cancer (Munoz-Galvan et al., \u003cspan class=\"CitationRef\"\u003e2019b\u003c/span\u003e), we wondered whether hypoxia could increase the expression of \u003cem\u003ePLD2\u003c/em\u003e in OC. For this, we selected SKOV3, OVCAR8 and ES-2 OC cell lines for our studies and analysed the expression levels of \u003cem\u003ePLD2\u003c/em\u003e under hypoxic conditions. We found that \u003cem\u003ePLD2\u003c/em\u003e expression was similar among the three cell lines and that oxygen levels of 3% led to a 2-fold increase in the expression of \u003cem\u003ePLD2\u003c/em\u003e in all of them, while the well-known hypoxia target genes \u003cem\u003eLDHA\u003c/em\u003e and \u003cem\u003eVEGFA\u003c/em\u003e showed a similar increase, validating the results (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eD). The hypoxia-induced increase in PLD2 expression was also observed at the protein level by performing immunofluorescence in the three OC cell lines (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eE). In addition, we further demonstrated these findings by using the HIF-hydroxylase inhibitor DMOG, which increases HIF levels, generating a hypoxia-like phenotype under normoxic conditions (Supplementary Figure S2A-B). Thus, we treated OC cells with DMOG and found that this treatment resulted in a similar increase in hypoxia marker gene expression as that induced by hypoxia, also leading to the observed increase in \u003cem\u003ePLD2\u003c/em\u003e expression at both the mRNA and protein levels (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eD-E). Therefore, we can conclude that \u003cem\u003ePLD2\u003c/em\u003e expression is promoted by hypoxia.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eHIF-1\u0026alpha; activates\u003c/strong\u003e \u003cstrong\u003ePLD2\u003c/strong\u003e \u003cstrong\u003etranscription through HREs at promoter and hypoxia-specific enhancer regions.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe aimed to understand how hypoxia promotes \u003cem\u003ePLD2\u003c/em\u003e expression. HIF-1\u0026alpha; is considered the master transcriptional regulator of the cellular response to hypoxia. It forms a heterodimer with ARNT that binds HREs to control the expression of hypoxia-response genes (Schito and Semenza, \u003cspan class=\"CitationRef\"\u003e2016\u003c/span\u003e). To address the possibility that HIF-1\u0026alpha; could regulate \u003cem\u003ePLD2\u003c/em\u003e expression at the transcriptional level, we first searched for possible cis-regulatory elements (CREs) near the \u003cem\u003ePLD2\u003c/em\u003e gene that may contain HREs. According to public 3D chromatin conformation experiments (micro-C) in human embryonic stem cells, the \u003cem\u003ePLD2\u003c/em\u003e gene is located within a topologically associating domain (TAD) of 190 kb with a high interaction frequency at the 3D level and is relatively isolated from the neighbouring regions (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eF). To identify CREs that may regulate \u003cem\u003ePLD2\u003c/em\u003e expression, we focused on a smaller region of 50 kb surrounding the \u003cem\u003ePLD2\u003c/em\u003e gene with a higher interaction frequency with the \u003cem\u003ePLD2\u003c/em\u003e promoter, which we called the \u003cem\u003ePLD2\u003c/em\u003e regulatory region. We scanned the CREs annotated by ENCODE within the \u003cem\u003ePLD2\u003c/em\u003e regulatory region for the presence of the DNA binding motif of HIF1A with a high score (\u0026gt;\u0026thinsp;90% relative score). We found 15 out of 35 CREs fulfilling this condition, some of which corresponded to gene promoters and others to enhancers, including the \u003cem\u003ePLD2\u003c/em\u003e promoter and an enhancer in \u003cem\u003ePLD2\u003c/em\u003e intron 12 (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eF). These CREs with high-score HIF1A motifs represent putative HREs.\u003c/p\u003e\n\u003cp\u003eTo assess the regulatory activity of these two CREs in the \u003cem\u003ePLD2\u003c/em\u003e gene (promoter and putative enhancer) containing the HIF1A motif, we cloned both genomic regions in promoter and enhancer reporter vectors controlling the expression of the luciferase gene. We transfected OC cells with these vectors and measured the luciferase activity under normoxia and hypoxia. We found that the \u003cem\u003ePLD2\u003c/em\u003e promoter was able to activate reporter expression in normoxia and that hypoxia led to a significant increase in luciferase activity (8-fold) (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eG), suggesting that the \u003cem\u003ePLD2\u003c/em\u003e promoter responds to hypoxic conditions by increasing \u003cem\u003ePLD2\u003c/em\u003e transcription. However, the enhancer contained within the \u003cem\u003ePLD2\u003c/em\u003e intron was unable to activate reporter expression in normoxia but led to a 22-fold increase in its expression in hypoxia (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eG), indicating that this enhancer acts as an HRE in OC cells.\u003c/p\u003e\n\u003cp\u003eNext, we wanted to validate functionally that the hypoxia-induced upregulation of \u003cem\u003ePLD2\u003c/em\u003e expression was indeed mediated by HIF-1\u0026alpha;. Therefore, we depleted \u003cem\u003eHIF1A\u003c/em\u003e using a small interfering RNA (siRNA) in the three OC cell lines and found that it suppressed the increase in the PLD2 protein levels induced by hypoxia (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eH). Additionally, we generated \u003cem\u003ehif1a\u003c/em\u003e mutant OC cell lines by transfecting cells with a mutant \u003cem\u003ehif1a\u003c/em\u003e allele (Hu et al., \u003cspan class=\"CitationRef\"\u003e2007\u003c/span\u003e) that is unable to be hydroxylated and, therefore, is constitutively active even under normoxic conditions. Interestingly, we observed that the PLD2 levels in \u003cem\u003ehif1a\u003c/em\u003e mutant cells were as high as those induced by hypoxia on a \u003cem\u003eHIF1\u0026alpha;\u003c/em\u003e wild-type background in both normoxia and hypoxia, confirming that HIF-1\u0026alpha; activates PLD2 expression (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eI-J). Altogether, these data indicate that hypoxia induces \u003cem\u003ePLD2\u003c/em\u003e expression in OC cells through transcriptional activation by HIF-1\u0026alpha; at HREs in the \u003cem\u003ePLD2\u003c/em\u003e gene.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eHypoxia alters the chromatin landscape of ovarian cancer cells in a PLD2-dependent manner.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe wondered whether \u003cem\u003ePLD2\u003c/em\u003e expression mediated by HIF-1\u0026alpha; could have an impact on hypoxia-induced gene regulation. For this, we first analyzed the effect of hypoxia and \u003cem\u003ePLD2\u003c/em\u003e expression in the epigenomic landscape of OC cells through ATAC-seq experiments in SKOV3 cells under normoxia and hypoxia and altered \u003cem\u003ePLD2\u003c/em\u003e expression under normoxia (\u003cem\u003ePLD2\u003c/em\u003e overexpression) and hypoxia (\u003cem\u003ePLD2\u003c/em\u003e depletion). We computationally called open chromatin regions (ATAC peaks) in normoxia and hypoxia and compared both conditions by a differential accessibility analysis; we detected 140 and 102 peaks with increased or decreased accessibility in hypoxia, respectively (Fig. \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eA). The heatmaps and aggregate profiles of these differentially accessible regions (DARs) showed that the peaks with increased accessibility in hypoxia were also more open upon \u003cem\u003ePLD2\u003c/em\u003e overexpression in normoxia, although to a lower extent, and vice versa, with the peaks showing decreased accessibility in hypoxia (Fig. \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eB), suggesting that \u003cem\u003ePLD2\u003c/em\u003e overexpression in normoxia has a similar effect on chromatin accessibility as hypoxia. Indeed, the DARs of EV- versus \u003cem\u003ePLD2\u003c/em\u003e-overexpressing cells under normoxia showed similar changes in accessibility under hypoxia, reinforcing the previous idea (Supplementary Figure S3A-B). Both DARs in normoxia versus hypoxia and control versus \u003cem\u003ePLD2\u003c/em\u003e overexpression were associated with genes enriched in Gene Ontology terms related to the response to hypoxia or well-known functions of \u003cem\u003ePLD2\u003c/em\u003e, respectively (Supplementary Figure S3C-D). Moreover, changes in accessibility induced by hypoxia were suppressed by \u003cem\u003ePLD2\u003c/em\u003e depletion (Fig. \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eB), suggesting that the effect of hypoxia in chromatin accessibility is mediated by PLD2. Similarly, changes in accessibility upon \u003cem\u003ePLD2\u003c/em\u003e overexpression in normoxia were absent in \u003cem\u003ePLD2\u003c/em\u003e-depleted cells in hypoxia (Supplementary Figure S3B), reinforcing the previous idea. Altogether, these data indicate that both hypoxia and \u003cem\u003ePLD2\u003c/em\u003e overexpression induce similar alterations in the chromatin landscape of OC cells and that the effect of hypoxia is mediated by PLD2.\u003c/p\u003e\n\u003cp\u003eNext, we sought to investigate the possible mechanisms driving the alterations in the chromatin accessibility landscape induced by hypoxia and \u003cem\u003ePLD2\u003c/em\u003e overexpression. For this, we first performed motif enrichment analyses of DARs. We found that DARs in hypoxia were enriched in the motifs of the AP-1, ETS and C2H2 zinc finger transcription factor families in the increased accessibility sites and the C2H2 zinc finger and fork head families in the decreased accessibility sites (Supplementary Figure S4A). Similar enrichments were found in DARs upon \u003cem\u003ePLD2\u003c/em\u003e overexpression (Supplementary Figure S4B), reinforcing the idea of a similar effect under both conditions. To further elucidate the possible TFs involved in hypoxia and PLD2-mediated epigenomic changes, we estimated differential TF binding among the conditions based on footprints in ATAC-seq data. Using this approach, we confirmed the increased chromatin binding of AP-1 family transcription factors, such as FOS and JUN, consistent with their implication in the response to hypoxia (Kunz and Ibrahim, \u003cspan class=\"CitationRef\"\u003e2003\u003c/span\u003e). (Fig. \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eC-D). We also found other TF families showing increased TF binding in hypoxia, such as homeobox, paired box or fork head TFs, and the C2H2 zinc finger TF ZBTB32, while TF families with decreased binding in hypoxia included members of the bHLH, NF-Y and ETS families, among others (Supplementary Figure S4C). When we compared these TFs with those showing increased binding upon \u003cem\u003ePLD2\u003c/em\u003e overexpression in normoxia, we found that most of these TFs overlapped (Supplementary Figure S4C). In particular, 25 TF motifs of the AP-1 family showed increased binding under both conditions. Interestingly, these AP-1 and 20 more TF motifs with increased binding in hypoxia showed decreased binding upon \u003cem\u003ePLD2\u003c/em\u003e depletion (Supplementary Figure S4C-E), indicating that they are dependent on \u003cem\u003ePLD2\u003c/em\u003e expression. Altogether, these results suggest a function of the AP-1 family of TFs mediating the alterations in the chromatin landscape mediated by hypoxia and PLD2.\u003c/p\u003e\n\u003cp\u003eSince we previously connected \u003cem\u003ePLD2\u003c/em\u003e overexpression with tumor stemness in colorectal cancer, we wondered whether alterations in the chromatin accessibility landscape of OC cells induced by hypoxia and PLD2 could result in increased expression of stemness genes. Therefore, we first selected ATAC peaks falling within the putative regulatory landscapes of the genes associated with stem cell maintenance and proliferation. Clustering of these 3,572 peaks revealed 4 groups with different accessibility levels and behaviours (Supplementary Figure S5A). Among them, Cluster 3 corresponded with peaks with increased accessibility in both hypoxic and \u003cem\u003ePLD2\u003c/em\u003e-overexpressing cells but decreased accessibility in \u003cem\u003ePLD2\u003c/em\u003e-depleted cells. Among the stemness genes associated with this cluster, we found \u003cem\u003eSOX9\u003c/em\u003e, \u003cem\u003ePROM1\u003c/em\u003e, \u003cem\u003eWNT7A\u003c/em\u003e or \u003cem\u003eJAG1\u003c/em\u003e (Fig. \u003cspan class=\"InternalRef\"\u003e2\u003c/span\u003eE; Supplementary Figure S5B). These results indicate that hypoxia promotes chromatin accessibility around stemness genes in a PLD2-dependent manner in OC cells, and suggests that hypoxia and PLD2 could be connected with tumor stemness in OC.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eHigh PLD2 expression in OC patients leads to the transcriptomic rewiring of stemness and hypoxia genes\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo determine whether \u003cem\u003ePLD2\u003c/em\u003e expression in OC patients is related to hypoxia and stemness, we analysed the expression of genes related to these functions in the three OC databases with available expression data from control individuals. First, we selected the genes annotated to the Gene Ontology (GO) term \u0026ldquo;Response to Hypoxia\u0026rdquo; and whose expression was significantly correlated with that of \u003cem\u003ePLD2\u003c/em\u003e in OC patients (p\u0026thinsp;\u0026lt;\u0026thinsp;0.05; r\u0026thinsp;\u0026gt;\u0026thinsp;0.2 or \u0026lt;-0.2). Then, we performed hierarchical clustering of patient and control individuals based on the expression levels of these genes (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eA; Supplementary Figure S6). In the GSE18520 database, the clustering clearly separated the control individuals (\u0026lsquo;non-tumoral\u0026rsquo;, NT) and a reduced group of patients that we termed \u0026lsquo;Tumoral Cluster 1\u0026rsquo; (T1) from most patients who clustered in what we termed \u0026lsquo;Tumoral Cluster 2\u0026rsquo; (T2) (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eA). The patients at T1 showed a transcriptional profile of hypoxia-related genes more similar to NT and clearly different from T2. In contrast, the clustering in the GSE4095 and GSE38666 databases clearly separated clusters of NT and tumoral (T) individuals, who showed distinct transcriptional profiles of hypoxia-related genes (Supplementary Figure S6). These results suggest that there is a switch in the expression of hypoxia-related genes in OC tumors compared with healthy ovaries.\u003c/p\u003e\n\u003cp\u003eNext, we repeated the hierarchical clustering with the genes annotated to the GO term \u0026ldquo;Stem cell maintenance\u0026rdquo; and whose expression was significantly correlated with that of \u003cem\u003ePLD2\u003c/em\u003e in OC patients (p\u0026thinsp;\u0026lt;\u0026thinsp;0.05; r\u0026thinsp;\u0026gt;\u0026thinsp;0.2 or \u0026lt;-0.2). Surprisingly, this clustering based on stem-related genes separated the patients and control individuals into the same clusters as the hypoxia-related genes (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eA; Supplementary Figure S6), suggesting a connection between both groups of genes that supports the model of CSC generation induced by hypoxia. Then, we plotted the expression levels of \u003cem\u003ePLD2\u003c/em\u003e in the 3 clusters obtained from the GSE18520 database and found that \u003cem\u003ePLD2\u003c/em\u003e exhibited significantly increased expression in Cluster T2 compared with that in Cluster NT, while Cluster T1 showed similar levels to NT (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eB). This observation suggests that a connection exists among \u003cem\u003ePLD2\u003c/em\u003e expression, the hypoxia response and stemness since \u003cem\u003ePLD2\u003c/em\u003e expression is misregulated only in patients showing transcriptional profiles highly different from healthy controls. In addition, we checked the expression of stem-associated and hypoxia-related genes in these clusters. As shown in Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eB, the expression of \u003cem\u003ePAX8\u003c/em\u003e, a well-known OC marker, was significantly increased in both patient Clusters T1 and T2, similar to other stemness genes, such as \u003cem\u003eSOX9\u003c/em\u003e, \u003cem\u003eSOX17\u003c/em\u003e, \u003cem\u003eEPCAM\u003c/em\u003e, \u003cem\u003ePROM1\u003c/em\u003e, \u003cem\u003eCD24\u003c/em\u003e, \u003cem\u003eNOTCH1\u003c/em\u003e and \u003cem\u003eWNT7A\u003c/em\u003e. Other genes in this group showed a significant\u003c/p\u003e\n\u003cp\u003eincrease in expression only in Cluster T2, coinciding with higher \u003cem\u003ePLD2\u003c/em\u003e levels, including \u003cem\u003ePAX2\u003c/em\u003e, \u003cem\u003ePOU5F1\u003c/em\u003e (\u003cem\u003eOCT4\u003c/em\u003e), \u003cem\u003eSOX5\u003c/em\u003e, \u003cem\u003eSOX11\u003c/em\u003e, \u003cem\u003eCD34\u003c/em\u003e and \u003cem\u003eTP63\u003c/em\u003e, while the others were unaffected or even significantly reduced, such as \u003cem\u003eKLF4\u003c/em\u003e, \u003cem\u003eNANOG\u003c/em\u003e or \u003cem\u003eSOX2\u003c/em\u003e, although the latter showed a nonsignificant increase in Cluster T2 (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eB). This finding suggests that there is an OC stemness signature in patients that may be stronger with higher \u003cem\u003ePLD2\u003c/em\u003e expression. In addition, some hypoxia-related genes showed a significant increase in expression in both Clusters T1 and T2, including \u003cem\u003eVEGFA\u003c/em\u003e, \u003cem\u003eSLC2A1\u003c/em\u003e (\u003cem\u003eGLUT1\u003c/em\u003e), \u003cem\u003eHK2\u003c/em\u003e and \u003cem\u003eNOX4\u003c/em\u003e, or only in T2, including \u003cem\u003eSLC2A4\u003c/em\u003e, \u003cem\u003eNOS1\u003c/em\u003e and \u003cem\u003eMMP14\u003c/em\u003e, and a nonsignificant increase in \u003cem\u003eSLC2A14\u003c/em\u003e (\u003cem\u003eGLUT14\u003c/em\u003e) (Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eB). Similar results were obtained in the NT and T clusters in the GSE4095 and GSE38666 databases (Supplementary Figure S6). Altogether, these results suggest that there is transcriptional rewiring of the expression of hypoxia- and stem-related genes in OC patients with \u003cem\u003ePLD2\u003c/em\u003e overexpression.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003ePLD2\u003c/strong\u003e \u003cstrong\u003epromotes tumorigenesis and CSC-like features in ovarian cancer cells.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe ability of hypoxia to induce a CSC-like phenotype in OC cells has been previously observed in several cancer types (Liang et al., \u003cspan class=\"CitationRef\"\u003e2012\u003c/span\u003e; Mathieu et al., \u003cspan class=\"CitationRef\"\u003e2011\u003c/span\u003e; Peng and Liu, \u003cspan class=\"CitationRef\"\u003e2015\u003c/span\u003e). We first aimed to validate these results in our OC cell lines SKOV3, OVCAR8 and ES-2 and found that hypoxia led to significant increases in the number of tumorspheres, which were generated by growing the cells under low-attachment conditions, and in the percentage of holoclones, both of which were used as a proxy for CSCs (Supplementary Figure S7A-B). Next, we analysed the expression of stem cell markers in OC cell lines grown under hypoxic conditions and detected an increase in the mRNA levels of \u003cem\u003eNANOG\u003c/em\u003e, \u003cem\u003eCD44\u003c/em\u003e, \u003cem\u003eSOX2\u003c/em\u003e and \u003cem\u003eEPCAM\u003c/em\u003e and the percentage of cells containing the surface CSC marker CD133 (Supplementary Figure S7C-D). These results confirm that hypoxia induces a CSC-like phenotype in OC cells.\u003c/p\u003e\n\u003cp\u003eNext, we wondered whether increased \u003cem\u003ePLD2\u003c/em\u003e expression led to an increase in the CSC population in OC cells, as suggested by the chromatin accessibility and gene expression data. Therefore, we first established OC cell lines expressing ectopic \u003cem\u003ePLD2\u003c/em\u003e cDNA or depleted of \u003cem\u003ePLD2\u003c/em\u003e using a short hairpin RNA (shRNA). The expression of \u003cem\u003ePLD2\u003c/em\u003e under these conditions was assessed at the mRNA and protein levels (Fig. \u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003eA-B and Supplementary Figure S8). We observed that the overexpression of \u003cem\u003ePLD2\u003c/em\u003e in OC cells led to a significant increase in the number of clones generated by the three cell lines, while a significant decrease was detected in the OVCAR8 cells upon \u003cem\u003ePLD2\u003c/em\u003e depletion (Fig. \u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003eC), suggesting that PLD2 promotes tumor growth. To address this question, we analysed the growth of these cell lines and found that the enhanced \u003cem\u003ePLD2\u003c/em\u003e expression led to a significant increase in proliferation, while the \u003cem\u003ePLD2\u003c/em\u003e depletion generated the opposite effect with statistical significance in all cell lines (Fig. \u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003eD). This effect was further confirmed \u003cem\u003ein vivo\u003c/em\u003e by generating xenograft models of OC cells overexpressing or depleted of \u003cem\u003ePLD2\u003c/em\u003e, showing an increase or decrease in the tumor volume, respectively, 50 days after transplantation (Fig. \u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003eE). Altogether, these results indicate that \u003cem\u003ePLD2\u003c/em\u003e expression promotes tumorigenesis.\u003c/p\u003e\n\u003cp\u003eNext, we wondered whether \u003cem\u003ePLD2\u003c/em\u003e expression was related to the formation of ovarian CSCs. First, we analysed the formation of different types of colonies, including holoclones, meroclones and paraclones, which are considered stem cells, transit-amplifying cells and differentiated cells, respectively (Barrandon and Green, \u003cspan class=\"CitationRef\"\u003e1987\u003c/span\u003e). We found a significant increase in the percentage of holoclones and a significant decrease in the percentage of paraclones in the three cell lines overexpressing \u003cem\u003ePLD2\u003c/em\u003e (Fig. \u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003eF). A significant decrease in holoclone formation was also observed in the three OC cell lines upon \u003cem\u003ePLD2\u003c/em\u003e depletion. Furthermore, we measured the formation of tumorspheres under low attachment conditions in OC cells overexpressing or depleted of \u003cem\u003ePLD2\u003c/em\u003e. We found that \u003cem\u003ePLD2\u003c/em\u003e overexpression led to a significant increase in the number of tumorspheres, while \u003cem\u003ePLD2\u003c/em\u003e depletion generated the opposite effect (Fig. \u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003eG), although we did not observe changes in the size of such tumorspheres. These data indicate that \u003cem\u003ePLD2\u003c/em\u003e expression promotes the formation of ovarian CSCs.\u003c/p\u003e\n\u003cp\u003eFinally, we analysed the expression of pluripotency and CSC marker genes in our OC cell lines overexpressing or depleted of \u003cem\u003ePLD2\u003c/em\u003e. We found that the expression of\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eSOX2\u003c/em\u003e, \u003cem\u003eCD44\u003c/em\u003e and \u003cem\u003eEPCAM\u003c/em\u003e was significantly increased in the ES-2 and SKOV3 cells overexpressing \u003cem\u003ePLD2\u003c/em\u003e, while \u003cem\u003eNANOG\u003c/em\u003e was only increased in ES-2 cells (Fig. \u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003eH). In OVCAR8 cells, the expression of these genes was increased in the same trend, although in a nonsignificant manner. Then, we measured the expression levels of these pluripotency genes in the tumorspheres extracts. First, we found that \u003cem\u003ePLD2\u003c/em\u003e was highly expressed in the tumorspheres compared with that in the total extracts transfected with only the empty vector (Fig. \u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003eH), confirming that CSCs indeed have higher expression levels of \u003cem\u003ePLD2\u003c/em\u003e. The expression of stemness genes was also increased in the tumorspheres compared with that in the total extracts, as expected, while they were further upregulated in most cases upon \u003cem\u003ePLD2\u003c/em\u003e overexpression, or downregulated upon \u003cem\u003ePLD2\u003c/em\u003e depletion (Fig. \u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003eH). Altogether, these data indicate that \u003cem\u003ePLD2\u003c/em\u003e, whose expression is induced by hypoxia, is an important oncogene in OC and that its overexpression leads to increased tumor stemness.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eThe hypoxia-induced stemness of ovarian cancer cells partially depends on\u003c/strong\u003e \u003cstrong\u003ePLD2\u003c/strong\u003e \u003cstrong\u003eexpression.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThus far, we have showed that both hypoxia and \u003cem\u003ePLD2\u003c/em\u003e expression led to an increase in CSCs in OC cells (Fig. \u003cspan class=\"InternalRef\"\u003e4\u003c/span\u003e, Supplementary Figure S7) and that \u003cem\u003ePLD2\u003c/em\u003e expression was increased under hypoxic conditions in a HIF-1\u0026alpha;-dependent manner (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003e). Thus, we wondered whether both phenomena were connected and whether the hypoxia-induced increase in CSC-like features was dependent on \u003cem\u003ePLD2\u003c/em\u003e expression. To address this, we first analysed the expression levels of stemness genes by RT‒qPCR using custom TaqMan Array plates containing probes against a selection of these genes in OC cells. We observed that either hypoxia or \u003cem\u003ePLD2\u003c/em\u003e overexpression in normoxia led to the increased expression of many stemness genes, while \u003cem\u003ePLD2\u003c/em\u003e depletion largely suppressed this increase (Fig. 5A). Indeed, hierarchical clustering of the four analysed conditions showed that the samples corresponding to EV hypoxic cells and \u003cem\u003ePLD2\u003c/em\u003e-overexpressing cells clustered together, while EV normoxic cells and \u003cem\u003ePLD2\u003c/em\u003e-depleted cells clustered separately (Fig. 5B). We confirmed these result by RT‒qPCR of individual representative genes, including \u003cem\u003eSOX2\u003c/em\u003e, \u003cem\u003eNANOG\u003c/em\u003e, \u003cem\u003eCD44\u003c/em\u003e and \u003cem\u003eEPCAM\u003c/em\u003e, showing that either hypoxia or \u003cem\u003ePLD2\u003c/em\u003e overexpression in normoxia lead to increased expression of stemness genes, while combination of both conditions further increased their expression (Fig. 5C). \u003cem\u003ePLD2\u003c/em\u003e depletion partially suppressed the hypoxia-induced enhancement of expression, with a lower non-statistically significant effect in normoxia, and rescue experiments confirmed the specificity of \u003cem\u003ePLD2\u003c/em\u003e depletion (Fig. 5C).\u003c/p\u003e\n\u003cp\u003eNext, we analyzed the formation of tumorspheres in normoxia and hypoxia with altered \u003cem\u003ePLD2\u003c/em\u003e expression. We observed that either \u003cem\u003ePLD2\u003c/em\u003e overexpression or hypoxia led to a similar increase in the formation of tumorspheres, with only a slightly higher nonsignificant increase when both conditions were combined (Fig. 5D). However, \u003cem\u003ePLD2\u003c/em\u003e depletion caused a partial suppression of the increase in tumorsphere formation under hypoxic conditions in SKOV3 and OVCAR8 cells, which was rescued by overexpressing back \u003cem\u003ePLD2\u003c/em\u003e in shRNA-transfected cells (Fig. 5D and Supplementary Figure S9A), suggesting that PLD2 is partially responsible for the hypoxia-induced stemness. Then, we measured the formation of holoclones, meroclones and paraclones under hypoxic conditions upon \u003cem\u003ePLD2\u003c/em\u003e overexpression or depletion and found that the increase in the percentage of holoclones induced by hypoxia was further enhanced by the \u003cem\u003ePLD2\u003c/em\u003e overexpression, while it was suppressed by \u003cem\u003ePLD2\u003c/em\u003e depletion and rescued back by expressing \u003cem\u003ePLD2\u003c/em\u003e after its depletion (Supplementary Figure S9B).\u003c/p\u003e\n\u003cp\u003eThen, we measured the protein levels of the pluripotency factors Sox2, Sox17, Sox9 and Notch1 (found to correlate with PLD2 in OC patients, Fig. \u003cspan class=\"InternalRef\"\u003e3\u003c/span\u003eB) by immunofluorescence in tumorspheres to determine whether PLD2 could influence their expression in CSCs. We validated PLD2 protein levels in tumorspheres (Supplementary Figure S9C) and found that either hypoxia or \u003cem\u003ePLD2\u003c/em\u003e overexpression led to an increase in the levels of Sox2, Sox9 and Notch1, while only hypoxia led to an increase in Sox17 protein levels (Fig. 5E). In addition, \u003cem\u003ePLD2\u003c/em\u003e depletion led to a partial suppression of the hypoxia-induced expression of Sox2, Sox9 and Notch1 that was rescued by expressing back \u003cem\u003ePLD2\u003c/em\u003e in these cells (Fig. 5E), suggesting that PLD2 plays a role in the generation of CSCs in hypoxia through these genes. These observations were confirmed at the mRNA level by RT‒qPCR (Fig. 5F). Altogether, these data indicate that PLD2 plays a major role in the induction of the CSC phenotype in hypoxia, promoting the expression of specific stem-related genes, such as \u003cem\u003eSOX2\u003c/em\u003e, \u003cem\u003eSOX9\u003c/em\u003e or \u003cem\u003eNOTCH1\u003c/em\u003e.\u003c/p\u003e\n\u003cp\u003eFinally, we extended the gene expression analyses to EMT genes using TaqMan Arrays to assess whether \u003cem\u003ePLD2\u003c/em\u003e expression may have a role in tumor invasion and metastasis. We found that either hypoxia or \u003cem\u003ePLD2\u003c/em\u003e overexpression in normoxia led to an increase in the expression of many of these genes, but \u003cem\u003ePLD2\u003c/em\u003e depletion was unable to suppress such increase (Supplementary Figure S9A). Indeed, hierarchical clustering of the samples did not exhibit the pattern observed in stemness genes (Supplementary Figure S9B), and results were further validated by RT‒qPCR of particular EMT genes (Supplementary Figure S9C). Consistently, invasiveness assays using Boyden\u0026rsquo;s chamber showed that both \u003cem\u003ePLD2\u003c/em\u003e overexpression and hypoxia were able to increase invasion, but \u003cem\u003ePLD2\u003c/em\u003e depletion did not have any effect (Supplementary Figure S9D). These results indicate that while the increased expression of stemness genes induced by hypoxia relies on \u003cem\u003ePLD2\u003c/em\u003e overexpression, this is not the case for EMT genes and suggests that PLD2 is a specific mediator of the increase in CSCs induced by hypoxia in OC cells.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eHypoxia-mediated reprogramming to induced pluripotent stem cells is dependent on PLD2.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe aimed to obtain additional evidence of the contribution of PLD2 to dedifferentiation or reprogramming events mediated by hypoxia that may generate ovarian CSCs from normal OC cells. Therefore, we performed reprogramming experiments of mouse embryonic fibroblasts (MEFs) to induced pluripotent stem cells (iPSCs) in normoxia and hypoxia and upon alteration of \u003cem\u003ePLD2\u003c/em\u003e expression levels (overexpression or depletion). We used a previously published protocol (Yoshida et al., \u003cspan class=\"CitationRef\"\u003e2009\u003c/span\u003e) in which MEFs were infected using a HEK293T cell-derived virus that provides OSKM genes and Nanog reporter retroviruses and then cocultured on SNL feeder cells that produce LIF. Then, the samples were incubated with or without hypoxia for 7 days, and the efficiency of iPSC generation was measured for additional 5 days. Cell reprogramming and the acquisition of pluripotency were assessed by colony morphology, alkaline phosphatase and \u003cem\u003enanog\u003c/em\u003e promoter-driven GFP expression analyses (Fig. \u003cspan class=\"InternalRef\"\u003e6\u003c/span\u003eA) to assess the effect of PLD2 and hypoxia on the efficiency of the reprogramming process and the acquisition of stem cell-like properties. Using this protocol, we found that, as expected, hypoxia led to a significant increase in the generation of iPSCs (Fig. \u003cspan class=\"InternalRef\"\u003e6\u003c/span\u003eB). Furthermore, we found that \u003cem\u003ePLD2\u003c/em\u003e overexpression in normoxia provoked a similar increase in iPSC generation, consistent with its effect on the generation of CSCs, and that the combination of both hypoxia and \u003cem\u003ePLD2\u003c/em\u003e overexpression further increased iPSC formation (Fig. \u003cspan class=\"InternalRef\"\u003e6\u003c/span\u003eB). This confirms that high \u003cem\u003ePLD2\u003c/em\u003e expression leads to dedifferentiation processes. Finally, \u003cem\u003ePLD2\u003c/em\u003e depletion in hypoxia suppressed the increased iPSC production induced by hypoxia, and this was recovered by expressing back \u003cem\u003ePLD2\u003c/em\u003e in \u003cem\u003ePLD2\u003c/em\u003e-depleted cells (Fig. \u003cspan class=\"InternalRef\"\u003e6\u003c/span\u003eB), suggesting that PLD2 is an important mediator in the activation of pluripotency by hypoxic conditions.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eOverexpression of\u003c/strong\u003e \u003cstrong\u003ePLD2\u003c/strong\u003e \u003cstrong\u003eleads to chemotherapy resistance in ovarian tumors.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eSince we showed that \u003cem\u003ePLD2\u003c/em\u003e overexpression leads to an increase in CSC-like cells in OC and CSCs were previously proposed to be responsible for chemotherapy resistance and tumor relapse, we wondered whether \u003cem\u003ePLD2\u003c/em\u003e overexpression could cause resistance to conventional therapy in ovarian tumors. Therefore, we first analysed the expression levels of \u003cem\u003ePLD2\u003c/em\u003e in our own cohort of OC patients. The immunohistochemistry analyses showed that \u003cem\u003ePLD2\u003c/em\u003e protein levels were higher in tumors than in healthy tissue (Fig. \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003eA), and RT‒qPCR revealed that \u003cem\u003ePLD2\u003c/em\u003e mRNA was significantly more abundant in OC patients than in control non-tumoral samples (Fig. \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003eB), thus confirming the results observed in the transcriptomic databases (Fig. \u003cspan class=\"InternalRef\"\u003e1\u003c/span\u003eA-B). Then, we separated our patient samples into those who were sensitive or resistant to platinum-based chemotherapy (without or with tumor relapse within the next 6 months after chemotherapy, respectively) and analysed \u003cem\u003ePLD2\u003c/em\u003e expression levels. Importantly, we found that the resistant patients showed significantly higher expression of \u003cem\u003ePLD2\u003c/em\u003e than the sensitive patients (Fig. \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003eC), suggesting that \u003cem\u003ePLD2\u003c/em\u003e overexpression may contribute to resistance to platinum-based therapy. Resistant patients in our cohort showed reduced OS and PFS (Fig. \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003eD-E), consistent with the reduced survival of patients with high \u003cem\u003ePLD2\u003c/em\u003e expression (Fig. \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003eC).\u003c/p\u003e\n\u003cp\u003eNext, we analysed the effect of \u003cem\u003ePLD2\u003c/em\u003e expression on resistance to platinum compounds in OC cells \u003cem\u003ein vitro\u003c/em\u003e. First, cells overexpressing or depleted of \u003cem\u003ePLD2\u003c/em\u003e were treated with increasing concentrations of cisplatin and carboplatin, and the IC\u003csub\u003e50\u003c/sub\u003e values were calculated in each case. We found that \u003cem\u003ePLD2\u003c/em\u003e overexpression led to a significant increase in the IC\u003csub\u003e50\u003c/sub\u003e values, while \u003cem\u003ePLD2\u003c/em\u003e depletion led to only a weak nonsignificant reduction (Fig. \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003eF). This finding suggests that higher \u003cem\u003ePLD2\u003c/em\u003e expression causes resistance to platinum-based compounds. We repeated these experiments under hypoxic conditions and found that hypoxia also led to increased IC\u003csub\u003e50\u003c/sub\u003e values (Fig. \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003eF). However, \u003cem\u003ePLD2\u003c/em\u003e depletion reduced the hypoxia-induced increase in IC\u003csub\u003e50\u003c/sub\u003e values, suggesting that enhanced resistance to cisplatin and carboplatin in OC cells under hypoxic conditions relies on \u003cem\u003ePLD2\u003c/em\u003e expression.\u003c/p\u003e\n\u003cp\u003eFinally, we performed \u003cem\u003ein vivo\u003c/em\u003e analyses to validate our findings by establishing xenograft models from SKOV3 and OVCAR8 cells overexpressing or depleted of \u003cem\u003ePLD2\u003c/em\u003e and analysing tumor growth upon treatment with cisplatin. The control tumors from cells transfected with the empty vector were sensitive to the cisplatin treatment, significantly reducing tumor growth in xenografts from both SKOV3 and OVCAR8 cells (Fig. \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003eH). However, tumors overexpressing \u003cem\u003ePLD2\u003c/em\u003e showed higher tumor growth that was not reduced upon cisplatin treatment, indicating both a higher aggressiveness of \u003cem\u003ePLD2\u003c/em\u003e-overexpressing tumors and the resistance of these tumors to cisplatin. However, \u003cem\u003ePLD2\u003c/em\u003e depletion resulted in significantly reduced tumor growth that was further reduced upon cisplatin treatment (Fig. \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003eH). Importantly, cisplatin treatment in control tumors led to increased survival, while \u003cem\u003ePLD2\u003c/em\u003e-overexpressing tumors did not exhibit improved survival, and \u003cem\u003ePLD2\u003c/em\u003e-depleted tumors exhibited increased survival independent of cisplatin treatment (Fig. \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003eI). These results indicate that the overexpression of \u003cem\u003ePLD2\u003c/em\u003e causes resistance to platinum-based chemotherapy in OC tumors.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCombination treatment with cisplatin and a PLD inhibitor suppresses chemotherapy resistance in ovarian cancer.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eFinally, we wondered whether the increased therapy resistance to platinum-based compounds induced by \u003cem\u003ePLD2\u003c/em\u003e overexpression and hypoxia could be suppressed by the pharmacological inhibition of PLD2. For this, we used the PLD inhibitor (PLDi) 5-Fluoro-2-indolyl des-chlorohalopemide (FIPI), which inhibits the catalytic activity of phospholipases D (Ganesan et al., \u003cspan class=\"CitationRef\"\u003e2015\u003c/span\u003e). First, we tested this possibility \u003cem\u003ein vitro\u003c/em\u003e by calculating the IC\u003csub\u003e50\u003c/sub\u003e in OC cells treated with cisplatin, PLDi and their combination in normoxia and hypoxia with altered levels of \u003cem\u003ePLD2\u003c/em\u003e. We found that the higher IC\u003csub\u003e50\u003c/sub\u003e to cisplatin in OC cells overexpressing \u003cem\u003ePLD2\u003c/em\u003e or in hypoxia was suppressed by the PLDi (Fig. \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003eG). Next, we validated these results \u003cem\u003ein vivo\u003c/em\u003e by establishing xenografts of OC cells expressing EV or \u003cem\u003ePLD2\u003c/em\u003e and treating mice with cisplatin, PLDi or their combination. We found that the increased tumor growth provoked by \u003cem\u003ePLD2\u003c/em\u003e overexpression was reduced upon treatment with PLDi (Fig. \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003eH). Moreover, although treatment with cisplatin did not reduce the higher tumor growth induced by \u003cem\u003ePLD2\u003c/em\u003e overexpression, its combination with PLDi led to a significant reduction in tumor growth that was stronger than that following treatment with PLDi alone (Fig. \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003eH). This finding was confirmed in xenografts from two OC cell lines (OVCAR8 and SKOV3) and led to an increase in survival (Fig. \u003cspan class=\"InternalRef\"\u003e7\u003c/span\u003eI). Altogether, these results indicate that chemotherapy resistance to cisplatin caused by \u003cem\u003ePLD2\u003c/em\u003e overexpression can be overcome by the pharmacological inhibition of PLD2, suggesting that combined treatment with cisplatin and PLDi is a promising alternative treatment for patients with high \u003cem\u003ePLD2\u003c/em\u003e expression levels.\u003c/p\u003e"},{"header":"DISCUSSION","content":"\u003cp\u003eWe show here that the HIF-1α-PLD2 axis is a major player in the chemoresistance of OC by promoting the generation of CSCs under hypoxic conditions. On the one hand, \u003cem\u003ePLD2\u003c/em\u003e expression is regulated by HIF-1α through HREs in its promoter and a hypoxia-specific enhancer; on the other hand, \u003cem\u003ePLD2\u003c/em\u003e expression is required for the full transcriptional and epigenomic rewiring promoted by hypoxia as evidenced using gene expression and chromatin accessibility analyses. In particular, the expression of stem-related genes, the opening of enhancers in the vicinity of these genes, the generation of CSCs induced by hypoxia, and the reprogramming of normal cells to iPSCs rely on normal \u003cem\u003ePLD2\u003c/em\u003e expression levels. These findings indicate that PLD2 is a major player in the response to hypoxia in cancer cells that leads to increased stem cell properties resulting in higher chemoresistance.\u003c/p\u003e \u003cp\u003eHypoxia is a feature of the tumor microenvironment in regions with low oxygen supply that is known to increase the stemness features of cancer cells in several types of cancer (Carnero and Lleonart, \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Li et al., 2009b; Mathieu et al., \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2011\u003c/span\u003e; Peng and Liu, \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e2015\u003c/span\u003e; Schwab et al., \u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Wang et al., \u003cspan citationid=\"CR71\" class=\"CitationRef\"\u003e2011\u003c/span\u003e), including OC in which hypoxia has been shown to increase the stem-like properties of cancer cells (Liang et al., \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). HIF-ARNT can regulate the expression of many genes that promote the hypoxic response (Kaelin and Ratcliffe, \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e2008\u003c/span\u003e; Schito and Semenza, \u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). Therefore, HIF factors may exert their effect of increasing CSCs via multiple mechanisms. In OC, HIF factors contribute to the upregulation of pluripotency factor genes, such as \u003cem\u003eSOX2\u003c/em\u003e or \u003cem\u003eOCT3/4\u003c/em\u003e (Liang et al., \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Seo et al., \u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e2016\u003c/span\u003e), proliferation pathways, such as Notch or Wnt (Chau et al., \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2013\u003c/span\u003e; Seo et al., \u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e2016\u003c/span\u003e), or epigenetic modulation by affecting chromatin modifiers, such as SIRT1 (Qin et al., \u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Here we show that the OC cell lines ES-2, SKOV3 and OVCAR8 under hypoxia show upregulated expression of \u003cem\u003ePLD2\u003c/em\u003e, which encodes phospholipase D2, in these cells in a Hif-1α-dependent manner (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e), as recently reported in colon cancer (Liu et al., \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Using DNA binding motif analyses and reporter assays, we also demonstrate that \u003cem\u003ePLD2\u003c/em\u003e expression is regulated at the transcriptional level by HREs located within the \u003cem\u003ePLD2\u003c/em\u003e promoter and a hypoxia-specific enhancer that activates \u003cem\u003ePLD2\u003c/em\u003e transcription in hypoxia, thus providing a mechanistic explanation of Hif-1α-mediated \u003cem\u003ePLD2\u003c/em\u003e overexpression. Therefore, it is likely that the HIf-1α-PLD2 axis works in other solid tumors under hypoxic conditions.\u003c/p\u003e \u003cp\u003eThe expression of \u003cem\u003ePLD2\u003c/em\u003e is elevated in several cancer types (Carnero et al., \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e1994\u003c/span\u003e; Frankel et al., \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e1999\u003c/span\u003e; Henkels et al., \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2013\u003c/span\u003e; Munoz-Galvan et al., \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2019b\u003c/span\u003e; Saito et al., \u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Song et al., \u003cspan citationid=\"CR64\" class=\"CitationRef\"\u003e1991\u003c/span\u003e; Zheng et al., \u003cspan citationid=\"CR75\" class=\"CitationRef\"\u003e2006\u003c/span\u003e), and here, using public patient databases and our cohort of patients, we show that this is also the case in OC in which it may be related to decreased OS (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). Furthermore, consistent with the results in OC cell lines, clustering analyses of OC patient gene expression revealed that \u003cem\u003ePLD2\u003c/em\u003e expression is correlated with the rewiring of transcriptomic programs of the response to hypoxia and stem cell maintenance (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). Indeed, we found two different clusters of patients, one of which showed gene expression patterns more different from healthy controls coinciding with higher \u003cem\u003ePLD2\u003c/em\u003e expression. Although the expression of the OC marker \u003cem\u003ePAX8\u003c/em\u003e and other stemness- and hypoxia-related genes was enhanced in both clusters, other markers showed increased expression only in the patient cluster with high \u003cem\u003ePLD2\u003c/em\u003e expression, suggesting that PLD2 may promote stemness through specific genes or pathways as we show in \u003cem\u003eSOX2\u003c/em\u003e, \u003cem\u003eSOX9\u003c/em\u003e and \u003cem\u003eNOTCH1\u003c/em\u003e, but not \u003cem\u003eSOX17\u003c/em\u003e. These data indicate that a correlation exists between \u003cem\u003ePLD2\u003c/em\u003e expression and highly altered transcriptomic programs of the response to hypoxia and stemness. In addition, we provide evidence that PLD2 is important for the alteration in the epigenomic landscape provoked by hypoxia since both hypoxia and \u003cem\u003ePLD2\u003c/em\u003e overexpression in normoxia lead to similar alterations in chromatin accessibility that are counteracted by \u003cem\u003ePLD2\u003c/em\u003e depletion, including the opening of enhancers in the proximity of genes related to the stem fate, which we show were upregulated (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). These changes connect hypoxia with stemness gene expression and likely occur through the activation of CREs bound by the AP-1 family of TFs. In a previous study analysing OC tumors, solid metastasis and effusions, higher \u003cem\u003ePLD2\u003c/em\u003e expression was described in effusions rather than solid tumors and metastasis (Harel-Dassa et al., \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). This finding is consistent with our findings since OC effusions, most of which are peritoneal, are characterized by low oxygen levels and a high content in CSCs (Munoz-Galvan and Carnero, \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Therefore, the hypoxia-PLD2 axis seems to play a major role in the generation of ovarian CSCs.\u003c/p\u003e \u003cp\u003eImportantly, \u003cem\u003ePLD2\u003c/em\u003e depletion partially suppresses the effect induced by hypoxia (Fig.\u0026nbsp;5), suggesting that PLD2 is an important mediator of hypoxia-induced stemness. This finding was also corroborated by reprogramming experiments of MEFs to iPSCs in which PLD2 was required for the increased reprogramming induced by hypoxia. Indeed, PLD2 might also potentiate the response to hypoxia by a positive feedback loop. This hypothesis can be supported since the combination of hypoxia and \u003cem\u003ePLD2\u003c/em\u003e overexpression additively enhanced the stem properties of OC cells. The activation of HIF-1α expression or activity by PLD2 has been reported in endothelial, glioma and renal cancer cells (Ghim et al., \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Han et al., \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Toschi et al., \u003cspan citationid=\"CR67\" class=\"CitationRef\"\u003e2008\u003c/span\u003e), but there is evidence of the opposite effect in HEK293 cells (Park et al., \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2015\u003c/span\u003e), suggesting that this feedback loop may work in specific contexts or cell types. Nevertheless, whether the effect of PLD2 on stemness is an autocrine or a paracrine effect remains to be elucidated. In this regard, we previously showed that exosomes from \u003cem\u003ePLD2\u003c/em\u003e-overexpressing colorectal cancer cells induced senescence in stromal fibroblasts (Munoz-Galvan et al., \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2019b\u003c/span\u003e), which, in turn, induced WNT pathway activation and increased stemness in tumor cells. In OC, hypoxia-induced exosomes have been involved in increased tumorigenic properties and chemoresistance by several mechanisms. These include exosome-containing oncogenic proteins, such as STAT3 and FAS (Dorayappan et al., \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2018\u003c/span\u003e), microRNAs that altered tumor-associated macrophages (Chen et al., \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e2018b\u003c/span\u003e; Zhu et al., \u003cspan citationid=\"CR76\" class=\"CitationRef\"\u003e2019\u003c/span\u003e), as well as plasma gelsolin, which induces the conversion of chemosensitive OC cells to chemoresistant cells (Asare-Werehene et al., \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Finally, a role of PLD1 and PLD2 inducing exosome secretion in OC cells has also been recently proposed (Onallah et al., \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e2022\u003c/span\u003e), suggesting that PLD2 may also influence the tumor microenvironment in OC.\u003c/p\u003e \u003cp\u003eSeveral mechanisms have been described to promote chemotherapy resistance under hypoxia in OC, including the upregulation of the \u003cem\u003eABCG2\u003c/em\u003e transporter gene, which increases drug efflux (He et al., \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2019\u003c/span\u003e; Wang et al., \u003cspan citationid=\"CR70\" class=\"CitationRef\"\u003e2016\u003c/span\u003e), c-\u003cem\u003eKIT\u003c/em\u003e overexpression (Chau et al., \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2013\u003c/span\u003e) and high cysteine levels (Nunes et al., \u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). CSCs are responsible for chemotherapy resistance (Batlle and Clevers, \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Beck and Blanpain, \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2013\u003c/span\u003e), and ovarian CSCs were identified sixteen years ago and reported to be chemoresistant (Bapat et al., \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2005\u003c/span\u003e; Hu et al., \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). In agreement with this idea, we previously found several markers linking ovarian CSCs and chemoresistance (Munoz-Galvan et al., \u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e2019a\u003c/span\u003e; Munoz-Galvan et al., \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Here, using OC cells and xenograft models, we show that the overexpression of \u003cem\u003ePLD2\u003c/em\u003e leads to resistance to platinum-derived compounds, including cisplatin and carboplatin (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e7\u003c/span\u003e). Moreover, \u003cem\u003ePLD2\u003c/em\u003e expression is higher in patients resistant to platinum-based chemotherapy than in sensitive patients, confirming our results in cell lines and mouse models. How PLD2 provokes such resistance is an intriguing issue, although it is likely that its enzymatic product PA, an important molecule acting as a second messenger in multiple cellular functions (Jang et al., \u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e2012\u003c/span\u003e), might play some relevant role in avoiding chemotherapy-induced cell damage. Indeed, we previously showed that PA administration has similar effects as \u003cem\u003ePLD2\u003c/em\u003e overexpression (Munoz-Galvan et al., \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2019b\u003c/span\u003e). This is particularly relevant in ovarian tumors, which are typically treated with platinum-based compounds and show high rates of chemoresistance, with frequent metastasis in hypoxic environments, such as abdominal ascites.\u003c/p\u003e"},{"header":"CONCLUSIONS","content":"\u003cp\u003eOur findings suggest a model in which hypoxia leads to the transcriptional overexpression of \u003cem\u003ePLD2\u003c/em\u003e in OC, which, in turn, generates PA and induces the generation of chemoresistant ovarian CSCs. Therefore, we propose an alternative treatment based on a combination of cisplatin and the pharmacological inhibition of PLD2. Our \u003cem\u003ein vitro\u003c/em\u003e and \u003cem\u003ein vivo\u003c/em\u003e results demonstrate that this combined treatment may be useful for patients with high \u003cem\u003ePLD2\u003c/em\u003e expression, who are resistant to conventional therapy with cisplatin alone.\u003c/p\u003e"},{"header":"METHODS","content":"\u003cp\u003e \u003cb\u003eCell culture.\u003c/b\u003e Cells were cultured according to the manufacturer's recommended procedures in McCoy (ES-2 line) or RPMI (SKOV3 and OVCAR8 lines) and incubated at 37\u0026deg;C in 5% CO\u003csub\u003e2\u003c/sub\u003e in a humidified atmosphere.\u003c/p\u003e \u003cp\u003e \u003cb\u003eGene transfer.\u003c/b\u003e The gene transfer was performed as previously described (Munoz-Galvan et al., \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). The PLD2 overexpression plasmid was described in (Munoz-Galvan et al., \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2019b\u003c/span\u003e). The shRNA against \u003cem\u003ePLD2\u003c/em\u003e was provided by Origene.\u003c/p\u003e \u003cp\u003e \u003cb\u003eProliferation assay.\u003c/b\u003e The proliferation assay was performed as previously described (Munoz-Galvan et al., \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2020\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003cb\u003eCytotoxic assay.\u003c/b\u003e ES-2, SKOV3 or OVCAR8 cells were seeded and then treated with platinum drugs and/or the PLD inhibitor 5-Fluoro-2-indolyl des-chlorohalopemide (FIPI) at 300nM concentration 24 hours later. After 96 hours, cells were stained with 0.5% crystal violet. Then, the crystal violet was solubilized in 20% acetic acid and quantified at 595 nm absorbance to measure the cell viability.\u003c/p\u003e \u003cp\u003e\u003cb\u003eMaintenance of mouse colonies.\u003c/b\u003e All experiments involving animals received expressed approval from the IBIS/HUVR Ethical Committee for the Care and Health of Animals. The mice were maintained in the IBIS animal facility according to the facility guidelines, which are based on the Real Decreto 53/2013 and were sacrificed by CO\u003csub\u003e2\u003c/sub\u003e inhalation either using a planned procedure or as a human endpoint when the animals showed significant signs of illness.\u003c/p\u003e \u003cp\u003e \u003cb\u003eColony formation assay and clonal heterogeneity analysis\u003c/b\u003e This analysis was performed as previously described in (Munoz-Galvan et al., \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Briefly, in total, 10\u003csup\u003e3\u003c/sup\u003e cells were seeded onto 10 cm plates, and every condition was evaluated in triplicate. The medium was replaced every 3 days for 12 days, and the colonies were fixed, stained and counted. The values are expressed as the number of observed colonies among the 10\u003csup\u003e3\u003c/sup\u003e seeded cells. To analyze the clonal heterogeneity, 10\u003csup\u003e2\u003c/sup\u003e random colonies were classified in triplicate as having the following phenotypes: holoclone, meroclone and paraclone (Li et al., \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e2008\u003c/span\u003e), which are considered stem cells, transit-amplifying cells and differentiated cells, respectively (Barrandon and Green, \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e1987\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003cb\u003eSphere-forming assay\u003c/b\u003e. In total, 2x 10\u003csup\u003e3\u003c/sup\u003e cells were resuspended in 1 ml of complete MammoCultTM Basal Medium (Stemcell Tech) and seeded in ultralow attachment plates. The cultures were imaged, the tumorspheres were counted, and their diameters were quantified using CellSenseDimension software on Days 2, 3 and 4.\u003c/p\u003e \u003cp\u003e \u003cb\u003eWestern blot analyses and immunofluorescence\u003c/b\u003e were performed according to standard procedures. Information of antibodies and dilutions is shown in Supplementary Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e.\u003c/p\u003e \u003cp\u003e \u003cb\u003eRT\u0026ndash;qPCR.\u003c/b\u003e The total RNA was isolated using an RNeasy kit (Qiagen), and cDNA was generated from 1 \u0026micro;g of RNA with MultiScribe Reverse Transcriptase (Applied Biosystems). qPCR was performed using a TaqMan Assay (Applied Biosystems) with probes. The relative mRNA expression was calculated as 2\u003csup\u003e\u0026minus;∆Ct\u003c/sup\u003e relative to the \u003cem\u003eACTB\u003c/em\u003e gene. Information of probes is shown in Supplementary Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e.\u003c/p\u003e \u003cp\u003e \u003cb\u003eFluorescence-activated cell sorting\u003c/b\u003e For FACS staining, live cells were incubated with antibodies for 30 minutes at dilutions specified in the manufacturer's protocols.\u003c/p\u003e \u003cp\u003e \u003cstrong\u003eiPSCs protocol\u003c/strong\u003e \u003cp\u003eBriefly, mouse embryonic fibroblasts (MEFs) were infected with 4 different retroviral vectors encoding the Yamanaka factors (pMXs-Oct3/4, pMXs-Sox2, pMXs-Klf4, and pMXs-cMyc), an additional lentiviral vector that expresses GFP in cells where the Nanog promoter/enhancer is active (mNanog-pGreenZeo), and the corresponding plasmid overexpressing PLD2, carrying shPLD2 or carrying Ev. Seventy-two hours after the infection, MEFs were seeded on top of an SNL feeder layer in the presence of ES media, and the media was renewed every 24 h. Cell reprogramming and the acquisition of pluripotency were assessed by colony morphology, GFP expression and alkaline phosphatase activity assays to assess the effect of the gene of interest/condition of interest* on the efficiency of the reprogramming process and the acquisition of stem cell-like properties.\u003c/p\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003eATAC-seq.\u003c/b\u003e\u0026nbsp;ATAC-seq assays were performed using standard protocols (Buenrostro et al., \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2013\u003c/span\u003e; Fernandez-Minan et al., \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2016\u003c/span\u003e), with minor modifications. Briefly, 70,000 ovarian cancer cells overexpressing PLD2 carrying shPLD2 or carrying Ev growing under normoxic or hypoxic conditions were collected by centrifugation for 5 min at 500 g 4\u0026deg;C. The supernatant was removed, and the cells were washed with PBS. Then, the cells were lysed in 50 \u0026micro;l of lysis buffer (10 mM Tris-HCl pH 7.4, 10 mM NaCl, 3 mM MgCl\u003csub\u003e2\u003c/sub\u003e, 0.1% NP-40, 1x Roche Complete protease inhibitors cocktail) by pipetting up and down. The whole cell lysate was used for TAGmentation, which was centrifuged for 10 min at 500 g 4\u0026deg;C, resuspended in 50 \u0026micro;l of the Transposition Reaction containing 2.5 \u0026micro;l of Tn5 enzyme and TAGmentation Buffer (10 mM Tris-HCl pH 8.0, 5 mM MgCl2, 10% w/v dimethylformamide), and incubated for 30 min at 37\u0026deg;C. Immediately after TAGmentation, DNA was purified using a Minelute PCR Purification Kit (Qiagen) and eluted in 20 \u0026micro;l. Libraries were generated by PCR amplification using NEBNext High-Fidelity 2X PCR Master Mix (NEB). The resulting libraries were multiplexed and sequenced in a HiSeq 4000 paired-end lane, producing 100 M 49-bp paired-end reads per sample.\u003c/p\u003e \u003cp\u003e \u003cb\u003eQuantification and statistical analysis.\u003c/b\u003e All statistical analyses were performed using GraphPad Prism 4. The distribution of the quantitative variables among different study groups was assessed using parametric (Student\u0026rsquo;s \u003cem\u003et test\u003c/em\u003e) or nonparametric (Kruskal\u0026ndash;Wallis or Mann\u0026ndash;Whitney) tests as appropriate. The experiments were performed a minimum of three times and were performed in independent triplicates each time. The survival data from the patient databases were analyzed by a log-rank Mantel‒Cox statistical test.\u003c/p\u003e \u003cp\u003e \u003cb\u003eAnalyses of cancer patient databases.\u003c/b\u003e We performed meta-analyses using the R2 Genomics analysis and visualization platform (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://hgserver1.amc.nl\u003c/span\u003e\u003cspan address=\"http://hgserver1.amc.nl\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) to analyze the \u003cem\u003ePLD2\u003c/em\u003e expression levels in tumor and nontumor ovarian samples from the databases. The statistical significance of the tumor versus normal samples was assessed (P\u0026thinsp;\u0026lt;\u0026thinsp;0.05). Patient survival was analyzed using an R2 Genomics analysis and visualization platform (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://hgserver1.amc.nl\u003c/span\u003e\u003cspan address=\"http://hgserver1.amc.nl\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e), which was developed by the Department of Oncogenomics of the Academic Medical Center (Amsterdam, Netherlands). Kaplan‒Meier plots showing patient survival were generated using the databases with available survival data with the scan method, which searches for the optimum survival cut-off based on statistical analyses (log-rank test), thereby identifying the most significant expression cut-off.\u003c/p\u003e \u003cp\u003e \u003cb\u003eATAC-seq data analyses.\u003c/b\u003e ATAC-seq reads were aligned to the GRCh38 (hg38) human genome assembly using Bowtie2 2.3.5 (Langmead and Salzberg, \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2012\u003c/span\u003e), and pairs separated by more than 2 kb were removed. For ATAC-seq, the Tn5 cutting site was determined as position \u0026minus;\u0026thinsp;4 (minus strand) or +\u0026thinsp;5 (plus strand) from each read start, and this position was extended 5 bp in both directions. Reads below 150 bp were considered nucleosome free. The conversion of the SAM alignment files to BAM was performed using SAMtools 1.9 (Li et al., 2009a). The conversion of BAM to BED files and peak analyses, such as overlaps or merges, were carried out using the Bedtools 2.29.2 suite (Quinlan and Hall, \u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e2010\u003c/span\u003e)⁠. The conversion of BED to BigWig files was performed using the genomecov tool from Bedtools and the wigToBigWig utility from UCSC (Haeussler et al., \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). For ATAC-seq, peaks were called using the MACS2 2.1.1.20160309 algorithm (Zhang et al., \u003cspan citationid=\"CR74\" class=\"CitationRef\"\u003e2008\u003c/span\u003e) with an FDR\u0026thinsp;\u0026lt;\u0026thinsp;0.05 for each replicate and merged into a single pool of peaks that was used to calculate the differentially accessible sites with the DESeq2 1.18.1 package in R 3.4.3 (Love et al., \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e2014\u003c/span\u003e); a value\u0026thinsp;\u0026lt;\u0026thinsp;0.05 was set as the cut-off for statistical significance of the differential accessibility. For visualization purposes, reads were extended 100 bp for ATAC-seq.\u0026nbsp;For the data comparison, all ATAC-seq experiments used were normalized using reads falling into peaks to counteract differences in background levels between experiments and replicates (Santos-Pereira et al., \u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e2019\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eHeatmaps, average profiles and k-means clustering of the ATAC-seq data were generated using the computeMatrix and plotHeatmap tools from the Deeptools 3.5 toolkit (Ramirez et al., \u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). The TF motif enrichment was calculated using XSTREME (Grant and Bailey, \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2021\u003c/span\u003e) with the standard parameters. For the gene assignment to ATAC peaks, we used the GREAT 3.0.0 tool (Hiller et al., \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2013\u003c/span\u003e)⁠, with the basal plus extension association rule and the standard parameters (5 kb upstream, 1 kb downstream, and 1 Mb maximum extension). For the footprinting analyses, we used TOBIAS 0.12.9 (Bentsen et al., \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). First, we performed bias correction using ATACorrect and calculated the footprint scores with ScoreBigwig, both with the standard parameters. Then, we used BINDetect to determine the differential TF binding of all vertebrate motifs in the JASPAR database (Fornes et al., \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). We considered motifs with a linear fold-change\u0026thinsp;\u0026ge;\u0026thinsp;15% between conditions differentially bound. Aggregated ATAC-seq signals in the footprints were visualized using PlotAggregate. ATAC-seq data generated in this study is available through the Gene Expression Omnibus (GEO) accession number GSE210599.\u003c/p\u003e \u003cp\u003e \u003cb\u003eIn vivo\u003c/b\u003e \u003cb\u003exenograft studies.\u003c/b\u003e Tumor growth was assayed following the subcutaneous injection of 4x10\u003csup\u003e6\u003c/sup\u003e SKOV3 or OVCAR8 cells that were transfected with a plasmid carrying PLD2 or shRNA against \u003cem\u003ePLD2\u003c/em\u003e in cohorts of five nude mice each that were analyzed weekly. The tumors were measured using callipers. All mice were sacrificed once the growth experiment was completed.\u003c/p\u003e \u003cp\u003e \u003cb\u003eIn vivo\u003c/b\u003e \u003cb\u003exenograft treatment.\u003c/b\u003e Tumors were harvested when they reached 1500 mm\u003csup\u003e3\u003c/sup\u003e, cut into 2\u0026times;2\u0026times;2 mm pieces and reimplanted. Mice were randomly allocated to the drug-treated and control-treated (solvent only) groups, and once the tumor reached 20 mm\u003csup\u003e3\u003c/sup\u003e, the mice received the appropriate treatment for 4 weeks (2 doses/week). The mice were monitored daily for signs of distress and weighed twice a week. The tumor size was measured, and the size was estimated according to the following equation: tumor volume = [length x width\u003csup\u003e2\u003c/sup\u003e]/2. The experiments were terminated when the tumor reached 350 mm\u003csup\u003e3\u003c/sup\u003e or when the clinical endpoint was reached. The drugs cisplatin and carboplatin were obtained from Pharmacy HUVR and were freshly prepared and administered by intraperitoneal injection. We used higher doses in mice, assuming a 70 kg average weight for humans (125 mg/dose in humans) (Munoz-Galvan et al., \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). We administered two doses per week as follows: 3.5 mg/kg of cisplatin (equivalent to 7 mg/kg, averaging 25 g body weights of each mouse) with or without 3 mg/kg of FIPI. We did not observe signs of toxicity.\u003c/p\u003e \u003cp\u003e \u003cb\u003eIn vivo\u003c/b\u003e \u003cb\u003exenografts from tumorspheres.\u003c/b\u003e This assay involved the subcutaneous injection of 1\u0026times;10\u003csup\u003e3\u003c/sup\u003e cells grown as tumorspheres into the hind legs of 4-week-old female athymic nude mice. The animals were treated as previously described, examined twice a week, incubated for 4 more weeks, and killed, and the tumors were extracted. The tumors were measured using callipers.\u003c/p\u003e \u003cp\u003e \u003cb\u003ePatient cohort.\u003c/b\u003e The entire procedure was approved by the local ethical committee of the HUVR (CEEA O309-N-15). A cohort of paraffin-embedded tissue samples from 25 patients with ovarian cancer was obtained from the biobank of Hospital Universitario Virgen del Roc\u0026iacute;o-Instituto de Biomedicina de Sevilla (Sevilla, Spain) for the RNA expression studies and the evaluation of the correlation of the clinicopathological features (see Supplementary Table S2). The samples were obtained from biopsies of patients who were subjected to platinum treatment and who were evaluated for their response according to the RECIST criteria; normal tissue, platinum-resistant tumor samples and platinum-sensitive tumor samples were obtained. The tumor samples were sent to the pathology laboratory for diagnosis and were prepared for storage with formalin fixation and paraffin embedding. The samples were stained with haematoxylin/eosin, and RNA was extracted from the tumor tissue.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate:\u0026nbsp;\u003c/strong\u003eAll methods were performed in accordance with the relevant guidelines and regulations of the Institute for Biomedical Research of Seville (IBIS) and University Hospital Virgen del Rocio (HUVR). All animal experiments and the entire procedures of the patient cohort were performed according to the experimental protocol approved by HUVR Animals Ethics (CEI 0309-N-15).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication:\u0026nbsp;\u003c/strong\u003eWritten consent for publication was obtained from all patients involved in our study\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and materials:\u0026nbsp;\u003c/strong\u003eThe datasets used and/or analysed during the current study are available from the corresponding author upon reasonable request.\u0026nbsp;ATAC-seq data generated\u0026nbsp;in this study is available through the Gene Expression Omnibus (GEO) accession number GSE210599.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests:\u0026nbsp;\u003c/strong\u003eThe authors declare that they have no other competing financial interests.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding:\u003c/strong\u003e This research was funded by Grants RTI2018-097455-B-I00 and PID2021-122629OB-I00 funded by MCIN/AEI/10.13039/501100011033 and by \u0026ldquo;ERDF A way of making Europe\u0026rdquo;, by the \u0026ldquo;European Union\u0026rdquo;. Additional grants from CIBER de C\u0026aacute;ncer (CB16/12/00275), from Consejeria de Salud (PI-0397-2017) and Project P18-RT-2501 from 2018 competitive research projects call within the scope of PAIDI 2020\u0026mdash;80% co-financed by the European Regional Development Fund (ERDF) from the Regional Ministry of Economic Transformation, Industry, Knowledge and Universities. Junta de Andaluc\u0026iacute;a. Special thanks to the AECC (Spanish Association of Cancer Research) Founding Ref. GC16173720CARR for supporting this work. SMG was funded by a grant from the Fundación AECC. EMV-S and JMS-P are funded by postdoctoral fellowships from Junta de Andaluc\u0026iacute;a (DOC_01655 and DOC_00512, respectively).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e: The authors thank the donors and the HUVR-IBiS Biobank (Andalusian Public Health System Biobank and ISCIII-Red de Biobancos PT17/0015/0041) for the human specimens that were used in this study.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor contributions\u003c/strong\u003e: SMG and AC conceived and designed this study. SMG and EMVS performed the experiments; JMSP analysed the NGS data; PEG collected the clinical data; and SMG and AC analysed and interpreted the data and drafted and edited the manuscript. All authors revised the manuscript\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n \u003cli\u003eAsare-Werehene, M., K. Nakka, A. Reunov, C.T. Chiu, W.T. Lee, M.R. Abedini, P.W. Wang, D.B. Shieh, F.J. Dilworth, E. Carmona, T. Le, A.M. Mes-Masson, D. Burger, and B.K. Tsang. 2020. The exosome-mediated autocrine and paracrine actions of plasma gelsolin in ovarian cancer chemoresistance. \u003cem\u003eOncogene\u003c/em\u003e 39:1600-1616.\u003c/li\u003e\n \u003cli\u003eBapat, S.A., A.M. Mali, C.B. Koppikar, and N.K. Kurrey. 2005. Stem and progenitor-like cells contribute to the aggressive behavior of human epithelial ovarian cancer. \u003cem\u003eCancer Res\u003c/em\u003e 65:3025-3029.\u003c/li\u003e\n \u003cli\u003eBarrandon, Y., and H. Green. 1987. Three clonal types of keratinocyte with different capacities for multiplication. \u003cem\u003eProc Natl Acad Sci U S A\u003c/em\u003e 84:2302-2306.\u003c/li\u003e\n \u003cli\u003eBatlle, E., and H. Clevers. 2017. Cancer stem cells revisited. \u003cem\u003eNat Med\u003c/em\u003e 23:1124-1134.\u003c/li\u003e\n \u003cli\u003eBeck, B., and C. Blanpain. 2013. Unravelling cancer stem cell potential. \u003cem\u003eNat Rev Cancer\u003c/em\u003e 13:727-738.\u003c/li\u003e\n \u003cli\u003eBentsen, M., P. Goymann, H. Schultheis, K. Klee, A. Petrova, R. Wiegandt, A. Fust, J. Preussner, C. Kuenne, T. Braun, J. Kim, and M. Looso. 2020. ATAC-seq footprinting unravels kinetics of transcription factor binding during zygotic genome activation. \u003cem\u003eNat Commun\u003c/em\u003e 11:4267.\u003c/li\u003e\n \u003cli\u003eBuenrostro, J.D., P.G. Giresi, L.C. Zaba, H.Y. Chang, and W.J. Greenleaf. 2013. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. \u003cem\u003eNat Methods\u003c/em\u003e 10:1213-1218.\u003c/li\u003e\n \u003cli\u003eCarnero, A., A. Cuadrado, L. del Peso, and J.C. Lacal. 1994. Activation of type D phospholipase by serum stimulation and ras-induced transformation in NIH3T3 cells. \u003cem\u003eOncogene\u003c/em\u003e 9:1387-1395.\u003c/li\u003e\n \u003cli\u003eCarnero, A., Y. Garcia-Mayea, C. Mir, J. Lorente, I.T. Rubio, and L.L. ME. 2016. The cancer stem-cell signaling network and resistance to therapy. \u003cem\u003eCancer Treat Rev\u003c/em\u003e 49:25-36.\u003c/li\u003e\n \u003cli\u003eCarnero, A., and M. Lleonart. 2016. The hypoxic microenvironment: A determinant of cancer stem cell evolution. \u003cem\u003eBioessays\u003c/em\u003e 38 Suppl 1:S65-74.\u003c/li\u003e\n \u003cli\u003eChau, W.K., C.K. Ip, A.S. Mak, H.C. Lai, and A.S. Wong. 2013. c-Kit mediates chemoresistance and tumor-initiating capacity of ovarian cancer cells through activation of Wnt/beta-catenin-ATP-binding cassette G2 signaling. \u003cem\u003eOncogene\u003c/em\u003e 32:2767-2781.\u003c/li\u003e\n \u003cli\u003eChen, A., J. Sceneay, N. Godde, T. Kinwel, S. Ham, E.W. Thompson, P.O. Humbert, and A. Moller. 2018a. Intermittent hypoxia induces a metastatic phenotype in breast cancer. \u003cem\u003eOncogene\u003c/em\u003e 37:4214-4225.\u003c/li\u003e\n \u003cli\u003eChen, X., J. Zhou, X. Li, X. Wang, Y. Lin, and X. Wang. 2018b. Exosomes derived from hypoxic epithelial ovarian cancer cells deliver microRNAs to macrophages and elicit a tumor-promoted phenotype. \u003cem\u003eCancer Lett\u003c/em\u003e 435:80-91.\u003c/li\u003e\n \u003cli\u003eDorayappan, K.D.P., R. Wanner, J.J. Wallbillich, U. Saini, R. Zingarelli, A.A. Suarez, D.E. Cohn, and K. Selvendiran. 2018. Hypoxia-induced exosomes contribute to a more aggressive and chemoresistant ovarian cancer phenotype: a novel mechanism linking STAT3/Rab proteins. \u003cem\u003eOncogene\u003c/em\u003e 37:3806-3821.\u003c/li\u003e\n \u003cli\u003eFernandez-Minan, A., J. Bessa, J.J. Tena, and J.L. Gomez-Skarmeta. 2016. Assay for transposase-accessible chromatin and circularized chromosome conformation capture, two methods to explore the regulatory landscapes of genes in zebrafish. \u003cem\u003eMethods Cell Biol\u003c/em\u003e 135:413-430.\u003c/li\u003e\n \u003cli\u003eFornes, O., J.A. Castro-Mondragon, A. Khan, R. van der Lee, X. Zhang, P.A. Richmond, B.P. Modi, S. Correard, M. Gheorghe, D. Baranasic, W. Santana-Garcia, G. Tan, J. Cheneby, B. Ballester, F. Parcy, A. Sandelin, B. Lenhard, W.W. Wasserman, and A. Mathelier. 2020. JASPAR 2020: update of the open-access database of transcription factor binding profiles. \u003cem\u003eNucleic Acids Res\u003c/em\u003e 48:D87-D92.\u003c/li\u003e\n \u003cli\u003eFrankel, P., M. Ramos, J. Flom, S. Bychenok, T. Joseph, E. Kerkhoff, U.R. Rapp, L.A. Feig, and D.A. Foster. 1999. Ral and Rho-dependent activation of phospholipase D in v-Raf-transformed cells. \u003cem\u003eBiochem Biophys Res Commun\u003c/em\u003e 255:502-507.\u003c/li\u003e\n \u003cli\u003eGanesan, R., M. Mahankali, G. Alter, and J. Gomez-Cambronero. 2015. Two sites of action for PLD2 inhibitors: The enzyme catalytic center and an allosteric, phosphoinositide biding pocket. \u003cem\u003eBiochim Biophys Acta\u003c/em\u003e 1851:261-272.\u003c/li\u003e\n \u003cli\u003eGhim, J., J.S. Moon, C.S. Lee, J. Lee, P. Song, A. Lee, J.H. Jang, D. Kim, J.H. Yoon, Y.J. Koh, C. Chelakkot, B.J. Kang, J.M. Kim, K.L. Kim, Y.R. Yang, Y. Kim, S.H. Kim, D. Hwang, P.G. Suh, G.Y. Koh, Y.Y. Kong, and S.H. Ryu. 2014. Endothelial deletion of phospholipase D2 reduces hypoxic response and pathological angiogenesis. \u003cem\u003eArterioscler Thromb Vasc Biol\u003c/em\u003e 34:1697-1703.\u003c/li\u003e\n \u003cli\u003eGrant, C.E., and T.L. Bailey. 2021. XSTREME: Comprehensive motif analysis of biological sequence datasets. \u003cem\u003ebioRxiv\u003c/em\u003e\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eHaeussler, M., A.S. Zweig, C. Tyner, M.L. Speir, K.R. Rosenbloom, B.J. Raney, C.M. Lee, B.T. Lee, A.S. Hinrichs, J.N. Gonzalez, D. Gibson, M. Diekhans, H. Clawson, J. Casper, G.P. Barber, D. Haussler, R.M. Kuhn, and W.J. Kent. 2019. The UCSC Genome Browser database: 2019 update. \u003cem\u003eNucleic Acids Res\u003c/em\u003e 47:D853-D858.\u003c/li\u003e\n \u003cli\u003eHan, S., J. Huh, W. Kim, S. Jeong, S. Min do, and Y. Jung. 2014. Phospholipase D activates HIF-1-VEGF pathway via phosphatidic acid. \u003cem\u003eExp Mol Med\u003c/em\u003e 46:e126.\u003c/li\u003e\n \u003cli\u003eHarel-Dassa, K., S. Yedgar, C.G. Trope, B. Davidson, and R. Reich. 2017. Phospholipase D messenger RNA expression and clinical role in high-grade serous carcinoma. \u003cem\u003eHum Pathol\u003c/em\u003e 62:115-121.\u003c/li\u003e\n \u003cli\u003eHarris, A.L. 2002. Hypoxia--a key regulatory factor in tumour growth. \u003cem\u003eNat Rev Cancer\u003c/em\u003e 2:38-47.\u003c/li\u003e\n \u003cli\u003eHe, M., H. Wu, Q. Jiang, Y. Liu, L. Han, Y. Yan, B. Wei, F. Liu, X. Deng, H. Chen, L. Zhao, M. Wang, X. Wu, W. Yao, H. Zhao, J. Chen, and M. Wei. 2019. Hypoxia-inducible factor-2alpha directly promotes BCRP expression and mediates the resistance of ovarian cancer stem cells to adriamycin. \u003cem\u003eMol Oncol\u003c/em\u003e 13:403-421.\u003c/li\u003e\n \u003cli\u003eHenkels, K.M., G.P. Boivin, E.S. Dudley, S.J. Berberich, and J. Gomez-Cambronero. 2013. Phospholipase D (PLD) drives cell invasion, tumor growth and metastasis in a human breast cancer xenograph model. \u003cem\u003eOncogene\u003c/em\u003e 32:5551-5562.\u003c/li\u003e\n \u003cli\u003eHiller, M., S. Agarwal, J.H. Notwell, R. Parikh, H. Guturu, A.M. Wenger, and G. Bejerano. 2013. Computational methods to detect conserved non-genic elements in phylogenetically isolated genomes: application to zebrafish. \u003cem\u003eNucleic Acids Res\u003c/em\u003e 41:e151.\u003c/li\u003e\n \u003cli\u003eHu, C.J., A. Sataur, L. Wang, H. Chen, and M.C. Simon. 2007. The N-terminal transactivation domain confers target gene specificity of hypoxia-inducible factors HIF-1alpha and HIF-2alpha. \u003cem\u003eMol Biol Cell\u003c/em\u003e 18:4528-4542.\u003c/li\u003e\n \u003cli\u003eHu, L., C. McArthur, and R.B. Jaffe. 2010. Ovarian cancer stem-like side-population cells are tumourigenic and chemoresistant. \u003cem\u003eBr J Cancer\u003c/em\u003e 102:1276-1283.\u003c/li\u003e\n \u003cli\u003eJang, J.H., C.S. Lee, D. Hwang, and S.H. Ryu. 2012. Understanding of the roles of phospholipase D and phosphatidic acid through their binding partners. \u003cem\u003eProg Lipid Res\u003c/em\u003e 51:71-81.\u003c/li\u003e\n \u003cli\u003eJing, X., F. Yang, C. Shao, K. Wei, M. Xie, H. Shen, and Y. Shu. 2019. Role of hypoxia in cancer therapy by regulating the tumor microenvironment. \u003cem\u003eMol Cancer\u003c/em\u003e 18:157.\u003c/li\u003e\n \u003cli\u003eKaelin, W.G., Jr., and P.J. Ratcliffe. 2008. Oxygen sensing by metazoans: the central role of the HIF hydroxylase pathway. \u003cem\u003eMol Cell\u003c/em\u003e 30:393-402.\u003c/li\u003e\n \u003cli\u003eKlemba, A., L. Bodnar, H. Was, K.K. Brodaczewska, G. Wcislo, C.A. Szczylik, and C. Kieda. 2020. Hypoxia-Mediated Decrease of Ovarian Cancer Cells Reaction to Treatment: Significance for Chemo- and Immunotherapies. \u003cem\u003eInt J Mol Sci\u003c/em\u003e 21:\u003c/li\u003e\n \u003cli\u003eKunz, M., and S.M. Ibrahim. 2003. Molecular responses to hypoxia in tumor cells. \u003cem\u003eMol Cancer\u003c/em\u003e 2:23.\u003c/li\u003e\n \u003cli\u003eLangmead, B., and S.L. Salzberg. 2012. Fast gapped-read alignment with Bowtie 2. \u003cem\u003eNat Methods\u003c/em\u003e 9:357-359.\u003c/li\u003e\n \u003cli\u003eLi, H., X. Chen, T. Calhoun-Davis, K. Claypool, and D.G. Tang. 2008. PC3 human prostate carcinoma cell holoclones contain self-renewing tumor-initiating cells. \u003cem\u003eCancer Res\u003c/em\u003e 68:1820-1825.\u003c/li\u003e\n \u003cli\u003eLi, H., B. Handsaker, A. Wysoker, T. Fennell, J. Ruan, N. Homer, G. Marth, G. Abecasis, R. Durbin, and S. Genome Project Data Processing. 2009a. The Sequence Alignment/Map format and SAMtools. \u003cem\u003eBioinformatics\u003c/em\u003e 25:2078-2079.\u003c/li\u003e\n \u003cli\u003eLi, Z., S. Bao, Q. Wu, H. Wang, C. Eyler, S. Sathornsumetee, Q. Shi, Y. Cao, J. Lathia, R.E. McLendon, A.B. Hjelmeland, and J.N. Rich. 2009b. Hypoxia-inducible factors regulate tumorigenic capacity of glioma stem cells. \u003cem\u003eCancer Cell\u003c/em\u003e 15:501-513.\u003c/li\u003e\n \u003cli\u003eLiang, D., Y. Ma, J. Liu, C.G. Trope, R. Holm, J.M. Nesland, and Z. Suo. 2012. The hypoxic microenvironment upgrades stem-like properties of ovarian cancer cells. \u003cem\u003eBMC Cancer\u003c/em\u003e 12:201.\u003c/li\u003e\n \u003cli\u003eLin, S.H., and A.C. Koong. 2018. Breathing New Life Into Hypoxia-Targeted Therapies for Non-Small Cell Lung Cancer. \u003cem\u003eJ Natl Cancer Inst\u003c/em\u003e 110:\u003c/li\u003e\n \u003cli\u003eLiu, M., K. Du, B. Jiang, and X. Wu. 2020. High Expression of PhospholipaseD2 Induced by Hypoxia Promotes Proliferation of Colon Cancer Cells through Activating NF- kappa Bp65 Signaling Pathway. \u003cem\u003ePathol Oncol Res\u003c/em\u003e 26:281-290.\u003c/li\u003e\n \u003cli\u003eLove, M.I., W. Huber, and S. Anders. 2014. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. \u003cem\u003eGenome Biol\u003c/em\u003e 15:550.\u003c/li\u003e\n \u003cli\u003eMathieu, J., Z. Zhang, W. Zhou, A.J. Wang, J.M. Heddleston, C.M. Pinna, A. Hubaud, B. Stadler, M. Choi, M. Bar, M. Tewari, A. Liu, R. Vessella, R. Rostomily, D. Born, M. Horwitz, C. Ware, C.A. Blau, M.A. Cleary, J.N. Rich, and H. Ruohola-Baker. 2011. HIF induces human embryonic stem cell markers in cancer cells. \u003cem\u003eCancer Res\u003c/em\u003e 71:4640-4652.\u003c/li\u003e\n \u003cli\u003eMaugeri-Sacca, M., P. Vigneri, and R. De Maria. 2011. Cancer stem cells and chemosensitivity. \u003cem\u003eClin Cancer Res\u003c/em\u003e 17:4942-4947.\u003c/li\u003e\n \u003cli\u003eMayer, A., M. Hockel, A. Wree, C. Leo, L.C. Horn, and P. Vaupel. 2008. Lack of hypoxic response in uterine leiomyomas despite severe tissue hypoxia. \u003cem\u003eCancer Res\u003c/em\u003e 68:4719-4726.\u003c/li\u003e\n \u003cli\u003eMunoz-Galvan, S., and A. Carnero. 2020. Targeting Cancer Stem Cells to Overcome Therapy Resistance in Ovarian Cancer. \u003cem\u003eCells\u003c/em\u003e 9:\u003c/li\u003e\n \u003cli\u003eMunoz-Galvan, S., B. Felipe-Abrio, M. Garcia-Carrasco, J. Dominguez-Pinol, E. Suarez-Martinez, E.M. Verdugo-Sivianes, A. Espinosa-Sanchez, L.E. Navas, D. Otero-Albiol, J.J. Marin, M.P. Jimenez-Garcia, J.M. Garcia-Heredia, A.G. Quiroga, P. Estevez-Garcia, and A. Carnero. 2019a. New markers for human ovarian cancer that link platinum resistance to the cancer stem cell phenotype and define new therapeutic combinations and diagnostic tools. \u003cem\u003eJ Exp Clin Cancer Res\u003c/em\u003e 38:234.\u003c/li\u003e\n \u003cli\u003eMunoz-Galvan, S., B. Felipe-Abrio, E.M. Verdugo-Sivianes, M. Perez, M.P. Jimenez-Garcia, E. Suarez-Martinez, P. Estevez-Garcia, and A. Carnero. 2020. Downregulation of MYPT1 increases tumor resistance in ovarian cancer by targeting the Hippo pathway and increasing the stemness. \u003cem\u003eMol Cancer\u003c/em\u003e 19:7.\u003c/li\u003e\n \u003cli\u003eMunoz-Galvan, S., A. Lucena-Cacace, M. Perez, D. Otero-Albiol, J. Gomez-Cambronero, and A. Carnero. 2019b. Tumor cell-secreted PLD increases tumor stemness by senescence-mediated communication with microenvironment. \u003cem\u003eOncogene\u003c/em\u003e 38:1309-1323.\u003c/li\u003e\n \u003cli\u003eNunes, S.C., C. Ramos, F. Lopes-Coelho, C.O. Sequeira, F. Silva, S. Gouveia-Fernandes, A. Rodrigues, A. Guimaraes, M. Silveira, S. Abreu, V.E. Santo, C. Brito, A. Felix, S.A. Pereira, and J. Serpa. 2018. Cysteine allows ovarian cancer cells to adapt to hypoxia and to escape from carboplatin cytotoxicity. \u003cem\u003eSci Rep\u003c/em\u003e 8:9513.\u003c/li\u003e\n \u003cli\u003eOnallah, H., S.T. Mannully, B. Davidson, and R. Reich. 2022. Exosome Secretion and Epithelial-Mesenchymal Transition in Ovarian Cancer Are Regulated by Phospholipase D. \u003cem\u003eInt J Mol Sci\u003c/em\u003e 23:\u003c/li\u003e\n \u003cli\u003ePark, M.H., S.S. Bae, K.Y. Choi, and S. Min do. 2015. Phospholipase D2 promotes degradation of hypoxia-inducible factor-1alpha independent of lipase activity. \u003cem\u003eExp Mol Med\u003c/em\u003e 47:e196.\u003c/li\u003e\n \u003cli\u003ePeng, G., and Y. Liu. 2015. Hypoxia-inducible factors in cancer stem cells and inflammation. \u003cem\u003eTrends Pharmacol Sci\u003c/em\u003e 36:374-383.\u003c/li\u003e\n \u003cli\u003eQin, J., Y. Liu, Y. Lu, M. Liu, M. Li, J. Li, and L. Wu. 2017. Hypoxia-inducible factor 1 alpha promotes cancer stem cells-like properties in human ovarian cancer cells by upregulating SIRT1 expression. \u003cem\u003eSci Rep\u003c/em\u003e 7:10592.\u003c/li\u003e\n \u003cli\u003eQuinlan, A.R., and I.M. Hall. 2010. BEDTools: a flexible suite of utilities for comparing genomic features. \u003cem\u003eBioinformatics\u003c/em\u003e 26:841-842.\u003c/li\u003e\n \u003cli\u003eRamirez, F., D.P. Ryan, B. Gruning, V. Bhardwaj, F. Kilpert, A.S. Richter, S. Heyne, F. Dundar, and T. Manke. 2016. deepTools2: a next generation web server for deep-sequencing data analysis. \u003cem\u003eNucleic Acids Res\u003c/em\u003e 44:W160-165.\u003c/li\u003e\n \u003cli\u003eSaito, M., M. Iwadate, M. Higashimoto, K. Ono, Y. Takebayashi, and S. Takenoshita. 2007. Expression of phospholipase D2 in human colorectal carcinoma. \u003cem\u003eOncol Rep\u003c/em\u003e 18:1329-1334.\u003c/li\u003e\n \u003cli\u003eSantos-Pereira, J.M., L. Gallardo-Fuentes, A. Neto, R.D. Acemel, and J.J. Tena. 2019. Pioneer and repressive functions of p63 during zebrafish embryonic ectoderm specification. \u003cem\u003eNat Commun\u003c/em\u003e 10:3049.\u003c/li\u003e\n \u003cli\u003eSchito, L., and G.L. Semenza. 2016. Hypoxia-Inducible Factors: Master Regulators of Cancer Progression. \u003cem\u003eTrends Cancer\u003c/em\u003e 2:758-770.\u003c/li\u003e\n \u003cli\u003eSchwab, L.P., D.L. Peacock, D. Majumdar, J.F. Ingels, L.C. Jensen, K.D. Smith, R.C. Cushing, and T.N. Seagroves. 2012. Hypoxia-inducible factor 1alpha promotes primary tumor growth and tumor-initiating cell activity in breast cancer. \u003cem\u003eBreast Cancer Res\u003c/em\u003e 14:R6.\u003c/li\u003e\n \u003cli\u003eSeo, E.J., D.K. Kim, I.H. Jang, E.J. Choi, S.H. Shin, S.I. Lee, S.M. Kwon, K.H. Kim, D.S. Suh, and J.H. Kim. 2016. Hypoxia-NOTCH1-SOX2 signaling is important for maintaining cancer stem cells in ovarian cancer. \u003cem\u003eOncotarget\u003c/em\u003e 7:55624-55638.\u003c/li\u003e\n \u003cli\u003eSiegel, R.L., K.D. Miller, and A. Jemal. 2020. Cancer statistics, 2020. \u003cem\u003eCA Cancer J Clin\u003c/em\u003e 70:7-30.\u003c/li\u003e\n \u003cli\u003eSimiantonaki, N., M. Taxeidis, C. Jayasinghe, U. Kurzik-Dumke, and C.J. Kirkpatrick. 2008. Hypoxia-inducible factor 1 alpha expression increases during colorectal carcinogenesis and tumor progression. \u003cem\u003eBMC Cancer\u003c/em\u003e 8:320.\u003c/li\u003e\n \u003cli\u003eSong, J.G., L.M. Pfeffer, and D.A. Foster. 1991. v-Src increases diacylglycerol levels via a type D phospholipase-mediated hydrolysis of phosphatidylcholine. \u003cem\u003eMol Cell Biol\u003c/em\u003e 11:4903-4908.\u003c/li\u003e\n \u003cli\u003eSun, H.C., Z.J. Qiu, J. Liu, J. Sun, T. Jiang, K.J. Huang, M. Yao, and C. Huang. 2007. Expression of hypoxia-inducible factor-1 alpha and associated proteins in pancreatic ductal adenocarcinoma and their impact on prognosis. \u003cem\u003eInt J Oncol\u003c/em\u003e 30:1359-1367.\u003c/li\u003e\n \u003cli\u003eTorre, L.A., B. Trabert, C.E. DeSantis, K.D. Miller, G. Samimi, C.D. Runowicz, M.M. Gaudet, A. Jemal, and R.L. Siegel. 2018. Ovarian cancer statistics, 2018. \u003cem\u003eCA Cancer J Clin\u003c/em\u003e 68:284-296.\u003c/li\u003e\n \u003cli\u003eToschi, A., J. Edelstein, P. Rockwell, M. Ohh, and D.A. Foster. 2008. HIF alpha expression in VHL-deficient renal cancer cells is dependent on phospholipase D. \u003cem\u003eOncogene\u003c/em\u003e 27:2746-2753.\u003c/li\u003e\n \u003cli\u003eUeda, S., T. Saeki, A. Osaki, T. Yamane, and I. Kuji. 2017. Bevacizumab Induces Acute Hypoxia and Cancer Progression in Patients with Refractory Breast Cancer: Multimodal Functional Imaging and Multiplex Cytokine Analysis. \u003cem\u003eClin Cancer Res\u003c/em\u003e 23:5769-5778.\u003c/li\u003e\n \u003cli\u003eUnruh, A., A. Ressel, H.G. Mohamed, R.S. Johnson, R. Nadrowitz, E. Richter, D.M. Katschinski, and R.H. Wenger. 2003. The hypoxia-inducible factor-1 alpha is a negative factor for tumor therapy. \u003cem\u003eOncogene\u003c/em\u003e 22:3213-3220.\u003c/li\u003e\n \u003cli\u003eWang, W.J., H. Sui, C. Qi, Q. Li, J. Zhang, S.F. Wu, M.Z. Mei, Y.Y. Lu, Y.T. Wan, H. Chang, and P.T. Guo. 2016. Ursolic acid inhibits proliferation and reverses drug resistance of ovarian cancer stem cells by downregulating ABCG2 through suppressing the expression of hypoxia-inducible factor-1alpha in vitro. \u003cem\u003eOncol Rep\u003c/em\u003e 36:428-440.\u003c/li\u003e\n \u003cli\u003eWang, Y., Y. Liu, S.N. Malek, P. Zheng, and Y. Liu. 2011. Targeting HIF1alpha eliminates cancer stem cells in hematological malignancies. \u003cem\u003eCell Stem Cell\u003c/em\u003e 8:399-411.\u003c/li\u003e\n \u003cli\u003eWilson, W.R., and M.P. Hay. 2011. Targeting hypoxia in cancer therapy. \u003cem\u003eNat Rev Cancer\u003c/em\u003e 11:393-410.\u003c/li\u003e\n \u003cli\u003eYoshida, Y., K. Takahashi, K. Okita, T. Ichisaka, and S. Yamanaka. 2009. Hypoxia enhances the generation of induced pluripotent stem cells. \u003cem\u003eCell Stem Cell\u003c/em\u003e 5:237-241.\u003c/li\u003e\n \u003cli\u003eZhang, Y., T. Liu, C.A. Meyer, J. Eeckhoute, D.S. Johnson, B.E. Bernstein, C. Nusbaum, R.M. Myers, M. Brown, W. Li, and X.S. Liu. 2008. Model-based analysis of ChIP-Seq (MACS). \u003cem\u003eGenome Biol\u003c/em\u003e 9:R137.\u003c/li\u003e\n \u003cli\u003eZheng, Y., V. Rodrik, A. Toschi, M. Shi, L. Hui, Y. Shen, and D.A. Foster. 2006. Phospholipase D couples survival and migration signals in stress response of human cancer cells. \u003cem\u003eJ Biol Chem\u003c/em\u003e 281:15862-15868.\u003c/li\u003e\n \u003cli\u003eZhu, X., H. Shen, X. Yin, M. Yang, H. Wei, Q. Chen, F. Feng, Y. Liu, W. Xu, and Y. Li. 2019. Macrophages derived exosomes deliver miR-223 to epithelial ovarian cancer cells to elicit a chemoresistant phenotype. \u003cem\u003eJ Exp Clin Cancer Res\u003c/em\u003e 38:81.\u003cstrong\u003e\u003c/strong\u003e\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":true,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"journal-of-experimental-and-clinical-cancer-research","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"jecc","sideBox":"Learn more about [Journal of Experimental \u0026 Clinical Cancer Research](http://jeccr.biomedcentral.com)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/jecc/default.aspx","title":"Journal of Experimental \u0026 Clinical Cancer Research","twitterHandle":"@OncoBioMed","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"BMC/SO AJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"ovarian cancer, hypoxia, phospholipase D, therapy resistance, stemness","lastPublishedDoi":"10.21203/rs.3.rs-3730407/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-3730407/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cstrong\u003eBackground:\u003c/strong\u003e Hypoxia in solid tumors is an important source of chemoresistance that can determine poor patient prognosis. Such chemoresistance relies on the presence of cancer stem cells (CSCs), and hypoxia promotes their generation through transcriptional activation by HIF transcription factors.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMethods.\u003c/strong\u003e We used OC cell lines, xenograft models, OC patient samples, transcriptional databases, iPSCs and ATAC-seq\u003cstrong\u003e.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eResults\u003c/strong\u003e Here, we show that hypoxia induces CSC formation and chemoresistance in ovarian cancer through transcriptional activation of the \u003cem\u003ePLD2\u003c/em\u003e gene. Mechanistically, HIF-1a activates \u003cem\u003ePLD2\u003c/em\u003etranscription through hypoxia response elements, and both hypoxia and \u003cem\u003ePLD2\u003c/em\u003eoverexpression lead to increased accessibility around stemness genes, detected by ATAC-seq, at sites bound by AP-1 transcription factors. This in turn provokes a rewiring of stemness genes, including the overexpression of \u003cem\u003eSOX2\u003c/em\u003e, \u003cem\u003eSOX9\u003c/em\u003e or \u003cem\u003eNOTCH1\u003c/em\u003e. \u003cem\u003ePLD2\u003c/em\u003e overexpression also leads to decreased patient survival, enhanced tumor growth and CSC formation, and increased iPSCs reprograming, confirming its role in dedifferentiation to a stem-like phenotype. Importantly, hypoxia-induced stemness is dependent on \u003cem\u003ePLD2\u003c/em\u003e expression, demonstrating that PLD2 is a major determinant of de-differentiation of ovarian cancer cells to stem-like cells in hypoxic conditions. Finally, we demonstrate that high \u003cem\u003ePLD2\u003c/em\u003eexpression increases chemoresistance to cisplatin and carboplatin treatments, both \u003cem\u003ein vitro\u003c/em\u003e and \u003cem\u003ein vivo\u003c/em\u003e, while its pharmacological inhibition restores sensitivity.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConclusions. \u003c/strong\u003eAltogether, our work highlights the importance of the HIF-1a-PLD2 axis for CSC generation and chemoresistance in OC and proposes an alternative treatment for patients with high \u003cem\u003ePLD2\u003c/em\u003e expression.\u003c/p\u003e","manuscriptTitle":"Suppression of hypoxia-induced stemness and chemoresistance in ovarian tumors","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2023-12-15 17:28:12","doi":"10.21203/rs.3.rs-3730407/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Major revision","date":"2023-12-22T02:09:27+00:00","index":"","fulltext":""},{"type":"reviewerAgreed","content":"","date":"2023-12-12T08:48:13+00:00","index":0,"fulltext":""},{"type":"reviewersInvited","content":"","date":"2023-12-12T07:25:47+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2023-12-12T01:04:11+00:00","index":"","fulltext":""},{"type":"submitted","content":"Journal of Experimental \u0026 Clinical Cancer Research","date":"2023-12-11T06:05:05+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"journal-of-experimental-and-clinical-cancer-research","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"jecc","sideBox":"Learn more about [Journal of Experimental \u0026 Clinical Cancer Research](http://jeccr.biomedcentral.com)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/jecc/default.aspx","title":"Journal of Experimental \u0026 Clinical Cancer Research","twitterHandle":"@OncoBioMed","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"BMC/SO AJ","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"fc95b69e-35b2-43f6-b6d9-aa50d757d115","owner":[],"postedDate":"December 15th, 2023","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2024-03-04T15:12:41+00:00","versionOfRecord":{"articleIdentity":"rs-3730407","link":"https://doi.org/10.1186/s13046-024-02988-y","journal":{"identity":"journal-of-experimental-and-clinical-cancer-research","isVorOnly":false,"title":"Journal of Experimental \u0026 Clinical Cancer Research"},"publishedOn":"2024-02-26 15:01:38","publishedOnDateReadable":"February 26th, 2024"},"versionCreatedAt":"2023-12-15 17:28:12","video":"","vorDoi":"10.1186/s13046-024-02988-y","vorDoiUrl":"https://doi.org/10.1186/s13046-024-02988-y","workflowStages":[]},"version":"v1","identity":"rs-3730407","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-3730407","identity":"rs-3730407","version":["v1"]},"buildId":"omnImTCwR2MFx8CMYfrG7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}