Brown, Kevin M.

Abstract Background: Genome-wide association studies (GWAS) have linked single nucleotide polymorphisms (SNPs) to pancreatic ductal adenocarcinoma (PDAC) risk at over 20 genomic loci. The Pancreatic Cancer Cohort Consortium and Pancreatic Cancer Case-Control Consortium are currently expanding on previous GWAS studies for PDAC from ∼9,000 cases and ∼12,000 controls to ∼36,000 cases and ∼800,000 controls in individuals of European and Asian ancestry. This ongoing GWAS phase is estimated to detect up to 50 new risk signals. Post-GWAS studies aim to fine-map and functionally assess risk variants in disease-relevant cell types and contexts to identify target genes and underlying mechanisms of risk. Massively parallel reporter assays (MPRAs) are a powerful method of assessing allele-specific gene regulatory effects. Integrating sequences into the genome using lentiviruses (Lenti-MPRA) enables the capture of variant function in the context of the native cellular chromatin. An added benefit is that Lenti-MPRAs can be used for difficult to transfect cells such as primary cells, that often reflect the disease relevant cell types. Methods: We have used a lentiviral MPRA to test 228 fine-mapped SNPs at the chr5p15.33 multi-cancer risk locus in PDAC cell lines (MIA PaCa-2, PANC-1) and the normal-derived ductal pancreas (HPDE) cell line. Additionally, we have optimized the transduction of large MPRA libraries into primary pancreas cells. Tunicamycin was used to induce ER stress in cell lines and primary cells, assessing downstream gene expression using RT-qPCR. Results: Our pilot MPRA demonstrated high reproducibility in PDAC cell lines (r2=0.93-0.97, MIA PaCa-2, PANC-1 cells) and normal-derived pancreatic ductal cells (r2=0.68-0.86, HPDE cells). We identified 33 SNPs out of 228 with significant allele-specific transcriptional activity in PDAC cell lines and 5 in HPDE cells (FDR ≤0.05). An additional 6-15% of variants exhibited significant allele-specific transcriptional activity when in the frames of the SNP-centered tested oligos were shifted left or right. We demonstrate the inducibility of ER-stress with tunicamycin treatment over 24 hours in cells lines and primary pancreatic cells. Additionally, we observed increased expression of CLPTM1IL at the chr5p15.33 locus after tunicamycin treatment over 24 hours in cell lines and primary pancreatic cells, suggesting ER-stress and inflammation may influence gene expression at this locus. Conclusions: Lentiviral MPRA assays in commonly used PDAC cell lines have demonstrated good reproducibility. The optimization studies and preliminary work in primary pancreatic cells provide a solid foundation for the next steps of assessing all GWAS loci in ER stress, an important disease relevant condition. Citation Format: Minal B Patel, Aidan O'Brien, Jun Zhong, Irene Collins, Jason W Hoskins, Stephen J Chanock, Samuel O Antwi, Rachel Z Stolzenberg-Solomon, Alison P Klein, Brian M Wolpin, Sara Lindström, H. Efsun Arda, Kevin M Brown, Katelyn E Connelly, Laufey T Amundadottir, Meagan Jezek. The functional characterization of pancreatic ductal adenocarcinoma GWAS risk variants in primary pancreatic cells – A pilot study [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Advances in Pancreatic Cancer Research; 2024 Sep 15-18; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2024;84(17 Suppl_2):Abstract nr A018.

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.