Telomere deregulation is a hallmark of cancer. Telomere length measured in lymphocytes (LTL) has been shown to be a risk marker for several cancers. For pancreatic ductal adenocarcinoma (PDAC) consensus is lacking whether risk is associated with long or short telomeres. Mendelian randomization approaches have shown that a score built from SNPs associated with LTL could be used as a robust risk marker. We explored this approach in a large scale study within the PANcreatic Disease ReseArch (PANDoRA) consortium. We analyzed 10 SNPs (ZNF676-rs409627, TERT-rs2736100, CTC1-rs3027234, DHX35-rs6028466, PXK-rs6772228, NAF1-rs7675998, ZNF208-rs8105767, OBFC1-rs9420907, ACYP2-rs11125529 and TERC-rs10936599) alone and combined in a LTL genetic score ("teloscore", which explains 2.2% of the telomere variability) in relation to PDAC risk in 2,374 cases and 4,326 controls. We identified several associations with PDAC risk, among which the strongest were with the TERT-rs2736100 SNP (OR = 1.54; 95%CI 1.35-1.76; p = 1.54 × 10-10 ) and a novel one with the NAF1-rs7675998 SNP (OR = 0.80; 95%CI 0.73-0.88; p = 1.87 × 10-6 , ptrend = 3.27 × 10-7 ). The association of short LTL, measured by the teloscore, with PDAC risk reached genome-wide significance (p = 2.98 × 10-9 for highest vs. lowest quintile; p = 1.82 × 10-10 as a continuous variable). In conclusion, we present a novel genome-wide candidate SNP for PDAC risk (TERT-rs2736100), a completely new signal (NAF1-rs7675998) approaching genome-wide significance and we report a strong association between the teloscore and risk of pancreatic cancer, suggesting that telomeres are a potential risk factor for pancreatic cancer.
The CDKN2A (p16) gene plays a key role in pancreatic cancer etiology. It is one of the most commonly somatically mutated genes in pancreatic cancer, rare germline mutations have been found to be associated with increased risk of developing familiar pancreatic cancer and CDKN2A promoter hyper-methylation has been suggested to play a critical role both in pancreatic cancer onset and prognosis. In addition several unrelated SNPs in the 9p21.3 region, that includes the CDNK2A, CDNK2B and the CDNK2B-AS1 genes, are associated with the development of cancer in various organs. However, association between the common genetic variability in this region and pancreatic cancer risk is not clearly understood. We sought to fill this gap in a case-control study genotyping 13 single nucleotide polymorphisms (SNPs) in 2,857 pancreatic ductal adenocarcinoma (PDAC) patients and 6,111 controls in the context of the Pancreatic Disease Research (PANDoRA) consortium. We found that the A allele of the rs3217992 SNP was associated with an increased pancreatic cancer risk (ORhet=1.14, 95% CI 1.01-1.27, p=0.026, ORhom=1.30, 95% CI 1.12-1.51, p=0.00049). This pleiotropic variant is reported to be a mir-SNP that, by changing the binding site of one or more miRNAs, could influence the normal cell cycle progression and in turn increase PDAC risk. In conclusion, we observed a novel association in a pleiotropic region that has been found to be of key relevance in the susceptibility to various types of cancer and diabetes suggesting that the CDKN2A/B locus could represent a genetic link between diabetes and pancreatic cancer risk.
Genome-wide association studies (GWAS) have identified common pancreatic cancer susceptibility variants at 13 chromosomal loci in individuals of European descent. To identify new susceptibility variants, we performed imputation based on 1000 Genomes (1000G) Project data and association analysis using 5,107 case and 8,845 control subjects from 27 cohort and case-control studies that participated in the PanScan I-III GWAS. This analysis, in combination with a two-staged replication in an additional 6,076 case and 7,555 control subjects from the PANcreatic Disease ReseArch (PANDoRA) and Pancreatic Cancer Case-Control (PanC4) Consortia uncovered 3 new pancreatic cancer risk signals marked by single nucleotide polymorphisms (SNPs) rs2816938 at chromosome 1q32.1 (per allele odds ratio (OR) = 1.20, P = 4.88x10 -15), rs10094872 at 8q24.21 (OR = 1.15, P = 3.22x10 -9) and rs35226131 at 5p15.33 (OR = 0.71, P = 1.70x10 -8). These SNPs represent independent risk variants at previously identified pancreatic cancer risk loci on chr1q32.1 ( NR5A2), chr8q24.21 ( MYC) and chr5p15.33 ( CLPTM1L- TERT) as per analyses conditioned on previously reported susceptibility variants. We assessed expression of candidate genes at the three risk loci in histologically normal ( n = 10) and tumor ( n = 8) derived pancreatic tissue samples and observed a marked reduction of NR5A2 expression (chr1q32.1) in the tumors (fold change -7.6, P = 5.7x10 -8). This finding was validated in a second set of paired ( n = 20) histologically normal and tumor derived pancreatic tissue samples (average fold change for three NR5A2 isoforms -31.3 to -95.7, P = 7.5x10 -4-2.0x10 -3). Our study has identified new susceptibility variants independently conferring pancreatic cancer risk that merit functional follow-up to identify target genes and explain the underlying biology.