<figure><img src="image_1.png" /></figure>

RESEARCH ARTICLE

<figure><img src="image_2.png" /></figure>

WILEY

# E2F transcription factors promote tumorigenicity in pancreatic ductal adenocarcinoma

Ludivine Bertonnier-Brouty<sup>1,2</sup> | Jonas Andersson<sup>2</sup> | Tuomas Kaprio<sup>3,4,5</sup> |

<figure><img src="image_3.png" /></figure>

Jaana Hagström<sup>3,4,5,6</sup> | Sara Bsharat<sup>1,2</sup> | Olof Asplund<sup>2</sup> | Gad Hatem<sup>2</sup> |

Caj Haglund<sup>3,4,5</sup> | Hanna Seppänen<sup>3,4,5</sup> | Rashmi B. Prasad<sup>2</sup> | Isabella Artner<sup>1,2</sup>

## Correspondence

Isabella Artner, Lund University Diabetes Center, Lund University, Jan Waldenströms gata 35, Malmö 21428, Sweden.

Email: isabella.artner@med.lu.se

## Funding information

Cancer and Allergy Foundation (2022, 2023); Novonordiskfonden, Grant/Award Number: NNF20OC0063485; The Swedish Cancer foundation (2022); Sigrid Juselius Foundation; Vetenskapsrådet, Grant/Award Number: 2020-0146; Fredrik and Ingrid Thuring foundation, Grant/Award Number: 2020-00596; The Royal Physiographic Society (2022)

## Abstract

**Background:** Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal cancers with limited treatment options, illustrating an urgent need to identify new drugable targets in PDACs.

**Objective:** Using the similarities between tumor development and normal embryonic development, which is accompanied by rapid cell expansion, we aimed to identify and characterize embryonic signaling pathways that were reinitiated during tumor formation and expansion.

**Methods and Results:** Here, we report that the transcription factors E2F1 and E2F8 are potential key regulators in PDAC. E2F1 and E2F8 RNA expression is mainly localized in proliferating cells in the developing pancreas and in malignant ductal cells in PDAC. Silencing of E2F1 and E2F8 in PANC-1 pancreatic tumor cells inhibited cell proliferation and impaired cell spreading and migration. Moreover, loss of E2F1 also affected cell viability and apoptosis with E2F expression in PDAC tissues correlating with expression of apoptosis and mitosis pathway genes, suggesting that E2F factors promote cell cycle regulation and tumorigenesis in PDAC cells.

**Conclusion:** Our findings illustrate that E2F1 and E2F8 transcription factors are expressed in pancreatic progenitor and PDAC cells, where they contribute to tumor cell expansion by regulation of cell proliferation, viability, and cell migration making these genes attractive therapeutic targets and potential prognostic markers for pancreatic cancer.

## KEYWORDS

E2F1, E2F8, fetal pancreas, pancreatic cancer, pancreatic ductal adenocarcinoma

<footnote n="1"><sup>1</sup>Lund Stem Cell Center, Lund University, Lund, Sweden</footnote>

<footnote n="2"><sup>2</sup>Lund University Diabetes Center, Lund University, Malmö, Sweden</footnote>

<footnote n="3"><sup>3</sup>Department of Surgery, Helsinki University Hospital, Helsinki, Finland</footnote>

<footnote n="4"><sup>4</sup>Translational Cancer Medicine Research Program, Faculty of Medicine, University of Helsinki, Helsinki, Finland</footnote>

<footnote n="5"><sup>5</sup>iCAN, Digital Cancer Precision Medicine, University of Helsinki and HUS Helsinki University Hospital, Helsinki, Finland</footnote>

<footnote n="6"><sup>6</sup>Department of Oral Pathology and Radiology, University of Turku, Turku, Finland</footnote>

<footnote>This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.</footnote>

<footnote>© 2024 The Authors. *Cancer Medicine* published by John Wiley & Sons Ltd.</footnote>
# 1 | INTRODUCTION

Pancreatic ductal adenocarcinoma (PDAC), the most prevalent pancreatic cancer,¹ has a poor prognosis and low survival rate (less than 10% 5-year survival rate). Incidences are expected to rise further in the near future with projections indicating a more than two-fold increase in the number of cases within the next decade.²⁻⁴ Low survival rates are caused by delayed diagnosis which results in initiation of treatment at an advanced stage when tumor cells have already started to invade surrounding tissues and metastasize.⁵,⁶ A prominent reason for delayed diagnostics is the lack of PDAC biomarkers, making it a priority to discover novel genetic and biological features of PDAC to facilitate diagnosis and development of novel treatments.

Recent studies relying on data derived from animal models have shown that embryonic pancreas and PDAC development share molecular similarities. Developmental regulators normally prevent apoptosis and ensure cell proliferation of progenitor cells, however, when mis-expressed in adult cells, may cause changes in cell plasticity and promote survival of tumor cells thereby contributing to tumor formation.⁷ For example, PDX1 is essential for pancreas formation⁸ and maintenance of adult endocrine cell function,⁹ however aberrant PDX1 expression in adult exocrine cells promotes formation of pancreatic neoplastic lesions and epithelial to mesenchymal transition.¹⁰ Another example is the developmental transcription factor HNF1A which is highly enriched in pancreatic cancer stem cells and HNF1A overexpression promotes tumor formation and invasiveness.¹¹ These examples illustrate that exocrine and ductal cells undergo dedifferentiation processes similar to development when becoming carcinogenic.

Human and mouse pancreas development differ significantly, restricting the possibilities of identifying common molecular pathways between human pancreas development and tumor formation using animal models. To identify novel signaling pathways in PDAC, a comparison of RNA sequencing data from human embryonic and adult pancreas with PDAC samples was performed. In this analysis, the transcriptional regulators E2F1 and E2F8 were found to be enriched in embryonic and PDAC transcriptomes.

Previous studies have linked the E2F family of transcription factors to cell cycle regulation where E2F1 mediates both cell proliferation and TP53/p53-dependent apoptosis,¹² whereas E2F1's role in PDAC has not been studied thoroughly, with only a link to PDAC chemoresistance being reported.¹³ Moreover, the function of E2F7 and E2F8 which represent a unique repressive arm of the E2F transcriptional network and control the E2F1-p53

apoptotic axis by directly binding to and repressing E2F1 transcription,¹⁴ has not been studied in PDAC so far.

Single cell sequencing analysis of human embryonic pancreas and PDAC showed that E2F1 and E2F8 were expressed in proliferating pancreatic progenitor cells and ductal cancer cells and co-expression analysis revealed that E2F expression was highly correlated with cell cycle regulation in both tissues. Functional analysis in cancer cells demonstrated that loss of E2F factors resulted in reduced colony formation, tumor cell spreading and migration, proliferation, and oncogene expression, supporting a role of these transcriptional regulators in contributing to carcinogenesis.

# 2 | MATERIALS AND METHODS

## 2.1 | Gene expression profiling in human fetal tissues

Embryonic pancreas was obtained from terminated fetuses (7–14 gestational weeks, *n* = 16). Informed consent was obtained from the participating women. Ethical permission has been obtained from the regional ethics committee in Lund (Dnr 2012/593, Dnr 2015/241, Dnr 2018-579).

Bulk RNAseq: RNA was extracted and libraries generated as previously described.¹⁵ For comparison with adult tissues, the RNA-Seq expression quantification pipeline from GTEX V8 (https://github.com/broadinstitute/gtex-pipeline/) was used.¹⁶,¹⁷ Paired-end 101 bp-long reads were aligned to the Reference Human Genome Build 38 using STAR v2.5.3a,¹⁸ annotation Gencode v26. After alignment and post-processing, expression quantification was performed using RSEM and RNASeqQC.¹⁹,²⁰ Normalization of library sizes was performed in edgeR by dividing counts by the library counts sum and multiplying the results by a million to obtain counts per million (CPM) values.

Fetal-specific expression was defined as ≥1CPM in >75% of the samples. If a gene showed expression in fetal but ≤1CPM in adult tissues, the gene was considered to have fetal specific expression.

Comparison with expression in adult pancreas: edge-R was used to perform differential expression analysis with age and sex as covariates. Batch correction was performed by considering platform and batches using COMBAT.

Single-cell RNAseq of embryonic pancreas has been performed as previously described.¹⁵ Unsupervised clustering of 3199 cells was performed and clusters were assigned using top 20 expressed genes (listed in Ref. [15]). Normalized counts were obtained, and correlation of E2Fs genes with all the genes expressed in at least three cells was performed using spearman correlation
implemented in R. Significant correlated genes with a false discovery rate (FDR) lower than 5% were then studied using Reactome analysis tool to identify relevant pathways.²¹

## 2.2 | Study population

This study is comprised of a cohort of 154 patients with pancreatic adenocarcinoma operated on without neoadjuvant therapy between 2000 and 2013 at the Department of Surgery, Helsinki University Hospital (HUH). As previously described,²² clinical and survival data were compiled from patient records and the Finnish Population Register Centre. Median age of patients at time of diagnosis was 64.8 (interquartile range 59.1–71.0) and median duration of overall survival was 2.0 years (range 0.0–13.1). The study protocol was approved by the Surgical Ethics Committee of HUH (Dnr HUH 226/E6/06, extension TKM02 §66 17.4.2013). Tissue archive samples examination was done with permission from the Finnish Medicines Agency (Dnr FIMEA/2021/006901 28.12.2021). Written informed consent was obtained upon inclusion in the study. Patient information, samples and data are handled and stored in accordance with the Declaration of Helsinki and other local regulations.

## 2.3 | Tumor tissue microarray immunohistochemistry

Tumor tissue microarrays (TMA) blocks were prepared, cut and stained as previously described.²² Once diagnosis confirmed by an experienced pathologist, immunohistochemical staining were performed on TMA slides as previously described.²² TMA slides were incubated with E2F1 (VD301790, Thermo Fischer, 1:100 diluted in Dako REAL Antibody Diluent S2022) or E2F8 antibody (ab185727, Abcam, 1:500 diluted in Dako REAL Antibody Diluent S2022) overnight or for 60 min, respectively. Brain and thyroid tissue served as positive control.

## 2.4 | Scoring of E2F staining in tumors

Immunoreactivity of tumor cells was evaluated independently by two researchers (T.K. and J.H) as previously described²² and scored separately for E2F1 and E2F8 as follows: 0, negative; 1, weakly positive; 2, moderately positive and 3, strongly positive staining. For statistical analysis, samples were grouped into low (negative—low), moderate and high. Survival analysis

was done using the Kaplan-Meier method and compared by log-rank test as previously described.²²

## 2.5 | Gene expression studies in PDAC samples

Sequencing data consisting of PDAC samples (n=24) and control samples (n=11) from Peng et al.²³ were used (accession number: CRA001160) and processed as closely as possible to the parameters described previously.²³ Ten cell types were identified and categorized among the PDAC samples separately according to the criteria used by Peng et al.²³ Clusters were assigned using top 100 differentially expressed genes.²³ These cells were re-clustered in Loupe browser with default settings to get a tSNE plot, which resulted in a total of 57,383 re-clustered cells (PDAC=41,862 and control=15,521). E2F expression was then filtered for, and plots were generated using the built-in function in Loupe browser. Expression correlation of E2F1 and E2F8 genes was performed as described in the “Gene expression profiling in human fetal tissues” paragraph.

## 2.6 | Target gene analyses

Potential E2F1 or E2F8 target genes were identified using ChIP-Atlas (http://chip-atlas.org<sup>24,25</sup>) with binding sites located around transcription start sites (±1kb) being considered. Potential target genes were classified according to their average binding scores of MACS2 and compared with the list of genes significantly co-expressed with E2F1 or E2F8 in PDAC tissues (“Gene expression studies in PDAC samples”). Overall survival and disease-free survival maps of identified target genes were generated using the Gene Expression Profiling Interactive Analysis (GEPIA) database (http://gepia2.cancer-pku.cn/index.html<sup>25</sup>).

## 2.7 | Cells

Human pancreatic cancer cell lines PANC-1<sup>26</sup> were obtained from the European Collection of Authenticated Cell Cultures (ECACC, 87092802) and were maintained in DMEM high glucose medium (Sigma, D6429) supplemented with 10% FBS and 1% penicillin-streptomycin. Cells were incubated at 37°C in a humidified incubator with 5% CO<sub>2</sub> and frequently tested for mycoplasma contamination using MycoAlert™ detection kit (Lonza, LT07-118).

Cells were transfected with negative control siRNA (Invitrogen, 4390843), E2F1 siRNA (Invitrogen, s4406) or
E2F8 siRNA (Invitrogen, s36210) at 10 nM with RNAiMax Lipofectamine (Invitrogen, 13778150) according to the manufacturer's guidelines and cultured for 48 h or 7 days.

## 2.8 | qPCR

Total RNA was extracted using RNeasy Qiagen kit and cDNA was generated. qPCR assays were performed using StepOnePlus Real-Time PCR System. Relative gene abundance was calculated using the ΔΔCt method with TBP, S18, and GAPDH as reference genes and expressed as fold change to control (Table S1).

## 2.9 | Western blot

Transfected PANC-1 cells were lysed using RIPA buffer with protease and phosphatase inhibitors. Whole-cell extracts (20 µg) were separated in Tris-glycine sodium dodecyl sulphate (SDS) gels and transferred onto polyvinylidenedifluoride (PVDF) membranes. E2F1 (1:300, Santa Cruz biotechnology, sc-56662), E2F8 (1:300, Santa Cruz biotechnology, sc-514064), Caspase 8 (1:1000, Cell signaling technology, CST-9746) and Caspase 9 (1:1000; Cell signaling technology, CST-9508) primary antibodies were used. Band density for target proteins was normalized against total protein load and analyzed using Image Lab software.

## 2.10 | Colony assay

To detect colony formation, cells were transfected in 24-well plates at a density of 600 cells/well. Medium was changed every 2 days. After 7 days of incubation, cells were fixed with paraformaldehyde (PFA) 4% for 15 min and stained with a 0.1% crystal solution in 20% EtOH for 30 min. Once dried, colonies were photographed and pictures were analyzed by ColonyArea plugin²⁷ with ImageJ software. Colony formation was quantified by determining the percentage of area covered by cell colonies and the cell density according to the intensity of the staining.

## 2.11 | Proliferation and apoptosis

Proliferation was quantified using Click-iT™ Plus EdU Cell Proliferation Kit for Imaging, Alexa Fluor™ 647 dye (Thermofisher, C10640). Cells were transfected with siRNAs in eight-well chamber slides at a density of 1 × 10⁴ cells per well (48 h) or 3.75 × 10³ cells per well (7 days). After

48 h or 7 days of incubation, cells were exposed to 10 µM of 5-ethynyl-2'-deoxyuridine (EdU) for 4 h at 37°C and then fixed with PFA 4%. EdU labelling was done according to manufacturer's instructions (Life Technologies, C10640) and nuclei were stained with Hoechst 33342. Images were acquired using a confocal microscope and quantification was done using the Cell Counter Plugin in ImageJ to manually tag and count stained cells for each color channel and apoptotic bodies.

## 2.12 | Triplex assay

PANC-1 cells were transfected with siRNAs in 96-well plate at a density of 2000 cells per well (48 h) or 500 cells per well (7 days). Cell viability, cytotoxicity and apoptosis were measured 48 h and 7 days after the knock-down using ApoTox-Glo triplex assay kit (Promega, G6320) according to manufacturer's instructions.

## 2.13 | Migration

Cell culture inserts in 24-well plate (8 µm pore size; Falcon, 353097) were used for the migration assay. Seven days after transfection, cells were trypsinized and counted. Ten percentage FBS medium was used as chemoattractant in the bottom chamber and 4 × 10⁴ cells diluted in 0.5% FBS medium were added in the top chamber. After 7 h of incubation at 37°C, cells were fixed and stained with eosin. Cells that didn't migrate were removed from the insert using a cotton swab. Insert membrane was cut and mounted and migrated cells were imaged and counted.

## 2.14 | Spreading

Seven days after transfection cells were trypsinized, and 4 × 10⁴ cells diluted in 0.5% FBS medium were cultured on coverslips coated with fibronectin (Sigma, F1141) in a 12-well plate as described by Pijuan.²⁸ One hour later, cells were fixed with PFA 4% and stained with a 0.1% crystal solution.

According to their shape, cells were counted as spread (visible cytoplasm surrounding the nucleus) or unspread (round cells). Ratio was calculated as the number of spread cells divided by the total number of cells on the coverslips.

## 2.15 | Statistics

Statistical analyses and graphs were performed using R software with ggplot2 and ggpubr packages. Each
experiment has been performed at least two times independently, with at least three replicates each time. Multiple comparisons between control and different siRNAs were analyzed using ANOVA with Tukey's post-hoc test. Normality and homoscedasticity of the variances were verified by checking graphically (i.e., Q-Q plot and residuals versus fits plot). Data were visualized using box plots, the central mark indicates the median and edges indicate 25th and 75th percentiles. Whiskers extend to the largest or smallest point comprised within 1.5× of the interquartile range from both edges.

3 | RESULTS

3.1 | E2F1 and E2F8 are expressed in the
developing human pancreas and PDAC

Expression of developmental regulators is re-initiated dur-
ing tumor formation⁷ and altered expression is associated
with malignancy. To identify key genes in PDAC develop-
ment, we compared expression of genes upregulated in
PDAC²⁹,³⁰ with expression data from human embryonic
pancreas we obtained previously.¹⁵ E2F1, E2F3, and E2F8
genes were expressed both in PDAC and in human embry-
onic pancreas anlagen (Figure 1A). From this list, E2F1
and E2F8 expression was significantly higher in embryonic
than in adult pancreas (Figure 1A). Thus, E2F1 and E2F8
expression are up-regulated in PDAC and human embry-
onic pancreas compared to adult tissues.

Using open-access databases,³¹⁻³³ we confirmed the higher expression levels of E2F1 and E2F8 in PDAC compared to adult pancreas (Table S2). Protein expression was detected only in PDAC samples. To complete these observations, we stained and scored E2F1 and E2F8 protein expression in 154 tumor samples. TMA analysis detected E2F1 and E2F8 protein expression in 83.8 and 90.3% of analyzed PDAC samples, respectively. Interestingly expression was detected mainly in the cytoplasm (E2F1) and both cytoplasm and nucleus (E2F8) of tumor cells (Figure 1B).

Kaplan–Meier analysis revealed that E2F8 expres-
sion associates significantly with DSS, with the poorer
outcome for patients with low E2F8 expression com-
pared to moderate (p = 0.015) or high E2F8 expression

(p = 0.004). No difference was seen between patients with moderate or high expression. Five-year DSS for patients with a low E2F8 expression was 9.1% (95% CI 0%–21.1%) compared to 21.2% (95% CI 11.8%–30.6%) among those with a moderate or 19.9% (95% CI 8.7%–31.1%) with a high expression. Low E2F1 expression associates with worse prognosis compared high expres- sion (p = 0.028), but no difference was seen compared to moderate expression (p = 0.424). No difference with statistical significance was seen between moderate and high expression (p = 0.066). Five-year DSS for patients with a low E2F1 expression was 16.7% (95% CI 7.3%–26.1%) compared to 19.0% (95% CI 9.0%–28.0%) among those with a moderate or 35.7% (95% CI 10.6%–60.8%) with a high expression. (Figure 1C,D). No significant association with patho-clinical parameters like age, sex, T-or N-status, stage, grade, lymph node ratio or peri- vascular/perineural invasion were detected (Table S3). These data suggest that variations in E2F expression in- fluence tumor malignancy.

To assess which type of embryonic and PDAC cells expressed E2F1 and E2F8, single cell RNA sequencing analysis was performed. In human embryonic pancreas E2F1 and E2F8 expression was predominantly detected in highly proliferating pancreatic progenitor cell populations (Figure 1E; Figure S1). In PDAC, E2F1, and E2F8 expression was detected in immune and type 2 ductal cells (Figure 1F; Figure S2A, S2B). Type 2 ductal cells are considered highly proliferative and malignant in contrast to type 1 cells²³ and in this cell population, only 9% of the E2F1⁺ or E2F8⁺ cells co-expressed E2F1 and E2F8 (Figure S2C). These findings show that E2F1 and E2F8 mRNA expression is present during normal pancreas development and re-initiated in malignant ductal cells with aberrant expression impacting patient survival suggesting an involvement in PDAC tumor formation and maintenance.

3.2 | Loss of E2F1 and E2F8 gene expression reduces clonogenic capacity and cell proliferation

To determine if E2F genes participate in tumor biology,
E2F1 and E2F8 expression was knocked down using

FIGURE 1 E2F gene expression in pancreas from fetal, adult, and cancer tissues. (A) Plot of E2F gene expression in fetal (7–14 weeks post-conception (PC)) and adult pancreas. Y-axis shows expression in LogCPM. Significant expression difference between fetal and adult tissues, fetal specificity (at least one CPM in at least 75% of the fetal samples) and E2F genes up regulated in PDAC (both in Mao's and Liu's studies<sup>27,28</sup>) are indicated for each E2F transcription factor. (B) Representative images of E2F1 and E2F8 immunohistochemical staining and distribution of TMA-samples according to staining intensity. Original magnification: 200×. Kaplan–Meier disease-specific survival (DSS) curves for PDACs according to E2F1 (C) or E2F8 (D) protein expression. (E, F) Associated expression patterns of E2F1 and E2F8 in developing pancreas from an 8-week PC embryo, in normal adult pancreas and PDAC tissues. Expression of E2F1 and E2F8 on the t-SNE embedding from Figures S1 (E) and S2 (F). Any cell with E2F1 or E2F8 detected was colored according to its expression level.
(A) E2F expression in fetal and adult pancreas

<figure><img src="image_4.png" /></figure>

(C) Survival curve - E2F1 protein expression

<figure><img src="image_5.png" /></figure>

(E) Fetal pancreas

<figure><img src="image_6.png" /></figure>

(B) E2F1 and E2F8 protein expression in PDAC

<figure><img src="image_7.png" /></figure>

<table><thead><tr><th>TMA-samples (n=154)</th><th>E2F1</th><th>E2F8</th></tr></thead><tbody><tr><td>negative</td><td>25 (16.2%)</td><td>15 (9.7%)</td></tr><tr><td>weak</td><td>36 (23.4%)</td><td>7 (4.5%)</td></tr><tr><td>moderate</td><td>79 (51.3%)</td><td>77 (50%)</td></tr><tr><td>strong</td><td>14 (9.1%)</td><td>55 (35.7%)</td></tr></tbody></table>

(D) Survival curve - E2F8 protein expression

<figure><img src="image_8.png" /></figure>

(F) Adult pancreas

<figure><img src="image_9.png" /></figure>

PDAC tissues

<figure><img src="image_10.png" /></figure>
siRNAs in the pancreatic cancer cell line PANC-1 (i.e., siE2F1 and siE2F8). This human pancreatic cancer cell line has been isolated from a pancreatic carcinoma of ductal cell origin and highly expresses E2F candidate genes.²³ Knock down (KD) efficiency of E2F1 and E2F8 siRNAs was confirmed by RT-qPCR and western blot. The transfection efficiencies at 48 h and 7 days after KDs were > 75% and 60% respectively (Figure S3).

To assess if siE2F1 and siE2F8 affect the capacity of tumor cells to grow we first examined the ability of KD cells to form colonies. Seven days after the transfection, the area and density of PANC-1 colonies were significantly decreased in siE2F1 and siE2F8 cells compared to the control. In mean, colony formation was reduced over 95% in siE2F1 and 70% in siE2F8 compared to control (area and density *p*-values <1e-8 for both KDs Figure 2A,B).

We next quantified the proliferative capacity of KD cells using an EdU incorporation assay (Figure 2C). Forty-eight hour after transfection, the proportion of EdU-positive cells was significantly decreased in siE2F1 (–20%, *p*-value 0.0052) and siE2F8 (–54%, *p*-value <1e-8) compared to the control (Figure 2D). Cell proliferation capacities were also tested 7 days post transfection to detect long term effects of siE2F1 and siE2F8. Seven days after the transfection, the proliferation capacity was still decreased in KD cells to 29% in siE2F1 (*p*-value 0.0055) and 26% in siE2F8 (*p*-value 0.019) compared to control cells. Taken together, these results suggest that E2F transcription factors regulate the cell cycle in PDAC cells.

## 3.3 | Loss of E2F1 expression affects cell viability and apoptosis

siE2F1 affected the clonogenic capacity to a greater extent than siE2F8 (*p*-value siE2F1 vs. siE2F8: 0.005). To assess if this is due to an effect of E2F function on cell viability, we determined apoptosis and cell viability rate of siE2F1 and siE2F8 using Triplex assay and quantifying the apoptotic bodies. Forty-eight hour after transfection, no significant differences in viability and apoptosis were observed between the KDs and control cells (Figure 3A–D). Seven-days post transfection, siE2F1 significantly reduced cell viability (–60%, *p*-value <1e-8), increased caspase activity (*p*-value 0.0001) and number of apoptotic bodies (*p*-value 3e-7) compared to control (Figure S4), while

siE2F8 did not affect cell caspase activity and cell viability (Figure 3D–G). Thus, low proliferation rate and increased apoptosis contribute to the decreased clonogenic capacity of siE2F1 cells. These results suggest that reduced E2F1 but not E2F8 expression promotes apoptosis.

## 3.4 | Loss of E2F1 and E2F8 impairs cell spreading and migration

Cell migration capacity is critical for tumor invasiveness. To determine if E2F1 and E2F8 regulate migration of tumor cells, cell spreading and migration assays were performed 7 days after transfection.

To study cell spreading, coverslips were coated with fibronectin, a main constituent of the tumor stroma³⁴ and a poor prognosis marker in PDAC.³⁵ We observed that cell spreading was reduced by 78% and 29% in siE2F1 and siE2F8 cells, respectively (Figure 4A,B, *p*-values E2F1 <1e-8, E2F8 2.67e-05) while 0.5% FBS medium was sufficient to promotes cell adhesion in control cells.

Transwell assays were conducted to determine the migration potentials of the transfected pancreatic cancer cells. After 7 h of incubation, control cells effectively migrated upon stimulation with chemo-attractants, while siE2F1 and siE2F8 cells had significantly reduced cell migration (Figure 4C,D, *p*-values siE2F1 <1e-8, siE2F8 0,001). Thus, reduced expression of E2F1 and E2F8 reduced migration and the metastatic potential of PANC-1 cells.

## 3.5 | E2F transcription factors are co-expressed with genes regulating the cell cycle

The functional studies described above indicate that E2F1 and E2F8 are critical for tumor cell survival, proliferation, and migration. To identify pathways potentially regulated by E2F transcription factors and involved in carcinogenesis, unbiased co-expression analyses of E2F1 and E2F8 in single PDAC, adult and embryonic pancreatic cells were performed (Tables S4–S8).

In fetal cells, E2F8 expression significantly correlates with 2297 genes, while only one gene correlates with E2F1 expression, probably due to the lower number of E2F1⁺ cells (Tables S4 and S5). In adult, 211 genes correlate with

**FIGURE 2** Colony formation and proliferation capacities of E2F transfected cells. Relative colony formation quantification (A) and morphological observation (B). (A) Results are shown as colony area or staining intensity percent 7 days after the knock-down compared to the negative control (transfected with scrambled RNA). (B) Colonies stained with crystal violet. (C) Dapi and EdU staining 48 h or 7 days after transfection. Cells were exposed to EdU for 4 h. (D) Results are shown as EdU positive cells to the total number of cells. (A, D) Measurements were derived from two independent experiments of three replicates per condition. Tukey's post-hoc test significances are indicated by stars compared to the control when significant, \* *p* < 0.05, \*\* *p* < 0.01, \*\*\* *p* < 0.001.
E2F1 expression (Table S6) and no gene expression correlations were obtained for E2F8 due to the low number of E2F8⁺ cells. Finally, in PDAC, unbiased co-expression

analysis showed that E2F1 and E2F8 expression significantly correlates with 4989 and 612 genes, respectively (Tables S7 and S8).

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>
<figure><img src="image_15.png" /><figcaption>FIGURE 3 Cell viability and apoptosis 48 h (A–D) and 7 days (D–G) after transfection. Cell viability 48 h (A) or 7 days (E) after transfection. Results are shown as viability/cytotoxicity ratio to the negative control. Relative caspase activity 48 h (B) or 7 days (F) after transfection. Results are shown as caspase-3/7 activity to the negative control. Quantification of the apoptotic bodies 48 h (C) or 7 days (G) after transfection. Results are shown as counted apoptotic bodies to the total number of cells. At least 1000 cells were quantified per condition. (D) Dapi staining 48 h or 7 days after transfection, red arrows show apoptotic bodies, scale bars: 50 µm. (A–G) Measurements were derived from two independent experiments of three replicates. (A–C, E–G) Tukey's post-hoc test significances are indicated by stars compared to the control when significant, *** p < 0.001.</figcaption></figure>

Only 17% of the genes significantly co-expressed with E2F1 in adult were also significantly correlated in PDAC, suggesting different target genes in normal pancreas and PDAC cells. In contrast, 51% of the genes that had significant correlation of gene expression with E2F8 in PDAC were also correlated in embryonic cells, and 56% of the genes that shared expression correlation to E2F8 in the embryo were significantly correlated to E2F1 in PDAC. These data suggest that E2F1 and E2F8 mainly act in the same gene

networks during embryogenesis and PDAC. Interestingly, E2F8 correlates with the atypical repressor E2F7 in the embryo (*p*-value 0.0016) whereas it correlates with the activators E2F1 and E2F2 in PDAC (*p*-values 0.0006 and 5e-06, respectively).

In PDAC, 83% of the genes with significant positive correlation with E2F8 expression were also positively correlated with E2F1 expression, which may suggest redundant function between E2F1 and E2F8. Using
<figure><img src="image_16.png" /><figcaption>**FIGURE 4** Spreading and migration capacities of PANC-1 cells 7 days after transfection. Plot showing the percentage of spreading cells (A). Cells spread on fibronectin coated plates (B). Spreading cells are indicated by a black arrow, round cells were considered to not have spread (B). Number of migrated cells per image after 7 h of incubation (C) and eosin saining of the migrated cells (D). (A, C) At least 1000 cells were quantified per condition. Scale bars: 50 µm. Measurements were derived from two independent experiments of three replicates per condition. Tukey's post-hoc test significances are indicated by stars compared to the control when significant, *** *p* < 0.001.</figcaption></figure>

Reactome analysis tools, we determined that genes co-expressed with E2F1 or E2F78 are involved in pathways regulating cell cycle and proliferation but also RNA and protein processing and translation (Table 1). These results suggest a central role of E2F1 and E2F8 in cancer cell biology.

In addition, E2F1 expression specifically correlates with genes involved in the regulation of apoptosis and mitochondrial translation pathways (Table 2), suggesting a more prominent role in PDAC regulation.

## 3.6 | E2F1 and E2F8 transcription factors can directly regulate known oncogenes

A comparison of genes co-expressed with E2F1 or E2F8 in PDAC (Tables S7 and S8) with published chromatin immunoprecipitation sequencing (ChiP-seq) results (ChIP-Atlas) was performed to obtain a list of potential direct E2F target genes in E2F+ PDAC cells. Once classified by Reactome pathways (Tables S9 and S10), we observed that most of these genes were involved in pathways regulating carcinogenesis like cell cycle, DNA replication, gene expression and programmed cell death.

A thorough literature review and GEPIA2 analysis showed that most of these genes were already known as potential oncogene and/or linked to poor survival in PDAC but few of them have been shown to be regulated by E2F transcription factors (Figure 5A; Tables S11 and S12).

To validate our analysis, we selected 25 oncogenes identified as potential E2F1 and/or E2F8 targets (Figure 5B) and assessed if their expression was impaired in E2F KD cells. Forty-eight hour post siE2F1 or siE2F8 transfection, all the selected genes were significantly downregulated in siE2F1 and/or siE2F8 cells compared to controls (Figure 5C). These results suggest that E2F1 and E2F8 directly regulate the expression of numerous known oncogenes in PDAC tumors and thus act as master regulators of pancreatic cancer biology.

## 4 | DISCUSSION

Great progress has been made in cancer treatment with overall survival rates increasing for many different cancer types.36 However, the diagnosis and treatment of pancreatic cancer has not been improved significantly during the past years. Immuno-and chemotherapy have limited success and only after total resection of tumor tissue. Thus, great efforts have been made to identify novel biomarkers and drugable targets. This may be facilitated by novel expression data sets from PDAC tissue. Embryonic progenitor and PDAC cells share many features like cell plasticity, survival, and proliferation potential suggesting that comparison between these cell types may lead to the identification of important genes for pancreatic tumorigenesis. Here we compared single and bulk gene expression from PDAC with expression data from human embryonic pancreas to identify novel PDAC genes that are present in embryonic pancreatic progenitor and tumor cells. Using this expression correlation approach, we identified E2F genes to be highly expressed in embryonic pancreas and PDAC.

In the developing pancreas, E2F1 and E2F8 are expressed in proliferating pancreatic progenitor cells
<table><caption>TABLE 1 Most relevant pathways identified by Reactome using significant co-expressed genes with E2F1 or E2F8 in PDAC or embryo development. All the entities have a false discovery rate (FDR) < 0.05.</caption><thead><tr><th rowspan="2">Most relevant pathways sorted by p-value</th><th colspan="6">Number of entities co-expressed with</th></tr><tr><th>E2F1 in PDAC</th><th>p value</th><th>E2F8 in PDAC</th><th>p value</th><th>E2F8 in embryo</th><th>p value</th><th>Entities in the pathway</th></tr></thead><tbody><tr><td>Formation of a pool of free 40S subunits</td><td>100</td><td>1.11e-16</td><td>53</td><td>1.11e-16</td><td>75</td><td>1.11e-16</td><td>106</td></tr><tr><td>GTP hydrolysis and joining of the 60S ribosomal subunit</td><td>110</td><td>1.11e-16</td><td>54</td><td>1.11e-16</td><td>79</td><td>1.11e-16</td><td>120</td></tr><tr><td>SRP-dependent cotranslational protein targeting to membrane</td><td>109</td><td>1.11e-16</td><td>54</td><td>1.11e-16</td><td>79</td><td>1.11e-16</td><td>119</td></tr><tr><td>L13a-mediated translational silencing of ceruloplasmin expression</td><td>109</td><td>1.11e-16</td><td>55</td><td>1.11e-16</td><td>78</td><td>1.11e-16</td><td>120</td></tr><tr><td>Nonsense-mediated decay (NMD)</td><td>110</td><td>1.11e-16</td><td>53</td><td>1.11e-16</td><td>80</td><td>1.11e-16</td><td>124</td></tr><tr><td>Nonsense mediated decay (NMD) enhanced by the Exon junction complex (EJC)</td><td>110</td><td>1.11e-16</td><td>53</td><td>1.11e-16</td><td>80</td><td>1.11e-16</td><td>124</td></tr><tr><td>Eukaryotic translation initiation</td><td>115</td><td>1.11e-16</td><td>55</td><td>1.11e-16</td><td>80</td><td>1.11e-16</td><td>130</td></tr><tr><td>Cap-dependent translation initiation</td><td>115</td><td>1.11e-16</td><td>55</td><td>1.11e-16</td><td>80</td><td>1.11e-16</td><td>130</td></tr><tr><td>DNA replication</td><td>137</td><td>1.11e-16</td><td>43</td><td>1.11e-16</td><td>80</td><td>2.66e-15</td><td>169</td></tr><tr><td>Major pathway of rRNA processing in the nucleolus and cytosol</td><td>153</td><td>1.11e-16</td><td>59</td><td>1.11e-16</td><td>89</td><td>1.11e-16</td><td>189</td></tr><tr><td>M phase</td><td>273</td><td>1.11e-16</td><td>69</td><td>1.11e-16</td><td>175</td><td>1.11e-16</td><td>416</td></tr><tr><td>Cell cycle mitotic</td><td>392</td><td>1.11e-16</td><td>129</td><td>1.11e-16</td><td>256</td><td>1.11e-16</td><td>596</td></tr><tr><td>G1/S transition</td><td>124</td><td>1.11e-16</td><td>42</td><td>1.11e-16</td><td>72</td><td>1.89e-15</td><td>150</td></tr><tr><td>Cell cycle</td><td>485</td><td>1.11e-16</td><td>148</td><td>1.11e-16</td><td>314</td><td>1.11e-16</td><td>734</td></tr><tr><td>Cell cycle checkpoints</td><td>214</td><td>1.11e-16</td><td>57</td><td>1.11e-16</td><td>140</td><td>1.11e-16</td><td>279</td></tr><tr><td>Metabolism of RNA</td><td>567</td><td>1.11e-16</td><td>130</td><td>1.11e-16</td><td>309</td><td>1.11e-16</td><td>829</td></tr><tr><td>mRNA splicing—major pathway</td><td>173</td><td>1.11e-16</td><td>39</td><td>2.86e-12</td><td>112</td><td>1.11e-16</td><td>213</td></tr><tr><td>Processing of capped intron containing pre-mRNA</td><td>237</td><td>1.11e-16</td><td>51</td><td>1.10e-14</td><td>139</td><td>1.11e-16</td><td>299</td></tr><tr><td>mRNA splicing</td><td>177</td><td>1.11e-16</td><td>40</td><td>3.03e-12</td><td>113</td><td>1.11e-16</td><td>224</td></tr></tbody></table>

(Continues)
<table><caption>TABLE 1 (Continued)</caption><thead><tr><th rowspan="2">Most relevant pathways sorted by p-value</th><th colspan="6">Number of entities co-expressed with</th></tr><tr><th>E2F1 in PDAC</th><th>p value</th><th>E2F8 in PDAC</th><th>p value</th><th>E2F8 in embryo</th><th>p value</th><th>Entities in the pathway</th></tr></thead><tbody><tr><td>rRNA processing in the nucleus and cytosol</td><td>162</td><td>1.11e-16</td><td>60</td><td>1.11e-16</td><td>92</td><td>1.11e-16</td><td>208</td></tr><tr><td>Mitotic metaphase and anaphase</td><td>177</td><td>1.11e-16</td><td>43</td><td>1.35e-12</td><td>106</td><td>1.11e-16</td><td>250</td></tr><tr><td>Mitotic anaphase</td><td>176</td><td>1.11e-16</td><td>42</td><td>1.35e-12</td><td>105</td><td>2.22e-16</td><td>249</td></tr><tr><td>rRNA processing</td><td>178</td><td>1.11e-16</td><td>64</td><td>1.11e-16</td><td>104</td><td>1.11e-16</td><td>246</td></tr><tr><td>Translation</td><td>261</td><td>1.11e-16</td><td>71</td><td>1.11e-16</td><td>124</td><td>5.55e-16</td><td>339</td></tr><tr><td>G2/M checkpoints</td><td>127</td><td>1.11e-16</td><td>32</td><td>2.85e-12</td><td>80</td><td>1.11e-16</td><td>153</td></tr><tr><td>Viral mRNA translation</td><td>88</td><td>2.97e-11</td><td>52</td><td>1.11e-16</td><td>70</td><td>1.11e-16</td><td>114</td></tr><tr><td>Nonsense mediated decay (NMD) independent of the exon junction complex (EJC)</td><td>94</td><td>4.44e-16</td><td>53</td><td>1.11e-16</td><td>71</td><td>1.11e-16</td><td>101</td></tr><tr><td>Translation initiation complex formation</td><td>56</td><td>5.50e-10</td><td>27</td><td>1.11e-16</td><td>39</td><td>4.93e-12</td><td>62</td></tr><tr><td>Mitotic G1 phase and G1/S transition</td><td>136</td><td>1.11e-16</td><td>51</td><td>1.11e-16</td><td>80</td><td>1.89e-15</td><td>174</td></tr><tr><td>Peptide chain elongation</td><td>89</td><td>2.44e-15</td><td>52</td><td>1.11e-16</td><td>70</td><td>1.11e-16</td><td>97</td></tr><tr><td>SARS-CoV-1 modulates host translation machinery</td><td>36</td><td>1.11e-06</td><td>25</td><td>1.11e-16</td><td>32</td><td>1.81e-12</td><td>41</td></tr><tr><td>Eukaryotic translation elongation</td><td>92</td><td>2.22e-15</td><td>53</td><td>1.11e-16</td><td>72</td><td>1.11e-16</td><td>102</td></tr><tr><td>rRNA processing in the nucleus and cytosol</td><td>162</td><td>1.11e-16</td><td>60</td><td>1.11e-16</td><td>92</td><td>1.11e-16</td><td>208</td></tr><tr><td>Eukaryotic translation termination</td><td>92</td><td>5.01e-14</td><td>52</td><td>1.11e-16</td><td>70</td><td>1.11e-16</td><td>106</td></tr><tr><td>Mitotic prometaphase</td><td>132</td><td>1.80e-09</td><td>38</td><td>9.00e-12</td><td>100</td><td>1.11e-16</td><td>211</td></tr></tbody></table>

suggesting a role in cell cycle regulation in human embryos. This is consistent with previous findings showing that E2F1 and E2F8 are present in developing Pdx1+ pancreatic progenitor cells³⁷,³⁸ and that E2F1 is essential for regulating proliferation in the mouse pancreas anlage.³⁸ E2F1 and E2F8 expression is also enriched in PDAC and the majority of E2F1 and E2F8 producing cells are malignant ductal cells suggesting that E2F transcription factors may also participate in

PDAC tumor formation similar to what has been observed in tumorigenesis and treatment resistance of colon, breast, and cervical cancer.³⁹⁻⁴¹ Previous reports showed that high expression of mainly E2F1 is predictive of a poor clinical outcome in PDAC patients,⁴²,⁴³ while our analysis of E2F gene expression in PDAC tissue showed that low E2F8 expression was associated with poor disease specific patient survival. Interestingly, E2F8 expression is mainly detected in the
<table><caption>TABLE 2 Most relevant pathways identified by Reactome using the genes significantly co-expressed only with E2F1 in PDAC.</caption><thead><tr><th rowspan="2">Most relevant pathways sorted by p-value identified using E2F1-only co-expressed entities in PDAC</th><th colspan="4">Number of entities co-expressed with</th><th rowspan="2">Entities in the pathway</th></tr><tr><th>E2F1</th><th>p value</th><th>E2F8</th><th>p value</th></tr></thead><tbody><tr><td>Processing of capped intron-containing pre-mRNA</td><td>237</td><td>1.11e-16</td><td>51</td><td>1.10e-14</td><td>299</td></tr><tr><td>Translation</td><td>261</td><td>1.11e-16</td><td>71</td><td>1.11e-16</td><td>339</td></tr><tr><td>Metabolism of RNA</td><td>567</td><td>1.11e-16</td><td>130</td><td>1.11e-16</td><td>829</td></tr><tr><td>Cell cycle</td><td>485</td><td>1.11e-16</td><td>148</td><td>1.11e-16</td><td>734</td></tr><tr><td>Cell cycle checkpoints</td><td>214</td><td>1.11e-16</td><td>57</td><td>1.11e-16</td><td>279</td></tr><tr><td>mRNA splicing—major pathway</td><td>173</td><td>1.11e-16</td><td>39</td><td>2.86e-12</td><td>213</td></tr><tr><td>mRNA splicing</td><td>177</td><td>1.11e-16</td><td>40</td><td>3.03e-12</td><td>224</td></tr><tr><td><strong>Mitochondrial translation elongation</strong></td><td><strong>81</strong></td><td><strong>9.80e-13</strong></td><td>10</td><td>0.0123<br>FDR 0.053</td><td>94</td></tr><tr><td>M phase</td><td>273</td><td>1.11e-16</td><td>69</td><td>1.11e-16</td><td>416</td></tr><tr><td><strong>Mitochondrial translation termination</strong></td><td><strong>80</strong></td><td><strong>2.45e-12</strong></td><td>10</td><td>0.0123<br>FDR 0.053</td><td>94</td></tr><tr><td><strong>Mitochondrial translation initiation</strong></td><td><strong>81</strong></td><td><strong>2.73e-12</strong></td><td>10</td><td>0.0141<br>FDR 0.056</td><td>96</td></tr><tr><td><strong>Mitochondrial translation</strong></td><td><strong>84</strong></td><td><strong>3.65e-12</strong></td><td>10</td><td>0.0204<br>FDR 0.067</td><td>102</td></tr><tr><td>Separation of sister chromatids</td><td>145</td><td>1.33e-15</td><td>33</td><td>1.13e-09</td><td>195</td></tr><tr><td>DNA replication pre-initiation</td><td>112</td><td>5.55e-16</td><td>29</td><td>4.04e-11</td><td>132</td></tr><tr><td>HIV infection</td><td>168</td><td>1.10e-13</td><td>28</td><td>3.78e-05</td><td>260</td></tr><tr><td>Mitotic anaphase</td><td>176</td><td>1.11e-16</td><td>42</td><td>4.63e-12</td><td>249</td></tr><tr><td>G2/M checkpoints</td><td>127</td><td>1.11e-16</td><td>32</td><td>2.85e-12</td><td>153</td></tr><tr><td>Mitotic metaphase and anaphase</td><td>177</td><td>1.11e-16</td><td>43</td><td>1.35e-12</td><td>250</td></tr><tr><td><strong>Regulation of apoptosis</strong></td><td><strong>51</strong></td><td><strong>1.34e-09</strong></td><td>5</td><td>0.1094<br>FDR 0.215</td><td>54</td></tr><tr><td>Activation of APC/C and APC/C:Cdc20 mediated degradation of mitotic proteins</td><td>69</td><td>2.20e-11</td><td>12</td><td>2.83e-04</td><td>77</td></tr><tr><td>APC/C:Cdc20 mediated degradation of securin</td><td>61</td><td>3.27e-10</td><td>9</td><td>0.0045</td><td>68</td></tr><tr><td>Cell cycle, mitotic</td><td>392</td><td>1.11e-16</td><td>129</td><td>1.11e-16</td><td>596</td></tr><tr><td>APC/C:Cdh1 mediated degradation of Cdc20 and other APC/C:Cdh1 targeted proteins in late mitosis/ early G1</td><td>66</td><td>7.42e-11</td><td>12</td><td>1.98e-04</td><td>74</td></tr><tr><td><strong>Regulation of activated PAK-2p34 by proteasome mediated degradation</strong></td><td><strong>48</strong></td><td><strong>1.37e-09</strong></td><td>5</td><td>0.0810<br>FDR 0.162</td><td>50</td></tr><tr><td>APC/C:Cdc20 mediated degradation of mitotic proteins</td><td>69</td><td>2.20e-11</td><td>12</td><td>2.52e-04</td><td>77</td></tr></tbody></table>

Note: Number of entities were then identified using the full significant gene list for E2F1 or E2F8 in PDAC. In bold, the E2F1-specific pathways. All the entities have a false discovery rate (FDR) < 0.05 if not indicated.

immunogenic PDAC subtype⁴⁴ and correlates with all the signature genes of the classical subtype²⁵,⁴⁵ but not with the genes linked to the basal subtype. These observations suggest that the low E2F8 expressing group

in our survival group may identify patients with basal PDAC which have a much shorter DSS rate, while patients with classical PDAC have moderate to high E2F8 expression.
<figure><img src="image_17.png" /><figcaption>FIGURE 5 Regulation of known oncogenes by E2F1 and E2F8. (A) Disease-free survival and overall survival heatmaps according to the expression of direct E2F target genes in PDAC. The red blocks denote higher risks. The rectangles with frames mean the significant unfavorable results in prognostic analyses. (B) For each oncogene, the lines indicate that the gene is co-expressed with E2F1 (yellow) or E2F8 (orange) in PDAC tissues and contains a direct binding site for E2F1 (yellow) or E2F8 (orange). (C) qPCR analyses showing the expression of 25 oncogenes 48 h after E2F1 or E2F8 KD. Significance was determined for each gene by a one-way ANOVA and Tukey's post-hoc test. NS *p* > 0.05, \* *p* < 0.05, \*\* *p* < 0.01, \*\*\* *p* < 0.001.</figcaption></figure>

E2F1 is a cell cycle regulator which has been associ-
ated indirectly with PDAC by studies describing potential
upstream$^{46,47}$ or downstream$^{48–50}$ target genes that play a
significant role in PDAC biology. By studying direct E2F1
down-regulation, we observed that E2F1 was not only regu-
lating proliferation as proposed previously, but also apoptosis
and cell migration. In normal tissues, a negative feedback
loop regulates E2F signaling, E2F1 inducing expression of
the repressor E2F8.$^{51}$ Few studies followed E2F8 function in
cancer, and so far, no potential E2F8 function in PDAC has
been reported. E2F8 has a dual role in cancer, as it can act
as tumor suppressor$^{52}$ or proto-oncogene$^{40,53}$ depending on
cancer type. Here we showed that in pancreas tumorigenesis,
E2F8 may regulate not only the cell cycle but also migration.

KD of E2F1 and E2F8 in the PDAC cell line PANC-1
resulted in reduced colony formation, with E2F1 having

a significantly stronger effect, suggesting a more promi-
nent role for E2F1 in PDAC. This is consistent with our
findings that E2F1, but not E2F8 KD induced tumor cell
apoptosis. This observation is in contrast to previous
studies which have shown that E2F1 is essential for in-
ducing apoptosis in response to genotoxic treatment in
early tumor cells,⁵⁴ however elevated E2F1 expression
results in resistance to apoptosis in late stage tumor cells
(reviewed in Ref.[55]) which suggests a dual role of this
gene in regulating cell viability/apoptosis of PDAC cells.
PANC-1 cell proliferation was impaired upon E2F1 and
E2F8 KDs suggesting that E2F factors promote cell cycle
regulation in PDAC cells. A finding that is further sup-
ported by expression correlation of E2F1 and E2F8 with
genes involved in cell cycle regulation and mitosis in
PDAC tissues.
No significant correlation between E2F1 or E2F8 with
epithelial–mesenchymal transition (EMT) marker genes
was found and E2F1 and E2F8 KD did not induce mesen-
chymal marker gene expression in PANC-1 cells (data were
not shown) suggesting that the reduced cell spreading/mi-
gration observed upon E2F1 and E2F8 KD is most likely
due to E2F regulating cell proliferation and survival.

E2F1 and E2F8 mainly positively correlate with the
same genes in PDAC tissues, suggesting similar and re-
dundant functions, as observed experimentally in pan-
creatic cancer cells. Our co-expression analyses showed
that the same pathways were correlating with E2F1 and
E2F8 expression in PDAC and embryonic tissues while it
is known that in healthy tissues, E2F8 and E2F1 have an-
tagonistic functions.¹⁴ Interestingly, in embryonic tissues,
E2F8 expression correlates with the atypical repressor
E2F7 whereas in PDAC E2F8 expression correlates with
the transcriptional activators E2F1 and E2F2. Morgunova
et al., showed that both atypical and typical E2Fs can bind
to similar DNA sequences.⁵⁶ Our data suggest that E2F8
may regulate the same pathways during PDAC and em-
bryonic development but acts differently according to its
co-factors.

Links between E2F genes and mitochondrial functions
are now well studied and revealed a role of E2F1 as regulator
of metabolism.⁵⁷⁻⁵⁹ In pancreatic cancer, it has been shown
that E2F1 may induces aerobic glycolysis,⁵⁹ here we link
E2F1 expression in PDAC with expression of genes involved
in the mitochondrial translation machinery. Interestingly,
inhibition of mitochondrial functions is explored as a prom-
ising cancer therapeutic strategy in different tumor types,
and some clinical trials are ongoing.⁶⁰

Using available ChiP-seq data, we identified genes
containing E2F binding sites that are co-expressed with
E2F1 or E2F8 in PDAC. Most of these genes are known
to be over-expressed in PDAC and act in numerous on-
cogenic pathways. For example, ATAD2 expression
is linked to cell invasion, migration, and gemcitabine
resistance⁶¹ and overexpression of TUBA1B, NDC80,
RCF4, and MCM proteins is associated with poorer out-
come.⁶²⁻⁶⁵ Using E2F1 and E2F8 KD cells, we confirmed
that loss of E2F impaired expression of the identified
targets genes. These results suggest that E2F1 and E2F8
over-expression in PDAC induce the expression of nu-
merous oncogenes involved in the regulation of the cell
cycle as PCNA⁶⁶ and CDK1⁶⁷ which affect the patient
prognosis.

5 | CONCLUSION

In summary our findings illustrate that E2F1 and E2F8
are expressed in pancreatic progenitor and PDAC cells,

where they contribute to tumor cell expansion by regu-
lation of cell proliferation, viability and cell migration
making these genes attractive therapeutic targets and
potential prognostic markers for pancreatic cancer.

AUTHOR CONTRIBUTIONS

**Ludivine Bertonnier-Brouty:** Conceptualization (equal); formal analysis (equal); funding acquisition (equal); investigation (equal); methodology (equal); visualization (equal); writing – original draft (equal). **Jonas Andersson:** Data curation (equal); methodology (equal); visualization (equal). **Tuomas Kaprio:** Investigation (equal); methodology (equal); visualization (equal). **Jaana Hagström:** Investigation (equal); methodology (equal). **Sara Bsharat:** Methodology (equal). **Olof Asplund:** Data curation (equal); methodology (equal). **Gad Hatem:** Data curation (equal); methodology (equal). **Caj Haglund:** Investigation (equal); methodology (equal). **Hanna Seppänen:** Investigation (equal); methodology (equal). **Rashmi B. Prasad:** Data curation (equal); methodology (equal); visualization (equal). **Isabella Artner:** Conceptualization (equal); funding acquisition (equal); investigation (equal); supervision (equal); writing – review and editing (equal).

FUNDING INFORMATION

This research was funded by the Fredrik and Ingrid Thuring
foundation (2020-00596), Swedish Research council (2020-
0146), Novonordisk Foundation (NNF20OC0063485), the
Swedish Cancer foundation (2022), the Royal Physiographic
Society (2022), the Sigrid Juselius Foundation, and the
Cancer and Allergy Foundation (2022, 2023).

CONFLICT OF INTEREST STATEMENT

The authors declare that the research was conducted in
the absence of any commercial or financial relationships
that could be construed as a potential conflict of interest.

DATA AVAILABILITY STATEMENT

The sequencing data generated for this study can be found
in the LUDC repository (www.ludc.lu.se/resources/
repository) under the following accession numbers and
are available upon reasonable request: LUDC2021.10.12
(bulk RNA sequencing data from fetal pancreata), and
LUDC2021.10.18 (single cell RNA sequencing data from
fetal pancreas).

ETHICS STATEMENT

Ethical permission has been obtained from the regional
ethics committee in Lund (Dnr 2012/593, Dnr 2015/241,
Dnr 2018-579). The Surgical Ethics Committee of HUH
approved the study protocol (Dnr HUH 226/E6/06, exten-
sion TKM02 §66 17.4.2013).
ORCID

Ludivine Bertonnier-Brouty https://orcid.org/0000-0003-1294-3244

REFERENCES

1. Kleeff J, Korc M, Apte M, et al. Pancreatic cancer. *Nat Rev Dis Primers*. 2016;2(1):1-22. doi:10.1038/nrdp.2016.22

2. Hu JX, Lin YY, Zhao CF, et al. Pancreatic cancer: a review of epidemiology, trend, and risk factors. *World J Gastroenterol*. 2021;27(27):4298-4321. doi:10.3748/WJG.V27.127.4298

3. Quante AS, Ming C, Rottmann M, et al. Projections of cancer incidence and cancer-related deaths in Germany by 2020 and 2030. *Cancer Med*. 2016;5(9):2649-2656. doi:10.1002/cam4.767

4. Rahib L, Smith BD, Aizenberg R, Rosenzweig AB, Fleshman JM, Matrisian LM. Projecting cancer incidence and deaths to 2030: the unexpected burden of thyroid, liver, and pancreas cancers in the United States. *Cancer Res*. 2014;74(11):2913-2921. doi:10.1158/0008-5472.CAN-14-0155

5. Kamisawa T, Wood LD, Itoi T, Takaori K. Pancreatic cancer. *Lancet*. 2016;388(10039):73-85. doi:10.1016/S0140-6736(16)00141-0

6. Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2021. *CA: Cancer J Clin*. 2021;71(1):7-33. doi:10.3322/caac.21654

7. Naqvi AAT, Hasan GM, Hassan MI. Investigating the role of transcription factors of pancreas development in pancreatic cancer. *Pancreatology*. 2018;18(2):184-190. doi:10.1016/j.pan.2017.12.013

8. Kaneto H, Matsuoka TA, Miyatsuka T, et al. Pdx-1 functions as a master factor in the pancreas. *Front Biosci*. 2008;13(16):6406-6420. doi:10.2741/3162/PDF

9. Zhu Y, Liu Q, Zhou Z, Ikeda Y. Pdx1, Neurogenin-3, and Mafa: critical transcription regulators for beta cell development and regeneration. *Stem Cell Res Ther*. 2017;8(1):1-7. doi:10.1186/s13287-017-0694-z

10. Roy N, Takeuchi KK, Ruggeri JM, et al. Pdx1 dynamically regulates pancreatic ductal adenocarcinoma initiation and maintenance. *Genes Dev*. 2016;30(24):2669-2683. doi:10.1101/gad.291021.116

11. Abel EV, Goto M, Magnuson B, et al. Hnf1a is a novel oncogene that regulates human pancreatic cancer stem cell properties. *eLife*. 2018;7:7. doi:10.7554/eLife.33947

12. Hershko T, Chaussepied M, Oren M, Ginsberg D. Novel link between E2f and P53: proapoptotic cofactors of P53 are transcriptionally upregulated by E2f. *Cell Death Differ*. 2005;12(4):377-383. doi:10.1038/sj.cdd.4401575

13. Ropolo A, Catrinacio C, Renna FJ, et al. A novel E2F1-EP300-VMP1 pathway mediates gemcitabine-induced autophagy in pancreatic cancer cells carrying oncogenic KRAS. *Front Endocrinol*. 2020;11:411. doi:10.3389/fendo.2020.00411

14. Li J, Ran C, Li E, et al. Synergistic function of E2F7 and E2F8 is essential for cell survival and embryonic development. *Dev Cell*. 2008;14(1):62-75. doi:10.1016/j.devcel.2007.10.017

15. Bsharat S, Monni E, Singh T, et al. Mafb-dependent neurotransmitter signaling promotes beta cell migration in the developing pancreas. *Development*. 2023;150(6): dev201009. doi:10.1242/dev.201009

16. GTEx Consortium; Laboratory, Data Analysis & Coordinating Center (LDACC)—Analysis Working Group; Statistical Methods groups—Analysis Working Group. Genetic effects on gene

expression across human tissues. *Nature*. 2017;550(7675):204-213. doi:10.1038/nature24277

17. Gamazon ER, Segre AV, van de Bunt M, et al. Using an atlas of gene regulation across 44 human tissues to inform complex disease- and trait-associated variation. *Nat Genet*. 2018;50(7):956-967. doi:10.1038/s41588-018-0154-4

18. Dobin A, Davis CA, Schlesinger F, et al. STAR: ultrafast universal RNA-Seq aligner. *Bioinformatics*. 2013;29(1):15-21. doi:10.1093/bioinformatics/bts635

19. DeLuca DS, Levin JZ, Sivachenko A, et al. RNA-SeQC: RNA-seq metrics for quality control and process optimization. *Bioinformatics*. 2012;28(11):1530-1532. doi:10.1093/bioinformatics/bts196

20. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. *BMC Bioinformatics*. 2011;12:323. doi:10.1186/1471-2105-12-323

21. Gillespie M, Jassal B, Stephan R, et al. The reactome pathway knowledgebase 2022. *Nucleic Acids Res*. 2022;50(D1):D687-D692. doi:10.1093/nar/gkab1028

22. Lemmetyinen TT, Viitala EW, Wartiovaara L, et al. Fibroblast-derived EGF ligand neuregulin 1 induces fetal-like repro-gramming of the intestinal epithelium without supporting tumorigenic growth. *Dis Model Mech*. 2023;16(4), dmm049692. doi:10.1242/dmm.049692

23. Peng J, Sun B-F, Chen C-Y, et al. Single-cell RNA-seq highlights intra-tumoral heterogeneity and malignant progression in pancreatic ductal adenocarcinoma. *Cell Res*. 2019;29(9):725-738. doi:10.1038/s41422-019-0195-y

24. Oki S, Ohta T, Shioi G, et al. ChIP-Atlas: a data-mining suite powered by full integration of public ChIP-seq data. *EMBO Rep*. 2018;19(12), e46255. doi:10.15252/embr.201846255

25. Tang Z, Kang B, Li C, Chen T, Zhang Z. GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. *Nucleic Acids Res*. 2019;47(W1):W556-W560. doi:10.1093/nar/gkz430

26. Lieber M, Mazzetta J, Nelson-Rees W, Kaplan M, Todaro G. Establishment of a continuous tumor-cell line (panc-1) from a human carcinoma of the exocrine pancreas. *Int J Cancer*. 1975;15(5):741-747. doi:10.1002/ijc.2910150505

27. Guzmán C, Bagga M, Kaur A, Westermarck J, Abankwa D. Colonyarea: an ImageJ plugin to automatically quantify colony formation in clonogenic assays. *PLoS One*. 2014;9(3):e92444. doi:10.1371/journal.pone.0092444

28. Pijuan J, Barceló C, Moreno DF, et al. In vitro cell migration, invasion, and adhesion assays: from cell imaging to data analysis. *Front Cell Dev Biol*. 2019;7:107. doi:10.3389/FCELL.2019.00107/BIBTEX

29. Mao Y, Shen J, Lu Y, et al. RNA sequencing analyses reveal novel differentially expressed genes and pathways in pancreatic cancer. *Oncotarget*. 2017;8(26):42537-42547. doi:10.18632/oncotarget.16451

30. Liu XS, Gao Y, Liu C, et al. Comprehensive analysis of prognostic and immune infiltrates for E2F transcription factors in human pancreatic adenocarcinoma. *Front Oncol*. 2021;10:3297. doi:10.3389/fonc.2020.606735

31. Cancer Genome Atlas Research N, Weinstein JN, Collisson EA, et al. The Cancer Genome Atlas Pan-Cancer analysis project. *Nat Genet*. 2013;45(10):1113-1120. doi:10.1038/ng.2764
32. Uhlen M, Zhang C, Lee S, et al. A pathology atlas of the human cancer transcriptome. *Science*. 2017;357(6352), ean2507.
doi:10.1126/science.aan2507

33. Lonsdale J, Thomas J, Salvatore M, et al. The Genotype-Tissue Expression (GTEx) project. *Nat Genet.* 2013;45:580-585.
doi:10.1038/ng.2653

34. Hu D, Ansari D, Zhou Q, Sasor A, Said Hilmersson K, Andersson R. Stromal fibronectin expression in patients with resected pancreatic ductal adenocarcinoma. *World J Surg Oncol.* 2019;17(1):29. doi:10.1186/s12957-019-1574-z

35. Amrutkar M, Aasrum M, Verbeke CS, Gladhaug IP. Secretion of fibronectin by human pancreatic stellate cells promotes chemoresistance to gemcitabine in pancreatic cancer cells. *BMC Cancer*. 2019;19(1):596. doi:10.1186/s12885-019-5803-1

36. Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2022. *CA Cancer J Clin.* 2022;72(1):7-33. doi:10.3322/CAAC.21708

37. Yu XX, Qiu WL, Yang L, et al. Defining multistep cell fate decision pathways during pancreatic development at single-cell resolution. *EMBO J.* 2019;38(8):e100164. doi:10.15252/embj.2018100164

38. Kim SY, Rane SG. The Cdk4-E2f1 pathway regulates early pancreas development by targeting Pdx1+ progenitors and Ngn3+ endocrine precursors. *Development*. 2011;138(10):1903-1912. doi:10.1242/dev.061481

39. Xu ZH, Qu H, Ren YY, Gong ZZ, Ri HJ, Chen X. An update on the potential roles of E2F family members in colorectal cancer. *Cancer Manag Res.* 2021;13:5509-5521. doi:10.2147/CMAR.S320193

40. Ye L, Guo L, He Z, et al. Upregulation of E2F8 promotes cell proliferation and tumorigenicity in breast cancer by modulating G1/S phase transition. *Oncotarget*. 2016;7(17):23757-23771. doi:10.18632/ONCOTARGET.8121

41. Eoh KJ, Kim HJ, Lee JW, et al. E2F8 induces cell proliferation and invasion through the epithelial–mesenchymal transition and notch signaling pathways in ovarian cancer. *Int J Mol Sci.* 2020;21(16):1-16. doi:10.3390/IJMS21165813

42. Luo L, Zhang G, Wu T, Wu G. Prognostic value of E2F transcription factor expression in pancreatic adenocarcinoma. *Med Sci Monit.* 2021;27:e933443. doi:10.12659/MSM.933443

43. Lan W, Bian B, Xia Y, et al. E2F signature is predictive for the pancreatic adenocarcinoma clinical outcome and sensitivity to E2Ff inhibitors, but not for the response to cytotoxic-based treatments. *Sci Rep.* 2018;8(1):8330. doi:10.1038/s41598-018-26,613-Z

44. Bailey P, Chang DK, Nones K, et al. Genomic analyses identify molecular subtypes of pancreatic cancer. *Nature*. 2016;531(7592):47-52. doi:10.1038/nature16965

45. Oh K, Yoo YJ, Torre-Healy LA, et al. Coordinated single-cell tumor microenvironment dynamics reinforce pancreatic cancer subtype. *Nat Commun.* 2023;14(1):5226. doi:10.1038/s41467-023-40,895-6

46. Ma L, Tian X, Guo H, et al. Long noncoding RNA H19 derived miR-675 regulates cell proliferation by down-regulating E2F-1 in human pancreatic ductal adenocarcinoma. *J Cancer*. 2018;9(2):389-399. doi:10.7150/jca.21347

47. Shao Z, Li C, Wu Q, et al. ZNF655 accelerates progression of pancreatic cancer by promoting the binding of

E2F1 and CDK1. *Oncogene*. 2022;11(1):44. doi:10.1038/s41389-022-00418-2

48. Zhu X, Shi C, Peng Y, et al. Thymidine kinase 1 silencing retards proliferative activity of pancreatic cancer cell via E2F1-TK1-P21 axis. *Cell Prolif.* 2018;51(3):e12428. doi:10.1111/cpr.12428

49. Ren X, Chen C, Luo Y, et al. lncRNA-PLACT1 sustains activation of Nf-kB pathway through a positive feedback loop with $\kappa B\alpha/E2F1$axis in pancreatic cancer. *Mol Cancer*. 2020;19(1):35. doi:10.1186/s12943-020-01153-1

50. Yang Q, Yang B, Chen M. Partner of NOB1 homolog transcriptionally activated by E2F transcription factor 1 promotes the malignant progression and inhibits ferroptosis of pancreatic cancer. *Chin J Physiol.* 2023;66(5):388-399. doi:10.4103/cjop.CJOP-D-23-00063

51. Lammens T, Li J, Leone G, De Veylder L. Atypical E2Fs: new players in the E2F transcription factor family. *Trends Cell Biol.* 2009;19(3):111-118. doi:10.1016/j.tcb.2009.01.002

52. Thurlings I, Martinez-Lopez LM, Westendorp B, et al. Synergistic functions of E2F7 and E2F8 are critical to suppress stress-induced skin cancer. *Oncogene*. 2017;36(6):829-839. doi:10.1038/onc.2016.251

53. Lee S, Park YR, Kim SH, et al. Geraniol suppresses prostate cancer growth through down-regulation of E2F8. *Cancer Med.* 2016;5(10):2899-2908. doi:10.1002/cam4.864

54. Ginsberg D. E2F1 pathways to apoptosis. *FEBS Lett.* 2002;529(1):122-125. doi:10.1016/S0014-5793(02)03270-2

55. Pützer BM, Engelmann D. E2F1 apoptosis counterattacked: evil strikes back. *Trends Mol Med.* 2013;19(2):89-98. doi:10.1016/J.MOLMED.2012.10.009

56. Morgunova E, Yin Y, Jolma A, et al. Structural insights into the DNA-binding specificity of E2F family transcription factors. *Nat Commun.* 2015;6:10050. doi:10.1038/ncomms10050

57. Lopez-Mejia IC, Fajas L. Cell cycle regulation of mitochondrial function. *Curr Opin Cell Biol.* 2015;33:19-25. doi:10.1016/j.ceb.2014.10.006

58. Benevolenskaya EV, Frolov MV. Emerging links between E2F control and mitochondrial function. *Cancer Res.* 2015;75(4):619-623. doi:10.1158/0008-5472.CAN-14-2173

59. Denechaud PD, Fajas L, Giralt A. E2F1, a novel regulator of metabolism. *Front Endocrinol (Lausanne).* 2017;8:311. doi:10.3389/fendo.2017.00311

60. Criscuolo D, Avolio R, Matassa DS, Esposito F. Targeting mitochondrial protein expression as a future approach for cancer therapy. *Front Oncol.* 2021;11:797265. doi:10.3389/fonc.2021.797265

61. Liu N, Funasaka K, Obayashi T, et al. ATAD2 is associated with malignant characteristics of pancreatic cancer cells. *Oncol Lett.* 2019;17(3):3489-3494. doi:10.3892/ol.2019.9960

62. Albahde MAH, Zhang P, Zhang Q, Li G, Wang W. Upregulated expression of TUBA1c predicts poor prognosis and promotes oncogenesis in pancreatic ductal adenocarcinoma via regulating the cell cycle. *Front Oncol.* 2020;10:49. doi:10.3389/fonc.2020.00049

63. Kandimalla R, Tomihara H, Banwait JK, et al. A 15-gene immune, stromal, and proliferation gene signature that significantly associates with poor survival in patients with pancreatic ductal adenocarcinoma. *Clin Cancer Res.* 2020;26(14):3641-3648. doi:10.1158/1078-0432.CCR-19-4044
64. Meng QC, Wang HC, Song ZL, et al. Overexpression of NDC80 is correlated with prognosis of pancreatic cancer and regulates cell proliferation. *Am J Cancer Res.* 2015;5(5):1730-1740.

65. Peng YP, Zhu Y, Yin LD, et al. The expression and prognostic roles of MCMs in pancreatic cancer. *PLoS One*. 2016;11(10):e0164150. doi:10.1371/journal.pone.0164150

66. Smith SJ, Li CM, Lingeman RG, et al. Molecular targeting of cancer-associated PCNA interactions in pancreatic ductal adenocarcinoma using a cell-penetrating peptide. *Mol Ther Oncolytics*. 2020;17:250-256. doi:10.1016/j.omto.2020.03.025

67. Wijnen R, Pecoraro C, Carbone D, et al. Cyclin dependent kinase-1 (CDK-1) inhibition as a novel therapeutic strategy against pancreatic ductal adenocarcinoma (PDAC). *Cancers (Basel)*. 2021;13(17), 4389. doi:10.3390/cancers13174389

## SUPPORTING INFORMATION

Additional supporting information can be found online in the Supporting Information section at the end of this article.

**How to cite this article:** Bertonnier-Brouty L, Andersson J, Kaprio T, et al. E2F transcription factors promote tumorigenicity in pancreatic ductal adenocarcinoma. *Cancer Med*. 2024;13:e7187. doi:10.1002/cam4.7187