J Cancer 2024; 15(18):6122-6134. doi:10.7150/jca.100651 This issue Cite

Research Paper

Development and Validation of a Nomogram for Axillary Lymph Node Metastasis Risk in Breast Cancer

Shijing Wang1#, He Zhang2#, Xin Wang3#, Juanhan Yu4, Qingfu Zhang4, Yiwen Zheng4, Tangbo Zhang4, Xiaoyun Mao1 Corresponding address

1. Department of Breast Surgery, The First Affiliated Hospital of China Medical University, Shenyang, Liaoning Province, 110001, China.
2. Department of Medical Imaging, Affiliated Hospital of Nantong University, Nantong, Jiangsu Province, 226001, China.
3. Department of Radiology, The First Affiliated Hospital of China Medical University, Shenyang, Liaoning Province, 110001, China.
4. Department of Pathology, The First Affiliated Hospital of China Medical University, Shenyang, Liaoning Province, 110001, China.
#Shijing Wang, He Zhang and Xin Wang contributed equally to this manuscript.

Citation:
Wang S, Zhang H, Wang X, Yu J, Zhang Q, Zheng Y, Zhang T, Mao X. Development and Validation of a Nomogram for Axillary Lymph Node Metastasis Risk in Breast Cancer. J Cancer 2024; 15(18):6122-6134. doi:10.7150/jca.100651. https://www.jcancer.org/v15p6122.htm
Other styles

File import instruction

Abstract

Graphic abstract

Purpose: Preoperative assessment of axillary lymph node (ALN) status is essential for breast cancer treatment planning. This study prospectively analyzed risk factors for ALN metastasis by comparing high-resolution computed tomography (HRCT) imaging with pathology and developed a nomogram to aid in diagnosis.

Methods: From April 2023 to May 2024, breast cancer patients confirmed by pathology participated in the study. All had chest HRCT before surgery, and ALN samples were anatomically matched to HRCT imaging and pathology. The least absolute shrinkage and selection operator (LASSO) regression helped refine metastasis features, and a nomogram was constructed using the final selected features determined by multivariate logistic regression. The nomogram's performance was evaluated with concordance index (C-index), calibration plot, and decision curve analysis, with internal validation through bootstrapping.

Results: A total of 302 ALN from 98 patients were included in this study. The predictors included in the nomogram encompassed the mean CT value, short diameter, border, and shape of ALN, as well as the Ki-67 status and histological grade of the primary tumor. The model exhibited satisfactory discrimination, with a C-index of 0.869 (95% CI: 0.826-0.912) and an AUC of 0.862 (95% CI, 0.815-0.909). The calibration curve demonstrated a high degree of concordance between the predicted and actual probabilities. The decision curve analysis demonstrated that the nomogram was clinically useful when the threshold for intervention was set at the metastasis possibility range of 1% to 86%.

Conclusion: The nomogram combined with preoperative pathology and HRCT imaging have the potential to improve the evaluation of ALN status.

Keywords: Breast cancer, Axillary lymph node, Metastasis, HRCT, Nomogram

Introduction

Breast cancer is the most prevalent and the most lethal form of cancer among women globally. The most recent global cancer burden data for 2022, released by the International Agency for Research on Cancer (IARC) of the World Health Organization, indicates that there were more than 2.3 million new cases of breast cancer and 670,000 deaths worldwide in 2022[1]. Status of axillary lymph node (ALN) is of considerable consequence for the staging, treatment planning, and prognostic assessment[2].

Pathological examination remains the most accurate method of assessing ALN status, it necessitates invasive procedures, which inevitably entails some risk[3]. Additionally, non-metastasis lymph node (NMLN) may be misclassified during surgery, leading to unnecessary dissections that don't benefit the patient and may cause complications. Furthermore, studies indicated that maintaining the integrity of NMLN can be advantageous for immunotherapy[4,5]. Consequently, there is a pressing clinical requirement for precise and non-invasive preoperative assessment method. Ultrasonography, while cost-effective and versatile, is heavily reliant on operator expertise and equipment, introducing subjectivity[6]. A limitation is its inability to differentiate between localized and advanced ALN metastasis when results are positive[7]. Additionally, detecting deep-seated or bone-adjacent ALN can be challenging due to anatomical constraints. Magnetic resonance imaging (MRI) offers superior soft-tissue resolution and can distinguish metastasis lymph node (MLN) from NMLN with contrast enhancement[8]. However, preoperative MRI is primarily employed to identify occult breast tumor and to plan breast-conserving surgery due to factors such as cost, examination time, and accessibility. Additionally, MRI has limited ability to display the entire axilla, and artifacts caused by cardiac motion can impede the detection of deep ALN[9].

Computed Tomography (CT) is commonly used to evaluate the progression and distant spread of breast cancer[10]. As the scan covers the axillary region, chest CT provides crucial information about ALN without exposing the patient to additional radiation. Furthermore, chest CT not only reveals deep ALN and Rotter's nodes but also offers a comprehensive view of both axillae, facilitating comparative analysis between the affected and unaffected sides[11]. Previous studies have indicated that standard CT scans based on size criteria can assist in preoperative ALN status assessment and enhance the accuracy of sentinel lymph node biopsy[12-14]. However, conventional CT has limitations in depicting small lymph nodes due to scan slice thickness and reconstruction algorithms[15]. Given the lack of pathology for every ALN, the majority of studies have employed an indirect approach to assess ALN status by sampling the entire specimen.

We prospectively integrated pathology with imaging by utilizing the high spatial resolution of HRCT technology combined with spiral CT. A comparative analysis was conducted on an individual ALN basis to develop a nomogram for predicting ALN status, and its performance was comprehensively evaluated.

Patients and methods

This study was approved by the Medical Ethics Committee of the First Affiliated Hospital of China Medical University (IRB-2024-2023-514). All patients provided informed consent. All methods were performed in accordance with the Declaration of Helsinki and relevant guidelines. All data generated or analyzed during this study are included in this article. No large language models were used in the preparation of this article.

Patients

From April 2023 to May 2024, 98 patients with breast cancer diagnosed via core needle biopsy (CNB) were prospectively enrolled. The inclusion criteria for the study were as follows: (1) complete pathology and HRCT imaging; (2) the HRCT examination is performed before CNB to avoid any impact on the imaging signs of the lymph nodes caused by the CNB procedure; (3) the interval between preoperative HRCT examination and surgery was less than one week, and no local or systemic treatment was received during this period. The exclusion criteria were as follows: (1) axillary fat or vascular artifacts that affected the assessment of ALN. (2) no ALN on preoperative HRCT imaging.

Image acquisition

All patients were scanned using a Siemens 64-slice 128-layer spiral CT machine in a supine position with arms raised, either head-first or feet-first, while maintaining a state of respiratory arrest at the end of inhalation. The scanning range extended from the lung apex to the diaphragm. The CT scanning parameters were set as follows: tube voltage of 120 kVp, a matrix of 512, a scanning slice thickness of 1 mm, and 1 mm thin-slice reconstruction. Upon completion of the scan, image reconstruction was performed using a high spatial frequency algorithm. Subsequently, the thin-slice reconstructed images were archived and transferred to the Picture Archive and Communication System (PACS).

Imaging-pathology comparison

The patient's basic information, HRCT imaging, and pathology (both the preoperative primary tumor pathology and the postoperative lymph node pathology) were recorded in accordance with data protection regulations. A radiologist with over 15 years of experience in breast imaging analyzed the HRCT images preoperatively, carefully selecting the region of interest (ROI) on the PACS while avoiding adjacent blood vessels and muscular tissues. The key parameters of the ALN were documented, including mean CT value, diameters, border, and shape (Figure 1). Moreover, the anatomical location of each ALN was recorded, utilizing Berg's localization criteria and referencing landmarks such as the pectoralis minor, adjacent vessels, and ribs. Anatomically, ALN are classified into three levels based on their relation to the pectoralis minor muscle. The origin of the thoracic acromial artery is located at the projected position of the inner edge of the pectoralis minor muscle, while the origin of the lateral thoracic artery is situated at the projected position of the outer edge of the pectoralis minor muscle. Level III (subclavian area) is situated medially to the inner edge of the pectoralis minor muscle and is aligned with the origin of the thoracic acromial artery. In contrast, Level I is situated lateral to the outer edge of the pectoralis minor and is marked by the starting point of the lateral thoracic artery. The area situated between the origins of the thoracic acromial and lateral thoracic arteries, which also encompasses the interpectoral (Rotter's) lymph nodes, is designated as Level II (Figure 2). All patients underwent sentinel lymph node biopsy or axillary lymph node dissection within one week of undergoing a chest HRCT examination. The fresh tissues were placed on a specimen measuring board and subsequently sent to the pathology department. A collaborative effort between pathologists, breast surgery experts, and radiology experts ensured that the specimens were matched and numbered individually based on the location and size of the lymph node corresponding to the HRCT imaging. All specimens obtained postoperatively underwent formalin fixation, paraffin embedding, and hematoxylin and eosin (HE) staining. Both micro- and macro-metastasis were considered MLN. Figure 3 depicts the complete process of comparison of imaging with pathology.

 Figure 1 

Key parameters of the HRCT imaging. The primary tumor and lymph nodes were evaluated in the transverse of the HRCT imaging, with a slice thickness of 1 mm. A. (D) diameter of the primary tumor. B. (L) long diameter of the lymph node; (S) short diameter of the lymph node. C. The mean CT value was calculated by delineating a region of interest that encompassed the entire lymph node, ensuring that the measurement line remained within the border of the lymph node. D. Reniform shape. E. Round shape. F. Other shape. G. Smooth border. H. Indistinct border. I. Fused border.

J Cancer Image
 Figure 2 

Anatomical location of axillary lymph node. A. Transverse anatomical schematic. MG: mammary gland, Pma: pectralis major muscle, Pmi: petcoralis minor muscle, LT: lateral thoracic artery (blue), TT: thoracic acromial artery (green). B. Anatomical location of axillary lymph node in a CT transverse view. The area where the level I axillary lymph node is located (red); The area where the level II axillary lymph node is located (yellow); The area where the level III axillary lymph node is located (blue). C. Coronal anatomical schematic. The pectoralis minor muscle (light green), level I axillary lymph node (red), level II axillary lymph node (yellow), level III axillary lymph node (blue). D. Anatomical location of axillary lymph node in a CT coronal view. The area where the level I axillary lymph node is located (red); The area where the level II axillary lymph node is located (yellow); The area where the level III axillary lymph node is located (blue).

J Cancer Image
 Figure 3 

Flow chart comparing imaging to pathology. CT detected a 2x1.5cm lymph node located laterally to the left pectoralis minor. The node was surgically located, and subsequent pathology confirmed breast cancer metastasis.

J Cancer Image

Statistical analysis

All analyses were performed using R version 4.4.0 (http:// www.Rproject.org). The normality of the continuous variables was evaluated using the Kolmogorov-Smirnov test. Data that adhere to a normal distribution was presented as the mean ± standard deviation (SD). For intergroup comparisons, the independent sample t test was employed. In instances where the data did not adhere to a normal distribution, they were expressed as the median accompanied by interquartile ranges (IQRs). Comparative analyses between the two groups were conducted using the Mann-whitney U test. Categorical variables were compared between groups using the chi-square test. The LASSO regression, which is particularly suitable for reducing high-dimensional data, was employed to identify the most predictive risk factors among MLN. Subsequently, a multivariate logistic regression analysis was conducted using variables that had been screened by LASSO regression to construct a predictive model. A nomogram was developed based on the logistic regression model. The performance of the nomogram was evaluated using the concordance index (C-index) and calibration plot. The receiver operator characteristic (ROC) curve and the area under the curve (AUC) were used to assess the diagnostic value. The nomogram was subjected to bootstrapping validation (1000 bootstrap resamples) to calculate a relatively corrected C-index and AUC. Decision curve analysis was conducted to determine the clinical usefulness of the nomogram by quantifying the net benefits at different threshold probabilities in the breast cancer cohort. All statistical significance levels were two-sided, with a p-value threshold of <0.05.

Results

Patient characteristics

From April 2023 to May 2024, this study prospectively enrolled 98 cases. Based on whether the patients developed ALN metastasis of breast cancer, these cases were divided into two groups: the metastasis group and the non-metastasis group. The metastasis group comprised 30 cases, while the non-metastasis group consisted of 68 cases. Table 1 presents a summary and comparison of the basic characteristics of the two groups.

Imaging features of axillary lymph node

For the purpose of the node-by-node assessment, only those ALN that could be identified were considered. HRCT identified 515 ALN, whereas pathological examination revealed 716. It is noteworthy that 302 ALN were successfully matched by both HRCT imaging and pathology. Among them, pathology confirmed 71 as MLN and 231 as NMLN (Figure 4). The mean CT value, short and long diameter of MLN significantly exceeded those of NMLN (both P<0.001). However, MLN exhibited a lower long diameter/short diameter (L/S) ratio compared to NMLN (P=0.003). The median short diameter of MLN was 6.1mm (IQRs: 3.3mm; range: 3.1-13mm), while that of NMLN was 5.0mm (IQRs: 1.8mm; range: 2-13.8mm). Similarly, the median long diameter of MLN was 8.9mm (IQRs: 4.1mm; range: 5.5-20mm) compared to 7.8mm (IQRs: 3mm; range: 4.4-21.6mm) for NMLN. The median L/S ratio for MLN was 1.45 (IQRs: 0.2; range: 1.06-2.40), whereas for NMLN, it was 1.61 (IQRs: 0.51; range: 1.03-3.20). Finally, the mean CT value for MLN was 38.99 HU (SD: 6.80HU), and for NMLN, it was 33.05 HU (SD: 6.32HU). Figure 5 depicts the histogram of lymph nodes distribution under different diameters, L/S ratio, and mean CT value. Notably, significant differences in border and shape were found between MLN and BLN (P=0.003 and P<0.001, respectively) (Table 2).

 Table 1 

Comparison of characteristics between the metastasis and non-metastasis group.

CharacteristicMetastasis group N, %Non-metastasis group N, %Total N, %
Age (years)
≤404 (13.3)7 (10.3)11 (11.2)
>4026 (86.7)61 (89.7)87 (88.8)
Menstrual State
Premenopausal11 (36.7)23 (33.8)34 (34.7)
Postmenopausal19 (63.3)45 (66.2)64 (65.3)
BMI
<18.52 (6.7)5 (7.4)7 (7.1)
18.5-2416 53.3)43 (63.2)59 (60.2)
>2412 (40)20 (29.4)32 (32.7)
Tumor diameter (cm)
≤219 (63.3)46 (67.6)65 (66.3)
>26 (20)11 (16.2)17 (17.3)
Unknown5 (16.7)11 (16.2)16 (16.4)
Histological grade
I2 (6.7)7 (10.3)9 (9.2)
II12 (40)44 (64.7)56 (57.1)
III16 (53.3)17 (25)33 (33.7)
Hormone receptor
Negative10 (33.3)16 (23.5)26 (26.5)
Positive20 (66.7)52 (76.5)72 (73.5)
HER-2
Negative23 (76.7)40 (58.8)63 (64.3)
Positive4 (13.3)20 (29.4)24 (24.5)
Unknown3 (10)8 (11.8)11 (11.2)
Ki-67 status
<20%6 (20)30 (44.1)36 (36.7)
≥20%24 (80)38 (55.9)62 (63.3)

Abbreviations: BMI-body mass index; Unknown-the information is unclear.

Feature selection

Based on the cohort comprising 98 patients with 302 lymph nodes, we initially identified 14 features spanning demographics, HRCT imaging, and pathology. As illustrated in Figure 6A and B, through the application of the LASSO regression, these features were further refined to identify 10 potential predictors. By leveraging the 1 Standard Error (SE) of the minimum criteria, we successfully streamlined the model to encompass six strongly correlated variables: mean CT value, short diameter, border, shape, Ki-67 status, and histological grade. This streamlining not only mitigates the risk of overfitting but also enhances computational efficiency, thereby bolstering the model's robustness and overall performance.

Development of an individualized prediction model

The results of the logistic regression analysis among the mean CT value, short diameter, border, shape, Ki-67 status, and histological grade are given in Table 3. The model that incorporated the above independent predictors was developed and presented as the nomogram (Figure 7).

 Figure 4 

Imaging versus the pathology in 716 lymph node of 98 patients. Abbreviations: ALN-axillary lymph node; MLN-metastasis lymph node; NMLN-non-metastasis lymph node.

J Cancer Image
 Figure 5 

The histogram of axillary lymph node distribution under different variables. A. Mean CT value; B. Short diameter; C. Long diameter; D. Long diameter/Short diameter. Abbreviations: NO.of LN- number of lymph nodes.

J Cancer Image
 Table 2 

HRCT imaging of axillary lymph node

FeatureMLN (n=71) N,%NMLN (n=231) N,%P value
Short diameter (mm), median (IQRs)6.1 (3.3)5.0 (1.8)<0.001**
Long diameter (mm), median (IQRs)8.9 (4.1)7.8 (3)<0.001**
L/S Ratio1.45 (0.2)1.61 (0.51)0.003**
Mean CT value (HU), mean ± SD38.99±6.8033.05±6.32<0.001***
Border0.003*
Smooth46 (64.8)187 (81)
Indistinct14 (19.7)34 (14.7)
Fused11 (15.5)10 (4.3)
Shape<0.001*
Reniform24 (33.8)143 (61.9)
Round29 (40.8)60 (26)
Other18 (25.4)28 (12.1)

Abbreviations: *chi-square test; **Mann-whitney U test; ***independent sample t test.

 Table 3 

Logistic multivariate analysis

VariableβStandard.ErrorZWald χ2Odds ratio (95% CI)P-value
Mean CT value0.10910.02654.117316.95251.115 (1.061-1.177)<0.001
Short diameter0.35330.08654.085816.69351.424 (1.207-1.697)<0.001
Border
SmoothReference
Indistinct0.58480.44921.3021.69511.795 (0.734-4.317)0.193
Fused1.66690.60732.74467.53295.296 (1.618-17.919)0.006
Shape
ReniformReference
Round0.89590.39412.27345.16842.45 (1.137-5.371)0.023
Other0.82570.45331.82143.31762.284 (0.932-5.563)0.069
Ki-67 status
Ki-67<20%Reference
Ki-67≥20%1.37020.4133.317811.00773.936 (1.800-9.148)<0.001
Histological grade
IReference
II0.42640.60110.70940.50331.532 (0.492-5.310)0.478
III2.00860.63413.167910.03537.453 (2.269-27.757)0.002
 Figure 6 

Selection of HRCT imaging and pathology features using the LASSO regression. A. Optimal lambda selection in the LASSO model involved fivefold cross-validation based on minimum criteria. The binomial deviance curve was plotted against log(lambda), with dotted lines indicating optimal values determined by both the minimum criteria and the 1-SE criteria (red). B. The LASSO coefficient profiles for 14 features were plotted versus log (lambda). A vertical line (red) indicates the optimal lambda chosen through fivefold cross-validation, resulting in six features with nonzero coefficients. Abbreviations: HRCT-high-resolution computed tomography; LASSO-least absolute shrinkage and selection operator; SE-standard error.

J Cancer Image
 Figure 7 

Developed prediction nomogram. The axillary lymph node metastasis nomogram was developed in the cohort, with the mean CT value, short diameter, border, shape, Ki-67 status, and the histological grade.

J Cancer Image
 Figure 8 

Calibration curves of the prediction nomogram in the cohort. The x-axis represents the predicted axillary lymph node metastasis risk. The y-axis represents the actual diagnosed axillary lymph node metastasis. The diagonal dotted line represents a perfect prediction by an ideal model. The blue dotted line represents the entire cohort, and the red solid line is bias-corrected by bootstrapping (B=1000 repetitions, boot), indicating nomogram performance.

J Cancer Image

Performance of the nomogram

The calibration curve of the nomogram for the prediction of ALN metastasis in breast cancer patients demonstrated good agreement in this cohort (Figure 8). The prediction nomogram for the cohort achieved a C-index of 0.869 (95% CI: 0.826-0.912), with bootstrapping validation confirming a score of 0.845. Likewise, the nomogram showed a diagnostic accuracy of 79.5%, an AUC of 0.862 (95% CI: 0.815-0.909), sensitivity of 78.9%, and specificity of 79.7% (Figure 9). Additional bootstrapping validation further supported an AUC of 0.857. Furthermore, the AUC determined using the mean CT value was 0.752 (95% CI, 0.684-0.811), while the AUC determined using the short diameter was 0.705 (95% CI, 0.636-0.770) (Figure 10). Taken together, these findings underscore the model's excellent discriminatory and predictive abilities.

Clinical use

Figure 11 presents the decision curve analysis for the prediction nomogram. The curve demonstrates that utilizing this nomogram to predict the risk of MLN provides greater benefit than the scheme, particularly when the threshold probability of MLN ranges from 1 to 86%. Within this range, net benefit was comparable with several overlaps, on the basis of the metastasis nomogram.

 Figure 9 

Receiver operating characteristic curve and prediction model based on nomogram. Abbreviations: AUC-area under the curve; CI-confidence interval.

J Cancer Image
 Figure 10 

Receiver operating characteristic curve and prediction model based on mean CT value, short diameter and nomogram. Abbreviations: AUC-area under the curve; CI-confidence interval.

J Cancer Image
 Figure 11 

Decision curve analysis for the prediction nomogram. The y-axis measures the net benefit. The red line represents the axillary lymph node metastasis risk nomogram. The thin solid line represents the assumption that all axillary lymph node are malignant. The thick solid line represents the assumption that no axillary lymph node are malignant.

J Cancer Image

Discussion

In clinical practice, preoperative non-invasive precise assessment of lymph node metastasis and accurate localization in breast cancer patients is an urgent issue that needs to be addressed. Preoperative imaging radiology is undoubtedly a very good tool, but it also presents significant challenges[16]. CT technology appeared in the 1970s but was rarely used for breast diseases[17-19]. However, in the 1990s, advanced spiral CT scanners improved imaging without extra radiation, aiding in assessing ALN status[20,21]. The integration of spiral CT scanners with HRCT technology offers the potential to elucidate the subtle structural and morphological alterations of lymph nodes, thereby improving the precision of ALN status assessment. In this study, 31% (71/231) of the ALN were identified as MLN. In the risk factor analysis, six characteristics were found to be associated with MLN, including the mean CT value, short diameter, border, and shape of the ALN, as well as the Ki-67 status and histological grade of the primary tumor. The prediction nomogram indicated that higher mean CT value, larger short diameter, fused border, round shape, Ki-67≥20%, and histological grade III may be the key individual factors associated with MLN.

Hounsfield units (HU) serve as a quantifiable measure of X-ray attenuation in CT imaging, offering an objective and reproducible metric that is less prone to subjective examiner interpretation[22,23]. It is generally believed that pathological parameters, including cell count, tumor stroma, and/or extracellular matrix, may have a correlation with CT value[24,25]. As previous study indicated that HU can accurately distinguished breast masses as cystic or solid with or without contrast with chest CT[26]. Urata et al. conducted a retrospective, dual-center study on breast cancer patients who underwent contrast-enhanced CT assessment. They discovered that CT value can predict ALN status with sensitivity, specificity, and accuracy levels of 79.6%, 80.5% and 80.2%, respectively[27]. Similarly, Zhang et al.'s research indicated that the venous phase slope of the spectral HU curve is the optimal single parameter for detecting metastatic sentinel lymph nodes[28]. In our study, the nomogram that we established indicated that the mean CT value carried the greatest weight, implying that lymph nodes with a higher mean CT value have a higher likelihood of metastasis (Figure 7). Furthermore, the AUC, determined using mean CT value, stood at 0.752 (95% CI, 0.684-0.811). The optimal cut-off value was identified as 41.2 HU, which yielded a sensitivity of 73.2%, a specificity of 71.4%, and an accuracy of 71.9% (Figure 10). Assessing lymph node status by CT value is advantageous due to its objectivity and availability during preoperative exams. However, variations in CT scanners, imaging parameters, and software across institutions can affect results. Thus, when predicting ALN status with CT value, the cut-off may only be accurate for patients examined at the same institution.

In clinical practice, the size of lymph nodes serves as a critical indicator of their status[29]. From a pathological perspective, MLN enlarge due to the infiltration of tumor cells and obstructed lymphatic drainage[30]. Numerous imaging studies have also demonstrated that the diameters of MLN are significantly larger than those of NMLN[31-33]. This study further supported these observations. Our findings showed that the median short diameter of MLN was 6.1mm (IQRs: 3.3mm), while that of NMLN was only 5.0mm (IQRs: 1.8mm) (P<0.001). Similarly, the median long diameter of MLN is 8.9mm (IQRs: 4.1mm), which is significantly higher than the 7.8mm (3.0mm) of NMLN (P<0.001). Lymph node diameters may vary in imaging due to angle changes relative to the imaging plane. Compared to the long diameter, the short diameter offers a more stable, accurate representation of the size and higher measurement repeatability, as it is less affected by angle changes[34-36]. The findings of this study revealed that the long diameter is not an independent predictor of ALN status. Instead, the short diameter proved superior to the long diameter in anticipating ALN status (AUC: 0.702 vs. 0.637). The established nomogram showed that the short diameter was included with a relatively high weight, which means that ALN with larger short diameters are more likely to metastasize. These results align with the Response Evaluation Criteria in Solid Tumors 1.1, which recommends using the short diameter to distinguish lymph nodes status[37].

However, there is still no consensus on the optimal cut-off value for the short diameter of MLN, which varies from 5mm to 10.4mm across different tumor types [12,29,38-42]. Despite data differences among study groups, the diagnostic trend was uniform: a larger short diameter of lymph nodes led to higher positive predictive value and specificity in metastasis diagnosis, but lower negative predictive value and sensitivity. In our study, the optimal cut-off value for the short diameter in predicting MLN was found to be 7.4mm, yielding a sensitivity of 70.4%, a specificity of 62.3%, and an accuracy of 64.2% (Figure 10). Furthermore, studies have demonstrated a correlation between lymph nodes (including ALN) with a L/S ratio less than 1.5 and the presence of metastasis[43-45]. In this study, we found that the L/S ratio of MLN was significantly lower than that of BLN, with 53.5% (38/71) of MLN exhibiting an L/S ratio below 1.5. It is noteworthy that while univariate analyses revealed significant differences between long diameter and L/S ratio, LASSO regression indicated that these two factors exhibited diminished significance when considered alongside other variables. This finding was further corroborated by subsequent multifactor logistic regression analysis. This indicated that the independent effects of long diameter and L/S ratio on the results are diminished when interactions between multiple variables are considered. They may be masked by other stronger influences or covariates that interact with other factors.

Besides CT value and size, the border and shape of lymph nodes are crucial for assessing lymph node status[46]. The NMLN typically exhibit a reniform shape, smooth border, an absence of focal prominence, and a clearly visible central fatty hilum[47,48]. The immune response triggered by cancer cells, as well as the process of their colonization and proliferation, may result in significant morphological and border changes in ALN[49]. With its high spatial resolution, HRCT is capable of accurately capturing these subtle changes, which in turn provides substantial support in assessing the status of the lymph nodes. Similar to previous studies, we found that MLN exhibited an odds ratio (OR) of 5.296 (95% CI: 1.618-17.919) for fused border and an OR of 2.45 (95% CI: 1.137-5.371) for a round shape in comparison to NMLN. Different from previous studies, indistinct border and other shape (such as lobulated or focally prominent) may not necessarily indicate MLN. Perhaps the border and shape of ALN are also affected by various factors, including non-neoplastic inflammation and fat content[50-54]. To optimize diagnostic accuracy, considering multiple factors like tumor metabolic intensity and inflammatory indicators is imperative.

Scholars believe that correlating the biological behavior of lymph node metastasis with the pathological characteristics of the primary tumor can not only improve the accuracy of judging lymph node metastasis but also provide clinical support for finding new targets for breast cancer treatment[55,56]. Previous studies have discovered a strong association between larger primary tumor diameter, vascular invasion, triple-negative breast cancer, and MLN[57-60]. However, most of these risk factors were based on postoperative pathological results, making them potentially unsuitable for the preoperative evaluation of ALN status. In our study, all histological features of the primary tumor were derived from preoperative CNB[61-63]. High histological grade (grade III) and Ki-67≥20% of the primary tumor were identified as independent predictors of MLN, both of which carried relatively high weights in the nomogram. However, characteristics such as HR, HER-2 and tumor diameter did not show significant differences. This finding may be due in part to the lack of HER2 information for some patients in this study (11/98), as well as the failure of CT scans to detect the primary tumor (16/98).

To the best of our knowledge, this is the first prospective study to create a nomogram using preoperative pathology and HRCT imaging for accurate preoperative ALN status evaluation. We developed and validated this prediction nomogram using six accessible variables to assess the status of each ALN in patients with breast cancer. The nomogram boasts user-friendly digital interface, enhanced accuracy, and prognostic reports that are simpler to comprehend, all aimed at facilitating improved clinical decision-making.

Limitations

The present study was limited by several factors. Primarily, despite extensive internal validation via bootstrap testing, our nomogram lacks external validation. Its applicability to other breast cancer populations remains uncertain, necessitating further external evaluation in broader contexts. Secondly, establishing a one-to-one correspondence between excised lymph nodes and HRCT imaging of ALN proved challenging, resulting in a relatively small number of matched ALN, particularly MLN. Thirdly, despite efforts to match all histologically identified lymph nodes, there was a notable discrepancy in the number of nodes within different short diameters, which could potentially impact the accuracy of the results.

Conclusions

This study provided a relatively accurate prediction tool of ALN status for breast cancer patients. However, these results should be interpreted cautiously as further extensive prospective studies are required to validate the diagnostic efficacy of the nomogram.

Acknowledgements

All authors have read and approved the manuscript.

Funding

This work was funded by the National Natural Science Foundation of China (grant numbers 12374413)

Author contributions

Shijing Wang, He Zhang, Xin Wang, Juanhan Yu, Qingfu Zhang, and Xiaoyun Mao wrote the main manuscript text, Shijing Wang and He Zhang counted the sample size, Shijing Wang, He Zhang, Xin Wang and Xiaoyun Mao participated in the design of the research scheme, Yiwen Zheng and Tangbo Zhang participated in reviewing and revising the language.

Ethics declarations

This study was approved by the Medical Ethics Committee of the First Affiliated Hospital of China Medical University (IRB-2024-2023-514). All patients provided informed consent. All methods were performed in accordance with the Declaration of Helsinki and relevant guidelines.

Data availability

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

Competing Interests

The authors have declared that no competing interest exists.

References

1. Bray F, Laversanne M, Sung H. et al. Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2024;74:229-263

2. Wang G, Zhang S, Wang M. et al. Prognostic significance of occult lymph node metastases in breast cancer: a meta-analysis. BMC Cancer. 2021;21:875

3. Pilger TL, Francisco DF, Candido Dos Reis FJ. Effect of sentinel lymph node biopsy on upper limb function in women with early breast cancer: A systematic review of clinical trials. Eur J Surg Oncol. 2021;47:1497-506

4. Rahim MK, Okholm TLH, Jones KB. et al. Dynamic CD8+ T cell responses to cancer immunotherapy in human regional lymph nodes are disrupted in metastatic lymph nodes. Cell. 2023;186:1127-1143.e18

5. Mj O, Mp M, Sn T. Tumor-draining lymph nodes are survival niches that support T cell priming against lymphatic transported tumor antigen and effects of immune checkpoint blockade in TNBC. Cancer immunology, immunotherapy : CII. 2021;70:2179-2185

6. Marino MA, Avendano D, Zapata P, Riedl CC, Pinker K. Lymph Node Imaging in Patients with Primary Breast Cancer: Concurrent Diagnostic Tools. Oncologist. 2020;25:e231-42

7. Schipper RJ, van Roozendaal LM, de Vries B. et al. Axillary ultrasound for preoperative nodal staging in breast cancer patients: is it of added value? Breast. 2013;22:1108-13

8. Scheenen TWJ, Zamecnik P. The Role of Magnetic Resonance Imaging in (Future) Cancer Staging: Note the Nodes. Invest Radiol. 2021;56:42-9

9. Zhou P, Wei Y, Chen G, Guo L, Yan D, Wang Y. Axillary lymph node metastasis detection by magnetic resonance imaging in patients with breast cancer: A meta-analysis. Thorac Cancer. 2018;9:989-96

10. Desperito E, Schwartz L, Capaccione KM. et al. Chest CT for Breast Cancer Diagnosis. Life (Basel). 2022;12:1699

11. Imai N, Kitayama M, Shibahara A. et al. Strategy for the accurate preoperative evaluation of the number of metastatic axillary lymph nodes in breast cancer. Asian J Surg. 2019;42:228-34

12. Shien T, Akashi-Tanaka S, Yoshida M. et al. Evaluation of axillary status in patients with breast cancer using thin-section CT. Int J Clin Oncol. 2008;13:314-9

13. Yasuta M, Sato S, Ishida T, Kiyohara T. Usefulness of CT-lymphography in sentinel lymph node navigation. Int J Clin Oncol. 2014;19:557-62

14. Fujita T, Miura H, Seino H. et al. Anatomical classification of breast sentinel lymph nodes using computed tomography-lymphography. Anat Sci Int. 2018;93:487-94

15. Ou X, Zhu J, Qu Y. et al. Imaging features of sentinel lymph node mapped by multidetector-row computed tomography lymphography in predicting axillary lymph node metastasis. BMC Med Imaging. 2021;21:193

16. Valente SA, Levine GM, Silverstein MJ. et al. Accuracy of predicting axillary lymph node positivity by physical examination, mammography, ultrasonography, and magnetic resonance imaging. Ann Surg Oncol. 2012;19:1825-30

17. McCullough EC, Payne JT. X-ray-transmission computed tomography. Med Phys. 1977;4:85-98

18. Chang CH, Sibala JL, Fritz SL, Dwyer SJ, Templeton AW. Specific value of computed tomographic breast scanner (CT/M) in diagnosis of breast diseases. Radiology. 1979;132:647-52

19. Chang CH, Sibala JL, Gallagher JH. et al. Computed tomography of the breast. A preliminary report. Radiology. 1977;124:827-9

20. Teifke A, Schweden F, Cagil H, Kauczor HU, Mohr W, Thelen M. [Spiral computerized tomography of the breast]. Rofo. 1994;161:495-500

21. Imhof H, Schibany N, Ba-Ssalamah A. et al. Spiral CT and radiation dose. Eur J Radiol. 2003;47:29-37

22. Zhou Z, Liu Y, Meng K. et al. Application of spectral CT imaging in evaluating lymph node metastasis in patients with gastric cancers: initial findings. Acta Radiol. 2019;60:415-24

23. Koezuka S, Sano A, Azuma Y. et al. Combination of mean CT value and maximum CT value as a novel predictor of lepidic predominant lesions in small lung adenocarcinoma presenting as solid nodules. Sci Rep. 2022;12:5450

24. Wienbeck S, Fischer U, Perske C. et al. Cone-beam Breast Computed Tomography: CT Density Does Not Reflect Proliferation Potential and Receptor Expression of Breast Carcinoma. Transl Oncol. 2017;10:599-603

25. Uhlig J, Fischer U, von Fintel E. et al. Contrast Enhancement on Cone-Beam Breast-CT for Discrimination of Breast Cancer Immunohistochemical Subtypes. Transl Oncol. 2017;10:904-10

26. Capaccione K, Desperito E, Asiimwe AC, Salvatore M. Can Hounsfield units on chest CT characterize breast nodules as cystic or solid? Transl Breast Cancer Res. 2024;5:6

27. Urata M, Kijima Y, Hirata M. et al. Computed tomography Hounsfield units can predict breast cancer metastasis to axillary lymph nodes. BMC Cancer. 2014;14:730

28. Zhang X, Zheng C, Yang Z. et al. Axillary Sentinel Lymph Nodes in Breast Cancer: Quantitative Evaluation at Dual-Energy CT. Radiology. 2018;289:337-46

29. Torabi M, Aquino SL, Harisinghani MG. Current concepts in lymph node imaging. J Nucl Med. 2004;45:1509-18

30. Yoshimura G, Sakurai T, Oura S. et al. Evaluation of Axillary Lymph Node Status in Breast Cancer with MRI. Breast Cancer. 1999;6:249-58

31. Iima M, Kataoka M, Okumura R, Togashi K. Detection of axillary lymph node metastasis with diffusion-weighted MR imaging. Clin Imaging. 2014;38:633-6

32. Fujioka T, Mori M, Kubota K. et al. Clinical Usefulness of Ultrasound-Guided Fine Needle Aspiration and Core Needle Biopsy for Patients with Axillary Lymphadenopathy. Medicina (Kaunas). 2021;57:722

33. Le-Petross HT, McCall LM, Hunt KK. et al. Axillary Ultrasound Identifies Residual Nodal Disease After Chemotherapy: Results From the American College of Surgeons Oncology Group Z1071 Trial (Alliance). AJR Am J Roentgenol. 2018;210:669-76

34. Li Y, Liu K, Ke Y. et al. Risk Factors Analysis of Pathologically Confirmed Cervical Lymph Nodes Metastasis in Oral Squamous Cell Carcinoma Patients with Clinically Negative Cervical Lymph Node: Results from a Cancer Center of Central China. J Cancer. 2019;10:3062-9

35. Wang D, Zhuang Z, Wu S. et al. A Dual-Energy CT Radiomics of the Regional Largest Short-Axis Lymph Node Can Improve the Prediction of Lymph Node Metastasis in Patients With Rectal Cancer. Front Oncol. 2022;12:846840

36. Guo Q, Dong Z, Jiang L, Zhang L, Li Z, Wang D. Assessing Whether Morphological Changes in Axillary Lymph Node Have Already Occurred Prior to Metastasis in Breast Cancer Patients by Ultrasound. Medicina (Kaunas). 2022;58:1674

37. Schwartz LH, Bogaerts J, Ford R. et al. Evaluation of lymph nodes with RECIST 1.1. Eur J Cancer. 2009;45:261-7

38. Lam TBL. Optimizing the diagnosis of pelvic lymph node metastasis in bladder cancer using computed tomography and magnetic resonance imaging. Cancer Commun (Lond). 2018;38:2

39. Li Y, Diao F, Shi S. et al. Computed tomography and magnetic resonance imaging evaluation of pelvic lymph node metastasis in bladder cancer. Chin J Cancer. 2018;37:3

40. Li ZX, Liu XB, Li Y. et al. [Application value of CT examination of lymph node short diameter in evaluating cardia-left gastric lymph node metastasis in thoracic esophageal squamous cell carcinoma]. Zhonghua Zhong Liu Za Zhi. 2023;45:962-6

41. Nagano H, Takumi K, Nakajo M. et al. Dual-Energy CT-Derived Electron Density for Diagnosing Metastatic Mediastinal Lymph Nodes in Non-Small Cell Lung Cancer: Comparison With Conventional CT and FDG PET/CT Findings. AJR Am J Roentgenol. 2022;218:66-74

42. Hong EK, Chalabi M, Landolfi F. et al. Colon cancer CT staging according to mismatch repair status: Comparison and suggestion of imaging features for high-risk colon cancer. Eur J Cancer. 2022;174:165-75

43. Feu J, Tresserra F, Fábregas R. et al. Metastatic breast carcinoma in axillary lymph nodes: in vitro US detection. Radiology. 1997;205:831-5

44. Liu Y, Liu SL, Jia WQ, Wang YY, Wang J. [Value of endobronchial ultrasound-guided transbronchial needle aspiration in diagnosis of the mediastinal enlarged lymph nodes]. Zhonghua Yi Xue Za Zhi. 2016;96:3067-71

45. Tohnosu N, Onoda S, Isono K. Ultrasonographic evaluation of cervical lymph node metastases in esophageal cancer with special reference to the relationship between the short to long axis ratio (S/L) and the cancer content. J Clin Ultrasound. 1989;17:101-6

46. Farrell TPJ, Adams NC, Stenson M. et al. The Z0011 Trial: Is this the end of axillary ultrasound in the pre-operative assessment of breast cancer patients? Eur Radiol. 2015;25:2682-7

47. Bedi DG, Krishnamurthy R, Krishnamurthy S. et al. Cortical morphologic features of axillary lymph nodes as a predictor of metastasis in breast cancer: in vitro sonographic study. AJR Am J Roentgenol. 2008;191:646-52

48. Yang WT, Chang J, Metreweli C. Patients with breast cancer: differences in color Doppler flow and gray-scale US features of benign and malignant axillary lymph nodes. Radiology. 2000;215:568-73

49. Cho P, Park CS, Park GE, Kim SH, Kim HS, Oh S-J. Diagnostic Usefulness of Diffusion-Weighted MRI for Axillary Lymph Node Evaluation in Patients with Breast Cancer. Diagnostics (Basel). 2023;13:513

50. van der Valk P, Meijer CJ. The histology of reactive lymph nodes. Am J Surg Pathol. 1987;11:866-82

51. diFlorio-Alexander RM, Song Q, Dwan D. et al. Fat-enlarged axillary lymph nodes are associated with node-positive breast cancer in obese patients. Breast Cancer Res Treat. 2021;189:257-67

52. Leong PW, Chotai NC, Kulkarni S. Imaging Features of Inflammatory Breast Disorders: A Pictorial Essay. Korean J Radiol. 2018;19:5-14

53. Thierry GR, Gentek R, Bajenoff M. Remodeling of reactive lymph nodes: Dynamics of stromal cells and underlying chemokine signaling. Immunol Rev. 2019;289:42-61

54. Kooreman LFS, Dieleman S, van Kuijk SMJ, Zur Hausen A, Smidt ML, Grabsch HI. The prognostic value of the histological shape of tumor negative sentinel nodes in breast cancer. Front Immunol. 2023;14:1258641

55. Ping J, Liu W, Chen Z, Li C. Lymph node metastases in breast cancer: Mechanisms and molecular imaging. Clin Imaging. 2023;103:109985

56. Solanki M, Visscher D. Pathology of breast cancer in the last half century. Hum Pathol. 2020;95:137-48

57. Xiong J, Zuo W, Wu Y. et al. Ultrasonography and clinicopathological features of breast cancer in predicting axillary lymph node metastases. BMC Cancer. 2022;22:1155

58. Li Y, Han D, Shen C, Duan X. Construction of a comprehensive predictive model for axillary lymph node metastasis in breast cancer: a retrospective study. BMC Cancer. 2023;23:1028

59. Chen W, Wang C, Fu F, Yang B, Chen C, Sun Y. A Model to Predict the Risk of Lymph Node Metastasis in Breast Cancer Based on Clinicopathological Characteristics. Cancer Manag Res. 2020;12:10439-47

60. Xie X, Xiong Z, Li X. et al. Nomogram to Predict Internal Mammary Lymph Nodes Metastasis in Patients With Breast Cancer. Front Oncol. 2019;9:1193

61. Cipolla C, Fricano S, Vieni S. et al. Validity of needle core biopsy in the histological characterisation of mammary lesions. Breast. 2006;15:76-80

62. Verkooijen HM, Core Biopsy After Radiological Localisation (COBRA) Study Group. Diagnostic accuracy of stereotactic large-core needle biopsy for nonpalpable breast disease: results of a multicenter prospective study with 95% surgical confirmation. Int J Cancer. 2002;99:853-9

63. Vieni S, Graceffa G, La Mendola R. et al. Application of a predictive model of axillary lymph node status in patients with sentinel node metastasis from breast cancer. A retrospective cohort study. Int J Surg. 2016;35:58-63

Author contact

Corresponding address Corresponding author: Xiaoyun Mao, https://orcid.org/0000-0003-4755-0301. Department of Breast Surgery, The First Affiliated Hospital of China Medical University, 155 Nanjing North Street, Heping District, Shenyang City, Liaoning Province, 110001, China. Email: xymaoedu.cn; Telephone number: 0086 24 83282618; Fax number: 0086 24 83282618.


Received 2024-7-8
Accepted 2024-8-31
Published 2024-10-7


Citation styles

APA
Wang, S., Zhang, H., Wang, X., Yu, J., Zhang, Q., Zheng, Y., Zhang, T., Mao, X. (2024). Development and Validation of a Nomogram for Axillary Lymph Node Metastasis Risk in Breast Cancer. Journal of Cancer, 15(18), 6122-6134. https://doi.org/10.7150/jca.100651.

ACS
Wang, S.; Zhang, H.; Wang, X.; Yu, J.; Zhang, Q.; Zheng, Y.; Zhang, T.; Mao, X. Development and Validation of a Nomogram for Axillary Lymph Node Metastasis Risk in Breast Cancer. J. Cancer 2024, 15 (18), 6122-6134. DOI: 10.7150/jca.100651.

NLM
Wang S, Zhang H, Wang X, Yu J, Zhang Q, Zheng Y, Zhang T, Mao X. Development and Validation of a Nomogram for Axillary Lymph Node Metastasis Risk in Breast Cancer. J Cancer 2024; 15(18):6122-6134. doi:10.7150/jca.100651. https://www.jcancer.org/v15p6122.htm

CSE
Wang S, Zhang H, Wang X, Yu J, Zhang Q, Zheng Y, Zhang T, Mao X. 2024. Development and Validation of a Nomogram for Axillary Lymph Node Metastasis Risk in Breast Cancer. J Cancer. 15(18):6122-6134.

This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions.
Popup Image