JUN mediates the senescence associated secretory phenotype and immune cell recruitment to prevent prostate cancer progression

Background Prostate cancer develops through malignant transformation of the prostate epithelium in a stepwise, mutation-driven process. Although activator protein-1 transcription factors such as JUN have been implicated as potential oncogenic drivers, the molecular programs contributing to prostate cancer progression are not fully understood. Methods We analyzed JUN expression in clinical prostate cancer samples across different stages and investigated its functional role in a Pten-deficient mouse model. We performed histopathological examinations, transcriptomic analyses and explored the senescence-associated secretory phenotype in the tumor microenvironment. Results Elevated JUN levels characterized early-stage prostate cancer and predicted improved survival in human and murine samples. Immune-phenotyping of Pten-deficient prostates revealed high accumulation of tumor-infiltrating leukocytes, particularly innate immune cells, neutrophils and macrophages as well as high levels of STAT3 activation and IL-1β production. Jun depletion in a Pten-deficient background prevented immune cell attraction which was accompanied by significant reduction of active STAT3 and IL-1β and accelerated prostate tumor growth. Comparative transcriptome profiling of prostate epithelial cells revealed a senescence-associated gene signature, upregulation of pro-inflammatory processes involved in immune cell attraction and of chemokines such as IL-1β, TNF-α, CCL3 and CCL8 in Pten-deficient prostates. Strikingly, JUN depletion reversed both the senescence-associated secretory phenotype and senescence-associated immune cell infiltration but had no impact on cell cycle arrest. As a result, JUN depletion in Pten-deficient prostates interfered with the senescence-associated immune clearance and accelerated tumor growth. Conclusions Our results suggest that JUN acts as tumor-suppressor and decelerates the progression of prostate cancer by transcriptional regulation of senescence- and inflammation-associated genes. This study opens avenues for novel treatment strategies that could impede disease progression and improve patient outcomes. Graphical Abstract Supplementary Information The online version contains supplementary material available at 10.1186/s12943-024-02022-x.


Background
Prostate cancer (PCa) is one of the most frequently diagnosed malignancies in men worldwide [1].Its significance lies not only in its prevalence but also in the potential to progress to aggressive forms that resist conventional treatments and lead to high mortality rates [2].The complex molecular programs that determine the routes of PCa progression are still incompletely understood.On the molecular level, the dysregulation of the phosphoinositide 3-kinase (Pl3K) and androgen receptor (AR) pathways has been implicated in the pathology of PCa [3].The constitutive activation of the Pl3K cascade, which is caused by mutations in the tumor-suppressor gene and Pl3K antagonist Phosphate and tensin homologue (PTEN), was identified in 20% of primary PCa tumors and represents a major oncogenic driver [4].The current standard treatment for primary advanced-stage PCa is the administration of anti-androgens to deprive the tumor of dihydrotestosterone.PCa inevitably escapes androgen deprivation by relapsing into castration resistant PCa (CRPC), which is associated with loss of PTEN tumor-suppressor activity in 50% of cases.The characteristic dissemination of CRPC into local and distant regions such as bone, is correlated with poor survival [3][4][5].
In a previously described mouse model, the abrogation of Pten in prostate epithelium (PE) caused activation of a p53-mediated senescence program [6][7][8].The emergence of senescence in cancer is considered a doubleedged sword: it either confers anti-tumorigenic effects when originating from tumor cells or results in protumorigenic outcomes when the tumor microenvironment (TME) is affected [9].This phenomenon is mainly attributed to the induction of a senescence-associated secretory phenotype (SASP), characterized by the secretion of soluble signaling factors, proteases and extracellular matrix proteins [10].In particular, pro-inflammatory cytokines such as IL-6, IL-1, TNF-α, CCL3 and CCL8 attract innate immune cells to the vicinity of the tumor site.As a collective, all components of SASP aid in creation of a pro-tumorigenic microenvironment and ultimately advance tumor progression depending on the tissue context.IL-6 and its downstream effector signal transducer and activator of transcription 3 (STAT3) are known to regulate apoptosis, angiogenesis, proliferation and differentiation, making them promising therapeutic targets in PCa [11].However, our group has recently challenged active IL-6/STAT3 signaling as a tumor driver in PCa, as loss of Stat3 unexpectedly resulted in increased tumor burden and was accompanied by a bypass of PTEN-loss induced cellular senescence (PICS) in a Pten-deficient PCa mouse model [12,13].
Besides the hyperactivation of PI3K/AKT and amplification of AR signaling, other mechanisms driving the progression of PCa include the activation of activator protein-1 (AP-1) mediated gene expression [14].AP-1 transcription factors (TF) such as JUN, were initially considered as proto-oncogenes [15] and deregulation of AP-1 family members was observed in several cancers [16].Previous studies have suggested that JUN modulates hepatocellular tumorigenesis as a regulator of cell cycle genes and has coactivator and repressor functions in the regulation of AR in the prostate [17][18][19].Recent evidence suggests tumorsuppressive functions for several members of the AP-1 TF family and their regulators [17,20].For example, the JUNactivating JUN N-terminal kinase (JNK) has previously been identified as a potent tumor-suppressor in a murine PCa model [21].JUNB, which is also activated by JNK has been associated with growth limiting properties in PCa and its activation may explain the mechanism of JNK's tumorsuppression [22].A recent study provides novel insights how the tumor-suppressive functions of AP-1 might be exerted, as JUN was particularly implicated as pioneering factor in bookmarking the enhancers of genes associated with the induction of the senescence program [23].
Here we investigated the role of Jun in a murine model of Pten-loss driven neoplasia of the PE and surveyed the consequence of JUN-deficiency in tumor development and senescence.

Histological staining
Hematoxylin and eosin (H&E), immunohistochemistry (IHC) and immunofluorescence (IF) stainings were performed on 2 µm sections of formalin-fixed paraffin embedded (FFPE) tissue.H&E staining was done according to routine diagnostic protocols.Details of IHC staining for the different markers are indicated in Supplementary Table 6 and all slides were counterstained with hematoxylin.
For the EpCAM IF staining, slides were dewaxed and heated in pH 6 citrate buffer.After blocking with 2% bovine serum albumin (Roth 8076.4), the slides were incubated in primary antibody (EpCAM, Elab Science, E-AB-70132, dilution 1:300) overnight.Next, slides were incubated for 1 h at room temperature in secondary antibody (Goat anti-Rabbit IgG Alexa Fluor 488, Dilution 1:500) and stained with DAPI.

Human tissue microarray analysis
The generation of human tissue microarrays (TMAs) of healthy and tumor prostate tissues was previously described [28].The TMAs were stained with an antibody for JUN (Supplementary Table 6) and analysed by trained pathologists.The JUN levels were determined by combining the staining intensity with the percentage of positive cells and graded into absent (0), low-grade (1), mediumgrade (2) and high-grade (3).We next stratified TMA samples according to Gleason scores, resulting in three groups (healthy: no Gleason score; low Gleason: Gleason score 5-6; and high Gleason: Gleason score 7-9) and analysed JUN expression for all groups.For the Kaplan-Meier analysis, patients were grouped into absent (0) (JUN absent ) and present (1-3) (JUN present ) JUN expression and correlated with biochemical recurrence (BCR) data.

Whole slide scan analysis
Analysis of IHC staining was performed with QuPath (version 0.3.2) [29].First, regions of interest were annotated, excluding non-prostate tissue such as urethra, seminal vesicles and ductus deferens.Cell detection was performed with the StarDist extension [30] for the NIMP-R14 staining and the built-in watershed cell detection plugin for F4/80, CD79b, JUN, Granzyme B and phosphorylated (p)STAT3.Parameters were chosen individually for each staining.Thereafter, smoothed features were calculated with a FWHM radius of 25 µm.The tissue was then classified into tumor/epithelium and stroma using an object classifier, trained individually for each staining.A threshold was set for the mean DAB optical density value, categorizing cells into positive or negative.For pSTAT3, multiple thresholds were set and cells were classified into 1*, 2* and 3* positive to calculate the H-score.The H-score was calculated by multiplying the percentage of cells by their respective intensity value and ranged from 1 to 300.Analysis was performed by a single investigator and evaluated by two independent pathologists.For quantification of Ki67 levels of tumor and non-tumor samples, we defined four circular regions of interest with a radius of 150 µm.Within each region, we manually counted the positive epithelial cells and used QuPath to detect the negative cells.Percentage of p21 CIP1/WAF1 positive epithelial cells was estimated by a blinded pathologist.For Galactosidase beta 1 (GLB1), regions of interest were annotated and categorized into positive and negative areas using a stringent pixel threshold for the DAB optical density value.The threshold was adjusted to detect the granular expression pattern.Results shown are from the anterior prostate.

Statistical analysis for immunohistochemistry
Measurements were exported as TSV files and imported into GraphPad PRISM (version 9.5.0).Significance was determined using an ordinary one-way ANOVA with Tukey's multiple comparisons tests for 3 or more groups.Graphs were created and formatted in GraphPad PRISM.

Magnetic cell sorting, library preparation and RNA sequencing
The preparation of sequencing libraries and subsequent RNA sequencing (RNA-seq) was performed as previously described [32].Briefly, prostates of 19-week-old mice were dissected, processed to yield a single cell suspension and EpCAM (CD326) positive cells were isolated by magnetic cell sorting (Magnisort ® , Thermo Fisher Scientific) using anti-CD326-biotin (13-5791-82, eBioscience).EpCAM positive cells were collected by centrifugation at 300 xg for 5 min at 4 °C and stored at -80 °C until further use.High-quality RNA, as assessed by 4200 TapeStation System (Agilent) was used for library preparation according to the manufacturer's instructions.

RNA sequencing data analysis
Single-end 75 bp reads sequencing of libraries was performed at CEITEC, Centre for Molecular Medicine (Brno, Czech Republic) as previously described [32].Genes with a false discovery rate (FDR) FDR-adjusted p-value < 0.05 and log2 fold change ≥ 1 or ≤ -1 were considered significantly up-or downregulated.

Kaplan-Meier survival analyses of public datasets
To assess whether expression levels of JUN, PTEN, IL1B, CCL3 and CCL8 affected survival capabilities of human PCa patients, we applied the KM plotter tool (https:// kmplot.com/ analy sis/) which computed probabilities of RFS based on the TCGA-PRAD study [33].Output data were used for re-plotting of survival curves and performing of cox-regression analyses with R packages "Survival" and "Survminer" and R-script "ggsurvplot".Combined KM plotter output was used for calculation of subgroups as stratified by expression levels of both genes, such as JUN and PTEN.Relapse-free survival (RFS) analysis was based on survival data of n = 333 PCa patients.Groups were automatically separated and the calculated, best performing and most significant threshold was used as a cut-off.Hazard ratios and p-values were retrieved from Cox-regression analyses.

JUN levels discriminate progression states in prostate cancer dependent on PTEN
To clarify the role of AP-1 TFs in PCa progression, we investigated the level of the master factor JUN in tissue microarrays (TMA) of low and high progressive human prostate tumors by immunohistochemistry (IHC).We performed semi-quantitative analysis and categorized each tumor based on JUN levels from 0 (absent), 1 (lowgrade), 2 (medium-grade) to 3 (high grade) (Fig. 1a).Patients were divided into present (n = 29 + 6 censored subjects) (JUN present ) and absent (n = 32 + 8 censored subjects) (JUN absent ) cohorts and correlated with biochemical recurrence (BCR) data (Supplementary Fig. 1a).PCa progression is marked by histological changes of the tumor architecture and is categorized by Gleason scoring [34].We observed a gradual decrease of JUN protein abundance from healthy tissue to primary tumors (low Gleason; Gleason score 5-6), reaching the lowest JUN expression state in advanced tumor stages (high Gleason; Gleason score 7-9) (Fig. 1a).The correlation between JUN protein and patient BCR status revealed a significantly (p = 1.8e-02) diminished BCR-free survival in patients with low JUN, whereas high JUN levels were associated with increased survival probability (Supplementary Fig. 1a).We next mined a publicly available transcriptome dataset ( [35]; n = 140) and stratified PCa patients into high-risk and low-risk groups as defined by the prognostic index and characterized by a significant difference in relapse-free survival (RFS) using the Surv-Express webtool [36] (p = 4e-04) (Supplementary Fig. 1b).We investigated JUN mRNA expression in the high-and low-risk groups and found significantly (p = 1.3e-30) higher JUN among low-risk patients compared to the high-risk group (Supplementary Fig. 1c).To explore JUN levels in advanced stages of PCa, we used the Taylor dataset [35], comprising primary tumors of different progression stages and Gleason scores (n = 131) as well as healthy prostate tissue (n = 29).Compared to healthy tissue, we observed higher levels of JUN in early disease stages with Gleason scores 5-6 and significantly decreased expression of JUN in high grade tumors (p = 3e-03; Gleason scores 7-9) (Fig. 1b).Concordantly, JUN was highly expressed in primary tumors (n = 131; n = 65) but significantly lower expressed in PCa metastases (n = 19; n = 25) as observed in two independent datasets (Fig. 1c-d; p = 1.3e-02; [35]; p = 5.3e-09; [37]).We next investigated levels of JUN, JUNB and FOS and observed a comparable regulation (Supplementary Fig. 1d-e).Metastatic CRPC and neuroendocrine PCa (NEPC) present aggressive tumor subtypes that emerge under androgen deprivation ) and low (Gleason score < 7) grade human prostate tumors.Data were retrieved from [35].Significance was determined by an unpaired, two-sided t-test or one-sided Anova.c High and low JUN levels significantly (p = 1.3e-02) discriminate primary prostate tumors (n = 131) (red) and metastases (n = 19) (blue).Data were retrieved from [35].d High and low JUN levels significantly (p = 5.3e-09) discriminate primary prostate tumors (n = 65) (red) and metastases (n = 25) (blue).Data were retrieved from [37].Significances in c-d were determined by an unpaired, two-sided t-test.e Principal component analysis (PCA) of prostate tumors of different developmental stages comprising normal prostate tissue, primary tumors and primary (p) and metastatic (m) CRPC and NEPC tumors.Datasets from [38].f Overlay of JUN expression with PCA clustering from e).JUN levels are color coded from high expression (yellow) to low expression (blue).g Kaplan-Meier survival analysis of TCGA-PRAD [33] tumors (n = 333) assessing levels of JUN and PTEN.Hazard ratios (HR) were determined by Cox-regression analysis: HR(JUN high vs.JUN low ) = 0.461, p = 3.8e-02 and HR(PTEN high vs.PTEN low ) = 0.307, p = 1.5e-03.Statistical testing was done with a logrank test.h Co-analysis between PTEN expression (RNA-Seq by Expectation-Maximization (RSEM) and PTEN protein level reverse-phase protein array (RPPA)).i Co-analysis between PTEN protein level (RRPA) and JUN expression (RSEM) (See figure on next page.)therapy and are associated with poor prognosis.We compared levels of JUN and its related TFs FOS and JUNB in primary (n = 715) and metastatic (n = 320) PCa [38], including CRPC and NEPC (Fig. 1e-f, Supplementary Fig. 1f ).The tumor-subtype and stage-dependent expression of JUN was highly significant when comparing healthy and primary (p = 2.8e-05), primary and metastatic CRPC (p = 2.6e-43) and primary and metastatic NEPC (p = 5.3e-04) (Supplementary Fig. 1g), suggesting JUN as a potential marker of progressive subtypes of PCa.In addition, our survey revealed higher levels of JUN in primary PCa associated with low Gleason scores than healthy prostates (Supplementary Fig. 1g-h), suggesting a gradual change of JUN levels in PCa development and progression.Our data implicate that JUN and other AP-1 factors except MAF and MAFB may act as suppressors rather than drivers of PCa which was reflected by hazard ratios (HR) calculated from RFS (Supplementary Fig. 1i).
Mutations in the tumor suppressor PTEN are considered as main drivers of oncogenic transformation and malignancy in PCa [4].As PTEN loss is highly correlated with increasing Gleason score and associated with activation of several downstream processes, primarily via hyperactivation of PI3K/AKT and inactivation of AR signaling [4], we next investigated synergistic effects of additional JUN alterations.We applied the KMplot tool to assess RFS of PCa patients ( [33]; n = 333) that were stratified into four risk groups JUN high /PTEN high , JUN high /PTEN low , JUN low /PTEN high and JUN low /PTEN low .We observed that patients featuring low levels of JUN and PTEN showed the lowest survival probability whereas patients with high JUN and PTEN expression presented with the most favorable prognosis.In contrast to singular JUN depletion, downregulation of PTEN alone resulted in intermediate survival probabilities confirming its role as main oncogenic driver in PCa (Fig. 1g).Finally, we surveyed reverse-phase protein array (RPPA) data of the TCGA-PRAD cohort [33] and observed that PTEN protein correlated well with PTEN mRNA levels whereas we identified an inverse relationship between PTEN and JUN levels (Fig. 1h-i).Although loss of JUN alone is not sufficient to cause significant changes in survival probability, our data suggest that the absence of PTEN promotes JUN to a survival-determining factor in PCa patients.
To investigate the morphological architecture of prostates upon Jun deletion in the PCa mouse model, we analyzed histological sections by hematoxilin and eosin (H&E) staining (Fig. 2c, top panel).Both wt and Jun PEΔ/Δ animals showed physiological growth patterns and morphology, characteristic for the respective prostate lobes.In Pten PEΔ/Δ and Jun PEΔ/Δ ; Pten PEΔ/Δ prostates, we observed hyperplastic epithelium growing in cribriform patterns into the lumen.Both groups showed anisocytosis, anisokaryosis and alterations in nucleus-to-cytoplasmic ratios, but largely without invasion of the stroma.
Next, we analyzed JUN levels in prostates of all genotypes.Supporting our immunoblot results, IHC revealed increased levels of total JUN predominantly in the PE of Pten PEΔ/Δ mice and absence in epithelial cells of Jun PEΔ/Δ and Jun PEΔ/Δ ; Pten PEΔ/Δ (Fig. 2c, bottom panel).We assessed the effects of Jun deficiency on tumor burden and survival by morphological and survival analyses.Macroscopically, prostates from Pten PEΔ/Δ and Jun PEΔ/Δ ; Pten PEΔ/Δ mice were notably enlarged as compared to wt or Jun PEΔ/Δ prostates (Fig. 2d).This finding was corroborated by prostate weight analysis (Fig. 2e).The additional deletion of Jun on the Pten-deficient background resulted in even higher prostate weights, hinting at JUN's potential function as a tumor-suppressor in murine PCa development.We performed a Kaplan-Meier survival analysis where overall survival or the occurrence of the discontinuation criteria according to the guidelines of the 3Rs principles were defined as the endpoint of the experiments (Fig. 2f ) [40].We observed comparable survival probabilities of wt and Jun PEΔ/Δ mice (p = 5.4e-01) but a significantly decreased survival of Pten PEΔ/Δ (mean survival 85.3 weeks, p = 6e-03) as compared to wt mice.Remarkably, the survival of Pten PEΔ/Δ mice was significantly (p < 1e-04) reduced by the additional deletion of Jun.Jun PEΔ/Δ ; Pten PEΔ/Δ mice exhibited a mean survival of 67.2 weeks.Despite the significantly reduced survival rates in Jun PEΔ/Δ ; Pten PEΔ/Δ mice, we did not detect metastatic lesions in the analysed genotypes (Fig. 2g).We therefore conclude that Jun-deficiency alone is not sufficient to induce prostate tumorigenesis, but causes a significant increase in tumor burden and a significant reduction in overall survival in combination with Pten knockout.The results of our murine PCa model reinforce our observations from human PCa samples, suggesting that JUN acts as a tumor-suppressor in PCa.
To determine whether aberrant cellular proliferation contributes to enhanced tumor growth in Jun PEΔ/Δ ; Pten PEΔ/Δ -deficient prostates, we assessed the number of Ki67 + epithelial cells by IHC.Although we noticed higher Ki67 levels in Jun PEΔ/Δ ; Pten PEΔ/Δ tumors by trend, the difference was not significant (p = 1.3e-01) when compared to Pten PEΔ/Δ prostates (Supplementary Fig. 2a).To investigate the effects of JUN ablation in vitro, we utilized the CRISPR/Cas9 technology in the human PCa cell lines DU145 (PTEN wildtype) and PC3 (PTEN mutated).We designed three individual guide RNAs for the JUN locus (Supplementary Fig. 2b) and used lentiviral transduction of empty vector (EV) and guide RNA (G1, G12, G14) plasmids.We identified varying efficiencies of JUN knockout in bulk cultures of DU145 and PC3 cell lines (Supplementary Fig. 2c, e) and no significant differences in cellular proliferation (Supplementary Fig. 2d, f ).We confirmed the results of unchanged proliferation in single clones of both cell lines which were selected according to complete loss of JUN protein (Supplementary Fig. 2g-j).The in vivo and in vitro results indicate that proliferation may not be the primary biological process influenced by JUN during PCa progression.

Transcriptome profiling reveals JUN-mediated alterations in senescence-associated secretion and immune response
To elucidate the tumor cell-specific molecular programs regulated by JUN in vivo, we performed transcriptome profiling of PE cells across all four experimental murine groups (Fig. 2a).To obtain a homogenous epithelial fraction, we enriched prostate lysates for the Epithelial cell adhesion molecule (EpCAM) showing a uniform expression in PE cells (Fig. 3a, Supplementary Fig. 3a) via magnetic cell separation [32] (Fig. 3b, Supplementary Fig. 3b).The correlation analysis revealed high congruence between Jun PEΔ/Δ ; Pten PEΔ/Δ and Pten PEΔ/Δ tumor and wt and Jun PEΔ/Δ samples (Fig. 3c).
We therefore investigated the enrichment level of different senescence signatures including "oncogene-induced senescence" (OIS), "SASP" signatures and the novel "Sen-Mayo" gene signature, consisting of 125 previously identified senescence/SASP-associated factors.SenMayo genes are transcriptionally regulated by senescence and allow identification of senescent cells across tissues [41].Sen-Mayo genes were significantly (qval = 2.40e-02) enriched in Pten PEΔ/Δ prostates and depleted (qval = 2.64e-02) in Jun PEΔ/Δ ; Pten PEΔ/Δ tumors (Fig. 3d).Among the depleted SenMayo genes in Jun-deficient Pten PEΔ/Δ prostates, we identified chemokines such as Ccl3, Ccl4 and Ccl8, along with pro-inflammatory cytokines such as Il1b and Tnfa (Fig. 3e).As these secreted cytokines and chemokines represent well described SASP factors, we next investigated a SASP core gene signature previously described in a Pten-deficient prostate model [42].Using GSEA, we indeed detected enrichment of the SASP core signature in Pten-deficient prostates which was reverted in Jun PEΔ/Δ ; Pten PEΔ/Δ animals (Fig. 3f ).To investigate further aspects of JUN-dependent regulation of senescence in Pten-deficient murine prostates, we stained formalin-fixed paraffin embedded (FFPE) material with the senescence markers p16 INK4A , p21 CIP1/WAF1 and Galactosidase beta 1 (GLB1) (Supplementary Fig. 3c).We did not observe differences in the amount of p16 INK4A positive cells between Pten PEΔ/Δ and Jun PEΔ/Δ ; Pten PEΔ/Δ tumors, but found significant changes in staining patterns.While we detected prominent nuclear staining in Pten PEΔ/Δ samples, Jun PEΔ/Δ ; Pten PEΔ/Δ revealed predominantly cytoplasmic localization, hinting at a potential inactivation of p16 INK4A via nuclear export [43].In wt and Jun PEΔ/Δ prostates, we observed a weak lobe-dependent expression pattern of p21 CIP1/WAF1 .Conversely, in Pten PEΔ/Δ and Jun PEΔ/Δ ; Pten PEΔ/Δ samples, p21 CIP1/WAF1 was expressed in each individual epithelial cell, with no discernible difference between the two groups.GLB1 staining displayed its characteristic granular expression pattern prompting us to quantify percentage of positive area however we found no significant difference between Pten PEΔ/Δ and Jun PEΔ/Δ ; Pten PEΔ/Δ groups.Apart from changes in the p16 staining pattern, we found no significant deregulation of the classic senescence-associated cell cycle markers, implicating that JUN affects the SASP but not senescence-associated cell cycle arrest.
As our results suggest JUN-dependent activation of the IL-6/STAT3 axis and our previous study connected loss of activated STAT3 in Pten-deficient PCa to increased tumor burden via disruption of senescence [12], we sought to analyze STAT3 tyrosine 705 (Y705) phosphorylation (pSTAT3 Y705 ) in the Jun-deficient background.We indeed detected reduced levels of pSTAT3 Y705 in both stroma (p = 5.0e-04) and epithelial cells (p < 1.0e-04) of Jun PEΔ/Δ ; Pten PEΔ/Δ compared to Pten PEΔ/Δ tumors (Supplementary Fig. 3d upper panel, Supplementary Fig. 3e) while total STAT3 levels remained constant (Supplementary Fig. 3d, lower panel).Our findings provide evidence that loss of JUN accompanied by reduced activation of STAT3 bypasses SASP and subsequently amplifies the tumor load in Jun PEΔ/Δ ; Pten PEΔ/Δ animals.We suggest an interplay of JUN and STAT3 mediating senescence-associated secretion of inflammatory factors in PCa in vivo, reinforcing JUN's proposed function as a pioneering factor of senescence [23].
Cells of the innate immune system, including neutrophil granulocytes, mast cells and macrophages serve as the primary defense against infections and consequently recruit T and B cells to infection sites [45].Among the DEGs of Pten PEΔ/Δ versus Jun PEΔ/Δ ; Pten PEΔ/Δ prostates, we identified neutrophil movement-specific gene signatures that play a crucial role in the recruitment of immune cells (Fig. 4b) [46].We observed that cytokines involved in chemotaxis of immune cells such as Ccl3, Ccl8 and Il1b were significantly deregulated between groups (Fig. 4c).To further dissect the potentially involved immune cell subsets, we conducted single sample GSEA using the M5 ontology gene sets signature from the molecular signature database (MsigDB).We identified enrichment of macrophage-and neutrophilspecific gene signatures characterized by cellular activities such as migration, activation/differentiation and enhanced expression indicating production of MIP1α/ CCL3 and GM-CSF.Moreover, single sample GSEA revealed processes related to other immune cell subsets such as mast cells, myeloid cells and CD8 + T cells that were significantly enriched in Pten PEΔ/Δ compared to wt prostates and depleted in Jun PEΔ/Δ ; Pten PEΔ/Δ (Fig. 4d).This implicates JUN in the control of inflammatory states during PCa progression.We validated the JUN-dependent regulation of IL-1β, TNF-α and NLRP3, all involved in the regulation of inflammatory response processes by immunoblot and cytokine analyses (Fig. 4e-f and Supplementary Fig. 4b).
To further examine the apparent shifts in immune system-related transcriptomic signatures, we assessed granulocytic or lymphocytic cell infiltrations based on microscopic characteristics in H&E staining of all four genotypes (Supplementary Fig. 4c).We detected no or low-grade infiltration by inflammatory cells in wt and Jun PEΔ/Δ specimens.In contrast, Pten PEΔ/Δ mouse prostates exhibited increased levels of high-and middlegrade infiltrations, which were significantly mitigated in Jun PEΔ/Δ ; Pten PEΔ/Δ prostates.Increased immune cell infiltration of Pten PEΔ/Δ prostates as identified by histopathological analysis therefore supported the results of transcriptome profiling.This highlights the importance of JUN in the regulation of inflammation by affecting the secretion of pro-inflammatory cytokines in Pten-deficient PCa.

Epithelial JUN deficiency modulates the migration of innate immune cells from the periphery
To investigate the distribution and abundance of infiltrating immune cells, we performed IHC stainings.Neutrophils and inflammatory monocytes were stained using the antibody clone NIMP-R14, which targets the specific cell surface markers and differentiation antigens Ly-6G and Ly-6C (Fig. 5a).In Pten PEΔ/Δ (See figure on next page.)Fig. 4 JUN expression determines the level of immune cell infiltration of Pten-loss driven tumors.a Heat map showing JUN-dependent regulation of genes related to innate immunity (upper panel) and inflammatory response (lower panel) in wt, Jun PEΔ/Δ , Pten PEΔ/Δ and Jun PEΔ/ Δ ;Pten PEΔ/Δ prostates.JUN-dependent core factors such as Il1b, Nlrp3 and Ccl5 are highlighted.b Heat map presenting the JUN-dependent regulation of genes involved in migration and chemotaxis of neutrophil granulocytes in Pten PEΔ/Δ and Jun PEΔ/Δ ;Pten PEΔ/Δ prostates.Genotypes and expression levels in a-b are color coded.c Expression levels (log2, FPKM) of Ccl3, Ccl8 and Il1b are significantly (Ccl3, p = 2.4e-04; Ccl8, p = 9.7e-05 and Il1b, p = 5.0e-03) reduced in EpCAM + cells of Jun PEΔ/Δ ;Pten PEΔ/Δ prostates.Significance was determined by an unpaired two-sided t-test.d Single-sample GSEA analysis using the M5 signature of Broad Institute's molecular signature database (MsigDB) revealing enrichment of macrophage-and neutrophil-associated properties in Pten PEΔ/Δ compared to Jun PEΔ/Δ ;Pten PEΔ/Δ prostates.e Western blot analysis of NLRP3 and non-cleaved Pro-IL-1β in all four experimental groups in biological replicates.β-ACTIN served as loading control.f Multiplex immunoassay of homogenized prostate samples of 19-week-old wt, Jun PEΔ/Δ , Pten PEΔ/Δ and Jun PEΔ/Δ ;Pten PEΔ/Δ animals for analysis of IL-1β levels in pico grams (pg)/ ml of indicated biological replicates.Statistical testing was done with one-way Anova, significant p-values are indicated prostates, we observed high numbers of neutrophils migrating from the blood vessels across the stroma into the epithelium, where they predominantly accumulated, and subsequently advanced into the lumen.In Jun PEΔ/Δ ; Pten PEΔ/Δ prostates, we detected significantly (p < 1.0e-04) less neutrophils in the stroma and epithelium, but the migration patterns remained consistent with Pten PEΔ/Δ tumors (Fig. 5b).In contrast, macrophages, stained by the marker F4/80 were primarily located in the stroma, with no significant differences between the groups (Fig. 5c-d).We observed significantly (p = 4.0e-04) less macrophages infiltrating the epithelium in prostates with additional deficiency of Jun.In conclusion, Pten PEΔ/Δ displayed a highly immune infiltrated phenotype, which was substantially reverted in prostates with additional deficiency of Jun.This observation suggests that JUN may be essential for tumor cell recognition by innate and consequently adaptive immune cells.
Neutrophils attract T cells to the site of inflammation via secretion of chemokines such as CCL2 and CCL5 [47,48].We utilized multiplex IHC to discern the T cell subsets, employing a marker panel consisting of CK/CD3/CD4/CD8/CD45/PD-1/DAPI.We observed various T cell subpopulations (T helper (CD4 + ) and cytotoxic T cells (CD8 + ), PD-1 positive and negative) mainly in the stroma and to a lesser degree in the epithelium (Supplementary Fig. 5a-b).To further investigate active cytotoxic T cells and natural killer cells, we performed Granzyme B IHC (Supplementary Fig. 5c-d).We did not observe a significant effect of Jun deficiency on any of the investigated populations.Additionally, we investigated the infiltration of B cells, stained by CD79b.B cells were found almost exclusively in the stroma, with significantly less infiltration in Jun PEΔ/Δ ; Pten PEΔ/Δ compared to Pten PEΔ/Δ prostates (Fig. 5e-f ).In summary, IHC validated the JUN-dependent modulation of the immune cell compartment, particularly affecting innate immune cells.This phenotype was likely provoked by a JUN-dependent regulation of neutrophil attracting chemokines such as IL-1β.

Increased expression of SASP factors is correlated with prolonged survival in prostate cancer
To translate our findings to the human disease, we investigated a potential association of JUN and SASP factors.We compared the levels of JUN, IL1B, CCL3 and CCL8 in patient data (TCGA-PRAD [33]) by PCA and found high IL1B, CCL3 and CCL8 mRNA in tumors expressing high levels of JUN.In contrast, tumors expressing low JUN levels revealed equally low amounts of IL1B, CCL3 and CCL8 (Fig. 6a-d).Upon separation of TCGA-PRAD tumors by using mean JUN expression as cut-off in JUN high and low expressing groups, we indeed confirmed enrichment of SASP factors and neutrophil marker ADGRE1 in JUN high subgroups (Fig. 6e).We next performed a Pearson correlation analysis and detected a significant but weak positive (R ≤ 0.64; p ≤ 0.05) association of IL1B, CCL3 and CCL8 with JUN and additional AP-1 factors such as JUND, JUNB, FOS and FOSB (Fig. 6f ).
Assuming that JUN mediates tumor-suppressor activity via positive regulation of SASP factors required for the recruitment of immune cells, we expected that high levels of SASP factors may be associated with favorable prognosis.Hence, we asked whether cytokine expression may act in concert with JUN to influence patient survival.We stratified patients according to their PTEN, JUN and IL1B (Fig. 6g) CCL3 (Fig. 6h) or CCL8 (Fig. 6i) levels and compared the RFS between groups.As expected, PTEN high groups (presented in gray shades) showed overall favorable outcomes and did not significantly differ from one another.In groups where all marker genes were lowly expressed we generally detected the worst prognosis while groups with a singular lowly expressed marker showed intermediate prognosis.CCL8 did not Fig. 5 Histological analysis of infiltrating immune cells reveals downregulated innate immune response in Jun PEΔ/Δ ;Pten PEΔ/Δ prostates.a Representative images of IHC stainings of NIMP-R14, a pan-marker of neutrophil granulocytes, indicating high neutrophil infiltration of Pten PEΔ/Δ prostates, reverted by the additional loss of Jun in Jun PEΔ/Δ ;Pten PEΔ/Δ prostates.Top row: 20.0 × magnification, scale bar represents 150 µm; Bottom row: 63.0 × magnification, scale bar represents 40 µm.b Quantification of NIMP-R14 + neutrophils in epithelium (left) and stroma (right).A significantly decreased (p < 1e-04) infiltration of neutrophils in tumors and adjacent stroma of Jun PEΔ/Δ ;Pten PEΔ/Δ prostates is evident.c Representative images of IHC stainings for the pan-marker of macrophages F4/80.A high infiltration of Pten PEΔ/Δ prostates and adjacent stroma by macrophages is evident and reverted by the additional loss of Jun in Jun PEΔ/Δ ;Pten PEΔ/Δ prostates.Top row: 40.0 × magnification, scale bar represents 60 µm; Bottom row: 100.0 × magnification, scale bar represents 30 µm.d Quantification of F4/80 + macrophages in epithelium (left) and stroma (right).A significantly decreased (p = 4e-04) infiltration of macrophages in tumors but not adjacent stroma (p = 8.3e-01) of Jun PEΔ/ Δ ;Pten PEΔ/Δ prostates is evident.e Representative images of IHC stainings of B cell infiltration using the pan-marker CD79b.A high infiltration of stroma adjacent to Pten PEΔ/Δ prostates by CD79b + B cells is evident and reverted by the additional loss of Jun in Jun PEΔ/Δ ;Pten PEΔ/Δ prostates.Top row: 40.0 × magnification, scale bar represents 60 µm; Bottom row: 100.0 × magnification, scale bar represents 30 µm.f Quantification of B cells in epithelium (left) and stroma (right).B cell infiltration as observed in the stroma of Pten PEΔ/Δ prostates was significantly decreased (p < 1e-04) in Jun PEΔ/Δ ;Pten PEΔ/Δ prostates.Statistical significance between Pten PEΔ/Δ and Jun PEΔ/Δ ;Pten PEΔ/Δ groups are indicated in b, d and f (See figure on next page.)follow the observed IL1B and CCL3 profiles suggesting that this chemokine does not act in concert with JUN in the absence of PTEN.Finally, we investigated a potential relationship between JUN and STAT3 activation and explored STAT3's role in immune modulation in TCGA-PRAD data [33].We correlated RPPA of pSTAT3 Y705 with levels of JUN and IL1B.We observed a weak but significant correlation of JUN (R = 0.47, p < 2.2e-16) and IL1B (R = 0.48, p = 2.2e-16) with pSTAT3 Y705 (Fig. 6j).As observed for JUN, PCa exhibiting a high (> 7) Gleason score showed reduced levels of pSTAT3 Y705 (p = 1.6e-02) when compared to low risk tumors (Gleason score ≤ 7) (Fig. 6k).To correlate the main oncogenic driver PTEN to pSTAT3 Y705 , we grouped patients according to varying PTEN levels and observed a dose dependent decrease of pSTAT3 Y705 (Fig. 6l).These results were analogous to our previous findings where we detected gradually reduced JUN expression in PTEN medium and high expressing tumors (Fig. 1i).These data might hint at interconnected mechanisms of both transcriptional regulators.In summary, we propose that levels of JUN and STAT3 potentially orchestrated via PTEN, determine progression stages of human prostate tumors by modulating the immune response through regulation of cytokines and interleukins as identified in the Jun-deficient murine PCa model.

Discussion
PCa is among the most frequently diagnosed malignancies in men worldwide and a significant number of patients progress to advanced and lethal stages.The mortality linked to metastatic PCa highlights the pressing need to elucidate its intricate mechanisms and pinpoint viable therapeutic interventions.Despite this urgency, the cellular mechanisms and environmental contexts that control PCa development and progression remain incompletely understood.Loss of PTEN is evident in 20% of primary human prostate carcinomas and escalates in 50% of metastatic CRPC [4].Comparable to the human situation, Pten loss leads to the formation of precancerous lesions in PE cells in mouse models [49,50].Aggressive carcinomas develop only in the presence of additional mutations [51], such as abnormal expression of ERG [52], loss of IL-6/STAT3 functionality [12,13], dysfunction of the methyltransferase Kmt2c [32] or activation of the RAS/MAPK cascade [51,53].While several studies indicate that augmented JUN expression drives PCa progression [14,54], the functional role of JUN and AP-1 TFs in PCa remains controversial.Intriguingly, genetic disruption of JunB in vivo accelerated the progressive phenotype of Pten-deficient PCa [22].A recent study using in vivo CRISPR to achieve combinatorial deletion of Pten and Fos in the PE led to increased tumorigenesis potentially via upregulation of Jun [55].In contrast, we did not observe compensatory functions of other AP-1 family members in our genetically engineered mouse model or the analysis of human patient data.We therefore speculate that Fos deletion might be more susceptible to influencing levels of other AP-1 members than Jun or that the upregulation of JUN might arise due to differential CRISPR targeting efficiencies of Pten and Fos in the chosen model.An alternative explanation for the increased tumorigenesis observed upon Fos knockout in PCa is provided by evidence that FOS stimulates the trans-activation properties of JUN but represses its AR co-activator function [56].Loss of Fos might therefore predominantly favor Jun's pro-proliferative co-activator function and weaken its anti-proliferative trans-activator function [57].Our results similarly point towards a context dependent tumor-suppressive role, rather than a driving function of JUN in PCa progression.
In the present study, our focus was to delineate the role of JUN in PCa.We first examined JUN levels in clinical PCa samples and analyzed JUN patterns across varying progression stages from three publicly available datasets [35,37,38].We found that JUN expression increased in tumors normal prostates, however did not observe the same effect when we compared healthy to low stages in TMA data.It remains to be elucidated whether this discrepancy can be explained by differences in protein versus RNA levels or the PTEN mutation status in early stages of malignant transformation.Importantly, the levels of JUN, FOS and JUNB and the levels of JUN in the TMA dataset were all significantly decreased with progression of PCa.This suggests that high JUN levels may protect from development of progressive disease, a hypothesis further supported by the increased survival rates of patients harboring high JUN expressing tumors.By co-integrating the tumor suppressor PTEN in our survival analysis, we found that levels of PTEN and JUN determine survival probabilities of PCa and revealed the worst prognosis in JUN low /PTENlow PCa but highest RFS in JUN high /PTEN high tumors.Encouraged by these findings, we studied the functional role of JUN in a murine PCa model, characterized by homozygous loss of Pten (Pten PEΔ/Δ ) [26,39].Mirroring  6 is TCGA-PRAD [33] the early prostatic intraepithelial neoplasia (PIN) stages of human PCa, JUN was significantly upregulated in Pten PEΔ/Δ prostates.Consistent with human patient data, depletion of Jun alone had no effect on the morphological architecture and growth of the prostate.Epithelial cells of Pten PEΔ/Δ prostates developed hyperplasia, subsequently forming prostate adenocarcinoma and rapidly progressing upon additional deletion of Jun.The aggressive phenotype observed in Jun PEΔ/Δ ; Pten PEΔ/Δ prostates resulted in decreased survival of mice and increased prostate weight and size.We did not detect signs of severe organ dysfunction, systemic inflammation or metastatic disease (Fig. 2g).
The TME is a dynamic system characterized by chronic inflammation and participation of diverse host components, but plays a pivotal role in cancer progression [58].Within the TME, immune cells such as tumor-associated macrophages (TAM) and tumor-associated neutrophils (TAN) both foster cancer progression or combat tumor cells, underscoring their dual roles in tumorigenesis [59][60][61].Central to this environment is the SASP, where senescent cells release a plethora of inflammatory mediators.SASP-driven effects often culminate in the immune-mediated clearance of potential tumorigenic cells, a process termed "senescence surveillance" [62,63].Our histopathologic examination of Pten PEΔ/Δ PCa samples revealed significant enrichment of neutrophils and macrophages that infiltrated the tumors and adjacent stroma.Concurrent deletion of Jun strikingly reduced tumor infiltration with neutrophils and macrophages and accelerated tumor growth.Transcriptomic analyses of Pten PEΔ/Δ and Jun PEΔ/Δ ; Pten PEΔ/Δ prostates revealed a JUN-dependent modulation of SASP-associated genes, but we did not identify compensatory upregulation of other AP-1 members as it has been described upon inactivation of FOS [64].To investigate further aspects of senescence and address senescence-associated cell cycle arrest, we conducted IHC stainings for p16 INK4a , p21 CIP/ WAF1 and GLB1.We did not observe quantitative differences in expression of these classical senescence markers indicating that JUN is involved in the regulation of SASP, but not senescence-associated cell cycle arrest.The regulation of SASP without affecting cell cycle arrest has been previously demonstrated.This was evidenced by the association of the chromatin reader BRD4 with recruitment to enhancer regions activating SASP genes in senescent cells [65].Recent findings have implicated AP-1 and in particular JUN as pioneering factors on a specific enhancer landscape essential for the execution of senescence-controlling programs [23].In line with these results, we propose that loss of Pten coupled with an increase in JUN likely instigates a JUN-driven SASP phenotype.
SASP involves the expression and secretion of inflammatory cytokines such as CCL3, CCL8, IL-1β and TNF-α [32,66] which subsequently recruit immune cells such as neutrophils, macrophages and T cells [10,67,68].As we observe downregulated secretion of IL-1β and TNF-α in Jun PEΔ/Δ ; Pten PEΔ/Δ prostates, we suggest that Jun depletion in Pten-deficient prostates disrupts SASP.This impedes the recruitment of neutrophils and macrophages, as well as tumor cell clearance by macrophages and dendritic cells [9,62].We thus propose JUN as a key regulator of SASP.Our results support a previous study where JUN depletion was linked to diminished inflammatory responses and reverting the senescent/SASP phenotype of RAS-OIS fibroblasts to a proliferating phenotype [23].Furthermore, GM-CSF, a direct JUN target has been shown to amplify macrophage and neutrophil immune responses [69] and modulate pro-inflammatory cytokine secretion such as TNFα and IL-6 [70].
Another intriguing mechanism of JUN-dependent modulation of the immune phenotype in PCa may depend on STAT3 levels.Our previous work identified activation of STAT3 and a p19 ARF -MDM2-p53 axis to induce senescence upon Pten depletion [12].Consistently, Jun loss was associated with decreased IL-6-JAK-STAT3 signaling, evidenced by significantly reduced pSTAT Y705 levels in Jun PEΔ/Δ ; Pten PEΔ/Δ prostates.ENCODE database exploration [71] revealed mutual promoter binding sites for JUN and STAT3 suggesting a potential JUN-STAT3 interplay in impacting senescence pathways in PCa (https:// maaya nlab.cloud/ Harmo nizome/ datas et/ CHEA+ Trans cript ion+ Factor+ Targe ts).This interplay is supported by results of a STAT3 binding analysis in CD4 + T cells, which suggests that STAT3 directly regulates the expression of Jun and Fos and may potentially function in a positive feedback loop [72].Therefore, therapeutic activation of STAT3 potentially causes SASP factor modulation and may elevate JUN levels in tumors, thereby restricting tumor progression and enhancing PCa patient survival.

Conclusions
In summary, our data suggest that JUN functions as a pivotal regulator of SASP and survival in PTEN-deficient PCa, orchestrating the recruitment dynamics of TAMs and TANs within the TME.Given the indispensable role of robust SASP in immune surveillance of preneoplastic anomalies, its therapeutic modulation presents intricate challenges.Our recent investigations have shown the potential of the antidiabetic agent metformin, which curtails multiple pro-inflammatory SASP components by inhibiting NF-κB nuclear translocation [73].Metformin increases STAT3 in advanced PCa cases, leading to significant tumor growth attenuation, underscored by reduced mTORC1/CREB and AR levels in a PCa murine model [13].The interplay between JUN, STAT3 and PTEN might represent a key mechanism that could be exploited for therapeutic advances.

Fig. 1
Fig. 1 JUN levels are correlated with prostate cancer progression stages.a Left panel: Representative immunohistochemistry (IHC) images of tissue microarrays (TMAs) investigating human prostate tumors (n = 60) with high or low Gleason scores stained for JUN protein.Scale bars indicate 150 µm (top row) and 30 µm (bottom row), images are presented in 16.8 × (top row) and 80.0 × magnification (bottom row).The area used for the higher magnification is indicated by the rectangle.Right panel: Violin plot showing JUN expression divided in absent (0), low-grade (1), medium-grade (2) and high-grade (3) in healthy (no Gleason score), low Gleason (Gleason score 5-6) and high Gleason (Gleason score 7-9) TMA samples.b JUN mRNA levels in high (Gleason score ≥ 7) and low (Gleason score < 7) grade human prostate tumors.Data were retrieved from[35].Significance was determined by an unpaired, two-sided t-test or one-sided Anova.c High and low JUN levels significantly (p = 1.3e-02) discriminate primary prostate tumors (n = 131) (red) and metastases (n = 19) (blue).Data were retrieved from[35].d High and low JUN levels significantly (p = 5.3e-09) discriminate primary prostate tumors (n = 65) (red) and metastases (n = 25) (blue).Data were retrieved from[37].Significances in c-d were determined by an unpaired, two-sided t-test.e Principal component analysis (PCA) of prostate tumors of different developmental stages comprising normal prostate tissue, primary tumors and primary (p) and metastatic (m) CRPC and NEPC tumors.Datasets from[38].f Overlay of JUN expression with PCA clustering from e).JUN levels are color coded from high expression (yellow) to low expression (blue).g Kaplan-Meier survival analysis of TCGA-PRAD[33] tumors (n = 333) assessing levels of JUN and PTEN.Hazard ratios (HR) were determined by Cox-regression analysis: HR(JUN high vs.JUN low ) = 0.461, p = 3.8e-02 and HR(PTEN high vs.PTEN low ) = 0.307, p = 1.5e-03.Statistical testing was done with a logrank test.h Co-analysis between PTEN expression (RNA-Seq by Expectation-Maximization (RSEM) and PTEN protein level reverse-phase protein array (RPPA)).i Co-analysis between PTEN protein level (RRPA) and JUN expression (RSEM)

Fig. 3
Fig.3 Transcriptome profiling of genetic models reveals a JUN-dependent regulation of innate immunity.a Representative immunofluorescence (IF) image of a wt murine prostate for the epithelial marker EpCAM (green).DAPI (blue) is shown as a nuclear stain.Top image: 40.0 × magnification, scale bar represents 60 µm; Bottom image: 147.5 × magnification, scale bar represents 20 µm.b Overview of sample preparation for transcriptome profiling of wt, Pten PEΔ/Δ , Jun PEΔ/Δ and Jun PEΔ/Δ ;Pten PEΔ/Δ prostate samples of 19-week-old animals.An antibody against the epithelial marker EpCAM was used to separate single cell suspensions of minced and digested prostates into EpCAM positive (pos) and negative (neg) fractions by magnetic cell sorting.EpCAM pos cells were used for RNA-seq expression profiling.c Heat map showing correlation analysis of tumor samples described in b) regarding global similarity of samples.The Pearson correlation coefficient (R) is shown (color coded).d Gene onthology (GO)-enrichment analysis of differentially expressed genes (DEGs) showing the top differentially regulated pathways between Pten PEΔ/Δ and Jun PEΔ/Δ ;Pten PEΔ/Δ .Significance as shown by FDR is color coded, enriched (positive normalized enrichment score (NES)) or depleted (negative NES) processes are indicated.Asterisk represents non-significant pathways (ns).e Heat map showing SenMayo genes most significantly (p ≤ 1e-02) regulated among Pten PEΔ/Δ and Jun PEΔ/Δ ; Pten PEΔ/Δ prostates.f GSEA enrichment analysis using the Guccini_core_SASP gene set in Pten PEΔ/Δ versus Jun PEΔ/Δ ;Pten PEΔ/Δ and Pten PEΔ/Δ versus wt animals.g Heat map representation of wt, Pten PEΔ/Δ , Jun PEΔ/Δ and Jun PEΔ/Δ ;Pten PEΔ/Δ samples showing DEGs."Innate immunity", FDR = 7.64e-05; "Immune system", FDR = 2.77e-04 and "Extracellular space", FDR = 6.60e-03 related processes most discriminated the groups.Genotypes and expression levels are color coded.h GO-enrichment analysis of DEGs showing the regulation of innate immune cells such as neutrophil granulocytes.Significance as shown by p-value is color coded, enriched (positive NES) or depleted (negative NES) processes are indicated.Shown are the signaling pathways enriched in Pten PEΔ/Δ tumors compared to wt (left side) and Jun PEΔ/Δ ;Pten PEΔ/Δ tumors compared to Pten PEΔ/Δ (right side)