Lambrechts, Diether

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 Background: CA125 is elevated in 80% of ovarian cancer cases and has proven utility in assessing response to therapy and prognosis. Unfortunately, CA125 is elevated in a variety of benign conditions and only in 50% of early stage ovarian cancers, resulting in a sensitivity and specificity unacceptable for population-based screening. Therefore, understanding factors that influence CA125 at presentation could provide important insights for interpreting CA125 values. Methods: Using pre-treatment CA125 and detailed epidemiologic data from 12 studies participating in the Ovarian Cancer Association Consortium (OCAC), we evaluated factors previously reported to influence CA125 levels, including age, race, oral contraceptive use, parity, tubal ligation, endometriosis, body mass index (BMI), personal history of breast cancer, or a family history of breast or ovarian cancer. We used linear regression to estimate the association between each variable and CA125 probit scores, which we used to standardize values between studies. Secondary analyses included adjustment for histologic subtype. We also estimated the associations within each study using log-transformed CA125 as the outcome and estimated summary measures using random effects meta-analyses. Results: Of the 4417 cases included in the analysis, 2918 (66%) were serous, 227 (5%) were mucinous, 484 (11%) were endometrioid, and 258 (6%) were clear cell carcinomas. Median CA125 values varied between studies with a high of 831 U/mL and a low of 271 U/mL. We observed no association between race, oral contraceptive use, tubal ligation, endometriosis, prior breast cancer, and family history of breast cancer and pre-treatment CA125. However, we observed increased pre-treatment CA125 levels with older age (>70 vs. < 50, p=0.03), parity (p=0.01), number of children (p=0.02), obesity (BMI>30 vs. 21-25, p=0.03), and family history of ovarian cancer (p=0.05). Results were similar but attenuated when we calculated summary estimates of the association using meta-analysis of study specific estimates of the association with log-transformed CA125 rather than probit scores. Age and BMI remained predictive of pre-treatment CA125 values after adjustment for histologic subtype. Our data are limited by between site variability in CA125 assays; therefore, validation of these findings is needed with adjustment for type of assay used or in a study population in which all cases had pre-treatment CA125 values measured with the same assay. Conclusions: Despite these limitations, results from this analysis of over 4000 ovarian cancer cases suggest that age and body size may influence pre-treatment CA125 values. Citation Format: KL Terry, A Babic, BY Karlan, MT Goodman, D Lambrechts, F Heitz, K Matsuo, I McNeish, T Pejovic, S Kruger Kjaer, PM Webb, E Hogdall, EL Goode, DW Cramer, for the Ovarian Cancer Association Consortium. Epidemiologic predictors of pre-treatment CA125 in women with ovarian cancer [abstract]. In: Proceedings of the 10th Biennial Ovarian Cancer Research Symposium; Sep 8-9, 2014; Seattle, WA. Philadelphia (PA): AACR; Clin Cancer Res 2015;21(16 Suppl):Abstract nr AS13.

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 < 0.05 and FDR < 0.05). These results were replicated (P < 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.

Studies assessing the effect of bevacizumab (BEV) on breast cancer (BC) outcome have shown different effects on progression‐free and overall survival, suggesting that a subgroup of patients may benefit from this treatment. Unfortunately, no biomarkers exist to identify these patients. Here, we investigate whether single nucleotide polymorphisms (SNPs) in VEGF pathway genes correlate with pathological complete response (pCR) in the neoadjuvant GeparQuinto trial. HER2‐negative patients were randomized into treatment arms receiving either BEV combined with standard chemotherapy or chemotherapy alone. In a pre‐planned biomarker study, DNA was collected from 729 and 724 patients, respectively from both treatment arms, and genotyped for 125 SNPs. Logistic regression assessed interaction between individual SNPs and both treatment arms to predict pCR. Five SNPs may be associated with a better response to BEV, but none of them remained significant after correction for multiple testing. The two SNPs most strongly associated, rs833058 and rs699947, were located upstream of the VEGF‐A promoter. Odds ratios for the homozygous common, heterozygous and homozygous rare rs833058 genotypes were 2.36 (95% CI, 1.49–3.75), 1.20 (95% CI, 0.88–1.64) and 0.61 (95% CI, 0.34–1.12). Notably, some SNPs in VEGF‐A exhibited a more pronounced effect in the triple‐negative subgroup. Several SNPs in VEGF‐A may be associated with improved pCR when receiving BEV in the neoadjuvant setting. Although none of the observed effects survived correction for multiple testing, our observations are consistent with previous studies on BEV efficacy in BC. Further research is warranted to clarify the predictive value of these markers.

Abstract Purpose: To analyze the effect of germline mutations in BRCA1 and BRCA2 on mortality in patients with ovarian cancer up to 10 years after diagnosis. Experimental Design: We used unpublished survival time data for 2,242 patients from two case–control studies and extended survival time data for 4,314 patients from previously reported studies. All participants had been screened for deleterious germline mutations in BRCA1 and BRCA2. Survival time was analyzed for the combined data using Cox proportional hazard models with BRCA1 and BRCA2 as time-varying covariates. Competing risks were analyzed using Fine and Gray model. Results: The combined 10-year overall survival rate was 30% [95% confidence interval (CI), 28%–31%] for non-carriers, 25% (95% CI, 22%–28%) for BRCA1 carriers, and 35% (95% CI, 30%–41%) for BRCA2 carriers. The HR for BRCA1 was 0.53 at time zero and increased over time becoming greater than one at 4.8 years. For BRCA2, the HR was 0.42 at time zero and increased over time (predicted to become greater than 1 at 10.5 years). The results were similar when restricted to 3,202 patients with high-grade serous tumors and to ovarian cancer–specific mortality. Conclusions: BRCA1/2 mutations are associated with better short-term survival, but this advantage decreases over time and in BRCA1 carriers is eventually reversed. This may have important implications for therapy of both primary and relapsed disease and for analysis of long-term survival in clinical trials of new agents, particularly those that are effective in BRCA1/2 mutation carriers. Clin Cancer Res; 21(3); 652–7. ©2014 AACR.

A large genotyping project within the Breast Cancer Association Consortium (BCAC) recently identified 41 associations between single nucleotide polymorphisms (SNPs) and overall breast cancer (BC) risk. We investigated whether the effects of these 41 SNPs, as well as six SNPs associated with estrogen receptor (ER) negative BC risk are modified by 13 environmental risk factors for BC. Data from 22 studies participating in BCAC were pooled, comprising up to 26,633 cases and 30,119 controls. Interactions between SNPs and environmental factors were evaluated using an empirical Bayes‐type shrinkage estimator. Six SNPs showed interactions with associated p ‐values ( p int ) <1.1 × 10 −3 . None of the observed interactions was significant after accounting for multiple testing. The Bayesian False Discovery Probability was used to rank the findings, which indicated three interactions as being noteworthy at 1% prior probability of interaction. SNP rs6828523 was associated with increased ER‐negative BC risk in women ≥170 cm (OR = 1.22, p = 0.017), but inversely associated with ER‐negative BC risk in women <160 cm (OR = 0.83, p = 0.039, p int = 1.9 × 10 −4 ). The inverse association between rs4808801 and overall BC risk was stronger for women who had had four or more pregnancies (OR = 0.85, p = 2.0 × 10 −4 ), and absent in women who had had just one (OR = 0.96, p = 0.19, p int = 6.1 × 10 −4 ). SNP rs11242675 was inversely associated with overall BC risk in never/former smokers (OR = 0.93, p = 2.8 × 10 −5 ), but no association was observed in current smokers (OR = 1.07, p = 0.14, p int = 3.4 × 10 −4 ). In conclusion, recently identified BC susceptibility loci are not strongly modified by established risk factors and the observed potential interactions require confirmation in independent studies.