Citations

Abstract The FKBP5 gene encodes the FKBP51 immunophilin, a negative regulator of AKT. Variability in level of FKBP5 expression is a major factor contributing to variation in response to several classes of chemotherapeutic agents including gemcitabine, first line treatment for pancreatic cancer. Genetic variation of FKBP5 could influence its expression, function, and, ultimately, treatment response of pancreatic cancer. We set out to comprehensively study the role of genetic variation in FKBP5 identified by Next Generation DNA resequencing in response to gemcitabine treatment of pancreatic cancer by utilizing both tumor and germline DNA samples from 43 pancreatic cancer patients, including 19 paired tumor-normal samples. Genotype-phenotype association studies with overall survival (in tumor only) as well as with FKBP5 gene expression were performed using the same samples in which Next Generation DNA resequencing had been performed, followed by functional genomic studies of top selected single nucleotide polymorphisms (SNPs). In-depth resequencing identified 404 FKBP5 SNPs in both normal and tumor DNA patient samples. SNPs with the smallest p-values for association with FKBP5 survival or expression were subjected to functional genomic studies. Electromobility shift assay (EMSA) showed that the rs73748206 “T” allele altered DNA-protein binding patterns, consistent with a significantly increased reporter gene activity. The effect of rs73748206 was also confirmed by mRNA expression studies using lymphoblastoid cell lines (LCL) derived from the same patients from whom DNA used for resequencing had been obtained. This comprehensive FKBP5 resequencing study provides insight into the role of inheritance in variation in response to gemcitabine therapy of pancreatic cancer. Citation Format: Katarzyna A. Ellsworth, Bruce W. Eckloff, Liang Li, Irene Moon, Brooke L. Fridley, Greg D. Jenkins, Erin Carlson, Abra Brisbin, Ryan Abo, William Bamlet, Gloria Petersen, Eric D. Wieben, Liewei Wang. Contribution of FKBP5 genetic variation to gemcitabine treatment and survival in pancreatic cancer. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 2209. doi:10.1158/1538-7445.AM2013-2209

Abstract Integrative genomics has the potential to uncover relevant loci, as clinical outcome and response to chemotherapies are most likely not due to a single gene (or data type) but rather a complex relationship involving genetic variation, mRNA, DNA methylation, and copy number variation. In addition to this complexity, many complex phenotypes are thought to be controlled by the interplay of multiple genes within the same molecular pathway or gene set (GS). To address these two challenges, we propose an integrative gene set analysis approach and apply this strategy to a cisplatin (CDDP) pharmacogenomics study involving lymphoblastoid cell lines for which genome-wide SNP and mRNA expression data was collected. Application of the integrative GS analysis implicated the role of the RNA binding and cytoskeletal part GSs. The genes LMNB1 and CENPF, within the cytoskeletal part GS, were functionally validated with siRNA knockdown experiments, where the knockdown of LMNB1 and CENPF resulted in CDDP resistance in multiple cancer cell lines. This study demonstrates the utility of an integrative GS analysis strategy for detecting novel genes associated with response to cancer therapies, moving closer to tailored therapy decisions for cancer patients.

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

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 DNA double-strand breaks (DSBs) are highly cytotoxic lesions caused by genotoxin such as ionizing radiation (IR). Following DSB, cells activate the DNA-damage response pathway by activating protein kinases such as ATM, ATR and DNA-PK, which then trigger histone H2AX phosphorylation and the accumulation of proteins such as Mediator of DNA damage checkpoint protein 1 (MDC1), 53BP1, and BRCA1 into IR-induced foci. This process initiates the repair of DSB. Previously we have identified a novel biomarker for radiation response, Bora, by using a genome-wide association approach using data generated with 300 human lymphoblastoid cell liens. Bora is known as a cofactor of Aurora-A, which is involved in mitosis regulation. During G2, Aurora-A, in complex with Bora, phosphorylates Plk1 at Thr-210. Activated Plk1 phosphorylates Cdc25C and Wee1, which induces activation of cyclin B-cdk1 complexes and promotes mitotic entry. In this study, we want to determine the underlying mechanisms by which Bora regulates response to IR. We found that knockdown of Bora in human cancer cell lines including HupT3, MIAPaCa-2 and Hela (siBora) exhibit increased clonogenic survival rate in response to IR, an observation that is opposite from the effect of Plk1 knockdown, suggesting a different mechanism involved. Knockdown of Bora also desensitize cancer cells to chemotherapeutic agents including gemcitabine and oxiplatin. This enhanced resistance to genotoxin is primarily due to increased DNA DSB repair. More than 90% of DSBs were repaired in siBora cells compared with 60% in control cells 8 hours after cell exposure to 2 Gy IR. Upon IR, more cells with Bora knockdown retained in G2-M phase compared with control. Furthermore, Bora directly interacts with BRCT domain of MDC1 via its C terminus, and over expression of Bora can abolish MDC1 foci formation induced by IR. All of these observations suggested a role of Bora in DNA damage response pathway and in turn, in chemo-radio resistance. Indeed, overexpression of the full length or C-terminus of Bora sensitized cancer cells to IR treatment, while the N-terminal Bora did not have significant effect on IR response. In summary, we have identified a novel function of Bora involved in DNA damage and repair pathway, which indicates a broad role of Bora in both tumorgenesis and chemo-radio resistance. 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 2660. doi:10.1158/1538-7445.AM2011-2660

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.