Simon A. Gayther

Patterns of Associations with Epidemiologic Factors by High-Grade Serous Ovarian Cancer Gene Expression Subtypes

Abstract Background: Ovarian high-grade serous carcinomas (HGSC) comprise four distinct molecular subtypes based on mRNA expression patterns, with differential survival. Understanding risk factor associations is important to elucidate the etiology of HGSC. We investigated associations between different epidemiologic risk factors and HGSC molecular subtypes. Methods: We pooled data from 11 case–control studies with epidemiologic and tumor gene expression data from custom NanoString CodeSets developed through a collaboration within the Ovarian Tumor Tissue Analysis consortium. The PrOTYPE-validated NanoString-based 55-gene classifier was used to assign HGSC gene expression subtypes. We examined associations between epidemiologic factors and HGSC subtypes in 2,070 cases and 16,633 controls using multivariable-adjusted polytomous regression models. Results: Among the 2,070 HGSC cases, 556 (27%) were classified as C1.MES, 340 (16%) as C5.PRO, 538 (26%) as C2.IMM, and 636 (31%) as C4.DIF. The key factors, including oral contraceptive use, parity, breastfeeding, and family history of ovarian cancer, were similarly associated with all subtypes. Heterogeneity was observed for several factors. Former smoking [OR = 1.25; 95% confidence interval (CI) = 1.03, 1.51] and genital powder use (OR = 1.42; 95% CI = 1.08, 1.86) were uniquely associated with C2.IMM. History of endometriosis was associated with C5.PRO (OR = 1.46; 95% CI = 0.98, 2.16) and C4.DIF (OR = 1.27; 95% CI = 0.94, 1.71) only. Family history of breast cancer (OR = 1.44; 95% CI = 1.16, 1.78) and current smoking (OR = 1.40; 95% CI = 1.11, 1.76) were associated with C4.DIF only. Conclusions: This study observed heterogeneous associations of epidemiologic and modifiable factors with HGSC molecular subtypes. Impact: The different patterns of associations may provide key information about the etiology of the four subtypes.

Rare germline genetic variation in PAX8 transcription factor binding sites and susceptibility to epithelial ovarian cancer

Abstract Common genetic variation throughout the genome and rare coding variants identified to date explain about half of the inherited genetic component of epithelial ovarian cancer risk. It is likely that rare variation in the noncoding genome will explain some of the unexplained heritability, but identifying such variants is challenging. The primary problem is a lack of statistical power to identify individual risk variants by association, as power is a function of sample size, effect size, and allele frequency. Power can be increased by using burden tests, which test for the association of carriers of any variant in a specified genomic region. This has the effect of increasing the putative effect allele frequency. PAX8 is a transcription factor that plays a critical role in tumor progression, migration, and invasion. Furthermore, regulatory elements proximal to target genes of PAX8 are enriched for common ovarian cancer risk variants. We hypothesized that rare variation in PAX8 binding sites is also associated with ovarian cancer risk but unlikely to be associated with risk of breast, colorectal, or endometrial cancer. We have used publicly available, whole-genome sequencing data from the UK 100,000 Genomes Project to evaluate the burden of rare variation in PAX8 binding sites across the genome. Data were available for 522 ovarian cancers, 2984 breast cancers, 2696 colorectal cancers, 836 endometrial cancers, and 2253 noncancer controls. Active binding sites were defined using data from multiple PAX8 and H3K27 chromatin immunoprecipitation sequencing experiments. We found no association between the burden of rare variation in PAX8 binding sites (defined in several ways) and risk of ovarian, breast, or endometrial cancer. An apparent association with colorectal cancer was likely to be a technical artifact as a similar association was also detected for rare variation in random regions of the genome. Despite the null result, this study provides a proof-of-principle for using burden testing to identify rare, noncoding germline genetic variation associated with disease. Larger sample sizes available from large-scale sequencing projects, together with improved understanding of the function of the noncoding genome, will increase the potential of similar studies in the future.

Predicted Proteome Association Studies of Breast, Prostate, Ovarian, and Endometrial Cancers Implicate Plasma Protein Regulation in Cancer Susceptibility

Abstract Background: Predicting protein levels from genotypes for proteome-wide association studies (PWAS) may provide insight into the mechanisms underlying cancer susceptibility. Methods: We performed PWAS of breast, endometrial, ovarian, and prostate cancers and their subtypes in several large European-ancestry discovery consortia (effective sample size: 237,483 cases/317,006 controls) and tested the results for replication in an independent European-ancestry GWAS (31,969 cases/410,350 controls). We performed PWAS using the cancer GWAS summary statistics and two sets of plasma protein prediction models, followed by colocalization analysis. Results: Using Atherosclerosis Risk in Communities (ARIC) models, we identified 93 protein–cancer associations [false discovery rate (FDR) < 0.05]. We then performed a meta-analysis of the discovery and replication PWAS, resulting in 61 significant protein–cancer associations (FDR < 0.05). Ten of 15 protein–cancer pairs that could be tested using Trans-Omics for Precision Medicine (TOPMed) protein prediction models replicated with the same directions of effect in both cancer GWAS (P < 0.05). To further support our results, we applied Bayesian colocalization analysis and found colocalized SNPs for SERPINA3 protein levels and prostate cancer (posterior probability, PP = 0.65) and SNUPN protein levels and breast cancer (PP = 0.62). Conclusions: We used PWAS to identify potential biomarkers of hormone-related cancer risk. SNPs in SERPINA3 and SNUPN did not reach genome-wide significance for cancer in the original GWAS, highlighting the power of PWAS for novel locus discovery, with the added advantage of providing directions of protein effect. Impact: PWAS and colocalization are promising methods to identify potential molecular mechanisms underlying complex traits.

Concurrent RB1 Loss and BRCA Deficiency Predicts Enhanced Immunologic Response and Long-term Survival in Tubo-ovarian High-grade Serous Carcinoma

Abstract Purpose: The purpose of this study was to evaluate RB1 expression and survival across ovarian carcinoma histotypes and how co-occurrence of BRCA1 or BRCA2 (BRCA) alterations and RB1 loss influences survival in tubo-ovarian high-grade serous carcinoma (HGSC). Experimental Design: RB1 protein expression was classified by immunohistochemistry in ovarian carcinomas of 7,436 patients from the Ovarian Tumor Tissue Analysis consortium. We examined RB1 expression and germline BRCA status in a subset of 1,134 HGSC, and related genotype to overall survival (OS), tumor-infiltrating CD8+ lymphocytes, and transcriptomic subtypes. Using CRISPR-Cas9, we deleted RB1 in HGSC cells with and without BRCA1 alterations to model co-loss with treatment response. We performed whole-genome and transcriptome data analyses on 126 patients with primary HGSC to characterize tumors with concurrent BRCA deficiency and RB1 loss. Results: RB1 loss was associated with longer OS in HGSC but with poorer prognosis in endometrioid ovarian carcinoma. Patients with HGSC harboring both RB1 loss and pathogenic germline BRCA variants had superior OS compared with patients with either alteration alone, and their median OS was three times longer than those without pathogenic BRCA variants and retained RB1 expression (9.3 vs. 3.1 years). Enhanced sensitivity to cisplatin and paclitaxel was seen in BRCA1-altered cells with RB1 knockout. Combined RB1 loss and BRCA deficiency correlated with transcriptional markers of enhanced IFN response, cell-cycle deregulation, and reduced epithelial–mesenchymal transition. CD8+ lymphocytes were most prevalent in BRCA-deficient HGSC with co-loss of RB1. Conclusions: Co-occurrence of RB1 loss and BRCA deficiency was associated with exceptionally long survival in patients with HGSC, potentially due to better treatment response and immune stimulation.