Abstract
Background: Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal malignancy. Next-generation sequencing (NGS) enables molecular characterization and may identify clinically actionable alterations; however, real-world multicenter data linking genomic subgroups to survival outcomes remain limited. We aimed to characterize the molecular landscape of NGS-tested PDAC in a Turkish multicenter cohort and evaluate the association of key molecular alterations, including KRAS status and KRAS variant subgroups, with survival outcomes. Methods: We conducted a multicenter retrospective cohort study including patients with pathologically diagnosed PDAC between 2017 and 2025 who underwent tumor-based NGS in routine clinical practice. Overall survival (OS) was calculated from the date of metastasis, defined as the date of diagnosis for de novo metastatic disease and the date of first documented distant recurrence for recurrent cases. Progression-free survival (PFS) was calculated from the initiation of first-line systemic therapy for metastatic disease to progression or death. Survival was estimated using the Kaplan-Meier method and compared using the log-rank test. Multivariable Cox proportional hazards models were constructed for OS and PFS using clinically relevant covariates selected a priori. Results: A total of 98 patients underwent molecular profiling, and survival analyses were performed in 92 patients with available OS/PFS data. KRAS mutations were detected in 83.7% (82/98) of patients, with predominant variants G12D (47.6%), G12V (30.5%), and G12R (12.2%). TP53 mutations were present in 59.2% (58/98) of tumors, and all tumors were microsatellite stable. Tumor mutational burden data were available for 72 patients; the median TMB was 3.83 mutations/Mb, and 15.3% of evaluable tumors had a TMB ≥ 10 mutations/Mb. Excluding KRAS, clinically actionable alterations were identified in 4.1% of patients, whereas an additional 32.7% harbored potentially actionable or investigational alterations. Median OS was 14.0 months (95% CI, 11.7-16.3), and median PFS was 6.0 months (95% CI, 4.3-7.7). In unadjusted analyses, OS and PFS did not differ significantly according to KRAS mutation status (OS, p = 0.967; PFS, p = 0.652), TP53 mutation status (OS, p = 0.404; PFS, p = 0.510), or KRAS variant subgroup (OS, p = 0.332; PFS, p = 0.194). In multivariable Cox analyses, KRAS mutation status was not independently associated with OS (aHR 1.13, 95% CI 0.56-2.28; p = 0.727) or PFS (aHR 1.09, 95% CI 0.59-2.01; p = 0.780), whereas ECOG performance status remained the strongest adverse clinical factor. Conclusions: In this multicenter real-world PDAC cohort, the molecular landscape was dominated by KRAS and TP53 alterations, whereas clinically actionable non-KRAS alterations were identified in only a minority of patients. After adjustment for major clinical covariates, KRAS mutation status was not independently associated with OS or PFS. Molecular profiling may still be useful for identifying uncommon potentially targetable alterations; however, larger clinically annotated multicenter studies are needed to better define its prognostic and treatment-directing value in routine practice.