Comparison of circulating oxidative damage products
In sensitivity analyses using heteroscedasticity-consistent standard errors and adjusting for the number of vaginal deliveries and duration of vulvar leukoplakia, the overall group term remained significant for all 23 circulating biomarker outcomes (all robust Wald p < 0.001). The adjusted group D versus group A contrasts were directionally consistent with the primary analyses for all outcomes. All seven oxidative-damage biomarkers differed across the recorded severity categories after false-discovery-rate correction. MDA increased from 3.38 ± 0.76 µmol/L in group A to 4.55 ± 1.08 µmol/L in group D, with Welch F(3,123.49) = 16.984, FDR-adjusted p < 0.001, and partial η2 = 0.292. The skewed distributions of 8-iso-PGF2α and lipid hydroperoxide were analyzed using the Kruskal–Wallis test, yielding H = 15.326, FDR-adjusted p = 0.002, ε2 = 0.055, and H = 14.500, FDR-adjusted p = 0.002, ε2 = 0.051, respectively. The remaining oxidative-damage biomarkers also differed significantly after correction, and all adjusted ordered-trend tests were significant (p < 0.001) (Figure 2A–G and Supplementary Table 1).
Comparison of circulating antioxidant defense indicators
Antioxidant defense indicators were generally lower in the higher-severity categories. SOD decreased from 90.60 ± 17.80 U/mL in group A to 75.30 ± 17.10 U/mL in group D (F = 8.769, FDR-adjusted p < 0.001, partial η2 = 0.105). The GSH/GSSG ratio differed across groups (H = 27.663, FDR-adjusted p < 0.001, ε2 = 0.110), as did NADPH oxidase activity (H = 13.838, FDR-adjusted p = 0.003, ε2 = 0.048). All eight antioxidant/redox outcomes remained significant after correction, and all adjusted ordered-trend tests were significant (p < 0.001) (Figure 3A–H and Supplementary Table 1).
Comparison of extracellular matrix degradation and remodeling indicators
All eight circulating ECM-related biomarkers differed across the recorded severity categories after false discovery rate correction. The MMP-9/TIMP-1 ratio increased from 1.88 ± 0.54 in group A to 3.12 ± 0.95 in group D, with Welch F(3,122.21) = 29.225, FDR-adjusted p < 0.001, and partial η2 = 0.418. MMP-9 also showed a substantial group difference (Welch F(3,123.70) = 20.396, FDR-adjusted p < 0.001, partial η2 = 0.331). TIMP-1 and TIMP-2 were lower in the higher recorded categories (FDR-adjusted p = 0.003 for both), whereas desmosine/isodesmosine differed by the Kruskal–Wallis test (H = 21.976, FDR-adjusted p < 0.001, ε2 = 0.085). All adjusted ordered-trend tests were significant (p < 0.001) (Figure 4A–H and Supplementary Table 1).
Comparison of clinical phenotypes of vulvar leukoplakia
Differences were observed among the four groups in vulvar pruritus visual analog scale (VAS) score, pain/burning score, leukoplakic lesion area score, skin atrophy score, and fissuring/sclerosis score (p < 0.05). Among them, the difference in leukoplakic lesion area score was the most significant (H = 14.295, p = 0.003) (Table 3).
Comparison of POP-Q-related anatomical indicators
Finally, the authors also observed differences among the four groups in Aa, Ba, C, Ap, Bp, gh, and pb (p < 0.001). The change in point C best reflected the difference in POP severity, changing from -6.38 ± 0.95 cm in the no-prolapse group to 2.88 ± 1.36 cm in the severe prolapse group (F = 822.788, p < 0.001). In addition, gh increased with increasing prolapse severity, whereas pb showed a decreasing trend. It is worth noting that total vaginal length (TVL) did not differ significantly among the four groups (p = 0.244) (Table 4). Based on the most distal values of Ba, Bp, and C, group A included 19 anterior, 9 posterior, and 29 co-dominant cases; group B included 26 anterior, 3 posterior, and 28 co-dominant cases; group C included 28 anterior, 4 posterior, 11 apical, and 14 co-dominant cases; and all 57 patients in group D had an apical most-distal compartment (Supplementary Table 2).
DATA AVAILABILITY:
The complete recalculated statistical results, including balance diagnostics, assumption checks, effect sizes, pairwise comparisons, and adjusted sensitivity analyses, are provided in Supplementary Table 3.
FIGURES AND TABLES LEGENDS:

Figure 1: Retrospective patient screening and propensity score-matched grouping workflow. The flowchart shows the initial screening of patients with vulvar leukoplakia, exclusion according to the eligibility criteria, group allocation before propensity score matching (PSM), and final balanced allocation into group A, group B, group C, and group D after PSM. Please click here to view a larger version of this figure.

Figure 2: Circulating oxidative damage products across different degrees of pelvic organ prolapse. (A) Malondialdehyde (MDA); (B) 8-hydroxy-2'-deoxyguanosine (8-OHdG); (C) 8-iso-prostaglandin F2α (8-iso-PGF2α); (D) advanced oxidation protein products (AOPP); (E) protein carbonyl; (F) lipid hydroperoxide; (G) advanced glycation end products (AGEs). Please click here to view a larger version of this figure.

Figure 3: Circulating antioxidant defense and redox balance indicators across different degrees of pelvic organ prolapse. (A) Superoxide dismutase (SOD); (B) glutathione peroxidase (GSH-Px); (C) catalase (CAT); (D) total antioxidant capacity (T-AOC); (E) glutathione (GSH); (F) oxidized glutathione (GSSG); (G) GSH/GSSG ratio; (H) nicotinamide adenine dinucleotide phosphate (NADPH) oxidase activity. Please click here to view a larger version of this figure.

Figure 4: Extracellular matrix degradation and remodeling indicators across different degrees of pelvic organ prolapse. (A) Matrix metalloproteinase-2 (MMP-2); (B) matrix metalloproteinase-9 (MMP-9); (C) tissue inhibitor of metalloproteinase-1 (TIMP-1); (D) tissue inhibitor of metalloproteinase-2 (TIMP-2); (E) MMP-9/TIMP-1 ratio; (F) type I collagen C-terminal telopeptide/cross-linked C-terminal telopeptide of type I collagen (CTX-I/ICTP); (G) procollagen type III N-terminal peptide (PIIINP); (H) desmosine/isodesmosine. Please click here to view a larger version of this figure.
Table 1: Baseline characteristics of patients with vulvar leukoplakia and different degrees of pelvic organ prolapse after propensity score matching. This table summarizes demographic, reproductive, metabolic, and treatment-related baseline characteristics across the four propensity score-matched groups, together with between-group statistics and standardized mean differences. Please click here to download this Table.
Table 2: Assay characteristics and quality control parameters for circulating oxidative stress, antioxidant defense, and extracellular matrix remodeling indicators. This table summarizes the sample type, detection method, analytical range, detection limit, calibration linearity, and intra- and inter-batch variability for the measured circulating biomarkers. Please click here to download this Table.
Table 3: Clinical phenotype scores of vulvar leukoplakia across different degrees of pelvic organ prolapse. This table compares vulvar pruritus, pain or burning, lesion area, skin atrophy, and fissuring or sclerosis scores across the four groups. Please click here to download this Table.
Table 4: Pelvic organ prolapses quantification-related anatomical indicators across different degrees of pelvic organ prolapse. This table compares pelvic organ prolapse quantification anatomical measurements across the four groups, covering anterior, posterior, and apical support as well as genital hiatus, perineal body, and total vaginal length. Please click here to download this Table.
Supplementary Table 1: Detailed statistical analysis of circulating biomarkers across recorded pelvic organ prolapse severity categories. Group-wise biomarker values, overall statistical tests, FDR-adjusted p values, effect sizes, variance assessments, and adjusted ordered-trend analyses are presented. Please click here to download this file.
Supplementary Table 2: Distribution of POP-Q severity ranges and most distal pelvic compartments across recorded prolapse categories. The table summarizes leading-edge ranges and the distribution of anterior, posterior, apical, and co-dominant most-distal compartments for Groups A–D. Please click here to download this file.
Supplementary Table 3: Participant-level analytical dataset for the propensity score-matched cohort. Individual demographic and clinical characteristics, circulating oxidative stress, antioxidant/redox and ECM biomarkers, vulvar phenotype scores, and POP-Q anatomical measurements are provided for the 228 matched participants. Please click here to download this file.