ORIGINAL RESEARCH article

Front. Oncol., 04 January 2024

Sec. Head and Neck Cancer

Volume 13 - 2023 | https://doi.org/10.3389/fonc.2023.1174457

Development and validation of a nomogram to predict cervical lymph node metastasis in head and neck squamous cell carcinoma

  • 1. Department of Radiation Oncology, Affiliated Hospital of Qingdao University, Qingdao, China

  • 2. Department of Oncology and Radiotherapy, Affiliated Hospital of Qingdao University, Qingdao, China

Article metrics

View details

1

Citations

1,9k

Views

879

Downloads

Abstract

Background:

Head and neck cancers are a heterogeneous, aggressive, and genetically complex collection of malignancies of the oral cavity, nasopharynx, oropharynx, hypopharynx, larynx, paranasal sinuses and salivary glands, which are difficult to treat. Regional lymph nodes metastasis is a significant poor prognosis factor for head and neck squamous cell carcinoma. Metastasis to the regional lymph nodes reduces the 5-year survival rate by 50% compared with that of patients with early-stage disease. Accurate evaluation of cervical lymph node is a vital component in the overall treatment plan for patients with squamous cell carcinoma of the head and neck. However, current models are struggle to accurately to predict cervical lymph node metastasis. Here, we analyzed the clinical, imaging, and pathological data of 272 patients with HNSCC confirmed by postoperative pathology and sought to develop and validate a nomogram for prediction of lymph node metastasis in patients with head and neck squamous cell carcinoma.

Methods:

We retrospectively analyzed the clinical, imaging, and pathological data of 272 patients with head and neck squamous cell carcinoma (HNSCC) confirmed by postoperative pathology at the Affiliated Hospital of Qingdao University from June 2017 to June 2021. Patients were randomly divided into the training and validation cohorts in a 3:1 ratio, and after screening risk factors by logistic regression, nomogram was developed for predicting lymph nodes metastasis, then the prediction model was verified by C-index, area under curve (AUC), and calibration curve.

Results:

Of the 272 patients, seven variables were screened to establish the predictive model, including the differentiation degree of the tumor [95% confidence interval(CI):1.224~6.735, P=0.015], long-to-short axis ratio of the lymph nodes (95%CI: 0.019~0.217, P<0.001), uneven/circular enhancement (95%CI: 1.476~16.715, P=0.010), aggregation of lymph nodes (95%CI:1.373~10.849, P=0.010), inhomogeneous echo (95%CI: 1.337~23.389, P=0.018), unclear/absent medulla of lymph nodes (95%CI: 2.514~43.989, P=0.001), and rich blood flow (95%CI: 1.952~85.632, P=0.008). The C-index was 0.910, areas under the curve of training cohort and verification cohort were 0.953 and 0.938 respectively, indicating the discriminative ability of this nomogram. The calibration curve showed a favorable compliance between the prediction of the model and actual observations. The clinical decision curve showed this model is clinically useful and had better discriminative ability between 0.25 and 0.9 for the probability of cervical LNs metastasis.

Conclusions:

We established a good prediction model for cervical lymph node metastasis in head and neck squamous cell carcinoma patients which can provide reference value and auxiliary diagnosis for clinicians in making neck management decisions of HNSCC patients.

1 Introduction

Lymph nodes(LNs) metastases is an important prognosis indicator for patients with head and neck squamous cell carcinoma(HNSCC) (1). According to the 2020 National Comprehensive Cancer Network guidelines, for HNSCC patients with LNs metastases, radical cervical LNs dissection and necessary postoperative chemoradiotherapy are required, while elective dissection or non-dissection of cervical LNs is more appropriate for patients without LNs metastases. Imaging is an effective tool for identifying lymph node metastases, but no single imaging technique can accurately diagnose, stage, and provide long-term monitoring of HNSCC. The combination of ultrasound and computed tomography(CT) may be helpful, but neither of them has 100% diagnostic sensitivity (2). Therefore, a comprehensive and systematic approach based on previous research is needed to maximize the benefit of preoperative imaging techniques. As a chart, nomogram could assess the incidence of clinical events quickly, visually, and accurately. Their validity has been demonstrated in assessing the probability of cervical LNs metastasis in some HNSCCs (including tongue, larynx, and hypopharynx) (35). However, they only evaluated the correlation between clinicopathological features and cervical LNs metastasis.

In this study, we constructed a novel predict model and validated it for predicting cervical LNs metastasis in patients with HNSCC based on the incorporated imaging features of cervical LNs and pathological features of primary tumor. Compared to other models, our study included more variables and exhibited superior discrimination capabilities. According to the risk stratification based on this nomogram, this prediction model can provide clinicians with aid in diagnosis, high-risk patients may profit from neck management, but low-risk patients may be conservative. We present this article in accordance with the TRIPOD reporting checklist.

2 Methods

2.1 Patient selection

This was a retrospective study compliant with the Declaration of Helsinki (as revised in 2013). All included patients gave their oral and written informed consent. This study was approved by the Institutional Review Board of the affiliated hospital of Qingdao University (No. QYFY WZLL 27167). We collected patients with HNSCC confirmed by postoperative pathology at the Affiliated Hospital of Qingdao University from June 2017 to June 2021. After acquiring complete information, following inclusion criteria were as follows: 1) patients with pathologically confirmed HNSCC who underwent LNs dissection; 2) age ≥18 years; 3) karnofsky performance status score ≥70; 4) patients with complete data who undergo contrast-enhanced CT and cervical ultrasound within 2 weeks before surgery in our hospital. The exclusion criteria included: 1) patients diagnosed with multiple primary malignant tumors; 2) patients who received any other anti-tumor treatments; 3) patients with a history of liver or kidney dysfunction; 4) patients whose CT or ultrasound presented poor image quality. All patients received tumor extended resection and elective neck dissection. The outcome indicator was cervical LNs metastasis reviewed by two independent pathologists.

2.2 Imaging and pathological analysis

A radiotherapy doctor, a radiologist, and a sonographer, each with over 5 years of experience in head and neck cancer imaging, respectively finished their training using randomly selected cases from 10 patients, to reach reasonable agreement in terms of measuring and reporting the imaging findings of the lymph nodes. They meticulously examined and documented the characteristics of the primary tumor, including its size. For metastatic lymph nodes, we gathered information from CT images, including the size of lymph nodes in terms of their long and short axes, the ratio of long to short axes, enhancement patterns, adjacent layer invasion, lymph node aggregation, and the side of lymph node metastasis. On ultrasound images, we assessed factors such as the presence of extracapsular invasion, lymph node echogenicity, hilus and medulla structure and lymph node blood flow of LN. The masses resected by surgery with routine paraffin sectioning and HE-staining were independently reviewed by two pathologists with more than 5 years of experience. Pathologists examine pathological sections of tumor tissues, including the microscopic observation and analysis of the morphological, structural, and histological features of tumor cells, to determine the tumor’s type and degree of differentiation.

CT scans were performed by 64-slice Optima CT660 (GE Healthcare, USA), 256-slice GE Revolution CT (GE Healthcare, USA), 64-slice Siemens Somatom Sensation CT (Siemens Medical Solutions, Forchheim, Germany), or 320-slice multidetector CT (Toshiba Medicals, Tokyo, Japan). Patients underwent the CT scans with slice thicknesses of 5mm and tube voltages of 120kV. Iohexol (350mg/ml) was administered intravenously (1ml/kg) using a flow of 4 ml/s. Ultrasound scans were performed on ultrasound scanners Philips Epiq 7, Philips iU22 or Philips EnVisor with 7 to 12 MHz linear probe.

2.3 Cohort definition and variables

These 272 patients were randomly divided into the training and validation cohorts in a 3:1 ratio. Based on the 8th AJCC, we collected the following potential predictors as study variables, which include baseline demographics at diagnosis (sex, age, height, weight, BMI), tumor features (type, size, degree of differentiation), radiological features of LNs (diameters of long and short axis, long-to-short axis ratio (L/S ratio), enhancement pattern, adjacent planes with infiltration, aggregation of LNs, unilateral or bilateral LNs), features in ultrasound images features(extranodal extension, echogenicity, hilus and medulla of LN, blood flow of LNs). LNs aggregation was defined as presence or clustering of more than three LNs in succession. Infiltration of adjacent plane was defined as indistinct LNs’ margins or perinodal soft tissue infiltration.

2.4 Statistical analysis

Statistical analysis of the data was performed using SPSS 26.0 software and R software (version 4.1.2, https://www.r-project.org/). First, the caret package was used to randomly divide the patients into the training cohort and validation cohort conforming to the ratio of 3:1. Second, SPSS was used to analyze the correlation and difference of different factors in the training cohort and validation cohort. The Mann–Whitney U test was used to compare continuous data, categorical data were tested by the chi-square test or Fisher’s exact test, Third, multivariate logistic regression models were developed to study effects of predictor variables. For the construction of the multivariable model, we used a forward model building approach. Fourth, we built a nomogram model by using the nomogram package. The nomogram was drawn step-by-step according to the method of Zhang et al (6). Subsequently, Model’s performance was evaluated by AUC and C-index using survival package and pROC package, and we also assessed discrimination, calibration, and clinical utility value by calibration curves and clinical decision curves using lrm package and rmda package. Two-sided P values < 0.05 were considered statistically significant.

3 Results

3.1 Clinical characteristics

The study was conducted as shown in the low chart (Figure 1). A total of 272 patients who met the criteria [mean age (SD), 62.99(4.89) years; 256 males and 16 females] were included in this study and their clinical, medical imaging and pathological data were collected for analysis. According to the postoperative pathological results of cervical LNs, the patients were divided into the two following groups: the nonmetastatic group (n=135) and the metastatic group (n=137). Comparing the data of the two groups (Table 1). There was a significant difference between the two groups in age (64.27 ± 9.04 vs. 62.20 ± 9.15, P=0.048), and no significant difference in sex(P=0.053), BMI(P=0.408), height(P=0.718), weight(P=0.396), and tumor type (P=0.113).

Figure 1

Table 1

VariablesNon-metastasisMetastasisp
(n=135)(n=137)
Age(year)64.27 ± 9.0462.20 ± 9.150.048*
Gender0.530
Female97
Male126130
BMI(kg/m2)23.04 ± 2.7922.82 ± 3.400.408
Height(cm)169.27 ± 2.39169.46 ± 6.760.718
Weight(kg)62.22 ± 10.2065.75 ± 12.010.396
Tumor type0.113
Larynx carcinoma6864
Oral carcinoma2323
Oral floor carcinoma103
Hypopharyngeal carcinoma3447

Baseline demographics and tumor type of the HNSCC patients.

Data are presented as mean ± standard deviation or n. *, statistical significance (P<0.05).

3.2 Grouping and comparison of the training and validation cohorts

The 272 patients were randomly divided into the training cohort(n=204) and validation cohort(n=68), and no significant differences were found between the two groups (P > 0.05), indicating the consistency of these two groups (Table 2).

Table 2

VariablesTraining cohortValidation cohortP
(n=204)(n=68)
Group0.73
Non-metastatic group10332
Metastatic group10136
Age(year)63.38 ± 8.5862.78 ± 10.670.54
Gender0.77
Female133
Male19165
BMI (kg/m2)22.79 ± 3.0922.37 ± 3.140.13
Height(cm)169.09 ± 6.37170.18 ± 7.120.19
Weight (Kg)62.22 ± 10.1267.96 ± 11.840.07
Tumor type0.34
Larynx carcinoma10131
Oral carcinoma3412
Oral floor carcinoma76
Hypopharyngeal carcinoma6219
Diameter of long axis(mm)12.51 ± 8.2814.95 ± 9.620.08
Diameter of short axis (mm)8.00 ± 6.259.77 ± 7.390.08
L/S ratio1.76 ± 0.481.71 ± 0.430.61
Adjacent plane with infiltration0.14
With94
Without19564
Aggregation of LNs0.12
With7533
Without12935
Enhancement pattern0.84
Uneven/circular3612
Obvious63
Mild/without16253
Side0.75
Unilateral15048
Bilateral5420
Capsule0.55
Clear/sharp16452
Extension4016
Medulla of LNs0.39
Clear12838
Unclear/absent7630
Hilus of LNs0.40
Clear15046
Unclear/absent5422
Blood flow0.24
Rich17152
Without3316
Echogenicity0.75
Homogeneous16148
Inhomogeneous4320
Differentiation degree0.96
Well5419
Moderately12139
Poor2910
Longitudinal diameter of tumor(cm)2.71 ± 1.192.87 ± 1.440.83

Comparison of the training and validation cohorts.

Data are presented as mean ± standard deviation or n.

Furthermore, patients in the training cohort were divided into the non-metastatic group (n=103) and the metastatic group (n=101) according to pathology. Comparing the data of these two groups, the following variables as follows showed statistically significant differences between these two groups (P < 0.001): the diameters of the long and short axis of the LNs, the L/S ratio, adjacent planes with infiltration, aggregation of LNs, enhancement pattern, extranodal extension, ultrasound echogenicity, hilus and medulla of LNs, blood flow, differentiation degree and longitudinal diameter of primary tumor (Table 3).

Table 3

VariablesTraining cohortValidation cohort
Non-metastatic group(n=103)Metastatic group(n=101)PNon-metastatic group(n=32)Metastatic group (n=36)P
Age64.39 ± 8.5562.46 ± 8.570.05863.90 ± 10.5861.78 ± 10.800.491
Gender0.8900.598
Female7621
Male96953035
BMI (kg/m2)22.79 ± 2.7622.77 ± 3.410.78423.83 ± 2.7922.95 ± 3.400.221
Height(cm)169.38 ± 6.13168.80 ± 6.620.685168.91 ± 7.75171.31 ± 6.910.193
Weight (Kg)62.22 ± 10.12169.09 ± 6.370.41668.23 ± 10.4567.72 ± 13.100.676
Diameter of LNs<0.001*<0.001*
Diameter of long axis(mm)8.44 ± 4.8416.65 ± 9.0110.01 ± 5.7819.32 ± 10.28
Diameter of short axis (mm)4.57 ± 3.3211.50 ± 6.605.27 ± 3.4813.76 ± 7.67
L/S ratio1.98 ± 0.361.53 ± 0.481.99 ± 0.281.46 ± 0.38
Adjacent planes of infiltration0.003*0.616
With1813
Without102933133
Aggregation of LNs<0.001*0.009*
With23521023
Without80492213
Enhancement<0.001*<0.001*
Uneven/circular234012
Obvious1512
Mild/without100623122
Side0.3420.287
Unilateral79712523
Bilateral2430713
Capsule<0.001*<0.001*
Complete96683121
Extension733115
Medulla0.001*<0.001*
Clear91372711
Unclear/absent1264525
Hilus of LNs0.002*0.002*
Clear94562818
Unclear/absent945418
Blood flow0.004*0.002*
Rich231214
Poor101703022
Ultrasound echo<0.001*0.001*
Homogeneous97642929
Inhomogeneous63733
Differentiation<0.001*0.383
Well4141118
Moderately55551821
Poor7737
Diameter of tumor(cm)2.46 ± 1.092.95 ± 1.240.002*2.40 ± 1.023.29 ± 1.620.012*

Statistical analysis of the training and validation cohorts.

Data are presented as mean ± standard deviation or n. *, statistical significance (P<0.05).

3.3 Single- and multifactor regression analysis

Differentiation degree and diameter of primary tumor, L/S ratio of LNs, enhancement pattern, aggregation of LNs, unilateral/bilateral, ultrasound echogenicity, hilus and medulla of LNs, blood flow and extranodal extension were included in the single regression analysis (Table 4). In the training cohort, unilateral/bilateral [95% confidence interval (CI): 0.744~2.599, P=0.301] was not identified as a risk factor for LNs metastasis, and the differences in other variables were statistically significant. We included them in the multifactorial regression model (Table 4), and the results showed that the differentiation degree of the tumor (OR:2.871, 95%CI: 1.224~6.735, P=0.015), L/S ratio of LNs (OR: 0.064, 95%CI:0.019~0.217, P<0.001), inhomogeneous enhancement (OR: 4.966, 95%CI: 1.476~16.715, P=0.010), LNs aggregation (OR: 3.860, 95%CI: 1.373~10.849, P=0.010), inhomogeneous echo (OR: 5.591, 95%CI: 1.337~23.389, P=0.018), rich blood flow (OR: 12.927, 95%CI: 1.952~85.632, P=0.008), and unclear/absent medulla (OR: 10.516, 95%CI: 2.514~43.989, P=0.001) were risk factors for cervical LNs metastasis.

Table 4

VariablesUnivariateMultivariate
OR (95%CI)POR (95%CI)P
Differentiation3.283 (1.962~5.495)<0.001*2.871 (1.224~6.735)0.015*
L/S ratio0.081 (0.036~0.185)<0.001*0.064 (0.019~0.217)<0.001*
Uneven/circular enhancement7.518 (2.971~19.024)<0.001*4.966 (1.476~16.715)0.010*
Aggregation of LNs3.691 (2.014~6.767)<0.001*3.860 (1.373~10.849)0.010*
Unilateral/bilateral1.391 (0.744~2.599)0.301
Diameter of long axis1.440 (1.122~1.849)0.004*1.275 (0.805~2.019)0.301
Inhomogeneous echo9.346 (3.730~23.418)<0.001*5.591 (1.337~23.389)0.018*
Unclear/absent medulla13.117 (6.351~27.094)<0.001*10.516 (2.514~43.989)0.001*
Unclear/absent hilus8.206 (3.726~18.076)<0.001*0.943 (0.182~4.892)0.945
Rich blood flow22.143 (5.131~95.553)<0.001*12.927 (1.952~85.632)0.008*
Capsule extension7.963 (3.161~20.059)<0.001*1.367 (0.271~6.900)0.705

Univariate and multivariate analysis of the training cohort.

*, statistical significance (P<0.05).

3.4 Nomogram development and assessment of consistency, calibration, and discrimination

We constructed a nomogram according to the variables screened. The total score could be obtained by summing the individual scores based on the nomogram which corresponds to the probability of cervical LNs metastasis in HNSCC patients (Figure 2), it also showed an example of using this to predict the probability of cervical LNs metastasis. The receiver operating characteristic (ROC) curve was plotted with area under curve (AUC) of 0.953 and a cut-off value of 0.478 was determined to differentiate between high and low risks of cervical LNs metastasis (Figure 3A). Risk stratification according to the nomogram shows that patients with a score > 80 are high-risk groups for cervical lymph node metastases, a score ≤ 80 are the low-risk groups for cervical lymph node metastasis. Corrected with the calibration curve (Figure 3B), the mean absolute error between the predicted and actual values was 0.021, which showed a good degree of compliance.

Figure 2

Figure 3

3.5 Validation cohort for evaluating model differentiation

The nomogram was validated with the validation cohort, and the ROC curve was plotted with an AUC of 0.938 (Figure 3C). Then, corrected with the calibration plot (Figure 3D), the mean absolute error between the predicted and actual values was 0.024. Predictive ROC analysis of individual metrics in the nomogram for lymph node metastasis is shown in Figure 4. The predicted risk was close to the actual risk, indicating good model calibration. After 400 bootstrap self-sampling internal tests, the C-index of the model was 0.910, indicating a good agreement between the predicted and the actual situation. The clinical decision curve was plotted (Figure 5). The nomogram predicted a net benefit for decision-making across a range of threshold probabilities between 0.25 and 0.9 for the probability of cervical LNs metastasis, with a good performance.

Figure 4

Figure 5

4 Discussion

HNSCC is the sixth most common tumor worldwide which is prone to LNs metastasis, there is a significant correlation between cervical LNs metastasis and poor prognosis (7, 8). Typical therapies for HNSCC are surgery and chemoradiotherapy, which have serious side effects. Surgery may cause significant tissue distortion and dysfunction. Chemoradiotherapy may also cause severe dysfunction, such as oral dryness, tissue necrosis, atrophy or fibrosis of bones and soft tissues, resulting in difficulty swallowing (9). Currently, the treatment modalities depends on the staging, zonation of the LNs and involvement of surrounding structures (10). Even in small aggressive primary tumors, there may be occult metastases in the cervical LNs (11). Therefore, we incorporated imaging features of the LNs and clinical characteristics of patients to construct a model to accurately predict LNs staging, which has reference value for clinicians in making neck management decision before initiating a definitive treatment.

Imaging can more accurately identify asymptomatic metastatic LNs and provide a morphologic assessment of such LNs, thus, it’s becoming an effective diagnostic tool for determining treatment strategies (12). Ultrasound and CT are the traditional imaging techniques commonly used for the clinical detection of cervical LNs, because they are easily accessible and non-invasive, but both also have limitations. CT is the often first-choice imaging technique for evaluating cervical LNs and as an objective technique, it’s less subjective. In addition, the use of standard iodine contrast in contrast-enhanced CT scans allows a more accurate evaluation of tumor morphology, including size, contour, margins, and internal structure. Ultrasound is a valuable diagnostic tool in tumor staging and follow-up. The advantages of ultrasound include wide availability, good tolerability, and no radiation exposure, with the added advantage of detecting tumor invasion of blood vessels, and can be used in conjunction with fine-needle aspiration cytology (FNAC) examination (13, 14). Studies have shown that ultrasound provides information on LNs’ location, medullary features, internal echoes, and blood flow distribution. Furthermore, it can be repeated but more subjective than CT scanning (15). Based on the advantages and disadvantages of these two imaging techniques, our study evaluated the accuracy of contrast-enhanced CT and ultrasound in diagnosing cervical LNs metastases in patients with HNSCC. A recent meta-analysis showed that the sensitivity of CT for diagnosing cervical LNs was 77-81% and ultrasound was 76-93% (16). The difference in diagnostic efficacy between the two methods was not statistically significant, neither could be used to accurately detect occult cervical LNs metastasis due to the low negative predictive value of both. Our findings are consistent with previously published results, the sensitivity of neither CT nor ultrasound is 100%, the rate of false-negative results of both methods is relatively high: 21.9% for CT and 17.5% for ultrasound, suggesting that they are not sufficiently reliable methods to exclude the presence of metastatic cervical LNs. So, in these patients especially in stage T1 and T2 tumors, we cannot decide not to perform elective neck dissection based on a “negative” preoperative staging result and should still follow treatment guidelines in making treatment decisions.

Multivariate analysis can obtain coefficients for risk factors and calculate specific risk values through model formulas, but it is difficult to incorporate the predicted values of these indicators (17, 18). Researchers are increasingly interested in nomogram, which is an intuitive graph based on logistic regression results (19). For clinical applications, one of the main advantages of nomograms is the ability to estimate individualized risks based on patients’ and disease’s characteristics. Our study constructed and validated a prediction model based on imaging and pathological parameters for predicting the probability of metastasis in cervical LNs in patients with HNSCC. Poorly differentiated tumors are widely considered to be highly malignant, and previous studies have shown that tumor differentiation is an important prognostic factor for HNSCC and an independent predictor of cervical LNs metastasis (2022). Consistent with the results of previous studies, our findings suggest that the degree of tumor differentiation is significantly associated with the risk of cervical LNs metastasis. The number, size and anatomical location of LNs were the only criteria used in the clinical staging system (23, 24). LNs’ size is the most discussed criterion for assessing cervical LNs metastasis, and there is a wide variation in the size of reactive and metastatic LNs; however, the challenge is the trade-off between sensitivity and specificity, which increases with decreasing size. Although various criteria have been proposed over the past few decades, however, none is perfect (25).

Therefore, investigators incorporated LNs shape into the imaging criteria to improve the sensitivity of diagnosis (26). Normal and reactive LNs tend to be oval or reniform. In contrast, metastatic LNs are round or spherical. The L/S ratio of LNs has been suggested to quantify LNs shape and predict LNs metastasis (27). Based on these studies, we evaluated the L/S ratio on CT images, it showed that the ratio was negatively correlated with cervical LNs metastasis. LNs aggregation was defined as three or more adjacent LNs. The combination of LNs aggregation and size criteria can improve the sensitivity of the diagnosis without compromising specificity. In our study, LNs aggregation was an independent risk factor for cervical LNs metastasis. In addition, abnormalities in the internal structure of LNs are the most reliable imaging finding in patients with HNSCC, regardless of node size (28, 29). When the tumor spreads to the cervical LNs, the internal structure of the LNs is disrupted, and the tumor cell areas, keratin pooling areas, necrotic areas and residual normal LNs tissues are diffusely distributed. Based on the altered internal structure, the heterogeneity of metastatic LNs leads to an inhomogeneous enhancement or circular enhancement on CT imaging, and inhomogeneous echo, loss of hilus and medulla on ultrasonic imaging (30, 31). Our results also showed the non-homogeneity of metastatic LNs in both CT and ultrasound imaging features. Furthermore, although studies have shown that the loss of hilus is considered as the diagnostic criteria for metastatic LNs on ultrasound (32), however, in the early stages of tumors, the hilus of LNs can still be observed because they have not been sufficiently disrupted (31). In our study, the absence of hilus differed between these two groups but not an independent risk factor for LNs metastasis, we speculate that it is because of the larger proportion of patients with T2 stage in our study population. Normal and reactive LNs usually show either a lymphatic portal vascular distribution or no vascular distribution. Tumors can infiltrate LNs and produce tumor angiogenic factors, leading to angiogenesis and peripheral vascular recruitment, lymphatic portal vessels are preserved until they are destroyed by tumor cells later, so mixed vascular distribution is seen in malignant LNs (33). Ultrasound assessment of the vascular pattern of cervical LNs has been reported to be highly reliable, with a reproducibility of 85% (34). In this study, blood flow characteristics of cervical LNs were evaluated, and the presence of irregular or abundant blood flow in LNs is a valid indicator of LNs metastasis. Individual imaging parameters may be nonspecific, but when appropriately combined and applied, they can improve the accuracy of LNs disease assessment. Our study combined imaging features of ultrasound and postcontrast CT with the degree of pathological tumor differentiation, making the nomogram become more accurate for LNs risk stratification. In conclusion, we believe that this nomogram to predict the risk of LNs metastasis has a theoretical and clinical basis, also can provide reference value and auxiliary diagnosis for clinicians in making neck management decisions of HNSCC patients.

Our study also has some limitations. The size of the tumor has been described in the literature as one of the most important factors affecting LNs metastasis (35). In this study, we observed that tumor length was not an independent risk factor for cervical LNs metastasis, which may be related to the high number of patients with early stage in our case. Early-stage tumors, due to their smaller size, may lack the statistical power to demonstrate independence in relation to LN metastasis. Thus, our conclusions might be limited by a selection bias towards early-stage cases. Future research should include a more diverse range of tumor stages to assess the relationship more accurately between tumor size and LN metastasis. In addition, in both CT and ultrasound imaging parameters, we did not find the LNs border is an independent risk factor for predicting cervical LNs metastasis, but earlier reports showed that metastatic LNs tend to have well-defined sharp borders. In contrast, benign LNs usually show indistinct borders (36). Such sharp borders are attributed to increased acoustic impedance differences between the metastatic node and the surrounding tissue due to tumor infiltration (31). Advanced metastatic LNs may exhibit ill-defined borders, indicating the presence of extranodal extension (ENE). Studies have shown that ENE, a prognostic factor for HNSCC, significantly reduces the 5-year survival of patients and is associated with higher rates of local recurrence and distant metastasis (37). Thus, we should expand the sample size to analyze the role of ENE in LNs metastasis.

Due to limited research resources, this study was a small-scale single-center retrospective study, selection bias and subjective differences in the review of images by ultrasound physicians and radiologists could not be completely excluded, the external validity of their results may be limited. While our study provides valuable insights and contributes to the existing body of knowledge, we acknowledge that its single-center nature constitutes a limitation in terms of broader applicability. Future multicenter studies are recommended to validate our findings across diverse populations and healthcare settings, thus enhancing the external validity and impact of the research. Moreover, we only performed a self-test with 400 resamples of the same population and validated data from the same institution, which may lead to overfitting of the model, so independent data should be selected for external validation at other institutions. Validation of multicenter or even prospective data with large sample sizes would better improve the predictive efficacy of this nomogram and provide a higher level of evidence for clinical application.

5 Conclusion

We established and verified the prediction model for cervical lymph node metastasis in HNSCC patients with good performance and can provide auxiliary diagnosis for clinicians.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by Ethics Committee of the affiliated hospital of Qingdao University. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

XHC and LZ contributed to conception and design of the study. XHC and LZ organized the database. XHC performed the statistical analysis. LZ wrote the first draft of the manuscript. YT, BL and HJL wrote sections of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.

Funding

The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1

    ChowLQM. Head and neck cancer. N Engl J Med (2020) 382:6072. doi: 10.1056/NEJMra1715715

  • 2

    HorváthAPrekoppPPolonyG. Accuracy of the preoperative diagnostic workup in patients with head and neck cancers undergoing neck dissection in terms of nodal metastases. Eur Arch Otorhinolaryngol (2021) 278:2041–6. doi: 10.1007/s00405-020-06324-w

  • 3

    JiangQTangALongS. Development and validation of a nomogram to predict the risk of occult cervical lymph node metastases in cN0 squamous cell carcinoma of the tongue. Br J Oral Maxillofac Surg (2019) 57:1092–7. doi: 10.1016/j.bjoms.2019.09.024

  • 4

    ChenLYWengWBWangW. Analyses of high-risk factors for cervical lymph node metastasis in laryngeal squamous cell carcinoma and establishment of nomogram prediction model. Ear Nose Throat J (2021) 100:657–62. doi: 10.1177/0145561320901613

  • 5

    YuanXHCYZhangM. A nomogram for predicting occult lymph node metastasis in early hypopharyngeal cancer with cN0. Eur Arch Otorhinolaryngol (2021) 278:3515–22. doi: 10.1007/s00405-021-06648-1

  • 6

    ZhangZGRBWKM. Nomogram for survival analysis in the presence of competing risks. Ann Transl Med (2017) 5:403. doi: 10.21037/atm.2017.07.27

  • 7

    BrayFFerlayJSoerjomataramI. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin (2018) 68:394424. doi: 10.3322/caac.21492

  • 8

    FerlayJColombetMSoerjomataramI. Estimating the global cancer incidence and mortality in 2018: GLOBOCAN sources and methods. Int J Cancer (2019) 144:1941–53. doi: 10.1002/ijc.31937

  • 9

    BhandariSSoniBWBahlAGhoshalS. Radiotherapy-induced oral morbidities in head and neck cancer patients. Spec Care Dentist (2020) 40:238–50. doi: 10.1111/scd.12469

  • 10

    KlausnerGTroussierIBlaisE. Neck management in head and neck squamous cell carcinomas: where do we stand? Med Oncol (2019) 36:40. doi: 10.1007/s12032-019-1265-1

  • 11

    DriessenDAJJDijkemaTWeijsWLJTakesRPPeggeSAHZámecnikPet al. Novel diagnostic approaches for assessment of the clinically negative neck in head and neck cancer patients. Front Oncol (2021) 10:637513. doi: 10.3389/fonc.2020.637513

  • 12

    JoynerDRizviTKaleliogluTJamesonMJMukherjeeS. Lymph node dissection: principles and postoperative imaging. Neuroimaging Clin N Am (2022) 32:1936. doi: 10.1016/j.nic.2021.09.001

  • 13

    AngelicoGSantoroAInzaniFSpadolaSFiorentinoVCianfriniFet al. Ultrasound-guided FNA cytology of groin lymph nodes improves the management of squamous cell carcinoma of the vulva: Results from a comparative cytohistological study. Cancer Cytopathol (2019) 127:514–20. doi: 10.1002/cncy.22154

  • 14

    HallurNSatharRSiddiquaAZakaullahSKothariC. Evaluation of metastatic lymph nodes in oral squamous cell carcinoma: A comparative study of clinical, FNAC, ultra sonography and computed tomography with post operative histopathology. Indian J Otolaryngol Head Neck Surg (2022) 74:5921–6. doi: 10.1007/s12070-021-02495-y

  • 15

    KHRCHD. Chapter 2 squamous cell carcinoma of the head and neck-imaging evaluation of regional lymph nodes and implications for management. Semin Ultrasound CT MR (2017) 38:466–78. doi: 10.1053/j.sult.2017.05.003

  • 16

    SunJLiBJLC. Computed tomography versus magnetic resonance imaging for diagnosing cervical lymph node metastasis of head and neck cancer: a systematic review and meta-analysis. Onco Targets Ther (2015) 8:1291–313. doi: 10.2147/OTT.S73924

  • 17

    MC-ECLJRTA. Clinical and histopathological risk factors for distant metastasis in head and neck cancer patients. Acta Otorhinolaryngol Ital (2021) 41:617. doi: 10.14639/0392-100X-N0879

  • 18

    RivièreDManciniJSantiniL. Lymph-node metastasis following total laryngectomy and total pharyngolaryngectomy for laryngeal and hypopharyngeal squamous cell carcinoma: Frequency, distribution and risk factors. Eur Ann Otorhinolaryngol Head Neck Dis (2018) 135:163–6. doi: 10.1016/j.anorl.2017.11.008

  • 19

    BVPMGJSJ. Nomograms in oncology: more than meets the eye[J. Lancet Oncol (2015) 16:173–80. doi: 10.1016/S1470-2045(14)71116-7

  • 20

    KMFMFOWSIR. Predictors for distant metastasis in head and neck cancer, with emphasis on age. Eur Arch Otorhinolaryngol (2021) 278:181–90. doi: 10.1007/s00405-020-06118-0

  • 21

    KMGEMAANS. Predictors of contralateral-bilateral nodal disease in oropharyngeal cancer: A National Cancer Data Base Study. Head Neck (2018) 40:338–48. doi: 10.1002/hed.24964

  • 22

    PanareseIAquinoGRonchiALongoFMontellaMCozzolinoIet al. Oral and Oropharyngeal squamous cell carcinoma: prognostic and predictive parameters in the etiopathogenetic route. Expert Rev Anticancer Ther (2019) 19:105–19. doi: 10.1080/14737140.2019.1561288

  • 23

    AMBGFLBES. The Eighth Edition AJCC Cancer Staging Manual: Continuing to build a bridge from a population-based to a more “personalized” approach to cancer staging. CA Cancer J Clin (2017) 67:93–9. doi: 10.3322/caac.21388

  • 24

    LWMPSG‘SullivanBO. Head and Neck cancers-major changes in the American Joint Committee on cancer eighth edition cancer staging manual. CA Cancer J Clin (2017) 67:122–37. doi: 10.3322/caac.21389

  • 25

    ZhangYYuDYangQLiW. Diagnostic efficacy of physical examination, preoperative ultrasound, and/or computed tomography in detecting lymph node metastasis: A single-center retrospective analysis of patients with squamous cell carcinoma of the head and neck. Oral Surg Oral Med Oral Pathol Oral Radiol (2022) 134:386–96. doi: 10.1016/j.oooo.2022.05.002

  • 26

    TreutleinCStollbergAScherlCAgaimyAEllmannSIroHet al. Diagnostic value of 3D dynamic contrast-enhanced magnetic resonance imaging in lymph node metastases of head and neck tumors: a correlation study with histology. Acta Radiol Open (2020) 9:205846012095196. doi: 10.1177/2058460120951966

  • 27

    WangJLuXZhengXXiaCLiP. Clinical value of preoperative ultrasound signs in evaluating axillary lymph node status in triple-negative breast cancer. J Oncol (2022) 2022:17. doi: 10.1155/2022/2590647

  • 28

    BloomerZWSnitchlerANShakirMKMHoangTD. Metastatic squamous cell carcinoma in a cervical lymph node. Clin Case Rep (2020) 8:2322–3. doi: 10.1002/ccr3.3326

  • 29

    MahmoodSMairMFagiryRAhmedMMMenonIIbrahimNet al. Diagnostic efficacy of combined CT and MRI in detecting nodal metastasis in patients with oral cancer. Oral Surg Oral Med Oral Pathol Oral Radiol (2022) 133:343–8. doi: 10.1016/j.oooo.2021.09.003

  • 30

    ForghaniRYuELeventalM. Imaging evaluation of lymphadenopathy and patterns of lymph node spread in head and neck cancer. Expert Rev Anticancer Ther (2015) 15:207–24. doi: 10.1586/14737140.2015.978862

  • 31

    PrativadiRDahiyaNKamayaA. Chapter 5 ultrasound characteristics of benign vs Malignant cervical lymph nodes. Semin Ultrasound CT MR (2017) 38:506–15. doi: 10.1053/j.sult.2017.05.005

  • 32

    WagnerJMAllemanAM. Ultrasonography of cervical lymph nodes. Radiol Clin North Am (2019) 57:485500. doi: 10.1016/j.rcl.2019.01.005

  • 33

    ChammasMCMacedoTAALoVWGomesACJulianoACerriGG. Predicting Malignant neck lymphadenopathy using color duplex sonography based on multivariate analysis. J Clin Ultrasound (2016) 44:587–94. doi: 10.1002/jcu.22380

  • 34

    De Koekkoek-DollPKRobertiSVan Den BrekelMWMaasMSmitLBeets-TanRet al. Value of assessing peripheral vascularization with micro-flow imaging, resistive index and absent hilum sign as predictor for Malignancy in lymph nodes in head and neck squamous cell carcinoma. Cancers (2021) 13:5071. doi: 10.3390/cancers13205071

  • 35

    OkabeRUekiYOhashiRTakeuchiMHashimotoSTakahashiTet al. Predicting cervical lymph node metastasis following endoscopic surgery in superficial head and neck carcinoma. Front Surg (2022) 8:813260. doi: 10.3389/fsurg.2021.813260

  • 36

    JayapalNRamSMMurthyVBasheerSShamsuddinSKhanA. Differentiation between benign and metastatic cervical lymph nodes using ultrasound. J Pharm Bioallied Sci (2019) 11:338. doi: 10.4103/JPBS.JPBS_26_19

  • 37

    MermodMTolstonogGSimonC. Extracapsular spread in head and neck squamous cell carcinoma: A systematic review and meta-analysis. Oral Oncol (2016) 62:6071. doi: 10.1016/j.oraloncology.2016.10.003

Summary

Keywords

head and neck squamous cell carcinoma (HNSCC), nomogram, pathologic, lymph node, image

Citation

Chen X, Zhang L, Lu H, Tan Y and Li B (2024) Development and validation of a nomogram to predict cervical lymph node metastasis in head and neck squamous cell carcinoma. Front. Oncol. 13:1174457. doi: 10.3389/fonc.2023.1174457

Received

26 February 2023

Accepted

12 December 2023

Published

04 January 2024

Volume

13 - 2023

Edited by

Olgun Elicin, University Hospital Bern, Switzerland

Reviewed by

Hao Chi, The Affiliated Hospital of Southwest Medical University, China

Wilfredo Alejandro González-Arriagada, University of the Andes, Chile

Updates

Copyright

*Correspondence: Haijun Lu,

†These authors have contributed equally to this work and share first authorship

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics