Search

The prevalence of childhood obesity is increasing worldwide, along with the associated common comorbidities of type 2 diabetes and cardiovascular disease in later life. Motivated by evidence for a strong genetic component, our prior genome-wide association study (GWAS) efforts for childhood obesity revealed 19 independent signals for the trait; however, the mechanism of action of these loci remains to be elucidated. To molecularly characterize these childhood obesity loci, we sought to determine the underlying causal variants and the corresponding effector genes within diverse cellular contexts. Integrating childhood obesity GWAS summary statistics with our existing 3D genomic datasets for 57 human cell types, consisting of high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq, and RNA-seq, we applied stratified LD score regression and calculated the proportion of genome-wide SNP heritability attributable to cell type-specific features, revealing pancreatic alpha cell enrichment as the most statistically significant. Subsequent chromatin contact-based fine-mapping was carried out for genome-wide significant childhood obesity loci and their linkage disequilibrium proxies to implicate effector genes, yielded the most abundant number of candidate variants and target genes at the BDNF , ADCY3 , TMEM18, and FTO loci in skeletal muscle myotubes and the pancreatic beta-cell line, EndoC-BH1. One novel implicated effector gene, ALKAL2 - an inflammation-responsive gene in nerve nociceptors - was observed at the key TMEM18 locus across multiple immune cell types. Interestingly, this observation was also supported through colocalization analysis using expression quantitative trait loci (eQTL) derived from the Genotype-Tissue Expression (GTEx) dataset, supporting an inflammatory and neurologic component to the pathogenesis of childhood obesity. Our comprehensive appraisal of 3D genomic datasets generated in a myriad of different cell types provides genomic insights into pediatric obesity pathogenesis., Competing Interests: KT, MP, JP, CS, SL, PS, NK, LG, NT, JO, YW, KH, AJ, JH, LA, MX, KB, SA, WY, PT, PS, KK, LC, ML, BZ, AC, AW, SG No competing interests declared

Identification of somatic driver mutations in the noncoding genome remains challenging. To comprehensively characterize noncoding driver mutations for pancreatic ductal adenocarcinoma (PDAC), we first created genome-scale maps of accessible chromatin regions (ACRs) and histone modification marks (HMMs) in pancreatic cell lines and purified pancreatic acinar and duct cells. Integration with whole-genome mutation calls from 506 PDACs revealed 314 ACRs/HMMs significantly enriched with 3,614 noncoding somatic mutations (NCSMs). Functional assessment using massively parallel reporter assays (MPRA) identified 178 NCSMs impacting reporter activity (19.45% of those tested). Focused luciferase validation confirmed negative effects on gene regulatory activity for NCSMs near CDKN2A and ZFP36L2 . For the latter, CRISPR interference (CRISPRi) further identified ZFP36L2 as a target gene (16.0 - 24.0% reduced expression, P = 0.023-0.0047) with disrupted KLF9 binding likely mediating the effect. Our integrative approach provides a catalog of potentially functional noncoding driver mutations and nominates ZFP36L2 as a PDAC driver gene., Competing Interests: Competing Interests statement The authors declare no competing interests.

Pancreatic Ductal Adenocarcinoma (PDAC) is the third leading cause of cancer-related deaths in the U.S. Both rare and common germline variants contribute to PDAC risk. Here, we fine-map and functionally characterize a common PDAC risk signal at 1p36.33 (tagged by rs13303010) identified through a genome wide association study (GWAS). One of the fine-mapped SNPs, rs13303160 (r 2 =0.93 in 1000G EUR samples, OR=1.23, P value=2.74x10 -9 ) demonstrated allele-preferential gene regulatory activity in vitro and allele-preferential binding of JunB and JunD in vitro and in vivo . Expression Quantitative Trait Locus (eQTL) analysis identified KLHL17 as a likely target gene underlying the signal. Proteomic analysis identified KLHL17 as a member of the Cullin-E3 ubiquitin ligase complex in PDAC-derived cells. In silico differential gene expression analysis of the GTExv8 pancreas data suggested an association between lower KLHL17 (risk associated) and pro-inflammatory pathways. We hypothesize that KLHL17 may mitigate inflammation by recruiting pro-inflammatory proteins for ubiquitination and degradation thereby influencing PDAC risk., Competing Interests: Competing Interests statement The authors declare no competing interests.

Genetic and epigenetic variations in regulatory enhancer elements increase susceptibility to a range of pathologies. Despite recent advances, linking enhancer elements to target genes and predicting transcriptional outcomes of enhancer dysfunction remain significant challenges. Using 3D chromatin conformation assays, we generated an extensive enhancer interaction dataset for the human pancreas, encompassing more than 20 donors and five major cell types, including both exocrine and endocrine compartments. We employed a network approach to parse chromatin interactions into enhancer-promoter tree models, facilitating a quantitative, genome-wide analysis of enhancer connectivity. With these tree models, we developed a machine learning algorithm to estimate the impact of enhancer perturbations on cell type-specific gene expression in the human pancreas. Orthogonal to our computational approach, we perturbed enhancer function in primary human pancreas cells using CRISPR interference and quantified the effects at the single-cell level through RNA FISH coupled with high-throughput imaging. Our enhancer tree models enabled the annotation of common germline risk variants associated with pancreas diseases, linking them to putative target genes in specific cell types. For pancreatic ductal adenocarcinoma, we found a stronger enrichment of disease susceptibility variants within acinar cell regulatory elements, despite ductal cells historically being assumed as the primary cell-of-origin. Our integrative approach-combining cell type-specific enhancer-promoter interaction mapping, computational models, and single-cell enhancer perturbation assays-produced a robust resource for studying the genetic basis of pancreas disorders.

Background: Pancreatic cancer is a leading cause of cancer-related death globally. Risk factors for pancreatic cancer include common genetic variants and potentially heavy alcohol consumption. We assessed if genetic variants modify the association between heavy alcohol consumption and pancreatic cancer risk., Methods: We conducted a genome-wide interaction analysis of single-nucleotide polymorphisms (SNP) by heavy alcohol consumption (more than three drinks per day) for pancreatic cancer in European ancestry populations from genome-wide association studies. Our analysis included 3,707 cases and 4,167 controls from case-control studies and 1,098 cases and 1,162 controls from cohort studies. Fixed-effect meta-analyses were conducted., Results: A potential novel region of association on 10p11.22, lead SNP rs7898449 (interaction P value (Pinteraction) = 5.1 × 10-8 in the meta-analysis; Pinteraction = 2.1 × 10-9 in the case-control studies; Pinteraction = 0.91 in the cohort studies), was identified. An SNP correlated with this lead SNP is an expression quantitative trait locus for the neuropilin 1 gene. Of the 17 genomic regions with genome-wide significant evidence of association with pancreatic cancer in prior studies, we observed suggestive evidence that heavy alcohol consumption modified the association for one SNP near LINC00673, rs11655237 on 17q25.1 (Pinteraction = 0.004)., Conclusions: We identified a novel genomic region that may be associated with pancreatic cancer risk in conjunction with heavy alcohol consumption located near an expression quantitative trait locus for neuropilin 1, a protein that plays an important role in the development and progression of pancreatic cancer., Impact: This work can provide insights into the etiology of pancreatic cancer, particularly in heavy drinkers., (©2024 American Association for Cancer Research.)

The prevalence of childhood obesity is increasing worldwide, along with the associated common comorbidities of type 2 diabetes and cardiovascular disease in later life. Motivated by evidence for a strong genetic component, our prior genome-wide association study (GWAS) efforts for childhood obesity revealed 19 independent signals for the trait; however, the mechanism of action of these loci remains to be elucidated. To molecularly characterize these childhood obesity loci we sought to determine the underlying causal variants and the corresponding effector genes within diverse cellular contexts. Integrating childhood obesity GWAS summary statistics with our existing 3D genomic datasets for 57 human cell types, consisting of high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq, and RNA-seq, we applied stratified LD score regression and calculated the proportion of genome-wide SNP heritability attributable to cell type-specific features, revealing pancreatic alpha cell enrichment as the most statistically significant. Subsequent chromatin contact-based fine-mapping was carried out for genome-wide significant childhood obesity loci and their linkage disequilibrium proxies to implicate effector genes, yielded the most abundant number of candidate variants and target genes at the BDNF , ADCY3 , TMEM18 and FTO loci in skeletal muscle myotubes and the pancreatic beta-cell line, EndoC-BH1. One novel implicated effector gene, ALKAL2 - an inflammation-responsive gene in nerve nociceptors - was observed at the key TMEM18 locus across multiple immune cell types. Interestingly, this observation was also supported through colocalization analysis using expression quantitative trait loci (eQTL) derived from the Genotype-Tissue Expression (GTEx) dataset, supporting an inflammatory and neurologic component to the pathogenesis of childhood obesity. Our comprehensive appraisal of 3D genomic datasets generated in a myriad of different cell types provides genomic insights into pediatric obesity pathogenesis.

A portion of the genetic basis for many common autoimmune disorders has been uncovered by genome-wide association studies (GWAS), but GWAS do not reveal causal variants, effector genes, or the cell types impacted by disease-associated variation. We have generated 3D genomic datasets consisting of promoter-focused Capture-C, Hi-C, ATAC-seq, and RNA-seq and integrated these data with GWAS of 16 autoimmune traits to physically map disease-associated variants to the effector genes they likely regulate in 57 human cell types. These 3D maps of gene cis -regulatory architecture are highly powered to identify the cell types most likely impacted by disease-associated genetic variation compared to 1D genomic features, and tend to implicate different effector genes than eQTL approaches in the same cell types. Most of the variants implicated by these cis -regulatory architectures are highly trait-specific, but nearly half of the target genes connected to these variants are shared across multiple autoimmune disorders in multiple cell types, suggesting a high level of genetic diversity and complexity among autoimmune diseases that nonetheless converge at the level of target gene and cell type. Substantial effector gene sharing led to the common enrichment of similar biological networks across disease and cell types. However, trait-specific pathways representing potential areas for disease-specific intervention were identified. To test this, we pharmacologically validated squalene synthase, a cholesterol biosynthetic enzyme encoded by the FDFT1 gene implicated by our approach in MS and SLE, as a novel immunomodulatory drug target controlling inflammatory cytokine production by human T cells. These data represent a comprehensive resource for basic discovery of gene cis -regulatory mechanisms, and the analyses reported reveal mechanisms by which autoimmune-associated variants act to regulate gene expression, function, and pathology across multiple, distinct tissues and cell types.

Objective: Genome wide association studies have identified an exon 6 CTRB2 deletion variant that associates with increased risk of pancreatic cancer. To acquire evidence on its causal role, we developed a new mouse strain carrying an equivalent variant in Ctrb1 , the mouse orthologue of CTRB2 ., Design: We used CRISPR/Cas9 to introduce a 707bp deletion in Ctrb1 encompassing exon 6 ( Ctrb1 Δexon6 ). This mutation closely mimics the human deletion variant. Mice carrying the mutant allele were extensively profiled at 3 months to assess their phenotype., Results: Ctrb1 Δexon6 mutant mice express a truncated CTRB1 that accumulates in the ER. The pancreas of homozygous mutant mice displays reduced chymotrypsin activity and total protein synthesis. The histological aspect of the pancreas is inconspicuous but ultrastructural analysis shows evidence of dramatic ER stress and cytoplasmic and nuclear inclusions. Transcriptomic analyses of the pancreas of mutant mice reveals acinar program down-regulation and increased activity of ER stress-related and inflammatory pathways. Heterozygous mice have an intermediate phenotype. Agr2 is one of the most up-regulated genes in mutant pancreata. Ctrb1 Δexon6 mice exhibit impaired recovery from acute caerulein-induced pancreatitis. Administration of TUDCA or sulindac partially alleviates the phenotype. A transcriptomic signature derived from the mutant pancreata is significantly enriched in normal human pancreas of CTRB2 exon 6 deletion variant carriers from the GTEx cohort., Conclusions: This mouse strain provides formal evidence that the Ctrb1 Δexon6 variant causes ER stress and inflammation in vivo , providing an excellent model to understand its contribution to pancreatic ductal adenocarcinoma development and to identify preventive strategies., Summary Box: What is already known about this subject?: - CTRB2 is one of the most abundant proteins produced by human pancreatic acinar cells. - A common exon 6 deletion variant in CTRB2 has been associated with an increased risk of pancreatic ductal adenocarcinoma. - Misfolding of digestive enzymes is associated with pancreatic pathology. What are the new findings?: - We developed a novel genetic model that recapitulates the human CTRB2 deletion variant in the mouse orthologue, Ctrb1 . - Truncated CTRB1 misfolds and accumulates in the ER; yet, mutant mice display a histologically normal pancreas at 3 months age.- CTRB1 and associated chaperones colocalize in the ER, the cytoplasm, and the nucleus of acinar cells.- Transcriptomics analysis reveals reduced activity of the acinar program and increased activity of pathways involved in ER stress, unfolded protein response, and inflammation.- Mutant mice are sensitized to pancreatic damage and do not recover properly from a mild caerulein-induced pancreatitis.- TUDCA administration partially relieves the ER stress in mutant mice. How might it impact on clinical practice in the foreseeable future?: - The new mouse model provides a tool to identify the mechanisms leading to increased pancreatic cancer risk in CTRB2 exon 6 carriers. - The findings suggest that drugs that cause ER stress relief and/or reduce inflammation might provide preventive opportunities.

Recent genome-wide association studies (GWAS) have revealed shared genetic components among alcohol, opioid, tobacco and cannabis use disorders. However, the extent of the underlying shared causal variants and effector genes, along with their cellular context, remain unclear. We leveraged our existing 3D genomic datasets comprising high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq and RNA-seq across >50 diverse human cell types to focus on genomic regions that coincide with GWAS loci. Using stratified LD regression, we determined the proportion of genomewide SNP heritability attributable to the features assayed across our cell types by integrating recent GWAS summary statistics for the relevant traits: alcohol use disorder (AUD), tobacco use disorder (TUD), opioid use disorder (OUD) and cannabis use disorder (CanUD). Statistically significant enrichments ( P <0.05) were observed in 14 specific cell types, with heritability reaching 9.2-fold for iPSC-derived cortical neurons and neural progenitors, confirming that they are crucial cell types for further functional exploration. Additionally, several pancreatic cell types, notably pancreatic beta cells, showed enrichment for TUD, with heritability enrichments up to 4.8-fold, suggesting genomic overlap with metabolic processes. Further investigation revealed significant positive genetic correlations between T2D with both TUD and CanUD (FDR<0.05) and a significant negative genetic correlation with AUD. Interestingly, after partitioning the heritability for each cell type's cis-regulatory elements, the correlation between T2D and TUD for pancreatic beta cells was greater (r=0.2) than the global genetic correlation value. Our study provides new genomic insights into substance use disorders and implicates cell types where functional follow-up studies could reveal causal variant-gene mechanisms underpinning these disorders., Competing Interests: COMPETING INTERESTS Dr. Kranzler is a member of advisory boards for Dicerna Pharmaceuticals, Sophrosyne Pharmaceuticals, Enthion Pharmaceuticals, and Clearmind Medicine; a consultant to Sobrera Pharmaceuticals and Altimmune; the recipient of research funding and medication supplies for an investigator-initiated study from Alkermes; a member of the American Society of Clinical Psychopharmacology’s Alcohol Clinical Trials Initiative, which was supported in the last three years by Alkermes, Dicerna, Ethypharm, Lundbeck, Mitsubishi, Otsuka, and Pear Therapeutics; and a holder of U.S. patent 10,900,082 titled: “Genotype-guided dosing of opioid agonists,” issued 26 January 2021. The remaining authors have nothing to disclose.

Genome-wide association studies (GWASs) have identified numerous lung cancer risk-associated loci. However, decoding molecular mechanisms of these associations is challenging since most of these genetic variants are non-protein-coding with unknown function. Here, we implemented massively parallel reporter assays (MPRAs) to simultaneously measure the allelic transcriptional activity of risk-associated variants. We tested 2,245 variants at 42 loci from 3 recent GWASs in East Asian and European populations in the context of two major lung cancer histological types and exposure to benzo(a)pyrene. This MPRA approach identified one or more variants (median 11 variants) with significant effects on transcriptional activity at 88% of GWAS loci. Multimodal integration of lung-specific epigenomic data demonstrated that 63% of the loci harbored multiple potentially functional variants in linkage disequilibrium. While 22% of the significant variants showed allelic effects in both A549 (adenocarcinoma) and H520 (squamous cell carcinoma) cell lines, a subset of the functional variants displayed a significant cell-type interaction. Transcription factor analyses nominated potential regulators of the functional variants, including those with cell-type-specific expression and those predicted to bind multiple potentially functional variants across the GWAS loci. Linking functional variants to target genes based on four complementary approaches identified candidate susceptibility genes, including those affecting lung cancer cell growth. CRISPR interference of the top functional variant at 20q13.33 validated variant-to-gene connections, including RTEL1, SOX18, and ARFRP1. Our data provide a comprehensive functional analysis of lung cancer GWAS loci and help elucidate the molecular basis of heterogeneity and polygenicity underlying lung cancer susceptibility., Competing Interests: Declaration of interests The authors declare no competing interests., (Published by Elsevier Inc.)

The complex relationship between chromatin accessibility, transcriptional regulation, and cancer transitions presents a daunting puzzle. Terekhanova et al. created a pan-cancer epigenetic and transcriptomic atlas at single-cell resolution, yielding important insights into the underlying chromatin architecture of cancer transitions and novel discoveries with the potential to advance precision medicine., Competing Interests: Declaration of interests The authors declare no competing interests., (Published by Elsevier Ltd.)

Here, we present a protocol to identify transcriptional regulators potentially mediating downstream biological effects of germline variants associated with complex traits of interest, which enables functional hypothesis generation independent of colocalizing expression quantitative trait loci (eQTLs). We describe steps for tissue-/cell-type-specific co-expression network modeling, expression regulator activity inference, and identification of representative phenotypic master regulators. Finally, we detail activity QTL and eQTL analyses. This protocol requires genotype, expression, and relevant covariables and phenotype data from existing eQTL datasets. For complete details on the use and execution of this protocol, please refer to Hoskins et al. 1 ., Competing Interests: Declaration of interests The authors declare no competing interests., (Published by Elsevier Inc.)

Background: In Western populations, pancreatic ductal adenocarcinoma (PDAC) risk has been found to be greater among individuals with non-O blood types than those with O blood type. However, the association has not been fully evaluated with respect to FUT2 (determining secretor status) and FUT3 (determining Lewis antigens) status, two biologically important genes in the expression of ABO blood groups with PDAC., Methods: We examined interactions in data from 8,027 cases and 11,362 controls in large pancreatic cancer consortia (PanScan I-III and PanC4) by using genetic variants to predict ABO blood groups (rs505922 and rs8176746), secretor status (rs601338), and Lewis antigens (rs812936, rs28362459, and rs3894326). Multivariable logistic regression was used to estimate ORs and 95% confidence intervals (CI) of the risk of PDAC adjusted for age and sex. We examined multiplicative interactions of ABO with secretor status and Lewis antigens by considering each product term between ABO and secretor and between ABO and Lewis antigens individually., Results: We found that the increased risk associated with non-O blood groups was somewhat stronger among secretors than nonsecretors [ORs, 1.28 (95% CI, 1.15-1.42) and 1.17 (95% CI, 1.03-1.32) respectively; Pinteraction = 0.002]. We did not find any interactions between ABO and Lewis antigens., Conclusions: Our large consortia data provide evidence of effect modification in the association between non-O blood type and pancreatic cancer risk by secretor status., Impact: Our results indicate that the association between ABO blood type and PDAC risk may vary by secretor status, but not by Lewis antigens., (©2023 American Association for Cancer Research.)

Background: There are conflicting data on whether nonalcoholic fatty liver disease (NAFLD) is associated with susceptibility to pancreatic cancer. Using Mendelian randomization (MR), we investigated the relationship between genetic predisposition to NAFLD and risk for pancreatic cancer., Methods: Data from genome-wide association studies (GWAS) within the Pancreatic Cancer Cohort Consortium (PanScan; cases n = 5,090, controls n = 8,733) and the Pancreatic Cancer Case Control Consortium (PanC4; cases n = 4,163, controls n = 3,792) were analyzed. We used data on 68 genetic variants with four different MR methods [inverse variance weighting (IVW), MR-Egger, simple median, and penalized weighted median] separately to predict genetic heritability of NAFLD. We then assessed the relationship between each of the four MR methods and pancreatic cancer risk, using logistic regression to calculate ORs and 95% confidence intervals (CI), adjusting for PC risk factors, including obesity and diabetes., Results: No association was found between genetically predicted NAFLD and pancreatic cancer risk in the PanScan or PanC4 samples [e.g., PanScan, IVW OR, 1.04; 95% confidence interval (CI), 0.88-1.22; MR-Egger OR, 0.89; 95% CI, 0.65-1.21; PanC4, IVW OR, 1.07; 95% CI, 0.90-1.27; MR-Egger OR, 0.93; 95% CI, 0.67-1.28]. None of the four MR methods indicated an association between genetically predicted NAFLD and pancreatic cancer risk in either sample., Conclusions: Genetic predisposition to NAFLD is not associated with pancreatic cancer risk., Impact: Given the close relationship between NAFLD and metabolic conditions, it is plausible that any association between NAFLD and pancreatic cancer might reflect host metabolic perturbations (e.g., obesity, diabetes, or metabolic syndrome) and does not necessarily reflect a causal relationship between NAFLD and pancreatic cancer., (©2023 American Association for Cancer Research.)

Background: The shared inherited genetic contribution to risk of different cancers is not fully known. In this study, we leverage results from 12 cancer genome-wide association studies (GWAS) to quantify pairwise genome-wide genetic correlations across cancers and identify novel cancer susceptibility loci., Methods: We collected GWAS summary statistics for 12 solid cancers based on 376 759 participants with cancer and 532 864 participants without cancer of European ancestry. The included cancer types were breast, colorectal, endometrial, esophageal, glioma, head and neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancers. We conducted cross-cancer GWAS and transcriptome-wide association studies to discover novel cancer susceptibility loci. Finally, we assessed the extent of variant-specific pleiotropy among cancers at known and newly identified cancer susceptibility loci., Results: We observed widespread but modest genome-wide genetic correlations across cancers. In cross-cancer GWAS and transcriptome-wide association studies, we identified 15 novel cancer susceptibility loci. Additionally, we identified multiple variants at 77 distinct loci with strong evidence of being associated with at least 2 cancer types by testing for pleiotropy at known cancer susceptibility loci., Conclusions: Overall, these results suggest that some genetic risk variants are shared among cancers, though much of cancer heritability is cancer-specific and thus tissue-specific. The increase in statistical power associated with larger sample sizes in cross-disease analysis allows for the identification of novel susceptibility regions. Future studies incorporating data on multiple cancer types are likely to identify additional regions associated with the risk of multiple cancer types., (© The Author(s) 2023. Published by Oxford University Press. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Competing Interests: The authors declare no competing interests.

Background: Pancreatic cancer presents as advanced disease in >80% of patients; yet, appropriate ages to consider prevention and early detection strategies are poorly defined. We investigated age-specific associations and attributable risks of pancreatic cancer for established modifiable and non-modifiable risk factors., Patients and Methods: We included 167 483 participants from two prospective US cohort studies with 1190 incident cases of pancreatic cancer during >30 years of follow-up; 5107 pancreatic cancer cases and 8845 control participants of European ancestry from a completed multicenter genome-wide association study (GWAS); and 248 893 pancreatic cancer cases documented in the US Surveillance, Epidemiology, and End Results (SEER) Program. Across different age categories, we investigated cigarette smoking, obesity, diabetes, height, and non-O blood group in the prospective cohorts; weighted polygenic risk score of 22 previously identified single nucleotide polymorphisms in the GWAS; and male sex and black race in the SEER Program., Results: In the prospective cohorts, all five risk factors were more strongly associated with pancreatic cancer risk among younger participants, with associations attenuated among those aged >70 years. The hazard ratios comparing participants with three to five risk factors with those with no risk factors were 9.24 [95% confidence interval (CI) 4.11-20.77] among those aged ≤60 years, 3.00 (95% CI 1.85-4.86) among those aged 61-70 years, and 1.46 (95% CI 1.10-1.94) among those aged >70 years (P heterogeneity  = 3×10 -5 ). These factors together were related to 65.6%, 49.7%, and 17.2% of incident pancreatic cancers in these age groups, respectively. In the GWAS and the SEER Program, the associations with the polygenic risk score, male sex, and black race were all stronger among younger individuals (P heterogeneity ≤0.01)., Conclusions: Established risk factors are more strongly associated with earlier-onset pancreatic cancer, emphasizing the importance of age at initiation for cancer prevention and control programs targeting this highly lethal malignancy., (Copyright © 2022 European Society for Medical Oncology. All rights reserved.)

Expression QTL (eQTL) analyses have suggested many genes mediating genome-wide association study (GWAS) signals but most GWAS signals still lack compelling explanatory genes. We have leveraged an adipose-specific gene regulatory network to infer expression regulator activities and phenotypic master regulators (MRs), which were used to detect activity QTLs (aQTLs) at cardiometabolic trait GWAS loci. Regulator activities were inferred with the VIPER algorithm that integrates enrichment of expected expression changes among a regulator's target genes with confidence in their regulator-target network interactions and target overlap between different regulators (i.e., pleiotropy). Phenotypic MRs were identified as those regulators whose activities were most important in predicting their respective phenotypes using random forest modeling. While eQTLs were typically more significant than aQTLs in cis, the opposite was true among candidate MRs in trans. Several GWAS loci colocalized with MR trans-eQTLs/aQTLs in the absence of colocalized cis-QTLs. Intriguingly, at the 1p36.1 BMI GWAS locus the EPHB2 cis-aQTL was stronger than its cis-eQTL and colocalized with the GWAS signal and 35 BMI MR trans-aQTLs, suggesting the GWAS signal may be mediated by effects on EPHB2 activity and its downstream effects on a network of BMI MRs. These MR and aQTL analyses represent systems genetic methods that may be broadly applied to supplement standard eQTL analyses for suggesting molecular effects mediating GWAS signals., Competing Interests: The authors have declared that no competing interests exist.

Genome-wide association studies (GWASs) have discovered 20 risk loci in the human genome where germline variants associate with risk of pancreatic ductal adenocarcinoma (PDAC) in populations of European ancestry. Here, we fine-mapped one such locus on chr16q23.1 (rs72802365, p = 2.51 × 10 -17 , OR = 1.36, 95% CI = 1.31-1.40) and identified colocalization (PP = 0.87) with aberrant exon 5-7 CTRB2 splicing in pancreatic tissues (p GTEx = 1.40 × 10 -69 , β GTEx = 1.99; p LTG = 1.02 × 10 -30 , β LTG = 1.99). Imputation of a 584 bp structural variant overlapping exon 6 of CTRB2 into the GWAS datasets resulted in a highly significant association with pancreatic cancer risk (p = 2.83 × 10 -16 , OR = 1.36, 95% CI = 1.31-1.42), indicating that it may underlie this signal. Exon skipping attributable to the deletion (risk) allele introduces a premature stop codon in exon 7 of CTRB2, yielding a truncated chymotrypsinogen B2 protein that lacks chymotrypsin activity, is poorly secreted, and accumulates intracellularly in the endoplasmic reticulum (ER). We propose that intracellular accumulation of a nonfunctional chymotrypsinogen B2 protein leads to ER stress and pancreatic inflammation, which may explain the increased pancreatic cancer risk in carriers of CTRB2 exon 6 deletion alleles., Competing Interests: Declaration of interests The authors declare no competing interests., (Published by Elsevier Inc.)

Background: Epidemiological studies have suggested positive associations for iron and red meat intake with risk of pancreatic ductal adenocarcinoma (PDAC). Inherited pathogenic variants in genes involved in the hepcidin-regulating iron metabolism pathway are known to cause iron overload and hemochromatosis., Objectives: The objective of this study was to determine whether common genetic variation in the hepcidin-regulating iron metabolism pathway is associated with PDAC., Methods: We conducted a pathway analysis of the hepcidin-regulating genes using single nucleotide polymorphism (SNP) summary statistics generated from 4 genome-wide association studies in 2 large consortium studies using the summary data-based adaptive rank truncated product method. Our population consisted of 9253 PDAC cases and 12,525 controls of European descent. Our analysis included 11 hepcidin-regulating genes [bone morphogenetic protein 2 (BMP2), bone morphogenetic protein 6 (BMP6), ferritin heavy chain 1 (FTH1), ferritin light chain (FTL), hepcidin (HAMP), homeostatic iron regulator (HFE), hemojuvelin (HJV), nuclear factor erythroid 2-related factor 2 (NRF2), ferroportin 1 (SLC40A1), transferrin receptor 1 (TFR1), and transferrin receptor 2 (TFR2)] and their surrounding genomic regions (±20 kb) for a total of 412 SNPs., Results: The hepcidin-regulating gene pathway was significantly associated with PDAC (P = 0.002), with the HJV, TFR2, TFR1, BMP6, and HAMP genes contributing the most to the association., Conclusions: Our results support that genetic susceptibility related to the hepcidin-regulating gene pathway is associated with PDAC risk and suggest a potential role of iron metabolism in pancreatic carcinogenesis. Further studies are needed to evaluate effect modification by intake of iron-rich foods on this association., (Published by Oxford University Press on behalf of the American Society for Nutrition 2021.)

Genome-wide association studies (GWASs) have identified thousands of cancer risk loci revealing many risk regions shared across multiple cancers. Characterizing the cross-cancer shared genetic basis can increase our understanding of global mechanisms of cancer development. In this study, we collected GWAS summary statistics based on up to 375,468 cancer cases and 530,521 controls for fourteen types of cancer, including breast (overall, estrogen receptor [ER]-positive, and ER-negative), colorectal, endometrial, esophageal, glioma, head/neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancer, to characterize the shared genetic basis of cancer risk. We identified thirteen pairs of cancers with statistically significant local genetic correlations across eight distinct genomic regions. Specifically, the 5p15.33 region, harboring the TERT and CLPTM1L genes, showed statistically significant local genetic correlations for multiple cancer pairs. We conducted a cross-cancer fine-mapping of the 5p15.33 region based on eight cancers that showed genome-wide significant associations in this region (ER-negative breast, colorectal, glioma, lung, melanoma, ovarian, pancreatic, and prostate cancer). We used an iterative analysis pipeline implementing a subset-based meta-analysis approach based on cancer-specific conditional analyses and identified ten independent cross-cancer associations within this region. For each signal, we conducted cross-cancer fine-mapping to prioritize the most plausible causal variants. Our findings provide a more in-depth understanding of the shared inherited basis across human cancers and expand our knowledge of the 5p15.33 region in carcinogenesis.

Germline variation and smoking are independently associated with pancreatic ductal adenocarcinoma (PDAC). We conducted genome-wide smoking interaction analysis of PDAC using genotype data from four previous genome-wide association studies in individuals of European ancestry (7,937 cases and 11,774 controls). Examination of expression quantitative trait loci data from the Genotype-Tissue Expression Project followed by colocalization analysis was conducted to determine whether there was support for common SNP(s) underlying the observed associations. Statistical tests were two sided and P < 5 × 10 -8 was considered statistically significant. Genome-wide significant evidence of qualitative interaction was identified on chr2q21.3 in intron 5 of the transmembrane protein 163 (TMEM163) and upstream of the cyclin T2 (CCNT2). The most significant SNP using the Empirical Bayes method, in this region that included 45 significantly associated SNPs, was rs1818613 [per allele OR in never smokers 0.87, 95% confidence interval (CI), 0.82-0.93; former smokers 1.00, 95% CI, 0.91-1.07; current smokers 1.25, 95% CI 1.12-1.40, P interaction = 3.08 × 10 -9 ). Examination of the Genotype-Tissue Expression Project data demonstrated an expression quantitative trait locus in this region for TMEM163 and CCNT2 in several tissue types. Colocalization analysis supported a shared SNP, rs842357, in high linkage disequilibrium with rs1818613 ( r 2 = 0. 94) driving both the observed interaction and the expression quantitative trait loci signals. Future studies are needed to confirm and understand the differential biologic mechanisms by smoking status that contribute to our PDAC findings. SIGNIFICANCE: This large genome-wide interaction study identifies a susceptibility locus on 2q21.3 that significantly modified PDAC risk by smoking status, providing insight into smoking-associated PDAC, with implications for prevention., (©2021 American Association for Cancer Research.)

Jews are estimated to be at increased risk of pancreatic cancer compared to non-Jews, but their observed 50-80% excess risk is not explained by known non-genetic or genetic risk factors. We conducted a GWAS in a case-control sample of American Jews, largely Ashkenazi, including 406 pancreatic cancer patients and 2332 controls, identified in the dbGaP, PanScan I/II, PanC4 and GERA data sets. We then examined resulting SNPs with P < 10 -7 in an expanded sample set, of 539 full- or part-Jewish pancreatic cancer patients and 4117 full- or part-Jewish controls from the same data sets. Jewish ancestries were genetically determined using seeded FastPCA. Among the full Jews, a novel genome-wide significant association was detected on chromosome 19p12 (rs66562280, per-allele OR = 1.55, 95% CI = 1.33-1.81, P = 10 -7.6 ). A suggestive relatively independent association was detected on chromosome 19p13.3 (rs2656937, OR = 1.53, 95% CI = 1.31-1.78, P = 10 -7.0 ). Similar associations were seen for these SNPs among the full and part Jews combined. This is the first GWAS conducted for pancreatic cancer in the increased-risk Jewish population. The SNPs rs66562280 and rs2656937 are located in introns of ZNF100-like and ARRDC5, respectively, and are known to alter regulatory motifs of genes that play integral roles in pancreatic carcinogenesis.

Background: Whether circulating polyunsaturated fatty acid (PUFA) levels are associated with pancreatic cancer risk is uncertain. Mendelian randomization (MR) represents a study design using genetic instruments to better characterize the relationship between exposure and outcome., Methods: We utilized data from genome-wide association studies within the Pancreatic Cancer Cohort Consortium and Pancreatic Cancer Case-Control Consortium, involving approximately 9,269 cases and 12,530 controls of European descent, to evaluate associations between pancreatic cancer risk and genetically predicted plasma n-6 PUFA levels. Conventional MR analyses were performed using individual-level and summary-level data., Results: Using genetic instruments, we did not find evidence of associations between genetically predicted plasma n-6 PUFA levels and pancreatic cancer risk [estimates per one SD increase in each PUFA-specific weighted genetic score using summary statistics: linoleic acid odds ratio (OR) = 1.00, 95% confidence interval (CI) = 0.98-1.02; arachidonic acid OR = 1.00, 95% CI = 0.99-1.01; and dihomo-gamma-linolenic acid OR = 0.95, 95% CI = 0.87-1.02]. The OR estimates remained virtually unchanged after adjustment for covariates, using individual-level data or summary statistics, or stratification by age and sex., Conclusions: Our results suggest that variations of genetically determined plasma n-6 PUFA levels are not associated with pancreatic cancer risk., Impact: These results suggest that modifying n-6 PUFA levels through food sources or supplementation may not influence risk of pancreatic cancer., (©2020 American Association for Cancer Research.)

Background: Although 20 pancreatic cancer susceptibility loci have been identified through genome-wide association studies in individuals of European ancestry, much of its heritability remains unexplained and the genes responsible largely unknown., Methods: To discover novel pancreatic cancer risk loci and possible causal genes, we performed a pancreatic cancer transcriptome-wide association study in Europeans using three approaches: FUSION, MetaXcan, and Summary-MulTiXcan. We integrated genome-wide association studies summary statistics from 9040 pancreatic cancer cases and 12 496 controls, with gene expression prediction models built using transcriptome data from histologically normal pancreatic tissue samples (NCI Laboratory of Translational Genomics [n = 95] and Genotype-Tissue Expression v7 [n = 174] datasets) and data from 48 different tissues (Genotype-Tissue Expression v7, n = 74-421 samples)., Results: We identified 25 genes whose genetically predicted expression was statistically significantly associated with pancreatic cancer risk (false discovery rate < .05), including 14 candidate genes at 11 novel loci (1p36.12: CELA3B; 9q31.1: SMC2, SMC2-AS1; 10q23.31: RP11-80H5.9; 12q13.13: SMUG1; 14q32.33: BTBD6; 15q23: HEXA; 15q26.1: RCCD1; 17q12: PNMT, CDK12, PGAP3; 17q22: SUPT4H1; 18q11.22: RP11-888D10.3; and 19p13.11: PGPEP1) and 11 at six known risk loci (5p15.33: TERT, CLPTM1L, ZDHHC11B; 7p14.1: INHBA; 9q34.2: ABO; 13q12.2: PDX1; 13q22.1: KLF5; and 16q23.1: WDR59, CFDP1, BCAR1, TMEM170A). The association for 12 of these genes (CELA3B, SMC2, and PNMT at novel risk loci and TERT, CLPTM1L, INHBA, ABO, PDX1, KLF5, WDR59, CFDP1, and BCAR1 at known loci) remained statistically significant after Bonferroni correction., Conclusions: By integrating gene expression and genotype data, we identified novel pancreatic cancer risk loci and candidate functional genes that warrant further investigation., (Published by Oxford University Press 2020. This work is written by US Government employees and is in the public domain in the US.)

Registry-based epidemiologic studies suggest associations between chronic inflammatory intestinal diseases and pancreatic ductal adenocarcinoma (PDAC). As genetic susceptibility contributes to a large proportion of chronic inflammatory intestinal diseases, we hypothesize that the genomic regions surrounding established genome-wide associated variants for these chronic inflammatory diseases are associated with PDAC. We examined the association between PDAC and genomic regions (±500 kb) surrounding established common susceptibility variants for ulcerative colitis, Crohn's disease, inflammatory bowel disease, celiac disease, chronic pancreatitis, and primary sclerosing cholangitis. We analyzed summary statistics from genome-wide association studies data for 8,384 cases and 11,955 controls of European descent from two large consortium studies using the summary data-based adaptive rank truncated product method to examine the overall association of combined genomic regions for each inflammatory disease group. Combined genomic susceptibility regions for ulcerative colitis, Crohn disease, inflammatory bowel disease, and chronic pancreatitis were associated with PDAC at P values < 0.05 (0.0040, 0.0057, 0.011, and 3.4 × 10 -6 , respectively). After excluding the 20 PDAC susceptibility regions (±500 kb) previously identified by GWAS, the genomic regions for ulcerative colitis, Crohn disease, and inflammatory bowel disease remained associated with PDAC ( P = 0.0029, 0.0057, and 0.0098, respectively). Genomic regions for celiac disease ( P = 0.22) and primary sclerosing cholangitis ( P = 0.078) were not associated with PDAC. Our results support the hypothesis that genomic regions surrounding variants associated with inflammatory intestinal diseases, particularly, ulcerative colitis, Crohn disease, inflammatory bowel disease, and chronic pancreatitis are associated with PDAC. SIGNIFICANCE: The joint effects of common variants in genomic regions containing susceptibility loci for inflammatory bowel disease and chronic pancreatitis are associated with PDAC and may provide insights to understanding pancreatic cancer etiology., (©2020 American Association for Cancer Research.)

Background: Obesity and diabetes are major modifiable risk factors for pancreatic cancer. Interactions between genetic variants and diabetes/obesity have not previously been comprehensively investigated in pancreatic cancer at the genome-wide level., Methods: We conducted a gene-environment interaction (GxE) analysis including 8,255 cases and 11,900 controls from four pancreatic cancer genome-wide association study (GWAS) datasets (Pancreatic Cancer Cohort Consortium I-III and Pancreatic Cancer Case Control Consortium). Obesity (body mass index ≥30 kg/m 2 ) and diabetes (duration ≥3 years) were the environmental variables of interest. Approximately 870,000 SNPs (minor allele frequency ≥0.005, genotyped in at least one dataset) were analyzed. Case-control (CC), case-only (CO), and joint-effect test methods were used for SNP-level GxE analysis. As a complementary approach, gene-based GxE analysis was also performed. Age, sex, study site, and principal components accounting for population substructure were included as covariates. Meta-analysis was applied to combine individual GWAS summary statistics., Results: No genome-wide significant interactions (departures from a log-additive odds model) with diabetes or obesity were detected at the SNP level by the CC or CO approaches. The joint-effect test detected numerous genome-wide significant GxE signals in the GWAS main effects top hit regions, but the significance diminished after adjusting for the GWAS top hits. In the gene-based analysis, a significant interaction of diabetes with variants in the FAM63A (family with sequence similarity 63 member A) gene (significance threshold P < 1.25 × 10 -6 ) was observed in the meta-analysis ( P GxE = 1.2 ×10 -6 , P Joint = 4.2 ×10 -7 )., Conclusions: This analysis did not find significant GxE interactions at the SNP level but found one significant interaction with diabetes at the gene level. A larger sample size might unveil additional genetic factors via GxE scans., Impact: This study may contribute to discovering the mechanism of diabetes-associated pancreatic cancer., (©2020 American Association for Cancer Research.)

Genome-wide association studies (GWAS) have led to the identification of hundreds of susceptibility loci across cancers, but the impact of further studies remains uncertain. Here we analyse summary-level data from GWAS of European ancestry across fourteen cancer sites to estimate the number of common susceptibility variants (polygenicity) and underlying effect-size distribution. All cancers show a high degree of polygenicity, involving at a minimum of thousands of loci. We project that sample sizes required to explain 80% of GWAS heritability vary from 60,000 cases for testicular to over 1,000,000 cases for lung cancer. The maximum relative risk achievable for subjects at the 99th risk percentile of underlying polygenic risk scores (PRS), compared to average risk, ranges from 12 for testicular to 2.5 for ovarian cancer. We show that PRS have potential for risk stratification for cancers of breast, colon and prostate, but less so for others because of modest heritability and lower incidence.

Pancreatic cancer is the fourth leading cause of cancer-related deaths in Japan. To identify risk loci, we perform a meta-analysis of three genome-wide association studies comprising 2,039 pancreatic cancer patients and 32,592 controls in the Japanese population. Here, we identify 3 (13q12.2, 13q22.1, and 16p12.3) genome-wide significant loci (P < 5.0 × 10 -8 ), of which 16p12.3 has not been reported in the Western population. The lead single nucleotide polymorphism (SNP) at 16p12.3 is rs78193826 (odds ratio = 1.46, 95% confidence interval = 1.29-1.66, P = 4.28 × 10 -9 ), an Asian-specific, nonsynonymous glycoprotein 2 (GP2) gene variant. Associations between selected GP2 gene variants and pancreatic cancer are replicated in 10,822 additional cases and controls of East Asian origin. Functional analyses using cell lines provide supporting evidence of the effect of rs78193826 on KRAS activity. These findings suggest that GP2 gene variants are probably associated with pancreatic cancer susceptibility in populations of East Asian ancestry.

Background: Pancreatic cancer is the third leading cause of cancer death in the United States, and 80% of patients present with advanced, incurable disease. Risk markers for pancreatic cancer have been characterized, but combined models are not used clinically to identify individuals at high risk for the disease., Methods: Within a nested case-control study of 500 pancreatic cancer cases diagnosed after blood collection and 1,091 matched controls enrolled in four U.S. prospective cohorts, we characterized absolute risk models that included clinical factors (e.g., body mass index, history of diabetes), germline genetic polymorphisms, and circulating biomarkers., Results: Model discrimination showed an area under ROC curve of 0.62 via cross-validation. Our final integrated model identified 3.7% of men and 2.6% of women who had at least 3 times greater than average risk in the ensuing 10 years. Individuals within the top risk percentile had a 4% risk of developing pancreatic cancer by age 80 years and 2% 10-year risk at age 70 years., Conclusions: Risk models that include established clinical, genetic, and circulating factors improved disease discrimination over models using clinical factors alone., Impact: Absolute risk models for pancreatic cancer may help identify individuals in the general population appropriate for disease interception., (©2020 American Association for Cancer Research.)

Background: Leukocyte telomere length has been associated with risk of subsequent pancreatic cancer. Few prospective studies have evaluated the association of prediagnostic leukocyte telomere length with pancreatic cancer survival., Methods: We prospectively examined the association of prediagnostic leukocyte telomere length with overall survival (OS) time among 423 participants diagnosed with pancreatic adenocarcinoma between 1984 and 2008 within the Health Professionals Follow-up Study, Nurses' Health Study, Physicians' Health Study, and Women's Health Initiative. We measured prediagnostic leukocyte telomere length in banked blood samples using quantitative PCR. Cox proportional hazards models were used to estimate HRs for OS with adjustment for potential confounders. We also evaluated 10 SNPs at the telomerase reverse transcriptase locus., Results: Shorter prediagnostic leukocyte telomere length was associated with reduced OS among patients with pancreatic cancer ( P trend = 0.04). The multivariable-adjusted HR for OS comparing the lowest with highest quintiles of leukocyte telomere length was 1.39 (95% confidence interval, 1.01-1.93), corresponding to a 3-month difference in median OS time. In an analysis excluding cases with blood collected within 2 years of cancer diagnosis, the association was moderately stronger (HR, 1.55; 95% confidence interval, 1.09-2.21; comparing the lowest with highest quintiles; P trend = 0.01). No prognostic association or effect modification for the prognostic association of prediagnostic leukocyte telomere length was noted in relation to the studied SNPs., Conclusions: Prediagnostic leukocyte telomere length was associated with pancreatic cancer survival., Impact: Prediagnostic leukocyte telomere length can be a prognostic biomarker in pancreatic cancer., (©2019 American Association for Cancer Research.)

Objective: The KRAS gene is the most frequently mutated gene in pancreatic cancer, and no successful anti-Ras therapy has been developed. Gastrin has been shown to stimulate pancreatic cancer in an autocrine fashion. We hypothesized that reactivation of the peptide gastrin collaborates with KRAS during pancreatic carcinogenesis., Methods: LSL-Kras; P48-Cre (KC) mutant KRAS transgenic mice were crossed with gastrin-KO (GKO) mice to develop GKO/KC mice. Pancreata were examined for 8 months for stage of pancreatic intraepithelial neoplasia lesions, inflammation, fibrosis, gastrin peptide, and microRNA expression. Pancreatic intraepithelial neoplasias from mice were collected by laser capture microdissection and subjected to reverse-phase protein microarray, for gastrin and protein kinases associated with signal transduction. Gastrin mRNA was measured by RNAseq in human pancreatic cancer tissues and compared to that in normal pancreas., Results: In the absence of gastrin, PanIN progression, inflammation, and fibrosis were significantly decreased and signal transduction was reversed to the canonical pathway with decreased KRAS. Gastrin re-expression in the PanINs was mediated by miR-27a. Gastrin mRNA expression was significantly increased in human pancreatic cancer samples compared to normal human pancreas controls., Conclusions: This study supports the mitogenic role of gastrin in activation of KRAS during pancreatic carcinogenesis.

Background: Genome-wide association studies (GWAS) identify associations of individual single-nucleotide polymorphisms (SNPs) with cancer risk but usually only explain a fraction of the inherited variability. Pathway analysis of genetic variants is a powerful tool to identify networks of susceptibility genes., Methods: We conducted a large agnostic pathway-based meta-analysis of GWAS data using the summary-based adaptive rank truncated product method to identify gene sets and pathways associated with pancreatic ductal adenocarcinoma (PDAC) in 9040 cases and 12 496 controls. We performed expression quantitative trait loci (eQTL) analysis and functional annotation of the top SNPs in genes contributing to the top associated pathways and gene sets. All statistical tests were two-sided., Results: We identified 14 pathways and gene sets associated with PDAC at a false discovery rate of less than 0.05. After Bonferroni correction (P ≤ 1.3 × 10-5), the strongest associations were detected in five pathways and gene sets, including maturity-onset diabetes of the young, regulation of beta-cell development, role of epidermal growth factor (EGF) receptor transactivation by G protein-coupled receptors in cardiac hypertrophy pathways, and the Nikolsky breast cancer chr17q11-q21 amplicon and Pujana ATM Pearson correlation coefficient (PCC) network gene sets. We identified and validated rs876493 and three correlating SNPs (PGAP3) and rs3124737 (CASP7) from the Pujana ATM PCC gene set as eQTLs in two normal derived pancreas tissue datasets., Conclusion: Our agnostic pathway and gene set analysis integrated with functional annotation and eQTL analysis provides insight into genes and pathways that may be biologically relevant for risk of PDAC, including those not previously identified., (Published by Oxford University Press 2018. This work is written by US Government employees and is in the public domain in the US.)

Most expression quantitative trait locus (eQTL) studies to date have been performed in heterogeneous tissues as opposed to specific cell types. To better understand the cell-type-specific regulatory landscape of human melanocytes, which give rise to melanoma but account for <5% of typical human skin biopsies, we performed an eQTL analysis in primary melanocyte cultures from 106 newborn males. We identified 597,335 cis -eQTL SNPs prior to linkage disequilibrium (LD) pruning and 4997 eGenes (FDR < 0.05). Melanocyte eQTLs differed considerably from those identified in the 44 GTEx tissue types, including skin. Over a third of melanocyte eGenes, including key genes in melanin synthesis pathways, were unique to melanocytes compared to those of GTEx skin tissues or TCGA melanomas. The melanocyte data set also identified trans -eQTLs, including those connecting a pigmentation-associated functional SNP with four genes, likely through cis -regulation of IRF4 Melanocyte eQTLs are enriched in cis -regulatory signatures found in melanocytes as well as in melanoma-associated variants identified through genome-wide association studies. Melanocyte eQTLs also colocalized with melanoma GWAS variants in five known loci. Finally, a transcriptome-wide association study using melanocyte eQTLs uncovered four novel susceptibility loci, where imputed expression levels of five genes ( ZFP90 , HEBP1 , MSC , CBWD1 , and RP11-383H13.1 ) were associated with melanoma at genome-wide significant P -values. Our data highlight the utility of lineage-specific eQTL resources for annotating GWAS findings, and present a robust database for genomic research of melanoma risk and melanocyte biology., (© 2018 Zhang et al.; Published by Cold Spring Harbor Laboratory Press.)

Diets with high inflammatory potential are suspected to increase risk for pancreatic cancer (PC). Using pooled analyses, we examined whether this association applies to populations from different geographic regions and population subgroups with varying risks for PC, including variation in ABO blood type. Data from six case-control studies (cases, n = 2414; controls, n = 4528) in the Pancreatic Cancer Case-Control Consortium (PanC4) were analyzed, followed by replication in five nested case-control studies (cases, n = 1268; controls, n = 4215) from the Pancreatic Cancer Cohort Consortium (PanScan). Two polymorphisms in the ABO locus (rs505922 and rs8176746) were used to infer participants' blood types. Dietary questionnaire-derived nutrient/food intake was used to compute energy-adjusted dietary inflammatory index (E-DII®) scores to assess inflammatory potential of diet. Pooled odds ratios (ORs) and 95% confidence intervals (CIs) were calculated using multivariable-adjusted logistic regression. Higher E-DII scores, reflecting greater inflammatory potential of diet, were associated with increased PC risk in PanC4 [ORQ5 versus Q1=2.20, 95% confidence interval (CI) = 1.85-2.61, Ptrend < 0.0001; ORcontinuous = 1.20, 95% CI = 1.17-1.24], and PanScan (ORQ5 versus Q1 = 1.23, 95% CI = 0.92-1.66, Ptrend = 0.008; ORcontinuous = 1.09, 95% CI = 1.02-1.15). As expected, genotype-derived non-O blood type was associated with increased PC risk in both the PanC4 and PanScan studies. Stratified analyses of associations between E-DII quintiles and PC by genotype-derived ABO blood type did not show interaction by blood type (Pinteraction = 0.10 in PanC4 and Pinteraction=0.13 in PanScan). The results show that consuming a pro-inflammatory diet and carrying non-O blood type are each individually, but not interactively, associated with increased PC risk., (© The Author(s) 2018. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

This corrects the article DOI: 10.1038/ncomms15034.

Chronic inflammation increases the risk of developing one of several types of cancer. Inflammatory responses are currently thought to be controlled by mechanisms that rely on transcriptional networks that are distinct from those involved in cell differentiation. The orphan nuclear receptor NR5A2 participates in a wide variety of processes, including cholesterol and glucose metabolism in the liver, resolution of endoplasmic reticulum stress, intestinal glucocorticoid production, pancreatic development and acinar differentiation. In genome-wide association studies, single nucleotide polymorphisms in the vicinity of NR5A2 have previously been associated with the risk of pancreatic adenocarcinoma. In mice, Nr5a2 heterozygosity sensitizes the pancreas to damage, impairs regeneration and cooperates with mutant Kras in tumour progression. Here, using a global transcriptomic analysis, we describe an epithelial-cell-autonomous basal pre-inflammatory state in the pancreas of Nr5a2 +/- mice that is reminiscent of the early stages of pancreatitis-induced inflammation and is conserved in histologically normal human pancreases with reduced expression of NR5A2 mRNA. In Nr5a2 +/- mice, NR5A2 undergoes a marked transcriptional switch, relocating from differentiation-specific to inflammatory genes and thereby promoting gene transcription that is dependent on the AP-1 transcription factor. Pancreatic deletion of Jun rescues the pre-inflammatory phenotype, as well as binding of NR5A2 to inflammatory gene promoters and the defective regenerative response to damage. These findings support the notion that, in the pancreas, the transcriptional networks involved in differentiation-specific functions also suppress inflammatory programmes. Under conditions of genetic or environmental constraint, these networks can be subverted to foster inflammation.

In 2020, 146,063 deaths due to pancreatic cancer are estimated to occur in Europe and the United States combined. To identify common susceptibility alleles, we performed the largest pancreatic cancer GWAS to date, including 9040 patients and 12,496 controls of European ancestry from the Pancreatic Cancer Cohort Consortium (PanScan) and the Pancreatic Cancer Case-Control Consortium (PanC4). Here, we find significant evidence of a novel association at rs78417682 (7p12/TNS3, P = 4.35 × 10 -8 ). Replication of 10 promising signals in up to 2737 patients and 4752 controls from the PANcreatic Disease ReseArch (PANDoRA) consortium yields new genome-wide significant loci: rs13303010 at 1p36.33 (NOC2L, P = 8.36 × 10 -14 ), rs2941471 at 8q21.11 (HNF4G, P = 6.60 × 10 -10 ), rs4795218 at 17q12 (HNF1B, P = 1.32 × 10 -8 ), and rs1517037 at 18q21.32 (GRP, P = 3.28 × 10 -8 ). rs78417682 is not statistically significantly associated with pancreatic cancer in PANDoRA. Expression quantitative trait locus analysis in three independent pancreatic data sets provides molecular support of NOC2L as a pancreatic cancer susceptibility gene.

Background: The higher risk of pancreatic cancer in Ashkenazi Jews compared with non-Jews is only partially explained by the increased frequency of BRCA1 and BRCA2 mutations in Ashkenazi Jews. Methods: We evaluated the impact of 16 established pancreatic cancer susceptibility loci in a case-control sample of American Jews, largely Ashkenazi, including 406 full-Jewish pancreatic cancer patients and 2,332 full-Jewish controls, genotyped as part of the Pancreatic Cancer Cohort and Case-Control Consortium I/II (PanScan I/II), Pancreatic Cancer Case-Control Consortium (PanC4), and Resource for Genetic Epidemiology Research on Adult Health and Aging (GERA) datasets. We compared risk in full-Jewish subjects with risk in part-Jewish; non-Jewish Southern European; and in the combined non-Jewish Eastern, Northern, Southern, and Western European (non-Jewish white European) subjects from the same datasets. Jewish ancestries were genetically identified using seeded Fast principal component analysis. Data were analyzed by unconditional logistic regression, and adjusted for age, sex, and principal components. Results: One SNP on chromosome 13q22.1 (rs9543325; OR, 1.36; 95% confidence interval, 1.16-1.58; P = 10 -4.1 ) was significant in full-Jews. Individual ORs and minor allele frequencies were similar between Jewish and non-Jewish white European subjects. The average ORs across the 16 pancreatic cancer susceptibility loci for full-Jewish, full- plus part-Jewish, non-Jewish Southern European, and non-Jewish white European subjects were 1.25, 1.30, 1.31, and 1.26, respectively. Conclusions: The 16 pancreatic cancer susceptibility loci similarly impact Jewish and non-Jewish white European subjects, both individually and as summary odds. Impact: These 16 pancreatic cancer susceptibility loci likely do not explain the higher risk seen in Ashkenazi Jews. Cancer Epidemiol Biomarkers Prev; 26(10); 1540-8. ©2017 AACR ., (©2017 American Association for Cancer Research.)

Background: Many cancers share specific genetic risk factors, including both rare high-penetrance mutations and common SNPs identified through genome-wide association studies (GWAS). However, little is known about the overall shared heritability across cancers. Quantifying the extent to which two distinct cancers share genetic origin will give insights to shared biological mechanisms underlying cancer and inform design for future genetic association studies. Methods: In this study, we estimated the pair-wise genetic correlation between six cancer types (breast, colorectal, lung, ovarian, pancreatic, and prostate) using cancer-specific GWAS summary statistics data based on 66,958 case and 70,665 control subjects of European ancestry. We also estimated genetic correlations between cancers and 14 noncancer diseases and traits. Results: After adjusting for 15 pair-wise genetic correlation tests between cancers, we found significant ( P < 0.003) genetic correlations between pancreatic and colorectal cancer (rg = 0.55, P = 0.003), lung and colorectal cancer (rg = 0.31, P = 0.001). We also found suggestive genetic correlations between lung and breast cancer (rg = 0.27, P = 0.009), and colorectal and breast cancer (rg = 0.22, P = 0.01). In contrast, we found no evidence that prostate cancer shared an appreciable proportion of heritability with other cancers. After adjusting for 84 tests studying genetic correlations between cancer types and other traits (Bonferroni-corrected P value: 0.0006), only the genetic correlation between lung cancer and smoking remained significant (rg = 0.41, P = 1.03 × 10 -6 ). We also observed nominally significant genetic correlations between body mass index and all cancers except ovarian cancer. Conclusions: Our results highlight novel genetic correlations and lend support to previous observational studies that have observed links between cancers and risk factors. Impact: This study demonstrates modest genetic correlations between cancers; in particular, breast, colorectal, and lung cancer share some degree of genetic basis. Cancer Epidemiol Biomarkers Prev; 26(9); 1427-35. ©2017 AACR ., (©2017 American Association for Cancer Research.)

Objective: Telomere shortening occurs as an early event in pancreatic tumorigenesis, and genetic variants at the telomerase reverse transcriptase ( TERT ) gene region have been associated with pancreatic cancer risk. However, it is unknown whether prediagnostic leucocyte telomere length is associated with subsequent risk of pancreatic cancer., Design: We measured prediagnostic leucocyte telomere length in 386 pancreatic cancer cases and 896 matched controls from five prospective US cohorts. ORs and 95% CIs were calculated using conditional logistic regression. Matching factors included year of birth, cohort (which also matches on sex), smoking status, fasting status and month/year of blood collection. We additionally examined single-nucleotide polymorphisms (SNPs) at the TERT region in relation to pancreatic cancer risk and leucocyte telomere length using logistic and linear regression, respectively., Results: Shorter prediagnostic leucocyte telomere length was associated with higher risk of pancreatic cancer (comparing extreme quintiles of telomere length, OR 1.72; 95% CI 1.07 to 2.78; p trend =0.048). Results remained unchanged after adjustment for diabetes, body mass index and physical activity. Three SNPs at TERT (linkage disequilibrium r 2 <0.25) were associated with pancreatic cancer risk, including rs401681 (per minor allele OR 1.33; 95% CI 1.12 to 1.59; p=0.002), rs2736100 (per minor allele OR 1.36; 95% CI 1.13 to 1.63; p=0.001) and rs2736098 (per minor allele OR 0.75; 95% CI 0.63 to 0.90; p=0.002). The minor allele for rs401681 was associated with shorter telomere length (p=0.023)., Conclusions: Prediagnostic leucocyte telomere length and genetic variants at the TERT gene region were associated with risk of pancreatic cancer., Competing Interests: Competing interests: None declared., (Published by the BMJ Publishing Group Limited. For permission to use (where not already granted under a licence) please go to http://www.bmj.com/company/products-services/rights-and-licensing/.)

Background: Chromosome 8q24 has emerged as an important genetic susceptibility region for several cancers, including prostate cancer; however, little is known about the contribution of DNA methylation in this region to risk., Methods: We prospectively evaluated DNA methylation at 8q24 in relation to prostate cancer using pre-diagnostic blood samples from 694 prostate cancer cases (including 172 aggressive cases) and 703 controls in the Prostate, Lung, Colorectal and Ovarian Cancer Screening Trial. We used logistic regression to estimate odds ratios and 95% confidence intervals., Results: Although none remained significant after adjustment for multiple testing (q>0.05), of the 50 CpG sites meeting quality control, we identified 8 sites that were nominally associated with prostate cancer (P trend <0.05), including 6 correlated (Spearman ρ: 0.20-0.52) sites in POU5F1B and 2 intergenic sites (most significant site: Chr8:128428897 in POU5F1B, P trend =0.01). We also identified two correlated (ρ=0.39) sites in MYC (Chr8:128753187 and Chr8:128753154) that were associated with aggressive (P trend =0.02 and 0.03), but not non-aggressive disease (P trend =0.70 and 0.20; P heterogeneity =0.01 and 4.6 × 10 -3 ). These findings persisted after adjustment for the top 8q24 prostate cancer variants in our study., Conclusions: Although requiring replication, our findings provide some evidence that 8q24 DNA methylation levels may be associated with prostate cancer risk.

Genome wide association studies (GWAS) have mapped multiple independent cancer susceptibility loci to chr5p15.33. Here, we show that fine-mapping of pancreatic and testicular cancer GWAS within one of these loci (Region 2 in CLPTM1L) focuses the signal to nine highly correlated SNPs. Of these, rs36115365-C associated with increased pancreatic and testicular but decreased lung cancer and melanoma risk, and exhibited preferred protein-binding and enhanced regulatory activity. Transcriptional gene silencing of this regulatory element repressed TERT expression in an allele-specific manner. Proteomic analysis identifies allele-preferred binding of Zinc finger protein 148 (ZNF148) to rs36115365-C, further supported by binding of purified recombinant ZNF148. Knockdown of ZNF148 results in reduced TERT expression, telomerase activity and telomere length. Our results indicate that the association with chr5p15.33-Region 2 may be explained by rs36115365, a variant influencing TERT expression via ZNF148 in a manner consistent with elevated TERT in carriers of the C allele.

Recent heritability analyses have indicated that genome-wide association studies (GWAS) have the potential to improve genetic risk prediction for complex diseases based on polygenic risk score (PRS), a simple modelling technique that can be implemented using summary-level data from the discovery samples. We herein propose modifications to improve the performance of PRS. We introduce threshold-dependent winner's-curse adjustments for marginal association coefficients that are used to weight the single-nucleotide polymorphisms (SNPs) in PRS. Further, as a way to incorporate external functional/annotation knowledge that could identify subsets of SNPs highly enriched for associations, we propose variable thresholds for SNPs selection. We applied our methods to GWAS summary-level data of 14 complex diseases. Across all diseases, a simple winner's curse correction uniformly led to enhancement of performance of the models, whereas incorporation of functional SNPs was beneficial only for selected diseases. Compared to the standard PRS algorithm, the proposed methods in combination led to notable gain in efficiency (25-50% increase in the prediction R2) for 5 of 14 diseases. As an example, for GWAS of type 2 diabetes, winner's curse correction improved prediction R2 from 2.29% based on the standard PRS to 3.10% (P = 0.0017) and incorporating functional annotation data further improved R2 to 3.53% (P = 2×10-5). Our simulation studies illustrate why differential treatment of certain categories of functional SNPs, even when shown to be highly enriched for GWAS-heritability, does not lead to proportionate improvement in genetic risk-prediction because of non-uniform linkage disequilibrium structure., Competing Interests: The authors have declared that no competing interests exist.

Leptin is an adipokine involved in regulating energy balance, which has been identified as a potential biologic link in the development of obesity-associated cancers, such as pancreatic cancer. In this prospective, nested case-control study of 470 cases and 1,094 controls from five U.S. cohorts, we used conditional logistic regression to evaluate pancreatic cancer risk by prediagnostic plasma leptin, adjusting for race/ethnicity, diabetes, body mass index, physical activity, plasma C-peptide, adiponectin, and 25-hydroxyvitamin D. Because of known differences in leptin levels by gender, analyses were conducted separately for men and women. We also evaluated associations between 32 tagging SNPs in the leptin receptor (LEPR) gene and pancreatic cancer risk. Leptin levels were higher in female versus male control participants (median, 20.8 vs. 6.7 ng/mL; P < 0.0001). Among men, plasma leptin was positively associated with pancreatic cancer risk and those in the top quintile had a multivariable-adjusted OR of 3.02 [95% confidence interval (CI), 1.27-7.16; P trend = 0.02] compared with men in the bottom quintile. Among women, circulating leptin was not associated with pancreatic cancer risk (P trend = 0.21). Results were similar across cohorts (P heterogeneity = 0.88 for two male cohorts and 0.35 for three female cohorts). In genetic analyses, rs10493380 in LEPR was associated with increased pancreatic cancer risk among women, with an OR per minor allele of 1.54 (95% CI, 1.18-2.02; multiple hypothesis-corrected P = 0.03). No SNPs were significantly associated with risk in men. In conclusion, higher prediagnostic levels of plasma leptin were associated with an elevated risk of pancreatic cancer among men, but not among women. Cancer Res; 76(24); 7160-7. ©2016 AACR., Competing Interests: of Potential Conflicts of Interest No potential conflicts of interest were disclosed, (©2016 American Association for Cancer Research.)

Genome-wide association studies (GWAS) have identified multiple common susceptibility loci for pancreatic cancer. Here we report fine-mapping and functional analysis of one such locus residing in a 610 kb gene desert on chr13q22.1 (marked by rs9543325). The closest candidate genes, KLF5, KLF12, PIBF1, DIS3 and BORA, range in distance from 265-586 kb. Sequencing three sub-regions containing the top ranked SNPs by imputation P-value revealed a 30 bp insertion/deletion (indel) variant that was significantly associated with pancreatic cancer risk (rs386772267, P = 2.30 × 10 -11 , OR = 1.22, 95% CI 1.15-1.28) and highly correlated to rs9543325 (r 2  =   0.97 in the 1000 Genomes EUR population). This indel was the most significant cis-eQTL variant in a set of 222 histologically normal pancreatic tissue samples (β = 0.26, P = 0.004), with the insertion (risk-increasing) allele associated with reduced DIS3 expression. DIS3 encodes a catalytic subunit of the nuclear RNA exosome complex that mediates RNA processing and decay, and is mutated in several cancers. Chromosome conformation capture revealed a long range (570 kb) physical interaction between a sub-region of the risk locus, containing rs386772267, and a region ∼6 kb upstream of DIS3 Finally, repressor regulatory activity and allele-specific protein binding by transcription factors of the TCF/LEF family were observed for the risk-increasing allele of rs386772267, indicating that expression regulation at this risk locus may be influenced by the Wnt signaling pathway. In conclusion, we have identified a putative functional indel variant at chr13q22.1 that associates with decreased DIS3 expression in carriers of pancreatic cancer risk-increasing alleles, and could therefore affect nuclear RNA processing and/or decay., (Published by Oxford University Press 2016. This work is written by US Government employees and is in the public domain in the US.)