Abstract
PURPOSE: To identify the clinical and genetic variables associated with rim enhancement of pancreatic ductal adenocarcinoma (PDAC) and to develop a dynamic contrast-enhanced (DCE) MRI-based radiomics model for predicting the genetic status from next-generation sequencing (NGS). MATERIALS AND METHODS: Patients with PDAC, who underwent pretreatment pancreatic DCE-MRI between November 2019 and July 2021, were eligible in this prospective study. Two radiologists evaluated presence of rim enhancement in PDAC, a known radiological prognostic indicator, on DCE MRI. NGS was conducted for the tissue from the lesion. The Mann-Whitney U and Chi-square tests were employed to identify clinical and genetic variables associated with rim enhancement in PDAC. For continuous variables predicting rim enhancement, the cutoff value was set based on the Youden's index from the receiver operating characteristic (ROC) curve. Radiomics features were extracted from a volume-of-interest of PDAC on four DCE maps (K(trans), K(ep), V(e), and iAUC). A random forest (RF) model was constructed using 10 selected radiomics features from a pool of 392 original features. This model aimed to predict the status of significant NGS variables associated with rim enhancement. The performance of the model was validated using test set. RESULTS: A total of 55 patients (32 men; median age 71 years) were randomly assigned to the training (n = 41) and test (n = 14) sets. In the training set, KRAS, TP53, CDKN2A, and SMAD4 mutation rates were 92.3%, 61.8%, 14.5%, and 9.1%, respectively. Tumor size and KRAS variant allele frequency (VAF) differed between rim-enhancing (n = 12) and nonrim-enhancing (n = 29) PDACs with a cutoff of 17.22%. The RF model's average AUC from 10-fold cross-validation for predicting KRAS VAF status was 0.698. In the test set comprising 6 tumors with low KRAS VAF and 8 with high KRAS VAF, the RF model's AUC reached 1.000, achieving a sensitivity of 75.0%, specificity of 100% and accuracy of 87.5%. CONCLUSION: Rim enhancement of PDAC is associated with KRAS VAF derived from NGS-based genetic information. For predicting the KRAS VAF status in PDAC, a radiomics model based on DCE maps showed promising results.