- Research
- Open access
- Published:
TP53 missense mutation reveals gain-of-function properties in small-sized KRAS transformed pancreatic ductal adenocarcinoma
Journal of Translational Medicine volume 21, Article number: 872 (2023)
Abstract
Background
Although the molecular features of pancreatic ductal adenocarcinoma (PDAC) have been well described, the impact of detailed gene mutation subtypes on disease progression remained unclear. This study aimed to evaluate the impact of different TP53 mutation subtypes on clinical characteristics and outcomes of patients with PDAC.
Methods
We included 639 patients treated with PDAC in Ruijin Hospital affiliated to Shanghai Jiaotong University School of Medicine between Jan 2019 and Jun 2021. The genomic alterations of PDAC were analyzed, and the association of TP53 mutation subtypes and other core gene pathway alterations with patients’ clinical characteristics were evaluated by Chi-squared test, Kaplan-Meier method and Cox regression model.
Results
TP53 missense mutation was significantly associated with poor differentiation in KRASmut PDAC (50.7% vs. 36.1%, P = 0.001). In small-sized (≤ 2 cm) KRASmut tumors, significantly higher LNs involvement (54.8% vs. 23.5%, P = 0.010) and distal metastic rate (20.5% vs. 2.9%, P = 0.030) were observed in those with TP53 missense mutation instead of truncating mutation. Compared with TP53 truncating mutation, missense mutation was significantly associated with reduced DFS (6.6 [5.6–7.6] vs. 9.2 [5.2–13.3] months, HR 0.368 [0.200–0.677], P = 0.005) and OS (9.6 [8.0-11.1] vs. 18.3 [6.7–30.0] months, HR 0.457 [0.248–0.842], P = 0.012) in patients who failed to receive chemotherapy, while higher OS (24.2 [20.8–27.7] vs. 23.8 [19.0–28.5] months, HR 1.461 [1.005–2.124], P = 0.047) was observed in TP53missense cases after chemotherapy.
Conclusions
TP53 missense mutation was associated with poor tumor differentiation, and revealed gain-of-function properties in small-sized KRAS transformed PDAC. Nonetheless, it was not associated with insensitivity to chemotherapy, highlighting the neoadjuvant therapy before surgery as the potential optimized strategy for the treatment of a subset of patients.
Introduction
Pancreatic ductal adenocarcinoma (PDAC), the classical and the most common subtype of pancreatic cancer, represents increased incidence and mortality rates worldwide [1]. Surgery is the only potential curative option for a relative small proportion of patients detected early before local progression and distal metastasis. However, most patients will experience a recurrence and die within 5 years after surgery [2].
Among the four most common driver genes of PDAC (KRAS, CDKN2A, TP53 and SMAD4), KRAS mutations are the most common presented in up to 90% of PDAC. Currently, the most widely recognized model of PDAC carcinogenesis is the stepwise model [3], in which firstly, KRAS mutations lead to low-grade dysplastic pancreatic intraepithelial neoplasias (PanINs). Secondly, high-grade PanINs and invasive adenocarcinomas are mostly driven by TP53 and/or CDKN2A and/or SMAD4 mutations. Other low-frequency genetic alterations include the alterations of the COMPASS and SWI/SNF family of Trithorax genes which may also interlink with TP53 and cell cycle pathways, and the alterations of other genes related to DNA damage repair pathway, mismatch repair pathway, RNA processing pathway, PI3K-Akt pathway, WNT pathway, NOTCH pathway, Hedgehog pathway and DNA modification pathway etc., [4] as listed in Additional file 4: Table S1_The Detailed Pathways and Genes Related to PDAC Carcinogenesis and Included in the Genetic Analyses Panel.
TP53 mutations occur in 60–80% of PDAC [5]. The majority of TP53 mutations are missense mutations clustered in the regions encoding its central DNA-binding domain, while a smaller but considerable subset result in truncated TP53 proteins expression [4, 6]. It was noted early on that KRASLSL.G12D/+, Trp53R172H/+, PdxCre mice (KPC mice) commonly showed an accelerated progression of PDAC with metastases observed in around 80% of the animals, while reduced invasiveness and decreased metastasis were observed in mice bearing a heterozygous Trp53 conditional knockout allele [7,8,9]. This has been explained by the potential gain-of-function (GOF) properties of TP53 missense mutations compared with other TP53 loss-of-function (LOF) mutations [10]. Indeed, recent studies also reinforced that TP53 missense mutations heightened the transformed KRAS function and they cooperated to drive PDAC invasiveness and metastasis [11, 12].
However, the GOF hypothesis also remains controversial now and then in the context of the large diversity of TP53 mutations [13]. On one hand, there is still a lack of clinical evidences concerning the differences between TP53 GOF and LOF mutations. On the other, it remains to be studied to which extent the GOF properties of TP53 missense mutations extend, because not all the TP53 missense mutations are the same. A recent study compared certain well-referenced TP53 GOF mutations (R175H, R248W, R248Q, R249S, R273H, R273L and R282W) with the rest in advanced and non-resectable PDAC, and found that they were associated with worse survival [14]. However, these subtypes of TP53 mutations only account for about 20% of patients with TP53 mutations, and about 15% of all the patients, while TP53 missense mutations account for 40–45% of all the PDAC patients. So we found it also interesting to study the GOF properties of TP53 missense mutations, even though not all the TP53 missense mutations are GOF mutations. In this study, we aimed to use a clinical database to evaluate the underlying GOF properties of TP53 missense mutations in the Chinese Asian population with PDAC.
Methods
Patient cohort
A total of 639 patients diagnosed with classical PDAC and undergoing genetic analysis in the pancreatic center of Shanghai Ruijin Hospital were enrolled in this study. All of them are Chinese Asians. There were 513 patients undergoing curative surgery, and the biopsy of the primary tumor or the metastatic site was performed in the rest 126 patients between Jan 1st, 2019 and Jun 30th, 2021. The rate of patients undergoing curative surgery was high because the genetic analysis was not done in all the patients with locally advanced or metastatic PDAC, and those late-staged patients without genetic analysis were not included. Patients diagnosed pathologically with other subtypes of pancreatic cancer (intraductal papillary mucinous neoplasm or mucinous cystic neoplasm associated with invasive carcinoma, colloid carcinoma, signet ring cell carcinoma, adenosquamous carcinoma, adenocarcinoma with squamous differentiation, squamous carcinoma, undifferentiated carcinoma, acinar cell carcinoma, ampulla carcinoma etc.) were excluded.
The postoperative treatment mainly referred to adjuvant chemotherapy, the schemes of which included FOLFIRINOX (leucovorin/FOL, fluorouracil/F, irinotecan/IRIN, oxaliplatin/OX) and AG (nab-paclitaxel/A, gemcitabine/G) consistent with the schemes of neoadjuvant therapy in PDAC, and included other derived schemes such as FOLFIRINOX plus PARP inhibitor, monotherapy by gemcitabine etc. [15, 16]. Patients undergoing neoadjuvant therapy before surgery were excluded because firstly, the genetic analysis in our center was mostly done after the neoadjuvant therapy and the surgery, instead of before the neoadjuvant therapy. Secondly, the schemes of neoadjuvant therapy included platinum-based FOLFIRINOX which might induce DNA damage and somatic mutations, and we found it hard to distinguish whether these mutations were carcinogenesis-associated or due to neoadjuvant therapy. And thirdly, the analysis in this study was already complex, so we didn’t include patients after neoadjuvant therapy in this study.
The clinical data were collected before and after surgery or biopsy. The tumor size was determined by CT and pathology in resectable cases, and by CT in metastatic cases. The follow-up data of patients undergoing curative surgery were used for survival analysis, and the last follow-up time was April 30th, 2023. The study was approved by the ethics committee of Shanghai Ruijin Hospital and all patients signed informed consent. The flowchart of the study was shown in Additional file 5: Fig. S1.
Genetic analysis
The genetic analyses were all performed in the clinical laboratory of Shanghai Ruijin Hospital using next-generation sequencing (NGS) as previously described [17]. Briefly, genomic DNA of the pancreatic tumor was extracted from the curative resection specimen or the biopsy specimen. Genomic DNA of peripheral blood lymphocyte was extracted from blood for somatic mutation calling and for germline variant calling if the consent of the latter was obtained.
As shown in Additional file 4: Table S1_The Detailed Pathways and Genes Related to PDAC Carcinogenesis and Included in the Genetic Analyses Panel, the panel included KRAS and other genes related to the RAS-MAPK pathway (ALK, AXL, DDR2, ERBB1-4, FGFR1-4, KDR, KIT, MET, NTRK1-3, PDGFRA and PDGFRB, RET, SYK, CIC, HRAS, NRAS, NF1, ARAF, BRAF, RAF1, MAP2K1-2, MAP3K1, JAK1-3), TP53, genes related to the cell cycle pathway (CDKN2A, CDKN2B, CDKN1B, RB1), the TGFβ pathway (SMAD4, SMAD2, SMAD3, TGFβR1, TGFβR2), the COMPASS family of the Trithorax genes (KDM6A, KMT2A, KMT2C, KMT2D), the SWI/SNF family of the Trithorax genes (ARID1A, ARID1B, ARID2, CEBPA, PBRM1, SMARCA4, SMARCB1), the homologous recombination pathway (ARID1A which also belongs to the SWI/SNF family, ATM, BAP1, BARD1, BLM, BRCA1, BRCA2, BRIP1, CHEK2, FANCA, FANCC, NBN, PALB2, RAD50, RAD51C, RAD51D), the mismatch repair pathway (MLH1, MSH2, MSH6, PMS2), the RNA processing pathway (RBM10, SF3B1, U2AF1), the PI3K-Akt pathway (PTEN, TSC1, TSC2), the WNT pathway (APC, AXIN1, CTNNB1, FAT1, SOX9, RNF43), the NOTCH pathway (NOTCH1, NOTCH2, NOTCH3, NOTCH4), the hedgehog pathway (PTCH1, SUFU), and the DNA modification pathway (MUTYH, DNMT3A, DNMT3B, TET2).
The missense mutation was determined to be deleterious or benign by Polymorphism Phenotyping v2 (PolyPhen-2) system [18]. Frameshift insertion or deletion, nonsense mutation, splice-site mutation and initiation codon mutation were considered to be truncating mutation. Other types of mutation included in-frame insertion or deletion, deletion-insertion mutation, fusion, silent mutation and readthrough mutation, and were considered to be deleterious in this study. The annotations of mutations were listed in Additional file 1: Raw Data 1_Annotations of mutations, and the clinical and genetic information of all the patients were recorded in Additional file 2: Raw Data 2_Information of patients.
Definition of mutational subtypes of PDAC
Among the 96 single-base substitution (SBS) COSMIC mutational signatures (Version 3.3) based on the prevalence of the six possible base substitutions, three types of signatures predominate in PDAC which are those related to age (SBS 1/clock-like) in the most of cases, to homologous recombination deficiency (HRD) in 10% (SBS 3) and to mismatch repair deficiency (MMRD) in 1–2% (SBS 6, 15, 21, 26 and 44) [3, 19, 20]. In this study, patients were classified as the MMRD subtype if deleterious mutations of the four MMRD genes existed, and as the HRD subtype if one was not the MMRD subtype and that deleterious mutations of the HRD genes were found. The rest of patients were considered as the age-related subtype.
Statistical analysis
SPSS version 22.0 (IBM) was used for statistical analysis and GraphPad Prism 5 was used for plotting survival curves. Differences between groups were evaluated by the Chi-squared test or Fisher’s exact test according to standard protocols. Briefly, the chi-squared test is considered if the sample size is large when expected counts all exceed 5. The Fisher’s Exact test is used for small samples when the use of chi-squared test is not appropriate, and where cells in the table have expected counts that are less than 5, and/or cell counts are smaller than 20, and/or the column or row marginal values are extremely uneven. The Kaplan–Meier plots and the univariate log-rank test were used for visually demonstrating disease-free survival (DFS) and overall survival (OS). Because the univariate log-rank test was not reliable enough to evaluate the survival differences, we used multivariate Cox regression model instead to compare the survival outcomes between different groups of the Kaplan–Meier plots, after adjusting for comutations and clinical covariates. The results of the multivariate Cox analysis were presented by hazard ratio (HR). A two-sided P value < 0.05 was considered statistically significant.
Cell culture, transfection and quantitative real-time PCR (qRT-PCR)
Human pancreatic cancer cell lines PANC-1 and CFPAC-1 were obtained from the Cell Bank of the Chinese Academy of Sciences. PANC-1 cells were cultured in DMEM and CFPAC-1 cells were cultured in IMDM supplemented with 10% fetal bovine serum (FBS).
For in vitro experiments, TP53 siRNA (si-TP53) and NC siRNA (si-NC) were synthesized by Bioegene (Shanghai, China). The transfections were performed using Hilymax (Dojindo). Cells were collected 48 h post-transfection. RNA was extracted using TRIzol reagent (Invitrogen, USA). The HiScript III RT SuperMix (TOYOBO, Japan) was used for reverse transcription, and the AceQ Universal SYBR qPCR Master Mix (AG, China) was used to detect the RNA expression levels, which were normalized to β-actin. The divergent and convergent primers of different types of mutant TP53 were used for RT-PCR. The siRNA sequences and the primer sequences are listed in Additional file 4: Table S2_The siRNA Sequences and the Primer Sequences Used in the Study.
Cell migration assays
The Transwell migration assays were performed as previously described [21]. Briefly, transfected pancreatic cancer cells (5 × 104 cells/100 µL) were suspended in serum-free medium and plated in the top chambers. The lower chambers were filled with 700 µL of DMEM or IMDM supplemented with 10% FBS. After 24 h, the pancreatic cancer cells that migrated from the upper chamber were fixed with 1% crystal violet stain solution for 20 min at room temperature and migrated cells were counted manually, and the relative cell migration percentage was calculated as the ratio of the number of migrated cells to that of the negative control.
For the wound healing assays, wound areas were made using 200 µL pipette tips after 48 h of transfection, and this time point was defined as 0 h. Afterwards, the cells were cultured for 36 h in serum-free medium after plating, and the wound areas were observed and photographed.
Results
Mutational landscape of the patient cohort
A total of 639 patients were included in this study and KRAS mutation was present in 595 (93.1%) patients, as shown in Additional file 4: Table S3_Patient Cohort and KRAS Status. 473 KRASmut patients and 40 KRASWT patients underwent curative surgery, whose overall clinical information were shown in Additional file 4: Table S4_Clinical Information of the Patient Cohort.
As shown in Additional file 4: Table S5_Core Gene Pathway Alterations in KRASmut and KRASWT PDAC, the 595 KRASmut PDAC included 491 (82.5%) age-related subtype, 92 (15.5%) HRD subtype and 12 (2.0%) MMRD subtype, while the 44 KRASWT PDAC included 26 (59.1%) age-related subtype, 17 (38.6%) HRD subtype and 1 (2.3%) MMRD subtype. Compared with KRASmut PDAC, KRASWT PDAC were more likely to be the HRD subtype (38.6% vs. 15.5%, P < 0.001), showed lower alteration frequencies in the three traditional pathways (TP53, cell cycle and TGFβ), comparable frequencies in Trithorax genes, and higher frequencies in NOTCH pathway, Hedgehog pathway and DNA modification pathway.
Core gene pathway alterations and association with tumor differentiation in KRASmut PDAC
The tumor differentiation status of each sample was judged by two independent pathologists from the Department of Pathology of Shanghai Ruijin Hospital with official pathology reports. The differentiation status of PDAC was according to defined WHO criteria, including the presence of tubular structures or solid growth, the presence of mucin, nuclear polymorphism and number of mitoses. Specifically, tumors with more than 30% of areas showing features of poor differentiation were categorized as poor differentiated PDAC. Examples of the histological features of well, moderate and poor differentiated PDAC are shown in Additional file 5: Fig. S2.
As shown in Table 1, among TP53 and other eleven pathways analyzed, TP53 was the only one related to poor tumor differentiation in KRASmut PDAC. Notably, the rate of TP53 missense mutation in poor differentiated PDAC was 50.7% (vs. 36.1% in moderate and well differentiated PDAC, P = 0.001), while the frequency of TP53 truncating mutation was comparable between the two groups (18.7% vs. 20.4%, P = 0.612). Same results were obtained in patients with age-related PDAC, where the rate of TP53 missense mutation was 50.9% in poor differentiated PDAC (vs. 40.0% in moderate and well differentiated PDAC, P = 0.021), as shown in Additional file 4: Table S6_Core Gene Pathway Alterations and Association with Tumor Differentiation in age-related KRASmut PDAC.
We also conducted multivariate Logistic regression to validate the influence of core gene pathway alterations on tumor differentiation in KRASmut PDAC. As shown in Additional file 4: Table S7_Multivariate Logistic Regression of Core Gene Pathway Alterations for Tumor Differentiation in KRASmut PDAC, TP53 missense mutation was the only alteration related to poor tumor differentiation in multivariate analysis (OR 1.848 [1.290–2.647], P = 0.001).
Core gene pathway alterations and association with lymph nodes involvement and distal metastasis according to tumor size in KRAS mut PDAC
At a median number of 13 lymph nodes examined per patient (range, 5–55), lymph nodes (LNs) involvement in resectable KRASmut PDAC only negatively correlated with WNT pathway alteration (LNs positive rate: 37.2% in WNTmut cases vs. 54.0% in WNTWT cases, P = 0.036) and Hedgehog pathway alteration (LNs positive rate: 16.7% in Hedgehogmut cases vs. 53.4% in HedgehogWT cases, P = 0.012), as shown in Table 2 and Additional file 4: Table S8_Other Gene Pathway Alterations and Association with Lymph Nodes Involvement. However, after stratification according to tumor size, we found that in the small-sized group (≤ 2 cm), the LNs involvement rate of patients with TP53 missense mutation remained as high as 54.8% compared with TP53WT cases (vs. 23.5%, P = 0.010). In cases with small-sized tumors, the LNs involvement rate of patients with TP53 truncating mutations (27.8% vs. 23.5% in TP53WT cases, P = 0.736), cell cycle pathway (50.0% vs. 34.7% in cell cycleWT cases, P = 0.310) and TGFβ pathway (41.7% vs. 36.1% in TGFβWT cases, P = 0.712) alterations was only slightly higher than relevant wild-type cases, without significant differences.
Interestingly, as shown in Table 2, the LNs involvement rate of large-sized (> 3 cm) resectable tumors with TP53 missense mutation was lower than those without TP53 mutation (LNs positive rate: 51.2% in TP53missense cases vs. 62.7% in TP53WT cases, P = 0.169), though without significant difference. It was also found that in resectable cases, the rate of TP53 missense mutation was higher in large-sized (> 3 cm) tumors (48.3% vs. 36.9% in small-sized tumors, P = 0.083), as shown in Additional file 4: Table S9_Frequency of Core Gene Pathway Alterations according to Tumor Size in KRASmut PDAC. We assumed that firstly, those with TP53 missense mutation that didn’t metastasize in the early stage were tumors with specifically lower invasiveness, therefore had lower LNs involvement rate. Secondly, if early metastasis hadn’t occurred, KRASmut PDAC with TP53 missense mutation tended to grow faster. Certainly, these two hypotheses need further validation.
With regard to distal metastasis, the metastatic rate of patients with TP53 mutations (20.4% vs. 17.5%, P = 0.389), cell cycle pathway (20.9% vs. 19.0%, P = 0.642) and Trithorax genes (22.0% vs. 18.6%, P = 0.384) alterations was only slightly higher than relevant wild-type cases, without significant differences, as shown in Table 3 and Additional file 4: Table S10_Other Gene Pathway Alterations and Association with Distal Metastasis. TGFβ pathway alteration was the only potential factor significantly associated with distal metastasis (metastatic rate: 26.8% in TGFβmut cases vs. 17.1% in TGFβWT cases, P = 0.011). The impact of TGFβ pathway alteration on distal metastasis seemed widely reported although its detailed characteristics remained controversial [22, 23]. After stratification according to tumor size, patients with TP53 missense mutation instead of truncating mutation showed significant higher metastatic rate than those without TP53 mutation (metastatic rate: 20.5% in KRASmutTP53missense PDAC vs. 2.9% in KRASmutTP53WT PDAC, P = 0.030), indicating that TP53 missense mutation probably led to early invasiveness and metastasis in KRASmut PDAC patients.
Similar results were obtained in age-related subtype of KRASmut PDAC. TP53 missense mutation in small-sized age-related cases was significantly associated with higher LNs involvement rate (LNs positive rate 53.6% in TP53missense cases vs. 25.8% in TP53WT cases, P = 0.029) and distal metastasis (metastatic rate: 17.6% in TP53missense cases vs. 0% in TP53WT cases, P = 0.025), as shown in Additional file 4: Table S11_Core Gene Pathway Alterations and Association with Lymph Nodes Involvement in Age-related PDAC, Additional file 4: Table S12_Frequency of Core Gene Pathway Alterations according to Tumor Size in Age-related KRASmut PDAC, and Additional file 4: Table S13_Core Gene Pathway Alterations and Association with Distal Metastasis in Age-related PDAC.
Core gene pathway alterations and association with patients’ survival in resectable KRASmut cases
The follow-up data of 458 patients with KRASmut PDAC and 38 patients with KRASWT PDAC who underwent curative surgery were used for survival analysis. At a median follow-up time of 22.4 months (range, 1.6–52.0 months), reduced DFS (12.8 vs. 16.1 months, HR 1.262 [0.829–1.922], P = 0.278) and OS (22.5 vs. 36.9 months, HR 1.545 [0.956–2.498], P = 0.076) were found in KRASmut cases compared with KRASWT cases, as shown in Fig. 1A and B.
As shown in Table 4, the multivariate survival analysis was conducted including 5 mutated pathways, sex, age, vascular invasion, tumor differentiation, tumor size, LNs involvement, resection margin status and adjuvant chemotherapy. It was demonstrated that in patients with KRASmut PDAC, KRASG12D mutation was associated with reduced DFS (11.8 vs. 14.0 months, P = 0.075) and OS (21.1 vs. 23.3 months, P = 0.139) compared with mutation of other KRAS codons without significant difference. TP53 mutation was associated with reduced DFS (10.6 months in TP53mut cases vs. 17.3 months in TP53WT cases, P = 0.005) and OS (21.0 months in TP53mut cases vs. 27.0 months in TP53WT cases, P = 0.047), and alterations of cell cycle pathway was also associated with reduced DFS (10.6 months in cell cyclemut cases vs. 13.9 months in cell cycleWT cases, P = 0.017) and OS (19.2 months in cell cyclemut cases vs. 23.3 months in cell cycleWT cases, P = 0.029) with significant difference. Alterations of TGFβ pathway was associated with reduced DFS and OS without significant difference, and that DFS and OS were comparable between patients with and without Trithorax genes alterations. The Kaplan–Meier curves for patients with KRASG12D mutation and with several pathway mutations were shown in Fig. 1C, D and E F, with HR calculated by Cox regression model after adjusting for comutations and clinical covariates.
The impact of TP53 missense mutation on patients’ survival
DFS or OS were comparable in patients with TP53 missense and truncating mutation (Fig. 2A and B). In patients with small-sized KRASmut PDAC, TP53 missense mutation was associated with reduced DFS (13.5 months in TP53missense cases vs. 21.8 months in TP53truncating cases, P = 0.046 by univariate log-rank test, HR 0.895 [0.311–2.576], P = 0.837 by Cox regression analysis) instead of OS (25.7 months in TP53missense cases vs. 28.4 months in TP53missense cases, HR 0.910 [0.300–2.761], P = 0.868), suggesting the potential early micrometastases in these cases (Fig. 2C and D).
However, the impact of TP53 mutation subtypes on patients’ survival remained controversial. Although TP53 missense mutation induced potential early metastases, TP53 truncating mutation was reported to more negatively correlated with patients’ survival than TP53 missense mutation, [24, 25] so we then evaluated whether this was due to tumor’s sensitivity to chemotherapy. Interestingly, TP53 missense mutation was significantly associated with reduced DFS (6.6 months in TP53missense cases vs. 9.2 months in TP53truncating cases, HR 0.368 [0.200–0.677], P = 0.005) and reduced OS (9.6 months in TP53missense cases vs. 18.3 months in TP53truncating cases, HR 0.457 [0.248–0.842], P = 0.012) in patients who failed to receive adjuvant chemotherapy (Fig. 2E, F). In patients who received chemotherapy, while DFS were comparable between the two groups (12.9 months in TP53missense cases vs. 13.2 months in TP53truncating cases, HR 1.231 [0.860–1.761], P = 0.257), TP53 missense mutation was even associated with slightly extended OS compared with tumors with TP53 truncating mutation (24.2 months in TP53missense cases vs. 23.8 months in TP53truncating cases, HR 1.461 [1.005–2.124], P = 0.047), as shown in Fig. 2G, H. The decreased survival differences between the missense and truncating group and the increased differences between the wild-type and truncating group in patients who received adjuvant chemotherapy induced us to hypothesize that tumors with TP53 missense mutation were not associated with insensitivity to chemotherapy, compared with tumors with TP53 truncating mutation.
The impact of the Trithorax gene alterations on patients’ survival according to TP53 mutation status
The Trithorax gene encodes a large family of proteins which serve as active epigenetic regulators counteracting the repressive gene expression programmes guided by the the Polycomb group of proteins [26]. The transcriptional repression of the locus INK4A/ARF downstream of TP53-Rb signaling axis is a potential target of the Polycomb group of proteins, thus alterations of the Trithorax genes might result in the repression of TP53 function, [27] as illustrated in Additional file 5:Â Fig. S3.
As the Trithorax genes were reported to be altered in less frequent but up to 10% of PDAC and that some of its components like KDM6A and ARID1A might influence the development, differentiation and metastasis of PDAC, [4, 28,29,30] we analyzed the relationship between the Trithorax genes alterations and the survival of patients with different TP53 mutation status. As shown in Fig. 3, although the survival outcomes between patients with and without Trithorax genes alterations were comparable (Fig. 3A and B), Trithorax genes alterations contributed to reduced DFS (12.9 months in Trxmut cases vs. 18.6 months in TrxWT cases, HR 1.580 [0.976–2.559], P = 0.063) and OS (24.0 months in Trxmut cases vs. 31.8 months in TrxWT cases, HR 1.254 [0.757–2.078], P = 0.379) in KRASmut, TP53WT PDAC patients (Fig. 3C and D), higher DFS (16.7 months in Trxmut cases vs. 11.6 months in TrxWT cases, HR 0.512 [0.267–0.982], P = 0.044) and OS (25.1 months in Trxmut cases vs. 19.9 months in TrxWT cases, HR 0.606 [0.309–1.190], P = 0.146) in KRASmut, TP53truncating PDAC patients (Fig. 3E, F), and slightly higher DFS (11.1 months in Trxmut cases vs. 9.5 months in TrxWT cases, HR 0.639 [0.426–0.956], P = 0.030) and OS (22.0 months in Trxmut cases vs. 20.6 months in TrxWT cases, HR 0.608 [0.385–0.960], P = 0.033) in KRASmut, TP53missense PDAC patients (Fig. 3G, H), indicating the different properties between TP53 missense and truncating mutations.
Validation by Independent TCGA cohort and functional experiments
We managed to validate our conclusions in two ways. Firstly, we used a cohort from the TCGA database including 126 pancreatic cancer patients to validate the association of TP53 mutation status with tumor differentiation. Secondly, we conducted functional experiments in vitro to validate the invasive properties of TP53 missense mutations.
We collected a cohort from the TCGA database with 126 pancreatic cancer patients who had recorded clinical information and SNV analysis, as recorded in Additional file 3: Raw Data 3_Information of the TCGA cohort. There are 93 KRASmut cases and 33 KRASWT cases. As shown in Additional file 4: Table S14_Core Gene Alterations and Association with Tumor Differentiation in KRASmut PDAC in the TCGA cohort, TP53 missense mutation was presented in 0 (0%) well differentiated tumors, while the rate was 54.0% in moderate differentiated tumors (P = 0.001) and 46.7% in poor differentiated tumors (P = 0.007). On the contrary, TP53 truncating mutation was presented in 72.7% well differentiated tumors vs. 20.0% in moderate differentiated tumors (P = 0.001) and 33.3% in poor differentiated tumors (P = 0.036), indicating TP53 missense mutation instead of truncating mutation was associated with poor tumor differentiation.
Notably, tumors with more than 30% of areas showing features of poor differentiation were categorized as poor differentiated PDAC in our study. As shown in Additional file 4: Table S14_Core Gene Alterations and Association with Tumor Differentiation in KRASmut PDAC in the TCGA cohort, there seemed to be more moderate differentiated tumors in the TCGA cohort, probably because that PDAC was defined as moderate differentiated in our hospital only when there was less than 30% of areas showing features of poor differentiation.
Since the volume of the TCGA cohort was not large enough to validate the invasive properties of TP53 missense mutations, we next used two pancreatic cancer cell lines both with KRAS and TP53 missense mutation to perform functional experiments. One cell line was PANC-1 carrying a heterozygous KRASG12D allele (c.35G > A) and a heterozygous TP53R273H allele (c.818G > A). The other cell line was CFPAC-1 carrying a heterozygous KRASG12V allele (c.35G > T) and a heterozygous TP53C242R allele (c.724T > C). A non-specific si-TP53 was transfected into the two cell lines to knock out the TP53 missense mutants, mimicking TP53 truncating mutations. As shown in Fig. 4, the Transwell migration assays (Fig. 4A) and the wound healing assays (Fig. 4B) showed that TP53 knock-out was associated with lower invasiveness compared with TP53 missense mutations. Using specific primers, qRT-PCR was performed to validate the knock-out of TP53 missense allele by non-specific si-TP53 (Fig. 4C).
Discussion
As shown inAdditional file 5:Â Fig. S4, the genomic alterations of PDAC involve different key pathways, which may also interplay one with another. It remained to be an intriguing question to which extent KRASmut and KRASWT PDAC varies one from the other. It was believed that RAS-MAPK pathway activation was still a crucial molecular driver of PDAC on condition that KRAS itself was not mutated [31]. Similarly in our study, 14 out of 44 (31.8%) KRAS WT PDAC bore a deleterious alteration of RAS-MAKP pathway other than KRAS. Alterations of other pathways also potentially contributed to the development of KRASWT PDAC, such as HRD pathway, NOTCH pathway, Hedgehog pathway and DNA modification pathway.
As for 595 KRASmut PDAC, 92 (15.5%) cases were classified as the HRD subtype. Some trivial facts reported by Golan T et al. were also found in our study [32]. Firstly, in KRASmut cases, TP53 was more likely to be intact in HRD subtype than in age-related subtype (mutation rate 52.2% vs. 66.6%, Additional file 4: Table S2). Secondly, TP53 mutation was more common in BRCA1 compared with BRCA2 inactivated cases (87.5% vs. 60.0%). Unlike the HRD subtype whose classification remained various and controversial, the MMRD subtype involved four genes and accounted for 2.0% of cases in our cohort, which was exactly the ratio reported by other groups [19, 33].
It was well established that PDAC could be classified into classical and basal-like/squamous subtypes according to the tumor’s transcriptome, while the latter tended to be more aggressive with increased metastatic potential [31, 34, 35]. Although not completely exclusive, the broadly defined TP53 mutation was to some extent associated with the basal-like subtype and tumor invasiveness [36, 37]. However, the outcomes of patients were usually more associated with tumors’ clinical and pathologic rather than molecular features [38]. In this study, we found that TP53 missense instead of truncating mutation was associated with poor differentiated PDAC, which was an important clinical prognostic factor. CDKN2A mutations have also been shown to promote undifferentiated PDAC in mouse models [34, 39]. Both in our cohort and in the TCGA cohort, however, we found that the alteration of the cell cycle pathway or CDKN2A itself was a less important factor than TP53 missense mutation to promote poor tumor differentiation, as shown in Table 1, Additional file 4: Table S7_Multivariate Logistic Regression of Core Gene Pathway Alterations for Tumor Differentiation in KRASmut PDAC, and Additional file 4: Table S14_Core Gene Alterations and Association with Tumor Differentiation in KRASmut PDAC in the TCGA cohort.
The cross-talk of different core gene pathways and its role in clinical outcomes, such as tumor differentiation, are complex in PDAC. For instance, the cell cycle pathway is also downstream of TP53 and Trithorax genes may interlink with TP53 and cell cycle pathways, as listed in Additional file 4: Table S1_The Detailed Pathways and Genes Related to PDAC Carcinogenesis and Included in the Genetic Analyses Panel. Form both univariate and multivariate analysis in our study, TP53 missense mutation was the only alteration significantly related to poor tumor differentiation. Interestingly, however, the two pathways which may interlink with TP53 were also related to poor tumor differentiation (cell cycle: 20.4% in poor differentiated PDAC vs. 16.2% in others; Trithorax: 24.6% in poor differentiated PDAC vs. 20.0% in others) yet without significant differences, as shown in Table 1 and Additional file 4: Table S7_Multivariate Logistic Regression of Core Gene Pathway Alterations for Tumor Differentiation in KRASmut PDAC. This might probably be explained by the cross-talk between TP53 and the two pathways. And for sure, the clinical outcomes of PDAC are results of the comprehensive interaction of different pathways. To be noted, potential alterations of GATA6 and MYC were not included in this study, because it is the copy number variation instead of the genetic mutation of these two genes that might influence the PDAC development and progression [31].
The GOF hypothesis involved the interaction between TP53, which mutated in its DNA binding domain, and other transcriptional regulators inducing gene expression modifications [40, 41]. It was recently further confirmed that compared with truncating mutations, TP53 missense mutations presented significantly greater cis-effect with higher TP53 protein expression and TP53-S315 phosphosite expression, while the trans-effects were similar between the two groups [42]. To our knowledge, we firstly described that under curative surgery along with adequate LNs resection and pathological examination, TP53 missense mutation was associated with higher LNs involvement rate in small-sized KRAS transformed PDAC. The distal metastatic rate was also significantly higher in small-sized tumors with TP53 missense mutation, which was in accordance with what was observed in genetically engineered mouse models [7,8,9]. Although it remains an open question whether and how LN metastasis plays an active role in shaping distant metastasis, [43] both of them indicate tumors’ invasiveness and the GOF properties of TP53 missense mutation. However, it was also true that when the tumor got bigger, the GOF properties of TP53 missense mutation became less obvious.
While the age-related subtype predominates in PDAC, the HRD subtype and the MMRD subtype are, to some extent, special PDAC subtypes [3]. Talia Golan et al. reported that the survival outcomes were comparable between the HRD subtype group and non-HRD group in resectable PDAC. However, the HRD status was predictive of platinum response and superior survival in advanced PDAC [32]. As for the MMRD subtype, Robert C Grant et al. found that MMRD-PDACs presented higher tumour mutational burden, were less likely to have mutations in KRAS and SMAD4, were more likely to have basal-like transcriptional programmes yet had longer OS after surgery [33]. As a result, the clinical manifestations of the MMRD subtype are different from the majority of PDAC, where basal-like tumors tend to be more aggressive with worse survival outcomes, as mentioned early in the discussion. Considering the potential differences among the three mutational subtypes of PDAC, we also investigated in this study whether same results were obtained in age-related PDAC to shed more solid evidence. Indeed, it was found that TP53 missense mutation was associated with poor differentiation in age-related KRASmut PDAC (Additional file 4: Table S6_Core Gene Pathway Alterations and Association with Tumor Differentiation in age-related KRASmut PDAC). And in small-sized (≤ 2 cm) age-related KRASmut tumors, TP53 missense mutation was associated with higher LNs involvement (Additional file 4: Table S11_Core Gene Pathway Alterations and Association with Lymph Nodes Involvement in Age-related PDAC) and distal metastic rate (Additional file 4: Table S13_Core Gene Pathway Alterations and Association with Distal Metastasis in Age-related PDAC).
The impact of the four driver genes’ mutation on patients’ outcomes have been studied since a decade ago [44, 45]. Qian ZR et al. reported the significant association of KRAS, especially KRASG12D mutation, and CDKN2A loss with DFS and OS, and of TP53 mutation with DFS only in resectable PDAC [24]. In this study, we evaluated the impact of genomic features on patients’ outcomes according to gene pathways, and found that in resectable KRASmut PDAC cases, KRASG12D and mutations of TP53, cell cycle and TGFβ pathways were all associated with poor survival, while the impact of TP53 and cell cycle pathway mutations turned out to be statistically significant after the multivariate analysis.
McIntyre CA et al. reported that TP53 truncating mutation predicted worse OS than TP53 missense mutation in resectable cases [25]. And a recent study showed that certain well-referenced TP53 GOF mutations (R175H, R248W, R248Q, R249S, R273H, R273L and R282W) were associated with worse survival in advanced and non-resectable PDAC [14]. In our study, we analyzed the survival difference between patients with missense mutations and truncating mutations. Instead of focusing on well-referenced TP53 GOF mutations, we mainly studied the GOF properties of all the patients with TP53 missense mutations. Firstly, the DFS and OS between the TP53 missense and truncating mutation group were comparable, except that TP53 missense mutation was associated with reduced DFS in small-sized PDAC, implying the potential early micrometastases in these cases. Secondly, TP53 missense mutation was associated with reduced DFS and OS than truncating mutation in patients who failed to receive adjuvant chemotherapy, probably due to its GOF properties. Thirdly, the impact of Trithorax genes alterations on patients’ outcomes varied one from another according to the TP53 mutation status. All these facts served as clinical evidences indicating the different properties between TP53 missense and truncating mutations, which might offer help to the potential targeted treatment against TP53 mutations in the future.
It was widely acknowledged that no matter what the mutation subtype is, TP53 mutations may lead to chemotherapy resistance compared with TP53WT cases, because the chemotherapy seems to be highly dependent on the apoptosis-inducing ability of TP53 [10]. Nonetheless, we still found in this study that the adjuvant chemotherapy extended DFS and OS in TP53missense cases (DFS from 6.6 to 12.9 months, OS from 9.6 to 24.2 months) and in TP53truncating cases (DFS from 9.2 to 13.2 months, OS from 18.3 to 23.8 months), indicating that the adjuvant chemotherapy may provide better survival both in TP53mut and TP53WT cases. Although after multivariate Cox analysis, TP53 missense mutation was an independent factor associated with better OS compared with TP53 truncating mutation in patients receiving adjuvant chemotherapy (24.2 months in TP53missense cases vs. 23.8 months in TP53truncating cases, HR 1.461 [1.005–2.124], P = 0.047), the OS of the TP53 missense group was only 0.4 months more than the truncating group. So we neither concluded that the chemotherapy would give better survival rate scenario in TP53missense group than in TP53truncating group, nor hypothesized that TP53 missense mutation was more sensitive to chemotherapy compared with TP53 truncating mutation. Instead, we found it more rational that TP53 missense mutation was not associated with insensitivity to chemotherapy compared with TP53 truncating mutation.
It was recently reported that the neoadjuvant therapy substantially ameliorated the 5-year OS rate compared to upfront surgery (20.5% vs. 6.5%) [46]. According to this study, PDAC with TP53 missense mutations tend to generate local LNs involvement and distal micrometastases in the early stage and are no less sensitive to chemotherapy. Besides, the majority of TP53 mutations in PDAC are the missense subtype [4, 10]. This might probably explain the significant higher 5-year OS rate in patients receiving neoadjuvant therapy, highlighting the chemotherapy before surgery as the potential optimized strategy for the treatment of a subset of patients with PDAC. In clinical practice, further study can be focused on the necessitation of neoadjuvant therapy in resectable PDAC according to patients’ mutational signatures, such as the different TP53 mutation subtypes. To achieve this, the biopsy of the tumor should be primarily done to perform genetic analysis, and relative randomized controlled trials can be undertaken in different patient subgroups with distinct TP53 mutation subtypes.
Some limitations existed in this study. Firstly, only a relative small subset of patients were operated with small-sized pancreatic tumor. Secondly, patients were only analysed based on the genomic features concerning core gene pathways instead of the multi-omic data, and the tumor heterogeneity was not taken into great consideration. Thirdly, the copy number variations including the loss of heterozygosity status of the TP53 mutations were not evaluated in this study.
Conclusions
TP53 missense mutation was associated with poor tumor differentiation, and revealed gain-of-function properties as it was predictive of early LNs involvement and distal metastasis. Nonetheless, it was not associated with insensitivity to chemotherapy, highlighting the neoadjuvant therapy before surgery as the potential optimized strategy for the treatment of a subset of patients.
Data availability
The data used and analysed in the study are available from the corresponding author on reasonable request.
References
Grossberg AJ, Chu LC, Deig CR, Fishman EK, Hwang WL, Maitra A, et al. Multidisciplinary standards of care and recent progress in pancreatic ductal adenocarcinoma. Cancer J Clin. 2020;70(5):375–403.
Sakamoto H, Attiyeh MA, Gerold JM, Makohon-Moore AP, Hayashi A, Hong J, et al. The Evolutionary origins of recurrent Pancreatic Cancer. Cancer Discov. 2020;10(6):792–805.
Connor AA, Gallinger S. Pancreatic cancer evolution and heterogeneity: integrating omics and clinical data. Nat Rev Cancer. 2022;22(3):131–42.
Hayashi A, Hong J, Iacobuzio-Donahue CA. The Pancreatic cancer genome revisited. Nat Rev Gastroenterol Hepatol. 2021;18(7):469–81.
Jones S, Zhang X, Parsons DW, Lin JC, Leary RJ, Angenendt P, et al. Core signaling pathways in human pancreatic cancers revealed by global genomic analyses. Science. 2008;321(5897):1801–6.
Bykov VJN, Eriksson SE, Bianchi J, Wiman KG. Targeting mutant p53 for efficient cancer therapy. Nat Rev Cancer. 2018;18(2):89–102.
Hingorani SR, Wang L, Multani AS, Combs C, Deramaudt TB, Hruban RH, et al. Trp53R172H and KrasG12D cooperate to promote chromosomal instability and widely metastatic pancreatic ductal adenocarcinoma in mice. Cancer Cell. 2005;7(5):469–83.
Muller PA, Caswell PT, Doyle B, Iwanicki MP, Tan EH, Karim S, et al. Mutant p53 drives invasion by promoting integrin recycling. Cell. 2009;139(7):1327–41.
Morton JP, Timpson P, Karim SA, Ridgway RA, Athineos D, Doyle B, et al. Mutant p53 drives Metastasis and overcomes growth arrest/senescence in Pancreatic cancer. Proc Natl Acad Sci USA. 2010;107(1):246–51.
Sabapathy K, Lane DP. Therapeutic targeting of p53: all mutants are equal, but some mutants are more equal than others. Nat Rev Clin Oncol. 2018;15(1):13–30.
Escobar-Hoyos LF, Penson A, Kannan R, Cho H, Pan CH, Singh RK, et al. Altered RNA splicing by mutant p53 activates oncogenic RAS Signaling in Pancreatic Cancer. Cancer Cell. 2020;38(2):198–211. e8.
Kim MP, Li X, Deng J, Zhang Y, Dai B, Allton KL, et al. Oncogenic KRAS Recruits an Expansive Transcriptional Network through Mutant p53 to Drive Pancreatic Cancer Metastasis. Cancer Discov. 2021;11(8):2094–111.
Kennedy MC, Lowe SW. Mutant p53: it’s not all one and the same. Cell Death Differ. 2022;29(5):983–7.
Pan M, Jiang C, Zhang Z, Achacoso N, Alexeeff S, Solorzano AV, et al. TP53 gain-of-function and non-gain-of-function mutations are Associated with Differential Prognosis in Advanced Pancreatic Ductal Adenocarcinoma. JCO Precis Oncol. 2023;7:e2200570.
Mizrahi JD, Surana R, Valle JW, Shroff RT. Pancreatic cancer. Lancet. 2020;395(10242):2008–20.
Springfeld C, Ferrone CR, Katz MHG, Philip PA, Hong TS, Hackert T, et al. Neoadjuvant therapy for Pancreatic cancer. Nat Rev Clin Oncol. 2023;20(5):318–37.
Zhan Q, Wen C, Zhao Y, Fang L, Jin Y, Zhang Z, et al. Identification of copy number variation-driven molecular subtypes informative for prognosis and treatment in pancreatic adenocarcinoma of a Chinese cohort. EBioMedicine. 2021;74:103716.
Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7(4):248–9.
Connor AA, Denroche RE, Jang GH, Timms L, Kalimuthu SN, Selander I, et al. Association of distinct mutational signatures with correlates of increased Immune Activity in Pancreatic Ductal Adenocarcinoma. JAMA Oncol. 2017;3(6):774–83.
Alexandrov LB, Kim J, Haradhvala NJ, Huang MN, Tian Ng AW, Wu Y, et al. The repertoire of mutational signatures in human cancer. Nature. 2020;578(7793):94–101.
Lin J, Zhai S, Zou S, Xu Z, Zhang J, Jiang L, et al. Positive feedback between lncRNA FLVCR1-AS1 and KLF10 may inhibit Pancreatic cancer progression via the PTEN/AKT pathway. J Exp Clin Cancer Res. 2021;40(1):316.
Iacobuzio-Donahue CA, Fu B, Yachida S, Luo M, Abe H, Henderson CM, et al. DPC4 gene status of the primary carcinoma correlates with patterns of failure in patients with Pancreatic cancer. J Clin Oncology: Official J Am Soc Clin Oncol. 2009;27(11):1806–13.
Winter JM, Tang LH, Klimstra DS, Liu W, Linkov I, Brennan MF, et al. Failure patterns in resected pancreas adenocarcinoma: lack of predicted benefit to SMAD4 expression. Ann Surg. 2013;258(2):331–5.
Qian ZR, Rubinson DA, Nowak JA, Morales-Oyarvide V, Dunne RF, Kozak MM, et al. Association of alterations in main driver genes with outcomes of patients with resected pancreatic ductal adenocarcinoma. JAMA Oncol. 2018;4(3):e173420.
McIntyre CA, Lawrence SA, Richards AL, Chou JF, Wong W, Capanu M, et al. Alterations in driver genes are predictive of survival in patients with resected pancreatic ductal adenocarcinoma. Cancer. 2020;126(17):3939–49.
Cenik BK, Shilatifard A. COMPASS and SWI/SNF complexes in development and Disease. Nat Rev Genet. 2021;22(1):38–58.
Piunti A, Shilatifard A. The roles of polycomb repressive complexes in mammalian development and cancer. Nat Rev Mol Cell Biol. 2021;22(5):326–45.
Yi Z, Wei S, Jin L, Jeyarajan S, Yang J, Gu Y, et al. KDM6A regulates cell plasticity and Pancreatic Cancer progression by Noncanonical Activin Pathway. Cell Mol Gastroenterol Hepatol. 2022;13(2):643–67.
Kimura Y, Fukuda A, Ogawa S, Maruno T, Takada Y, Tsuda M, et al. ARID1A maintains differentiation of pancreatic ductal cells and inhibits development of pancreatic ductal adenocarcinoma in mice. Gastroenterology. 2018;155(1):194–209e2.
Wang SC, Nassour I, Xiao S, Zhang S, Luo X, Lee J, et al. SWI/SNF component ARID1A restrains pancreatic neoplasia formation. Gut. 2019;68(7):1259–70.
Cancer Genome Atlas Research Network. Electronic address aadhe, Cancer Genome Atlas Research N. Integrated Genomic Characterization of Pancreatic Ductal Adenocarcinoma. Cancer Cell. 2017;32(2):185–203. e13.
Golan T, O’Kane GM, Denroche RE, Raitses-Gurevich M, Grant RC, Holter S, et al. Genomic features and classification of homologous recombination deficient pancreatic ductal adenocarcinoma. Gastroenterology. 2021;160(6):2119–32e9.
Grant RC, Denroche R, Jang GH, Nowak KM, Zhang A, Borgida A, et al. Clinical and genomic characterisation of mismatch repair deficient pancreatic adenocarcinoma. Gut. 2021;70(10):1894–903.
Collisson EA, Sadanandam A, Olson P, Gibb WJ, Truitt M, Gu S, et al. Subtypes of pancreatic ductal adenocarcinoma and their differing responses to therapy. Nat Med. 2011;17(4):500–3.
Moffitt RA, Marayati R, Flate EL, Volmar KE, Loeza SG, Hoadley KA, et al. Virtual microdissection identifies distinct tumor- and stroma-specific subtypes of pancreatic ductal adenocarcinoma. Nat Genet. 2015;47(10):1168–78.
Bailey P, Chang DK, Nones K, Johns AL, Patch AM, Gingras MC, et al. Genomic analyses identify molecular subtypes of Pancreatic cancer. Nature. 2016;531(7592):47–52.
Wilson SNK, Grant GW, Seto RC, O’Kane M, Vajpeyi G. Morphological classification of pancreatic ductal adenocarcinoma that predicts molecular subtypes and correlates with clinical outcome. Gut. 2020;69(2):317–28.
Connor AA, Denroche RE, Jang GH, Lemire M, Zhang A, Chan-Seng-Yue M, et al. Integration of genomic and transcriptional features in Pancreatic Cancer reveals increased cell cycle progression in Metastases. Cancer Cell. 2019;35(2):267–82. e7.
Bardeesy N, Aguirre AJ, Chu GC, Cheng KH, Lopez LV, Hezel AF, et al. Both p16(Ink4a) and the p19(Arf)-p53 pathway constrain progression of pancreatic adenocarcinoma in the mouse. Proc Natl Acad Sci U S A. 2006;103(15):5947–52.
Weissmueller S, Manchado E, Saborowski M, Morris JPt, Wagenblast E, Davis CA, et al. Mutant p53 drives Pancreatic cancer Metastasis through cell-autonomous PDGF receptor beta signaling. Cell. 2014;157(2):382–94.
Pfister NT, Prives C. Transcriptional Regulation by Wild-Type and Cancer-Related Mutant Forms of p53. Cold Spring Harb perspect Med. 2017;7(2):a026054.
Cao L, Huang C, Cui Zhou D, Hu Y, Lih TM, Savage SR, et al. Proteogenomic characterization of pancreatic ductal adenocarcinoma. Cell. 2021;184(19):5031–52e26.
Reticker-Flynn NE, Zhang W, Belk JA, Basto PA, Escalante NK, Pilarowski GOW, et al. Lymph node colonization induces tumor-immune tolerance to promote distant Metastasis. Cell. 2022;185(11):1924–42e23.
Blackford A, Serrano OK, Wolfgang CL, Parmigiani G, Jones S, Zhang X, et al. SMAD4 gene mutations are associated with poor prognosis in pancreatic cancer. Clin cancer Research. 2009;15(14):4674–9.
Oshima M, Okano K, Muraki S, Haba R, Maeba T, Suzuki Y, et al. Immunohistochemically detected expression of 3 major genes (CDKN2A/p16, TP53, and SMAD4/DPC4) strongly predicts survival in patients with resectable Pancreatic cancer. Ann Surg. 2013;258(2):336–46.
Versteijne E, van Dam JL, Suker M, Janssen QP, Groothuis K, Akkermans-Vogelaar JM, et al. Neoadjuvant Chemoradiotherapy Versus Upfront Surgery for Resectable and Borderline Resectable Pancreatic Cancer: long-term results of the Dutch Randomized PREOPANC Trial. J Clin Oncology: Official J Am Soc Clin Oncol. 2022;40(11):1220–30.
Acknowledgements
Not applicable.
Funding
Not applicable.
Author information
Authors and Affiliations
Contributions
YZ, JJ and JZ designed the study, collected the patient data, and drafted the paper. YZ. and YJ conducted the functional experiments. DC, XD and BS participated in the design of the study and edited the final paper. YJ, NF and MC collected the patient data. JW, KQ and YJ participated in the statistical analysis and data interpretation.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The study was approved by the ethics committee of Shanghai Ruijin Hospital and all patients signed informed consent.
Consent for publication
All authors agree to publish.
Competing interests
All the authors have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1.
Raw Data 1_Annotations of mutations.
Additional file 2.
Raw Data 2_Information of patients.
Additional file 3.
Raw Data 3_Information of the TCGA cohort.
Additional file 4: Table S1–S14.
Table S1. The Detailed Pathways and Genes Related to PDAC Carcinogenesis and Included in the Genetic Analyses Panel. Table S2. The siRNA Sequences and the Primer Sequences Used in the Study. Table S3. Patient Cohort and KRAS Status. Table S4. Clinical Information of the Patient Cohort. Table S5. Core Gene Pathway Alterations in KRASmut and KRASWT PDAC. Table S6. Core Gene Pathway Alterations and Association with Tumor Differentiation in age-related KRASmut PDAC. Table S7. Multivariate Logistic Regression of Core Gene Pathway Alterations for Tumor Differentiation in KRASmut PDAC. Table S8. Other Gene Pathway Alterations and Association with Lymph Nodes Involvement. Table S9. Frequency of Core Gene Pathway Alterations according to Tumor Size in KRASmut PDAC. Table S10. Other Gene Pathway Alterations and Association with Distal Metastasis. Table S11. Core Gene Pathway Alterations and Association with Lymph Nodes Involvement in Age-related PDAC. Table S12. Frequency of Core Gene Pathway Alterations according to Tumor Size in Age-related KRASmut PDAC. Table S13. Core Gene Pathway Alterations and Association with Distal Metastasis in Age-related PDAC. Table S14. Core Gene Alterations and Association with Tumor Differentiation in KRASmut PDAC in the TCGA cohort.
Additional file 5:Â Figure S1.
Study design and flowchart of the study. Figure S2. Examples of the histological features of A well, B moderate and C poor differentiated PDAC. Figure S3. Illustration of how alterations of the Trithorax genes might result in the repression of TP53 function. Figure S4. The stepwise mode of PDAC carcinogenesis and relative genetic alterations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Zhou, Y., Jin, J., Ji, Y. et al. TP53 missense mutation reveals gain-of-function properties in small-sized KRAS transformed pancreatic ductal adenocarcinoma. J Transl Med 21, 872 (2023). https://doi.org/10.1186/s12967-023-04742-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12967-023-04742-y