ORIGINAL RESEARCH article

Front. Oncol., 28 April 2023

Sec. Breast Cancer

Volume 13 - 2023 | https://doi.org/10.3389/fonc.2023.1139189

Histogram analysis of multi-model high-resolution diffusion-weighted MRI in breast cancer: correlations with molecular prognostic factors and subtypes

  • 1. Department of Radiology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China

  • 2. Department of Radiology, Xiangyang Central Hospital, Affiliated Hospital of Hubei University of Arts and Science, Xiangyang, China

  • 3. Department of Radiology, Xiantao First People’s Hospital Affiliated to Yangtze University, Xiantao, China

  • 4. Magnetic Resonance (MR) Scientific Marketing, Siemens Healthineers Ltd., Wuhan, China

  • 5. Magnetic Resonance (MR) Collaborations, Siemens Healthineers Ltd., Chengdu, China

Abstract

Objective:

To investigate the correlations between quantitative diffusion parameters and prognostic factors and molecular subtypes of breast cancer, based on a single fast high-resolution diffusion-weighted imaging (DWI) sequence with mono-exponential (Mono), intravoxel incoherent motion (IVIM), diffusion kurtosis imaging (DKI) models.

Materials and Methods:

A total of 143 patients with histopathologically verified breast cancer were included in this retrospective study. The multi-model DWI-derived parameters were quantitatively measured, including Mono-ADC, IVIM-D, IVIM-D*, IVIM-f, DKI-Dapp, and DKI-Kapp. In addition, the morphologic characteristics of the lesions (shape, margin, and internal signal characteristics) were visually assessed on DWI images. Next, Kolmogorov–Smirnov test, Mann-Whitney U test, Spearman’s rank correlation, logistic regression, receiver operating characteristic (ROC) curve, and Chi-squared test were utilized for statistical evaluations.

Results:

The histogram metrics of Mono-ADC, IVIM-D, DKI-Dapp, and DKI-Kapp were significantly different between estrogen receptor (ER)-positive vs. ER-negative groups, progesterone receptor (PR)-positive vs. PR-negative groups, Luminal vs. non-Luminal subtypes, and human epidermal receptor factor-2 (HER2)-positive vs. non-HER2-positive subtypes. The histogram metrics of Mono-ADC, DKI-Dapp, and DKI-Kapp were also significantly different between triple-negative (TN) vs. non-TN subtypes. The ROC analysis revealed that the area under the curve considerably improved when the three diffusion models were combined compared with every single model, except for distinguishing lymph node metastasis (LNM) status. For the morphologic characteristics of the tumor, the margin showed substantial differences between ER-positive and ER-negative groups.

Conclusions:

Quantitative multi-model analysis of DWI showed improved diagnostic performance for determining the prognostic factors and molecular subtypes of breast lesions. The morphologic characteristics obtained from high-resolution DWI can be identifying ER statuses of breast cancer.

Introduction

Diffusion-weighted imaging (DWI) has been proven to be a potential diagnostic tool for the evaluation of breast lesions (1). In clinical practice, single-shot echo planar imaging (ss-EPI) sequence is commonly utilized in breast DWI protocols due to fast acquisition time (2, 3). However, previous studies have seldom evaluated the morphologic analysis of lesions on ss-EPI images due to geometric distortion and poor spatial resolution (4, 5). Readout-segmented EPI (rs-EPI) can reduce distortions and maintain high resolution, but it was limited by long scan times (6). Recently, the simultaneous multislice (SMS) technique, which allows the exciting acquire multiple slices at once, has been introduced to reduce the scan time (7, 8). The SMS technique was combined with rs-EPI to generate images with less image distortion and higher spatial resolution for breast lesions in a clinically acceptable scanning duration (9, 10). The application of SMS rs-EPI makes it feasible to qualitatively assess the morphologic characteristics of breast lesions.

In addition to the qualitative analysis of SMS rs-EPI images, our study also focused on the quantitative analysis of multi-model DWI. Conventional DWI is based on a mono-exponential (Mono) model that was first proposed to reflect the random Brownian motion of water molecules diffusing into biological media by quantifying apparent diffusion coefficient (ADC) values (11). However, water diffusion in complex biological media may be influenced by the blood microcirculation in capillaries, leading to a non-Gaussian distribution (12, 13). To address this, advanced diffusion models, including intravoxel incoherent motion (IVIM) and diffusion kurtosis imaging (DKI), have been developed to reflect the diffusion behavior of water molecules in tumors more accurately (14). A few studies have investigated the correlations between IVIM- or DKI-derived parameters with several clinical prognostic factors and molecular subtypes. However, the conclusions have still not reached a consensus (15, 16). Furthermore, most studies have reported that more information can be parsed from histogram analysis, which can reflect the microstructures and heterogeneity of breast cancer (17–20).

Therefore, the aim of this study was to apply three diffusion models (Mono, IVIM, and DKI) to determine the parameters valuable for differentiating between prognostic factor statuses and molecular subtypes, as well as to assess the correlations of morphologic characteristics with prognostic factors and molecular subtypes.

Materials and methods

Patients

This retrospective study was approved by our institutional review board, and informed consent was obtained. From September 2020 to May 2021, 216 female patients who underwent breast MRI in our hospital and fulfilled the following criteria were selected: (1) the patients did not undergo chemotherapy, or any other interventions before they were examined by MRI; (2) the pathologic type of breast lesions was confirmed by surgery or biopsy; and (3) relevant pathologic data of patients were complete. The exclusion criteria included: 1) non-mass like enhancement lesions detected on dynamic contrast-enhanced (DCE)-MRI (n = 34); 2) the max diameter of mass lesions< 1cm (n = 16); 3) poor DWI image quality due to patient motion or susceptibility artifact (n = 23). Only the largest lesion was analyzed when multiple lesions were detected in the bilateral breast. Finally, 143 patients (mean age, 48.57 ± 12.01 years, range, 26 – 81 years) with 143 mass lesions (mean diameter, 2.48 ± 0.95 cm) were included in the study. More detailed characteristics of the 143 patients are summarized in Table 1.

Table 1

CharacteristicsN (%)
 Age at diagnosis≤ 5078 (54.5)
> 5065 (45.5)
 Long diameter (cm)> 291 (63.6)
≤ 252 (36.4)
 SideRight68 (47.6)
Left75 (52.4)
 Menopausal statusPremenopausal74 (51.7)
Postmenopausal69 (48.3)
 Histological typeIDC109 (76.2)
ILC11 (7.7)
Papillary carcinoma8 (5.6)
DCIS15 (10.5)
 ERPositive80 (55.9)
Negative63 (44.1)
 PRPositive76 (53.1)
Negative67 (46.9)
 HER2Positive48 (33.6)
Negative95 (66.4)
 Ki-67Positive79 (55.2)
Negative64 (44.8)
 LNMPositive49 (34.3)
Negative94 (65.7)
Molecular subtypes
 Luminal A45 (31.5)
 Luminal B41 (28.7)
 HER2-positive25 (17.5)
 Triple-negative32 (22.3)
Morphological features
 ShapeRound24 (16.8)
Oval49 (34.3)
Irregular70 (48.9)
 MarginSmooth76 (53.1)
Spiculated16 (11.2)
Irregular51 (35.7)
 SignalHomogenous55 (38.4)
Heterogenous64 (44.8)
Rim24 (16.8)

Study population and histopathological characteristics.

DCIS, Ductal carcinoma in situ; IDC, Invasive ductal carcinoma; ILC, Invasive lobular carcinoma; ER, Estrogen receptor; PR, Progesterone receptor; HER2, Human epidermal growth factor receptor 2; LNM, Lymph nodes metastasis; TN, Triple-negative.

MRI scans

Breast MRI was performed on a 3T MRI scanner (MAGNETOM Skyra, Siemens Healthcare, Erlangen, Germany) using a dedicated 16-channel phased-array bilateral breast surface coil. The breast MRI protocol included the following sequences: axial fat-saturated T2-weighted imaging, T1-weighted DCE-MRI with the time-resolved angiography (TWIST) with a volumetric interpolated breath-hold examination (VIBE) technique, and SMS rs-EPI sequence. Detailed imaging parameters are provided in Table 2.

Table 2

ParametersT2WISMS rs-EPIDCE-MRI
Repetition time (ms)370023505.24
Echo time (ms)101722.46
Field of view (mm2)320 x 320280 x 280320 x 320
Matrix224 x 320122 x 188182 x 320
Slice thickness (mm)451.5
Pixel bandwidth (Hz/Px)347887780
Parallel imagingGRAPPA (x2)GRAPPA (x2)CAIPIRINHA (x4)
b-values (sec/mm2)/0, 50, 100, 200, 400, 800, 1000, 2000/
Readout segment/5/
Multi-slice mode/Slice acceleration (x2)/
Temporal resolution (sec/phase)//5.74
Acquisition time (min:sec)2:064:395:57

Sequence parameters for T2-weighted imaging, SMS rs-EPI, and DCE-MRI.

SMS, Simultaneous multi-slice; DCE, Dynamic contrast-enhanced; GRAPPA, Generalized autocalibrating partially parallel acquisition; CAIPIRINHA, Controlled aliasing in parallel imaging results in higher acceleration.

/ indicates Non-applicable.

Image analysis

The images were independently analyzed by two breast readers (with 3 and 5 years of experience, respectively) using an in-house-developed DKI tool software. Both readers were informed that the patients had breast cancer but were blinded to the detailed pathologic data. Two-dimensional (2D) region of interest (ROI) were manually delineated, which excluded the cystic or necrotic portions of the tumor, on high-b-value (b=1000 s/mm2) SMS rs-EPI images, with the reference of the corresponding T2-weighted and DCE-MRI images. The ROI was then copied to other parametric maps [including ADC, pure diffusion (D), pseudo-diffusion coefficient (D*), perfusion fraction (f), apparent diffusional kurtosis (Kapp), and apparent diffusion coefficient (Dapp) maps] using the DKI tool software. Finally, the histogram information of each ROI map was generated, including mean, median, percentile values (25th and 75th), kurtosis, and skewness. For example, the mean and 75th percentile metrics of Mono-ADC were presented as Mono-ADCmean and Mono-ADC75th, respectively.

The corresponding mathematical expressions were as follows:

1. Mono-exponential model (16):

where Sb is the signal intensity on the DWI image at a certain b value (800 sec/mm2) and S0 is the signal intensity value in the voxels with b values of 0.

2.DKI model (21):

where Sb is he signal intensity on the DWI image according to all b-values (0, 50, 100, 200, 400, 800, 1000, and 2000 sec/mm2). Dapp represents the non-Gaussian diffusion coefficient and Kapp represents the apparent kurtosis coefficient without unit.

3. IVIM model (22):

where Sb is he signal intensity on the DWI image according to the b-value (0, 50, 100, 200, and 400 sec/mm2). D is the true diffusion coefficient representing the simple movement of water molecules in the tissue (unit: mm2/s), D* is the pseudo-diffusion coefficient representing perfusion-related diffusion (unit: mm2/s), and f is the fraction of fast diffusion representing the diffusion linked to microcirculation (0 ≤ f ≤ 1).

Morphologic analysis

Two experienced readers independently assessed several morphologic characteristics on SMS rs-EPI images with b = 1000 mm2/s according to the Breast Imaging Reporting and Data System lexicon (BI-RADS edition 2013). Since enhancement is mostly used to evaluate breast lesions on DCE-MRI, two readers evaluated breast lesions on DWI images using internal signal characteristics, which were defined as lesions with homogeneous signal, heterogeneous signal, or only high-signal at the rim (23, 24). Each morphological characteristic was specifically evaluated as follows:

  • Lesion shape: 1 for round, 2 for oval, 3 for irregular.

  • Lesion margin: 1 for smooth, 2 for spiculated, 3 for irregular.

  • Lesion internal signal characteristics: 1 for heterogeneous, 2 for homogeneous, 3 for rim.

  • The max diameter of lesion was measured on the largest tumor section.

Histopathologic assessment

Histopathologic results were obtained from the electronic medical records of each patient in our hospital. Estrogen receptor (ER) positivity and progesterone receptor (PR) positivity were defined as the presence of 1% or more positively stained nuclei in 10 high-power fields (25). Human epidermal growth factor receptor 2 (HER2) was considered positive if it was scored 3+ for immunohistochemically stained tissue, or gene amplification was observed with fluorescence in situ hybridization (FISH) (26). More than 20% of cancer nuclei were positively stained for Ki-67 (12). Lymph node metastasis (LNM) was confirmed by the clinician performing the histopathologic examination (13). According to the statuses of ER, PR, HER2, and Ki-67, the breast tumors were further classified as Luminal A, Luminal B, HER2-positive, and triple-negative (TN) (27).

Statistical analysis

All statistical analyses were performed using MedCalc software (version 15.0, Ostend, Belgium) and SPSS software (version 26.0, IL, USA). The inter-reader agreement for diffusion parameters and morphological characteristics was assessed by using the intraclass correlation coefficient (ICC): ICC ≤ 0.40, poor agreement; 0.40 - 0.59, fair agreement; 0.60 - 0.74, good agreement; 0.75 - 1.00, excellent agreement. The categorical variables were as follows: prognostic factors including ER, PR, HER2, Ki-67, and LNM (positive vs. negative) and molecular subtypes (Luminal type vs. non-Luminal type, TN type vs. non-TN type, and HER2-positive type vs. non-HER2-positive type). All data were tested first with the Kolmogorov–Smirnov test for normality analysis. The quantitative diffusion parameters and max diameter of lesions between different subgroups were compared using the Mann–Whitney U test. Spearman correlations were used to characterize the correlations of multi-model-derived histogram metrics with prognostic factors and molecular subtypes. With pathologic results as the gold standard, the receiver operating characteristic (ROC) curve analysis was used to assess the diagnostic efficacy of each parameter or each model, and the area under the ROC curve (AUC) was calculated. Then, the largest AUC of each parameter was selected to establish the IVIM model (D, D*, and f), the DKI model (Kapp and Dapp), and the combined three diffusion models (Mono, IVIM, and DKI) using logistic regression. The AUC comparisons were performed using the DeLong test. The morphologic characteristics were compared using the Chi-squared test. For all tests, the significance was set at p< 0.05/8 = 0.00625 (control for multiple comparisons across five prognostic factors and three molecular subtypes).

Results

Inter-reader agreement

As shown in Table 3, there was an excellent agreement between two readers regarding the representative mean and median metrics of diffusion parameters (range of ICCs, 0.827 – 0.939) and morphological characteristics including the shape, margin, and internal signal (range of ICCs, 0.857 – 0.890).

Table 3

ParametersMetricsICC95% Confidence Interval
Mono-ADCmean0.8930.854 – 0.922
median0.8820.839 – 0.914
IVIM-Dmean0.8610.807 – 0.900
median0.8270.759 – 0.875
IVIM-D*mean0.9390.915 – 0.956
median0.9180.885 – 0.941
IVIM-fmean0.8320.766 – 0.879
median0.8710.820 – 0.907
DKI-Kappmean0.8890.849 – 0.919
median0.9330.908 – 0.951
DKI-Dappmean0.9270.900 – 0.947
median0.9180.888 – 0.940
Morphological characteristics
Shape0.8570.807 – 0.895
Margin0.8670.819 – 0.902
Internal signal0.8900.851 – 0.920

Interobserver agreement for diffusion parameters and morphological characteristics by two readers.

ICC, intraclass correlation coefficient. D* is pseudo-diffusion coefficient.

Relationship of diffusion parameters with prognostic factors and molecular subtypes

The histogram metrics of various diffusion parameters among prognostic factors and molecular subtypes of breast cancer are displayed in Table 4. For Mono-ADC, IVIM-D, and DKI-Dapp, all histogram metrics (mean, median, 25th, and 75th percentile) were significantly lower while DKI-Kapp histogram metrics were significantly higher in ER-positive groups compared to those in ER-negative groups (all p< 0.0625), the same trend was found in PR-positive groups compared with the PR-negative groups (all p< 0.0625).

Table 4

ERPRHER2Ki-67
ParametersHistogram metricsNegativePositiveNegativePositiveNegativePositiveNegativePositive
Mono-ADC25th0.869 ± 0.1470.739 ± 0.1180.853 ± 0.1490.746 ± 0.1250.789 ± 0.1480.809 ± 0.1440.788 ± 0.1350.803 ± 0.156
Median0.963 ± 0.1630.812 ± 0.1220.948 ± 0.1630.817 ± 0.1300.867 ± 0.1600.902 ± 0.1590.868 ± 0.1440.888 ± 0.172
Mean0.979 ± 0.1630.827 ± 0.1260.964 ± 0.1650.832 ± 0.1310.880 ± 0.1610.921 ± 0.1600.881 ± 0.1450.904 ± 0.174
75th1.073 ± 0.1920.906 ± 0.1411.060 ± 0.1920.909 ± 0.1460.960 ± 0.1811.018 ± 0.1870.964 ± 0.1680.991 ± 0.198
Kurtosis3.739 ± 1.8153.650 ± 1.8593.615 ± 1.7833.755 ± 1.8873.903 ± 2.0013.266 ± 1.3703.597 ± 1.7083.764 ± 1.937
Skewness0.612 ± 0.6760.509 ± 0.7170.615 ± 0.6330.500 ± 0.7520.582 ± 0.7410.499 ± 0.6100.453 ± 0.7640.636 ± 0.634
IVIM-D25th0.990 ± 0.2340.857 ± 0.2050.971 ± 0.2310.867 ± 0.2140.914 ± 0.2290.919 ± 0.2280.909 ± 0.2300.921 ± 0.227
Median1.128 ± 0.2260.994 ± 0.1741.113 ± 0.2261.000 ± 0.1781.044 ± 0.2101.072 ± 0.2071.041 ± 0.1891.063 ± 0.224
Mean1.147 ± 0.2231.008 ± 0.1681.131 ± 0.2241.016 ± 0.1721.059 ± 0.2051.091 ± 0.2061.056 ± 0.1811.080 ± 0.224
75th1.280 ± 0.2381.136 ± 0.1851.266 ± 0.2411.141 ± 0.1841.182 ± 0.2141.235 ± 0.2321.183 ± 0.1901.213 ± 0.243
Kurtosis3.440 ± 1.5533.349 ± 1.3663.376 ± 1.5193.401 ± 1.3903.585 ± 1.6133.002 ± 0.9443.416 ± 1.5833.367 ± 1.337
Skewness0.459 ± 0.6280.390 ± 0.6080.481 ± 0.5970.366 ± 0.6310.465 ± 0.6500.331 ± 0.5370.328 ± 0.6820.495 ± 0.550
IVIM-D*25th3.211 ± 4.4322.775 ± 3.9922.884 ± 3.9613.040 ± 4.3933.188 ± 4.4772.530 ± 3.5313.292 ± 4.5432.704 ± 3.875
Median8.321 ± 5.9007.307 ± 5.0407.857 ± 5.5367.663 ± 5.3877.818 ± 5.4627.627 ± 5.4488.050 ± 5.5017.513 ± 5.412
Mean10.022 ± 4.8719.096 ± 4.0409.533 ± 4.4609.478 ± 4.4399.583 ± 4.4359.347 ± 4.4729.741 ± 4.7199.312 ± 4.208
75th14.621 ± 6.30413.424 ± 5.33414.180 ± 6.01713.749 ± 5.61714.020 ± 5.38713.816 ± 6.57613.963 ± 6.02713.941 ± 5.632
Kurtosis5.959 ± 9.4925.912 ± 4.5165.874 ± 5.4175.984 ± 4.5385.868 ± 5.1746.060 ± 4.5275.445 ± 3.7036.327 ± 5.762
Skewness1.282 ± 1.0011.278 ± 0.9071.256 ± 1.0041.301 ± 0.8981.236 ± 0.9611.367 ± 0.9216.327 ± 5.7621.333 ± 1.029
IVIM-f25th0.027 ± 0.0250.025 ± 0.0230.026 ± 0.0250.026 ± 0.0250.027 ± 0.0250.024 ± 0.0210.026 ± 0.0240.026 ± 0.023
Median0.052 ± 0.0290.051 ± 0.0310.052 ± 0.0290.052 ± 0.0320.052 ± 0.0320.052 ± 0.0270.055 ± 0.0330.050 ± 0.029
Mean0.058 ± 0.0280.057 ± 0.0300.057 ± 0.0280.058 ± 0.0310.057 ± 0.0310.059 ± 0.0270.061 ± 0.0320.055 ± 0.027
75th0.082 ± 0.0410.081 ± 0.0460.082 ± 0.0400.082 ± 0.0470.080 ± 0.0460.084 ± 0.0410.087 ± 0.0490.077 ± 0.039
Kurtosis3.695 ± 2.3323.597 ± 2.9703.533 ± 2.2553.735 ± 3.0493.814 ± 3.1683.296 ± 1.3273.388 ± 1.7203.845 ± 3.283
Skewness0.626 ± 0.7090.624 ± 0.7580.591 ± 0.7020.654 ± 0.7650.662 ± 0.7860.551 ± 0.6190.611 ± 0.6430.635 ± 0.804
DKI-Kapp25th0.796 ± 0.1280.901 ± 0.1580.811 ± 0.1260.893 ± 0.1670.872 ± 0.1590.820 ± 0.1390.867 ± 0.1420.844 ± 0.163
Median0.883 ± 0.1250.993 ± 0.1650.895 ± 0.1260.988 ± 0.1710.956 ±0.1680.921 ± 0.1360.966 ± 0.1370.927 ± 0.173
Mean0.878 ± 0.1240.991 ± 0.1480.890 ± 0.1220.986 ± 0.1550.956 ± 0.1510.911 ± 0.1400.957 ± 0.1440.928 ± 0.151
75th0.963 ± 0.1331.081 ± 0.1600.973 ± 0.1281.078 ± 0.1691.037 ± 0.1641.013 ± 0.1491.051 ± 0.1541.011 ± 0.162
Kurtosis3.594 ± 2.0313.653 ± 1.8853.618 ± 2.1533.635 ± 1.7533.751 ± 1.9093.381 ± 2.0093.567 ± 1.7653.675 ± 2.087
Skewness-0.124 ± 0.7520.005 ± 0.812-0.136 ± 0.7420.023 ± 0.8210.021 ± 0.804-0.195 ± 0.738-0.046 ± 0.813-0.056 ± 0.769
DKI-Dapp25th1.179 ± 0.2091.003 ± 0.2081.148 ± 0.2111.014 ± 0.2171.072 ± 0.2381.087 ± 0.1941.079 ± 0.2061.075 ± 0.239
Median1.318 ± 0.2271.124 ± 0.2031.298 ± 0.2271.131 ± 0.2131.195 ± 0.2411.238 ± 0.2191.205 ± 0.2171.213 ± 0.249
Mean1.335 ± 0.2231.146 ± 0.2051.316 ± 0.2251.153 ± 0.2111.215 ± 0.2411.259 ± 0.2111.222 ± 0.2171.236 ± 0.245
75th1.481 ± 0.2651.274 ± 0.2201.467 ± 0.2651.276 ± 0.2241.339 ± 0.2621.417 ± 0.2561.352 ± 0.2451.377 ± 0.275
Kurtosis3.469 ± 1.6413.311 ± 1.1603.314 ± 1.5753.439 ± 1.2103.552 ± 1.5283.041 ± 0.9923.419 ± 1.2433.349 ± 1.505
Skewness0.503 ± 0.6710.471 ± 0.6040.507 ± 0.6020.466 ± 0.6610.540 ± 0.6570.377 ± 0.5710.387 ± 0.6730.562 ± 0.590
LNMLuminalTNHER2
ParametersHistogram metricsNegativePositiveNon-LuminalLuminalNon-TNTNNon-HER2-positiveHER2-positive
Mono-ADC25th0.805 ± 0.1510.780 ± 0.1370.878 ± 0.1420.742 ± 0.1230.772 ± 0.1390.878 ± 0.1430.779 ± 0.1420.877 ± 0.144
Median0.885 ± 0.1650.867 ± 0.1500.975 ± 0.1590.815 ± 0.1260.853 ± 0.1500.967 ± 0.1630.856 ± 0.1520.984 ± 0.155
Mean0.903 ± 0.1670.877 ± 0.1500.991 ± 0.1600.830 ± 0.1280.869 ± 0.1520.982 ± 0.1650.871 ± 0.1541.002 ± 0.155
75th0.986 ± 0.1920.966 ± 0.1721.087 ± 0.1900.908 ± 0.1430.953 ± 0.1741.071 ± 0.1960.952 ± 0.1741.107 ± 0.184
Kurtosis3.726 ± 1.9544.084 ± 2.5613.733 ± 1.8973.660 ± 1.8023.645 ± 1.7803.842 ± 2.0323.709 ± 1.8603.593 ± 1.738
Skewness0.645 ± 0.6350.380 ± 0.7850.616 ± 0.6670.513 ± 0.7200.525 ± 0.7040.657 ± 0.6800.552 ± 0.7090.564 ± 0.661
IVIM-D25th0.920 ± 0.2490.907 ± 0.1820.996 ± 0.2360.862 ± 0.2060.897 ± 0.2100.981 ± 0.2730.895 ± 0.2311.016 ± 0.181
Median1.057 ± 0.2251.046 ± 0.1741.138 ± 0.2270.997 ± 0.1751.034 ± 0.1951.120 ± 0.2401.030 ± 0.2011.160 ± 0.212
Mean1.073 ± 0.2201.063 ± 0.1771.158 ± 0.2241.011 ± 0.1691.049 ± 0.1891.142 ± 0.2421.047 ± 0.1991.178 ± 0.202
75th1.201 ± 0.2321.198 ± 0.1991.294 ± 0.2391.137 ± 0.1841.178 ± 0.2091.276 ± 0.2481.174 ± 0.2111.319 ± 0.231
Kurtosis3.437 ± 1.5403.296 ± 1.2593.368 ± 1.5213.403 ± 1.4053.323 ± 1.3543.617 ± 1.7373.461 ± 1.4973.049 ± 1.145
Skewness0.429 ± 0.6140.403 ± 0.6250.470 ± 0.5900.387 ± 0.6330.394 ± 0.6130.513 ± 0.6280.421 ± 0.6320.416 ± 0.547
IVIM-D*25th3.278 ± 4.3412.371 ± 3.8313.033 ± 4.0762.923 ± 4.2752.779 ± 4.0283.619 ± 4.6893.112 ± 4.3812.282 ± 3.055
Median8.010 ± 5.5517.262 ± 5.2388.117 ± 5.7477.513 ± 5.2467.642 ± 5.2818.140 ± 6.0297.683 ± 5.4508.088 ± 5.488
Mean9.706 ± 4.5169.116 ± 4.2899.726 ± 4.6879.356 ± 4.2789.356 ± 4.29810.015 ± 4.9139.535 ± 4.4479.357 ± 4.454
75th14.146 ± 5.90913.578 ± 5.59914.291 ± 6.28413.736 ± 5.46613.732 ± 5.71014.712 ± 6.09513.994 ± 5.63413.752 ± 6.604
Kurtosis6.106 ± 5.2965.600 ± 4.2435.984 ± 5.7635.899 ± 4.3676.014 ± 4.6685.648 ± 5.9045.831 ± 4.8066.412 ± 5.668
Skewness1.287 ± 0.9981.267 ± 0.8471.270 ± 1.0491.287 ± 0.8781.306 ± 0.9331.191 ± 1.0021.261 ± 0.9101.371 ± 1.118
IVIM-f25th0.027 ± 0.0230.023 ± 0.0250.027 ± 0.0250.025 ± 0.0230.025 ± 0.0230.029 ± 0.0280.026 ± 0.0240.024 ± 0.022
Median0.053 ± 0.0300.049 ± 0.0320.053 ± 0.0300.051 ± 0.0310.051 ± 0.0300.055 ± 0.0340.052 ± 0.0320.050 ± 0.025
Mean0.059 ± 0.0300.055 ± 0.0290.059 ± 0.0290.057 ± 0.0300.056 ± 0.0280.062 ± 0.0350.058 ± 0.0310.054 ± 0.019
75th0.083 ± 0.0450.079 ± 0.0420.084 ± 0.0420.080 ± 0.0450.080 ± 0.0410.088 ± 0.0520.082 ± .0470.078 ± 0.024
Kurtosis3.719 ± 3.0643.490 ± 1.8253.597 ± 2.3953.669 ± 2.8963.570 ± 2.6283.882 ± 2.9653.727 ± 2.9043.231 ± 1.335
Skewness0.618 ± 0.7850.638 ± 0.6320.591 ± 0.7440.647 ± 0.7310.607 ± 0.7130.684 ± 0.8140.657 ± 0.7510.471 ± 0.641
DKI-Kapp25th0.846 ± 0.1690.871 ± 0.1190.794 ± 0.1240.895 ± 0.1590.867 ± 0.1640.813 ± 0.1060.873 ± 0.1510.769 ± 0.141
Median0.936 ± 0.1700.960 ± 0.1330.879 ± 0.1250.988 ± 0.1640.961 ± 0.1640.887 ± 0.1210.961 ± 0.1590.867 ± 0.131
Mean0.932 ± 0.1560.957 ± 0.1320.873 ± 0.1190.986 ± 0.1490.957 ± 0.1550.886 ± 0.1050.959 ± 0.1450.856 ± 0.136
75th1.020 ± 0.1671.047 ± 0.1440.956 ± 0.1271.077 ± 0.1611.049 ± 0.1650.961 ± 0.1171.046 ± 0.1580.950 ± 0.141
Kurtosis3.388 ± 1.4894.084 ± 2.5613.738 ± 2.2983.553 ± 1.6793.629 ± 1.9513.619 ± 1.9513.571 ± 1.7493.891 ± 2.713
Skewness-0.068 ± 0.708-0.020 ± 0.926-0.151 ± 0.7890.015 ± 0.782-0.060 ± 0.801-0.023 ± 0.7440.001 ± 0.769-0.315 ± 0.830
DKI-Dapp25th1.084 ± 0.2351.062 ± 0.2011.185 ±0.1971.005 ± 0.2121.044 ± 0.2171.192 ± 0.2111.056 ± 0.2271.176 ± 0.182
Median1.214 ±0.2431.201 ± 0.2191.335 ±0.2181.126 ± 0.2071.176 ± 0.2251.326 ± 0.2301.180 ± 0.2301.345 ± 0.205
Mean1.238 ± 0.2421.214 ± 0.2131.352 ± 0.2161.148 ± 0.2061.196 ± 0.2211.347 ± 0.2341.202 ± 0.2311.359 ± 0.195
75th1.371 ± 0.2701.354 ± 0.2461.503 ± 0.2611.275 ±0.2201.331 ± 0.2491.483 ± 0.2721.331 ± 0.2521.527 ± 0.250
Kurtosis3.458 ± 1.4833.231 ± 1.1893.414 ± 1.6743.358 ± 1.1743.290 ± 1.1823.694 ± 1.9373.449 ± 1.4203.057 ± 1.026
Skewness0.563 ± 0.6220.337 ± 0.6310.500 ± 0.6360.476 ± 0.6330.444 ± 0.6230.628 ± 0.6540.517 ± 0.6390.336 ± 0.585

Comparisons of mono, IVIM and DKI histogram metrics between different groups with molecular prognostic factors and subtypes.

The data for significance is shown in bold (p < 0.0625). ER, Estrogen receptor; PR, Progesterone receptor; HER2, Human epidermal growth factor receptor 2; LNM, Lymph nodes metastasis; TN, Triple-negative.

Luminal type vs. non-Luminal type revealed that considerable differences originated from histogram metrics (mean, median, 25th, and 75th percentile) of Mono-ADC, IVIM-D, DKI-Kapp, and DKI-Dapp (all p< 0.0625). Significantly higher histogram metrics (mean, median, 25th, and 75th percentile) of Mono-ADC and DKI-Dapp while lower histogram metrics (mean, median, 25th, and 75th percentile) of DKI-Kapp were found in the TN type than in the non-TN type (all p< 0.0625). Additionally, the Mono-ADC (mean, median, 25th, and 75th percentile), IVIM-D(mean, median, and 75th percentile), and DKI-Dapp (mean, median, 25th, and 75th percentile) values were significantly higher and the DKI-Kapp (mean, median, and 25th percentile) values were significantly lower in the HER2-positive type than in the non-HER2-positive type (all p< 0.0625). No statistically significant difference was observed in the negative and positive groups between HER2, Ki-67, and LNM (all p > 0.00625).

Considerable correlations were observed between ER and PR groups as well as Luminal, TN, and HER2-positive types. Diffusion parameters (Mono-ADC, IVIM-D, DKI-Dapp, and DKI-Kapp) largely involved the histogram metrics (mean, median, 25th, and 75th percentile). When including all parameters in three diffusion models, 74 correlations were remarkable (Figure 1).

Figure 1

Among single model parameters, Mono-ADCmedian and Mono-ADCmean generated the best AUC in the positive and negative groups between ER (AUC = 0.766, p< 0.001) and PR (AUC = 0.735, p< 0.001), respectively. Meanwhile, DKI-Kappkurtosis, IVIM-Dskewness, and DKI-Dappskewness generated the best AUC in the positive and negative groups between HER2 (AUC = 0.632, p = 0.010), Ki-67 (AUC = 0.572, p = 0.049), and LNM (AUC = 0.603, p = 0.044), respectively. Regarding the differentiation of Luminal type vs. non-Luminal type, TN type vs. non-TN type, as well as HER2-positive type vs. non-HER2-positive type, the best AUC was derived from the Mono-ADCmean (AUC = 0.785, p< 0.001), Mono-ADC25th (AUC = 0.719, p< 0.001), and Mono-ADC75th (AUC = 0.738, p< 0.001), respectively (Table 5).

Table 5

ADCDD*fKappDappADCDD*fKappDapp
MetricsER-positive vs.ER-negativePR-positive vs. PR-negative
25th0.7520.6870.5130.5230.7370.7520.7050.6480.5110.5080.6920.679
Median0.7660.6930.5510.5120.7530.7390.7300.6620.5100.5120.7070.704
Mean0.7650.7020.5550.5110.7460.7340.7350.6710.5060.5070.7040.699
75th0.7500.6890.5850.5260.7420.7170.7280.6650.5340.5310.6980.701
Kurtosis0.5040.5210.5230.5150.550.5310.4570.5410.5590.5330.5750.596
Skewness0.5180.5370.5350.5140.5380.5070.5160.5490.5560.5270.5580.514
Single model0.7660.7060.7550.7350.6710.719
Three models0.7840.747
MetricsHER2-positive vs. HER2-negativeKi-67-positive vs. Ki-67-negative
25th0.5370.5180.5290.5230.6070.5170.5150.5040.5390.5020.5260.506
Median0.5650.5410.5000.5400.5810.5570.5180.5080.5240.5570.5440.505
Mean0.5840.5520.5030.5590.5810.5660.5240.5160.5220.5630.5350.509
75th0.5930.5780.5090.5750.5470.6020.5250.5190.5020.5580.5420.519
Kurtosis0.6020.6170.5270.5040.6320.5930.5140.5030.5010.5390.5270.538
Skewness0.5540.5790.5450.5370.5720.5880.5690.5720.5080.5080.5080.560
Single model0.6020.6200.6220.5690.5880.611
Three models0.6590.630
MetricsLNM-positive vs. LNM-negativeLuminal vs. Non-Luminal
25th0.5500.5400.5630.5770.5460.5430.7660.6970.5010.5210.7330.746
Median0.5340.5140.5410.5680.5440.5210.78107040.52340.5260.7410.758
Mean0.5460.5070.5430.5480.5540.5350.7850.7180.5280.5270.7380.752
75th0.5290.5100.5360.5400.5570.5210.7720.7060.51510.5450.7370.741
Kurtosis0.5050.5130.5250.5020.5590.5170.5190.5450.5560.5230.5530.565
Skewness0.5920.5010.5140.5140.5250.6030.5030.5290.5610.5290.5620.532
Single model0.5920.5800.6170.7850.7220.773
Three models0.6160.796
MetricsTN vs. Non-TNHER2-positive vs. Non-HER2-positive
25th0.7190.6300.5310.5410.6470.7040.6790.6720.5350.5140.6960.663
Median0.7010.6250.5250.5130.6730.6840.7260.6890.5270.5260.6960.707
Mean0.6970.6360.5320.5290.6620.6780.7370.6980.5080.5110.7010.705
75th0.6770.6200.5490.5260.6850.6480.7380.6990.5260.5430.6720.723
Kurtosis0.5080.5310.5730.5140.5130.5030.5410.5120.5050.5220.5720.612
Skewness0.5310.5690.5770.5260.5190.5550.5320.6240.5090.5790.6250.620
Single model0.7190.6550.7140.7380.7140.738
Three models0.7360.747

AUC of histogram metrics derived from mono, IVIM, and DKI models to predict molecular prognostic factors and subtypes.

The best AUC of every diffusion parameter is shown in bold. ER, Estrogen receptor; PR, Progesterone receptor; HER2, Human epidermal growth factor receptor 2; LNM, Lymph nodes metastasis; TN, Triple-negative; D*, pseudo-diffusion coefficient.

Among single models, the DKI model generated the best AUC in the HER2-positive and HER2-negative groups (AUC = 0.622, p = 0.017), Ki-67-negative and Ki-67-positive groups (AUC = 0.611, p = 0.022), and LNM-positive and LNM-negative groups (AUC = 0.617, p = 0.022). The Mono model generated the best AUC in the ER-positive and ER-negative groups (AUC = 0.766, p< 0.001), PR-positive and PR-negative groups (AUC = 0.735, p< 0.001), Luminal type vs. non-Luminal type (AUC = 0.785, p< 0.001), as well as TN type vs. non-TN type (AUC = 0.719, p< 0.001). Both Mono and DKI models generated the best AUC in the HER2-positive type vs. non-HER2-positive type (AUC = 0.738, p< 0.001) (Table 5).

Regarding the differentiation of positive and negative groups between ER, PR, HER2, and Ki-67, the AUC of the combination of Mono, IVIM, and DKI resulted in the best discriminatory power compared with either model alone. The comparisons of Luminal type versus non-Luminal type, TN type versus non-TN type, and HER2-positive type versus non-HER2-positive type revealed that AUC considerably improved when the combination of Mono, IVIM, and DKI was used compared with either model alone (Table 5).

Comparison of morphologic characteristics between the groups of molecular prognostic factors and subtypes

As summarized in Table 6, the results demonstrated that the margin of breast cancer had significant differences between the ER-positive and ER-negative groups (p = 0.002). No significant differences were observed in residual groups (all p > 0.00625). Two representative cases are shown in Figures 2; 3.

Table 6

Max diameterShapeMarginInternal signal
groupsp-valueRoundOvalIrregularp-valuesmoothspiculatedirregularp-valuehomogeneousheterogeneousrimp-value
ERPositive2.40 ± 0.960.2691137320.0925213150.002*3927140.040
Negative2.59 ± 1.0013183231428193410
PRPositive2.39 ± 0.900.2091134310.253499180.1872134120.097
Negative2.59 ± 1.011321333482591114
HER2Positive2.44 ± 0.900.7881016220.544276150.953201990.840
Negative2.51 ± 0.99143942561128384215
Ki-67Positive2.52 ± 0.940.4091730320.214517210.1863034150.668
Negative2.43 ± 0.977253232102228279
LNMPositive2.48 ± 0.950.992618250.443266170.652202180.994
Negative2.48 ± 0.96123935571126384016
Luminal vs.2.40 ± 0.930.2301216290.1065513180.0104230140.034
Non-Luminal type2.60 ± 0.105161028425163110
TN vs.2.62 ± 0.990.301610160.636163130.33281950.074
Non-TN type2.45 ± 0.95184548671430504219
HER2-positive vs.2.58 ± 1.030.68066130.230121120.06981250.626
Non-HER2-poitive type2.46 ± 0.84184951711631504919

Magnetic resonance imaging morphological characteristics of molecular prognostic factors and subtypes.

*indicates that the correlation is significant at the level of 0.00625 (double-tailed). ER, Estrogen receptor; PR, Progesterone receptor; HER2, Human epidermal growth factor receptor 2; LNM, Lymph nodes metastasis; TN, Triple-negative.

Figure 2

Figure 3

Discussion

In this study, we evaluated the correlation of Mono, IVIM, and DKI parameters with prognostic factors and molecular subtypes of breast cancer using histogram analysis. The Mono and DKI models yielded greater AUC to discriminate prognostic factors and molecular subtypes compared with the IVIM model. The AUC significantly improved when the combination of the three diffusion models was used compared with either model alone except for discriminating LNM-positive and negative. Additionally, the qualitative DWI analysis based on the morphologic characteristics could distinguish between ER-positive and -negative groups.

Previous studies have demonstrated the correlations of diffusion parameters derived from Mono, IVIM, and DKI models with breast cancer prognostic factors (18, 26, 28, 29). ER overexpression could inhibit angiogenesis to reduce perfusion contribution as well as increase cellularity to restrict water diffusion (11, 12, 29). Low perfusion contribution and high cellularity could both result in decreased histogram metrics of Mono-ADC, DKI-Dapp, and IVIM-D and increased histogram metrics of DKI-Kapp in the ER-positive group. The higher DKI-Kappmean in ER-positive tumors was consistent with the result of Yang et al. (16). Due to similarities in hormone receptor effects, PR-positive tumors also have same trend as ER-positive tumors. In our study, the histogram metrics of various diffusion parameters failed to reveal a remarkable difference between the statuses of HER2, Ki-67, and LNM. We speculated that this difference might be related to the inclusion of lesions, the selection of the b values, and the delineation of the ROI.

In terms of molecular subtypes, we analyzed them statistically in the form of binary classification. Previous studies demonstrated that IVIM-D75th was lower in the Luminal type than in the HER2-positive type, and higher IVIM-D and lower IVIM-D* in Luminal A type than in the other subtypes (25, 30). These results were not entirely consistent with our study. Due to the Luminal type being defined as ER and/or PR positive, histogram metrics of Mono-ADC, IVIM-D, DKI-Dapp, and DKI-Kapp can be used to distinguish Luminal type from non-Luminal type, as similar to distinguishing ER and PR status. You et al. revealed that DKI-Kapp entropy value could identify the HER2-positive type and non-HER2-positive type (20). Our study also showed DKI-Kapp histogram metrics, particularly mean, median, and 25th percentile, could differentiate HER2-positive type and non-HER2-positive type. Suo et al. have demonstrated higher Mono-ADC values in the TN subtype than in other subtypes (12); this tendency was also observed in our study with higher Mono-ADC, IVIM-D, and DKI-Dapp histogram metrics in the TN type than those in the non-TN type. The reason may be that the TN type shows a decrease in tumor cellularity with an associated increase in diffusion (31, 32). In summary, various diffusion parameters can quantify tissue cell density, perfusion contribution, and water motion in vivo and may serve as a potential biomarker for differentiating molecular subtypes.

Besides comparing individual parameters, the ROC of various models was also compared. The present study revealed that the AUC of the Mono or DKI model was higher than that of the IVIM model. That is, the Mono or DKI model was superior to the IVIM model in evaluating the correlations of prognostic factors and molecular subtypes of breast cancers. Yang et al. demonstrated that the DKI model was not superior to the Mono model in reflecting the prognostic information of breast cancer (16). Cho et al. demonstrated that the AUC of the IVIM model was higher than that of the Mono model, whereas Feng reported that the AUC of the IVIM model was lower than that of the Mono model (15, 17). The contradictory results might have resulted from the distinct choices of multi-b values and poor repeatability of multi-models. Therefore, the diagnostic value of the three models with various ranges of multi-b values needs further exploration.

Kul et al. reported that the morphology evaluated on DWI provided 83%-84% accuracy in distinguishing between benign and malignant lesions (33). However, Kang et al. reported that the specificity of the high-signal rim in DWI was higher than that of the ADCmean value (80.6% vs. 63.9%) (34). Related studies include one by Cho, who showed that ER-positive tumor tends to show a not-circumscribed margin in mammography compared to ER-negative tumors (35). Different from our present study, the characteristic of smooth margin was more frequently observed in ER-positive tumors. Another study by Yuan et al, reported that the rate of burr sign in ER-positive in DCE-MRI was higher than that in negative groups (36). The trend was also observed in our study but was not significant. Although this study was a preliminary work, the morphologic characteristics assessed using SMS rs-EPI might provide a noninvasive tool for assessing the biologic characteristics and heterogeneity of breast cancers.

The present study had several limitations. First, the patient population was relatively small, and hence a selection bias might exist. Second, 2D ROI was manually drawn on the slice with the largest tumor diameter. This method did not reflect the overall tumor heterogeneity. Third, all MRI data were obtained in a single institution. Further studies are needed to verify the generalizability and reproducibility of our results.

In conclusion, the histogram metrics of multiparametric DWI and morphologic characteristics might be of use in providing prognostic information regarding breast cancer, thus potentially contributing to individualized treatment plans for patients with breast cancer.

Statements

Data availability statement

The datasets presented in this article are not readily available because the datasets generated or analyzed during the study are available from the corresponding author on reasonable request. Requests to access the datasets should be directed to YQ, .

Author contributions

Conceptualization: YQ, FW, and TA. Data curation: YQ, CT, QH, LH. Formal analysis: YQ, QH, FW, and LH. Investigation: YQ, JY, and QH. Methodology: YQ, FW, CT, and MH. Project administration: TA. Software: TY and HZ. Supervision: TA. Visualization: YQ, FW, and MH. Writing - original draft: YQ and FW. Writing - review and editing: TA and TY. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by the Keypoint Research and Development Program of Hubei Province (Grant Number: 2022BCE019).

Conflict of interest

Authors HZ and TY were employed by the company Siemens Healthineers Ltd.

The remaining 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.

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.

References

  • 1

    LuoHBDuMYLiuYYWangMQingHMWenZPet al. Differentiation between luminal a and b molecular subtypes of breast cancer using pharmacokinetic quantitative parameters with histogram and texture features on preoperative dynamic contrast-enhanced magnetic resonance imaging. Acad Radiol (2020) 27(3):e35–44. doi: 10.1016/j.acra.2019.05.002

  • 2

    WeissJMartirosianPTaronJOthmanAEKuestnerTErbMet al. Feasibility of accelerated simultaneous multislice diffusion-weighted MRI of the prostate. J Magn Reson Imaging (2017) 46(5):1507–15. doi: 10.1002/jmri.25665

  • 3

    HuYZhanCYangZZhangXZhangHLiuWet al. Accelerating acquisition of readout-segmented echo planar imaging with a simultaneous multi-slice (SMS) technique for diagnosing breast lesions. Eur Radiol (2021) 31(5):2667–76. doi: 10.1007/s00330-020-07393-5

  • 4

    TangCQinYHuQAiT. Diagnostic value of multi-model high-resolution diffusion-weighted MR imaging in breast lesions: based on simultaneous multi-slice readout-segmented echo-planar imaging. Eur J Radiol (2022) 154:110439. doi: 10.1016/j.ejrad.2022.110439

  • 5

    FilliLGhafoorSKenkelDLiuWWeilandEAndreisekGet al. Simultaneous multi-slice readout-segmented echo planar imaging for accelerated diffusion-weighted imaging of the breast. Eur J Radiol (2016) 85(1):274–8. doi: 10.1016/j.ejrad.2015.10.009

  • 6

    JiangJSZhuLNWuQSunYLiuWXuXQet al. Feasibility study of using simultaneous multi-slice RESOLVE diffusion weighted imaging to assess parotid gland tumors: comparison with conventional RESOLVE diffusion weighted imaging. BMC Med Imaging (2020) 20(1):93. doi: 10.1186/s12880-020-00492-1

  • 7

    McKayJAChurchALRubinNEmoryTHHovenNFKuehn-HajderJEet al. A comparison of methods for high-Spatial-Resolution diffusion-weighted imaging in breast MRI. Radiology (2020) 297(2):304–12. doi: 10.1148/radiol.2020200221

  • 8

    Phi VanVDBeckerASCiritsisAReinerCSBossA. Intravoxel incoherent motion analysis of abdominal organs: application of simultaneous multislice acquisition. Invest Radiol (2018) 53(3):179–85. doi: 10.1097/RLI.0000000000000426

  • 9

    BarthMBreuerFKoopmansPJNorrisDGPoserBA. Simultaneous multislice (SMS) imaging techniques. Magn Reson Med (2016) 75(1):63–81. doi: 10.1002/mrm.25897

  • 10

    TavakoliAAttenbergerUIBudjanJStemmerANickelDKannengiesserSet al. Improved liver diffusion-weighted imaging at 3 T using respiratory triggering in combination with simultaneous multislice acceleration. Invest Radiol (2019) 54(12):744–51. doi: 10.1097/RLI.0000000000000594

  • 11

    MaricJBobanJIvkovic-KapiclTDjilasDVucaj-CirilovicVBogdanovic-StojanovicD. Differentiation of breast lesions and distinguishing their histological subtypes using diffusion-weighted imaging and ADC values. Front Oncol (2020) 10:332. doi: 10.3389/fonc.2020.00332

  • 12

    SuoSZhangDChengFCaoMHuaJLuJet al. Added value of mean and entropy of apparent diffusion coefficient values for evaluating histologic phenotypes of invasive ductal breast cancer with MR imaging. Eur Radiol (2019) 29(3):1425–34. doi: 10.1007/s00330-018-5667-9

  • 13

    HuangYLinYHuWMaCLinWWangZet al. Diffusion kurtosis at 3.0T as an in vivo imaging marker for breast cancer characterization: correlation with prognostic factors. J Magn Reson Imaging (2019) 49(3):845–56. doi: 1002/jmri.26249

  • 14

    HeMSongYLiHLuJLiYDuanSet al. Histogram analysis comparison of monoexponential, advanced diffusion-weighted imaging, and dynamic contrast-enhanced MRI for differentiating borderline from malignant epithelial ovarian tumors. J Magn Reson Imaging (2020) 52(1):257–68. doi: 10.1002/jmri.27037

  • 15

    ChoGYMoyLKimSGBaeteSHMoccaldiMBabbJSet al. Evaluation of breast cancer using intravoxel incoherent motion (IVIM) histogram analysis: comparison with malignant status, histological subtype, and molecular prognostic factors. Eur Radiol (2016) 26(8):2547–58. doi: 10.1007/s00330-015-4087-3

  • 16

    YangZLLiYZhanCAHuYQGuoYHXiaLMet al. Evaluation of suspicious breast lesions with diffusion kurtosis MR imaging and connection with prognostic factors. Eur J Radiol (2021) 145:110014. doi: 10.1016/j.ejrad.2021.110014

  • 17

    FengWGaoYLuXRXuYSGuoZZLeiJQ. Correlation between molecular prognostic factors and magnetic resonance imaging intravoxel incoherent motion histogram parameters in breast cancer. Magn Reson Imaging (2022) 85:262–70. doi: 10.1016/j.mri.2021.10.027

  • 18

    KimEJKimSHParkGEKangBJSongBJKimYJet al. Histogram analysis of apparent diffusion coefficient at 3.0t: correlation with prognostic factors and subtypes of invasive ductal carcinoma. J Magn Reson Imaging (2015) 42(6):1666–78. doi: 10.1002/jmri.24934

  • 19

    ChoiYKimSHYounIKKangBJParkWCLeeAet al. Rim sign and histogram analysis of apparent diffusion coefficient values on diffusion-weighted MRI in triple-negative breast cancer: comparison with ER-positive subtype. PloS One (2017) 12(5):e0177903. doi: 10.1371/journal.pone.0177903

  • 20

    YouCLiJZhiWChenYYangWGuYet al. The volumetric-tumour histogram-based analysis of intravoxel incoherent motion and non-Gaussian diffusion MRI: association with prognostic factors in HER2-positive breast cancer. J Transl Med (2019) 17(1):182. doi: 10.1186/s12967-019-1911-6

  • 21

    KangHSKimJYKimJJKimSLeeNKLeeJWet al. Diffusion kurtosis MR imaging of invasive breast cancer: correlations with prognostic factors and molecular subtypes. J Magn Reson Imaging (2022) 56(1):110–20. doi: 10.1002/jmri.27999

  • 22

    Le BihanD. What can we see with IVIM MRI? Neuroimage (2019) 187:56–67. doi: 10.1016/j.neuroimage.2017.12.062

  • 23

    KishimotoAOKataokaMIimaMHondaMMiyakeKKOhashiAet al. The comparison of high-resolution diffusion weighted imaging (DWI) with high-resolution contrast-enhanced MRI in the evaluation of breast cancers. Magn Reson Imaging (2020) 71:161–9. doi: 10.1016/j.mri.2020.03.007

  • 24

    BarentszMWTavianiVChangJMIkedaDMMiyakeKKBanerjeeSet al. Assessment of tumor morphology on diffusion-weighted (DWI) breast MRI: diagnostic value of reduced field of view DWI. J Magn Reson Imaging (2015) 42(6):1656–65. doi: 10.1002/jmri.24929

  • 25

    LeeYJKimSHKangBJKangYJYooHYooJet al. Intravoxel incoherent motion (IVIM)-derived parameters in diffusion-weighted MRI: associations with prognostic factors in invasive ductal carcinoma. J Magn Reson Imaging (2017) 45(5):1394–406. doi: 10.1002/jmri.25514

  • 26

    SuoSChengFCaoMKangJWangMHuaJet al. Multiparametric diffusion-weighted imaging in breast lesions: association with pathologic diagnosis and prognostic factors. J Magn Reson Imaging (2017) 46(3):740–50. doi: 10.1002/jmri.25612

  • 27

    BoriaFTagliatiCBaldassarreSErcolaniPMarconiESimonettiBFet al. Morphological MR features and quantitative ADC evaluation in invasive breast cancer: correlation with prognostic factors. Clin Imaging (2018) 50:141–6. doi: 10.1016/j.clinimag.2018.02.011

  • 28

    HorvatJVIyerAMorrisEAApteABernard-DavilaBMartinezDFet al. Histogram analysis and visual heterogeneity of diffusion-weighted imaging with apparent diffusion coefficient mapping in the prediction of molecular subtypes of invasive breast cancers. Contrast Media Mol Imaging (2019) 2019:2972189. doi: 10.1155/2019/2972189

  • 29

    HorvatJVBernard-DavilaBHelbichTHZhangMMorrisEAThakurSBet al. Diffusion-weighted imaging (DWI) with apparent diffusion coefficient (ADC) mapping as a quantitative imaging biomarker for prediction of immunohistochemical receptor status, proliferation rate, and molecular subtypes of breast cancer. J Magn Reson Imaging (2019) 50(3):836–46. doi: 10.1002/jmri.26697

  • 30

    KawashimaHMiyatiTOhnoNOhnoMInokuchiMIkedaHet al. Differentiation between luminal-a and luminal-b breast cancer using intravoxel incoherent motion and dynamic contrast-enhanced magnetic resonance imaging. Acad Radiol (2017) 24(12):1575–81. doi: 10.1016/j.acra.2017.06.016

  • 31

    LiQXiaoQYangMChaiQHuangYWuPYet al. Histogram analysis of quantitative parameters from synthetic MRI: correlations with prognostic factors and molecular subtypes in invasive ductal breast cancer. Eur J Radiol (2021) 139:109697. doi: 10.1016/j.ejrad.2021.109697

  • 32

    YoukJHSonEJChungJKimJAKimEK. Triple-negative invasive breast cancer on dynamic contrast-enhanced and diffusion-weighted MR imaging: comparison with other breast cancer subtypes. Eur Radiol (2012) 22(8):1724–34. doi: 10.1007/s00330-012-2425-2

  • 33

    KulSMetinYKulMMetinNEyubogluIOzdemirOet al. Assessment of breast mass morphology with diffusion-weighted MRI: beyond apparent diffusion coefficient. J Magn Reson Imaging (2018) 48(6):1668–77. doi: 10.1002/jmri.26175

  • 34

    KangBJLipsonJAPlaneyKRZackrissonSIkedaDMKaoJet al. Rim sign in breast lesions on diffusion-weighted magnetic resonance imaging: diagnostic accuracy and clinical usefulness. J Magn Reson Imaging (2015) 41(3):616–23. doi: 10.1002/jmri.24617

  • 35

    ChoN. Imaging features of breast cancer molecular subtypes: state of the art. J Pathol Transl Med (2021) 55(1):16–25. doi: 10.4132/jptm.2020.09.03

  • 36

    YuanCJinFGuoXZhaoSLiWGuoHet al. Correlation analysis of breast cancer DWI combined with DCE-MRI imaging features with molecular subtypes and prognostic factors. J Med Syst (2019) 43(4):83. doi: 10.1007/s10916-019-1197-5

Summary

Keywords

diffusion weight imaging, diffusion kurtosis imaging, intravoxel incoherent motion, breast cancer, prognosis, molecular subtypes

Citation

Qin Y, Wu F, Hu Q, He L, Huo M, Tang C, Yi J, Zhang H, Yin T and Ai T (2023) Histogram analysis of multi-model high-resolution diffusion-weighted MRI in breast cancer: correlations with molecular prognostic factors and subtypes. Front. Oncol. 13:1139189. doi: 10.3389/fonc.2023.1139189

Received

06 January 2023

Accepted

17 April 2023

Published

28 April 2023

Volume

13 - 2023

Edited by

Maryam Afzali, University of Leeds, United Kingdom

Reviewed by

Mami Iima, Kyoto University, Japan; Yan Lin, Second Affiliated Hospital of Shantou University Medical College, China

Updates

Copyright

*Correspondence: Tao Ai,

†These authors have contributed equally to this work

Disclaimer

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics