ABSTRACT
PURPOSE
This study aimed to evaluate the prognostic significance of tumoral apparent diffusion coefficient (tADC), peritumoral ADC (pADC), and the peritumoral-to-tumoral ADC (p/tADC) ratio in invasive breast cancers using high-resolution diffusion-weighted imaging.
METHODS
A retrospective cohort of 149 patients with invasive breast cancer was analyzed. pADC, tADC, and p/tADC values were independently measured by two experienced radiologists. The associations between these parameters and histological grade, receptor status, Ki-67 index, and molecular subtypes were assessed. Interobserver agreement was evaluated using the intraclass correlation coefficient (ICC), and multivariate logistic regression was performed to identify independent predictors.
RESULTS
Although tADC was significantly associated only with tumor size (P < 0.050), both pADC and p/tADC were significantly associated with high histological grade, estrogen receptor/progesterone receptor negativity, human epidermal growth factor receptor 2 positivity, and high Ki-67 index(P < 0.01). The p/tADC ratio emerged as an independent predictor of high-grade tumors and Ki-67 ≥ 20%. Excellent interobserver agreement was observed for all ADC measurements (ICC > 0.90).
CONCLUSION
The p/tADC ratio reflects stromal and microenvironmental alterations associated with tumor aggressiveness and demonstrates stronger prognostic relevance than tADC alone. Incorporating the p/tADC ratio into routine magnetic resonance imaging (MRI) interpretation may improve preoperative risk stratification and support personalized treatment planning in invasive breast cancer.
CLINICAL SIGNIFICANCE
The p/tADC ratio reflects stromal and microenvironmental changes associated with tumor aggressiveness and may serve as a noninvasive imaging biomarker. Its integration into preoperative MRI protocols may facilitate risk stratification and personalized treatment planning.
Main points
• The peritumoral-to-tumoral apparent diffusion coefficient (p/tADC) ratio is significantly elevated in tumors with aggressive biological features.
• Unlike tumoral ADC, peritumoral ADC and p/tADC show robust correlations with multiple prognostic markers.
• ADC measurements demonstrated excellent interobserver reproducibility, supporting their clinical reliability.
Breast cancer is among the most frequently encountered malignancies worldwide, with an annual incidence of approximately 2.3 million cases and over 666,000 deaths, ranking as the fourth leading cause of cancer-related mortality.1 To address this burden, global efforts have intensified through early detection strategies, including national screening programs and awareness campaigns.
Tumor microenvironment
Breast cancers exhibit remarkable heterogeneity in both genetic and cellular behavior, influencing disease progression and treatment response.2, 3 The tumor microenvironment (TME)—comprising immune cells, fibroblasts, and the extracellular matrix—plays a pivotal role in tumor development.4-9 In invasive breast cancers, phenomena such as lymphatic obstruction and matrix remodeling around tumors contribute to peritumoral edema, a known feature of malignancy.5 This edema, linked to increased vascular permeability via proteolytic enzymes and angiogenesis, has been associated with more aggressive tumor behavior and poorer prognosis.10-13 Notably, the peritumoral stroma harbors genetic distinctions from intratumoral regions, suggesting it provides biologically distinct and prognostically relevant information.14-16
Magnetic resonance imaging and diffusion-weighted imaging
Magnetic resonance imaging (MRI) offers superior sensitivity for breast cancer detection, enabling detailed anatomical and functional assessments.17 However, its specificity remains limited.18 Diffusion-weighted imaging (DWI), which measures water-molecule movement, enhances diagnostic accuracy by providing apparent diffusion coefficient (ADC) values, which are typically lower in malignant tissues due to higher cellularity.19-22
Several studies have explored correlations between ADC values and histopathologic features, including tumor size, grade, lymphovascular invasion, nodal status, and molecular subtype; however, results have been inconsistent due to methodological variability.23-25 Peritumoral edema has emerged as a promising imaging marker, but limited data exist on its association with histopathologic aggressiveness or DWI parameters.26, 27
In this study, the term peritumoral region specifically refers to the stromal tissue surrounding the tumor rather than to peritumoral edema alone. Although peritumoral edema seen on T2-weighted imaging may be a visual manifestation of stromal changes, our diffusion-weighted analysis focused on quantitative ADC measurements from the biologically active peritumoral stroma, excluding cystic, necrotic, or purely edematous regions.
Aim
This study investigates the relationship between MRI features and histopathological characteristics of breast cancer, with a specific focus on the prognostic significance of the peritumoral-to-tumoral ADC (p/tADC) ratio.
Methods
Study group
This retrospective study was conducted with the approval of the Ethics Committee of Ankara Bilkent City Hospital (decision number: E1-23-3199, date: 11.01.2023). Between April 2020 and July 2023, 367 patients who underwent high-resolution DWI with a narrow-field-of-view protocol on a 3-Tesla (T) MRI scanner were reviewed. Among these, 149 women (mean age: 53 ± 10.5 years; range 33–86 years) who had preoperative breast MRI- and biopsy-confirmed invasive breast cancer were included in the study. Patients with pure in situ carcinoma, a prior history of neoadjuvant chemotherapy, or the absence of surgery following breast MRI were excluded from the study.
Demographic data and the pathology reports of the patients were retrieved from the hospital information management system.
The patient selection process is summarized in Figure 1.
Magnetic resonance imaging protocol
All breast MRI examinations were performed on a 3.0-T scanner (Signa Pioneer, General Electric Medical Systems, Milwaukee, WI, USA) in the prone position using a dedicated 16-channel phased-array bilateral breast coil (NeoCoil 3.0 T, General Electric Medical Systems, Milwaukee, WI, USA).
The imaging protocol included the following sequences:
• Axial T1-weighted fast spin-echo [repetition time (TR)/echo time (TE): 484/8.5 ms; slice thickness: 5 mm]
• Axial and sagittal T2-weighted fat-suppressed fast spin-echo (TR/TE: 4,786/82.5 ms and 5,116/83.5 ms, respectively; slice thickness:
5 mm)
• Dynamic contrast-enhanced (DCE) axial T1-weighted fat-suppressed sequence (TR/TE: 6.1/1.7 ms; slice thickness: 2.2 mm), acquired before and after intravenous administration of a gadolinium-based contrast agent (0.1 mmol/kg), followed by a 20-mL saline flush
• DWI using single-shot echo-planar imaging with b values of 50 and 800 s/mm². ADC maps were automatically generated on the workstation.
Histopathological evaluation
Histopathological specimens obtained from surgical samples were evaluated by the Pathology Department of Ankara Bilkent City Hospital. Pathological data included information on estrogen receptor (ER) status, progesterone receptor (PR) status, human epidermal growth factor receptor 2 (HER2) status, Ki-67 index, and tumor size.
Apparent diffusion coefficient measurements
ADC measurements and T1-weighted, T2-weighted, and DCE images were used for anatomical reference. On the ADC map, the slice with the largest lesion cross-sectional diameter was selected. For tADC measurement, a circular region of interest (ROI) with a fixed diameter of 10 pixels (approximately 10 mm²) was manually placed within the most restricted (visually darkest) solid portion of the tumor, avoiding necrotic, cystic, hemorrhagic, and fatty areas. The lowest ADC value obtained from three measurements was recorded as the representative tADC. For peritumoral ADC (pADC) assessment, three ROIs of identical size and shape (10 pixels in diameter) were placed in the adjacent breast parenchyma within a 5–10 mm distance from the tumor border, excluding fibroglandular tissue, large vessels, ducts, and artifacts. These ROIs were placed in areas with the highest signal intensity on the ADC map, consistent with increased diffusivity. Among these, the highest ADC value was recorded as the pADC, in line with previously published methods.28 The p/tADC ratio was calculated by dividing the pADC by the tADC. Sample images of tumors with high and low pADC values were shown in Figures 2 and 3, respectively. All measurements were performed by two radiologists with over 5 years of experience in breast imaging. Each measurement was performed independently, and the mean value was used for analysis. The interobserver agreement was assessed using the intraclass correlation coefficient (ICC).
Pathological assessment
All patients underwent pathological evaluation of surgical specimens. Parameters such as histological grade, ER status, PR status, HER2 status, Ki-67 index, and axillary lymph node status were analyzed as binary variables as follows:
• ER/PR status was evaluated using the Allred scoring system, and a total score >2 was considered positive.
• The HER2 score was evaluated immunohistochemically on a scale of 0–3; scores of 0–1 were considered negative, and a score of 3 was considered positive. Cases with a score of 2 underwent further evaluation with fluorescent in situ hybridization testing.
• Ki-67 values < 20% were classified as negative and values ≥ 20% as positive.
Statistical analysis
All statistical analyses were performed using SPSS version 30.0 (IBM, Chicago, IL, USA). ADC values were considered continuous dependent variables. Histopathological parameters were treated as binary independent variables, including the following:
• Histological grade: Grade 1–2 vs. Grade 3
• ER and PR status: Negative vs. Positive
• HER2 status: Negative vs. Positive
• Ki-67: < 20% vs. ≥ 20%
• Axillary lymph node metastasis: Negative vs. Positive
• Tumor type: Ductal vs. Lobular
• Molecular subtype: Luminal (luminal A, luminal B and luminal B-HER2+) vs. non luminal (HER2-enriched, basal-like)
The normality of data distribution was assessed using the Shapiro–Wilk test. For normally distributed variables, comparisons between two groups were conducted using the independent samples Student’s t-test, and results were expressed as mean ± standard deviation (SD).
For comparisons among luminal molecular subtypes (luminal A, luminal B, luminal B-HER2+, HER2-enriched, and basal-like), one-way analysis of variance was performed. When overall significance was observed, post-hoc pairwise comparisons were conducted using the Tukey test to identify intergroup differences.
Two-tailed tests were used for all analyses, and a P value < 0.05 was considered statistically significant.
Receiver operating characteristic (ROC) curve analysis was performed to determine optimal cut-off values for the p/tADC ratio to predict aggressive tumor characteristics. For each cut-off, sensitivity, specificity, and positive likelihood ratios (LR+) were calculated to assess discriminative performance.
Variables with P < 0.05 in univariate analysis were entered into a multivariate binary logistic regression model to identify independent predictors of a high p/tADC ratio. The p/tADC ratio was dichotomized at 2.3 based on ROC analysis; ROC curves were generated to evaluate the ability of the p/tADC ratio to distinguish high-grade tumors (grade 3) and a high Ki-67 index (≥ 20%) from lower-risk tumors. The cut-off value was selected by maximizing the Youden index across these prognostic measures. Model calibration and performance were assessed using the Hosmer–Lemeshow goodness-of-fit test and Nagelkerke R².
Interobserver agreement
To assess measurement reproducibility between the two radiologists, interobserver agreement was evaluated for tADC, pADC, and p/tADC values in a randomly selected subset of 40 lesions (approximately 25% of the sample). The ICC (two-way random-effects model, absolute agreement) was calculated. Agreement was interpreted as follows: poor (< 0.50), moderate (0.50–0.75), good (0.75–0.90), and excellent (> 0.90).
Additionally, Bland–Altman plots were generated to visualize systematic bias and 95% limits of agreement between the two observers for each parameter.
Results
Between April 2020 and July 2023, a total of 149 women (mean age: 53 ± 10.5 years; range: 33–86 years) with invasive breast cancer were included in the study.
Interobserver agreement
Excellent interobserver agreement was observed for all ADC-based measurements. The ICC values for tADC, pADC, and the p/tADC ratio were 0.94 [95% confidence interval (CI): 0.90–0.97], 0.91 (95% CI: 0.86–0.95), and 0.93 (95% CI: 0.88–0.96), respectively (all P < 0.001).
Bland–Altman analysis demonstrated a mean difference of 0.02 × 10-³ mm²/s for tADC and 0.05 × 10-³ mm²/s for pADC, with no significant proportional bias. All differences lay within ± 1.96 SD, confirming acceptable agreement between readers (Table 1).
Demographic data and tumor characteristics are summarized in Table 2.
tADC values were significantly associated with tumor size (P = 0.005). However, no significant associations were found with histological grade (P = 0.252), Ki-67 index (P = 0.635), ER (P = 0.562), PR (P = 0.652), HER2 (P = 0.556), axillary lymph node status (P = 0.957), or molecular subtype (P = 0.703).
pADC was significantly associated with tumor size (P = 0.009), axillary lymph node positivity (P = 0.02), molecular subtype (P < 0.001), HER2 (P = 0.002), ER (P = 0.001), PR (P = 0.002), Ki-67 (P < 0.001), and histological grade (P = 0.010).
The p/tADC ratio showed significant associations with tumor size (P < 0.001), histological grade (P = 0.004), Ki-67 index (P < 0.001), ER (P = 0.006), PR (P = 0.015), and HER2 (P = 0.010) but not with lymph node positivity (P = 0.08).
Associations between pADC and tADC values, p/tADC ratios, histopathological parameters, and lymph nodes status are shown in Table 3.
Comparison by molecular subtype
No statistically significant difference in tADC values was observed between luminal and non-luminal subgroups (P = 0.703). By contrast, both pADC and p/tADC values demonstrated a statistically significant difference among the luminal subtypes (P < 0.001 and P = 0.002, respectively).
No statistically significant differences in pADC (P = 0.085) or p/tADC ratios (P = 0.327) were observed between luminal A and luminal B groups. However, the luminal A group demonstrated significantly lower pADC and p/tADC ratios than the luminal B-HER2+ (P = 0.003, P = 0.007, respectively), HER2-enriched (P < 0.001, P = 0.035, respectively), and basal-like tumor groups (P = 0.001, P = 0.003, respectively).
By contrast, no significant differences were found among the luminal B, luminal B-HER2+, HER2-enriched, and basal-like groups in terms of pADC and p/tADC ratios. Additionally, tADC values did not differ significantly across molecular subtypes.
Prognostic performance of peritumoral-to-tumoral apparent diffusion coefficient ratios
Analysis of the prognostic predictive power of the p/tADC ratio revealed the highest discriminative performance for tumor size (LR+: 1.99), whereas the lowest was observed for HER2 status (LR+: 1.56). Cut-off values, sensitivity, specificity, and LR+ are summarized in Table 4.
Multivariate logistic regression
In the multivariate logistic regression analysis, the model demonstrated acceptable fit (Hosmer–Lemeshow χ²: 11.55, df: 8, P = 0.17) and accounted for approximately 24% of the variance in the high p/tADC ratio (Nagelkerke R²: 0.24). High histological grade and a Ki-67 index ≥ 20% were identified as independent predictors of a high p/tADC ratio (≥ 2.3), whereas ER, PR, and HER2 status did not remain significant in the multivariate model (Table 5). The ROC analysis confirmed that a p/tADC threshold of 2.3 provided the optimal balance of sensitivity and specificity for differentiating high-grade tumors and elevated Ki-67 expression.
Discussion
This study underscores the prognostic potential of the p/tADC ratio as a non-invasive imaging biomarker in invasive breast cancer. We demonstrated that although tADC was significantly associated only with tumor size, the p/tADC ratio showed independent associations with high histological grade and elevated Ki-67 index and was significantly different across molecular subtypes. These findings suggest that alterations in diffusion within the peritumoral microenvironment provide clinically relevant biological information beyond what tADC alone can reveal.
The inverse association between ADC values and cellularity is well established, as lower tADC values often indicate high cellular density and reduced extracellular space.23, 29 However, prior studies have reported inconsistent correlations between ADC and key pathological markers such as receptor status, Ki-67, and tumor grade.20, 30 Our findings are consistent with this inconsistency, as tADC did not correlate significantly with ER, PR, HER2, or Ki-67 in our cohort.
Interobserver agreement for all diffusion metrics was excellent (ICC > 0.90), confirming the reproducibility of our measurement protocol.
In contrast to tADC, both pADC and the p/tADC ratio showed significant associations with histological grade, hormone receptor status, HER2 expression, and Ki-67 index. This supports the notion that the peritumoral stroma is an active component of the TME and may mirror stromal remodeling, angiogenesis, immune suppression, and extracellular matrix changes.6, 16, 31, 32
Several histopathological and imaging-based studies have identified edema, inflammation, and stromal activation in peritumoral tissues of high-grade tumors.11, 33 Zhang et al.33 showed that peritumoral edema is more pronounced in HER2+ and triple-negative tumors, reflecting elevated vascular permeability and cytokine-mediated stromal stress.34 These findings are further supported by DWI-based studies linking high pADC values to HER2 overexpression and high Ki-67,12, 35 which align with our observation that aggressive tumors exhibit higher pADC and p/tADC ratios.
From a molecular perspective, our analysis revealed that HER2-enriched, basal-like, and triple-negative tumors had significantly higher pADC and p/tADC values than luminal A subtypes. These subtypes are known to exhibit more intense stromal remodeling and angiogenesis.36-39
Notably, tADC did not differentiate tumor grade, whereas the p/tADC ratio did (P = 0.004 vs. P = 0.252). This supports the findings of Doğan et al.,40 who argued that peritumoral diffusion metrics better reflect biological aggressiveness, and those of Choi et al.,41 who suggested that the p/tADC ratio mitigates variability introduced by tumor heterogeneity and imaging artifacts
We also calculated cut-off values for predicting high-grade features. The optimal p/tADC threshold of 2.3, determined via ROC analysis, yielded the best balance between sensitivity and specificity. The p/tADC ratio was moderately predictive across different features, with the highest LR+ for tumor size (LR+: 1.99) and the lowest for HER2 status (LR+: 1.56). Although not highly diagnostic on their own, these values offer quantifiable markers that can assist in clinical risk stratification.
In multivariate analysis, only high histological grade and Ki-67 ≥ 20% retained statistical significance, emphasizing the p/tADC ratio’s closer relationship with proliferative activity and stromal remodeling rather than receptor expression per se.
Our findings are consistent with those of Okuma et al.,28 who also identified associations between p/tADC and tumor size and grade, Ki-67, and lymph node status. Similarly, El-Metwally et al.42 supported a link between ADC metrics and Ki-67, whereas associations with hormone receptors were weaker.
Incorporating peritumoral diffusion metrics into breast MRI has been supported by a recent meta-analysis published in 2024 (PMID: 38334760).43 However, many prior studies lacked interobserver reproducibility assessments or multivariate control. Our study addresses these gaps by demonstrating that the p/tADC ratio is not only reproducible (ICC > 0.90) but also an independent predictor of histopathological aggressiveness.
We also contributed to the literature by offering molecular subtype-specific analysis and defined p/tADC cut-offs for biologically aggressive features. Together, these findings position the p/tADC ratio as a promising adjunct to conventional MRI parameters for risk stratification in invasive breast cancer.
Our findings highlight the clinical value of the p/tADC ratio as a non-invasive imaging biomarker of tumor aggressiveness. The p/tADC ratio was independently associated with high histological grade and Ki-67 index and showed distinct patterns across molecular subtypes. By using a standardized ROI placement and defining specific cut-off values, this metric may aid in preoperative risk stratification. Integrating the p/tADC ratio into routine MRI protocols could support more personalized treatment planning. Future studies should also explore advanced diffusion techniques, such as diffusion tensor imaging, to better characterize peritumoral tissue heterogeneity.


