Gronwald, Jacek

Abstract Background: The lack of consensus on whether bilateral oophorectomy impacts risk of developing breast cancer among BRCA1 mutation carriers might be attributed to various biases, specifically, cancer-induced testing bias due to inclusion of prevalent cases. We conducted two complementary matched case–control analyses to evaluate the association of oophorectomy and BRCA1 breast cancer. Methods: A research questionnaire was administered every two years to collect information on exposures and disease. In the first analysis, we limited the study to prevalent breast cancer cases (diagnosed prior to study entry; n = 2,962) who were matched to controls on year of birth and country of residence (n = 4,358). In the second approach, we limited to 330 incident cases (diagnosed in the follow-up period) and 1,548 matched controls. Conditional logistic regression was used to estimate the adjusted odds ratios (OR) and 95% confidence intervals (CI) of invasive breast cancer. Results: In the first approach, there was a significant inverse association between oophorectomy and the risk of developing breast cancer [OR = 0.43; 95% confidence interval (CI), 0.34–0.55; P < 00001]. In the second approach, there was no association between oophorectomy and risk (OR = 1.21; 95% CI, 0.87–1.70; P = 0.26). Conclusions: The inclusion of women with a personal history of breast cancer prior to ascertainment likely impacts upon the association of oophorectomy and BRCA1 breast cancer risk. Impact: Oophorectomy is unlikely a determinant of breast cancer risk in BRCA1 mutation carriers but should be offered at age 35 to reduce the risk of ovarian and fallopian tube cancer.

Abstract ARID1A (BAF250a) is a component of the SWI/SNF chromatin modifying complex, plays an important tumour suppressor role, and is considered prognostic in several malignancies. However, in ovarian carcinomas there are contradictory reports on its relationship to outcome, immune response, and correlation with clinicopathological features. We assembled a series of 1623 endometriosis‐associated ovarian carcinomas, including 1078 endometrioid (ENOC) and 545 clear cell (CCOC) ovarian carcinomas, through combining resources of the Ovarian Tumor Tissue Analysis (OTTA) Consortium, the Canadian Ovarian Unified Experimental Resource (COEUR), local, and collaborative networks. Validated immunohistochemical surrogate assays for ARID1A mutations were applied to all samples. We investigated associations between ARID1A loss/mutation, clinical features, outcome, CD8 + tumour‐infiltrating lymphocytes (CD8 + TILs), and DNA mismatch repair deficiency (MMRd). ARID1A loss was observed in 42% of CCOCs and 25% of ENOCs. We found no associations between ARID1A loss and outcomes, stage, age, or CD8 + TIL status in CCOC. Similarly, we found no association with outcome or stage in endometrioid cases. In ENOC, ARID1A loss was more prevalent in younger patients ( p = 0.012) and was associated with MMRd ( p < 0.001) and the presence of CD8 + TILs ( p = 0.008). Consistent with MMRd being causative of ARID1A mutations, in a subset of ENOCs we also observed an association with ARID1A loss‐of‐function mutation as a result of small indels ( p = 0.035, versus single nucleotide variants). In ENOC, the association with ARID1A loss, CD8 + TILs, and age appears confounded by MMRd status. Although this observation does not explicitly rule out a role for ARID1A influence on CD8 + TIL infiltration in ENOC, given current knowledge regarding MMRd, it seems more likely that effects are dominated by the hypermutation phenotype. This large dataset with consistently applied biomarker assessment now provides a benchmark for the prevalence of ARID1A loss‐of‐function mutations in endometriosis‐associated ovarian cancers and brings clarity to the prognostic significance. © 2021 The Pathological Society of Great Britain and Ireland.

Abstract Background: Accumulating evidence suggests a relationship between endometrial cancer and ovarian cancer. Independent genome-wide association studies (GWAS) for endometrial cancer and ovarian cancer have identified 16 and 27 risk regions, respectively, four of which overlap between the two cancers. We aimed to identify joint endometrial and ovarian cancer risk loci by performing a meta-analysis of GWAS summary statistics from these two cancers. Methods: Using LDScore regression, we explored the genetic correlation between endometrial cancer and ovarian cancer. To identify loci associated with the risk of both cancers, we implemented a pipeline of statistical genetic analyses (i.e., inverse-variance meta-analysis, colocalization, and M-values) and performed analyses stratified by subtype. Candidate target genes were then prioritized using functional genomic data. Results: Genetic correlation analysis revealed significant genetic correlation between the two cancers (rG = 0.43, P = 2.66 × 10−5). We found seven loci associated with risk for both cancers (PBonferroni &lt; 2.4 × 10−9). In addition, four novel subgenome-wide regions at 7p22.2, 7q22.1, 9p12, and 11q13.3 were identified (P &lt; 5 × 10−7). Promoter-associated HiChIP chromatin loops from immortalized endometrium and ovarian cell lines and expression quantitative trait loci data highlighted candidate target genes for further investigation. Conclusions: Using cross-cancer GWAS meta-analysis, we have identified several joint endometrial and ovarian cancer risk loci and candidate target genes for future functional analysis. Impact: Our research highlights the shared genetic relationship between endometrial cancer and ovarian cancer. Further studies in larger sample sets are required to confirm our findings.

Abstract Purpose: Gene expression–based molecular subtypes of high-grade serous tubo-ovarian cancer (HGSOC), demonstrated across multiple studies, may provide improved stratification for molecularly targeted trials. However, evaluation of clinical utility has been hindered by nonstandardized methods, which are not applicable in a clinical setting. We sought to generate a clinical grade minimal gene set assay for classification of individual tumor specimens into HGSOC subtypes and confirm previously reported subtype-associated features. Experimental Design: Adopting two independent approaches, we derived and internally validated algorithms for subtype prediction using published gene expression data from 1,650 tumors. We applied resulting models to NanoString data on 3,829 HGSOCs from the Ovarian Tumor Tissue Analysis consortium. We further developed, confirmed, and validated a reduced, minimal gene set predictor, with methods suitable for a single-patient setting. Results: Gene expression data were used to derive the predictor of high-grade serous ovarian carcinoma molecular subtype (PrOTYPE) assay. We established a de facto standard as a consensus of two parallel approaches. PrOTYPE subtypes are significantly associated with age, stage, residual disease, tumor-infiltrating lymphocytes, and outcome. The locked-down clinical grade PrOTYPE test includes a model with 55 genes that predicted gene expression subtype with &gt;95% accuracy that was maintained in all analytic and biological validations. Conclusions: We validated the PrOTYPE assay following the Institute of Medicine guidelines for the development of omics-based tests. This fully defined and locked-down clinical grade assay will enable trial design with molecular subtype stratification and allow for objective assessment of the predictive value of HGSOC molecular subtypes in precision medicine applications. See related commentary by McMullen et al., p. 5271

Purpose The known epithelial ovarian cancer (EOC) susceptibility genes account for less than 50% of the heritable risk of ovarian cancer suggesting that other susceptibility genes exist. The aim of this study was to evaluate the contribution to ovarian cancer susceptibility of rare deleterious germline variants in a set of candidate genes. Methods We sequenced the coding region of 54 candidate genes in 6385 invasive EOC cases and 6115 controls of broad European ancestry. Genes with an increased frequency of putative deleterious variants in cases versus controls were further examined in an independent set of 14 135 EOC cases and 28 655 controls from the Ovarian Cancer Association Consortium and the UK Biobank. For each gene, we estimated the EOC risks and evaluated associations between germline variant status and clinical characteristics. Results The ORs associated for high-grade serous ovarian cancer were 3.01 for PALB2 (95% CI 1.59 to 5.68; p=0.00068), 1.99 for POLK (95% CI 1.15 to 3.43; p=0.014) and 4.07 for SLX4 (95% CI 1.34 to 12.4; p=0.013). Deleterious mutations in FBXO10 were associated with a reduced risk of disease (OR 0.27, 95% CI 0.07 to 1.00, p=0.049). However, based on the Bayes false discovery probability, only the association for PALB2 in high-grade serous ovarian cancer is likely to represent a true positive. Conclusions We have found strong evidence that carriers of PALB2 deleterious mutations are at increased risk of high-grade serous ovarian cancer. Whether the magnitude of risk is sufficiently high to warrant the inclusion of PALB2 in cancer gene panels for ovarian cancer risk testing is unclear; much larger sample sizes will be needed to provide sufficiently precise estimates for clinical counselling.

AbstractPrevious transcriptome‐wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome‐wide association studies (GWAS), but analyses of breast cancer subtype‐specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta‐analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER−). We further compared associations with ER+ and ER− subtypes, using a case‐only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER– breast cancer. We further identified 30 TWAS‐significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast‐cancer gene in three of six regions harboring multiple TWAS‐significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER− breast cancer.

Abstract Background The purpose of this study was to estimate precise age-specific tubo-ovarian carcinoma (TOC) and breast cancer (BC) risks for carriers of pathogenic variants in RAD51C and RAD51D. Methods We analyzed data from 6178 families, 125 with pathogenic variants in RAD51C, and 6690 families, 60 with pathogenic variants in RAD51D. TOC and BC relative and cumulative risks were estimated using complex segregation analysis to model the cancer inheritance patterns in families while adjusting for the mode of ascertainment of each family. All statistical tests were two-sided. Results Pathogenic variants in both RAD51C and RAD51D were associated with TOC (RAD51C: relative risk [RR] = 7.55, 95% confidence interval [CI] = 5.60 to 10.19; P = 5 × 10-40; RAD51D: RR = 7.60, 95% CI = 5.61 to 10.30; P = 5 × 10-39) and BC (RAD51C: RR = 1.99, 95% CI = 1.39 to 2.85; P = 1.55 × 10-4; RAD51D: RR = 1.83, 95% CI = 1.24 to 2.72; P = .002). For both RAD51C and RAD51D, there was a suggestion that the TOC relative risks increased with age until around age 60 years and decreased thereafter. The estimated cumulative risks of developing TOC to age 80 years were 11% (95% CI = 6% to 21%) for RAD51C and 13% (95% CI = 7% to 23%) for RAD51D pathogenic variant carriers. The estimated cumulative risks of developing BC to 80 years were 21% (95% CI = 15% to 29%) for RAD51C and 20% (95% CI = 14% to 28%) for RAD51D pathogenic variant carriers. Both TOC and BC risks for RAD51C and RAD51D pathogenic variant carriers varied by cancer family history and could be as high as 32–36% for TOC, for carriers with two first-degree relatives diagnosed with TOC, or 44–46% for BC, for carriers with two first-degree relatives diagnosed with BC. Conclusions These estimates will facilitate the genetic counseling of RAD51C and RAD51D pathogenic variant carriers and justify the incorporation of RAD51C and RAD51D into cancer risk prediction models.

Abstract 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 &lt; 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.

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.

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