Batzler, Anthony

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: 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 Purpose: Inherited variability in the prognosis of lung cancer patients treated with platinum-based chemotherapy has been widely investigated. However, the overall contribution of genetic variation to platinum response is not well established. To identify novel candidate single nucleotide polymorphisms (SNP)/genes, we carried out a genome-wide association study (GWAS) for cisplatin cytotoxicity by using lymphoblastoid cell lines (LCL), followed by an association study of selected SNPs from the GWAS with overall survival (OS) in lung cancer patients. Experimental Design: A GWAS for cisplatin was conducted with 283 ethnically diverse LCLs. A total of 168 top SNPs were genotyped in 222 small cell lung cancer (SCLC) and 961 non-SCLC (NSCLC) patients treated with platinum-based therapy. Association of the SNPs with OS was determined by using the Cox regression model. Selected candidate genes were functionally validated by siRNA knockdown in human lung cancer cells. Results: Among 157 successfully genotyped SNPs, 9 and 10 SNPs were top SNPs associated with OS for patients with NSCLC and SCLC, respectively, although they were not significant after adjusting for multiple testing. Fifteen genes, including 7 located within 200 kb up or downstream of the 4 top SNPs and 8 genes for which expression was correlated with 3 SNPs in LCLs were selected for siRNA screening. Knockdown of DAPK3 and METTL6, for which expression levels were correlated with the rs11169748 and rs2440915 SNPs, significantly decreased cisplatin sensitivity in lung cancer cells. Conclusions: This series of clinical and complementary laboratory-based functional studies identified several candidate genes/SNPs that might help predict treatment outcomes for platinum-based therapy of lung cancer. Clin Cancer Res; 17(17); 5801–11. ©2011 AACR.

Abstract Background: NSABP P1 and P2 are the largest prevention trials of tamoxifen and raloxifene in women at high risk of developing breast cancer. Both agents have been shown to decrease, but not eliminate breast cancer risk in these women. The goal of our study was to identify single nucleotide polymorphisms (SNPs) associated with development of breast cancer in high-risk women treated with tamoxifen or raloxifene. Methods: We performed a nested matched case-control genome-wide association study (GWAS) in subjects entered on the P2, which included only postmenopausal women, and subjects at least 50 years old entered on the P1 phase III prevention trials who were treated with tamoxifen or raloxifene. Cases were defined as subjects who developed invasive breast cancer or ductal carcinoma in situ and each was matched to two controls based on trial (P1, P2), age at trial entry, Gail risk score, history of lobular carcinoma in situ, and time on study. The study was restricted to the 94% of subjects self-identified as Caucasian. Genotyping was performed with the Illumina Human610-Quad chip. Eigenstrat analyses were performed to control for population stratification. SNPs were imputed within 200 Kb on either side of regions containing SNPs with p ≤ 3E-05. CYP2D6 genotyping was performed on all subjects. Results: The GWAS included 592 cases and 1171 controls. 547,356 SNPs were used in the analyses after removing 14,019 for call rates less than 98%, 30,560 for minor allele frequencies <1%, and 301 with Hardy-Weinberg equilibrium p-value <1E-06. Eigenvectors did not alter the results when used as covariates. Eleven genotyped SNPs with p-values ranging from 2.12E-06 to 9.98E-06 (conditional logistic regression with eigenvector adjustment) were identified; and 22 SNPs with p < 4E-05 were identified by imputation around these SNPs on chromosomes (Chr) 4, 8, 9, 13, 16 and 22, which were then genotyped with the Invader platform. Following imputation and genotyping, initial functional genomic studies have focused on 8 SNPs on Chr 16 (p-values 1.81E-06 to 9.55E-06, with the lowest p-value for an imputed and fine mapped SNP), all of which were in ZNF423, a gene known to encode a protein that associates with SMAD family members. In lymphoblastoid cell lines genotyped for the 2 top SNPs in ZNF423 (rs8060157 and rs11076499) and stably transfected with estrogen receptor-α, the variant and wild type SNPs displayed striking differences in estradiol-induced expression of ZNF423, SMAD3, BRCA1, and BRCA2. Conclusions: This GWAS revealed SNPs associated with the development of breast cancer in high-risk women treated with tamoxifen or raloxifene. Genome-wide significance was not achieved, but SNPs in ZNF423 identified on Chr 16 were associated with differential expression of ZNF423, SMAD3, BRCA1, and BRCA2 in response to estradiol in vitro, providing direction for future research. (Supported in part by NIH grants U01GM61388, U01GM63173, P50CA116201, U10CA77202, U10CA37377, U10CA69974, U24CA114732, and the Biobank Japan Project funded by the Ministry of Education, Culture, Sports, Science and Technology, Japan) Citation Information: Cancer Res 2010;70(24 Suppl):Abstract nr PD05-02.

Abstract Background Radiation therapy is used to treat half of the all cancer patients. Response to radiation therapy varies widely, in terms of both response and toxicity related events. Previous studies suggested almost 80% of the variation in normal tissue response to radiation treatment was likely to be due to genetic factors. Therefore, we performed a genome-wide association study (GWAS) to identify biomarkers that might help predict radiation response using a the “Human Variation Panel” lymphoblastoid cell lines (LCLs) that consists 287 ethnically defined cell lines. Method Basal gene expression level, 1.3 million genome-wide SNP data were obtained for 287 human LCLs using Affymetrix U133 plus 2.0 GeneChip, Illumina HumanHap 550K, 510S BeadChips and Affymetix SNP Array 6.0. Radiation cytotoxicity assay was also performed to obtain the area under radiation dose response curve (AUC), a summary measure for cytotoxicity for the same cell lines. Functional validation of candidate genes, selected based on the integrated analysis of SNP, expression and AUC data, was performed in multiple cancer cell lines using specific siRNAs followed by MTS and colony formation assays. Results 27 loci, containing at least 2 SNPs of p values <10−4 within 50kb were associated with radiation AUC. 270 expression probe sets were associated with radiation AUC (p< 10−3). The integrated analysis of 27 loci, 54,000 basal expression probe sets and radiation AUC identified 50 unique SNPs associated with AUC were also associated with gene expression of 39 unique annotated genes. In addition, these 39 genes were associated with radiation AUC with p values <10−3. Based on the GWAS and integrated analysis of SNP, expression array and radiation AUC data, 22 genes were selected for functional validation studies using siRNA knock down followed by MTS and colony formation assays in multiple tumor cell lines. C13orf34 (Bora), MAD2L1, PLK4, TPD52 significantly altered the radiation sensitivity in at least 2 cancer cell lines for both MTS and colony formation assays (p<0.05). Knock down of DEPDC1B, a gene containing 2 SNPs at the 3′-end that was significantly related with Bora gene expression level (p<10−4) significantly desensitized cancer cells to radiation treatment. Both microarray data and QRT-PCR analysis indicated that expression levels for those two genes were significantly correlated. Conclusion Our results using the lymphoblastoid cell lines might help identify novel biomarkers that contribute to variation in response to radiation therapy and help understand the mechanisms underlying the association. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 3602.

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.

Purpose We performed a case-control genome-wide association study (GWAS) to identify single nucleotide polymorphisms (SNPs) associated with musculoskeletal adverse events (MS-AEs) in women treated with aromatase inhibitors (AIs) for early breast cancer. Patients and Methods A nested case-control design was used to select patients enrolled onto the MA.27 phase III trial comparing anastrozole with exemestane. Cases were matched to two controls and were defined as patients with grade 3 or 4 MS-AEs (according to the National Cancer Institute's Common Terminology Criteria for Adverse Events v3.0) or those who discontinued treatment for any grade of MS-AE within the first 2 years. Genotyping was performed with the Illumina Human610-Quad BeadChip. Results The GWAS included 293 cases and 585 controls. A total of 551,358 SNPs were analyzed, followed by imputation and fine mapping of a region of interest on chromosome 14. Four SNPs on chromosome 14 had the lowest P values (2.23E-06 to 6.67E-07). T-cell leukemia 1A (TCL1A) was the gene closest (926-7000 bp) to the four SNPs. Functional genomic studies revealed that one of these SNPs (rs11849538) created an estrogen response element and that TCL1A expression was estrogen dependent, was associated with the variant SNP genotypes in estradiol-treated lymphoblastoid cells transfected with estrogen receptor alpha and was directly related to interleukin 17 receptor A (IL17RA) expression. Conclusion This GWAS identified SNPs associated with MS-AEs in women treated with AIs and with a gene (TCL1A) which, in turn, was related to a cytokine (IL17). These findings provide a focus for further research to identify patients at risk for MS-AEs and to explore the mechanisms for these adverse events.

Abstract Aromatase inhibitors play a prominent role in the management of postmenopausal women with endocrine-sensitive breast cancer, but there is large variability in both efficacy and tolerability. The purpose of our study was to define interindividual variation in anastrozole metabolism and pharmacodynamics among patients treated with the approved daily dose of 1 mg in a standard practice setting as adjuvant therapy for resected early breast cancer. This study was performed in 191 women in whom pretreatment and during anastrozole plasma concentrations of estrone (E1), estradiol (E2), estrone conjugates, androstenedione, and testosterone were determined and correlated with plasma concentrations of anastrozole and anastrozole metabolites. There were large interindividual variations in plasma anastrozole and anastrozole metabolite concentrations, as well as pretreatment and postdrug plasma E1, E2, and E1 conjugate and estrogen precursor (androstenedione and testosterone) concentrations. E1 and E2 concentrations were below the lower limit of quantitation (LLQ) in most patients after anastrozole therapy (83% for both), but those with detectable concentrations had a broad range (1.58–45.2 and 0.635–97.0 pg/mL, respectively). E1 conjugates after anastrozole therapy were above the LLQ in most patients (93%), with wide interpatient variability (3.50–2,990 pg/mL). Two patients seemed to extensively metabolize anastrozole and failed to display substantial decreases in estrogens. Acknowledging the potential factor of variable compliance, our results showed large interindividual variation in anastrozole metabolism and its effect on circulating estrogens in postmenopausal patients. These findings may have implications with regard to efficacy and adverse events and may indicate the need to “individualize” therapy with this drug. Cancer Res; 70(8); 3278–86. ©2010 AACR.

Abstract Background: Lung cancer is commonly treated with platinum compounds. The “glutathione pathway” participates in the metabolism of platinum compounds. We set out to test the hypotheses that single nucleotide polymorphisms (SNPs) or copy number polymorphisms for genes within the glutathione pathway might influence survival in lung cancer patients treated with these drugs. Methods: Germline DNA samples from 973 lung cancer patients were genotyped for 290 glutathione pathway SNPs. GSTT1 copy number was also assayed. We determined the association of these polymorphisms with survival for lung cancer patients, followed by functional genomic validation. Results: We observed suggestive associations between survival and GSTT1 copy number (P = 0.017), and GSTA5, GSTM4, and ABCC4 SNPs, adjusted for covariates (P = 0.018, 0.002, and 0.002, respectively) or not (P = 0.005, 0.011, and 0.002). One hundred lymphoblastoid cell lines were then treated with cisplatin, and IC50 values were significantly associated with the GSTM4 SNP (P = 0.019). Furthermore, GSTM4, GSTT1, and ABCC4 overexpression significantly decreased cisplatin sensitivity in lung cancer and HEK293T cell lines. Conclusions: These results suggest that GSTM4 polymorphisms are biomarkers for the prediction of cisplatin response. ABCC4 polymorphisms, as well as GSTT1 copy number, may also help to predict cisplatin response, but further validation is required. These results represent a step toward the individualized chemotherapy of lung cancer. Cancer Epidemiol Biomarkers Prev; 19(3); 811–21

Abstract Aromatase (CYP19) is a critical enzyme in estrogen biosynthesis and aromatase inhibitors (AI) are employed widely for endocrine therapy in postmenopausal women with breast cancer. We hypothesized that single nucleotide polymorphisms (SNPs) in the CYP19 gene may alter the effectiveness of AI therapy in the neoadjuvant setting. Genomic DNA was obtained for sequencing from 52 women pre-AI and post-AI treatment in this setting. Additionally, genomic DNA obtained from 82 samples of breast cancer and 19 samples of normal breast tissue was subjected to resequencing. No differences in CYP19 sequence were observed between tumor and germ-line DNA in the same patient. A total of 48 SNPs were identified including 4 novel SNPs when compared with previous resequencing data. For genotype-phenotype association studies, we determined the levels of aromatase activity, estrone, estradiol, and tumor size in patients pre-AI and post-AI treatment. We defined two tightly linked SNPs (rs6493497 and rs7176005 in the 5'-flanking region of CYP19 exon 1.1) that were significantly associated with a greater change in aromatase activity after AI treatment. In a follow-up study of 200 women with early-stage breast cancer who were treated with adjuvant anastrozole, these same two SNPs were also associated with higher plasma estradiol levels in patients pre-AI and post-AI treatment. Electrophoretic mobility shift and reporter gene assays confirmed likely functional effects of these two SNPs on transcription of CYP19. Our findings indicate that two common genetic polymorphisms in the aromatase gene CYP19 vary the response of breast cancer patients to aromatase inhibitors. Cancer Res; 70(1); 319–28