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 tumor-derived tissue samples (n=115) from The Cancer Genome Atlas (TCGA) was 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 tumor-derived tissues (FDR<0.1), with the strongest effects seen near transcriptional start sites (TSS). Approximately 23% and 42% of genes with significant cis-eQTLs appeared to be specific for tumor and normal derived tissues, respectively. Significant enrichment of cis-eQTL variants was noted in noncoding regulatory regions, in particular for pancreatic tissues (1.53–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 tumor-derived (P=8.3×10−5) tissues. The high linkage disequilibrium (LD) 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 dataset for further studies on gene expression and its regulation in pancreatic tissues.