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.

Discovery and validation of serum glycoprotein biomarkers for high grade serous ovarian cancer

PURPOSE

This study aimed to identify serum glycoprotein biomarkers for early detection of high-grade serous ovarian cancer (HGSOC), the most common and aggressive histotype of ovarian cancer.

EXPERIMENTAL DESIGN

The glycoproteomics pipeline lectin magnetic bead array (LeMBA)-mass spectrometry (MS) was used in age-matched case-control serum samples. Clinical samples collected at diagnosis were divided into discovery (n = 30) and validation (n = 98) sets. We also analysed a set of preclinical sera (n = 30) collected prior to HGSOC diagnosis in the UK Collaborative Trial of Ovarian Cancer Screening.

RESULTS

A 7-lectin LeMBA-MS/MS discovery screen shortlisted 59 candidate proteins and three lectins. Validation analysis using 3-lectin LeMBA-multiple reaction monitoring (MRM) confirmed elevated A1AT, AACT, CO9, HPT and ITIH3 and reduced A2MG, ALS, IBP3 and PON1 glycoforms in HGSOC. The best performing multimarker signature had 87.7% area under the receiver operating curve, 90.7% specificity and 70.4% sensitivity for distinguishing HGSOC from benign and healthy groups. In the preclinical set, CO9, ITIH3 and A2MG glycoforms were altered in samples collected 11.1 ± 5.1 months prior to HGSOC diagnosis, suggesting potential for early detection.

CONCLUSIONS AND CLINICAL RELEVANCE

Our findings provide evidence of candidate early HGSOC serum glycoprotein biomarkers, laying the foundation for further study in larger cohorts.

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.

Methylation Signature Implicated in Immuno-Suppressive Activities in Tubo-Ovarian High-Grade Serous Carcinoma

BACKGROUND

Better understanding of prognostic factors in tubo-ovarian high-grade serous carcinoma (HGSC) is critical, as diagnosis confers an aggressive disease course. Variation in tumor DNA methylation shows promise predicting outcome, yet prior studies were largely platform-specific and unable to evaluate multiple molecular features.

METHODS

We analyzed genome-wide DNA methylation in 1,040 frozen HGSC, including 325 previously reported upon, seeking a multi-platform quantitative methylation signature that we evaluated in relation to clinical features, tumor characteristics, time to recurrence/death, extent of CD8+ tumor-infiltrating lymphocytes (TIL), gene expression molecular subtypes, and gene expression of the ATP-binding cassette transporter TAP1.

RESULTS

Methylation signature was associated with shorter time to recurrence, independent of clinical factors (N = 715 new set, hazard ratio (HR), 1.65; 95% confidence interval (CI), 1.10-2.46; P = 0.015; N = 325 published set HR, 2.87; 95% CI, 2.17-3.81; P = 2.2 × 10-13) and remained prognostic after adjustment for gene expression molecular subtype and TAP1 expression (N = 599; HR, 2.22; 95% CI, 1.66-2.95; P = 4.1 × 10-8). Methylation signature was inversely related to CD8+ TIL levels (P = 2.4 × 10-7) and TAP1 expression (P = 0.0011) and was associated with gene expression molecular subtype (P = 5.9 × 10-4) in covariate-adjusted analysis.

CONCLUSIONS

Multi-center analysis identified a novel quantitative tumor methylation signature of HGSC applicable to numerous commercially available platforms indicative of shorter time to recurrence/death, adjusting for other factors. Along with immune cell composition analysis, these results suggest a role for DNA methylation in the immunosuppressive microenvironment.

IMPACT

This work aids in identification of targetable epigenome processes and stratification of patients for whom tailored treatment may be most beneficial.

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.

Association of p16 expression with prognosis varies across ovarian carcinoma histotypes: an Ovarian Tumor Tissue Analysis consortium study

We aimed to validate the prognostic association of p16 expression in ovarian high-grade serous carcinomas (HGSC) and to explore it in other ovarian carcinoma histotypes. p16 protein expression was assessed by clinical-grade immunohistochemistry in 6525 ovarian carcinomas including 4334 HGSC using tissue microarrays from 24 studies participating in the Ovarian Tumor Tissue Analysis consortium. p16 expression patterns were interpreted as abnormal (either overexpression referred to as block expression or absence) or normal (heterogeneous). CDKN2A (which encodes p16) mRNA expression was also analyzed in a subset (n = 2280) mostly representing HGSC (n = 2010). Association of p16 expression with overall survival (OS) was determined within histotypes as was CDKN2A expression for HGSC only. p16 block expression was most frequent in HGSC (56%) but neither protein nor mRNA expression was associated with OS. However, relative to heterogeneous expression, block expression was associated with shorter OS in endometriosis-associated carcinomas, clear cell [hazard ratio (HR): 2.02, 95% confidence (CI) 1.47-2.77, p < 0.001] and endometrioid (HR: 1.88, 95% CI 1.30-2.75, p = 0.004), while absence was associated with shorter OS in low-grade serous carcinomas (HR: 2.95, 95% CI 1.61-5.38, p = 0.001). Absence was most frequent in mucinous carcinoma (50%), and was not associated with OS in this histotype. The prognostic value of p16 expression is histotype-specific and pattern dependent. We provide definitive evidence against an association of p16 expression with survival in ovarian HGSC as previously suggested. Block expression of p16 in clear cell and endometrioid carcinoma should be further validated as a prognostic marker, and absence in low-grade serous carcinoma justifies CDK4 inhibition.

Racial/ethnic differences in the epidemiology of ovarian cancer: a pooled analysis of 12 case-control studies

Background

Ovarian cancer incidence differs substantially by race/ethnicity, but the reasons for this are not well understood. Data were pooled from the African American Cancer Epidemiology Study (AACES) and 11 case-control studies in the Ovarian Cancer Association Consortium (OCAC) to examine racial/ethnic differences in epidemiological characteristics with suspected involvement in epithelial ovarian cancer (EOC) aetiology.

Methods

We used multivariable logistic regression to estimate associations for 17 reproductive, hormonal and lifestyle characteristics and EOC risk by race/ethnicity among 10 924 women with invasive EOC (8918 Non-Hispanic Whites, 433 Hispanics, 911 Blacks, 662 Asian/Pacific Islanders) and 16 150 controls (13 619 Non-Hispanic Whites, 533 Hispanics, 1233 Blacks, 765 Asian/Pacific Islanders). Likelihood ratio tests were used to evaluate heterogeneity in the risk factor associations by race/ethnicity.

Results

We observed statistically significant racial/ethnic heterogeneity for hysterectomy and EOC risk (P = 0.008), where the largest odds ratio (OR) was observed in Black women [OR = 1.64, 95% confidence interval (CI) = 1.34-2.02] compared with other racial/ethnic groups. Although not statistically significant, the associations for parity, first-degree family history of ovarian or breast cancer, and endometriosis varied by race/ethnicity. Asian/Pacific Islanders had the greatest magnitude of association for parity (≥3 births: OR = 0.38, 95% CI = 0.28-0.54), and Black women had the largest ORs for family history (OR = 1.77, 95% CI = 1.42-2.21) and endometriosis (OR = 2.42, 95% CI = 1.65-3.55).

Conclusions

Although racial/ethnic heterogeneity was observed for hysterectomy, our findings support the validity of EOC risk factors across all racial/ethnic groups, and further suggest that any racial/ethnic population with a higher prevalence of a modifiable risk factor should be targeted to disseminate information about prevention.

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.

Abstract 2787: Identifying predictive markers of endocrine response in high-grade serous ovarian cancer using RNA sequencing

Epithelial ovarian cancer is a complex disease and patients vary considerably in response to treatment. A subset of patients responds well to endocrine therapy, however hormone receptor positivity does not predict response and there are no biomarkers in clinical use to help select patients that would benefit from endocrine agents. We aimed to identify markers of endocrine response in high-grade serous ovarian cancer (HGSOC) by analysis of genes that are differentially expressed in tumors from patients that are endocrine responders compared with non-responders, combined with analysis of estrogen (E2) regulated genes that differ between E2 sensitive and E2 insensitive HGSOC cell lines.

Methods: Among women recruited to the Australian Ovarian Cancer Study, we identified 10 HGSOC patients treated with endocrine therapy (5 responders and 5 non-responders, based on GCIG CA125 response criteria) and performed RNAseq on cryopreserved tumor tissue. We also utilized paired E2 receptor (ER) positive HGSOC cell lines, derived from the same patient, with differential sensitivity to the growth effect of E2, PEO1 (E2 insensitive) and PEO4 (E2 sensitive). RNA was extracted following treatment with 0.1 nM E2 (or vehicle) for 24 hrs and gene expression changes in response to E2 were determined using RNAseq. Expression of selected ER responsive genes was validated using PCR array. Differential gene expression was determined using EdgeR. Pathway and Gene Ontology enrichment analysis of differentially expressed genes were performed using Metacore.

Results: Between endocrine responders and non-responders, 27 genes were significantly differentially expressed (False Discovery Rate (FDR)&lt;0.05, absolute log2 Fold-Change (FC)≥1). In the paired HGSOC cell lines, significantly more genes were regulated by E2 in E2-sensitive PEO4, compared with E2-insensitive PEO1 (896 genes compared to 56 genes, respectively, FDR&lt;0.05, absolute log2FC≥1). In response to E2, up-regulation of NRP1 in PEO1, up-regulation of LTBP1, MYC, PDZK1 and down-regulation of S100A6 and CYP1A1 in PEO4 were validated by PCR array profiling of ER-related genes. Genes regulated by E2 in cell lines and also differentially expressed between endocrine responders and non-responders included anoctamin 1 (ANO1) and NOTCH-Regulated Ankyrin Repeat Protein (NRARP). Both were up-regulated in endocrine responders and also up-regulated by E2 in PEO4.

Conclusion: The use of RNAseq in both tumors and cell line models was able to identify plausible sets of genes with prior evidence of association with endocrine sensitivity. In breast cancer, ANO1 has been associated with good prognosis following tamoxifen treatment, while NOTCH signaling, implicated in endocrine resistance, is negatively regulated by NRARP. Further validation is required to test the utility of these genes as predictive biomarkers that may aid in identifying HGSOC patients likely to respond to endocrine treatment.

Citation Format: Cristina Mapagu, Sian Fereday, Australian Ovarian Cancer Study Group, David D. Bowtell, Paul R. Harnett, Anna deFazio. Identifying predictive markers of endocrine response in high-grade serous ovarian cancer using RNA sequencing [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 2787. doi:10.1158/1538-7445.AM2017-2787

Germline polymorphisms in an enhancer of PSIP1 are associated with progression-free survival in epithelial ovarian cancer

Women with epithelial ovarian cancer (EOC) are usually treated with platinum/taxane therapy after cytoreductive surgery but there is considerable inter-individual variation in response. To identify germline single-nucleotide polymorphisms (SNPs) that contribute to variations in individual responses to chemotherapy, we carried out a multi-phase genome-wide association study (GWAS) in 1,244 women diagnosed with serous EOC who were treated with the same first-line chemotherapy, carboplatin and paclitaxel. We identified two SNPs (rs7874043 and rs72700653) in TTC39B (best P=7x10-5, HR=1.90, for rs7874043) associated with progression-free survival (PFS). Functional analyses show that both SNPs lie in a putative regulatory element (PRE) that physically interacts with the promoters of PSIP1, CCDC171 and an alternative promoter of TTC39B. The C allele of rs7874043 is associated with poor PFS and showed increased binding of the Sp1 transcription factor, which is critical for chromatin interactions with PSIP1. Silencing of PSIP1 significantly impaired DNA damage-induced Rad51 nuclear foci and reduced cell viability in ovarian cancer lines. PSIP1 (PC4 and SFRS1 Interacting Protein 1) is known to protect cells from stress-induced apoptosis, and high expression is associated with poor PFS in EOC patients. We therefore suggest that the minor allele of rs7874043 confers poor PFS by increasing PSIP1 expression.

Comprehensive genetic assessment of the ESR1 locus identifies a risk region for endometrial cancer

Excessive exposure to estrogen is a well-established risk factor for endometrial cancer (EC), particularly for cancers of endometrioid histology. The physiological function of estrogen is primarily mediated by estrogen receptor alpha, encoded by ESR1. Consequently, several studies have investigated whether variation at the ESR1 locus is associated with risk of EC, with conflicting results. We performed comprehensive fine-mapping analyses of 3633 genotyped and imputed single nucleotide polymorphisms (SNPs) in 6607 EC cases and 37 925 controls. There was evidence of an EC risk signal located at a potential alternative promoter of the ESR1 gene (lead SNP rs79575945, P=1.86×10(-5)), which was stronger for cancers of endometrioid subtype (P=3.76×10(-6)). Bioinformatic analysis suggests that this risk signal is in a functionally important region targeting ESR1, and eQTL analysis found that rs79575945 was associated with expression of SYNE1, a neighbouring gene. In summary, we have identified a single EC risk signal located at ESR1, at study-wide significance. Given SNPs located at this locus have been associated with risk for breast cancer, also a hormonally driven cancer, this study adds weight to the rationale for performing informed candidate fine-scale genetic studies across cancer types.

Investigation of gene-environment interactions between 47 newly identified breast cancer susceptibility loci and environmental risk factors

A large genotyping project within the Breast Cancer Association Consortium (BCAC) recently identified 41 associations between single nucleotide polymorphisms (SNPs) and overall breast cancer (BC) risk. We investigated whether the effects of these 41 SNPs, as well as six SNPs associated with estrogen receptor (ER) negative BC risk are modified by 13 environmental risk factors for BC. Data from 22 studies participating in BCAC were pooled, comprising up to 26,633 cases and 30,119 controls. Interactions between SNPs and environmental factors were evaluated using an empirical Bayes-type shrinkage estimator. Six SNPs showed interactions with associated p-values (pint ) <1.1 × 10(-3) . None of the observed interactions was significant after accounting for multiple testing. The Bayesian False Discovery Probability was used to rank the findings, which indicated three interactions as being noteworthy at 1% prior probability of interaction. SNP rs6828523 was associated with increased ER-negative BC risk in women ≥170 cm (OR = 1.22, p = 0.017), but inversely associated with ER-negative BC risk in women <160 cm (OR = 0.83, p = 0.039, pint = 1.9 × 10(-4) ). The inverse association between rs4808801 and overall BC risk was stronger for women who had had four or more pregnancies (OR = 0.85, p = 2.0 × 10(-4) ), and absent in women who had had just one (OR = 0.96, p = 0.19, pint = 6.1 × 10(-4) ). SNP rs11242675 was inversely associated with overall BC risk in never/former smokers (OR = 0.93, p = 2.8 × 10(-5) ), but no association was observed in current smokers (OR = 1.07, p = 0.14, pint = 3.4 × 10(-4) ). In conclusion, recently identified BC susceptibility loci are not strongly modified by established risk factors and the observed potential interactions require confirmation in independent studies.

MicroRNA related polymorphisms and breast cancer risk

Genetic variations, such as single nucleotide polymorphisms (SNPs) in microRNAs (miRNA) or in the miRNA binding sites may affect the miRNA dependent gene expression regulation, which has been implicated in various cancers, including breast cancer, and may alter individual susceptibility to cancer. We investigated associations between miRNA related SNPs and breast cancer risk. First we evaluated 2,196 SNPs in a case-control study combining nine genome wide association studies (GWAS). Second, we further investigated 42 SNPs with suggestive evidence for association using 41,785 cases and 41,880 controls from 41 studies included in the Breast Cancer Association Consortium (BCAC). Combining the GWAS and BCAC data within a meta-analysis, we estimated main effects on breast cancer risk as well as risks for estrogen receptor (ER) and age defined subgroups. Five miRNA binding site SNPs associated significantly with breast cancer risk: rs1045494 (odds ratio (OR) 0.92; 95% confidence interval (CI): 0.88-0.96), rs1052532 (OR 0.97; 95% CI: 0.95-0.99), rs10719 (OR 0.97; 95% CI: 0.94-0.99), rs4687554 (OR 0.97; 95% CI: 0.95-0.99, and rs3134615 (OR 1.03; 95% CI: 1.01-1.05) located in the 3' UTR of CASP8, HDDC3, DROSHA, MUSTN1, and MYCL1, respectively. DROSHA belongs to miRNA machinery genes and has a central role in initial miRNA processing. The remaining genes are involved in different molecular functions, including apoptosis and gene expression regulation. Further studies are warranted to elucidate whether the miRNA binding site SNPs are the causative variants for the observed risk effects.

A large-scale assessment of two-way SNP interactions in breast cancer susceptibility using 46,450 cases and 42,461 controls from the breast cancer association consortium

Part of the substantial unexplained familial aggregation of breast cancer may be due to interactions between common variants, but few studies have had adequate statistical power to detect interactions of realistic magnitude. We aimed to assess all two-way interactions in breast cancer susceptibility between 70,917 single nucleotide polymorphisms (SNPs) selected primarily based on prior evidence of a marginal effect. Thirty-eight international studies contributed data for 46,450 breast cancer cases and 42,461 controls of European origin as part of a multi-consortium project (COGS). First, SNPs were preselected based on evidence (P < 0.01) of a per-allele main effect, and all two-way combinations of those were evaluated by a per-allele (1 d.f.) test for interaction using logistic regression. Second, all 2.5 billion possible two-SNP combinations were evaluated using Boolean operation-based screening and testing, and SNP pairs with the strongest evidence of interaction (P < 10(-4)) were selected for more careful assessment by logistic regression. Under the first approach, 3277 SNPs were preselected, but an evaluation of all possible two-SNP combinations (1 d.f.) identified no interactions at P < 10(-8). Results from the second analytic approach were consistent with those from the first (P > 10(-10)). In summary, we observed little evidence of two-way SNP interactions in breast cancer susceptibility, despite the large number of SNPs with potential marginal effects considered and the very large sample size. This finding may have important implications for risk prediction, simplifying the modelling required. Further comprehensive, large-scale genome-wide interaction studies may identify novel interacting loci if the inherent logistic and computational challenges can be overcome.

ABCB1 (MDR1) polymorphisms and ovarian cancer progression and survival: a comprehensive analysis from the Ovarian Cancer Association Consortium and The Cancer Genome Atlas

OBJECTIVE

ABCB1 encodes the multi-drug efflux pump P-glycoprotein (P-gp) and has been implicated in multi-drug resistance. We comprehensively evaluated this gene and flanking regions for an association with clinical outcome in epithelial ovarian cancer (EOC).

METHODS

The best candidates from fine-mapping analysis of 21 ABCB1 SNPs tagging C1236T (rs1128503), G2677T/A (rs2032582), and C3435T (rs1045642) were analysed in 4616 European invasive EOC patients from thirteen Ovarian Cancer Association Consortium (OCAC) studies and The Cancer Genome Atlas (TCGA). Additionally we analysed 1,562 imputed SNPs around ABCB1 in patients receiving cytoreductive surgery and either 'standard' first-line paclitaxel-carboplatin chemotherapy (n=1158) or any first-line chemotherapy regimen (n=2867). We also evaluated ABCB1 expression in primary tumours from 143 EOC patients.

RESULT

Fine-mapping revealed that rs1128503, rs2032582, and rs1045642 were the best candidates in optimally debulked patients. However, we observed no significant association between any SNP and either progression-free survival or overall survival in analysis of data from 14 studies. There was a marginal association between rs1128503 and overall survival in patients with nil residual disease (HR 0.88, 95% CI 0.77-1.01; p=0.07). In contrast, ABCB1 expression in the primary tumour may confer worse prognosis in patients with sub-optimally debulked tumours.

CONCLUSION

Our study represents the largest analysis of ABCB1 SNPs and EOC progression and survival to date, but has not identified additional signals, or validated reported associations with progression-free survival for rs1128503, rs2032582, and rs1045642. However, we cannot rule out the possibility of a subtle effect of rs1128503, or other SNPs linked to it, on overall survival.

Estimating single nucleotide polymorphism associations using pedigree data: applications to breast cancer

BACKGROUND

Pedigrees with multiple genotyped family members have been underutilised in breast cancer (BC) genetic-association studies. We developed a pedigree-based analytical framework to characterise single-nucleotide polymorphism (SNP) associations with BC risk using data from 736 BC families ascertained through multiple affected individuals. On average, eight family members had been genotyped for 24 SNPs previously associated with BC.

METHODS

Breast cancer incidence was modelled on the basis of SNP effects and residual polygenic effects. Relative risk (RR) estimates were obtained by maximising the retrospective likelihood (RL) of observing the family genotypes conditional on all disease phenotypes. Models were extended to assess parent-of-origin effects (POEs).

RESULTS

Thirteen SNPs were significantly associated with BC under the pedigree RL approach. This approach yielded estimates consistent with those from large population-based studies. Logistic regression models ignoring pedigree structure generally gave larger RRs and association P-values. SNP rs3817198 in LSP1, previously shown to exhibit POE, yielded maternal and paternal RR estimates that were similar to those previously reported (paternal RR=1.12 (95% confidence interval (CI): 0.99-1.27), P=0.081, one-sided P=0.04; maternal RR=0.94 (95% CI: 0.84-1.06), P=0.33). No other SNP exhibited POE.

CONCLUSION

Our pedigree-based methods provide a valuable and efficient tool for characterising genetic associations with BC risk or other diseases and can complement population-based studies.

Cigarette smoking and risk of ovarian cancer: a pooled analysis of 21 case-control studies

PURPOSE

The majority of previous studies have observed an increased risk of mucinous ovarian tumors associated with cigarette smoking, but the association with other histological types is unclear. In a large pooled analysis, we examined the risk of epithelial ovarian cancer associated with multiple measures of cigarette smoking with a focus on characterizing risks according to tumor behavior and histology.

METHODS

We used data from 21 case-control studies of ovarian cancer (19,066 controls, 11,972 invasive and 2,752 borderline cases). Study-specific odds ratios (OR) and 95 % confidence intervals (CI) were obtained from logistic regression models and combined into a pooled odds ratio using a random effects model.

RESULTS

Current cigarette smoking increased the risk of invasive mucinous (OR = 1.31; 95 % CI: 1.03-1.65) and borderline mucinous ovarian tumors (OR = 1.83; 95 % CI: 1.39-2.41), while former smoking increased the risk of borderline serous ovarian tumors (OR = 1.30; 95 % CI: 1.12-1.50). For these histological types, consistent dose-response associations were observed. No convincing associations between smoking and risk of invasive serous and endometrioid ovarian cancer were observed, while our results provided some evidence of a decreased risk of invasive clear cell ovarian cancer.

CONCLUSIONS

Our results revealed marked differences in the risk profiles of histological types of ovarian cancer with regard to cigarette smoking, although the magnitude of the observed associations was modest. Our findings, which may reflect different etiologies of the histological types, add to the fact that ovarian cancer is a heterogeneous disease.

Epigenetic analysis leads to identification of HNF1B as a subtype-specific susceptibility gene for ovarian cancer

HNF1B is overexpressed in clear cell epithelial ovarian cancer, and we observed epigenetic silencing in serous epithelial ovarian cancer, leading us to hypothesize that variation in this gene differentially associates with epithelial ovarian cancer risk according to histological subtype. Here we comprehensively map variation in HNF1B with respect to epithelial ovarian cancer risk and analyse DNA methylation and expression profiles across histological subtypes. Different single-nucleotide polymorphisms associate with invasive serous (rs7405776 odds ratio (OR)=1.13, P=3.1 × 10(-10)) and clear cell (rs11651755 OR=0.77, P=1.6 × 10(-8)) epithelial ovarian cancer. Risk alleles for the serous subtype associate with higher HNF1B-promoter methylation in these tumours. Unmethylated, expressed HNF1B, primarily present in clear cell tumours, coincides with a CpG island methylator phenotype affecting numerous other promoters throughout the genome. Different variants in HNF1B associate with risk of serous and clear cell epithelial ovarian cancer; DNA methylation and expression patterns are also notably distinct between these subtypes. These findings underscore distinct mechanisms driving different epithelial ovarian cancer histological subtypes.

Comparison of 6q25 breast cancer hits from Asian and European Genome Wide Association Studies in the Breast Cancer Association Consortium (BCAC)

The 6q25.1 locus was first identified via a genome-wide association study (GWAS) in Chinese women and marked by single nucleotide polymorphism (SNP) rs2046210, approximately 180 Kb upstream of ESR1. There have been conflicting reports about the association of this locus with breast cancer in Europeans, and a GWAS in Europeans identified a different SNP, tagged here by rs12662670. We examined the associations of both SNPs in up to 61,689 cases and 58,822 controls from forty-four studies collaborating in the Breast Cancer Association Consortium, of which four studies were of Asian and 39 of European descent. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI). Case-only analyses were used to compare SNP effects in Estrogen Receptor positive (ER+) versus negative (ER-) tumours. Models including both SNPs were fitted to investigate whether the SNP effects were independent. Both SNPs are significantly associated with breast cancer risk in both ethnic groups. Per-allele ORs are higher in Asian than in European studies [rs2046210: OR (A/G) = 1.36 (95% CI 1.26-1.48), p = 7.6 × 10(-14) in Asians and 1.09 (95% CI 1.07-1.11), p = 6.8 × 10(-18) in Europeans. rs12662670: OR (G/T) = 1.29 (95% CI 1.19-1.41), p = 1.2 × 10(-9) in Asians and 1.12 (95% CI 1.08-1.17), p = 3.8 × 10(-9) in Europeans]. SNP rs2046210 is associated with a significantly greater risk of ER- than ER+ tumours in Europeans [OR (ER-) = 1.20 (95% CI 1.15-1.25), p = 1.8 × 10(-17) versus OR (ER+) = 1.07 (95% CI 1.04-1.1), p = 1.3 × 10(-7), p(heterogeneity) = 5.1 × 10(-6)]. In these Asian studies, by contrast, there is no clear evidence of a differential association by tumour receptor status. Each SNP is associated with risk after adjustment for the other SNP. These results suggest the presence of two variants at 6q25.1 each independently associated with breast cancer risk in Asians and in Europeans. Of these two, the one tagged by rs2046210 is associated with a greater risk of ER- tumours.

Tandem duplication of chromosomal segments is common in ovarian and breast cancer genomes

The application of paired-end next generation sequencing approaches has made it possible to systematically characterize rearrangements of the cancer genome to base-pair level. Utilizing this approach, we report the first detailed analysis of ovarian cancer rearrangements, comparing high-grade serous and clear cell cancers, and these histotypes with other solid cancers. Somatic rearrangements were systematically characterized in eight high-grade serous and five clear cell ovarian cancer genomes and we report here the identification of > 600 somatic rearrangements. Recurrent rearrangements of the transcriptional regulator gene, TSHZ3, were found in three of eight serous cases. Comparison to breast, pancreatic and prostate cancer genomes revealed that a subset of ovarian cancers share a marked tandem duplication phenotype with triple-negative breast cancers. The tandem duplication phenotype was not linked to BRCA1/2 mutation, suggesting that other common mechanisms or carcinogenic exposures are operative. High-grade serous cancers arising in women with germline BRCA1 or BRCA2 mutation showed a high frequency of small chromosomal deletions. These findings indicate that BRCA1/2 germline mutation may contribute to widespread structural change and that other undefined mechanism(s), which are potentially shared with triple-negative breast cancer, promote tandem chromosomal duplications that sculpt the ovarian cancer genome.

MicroRNA genes and their target 3'-untranslated regions are infrequently somatically mutated in ovarian cancers

MicroRNAs are key regulators of gene expression and have been shown to have altered expression in a variety of cancer types, including epithelial ovarian cancer. MiRNA function is most often achieved through binding to the 3′-untranslated region of the target protein coding gene. Mutation screening using massively-parallel sequencing of 712 miRNA genes in 86 ovarian cancer cases identified only 5 mutated miRNA genes, each in a different case. One mutation was located in the mature miRNA, and three mutations were predicted to alter the secondary structure of the miRNA transcript. Screening of the 3′-untranslated region of 18 candidate cancer genes identified one mutation in each of AKT2, EGFR, ERRB2 and CTNNB1. The functional effect of these mutations is unclear, as expression data available for AKT2 and EGFR showed no increase in gene transcript. Mutations in miRNA genes and 3′-untranslated regions are thus uncommon in ovarian cancer.

Polymorphisms in stromal genes and susceptibility to serous epithelial ovarian cancer: a report from the Ovarian Cancer Association Consortium

Alterations in stromal tissue components can inhibit or promote epithelial tumorigenesis. Decorin (DCN) and lumican (LUM) show reduced stromal expression in serous epithelial ovarian cancer (sEOC). We hypothesized that common variants in these genes associate with risk. Associations with sEOC among Caucasians were estimated with odds ratios (OR) among 397 cases and 920 controls in two U.S.-based studies (discovery set), 436 cases and 1,098 controls in Australia (replication set 1) and a consortium of 15 studies comprising 1,668 cases and 4,249 controls (replication set 2). The discovery set and replication set 1 (833 cases and 2,013 controls) showed statistically homogeneous (P(heterogeneity)≥0.48) decreased risks of sEOC at four variants: DCN rs3138165, rs13312816 and rs516115, and LUM rs17018765 (OR = 0.6 to 0.9; P(trend) = 0.001 to 0.03). Results from replication set 2 were statistically homogeneous (P(heterogeneity)≥0.13) and associated with increased risks at DCN rs3138165 and rs13312816, and LUM rs17018765: all ORs = 1.2; P(trend)≤0.02. The ORs at the four variants were statistically heterogeneous across all 18 studies (P(heterogeneity)≤0.03), which precluded combining. In post-hoc analyses, interactions were observed between each variant and recruitment period (P(interaction)≤0.003), age at diagnosis (P(interaction) = 0.04), and year of diagnosis (P(interaction) = 0.05) in the five studies with available information (1,044 cases, 2,469 controls). We conclude that variants in DCN and LUM are not directly associated with sEOC, and that confirmation of possible effect modification of the variants by non-genetic factors is required.

Deregulation of MYCN, LIN28B and LET7 in a molecular subtype of aggressive high-grade serous ovarian cancers

Molecular subtypes of serous ovarian cancer have been recently described. Using data from independent datasets including over 900 primary tumour samples, we show that deregulation of the Let-7 pathway is specifically associated with the C5 molecular subtype of serous ovarian cancer. DNA copy number and gene expression of HMGA2, alleles of Let-7, LIN28, LIN28B, MYC, MYCN, DICER1, and RNASEN were measured using microarray and quantitative reverse transcriptase PCR. Immunohistochemistry was performed on 127 samples using tissue microarrays and anti-HMGA2 antibodies. Fluorescence in situ hybridisation of bacterial artificial chromosomes hybridized to 239 ovarian tumours was used to measure translocation at the LIN28B locus. Short interfering RNA knockdown in ovarian cell lines was used to test the functionality of associations observed. Four molecular subtypes (C1, C2, C4, C5) of high-grade serous ovarian cancers were robustly represented in each dataset and showed similar pattern of patient survival. We found highly specific activation of a pathway involving MYCN, LIN28B, Let-7 and HMGA2 in the C5 molecular subtype defined by MYCN amplification and over-expression, over-expression of MYCN targets including the Let-7 repressor LIN28B, loss of Let-7 expression and HMGA2 amplification and over-expression. DICER1, a known Let-7 target, and RNASEN were over-expressed in C5 tumours. We saw no evidence of translocation at the LIN28B locus in C5 tumours. The reported interaction between LIN28B and Let-7 was recapitulated by siRNA knockdown in ovarian cancer cell lines. Our results associate deregulation of MYCN and downstream targets, including Let-7 and oncofetal genes, with serous ovarian cancer. We define for the first time how elements of an oncogenic pathway, involving multiple genes that contribute to stem cell renewal, is specifically altered in a molecular subtype of serous ovarian cancer. By defining the drivers of a molecular subtype of serous ovarian cancers we provide a novel strategy for targeted therapeutic intervention.

Amplicon-dependent CCNE1 expression is critical for clonogenic survival after cisplatin treatment and is correlated with 20q11 gain in ovarian cancer

Genomic amplification of 19q12 occurs in several cancer types including ovarian cancer where it is associated with primary treatment failure. We systematically attenuated expression of genes within the minimally defined 19q12 region in ovarian cell lines using short-interfering RNAs (siRNA) to identify driver oncogene(s) within the amplicon. Knockdown of CCNE1 resulted in G1/S phase arrest, reduced cell viability and apoptosis only in amplification-carrying cells. Although CCNE1 knockdown increased cisplatin resistance in short-term assays, clonogenic survival was inhibited after treatment. Gain of 20q11 was highly correlated with 19q12 amplification and spanned a 2.5 Mb region including TPX2, a centromeric protein required for mitotic spindle function. Expression of TPX2 was highly correlated with gene amplification and with CCNE1 expression in primary tumors. siRNA inhibition of TPX2 reduced cell viability but this effect was not amplicon-dependent. These findings demonstrate that CCNE1 is a key driver in the 19q12 amplicon required for survival and clonogenicity in cells with locus amplification. Co-amplification at 19q12 and 20q11 implies the presence of a cooperative mutational network. These observations have implications for the application of targeted therapies in CCNE1 dependent ovarian cancers.

Abstract 4727: TERT SNPs and ovarian cancer risk

There is increasing evidence that neighboring stromal cells are not uninvolved bystanders in the pathway to cancer. We hypothesized that SNPs in genes expressed as a consequence of interactions between ovarian cancer cells and the host micro-environment could contribute to ovarian cancer susceptibility. We therefore used a two-stage approach to comprehensively evaluate common single nucleotide polymorphisms (SNPs) in the Ovarian Cancer Association Consortium (OCAC) for 173 genes involved in stromal-epithelial interactions. Cases with epithelial ovarian cancer (n=675) and controls (n=1,162) were genotyped at 1536 SNPs using an Illumina GoldenGate™ BeadArray array. Logistic regression was used to estimate odds ratios (ORs) and adjust for potential confounders. Based on Positive Predictive Value estimates, three SNPs - PODXL rs1013368, ITGA6 rs13027811 and MMP3 rs522616 - were selected for replication using TaqMan genotyping in 3,734 serous invasive cases and 10,067 controls from 16 OCAC case-control studies. An additional 18 SNPs with Pper-allele &lt;0.05 in the test set, and which fitted into an iPLEX design, were selected for replication by a subset of five OCAC studies (n=1771 serous invasive cases; n=4022 controls). The associations observed in the test set for the PODXL, ITGA6 and MMP3 SNPs were completely attenuated in the replication set (adj. Pper-allele ≥0.5). However, TERT rs7726159 showed evidence of a log additive association with ovarian cancer risk in the replication study (Pper-allele =0.03). Combined analysis of the test and replication sets for this SNP showed a significant association with risk of serous ovarian cancer among non-Hispanic whites [adj. ORper-allele 1.14 (1.04 - 1.24) p = 0.003]. A stronger association was observed for serous invasive cases when data for all ethnicities were included [adj. ORper-allele 1.18 (1.09 - 1.28) p = 2.46 × 10−5]. Stratification by histological subtype and malignant potential showed that TERT rs7726159 was also associated with risk of serous borderline tumors, but not with any other subtypes. Although TERT variants have not been previously reported to be associated with ovarian cancer, multiple genome-wide association studies have reported associations with this locus and risk of glioma, basal cell carcinoma, and several cancer types that are known to originate in the epithelium, including lung, bladder, prostate and cervical cancer. Further analysis of this locus is necessary to definitively identify the causal SNP(s). Our study adds to the growing evidence that, as well as 8q24, the TERT locus at 5p15.33, is a general cancer susceptibility locus.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 4727.

Polymorphism in the IL18 gene and epithelial ovarian cancer in non-Hispanic white women

Over 22,000 cases of ovarian cancer were diagnosed in 2007 in the United States, but only a fraction of them can be attributed to mutations in highly penetrant genes such as BRCA1. To determine whether low-penetrance genetic variants contribute to ovarian cancer risk, we genotyped 1,536 single nucleotide polymorphisms (SNP) in several candidate gene pathways in 848 epithelial ovarian cancer cases and 798 controls in the North Carolina Ovarian Cancer Study (NCO) using a customized Illumina array. The inflammation gene interleukin-18 (IL18) showed the strongest evidence for association with epithelial ovarian cancer in a gene-by-gene analysis (P = 0.002) with a <25% chance of being a false-positive finding (q value = 0.240). Using a multivariate model search algorithm over 11 IL18 tagging SNPs, we found that the association was best modeled by rs1834481. Further, this SNP uniquely tagged a significantly associated IL18 haplotype and there was an increased risk of epithelial ovarian cancer per rs1834481 allele (odds ratio, 1.24; 95% confidence interval, 1.06-1.45). In a replication stage, 12 independent studies from the Ovarian Cancer Association Consortium (OCAC) genotyped rs1834481 in an additional 5,877 cases and 7,791 controls. The fixed effects estimate per rs1834481 allele was null (odds ratio, 0.99; 95% confidence interval, 0.94-1.05) when data from the 12 OCAC studies were combined. The effect estimate remained unchanged with the addition of the initial North Carolina Ovarian Cancer Study data. This analysis shows the importance of consortia, like the OCAC, in either confirming or refuting the validity of putative findings in studies with smaller sample sizes. (Cancer Epidemiol Biomarkers Prev 2008;17(12):3567-72).

High resolution melting for mutation scanning of TP53 exons 5-8

BACKGROUND

p53 is commonly inactivated by mutations in the DNA-binding domain in a wide range of cancers. As mutant p53 often influences response to therapy, effective and rapid methods to scan for mutations in TP53 are likely to be of clinical value. We therefore evaluated the use of high resolution melting (HRM) as a rapid mutation scanning tool for TP53 in tumour samples.

METHODS

We designed PCR amplicons for HRM mutation scanning of TP53 exons 5 to 8 and tested them with DNA from cell lines hemizygous or homozygous for known mutations. We assessed the sensitivity of each PCR amplicon using dilutions of cell line DNA in normal wild-type DNA. We then performed a blinded assessment on ovarian tumour DNA samples that had been previously sequenced for mutations in TP53 to assess the sensitivity and positive predictive value of the HRM technique. We also performed HRM analysis on breast tumour DNA samples with unknown TP53 mutation status.

RESULTS

One cell line mutation was not readily observed when exon 5 was amplified. As exon 5 contained multiple melting domains, we divided the exon into two amplicons for further screening. Sequence changes were also introduced into some of the primers to improve the melting characteristics of the amplicon. Aberrant HRM curves indicative of TP53 mutations were observed for each of the samples in the ovarian tumour DNA panel. Comparison of the HRM results with the sequencing results revealed that each mutation was detected by HRM in the correct exon. For the breast tumour panel, we detected seven aberrant melt profiles by HRM and subsequent sequencing confirmed the presence of these and no other mutations in the predicted exons.

CONCLUSION

HRM is an effective technique for simple and rapid scanning of TP53 mutations that can markedly reduce the amount of sequencing required in mutational studies of TP53.