This retrospective multicenter cohort study received ethical approval from the Institutional Review Board of the Obstetrics and Gynecology Hospital of Fudan University (Shanghai, China; approval no. 2019‐87) and registered with the Chinese Clinical Trial Registry (ChiCTR1900028702). Oral informed consent was obtained from all participants via telephone follow-up in accordance with institutional ethical standards. All personally identifiable information was permanently removed prior to analysis. Data were coded using unique identifiers and stored securely on encrypted, password-protected servers with role-based access controls to ensure confidentiality. No compensation was provided to participants.

A total of 1076 patients with CC who underwent surgical resection between January 2014 and December 2022 were identified from 2 tertiary hospitals in China: the Obstetrics and Gynecology Hospital affiliated with Fudan University and Shanghai First Maternity and Infant Hospital. The inclusion criteria were as follows: (1) pathologically confirmed FIGO stage IA1 with positive lymph-vascular space invasion (LVSI) to stage IIB CC; (2) radical hysterectomy performed according to National Comprehensive Cancer Network (NCCN) guidelines appropriate to the disease stage at the time [19-21]; (3) availability of high-quality preoperative pelvic contrast-enhanced MR images, including the coronal T2-weighted imaging (T2WI) fat-suppressed sequence; and (4) at least 3 years of follow-up data. The exclusion criteria were as follows: (1) receipt of chemotherapy or radiotherapy before surgery (n=23); (2) incomplete medical records (n=61); (3) absence of preoperative MR imaging (MRI) or MRI performed at another institution (n=536); (4) unsatisfactory MR image quality (n=98); and (5) loss to follow-up (n=108). The final study population comprised 250 patients (Figure 1).

For eligible patients, demographic data, laboratory test results, treatment details, and tumor characteristics were collected from medical records. All records were reviewed concurrently by 3 experts and independently verified for accuracy by 2 additional reviewers. Following diagnosis, demographic variables—including age and comorbidity (hypertension or diabetes)—as well as laboratory data, including squamous cell carcinoma antigen levels and human papillomavirus (HPV) infection status, were recorded. The history of the loop electrosurgical excision procedure was also noted.

All surgical procedures during the study period were performed by faculty members with completed fellowship training in gynecologic oncology. Surgical data included surgical approach, operative time, estimated blood loss, and use of blood transfusion. Tumor characteristics included FIGO stage, tumor size, histologic type, depth of stromal invasion (DSI), LVSI, surgical margin status, parametrial involvement, lymph node metastasis, keratinization, degree of differentiation, and expression of P53, P16, and Ki-67. According to the NCCN guidelines [3], all patients received adjuvant therapy if they met one of the following criteria: (1) presence of any high-risk factor, including positive surgical margins, parametrial involvement, or lymph node metastasis; or (2) fulfillment of any Sedlis criteria [22] for intermediate-risk factors, including tumor size, LVSI, and DSI.

Patients were followed according to the 2022 NCCN guidelines [3] after discharge and completion of initial treatment. HPV testing, liquid-based cytology, tumor marker evaluation, and ultrasonography were conducted every 3 months for the first 2 years, every 6 months for the next 2 years, and annually thereafter. Chest computed tomography (CT), contrast-enhanced upper abdominal CT, and pelvic MRI were performed annually. Telephone follow-ups were also conducted, and patients with symptoms or abnormal findings suggestive of recurrence were advised to undergo the aforementioned tests. In cases of suspected organ or lymph node metastasis, needle aspiration biopsy was performed when clinically indicated.

Coronal preoperative T2WI fat-suppressed pelvic contrast-enhanced MR images were collected for all 250 patients. This included both the original MR sequence source images and lesion segmentation images annotated by radiologists. The MR source images were first normalized before further processing. To standardize the 2 types of MR source image resolutions (256×256 pixels and 320×320 pixels), images of 320×320 pixels were center-cropped to 256×256 pixels. Meanwhile, the pixel intensity values, which originally ranged from 0 to over 1000, were rescaled to the range [0, 1]. After obtaining standardized MR source images, 2 experienced radiologists performed lesion annotations. The region of interest for CC was manually delineated on T2WI images using ITK-SNAP (version 3.8.0). Each radiologist delineated half of the cases and independently reviewed the other half. In cases of disagreement, the final decision was reached through discussion or by consulting a third radiologist. During the segmentation process, the radiologists were blinded to the patients’ clinical information.

After obtaining standardized MR source images, 2 experienced radiologists performed lesion annotations. The region of interest for CC was manually delineated on T2WI images using ITK-SNAP (version 3.8.0). Each radiologist delineated half of the cases and independently reviewed the other half. In cases of disagreement, the final decision was reached through discussion or by consulting a third radiologist. During the segmentation process, the radiologists were blinded to the patients’ clinical information.

Subsequently, cervical lesion identification was modeled using the radiologist-annotated MR images. An optimized algorithm based on the 3D U-Net architecture [23], embedded with a squeeze-and-excitation layer [24], was adopted to segment the lesion pixels. Following multiple preprocessing steps (see Supplementary Materials in Multimedia Appendix 1), the input data (C×W×H×D) were compressed and normalized using the Sigmoid activation function [25], then reshaped to their original dimensions. The grayscale output was segmented by computing the probability of each pixel being classified as a positive sample through the Sigmoid function. To address the issue of sample imbalance—where MR image sequences contained significantly more frames with lesion markers than without—an improved Focal Loss function [26] was implemented. A weighting factor (α) was applied to balance positive and negative samples, while a power function was used to reduce the loss contribution of easy samples. The optimal parameters used in this study were α=.25 and γ=2.

Parametrial invasion and lymphatic metastasis were 2 key pathological indicators derived from the surgical specimens of cervical carcinoma. Parametrial invasion refers to the infiltration of neoplastic cells beyond the cervix into the surrounding connective tissue (parametrium), while lymphatic metastasis denotes the spread of malignant cells to regional lymph nodes.

For survival outcomes, recurrence-free survival (RFS) was defined as the interval from the initial CC diagnosis to either the first documented recurrence or the last follow-up. Overall survival (OS) was defined as the interval from diagnosis to CC–related death or the last follow-up. Regarding recurrence classification, local recurrence was defined as the pathologically confirmed first reappearance of cancer in the cervix or vagina following complete treatment, confined to the pelvic region. Distant recurrence was defined as the first pathologically confirmed relapse beyond the pelvis, including peritoneal dissemination or metastasis to distant organs.

After selecting MR source images that corresponded to lesion-annotated segmentation files, the filtered original images (256×256 pixels) were flattened into 65,536-dimensional vectors. The clinical dataset consisted of 22 variables (Table S1 in Multimedia Appendix 1), including both continuous and categorical features. Continuous variables were normalized using min-max scaling, rescaling each feature into the [0, 1] range using the following formula:

where Xmin and Xmax represent the minimum and maximum observed values of that feature. For example, for a variable ranging from 2 to 10, a value of 5 would be transformed as: (5–2)/(10–2)=0.375. Categorical variables were encoded using one-hot encoding, which converts nominal variables into orthogonal binary vectors. For a categorical feature with 3 classes (1, 2, 3), the transformation was as follows: Class 1 → [1, 0, 0]; Class 2 → [0, 1, 0]; and Class 3 → [0, 0, 1].

After preprocessing, categorical and continuous variables were concatenated into 5-dimensional (preoperative) and 57-dimensional (postoperative) clinical feature vectors. These clinical vectors were then fused with the 65,536-dimensional image vectors, resulting in integrated feature vectors of 65,593 dimensions.

Two integrated models were developed: (1) The preoperative recognition model integrated preoperative clinical parameters (age, comorbidity, HPV status, squamous cell antigen carcinoma level, and loop electrosurgical excision procedure history) and MR images to generate binary classifications for parametrial invasion and lymph node metastasis; and (2) the postoperative prognostic model combined 21 postoperative clinical parameters (Table S1 in Multimedia Appendix 1) and MR images to predict recurrence, mortality, and individualized RFS or OS times.

A total of 7 ML algorithms were used for model construction. These included classical ML algorithms such as K-nearest neighbor (KNN), support vector machine (SVM), decision tree (DT), and random forest (RF), as well as their imbalanced-data variants—weighted KNN, balanced RF, and weighted DT. Weighted KNN adjusted sample weights inversely proportional to class frequency. Specifically, minority class samples (in this study, positive cases) were assigned higher weights in distance calculations, forcing the algorithm to prioritize minority class neighborhoods during prediction. Balanced RF used under-sampling by randomly removing majority class samples (negative cases) in each bootstrap iteration to create balanced subsets for training individual trees. Weighted DT modified the splitting criterion by incorporating a class-weighted penalty: a weight of misclassifying a positive sample was set to amplify the cost of false negatives during node splitting. Besides, to mitigate the influence of class imbalance during model development, negative samples were down-sampled to achieve a 1:1 ratio with positive cases during training.

The study process consisted of 4 main stages: patient data input, data preprocessing, model development, and model evaluation (Figure 1). Stratified random sampling was used to divide the complete dataset from the 2 hospitals into training and test sets at a ratio of 8:2. Model training was conducted using the training set, while final validation was performed on the test set. To minimize overfitting and reduce bias, 5-fold cross-validation was used on the training set for hyperparameter tuning.

Model performance for classification tasks was evaluated using sensitivity, specificity, accuracy, precision, F1-score, weighted accuracy, and area under the receiver operating characteristic curve (AUC). The mean absolute error (MAE) and concordance index (C-index) were used to assess the accuracy of individual-specific RFS and OS time predictions.

To facilitate clinical application and enhance accessibility for physicians, a web-based predictive diagnostic support tool was developed using Python, enabling DICOM upload and automated prediction.

Continuous variables were reported as means with standard deviations (SDs) for normally distributed data and as medians with interquartile ranges (IQRs) for non-normally distributed data, while categorical variables were summarized as counts and percentages. No significant interactions were observed among variables based on correlation matrix analysis. Concordance of continuous variables was assessed using intraclass correlation coefficients to evaluate consistency between pathology and imaging reports. Categorical variables were analyzed using the χ² test and the Kappa coefficient.

All statistical analyses were performed using R statistical software (version 4.1.1; R Foundation for Statistical Computing) and Python programming software (version 3.10.4; Python Software Foundation). All tests were 2-sided, and a P value of less than .05 was considered statistically significant unless otherwise stated.

A total of 250 patients with FIGO stage IA1 (LVSI+) to IIB CC who underwent radical hysterectomy between 2014 and 2022 were included in the study population (Table 1). The median age was 48.8 years. Most patients were at stage I (n=182, 72.8%) and had a squamous histologic type (n=200, 80%). More than half of the patients (n=191, 76.4%) received adjuvant therapy. In total, 24 women (9.6%) experienced recurrence, and 11 (4.4%) died during the follow-up period. Among the recurrence cases, 9 (37.5%) were local recurrences and 15 (62.5%) were distant. Among the distant recurrences, 4 (16.7%) occurred in the thoracic region, 5 (20.8%) in the abdominal region, and 6 (25%) in bone. The median RFS was 33.8 (IQR 24.9‐42.4) months, and the median OS was 34.6 (IQR 26.0‐42.8) months. The 1-year RFS and OS rates were 94.4% and 99.2%, respectively, while the 3-year RFS and OS rates were 90.1% and 95.4%, respectively (Figure S1 in Multimedia Appendix 1).

Integrated models combining MR images with 5 clinical parameters demonstrated variable performance across 7 ML algorithms (Table 2).

For parametrial invasion diagnosis, classical ML models generally exhibited poor sensitivity (0.00‐0.57). Specifically, RF and KNN failed to detect any positive cases (sensitivity=0.00), while SVM and DT showed low sensitivity (0.57 and 0.50) and specificity (0.56 and 0.59). On the contrary, imbalanced-data variants exhibited improved performance: (1) Balanced RF achieved balanced metrics (sensitivity=0.81, specificity=0.85, F1-score=0.64); (2) weighted KNN attained high sensitivity (0.98), while weighted DT prioritized specificity (0.93).

For lymph node metastasis detection, classical ML models showed limited sensitivity (0.31‐0.66), with RF achieving high specificity (0.87) and SVM delivering the best traditional performance (sensitivity=0.66, specificity=0.52, F1-score=0.54). Imbalanced-data variants exhibited enhancements: (1) Weighted KNN demonstrated optimal overall performance (sensitivity=0.67, specificity +=0.61, F1-score=0.58); (2) balanced RF achieved high specificity (0.87) and precision (0.68); and (3) weighted DT showed marginal improvement over classical DT (sensitivity=0.38 vs 0.53, F1-score=0.42 vs 0.49).

For prognostic prediction (Table 3), all classical ML models (KNN, SVM, DT, and RF) failed to identify recurrence or mortality cases (sensitivity=0.00). Imbalanced data variants exhibited divergent outcomes: For recurrence prediction, (1) weighted KNN achieved clinically actionable performance (sensitivit=0.80, specificity=0.96, F1-score=0.73); (2) balanced RF detected recurrences effectively (sensitivity=0.80) but exhibited poor specificity (0.37); and (3) weighted DT failed to identify recurrence cases (sensitivity=0.00). For mortality prediction, (1) weighted KNN showed perfect specificity and precision (specificity=0.99, precision=0.99) but demonstrated low sensitivity (0.33); (2) balanced RF showed limited detection capability (sensitivity=0.33, F1-score=0.11); and (3) weighted DT failed to identify mortality cases (sensitivity=0.00). Among all models, weighted KNN yielded the highest AUC for recurrence (0.861) and mortality (0.765) (Figure 2). Given the importance of sensitivity in guiding postoperative treatment and follow-up, weighted accuracy was introduced as an evaluation metric, and multiple weighting strategies were tested (Table S2 in Multimedia Appendix 1). Results revealed that increasing the weight of sensitivity led to improved detection accuracy for potential recurrence and mortality cases.

For individualized survival time prediction (Table 4), weighted KNN achieved the best performance for recurrence with the lowest MAE (8.53 mo) and highest C-index (0.83), while regression tree (RT) achieved the best performance for mortality (MAE=4.36 mo; C-index=0.99).

To enhance the clinical applicability of the prediction model and improve physicians’ access, usability, and integration into practice, the optimal artificial intelligence (AI) prediction models were embedded into a web-based software platform designed to assist in preoperative evaluation and prognosis prediction on MR images (Figure 3) [27]. By inputting key preoperative clinical parameters and uploading DICOM (DCM) files of preoperative MR images, the system can automatically identify possible lesions and generate predictions regarding the patient’s risk of parametrial invasion and lymphatic metastasis. Similarly, by entering relevant postoperative clinical, surgical, and pathological data along with uploading postoperative MR images, the system can estimate the patient’s risk of recurrence and mortality, as well as predict individualized durations of RFS and OS.

This study developed integrated ML models combining clinicopathological characteristics with MRI data for preoperative assessment of parametrial invasion and lymph node metastasis, alongside postoperative prognosis prediction in early-stage CC. Our key findings demonstrate the following: (1) Integrated models improved diagnosing sensitivity of parametrial invasion and lymph node metastasis, enhancing preoperative staging accuracy to better guide surgical plan; (2) postoperative prognostic models achieved robust performance in individualized recurrence (AUC 0.861) and survival prediction (AUC 0.765), providing a novel tool for precise adjuvant treatment and individualized follow-up strategy; and (3) a clinically deployable AI platform operationalizes these models for clinical workflow integration, enabling automated risk stratification.

The core strength of our study lies in developing a comprehensive AI-driven clinical support system that uniquely bridges preoperative staging and postoperative prognosis. By integrating MRI with clinicopathological data through imbalanced-data algorithms (eg, balanced RF and weighted KNN), we significantly enhanced sensitivity for detecting parametrial invasion and lymph node metastasis preoperatively alongside postoperative recurrence and mortality. Crucially, this framework was operationalized through a clinically deployable platform enabling automated DICOM analysis, risk stratification, and individualized survival time prediction, thus forming a closed-loop decision pathway from diagnosis to follow-up.

Despite these advances, several limitations merit careful consideration. Foremost among these is the limited population and tertiary-center recruitment bias. The patient population in this study was limited and may not adequately represent the broader clinical population, particularly individuals in primary care settings or from different geographic regions. Patients were recruited from 2 specialized hospitals, which may have introduced selection bias due to potential differences in disease severity or demographic characteristics. Crucially, only 23.2% of the patients with CC initially identified possessed the necessary MRI data for modeling inclusion. This limitation stemmed primarily from the requirement for locally archived DICOM files. Importantly, most excluded patients underwent initial MRI scans at external institutions—predominantly primary care hospitals where only non-digitalized paper reports were available, precluding DICOM file retrieval. While subtle technical variations in MRI scanners or protocols across institutions may exist, their impact on model generalizability is likely minor compared to the dominant limitation of sample size scarcity. Second, the relatively short follow-up period (<5 y) restricted our ability to evaluate long-term outcomes. This limitation may hinder the detection of delayed disease recurrence and treatment-related effects, especially in chronic disease contexts where complications can take several years to manifest. Third, although DL models are widely recognized for their superior performance when applied to large datasets, this study exhibited suboptimal performance due to the limited sample size and parameter constraints.

To bridge these gaps and advance clinical translation, 3 strategic priorities emerge. First, we will expand the cohort through multicenter collaborations targeting 1200+ patients—deliberately enriching underrepresented positive cases—while extending follow-up beyond 5 years to capture long-term outcomes. Second, transfer learning will be implemented to develop advanced 3D CNN architectures using volumetric DICOM data, overcoming current sample size barriers and unlocking DL’s latent potential. Third, recognizing that coronal T2WI—while optimal for bilateral tumor assessment—represents only one facet of MRI’s diagnostic capability, we will implement multidimensional sequence analysis. Concurrently, we will systematically integrate axial/sagittal planes and functional sequences (DWI, DCE-MRI) to refine microinvasion detection, ultimately validating this optimized system through prospective trials assessing its impact on surgical planning and survival outcomes.

In conclusion, this study used ML techniques to develop diagnostic and prognostic models using clinical and MRI data for patients with FIGO stage IA1 (LVSI+) to IIB CC. These models were designed to detect preoperative parametrial invasion and lymph node metastasis and to predict postoperative survival and recurrence. We trained and externally validated the models using data from 250 patients across 2 tertiary hospitals. In both diagnostic and prognostic applications, integrated ML models, particularly weighted KNN, demonstrated favorable performance with clinically applicable sensitivities. These findings suggest that integrated ML models may offer a valuable tool for improving preoperative staging and individualized prognosis prediction in CC, supporting more personalized and precise treatment strategies.