ORIGINAL RESEARCH article

Front. Oncol., 24 June 2021

Sec. Women's Cancer

Volume 11 - 2021 | https://doi.org/10.3389/fonc.2021.628814

Mutant p53 and Twist1 Co-Expression Predicts Poor Prognosis and Is an Independent Prognostic Factor in Breast Cancer

  • 1. Department of Breast-Thyroid-Surgery and Cancer Research Center, Xiang’an Hospital of Xiamen University, Xiamen, China

  • 2. Clinical Central Research Core, School of Medicine, Xiang’an Hospital of Xiamen University, Xiamen, China

  • 3. Key Laboratory for Endocrine-Related Cancer Precision Medicine of Xiamen, Xiang’an Hospital of Xiamen University, Xiamen, China

  • 4. Cancer Research Center, School of Medicine, Xiamen University, Xiamen, China

  • 5. Guangdong Provincial Key Laboratory for Breast Cancer Diagnosis and Treatment, Cancer Hospital of Shantou University Medical College, Shantou, China

  • 6. Department of Pathology, Cancer Hospital of Shantou University Medical College, Shantou, China

  • 7. Department of Breast Center, Cancer Hospital of Shantou University Medical College, Shantou, China

  • 8. Department of Medical Oncology, Xiang’an Hospital of Xiamen University, Xiamen, China

Abstract

Purpose:

The basic helix-loop-helix transcription factor (bHLH) transcription factor Twist1 plays a key role in embryonic development and tumorigenesis. p53 is a frequently mutated tumor suppressor in cancer. Both proteins play a key and significant role in breast cancer tumorigenesis. However, the regulatory mechanism and clinical significance of their co-expression in this disease remain unclear. The purpose of this study was to analyze the expression patterns of p53 and Twist1 and determine their association with patient prognosis in breast cancer. We also investigated whether their co-expression could be a potential marker for predicting patient prognosis in this disease.

Methods:

Twist1 and mutant p53 expression in 408 breast cancer patient samples were evaluated by immunohistochemistry. Kaplan-Meier Plotter was used to analyze the correlation between co-expression of Twist1 and wild-type or mutant p53 and prognosis for recurrence-free survival (RFS) and overall survival (OS). Univariate analysis, multivariate analysis, and nomograms were used to explore the independent prognostic factors in disease-free survival (DFS) and OS in this cohort.

Results:

Of the 408 patients enrolled, 237 (58%) had high mutant p53 expression. Two-hundred twenty patients (53.9%) stained positive for Twist1, and 188 cases were Twist1-negative. Furthermore, patients that co-expressed Twist1 and mutant p53 (T+P+) had significantly advanced-stage breast cancer [stage III, 61/89 T+P+ (68.5%) vs. 28/89 T-P- (31.5%); stage II, 63/104 T+P+ (60.6%)vs. 41/104 T-P- (39.4%)]. Co-expression was negatively related to early clinical stage (i.e., stages 0 and I; P = 0.039). T+P+ breast cancer patients also had worse DFS (95% CI = 1.217–7.499, P = 0.017) and OS (95% CI = 1.009–9.272, P = 0.048). Elevated Twist1 and mutant p53 expression predicted shorter RFS in basal-like patients. Univariate and multivariate analysis identified three variables (i.e., lymph node involvement, larger tumor, and T+P+) as independent prognostic factors for DFS. Lymph node involvement and T+P+ were also independent factors for OS in this cohort. The total risk scores and nomograms were reliable for predicting DFS and OS in breast cancer patients.

Conclusions:

Our results revealed that co-expression of mutant p53 and Twist1 was associated with advanced clinical stage, triple negative breast cancer (TNBC) subtype, distant metastasis, and shorter DFS and OS in breast cancer patients. Furthermore, lymph nodes status and co-expression of Twist1 and mutant p53 were classified as independent factors for DFS and OS in this cohort. Co-evaluation of mutant p53 and Twist1 might be an appropriate tool for predicting breast cancer patient outcome.

Introduction

Breast cancer (BC) is a heterogeneous disease that is classified into different subtypes depending on the expression of the estrogen and progesterone receptors and the amplification of human epidermal growth factor receptor 2 (HER2/ERBB2) (). In particular, molecular profiling has divided BC into five subtypes: luminal A, luminal B, HER2+, and basal-like and triple-negative breast cancer (TNBC) (). There are no effective targeted drugs for TNBC due to the lack of estrogen receptor (ERα) or HER2 overexpression. TNBC is one of the most heterogeneous breast cancer subtypes, accounting for 15–20% of all breast cancer. It has poor clinical and pathological features and is considered more aggressive than other BC subtypes, such as luminal A and luminal B (). Therefore, it is important to identify the molecules and signaling pathways related to TNBC progression. Recent studies demonstrated that the epithelial-to-mesenchymal transition (EMT) process promotes tumor invasion and metastasis in different solid tumor types (). Thus, it is important to identify the transcription factors involved in EMT of TNBC progression (). EMT is a process characterized by the loss of cell-to-cell adhesion and tight cell junctions, and increased motility and invasion. Twist1 is a major regulator of EMT in breast cancer (). We previously revealed that Twist1 was overexpressed in breast cancers with higher node status and clinical stage and positively associated with EMT in breast cancer ().

Mutation of the tumor suppressor gene p53 is found in various cancers, including breast, melanoma, bladder, and colon, particularly in tumors resistant to endocrine therapy or radiotherapy (). Altered p53 protein is prevalently associated with TNBC, and loss of p53 function has been linked to the induction of EMT. Recently, several studies demonstrated that p53-regulated EMT is a critical process in cancer metastasis and chemotherapy resistance. Chang et al. () showed that p53 regulates EMT by activating the microRNA miR-200c. Kim et al. () demonstrated that p53 loss-of-function or mutation promotes EMT by derepressing Snail. Most noteworthy, van Nes et al. () found that Twist1 overexpression is correlated with invasive lobular carcinoma, which is consistent with p53 changes in breast cancer and establishes a mechanistic link between Twist1, p53, and tumor progression. However, studies on the relationship between mutant p53 and Twist1 and their combined role in determining the prognosis and survival of breast cancer patients are scarce. In the present study, we investigated the clinical significance and the combined detection value of mutant p53 and Twist1 in breast cancer patients. We report that mutant p53 and Twist1 co-expression could represent a feasible and effective marker for predicting patient prognosis and survival in breast cancer.

Materials and Methods

Patients

A total of 408 breast cancer patients treated at the Cancer Hospital of Shantou University Medical College (SUMC) between 2011 and 2013 were enrolled in this study. None of the patients received any treatment before surgery. Clinicopathological data, including age, menopausal status, tumor size, nodal status, TNM stage, histology, and morphology, were collected. All samples were fixed in buffered formalin and embedded in paraffin for immunohistochemistry (IHC). The DFS and the OS of the patients were quantified. This study was approved by the Medical Ethics Committee of the Cancer Hospital of Shantou University Medical College.

Immunohistochemistry

All samples were fixed in 10% neutral buffered formalin for 8 to 48 h, followed by dehydration with alcohol and xylene. The dehydrated samples were embedded in paraffin. For immunohistochemical staining, 4-μm-thick paraffin-embedded sections were obtained. The sections were treated with 0.3% H2O2 for 10 min at room temperature. For antigen retrieval, the sections were heated in 0.01 mmol/L sodium citrate (pH 6.0) in a microwave oven. Sections were then incubated with primary antibodies at 4°C overnight in a humid environment. The stained slides were rinsed thrice in 0.1 mmol/L PBS for 2 min and then incubated for 30 min at room temperature with horseradish peroxidase-conjugated goat anti-rabbit or mouse secondary antibody. Visualization of the stained bands was performed with 3’ 3-diaminobenzidine. Finally, the sections were stained for 5 min with 3,3′-diaminobenzidine (DAB; Fuzhou Maixin Biotechnology Development Co., Ltd., China; DAB-0031) and counterstained with hematoxylin for 1 min. The IHC antibodies were as follows: mutant p53 (Cat. No. 0010-2, Fuzhou Maixin Biotech, Fuzhou, China) and Twist1 (Cat. No. ab50581, Abcam, Cambridge, UK). All immunostained sections were evaluated in a coded manner by the principal author, who was blinded to the patients’ clinicopathologic data. Twist1 expression was detected mainly in the cell membrane and cytoplasm. The expressions were semi-quantitatively determined according to the percentage of positive cancer cells. Staining intensity was classified as four grades: none (0), weak (1), moderate (2) and strong (3). The percentage of positive cancer cells was classified as 4 grades: 0 (0%), 1 (1%-10%), 2 (11%-49%) and 3 (50%-100%). The total score was a product of two scores and the final score of one sample was the mean of 10 microscopic fields. The median score was determined, according to which cancers were categorized into Negative- (score 0-3) and Positive-expression (score 4-6) cancers ().Brown stained cells in the nucleus were defined as mutant p53 positive cells and the specimens were divided into the following categories: Mutant p53-negative (-), no brown stained; weakly positive (+), <30% cells were mutant p53-positive; moderately positive (++), 30-70% cells were mutant p53-positive; and strongly positive (+++), >70% cells were mutant p53-positive. “-” and “+” were defined as Mutant p53-negative, “++” and “+++” were defined as Mutant p53-positive (). The following criteria are the definition of molecular typing of breast cancer:TNBC: ER-/PR-/HER2-;HER-2 enriched: ER-/PR-/HER2+;Luminal A: ER+ and/or PR+, HER2-, Ki-67 low expression(≤ 14%);Luminal B: ER+ and/or PR+, HER2-, Ki-67 high expression (>14%) (, ).

Kaplan-Meier Analysis

The free online clinical database, Kaplan-Meier Plotter (http://kmplot.com), was used to analyze the correlation between Twist1 co-expression with wild-type or mutant p53 and prognosis for recurrence-free survival in patients suffering from breast cancer.

Statistical Analysis

All data were analyzed using SPSS 13.0 statistical software (SPSS Inc., Chicago, IL). The counting data are expressed as the rate using Pearson’s χ2 test. Spearman was used for correlation analysis. Survival curves were assessed using the Kaplan-Meier method, and the data were compared using the log-rank test. Univariate and multivariate analyses were used to determine the impact of the variables on patient survival. Two-sided P-values < 0.05 were considered statistically significant.

Results

Intersection of Mutant p53 and Twist1 in Clinical and Pathological Features of Breast Cancer

We stained 408 breast cancer samples to determine the expression of mutant p53 and Twist1 (Figures 1A–D) and IgG was used as negative control (Figure 1E). Mutant p53 expression was positively associated with invasive ductal carcinoma (IDC) (213/355, 60%, P = 0.043). A significant positive relationship was observed between mutant p53 expression and HER2 overexpression (81/119, 68.1%, P = 0.009), the luminal B subtype (103/180, 57.2%) and followed by TNBC (36/63, 57.1%) (P = 0.004). There were no significant differences with age, menopausal status, lymph node status, TNM stage, histological grade, or ER or PR status. Of the samples analyzed, 220/408 (53.9%) were Twist1-positive. The Twist1-positive rate was significantly higher in ER- (113/147; 76.9%) than in ER+ (107/261; 41.0%) breast cancer (P < 0.001). Similarly, it was higher in PR- (127/180; 70.6%) than in PR+ (93/228; 40.8%) breast cancer (P < 0.001). Thus, Twist1 was negatively correlated with ER and PR expression. Twist1 expression was significantly higher in TNBC (55/63, 87.3%), follow by HER2-enriched (43/119, 63.9%) and the luminal B (79/180, 43.9%) and luminal A (10/46, 21.7%) subtypes (P < 0.001). In contrast, no significant differences were seen with lymph nodal status, TNM stage, histological type, grade, age, or menopausal status. However, Twist1 was positively correlated with mutant p53, E-cadherin, and Ki67 (P < 0.001). The relationships between Twist1, mutant p53, E-cadherin, and other biomarkers are shown in Tables 1 and 2.

Figure 1

Table 1

FeaturesP53 (%)PTwist1 (%)P
NegativePositiveNegativePositive
Age
≤ 60133 (40.5)195(59.5)0.259154 (47.0)174(53.0)0.474
> 6038 (47.5)42 (52.5)34(42.5)46(57.5)
Menopausal status
Pre89 (41.4)126(58.6)0.82396(44.7)119(55.3)0.542
Post82 (42.5)111(57.5)92 (47.7)101(52.3)
Tumor size
Tis-T152 (45.2)63 (54.8)0.48363 (54.8)52 (45.2)0.073
T2-T3102(39.7)155(60.3)108 (42.0)149 (58)
T417 (47.2)19 (52.8)17 (47.2)19 (52.8)
Nodal status
N088 (45.1)107(54.9)0.64699(50.8)96 (49.2)0.294
N137 (38.9)58 (61.1)42(44.2)53 (55.8)
N228 (40.0)42 (60.0)28 (40.0)42 (60.0)
N318 (37.5)30 (62.5)19(39.6)29 (60.4)
Histological grade
I23 (36.5)40 (63.5)0.42631 (49.2)32 (50.8)0.285
II72 (45.6)86 (54.4)79 (50.0)79 (50.0)
III76 (40.9)110(59.1)78 (41.9)108(58.1)
Histological type
IDC142(40.0)213(60.0)0.043161 (45.4)194(54.6)0.446
Non IDC29 (54.7)24 (45.3)27 (50.9)26(49.1)
ER
Negative57 (38.8)90 (61.2)0.33534 (23.1)113(76.9)<0.001
Positive114(43.7)147(56.3)154 (59.0)107(41.0)
PR
Negative77 (42.8)103(57.2)0.75353 (29.4)127(70.6)<0.001
Positive94 (41.2)134(58.8)135 (59.2)93 (40.8)
HER2
Negative133(46.0)156(54.0)0.009145 (50.2)144(49.8)0.010
Positive38 (31.9)81 (68.1)43 (36.1)76 (63.9)
Molecular subtypes
Luminal A29 (63.0)17 (37.0)0.00436 (78.3)10(21.7)<0.001
Luminal B77 (42.8)103(57.2)101(56.1)79 (43.9)
HER2-enriched38 (31.9)81 (68.1)43 (36.1)76 (63.9)
TNBC27 (42.9)36 (57.1)8 (12.7)55 (87.3)

Association between mutant p53 and Twist1 expression and clinicopathological characteristics in breast cancer patients (n = 408).

Statistically significant values (P < 0.05) are in bold.

Table 2

BiomarkerTwist1 (%)χ2P
Negative (n = 188)Positive (n = 220)
P53
Negative99(57.9)72(42.1)16.543<0.001
Positive89(37.6)148(62.4)
Ki67
≤ 14%45(78.9)12(21.1)28.810<0.001
> 14%143(40.7)208(59.3)
E-cadherin
Negative9(11.2)71(88.8)48.583<0.001
Positive179(54.6)149(45.4).

Association between Twist1 expression and other molecules in breast cancer patients.

Statistically significant values (P < 0.05) are in bold.

Combined Analysis of the Clinical and Pathological Features of Breast Cancer With Mutant p53 and Twist1

Because mutant p53 expression was positively related to Twist1, we classified the study patients into four groups: Twist1- and mutant p53- (T-P-), Twist1+ and mutant p53- (T+P-), Twist1- and mutant p53+ (T-P+), and Twist1+ and mutant p53+ (T+P+). Of the 408 patients, 148 patients belonged to the T+P+ group, while 99 patients were assigned to the T-P- group. It is of note that 37.5% of stage 0 patients (3/8), 45.7% of stage I patients (21/46), 60.6% of stage II patients (63/104), and 68.5% of stage III patients (61/89) were classified as T+P+, indicating that the T+P+ frequency was significantly higher in advanced clinical stages (P = 0.039, Table 3).

Table 3

FeaturesCase (%)χ2PCase (%)χ2P
T-P-T+P+T+P-T-P+
(n = 99)(n = 148)(n = 72)(n = 89)
Age
≤6079(39.7)120(60.3)0.0620.80354(41.9)75(58.1)2.1480.143
>6020(41.7)28(58.3)18(56.3)14(43.8)
Menopausal status
Pre50(38.5)80(61.5)0.3000.58439(45.9)46(54.1)0.0980.754
Post49(41.9)68(58.1)33(43.4)43(56.6)
Tumor size
Tis-T134 (50.0)34(50.0)4.5520.10318(38.3)29(61.7)1.1600.560
T2-T356(35.2)103(64.8)46(46.9)52(53.1)
T49(45.0)11(55.0)8(50.0)8(50.0)
Nodal status
N054(46.6)62(53.4)4.1990.24134(43.0)45(57.0)0.5040.918
N122(36.7)38(63.3)15(42.9)20(57.1)
N215(34.1)29(65.9)13(50.0)13(50.0)
N38(29.6)19(70.4)10(47.6)11(52.4)
Clinical stage
05(62.5)3(37.5)8.3450.0391(33.3)2(66.7)0.3440.076
I25(54.3)21(45.7)18(46.2)21(53.8)
II41(39.4)63(60.6)27(42.9)36(57.1)
III28(31.5)61(68.5)26(46.4)30(53.6)
Histological grade
I15(38.5)24(61.5)2.3080.3158(33.3)16(66.7)1.7000.427
II43(46.2)50(53.8)29(44.6)36(55.4)
III41(36.0)73(64.0)35(48.6)37(51.4)
Histological type
IDC17(54.8)14(45.2)3.2150.07312(54.5)10(45.5)0.9950.319
Non-IDC82(38.0)134(62.0)60(43.2)79(56.8)
ER
Negative19(20.2)75(79.8)24.943<0.00138(71.7)15(28.3)23.260<0.001
Positive80(52.3)73(47.7)34(31.5)74(68.5)
PR
Negative30(27.3)80(72.7)13.547<0.00147(67.1)23(32.9)25.187<0.001
Positive69(50.4)68(49.6)25(27.5)66(72.5)
HER2
Negative84(46.9)95(53.1)12.69<0.00149(44.5)61(55.5)0.0040.948
Positive15(22.1)53(77.9)23(45.1)28(54.9)
Molecular subtypes
Luminal A24(82.8)5(17.2)44.415<0.0015(29.4)12(70.6)27.07<0.001
Luminal B53(49.1)55(50.9)24(33.3)48(66.7)
HER2-enriched15(22.1)53(77.9)23(45.1)28(54.9)
TNBC7(16.7)35(83.3)20(95.2)1(4.8)

Correlations between mutant p53 and Twist1 expression and clinicopathologic parameters.

Statistically significant values (P < 0.05) are in bold.

A comparison between T+P+ and T-P- groups showed that more ER- patients (75/94, 79.8%) and PR- patients (80/110, 72.7%) were T+P+ (P < 0.001). Similarly, significantly higher mutant p53 and Twist1 expression levels were observed with the HER2+ (77.9% T+P+ vs. 22.1% T-P-, P < 0.001) and TNBC (83.3% T+P+ vs. 16.7% T-P-, P < 0.001) subtypes.

As shown in Table 4, patients with T+P+ had higher Ki67, E-cadherin, and VEGF-C expression than the T-P- group (P < 0.001). There were no statistically significant differences between the T+P- and T-P+ groups in terms of age, menopausal status, tumor size or clinical stage.

Table 4

BiomarkerCase (%)χ2PCase (%)χ2P
T-P-T+P+T+P-T-P+
(n = 72)(n = 172)(n = 48)(n = 116)
Ki67
≤ 14%30(83.3)6(16.7)32.827<0.0016(28.6)15(71.4)2.5480.110
> 14%69(32.7)142(67.3)66(47.1)74(52.9)
E-cadherin
Negative4(8.7)42(91.3)23.185<0.00129(85.3)5(14.7)28.702<0.001
Positive95(47.3)106(52.7)43(33.9)84(66.1)
VEGF-C
Negative45(57.0)34(43.0)13.782<0.00114(34.1)27(65.9)2.4880.115
Positive54(32.1)114(67.9)58(48.3)62(51.7)

Correlations of mutant p53 and Twist1 expression with other molecules in breast cancer patients.

Statistically significant values (P < 0.05) are in bold.

Univariate and Multivariate Survival Analysis for Predicting DFS and OS

We next assessed the prognostic value of Twist1 with or without p53 mutation in breast cancer using the Kaplan-Meier plotter. The results demonstrated that elevated Twist1 expression combined with mutant p53 predicted shorter RFS for basal-like breast cancer patients (P = 0.027) (Figure 2G); however, there was no difference for the OS (P = 0.18) (Figure 3G). There were no statistical differences in RFS and OS for the other groups of patients (Figures 2 and 3).

Figure 2

Figure 3

To investigate which factors could predict the clinical outcome of breast cancer patients, we performed univariate and multivariate survival analysis. DFS was followed up for 2 to 74 months and OS for 3 to 74 months. In our cohort, 148 patients belonged to the T+P+ group, while 99 patients were T-P-. The variables included in the analysis were tumor size, node stage, molecular subtype, and the different combinations of Twist1 and mutant p53 (i.e., T-P-, T+P-, T-P+, and T+P+). Tables 5 and 6 show the main results of the survival analysis. By univariate analysis, we found the node status was positively related with DFS (N1 vs. N0, HR = 3.749, 95% CI = 1.476–9.524, P = 0.005; N2 vs. N0, HR = 7.608, 95% CI = 3.154–18.350, P < 0.001; N3 vs. N0, HR = 13.977 95% CI = 5.871–33.277, P < 0.001). The multivariate analysis demonstrated that node status was an independent factor that influenced DFS (N1 vs. N0, HR = 2.757, 95% CI = 1.069–7.112, P = 0.036; N2 vs. N0, HR = 5.386, 95%CI = 2.134–13.595, P < 0.001; N3 vs. N0, HR = 10.757, 95% CI = 4.286–27.000, P < 0.001). A larger tumor size was another independent risk factor for DFS in breast cancer (Table 5). Multivariate analysis demonstrated that node status was also an independent factor influencing OS (N1 vs. N0, HR = 8.522 95% CI = 1.789–40.606, P = 0.007; N2 vs. N0, HR = 13.717, 95% CI = 2.780–67.686, P = 0.001; N3 vs. N0, HR = 41.619, 95% CI = 8.943–193.679, P < 0.001). However, there were no significant differences between tumor size, molecular typing, and OS by multivariate analysis. Patients with co-expression of mutant p53 and Twist1 (T+P+) had worse DFS (P = 0.017) and OS (P = 0.018) compared to patients in the other groups (i.e., T-P-, T+P-, and T-P+) (Figures 4A, B).

Table 5

FeaturesUnivariateMultivariate
Hazard ratio95% CIPHazard ratio95% CIP
Tumor Size
T2-T3 vs.T0-T19.8432.381,40.6970.0024.3991.028, 18.8250.046
T4 vs. T0-T121.5164.815,96.151<0.0017.3661.547,35.0670.012
Nodal status
N1 vs. N03.7491.476, 9.5240.0052.7571.069,7.1120.036
N2 vs. N07.6083.154,18.350<0.0015.3862.134,13.595<0.001
N3 vs. N013.9775.871,33.277<0.00110.7574.286,27.000<0.001
Molecular subtype
Luminal B vs. luminal A1.8400.549, 6.1720.3231.1070.318,3.8560.873
Her2 vs. Luminal A2.7910.829, 9.3960.0971.9400.548,6.8690.304
TNBC vs. Luminal A2.8990.808,10.3940.1022.3120.604,8.8520.221
Case
T+P- vs. T-P-2.6120.966,7.0620.0591.7240.627,4.7390.291
T-P+ vs. T-P-1.9070.693,5.2490.2111.6960.613,4.6870.309
T+P+ vs. T-P-3.5771.480,8.6480.0053.0201.217,7.4990.017

Disease-free survival of patients diagnosed with breast cancer.

Statistically significant values (P < 0.05) are in bold.

Table 6

FeaturesUnivariateMultivariate
Hazard ratio95% CIPHazard ratio95% CIP
Tumor Size
T2-T3 vs. T0-T16.3541.511, 26.7200.0121.8680.422, 8.2730.410
T4 vs. T0-T114.9613.232,69.2590.0013.0780.606,15.6370.175
Node status
N1 vs. N09.9302.145,45.9710.0038.5221.789,40.6060.007
N2 vs. N015.0123.289,68.524<0.00113.7172.780,67.6860.001
N3 vs. N040.3189.308,174.632<0.00141.6198.943,193.679<0.001
molecular subtype
luminal B vs. luminal A1.5710.352,7.0220.5540.8250.174,3.9170.809
Her2 vs. luminal A2.8390.645,12.4940.1681.9650.416,9.2810.394
TNBC vs. luminal A3.9150.858,17.8760.0783.0650.612,15.3440.173
Case
T+P- vs. T-P-2.7790.837,9.2280.0951.6940.500,5.7420.398
T-P+ vs. T-P-1.4040.377,5.2280.6131.1940.319,4.4730.793
T+P+ vs. T-P-3.9141.343,11.4120.0123.0591.009,9.2720.048

Overall survival of patients diagnosed with breast cancer.

Statistically significant values (P < 0.05) are in bold.

Figure 4

Importantly, the multivariate analysis demonstrated that the T+P+ group had a worse DFS (95% CI = 1.217–7.499, P = 0.017) and OS (95% CI = 1.009–9.272, P = 0.048) than the T-P-group, but the T+P- and T-P+ groups were not statistically different from the T-P- group. Univariate and multivariate analysis identified three variables (high node stage, larger tumor size, and Twist1 and mutant p53 co-expression) as independent risk factors for DFS and node involvement and T+P+ as independent risk factors for OS in breast cancer. Therefore, total risk scores and nomograms were considered suitable for predicting the DFS and OS of breast cancer patients (Figure 5).

Figure 5

Discussion

Twist1 is a master transcription factor that induces migration and invasion and promotes EMT in various cancer cells (). Wild-type p53 is a transcription factor that promotes cell cycle arrest, DNA repair, cellular senescence, and apoptosis (). p53 inhibits cancer cell proliferation and metastasis; however, mutant TP53 has a transforming function. Previous studies have demonstrated that mutant p53 is overexpressed and accumulates to high levels in cancer cells, and is highly correlated with EMT (). Moreover, Twist1 overexpression can inhibit oncogene-induced premature senescence by blocking key regulators of the p53 and Rb-dependent pathways (). Numerous studies have shown that there is a regulatory relationship between Twist1 and p53. For instance, Shiota et al. (26) suggested that Twist1 suppresses the DNA-binding activity of p53. Moreover, Kogan-Sakin et al. (27) revealed that Twist1 might be upregulated following p53 mutation in cancer cells. However, there are no reports on the combined detection of Twist1 and mutant p53 at the protein level and its prognostic impact on breast cancer patients.

We analyzed the expression patterns of mutant p53 and Twist1 in breast cancer to explore the role of these molecular markers in breast carcinogenesis, their influence on prognosis, and clinical practicability in breast cancer. Our study showed that mutant p53 expression was positively associated with invasive ductal carcinoma (IDC). Mutant p53 was highest in the HER2+ breast cancer subtype, followed by the luminal B and TNBC subtypes. Twist1 positivity was significantly higher in TNBC patients. Consistent with previous studies (, 28), we further confirmed that high Twist1expression was associated with larger tumor size, higher node involvement, and shorter DFS. Most importantly, we found, for the first time, that Twist1 was positively correlated with mutant p53, E-cadherin, and Ki67. The association of mutant p53 with Twist1 in clinical tissues was supported by the work of Yang-Hartwich et al. (29), who showed that wild-type p53 promoted Twist1 protein degradation and inhibited EMT, maintaining the epithelial phenotype. Interestingly, a preclinical study by Li et al. (30) found that Twist1 had noteworthy impacts on p53 and p21 induction, which were similar to our current results. Maestro et al. (31) found that Twist regulated p53 indirectly by modulating the ARF/MDM2/p53 pathway. Similarly, Piccinin et al. (32) found that Twist1 binds to the p53 C terminus through the Twist box to inactivate p53 in mesenchymal tumors. Conversely, Twist promotes the reprogramming of glucose metabolism in MCF10A-Twist cells and Twist-positive breast cancer cells by inhibiting the p53 pathway (33). Taken together, these findings suggest that a mutual regulation mechanism might exist between mutant p53 and Twist1. However, the clinical significance of co-expression of two proteins has not been clearly defined.

In this study, we also investigated the relationship and clinical significance of mutant p53 and Twist1. Our results suggested that co-overexpression of mutant p53 and Twist1 in breast cancer was significantly associated with larger tumor size, greater lymph node involvement, and more advanced TNM stage. Co-expression of mutant p53 and Twist1 (T+P+ group) was significantly higher in TNBC, followed by HER2-enriched, suggesting the aggressiveness of T+P+ tumors. Previously, Twist 1 overexpression was shown to be associated with poor prognosis in breast cancer (, 34). Similarly, abnormal p53 expression has been widely accepted as a poor prognostic factor in breast cancer (35). Moreover, Sanambar et al. (36) found that p53 mutation detected by immunohistochemistry could effectively predict the prognosis of breast cancer patients. In our study, regression analysis was used to assess the role of mutant p53 by immunohistochemistry. As expected, the patients in the T+P+ group presented the worst DFS and OS by univariate analysis. Furthermore, our results demonstrated that co-expression of Twist1 and mutant p53 (T+P+) is an independent prognostic factor for both DFS and OS in breast cancer patients using multivariate analysis. Similarly, Kaplan-Meier analysis using the Kaplan-Meier plotter database demonstrated that co-expression of Twist1 and mutant p53 predicted shorter RFS in patients with basal-like breast cancer. However, due to the small sample size of the database, we did not find a statistical difference between Twist1 with mutant or wild-type p53 and RFS and OS in other subtypes.

A nomograph was established based on the survival analysis, and risk scores for DFS and OS time were calculated according to the Cox regression coefficient. Internal validation patients were randomly selected from the total population, and the verification confirmed that the risk score was related to DFS and OS time, showing the nomogram’s reliability. However, the weakness is the lack of validation queues in our research. Future research should enroll more cases to further verify the nomogram.

This study is the first report showing that co-expression of Twist1 and mutant p53 could be used to evaluate treatment efficacy and prognosis in breast cancer patients. The relationship between the two proteins needs further study, such as the relationship between the specific mutation site of the p53 gene and breast cancer and the interactive regulatory mechanism of mutant p53 and Twist1. Using small molecular inhibitors that inhibit the expression of Twist1 and mutant p53 may be one strategy for targeted treatment of TNBC (37).

Funding

This work was partly supported by the National Natural Science Foundation of China (No. 81602345 and No. 91859120), the Science and Technology Planning Project of Guangdong Province (No. 2016A020215145), the Science and Technology Planning Project of Shantou (No. 190917085269842), the Start-Up Funds from Xiamen University, the Science and Technology Project of Xiamen Municipal Bureau of Science and Technology (No.3502Z20199047), the Natural Science Foundation of Fujian Province of China (No. 2020J01015), Fujian Major Scientific and Technological Special Project for “Social Development” (No.2020YZ016002).

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by Medical Ethics Committee of the Cancer Hospital of Shantou University Medical College. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author contributions

Conception and design: G-JZ and X-LW. Administrative support: G-JZ. Provision of study materials or patients: G-JZ. Collection and assembly of data: Y-QZ. Data analysis and interpretation: Y-QZ and FZ. All authors contributed to the article and approved the submitted version.

Acknowledgments

The authors would like to thank the BioMed Proofreading® LLC (USA) for language editing.

Conflict of interest

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

References

  • 1

    RussnesHGLingjærdeOCBørresen-DaleALCaldasC. Breast Cancer Molecular Stratification: From Intrinsic Subtypes to Integrative Clusters. Am J Pathol (2017) 187(10):2152–62. doi: 10.1016/j.ajpath.2017.04.022

  • 2

    TurkozFPSolakMPetekkayaIKeskinOKertmenNSariciFet al. Association Between Common Risk Factors and Molecular Subtypes in Breast Cancer Patients. Breast (2013) 22(3):344–50. doi: 10.1016/j.breast.2012.08.005

  • 3

    PodoFBuydensLMCDeganiHHilhorstRKlippEGribbestadISet al. Triple-Negative Breast Cancer: Present Challenges and New Perspectives. Mol Oncol (2010) 4(3):209–29. doi: 10.1016/j.molonc.2010.04.006

  • 4

    DongreARashidianMReinhardtFBagnatoAKeckesovaZPloeghHLet al. Epithelial-to-Mesenchymal Transition Contributes to Immunosuppression in Breast Carcinomas. Cancer Res (2017) 77(15):3982–9. doi: 10.1158/0008-5472.CAN-16-3292

  • 5

    ZahediAPhandthong RChailiARemarkGTalbotPet al. Epithelial-to-Mesenchymal Transition of A549 Lung Cancer Cells Exposed to Electronic Cigarettes. Lung Cancer (2018) 122:224–33. doi: 10.1016/j.lungcan.2018.06.010

  • 6

    ZhouPLiBLiuFZhangMWangQLiuYet al. The Epithelial to Mesenchymal Transition (EMT) and Cancer Stem Cells: Implication for Treatment Resistance in Pancreatic Cancer. Mol Cancer (2017) 16(1):52. doi: 10.1186/s12943-017-0624-9

  • 7

    KhaledNBidetY. New Insights Into the Implication of Epigenetic Alterations in the EMT of Triple Negative Breast Cancer. Cancers (Basel) (2019) 11(4):559. doi: 10.3390/cancers11040559

  • 8

    TaubeJHHerschkowitzJIKomurovKZhouAYGuptaSYangJet al. Core Epithelial-to-Mesenchymal Transition Interactome Gene-Expression Signature is Associated With Claudin-Low and Metaplastic Breast Cancer Subtypes. Proc Natl Acad Sci U S A (2010) 107(35):15449–54. doi: 10.1073/pnas.1015095107

  • 9

    ZhangY-QWeiX-LLiangY-KChenW-LZhangFBaiJ-Wet al. Over-Expressed Twist Associates With Markers of Epithelial Mesenchymal Transition and Predicts Poor Prognosis in Breast Cancers Via ERK and Akt Activation. PloS One (2015) 10(8):e0135851. doi: 10.1371/journal.pone.0135851

  • 10

    AgarwalAKasinathanAGanesanRBalasubramanianABhaskaranJSureshSet al. Curcumin Induces Apoptosis and Cell Cycle Arrest Via the Activation of Reactive Oxygen Species-Independent Mitochondrial Apoptotic Pathway in Smad4 and p53 Mutated Colon Adenocarcinoma HT29 Cells. Nutr Res (2018) 51:6781. doi: 10.1016/j.nutres.2017.12.011

  • 11

    Shattuck-BrandtRLChenS-CMurrayEJohnsonCACrandallHO'NeaJFet al. Metastatic Melanoma Patient-Derived Xenografts Respond to MDM2 Inhibition as a Single Agent or in Combination With BRAF/MEK Inhibition. Clin Cancer Res (2020) 26(14):3803–18. doi: 10.1158/1078-0432.CCR-19-1895

  • 12

    LuJMcEachernDLiSEllisMJWangS. Reactivation of p53 by MDM2 Inhibitor MI-77301 for the Treatment of Endocrine-Resistant Breast Cancer. Mol Cancer Ther (2016) 15(12):2887–93. doi: 10.1158/1535-7163.MCT-16-0028

  • 13

    DokumcuKSimonianMFarahaniRM. miR4673 Improves Fitness Profile of Neoplastic Cells by Induction of Autophagy. Cell Death Dis (2018) 9(11):1068. doi: 10.1038/s41419-018-1088-6

  • 14

    ChangC-JChaoC-HXiaWYangJ-YXiongYLiC-Wet al. p53 Regulates Epithelial-Mesenchymal Transition and Stem Cell Properties Through Modulating Mirnas. Nat Cell Biol (2011) 13(3):317–23. doi: 10.1038/ncb2173

  • 15

    KimNHKimHSLiX-YLeeIChoiH-SKangSEet al. A P53/miRNA-34 Axis Regulates Snail1-dependent Cancer Cell Epithelial-Mesenchymal Transition. J Cell Biol (2011) 195(3):417–33. doi: 10.1083/jcb.201103097

  • 16

    van NesJGHde KruijfEMPutterHFaratianDMunroACampbell Fet al. Co-Expression of SNAIL and TWIST Determines Prognosis in Estrogen Receptor-Positive Early Breast Cancer Patients. Breast Cancer Res Treat (2012) 133(1):4959. doi: 10.1007/s10549-011-1684-y

  • 17

    LiYZhangJ. Expression of Mutant p53 in Oral Squamous Cell Carcinoma is Correlated With the Effectiveness of Intra-Arterial Chemotherapy. Oncol Lett (2015) 10(5):2883–7. doi: 10.3892/ol.2015.3651

  • 18

    GoldhirschAWinerEPCoatesASGelberRDPiccart-GebhartMThürlimannBet al. Personalizing the Treatment of Women With Early Breast Cancer: Highlights of the St Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2013. Ann Oncol (2013) 24(9):2206–23. doi: 10.1093/annonc/mdt303

  • 19

    ZhangL-YZhangY-QZengY-ZZhuJ-LChenHCWeiX-Let al. TRPC1 Inhibits the Proliferation and Migration of Estrogen Receptor-Positive Breast Cancer and Gives a Better Prognosis by Inhibiting the PI3K/AKT Pathway. Breast Cancer Res Treat (2020) 182(1):2133. doi: 10.1007/s10549-020-05673-8

  • 20

    XuYQinLSunTWuHHeTYangZet al. Twist1 Promotes Breast Cancer Invasion and Metastasis by Silencing Foxa1 Expression. Oncogene (2017) 36(8):1157–66. doi: 10.1038/onc.2016.286

  • 21

    HartsoughEJWeissMBHeilmanSAPurwinTJKugelCH3rdRosenbaumSRet al. CADM1 is a TWIST1-regulated Suppressor of Invasion and Survival. Cell Death Dis (2019) 10(4):281. doi: 10.1038/s41419-019-1515-3

  • 22

    WangSZhongLLiYXiaoDZhangRLiaoDet al. Up-Regulation of PCOLCE by TWIST1 Promotes Metastasis in Osteosarcoma. Theranostics (2019) 9(15):4342–53. doi: 10.7150/thno.34090

  • 23

    FridmanJSLoweSW. Control of Apoptosis by P53. Oncogene (2003) 22(56):9030–40. doi: 10.1038/sj.onc.1207116

  • 24

    SigalARotterV. Oncogenic Mutations of the p53 Tumor Suppressor: The Demons of the Guardian of the Genome. Cancer Res (2000) 60(24):6788–93.

  • 25

    AnsieauSBastidJDoreauAMorelA-PBouchetBPThomasCet al. Induction of EMT by Twist Proteins as a Collateral Effect of Tumor-Promoting Inactivation of Premature Senescence. Cancer Cell (2008) 14(1):7989. doi: 10.1016/j.ccr.2008.06.005

  • 26

    ShiotaMIzumiHOnitsukaTMiyamotoNKashiwagiEKidaniAet al. Twist and p53 Reciprocally Regulate Target Genes Via Direct Interaction. Oncogene (2008) 27(42):5543–53. doi: 10.1038/onc.2008.176

  • 27

    Kogan-SakinITabachYBuganimYMolchadskyASolomonHMadarSet al. Mutant p53(R175H) Upregulates Twist1 Expression and Promotes Epithelial-Mesenchymal Transition in Immortalized Prostate Cells. Cell Death Differ (2011) 18(2):271–81. doi: 10.1038/cdd.2010.94

  • 28

    SchirosiLDe SummaSTommasiSParadisoAGaspariniGPopescuOet al. VEGF and TWIST1 in a 16-Biomarker Immunoprofile Useful for Prognosis of Breast Cancer Patients. Int J Cancer (2017) 141(9):1901–11. doi: 10.1002/ijc.30868

  • 29

    Yang-HartwichYTedjaYYRobertsCMGoodner-BinghamJCardenasCGureaMet al. p53-Pirh2 Complex Promotes Twist1 Degradation and Inhibits EMT. Mol Cancer Res (2019) 17(1):153–64. doi: 10.1158/1541-7786.MCR-18-0238

  • 30

    LiQ-QXuJ-DWangW-JCaoX-XChenQTangFet al. Twist1-Mediated Adriamycin-Induced Epithelial-Mesenchymal Transition Relates to Multidrug Resistance and Invasive Potential in Breast Cancer Cells. Clin Cancer Res (2009) 15(8):2657–65. doi: 10.1158/1078-0432.CCR-08-2372

  • 31

    MaestroRDei TosAPHamamoriYKrasnokutskySSartorelliVKedesLet al. Twist is a Potential Oncogene That Inhibits Apoptosis. Genes Dev (1999) 13(17):2207–17. doi: 10.1101/gad.13.17.2207

  • 32

    PiccininSToninESessaSDemontisSRossiSPecciariniLet al. A “Twist Box” Code of p53 Inactivation: Twist Box: p53 Interaction Promotes p53 Degradation. Cancer Cell (2012) 22(3):404–15. doi: 10.1016/j.ccr.2012.08.003

  • 33

    YangLHouYYuanJTangSZhangHZhuQet al. Twist Promotes Reprogramming of Glucose Metabolism in Breast Cancer Cells Through PI3K/AKT and p53 Signaling Pathways. Oncotarget (2015) 6(28):25755–69. doi: 10.18632/oncotarget.4697

  • 34

    StratiANikolaouMGeorgouliasVLianidouES. Prognostic Significance of TWIST1, CD24, CD44, and ALDH1 Transcript Quantification in EpCAM-Positive Circulating Tumor Cells From Early Stage Breast Cancer Patients. Cells (2019) 8(7):652. doi: 10.3390/cells8070652

  • 35

    ConstantinouCPapadopoulosSKarydaEAlexopoulosAAgnantiNBatistatouA. Expression and Clinical Significance of Claudin-7, PDL-1, PTEN, c-Kit, c-Met, c-Myc, ALK, CK5/6, CK17, p53, EGFR, Ki67, p63 in Triple-Negative Breast Cancer-A Single Centre Prospective Observational Study. In Vivo (2018) 32(2):303–11. doi: 10.21873/invivo.11238

  • 36

    SadighiSZokaasadiMKasaeianAMaghsudiSJahanzadIFumaniHK. The Effect of Immunohistochemically Detected p53 Accumulation in Prognosis of Breast Cancer; A Retrospective Survey of Outcome. PloS One (2017) 12(8):e0182444. doi: 10.1371/journal.pone.0182444

  • 37

    ParralesAIwakumaT. Targeting Oncogenic Mutant p53 for Cancer Therapy. Front Oncol (2015) 5:288. doi: 10.3389/fonc.2015.00288

Summary

Keywords

mutant p53, Twist1, breast cancer, independent prognostic factor, poor prognosis

Citation

Zhang Y-Q, Zhang F, Zeng Y-Z, Chen M, Huang W-H, Wu J-D, Chen W-L, Gao W-L, Bai J-W, Yang R-Q, Zeng H-C, Wei X-L and Zhang G-J (2021) Mutant p53 and Twist1 Co-Expression Predicts Poor Prognosis and Is an Independent Prognostic Factor in Breast Cancer. Front. Oncol. 11:628814. doi: 10.3389/fonc.2021.628814

Received

13 November 2020

Accepted

04 June 2021

Published

24 June 2021

Volume

11 - 2021

Edited by

Yi-Zhou Jiang, Fudan University, China

Reviewed by

Hailin Tang, Sun Yat-sen University Cancer Center (SYSUCC), China; Jin-Ming Yang, University of Kentucky, United States; Na Niu, University of Texas MD Anderson Cancer Center, United States; Donghong Zhang, Georgia State University, United States

Updates

Copyright

*Correspondence: Guo-Jun Zhang, ; Xiao-Long Wei,

†These authors have contributed equally to this work

This article was submitted to Women’s Cancer, a section of the journal Frontiers in Oncology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics