Citations

Abstract Selective Estrogen Receptor Modulators (SERMs) such as tamoxifen and raloxifene can reduce the occurrence of breast cancer in high risk women by approximately 50%. This presentation will describe the use of a genome-wide association study (GWAS) to identify novel genes associated with breast cancer risk during SERM prevention therapy. Those results were pursued by functional genomic experiments that defined mechanisms responsible for the associations observed during the GWAS. These observations not only moved “beyond biomarkers” to function and underlying mechanism, but they also provided novel insight into mechanisms responsible for individual variation in SERM response phenotypes. Citation Format: Richard Weinshilboum. SERM breast cancer pharmaco-genomics: Beyond biomarkers. [abstract]. In: Proceedings of the Eleventh Annual AACR International Conference on Frontiers in Cancer Prevention Research; 2012 Oct 16-19; Anaheim, CA. Philadelphia (PA): AACR; Cancer Prev Res 2012;5(11 Suppl):Abstract nr PL04-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

Background— The primary role of natriuretic peptide receptor-3 (NPR3) or NPR-C is in the clearance of natriuretic peptides that play an important role in modulating intravascular volume and vascular tone. Genetic variation in NPR3 has been associated with variation in blood pressure and obesity. Despite the importance of NPR3, sequence variation in the gene has not been addressed using DNA from different ethnic populations. We set out to identify and functionally characterize genetic variation in NPR3 in 3 ethnic groups. Methods and Results— DNA samples from 96 European American, 96 African American, and 96 Han Chinese American healthy subjects were used to resequence NPR3 exons, splice junctions, and flanking regions. We identified 105 polymorphisms, 50 of which were novel, including 8 nonsynonymous single-nucleotide polymorphisms, 7 were novel. Expression constructs were created for the nonsynonymous single-nucleotide polymorphisms. HEK293 cells were transfected with constructs for wild type and variant allozymes; and recombinant proteins were measured by quantitative Western blot analysis. The most significant change in NPR3 protein was observed for the Arg146 variant allozyme, with 20% of wild-type protein, primarily because of autophagy-dependent degradation. NPR3 structural modeling confirmed that the Arg146 variant protein was not compatible with wild-type conformation and could result in protein misfolding or instability. Conclusions— Multiple novel NPR3 genetic polymorphisms were identified in 3 ethnic groups. The Arg146 allozyme displayed a significant decrease in protein quantity because of degradation mediated predominantly by autophagy. This genetic variation could have a significant effect on the metabolism of natriuretic peptides with potential clinical implications.

Abstract Purpose: Mammographic breast density (MBD) is decreased by tamoxifen, but the effect of aromatase inhibitors is less clear. Experimental Design: We enrolled early-stage postmenopausal patients with breast cancer initiating adjuvant aromatase inhibitor therapy and ascertained mammograms before and at an average 10 months of aromatase inhibitor therapy. We matched cases to healthy postmenopausal women (controls) from a large mammography screening cohort on age, baseline body mass index, baseline MBD, and interval between mammograms. We estimated change in MBD using a computer-assisted thresholding program (Cumulus) and compared differences between cases and matched controls. Results: In predominantly White women (96%), we found 14% of the 387 eligible cases had a MBD reduction of at least 5% after an average of 10 months of aromatase inhibitor therapy. MBD reductions were associated with higher baseline MBD, aromatase inhibitor use for more than 12 months, and prior postmenopausal hormone use. Comparing each case with her matched control, there was no evidence of an association of change in MBD with aromatase inhibitor therapy [median case–control difference among 369 pairs was −0.1% (10th and 90th percentile: −5.9%, 5.2%) P = 0.51]. Case–control differences were similar by type of aromatase inhibitor (P's 0.41 and 0.56); prior use of postmenopausal hormones (P = 0.85); baseline MBD (P = 0.55); and length of aromatase inhibitor therapy (P = 0.08). Conclusions: In postmenopausal women treated with aromatase inhibitors, 14% of cases had a MBD reduction of more than 5%, but these decreases did not differ from matched controls. These data suggest that MBD is not a clinically useful biomarker for predicting the value of aromatase inhibitor therapy in White postmenopausal women. Clin Cancer Res; 19(8); 2144–53. ©2013 AACR.

Genome‐wide association studies (GWAS) of complex traits have generated many association signals for single nucleotide polymorphisms (SNPs). To understand the underlying causal genetic variant(s), focused DNA resequencing of targeted genomic regions is commonly used, yet the current cost of resequencing limits sample sizes for resequencing studies. Information from the large GWAS can be used to guide choice of samples for resequencing, such as the SNP genotypes in the targeted genomic region. Viewing the GWAS tag‐SNPs as imperfect surrogates for the underlying causal variants, yet expecting that the tag‐SNPs are correlated with the causal variants, a reasonable approach is a two‐phase case‐control design, with the GWAS serving as the first‐phase and the resequencing study serving as the second‐phase. Using stratified sampling based on both tag‐SNP genotypes and case‐control status, we explore the gains in power of a two‐phase design relative to randomly sampling cases and controls for resequencing (i.e., ignoring tag‐SNP genotypes). Simulation results show that stratified sampling based on both tag‐SNP genotypes and case‐control status is not likely to have lower power than stratified sampling based only on case‐control status, and can sometimes have substantially greater power. The gain in power depends on the amount of linkage disequilibrium between the tag‐SNP and causal variant alleles, as well as the effect size of the causal variant. Hence, the two‐phase design provides an efficient approach to follow‐up GWAS signals with DNA resequencing.

Abstract Purpose: Controversy exists about CYP2D6 genotype and tamoxifen efficacy. Experimental Design: A matched case–control study was conducted using the Austrian Breast and Colorectal Cancer Study Group Trial 8 (ABCSG8) that randomized postmenopausal women with estrogen receptor (ER)-positive breast cancer to tamoxifen for 5 years (arm A) or tamoxifen for 2 years followed by anastrozole for 3 years (arm B). Cases had disease recurrence, contralateral breast cancer, second non–breast cancer, or died. For each case, controls were identified from the same treatment arm of similar age, surgery/radiation, and tumor–node—metastasis (TNM) stage. Genotyping was conducted for alleles associated with no (PM; *3, *4, *6), reduced (IM; *10, and *41), and extensive (EM: absence of these alleles) CYP2D6 metabolism. Results: The common CYP2D6*4 allele was in Hardy–Weinberg equilibrium. In arm A during the first 5 years of therapy, women with two poor alleles [PM/PM: OR, 2.45; 95% confidence interval (CI), 1.05–5.73, P = 0.04] and women with one poor allele (PM/IM or PM/EM: OR, 1.67; 95% CI, 0.95–2.93; P = 0.07) had a higher likelihood of an event than women with two extensive alleles (EM/EM). In years 3 to 5 when patients remained on tamoxifen (arm A) or switched to anastrozole (arm B), PM/PM tended toward a higher likelihood of a disease event relative to EM/EM (OR, 2.40; 95% CI, 0.86–6.66; P = 0.09) among women on arm A but not among women on arm B (OR, 0.28; 95% CI, 0.03–2.30). Conclusion: In ABCSG8, the negative effects of reduced CYP2D6 metabolism were observed only during the period of tamoxifen administration and not after switching to anastrozole. Clin Cancer Res; 19(2); 500–7. ©2012 AACR.

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.