$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
As shown in Table 1, there were no significant differences in demographic characteristics or biometric parameters between the MG and control groups. Crucially, the Standardized Mean Difference (SMD) for essential baseline covariates, such as axial length (AL) (SMD = 0.232) and age (SMD = 0.298), was well below the 0.3 threshold, establishing robust balance in these measured baseline variables between the cohorts. However, significant anatomical differences were observed in the UBM parameters (Figure 1). The MG group presented significantly decreased sCLD (unadjusted P < 0.001, SMD = 2.364), increased LV (unadjusted P < 0.001, SMD = 1.596), shallower ACD (unadjusted P = 0.0007, SMD = 1.050), and smaller TCPA (unadjusted P = 0.0037, SMD = 0.882). For multiple comparisons of UBM parameters, a False Discovery Rate (FDR) correction was applied, and the significance of sCLD, LV, ACD, and TCPA remained stable (all adjusted P < 0.05). The SMD values for all these parameters exceeded the threshold of 0.8, indicating large effect sizes and clinically meaningful differences between groups.
To rigorously account for parameter interdependence and isolate independent anatomical predictors, multivariable logistic regression modeling adjusting for baseline AL was performed (Table 2). The independent statistical significance of a shortened sCLD (P < 0.001), an increased LV (P < 0.001), and a narrowed TCPA (adjusted P = 0.0099) remained highly robust, demonstrating that these localized configurations are associated with MG development independently after statistical adjustment for axial length. Consequently, AL was included along with sCLD, LV, ACD, and TCPA in subsequent analyses to observe interrelationships.
To assess intergroup differences in ocular parameter coordination, correlation analyses were performed. Although none of these correlation differences reached statistical significance after FDR correction for multiple comparisons (all adjusted P > 0.05), a consistent trend of weakened correlations was observed in the MG group. Notable changes included: AL-ACD (control r = 0.718 vs. MG r = 0.334; r = -0.384); TCPA-LV (control r = -0.609 vs. MG r = -0.245; r = +0.364); TCPA-AL (control r = 0.531 vs. MG r = 0.289; r = -0.242); and sCLD-TCPA (control r = 0.224 vs. MG r = 0.017; r = -0.207). Scatter plots (Figure 2) visually confirmed this structural discoordination, showing greater data dispersion in the MG group, indicating a trend of preoperative weakening of physiological coupling among key ocular parameters in eyes that subsequently developed MG.
PCA was employed as an exploratory tool to identify multivariate anatomical patterns. The first two principal components (PCs) accounted for 79.0% of the total variance (PC1: 62.2%; PC2: 16.8%). PC1, which accounted for the largest proportion of variance (62.2%), was characterized by high negative loadings from LV and positive loadings from ACD and sCLD, suggesting a pattern interpretable as anterior segment crowding. PC2 was primarily defined by AL (loading = 0.741). The score plot (Figure 3A) demonstrated that the MG group clustered significantly toward the negative end of PC1 (t = 5.563, P < 0.001, Cohen’s d = 1.606) and exhibited higher scores on PC2 (t = -2.706, P = 0.010, Cohen’s d = 0.781) compared to controls. The loading plot (Figure 3B) visually illustrates the opposing contributions of LV, ACD, and sCLD to PC1. Overall, exploratory principal component analysis supported a dominant pattern of anterior segment crowding configuration in this malignant glaucoma cohort, characterized by a synergistic pattern of elevated LV with reduced ACD and sCLD.
Based on the hypothesized ciliary block mechanism and the anatomical patterns revealed by PCA, three exploratory composite risk scores were constructed. The ciliary block score integrates the risk values of sCLD, TCPA, and LV, reflecting the spatial relationship between the ciliary body and lens, representing a hypothesized anatomical risk substrate rather than a proven causal mechanism. The anterior segment crowding score combines the risk values of the sCLD, LV, and ACD, reflecting the configuration feature of anterior segment space crowding. The comprehensive risk score was the sum of the risk values of all five parameters (sCLD, TCPA, LV, ACD, and AL). Compared with the comprehensive risk score (Cohen’s d = 1.525) and the anterior segment crowding score (Cohen’s d = 1.856), the ciliary block score demonstrated the largest effect size (Cohen’s d = 1.995, P < 0.001) for group separation in this cohort. It is emphasized that these composite scores are sample-specific, unweighted exploratory constructs derived from z-scores and explicitly lack external clinical validation. These differential effect sizes were visually corroborated by the parallel coordinates plot (Figure 4), which visualized the multidimensional profiles of individual eyes across the five standardized parameters. The trajectories of the MG group (orange lines) consistently exhibited a pattern of elevated LV coupled with reduced ACD, TCPA, and sCLD, graphically illustrating the strong effects captured by both the ciliary block and anterior segment crowding scores. In contrast, the weaker discriminative power of the comprehensive risk score aligns with the less distinct separation in axial length observed between the groups in the plot.
In a highly exploratory, post-hoc analysis, five extreme high-risk cases with the highest ciliary block scores from the MG group were selected and compared with five control cases exhibiting the lowest scores. As detailed in Table 3, all high-risk MG cases demonstrated characteristic ciliary body-lens crowding, featuring shortened sCLD (-0.199 to -0.095 mm), narrowed TCPA (34.1–42.4°), and elevated LV (0.98–1.35 mm), contrasting sharply with normal parameters in low-risk controls. Analysis of these extreme cases revealed potential anatomical phenotypes (Table 3), predominantly categorized as a lens anterior displacement dominant pattern (n = 4), which was characterized by markedly elevated LV (1.33 mm ± 0.03 mm) serving as the main risk driver, alongside a single unique case (n = 1) exhibiting an extreme ciliary body proximity pattern. Across all five extreme high-risk cases, the revised overall composite score was 4.20 ± 0.79. Because this subgroup is too small to support broader phenotype claims, any conclusions drawn from this 5-case extreme subset remain highly provisional and strictly localized to this sample.
Using a computer-based randomization procedure, UBM images from 10 patients were randomly selected without replacement from the total study cohort via simple random sampling. To ensure randomization of the measured quadrants, this study independently assigned a random measurement quadrant for CBT0, CBT1000, CBTmax, APCB, TCA, and sCLD. Following a dual-stratified framework, global ocular metrics for general inter-group comparisons were derived by averaging measurements acquired across all four anatomical quadrants, whereas for the reliability analysis, a single quadrant was randomly allocated. An experienced glaucoma specialist, blinded to patient group assignment and initial measurement results, performed repeated measurements on the selected images. As shown in Table 4, the evaluated parameters included ACD, CCD, LV, CBT0 in the nasal quadrant (CBT0N), CBT1000 in the inferior quadrant (CBT1000I), CBTmax in the inferior quadrant (CBTmaxI), APCB in the inferior quadrant (APCBI), TCPA in the temporal quadrant (TCPAT), and sCLD in the temporal quadrant (sCLDT). The ICC values for all the parameters were above 0.90, demonstrating excellent interobserver reliability and reproducibility of the measurement protocol.
In summary, preoperative metrics indicate that eyes exhibiting an association with post-surgical malignant glaucoma possess distinct morphological profiles characterized by an independent cluster of narrow ciliary-lenticular space (sCLD), elevated lens protrusion (LV), and a compressed ciliary complex (TCPA), accompanied by an exploratory trend of structural discoordination among biometric variables.
DATA AVAILABILITY:
The fully de-identified raw datasets are provided in Supplementary File 1.

Figure 1: UBM measurement parameters of the anterior segment. (A) illustrates the measurements of anterior chamber depth (ACD), lens vault (LV), and ciliary process-ciliary process distance (CCD) on a representative horizontal panoramic scan matrix. (B) quantifies various localized ciliary body parameters, including ciliary body thickness at the scleral spur (CBT0), maximum ciliary body thickness (CBTmax), ciliary body thickness 1000 m posterior to the scleral spur (CBT1000), anterior placement of the ciliary body (APCB), trabecular-ciliary process angle (TCPA), and simulated cilio-lenticular distance (sCLD) on a representative radial quadrant scan recess. Scale bar = 1 mm. ACD, anterior chamber depth; LV, lens vault; CCD, ciliary process-ciliary process distance; CBT0, ciliary body thickness at the scleral spur; CBTmax, maximum ciliary body thickness; CBT1000, ciliary body thickness 1000 m posterior to the scleral spur; APCB, anterior placement of ciliary body; TCPA, trabecular-ciliary process angle; sCLD, simulated cilio-lenticular distance; SS, scleral spur. Scale bars: 1 mm. Please click here to view a larger version of this figure.

Figure 2: Comparison of correlation patterns between the MG and control groups for key ocular parameter pairs. Scatter plots illustrate the linear relationships between selected ocular parameters with fitted regression lines and 95% confidence intervals: (A) AL vs. ACD, (B) TCPA vs. LV, (C) AL vs. TCPA, and (D) sCLD vs. TCPA. Control group data (blue) demonstrated correlation coefficients of r = 0.718, –0.609, 0.531, and 0.224 for panels A–D, respectively, whereas the MG group (red) showed corresponding r values of 0.334, –0.245, 0.289, and 0.017. The between-group differences in correlation coefficients (Δr = –0.207 to –0.385) suggest altered inter-parameter associations in MG eyes. MG, malignant glaucoma; AL, axial length; ACD, anterior chamber depth; TCPA, trabecular-ciliary process angle; LV, lens vault; sCLD, simulated cilio-lenticular distance. Please click here to view a larger version of this figure.

Figure 3: Principal component analysis of ocular parameters in PACG patients with postoperative malignant glaucoma. (A) PCA score plot showing the distribution of MG (purple) and control (blue) groups in the principal component space. The MG group was significantly distributed toward the negative region of PC1 (p < 0.001) and showed significantly higher scores on PC2 (p = 0.010) compared to controls. (B) PCA loading plot illustrating the contributions of key ocular parameters to the first two principal components. Point size indicates representation quality (cos2), and color denotes contribution level (gray: <20%; orange: 20-40%; dark orange: >40%). PC1 and PC2 collectively explain 79.0% of total variance (PC1: 62.2%; PC2: 16.8%), with PC1 reflecting anterior segment crowding and PC2 primarily associated with axial length variation. MG, malignant glaucoma; LV, lens vault; ACD, anterior chamber depth; sCLD, simulated cilio-lenticular distance; TCPA, trabecular-ciliary process angle; AL, axial length. Please click here to view a larger version of this figure.

Figure 4: Parallel coordinates plot showcasing multidimensional ocular profiles. This chart visualizes the comprehensive geometric trajectories of individual eyes from the Control (blue lines, n=24) and Malignant Glaucoma (MG, orange lines, n = 24) cohorts across five biometric parameters (AL, ACD, LV, TCPA, and sCLD). Crucially, the "Standardized Value" on the Y-axis represents the sample-specific min-max scaled values [0, 1] derived directly from standardized z-scores to eliminate scale-dependent variance. The trajectories graphically demonstrate that the pre-surgical risk substrate of the MG group is fundamentally characterized by an independent cluster of elevated LV tightly coupled with compressed axial, angular, and cilio-lenticular spaces. Please click here to view a larger version of this figure.
| Characteristic | Control (n=24) | Malignant Glaucoma (n=24) | P-value | SMD |
| Age, years (Median [IQR]) | 65.00 [60.00, 69.00] | 61.00 [54.75, 66.00] | 0.0984 | 0.298 |
| Sex, female, n (%) | 19 (79.2%) | 22 (91.7%) | 0.4158 | 0.36 |
| Eye, OD, n (%) | 12 (50.0%) | 15 (62.5%) | 0.5612 | 0.254 |
| PAS, <180°, n (%) | 11 (45.8%) | 9 (37.5%) | 0.7702 | 0.17 |
| IOP, mmHg (Mean ± SD) | 24.14 ± 4.90 | 24.46 ± 4.87 | 0.8189 | 0.066 |
| AL, mm (Mean ± SD) | 21.67 ± 0.61 | 21.52 ± 0.69 | 0.4257 | 0.232 |
| ACD, mm (Mean ± SD) | 2.03 ± 0.26 | 1.75 ± 0.28 | <0.001 | 1.05 |
Table 1: Baseline characteristics of study participants. Data are presented as mean ± standard deviation, median (interquartile range), or number (percentage) as appropriate. P-values were calculated using Student’s t-test for normally distributed continuous variables, Mann-Whitney U test for non-normally distributed continuous variables, and Chi-square test or Fisher’s exact test for categorical variables. To account for multiple comparisons across the 15 baseline characteristics, Benjamini-Hochberg false discovery rate (FDR) correction was applied to the raw P-values; adjusted P-values are reported in the “P. adjust” column, with P < 0.05 considered statistically significant. SMD, standardized mean difference; OD, oculus dexter (right eye); PAS, peripheral anterior synechiae; AL, axial length; IOP, intraocular pressure; ACD, anterior chamber depth; CBT0, ciliary body thickness at the scleral spur; CBTmax, maximum ciliary body thickness; CBT1000, ciliary body thickness 1000 µm posterior to the scleral spur; APCB, anterior placement of ciliary body; TCPA, trabecular-ciliary process angle; sCLD, simulated cilio-lenticular distance; LV, lens vault; CCD, ciliary process-ciliary process distance. Significant differences (P. adjust < 0.05, SMD > 0.8) were observed in ACD, TCPA, sCLD, and LV between groups.
| Parameter | Control (n=24) | Malignant Glaucoma (n=24) | Unadjusted P-value | Adjusted P-value (vs AL) |
| LV, mm | 0.834 ± 0.151 | 1.081 ± 0.158 | <0.001 | <0.001 |
| CCD, mm | 9.474 ± 0.316 | 9.612 ± 0.310 | 0.1344 | 0.0877 |
| CBT0, mm | 0.851 ± 0.078 | 0.854 ± 0.089 | 0.9002 | 0.7691 |
| CBTmax, mm | 0.938 ± 0.089 | 0.926 ± 0.089 | 0.6516 | 0.7637 |
| CBT1000, mm | 0.464 ± 0.082 | 0.471 ± 0.067 | 0.7494 | 0.5807 |
| APCB, mm | 0.629 ± 0.195 | 0.631 ± 0.133 | 0.9636 | 0.8426 |
| TCPA, ° | 53.822 ± 9.442 | 45.824 ± 8.685 | 0.0037 | 0.0099 |
| sCLD, mm | 0.118 ± 0.063 | -0.052 ± 0.080 | <0.001 | <0.001 |
| LV/ACD | 0.422 ± 0.120 | 0.645 ± 0.192 | <0.001 | <0.001 |
| ACD/AL | 0.094 ± 0.010 | 0.081 ± 0.012 | <0.001 | 0.0039 |
Table 2: Two-Sample t-test and multivariable logistic regression analysis of preoperative Ultrasound Biomicroscopy (UBM) parameters. This table presents the unadjusted baseline comparisons alongside the independent statistical parameters derived after multivariable adjustment for axial length (AL). Data are presented as mean ± standard deviation. sCLD, simulated cilio-lenticular distance; LV, lens vault; ACD, anterior chamber depth; TCPA, trabecular-ciliary process angle; CCD, ciliary process-ciliary process distance; CBT0, ciliary body thickness at the scleral spur; CBTmax, maximum ciliary body thickness; CBT1000, ciliary body thickness 1000 m posterior to the scleral spur; APCB, anterior placement of the ciliary body.
| Case ID | Risk Pattern | sCLD (mm) [Risk Score] | TCPA (°) [Risk Score] | LV (mm) [Risk Score] | Exploratory Composite Score | Dominant Risk Factors |
| MG_HighRisk_1 | 1 | -0.109 [1.27] | 34.1 [1.60] | 1.35 [1.98] | 4.85 | LV risk (1.98) + TCPA risk (1.60) |
| MG_HighRisk_2 | 1 | -0.095 [1.15] | 35.9 [1.41] | 1.35 [1.98] | 4.54 | LV risk (1.98) + TCPA risk (1.41) |
| MG_HighRisk_3 | 1 | -0.132 [1.48] | 36.7 [1.33] | 1.29 [1.69] | 4.51 | LV risk (1.69) + sCLD risk (1.48) |
| MG_HighRisk_4 | 1 | -0.126 [1.42] | 39.6 [1.04] | 1.32 [1.82] | 4.27 | LV risk (1.82) + sCLD risk (1.42) |
| MG_HighRisk_5 | 2 | -0.142 [1.56] | 42.4 [0.78] | 0.98 [0.49] | 2.83 | sCLD risk (1.56) + TCPA risk (0.78) |
| Mean ± SD | — | -0.121 ± 0.019 | 37.7 ± 3.4 | 1.26 ± 0.16 | 4.20 ± 0.79 | (Overall Profile) |
Table 3: Comparative analysis of anatomical parameters and standardized risk scores in extreme high-risk and low-risk cases. This table presents a comprehensive comparison of anatomical parameters and standardized risk scores for extreme cases from the MG (high-risk) and control (low-risk) groups. Values in brackets represent standardized risk scores, calculated by Z-score transformation with direction adjustment so that positive values indicate higher MG risk. The Ciliary Block Score is the sum of sCLD, TCPA, and LV risk scores. Note the consistent pattern: high-risk cases show positive risk scores and abnormal anatomical parameters (negative sCLD, narrow TCPA, elevated LV), while low-risk controls show negative risk scores and normal anatomical parameters. MG, malignant glaucoma; sCLD, simulated cilio-lenticular distance; TCPA, trabecular-ciliary process angle; LV, lens vault.
| Parameter | ICC | Mean 1 | Mean 2 | Difference | 95% CI of Difference |
| ACD | 0.977 | 2.063 | 2.019 | 0.044 | (0.013, 0.075) |
| CCD | 0.908 | 9.845 | 9.763 | 0.082 | (-0.010, 0.175) |
| LV | 0.969 | 0.89 | 0.874 | 0.016 | (-0.020, 0.052) |
| CBT0N | 0.93 | 0.808 | 0.8 | 0.008 | (-0.006, 0.022) |
| CBT1000I | 0.96 | 0.515 | 0.523 | -0.008 | (-0.027, 0.011) |
| CBTmaxI | 0.946 | 0.933 | 0.954 | -0.021 | (-0.052, 0.010) |
| APCBI | 0.959 | 0.695 | 0.685 | 0.01 | (-0.043, 0.062) |
| TCPAT | 0.992 | 52.163 | 50.886 | 1.277 | (0.512, 2.042) |
| sCLDT | 0.974 | 0.098 | 0.094 | 0.004 | (-0.013, 0.021) |
Table 4: Interobserver reliability of parameters measured by UBM. ICC values > 0.90 are considered to represent excellent reliability. Mean 1 and Mean 2 represent the average values from the first and second observers, respectively. Difference represents the mean difference between observers (Mean 1 – Mean 2). ICC, intraclass correlation coefficient; CI, confidence interval; ACD, anterior chamber depth; CCD, ciliary process-ciliary process distance; LV, lens vault; CBT0N, ciliary body thickness at the scleral spur in the nasal quadrant; CBT1000I, ciliary body thickness at 1000 µm from the scleral spur in the inferior quadrant; CBTmaxI, maximum ciliary body thickness in the inferior quadrant; APCBI, anterior placement of the ciliary body in the inferior quadrant; TCPAT, trabecular-ciliary process angle in the temporal quadrant; sCLDT, simulated cilio-lenticular distance in the temporal quadrant; UBM, ultrasound biomicroscopy.
Supplementary File 1: The fully de-identified raw datasets. Please click here to download this file.