Coffee consumption and risk of endometrial cancer: a pooled analysis of individual participant data in the Epidemiology of Endometrial Cancer Consortium (E2C2)

BACKGROUND

Epidemiologic studies suggest that coffee consumption may be inversely associated with risk of endometrial cancer (EC), the most common gynecological malignancy in developed countries. Furthermore, coffee consumption may lower circulating concentrations of estrogen and insulin, hormones implicated in endometrial carcinogenesis. Antioxidants and other chemopreventive compounds in coffee may have anticarcinogenic effects. Based on available meta-analyses, the World Cancer Research Fund (WCRF) concluded that consumption of coffee probably protects against EC.

OBJECTIVES

Our main aim was to examine the association between coffee consumption and EC risk by combining individual-level data in a pooled analysis. We also sought to evaluate potential effect modification by other risk factors for EC.

METHODS

We combined individual-level data from 19 epidemiologic studies (6 cohort, 13 case-control) of 12,159 EC cases and 27,479 controls from the Epidemiology of Endometrial Cancer Consortium (E2C2). Logistic regression was used to calculate ORs and their corresponding 95% CIs. All models were adjusted for potential confounders including age, race, BMI, smoking status, diabetes status, study design, and study site.

RESULTS

Coffee drinkers had a lower risk of EC than non-coffee drinkers (multiadjusted OR: 0.87; 95% CI: 0.79, 0.95). There was a dose-response relation between higher coffee consumption and lower risk of EC: compared with non-coffee drinkers, the adjusted pooled ORs for those who drank 1, 2-3, and >4 cups/d were 0.90 (95% CI: 0.82, 1.00), 0.86 (95% CI: 0.78, 0.95), and 0.76 (95% CI: 0.66, 0.87), respectively (P-trend < 0.001). The inverse association between coffee consumption and EC risk was stronger in participants with BMI > 25 kg/m2.

CONCLUSIONS

The results of the largest analysis to date pooling individual-level data further support the potentially beneficial health effects of coffee consumption in relation to EC, especially among females with higher BMI.

Polygenic risk modeling for prediction of epithelial ovarian cancer risk

Polygenic risk scores (PRS) for epithelial ovarian cancer (EOC) have the potential to improve risk stratification. Joint estimation of Single Nucleotide Polymorphism (SNP) effects in models could improve predictive performance over standard approaches of PRS construction. Here, we implemented computationally efficient, penalized, logistic regression models (lasso, elastic net, stepwise) to individual level genotype data and a Bayesian framework with continuous shrinkage, "select and shrink for summary statistics" (S4), to summary level data for epithelial non-mucinous ovarian cancer risk prediction. We developed the models in a dataset consisting of 23,564 non-mucinous EOC cases and 40,138 controls participating in the Ovarian Cancer Association Consortium (OCAC) and validated the best models in three populations of different ancestries: prospective data from 198,101 women of European ancestries; 7,669 women of East Asian ancestries; 1,072 women of African ancestries, and in 18,915 BRCA1 and 12,337 BRCA2 pathogenic variant carriers of European ancestries. In the external validation data, the model with the strongest association for non-mucinous EOC risk derived from the OCAC model development data was the S4 model (27,240 SNPs) with odds ratios (OR) of 1.38 (95% CI: 1.28-1.48, AUC: 0.588) per unit standard deviation, in women of European ancestries; 1.14 (95% CI: 1.08-1.19, AUC: 0.538) in women of East Asian ancestries; 1.38 (95% CI: 1.21-1.58, AUC: 0.593) in women of African ancestries; hazard ratios of 1.36 (95% CI: 1.29-1.43, AUC: 0.592) in BRCA1 pathogenic variant carriers and 1.49 (95% CI: 1.35-1.64, AUC: 0.624) in BRCA2 pathogenic variant carriers. Incorporation of the S4 PRS in risk prediction models for ovarian cancer may have clinical utility in ovarian cancer prevention programs.

Pregnancy outcomes and risk of endometrial cancer: A pooled analysis of individual participant data in the Epidemiology of Endometrial Cancer Consortium

A full-term pregnancy is associated with reduced endometrial cancer risk; however, whether the effect of additional pregnancies is independent of age at last pregnancy is unknown. The associations between other pregnancy-related factors and endometrial cancer risk are less clear. We pooled individual participant data from 11 cohort and 19 case-control studies participating in the Epidemiology of Endometrial Cancer Consortium (E2C2) including 16 986 women with endometrial cancer and 39 538 control women. We used one- and two-stage meta-analytic approaches to estimate pooled odds ratios (ORs) for the association between exposures and endometrial cancer risk. Ever having a full-term pregnancy was associated with a 41% reduction in risk of endometrial cancer compared to never having a full-term pregnancy (OR = 0.59, 95% confidence interval [CI] 0.56-0.63). The risk reduction appeared the greatest for the first full-term pregnancy (OR = 0.78, 95% CI 0.72-0.84), with a further ~15% reduction per pregnancy up to eight pregnancies (OR = 0.20, 95% CI 0.14-0.28) that was independent of age at last full-term pregnancy. Incomplete pregnancy was also associated with decreased endometrial cancer risk (7%-9% reduction per pregnancy). Twin births appeared to have the same effect as singleton pregnancies. Our pooled analysis shows that, while the magnitude of the risk reduction is greater for a full-term pregnancy than an incomplete pregnancy, each additional pregnancy is associated with further reduction in endometrial cancer risk, independent of age at last full-term pregnancy. These results suggest that the very high progesterone level in the last trimester of pregnancy is not the sole explanation for the protective effect of pregnancy.

Lack of association between modifiable exposures and glioma risk: a Mendelian randomization analysis

Background

The etiological basis of glioma is poorly understood. We have used genetic markers in a Mendelian randomization (MR) framework to examine if lifestyle, cardiometabolic, and inflammatory factors influence the risk of glioma. This methodology reduces bias from confounding and is not affected by reverse causation.

Methods

We identified genetic instruments for 37 potentially modifiable risk factors and evaluated their association with glioma risk using data from a genome-wide association study of 12 488 glioma patients and 18 169 controls. We used the estimated odds ratio of glioma associated with each of the genetically defined traits to infer evidence for a causal relationship with the following exposures:

Lifestyle and dietary factors—height, plasma insulin-like growth factor 1, blood carnitine, blood methionine, blood selenium, blood zinc, circulating adiponectin, circulating carotenoids, iron status, serum calcium, vitamins (A1, B12, B6, E, and 25-hydroxyvitamin D), fatty acid levels (monounsaturated, omega-3, and omega-6) and circulating fetuin-A;

Cardiometabolic factors—birth weight, high density lipoprotein cholesterol, low density lipoprotein cholesterol, total cholesterol, total triglycerides, basal metabolic rate, body fat percentage, body mass index, fasting glucose, fasting proinsulin, glycated hemoglobin levels, diastolic and systolic blood pressure, waist circumference, waist-to-hip ratio; and

Inflammatory factors— C-reactive protein, plasma interleukin-6 receptor subunit alpha and serum immunoglobulin E.

Results

After correction for the testing of multiple potential risk factors and excluding associations driven by one single nucleotide polymorphism, no significant association with glioma risk was observed (ie, P_Corrected > 0.05).

Conclusions

This study did not provide evidence supporting any of the 37 factors examined as having a significant influence on glioma risk.

The p.Ser64Leu and p.Pro104Leu missense variants of PALB2 identified in familial pancreatic cancer patients compromise the DNA damage response

PALB2 has been identified as a breast and pancreatic cancer susceptibility gene. Utilizing a targeted sequencing approach, we discovered two novel germline missense PALB2 variants c.191C>T and c.311C>T, encoding p.Ser64Leu and p.Pro104Leu, respectively, in individuals in a pancreatic cancer registry. No missense PALB2 variants from familial pancreatic cancer patients, and few PALB2 variants overall, have been functionally characterized. Given the known role of PALB2, we tested the impact of p.Ser64Leu and p.Pro104Leu variants on DNA damage responses. Neither p.Ser64Leu nor p.Pro104Leu have clear effects on interactions with BRCA1 and KEAP1, which are mediated by adjacent motifs in PALB2. However, both variants are associated with defective recruitment of PALB2, and the RAD51 recombinase downstream, to DNA damage foci. Furthermore, p.Ser64Leu and p.Pro104Leu both largely compromise DNA double-strand break-initiated homologous recombination, and confer increased cellular sensitivity to ionizing radiation (IR) and the poly (ADP-ribose) polymerase (PARP) inhibitor Olaparib. Taken together, our results represent the first demonstration of functionally deleterious PALB2 missense variants associated with familial pancreatic cancer and of deleterious variants in the N-terminus outside of the coiled-coil domain. Furthermore, our results suggest the possibility of personalized treatments, using IR or PARP inhibitor, of pancreatic and other cancers that carry a deleterious PALB2 variant.

Incorporating multiple sets of eQTL weights into gene-by-environment interaction analysis identifies novel susceptibility loci for pancreatic cancer

It is of great scientific interest to identify interactions between genetic variants and environmental exposures that may modify the risk of complex diseases. However, larger sample sizes are usually required to detect gene-by-environment interaction (G × E) than required to detect genetic main association effects. To boost the statistical power and improve the understanding of the underlying molecular mechanisms, we incorporate functional genomics information, specifically, expression quantitative trait loci (eQTLs), into a data-adaptive G × E test, called aGEw. This test adaptively chooses the best eQTL weights from multiple tissues and provides an extra layer of weighting at the genetic variant level. Extensive simulations show that the aGEw test can control the Type 1 error rate, and the power is resilient to the inclusion of neutral variants and noninformative external weights. We applied the proposed aGEw test to the Pancreatic Cancer Case-Control Consortium (discovery cohort of 3,585 cases and 3,482 controls) and the PanScan II genome-wide association study data (replication cohort of 2,021 cases and 2,105 controls) with smoking as the exposure of interest. Two novel putative smoking-related pancreatic cancer susceptibility genes, TRIP10 and KDM3A, were identified. The aGEw test is implemented in an R package aGE.

Gallbladder disease, cholecystectomy, and pancreatic cancer risk in the International Pancreatic Cancer Case-Control Consortium (PanC4)

BACKGROUND

The association among gallbladder disease, cholecystectomy, and pancreatic cancer is unclear. Moreover, time interval between gallbladder disease or cholecystectomy and pancreatic cancer diagnosis is not considered in most previous studies.

AIM

To quantify the association among gallbladder disease, cholecystectomy, and pancreatic cancer, considering time since first diagnosis of gallbladder disease or cholecystectomy.

METHODS

We used data from nine case-control studies within the Pancreatic Cancer Case-Control Consortium, including 5760 cases of adenocarcinoma of the exocrine pancreas and 8437 controls. We estimated pooled odds ratios and the corresponding 95% confidence intervals by estimating study-specific odds ratios through multivariable unconditional logistic regression models, and then pooling the obtained estimates using fixed-effects models.

RESULTS

Compared with patients with no history of gallbladder disease, the pooled odds ratio of pancreatic cancer was 1.69 (95% confidence interval, 1.51-1.88) for patients reporting a history of gallbladder disease. The odds ratio was 4.90 (95% confidence interval, 3.45-6.97) for gallbladder disease diagnosed <2 years before pancreatic cancer diagnosis and 1.11 (95% confidence interval, 0.96-1.29) when ≥2 years elapsed. The pooled odds ratio was 1.64 (95% confidence interval, 1.43-1.89) for patients who underwent cholecystectomy, as compared to those without cholecystectomy. The odds ratio was 7.00 (95% confidence interval, 4.13-11.86) for a surgery <2 years before pancreatic cancer diagnosis and 1.28 (95% confidence interval, 1.08-1.53) for a surgery ≥2 years before.

CONCLUSIONS

There appears to be no long-term effect of gallbladder disease on pancreatic cancer risk, and at most a modest one for cholecystectomy. The strong short-term association can be explained by diagnostic bias and reverse causation.

A pooled genome-wide association study identifies pancreatic cancer susceptibility loci on chromosome 19p12 and 19p13.3 in the full-Jewish population

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.

Detection of first gammaherpesvirus sequences in Central African bats

Herpesviruses have been identified in many species; however, relatively few bat herpesvirus are known, considering the enormous diversity of bats. We used consensus PCR to test bats from the Republic of the Congo and found DNA of two different novel bat herpesviruses. One was detected in a Pipistrellus nanulus, the other in a Triaenops persicus bat and both resemble gammaherpesviruses. On the amino acid level, the amplified sequences differ by 55% from each other, and by 27% and 25% from the next closest known viruses. The findings point towards the diversity of herpesviruses in Central African bats.

A region-based gene association study combined with a leave-one-out sensitivity analysis identifies SMG1 as a pancreatic cancer susceptibility gene

Pancreatic adenocarcinoma (PC) is a lethal malignancy that is familial or associated with genetic syndromes in 10% of cases. Gene-based surveillance strategies for at-risk individuals may improve clinical outcomes. However, familial PC (FPC) is plagued by genetic heterogeneity and the genetic basis for the majority of FPC remains elusive, hampering the development of gene-based surveillance programs. The study was powered to identify genes with a cumulative pathogenic variant prevalence of at least 3%, which includes the most prevalent PC susceptibility gene, BRCA2. Since the majority of known PC susceptibility genes are involved in DNA repair, we focused on genes implicated in these pathways. We performed a region-based association study using the Mixed-Effects Score Test, followed by leave-one-out characterization of PC-associated gene regions and variants to identify the genes and variants driving risk associations. We evaluated 398 cases from two case series and 987 controls without a personal history of cancer. The first case series consisted of 109 patients with either FPC (n = 101) or PC at ≤50 years of age (n = 8). The second case series was composed of 289 unselected PC cases. We validated this discovery strategy by identifying known pathogenic BRCA2 variants, and also identified SMG1, encoding a serine/threonine protein kinase, to be significantly associated with PC following correction for multiple testing (p = 3.22x10-7). The SMG1 association was validated in a second independent series of 532 FPC cases and 753 controls (p<0.0062, OR = 1.88, 95%CI 1.17-3.03). We showed segregation of the c.4249A>G SMG1 variant in 3 affected relatives in a FPC kindred, and we found c.103G>A to be a recurrent SMG1 variant associating with PC in both the discovery and validation series. These results suggest that SMG1 is a novel PC susceptibility gene, and we identified specific SMG1 gene variants associated with PC risk.

Outcome of Pancreatic Cancer Surveillance Among High-Risk Individuals Tested for Germline Mutations in BRCA1 and BRCA2

Germline mutations in BRCA1/2 are risk factors for pancreatic ductal adenocarcinoma (PDAC). The aim of this study was to evaluate whether results of surveillance for PDAC in high risk individuals (HRI) differ between those with and without a pathogenic BRCA1/2 mutation. This prospective study was conducted within the Pancreatic Tumor Registry at a major cancer center. There were 83 HRIs with ≥1 first-degree relative with PDAC who underwent surveillance and testing for pathogenic germline mutations in BRCA1/2 A secondary analysis includes 18 HRIs with known mutations in BRCA1/2 but with weaker family history. HRIs were evaluated over time using magnetic resonance cholangiopancreatography (MRCP) and endoscopic ultrasound when indicated by MRCP findings. We reviewed imaging results, blinded to mutation status. Demographic information was obtained from interviewer-administered questionnaires. The outcome was the proportion with any pancreatic abnormality identified at initial or follow-up surveillance. Among the 83 HRIs in the main analysis, 48 had a mutation in BRCA1/2 and 35 did not. Overall, 16 of 48 (33%) BRCA1/2-positive and 13 of 35 (37%) BRCA1/2-negative participants had pancreatic abnormalities on imaging; in each group, all but one finding was an intraductal papillary mucinous neoplasm. Among those with pathogenic mutations but weaker family history, results were similar: 7 of 18 (39%) with pancreatic abnormalities. Results of surveillance for pancreatic abnormalities on imaging are similar regardless of BRCA1/2 mutation status. While the results from this small study need confirmation in other studies, at present there does not appear to be increased yield from targeting individuals with BRCA1/2 mutations for surveillance.

Analysis of Heritability and Genetic Architecture of Pancreatic Cancer: A PanC4 Study

BACKGROUND

Pancreatic cancer is the fourth-leading cause of cancer death in both men and women in the United States. The currently identified common susceptibility loci account for a small fraction of estimated heritability. We sought to estimate overall heritability of pancreatic cancer and partition the heritability by variant frequencies and functional annotations.

METHODS

Analysis using the genome-based restricted maximum likelihood method (GREML) was conducted on Pancreatic Cancer Case-Control Consortium (PanC4) genome-wide association study (GWAS) data from 3,568 pancreatic cancer cases and 3,363 controls of European Ancestry.

RESULTS

Applying linkage disequilibrium- and minor allele frequency-stratified GREML (GREML-LDMS) method to imputed GWAS data, we estimated the overall heritability of pancreatic cancer to be 21.2% (SE = 4.8%). Across the functional groups (intronic, intergenic, coding, and regulatory variants), intronic variants account for most of the estimated heritability (12.4%). Previously identified GWAS loci explained 4.1% of the total phenotypic variation of pancreatic cancer. Mutations in hereditary pancreatic cancer susceptibility genes are present in 4% to 10% of patients with pancreatic cancer, yet our GREML-LDMS results suggested these regions explain only 0.4% of total phenotypic variance for pancreatic cancer.

CONCLUSIONS

Although higher than previous studies, our estimated 21.2% overall heritability may still be downwardly biased due to the inherent limitation that the contribution of rare variants in genes with a substantive overall impact on disease are not captured when applying these commonly used methods to imputed GWAS data.

IMPACT

Our work demonstrated the importance of rare and common variants in pancreatic cancer risk.

Abstract 4173: Previously identified common glioma risk SNPs are associated with familial glioma

BACKGROUND: Approximately 5% of gliomas occur in individuals with a family history of glioma, and first-degree relatives of brain tumor cases have a two-fold increase in risk of brain tumor. Recent somatic characterization has shown that tumors from familial cases are indistinguishable from sporadic cases, suggesting that familial cases may arise through similar mechanisms of gliomagenesis, and therefore may be associated with common variants as well as rare mutations. In this analysis, we assessed whether previously identified common risk variants are associated with familial glioma.

METHODS: Data were obtained from the Glioma International Case Control (GICC) and Gliogene studies for 448 cases with reported family history, 4,405 cases without reported family history, and 3,288 controls. We assessed 25 risk loci previously identified by glioma GWAS, and odds ratios (OR) and 95% confidence intervals (95%CI) were calculated using an additive genetic logistic regression model adjusted for age, sex, and the first two principal components for familial cases versus unaffected controls, and non-familial cases versus controls. Results were considered significant at p&lt;0.002 (Bonferroni correction for 25 tests).

RESULTS: Significant associations were detected at 5/25 loci, including TERT, EGFR, CCDC26, CDKN2B, and RTEL1. The strongest association was at rs55705857 (CCDC26, OR=2.7, p=7.49x10-17). For GBM (222 familial cases), significant associations were detected at 6/26 loci (TERT, EGFR, CDKN2B, TP53 and RTEL1), while in non-GBM (205 familial cases) significant associations were detected at 3/25 loci (LRIG1, CCDC26, PHLDB1). These SNPs were further examined using a case-only approach comparing familial to non-familial cases, and there was no significant difference in allele frequencies by family history status. There was a strong correlation between log(OR) for familial cases only versus non-familial cases (adjusted R2=0.88).

CONCLUSIONS: In this analysis we identified a significant association between familial glioma and five common risk loci previously identified by glioma GWAS. This provides further evidence of shared pathways of genetic risk and gliomagenesis between familial and non-familial glioma. Further exploration is necessary to determine the overall contribution of common genetic variation to risk of familial glioma.

Citation Format: Quinn T. Ostrom, Georgina Armstrong, Christopher I. Amos, Jonine L. Bernstein, Elizabeth B. Claus, Jeanette E. Eckel-Passow, Dora Il'yasova, Christoffer Johansen, Daniel H. Lachance, Rose K. Lai, Ryan T. Merrell, Sara H. Olson, Joellen H. Schildkraut, Sanjay S. Shete, Richard S. Houlston, Robert B. Jenkins, Margaret R. Wrensch, Beatrice Melin, Jill S. Barnholtz-Sloan, Melissa L. Bondy. Previously identified common glioma risk SNPs are associated with familial glioma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 4173.

Agnostic Pathway/Gene Set Analysis of Genome-Wide Association Data Identifies Associations for Pancreatic Cancer

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.

A comprehensive gene-environment interaction analysis in Ovarian Cancer using genome-wide significant common variants

As a follow-up to genome-wide association analysis of common variants associated with ovarian carcinoma (cancer), our study considers seven well-known ovarian cancer risk factors and their interactions with 28 genome-wide significant common genetic variants. The interaction analyses were based on data from 9971 ovarian cancer cases and 15,566 controls from 17 case-control studies. Likelihood ratio and Wald tests for multiplicative interaction and for relative excess risk due to additive interaction were used. The top multiplicative interaction was noted between oral contraceptive pill (OCP) use (ever vs. never) and rs13255292 (p value = 3.48 × 10-4 ). Among women with the TT genotype for this variant, the odds ratio for OCP use was 0.53 (95% CI = 0.46-0.60) compared to 0.71 (95%CI = 0.66-0.77) for women with the CC genotype. When stratified by duration of OCP use, women with 1-5 years of OCP use exhibited differential protective benefit across genotypes. However, no interaction on either the multiplicative or additive scale was found to be statistically significant after multiple testing correction. The results suggest that OCP use may offer increased benefit for women who are carriers of the T allele in rs13255292. On the other hand, for women carrying the C allele in this variant, longer (5+ years) use of OCP may reduce the impact of carrying the risk allele of this SNP. Replication of this finding is needed. The study presents a comprehensive analytic framework for conducting gene-environment analysis in ovarian cancer.

Glioma risk associated with extent of estimated European genetic ancestry in African Americans and Hispanics

Glioma incidence is highest in non-Hispanic Whites, and to date, glioma genome-wide association studies (GWAS) to date have only included European ancestry (EA) populations. African Americans and Hispanics in the US have varying proportions of EA, African (AA) and Native American ancestries (NAA). It is unknown if identified GWAS loci or increased EA is associated with increased glioma risk. We assessed whether EA was associated with glioma in African Americans and Hispanics. Data were obtained for 832 cases and 675 controls from the Glioma International Case-Control Study and GliomaSE Case-Control Study previously estimated to have <80% EA, or self-identify as non-White. We estimated global and local ancestry using fastStructure and RFMix, respectively, using 1,000 genomes project reference populations. Within groups with ≥40% AA (AFR_≥0.4 ), and ≥15% NAA (AMR_≥0.15 ), genome-wide association between local EA and glioma was evaluated using logistic regression conditioned on global EA for all gliomas. We identified two regions (7q21.11, p = 6.36 × 10^-4 ; 11p11.12, p = 7.0 × 10^-4 ) associated with increased EA, and one associated with decreased EA (20p12.13, p = 0.0026) in AFR_≥0.4 . In addition, we identified a peak at rs1620291 (p = 4.36 × 10^-6 ) in 7q21.3. Among AMR_≥0.15 , we found an association between increased EA in one region (12q24.21, p = 8.38 × 10^-4 ), and decreased EA in two regions (8q24.21, p = 0. 0010; 20q13.33, p = 6.36 × 10^-4 ). No other significant associations were identified. This analysis identified an association between glioma and two regions previously identified in EA populations (8q24.21, 20q13.33) and four novel regions (7q21.11, 11p11.12, 12q24.21 and 20p12.13). The identifications of novel association with EA suggest regions to target for future genetic association studies.

Aspirin, NSAIDs, and Glioma Risk: Original Data from the Glioma International Case-Control Study and a Meta-analysis

BACKGROUND

There have been few studies of sufficient size to address the relationship between glioma risk and the use of aspirin or NSAIDs, and results have been conflicting. The purpose of this study was to examine the associations between glioma and aspirin/NSAID use, and to aggregate these findings with prior published studies using meta-analysis.

METHODS

The Glioma International Case-Control Study (GICC) consists of 4,533 glioma cases and 4,171 controls recruited from 2010 to 2013. Interviews were conducted using a standardized questionnaire to obtain information on aspirin/NSAID use. We examined history of regular use for ≥6 months and duration-response. Restricted maximum likelihood meta-regression models were used to aggregate site-specific estimates, and to combine GICC estimates with previously published studies.

RESULTS

A history of daily aspirin use for ≥6 months was associated with a 38% lower glioma risk, compared with not having a history of daily use [adjusted meta-OR = 0.62; 95% confidence interval (CI), 0.54-0.70]. There was a significant duration-response trend (P = 1.67 × 10^-17), with lower ORs for increasing duration of aspirin use. Duration-response trends were not observed for NSAID use. In the meta-analysis aggregating GICC data with five previous studies, there was a marginally significant association between use of aspirin and glioma (mOR = 0.84; 95% CI, 0.70-1.02), but no association for NSAID use.

CONCLUSIONS

Our study suggests that aspirin may be associated with a reduced risk of glioma.

IMPACT

These results imply that aspirin use may be associated with decreased glioma risk. Further research examining the association between aspirin use and glioma risk is warranted.

Transcriptome-Wide Association Study Identifies New Candidate Susceptibility Genes for Glioma

Genome-wide association studies (GWAS) have so far identified 25 loci associated with glioma risk, with most showing specificity for either glioblastoma (GBM) or non-GBM tumors. The majority of these GWAS susceptibility variants reside in noncoding regions and the causal genes underlying the associations are largely unknown. Here we performed a transcriptome-wide association study to search for novel risk loci and candidate causal genes at known GWAS loci using Genotype-Tissue Expression Project (GTEx) data to predict cis-predicted gene expression in relation to GBM and non-GBM risk in conjunction with GWAS summary statistics on 12,488 glioma cases (6,183 GBM and 5,820 non-GBM) and 18,169 controls. Imposing a Bonferroni-corrected significance level of P < 5.69 × 10^-6, we identified 31 genes, including GALNT6 at 12q13.33, as a candidate novel risk locus for GBM (mean Z = 4.43; P = 5.68 × 10^-6). GALNT6 resides at least 55 Mb away from any previously identified glioma risk variant, while all other 30 significantly associated genes were located within 1 Mb of known GWAS-identified loci and were not significant after conditioning on the known GWAS-identified variants. These data identify a novel locus (GALNT6 at 12q13.33) and 30 genes at 12 known glioma risk loci associated with glioma risk, providing further insights into glioma tumorigenesis. SIGNIFICANCE: This study identifies new genes associated with glioma risk, increasing understanding of how these tumors develop.

Sex-specific gene and pathway modeling of inherited glioma risk

Background

To date, genome-wide association studies (GWAS) have identified 25 risk variants for glioma, explaining 30% of heritable risk. Most histologies occur with significantly higher incidence in males, and this difference is not explained by currently known risk factors. A previous GWAS identified sex-specific glioma risk variants, and this analysis aims to further elucidate risk variation by sex using gene- and pathway-based approaches.

Methods

Results from the Glioma International Case-Control Study were used as a testing set, and results from 3 GWAS were combined via meta-analysis and used as a validation set. Using summary statistics for nominally significant autosomal SNPs (P < 0.01 in a previous meta-analysis) and nominally significant X-chromosome SNPs (P < 0.01), 3 algorithms (Pascal, BimBam, and GATES) were used to generate gene scores, and Pascal was used to generate pathway scores. Results were considered statistically significant in the discovery set when P < 3.3 × 10-6 and in the validation set when P < 0.001 in 2 of 3 algorithms.

Results

Twenty-five genes within 5 regions and 19 genes within 6 regions reached statistical significance in at least 2 of 3 algorithms in males and females, respectively. EGFR was significantly associated with all glioma and glioblastoma in males only and a female-specific association in TERT, all of which remained nominally significant after conditioning on known risk loci. There were nominal associations with the BioCarta telomeres pathway in both males and females.

Conclusions

These results provide additional evidence that there may be differences by sex in genetic risk for glioma. Additional analyses may further elucidate the biological processes through which this risk is conferred.

Age-specific genome-wide association study in glioblastoma identifies increased proportion of 'lower grade glioma'-like features associated with younger age

Glioblastoma (GBM) is the most common malignant brain tumor in the United States. Incidence of GBM increases with age, and younger age-at-diagnosis is significantly associated with improved prognosis. While the relationship between candidate GBM risk SNPs and age-at-diagnosis has been explored, genome-wide association studies (GWAS) have not previously been stratified by age. Potential age-specific genetic effects were assessed in autosomal SNPs for GBM patients using data from four previous GWAS. Using age distribution tertiles (18-53, 54-64, 65+) datasets were analyzed using age-stratified logistic regression to generate p values, odds ratios (OR), and 95% confidence intervals (95%CI), and then combined using meta-analysis. There were 4,512 total GBM cases, and 10,582 controls used for analysis. Significant associations were detected at two previously identified SNPs in 7p11.2 (rs723527 [p54-63 = 1.50x10-9 , OR54-63 = 1.28, 95%CI54-63 = 1.18-1.39; p64+ = 2.14x10-11 , OR64+ = 1.32, 95%CI64+ = 1.21-1.43] and rs11979158 [p54-63 = 6.13x10-8 , OR54-63 = 1.35, 95%CI54-63 = 1.21-1.50; p64+ = 2.18x10-10 , OR64+ = 1.42, 95%CI64+ = 1.27-1.58]) but only in persons >54. There was also a significant association at the previously identified lower grade glioma (LGG) risk locus at 8q24.21 (rs55705857) in persons ages 18-53 (p18-53 = 9.30 × 10-11 , OR18-53 = 1.76, 95%CI18-53 = 1.49-2.10). Within The Cancer Genome Atlas (TCGA) there was higher prevalence of 'LGG'-like tumor characteristics in GBM samples in those 18-53, with IDH1/2 mutation frequency of 15%, as compared to 2.1% [54-63] and 0.8% [64+] (p = 0.0005). Age-specific differences in cancer susceptibility can provide important clues to etiology. The association of a SNP known to confer risk for IDH1/2 mutant glioma and higher prevalence of IDH1/2 mutation within younger individuals 18-53 suggests that more younger individuals may present initially with 'secondary glioblastoma.'

A Transcriptome-Wide Association Study Among 97,898 Women to Identify Candidate Susceptibility Genes for Epithelial Ovarian Cancer Risk

Large-scale genome-wide association studies (GWAS) have identified approximately 35 loci associated with epithelial ovarian cancer (EOC) risk. The majority of GWAS-identified disease susceptibility variants are located in noncoding regions, and causal genes underlying these associations remain largely unknown. Here, we performed a transcriptome-wide association study to search for novel genetic loci and plausible causal genes at known GWAS loci. We used RNA sequencing data (68 normal ovarian tissue samples from 68 individuals and 6,124 cross-tissue samples from 369 individuals) and high-density genotyping data from European descendants of the Genotype-Tissue Expression (GTEx V6) project to build ovarian and cross-tissue models of genetically regulated expression using elastic net methods. We evaluated 17,121 genes for their cis-predicted gene expression in relation to EOC risk using summary statistics data from GWAS of 97,898 women, including 29,396 EOC cases. With a Bonferroni-corrected significance level of P < 2.2 × 10-6, we identified 35 genes, including FZD4 at 11q14.2 (Z = 5.08, P = 3.83 × 10-7, the cross-tissue model; 1 Mb away from any GWAS-identified EOC risk variant), a potential novel locus for EOC risk. All other 34 significantly associated genes were located within 1 Mb of known GWAS-identified loci, including 23 genes at 6 loci not previously linked to EOC risk. Upon conditioning on nearby known EOC GWAS-identified variants, the associations for 31 genes disappeared and three genes remained (P < 1.47 × 10-3). These data identify one novel locus (FZD4) and 34 genes at 13 known EOC risk loci associated with EOC risk, providing new insights into EOC carcinogenesis.Significance: Transcriptomic analysis of a large cohort confirms earlier GWAS loci and reveals FZD4 as a novel locus associated with EOC risk. Cancer Res; 78(18); 5419-30. ©2018 AACR.

rs495139 in the TYMS-ENOSF1 Region and Risk of Ovarian Carcinoma of Mucinous Histology

Thymidylate synthase (TYMS) is a crucial enzyme for DNA synthesis. TYMS expression is regulated by its antisense mRNA, ENOSF1. Disrupted regulation may promote uncontrolled DNA synthesis and tumor growth. We sought to replicate our previously reported association between rs495139 in the TYMS-ENOSF1 3' gene region and increased risk of mucinous ovarian carcinoma (MOC) in an independent sample. Genotypes from 24,351 controls to 15,000 women with invasive OC, including 665 MOC, were available. We estimated per-allele odds ratios (OR) and 95% confidence intervals (CI) using unconditional logistic regression, and meta-analysis when combining these data with our previous report. The association between rs495139 and MOC was not significant in the independent sample (OR = 1.09; 95% CI = 0.97⁻1.22; p = 0.15; N = 665 cases). Meta-analysis suggested a weak association (OR = 1.13; 95% CI = 1.03⁻1.24; p = 0.01; N = 1019 cases). No significant association with risk of other OC histologic types was observed (p = 0.05 for tumor heterogeneity). In expression quantitative trait locus (eQTL) analysis, the rs495139 allele was positively associated with ENOSF1 mRNA expression in normal tissues of the gastrointestinal system, particularly esophageal mucosa (r = 0.51, p = 1.7 × 10-28), and nonsignificantly in five MOC tumors. The association results, along with inconclusive tumor eQTL findings, suggest that a true effect of rs495139 might be small.

Variants in genes encoding small GTPases and association with epithelial ovarian cancer susceptibility

Epithelial ovarian cancer (EOC) is the fifth leading cause of cancer mortality in American women. Normal ovarian physiology is intricately connected to small GTP binding proteins of the Ras superfamily (Ras, Rho, Rab, Arf, and Ran) which govern processes such as signal transduction, cell proliferation, cell motility, and vesicle transport. We hypothesized that common germline variation in genes encoding small GTPases is associated with EOC risk. We investigated 322 variants in 88 small GTPase genes in germline DNA of 18,736 EOC patients and 26,138 controls of European ancestry using a custom genotype array and logistic regression fitting log-additive models. Functional annotation was used to identify biofeatures and expression quantitative trait loci that intersect with risk variants. One variant, ARHGEF10L (Rho guanine nucleotide exchange factor 10 like) rs2256787, was associated with increased endometrioid EOC risk (OR = 1.33, p = 4.46 x 10-6). Other variants of interest included another in ARHGEF10L, rs10788679, which was associated with invasive serous EOC risk (OR = 1.07, p = 0.00026) and two variants in AKAP6 (A-kinase anchoring protein 6) which were associated with risk of invasive EOC (rs1955513, OR = 0.90, p = 0.00033; rs927062, OR = 0.94, p = 0.00059). Functional annotation revealed that the two ARHGEF10L variants were located in super-enhancer regions and that AKAP6 rs927062 was associated with expression of GTPase gene ARHGAP5 (Rho GTPase activating protein 5). Inherited variants in ARHGEF10L and AKAP6, with potential transcriptional regulatory function and association with EOC risk, warrant investigation in independent EOC study populations.

Abstract 233: Evaluating glioma risk associated with extent of European admixture in African-Americans and Latinos

Glioma incidence is highest in non-Hispanic Whites, where it occurs ~2x as frequently compared with other race/ethnicity groups. Glioma GWAS to date have included European ancestry populations only, and it is unknown whether variants identified by these analyses are associated with glioma in non- European ancestry populations. African Americans and Hispanics are admixed populations with varying proportions of European ancestry. While global ancestry may be similar within admixed groups, the proportion of European ancestry at each allele can vary across the genome. As glioma is more common in European ancestry populations, the presence of increased local European ancestry in these admixed populations could be used to identify glioma risk loci. Here we assessed whether excess European ancestry at established risk loci (Melin et al, Nature Genetics, 2017) was associated with glioma risk in non-European ancestry populations. Global ancestry was estimated using fastStructure, and local ancestry was estimated using RFMix. Both methods used 1,000 genomes project reference populations (African: YRI; European: CEU; East Asian: CHB/JPT; and Native American: CLM/PEL/MXL). We evaluated differences in local European ancestry between cases and controls using logistic regression conditioned on global European ancestry within 500kb of 25 previously identified risk variants among individuals with ≥50% African ancestry, and ≥30% Native American ancestry for all gliomas, and for grade IV glioblastoma (GBM) and grade II-III non-GBM. There were 347 individuals (184 cases and 163 controls) with ≥50% global African ancestry, and 277 individuals (153 cases and 124 controls) with ≥30% global American ancestry. There was no significant difference in proportion of global European ancestry between cases and controls with ≥50% global African ancestry (cases: 18.2%, controls: 17.7%, p=0.6834), and no significant difference in proportion of global European ancestry between cases and controls with ≥30% global American ancestry (cases: 51.1%, controls: 49.0%, p=0.2123). Among individuals with &gt;50% African ancestry, we observed a nominally significant association between all glioma and increased local European ancestry at 7p11.2 (EGFR, pmin=0.0070) and between GBM and increased local European ancestry at 22q13.1 (CSNK1E, pmin=0.0098), both near SNPs previously associated with glioblastoma in majority European-ancestry populations. The dataset used for this analysis represents the largest collection of genotyped non-European glioma cases. These results suggest that glioma risk in African Americans may be associated with an increased local European ancestry variants at glioma risk loci previously identified in majority European ancestry populations (7p11.2 and 22q13.1).

Citation Format: Quinn T. Ostrom, Kathleen M. Egan, L. Burt Nabors, Travis Gerke, Reid C. Thompson, Jeffrey J. Olson, Renato LaRocca, Sajeel Chowdhary, Jeanette E. Eckel-Passow, Georgina Armstrong, John K. Wiencke, Christopher I. Amos, Jonine L. Bernstein, Elizabeth B. Claus, Dora Il'yasova, Christoffer Johansen, Daniel H. Lachance, Rose Lai, Ryan T. Merrell, Sara H. Olson, Siegal Sadetzki, Joellen Schildkraut, Sanjay Shete, Richard S. Houlston, Robert B. Jenkins, Beatrice Melin, Melissa L. Bondy, Jill S. Barnholtz-Sloan. Evaluating glioma risk associated with extent of European admixture in African-Americans and Latinos [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 233.

Sex-specific glioma genome-wide association study identifies new risk locus at 3p21.31 in females, and finds sex-differences in risk at 8q24.21

Incidence of glioma is approximately 50% higher in males. Previous analyses have examined exposures related to sex hormones in women as potential protective factors for these tumors, with inconsistent results. Previous glioma genome-wide association studies (GWAS) have not stratified by sex. Potential sex-specific genetic effects were assessed in autosomal SNPs and sex chromosome variants for all glioma, GBM and non-GBM patients using data from four previous glioma GWAS. Datasets were analyzed using sex-stratified logistic regression models and combined using meta-analysis. There were 4,831 male cases, 5,216 male controls, 3,206 female cases and 5,470 female controls. A significant association was detected at rs11979158 (7p11.2) in males only. Association at rs55705857 (8q24.21) was stronger in females than in males. A large region on 3p21.31 was identified with significant association in females only. The identified differences in effect of risk variants do not fully explain the observed incidence difference in glioma by sex.

Glioma-related seizures in relation to histopathological subtypes: a report from the glioma international case-control study

BACKGROUND

The purpose of this study was to evaluate the distribution of glioma-related seizures and seizure control at the time of tumor diagnosis with respect to tumor histologic subtypes, tumor treatment and patient characteristics, and to compare seizure history preceding tumor diagnosis (or study enrollment) between glioma patients and healthy controls.

METHODS

The Glioma International Case Control study (GICC) risk factor questionnaire collected information on demographics, past medical/medication history, and occupational history. Cases from eight centers were also asked detailed questions on seizures in relation to glioma diagnosis; cases (n = 4533) and controls (n = 4171) were also asked about seizures less than 2 years from diagnosis and previous seizure history more than 2 years prior to tumor diagnosis, including childhood seizures.

RESULTS

Low-grade gliomas (LGGs), particularly oligodendrogliomas/oligoastrocytomas, had the highest proportion of glioma-related seizures. Patients with low-grade astrocytoma demonstrated the most medically refractory seizures. A total of 83% of patients were using only one antiepileptic drug (AED), which was levetiracetam in 71% of cases. Gross total resection was strongly associated with reduced seizure frequency (p < 0.009). No significant difference was found between glioma cases and controls in terms of seizure occurring more than 2 years before diagnosis or during childhood.

CONCLUSIONS

Our study showed that glioma-related seizures were most common in low-grade gliomas. Gross total resection was associated with lower seizure frequency. Additionally, having a history of childhood seizures is not a risk factor ***for developing glioma-related seizures or glioma.

Adult height is associated with increased risk of ovarian cancer: a Mendelian randomisation study

BACKGROUND

Observational studies suggest greater height is associated with increased ovarian cancer risk, but cannot exclude bias and/or confounding as explanations for this. Mendelian randomisation (MR) can provide evidence which may be less prone to bias.

METHODS

We pooled data from 39 Ovarian Cancer Association Consortium studies (16,395 cases; 23,003 controls). We applied two-stage predictor-substitution MR, using a weighted genetic risk score combining 609 single-nucleotide polymorphisms. Study-specific odds ratios (OR) and 95% confidence intervals (CI) for the association between genetically predicted height and risk were pooled using random-effects meta-analysis.

RESULTS

Greater genetically predicted height was associated with increased ovarian cancer risk overall (pooled-OR (pOR) = 1.06; 95% CI: 1.01-1.11 per 5 cm increase in height), and separately for invasive (pOR = 1.06; 95% CI: 1.01-1.11) and borderline (pOR = 1.15; 95% CI: 1.02-1.29) tumours.

CONCLUSIONS

Women with a genetic propensity to being taller have increased risk of ovarian cancer. This suggests genes influencing height are involved in pathways promoting ovarian carcinogenesis.

Characterising cis-regulatory variation in the transcriptome of histologically normal and tumour-derived pancreatic tissues

OBJECTIVE

To elucidate the genetic architecture of gene expression in pancreatic tissues.

DESIGN

We performed expression quantitative trait locus (eQTL) analysis in histologically normal pancreatic tissue samples (n=95) using RNA sequencing and the corresponding 1000 genomes imputed germline genotypes. Data from pancreatic tumour-derived tissue samples (n=115) from The Cancer Genome Atlas were included for comparison.

RESULTS

We identified 38 615 cis-eQTLs (in 484 genes) in histologically normal tissues and 39 713 cis-eQTL (in 237 genes) in tumour-derived tissues (false discovery rate <0.1), with the strongest effects seen near transcriptional start sites. Approximately 23% and 42% of genes with significant cis-eQTLs appeared to be specific for tumour-derived and normal-derived tissues, respectively. Significant enrichment of cis-eQTL variants was noted in non-coding regulatory regions, in particular for pancreatic tissues (1.53-fold to 3.12-fold, p≤0.0001), indicating tissue-specific functional relevance. A common pancreatic cancer risk locus on 9q34.2 (rs687289) was associated with ABO expression in histologically normal (p=5.8×10^-8) and tumour-derived (p=8.3×10^-5) tissues. The high linkage disequilibrium between this variant and the O blood group generating deletion variant in ABO (exon 6) suggested that nonsense-mediated decay (NMD) of the 'O' mRNA might explain this finding. However, knockdown of crucial NMD regulators did not influence decay of the ABO 'O' mRNA, indicating that a gene regulatory element influenced by pancreatic cancer risk alleles may underlie the eQTL.

CONCLUSIONS

We have identified cis-eQTLs representing potential functional regulatory variants in the pancreas and generated a rich data set for further studies on gene expression and its regulation in pancreatic tissues.

Transcriptional regulation by NR5A2 links differentiation and inflammation in the pancreas

Chronic inflammation increases the risk of several cancer types. The current notion is that the control of inflammatory responses relies on transcriptional networks distinct from those involved in cell differentiation ^1–3. The orphan nuclear receptor NR5A2 participates in a wide variety of processes including cholesterol and glucose metabolism in the liver, resolution of ER stress, intestinal glucocorticoid production, pancreatic development, and acinar differentiation ^4–8. Single nucleotide polymorphisms (SNPs) in the vicinity of NR5A2 have been associated with the risk of pancreatic adenocarcinoma (PDAC) through genome wide association studies ^9,10. In mice, Nr5a2 heterozygosity sensitizes the pancreas to damage, impairs regeneration, and cooperates with mutant KRas in tumor progression ¹¹. Through global transcriptomic analysis, we describe here an epithelial cell-autonomous basal pre-inflammatory state in the pancreas of Nr5a2^+/− mice that is reminiscent of early stages of pancreatitis-induced inflammation and is conserved in histologically normal human pancreata with reduced NR5A2 mRNA expression. In Nr5a2^+/− mice, Nr5a2 undergoes a dramatic transcriptional switch relocating from differentiation-specific to inflammatory genes thereby promoting AP-1-dependent gene transcription. Pancreatic deletion of c-Jun rescues the pre-inflammatory phenotype, Nr5a2 binding to inflammatory gene promoters, and the defective regenerative response to damage. These findings support the notion that, in the pancreas, the same transcriptional networks involved in differentiation-specific functions suppress inflammatory programmes. These networks can be subverted to foster inflammation upon genetic or environmental constraints.

Genome-wide meta-analysis identifies five new susceptibility loci for pancreatic cancer

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.

Mendelian randomisation study of the relationship between vitamin D and risk of glioma

To examine for a causal relationship between vitamin D and glioma risk we performed an analysis of genetic variants associated with serum 25-hydroxyvitamin D (25(OH)D) levels using Mendelian randomisation (MR), an approach unaffected by biases from confounding. Two-sample MR was undertaken using genome-wide association study data. Single nucleotide polymorphisms (SNPs) associated with 25(OH)D levels were used as instrumental variables (IVs). We calculated MR estimates for the odds ratio (OR) for 25(OH)D levels with glioma using SNP-glioma estimates from 12,488 cases and 18,169 controls, using inverse-variance weighted (IVW) and maximum likelihood estimation (MLE) methods. A non-significant association between 25(OH)D levels and glioma risk was shown using both the IVW (OR = 1.21, 95% confidence interval [CI] = 0.90–1.62, P = 0.201) and MLE (OR = 1.20, 95% CI = 0.98–1.48, P = 0.083) methods. In an exploratory analysis of tumour subtype, an inverse relationship between 25(OH)D levels and glioblastoma (GBM) risk was identified using the MLE method (OR = 0.62, 95% CI = 0.43–0.89, P = 0.010), but not the IVW method (OR = 0.62, 95% CI = 0.37–1.04, P = 0.070). No statistically significant association was shown between 25(OH)D levels and non-GBM glioma. Our results do not provide evidence for a causal relationship between 25(OH)D levels and all forms of glioma risk. More evidence is required to explore the relationship between 25(OH)D levels and risk of GBM.

The oral microbiota in patients with pancreatic cancer, patients with IPMNs, and controls: a pilot study

PURPOSE

Poor oral health appears to be a risk factor for pancreatic cancer, possibly implicating the oral microbiota. In this pilot study, we evaluated the characteristics of the oral microbiota in patients with pancreatic ductal adenocarcinoma (PDAC), intraductal papillary mucinous neoplasms (IPMN), and healthy controls.

METHODS

Forty newly diagnosed PDAC patients, 39 IPMN patients, and 58 controls, excluding current smokers and users of antibiotics, provided saliva samples. Common oral bacterial species were comprehensively surveyed by sequencing of the 16S rRNA microbial genes. We obtained measures of diversity and the mean relative proportions of individual taxa. We explored the degree to which these measures differed according to respondent characteristics based on individual interviews.

RESULTS

PDAC cases did not differ in diversity measures from either controls or IPMN cases. PDAC cases had higher mean relative proportions of Firmicutes and related taxa, while controls had higher mean relative proportions of Proteobacteria and related taxa. Results were generally similar when comparing PDAC to IPMN cases. Among IPMNs and controls combined, younger individuals had higher levels of several taxa within the Proteobacteria. The only other variable consistently related to mean relative proportions was mouthwash use, with taxa within Firmicutes more common among users.

CONCLUSIONS

While there were no differences in diversity of the oral microbiota among these groups, there were differences in the mean relative proportions of some taxa. Characteristics of the oral microbiota are not associated with most measures of oral health.

Impact of Sixteen Established Pancreatic Cancer Susceptibility Loci in American Jews

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.

Abstract 1315: Estimating sex-specific effects of genetic loci associated with glioma risk

BACKGROUND: Incidence of glioma varies significantly by sex, and most glioma histologies occur with greater incidence in males. Previous analyses have examined the impact of estrogen exposure as a risk factor for these tumors, but have found results of varying significance and low effect size. There may be differences in effect of previously discovered risk alleles that contribute to sex differences.

METHODS: Using data collected for three previous glioma GWAS in European-ancestry populations (MD Anderson Cancer Center, the San Francisco Adult Glioma Study, and the Glioma International Case Control Study) we assessed sex-specific effects for 14 previously identified and 13 newly identified glioma risk SNPs (27 total) overall and for glioblastoma (GBM) and non-GBM tumors separately. There were 3,892 male cases (59% GBM), 4,522 male controls, 2,500 female cases (52% GBM) and 4,940 female controls. Sex-specific odds ratios (ORM and ORF), 95% confidence intervals (95% CI) and p values (pM and pF) were generated using stratified logistic regression models. Data from each study were analyzed separately and combined using inverse variance weighted meta-analysis. Results were considered statistically significant at p&lt;6.2x10-4.

RESULTS: In GBM, rs11979158 (7p11.2, pM=1.01x10-10, ORM=1.43 [95% CI: 1.28-1.59]; pF=3.43x10-3, ORF=1.22 [95% CI: 1.07-1.39]) and rs2562152 (16p13.3, pM=5.59x10-4, ORM=1.23 [95% CI: 1.09-1.39]; pF=1.22x10-1, ORF=1.12 [95% CI: 0.97-1.29]) had significant effect in males only. In non-GBM gliomas, rs12076373 (1q44, pM=3.56x10-7, ORM=1.41, 95% CI: 1.23-1.61; pF=1.76x10-2, ORF=1.20, 95% CI: 1.03-1.38), rs11979158 (7p11.2, pM=1.77x10-5, ORM=1.32 [95% CI: 1.16-1.49]; pF=2.73x10-1, ORF=1.08 [95% CI: 0.94-1.24]), and rs3751667 (16p13.3, pM=9.44x10-7, ORM=1.30 [95% CI: 1.17-1.44]; pF=5.18x10-2, ORF=1.13 [95% CI: 1.00-1.28]) had significant effect in males only. Effect size for rs55705857 (8q24.21) varied significantly by sex, with ORM=2.63 (95% CI: 2.24-3.09, pM=4.42x10-32), as compared to ORF=3.95 (95% CI: 3.28-4.76, pF=1.82x10-47). A sensitivity analysis was performed due to allele frequency heterogeneity by study and results did not change.

CONCLUSIONS: Sex differences and other demographic differences in cancer susceptibility can provide important clues to etiology, and these differences can be leveraged for discovery in genetic association studies. Significant differences in effect size may suggest variation in genetic effect of risk alleles or in unidentified risk factors that vary in prevalence or effect by sex. There may also be differences in the distribution of molecular subtypes within each histology by sex. Further investigation using an agnostic approach may further elucidate the relationship between effect of risk alleles and sex.

Citation Format: Quinn T. Ostrom, Ben Kinnersley, Margaret Wrensch, Jeanette E. Eckel-Passow, Georgina Armstrong, Terri Rice, Yanwen Chen, John Wiencke, Lucie McCoy, Helen Hansen, Christopher Amos, Jonine L. Bernstein, Elizabeth B. Claus, Dora Il'yasova, Christoffer Johansen, Daniel Lachance, Rose Lai, Ching C. Lau, Ryan T. Merrell, Sara H. Olson, Siegal Sadetzki, Joellen Schildkraut, Sanjay Shete, Richard S. Houlston, Robert B. Jenkins, Beatrice Melin, Melissa Bondy, Jill S. Barnholtz-Sloan. Estimating sex-specific effects of genetic loci associated with glioma risk [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1315. doi:10.1158/1538-7445.AM2017-1315

Abstract 1302: Genome-wide association study of glioma reveals specific differences in genetic susceptibility to glioblastoma and non-glioblastoma

BACKGROUND: Glioma accounts for ~27% of all primary brain tumors and is responsible for ~13,000 cancer-related deaths in the US each year. Glioma tumors can be broadly classified into glioblastoma (GBM) and lower-grade non-GBM. Typically gliomas have a poor prognosis irrespective of medical care, with the most common form, GBM, having a five-year survival rate of only 5%. While genome-wide association studies (GWAS) have transformed our understanding of glioma susceptibility, individual studies have had limited power to identify risk loci.

METHODS: We performed the largest glioma GWAS to date, comprising a meta-analysis of six existing GWAS (6,405 cases, 14,100 controls) as well as new GWAS from the Glioma International Case Control Consortium (GICC; 4,572 cases and 3,286 controls) and University of California, San Francisco (UCSF)-Mayo (1,519 cases, 804 controls), totaling 12,496 cases (6,191 classified as GBM, 5,819 as non-GBM) and 18,190 controls.

RESULTS: We identified five new risk loci for GBM at 1p31.3 (rs12752552; near JAK1, P=2.04×10-9, odds ratio (OR)=1.22), 11q14.1 (rs11233250; P=9.95×10-10, OR=1.24), 16p13.3 (rs2562152; near MPG, P=1.93x10-8, OR=1.21), 16q12.1 (rs10852606; HEATR3, P=1.29×10-11, OR=1.18), 22q13.1 (rs2235573; P=1.76×10-10, OR=1.15) and eight for non-GBM at 1q32.1 (rs4252707; MDM4, P=3.34×10-9, OR=1.19), 1q44 (rs12076373; AKT3, P=2.63×10-10, OR=1.23), 2q33.3 (rs7572263; near IDH1, P=2.18×10-10, OR=1.20), 3p14.1 (rs11706832; LRIG1, P=7.66×10-9, OR=1.15), 10q24.33 (rs11598018; OBFC1, P=3.39×10-8, OR=1.14), 11q21 (rs7107785; P=3.87×10-10, OR=1.16), 14q12 (rs10131032; P=5.07x10-11, OR=1.33) and 16p13.3 (rs3751667; P=2.61×10-9, OR=1.18). Case-only analyses confirmed the specificity of 11q14.1, 16p13.3 and 22q13.1 associations for GBM and 1q44, 2q33.3, 3p14.1, 11q21 and 14q12 for non-GBM tumors. In the combined meta-analysis, among previously published glioma risk SNPs, those for all glioma at 17p13.1 (TP53), GBM at 5p15.33 (TERT), 7p11.2 (EGFR), 9p21.3 (CDKN2B-AS1) and 20q13.33 (RTEL1) and for non-GBM at 8q24.21 (CCDC26), 11q23.2, 11q23.3 (PHLDB1) and 15q24.2 (ETFA) showed even greater evidence for association. SNPs at 10q25.2 and 12q12.1 for non-GBM tumors retained genome-wide significance (i.e. P&lt;5.0x10-8). Associations at the previously reported loci for GBM at 3q26.2 (near TERC) and 12q23.33 (POLR3B) did not retain statistical significance.

CONCLUSIONS: Our findings substantiate genetic susceptibility to GBM and non-GBM glioma being highly distinct, consistent with their distinctive molecular profiles presumably resulting from different etiological pathways. Functional analyses should lead to further insights into the biological basis of the different glioma histologies. Such information can inform gene discovery initiatives and therefore have a measurable impact on the successful development of new therapeutic agents.

Citation Format: Ben Kinnersley, Beatrice S. Melin, Jill S. Barnholtz-Sloan, Margaret R. Wrensch, Christoffer Johansen, Dora Il’yasova, Quinn Ostrom, and members of GICC, Karim Labreche, Jeanette E. Eckel-Passow, Paul A. Decker, Marianne Labussière, Ahmed Idbaih, Khe Hoang-Xuan, Anna-Luisa Di Stefano, Karima Mokhtari, Jean-Yves Delattre, Pilar Galan, Konstantinos Gousias, Johannes Schramm, Minouk J. Schoemaker, Sarah J. Fleming, Stefan Herms, Stefanie Heilmann, Marcus M. Nöthen, Heinz-Erich Wichmann, Stefan Schreiber, Anthony Swerdlow, Mark Lathrop, Matthias Simon, Marc Sanson, Preetha Rajaraman, Stephen Chanock, Martha Linet, Zhaoming Wang, Meredith Yeager, Rose K. Lai, Elizabeth B. Claus, Sara H. Olson, Robert B. Jenkins, Richard S. Houlston, Melissa L. Bondy. Genome-wide association study of glioma reveals specific differences in genetic susceptibility to glioblastoma and non-glioblastoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1302. doi:10.1158/1538-7445.AM2017-1302

Abstract 1326: A pooled genome-wide association study of pancreatic cancer susceptibility loci in American Jews

Jews are estimated to be at increased risk of pancreatic cancer compared to non-Jews. The higher risk of pancreatic cancer in Jews can be partially explained by the increased frequency of BRCA1 and BRCA2 mutations in the Ashkenazi Jewish population; however, the remaining 40-70% excess risk of pancreatic cancer in Jews is not explained by other established non-genetic and genetic risk factors. The genetic origins of modern Ashkenazi Jews can be traced back to a founder population bottleneck between the 11th and 15th centuries follow by a rapid endogamous population expansion, making the Ashkenazi Jewish population a genetic isolate ideal for exploring genetic contributions to pancreatic cancer. We conducted a genome-wide association study (GWAS) in a case-control sample of American Jews, largely Ashkenazi, including 406 pancreatic cancer patients and 2,332 controls, identified in the database of Genotypes and Phenotypes (dbGaP) Pancreatic Cancer Cohort and Case-Control Consortium I/II (PanScan I/II), Pancreatic Cancer Case-Control Consortium (PanC4), and Genetic Epidemiology Research on Adult Health and Aging (GERA) data sets. We then examined those single nucleotide polymorphisms (SNPs) with P&lt;10-7 in an expanded sample set, of 539 full- plus part-Jewish pancreatic cancer patients and 4,117 full- plus part-Jewish controls from the same data sets. Jewish ancestry was genetically determined using seeded fast principal component analysis. In our data set, the part-Jewish subjects had mixed genetic ancestry both from Jews and from non-Jewish white-Europeans. Among the full Jews, a novel genome-wide significant association was detected on chromosome 19p12 (rs66562280, OR=1.55, 95% CI=1.33-1.81, P=2.30x10-8). A suggestive association was detected on chromosome 19p13.3 (rs2656937, OR=1.53, 95% CI=1.31-1.78, P=9.96x10-8). Similar associations were seen for these SNPs among the full- plus part-Jews. No SNPs meeting the quality control filter inclusion criteria were found in the full- plus part-Jewish subjects to be significant at the P=10-7 level. Overall, we identified in Jews one novel susceptibility locus and one suggestive novel susceptibility locus for pancreatic cancer that warrant follow-up work. This is the first GWAS conducted for pancreatic cancer in the increased-risk Jewish population. The novel susceptibility locus discovered on chromosome 19p12 could explain up to ~15%, and the novel suggestive susceptibility locus on chromosome 19p13.3 could explain up to an additional ~15% of the increased risk for pancreatic cancer in the Jewish population.

Citation Format: Samantha A. Streicher, Alison P. Klein, Sara H. Olson, Robert C. Kurtz, Andrew T. DeWan, Hongyu Zhao, Harvey A. Risch. A pooled genome-wide association study of pancreatic cancer susceptibility loci in American Jews [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1326. doi:10.1158/1538-7445.AM2017-1326

Epidemiology of pancreatic adenocarcinoma

Pancreatic ductal adenocarcinoma is one of the most lethal cancers worldwide. The highest incidence rates her found are in North America and in Western Europe while lower rates in Asian Africa, with age standard incidence rates of 7.2 and 2.8 per 100,000 populations. Unfortunately the vast majority of individuals with pancreatic cancer present with symptoms, and once symptoms develop the chance for surgery is only about 20%. Additionally he incidence rate and mortality from pancreatic ductal adenocarcinoma in the United States shows a very close association suggesting that her earlier detection and treatment does little to change the outcome from this disease. Multiple of environmental and genetic risk factors have been implicated in the development of pancreatic ductal adenocarcinoma. We have reviewed those risk factors we believe are most important the development of this lethal disease. It is hoped that in the future, identification of biomarkers unique in the development of pancreatic ductal adenocarcinoma will be identified leading to earlier detection and a greater frequency of potential cure of this disease.

Proportion of cancer in a Middle eastern country attributable to established risk factors

Background

Providing an estimate of the percentage of cancer in Lebanon by 2018 that is due to the exposure to risk factors in 2008. Factors include: smoking, body mass index (BMI), physical inactivity, dietary factors, alcohol consumption, infections, and air pollution in adults.

Method

Population Attributable Fraction (PAF) was calculated using the proportion of the population exposed and relative risks for each risk factor from meta-analyses. The PAF estimates the proportion of cases in which exposure may have played a causal role.

Results

Smoking caused most cancer cases, and it will further add a total of 1800 new cases by 2018. Among many other cancers, lung cancer had the largest proportion attributable of around 75%. BMI is expected to increase colorectal, liver and gastric cardia carcinoma specifically in males. High physical activity has a an average of 15% protection rate on cancer on colorectal cancer. Minimal adherence to Mediterranean diet will affect gastric cancer incidence by 7%. Cases of oropharyngeal and esophageal cancer will be the result of alcohol consumption mainly in males. H.Pylori infection is expected to result in half of the gastric cases by 2018. The high exposure to air pollution is expected to contribute by 13% to lung cancer cases in 2018.

Conclusion

The highest benefits can be achieved by controlling tobacco smoking. Interrelated and small changes in weight, physical activity and healthy diet with limited alcohol consumption can protect against several GI cancers in the long run. These results can be used to determine public health interventions that target important risk factors in the general population.

Genome-wide association study of glioma subtypes identifies specific differences in genetic susceptibility to glioblastoma and non-glioblastoma tumors

Genome-wide association studies (GWAS) have transformed our understanding of glioma susceptibility, but individual studies have had limited power to identify risk loci. We performed a meta-analysis of existing GWAS and two new GWAS, which totaled 12,496 cases and 18,190 controls. We identified five new loci for glioblastoma (GBM) at 1p31.3 (rs12752552; P = 2.04 × 10-9, odds ratio (OR) = 1.22), 11q14.1 (rs11233250; P = 9.95 × 10-10, OR = 1.24), 16p13.3 (rs2562152; P = 1.93 × 10-8, OR = 1.21), 16q12.1 (rs10852606; P = 1.29 × 10-11, OR = 1.18) and 22q13.1 (rs2235573; P = 1.76 × 10-10, OR = 1.15), as well as eight loci for non-GBM tumors at 1q32.1 (rs4252707; P = 3.34 × 10-9, OR = 1.19), 1q44 (rs12076373; P = 2.63 × 10-10, OR = 1.23), 2q33.3 (rs7572263; P = 2.18 × 10-10, OR = 1.20), 3p14.1 (rs11706832; P = 7.66 × 10-9, OR = 1.15), 10q24.33 (rs11598018; P = 3.39 × 10-8, OR = 1.14), 11q21 (rs7107785; P = 3.87 × 10-10, OR = 1.16), 14q12 (rs10131032; P = 5.07 × 10-11, OR = 1.33) and 16p13.3 (rs3751667; P = 2.61 × 10-9, OR = 1.18). These data substantiate that genetic susceptibility to GBM and non-GBM tumors are highly distinct, which likely reflects different etiology.

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.

Pelvic Inflammatory Disease and the Risk of Ovarian Cancer and Borderline Ovarian Tumors: A Pooled Analysis of 13 Case-Control Studies

Inflammation has been implicated in ovarian carcinogenesis. However, studies investigating the association between pelvic inflammatory disease (PID) and ovarian cancer risk are few and inconsistent. We investigated the association between PID and the risk of epithelial ovarian cancer according to tumor behavior and histotype. We pooled data from 13 case-control studies, conducted between 1989 and 2009, from the Ovarian Cancer Association Consortium (OCAC), including 9,162 women with ovarian cancers, 2,354 women with borderline tumors, and 14,736 control participants. Study-specific odds ratios were estimated and subsequently combined into a pooled odds ratio using a random-effects model. A history of PID was associated with an increased risk of borderline tumors (pooled odds ratio (pOR) = 1.32, 95% confidence interval (CI): 1.10, 1.58). Women with at least 2 episodes of PID had a 2-fold increased risk of borderline tumors (pOR = 2.14, 95% CI: 1.08, 4.24). No association was observed between PID and ovarian cancer risk overall (pOR = 0.99, 95% CI: 0.83, 1.19); however, a statistically nonsignificantly increased risk of low-grade serous tumors (pOR = 1.48, 95% CI: 0.92, 2.38) was noted. In conclusion, PID was associated with an increased risk of borderline ovarian tumors, particularly among women who had had multiple episodes of PID. Although our results indicated a histotype-specific association with PID, the association of PID with ovarian cancer risk is still somewhat uncertain and requires further investigation.

A targeted genetic association study of epithelial ovarian cancer susceptibility

BACKGROUND

Genome-wide association studies have identified several common susceptibility alleles for epithelial ovarian cancer (EOC). To further understand EOC susceptibility, we examined previously ungenotyped candidate variants, including uncommon variants and those residing within known susceptibility loci.

RESULTS

At nine of eleven previously published EOC susceptibility regions (2q31, 3q25, 5p15, 8q21, 8q24, 10p12, 17q12, 17q21.31, and 19p13), novel variants were identified that were more strongly associated with risk than previously reported variants. Beyond known susceptibility regions, no variants were found to be associated with EOC risk at genome-wide statistical significance (p <5x10(-8)), nor were any significant after Bonferroni correction for 17,000 variants (p< 3x10-6).

METHODS

A customized genotyping array was used to assess over 17,000 variants in coding, non-coding, regulatory, and known susceptibility regions in 4,973 EOC cases and 5,640 controls from 13 independent studies. Susceptibility for EOC overall and for select histotypes was evaluated using logistic regression adjusted for age, study site, and population substructure.

CONCLUSION

Given the novel variants identified within the 2q31, 3q25, 5p15, 8q21, 8q24, 10p12, 17q12, 17q21.31, and 19p13 regions, larger follow-up genotyping studies, using imputation where necessary, are needed for fine-mapping and confirmation of low frequency variants that fall below statistical significance.

PALB2, CHEK2 and ATM rare variants and cancer risk: data from COGS

BACKGROUND

The rarity of mutations in PALB2, CHEK2 and ATM make it difficult to estimate precisely associated cancer risks. Population-based family studies have provided evidence that at least some of these mutations are associated with breast cancer risk as high as those associated with rare BRCA2 mutations. We aimed to estimate the relative risks associated with specific rare variants in PALB2, CHEK2 and ATM via a multicentre case-control study.

METHODS

We genotyped 10 rare mutations using the custom iCOGS array: PALB2 c.1592delT, c.2816T>G and c.3113G>A, CHEK2 c.349A>G, c.538C>T, c.715G>A, c.1036C>T, c.1312G>T, and c.1343T>G and ATM c.7271T>G. We assessed associations with breast cancer risk (42 671 cases and 42 164 controls), as well as prostate (22 301 cases and 22 320 controls) and ovarian (14 542 cases and 23 491 controls) cancer risk, for each variant.

RESULTS

For European women, strong evidence of association with breast cancer risk was observed for PALB2 c.1592delT OR 3.44 (95% CI 1.39 to 8.52, p=7.1×10-5), PALB2 c.3113G>A OR 4.21 (95% CI 1.84 to 9.60, p=6.9×10-8) and ATM c.7271T>G OR 11.0 (95% CI 1.42 to 85.7, p=0.0012). We also found evidence of association with breast cancer risk for three variants in CHEK2, c.349A>G OR 2.26 (95% CI 1.29 to 3.95), c.1036C>T OR 5.06 (95% CI 1.09 to 23.5) and c.538C>T OR 1.33 (95% CI 1.05 to 1.67) (p≤0.017). Evidence for prostate cancer risk was observed for CHEK2 c.1343T>G OR 3.03 (95% CI 1.53 to 6.03, p=0.0006) for African men and CHEK2 c.1312G>T OR 2.21 (95% CI 1.06 to 4.63, p=0.030) for European men. No evidence of association with ovarian cancer was found for any of these variants.

CONCLUSIONS

This report adds to accumulating evidence that at least some variants in these genes are associated with an increased risk of breast cancer that is clinically important.

Inherited variants affecting RNA editing may contribute to ovarian cancer susceptibility: results from a large-scale collaboration

RNA editing in mammals is a form of post-transcriptional modification in which adenosine is converted to inosine by the adenosine deaminases acting on RNA (ADAR) family of enzymes. Based on evidence of altered ADAR expression in epithelial ovarian cancers (EOC), we hypothesized that single nucleotide polymorphisms (SNPs) in ADAR genes modify EOC susceptibility, potentially by altering ovarian tissue gene expression. Using directly genotyped and imputed data from 10,891 invasive EOC cases and 21,693 controls, we evaluated the associations of 5,303 SNPs in ADAD1, ADAR, ADAR2, ADAR3, and SND1. Unconditional logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI), with adjustment for European ancestry. We conducted gene-level analyses using the Admixture Maximum Likelihood (AML) test and the Sequence-Kernel Association test for common and rare variants (SKAT-CR). Association analysis revealed top risk-associated SNP rs77027562 (OR (95% CI)= 1.39 (1.17-1.64), P=1.0x10-4) in ADAR3 and rs185455523 in SND1 (OR (95% CI)= 0.68 (0.56-0.83), P=2.0x10-4). When restricting to serous histology (n=6,500), the magnitude of association strengthened for rs185455523 (OR=0.60, P=1.0x10-4). Gene-level analyses revealed that variation in ADAR was associated (P<0.05) with EOC susceptibility, with PAML=0.022 and PSKAT-CR=0.020. Expression quantitative trait locus analysis in EOC tissue revealed significant associations (P<0.05) with ADAR expression for several SNPs in ADAR, including rs1127313 (G/A), a SNP in the 3' untranslated region. In summary, germline variation involving RNA editing genes may influence EOC susceptibility, warranting further investigation of inherited and acquired alterations affecting RNA editing.

Functional characterization of a chr13q22.1 pancreatic cancer risk locus reveals long-range interaction and allele-specific effects on DIS3 expression

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 (r2 = 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.

Menstrual and Reproductive Factors, Hormone Use, and Risk of Pancreatic Cancer: Analysis From the International Pancreatic Cancer Case-Control Consortium (PanC4)

OBJECTIVES

We aimed to evaluate the relation between menstrual and reproductive factors, exogenous hormones, and risk of pancreatic cancer (PC).

METHODS

Eleven case-control studies within the International Pancreatic Cancer Case-control Consortium took part in the present study, including in total 2838 case and 4748 control women. Pooled estimates of odds ratios (ORs) and their 95% confidence intervals (CIs) were calculated using a 2-step logistic regression model and adjusting for relevant covariates.

RESULTS

An inverse OR was observed in women who reported having had hysterectomy (ORyesvs.no, 0.78; 95% CI, 0.67-0.91), remaining significant in postmenopausal women and never-smoking women, adjusted for potential PC confounders. A mutually adjusted model with the joint effect for hormone replacement therapy (HRT) and hysterectomy showed significant inverse associations with PC in women who reported having had hysterectomy with HRT use (OR, 0.64; 95% CI, 0.48-0.84).

CONCLUSIONS

Our large pooled analysis suggests that women who have had a hysterectomy may have reduced risk of PC. However, we cannot rule out that the reduced risk could be due to factors or indications for having had a hysterectomy. Further investigation of risk according to HRT use and reason for hysterectomy may be necessary.

Assessment of variation in immunosuppressive pathway genes reveals TGFBR2 to be associated with risk of clear cell ovarian cancer

BACKGROUND

Regulatory T (Treg) cells, a subset of CD4+ T lymphocytes, are mediators of immunosuppression in cancer, and, thus, variants in genes encoding Treg cell immune molecules could be associated with ovarian cancer.

METHODS

In a population of 15,596 epithelial ovarian cancer (EOC) cases and 23,236 controls, we measured genetic associations of 1,351 SNPs in Treg cell pathway genes with odds of ovarian cancer and tested pathway and gene-level associations, overall and by histotype, for the 25 genes, using the admixture likelihood (AML) method. The most significant single SNP associations were tested for correlation with expression levels in 44 ovarian cancer patients.

RESULTS

The most significant global associations for all genes in the pathway were seen in endometrioid ( p = 0.082) and clear cell ( p = 0.083), with the most significant gene level association seen with TGFBR2 ( p = 0.001) and clear cell EOC. Gene associations with histotypes at p < 0.05 included: IL12 ( p = 0.005 and p = 0.008, serous and high-grade serous, respectively), IL8RA ( p = 0.035, endometrioid and mucinous), LGALS1 ( p = 0.03, mucinous), STAT5B ( p = 0.022, clear cell), TGFBR1 ( p = 0.021 endometrioid) and TGFBR2 ( p = 0.017 and p = 0.025, endometrioid and mucinous, respectively).

CONCLUSIONS

Common inherited gene variation in Treg cell pathways shows some evidence of germline genetic contribution to odds of EOC that varies by histologic subtype and may be associated with mRNA expression of immune-complex receptor in EOC patients.

Risk Prediction for Epithelial Ovarian Cancer in 11 United States-Based Case-Control Studies: Incorporation of Epidemiologic Risk Factors and 17 Confirmed Genetic Loci

Previously developed models for predicting absolute risk of invasive epithelial ovarian cancer have included a limited number of risk factors and have had low discriminatory power (area under the receiver operating characteristic curve (AUC) < 0.60). Because of this, we developed and internally validated a relative risk prediction model that incorporates 17 established epidemiologic risk factors and 17 genome-wide significant single nucleotide polymorphisms (SNPs) using data from 11 case-control studies in the United States (5,793 cases; 9,512 controls) from the Ovarian Cancer Association Consortium (data accrued from 1992 to 2010). We developed a hierarchical logistic regression model for predicting case-control status that included imputation of missing data. We randomly divided the data into an 80% training sample and used the remaining 20% for model evaluation. The AUC for the full model was 0.664. A reduced model without SNPs performed similarly (AUC = 0.649). Both models performed better than a baseline model that included age and study site only (AUC = 0.563). The best predictive power was obtained in the full model among women younger than 50 years of age (AUC = 0.714); however, the addition of SNPs increased the AUC the most for women older than 50 years of age (AUC = 0.638 vs. 0.616). Adapting this improved model to estimate absolute risk and evaluating it in prospective data sets is warranted.