Robust Tests for Additive Gene-Environment Interaction in Case-Control Studies Using Gene-Environment Independence

There have been recent proposals advocating the use of additive gene-environment interaction instead of the widely used multiplicative scale, as a more relevant public health measure. Using gene-environment independence enhances statistical power for testing multiplicative interaction in case-control studies. However, under departure from this assumption, substantial bias in the estimates and inflated type I error in the corresponding tests can occur. In this paper, we extend the empirical Bayes (EB) approach previously developed for multiplicative interaction, which trades off between bias and efficiency in a data-adaptive way, to the additive scale. An EB estimator of the relative excess risk due to interaction is derived, and the corresponding Wald test is proposed with a general regression setting under a retrospective likelihood framework. We study the impact of gene-environment association on the resultant test with case-control data. Our simulation studies suggest that the EB approach uses the gene-environment independence assumption in a data-adaptive way and provides a gain in power compared with the standard logistic regression analysis and better control of type I error when compared with the analysis assuming gene-environment independence. We illustrate the methods with data from the Ovarian Cancer Association Consortium.Multiple funders listed on paper

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

Population-based targeted sequencing of 54 candidate genes identifies PALB2 as a susceptibility gene for high-grade serous ovarian cancer

Purpose: The known epithelial ovarian cancer (EOC) susceptibility genes account for less than 50% of the heritable risk of ovarian cancer suggesting that other susceptibility genes exist. The aim of this study was to evaluate the contribution to ovarian cancer susceptibility of rare deleterious germline variants in a set of candidate genes. Methods: We sequenced the coding region of 54 candidate genes in 6385 invasive EOC cases and 6115 controls of broad European ancestry. Genes with an increased frequency of putative deleterious variants in cases versus controls were further examined in an independent set of 14 135 EOC cases and 28 655 controls from the Ovarian Cancer Association Consortium and the UK Biobank. For each gene, we estimated the EOC risks and evaluated associations between germline variant status and clinical characteristics. Results: The ORs associated for high-grade serous ovarian cancer were 3.01 for PALB2 (95% CI 1.59 to 5.68; p=0.00068), 1.99 for POLK (95% CI 1.15 to 3.43; p=0.014) and 4.07 for SLX4 (95% CI 1.34 to 12.4; p=0.013). Deleterious mutations in FBXO10 were associated with a reduced risk of disease (OR 0.27, 95% CI 0.07 to 1.00, p=0.049). However, based on the Bayes false discovery probability, only the association for PALB2 in high-grade serous ovarian cancer is likely to represent a true positive. Conclusions: We have found strong evidence that carriers of PALB2 deleterious mutations are at increased risk of high-grade serous ovarian cancer. Whether the magnitude of risk is sufficiently high to warrant the inclusion of PALB2 in cancer gene panels for ovarian cancer risk testing is unclear; much larger sample sizes will be needed to provide sufficiently precise estimates for clinical counselling

Enrichment of putative PAX8 target genes at serous epithelial ovarian cancer susceptibility loci

Background: Genome-wide association studies (GWAS) have identified 18 loci associated with serous ovarian cancer (SOC) susceptibility but the biological mechanisms driving these findings remain poorly characterised. Germline cancer risk loci may be enriched for target genes of transcription factors (TFs) critical to somatic tumorigenesis. Methods: All 615 TF-target sets from the Molecular Signatures Database were evaluated using gene set enrichment analysis (GSEA) and three GWAS for SOC risk: discovery (2196 cases/4396 controls), replication (7035 cases/21 693 controls; independent from discovery), and combined (9627 cases/30 845 controls; including additional individuals). Results: The PAX8-target gene set was ranked 1/615 in the discovery (PGSEA<0.001; FDR=0.21), 7/615 in the replication (PGSEA=0.004; FDR=0.37), and 1/615 in the combined (PGSEA<0.001; FDR=0.21) studies. Adding other genes reported to interact with PAX8 in the literature to the PAX8-target set and applying an alternative to GSEA, interval enrichment, further confirmed this association (P=0.006). Fifteen of the 157 genes from this expanded PAX8 pathway were near eight loci associated with SOC risk at P<10−5 (including six with P<5 × 10−8). The pathway was also associated with differential gene expression after shRNA-mediated silencing of PAX8 in HeyA8 (PGSEA=0.025) and IGROV1 (PGSEA=0.004) SOC cells and several PAX8 targets near SOC risk loci demonstrated in vitro transcriptomic perturbation. Conclusions: Putative PAX8 target genes are enriched for common SOC risk variants. This finding from our agnostic evaluation is of particular interest given that PAX8 is well-established as a specific marker for the cell of origin of SOC

Assessing the genetic architecture of epithelial ovarian cancer histological subtypes.

Epithelial ovarian cancer (EOC) is one of the deadliest common cancers. The five most common types of disease are high-grade and low-grade serous, endometrioid, mucinous and clear cell carcinoma. Each of these subtypes present distinct molecular pathogeneses and sensitivities to treatments. Recent studies show that certain genetic variants confer susceptibility to all subtypes while other variants are subtype-specific. Here, we perform an extensive analysis of the genetic architecture of EOC subtypes. To this end, we used data of 10,014 invasive EOC patients and 21,233 controls from the Ovarian Cancer Association Consortium genotyped in the iCOGS array (211,155 SNPs). We estimate the array heritability (attributable to variants tagged on arrays) of each subtype and their genetic correlations. We also look for genetic overlaps with factors such as obesity, smoking behaviors, diabetes, age at menarche and height. We estimated the array heritabilities of high-grade serous disease ([Formula: see text] = 8.8 ± 1.1 %), endometrioid ([Formula: see text] = 3.2 ± 1.6 %), clear cell ([Formula: see text] = 6.7 ± 3.3 %) and all EOC ([Formula: see text] = 5.6 ± 0.6 %). Known associated loci contributed approximately 40 % of the total array heritability for each subtype. The contribution of each chromosome to the total heritability was not proportional to chromosome size. Through bivariate and cross-trait LD score regression, we found evidence of shared genetic backgrounds between the three high-grade subtypes: serous, endometrioid and undifferentiated. Finally, we found significant genetic correlations of all EOC with diabetes and obesity using a polygenic prediction approach.The Ovarian Cancer Association Consortium is supported by a grant from the Ovarian Cancer Research Fund thanks to donations by the family and friends of Kathryn Sladek Smith (PPD/RPCI.07). The Nurses’ Health Studies would like to thank the participants and staff of the Nurses' Health Study and Nurses' Health Study II for their valuable contributions as well as the following state cancer registries for their help: AL, AZ, AR, CA, CO, CT, DE, FL, GA, ID, IL, IN, IA, KY, LA, ME, MD, MA, MI, NE, NH, NJ, NY, NC, ND, OH, OK, OR, PA, RI, SC, TN, TX, VA, WA, WY. The authors assume full responsibility for analyses and interpretation of these data. Funding of the constituent studies was provided by the California Cancer Research Program (00-01389V-20170, N01-CN25403, 2II0200); the Canadian Institutes of Health Research (MOP-86727); Cancer Australia; Cancer Council Victoria; Cancer Council Queensland; Cancer Council New South Wales; Cancer Council South Australia; Cancer Council Tasmania; Cancer Foundation of Western Australia; the Cancer Institute of New Jersey; Cancer Research UK (C490/A6187, C490/A10119, C490/A10124); the Danish Cancer Society (94-222-52); the ELAN Program of the University of Erlangen-Nuremberg; the Eve Appeal; the Helsinki University Central Hospital Research Fund; Helse Vest; the Norwegian Cancer Society; the Norwegian Research Council; the Ovarian Cancer Research Fund; Nationaal Kankerplan of Belgium; the L & S Milken Foundation; the Polish Ministry of Science and Higher Education (4 PO5C 028 14, 2 PO5A 068 27); the Roswell Park Cancer Institute Alliance Foundation; the US National Cancer Institute (K07-CA095666, K07-CA80668, K07-CA143047, K22-CA138563, N01-CN55424, N01-PC67001, N01-PC067010, N01-PC035137, P01-CA017054, P01-CA087696, P30-CA072720, P30-CA15083, P30-CA008748, P50-CA159981, P50-CA105009, P50-CA136393, R01-CA149429, R01-CA014089, R01-CA016056, R01-CA017054, R01-CA049449, R01-CA050385, R01-CA054419, R01-CA058598, R01-CA058860, R01-CA061107, R01-CA061132, R01-CA063678, R01-CA063682, R01-CA067262, R01-CA071766, R01-CA074850, R01-CA080978, R01-CA083918, R01-CA087538, R01-CA092044, R01-CA095023, R01-CA122443, R01-CA112523, R01-CA114343, R01-CA126841, R01-CA136924, R03-CA113148, R03-CA115195, U01-CA069417, U01-CA071966, UM1-CA186107, UM1-CA176726 and Intramural research funds); the NIH/National Center for Research Resources/General Clinical Research Center (MO1-RR000056); the US Army Medical Research and Material Command (DAMD17-01-1-0729, DAMD17-02-1-0666, DAMD17-02-1-0669, W81XWH-07-0449, W81XWH-10-1-02802); the US Public Health Service (PSA-042205); the National Health and Medical Research Council of Australia (199600 and 400281); the German Federal Ministry of Education and Research of Germany Programme of Clinical Biomedical Research (01GB 9401); the State of Baden-Wurttemberg through Medical Faculty of the University of Ulm (P.685); the German Cancer Research Center; the Minnesota Ovarian Cancer Alliance; the Mayo Foundation; the Fred C. and Katherine B. Andersen Foundation; the Lon V. Smith Foundation (LVS-39420); the Oak Foundation; Eve Appeal; the OHSU Foundation; the Mermaid I project; the Rudolf-Bartling Foundation; the UK National Institute for Health Research Biomedical Research Centres at the University of Cambridge, Imperial College London, University College Hospital ‘Womens Health Theme’ and the Royal Marsden Hospital; and WorkSafeBC 14. Investigator-specific funding: G.C.P receives scholarship support from the University of Queensland and QIMR Berghofer. Y.L. was supported by the NHMRC Early Career Fellowship. G.C.T. is supported by the National Health and Medical Research Council. S.M. was supported by an ARC Future Fellowship

Shared genetics underlying epidemiological association between endometriosis and ovarian cancer

Epidemiological studies have demonstrated associations between endometriosis and certain histotypes of ovarian cancer, including clear cell, low-grade serous and endometrioid carcinomas. We aimed to determine whether the observed associations might be due to shared genetic aetiology. To address this, we used two endometriosis datasets genotyped on common arrays with full-genome coverage (3194 cases and 7060 controls) and a large ovarian cancer dataset genotyped on the customized Illumina Infinium iSelect (iCOGS) arrays (10 065 cases and 21 663 controls). Previous work has suggested that a large number of genetic variants contribute to endometriosis and ovarian cancer (all histotypes combined) susceptibility. Here, using the iCOGS data, we confirmed polygenic architecture for most histotypes of ovarian cancer. This led us to evaluate if the polygenic effects are shared across diseases. We found evidence for shared genetic risks between endometriosis and all histotypes of ovarian cancer, except for the intestinal mucinous type. Clear cell carcinoma showed the strongest genetic correlation with endometriosis (0.51, 95% CI = 0.18-0.84). Endometrioid and low-grade serous carcinomas had similar correlation coefficients (0.48, 95% CI = 0.07-0.89 and 0.40, 95% CI = 0.05-0.75, respectively). High-grade serous carcinoma, which often arises from the fallopian tubes, showed a weaker genetic correlation with endometriosis (0.25, 95% CI = 0.11-0.39), despite the absence of a known epidemiological association. These results suggest that the epidemiological association between endometriosis and ovarian adenocarcinoma may be attributable to shared genetic susceptibility loci.Other Research Uni

Evaluation of Candidate Stromal Epithelial Cross-Talk Genes Identifies Association between Risk of Serous Ovarian Cancer and TERT, a Cancer Susceptibility “Hot-Spot”

We hypothesized that variants in genes expressed as a consequence of interactions between ovarian cancer cells and the host micro-environment could contribute to cancer susceptibility. We therefore used a two-stage approach to evaluate common single nucleotide polymorphisms (SNPs) in 173 genes involved in stromal epithelial interactions in the Ovarian Cancer Association Consortium (OCAC). In the discovery stage, cases with epithelial ovarian cancer (n = 675) and controls (n = 1,162) were genotyped at 1,536 SNPs using an Illumina GoldenGate assay. Based on Positive Predictive Value estimates, three SNPs—PODXL rs1013368, ITGA6 rs13027811, and MMP3 rs522616—were selected for replication using TaqMan genotyping in up to 3,059 serous invasive cases and 8,905 controls from 16 OCAC case-control studies. An additional 18 SNPs with Pper-allele<0.05 in the discovery stage were selected for replication in a subset of five OCAC studies (n = 1,233 serous invasive cases; n = 3,364 controls). The discovery stage associations in PODXL, ITGA6, and MMP3 were attenuated in the larger replication set (adj. Pper-allele≥0.5). However genotypes at TERT rs7726159 were associated with ovarian cancer risk in the smaller, five-study replication study (Pper-allele = 0.03). Combined analysis of the discovery and replication sets for this TERT SNP showed an increased risk of serous ovarian cancer among non-Hispanic whites [adj. ORper-allele 1.14 (1.04–1.24) p = 0.003]. Our study adds to the growing evidence that, like the 8q24 locus, the telomerase reverse transcriptase locus at 5p15.33, is a general cancer susceptibility locus

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

Common variants at theCHEK2gene locus and risk of epithelial ovarian cancer

Genome-wide association studies have identified 20 genomic regions associated with risk of epithelial ovarian cancer (EOC), but many additional risk variants may exist. Here, we evaluated associations between common genetic variants [single nucleotide polymorphisms (SNPs) and indels] in DNA repair genes and EOC risk. We genotyped 2896 common variants at 143 gene loci in DNA samples from 15 397 patients with invasive EOC and controls. We found evidence of associations with EOC risk for variants at FANCA, EXO1, E2F4, E2F2, CREB5 and CHEK2 genes (P ≤ 0.001). The strongest risk association was for CHEK2 SNP rs17507066 with serous EOC (P = 4.74 x 10(-7)). Additional genotyping and imputation of genotypes from the 1000 genomes project identified a slightly more significant association for CHEK2 SNP rs6005807 (r (2) with rs17507066 = 0.84, odds ratio (OR) 1.17, 95% CI 1.11-1.24, P = 1.1×10(-7)). We identified 293 variants in the region with likelihood ratios of less than 1:100 for representing the causal variant. Functional annotation identified 25 candidate SNPs that alter transcription factor binding sites within regulatory elements active in EOC precursor tissues. In The Cancer Genome Atlas dataset, CHEK2 gene expression was significantly higher in primary EOCs compared to normal fallopian tube tissues (P = 3.72×10(-8)). We also identified an association between genotypes of the candidate causal SNP rs12166475 (r (2) = 0.99 with rs6005807) and CHEK2 expression (P = 2.70×10(-8)). These data suggest that common variants at 22q12.1 are associated with risk of serous EOC and CHEK2 as a plausible target susceptibility gene.Other Research Uni