Abstract
BACKGROUND: With poor prognosis and high mortality, pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal malignancies. Standard of care therapies for PDAC have included gemcitabine for the past three decades, although resistance often develops within weeks of chemotherapy initiation through an array of possible mechanisms. METHODS: We reanalyzed publicly available RNA-seq gene expression profiles of 28 PDAC patient-derived xenograft (PDX) models before and after a 21-day gemcitabine treatment using our validated analysis pipeline to identify molecular markers of intrinsic and acquired resistance. RESULTS: Using normalized RNA-seq quantification measurements, we first identified oxidative phosphorylation and interferon alpha pathways as the two most enriched cancer hallmark gene sets in the baseline gene expression profile associated with intrinsic gemcitabine resistance and sensitivity, respectively. Furthermore, we discovered strong correlations between drug-induced expression changes in glycolysis and oxidative phosphorylation genes and response to gemcitabine, which suggests that these pathways may be associated with acquired gemcitabine resistance mechanisms. Thus, we developed prediction models using baseline gene expression profiles in those pathways and validated them in another dataset of 12 PDAC models from Novartis. We also developed prediction models based on drug-induced expression changes in genes from the Molecular Signatures Database (MSigDB)'s curated 50 cancer hallmark gene sets. Finally, pathogenic TP53 mutations correlated with treatment resistance. CONCLUSION: Our results demonstrate that concurrent upregulation of both glycolysis and oxidative phosphorylation pathways occurs in vivo in PDAC PDXs following gemcitabine treatment and that pathogenic TP53 status had association with gemcitabine resistance in these models. Our findings may elucidate the molecular basis for gemcitabine resistance and provide insights for effective drug combination in PDAC chemotherapy.