Front. Oncol.Frontiers in OncologyFront. Oncol.2234-943XFrontiers Media S.A.10.3389/fonc.2025.1539748OncologyOriginal ResearchPrediction of lymphovascular invasion in non-mass enhancement breast cancer using DCE-MRI and clinical-pathological featuresGuoShiqi1‡ZhaoKai1†‡XuHu1‡XieYujiao12LiQingyang1LiangJunqing3ChenSiyi1SunJiahong1GaoZhaofeng12ZhuLi1*WangJiandong1*1Department of General Surgery, The First Medical Center of Chinese People's Liberation Army General Hospital, Beijing, China2School of Medicine Nankai University, Tianjin, China3Department of Breast Surgery, Affiliated Cancer Hospital of Inner Mongolia Medical University, Hohhot, Inner Mongolia, China
Edited by: Sharon R. Pine, University of Colorado Anschutz Medical Campus, United States
Reviewed by: Bilgin Kadri Aribas, Bülent Ecevit University, Türkiye
Aya El Rawy, Suez Canal University, Egypt
*Correspondence: Li Zhu, zhuxiaoli0430@163.com; Jiandong Wang, vicky1968@163.com
†Present address: Kai Zhao, Beijing Institute for Brain Research, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing, China
‡These authors have contributed equally to this work and share first authorship
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Objective
The present study explores the relationship between the distribution patterns of non-mass enhancement (NME) type invasive breast cancer in dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) and lymphovascular invasion (LVI).
Methods
A retrospective analysis was conducted on 192 female patients with NME-type breast cancer who underwent DCE-MRI between January 2019 and December 2023. Based on postoperative pathological results, the patients were divided into two groups: LVI-positive [LVI(+)] (N = 50) and LVI-negative [LVI(-)] (N = 142). A generalized variance inflation factor (GVIF) analysis was used to identify variables with high multicollinearity. Multivariate logistic regression was used to analyze the risk factors associated with LVI. The performance of the Distribution of NME + ADC + Molecular subtype was evaluated using receiver operating characteristic (ROC) curves and the areas under the curve (AUC). A nomogram was built based on the predictive factors and internally evaluated using a bootstrap resampling method (1000 bootstrap resamples). The performance of the predictive model was evaluated by calibration curve and decision curve analysis (DCA). The DeLong test was applied to compare differences between AUC values, while net reclassification improvement (NRI) and integrated discrimination improvement (IDI) were used to assess the predictive ability of adding the Distribution of NME to the basic model [apparent diffusion coefficient (ADC) + Molecular subtype)].
Results
Compared to focal distribution, patients with linear distribution of NME had a higher risk of LVI positivity (P = 0.030). Distribution of NME + ADC + Molecular subtype demonstrated a relatively strong ability to predict LVI status, with an AUC of 0.723. Compared to the performance of each risk factor alone in predicting LVI, the differences in AUC were statistically significant (P = 0.008, P = 0.006, P = 0.012, DeLong test). Additionally, the inclusion of Distribution of NME could effectively improve the ability of basic model (ADC + Molecular subtype) to predict LVI, its NRI value was 0.389 (P = 0.013) and its IDI value was 0.047 (P = 0.008).
Conclusion
Distribution of NME + ADC + Molecular subtype was effective in predicting LVI status, with an AUC of 0.723. The inclusion of Distribution of NME significantly improved its predictive ability for LVI.
lymphovascular invasionmolecular subtypenon-mass enhancementDCE-MRIinvasive breast cancerNational Natural Science Foundation of China10.13039/501100001809National Key Laboratory Foundation of China10.13039/501100014206section-in-acceptanceBreast CancerIntroduction
According to GLOBOCAN 2022, breast cancer remains the most common cancer among women worldwide, and imaging examinations play a crucial role in the early detection, diagnosis, and treatment of breast cancer. Dynamic Contrast-Enhanced Magnetic Resonance Imaging (DCE-MRI) has the highest diagnostic accuracy among breast imaging techniques and can be used for preoperative evaluation, treatment planning, and assessing chemotherapy response in breast lesions (1). According to the fifth edition of the Breast Imaging Reporting and Data System (BI-RADS) lexicon, lesions on DCE-MRI are classified into mass enhancement (ME) and non-mass enhancement (NME) (2). NME is a unique imaging feature, characterized by enhancement within fibro glandular tissue without a mass effect, and occurs less frequently than ME.
Lymphovascular invasion (LVI) is a precursor of regional lymph node metastasis in invasive breast cancer and an independent prognostic factor for distant metastasis (3–5). It is significantly associated with poor prognosis and thus plays a crucial role in influencing treatment decisions for breast cancer. LVI is typically detected in postoperative pathological specimens, but preoperative treatments can affect the evaluation of it (6).
Currently, the vast majority of research on NME lesions focuses on predicting malignancy or exploring ways to reduce unnecessary biopsies, there are still few studies examining the factors influencing LVI in NME-type breast cancer (7). Conventional LVI prediction models (e.g., radiomic or clinicopathological frameworks) primarily rely on morphological features of mass-type lesions—such as spiculated margins and enhancement heterogeneity. These characteristics manifest as regional distribution patterns rather than focal masses in NME-type lesions, resulting in significantly diminished predictive performance of traditional models (8, 9). Standard MRI radiomic approaches depend on texture features (e.g., gray-level co-occurrence matrices), yet the diffuse distribution of NME lesions compromises texture feature extraction due to low signal-to-noise ratios and substantial interference from background parenchymal enhancement (BPE) (10).
Therefore, our study will investigate the associations between DCE-MRI imaging features, clinicopathological features, and LVI in NME-type invasive breast cancer patients. We will establish a dedicated NME-specific LVI prediction framework, aiming to provide clinical decision support tools for personalized surgical navigation and evidence-based treatment optimization.
Materials and methodsPatients
This study was a retrospective case-control study, involving 192 patients with non-mass enhancement (NME) invasive breast cancer who underwent DCE-MRI between January 2019 and December 2023. Based on postoperative pathological results, the patients were divided into two groups: LVI-positive [LVI(+)] (N = 50) and LVI-negative [LVI(-)] (N = 142). Inclusion criteria included: a) Female patients aged ≥18 years with pathologically confirmed invasive breast cancer after surgery; b) Patients who underwent DCE-MRI at our center prior to surgery, and whose imaging showed non-mass enhancement; c) Patients who had not received systemic or local treatments, including neoadjuvant therapy, radiotherapy, biopsy or surgery, within two years before the DCE-MRI examination. Exclusion criteria included: a) Patients with inflammatory breast cancer; b) Patients lacking key imaging parameters, clinical-pathological information, or follow-up data; c) Patients with concurrent other malignancies; d) Patients with pregnancy or a history of metformin use. This study was approved by the Ethics Review Committee of the hospital.
MRI examinations
A 3.0-Tesla MRI scanner with an 8-channel phased array breast coil was used to conduct the MRI examinations (Discovery 750, GE HealthCare, Chicago, IL, USA). The imaging protocol, lasting for 18 minutes, included 4 pulse sequences: diffusion-weighted imaging (DWI), T2-weighted imaging (T2WI), T1-weighted imaging (T1WI), and DCE. All sequences were spatially matched in axial view, with a field of view of 320 mm × 320 mm and 190 mm of z-axis coverage. The b value of DWI was 0 and 1,000 s/mm2 in 3 orthogonal diffusion gradients, the inversion recovery (IR) was 250 ms for fat suppression, the repetition time (TR) was 5,400 ms, minimum echo time (TE), and the matrix size was 128×128. The T2WI used iterative decomposition of water and fat with echo asymmetric and least squares estimation (IDEAL) for fat suppression, with a TR of 5,000 ms, a TE of 68 ms, and a matrix size of 320×256. The T1WI and DCE had the same geometric location and were performed with the following parameters: an isotropic spatial resolution of 1.0 mm × 1.0 mm × 1.0 mm, 192 partitions in the axial view, a minimum TR/TE, and a flip angle of 120°. The DCE scan repeated 6 continuous phases without interruption, each of which lasted 120 seconds (s). Following the pre-contrast phase, 0.5 M of gadopentetic acid (Gd-DTPA) was administered into the antecubital vein at a rate of 2 mL/s and a dose of 0.1 mmol/kg body weight and was followed by a 20-mL flush of saline. The time-intensity curve (TIC) was classified using the enhancement ratio difference between the first-enhanced phase and the last-enhanced phase (8–10 minutes after contrast delivery). The persistent, plateau and washout TICs were defined as enhancement ratios of <−10%, −10% to 10%, and >10%, respectively.
DCE-MRI imaging features
The study used a blinded diagnostic approach. Each patient’s imaging data were independently reviewed and evaluated by two radiologists specializing in breast diseases, with 8 and 12 years of experience, respectively, and formal professional training. Both radiologists independently interpreted and assessed the DCE-MRI images without knowledge of the patient’s clinical data, other test results, or treatment measures.
The description of imaging features followed the fifth edition of the BI-RADS lexicon (2).NME refers to enhancement detected in fibroglandular tissue regions on DCE-MRI without a mass effect. NME is categorized based on internal enhancement patterns (homogeneous, heterogeneous, clustered ring, and clumped enhancement) and distribution patterns (focal, linear, segmental, regional, multiple regions, and diffuse enhancement). Background parenchymal enhancement (BPE) refers to the enhancement of normal breast parenchyma and is classified as mild, moderate, or marked. The apparent diffusion coefficient (ADC) is a parameter describing the rate and extent of water molecule diffusion in DWI. The rate of early enhancement is calculated using the formula: (post-enhancement signal intensity - pre-enhancement signal intensity) × 100%/pre-enhancement signal intensity.
Histopathological evaluation
All histopathological assessments were performed on surgical specimens. Rosen’s histopathological criteria defined LVI as the presence of tumor emboli in the lymphatic and vascular systems. Immunohistochemistry (IHC) Staining: Estrogen receptor (ER) high expression is defined as ≥10% of tumor cell nuclei showing positive staining, while low or no expression is defined as <10% of nuclei staining positive. Progesterone receptor (PR) high expression is similarly defined as ≥10% positive staining in tumor cell nuclei, with low or no expression defined as <10%. Ki-67 is considered highly expressed if >14%, and low expression if ≤14%.
Human epidermal growth factor receptor 2 (HER2) expression is classified based on the proportion and intensity of membrane staining: “Negative” (“-”) refers to no staining or ≤10% of invasive cancer cells showing incomplete and faint membrane staining; “1+” refers to >10% of invasive cancer cells with incomplete and weak membrane staining; “2+” indicates >10% of invasive cancer cells with weak to moderate complete membrane staining or ≤10% of cells with strong complete membrane staining. For patients with “2+” HER2 expression, fluorescence in situ hybridization (FISH) is performed to determine HER2 gene amplification; “3+” indicates >10% of invasive cancer cells with strong, complete, and uniform membrane staining. HER2 status is defined as follows (11):
HER2-0: IHC score of “0”
HER2-low: IHC score of “1+” or “2+” with a negative FISH result
HER2-positive: IHC score of “3+” or “2+” with a positive FISH result
According to the 2017 St. Gallen International Breast Cancer Conference consensus (12), breast cancer is classified into four subtypes. Luminal A: ER/PR positive, HER2 negative, and Ki-67 ≤14%; Luminal B: Includes HER2 negative type (ER/PR positive, HER2 negative, and Ki-67 >14%) and HER2 positive type (ER/PR positive, HER2 positive, and any Ki-67 level); HER2 overexpression: ER/PR negative, HER2 positive, and any Ki-67 level; Triple-negative: ER, PR, and HER2 all negative, with any Ki-67 level.
Statistical analysis
Data analysis was performed using R version 4.4.0. All statistical tests were two-sided. Continuous variables were expressed as mean ± standard deviation (SD) or as median and interquartile range (IQR), depending on data distribution. Comparisons of continuous variables between the two groups were conducted using the Student’s t-test or the Mann-Whitney U test. Categorical variables were presented as counts and percentages, and comparisons between the two groups were made using the chi-square (χ2) test or Fisher’s exact test. A generalized variance inflation factor (GVIF) analysis was used to identify variables with high multicollinearity. Risk factors for LVI formation in patients with invasive breast cancer were analyzed using multivariate logistic regression. The diagnostic efficacy of each risk factor was evaluated by receiver operating characteristic (ROC) curves and the area under the curve (AUC). A nomogram was built based on the predictive factors and internally evaluated using a bootstrap resampling method (1000 bootstrap resamples). The performance of the predictive model was evaluated by calibration curve. Applying decision curve analysis (DCA) to evaluate the clinical usefulness of the model. DeLong test was used to compare differences between AUC values. Net reclassification improvement (NRI) and integrated discrimination improvement (IDI) were used to assess the predictive ability and clinical benefit of adding the Distribution of NME to the basic model (ADC + Molecular subtype). P-value below 0.05 was considered statistically significant.
ResultsBaseline characteristics
A total of 192 breast cancer lesions from 192 patients were included in this study. Of the patients, 73 (38%) were younger than 45 years, and 177 (92.2%) had invasive ductal carcinoma as their pathological type, other non-ductal forms of invasive breast cancer were classified as papillary carcinoma (10 cases), mucinous carcinoma (3 cases), and invasive lobular carcinoma (2 cases). Due to the limited overall number of cases for these subtypes, they were collectively grouped under the category “Others” for statistical analysis. The median size of the invasive cancers in all patients was 2.0 cm (range: 1.5–3.0 cm). Baseline characteristics are shown in Table 1.
Clinical-pathological and DCE-MRI features of patients with breast cancer between LVI(+) and LVI(-) (n, %).
Features
LVI(-) N=142
LVI(+) N=50
All N=192
P value
Clinicopathological features
Age
0.028
<45 yrs
47 (33.1)
26 (52.0)
73 (38.0)
≥45 yrs
95 (66.9)
24 (48.0)
119 (62.0)
Tumour type
1.000
IDC
131 (92.3)
46 (92.0)
177 (92.2)
Others
11 (7.7)
4 (8.0)
15 (7.8)
Tumor site
0.327
Left
73 (51.41%)
21 (42.00%)
94 (48.96%)
Right
69 (48.59%)
29 (58.00%)
98 (51.04%)
Invasive carcinoma size [M(Q1~Q3), cm]
2.0 (1.4~3.0)
2.7 (1.5~4.0)
2.0 (1.5~3.0)
0.007
pT (TNM)
0.030
1
81 (57.1)
21 (42.0)
102 (53.1)
2
53 (37.3)
20 (40.0)
73 (38.0)
3
6 (4.2)
8 (16.0)
14 (7.3)
4
2 (1.4)
1 (2.0)
3 (1.6)
ER status
0.067
Low
33 (23.2)
19 (38.0)
52 (27.1)
High
109 (76.8)
31 (62.0)
140 (72.9)
PR status
0.142
Low
29 (20.4)
16 (32.0)
45 (23.4)
High
113 (79.6)
34 (68.0)
147 (76.6)
Ki-67 status
0.013
Low
29 (20.4)
2 (4.0)
31 (16.2)
High
113 (79.6)
48 (96.0)
161 (83.8)
HER2 status
0.649
HER2-0
14 (9.8)
3 (6.0)
17 (8.8)
HER2-low
88 (62.0)
30 (60.0)
118 (61.5)
HER2-positive
40 (28.2)
17 (34.0)
57 (29.7)
Molecular subtype
0.019
Luminal A
28 (19.7)
2 (4.0)
30 (15.6)
Luminal B
88 (62.0)
33 (66.0)
121 (63.0)
HER2-overexpression
18 (12.7)
9 (18.0)
27 (14.1)
Triple-negative
8 (5.6)
6 (12.0)
14 (7.3)
Nerve invasion
1.000
Negative
136 (95.8)
48 (96.0)
184 (95.8)
Positive
6 (4.2)
2 (4.0)
8 (4.2)
SBR grade
0.696
1
13 (9.2)
4 (8.0)
17 (8.8)
2
113 (79.6)
38 (76.0)
151 (78.7)
3
16 (11.2)
8 (16.0)
24 (12.5)
ALN status
<0.001
Negative
87 (61.3)
10 (20.0)
97 (50.5)
Positive
55 (38.7)
40 (80.0)
95 (49.5)
pN (TNM)
<0.001
0
87 (61.3)
10 (20.0)
97 (50.5)
1
32 (22.5)
14 (28.0)
46 (24.0)
2
12 (8.4)
9 (18.0)
21 (10.9)
3
11 (7.8)
17 (34.0)
28 (14.6)
DCE-MRI features
Internal enhancement
Homogeneous
2 (1.4)
0 (0.0)
2 (1.0)
1.000
Heterogeneous
48 (33.8)
23 (46.0)
71 (37.0)
0.172
Consecutive nodular
32 (22.5)
12 (24.0)
44 (22.9)
0.987
Agglomerated annularly
89 (62.7)
28 (56.0)
117 (60.9)
0.507
Distribution
0.070
Focal
16 (11.3)
2 (4.0)
18 (9.4)
Linear
4 (2.8)
3 (6.0)
7 (3.7)
Segmental
83 (58.5)
24 (48.0)
107 (55.7)
Regional
18 (12.7)
8 (16.0)
26 (13.5)
Multiple regional
11 (7.8)
3 (6.0)
14 (7.3)
Diffuse
10 (7.0)
10 (20.0)
20 (10.4)
BPE
0.697
Minimal
28 (19.7)
11 (22.0)
39 (20.3)
Mild
65 (45.8)
19 (38.0)
84 (43.8)
Moderate
29 (20.4)
10 (20.0)
39 (20.3)
Marked
20 (14.1)
10 (20.0)
30 (15.6)
ADC value[M(Q1~Q3),×10–3 mm2/s]
0.99 (0.87~1.09)
0.92 (0.79~1.03)
0.97 (0.84~1.06)
0.027
TIC type
0.532
Persistent
6 (4.2)
1 (2.0)
7 (3.6)
Plateau
42 (29.6)
19 (38.0)
61 (31.8)
Washout
94 (66.2)
30 (60.0)
124 (64.6)
Rate of early enhancement
0.902
>120%
99 (69.7)
36 (72.0)
135 (70.3)
≤120%
43 (30.3)
14 (28.0)
57 (29.7)
LVI, lymphovascular invasion; IDC, invasive ductal carcinoma; pT, Pathological T stage of TNM; ER, estrogen receptor; PR, progesterone receptor; HER2, human epidermal growth factor receptor 2; SBR, Scarff-Bloom-Richardson; ALN, axillary lymph node; pN, Pathological N stage of TNM; BPE, background parenchyma enhancement; ADC, apparent diffusion coefficient; TIC, time-intensity curve.
Compared to the LVI(-) group, the LVI(+) group had a higher proportion of patients under 45 years of age (52% vs. 33.1%, P = 0.028), larger tumors (2.7 cm vs. 2.0 cm, P = 0.007), and a statistically significant difference in pathological T stage (P = 0.030). In terms of IHC, the proportion of high Ki-67 expression was higher in the LVI(+) group (96.0% vs. 79.6%, P = 0.013), while no statistically significant differences were observed in ER status (P = 0.067), PR status (P = 0.142), or HER2 status (P = 0.649). When comparing molecular subtypes, the proportion of Luminal A subtype was lower in the LVI(+) group (4.0% vs. 19.7%, P < 0.019), while the proportion of triple-negative subtype was higher (12.0% vs. 5.6%, P < 0.019). Additionally, the LVI(+) group had a significantly higher rate of axillary lymph node positivity (80.0% vs. 38.7%, P < 0.001), and there was a statistically significant difference in pathological N stage (P < 0.001). No statistically significant differences were observed between the two groups in terms of neural invasion or SBR grading (P = 1.000; P = 0.696) (Table 1).
DCE-MRI features analysis
Compared to the LVI(-) group, the LVI(+) group had a lower median ADC value (0.92×10−³ mm²/s vs. 0.99×10−³ mm²/s, P = 0.027). No statistically significant differences were observed between the two groups in terms of internal enhancement patterns (P = 1.000; P = 0.172; P = 0.987; P = 0.507), distribution patterns (P = 0.070), BPE (P = 0.697), TIC types (P = 0.532), or early enhancement rate (P = 0.902) (Table 1).
Risk factors of LVI in patients with invasive breast cancer
Table 2 presents the GVIF analysis used to assess multicollinearity. The regression model demonstrated an absence of significant multicollinearity, as all adjusted GVIF values were below 2.24 except for PR status. Although the GVIF for PR status exceeded 2.24, it remained below the 3.16 tolerance threshold, indicating acceptable collinearity levels. After controlling for confounding factors through multivariate logistic regression analysis, four risk factors associated with LVI formation were identified (Table 3). Patients aged ≥45 years (OR = 0.406, 95% CI: 0.191-0.844, P = 0.017) and those with higher ADC values (OR = 0.133, 95% CI: 0.017-0.825, P = 0.041) had a lower risk of LVI formation. Among the NME distribution patterns, patients with linear distribution had a significantly higher risk of LVI positivity compared to those with focal distribution (OR = 13.540, 95% CI: 1.390-172.644, P = 0.030) (Figures 1, 2), while no statistically significant differences were found for segmental, regional, multiple region, or diffuse enhancement patterns (all P > 0.05). Regarding molecular subtypes, compared to the Luminal A subtype, patients with Luminal B (OR = 5.081, 95% CI: 1.326-33.859, P = 0.039), HER2 overexpression (OR = 9.378, 95% CI: 1.922-71.460, P = 0.012), and triple-negative subtype (OR = 11.599, 95% CI: 2.043-96.710, P = 0.010) had a significantly higher risk of LVI formation.
DCE-MRI images of a 47-year-old woman with invasive ductal carcinoma of the right breast. Axial contrast-enhanced image (A) and sagittal contrast-enhanced image (B) show non-mass-like enhancement with the linear distribution. Lymphovascular invasion was confirmed at the histopathological examination.
Magnetic resonance images of a breast showing invasive cancer in two views. Image (A) displays an axial-sectional view with an arrow pointing to a specific area. Image (B) shows a side view with an arrow indicating the linear distribution.
DCE-MRI images of a 65-year-old woman with invasive ductal carcinoma of the left breast. Axial contrast-enhanced image (A) and sagittal contrast-enhanced image (B) show non-mass-like enhancement with the focal distribution. Lymphovascular invasion was negative at the histopathological examination.
MRI scans of breast tissue. Image (A) displays an axial-sectional view and Image (B) presents a sagittal view，they show the non-mass-like enhancement with the focal distribution.Diagnostic performance of each risk factor for predicting LVI
We evaluated the diagnostic performance of each risk factor for predicting LVI using ROC curves, and the DeLong test was applied to compare differences in AUC values. Figure 3 and Table 4 present the diagnostic performance of Distribution of NME + ADC + Molecular subtype and other individual risk factors for predicting LVI.
ROC curve of each risk factor predicting LVI.
ROC curve graph with sensitivity on the y-axis and 1-specificity on the x-axis. Four lines represent: Distribution of NME+ADC+Molecular Subtype (red), Distribution of NME (green), Molecular Subtype (blue), and ADC (purple). The diagonal line indicates random classification.
Comparison of AUC values of each risk factor for predicting LVI.
Risk factors
AUC
95% CI
P value
Distribution of NME + ADC + Molecular subtype
0.723
0.645-0.802
Ref.
Distribution of NME
0.625
0.540-0.709
0.008
Molecular subtype
0.614
0.540-0.687
0.006
ADC
0.606
0.513-0.698
0.012
LVI, lymphovascular invasion; AUC, Area under the curve; NME, non-mass enhancement; ADC, apparent diffusion coefficient.
As shown in Table 4, when using only the Distribution of NME to predict LVI, the AUC was 0.625. However, when ADC and Molecular subtype were added, the AUC significantly increased to 0.723 (P = 0.008, DeLong test). Similarly, the AUCs for Molecular subtype and ADC alone were 0.614 and 0.606, respectively. When compared to the AUC of Distribution of NME + ADC + Molecular subtype, the differences were statistically significant (P = 0.006, P = 0.012, DeLong test). To control Type I error across the 3 comparisons, Bonferroni correction was applied with an adjusted significance threshold of α = 0.0167. All post-comparison P values remained below this adjusted α, indicating statistically significant differences after correction. The bootstrap-derived AUC (1,000 iterations) showed high concordance with the apparent AUC upon internal validation (Figure 4). The calibration plot revealed good predictive accuracy between the actual probability and predicted probability (Figure 5). The DCA demonstrates that the predictive model yields net benefit across clinically relevant threshold probabilities (Figure 6). These indicate the potential utility of Distribution of NME + ADC + Molecular subtype in predicting LVI. Additionally, compared to the basic model (ADC + Molecular subtype), the inclusion of Distribution of NME improved the model’s reclassification ability, with a NRI of 0.389 (95% CI: 0.083-0.696) and an IDI of 0.047 (95% CI: 0.012-0.083), as shown in Table 5. These findings demonstrate that the inclusion of Distribution of NME significantly enhances the predictive ability of the model for LVI.
Resampling-corrected ROC for overfitting.
Bootstrap ROC curve showing the true positive rate versus false positive rate. The original ROC curve is in blue with an AUC of 0.723 and a 95% confidence interval of 0.645 to 0.802. Multiple gray lines represent bootstrap ROC iterations with a mean AUC of 0.726 and a standard deviation of 0.039.
Calibration curves for the model integrating Distribution of NME + ADC + Molecular subtype (A) and resampling-corrected calibration curve for overfitting (B).
Graph A displays a calibration plot of actual vs. predicted probabilities with a diagonal line for ideal calibration, a solid line for logistic calibration, and a dotted line for nonparametric. Statistical metrics like Dxy, C (ROC), and Brier score are noted. Graph B shows a calibration plot of observed vs. predicted probabilities with a dashed line for ideal, a solid line for bias-corrected, and a dotted line for apparent calibration. B equals one thousand repetitions, boot, with a mean absolute error of 0.027 for n equals 192.
DCA of the Distribution of NME + ADC + Molecular subtype prediction model. The y-axis represents net benefit. The solid red line denotes the predictive model incorporating Distribution of NME + ADC + Molecular subtype. The solid gray line indicates the ‘treat-all’ strategy (assuming all patients belong to the LVI-positive cohort). The solid black line represents the ‘treat-none’ strategy (assuming no patients are included in the LVI-positive cohort). The x-axis displays threshold probability, defined as the probability at which the expected benefit of treatment equals the expected benefit of avoiding treatment. The DCA demonstrates that the predictive model yields net benefit across clinically relevant threshold probabilities.
Decision curve analysis graph showing net benefit on the y-axis and high risk threshold on the x-axis, ranging from zero to one. The red line represents the model, the gray line represents none, and the black horizontal line represents all. The model demonstrates that the predictive model yields net benefit across clinically relevant threshold probabilities.
Model performance after the addition of Distribution of NME to the basic model (ADC + Molecular subtype) for predicting LVI.
Models
NRI
95% CI
P value
IDI
95% CI
P value
ADC + Molecular subtype
–
–
–
–
–
–
ADC + Molecular subtype +Distribution of NME
0.389
0.083-0.696
0.013
0.047
0.012-0.083
0.008
NRI, net classification improvement; IDI, integrated discrimination improvement; CI, confidence interval.
Discussion
The lesions of NME-type breast cancer are loose, often interspersed with normal fibroglandular tissue, making clinical palpation challenging. The pathology of NME-type breast cancer is diverse, and it lacks typical imaging features. Mammography and ultrasound have low sensitivity in detecting NME-type lesions, whereas DCE-MRI, with its high soft tissue and spatial resolution, allows for multiparametric and multidirectional imaging. This provides a detailed view of the morphological and hemodynamic characteristics of NME-type breast cancer, offering higher sensitivity. Due to the scattered nature of NME lesions, clinicians often approach breast-conserving surgery for NME-type breast cancer with caution, as these lesions may increase the risk of positive surgical margins, potentially affecting patient prognosis (13). Moreover, the St. Gallen Consensus reported that breast cancer patients with LVI positivity have a higher rate of local recurrence after breast-conserving surgery (14). Similarly, a study by LIU J et al. (13) found that LVI positivity is an independent predictor of positive surgical margins and an increased risk of local recurrence in breast-conserving surgery. Therefore, we speculate whether there are associated factors that may contribute to the increased risk of LVI in patients with NME-type breast cancer.
At present, the vast majority of research on NME lesions focuses on predicting malignancy or exploring ways to reduce unnecessary biopsies. As a highly specific form of breast cancer on DCE-MRI, the prognostic significance of NME-type breast cancer requires further investigation. LVI is a known risk factor for both local recurrence and distant metastasis in breast cancer patients, and it is also associated with poorer disease-free survival (DFS) and overall survival (OS) (15). Many studies have also identified other factors that influence breast cancer prognosis, such as age, axillary lymph node involvement, tumor histological grade, hormone receptor (HR) status, HER2 expression, and Ki-67 expression (16, 17). Some studies have reported correlations between LVI status and these established prognostic factors in breast cancer patients (18), however, results across different studies are not entirely consistent (19–21). Currently, there is limited research investigating the factors influencing LVI in NME-type breast cancer. This may be due to the relatively low incidence of NME, accounting for approximately 21% of invasive breast cancer cases (22). Our study focuses on this population, specifically exploring the relationship between DCE-MRI imaging features, clinical-pathological features, and LVI in patients with NME-type breast cancer.
In addition to local invasion, one of the key characteristics of cancer is its ability to spread, leading to the metastasis of cancer cells. Metastatic disease is the primary cause of death in cancer patients (23). Before cancer cells metastasize to secondary sites, they often first enter and spread through the lymphatic system. Therefore, LVI and minor perineural and neural invasion, is one of the biological prerequisites for systemic spread and metastasis (15). A study by Liu G et al. (24) reviewed the medical records of 14,782 unilateral breast cancer patients who underwent surgery at the Cancer Hospital of the Chinese Academy of Medical Sciences (CHCAMS) between January 2009 and September 2017. The study found that younger age was associated with poorer overall survival (OS), though the association with disease-free survival (DFS) was not significant. However, other studies have suggested that younger age is related to worse DFS (25, 26). These inconsistencies may be due to differences in histopathological types and staging. In our study, a higher proportion of patients in the LVI(+) group were younger than 45 years (52.0% vs. 33.1%, P = 0.028). Logistic multivariate regression analysis revealed that patients aged ≥45 years had a lower risk of LVI formation (OR = 0.406, 95%CI: 0.191-0.844, P = 0.017), which is consistent with previous studies (27, 28). It has been well-established that LVI positivity is a precursor for regional lymph node metastasis and an independent prognostic factor for distant metastasis in invasive breast cancer (3–5). Moreover, a study by Lee SJ et al. (29) demonstrated that LVI is also associated with higher recurrence rates and lower survival after surgery. Taken together with previous research and the findings from our study, these results suggest that the presence of LVI in younger breast cancer patients may be one of the factors contributing to their relatively poorer prognosis.
In terms of clinical-pathological features, CHOI BB’s study (30) suggested that histological grade, Ki-67 expression, and HR status were not correlated with LVI. However, in our study, the proportion of high Ki-67 expression was significantly higher in the LVI(+) group compared to the LVI(-) group (96.0% vs. 79.6%, P = 0.013). This discrepancy may be due to the smaller sample size in CHOI BB’s study, which included only 132 patients. On the other hand, a study by ABE H et al. (31) yielded results similar to ours, finding that LVI was associated with high tumor cell proliferation and high Ki-67 expression. Our study also found that, compared to the Luminal A subtype, patients with Luminal B, HER2 overexpression and triple-negative breast cancer had a higher risk of LVI formation. Previous studies have reported that LVI is more common in HR-negative, HER2 positive, and high Ki-67 expression tumors (32–34). According to the major intrinsic biological subtype of breast cancer, LVI is more common in HER2-overexpression and Luminal B subtypes (32, 34), which is similar to the results of our study. Furthermore, Van den Eynden et al. (35) found that LVI was significantly associated with lymph node metastasis, a conclusion also supported by SHEN SD et al. (18). In our study, the LVI(+) group had a higher rate of axillary lymph node positivity (80.0% vs. 38.7%, P < 0.001). Axillary lymph node metastasis often indicates an increased risk of distant metastasis, which helps explain why LVI is a risk factor for distant metastasis and is associated with poorer DFS and OS outcomes (13). Therefore, for patients diagnosed with invasive breast cancer via core needle biopsy, assessing HR, HER2 status, and Ki-67 expression in the biopsy tissue can help clinicians preliminarily determine the molecular subtype. This enables clinicians to be vigilant about the potential presence of LVI and select the more appropriate treatment strategy.
Most current research on NME detected by MRI focuses on distinguishing between benign and malignant lesions. Malignant breast tumors typically exhibit low ADC values, which is attributed to the increased number of cells, larger nuclei, abundant macromolecular proteins, and reduced extracellular space in malignant tumors, all of which significantly reduce the water molecule diffusion capacity. Additionally, ADC values may be affected by factors such as hemorrhage within the lesion, highly dense fluids, or enlarged lymph nodes in the breast (36). Considering the intermingling of tumors and normal tissues in NME lesions, the ADC value of lesions was obtained by selecting multiple points for the measurement and taking the average value to avoid bias in the data as much as possible. Some studies have shown that ADC values are related to tumor invasiveness in breast cancer patients and are associated with prognosis (37–39), other researches have demonstrated a negative correlation between ADC values and LVI (21, 40, 41). Retrospective cohort studies by DURANDO M et al. (42) and IGARASHI T et al. (43) both found that ADC values were significantly lower in the LVI-positive group compared to the LVI-negative group, which is consistent with our findings.
A study by CHIKARMANE SA et al. (44) found that the most common distribution pattern of NME-type breast cancer on DCE-MRI was focal (49.7%; 102/205), followed by linear (25.6%; 53/205). However, other studies have indicated that segmental distribution was more common (45, 46). The discrepancy in findings may be attributed to differences in study populations and sample sizes. In our study, regardless of LVI status, the most common distribution pattern was segmental. The reason for this result may be related to the growth of the tumor along the duct. After multivariate adjustment, we found that patients with linear distribution of NME had a significantly higher risk of LVI positivity compared to those with focal distribution (P = 0.030). Although the total number of cases with linear distribution in our data was small (N = 7), this result suggests a possible correlation between lesion distribution patterns and LVI formation in NME patients. In the future, we plan to collect more samples to further validate the reliability of this conclusion. To further investigate whether distribution patterns play a role in predicting LVI, we applied NRI and IDI for evaluation. The results demonstrated that the inclusion of Distribution of NME significantly improved the predictive ability of the basic model (ADC + Molecular subtype) for LVI [NRI (0.389; 95% CI: 0.083-0.696) and IDI (0.047; 95% CI: 0.012-0.083)]. Compared to the evaluation of LVI by individual risk factors, the AUC increased significantly to 0.723. Therefore, we believe that Distribution of NME may have potential utility in predicting LVI. Currently, no other studies have explored the relationship between distribution patterns of NME-type breast cancer and LVI, making our study somewhat novel in this regard. This may further alert clinicians to pay closer attention to NME-type breast cancer in order to develop more appropriate treatment strategies for patients at risk of LVI.
This study has several limitations. First, as a retrospective study, future prospective studies with larger sample sizes are needed to further validate our findings. Second, since this is a single-center study, there is a certain degree of selection bias in the enrolled patient population. Third, due to the relatively short follow-up period, we did not explore the relationship between NME and the prognosis of patients with invasive breast cancer. Lastly, the pathological types included in our study were limited to invasive carcinomas, meaning our sample does not fully represent the histological diversity of the general population.
Conclusions
Our study found that age ≥ 45 years and higher ADC values were associated with a lower risk of LVI formation. Among the NME distribution patterns, patients with linear distribution had a higher risk of LVI positivity compared to those with focal distribution. Regarding molecular subtypes, patients with Luminal B, HER2-overexpression, and triple-negative subtypes had a higher risk of LVI formation compared to those with Luminal A. Distribution of NME + ADC + Molecular subtype demonstrated relatively good predictive ability for LVI, with an AUC of 0.723, and the inclusion of Distribution of NME significantly improved its predictive capability. This study provides a reference for evaluating LVI during the preoperative imaging stage, contributing to clinical decision-making and personalized treatment.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving humans were approved by Ethics Committee Of Chinese PLA General Hospital. The studies were conducted in accordance with the local legislation and institutional requirements. The human samples used in this study were acquired from a by- product of routine care or industry. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.
Author contributions
SG: Data curation, Writing – original draft, Writing – review & editing, Formal Analysis, Methodology, Conceptualization, Investigation. KZ: Investigation, Writing – original draft, Writing – review & editing, Data curation, Formal Analysis, Methodology, Validation. HX: Data curation, Formal Analysis, Investigation, Methodology, Validation, Writing – original draft, Writing – review & editing. YX: Data curation, Methodology, Writing – original draft. QL: Data curation, Methodology, Writing – original draft. JL: Investigation, Methodology, Writing – original draft. SC: Investigation, Methodology, Writing – original draft. JS: Investigation, Methodology, Writing – original draft. ZG: Investigation, Methodology, Writing – original draft. LZ: Conceptualization, Methodology, Writing – review & editing. JW: Conceptualization, Investigation, Project administration, Resources, Supervision, Validation, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. This work was supported by the National Natural Science Foundation of China (No. 62476285) and Major Project of the Open Fund of Chinese State Key Laboratory (No. EMC2024N001).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
ReferencesMannRMChoNMoyL. Breast MRI: state of the art. Radiology. (2019) 292:520–36. doi: 10.1148/radiol.2019182947, PMID: 31361209SpakDAPlaxcoJSSantiagoLDrydenMJDoganBE. BI-RADS(®) fifth edition: A summary of changes. Diagn Interv Imaging. (2017) 98:179–90. doi: 10.1016/j.diii.2017.01.001, PMID: 28131457RagageFDebledMMacGroganGBrousteVDesrousseauxMSoubeyranI. Is it useful to detect lymphovascular invasion in lymph node-positive patients with primary operable breast cancer? Cancer. (2010) 116:3093–101. doi: 10.1002/cncr.25137, PMID: 20564641MaQDieterichLCDetmarM. Multiple roles of lymphatic vessels in tumor progression. Curr Opin Immunol. (2018) 53:7–12. doi: 10.1016/j.coi.2018.03.018, PMID: 29605736ZhangSZhangDGongMWenLLiaoCZouL. High lymphatic vessel density and presence of lymphovascular invasion both predict poor prognosis in breast cancer. BMC Cancer. (2017) 17:335. doi: 10.1186/s12885-017-3338-x, PMID: 28514957KayadibiYKocakBUcarNAkanYNYildirimEBektasS. MRI radiomics of breast cancer: machine learning-based prediction of lymphovascular invasion status. Acad Radiol. (2022) 29:S126–S34. doi: 10.1016/j.acra.2021.10.026, PMID: 34876340ZhangJLiLZhangLZheXTangMLeiX. Meta-analysis of dynamic contrast enhancement and diffusion-weighted MRI for differentiation of benign from Malignant non-mass enhancement breast lesions. Front Oncol. (2024) 14. doi: 10.3389/fonc.2024.1332783, PMID: 38544833YangXFanXLinSZhouYLiuHWangX. Assessment of lymphovascular invasion in breast cancer using a combined MRI morphological features, radiomics, and deep learning approach based on dynamic contrast-enhanced MRI. J Magn Reson Imaging. (2024) 59:2238–49. doi: 10.1002/jmri.29060, PMID: 37855421LiuZFengBLiCChenYChenQLiX. Preoperative prediction of lymphovascular invasion in invasive breast cancer with dynamic contrast-enhanced-MRI-based radiomics. J Magn Reson Imaging. (2019) 50:847–57. doi: 10.1002/jmri.26688, PMID: 30773770YangXWangXZuoZZengWLiuHZhouL. Radiomics-based analysis of dynamic contrast-enhanced magnetic resonance image: A prediction nomogram for lymphovascular invasion in breast cancer. Magn Reson Imaging. (2024) 112:89–99. doi: 10.1016/j.mri.2024.07.001, PMID: 38971267TarantinoPVialeGPressMFHuXPenault-LlorcaFBardiaA. ESMO expert consensus statements (ECS) on the definition, diagnosis, and management of HER2-low breast cancer. Ann Oncol. (2023) 34:645–59. doi: 10.1016/j.annonc.2023.05.008, PMID: 37269905CuriglianoGBursteinHJWinerEPGnantMDubskyPLoiblS. De-escalating and escalating treatments for early-stage breast cancer: the St. Gallen International Expert Consensus Conference on the Primary Therapy of Early Breast Cancer 2017. Ann Oncol. (2017) 28:1700–12. doi: 10.1093/annonc/mdx308, PMID: 28838210LiuJZhouDPanTCaoLYangQZhangJ. Socioeconomic and clinical factors affecting the proportion of breast conserving surgery in Chinese women with breast cancer. Gland Surg. (2022) 11:341–51. doi: 10.21037/gs-22-25, PMID: 35284315IgnatiadisMBuyseMSotiriouC. St Gallen International Expert Consensus on the primary therapy of early breast cancer: an invaluable tool for physicians and scientists. Ann Oncol. (2015) 26:1519–20. doi: 10.1093/annonc/mdv259, PMID: 26063634KuhnEGambiniDDespiniLAsnaghiDRunzaLFerreroS. Updates on lymphovascular invasion in breast cancer. Biomedicines. (2023) 11(3):968. doi: 10.3390/biomedicines11030968, PMID: 36979946SchwartzAMHensonDEChenDRajamarthandanS. Histologic grade remains a prognostic factor for breast cancer regardless of the number of positive lymph nodes and tumor size: a study of 161 708 cases of breast cancer from the SEER Program. Arch Pathol Lab Med. (2014) 138:1048–52. doi: 10.5858/arpa.2013-0435-OA, PMID: 25076293CuriglianoGBursteinHJGnantMLoiblSCameronDReganMM. Understanding breast cancer complexity to improve patient outcomes: The St Gallen International Consensus Conference for the Primary Therapy of Individuals with Early Breast Cancer 2023. Ann Oncol. (2023) 34:970–86. doi: 10.1016/j.annonc.2023.08.017, PMID: 37683978ShenSDZhongSZWangCZHuangWH. Correlation of lymphovascular invasion with clinicopathological factors in invasive breast cancer: a meta-analysis. Int J Clin Exp Med. (2015) 8:17789–95.UematsuT. Focal breast edema associated with Malignancy on T2-weighted images of breast MRI: peritumoral edema, prepectoral edema, and subcutaneous edema. Breast Cancer. (2014) 22:66–70. doi: 10.1007/s12282-014-0572-9, PMID: 25336185UematsuTKasamiMWatanabeJ. Is evaluation of the presence of prepectoral edema on T2-weighted with fat-suppression 3 T breast MRI a simple and readily available noninvasive technique for estimation of prognosis in patients with breast cancer? Breast Cancer. (2014) 21(6):684–92. doi: 10.1007/s12282-013-0440-z, PMID: 23341125MoriNMugikuraSTakasawaCMiyashitaMShimauchiAOtaH. Peritumoral apparent diffusion coefficients for prediction of lymphovascular invasion in clinically node-negative invasive breast cancer. Eur Radiol. (2015) 26:331–9. doi: 10.1007/s00330-015-3847-4, PMID: 26024846LeeSMNamKJChooKSKimJYJeongDWKimHY. Patterns of Malignant non-mass enhancement on 3-T breast MRI help predict invasiveness: using the BI-RADS lexicon fifth edition. Acta Radiol. (2018) 59:1292–9. doi: 10.1177/0284185118759139, PMID: 29758996MohammedRAMartinSGMahmmodAMMacmillanRDGreenARPaishEC. Objective assessment of lymphatic and blood vascular invasion in lymph node-negative breast carcinoma: findings from a large case series with long-term follow-up. J Pathol. (2011) 223:358–65. doi: 10.1002/path.2810, PMID: 21171081LiuGKongXDaiQChengHWangJGaoJ. Clinical features and prognoses of patients with breast cancer who underwent surgery. JAMA Netw Open. (2023) 6:e2331078. doi: 10.1001/jamanetworkopen.2023.31078, PMID: 37624596HanWKimSWParkIAKangDKimSWYounYK. Young age: an independent risk factor for disease-free survival in women with operable breast cancer. BMC Cancer. (2004) 4:82. doi: 10.1186/1471-2407-4-82, PMID: 15546499ChenBFTsaiYFHuangCCHsuCYLienPJWangYL. Clinicopathological characteristics and treatment outcomes of luminal B1 breast cancer in Taiwan. J Chin Med Assoc. (2022) 85:190–7. doi: 10.1097/JCMA.0000000000000632, PMID: 34643617KolliasJElstonCWEllisIORobertsonJFBlameyRW. Early-onset breast cancer–histopathological and prognostic considerations. Br J Cancer. (1997) 75:1318–23. doi: 10.1038/bjc.1997.223, PMID: 9155052ShannonCSmithIE. Breast cancer in adolescents and young women. Eur J Cancer. (2003) 39:2632–42. doi: 10.1016/S0959-8049(03)00669-5, PMID: 14642925LeeSJGoJAhnBSAhnJHKimJYParkHS. Lymphovascular invasion is an independent prognostic factor in breast cancer irrespective of axillary node metastasis and molecular subtypes. Front Oncol. (2023) 13:1269971. doi: 10.3389/fonc.2023.1269971, PMID: 38053656ChoiBB. Dynamic contrast enhanced-MRI and diffusion-weighted image as predictors of lymphovascular invasion in node-negative invasive breast cancer. World J Surg Oncol. (2021) 19:76. doi: 10.1186/s12957-021-02189-3, PMID: 33722246AbeHSchmidtRAKulkarniKSennettCAMuellerJSNewsteadGM. Axillary lymph nodes suspicious for breast cancer metastasis: sampling with US-guided 14-gauge core-needle biopsy–clinical experience in 100 patients. Radiology. (2009) 250:41–9. doi: 10.1148/radiol.2493071483, PMID: 18955508HouvenaeghelGCohenMClasseJMReyalFMazouniCChopinN. Lymphovascular invasion has a significant prognostic impact in patients with early breast cancer, results from a large, national, multicenter, retrospective cohort study. ESMO Open. (2021) 6:100316. doi: 10.1016/j.esmoop.2021.100316, PMID: 34864349RakhaEAMartinSLeeAHMorganDPharoahPDHodiZ. The prognostic significance of lymphovascular invasion in invasive breast carcinoma. Cancer. (2012) 118:3670–80. doi: 10.1002/cncr.26711, PMID: 22180017NishimuraROsakoTOkumuraYNakanoMOhtsukaHFujisueM. An evaluation of lymphovascular invasion in relation to biology and prognosis according to subtypes in invasive breast cancer. Oncol Lett. (2022) 24:245. doi: 10.3892/ol.2022.13366, PMID: 35761943Van den EyndenGGvan der AuweraIVan LaereSJColpaertCGvan DamPDirixLY. Distinguishing blood and lymph vessel invasion in breast cancer: a prospective immunohistochemical study. Br J Cancer. (2006) 94:1643–9. doi: 10.1038/sj.bjc.6603152, PMID: 16670715ImamuraTIsomotoISueyoshiEYanoHUgaTAbeK. Diagnostic performance of ADC for Non-mass-like breast lesions on MR imaging. Magn Reson Med Sci. (2010) 9:217–25. doi: 10.2463/mrms.9.217, PMID: 21187691MartincichLDeantoniVBertottoIRedanaSKubatzkiFSarottoI. Correlations between diffusion-weighted imaging and breast cancer biomarkers. Eur Radiol. (2012) 22:1519–28. doi: 10.1007/s00330-012-2403-8, PMID: 22411304ParkSHChoiHYHahnSY. Correlations between apparent diffusion coefficient values of invasive ductal carcinoma and pathologic factors on diffusion-weighted MRI at 3.0 Tesla. J Magn Reson Imaging. (2013) 41(1):175–82. doi: 10.1002/jmri.24519, PMID: 24353241ThakurSBDurandoMMilansSChoGYGennaroLSuttonEJ. Apparent diffusion coefficient in estrogen receptor-positive and lymph node-negative invasive breast cancers at 3.0T DW-MRI: A potential predictor for an oncotype Dx test recurrence score. J Magn Reson Imaging. (2017) 47(2):401–9. doi: 10.1002/jmri.25796, PMID: 28640531KaranBPourbagherATorunN. Diffusion-weighted imaging and 18F-fluorodeoxyglucose positron emission tomography/computed tomography in breast cancer: Correlation of the apparent diffusion coefficient and maximum standardized uptake values with prognostic factors. J Magn Reson Imaging. (2015) 43(6):1434–44. doi: 10.1002/jmri.25112, PMID: 26663655YangWQiangJWTianHPChenBWangAJZhaoJG. Minimum apparent diffusion coefficient for predicting lymphovascular invasion in invasive cervical cancer. J Magnetic Resonance Imaging. (2016) 45:1771–9. doi: 10.1002/jmri.25542, PMID: 27859961DurandoMGennaroLChoGYGiriDDGnanasigamaniMMPatilS. Quantitative apparent diffusion coefficient measurement obtained by 3.0Tesla MRI as a potential noninvasive marker of tumor aggressiveness in breast cancer. Eur J Radiol. (2016) 85:1651–8. doi: 10.1016/j.ejrad.2016.06.019, PMID: 27501902IgarashiTFurubeHAshidaHOjiriH. Breast MRI for prediction of lymphovascular invasion in breast cancer patients with clinically negative axillary lymph nodes. Eur J Radiol. (2018) 107:111–8. doi: 10.1016/j.ejrad.2018.08.024, PMID: 30292254ChikarmaneSAMichaelsAYGiessCS. Revisiting nonmass enhancement in breast MRI: analysis of outcomes and follow-up using the updated BI-RADS atlas. Am J Roentgenol. (2017) 209:1178–84. doi: 10.2214/AJR.17.18086, PMID: 28834447AydinH. The MRI characteristics of non-mass enhancement lesions of the breast: associations with Malignancy. Br J Radiol. (2019) 92:20180464s. doi: 10.1259/bjr.20180464, PMID: 30673299de Faria Castro FleuryECastroCdo AmaralMSCRoveda JuniorD. Management of non-mass enhancement at breast magnetic resonance in screening settings referred for magnetic resonance-guided biopsy. Breast Cancer (Auckl). (2022) 16:11782234221095897. doi: 10.1177/11782234221095897, PMID: 35602239