Winham, Stacey J

Abstract Nonsteroidal anti-inflammatory agents (NSAID) are associated with modest inconsistent reductions in breast cancer risk in population-based cohorts, whereas two focused studies of patients with benign breast disease (BBD) have found lower risk with NSAID use. Given that BBD includes fibroinflammatory lesions linked to elevated breast cancer risk, we assessed whether NSAID use was associated with lower breast cancer risk among patients with BBD. Participants were postmenopausal women in the Cancer Prevention Study-II (CPS-II), a prospective study of cancer incidence and mortality, who completed follow-up surveys in 1997 with follow-up through June 30, 2015. History of BBD, NSAID use, and covariate data were updated biennially. This analysis included 23,615 patients with BBD and 36,751 patients with non-BBD, including 3,896 incident breast cancers over an average of 12.72 years of follow-up among participants. NSAID use, overall and by formulation, recency, duration, and pills per month was analyzed versus breast cancer risk overall and by BBD status using multivariable-adjusted Cox models; BBD status and NSAID use were modeled as time-dependent exposures. Patients with BBD who reported using NSAIDs experienced lower breast cancer risk (HR, 0.87; 95% CI, 0.78–0.97), with similar effects for estrogen receptor (ER)-positive breast cancers [HR, 0.85; 95% confidence interval (CI), 0.74–0.97] and ER-negative breast cancers (HR, 0.87; 95% CI, 0.59–1.29); among women without BBD, NSAID use was unrelated to risk (HR, 1.02; 95% CI, 0.92–1.13; Pinteraction = 0.04). Associations stratified by age, obesity, menopausal hormone use, and cardiovascular disease were similar. Among patients with BBD, NSAID use appears linked to lower breast cancer risk. Further studies to assess the value of NSAID use among patients with BBD are warranted. Prevention Relevance: We examined whether NSAID use, a modifiable exposure, is associated with breast cancer risk in postmenopausal women from the Cancer Prevention Study-II with self-reported benign breast disease, an often inflammatory condition associated with higher rates of breast cancer.

Abstract Purpose: To identify molecular subclasses of clear cell ovarian carcinoma (CCOC) and assess their impact on clinical presentation and outcomes. Experimental Design: We profiled 421 primary CCOCs that passed quality control using a targeted deep sequencing panel of 163 putative CCOC driver genes and whole transcriptome sequencing of 211 of these tumors. Molecularly defined subgroups were identified and tested for association with clinical characteristics and overall survival. Results: We detected a putative somatic driver mutation in at least one candidate gene in 95% (401/421) of CCOC tumors including ARID1A (in 49% of tumors), PIK3CA (49%), TERT (20%), and TP53 (16%). Clustering of cancer driver mutations and RNA expression converged upon two distinct subclasses of CCOC. The first was dominated by ARID1A-mutated tumors with enriched expression of canonical CCOC genes and markers of platinum resistance; the second was largely comprised of tumors with TP53 mutations and enriched for the expression of genes involved in extracellular matrix organization and mesenchymal differentiation. Compared with the ARID1A-mutated group, women with TP53-mutated tumors were more likely to have advanced-stage disease, no antecedent history of endometriosis, and poorer survival, driven by their advanced stage at presentation. In women with ARID1A-mutated tumors, there was a trend toward a lower rate of response to first-line platinum-based therapy. Conclusions: Our study suggests that CCOC consists of two distinct molecular subclasses with distinct clinical presentation and outcomes, with potential relevance to both traditional and experimental therapy responsiveness. See related commentary by Lheureux, p. 4838

Abstract Background: There is suggestive evidence that inflammation is related to ovarian cancer survival. However, more research is needed to identify inflammation-related factors that are associated with ovarian cancer survival and to determine their combined effects. Methods: This analysis used pooled data on 8,147 women with invasive epithelial ovarian cancer from the Ovarian Cancer Association Consortium. The prediagnosis inflammation-related exposures of interest included alcohol use; aspirin use; other nonsteroidal anti-inflammatory drug use; body mass index; environmental tobacco smoke exposure; history of pelvic inflammatory disease, polycystic ovarian syndrome, and endometriosis; menopausal hormone therapy use; physical inactivity; smoking status; and talc use. Using Cox proportional hazards models, the relationship between each exposure and survival was assessed in 50% of the data. A weighted inflammation-related risk score (IRRS) was developed, and its association with survival was assessed using Cox proportional hazards models in the remaining 50% of the data. Results: There was a statistically significant trend of increasing risk of death per quartile of the IRRS [HR = 1.09; 95% confidence interval (CI), 1.03–1.14]. Women in the upper quartile of the IRRS had a 31% higher death rate compared with the lowest quartile (95% CI, 1.11–1.54). Conclusions: A higher prediagnosis IRRS was associated with an increased mortality risk after an ovarian cancer diagnosis. Further investigation is warranted to evaluate whether postdiagnosis exposures are also associated with survival. Impact: Given that pre- and postdiagnosis exposures are often correlated and many are modifiable, our study results can ultimately motivate the development of behavioral recommendations to enhance survival among patients with ovarian cancer.

Abstract Background: Ovarian clear cell carcinoma (OCCC) is a rare ovarian cancer histotype that tends to be resistant to standard platinum-based chemotherapeutics. We sought to better understand the role of DNA methylation in clinical and biological subclassification of OCCC. Methods: We interrogated genome-wide methylation using DNA from fresh frozen tumors from 271 cases, applied nonsmooth nonnegative matrix factorization (nsNMF) clustering, and evaluated clinical associations and biological pathways. Results: Two approximately equally sized clusters that associated with several clinical features were identified. Compared with Cluster 2 (N = 137), Cluster 1 cases (N = 134) presented at a more advanced stage, were less likely to be of Asian ancestry, and tended to have poorer outcomes including macroscopic residual disease following primary debulking surgery (P < 0.10). Subset analyses of targeted tumor sequencing and IHC data revealed that Cluster 1 tumors showed TP53 mutation and abnormal p53 expression, and Cluster 2 tumors showed aneuploidy and ARID1A/PIK3CA mutation (P < 0.05). Cluster-defining CpGs included 1,388 CpGs residing within 200 bp of the transcription start sites of 977 genes; 38% of these genes (N = 369 genes) were differentially expressed across cluster in transcriptomic subset analysis (P < 10−4). Differentially expressed genes were enriched for six immune-related pathways, including IFNα and IFNγ responses (P < 10−6). Conclusions: DNA methylation clusters in OCCC correlate with disease features and gene expression patterns among immune pathways. Impact: This work serves as a foundation for integrative analyses that better understand the complex biology of OCCC in an effort to improve potential for development of targeted therapeutics.

Abstract Background: Many loci have been found to be associated with risk of epithelial ovarian cancer (EOC). However, although there is considerable variation in progression-free survival (PFS), no loci have been found to be associated with outcome at genome-wide levels of significance. Methods: We carried out a genome-wide association study (GWAS) of PFS in 2,352 women with EOC who had undergone cytoreductive surgery and standard carboplatin/paclitaxel chemotherapy. Results: We found seven SNPs at 12q24.33 associated with PFS (P < 5 × 10–8), the top SNP being rs10794418 (HR = 1.24; 95% CI, 1.15–1.34; P = 1.47 × 10–8). High expression of a nearby gene, ULK1, is associated with shorter PFS in EOC, and with poor prognosis in other cancers. SNP rs10794418 is also associated with expression of ULK1 in ovarian tumors, with the allele associated with shorter PFS being associated with higher expression, and chromatin interactions were detected between the ULK1 promoter and associated SNPs in serous and endometrioid EOC cell lines. ULK1 knockout ovarian cancer cell lines showed significantly increased sensitivity to carboplatin in vitro. Conclusions: The locus at 12q24.33 represents one of the first genome-wide significant loci for survival for any cancer. ULK1 is a plausible candidate for the target of this association. Impact: This finding provides insight into genetic markers associated with EOC outcome and potential treatment options. See related commentary by Peres and Monteiro, p. 1604

AbstractInterest in analyzing X chromosome single nucleotide polymorphisms (SNPs) is growing and several approaches have been proposed. Prior studies have compared power of different approaches, but bias and interpretation of coefficients have received less attention. We performed simulations to demonstrate the impact of X chromosome model assumptions on effect estimates. We investigated the coefficient biases of SNP and sex effects with commonly used models for X chromosome SNPs, including models with and without assumptions of X chromosome inactivation (XCI), and with and without SNP–sex interaction terms. Sex and SNP coefficient biases were observed when assumptions made about XCI and sex differences in SNP effect in the analysis model were inconsistent with the data‐generating model. However, including a SNP–sex interaction term often eliminated these biases. To illustrate these findings, estimates under different genetic model assumptions are compared and interpreted in a real data example. Models to analyze X chromosome SNPs make assumptions beyond those made in autosomal variant analysis. Assumptions made about X chromosome SNP effects should be stated clearly when reporting and interpreting X chromosome associations. Fitting models with SNP × Sex interaction terms can avoid reliance on assumptions, eliminating coefficient bias even in the absence of sex differences in SNP effect.

Abstract Background: Benign breast disease (BBD) is a risk factor for breast cancer (BC); however, little is known about the genetic alterations present at the time of BBD diagnosis, and how these relate to risk of incident BC. Methods: A subset of a long-term BBD cohort was selected to examine DNA variation across three BBD groups (42 future estrogen receptor-positive (ER+) BC, 36 future estrogen receptor-negative (ER-) BC and 42 controls cancer-free for at least 16 years post-BBD). DNA extracted from archival formalin fixed, paraffin-embedded (FFPE) tissue blocks was tested and analyzed for presence of DNA alterations using a targeted panel of 93 BC-associated genes. To address artifacts frequently observed in FFPE tissues (e.g., C>T changes), we applied three filtering strategies based on alternative allele frequencies and nucleotide substitution context. Gene-level associations were performed using two burden tests, SKAT-O and logistic regression, and adjusted for clinical and technical covariates including epithelial percentage, histologic impression, patient age and year of BBD biopsy. Results: After filtering, the variant frequency of SNPs in our samples was highly consistent with population allele frequencies reported in 1KG/ExAC (0.986, p <1e-16). Ten genes were found to be associated with later cancer status by four of 12 association methods, with nominal p-values <0.05: MED12, MSH2, BRIP1, PMS1, GATA3, MUC16, FAM175A, EXT2, MLH1 and TGFB1. Additionally, analysis revealed 11 and 19 gene-level associations with p-values <0.05 and OR<1 for all cancer cases and ER- cases, respectively. Comparing between the three case groups, BBD ER+ cases were closer to controls in mutation profile, while BBD ER- cases were distinct. Notably, the variant burden was significantly higher in controls than in either ER+ or ER- cases. CD45 expression was associated with mutational burden (p=0.003).Conclusions: Somatic mutations were more frequent in benign breast tissue from women who did not develop cancer, opening questions of clonal diversity or immune-mediated restraint on future cancer development. CD45 expression was positively associated with mutational burden, most strongly in controls. Further studies in both normal and premalignant tissues are needed to better understand the role of somatic gene mutations and their contribution to future cancer development.

Abstract Background: Clinical use of breast cancer risk prediction requires simplified models. We evaluate a simplified version of the validated Rosner–Colditz model and add percent mammographic density (MD) and polygenic risk score (PRS), to assess performance from ages 45–74. We validate using the Mayo Mammography Health Study (MMHS). Methods: We derived the model in the Nurses' Health Study (NHS) based on: MD, 77 SNP PRS and a questionnaire score (QS; lifestyle and reproductive factors). A total of 2,799 invasive breast cancer cases were diagnosed from 1990–2000. MD (using Cumulus software) and PRS were assessed in a nested case–control study. We assess model performance using this case–control dataset and evaluate 10-year absolute breast cancer risk. The prospective MMHS validation dataset includes 21.8% of women age &lt;50, and 434 incident cases identified over 10 years of follow-up. Results: In the NHS, MD has the highest odds ratio (OR) for 10-year risk prediction: ORper SD = 1.48 [95% confidence interval (CI): 1.31–1.68], followed by PRS, ORper SD = 1.37 (95% CI: 1.21–1.55) and QS, ORper SD = 1.25 (95% CI: 1.11–1.41). In MMHS, the AUC adjusted for age + MD + QS 0.650; for age + MD + QS + PRS 0.687, and the NRI was 6% in cases and 16% in controls. Conclusion: A simplified assessment of QS, MD, and PRS performs consistently to discriminate those at high 10-year breast cancer risk. Impact: This simplified model provides accurate estimation of 10-year risk of invasive breast cancer that can be used in a clinical setting to identify women who may benefit from chemopreventive intervention. See related commentary by Tehranifar et al., p. 587

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.