The value of multi-parametric ultrasonography nomogram in predicting nipple-areolar complex involvement in breast cancer
Introduction
Breast cancer (BC) is the most commonly diagnosed malignancy in women worldwide, posing a significant threat to global public health (1). Increasingly preferred in BC surgery, nipple-sparing mastectomy (NSM) removes breast tissue while preserving the skin and nipple-areolar complex (NAC), offering patients the combination of excellent cosmetic and oncological outcomes (2,3). Previous studies have reported that, despite a clinically normal NAC on physical examination, NAC involvement can be identified in approximately 16.7% of cases by postoperative pathology (4). Consequently, the primary oncological risk associated with NSM is the potential retention of residual tumor cells within the preserved NAC, which may subsequently lead to local recurrence (5). The literature reports a low incidence of local cancer recurrence after NSM, provided that careful patient selection criteria are applied (6). Therefore, precise assessment of NAC involvement in BC is crucial for appropriate patient management.
Breast magnetic resonance imaging (MRI) is routinely employed for the preoperative evaluation of patients being considered for NSM with higher sensitivity than mammography and conventional ultrasound (CUS) (3). MRI features that have been shown to significantly correlate with NAC include tumor size, multifocality/multicentricity, central location, short tumor-to-nipple distance (TND), ductal enhancement extending to the nipple (DEEN), unilateral NAC enhancement, nipple retraction, and thickening of the NAC (7,8). However, MRI is not suitable for patients with contraindications, such as renal impairment, metal implants, and claustrophobia.
CUS, a common imaging modality, plays an important role in the screening and diagnosis of BC due to its non-invasive, radiation-free nature and low cost (9). Strain elastography (SE) and contrast-enhanced ultrasound (CEUS) have emerged as valuable complements to clinical breast ultrasound. SE identifies pathologies by detecting increased stiffness resulting from calcification, fibrous tissue, or collagen components through variations in tissue elasticity (10). Beyond offering a qualitative visual map of tissue hardness, SE provides a quantitative evaluation using strain rate values (11,12). CEUS provides excellent spatial and temporal resolution to display microcirculatory perfusion, which aids in distinguishing mass lesions, staging invasive cancer, detecting tumor recurrence, and assessing response to neoadjuvant chemotherapy (13,14). Previously, we developed a nomogram combining CUS and CEUS for preoperative prediction of NAC involvement in BC, which show a strong predictive performance, with an area under the curve (AUC) of 0.891 and an overall accuracy of 85.5% (15). Although the combination of SE and CEUS has been shown to be highly valuable in BC diagnosis and therapeutic evaluation, its potential for evaluating malignant involvement of the NAC remains largely unexplored.
In this study, we aimed to compare and evaluate the diagnostic value of combining CUS, SE, and CEUS in predicting NAC involvement in BC, and to establish a relatively objective and clinically practical multi-parametric ultrasonography nomogram for surgical planning. We present this article in accordance with the STARD reporting checklist (available at https://qims.amegroups.com/article/view/10.21037/qims-2026-1-0142/rc).
Methods
This prospective study was registered with the Chinese Clinical Trial Registry (registration no. ChiCTR2000034240) and approved by the Ethics Committee of The First Affiliated Hospital of Shenzhen University Health Science Center, Shenzhen Second People’s Hospital (No. 20200601005). Written informed consent was provided by all study participants. The study was conducted in accordance with the Declaration of Helsinki and its subsequent amendments.
Study participants
Consecutive female BC patients confirmed by needle biopsy between January 2020 and December 2024 were prospectively enrolled. No patient was preselected based on NAC-related symptoms or signs. The inclusion criteria were as follows: (I) newly diagnosed primary BC confirmed by needle biopsy; (II) willingness to undergo CUS, SE, and CEUS examination before surgery and no contraindications for CEUS; (III) subareolar intraoperative biopsy and postoperative immunohistochemical staining were used to assess NAC involvement; and (IV) duration of ≤1 month between ultrasound examination and surgery. The exclusion criteria were as follows: (I) a history of previous surgery, neoadjuvant chemotherapy, radiation therapy, or endocrine therapy involving the NAC; (II) incomplete and low-quality CUS, SE, and CEUS data; and (III) incomplete operative and postoperative pathological information. The patient selection process, including inclusion and exclusion criteria, is illustrated in Figure 1.
CUS, SE, and CEUS examinations
All examinations were conducted following standardized protocols by two radiologists, each having over a decade of experience in breast ultrasound imaging and specific expertise in SE and CEUS. We conducted all ultrasound examinations (CUS, SE, and CEUS) using a GE LOGIQ E9 system (GE Healthcare, Milwaukee, WI, USA), which was equipped with L6–15 (6–15 MHz) and L3–9 (3–9 MHz) linear array transducers. The contrast agent used was Sono Vue (Bracco, Milan, Italy).
Patients were placed in a supine position with both arms elevated and the hands placed beneath the head. The examination was performed with the mechanical index set between 0.18 and 0.22, and with the focus, gain, and depth appropriately adjusted. A thorough CUS scan was performed to image the lesion and nipple simultaneously in the radial plane. For color Doppler assessment, minimal compression was applied while performing multi-planar and multi-angle scanning. Both grayscale and color Doppler characteristics were recorded. For patients with multiple lesions, the lesion closest to the NAC was selected for analysis; for those with non-mass lesions, the area displaying the most conspicuous diffuse calcifications or the most marked architectural distortion was deemed representative and selected.
Furthermore, elasticity images were acquired using SE in dual-display mode. Sampling boxes were placed to contain the nipple, lesions, and sufficient surrounding normal parenchyma. With stable pre-compression, the transducer was oriented perpendicularly and then gently oscillated to produce subtle tissue deformation. A relatively stable elastogram was captured when the color frame displayed green and full, indicating a stable pressure curve.
Lastly, following the acquisition of clear CUS and SE images, the dual-screen CEUS mode was activated to capture the enhancement characteristics. A radial imaging plane simultaneously encompassing both the lesion and the nipple was selected as the target section. Standardized settings included a low mechanical index of 0.06 for all participants to avoid microbubble disruption and to maintain stable contrast enhancement for low-velocity blood flow detection, with the focus placed at the deepest margin of the lesion. Minimal transducer pressure was maintained to avoid vascular compression. A 4.8 mL bolus of SonoVue contrast agent was administered intravenously via the antecubital vein, followed by a 5 mL saline flush. Dynamic imaging commenced at the start of injection and was continuously recorded for 120 seconds for subsequent analysis.
Image interpretation
The CUS, SE, and CEUS images were reviewed independently by two senior radiologists who were blinded to all clinical and pathological data. Any discrepancies between their assessments were resolved through consultation with a third senior radiologist to reach a consensus. CUS images were read according to the 5th edition of the American College of Radiology Breast Imaging Reporting and Data System ultrasound lexicon (16). The following CUS characteristics were recorded: tumor size, TND, location, lesion type, shape, orientation, margin, echo pattern, internal homogeneity, microcalcifications, posterior features, areolar skin thickness, nipple inversion, vascularity, and ductal echo extending from the lesion (defined as the presence of abnormal band-shaped hypoechoic ducts between the lesion and the nipple).
SE images were assessed using the 5-point elasticity scoring system established by Itoh et al. (17). According to this scale, a score of 1 was assigned when the entire lesion appeared green; a score of 2 indicated a lesion with mixed blue and green coloration; a score of 3 was given when the lesion displayed a blue center; a score of 4 corresponded to an entirely blue lesion; and a score of 5 was assigned when both the lesion and the surrounding tissue appeared blue. After scoring, the measurement software of strain rate ratio of elastography was used. The region of interest (ROI) was respectively placed in the lesion, surrounding gland tissues, adjacent fat tissue, and nipple. The strain value was recorded three times for each location, yielding the mean elastography strain value of the lesion (Em-lesion), gland (Em-gland), fat (Em-fat), and nipple (Em-nipple). Then, the lesion-to-gland strain ratio (LGR), lesion-to-fat strain ratio (LFR), and lesion-to-nipple strain ratio (LNR) were calculated.
CEUS features were extracted through frame-by-frame playback analysis of the CEUS images. The following CEUS characteristics were recorded: CEUS tumor size, enhancement homogeneity, enhancement margin, enhancement shape, enhancement order, perfusion defects, enhancement scope, perforator vessels, DEEN, and focal nipple enhancement (FNE). DEEN was defined as ductal enhancement extending from the tumor to the nipple. FNE was defined as focal enhancement in the nipple.
Pathologic evaluation
Pathological evaluation of NAC involvement was performed by an experienced pathologist. When the NAC was surgically removed by total mastectomy, skin-sparing mastectomy, or central lumpectomy, the NAC specimens were assessed according to standard protocols. Specifically, the NAC was dissected from the tissue sample, sectioned vertically along the midline, and sliced at 2–3 mm intervals. Diagnosis was based on hematoxylin and eosin staining and immunohistochemical staining. In cases where the NAC was preserved (e.g., NSM or central lumpectomy without imaging suspicion of involvement), intraoperative histological evaluation was performed on frozen sections of the retroareolar margin. Pathological findings of ductal carcinoma in situ, lobular carcinoma, invasive ductal carcinoma, lymphovascular invasion, or Paget’s disease were classified as positive for NAC involvement; benign lesions such as intraductal papilloma, atypical hyperplasia, or mastopathy were deemed NAC without involvement.
Statistical analysis
The data were analyzed using Empower (R) (X&Y Solutions, Inc., Boston, MA, USA), R software (version 4.0, http://www.empowerstats.net/analysis/index4.php), and SPSS 22.0 software (IBM Corp., Armonk, NY, USA). The normality of continuous variables was examined by the Shapiro-Wilk test. Data following a normal distribution are reported as mean ± standard deviation (SD), whereas non-normally distributed data are summarized as median (25th–75th percentiles). Categorical variables are reported as counts and percentages. Univariate analysis between groups were evaluated with independent-samples t-tests or Mann-Whitney U tests for continuous variables, and chi-square or McNemar tests for categorical variables. Least absolute shrinkage and selection operator (LASSO) logistic regression was employed to identify the independent risk factors among CUS, SE, and CEUS features associated with NAC involvement. Receiver operating characteristic (ROC) curves were generated for each predictive model based on the independent risk factors. The sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), accuracy, and AUC were calculated, and the statistical significance of differences in AUC was evaluated using the Delong test. A multi-parametric ultrasonography nomogram integrating CUS, SE, and CEUS features was developed based on the selected predictors. Internal validation was carried out via 500 bootstrap resamples; calibration curve analysis assessed predictive accuracy, and decision curve analysis (DCA) evaluated clinical utility. Statistical significance was defined as a two-sided P value of less than 0.05.
Results
Clinicopathologic information of enrolled patients
A total of 166 lesions (from 164 patients) met the inclusion criteria and were analyzed (Figure 1). All participants were female with the median age of 50 years (range, 28–82 years). According to the pathological results, approximately 19.87% (33/166) of the patients had confirmed NAC with involvement and 80.12% (133/166) patients were NAC with non-involvement. Comparison of clinical characteristics revealed a statistically significant difference in the prevalence of nipple discharge between patients with and without pathological NAC involvement (P<0.05, Table 1). However, no significant differences were observed between the two groups regarding age, menopausal status, palpable mass, or pain (all P>0.05).
Table 1
| Characteristics | Pathological NAC Involvement | P value | |
|---|---|---|---|
| Involvement (N=33) | Non-involvement (N=133) | ||
| Age (years) | 50.6±11.2 | 49.2±11.8 | 0.508 |
| Menopause | |||
| Yes | 16 (48.5) | 62 (47.0) | 0.876 |
| No | 17 (51.5) | 70 (53.0) | |
| Palpable breast mass | 0.674 | ||
| Yes | 28 (90.3) | 113 (87.6) | |
| No | 3 (9.7) | 16 (12.4) | |
| Pain | 0.480 | ||
| Yes | 6 (19.4) | 19 (14.3) | |
| No | 25 (80.6) | 114 (85.7) | |
| Nipple discharge | 0.049 | ||
| Yes | 5 (15.2) | 7 (5.3) | |
| No | 28 (84.8) | 126 (94.7) | |
Data are reported as mean ± standard deviation or n (%). NAC, nipple-areolar complex.
Ultrasonographic features of patients with NAC involvement and non-involvement
Various CUS and SE parameters, as well as CEUS parameters, exhibited statistically significant differences between the patients with NAC involvement and non-involvement groups (Tables 2-4). On CUS, these were a shorter TND (6.5±8.5 vs. 14.0±9.8 mm; P<0.001), thicker areolar skin (3.4±0.7 vs. 3.1±0.5 mm; P=0.003), and higher percentage of nipple retraction (18.2% vs. 3.8%; P=0.003) in the patients with NAC involvement group. SE findings included higher Em-nipple (3.8±1.5 vs. 2.2±1.1; P<0.001) and a lower LNR (1.4±0.8 vs. 2.3±1.2; P<0.001). On CEUS, DEEN (84.8% vs. 16.5%; P<0.001) and FNE (75.8% vs. 12.8%; P<0.001) were the most prominent features, with peripheral claw-like enhancement also being more common (90.9% vs. 71.4%; P=0.020).
Table 2
| Parameter | Pathological NAC involvement | P value | |
|---|---|---|---|
| Involvement (N=33) | Non-involvement (N=133) | ||
| 2-D tumor size (mm) | 25.1±14.0 | 21.9±12.0 | 0.176 |
| 2-D TND (mm) | 6.5±8.5 | 14.0±9.8 | <0.001 |
| Areola skin thickness (mm) | 3.4±0.7 | 3.1±0.5 | 0.003 |
| Lesion type | 0.087 | ||
| Non-mass | 17 (51.5) | 47 (35.3) | |
| Mass | 16 (48.5) | 86 (64.7) | |
| Nipple retraction | 0.003 | ||
| Present | 6 (18.2) | 5 (3.8) | |
| Absent | 27 (81.8) | 128 (96.2) | |
| Shape | 0.884 | ||
| Round/oval | 2 (6.1) | 9 (6.8) | |
| Irregular | 31 (93.9) | 124 (93.2) | |
| Orientation | 0.739 | ||
| Parallel | 27 (81.8) | 112 (84.2) | |
| Not parallel | 6 (18.2) | 21 (15.8) | |
| Margin | 0.321 | ||
| Circumscribed | 29 (87.9) | 107 (80.5) | |
| Non-circumscribed | 4 (12.1) | 26 (19.5) | |
| Echo pattern | 0.481 | ||
| Hypoechoic | 22 (66.7) | 92 (69.2) | |
| Complex cystic and solid | 3 (9.1) | 11 (8.3) | |
| Hyperechoic/isoechoic | 0 (0.0) | 7 (5.3) | |
| Heterogeneous | 8 (24.2) | 23 (17.3) | |
| Internal homogeneity | 0.992 | ||
| Homogeneous | 2 (6.1) | 8 (6.0) | |
| Heterogeneous | 31 (93.9) | 125 (94.0) | |
| Microcalcifications | 0.173 | ||
| Absent | 13 (39.4) | 70 (52.6) | |
| Present | 20 (60.6) | 63 (47.4) | |
| Posterior features | 0.794 | ||
| Enhancement | 9 (27.3) | 48 (36.1) | |
| Shadowing | 9 (27.3) | 29 (21.8) | |
| Combined pattern | 3 (9.1) | 12 (9.0) | |
| No posterior acoustic features | 12 (36.4) | 44 (33.1) | |
| Echogenic halo | 0.321 | ||
| Present | 11 (33.3) | 33 (24.8) | |
| Absent | 22 (66.7) | 100 (75.2) | |
| Ductal echo extending from lesion | <0.001 | ||
| Present | 24 (72.7) | 35 (26.3) | |
| Absent | 9 (27.3) | 98 (73.7) | |
| Adler grade | 0.568 | ||
| 0 | 4 (12.1) | 16 (12.0) | |
| 1 | 15 (45.5) | 74 (55.6) | |
| 2 | 11 (33.3) | 29 (21.8) | |
| 3 | 3 (9.1) | 14 (10.5) | |
Data are reported as mean ± standard deviation or n (%). CUS, conventional ultrasound; NAC, nipple-areolar complex; TND, tumor-nipple distance.
Table 3
| Parameter | Pathological NAC involvement | P value | |
|---|---|---|---|
| Involvement (N=33) | Non-involvement (N=133) | ||
| Tsukuba score | 0.641 | ||
| 2 | 1 (3.0) | 6 (4.5) | |
| 3 | 4 (12.1) | 10 (7.5) | |
| 4 | 7 (21.2) | 40 (30.1) | |
| 5 | 21 (63.6) | 77 (57.9) | |
| Em-lesion | 4.5±1.1 | 4.3±0.9 | 0.331 |
| Em-gland | 1.0±0.3 | 1.0±0.4 | 0.994 |
| Em-fat continuous | 0.8±0.4 | 0.9±0.7 | 0.633 |
| Em-nipple | 3.8±1.5 | 2.2±1.1 | <0.001 |
| LGR | 6.8±3.4 | 6.2±2.9 | 0.307 |
| LFR | 6.8±3.4 | 6.2±2.9 | 0.307 |
| LNR | 1.4±0.8 | 2.3±1.2 | <0.001 |
Data are reported as mean ± standard deviation or n (%). Em, the mean elastography strain value; LFR, lesion-to-fat strain ratio; LGR, lesion-to-gland strain ratio; LNR, lesion-to-nipple strain ratio; NAC, nipple-areolar complex; SE, strain elastography.
Table 4
| Parameter | Pathological NAC involvement | P value | |
|---|---|---|---|
| Involvement (N=33) | Non-involvement (N=133) | ||
| CEUS tumor size, mm | 32.1±13.2 | 28.4±12.3 | 0.122 |
| Enhancement direction | 0.790 | ||
| Centripetal | 19 (57.6) | 73 (54.9) | |
| Centrifugal | 3 (9.1) | 18 (13.5) | |
| Diffuse | 11 (33.3) | 42 (31.6) | |
| Enhancement intensity | 0.396 | ||
| Hyperenhancement | 32 (97.0) | 125 (94.0) | |
| Hypoenhancement | 0 (0.0) | 6 (4.5) | |
| Isoenhancement | 1 (3.0) | 2 (1.5) | |
| Enhancement mode | 0.795 | ||
| Homogeneous | 1 (3.0) | 3 (2.3) | |
| Heterogeneous | 32 (97.0) | 130 (97.7) | |
| Enhancement shape | 0.879 | ||
| Regular | 3 (9.1) | 11 (8.3) | |
| Irregular | 30 (90.9) | 122 (91.7) | |
| Enhancement boundary | 0.134 | ||
| Clear | 3 (9.1) | 27 (20.3) | |
| Not clear | 30 (90.9) | 106 (79.7) | |
| Enhancement scope | 0.769 | ||
| Enlarged | 28 (84.8) | 110 (82.7) | |
| Not enlarged | 5 (15.2) | 23 (17.3) | |
| Perfusion defect | 0.441 | ||
| Present | 7 (21.2) | 37 (27.8) | |
| Absent | 26 (78.8) | 96 (72.2) | |
| Peripheral claw-like enhancement | 0.020 | ||
| Present | 30 (90.9) | 95 (71.4) | |
| Absent | 3 (9.1) | 38 (28.6) | |
| DEEN | <0.001 | ||
| Present | 28 (84.8) | 22 (16.5) | |
| Absent | 5 (15.2) | 111 (83.5) | |
| FNE | <0.001 | ||
| Present | 25 (75.8) | 17 (12.8) | |
| Absent | 8 (24.2) | 116 (87.2) | |
Data are reported as mean ± standard deviation or n (%). CEUS, contrast-enhanced ultrasound; DEEN, ductal enhancement extending to the nipple; FNE, focal nipple enhancement; NAC, nipple-areolar complex.
Selection of predictive features on multi-parametric ultrasonography
LASSO regression analysis selected four key predictors of NAC involvement from 35 initial sonographic features (Figure 2): ductal echo extending from the lesion on CUS, the Em-nipple value on SE, and both DEEN and FNE on CEUS. Figure 3 illustrates the representative sonographic features that were selected.
Predictive efficacy and comparison of different models
Diagnostic performance for predicting NAC involvement improved with the sequential integration of sonographic modalities (Figure 4). As detailed in Table 5, the CUS + SE model demonstrated improved specificity (0.820) and a higher AUC (0.859) than CUS alone (AUC: 0.732); the combined CUS + SE + CEUS model yielded the highest diagnostic efficacy, achieving excellent sensitivity (0.909), a very high NPV (0.973), and an AUC of 0.917, with all pairwise AUC comparisons reaching statistical significance (P<0.05).
Table 5
| Model | Specificity | Sensitivity | PPV | NPV | Accuracy | AUC | P† |
|---|---|---|---|---|---|---|---|
| CUS | 0.737 | 0.727 | 0.407 | 0.916 | 0.735 | 0.732 | 0.001a |
| CUS + SE | 0.820 | 0.758 | 0.510 | 0.932 | 0.807 | 0.859 | 0.009b |
| CUS + SE + CEUS | 0.820 | 0.909 | 0.556 | 0.973 | 0.837 | 0.917 | <0.001c |
†, the DeLong test was adopted for pairwise comparison of AUCs; a, CUS vs. CUS + SE; b, CUS + SE vs. CUS + SE + CEUS; c, CUS vs. CUS + SE + CEUS. AUC, area under the curve; CEUS, contrast-enhanced ultrasound; CUS, conventional ultrasound; NAC, nipple–areolar complex; NPV, negative predictive value; PPV, positive predictive value; SE, strain elastography.
Development, performance, and validation of the nomogram
The multi-parametric ultrasonography nomogram was constructed based on the four predictors of NAC involvement identified through LASSO regression (Figure 5). The nomogram was verified internally using 500 bootstrap samples to reduce overfitting bias. The internal validation ROC curve (Figure 6A) depicted the bootstrap 95% confidence interval (CI) for the nomogram with the AUC of 0.916. The calibration plot (Figure 6B) illustrated a good fit between the model’s predictions and the actual probabilities. The model’s decision curve (Figure 6C) charted above the “treat all” and “treat none” reference lines, indicating its robust clinical utility for predicting NAC involvement and supporting its value in preoperative decision-making.
Discussion
Driven by a heightened focus on aesthetic outcomes and quality of life post-treatment, NSM has seen broadened adoption among BC patients opting for mastectomy with reconstruction (2). Nonetheless, the oncological safety of NSM is a subject of ongoing concern, as the preserved NAC may serve as a potential reservoir for occult tumor cells (18). Therefore, prudent patient selection is imperative to ensure the oncological safety of NSM. In this study, we developed a multi-parametric ultrasonography nomogram integrating CUS, SE, and CEUS for NSM patient selection based on the ductal echo extending from the lesion, Em-nipple, DEEN, and FNE. With an AUC of 0.917, sensitivity of 90.9%, and NPV of 97.3%, this comprehensive nomogram provides reliable preoperative assessment to guide surgical planning. In our study, pathological involvement of the NAC was confirmed in 19.9% (33/166) of the original cohort. This incidence is consistent with rates reported in the literature, which typically range from 8% to 27.7% (19-21). Our analysis revealed that patients with nipple discharge were more prone to NAC involvement—a finding that aligns with the work of Billar et al., in which abnormal nipple signs or symptoms showed 61% sensitivity, 86% specificity, 45% PPV, and 92% NPV for diagnosing nipple pathology (22). Nipple discharge may result from the involvement of NAC, which occurs as cancer cells extend along the mammary ductal system toward the nipple (23).
Initial univariate analysis revealed a significant association between multiple ultrasound factors and NAC involvement, which included 2-D TND, areolar skin thickness, nipple retraction, ductal echo extending from the lesion, Em-nipple, LNR, peripheral claw-like enhancement, DEEN, and FNE. MRI features significantly correlated with NAC in previous studies (8,24-28), such as TND, nipple retraction, and thickening of the NAC, DEEN, and abnormal nipple enhancement, which are largely in line with the findings in our study. Owing to the disparity between the number of candidate ultrasound predictors and the available sample size, a penalized LASSO logistic regression model was utilized to select an optimal subset of variables most strongly linked to NAC involvement. Following feature selection, four predictors were retained: ductal echo extending from the lesion, the Em-nipple value, DEEN, and FNE. In our previous study, ductal echo extending from the lesion was classified as absent, segmental, or full course (15). As the sample size increased in this study, the frequent presentation of full course of ductal echo extending from the lesion in the NAC involvement group prompted us to dichotomize this feature into a binary variable (present or absent) for simpler recognition. Ductal echo extending from the lesion proximally along the duct may correlate with tumor invasion toward the NAC region. Our results revealed that DEEN was present in 30.1% of cases, with NAC involvement identified in 56.7% of these. A similar proportion (32.1%) of cases showing continuous enhancement from the lesion to the nipple was reported by Sakamoto et al. (26), and NAC involvement was subsequently confirmed in 62% of these cases. As for FNE, previous MRI studies have linked asymmetric nipple enhancement with occult NAC involvement. Although we did not perform contralateral CEUS evaluation of the NAC, we interpret FNE as a manifestation of asymmetric enhancement, a concept supported by its significantly higher incidence in patients with NAC involvement (75.8% vs. 12.8%). The presence of DEEN and FNE is likely due to malignant invasion, which leads to increased vessel density and permeability, thereby manifesting as hyperenhancement on CEUS. To our knowledge, this is the first study to utilize SE for nipple assessment to predict NAC involvement. Our results demonstrated that both Em-lesion and LNR differed significantly between BC patients with and without NAC involvement. Evidence indicates that tumor stiffness quantified via elastography is significantly associated with known indicators of poor prognosis, such as lymph node metastasis, lymphovascular invasion, and immunohistochemical biomarkers (12,29,30).
In our study, the combined CUS + SE + CEUS model yielded the highest diagnostic efficacy (AUC =0.917), which was superior to CUS (AUC=0.732) and CUS + SE (AUC =0.859) models. The value of this combined approach has also been demonstrated in other domains of breast research, including but not limited to differentiating benign from malignant lesions, predicting response to neoadjuvant therapy, and assessing distant metastasis (30-32). Niu et al. developed a nomogram based on a multi-parametric ultrasound approach (incorporating gray-scale, Doppler, shear-wave elastography, and CEUS) to predict ductal carcinoma in situ invasion, which achieved an AUC of 0.903 with 89.7% sensitivity, 73.8% specificity, and 84.3% overall accuracy, outperforming CUS, CUS + SWE, and CUS + CEUS (33). Another study by Liu et al. developed a nomogram integrating SE, CEUS features, and Ki-67 index, which demonstrated excellent performance (AUC =0.986) in predicting pathological complete response to neoadjuvant chemotherapy in BC patients (34). Therefore, developing comprehensive ultrasound models represents a clear and promising direction for future research in breast imaging and personalized treatment strategy formulation.
In terms of the predictive model in our study, the proposed multi-parametric ultrasonography nomogram based on the features selected by LASSO logistic regression demonstrated good discriminatory capability for predicting NAC involvement. Its performance was evidenced by an AUC of 0.916 in the primary analysis and further supported by an internally validated AUC of 0.917. Notwithstanding the existence of several high-performing preoperative predictive models (AUC =0.792–0.984) (8,35,36), most are limited to predictors from breast dynamic contrast-enhanced MRI, contrast-enhanced cone-beam breast computed tomography, clinical data, and their feature selection entails conventional univariate or multivariate regression. The nomogram was developed by integrating and linearly weighting each diagnostic factor to generate an individualized probability. This tool is designed to enable physicians to assess NAC involvement risk rapidly and accurately, thereby facilitating the identification of suitable NSM candidates. Owing to its parsimonious nature, the multi-parametric ultrasonography nomogram offers a practical and efficient diagnostic tool, particularly valuable for radiologists or clinicians with limited experience.
This study has several limitations. Its single-center design and lack of external validation may affect the generalizability of our findings. Furthermore, the limited sample size restricted more in-depth subgroup analyses. Finally, our model did not incorporate quantitative CEUS parameters, which might offer additional diagnostic insights. Based on these limitations, future research should be directed toward multicenter, large-scale prospective studies to externally validate and refine our model. Expanding the cohort will also enable robust subgroup analyses based on lesion characteristics and molecular subtypes. Moreover, integrating quantitative CEUS perfusion parameters represents a promising avenue to further enhance the model’s objectivity and predictive power.
Conclusions
The multi-parametric ultrasonography nomogram integrating CUS, SE, and CEUS demonstrates excellent diagnostic performance for predicting NAC involvement in BC patients and is particularly valuable and convenient for identifying appropriate candidates for NSM, potentially preserving cosmetic outcomes without compromising oncologic safety.
Acknowledgments
None.
Footnote
Reporting Checklist: The authors have completed the STARD reporting checklist. Available at https://qims.amegroups.com/article/view/10.21037/qims-2026-1-0142/rc
Data Sharing Statement: Available at https://qims.amegroups.com/article/view/10.21037/qims-2026-1-0142/dss
Funding: This study 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-2026-1-0142/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. This prospective study was approved by the Ethics Committee of The First Affiliated Hospital of Shenzhen University Health Science Center, Shenzhen Second People’s Hospital (No. 20200601005), and informed consent was obtained from all patients prior to the procedure.
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
- Bray F, Laversanne M, Sung H, Ferlay J, Siegel RL, Soerjomataram I, Jemal A. Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin 2024;74:229-63. [Crossref] [PubMed]
- Rancati AO, Nahabedian MY, Angrigiani C, Irigo M, Dorr J, Acquaviva J, Rancati A. Revascularization of the Nipple-Areola Complex following Nipple-Sparing Mastectomy. Plast Reconstr Surg 2023;151:254-62. [Crossref] [PubMed]
- Lee SC, Mendez-Broomberg K, Eacobacci K, Vincoff NS, Gupta E, McElligott SE. Nipple-sparing Mastectomy: What the Radiologist Should Know. Radiographics 2022;42:321-39. [Crossref] [PubMed]
- Blanckaert M, Vranckx J. Oncological safety of therapeutic 'nipple-sparing mastectomy' followed by reconstruction: a systematic review. Acta Chir Belg 2021;121:155-63. [Crossref] [PubMed]
- Galimberti V, Vicini E, Corso G, Morigi C, Fontana S, Sacchini V, Veronesi P. Nipple-sparing and skin-sparing mastectomy: Review of aims, oncological safety and contraindications. Breast 2017;34:S82-4. [Crossref] [PubMed]
- Clarijs ME, Peeters NJMCV, van Dongen SAF, Koppert LB, Pusic AL, Mureau MAM, Rijken BFM. Quality of Life and Complications after Nipple- versus Skin-Sparing Mastectomy followed by Immediate Breast Reconstruction: A Systematic Review and Meta-Analysis. Plast Reconstr Surg 2023;152:12e-24e. [Crossref] [PubMed]
- Byon JH, Hwang S, Choi H, Choi EJ. Diagnostic Accuracy of Magnetic Resonance Imaging Features and Tumor-to-Nipple Distance for the Nipple-Areolar Complex Involvement of Breast Cancer: A Systematic Review and Meta-Analysis. Korean J Radiol 2023;24:739-51. [Crossref] [PubMed]
- Seki H, Sakurai T, Ishiguro Y, Kanno M, Ikebata A, Katsuki Y, Kaburagi T, Okumura T, Karahashi T, Nakajima K, Shimizu K. A novel MRI-based predictive index can identify patients suitable for preservation of the nipple-areola complex in breast reconstructive surgery. Eur J Surg Oncol 2021;47:225-31. [Crossref] [PubMed]
- Berg WA, Bandos AI, Mendelson EB, Lehrer D, Jong RA, Pisano ED. Ultrasound as the Primary Screening Test for Breast Cancer: Analysis From ACRIN 6666. J Natl Cancer Inst 2016;108:djv367. [Crossref] [PubMed]
- Sigrist RMS, Liau J, Kaffas AE, Chammas MC, Willmann JK. Ultrasound Elastography: Review of Techniques and Clinical Applications. Theranostics 2017;7:1303-29. [Crossref] [PubMed]
- Sinha D, Sharma S, Kundaragi NG, Kale SK. Added value of strain elastography in the characterisation of breast lesions: A prospective study. Ultrasound 2020;28:164-73. [Crossref] [PubMed]
- Mutala TM, Mwango GN, Aywak A, Cioni D, Neri E. Determining the elastography strain ratio cut off value for differentiating benign from malignant breast lesions: systematic review and meta-analysis. Cancer Imaging 2022;22:12. [Crossref] [PubMed]
- Hu W, Dong Y, Zhang X, Zhang H, Li F, Bai M. The clinical value of Arrival-time Parametric Imaging using contrast-enhanced ultrasonography in differentiating benign and malignant breast lesions. Clin Hemorheol Microcirc 2020;75:369-82. [Crossref] [PubMed]
- Balleyguier C, Opolon P, Mathieu MC, Athanasiou A, Garbay JR, Delaloge S, Dromain C. New potential and applications of contrast-enhanced ultrasound of the breast: Own investigations and review of the literature. Eur J Radiol 2009;69:14-23. [Crossref] [PubMed]
- Fang J, Deng Q, Zhang J, Ma Y, Jin C, Lu J, Hao Y, Ma Y, Liu W, Li Z, Liu G, Mo Y, Xiao Y, Zheng C, Yang Y, Wu T, Zhao C, Zhou X, Zhou P. A Nomogram Based on Conventional and Contrast-Enhanced Ultrasound for Pre-operative Prediction of Nipple-Areola Complex Involvement in Breast Cancer: A Prospective Study. Ultrasound Med Biol 2023;49:2565-72. [Crossref] [PubMed]
- Spak DA, Plaxco JS, Santiago L, Dryden MJ, Dogan BE. BI-RADS® fifth edition: A summary of changes. Diagn Interv Imaging 2017;98:179-90.
- Itoh A, Ueno E, Tohno E, Kamma H, Takahashi H, Shiina T, Yamakawa M, Matsumura T. Breast disease: clinical application of US elastography for diagnosis. Radiology 2006;239:341-50. [Crossref] [PubMed]
- Wu ZY, Kim HJ, Lee JW, Chung IY, Kim JS, Lee SB, Son BH, Eom JS, Kim SB, Gong GY, Kim HH, Ahn SH, Ko B. Breast Cancer Recurrence in the Nipple-Areola Complex After Nipple-Sparing Mastectomy With Immediate Breast Reconstruction for Invasive Breast Cancer. JAMA Surg 2019;154:1030-7. [Crossref] [PubMed]
- Koh J, Park AY, Ko KH, Jung HK. MRI diagnostic features for predicting nipple-areolar-complex involvement in breast cancer. Eur J Radiol 2020;122:108754. [Crossref] [PubMed]
- Piato JR, de Andrade RD, Chala LF, de Barros N, Mano MS, Melitto AS, Goncalves R, Soares JM Junior, Baracat EC, Filassi JR. MRI to Predict Nipple Involvement in Breast Cancer Patients. AJR Am J Roentgenol 2016;206:1124-30. [Crossref] [PubMed]
- Wang J, Xiao X, Wang J, Iqbal N, Baxter L, Skinner KA, Hicks DG, Hajdu SI, Tang P. Predictors of nipple-areolar complex involvement by breast carcinoma: histopathologic analysis of 787 consecutive therapeutic mastectomy specimens. Ann Surg Oncol 2012;19:1174-80. [Crossref] [PubMed]
- Billar JA, Dueck AC, Gray RJ, Wasif N, Pockaj BA. Preoperative predictors of nipple-areola complex involvement for patients undergoing mastectomy for breast cancer. Ann Surg Oncol 2011;18:3123-8. [Crossref] [PubMed]
- Morimoto T, Komaki K, Inui K, Umemoto A, Yamamoto H, Harada K, Inoue K. Involvement of nipple and areola in early breast cancer. Cancer 1985;55:2459-63. [Crossref] [PubMed]
- Seki H, Sakurai T, Mizuno S, Tokuda T, Kaburagi T, Seki M, Karahashi T, Nakajima K, Shimizu K, Jinno H. A novel nipple-areola complex involvement predictive index for indicating nipple-sparing mastectomy in breast cancer patients. Breast Cancer 2019;26:808-16. [Crossref] [PubMed]
- Cho J, Chung J, Cha ES, Lee JE, Kim JH. Can preoperative 3-T MRI predict nipple-areolar complex involvement in patients with breast cancer? Clin Imaging 2016;40:119-24. [Crossref] [PubMed]
- Sakamoto N, Tozaki M, Hoshi K, Fukuma E. Is MRI useful for the prediction of nipple involvement? Breast Cancer 2013;20:316-22. [Crossref] [PubMed]
- Byon W, Kim E, Kwon J, Park YL, Park C. Magnetic resonance imaging and clinicopathological factors for the detection of occult nipple involvement in breast cancer patients. J Breast Cancer 2014;17:386-92. [Crossref] [PubMed]
- Liu Z, Li X, Feng B, Li C, Chen Y, Yi L, Li Z, Li R, Long W. MIP image derived from abbreviated breast MRI: potential to reduce unnecessary sub-nipple biopsies during nipple-sparing mastectomy for breast cancer. Eur Radiol 2021;31:3683-92. [Crossref] [PubMed]
- Evans A, Whelehan P, Thomson K, McLean D, Brauer K, Purdie C, Baker L, Jordan L, Rauchhaus P, Thompson A. Invasive breast cancer: relationship between shear-wave elastographic findings and histologic prognostic factors. Radiology 2012;263:673-7. [Crossref] [PubMed]
- Wan C, Zhou L, Jin Y, Li F, Wang L, Yin W, Wang Y, Li H, Jiang L, Lu J. Strain ultrasonic elastography imaging features of locally advanced breast cancer: association with response to neoadjuvant chemotherapy and recurrence-free survival. BMC Med Imaging 2023;23:216. [Crossref] [PubMed]
- Zhu Y, Fei XC, Zhou QF, Jia XH, Zhou JQ. Prediction of the Level of Tumor-infiltrating Lymphocytes in Breast Cancer Patients with a Nomogram Based on Clinicopathologic Variables, Conventional Ultrasound and Dual-mode Elastography. Acad Radiol 2025;32:5703-12. [Crossref] [PubMed]
- Gu LS, Zhang R, Wang Y, Liu XM, Ma F, Wang JY, Sun XY, Liu MJ, Wang B, Zou SM. Characteristics of contrast-enhanced ultrasonography and strain elastography of locally advanced breast cancer. J Thorac Dis 2019;11:5274-89. [Crossref] [PubMed]
- Niu Q, Li H, Du L, Wang R, Lin J, Chen A, Jia C, Jin L, Li F. Development of a Multi-Parametric ultrasonography nomogram for prediction of invasiveness in ductal carcinoma in situ. Eur J Radiol 2024;175:111415. [Crossref] [PubMed]
- Liu Q, Tang L, Chen M. Ultrasound Strain Elastography and Contrast-Enhanced Ultrasound in Predicting the Efficacy of Neoadjuvant Chemotherapy for Breast Cancer: A Nomogram Integrating Ki-67 and Ultrasound Features. J Ultrasound Med 2022;41:2191-201. [Crossref] [PubMed]
- Huang W, Qiu Z, Mu T, Li X. A nomogram based on clinical factors for preoperative prediction of nipple involvement in breast cancer. Front Surg 2022;9:923554. [Crossref] [PubMed]
- Huang J, He N, Li J, Chen J, Guan C, Wu Y, Lu Q. Performance of contrast-enhanced cone-beam breast CT to predict nipple-areolar complex involvement in early-stage breast cancer. Eur Radiol 2026;36:45-54. [Crossref] [PubMed]

