Citations

Abstract Background: Candidate gene and candidate pathway-based pharmacogenomic studies have been used to investigate the inherited variability of prognosis of lung cancer treated with platinum-based chemoptherapy. In order to make it possible to conduct an ‘unbiased’ query across the entire genome, a genome-wide association study (GWAS) in 283 lymphoblastoid cell lines (LCLs) was carried out to help identify candidate SNPs associated with variation in cisplatin cytotoxicity. Furthermore, we set out to test the hypothesis that genetic variation resulting in differences in cisplatin cytotoxicity might influence the overall survival of patients with lung cancer treated with platinum-based antineoplastic agents. Material and Methods: We applied ethnically diverse Coriell “Human Variation Panel” LCLs for GWAS. Then, we picked 168 top significant SNPs selected from GWAS in 283 LCLs and genotyped these SNPs in 1183 lung cancer patients (222 SCLC and 961 NSCLC) who received platinum-based chemotherapy. Genotyping was performed using the Illumina Golden Gate platform. The association of SNPs with lung cancer overall survival was analyzed by Cox regression model after adjusting for disease stage. Functional validation of candidate genes was performed by siRNA screening in NSCLC cell lines CRL5872, H1299 and SCLC cell line CRL5823. Results: A total of 157 SNPs were successfully genotyped for 1183 lung cancer patient samples. The Cox regression analysis indicated that 9 and 10 different SNPs were associated with overall survival for NSCLC and SCLC patients, respectively (p value <0.05). Fifteen genes including 7 genes located in the region of 200kb up or downstream of 4 SNPs (rs11169748, rs2440915, rs5952066 and rs7620841) and 8 genes for which the expression was correlated with three SNPs (rs11169748, rs2440915 and rs5952066) were selected for siRNA screening. Knockdown of DAPK3 and METTL6, for which the expression levels were correlated with the SNPs rs11169748 and rs2440915, significantly decreased cisplatin sensitivity in lung cancer cell lines CRL5872, H1299 and CRL5823. Conclusions: This series of clinical and complementary laboratory-based functional studies identified several candidate genes and SNPs that might help predict treatment outcomes of platinum-based therapy in lung cancer. 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 1285. doi:10.1158/1538-7445.AM2011-1285

Abstract Variation in drug response results from a combination of factors that include differences in gender, ethnicity, and environment, as well as genetic variation that may result in differences in mRNA and protein expression. This article presents two integrative analytic approaches that make use of both genome-wide SNP and mRNA expression data available on the same set of subjects: a step-wise integrative approach and a comprehensive analysis using sparse canonical correlation analysis (SCCA). In addition to applying standard SCCA, we present a novel modification of SCCA which allows different weighting for the various pair-wise relationships in the SCCA. These integrative approaches are illustrated with both simulated data and data from a pharmacogenomic study of the drug gemcitabine. Results from these analyses found little overlap in terms of genes detected, possibly detecting different biological mechanisms. In addition, we found the proposed weighted SCCA to outperform its unweighted counterpart in detecting associations between the genomic features and phenotype. Further research is needed to develop and assess new integrative methods for pharmacogenomic studies, as these types of analyses may uncover novel insights into the relationship between genomic variation and drug response.

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 Taxane, a class of common anticancer agent, is widely used to treat many cancers. Among the taxanes, paclitaxel and docetaxel are the most commonly used agents. The mechanism of action for taxane drugs is through stabilizing microtubules by enhancing tubulin assembly and against depolymerization, thereby disrupting the dynamics of the microtubule. Large inter-individual variation has been observed in response to taxane treatment. Therefore, in the present study, we conducted a genome-wide association (GWA) study to identify genetic polymorphisms that might contribute to variation in response to taxanes drugs using the “Human Variation Panel” lymphoblastoid cell line (LCL) model system that consists of 287 ethnically defined cell lines. Method Genome-wide SNP and basal gene expression data were obtained for 287 LCLs using Illumina HumanHap 550K, 510S BeadChips, Affymetix SNP Array 6.0 and Affymetrix U133 plus 2.0 GeneChip. Cytotoxicity assay was also performed in the same cell lines for both paclitaxel and docetaxel. Candidate genes were selected based on an integrated analysis of genome-wide SNPs, expression array data and IC50 values (a cytotoxicity phenotype indicating the drug concentration required for 50% inhibition of maximal cell growth) for paclitaxel and docetaxel. SiRNA screening was performed to functionally validate those top candidates, followed by MTS assay in a breast cancer cell line, MDA-MB231. Results Of the 54,000 probe sets considered, 390 probe sets, representing 272 unique annotated genes, were found to be associated with docetaxel cytotoxicity (IC50 value) with p value <10−3, while 176 probe sets, representing 124 unique genes, were associated with paclitaxel IC50. Of these probe sets, 50 probe sets were associated with IC50 for both drugs. Analysis of 1.3 million SNPs identified 11 loci, defined as a region containing at least 1 SNP of p value <10−5 and 1 SNP of p value <10−4 within 50kb, associated with docetaxel IC50, while 7 loci were associated with paclitaxel IC50. One locus was found to be associated with both drugs. Previous studies suggested that gene expression is regulated in a cell-type specific manner. Therefore, based on our analysis and biological implication, we selected 16 top candidate genes that were expressed in our LCLs to perform functional validation studies using siRNA, followed by MTS assay in a breast cancer cell line, MDA-MB231. Functional studies showed that knockdown of SLC44A1 and HSPA5 significantly desensitized the cancer cell to both taxanes treatment. Conclusion GWAS using LCL model system could provide potential pharmacogenomic candidates for future clinical translational studies. Functional validation in tumor cell lines may help determine the genomic mechanism contributing to the inter-individual variation of taxane response. 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 2599.

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.