The value of multi-parametric ultrasonography nomogram in predicting nipple-areolar complex involvement in breast cancer
Original Article

The value of multi-parametric ultrasonography nomogram in predicting nipple-areolar complex involvement in breast cancer

Yuqin Ma#, Siqi Tang#, Jingwen Zhang#, Yuanyuan Ma, Qingshu Lian, Chunchun Jin, Jianghao Lu, Peng Zhou

Department of Ultrasound, The First Affiliated Hospital of Shenzhen University Health Science Center, Shenzhen Second People’s Hospital, Shenzhen, China

Contributions: (I) Conception and design: Yuqin Ma, S Tang, J Zhang; (II) Administrative support: P Zhou; (III) Provision of study materials or patients: Yuanyuan Ma, Q Lian; (IV) Collection and assembly of data: Yuqin Ma, S Tang, J Lu; (V) Data analysis and interpretation: Yuqin Ma, C Jin; (VI) Manuscript writing: All authors; (VII) Final approval of manuscript: All authors.

#These authors contributed equally to this work.

Correspondence to: Peng Zhou, MD. Department of Ultrasound, The First Affiliated Hospital of Shenzhen University Health Science Center, Shenzhen Second People’s Hospital, No. 3002, Sungang West Road, Futian District, Shenzhen 518000, China. Email: zpdahai@163.com.

Background: Accurate assessment of nipple-areolar complex (NAC) involvement in breast cancer (BC) plays a key role in guiding preoperative surgical strategy. This study developed and evaluated a multi-parametric ultrasonography nomogram, incorporating conventional ultrasound (CUS), strain elastography (SE), and contrast-enhanced ultrasound (CEUS) to predict NAC involvement in BC.

Methods: In this prospective study, consecutive patients were enrolled if they had a preoperative biopsy confirming primary BC and were scheduled for surgical treatment. All patients underwent CUS, SE, and CEUS examinations preoperatively. Histopathological findings of NAC after surgical resection served as the gold standard. Least absolute shrinkage and selection operator (LASSO) regression was applied to identify the predictors most strongly associated with NAC involvement. The diagnostic performance of different model combinations (CUS, CUS + SE, CUS + SE + CEUS) was compared using receiver operating characteristic (ROC) analysis with the DeLong test. A multi-parametric ultrasonography nomogram was developed and internally validated by 500 bootstrap resamples, with calibration curve analysis assessing its predictive accuracy and decision curve analysis (DCA) evaluating its clinical utility.

Results: A total of 166 lesions were enrolled from 164 patients, comprising 33 (19.9%) lesions with NAC involvement and 133 (80.1%) lesions with NAC non-involvement. LASSO regression identified four independent predictors: ductal echo extending from the lesion on CUS, the mean elastography strain value of nipple on SE, and ductal enhancement extending to the nipple and focal nipple enhancement on CEUS. The combined CUS + SE + CEUS model demonstrated the highest diagnostic efficacy with an area under the curve (AUC) of 0.917, a sensitivity of 90.91%, and a negative predictive value (NPV) of 97.32%. Its performance was significantly better than that of CUS alone (AUC =0.732, P<0.001) and CUS + SE (AUC =0.859, P=0.009). The nomogram based on multi-parametric ultrasonography showed good calibration and internal validation (bias corrected AUC =0.916). DCA confirmed net benefit across clinically relevant thresholds.

Conclusions: The multi-parametric ultrasonography nomogram demonstrated excellent performance in predicting NAC involvement and offers a valuable, non-invasive tool for preoperative surgical planning.

Keywords: Breast cancer (BC); strain elastography (SE); contrast-enhanced ultrasound (CEUS); nipple-areolar complex involvement (NAC involvement); nomogram


Submitted Jan 19, 2026. Accepted for publication May 25, 2026. Published online Jun 16, 2026.

doi: 10.21037/qims-2026-1-0142


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.

Figure 1 Flowchart of patient selection process. CEUS, contrast-enhanced ultrasound; CUS, conventional ultrasound; NAC, nipple-areolar complex; SE, strain elastography.

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

Clinical characteristics of patients

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

Sonographic features on CUS between patients with NAC involvement and non-involvement

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

Sonographic features on SE between patients with NAC involvement and non-involvement

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

Sonographic features on CEUS between patients with NAC involvement and non-involvement

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.

Figure 2 Feature selection using LASSO regression. (A) LASSO coefficient profiles of the sonographic features. The vertical dashed line indicates the optimal value of log(λ) selected by 10-fold cross-validation, which resulted in 4 features with nonzero coefficients. (B) Binomial deviance across log(λ) values. The vertical dashed lines correspond to the value that minimized the deviance, confirming the selection of 4 diagnostic features. LASSO, least absolute shrinkage and selection operator.
Figure 3 A 54-year-old woman with a mass located in the left upper outer quadrant of the breast, with invasive ductal carcinoma and NAC involvement confirmed by pathology. (A) Gray-scale ultrasonography demonstrated the ductal echo extending from the lesion (arrow). (B) Color Doppler revealed that the ductal echo extending from the lesion (arrow) appeared to have punctate blood flow signals. (C) SE demonstrated that the local lesion and the nipple (arrow) presented firm texture. (D) On arterial-phase CEUS, focal nipple enhancement (white arrow) and ductal enhancement extending to the nipple (yellow arrow) were detected. CEUS, contrast-enhanced ultrasound; SE, strain elastography.

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).

Figure 4 ROC curve of different models for predicting NAC involvement. AUC, area under the curve; CEUS, contrast-enhanced ultrasound; CUS, conventional ultrasound; NAC, nipple-areolar complex; ROC, receiver operating characteristic; SE, strain elastography.

Table 5

Diagnostic efficacy of CUS, CUS + SE, and CUS + SE + CEUS models in predicting NAC involvement

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.

Figure 5 Nomogram with ductal echo extending from the lesion, Em-nipple, DEEN, and FNE. DEEN, ductal enhancement extending to the nipple; Em, the mean elastography strain value; FNE, focal nipple enhancement; NAC, nipple-areolar complex.
Figure 6 Performance evaluation of the nomogram. (A) Internal validation ROC curve. The shaded blue area represented the bootstrap 95% confidence interval for the AUC. (B) Calibration plot showed good consistency between the predicted probability and the actual observed frequency. (C) DCA for the nomogram demonstrated a superior net benefit across a wide range of clinically relevant threshold probabilities. AUC, area under the curve; DCA, decision curve analysis; ROC, receiver operating characteristic.

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 the Shenzhen Second People’s Hospital Clinical Research Fund of Shenzhen High-level Hospital Construction Project (Nos. 20203357001 and 2023yjlcyj003), the Shenzhen Health Economics Association Scientific Research Fund (No. 2023107), and the Shenzhen Key Medical Discipline Construction Fund (No. SZXK052).

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

  1. 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]
  2. 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]
  3. 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]
  4. 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]
  5. 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]
  6. 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]
  7. 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]
  8. 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]
  9. 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]
  10. 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]
  11. 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]
  12. 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]
  13. 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]
  14. 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]
  15. 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]
  16. 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.
  17. 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]
  18. 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]
  19. 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]
  20. 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]
  21. 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]
  22. 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]
  23. 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]
  24. 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]
  25. 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]
  26. 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]
  27. 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]
  28. 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]
  29. 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]
  30. 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]
  31. 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]
  32. 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]
  33. 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]
  34. 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]
  35. 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]
  36. 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]
Cite this article as: Ma Y, Tang S, Zhang J, Ma Y, Lian Q, Jin C, Lu J, Zhou P. The value of multi-parametric ultrasonography nomogram in predicting nipple-areolar complex involvement in breast cancer. Quant Imaging Med Surg 2026;16(7):556. doi: 10.21037/qims-2026-1-0142

Download Citation