Genome-wide association analyses of ovarian cancer patients undergoing primary debulking surgery identify candidate genes for residual disease

Survival from ovarian cancer depends on the resection status after primary surgery. We performed genome-wide association analyses for resection status of 7705 ovarian cancer patients, including 4954 with high-grade serous carcinoma (HGSOC), to identify variants associated with residual disease. The most significant association with resection status was observed for rs72845444, upstream of MGMT, in HGSOC (p = 3.9 × 10-8). In gene-based analyses, PPP2R5C was the most strongly associated gene in HGSOC after stage adjustment. In an independent set of 378 ovarian tumours from the AGO-OVAR 11 study, variants near MGMT and PPP2R5C correlated with methylation and transcript levels, and PPP2R5C mRNA levels predicted progression-free survival in patients with residual disease. MGMT encodes a DNA repair enzyme, and PPP2R5C encodes the B56γ subunit of the PP2A tumour suppressor. Our results link heritable variation at these two loci with resection status in HGSOC.

Large-scale genome-wide association study of 398,238 women unveils seven novel loci associated with high-grade serous epithelial ovarian cancer risk

Background

Nineteen genomic regions have been associated with high-grade serous ovarian cancer (HGSOC). We used data from the Ovarian Cancer Association Consortium (OCAC), Consortium of Investigators of Modifiers of BRCA1/BRCA2 (CIMBA), UK Biobank (UKBB), and FinnGen to identify novel HGSOC susceptibility loci and develop polygenic scores (PGS).

Methods

We analyzed >22 million variants for 398,238 women. Associations were assessed separately by consortium and meta-analysed. OCAC and CIMBA data were used to develop PGS which were trained on FinnGen data and validated in UKBB and BioBank Japan.

Results

Eight novel variants were associated with HGSOC risk. An interesting discovery biologically was finding that TP53 3'-UTR SNP rs78378222 was associated with HGSOC (per T allele relative risk (RR)=1.44, 95%CI:1.28-1.62, P=1.76×10-9). The optimal PGS included 64,518 variants and was associated with an odds ratio of 1.46 (95%CI:1.37-1.54) per standard deviation in the UKBB validation (AUROC curve=0.61, 95%CI:0.59-0.62).

Conclusions

This study represents the largest GWAS for HGSOC to date. The results highlight that improvements in imputation reference panels and increased sample sizes can identify HGSOC associated variants that previously went undetected, resulting in improved PGS. The use of updated PGS in cancer risk prediction algorithms will then improve personalized risk prediction for HGSOC.

Exome sequencing identifies breast cancer susceptibility genes and defines the contribution of coding variants to breast cancer risk

Linkage and candidate gene studies have identified several breast cancer susceptibility genes, but the overall contribution of coding variation to breast cancer is unclear. To evaluate the role of rare coding variants more comprehensively, we performed a meta-analysis across three large whole-exome sequencing datasets, containing 26,368 female cases and 217,673 female controls. Burden tests were performed for protein-truncating and rare missense variants in 15,616 and 18,601 genes, respectively. Associations between protein-truncating variants and breast cancer were identified for the following six genes at exome-wide significance (P < 2.5 × 10-6): the five known susceptibility genes ATM, BRCA1, BRCA2, CHEK2 and PALB2, together with MAP3K1. Associations were also observed for LZTR1, ATR and BARD1 with P < 1 × 10-4. Associations between predicted deleterious rare missense or protein-truncating variants and breast cancer were additionally identified for CDKN2A at exome-wide significance. The overall contribution of coding variants in genes beyond the previously known genes is estimated to be small.

Folate Intake and Ovarian Cancer Risk among Women with Endometriosis: A Case-Control Study from the Ovarian Cancer Association Consortium

BACKGROUND

Although folate intake has not been associated with an increased risk of ovarian cancer overall, studies of other cancer types have suggested that high folate intake may promote carcinogenesis in precancerous lesions. Women with endometriosis (a potential precancerous lesion) have an increased risk of developing ovarian cancer; however, whether high folate intake increases risk in this group is unknown.

METHODS

We conducted a pooled analysis of six case-control studies from the Ovarian Cancer Association Consortium to investigate the association between folate intake and risk of ovarian cancer among women with and without self-reported endometriosis. We included 570 cases/558 controls with and 5,171/7,559 without endometriosis. We used logistic regression to estimate odds ratios (OR) and 95% confidence intervals for the association between folate intake (dietary, supplemental, and total) and ovarian cancer risk. Finally, we used Mendelian randomization (MR) to evaluate our results using genetic markers as a proxy for folate status.

RESULTS

Higher dietary folate intake was associated with an increased risk of ovarian cancer for women with endometriosis [OR, 1.37 (1.01-1.86)] but not for women without endometriosis. There was no association between supplemental folate intake and ovarian cancer risk for women with or without endometriosis. A similar pattern was seen using MR.

CONCLUSIONS

High dietary folate intake may be associated with an increased risk of ovarian cancer among women with endometriosis.

IMPACT

Women with endometriosis with high folate diets may be at increased risk of ovarian cancer. Further research is needed on the potential cancer-promoting effects of folate in this group.

Copy Number Variants Are Ovarian Cancer Risk Alleles at Known and Novel Risk Loci

BACKGROUND

Known risk alleles for epithelial ovarian cancer (EOC) account for approximately 40% of the heritability for EOC. Copy number variants (CNVs) have not been investigated as EOC risk alleles in a large population cohort.

METHODS

Single nucleotide polymorphism array data from 13 071 EOC cases and 17 306 controls of White European ancestry were used to identify CNVs associated with EOC risk using a rare admixture maximum likelihood test for gene burden and a by-probe ratio test. We performed enrichment analysis of CNVs at known EOC risk loci and functional biofeatures in ovarian cancer-related cell types.

RESULTS

We identified statistically significant risk associations with CNVs at known EOC risk genes; BRCA1 (PEOC = 1.60E-21; OREOC = 8.24), RAD51C (Phigh-grade serous ovarian cancer [HGSOC] = 5.5E-4; odds ratio [OR]HGSOC = 5.74 del), and BRCA2 (PHGSOC = 7.0E-4; ORHGSOC = 3.31 deletion). Four suggestive associations (P < .001) were identified for rare CNVs. Risk-associated CNVs were enriched (P < .05) at known EOC risk loci identified by genome-wide association study. Noncoding CNVs were enriched in active promoters and insulators in EOC-related cell types.

CONCLUSIONS

CNVs in BRCA1 have been previously reported in smaller studies, but their observed frequency in this large population-based cohort, along with the CNVs observed at BRCA2 and RAD51C gene loci in EOC cases, suggests that these CNVs are potentially pathogenic and may contribute to the spectrum of disease-causing mutations in these genes. CNVs are likely to occur in a wider set of susceptibility regions, with potential implications for clinical genetic testing and disease prevention.

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.

Rare germline copy number variants (CNVs) and breast cancer risk

Germline copy number variants (CNVs) are pervasive in the human genome but potential disease associations with rare CNVs have not been comprehensively assessed in large datasets. We analysed rare CNVs in genes and non-coding regions for 86,788 breast cancer cases and 76,122 controls of European ancestry with genome-wide array data. Gene burden tests detected the strongest association for deletions in BRCA1 (P = 3.7E-18). Nine other genes were associated with a p-value < 0.01 including known susceptibility genes CHEK2 (P = 0.0008), ATM (P = 0.002) and BRCA2 (P = 0.008). Outside the known genes we detected associations with p-values < 0.001 for either overall or subtype-specific breast cancer at nine deletion regions and four duplication regions. Three of the deletion regions were in established common susceptibility loci. To the best of our knowledge, this is the first genome-wide analysis of rare CNVs in a large breast cancer case-control dataset. We detected associations with exonic deletions in established breast cancer susceptibility genes. We also detected suggestive associations with non-coding CNVs in known and novel loci with large effects sizes. Larger sample sizes will be required to reach robust levels of statistical significance.

Identification of a Locus Near ULK1 Associated With Progression-Free Survival in Ovarian Cancer

BACKGROUND

Many loci have been found to be associated with risk of epithelial ovarian cancer (EOC). However, although there is considerable variation in progression-free survival (PFS), no loci have been found to be associated with outcome at genome-wide levels of significance.

METHODS

We carried out a genome-wide association study (GWAS) of PFS in 2,352 women with EOC who had undergone cytoreductive surgery and standard carboplatin/paclitaxel chemotherapy.

RESULTS

We found seven SNPs at 12q24.33 associated with PFS (P < 5 × 10-8), the top SNP being rs10794418 (HR = 1.24; 95% CI, 1.15-1.34; P = 1.47 × 10-8). High expression of a nearby gene, ULK1, is associated with shorter PFS in EOC, and with poor prognosis in other cancers. SNP rs10794418 is also associated with expression of ULK1 in ovarian tumors, with the allele associated with shorter PFS being associated with higher expression, and chromatin interactions were detected between the ULK1 promoter and associated SNPs in serous and endometrioid EOC cell lines. ULK1 knockout ovarian cancer cell lines showed significantly increased sensitivity to carboplatin in vitro.

CONCLUSIONS

The locus at 12q24.33 represents one of the first genome-wide significant loci for survival for any cancer. ULK1 is a plausible candidate for the target of this association.

IMPACT

This finding provides insight into genetic markers associated with EOC outcome and potential treatment options.See related commentary by Peres and Monteiro, p. 1604.

Pleiotropy-guided transcriptome imputation from normal and tumor tissues identifies candidate susceptibility genes for breast and ovarian cancer

Familial, sequencing, and genome-wide association studies (GWASs) and genetic correlation analyses have progressively unraveled the shared or pleiotropic germline genetics of breast and ovarian cancer. In this study, we aimed to leverage this shared germline genetics to improve the power of transcriptome-wide association studies (TWASs) to identify candidate breast cancer and ovarian cancer susceptibility genes. We built gene expression prediction models using the PrediXcan method in 681 breast and 295 ovarian tumors from The Cancer Genome Atlas and 211 breast and 99 ovarian normal tissue samples from the Genotype-Tissue Expression project and integrated these with GWAS meta-analysis data from the Breast Cancer Association Consortium (122,977 cases/105,974 controls) and the Ovarian Cancer Association Consortium (22,406 cases/40,941 controls). The integration was achieved through application of a pleiotropy-guided conditional/conjunction false discovery rate (FDR) approach in the setting of a TWASs. This identified 14 candidate breast cancer susceptibility genes spanning 11 genomic regions and 8 candidate ovarian cancer susceptibility genes spanning 5 genomic regions at conjunction FDR < 0.05 that were >1 Mb away from known breast and/or ovarian cancer susceptibility loci. We also identified 38 candidate breast cancer susceptibility genes and 17 candidate ovarian cancer susceptibility genes at conjunction FDR < 0.05 at known breast and/or ovarian susceptibility loci. The 22 genes identified by our cross-cancer analysis represent promising candidates that further elucidate the role of the transcriptome in mediating germline breast and ovarian cancer risk.

Author Correction: A case-only study to identify genetic modifiers of breast cancer risk for BRCA1/BRCA2 mutation carriers

A Correction to this paper has been published: https://doi.org/10.1038/s41467-021-23162-4

A case-only study to identify genetic modifiers of breast cancer risk for BRCA1/BRCA2 mutation carriers

Breast cancer (BC) risk for BRCA1 and BRCA2 mutation carriers varies by genetic and familial factors. About 50 common variants have been shown to modify BC risk for mutation carriers. All but three, were identified in general population studies. Other mutation carrier-specific susceptibility variants may exist but studies of mutation carriers have so far been underpowered. We conduct a novel case-only genome-wide association study comparing genotype frequencies between 60,212 general population BC cases and 13,007 cases with BRCA1 or BRCA2 mutations. We identify robust novel associations for 2 variants with BC for BRCA1 and 3 for BRCA2 mutation carriers, P < 10-8, at 5 loci, which are not associated with risk in the general population. They include rs60882887 at 11p11.2 where MADD, SP11 and EIF1, genes previously implicated in BC biology, are predicted as potential targets. These findings will contribute towards customising BC polygenic risk scores for BRCA1 and BRCA2 mutation carriers.

Mendelian randomization analyses suggest a role for cholesterol in the development of endometrial cancer

Blood lipids have been associated with the development of a range of cancers, including breast, lung and colorectal cancer. For endometrial cancer, observational studies have reported inconsistent associations between blood lipids and cancer risk. To reduce biases from unmeasured confounding, we performed a bidirectional, two-sample Mendelian randomization analysis to investigate the relationship between levels of three blood lipids (low-density lipoprotein [LDL] and high-density lipoprotein [HDL] cholesterol, and triglycerides) and endometrial cancer risk. Genetic variants associated with each of these blood lipid levels (P < 5 × 10-8 ) were identified as instrumental variables, and assessed using genome-wide association study data from the Endometrial Cancer Association Consortium (12 906 cases and 108 979 controls) and the Global Lipids Genetic Consortium (n = 188 578). Mendelian randomization analyses found genetically raised LDL cholesterol levels to be associated with lower risks of endometrial cancer of all histologies combined, and of endometrioid and non-endometrioid subtypes. Conversely, higher genetically predicted HDL cholesterol levels were associated with increased risk of non-endometrioid endometrial cancer. After accounting for the potential confounding role of obesity (as measured by genetic variants associated with body mass index), the association between genetically predicted increased LDL cholesterol levels and lower endometrial cancer risk remained significant, especially for non-endometrioid endometrial cancer. There was no evidence to support a role for triglycerides in endometrial cancer development. Our study supports a role for LDL and HDL cholesterol in the development of non-endometrioid endometrial cancer. Further studies are required to understand the mechanisms underlying these findings.

European polygenic risk score for prediction of breast cancer shows similar performance in Asian women

Polygenic risk scores (PRS) have been shown to predict breast cancer risk in European women, but their utility in Asian women is unclear. Here we evaluate the best performing PRSs for European-ancestry women using data from 17,262 breast cancer cases and 17,695 controls of Asian ancestry from 13 case-control studies, and 10,255 Chinese women from a prospective cohort (413 incident breast cancers). Compared to women in the middle quintile of the risk distribution, women in the highest 1% of PRS distribution have a ~2.7-fold risk and women in the lowest 1% of PRS distribution has ~0.4-fold risk of developing breast cancer. There is no evidence of heterogeneity in PRS performance in Chinese, Malay and Indian women. A PRS developed for European-ancestry women is also predictive of breast cancer risk in Asian women and can help in developing risk-stratified screening programmes in Asia.

Fine-mapping of 150 breast cancer risk regions identifies 191 likely target genes

Genome-wide association studies have identified breast cancer risk variants in over 150 genomic regions, but the mechanisms underlying risk remain largely unknown. These regions were explored by combining association analysis with in silico genomic feature annotations. We defined 205 independent risk-associated signals with the set of credible causal variants in each one. In parallel, we used a Bayesian approach (PAINTOR) that combines genetic association, linkage disequilibrium and enriched genomic features to determine variants with high posterior probabilities of being causal. Potentially causal variants were significantly over-represented in active gene regulatory regions and transcription factor binding sites. We applied our INQUSIT pipeline for prioritizing genes as targets of those potentially causal variants, using gene expression (expression quantitative trait loci), chromatin interaction and functional annotations. Known cancer drivers, transcription factors and genes in the developmental, apoptosis, immune system and DNA integrity checkpoint gene ontology pathways were over-represented among the highest-confidence target genes.

Publisher Correction: Shared heritability and functional enrichment across six solid cancers

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

An amendment to this paper has been published and can be accessed via a link at the top of the 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.

Genetic Data from Nearly 63,000 Women of European Descent Predicts DNA Methylation Biomarkers and Epithelial Ovarian Cancer Risk

DNA methylation is instrumental for gene regulation. Global changes in the epigenetic landscape have been recognized as a hallmark of cancer. However, the role of DNA methylation in epithelial ovarian cancer (EOC) remains unclear. In this study, high-density genetic and DNA methylation data in white blood cells from the Framingham Heart Study (N = 1,595) were used to build genetic models to predict DNA methylation levels. These prediction models were then applied to the summary statistics of a genome-wide association study (GWAS) of ovarian cancer including 22,406 EOC cases and 40,941 controls to investigate genetically predicted DNA methylation levels in association with EOC risk. Among 62,938 CpG sites investigated, genetically predicted methylation levels at 89 CpG were significantly associated with EOC risk at a Bonferroni-corrected threshold of P < 7.94 × 10-7. Of them, 87 were located at GWAS-identified EOC susceptibility regions and two resided in a genomic region not previously reported to be associated with EOC risk. Integrative analyses of genetic, methylation, and gene expression data identified consistent directions of associations across 12 CpG, five genes, and EOC risk, suggesting that methylation at these 12 CpG may influence EOC risk by regulating expression of these five genes, namely MAPT, HOXB3, ABHD8, ARHGAP27, and SKAP1. We identified novel DNA methylation markers associated with EOC risk and propose that methylation at multiple CpG may affect EOC risk via regulation of gene expression. SIGNIFICANCE: Identification of novel DNA methylation markers associated with EOC risk suggests that methylation at multiple CpG may affect EOC risk through regulation of gene expression.

Functional Analysis and Fine Mapping of the 9p22.2 Ovarian Cancer Susceptibility Locus

Genome-wide association studies have identified 40 ovarian cancer risk loci. However, the mechanisms underlying these associations remain elusive. In this study, we conducted a two-pronged approach to identify candidate causal SNPs and assess underlying biological mechanisms at chromosome 9p22.2, the first and most statistically significant associated locus for ovarian cancer susceptibility. Three transcriptional regulatory elements with allele-specific effects and a scaffold/matrix attachment region were characterized and, through physical DNA interactions, BNC2 was established as the most likely target gene. We determined the consensus binding sequence for BNC2 in vitro, verified its enrichment in BNC2 ChIP-seq regions, and validated a set of its downstream target genes. Fine-mapping by dense regional genotyping in over 15,000 ovarian cancer cases and 30,000 controls identified SNPs in the scaffold/matrix attachment region as among the most likely causal variants. This study reveals a comprehensive regulatory landscape at 9p22.2 and proposes a likely mechanism of susceptibility to ovarian cancer. SIGNIFICANCE: Mapping the 9p22.2 ovarian cancer risk locus identifies BNC2 as an ovarian cancer risk gene.See related commentary by Choi and Brown, p. 439.

Shared heritability and functional enrichment across six solid cancers

Quantifying the genetic correlation between cancers can provide important insights into the mechanisms driving cancer etiology. Using genome-wide association study summary statistics across six cancer types based on a total of 296,215 cases and 301,319 controls of European ancestry, here we estimate the pair-wise genetic correlations between breast, colorectal, head/neck, lung, ovary and prostate cancer, and between cancers and 38 other diseases. We observed statistically significant genetic correlations between lung and head/neck cancer (rg = 0.57, p = 4.6 × 10-8), breast and ovarian cancer (rg = 0.24, p = 7 × 10-5), breast and lung cancer (rg = 0.18, p =1.5 × 10-6) and breast and colorectal cancer (rg = 0.15, p = 1.1 × 10-4). We also found that multiple cancers are genetically correlated with non-cancer traits including smoking, psychiatric diseases and metabolic characteristics. Functional enrichment analysis revealed a significant excess contribution of conserved and regulatory regions to cancer heritability. Our comprehensive analysis of cross-cancer heritability suggests that solid tumors arising across tissues share in part a common germline genetic basis.

Polygenic Risk Scores for Prediction of Breast Cancer and Breast Cancer Subtypes

Stratification of women according to their risk of breast cancer based on polygenic risk scores (PRSs) could improve screening and prevention strategies. Our aim was to develop PRSs, optimized for prediction of estrogen receptor (ER)-specific disease, from the largest available genome-wide association dataset and to empirically validate the PRSs in prospective studies. The development dataset comprised 94,075 case subjects and 75,017 control subjects of European ancestry from 69 studies, divided into training and validation sets. Samples were genotyped using genome-wide arrays, and single-nucleotide polymorphisms (SNPs) were selected by stepwise regression or lasso penalized regression. The best performing PRSs were validated in an independent test set comprising 11,428 case subjects and 18,323 control subjects from 10 prospective studies and 190,040 women from UK Biobank (3,215 incident breast cancers). For the best PRSs (313 SNPs), the odds ratio for overall disease per 1 standard deviation in ten prospective studies was 1.61 (95%CI: 1.57-1.65) with area under receiver-operator curve (AUC) = 0.630 (95%CI: 0.628-0.651). The lifetime risk of overall breast cancer in the top centile of the PRSs was 32.6%. Compared with women in the middle quintile, those in the highest 1% of risk had 4.37- and 2.78-fold risks, and those in the lowest 1% of risk had 0.16- and 0.27-fold risks, of developing ER-positive and ER-negative disease, respectively. Goodness-of-fit tests indicated that this PRS was well calibrated and predicts disease risk accurately in the tails of the distribution. This PRS is a powerful and reliable predictor of breast cancer risk that may improve breast cancer prevention programs.

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.

Identification of nine new susceptibility loci for endometrial cancer

Endometrial cancer is the most commonly diagnosed cancer of the female reproductive tract in developed countries. Through genome-wide association studies (GWAS), we have previously identified eight risk loci for endometrial cancer. Here, we present an expanded meta-analysis of 12,906 endometrial cancer cases and 108,979 controls (including new genotype data for 5624 cases) and identify nine novel genome-wide significant loci, including a locus on 12q24.12 previously identified by meta-GWAS of endometrial and colorectal cancer. At five loci, expression quantitative trait locus (eQTL) analyses identify candidate causal genes; risk alleles at two of these loci associate with decreased expression of genes, which encode negative regulators of oncogenic signal transduction proteins (SH2B3 (12q24.12) and NF1 (17q11.2)). In summary, this study has doubled the number of known endometrial cancer risk loci and revealed candidate causal genes for future study.

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.

Genetic overlap between endometriosis and endometrial cancer: evidence from cross-disease genetic correlation and GWAS meta-analyses

Epidemiological, biological, and molecular data suggest links between endometriosis and endometrial cancer, with recent epidemiological studies providing evidence for an association between a previous diagnosis of endometriosis and risk of endometrial cancer. We used genetic data as an alternative approach to investigate shared biological etiology of these two diseases. Genetic correlation analysis of summary level statistics from genomewide association studies (GWAS) using LD Score regression revealed moderate but significant genetic correlation (rg  = 0.23, P = 9.3 × 10-3 ), and SNP effect concordance analysis provided evidence for significant SNP pleiotropy (P = 6.0 × 10-3 ) and concordance in effect direction (P = 2.0 × 10-3 ) between the two diseases. Cross-disease GWAS meta-analysis highlighted 13 distinct loci associated at P ≤ 10-5 with both endometriosis and endometrial cancer, with one locus (SNP rs2475335) located within PTPRD associated at a genomewide significant level (P = 4.9 × 10-8 , OR = 1.11, 95% CI = 1.07-1.15). PTPRD acts in the STAT3 pathway, which has been implicated in both endometriosis and endometrial cancer. This study demonstrates the value of cross-disease genetic analysis to support epidemiological observations and to identify biological pathways of relevance to multiple diseases.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

In the version of this article originally published, one of the two authors with the name Wei Zhao was omitted from the author list and the affiliations for both authors were assigned to the single Wei Zhao in the author list. In addition, the ORCID for Wei Zhao (Department of Biostatistics and Epidemiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA) was incorrectly assigned to author Wei Zhou. The errors have been corrected in the HTML and PDF versions of the article.

Publisher Correction: Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

In the published version of this paper, the name of author Emanuele Di Angelantonio was misspelled. This error has now been corrected in the HTML and PDF versions of the article.

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.

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.

Protein-altering variants associated with body mass index implicate pathways that control energy intake and expenditure in obesity

Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity.

Analyses of germline variants associated with ovarian cancer survival identify functional candidates at the 1q22 and 19p12 outcome loci

We previously identified associations with ovarian cancer outcome at five genetic loci. To identify putatively causal genetic variants and target genes, we prioritized two ovarian outcome loci (1q22 and 19p12) for further study. Bioinformatic and functional genetic analyses indicated that MEF2D and ZNF100 are targets of candidate outcome variants at 1q22 and 19p12, respectively. At 19p12, the chromatin interaction of a putative regulatory element with the ZNF100 promoter region correlated with candidate outcome variants. At 1q22, putative regulatory elements enhanced MEF2D promoter activity and haplotypes containing candidate outcome variants modulated these effects. In a public dataset, MEF2D and ZNF100 expression were both associated with ovarian cancer progression-free or overall survival time. In an extended set of 6,162 epithelial ovarian cancer patients, we found that functional candidates at the 1q22 and 19p12 loci, as well as other regional variants, were nominally associated with patient outcome; however, no associations reached our threshold for statistical significance (p<1×10-5). Larger patient numbers will be needed to convincingly identify any true associations at these loci.

Germline whole exome sequencing and large-scale replication identifies FANCM as a likely high grade serous ovarian cancer susceptibility gene

We analyzed whole exome sequencing data in germline DNA from 412 high grade serous ovarian cancer (HGSOC) cases from The Cancer Genome Atlas Project and identified 5,517 genes harboring a predicted deleterious germline coding mutation in at least one HGSOC case. Gene-set enrichment analysis showed enrichment for genes involved in DNA repair (p = 1.8×10-3). Twelve DNA repair genes - APEX1, APLF, ATX, EME1, FANCL, FANCM, MAD2L2, PARP2, PARP3, POLN, RAD54L and SMUG1 - were prioritized for targeted sequencing in up to 3,107 HGSOC cases, 1,491 cases of other epithelial ovarian cancer (EOC) subtypes and 3,368 unaffected controls of European origin. We estimated mutation prevalence for each gene and tested for associations with disease risk. Mutations were identified in both cases and controls in all genes except MAD2L2, where we found no evidence of mutations in controls. In FANCM we observed a higher mutation frequency in HGSOC cases compared to controls (29/3,107 cases, 0.96 percent; 13/3,368 controls, 0.38 percent; P=0.008) with little evidence for association with other subtypes (6/1,491, 0.40 percent; P=0.82). The relative risk of HGSOC associated with deleterious FANCM mutations was estimated to be 2.5 (95% CI 1.3 - 5.0; P=0.006). In summary, whole exome sequencing of EOC cases with large-scale replication in case-control studies has identified FANCM as a likely novel susceptibility gene for HGSOC, with mutations associated with a moderate increase in risk. These data may have clinical implications for risk prediction and prevention approaches for high-grade serous ovarian cancer in the future and a significant impact on reducing disease mortality.

Prediction of Breast and Prostate Cancer Risks in Male BRCA1 and BRCA2 Mutation Carriers Using Polygenic Risk Scores

Purpose BRCA1/2 mutations increase the risk of breast and prostate cancer in men. Common genetic variants modify cancer risks for female carriers of BRCA1/2 mutations. We investigated-for the first time to our knowledge-associations of common genetic variants with breast and prostate cancer risks for male carriers of BRCA1/ 2 mutations and implications for cancer risk prediction. Materials and Methods We genotyped 1,802 male carriers of BRCA1/2 mutations from the Consortium of Investigators of Modifiers of BRCA1/2 by using the custom Illumina OncoArray. We investigated the combined effects of established breast and prostate cancer susceptibility variants on cancer risks for male carriers of BRCA1/2 mutations by constructing weighted polygenic risk scores (PRSs) using published effect estimates as weights. Results In male carriers of BRCA1/2 mutations, PRS that was based on 88 female breast cancer susceptibility variants was associated with breast cancer risk (odds ratio per standard deviation of PRS, 1.36; 95% CI, 1.19 to 1.56; P = 8.6 × 10-6). Similarly, PRS that was based on 103 prostate cancer susceptibility variants was associated with prostate cancer risk (odds ratio per SD of PRS, 1.56; 95% CI, 1.35 to 1.81; P = 3.2 × 10-9). Large differences in absolute cancer risks were observed at the extremes of the PRS distribution. For example, prostate cancer risk by age 80 years at the 5th and 95th percentiles of the PRS varies from 7% to 26% for carriers of BRCA1 mutations and from 19% to 61% for carriers of BRCA2 mutations, respectively. Conclusion PRSs may provide informative cancer risk stratification for male carriers of BRCA1/2 mutations that might enable these men and their physicians to make informed decisions on the type and timing of breast and prostate cancer risk management.

Genomic analyses identify hundreds of variants associated with age at menarche and support a role for puberty timing in cancer risk

The timing of puberty is a highly polygenic childhood trait that is epidemiologically associated with various adult diseases. Using 1000 Genomes Project-imputed genotype data in up to ∼370,000 women, we identify 389 independent signals (P < 5 × 10-8) for age at menarche, a milestone in female pubertal development. In Icelandic data, these signals explain ∼7.4% of the population variance in age at menarche, corresponding to ∼25% of the estimated heritability. We implicate ∼250 genes via coding variation or associated expression, demonstrating significant enrichment in neural tissues. Rare variants near the imprinted genes MKRN3 and DLK1 were identified, exhibiting large effects when paternally inherited. Mendelian randomization analyses suggest causal inverse associations, independent of body mass index (BMI), between puberty timing and risks for breast and endometrial cancers in women and prostate cancer in men. In aggregate, our findings highlight the complexity of the genetic regulation of puberty timing and support causal links with cancer susceptibility.

Identification of 12 new susceptibility loci for different histotypes of epithelial ovarian cancer

To identify common alleles associated with different histotypes of epithelial ovarian cancer (EOC), we pooled data from multiple genome-wide genotyping projects totaling 25,509 EOC cases and 40,941 controls. We identified nine new susceptibility loci for different EOC histotypes: six for serous EOC histotypes (3q28, 4q32.3, 8q21.11, 10q24.33, 18q11.2 and 22q12.1), two for mucinous EOC (3q22.3 and 9q31.1) and one for endometrioid EOC (5q12.3). We then performed meta-analysis on the results for high-grade serous ovarian cancer with the results from analysis of 31,448 BRCA1 and BRCA2 mutation carriers, including 3,887 mutation carriers with EOC. This identified three additional susceptibility loci at 2q13, 8q24.1 and 12q24.31. Integrated analyses of genes and regulatory biofeatures at each locus predicted candidate susceptibility genes, including OBFC1, a new candidate susceptibility gene for low-grade and borderline serous EOC.

Rare and low-frequency coding variants alter human adult height

Height is a highly heritable, classic polygenic trait with approximately 700 common associated variants identified through genome-wide association studies so far. Here, we report 83 height-associated coding variants with lower minor-allele frequencies (in the range of 0.1-4.8%) and effects of up to 2 centimetres per allele (such as those in IHH, STC2, AR and CRISPLD2), greater than ten times the average effect of common variants. In functional follow-up studies, rare height-increasing alleles of STC2 (giving an increase of 1-2 centimetres per allele) compromised proteolytic inhibition of PAPP-A and increased cleavage of IGFBP-4 in vitro, resulting in higher bioavailability of insulin-like growth factors. These 83 height-associated variants overlap genes that are mutated in monogenic growth disorders and highlight new biological candidates (such as ADAMTS3, IL11RA and NOX4) and pathways (such as proteoglycan and glycosaminoglycan synthesis) involved in growth. Our results demonstrate that sufficiently large sample sizes can uncover rare and low-frequency variants of moderate-to-large effect associated with polygenic human phenotypes, and that these variants implicate relevant genes and pathways.

Genetic Risk Score Mendelian Randomization Shows that Obesity Measured as Body Mass Index, but not Waist:Hip Ratio, Is Causal for Endometrial Cancer

BACKGROUND

The strongest known risk factor for endometrial cancer is obesity. To determine whether SNPs associated with increased body mass index (BMI) or waist-hip ratio (WHR) are associated with endometrial cancer risk, independent of measured BMI, we investigated relationships between 77 BMI and 47 WHR SNPs and endometrial cancer in 6,609 cases and 37,926 country-matched controls.

METHODS

Logistic regression analysis and fixed effects meta-analysis were used to test for associations between endometrial cancer risk and (i) individual BMI or WHR SNPs, (ii) a combined weighted genetic risk score (wGRS) for BMI or WHR. Causality of BMI for endometrial cancer was assessed using Mendelian randomization, with BMIwGRS as instrumental variable.

RESULTS

The BMIwGRS was significantly associated with endometrial cancer risk (P = 3.4 × 10-17). Scaling the effect of the BMIwGRS on endometrial cancer risk by its effect on BMI, the endometrial cancer OR per 5 kg/m2 of genetically predicted BMI was 2.06 [95% confidence interval (CI), 1.89-2.21], larger than the observed effect of BMI on endometrial cancer risk (OR = 1.55; 95% CI, 1.44-1.68, per 5 kg/m2). The association attenuated but remained significant after adjusting for BMI (OR = 1.22; 95% CI, 1.10-1.39; P = 5.3 × 10-4). There was evidence of directional pleiotropy (P = 1.5 × 10-4). BMI SNP rs2075650 was associated with endometrial cancer at study-wide significance (P < 4.0 × 10-4), independent of BMI. Endometrial cancer was not significantly associated with individual WHR SNPs or the WHRwGRS.

CONCLUSIONS

BMI, but not WHR, is causally associated with endometrial cancer risk, with evidence that some BMI-associated SNPs alter endometrial cancer risk via mechanisms other than measurable BMI.

IMPACT

The causal association between BMI SNPs and endometrial cancer has possible implications for endometrial cancer risk modeling. Cancer Epidemiol Biomarkers Prev; 25(11); 1503-10. ©2016 AACR.

Association of vitamin D levels and risk of ovarian cancer: a Mendelian randomization study

BACKGROUND

In vitro and observational epidemiological studies suggest that vitamin D may play a role in cancer prevention. However, the relationship between vitamin D and ovarian cancer is uncertain, with observational studies generating conflicting findings. A potential limitation of observational studies is inadequate control of confounding. To overcome this problem, we used Mendelian randomization (MR) to evaluate the association between single nucleotide polymorphisms (SNPs) associated with circulating 25-hydroxyvitamin D [25(OH)D] concentration and risk of ovarian cancer.

METHODS

We employed SNPs with well-established associations with 25(OH)D concentration as instrumental variables for MR: rs7944926 (DHCR7), rs12794714 (CYP2R1) and rs2282679 (GC). We included 31 719 women of European ancestry (10 065 cases, 21 654 controls) from the Ovarian Cancer Association Consortium, who were genotyped using customized Illumina Infinium iSelect (iCOGS) arrays. A two-sample (summary data) MR approach was used and analyses were performed separately for all ovarian cancer (10 065 cases) and for high-grade serous ovarian cancer (4121 cases).

RESULTS

The odds ratio for epithelial ovarian cancer risk (10 065 cases) estimated by combining the individual SNP associations using inverse variance weighting was 1.27 (95% confidence interval: 1.06 to 1.51) per 20 nmol/L decrease in 25(OH)D concentration. The estimated odds ratio for high-grade serous epithelial ovarian cancer (4121 cases) was 1.54 (1.19, 2.01).

CONCLUSIONS

Genetically lowered 25-hydroxyvitamin D concentrations were associated with higher ovarian cancer susceptibility in Europeans. These findings suggest that increasing plasma vitamin D levels may reduce risk of ovarian cancer.

Cross-Cancer Genome-Wide Analysis of Lung, Ovary, Breast, Prostate, and Colorectal Cancer Reveals Novel Pleiotropic Associations

Identifying genetic variants with pleiotropic associations can uncover common pathways influencing multiple cancers. We took a two-stage approach to conduct genome-wide association studies for lung, ovary, breast, prostate, and colorectal cancer from the GAME-ON/GECCO Network (61,851 cases, 61,820 controls) to identify pleiotropic loci. Findings were replicated in independent association studies (55,789 cases, 330,490 controls). We identified a novel pleiotropic association at 1q22 involving breast and lung squamous cell carcinoma, with eQTL analysis showing an association with ADAM15/THBS3 gene expression in lung. We also identified a known breast cancer locus CASP8/ALS2CR12 associated with prostate cancer, a known cancer locus at CDKN2B-AS1 with different variants associated with lung adenocarcinoma and prostate cancer, and confirmed the associations of a breast BRCA2 locus with lung and serous ovarian cancer. This is the largest study to date examining pleiotropy across multiple cancer-associated loci, identifying common mechanisms of cancer development and progression. Cancer Res; 76(17); 5103-14. ©2016 AACR.

Meta-analysis of Genome Wide Association Studies Identifies Genetic Markers of Late Toxicity Following Radiotherapy for Prostate Cancer

Nearly 50% of cancer patients undergo radiotherapy. Late radiotherapy toxicity affects quality-of-life in long-term cancer survivors and risk of side-effects in a minority limits doses prescribed to the majority of patients. Development of a test predicting risk of toxicity could benefit many cancer patients. We aimed to meta-analyze individual level data from four genome-wide association studies from prostate cancer radiotherapy cohorts including 1564 men to identify genetic markers of toxicity. Prospectively assessed two-year toxicity endpoints (urinary frequency, decreased urine stream, rectal bleeding, overall toxicity) and single nucleotide polymorphism (SNP) associations were tested using multivariable regression, adjusting for clinical and patient-related risk factors. A fixed-effects meta-analysis identified two SNPs: rs17599026 on 5q31.2 with urinary frequency (odds ratio [OR] 3.12, 95% confidence interval [CI] 2.08-4.69, p-value 4.16×10(-8)) and rs7720298 on 5p15.2 with decreased urine stream (OR 2.71, 95% CI 1.90-3.86, p-value=3.21×10(-8)). These SNPs lie within genes that are expressed in tissues adversely affected by pelvic radiotherapy including bladder, kidney, rectum and small intestine. The results show that heterogeneous radiotherapy cohorts can be combined to identify new moderate-penetrance genetic variants associated with radiotherapy toxicity. The work provides a basis for larger collaborative efforts to identify enough variants for a future test involving polygenic risk profiling.

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.

Correction: The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study

[This corrects the article DOI: 10.1371/journal.pgen.1005378.].

Adult body mass index and risk of ovarian cancer by subtype: a Mendelian randomization study

BACKGROUND

Observational studies have reported a positive association between body mass index (BMI) and ovarian cancer risk. However, questions remain as to whether this represents a causal effect, or holds for all histological subtypes. The lack of association observed for serous cancers may, for instance, be due to disease-associated weight loss. Mendelian randomization (MR) uses genetic markers as proxies for risk factors to overcome limitations of observational studies. We used MR to elucidate the relationship between BMI and ovarian cancer, hypothesizing that genetically predicted BMI would be associated with increased risk of non-high grade serous ovarian cancers (non-HGSC) but not HGSC.

METHODS

We pooled data from 39 studies (14 047 cases, 23 003 controls) in the Ovarian Cancer Association Consortium. We constructed a weighted genetic risk score (GRS, partial F-statistic = 172), summing alleles at 87 single nucleotide polymorphisms previously associated with BMI, weighting by their published strength of association with BMI. Applying two-stage predictor-substitution MR, we used logistic regression to estimate study-specific odds ratios (OR) and 95% confidence intervals (CI) for the association between genetically predicted BMI and risk, and pooled these using random-effects meta-analysis.

RESULTS

Higher genetically predicted BMI was associated with increased risk of non-HGSC (pooled OR = 1.29, 95% CI 1.03-1.61 per 5 units BMI) but not HGSC (pooled OR = 1.06, 95% CI 0.88-1.27). Secondary analyses stratified by behaviour/subtype suggested that, consistent with observational data, the association was strongest for low-grade/borderline serous cancers (OR = 1.93, 95% CI 1.33-2.81).

CONCLUSIONS

Our data suggest that higher BMI increases risk of non-HGSC, but not the more common and aggressive HGSC subtype, confirming the observational evidence.

Five endometrial cancer risk loci identified through genome-wide association analysis

We conducted a meta-analysis of three endometrial cancer genome-wide association studies (GWAS) and two follow-up phases totaling 7,737 endometrial cancer cases and 37,144 controls of European ancestry. Genome-wide imputation and meta-analysis identified five new risk loci of genome-wide significance at likely regulatory regions on chromosomes 13q22.1 (rs11841589, near KLF5), 6q22.31 (rs13328298, in LOC643623 and near HEY2 and NCOA7), 8q24.21 (rs4733613, telomeric to MYC), 15q15.1 (rs937213, in EIF2AK4, near BMF) and 14q32.33 (rs2498796, in AKT1, near SIVA1). We also found a second independent 8q24.21 signal (rs17232730). Functional studies of the 13q22.1 locus showed that rs9600103 (pairwise r(2) = 0.98 with rs11841589) is located in a region of active chromatin that interacts with the KLF5 promoter region. The rs9600103[T] allele that is protective in endometrial cancer suppressed gene expression in vitro, suggesting that regulation of the expression of KLF5, a gene linked to uterine development, is implicated in tumorigenesis. These findings provide enhanced insight into the genetic and biological basis of endometrial cancer.

No clinical utility of KRAS variant rs61764370 for ovarian or breast cancer

OBJECTIVE

Clinical genetic testing is commercially available for rs61764370, an inherited variant residing in a KRAS 3' UTR microRNA binding site, based on suggested associations with increased ovarian and breast cancer risk as well as with survival time. However, prior studies, emphasizing particular subgroups, were relatively small. Therefore, we comprehensively evaluated ovarian and breast cancer risks as well as clinical outcome associated with rs61764370.

METHODS

Centralized genotyping and analysis were performed for 140,012 women enrolled in the Ovarian Cancer Association Consortium (15,357 ovarian cancer patients; 30,816 controls), the Breast Cancer Association Consortium (33,530 breast cancer patients; 37,640 controls), and the Consortium of Modifiers of BRCA1 and BRCA2 (14,765 BRCA1 and 7904 BRCA2 mutation carriers).

RESULTS

We found no association with risk of ovarian cancer (OR=0.99, 95% CI 0.94-1.04, p=0.74) or breast cancer (OR=0.98, 95% CI 0.94-1.01, p=0.19) and results were consistent among mutation carriers (BRCA1, ovarian cancer HR=1.09, 95% CI 0.97-1.23, p=0.14, breast cancer HR=1.04, 95% CI 0.97-1.12, p=0.27; BRCA2, ovarian cancer HR=0.89, 95% CI 0.71-1.13, p=0.34, breast cancer HR=1.06, 95% CI 0.94-1.19, p=0.35). Null results were also obtained for associations with overall survival following ovarian cancer (HR=0.94, 95% CI 0.83-1.07, p=0.38), breast cancer (HR=0.96, 95% CI 0.87-1.06, p=0.38), and all other previously-reported associations.

CONCLUSIONS

rs61764370 is not associated with risk of ovarian or breast cancer nor with clinical outcome for patients with these cancers. Therefore, genotyping this variant has no clinical utility related to the prediction or management of these cancers.

Assessment of Multifactor Gene-Environment Interactions and Ovarian Cancer Risk: Candidate Genes, Obesity, and Hormone-Related Risk Factors

BACKGROUND

Many epithelial ovarian cancer (EOC) risk factors relate to hormone exposure and elevated estrogen levels are associated with obesity in postmenopausal women. Therefore, we hypothesized that gene-environment interactions related to hormone-related risk factors could differ between obese and non-obese women.

METHODS

We considered interactions between 11,441 SNPs within 80 candidate genes related to hormone biosynthesis and metabolism and insulin-like growth factors with six hormone-related factors (oral contraceptive use, parity, endometriosis, tubal ligation, hormone replacement therapy, and estrogen use) and assessed whether these interactions differed between obese and non-obese women. Interactions were assessed using logistic regression models and data from 14 case-control studies (6,247 cases; 10,379 controls). Histotype-specific analyses were also completed.

RESULTS

SNPs in the following candidate genes showed notable interaction: IGF1R (rs41497346, estrogen plus progesterone hormone therapy, histology = all, P = 4.9 × 10(-6)) and ESR1 (rs12661437, endometriosis, histology = all, P = 1.5 × 10(-5)). The most notable obesity-gene-hormone risk factor interaction was within INSR (rs113759408, parity, histology = endometrioid, P = 8.8 × 10(-6)).

CONCLUSIONS

We have demonstrated the feasibility of assessing multifactor interactions in large genetic epidemiology studies. Follow-up studies are necessary to assess the robustness of our findings for ESR1, CYP11A1, IGF1R, CYP11B1, INSR, and IGFBP2 Future work is needed to develop powerful statistical methods able to detect these complex interactions.

IMPACT

Assessment of multifactor interaction is feasible, and, here, suggests that the relationship between genetic variants within candidate genes and hormone-related risk factors may vary EOC susceptibility. Cancer Epidemiol Biomarkers Prev; 25(5); 780-90. ©2016 AACR.

Investigation of Exomic Variants Associated with Overall Survival in Ovarian Cancer

BACKGROUND

While numerous susceptibility loci for epithelial ovarian cancer (EOC) have been identified, few associations have been reported with overall survival. In the absence of common prognostic genetic markers, we hypothesize that rare coding variants may be associated with overall EOC survival and assessed their contribution in two exome-based genotyping projects of the Ovarian Cancer Association Consortium (OCAC).

METHODS

The primary patient set (Set 1) included 14 independent EOC studies (4,293 patients) and 227,892 variants, and a secondary patient set (Set 2) included six additional EOC studies (1,744 patients) and 114,620 variants. Because power to detect rare variants individually is reduced, gene-level tests were conducted. Sets were analyzed separately at individual variants and by gene, and then combined with meta-analyses (73,203 variants and 13,163 genes overlapped).

RESULTS

No individual variant reached genome-wide statistical significance. A SNP previously implicated to be associated with EOC risk and, to a lesser extent, survival, rs8170, showed the strongest evidence of association with survival and similar effect size estimates across sets (Pmeta = 1.1E-6, HRSet1 = 1.17, HRSet2 = 1.14). Rare variants in ATG2B, an autophagy gene important for apoptosis, were significantly associated with survival after multiple testing correction (Pmeta = 1.1E-6; Pcorrected = 0.01).

CONCLUSIONS

Common variant rs8170 and rare variants in ATG2B may be associated with EOC overall survival, although further study is needed.

IMPACT

This study represents the first exome-wide association study of EOC survival to include rare variant analyses, and suggests that complementary single variant and gene-level analyses in large studies are needed to identify rare variants that warrant follow-up study. Cancer Epidemiol Biomarkers Prev; 25(3); 446-54. ©2016 AACR.

CYP19A1 fine-mapping and Mendelian randomization: estradiol is causal for endometrial cancer

Candidate gene studies have reported CYP19A1 variants to be associated with endometrial cancer and with estradiol (E2) concentrations. We analyzed 2937 single nucleotide polymorphisms (SNPs) in 6608 endometrial cancer cases and 37 925 controls and report the first genome wide-significant association between endometrial cancer and a CYP19A1 SNP (rs727479 in intron 2, P=4.8×10(-11)). SNP rs727479 was also among those most strongly associated with circulating E2 concentrations in 2767 post-menopausal controls (P=7.4×10(-8)). The observed endometrial cancer odds ratio per rs727479 A-allele (1.15, CI=1.11-1.21) is compatible with that predicted by the observed effect on E2 concentrations (1.09, CI=1.03-1.21), consistent with the hypothesis that endometrial cancer risk is driven by E2. From 28 candidate-causal SNPs, 12 co-located with three putative gene-regulatory elements and their risk alleles associated with higher CYP19A1 expression in bioinformatical analyses. For both phenotypes, the associations with rs727479 were stronger among women with a higher BMI (Pinteraction=0.034 and 0.066 respectively), suggesting a biologically plausible gene-environment interaction.