Abstract

Objectives. To investigate the association between quantitative parameters generated using synthetic magnetic resonance imaging (SyMRI) and diffusion-weighted imaging (DWI) and Ki-67 expression level in patients with invasive ductal breast cancer (IDC). Method. We retrospectively reviewed the records of patients with IDC who underwent SyMRI and DWI before treatment. Precontrast and postcontrast relaxation times (T1, longitudinal; T2, transverse), proton density (PD) parameters, and apparent diffusion coefficient (ADC) values were measured in breast lesions. Univariate and multivariate regression analyses were performed to screen for statistically significant variables to differentiate the high (≥30%) and low (<30%) Ki-67 expression groups. Their performance was evaluated by receiver operating characteristic (ROC) curve analysis. Results. We analyzed 97 patients. Multivariate regression analysis revealed that the high Ki-67 expression group (n = 57) had significantly higher parameters generated using SyMRI (pre-T1, ) and lower ADC values () compared with the low Ki-67 expression group (n = 40). Pre-T1 showed the best diagnostic performance for predicting the Ki-67 expression level in patients with invasive ductal breast cancer (areas under the ROC curve (AUC), 0.711; 95% confidence interval (CI), 0.609–0.813). Conclusions. Pre-T1 could be used to predict the pretreatment Ki-67 expression level in invasive ductal breast cancer.

1. Introduction

Breast cancer is a highly heterogeneous disease and the most frequently diagnosed cancer in almost all regions of the world [1, 2]. Because the various molecular subtypes of breast cancer lead to different clinical outcomes, molecular typing is critical for the accurate diagnosis and treatment of affected patients. According to the St. Gallen 2013 immunohistochemical (IHC) classification, IHC biomarkers, including estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2), and Ki-67 can be used as molecular subtypes in breast cancer [3].

Ki-67 is a nuclear DNA-binding protein highly expressed in the G1, S, and G2 phases of the cell cycle but not in the quiescent G0 phase. This protein is a clinically important proliferation marker used for grading various cancers [4, 5]. Clinical studies demonstrated that patients with a high Ki-67 index showed a favorable response to chemotherapy but had a relatively poor prognosis [4, 6]. Thus, timely identification of this population is important for appropriate treatment planning. Currently, we determine the preoperative Ki-67 proliferation index using IHC, which requires sufficient tumor tissue, typically obtained by core needle biopsy. However, biopsies are invasive and cannot be used to evaluate the status of the entire lesion or follow changes in the tumor microenvironment after neoadjuvant therapy [5, 7]. Furthermore, it is impossible to obtain tumor tissue specimens in many cases.

Magnetic resonance imaging (MRI) is a noninvasive tool that can assess multiple biomarkers in patients with breast cancer [8, 9]. The use of several advanced sequences in probing breast cancer biology, including diffusion-weighted imaging (DWI) and synthetic MRI (SyMRI) is still under investigation.

DWI is a noncontrast imaging technique that characterizes tissues by their random movement of water molecules within them [1012]. The apparent diffusion coefficient (ADC) derived from DWI is a very promising quantitative parameter for the differential diagnosis of breast lesions and monitoring the effect of neoadjuvant treatment [13, 14]. However, while some studies analyzed associations between DWI and histopathological features in breast cancer, including between ADC and the expression of Ki-67, their findings were quite divergent. For example, some authors [5, 15] reported significant correlations between ADC and Ki-67 status, but others [1618] did not identify this relationship.

SyMRI is a compilation pulse sequence that uses a multidynamic and multiecho (MDME) acquisition method. SyMRI provides absolute values for tissue properties, such as T1 and T2 relaxation times and proton density (PD), from a single acquisition with good accuracy and reproducibility, even across instruments from different vendors [1923]. Thus far, some studies have reported on the value of SyMRI in discriminating breast lesions [2427]. However, findings describing the relationships between the parameters generated using SyMRI and Ki-67 status are scarce and inconsistent.

Given the above facts, further assessment of these potential imaging biomarkers is warranted. Therefore, we investigated, for the first time, the ability of quantitative parameters derived from DWI and SyMRI to accurately predict the pretreatment Ki-67 expression level in invasive ductal breast cancer (IDC).

2. Material and Methods

2.1. Patients

Our institutional review board approved this retrospective study and waived the requirement for informed consent. We enrolled consecutive female patients who underwent breast MRI (including SyMRI and DWI) for a breast lesion that was subsequently diagnosed as histopathological IDC after excision surgery or core needle biopsy between January 2020 and August 2022. The exclusion criteria were as follows: (1) neoadjuvant treatment before breast MRI, (2) incomplete scans, (3) nonmass enhancement lesions on dynamic contrast-enhanced MRI (DCE-MRI), and (4) lesions too small to be identified by the SyMRI or DWI sequence. For patients with multiple lesions, only the largest mass on the DCE-MRI was used for analysis. The patient selection is summarized in Figure 1.

2.2. MRI Protocol

All breast MR images were acquired using a 3T MRI scanner (Signa Pioneer, GE Healthcare, Chicago, IL, US) with an 8-channel phased-array breast coil. The patients entered the scanner feet first in the prone position. We used a power injector to administer a body weight-adjusted dose of Gd-DTPA (0.1 mmol/kg) intravenously at a rate of 3.0 ml/s to each patient, followed by 20 ml of saline flush at the same rate.

The scanning protocol consisted of conventional and quantitative scanning. T1- and T2-weighted images and DCE-MRI images were acquired. DCE-MRI was obtained using 3-dimensional (3D) differential subsampling with Cartesian ordering (DISCO) technology. Quantitative MR images included those obtained using SyMRI and DWI. SyMRI was performed before and after enhancement using a 2D fast spin-echo MDME sequence. The scan parameters are listed in Table 1.

2.3. Image Analysis

Two radiologists, each with more than seven years of breast MRI experience, performed a consensus review of the MRI findings. Lesion features, including size, shape, margin, enhancement pattern, and kinetic curves, were evaluated according to the Breast Imaging Reporting and Data System (BI-RADS) MRI lexicon.

Synthetic MRI data were postprocessed using SyMRI 8.0 software (Synthetic MR, Linköping, Sweden). Using DCE-MRI for lesion localization, a region of interest (ROI) was manually drawn to cover the enhancing solid portion of the lesion on approximately the same slice on the synthetic images and ADC maps. For each radiologist, the best quality image among the synthetic images was chosen for analysis, and values for T1, T2 relaxation time, and PD were simultaneously calculated within an ROI. All of the quantitative parameters (T1, T2, PD, and ADC values) were automatically produced in the workstation. The average measurements of the two radiologists were used for analysis. We recorded preparameters (pre-T1, pre-T2, and pre-PD) and postparameters (post-T1, post-T2, and post-PD generated from SyMRI.

2.4. Histopathologic Analysis

Pathologic reports were reviewed to identify lymph node metastasis and histologic grade. The histological grade was assessed using the Elston–Ellis system. ER, PR, and HER2 and the Ki-67 expression level were evaluated by IHC. ER and PR studies were considered positive when at least 1% of the tumor cells showed positive nuclear staining. The HER2 scanning intensity was scored as 0, 1+, 2+, or 3+; we considered 0 and 1+ negative, 2+ equivocal, and 3+ positive. An immunohistochemistry HER2 score of 2+ was further explored by in situ hybridization to determine HER2 gene amplification. A Ki-67 level ≥30% was considered high and a level <30% was considered low [28].

2.5. Statistical Analysis

Statistical analyses were performed using R language, version 4.1.0 (R Foundation for Statistical Computing, Vienna, Austria). The R packages doBy and ggplot2 were applied. All continuous variables are expressed as means ± standard deviation (SD), while categorical variables are shown as totals and proportions. The interobserver consistencies for all quantitative parameters between the two radiologists were evaluated with the intraclass correlation coefficient (ICC) analysis. Agreement was defined as good (ICC > 0.75), moderate (ICC = 0.5–0.75), or poor (ICC < 0.5). Clinicopathological and DCE-MRI features were compared using the Mann–Whitney U and chi-square tests. Quantitative parameters (T1, T2, PD, and ADC values) were performed both univariate and multivariate logistic regression analyses with a variable selection criterion of . We estimated the area under the receiver operating characteristic curve (AUC) to evaluate the predictive ability of the quantitative parameters. For all tests, was considered statistically significant.

3. Results

3.1. Patients

We enrolled 136 patients. After applying our exclusion criteria, 97 patients (age range, 30–71 years; mean age, 49.5 years) were included in the study.

3.2. Clinicopathological and DCE-MRI Features

The Ki-67 proliferation indices obtained by analyses of core needle biopsies or surgically excised specimens ranged from 8% to 90%. We categorized 40 tumors (41.2%) as showing low proliferation (the low proliferation group) and the remaining 57 (58.8%) as demonstrating high proliferation (the high proliferation group) using the 30% cutoff. Our comparison of clinicopathological and DCE-MRI features between the two groups is shown in Table 2.

3.3. Interobserver Agreement on Quantitative Parameters

All ICCs between the two radiologists for the quantitative parameters generated from SyMRI and ADC maps were all greater than 0.75, indicating good agreements (range, 0.782–0.890; ).

3.4. Quantitative Parameters

Histopathological and imaging findings of tumors and their associations with the quantitative parameters generated from SyMRI and DWI are summarized in Table 3. The post-T1 values were significantly lower in the noncircumscribed lesions than those in the circumscribed lesions (). Furthermore, pre-T1, post-T1, and post-T2 in the high proliferation group were significantly higher than those in the low proliferation group (pre-T1, ; post-T1, ; post-T2, ), and ADC values in the high proliferation group were significantly lower than those in the low proliferation group (). Among these, multiple logistic regression analysis showed that pre-T1 (OR = 1.003; ) and ADC (OR = 0.035; ) values were statistically significant parameters in predicting the Ki-67 expression level (Table 4). Representative images from tumors with a low and high Ki-67 proliferation index are presented in Figures 2 and 3, respectively. Table 5 shows the ROC analysis of the quantitative parameters. The area under the curve (AUC) of pre-T1 had the greatest discriminative ability with the highest AUC (AUC = 0.711; 95% confidence interval (CI), 0.609–0.813) (Figure 4).

4. Discussion

We compared multiple quantitative parameters (T1, T2, PD, and ADC values) to determine the expression level of Ki-67 in breast IDC. The results showed that the pre-T1 values were a significant predictor for the pretreatment Ki-67 expression level.

Matsuda et al. [29] reported that only the T1-Gd SD obtained from SyMRI was useful to predict Ki-67 status. However, their study only included patients with ER + breast cancer, and they used 14% as a cutoff value to define their low and high Ki-67 proliferation groups. According to the International Ki-67 in Breast Cancer Working Group (IKWG) Consensus Meeting Guideline (2019) [28], Ki-67 analysis should only be used to guide clinical decisions with 5% or less or 30% or more. Therefore, a reevaluation of the previous studies on associations between imaging parameters and Ki-67 expression was needed. We classified IDC lesions into two groups: low Ki-67 group (Ki-67 < 30%) and high Ki-67 group (Ki-67 ≥ 30%) in this study.

T1 values have been suggested to reflect extracellular expansion and underlying pathophysiological processes, using intrinsic and fundamental tissue property change [30, 31]. A previous research showed that hyperproliferative cancers (i.e., with high Ki-67 expression) might outgrow the oxygen supply of their vascular system, resulting in cell necrosis [32]. Therefore, we infer that there may be different native T1 values between the Ki-67 low and high groups. In the present study, tumors with high Ki-67 expression had higher native T1 values than those with low Ki-67 expression, and pre-T1 values differed significantly between the groups. The ROC analysis showed that a possible use of pre-T1 for discrimination of Ki-67 expression level in IDC has a relatively high AUC (AUC = 0.711). These results suggest the potential of pre-T1 values as a noninvasive surrogate biomarker for estimation of Ki-67 expression level.

Du et al. [33] showed that pre-T2 values had a statistically significant difference between the groups, and they used 14% as a Ki-67 cutoff value. In contrast to their results, we did not find a statistically significant association between pre-T2 values and Ki-67 expression level using 30% as a Ki-67 cutoff value.

Furthermore, our results showed that the post-T1 values were significantly lower in the noncircumscribed lesions than in the circumscribed lesions (). Further studies are required to address the meaning of this result. Nevertheless, we did not find a correlation between the quantitative parameters generated from SyMRI after contrast agent injection and the Ki-67 expression level. This may be because our study did not include all the different molecular subtypes of breast cancer.

Another interesting aspect of this study is the fact that ADC values cannot be used as a surrogate marker for proliferation activity in IDC. Our results showed that ADC values in the high proliferation group were significantly lower than those in the low proliferation group (). In spite of this, the ROC curve was not found to be useful in diagnostics (AUC = 0.683). The results are consistent with those reported previously [15, 34].

Our study had several limitations. First, patients were enrolled from a single institution with a limited number of relevant cases. Therefore, further studies with a large sample size are required to validate our findings. Second, patients with small lesions that could not be accurately located based on SyMRI were excluded from enrollment, possibly causing sampling bias. Third, we did not apply automatic image registration before and after enhancement to correct for the offset position between SyMRI and DWI, which could have reduced the stability of the analyzed data. Further study should be performed to confirm the accuracy of our findings.

5. Conclusions

In summary, synthetic MRI maybe a useful tool that can be used to predict the Ki-67 expression level in patients with IDC.

Data Availability

All data generated or analyzed during this study are included in this article.

Ethical Approval

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 (A research study was carried out in accordance with the Declaration of Helsinki (as revised in 2013).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the Medical Science and Technology Plan Joint-Project of Henan Province (grant number LHGJ20200468).