Brinton, Louise A.

Previous genome‐wide association studies among women of European ancestry identified two independent breast cancer susceptibility loci represented by single nucleotide polymorphisms (SNPs) rs13281615 and rs11780156 at 8q24. A fine‐mapping study across 2.06 Mb (chr8:127,561,724–129,624,067, hg19) in 55,540 breast cancer cases and 51,168 controls within the Breast Cancer Association Consortium was conducted. Three additional independent association signals in women of European ancestry, represented by rs35961416 (OR = 0.95, 95% CI = 0.93–0.97, conditional p = 5.8 × 10−6), rs7815245 (OR = 0.94, 95% CI = 0.91–0.96, conditional p = 1.1 × 10−6) and rs2033101 (OR = 1.05, 95% CI = 1.02–1.07, conditional p = 1.1 × 10−4) were found. Integrative analysis using functional genomic data from the Roadmap Epigenomics, the Encyclopedia of DNA Elements project, the Cancer Genome Atlas and other public resources implied that SNPs rs7815245 in Signal 3, and rs1121948 in Signal 5 (in linkage disequilibrium with rs11780156, r2 = 0.77), were putatively functional variants for two of the five independent association signals. The results highlighted multiple 8q24 variants associated with breast cancer susceptibility in women of European ancestry.

Abstract Background: Many epithelial ovarian cancer (EOC) risk factors relate to hormone exposure and elevated estrogen levels are associated with obesity in postmenopausal women. Therefore, we hypothesized that gene–environment interactions related to hormone-related risk factors could differ between obese and non-obese women. Methods: We considered interactions between 11,441 SNPs within 80 candidate genes related to hormone biosynthesis and metabolism and insulin-like growth factors with six hormone-related factors (oral contraceptive use, parity, endometriosis, tubal ligation, hormone replacement therapy, and estrogen use) and assessed whether these interactions differed between obese and non-obese women. Interactions were assessed using logistic regression models and data from 14 case–control studies (6,247 cases; 10,379 controls). Histotype-specific analyses were also completed. Results: SNPs in the following candidate genes showed notable interaction: IGF1R (rs41497346, estrogen plus progesterone hormone therapy, histology = all, P = 4.9 × 10−6) and ESR1 (rs12661437, endometriosis, histology = all, P = 1.5 × 10−5). The most notable obesity–gene–hormone risk factor interaction was within INSR (rs113759408, parity, histology = endometrioid, P = 8.8 × 10−6). Conclusions: We have demonstrated the feasibility of assessing multifactor interactions in large genetic epidemiology studies. Follow-up studies are necessary to assess the robustness of our findings for ESR1, CYP11A1, IGF1R, CYP11B1, INSR, and IGFBP2. Future work is needed to develop powerful statistical methods able to detect these complex interactions. Impact: Assessment of multifactor interaction is feasible, and, here, suggests that the relationship between genetic variants within candidate genes and hormone-related risk factors may vary EOC susceptibility. Cancer Epidemiol Biomarkers Prev; 25(5); 780–90. ©2016 AACR.

Abstract Background: While numerous susceptibility loci for epithelial ovarian cancer (EOC) have been identified, few associations have been reported with overall survival. In the absence of common prognostic genetic markers, we hypothesize that rare coding variants may be associated with overall EOC survival and assessed their contribution in two exome-based genotyping projects of the Ovarian Cancer Association Consortium (OCAC). Methods: The primary patient set (Set 1) included 14 independent EOC studies (4,293 patients) and 227,892 variants, and a secondary patient set (Set 2) included six additional EOC studies (1,744 patients) and 114,620 variants. Because power to detect rare variants individually is reduced, gene-level tests were conducted. Sets were analyzed separately at individual variants and by gene, and then combined with meta-analyses (73,203 variants and 13,163 genes overlapped). Results: No individual variant reached genome-wide statistical significance. A SNP previously implicated to be associated with EOC risk and, to a lesser extent, survival, rs8170, showed the strongest evidence of association with survival and similar effect size estimates across sets (Pmeta = 1.1E−6, HRSet1 = 1.17, HRSet2 = 1.14). Rare variants in ATG2B, an autophagy gene important for apoptosis, were significantly associated with survival after multiple testing correction (Pmeta = 1.1E−6; Pcorrected = 0.01). Conclusions: Common variant rs8170 and rare variants in ATG2B may be associated with EOC overall survival, although further study is needed. Impact: This study represents the first exome-wide association study of EOC survival to include rare variant analyses, and suggests that complementary single variant and gene-level analyses in large studies are needed to identify rare variants that warrant follow-up study. Cancer Epidemiol Biomarkers Prev; 25(3); 446–54. ©2016 AACR.

ABSTRACTEpithelial‐mesenchymal transition (EMT) is a process whereby epithelial cells assume mesenchymal characteristics to facilitate cancer metastasis. However, EMT also contributes to the initiation and development of primary tumors. Prior studies that explored the hypothesis that EMT gene variants contribute to epithelial ovarian carcinoma (EOC) risk have been based on small sample sizes and none have sought replication in an independent population. We screened 15,816 single‐nucleotide polymorphisms (SNPs) in 296 genes in a discovery phase using data from a genome‐wide association study of EOC among women of European ancestry (1,947 cases and 2,009 controls) and identified 793 variants in 278 EMT‐related genes that were nominally (P < 0.05) associated with invasive EOC. These SNPs were then genotyped in a larger study of 14,525 invasive‐cancer patients and 23,447 controls. A P‐value <0.05 and a false discovery rate (FDR) <0.2 were considered statistically significant. In the larger dataset, GPC6/GPC5 rs17702471 was associated with the endometrioid subtype among Caucasians (odds ratio (OR) = 1.16, 95% CI = 1.07–1.25, P = 0.0003, FDR = 0.19), whereas F8 rs7053448 (OR = 1.69, 95% CI = 1.27–2.24, P = 0.0003, FDR = 0.12), F8 rs7058826 (OR = 1.69, 95% CI = 1.27–2.24, P = 0.0003, FDR = 0.12), and CAPN13 rs1983383 (OR = 0.79, 95% CI = 0.69–0.90, P = 0.0005, FDR = 0.12) were associated with combined invasive EOC among Asians. In silico functional analyses revealed that GPC6/GPC5 rs17702471 coincided with DNA regulatory elements. These results suggest that EMT gene variants do not appear to play a significant role in the susceptibility to EOC.

Abstract Background: A recent association study identified a common variant (rs9790517) at 4q24 to be associated with breast cancer risk. Independent association signals and potential functional variants in this locus have not been explored. Methods: We conducted a fine-mapping analysis in 55,540 breast cancer cases and 51,168 controls from the Breast Cancer Association Consortium. Results: Conditional analyses identified two independent association signals among women of European ancestry, represented by rs9790517 [conditional P = 2.51 × 10−4; OR, 1.04; 95% confidence interval (CI), 1.02–1.07] and rs77928427 (P = 1.86 × 10−4; OR, 1.04; 95% CI, 1.02–1.07). Functional annotation using data from the Encyclopedia of DNA Elements (ENCODE) project revealed two putative functional variants, rs62331150 and rs73838678 in linkage disequilibrium (LD) with rs9790517 (r2 ≥ 0.90) residing in the active promoter or enhancer, respectively, of the nearest gene, TET2. Both variants are located in DNase I hypersensitivity and transcription factor–binding sites. Using data from both The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC), we showed that rs62331150 was associated with level of expression of TET2 in breast normal and tumor tissue. Conclusion: Our study identified two independent association signals at 4q24 in relation to breast cancer risk and suggested that observed association in this locus may be mediated through the regulation of TET2. Impact: Fine-mapping study with large sample size warranted for identification of independent loci for breast cancer risk. Cancer Epidemiol Biomarkers Prev; 24(11); 1680–91. ©2015 AACR.

Abstract Purpose: Breast biopsies with benign findings (benign breast disease [BBD]) outnumber biopsies revealing in-situ or invasive cancer 4:1, providing an opportunity to identify histopathological features that predict breast cancer risk. Greater degrees of involution of terminal duct lobular units (TDLUs), the structures within the breast that produce milk and the primary source of breast cancer precursors, have been inversely associated with breast cancer risk among women with BBD. We recently developed measures to quantitate levels of TDLU involution, which demonstrate high inter/intra pathologist reproducibility. Here we assessed whether TDLU counts/100mm2 and median TDLU span (microns), two measures inversely related to degree of TDLU involution, are associated with subsequent breast cancer risk among women with BBD. Methods: From the Mayo BBD cohort (n = 9,376), we evaluated benign biopsies from 99 women who later developed breast cancer (cases) and 145 age-matched controls who did not develop breast cancer. Digitized images of biopsy sections were reviewed to enumerate TDLUs/mm2 and measure median TDLU span (microns) for up to ten normal TDLUs. Breast cancer risk factors were available from questionnaires or medical records. To assess associations with breast cancer risk, subjects were categorized into quartile levels of TDLU counts and median TDLU spans, based on data from controls. Odds ratios (ORs) and 95% confidence intervals (CI) were estimated using multivariable conditional logistic regression models, adjusted for histologic type of BBD, family history of breast cancer and menopausal hormone use. Results: Compared to controls, women who developed breast cancer had higher median number of TDLU counts/100mm2 (28 vs 20, p = 0.03) and larger TDLU spans (300 vs 267 microns, p = 0.14). In multivariable models higher TDLU counts (quartile4 vs. quartile1, OR = 2.44, 95%CI = 0.96-6.19, ptrend = 0.02) and larger TDLU span measures (quartile4 vs. quartile1, OR = 2.83, 95%CI = 1.13-7.06, ptrend = 0.03) were associated with subsequent diagnosis of breast cancer. Combinatorial metrics of TDLU counts with median TDLU span measures identified women at higher risk; specifically, women above the median for both TDLU span and TDLU counts had an OR = 3.75 (95%CI = 1.40-10.00, ptrend = 0.008), compared with women below the median for TDLU span and TDLU counts. Conclusion: These data show that lack of TDLU involution, as measured by increased persistence of TDLU counts, larger median TDLU spans and cross-classification using these measures was associated with increased breast cancer risk among women with BBD, extending prior work in this cohort based on other TDLU involution metrics. Future studies to identify determinants of TDLU involution, its association with breast cancer risk and its potential as an intermediate endpoint in prevention studies warrant consideration. Citation Format: Jonine D. Figueroa, Ruth Pfeiffer, Maya Palakal, Amy C. Degnim, Derek Radisky, Lynn C. Hartmann, Marlene Frost, Melody L. Stallings Mann, Louise A. Brinton, Daphne Papathomas, Daniel Visscher, Mark E. Sherman. Standardized measures of lobular involution and subsequent breast cancer risk among women with benign breast disease. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4682. doi:10.1158/1538-7445.AM2015-4682

Abstract Introduction: Increased levels of circulating estrogens are associated with elevated risks of postmenopausal breast cancer. Previously, we showed that levels of sex steroid hormones in breast adipose tissue are positively correlated with circulating hormones. Further evaluation of hormone levels in breast tissues may help to explain the relationship between estrogens and breast cancer risk. Specifically, data suggest that focal breast lesions consisting of macrophages encircling dying adipocytes, referred to as “crown-like structures (CLS)”, are associated with increased aromatase activity, which could increase conversion of androgens to estrogens locally in the breast. Therefore, we assessed whether the presence of CLS affects the correlations between blood and breast adipose tissue hormone levels. Methods: Participants included 84 postmenopausal women with invasive breast cancer in the Polish Breast Cancer Study, a population based case-control study conducted from 2000 to 2003. We previously reported sex steroid hormone levels in cancer-free breast adipose tissue removed as part of breast excisions performed for cancer. Here, using immunohistochemical stains, we evaluated the number of CD68 positive macrophages and the presence of full CLS (i.e., complete adipocyte encirclement) in benign breast tissue sections prepared from surgical pathology specimens containing breast cancer. The relationship between the number of CD68 positive macrophages per unit area of fat and hormone levels was examined using Spearman rank correlations (rs). We assessed potential effect modification of the association between blood and tissue hormone levels by presence or absence of CLS by including an interaction term in linear regression models. Results: The number of CD68 positive macrophages per unit area of fat was not significantly correlated with levels of estrone, estradiol, androstenedione or testosterone measured in blood (rs = 0.001, -0.04, -0.07, and -0.09, respectively; p&gt;0.05) or breast adipose tissue (rs = -0.15, 0.06, -0.06, and -0.08, respectively; p&gt;0.05). Correlations between serum and adipose tissue hormones did not statistically significantly differ between women with CLS (n = 30) compared to those without CLS (n = 54). CLS were found more frequently in breast adipose tissue from overweight (i.e., women with a body mass index (BMI) between 25 and 29.9) (43%) and obese women (BMI &gt;30) (30%) versus normal weight women, with (BMI between 18.5 and 24.9), (27%); (p = 0.04, χ2 test). The number of CD68 positive macrophages per unit area of fat was not significantly related to BMI (p = 0.31). Conclusion: Among postmenopausal women with invasive breast cancer, the presence of CLS was more common among overweight and obese women, but CLS did not influence the correlation between circulating and breast adipose tissue sex steroid hormone levels. Citation Format: Maeve Mullooly, Hannah P. Yang, Roni T. Falk, Sarah Nyante, Renata Cora, Ruth M. Pfeiffer, Derek C. Radisky, Daniel W. Visscher, Lynn C. Hartmann, Amy C. Degnim, Frank Z. Stanczyk, Jonine F. Figueroa, Montserrat Garcia-Closas, Jolanta Lissowska, Melissa T. Troester, Louise A. Brinton, Mark E. Sherman, Gretchen L. Gierach. Investigation of the relationship between crown-like structures and adipose tissue hormone levels among postmenopausal women with breast cancer. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2767. doi:10.1158/1538-7445.AM2015-2767

Abstract Background: Genome-wide association studies (GWAS) have so far reported 12 loci associated with serous epithelial ovarian cancer (EOC) risk. We hypothesized that some of these loci function through nearby transcription factor (TF) genes and that putative target genes of these TFs as identified by coexpression may also be enriched for additional EOC risk associations. Methods: We selected TF genes within 1 Mb of the top signal at the 12 genome-wide significant risk loci. Mutual information, a form of correlation, was used to build networks of genes strongly coexpressed with each selected TF gene in the unified microarray dataset of 489 serous EOC tumors from The Cancer Genome Atlas. Genes represented in this dataset were subsequently ranked using a gene-level test based on results for germline SNPs from a serous EOC GWAS meta-analysis (2,196 cases/4,396 controls). Results: Gene set enrichment analysis identified six networks centered on TF genes (HOXB2, HOXB5, HOXB6, HOXB7 at 17q21.32 and HOXD1, HOXD3 at 2q31) that were significantly enriched for genes from the risk-associated end of the ranked list (P &lt; 0.05 and FDR &lt; 0.05). These results were replicated (P &lt; 0.05) using an independent association study (7,035 cases/21,693 controls). Genes underlying enrichment in the six networks were pooled into a combined network. Conclusion: We identified a HOX-centric network associated with serous EOC risk containing several genes with known or emerging roles in serous EOC development. Impact: Network analysis integrating large, context-specific datasets has the potential to offer mechanistic insights into cancer susceptibility and prioritize genes for experimental characterization. Cancer Epidemiol Biomarkers Prev; 24(10); 1574–84. ©2015 AACR.

Abstract Purpose: Chemotherapy resistance remains a major challenge in the treatment of ovarian cancer. We hypothesize that germline polymorphisms might be associated with clinical outcome. Experimental Design: We analyzed approximately 2.8 million genotyped and imputed SNPs from the iCOGS experiment for progression-free survival (PFS) and overall survival (OS) in 2,901 European epithelial ovarian cancer (EOC) patients who underwent first-line treatment of cytoreductive surgery and chemotherapy regardless of regimen, and in a subset of 1,098 patients treated with ≥4 cycles of paclitaxel and carboplatin at standard doses. We evaluated the top SNPs in 4,434 EOC patients, including patients from The Cancer Genome Atlas. In addition, we conducted pathway analysis of all intragenic SNPs and tested their association with PFS and OS using gene set enrichment analysis. Results: Five SNPs were significantly associated (P ≤ 1.0 × 10−5) with poorer outcomes in at least one of the four analyses, three of which, rs4910232 (11p15.3), rs2549714 (16q23), and rs6674079 (1q22), were located in long noncoding RNAs (lncRNAs) RP11-179A10.1, RP11-314O13.1, and RP11-284F21.8, respectively (P ≤ 7.1 × 10−6). ENCODE ChIP-seq data at 1q22 for normal ovary show evidence of histone modification around RP11-284F21.8, and rs6674079 is perfectly correlated with another SNP within the super-enhancer MEF2D, expression levels of which were reportedly associated with prognosis in another solid tumor. YAP1- and WWTR1 (TAZ)-stimulated gene expression and high-density lipoprotein (HDL)-mediated lipid transport pathways were associated with PFS and OS, respectively, in the cohort who had standard chemotherapy (pGSEA ≤6 × 10−3). Conclusions: We have identified SNPs in three lncRNAs that might be important targets for novel EOC therapies. Clin Cancer Res; 21(23); 5264–76. ©2015 AACR.