Diagnostic efficacy of dual-energy computed tomography-based fractal analysis for assessing extramural venous invasion/tumor deposits and peripheral nerve invasion in rectal cancer
Introduction
In recent years, with the popularization of rectal cancer (RC) screening modalities, such as computed tomography (CT) and magnetic resonance imaging (MRI), and the progress of diagnostic technology, the prevalence of RC has gradually increased, becoming one the most prevalent malignancies globally (1,2). Although the prognosis of RC is informed by tumor-node-metastasis (TNM) staging (3,4), several pathological features independent of the TNM system—including status of extramural venous invasion (EMVI), tumor deposits (TDs), and peripheral nerve invasion (PNI)—have garnered increased attention in guiding treatment strategies for patients with RC.
EMVI, defined histologically as the involvement of veins beyond the muscularis propria, may promote tumor dissemination along vascular structures and is recognized as an independent prognostic factor associated with an increased risk of local recurrence, distant metastasis, and reduced overall survival (5,6). TDs, defined as discrete tumor nodules lacking residual lymph node tissue, are also associated with poor prognosis in colorectal cancer, and this adverse prognostic effect persists independently of lymph node metastasis (7,8). PNI is the infiltration of tumor cells into and along nerve sheaths and is significantly associated with higher postoperative recurrence rates and poorer prognosis in patients with RC (9,10). Preoperative evaluation of EMVI, TD, and PNI allows for the more accurate stratification of patients with RC. This comprehensive approach enhances risk assessment, supports the development of individualized treatment strategies, and ultimately contributes to improved clinical outcomes.
MRI and CT are both recommended as preoperative imaging modalities in current clinical guidelines and are commonly used in combination with TNM staging. Although MRI is considered the gold standard for evaluating high-risk pathological features in patients with RC, its clinical utility is limited by contraindications and its inability to comprehensively assess distant metastases. Contrast-enhanced CT is a more readily applicable imaging modality in current clinical practice, as recommended by existing guidelines. Nevertheless, its inadequate soft tissue resolution complicates the assessment of risk factors, such as PNI, EMVI, and TD in the tumor periphery. Compared with conventional CT, dual-energy CT (DECT) enhances diagnostic accuracy and offers quantitative data to evaluate the biological behavior of malignancies (11). DECT has demonstrated significant applicability in the assessment of colorectal cancer, particularly in refining T staging and forecasting EMVI and lymph node metastases (12-14).
Fractal analysis offers a reliable method for quantifying tumor heterogeneity via fractal dimension (FD) (15,16). FD is a non-integer metric that characterizes the inherent morphology of the tissue, with its magnitude indicating the intricacy of tumor heterogeneity (17). Elevated FD typically indicates increased tumor aggressiveness, and hence, fractal analysis is extensively employed in tumor risk assessment. Fractal analysis has been shown to predict the effectiveness of neoadjuvant radiation in locally advanced RC (18).
We hypothesized that fractal analysis of DECT can elucidate the biological attributes of high-risk factors associated with RC pathology. This study aimed to assess the diagnostic significance of quantitative parameters derived from fractal analysis with DECT in the preoperative prediction of EMVI/TD and PNI status. We present this article in accordance with the STARD reporting checklist (available at https://qims.amegroups.com/article/view/10.21037/qims-2025-862/rc).
Methods
Patient population
This study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments and was approved by the Ethics Committees of the Second Affiliated Hospital of Chongqing Medical University (No. 2020-366). Written informed consent was obtained from all participants. This study comprised a retrospective development phase and a prospective validation phase. Patients with histologically confirmed RC who underwent preoperative DECT between April 2021 and November 2024 were consecutively enrolled. Participants were divided into a retrospective development cohort (n=85) enrolled between April 2021 and October 2023 and a prospective independent validation cohort (n=45) enrolled between October 2023 and November 2024. The same exclusion criteria were applied to both cohorts: (I) pathologically confirmed benign lesions or non-adenocarcinomas (n=28); (II) administration of anticancer treatment prior to CT examination (n=56); (III) inadequate visualization of tumor lesions (n=12); (IV) presence of distant metastases (n=8); and (V) absence of surgical or histological examination within 2 weeks after CT examination (n=11) (Figure 1).
Histological examination
A board-certified gastrointestinal pathologist (>5 years of experience) diagnosed EMVI [tumor within blood vessels beyond muscularis propria confirmed by immunohistochemical (IHC) tests for CD31/CD34], TD (discrete tumor foci in the perirectal fat without regional lymph node or identifiable vascular or neural structures), and PNI (tumor encircling ≥33% of nerve circumference verified by IHC test for S100) according to the International Collaboration on Cancer Reporting (ICCR) guidelines (19). All specimens underwent 10% formalin fixation (24–48 h), 4-µm sectioning, and hematoxylin and eosin staining, with equivocal cases resolved through IHC and senior pathologist review.
DECT protocol and reconstruction
All CT exams were conducted with a DECT SOMATOM Force system (Siemens Healthineers, Erlangen, Germany). The scanning parameters included tube voltages of 100 and 140 kV, a tube current of 350 mA, a pitch of 1, and a gantry rotation time of 0.5 seconds. The contrast agent (370 mg/mL; Ultravist, Schering, Berlin, Germany) was administered into an elbow vein with a dual-head power injector at a prescribed dosage of 1.5 mL/kg and a flow rate of 3 mL/s. The scan commenced automatically in the arterial phase after a 15-second wait, which was triggered by the abdominal aortic threshold reaching 120 Hounsfield units (HU).
The DECT data underwent postprocessing on a commercial 3D multimodality workstation (syngo.via v. VA30A, Siemens Healthineers) with a customized soft tissue convolution kernel (Qr40, Siemens Healthineers) and an iterative reconstruction technique at a strength level of 3 [advanced modeled iterative reconstruction (ADMIRE); Siemens Healthineers]. Standard linear mixed pictures were thereafter created automatically with a mixing ratio of 0.6 (M0.6), with the integration of 60% of the data from a low tube voltage of 100 keV and 40% from a high tube voltage of 140 keV (20). Arterial-phase images with a 1-mm slice thickness were imported into the workstation and displayed in dual-energy mode to automatically generate conventional mixed-energy images, virtual monoenergetic maps, iodine concentration maps, and effective atomic number (Zeff) maps. This study substituted conventional mixed-energy CT pictures with regular linear blend images, as indicated by previous research (21).
Fractal analysis
Fractal analysis was performed with ImageJ software v. 1.54d (National Institutes of Health, Bethesda, MD, USA; http://rsb.info.nih.gov/ij) along with the FracLac plugin version 2.5 (http://rsb.info.nih.gov/ij/plugins/fraclac/FLHelp/Introduction.htm). FracLac computed FD by employing the conventional box-counting technique on images obtained from 12 distinct grid positions. Fractal analyses were conducted by two radiologists (C.Z. and H.W., with 5 and 3 years of abdominal CT experience, respectively) who were blinded to each other’s data and any additional clinical information. C.Z. conducted the fractal analysis of the total datasets, whereas H.W. analyzed the development cohort. The tumor was identified on 40-keV monoenergetic pictures, with the maximal tumor level assessed in the sagittal orientation. Region of interest (ROI) delineation was conducted in the axial position, with the focus on areas of notable lesion enhancement and with prominent blood vessels and the intestinal lumen being excluded (Figure 2). For the standard linear blend image (M0.6), the window width and window level were configured to 400 and 40 HU, respectively. These were subsequently imported to ImageJ software and converted to the 8-bit type image. The window width and window level for the iodine concentration images and the Zeff image were configured to 400 and 40 HU, respectively. The color contrast for the iodine maps (window width, 200; window level, 0) and the Zeff image (window width, 8; window level, 7) was adjusted, and the images were then imported to ImageJ to change the image type to 8-bit color. The minimal box size was established at 2 pixels and progressively augmented throughout the sample duration until it attained 45% of the maximum total area specified (22). After window-width/level optimization, the ROI delineated on the 40-keV monoenergetic image was propagated to the conventional mixed-energy images, iodine concentration images, and Zeff images for the extraction of FD parameters (FD-Con, FD-IC, and FD-Zeff, respectively).
Statistical analysis
The Chi-squared test or Fisher exact test was employed for categorical variables. Normally distributed continuous data are presented as the mean ± standard deviation and were analyzed with the t-test; nonnormally distributed continuous variables are presented as medians and quartiles [median (interquartile range)] and were analyzed with the Mann-Whitney test. Bland-Altman plots and the intraclass correlation coefficient (ICC) were employed to assess the concordance between two radiologists evaluating FD. Receiver operating characteristic (ROC) curves were constructed, and the area under the curve (AUC) and Youden J index were computed. The DeLong test was employed to compare the AUCs of various approaches. Univariate and multivariate logistic regression analyses were employed to identify independent risk factors for PNI and EMVI/TD status. Statistical analysis was conducted via MedCalc v. 19.2 (MedCalc Software, Ostend, Belgium) and GraphPad Prism v. 9.4.1 (Dotmatics, Boston, MA, USA), with P<0.05 indicating statistical significance.
Results
Clinicopathological characteristics
The study included 130 patients, with 85 and 45 in the development and validation cohorts, respectively. Furthermore, in the development cohort, the PNI subgroup had 37 positive and 48 negative individuals, whereas the EMVI/TD subgroup included 19 positive and 66 negative patients. In the validation cohort, the PNI subgroup had 22 positive and 23 negative patients, whereas the EMVI/TD subgroup included 16 positive and 29 negative patients. In the EMVI/TD subgroup, the percentage of patients with carcinoembryonic antigen levels beyond 5 ng/mL was significantly greater in the positive group than in negative group in both the development and validation cohorts (P=0.038 and P=0.007, respectively). No significant difference was observed between the PNI-positive and PNI-negative groups (Table 1).
Table 1
| Characteristics | Development cohort (n=85) | Validation cohort (n=45) | |||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| PNI | EMVI/TD | PNI | EMVI/TD | ||||||||||||
| Positive (n=37) | Negative (n=48) | P | Positive (n=19) | Negative (n=66) | P | Positive (n=22) | Negative (n=23) | P | Positive (n=16) | Negative (n=29) | P | ||||
| Age (years) | 65 [57, 70] | 65 [53, 68] | 0.385 | 67 [60, 70] | 64 [53, 69] | 0.415 | 66 [53, 71] | 61 [56, 64] | 0.242 | 60±11 | 65±11 | 0.208 | |||
| Gender | 0.947 | 0.932 | 0.601 | 0.578 | |||||||||||
| Male | 19 | 25 | 10 | 34 | 13 | 12 | 8 | 17 | |||||||
| Female | 18 | 23 | 9 | 32 | 9 | 11 | 8 | 12 | |||||||
| Clinical T stage | 0.257 | 0.097 | 0.609 | 0341 | |||||||||||
| T1–T2 | 11 | 20 | 10 | 21 | 7 | 9 | 4 | 12 | |||||||
| T3–T4 | 26 | 28 | 9 | 45 | 15 | 14 | 12 | 17 | |||||||
| Clinical N stage | 0.452 | 0.840 | 0.295 | 0.256 | |||||||||||
| N0 | 17 | 26 | 9 | 33 | 9 | 13 | 6 | 16 | |||||||
| N1–N2 | 20 | 22 | 10 | 33 | 13 | 10 | 10 | 13 | |||||||
| Tumor length (cm) | 4.544±2.080 | 4.673±1.738 | 0.757 | 3.934±1.115 | 4.813±2.016 | 0.073 | 4.275±1.728 | 4.363±1.943 | 0.873 | 4.767±1.840 | 4.074±1.793 | 0.226 | |||
| Tumor thickness (cm) | 1.288±0.517 | 1.199±0.418 | 0.380 | 1.361±0.362 | 1.203±0.485 | 0.192 | 1.322±0.479 | 1.201±0.383 | 0.352 | 1.230 [0.925, 1.875] | 1.210±0.407 | 0.356 | |||
| CEA (ng/mL) | 0.144 | 0.038 | 0.302 | 0.007 | |||||||||||
| ≤5 | 24 | 38 | 10 | 52 | 11 | 15 | 5 | 21 | |||||||
| >5 | 13 | 10 | 9 | 14 | 11 | 8 | 11 | 8 | |||||||
| CA19-9 (ng/mL) | 0.396 | 0.122 | >0.99 | 0.692 | |||||||||||
| ≤37 | 33 | 46 | 16 | 63 | 18 | 19 | 14 | 23 | |||||||
| >37 | 4 | 2 | 3 | 3 | 4 | 4 | 2 | 6 | |||||||
Data are presented as mean ± standard deviation, number, or median [interquartile range]. CA19-9, carbohydrate antigen 19-9; CEA carcinoembryonic antigen; EMVI, extramural venous invasion; N, node; PNI, peripheral nerve invasion; T, tumor; TD, tumor deposit.
Interobserver concordance of FD parameters
Excellent interobserver agreement was confirmed for all FD parameters in the development cohort, which was supported by both the ICC and Bland-Altman analyses. The ICCs for FD-Con, FD-IC, and FD-Zeff were 0.973 [95% confidence interval (CI): 0.959–0.982], 0.971 (95% CI: 0.955–0.981), and 0.965 (95% CI: 0.947–0.971), respectively. Concordantly, the Bland-Altman plots revealed minimal bias with mean differences of 0.004, −0.002, and −0.003 and narrow limits of agreement for FD-Con, FD-IC, and FD-Zeff, respectively (Figure 3).
Analysis of fractal parameters
Table 2 summarizes the distribution of FD values in the two subgroups. In the PNI subgroup, the positive group, as compared to the negative group, had a significantly higher mean FD-IC [development cohort: 0.986±0.118 vs. 0.865 (0.796, 1.015), P=0.001; validation cohort: 0.985±0.105 vs. 0.914±0.093, P=0.019] and FD-Zeff (development cohort: 0.947±0.125 vs. 0.868±0.103, P=0.002; validation cohort: 1.002±0.110 vs. 0.890±0.130, P=0.003), while the FD-Con was not significantly different. In the EMVI/TD subgroup, the negative group, as compared to the positive group, had a significantly higher mean FD-IC (development cohort: 1.013±0.097 vs. 0.908±0.136, P=0.002; validation cohort: 1.014±0.087 vs. 0.913±0.096, P=0.001) and FD-Zeff (development cohort: 1.006±0.080 vs. 0.873±0.112, P<0.001; validation cohort: 1.048±0.097 vs. 0.887±0.113, P=0.001), while the FD-Con was not significantly different.
Table 2
| Parameter | PNI | EMVI/TD | ||||||
|---|---|---|---|---|---|---|---|---|
| Positive | Negative | P | Positive | Negative | P | |||
| Development cohort (n=85) | ||||||||
| FD-Con | 0.944±0.121 | 0.900±0.099 | 0.067 | 0.932±0.098 | 0.915±0.115 | 0.556 | ||
| FD-IC | 0.986±0.118 | 0.865 (0.796, 1.015) | 0.001 | 1.013±0.097 | 0.908±0.136 | 0.002 | ||
| FD-Zeff | 0.947±0.125 | 0.868±0.103 | 0.002 | 1.006±0.080 | 0.873±0.112 | <0.001 | ||
| Validation cohort (n=45) | ||||||||
| FD-Con | 0.972±0.097 | 0.946±0.115 | 0.415 | 0.971±0.089 | 0.952±0.115 | 0.571 | ||
| FD-IC | 0.985±0.105 | 0.914±0.093 | 0.019 | 1.014±0.087 | 0.913±0.096 | 0.001 | ||
| FD-Zeff | 1.002±0.110 | 0.890±0.130 | 0.003 | 1.048±0.097 | 0.887±0.113 | 0.001 | ||
Data are presented as mean ± standard deviation or median (interquartile range). EMVI, extramural venous invasion; FD, fractal dimension; FD-Con, FD from conventional mixed-energy images; FD-IC, FD from iodine concentration images; FD-Zeff, FD from effective atomic number images; PNI, peripheral nerve invasion; TD, tumor deposit.
Diagnostic efficacy of FD parameters
Table 3 and Figure 4 show the diagnostic performance of significant FD characteristics. In the PNI subgroup, the AUC with 95% CI, sensitivity, and specificity of FD-IC were 0.706 (95% CI: 0.595–0.816), 0.892, and 0.521 for the development cohort, respectively, while for the validation cohort, they were 0.732 (95% CI: 0.579–0.885), 0.727, and 0.783, respectively; meanwhile, these parameters for FD-Zeff were 0.710 (95% CI: 0.594–0.852), 0.649, and 0.750 in the development cohort, respectively, while they were 0.752 (95% CI: 0.630–0.912), 0.682, and 0.826 in the validation cohort, respectively. In the EMVI/TD subgroup, the AUC for FD-IC was 0.739 (95% CI: 0.629–0.849) in the development cohort, and the sensitivity and specificity were 0.895 and 0.561, respectively; while in the validation cohort, the AUC was 0.795 (95% CI: 0.662–0.929), with a sensitivity and specificity of 0.938 and 0.586, respectively; meanwhile, these parameters for FD-Zeff were 0.845 (95% CI: 0.758–0.931), 0.895, and 0.727 for the development cohort, respectively, while they were 0.855 (95% CI: 0.740–0.969), 0.750, and 0.862 for the validation cohort, respectively. Table 4 presents the variations in AUC observed for the FD parameters. In the PNI subgroup, the AUC of FD-Zeff was significantly higher than that of FD-Con in the validation cohort (P=0.044). In the EMVI/TD subgroup, both FD-IC and FD-Zeff exhibited significantly elevated AUCs compared to FD-Con (development cohort: P=0.008 and P<0.001; validation cohort: P=0.030 and P=0.002).
Table 3
| Parameter | PNI | EMVI/TD | |||||||
|---|---|---|---|---|---|---|---|---|---|
| AUC (95% CI) | Threshold | Sensitivity | Specificity | AUC (95% CI) | Threshold | Sensitivity | Specificity | ||
| Development cohort (n=85) | |||||||||
| FD-IC | 0.706 (0.595, 0.816) | 0.869 | 0.892 | 0.521 | 0.739 (0.629, 0.849) | 0.911 | 0.895 | 0.561 | |
| FD-Zeff | 0.710 (0.594, 0.825) | 0.916 | 0.649 | 0.750 | 0.845 (0.758, 0.931) | 0.918 | 0.895 | 0.727 | |
| Combined | 0.744 (0.638, 0.849) | – | 0.811 | 0.625 | 0.848 (0.766, 0.929) | – | 1.000 | 0.606 | |
| Validation cohort (n=45) | |||||||||
| FD-IC | 0.732 (0.579, 0.885) | – | 0.727 | 0.783 | 0.795 (0.662, 0.929) | – | 0.938 | 0.586 | |
| FD-Zeff | 0.752 (0.603, 0.902) | – | 0.818 | 0.696 | 0.855 (0.740, 0.969) | – | 0.750 | 0.862 | |
| Combined | 0.771 (0.630, 0.912) | – | 0.682 | 0.826 | 0.879 (0.780, 0.978) | – | 0.862 | 0.750 | |
AUC, area under the curve; CI, confidence interval; EMVI, extramural venous invasion; FD, fractal dimension; FD-IC, FD from iodine concentration images; FD-Zeff, FD from effective atomic number images; PNI, peripheral nerve invasion; TD, tumor deposit.
Table 4
| Parameter | PNI | EMVI/TD | |||||
|---|---|---|---|---|---|---|---|
| AUC difference | Z | P | AUC difference | Z | P | ||
| Development cohort (n=85) | |||||||
| FD-Con vs. FD-IC | 0.072 | 1.064 | 0.287 | 0.188 | 2.664 | 0.008 | |
| FD-IC vs. FD-Zeff | 0.004 | 0.060 | 0.952 | 0.106 | 1.862 | 0.063 | |
| FD-Zeff vs. FD-Con | 0.076 | 1.086 | 0.277 | 0.294 | 4.044 | <0.001 | |
| Validation cohort (n=45) | |||||||
| FD-Con vs. FD-IC | 0.189 | 1.716 | 0.086 | 0.237 | 2.173 | 0.030 | |
| FD-IC vs. FD-Zeff | 0.020 | 0.245 | 0.806 | 0.059 | 0.850 | 0.395 | |
| FD-Zeff vs. FD-Con | 0.208 | 2.016 | 0.044 | 0.296 | 3.035 | 0.002 | |
AUC, area under the curve; EMVI, extramural venous invasion; FD, fractal dimension; FD-Con, FD from conventional mixed-energy images; FD-IC, FD from iodine concentration images; FD-Zeff, FD from effective atomic number images; PNI, peripheral nerve invasion; TD, tumor deposit.
Univariate and multivariate analyses
Table 5 presents the results of the univariate and multivariate logistic regression analyses for the development cohort. FD-IC was identified as an independent risk factor for predicting PNI [odds ratio (OR) 77.873; 95% CI: 1.362–5.447×103; P=0.039], while FD-Zeff was as an independent risk factor for EMVI/TD status (OR 1.109×104, 95% CI: 25.550–1.907×107; P=0.006).
Table 5
| Parameter | Univariate analysis | Multivariate analysis | |||||
|---|---|---|---|---|---|---|---|
| Odds ratio | 95% CI | P | Odds ratio | 95% CI | P | ||
| PNI | |||||||
| Age | 1.025 | 0.985, 1.070 | 0.227 | ||||
| Gender | 0.754 | 0.316, 1.788 | 0.521 | ||||
| Clinical T stage | 1.668 | 0.316, 1.788 | 0.259 | ||||
| Clinical N stage | 1.390 | 0.589, 3.316 | 0.453 | ||||
| Tumor length | 0.964 | 0.763, 1.213 | 0.754 | ||||
| Tumor thickness | 1.529 | 0.600, 4.053 | 0.377 | ||||
| CEA | 2.058 | 0.785, 5.546 | 0.143 | ||||
| CA19-9 | 2.788 | 0.513, 2.096 | 0.252 | ||||
| FD-Con | 44.410 | 0.817, 3.321×103 | 0.071 | 2.253 | 0.021, 2.473×102 | 0.730 | |
| FD-IC | 3.803×102 | 11.010, 2.019×104 | 0.002 | 73.870 | 1.362, 5.447×103 | 0.039 | |
| FD-Zeff | 5.016×102 | 9.357, 4.454×104 | 0.004 | 44.770 | 0.523, 5.963×103 | 0.100 | |
| EMVI/TD | |||||||
| Age | 1.025 | 0.978, 1.082 | 0.324 | ||||
| Gender | 0.871 | 0.311, 2.461 | 0.791 | ||||
| Clinical T stage | 0.420 | 0.146, 1.188 | 0.102 | ||||
| Clinical N stage | 1.181 | 0.423, 3.335 | 0.750 | ||||
| Tumor length | 0.766 | 0.559, 1.018 | 0.077 | 0.756 | 0.508, 1.077 | 0.138 | |
| Tumor thickness | 2.083 | 0.692, 6.507 | 0.193 | ||||
| CEA | 3.343 | 1.132, 9.959 | 0.028 | 2.300 | 0.635, 8.361 | 0.199 | |
| CA19-9 | 3.938 | 0.675, 23.100 | 0.112 | ||||
| FD-Con | 4.125 | 0.040, 4.891×102 | 0.551 | ||||
| FD-IC | 5.869×102 | 9.111, 6.998×104 | 0.005 | 7.346 | 0.026, 2.355×103 | 0.489 | |
| FD-Zeff | 7.811×104 | 3.732×102, 4.939×107 | <0.001 | 1.109×104 | 25.550, 1.907×107 | 0.006 | |
CA19-9, carbohydrate antigen 19-9; CEA carcinoembryonic antigen; CI, confidence interval; EMVI, extramural venous invasion; FD, fractal dimension; FD-Con, FD from conventional mixed-energy images; FD-IC, FD from iodine concentration images; FD-Zeff, FD from effective atomic number images; N, node; PNI, peripheral nerve invasion; T, tumor; TD, tumor deposit.
Discussion
The preoperative diagnosis of PNI and EMVI/TD is particularly valuable for personalized treatment planning and improving the prognosis of patients with RC (23,24). This study primarily assessed the diagnostic efficacy of FD parameters derived from conventional and functional DECT imaging of RC. FD-IC and FD-Zeff were elevated in positive patients compared to negative patients in both the PNI and EMVI/TD groups. Multivariate logistic regression analysis indicated that FD-IC was an independent predictor of PNI, while FD-Zeff served as an independent predictor of EMVI/TD. Additionally, we employed these factors to develop integrated models for PNI and EMVI/TD. In the validation cohort, the AUC, sensitivity, and specificity of the PNI combination model were 0.771, 0.682, and 0.826, respectively, whereas the corresponding parameters for the EMVI/TD combination model were 0.879, 0.862, and 0.750, respectively. The diagnostic accuracy of PNI and EMVI/TD was significantly enhanced with the combination of multiparameter risk variables.
Disordered angiogenesis is a fundamental characteristic of malignant tumors, being directly associated with the proliferation and metabolism of cancer cells. The irregular blood supply distribution within the tumor microenvironment, along with the influence of various factors, alters the invasiveness and self-renewal capacity of cancer cells at the tumor periphery. This leads to disordered vascular growth patterns and induces hypoxia and necrosis within the tumor, which exhibit distinct fractal characteristics (25). Tumor heterogeneity is associated with aggressiveness and affects treatment choices and prognosis (26). Fractal analysis has been shown to be capable of efficiently evaluating tumor heterogeneity (27). Previous research on fractal analysis has predominantly used conventional CT images (28); however, DECT offers enhanced information relative to traditional CT, including perfusion parameters (29), which might directly indicate the angiogenesis and blood supply of tumor tissue (30). This study used a multidimensional approach to examine fractal analysis with iodine and Zeff maps to predict EMVI/TD and PNI status in patients with RC. Our findings indicate that FD-IC and FD-Zeff may effectively evaluate PNI and EMVI/TD status, with both the PNI- and EMVI/TD-positive groups exhibiting higher FD than their negative counterparts. These results aligned with a previous study that used fractal analysis to predict advanced T stage and pathological stage in gastric cancer, with FD serving as an independent predictive factor in the lymph node-positive subgroup of gastric cancer (31). Highly invasive tumor tissues exhibit elevated microvessel density, which augments the proliferative capacity and disarray of cancer cells (21). Meanwhile, elevated tumor heterogeneity is induced by factors such as hypoxia and necrosis, ultimately resulting in a rise in FD (18). The data of our study indicated that a high FD value correlate with an increased likelihood for tumor invasion.
We subsequently compared the diagnostic efficacy of FD-Con with that of FD-IC and FD-Zeff parameters in the study groups. Within the PNI subgroup, both FD-IC and FD-Zeff yielded a higher AUC than did FD-Con; however, significant differences were observed between FD-Zeff and FD-Con only in the validation cohort (P=0.044). In the EMVI/TD subgroup, we observed that FD-IC and FD-Zeff exhibited considerably superior diagnostic accuracy for the evaluation of EMVI/TD status compared to FD-Con in both the development and validation cohorts (all P values <0.01). Given the significant interobserver variability of FD-Con, we hypothesize that FD-IC and FD-Zeff offer more benefits in evaluating these high-risk factors.
In one study, an iodine map–based radiomics model exhibited superior diagnostic efficacy in predicting peritoneal metastasis in patients with gastric cancer compared to the standard CT-based radiomics model, which is in line with our findings (32). In contrast to gross pathology and multiparametric MRI, standard CT may inadequately represent the tumor size, hence diminishing the precision of the real FD measurement of the tumor (33). DECT enhances tumor identification accuracy by reflecting the internal perfusion heterogeneity of tumor tissue, thereby improving the diagnostic utility of fractal analysis in the assessment of PNI and EMVI/TD status in patients with RC as compared to conventional CT (11,34,35).
At present, a dependable noninvasive method for predicting PNI and EMVI/TD status in patients with RC is absent for clinical settings, and fractal analysis with DECT may represent a viable option. From a clinical standpoint, compared with radiomics and deep learning approaches, fractal analysis demonstrates superior interpretability and reduced computational demands in assessing tumor aggressiveness, thereby offering a novel strategy for the pretherapeutic evaluation of RC (36,37).
This study involved several limitations that should be addressed. First, owing to the modest cohort size and the low counts of EMVI-positive and TD-positive patients, EMVI and TD were analyzed in a combined subgroup. Second, patients with distant metastases were excluded to avoid confounding effects on tumor vascular invasion patterns. Third, all data were acquired with a fixed CT scanner and imaging protocol at a single center, and external validation and protocol refinements will be implemented in subsequent studies. Finally, regions of interest were delineated according to the maximum axial area image. However, malignancies can be more precisely recognized in a sequence of consecutive cross-axis images in 3D analysis than they can be in a single cross-axis image in 2D analysis (38).
Conclusions
DECT-based fractal analysis can effectively assess PNI and EMVI/TD status in patients with RC, complementing existing clinical guidelines.
Acknowledgments
None.
Footnote
Reporting Checklist: The authors have completed the STARD reporting checklist. Available at https://qims.amegroups.com/article/view/10.21037/qims-2025-862/rc
Data Sharing Statement: Available at https://qims.amegroups.com/article/view/10.21037/qims-2025-862/dss
Funding: None.
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://qims.amegroups.com/article/view/10.21037/qims-2025-862/coif). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. This study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments. Written informed consent was obtained from all participants. The study protocol was approved by the Ethics Committees of the Second Affiliated Hospital of Chongqing Medical University (No. 2020-366).
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Benson AB, Venook AP, Al-Hawary MM, Azad N, Chen YJ, Ciombor KK, et al. Rectal Cancer, Version 2.2022, NCCN Clinical Practice Guidelines in Oncology. J Natl Compr Canc Netw 2022;20:1139-67. [Crossref] [PubMed]
- Siegel RL, Giaquinto AN, Jemal A. Cancer statistics, 2024. CA Cancer J Clin 2024;74:12-49. [Crossref] [PubMed]
- Keller DS, Berho M, Perez RO, Wexner SD, Chand M. The multidisciplinary management of rectal cancer. Nat Rev Gastroenterol Hepatol 2020;17:414-29. [Crossref] [PubMed]
- Riffet M, Dupont B, Faisant M, Cerasuolo D, Menahem B, Alves A, Dubois F, Levallet G, Bazille C. New Histoprognostic Factors to Consider for the Staging of Colon Cancers: Tumor Deposits, Invasive Tumor Infiltration and High-Grade Budding. Int J Mol Sci 2023;24:3573. [Crossref] [PubMed]
- Sohn B, Lim JS, Kim H, Myoung S, Choi J, Kim NK, Kim MJ. MRI-detected extramural vascular invasion is an independent prognostic factor for synchronous metastasis in patients with rectal cancer. Eur Radiol 2015;25:1347-55. [Crossref] [PubMed]
- Smith NJ, Barbachano Y, Norman AR, Swift RI, Abulafi AM, Brown G. Prognostic significance of magnetic resonance imaging-detected extramural vascular invasion in rectal cancer. Br J Surg 2008;95:229-36. [Crossref] [PubMed]
- Lord AC, D'Souza N, Pucher PH, Moran BJ, Abulafi AM, Wotherspoon A, Rasheed S, Brown G. Significance of extranodal tumour deposits in colorectal cancer: A systematic review and meta-analysis. Eur J Cancer 2017;82:92-102. [Crossref] [PubMed]
- Lord AC, D'Souza N, Shaw A, Rokan Z, Moran B, Abulafi M, Rasheed S, Chandramohan A, Corr A, Chau I, Brown G. MRI-Diagnosed Tumor Deposits and EMVI Status Have Superior Prognostic Accuracy to Current Clinical TNM Staging in Rectal Cancer. Ann Surg 2022;276:334-44. [Crossref] [PubMed]
- van Wyk HC, Going J, Horgan P, McMillan DC. The role of perineural invasion in predicting survival in patients with primary operable colorectal cancer: A systematic review. Crit Rev Oncol Hematol 2017;112:11-20. [Crossref] [PubMed]
- Knijn N, Mogk SC, Teerenstra S, Simmer F, Nagtegaal ID. Perineural Invasion is a Strong Prognostic Factor in Colorectal Cancer: A Systematic Review. Am J Surg Pathol 2016;40:103-12. [Crossref] [PubMed]
- Arico' FM, Trimarchi R, Portaluri A, Barilla' C, Migliaccio N, Bucolo GM, Cicero G, Sofia C, Booz C, Vogl TJ, Marino MA, Ascenti V, D'Angelo T, Mazziotti S, Ascenti G. Virtual monoenergetic dual-layer dual-energy CT images in colorectal cancer: CT diagnosis could be improved? Radiol Med 2023;128:891-9. [Crossref] [PubMed]
- Zhang J, Qi H, Yang C, Liu L, Wang Y, Li W. Preoperative prediction of lymphovascular invasion and T-staging of rectal cancer via a dual-energy computed tomography iodine map: a feasibility study. Diagn Interv Radiol 2025;31:1-9. [Crossref] [PubMed]
- Gao W, Zhang Y, Dou Y, Zhao L, Wu H, Yang Z, Liu A, Zhu L, Hao F. Association between extramural vascular invasion and iodine quantification using dual-energy computed tomography of rectal cancer: a preliminary study. Eur J Radiol 2023;158:110618. [Crossref] [PubMed]
- Wang D, Zhuang Z, Wu S, Chen J, Fan X, Liu M, Zhu H, Wang M, Zou J, Zhou Q, Zhou P, Xue J, Meng X, Ju S, Zhang L. A Dual-Energy CT Radiomics of the Regional Largest Short-Axis Lymph Node Can Improve the Prediction of Lymph Node Metastasis in Patients With Rectal Cancer. Front Oncol 2022;12:846840. [Crossref] [PubMed]
- Park YW, Kim S, Ahn SS, Han K, Kang SG, Chang JH, Kim SH, Lee SK, Park SH. Magnetic resonance imaging-based 3-dimensional fractal dimension and lacunarity analyses may predict the meningioma grade. Eur Radiol 2020;30:4615-22. [Crossref] [PubMed]
- Castello A, Russo C, Grizzi F, Qehajaj D, Lopci E. Prognostic Impact of Intratumoral Heterogeneity Based on Fractal Geometry Analysis in Operated NSCLC Patients. Mol Imaging Biol 2019;21:965-72. [Crossref] [PubMed]
- Lennon FE, Cianci GC, Cipriani NA, Hensing TA, Zhang HJ, Chen CT, Murgu SD, Vokes EE, Vannier MW, Salgia R. Lung cancer-a fractal viewpoint. Nat Rev Clin Oncol 2015;12:664-75. [Crossref] [PubMed]
- Tochigi T, Kamran SC, Parakh A, Noda Y, Ganeshan B, Blaszkowsky LS, Ryan DP, Allen JN, Berger DL, Wo JY, Hong TS, Kambadakone A. Response prediction of neoadjuvant chemoradiation therapy in locally advanced rectal cancer using CT-based fractal dimension analysis. Eur Radiol 2022;32:2426-36. [Crossref] [PubMed]
- Loughrey MB, Webster F, Arends MJ, Brown I, Burgart LJ, Cunningham C, Flejou JF, Kakar S, Kirsch R, Kojima M, Lugli A, Rosty C, Sheahan K, West NP, Wilson RH, Nagtegaal ID. Dataset for Pathology Reporting of Colorectal Cancer: Recommendations From the International Collaboration on Cancer Reporting (ICCR). Ann Surg 2022;275:e549-61. [Crossref] [PubMed]
- Albrecht MH, Vogl TJ, Martin SS, Nance JW, Duguay TM, Wichmann JL, De Cecco CN, Varga-Szemes A, van Assen M, Tesche C, Schoepf UJ. Review of Clinical Applications for Virtual Monoenergetic Dual-Energy CT. Radiology 2019;293:260-71. [Crossref] [PubMed]
- Chen JF, Yang J, Chen WJ, Wei X, Yu XL, Huang DD, Deng H, Luo YD, Liu XJ. Mono+ algorithm assessment of the diagnostic value of dual-energy CT for high-risk factors for colorectal cancer: a preliminary study. Quant Imaging Med Surg 2024;14:432-46. [Crossref] [PubMed]
- Captur G, Muthurangu V, Cook C, Flett AS, Wilson R, Barison A, Sado DM, Anderson S, McKenna WJ, Mohun TJ, Elliott PM, Moon JC. Quantification of left ventricular trabeculae using fractal analysis. J Cardiovasc Magn Reson 2013;15:36. [Crossref] [PubMed]
- Zhang XY, Wang S, Li XT, Wang YP, Shi YJ, Wang L, Wu AW, Sun YS. MRI of Extramural Venous Invasion in Locally Advanced Rectal Cancer: Relationship to Tumor Recurrence and Overall Survival. Radiology 2018;289:677-85. [Crossref] [PubMed]
- Uccello TP, Lesch ML, Ullman NA, Kintzel SA, Gradzewicz LB, Velagaleti T, Fleming FJ, Mills BN, Murphy JD, Garrett-Larsen J, Qiu H, Drage MG, Ye J, Gavras NW, Johnston CJ, Love TMT, Repasky EA, Linehan DC, Lord EM, Gerber SA. Radiation Therapy Exacerbates Tumor-Promoting Innervation and Nerve Signaling in Rectal Cancer. Int J Radiat Oncol Biol Phys 2023;115:733-45. [Crossref] [PubMed]
- Tauriello DV, Calon A, Lonardo E, Batlle E. Determinants of metastatic competency in colorectal cancer. Mol Oncol 2017;11:97-119. [Crossref] [PubMed]
- Cusumano D, Dinapoli N, Boldrini L, Chiloiro G, Gatta R, Masciocchi C, Lenkowicz J, Casà C, Damiani A, Azario L, Van Soest J, Dekker A, Lambin P, De Spirito M, Valentini V. Fractal-based radiomic approach to predict complete pathological response after chemo-radiotherapy in rectal cancer. Radiol Med 2018;123:286-95. [Crossref] [PubMed]
- Hayano K, Lee SH, Yoshida H, Zhu AX, Sahani DV. Fractal analysis of CT perfusion images for evaluation of antiangiogenic treatment and survival in hepatocellular carcinoma. Acad Radiol 2014;21:654-60. [Crossref] [PubMed]
- Hayano K, Yoshida H, Zhu AX, Sahani DV. Fractal analysis of contrast-enhanced CT images to predict survival of patients with hepatocellular carcinoma treated with sunitinib. Dig Dis Sci 2014;59:1996-2003. [Crossref] [PubMed]
- Zhang LJ, Wu S, Wang M, Lu L, Chen B, Jin L, Wang J, Larson AC, Lu GM. Quantitative dual energy CT measurements in rabbit VX2 liver tumors: Comparison to perfusion CT measurements and histopathological findings. Eur J Radiol 2012;81:1766-75. [Crossref] [PubMed]
- Moding EJ, Clark DP, Qi Y, Li Y, Ma Y, Ghaghada K, Johnson GA, Kirsch DG, Badea CT. Dual-energy micro-computed tomography imaging of radiation-induced vascular changes in primary mouse sarcomas. Int J Radiat Oncol Biol Phys 2013;85:1353-9. [Crossref] [PubMed]
- Watanabe H, Hayano K, Ohira G, Imanishi S, Hanaoka T, Hirata A, Kano M, Matsubara H. Quantification of Structural Heterogeneity Using Fractal Analysis of Contrast-Enhanced CT Image to Predict Survival in Gastric Cancer Patients. Dig Dis Sci 2021;66:2069-74. [Crossref] [PubMed]
- Chen Y, Xi W, Yao W, Wang L, Xu Z, Wels M, Yuan F, Yan C, Zhang H. Dual-Energy Computed Tomography-Based Radiomics to Predict Peritoneal Metastasis in Gastric Cancer. Front Oncol 2021;11:659981. [Crossref] [PubMed]
- Michallek F, Haouari MA, Dana O, Perrot A, Silvera S, Dallongeville A, Dewey M, Zins M. Fractal analysis improves tumour size measurement on computed tomography in pancreatic ductal adenocarcinoma: comparison with gross pathology and multi-parametric MRI. Eur Radiol 2022;32:5053-63. [Crossref] [PubMed]
- Wu J, Zhang Q, Zhao Y, Liu Y, Chen A, Li X, Wu T, Li J, Guo Y, Liu A. Radiomics Analysis of Iodine-Based Material Decomposition Images With Dual-Energy Computed Tomography Imaging for Preoperatively Predicting Microsatellite Instability Status in Colorectal Cancer. Front Oncol 2019;9:1250. [Crossref] [PubMed]
- Wang G, Liu X, Zhou J. Differentiating gastric schwannoma from gastric stromal tumor (≤5 cm) by histogram analysis based on iodine-based material decomposition images: a preliminary study. Front Oncol 2023;13:1243300. [Crossref] [PubMed]
- Michallek F, Dewey M. Fractal analysis in radiological and nuclear medicine perfusion imaging: a systematic review. Eur Radiol 2014;24:60-9. [Crossref] [PubMed]
- Stanzione A, Verde F, Romeo V, Boccadifuoco F, Mainenti PP, Maurea S. Radiomics and machine learning applications in rectal cancer: Current update and future perspectives. World J Gastroenterol 2021;27:5306-21. [Crossref] [PubMed]
- Sanghera B, Banerjee D, Khan A, Simcock I, Stirling JJ, Glynne-Jones R, Goh V. Reproducibility of 2D and 3D fractal analysis techniques for the assessment of spatial heterogeneity of regional blood flow in rectal cancer. Radiology 2012;263:865-73. [Crossref] [PubMed]

