Abstract 1194: Cross-cancer GWAS meta-analysis of more than 370,000 cases and 530,000 controls identifies multiple novel cancer risk regions

Genome-wide association studies (GWAS) have identified hundreds of common, low-penetrance alleles associated with cancer risk. However, known rare and common risk alleles only explain between 10% and 30% of the familial relative risk for different cancers and multiple lines of evidence indicate that many more risk alleles remain to be discovered. We have demonstrated genetic correlations between cancers, reflecting a shared genetic origin for solid tumors. These results suggest that jointly analyzing multiple cancer sites will lead to the discovery of novel risk regions.

We conducted a cross-cancer GWAS meta-analysis by leveraging GWAS summary statistics from 12 solid cancers (breast, colorectal, endometrial, esophageal, glioma, head and neck, lung, melanoma, ovarian, pancreatic, prostate and renal cancers) with a total of 373,818 cases and 532,382 controls of European ancestry. All studies had been imputed to either 1,000 Genomes or the Haplotype Reference Consortium panel. We conducted four meta-analysis using (1) fixed-effect, (2) random-effect, (3) one-sided subset (ASSET) and (4) two-sided subset (ASSET) models. The subset analysis were conducted assuming either the same direction of effects across cancers (one-sided ASSET) or allowed for opposite direction of effects across cancers (two-sided ASSET). In all analyses, we used tetrachoric correlations to account for sample overlap across cancer sites. In total, we tested 10,223,013 variants for association. We considered regions with a p-value<1.25 × 10−8 in at least one of the four meta-analysis approaches and located at least 500kb away from known cancer risk SNPs as novel.

We identified eight novel regions that reached genome-wide significance. Of those eight regions, two were identified from fixed-effects meta-analysis, three from random effects meta-analysis, one from the two-sided subset analysis, and two regions (15.q15.3 and 21q22.3) were identified at p<1.25 × 10−8 in three of the meta-analysis approaches. Among novel findings is a deleterious missense variant located in RREB1 previously associated with type 2 diabetes, a deleterious missense variant located in DSTYK previously associated with waist-to-hip ratio and triglycerides, and an intergenic variant in linkage disequilibrium (LD) with variants in TMEM18, previously associated with body mass index. Other potential target genes among the newly discovered regions include TP53BP1 and PCNT, both previously implicated in carcinogenesis.

Citation Format: Sara Lindström, Siddhartha Kar, Lu Wang, Constance Turman, James MacDonald, Theo Bammler, BCAC, OCAC, PRACTICAL, Jeroen Huyghe, Stephanie Schmit, Tracy A. O'Mara, Deborah J. Thompson, Puya Gharahkhani, Stuart MacGregor, Paul Brennan, Richard S. Houlston, Beatrice S. Melin, Christopher I. Amos, James McKay, Mark M. Iles, Matthew H. Law, Alison Klein, Laufey Amundadottir, Bogdan Pasaniuc, Paul Pharoah, Rayjean J. Hung, Peter Kraft. Cross-cancer GWAS meta-analysis of more than 370,000 cases and 530,000 controls identifies multiple novel cancer risk regions [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1194.

Hepcidin-regulating Iron-metabolism Genes and Pancreatic Ductal Adenocarcinoma: A Pathway Analysis of Genome-wide Association Studies

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal cancer, and epidemiological studies have suggested positive associations with iron and red meat intake. Rare mutations in genes involved in the hepcidin-regulating pathway are known to cause iron overload and hemochromatosis. We hypothesize that the hepcidin-regulating pathway as characterized by common variants from genome-wide association studies will be associated with PDAC. Methods: We conducted a large pathway-based meta-analysis of the hepcidin-regulating genes using the summary based adaptive rank truncated product (sARTP) method in 9,253 PDAC cases and 12,525 controls of European descent from the Pancreatic Cancer Cohort (PanScan) and the Pancreatic Cancer Case-Control (PANC4) consortia. Our analysis included 11 hepcidin-regulating genes (BMP2, BMP6, FTH1, FTL, HAMP, HFE, HJV, NRF2, SLC40A1, TFR1, TFR2) and adjacent genomic regions (20 kb upstream and downstream) with a total of 412 single-nucleotide polymorphisms (SNPs). We also conducted the sARTP with four iron status biomarkers (serum iron, transferrin, transferrin saturation, and ferritin, n = 23,986) using summary statistics from previous GWAS studies (Benyamin, et al. 2014) to examine if the hepcidin-regulating genes were also associated with these iron traits. The sARTP method combines SNP-level associations across variants in a gene or a pathway. Signals from up to five of the most associated SNPs for each gene studied were accumulated. Results: The hepcidin-regulating pathway was significantly associated with PDAC (P-value = 0.002) with the HJV, TFR2, and TFR1 genes contributing the most to the association (gene level P-values = 0.001, 0.014, and 0.019, respectively). The pathway associations were more significant in women than men. This pathway was also significantly associated with the four biomarkers of iron metabolism (P-values <1.5 × 10–7). Conclusions: Our results support that genetic susceptibility related to the hepcidin-regulating pathway is associated with PDAC and a potential role of iron metabolism in pancreatic carcinogenesis. Further studies are needed to evaluate the modifying effect of iron-rich foods and genetic susceptibility of this pathway and PDAC risk.

A Pathway Analysis of Hereditary Hemochromatosis-related Genes and Pancreatic Ductal Adenocarcinoma Risk (FS11-05-19)

Objectives

Hereditary primary hemochromatosis is characterized by dysregulation of iron homeostasis and is caused by a genetic predisposition to absorb too much iron from foods. Hemochromatosis has been associated with some chronic diseases, including hepatocellular carcinoma and type 2 diabetes mellitus. Type 2 diabetes is an established risk factor and high red meat intake has been associated with pancreatic ductal adenocarcinoma (PDAC). We hypothesize that genetic susceptibility to hemochromatosis as determined by known hereditary hemochromatosis-related genes will be associated with PDAC.

Methods

We conducted a pathway analysis of genes known to contribute to hereditary hemochromatosis using the summary-based adaptive rank truncated product (sARTP) method on GWAS summary statistics derived from 9038 PDAC cases and 12,389 controls of European descent collected by the Pancreatic Cancer Cohort Consortium (PanScan) and the Pancreatic Cancer Case-Control Consortium (PANC4).Our analysis included 7 hereditary hemochromatosis genes (HFE, BMP2, HJV, HAMP, TFR2, SLC40A1, and FTH1) and close genomic regions (20 kb upstream and 20 kb downstream) with a total of 176 single nucleotide polymorphisms (SNPs). The sARTP method combines SNP-level associations across SNPs in a gene or a pathway.

Results

The hereditary hemochromatosis pathway was significantly associated with PDAC (P-value = 0.011). HJV and TFR2 genes contributed the most to the association with PDAC risk (gene level P-values = 0.003 and 0.013, respectively).

Conclusions

This study supports the hypothesis that genetic susceptibility related to hereditary hemochromatosis genes are associated with PDAC. Further studies should evaluate the modifying effect of iron-rich foods and genetic susceptibility of hemochromatosis and PDAC risk.

Funding Sources

This work was supported by the Intramural Research Program, Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health.

CLPTM1L/CRR9 ectodomain interaction with GRP78 at the cell surface signals for survival and chemoresistance upon ER stress in pancreatic adenocarcinoma cells

Altered regulation of endoplasmic reticulum (ER) homeostasis has been implicated in many cancers and has recently become a therapeutic and chemosensitization target of interest. We have identified Cleft Lip and Palate Transmembrane 1-Like (CLPTM1L)/Cisplatin Resistance Related Protein 9 (CRR9) as an ER stress related mediator of cytoprotection in pancreatic cancer. We recently demonstrated that CLPTM1L is highly expressed in pancreatic ductal adenocarcinoma and associated with poor outcome. Furthermore, we have discovered that CLPTM1L interacts with phosphoinositol-3-kinase-alpha at the tumor cell surface and causes up-regulation of Bcl-xL and pAkt mediated survival signaling. Here, we demonstrate surface relocalization and survival signaling by CLPTM1L triggered by endoplasmic reticular (ER) stress. We demonstrate the interaction of CLPTM1L with the central ER stress survival mediator, Glucose Regulated Protein 78 (GRP78)/Binding Immunoglobulin Protein (BiP) and PI3K-alpha /p110α. This interaction and surface relocalization of CLPTM1L and GRP78 is induced by ER stress, including that caused by treatment with gemcitabine. We demonstrate that the extracellular loop of CLPTM1L is required for gemcitabine resistance and interaction with GRP78. This interaction and the chemoresistance effect conferred by this pathway is targetable with our recently developed inhibitory CLPTM1L antibodies, which may represent novel modalities of chemosensitization and treatment of pancreatic adenocarcinoma. Anchorage independent growth, GRP78-mediated chemoresistance, and Akt phosphorylation were abrogated by inhibition of CLPTM1L. These findings demonstrate a novel and potentially targetable mechanism of cytoprotection and chemoresistance in pancreatic tumors.

Association Between Telomere Length and Risk of Cancer and Non-Neoplastic Diseases: A Mendelian Randomization Study

IMPORTANCE

The causal direction and magnitude of the association between telomere length and incidence of cancer and non-neoplastic diseases is uncertain owing to the susceptibility of observational studies to confounding and reverse causation.

OBJECTIVE

To conduct a Mendelian randomization study, using germline genetic variants as instrumental variables, to appraise the causal relevance of telomere length for risk of cancer and non-neoplastic diseases.

DATA SOURCES

Genomewide association studies (GWAS) published up to January 15, 2015.

STUDY SELECTION

GWAS of noncommunicable diseases that assayed germline genetic variation and did not select cohort or control participants on the basis of preexisting diseases. Of 163 GWAS of noncommunicable diseases identified, summary data from 103 were available.

DATA EXTRACTION AND SYNTHESIS

Summary association statistics for single nucleotide polymorphisms (SNPs) that are strongly associated with telomere length in the general population.

MAIN OUTCOMES AND MEASURES

Odds ratios (ORs) and 95% confidence intervals (CIs) for disease per standard deviation (SD) higher telomere length due to germline genetic variation.

RESULTS

Summary data were available for 35 cancers and 48 non-neoplastic diseases, corresponding to 420 081 cases (median cases, 2526 per disease) and 1 093 105 controls (median, 6789 per disease). Increased telomere length due to germline genetic variation was generally associated with increased risk for site-specific cancers. The strongest associations (ORs [95% CIs] per 1-SD change in genetically increased telomere length) were observed for glioma, 5.27 (3.15-8.81); serous low-malignant-potential ovarian cancer, 4.35 (2.39-7.94); lung adenocarcinoma, 3.19 (2.40-4.22); neuroblastoma, 2.98 (1.92-4.62); bladder cancer, 2.19 (1.32-3.66); melanoma, 1.87 (1.55-2.26); testicular cancer, 1.76 (1.02-3.04); kidney cancer, 1.55 (1.08-2.23); and endometrial cancer, 1.31 (1.07-1.61). Associations were stronger for rarer cancers and at tissue sites with lower rates of stem cell division. There was generally little evidence of association between genetically increased telomere length and risk of psychiatric, autoimmune, inflammatory, diabetic, and other non-neoplastic diseases, except for coronary heart disease (OR, 0.78 [95% CI, 0.67-0.90]), abdominal aortic aneurysm (OR, 0.63 [95% CI, 0.49-0.81]), celiac disease (OR, 0.42 [95% CI, 0.28-0.61]) and interstitial lung disease (OR, 0.09 [95% CI, 0.05-0.15]).

CONCLUSIONS AND RELEVANCE

It is likely that longer telomeres increase risk for several cancers but reduce risk for some non-neoplastic diseases, including cardiovascular diseases.

Abstract 4604: DNA methylation at chromosome 8q24 in peripheral blood and prostate cancer risk

Chromosome 8q24 has emerged as an important region for the genetic susceptibility to several cancers, including prostate cancer; however, little is known about the contribution of DNA methylation in this region to risk. To evaluate DNA methylation levels at 8q24 in relation to the risk of prostate cancer, we conducted a nested case-control study using pre-diagnostic peripheral blood samples from 694 prostate cancer cases and 703 controls in the Prostate, Lung, Colorectal, and Ovarian (PLCO) Cancer Screening Trial. We quantified DNA methylation levels at 63 specific CpG sites nearby cancer susceptibility single nucleotide polymorphisms (SNPs) at 8q24 or the MYC oncogene using pyrosequencing with bisulfite-treated DNA. For each of the 50 CpG sites meeting our quality control requirements, we used logistic regression to estimate odds ratios (OR) and 95% confidence intervals (CI) for the risk of prostate cancer associated with each unit increase in percent DNA methylation. We identified 7 CpG sites whose DNA methylation levels were associated with prostate cancer risk (p less than 0.05), including 5 CpG sites in POU5F1B that were moderately correlated with one another (Spearman rho: 0.20-0.52), as well as 2 intergenic CpG sites (Chr8:128081593 and Chr8:128513944). The most significant association was observed for Chr8:128498079 in POU5F1B (OR = 0.97, 95% CI: 0.95-0.99, p-value = 0.01). The five POU5F1B CpG sites found to be associated with prostate cancer risk were significantly correlated with several established cancer susceptibility SNPs at 8q24 from genome-wide association studies, including rs6983267 (rho: 0.14-0.36), which has been previously associated with prostate cancer, among other cancers. Rs6983267 was only marginally associated with prostate cancer risk in this study (p = 0.10). When rs6983267 was included in the models for the five correlated CpG sites at POU5F1B, the associations for these CpG sites with prostate cancer were attenuated, suggesting that the POU5F1B CpG sites may be part of the same biological pathway conferring increased prostate cancer risk. Our findings suggest that DNA methylation at 8q24 may contribute to the risk of prostate cancer, and epigenetic mechanisms may help explain some of the SNP-associated risk in this important cancer susceptibility region.

Citation Format: Kathryn Hughes Barry, Lee E. Moore, Joshua Sampson, Stella Koutros, Liying Yan, Ann Meyer, Mahitha Reddy, Michael B. Cook, Joseph F. Fraumeni, Meredith Yeager, Laufey Amundadottir, Sonja I. Berndt. DNA methylation at chromosome 8q24 in peripheral blood and prostate cancer risk. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4604. doi:10.1158/1538-7445.AM2015-4604

Two susceptibility loci identified for prostate cancer aggressiveness

Most men diagnosed with prostate cancer will experience indolent disease; hence, discovering genetic variants that distinguish aggressive from nonaggressive prostate cancer is of critical clinical importance for disease prevention and treatment. In a multistage, case-only genome-wide association study of 12,518 prostate cancer cases, we identify two loci associated with Gleason score, a pathological measure of disease aggressiveness: rs35148638 at 5q14.3 (RASA1, P=6.49 × 10(-9)) and rs78943174 at 3q26.31 (NAALADL2, P=4.18 × 10(-8)). In a stratified case-control analysis, the SNP at 5q14.3 appears specific for aggressive prostate cancer (P=8.85 × 10(-5)) with no association for nonaggressive prostate cancer compared with controls (P=0.57). The proximity of these loci to genes involved in vascular disease suggests potential biological mechanisms worthy of further investigation.

DNA methylation levels at chromosome 8q24 in peripheral blood are associated with 8q24 cancer susceptibility loci

Chromosome 8q24 has emerged as an important region for genetic susceptibility to various cancers, but little is known about the contribution of DNA methylation at 8q24. To evaluate variability in DNA methylation levels at 8q24 and the relationship with cancer susceptibility single nucleotide polymorphisms (SNPs) in this region, we quantified DNA methylation levels in peripheral blood at 145 CpG sites nearby 8q24 cancer susceptibility SNPs or MYC using pyrosequencing among 80 Caucasian men in the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial. For the 60 CpG sites meeting quality control, which also demonstrated temporal stability over a 5-year period, we calculated pairwise Spearman correlations for DNA methylation levels at each CpG site with 42 8q24 cancer susceptibility SNPs. To account for multiple testing, we adjusted P values into q values reflecting the false discovery rate (FDR). In contrast to the MYC CpG sites, most sites nearby the SNPs demonstrated good reproducibility, high methylation levels, and moderate-high between-individual variation. We observed 10 statistically significant (FDR < 0.05) CpG site-SNP correlations. These included correlations between an intergenic CpG site at Chr8:128393157 and the prostate cancer SNP rs16902094 (ρ = -0.54; P = 9.7 × 10(-7); q = 0.002), a PRNCR1 CpG site at Chr8:128167809 and the prostate cancer SNP rs1456315 (ρ = 0.52; P = 1.4 × 10(-6); q = 0.002), and two POU5F1B CpG sites and several prostate/colorectal cancer SNPs (for Chr8:128498051 and rs6983267, ρ = 0.46; P = 2.0 × 10(-5); q = 0.01). This is the first report of correlations between blood DNA methylation levels and cancer susceptibility SNPs at 8q24, suggesting that DNA methylation at this important susceptibility locus may contribute to cancer risk.

Abstract 288: DNA methylation patterns in peripheral blood and the relationship with cancer susceptibility loci at chromosome 8q24

Chromosome 8q24 has emerged as an important region for genetic susceptibility to several cancers, but little is known about the contribution of DNA methylation in this region. To explore the extent of variation in DNA methylation at 8q24 in peripheral blood and its relationship with genetic variation in the region, we conducted a cross-sectional study using blood samples from 80 non-Hispanic Caucasian males in the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial. Specifically, we aimed to evaluate between-individual variation in DNA methylation levels at specific CpG sites at 8q24 and to investigate the underlying genetic structure in the region by examining correlations for DNA methylation levels with each other and with the established cancer susceptibility SNPs at 8q24. We quantified DNA methylation levels at 145 CpG sites nearby cancer susceptibility single nucleotide polymorphisms (SNPs) at 8q24 or the MYC oncogene using pyrosequencing of bisulfite-treated DNA, which is considered a highly sensitive method to detect differences in DNA methylation levels between individuals. We calculated pairwise Spearman correlations (rho), adjusting for multiple testing using the False Discovery Rate (FDR) method. We identified a large number of CpG sites that were reproducible and demonstrated moderate to high between-individual variation in our study population. Among these CpG sites, some sites within or nearby MYC and POU5F1B were strongly correlated with one another (highest rho=0.74), suggesting a coordination of DNA methylation levels in gene regions. Moreover, we observed strong correlations between several CpG sites and some of the known cancer susceptibility SNPs at 8q24. Some of the correlations remained statistically significant after adjustment for multiple comparisons, including a CpG site and an established prostate cancer SNP in the long non-coding RNA PRNCR1 (Chr8:128167809 and rs1456315, rho=0.52; p-value=1.4x10-6; FDR-adjusted p-value=0.002) and a CpG site (Chr8:128498051) in POU5F1B and the known prostate/colorectal cancer SNP rs6983267 (rho=0.46; p-value=2.0x10-5; FDR-adjusted p-value=0.01). This is the first study to report associations between DNA methylation levels at 8q24 in peripheral blood and 8q24 cancer susceptibility SNPs, suggesting that DNA methylation at this important susceptibility locus may contribute to risk. Additional studies are needed to clarify the relationship between genetic and epigenetic variation at 8q24 and cancer risk.

Citation Format: Kathryn Hughes Barry, Lee Moore, Joshua Sampson, Liying Yan, Ann Meyer, Charles C. Chung, Meredith Yeager, Laufey Amundadottir, Sonja I. Berndt. DNA methylation patterns in peripheral blood and the relationship with cancer susceptibility loci at chromosome 8q24. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 288. doi:10.1158/1538-7445.AM2014-288

Abstract 287: Characterizing the genetic basis of methylome diversity in histologically normal human lung tissue

The genetic regulation of the human epigenome is not fully appreciated. Here we describe the effects of genetic variants on the DNA methylome in human lung based on methylation-quantitative trait loci (meQTL) analyses. We report 34,304 cis- and 585 trans-meQTLs, a genetic-epigenetic interaction of surprising magnitude, including a regulatory hotspot. These findings are replicated in both breast and kidney tissues and show distinct patterns: cis-meQTLs mostly localize to CpG sites outside of genes, promoters, and CpG islands (CGIs) while trans-meQTLs are over-represented in promoter CGIs. meQTL SNPs are enriched in CTCF binding sites, DNaseI hypersensitivity regions and histone marks. Importantly, 4 of the 5 established lung cancer risk loci in European ancestry are cis-meQTLs and, in aggregate, cis-meQTLs are enriched for lung cancer risk in a genome-wide analysis of 11,587 subjects. Thus, inherited genetic variation may affect lung carcinogenesis by regulating the human methylome.

Citation Format: Jianxin Shi, Crystal Marconett, Jubao Duan, Paula Hyland, Peng Li, Zhaoming Wang, William Wheeler, Mihaela Campan, Jing Huang, Weiyin Zhou, Tim Triche, Laufey Amundadottir, Amy Hutchinson, Po-Han Chen, Beiyun Zhou, Brian Chung, Pier Alberto Bertazzi, Andrew W. Bergen, Mathew Freedman, Diane Lee, Kim Siegmund, Andrew Clay Warner, Ben Berman, Angela C. Pesatori, Zea Borok, Dario Consonni, Nilanjan Chatterjee, Margaret Tucker, Neil Caporaso, Stephen J. Chanock, Ite A. Laird-Offringa, Maria Teresa Landi. Characterizing the genetic basis of methylome diversity in histologically normal human lung tissue. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 287. doi:10.1158/1538-7445.AM2014-287

Genetic polymorphisms in the 9p21 region associated with risk of multiple cancers

The chromosome 9p21 region has been implicated in the pathogenesis of multiple cancers. We analyzed 9p21 single nucleotide polymorphisms (SNPs) from eight genome-wide association studies (GWAS) with data deposited in dbGaP, including studies of esophageal squamous cell carcinoma (ESCC), gastric cancer (GC), pancreatic cancer, renal cell carcinoma (RCC), lung cancer (LC), breast cancer (BrC), bladder cancer (BC) and prostate cancer (PrC). The number of subjects ranged from 2252 (PrC) to 7619 (LC). SNP-level analyses for each cancer were conducted by logistic regression or random-effects meta-analysis. A subset-based statistical approach (ASSET) was performed to combine SNP-level P values across multiple cancers. We calculated gene-level P values using the adaptive rank truncated product method. We identified that rs1063192 and rs2157719 in the CDKN2A/2B region were significantly associated with ESCC and rs2764736 (3' of TUSC1) was associated with BC (P ≤ 2.59 × 10(-6)). ASSET analyses identified four SNPs significantly associated with multiple cancers: rs3731239 (CDKN2A intronic) with ESCC, GC and BC (P = 3.96 × 10(-) (4)); rs10811474 (3' of IFNW1) with RCC and BrC (P = 0.001); rs12683422 (LINGO2 intronic) with RCC and BC (P = 5.93 × 10(-) (4)) and rs10511729 (3' of ELAVL2) with LC and BrC (P = 8.63 × 10(-) (4)). At gene level, CDKN2B, CDKN2A and CDKN2B-AS1 were significantly associated with ESCC (P ≤ 4.70 × 10(-) (5)). Rs10511729 and rs10811474 were associated with cis-expression of 9p21 genes in corresponding cancer tissues in the expression quantitative trait loci analysis. In conclusion, we identified several genetic variants in the 9p21 region associated with the risk of multiple cancers, suggesting that this region may contribute to a shared susceptibility across different cancer types.

Genome-wide association study of survival in patients with pancreatic adenocarcinoma

BACKGROUND AND OBJECTIVE

Survival of patients with pancreatic adenocarcinoma is limited and few prognostic factors are known. We conducted a two-stage genome-wide association study (GWAS) to identify germline variants associated with survival in patients with pancreatic adenocarcinoma.

METHODS

We analysed overall survival in relation to single nucleotide polymorphisms (SNPs) among 1005 patients from two large GWAS datasets, PanScan I and ChinaPC. Cox proportional hazards regression was used in an additive genetic model with adjustment for age, sex, clinical stage and the top four principal components of population stratification. The first stage included 642 cases of European ancestry (PanScan), from which the top SNPs (p≤10(-5)) were advanced to a joint analysis with 363 additional patients from China (ChinaPC).

RESULTS

In the first stage of cases of European descent, the top-ranked loci were at chromosomes 11p15.4, 18p11.21 and 1p36.13, tagged by rs12362504 (p=1.63×10(-7)), rs981621 (p=1.65×10(-7)) and rs16861827 (p=3.75×10(-7)), respectively. 131 SNPs with p≤10(-5) were advanced to a joint analysis with cases from the ChinaPC study. In the joint analysis, the top-ranked SNP was rs10500715 (minor allele frequency, 0.37; p=1.72×10(-7)) on chromosome 11p15.4, which is intronic to the SET binding factor 2 (SBF2) gene. The HR (95% CI) for death was 0.74 (0.66 to 0.84) in PanScan I, 0.79 (0.65 to 0.97) in ChinaPC and 0.76 (0.68 to 0.84) in the joint analysis.

CONCLUSIONS

Germline genetic variation in the SBF2 locus was associated with overall survival in patients with pancreatic adenocarcinoma of European and Asian ancestry. This association should be investigated in additional large patient cohorts.

Integrating gene expression and epidemiological data for the discovery of genetic interactions associated with cancer risk

Dozens of common genetic variants associated with cancer risk have been identified through genome-wide association studies (GWASs). However, these variants only explain a modest fraction of the heritability of disease. The missing heritability has been attributed to several factors, among them the existence of genetic interactions (G × G). Systematic screens for G × G in model organisms have revealed their fundamental influence in complex phenotypes. In this scenario, G × G overlap significantly with other types of gene and/or protein relationships. Here, by integrating predicted G × G from GWAS data and complex- and context-defined gene coexpression profiles, we provide evidence for G × G associated with cancer risk. G × G predicted from a breast cancer GWAS dataset identified significant overlaps [relative enrichments (REs) of 8-36%, empirical P values < 0.05 to 10(-4)] with complex (non-linear) gene coexpression in breast tumors. The use of gene or protein data not specific for breast cancer did not reveal overlaps. According to the predicted G × G, experimental assays demonstrated functional interplay between lipoma-preferred partner and transforming growth factor-β signaling in the MCF10A non-tumorigenic mammary epithelial cell model. Next, integration of pancreatic tumor gene expression profiles with pancreatic cancer G × G predicted from a GWAS corroborated the observations made for breast cancer risk (REs of 25-59%). The method presented here can potentially support the identification of genetic interactions associated with cancer risk, providing novel mechanistic hypotheses for carcinogenesis.

Abstract 3137: Functional analysis of the chr13q22.1 pancreatic cancer risk locus suggests allele-specific effects on DIS3 expression

Pancreatic cancer is the 10th most common cancer and 4th most common cause of cancer mortality in the United States. A genome wide association study has revealed pancreatic cancer susceptibility loci on chromosomes 13q22.1, 1q32.1 and 5p15.33. The chr13q22.1 region contains the most significant of all risk SNPs, rs9543325, (P= 3.27 x 10−11) which is located in a gene desert. The nearest genes are KLF5, KLF12, PIBF1, DIS3 and BORA, which range in distance from 265kb to 586kb, respectively, from the most significant variants. Given the high linkage disequilibrium across this risk locus, there are many candidate functional variants to consider. Imputation in the region did not improve the signal but gave a set of highly correlated SNPs which likely includes the functional variant(s). In an effort to identify the functional variant(s) we performed eQTL analyses to test the association between the genotypes of these candidate SNPs and expression of nearby genes. Among 100 normal pancreatic tissue samples, DIS3 showed the strongest association with SNPs in our risk locus (P-values as low as 0.0004). Mutations in DIS3 have been identified in acute myeloid leukemia and multiple myeloma, and its expression has been correlated with metastatic potential in colorectal cancer, suggesting this gene could be important in pancreatic cancer biology. Chromosome Conformation Capture (3C) was performed to test for physical interactions between the risk locus and nearby genes. This assay confirmed the three dimensional proximity of the risk locus and DIS3 promoter. Sub-regions of the risk locus were then cloned upstream of a minimal promoter controlling luciferase expression to assay for potential allele specific enhancer/silencer activity. This assay revealed a sub-region downstream of rs9543325 that causes allele-specific silencing. These results suggest that a sub-region of the chr13q22.1 risk locus has allele-specific effects on DIS3 expression, which may have implications in the susceptibility to pancreatic cancer.

Citation Format: Jason Hoskins, Abdisamad Ibrahim, Jinping Jia, Irene Collins, Hemang Parikh, Gloria M. Petersen, Laufey Amundadottir. Functional analysis of the chr13q22.1 pancreatic cancer risk locus suggests allele-specific effects on DIS3 expression. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 3137. doi:10.1158/1538-7445.AM2013-3137

Abstract 2556: Functional characterization of the pancreatic cancer TERT-CLPTM1L risk locus on chr5p15.33

Pancreatic cancer is a highly lethal cancer with few well established risk factors. PanScan, a genome wide association study (GWAS) of pancreatic cancer has identified four pancreatic cancer susceptibility loci in populations of European ancestry. One of these is in a multi cancer locus on chr5p15.33 where the most significant SNP from the GWAS (rs401681, P=3.7x10-7, ORAllele=1.19) lies in a region containing two genes, TERT and CLPTM1L. The TERT gene encodes the catalytic subunit of telomerase, well known for its essential role in maintaining telomere ends. The function of CLPTM1L is not as clear, although it has been proposed to play a role in apoptosis. It is predicted to encode a protein with 6 transmembrane (TM) domains and two large hydrophilic domains: a loop of 253 aa between the first and second TM domains, and a C-terminal tail of 89 aa.

We have performed imputation to fine-map the signal to a SNP three orders of magnitude more significant than the GWAS SNP (PImputed=1.4x10-10, ORAllele=1.30). As this SNP is located in the CLPTM1L gene we have performed a series of experiments to investigate the function of the CLPTM1L gene and its encoded protein. Immunofluorescence analysis in pancreatic cancer cells (PANC-1) indicates that it localizes to the endoplasmic reticulum. Affinity purification and mass spectrometry (HEK-293T, hTERT-HPNE and PANC-1 cells) identified MYH9, a non-muscle heavy chain myosin, as a potential interacting protein. The interaction has been validated by co-immunoprecipitation and co-localization experiments. To examine if CLPTM1L plays a role in growth control, we created stable PANC-1 cell lines overexpressing the full length CLPTM1L gene as well as two deletions, a C-terminal deletion and a loop deletion, and assayed growth in vitro and in vivo. Cell lines overexpressing full length CLPTM1L grow faster in vitro and in vivo as compared to cells containing empty vector. Interestingly, the two CLPTM1L mutants abolish this effect. Furthermore, we have shown by RNA-seq that the CLPTM1L gene is overexpressed in pancreatic tumors as compared to normal pancreatic tissues. Our results indicate that CLPTM1L may play a role in the control of cell growth and oncogenesis in the pancreas. Our current efforts aim at further characterizing the function of CLPTM1L and to correlate pancreatic cancer risk variants on 5p15.33 to molecular phenotypes to attempt to explain the underlying biology of the risk.

Citation Format: Jinping Jia, Irene Collins, Marta Dzyadyk, Abbey Thompson, Adam Cheuk, Hemang Parikh, Zhaoming Wang, Chris Westlake, Allen Bosley, Gloria Petersen, Thorkell Andresson, Laufey Amundadottir. Functional characterization of the pancreatic cancer TERT-CLPTM1L risk locus on chr5p15.33. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 2556. doi:10.1158/1538-7445.AM2013-2556

Note: This abstract was not presented at the AACR Annual Meeting 2013 because the presenter was unable to attend.

Characterization of SNPs associated with prostate cancer in men of Ashkenazic descent from the set of GWAS identified SNPs: impact of cancer family history and cumulative SNP risk prediction

BACKGROUND

Genome-wide association studies (GWAS) have identified multiple SNPs associated with prostate cancer (PrCa). Population isolates may have different sets of risk alleles for PrCa constituting unique population and individual risk profiles.

METHODS

To test this hypothesis, associations between 31 GWAS SNPs of PrCa were examined among 979 PrCa cases and 1,251 controls of Ashkenazic descent using logistic regression. We also investigated risks by age at diagnosis, pathological features of PrCa, and family history of cancer. Moreover, we examined associations between cumulative number of risk alleles and PrCa and assessed the utility of risk alleles in PrCa risk prediction by comparing the area under the curve (AUC) for different logistic models.

RESULTS

Of the 31 genotyped SNPs, 8 were associated with PrCa at p ≤ 0.002 (corrected p-value threshold) with odds ratios (ORs) ranging from 1.22 to 1.42 per risk allele. Four SNPs were associated with aggressive PrCa, while three other SNPs showed potential interactions for PrCa by family history of PrCa (rs8102476; 19q13), lung cancer (rs17021918; 4q22), and breast cancer (rs10896449; 11q13). Men in the highest vs. lowest quartile of cumulative number of risk alleles had ORs of 3.70 (95% CI 2.76-4.97); 3.76 (95% CI 2.57-5.50), and 5.20 (95% CI 2.94-9.19) for overall PrCa, aggressive cancer and younger age at diagnosis, respectively. The addition of cumulative risk alleles to the model containing age at diagnosis and family history of PrCa yielded a slightly higher AUC (0.69 vs. 0.64).

CONCLUSION

These data define a set of risk alleles associated with PrCa in men of Ashkenazic descent and indicate possible genetic differences for PrCa between populations of European and Ashkenazic ancestry. Use of genetic markers might provide an opportunity to identify men at highest risk for younger age of onset PrCa; however, their clinical utility in identifying men at highest risk for aggressive cancer remains limited.

Polymorphisms in genes related to one-carbon metabolism are not related to pancreatic cancer in PanScan and PanC4

PURPOSE

The evidence of a relation between folate intake and one-carbon metabolism (OCM) with pancreatic cancer (PanCa) is inconsistent. In this study, the association between genes and single-nucleotide polymorphisms (SNPs) related to OCM and PanCa was assessed.

METHODS

Using biochemical knowledge of the OCM pathway, we identified thirty-seven genes and 834 SNPs to examine in association with PanCa. Our study included 1,408 cases and 1,463 controls nested within twelve cohorts (PanScan). The ten SNPs and five genes with lowest p values (<0.02) were followed up in 2,323 cases and 2,340 controls from eight case-control studies (PanC4) that participated in PanScan2. The correlation of SNPs with metabolite levels was assessed for 649 controls from the European Prospective Investigation into Cancer and Nutrition.

RESULTS

When both stages were combined, we observed suggestive associations with PanCa for rs10887710 (MAT1A) (OR 1.13, 95 %CI 1.04-1.23), rs1552462 (SYT9) (OR 1.27, 95 %CI 1.02-1.59), and rs7074891 (CUBN) (OR 1.91, 95 %CI 1.12-3.26). After correcting for multiple comparisons, no significant associations were observed in either the first or second stage. The three suggested SNPs showed no correlations with one-carbon biomarkers.

CONCLUSIONS

This is the largest genetic study to date to examine the relation between germline variations in OCM-related genes polymorphisms and the risk of PanCa. Suggestive evidence for an association between polymorphisms and PanCa was observed among the cohort-nested studies, but this did not replicate in the case-control studies. Our results do not strongly support the hypothesis that genes related to OCM play a role in pancreatic carcinogenesis.

Diabetes and risk of pancreatic cancer: a pooled analysis from the pancreatic cancer cohort consortium

PURPOSE

Diabetes is a suspected risk factor for pancreatic cancer, but questions remain about whether it is a risk factor or a result of the disease. This study prospectively examined the association between diabetes and the risk of pancreatic adenocarcinoma in pooled data from the NCI pancreatic cancer cohort consortium (PanScan).

METHODS

The pooled data included 1,621 pancreatic adenocarcinoma cases and 1,719 matched controls from twelve cohorts using a nested case-control study design. Subjects who were diagnosed with diabetes near the time (<2 years) of pancreatic cancer diagnosis were excluded from all analyses. All analyses were adjusted for age, race, gender, study, alcohol use, smoking, BMI, and family history of pancreatic cancer.

RESULTS

Self-reported diabetes was associated with a forty percent increased risk of pancreatic cancer (OR = 1.40, 95 % CI: 1.07, 1.84). The association differed by duration of diabetes; risk was highest for those with a duration of 2-8 years (OR = 1.79, 95 % CI: 1.25, 2.55); there was no association for those with 9+ years of diabetes (OR = 1.02, 95 % CI: 0.68, 1.52).

CONCLUSIONS

These findings provide support for a relationship between diabetes and pancreatic cancer risk. The absence of association in those with the longest duration of diabetes may reflect hypoinsulinemia and warrants further investigation.

Abstract 28: Characterization of SNPs associated with prostate cancer in men of Ashkenazic descent from the set of GWAS identified SNPs: Impact of cancer family history

Background: Genome-wide association studies (GWAS) have identified several SNPs that are independently associated with prostate cancer (PrCa) in Caucasian men. Population isolates with unique allele frequencies are likely to have different sets of risk alleles for PrCa constituting personalized population risk stratification in addition to individual risk profiles.

Methods: To test this hypothesis, we examined associations between 31 GWAS SNPs and PrCa among 979 cases and 1,251 controls of Ashkenazic descent. Odds ratios (OR) and 95% confidence intervals (CIs) were estimated for each SNP using logistic regression models adjusted for age. We also investigated risk in strata by age at diagnosis, pathological features of PrCa, as well as family history of PrCa and other cancers.

Results: Of 31 SNPs, 19 (62%) were associated with PrCa at p&lt;0.05 in one or more genetic risk models. The strongest associations with PrCa were for rs6983267 at 8q24 (OR=1.78, p=5.7 x 10−7;comparison of GG vs. TT genotype), and rs1465618 at 2p21 (OR=2.19; p=0.003; AA vs. AG/GG). Six SNPs were associated with aggressive PrCa, while six others with less aggressive cancer. Four SNPs showed statistically significant interactions between PrCa and family history of PrCa (rs8102476 at 19q13), lung cancer (rs17021918 at 4q22), colorectal cancer (rs6983267 at 8q24), and breast cancer (rs10896449 at 11q13).

Conclusion: These data define a specific set of risk alleles associated with PrCa in Ashkenazic men. Results demonstrate differences between men of European and Ashkenazic ancestry in relation to genetic susceptibility of PrCa. Interactions of genetic variants with family history of specific cancers indicate a complex network of inherited cancer risk syndromes still to be defined. Use of genetic markers to identify men at higher risk of PrCa will provide an opportunity for personalized care of PrCa prevention through screening of high-risk individuals, independent of family history information.

Citation Format: Ilir Agalliu, Zhaoming Wang, Anne Dunn, Hemang Parikh, Timothy Myers, Robert Burk, Laufey Amundadottir. Characterization of SNPs associated with prostate cancer in men of Ashkenazic descent from the set of GWAS identified SNPs: Impact of cancer family history. [abstract]. In: Proceedings of the AACR Special Conference on Post-GWAS Horizons in Molecular Epidemiology: Digging Deeper into the Environment; 2012 Nov 11-14; Hollywood, FL. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2012;21(11 Suppl):Abstract nr 28.

Pathway analysis of genome-wide association study data highlights pancreatic development genes as susceptibility factors for pancreatic cancer

Four loci have been associated with pancreatic cancer through genome-wide association studies (GWAS). Pathway-based analysis of GWAS data is a complementary approach to identify groups of genes or biological pathways enriched with disease-associated single-nucleotide polymorphisms (SNPs) whose individual effect sizes may be too small to be detected by standard single-locus methods. We used the adaptive rank truncated product method in a pathway-based analysis of GWAS data from 3851 pancreatic cancer cases and 3934 control participants pooled from 12 cohort studies and 8 case-control studies (PanScan). We compiled 23 biological pathways hypothesized to be relevant to pancreatic cancer and observed a nominal association between pancreatic cancer and five pathways (P < 0.05), i.e. pancreatic development, Helicobacter pylori lacto/neolacto, hedgehog, Th1/Th2 immune response and apoptosis (P = 2.0 × 10(-6), 1.6 × 10(-5), 0.0019, 0.019 and 0.023, respectively). After excluding previously identified genes from the original GWAS in three pathways (NR5A2, ABO and SHH), the pancreatic development pathway remained significant (P = 8.3 × 10(-5)), whereas the others did not. The most significant genes (P < 0.01) in the five pathways were NR5A2, HNF1A, HNF4G and PDX1 for pancreatic development; ABO for H.pylori lacto/neolacto; SHH for hedgehog; TGFBR2 and CCL18 for Th1/Th2 immune response and MAPK8 and BCL2L11 for apoptosis. Our results provide a link between inherited variation in genes important for pancreatic development and cancer and show that pathway-based approaches to analysis of GWAS data can yield important insights into the collective role of genetic risk variants in cancer.

Y chromosome haplogroups and prostate cancer in populations of European and Ashkenazi Jewish ancestry

Genetic variation on the Y chromosome has not been convincingly implicated in prostate cancer risk. To comprehensively analyze the role of inherited Y chromosome variation in prostate cancer risk in individuals of European ancestry, we genotyped 34 binary Y chromosome markers in 3,995 prostate cancer cases and 3,815 control subjects drawn from four studies. In this set, we identified nominally significant association between a rare haplogroup, E1b1b1c, and prostate cancer in stage I (P = 0.012, OR = 0.51; 95% confidence interval 0.30–0.87). Population substructure of E1b1b1c carriers suggested Ashkenazi Jewish ancestry, prompting a replication phase in individuals of both European and Ashkenazi Jewish ancestry. The association was not significant for prostate cancer overall in studies of either Ashkenazi Jewish (1,686 cases and 1,597 control subjects) or European (686 cases and 734 control subjects) ancestry (P_meta = 0.078), but a meta-analysis of stage I and II studies revealed a nominally significant association with prostate cancer risk (P_meta = 0.010, OR = 0.77; 95% confidence interval 0.62–0.94). Comparing haplogroup frequencies between studies, we noted strong similarities between those conducted in the US and France, in which the majority of men carried R1 haplogroups, resembling Northwestern European populations. On the other hand, Finns had a remarkably different haplogroup distribution with a preponderance of N1c and I1 haplogroups. In summary, our results suggest that inherited Y chromosome variation plays a limited role in prostate cancer etiology in European populations but warrant follow-up in additional large and well characterized studies of multiple ethnic backgrounds.

Abstract 2778: Fine mapping the KLK3 locus on chromosome 19q13.33 associated with prostate cancer susceptibility and PSA levels

Prostate cancer is the most commonly diagnosed non-cutaneous cancer in men in the United States with over 200,000 new cases and 30,000 deaths estimated in 2010. Measurements of serum prostate specific antigen (PSA) protein levels form the basis for a widely used test to screen men for prostate cancer. Germline variants in the gene that encodes the PSA protein (KLK3) have been shown to be associated with serum PSA levels. In genome-wide association studies (GWAS) of prostate cancer, KLK3 variants have been suggested to be prostate cancer risk factors. Based on a resequence analysis of a 56 kb region on chromosome 19q13.33, centered on the KLK3 gene, that characterized common germline variations in the region, we fine mapped this locus by genotyping tag SNPs in 3,522 prostate cancer cases and 3,338 controls from five cohort and case-control studies. The cohort studies are: the Prostate, Lung, Colon and Ovarian (PLCO) Cancer Screening Trial, the Alpha-Tocopherol, Beta-Carotene Cancer Prevention Study (ATBC), the American Cancer Society Prevention Study II Nutrition Cohort (CPS-II) and the Cohort of Norway (CONOR). The case-control study is the French Prostate Case-Control Study (CeRePP) which is a hospital based case-control study. This study has &gt; 80% power to detect an association with an odds ratio of 1.25 (assuming a MAF of 0.05, a prevalence of prostate cancer = 1.5067% and alpha = 0.05). We did not observe a strong association for KLK3 variants, previously reported to confer risk for prostate cancer (rs2735839; P = 0.20). However, three highly correlated SNPs (rs17632542, rs62113212 and rs62113214) demonstrated an association with prostate cancer (P = 3.41×10−4, per-allele odds ratio (OR) = 0.77, 95% CI = 0.67-0.89). In particular, the signal is apparent in nonaggressive prostate cancer cases with Gleason score &lt; 7 and disease stage &lt; III (P = 2.62×10−5, per-allele trend OR = 0.67, 95% CI = 0.56-0.81) but not in advanced cases with Gleason score &gt; 8 or stage ≥ III (P = 0.31, per-allele OR = 1.12, 95% CI = 0.90-1.40). One of the three highly correlated SNPs, rs17632542, introduces a non-synonymous amino acid change in the KLK3 protein with a predicted benign or neutral functional impact. Baseline PSA levels were lower in control subjects with one or more minor alleles at any one of the three SNPs (1.12 ng/ml, 95% CI = 0.96-1.28) as compared to those with no minor alleles (1.61 ng/ml, 95% CI = 1.49-1.72) (P = 9.70×10−5). Together our results suggest that germline KLK3 variants could influence the diagnosis of nonaggressive prostate cancer by influencing the likelihood of biopsy. It is possible that the KLK3 locus also contributes to prostate cancer risk, but additional studies will be needed to dissect the contribution of KLK3 to PSA and prostate risk separately.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 2778. doi:10.1158/1538-7445.AM2011-2778

Abstract 1170: MicroRNA gene polymorphisms in cancer susceptibility loci 5p12, 8q24 and 11q13

Using Genome Wide Association Studies (GWAS), we and others have previously identified and replicated genetic variations in chromosome 5p12, 8q24 and 11q13. Interestingly, the majority of risk-associated variations map to non-coding regions of the genome and thus, the molecular mechanism of single nucleotide polymorphisms (SNPs) underlying increased cancer risk remains unknown. To assess whether these SNPs could affect the biogenesis and function of non-coding RNA species, such as microRNAs (miRNA), here we assessed miRNA gene polymorphisms in these cancer susceptibility loci. We catalogued two miRNA genes in 5p12, sixteen miRNA genes in 8q24, and seven miRNA genes in 11q13. We re-sequenced the pri- and pre-miRNA genomic sequences in these loci using HapMap DNA samples. As a result, in 8q24 miRNA genes, we identified three new SNPs and confirmed two previously reported SNPs. In 11q13, we identified four genetic variations within different miRNA genes. However, linkage disequilibrium plot analysis of both loci didn't show any statistically significant association of any of these SNPs with their corresponding GWAS signals. Thus, we consider it unlikely that miRNA in these regions are direct downstream effectors of cancer-associated genetic variations. Nevertheless, because these risk loci are prone to genomic instability, we assessed whether the SNPs could affect miRNA gene expression since indirect associations with increased cancer susceptibility can't be excluded. To address this question, we cloned pre-miRNA genes from these loci, with or without genetic variations, into miRNA expression vectors and determined their relative expression in transfected 293T cells. We tested for an 8q24-encoded miRNA (miR-1206) or an 11q13 encoded-miRNA (miR-612) expression using miRNA TaqMan assays, and we found that the G allele polymorphism (rs2114358) in miR-1206 significantly downregulated miRNA expression levels compared to the reference A allele. On the other hand, an rs550894 (C/A) in miR-612 significantly upregulated miRNA expression levels compared to that of the reference C allele. These data indicate that the miRNA SNPs in these loci have indeed functional consequences that lead to dysregulation of miRNA biogenesis. Whether such changes in miRNA expression could also result in increased cancer risk needs to be assessed in future studies by identifying the downstream target genes of these miRNAs.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 1170. doi:10.1158/1538-7445.AM2011-1170

Large-scale exploration of gene-gene interactions in prostate cancer using a multistage genome-wide association study

Recent genome-wide association studies have identified independent susceptibility loci for prostate cancer that could influence risk through interaction with other, possibly undetected, susceptibility loci. We explored evidence of interaction between pairs of 13 known susceptibility loci and single nucleotide polymorphisms (SNP) across the genome to generate hypotheses about the functionality of prostate cancer susceptibility regions. We used data from Cancer Genetic Markers of Susceptibility: Stage I included 523,841 SNPs in 1,175 cases and 1,100 controls; Stage II included 27,383 SNPs in an additional 3,941 cases and 3,964 controls. Power calculations assessed the magnitude of interactions our study is likely to detect. Logistic regression was used with alternative methods that exploit constraints of gene-gene independence between unlinked loci to increase power. Our empirical evaluation demonstrated that an empirical Bayes (EB) technique is powerful and robust to possible violation of the independence assumption. Our EB analysis identified several noteworthy interacting SNP pairs, although none reached genome-wide significance. We highlight a Stage II interaction between the major prostate cancer susceptibility locus in the subregion of 8q24 that contains POU5F1B and an intronic SNP in the transcription factor EPAS1, which has potentially important functional implications for 8q24. Another noteworthy result involves interaction of a known prostate cancer susceptibility marker near the prostate protease genes KLK2 and KLK3 with an intronic SNP in PRXX2. Overall, the interactions we have identified merit follow-up study, particularly the EPAS1 interaction, which has implications not only in prostate cancer but also in other epithelial cancers that are associated with the 8q24 locus.

Variant ABO blood group alleles, secretor status, and risk of pancreatic cancer: results from the pancreatic cancer cohort consortium

BACKGROUND

Subjects with non-O ABO blood group alleles have increased risk of pancreatic cancer. Glycosyltransferase activity is greater for the A(1) versus A(2) variant, whereas O01 and O02 variants are nonfunctioning. We hypothesized: 1) A(1) allele would confer greater risk than A(2) allele, 2) protective effect of the O allele would be equivalent for O01 and O02 variants, 3) secretor phenotype would modify the association with risk.

METHODS

We determined ABO variants and secretor phenotype from single nucleotide polymorphisms in ABO and FUT2 genes in 1,533 cases and 1,582 controls from 12 prospective cohort studies. Adjusted odds ratios (OR) for pancreatic cancer were calculated using logistic regression.

RESULTS

An increased risk was observed in participants with A(1) but not A(2) alleles. Compared with subjects with genotype O/O, genotypes A(2)/O, A(2)/A(1), A(1)/O, and A(1)/A(1) had ORs of 0.96 (95% CI, 0.72-1.26), 1.46 (95% CI, 0.98-2.17), 1.48 (95% CI, 1.23-1.78), and 1.71 (95% CI, 1.18-2.47). Risk was similar for O01 and O02 variant O alleles. Compared with O01/O01, the ORs for each additional allele of O02, A(1), and A(2) were 1.00 (95% CI, 0.87-1.14), 1.38 (95% CI, 1.20-1.58), and 0.96 (95% CI, 0.77-1.20); P, O01 versus O02 = 0.94, A(1) versus A(2) = 0.004. Secretor phenotype was not an effect modifier (P-interaction = 0.63).

CONCLUSIONS

Among participants in a large prospective cohort consortium, ABO allele subtypes corresponding to increased glycosyltransferase activity were associated with increased pancreatic cancer risk.

IMPACT

These data support the hypothesis that ABO glycosyltransferase activity influences pancreatic cancer risk rather than actions of other nearby genes on chromosome 9q34.

Family history of cancer and risk of pancreatic cancer: a pooled analysis from the Pancreatic Cancer Cohort Consortium (PanScan)

A family history of pancreatic cancer has consistently been associated with increased risk of pancreatic cancer. However, uncertainty remains about the strength of this association. Results from previous studies suggest a family history of select cancers (i.e., ovarian, breast and colorectal) could also be associated, although not as strongly, with increased risk of pancreatic cancer. We examined the association between a family history of 5 types of cancer (pancreas, prostate, ovarian, breast and colorectal) and risk of pancreatic cancer using data from a collaborative nested case-control study conducted by the Pancreatic Cancer Cohort Consortium. Cases and controls were from cohort studies from the United States, Europe and China, and a case-control study from the Mayo Clinic. Analyses of family history of pancreatic cancer included 1,183 cases and 1,205 controls. A family history of pancreatic cancer in a parent, sibling or child was associated with increased risk of pancreatic cancer [multivariate-adjusted odds ratios (ORs) = 1.76, 95% confidence interval (CI) = 1.19-2.61]. A family history of prostate cancer was also associated with increased risk (OR = 1.45, 95% CI = 1.12-1.89). There were no statistically significant associations with a family history of ovarian cancer (OR = 0.82, 95% CI = 0.52-1.31), breast cancer (OR = 1.21, 95% CI = 0.97-1.51) or colorectal cancer (OR = 1.17, 95% CI = 0.93-1.47). Our results confirm a moderate sized association between a family history of pancreatic cancer and risk of pancreatic cancer and also provide evidence for an association with a family history of prostate cancer worth further study.

Anthropometric measures, body mass index, and pancreatic cancer: a pooled analysis from the Pancreatic Cancer Cohort Consortium (PanScan)

BACKGROUND

Obesity has been proposed as a risk factor for pancreatic cancer.

METHODS

Pooled data were analyzed from the National Cancer Institute Pancreatic Cancer Cohort Consortium (PanScan) to study the association between prediagnostic anthropometric measures and risk of pancreatic cancer. PanScan applied a nested case-control study design and included 2170 cases and 2209 control subjects. Odds ratios (ORs) and 95% confidence intervals (CIs) were estimated using unconditional logistic regression for cohort-specific quartiles of body mass index (BMI [calculated as weight in kilograms divided by height in meters squared]), weight, height, waist circumference, and waist to hip ratio as well as conventional BMI categories (underweight, <18.5; normal weight, 18.5-24.9; overweight, 25.0-29.9; obese, 30.0-34.9; and severely obese, > or = 35.0). Models were adjusted for potential confounders.

RESULTS

In all of the participants, a positive association between increasing BMI and risk of pancreatic cancer was observed (adjusted OR for the highest vs lowest BMI quartile, 1.33; 95% CI, 1.12-1.58; P(trend) < .001). In men, the adjusted OR for pancreatic cancer for the highest vs lowest quartile of BMI was 1.33 (95% CI, 1.04-1.69; P(trend) < .03), and in women it was 1.34 (95% CI, 1.05-1.70; P(trend) = .01). Increased waist to hip ratio was associated with increased risk of pancreatic cancer in women (adjusted OR for the highest vs lowest quartile, 1.87; 95% CI, 1.31-2.69; P(trend) = .003) but less so in men.

CONCLUSIONS

These findings provide strong support for a positive association between BMI and pancreatic cancer risk. In addition, centralized fat distribution may increase pancreatic cancer risk, especially in women.

Evaluation of association of HNF1B variants with diverse cancers: collaborative analysis of data from 19 genome-wide association studies

BACKGROUND

Genome-wide association studies have found type 2 diabetes-associated variants in the HNF1B gene to exhibit reciprocal associations with prostate cancer risk. We aimed to identify whether these variants may have an effect on cancer risk in general versus a specific effect on prostate cancer only.

METHODOLOGY/PRINCIPAL FINDINGS

In a collaborative analysis, we collected data from GWAS of cancer phenotypes for the frequently reported variants of HNF1B, rs4430796 and rs7501939, which are in linkage disequilibrium (r(2) = 0.76, HapMap CEU). Overall, the analysis included 16 datasets on rs4430796 with 19,640 cancer cases and 21,929 controls; and 21 datasets on rs7501939 with 26,923 cases and 49,085 controls. Malignancies other than prostate cancer included colorectal, breast, lung and pancreatic cancers, and melanoma. Meta-analysis showed large between-dataset heterogeneity that was driven by different effects in prostate cancer and other cancers. The per-T2D-risk-allele odds ratios (95% confidence intervals) for rs4430796 were 0.79 (0.76, 0.83)] per G allele for prostate cancer (p<10(-15) for both); and 1.03 (0.99, 1.07) for all other cancers. Similarly for rs7501939 the per-T2D-risk-allele odds ratios (95% confidence intervals) were 0.80 (0.77, 0.83) per T allele for prostate cancer (p<10(-15) for both); and 1.00 (0.97, 1.04) for all other cancers. No malignancy other than prostate cancer had a nominally statistically significant association.

CONCLUSIONS/SIGNIFICANCE

The examined HNF1B variants have a highly specific effect on prostate cancer risk with no apparent association with any of the other studied cancer types.

Abstract 164: High-throughput epigenetic analysis of susceptibility loci identified by GWAS in pancreatic cancer

Pancreatic cancer is a highly lethal cancer with few well established risk factors. A genome wide association study (GWAS) of pancreatic cancer (PanScan) is being conducted within the framework of the NCI-sponsored Cohort Consortium and the Pancreatic Cancer Case-Control Consortium (PANC4). Susceptibility loci discovered to date in PanScan are non-coding variants that lie in intronic or intergenic regions, suggesting that the underlying signals may function through regulatory mechanisms that influence gene expression or splicing. Alternatively, these may lie in unannotated transcripts and directly affect their function.

Epigenetic mechanisms, such as DNA methylation, can affect the regulation of gene expression and plays a critical role in the development of many human diseases including cancer. Powerful methods exist to analyze DNA methylation patterns in higher eukaryotes including methylated DNA immunoprecipitation (MeDIP), an affinity based approach to enrich methylated DNA regions from genomic DNA, which can be combined with microarrays to profile genomic DNA methylation patterns. We created a new semi-custom MeDIP-optimized array design based on our Human DNA Methylation 3×720K CpG Island Plus RefSeq Promoter array by tiling additional regions associated with pancreatic cancer that were identified in the PanScan study. To complement DNA methylation data a genome-wide transcriptome analysis was performed with RNA-sequencing (RNA-seq). Here we describe a comprehensive genome wide analysis of 8 pancreatic cancer cell lines to examine methylation patterns of Refseq promoters and annotated CpG islands as well as transcribed sequences.

One of the susceptibility loci from PanScan is in the vicinity of the ABO gene on Chr9q34 where four SNPs (rs505922, rs495828, rs657152 and rs630014) are associated with a significantly increased risk of pancreatic cancer. As a pilot study, DNA methylation patterns and expressed sequences in this locus were investigated in cell lines derived from pancreatic tumors and normal pancreatic tissues. Our genome wide DNA methylation and RNA-seq analysis aims at establishing a comprehensive catalog of epigenetic patterns in pancreatic cell lines that could provide plausibility for the association signal in the ABO gene and other GWAS regions and thereby, initiate the characterization of the molecular phenotype of the susceptibility alleles for pancreatic cancer risk.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 164.

Pancreatic cancer risk and ABO blood group alleles: results from the pancreatic cancer cohort consortium

A recent genome-wide association study (PanScan) identified significant associations at the ABO gene locus with risk of pancreatic cancer, but the influence of specific ABO genotypes remains unknown. We determined ABO genotypes (OO, AO, AA, AB, BO, and BB) in 1,534 cases and 1,583 controls from 12 prospective cohorts in PanScan, grouping participants by genotype-derived serologic blood type (O, A, AB, and B). Adjusted odds ratios (ORs) for pancreatic cancer by ABO alleles were calculated using logistic regression. Compared with blood type O, the ORs for pancreatic cancer in subjects with types A, AB, and B were 1.38 [95% confidence interval (95% CI), 1.18-1.62], 1.47 (95% CI, 1.07-2.02), and 1.53 (95% CI, 1.21-1.92), respectively. The incidence rates for blood types O, A, AB, and B were 28.9, 39.9, 41.8, and 44.5 cases per 100,000 subjects per year. An increase in risk was noted with the addition of each non-O allele. Compared with OO genotype, subjects with AO and AA genotype had ORs of 1.33 (95% CI, 1.13-1.58) and 1.61 (95% CI, 1.22-2.18), whereas subjects with BO and BB genotypes had ORs of 1.45 (95% CI, 1.14-1.85) and 2.42 (1.28-4.57). The population attributable fraction for non-O blood type was 19.5%. In a joint model with smoking, current smokers with non-O blood type had an adjusted OR of 2.68 (95% CI, 2.03-3.54) compared with nonsmokers of blood type O. We concluded that ABO genotypes were significantly associated with pancreatic cancer risk.

Identification of a new prostate cancer susceptibility locus on chromosome 8q24

We report a genome-wide association study in 10,286 cases and 9,135 controls of European ancestry in the Cancer Genetic Markers of Susceptibility (CGEMS) initiative. We identify a new association with prostate cancer risk on chromosome 8q24 (rs620861, P = 1.3 x 10(-10), heterozygote OR = 1.17, 95% CI 1.10-1.24; homozygote OR = 1.33, 95% CI 1.21-1.45). This defines a new locus associated with prostate cancer susceptibility on 8q24.

Cigarette smoking and pancreatic cancer: a pooled analysis from the pancreatic cancer cohort consortium

Smoking is an established risk factor for pancreatic cancer; however, detailed examination of the association of smoking intensity, smoking duration, and cumulative smoking dose with pancreatic cancer is limited. The authors analyzed pooled data from the international Pancreatic Cancer Cohort Consortium nested case-control study (1,481 cases, 1,539 controls). Odds ratios and 95% confidence intervals were calculated by using unconditional logistic regression. Smoking intensity effects were examined with an excess odds ratio model that was linear in pack-years and exponential in cigarettes smoked per day and its square. When compared with never smokers, current smokers had a significantly elevated risk (odds ratio (OR) = 1.77, 95% confidence interval (CI): 1.38, 2.26). Risk increased significantly with greater intensity (> or =30 cigarettes/day: OR = 1.75, 95% CI: 1.27, 2.42), duration (> or =50 years: OR = 2.13, 95% CI: 1.25, 3.62), and cumulative smoking dose (> or =40 pack-years: OR = 1.78, 95% CI: 1.35, 2.34). Risk more than 15 years after smoking cessation was similar to that for never smokers. Estimates of excess odds ratio per pack-year declined with increasing intensity, suggesting greater risk for total exposure delivered at lower intensity for longer duration than for higher intensity for shorter duration. This finding and the decline in risk after smoking cessation suggest that smoking has a late-stage effect on pancreatic carcinogenesis.

Genome-wide association study identifies variants in the ABO locus associated with susceptibility to pancreatic cancer

We conducted a two-stage genome-wide association study of pancreatic cancer, a cancer with one of the lowest survival rates worldwide. We genotyped 558,542 SNPs in 1,896 individuals with pancreatic cancer and 1,939 controls drawn from 12 prospective cohorts plus one hospital-based case-control study. We conducted a combined analysis of these groups plus an additional 2,457 affected individuals and 2,654 controls from eight case-control studies, adjusting for study, sex, ancestry and five principal components. We identified an association between a locus on 9q34 and pancreatic cancer marked by the SNP rs505922 (combined P = 5.37 x 10(-8); multiplicative per-allele odds ratio 1.20; 95% confidence interval 1.12-1.28). This SNP maps to the first intron of the ABO blood group gene. Our results are consistent with earlier epidemiologic evidence suggesting that people with blood group O may have a lower risk of pancreatic cancer than those with groups A or B.

Comprehensive resequence analysis of a 136 kb region of human chromosome 8q24 associated with prostate and colon cancers

Recently, genome-wide association studies have identified loci across a segment of chromosome 8q24 (128,100,000–128,700,000) associated with the risk of breast, colon and prostate cancers. At least three regions of 8q24 have been independently associated with prostate cancer risk; the most centromeric of which appears to be population specific. Haplotypes in two contiguous but independent loci, marked by rs6983267 and rs1447295, have been identified in the Cancer Genetic Markers of Susceptibility project (http://cgems.cancer.gov), which genotyped more than 5,000 prostate cancer cases and 5,000 controls of European origin. The rs6983267 locus is also strongly associated with colorectal cancer. To ascertain a comprehensive catalog of common single-nucleotide polymorphisms (SNPs) across the two regions, we conducted a resequence analysis of 136 kb (chr8: 128,473,000–128,609,802) using the Roche/454 next-generation sequencing technology in 39 prostate cancer cases and 40 controls of European origin. We have characterized a comprehensive catalog of common (MAF > 1%) SNPs within this region, including 442 novel SNPs and have determined the pattern of linkage disequilibrium across the region. Our study has generated a detailed map of genetic variation across the region, which should be useful for choosing SNPs for fine mapping of association signals in 8q24 and investigations of the functional consequences of select common variants.

Testicular germ cell tumours in Iceland

The purpose of this study was to examine the pathology of all germ cell tumours of the testis diagnosed in Iceland 1955-2002. A total of 214 patients were included in the study. The current age-standardized incidence was found to be 6.1 per 100,000 and had increased almost fourfold during the study period. Seminoma was diagnosed in 55% of cases. Non-seminomas were diagnosed in 45%, and these were further classified as mixed germ cell tumours (33%), embryonal carcinoma (8%), teratoma (3%), and yolk sac tumour (n=1). The mean age at diagnosis was significantly higher for the seminomas than the non-seminomas (38 years versus 29 years) (p<0.001) and the non-seminomas were diagnosed at a significantly higher stage than the seminomas (p<0.001). Thus, in seminoma patients the tumour was localized to the testis (stage I) in 81% of cases, in 17% of patients the tumour had spread to the lymph nodes (stage II or III), and only 2% had extranodal metastasis at diagnosis (stage IV). In contrast, in the non-seminoma patients, the tumours were found to be stage I in 56%, stage II or III in 24%, and stage IV in 20% of cases. No significant difference in staging was found between non-seminoma subtypes. Identification of necrosis or vascular invasion was significantly associated with metastatic disease at diagnosis (p=0.002). During the study period a significant increase in stage I tumours was found as well as a decrease in the size of the tumours.

Mechanism of inhibition of the Drosophila and mammalian EGF receptors by the transmembrane protein Kekkon 1

The transmembrane protein Kekkon 1 (Kek1) has previously been shown to act in a negative feedback loop to downregulate the Drosophila Epidermal Growth Factor Receptor (DER) during oogenesis. We show that this protein plays a similar role in other DER-mediated developmental processes. Structure-function analysis reveals that the extracellular Leucine-Rich Repeat (LRR) domains of Kek1 are critical for its function through direct association with DER, whereas its cytoplasmic domain is required for apical subcellular localization. In addition, the use of chimeric proteins between Kek1 extracellular and transmembrane domains fused to DER intracellular domain indicates that Kek1 forms an heterodimer with DER in vivo. To characterize more precisely the mechanism underlying the Kek1/DER interaction, we used mammalian ErbB/EGFR cell-based assays. We show that Kek1 is capable of physically interacting with each of the known members of the mammalian ErbB receptor family and that the Kek1/EGFR interaction inhibits growth factor binding, receptor autophosphorylation and Erk1/2 activation in response to EGF. Finally, in vivo experiments show that Kek1 expression potently suppresses the growth of mouse mammary tumor cells derived from aberrant ErbB receptors activation, but does not interfere with the growth of tumor cells derived from activated Ras. Our results underscore the possibility that Kek1 may be used experimentally to inhibit ErbB receptors and point to the possibility that, as yet uncharacterized, mammalian transmembrane LRR proteins might act as modulators of growth factor signalling.