Recent studies have offered ample insight into genome-wide expression patterns to define pancreatic ductal adenocarcinoma (PDAC) subtypes, although there remains a lack of knowledge regarding the underlying epigenomics of PDAC. Here we perform multi-parametric integrative analyses of chromatin immunoprecipitation-sequencing (ChIP-seq) on multiple histone modifications, RNA-sequencing (RNA-seq), and DNA methylation to define epigenomic landscapes for PDAC subtypes, which can predict their relative aggressiveness and survival. Moreover, we describe the state of promoters, enhancers, super-enhancers, euchromatic, and heterochromatic regions for each subtype. Further analyses indicate that the distinct epigenomic landscapes are regulated by different membrane-to-nucleus pathways. Inactivation of a basal-specific super-enhancer associated pathway reveals the existence of plasticity between subtypes. Thus, our study provides new insight into the epigenetic landscapes associated with the heterogeneity of PDAC, thereby increasing our mechanistic understanding of this disease, as well as offering potential new markers and therapeutic targets.

Preclinical models based on patient-derived xenografts have remarkable specificity in distinguishing transformed human tumor cells from non-transformed murine stromal cells computationally. We obtained 29 pancreatic ductal adenocarcinoma (PDAC) xenografts from either resectable or non-resectable patients (surgery and endoscopic ultrasound-guided fine-needle aspirate, respectively). Extensive multiomic profiling revealed two subtypes with distinct clinical outcomes. These subtypes uncovered specific alterations in DNA methylation and transcription as well as in signaling pathways involved in tumor-stromal cross-talk. The analysis of these pathways indicates therapeutic opportunities for targeting both compartments and their interactions. In particular, we show that inhibiting NPC1L1 with Ezetimibe, a clinically available drug, might be an efficient approach for treating pancreatic cancers. These findings uncover the complex and diverse interplay between PDAC tumors and the stroma and demonstrate the pivotal role of xenografts for drug discovery and relevance to PDAC.

Pancreatic ductal adenocarcinoma (PDAC) is composed of stromal, immune, and cancerous epithelial cells. Transcriptomic analysis of the epithelial compartment allows classification into different phenotypic subtypes as classical and basal‐like. However, little is known about the intra‐tumor heterogeneity particularly in the epithelial compartment. Growing evidences suggest that this phenotypic segregation is not so precise and different cancerous cell types may coexist in a single tumor. To test this hypothesis, we performed single‐cell transcriptomic analyses using combinational barcoding exclusively on epithelial cells from six different classical PDAC patients obtained by Endoscopic Ultrasound (EUS) with Fine Needle Aspiration (FNA). To purify the epithelial compartment, PDAC were grown as biopsy‐derived pancreatic cancer organoids. Single‐cell transcriptomic analysis allowed the identification of four main cell clusters present in different proportions in all tumors. Remarkably, although all these tumors were classified as classical, one cluster present in all corresponded to a basal‐like phenotype. These results reveal an unanticipated high heterogeneity of pancreatic cancers and demonstrate that basal‐like cells, which have a highly aggressive phenotype, are more widespread than expected.

Growing evidence indicates that Rab GTPases, key regulators of intracellular transport in eukaryotic cells, play an important role in cancer. We analysed the deregulation at the transcriptional level of the genes encoding Rab proteins and Rab-interacting proteins in bladder cancer pathogenesis, distinguishing between the two main progression pathways so far identified in bladder cancer: the Ta pathway characterized by a high frequency of FGFR3 mutation and the carcinoma in situ pathway where no or infrequent FGFR3 mutations have been identified. A systematic literature search identified 61 genes encoding Rab proteins and 223 genes encoding Rab-interacting proteins. Transcriptomic data were obtained for normal urothelium samples and for two independent bladder cancer data sets corresponding to 152 and 75 tumors. Gene deregulation was analysed with the SAM (significant analysis of microarray) test or the binomial test. Overall, 30 genes were down-regulated, and 13 were up-regulated in the tumor samples. Five of these deregulated genes (LEPRE1, MICAL2, RAB23, STXBP1, SYTL1) were specifically deregulated in FGFR3-non-mutated muscle-invasive tumors. No gene encoding a Rab or Rab-interacting protein was found to be specifically deregulated in FGFR3-mutated tumors. Cluster analysis showed that the RAB27 gene cluster (comprising the genes encoding RAB27 and its interacting partners) was deregulated and that this deregulation was associated with both pathways of bladder cancer pathogenesis. Finally, we found that the expression of KIF20A and ZWINT was associated with that of proliferation markers and that the expression of MLPH, MYO5B, RAB11A, RAB11FIP1, RAB20 and SYTL2 was associated with that of urothelial cell differentiation markers. This systematic analysis of Rab and Rab effector gene deregulation in bladder cancer, taking relevant tumor subgroups into account, provides insight into the possible roles of Rab proteins and their effectors in bladder cancer pathogenesis. This approach is applicable to other group of genes and types of cancer.

Chondrosarcomas are primary cancers of cartilaginous tissue with highly contrasting prognoses. These tumors are defined by recurrent mutations in the IDH genes and other genetic alterations including inactivation of CDKN2A and COL2A1; however, these have no clinical value. Here we use multi-omics molecular profiles from a series of cartilage tumors and find an mRNA classification that identifies two subtypes of chondrosarcomas defined by a balance in tumor differentiation and cell cycle activation. The microRNA classification reveals the importance of the loss of expression of the 14q32 locus in defining the level of malignancy. Finally, DNA methylation is associated with IDH mutations. We can use the multi-omics classifications to predict outcome. We propose an mRNA-only classifier to reproduce the integrated multi-omics classification, and its application to relapsed tumor samples shows the progressive nature of the classification. Thus, it may be possible to use mRNA-based signatures to detect patients with high-risk chondrosarcomas.