BackgroundTo provide a preoperative prediction model for lymph node metastasis in pancreatic cancer patients and provide molecular information of key radiomic features.MethodsTwo cohorts comprising 151 and 54 pancreatic cancer patients were included in the analysis. Radiomic features from the tumor region of interests were extracted by using PyRadiomics software. We used a framework that incorporated 10 machine learning algorithms and generated 77 combinations to construct radiomics-based models for lymph node metastasis prediction. Weighted gene coexpression network analysis (WGCNA) was subsequently performed to determine the relationships between gene expression levels and radiomic features. Molecular pathways enrichment analysis was performed to uncover the underlying molecular features.ResultsPatients in the in-house cohort (mean age, 61.3 years ± 9.6 [SD]; 91 men [60%]) were separated into training (n = 105, 70%) and validation (n = 46, 30%) cohorts. A total of 1,239 features were extracted and subjected to machine learning algorithms. The 77 radiomic models showed moderate performance for predicting lymph node metastasis, and the combination of the StepGBM and Enet algorithms had the best performance in the training (AUC = 0.84, 95% CI = 0.77–0.91) and validation (AUC = 0.85, 95% CI = 0.73–0.98) cohorts. We determined that 15 features were core variables for lymph node metastasis. Proliferation-related processes may respond to the main molecular alterations underlying these features.ConclusionsMachine learning-based radiomics could predict the status of lymph node metastasis in pancreatic cancer, which is associated with proliferation-related alterations.