Cytidine deaminases APOBEC3C and APOBEC3D promote DNA replication stress resistance in pancreatic cancer cells

Gemcitabine is a potent inhibitor of DNA replication and is a mainstay therapeutic for diverse cancers, particularly pancreatic ductal adenocarcinoma (PDAC). However, most tumors remain refractory to gemcitabine therapies. Here, to define the cancer cell response to gemcitabine, we performed genome-scale CRISPR-Cas9 chemical-genetic screens in PDAC cells and found selective loss of cell fitness upon disruption of the cytidine deaminases APOBEC3C and APOBEC3D. Following gemcitabine treatment, APOBEC3C and APOBEC3D promote DNA replication stress resistance and cell survival by deaminating cytidines in the nuclear genome to ensure DNA replication fork restart and repair in PDAC cells. We provide evidence that the chemical-genetic interaction between APOBEC3C or APOBEC3D and gemcitabine is absent in nontransformed cells but is recapitulated across different PDAC cell lines, in PDAC organoids and in PDAC xenografts. Thus, we uncover roles for APOBEC3C and APOBEC3D in DNA replication stress resistance and offer plausible targets for improving gemcitabine-based therapies for PDAC.

Differences in the transcriptome response in the gills of sea lamprey acutely exposed to 3-trifluoromethyl-4-nitrophenol (TFM), niclosamide or a TFM:niclosamide mixture

Sea lamprey (Petromyzon marinus) control in the Laurentian Great Lakes of North America makes use of two pesticides: 3-trifluoromethyl-4-nitrophenol (TFM) and niclosamide, which are often co-applied. Sea lamprey appear to be vulnerable to these agents resulting from a lack of detoxification responses with evidence suggesting that lampricide mixtures produce a synergistic effect. However, there is a lack of information pertaining to the physiological responses of sea lamprey to niclosamide and TFM:niclosamide mixtures. Here, we characterized the transcriptomic responses of the sea lamprey to TFM, niclosamide, and a TFM:niclosamide (1.5 %) mixture in the gill. Along with a control, larval sea lamprey were exposed to each treatment for 6 h, after which gill tissues were extracted for measuring whole-transcriptome responses using RNA sequencing. Differential gene expression patterns were summarized, which included identifying the broad roles of genes and common expression patterns among the treatments. While niclosamide treatment resulted in no differentially expressed genes, TFM- and mixture-treated fish had several differentially expressed genes that were associated with the cell cycle, DNA damage, metabolism, immune function, and detoxification. However, there was no common differential expression among treatments. For the first time, we characterized the transcriptomic response of sea lamprey to niclosamide and a TFM:niclosamide mixture and identified that these agents impact mRNA transcript abundance of genes associated with the cell cycle and cellular death, and immune function, which are likely mediated through mitochondrial dysregulation. These results may help to inform the production of more targeted and effective lampricides in sea lamprey control efforts.

Chemokine expression predicts T cell-inflammation and improved survival with checkpoint inhibition across solid cancers

Immune checkpoint inhibitors (ICI) are highly effective in specific cancers where canonical markers of antitumor immunity are used for patient selection. Improved predictors of T cell-inflammation are needed to identify ICI-responsive tumor subsets in additional cancer types. We investigated associations of a 4-chemokine expression signature (c-Score: CCL4, CCL5, CXCL9, CXCL10) with metrics of antitumor immunity across tumor types. Across cancer entities from The Cancer Genome Atlas, subgroups of tumors displayed high expression of the c-Score (c-Scorehi) with increased expression of immune checkpoint (IC) genes and transcriptional hallmarks of the cancer-immunity cycle. There was an incomplete association of the c-Score with high tumor mutation burden (TMB), with only 15% of c-Scorehi tumors displaying ≥10 mutations per megabase. In a heterogeneous pan-cancer cohort of 82 patients, with advanced and previously treated solid cancers, c-Scorehi tumors had a longer median time to progression (103 versus 72 days, P = 0.012) and overall survival (382 versus 196 days, P = 0.038) following ICI therapy initiation, compared to patients with low c-Score expression. We also found c-Score stratification to outperform TMB assignment for overall survival prediction (HR = 0.42 [0.22-0.79], P = 0.008 versus HR = 0.60 [0.29-1.27], P = 0.18, respectively). Assessment of the c-Score using the TIDE and PredictIO databases, which include ICI treatment outcomes from 10 tumor types, provided further support for the c-Score as a predictive ICI therapeutic biomarker. In summary, the c-Score identifies patients with hallmarks of T cell-inflammation and potential response to ICI treatment across cancer types, which is missed by TMB assignment.

Spectrum of Response to Platinum and PARP Inhibitors in Germline BRCA-Associated Pancreatic Cancer in the Clinical and Preclinical Setting

Germline BRCA-associated pancreatic ductal adenocarcinoma (glBRCA PDAC) tumors are susceptible to platinum and PARP inhibition. The clinical outcomes of 125 patients with glBRCA PDAC were stratified based on the spectrum of response to platinum/PARP inhibition: (i) refractory [overall survival (OS) <6 months], (ii) durable response followed by acquired resistance (OS <36 months), and (iii) long-term responders (OS >36 months). Patient-derived xenografts (PDX) were generated from 25 patients with glBRCA PDAC at different clinical time points. Response to platinum/PARP inhibition in vivo and ex vivo culture (EVOC) correlated with clinical response. We deciphered the mechanisms of resistance in glBRCA PDAC and identified homologous recombination (HR) proficiency and secondary mutations restoring partial functionality as the most dominant resistant mechanism. Yet, a subset of HR-deficient (HRD) patients demonstrated clinical resistance. Their tumors displayed basal-like molecular subtype and were more aneuploid. Tumor mutational burden was high in HRD PDAC and significantly higher in tumors with secondary mutations. Anti-PD-1 attenuated tumor growth in a novel humanized glBRCA PDAC PDX model. This work demonstrates the utility of preclinical models, including EVOC, to predict the response of glBRCA PDAC to treatment, which has the potential to inform time-sensitive medical decisions.

SIGNIFICANCE

glBRCA PDAC has a favorable response to platinum/PARP inhibition. However, most patients develop resistance. Additional treatment options for this unique subpopulation are needed. We generated model systems in PDXs and an ex vivo system (EVOC) that faithfully recapitulate these specific clinical scenarios as a platform to investigate the mechanisms of resistance for further drug development.

Systemic inflammatory prognostic scores in advanced pancreatic adenocarcinoma

BACKGROUND

Systemic inflammatory scores may aid prognostication and patient selection for trials. We compared five scores in advanced pancreatic adenocarcinoma (PDAC).

METHODS

Unresectable/metastatic PDAC patients enrolled in the Comprehensive Molecular Characterisation of Advanced Pancreatic Ductal Adenocarcinoma for Better Treatment Selection trial (NCT02750657) were included. Patients had pre-treatment biopsies for whole genome and RNA sequencing. CD8 immunohistochemistry was available in a subset. The neutrophil-to-lymphocyte ratio (NLR), platelet-to-lymphocyte ratio, Prognostic Nutritional Index, Gustave Roussy Immune Score (GRIm-S), and Memorial Sloan Kettering Prognostic Score (MPS) were calculated. Overall survival (OS) was estimated using Kaplan-Meier methods. Associations between inflammatory scores, clinical/genomic characteristics, and OS were analysed.

RESULTS

We analysed 263 patients. High-risk NLR, GRIm-S and MPS were poorly prognostic. The GRIm-S had the highest predictive ability: median OS 6.4 vs. 10 months for high risk vs. low-risk (P < 0.001); HR 2.26 (P < 0.001). ECOG ≥ 1, the basal-like subtype, and low-HRDetect were additional poor prognostic factors (P < 0.01). Inflammatory scores did not associate with RNA-based classifiers or homologous recombination repair deficiency genotypes. High-risk MPS (P = 0.04) and GRIm-S (P = 0.02) patients had lower median CD8 + tumour-infiltrating lymphocytes.

CONCLUSIONS

Inflammatory scores incorporating NLR have prognostic value in advanced PDAC. Understanding immunophenotypes of poor-risk patients and using these scores in trials will advance the field.

Molecular characterisation of pancreatic ductal adenocarcinoma with NTRK fusions and review of the literature

AIMS

The majority of pancreatic ductal adenocarcinomas (PDACs) harbour oncogenic mutations in KRAS with variants in TP53, CDKN2A and SMAD4 also prevalent. The presence of oncogenic fusions including NTRK fusions are rare but important to identify. Here we ascertain the prevalence of NTRK fusions and document their genomic characteristics in a large series of PDAC.

METHODS

Whole genome sequencing and RNAseq were performed on a series of patients with resected or locally advanced/metastatic PDAC collected between 2008 and 2020 at a single institution. A subset of specimens underwent immunohistochemistry (IHC) analysis. Clinical and molecular characterisation and IHC sensitivity and specificity were evaluated.

RESULTS

400 patients were included (resected n=167; locally advanced/metastatic n=233). Three patients were identified as harbouring an NTRK fusion, two EML4-NTRK3 (KRAS-WT) and a single novel KANK1-NTRK3 fusion. The latter occurring in the presence of a subclonal KRAS mutation. Typical PDAC drivers were present including mutations in TP53 and CDKN2A. Substitution base signatures and tumour mutational burden were similar to typical PDAC. The prevalence of NTRK fusions was 0.8% (3/400), while in KRAS wild-type tumours, it was 6.25% (2/32). DNA prediction alone documented six false-positive cases. RNA analysis correctly identified the in-frame fusion transcripts. IHC analysis was negative in the KANK1-NTRK3 fusion but positive in a EML4-NTRK3 case, highlighting lower sensitivity of IHC.

CONCLUSION

NTRK fusions are rare; however, with emerging therapeutic options targeting these fusions, detection is vital. Reflex testing for KRAS mutations and subsequent RNA-based screening could help identify these cases in PDAC.

Abstract A036: Plasma extracellular RNA for the detection of pancreatic cancer

Pancreatic ductal adenocarcinoma (PDA) is expected to become the second leading cause of cancer-related death by 2030. Patients are often diagnosed at the metastatic stage, at which point treatments are limited and have non-uniform efficacy in patients. Advances in transcriptomic subtyping of tumor tissue can segregate patients into two broad subtypes, Classical and Basal-like, which are predictive of survival and response to first-line chemotherapy. However, tumor tissue is often difficult or impossible to obtain, particularly in patients with advanced disease, hindering the clinical implementation of transcriptomic subtyping. Profiling the extracellular RNA (exRNA) content in plasma provides an opportunity for a non-invasive readout of the tumor transcriptome. Here, we aim to establish the feasibility of liquid biopsy—particularly, plasma exRNA-profiling—for the diagnosis and subtyping of PDA. Extracellular vesicle (EV) encapsulated exRNA was isolated from 1mL of frozen plasma collected between 2017-2020 from a pilot cohort of 8 patients with stage III-IV PDA and processed alongside plasma from 6 healthy volunteers. UMI-tagged libraries were constructed from low-input exRNA to generate sequencing data. Gene body coverage analysis indicated that exRNAs consisted of intact full-length transcripts. Of these transcripts, we found that there was a median of 34,975 genes encoded by exRNAs, deriving mostly from protein-coding regions. Most genes captured by exRNA sequencing had low expression. Comparing between groups, we detected a significantly greater number of genes and novel splicing junctions in PDA patients. Of the cancer-specific transcripts, differential expression analysis revealed significant upregulation of membrane organelle and RNA-processing related pathways. PDA exRNA profiles did not completely recapitulate the bulk tumor transcriptome, suggesting specific mechanisms for EV encapsulation of RNAs and exRNA release. To classify tumor subtypes by plasma exRNA sequencing, we identified an exRNA-specific gene panel from bulk tumor transcriptome signatures. Using this gene panel, tumor subtype was correctly inferred from 6/7 PDA patients. A combined gene panel derived from all tumor-associated exRNA signatures distinguished 7/8 PDA patients from healthy controls. In conclusion, plasma exRNA recapitulated elements of the tumor transcriptome and has utility in PDA diagnosis. The establishment of diagnostic markers will be further investigated in the expanded study cohort comprised of 91 unique stage III-IV PDA donors with plasma collected at baseline, 150 stage III-IV PDA donors with plasma banked at timepoints throughout chemotherapy treatment, and 100 age- and sex-matched non-PDA controls donors. Leveraging the study cohort, future work will also investigate plasma exRNA signatures as prognostic markers to predict treatment outcomes and monitor treatment responses in patients with advanced PDA.

Citation Format: Linsey Gong, Ling Xi, Yuanchang Fang, Karen Ng, Gun Ho Jang, Amy X. Zhang, Julie M. Wilson, Stephanie Ramotar, Jennifer J. Knox, Grainne M. O'Kane, Steven Gallinger, Hansen H. He, Faiyaz Notta. Plasma extracellular RNA for the detection of pancreatic cancer [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer; 2022 Sep 13-16; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2022;82(22 Suppl):Abstract nr A036.

Abstract B031: Loss of GATA6 leads to reduced cell fitness and PDA molecular subtype changes

Background: Pancreatic ductal adenocarcinoma (PDA) is defined by two main transcriptional subtypes, Basal-like & Classical. Classical tumors are more frequent and tend to respond better to chemotherapy with higher expression of genes related to pancreatic lineage development, the epithelial compartment, and adhesion. Conversely, Basal-like tumors are more aggressive, therapy resistant and have higher expression of laminins and keratins. The Classical subtype correlates with GATA6 expression, and ~16% pancreatic tumors show amplification in the transcription factor GATA6 suggesting its importance in the Classical subtype and in tumor survival (Chan-Seng-Yue et al., 2020). The goal of this project is to elucidate the role of GATA6 expression in PDA tumors showing that GATA6 is a key molecular regulator of the Classical subtype influencing tumor growth, maintenance, and phenotype. Methods: Patient derived organoids (PDOs) spanning the spectrum of GATA6 expression were used with CRISPR/Cas9 to establish GATA6 KO lines (n=8). These were then subject to a growth competition assay between cells with and without GATA6 KO and the amount of knockout over time was monitored by Sanger sequencing. Tumor growth curves were generated in a subcutaneous xenograft mouse model over 1.5 months to assess the in vivo growth potential of GATA6 KO cells in a sample with high GATA6 expression. 3’RNA-seq was used to assess gene expression changes in subtype, and GSEA was used to assess pathway enrichment. This was done alongside IHC, capillary Westerns, and flow cytometry to assess changes that GATA6 KO incurs on a functional level. Results: GATA6 KO tended to result in reduced fitness (dropout) in cells when ablated, although the rate and level of dropout varied, suggesting different levels of dependence on GATA6 expression for survival/growth. This was confirmed in vivo when the GATA6 KO failed to grow as robustly as control groups in a sample with high native GATA6 expression. Molecular analysis of all PDO samples showed shifts in expression from Classical to Basal-like. This shift was denoted in changes in subtype score, changes in marker gene expression (i.e., ANXA10 & TP63), and in changes in key pathways associated with each subtype. This shows that GATA6 not only correlates with the Classical subtype but plays a causal role in maintaining it. Even though there was a shift towards Basal-like expression several key characteristics of the subtype were notably absent including enrichment for epithelial to mesenchymal transition, down regulation of MYC and its targets, and the lack of aggressive growth in in vitro and in vivo models (dropout of GATA6 KO cells, low Ki67 expression, and upregulation of apoptosis). This suggests that KO of GATA6 is not sufficient in of itself to cause a true Basal-like phenotype and that instead GATA6 KO might lead to an in-between state that does not quite offer the same growth advantages and may open GATA6 and its pathways to be a potential therapeutic target.

Citation Format: Tristan S.D. Woo, Kena F. Figueroa, Sabrina Ge, Irene Xie, Karen Ng, Pinjiang Cao, Julie M. Wilson, Stephanie Ramotar, Jennifer J. Knox, Grainne M. O'Kane, Steven Gallinger, Faiyaz Notta. Loss of GATA6 leads to reduced cell fitness and PDA molecular subtype changes [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer; 2022 Sep 13-16; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2022;82(22 Suppl):Abstract nr B031.

Abstract A003: Developmental parallels in pancreatic cancer subtypes revealed by single cell transcriptomics

Recent research into the transcriptional subtypes in pancreatic ductal adenocarcinoma (PDA) has demonstrated their importance in predicting survival and therapeutic response. The Basal-like subtype is enriched in metastases and generally predicts poorer prognosis, while the Classical subtype is defined by gene expression reflective of pancreatic progenitors. The two-subtype scheme can be further subdivided into five subtypes that help better reflect the vast transcriptional heterogeneity found in PDA tumors. Basal-like A tumors fare worse than Basal-like B tumors. KRAS wildtype tumors are enriched in Classical B, which have the best prognosis, while Classical A is more reflective of the two-subtype Classical program. Hybrid tumors, generally a mix of Basal-like B and Classical A, are intermediate in prognosis between the Classical and Basal-like subtypes. Despite the importance of the subtypes, we have yet to fully understand how subtype programs within single cells of a tumor are intrinsically regulated, and how they contribute to tumor-level transcriptional observations. To investigate, we generated and leveraged a single cell RNA sequencing dataset enriched for malignant tumor cells spanning both primary and metastatic disease. Along with the richly annotated COMPASS dataset featuring bulk RNA sequencing of over 300 laser captured samples, we developed a single sample classifier utilizing gene pairs for assigning Basal-like A, Basal-like B, Hybrid, Classical A or Classical B to bulk tumors as well as individual cells, allowing interrogation of the subtypes across modalities. We find that the bulk subtype of the tumor is generally reflected in the dominant cell subtype. However, there are frequently cells of other subtypes within individual tumors, and substantial differences in cell subtype composition across tumors of the same subtype. In addition, we find that the proportion of Basal-like cells is a significant risk prognosticator. Investigating deeper into the transcriptional drivers of the subtypes reveals differences in expression of programs linked to endodermal development. Transcription factors such as FOXC1 and IRX3 generally related to foregut development are enriched in Basal-like A, forming a gradient of high expression to lower expression from Basal-like A, Basal-like B, Hybrid, Classical A to Classical B. Interestingly, midgut transcription factors such as CDX2 and ELF3 are expressed in the reverse order. Integrating data from single cell endoderm RNA sequencing further supports this finding, suggesting that the subtypes may arise from tumor cells exploiting pathways normally expressed early in development. This work provides a stepping stone to better understanding the mechanisms behind transcriptional subtypes in PDA.

Citation Format: Sabrina Ge, Gun Ho Jang, Michelle Chan-Seng-Yue, Eugenia Flores Figueroa, Karen Ng, Stephanie Ramotar, Anna Dodd, Grainne M. O'Kane, Julie M. Wilson, Jennifer J. Knox, Sandra E. Fischer, Steven Gallinger, Faiyaz Notta. Developmental parallels in pancreatic cancer subtypes revealed by single cell transcriptomics [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer; 2022 Sep 13-16; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2022;82(22 Suppl):Abstract nr A003.

Integrative analysis of KRAS wildtype metastatic pancreatic ductal adenocarcinoma reveals mutation and expression-based similarities to cholangiocarcinoma

Oncogenic KRAS mutations are absent in approximately 10% of patients with metastatic pancreatic ductal adenocarcinoma (mPDAC) and may represent a subgroup of mPDAC with therapeutic options beyond standard-of-care cytotoxic chemotherapy. While distinct gene fusions have been implicated in KRAS wildtype mPDAC, information regarding other types of mutations remain limited, and gene expression patterns associated with KRAS wildtype mPDAC have not been reported. Here, we leverage sequencing data from the PanGen trial to perform comprehensive characterization of the molecular landscape of KRAS wildtype mPDAC and reveal increased frequency of chr1q amplification encompassing transcription factors PROX1 and NR5A2. By leveraging data from colorectal adenocarcinoma and cholangiocarcinoma samples, we highlight similarities between cholangiocarcinoma and KRAS wildtype mPDAC involving both mutation and expression-based signatures and validate these findings using an independent dataset. These data further establish KRAS wildtype mPDAC as a unique molecular entity, with therapeutic opportunities extending beyond gene fusion events.

KRAS wildtype metastatic pancreatic ductal adenocarcinoma (mPDAC) could represent a distinct molecular entity from other PDACs. Here, the authors analyse KRAS wildtype mPDAC tumours using genomics and transcriptomics and find molecular similarities with cholangiocarcinomas.

Variation in the Transcriptome Response and Detoxification Gene Diversity Drives Pesticide Tolerance in Fishes

Pesticides are critical for invasive species management but often have negative effects on nontarget native biota. Tolerance to pesticides should have an evolutionary basis, but this is poorly understood. Invasive sea lamprey (Petromyzon marinus) populations in North America have been controlled with a pesticide lethal to them at lower concentrations than native fishes. We addressed how interspecific variation in gene expression and detoxification gene diversity confer differential pesticide sensitivity in two fish species. We exposed sea lamprey and bluegill (Lepomis macrochirus), a tolerant native species, to 3-trifluoromethyl-4-nitrophenol (TFM), a pesticide commonly used in sea lamprey control. We then used whole-transcriptome sequencing of gill and liver to characterize the cellular response in both species. Comparatively, bluegill exhibited a larger number of detoxification genes expressed and a larger number of responsive transcripts overall, which likely contributes to greater tolerance to TFM. Understanding the genetic and physiological basis for pesticide tolerance is crucial for managing invasive species.

Characterization of 3 different types of aquaporins in Carcinus maenas and their potential role in osmoregulation

Intertidal crustaceans like Carcinus maenas shift between an osmoconforming and osmoregulating state when inhabiting full-strength seawater and dilute environments, respectively. While the bodily fluids and environment of marine osmoconformers are approximately isosmotic, osmoregulating crabs inhabiting dilute environments maintain their bodily fluid osmolality above that of their environment by actively absorbing and retaining osmolytes (e.g., Na⁺, Cl^-, urea) while eliminating excess water. Few studies have investigated the role of aquaporins (AQPs) in the osmoregulatory organs of crustaceans, especially within brachyuran species. In the current study, three different aquaporins were identified within a transcriptome of C. maenas, including a classical AQP (CmAQP1), an aquaglyceroporin (CmGLP1), and a big-brain protein (CmBIB1), all of which are expressed in the gills and the antennal glands. Functional expression of these aquaporins confirmed water transport capabilities for CmAQP1, CmGLP1, but not for CmBIB1, while CmGLP1 also transported urea. Higher relative CmAQP1 mRNA expression within tissues of osmoconforming crabs suggests the apical/sub-apically localized channel attenuates osmotic gradients created by non-osmoregulatory processes while its downregulation in dilute media reduces the water permeability of tissues to facilitate osmoregulation. Although hemolymph urea concentrations rose upon exposure to brackish water, urea was not detected in the final urine. Due to its urea-transport capabilities, CmGLP1 is hypothesized to be involved in a urea retention mechanism believed to be involved in the production of diluted urine. Overall, these results suggest that AQPs are involved in osmoregulation and provide a basis for future mechanistic studies investigating the role of AQPs in volume regulation in crustaceans.

Distinct mechanisms of innate and adaptive immune regulation underlie poor oncologic outcomes associated with KRAS-TP53 co-alteration in pancreatic cancer

Co-occurrent KRAS and TP53 mutations define a majority of patients with pancreatic ductal adenocarcinoma (PDAC) and define its pro-metastatic proclivity. Here, we demonstrate that KRAS-TP53 co-alteration is associated with worse survival compared with either KRAS-alone or TP53-alone altered PDAC in 245 patients with metastatic disease treated at a tertiary referral cancer center, and validate this observation in two independent molecularly annotated datasets. Compared with non-TP53 mutated KRAS-altered tumors, KRAS-TP53 co-alteration engenders disproportionately innate immune-enriched and CD8+ T-cell-excluded immune signatures. Leveraging in silico, in vitro, and in vivo models of human and murine PDAC, we discover a novel intersection between KRAS-TP53 co-altered transcriptomes, TP63-defined squamous trans-differentiation, and myeloid-cell migration into the tumor microenvironment. Comparison of single-cell transcriptomes between KRAS-TP53 co-altered and KRAS-altered/TP53WT tumors revealed cancer cell-autonomous transcriptional programs that orchestrate innate immune trafficking and function. Moreover, we uncover granulocyte-derived inflammasome activation and TNF signaling as putative paracrine mediators of innate immunoregulatory transcriptional programs in KRAS-TP53 co-altered PDAC. Immune subtyping of KRAS-TP53 co-altered PDAC reveals conflation of intratumor heterogeneity with progenitor-like stemness properties. Coalescing these distinct molecular characteristics into a KRAS-TP53 co-altered "immunoregulatory program" predicts chemoresistance in metastatic PDAC patients enrolled in the COMPASS trial, as well as worse overall survival.

Is Biannual Surveillance for Pancreatic Cancer Sufficient in Individuals With Genetic Syndromes or Familial Pancreatic Cancer?

BACKGROUND

Individuals with a family history of pancreatic adenocarcinoma (PC) or with a germline mutation in a PC susceptibility gene are at increased risk of developing PC. These high-risk individuals (HRIs) may benefit from PC surveillance.

METHODS

A PC surveillance program was developed to evaluate the detection of premalignant lesions and early-stage PCs using biannual imaging and to determine whether locally advanced or metastatic PCs develop despite biannual surveillance. From January 2013 to April 2020, asymptomatic HRIs were enrolled and followed with alternating MRI and endoscopic ultrasound every 6 months.

RESULTS

Of 75 HRIs, 43 (57.3%) had a germline mutation in a PC susceptibility gene and 32 (42.7%) had a familial pancreatic cancer (FPC) pedigree. Branch-duct intraductal papillary mucinous neoplasms (BD-IPMNs) were identified in 26 individuals (34.7%), but only 2 developed progressive lesions. One patient with Peutz-Jeghers syndrome (PJS) developed locally advanced PC arising from a BD-IPMN. Whole-genome sequencing of this patient's PC and of a second patient with PJS-associated PC from the same kindred revealed biallelic inactivation of STK11 in a KRAS-independent manner. A review of 3,853 patients from 2 PC registries identified an additional patient with PJS-associated PC. All 3 patients with PJS developed advanced PC consistent with the malignant transformation of an underlying BD-IPMN in <6 months. The other surveillance patient with a progressive lesion had FPC and underwent resection of a mixed-type IPMN that harbored polyclonal KRAS mutations.

CONCLUSIONS

PC surveillance identifies a high prevalence of BD-IPMNs in HRIs. Patients with PJS with BD-IPMNs may be at risk for accelerated malignant transformation.

Intra-Tumoral CD8+ T-Cell Infiltration and PD-L1 Positivity in Homologous Recombination Deficient Pancreatic Ductal Adenocarcinoma

The immune contexture of pancreatic ductal adenocarcinoma (PDAC) is generally immunosuppressive. A role for immune checkpoint inhibitors (ICIs) in PDAC has only been demonstrated for the rare and hypermutated mismatch repair (MMR) deficient (MMR-d) subtype. Homologous recombination repair (HR) deficient (HR-d) PDAC is more prevalent and may encompass up to 20% of PDAC. Its genomic instability may promote a T-cell mediated anti-tumor response with therapeutic sensitivity to ICIs. To investigate the immunogenicity of HR-d PDAC, we used multiplex immunohistochemistry (IHC) to compare the density and spatial distribution of CD8+ cytotoxic T-cells, FOXP3+ regulatory T-cells (Tregs), and CD68+ tumor-associated macrophages (TAMs) in HR-d versus HR/MMR-intact PDAC. We also evaluated the IHC positivity of programmed death-ligand 1 (PD-L1) across the subgroups. 192 tumors were evaluated and classified as HR/MMR-intact (n=166), HR-d (n=25) or MMR-d (n=1) based on germline testing and tumor molecular hallmarks. Intra-tumoral CD8+ T-cell infiltration was higher in HR-d versus HR/MMR-intact PDAC (p<0.0001), while CD8+ T-cell densities in the peri-tumoral and stromal regions were similar in both groups. HR-d PDAC also displayed increased intra-tumoral FOXP3+ Tregs (p=0.049) and had a higher CD8+:FOXP3+ ratio (p=0.023). CD68+ TAM expression was similar in HR-d and HR/MMR-intact PDAC. Finally, 6 of the 25 HR-d cases showed a PD-L1 Combined Positive Score of >=1, whereas none of the HR/MMR-intact cases met this threshold (p<0.00001). These results provide immunohistochemical evidence for intra-tumoral CD8+ T-cell enrichment and PD-L1 positivity in HR-d PDAC, suggesting that HR-d PDAC may be amenable to ICI treatment strategies.

Abstract PO-077: Establishment of a novel living biobank of patient-derived pancreatic cancer organoids with genomic and drug response characterization

Advanced pancreatic cancer has a dismal prognosis and current treatment options (FOLFIRINOX, Gemcitabine/nab-paclitaxel [GnP]) are associated with toxicity. Although some patients achieve partial responses, most progress rapidly and become chemorefractory. While RNA subtypes, genomic alterations, and protein biomarkers have prognostic value, predictive biomarkers to guide therapy are needed. Patient-derived organoids (PDOs) are an increasingly popular model for predicting patient responses to standard-of-care therapy and investigating personalized therapy options. We present a novel biobank of 42 PDOs and drug profiling data with 5 standard of care agents and 3 kinase inhibitors. Tissue was processed from n=103 biopsies from 97 patients with a confirmed pathologic diagnosis of advanced (Stage III-IV) pancreatic ductal adenocarcinoma who presented to a single Canadian tertiary care centre between 2017-2020. Matched WGS was available in all cases. Our PDO generation success was 42/103 (41%). We observed a trend towards decreased establishment in tumors that were KRAS WT, TP53 WT, or had higher HRDetect scores. Conversely, polyploidy, SMAD4 WT, and major imbalances in mutant KRAS were associated with successful PDO establishment. These associations were not statistically significant after multiple comparisons correction, but suggest selection for success with more aggressive tumors. Drug profiling was performed on all 42 PDOs with the individual agents of FOLFIRINOX (5-FU, irinotecan, oxaliplatin), GnP (gemcitabine, paclitaxel), and three targeted agents (afatinib, trametinib, and talazoparib). Combination testing was also performed for gemcitabine + paclitaxel. Drug responses were measured through both viability and growth rate (GRMetrics). We found that GRMetrics minimized effects from different PDO growth rates. Matched clinical data were available for 23 patients who received FOLFIRINOX, 11 patients who received GnP, and one patient who received gemcitabine monotherapy. Similar to previous studies, we found that in vitro PDO responses to 5-FU, irinotecan, and GnP were correlated with patient responses based on RECIST criteria. Interestingly, and similar to previous reports in colorectal cancer PDOs, we found that oxaliplatin responses were not predictive of RECIST response. As expected, PDOs were resistant to afatinib (EGFRi), which reflects negative clinical trials, and may also be masked by use of EGF in growth media. A range of responses to trametinib (MEKi) were seen but were not correlated with KRAS allelic dosage. A similar range of response was seen to talazoparib (PARPi), but did not correlate with oxaliplatin response or HRDetect scores. In summary, we have established a novel biobank of PDOs from advanced pancreatic cancer patients. Notably, PDOs were less likely to establish from tumors that were KRAS WT or HR-deficient, even though these patients are likely to benefit from targeted approaches. Further investigation is required to develop PDO use in clinical drug prediction and drug discovery.

Citation Format: Irene Y. Xie, Laura Tamblyn, Karen Ng, Eugenia Flores-Figueroa, Julie M. Wilson, Gun Ho Jang, Amy X. Zhang, Stephanie Ramotar, Anna Dodd, Nikolina Radulovich, Jennifer J. Knox, Grainne M. O'Kane, Steven Gallinger, Faiyaz Notta. Establishment of a novel living biobank of patient-derived pancreatic cancer organoids with genomic and drug response characterization [abstract]. In: Proceedings of the AACR Virtual Special Conference on Pancreatic Cancer; 2021 Sep 29-30. Philadelphia (PA): AACR; Cancer Res 2021;81(22 Suppl):Abstract nr PO-077.

The value of GATA6 immunohistochemistry and computer-assisted diagnosis to predict clinical outcome in advanced pancreatic cancer

Combination chemotherapy, either modified FOLFIRINOX (mFFX) or gemcitabine-nabpaclitaxel, are used in the treatment of most patients with advanced pancreatic ductal adenocarcinoma (PDAC), yet robust biomarkers of outcome are currently lacking to guide regimen selection. Here, we tested GATA6 immunohistochemistry (IHC) as a putative biomarker in advanced PDAC. GATA6 is a transcription factor in normal pancreas development. Two pathologists, blinded to clinical and molecular data, independently assessed GATA6 IHC in biopsy specimens of 130 patients with advanced PDAC, in 2 distinct phases (without and with computer assistance using the open source software QuPath). Low GATA6 IHC expression was associated with shorter overall survival [median OS 6.2 months for patients with GATA6 low tumors vs. 11.5 months for patients with GATA6 high tumors, HR 1.66 (95% CI 1.15-2.40), P = 0.007]. Progression appears to be higher in GATA6-low tumors compared to GATA6-high tumors in patients treated with mFFX (P = 0.024) but not in patients treated with gemcitabine regimens. GATA6 IHC expression was significantly associated with molecular subtypes (P = 0.0003). Digital assistance markedly improved interrater concordance (Cohen's kappa scores of 0.32 vs. 0.95). Our results provide strong evidence that GATA6 IHC can be used as a single biomarker in the clinic to predict clinical outcome in advanced PDAC, warranting further investigation in prospective clinical trials. These results provide the basis for an improved classification of PDAC and future biomarker design using digital pathology workflow.

Patient-derived tumor xenograft and organoid models established from resected pancreatic, duodenal and biliary cancers

Patient-derived xenograft (PDX) and their xenograft-derived organoid (XDO) models that recapitulate the genotypic and phenotypic landscape of patient cancers could help to advance research and lead to improved clinical management. PDX models were established from 276 pancreato-duodenal and biliary cancer resections. Initial, passage 0 (P0) engraftment rates were 59% (118/199) for pancreatic, 86% (25/29) for duodenal, and 35% (17/48) for biliary ductal tumors. Pancreatic ductal adenocarcinoma (PDAC), had a P0 engraftment rate of 62% (105/169). KRAS mutant and wild-type PDAC models were molecularly profiled, and XDO models were generated to perform initial drug response evaluations. Subsets of PDAC PDX models showed global copy number variants and gene expression profiles that were retained with serial passaging, and they showed a spectrum of somatic mutations represented in patient tumors. PDAC XDO models were established, with a success rate of 71% (10/14). Pathway activation of KRAS-MAPK in PDXs was independent of KRAS mutational status. Four wild-type KRAS models were characterized by one with EGFR (L747-P753 del), two with BRAF alterations (N486_P490del or V600E), and one with triple negative KRAS/EGFR/BRAF. Model OCIP256, characterized by BRAF (N486-P490 del), had activated phospho-ERK. A combination treatment of a pan-RAF inhibitor (LY3009120) and a MEK inhibitor (trametinib) effectively suppressed phospho-ERK and inhibited growth of OCIP256 XDO and PDX models. PDAC/duodenal adenocarcinoma have high success rates forming PDX/organoid and retaining their phenotypic and genotypic features. These models may be effective tools to evaluate novel drug combination therapies.

Genomic Features and Classification of Homologous Recombination Deficient Pancreatic Ductal Adenocarcinoma

BACKGROUND AND AIMS

Homologous recombination deficiency (HRD) in pancreatic ductal adenocarcinoma (PDAC), remains poorly defined beyond germline (g) alterations in BRCA1, BRCA2, and PALB2.

METHODS

We interrogated whole genome sequencing (WGS) data on 391 patients, including 49 carriers of pathogenic variants (PVs) in gBRCA and PALB2. HRD classifiers were applied to the dataset and included (1) the genomic instability score (GIS) used by Myriad's MyChoice HRD assay; (2) substitution base signature 3 (SBS3); (3) HRDetect; and (4) structural variant (SV) burden. Clinical outcomes and responses to chemotherapy were correlated with HRD status.

RESULTS

Biallelic tumor inactivation of gBRCA or PALB2 was evident in 43 of 49 germline carriers identifying HRD-PDAC. HRDetect (score ≥0.7) predicted gBRCA1/PALB2 deficiency with highest sensitivity (98%) and specificity (100%). HRD genomic tumor classifiers suggested that 7% to 10% of PDACs that do not harbor gBRCA/PALB2 have features of HRD. Of the somatic HRDetect^hi cases, 69% were attributed to alterations in BRCA1/2, PALB2, RAD51C/D, and XRCC2, and a tandem duplicator phenotype. TP53 loss was more common in BRCA1- compared with BRCA2-associated HRD-PDAC. HRD status was not prognostic in resected PDAC; however in advanced disease the GIS (P = .02), SBS3 (P = .03), and HRDetect score (P = .005) were predictive of platinum response and superior survival. PVs in gATM (n = 6) or gCHEK2 (n = 2) did not result in HRD-PDAC by any of the classifiers. In 4 patients, BRCA2 reversion mutations associated with platinum resistance.

CONCLUSIONS

Germline and parallel somatic profiling of PDAC outperforms germline testing alone in identifying HRD-PDAC. An additional 7% to 10% of patients without gBRCA/PALB2 mutations may benefit from DNA damage response agents.

Vaccination with Detoxified Leukocidin AB Reduces Bacterial Load in a Staphylococcus aureus Minipig Deep Surgical Wound Infection Model

Vaccines against Staphylococcus aureus have eluded researchers for >3 decades while the burden of staphylococcal diseases has increased. Early vaccine attempts mainly used rodents to characterize preclinical efficacy, and all subsequently failed in human clinical efficacy trials. More recently, leukocidin AB (LukAB) has gained interest as a vaccine antigen. We developed a minipig deep surgical wound infection model offering 3 independent efficacy readouts: bacterial load at the superficial and at the deep-seated surgical site, and dissemination of bacteria. Due to similarities with humans, minipigs are an attractive option to study novel vaccine candidates. With this model, we characterized the efficacy of a LukAB toxoid as vaccine candidate. Compared to control animals, a 3-log reduction of bacteria at the deep-seated surgical site was observed in LukAB-treated minipigs and dissemination of bacteria was dramatically reduced. Therefore, LukAB toxoids may be a useful addition to S. aureus vaccines and warrant further study.

Coronary artery calcium scoring: an evidence-based guide for primary care physicians

Primary care physicians often must decide whether statin therapy would be appropriate (in addition to lifestyle modification) for managing asymptomatic individuals with borderline or intermediate risk for developing atherosclerotic cardiovascular disease (ASCVD), as assessed on the basis of traditional risk factors. In appropriate subjects, a simple, noninvasive measurement of coronary artery calcium can help clarify risk. Coronary atherosclerosis is a chronic inflammatory disease, with atherosclerotic plaque formation involving intimal inflammation and repeated cycles of erosion and fibrosis, healing and calcification. Atherosclerotic plaque formation represents the prognostic link between risk factors and future clinical events. The presence of coronary artery calcification is almost exclusively an indication of coronary artery disease, except in certain metabolic conditions. Coronary artery calcification can be detected and quantified in a matter of seconds by noncontrast electrocardiogram-gated low-dose X-ray computed tomography (coronary artery calcium scoring [CACS]). Since the publication of the seminal work by Dr. Arthur Agatston in 1990, a wealth of CACS-based prognostic data has been reported. In addition, recent guidelines from various professional societies conclude that CACS may be considered as a tool for reclassifying risk for atherosclerotic cardiovascular disease in patients otherwise assessed to have intermediate risk, so as to more accurately inform decisions about possible statin therapy in addition to lifestyle modification as primary preventive therapy. In this review, we provide an overview of CACS, from acquisition to interpretation, and summarize the scientific evidence for and the appropriate use of CACS as put forth in current clinical guidelines.

Subtype-Discordant Pancreatic Ductal Adenocarcinoma Tumors Show Intermediate Clinical and Molecular Characteristics

PURPOSE

RNA-sequencing-based subtyping of pancreatic ductal adenocarcinoma (PDAC) has been reported by multiple research groups, each using different methodologies and patient cohorts. "Classical" and "basal-like" PDAC subtypes are associated with survival differences, with basal-like tumors associated with worse prognosis. We amalgamated various PDAC subtyping tools to evaluate the potential of such tools to be reliable in clinical practice.

EXPERIMENTAL DESIGN

Sequencing data for 574 PDAC tumors was obtained from prospective trials and retrospective public databases. Six published PDAC subtyping strategies (Moffitt regression tools, clustering-based Moffitt, Collisson, Bailey, and Karasinska subtypes) were used on each sample, and results were tested for subtype call consistency and association with survival.

RESULTS

Basal-like and classical subtype calls were concordant in 88% of patient samples, and survival outcomes were significantly different (P < 0.05) between prognostic subtypes. Twelve percent of tumors had subtype-discordant calls across the different methods, showing intermediate survival in univariate and multivariate survival analyses. Transcriptional profiles compatible with that of a hybrid subtype signature were observed for subtype-discordant tumors, in which classical and basal-like genes were concomitantly expressed. Subtype-discordant tumors showed intermediate molecular characteristics, including subtyping gene expression (P < 0.0001) and mutant KRAS allelic imbalance (P < 0.001).

CONCLUSIONS

Nearly 1 in 6 patients with PDAC have tumors that fail to reliably fall into the classical or basal-like PDAC subtype categories, based on two regression tools aimed toward clinical practice. Rather, these patient tumors show intermediate prognostic and molecular traits. We propose close consideration of the non-binary nature of PDAC subtypes for future incorporation of subtyping into clinical practice.

Validation of Prognostic Radiomic Features From Resectable Pancreatic Ductal Adenocarcinoma in Patients With Advanced Disease Undergoing Chemotherapy

BACKGROUND

Radiomic features in pancreatic ductal adenocarcinoma (PDAC) often lack validation in independent test sets or are limited to early or late stage disease. Given the lethal nature of PDAC it is possible that there are similarities in radiomic features of both early and advanced disease reflective of aggressive biology.

PURPOSE

To assess the performance of prognostic radiomic features previously published in patients with resectable PDAC in a test set of patients with unresectable PDAC undergoing chemotherapy.

METHODS

The pre-treatment CT of 108 patients enrolled in a prospective chemotherapy trial were used as a test cohort for 2 previously published prognostic radiomic features in resectable PDAC (Sum Entropy and Cluster Tendency with square-root filter[Sqrt]). We assessed the performance of these 2 radiomic features for the prediction of overall survival (OS) and time to progression (TTP) using Cox proportional-hazard models.

RESULTS

Sqrt Cluster Tendency was significantly associated with outcome with a hazard ratio (HR) of 1.27(for primary pancreatic tumor plus local nodes), (Confidence Interval(CI):1.01 -1.6, P-value = 0.039) for OS and a HR of 1.25(CI:1.00 -1.55, P-value = 0.047) for TTP. Sum entropy was not associated with outcomes. Sqrt Cluster Tendency remained significant in multivariate analysis.

CONCLUSION

The CT radiomic feature Sqrt Cluster Tendency, previously demonstrated to be prognostic in resectable PDAC, remained a significant prognostic factor for OS and TTP in a test set of unresectable PDAC patients. This radiomic feature warrants further investigation to understand its biologic correlates and CT applicability in PDAC patients.

A Preclinical Trial and Molecularly Annotated Patient Cohort Identify Predictive Biomarkers in Homologous Recombination–deficient Pancreatic Cancer

Purpose:

Pancreatic ductal adenocarcinoma (PDAC) arising in patients with a germline BRCA1 or BRCA2 (gBRCA) mutation may be sensitive to platinum and PARP inhibitors (PARPi). However, treatment stratification based on gBRCA mutational status alone is associated with heterogeneous responses.

Experimental Design:

We performed a seven-arm preclinical trial consisting of 471 mice, representing 12 unique PDAC patient-derived xenografts, of which nine were gBRCA mutated. From 179 patients whose PDAC was whole-genome and transcriptome sequenced, we identified 21 cases with homologous recombination deficiency (HRD), and investigated prognostic biomarkers.

Results:

We found that biallelic inactivation of BRCA1/BRCA2 is associated with genomic hallmarks of HRD and required for cisplatin and talazoparib (PARPi) sensitivity. However, HRD genomic hallmarks persisted in xenografts despite the emergence of therapy resistance, indicating the presence of a genomic scar. We identified tumor polyploidy and a low Ki67 index as predictors of poor cisplatin and talazoparib response. In patients with HRD PDAC, tumor polyploidy and a basal-like transcriptomic subtype were independent predictors of shorter survival. To facilitate clinical assignment of transcriptomic subtype, we developed a novel pragmatic two-marker assay (GATA6:KRT17).

Conclusions:

In summary, we propose a predictive and prognostic model of gBRCA-mutated PDAC on the basis of HRD genomic hallmarks, Ki67 index, tumor ploidy, and transcriptomic subtype.

Delving into Early-onset Pancreatic Ductal Adenocarcinoma: How Does Age Fit In?

PURPOSE

With the rising incidence of early-onset pancreatic cancer (EOPC), molecular characteristics that distinguish early-onset pancreatic ductal adenocarcinoma (PDAC) tumors from those arising at a later age are not well understood.

EXPERIMENTAL DESIGN

We performed bioinformatic analysis of genomic and transcriptomic data generated from 269 advanced (metastatic or locally advanced) and 277 resectable PDAC tumor samples. Patient samples were stratified into EOPC (age of onset ≤55 years; n = 117), intermediate (age of onset 55-70 years; n = 264), and average (age of onset ≥70 years; n = 165) groups. Frequency of somatic mutations affecting genes commonly implicated in PDAC, as well as gene expression patterns, were compared between EOPC and all other groups.

RESULTS

EOPC tumors showed significantly lower frequency of somatic single-nucleotide variant (SNV)/insertions/deletions (indel) in CDKN2A (P = 0.0017), and were more likely to achieve biallelic mutation of CDKN2A through homozygous copy loss as opposed to heterozygous copy loss coupled with a loss-of-function SNV/indel mutation, the latter of which was more common for tumors with later ages of onset (P = 1.5e-4). Transcription factor forkhead box protein C2 (FOXC2) was significantly upregulated in EOPC tumors (P = 0.032). Genes significantly correlated with FOXC2 in PDAC samples were enriched for gene sets related to epithelial-to-mesenchymal transition (EMT) and included VIM (P = 1.8e-8), CDH11 (P = 6.5e-5), and CDH2 (P = 2.4e-2).

CONCLUSIONS

Our comprehensive analysis of sequencing data generated from a large cohort of PDAC patient samples highlights a distinctive pattern of biallelic CDKN2A mutation in EOPC tumors. Increased expression of FOXC2 in EOPC, with the correlation between FOXC2 and EMT pathways, represents novel molecular characteristics of EOPC.See related commentary by Lou, p. 8.

Author Correction: Transcription phenotypes of pancreatic cancer are driven by genomic events during tumor evolution

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

A Four-Chemokine Signature Is Associated with a T-cell-Inflamed Phenotype in Primary and Metastatic Pancreatic Cancer

PURPOSE

The molecular drivers of antitumor immunity in pancreatic ductal adenocarcinoma (PDAC) are poorly understood, posing a major obstacle for the identification of patients potentially amenable for immune-checkpoint blockade or other novel strategies. Here, we explore the association of chemokine expression with effector T-cell infiltration in PDAC.

EXPERIMENTAL DESIGN

Discovery cohorts comprised 113 primary resected PDAC and 107 PDAC liver metastases. Validation cohorts comprised 182 PDAC from The Cancer Genome Atlas and 92 PDACs from the Australian International Cancer Genome Consortium. We explored associations between immune cell counts by immunohistochemistry, chemokine expression, and transcriptional hallmarks of antitumor immunity by RNA sequencing (RNA-seq), and mutational burden by whole-genome sequencing.

RESULTS

Among all known human chemokines, a coregulated set of four (CCL4, CCL5, CXCL9, and CXCL10) was strongly associated with CD8⁺ T-cell infiltration (P < 0.001). Expression of this "4-chemokine signature" positively correlated with transcriptional metrics of T-cell activation (ZAP70, ITK, and IL2RB), cytolytic activity (GZMA and PRF1), and immunosuppression (PDL1, PD1, CTLA4, TIM3, TIGIT, LAG3, FASLG, and IDO1). Furthermore, the 4-chemokine signature marked tumors with increased T-cell activation scores (MHC I presentation, T-cell/APC costimulation) and elevated expression of innate immune sensing pathways involved in T-cell priming (STING and NLRP3 inflammasome pathways, BATF3-driven dendritic cells). Importantly, expression of this 4-chemokine signature was consistently indicative of a T-cell-inflamed phenotype across primary PDAC and PDAC liver metastases.

CONCLUSIONS

A conserved 4-chemokine signature marks resectable and metastatic PDAC tumors with an active antitumor phenotype. This could have implications for the appropriate selection of PDAC patients in immunotherapy trials.

Transcription phenotypes of pancreatic cancer are driven by genomic events during tumor evolution

Pancreatic adenocarcinoma presents as a spectrum of a highly aggressive disease in patients. The basis of this disease heterogeneity has proved difficult to resolve due to poor tumor cellularity and extensive genomic instability. To address this, a dataset of whole genomes and transcriptomes was generated from purified epithelium of primary and metastatic tumors. Transcriptome analysis demonstrated that molecular subtypes are a product of a gene expression continuum driven by a mixture of intratumoral subpopulations, which was confirmed by single-cell analysis. Integrated whole-genome analysis uncovered that molecular subtypes are linked to specific copy number aberrations in genes such as mutant KRAS and GATA6. By mapping tumor genetic histories, tetraploidization emerged as a key mutational process behind these events. Taken together, these data support the premise that the constellation of genomic aberrations in the tumor gives rise to the molecular subtype, and that disease heterogeneity is due to ongoing genomic instability during progression.

Application of multi-criteria decision analysis to assess sampling strategies in eutrophic urbanized waterbodies

Different water quality sampling practices such as location selection or frequency can inform future watershed management strategies. The objective of this work was to compare water quality sampling strategies based on different weighted criteria to determine the optimal sampling frequency and sampling location for an urbanized, eutrophic, freshwater system. Weekly water sampling was conducted over a 2-year period at five locations for six water quality parameters. This high frequency (HF) dataset was then deconstructed into a lower frequency (LF) dataset to simulate a monthly sampling strategy. Statistical analyses conducted showed that for all sampling locations the LF datasets were not significantly different from the HF datasets, suggesting monthly sampling is sufficient to capture the overall water quality conditions in this system. A multi-criteria decision analysis was constructed for statistical and operational criteria to determine the optimal sampling locations given different criteria weights. Results showed that the optimal sampling location changed depending on the criteria weighting, suggesting that statistical analyses alone would not be sufficient to determine optimal sampling locations in this system. This analysis was then used if optimal sampling location depended on specific water quality monitoring goals. Results showed that the optimal location depends on the particular water quality monitoring goals and that this effect should also be considered in the design of future sampling programs.

Integration of Genomic and Transcriptional Features in Pancreatic Cancer Reveals Increased Cell Cycle Progression in Metastases

We integrated clinical, genomic, and transcriptomic data from 224 primaries and 95 metastases from 289 patients to characterize progression of pancreatic ductal adenocarcinoma (PDAC). Driver gene alterations and mutational and expression-based signatures were preserved, with truncations, inversions, and translocations most conserved. Cell cycle progression (CCP) increased with sequential inactivation of tumor suppressors, yet remained higher in metastases, perhaps driven by cell cycle regulatory gene variants. Half of the cases were hypoxic by expression markers, overlapping with molecular subtypes. Paired tumor heterogeneity showed cancer cell migration by Halstedian progression. Multiple PDACs arising synchronously and metachronously in the same pancreas were actually intra-parenchymal metastases, not independent primary tumors. Established clinical co-variates dominated survival analyses, although CCP and hypoxia may inform clinical practice.

Whole genomes define concordance of matched primary, xenograft, and organoid models of pancreas cancer

Pancreatic ductal adenocarcinoma (PDAC) has the worst prognosis among solid malignancies and improved therapeutic strategies are needed to improve outcomes. Patient-derived xenografts (PDX) and patient-derived organoids (PDO) serve as promising tools to identify new drugs with therapeutic potential in PDAC. For these preclinical disease models to be effective, they should both recapitulate the molecular heterogeneity of PDAC and validate patient-specific therapeutic sensitivities. To date however, deep characterization of the molecular heterogeneity of PDAC PDX and PDO models and comparison with matched human tumour remains largely unaddressed at the whole genome level. We conducted a comprehensive assessment of the genetic landscape of 16 whole-genome pairs of tumours and matched PDX, from primary PDAC and liver metastasis, including a unique cohort of 5 'trios' of matched primary tumour, PDX, and PDO. We developed a pipeline to score concordance between PDAC models and their paired human tumours for genomic events, including mutations, structural variations, and copy number variations. Tumour-model comparisons of mutations displayed single-gene concordance across major PDAC driver genes, but relatively poor agreement across the greater mutational load. Genome-wide and chromosome-centric analysis of structural variation (SV) events highlights previously unrecognized concordance across chromosomes that demonstrate clustered SV events. We found that polyploidy presented a major challenge when assessing copy number changes; however, ploidy-corrected copy number states suggest good agreement between donor-model pairs. Collectively, our investigations highlight that while PDXs and PDOs may serve as tractable and transplantable systems for probing the molecular properties of PDAC, these models may best serve selective analyses across different levels of genomic complexity.

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.

Mutations in Mitochondrial DNA From Pancreatic Ductal Adenocarcinomas Associate With Survival Times of Patients and Accumulate as Tumors Progress

Somatic mutations have been found in the mitochondria in different types of cancer cells, but it is not clear whether these affect tumorigenesis or tumor progression. We analyzed mitochondrial genomes of 268 early-stage, resected pancreatic ductal adenocarcinoma tissues and paired non-tumor tissues. We defined a mitochondrial somatic mutation (mtSNV) as a position where the difference in heteroplasmy fraction between tumor and normal sample was ≥0.2. Our analysis identified 304 mtSNVs, with at least 1 mtSNV in 61% (164 of 268) of tumor samples. The noncoding control region had the greatest proportion of mtSNVs (60 of 304 mutations); this region contains sites that regulate mitochondrial DNA transcription and replication. Frequently mutated genes included ND5, RNR2, and CO1, plus 29 mutations in transfer RNA genes. mtSNVs in 2 separate mitochondrial genes (ND4 and ND6) were associated with shorter overall survival time. This association appeared to depend on the level of mtSNV heteroplasmy. Non-random co-occurrence between mtSNVs and mutations in nuclear genes indicates interactions between nuclear and mitochondrial DNA. In an analysis of primary tumors and metastases from 6 patients, we found tumors to accumulate mitochondrial mutational mutations as they progress.

Recapitulating the clinical scenario of BRCA-associated pancreatic cancer in pre-clinical models

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal malignancies. BRCA-associated PDAC comprises a clinically relevant subtype. A portion of these patients are highly susceptible to DNA damaging therapeutics, however, responses are heterogeneous and clinical resistance evolves. We have developed unique patient-derived xenograft (PDX) models from metastatic lesions of germline BRCA-mutated patients obtained at distinct time points; before treatment and at progression. Thus, closely mimicking clinical scenarios, to further investigate treatment naïve and resistant patients. DNA was isolated from six BRCA-mutated PDXs and classified by whole-genome sequencing to stable-genome or homologous recombination deficient (HRD)-genome. The sensitivity to DNA-damaging agents was evaluated in vivo in three BRCA-associated PDAC PDXs models: (1) HRD-genome naïve to treatments; (2) stable-genome naïve to treatment; (3) HRD-genome resistant to treatment. Correlation between disease course at tissue acquisition and response to PARP inhibitor (PARPi)/platinum was demonstrated in PDXs in vivo. Only the HRD-genome PDX, naïve to treatment, was sensitive to PARP inhibitor/cisplatin treatments. Our results demonstrate heterogeneous responses to DNA damaging agents/PARPi in BRCA-associated PDX thus reflecting the wide clinical spectrum. An HRD-genome PDX generated from a naïve to treatment biopsy was sensitive to platinum/PARPi whereas no benefit was observed in treating a HRD-genome PDXs generated from a patient that had acquired resistance nor stable-genome PDXs.

Case Report: First Reported Combined Heart-Liver Transplant in a Patient With a Congenital Solitary Kidney

We report a case of successful combined heart liver transplant in a patient with a congenital solitary kidney. The patient had normal renal function before combined heart-liver transplantation and developed acute kidney injury requiring slow continuous dialysis and subsequent intermittent dialysis for almost 8 weeks post transplantation. Her renal function recovered and she remains off dialysis now 7 months post transplantation. She only currently has mild chronic renal insufficiency. We believe this is the first reported case of successful heart liver transplant in a patient with a congenital solitary kidney.

Exome-Wide Association Study of Pancreatic Cancer Risk

We conducted a case-control exome-wide association study to discover germline variants in coding regions that affect risk for pancreatic cancer, combining data from 5 studies. We analyzed exome and genome sequencing data from 437 patients with pancreatic cancer (cases) and 1922 individuals not known to have cancer (controls). In the primary analysis, BRCA2 had the strongest enrichment for rare inactivating variants (17/437 cases vs 3/1922 controls) (P = 3.27x10-6; exome-wide statistical significance threshold P < 2.5x10-6). Cases had more rare inactivating variants in DNA repair genes than controls, even after excluding 13 genes known to predispose to pancreatic cancer (adjusted odds ratio, 1.35; P = .045). At the suggestive threshold (P < .001), 6 genes were enriched for rare damaging variants (UHMK1, AP1G2, DNTA, CHST6, FGFR3, and EPHA1) and 7 genes had associations with pancreatic cancer risk, based on the sequence-kernel association test. We confirmed variants in BRCA2 as the most common high-penetrant genetic factor associated with pancreatic cancer and we also identified candidate pancreatic cancer genes. Large collaborations and novel approaches are needed to overcome the genetic heterogeneity of pancreatic cancer predisposition.

Genomics-Driven Precision Medicine for Advanced Pancreatic Cancer: Early Results from the COMPASS Trial

Purpose: 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.Experimental Design: 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.Results: 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 (P = 0.004). The best progression-free survival was observed in those with classical subtype treated with m-FOLFIRINOX. GATA6 expression in tumor measured by RNA in situ 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.Conclusions: Prospective genomic profiling of advanced PDAC is feasible, and our early data indicate that chemotherapy response differs among patients with different genomic/transcriptomic subtypes. Clin Cancer Res; 24(6); 1344-54. ©2017 AACR.

Cortisol plays a role in the high environmental ammonia associated suppression of the immune response in zebrafish

Exposure to high environmental ammonia (HEA) levels increases the vulnerability of fishes to parasitic, viral and bacterial diseases. We tested the hypothesis that elevated plasma cortisol levels play a role in the HEA-mediated immunosuppression in fishes. To this end, we tested the effect of exogenous cortisol treatment on the lipopolysaccharide (LPS)-induced immune response in zebrafish (Danio rerio). Also, to test whether glucocorticoid receptor (GR) signaling is involved in HEA-mediated immunosuppression, zebrafish were treated with mifepristone, a GR antagonist, and the LPS-induced immune response assessed after HEA exposure. We evaluated a panel of important immunity-related genes including interleukin 1β (il1b) and suppressor of cytokine signaling (socs-1a, 2, 3) and acute phase response genes [serum amyloid A (saa), transferrin (tfa), leukocyte cell-derived chemotaxin 2-like (lect2l), haptoglobin (hp), hepcidin (=hepatic anti-microbial peptide hamp), and complement component 3b (c3b)] by real-time quantitative PCR. Our results demonstrate that exogenous cortisol administration as well as elevated cortisol levels in response to HEA exposure modulate mRNA transcript levels of key mediators of the innate immune response in zebrafish. Mifepristone treatment reduced whole body cortisol levels and eliminated the HEA-mediated changes in transcript abundance of socs1a, il1b, as well as APR genes. Together, these results suggest that the HEA effect on the innate immune response is in part mediated by cortisol signaling, while the mode of action, including the receptors involved remains to be elucidated.

Correlation of 18F-Fluorodeoxyglucose Positron Emission Tomography Parameters with Patterns of Disease Progression in Locally Advanced Pancreatic Cancer after Definitive Chemoradiotherapy

AIMS

A proportion of patients with pancreatic cancer never develop metastatic disease. We evaluated a role for 18F-fluorodeoxyglucose positron emission tomography (FDG-PET) in identifying a subset of patients with locally advanced pancreatic cancer (LAPC) who never develop metastatic disease and only experience local disease and may therefore benefit from local treatment intensification.

MATERIAL AND METHODS

Patients with histologically confirmed LAPC entered a single-centre phase II study of definitive upfront chemoradiotherapy (CRT). All patients underwent FDG-PET/CT before and 6 weeks after CRT. Tumour volume, standardised uptake values (SUVmax, SUVpeak, SUVmean, SUVmedian) and total lesion glycolysis (TLG) were measured on each scan and the response in each parameter was evaluated. The presence or absence of metastatic disease was noted on contrast-enhanced CT carried out every 3 months for 1 year and then at clinician discretion.

RESULTS

Twenty-three patients with LAPC were recruited; 17/23 completed treatment and had interpretable sequential imaging. Twenty-four per cent of patients only ever experienced local disease. Median pre-CRT FDG-PET parameters were significantly lower in patients with local disease only during follow-up compared with those who developed metastatic disease: SUVmax 3.8 versus 8.6 (P=0.006), SUVpeak 2.5 versus 7.5 (P=0.002), SUVmean 1.8 versus 3.3 (P=0.001), SUVmedian 1.7 versus 3.0 (P=0.002), TLG 26.9 versus 115.9 (P=0.006). Tumour volume, post-CRT FDG-PET values and their relative change were not statistically different between local disease and metastatic disease groups. Receiver operating characteristic curves for pre-CRT FDG-PET parameters to predict those who never develop metastatic disease all had areas under the curve (AUCs) ≥ 0.932. Pre-CRT FDG-PET SUVmax < 6.2 predicted patients with local disease only during follow-up with 100.0% sensitivity and 92.3% specificity, 80.0% positive predictive value and 100% negative predictive value.

CONCLUSIONS

Our findings suggest that patients with less FDG-avid tumours are less likely to metastasise and may therefore benefit from upfront local treatment intensification.

Association of Distinct Mutational Signatures With Correlates of Increased Immune Activity in Pancreatic Ductal Adenocarcinoma

IMPORTANCE

Outcomes for patients with pancreatic ductal adenocarcinoma (PDAC) remain poor. Advances in next-generation sequencing provide a route to therapeutic approaches, and integrating DNA and RNA analysis with clinicopathologic data may be a crucial step toward personalized treatment strategies for this disease.

OBJECTIVE

To classify PDAC according to distinct mutational processes, and explore their clinical significance.

DESIGN, SETTING, AND PARTICIPANTS

We performed a retrospective cohort study of resected PDAC, using cases collected between 2008 and 2015 as part of the International Cancer Genome Consortium. The discovery cohort comprised 160 PDAC cases from 154 patients (148 primary; 12 metastases) that underwent tumor enrichment prior to whole-genome and RNA sequencing. The replication cohort comprised 95 primary PDAC cases that underwent whole-genome sequencing and expression microarray on bulk biospecimens.

MAIN OUTCOMES AND MEASURES

Somatic mutations accumulate from sequence-specific processes creating signatures detectable by DNA sequencing. Using nonnegative matrix factorization, we measured the contribution of each signature to carcinogenesis, and used hierarchical clustering to subtype each cohort. We examined expression of antitumor immunity genes across subtypes to uncover biomarkers predictive of response to systemic therapies.

RESULTS

The discovery cohort was 53% male (n = 79) and had a median age of 67 (interquartile range, 58-74) years. The replication cohort was 50% male (n = 48) and had a median age of 68 (interquartile range, 60-75) years. Five predominant mutational subtypes were identified that clustered PDAC into 4 major subtypes: age related, double-strand break repair, mismatch repair, and 1 with unknown etiology (signature 8). These were replicated and validated. Signatures were faithfully propagated from primaries to matched metastases, implying their stability during carcinogenesis. Twelve of 27 (45%) double-strand break repair cases lacked germline or somatic events in canonical homologous recombination genes-BRCA1, BRCA2, or PALB2. Double-strand break repair and mismatch repair subtypes were associated with increased expression of antitumor immunity, including activation of CD8-positive T lymphocytes (GZMA and PRF1) and overexpression of regulatory molecules (cytotoxic T-lymphocyte antigen 4, programmed cell death 1, and indolamine 2,3-dioxygenase 1), corresponding to higher frequency of somatic mutations and tumor-specific neoantigens.

CONCLUSIONS AND RELEVANCE

Signature-based subtyping may guide personalized therapy of PDAC in the context of biomarker-driven prospective trials.

A renewed model of pancreatic cancer evolution based on genomic rearrangement patterns

Pancreatic cancer, a highly aggressive tumour type with uniformly poor prognosis, exemplifies the classically held view of stepwise cancer development. The current model of tumorigenesis, based on analyses of precursor lesions, termed pancreatic intraepithelial neoplasm (PanINs) lesions, makes two predictions: first, that pancreatic cancer develops through a particular sequence of genetic alterations (KRAS, followed by CDKN2A, then TP53 and SMAD4); and second, that the evolutionary trajectory of pancreatic cancer progression is gradual because each alteration is acquired independently. A shortcoming of this model is that clonally expanded precursor lesions do not always belong to the tumour lineage, indicating that the evolutionary trajectory of the tumour lineage and precursor lesions can be divergent. This prevailing model of tumorigenesis has contributed to the clinical notion that pancreatic cancer evolves slowly and presents at a late stage. However, the propensity for this disease to rapidly metastasize and the inability to improve patient outcomes, despite efforts aimed at early detection, suggest that pancreatic cancer progression is not gradual. Here, using newly developed informatics tools, we tracked changes in DNA copy number and their associated rearrangements in tumour-enriched genomes and found that pancreatic cancer tumorigenesis is neither gradual nor follows the accepted mutation order. Two-thirds of tumours harbour complex rearrangement patterns associated with mitotic errors, consistent with punctuated equilibrium as the principal evolutionary trajectory. In a subset of cases, the consequence of such errors is the simultaneous, rather than sequential, knockout of canonical preneoplastic genetic drivers that are likely to set-off invasive cancer growth. These findings challenge the current progression model of pancreatic cancer and provide insights into the mutational processes that give rise to these aggressive tumours.

Abstract 3124: Prognostic biomarkers for pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) has the lowest 5-year-survival rate of common cancers (&lt; 7%), but despite intensive research efforts over the past several decades, its dismal prognosis has barely improved. Recent genome-wide association studies (GWAS) have linked several genetic factors with pancreatic cancer, but there is generally poor overlap in the results from such studies, presumably due to population heterogeneity. In this study, we collected previously published pancreatic risk loci from multiple GWAS studies, and correlated these loci with clinical information gathered on a well-characterized cohort of 148 PDAC patients from the University Health Network in Toronto Canada, and several other North American hospitals, which have undergone whole-genome sequencing. Using this cohort, we attempted to identify prognostic genomic biomarkers by correlating the germline genomic alterations observed in our cohort to overall survival (OS). By understanding the mechanism of germline variants that alter OS, we hope to develop insights that will lead to improved detection and therapy for PDAC patients.

Among the 67 published risk loci we tested using a multivariate Cox proportional hazard model, we found a strong positive association of the single nucleotide polymorphism rs4785367 (RefSNP alleles: C/T on forward strand; MAF = 0.474) with overall survival of PDAC donors: HR = 0.426; CI = 0.268 - 0.686; p-value = 0.00029. A more detailed analysis at the genotype level revealed that the presence of the homozygous minor allele has a stronger effect than either the heterozygous or homozygous major allele. The SNP falls within the intergenic region between the ZNF423 and TMEM188 genes, within the exon 2 of lncRNA RP11-305A4.3 and overlaps a CTCF regulatory domain. Preliminary gene expression analysis from RNA sequencing data on a subset of PDAC donors (n = 28) shows that patients carrying the minor allele have significant higher TMEM188 expression than of the major type (p-value = 0.012), suggesting that this allele may influence the course of PDAC via TMEM188 activity. A recent study linked the gene product to activation of NK cells, which in turns increases the defense mechanism against the pathogens, infections and transformed tumors. These findings suggest a possible molecular mechanism influencing the course of PDAC. We are also exploring the effect of the presence of the minor allele on the regulatory CTCF region, by applying an integrative pipeline for risk SNP analysis to pancreatic cancer. This will possibly detect the effect on the NANOG motif binding and/or on CTCF looping.

In summary, the present study detected the rs4785367 as a prognostic biomarker for pancreatic cancer, with the novelty of increased TMEM188 gene expression being linked to the presence of the alternate allele in PDAC patients. Further investigations on this and on assessing the effect of the polymorphism on the regulatory CTCF feature are in progress.

Citation Format: Cristina Baciu, Robert Grant, Hansen He, Musaddeque Ahmed, Robert E. Denroche, Lee Timms, Gun Ho Jang, Ayelet Borgida, Xihui Lin, Paul C. Boutros, Dianne Chadwich, Sheng-Ben Liang, Sagedeh Shahabi, Michael H.A. Roehrl, Sean Cleary, Julie M. Wilson, John D. McPherson, Lincoln Stein, Steven Gallinger. Prognostic biomarkers for pancreatic cancer. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 3124.