Abstract
Ovarian cancer (OC) cells frequently metastasize to the omentum, and adipocytes play a significant role in ovarian tumor progression. Therapeutic interventions targeting aberrant DNA methylation in ovarian tumors have shown promise in the clinic, but the effects of epigenetic therapy on the tumor microenvironment are understudied. Here, we examined the effect of adipocytes on OC cell behavior in culture and impact of targeting DNA methylation in adipocytes on OC metastasis. The presence of adipocytes increased OC cell migration and invasion, and proximal and direct coculture of adipocytes increased OC proliferation alone or after treatment with carboplatin. Treatment of adipocytes with hypomethylating agent guadecitabine decreased migration and invasion of OC cells toward adipocytes. Subcellular protein fractionation of adipocytes treated with guadecitabine revealed decreased DNA methyltransferase 1 (DNMT1) levels even in the presence of DNA synthesis inhibitor, aphidicolin. Methyl-Capture- and RNA-sequencing analysis of guadecitabine-treated adipocytes revealed derepression of tumor-suppressor genes and epithelial–mesenchymal transition inhibitors. SUSD2, a secreted tumor suppressor downregulated by promoter CpG island methylation in adipocytes, was upregulated after guadecitabine treatment, and recombinant SUSD2 decreased OC cell migration and invasion. Integrated analysis of the methylomic and transcriptomic data identified pathways associated with inhibition of matrix metalloproteases and fatty acid α-oxidation, suggesting a possible mechanism of how epigenetic therapy of adipocytes decreases metastasis. In conclusion, the effect of DNMT inhibitor on fully differentiated adipocytes suggests that hypomethylating agents may affect the tumor microenvironment to decrease cancer cell metastasis.
Implications: Epigenetic targeting of tumor microenvironment can affect metastatic behavior of ovarian cancer cells. Mol Cancer Res; 16(8); 1226–40. ©2018 AACR.
This article is featured in Highlights of This Issue, p. 1207
Introduction
Increased adiposity is associated with increased cancer incidence, morbidity, and mortality (1). Cancer-associated adipocytes support the growth and metastasis of cancer cells, contributing to an invasive and aggressive phenotype (1). Approximately 22,440 women were diagnosed with epithelial ovarian cancer (EOC) in 2016, and the 5-year survival rate was only 28.3%, making EOC the deadliest gynecologic disease (2). EOC patients are often diagnosed at late stage with 80% of metastasis commonly found on the omentum (3), a membranous fatty tissue layer in the lower abdominal area, and a recent meta-analysis that included 25,157 women with EOC from 47 epidemiologic studies demonstrated that EOC risk is associated with higher body mass index and weight, suggesting that increased adiposity is correlated with EOC risk (4).
The degree of adipose tissue invasion is indicative of the aggressiveness of the tumor (3), and obesity contributes to OC metastatic success through increased lipogenesis and lipid transport, demonstrating the importance of the role adipocytes play in tumor progression as a rich energy source of free fatty acids (5). EOC-associated adipocytes have been shown to transfer lipids to OC cells through upregulation of FABP4 (3), and overexpression of SIK2 in metastatic ovarian tumor samples from adipocyte-rich areas promotes OC cell fatty acid oxidation and activates the PI3K/AKT pathway leading to cancer cell growth (6). Due to the supportive role adipocytes play in OC tumor growth, adipocytes might serve as an additional therapeutic target for OC metastasis.
Obesity-induced, proinflammatory cytokines promote DNMT1 expression and its enzymatic activity, leading to greater global DNA methylation and epigenetic silencing of genes in adipocytes (7, 8). The presence of cancer cells may epigenetically modify adipocytes, lead to disruption of cell signaling between adipocytes and OC cells, and promote OC metastasis. Inhibiting DNMT1 in adipocytes was previously shown to demethylate and reexpress adiponectin, an important regulator of glucose and lipid metabolism (7), providing proof of concept that epigenetic inhibition in adipocytes can result in reexpression of hypermethylated genes. Adipocytes secrete cytokines that increase OC cell homing, and once adjacent to the adipocyte, OC cells alter metabolic processes of adipocytes to undergo lipolysis and transfer of lipids (3). The epigenetic alterations that occur in adipocytes after OC cell interaction remain unexplored.
Adipose tissue responds to changing metabolic settings such as obesity by epigenetic regulation of metabolic genes (8) like adiponectin DNA methylation, which can be reversed by hypomethylating therapy (7). However, epigenetic targeting of adipocytes has not been studied in the context of cancer progression where DNA methylation alterations may have occurred due to OC cell interaction. We have demonstrated that the next-generation DNMT inhibitor guadecitabine (SGI-110) has significant activity against OC (9, 10), but effects of the hypomethylating agent on the microenvironment in general and on adipocytes in particular are not well understood. Here, we show that guadecitabine has direct effects on adipocytes. Guadecitabine treatment decreased DNMT1 levels and global demethylation, corresponding to gene expression changes in adipocytes. Furthermore, epigenetic reexpression of tumor-suppressor genes, including SUSD2, suppressed OC migration and invasion. We demonstrate that guadecitabine-treated adipocytes inhibited OC cell migration and invasion, suggesting that in addition to direct effects of DNMT inhibitors on tumor cells, epigenetic-induced alterations in the tumor microenvironment may contribute to decreased OC cell metastasis and has implications that adipocytes may be a therapeutic target for inhibiting OC cell metastasis.
Materials and Methods
Cell culture
High-grade serous (HGS) EOC cell lines OVCAR4, OVCAR5, OVCAR8, and OVCAR8-RFP were grown in DMEM media. OVCAR8 cells were transduced to express red fluorescent protein (RFP) using commercially procured lentiviral particles expressing RFP and blasticidin gene selection marker (11). Kuramochi (HGSOC) and SKOV3 (endometrioid OC) cells were grown in RPMI 1640 media and McCoy's 5A media, respectively. Normal, human primary subcutaneous preadipocytes (ATCC PCS-210-010) were grown in fibroblast basal medium (ATCC PCS-201-030) supplemented with phenol red. Simpson–Golabi–Behman syndrome (SGBS) human adipocyte cells were (provided by Dr. Martin Wabitsch) differentiated as described (12). All media were supplemented with 10% FBS, antibiotic–antimycotic, 25 mmol/L HEPES, MEM vitamins, and MEM nonessential amino acids. Cells were treated with guadecitabine (Astex Pharmaceuticals, Inc.) dissolved in DMSO or carboplatin (Sigma) dissolved in PBS. Cells were treated with 100 nmol/L guadecitabine and refreshed every day for 3 days.
Adipogenic induction of preadipocytes
Preadipocytes were differentiated following an adipogenic differentiation protocol (13). Briefly, preadipocytes were seeded at 18,000 cells/cm2. Initiation of adipogenesis was induced by replacing the fibroblast basal medium with adipogenic induction medium for 15 days. Oil Red O staining was performed to confirm adipocyte differentiation.
Western blot
Adipocyte-conditioned media (ACM) were collected and concentrated (Amicon Ultra-0.5 Centrifugal Filter Units 10 kDa or 30 kDa; Millipore) and proteins were extracted with RIPA, as described previously (9). Primary antibodies were DNMT1 (Cell Signaling Technology; #5032S), histone H3 (Cell Signaling Technology; #4499S), GAPDH (Santa Cruz Biotechnology; #sc-25778), PPARγ (Cell Signaling Technology; #2435P), and SUSD2 (Prestige Antibodies; #HPA004117). After incubation with peroxidase-labeled secondary anti-rabbit (KPL #374-1506) or anti-mouse antibodies (KPL #474-1806), protein bands were visualized using ECL reagents.
Propidium iodide cell-cycle and cell proliferation analysis using flow cytometry
Cell-cycle analysis of preadipocyte and adipocytes was performed according to Kim and colleagues (14). Cells were treated with aphidicolin, ethanol fixed, stained with propidium iodide, and analyzed on BD LSR II flow cytometer. For cell proliferation analysis, OVCAR8-RFP cells were seeded directly onto adipocytes (1:3) treated with guadecitabine. Then, OVCAR8-RFP and adipocytes were treated with carboplatin for 1 week, imaged using EVOS FL Auto Cell Imaging System, and analyzed on BD LSR II flow cytometer.
Caspase-Glo 3/7 apoptosis assay and trypan blue assay
Caspase-Glo 3/7 assay (Promega) was performed according to the manufacturer's protocol on preadipocytes or adipocytes seeded in a 96-well plate and treated with aphidicolin (1–4 μg/mL, 4X) or guadecitabine, respectively. Trypan Blue assay was performed at 0.4% trypan blue.
Subcellular fractionation
Cytoplasm and nuclear extraction protocol was adapted from O'Hagan and colleagues (15). Cell pellets were sequentially washed in CEBN buffer to obtain the cytoplasm, CEB buffer (CEBN without IGEPAL CA-630) to obtain the total nuclear fraction, and soluble nuclear buffer to obtain the chromatin and soluble nuclear fraction.
qRT-PCR
RNA was extracted with an RNeasy kit (Qiagen; #74104), cDNA was prepared using M-MLV RT system (Promega), and qPCR was performed with cDNA and primers. Data were normalized to EEF control gene. Fold change was calculated using 2−Δ (ΔCT); where ΔCT = CT, target – CT, GAPDH and Δ(ΔCT) = ΔCT, stimulated – ΔCT, control. Primer sequences were listed in Supplementary Table S1. TSH2B and GAPDH primers were control primers provided in a Diagenode MethylCap kit.
MTT cell proliferation assay
MTT assay was performed as described (9). OC cells in 96-well plates were treated with centrifuged and clarified ACM or recombinant SUSD2 (100, 500 ng/mL) with carboplatin (7 days). OC cells were seeded in the top compartment of 0.45 μm Boyden chamber assay (Corning #353095) with proximally cocultured adipocytes in the bottom of the 24-well plate in serum-free media.
Migration and invasion assay
Adipocytes seeded in 24-well plates were treated with guadecitabine (100 nmol/L, 3X daily). After drug treatment, to remove any direct effects of guadecitabine on tumor cells, old media were diluted (75% with serum-free media), and this dilution approach was performed 4X before replacing the media with 500 μL of fresh serum-free media prior to the addition of cancer cells 24 hours later. Adipocytes cannot be washed due to sensitivity to air and agitation; therefore, dilution was utilized to remove guadecitabine from the media. To test whether residual guadecitabine has any effect on cell migration and invasion, we repeated this experiment as stated above with guadecitabine in 10% FBS serum media without adipocytes. Clarified ACM (serum-free) were collected after guadecitabine treatment was removed, diluted 4X, and replaced with fresh media (24 hours) to be used as chemoattractant. Recombinant, human SUSD2 (Abnova; #H00056241-P01, 100 or 500 ng/mL) in 10% FBS DMEM was also used as a chemoattractant in the bottom of the well to test the effect of SUSD2 on migration and invasion. Cells (200,000) were seeded in migration (8 μm pore; Corning; #354578) and matrigel invasion Boyden chamber inserts (8 μm pore; Corning; #354480) with serum-free media and allowed to migrate or invade for 8 and 16 hours, respectively. The HEMA3 stain set (Protocol) was used to stain inserts. Cells were imaged using the 5X objective (5 fields per insert) and counted using ImageJ.
Bromodeoxyuridine cell proliferation assay through ELISA and fluorescence microscopy
Bromodeoxyuridine (BrdUrd) cell proliferation assay was performed (Cell Signaling Technology) according to the manufacturer's protocol on preadipocytes treated with aphidicolin daily (1–4 μg/mL, 4X). Also, adipocytes were differentiated on glass coverslips in 6-well plates and treated with aphidicolin. Cells were stained with 10 μmol/L BrdUrd (2 hours), fixed in 10% formalin, permeabilized (0.1% Triton X-100 buffer), acid-treated (HCl and phosphate/citric acid buffer), and stained (BrdUrd primary antibody 1:10,000, Cell Signaling Technology #5292S; Alexa Fluor 568 secondary antibody 1:10,000, Invitrogen #A10042). Slides were mounted using Fluoroshield (Sigma) and viewed on Nikon NiE microscope.
SUSD2 ELISA
Adipocytes were seeded in 6-well plate in serum-free DMEM:F12 and treated with guadecitabine (100 nmol/L, 3X daily); on day 4, ACM were collected, centrifuged to clarify, and used in SUSD2 ELISA assay (LifeSpan BioSciences, Inc.) according to the manufacturer's protocol.
Indirect flow cytometry of SUSD2
Adipocytes (∼1 × 106 cells) were harvested, not permeabilized, and incubated with primary SUSD2 antibody (1:1,000; Prestige Antibodies; #HPA004117) and 10% FBS overnight at 4°C. Cells were washed with PBS (3X), incubated with secondary Alexa Fluor 488 antibody (1:5,000; Invitrogen; #A11034) for 1 hour at room temperature in the dark, washed with PBS (3X), suspended in 0.5 mL PBS, 3% BSA, filtered, and analyzed on BD LSR II flow cytometer.
Next-generation sequencing analyses
For Methyl-Capture-sequencing (MC-seq) library preparation, DNA was extracted with an AllPrep kit (Qiagen; #80204) from adipocytes, sonicated (180–200 bp), and used as input for a methylated DNA capture kit (Diagenode) according to the manufacturer's protocol. MC- and RNA-seq libraries were prepared with a NEXTflex Methyl-Seq1 kit (Bioo Scientific) or a TruSeq Stranded mRNA Library Prep kit (Illumina #RS-122-2101), respectively, according to the manufacturer's protocol. Libraries were verified using Qubit and Bioanalyzer, prior to sequencing on Illumina NextSeq. NextSeq read sequences were cleaned using Trimmomatic version 0.32 with the parameters “2:20:5 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 MINLEN:17″ to remove adapter sequences and perform quality trimming. Custom perl scripts and FASTX-Toolkit version 0.0.13.2 (-l 17 -v -M 1 –i) were used to identify partial adapters.
For MC-seq data, reads were mapped against GRCh38.p5 using Bowtie version 2.2.9 (-k 20). A custom perl script was used to obtain only the mapping results with the best alignment score. Initial differential methylation analysis was carried out using the MethylAction package (version 2.0) in R (version 3.3.1) and Gencode v24 as the annotation (16). Window sizes of 50 base pairs were used in the getCounts function of the package. The Goldmine package (version 1.0) in R was used to assign biological annotation to the joined window results produced by MethylAction (16).
For RNA-seq data, the resulting reads were mapped against GRCh38.p5 using TopHat2 version 2.1.1 with parameters “–b2-very-sensitive –read-edit-dist 2 –max-multihits 100.″ Read counts for each gene were created using custom perl scripts and Gencode v24 as the annotation. Custom perl scripts were used for estimation of transcript abundances based on Fragments Per Kilobase of exon per Million fragments mapped. Initial differential expression analysis was carried out using the DESeq2 package (version 1.12.4) in R/Bioconductor (R version 3.3.1; ref. 17). Pathway analysis was performed on both MC- and RNA-seq data using Ingenuity Pathway Analysis (IPA; Qiagen).
Statistical analysis
All data are presented as mean values ± SD of at least three biological experiments unless otherwise indicated. IC50 values were determined by Prism 6 (GraphPad Software), using logarithm-normalized sigmoidal dose curve fitting. The estimate variation within each group was similar; therefore, the Student t test was used to statistically analyze the significant difference among different groups by using Prism 4.0 (GraphPad Software).
Results
Epigenetic inhibition of adipocytes decreases OC cell migration and invasion
Adipocytes of the omentum play a key role in promoting the growth of metastatic lesions (18). To examine the effects of adipocytes on migration and invasion, EOC cells were seeded in the top compartment of the Boyden chamber and allowed to migrate (8 hours) and invade (16 hours) toward adipocytes in the bottom of the well. The presence of adipocytes increased (P < 0.05) migration (>20-fold; Fig. 1A) and invasion (>10-fold; Fig. 1B) of all EOC cells lines examined (Supplementary Fig. S1A and S1B). Based on previous evidence of obesity-induced methylation changes in adipose tissue (8) and altered methylation in adipocytes by DNMT inhibition (7), we hypothesized that epigenetic inhibition of adipocytes may demethylate and reexpress factors that could affect migration and invasion of OC cells. Guadecitabine pretreatment of adipocytes (100 nmol/L, 3X daily) decreased (P < 0.05) EOC cell migration (>26%; Fig. 1C) and invasion (>27%; Fig. 1D) toward adipocytes with the exception of SKOV3 (Supplementary Fig. S1C and S1D). To test whether residual guadecitabine has any effect on cell migration and invasion, we repeated this experiment as before with guadecitabine (100 nmol/L, 3X daily) in 10% FBS serum media without adipocytes. There was no change in migration and invasion in the presence of residual guadecitabine (Supplementary Fig. S2A and S2B), suggesting that decreased cell migration and invasion toward adipocytes were a result of guadecitabine treatment of adipocytes and not residual drug concentration. To assess whether guadecitabine has an effect on adipocyte viability, adipocytes were treated with guadecitabine (100 nmol/L, daily 3X), and caspase assay was performed. No effect on caspase levels in adipocytes was observed after guadecitabine treatment (Supplementary Fig. S2G), demonstrating that a change in adipocyte population was not responsible for the decreased migration and invasion after guadecitabine treatment.
Adipocyte cancer cell's two-way communication occurs through secreted factors in the Boyden chamber assay. To assess if similar results could be obtained using ACM, migration and invasion assays were repeated with Kuramochi and OVCAR4 cells. A decrease (P < 0.001) in migration (Supplementary Fig. S2C and S2D) and invasion (Supplementary Fig. S2E and S2F) toward ACM compared with adipocytes was observed, and no effect of ACM from guadecitabine-treated adipocytes on either migration or invasion was seen (Supplementary Fig. S2C–S2F). One explanation for these observations is that cytokines present in ACM may be unstable, and without constant cytokine secretion by adipocytes, migration and invasion were not induced. However, as initial Boyden chamber assays of adipocytes and OC cells were indirectly cocultured, secreted factors were likely responsible for the observed effects.
Decreased EOC cell migration and invasion after guadecitabine treatment of adipocytes suggested an alteration of secreted factors into the media. ACM collected from adipocytes treated with or without guadecitabine was used to culture OVCAR8 and OVCAR4 cells, and expressions of known epithelial–mesenchymal transition (EMT) genes (SLUG, FN1, TWIST1, CDH1, CDH2, VIM, and ITGB1) were examined relative to basal EMT marker levels in OVCAR4/8 cells grown in regular media. Guadecitabine ACM increased (P < 0.01) CDH1 expression in OVCAR8 (Fig. 1E) and decreased (P < 0.05) expression of SLUG, FN1, and CDH2 in OVCAR4 (Fig. 1F), suggesting that guadecitabine altered secretion of adipokine(s) which modified EMT marker expression in EOC.
Guadecitabine induces decreased DNMT1 levels in adipocytes
The known molecular mechanism of guadecitabine is incorporation of decitabine (5-aza-CdR) into replicating DNA during cell division which covalently binds DNMT on carbon 6 of decitabine and leads to DNMT depletion (19, 20). However, as mature adipocytes do not divide (21), and the cultures represented a heterogeneous mixture of mature and preadipocytes, we tested whether undifferentiated adipocytes may have contributed to the observed response to guadecitabine. Expression of the preadipocyte marker DLK1 (22) was low (P < 0.001) in adipocytes (Supplementary Fig. S3A) compared with mature adipocyte marker PPARγ (Supplementary Fig. S3B). Furthermore, after adipocyte differentiation (i.e., mature), increased (P < 0.05) expression of FABP4 was observed (Supplementary Fig. S3A). Because of the known mechanism of action of guadecitabine and the absence of preadipocyte markers in the differentiated cultures, it was of interest to determine whether or not any cells in the differentiated culture were undergoing DNA synthesis. Optimal concentration of aphidicolin to inhibit DNA synthesis (1 μg/mL aphidicolin, daily 4X) was determined by BrdUrd ELISA (Supplementary Fig. S3C). No cell toxicity was detected by caspase 3/7 (Supplementary Fig. S3D) and trypan blue assays (Supplementary Fig. S3E). BrdUrd staining in preadipocytes was reduced after aphidicolin treatment (P < 0.05), whereas aphidicolin treatment of adipocytes resulted in negative BrdUrd staining similar to control (Fig. 2A). In addition, aphidicolin treatment decreased (P < 0.05) preadipocyte cell proliferation but had no effect on adipocyte proliferation (Fig. 2B). Furthermore, cell-cycle analysis revealed that 6.78% of preadipocytes were in S phase, whereas only 1.80% of adipocytes were in S phase with no difference after aphidicolin treatment (Fig. 2C).
To test whether guadecitabine depleted soluble DNMT1 and/or chromatin-bound DNMT1 (23), subcellular fractionation of adipocytes was performed. DNMT1 was not detectable in cytoplasm; however, chromatin-bound DNMT1 (Fig. 2D) and whole-cell DNMT1 (Supplementary Fig. S3F) were decreased after guadecitabine treatment, despite the presence of DNA synthesis inhibitor. GAPDH and histone H3 serve as cytoplasmic and nuclear controls, respectively, for the cellular fractionation. Furthermore, DNMT1 mRNA expression in adipocytes was unchanged after guadecitabine treatment (Supplementary Fig. S3G), suggesting that the molecular mechanism of guadecitabine in adipocytes is replication-independent because DNMT1 mRNA levels positively correlated with replication. Taken together, these results demonstrated guadecitabine-induced decreased DNMT1 levels in adipocytes in the absence of cell division or DNA synthesis.
Adipocytes increase HGSOC cell proliferation
Previous studies indicated a role of cancer associated adipocytes in proliferation and carboplatin resistance (24) and that hypomethylating agents can reverse platinum resistance in EOC cells (9, 25). Therefore, we investigated if adipocytes increased OC cell proliferation and whether DNMT inhibition of adipocytes could reverse proliferation of EOC cells. Unexpectedly, ACM alone had no effect on HGSOC cell proliferation (Supplementary Fig. S4A–S4D); however, indirect coculture of HGSOC cells with adipocytes increased (P < 0.05) OC cell proliferation, and pretreatment with guadecitabine resulted in no decrease (Fig. 3A–C) except in Kuramochi cells (Fig. 3D). To test the effects of direct coculture, adipocytes were treated with guadecitabine, the DNMT inhibitor was removed, OVCAR8-RFP cells were then seeded directly onto adipocytes, and both cells were treated with increasing doses of carboplatin. The presence of adipocytes increased (P < 0.05) OVCAR8-RFP proliferation in a carboplatin dose-dependent manner (Fig. 3F) compared with control (Fig. 3E), and guadecitabine treatment of adipocytes resulted in no change in OVCAR8-RFP proliferation (Fig. 3G), compared with carboplatin alone. OVCAR8-RFP proliferation with adipocytes was quantified using flow cytometry (Fig. 3H). These results demonstrated adipocytes increased OVCAR8-RFP proliferation when treated with carboplatin in direct coculture with adipocytes.
Guadecitabine induces DNA demethylation and gene expression in adipocytes
To assay global methylation changes, Methyl-Capture-seq (MC-seq) was performed on adipocytes after low-dose (100 nmol/L, daily 3X) guadecitabine treatment. The methyl-capture DNA library was validated prior to sequencing by qRT-PCR of TSH2B (hypermethylated gene control) and GAPDH (hypomethylated gene control; Supplementary Fig. S5A). MethylAction (16) analysis resulted in 3,246 hypomethylated and 309 hypermethylated genes compared with control (FC > 10; FDR < 0.05; Fig. 4A; Supplementary Table S2). Goldmine was used to reveal the genomic context of differentially methylated regions (DMR) of guadecitabine-treated adipocytes. The distribution of DMRs after guadecitabine treatment showed that hypomethylation in promoter (>0.03), exon (>0.02), intron (>0.15), and 3′ end (>0.02) regions and hypermethylation in intergenic region (>0.20) was increased relative to control adipocytes (Fig. 4B). IPA of the demethylated genes generated 52 pathways enriched (FDR < 0.05) in cell signaling, metabolic, and immunogenic pathways (Fig. 4C–E; Supplementary Table S3 for full list of pathways).
To investigate whether the effects of demethylation in adipocytes led to reexpression of genes after guadecitabine treatment, poly(A) tail RNA-seq was performed on the same adipocytes used for MC-seq. Although over 3,000 genes were demethylated, only 753 genes were upregulated after guadecitabine treatment (fold change > 1.2, P < 0.05) and 712 genes were downregulated, as seen in the volcano plot (Fig. 5A) and heatmap (Fig. 5B; Supplementary Table S4 for full list of differentially expressed genes). IPA generated 57 pathways enriched (FDR < 0.05) in cell signaling, metabolic, and immunogenic pathways (Fig. 5C–E; Supplementary Table S3 for full list of pathways generated by 286 differentially expressed genes), similar to the pathways enriched by the differentially methylated genes. Six of 7 genes in IL8 signaling pathway (PGF, RAC2, GNG2, BCL2L1, EGFR, RHOJ, and VEGFA) were validated using qRT-PCR (Supplementary Fig. S6B) compared with RNA-seq data (Supplementary Fig. S6A).
Integrated analysis of MC- and RNA-seq data generated 85 genes that remained hypermethylated and downregulated and 279 genes that were hypomethylated and upregulated (Supplementary Table S5) by guadecitabine. IPA analysis of the hypomethylated and upregulated genes generated 74 pathways again enriched (FDR < 0.05) for cell signaling, metabolic, and immunogenic pathways (Supplementary Fig. S5B–S5D; Supplementary Table S3 for full list of integrated analysis pathways).
To identify potential changes in methylation and gene expression due to the two-way communication between adipocytes and OC cells, we indirectly cocultured OVCAR5-RFP cells using a 0.45 μm Boyden chamber with adipocytes for 4 days, collected the RNA from adipocytes, and performed RNA-seq. The top 5 altered pathways in adipocytes revealed that coculture with OC cells suppressed LXR/RXR activation pathway and activated IL6 signaling, leukocyte extravasation signaling, dendritic cell maturation, and IL8 signaling compared with pathway alterations in adipocytes treated with guadecitabine (FC > 1.2; FDR < 0.05; Supplementary Fig. S6D; full list in Supplementary Table S6). These results demonstrate that communication between OC cells and adipocytes may lead to changes in pathways in adipocytes.
SUSD2 upregulation and secretion decrease migration and invasion
It was of interest to examine the above results for genes associated with cytokines or other secretory factors known to inhibit OC cell progression. Demethylated and reexpressed genes of interest (FDR < 0.05) included SUSD2, TFP12, GREM1, NBL1, TMBIM1, LY6K, ISG15, PLA2G4A, TRIM29, PSG4, PRDX2, ACOT7, GPAT2, PLA2G5, and SLC25A1, and these were validated by qRT-PCR in the same RNA samples used for sequencing (Fig. 6A). Separate adipocyte samples treated with guadecitabine were tested, and only SUSD2, TFPI2, GREM1, PLA2G4A, and TRIM29 corresponded with increased (P < 0.05) expression (Fig. 6B). Of these genes, significant epigenetic reexpression (fold change = 3.15; Supplementary Table S4) of SUSD2 (sushi domain containing 2) and demethylation of SUSD2 in the promoter region with reexpression in the gene body was observed (Fig. 6C; Integrated Genome Viewer).
Previously described as a secreted tumor suppressor in HGSOC (26), it was of interest to test whether guadecitabine treatment altered secretion of SUSD2 in the adipocyte culture media. SUSD2 ELISA was performed on ACM (control vs. guadecitabine), and results demonstrated that guadecitabine treatment increased (P < 0.01) SUSD2 secretion (1.48-fold) from 46.3 ± 3.37ng/mL to 67.6 ± 3.91 ng/mL (Fig. 7A). Also, flow cytometry analysis of surface level SUSD2 on adipocytes demonstrated an increase (P < 0.05) after guadecitabine treatment (Fig. 7B). Furthermore, treatment with recombinant SUSD2 (100 ng/mL, 500 ng/mL placed in the bottom of the Boyden chamber) decreased (P < 0.05) OVCAR8 (fold change > 81%) and OVCAR4 (fold change > 30%) cell migration and OVCAR8 (fold change > 42%) and OVCAR4 (fold change > 41%) cell invasion (Fig. 7C). However, recombinant SUSD2 had no effect on OVCAR4 (Fig. 7D) and OVCAR8 proliferation (Fig. 7E) compared with carboplatin alone. The use of 100 ng/mL SUSD2 used in the migration, invasion, and proliferation assays is higher than levels determined by ELISA (46–68 ng/mL) but not unreasonable. Furthermore, higher dose SUSD2 (500 ng/mL) had no additional effect in migration, invasion, and proliferation, suggesting a level of saturation.
To note, SUSD2 mRNA levels of the cell lines were calculated using the slope generated by linear regression analysis of a SUSD2 primer efficiency curve generated by qRT-PCR (Supplementary Fig. S6C). SUSD2 concentrations for OVCAR4, OVCAR5, OVCAR8, Kuramochi, and SKOV3 cells were 39.77, 103.51, 10.67, 3.11, and 0.95 ng/μL, respectively. Collectively, these results demonstrated that DNMT inhibition of adipocytes increased SUSD2 secretion resulting in decreased HGSOC cell migration and invasion toward adipocytes.
Discussion
The majority of OC patients present at late stage with extensive metastases to the omentum. Adipocytes have been associated with metastasis of ovarian and other cancers (1). Adipose tissue represents a hypoxic area with elevated levels of proinflammatory factors and increased macrophage invasion (27). Mice fed a “Western” diet revealed increased adhesion to mesothelial explants and ovarian tumor burden compared with control (5), demonstrating that the adipocyte-rich metastatic niche plays a key role in OC progression. Despite the prominent role adipocytes play in OC, the molecular mechanism involved in increased metastasis remains incompletely understood. In this study, we demonstrate adipocytes increase migration and invasion of OC cells, and that epigenetic drug treatment of adipocytes reverses this trend by upregulating SUSD2 and other tumor-suppressor genes.
The omentum is a complex organ with endocrine, metabolic, and immune-regulatory roles. Recent findings have shown cancer-associated adipocytes secrete adipokines that stimulate adhesion, migration, and invasion of cancer cells (1). Here, we demonstrate that adipocytes increase migration and invasion of OC cells. Once associated with adipocytes, OC cells are able to engage in exchange of metabolites that promote tumorigenesis, including direct transfer of lipids through FABP4 (3) and adipocyte activated SIK2 regulation of fatty acid oxidation in OC cells (6). Adipocytes have been demonstrated to protect cancer cells from therapies by enhancing EMT process, upregulating prosurvival genes, increasing oxidative stress response, activating autophagy, increasing adipocyte-secreted glutamine, and dysregulating of genes involved in gemcitabine transport and metabolism (28). Our findings demonstrate that direct coculture of carboplatin-treated OVCAR8-RFP cells and adipocytes increases OVCAR8-RFP cell proliferation. However, there was no change in cell proliferation after guadecitabine treatment compared with carboplatin alone, suggesting that treatment of adipocytes with guadecitabine primarily affects migration and invasion. These results further indicate the importance of adipocytes in OC progression and suggest targeting adipocytes as a potential therapeutic strategy in OC.
DNA methylation changes after DNMT1 inhibition in adipose-derived stem cells (29), differentiating human adipocytes (30), and 3T3-L1 preadipocyte differentiation (31) have been reported, and epigenetic alterations in cancer cells in the context of metabolic remodeling caused by obesity have been examined (32). Although altered DNA methylation of important metabolic regulators in adipocytes of obese subjects has been described, including hypermethylation of adiponectin (7) and hypomethylation of leptin (33), epigenetic alterations of adipocytes associated with cancer have not been examined. In this study, we demonstrate that by altering DNA methylation and upregulating adipokines involved in negative regulation of OC cell homing, guadecitabine treatment of adipocytes decreases migration and invasion of OC cells. In agreement with observations of others demonstrating that adipocytes induce EMT of OC cells (3), our data show that guadecitabine treatment of adipocytes alters adipokine secretion and subsequently altered EMT markers in OC cells, including CDH1, SLUG, FN1, and CDH2, further suggesting that epigenetic therapy of adipocytes can influence proximal cancer cells.
The known molecular mechanism of cytidine analogue DNMT inhibitors is replication-dependent (34); however, mature adipocytes do not divide. We took several approaches (BrdUrd immunohistochemical staining, MTT proliferation assay, and propidium iodide cell-cycle analysis) to confirm that the adipocytes used in this study were neither dividing nor synthesizing new DNA. In support of our results, 5-azacitidine induces TNFα demethylation in mature 3T3-L1 cells (35), and treatment of 3T3-L1 cells with 5-aza-dC at the late stage of differentiation markedly enhances expression of adipocyte lipogenic markers, such as Fas and Srebp1c (31). Furthermore, we show that guadecitabine decreases DNMT1 levels in both soluble and chromatin subcellular fractions of mature adipocytes despite the presence of DNA synthesis inhibitor, and previous studies report DNMT1 degradation following 5-aza-deoxycytidine in nondividing cancer cells (23, 36). Interestingly, covalent bond formation between DNMT1 and 5-aza-deoxycytidine–incorporated DNA is not essential for DNMT1 enzyme degradation (36), and 5-aza-deoxycytidine–mediated DNMT1 degradation is not absolutely associated with cell-cycle progression in HCT116 colon cancer cells but is associated with their p53 genotype (23). These studies suggest an alternative mechanism of cytosine analogue-induced DNMT1 degradation that is not strictly dependent on cell division. Our data suggest that DNMT inhibition can occur in nondividing adipocytes and elicit antimetastatic effects on OC cells.
The majority of studies regarding adipose tissue genome-wide DNA methylation were performed in the context of body mass index, exercise, or adipogenesis (37–39). However, genome-wide DNA methylation alterations after DNMT inhibition of mature adipocytes remain relatively unexplored. We show that guadecitabine treatment of adipocytes in vitro may alter cell signaling pathways including ERK/MAPK and cAMP-mediated signaling, fatty acid activation metabolic signaling, and toll-like receptor signaling, all of which have been implicated in regulating metastasis and cancer metabolism (40, 41). Our future studies will include validation of IL8 pathway–related genes beyond mRNA expression. In addition, we observed upregulation of SUSD2 (adhesion molecule), TFPI2 (serine protease inhibitor), and GREM1 (BMP antagonist) in adipocytes treated with guadecitabine, and pathway analysis demonstrates inhibition of matrix metalloprotease and EMT pathways and suggests that hypomethylating agents have an indirect effect on OC cells through altered adipokine secretion.
DNMT inhibitor treatment has previously been demonstrated to induce interferon response by derepressing hypermethylated endogenous retrovirus genes and subsequently sensitizing OC cells to immune checkpoint therapy (42). Our observation is that guadecitabine upregulates PLA2G4A and TRIM29, both of which are involved in inflammatory response (43) and IFNγ signaling (44), respectively. RNA-seq of adipocytes cocultured with OC cells demonstrates that two-way communication between OC cells and adipocytes also affects adipocyte gene expression and pathways. IL6 and IL8 pathways in adipocytes are upregulated after coculture with OC cells, demonstrating an increase in inflammatory cytokines that are known to increase EMT. Also, LXR/RXR activation pathway, which plays a major role in lipid homeostasis, is downregulated (FC > 1.2, FDR < 0.05), suggesting an alteration in lipid transport in adipocytes when cultured in proximity with OC cells. Based on these collective results, we suggest that epigenetic therapy can indirectly inhibit OC metastasis by modifying signaling pathways in adipocytes and possibly alter the immune response by inducing inflammatory response genes and lipid transport in adipocytes in the tumor microenvironment.
SUSD2 has been shown to be a tumor suppressor in non–small cell lung (45), colon cancer (46), renal cell (47), hepatocellular (48), and ovarian cancer (26), and high SUSD2 levels associate with improved survival in OC (26). SUSD2 reexpression in renal cell carcinoma and lung cancer inhibited clonogenicity and proliferation (47), and SUSD2 knockdown in OVCAR3, OVSAHO, and Kuramochi cell lines increases migration and mesothelial clearance ability (26). Our data show guadecitabine-treated adipocytes upregulate SUSD2 secretion into the media and recombinant SUSD2 decreases HGSOC cell migration and invasion. SUSD2 is localized on the plasma membrane where the N-terminus is cleaved with only the transmembrane C-terminus domain remaining (49). Our study demonstrates that soluble SUSD2 interacts with OC cells to decrease metastatic behavior. Only two reports detail SUSD2-binding partners, galectin-1 (Gal-1; binds β-galactoside residues) in breast cancer (50) and colon-derived SUSD2-binding factor (CSBF) in colon cancer (46). From seven different publicly available protein–protein interaction datasets, several SUSD2-binding partners were predicted (Supplementary Table S1). DLG4 is the only membrane-associated protein that could be a SUSD2 receptor. However, DLG4–SUSD2 interaction has yet to be demonstrated. Further studies are needed to define the exact binding partner of SUSD2 in the context of OC and adipocytes. In addition, we show that recombinant SUSD2 interacts with OC cells to decrease migration and invasion, warranting further investigation into the intracellular mechanism.
Our study is the first to examine the effects of a DNMT inhibitor on adipocytes in the context of the tumor microenvironment. Guadecitabine treatment of adipocytes decreases HGSOC migration and invasion at least partly through SUSD2 upregulation. These findings suggest that in addition to direct effects on cancer cells, epigenetic therapy may affect the tumor microenvironment through altered cytokine secretion to prevent OC metastasis.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Authors' Contributions
Conception and design: J. Tang, A. Özeş, M.S. Stack, K.P. Nephew
Development of methodology: J. Tang, A. Özeş, H. O'Hagan
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): J. Tang, N. Pulliam, K.P. Nephew
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): J. Tang, A. Buechlein, D. Rusch, K.P. Nephew
Writing, review, and/or revision of the manuscript: J. Tang, H. Keer, H. O'Hagan, M.S. Stack, K.P. Nephew
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): J. Tang, N. Ding, H. Keer, K.P. Nephew
Study supervision: K.P. Nephew
Acknowledgments
This work was funded by the NCI Award CA182832-01, the Ovarian Cancer Research Fund Alliance (PPDIU01.2011), National Institute of Environmental Health Sciences (R01ES023183; to H. O'Hagan), the V-Foundation, Indiana Clinical and Translational Sciences Institute, and the Medical Sciences Doane and Eunice Wright Memorial Fellowship.
The authors thank Dr. F. Fang and Vaishnavi Murali for technical assistance and Drs. Mohammad Azab and Pietro Taverna (Astex Pharmaceuticals, Inc.) for providing guadecitabine. We thank Dr. Martin Wabitsch (University Medical Center Ulm, Germany) for SGBS cells. We thank Dr. Anirban Mitra for helpful discussions.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.