Nakamura, Yusuke

Abstract Purpose: Germline pathogenic variants in CHEK2 confer moderately elevated breast cancer risk (odds ratio, OR ∼ 2.5), qualifying carriers for enhanced breast cancer screening. Besides pathogenic variants, dozens of missense CHEK2 variants of uncertain significance (VUS) have been identified, hampering the clinical utility of germline genetic testing (GGT). Experimental Design: We collected 460 CHEK2 missense VUS identified by the ENIGMA consortium in 15 countries. Their functional characterization was performed using CHEK2-complementation assays quantifying KAP1 phosphorylation and CHK2 autophosphorylation in human RPE1–CHEK2-knockout cells. Concordant results in both functional assays were used to categorize CHEK2 VUS from 12 ENIGMA case–control datasets, including 73,048 female patients with breast cancer and 88,658 ethnicity-matched controls. Results: A total of 430/460 VUS were successfully analyzed, of which 340 (79.1%) were concordant in both functional assays and categorized as functionally impaired (N = 102), functionally intermediate (N = 12), or functionally wild-type (WT)–like (N = 226). We then examined their association with breast cancer risk in the case–control analysis. The OR and 95% CI (confidence intervals) for carriers of functionally impaired, intermediate, and WT-like variants were 2.83 (95% CI, 2.35–3.41), 1.57 (95% CI, 1.41–1.75), and 1.19 (95% CI, 1.08–1.31), respectively. The meta-analysis of population-specific datasets showed similar results. Conclusions: We determined the functional consequences for the majority of CHEK2 missense VUS found in patients with breast cancer (3,660/4,436; 82.5%). Carriers of functionally impaired missense variants accounted for 0.5% of patients with breast cancer and were associated with a moderate risk similar to that of truncating CHEK2 variants. In contrast, 2.2% of all patients with breast cancer carried functionally wild-type/intermediate missense variants with no clinically relevant breast cancer risk in heterozygous carriers.

Abstract Purpose: Tamoxifen and raloxifene, are primary prevention strategies for women at high risk of breast cancer (Visvanathan, 2013). Recent advances in genetic studies of breast cancer risk have identified common susceptibility loci that explain 14% of familial risk for breast cancer in the general population (Michailidou, 2013). However, it is not known if these loci are risk factors for breast cancer among high-risk women treated with SERMs for breast cancer prevention. We hypothesized that the large risk reduction associated with SERMs, coupled with the fact that several breast cancer loci correlate with family history, may limit the contribution of these common genetic loci to breast cancer in this high risk population. We present the first report to evaluate 75 established breast cancer susceptibility loci, in the context of a polygenic risk score (PRS), as a risk factor for breast cancer among high risk women taking raloxifene and tamoxifen for breast cancer prevention. Methods: We conducted a matched case-control study of 594 cases (i.e., participants who developed breast cancer while on SERM therapy) and 1,171 matched controls selected from the 33,000 participants enrolled in the NSABP P-1 and P-2 breast cancer prevention trials. Genotypes of 75 single nucleotide polymorphisms (SNPs) were available from a genome-wide association study conducted at the RIKEN Center for Genomic Medicine. We formed a quantitative PRS from reported per-SNP odds ratios (OR) for the 75 susceptibility loci. Conditional logistic regression was used to examine the PRS as a risk factor for breast cancer and to assess whether the PRS and breast cancer association differed by treatment type, family history, or other clinical characteristics. Analyses also examined associations of PRS with invasive vs. in situ cancer and ER-positive vs. ER-negative cancer. Results: The PRS ranged from 3.98 to 7.74, and a one unit change in PRS was associated with a 42% increase in breast cancer (OR=1.42; 95% CI: 1.18-1.70; P = 0.0002). There was evidence of a stronger association of PRS with breast cancer among women with no first-degree family history (OR=1.62, 95% CI: 1.18-2.21) compared to those with a positive family history (OR=1.32, 95% CI: 1.06-1.66) (Pintx<0.05). The PRS also appeared a stronger risk factor for ER-positive (OR=1.59, 95% CI: 1.25-2.02, P < 0.0002) vs. ER-negative (OR=1.05, 95% CI: 0.68-1.62, P=0.84) breast cancer, although differences did not reach statistical significance (Pintx=0.10). PRS and breast cancer associations were similar across tamoxifen and raloxifene treatments, age at trial entry, 5-year predicted Gail model risk, hysterectomy status, BMI, presence of atypical hyperplasia and invasive vs. in situ cancer. Conclusion: A polygenic risk score composed of 75 loci was a risk factor for ER-positive breast cancer, especially in the absence of a first-degree family history of breast cancer. Further, the PRS associations with breast cancer were similar for women taking tamoxifen or raloxifene for prevention. These data suggest that common genetic variation adds information on risk of ER-positive breast cancer in a high-risk population receiving SERMs. Citation Format: Celine M Vachon, Daniel J Schaid, James N Ingle, Matthew P Goetz, Donald L Wickerham, Michiaki Kubo, Erin E Carlson, Soonmyung Paik, Norman Wolmark, Yusuke Nakamura, Liewei Wang, Richard M Weinshilboum, Fergus J Couch. The contribution of common genetic variation to breast cancer risk among women receiving tamoxifen or raloxifene within the National Surgical Adjuvant Breast and Bowel Project (NSABP) P-1 and P-2 trials [abstract]. In: Proceedings of the Thirty-Seventh Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2014 Dec 9-13; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2015;75(9 Suppl):Abstract nr P6-10-03.

Abstract The selective estrogen receptor modulators (SERM) tamoxifen and raloxifene can reduce the occurrence of breast cancer in high-risk women by 50%, but this U.S. Food and Drug Administration-approved prevention therapy is not often used. We attempted to identify genetic factors that contribute to variation in SERM breast cancer prevention, using DNA from the NSABP P-1 and P-2 breast cancer prevention trials. An initial discovery genome-wide association study identified common single-nucleotide polymorphisms (SNP) in or near the ZNF423 and CTSO genes that were associated with breast cancer risk during SERM therapy. We then showed that both ZNF423 and CTSO participated in the estrogen-dependent induction of BRCA1 expression, in both cases with SNP-dependent variation in induction. ZNF423 appeared to be an estrogen-inducible BRCA1 transcription factor. The OR for differences in breast cancer risk during SERM therapy for subjects homozygous for both protective or both risk alleles for ZNF423 and CTSO was 5.71. Significance: This study identified novel, functionally polymorphic genes involved in the estrogen-dependent regulation of BRCA1 expression, as well as a novel mechanism for genetic variation in SERM therapeutic effect. These observations, and definition of their underlying mechanisms, represent steps toward pharmacogenomically individualized SERM breast cancer prevention. Cancer Discov; 3(7); 812–25. ©2013 AACR. See related commentary by Machiela and Chanock, p. 728 This article is highlighted in the In This Issue feature, p. 705

Pharmacogenomics research provides insights into a personalized approach to prescribing and dosing medications.

Abstract Introduction: The NSABP P-1 and P-2 breast cancer prevention trials have demonstrated that the SERMs tamoxifen and raloxifene can reduce the risk of breast cancer in high risk women. Methods: We performed a nested matched case-control GWAS utilizing DNA from participants enrolled in P-1 and P-2. 592 participants developed invasive breast cancer or ductal carcinoma in situ, and were matched with 1171 controls who did not. Genome-wide genotyping with the Illumina 610-Quad chip and CYP2D6 genotyping for pertinent alleles were performed. Functional genomic studies involved the use of cultured cells, with siRNA knockdown, followed by qRT-PCR and Western blot analyses. Results: Eleven SNPs with p-values < 3E-05 were identified. Twenty one additional SNPs with p < 4E-05 were identified by imputation around the genotyped SNPs on chromosomes (Chrs) 4, 8, 9, 13, 16 and 22, which were then genotyped with the Invader platform. Initial functional genomic studies focused on 6 SNPs on Chr16 (p-values 1.81-9.55E-06), all of which were in ZNF423, a gene encoding a putative zinc-finger protein. The Chr16 minor variant SNP sequences had less risk than the common variant (OR=0.7). Our functional studies showed that incubation of Hs578T cells stably transfected with estrogen receptor (ER) α with 0.1 nM estradiol (E2) induced expression of ZNF423, BRCA1/2 for WT but not variant SNP sequences. There was also a striking difference between E2-induced ZNF423, BRCA1/2 expression for variant and WT SNPs in lymphoblastoid cell lines (LCLs) stably transfected with ERα, with only WT showing induction. ChIP assay showed that ZNF423 could bind to the 5’-flanking region of BRCA1, and reporter gene assays showed that ZNF423 could regulate BRCA1/2 transcription. In LCLs stably transfected with ERα and genotyped for ZNF423 SNPs, estrogen-dependent expression of all 3 genes occurred only in the presence of the WT, but not variant ZNF423 SNP sequences. Blockade of ERα with 4-OH tamoxifen prevented induction of expression in cells with WT sequences, but resulted in gene dose-dependent increases in BRCA1/2 expression in cells with variant SNP sequences-suggested that the minor variant has a “protective” effect during clinical SERM therapy and verified by ChIP assays. Conclusions: The SNPs in ZNF423 identified during GWAS were associated with differential E2-dependent induction of ZNF423, BRCA1 and BRCA2-implying that ZNF423 is “upstream” for the estrogen-dependent induction of BRCA1/2 expression. However, LCLs with variant SNP sequences showed greatly enhanced expression of ZNF423, BRCA1/2 during 4-OH tamoxifen exposure, while cells with WT sequences did not. These results may reveal a novel mechanism for SERM-dependent regulation of BRCA1/2 expression and may have implications for patient selection for SERM therapy. Funded by U10-CA37377, U10-CA69974 and U14-GM61388. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 4727. doi:10.1158/1538-7445.AM2011-4727

Abstract Introduction: Androstenedione is a precursor for estrogen biosynthesis catalyzed by aromatase. Plasma androstenedione concentrations in postmenopausal women with early stage breast cancer display large interindividual variation (Ingle et al., Cancer Res. 70:3278–86, 2010). We performed a GWAS to identify single nucleotide polymorphisms (SNPs) associated with plasma levels of androstenedione in 776 postmenopausal women with resected early stage breast cancer. Methods: Plasma androstenedione concentrations were measured by GC-MS/MS. GWAS was performed using Illumina Human610-Quad BeadChips. The GWAS results were controlled for population stratification (PS) by Eigenstrat analyses. Plasma androstenedione concentrations were associated with genotype by using a quasi-likelihood model with adjustment for PS, BMI, age and other clinical variables associated with androstenedione concentrations (p<0.01). Imputation using HapMap2 data was performed within 200 kb on either side of regions containing SNPs with p<E-05. Results: 563,945 SNPs were analyzed, followed by imputation. The genotyped SNP with the lowest p value, rs4736317 (p=4.61E-06), was on chromosome 8 near the CYP11B1 and CYP11B2 genes. This SNP had a multiplicative effect of 1.128 (95% Cl: 1.07, 1.19) and a MAF of 0.40. After imputation, 20 SNPs with apparent p<E-06 were observed in this same region, with the lowest p value for rs9773022 (p=8.1E-08). We then performed functional genomic studies, and observed that two of the SNPs with low p values (rs4736312, p=1.49E-05 and rs4736350, p=6.34E-06) created apparent Oct-1 transcription factor binding sites (TFBS). Differential nuclear protein binding of Oct-1 to probes for these SNPs was verified with EMS and supershift assays using nuclear extract from A549 and NCI-H295R cells. Treatment of these two cell lines with 10 nM androstenedione led to the disruption of Oct-1 binding with the A (wild type) allele on rs4736312. rs4736312 is located in the 3'-UTR of CYP11B1 15 bp distant from response elements for the progesterone receptor and the glucocorticoid receptor. We also demonstrated that CYP11B1 and CYP11B2 were induced by estradiol in U2OS cell lines that express ER receptor and by androstenedione in H295R cells that express the androgen receptor. Transfection of a reporter gene construct with the rs4736312 SNP sequences into COS-1 and H295R cells suppressed reporter gene activity for the variant SNP sequence as compared to wild type, an observation that may correlate with the higher plasma levels of androstenedione in women (549× 328 vs 410×208 mean×SD, M/ml) with the variant allele. Conclusion: We have identified a series of SNPs in a region of chromosome 8 in and around CYP11B1 and CYP11B2 that are associated with elevated plasma androstenedione concentrations in postmenopausal women with resected early stage breast cancer. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference: Molecular Targets and Cancer Therapeutics; 2011 Nov 12-16; San Francisco, CA. Philadelphia (PA): AACR; Mol Cancer Ther 2011;10(11 Suppl):Abstract nr C114.

Abstract Background and Aims: Cancer and Leukemia Group B 80303 was a randomized, phase III study in patients with advanced pancreatic cancer treated with gemcitabine plus either bevacizumab or placebo. We prospectively collected germline DNA and conducted a genome-wide association study (GWAS) using overall survival (OS) as the endpoint. Experimental Design: DNA from 351 patients was genotyped for more than 550,000 single-nucleotide polymorphisms (SNP). Associations between OS and SNPs were investigated using the log-linear 2-way multiplicative Cox proportional hazards model. The subset of 294 genetically European patients was used for the primary analysis. Results: A nonsynonymous SNP in interleukin (IL)17F (rs763780, H161R) and an intronic SNP in strong linkage disequilibrium (rs7771466) were associated with OS using genome-wide criteria (P ≤ 10−7). Median OS was significantly shorter (P = 2.61 × 10−8) for the rs763780 heterozygotes [3.1 months; 95% confidence interval (CI), 2.3–4.3] than for the patients without this variant (6.8 months; 95% CI, 5.8–7.3). After adjustment by stratification factors, the P value for the association was 9.51 × 10−7. Conclusions: The variant 161R form of IL-17F is a natural antagonist of the antiangiogenic effects of wild-type 161H IL-17F, and angiogenesis may play an important role in the metastatic spread of pancreatic cancer. In this preliminary study, we hypothesize that the angiogenetic potential of pancreatic cancers in patients with variant IL-17F is higher than that of tumors in patients with wild-type IL-17F, conferring worse prognosis. This exploratory GWAS may provide the foundation for testing the biology and clinical effects of novel genes and their heritable variants through mechanistic and confirmatory studies in pancreatic cancer. Clin Cancer Res; 18(2); 577–84. ©2011 AACR.

Abstract Background: Controversy exists regarding the association between CYP2D6 enzyme activity and tamoxifen effectiveness in the adjuvant treatment of invasive breast cancer; however, this association in the primary prevention of breast cancer is unknown. Methods: We conducted a nested case–control study in the context of the NSABP P1 and P2 prevention trials to determine the impact of CYP2D6 genotype, CYP2D6 inhibitor use, and metabolizer status (CYP2D6 genotype combined with CYP2D6 inhibitor use), on breast cancer events. Women who developed breast cancer (both noninvasive and invasive) while on 5 years of selective estrogen receptor modulators therapy (cases) were matched to controls free of breast cancer. Comprehensive CYP2D6 genotyping was conducted for alleles associated with absent (*3, *4, *5, and *6), reduced (*10, *17, and *41), and increased (*1XN and *2XN) enzyme activity. Information regarding the use of CYP2D6 inhibitors was recorded. Results: A total of 591 cases were matched to 1,126 controls and DNA was genotyped in more than 97%. In patients treated with tamoxifen, there was no association of CYP2D6 genotype [OR (extensive/poor metabolizer): 0.90; 95% CI: 0.46–1.74, P = 0.74), use of a potent CYP2D6 inhibitor (OR 0.92; 95% CI: 0.575–1.486), or CYP2D6 metabolizer status (OR 1.03; 95% CI: 0.669–1.607) with breast cancer occurrence. Likewise, there was no association between any CYP2D6 metabolism parameter with breast cancer events in raloxifene-treated patients. Conclusion: In the NSABP P1 and P2 clinical trials, alterations in CYP2D6 metabolism are not associated with either tamoxifen or raloxifene efficacy. Clin Cancer Res; 17(21); 6944–51. ©2011 AACR.

Abstract Background: The third-generation aromatase inhibitors (AIs) are established adjuvant therapy for postmenopausal women with early stage breast cancer. Musculoskeletal adverse events (MS AE) have become increasingly recognized as an issue for many women receiving AIs resulting in decreased quality of life and, in some instances, discontinuation of the AI. MA.27 is an ongoing randomized clinical trial comparing anastrozole with exemestane that accrued 7576 patients (pts), including 6827 from North America, about 79% of whom provided DNA and consent for genetic testing. A review of withdrawals related to adverse events revealed that the most common cause was MS AE. The hypothesis for our study was that a genome-wide association (GWA) case control study would identify single nucleotide polymorphisms (SNPs) associated with MS AE.Methods: Cases were defined as pts with grade 3 or 4 MS AE (NCI CTCAE v3.0) or who went off treatment for any grade MS AE within the first 2 years and were matched to 2 controls based on treatment arm, length of follow-up, prior chemotherapy, age and prior celecoxib (at least 3 months prior to MS AE). The study was restricted to the 94% of pts self-identified as Caucasian. Genotyping was performed with the Illumina Human610-Quad. Genotype quality control prior to analysis required call rates ≥98%. Eigenstrat analyses were performed to control for population stratification.Results: The GWA study included 293 cases and 586 controls with 6 pts (0.7%) (4 cases, 2 controls) excluded for call rates <98%. 551,358 SNPs were used in the analyses after removing 29,439 for a minor allele frequency <1% and 82 with Hardy Weinberg Equilibrium p-values <10-6. Eigenvectors did not impact the results when used as covariates. Four SNPs with the lowest p-values (3.3-8.2 x 10-6, Armitage test), 3 in high LD (r2>.8) on chromosome (ch) 14 (rs7158782, rs7159713, rs2369049) and 1 on the X ch (rs6637820) were identified. The gene closest (4000-7000 bp) to the 3 SNPs on ch 14 was T-Cell Leukemia 1A (TCL1A), encoding a protein known to augment AKT kinase activity. Expression array data from lymphoblastoid cell lines from 300 subjects of 3 ethnicities revealed that all 3 SNPs were associated with decreased TCL1A expression after adjusting for ethnicity using 2 probe sets (p=0.006 to 0.03). Electromobility shift assay (EMSA) showed that 2 of these SNPs (rs7158782, rs7159713) displayed a shift and the former SNP showed substantially less binding for variant than wild type. The gene closest to the X ch SNP was Immunoglobulin Superfamily 1 (IGSF1), but EMSA did not show a shift for this SNP.Conclusions: MS AEs represent a major impediment to optimal use of AIs in women with breast cancer and this GWA study identified SNPs on ch 14 that provide a focus for further research to identify pts at risk for, and means to ameliorate, this adverse event.(Supported in part by NIH grants U01GM61388, U01GM63173, P50CA116201 and U10CA77202) Citation Information: Cancer Res 2009;69(24 Suppl):Abstract nr 15.