Mesenchymal-like pancreatic cancer cells harbor specific genomic alterations more frequently than their epithelial-like counterparts

Establishment and Molecular Characterization of Two Patient-Derived Pancreatic Ductal Adenocarcinoma Cell Lines as Preclinical Models for Treatment Response

The prognosis of pancreatic ductal adenocarcinoma (PDAC) is exceedingly poor. Although surgical resection is the only curative treatment option, multimodal treatment is of the utmost importance, as only about 20% of tumors are primarily resectable at the time of diagnosis. The choice of chemotherapeutic treatment regimens involving gemcitabine and FOLFIRINOX is currently solely based on the patient's performance status, but, ideally, it should be based on the tumors' individual biology. We established two novel patient-derived primary cell lines from surgical PDAC specimens. LuPanc-1 and LuPanc-2 were derived from a pT3, pN1, G2 and a pT3, pN2, G3 tumor, respectively, and the clinical follow-up was fully annotated. STR-genotyping revealed a unique profile for both cell lines. The population doubling time of LuPanc-2 was substantially longer than that of LuPanc-1 (84 vs. 44 h). Both cell lines exhibited a typical epithelial morphology and expressed moderate levels of CK7 and E-cadherin. LuPanc-1, but not LuPanc-2, co-expressed E-cadherin and vimentin at the single-cell level, suggesting a mixed epithelial-mesenchymal differentiation. LuPanc-1 had a missense mutation (p.R282W) and LuPanc-2 had a frameshift deletion (p.P89X) in TP53. BRCA2 was nonsense-mutated (p.Q780*) and CREBBP was missense-mutated (p.P279R) in LuPanc-1. CDKN2A was missense-mutated (p.H83Y) in LuPanc-2. Notably, only LuPanc-2 harbored a partial or complete deletion of DPC4. LuPanc-1 cells exhibited high basal and transforming growth factor (TGF)-β1-induced migratory activity in real-time cell migration assays, while LuPanc-2 was refractory. Both LuPanc-1 and LuPanc-2 cells responded to treatment with TGF-β1 with the activation of SMAD2; however, only LuPanc-1 cells were able to induce TGF-β1 target genes, which is consistent with the absence of DPC4 in LuPanc-2 cells. Both cell lines were able to form spheres in a semi-solid medium and in cell viability assays, LuPanc-1 cells were more sensitive than LuPanc-2 cells to treatment with gemcitabine and FOLFIRINOX. In summary, both patient-derived cell lines show distinct molecular phenotypes reflecting their individual tumor biology, with a unique clinical annotation of the respective patients. These preclinical ex vivo models can be further explored for potential new treatment strategies and might help in developing personalized (targeted) therapy regimens.

Examining the efficacy of localised gemcitabine therapy for the treatment of pancreatic cancer using a hybrid agent-based model

The prognosis for pancreatic ductal adenocarcinoma (PDAC) patients has not significantly improved in the past 3 decades, highlighting the need for more effective treatment approaches. Poor patient outcomes and lack of response to therapy can be attributed, in part, to a lack of uptake of perfusion of systemically administered chemotherapeutic drugs into the tumour. Wet-spun alginate fibres loaded with the chemotherapeutic agent gemcitabine have been developed as a potential tool for overcoming the barriers in delivery of systemically administrated drugs to the PDAC tumour microenvironment by delivering high concentrations of drug to the tumour directly over an extended period. While exciting, the practicality, safety, and effectiveness of these devices in a clinical setting requires further investigation. Furthermore, an in-depth assessment of the drug-release rate from these devices needs to be undertaken to determine whether an optimal release profile exists. Using a hybrid computational model (agent-based model and partial differential equation system), we developed a simulation of pancreatic tumour growth and response to treatment with gemcitabine loaded alginate fibres. The model was calibrated using in vitro and in vivo data and simulated using a finite volume method discretisation. We then used the model to compare different intratumoural implantation protocols and gemcitabine-release rates. In our model, the primary driver of pancreatic tumour growth was the rate of tumour cell division. We were able to demonstrate that intratumoural placement of gemcitabine loaded fibres was more effective than peritumoural placement. Additionally, we quantified the efficacy of different release profiles from the implanted fibres that have not yet been tested experimentally. Altogether, the model developed here is a tool that can be used to investigate other drug delivery devices to improve the arsenal of treatments available for PDAC and other difficult-to-treat cancers in the future.

Occult polyclonality of preclinical pancreatic cancer models drives in vitro evolution

Heterogeneity is a hallmark of cancer. The advent of single-cell technologies has helped uncover heterogeneity in a high-throughput manner in different cancers across varied contexts. Here we apply single-cell sequencing technologies to reveal inherent heterogeneity in assumptively monoclonal pancreatic cancer (PDAC) cell lines and patient-derived organoids (PDOs). Our findings reveal a high degree of both genomic and transcriptomic polyclonality in monolayer PDAC cell lines, custodial variation induced by growing apparently identical cell lines in different laboratories, and transcriptomic shifts in transitioning from 2D to 3D spheroid growth models. Our findings also call into question the validity of widely available immortalized, non-transformed pancreatic lines as contemporaneous "control" lines in experiments. We confirm these findings using a variety of independent assays, including but not limited to whole exome sequencing, single-cell copy number variation sequencing (scCNVseq), single-nuclei assay for transposase-accessible chromatin with sequencing, fluorescence in-situ hybridization, and single-cell RNA sequencing (scRNAseq). We map scRNA expression data to unique genomic clones identified by orthogonally-gathered scCNVseq data of these same PDAC cell lines. Further, while PDOs are known to reflect the cognate in vivo biology of the parental tumor, we identify transcriptomic shifts during ex vivo passage that might hamper their predictive abilities over time. The impact of these findings on rigor and reproducibility of experimental data generated using established preclinical PDAC models between and across laboratories is uncertain, but a matter of concern.

Identification of the Neurokinin-1 Receptor as Targetable Stratification Factor for Drug Repurposing in Pancreatic Cancer

The SP/NK1R-complex plays an important role in tumor proliferation. Targeting of the neurokinin-1 receptor in previous studies with its antagonist aprepitant (AP) resulted in anti-tumoral effects in colorectal cancer and hepatoblastoma. However, there is still a lack of knowledge regarding its effects on pancreatic cancer. Therefore, we treated human pancreatic ductal adenocarcinoma (PDAC) cell lines (Capan-1, DanG, HuP-T3, Panc-1, and MIA PaCa-2) and their cancer stem cell-like cells (CSCs) with AP and analyzed functional effects by MTT-, colony, and sphere formation assays, respectively; moreover, we monitored downstream mechanisms by flow cytometry. NK1R inhibition resulted in dose-dependent growth reduction in both CSCs and non-CSCs without induction of apoptosis in most PDAC cell lines. More importantly, we identified striking AP dependent cell cycle arrest in all parental cells. Furthermore, gene expression and the importance of key genes in PDAC tumorigenesis were analyzed combining RT-qPCR in eight PDAC cell lines with publicly available datasets (TCGA, GEO, CCLE). Surprisingly, we found a better overall survival in patients with high NK1R levels, while at the same time, NK1R was significantly decreased in PDAC tissue compared to normal tissue. Interestingly, there is currently no differentiation between the isoforms of NK1R (truncated and full; NK1R-tr and -fl) in any of the indicated public transcriptomic records, although many publications already emphasize on important regulatory differences between the two isoforms of NK1R in many cancer entities. In conclusion, analysis of splice variants might potentially lead to a stratification of PDAC patients for NK1R-directed therapies. Furthermore, we presume PDAC patients with high expressions of NK1R-tr might benefit from treatment with AP to improve chemoresistance. Therefore, analysis of splice variants might potentially lead to a stratification of PDAC patients for NK1R-directed therapies.

Unbiased in vivo preclinical evaluation of anticancer drugs identifies effective therapy for the treatment of pancreatic adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is typically diagnosed at an advanced stage, which limits surgical options and portends a dismal prognosis. Current oncologic PDAC therapies confer marginal benefit and, thus, a significant unmet clinical need exists for new therapeutic strategies. To identify effective PDAC therapies, we leveraged a syngeneic orthotopic PDAC transplant mouse model to perform a large-scale, in vivo screen of 16 single-agent and 41 two-drug targeted therapy combinations in mice. Among 57 drug conditions screened, combined inhibition of heat shock protein (Hsp)-90 and MEK was found to produce robust suppression of tumor growth, leading to an 80% increase in the survival of PDAC-bearing mice with no significant toxicity. Mechanistically, we observed that single-agent MEK inhibition led to compensatory activation of resistance pathways, including components of the PI3K/AKT/mTOR signaling axis, which was overcome with the addition of HSP90 inhibition. The combination of HSP90(i) + MEK(i) was also active in vitro in established human PDAC cell lines and in vivo in patient-derived organoid PDAC transplant models. These findings encourage the clinical development of HSP90(i) + MEK(i) combination therapy and highlight the power of clinically relevant in vivo model systems for identifying cancer therapies.

Mesenchymal-like pancreatic cancer cells harbor specific genomic alterations more frequently than their epithelial-like counterparts

The aggressiveness of pancreatic cancer is associated with the acquisition of mesenchymal characteristics by a subset of pancreatic cancer cells. The factors driving the development of this subset are not well understood. In this study, we tested the hypothesis that acquisition of a mesenchymal phenotype occurs selectively in tumor cells that harbor specific enabling genetic alterations. We obtained whole-genome comparative genomic hybridization (CGH) measurements on pancreatic cancer cell lines that have either an epithelial-like (17 cell lines) or a mesenchymal-like (9 cell lines) phenotype in vitro. The total amounts of amplifications and deletions were equivalent between the epithelial and mesenchymal groups, but 20 genes showed a major difference between the groups in prevalence of alterations. All 20 alterations (18 deletions and 2 amplifications) were more prevalent in the mesenchymal group, confirming the advanced nature of this cellular subtype. CDKN2A was altered in more than 50% of both groups, but co-deletions in neighboring genes, and concomitant loss of gene expression, were more prevalent in the mesenchymal group, suggesting that the size of the loss around CDKN2A affects cell phenotype. Whole-genome CGH on 11 primary cancer tissues revealed that the 20 genes were altered at a higher prevalence (up to 55% of the cases for certain genes) than randomly selected sets of 20 genes, with the same direction of alteration as in the cell lines. These findings support the concept that specific genetic alterations enable phenotype plasticity and provide promising candidate genes for further research.