Next Article in Journal
Molecular Pincers Using a Combination of N-H and C-H Donors for Anion Binding
Previous Article in Journal
Scion-to-Rootstock Mobile Transcription Factor CmHY5 Positively Modulates the Nitrate Uptake Capacity of Melon Scion Grafted on Squash Rootstock
Previous Article in Special Issue
A Novel L-Phenylalanine Dipeptide Inhibits the Growth and Metastasis of Prostate Cancer Cells via Targeting DUSP1 and TNFSF9
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 1
10.3390/ijms24010164
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Age- and Stage-Dependent Prostate Cancer Aggressiveness Associated with Differential Notch Signaling

by
Magdalena Julita Orzechowska
*,†,
Dorota Anusewicz
† and
Andrzej K. Bednarek
Department of Molecular Carcinogenesis, Medical University of Lodz, 90-752 Lodz, Poland
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2023, 24(1), 164; https://doi.org/10.3390/ijms24010164
Submission received: 13 October 2022 / Revised: 14 December 2022 / Accepted: 17 December 2022 / Published: 22 December 2022
(This article belongs to the Special Issue Cellular and Molecular Mechanisms Involved in Prostate Cancer Cell Migration and Invasion)
Browse Figures
Review Reports Versions Notes

Abstract

:
Prostate cancer (PC) remains a worldwide challenge, as does the question of how to distinguish its indolent from its aggressive form to reconcile proper management of the disease with age-related life expectations. This study aimed to differentiate the Notch-driven course of PC regarding patients’ ages and stage of their disease. We analyzed 397 PC samples split into age subgroups of ≦55, 60–70, and >70 years old, as well as early vs. late stage. The clinical association of Notch signaling was evaluated by DFS and UpSet analyses. The clustering of downstream effectors was performed with ExpressCluster. Finally, for the most relevant findings, functional networks were constructed with MCODE and stringApp. The results have been validated with an independent cohort. We identified specific patterns of Notch expression associated with unfavorable outcomes, which were reflected by entering into a hybrid epithelial/mesenchymal state and thus reaching tumor plasticity with its all consequences. We characterized the molecular determinants of the age-related clinical behavior of prostate tumors that stem from different invasive properties depending on the route of the EMT program. Of the utmost relevance is the discovery of age- and stage-specific combinations of the Notch molecules predicting unfavorable outcomes and constituting a new prognostic and therapeutic approach for PCs.

1. Introduction

Prostate cancer (PC) is the most frequently diagnosed cancer type in men and leads the most cancer-related mortality rates right after lung cancer, globally [1]. Elderly men comprise the majority of PC cases; however, recently, an increasing trend is reported among young males [2]. The course of PC is very heterogeneous, posing major challenges in proper management of the disease and accurate risk stratification. The former is strongly related to the age of PC onset and requires individualized decisions determined by balancing the benefits and risks of active treatment while considering comorbid conditions, life expectancy, and tumor characteristics [3]. In turn, risk stratification aims to distinguish indolent from aggressive PCs of the advanced stage. The latter is much more frequently diagnosed among younger males and reaches the highest mortality among all age groups, hence early-onset PC having been already determined to be a distinct clinical–pathological phenotype with a poor prognosis [4]. In contrast, older men tend to be diagnosed with tumors of a more advanced stage, but, simultaneously, they also appear to die from other comorbid causes. Many reports show that the rate of disease is associated with increasing age, as a small microscopic focus, and with reduced aggressiveness [5]. The difficulties in adequate PC management to either reduce unnecessary overtreatment of more benign cases or apply efficient strategies for metastatic disease prompt an expansion of the current knowledge, with molecular mechanisms and hallmarks reflecting the age-specific course of prostate tumorigenesis ultimately aiding the stratification of high- from low-risk tumors. Such distinction lays in functional alterations of regulatory pathways driving the cancer cell behavior that is manifested in the clinical observations among individuals. Currently, the early detection of PC, as well as clinical decisions on the management, are mostly based on prostate-specific antigen (PSA) screening. However, PSA is rather a prostate- than prostate cancer-specific biomarker; thereby, the information provided by its measurement is often ambiguous. On the other hand, PSA has been considered a more accurate secondary indicator to monitor PC progression among already diagnosed PC patients [6].
Some of the pathways carry external signals from the environment or neighboring cells, others trigger the internal programs of cellular death and proliferation. Notch is one of the major developmental pathways orchestrating the whole cellular machinery during pre- and postnatal life: proliferation, differentiation, and survival, to determine the fate of the cells and tissue architecture [7]. Unsurprisingly, due to its roles, an overwhelming body of evidence has been collected on cancerogenic consequences of Notch deregulation [8], including on PC [9]. To date, expression profiling studies reveal that the core components of Notch signaling may distinguish aggressive PCs of high Gleason grade and indicate that the pathway is functional especially in progressing metastatic disease [10]. Although most of the Notch-dedicated studies consider the tumorigenic effects of the mutations occurring within the pathway or single-gene alterations, the aberrations of core signaling resulting in cancer progression due to affected downstream output without overt mutations in the core members of the Notch pathway seem of the greatest interest.
In this study, we attempted to differentiate the Notch-driven course of tumorigenesis and disease progression of PC regarding the age of onset as well as the stage of the disease. We aimed to construct a global view of Notch-related expression alterations that occur in the core components, but, more importantly, we also took the downstream consequences of the primary Notch deregulation into consideration to reveal the complex background of age-associated differences in PC course. Finally, we stressed the potential of Notch profiling exploitation in stratifying patients who truly need active treatment compared to those who can remain under surveillance.

2. Results

2.1. Decreased Notch Activity Restrains PC Progression in Older Patients

The expression profiling of the core Notch components revealed differences between the age groups of the same tumor stage as well as within the same age groups between the tumor stages (Figure 1). The age group of above 70 years old was distinct from the age groups of below 55 or 60–70 years old, which were similar to each other to some extent, especially regarding the late stage of the disease. A more detailed comparison of the Notch expression profiles between the tumor stages revealed the spectrum of Notch activity reaching a peak at 60–70 years old. At the early stage of PC among patients, the age group of under 55 years old expressed ligands (DLL1, DLL3, JAG1), Notch-specific TFs (HES1 and HEY1), and major components of the γ-secretase complex (PSEN1, PSEN2, PSENEN, NCSTN). We also observed increased activity of epigenetic regulators such as histone deacetylases (HDAC1, HDAC2), as well as cross-talk with the Wnt pathway through DVL3. At the late stage of PC in patients below 55 years of age, similarly to early and then late PC in 60–70-year-olds, we observed heightening of Notch activity through triggering the expression of additional ligands (DLL4, JAG2, DLK1), Notch receptors, TFs (HES4, HES5, HEY2, HEYL, PTCRA), and intensified cross-talk with the Wnt pathway through DVLs (DVL1, DVL2) until its activity extinction in elderly men, especially at the late stage of PC. However, we also observed disparate epigenetic regulation through histone deacetylases (HDAC1, HDAC2), CREBBP, and EP300, which switched on and off in different PC subgroups, accordingly (Table 1). Detailed results of the comparison may be found in the Supplementary Table S1 and GitHub repository (https://github.com/orzechmag/notch-age-pc).
Moreover, the core Notch signaling demonstrated the specific age- and tumor stage-related expression signatures enclosed in Supplementary Table S2. We identified a unique expression profile of 13 Notch genes among patients below 55 years old, comprising, among others, DLL3 and HEY1 at the early stage, whereas JAG2, DLL4, HEYL, NOTCH3, NOTCH4, and DVL1 were relevant at the late stage. Furthermore, 18 genes were characterizing tumors of 60–70-year-olds, such as HES1, HES5, and HEY2 at the early stage and NOTCH1, JAG1, RBPJ, NCSTN at the late stage. Finally, 37 genes were unique for tumors in elderly men, including JAG1, DLL4, DVL1, NOTCH1, NOTCH2, NOTCH3, PTCRA, NCSTN, and RBPJ at the early stage, while DLL3, DLK1 were relevant for the late stage of the disease. Significantly, 13 genes overlapped between groups of patients below 55 and 60–70 years old (early vs. late tumor stage DOWN-UP: KAT2B, NUMB, DVL2, MAML1, PTCRA, CREBBP, MAML2, NOTCH2, APH1B, ADAM17, RBPJL, UP-DOWN: SNW1, LFNG), while none of the genes intersected between these groups and the elderly men. The summary of the above comparisons is shown in Table 2.

2.2. The Specific Combination of Jagged Ligands and Notch Receptors Drives the Metastatic Potential of PC

We aimed to evaluate the effects of the particular Notch members profiles on DFS analysis among age groups, and stage of the disease, accordingly. Detailed results are shown in Table 3. Among the major findings, we observed contrasting effects of JAG1 and JAG2 on DFS with increasing age, although independently of the disease stage (JAG2: HR≦55 = 3.03, p = 0.028; HR60–70 = 6.12, p = 0.044; late-stage HR = 1.94, p = 0.043; HR>70 = 0.069, p = 0.012; JAG1: HR>70 = 0.079, p = 0.007; late-stage HR>70 = 0.093, p = 0.012). The same trend could be also observed among Notch receptors: NOTCH4 significantly altered the prognosis among the group below 55 years old (HR = 3.34, p = 0.023) and in patients above 70 years old (HR = 0.069, p = 0.012). In turn, NOTCH2 and NOTCH3 were associated with DFS in patients aged 60–70 years old (NOTCH2: HR = 1.89, p = 0.041; NOTCH3: HR = 2.37, p = 0.007; early-stage HR = 24.9, p = 0.0002), but they were not significant in the other age or stage groups. Furthermore, some of the Notch genes differentiated DFS with tumor stage within patients aged 60–70 years old, such as ADAM10, DTX3L, DVL3, and HDAC1, which were significantly lowered at the early stage (ADAM10: HR = 14.4, p = 0.0043; DTX3L: HR = 11.2, p = 0.016; DVL3: HR = 9.73, p = 0.023; HDAC1: HR = 9.4, p = 0.025) and heightened at the late stage (ADAM10: HR = 0.373, p = 0.011; DTX3L: HR = 0.515, p = 0.034; DVL3: HR = 0.299, p = 0.0019; HDAC1: HR = 0.398, p = 0.037) of PC. We also performed the UpSet analysis (form of Venn diagrams) to determine whether some of the PC patients bear a specific Notch signature, i.e., the intersect of the Notch core members profiles, that is associated with DFS prognosis (Supplementary File S3). Among 53 patients below 55 years old at the early stage of the disease, we identified a subgroup of 20 patients (37.7%) with a Notch signature consisting of RBPJ, ADAM17, RBPJL, and ATXN1; however, there was no significant difference in DFS between the possessors and non-possessors of the signature (HR < 0.001, p = 0.23, Figure 2). At the late stage of the disease, half of the group aged below 55 years old showed the signature composed of KAT2A, LFNG, HDAC2, SNW1, and PSEN2 corresponding to significantly better DFS rates (HR = 0.081, p = 0.0028, Figure 2). In the group of PC patients aged 60–70 years old at the early stage of the disease, 58 patients (69%) had a common Notch profile of NOTCH3, APH1A, HDAC2, MFNG, and DTX4, improving the DFS prognosis (HR < 0.001, p = 0.021, Figure 2), while 124 (73.4%) patients at the late stage of PC showed the profile of HES5, ADAM17, CREBBP, DVL3, and ADAM10 (HR = 0.35, p = 0.0016, Figure 2). Finally, we found a favorable Notch signature among 6 (27.3%) late-stage PC patients aged above 70 years old, consisting of JAG1, PTCRA, HEYL, DTX3, and RFNG (HR < 0.001, p = 0.013, Figure 2). The detailed statistics of the DFS analysis for carriers and non-carriers of the identified Notch signature are shown in Table 4.

2.3. Notch Signaling Targets Tumor Properties Such as Aggressiveness and Invasive Potential

Alterations within the core Notch signaling entail deregulation of the downstream biological processes executed by the targets of Notch-specific TFs, such as the HES and HEY families manifesting in the clinical characteristics of the tumor. By the clustering of the Notch effectors, we identified 188 unique profiles therein, within which we determined specific clusters of high biological relevance, distinctive between the age groups and stage of the disease (Table 5, Supplementary File S3). Then, we constructed a functional network of high confidence subclusters based on the topology to find densely interconnected regions representing intrinsic/entangled pathways or processes. It led us to unravel that Notch signaling targets PC properties such as aggressiveness and invasive potential, which are modified along with progressing age and stage of the disease. Though the genes predominating within the network were associated mainly with the cytoskeleton (actin microfilaments, intermediate filaments, microtubules, and molecular motors), adhesion, extracellular matrix (ECM), establishing cell polarity, histone modifications, and epithelial-to-mesenchymal transition (EMT)-specific transcription factors, as well as angiogenesis, apoptosis, and cell cycle, our investigations revealed four major subclusters of Notch—EMT regulatory loop, kinesin motors, integrin signaling, and cellular migration—as the highly interconnected players driving PC progression (Figure 3).

2.3.1. Reshaping of Intercellular Communication—Cytoskeleton, Basement Membrane, and Extracellular Matrix

Alterations within the architecture of actin microfilaments were mostly specific for the age group below 55 years old and concurrently contrasted to mutually similar groups of 60–70 and above 70 years old. In addition, these profiles reverted at the late stage vs. early stage of PC. In turn, the profiles of the components of the intermediate filaments, such as FLG (filaggrin) and VIM (vimentin), distinguished the group of 60–70-year-olds from all groups at the early stage of the disease, with lowered FLG expression and heightened VIM expression therein; whereas, with PC progression, elevated FLG expression was kept only among patients above 70 years old, and significantly, this group acquired VIM expression while the others maintained the early stage profile. Furthermore, the PC early-stage expression of the microtubules constituents was mostly increased in patients above 70 years old, whereas higher expression of the molecular motors-related genes distinguished younger patients below 55 years old. In contrast, at the late stage of PC, the microtubules and molecular motors genes showed unified patterns of expression. Particularly, patients below 55 years old demonstrated a decrease in expression, while the groups aged 60–70 and above 70 years old presented the rise. Among adhesion-related molecules, the expression of CDH1 (E-cadherin) was higher among patients below 55 years and above 70 years at the early stage of the disease and, significantly, the CDH1 expression was lost among patients below 55 years old with stage progression. In turn, CTNNB1 (β-catenin) and CTNNA2 (α-catenin) showed opposite profiles differentiating patients at both stages of the disease, although shifting in the patterns. Specifically, we observed an increase in the CTNNB1 expression among patients below 55 years, in contrast to the groups of 60–70 and above 70 years old, while CTNNA2 was decreased in the former and elevated in the latter at the early stage of PC. With disease progression, CTNNB1 and CTNNA2 switched their expression patterns to the opposite in comparison to the early PC stage (CTNNB1: downregulation in the group below 55 years old vs. upregulation in the groups of 60–70 and above 70 years old; CTNNA2: upregulation in the groups of below 55 and 60–70 years old vs. downregulation in the group below 55 years old).
Among ECM-related genes, collagens (COL4A4, COL6A2, COL8A1, COL14A1), integrins (ITGA1, ITGA6, ITGA8, ITGA10, ITGAV, ITGB2, ITGB8), as well as MMP7, were mostly overexpressed in patients below 55 years old, rather than in the groups of 60–70 and above 70 years old at the early stage of PC, and the profile switched to the opposite with stage progression.

2.3.2. Complementary Alterations—The Important Players

Of the EMT-specific transcription factors, SNAI2 maintained its expression profile along with stage progression, although it altered with age: it showed downregulation in patients below 55 years old, whereas it was increased in the groups of 60–70 and above 70 years old. TWIST1 expression was lowered in the groups below 55 and above 70 years old at the early stage of PC, and the former group additionally gained the TWIST1 expression at the late stage. ZEB2, in turn, demonstrated opposite profiles with the progression of PC, as it was upregulated in patients below 55 years old in contrast to the groups of 60–70 and above 70 years old at the early stage, but, at the late stage, the expression dropped in the former and heightened in the latter groups. Additionally, we observed similar patterns in the expression of antiapoptotic BCL2, cell cycle regulator and suppressor RB1, epigenetic modificators (CREBBP, EP300), and hypoxia-related transcription factor HIF1A. Finally, we found several other notable alterations of a different character. Particularly, CDK12 expression was heightened only among patients below 55 years old at the early stage and above 70 years old at the late stage of PC. In turn, CDKN1A was lowered among patients above 70 years old at the early stage. In contrast, the pattern switched into overexpression among patients below 55 years old and downregulation in the older groups with PC progression. Interestingly, TP63 showed the same pattern of expression independently of the disease stage, involving downregulation among patients below 55 years old compared to older patients. Similarly, COCH expression was heightened only among patients above 70 years old, independently of disease progression. ESR1 expression was, in turn, decreased in this group, compared with the younger patients at the early PC stage, although the profile switched to the opposite with stage progression.

2.3.3. Specified Pattern of Notch Signaling Escalates PC Invasiveness through Affecting Plasticity

At the end, we attempted to evaluate the tumor properties associated with aggressiveness and invasive potential among the patients bearing the specific Notch signatures identified in the age and stage subgroups via a functional network consisting of cytoskeleton, adhesion, and ECM-related molecules, as well as key participants of EMT (Figure 4, Table 6). The general finding is that in the PC subgroups with specific Notch signatures associated with unfavorable DFS, these components were upregulated at the early stage of the disease and downregulated at the late stage, unrelated to age. Specifically, at the early stage of the disease, among patients aged below 55 years old that shared the unfavorable profile of the Notch signature (RBPJ, RBPJL, ADAM17, and ATXN1), we observed mainly alterations involving the remodeling of cytoskeleton components: actin microfilaments ACTA2 (actin α2) and ANXA1 (annexin A1) were downregulated in favor of elevation of DST (dystonin), FLNB (filamin B), FN1 (fibronectin 1), MYO1F (myosin IF), and MYO7A (myosin VIIA); intermediate filaments: FLG (filaggrin), cytokeratins KRT19 and KRT5, as well as VIM (vimentin); microtubules: tubulins TUBGCP4 and TUBGCP5; molecular motors: dynein DYNC1H1 and major kinesins KIF1B, KIF21A, KIF21B, KIF26A, KIF2A, KIF3A, KIF3B, KIF5A, and KIFAP3. Regarding adhesion, we observed (high CDH2 (N-cadherin), CDH22, and low CDH1 (E-cadherin); increases in CTNNB1 (β-catenin) and CTNNA2 (α-catenin), GJA1 (gap junction protein a 1; connexin 43), and SDC1 (syndecan 1), decreases in GSC (goosecoid homeobox), and MUC1 (mucin 1). Regarding ECM, we observed increases in COL8A1, DDR2 (discoidin domain RTK2), FBN1 (fibrillin 1), FRAS1, integrins such as ITGA1, ITGA6, ITGA8, ITGAV, ITGB2, and ITGB6, laminin LAMA3, LGALS3 (galectin 3), matrix metalloprotease MMP3, NTN4 (netrin 4), and TNN (tenascin N), and a decrease in VTN (vitronectin), GDF5, ITGA10, and ITGA5. Finally, we identified alterations among the major TFs driving the EMT program, such as a 1.5-fold rise in SNAI2, TCF4, and ZEB2 expression, whereas FOXC2 and TWIST1 dropped. At the late stage of the disease, among patients bearing the unfavorable profile of the Notch signature (KAT2A, LFNG, HDAC2, SNW1, PSEN2), we observed that the major components of actin microfilaments were significantly downregulated (ACTA2, ANK2, DST, FN1), the same as myosins (MYO1F and MYO7A) and intermediate filaments including FLG, KRT5, KRT19, and VIM. Alterations in adhesion comprised decreases in E-cadherin (CDH1) and CDH11 and a striking drop in CDH22. In turn, α- and β-catenins (CTNNA2, CTNNB1) decreased compared to early-stage PC. On the other hand, we observed diminishing of the expression of ECM molecules (collagens, integrins, laminins, matrix metalloproteases) by at least 1.5 times at the late stage of the disease in patients bearing the unfavorable profile of the Notch signature. Significantly, ITGA10 and VTN were overexpressed, as opposed to their profiles from the early stage of the disease. We also observed noticeable decreases in expression of all major EMT-related TFs, especially SNAI1 and SNAI2. In patients agedbetween 60 and 70 years old at the early stage of the disease, the profiles of diverging subgroups of the DFS-related Notch signature (NOTCH3, APH1A, HDAC2, MFNG, DTX4) were more diverse than in patients under 50 years old. The components of actin microfilaments were mostly upregulated among patients sharing the unfavorable profile of the Notch signature, but to a lesser extent than in younger patients (e.g., FLNB, FN1, MYO7A). In turn, ACTA2 showed the same trend towards a decrease in expression. Opposite trends were observed for ANK2 and DST, which were noticeably downregulated, whereas both genes were overexpressed among younger patients (Table 6). Significantly, intermediate filaments (FLG, KRT19, KRT5, KRT8) were all downregulated, and, compared to younger patients, they switched to the adversative profiles. Among adhesion molecules, we observed a rise in CDH1 expression with a simultaneous decrease in the level of CDH2. Additionally, CDH15 and CDH22 lessened their expression in contrast to younger patients with the unfavorable profile of the Notch signature. Also, CTNNA2 showed downregulation accompanied by CTNNB1 elevation. Molecules such as DSP, GSC, ICAM1, ICAM5, MUC1, OCLN, and TJP1 showed overexpression and reached higher levels than in younger patients, in contrast to GJA1, NID1, and SDC1, which significantly dropped. The expression of the ECM components was more diverse, as some of them gained and some of them lost their expression. In general, collagens showed downregulation, unlike in younger patients. Genes such as DDR2, FBN1, and FRAS1 decreased their expression compared to younger patients. Integrins, laminins, matrix metalloproteases, NTN4, TNN, and TNXB decreased their expression, whereas ITGA10, LAMC2, and VTN increased contrastingly to patients under 55 years old. Finally, among TFs, SNAI2 showed a marked decrease in expression, whereas the remaining ones remained overexpressed. In the late stage of the disease, the carriage of the unfavorable profile of the Notch signature consisting of HES5, ADAM17, CREBBP, DVL3, ADAM10 resulted in downregulation of major adhesion and cytoskeleton-related molecules. Although all actin microfilaments were decreased, ANK2, CALD1, and FN1 dropped in the expression, thus resembling the late-stage PC in patients under 55 years old and half-reciprocal to the early-stage PC among 60–70-year-old patients. The intermediate filaments were convergent with both late-stage PC in patients below 55 years old and early-stage PC of 60–70year-old patients, showing downregulation. Besides, microtubules and molecular motors were downregulated and contrasted to early-stage PC among 60–70-year-old patients, as well as late-stage PC in patients below 55 years old. Among adhesion and ECM molecules, we also observed a general decrease in expression, hence demonstrating a unique pattern dissimilar to the comparative age and stage groups. Unlike early-stage PC among 60–70-year-old patients and late-stage PC in patients agedbelow 55 years old, CDH2 and CDH4, as well as LIN7A, COL8A1, DSP, GSC, ICAM1, ICAM5, FRAS1, ITGA10, ITGA8, ITGAV, OCLN, PLEC, TJP1, and NPNT, showed a drop in expression, whereas CDH15, CDH26, and CTNNA2 were overexpressed, akin to late-stage PC in patients below 55 years old. Regarding TFs, only TWIST1 was increased. Such a pattern was also revealed for late-stage PC in patients below 55 years old. Among patients over 70 years old, it was only possible to analyze the late-stage disease, where we identified the unfavorable Notch signature of JAG1, PTCRA, HEYL, DTX3, and RFNG. The carriage of the unfavorable profile resulted in lowered expression of actin microfilaments, whereas MYO1F was elevated. The same trend was observable in intermediate filaments, especially for FLG, as well as for microtubules and molecular motors. In addition, adhesion molecules were lowered, excepting CDH11, DSP, GSC, MUC1, S100A4, and SDC1, similarly to ECM molecules, among which only ITGB2, LAMA1, LAMC3, MMP9, and TNN showed overexpression. Finally, all TFs were downregulated; however, TWIST1 was higher, resembling the late-stage PC of the 60–70-year-old patients.

2.4. Validation of the Findings

For the cross-validation of our findings, we employed the independent microarray study of Rubicz et al., which resembles the design of our investigations with PC staging (local and regional disease) as well as the patients’ age (≤55, 60–70, and 70> years old). We performed the validation at two levels. Firstly, the comparative profiling of the core Notch members was performed and revealed approx. 47% similarity between the expression patterns in the primary and cross-validation studies. Importantly, we also confirmed shifts in the expression of particular Notch members between disease stages of two (SNW1, DLL3), three (PSENEN, NUMBL, DTX2), and two (NCSTN, EP300) losses among patients below 55 years old, 60–70 years old, and above 70 years old, respectively, whereas five (KAT2B, JAG2, HEYL, MAML2, NOTCH2), two (NUMB, NCSTN), and four (PSENEN, SNW1, DLL3, KAT2A) genes gained the expression with disease progression, respectively. The comparison of the general trends in expression patterns is presented in Table 2 and Supplementary File S2. Subsequently, we conducted comparative clustering of the Notch targets, and, as shown in the right section of Table 5 and Supplementary File S3, we successfully confirmed many of the primary findings, but there were also some discrepancies.

3. Discussion

A well-known and widely used prostate cancer Gleason diagnosis system reflects pathomorphological differentiation of cancer tissue. Such microscopically classified cellular clonality is directly associated with changes in cell molecular activities, meaning differential regulation of gene expression. The Notch signaling pathway is one out of several main regulators of cellular functional differentiation and tissue behavior. In this study, we demonstrated transcriptomic representation of the effects of Notch signaling on PC course and progression according to age of disease onset and tumor stage. We observed that the activity of the core Notch pathway is significantly modulated by both factors, following the natural route of cellular aging, which may be compared to the light dimmer. On one hand, it leads to progressive loss of organic functions and self-renewal capacity in humans [11]. On the other hand, carcinogenic transformation and progression requires uncontrolled proliferation, restrained apoptosis, and the gaining of new functions [12,13]. These seeming dissimilarities are, in fact, two ends of the same stick of common molecular pathways that cause, on one side, loss of organ functionality and, on the other side, confer pro-survival functions at earlier stages of life [14]. In line with the above, we revealed an increase in Notch signaling activity among patients under 55 years of age that intensified with more advanced-stage PC involving the most essential Notch members, such as receptors (NOTCH2-4), ligands (JAG2, DLL4), and HEYL, a Notch-specific TF (Figure 1). Among patients aged 60–70 years old, the core Notch pathway was still upregulated, although of a different pattern including NOTCH1-2 receptors, JAG1 ligand, NCSTN regulator, as well as Notch-related epigenetic transcription modulators such as HDAC1, CREBBP, and EP300 (Figure 1). Among the oldest patients, we revealed dimming of the Notch activity through the decrease of expression in its main elements (JAG1, NOTCH1-3, CREBBP, EP300, γ-secretase complex) in favor of upregulation of Notch modulators (RFNG, LFNG, HDAC2) (Figure 1). In turn, late-stage PC did not reveal any significant Notch activity, which resembles the Goldilocks Principle, encompassing neither too low, nor too high activity of the pathway to maintain homeostasis [15]. These findings were supported with specific Notch profiles affecting DFS that, importantly, switched with progressing age and stage of the PC patients (Figure 2, Table 3). These observations are, to a large extent, consistent with the current state of knowledge, apart from influence of age and stage, which, as yet, has not been widely considered. For instance, in vitro experiments using PC cell lines revealed that NOTCH1 activity was promoting cancer cell migration and invasion as well as augmenting aggressiveness of the tumor [16,17]. In humans, a raise in JAG1 was associated with metastatic PC lesions [18,19], similarly to NOTCH3, the expression of which was inversely correlated with survival [20]. Fundamental Notch members were also upregulated in PCs of Gleason 8 compared to the tumors of lower grades [10]. It is worth emphasizing that, among younger patients, we observed an asymmetry between Jagged and Delta ligands, giving cells a rise to a hybrid state described in a series of articles by Boareto et al. In brief, due to Delta–Notch signaling, cells may acquire one of two possible states: Sender (high expression of ligands) or Receiver (high expression of receptors) [21]. However, Delta–Jagged asymmetry enables the cells to acquire a hybrid Sender/Receiver state that implicates a worse clinical outcome. It was shown that during angiogenesis, the asymmetry leads to the hybrid tip/stalk phenotype allowing chaotic and fast sprouting to ensure efficient oxygen supply to rapidly growing tumors [22]. Moreover, it is believed that the hybrid phenotype perpetuates the meta-stable hybrid epithelial/mesenchymal phenotype linked with more aggressive tumor behavior and stem cell-like features [23]. Remarkably, this evidence is a molecular substantiation of the observations of increased PC aggressiveness in younger males and is a proof of being a separate biological entity.
The phenotypic changes of PC induced by Notch are mediated by a transcriptional program of its numerous downstream effectors, including genes linked to invasive properties of cancer that we presented in a previous study [24]. The remarkable changes associating disease progression involve remodeling of cell–cell and cell–ECM interactions such as reshaped adhesion, cell junctions, cytoskeleton components (actin microfilaments, intermediate filaments, microtubules, and molecular motors), protease activity, constituents of the ECM, and differentiation factors determining the fate of the cell and its stem state-related properties [25]. Currently, rapidly growing data confirm the clinical premises of the more aggressive course of early-onset (herein diagnosed at below 55 years of age) PC compared to that in elderly men, which apparently arise from the distinct biology of the tumors [26,27,28,29]. Prior work has implicated age-related specific profiles of EMT, which, at molecular scale, reflected a Gleason grade [30]. In the current research, we performed a much wider description of how age alters tumor invasiveness at cellular resolution; in particular, we emphasized which of these changes are orchestrated by the Notch activity. Each age group showed specific profiles determining the course of PC, as well as a disparate trajectory of evolution associated with stage progression, which, in network analysis, emerged as the four most tightly interconnected clusters (Figure 3). In fact, they could be identified with biological groups because they belong to such as the EMT program, molecular motors, and parts of cell–ECM adhesion (Figure 3).
Several attempts have been made to unveil genomic landscapes of PC in younger and elderly men, although they conferred mostly mutations or signature genes [31,32]. There are well-known risk factors of PC identified as mutations of genes BRCA1, BRCA2, mismatch repair genes (MSH2, MSH6, MLH1, PMS2, EPCAM) and HOXB13 gene coding transcription factor. In our studied cohort of PC patients, the highest mutation frequency has BRCA2 (5%), BRCA1 (2.2%), HOXB13 (1.6%) and PMS2 (1%), and all others have below 1% mutations. However, these mutation profiles were not related to the molecular and age groups we identified. Conversely, as shown in Figure 5, we deciphered that aggressiveness of early-onset PC originates from reshaping of the cytoskeleton, adhesion, and ECM accompanying EMT-related plasticity of the tumor cells, which was modulated across the progressing stage of PC. The above supports gain in actin microfilaments such as caldesmon (CALD1) and myosins (MYO1D, MYO1F, MYO7A) participating in, e.g., assembling the stress fibers and motility-enabling protrusions [33]. Younger patients also showed high activity of microtubule-related molecular motors, i.e., kinesins and dynein, likely arising from uncontrolled cell division, as well as hypoxia-inducible factor 1-α (HIF1A). To date, several members of the kinesins family (KIF11, KIF15, KIF18B) have been correlated with poorer prognosis and aggressive forms of PC [34,35,36]. In turn, hypoxia in tumors arises from the metabolic requirements of rapidly dividing cells, which cannot be met due to ineffective neovascularization. However, cancer cells tend to preferentially metabolize glucose over lactate, despite aerobic conditions, known as the Warburg effect. Under hypoxic conditions, HIF1A induces, among others, angiogenesis, the program of EMT, and proliferation, as well as the migration of tumor cells. It is also a key regulator of enzymes involved in aerobic glycolysis [37]. The overexpression of HIF1A was, for instance, implicated in the biochemical progression of PC [38], as well as in developing resistance to hormonal therapies in castrate-resistant PC [39].
EMT has been repeatedly demonstrated in the progression and metastasis of many human malignancies [40]. During EMT, epithelial cells gain migratory capacities due to loosening junctions and cell–cell adhesion, as well as remodeling of the cytoskeleton to establish spindle-like morphology. EMT is mainly identified by a loss of major epithelial markers (cytokeratins, E-cadherin) to the detriment of acquiring mesenchymal markers (N-cadherin, vimentin, fibronectin). Importantly, EMT is more and more often considered a hybrid spectrum of epithelial and mesenchymal states, in fact translating to the plasticity that leads to stem cell-like, aggressive cancer phenotypes [41,42,43]. Such cells attaining a hybrid epithelial/mesenchymal phenotype exhibit both epithelial and mesenchymal features (i.e., they are capable of migrating with restricted cell–cell adhesion), which might be a determinant of collective cell migration of multicellular aggregates in the ECM or clustered circulating tumor cells found in the bloodstream of PC patients [44]. Consequently, cells expressing markers of both states have much higher metastatic potential due to enhanced resistance to cell death or more efficient extravasation rates. Thus, this hybrid state resembles the “Mr. Hyde” personality that constitutes much higher metastatic risk in patients than a complete EMT phenotype or “Dr. Jekyll” side [41]. It was further supported by several studies reporting, among others, co-expression of both E-cadherin and vimentin in, e.g., invasive breast cancer, that when exhibited by the tumors concomitantly conferred the worst DFS and OS outcomes [45]. The process of EMT itself may be regulated by many mechanisms and signaling pathways, such as Notch signaling [46], which modulate the expression of EMT-triggering TFs. To date, various members of the EMT program have been reported in PC research in association with phenotypic disparities of the tumor [47]. For instance, Cheaito et al. demonstrated that co-expression of cytokeratin 8 and vimentin correlated with higher Gleason grades and worse prognosis of biochemical recurrence-free survival [48]. In the review devoted to PC plasticity, Papanikolaou and collaborators demonstrated similar findings that were reported on N-cadherin, thus guiding the cross-talk between the tumor stroma and its epithelial cells. EMT has been also associated with the androgen-signaling axis, which is of very high importance in PC, specifically. AR has been considered a transcriptional repressor of E-cadherin in a manner comparable to Snail and Twist TFs, as well as an activator of β-catenin; hence, AR attributed to more mesenchymal phenotypes of PC cells [49].
Figure 5 presents the major findings and how the observed alterations translate into biological consequences for the cell (dark pink elements). In the present study, we noted concurrent overexpression of E-cadherin (CDH1) and β-catenin (CTNNB1) at the early stage of PC among younger males, which remains in line with our previous studies [30]. It also suggests the phenotype of tumor cells that only partially execute the EMT program, i.e., contemporary expression of both epithelial (CDH1, CTNNA2, LAMA1) and mesenchymal (FN1, CTNNB1, MMP7, ZEB2) markers. At the late stage, the cells lost the aforementioned characteristics, thus suggesting a decrease in the invasive potential compared to early-stage tumor cells in men below 55 years of age. In addition, at particular stages of the disease, patients 60–70 and above 70 years old showed very similar features, mostly contrasting to those of patients below 55 years old. Importantly, late-stage tumors of older patients resembled a more aggressive form of PC than in younger patients. Specifically, patients 60–70 years old showed high levels of fibronectin (FN1), β-catenin (CTNNB1), matrix metalloproteases (MMP2, MMP7), vimentin (VIM), microtubules and molecular motors, and EMT-triggering TFs such as SNAI2, ZEB2, TWIST1, HIF1A, and SMAD3, whereas the expression of E-cadherin (CDH1) was lost. Remarkably, these are the evidence for the predominant mesenchymal phenotype of the tumor cells that executed EMT and acquired migratory potential. Among the oldest patients, we noted heightened E-cadherin (CDH1) with concurrent β-catenin (CTNNB1), fibronectin (FN1), vimentin (VIM), microtubules and associated molecular motors, and EMT-related TFs such as SNAI2, ZEB2, HIF1A, and SMAD3 (detailed results are enclosed in Table 5). Notably, the similar trends observed in the cross-validation study comparing local and regional PCs emphasize the importance of the findings that clearly suggest disparate molecular mechanisms of prostate tumor invasion and progression in association with the age of PC onset. These observations have a number of implications for understanding the clinical behavior of PC diagnosed at different ages and stages, as well as highlight the role of Notch in influencing the plasticity of prostate tumors through modulation of its downstream effectors. These properties are an elicitation of disease progression through executing the sequential program, enabling the cancer cells’ detachment from the primary mass [50].
Much evidence indicates Notch signaling as a potential therapeutical target, and it is currently under clinical trials in some malignancies [51]. However, the reported Notch deregulation and its associations with various aspects of prostate tumorigenesis boils down to specific interactions and single members of the core signaling. Formerly, we confirmed that the carcinogenic effects of Notch are context-specific, i.e., they depend on the intensity of signaling, the pattern of the pathway members expression, and tissue type. The most important finding revealed that the characteristics of each solid tumor, including PC, are represented by specific profiles of the Notch pathway at two levels of signaling: the core and expression of the downstream effectors orchestrated by Notch-specific TFs (HES/HEY families) [52]. In the current study, we correlated the expression of the individual Notch participants with the DFS prognosis among specific age groups, as well as in combination with the PC stage, followed by the evaluation of the downstream effects among Notch targets. Though the current state of knowledge was lacking such findings, of the highest interest was whether the expressions of Notch participants combine into a specific pattern of activation translating into aberrant downstream mechanisms driving aggressive clinical behavior and thus could be identified as a PC-related prognostic factor of unfavorable prognosis. Ultimately, among each group of PC patients, we identified a repeating expression pattern of Notch, translating into significant differences in DFS that were detected in approx. half of the group (Figure 2, Table 3, Supplementary File S3). However, regardless of the known role that each molecule plays during transducing Notch signals, very little is known about their contribution in prostate tumorigenesis and progression. The closest idea to our study was research conducted by Kwon and collaborators that generated a Notch signature score involving the expression of receptors (NOTCH1-4), ligands (JAG1-2, DLL1, DLL3, DLL4), DTX1, and Notch-specific TFs (HES1-2, HEY1-2). It was demonstrated that the lower the Notch signature, the more favorable the prognosis of disease recurrence (i.e., PSA recurrence-free survival) [53]. In contrast, our signatures included specific mixtures of expression patterns that correlated with better DFS, i.e., RBPJ, RBPJL, ADAM17, ATNX1, and KAT2A, LFNG, HDAC2, SNW1, PSEN2 for patients below 55 years old at the early and late stages of PC, respectively; NOTCH3, APH1A, HDAC2, MFNG, DTX4 and HES5, ADAM17, CREBBP, DVL3, ADAM10 for patients 60–70 years old at the early and late stages of PC, respectively; finally, JAG1, PTCRA, HEYL, DTX3, RFNG for patients above 70 years old at the late stage of the disease (Figure 2, Table 3, Supplementary File S3). Hence, it suggests age-diversified mechanisms of Notch activation drive prostate progressions that reach beyond the fundamental transducers of Notch signaling and shed light on its modulators, affecting the strength of signaling. Following this, we observed that the unfavorable profiles of each signature involved escalated forms of downstream mechanisms of invasion and tumor plasticity. Importantly, they constituted their adversities, along with the disease stage, at the same age of PC onset (Figure 4) and remain consistent with the previous findings of cluster analysis.

4. Materials and Methods

For the study, we obtained the RNA-seq expression data with matched clinical information of 502 prostate adenocarcinoma (PC) patients provided by The Cancer Genome Atlas (TCGA). The data were downloaded via the NCBI Gene Expression Omnibus (GEO) series GSE62944 of Rahman et al. that published alternatively reprocessed and compiled RNA-seq with corresponding clinical data for the TCGA samples, with the Rsubread package offering improved analytic performance (normalized TPM values for tumor samples; data status of 27 January 2015) [54]. Then, we restricted the cohort to the specific age groups of ≤55, 60–70, and 70> years old, constituting a total of 397 patients. Additionally, our considerations were focused on the comparison between early-stage tumors restricted to the prostate gland, and late, locally advanced, or metastasizing stage tumors, which were classified with the American Joint Committee on Cancer (AJCC) 8th Edition TNM staging system [55] and subsequently grouped into early (stages I–II) or late (stages III–IV) stage subgroups. The detailed clinical characteristics of the study cohort are shown in Table 7.
The Notch signaling pathway and its core participants in humans were determined according to the KEGG database (hsa04330) [56] and MSigDB v7.2 [57]. The downstream effectors of the Notch signaling comprised a total of 10,740 targets of Notch-specific transcription factors (TFs) belonging to the HES and HEY families (HES1, HES4, HES5, HEY1, HEY2, HEYL) identified through the GTRD database v19.10 [58,59]. The identified target genes haven been annotated with their biological roles and divided into functional groups according to the current knowledge and literature, as well as our previous studies [24,52].
The expression profiling of the Notch core signaling between the groups of patients, according to the combination of age and tumor stage, was performed based on hierarchical clustering of the median gene expression in particular groups with the Pearson distance metric and the complete linkage method. The clustering was performed and visualized with the NMF R package with aheatmap() function. Additionally, to support the hypothesis that Notch signaling differentiates the clinical course of PC among the age groups and in association with the tumor stage, we performed disease-free survival (DFS) analysis regarding the effects of the expression of particular Notch members. The analysis was performed with the EvaluateCutpoints RShiny app through the maxstat algorithm [60]. In addition, we aimed to identify the subgroups of PC patients bearing specific combinations of the Notch core gene profiles associated with DFS prognosis by applying the UpSetR algorithm (UpSetR R package) [61]. For that purpose, we aggregated patients based on intersections of the dummy-encoded expression of particular Notch genes related to DFS outcomes that revealed specific patterns among PC patients. General DFS analysis was subsequently performed to compare the survival of PC patients with and without an identified Notch signature (log-rank test, p < 0.05; survminer R package). Furthermore, ExpressCluster software (http://cbdm.hms.harvard.edu/) was used to find common and unique expression profiles of the downstream Notch effectors among the age groups combined with stage of the disease. Clustering was performed by applying the K-means++ algorithm, z-norm (mean = 0, var = 1) signal transformation, rank correlation distance metric, 1000 iterations, and 400 clusters (K), as recommended for six class comparison. Profiles indicating unique contrasts between age groups within the tumor stage, as well as between tumor stages within the age group, were considered as significant. Finally, by using the relevant clusters in terms of functional association with age-related PC course and progression, we constructed the biological networks of interactions with stringApp in Cytoscape [62]. The most pertinent network that we successfully constructed was followed by MCODE identification of highly interconnected regions therein [63]. Moreover, the downstream biological differentiation between the subgroups of patients with specific patterns of Notch signaling was demonstrated by using biological networks from the previous step, with plotted changes in the expression between stage-related groups within the age groups of PC patients.
The findings of the present study have been cross-validated with an independent microarray study deposited in the NCBI GEO (GSE141551) aiming to evaluate the transcriptome-based genetic profiles of 503 localized prostate cancers (PC) associated with the aggressive course concerning age and stage (local vs. regional) of the disease [64]. Due to sample restriction according to the age groups of interest (≤55, 60–70, and 70> years old) we finally employed 380 individuals for validation purpose of the primary findings (local stage, ≤55 years old: 108; regional stage, ≤55 years old: 49; local stage, 60–70 years old: 129; regional stage, 60–70 years old: 69; local stage, 70> years old: 17; regional stage, 70> years old: 8). The validation was performed at two levels: profiling of the Notch core and ExpressCluster analysis of expression profiles of the downstream Notch effectors. Both analyses were performed analogously to the primary ones.

5. Conclusions

In a conclusion, PC is a worldwide medical challenge. The heterogeneous course of disease is conditioned by multiple factors, including the age of onset. This factor requires special attention, as aging is considered the most important constituent for cancer development and, importantly, determines decisions on the treatment strategy applied. While distinguishing the aggressive from the indolent form of PC is problematic, we characterized the molecular determinants of the age-related clinical behavior of the tumors. Primarily, it stems from different invasive properties depending on the route of the EMT program and adherent remodeling of cellular architecture. The Notch pathway, one of the most conservative mechanisms of signaling, through specific patterns of activation, orchestrates transcriptional programs corresponding to PC outcomes. While excessive Notch signaling in younger patients predestines more aggressive forms of PC, insufficiency in the signaling worsens the prognosis in elderly men. Of utmost relevance is the discovery of age- and stage-specific combinations of the Notch molecules, herein described as Notch signatures, which predict unfavorable outcomes arising from escalated invasive features, indicating an undifferentiated, stem cell-like state of the tumor cells and could constitute a new prognostic and therapeutic approach for prostate cancers.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms24010164/s1.

Author Contributions

M.J.O.: conceptualization, data curation, formal analysis, funding acquisition, investigation, methodology, software, validation, visualization, writing-original draft; D.A.: conceptualization, data curation, formal analysis, investigation, methodology, software, validation, visualization, writing—original draft; A.K.B.: conceptualization, supervision, writing—review & editing. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the National Science Center, Poland 2016/23/N/NZ2/02372.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets used for this study can be found in the Gene Expression Omnibus (accession number: GSE62944 and GSE141551; https://www.ncbi.nlm.nih.gov/geo/; accessed on 10 March 2021). Supplementary Materials supporting the presented findings are available in the GitHub repository (https://github.com/orzechmag/notch-age-pc).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef]
  2. Zhou, C.K.; Check, D.P.; Lortet-Tieulent, J.; Laversanne, M.; Jemal, A.; Ferlay, J.; Bray, F.; Cook, M.B.; Devesa, S.S. Prostate Cancer Incidence in 43 Populations Worldwide: An Analysis of Time Trends Overall and by Age Group. Int. J. Cancer 2016, 138, 1388–1400. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Stangelberger, A.; Waldert, M.; Djavan, B. Prostate Cancer in Elderly Men. Rev. Urol. 2008, 10, 111–119. [Google Scholar] [PubMed]
  4. Bleyer, A.; Spreafico, F.; Barr, R. Prostate Cancer in Young Men: An Emerging Young Adult and Older Adolescent Challenge. Cancer 2020, 126, 46–57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Scardino, P.T. Early Detection of Prostate Cancer. Urol. Clin. N. Am. 1989, 16, 635–655. [Google Scholar] [CrossRef]
  6. Ferraro, S.; Bussetti, M.; Panteghini, M. Serum Prostate-Specific Antigen Testing for Early Detection of Prostate Cancer: Managing the Gap between Clinical and Laboratory Practice. Clin. Chem. 2021, 67, 602–609. [Google Scholar] [CrossRef]
  7. Siebel, C.; Lendahl, U. Notch Signaling in Development, Tissue Homeostasis, and Disease. Physiol. Rev. 2017, 97, 1235–1294. [Google Scholar] [CrossRef] [Green Version]
  8. Aster, J.C.; Pear, W.S.; Blacklow, S.C. The Varied Roles of Notch in Cancer. Annu. Rev. Pathol. 2017, 12, 245–275. [Google Scholar] [CrossRef] [Green Version]
  9. Carvalho, F.L.F.; Simons, B.W.; Eberhart, C.G.; Berman, D.M. Notch Signaling in Prostate Cancer: A Moving Target: Notch Signaling in Prostate Cancer. Prostate 2014, 74, 933–945. [Google Scholar] [CrossRef] [Green Version]
  10. Ross, A.E.; Marchionni, L.; Vuica-Ross, M.; Cheadle, C.; Fan, J.; Berman, D.M.; Schaeffer, E.M. Gene Expression Pathways of High Grade Localized Prostate Cancer: Pathways of High Grade Prostate Cancer. Prostate 2011, 71, 1568–1577. [Google Scholar] [CrossRef]
  11. López-Otín, C.; Blasco, M.A.; Partridge, L.; Serrano, M.; Kroemer, G. The Hallmarks of Aging. Cell 2013, 153, 1194–1217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Hanahan, D.; Weinberg, R.A. The Hallmarks of Cancer. Cell 2000, 100, 57–70. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Hanahan, D.; Weinberg, R.A. Hallmarks of Cancer: The next Generation. Cell 2011, 144, 646–674. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Zinger, A.; Cho, W.C.; Ben-Yehuda, A. Cancer and Aging—The Inflammatory Connection. Aging Dis. 2017, 8, 611–627. [Google Scholar] [CrossRef] [Green Version]
  15. Braune, E.-B.; Lendahl, U. Notch—A Goldilocks Signaling Pathway in Disease and Cancer Therapy. Discov. Med. 2016, 21, 189–196. [Google Scholar]
  16. Bin Hafeez, B.; Adhami, V.M.; Asim, M.; Siddiqui, I.A.; Bhat, K.M.; Zhong, W.; Saleem, M.; Din, M.; Setaluri, V.; Mukhtar, H. Targeted Knockdown of Notch1 Inhibits Invasion of Human Prostate Cancer Cells Concomitant with Inhibition of Matrix Metalloproteinase-9 and Urokinase Plasminogen Activator. Clin. Cancer Res. 2009, 15, 452–459. [Google Scholar] [CrossRef] [Green Version]
  17. Kashat, M.; Azzouz, L.; Sarkar, S.H.; Kong, D.; Li, Y.; Sarkar, F.H. Inactivation of AR and Notch-1 Signaling by MiR-34a Attenuates Prostate Cancer Aggressiveness. Am. J. Transl. Res. 2012, 4, 432–442. [Google Scholar]
  18. LaTulippe, E.; Satagopan, J.; Smith, A.; Scher, H.; Scardino, P.; Reuter, V.; Gerald, W.L. Comprehensive Gene Expression Analysis of Prostate Cancer Reveals Distinct Transcriptional Programs Associated with Metastatic Disease. Cancer Res. 2002, 62, 4499–4506. [Google Scholar]
  19. Yu, Y.; Zhang, Y.; Guan, W.; Huang, T.; Kang, J.; Sheng, X.; Qi, J. Androgen Receptor Promotes the Oncogenic Function of Overexpressed Jagged1 in Prostate Cancer by Enhancing Cyclin B1 Expression via Akt Phosphorylation. Mol. Cancer Res. 2014, 12, 830–842. [Google Scholar] [CrossRef] [Green Version]
  20. Hudson, R.S.; Yi, M.; Esposito, D.; Watkins, S.K.; Hurwitz, A.A.; Yfantis, H.G.; Lee, D.H.; Borin, J.F.; Naslund, M.J.; Alexander, R.B.; et al. MicroRNA-1 Is a Candidate Tumor Suppressor and Prognostic Marker in Human Prostate Cancer. Nucleic Acids Res. 2012, 40, 3689–3703. [Google Scholar] [CrossRef]
  21. Boareto, M.; Jolly, M.K.; Lu, M.; Onuchic, J.N.; Clementi, C.; Ben-Jacob, E. Jagged-Delta Asymmetry in Notch Signaling Can Give Rise to a Sender/Receiver Hybrid Phenotype. Proc. Natl. Acad. Sci. USA 2015, 112, E402–E409. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Boareto, M.; Jolly, M.K.; Ben-Jacob, E.; Onuchic, J.N. Jagged Mediates Differences in Normal and Tumor Angiogenesis by Affecting Tip-Stalk Fate Decision. Proc. Natl. Acad. Sci. USA 2015, 112, E3836–E3844. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Bocci, F.; Gearhart-Serna, L.; Boareto, M.; Ribeiro, M.; Ben-Jacob, E.; Devi, G.R.; Levine, H.; Onuchic, J.N.; Jolly, M.K. Toward Understanding Cancer Stem Cell Heterogeneity in the Tumor Microenvironment. Proc. Natl. Acad. Sci. USA 2019, 116, 148–157. [Google Scholar] [CrossRef] [PubMed]
  24. Orzechowska, M.; Anusewicz, D.; Bednarek, A.K. Functional Gene Expression Differentiation of the Notch Signaling Pathway in Female Reproductive Tract Tissues-A Comprehensive Review With Analysis. Front. Cell Dev. Biol. 2020, 8, 592616. [Google Scholar] [CrossRef]
  25. Martin, T.A.; Ye, L.; Sanders, A.J.; Lane, J.; Jiang, W.G. Cancer Invasion and Metastasis: Molecular and Cellular Perspective; Landes Bioscience: Austin, TX, USA, 2013. [Google Scholar]
  26. Salinas, C.A.; Tsodikov, A.; Ishak-Howard, M.; Cooney, K.A. Prostate Cancer in Young Men: An Important Clinical Entity. Nat. Rev. Urol. 2014, 11, 317–323. [Google Scholar] [CrossRef] [Green Version]
  27. Zheng, Y.; Lin, S.X.; Wu, S.; Dahl, D.M.; Blute, M.L.; Zhong, W.-D.; Zhou, X.; Wu, C.-L. Clinicopathological Characteristics of Localized Prostate Cancer in Younger Men Aged ≤ 50 Years Treated with Radical Prostatectomy in the PSA Era: A Systematic Review and Meta-Analysis. Cancer Med. 2020, 9, 6473–6484. [Google Scholar] [CrossRef]
  28. Hussein, S.; Satturwar, S.; Van der Kwast, T. Young-Age Prostate Cancer. J. Clin. Pathol. 2015, 68, 511–515. [Google Scholar] [CrossRef] [Green Version]
  29. Chalmers, Z.R.; Burns, M.C.; Ebot, E.M.; Frampton, G.M.; Ross, J.S.; Hussain, M.H.A.; Abdulkadir, S.A. Early-Onset Metastatic and Clinically Advanced Prostate Cancer Is a Distinct Clinical and Molecular Entity Characterized by Increased TMPRSS2-ERG Fusions. Prostate Cancer Prostatic. Dis. 2021, 24, 558–566. [Google Scholar] [CrossRef]
  30. Jedroszka, D.; Orzechowska, M.; Hamouz, R.; Gorniak, K.; Bednarek, A.K. Markers of Epithelial-to-Mesenchymal Transition Reflect Tumor Biology According to Patient Age and Gleason Score in Prostate Cancer. PLoS ONE 2017, 12, e0188842. [Google Scholar] [CrossRef] [Green Version]
  31. Gerhauser, C.; Favero, F.; Risch, T.; Simon, R.; Feuerbach, L.; Assenov, Y.; Heckmann, D.; Sidiropoulos, N.; Waszak, S.M.; Hübschmann, D.; et al. Molecular Evolution of Early-Onset Prostate Cancer Identifies Molecular Risk Markers and Clinical Trajectories. Cancer Cell 2018, 34, 996–1011.e8. [Google Scholar] [CrossRef] [Green Version]
  32. Sharad, S.; Allemang, T.C.; Li, H.; Nousome, D.; Ku, A.T.; Whitlock, N.C.; Sowalsky, A.G.; Cullen, J.; Sesterhenn, I.A.; McLeod, D.G.; et al. Age and Tumor Differentiation-Associated Gene Expression Based Analysis of Non-Familial Prostate Cancers. Front. Oncol. 2021, 10, 584280. [Google Scholar] [CrossRef] [PubMed]
  33. Tojkander, S.; Gateva, G.; Lappalainen, P. Actin Stress Fibers—Assembly, Dynamics and Biological Roles. J. Cell Sci. 2012, 125, 1855–1864. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Wu, Y.-P.; Ke, Z.-B.; Zheng, W.-C.; Chen, Y.-H.; Zhu, J.-M.; Lin, F.; Li, X.-D.; Chen, S.-H.; Cai, H.; Zheng, Q.-S.; et al. Kinesin Family Member 18B Regulates the Proliferation and Invasion of Human Prostate Cancer Cells. Cell Death Dis. 2021, 12, 302. [Google Scholar] [CrossRef] [PubMed]
  35. Piao, X.-M.; Byun, Y.J.; Jeong, P.; Ha, Y.-S.; Yoo, E.S.; Yun, S.J.; Kim, W.-J. Kinesin Family Member 11 MRNA Expression Predicts Prostate Cancer Aggressiveness. Clin. Genitourin. Cancer 2017, 15, 450–454. [Google Scholar] [CrossRef]
  36. Qureshi, Z.; Ahmad, M.; Yang, W.-X.; Tan, F.-Q. Kinesin 12 (KIF15) Contributes to the Development and Tumorigenicity of Prostate Cancer. Biochem. Biophys. Res. Commun. 2021, 576, 7–14. [Google Scholar] [CrossRef]
  37. Courtnay, R.; Ngo, D.C.; Malik, N.; Ververis, K.; Tortorella, S.M.; Karagiannis, T.C. Cancer Metabolism and the Warburg Effect: The Role of HIF-1 and PI3K. Mol. Biol. Rep. 2015, 42, 841–851. [Google Scholar] [CrossRef]
  38. Vergis, R.; Corbishley, C.M.; Norman, A.R.; Bartlett, J.; Jhavar, S.; Borre, M.; Heeboll, S.; Horwich, A.; Huddart, R.; Khoo, V.; et al. Intrinsic Markers of Tumour Hypoxia and Angiogenesis in Localised Prostate Cancer and Outcome of Radical Treatment: A Retrospective Analysis of Two Randomised Radiotherapy Trials and One Surgical Cohort Study. Lancet Oncol. 2008, 9, 342–351. [Google Scholar] [CrossRef]
  39. Tran, M.G.B.; Bibby, B.A.S.; Yang, L.; Lo, F.; Warren, A.Y.; Shukla, D.; Osborne, M.; Hadfield, J.; Carroll, T.; Stark, R.; et al. Independence of HIF1a and Androgen Signaling Pathways in Prostate Cancer. BMC Cancer 2020, 20, 469. [Google Scholar] [CrossRef]
  40. Ribatti, D.; Tamma, R.; Annese, T. Epithelial-Mesenchymal Transition in Cancer: A Historical Overview. Transl. Oncol. 2020, 13, 100773. [Google Scholar] [CrossRef]
  41. Jolly, M.K.; Boareto, M.; Huang, B.; Jia, D.; Lu, M.; Ben-Jacob, E.; Onuchic, J.N.; Levine, H. Implications of the Hybrid Epithelial/Mesenchymal Phenotype in Metastasis. Front. Oncol. 2015, 5, 155. [Google Scholar] [CrossRef] [Green Version]
  42. Jolly, M.K.; Tripathi, S.C.; Jia, D.; Mooney, S.M.; Celiktas, M.; Hanash, S.M.; Mani, S.A.; Pienta, K.J.; Ben-Jacob, E.; Levine, H. Stability of the Hybrid Epithelial/Mesenchymal Phenotype. Oncotarget 2016, 7, 27067–27084. [Google Scholar] [CrossRef] [Green Version]
  43. Nieto, M.A.; Huang, R.Y.-J.; Jackson, R.A.; Thiery, J.P. EMT: 2016. Cell 2016, 166, 21–45. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Gupta, P.B.; Pastushenko, I.; Skibinski, A.; Blanpain, C.; Kuperwasser, C. Phenotypic Plasticity: Driver of Cancer Initiation, Progression, and Therapy Resistance. Cell Stem Cell 2019, 24, 65–78. [Google Scholar] [CrossRef] [PubMed]
  45. Yamashita, N.; Tokunaga, E.; Iimori, M.; Inoue, Y.; Tanaka, K.; Kitao, H.; Saeki, H.; Oki, E.; Maehara, Y. Epithelial Paradox: Clinical Significance of Coexpression of E-Cadherin and Vimentin With Regard to Invasion and Metastasis of Breast Cancer. Clin. Breast Cancer 2018, 18, e1003–e1009. [Google Scholar] [CrossRef] [PubMed]
  46. Stylianou, S.; Clarke, R.B.; Brennan, K. Aberrant Activation of Notch Signaling in Human Breast Cancer. Cancer Res. 2006, 66, 1517–1525. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Odero-Marah, V.; Hawsawi, O.; Henderson, V.; Sweeney, J. Epithelial-Mesenchymal Transition (EMT) and Prostate Cancer. Adv. Exp. Med. Biol. 2018, 1095, 101–110. [Google Scholar] [CrossRef] [PubMed]
  48. Cheaito, K.A.; Bahmad, H.F.; Hadadeh, O.; Saleh, E.; Dagher, C.; Hammoud, M.S.; Shahait, M.; Mrad, Z.A.; Nassif, S.; Tawil, A.; et al. EMT Markers in Locally-Advanced Prostate Cancer: Predicting Recurrence? Front. Oncol. 2019, 9, 131. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Papanikolaou, S.; Vourda, A.; Syggelos, S.; Gyftopoulos, K. Cell Plasticity and Prostate Cancer: The Role of Epithelial–Mesenchymal Transition in Tumor Progression, Invasion, Metastasis and Cancer Therapy Resistance. Cancers 2021, 13, 2795. [Google Scholar] [CrossRef]
  50. Wu, J.-S.; Jiang, J.; Chen, B.-J.; Wang, K.; Tang, Y.-L.; Liang, X.-H. Plasticity of Cancer Cell Invasion: Patterns and Mechanisms. Transl. Oncol. 2021, 14, 100899. [Google Scholar] [CrossRef]
  51. Allen, F.; Maillard, I. Therapeutic Targeting of Notch Signaling: From Cancer to Inflammatory Disorders. Front. Cell Dev. Biol. 2021, 9, 649205. [Google Scholar] [CrossRef]
  52. Anusewicz, D.; Orzechowska, M.; Bednarek, A.K. Notch Signaling Pathway in Cancer—Review with Bioinformatic Analysis. Cancers 2021, 13, 768. [Google Scholar] [CrossRef] [PubMed]
  53. Kwon, O.-J.; Zhang, L.; Wang, J.; Su, Q.; Feng, Q.; Zhang, X.H.F.; Mani, S.A.; Paulter, R.; Creighton, C.J.; Ittmann, M.M.; et al. Notch Promotes Tumor Metastasis in a Prostate-Specific Pten-Null Mouse Model. J. Clin. Investig. 2016, 126, 2626–2641. [Google Scholar] [CrossRef] [Green Version]
  54. Rahman, M.; Jackson, L.K.; Johnson, W.E.; Li, D.Y.; Bild, A.H.; Piccolo, S.R. Alternative Preprocessing of RNA-Sequencing Data in The Cancer Genome Atlas Leads to Improved Analysis Results. Bioinformatics 2015, 31, 3666–3672. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Paner, G.P.; Stadler, W.M.; Hansel, D.E.; Montironi, R.; Lin, D.W.; Amin, M.B. Updates in the Eighth Edition of the Tumor-Node-Metastasis Staging Classification for Urologic Cancers. Eur. Urol. 2018, 73, 560–569. [Google Scholar] [CrossRef] [PubMed]
  56. Kanehisa, M.; Goto, S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  57. Liberzon, A.; Birger, C.; Thorvaldsdottir, H.; Ghandi, M.; Mesirov, J.P.; Tamayo, P. The Molecular Signatures Database (MSigDB) Hallmark Gene Set Collection. Cell Syst. 2015, 1, 417–425. [Google Scholar] [CrossRef] [Green Version]
  58. Yevshin, I.; Sharipov, R.; Valeev, T.; Kel, A.; Kolpakov, F. GTRD: A Database of Transcription Factor Binding Sites Identified by ChIP-Seq Experiments. Nucleic Acids Res. 2017, 45, D61–D67. [Google Scholar] [CrossRef]
  59. Yevshin, I.; Sharipov, R.; Kolmykov, S.; Kondrakhin, Y.; Kolpakov, F. GTRD: A Database on Gene Transcription Regulation-2019 Update. Nucleic Acids Res. 2019, 47, D100–D105. [Google Scholar] [CrossRef] [Green Version]
  60. Ogluszka, M.; Orzechowska, M.; Jedroszka, D.; Witas, P.; Bednarek, A.K. Evaluate Cutpoints: Adaptable Continuous Data Distribution System for Determining Survival in Kaplan-Meier Estimator. Comput. Methods Programs Biomed. 2019, 177, 133–139. [Google Scholar] [CrossRef]
  61. Conway, J.R.; Lex, A.; Gehlenborg, N. UpSetR: An R Package for the Visualization of Intersecting Sets and Their Properties. Bioinformatics 2017, 33, 2938–2940. [Google Scholar] [CrossRef] [Green Version]
  62. Doncheva, N.T.; Morris, J.H.; Gorodkin, J.; Jensen, L.J. Cytoscape StringApp: Network Analysis and Visualization of Proteomics Data. J. Proteome Res. 2019, 18, 623–632. [Google Scholar] [CrossRef] [PubMed]
  63. Bader, G.D.; Hogue, C.W.V. An Automated Method for Finding Molecular Complexes in Large Protein Interaction Networks. BMC Bioinform. 2003, 4, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Rubicz, R.; Zhao, S.; April, C.; Wright, J.L.; Kolb, S.; Coleman, I.; Lin, D.W.; Nelson, P.S.; Ostrander, E.A.; Feng, Z.; et al. Expression of Cell Cycle-Regulated Genes and Prostate Cancer Prognosis in a Population-Based Cohort. Prostate 2015, 75, 1354–1362. [Google Scholar] [CrossRef] [PubMed]
Ijms 24 00164 g001 550
Figure 1. The profiles of core Notch signaling differentiate age subgroups of PC.
Figure 1. The profiles of core Notch signaling differentiate age subgroups of PC.
Ijms 24 00164 g001
Ijms 24 00164 g002 550
Figure 2. The specified Notch patterns of expression altered the DFS among subgroups of PC patients through the modulation of reshaping the cellular architecture. The unfavorable prognosis of signatures consisting of RBPJ, ADAM17, RBPJL, ATXN1 in early-stage PC below 55 years old, KAT2A, LFNG, HDAC2, SNW1, PSEN2 in late-stage PC below 55 years old, NOTCH3, APH1A, HDAC2, MFNG, DTX4 in early-stage PC among 60–70-year-olds, HES5, ADAM17, CREBBP, DVL3, ADAM10 in late-stage PC among 60–70-year-olds, and JAG1, PTCRA, HEYL, DTX3, RFNG in late-stage PC above 70 years old. Firstly, the DFS analysis was performed regarding the effects of the expression of particular Notch members through the maxstat algorithm, determining the optimal expression cutpoint stratifying the survival in the most significant manner. Next, by applying the UpSetR algorithm, we identified the specific subgroups of PC patients bearing a particular combination of the Notch core gene profiles associated with the DFS prognosis.
Figure 2. The specified Notch patterns of expression altered the DFS among subgroups of PC patients through the modulation of reshaping the cellular architecture. The unfavorable prognosis of signatures consisting of RBPJ, ADAM17, RBPJL, ATXN1 in early-stage PC below 55 years old, KAT2A, LFNG, HDAC2, SNW1, PSEN2 in late-stage PC below 55 years old, NOTCH3, APH1A, HDAC2, MFNG, DTX4 in early-stage PC among 60–70-year-olds, HES5, ADAM17, CREBBP, DVL3, ADAM10 in late-stage PC among 60–70-year-olds, and JAG1, PTCRA, HEYL, DTX3, RFNG in late-stage PC above 70 years old. Firstly, the DFS analysis was performed regarding the effects of the expression of particular Notch members through the maxstat algorithm, determining the optimal expression cutpoint stratifying the survival in the most significant manner. Next, by applying the UpSetR algorithm, we identified the specific subgroups of PC patients bearing a particular combination of the Notch core gene profiles associated with the DFS prognosis.
Ijms 24 00164 g002
Ijms 24 00164 g003 550
Figure 3. The functional network of the major molecules underlying invasive properties of PC with highly interconnected compartments (A). MCODE clustering indicated four of the most significant biological subgroups underlying age-associated course of PC, i.e., (B) the key participants of EMT program, (C) microtubule-related molecular motors, (D,E) parts of cell-ECM adhesion and interactions. We applied the clustering (ExpressCluster) to find common and unique expression profiles of the downstream Notch effectors among the age groups combined with the stage of the disease. Furthermore, by using the relevant clusters in terms of functional association with age-related PC course and progression, we constructed the biological networks of interactions, followed by MCODE identification of highly interconnected regions therein.
Figure 3. The functional network of the major molecules underlying invasive properties of PC with highly interconnected compartments (A). MCODE clustering indicated four of the most significant biological subgroups underlying age-associated course of PC, i.e., (B) the key participants of EMT program, (C) microtubule-related molecular motors, (D,E) parts of cell-ECM adhesion and interactions. We applied the clustering (ExpressCluster) to find common and unique expression profiles of the downstream Notch effectors among the age groups combined with the stage of the disease. Furthermore, by using the relevant clusters in terms of functional association with age-related PC course and progression, we constructed the biological networks of interactions, followed by MCODE identification of highly interconnected regions therein.
Ijms 24 00164 g003
Ijms 24 00164 g004 550
Figure 4. The functional network comprising the most essential invasion-related molecules shows expression profiles shared by (A) early-onset PC, (B) 60–70-year-old PC patients, and (C) elderly PC patients, reflecting the DFS prognosis according to the specific Notch signature. The adjacent Table 6 presents fold change (FC) of gene expression across different subgroups of patients.
Figure 4. The functional network comprising the most essential invasion-related molecules shows expression profiles shared by (A) early-onset PC, (B) 60–70-year-old PC patients, and (C) elderly PC patients, reflecting the DFS prognosis according to the specific Notch signature. The adjacent Table 6 presents fold change (FC) of gene expression across different subgroups of patients.
Ijms 24 00164 g004
Ijms 24 00164 g005 550
Figure 5. PC occurring in males below 55 years old showed intense activation of Notch signaling of Sender phenotype at early stage of disease, and, especially, at the late stage of disease through acquisition of hybrid Sender/Receiver phenotype (high expression of both Notch ligands and receptors). In turn, among males over 70 years old, Notch signaling has been diminished and turned into Receiver phenotype of PC cells at early stage of disease, whereas at the late-stage Notch has been inactivated. Downstream effects governed by Notch-specific TFs (HES/HEY) involved a change of cell morphology through reorganization of cellular architecture. The sedentary epithelial cells have apical/basal polarity, the cytoskeleton comprises the actin cortex, cytokeratin intermediate filaments, and microtubules with dynein and kinesin motors. The cell is anchored to the ECM by actin-linked laminin junctions, and the between-cells adhesion is asserted by E-cadherin bridges. Due to Notch signaling in younger males, PC cells undergo partial EMT triggered by ZEB2 concomitantly expressing epithelial (CDH1) and mesenchymal (CTNNB1, FN1) markers. As demonstrated by the elements in dark pink, consequently, the cells start to migrate, pushed forward by contractions via stress fibers and detachment from the matrix eased by the activity of MMP7. Finally, at the leading edge, the invading cells adhere by focal adhesion (involving LAMA1). Many of these features have been lost with disease progression, and tough gain in VIM has been noticed. The PC cells in older males demonstrated, in turn, inverse profiles to those of patients below 55 years old. Moreover, we observed a much-intensified activation of the EMT program through involvement of all specific TFs (SNAI2, TWIST1, ZEB2, HIF1A, SMAD3) at both stages of PC. Late-stage tumors restored the expression of CDH1, simultaneously maintaining CTNNB1 and VIM independent of PC stage. Tumors of both stages kept adhesive molecules like integrins and collagens. Importantly, high activity of microtubules and kinesin/dynein was kept.
Figure 5. PC occurring in males below 55 years old showed intense activation of Notch signaling of Sender phenotype at early stage of disease, and, especially, at the late stage of disease through acquisition of hybrid Sender/Receiver phenotype (high expression of both Notch ligands and receptors). In turn, among males over 70 years old, Notch signaling has been diminished and turned into Receiver phenotype of PC cells at early stage of disease, whereas at the late-stage Notch has been inactivated. Downstream effects governed by Notch-specific TFs (HES/HEY) involved a change of cell morphology through reorganization of cellular architecture. The sedentary epithelial cells have apical/basal polarity, the cytoskeleton comprises the actin cortex, cytokeratin intermediate filaments, and microtubules with dynein and kinesin motors. The cell is anchored to the ECM by actin-linked laminin junctions, and the between-cells adhesion is asserted by E-cadherin bridges. Due to Notch signaling in younger males, PC cells undergo partial EMT triggered by ZEB2 concomitantly expressing epithelial (CDH1) and mesenchymal (CTNNB1, FN1) markers. As demonstrated by the elements in dark pink, consequently, the cells start to migrate, pushed forward by contractions via stress fibers and detachment from the matrix eased by the activity of MMP7. Finally, at the leading edge, the invading cells adhere by focal adhesion (involving LAMA1). Many of these features have been lost with disease progression, and tough gain in VIM has been noticed. The PC cells in older males demonstrated, in turn, inverse profiles to those of patients below 55 years old. Moreover, we observed a much-intensified activation of the EMT program through involvement of all specific TFs (SNAI2, TWIST1, ZEB2, HIF1A, SMAD3) at both stages of PC. Late-stage tumors restored the expression of CDH1, simultaneously maintaining CTNNB1 and VIM independent of PC stage. Tumors of both stages kept adhesive molecules like integrins and collagens. Importantly, high activity of microtubules and kinesin/dynein was kept.
Ijms 24 00164 g005
Table
Table 1. A detailed comparison of the Notch expression between tumor stages according to the age of PC onset.
Table 1. A detailed comparison of the Notch expression between tumor stages according to the age of PC onset.
Age Group
≦5560–70>70
EarlyLateEarlyLateEarlyLate
UPcoreligands: DLL1, DLL3, JAG1, TFs: HEY1, HES1, RBPJ, modulators: PSENEN, PSEN2, LFNG, APH1A, NCSTN, ADAM10, PSEN1, MAML3ligands: DLL1, DLL4, JAG1, JAG2, receptors: NOTCH2, NOTCH3, NOTCH4, TFs: RBPJ, RBPJL, HEYL, PTCRA, HES1, modulators: APH1A, NCSTN, MAML2, NUMB, PSEN1, APH1B, PSEN2, MAML3, ADAM10, ADAM17, MAML1, PSENENligands: DLL1, DLL3, DLL4, JAG2, DLK1,
receptors: NOTCH3, NOTCH4,
TFs: HEY1, HES1, HEYL, HES4, HEY2, HES5
modulators: PSENEN, PSEN2, LFNG, RFNG, NUMBL		ligands: DLL4, DLL1,
JAG1, JAG2, DLL3, DLK1,
receptors: NOTCH1, NOTCH2, NOTCH3, NOTCH4,
TFs: RBPJ, RBPJL, HEYL, PTCRA, HES4, HEY1,
modulators: APH1A, NCSTN, MAML2, NUMB, PSEN1, APH1B,
PSEN2, MAML3, ADAM10, ADAM17,
MAML1,
RFNG,
MFNG		ligands: JAG1, DLL4,
receptors: NOTCH1, NOTCH2, NOTCH3,
TFs: RBPJ, PTCRA,
RBPJL, HEYL, HES4, HEY2, HES5,
modulators:
NCSTN, ADAM10, PSEN1, MAML3, NUMB, MAML2, APH1B, MAML1, ADAM17, MFNG, NUMBL 		ligands: DLK1, DLL3,
receptors:
TFs: HEY2,
HES4, HES5,
modulators: NUMBL, LFNG, PSENEN, RFNG
cross-talkepigenetic modulators: HDAC1, HDAC2, EP300,
regulators: CTBP1, CTBP2, KAT2A, SNW1, DTX3, CIR1, ATXN1L,
Wnt cross-talk: DVL3		epigenetic modulators: CREBBP, EP300, HDAC1, HDAC2,
regulators: ATXN1L, ATXN1, CIR1, DTX1, DTX3, DTX4, KAT2B, NCOR2, CTBP1,
Wnt cross-talk: DVL1, DVL2, DVL3 		regulators: CTBP1, KAT2A, SNW1, DTX3, CIR1, NCOR2, DTX4, DTX2,
Wnt cross-talk: DVL1, DVL3		epigenetic modulators: CREBBP, EP300,
HDAC1,
regulators: ATXN1L, ATXN1, CIR1, DTX4,
NCOR2, KAT2B, CTBP1, DTX3, KAT2A,
Wnt cross-talk: DVL3, DVL2, DVL1		epigenetic modulators: EP300, CREBBP,
regulators: ATXN1L, ATXN1L,
DTX3L, DTX1,
KAT2B, NCOR2, DTX4, DTX2,
Wnt cross-talk: DVL2, DVL1		epigenetic modulators: HDAC2,
regulators: SNW1, DTX3L,
CTBP2, DTX2, KAT2A
DOWNcoreligands: JAG2, DLL4, DLK1
receptors: NOTCH4, NOTCH1, NOTCH2, NOTCH3,
TFs: PTCRA, RBPJL, HEYL, HES4, HEY2, HES5,
modulators: RFNG, NUMB, MAML2, APH1B, MAML1, ADAM17, MFNG, NUMBL		ligands: DLK1, DLL3,
receptors: NOTCH1
TFs: HEY2, HES4, HES5, HEY1,
modulators: NUMBL, LFNG, RFNG, MFNG		ligands: JAG1,
receptors: NOTCH1, NOTCH2,
TFs: RBPJ, PTCRA, RBPJL,
modulators: APH1A, NCSTN, ADAM10, PSEN1, MAML3, NUMB, MAML2, APH1B, MAML1, ADAM17, MFNG		TFs: HES1, HEY2, HES5,
modulators: NUMBL, LFNG, PSENEN		ligands: DLL1, DLL3, JAG2, DLK1,
receptors: NOTCH4,
TFs: HEY1, HES1,
modulators: PSENEN, PSEN2, LFNG, APH1A, RFNG		ligands: DLL4, DLL1, JAG1, JAG2,
receptors: NOTCH2, NOTCH3, NOTCH4, NOTCH1,
TFs: RBPJ, RBPJL, HEYL, PTCRA, HES1, HEY1,
modulators: APH1A, NCSTN, MAML2, NUMB, PSEN1, APH1B, PSEN2, MAML3, ADAM10, ADAM17, MAML1, MFNG
cross-talkepigenetic modulators: CREBBP,
regulators: DTX3L, DTX1, KAT2B, NCOR2, DTX4, DTX2,
Wnt cross-talk: DVL2, DVL1, 		regulators: SNW1, DTX3L, CTBP2, DTX2, KAT2Aepigenetic modulators: HDAC1, HDAC2, EP300, CREBBP,
regulators: CTBP2, ATXN1L, ATXN1L, DTX3L, DTX1, KAT2B,
Wnt cross-talk: DVL2		epigenetic modulators: HDAC2,
regulators: DTX1, SNW1, DTX3L, CTBP2, DTX2		epigenetic modulators: HDAC1, HDAC2,
regulators: CTBP1, KAT2A, SNW1, DTX3, CIR1, CTBP2,
Wnt cross-talk: DVL3		epigenetic modulators: CREBBP, EP300, HDAC1,
regulators: ATXN1L, ATXN1, CIR1, DTX4, NCOR2, KAT2B, CTBP1, DTX3, DTX1,
Wnt cross-talk: DVL3, DVL2, DVL1
Table
Table 2. Corroboration of the primary findings with validation cohort (detailed lists of genes are available in Supplementary Table S2).
Table 2. Corroboration of the primary findings with validation cohort (detailed lists of genes are available in Supplementary Table S2).
Disease StageExpression PatternInAge Group
≤5560–7070>
General Trends
early/localUPTCGA only7825
Rubicz et al. only14198
Overlapping trend182010
DOWNTCGA only182310
Rubicz et al. only6521
Overlapping trend12611
late/regionalUPTCGA only19268
Rubicz et al. only8720
Overlapping trend21187
DOWNTCGA only10825
Rubicz et al. only15217
Overlapping trend6416
Shift in Expression between Stages
early -> late
local -> regional		UP -> DOWNTCGA only6828
Rubicz et al. only13187
Overlapping trend232
DOWN -> UPTCGA only20239
Rubicz et al. only10516
Overlapping trend524
Table
Table 3. The core participants of the Notch pathway differentiate DFS according to patients’ age and stage of the disease. The number represents the hazard ratio (HR) and the color indicates the trend–the blue color indicates that lower expression (i.e., below the determined cutpoint) was more favorable, and the red indicates that higher expression (i.e., above the determined cutpoint) was more favorable regarding DFS prognosis. Legend: 1 All patients were tumor-free; * p < 0.05, ** p < 0.01, *** p < 0.001.
Table 3. The core participants of the Notch pathway differentiate DFS according to patients’ age and stage of the disease. The number represents the hazard ratio (HR) and the color indicates the trend–the blue color indicates that lower expression (i.e., below the determined cutpoint) was more favorable, and the red indicates that higher expression (i.e., above the determined cutpoint) was more favorable regarding DFS prognosis. Legend: 1 All patients were tumor-free; * p < 0.05, ** p < 0.01, *** p < 0.001.
≤5560–7070>
General
HR		Early
HR		Late
HR		General
HR		Early
HR		Late
HR		General
HR		Early
HR 1		Late
HR
ADAM10 5.22 * 14.4 **0.373 *
ADAM173.16 *100> *4.46 * 0.455 *
APH1A 2.75 ***15.1 **1.9 *
APH1B0.29 * 0.069 *
ATXN1 100> * <0.001 ** <0.001 *
ATXN1L 0.14 *
CIR13.01 * 0.112 *
CREBBP 0.402 *
CTBP1 2.97 *
CTBP2 0.251 * <0.001 **
DLK1 2.15 * 2.26 *
DLL1 0.35 * 0.32 **
DLL3 0.482 *0.12 *
DLL4 0.069 *
DTX15.64 * <0.001 *
DTX2 0.347 *<0.001 **
DTX3 1.91 * 0.0874 * 0.107 *
DTX3L 2.85 **11.2 *0.515 *
DTX4 0.0873 *
DVL1
DVL2
DVL34.91 * 2.18 **9.73 *0.299 **9.54 *
EP300 <0.001 *
HDAC1 9.4 *0.398 *
HDAC2 3.19 * 12.7 **
HES1 0.451 *
HES4
HES5 0.339 **100> **
HEY10.296 * 0.157 *
HEY2
HEYL 0.284 *1.86 * 100> **
JAG1 0.079 ** 0.093 *
JAG23.03 * 6.12 * 1.94 *0.069 *
KAT2A 100> *3.58 * 100> *
KAT2B 0.382 **0.107 *
LFNG 6.77 ** <0.001 ***
MAML1
MAML2
MAML3 <0.001 *
MFNG 0.0539 ***
NCOR2
NCSTN 0.434 **0.0954 **
NOTCH1
NOTCH2 1.89 *
NOTCH3 2.37 **24.9 ***
NOTCH43.34 * 0.069 *
NUMB3.25 * 5.65 *
NUMBL
PSEN13.48 * 3.87 * 0.44 **
PSEN20.234 * 0.226 *1.91 *<0.001 * <0.001 * <0.001 **
PSENEN6.76 * 0.099 ** 0.251 *
PTCRA <0.001 * 12.3 **
RBPJ 100> * <0.001 *
RBPJL 100> *
RFNG 0.0874 * 0.107 *
SNW14.21 * 7.62 ***1.95 *14.1 **
Table
Table 4. Cluster analysis revealed molecular alterations affecting invasive properties of PC driven by the Notch pathway. Blue color and “↓” indicates downregulation of the expression, red color and “↑” indicates upregulation of the expression, and NA indicates not available expression profile of a particular gene. All the resulting clusters are available in the Supplementary File S3.
Table 4. Cluster analysis revealed molecular alterations affecting invasive properties of PC driven by the Notch pathway. Blue color and “↓” indicates downregulation of the expression, red color and “↑” indicates upregulation of the expression, and NA indicates not available expression profile of a particular gene. All the resulting clusters are available in the Supplementary File S3.
RNA-SeqMicroarrays
EarlyLateLocalRegional
≤5560–7070>≤5560–7070>≤5560–7070>≤5560–7070>
ANK2actin microfilaments↓↑↑↓↑↑↓↑↓↓↑↓
ANXA1↓↑↑↓↓↑↓↑↓↓↑↓
CALD1↑↓↑↓↓↑↑↑↑↓↓↓
CRB2↑↓↑↑↓↑↑↑↓↑↑↓
DST↑↓↓↓↑↑↑↑↓↓↓↑
FLNB↑↓↓↓↑↑↑↑↓↑↑↓
FN1↑↓↓↓↑↑↓↑↓↓↑↓
HIP1↑↑↓↑↓↓↓↓↓↑↑↑
MYO1D↑↓↓↓↑↑↑↑↓↑↑↓
MYO1F↑↑↓↑↑↑↓↓↑↓↓↑
MYO7A↓↓↓↓↓↑↑↓↓↑↓↑
CDH1adhesion↑↓↑↓↓↑↑↓↓↑↓↓
CDH11↓↓↓↓↑↑↓↓↑↓↓↑
CDH15↓↓↑↑↑↓↑↑↓↑↓↓
CDH22↑↑↓↓↓↓NA
CDH26↓↑↑↓↑↑↑↑↑↑↑↓
CDH4↓↓↑↑↑↓↓↓↓↑↑↓
CDHR2↑↓↑↑↓↓↓↑↓↓↓↑
CDHR5↑↓↓↑↓↓↑↑↓↑↑↓
CLEC16A↑↓↓↓↑↑↓↑↑↑↓↓
CTNNA2↓↑↑↑↑↓↑↑↑↓↑↓
CTNNB1↑↓↓↓↑↑↑↑↓↓↑↑
GJA1↓↑↓↓↑↑↓↓↑↓↓↑
ICAM1↑↓↓↓↓↑↓↑↓↓↑↑
ICAM5↑↑↓↑↑↓↓↑↓↓↑↑
PLEC↑↓↓↓↑↑↑↑↓↑↓↓
TJP3↑↓↓↓↓↑↓↓↑↓↓↑
VEGFAangiogenesis↑↓↑↓↓↑↓↓↑↓↓↑
VEGFC↑↓↓↑↓↓↓↑↓↑↑↑
BCL2apoptosis↑↓↓↓↑↑↓↓↑↓↓↑
CDK12cell cycle↑↓↓↓↓↑↑↑↓↑↑↓
CDK19↓↓↑↓↑↑↓↓↑↑↑↓
CDKN1A↑↑↓↑↓↓↓↑↓↓↑↓
RB1↑↓↓↓↑↑↓↓↑↓↓↑
TP63↓↑↑↓↑↑↓↓↑↓↓↑
LIN7Acell polarity↓↑↓↑↑↓↓↓↓↓↑↑
COCHECM↓↓↑↓↓↑↓↓↑↓↓↑
COL14A1↑↓↓↓↑↑↓↑↓↓↑↓
COL4A4↑↓↓↓↓↑↑↑↑↓↓↓
COL6A2↑↑↓↓↑↑↑↑↓↓↓↓
COL8A1↑↓↓↓↑↑↓↓↑↓↓↑
FBN1↑↓↓↓↑↑↑↑↓↑↓↑
FRAS1↓↑↓↓↑↑↓↑↓↑↑↑
GDF5↑↓↑↑↓↑↑↓↑↓↑↑
ITGA1↑↓↓↓↑↑↓↓↑↓↓↑
ITGA10↓↑↑↓↑↑↓↓↑↓↓↑
ITGA6↑↓↑↓↑↑↓↓↑↓↓↑
ITGA8↑↓↓↓↑↑↓↓↑↓↓↑
ITGAV↑↓↓↓↑↑↓↓↑↑↓↑
ITGB2↑↓↓↓↑↑↓↓↑↓↓↑
ITGB8↓↑↑↓↓↑↓↓↑↑↓↓
LAMA1↑↓↓↑↓↓↓↓↑↑↓↑
LAMC2↓↑↑↓↓↑↑↑↓↑↑↓
LAMC3↓↑↓↓↑↑↓↓↑↓↓↑
LGALS3↑↓↓↓↓↓↑↑↑↑↓↓
MATN1↑↓↑↑↓↑↓↓↓↓↑↑
MATN2↓↑↑↓↑↑↓↓↑↑↓↑
MMP7↑↓↓↓↑↑↓↑↑↓↓↓
NPNT↑↓↓↑↓↓↑↑↓↓↑↓
NTN4↓↑↑↓↑↑↓↓↑↑↓↓
TNN↓↓↓↑↓↑↑↑↓↑↑↓
TNXB↑↓↓↓↑↑NA
VTN↓↑↓↑↑↓NA
MMP2↓↑↓↓↑↑↓↓↓↓↓↑
EGFREGFR signaling↑↓↓↓↑↑↓↓↑↓↓↑
ERBB2↓↑↑↓↑↑↑↑↓↑↑↓
FCHO2endocytosis↑↓↑↓↑↑↓↓↑↓↓↑
ESR1estrogen↑↑↓↓↑↑↑↑↓↓↓↓
FGF1FGF signaling↑↑↓↓↑↑↓↑↓↓↓↑
FGF12↓↑↑↓↑↑↓↓↑↑↑↓
FGF7↑↓↓↑↓↓NA
FGF9↑↓↓↑↓↑↓↑↑↓↑↓
CREBBPepigenetic regulators↑↓↓↓↑↑↑↑↓↓↑↓
DNMT1↑↓↓↑↑↑↓↓↑↓↓↑
EP300↑↓↓↓↑↑↓↓↑↑↓↑
HDAC1↓↓↑↓↑↑↑↑↓↑↑↓
HDAC6↑↓↑↓↑↑↑↑↓↓↓↓
FLGintermediate filaments↑↓↑↓↓↑↑↑↓↑↑↓
VIM↓↑↓↓↑↑↑↑↓↓↓↑
MAP2microtubules↑↓↓↓↓↓↓↓↑↓↓↑
MAP7↓↓↑↓↑↑↓↓↑↑↓↑
MAPRE3↓↓↑↓↑↑↓↑↓↓↓↑
MZT1↓↓↑↓↑↑↓↓↑↑↑↓
TUBA4A↓↑↑↓↑↑↓↓↑↓↓↑
TUBGCP4↑↓↑↓↑↑NA
TUBGCP5↑↓↑↓↓↑NA
DYNC1H1molecular motors↑↓↓↓↑↑↓↑↓↓↑↓
KIF1B↑↓↓↓↑↑↑↑↓↑↓↑
KIF21A↑↓↓↓↑↑↑↑↓↑↑↓
KIF21B↑↓↓↓↑↑↓↓↑↓↓↑
KIF26A↑↑↓↓↑↑NA
KIF2A↑↓↓↓↑↑NA
KIF3A↑↓↓↓↑↑↓↑↓↓↑↓
KIF3B↑↓↓↓↓↓↓↑↓↑↑↑
KIF5A↓↓↓↓↓↑↓↓↓↓↑↑
KIFAP3↑↓↓↓↑↑↓↓↑↓↓↑
HES1Notch↓↓↑↑↑↑↓↓↑↑↓↓
HES4↓↑↓↑↑↓↓↑↓↓↑↓
HES5↓↑↑↑↑↑↓↑↑↑↓↓
HEY1↓↑↓↑↑↓↓↑↓↓↑↑
HEY2↓↓↓↑↑↓↑↑↓↑↑↓
HEYL↓↑↓↑↑↑↓↑↓↑↓↑
NOTCH1↑↑↓↓↑↑↓↓↑↓↓↑
ERGTF↓↓↑↓↑↑↑↑↓↑↑↓
GATA1↑↓↓↓↓↑↑↑↓↑↑↓
MYC↑↓↓↑↑↓↓↑↑↓↑↓
TCF4↑↓↓↓↑↑↓↓↑↓↓↑
TFAP2A↓↑↑↓↑↑↑↓↓↑↑↓
SNAI2TF/EMT↓↑↑↓↑↑↓↓↑↓↓↑
TWIST1↓↑↓↑↑↓↑↓↓↑↑↑
ZEB2↑↓↓↓↑↑↓↓↑↓↓↑
HIF1ATF/hypoxia↑↓↓↓↑↑↓↑↓↓↑↓
SMAD3TGFβ signaling/EMT↓↑↓↓↑↑↑↑↑↓↓↓
Table
Table 5. The Notch signature specific for each PC subgroup significantly differentiates DFS.
Table 5. The Notch signature specific for each PC subgroup significantly differentiates DFS.
Notch SignaturePrognosisnrmean *
Survival		SE(rmean) *Median Survival
(95% CI *)		HR*
(95% CI)		p
PatientsEvents
≦55 years oldearly stageRBPJ, RBPJL, ADAM17, ATXN1unfavorable2533734.3415.83524
(3524-NA)		2.23 × 10−09
(0-Inf)		0.23
favorable19025760NA0.08
(0.01–0.64)
late stageKAT2A, LFNG, HDAC2, SNW1, PSEN2unfavorable24101366.8185.41855
(791-NA)		8 × 10−10
(0-Inf)		0.0028
favorable2522122.4137.9NA
(2051-NA)		0.35
(0.18–0.69)
60–70 years oldearly stageNOTCH3, APH1A, HDAC2, MFNG, DTX4unfavorable2332819375.4NA3.5 × 10−10
(0-Inf)		0.021
favorable43027530NA2.23 × 10−09
(0-Inf)
late stageHES5, ADAM17, CREBBP, DVL3, ADAM10unfavorable3416952.9115.2717
(602-NA)		0.08
(0.01–0.64)		0.0016
favorable952519332101890
(1395-NA)		8 × 10−10
(0-Inf)
>70 years oldlate stageJAG1, PTCRA, HEYL, DTX3, RFNGunfavorable114467.294.4485
(244-NA)		0.35
(0.18–0.69)		0.013
favorable52104954.41049
(972-NA)		3.5 × 10−10
(0-Inf)
* rmean survival—restricted mean survival time; SE(rmean)—standard error of restricted mean survival time; 95% CI—95% confidence intervals; HR—hazard ratio.
Table
Table 6. The complementary data underlying networks in Figure 4 presenting expression profiles of the most essential invasion-related molecules that reflect the DFS prognosis according to the specific Notch signature. The values represent fold changes (FC) between unfavorable vs. favorable prognosis according to the specified Notch patterns.
Table 6. The complementary data underlying networks in Figure 4 presenting expression profiles of the most essential invasion-related molecules that reflect the DFS prognosis according to the specific Notch signature. The values represent fold changes (FC) between unfavorable vs. favorable prognosis according to the specified Notch patterns.
≦5560–70>70
EarlyLateEarlyLateLate
actin microfilamentACTA20.840.380.480.880.32
ANK21.170.560.860.490.38
ANXA10.760.740.910.890.61
CALD11.320.761.070.540.34
CRB21.071.251.230.770.92
DST1.820.810.800.600.87
FLNB1.571.061.390.901.03
FN11.420.721.110.470.50
HIP11.281.141.180.690.69
MYO1D1.280.861.250.850.76
MYO1F1.730.780.940.651.21
MYO7A1.630.661.320.840.44
adhesionCDH11.120.921.291.020.81
CDH111.350.861.020.441.25
CDH151.311.640.581.710.85
CDH21.291.070.920.550.56
CDH225.320.420.520.800.71
CDH261.611.190.901.430.90
CDH40.861.010.930.560.63
CDHR21.402.251.830.670.81
CDHR51.001.381.100.670.42
CLEC16A1.371.021.220.780.76
CTNNA21.241.000.821.400.39
CTNNB11.441.081.340.810.75
DSP1.361.171.500.771.18
GJA11.940.550.760.670.43
GSC0.871.211.300.931.24
ICAM11.300.981.470.790.80
ICAM51.001.211.550.790.70
MUC10.860.811.221.071.13
NID11.320.730.790.480.41
OCLN1.301.491.880.560.71
PLEC1.251.061.100.760.59
S100A41.330.770.950.851.52
SDC11.950.660.860.872.34
SERPINH11.130.911.071.000.66
TJP11.231.121.550.710.41
TJP31.030.691.061.020.95
cell polarityLIN7A1.281.231.220.850.61
ECMCOCH1.420.601.250.800.31
COL14A11.160.620.980.430.40
COL4A11.170.740.950.490.44
COL4A41.300.720.710.460.66
COL6A21.090.640.650.850.72
COL8A11.400.911.160.261.16
DDR22.550.881.400.280.47
FBN11.510.680.990.420.66
FRAS12.101.061.850.611.11
GDF50.750.600.530.940.35
ITGA11.640.641.030.420.42
ITGA100.881.391.310.760.91
ITGA50.890.580.840.720.52
ITGA61.650.831.380.750.78
ITGA81.501.061.020.360.89
ITGAV1.580.911.540.700.68
ITGB21.710.810.830.611.52
ITGB64.020.591.110.360.65
ITGB81.270.851.060.550.75
LAMA11.130.880.680.913.43
LAMA21.400.870.850.491.06
LAMA31.730.541.070.680.44
LAMA41.240.750.890.530.46
LAMA51.200.761.140.860.58
LAMC21.370.422.120.770.31
LAMC31.220.880.780.641.68
LGALS32.250.510.881.100.97
MATN11.080.901.060.890.59
MATN21.440.660.900.540.63
MMP21.060.570.870.590.72
MMP31.801.090.820.641.04
MMP71.310.880.790.580.99
MMP91.050.790.720.631.72
NPNT1.381.120.970.791.26
NTN42.220.791.090.680.82
TNN1.660.710.590.521.75
TNXB1.230.600.830.590.70
VTN0.650.980.880.770.63
intermediate filamentsFLG1.580.650.930.430.19
KRT180.890.860.801.440.94
KRT191.570.480.881.051.92
KRT52.900.470.820.840.77
KRT80.960.890.851.421.03
VIM1.140.790.990.740.69
microtubulesMAP21.311.201.230.750.74
MAP71.030.991.610.850.64
MAPRE30.951.021.170.970.46
MZT10.941.111.130.780.63
TUBA4A1.020.530.811.030.74
TUBGCP41.491.051.270.730.64
TUBGCP51.301.111.210.900.76
molecular motorsDYNC1H11.321.141.270.880.66
KIF1B1.621.021.580.850.60
KIF21A1.731.101.910.680.76
KIF21B1.411.171.280.580.82
KIF26A1.481.301.110.540.77
KIF2A1.521.061.530.620.66
KIF3A1.271.081.480.690.78
KIF3B1.291.181.250.870.78
KIF5A1.330.671.330.520.25
KIFAP31.271.001.180.830.76
TFFOXC20.780.810.750.440.38
LEF11.050.890.630.790.45
SMAD21.190.901.300.840.54
SMAD31.280.951.200.750.79
SNAI11.080.631.020.730.62
SNAI22.650.600.710.720.47
TCF31.101.001.290.890.79
TCF41.560.761.300.600.72
TWIST10.780.770.981.231.20
ZEB11.310.751.060.550.41
ZEB21.530.821.070.430.87
Table
Table 7. Detailed characteristics of the studied cohort.
Table 7. Detailed characteristics of the studied cohort.
All
(n = 397)		≤55 Years Old (n = 109)60–70 Years Old (n = 253)70> Years Old (n = 35)
Early
(n = 53)		Late
(n = 56)		Early
(n = 84)		Late
(n = 169)		Early
(n = 13)		Late
(n = 22)
Age63 (41–78)515364657272
-median
Race
-Asian2-1-1--
-Black or African American612-21-
-White1251518324974
-NA264373552117518
Tumor status250413763881110
-tumor free65312341-6
-with tumor8297184026
-NA
Vital status3915354831661322
-alive6-213--
-dead
Adjuvant radiation treatment19130234279512
-no32-5126--
-yes17423283164810
-NA
Treatment outcome of first course1241814304958
-complete remission/response23-4-14-4
-partial remission/response19--217--
-stable disease12-4-8--
-progressive disease21935345081810
-NA
Laterality3434352711461120
-bilateral17321101-
-left316112912
-right611-4--
-NA
Residual tumor254432871921010
-R0117722768211
-R141-11-1
-R212-237--
-RX1024211-
-NA
New tumor event19534254873510
-no53-9437-3
-yes1481922325989
-NA
Targeted molecular therapy 214384510
-no193306-21-2
-yes29-294164810
-NA17523
Stage
T [clinical/pathologic]
-T1a1/----1/---
-T1b1/----1/---
-T1c141/-30/-17/-45/-43/-3/-3/-
-T211/-2/-1/-2/-3/--3/-
-T2a44/126/47/-9/618/-1/23/-
-T2b45/113/49/-7/622/-1/13/-
-T2c37/1286/454/28/7118/-1/10-
-T3a28/127-7/282/-15/91-4/8
-T3b14/110-3/25-10/73-1/12
-T48-Feb-1/1-1/5--/2
-NA73/16/-7/-11/137/-7/-5/-
N [pathologic]
-N02823534681211014
-N160-16-37-7
-NA55186161131
M [clinical]
-M03624653801521021
-M1a1---1--
-M1c1-1----
-NA337241631
Zone of origin
-central zone41-12--
-overlapping/multiple zones987131055112
-peripheral zone1121513255234
-transition zone6131--1
-NA1772927476095
Gleason score
-6331418631
-72033529636277
-8494862632
-9109-17773-12
-103-1-2--
Primary pattern
-21-----1
-31494414502984
-420893433116511
-539-8124-6
Secondary pattern
-31272118305053
-418530264761813
-585212758-6
Biochemical recurrence indicator
-no2954341681141019
-yes4618433--
-NA5697122233
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Orzechowska, M.J.; Anusewicz, D.; Bednarek, A.K. Age- and Stage-Dependent Prostate Cancer Aggressiveness Associated with Differential Notch Signaling. Int. J. Mol. Sci. 2023, 24, 164. https://doi.org/10.3390/ijms24010164

AMA Style

Orzechowska MJ, Anusewicz D, Bednarek AK. Age- and Stage-Dependent Prostate Cancer Aggressiveness Associated with Differential Notch Signaling. International Journal of Molecular Sciences. 2023; 24(1):164. https://doi.org/10.3390/ijms24010164

Chicago/Turabian Style

Orzechowska, Magdalena Julita, Dorota Anusewicz, and Andrzej K. Bednarek. 2023. "Age- and Stage-Dependent Prostate Cancer Aggressiveness Associated with Differential Notch Signaling" International Journal of Molecular Sciences 24, no. 1: 164. https://doi.org/10.3390/ijms24010164

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop