Jakubowska, Anna

Background: Breast cancer polygenic risk scores (PRS) and traditional risk models (e.g., the Gail model [Gail]) are known to contribute largely independent information, but it is unclear how the overlap varies by ancestry, age, disease type (invasive breast cancer, DCIS), and risk threshold. Methods: In a retrospective case–control study, we evaluated risk prediction performance in 180,398 women (161,849 of European ancestry; 18,549 of Asian ancestry). Odds ratios (ORs) from logistic regression models and the area under the receiver operating characteristic curve (AUC) were estimated. Results: PRS for invasive disease showed a stronger association in younger (<50 years) women (OR = 2.51, AUC = 0.622) than in women ≥ 50 years (OR = 2.06, AUC = 0.653) of European ancestry. PRS performance in Asians was lower (OR range = 1.62–1.64, AUC = 0.551–0.600). Gail performance was modest across groups and poor in younger Asian women (OR = 0.94–0.99, AUC = 0.523–0.533). Age interactions were observed for both PRS (p < 0.001) and Gail (p < 0.001) in Europeans, whereas in Asians, age interaction was observed only for Gail (invasive: p < 0.001; DCIS: p = 0.002). PRS identified more high-risk individuals than Gail in Asian populations, especially ≥50 years, while Gail identified more in Europeans. Overlap between PRS, Gail, and family history was limited at higher thresholds. Calibration analysis, comparing empirical and model-based ROC curves, showed divergence for both PRS and Gail (p < 0.001), which indicates miscalibration. In Europeans, family history and prior biopsies drove Gail discrimination. In younger Asians, age at first live birth was influential. Conclusions: PRS adds value to risk stratification beyond traditional tools, especially in younger women and Asian ancestry populations.

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.

Abstract Background: It is not known whether modifiable lifestyle factors that predict survival after invasive breast cancer differ by subtype. Methods: We analyzed data for 121,435 women diagnosed with breast cancer from 67 studies in the Breast Cancer Association Consortium with 16,890 deaths (8,554 breast cancer specific) over 10 years. Cox regression was used to estimate associations between risk factors and 10-year all-cause mortality and breast cancer–specific mortality overall, by estrogen receptor (ER) status, and by intrinsic-like subtype. Results: There was no evidence of heterogeneous associations between risk factors and mortality by subtype (Padj &gt; 0.30). The strongest associations were between all-cause mortality and BMI ≥30 versus 18.5–25 kg/m2 [HR (95% confidence interval (CI), 1.19 (1.06–1.34)]; current versus never smoking [1.37 (1.27–1.47)], high versus low physical activity [0.43 (0.21–0.86)], age ≥30 years versus &lt;20 years at first pregnancy [0.79 (0.72–0.86)]; &gt;0–&lt;5 years versus ≥10 years since last full-term birth [1.31 (1.11–1.55)]; ever versus never use of oral contraceptives [0.91 (0.87–0.96)]; ever versus never use of menopausal hormone therapy, including current estrogen–progestin therapy [0.61 (0.54–0.69)]. Similar associations with breast cancer mortality were weaker; for example, 1.11 (1.02–1.21) for current versus never smoking. Conclusions: We confirm associations between modifiable lifestyle factors and 10-year all-cause mortality. There was no strong evidence that associations differed by ER status or intrinsic-like subtype. Impact: Given the large dataset and lack of evidence that associations between modifiable risk factors and 10-year mortality differed by subtype, these associations could be cautiously used in prognostication models to inform patient-centered care.

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.

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 Genome-wide association studies have identified 40 ovarian cancer risk loci. However, the mechanisms underlying these associations remain elusive. In this study, we conducted a two-pronged approach to identify candidate causal SNPs and assess underlying biological mechanisms at chromosome 9p22.2, the first and most statistically significant associated locus for ovarian cancer susceptibility. Three transcriptional regulatory elements with allele-specific effects and a scaffold/matrix attachment region were characterized and, through physical DNA interactions, BNC2 was established as the most likely target gene. We determined the consensus binding sequence for BNC2 in vitro, verified its enrichment in BNC2 ChIP-seq regions, and validated a set of its downstream target genes. Fine-mapping by dense regional genotyping in over 15,000 ovarian cancer cases and 30,000 controls identified SNPs in the scaffold/matrix attachment region as among the most likely causal variants. This study reveals a comprehensive regulatory landscape at 9p22.2 and proposes a likely mechanism of susceptibility to ovarian cancer. Significance: Mapping the 9p22.2 ovarian cancer risk locus identifies BNC2 as an ovarian cancer risk gene. See related commentary by Choi and Brown, p. 439

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.

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