Clinical application of gadoterate meglumine in abbreviated dynamic contrast-enhanced magnetic resonance imaging of the breast with ultrafast imaging: a single-center retrospective study
Introduction
Globally, breast cancer is the most commonly diagnosed cancer in women, and in the United States, it ranks as the second most common cause of cancer-related death among females, following lung cancer. Advances in understanding the diverse biological profiles of breast tumors and their varied clinical outcomes have significantly increased the complexity and sophistication required for accurate diagnosis. Digital mammography, also called full-field digital mammography (FFDM), has been the most widely available imaging modality for screening breast cancer. The major drawback of FFDM is the effect of summation of breast tissue which can obscure some breast lesions. Mammography remains an economically efficient modality and has played a crucial role in facilitating earlier detection of breast malignancies; however, its sensitivity is notably limited, dropping below 50% in women with dense breast tissue and reaching only around 60% in those with non-dense yet structurally complex breasts (1).
Currently, dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) is considered the most sensitive imaging technique for breast cancer detection, providing a significant diagnostic benefit compared to mammography (2). Magnetic resonance imaging (MRI) identifies approximately 14.7 to 15.5 cases of breast cancer per 1,000 examinations in both high-risk and average-risk populations, outperforming mammography and ultrasound (US) in detection yield (2,3). MRI is also substantially more effective than that of digital breast tomosynthesis (DBT) and US, which had a supplement cancer detection rate of 1.2 per 1,000 cases for DBT and 3.5–4.4 per 1,000 cases for US, respectively (4,5).
The interval cancer rate associated with MRI is markedly lower, ranging from 0% to 11% (2,6-9), as compared with other modalities (10,11). Moreover, MRI contributes significantly to improved survival outcomes by enabling the early detection of biologically aggressive and high-grade breast tumors at smaller sizes and earlier stages, often before axillary lymph node involvement, outperforming other imaging modalities (10,11). Given these cumulative advantages, broadening DCE-MRI’s involvement as a stand-alone screening tool is worthwhile.
Gadolinium-based contrast agents (GBCAs) exhibit diverse physicochemical characteristics, among which T1 relaxivity is particularly critical; higher T1 relaxivity enhances signal intensity, thereby improving the detection and delineation of lesions in the brain and breast (12-15). The ionic nature of GBCAs represents another key physicochemical property. These agents may exist in either ionic or non-ionic forms, and this variation in molecular charge can influence tissue contrast uptake, particularly in structures rich in negatively charged elements such as mucopolysaccharides (16,17). It is well recognized that malignant breast lesions are abundant in mucopolysaccharide acids (18-20).
Gadoterate meglumine (Dotarem®, Guerbet, Raleigh, NC, USA), a macrocyclic and ionic GBCA, exhibits lower relaxivity, osmolality, and viscosity relative to gadobenate dimeglumine (21) and is utilized for breast MRI in certain clinical settings (22). Macrocyclic GBCAs feature a pre-formed, rigid ligand architecture that tightly encapsulates the Gd3+ ion, resulting in enhanced thermodynamic and kinetic stability. Compared to their linear counterparts, macrocyclic GBCAs demonstrate superior resistance to dechelation, thereby reducing the likelihood of gadolinium release within the body.
Dotarem is employed to enhance the visibility of vascular structures and soft tissues during MRI examinations. In breast imaging, Dotarem facilitates a more detailed and comprehensive visualization of lesions, aiding radiologists in the early detection and characterization of breast abnormalities by clearer delineation of lesions’ conspicuity (22-25). Intravenous gadoteric acid has been shown in this noninterventional surveillance research to be a safe and efficient contrast agent for use in magnetic resonance (MR) mammography (22). Some investigations using Dotarem for breast MRI imaging have been conducted (22,26-28). However, there has been paucity of quantitative kinetic analysis of signal enhancement due to Dotarem in breast imaging. There have been no prior studies evaluating the effectiveness of Dotarem with high temporal resolution sampling—ultrafast DCE-MRI. The present study provides valuable new information concerning use of Dotarem pharmacokinetics to improve breast cancer. The aim of this study was to evaluate the feasibility and clinical application of Dotarem in DCE-MRI of the breast for breast cancer detection, focusing on quantitative pharmacokinetic analysis with ultrafast DCE-MRI. We present this article in accordance with the STARD reporting checklist (available at https://qims.amegroups.com/article/view/10.21037/qims-2025-653/rc).
Methods
This was a retrospective observational cross-sectional study performed at the Breast Imaging Center and MRI Research Center (MRIRC), Biological Science Division, Department of Radiology, University of Chicago, conducted between November 7th, 2020, and May 30th, 2023. The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments, also in compliance with the Health Insurance Portability and Accountability Act. The study was approved by the Institutional Review Board of the University of Chicago (No. IRB19-2074), and the requirement for written informed consent was waived by institutional policy and the retrospective study design.
Patient recruitment
The patient cohort consisted of women who underwent breast ultrafast DCE-MRI with Dotarem at Department of Radiology, University of Chicago during the study period for various diagnostic indications, including breast pain, breast lumpiness, palpable breast masses, or nipple discharge. Women included in this study were those who underwent ultrafast DCE-MRI of the breast due to either suspicious imaging findings suggestive of malignancy—classified as Breast Imaging Reporting and Data System (BI-RADS) category 4 or 5—or a confirmed diagnosis of breast cancer (BI-RADS category 6). As defined by the BI-RADS lexicon, category 4 lesions present an estimated malignancy probability ranging from 2% to 95%, whereas category 5 lesions are associated with a likelihood of malignancy exceeding 95%. A total of 50 women were initially recruited for this study to evaluate breast lesions and calculate semi-quantitative MRI parameters. One patient was excluded from this study due to a spine condition that impaired her ability to maintain a stable prone position, resulting in suboptimal image quality. This yielded a final cohort of 49 eligible patients. All index lesions were confirmed by core needle biopsy for histopathologic diagnosis.
MRI acquisition
All MRI scans were conducted with patients positioned prone, utilizing a 3.0 Tesla Philips TX platform (Achieva or Ingenia; Philips Healthcare, Best, The Netherlands) equipped with a specialized 16-channel bilateral breast coil (MammoTrak, Philips Healthcare). A standard dose of Dotarem (0.1 mM/kg) was administered intravenously through the antecubital vein using an automated power injector (Spectris Solaris EP, Medrad, PA, USA) at a flow rate of 2.5 mL/s, followed by a 20 mL saline flush. Abbreviated DCE-MRI that includes an ultrafast DCE-MRI sequence in the axial orientation. The acquisition parameters for the ultrafast and high spatial resolution sequence are summarized in Table 1. The ultrafast DCE-MRI protocol included five pre-contrast acquisitions followed by nineteen post-contrast dynamic phases.
Table 1
| Parameters | Ultrafast DCE-MRI protocol | High resolution protocol |
|---|---|---|
| Repetition time/echo time (ms) | 3.2/1.6 | 4.8/2.4 |
| Voxel dimensions (mm3) | 1.5×1.5×4.0 | 0.8×0.8×1.6 |
| SENSE factor (right-left) | 3 | 2.5 |
| SENSE factor (foot-head) | 2 | 2 |
| Partial Fourier acquisition | ky: 0.65; kz: 0.70 | ky: 0.85; kz: 1.00 |
| Temporal resolution (seconds) | 3.5 to 4.6 | 61 to 79.5 |
| Slice count range | 80 to 110 | 220 to 275 |
| Inter-slice gap (mm) | 2.0 | 0.8 |
| Flip angle (degrees) | 10 | 10 |
| Imaging field (AP × RL) (mm2) | (320–360)×(320–360) | (320–360)×(320–360) |
| Fat signal suppression | SPAIR | SPAIR |
AP, antero-posterior; DCE-MRI, dynamic contrast-enhanced magnetic resonance imaging; RL, right-left; SENSE, sensitivity encoding; SPAIR, spectral attenuated inversion recovery.
Image analysis
Image data were processed using a custom-developed MATLAB-based platform in conjunction with the open-source software 3D Slicer (29). Motion correction for ultrafast DCE-MRI was achieved through a non-rigid image registration technique (30). A semi-automatic volumetric segmentation of the lesions was conducted using the ‘Level Tracing’ algorithm in 3D Slicer, applied to subtraction ultrafast DCE-MRI (the last post-contrast image-pre-contrast image) by a radiologist with 5 years of experience of breast imaging. The ‘Level Tracing’ algorithm delineated lesion boundaries by automatically identifying voxels of similar intensity within the same imaging slice. The volumetric tumor region was then captured using a seeded region growing algorithm, which expanded the initial seed area into a three-dimensional (3D) tumor volume (31). A slice of normal background parenchymal enhancement (BPE) at the central part of the contralateral breast was selected randomly without any suspicious lesion or benign pathology and then segmented. Then a region of interest (ROI)-based data analysis was used for the quantitative analysis.
Average signal enhancement in each lesion and the selected normal parenchymal ROIs was calculated. The signal enhancement ratio (SER) over time per ROI, , was fitted with an empirical mathematical model (EMM):
where denotes the average signal intensity over the ROI at time , and is the average signal intensity over ROI from the pre-contrast image, A is the upper limit of initial contrast agent uptake, is the rate of contrast agent uptake, is the bolus arrival time (BAT, s). Then, the maximum enhancement slope is calculated as .
The secondary parameter, area under the enhancement-time curve over 30 seconds (AUC30) for , was calculated by integrating Eq. [1], i.e.,
where t30 =30 s.
The time to reach 90% of the maximum enhancement (T90) was calculated and used as a surrogate metric for time-to-peak enhancement. In additional to EMM, Tofts model was used to obtain quantitative pharmacokinetic parameters (32). Contrast agent concentration curve (in minutes) in lesions and normal parenchyma ROI was calculated using a previously published method (33) and fitted with the Tofts model:
where Ktrans (min−1) refers to the transfer constant, ve denotes the fractional volume of the extravascular extracellular space, and is the arterial input function (AIF) that was traced over an ROI within the central area of the descending thoracic aorta and average from three adjacent slices for individual patient. BAT in the aorta was calculated using Eq. [1]. BATs in lesions and parenchyma were calculated relative to the BAT in the thoracic aorta.
Statistical analysis
All statistical computations were carried out using STATA statistical software version 16 (Stata Corp, College Station, TX, USA). The threshold for statistical significance was set at a P value of less than 0.05. Continuous variables were tested for normality through visual inspection of a histogram and a normal density plot. Data that are normally distributed are reported as mean and standard deviation (SD), while non-normally distributed data are summarized using median and interquartile range (IQR).
Percentage values, mean, and SD of clinicopathological data were calculated for patients’ variables. The mean and SD of all kinetic parameters for each lesion dataset were calculated over ROI in the lesion region. To compare mean values of various kinetic parameters between malignant tumors and background parenchyma, as well as average ratios between malignant lesions and benign lesions or BPE, appropriate unpaired statistical tests were applied based on data distribution and variance. Welch’s t-test was used for normally distributed variables due to its robustness in handling unequal variances and sample sizes. For non-normally distributed data, the Mann-Whitney U test was employed as a nonparametric alternative for comparing independent groups.
Sensitivity and specificity were calculated at the optimal cutoff point, as determined by the maximum Youden’s index, along with their corresponding 95% confidence intervals (CIs).
To address the limitations associated with a small sample size and class imbalance—particularly the underrepresentation of benign lesions—we first calculated the area under the receiver operating characteristic (ROC) curve (AUC) values for each kinetic parameter to assess the discriminating performance of malignant lesions from background parenchymal enhancement (BPE) and benign lesions. Then, we applied nonparametric bootstrap resampling (1,000 iterations) to compute their corresponding 95% CIs. This resampling technique enables the estimation of variability in classifier performance under conditions of data imbalance and enhances the robustness of AUC interpretation by providing a distribution-based measure of uncertainty.
Results
In total, 49 eligible patients presenting with pathologically proven breast lesions, including suspicious abnormalities and known cases of breast malignancy, were enrolled in this study (Figure 1). There are 47 patients with 60 malignant lesions (cancer cohort) and 2 patients with benign lesions (one is chronic mastitis; the other is a phyllodes tumor, borderline subtype). The age of all patients ranged from 34 to 81 years, with the median of all patients being 49 years and the median in the cancer cohort being 52 years. Patient demographic data and lesion characteristics are summarized in Table 2.
Table 2
| Characteristics | Data |
|---|---|
| Number of patients | 49 |
| Age (years) | 49.00±11.89 |
| Age (years) in malignant group | 52.02±11.84 |
| Menopause status | |
| Postmenopausal | 24 (49.0) |
| Perimenopausal | 4 (8.2) |
| Premenopausal | 21 (42.8) |
| Breast density | |
| Dense breast | 30 (61.2) |
| Non-dense breast | 19 (38.8) |
| Qualitative BPE by radiologist | |
| Minimal | 15 |
| Mild | 29 |
| Moderate | 3 |
| Marked | 2 |
| Lesion volume (mm3) | |
| AP dimension | 36.96±25.41 |
| ML dimension | 26.69±18.38 |
| CC dimension | 20.62±14.03 |
| Lateralization | |
| Right | 24 |
| Left | 22 |
| Bilateral | 3 |
| Pathological data among total lesions (n=62) | |
| Malignant pathology of index lesion (n=60) | |
| Subtype of cancer | |
| DCIS | 3 (5.0) |
| IDC | |
| Grade 1 | 5 (8.3) |
| Grade 2 | 29 (48.3) |
| Grade 3 | 17 (28.3) |
| ILC | 2 (3.3) |
| Mixed IDC and ILC | 4 (6.8) |
| Hormone receptor | |
| ER+ and/or PR+ | 52 |
| HER2/neu+ | 1 |
| Triple negative receptor | 8 |
| Lymph node involvement | |
| Yes | 14 (23.3) |
| No | 46 (76.7) |
| Benign pathology result (n=2) | |
| Phyllodes | 1 |
| Chronic mastitis | 1 |
Data are presented as number, mean ± SD, or number (%). AP, antero-posterior; BPE, background parenchymal enhancement; CC, craniocaudal; DCIS, ductal carcinoma in situ; ER, estrogen receptor; HER2, human epidermal growth factor receptor; IDC, invasive ductal carcinoma; ILC, invasive lobular carcinoma; ML, medio-lateral; PR, progesterone receptor; SD, standard deviation.
The average and SD of kinetic parameters of 60 malignant lesions and 2 benign lesions with 47 BPE and average ratio of malignant lesions to BPE were summarized in Table 3. Two cases were excluded from BPE cohort due to silicone implant insertions. All kinetic parameters except ve are significantly different in malignant lesions versus BPE based on the Welch’s t-test and the Mann-Whitney U test (P<0.05).
Table 3
| Parameters | Malignant lesions (n=60) | BPE with benign lesions (n=49) | P value | Average ratio of malignant: BPE |
|---|---|---|---|---|
| A (%) | 87.02±31.04 | 21.39±14.30 | <0.001* | 4.96±3.30 |
| α (%/s) | 9.12±4.73 | 3.78±1.86 | <0.001* | 2.48±1.74 |
| A*α (s−1) | 0.09±0.06 | 0.01±0.01 | <0.001* | 13.13±12.38 |
| AUC30 | 13.96±7.65 | 1.57±1.37 | <0.001* | 12.16±11.11 |
| BAT (s) | 11.34±6.55 | 9.59±6.90 | 0.132 | 2.45±1.90 |
| T90 (s) | 31.37±10.64 | 44.12±13.01 | <0.001* | 0.67±0.25 |
| Ktrans (min−1) | 0.41±1.05 | 0.03±0.06 | <0.001* | 15.58±14.35 |
| ve | 0.22±0.29 | 0.26±0.35 | 0.055 | 4.62±11.02 |
Data are presented as mean ± SD. *, P<0.05. P value of A and α was calculated from Welch’s t-test; P value for the remaining parameters was calculated from the Mann-Whitney U test. A, upper limit of initial contrast agent uptake; α, rate of contrast agent uptake; A*α, maximum enhancement slope; AUC30, area under the enhancement-time curve over 30 seconds; BAT, bolus arrival time; BPE, background parenchymal enhancement; EMM, empirical mathematical model; Ktrans, transfer constant; SD, standard deviation; T90, time to reach 90% of the maximum enhancement; ve, fractional volume of the extravascular extracellular space.
A scatter plot of the average percentage in signal enhancement from the ROIs for all malignant lesions and 47 ROIs for BPE with respective EMM fits is presented in Figure 2. Most malignant lesions were more strongly enhanced than parenchyma, especially within the first minute after contrast injection.
Diagnostic accuracy in terms of sensitivity and specificity with a 95% CI of each kinetic parameter for discriminating breast malignant lesions from background parenchyma and benign lesions is demonstrated in Table 4.
Table 4
| Parameters | Sensitivity (95% CI), % | Specificity (95% CI), % | Optimal cutoff point |
|---|---|---|---|
| A (%) | 96.7 (88.5, 99.6) | 85.7 (72.8, 94.1) | 0.356 |
| α (%/s) | 71.7 (58.6, 82.5) | 89.8 (77.8, 96.6) | 0.065 |
| A*α (s−1) | 95.0 (86.1,99.0) | 91.8 (80.4, 97.7) | 0.017 |
| AUC30 | 96.7 (88.5, 99.6) | 91.8 (80.4, 97.7) | 3.508 |
| BAT (s) | 41.7 (29.1, 55.1) | 75.5 (61.1, 86.7) | 13.796 |
| T90 (s) | 85.7 (72.8, 94.1) | 70.0 (56.8, 81.2) | 34.000 |
| Ktrans (min−1) | 98.3 (91.1, 100) | 77.6 (63.4, 88.2) | 0.029 |
| ve | 95.0 (86.1, 99.0) | 42.9 (28.8, 57.8) | 0.043 |
A, upper limit of initial contrast agent uptake; α, rate of contrast agent uptake; A*α, maximum enhancement slope; AUC, area under receiver operating characteristic curve; AUC30, area under the enhancement-time curve over 30 seconds; BAT, bolus arrival time; BPE, background parenchymal enhancement; CI, confidence interval; Ktrans, transfer constant; T90, time to reach 90% of the maximum enhancement; ve, fractional volume of the extravascular extracellular space.
Figure 3 showed the ROC curves for the evaluated kinetic parameters: A, α, A*α, AUC30, BAT, Ktrans, and ve. The highest discriminative performance of malignant lesions from BPE and benign lesions with AUC values were observed for A*α and AUC30, both at 0.963, with slightly different 95% CI (0.918–0.995 and 0.917–0.996, respectively), indicating comparable but not statistically identical performance. A and Ktrans also showed strong diagnostic accuracy, with AUCs of 0.955 (95% CI: 0.915–0.985) and 0.938 (95% CI: 0.885–0.977), respectively. Rate of contrast agent uptake (α) demonstrated moderate performance (AUC =0.868; 95% CI: 0.793–0.928), while BAT and ve performed poorly (AUCs =0.580 and 0.579). Optimal thresholds for each parameter are indicated on the curves. These results highlight A*α, AUC30, A, and Ktrans as the most effective diagnostic markers in this analysis, as evidenced by their high AUCs and narrow CIs.
An example case of the patients with malignant breast lesions was demonstrated in Figure 4, which shows the differentiation of breast tumors from BPE.
Discussion
Kinetic parameters for detecting breast lesions in MRI images
Assuming that breast cancers are highly permeable, the mean Ktrans (0.41 min−1) is equivalent to blood flow of 40 mLs/100 gms of tissue/minute—which is similar to—but somewhat higher than measurements of breast cancer blood flow using other methods (34-37). Other groups have reported Ktrans either much higher (36) or much lower (34,35,37) than measurements obtained using other methods. In addition, the large range of Ktrans values are consistent with the inter-tumor variation in aggressiveness and angiogenesis.
Numerous prior studies have proposed that quantitative kinetic analysis holds promise for improving the detection of malignant lesions. Xu et al. (36) studied the difference in mass performance in DCE-MRI [gadopentetate dimeglumine (Gd-DTPA); Magnevist] characteristics between low-risk and non-low-risk breast cancer recurrences, revealing that low-risk tumors showed significantly lower Ktrans than non-low-risk breast cancer recurrences. The mean Ktrans of low-risk invasive ductal carcinoma (IDC) tumors is 0.928 (SD =0.630), versus that of non-low-risk tumors, which is 1.275 (SD =0.665, P=0.043).
In 2019, Wu et al. (34) conducted a quantitative evaluation of vascular characteristics obtained from ultrafast DCE-MRI using MultiHance to differentiate malignant from benign breast tumors. Applying the Tofts-Kety model pharmacokinetics model, they reported a mean Ktrans value of 0.073 for malignant lesions, with a range of 0.013 to 0.14.
In 2023, DiCarlo et al. (35) designed a method to select an optimal time point to measure DCE-MRI signal intensity and analyzed the pharmacokinetic parameter Ktrans with the Kety-Tofts model. Their analysis reported a median Ktrans for low-risk IDC of 0.07 with an IQR of 0.04.
There are very few studies measured the contrast pharmacokinetic analysis for breast lesion and BPE. In 2017, Kim et al. (37) assessed separation of benign and malignant breast lesions based on routine diagnostic breast MRI exams and with a single dose of Gd-DTPA (Magnevist) and MRI-guided biopsy when necessary. The extended Tofts model was utilized to derive kinetic parameters, and dynamic enhancement in the background parenchyma was assessed using a fibroglandular tissue mask created through principal component analysis (PCA). This analysis revealed that patients with malignancy had significantly higher median Ktrans values in both the malignant lesion (0.081 min−1) and background parenchyma (0.032 min−1) compared to those without malignancy, who demonstrated median values of 0.056 and 0.017 min−1, respectively.
We used the AIF of individual patients and calculated the concentration as a function of time of the thoracic aorta. Use of a rapidly sampled AIF for each patient corrects of interpatient variability and increases sensitivity to rapid enhancement in cancers. These could be one reason why there has been a very large difference between cancer and parenchyma for Ktrans compared to the previous studies in Kim et al. (37,38).
Moreover, according to the AUC of each kinetic parameter for discriminating malignancies from BPE and benign lesions, it was found that A*α, AUC30, A, and Ktrans demonstrated excellent discriminating performance for cancer detection. Furthermore, α is good level for differentiation between malignancies from BPE and benign lesions. BAT and ve do not provide reliable discrimination with the current data acquisition methods. BAT is noise-sensitive and dependent on temporal resolution. The accuracy of ve may be impacted by the limited contrast washout period, given the postcontrast ultrafast DCE-MRI series spans only 65 to 85 seconds. In the future validation, a unidirectional model may be more appropriate for a truncated ultrafast DCE-MRI series, particularly if additional time points can be sampled within the first 15 second post-contrast, when it is reasonable to assume negligible water exchange from the extravascular extracellular space into the blood plasma compartment (39). Improved methods for data acquisition may result in more accurate measurements of BAT and ve.
Ultrafast and BPE
In this study, the use of Dotarem resulted in a lower mean percentage enhancement of background parenchymal tissue compared to values reported by Pineda et al. (38), who employed gadobenate dimeglumine (MultiHance) as the contrast agent. This comparison was derived from scatterplot evaluations of average signal enhancement in both malignant lesions and BPE, along with their respective fits based on the EMM. Given the observation that godoterate meglumine (Dotarem) had lower BPE signal enhancement, it may be advantageous for identifying suspicious breast lesions during breast MRI screening settings, particularly in patients who are young or in premenopausal status and may have fluctuating BPE during the menstrual cycle (40). For low-grade ductal carcinoma in situ (DCIS), which inherently had little to mild enhancement, there might also be advantages.
Ultrafast usefulness and future application
Ultrafast sampling for DCE-MRI is needed to measure blood flow in tumors accurately, differentiate between more and less aggressive cancers, and enhance the contrast between cancer and benign (34,38). For the average breast cancer in the present cohort, the initial 30% of enhancement (which contains critical diagnostic information) occurs within 5–7 seconds. Our current sampling rate, 3.5 seconds per image, gives us 1–2 data points during the important phase of enhancement. This is helpful, but we hope to increase to even higher temporal resolution to obtain over twice as many data points during the first 5 seconds of enhancement. In addition, the average width of the AIF is about 7 seconds, and our current sampling rate gives us about 2 points to sample most of the AIF. Again, this is useful, but significantly higher temporal resolution would be beneficial.
The combination of early kinetic parameters derived from ultrafast DCE-MRI and morphologic detail obtained from high-spatial-resolution post-contrast imaging offers a complementary approach that may improve the dynamic characterization and diagnostic performance of breast lesion assessment on MRI (38).
The current study has some limitations. First, the dataset in this study is relatively small. A larger sample size should be explored in future studies. Second, not every case had a baseline MRI scan before receiving treatment. In certain cases, patients had already undergone neoadjuvant chemotherapy, which may have impacted the kinetic parameters under analysis. Third, this patient cohort has a limited number of benign lesions. In addition, there is a substantial class imbalance of malignant and benign lesions. This imbalance reflects the clinical nature of our retrospective diagnostic cohort, which included patients referred for MRI due to suspicious findings rather than routine screening. As a result, the robustness of specificity estimates and the generalizability of AUC values are limited. A larger benign cohort would provide better clarity in assessing the discriminating power between malignant and benign lesions. Furthermore, AIFs were extracted on a subject-specific basis from thoracic arteries, whose dimensions were smaller than the 3×3 voxel regions employed to calculate the time course of average signal intensity. The accuracy of AIF estimation may be compromised by partial volume effects and motion artifacts, potentially leading to reduced measurements of both the enhancement rate and peak amplitude. This measurement inaccuracy could lead to overestimated values of Ktrans and ve. Additionally, a larger dataset with a diverse range of pathologic subtypes would facilitate the examination of correlations between lesion characteristics (such as cancer classification or grading) and could offer new insights into kinetic parameters for distinguishing between cancer subtypes.
Conclusions
Based on quantitative pharmacokinetic analysis in ultrafast DCE-MRI, the study demonstrated the viability of utilizing Dotarem as an intravenous GBCA for DCE-MRI of the breast with preserved clinical success in breast cancer detection.
Acknowledgments
The authors thank to Guerbet, USA for support and Brenda Gonzales for assistance as research coordinator.
Footnote
Reporting Checklist: The authors have completed the STARD reporting checklist. Available at https://qims.amegroups.com/article/view/10.21037/qims-2025-653/rc
Data Sharing Statement: Available at https://qims.amegroups.com/article/view/10.21037/qims-2025-653/dss
Funding: This work was supported by
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://qims.amegroups.com/article/view/10.21037/qims-2025-653/coif). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments, also in compliance with the Health Insurance Portability and Accountability Act. The study was approved by the Institutional Review Board of the University of Chicago (No. IRB19-2074), and the requirement for written informed consent was waived by institutional policy and the retrospective study design.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- von Euler-Chelpin M, Lillholm M, Vejborg I, Nielsen M, Lynge E. Sensitivity of screening mammography by density and texture: a cohort study from a population-based screening program in Denmark. Breast Cancer Res 2019;21:111. [Crossref] [PubMed]
- Kuhl CK, Strobel K, Bieling H, Leutner C, Schild HH, Schrading S. Supplemental Breast MR Imaging Screening of Women with Average Risk of Breast Cancer. Radiology 2017;283:361-70. [Crossref] [PubMed]
- Berg WA, Zhang Z, Lehrer D, Jong RA, Pisano ED, Barr RG, Böhm-Vélez M, Mahoney MC, Evans WP 3rd, Larsen LH, Morton MJ, Mendelson EB, Farria DM, Cormack JB, Marques HS, Adams A, Yeh NM, Gabrielli G. Detection of breast cancer with addition of annual screening ultrasound or a single screening MRI to mammography in women with elevated breast cancer risk. JAMA 2012;307:1394-404. [Crossref] [PubMed]
- Friedewald SM, Rafferty EA, Rose SL, Durand MA, Plecha DM, Greenberg JS, Hayes MK, Copit DS, Carlson KL, Cink TM, Barke LD, Greer LN, Miller DP, Conant EF. Breast cancer screening using tomosynthesis in combination with digital mammography. JAMA 2014;311:2499-507. [Crossref] [PubMed]
- Scheel JR, Lee JM, Sprague BL, Lee CI, Lehman CD. Screening ultrasound as an adjunct to mammography in women with mammographically dense breasts. Am J Obstet Gynecol 2015;212:9-17. [Crossref] [PubMed]
- Riedl CC, Luft N, Bernhart C, Weber M, Bernathova M, Tea MK, Rudas M, Singer CF, Helbich TH. Triple-modality screening trial for familial breast cancer underlines the importance of magnetic resonance imaging and questions the role of mammography and ultrasound regardless of patient mutation status, age, and breast density. J Clin Oncol 2015;33:1128-35. [Crossref] [PubMed]
- Kuhl C, Weigel S, Schrading S, Arand B, Bieling H, König R, Tombach B, Leutner C, Rieber-Brambs A, Nordhoff D, Heindel W, Reiser M, Schild HH. Prospective multicenter cohort study to refine management recommendations for women at elevated familial risk of breast cancer: the EVA trial. J Clin Oncol 2010;28:1450-7. [Crossref] [PubMed]
- Vreemann S, Gubern-Mérida A, Schlooz-Vries MS, Bult P, van Gils CH, Hoogerbrugge N, Karssemeijer N, Mann RM. Influence of Risk Category and Screening Round on the Performance of an MR Imaging and Mammography Screening Program in Carriers of the BRCA Mutation and Other Women at Increased Risk. Radiology 2018;286:443-51. [Crossref] [PubMed]
- Luiten JD, Voogd AC, Luiten EJT, Duijm LEM. Trends in incidence and tumour grade in screen-detected ductal carcinoma in situ and invasive breast cancer. Breast Cancer Res Treat 2017;166:307-14. [Crossref] [PubMed]
- Sardanelli F, Podo F. Breast MR imaging in women at high-risk of breast cancer. Is something changing in early breast cancer detection? Eur Radiol 2007;17:873-87. [Crossref] [PubMed]
- Podo F, Santoro F, Di Leo G, Manoukian S, de Giacomi C, Corcione S, Cortesi L, Carbonaro LA, Trimboli RM, Cilotti A, Preda L, Bonanni B, Pensabene M, Martincich L, Savarese A, Contegiacomo A, Sardanelli F. Triple-Negative versus Non-Triple-Negative Breast Cancers in High-Risk Women: Phenotype Features and Survival from the HIBCRIT-1 MRI-Including Screening Study. Clin Cancer Res 2016;22:895-904. [Crossref] [PubMed]
- Attenberger UI, Runge VM, Morelli JN, Williams J, Jackson CB, Michaely HJ. Evaluation of gadobutrol, a macrocyclic, nonionic gadolinium chelate in a brain glioma model: comparison with gadoterate meglumine and gadopentetate dimeglumine at 1.5 T, combined with an assessment of field strength dependence, specifically 1.5 versus 3 T. J Magn Reson Imaging 2010;31:549-55. [Crossref] [PubMed]
- Knopp MV, Bourne MW, Sardanelli F, Wasser MN, Bonomo L, Boetes C, et al. Gadobenate dimeglumine-enhanced MRI of the breast: analysis of dose response and comparison with gadopentetate dimeglumine. AJR Am J Roentgenol 2003;181:663-76. [Crossref] [PubMed]
- Martincich L, Faivre-Pierret M, Zechmann CM, Corcione S, van den Bosch HC, Peng WJ, et al. Multicenter, double-blind, randomized, intraindividual crossover comparison of gadobenate dimeglumine and gadopentetate dimeglumine for Breast MR imaging (DETECT Trial). Radiology 2011;258:396-408. [Crossref] [PubMed]
- Pediconi F, Catalano C, Padula S, Roselli A, Dominelli V, Cagioli S, Kirchin MA, Pirovano G, Passariello R. Contrast-enhanced MR mammography: improved lesion detection and differentiation with gadobenate dimeglumine. AJR Am J Roentgenol 2008;191:1339-46. [Crossref] [PubMed]
- Gillis A, Gray M, Burstein D. Relaxivity and diffusion of gadolinium agents in cartilage. Magn Reson Med 2002;48:1068-71. [Crossref] [PubMed]
- Wiener E, Woertler K, Weirich G, Rummeny EJ, Settles M. Contrast enhanced cartilage imaging: Comparison of ionic and non-ionic contrast agents. Eur J Radiol 2007;63:110-9. [Crossref] [PubMed]
- Marotta M, D'Armiento FP, Martino G, Donato G, Nazzaro A, Vecchione R, Rosati P. Glycosaminoglycans in human breast cancer: morphological and biochemical study. Appl Pathol 1985;3:164-9. [PubMed]
- Olsen EB, Trier K, Eldov K, Ammitzbøll T. Glycosaminoglycans in human breast cancer. Acta Obstet Gynecol Scand 1988;67:539-42. [Crossref] [PubMed]
- Takeuchi J, Sobue M, Sato E, Shamoto M, Miura K. Variation in glycosaminoglycan components of breast tumors. Cancer Res 1976;36:2133-9. [PubMed]
- Kanal E, Maravilla K, Rowley HA. Gadolinium contrast agents for CNS imaging: current concepts and clinical evidence. AJNR Am J Neuroradiol 2014;35:2215-26. [Crossref] [PubMed]
- Seithe T, Braun J, Wolf M, Vahldiek J, Wolny D, Auer J, Pociej J, Heine O, Hamm B, de Bucourt M. Diagnostic efficacy and safety of gadoteric acid MR mammography in 1537 patients. Eur J Radiol 2016;85:2281-7. [Crossref] [PubMed]
- Shin HJ, Chae EY, Choi WJ, Ha SM, Park JY, Shin KC, Cha JH, Kim HH. Diagnostic Performance of Fused Diffusion-Weighted Imaging Using Unenhanced or Postcontrast T1-Weighted MR Imaging in Patients With Breast Cancer. Medicine (Baltimore) 2016;95:e3502. [Crossref] [PubMed]
- Henderson SA, Muhammad Gowdh N, Purdie CA, Jordan LB, Evans A, Brunton T, Thompson AM, Vinnicombe S. Breast cancer: influence of tumour volume estimation method at MRI on prediction of pathological response to neoadjuvant chemotherapy. Br J Radiol 2018;91:20180123. [Crossref] [PubMed]
- Muttalib M, Ibrahem R, Khashan AS, Hajaj M. Prospective MRI assessment for invasive lobular breast cancer. Correlation with tumour size at histopathology and influence on surgical management. Clin Radiol 2014;69:23-8. [Crossref] [PubMed]
- Clauser P, Helbich TH, Kapetas P, Pinker K, Bernathova M, Woitek R, Kaneider A, Baltzer PAT. Breast lesion detection and characterization with contrast-enhanced magnetic resonance imaging: Prospective randomized intraindividual comparison of gadoterate meglumine (0.15 mmol/kg) and gadobenate dimeglumine (0.075 mmol/kg) at 3T. J Magn Reson Imaging 2019;49:1157-65. [Crossref] [PubMed]
- Fallenberg EM, Renz DM, Karle B, Schwenke C, Ingod-Heppner B, Reles A, Engelken FJ, Huppertz A, Hamm B, Taupitz M. Intraindividual, randomized comparison of the macrocyclic contrast agents gadobutrol and gadoterate meglumine in breast magnetic resonance imaging. Eur Radiol 2015;25:837-49. [Crossref] [PubMed]
- Arponen O, Sudah M, Sutela A, Taina M, Masarwah A, Liimatainen T, Vanninen R. Gadoterate meglumine decreases ADC values of breast lesions depending on the b value combination. Sci Rep 2018;8:87. [Crossref] [PubMed]
- Fedorov A, Beichel R, Kalpathy-Cramer J, Finet J, Fillion-Robin JC, Pujol S, Bauer C, Jennings D, Fennessy F, Sonka M, Buatti J, Aylward S, Miller JV, Pieper S, Kikinis R. 3D Slicer as an image computing platform for the Quantitative Imaging Network. Magn Reson Imaging 2012;30:1323-41. [Crossref] [PubMed]
- Vercauteren T, Pennec X, Perchant A, Ayache N. Diffeomorphic demons: efficient non-parametric image registration. Neuroimage 2009;45:S61-72. [Crossref] [PubMed]
- Ren Z, Pineda FD, Howard FM, Hill E, Szasz T, Safi R, Medved M, Nanda R, Yankeelov TE, Abe H, Karczmar GS. Differences Between Ipsilateral and Contralateral Early Parenchymal Enhancement Kinetics Predict Response of Breast Cancer to Neoadjuvant Therapy. Acad Radiol 2022;29:1469-79. [Crossref] [PubMed]
- Tofts PS, Brix G, Buckley DL, Evelhoch JL, Henderson E, Knopp MV, Larsson HB, Lee TY, Mayr NA, Parker GJ, Port RE, Taylor J, Weisskoff RM. Estimating kinetic parameters from dynamic contrast-enhanced T(1)-weighted MRI of a diffusable tracer: standardized quantities and symbols. J Magn Reson Imaging 1999;10:223-32. [Crossref] [PubMed]
- Ren Z, Pineda FD, Howard FM, Fan X, Nanda R, Abe H, Kulkarni K, Karczmar GS. Bilateral asymmetry of quantitative parenchymal kinetics at ultrafast DCE-MRI predict response to neoadjuvant chemotherapy in patients with HER2+ breast cancer. Magn Reson Imaging 2023;104:9-15. [Crossref] [PubMed]
- Wu C, Pineda F, Hormuth DA 2nd, Karczmar GS, Yankeelov TE. Quantitative analysis of vascular properties derived from ultrafast DCE-MRI to discriminate malignant and benign breast tumors. Magn Reson Med 2019;81:2147-60. [Crossref] [PubMed]
- DiCarlo JC, Jarrett AM, Kazerouni AS, Virostko J, Sorace A, Slavkova KP, Woodard S, Avery S, Patt D, Goodgame B, Yankeelov TE. Analysis of simplicial complexes to determine when to sample for quantitative DCE MRI of the breast. Magn Reson Med 2023;89:1134-50. [Crossref] [PubMed]
- Xu T, Zhang L, Xu H, Kang S, Xu Y, Luo X, Hua T, Tang G. Prediction of low-risk breast cancer using quantitative DCE-MRI and its pathological basis. Oncotarget 2017;8:114360-70. [Crossref] [PubMed]
- Kim SG, Freed M, Leite APK, Zhang J, Seuss C, Moy L. Separation of benign and malignant breast lesions using dynamic contrast enhanced MRI in a biopsy cohort. J Magn Reson Imaging 2017;45:1385-93. [Crossref] [PubMed]
- Pineda FD, Medved M, Wang S, Fan X, Schacht DV, Sennett C, Oto A, Newstead GM, Abe H, Karczmar GS. Ultrafast Bilateral DCE-MRI of the Breast with Conventional Fourier Sampling: Preliminary Evaluation of Semi-quantitative Analysis. Acad Radiol 2016;23:1137-44. [Crossref] [PubMed]
- Heye AK, Thrippleton MJ, Armitage PA, Valdés Hernández MDC, Makin SD, Glatz A, Sakka E, Wardlaw JM. Tracer kinetic modelling for DCE-MRI quantification of subtle blood-brain barrier permeability. Neuroimage 2016;125:446-55. [Crossref] [PubMed]
- Jung Y, Jeong SK, Kang DK, Moon Y, Kim TH. Quantitative analysis of background parenchymal enhancement in whole breast on MRI: Influence of menstrual cycle and comparison with a qualitative analysis. Eur J Radiol 2018;103:84-9. [Crossref] [PubMed]

