Abstract
Objective:
To study and predict the risk of ovarian metastasis (OM) in patients with cervical adenocarcinoma (ADC).
Methods:
Patients with ADC who received surgical treatment from January 2015 to December 2022 in the Obstetrics and Gynecology Hospital of Fudan University were included in the study. Patients were further divided into OP (ovaries were preserved in surgery) and BSO (bilateral salpingo-oophorectomy) groups. For the patients who accepted BSO, 60% of the patients were randomly grouped into a training cohort, and predictive prognostic models were constructed with 10-fold cross-validation through random forest, LASSO, stepwise, and optimum subset analysis. The model with the highest area under receiver operator curve (AUC) was screened out in the testing cohort. The nomogram and its calibration curve presented the possibility of OM in future patients. The prognoses between the different populations were compared using Kaplan–Meier analysis. All data processing was carried out in R 4.2.0 software.
Results:
A total of 934 patients were enrolled in our cohort; 266 patients had their ovaries preserved and 668 patients had BSO according to the previous criteria reported The clinical safety with these criteria was secured, while the 5-year overall survival had no significant difference between the BSO and OP groups (p = 0.069), which suggested that the current criteria could be extended and are more precise. Four predictive models for ovarian metastasis by machine learning were constructed in our study, and the random forest model that obtained the highest AUC in both training and testing sets (0.971 for training and 0.962 for testing set) was taken as the best model. The optimal cut-off point of the ROC curve (specificity 99.5% and 90% sensitivity) was utilized to stratify the patients into high- and low-risk OM. Further comparing the survival curves of the high and low-OM risk groups, it was found that both DFS and OS were significantly prolonged in the low-risk group (p < 0.01). On the basis of this random forest model, a nomogram was used to calculate the OM risk, and the results were validated with calibration. The predictive model was further applied to the whole cohort (934 patients), and we identified the OM low-risk population (n = 822) and the patients with high risk who should be recommended for BSO (n = 112). No significant difference was found in the 5-year survival between the low-risk group with our model and the patients who already had ovaries preserved according to the previous criteria (n = 266).
Conclusion:
The predictive model constructed in our study could identify the low-risk population of OM in patients with ADC, which remarkably extended the number with the previous criteria, for whom we could potentially preserve ovaries to help improve their life quality.
Introduction
Although screening with cytology and human papillomavirus (HPV) have undoubtedly led to a remarkable decline in cervical cancer burden, cervical cancer remains the fourth most common cancer among women worldwide. In young female of reproductive age, cervical cancer is consistently the second leading cause of cancer death in women aged 20–39 years according to Cancer Statistics, 2023 (). The emergence of HPV vaccination could virtually eliminate HPV-related cervical cancer, most of which were squamous cervical cancer, and cervical adenocarcinoma (ADC) accounted for approximately 25% of all cervical cancer cases showing a significant increase, particularly among young women as to the rapidly growing use of HPV vaccine (–). ADC patients have considerable life expectancy as the 5-year overall survival (OS) of early-stage ADC is above 80% (). Thus, ovarian preservation (OP) during surgery seems crucial for younger patients to conserve ovarian endocrine function and fertility. However, considering the medical safety, OP during surgery has been always controversial for ADC patients. Some previous studies recommended that ADC patients should be routinely ovariectomized as the incidence of OM in ADC was reported to be above 10% (, ) Yet, more evidence arose to support OP in selective ADC patients (–) and the selection criteria for OP reported in 2015 (Table 1) (). OP surgery for ADC patients has been widely performed in accordance with the criteria in our center since 2015. Subsequent pathological analysis of the removed ovaries from patients who did not meet the criteria revealed that many of them have no OM, which might indicate that the criteria should be modified, and a predictive model is needed to extend the range involving more ADC patients to preserve ovaries.
Table 1
| Preoperative patients’ characteristics |
|---|
| –Age ≤45 years |
| –Patients who desire to retain ovarian function |
| –No familial predisposition to ovarian cancer |
| Tumor characteristics based on clinical and imaging evidence |
| –FIGO stage ≤IB2 |
| –Size of the tumor ≤4 cm |
| –No parametrial invasion |
| –No corpus invasion |
| –No deep stromal invasion |
| –No evidence of lymph node metastasis (CT scan or PET scan) |
| –No lymphatic vascular space invasion (LVSI) |
| Intraoperative findings |
| –No extra-uterine spread |
| –No lymph node metastasis |
| –Normal appearance of ovaries |
Detailed selection criteria for OP reported ().
However, studies on OP in cervical adenocarcinoma since 2015 have focused exclusively upon the clinical safety and potential high-risk factors causing ovarian metastasis, while no new criteria or predictive models were constructed to our knowledge (, –). Therefore, we collected and studied the characteristics of ADC patients treated in the past 7 years in our center to build and validate a reliable model for OM in cervical adenocarcinoma, subsequently to more accurately preserve ovaries in ADC patients and largely improve their life quality. This was the first time a model for predicting ovarian metastasis of cervical adenocarcinoma was constructed using the machine learning method, which can be practically utilized in real clinic scenario.
Materials and methods
Study setting and design
This was a retrospective study conducted from January 2015 to December 2022 in the Obstetrics and Gynecology Hospital of Fudan University in Shanghai, China, which is the largest medical center for cervical cancer diagnosis and treatment in China with our patients covering 34 states all over China.
Study population
Cervical cancer patients at early stage [stages IA–IIA, restaged according to FIGO 2019 ()], surgically treated between January 2015 and December 2022, were reviewed with informed consent. The study was approved by the medical ethics committee of the Obstetrics and Gynecology Hospital of Fudan University (Ob&Gyn Hospital) (No. 2022-123) on 31 August 2022 and conducted in accordance with the Declaration of Helsinki. Patients with stage IA1 disease were treated through hysterectomy. Modified radical hysterectomy () was performed on patients with FIGO stage IA1 coupled with positive margin or LVSI and FIGO stage IA2. Radical hysterectomy was performed on patients with FIGO stages IB to IIA2 cervical cancer. Pelvic lymphadenectomy ± para-aortic lymphadenectomy was performed on all except for those with FIGO stage IA1 disease. Fallopian tubes were routinely removed and whether or not to remove the ovaries was decided according to the criteria shown in Table 1.
Explanatory and outcome variables
Both clinical and pathological data were obtained from medical records. The explanatory clinicopathological variables included patient age, pre- and post-surgery pathological diagnosis (WHO 2014), tumor size, lymph node metastasis (LNM), location of LNM, depth of stromal invasion, LVSI status, parametrial invasion, uterine corpus invasion (UCI), surgical margin, vaginal involvement, fallopian tube metastasis, OM, FIGO stage (2019), and treatment modalities. The outcome variables included OM according to post-surgery pathology, as well as overall survival (OS) and disease-free survival (DFS) measured from the date of surgery to recurrence/death or the last time patients were censored at follow-up. The end point of follow-up was December 2023, and the date of recurrence/death was obtained from our hospital’s follow-up records.
Descriptive statistics and variable selection
For continuous covariates, we employed receiver operating characteristic (ROC) curves to determine cut-off points for classifying patients into high- and low-risk groups. This approach, facilitated by the “ROCR” R package, allowed us to identify thresholds that optimally discriminated between patients with and without ovarian metastasis. Categorical variables were compared using Chi-squared and Fisher’s exact tests, implemented through the “tableone” R package, to ensure that our baseline characteristics were balanced across risk groups.
Model construction, selection, and calibration
The 668 patients who underwent bilateral salpingo-oophorectomy (BSO) were randomly divided into training and testing cohorts (including 401 and 267 patients, respectively) using “caret” R package. To identify potentially significant factors affecting ovarian metastasis risks, we performed univariate logistic regression analyses on the training cohort. Factors with p-values < 0.1 were considered candidates for inclusion in the multivariate models. Then, we employed multiple methods to construct risk scores in the training cohort, including Least Absolute Shrinkage and Selection Operator (LASSO) regression, stepwise multivariate logistic analysis, optimum subset logistic regression, and random survival forest analysis. Each method was chosen for its unique strengths in handling high-dimensional data, avoiding overfitting, and capturing complex interactions. The method that yielded the highest area under the curve (AUC) of the ROC in the testing cohort was selected as the final model for ovarian metastasis risk prediction. The cutoff point determined by the ROC was used to evaluate whether our patients could preserve ovaries or not in our predictive model. After that, the nomogram and related calibration curves were established based on a new model to present the probability of ovarian metastasis in future patients.
Software and tools
All data processing and analyses were conducted using R 4.2.0 software leveraging a range of specialized packages including “glmnet,” “caret,” “randomForest,” “My.stepwise,” “forestploter,” “forestplot,” “bestglm,” “leaps,” “genefilter,” “Hmisc,” “ISLR,” “rms,” “regplot,” and “ROCR.” These packages offer powerful and flexible tools for statistical modeling, visualization, and interpretation enabling us to conduct our analyses with precision and efficiency.
Results
Patient characteristics
A total of 934 patients who were screened with inclusion criteria were finally enrolled in our study (Figure 1), with an average age of 46.32 ± 9.53 years (median: 45.50 years; range: 22–74 years), among whom BSO was performed in 668 cases and OP in 266 patients based on the criteria (Table 1) (). Patients who had OP were younger than those who had their ovaries removed (37.42 ± 5.79 years vs. 50.04 ± 8.23 years). Patients in the BSO group possessed more risk factors as expected, such as advanced stage, LNM, LVSI, DMI, parametric invasion, and UCI. The results of the univariate analysis in Table 2 also confirmed that ovarian metastasis should be associated with these risk factors, further supporting the validity of the selection criteria. Besides, different pathology types may have different risks of ovarian metastasis, while histopathology type was not included in the previous criteria. The association of ovarian metastasis with different histopathology types was analyzed (Table 3) and further classified into high- and low-risk histopathology types with a cutoff risk of 10%. Among the high-risk pathologic type patients, only 4/74 had ovaries preserved (Table 2), and the rate of ovarian metastasis for those who had their ovaries removed was as high as 17.14% (Table 4).
Figure 1
Table 2
| Total | BSO | OP | p-Value | |
|---|---|---|---|---|
| 668 | 266 | |||
| Age (n, %) | <0.001 | |||
| ≤50 | 640 (68.5) | 376 (56.3) | 264 (99.2) | |
| >50 | 294 (31.5) | 292 (43.7) | 2 (0.8) | |
| Histological classification (n, %) | <0.001 | |||
| Low risk | 858 (92.1) | 598 (89.5) | 260 (98.5) | |
| High risk | 74 (7.9) | 70 (10.5) | 4 (1.5) | |
| FIGO stage (n, %) | 0.002 | |||
| I | 679 (72.7) | 462 (69.2) | 217 (81.6) | |
| II | 85 (9.1) | 67 (10.0) | 18 (6.8) | |
| III | 163 (17.5) | 133 (19.9) | 30 (11.3) | |
| IV | 7 (0.7) | 6 (0.9) | 1 (0.4) | |
| Tumor size (n, %) | 0.075 | |||
| ≤4 cm | 717 (83.1) | 510 (81.6) | 207 (87.0) | |
| >4 cm | 146 (16.9) | 115 (18.4) | 31 (13.0) | |
| LNM (n, %) | 0.006 | |||
| Negative | 726 (80.8) | 508 (78.4) | 218 (87.2) | |
| Pelvic LN | 104 (11.6) | 82 (12.7) | 22 (8.8) | |
| Common iliac LN | 50 (5.6) | 40 (6.2) | 10 (4.0) | |
| Paraaortic LN | 18 (2.0) | 18 (2.8) | 0 (0.0) | |
| DMI (n, %) | <0.001 | |||
| Superficial 1/3 | 482 (51.6) | 304 (45.5) | 178 (66.9) | |
| Middle 1/3 | 29 (3.1) | 20 (3.0) | 9 (3.4) | |
| Deep 1/3 | 389 (41.6) | 315 (47.2) | 74 (27.8) | |
| Whole cervical layer | 34 (3.6) | 29 (4.3) | 5 (1.9) | |
| LUSI (n, %) | <0.001 | |||
| Negative | 727 (77.8) | 497 (74.4) | 230 (86.5) | |
| Positive | 207 (22.2) | 171 (25.6) | 36 (13.5) | |
| LVSI (n, %) | 0.001 | |||
| Negative | 583 (62.4) | 392 (58.7) | 191 (71.8) | |
| Sporadic | 334 (35.8) | 263 (39.4) | 71 (26.7) | |
| Extensive | 17 (1.8) | 13 (1.9) | 4 (1.5) | |
| Parametrial involvement (n, %) | 0.002 | |||
| Negative | 818 (91.9) | 579 (90.0) | 239 (96.8) | |
| Positive | 72 (8.1) | 64 (10.0) | 8 (3.2) | |
| Incision margin (n, %) | 0.006 | |||
| Negative | 895 (95.8) | 632 (94.6) | 263 (98.9) | |
| Positive | 39 (4.2) | 36 (5.4) | 3 (1.1) | |
| Vaginal invasion (n, %) | <0.001 | |||
| Negative | 832 (89.1) | 577 (86.4) | 255 (95.9) | |
| Positive | 102 (10.9) | 91 (13.6) | 11 (4.1) | |
| Fallopian metastasis (n, %) | 0.033 | |||
| Negative | 825 (96.7) | 641 (96.0) | 184 (99.5) | |
| Positive | 28 (3.3) | 27 (4.0) | 1 (0.5) | |
| Adjuvant chemotherapy (n, %) | <0.001 | |||
| No | 413 (46.4) | 264 (41.6) | 149 (58.4) | |
| Yes | 477 (53.6) | 371 (58.4) | 106 (41.6) | |
| Adjuvant radiotherapy (n, %) | <0.001 | |||
| No | 457 (51.5) | 294 (46.3) | 163 (64.7) | |
| Yes | 430 (48.5) | 341 (53.7) | 89 (35.3) | |
Clinicopathologic comparison in patients with ovarian preservation (OP) or bilateral salpingo-oophorectomy (BSO) during surgical treatment.
FIGO, International Federation of Gynecology and Obstetrics; LNM, lymph node metastasis; LUSI, lower uterine segment involvement; DMI, depth of myometrial invasion; LVSI, lymph vascular space invasion.
Table 3
| Pathological type (WHO 2014) | OM | |
|---|---|---|
| Yes | No | |
| Low risk | ||
| Villoglandular | 0 (0.0) | 3 (100.0) |
| Endometrioid | 0 (0.0) | 5 (100.0) |
| Intestinal type | 0 (0.0) | 21 (100.0) |
| Clear cell type | 0 (0.0) | 8 (100.0) |
| Mixed with neuroendocrine | 0 (0.0) | 15 (100.0) |
| Usual adenocarcinoma type | 3 (1.18) | 252 (98.82) |
| Mucinous, NOS | 1 (1.30) | 76 (98.70) |
| Adenosquamous | 4 (1.87) | 210 (98.13) |
| High risk | ||
| Gastric type | 5 (11.90) | 37 (88.10) |
| Serous | 2 (15.38) | 11 (84.62) |
| MDA | 4 (28.57) | 10 (71.43) |
| Signet ring cell type | 1 (100.0) | 0 (0.0) |
The OM distribution in different pathological subtypes of ADC.
Table 4
| 668 | BSO group | |||
|---|---|---|---|---|
| No OM | OM | p-Value | ||
| 647 | 21 | |||
| Age (n, %) | 0.138 | |||
| ≤50 | 376 (56.3) | 368 (56.9) | 8 (38.1) | |
| >50 | 292 (43.7) | 279 (43.1) | 13 (61.9) | |
| Histological classification (n, %) | <0.001 | |||
| Low risk | 598 (89.5) | 589 (91.0) | 9 (42.9) | |
| High risk | 70 (10.5) | 58 (9.0) | 12 (57.1) | |
| FIGO stage (n, %) | <0.001 | |||
| I | 462 (69.2) | 457 (70.6) | 5 (23.8) | |
| II | 67 (10.0) | 65 (10.0) | 2 (9.5) | |
| III | 133 (19.9) | 123 (19.0) | 10 (47.6) | |
| IV | 6 (0.9) | 2 (0.3) | 4 (19.0) | |
| Tumor size (n, %) | 0.025 | |||
| ≤4 cm | 510 (81.6) | 498 (82.3) | 12 (60.0) | |
| >4 cm | 115 (18.4) | 107 (17.7) | 8 (40.0) | |
| LNM (n, %) | <0.001 | |||
| Negative | 508 (78.4) | 502 (79.8) | 6 (31.6) | |
| Pelvic LN | 82 (12.7) | 77 (12.2) | 5 (26.3) | |
| Iliac LN | 40 (6.2) | 36 (5.7) | 4 (21.1) | |
| Paraaortic LN | 18 (2.8) | 14 (2.2) | 4 (21.1) | |
| DMI (n, %) | <0.001 | |||
| Superficial 1/3 | 304 (45.5) | 302 (46.7) | 2 (9.5) | |
| Middle 1/3 | 20 (3.0) | 20 (3.1) | 0 (0.0) | |
| Deep 1/3 | 315 (47.2) | 301 (46.5) | 14 (66.7) | |
| Whole cervical layer | 29 (4.3) | 24 (3.7) | 5 (23.8) | |
| LUSI (n, %) | <0.001 | |||
| Negative | 497 (74.4) | 493 (76.2) | 4 (19.0) | |
| Positive | 171 (25.6) | 154 (23.8) | 17 (81.0) | |
| LVSI (n, %) | <0.001 | |||
| Negative | 392 (58.7) | 385 (59.5) | 7 (33.3) | |
| Sporadic | 263 (39.4) | 255 (39.4) | 8 (38.1) | |
| Extensive | 13 (1.9) | 7 (1.1) | 6 (28.6) | |
| Parametrial involvement (n, %) | <0.001 | |||
| Negative | 579 (90.0) | 569 (91.0) | 10 (55.6) | |
| Positive | 64 (10.0) | 56 (9.0) | 8 (44.4) | |
| Incision margin (n, %) | <0.001 | |||
| Negative | 632 (94.6) | 621 (96.0) | 11 (52.4) | |
| Positive | 36 (5.4) | 26 (4.0) | 10 (47.6) | |
| Vaginal invasion (n, %) | <0.001 | |||
| Negative | 577 (86.4) | 569 (87.9) | 8 (38.1) | |
| Positive | 91 (13.6) | 78 (12.1) | 13 (61.9) | |
| Fallopian metastasis (n, %) | <0.001 | |||
| Negative | 641 (96.0) | 631 (97.5) | 10 (47.6) | |
| Positive | 27 (4.0) | 16 (2.5) | 11 (52.4) | |
| Adjuvant chemotherapy (n, %) | 0.001 | |||
| No | 264 (41.6) | 264 (42.7) | 0 (0.0) | |
| Yes | 371 (58.4) | 354 (57.3) | 17 (100.0) | |
| Adjuvant radiotherapy (n, %) | 0.031 | |||
| No | 294 (46.3) | 291 (47.1) | 3 (17.6) | |
| Yes | 341 (53.7) | 327 (52.9) | 14 (82.4) | |
Clinicopathologic characteristics associated with OM in the BSO group.
OM and survival analysis in patients with ADC
The overall 5-year OS among our cohort with 934 patients was 88.8%. In patients with adnexectomy, the total rate of OM was 3.14%. For patients with OM, all of them received adjuvant chemotherapy, and 82.4% had radiation therapy (RT) (Table 4). Compared to patients without OM, the 5-year DFS (Figure 2C) and OS (Figure 2D) of patients with OM decrease significantly, which indicated that surgery plus adjuvant chemoradiotherapy might not improve the outcome of patients with OM. On the other hand, there was no significant difference in 5-year DFS and OS between the BSO and OP groups (p = 0.09 and 0.069, Figures 2A, B).
Figure 2
Construction of the predictive model for OM in patients with ADC
Patients who underwent BSO were divided it into training and test sets for machine learning of ovarian metastasis at a ratio of 6:4. No difference was found between these two sets regarding the baseline clinicopathological characteristics (Table 5). Clinicopathological variables were analyzed for association with OM (Figure 3), and cutoff value for age and tumor size as continuous variables were defined by ROC curves at 50-year old and 4 cm, respectively (Figure 3A). Variables significantly associated with OM screened by univariate regression analysis in the training set included fallopian tube metastasis, LVSI (categorized as extensive and non-extensive as no significant difference was found between no LVSI and sporadic LVSI), DMI, LNM metastasis, incision margin status, UCI, pathological type, and tumor size. The two critical variables associated with OM fallopian tube metastasis and extensive LVSI status were proven by the OR value generated by univariate analysis (Figure 3B) and random forest model analysis.
Figure 3
Table 5
| n = 668 | Cohort | |||
|---|---|---|---|---|
| Training | Testing | p-Value | ||
| 401 | 267 | |||
| Age (n, %) | 0.847 | |||
| ≤50 | 376 (56.3) | 224 (55.9) | 152 (56.9) | |
| >50 | 292 (43.7) | 177 (44.1) | 115 (43.1) | |
| Histological classification (n, %) | 0.902 | |||
| Low risk | 598 (89.5) | 358 (89.3) | 240 (89.9) | |
| High risk | 70 (10.5) | 43 (10.7) | 27 (10.1) | |
| Tumor size (n, %) | 0.646 | |||
| ≤4 cm | 449 (71.8) | 271 (72.7) | 178 (70.6) | |
| >4 cm | 176 (28.2) | 102 (27.3) | 74 (29.4) | |
| LNM (n, %) | 0.574 | |||
| No | 508 (78.4) | 309 (79.2) | 199 (77.1) | |
| Pelvic | 82 (12.7) | 51 (13.1) | 31 (12.0) | |
| Common iliac | 40 (6.2) | 21 (5.4) | 19 (7.4) | |
| Paraaortic | 18 (2.8) | 9 (2.3) | 9 (3.5) | |
| DMI (n, %) | 0.604 | |||
| Superficial and middle | 324 (48.5) | 190 (47.4) | 134 (50.2) | |
| Deep | 315 (47.2) | 195 (48.6) | 120 (44.9) | |
| Whole | 29 (4.3) | 16 (4.0) | 13 (4.9) | |
| LUSI (n, %) | 0.835 | |||
| Negative | 497 (74.4) | 300 (74.8) | 197 (73.8) | |
| Positive | 171 (25.6) | 101 (25.2) | 70 (26.2) | |
| LVSI (n, %) | 0.456 | |||
| Non-extensive | 655 (98.1) | 395 (98.5) | 260 (97.4) | |
| Extensive | 13 (1.9) | 6 (1.5) | 7 (2.6) | |
| Parametrial involvement (n, %) | 0.558 | |||
| Negative | 579 (90.0) | 344 (89.4) | 235 (91.1) | |
| Positive | 64 (10.0) | 41 (10.6) | 23 (8.9) | |
| Incision margin (n, %) | 0.698 | |||
| Negative | 632 (94.6) | 381 (95.0) | 251 (94.0) | |
| Positive | 36 (5.4) | 20 (5.0) | 16 (6.0) | |
| Vaginal invasion (n, %) | 0.508 | |||
| Negative | 577 (86.4) | 343 (85.5) | 234 (87.6) | |
| Positive | 91 (13.6) | 58 (14.5) | 33 (12.4) | |
| Fallopian metastasis (%) | 0.085 | |||
| Negative | 641 (96.0) | 380 (94.8) | 261 (97.8) | |
| Positive | 27 (4.0) | 21 (5.2) | 6 (2.2) | |
| Ovarian metastasis (%) | 0.34 | |||
| Negative | 647 (96.9) | 391 (97.5) | 256 (95.9) | |
| Positive | 21 (3.1) | 10 (2.5) | 11 (4.1) | |
Clinicopathologic characteristics of the training and testing population in 668 patients who accepted BSO.
Our data showed that no significant difference was found between no LVSI and sporadic LVSI affecting OM; also, no difference was found between the superficial infiltration and middle one-third invasion of the cervix in affecting OM, so we categorized these two parameters into extensive LVSI vs. non-extensive LVSI, and deep + whole muscular invasion vs. superficial + middle one-third muscular infiltration of the cervix for further analysis.
Machine learning was performed with the following four methods: random forest, stepwise regression, LASSO, and optimal subset. The random forest model obtained the highest AUC in both the training and testing sets (0.971 for training set and 0.962 for testing set), considered to be the optimal predictive model. Also, the OM high- and low-risk subgroups were divided according to the optimal cutoff point of the ROC curve with its specificity at 99.5% and sensitivity at 90% (Figures 4A, B). Both DFS and OS differed significantly with this categorization (Figures 4C, D), which was further validated in the testing set (Figures 4E, F) with ovarian metastasis rate as high as 25.00% in the high-risk group and 0.44% in the low-risk group. On the basis of this random forest model, we further utilized nomogram to predict OM risk and validated the model with calibration curves (Figure 5). Details of the other three methods (stepwise regression, LASSO, and optimal subset) are displayed in Supplementary Figure S1.
Figure 4
Figure 5
Security verification of the predictive model
Considering the extremely poor prognosis of the occurring ovarian metastases, the low-risk group, predicted by the constructed model, who were recommended to preserve their ovaries, must be guaranteed to be secure. After redefining the whole population into low- and high-risk groups with the predictive model, the low-risk group included 822 patients who were recommended to preserve their ovaries, increasing the OP population screened by the current criteria by 3.21-fold (Table 6). The 5-year DFS and OS of the low-risk group had no difference compared with the patients with ovary preservation screened by the current criteria (Figure 6), with OM in the low-risk group at only 0.24% and the high-risk group at 16.96% in the total population, which proved the effectiveness and security of this predictive model.
Table 6
| OM | OP | |||||
|---|---|---|---|---|---|---|
| No | Yes | p-Value | No | Yes | p-Value | |
| High risk | 93 (10.2) | 19 (90.5) | <0.001 | 102 (15.3) | 10 (3.8) | <0.001 |
| Low risk | 820 (89.8) | 2 (9.5) | 566 (84.7) | 256 (96.2) | ||
| Total | 913 | 21 | 668 | 266 | ||
The evaluation of the total population with the predictive model.
Figure 6
Discussion
Patients with early-stage ADC usually have optimistic oncological outcomes, with the incidence of OM at 0.5%–11% (
OP could remarkably benefit young patients not only to retain fertility but also improve their quality of life through saving their ovarian endocrine function. Once ovaries are removed, patients sharply seek iatrogenic menopause, which could lead to severe menopausal symptoms, including early hot flashes, vaginal atrophy, an increased risk of cardiovascular disease, osteoporosis, hip fracture, Alzheimer’s disease, and mental health disorders (
Most of the parameters in the previous selection criteria mainly depend upon biopsy pathological results and imaging evidence. LNM status mainly relied on the pre-operative imaging test; however, it was reported that the preoperative predictive accuracy of LNM by PET-CT was only approximately 70% (
Table 7
| Preoperative patients’ characteristics |
|---|
| –Premenopausal women |
| –Patients who desire to retain ovarian function |
| –No familial predisposition to ovarian cancer |
| –OM rate calculated based on preoperative examination ≤ cutoff value* |
| Intraoperative findings |
| –No extra-uterine spread except LNM |
| –Normal appearance of ovaries |
| –OM rate calculated based on pre- and intra-operative findings ≤ cutoff value |
| Postoperative pathology |
| –OM rate ≤ cutoff value, close surveillance according to the routine of cervical cancer follow-up |
Updated selection criteria for OP recommended.
*The threshold for ovarian metastasis rate to have BSO performed could be individualized and decided with the gynecology oncologist and the patient together.
A multidisciplinary team should be summoned to implement this predictive model; especially, the imaging and pathology departments, which might be limited in resource-poor areas, should be consulted when making decisions. Additionally, it should be emphasized that patients could be assessed as low-risk group pre- and intra-operation, while upgraded as high-risk for OM post-operation, which required another operation to remove the ovaries. These risks should be carefully communicated with patients before OP surgery using this predictive model.
Our study constructed a prediction model for ovarian metastasis in patients with cervical adenocarcinoma/adenosquamous carcinoma for the first time to the best of our knowledge. Current clinical researches are more concerned on how to identify the population with a high risk of OM and if the ovaries preserved have a higher potential of subsequent recurrence than other organs. The OM risk reported previously has been inconsistent. Ritsu Yamamoto et al. found that the metastatic rate of adenocarcinoma was as high as 10.2% (10/98 = 10.2%) (
Cheng et al. reviewed 10 studies on ADC, staged at I–II with surgical treatment, by meta-analysis, with sample sizes of 5,075, and 1,199 who had their ovaries preserved, and 3,867 had their ovaries removed. The overall metastatic rate was 3.61% in patients with ovaries removed, and there was no significant difference in prognosis between patients with preserved ovaries and those with ovaries removed [OS (OR 1.00, 95% CI 0.64–1.56, I2 = 25.7%), PFS (OR 0.98, 95% CI 0.57–1.66, I2 = 0%)] (
For the second concern raised here, Giovanni Scambia et al. concluded that the rate of ovarian recurrence in cervical cancer patients stayed at a relatively low level (0%–1.3%) through a systematic literature review (
Although our predictive model is optimized through machine learning, it still has some limitation. First, it is based on a single-center, retrospective study, and a further multi-center prospective randomized control trial is needed to help validate the model. Second, some parameters in our model require high-quality intraoperative frozen pathological reports in clinical setting. Risk factors, such as fallopian tube metastasis, UCI, and LNM (lymph nodes with suspicious preoperative imaging or intraoperative microscopic metastasis), need intraoperative frozen pathology to detect and help make a clinical decision if OP could safely be performed. Additionally, our model is based on the WHO 2014 pathological classification, while a new classification system (International Endocervical Adenocarcinoma Criteria and Classification, IECC) was established in 2018 (
Conclusion
Ovarian preservation is remarkably important for young patients with ADC to maintain their life quality, and concurrently, patients with OM have extremely poor prognosis, so the individual decision of whether to preserve the ovaries before surgical treatment should be accurately made. Our study constructed a precise predictive model through machine learning with a decent sample size of 934 patients with ADC, which could potentially reduce 60% of young patients with ADC to safely preserve their ovaries compared with the current selective criteria.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Ethics statement
This study was approved by the Ethic Committee of The Obstetrics and Gynecology Hospital of Fudan University (No. 2021-34) on 18 January 2021. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
YZ: Writing – original draft, Writing – review & editing, Conceptualization. YS: Writing – original draft, Writing – review & editing, Supervision. XX: Data curation, Formal Analysis, Methodology, Software, Writing – original draft. JD: Conceptualization, Writing – review & editing. KH: Conceptualization, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research project was funded by Shanghai Municipal Health Commission (no. 20214Y0021).
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.
Publisher’s note
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2024.1464565/full#supplementary-material
Supplementary Figure 1(A) LASSO coefficient profiles of the 13 key clinicopathological variables for the prediction of ovarian metastasis in the training cohort. (B) Tuning parameter selection by 10-fold cross-validation in the LASSO model of training cohort. The partial likelihood deviance was plotted against log (lambda/λ), and λ was the tuning parameter. The partial likelihood deviance values are shown, and error bars represented s.e. The dotted vertical lines showed the optimal values through minimum criteria and 1−s.e. criteria. (C) Tenfold cross-validation showed that three variables could minimize the error of optimum subset regression analysis in the training cohort. The “bestglm” R package was used.
References
1
SiegelRLMillerKDWagleNSJemalA. Cancer statistics, 2023. CA: Cancer J Clin. (2023) 73:17–48. doi: 10.3322/caac.21763
2
PeiJShuTWuCLiMXuMJiangMet al. Impact of human papillomavirus vaccine on cervical cancer epidemic: evidence from the surveillance, epidemiology, and end results program. Front Public Health. (2022) 10:998174. doi: 10.3389/fpubh.2022.998174
3
SmithAJBBeavisALRositchAFLevinsonK. Disparities in diagnosis and treatment of cervical adenocarcinoma compared with squamous cell carcinoma: an analysis of the national cancer database, 2004-2017. J lower genital tract Dis. (2023) 27:29–34. doi: 10.1097/lgt.0000000000000702
4
LoureiroJOlivaE. The spectrum of cervical glandular neoplasia and issues in differential diagnosis. Arch Pathol Lab Med. (2014) 138:453–83. doi: 10.5858/arpa.2012-0493-RA
5
BrayFCarstensenBMøllerHZappaMZakeljMPLawrenceGet al. Incidence trends of adenocarcinoma of the cervix in 13 European countries. Cancer epidemiology Biomarkers Prev. (2005) 14:2191–9. doi: 10.1158/1055-9965.epi-05-0231
6
LandoniFZanagnoloVLovato-DiazLManeoARossiRGadducciAet al. Ovarian metastases in early-stage cervical cancer (Ia2-iia): A multicenter retrospective study of 1965 patients (a cooperative task force study). Int J gynecological Cancer. (2007) 17:623–8. doi: 10.1111/j.1525-1438.2006.00854.x
7
YamamotoROkamotoKYukiharuTKaneuchiMNegishiHSakuragiNet al. A study of risk factors for ovarian metastases in stage ib-iiib cervical carcinoma and analysis of ovarian function after a transposition. Gynecologic Oncol. (2001) 82:312–6. doi: 10.1006/gyno.2001.6277
8
ChengHHuoLZongLKongYYangJXiangY. Oncological outcomes and safety of ovarian preservation for early stage adenocarcinoma of cervix: A systematic review and meta-analysis. Front Oncol. (2019) 9:777. doi: 10.3389/fonc.2019.00777
9
MatsuoKShimadaMYamaguchiSKanaoHNakanishiTSaitoTet al. Identifying a candidate population for ovarian conservation in young women with clinical stage ib-iib cervical cancer. Int J Cancer. (2018) 142:1022–32. doi: 10.1002/ijc.31084
10
ZhouJChenYZhangPLouH. Ovarian preservation in adenocarcinoma of the uterine cervix. J Ovarian Res. (2017) 10:48. doi: 10.1186/s13048-017-0339-y
11
ZhouLSunCTLinLXieYHuangYLiQet al. Independent risk factors for ovarian metastases in stage ia-iib cervical carcinoma. Acta obstetricia gynecologica Scandinavica. (2019) 98:18–23. doi: 10.1111/aogs.13442
12
TouhamiOPlanteM. Should ovaries be removed or not in (Early-stage) adenocarcinoma of the uterine cervix: A review. Gynecologic Oncol. (2015) 136:384–8. doi: 10.1016/j.ygyno.2014.12.011
13
BizzarriNPavoneMLoverroMQuerleuDFagottiAScambiaG. Ovarian preservation in gynecologic oncology: current indications and techniques. Curr Opin Oncol. (2023) 35:401–11. doi: 10.1097/cco.0000000000000969
14
ShengJYiYNieJ. Influence of ovarian-sparing surgery and ovariectomy on prognosis in early cervical adenocarcinoma: A systematic review and meta-analysis. Gland Surg. (2022) 11:1103–10. doi: 10.21037/gs-22-310
15
TheplibAHanprasertpongJLeetanapornK. Safety and prognostic impacts of ovarian preservation during radical hysterectomy for early-stage adenocarcinoma and adenosquamous cervical cancer. BioMed Res Int. (2020) 2020:5791381. doi: 10.1155/2020/5791381
16
HuJJiaoXYangZCuiHGuoHWuYet al. Should ovaries be removed or not in early-stage cervical adenocarcinoma: A multicenter retrospective study of 105 patients. J Obstet Gynaecol. (2017) 37:1065–9. doi: 10.1080/01443615.2017.1323198
17
BhatlaNBerekJSCuello FredesMDennyLAGrenmanSKarunaratneKet al. Corrigendum to "Revised figo staging for carcinoma of the cervix uteri" [Int J Gynecol Obstet 145(2019) 129-135. Int J gynaecology obstetrics. (2019) 147:279–80. doi: 10.1002/ijgo.12969
18
NezhatCNezhatFBurrellMOBenignoBWelanderCE. Laparoscopic radical hysterectomy with paraaortic and pelvic node dissection. Am J obstetrics gynecology. (1994) 170:699. doi: 10.1016/s0002-9378(94)70251-9
19
ShimadaMKigawaJNishimuraRYamaguchiSKuzuyaKNakanishiTet al. Ovarian metastasis in carcinoma of the uterine cervix. Gynecologic Oncol. (2006) 101:234–7. doi: 10.1016/j.ygyno.2005.10.004
20
RoccaWAGrossardtBRde AndradeMMalkasianGDMeltonLJ3rd. Survival patterns after oophorectomy in premenopausal women: A population-based cohort study. Lancet Oncol. (2006) 7:821–8. doi: 10.1016/s1470-2045(06)70869-5
21
MunkarahARJhingranAIyerRBWallaceSEifelPJGershensonDet al. Utility of lymphangiography in the prediction of lymph node metastases in patients with cervical cancer. Int J gynecological Cancer. (2002) 12:755–9. doi: 10.1046/j.1525-1438.2002.t01-1-01150.x
22
LevinsonKBeavisALPurdyCRositchAFViswanathanAWolfsonAHet al. Beyond sedlis-a novel histology-specific nomogram for predicting cervical cancer recurrence risk: an Nrg/Gog ancillary analysis. Gynecologic Oncol. (2021) 162:532–8. doi: 10.1016/j.ygyno.2021.06.017
23
StolnicuSBarsanIHoangLPatelPTerinteCPesciAet al. International endocervical adenocarcinoma criteria and classification (Iecc): A new pathogenetic classification for invasive adenocarcinomas of the endocervix. Am J Surg Pathol. (2018) 42:214–26. doi: 10.1097/pas.0000000000000986
Summary
Keywords
cervical adenocarcinoma, ovarian metastasis, ovarian preservation, predictive model and machine learning, random forest
Citation
Zhang Y, Shi Y, Xiang X, Ding J and Hua K (2024) How do we safely preserve ovaries in patients with cervical adenocarcinoma: risk factors and predictive models. Front. Oncol. 14:1464565. doi: 10.3389/fonc.2024.1464565
Received
14 July 2024
Accepted
08 October 2024
Published
30 October 2024
Volume
14 - 2024
Edited by
Kai Ding, University of Oklahoma Health Sciences Center, United States
Reviewed by
Mikel Gorostidi, University of the Basque Country, Spain
Monica Ewomazino Akokuwebe, University of the Witwatersrand, South Africa
Mirella Fortunato, Azienda Sanitaria Ospedaliera S.Croce e Carle Cuneo, Italy
Xutong Zheng, China Medical University, China
Updates

Check for updates
Copyright
© 2024 Zhang, Shi, Xiang, Ding and Hua.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Keqin Hua, huakeqinjiaoshou@163.com; Jingxin Ding, djxdd@sina.com
†These authors have contributed equally to this work
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.