Concurrent RB1 loss and BRCA-deficiency predicts enhanced immunological response and long-term survival in tubo-ovarian high-grade serous carcinoma

Background

Somatic loss of the tumour suppressor RB1 is a common event in tubo-ovarian high-grade serous carcinoma (HGSC), which frequently co-occurs with alterations in homologous recombination DNA repair genes including BRCA1 and BRCA2 (BRCA). We examined whether tumour expression of RB1 was associated with survival across ovarian cancer histotypes (HGSC, endometrioid (ENOC), clear cell (CCOC), mucinous (MOC), low-grade serous carcinoma (LGSC)), and how co-occurrence of germline BRCA pathogenic variants and RB1 loss influences long-term survival in a large series of HGSC.

Patients and methods

RB1 protein expression patterns were classified by immunohistochemistry in epithelial ovarian carcinomas of 7436 patients from 20 studies participating in the Ovarian Tumor Tissue Analysis consortium and assessed for associations with overall survival (OS), accounting for patient age at diagnosis and FIGO stage. We examined RB1 expression and germline BRCA status in a subset of 1134 HGSC, and related genotype to survival, tumour infiltrating CD8+ lymphocyte counts and transcriptomic subtypes. Using CRISPR-Cas9, we deleted RB1 in HGSC cell lines with and without BRCA1 mutations to model co-loss with treatment response. We also performed genomic analyses on 126 primary HGSC to explore the molecular characteristics of concurrent homologous recombination deficiency and RB1 loss.

Results

RB1 protein loss was most frequent in HGSC (16.4%) and was highly correlated with RB1 mRNA expression. RB1 loss was associated with longer OS in HGSC (hazard ratio [HR] 0.74, 95% confidence interval [CI] 0.66-0.83, P = 6.8 ×10-7), but with poorer prognosis in ENOC (HR 2.17, 95% CI 1.17-4.03, P = 0.0140). Germline BRCA mutations and RB1 loss co-occurred in HGSC (P < 0.0001). Patients with both RB1 loss and germline BRCA mutations had a superior OS (HR 0.38, 95% CI 0.25-0.58, P = 5.2 ×10-6) compared to patients with either alteration alone, and their median OS was three times longer than non-carriers whose tumours retained RB1 expression (9.3 years vs. 3.1 years). Enhanced sensitivity to cisplatin (P < 0.01) and paclitaxel (P < 0.05) was seen in BRCA1 mutated cell lines with RB1 knockout. Among 126 patients with whole-genome and transcriptome sequence data, combined RB1 loss and genomic evidence of homologous recombination deficiency was correlated with transcriptional markers of enhanced interferon response, cell cycle deregulation, and reduced epithelial-mesenchymal transition in primary HGSC. CD8+ lymphocytes were most prevalent in BRCA-deficient HGSC with co-loss of RB1.

Conclusions

Co-occurrence of RB1 loss and BRCA mutation was associated with exceptionally long survival in patients with HGSC, potentially due to better treatment response and immune stimulation.

Lifetime ovulatory years and risk of epithelial ovarian cancer: a multinational pooled analysis

BACKGROUND

The role of ovulation in epithelial ovarian cancer (EOC) is supported by the consistent protective effects of parity and oral contraceptive use. Whether these factors protect through anovulation alone remains unclear. We explored the association between lifetime ovulatory years (LOY) and EOC.

METHODS

LOY was calculated using 12 algorithms. Odds ratios (ORs) and 95% confidence intervals (CIs) estimated the association between LOY or LOY components and EOC among 26 204 control participants and 21 267 case patients from 25 studies. To assess whether LOY components act through ovulation suppression alone, we compared beta coefficients obtained from regression models with expected estimates assuming 1 year of ovulation suppression has the same effect regardless of source.

RESULTS

LOY was associated with increased EOC risk (OR per year increase = 1.014, 95% CI = 1.009 to 1.020 to OR per year increase = 1.044, 95% CI = 1.041 to 1.048). Individual LOY components, except age at menarche, also associated with EOC. The estimated model coefficient for oral contraceptive use and pregnancies were 4.45 times and 12- to 15-fold greater than expected, respectively. LOY was associated with high-grade serous, low-grade serous, endometrioid, and clear cell histotypes (ORs per year increase = 1.054, 1.040, 1.065, and 1.098, respectively) but not mucinous tumors. Estimated coefficients of LOY components were close to expected estimates for high-grade serous but larger than expected for low-grade serous, endometrioid, and clear cell histotypes.

CONCLUSIONS

LOY is positively associated with nonmucinous EOC. Differences between estimated and expected model coefficients for LOY components suggest factors beyond ovulation underlie the associations between LOY components and EOC in general and for non-HGSOC.

p53 and ovarian carcinoma survival: an Ovarian Tumor Tissue Analysis consortium study

Our objective was to test whether p53 expression status is associated with survival for women diagnosed with the most common ovarian carcinoma histotypes (high-grade serous carcinoma [HGSC], endometrioid carcinoma [EC], and clear cell carcinoma [CCC]) using a large multi-institutional cohort from the Ovarian Tumor Tissue Analysis (OTTA) consortium. p53 expression was assessed on 6,678 cases represented on tissue microarrays from 25 participating OTTA study sites using a previously validated immunohistochemical (IHC) assay as a surrogate for the presence and functional effect of TP53 mutations. Three abnormal expression patterns (overexpression, complete absence, and cytoplasmic) and the normal (wild type) pattern were recorded. Survival analyses were performed by histotype. The frequency of abnormal p53 expression was 93.4% (4,630/4,957) in HGSC compared to 11.9% (116/973) in EC and 11.5% (86/748) in CCC. In HGSC, there were no differences in overall survival across the abnormal p53 expression patterns. However, in EC and CCC, abnormal p53 expression was associated with an increased risk of death for women diagnosed with EC in multivariate analysis compared to normal p53 as the reference (hazard ratio [HR] = 2.18, 95% confidence interval [CI] 1.36-3.47, p = 0.0011) and with CCC (HR = 1.57, 95% CI 1.11-2.22, p = 0.012). Abnormal p53 was also associated with shorter overall survival in The International Federation of Gynecology and Obstetrics stage I/II EC and CCC. Our study provides further evidence that functional groups of TP53 mutations assessed by abnormal surrogate p53 IHC patterns are not associated with survival in HGSC. In contrast, we validate that abnormal p53 IHC is a strong independent prognostic marker for EC and demonstrate for the first time an independent prognostic association of abnormal p53 IHC with overall survival in patients with CCC.

CCNE1 and survival of patients with tubo-ovarian high-grade serous carcinoma: An Ovarian Tumor Tissue Analysis consortium study

BACKGROUND

Cyclin E1 (CCNE1) is a potential predictive marker and therapeutic target in tubo-ovarian high-grade serous carcinoma (HGSC). Smaller studies have revealed unfavorable associations for CCNE1 amplification and CCNE1 overexpression with survival, but to date no large-scale, histotype-specific validation has been performed. The hypothesis was that high-level amplification of CCNE1 and CCNE1 overexpression, as well as a combination of the two, are linked to shorter overall survival in HGSC.

METHODS

Within the Ovarian Tumor Tissue Analysis consortium, amplification status and protein level in 3029 HGSC cases and mRNA expression in 2419 samples were investigated.

RESULTS

High-level amplification (>8 copies by chromogenic in situ hybridization) was found in 8.6% of HGSC and overexpression (>60% with at least 5% demonstrating strong intensity by immunohistochemistry) was found in 22.4%. CCNE1 high-level amplification and overexpression both were linked to shorter overall survival in multivariate survival analysis adjusted for age and stage, with hazard stratification by study (hazard ratio [HR], 1.26; 95% CI, 1.08-1.47, p = .034, and HR, 1.18; 95% CI, 1.05-1.32, p = .015, respectively). This was also true for cases with combined high-level amplification/overexpression (HR, 1.26; 95% CI, 1.09-1.47, p = .033). CCNE1 mRNA expression was not associated with overall survival (HR, 1.00 per 1-SD increase; 95% CI, 0.94-1.06; p = .58). CCNE1 high-level amplification is mutually exclusive with the presence of germline BRCA1/2 pathogenic variants and shows an inverse association to RB1 loss.

CONCLUSION

This study provides large-scale validation that CCNE1 high-level amplification is associated with shorter survival, supporting its utility as a prognostic biomarker in HGSC.

Increased FOXJ1 protein expression is associated with improved overall survival in high-grade serous ovarian carcinoma: an Ovarian Tumor Tissue Analysis Consortium Study

BACKGROUND

Recently, we showed a >60% difference in 5-year survival for patients with tubo-ovarian high-grade serous carcinoma (HGSC) when stratified by a 101-gene mRNA expression prognostic signature. Given the varied patient outcomes, this study aimed to translate prognostic mRNA markers into protein expression assays by immunohistochemistry and validate their survival association in HGSC.

METHODS

Two prognostic genes, FOXJ1 and GMNN, were selected based on high-quality antibodies, correlation with protein expression and variation in immunohistochemical scores in a preliminary cohort (n = 134 and n = 80, respectively). Six thousand four hundred and thirty-four (FOXJ1) and 5470 (GMNN) formalin-fixed, paraffin-embedded ovarian neoplasms (4634 and 4185 HGSC, respectively) represented on tissue microarrays from the Ovarian Tumor Tissue Analysis consortium underwent immunohistochemical staining and scoring, then univariate and multivariate survival analysis.

RESULTS

Consistent with mRNA, FOXJ1 protein expression exhibited a linear, increasing association with improved overall survival in HGSC patients. Women with >50% expression had the most favourable outcomes (HR = 0.78, 95% CI 0.67-0.91, p < 0.0001). GMNN protein expression was not significantly associated with overall HSGC patient survival. However, HGSCs with >35% GMNN expression showed a trend for better outcomes, though this was not significant.

CONCLUSION

We provide foundational evidence for the prognostic value of FOXJ1 in HGSC, validating the prior mRNA-based prognostic association by immunohistochemistry.

Profiling the immune landscape in mucinous ovarian carcinoma

OBJECTIVE

Mucinous ovarian carcinoma (MOC) is a rare histotype of ovarian cancer, with low response rates to standard chemotherapy, and very poor survival for patients diagnosed at advanced stage. There is a limited understanding of the MOC immune landscape, and consequently whether immune checkpoint inhibitors could be considered for a subset of patients.

METHODS

We performed multicolor immunohistochemistry (IHC) and immunofluorescence (IF) on tissue microarrays in a cohort of 126 MOC patients. Cell densities were calculated in the epithelial and stromal components for tumor-associated macrophages (CD68+/PD-L1+, CD68+/PD-L1-), T cells (CD3+/CD8-, CD3+/CD8+), putative T-regulatory cells (Tregs, FOXP3+), B cells (CD20+/CD79A+), plasma cells (CD20-/CD79a+), and PD-L1+ and PD-1+ cells, and compared these values with clinical factors. Univariate and multivariable Cox Proportional Hazards assessed overall survival. Unsupervised k-means clustering identified patient subsets with common patterns of immune cell infiltration.

RESULTS

Mean densities of PD1+ cells, PD-L1- macrophages, CD4+ and CD8+ T cells, and FOXP3+ Tregs were higher in the stroma compared to the epithelium. Tumors from advanced (Stage III/IV) MOC had greater epithelial infiltration of PD-L1- macrophages, and fewer PD-L1+ macrophages compared with Stage I/II cancers (p = 0.004 and p = 0.014 respectively). Patients with high epithelial density of FOXP3+ cells, CD8+/FOXP3+ cells, or PD-L1- macrophages, had poorer survival, and high epithelial CD79a + plasma cells conferred better survival, all upon univariate analysis only. Clustering showed that most MOC (86%) had an immune depleted (cold) phenotype, with only a small proportion (11/76,14%) considered immune inflamed (hot) based on T cell and PD-L1 infiltrates.

CONCLUSION

In summary, MOCs are mostly immunogenically 'cold', suggesting they may have limited response to current immunotherapies.

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.

MCM3 is a novel proliferation marker associated with longer survival for patients with tubo-ovarian high-grade serous carcinoma

Tubo-ovarian high-grade serous carcinomas (HGSC) are highly proliferative neoplasms that generally respond well to platinum/taxane chemotherapy. We recently identified minichromosome maintenance complex component 3 (MCM3), which is involved in the initiation of DNA replication and proliferation, as a favorable prognostic marker in HGSC. Our objective was to further validate whether MCM3 mRNA expression and possibly MCM3 protein levels are associated with survival in patients with HGSC. MCM3 mRNA expression was measured using NanoString expression profiling on formalin-fixed and paraffin-embedded tissue (N = 2355 HGSC) and MCM3 protein expression was assessed by immunohistochemistry (N = 522 HGSC) and compared with Ki-67. Kaplan-Meier curves and the Cox proportional hazards model were used to estimate associations with survival. Among chemotherapy-naïve HGSC, higher MCM3 mRNA expression (one standard deviation increase in the score) was associated with longer overall survival (HR = 0.87, 95% CI 0.81-0.92, p < 0.0001, N = 1840) in multivariable analysis. MCM3 mRNA expression was highest in the HGSC C5.PRO molecular subtype, although no interaction was observed between MCM3, survival and molecular subtypes. MCM3 and Ki-67 protein levels were significantly lower after exposure to neoadjuvant chemotherapy compared to chemotherapy-naïve tumors: 37.0% versus 46.4% and 22.9% versus 34.2%, respectively. Among chemotherapy-naïve HGSC, high MCM3 protein levels were also associated with significantly longer disease-specific survival (HR = 0.52, 95% CI 0.36-0.74, p = 0.0003, N = 392) compared to cases with low MCM3 protein levels in multivariable analysis. MCM3 immunohistochemistry is a promising surrogate marker of proliferation in HGSC.

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.

Adjuvant chemotherapy and radiation for patients with high-risk stage I endometrial cancer treated with curative intent surgery: impact on recurrence and survival

BACKGROUND

Survival benefits of post-operative systemic and radiation therapy in high-risk stage I endometrial cancer are uncertain.

OBJECTIVE

To compare recurrence patterns and survival outcomes of post-surgical treatment in patients with high-risk stage I endometrial cancer and to determine whether adjuvant therapy significantly improves outcomes.

METHODS

High-risk stage I endometrial cancer was defined as either stage IB grade 3 endometrioid histology or myoinvasive non-endometrioid histology. Consecutive patients diagnosed between January 2000 and December 2010 in eight cancer centers were included. Patients, disease, and treatment characteristics were summarized by descriptive statistics. Overall survival, disease-specific survival, and relapse-free survival were examined using Cox's proportional hazards regression and log-rank test. Survival curves were estimated using the Kaplan-Meier method.

RESULTS

Of 2317 patients with stage I endometrial cancer, 414 patients had high-risk disease. Use of chemotherapy did not improve overall survival (relative risk (RR) 0.70, 95% CI 0.46 to 1.14, p=0.13) or disease-specific survival (RR 1.06, 95% CI 0.61 to 1.85, p=0.84). Significant improvement in recurrence-free survival was observed in patients who received chemotherapy (RR 0.61, 95% CI 0.39 to 0.95, p=0.03). Use of radiation therapy did not improve overall survival, recurrence-free survival, or disease-specific survival. Patients who received four cycles or fewer of chemotherapy versus five to six cycles had similar overall survival, disease-specific survival, and recurrence-free survival.

CONCLUSIONS

Post-operative chemotherapy or radiation in stage I high-risk endometrial cancer is not associated with improved cancer-specific or overall survival. More than four cycles of chemotherapy did not improve survival compared with four cycles or fewer.

Accurate Distinction of Ovarian Clear Cell From Endometrioid Carcinoma Requires Integration of Phenotype, Immunohistochemical Predictions, and Genotype: Implications for Lynch Syndrome Screening

Ovarian clear cell carcinoma (OCCC) and ovarian endometrioid carcinoma (OEC) are both associated with endometriosis but differ in histologic phenotype, biomarker profile, and survival. Our objectives were to refine immunohistochemical (IHC) panels that help distinguish the histotypes and reassess the prevalence of mismatch repair deficiency (MMRd) in immunohistochemically confirmed OCCC. We selected 8 candidate IHC markers to develop first-line and second-line panels in a training set of 344 OCCC/OEC cases. Interobserver reproducibility of histotype diagnosis was assessed in an independent testing cohort of 100 OCC/OEC initially without and subsequently with IHC. The prevalence of MMRd was evaluated using the testing cohort and an expansion set of 844 ovarian carcinomas. The 2 prototypical combinations (OCCC: Napsin A+/HNF1B diffusely+/PR-; OEC: Napsin A-/HNF1B nondiffuse/PR+) occurred in 75% of cases and were 100% specific. A second-line panel (ELAPOR1, AMACR, CDX2) predicted the remaining cases with 83% accuracy. Integration of IHC improved interobserver reproducibility (κ=0.778 vs. 0.882, P<0.0001). The prevalence of MMRd was highest in OEC (11.5%, 44/383), lower in OCCC (1.7%, 5/297), and high-grade serous carcinomas (0.7%, 5/699), and absent in mucinous (0/126) and low-grade serous carcinomas (0/50). All 5 MMRd OCCC were probable Lynch syndrome cases with prototypical IHC profile but ambiguous morphologic features: 3/5 with microcystic architecture and 2/5 with intratumoral stromal inflammation. Integration of first-line and second-line IHC panels increases diagnostic precision and enhances prognostication and triaging for predisposing/predictive molecular biomarker testing. Our data support universal Lynch syndrome screening in all patients with OEC when the diagnosis of other histotypes has been vigorously excluded.

Refined cut-off for TP53 immunohistochemistry improves prediction of TP53 mutation status in ovarian mucinous tumors: implications for outcome analyses

TP53 mutations are implicated in the progression of mucinous borderline tumors (MBOT) to mucinous ovarian carcinomas (MOC). Optimized immunohistochemistry (IHC) for TP53 has been established as a proxy for the TP53 mutation status in other ovarian tumor types. We aimed to confirm the ability of TP53 IHC to predict TP53 mutation status in ovarian mucinous tumors and to evaluate the association of TP53 mutation status with survival among patients with MBOT and MOC. Tumor tissue from an initial cohort of 113 women with MBOT/MOC was stained with optimized IHC for TP53 using tissue microarrays (75.2%) or full sections (24.8%) and interpreted using established criteria as normal or abnormal (overexpression, complete absence, or cytoplasmic). Cases were considered concordant if abnormal IHC staining predicted deleterious TP53 mutations. Discordant tissue microarray cases were re-evaluated on full sections and interpretational criteria were refined. The initial cohort was expanded to a total of 165 MBOT and 424 MOC for the examination of the association of survival with TP53 mutation status, assessed either by TP53 IHC and/or sequencing. Initially, 82/113 (72.6%) cases were concordant using the established criteria. Refined criteria for overexpression to account for intratumoral heterogeneity and terminal differentiation improved concordance to 93.8% (106/113). In the expanded cohort, 19.4% (32/165) of MBOT showed evidence for TP53 mutation and this was associated with a higher risk of recurrence, disease-specific death, and all-cause mortality (overall survival: HR = 4.6, 95% CI 1.5-14.3, p = 0.0087). Within MOC, 61.1% (259/424) harbored a TP53 mutation, but this was not associated with survival (overall survival, p = 0.77). TP53 IHC is an accurate proxy for TP53 mutation status with refined interpretation criteria accounting for intratumoral heterogeneity and terminal differentiation in ovarian mucinous tumors. TP53 mutation status is an important biomarker to identify MBOT with a higher risk of mortality.

Tattoos and Hematologic Malignancies in British Columbia, Canada

BACKGROUND

Tattoos may cause a variety of adverse reactions in the body, including immune reactions and infections. However, it is unknown whether tattoos may increase the risk of lymphatic cancers such as non-Hodgkin lymphoma (NHL) and multiple myeloma.

METHODS

Participants from two population-based case-control studies were included in logistic regression models to examine the association between tattoos and risk of NHL and multiple myeloma.

RESULTS

A total of 1,518 participants from the NHL study (737 cases) and 742 participants from the multiple myeloma study (373 cases) were included in the analyses. No statistically significant associations were found between tattoos and risk of NHL or multiple myeloma after adjusting for age, sex, ethnicity, education, body mass index, and family history.

CONCLUSIONS

We did not identify any significant associations between tattoos and risk of multiple myeloma, NHL, or NHL subtypes in these studies.

IMPACT

Though biologically plausible, tattoos were not associated with increased risk of NHL or multiple myeloma in this study. Future studies with greater detail regarding tattoo exposure may provide further insights.

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.

Pathological chemotherapy response score is prognostic in tubo-ovarian high-grade serous carcinoma: A systematic review and meta-analysis of individual patient data

OBJECTIVE

There is a need to develop and validate biomarkers for treatment response and survival in tubo-ovarian high-grade serous carcinoma (HGSC). The chemotherapy response score (CRS) stratifies patients into complete/near-complete (CRS3), partial (CRS2), and no/minimal (CRS1) response after neoadjuvant chemotherapy (NACT). Our aim was to review current evidence to determine whether the CRS is prognostic in women with tubo-ovarian HGSC treated with NACT.

METHODS

We established an international collaboration to conduct a systematic review and meta-analysis, pooling individual patient data from 16 sites in 11 countries. Patients had stage IIIC/IV HGSC, 3-4 NACT cycles and >6-months follow-up. Random effects models were used to derive combined odds ratios in the pooled population to investigate associations between CRS and progression free and overall survival (PFS and OS).

RESULTS

877 patients were included from published and unpublished studies. Median PFS and OS were 15 months (IQR 5-65) and 28 months (IQR 7-92) respectively. CRS3 was seen in 249 patients (28%). The pooled hazard ratios (HR) for PFS and OS for CRS3 versus CRS1/CRS2 were 0·55 (95% CI, 0·45-0·66; P < 0·001) and 0·65 (95% CI 0·50-0·85, P = 0·002) respectively; no heterogeneity was identified (PFS: Q = 6·42, P = 0·698, I2 = 0·0%; OS: Q = 6·89, P = 0·648, I2 = 0·0%). CRS was significantly associated with PFS and OS in multivariate models adjusting for age and stage. Of 306 patients with known germline BRCA1/2 status, those with BRCA1/2 mutations (n = 80) were more likely to achieve CRS3 (P = 0·027).

CONCLUSIONS

CRS3 was significantly associated with improved PFS and OS compared to CRS1/2. This validation of CRS in a real-world setting demonstrates it to be a robust and reproducible biomarker with potential to be incorporated into therapeutic decision-making and clinical trial design.

Genome-wide association studies identify susceptibility loci for epithelial ovarian cancer in east Asian women

OBJECTIVE

Genome-wide association studies (GWASs) for epithelial ovarian cancer (EOC) have focused largely on populations of European ancestry. We aimed to identify common germline variants associated with EOC risk in Asian women.

METHODS

Genotyping was performed as part of the OncoArray project. Samples with >60% Asian ancestry were included in the analysis. Genotyping was performed on 533,631 SNPs in 3238 Asian subjects diagnosed with invasive or borderline EOC and 4083 unaffected controls. After imputation, genotypes were available for 11,595,112 SNPs to identify associations.

RESULTS

At chromosome 6p25.2, SNP rs7748275 was associated with risk of serous EOC (odds ratio [OR] = 1.34, P = 8.7 × 10-9) and high-grade serous EOC (HGSOC) (OR = 1.34, P = 4.3 × 10-9). SNP rs6902488 at 6p25.2 (r2 = 0.97 with rs7748275) lies in an active enhancer and is predicted to impact binding of STAT3, P300 and ELF1. We identified additional risk loci with low Bayesian false discovery probability (BFDP) scores, indicating they are likely to be true risk associations (BFDP <10%). At chromosome 20q11.22, rs74272064 was associated with HGSOC risk (OR = 1.27, P = 9.0 × 10-8). Overall EOC risk was associated with rs10260419 at chromosome 7p21.3 (OR = 1.33, P = 1.2 × 10-7) and rs74917072 at chromosome 2q37.3 (OR = 1.25, P = 4.7 × 10-7). At 2q37.3, expression quantitative trait locus analysis in 404 HGSOC tissues identified ESPNL as a putative candidate susceptibility gene (P = 1.2 × 10-7).

CONCLUSION

While some risk loci were shared between East Asian and European populations, others were population-specific, indicating that the landscape of EOC risk in Asian women has both shared and unique features compared to women of European ancestry.

Evaluation of vitamin D biosynthesis and pathway target genes reveals UGT2A1/2 and EGFR polymorphisms associated with epithelial ovarian cancer in African American Women

An association between genetic variants in the vitamin D receptor (VDR) gene and epithelial ovarian cancer (EOC) was previously reported in women of African ancestry (AA). We sought to examine associations between genetic variants in VDR and additional genes from vitamin D biosynthesis and pathway targets (EGFR, UGT1A, UGT2A1/2, UGT2B, CYP3A4/5, CYP2R1, CYP27B1, CYP24A1, CYP11A1, and GC). Genotyping was performed using the custom-designed 533,631 SNP Illumina OncoArray with imputation to the 1,000 Genomes Phase 3 v5 reference set in 755 EOC cases, including 537 high-grade serous (HGSOC), and 1,235 controls. All subjects are of African ancestry (AA). Logistic regression was performed to estimate odds ratios (OR) and 95% confidence intervals (CI). We further evaluated statistical significance of selected SNPs using the Bayesian False Discovery Probability (BFDP). A significant association with EOC was identified in the UGT2A1/2 region for the SNP rs10017134 (per allele OR = 1.4, 95% CI = 1.2-1.7, P = 1.2 × 10-6 , BFDP = 0.02); and an association with HGSOC was identified in the EGFR region for the SNP rs114972508 (per allele OR = 2.3, 95% CI = 1.6-3.4, P = 1.6 × 10-5 , BFDP = 0.29) and in the UGT2A1/2 region again for rs1017134 (per allele OR = 1.4, 95% CI = 1.2-1.7, P = 2.3 × 10-5 , BFDP = 0.23). Genetic variants in the EGFR and UGT2A1/2 may increase susceptibility of EOC in AA women. Future studies to validate these findings are warranted. Alterations in EGFR and UGT2A1/2 could perturb enzyme efficacy, proliferation in ovaries, impact and mark susceptibility to EOC.

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.

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.

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

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

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

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

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

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

The Super-Seniors Study: Phenotypic characterization of a healthy 85+ population

Background

To understand why some people live to advanced age in good health and others do not, it is important to study not only disease, but also long-term good health. The Super-Seniors Study aims to identify factors associated with healthy aging.

Methods

480 healthy oldest-old ‘Super-Seniors’ aged 85 to 105 years and never diagnosed with cancer, cardiovascular disease, diabetes, dementia, or major pulmonary disease, were compared to 545 mid-life controls aged 41–54, who represent a group that is unselected for survival from late-life diseases. Health and lifestyle information, personal and family medical history, and blood samples were collected from all participants. Super-Seniors also underwent four geriatric tests.

Results

Super-Seniors showed high cognitive (Mini-Mental State Exam mean = 28.3) and functional capacity (Instrumental Activities of Daily Living Scale mean = 21.4), as well as high physical function (Timed Up and Go mean = 12.3 seconds) and low levels of depression (Geriatric Depression Scale mean = 1.5). Super-Seniors were less likely to be current smokers than controls, but the frequency of drinking alcohol was the same in both groups. Super-Seniors were more likely to have 4 or more offspring; controls were more likely to have no children. Female Super-Seniors had a mean age of last fertility 1.9 years older than controls, and were 2.3 times more likely to have had a child at ≥ 40 years. The parents of Super-Seniors had mean ages of deaths of 79.3 years for mothers, and 74.5 years for fathers, each exceeding the life expectancy for their era by a decade.

Conclusions

Super-Seniors are cognitively and physically high functioning individuals who have evaded major age-related chronic diseases into old age, representing the approximately top 1% for healthspan. The familiality of long lifespan of the parents of Super-Seniors supports the hypothesis that heritable factors contribute to this desirable phenotype.

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.

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.

A targeted genetic association study of epithelial ovarian cancer susceptibility

BACKGROUND

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

RESULTS

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

METHODS

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

CONCLUSION

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Tea, coffee, and caffeinated beverage consumption and risk of epithelial ovarian cancers

BACKGROUND

The risk for epithelial ovarian cancer associated with the consumption of caffeinated beverages (tea, coffee, and soft drinks) and green tea is inconclusive. However, few studies have investigated the type of caffeinated beverage or the type of tea.

OBJECTIVE

We assessed consumption of tea (black/caffeinated tea and green tea separately), coffee, and caffeinated soft drinks, as well as level of consumption, and the risk for epithelial ovarian cancer and its histotypes.

STUDY DESIGN

This study was conducted within a population-based case-control study in Alberta and British Columbia, Canada from 2001 to 2012. After restricting to cases of epithelial invasive cancers and controls aged 40-79 years who completed an interview that included coffee, soft drink, and tea consumption (ascertained starting in 2005 in British Columbia and 2008 in Alberta), there were a total of 524 cases and 1587 controls. Those that did not meet the threshold for beverage consumption (at least once per month for 6 months or more) were classified as non-drinkers. Adult lifetime cumulative consumption (cup-years=cups/day*years) was calculated. Unconditional logistic regression was used to estimate adjusted odds ratios (aOR) and 95% confidence intervals (CI) to describe the association between the relevant drink consumption and risk.

RESULTS

No excess risk was seen for coffee or caffeinated soft drinks. Similarly, any tea consumption was not associated with risk, but when stratified by the type of tea, there was an increase in risk in black tea only drinkers (aOR=1.56; 95% CI:1.07-2.28 for >40 cup-years), but no excess risk for the exclusive green tea drinkers. Similar findings were observed for post-menopausal women. The association for black tea only consumption was mainly seen in the endometrioid histotype (aOR=3.19; 95% CI: 1.32-7.69).

CONCLUSION

Black tea consumption may be associated with an increased risk epithelial ovarian carcinoma. The excess risk is seen only in the endometrioid histotype but not in serous or clear cell. Further studies are required to confirm these findings and identify the constituents in black tea that may increase the risk.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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.

Genome-Wide Meta-Analyses of Breast, Ovarian, and Prostate Cancer Association Studies Identify Multiple New Susceptibility Loci Shared by at Least Two Cancer Types

Breast, ovarian, and prostate cancers are hormone-related and may have a shared genetic basis, but this has not been investigated systematically by genome-wide association (GWA) studies. Meta-analyses combining the largest GWA meta-analysis data sets for these cancers totaling 112,349 cases and 116,421 controls of European ancestry, all together and in pairs, identified at P < 10(-8) seven new cross-cancer loci: three associated with susceptibility to all three cancers (rs17041869/2q13/BCL2L11; rs7937840/11q12/INCENP; rs1469713/19p13/GATAD2A), two breast and ovarian cancer risk loci (rs200182588/9q31/SMC2; rs8037137/15q26/RCCD1), and two breast and prostate cancer risk loci (rs5013329/1p34/NSUN4; rs9375701/6q23/L3MBTL3). Index variants in five additional regions previously associated with only one cancer also showed clear association with a second cancer type. Cell-type-specific expression quantitative trait locus and enhancer-gene interaction annotations suggested target genes with potential cross-cancer roles at the new loci. Pathway analysis revealed significant enrichment of death receptor signaling genes near loci with P < 10(-5) in the three-cancer meta-analysis.

SIGNIFICANCE

We demonstrate that combining large-scale GWA meta-analysis findings across cancer types can identify completely new risk loci common to breast, ovarian, and prostate cancers. We show that the identification of such cross-cancer risk loci has the potential to shed new light on the shared biology underlying these hormone-related cancers. Cancer Discov; 6(9); 1052-67. ©2016 AACR.This article is highlighted in the In This Issue feature, p. 932.

Exome genotyping arrays to identify rare and low frequency variants associated with epithelial ovarian cancer risk

Rare and low frequency variants are not well covered in most germline genotyping arrays and are understudied in relation to epithelial ovarian cancer (EOC) risk. To address this gap, we used genotyping arrays targeting rarer protein-coding variation in 8,165 EOC cases and 11,619 controls from the international Ovarian Cancer Association Consortium (OCAC). Pooled association analyses were conducted at the variant and gene level for 98,543 variants directly genotyped through two exome genotyping projects. Only common variants that represent or are in strong linkage disequilibrium (LD) with previously-identified signals at established loci reached traditional thresholds for exome-wide significance (P < 5.0 × 10 -  7). One of the most significant signals (Pall histologies = 1.01 × 10 -  13;Pserous = 3.54 × 10 -  14) occurred at 3q25.31 for rs62273959, a missense variant mapping to the LEKR1 gene that is in LD (r2 = 0.90) with a previously identified 'best hit' (rs7651446) mapping to an intron of TIPARP. Suggestive associations (5.0 × 10 -  5 > P≥5.0 ×10 -  7) were detected for rare and low-frequency variants at 16 novel loci. Four rare missense variants were identified (ACTBL2 rs73757391 (5q11.2), BTD rs200337373 (3p25.1), KRT13 rs150321809 (17q21.2) and MC2R rs104894658 (18p11.21)), but only MC2R rs104894668 had a large effect size (OR = 9.66). Genes most strongly associated with EOC risk included ACTBL2 (PAML = 3.23 × 10 -  5; PSKAT-o = 9.23 × 10 -  4) and KRT13 (PAML = 1.67 × 10 -  4; PSKAT-o = 1.07 × 10 -  5), reaffirming variant-level analysis. In summary, this large study identified several rare and low-frequency variants and genes that may contribute to EOC susceptibility, albeit with possible small effects. Future studies that integrate epidemiology, sequencing, and functional assays are needed to further unravel the unexplained heritability and biology of this disease.

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.

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.

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.

Hormone Contraception before the First Birth and Ovarian Cancer Risk

Introduction: Combined oral contraceptive (OC) use strongly and consistently reduces the risk for epithelial ovarian cancer (EOC); longer durations of use and more recent use are associated with the strongest reductions in risk. However, it is unknown if exclusive OC use before the first birth is associated with a reduction in EOC risk many years later. Therefore, we investigated the risk for EOC among parous women associated with exclusive OC use before the first birth. Methods: From a population-based case-control study in Alberta and British Columbia, Canada, 2001–2011, we included 1144 invasive EOC cases and 2513 controls who were &gt;40 years at diagnosis/reference date. Participants reported OC use and all pregnancies via a telephone interview or self-administered questionnaire (in the early years of the study). Duration of OC use was evaluated as a continuous variable and by categories: non-users (never use or &lt;6 months of use), &lt;5, 5–&lt;10, &gt;10 years, unknown. Using logistic regression, we estimated adjusted odds ratios (aORs) and 95% confidence intervals (CI), controlling for study site, age, parity, breastfeeding, first degree family history of breast/ovarian cancer, tubal ligation, and BMI. Results: OC use at any time during reproductive life was associated with a 42% reduced risk for EOC relative to non-users (aOR = 0.58, 95% CI, 0.49–0.69). Among parous women, each additional year of exclusive OC use before the first birth conferred an 11% risk reduction relative to non-users (aOR = 0.89 95% CI, 0.86–0.94, linear trend p-value &lt;0.01). Results were similar when we restricted to cases with high grade serous cancers (aOR = 0.90 95% CI, 0.84–0.95, linear trend p-value &lt;0.01) and for cases with endometrioid/clear cell cancer (aOR 0.88 95% CI, 0.80–0.95, linear trend p-value &lt; 0.01). Discussion: Among parous women, exclusive use of OCs before the first birth was associated with a strong reduction in EOC risk many years later. Because OCs stop ovulation, this reduced risk may be due to a reduction in lifetime ovulatory cycles. However, it is also possible that OC use at younger ages, before the first birth, represents a window of opportunity to have a substantial impact on reducing risk that remains for many years, informing possible prevention strategies.

Lifetime physical activity and the risk of non-Hodgkin lymphoma

Research regarding the association between physical activity and the risk of non-Hodgkin lymphoma (NHL) is limited and inconsistent, and few studies have investigated whether the intensity and timing of physical activity influence the association. A case-control study of NHL was conducted in British Columbia, Canada, in 2000 to 2004. Data were collected on various NHL risk factors, including moderate-intensity and vigorous-intensity physical activity performed over the lifetime. Logistic regression was used to estimate the association between physical activity and the risk of NHL. This analysis included 818 controls and 749 cases. Lifetime vigorous-intensity physical activity was inversely associated with NHL risk. Participants in the second, third, and fourth quartiles of lifetime vigorous-intensity physical activity had an approximately 25% to 30% lower risk of NHL than those in the lowest quartile [adjusted odds ratios, 0.69 (95% confidence interval [CI], 0.52-0.93); 0.68 (95% CI, 0.50-0.92); and 0.75 (95% CI, 0.55-1.01), respectively]. No consistent associations were observed for total or moderate-intensity physical activity. There were no apparent age periods in which physical activity appeared to confer a greater risk reduction. In this study, we found that lifetime vigorous-intensity physical activity was associated with a significantly reduced risk of NHL. Given this finding, more research on physical activity intensity and timing in relation to NHL risk is warranted.

Adult lifetime alcohol consumption and invasive epithelial ovarian cancer risk in a population-based case-control study

OBJECTIVE

Meta-analyses report a null association between recent alcohol consumption and ovarian cancer risk. However, because few studies investigated different types of alcohol over adult ages, we investigated adult lifetime and type (beer, wine, spirits) of consumption and risk.

METHODS

Consumption after age 20years was ascertained in 1144 invasive epithelial ovarian cancer cases and 2513 controls in a population-based case-control study (Alberta and British Columbia, Canada, 2001-2012). Non-drinkers consumed any types of alcohol <12 times per year on average. Logistic regression was use to estimate adjusted odds ratios [aOR] and 95% confidence intervals [CIs].

RESULTS

Wine consumption was associated with a risk reduction (aOR=0.67, 95% CI: 0.50-0.88) relative to non-drinkers, but not beer (aOR=1.06, 95% CI: 0.71-1.58) or spirits (aOR=0.98, 95% CI: 0.69-1.39). The reduced risk was stronger for exclusive red wine drinkers (aOR=0.44, 95% CI: 0.19-0.92) than white wine drinkers (aOR=0.79, 95% CI: 0.46-1.34), although most women drank both types of wine. Risk decreased with increasing cumulative consumption of any wine (P-trend<0.05) and was evident for the serous histotype. Wine consumption initiated prior to age 50 was associated with a risk reduction (e.g., at 40-49years, aOR=0.58, 95% CI: 0.42-0.78), but not drinking initiated after 50years of age. For any type, level, or age at initiation of alcohol consumption, we found no increased risks.

CONCLUSIONS

For the moderate consumption in this study, higher levels of wine consumption were generally associated with risk reductions; reductions may be stronger for red wine. Our results suggest that alcohol consumption that is guideline concordant will not increase epithelial ovarian cancer risk.

Epithelial-Mesenchymal Transition (EMT) Gene Variants and Epithelial Ovarian Cancer (EOC) Risk

Epithelial-mesenchymal transition (EMT) is a process whereby epithelial cells assume mesenchymal characteristics to facilitate cancer metastasis. However, EMT also contributes to the initiation and development of primary tumors. Prior studies that explored the hypothesis that EMT gene variants contribute to epithelial ovarian carcinoma (EOC) risk have been based on small sample sizes and none have sought replication in an independent population. We screened 15,816 single-nucleotide polymorphisms (SNPs) in 296 genes in a discovery phase using data from a genome-wide association study of EOC among women of European ancestry (1,947 cases and 2,009 controls) and identified 793 variants in 278 EMT-related genes that were nominally (P < 0.05) associated with invasive EOC. These SNPs were then genotyped in a larger study of 14,525 invasive-cancer patients and 23,447 controls. A P-value <0.05 and a false discovery rate (FDR) <0.2 were considered statistically significant. In the larger dataset, GPC6/GPC5 rs17702471 was associated with the endometrioid subtype among Caucasians (odds ratio (OR) = 1.16, 95% CI = 1.07-1.25, P = 0.0003, FDR = 0.19), whereas F8 rs7053448 (OR = 1.69, 95% CI = 1.27-2.24, P = 0.0003, FDR = 0.12), F8 rs7058826 (OR = 1.69, 95% CI = 1.27-2.24, P = 0.0003, FDR = 0.12), and CAPN13 rs1983383 (OR = 0.79, 95% CI = 0.69-0.90, P = 0.0005, FDR = 0.12) were associated with combined invasive EOC among Asians. In silico functional analyses revealed that GPC6/GPC5 rs17702471 coincided with DNA regulatory elements. These results suggest that EMT gene variants do not appear to play a significant role in the susceptibility to EOC.

Common variants at the CHEK2 gene 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.

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.

Network-Based Integration of GWAS and Gene Expression Identifies a HOX-Centric Network Associated with Serous Ovarian Cancer Risk

BACKGROUND

Genome-wide association studies (GWAS) have so far reported 12 loci associated with serous epithelial ovarian cancer (EOC) risk. We hypothesized that some of these loci function through nearby transcription factor (TF) genes and that putative target genes of these TFs as identified by coexpression may also be enriched for additional EOC risk associations.

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.

IMPACT

Network analysis integrating large, context-specific datasets has the potential to offer mechanistic insights into cancer susceptibility and prioritize genes for experimental characterization.

Cis-eQTL analysis and functional validation of candidate susceptibility genes for high-grade serous ovarian cancer

Genome-wide association studies have reported 11 regions conferring risk of high-grade serous epithelial ovarian cancer (HGSOC). Expression quantitative trait locus (eQTL) analyses can identify candidate susceptibility genes at risk loci. Here we evaluate cis-eQTL associations at 47 regions associated with HGSOC risk (P≤10(-5)). For three cis-eQTL associations (P<1.4 × 10(-3), FDR<0.05) at 1p36 (CDC42), 1p34 (CDCA8) and 2q31 (HOXD9), we evaluate the functional role of each candidate by perturbing expression of each gene in HGSOC precursor cells. Overexpression of HOXD9 increases anchorage-independent growth, shortens population-doubling time and reduces contact inhibition. Chromosome conformation capture identifies an interaction between rs2857532 and the HOXD9 promoter, suggesting this SNP is a leading causal variant. Transcriptomic profiling after HOXD9 overexpression reveals enrichment of HGSOC risk variants within HOXD9 target genes (P=6 × 10(-10) for risk variants (P<10(-4)) within 10 kb of a HOXD9 target gene in ovarian cells), suggesting a broader role for this network in genetic susceptibility to HGSOC.

Common Genetic Variation in Circadian Rhythm Genes and Risk of Epithelial Ovarian Cancer (EOC)

Disruption in circadian gene expression, whether due to genetic variation or environmental factors (e.g., light at night, shiftwork), is associated with increased incidence of breast, prostate, gastrointestinal and hematologic cancers and gliomas. Circadian genes are highly expressed in the ovaries where they regulate ovulation; circadian disruption is associated with several ovarian cancer risk factors (e.g., endometriosis). However, no studies have examined variation in germline circadian genes as predictors of ovarian cancer risk and invasiveness. The goal of the current study was to examine single nucleotide polymorphisms (SNPs) in circadian genes BMAL1, CRY2, CSNK1E, NPAS2, PER3, REV1 and TIMELESS and downstream transcription factors KLF10 and SENP3 as predictors of risk of epithelial ovarian cancer (EOC) and histopathologic subtypes. The study included a test set of 3,761 EOC cases and 2,722 controls and a validation set of 44,308 samples including 18,174 (10,316 serous) cases and 26,134 controls from 43 studies participating in the Ovarian Cancer Association Consortium (OCAC). Analysis of genotype data from 36 genotyped SNPs and 4600 imputed SNPs indicated that the most significant association was rs117104877 in BMAL1 (OR = 0.79, 95% CI = 0.68-0.90, p = 5.59 × 10^-4]. Functional analysis revealed a significant down regulation of BMAL1 expression following cMYC overexpression and increasing transformation in ovarian surface epithelial (OSE) cells as well as alternative splicing of BMAL1 exons in ovarian and granulosa cells. These results suggest that variation in circadian genes, and specifically BMAL1, may be associated with risk of ovarian cancer, likely through disruption of hormonal pathways.