1. Introduction
Pancreatic ductal adenocarcinoma (PDAC) currently accounts for approximately 3% of all new cancer diagnoses but for 8% of all cancer deaths, making it the third leading cause of cancer-related death in the United States [
1]. The aggressive nature of PDAC is multifactorial and can be attributed to advanced stage at diagnosis and resistance to radiation and chemotherapy. Overall PDAC patients have a 5-year relative survival of about 11%; patients with metastatic disease have only a 3% 5-year relative survival, an especially grim prognosis considering nearly 50% of patients are found to have metastatic disease at their initial diagnosis [
1]. Current chemotherapy regimens such as FOLFIRINOX (5-fluorouracil, oxaliplatin, irinotecan, and leucovorin) and gemcitabine plus nab-paclitaxel have offered modest improvements in survival [
2]. These improvements have yielded median survival times of just 8-11 months in patients with metastatic disease and less than 2 years in patients with locally advanced disease [
3,
4]. Thus, new treatment modalities are critically important in improving survival among PDAC patients.
Pharmacologic ascorbate (P-AscH
-, high-dose intravenous vitamin C) was first proposed as a cancer treatment in the 1970s but interest soon waned. However, it has experienced a resurgence in interest due to multiple clinical trials with encouraging results [
5]. P-AscH
- acts as a pro-oxidant, delivering high concentrations of extracellular hydrogen peroxide (H
2O
2) to the cancer cell microenvironment [
6,
7,
8]. H
2O
2 can easily cross the cell membrane and form the highly reactive hydroxyl radical (HO
∙), resulting in protein and DNA damage and ultimately cell death [
6,
7,
8,
9]. This H
2O
2-induced cytotoxicity has been shown to be selective to cancer cells due to their overall higher endogenous levels of reactive oxygen species (ROS) as well as lower levels of peroxide-removing enzymes such as catalase and glutathione peroxidase when compared to normal cells [
7,
8,
10,
11,
12]. P-AscH
- has shown promise as an effective adjunct to standard of care chemotherapy in phase I and II clinical trials in PDAC, lung cancer, and glioblastoma in recent years [
9,
13,
14,
15].
In addition to its direct cytotoxic effects at pharmacologic doses, ascorbate has also been shown to act as an important cofactor for the activity of the ten-eleven translocation (TET) methylcytosine dioxygenase family of enzymes. These enzymes are responsible for removing methyl groups on cytosine bases in DNA, converting 5-methylcytosine (5-mC) to 5-hydroxymethylcytosine (5-hmC) [
16,
17]. Conversely, DNA methyltransferases (DNMT1, DNMT3A, DNMT3B) add methyl groups to cytosine bases [
18,
19]. The methylation status of 5’ cytosine bases within promoter regions and CpG islands influences gene expression and can therefore have a significant impact on the phenotype of tissues and organisms. Hypermethylation of promoter regions decreases downstream gene expression likely by altering chromatin organization and function as well as interaction with transcriptional repressors (methyl-CpG-binding proteins, MeCP1 and MeCP2) [
19,
20,
21]. Hypermethylation within promoter regions of tumor suppressor and DNA repair genes has been identified in every major tumor type, including PDAC [
22,
23,
24]. Relatedly, overexpression of DNMT has been associated with more aggressive forms of cancer, suggesting that DNMT overexpression leading to a hypermethylated state results in aberrant gene expression and decreased expression of tumor suppressor genes [
25,
26,
27]. Studies have shown that PDAC is hypermethylated and exhibits DNMT1 overexpression, leading to reduced expression of many tumor suppressor genes including p16 and APC [
23,
28]. Previous studies have also demonstrated that DNMT1 inhibition in PDAC reduces cancer cell growth both
in vitro and
in vivo [
29,
30]. Inhibitors of DNMT such as 5-azacytidine and 5-aza-2’-deoxycytidine are frequently used in hematologic malignancies, such as myelodysplastic syndrome and AML, but their use is being increasingly investigated in solid organ tumors with promising results in early studies in colorectal cancer, ovarian cancer, and PDAC [
28,
31,
32,
33]. Reducing the degree of DNA methylation in PDAC may result in increased expression of epigenetically silenced tumor suppressor genes and help to enhance the cytotoxicity of other chemotherapeutic agents.
Two genes that are silenced in PDAC are dual oxidase 1 and 2 (DUOX 1 and 2), members of the NADPH oxidase family of enzymes that produce extracellular hydrogen peroxide (H
2O
2) in normal tissues throughout the body [
34,
35]. These genes have been shown to be epigenetically silenced in lung cancer, hepatocellular carcinoma, and PDAC, and their expression in PDAC is increased for up to 72 h following treatment with P-AscH
-, leading to sustained oxidative stress and cytotoxicity [
34,
36,
37]. We hypothesized that P-AscH
- and DNMT inhibitors may act synergistically to increase DUOX expression, leading to increased H
2O
2-induced cytotoxicity in PDAC. Here we demonstrate that P-AscH
- and DNMT inhibitors act synergistically to reduce PDAC cell growth in a dose-dependent manner through an H
2O
2-dependent pathway secondary to increased DUOX expression. We also demonstrate the efficacy of P-AscH
- and DNMT inhibitors
in vivo, where tumor volume was significantly reduced and DUOX1 expression was increased in a xenograft model. Combination treatment with P-AscH
- and DNMT inhibitors may offer an epigenetic approach to the treatment of PDAC.
2. Materials and Methods
2.1. Survival Curve Analysis
The NIH Genomic Data Commons (GDC) Cancer Genome Atlas (TCGA) Pancreatic Cancer Database (PAAD) (
n = 223) was accessed
via the University of California Santa Cruz (UCSC) Xena Functional Genomics Explorer [
38]. Compared variables were “DNA methylation – Illumina Human Methylation 450” and “OS” (overall survival). The dataset was filtered using “sample_type.samples,” “primary tumor,” “morphology.diagnoses,” and “8500/3” to only include primary tumor samples from patients with invasive ductal adenocarcinoma. Adenocarcinoma not otherwise specified (
n = 24), neuroendocrine carcinoma (
n = 9), mucinous adenocarcinoma (
n = 6), adenocarcinoma with mixed subtypes (
n = 1), and undifferentiated carcinoma (
n = 1) were excluded. After excluding normal solid tissue (
n = 37) and metastatic (
n = 1) samples, the final number of samples for evaluation of DUOX1 methylation status was
n = 148. DUOX1 methylation status was stratified into two groups by dividing the filtered data set at the median into a low methylation and a high methylation group. Both groups were plotted on a Kaplan Meier curve and compared using the Gehan-Breslow-Wilcoxon test.
2.2. Cell Culture
Human PDAC cell lines MIA PaCa-2 and PANC-1 were cultured in DMEM (Gibco, 11965) supplemented with 10% FBS (Gibco, 26140) and 1% penicillin–streptomycin antibiotic (Gibco, 15140). The human patient-derived PDAC cell line PDX-339 was cultured in DMEM/F-12 media (Gibco, 11320) supplemented with 10% FBS, insulin (Gibco, 12585), EGF (Gibco, PHG0311), hydrocortisone (Sigma, H0888), bovine pituitary extract (Gibco, 13028), and 1% penicillin–streptomycin antibiotic. Human cell lines (MIA PaCa-2 and PANC-1) were purchased directly from the ATCC and were passaged for fewer than 6 months after receipt. No additional authentication was performed. The patient-derived cell line (PDX-339) was obtained from the Medical College of Wisconsin surgical oncology tissue bank. Regardless of varying cell type and media components, all cells were treated with ascorbate in fresh 10% DMEM media for 1 h at 37 °C. Media was then replaced, ascorbate was removed, and cells were allowed to incubate for 48 h. Ascorbate came from a stock solution of 1 mol/L (pH 7) made under argon and stored with a tight-fitting stopper at 4 °C. Ascorbate concentration was verified at 265 nm, ε
265 = 14,500 M
−1 cm
−1 [
39]. To enhance rigor and reproducibility, final concentrations were calculated in units of moles-per-cell to account for variation in media, cell density, and cellular metabolism [
40]. A 1 M 5-azacytidine (AZC) (Sigma, A2385) solution was made in distilled H
2O, passed through a 0.22 µm filter, aliquoted and stored at
-80 °C. For
in vitro use, a 100 µM 5-aza-2’-deoxycytidine (AZD) (Sigma, A3656) solution was made in distilled H
2O, passed through a 0.22 µm filter, aliquoted and stored at -80 °C. For
in vivo use, a separate 438 µM AZD (Tocris, 2624) solution was made in PBS -/- (Gibco, 14190144), passed through a 0.22 µm filter, aliquoted and stored at -80 °C.
qRT-PCR-1 × 105 cells were seeded for 5-7 d and then treated with 0.5-2 µM AZC for 5 d or 0.1-1 µM AZD for 3 d with fresh media and AZC or AZD replaced daily, with and without 10-20 pmol/cell P-AscH- for 1 h. Following 5 d of AZC treatment or 3 d of AZD treatment and 48 h after P-AscH-, total cellular RNA was isolated using Trizol reagent (Invitrogen, 15596026) in conjunction with the Direct-zol RNA Miniprep Kit (Zymo Research) according to the manufacturer's protocol. RNA was quantified using a ND1000 Nano-Drop spectrophotometer (Nano-drop). cDNA was synthesized using the high-capacity cDNA archive kit (Applied Biosystems). qRT-PCR assays were performed using 2x Power SYBR Green real-time master mix (Applied Biosystems, 4368702) under the following set up: 95 °C for 10 min, followed by 40 cycles of 95 °C for 15 s and 60 °C for 1 min (StepOne plus Sequence Detection System, Applied Biosystems). Primers were ordered from Integrated DNA Technology (DUOX1 gene ID: 53905, DUOX2 gene ID: 50506).
Primer sequences:
DUOX1 Forward-GGATGCTGAGATCCTTCATCGAGA.
DUOX1 Reverse-ACCTCCACCCCTTTGACACAGAG.
DUOX2 Forward-TGTGTATGAGTGGCTGCCCAGC.
DUOX2 Reverse-ACTGCTCAGAGGCCACCACAAA.
Western Blotting-Protein samples were isolated and prepared in phosphosafe buffer (EDM Millipore, 71296) containing protease inhibitor cocktail (Sigma). Protein samples cells were isolated in 1x RIPA buffer with PhosSTOP phosphatase inhibitor cocktail (Sigma, 04906845001) and cOmplete Mini protease inhibitor cocktail (Sigma, 11836153001). Protein concentration was measured using a Bradford protein assay. Total protein (30-40 µg) was loaded on a 4-20% SDS-PAGE gradient gel (Bio-Rad). Membranes were blocked in 5% BSA in TBS-T. Primary antibodies included: DUOX1 (Santa Cruz Biotechnology, H-9, 1:1000), DUOX2 (Santa Cruz Biotechnology, E-8, 1:1000), DNMT1 (Santa Cruz Biotechnology, H-12, 1:1000), Tubulin (1:500-1000, University of Iowa Developmental Studies Hybridoma Bank, RRID: AB_528499), and GAPDH (Cell Signaling, D16H11, 1:1000). Appropriate horseradish peroxidase-linked secondary antibodies were used at a concentration of 1:20,000 to 1:50,000. Blots were visualized using SuperSignal West Pico PLUS substrate (Thermo Scientific, 34580) on X-ray film.
2.3. Measurement of Pro-Oxidants
Bovine catalase (Sigma, C40) was utilized in the media at 100 mg/mL for 5 d along with 2 µM AZC replaced daily and immediately prior to P-AscH- treatment (PDX-339, 1 mmol/L which is equivalent to 10 pmol cell−1) for 1 h and fresh media was replaced. Cells were then allowed to proliferate for 48 h prior to the DCFH-DA [5-(and-6)-carboxy-20 ,70 -dichlorodihydrofluorescein diacetate; Molecular Probes, C400] assay. Cells were incubated at 37 °C and protected from light in 15 mmol/L DCFH-DA for 30 minutes in PBS. They were then harvested and filtered for flow cytometry analysis. A CBD LSRII cytometer (BD Biosciences) was used to measure DCFH-DA oxidation at 504/529 nm. Data were analyzed using FlowJo Software.
2.4. Clonogenic Survival
1 x 105 cells were seeded 5-7 d prior to assay. Cells were treated alone or in combination with 0.5-2 µM AZC for 5 d, 0.1-0.3 µM AZD for 3 d, 10-20 pmol/cell P-AscH- for 1 h, and bovine catalase (Sigma, C40) at 100 mg/mL. Following 5 d of AZC treatment or 3 d of AZD treatment and immediately after P-AscH- treatment, cells were trypsinized with TrypLE Express (Gibco, 12604) to form a single-cell suspension and counted using a Countess II automated cell counter (Thermo Fisher) to determine the number of cells plated into each well. Cells were cultured for 10-14 days before being fixed and stained for analysis of surviving fraction. Colonies containing ≥ 50 cells were scored. Surviving fraction was defined as number of colonies counted/number of cells seeded. Each experimental condition was normalized to their own control to determine a normalized surviving fraction. Each condition was done in triplicate, and experiments were performed at no less than n = 3.
2.5. In Vivo Experiments
The animal protocols were reviewed and approved by the Animal Care and Use Committee of The University of Iowa. Thirty-day old athymic nude mice (Foxn1nu) were obtained from Envigo and animals were allowed to acclimate in the unit for 1 week prior to any manipulation. MIA PaCa-2 cells (2 x 106) were injected subcutaneously into each flank region of nude mice with a 1 mL tuberculin syringe equipped with a 25-gauge needle. Tumors grew to approximately 5 mm in diameter before experimental treatment began. Mice were divided into four treatment groups with equal tumor volume. Mice were treated once daily with I.P. saline (1 M), P-AscH- (4 g/kg), AZD (1 g/kg), or a combination of P-AscH- and AZD for 21 days. Tumor volume was measured twice weekly using handheld calipers. Mice were euthanized and tumors were harvested and processed for experimental analyses.
Immunofluorescent Staining-Mouse tumor tissue was fixed in paraformaldehyde and imbedded in paraffin. Tissues were sectioned at 6–10 µm and placed on SuperFrost Plus slides and baked overnight at 42-45 °C. Slides were de-paraffinized and were further processed for immunohistochemistry. Immunofluorescent samples were subject to antigen unmasking protocol at 95 °C for 15 min, cooled to room temperature, washed with double distilled water, and allowed to dry. Samples were blocked with 5% normal goat serum prior to incubation with DUOX1 (1:50) and DNMT1 (1:10) primary antibodies in 1% normal goat serum at 4 °C for 24 h. Slides were washed 3 times in PBS for 5 min. Goat anti-mouse secondary (1:400) conjugated to FITC in 1% normal goat serum at 4 °C overnight. Slides were washed 3 times with PBS for 5 min. Nuclear staining was performed using Topoisomerase-3 (1:1000) for 15 min and samples were washed once with PBS. Coverslips were mounted on slides with immunofluorescence mounting medium and visualized using a Zeiss Confocal Microscope 40x oil objective. Quantification was performed using only DUOX1 or DNMT1 staining images. ImageJ was used to quantify the fluorescence for each image and values were normalized to each image’s nuclear content.
2.6. Statistical Methods
Data are presented as the mean ± SEM. For statistical analyses of two groups, unpaired 2-tailed Student’s t-tests were utilized. To study statistical differences between multiple comparisons, significance was determined using one-way ANOVA analysis with Tukey’s multiple-comparisons test. For survival analysis, Gehan-Breslow-Wilcoxon test was used to compare overall survival. All analyses were performed in GraphPad Prism 9.0 (GraphPad Software, Inc.).
2.7. Data Availability
All data are contained within the manuscript.
4. Discussion
Even with advancements in chemotherapy regimens and radiation, median survival for PDAC is still extremely low, with a 5-year relative survival of just 11% across all stages and median survival of approximately 18 months for locally advanced disease [
1,
2,
3,
4]. P-AscH
- has shown promise as an adjunct to current chemo- and radiotherapy in recent phase I and II clinical trials [
13,
14,
15]. Identifying therapies that can act synergistically with P-AscH
- would offer a novel approach to improve cancer-specific cytotoxicity. Our results demonstrate that P-AscH
- and DNMT inhibitors have at least additive and potentially synergistic effects on increasing cytotoxicity in pancreatic cancer cells.
Epigenetic modifications of tumor suppressor genes have been studied extensively. In PDAC specifically, the NADPH oxidase enzymes DUOX1 and DUOX2 are downregulated, but their expression is increased following exposure to H
2O
2 or P-AscH
- [
34,
43]. Available data suggest that DUOX1 is also often suppressed in hepatocellular carcinoma, a carcinoma of epithelial origin similar to PDAC [
37,
46]. In hepatocellular carcinoma, DUOX1 expression has been found to correlate with improved overall survival [
46]. As seen in
Figure 1A, we demonstrate the association of DUOX1 hypermethylation with decreased overall survival in PDAC. Thus, therapies that can act with P-AscH
- to increase DUOX expression and thereby increase downstream H
2O
2 production may benefit PDAC treatment.
Recent studies demonstrate that DNMT1 is overexpressed in PDAC, and its expression increases with progressing malignant potential,
i.e. from pre-neoplastic lesions to invasive cancer [
47]. This suggests that DNMT inhibitors may be effective adjuncts for PDAC treatment and offer an effective pathway to increase the expression of epigenetically silenced tumor suppressor genes, including DUOX. We demonstrate that commonly utilized DNMT inhibitors (5-azacytidine [AZC] and 5-aza-2’-deoxycytidine [AZD]) produce sustained increases in DUOX expression, leading to increased H
2O
2 production and H
2O
2-induced cytotoxicity. Our present study also demonstrates that modest concentrations of AZC (1 µM) and AZD (0.5-1 µM) produce significant increases in DUOX1 and DUOX2 expression across multiple PDAC cell lines compared to control [
41,
44]. These increases in DUOX expression are dose-dependent and significantly increased with the addition of P-AscH
-, with DUOX1 mRNA expression increased more than 80-fold following treatment with AZD and P-AscH
-. Furthermore, the increase in DUOX1 expression is sustained for up to 72 h in the absence of treatment, which was also seen
in vivo, where DUOX1 expression in the combination treatment group is increased compared to either treatment alone. The pathway for this large increase in enzyme expression is likely multifaceted and includes decreased genomic methylation secondary to decreased DNMT activity and increased degradation by DNMT inhibitors, decreased genomic methylation secondary to increased TET activity by P-AscH
-, and the H
2O
2-induced increased expression and activity of DUOX enzymes from P-AscH
- and the endogenously produced H
2O
2 [
7,
8,
16,
17]. This robust pathway to increased DUOX expression may offer a consistent, sustained mechanism for increasing H
2O
2 production in PDAC.
Increased DUOX expression correlates with increased production of H
2O
2. The PDX-339 cell line was chosen for this experiment as PDX cell lines have been shown to retain native tumor heterogeneity
in vitro and may more accurately reproduce clinical response [
48]. Multiple previous studies have demonstrated the cancer-specific cytotoxic properties of H
2O
2 in PDAC [
7,
8,
10,
11]. Across all PDAC cell lines, we demonstrate significantly increased cytotoxicity and decreased clonogenic survival when DNMT inhibitors are combined with P-AscH
- compared with P-AscH
- alone, with surviving fractions less than 10% at conservative treatment doses. Similar to DUOX expression, these effects are dose-dependent. With the addition of catalase to the media prior to AZD or P-AscH
- treatment, the decreased clonogenic survival of PDAC cells was reversed, indicating a H
2O
2-mediated mechanism of cytotoxicity following DNMT inhibitor treatment. These results are replicated
in vivo, where the combination of AZD and P-AscH
- demonstrates inhibited tumor growth in a pancreatic cancer xenograft model. While P-AscH
- and DNMT inhibitors have been utilized independently in
in vivo models of PDAC, this is the first study to demonstrate the combination leading to tumor growth inhibition [
14,
30,
33,
47].
Hypoxia is a well-known modulator of the PDAC microenvironment. Low oxygen content can impact enzyme function, cellular stability, DNA repair mechanisms, redox states, sensitivity to therapies, and gene expression [
49,
50]. We demonstrate that hypoxia increases DNMT1 protein expression in PDAC and that P-AscH
- reverses this effect. Similar results were seen in the
in vivo tumor xenografts suggesting that hypoxia could be another factor leading to DNMT1 overexpression and hypermethylation in PDAC.
Treatments for PDAC must continue to evolve. Epigenetic alterations provide new opportunities for PDAC treatment. We have shown that sustained increases in DUOX expression can be achieved with DNMT inhibitors and that P-AscH- enhances the cytotoxic and epigenetic effects of DNMT inhibitors through a peroxide-mediated mechanism. P-AscH- decreases DNMT1 expression in hypoxia in vitro and in vivo, providing an additional mechanism to limit DNMT1 activity and DNA methylation. Finally, the combination of a DNMT inhibitor and P-AscH- in vivo increases DUOX1 expression and decreases tumor volume compared to either treatment alone. Phase I and II clinical trials of DNMT inhibitors in PDAC are currently underway, where these therapies in combination with P-AscH- may provide new avenues for PDAC treatment through an epigenetic mechanism.