$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
The final analytical cohort comprised 36 patients with persistent high-risk HPV infection following cervical cancer surgery, while baseline characteristics demonstrated general comparability between treatment groups (Table 1), a statistically significant age difference was noted (Interferon group: 56.0 ± 5.6 years vs. Paiteling group: 50.0 ± 5.0 years; p. = 0.005). As illustrated in Figure 2, participants were allocated to either Paiteling (n = 17) or interferon (n = 19) treatment groups through a standardized screening and enrollment process. The comprehensive treatment protocol for Paiteling administration across the three therapeutic stages is detailed in Figure 3A, while Figure 3B presents the comparative timeline of both treatment regimens with key assessment points.
Treatment efficacy demonstrated a striking temporal pattern, with a complete reversal of outcomes between the two groups during the follow-up period (Figure 4A). At the 6-month assessment, the Interferon group showed significantly higher cure rates (68.4%, 13/19) compared to the Paiteling group (17.6%, 3/17; 2= 7.646, p = 0.006). However, by the 12-month follow-up, this relationship reversed dramatically, with the Paiteling group achieving 94.1% clearance (16/17) versus 63.2% (12/19) in the Interferon group (2= 4.976, p. = 0.044). The number needed to treat analysis indicated that only 3.2 patients required Paiteling treatment to achieve one additional clearance compared to interferon. The viral load dynamics revealed distinct patterns between treatments (Figure 4B), with the Paiteling group exhibiting gradual but sustained viral suppression throughout the 12 months, while the Interferon group demonstrated an initial rapid response followed by viral rebound in 36.8% of cases. This pattern provides a mechanistic explanation for the observed long-term efficacy advantage of Paiteling treatment. Subgroup analysis revealed consistent treatment effects favoring Paiteling across all HPV genotypes (Figure 4C), with particularly strong effects observed for HPV 16 (100% vs 57.1% clearance) and HPV 18 (100% vs 66.7% clearance). The genotype-specific treatment effects were further supported by Firth’s penalized likelihood method, which demonstrated the consistent benefit of Paiteling across all subgroups.
Statistical analyses robustly supported these clinical findings. Univariable logistic regression (Table 2) identified treatment group as a significant predictor of HPV clearance (OR = 9.33, 95% CI: 1.40–187.02, p = 0.019), while multivariable analysis adjusting for age, baseline viral load, and HPV genotype (Table 3) maintained this strong association (adjusted OR = 9.52, 95% CI: 1.11–220.32, p. = 0.039). Importantly, by incorporating age into this multivariable model, it was confirmed that the superior therapeutic effect of Paiteling remained significant independent of the baseline age disparity. None of the covariates reached statistical significance in the adjusted model, confirming treatment group as the primary determinant of therapeutic outcome.
The multivariable prediction model for HPV clearance demonstrated moderate discriminative ability with an area under the ROC curve of 0.583 (95% CI: 0.35–0.82) (Figure 5A). The model achieved 63.6% overall accuracy, 75.0% sensitivity, and 33.3% specificity. Variable importance analysis (Figure 5B) confirmed treatment group as the strongest predictor in the model, consistent with the primary efficacy analyses. The calibration plot (Figure 5C) showed reasonable agreement between predicted probabilities and observed outcomes across the risk spectrum.
Safety monitoring documented a favorable tolerability profile for both interventions (Table 4). Adverse events occurred in 29.4% (5/17) of Paiteling recipients and 10.5% (2/19) of Interferon patients (p. = 0.226), with all events being mild to moderate in severity and resolving completely without treatment discontinuations. The most common adverse events included a local burning sensation, pruritus, and local pain, all of which were managed conservatively.
Comprehensive data quality assessment confirmed complete follow-up for all 36 enrolled patients, with no missing values for primary efficacy endpoints. The statistical approaches employed, including profile likelihood estimation and Firth’s correction for small sample limitations, ensured robust and reliable estimation of treatment effects despite the modest sample size. All analyses were conducted according to pre-specified statistical plans, maintaining methodological rigor and transparency throughout the reporting of results.
DATA AVAILABILITY:
The raw data supporting the conclusions of this study are provided in Supplementary File 1.

Figure 1: Representative Images of cervical surface following iodine staining. Colposcopic images demonstrating the appearance of cervical epithelium after iodine solution application. Yellow/White areas indicate typical iodine-negative regions corresponding to potential lesions. Please click here to view a larger version of this figure.

Figure 2: Schematic representation of the retrospective cohort design illustrating participant flow from screening to final analysis, including treatment allocation and follow-up assessments. Please click here to view a larger version of this figure.

Figure 3: Treatment protocol overview. (A) Schematic representation of the three-stage Paiteling administration protocol, including timing, concentration specifications, and application methodology. (B) Comparative treatment timeline of the two intervention regimens, highlighting key efficacy evaluation time points at 6 and 12 months. Please click here to view a larger version of this figure.

Figure 4: Comprehensive analysis of treatment efficacy and subgroup results. (A) Temporal patterns of HPV clearance rates demonstrating the significant efficacy reversal between Paiteling and Interferon groups from 6 to 12 months follow-up. (B) Longitudinal viral load dynamics showing distinct response trajectories, with Paiteling exhibiting sustained viral suppression and interferon showing initial response followed by viral rebound in 36.8% of cases. (C) Forest plot of genotype-specific subgroup analysis using Firth’s penalized likelihood method, illustrating consistent treatment benefit of Paiteling across all HPV genotypes with corresponding odds ratios and confidence intervals. Data are expressed as frequencies (percentages) for categorical outcomes. Where applicable, error bars represent the standard error of the mean (SEM) or 95% confidence intervals. Statistical significance is denoted as *p < 0.05 or p. < 0.01. Please click here to view a larger version of this figure.

Figure 5: Statistical modeling and predictive performance evaluation. (A) ROC curve analysis of the multivariable prediction model for HPV clearance showing moderate discriminative ability. (B) Variable importance plot from the multivariable logistic regression model demonstrating the relative contribution of each predictor variable, with treatment group identified as the strongest predictor. (C) Calibration plot showing the relationship between predicted probabilities and observed outcomes, indicating reasonable model calibration across the risk spectrum. Error bars in the calibration plot represent 95% confidence intervals. Please click here to view a larger version of this figure.
| Characteristic | Overall (N=36) | Interferon (N=19) | Paiteling Group (N=17) | p-value |
| Age (years) | 53.2 ± 6.1 | 56.0 ± 5.6 | 50.0 ± 5.0 | 0.005 |
| Baseline Viral Load (pg/mL) | 23.7 ± 6.4 | 24.1 ± 6.1 | 23.2 ± 6.8 | 0.7 |
| HPV Genotype | | | >0.9 |
| 16 | 13 (36%) | 7 (37%) | 6 (35%) | |
| 18 | 11 (31%) | 6 (32%) | 5 (29%) | |
| Other | 12 (33%) | 6 (32%) | 6 (35%) | |
| 12-month HPV Clearance | 28 (78%) | 12 (63%) | 16 (94%) | 0.044 |
Table 1: Baseline characteristics of study participants. Demographic and clinical characteristics of the 36 patients with persistent high-risk HPV infection following cervical cancer surgery, stratified by treatment group. Continuous variables are presented as mean ± standard deviation (SD); categorical variables as counts (percentages). P-values were calculated using t-tests for continuous variables and 2 tests or Fisher’s exact tests for categorical variables.
| Variable | Term | Odds Ratio | 95% CI | p-value |
| Treatment | Paiteling Group | 9.33 | 1.40-187.02 | 0.019 |
| Age | Per year increase | 0.92 | 0.79-1.05 | 0.223 |
| Baseline Viral Load | Per unit increase | 1.02 | 0.90-1.16 | 0.796 |
| HPV Genotype 18 | vs HPV 16 | 1.35 | 0.18-12.09 | 0.769 |
| HPV Genotype Other | vs HPV 16 | 0.9 | 0.13-5.99 | 0.91 |
Table 2: Univariable analysis for predictors of HPV clearance. Unadjusted logistic regression analysis examining the association between each predictor variable and HPV clearance at the 12-month follow-up. Odds ratios (OR) with 95% confidence intervals (CI) and p-values are presented for the treatment group and all covariates, demonstrating the treatment group as a significant predictor. To ensure robust estimation in this modest sample size, 95% CIs were estimated using the profile likelihood method, and corresponding p-values were derived from likelihood ratio tests.
| Variable | Term | Adjusted Odds Ratio | 95% CI | p-value |
| Treatment | Paiteling Group | 9.52 | 1.11-220.32 | 0.039 |
| Age | Per year increase | 0.98 | 0.82-1.18 | 0.866 |
| Baseline Viral Load | Per unit increase | 1.04 | 0.89-1.23 | 0.581 |
| HPV Genotype 18 | vs HPV 16 | 1.66 | 0.18-18.97 | 0.655 |
| HPV Genotype Other | vs HPV 16 | 0.96 | 0.10-9.42 | 0.973 |
Table 3: Multivariable analysis for predictors of HPV clearance. Adjusted logistic regression analysis examining the independent association between treatment group and HPV clearance after controlling for age, baseline viral load, and HPV genotype. Adjusted odds ratios (aOR) with 95% confidence intervals and p-values demonstrate the robustness of the treatment effect. Consistent with univariable analyses, 95% CIs were calculated via. the profile likelihood method, with p-values obtained from likelihood ratio tests.
| Adverse Reactions | Paiteling Group (n=17) | Interferon Group (n=19) | χ² | p-value |
| Local Burning Sensation | 1 (5.9%) | 0 (0%) | 1.15 | 0.284 |
| Itching | 1 (5.9%) | 1 (5.3%) | 0.007 | 0.935 |
| Local Pain | 2 (11.8%) | 1 (5.3%) | 0.496 | 0.479 |
| Mild Joint Aches | 1 (5.9%) | 0 (0%) | 1.15 | 0.284 |
| Total | 5 (29.4%) | 2 (10.5%) | 2.043 | 0.153 |
Table 4: Safety profile and adverse events comparison. Safety assessment documenting the frequency, types, and severity of adverse events in both treatment groups. All reported adverse events were mild to moderate in severity, with no serious adverse events or treatment discontinuations observed. Statistical comparisons were performed using2 tests or Fisher’s exact tests.
Supplementary File 1: Raw data supporting the conclusions of this study.Please click here to download this file.