Amundadottir, Laufey

Objective To elucidate the genetic architecture of gene expression in pancreatic tissues. Design We performed expression quantitative trait locus (eQTL) analysis in histologically normal pancreatic tissue samples (n=95) using RNA sequencing and the corresponding 1000 genomes imputed germline genotypes. Data from pancreatic tumour-derived tissue samples (n=115) from The Cancer Genome Atlas were included for comparison. Results We identified 38 615 cis -eQTLs (in 484 genes) in histologically normal tissues and 39 713 cis -eQTL (in 237 genes) in tumour-derived tissues (false discovery rate <0.1), with the strongest effects seen near transcriptional start sites. Approximately 23% and 42% of genes with significant cis -eQTLs appeared to be specific for tumour-derived and normal-derived tissues, respectively. Significant enrichment of cis -eQTL variants was noted in non-coding regulatory regions, in particular for pancreatic tissues (1.53-fold to 3.12-fold, p≤0.0001), indicating tissue-specific functional relevance. A common pancreatic cancer risk locus on 9q34.2 (rs687289) was associated with ABO expression in histologically normal (p=5.8×10 −8 ) and tumour-derived (p=8.3×10 −5 ) tissues. The high linkage disequilibrium between this variant and the O blood group generating deletion variant in ABO (exon 6) suggested that nonsense-mediated decay (NMD) of the ‘O’ mRNA might explain this finding. However, knockdown of crucial NMD regulators did not influence decay of the ABO ‘O’ mRNA, indicating that a gene regulatory element influenced by pancreatic cancer risk alleles may underlie the eQTL. Conclusions We have identified cis -eQTLs representing potential functional regulatory variants in the pancreas and generated a rich data set for further studies on gene expression and its regulation in pancreatic tissues.

Abstract Individuals from pancreatic cancer families are at increased risk, not only of pancreatic cancer, but also of melanoma, breast, ovarian, and colon cancers. While some of the increased risk may be due to mutations in high-penetrance genes (i.e., BRCA2, PALB2, ATM, p16/CDKN2A or DNA mismatch repair genes), common genetic variants may also be involved. In a high-risk population of cases with either a family history of pancreatic cancer or early-onset pancreatic cancer (diagnosis before the age of 50 years), we examined the role of genetic variants previously associated with risk of pancreatic, breast, ovarian, or prostate cancer. We genotyped 985 cases (79 early-onset cases, 906 cases with a family history of pancreatic cancer) and 877 controls for 215,389 SNPs using the iSelect Collaborative Oncological Gene-Environment Study (iCOGS) array with custom content. Logistic regression was performed using a log-linear additive model. We replicated several previously reported pancreatic cancer susceptibility loci, including recently identified variants on 2p13.3 and 7p13 (2p13.3, rs1486134: OR = 1.36; 95% CI, 1.13–1.63; P = 9.29 × 10−4; 7p13, rs17688601: OR = 0.76; 95% CI, 0.63–0.93; P = 6.59 × 10−3). For the replicated loci, the magnitude of association observed in these high-risk patients was similar to that observed in studies of unselected patients. In addition to the established pancreatic cancer loci, we also found suggestive evidence of association (P &lt; 5 × 10−5) to pancreatic cancer for SNPs at HDAC9 (7p21.1) and COL6A2 (21q22.3). Even in high-risk populations, common variants influence pancreatic cancer susceptibility. Cancer Epidemiol Biomarkers Prev; 25(7); 1185–91. ©2016 AACR.

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 away). 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 64 normal derived pancreatic tissue samples, DIS3 expression showed the strongest association with a 30 bp indel variant in the risk locus (P = 4.8 × 10-4), 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, Sarah Manmiller, Jinping Jia, Hemang Parikh, 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 Translation of the Cancer Genome; Feb 7-9, 2015; San Francisco, CA. Philadelphia (PA): AACR; Cancer Res 2015;75(22 Suppl 1):Abstract nr A1-09.

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 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 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 &lt; 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.