METHODS AND FINDINGS: Genetic instruments to proxy 12 risk factors were constructed by identifying single nucleotide polymorphisms (SNPs) that were robustly (P < 5 × 10-8) and independently associated with each respective risk factor in previously reported genome-wide association studies. These risk factors included genetic liability to 3 factors (endometriosis, polycystic ovary syndrome, type 2 diabetes) scaled to reflect a 50% higher odds liability to disease. We obtained summary statistics for the association of these SNPs with risk of overall and histotype-specific invasive epithelial ovarian cancer (22,406 cases; 40,941 controls) and low malignant potential tumours (3,103 cases; 40,941 controls) from the Ovarian Cancer Association Consortium (OCAC). The OCAC dataset comprises 63 genotyping project/case-control sets with participants of European ancestry recruited from 14 countries (US, Australia, Belarus, Germany, Belgium, Denmark, Finland, Norway, Canada, Poland, UK, Spain, Netherlands, and Sweden). SNPs were combined into multi-allelic inverse-variance-weighted fixed or random effects models to generate effect estimates and 95% confidence intervals (CIs). Three complementary sensitivity analyses were performed to examine violations of MR assumptions: MR-Egger regression and weighted median and mode estimators. A Bonferroni-corrected P value threshold was used to establish strong evidence (P < 0.0042) and suggestive evidence (0.0042 < P < 0.05) for associations. In MR analyses, there was strong or suggestive evidence that 2 of the 12 risk factors were associated with invasive epithelial ovarian cancer and 8 of the 12 were associated with 1 or more invasive epithelial ovarian cancer histotypes. There was strong evidence that genetic liability to endometriosis was associated with an increased risk of invasive epithelial ovarian cancer (odds ratio [OR] per 50% higher odds liability: 1.10, 95% CI 1.06-1.15; P = 6.94 × 10-7) and suggestive evidence that lifetime smoking exposure was associated with an increased risk of invasive epithelial ovarian cancer (OR per unit increase in smoking score: 1.36, 95% CI 1.04-1.78; P = 0.02). In analyses examining histotypes and low malignant potential tumours, the strongest associations found were between height and clear cell carcinoma (OR per SD increase: 1.36, 95% CI 1.15-1.61; P = 0.0003); age at natural menopause and endometrioid carcinoma (OR per year later onset: 1.09, 95% CI 1.02-1.16; P = 0.007); and genetic liability to polycystic ovary syndrome and endometrioid carcinoma (OR per 50% higher odds liability: 0.89, 95% CI 0.82-0.96; P = 0.002). There was little evidence for an association of genetic liability to type 2 diabetes, parity, or circulating levels of 25-hydroxyvitamin D and sex hormone binding globulin with ovarian cancer or its subtypes. The primary limitations of this analysis include the modest statistical power for analyses of risk factors in relation to some less common ovarian cancer histotypes (low grade serous, mucinous, and clear cell carcinomas), the inability to directly examine the association of some ovarian cancer risk factors that did not have robust genetic variants available to serve as proxies (e.g., oral contraceptive use, hormone replacement therapy), and the assumption of linear relationships between risk factors and ovarian cancer risk.

METHODS: In total, DNA samples were obtained from 14,525 case subjects with invasive EOC and from 23,447 controls from 43 sites in the Ovarian Cancer Association Consortium (OCAC). Two hundred seventy nine SNPs, representing 131 genes, were genotyped using an Illumina Infinium iSelect BeadChip as part of the Collaborative Oncological Gene-environment Study (COGS). SNP analyses were conducted using unconditional logistic regression under a log-additive model, and the FDR q<0.2 was applied to adjust for multiple comparisons.

RESULTS: The most significant evidence of an association for all invasive cancers combined and for the serous subtype was observed for SNP rs17216603 in the iron transporter gene HEPH (invasive: OR = 0.85, P = 0.00026; serous: OR = 0.81, P = 0.00020); this SNP was also associated with the borderline/low malignant potential (LMP) tumors (P = 0.021). Other genes significantly associated with EOC histological subtypes (p<0.05) included the UGT1A (endometrioid), SLC25A45 (mucinous), SLC39A11 (low malignant potential), and SERPINA7 (clear cell carcinoma). In addition, 1785 SNPs in six genes (HEPH, MGST1, SERPINA, SLC25A45, SLC39A11 and UGT1A) were imputed from the 1000 Genomes Project and examined for association with INV EOC in white-European subjects. The most significant imputed SNP was rs117729793 in SLC39A11 (per allele, OR = 2.55, 95% CI = 1.5-4.35, p = 5.66x10-4).

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).

METHODS: We selected TF genes within 1 Mb of the top signal at the 12 genome-wide significant risk loci. Mutual information, a form of correlation, was used to build networks of genes strongly coexpressed with each selected TF gene in the unified microarray dataset of 489 serous EOC tumors from The Cancer Genome Atlas. Genes represented in this dataset were subsequently ranked using a gene-level test based on results for germline SNPs from a serous EOC GWAS meta-analysis (2,196 cases/4,396 controls).

RESULTS: Gene set enrichment analysis identified six networks centered on TF genes (HOXB2, HOXB5, HOXB6, HOXB7 at 17q21.32 and HOXD1, HOXD3 at 2q31) that were significantly enriched for genes from the risk-associated end of the ranked list (P < 0.05 and FDR < 0.05). These results were replicated (P < 0.05) using an independent association study (7,035 cases/21,693 controls). Genes underlying enrichment in the six networks were pooled into a combined network.

CONCLUSION: We identified a HOX-centric network associated with serous EOC risk containing several genes with known or emerging roles in serous EOC development.

METHODS: Using LDScore regression, we explored the genetic correlation between endometrial cancer and ovarian cancer. To identify loci associated with the risk of both cancers, we implemented a pipeline of statistical genetic analyses (i.e., inverse-variance meta-analysis, colocalization, and M-values) and performed analyses stratified by subtype. Candidate target genes were then prioritized using functional genomic data.

RESULTS: Genetic correlation analysis revealed significant genetic correlation between the two cancers (rG = 0.43, P = 2.66 × 10-5). We found seven loci associated with risk for both cancers (P Bonferroni < 2.4 × 10-9). In addition, four novel subgenome-wide regions at 7p22.2, 7q22.1, 9p12, and 11q13.3 were identified (P < 5 × 10-7). Promoter-associated HiChIP chromatin loops from immortalized endometrium and ovarian cell lines and expression quantitative trait loci data highlighted candidate target genes for further investigation.

CONCLUSIONS: Using cross-cancer GWAS meta-analysis, we have identified several joint endometrial and ovarian cancer risk loci and candidate target genes for future functional analysis.

METHODS: We evaluated the potential of gene-based aggregation in the Breast Cancer Association Consortium cohorts including 83,471 cases and 59,199 controls. Low-frequency variants were aggregated for individual genes' coding and regulatory regions. Association results in European ancestry samples were compared to single-marker association results in the same cohort. Gene-based associations were also combined in meta-analysis across individuals with European, Asian, African, and Latin American and Hispanic ancestry.