Wang, Liewei

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.

Abstract We have recently showed that the immunophilin FKBP5 (also known as FKBP5) as a scaffolding protein that can enhance PHLPP-AKT interaction and facilitate PHLPP-mediated dephosphorylation of AKT-Ser473, negatively regulating AKT activation in vitro. Therefore, FKBP5 might function as a tumor suppressor and levels of FKBP5 would affect cell response to chemotherapy. In the current study, we took a step forward by using pancreatic caner xenograft mice model to show that FKBP5−/− promotes tumor growth and renders resistance to gemcitabine and addition of AKT inhibitors have more significant sensitization effect on FKBP5−/− mice to gemcitabine treatment. The combination of triciribine and gemcitabine resulted in cooperative inhibition of pancreatic cancer cell proliferation and the growth of xenografts in vivo without increasing toxicity, as well as cooperative inhibition of AKT signaling pathways at concentrations that were ineffective as individual agents. These results provide evidence for future clinical studies designed to tailor the therapy based on the biomarker. 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 C15.

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.