Background: There is still no reasonably accurate method of preoperatively predicting central lymph node metastasis (LNM), and it is essential to develop an effective evaluation model for predicting LNM in papillary thyroid carcinoma (PTC) patients. Methods: PTC samples were collected from The Cancer Genome Atlas database. Candidate genes were identified as continuously upregulated or downregulated genes in the process of N0 to N1a and N1a to N1b. The least absolute shrinkage and selection operator (LASSO) regression analysis was used to construct the predictive model for LNM. Multivariate logistic regression analysis was performed to screen the potential factors related to LNM, and a nomogram was established. The risk score of the gene signature model for predicting disease-free survival (DFS) was evaluated by Kaplan-Meier analysis. Results: A 14-gene signature was developed by LASSO regression for predicting LNM based on 69 differential expression genes (DEGs) that were continuously upregulated or downregulated in the progress of PTC. The receiver operating characteristic (ROC) curves of the 14-gene signature predicting LNM, central LNM and lateral LNM were generated. The area under the ROC (AUC) values were 0.806 [95% confidence interval (CI): 0.7608-0.8815], 0.755 (95% CI: 0.6839-0.8263) and 0.821 (95% CI: 0.7608-0.8815). The nomogram's C-index value, including the 14-gene signature and other potential risk factors, was 0.786 (95% CI: 0.7296-0.8425), and the calibration exhibited fairly good consistency with the perfect prediction. Based on the 14-gene risk score, high-risk PTC patients had a worse DFS. Conclusions: A novel 14-gene signature was developed for predicting LNM in PTC patients. The risk score also correlated with DFS in PTC patients.