A mucus production programme promotes classical pancreatic ductal adenocarcinoma.

OBJECTIVE: The optimal therapeutic response in cancer patients is highly dependent upon the differentiation state of their tumours. Pancreatic ductal adenocarcinoma (PDA) is a lethal cancer that harbours distinct phenotypic subtypes with preferential sensitivities to standard therapies. This study aimed to investigate intratumour heterogeneity and plasticity of cancer cell states in PDA in order to reveal cell state-specific regulators. DESIGN: We analysed single-cell expression profiling of mouse PDAs, revealing intratumour heterogeneity and cell plasticity and identified pathways activated in the different cell states. We performed comparative analysis of murine and human expression states and confirmed their phenotypic diversity in specimens by immunolabeling. We assessed the function of phenotypic regulators using mouse models of PDA, organoids, cell lines and orthotopically grafted tumour models. RESULTS: Our expression analysis and immunolabeling analysis show that a mucus production programme regulated by the transcription factor SPDEF is highly active in precancerous lesions and the classical subtype of PDA - the most common differentiation state. SPDEF maintains the classical differentiation and supports PDA transformation in vivo. The SPDEF tumour-promoting function is mediated by its target genes AGR2 and ERN2/IRE1ß that regulate mucus production, and inactivation of the SPDEF programme impairs tumour growth and facilitates subtype interconversion from classical towards basal-like differentiation. CONCLUSIONS: Our findings expand our understanding of the transcriptional programmes active in precancerous lesions and PDAs of classical differentiation, determine the regulators of mucus production as specific vulnerabilities in these cell states and reveal phenotype switching as a response mechanism to inactivation of differentiation states determinants.

Genetic Ancestry Inference from Cancer-Derived Molecular Data across Genomic and Transcriptomic Platforms.

Genetic ancestry-oriented cancer research requires the ability to perform accurate and robust genetic ancestry inference from existing cancer-derived data, including whole-exome sequencing, transcriptome sequencing, and targeted gene panels, very often in the absence of matching cancer-free genomic data. Here we examined the feasibility and accuracy of computational inference of genetic ancestry relying exclusively on cancer-derived data. A data synthesis framework was developed to optimize and assess the performance of the ancestry inference for any given input cancer-derived molecular profile. In its core procedure, the ancestral background of the profiled patient is replaced with one of any number of individuals with known ancestry. The data synthesis framework is applicable to multiple profiling platforms, making it possible to assess the performance of inference specifically for a given molecular profile and separately for each continental-level ancestry; this ability extends to all ancestries, including those without statistically sufficient representation in the existing cancer data. The inference procedure was demonstrated to be accurate and robust in a wide range of sequencing depths. Testing of the approach in four representative cancer types and across three molecular profiling modalities showed that continental-level ancestry of patients can be inferred with high accuracy, as quantified by its agreement with the gold standard of deriving ancestry from matching cancer-free molecular data. This study demonstrates that vast amounts of existing cancer-derived molecular data are potentially amenable to ancestry-oriented studies of the disease without requiring matching cancer-free genomes or patient self-reported ancestry. SIGNIFICANCE: The development of a computational approach that enables accurate and robust ancestry inference from cancer-derived molecular profiles without matching cancer-free data provides a valuable methodology for genetic ancestry-oriented cancer research.

Early-Life Exposure to Environmental Contaminants Perturbs the Sperm Epigenome and Induces Negative Pregnancy Outcomes for Three Generations via the Paternal Lineage.

Due to the grasshopper effect, the Arctic food chain in Canada is contaminated with persistent organic pollutants (POPs) of industrial origin, including polychlorinated biphenyls and organochlorine pesticides. Exposure to POPs may be a contributor to the greater incidence of poor fetal growth, placental abnormalities, stillbirths, congenital defects and shortened lifespan in the Inuit population compared to non-Aboriginal Canadians. Although maternal exposure to POPs is well established to harm pregnancy outcomes, paternal transmission of the effects of POPs is a possibility that has not been well investigated. We used a rat model to test the hypothesis that exposure to POPs during gestation and suckling leads to developmental defects that are transmitted to subsequent generations via the male lineage. Indeed, developmental exposure to an environmentally relevant Arctic POPs mixture impaired sperm quality and pregnancy outcomes across two subsequent, unexposed generations and altered sperm DNA methylation, some of which are also observed for two additional generations. Genes corresponding to the altered sperm methylome correspond to health problems encountered in the Inuit population. These findings demonstrate that the paternal methylome is sensitive to the environment and that some perturbations persist for at least two subsequent generations. In conclusion, although many factors influence health, paternal exposure to contaminants plays a heretofore-underappreciated role with sperm DNA methylation contributing to the molecular underpinnings involved.

Intraductal Transplantation Models of Human Pancreatic Ductal Adenocarcinoma Reveal Progressive Transition of Molecular Subtypes.

Pancreatic ductal adenocarcinoma (PDAC) is the most lethal common malignancy, with little improvement in patient outcomes over the past decades. Recently, subtypes of pancreatic cancer with different prognoses have been elaborated; however, the inability to model these subtypes has precluded mechanistic investigation of their origins. Here, we present a xenotransplantation model of PDAC in which neoplasms originate from patient-derived organoids injected directly into murine pancreatic ducts. Our model enables distinction of the two main PDAC subtypes: intraepithelial neoplasms from this model progress in an indolent or invasive manner representing the classical or basal-like subtypes of PDAC, respectively. Parameters that influence PDAC subtype specification in this intraductal model include cell plasticity and hyperactivation of the RAS pathway. Finally, through intratumoral dissection and the direct manipulation of RAS gene dosage, we identify a suite of RAS-regulated secreted and membrane-bound proteins that may represent potential candidates for therapeutic intervention in patients with PDAC. SIGNIFICANCE: Accurate modeling of the molecular subtypes of pancreatic cancer is crucial to facilitate the generation of effective therapies. We report the development of an intraductal organoid transplantation model of pancreatic cancer that models the progressive switching of subtypes, and identify stochastic and RAS-driven mechanisms that determine subtype specification.See related commentary by Pickering and Morton, p. 1448.This article is highlighted in the In This Issue feature, p. 1426.

SOAT1 promotes mevalonate pathway dependency in pancreatic cancer.

Pancreatic ductal adenocarcinoma (PDAC) has a dismal prognosis, and new therapies are needed. Altered metabolism is a cancer vulnerability, and several metabolic pathways have been shown to promote PDAC. However, the changes in cholesterol metabolism and their role during PDAC progression remain largely unknown. Here we used organoid and mouse models to determine the drivers of altered cholesterol metabolism in PDAC and the consequences of its disruption on tumor progression. We identified sterol O-acyltransferase 1 (SOAT1) as a key player in sustaining the mevalonate pathway by converting cholesterol to inert cholesterol esters, thereby preventing the negative feedback elicited by unesterified cholesterol. Genetic targeting of Soat1 impairs cell proliferation in vitro and tumor progression in vivo and reveals a mevalonate pathway dependency in p53 mutant PDAC cells that have undergone p53 loss of heterozygosity (LOH). In contrast, pancreatic organoids lacking p53 mutation and p53 LOH are insensitive to SOAT1 loss, indicating a potential therapeutic window for inhibiting SOAT1 in PDAC.

Bioactivation of Napabucasin Triggers Reactive Oxygen Species-Mediated Cancer Cell Death.

PURPOSE: Napabucasin (2-acetylfuro-1,4-naphthoquinone or BBI-608) is a small molecule currently being clinically evaluated in various cancer types. It has mostly been recognized for its ability to inhibit STAT3 signaling. However, based on its chemical structure, we hypothesized that napabucasin is a substrate for intracellular oxidoreductases and therefore may exert its anticancer effect through redox cycling, resulting in reactive oxygen species (ROS) production and cell death. EXPERIMENTAL DESIGN: Binding of napabucasin to NAD(P)H:quinone oxidoreductase-1 (NQO1), and other oxidoreductases, was measured. Pancreatic cancer cell lines were treated with napabucasin, and cell survival, ROS generation, DNA damage, transcriptomic changes, and alterations in STAT3 activation were assayed in vitro and in vivo. Genetic knockout or pharmacologic inhibition with dicoumarol was used to evaluate the dependency on NQO1. RESULTS: Napabucasin was found to bind with high affinity to NQO1 and to a lesser degree to cytochrome P450 oxidoreductase (POR). Treatment resulted in marked induction of ROS and DNA damage with an NQO1- and ROS-dependent decrease in STAT3 phosphorylation. Differential cytotoxic effects were observed, where NQO1-expressing cells generating cytotoxic levels of ROS at low napabucasin concentrations were more sensitive. Cells with low or no baseline NQO1 expression also produced ROS in response to napabucasin, albeit to a lesser extent, through the one-electron reductase POR. CONCLUSIONS: Napabucasin is bioactivated by NQO1, and to a lesser degree by POR, resulting in futile redox cycling and ROS generation. The increased ROS levels result in DNA damage and multiple intracellular changes, one of which is a reduction in STAT3 phosphorylation.

Organoid Profiling Identifies Common Responders to Chemotherapy in Pancreatic Cancer.

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.Significance: 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. Cancer Discov; 8(9); 1112-29. ©2018 AACR.See related commentary by Collisson, p. 1062This article is highlighted in the In This Issue feature, p. 1047.

Inferring and modeling inheritance of differentially methylated changes across multiple generations.

High-throughput methylation sequencing enables genome-wide detection of differentially methylated sites (DMS) or regions (DMR). Increasing evidence suggests that treatment-induced DMS can be transmitted across generations, but the analysis of induced methylation changes across multiple generations is complicated by the lack of sound statistical methods to evaluate significance levels. Due to software design, DMS detection was usually made on each generation separately, thus disregarding stochastic effects expected when a large number of DMS is detected in each generation. Here, we present a novel method based on Monte Carlo sampling, methylInheritance, to evaluate that the number of conserved DMS between several generations is associated to an effect inherited from a treatment and not randomness. Moreover, we developed an inheritance simulation package, methInheritSim, to demonstrate the performance of the methylInheritance method and to evaluate the power of different experimental designs. Finally, we applied methylInheritance to a DNA methylation dataset obtained from early-life persistent organic pollutants (POPs) exposed Sprague-Dawley female rats and their descendants through a paternal transmission. The results show that metylInheritance can efficiently identify treatment-induced inherited methylation changes. Specifically, we identified two intergenerationally conserved DMS at transcription start site (TSS); one of those persisted transgenerationally. Three transgenerationally conserved DMR were found at intra or integenic regions.

metagene Profiles Analyses Reveal Regulatory Element's Factor-Specific Recruitment Patterns.

ChIP-Sequencing (ChIP-Seq) provides a vast amount of information regarding the localization of proteins across the genome. The aggregation of ChIP-Seq enrichment signal in a metagene plot is an approach commonly used to summarize data complexity and to obtain a high level visual representation of the general occupancy pattern of a protein. Here we present the R package metagene, the graphical interface Imetagene and the companion package similaRpeak. Together, they provide a framework to integrate, summarize and compare the ChIP-Seq enrichment signal from complex experimental designs. Those packages identify and quantify similarities or dissimilarities in patterns between large numbers of ChIP-Seq profiles. We used metagene to investigate the differential occupancy of regulatory factors at noncoding regulatory regions (promoters and enhancers) in relation to transcriptional activity in GM12878 B-lymphocytes. The relationships between occupancy patterns and transcriptional activity suggest two different mechanisms of action for transcriptional control: i) a "gradient effect" where the regulatory factor occupancy levels follow transcription and ii) a "threshold effect" where the regulatory factor occupancy levels max out prior to reaching maximal transcription. metagene, Imetagene and similaRpeak are implemented in R under the Artistic license 2.0 and are available on Bioconductor.

Development of highly reliable in silico SNP resource and genotyping assay from exome capture and sequencing: an example from black spruce (Picea mariana).

Picea mariana is a widely distributed boreal conifer across Canada and the subject of advanced breeding programmes for which population genomics and genomic selection approaches are being developed. Targeted sequencing was achieved after capturing P. mariana exome with probes designed from the sequenced transcriptome of Picea glauca, a distant relative. A high capture efficiency of 75.9% was reached although spruce has a complex and large genome including gene sequences interspersed by some long introns. The results confirmed the relevance of using probes from congeneric species to perform successfully interspecific exome capture in the genus Picea. A bioinformatics pipeline was developed including stringent criteria that helped detect a set of 97,075 highly reliable in silico SNPs. These SNPs were distributed across 14,909 genes. Part of an Infinium iSelect array was used to estimate the rate of true positives by validating 4267 of the predicted in silico SNPs by genotyping trees from P. mariana populations. The true positive rate was 96.2% for in silico SNPs, compared to a genotyping success rate of 96.7% for a set 1115 P. mariana control SNPs recycled from previous genotyping arrays. These results indicate the high success rate of the genotyping array and the relevance of the selection criteria used to delineate the new P. mariana in silico SNP resource. Furthermore, in silico SNPs were generally of medium to high frequency in natural populations, thus providing high informative value for future population genomics applications.

Using informative Multinomial-Dirichlet prior in a t-mixture with reversible jump estimation of nucleosome positions for genome-wide profiling.

Genome-wide mapping of nucleosomes has revealed a great deal about the relationships between chromatin structure and control of gene expression. Recent next generation CHIP-chip and CHIP-Seq technologies have accelerated our understanding of basic principles of chromatin organization. These technologies have taught us that nucleosomes play a crucial role in gene regulation by allowing physical access to transcription factors. Recent methods and experimental advancements allow the determination of nucleosome positions for a given genome area. However, most of these methods estimate the number of nucleosomes either by an EM algorithm using a BIC criterion or an effective heuristic strategy. Here, we introduce a Bayesian method for identifying nucleosome positions. The proposed model is based on a Multinomial-Dirichlet classification and a hierarchical mixture distributions. The number and the positions of nucleosomes are estimated using a reversible jump Markov chain Monte Carlo simulation technique. We compare the performance of our method on simulated data and MNase-Seq data from Saccharomyces cerevisiae against PING and NOrMAL methods.

The landscape of nucleotide polymorphism among 13,500 genes of the conifer picea glauca, relationships with functions, and comparison with medicago truncatula.

Gene families differ in composition, expression, and chromosomal organization between conifers and angiosperms, but little is known regarding nucleotide polymorphism. Using various sequencing strategies, an atlas of 212k high-confidence single nucleotide polymorphisms (SNPs) with a validation rate of more than 92% was developed for the conifer white spruce (Picea glauca). Nonsynonymous and synonymous SNPs were annotated over the corresponding 13,498 white spruce genes representative of 2,457 known gene families. Patterns of nucleotide polymorphisms were analyzed by estimating the ratio of nonsynonymous to synonymous numbers of substitutions per site (A/S). A general excess of synonymous SNPs was expected and observed. However, the analysis from several perspectives enabled to identify groups of genes harboring an excess of nonsynonymous SNPs, thus potentially under positive selection. Four known gene families harbored such an excess: dehydrins, ankyrin-repeats, AP2/DREB, and leucine-rich repeat. Conifer-specific sequences were also generally associated with the highest A/S ratios. A/S values were also distributed asymmetrically across genes specifically expressed in megagametophytes, roots, or in both, harboring on average an excess of nonsynonymous SNPs. These patterns confirm that the breadth of gene expression is a contributing factor to the evolution of nucleotide polymorphism. The A/S ratios of Medicago truncatula genes were also analyzed: several gene families shared between P. glauca and M. truncatula data sets had similar excess of synonymous or nonsynonymous SNPs. However, a number of families with high A/S ratios were found specific to P. glauca, suggesting cases of divergent evolution at the functional level.

Development of high-density SNP genotyping arrays for white spruce (Picea glauca) and transferability to subtropical and nordic congeners.

High-density SNP genotyping arrays can be designed for any species given sufficient sequence information of high quality. Two high-density SNP arrays relying on the Infinium iSelect technology (Illumina) were designed for use in the conifer white spruce (Picea glauca). One array contained 7338 segregating SNPs representative of 2814 genes of various molecular functional classes for main uses in genetic association and population genetics studies. The other one contained 9559 segregating SNPs representative of 9543 genes for main uses in population genetics, linkage mapping of the genome and genomic prediction. The SNPs assayed were discovered from various sources of gene resequencing data. SNPs predicted from high-quality sequences derived from genomic DNA reached a genotyping success rate of 64.7%. Nonsingleton in silico SNPs (i.e. a sequence polymorphism present in at least two reads) predicted from expressed sequenced tags obtained with the Roche 454 technology and Illumina GAII analyser resulted in a similar genotyping success rate of 71.6% when the deepest alignment was used and the most favourable SNP probe per gene was selected. A variable proportion of these SNPs was shared by other nordic and subtropical spruce species from North America and Europe. The number of shared SNPs was inversely proportional to phylogenetic divergence and standing genetic variation in the recipient species, but positively related to allele frequency in P. glauca natural populations. These validated SNP resources should open up new avenues for population genetics and comparative genetic mapping at a genomic scale in spruce species.