Ámundadóttir, Laufey T.

Abstract Pancreatic ductal adenocarcinoma (PDAC), as with all cancer, is driven by dysregulation of genetic programs due to accumulation of somatic mutations, epigenetic modifications and changes in the micro-environment. Holistic views of perturbations in gene expression networks could illuminate key regulators and pathways for this highly lethal cancer. Toward this end, we performed massively-parallel mRNA-sequencing in normal (n = 10) and tumor (n = 8) derived pancreatic tissue samples as well as pancreatic cancer cell lines (n = 9), and determined differential gene expression (DE) patterns. Sub-network enrichment analyses of all expressed genes based on magnitude of DE identified HNF1A as the regulator of the most significantly and consistently dysregulated expression sub-network in our pancreatic tumor samples (median P = 7.56 x 10-7, median rank = 1, range = 1-25). To explore the effect of HNF1A expression in pancreatic tumor-derived cells we generated stable HNF1A-inducible clones in two pancreatic cancer cell lines. Induction of HNF1A overexpression caused severe growth inhibition (5.3-fold, P = 4.5x10-5 for MIA PaCa-2 clones; 7.2-fold, P = 2.2x10-5 for PANC-1 clones), and G0/G1 cell cycle arrest. This was accompanied by down-regulation of 49 out of 84 assayed cell cycle genes, while only 1 displayed increased expression. These data, combined with HNF1A’s direct regulation of pancreatic development and homeostasis genes (i.e. PDX1, PTF1A and NR5A2), suggest it may be an important tumor suppressor in pancreatic cells. Citation Format: Jason W. Hoskins, Marta Flandez, Jinping Jia, Hemang Parikh, Irene Collins, Mickey Emmanuel, Abdisamad Ibrahim, Wenming Xiao, John Powell, Nuria Malats, Gloria M. Petersen, Fransisco X. Real, Laufey T. Amundadottir. Transcriptome analysis in pancreatic cancer reveals a tumor suppressor function for HNF1A. [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer: Innovations in Research and Treatment; May 18-21, 2014; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2015;75(13 Suppl):Abstract nr B112.

Abstract A genome-wide association study (GWAS) of pancreatic cancer conducted within the NCI Cohort Consortium (PanScan I and II) identified pancreatic cancer susceptibility loci on chromosomes 1q32.1/NR5A2, 5p15.33/CLPTM1L/TERT, 9q34.2/ABO, and 13q22.1. The most significant single-nucleotide polymorphism (SNP) identified on 13q22.1, rs9543325, lies in a 600 kb gene desert; the nearest genes are KLF5, KLF12, PIBF1, DIS3, and BORA (265–586 kb). Imputation using the 1000 Genomes and DCEG reference datasets did not improve the GWAS signal, but produced a set of highly correlated SNPs for functional follow-up. We performed eQTL analysis to test for association between the genotypes of these functional candidate variants and expression of nearby genes. Among 100 normal derived pancreatic tissue samples, DIS3 expression showed the strongest association with a novel 30 bp indel variant in the risk locus (P = 4.0 × 104), indicating risk alleles associate with reduced DIS3 expression. Mutations in DIS3 have previously been identified in acute myeloid leukemia and multiple myeloma, and its expression has been correlated with metastatic potential in colorectal cancer, suggesting DIS3 is relevant to cancer biology. Chromosome conformation capture identified a physical interaction between the indel-containing locus and a region near the DIS3 promoter. Luciferase assay for regulatory function of this indel-containing locus revealed allele-specific silencer activity for the insertion allele. Supershift electromobility shift assay (EMSA) demonstrated binding of LEF1 specifically to the insertion allele of the indel, which contains two in silico predicted LEF1 binding elements. Finally, through immunohistochemical analysis, high DIS3 protein levels associated with better survival for pancreatic cancer patients (hazard ratio = 2.87, 95% CI = 1.49–5.53, P = 0.001). Our results suggest that at least one target gene for the pancreatic cancer risk variants on chr13q22.1 may be DIS3, and that the underlying biology may be mediated by the novel indel through a long-range repressive effect on DIS3 expression. Citation Format: Jason W. Hoskins, Abdisamad Ibrahim, Mickey Emmanuel, Hemang Parikh, Jinping Jia, Irene Collins, Kris Ylaya, Sean F. Altekruse, Stephen M. Hewitt, Gloria M. Petersen, Laufey T. Amundadottir. Functional analysis of the chr13q22.1 pancreatic cancer risk locus suggests allele-specific effects on DIS3 expression with prognostic implications. [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer: Innovations in Research and Treatment; May 18-21, 2014; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2015;75(13 Suppl):Abstract nr B111.

Abstract Genome wide association studies (GWAS) have mapped multiple independent cancer risk loci (n = 6) to a small region on chr5p15.33 for at least ten distinct cancers, including bladder, breast, glioma, lung, melanoma, non-melanoma skin, ovarian, pancreas, prostate, and testicular germ cell cancer. This region harbors two plausible target genes, TERT which encodes the catalytic subunit of telomerase reverse transcriptase which maintains chromosome ends by adding telomeres repeats, and CLPTM1L which encodes the cleft lip and palate transmembrane protein 1-like protein which promotes cancer cell growth, protect cells from apoptosis, and can induce abnormal cytokinesis. The most significant pancreatic cancer GWAS SNP on chr5p15.33 was rs401681 (P = 3.7×10−7, OR = 1.19), located in the 13th intron of the CLPTM1L gene and within the independent multi-cancer risk locus in this region most distant from TERT. Imputation and fine mapping in pancreatic cancer resolved this signal by three orders of magnitude (P = 1.4×10−10, OR = 1.30) to a set of seven highly correlated SNPs. Electrophoretic mobility shift assays (EMSA) and luciferase assays across these SNPs led us to a single SNP, rs36115365, that showed allele-specific effects in both assays and overlapped with prominent ENCODE marks indicating regulatory potential. Examining this SNP across cell lines derived from multiple types of cancer linked to this locus, specifically pancreatic, testicular, lung, and melanoma, established a consistent pattern of a specific protein binding and enhanced regulatory activity for the risk allele across eight cancer cell lines (two cell lines per cancer type). siRNA-based targeting of the putative gene regulatory-region harboring rs36115365 resulted in a strong inhibition of TERT expression (average 60%, range 45-85%) but no effect on CLPTM1L expression. This inhibition of TERT expression was stronger from the pancreatic cancer risk as compared to the protective allele. Quantitative mass spectrometry identified Zinc finger protein 148 (ZNF148) as the protein specifically binding the risk allele, further supported by direct and allele-specific binding of purified recombinant ZNF148 in EMSA assays. Lastly, siRNA-mediated knockdown of ZNF148 resulted in significant reduction of TERT but not CLPTM1L expression across multiple cancer cell lines. Our results indicate that the pancreatic cancer GWAS signal at the single independent chr5p15.33 multi-cancer risk locus most distant from TERT may be explained by a single SNP that influences TERT expression via ZNF148. Citation Format: Kevin M. Brown, Jun Fang, Jinping Jia, Zhaoming Wang, Matthew Makowski, Tongwu Zhang, Jason Hoskins, Jiyeon Choi, Younghun Han, Mingfeng Zhang, Mai Xu, Peter Kanetsky, Andresson Thorkell, Gloria M. Petersen, Katherine L. Nathanson, Christopher I. Amos, Maria T. Landi, Stephen J. Chanock, Michiel Vermeulen, Laufey T. Amundadottir. Functional characterization of a multicancer risk locus on chr5p15.33 reveals regulation of TERT by ZNF148. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4610. doi:10.1158/1538-7445.AM2015-4610

A small number of common susceptibility loci have been identified for pancreatic cancer, one of which is marked by rs401681 in the TERT–CLPTM1L gene region on chromosome 5p15.33. Because this region is characterized by low linkage disequilibrium, we sought to identify whether additional single nucleotide polymorphisms (SNPs) could be related to pancreatic cancer risk, independently of rs401681. We performed an in‐depth analysis of genetic variability of the telomerase reverse transcriptase (TERT) and the telomerase RNA component (TERC) genes, in 5,550 subjects with pancreatic cancer and 7,585 controls from the PANcreatic Disease ReseArch (PANDoRA) and the PanScan consortia. We identified a significant association between a variant in TERT and pancreatic cancer risk (rs2853677, odds ratio = 0.85; 95% confidence interval = 0.80–0.90, p = 8.3 × 10−8). Additional analysis adjusting rs2853677 for rs401681 indicated that the two SNPs are independently associated with pancreatic cancer risk, as suggested by the low linkage disequilibrium between them (r2 = 0.07, D′ = 0.28). Three additional SNPs in TERT reached statistical significance after correction for multiple testing: rs2736100 (p = 3.0 × 10−5), rs4583925 (p = 4.0 × 10−5) and rs2735948 (p = 5.0 × 10−5). In conclusion, we confirmed that the TERT locus is associated with pancreatic cancer risk, possibly through several independent variants.

Abstract Pancreatic cancer is the 10th most common cancer and 4th most common cause of cancer mortality in the United States. A genome wide association study has revealed pancreatic cancer susceptibility loci on chromosomes 13q22.1, 1q32.1 and 5p15.33. The chr13q22.1 region contains the most significant of all risk SNPs, rs9543325, (P= 3.27 x 10−11) which is located in a gene desert. The nearest genes are KLF5, KLF12, PIBF1, DIS3 and BORA, which range in distance from 265kb to 586kb, respectively, from the most significant variants. Given the high linkage disequilibrium across this risk locus, there are many candidate functional variants to consider. Imputation in the region did not improve the signal but gave a set of highly correlated SNPs which likely includes the functional variant(s). In an effort to identify the functional variant(s) we performed eQTL analyses to test the association between the genotypes of these candidate SNPs and expression of nearby genes. Among 100 normal pancreatic tissue samples, DIS3 showed the strongest association with SNPs in our risk locus (P-values as low as 0.0004). Mutations in DIS3 have been identified in acute myeloid leukemia and multiple myeloma, and its expression has been correlated with metastatic potential in colorectal cancer, suggesting this gene could be important in pancreatic cancer biology. Chromosome Conformation Capture (3C) was performed to test for physical interactions between the risk locus and nearby genes. This assay confirmed the three dimensional proximity of the risk locus and DIS3 promoter. Sub-regions of the risk locus were then cloned upstream of a minimal promoter controlling luciferase expression to assay for potential allele specific enhancer/silencer activity. This assay revealed a sub-region downstream of rs9543325 that causes allele-specific silencing. These results suggest that a sub-region of the chr13q22.1 risk locus has allele-specific effects on DIS3 expression, which may have implications in the susceptibility to pancreatic cancer. Citation Format: Jason Hoskins, Abdisamad Ibrahim, Jinping Jia, Irene Collins, Hemang Parikh, Gloria M. Petersen, Laufey Amundadottir. Functional analysis of the chr13q22.1 pancreatic cancer risk locus suggests allele-specific effects on DIS3 expression. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 3137. doi:10.1158/1538-7445.AM2013-3137

Abstract Genome-wide association studies (GWAS) of 10 different cancers have identified pleiotropic cancer predisposition loci across a region of chromosome 5p15.33 that includes the TERT and CLPTM1L genes. Of these, susceptibility alleles for pancreatic cancer have mapped to the CLPTM1L gene, thus prompting an investigation of the function of CLPTM1L in the pancreas. Immunofluorescence analysis indicated that CLPTM1L localized to the endoplasmic reticulum where it is likely embedded in the membrane, in accord with multiple predicted transmembrane domains. Overexpression of CLPTM1L enhanced growth of pancreatic cancer cells in vitro (1.3–1.5–fold; PDAY7 < 0.003) and in vivo (3.46-fold; PDAY68 = 0.039), suggesting a role in tumor growth; this effect was abrogated by deletion of two hydrophilic domains. Affinity purification followed by mass spectrometry identified an interaction between CLPTM1L and non-muscle myosin II (NMM-II), a protein involved in maintaining cell shape, migration, and cytokinesis. The two proteins colocalized in the cytoplasm and, after treatment with a DNA-damaging agent, at the centrosomes. Overexpression of CLPTM1L and depletion of NMM-II induced aneuploidy, indicating that CLPTM1L may interfere with normal NMM-II function in regulating cytokinesis. Immunohistochemical analysis revealed enhanced staining of CLPTM1L in human pancreatic ductal adenocarcinoma (n = 378) as compared with normal pancreatic tissue samples (n = 17; P = 1.7 × 10−4). Our results suggest that CLPTM1L functions as a growth-promoting gene in the pancreas and that overexpression may lead to an abrogation of normal cytokinesis, indicating that it should be considered as a plausible candidate gene that could explain the effect of pancreatic cancer susceptibility alleles on chr5p15.33. Cancer Res; 74(10); 2785–95. ©2014 AACR.

Abstract Background: Melanoma cases may exist in pancreatic cancer kindreds, whereas there is increased risk of pancreatic cancer in familial melanoma. The two cancers may share genetic susceptibility variants in common. Methods: Three dbGaP (datasets in Genotypes and Phenotypes)-deposited GWAS (genome-wide association study) datasets (MD Anderson melanoma, PanScan 1, and PanScan 2 for pancreatic cancer) were used. Thirty-seven melanoma susceptibility variants in 22 genomic regions from published GWAS, plus melanoma-related genes and pathways were examined for pancreatic cancer risk in the PanScan datasets. Conversely, nine known pancreatic cancer susceptibility variants were examined for melanoma risk in the MD Anderson dataset. Results: In the PanScan data, initial associations were found with melanoma susceptibility variants in NCOA6 [rs4911442; OR, 1.32; 95% confidence interval (CI), 1.03–1.70; P = 0.03], YWHAZP5 (rs17119461; OR, 2.62; 95% CI, 1.08–6.35; P = 0.03), and YWHAZP5 (rs17119490; OR, 2.62; 95% CI, 1.08–6.34; P = 0.03), TYRP1 (P = 0.04), and IFNA13 (P = 0.04). In the melanoma dataset, two pancreatic cancer susceptibility variants were associated: NR5A2 (rs12029406; OR, 1.39; 95% CI, 1.01–1.92; P = 0.04) and CLPTM1L-TERT (rs401681; OR, 1.16; 95% CI, 1.01–1.34; P = 0.04). None of these associations remained significant after correcting for multiple comparisons. Conclusion: Reported variants of melanoma genes and pathways do not play a role in pancreatic cancer predisposition. Reciprocally, pancreatic cancer susceptibility variants are not associated with melanoma risk. Impact: Known melanoma-related genes and pathways, as well as GWAS-derived susceptibility variants of melanoma and pancreatic cancer, do not explain the shared genetic etiology of these two cancers. Cancer Epidemiol Biomarkers Prev; 23(6); 1121–4. ©2014 AACR.

Background and objective Survival of patients with pancreatic adenocarcinoma is limited and few prognostic factors are known. We conducted a two-stage genome-wide association study (GWAS) to identify germline variants associated with survival in patients with pancreatic adenocarcinoma. Methods We analysed overall survival in relation to single nucleotide polymorphisms (SNPs) among 1005 patients from two large GWAS datasets, PanScan I and ChinaPC. Cox proportional hazards regression was used in an additive genetic model with adjustment for age, sex, clinical stage and the top four principal components of population stratification. The first stage included 642 cases of European ancestry (PanScan), from which the top SNPs (p≤10 −5 ) were advanced to a joint analysis with 363 additional patients from China (ChinaPC). Results In the first stage of cases of European descent, the top-ranked loci were at chromosomes 11p15.4, 18p11.21 and 1p36.13, tagged by rs12362504 (p=1.63×10 −7 ), rs981621 (p=1.65×10 −7 ) and rs16861827 (p=3.75×10 −7 ), respectively. 131 SNPs with p≤10 −5 were advanced to a joint analysis with cases from the ChinaPC study. In the joint analysis, the top-ranked SNP was rs10500715 (minor allele frequency, 0.37; p=1.72×10 −7 ) on chromosome 11p15.4, which is intronic to the SET binding factor 2 (SBF2) gene. The HR (95% CI) for death was 0.74 (0.66 to 0.84) in PanScan I, 0.79 (0.65 to 0.97) in ChinaPC and 0.76 (0.68 to 0.84) in the joint analysis. Conclusions Germline genetic variation in the SBF2 locus was associated with overall survival in patients with pancreatic adenocarcinoma of European and Asian ancestry. This association should be investigated in additional large patient cohorts.