Integrating pathology, chromosomal instability and mutations for risk stratification in early-stage endometrioid endometrial carcinoma
Cell & Bioscience volume 10, Article number: 122 (2020)
Risk stratifications for endometrial carcinoma (EC) depend on histopathology and molecular pathology. Histopathological risk stratification lacks reproducibility, neglects heterogeneity and contributes little to surgical procedures. Existing molecular stratification is useless in patients with specific pathological or molecular characteristics and cannot guide postoperative adjuvant radiotherapies. Chromosomal instability (CIN), the numerical and structural alterations of chromosomes resulting from ongoing errors of chromosome segregation, is an intrinsic biological mechanism for the evolution of different prognostic factors of histopathology and molecular pathology and may be applicable to the risk stratification of EC.
By analyzing CIN25 and CIN70, two reliable gene expression signatures for CIN, we found that EC with unfavorable prognostic factors of histopathology or molecular pathology had serious CIN. However, the POLE mutant, as a favorable prognostic factor, had elevated CIN signatures, and the CTNNB1 mutant, as an unfavorable prognostic factor, had decreased CIN signatures. Only if these two mutations were excluded were CIN signatures strongly prognostic for outcomes in different adjuvant radiotherapy subgroups. Integrating pathology, CIN signatures and POLE/CTNNB1 mutation stratified stageIendometrioid EC into four groups with improved risk prognostication and treatment recommendations.
We revealed the possibility of integrating histopathology and molecular pathology by CIN for risk stratification in early-stage EC. Our integrated risk model deserves further improvement and validation.
Endometrial carcinoma (EC) is the sixth most common malignant tumor in females worldwide and the second most common in the female reproductive system . The risk stratification of EC is the prerequisite for the accurate evaluation of prognosis, and its ultimate goal is to improve the outcome of patients through the optimization of treatment guidelines. There are currently two kinds of stratification systems, conventional pathology assignment in the guidelines and emerging molecular classification proposed by The Cancer Genome Atlas (TCGA) [2, 3].
In the former system, prognostic factors of histopathology, such as histopathological type, grade, stage, myometrial invasion (MI) and lymphovascular space invasion (LVSI), constitute indications for risk assessment and adjuvant radiotherapy . Numerous retro- and prospective clinical studies have demonstrated that the number and severity of prognostic factors of histopathology positively correlate with the risk of recurrence and the extent of adjuvant therapy in EC . Nevertheless, the lack of consensus among pathologists on the histopathological type and tumor grade assignment has resulted in the same woman receiving different classifications, treatments, and clinical outcomes . In addition to this poor reproducibility of prognostic factors, tremendous diversity in clinical outcomes of patients with the same clinicopathological features suggests that the heterogeneity of EC is ignored in this traditional system . Since most of the prognostic factors of histopathology used for risk stratification are only available after surgery, such as MI and LVSI, this risk model contributes little to decisions regarding surgical procedures.
Existing molecular prognostic factors, such as POLE mutation, copy number variation (CNV) and abnormal expression of mismatch repair proteins, classify EC into four molecular subtypes: POLE-mutant, microsatellite instability (MSI), low copy number variation (CNV-L) and high copy number variation (CNV-H) . In addition, CTNNB1 mutation and L1CAM expression are two independent unfavorable prognostic factors [5,6,7]. The accurate and objective detection of all these molecular features makes up for the defects of histopathology mentioned above and improves the risk assessment of EC [5, 7]. However, this prognostic refinement, which only exists in patients categorized as “high-intermediate-risk” by the guidelines , is not conclusive in “high-risk” EC and is utterly ineffective in “low-risk” disease [8, 9]. In addition to being very expensive and complicated, multiplatform and multimolecular detections also generate some “multiple classifiers” that cannot be stratified accurately and reasonably due to the multiple molecular features in the same patient [9, 10]. More importantly, adjuvant radiotherapy recommendations for patients with specific molecular abnormalities still come from guidelines based on histopathology, and no targeted indication can be used as a Ref. . Therefore, both histopathological and existing molecular stratifications have advantages and disadvantages. We envisioned whether there were more suitable biomarkers and strategies to integrate histopathology and molecular pathology in clinical practice.
Chromosomal instability (CIN), which originates from ongoing errors of chromosome segregation and eventually manifests as both numerical and structural aberrations of chromosomes, including aneuploidy, polyploidy, and CNV [11, 12], exists in approximately 60%–80% of tumors . On the one hand, CIN contributes to adverse phenotypes of tumors, including malignant transformation, poor differentiation, invasion, metastasis, immune evasion and treatment resistance [14,15,16,17,18]. On the other hand, it is the end result of a number of molecular processes, such as mutations in DNA checkpoint genes, microtubule spindle defects, telomere dysfunction and even MSI [19,20,21]. As a common hallmark and mechanism underlying different phenotypes and molecular features of tumors, CIN may be a common entry point to explore different prognostic factors of histopathology and molecular pathology in EC. Although the respective roles of chromosomal content and chromatin structure in EC have been associated with histopathology and molecular pathology [22,23,24,25], the overall impact of the numerical and structural aberrations of chromosomes, which is the significance of CIN, is unclear. Since there is no CIN-specific biomarker for EC, we selected the CIN25 and CIN70 signatures from a pan-cancer genomic instability study to measure the CIN status . Based on the top 25 and 70 genes that have correlations with “total functional aneuploidy” in solid tumors, CIN25 and CIN70 signatures have been proven to fully reflect the numerical and structural complexities of chromosomes and have been successfully used in a broad variety of cancer types and research fields [14, 15, 26,27,28]. In the present study, our aims were, first, to investigate the interrelationships between the CIN signature and prognostic factors of histopathology or molecular pathology in EC; and second, relying on the integration of the CIN signature and existing stratification systems, to design a novel risk stratification model for improved prognostic refinement and better management of EC.
Relationships between CIN and prognostic factors of histopathology in EC
To investigate the CIN reflected by CIN signatures in EC, we first confirmed the difference in CIN signatures between benign and malignant endometria. In the TCGA Uterine Corpus Endometrial Carcinoma (UCEC) cohort, 23 cancer samples had notably increased CIN25 and CIN70 expression levels compared to matched adjacent normal tissues (CIN25: p < 0.001; CIN70: p < 0.001; Additional file 1: Figure S1a). Analysis in the GSE63678 dataset, which contained endometrioid EC (EEC) and four rare pathological types (mixed carcinoma with villoglandular, squamous differentiation, clear cell or papillary serous) gave similar results (CIN25: p = 0.003; CIN70: p = 0.003; Additional file 1: Figure S1b). Additionally, in the GSE17025 dataset, ECs had significantly increased CIN25 and CIN70 compared with benign lesions of the endometrium, including polyps and atrophic, inactive or cystic endometria (CIN25: p < 0.001; CIN70: p < 0.001; Additional file 1: Figure S1c).
The nearly identical outcomes of these detections indicated that abnormal chromosomal stability represented by elevated CIN signatures was a dominant feature of EC. For further exploration of CIN in EC, we then compared CIN signatures among prognostic factors of histopathology.
First, in the TCGA UCEC cohort, the highest, intermediate and lowest CIN25 and CIN70 values were found in Grade 3, Grade 2 and Grade 1 patients, respectively (CIN25: p < 0.001; CIN70: p < 0.001; Fig. 1a). Meta-analysis including 483 Grade 1 & 2 and 478 Grade 3 patients from 9 EC datasets confirmed the aggravated CIN in Grade 3 (CIN25 SMD: 0.985, 95% confidence interval (CI): 0.85 to 1.13, p = 0.000; CIN70 SMD: 1.009, 95% CI: 0.87 to 1.15, p = 0.000; Fig. 1b and Additional file 1: Figure S1d). This finding suggested that the more serious the CIN, the poorer the tumor differentiation. Second, we observed low-level expression of CIN signatures in EECs and high-level expression in non-EECs from TCGA (CIN25: p < 0.001; CIN70: p < 0.001; Fig. 1c). That unfavorable histopathological type of EC tended to have severe CIN was further demonstrated by meta-analysis comprising 710 EEC and 253 non-EEC patients from 9 datasets (CIN25 SMD: 0.69, 95% CI 0.54 to 0.84, p = 0.000; CIN70 SMD: 0.63, 95% CI 0.48 to 0.78, p = 0.000; Fig. 1d and Additional file 1: Figure S1e). Third, a meta-analysis of 447 Stage I & II patients versus 184 Stage III & IV patients from 4 EC datasets showed that Stage III & IV patients had obviously increased CIN signatures (CIN25 SMD: 0.602, 95% CI 0.43 to 0.78, p = 0.000; CIN70 SMD: 0.592, 95% CI 0.42 to 0.77, p = 0.000; Fig. 1e and Additional file 1: Figure S1f). Furthermore, patients with longer distances of lymph node metastasis (aortic) or deeper MI (MI > 50%) in TCGA had much higher CIN25 and CIN70 (CIN25: p < 0.05; CIN70: p < 0.05; Fig. 1f, g). Thus, the variation in CIN was also an important characteristic of EC progression. Finally, we detected significantly positive correlations between diagnosis age and CIN signatures in EEC patients with Stage I, Grades 1 & 2 and MI < 50% from TCGA (CIN25: R2 = 0.05, p = 0.010; CIN70: R2 = 0.05, p = 0.011; Fig. 1h). Analysis of the GSE17025 dataset yielded similar results (CIN25: R2 = 0.20, p = 0.050; CIN70: R2 = 0.25, p = 0.025; Fig. 1i). Patients older than 60 tended to have elevated CIN25 and CIN70 compared with younger patients in TCGA (CIN25: p = 0.0050; CIN70: p = 0.0054; Fig. 2j left). This trend between the two age subgroups did not reach a level of statistical significance in GSE17025, possibly due to insufficient samples (Fig. 2j left).
Relationships between CIN and prognostic factors of molecular pathology in EC
As all unfavorable prognostic factors of histopathology are tightly associated with aggravated CIN, we speculated whether CIN signatures could be used to conduct risk assessments for different patients in the same adjuvant radiotherapy subgroup classified by the guidelines (observation (OB) subgroup, vaginal brachytherapy (VBT) subgroup and external beam radiation therapy (EBRT) subgroup; "Materials and methods" section and Table 1), thus providing some opportunities to further optimize indications for postoperative adjuvant therapy. Although patients with a high risk of recurrence or progression tended to have high CIN signatures, the areas under the curve (AUCs) for 5-year disease-free survival (DFS) of the OB, VBT, EBRT and EBRT EEC subgroups were not more than 0.67 (Fig. 2a), which was unsatisfactory and prompted us to investigate possible factors for weakening the predictive power of CIN signatures.
Prognostic factors of molecular pathology became the focus of our investigation. Among the TCGA molecular subtypes of EC except POLE-mutant, CNV-L, MSI and CNV-H had the lowest, intermediate and highest risks of recurrence, respectively, and correspondingly had the lowest, intermediate and highest CIN25 and CIN70 (CIN25: p < 0.001; CIN70: p < 0.001; Fig. 2b, c) [5, 29, 30], which once again implied that CIN might positively correlate with the risk of recurrence in EC. The only exceptional subtype was POLE-mutant, whose prognosis was the best among the four TCGA molecular subtypes, but its CIN signature expression was comparable to that of CNV-H, which had the worst outcome (CIN25: p > 0.05; CIN70: p > 0.05; Fig. 2b, c) [5, 29, 30]. This phenomenon inspired us to explore whether other mutations with prognostic value also had special CIN signatures and in which adjuvant radiotherapy subgroup these special CIN signatures existed. To this end, we compared CIN signatures in wild-type patients with those in POLE, CTNNB1, PTEN, PIK3CA, FGFR2 and PPP2R1A mutant patients from subgroups of OB, VBT, EBRT and ICGC PanCancer Analysis of Whole Genomes (PCAWG) (Fig. 2d and Additional file 2: Figure S2). POLE mutant patients in the OB and VBT subgroups did not relapse or die (Fig. 2e) but had higher expression of CIN25 and CIN70 compared with wild-type patients (CIN25: p < 0.05; CIN70: p < 0.05; Fig. 2d), which might interfere with the risk assessment of CIN signatures. In the OB and EBRT subgroups, the CTNNB1 mutation was another special mutation that had much lower CIN signatures (CIN25: p < 0.05; CIN70: p < 0.05; Fig. 2d and Additional file 2: Figure S2e) but had a much worse prognosis than the wide type (Fig. 2f) [5, 31]. Multivariable Cox models further demonstrated that CTNNB1 mutation was an unfavorable prognostic factor independent of CIN signatures in the OB and EBRT subgroups (Tables 2 and 3). However, this conclusion did not hold in the VBT subgroup, whose CIN signature expression was exactly similar between the CTNNB1 mutant and the wild type (Fig. 2d and Table 2).
CIN signatures were prognostic in different adjuvant radiotherapy subgroups of EC
Consequently, we tested the prognostic value of CIN signatures in different adjuvant radiotherapy subgroups excluding different special mutations. For the OB subgroup without POLE and CTNNB1 mutations, the AUC based on CIN70 was 0.76 (Fig. 3a), and the CIN70 High group predicted worse DFS than the CIN70 Low group (Fig. 3b). For POLE wild types from the VBT subgroup, the AUC based on CIN25 was 0.71 (Fig. 3a), and the CIN25 High group had a much lower 5-year DFS rate than the CIN25 Low group (Fig. 3c). For CTNNB1 wild types from EBRT and EBRT EEC patients, the AUCs based on CIN25 were 0.62 and 0.72 (Fig. 3a), and the outcomes of the CIN25 High group were much worse than those of the CIN25 Low group (Fig. 3d, e). The predictive powers of the Fraction Genome Altered (FGA) and Aneuploidy Score, two signatures that only evaluate chromosomal content, were far inferior to that of CIN signatures (Fig. 3a). Recurrent patients belonging to different histopathological types or TCGA molecular subtypes can be effectively evaluated by CIN signatures in different adjuvant radiotherapy subgroups (Fig. 3f).
Since the AUCs based on CIN70 for DFS and OS of CTNNB1-mutant patients from the OB subgroup were 0.71 and 0.72 (Fig. 3g), we were curious whether CIN could also play a role in the risk assessment of these patients. Although no statistically significant association between the CIN70 Low group and the CIN70 High group was observed, patients with sufficiently long follow-up in the CIN70 High group exhibited a trend toward worse 5-year DFS (Fig. 3h left). We extended our analysis to 10-year OS and found that the outcome of the CIN70 High group was much worse than that of the CIN70 Low group (Fig. 3h right). We therefore reasoned that the CIN signature could and should be used to stratify the CTNNB1-mutant patients from the OB subgroup.
Integrated risk assessment for Stage I EEC from TCGA
According to the different effects of CIN signatures, mutations and pathology, a risk assessment model integrating all these factors is proposed in Fig. 4a for Stage I EEC. In this model, four risk profiles (low, intermediate, high and ultrahigh risk) with different prognoses were considered suitable to receive OB, VBT, EBRT and radiotherapy in combination with systemic therapy after surgery. Among the different existing risk stratification systems, our integrated risk model had the highest AUCs for both DFS and OS (AUC for DFS = 0.75, AUC for OS = 0.76; Fig. 4b) and was the only system that had significant prognostic value for both DFS and OS (Fig. 4c, d; Additional file 3: Figure S3).
Through the comparison and meta-analysis of CIN signatures in multiple EC datasets, our study demonstrated that unfavorable prognostic factors of histopathology and molecular pathology, including poor differentiation, non-EEC, advanced disease, deep MI, advanced age, MSI and CNV-H, usually aggravated CIN. A favorable prognostic factor, POLE mutation, and an unfavorable prognostic factor, CTNNB1 mutation, did not follow the above trend. The prognostic value of CIN signatures was well established in different adjuvant radiotherapy subgroups without POLE/CTNNB1 mutations and in CTNNB1 mutant patients from the OB subgroup. An integrated risk model that combines pathology, CIN signatures and mutations was defined for improved prognostic refinement and better management of Stage I EEC.
Most non-EECs are serous and high-grade cancers that exactly have complex aneuploidies and polyploidy ; hence, CIN showed consistent changes in fields of histopathological type and tumor differentiation of EC (Fig. 1a–d). At least three potential mechanisms generated by CIN, including the induction of mesenchymal transition, the activation of the STING pathway and immune evasion, may contribute to invasion and metastasis , which may explain the high CIN25 and CIN70 in Stage III & IV patients and in patients with deep MI or aortic lymph node metastasis (Fig. 1e–g). Although we cannot verify the CIN status in LVSI-positive patients due to a lack of sufficient pathological information, we speculate that CIN may also increase in LVSI-positive cases since aneuploidy has been correlated with the LVSI of EC . Given the propensity for aging somatic cells to generate unstable chromosomes resulting from gene misexpression, telomeric attrition and senescence failure [33,34,35], older EC patients were more prone to CIN enrichment (Fig. 1h–j).
Several well-recognized molecular features of EC also have characteristic CIN. One of the final results of CIN is CNV . Therefore, we observed the lowest CIN signature expression in CNV-L and the highest expression in CNV-H (Fig. 2c). The fact that MSI causes some degree of genomic instability and the tendency for MSI to have aggressive phenotypes are two possible reasons for the moderate exacerbation of CIN in MSI patients [18,19,20, 29, 30]. From CNV-L to MSI and then to CNV-H, as the CIN gradually becomes serious, the risk of recurrence gradually increases (Fig. 2b). In terms of the MSI subtype itself, high CIN signatures were unfavorable prognostic factors . These two pieces of evidence, combined with the fact that CIN signatures did identify recurrent patients who belonged to different TCGA molecular subtypes in each adjuvant radiotherapy subgroup (Fig. 3b–f), implies that CNV-L, MSI, and CNV-H may be pooled together for prognosis evaluation by CIN.
Mutation of POLE causes impaired proofreading activity and DNA repair ability, followed by poor fidelity of DNA replication and severe genomic instability [36, 37]. This makes the CIN of the POLE-mutant subtype roughly the same as that of CNV-H (Fig. 2c). Unlike POLE mutation, however, why the mutation of CTNNB1 is associated with a more stable chromosome status is not clear (Fig. 2d). The aberrant WNT/CTNNB1 pathway in colon cancer always induces CIN [38, 39], so the complete opposite relationship between CTNNB1 mutation and CIN in EC is confusing and interesting. Considering that patients with unstable chromosomes usually have poor clinical outcomes , how aggravated CIN produces an excellent prognosis in POLE-mutant patients and how alleviated CIN leads to poor outcomes in CTNNB1-mutant patients is another important issue worthy of further research (Fig. 2d–f; Tables 2 and 3). Serious CIN allows tumors to have different clonal selections in response to various biological stimuli and environmental stresses. However, this selective advantage also has a fitness cost for CIN because the extremely excessive instability of chromosomes is not conducive to the stable survival of the tumor cell itself [11, 14, 27, 40, 41]. For this reason, in addition to the immune activation triggered by POLE-related mutations , severe CIN may contribute to the excellent prognosis of POLE-mutant cases. Similarly, CTNNB1-mutant cases, which benefit from the progression and proliferation caused by the activation of WNT/CTNNB1 signaling [6, 43], may protect cells from the adverse effects of this pathway activation with the help of the alleviated CIN. Although this conjecture is still to be confirmed by molecular biology, it may provide CIN-targeted therapeutic strategies for mutation-specific EC.
Based on these data and references, the inherent biological connections between CIN and different prognostic factors of EC suggest that CIN may be a common hallmark in the evolution of different clinicopathological and molecular features, which is the root cause for the success of our integrated risk model (Figs. 3 and 4). From the perspective of risk assessment, the CIN signature, on the one hand, properly addressed the problems of heterogeneity and reproducibility in the conventional pathology system by the precise quantification of CIN status, thereby achieving prognostic refinement. “Multiple classifiers” that cannot be stratified by TCGA subtypes can also obtain accurate and reasonable risk assessments. On the other hand, the prognostic refinement achieved by CIN signatures existed in all adjuvant radiotherapy subgroups in the guidelines, which means that CIN may have more universal applications compared to other risk stratification systems such as TCGA subtypes, FGA and Aneuploidy Score. From a therapeutic point of view, high concordance of molecular alterations between curettage samples and hysterectomy specimens from EC suggested the potential for CIN signatures to guide surgical management [24, 44]. More importantly, because the accurate risk stratification accomplished by the CIN signature presupposed the adjuvant radiotherapy classification based on the guidelines, the treatment recommendations obtained from our integrated risk model may be an intact inheritance of and effective supplement to the indications for postoperative radiotherapy in the guidelines. In summary, the intrinsic relationships between CIN and clinicopathological or molecular features make CIN a bridge that comprehensively integrates histopathology and molecular pathology, which is difficult for other biomarkers to achieve.
To further refine our integrated risk model, we face two outstanding challenges. First, the CIN signature should be optimized on the basis of CIN25 and CIN70. Different adjuvant radiotherapy subgroups in the guidelines have different clinicopathological and molecular features (Table 1), which leads to the same CIN signature having different capabilities of risk assessment in different adjuvant radiotherapy subgroups. For the same adjuvant radiotherapy subgroup, different CIN signatures also have different risk assessment capabilities. Consequently, we believe that it is necessary to improve the respective CIN signatures for the OB, OB & CTNNB1 mutation, VBT and EBRT EEC patients classified by the guidelines to realize the full potential of CIN. Second, the relationships between CIN and LVSI or several EC prognostic factors assessed by immunohistochemistry , such as L1CAM, ER and PR, remain to be explored. It is unclear whether these features are still independent prognostic factors in our integrated model. We look forward to high-quality retrospective studies with mature long-term follow-up data and large sample sizes that will meet these two challenges and provide a solid foundation for future clinical applications.
Overall, except for POLE and CTNNB1 mutations, serious CIN represented by increased CIN25 and CIN70 are characteristic of unfavorable prognostic factors in EC. Integration of pathology, CIN signatures and mutation of POLE/CTNNB1 in Stage I EEC leads to improved prognostic refinement with potential clinical utility. Our integrated risk model holds promise to reduce both overtreatment and undertreatment and deserves further validation and improvement.
Materials and methods
Clinical information, gene expression (Z-score), mutation profiles, FGA and Aneuploidy Score of the TCGA UCEC cohort are available at cBioPortal for Cancer Genomics . RNAseq data (FPKM-UQ) of EC samples and matched adjacent normal tissues were downloaded from the TCGA data portal . Mutation profiles and RNAseq data (FPKM-UQ) of EEC samples in the ICGC PCAWG project were downloaded from the ICGC data portal [48, 49]. E-MTAB-1358 and E-MTAB-5018 were downloaded from the ArrayExpress database [50,51,52]. GSE2109, GSE23518, GSE56026, GSE33723, GSE32507, GSE17025, GSE24537, GSE21882 and GSE63678 were downloaded from the Gene Expression Omnibus (GEO) database [53,54,55,56,57,58,59,60,61]. Each dataset downloaded from ArrayExpress and GEO databases was standardized using the Z-score transformation before calculations of CIN signatures. CIN25 and CIN70 signatures for each sample were the average expressions of the 25 and 70 genes identified by Cart et al. (Additional file 4: Table S1)  and were compared among prognostic factors of histopathology or molecular pathology. The clinicopathological information for each dataset is shown in the meta-analysis.
Adjuvant radiotherapy classification for StageI patients in the TCGA UCEC cohort
There were three adjuvant therapeutic strategies after surgery for stage I EC patients, namely, observation (OB), vaginal brachytherapy (VBT), and external beam radiation therapy (EBRT). Indications for the three adjuvant radiotherapies in the guidelines of ESMO-ESGO-ESTRO were based on six established clinicopathological risk factors, including age, histologic type, grade, stage, MI, and LVSI . LVSI was missing in the TCGA UCEC cohort; therefore, we had to conduct the classification with the other five risk factors. The OB subgroup in the guidelines was defined as a) stage IA EEC with Grades 1 & 2 and b) stage IB EEC with Grades 1 & 2 and less than 60 years old. Patients in the EBRT subgroup followed the following criteria: a) stage IB EEC with Grade 3; b) stage I non-EEC. The VBT subgroup consisted of the remaining patients, including a) stage IB EEC with Grades 1 & 2 and age > 60; b) stage IA EEC with Grade 3. Patients who did not have complete or accurate information for classification and survival analysis, who had other malignancies or who had positive surgical margins were excluded. Ultimately, there were 123, 92 and 79 patients in the OB, VBT and EBRT subgroups, respectively. Detailed information is presented in Table 1.
Survival analysis for different adjuvant radiotherapy subgroups
In the OB, VBT, EBRT and EBRT EEC subgroups, AUC and optimal cutoff values based on CIN25 and CIN70 signatures, FGA and Aneuploidy Score were determined by the time-dependent receiver operating curve using the “survivalROC” package on the R platform. Kaplan–Meier curves and log-rank tests were carried out to predict 5-year DFS and 10-year OS based on the optimal cutoff values or mutation status of different subgroups. Cox proportional hazards models were used to evaluate the prognostic value of mutations and CIN signatures. Covariates violating the proportional hazards assumption were added as time-dependent covariates in the Cox regression models.
If not specified otherwise, comparisons among clinicopathological features and among molecular alterations were tested using the Chi square test for categorical variables and using the Mann–Whitney test or one-way analysis of variance (ANOVA) in continuous cases. SPSS 21.0 and GraphPad Prism 8 software programs were used to perform statistical and survival analyses and to plot graphs. Meta-analysis was performed using STATA 12.0. Data in this article are presented as the mean ± S.D., and p values were based on two-sided tests with < 0.05 considered statistically significant.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its supplementary information files.
The Cancer Genome Atlas
Lymphovascular space invasion
Copy number variation
Low copy number variation
High copy number variation
Uterine Corpus Endometrial Carcinoma
Endometrioid endometrial carcinoma
External beam radiation therapy
area under the curve
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Colombo N, Creutzberg C, Amant F, Bosse T, Gonzalez-Martin A, Ledermann J, et al. ESMO-ESGO-ESTRO Consensus Conference on Endometrial Cancer: diagnosis, treatment and follow-up. Ann Oncol. 2016;27(1):16-41.
Cancer Genome Atlas Research N, Kandoth C, Schultz N, Cherniack AD, Akbani R, Liu Y, et al. Integrated genomic characterization of endometrial carcinoma. Nature. 2013. 497(7447):67-73.
Clarke BA, Gilks CB. Endometrial carcinoma: controversies in histopathological assessment of grade and tumour cell type. J Clin Pathol. 2010;63(5):410–5.
Stelloo E, Nout RA, Osse EM, Jurgenliemk-Schulz IJ, Jobsen JJ, Lutgens LC, et al. Improved risk assessment by integrating molecular and clinicopathological factors in early-stage endometrial cancer-combined analysis of the PORTEC cohorts. Clin Cancer Res. 2016;22(16):4215–24.
Liu Y, Patel L, Mills GB, Lu KH, Sood AK, Ding L, et al. Clinical significance of CTNNB1 mutation and Wnt pathway activation in endometrioid endometrial carcinoma. J Natl Cancer Inst. 2014. 106(9).
Zeimet AG, Reimer D, Huszar M, Winterhoff B, Puistola U, Azim SA, et al. L1CAM in early-stage type I endometrial cancer: results of a large multicenter evaluation. J Natl Cancer Inst. 2013;105(15):1142–50.
Stelloo E, Bosse T, Nout RA, MacKay HJ, Church DN, Nijman HW, et al. Refining prognosis and identifying targetable pathways for high-risk endometrial cancer; a TransPORTEC initiative. Mod Pathol. 2015;28(6):836–44.
McAlpine J, Leon-Castillo A, Bosse T. The rise of a novel classification system for endometrial carcinoma; integration of molecular subclasses. J Pathol. 2018;244(5):538–49.
Schultheis AM, Martelotto LG, De Filippo MR, Piscuglio S, Ng CK, Hussein YR, et al. TP53 Mutational Spectrum in Endometrioid and Serous Endometrial Cancers. Int J Gynecol Pathol. 2016;35(4):289–300.
Turajlic S, Sottoriva A, Graham T, Swanton C. Resolving genetic heterogeneity in cancer. Nat Rev Genet. 2019;20(7):404–16.
Sansregret L, Vanhaesebroeck B, Swanton C. Determinants and clinical implications of chromosomal instability in cancer. Nat Rev Clin Oncol. 2018;15(3):139–50.
Carter SL, Cibulskis K, Helman E, McKenna A, Shen H, Zack T, et al. Absolute quantification of somatic DNA alterations in human cancer. Nat Biotechnol. 2012;30(5):413–21.
Endesfelder D, McGranahan N, Birkbak NJ, Szallasi Z, Kschischo M, Graham TA, et al. A breast cancer meta-analysis of two expression measures of chromosomal instability reveals a relationship with younger age at diagnosis and high risk histopathological variables. Oncotarget. 2011;2(7):529–37.
Swanton C, Nicke B, Schuett M, Eklund AC, Ng C, Li Q, et al. Chromosomal instability determines taxane response. Proc Natl Acad Sci USA. 2009;106(21):8671–6.
Turajlic S, Swanton C. Metastasis as an evolutionary process. Science. 2016;352(6282):169–75.
Bakhoum SF, Cantley LC. The multifaceted role of chromosomal instability in cancer and its microenvironment. Cell. 2018;174(6):1347–60.
Rajagopalan H, Nowak MA, Vogelstein B, Lengauer C. The significance of unstable chromosomes in colorectal cancer. Nat Rev Cancer. 2003;3(9):695–701.
Li LS, Kim NG, Kim SH, Park C, Kim H, Kang HJ, et al. Chromosomal imbalances in the colorectal carcinomas with microsatellite instability. Am J Pathol. 2003;163(4):1429–36.
Trautmann K, Terdiman JP, French AJ, Roydasgupta R, Sein N, Kakar S, et al. Chromosomal instability in microsatellite-unstable and stable colon cancer. Clin Cancer Res. 2006;12(21):6379–85.
Grady WM. Genomic instability and colon cancer. Cancer Metastasis Rev. 2004;23(1–2):11–27.
Proctor L, Pradhan M, Leung S, Cheng A, Lee CH, Soslow RA, et al. Assessment of DNA Ploidy in the ProMisE molecular subgroups of endometrial cancer. Gynecol Oncol. 2017;146(3):596–602.
Hveem TS, Njolstad TS, Nielsen B, Syvertsen RA, Nesheim JA, Kjaereng ML, et al. Changes in chromatin structure in curettage specimens identifies high-risk patients in endometrial cancer. Cancer Epidemiol Biomarkers Prev. 2017;26(1):61–7.
Njolstad TS, Trovik J, Hveem TS, Kjaereng ML, Kildal W, Pradhan M, et al. DNA ploidy in curettage specimens identifies high-risk patients and lymph node metastasis in endometrial cancer. Br J Cancer. 2015;112(10):1656–64.
Song T, Lee JW, Choi CH, Kim TJ, Bae DS, Sung CO, et al. Ploidy and S-phase fraction are correlated with lymphovascular space invasion that is predictive of outcomes in endometrial cancer. Int J Clin Oncol. 2012;17(6):590–7.
Carter SL, Eklund AC, Kohane IS, Harris LN, Szallasi Z. A signature of chromosomal instability inferred from gene expression profiles predicts clinical outcome in multiple human cancers. Nat Genet. 2006;38(9):1043–8.
Birkbak NJ, Eklund AC, Li Q, McClelland SE, Endesfelder D, Tan P, et al. Paradoxical relationship between chromosomal instability and survival outcome in cancer. Cancer Res. 2011;71(10):3447–52.
Decaux O, Lode L, Magrangeas F, Charbonnel C, Gouraud W, Jezequel P, et al. Prediction of survival in multiple myeloma based on gene expression profiles reveals cell cycle and chromosomal instability signatures in high-risk patients and hyperdiploid signatures in low-risk patients: a study of the Intergroupe Francophone du Myelome. J Clin Oncol. 2008;26(29):4798–805.
Talhouk A, McConechy MK, Leung S, Yang W, Lum A, Senz J, et al. Confirmation of ProMisE: A simple, genomics-based clinical classifier for endometrial cancer. Cancer. 2017;123(5):802–13.
Kommoss S, McConechy MK, Kommoss F, Leung S, Bunz A, Magrill J, et al. Final validation of the ProMisE molecular classifier for endometrial carcinoma in a large population-based case series. Ann Oncol. 2018;29(5):1180–8.
Wortman BG, Bosse T, Nout RA, Lutgens L, van der Steen-Banasik EM, Westerveld H, et al. Molecular-integrated risk profile to determine adjuvant radiotherapy in endometrial cancer: evaluation of the pilot phase of the PORTEC-4a trial. Gynecol Oncol. 2018;151(1):69–75.
Taylor AM, Shih J, Ha G, Gao GF, Zhang X, Berger AC, et al. Genomic and Functional Approaches to Understanding Cancer Aneuploidy. Cancer Cell. 2018. 33(4):676-89 e3.
Geigl JB, Langer S, Barwisch S, Pfleghaar K, Lederer G, Speicher MR. Analysis of gene expression patterns and chromosomal changes associated with aging. Cancer Res. 2004;64(23):8550–7.
Hackett JA, Feldser DM, Greider CW. Telomere dysfunction increases mutation rate and genomic instability. Cell. 2001;106(3):275–86.
Artandi SE, Chang S, Lee SL, Alson S, Gottlieb GJ, Chin L, et al. Telomere dysfunction promotes non-reciprocal translocations and epithelial cancers in mice. Nature. 2000;406(6796):641–5.
Henninger EE, Pursell ZF. DNA polymerase epsilon and its roles in genome stability. IUBMB Life. 2014;66(5):339–51.
Tubbs A, Nussenzweig A. Endogenous DNA damage as a source of genomic instability in cancer. Cell. 2017;168(4):644–56.
Hadjihannas MV, Bruckner M, Jerchow B, Birchmeier W, Dietmaier W, Behrens J. Aberrant Wnt/beta-catenin signaling can induce chromosomal instability in colon cancer. Proc Natl Acad Sci USA. 2006;103(28):10747–52.
Hadjihannas MV, Behrens J. CIN By WNT: growth pathways, mitotic control and chromosomal instability in cancer. Cell Cycle. 2006;5(18):2077–81.
Jamal-Hanjani M, A’Hern R, Birkbak NJ, Gorman P, Gronroos E, Ngang S, et al. Extreme chromosomal instability forecasts improved outcome in ER-negative breast cancer: a prospective validation cohort study from the TACT trial. Ann Oncol. 2015;26(7):1340–6.
Roylance R, Endesfelder D, Gorman P, Burrell RA, Sander J, Tomlinson I, et al. Relationship of extreme chromosomal instability with long-term survival in a retrospective analysis of primary breast cancer. Cancer Epidemiol Biomarkers Prev. 2011;20(10):2183–94.
van Gool IC, Bosse T, Church DN. POLE proofreading mutation, immune response and prognosis in endometrial cancer. Oncoimmunology. 2016;5(3):e1072675.
Klaus A, Birchmeier W. Wnt signalling and its impact on development and cancer. Nat Rev Cancer. 2008;8(5):387–98.
Stelloo E, Nout RA, Naves LC, Ter Haar NT, Creutzberg CL, Smit VT, et al. High concordance of molecular tumor alterations between pre-operative curettage and hysterectomy specimens in patients with endometrial carcinoma. Gynecol Oncol. 2014;133(2):197–204.
Karnezis AN, Leung S, Magrill J, McConechy MK, Yang W, Chow C, et al. Evaluation of endometrial carcinoma prognostic immunohistochemistry markers in the context of molecular classification. J Pathol Clin Res. 2017;3(4):279–93.
cBioPortal for Cancer Genomics Available from: https://www.cbioportal.org/.
TCGA data portal Available from: www.portal.gdc.cancer.gov/.
Consortium ITP-CAoWG. Pan-cancer analysis of whole genomes. Nature. 2020. 578(7793):82-93.
ICGC data portal https://dcc.icgc.org/.
Krakstad C, Birkeland E, Seidel D, Kusonmano K, Petersen K, Mjos S, et al. High-throughput mutation profiling of primary and metastatic endometrial cancers identifies KRAS, FGFR2 and PIK3CA to be frequently mutated. PLoS ONE. 2012;7(12):e52795.
Kusonmano K, Halle MK, Wik E, Hoivik EA, Krakstad C, Mauland KK, et al. Identification of highly connected and differentially expressed gene subnetworks in metastasizing endometrial cancer. PLoS ONE. 2018;13(11):e0206665.
Chiyoda T, Tsuda H, Tanaka H, Kataoka F, Nomura H, Nishimura S, et al. Expression profiles of carcinosarcoma of the uterine corpus-are these similar to carcinoma or sarcoma? Genes Chromosomes Cancer. 2012;51(3):229–39.
Risinger JI, Allard J, Chandran U, Day R, Chandramouli GV, Miller C, et al. Gene expression analysis of early stage endometrial cancers reveals unique transcripts associated with grade and histology but not depth of invasion. Front Oncol. 2013;3:139.
Mhawech-Fauceglia P, Wang D, Kesterson J, Syriac S, Clark K, Frederick PJ, et al. Gene expression profiles in stage I uterine serous carcinoma in comparison to grade 3 and grade 1 stage I endometrioid adenocarcinoma. PLoS ONE. 2011;6(3):e18066.
Levan K, Partheen K, Osterberg L, Olsson B, Delle U, Eklind S, et al. Identification of a gene expression signature for survival prediction in type I endometrial carcinoma. Gene Expr. 2010;14(6):361–70.
Mhawech-Fauceglia P, Wang D, Kesterson J, Clark K, Monhollen L, Odunsi K, et al. Microarray analysis reveals distinct gene expression profiles among different tumor histology, stage and disease outcomes in endometrial adenocarcinoma. PLoS ONE. 2010;5(11):e15415.
Romero-Perez L, Castilla MA, Lopez-Garcia MA, Diaz-Martin J, Biscuola M, Ramiro-Fuentes S, et al. Molecular events in endometrial carcinosarcomas and the role of high mobility group AT-hook 2 in endometrial carcinogenesis. Hum Pathol. 2013;44(2):244–54.
Pappa KI, Polyzos A, Jacob-Hirsch J, Amariglio N, Vlachos GD, Loutradis D, et al. Profiling of discrete gynecological cancers reveals novel transcriptional modules and common features shared by other cancer types and embryonic stem cells. PLoS ONE. 2015;10(11):e0142229.
Kharma B, Baba T, Matsumura N, Kang HS, Hamanishi J, Murakami R, et al. STAT1 drives tumor progression in serous papillary endometrial cancer. Cancer Res. 2014;74(22):6519–30.
Gene Expression Omnibus Available from: https://www.ncbi.nlm.nih.gov/gds.
This research was funded by the National Key Research and Development Program, Grant Number 2016YFC1303012 and was funded by the Nature and Science Foundation of China, Grant Number 81402163, 81402164, 81572569, 81874106.
Ethics approval and consent to participate
Consent for publication
All authors consent for publication.
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Comparison of CIN signatures among prognostic factors of histopathology. a Paired scatter plots of CIN25 (left) and CIN70 (right) expression profiles comparing endometrial carcinomas (red) with paired adjacent normal tissues (blue) from TCGA dataset (n = 23). Each pair of dots indicates the amount of CIN25 or CIN70 expression for a particular patient. P values represent paired t test calculations. b Boxplots of CIN25 (left) and CIN70 (right) expression in normal endometria (blue) and endometrial carcinomas (red). Each dot indicates the amount of CIN25 or CIN70 expression for a particular sample in the GSE63678 dataset. (c) Same as (b) but comparing benign lesions (blue) with endometrial carcinomas (red) from the GSE17025 dataset. P values presented in (b) and (c) are Mann–Whitney test calculations. d Forest plot comparing CIN70 expression in Grade 1 & 2 versus Grade 3 patients. e Forest plot comparing CIN70 expression in EEC versus non-EEC samples. f Forest plot comparing CIN70 expression in Stage I & II versus Stage III & IV patients. In (d), (e) and (f), an inverse variance (IV) fixed effects method was used to meta-analyze the data; squares (blue) represent standardized mean difference (SMD); square size is proportional to weights used in the analysis; bars (gray) represent 95% confidence intervals (CI); diamonds (yellow) represent overall SMD with associated 95% CI (lateral tips). P values: *p < 0.05, **p < 0.01, ***p < 0.001.
Relationships between CIN and molecular prognostic factors. a Boxplot of CIN25 and CIN70 expression in PTEN wild-type (blue) versus mutant patients (red) from the OB, VBT and EBRT groups. b Boxplot of CIN25 and CIN70 expression in PIK3CA wild-type (blue) versus mutant patients (red) from the OB, VBT and EBRT groups. c Boxplot of CIN25 and CIN70 expression in FGFR2 wild-type (blue) versus mutant patients (red) from the OB, VBT and EBRT groups. d Boxplot of CIN25 and CIN70 expression in PPP2R1A wild-type (blue) versus mutant patients (red) from the OB, VBT and EBRT groups. e Boxplot of CIN25 and CIN70 expression in wild-type POLE, CTNNB1, PTEN, PIK3CA, FGFR2 and PPP2R1A (blue) versus mutant patients (red) from the PCAWG EEC cohort. P values in (a)–(e) are Mann–Whitney test calculations (*p < 0.05, **p < 0.01, ***p < 0.001, n.s. not significantly different).
Kaplan–Meier plot for DFS and OS in Stage I EEC from TCGA based on Guidelines, FGA, Aneuploidy Score and TCGA subtypes. a 5-year PFS. b 10-year OS. Patients are grouped by quartiles of FGA and Aneuploidy Score. P values for (a) and (b) are calculations of the log-rank test.
Genes in CIN25 and CIN70.
About this article
Cite this article
Li, Y., Li, J., Guo, E. et al. Integrating pathology, chromosomal instability and mutations for risk stratification in early-stage endometrioid endometrial carcinoma. Cell Biosci 10, 122 (2020). https://doi.org/10.1186/s13578-020-00486-0