Polymorphisms in transcription factor binding sites and enhancer regions and pancreatic ductal adenocarcinoma risk.

Genome-wide association studies (GWAS) are a powerful tool for detecting variants associated with complex traits and can help risk stratification and prevention strategies against pancreatic ductal adenocarcinoma (PDAC). However, the strict significance threshold commonly used makes it likely that many true risk loci are missed. Functional annotation of GWAS polymorphisms is a proven strategy to identify additional risk loci. We aimed to investigate single-nucleotide polymorphisms (SNP) in regulatory regions [transcription factor binding sites (TFBSs) and enhancers] that could change the expression profile of multiple genes they act upon and thereby modify PDAC risk. We analyzed a total of 12,636 PDAC cases and 43,443 controls from PanScan/PanC4 and the East Asian GWAS (discovery populations), and the PANDoRA consortium (replication population). We identified four associations that reached study-wide statistical significance in the overall meta-analysis: rs2472632(A) (enhancer variant, OR 1.10, 95%CI 1.06,1.13, p = 5.5 × 10-8), rs17358295(G) (enhancer variant, OR 1.16, 95%CI 1.10,1.22, p = 6.1 × 10-7), rs2232079(T) (TFBS variant, OR 0.88, 95%CI 0.83,0.93, p = 6.4 × 10-6) and rs10025845(A) (TFBS variant, OR 1.88, 95%CI 1.50,1.12, p = 1.32 × 10-5). The SNP with the most significant association, rs2472632, is located in an enhancer predicted to target the coiled-coil domain containing 34 oncogene. Our results provide new insights into genetic risk factors for PDAC by a focused analysis of polymorphisms in regulatory regions and demonstrating the usefulness of functional prioritization to identify loci associated with PDAC risk.

Potential association between PSCA rs2976395 functional variant and pancreatic cancer risk.

Correlated regions of systemic interindividual variation (CoRSIV) represent a small proportion of the human genome showing DNA methylation patterns that are the same in all human tissues, are different among individuals, and are partially regulated by genetic variants in cis. In this study we aimed at investigating single-nucleotide polymorphisms (SNPs) within CoRSIVs and their involvement with pancreatic ductal adenocarcinoma (PDAC) risk. We analyzed 29,099 CoRSIV-SNPs and 133,615 CoRSIV-mQTLs in 14,394 cases and 247,022 controls of European and Asian descent. We observed that the A allele of the rs2976395 SNP was associated with increased PDAC risk in Europeans (p = 2.81 × 10-5). This SNP lies in the prostate stem cell antigen gene and is in perfect linkage disequilibrium with a variant (rs2294008) that has been reported to be associated with risk of many other cancer types. The A allele is associated with the DNA methylation level of the gene according to the PanCan-meQTL database and with overexpression according to QTLbase. The expression of the gene has been observed to be deregulated in many tumors of the gastrointestinal tract including pancreatic cancer; however, functional studies are needed to elucidate the function relevance of the association.

Polymorphic variants involved in methylation regulation: a strategy to discover risk loci for pancreatic ductal adenocarcinoma.

INTRODUCTION: Only a small number of risk factors for pancreatic ductal adenocarcinoma (PDAC) has been established. Several studies identified a role of epigenetics and of deregulation of DNA methylation. DNA methylation is variable across a lifetime and in different tissues; nevertheless, its levels can be regulated by genetic variants like methylation quantitative trait loci (mQTLs), which can be used as a surrogate. MATERIALS AND METHODS: We scanned the whole genome for mQTLs and performed an association study in 14 705 PDAC cases and 246 921 controls. The methylation data were obtained from whole blood and pancreatic cancer tissue through online databases. We used the Pancreatic Cancer Cohort Consortium and the Pancreatic Cancer Case-Control Consortium genome-wide association study (GWAS) data as discovery phase and the Pancreatic Disease Research consortium, the FinnGen project and the Japan Pancreatic Cancer Research consortium GWAS as replication phase. RESULTS: The C allele of 15q26.1-rs12905855 showed an association with a decreased risk of PDAC (OR=0.90, 95% CI 0.87 to 0.94, p=4.93×10-8 in the overall meta-analysis), reaching genome-level statistical significance. 15q26.1-rs12905855 decreases the methylation of a 'C-phosphate-G' (CpG) site located in the promoter region of the RCCD1 antisense (RCCD1-AS1) gene which, when expressed, decreases the expression of the RCC1 domain-containing (RCCD1) gene (part of a histone demethylase complex). Thus, it is possible that the rs12905855 C-allele has a protective role in PDAC development through an increase of RCCD1 gene expression, made possible by the inactivity of RCCD1-AS1. CONCLUSION: We identified a novel PDAC risk locus which modulates cancer risk by controlling gene expression through DNA methylation.

Duodenal Fluid Analysis as a Rewarding Approach to Detect Low-Abundance Mutations in Biliopancreatic Cancers.

Diagnosis of biliopancreatic cancers by the available serum tumor markers, imaging, and histopathological tissue specimen examination remains a challenge. Circulating cell-free DNA derived from matched pairs of secretin-stimulated duodenal fluid (DF) and plasma from 10 patients with biliopancreatic diseases and 8 control subjects was analyzed using AmpliSeq™ HD technology for Ion Torrent Next-Generation Sequencing to evaluate the potential of liquid biopsy with DF in biliopancreatic cancers. The median cfDNA concentration was greater in DF-derived than in plasma-derived samples. A total of 13 variants were detected: 11 vs. 1 were exclusive for DF relative to the plasma source, and 1 was shared between the two body fluids. According to the four-tier systems, 10 clinical tier-I-II (76.9%), 1 tier-III (7.7%), and 2 tier-IV (15.4%) variants were identified. Notably, the 11 tier-I-III variants were exclusively found in DF-derived cfDNA from five patients with biliopancreatic cancers, and were detected in seven genes (KRAS, TP53, BRAF, CDKN2A, RNF43, GNAS, and PIK3CA); 82% of the tier-I-III variants had a low abundance, with a VAF < 6%. The mutational profiling of DF seems to be a reliable and promising tool for identifying cancer-associated alterations in malignant cancers of the biliopancreatic tract.

Focus on Achalasia in the Omics Era.

Achalasia is a rare and complex esophageal disease of unknown etiology characterized by difficulty in swallowing due to the lack of opening of the lower esophageal sphincter and the absence of esophageal peristalsis. Recent advancements in technology for analyzing DNA, RNA and biomolecules in high-throughput techniques are offering new opportunities to better understand the etiology and the pathogenetic mechanisms underlying achalasia. Through this narrative review of the scientific literature, we aim to provide a comprehensive assessment of the state-of-the-art knowledge on omics of achalasia, with particular attention to those considered relevant to the pathogenesis of the disease. The notion and importance of the multi-omics approach, its limitations and future directions are also introduced, and it is highlighted how the integration of single omics data will lead to new insights into the development of achalasia and offer clinical tools which will allow early diagnosis and better patient management.

A scan of all coding region variants of the human genome, identifies 13q12.2-rs9579139 and 15q24.1-rs2277598 as novel risk loci for pancreatic ductal adenocarcinoma.

Coding sequence variants comprise a small fraction of the germline genetic variability of the human genome. However, they often cause deleterious change in protein function and are therefore associated with pathogenic phenotypes. To identify novel pancreatic ductal adenocarcinoma (PDAC) risk loci, we carried out a complete scan of all common missense and synonymous SNPs and analysed them in a case-control study comprising four different populations, for a total of 14 538 PDAC cases and 190 657 controls. We observed a statistically significant association between 13q12.2-rs9581957-T and PDAC risk (P = 2.46 × 10-9), that is in linkage disequilibrium (LD) with a deleterious missense variant (rs9579139) of the URAD gene. Recent findings suggest that this gene is active in peroxisomes. Considering that peroxisomes have a key role as molecular scavengers, especially in eliminating reactive oxygen species, a malfunctioning URAD protein might expose the cell to a higher load of potentially DNA damaging molecules and therefore increase PDAC risk. The association was observed in individuals of European and Asian ethnicity. We also observed the association of the missense variant 15q24.1-rs2277598-T, that belongs to BBS4 gene, with increased PDAC risk (P = 1.53 × 10-6). rs2277598 is associated with body mass index and is in LD with diabetes susceptibility loci. In conclusion, we identified two missense variants associated with the risk of developing PDAC independently from the ethnicity highlighting the importance of conducting reanalysis of genome-wide association studies (GWASs) in light of functional data.

Association between a polymorphic variant in the CDKN2B-AS1/ANRIL gene and pancreatic cancer risk.

Genes carrying high-penetrance germline mutations may also be associated with cancer susceptibility through common low-penetrance genetic variants. To increase the knowledge on genetic pancreatic ductal adenocarcinoma (PDAC) aetiology, the common genetic variability of PDAC familial genes was analysed in our study. We conducted a multiphase study analysing 7745 single nucleotide polymorphisms (SNPs) from 29 genes reported to harbour a high-penetrance PDAC-associated mutation in at least one published study. To assess the effect of the SNPs on PDAC risk, a total of 14 666 PDAC cases and 221 897 controls across five different studies were analysed. The T allele of the rs1412832 polymorphism, that is situated in the CDKN2B-AS1/ANRIL, showed a genome-wide significant association with increased risk of developing PDAC (OR = 1.11, 95% CI = 1.07-1.15, P = 5.25 × 10-9 ). CDKN2B-AS1/ANRIL is a long noncoding RNA, situated in 9p21.3, and regulates many target genes, among which CDKN2A (p16) that frequently shows deleterious somatic and germline mutations and deregulation in PDAC. Our results strongly support the role of the genetic variability of the 9p21.3 region in PDAC aetiopathogenesis and highlight the importance of secondary analysis as a tool for discovering new risk loci in complex human diseases.

Evaluation of inherited germline mutations in cancer susceptibility genes among pancreatic cancer patients: a single-center study.

BACKGROUND: Germline mutations in cancer susceptibility genes were identified in pancreatic cancer (PanC) patients with a sporadic disease and in those unselected for family cancer history. METHODS: With the aim to determine the prevalence of germline predisposition genes mutations in PanC, and to evaluate whether they were associated with the presence of PanC, we profiled a custom AmpliSeq panel of 27 cancer susceptibility genes in 47 PanC patients and 51 control subjects by using the Ion Torrent PGM system. RESULTS: Multigene panel testing identified a total of 31 variants in 27 PanC (57.4%), including variants with pathogenic/likely pathogenic effect, those of uncertain significance, and variants whose clinical significance remains currently undefined. Five patients carried more than one variant in the same gene or in different genes. Eight patients (17.0%) had at least one pathogenic/likely pathogenic variant in four main genes: CFTR (10.6%), BRCA2 (8.5%), ATM and CHEK2 (2.1%). Pathogenic/likely pathogenic mutation were identified in patients with positive PanC family history (20%) or in patients without first-degree relatives affected by PanC (13.6%). All the BRCA2 mutation carriers were unselected PanC patients. The presence of mutations in BRCA2 was significantly associated with an increased occurrence of PanC and with positive family history for endometrial cancer (p = 0.018). CONCLUSIONS: This study confirmed the potential remarkable contribution of BRCA2 in assessing the presence of PanC. Overall our findings supported the recommendation of offering the germline testing to all the PanC patients with the intent to reduce the number of underdiagnosed carriers of mutations in predisposition genes, and not to preclude their relatives from the opportunity to benefit from surveillance programs.

Exploring the Neandertal legacy of pancreatic ductal adenocarcinoma risk in Eurasians.

BACKGROUND: The genomes of present-day non-Africans are composed of 1-3% of Neandertal-derived DNA as a consequence of admixture events between Neandertals and anatomically modern humans about 50-60 thousand years ago. Neandertal-introgressed single nucleotide polymorphisms (aSNPs) have been associated with modern human disease-related traits, which are risk factors for pancreatic ductal adenocarcinoma (PDAC), such as obesity, type 2 diabetes, and inflammation. In this study, we aimed at investigating the role of aSNPs in PDAC in three Eurasian populations. RESULTS: The high-coverage Vindija Neandertal genome was used to select aSNPs in non-African populations from 1000 Genomes project phase 3 data. Then, the association between aSNPs and PDAC risk was tested independently in Europeans and East Asians, using existing GWAS data on more than 200 000 individuals. We did not find any significant associations between aSNPs and PDAC in samples of European descent, whereas, in East Asians, we observed that the Chr10p12.1-rs117585753-T allele (MAF = 10%) increased the risk to develop PDAC (OR = 1.35, 95%CI 1.19-1.54, P = 3.59 × 10-6), with a P-value close to a threshold that takes into account multiple testing. CONCLUSIONS: Our results show only a minimal contribution of Neandertal SNPs to PDAC risk.

Genetically Determined Telomere Length Is Associated with Pancreatic Neuroendocrine Neoplasms Onset.

INTRODUCTION: Telomere length (TL) is a potential indicator of cancer predisposition; however, the multitude of techniques used to measure it causes the results to be heterogeneous and, in some cases, controversial. In the last years, several studies adopted a strategy based on TL-associated genetic variants to generate a polygenic score, often referred as teloscore, used in lieu of direct TL measurement. For pancreatic neuroendocrine neoplasms (PanNEN), this strategy has not been attempted yet. METHODS: A teloscore was generated using 11 SNPs (NAF1-rs7675998, ZNF676-rs409627, TERC-rs10936599, CTC1-rs3027234, PXK-rs6772228, DHX35-rs6028466, OBFC1-rs9420907, ZNF208-rs8105767, ACYP2-rs11125529, TERT-rs2736100, and ZBTB46-rs755017), and 291 PanNEN cases and 1,686 controls collected by the PANcreatic Disease ReseArch (PANDoRA) consortium were genotyped to analyse the association of the teloscore and its individual SNPs with the risk of developing PanNEN. RESULTS: An association between genetically determined long telomeres and the risk of developing PanNEN (OR = 1.99, CI: 1.33-2.98, p = 0.0008) for highest versus median (third) quintile was observed. In addition, two novel SNPs associated with PanNEN risk were identified: ZNF676-rs409627 (ORC/C_vs_G/G = 2.27, CI: 1.58-3.27, p = 8.80 × 10-6) and TERT-rs2736100 (ORC/A_vs_C/C = 2.03, CI: 1.42-2.91, p = 1.06 × 10-4). CONCLUSION: In conclusion, this study provides for the first time a clear indication of the association between long genetically determined telomeres and increased risk of developing PanNEN.

Polygenic and multifactorial scores for pancreatic ductal adenocarcinoma risk prediction.

BACKGROUND: Most cases of pancreatic ductal adenocarcinoma (PDAC) are asymptomatic in early stages, and the disease is typically diagnosed in advanced phases, resulting in very high mortality. Tools to identify individuals at high risk of developing PDAC would be useful to improve chances of early detection. OBJECTIVE: We generated a polygenic risk score (PRS) for PDAC risk prediction, combining the effect of known risk SNPs, and carried out an exploratory analysis of a multifactorial score. METHODS: We tested the associations of the individual known risk SNPs on up to 2851 PDAC cases and 4810 controls of European origin from the PANcreatic Disease ReseArch (PANDoRA) consortium. Thirty risk SNPs were included in a PRS, which was computed on the subset of subjects that had 100% call rate, consisting of 839 cases and 2040 controls in PANDoRA and 6420 cases and 4889 controls from the previously published Pancreatic Cancer Cohort Consortium I-III and Pancreatic Cancer Case-Control Consortium genome-wide association studies. Additional exploratory multifactorial scores were constructed by complementing the genetic score with smoking and diabetes. RESULTS: The scores were associated with increased PDAC risk and reached high statistical significance (OR=2.70, 95% CI 1.99 to 3.68, p=2.54×10-10 highest vs lowest quintile of the weighted PRS, and OR=14.37, 95% CI 5.57 to 37.09, p=3.64×10-8, highest vs lowest quintile of the weighted multifactorial score). CONCLUSION: We found a highly significant association between a PRS and PDAC risk, which explains more than individual SNPs and is a step forward in the direction of the construction of a tool for risk stratification in the population.

Downexpression of miR-200c-3p Contributes to Achalasia Disease by Targeting the PRKG1 Gene.

Achalasia is an esophageal smooth muscle motility disorder with unknown pathogenesis. Taking into account our previous results on the downexpression of miR-200c-3p in tissues of patients with achalasia correlated with an increased expression of PRKG1, SULF1, and SYDE1 genes, our aim was to explore the unknown biological interaction between these genes and human miR-200c-3p and if this relation could unravel their functional role in the etiology of achalasia. To search for putative miR-200c-3p binding sites in the 3'-UTR of PRKG1, SULF1 and SYDE1, a bioinformatics tool was used. To test whether PRKG1, SULF1, and SYDE1 are targeted by miR-200c-3p, a dual-luciferase reporter assay and quantitative PCR on HEK293 and fibroblast cell lines were performed. To explore the biological correlation between PRKG1 and miR-200c-3p, an immunoblot analysis was carried out. The overexpression of miR-200c-3p reduced the luciferase activity in cells transfected with a luciferase reporter containing a fragment of the 3'-UTR regions of PRKG1, SULF1, and SYDE1 which included the miR-200c-3p seed sequence. The deletion of the miR-200c-3p seed sequence from the 3'-UTR fragments abrogated this reduction. A negative correlation between miR-200c-3p and PRKG1, SULF1, and SYDE1 expression levels was observed. Finally, a reduction of the endogenous level of PRKG1 in cells overexpressing miR-200c-3p was detected. Our study provides, for the first time, functional evidence about the PRKG1 gene as a direct target and SULF1 and SYDE1 as potential indirect substrates of miR-200c-3p and suggests the involvement of NO/cGMP/PKG signaling in the pathogenesis of achalasia.

Associations between pancreatic expression quantitative traits and risk of pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is among the most lethal cancers. Its poor prognosis is predominantly due to the fact that most patients remain asymptomatic until the disease reaches an advanced stage, alongside the lack of early markers and screening strategies. A better understanding of PDAC risk factors is essential for the identification of groups at high risk in the population. Genome-wide association studies (GWAS) have been a powerful tool for detecting genetic variants associated with complex traits, including pancreatic cancer. By exploiting functional and GWAS data, we investigated the associations between polymorphisms affecting gene function in the pancreas (expression quantitative trait loci, eQTLs) and PDAC risk. In a two-phase approach, we analysed 13 713 PDAC cases and 43 784 controls and identified a genome-wide significant association between the A allele of the rs2035875 polymorphism and increased PDAC risk (P = 7.14 × 10-10). This allele is known to be associated with increased expression in the pancreas of the keratin genes KRT8 and KRT18, whose increased levels have been reported to correlate with various tumour cell characteristics. Additionally, the A allele of the rs789744 variant was associated with decreased risk of developing PDAC (P = 3.56 × 10-6). This single nucleotide polymorphism is situated in the SRGAP1 gene and the A allele is associated with higher expression of the gene, which in turn inactivates the cyclin-dependent protein 42 (CDC42) gene expression, thus decreasing the risk of PDAC. In conclusion, we present here a functional-based novel PDAC risk locus and an additional strong candidate supported by significant associations and plausible biological mechanisms.

Do myeloproliferative neoplasms and multiple myeloma share the same genetic susceptibility loci?

Myeloproliferative neoplasms (MPNs) are a group of diseases that cause myeloid hematopoietic cells to overproliferate. Epidemiological and familial studies suggest that genetic factors contribute to the risk of developing MPN, but the genetic susceptibility of MPN is still not well known. Indeed, only few loci are known to have a clear role in the predisposition to this disease. Some studies reported a diagnosis of MPNs and multiple myeloma (MM) in the same patients, but the biological causes are still unclear. We tested the hypothesis that the two diseases share at least partly the same genetic risk loci. In the context of a European multicenter study with 460 cases and 880 controls, we analyzed the effect of the known MM risk loci, individually and in a polygenic risk score (PRS). The most significant result was obtained among patients with chronic myeloid leukemia (CML) for PS0RS1C1-rs2285803, which showed to be associated with an increased risk (OR = 3.28, 95% CI 1.79-6.02, P = .00012, P = .00276 when taking into account multiple testing). Additionally, the PRS showed an association with MPN risk when comparing the last with the first quartile of the PRS (OR = 2.39, 95% CI 1.64-3.48, P = 5.98 × 10-6 ). In conclusion, our results suggest a potential common genetic background between MPN and MM, which needs to be further investigated.

Genome-wide scan of long noncoding RNA single nucleotide polymorphisms and pancreatic cancer susceptibility.

Pancreatic ductal adenocarcinoma (PDAC) is projected to become the second cancer-related cause of death by 2030. Identifying novel risk factors, including genetic risk loci, could be instrumental in risk stratification and implementation of prevention strategies. Long noncoding RNAs (lncRNAs) are involved in regulation of key biological processes, and the possible role of their genetic variability has been unexplored so far. Combining genome wide association studies and functional data, we investigated the genetic variability in all lncRNAs. We analyzed 9893 PDAC cases and 9969 controls and identified a genome-wide significant association between the rs7046076 SNP and risk of developing PDAC (P = 9.73 × 10-9 ). This SNP is located in the NONHSAG053086.2 (lnc-SMC2-1) gene and the risk allele is predicted to disrupt the binding of the lncRNA with the micro-RNA (miRNA) hsa-mir-1256 that regulates several genes involved in cell cycle, such as CDKN2B. The CDKN2B region is pleiotropic and its genetic variants have been associated with several human diseases, possibly though an imperfect interaction between lncRNA and miRNA. We present a novel PDAC risk locus, supported by a genome-wide statistical significance and a plausible biological mechanism.

Genome-wide association study identifies an early onset pancreatic cancer risk locus.

Early onset pancreatic cancer (EOPC) is a rare disease with a very high mortality rate. Almost nothing is known on the genetic susceptibility of EOPC, therefore, we performed a genome-wide association study (GWAS) to identify novel genetic variants specific for patients diagnosed with pancreatic ductal adenocarcinoma (PDAC) at younger ages. In the first phase, conducted on 821 cases with age of onset ≤60 years, of whom 198 with age of onset ≤50, and 3227 controls from PanScan I-II, we observed four SNPs (rs7155613, rs2328991, rs4891017 and rs12610094) showing an association with EOPC risk (P < 1 × 10-4 ). We replicated these SNPs in the PANcreatic Disease ReseArch (PANDoRA) consortium and used additional in silico data from PanScan III and PanC4. Among these four variants rs2328991 was significant in an independent set of 855 cases with age of onset ≤60 years, of whom 265 with age of onset ≤50, and 4142 controls from the PANDoRA consortium while in the in silico data, we observed no statistically significant association. However, the resulting meta-analysis supported the association (P = 1.15 × 10-4 ). In conclusion, we propose a novel variant rs2328991 to be involved in EOPC risk. Even though it was not possible to find a mechanistic link between the variant and the function, the association is supported by a solid statistical significance obtained in the largest study on EOPC genetics present so far in the literature.

Genetic variability of the ABCC2 gene and clinical outcomes in pancreatic cancer patients.

Pancreatic ductal adenocarcinoma (PDAC) has an extremely poor prognosis, caused by various factors, such as the aggressiveness of the disease, the limited therapeutic options and the lack of early detection and risk markers. The ATP binding cassette subfamily C member 2 (ABCC2) protein plays a critical role in response to various drugs and is differentially expressed in gemcitabine sensitive and resistant cells. Moreover, single nucleotide polymorphisms (SNPs) in the gene have been associated with differential outcomes and prognosis in several tumour types. The aim of this study was to investigate the possible association between SNPs in the ABCC2 gene and overall survival (OS) in PDAC patients. We analysed 12 polymorphisms, including tagging-SNPs covering all the genetic variability of the ABCC2 gene and genotyped them in 1415 PDAC patients collected within the Pancreatic Disease ReseArch (PANDoRA) consortium. We tested the association between ABCC2 SNPs and PDAC OS using Cox proportional hazard models. We analysed PDAC patients dividing them by stage and observed that the minor alleles of three SNPs showed an association with worse OS [rs3740067: hazard ratio (HR) = 3.29, 95% confidence interval (CI) = 1.56-6.97, P = 0.002; rs3740073: HR = 3.11, 95% CI = 1.52-6.38, P = 0.002 and rs717620: HR = 2.90, 95% CI = 1.41-5.95, P = 0.004, respectively] in stage I patients. In patients with more advanced PDAC, we did not observe any statistically significant association. Our results suggest that rs3740067, rs3740073 and rs717620 could be promising prognostic markers in stage I PDAC patients.

Genetic determinants of telomere length and risk of pancreatic cancer: A PANDoRA study.

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.

High methylation levels of PCDH10 predict poor prognosis in patients with pancreatic ductal adenocarcinoma.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal malignancies and is not a clinically homogeneous disease, but subsets of patients with distinct prognosis and response to therapy can be identified by genome-wide analyses. Mutations in major PDAC driver genes were associated with poor survival. By bioinformatics analysis, we identified protocadherins among the most frequently mutated genes in PDAC suggesting an important role of these genes in the biology of this tumor. Promoter methylation of protocadherins has been suggested as a prognostic marker in different tumors, but in PDAC this epigenetic modification has not been extensively studied. Thus, we evaluated whether promoter methylation of three frequently mutated protocadherins, PCDHAC2, PCDHGC5 and PCDH10 could be used as survival predictors in PDAC patients. METHODS: DNA extracted from 23 PDACs and adjacent non-neoplastic pancreatic tissues were bisulfite treated. Combined Bisulfite Restriction Analysis (COBRA) coupled to denaturing high-performance liquid chromatography (dHPLC) detection and bisulfite genomic sequencing (BGS) were used to determine the presence of methylated CpG dinucleotides in the promoter amplicons analyzed. RESULTS: In an exploratory analysis, two protocadherins showed the same pattern of CpG methylation in PDAC and adjacent non-neoplastic pancreatic tissues: lack of methylation for PCDHAC2, complete methylation for PCDHGC5. Conversely, the third protocadherin analyzed, PCDH10, showed a variable degree of CpG methylation in PDAC and absence of methylation in adjacent non-neoplastic pancreatic tissues. At Kaplan-Meier analysis, high levels of PCDH10 methylation defined according to the receiver operating characteristic (ROC) curve analysis were significantly associated with worse progression-free survival (PFS) rates (P = 0.008), but not with overall survival (OS). High levels of PCDH10 methylation were a prognostic factor influencing PFS (HR = 4.0: 95% CI, 1.3-12.3; P = 0.016), but not the OS. CONCLUSIONS: In this study, we show for the first time that the methylation status of PCDH10 can predict prognosis in PDAC patients with a significant impact on the outcome in terms of progression-free survival. High levels of PCDH10 promoter methylation could be useful to identify patients at high risk of disease progression, contributing to a more accurate stratification of PDAC patients for personalized clinical management.

Common genetic variants associated with pancreatic adenocarcinoma may also modify risk of pancreatic neuroendocrine neoplasms.

Pancreatic neuroendocrine neoplasms (pNEN) account for less than 5% of all pancreatic neoplasms and genetic association studies on susceptibility to the disease are limited. We sought to identify possible overlap of genetic susceptibility loci between pancreatic ductal adenocarcinoma (PDAC) and pNEN; therefore, PDAC susceptibility variants (n = 23) from Caucasian genome-wide association studies (GWAS) were genotyped in 369 pNEN cases and 3277 controls from the PANcreatic Disease ReseArch (PANDoRA) consortium to evaluate the odds associated with pNEN risk, disease onset and tumor characteristics. Main effect analyses showed four PDAC susceptibility variants-rs9854771, rs1561927, rs9543325 and rs10919791 to be associated with pNEN risk. Subsequently, only associations with rs9543325, rs10919791 and rs1561927 were noteworthy with false positive report probability (FPRP) tests. Stratified analyses considering age at onset (50-year threshold), showed rs2736098, rs16986825 and rs9854771 to be associated with risk of developing pNEN at a younger age. Stratified analyses also showed some single nucleotide polymorphisms to be associated with different degrees of tumor grade, metastatic potential and functionality. Our results identify known GWAS PDAC susceptibility loci, which may also be involved in sporadic pNEN etiology and suggest that some genetic mechanisms governing pathogenesis of these two entities may be similar, with few of these loci being more influential in younger cases or tumor subtypes.