Pancreatic ductal adenocarcinoma (PDAC) remains a lethal disease with a 5-year survival of 4%. A key hallmark of PDAC is extensive stromal involvement, which makes capturing precise tumor-specific molecular information difficult. Here, we have overcome this problem by applying blind source separation to a diverse collection of PDAC gene expression microarray data, which includes primary, metastatic, and normal samples. By digitally separating tumor, stroma, and normal gene expression, we have identified and validated two tumor-specific subtypes including a “basal-like” subtype which has worse outcome, and is molecularly similar to basal tumors in bladder and breast cancer. Furthermore, we define “normal” and “activated” stromal subtypes which are independently prognostic. Our results provide new insight into the molecular composition of PDAC which may be used to tailor therapies or provide decision support in a clinical setting where the choice and timing of therapies is critical.

SUMMARY We performed integrated genomic, transcriptomic and proteomic profiling of 150 pancreatic ductal adenocarcinoma (PDAC) specimens, including samples with characteristic low neoplastic cellularity. Deep whole-exome sequencing revealed recurrent somatic mutations in KRAS, TP53, CDKN2A, SMAD4, RNF43, ARID1A, TGFβR2, GNAS, RREB1 and PBRM1. KRAS wild-type tumors harbored alterations in other oncogenic drivers, including GNAS, BRAF, CTNNB1 and additional RAS pathway genes. A subset of tumors harbored multiple KRAS mutations, with some showing evidence of biallelic mutations. Protein profiling identified a favorable prognosis subset with low epithelial-mesenchymal transition and high MTOR pathway scores. Associations of non-coding RNAs with tumor-specific mRNA subtypes were also identified. Our integrated multi-platform analysis reveals a complex molecular landscape of PDAC and provides a roadmap for precision medicine.

Pancreatic cancer is the most lethal common solid malignancy. Systemic therapies are often ineffective, and predictive biomarkers to guide treatment are urgently needed. We generated a pancreatic cancer patient-derived organoid (PDO) library that recapitulates the mutational spectrum and transcriptional subtypes of primary pancreatic cancer. New driver oncogenes were nominated and transcriptomic analyses revealed unique clusters. PDOs exhibited heterogeneous responses to standard-of-care chemotherapeutics and investigational agents. In a case study manner, we found that PDO therapeutic profiles paralleled patient outcomes and that PDOs enabled longitudinal assessment of chemosensitivity and evaluation of synchronous metastases. We derived organoid-based gene expression signatures of chemosensitivity that predicted improved responses for many patients to chemotherapy in both the adjuvant and advanced disease settings. Finally, we nominated alternative treatment strategies for chemorefractory PDOs using targeted agent therapeutic profiling. We propose that combined molecular and therapeutic profiling of PDOs may predict clinical response and enable prospective therapeutic selection. New approaches to prioritize treatment strategies are urgently needed to improve survival and quality of life for patients with pancreatic cancer. Combined genomic, transcriptomic, and therapeutic profiling of PDOs can identify molecular and functional subtypes of pancreatic cancer, predict therapeutic responses, and facilitate precision medicine for patients with pancreatic cancer. .

To perform real-time whole genome sequencing (WGS) and RNA sequencing (RNASeq) of advanced pancreatic ductal adenocarcinoma (PDAC) to identify predictive mutational and transcriptional features for better treatment selection. Patients with advanced PDAC were prospectively recruited prior to first-line combination chemotherapy. Fresh tumor tissue was acquired by image-guided percutaneous core biopsy for WGS and RNASeq. Laser capture microdissection was performed for all cases. Primary endpoint was feasibility to report WGS results prior to first disease assessment CT scan at 8 weeks. The main secondary endpoint was discovery of patient subsets with predictive mutational and transcriptional signatures. Sixty-three patients underwent a tumor biopsy between December 2015 and June 2017. WGS and RNASeq were successful in 62 (98%) and 60 (95%), respectively. Genomic results were reported at a median of 35 days (range, 19-52 days) from biopsy, meeting the primary feasibility endpoint. Objective responses to first-line chemotherapy were significantly better in patients with the classical PDAC RNA subtype compared with those with the basal-like subtype ( = 0.004). The best progression-free survival was observed in those with classical subtype treated with m-FOLFIRINOX. expression in tumor measured by RNA hybridization was found to be a robust surrogate biomarker for differentiating classical and basal-like PDAC subtypes. Potentially actionable genetic alterations were found in 30% of patients. Prospective genomic profiling of advanced PDAC is feasible, and our early data indicate that chemotherapy response differs among patients with different genomic/transcriptomic subtypes. .