Abstract
BACKGROUND: Identifying occult central lymph node metastasis (CLNM) is essential for guiding prophylactic lymph node dissection (PLND) in patients with cN0 stage papillary thyroid microcarcinoma (PTMC). This study aimed to identify molecular prognostic biomarkers associated with PTMC and develop a clinical-molecular prediction model for CLNM. METHODS: Differentially expressed genes (DEGs) in PTMC were identified through bioinformatics analysis of the TCGA database. Prognostic DEGs were selected using Cox and LASSO regression analyses, and a risk-scoring model was constructed based on these genes. The prognostic value of the model was validated using Kaplan-Meier survival analysis and ROC curves. DEG expression levels were compared between patients with CLNM and those without (NCLNM). Clinical data and surgical specimens were collected from 404 patients with cN0 stage PTMC treated at the First Affiliated Hospital of Ningbo University in 2022. The cohort was randomly divided into a derivation cohort (n = 323) and a validation cohort (n = 81). DEG expression was quantified using RT-qPCR. Univariate and multivariate logistic regression analyses were conducted in the derivation cohort to identify predictors of CLNM and develop a predictive model. The model's performance was evaluated using the Hosmer-Lemeshow test, ROC curves, calibration curves, and decision curve analysis (DCA). RESULTS: In the TCGA database, FN1, MT-1 F, and TFF3 were identified as prognostic biomarkers. Risk scores based on these genes achieved AUCs of 0.789 (5 years) and 0.674 (10 years) for predicting disease-free survival. Furthermore, FN1, MT-1 F, and TFF3 expression levels were significantly higher in the CLNM group compared to the NCLNM group. Among the 404 PTMC patients, the incidence of CLNM was 42.6% (n = 172). RT-qPCR analysis demonstrated significantly elevated expression of FN1 in both PTMC tissues compared to normal tissues and in the CLNM group relative to the NCLNM group, while MT-1 F and TFF3 exhibited markedly reduced expression levels. In the derivation cohort, FN1, MT-1 F, TFF3, tumor size ≥5 mm, calcification, multifocality, and extrathyroidal extension were independent predictors of CLNM. The prediction model based on these factors showed AUCs of 0.736 (derivation cohort) and 0.813 (validation cohort). Moreover, calibration curves, the Hosmer-Lemeshow test (χ² = 2.411, P = 0.966), and DCA confirmed the model's robust performance and clinical utility. CONCLUSION: FN1, MT-1 F, and TFF3 are valuable prognostic biomarkers for PTMC. The clinical-molecular prediction model incorporating these genes provides a basis for personalized PLND decision-making in cN0 stage PTMC patients. TRIAL REGISTRATION NUMBER: Not applicable.