Two cytidine analogues, gemcitabine (dFdC) and 1-β-d-arabinofuranosylcytosine (AraC), show significant therapeutic effect in a variety of cancers. However, response to these drugs varies widely. Evidence from tumor biopsy samples shows that expression levels for genes involved in the cytidine transport, metabolism, and bioactivation pathway contribute to this variation in response. In the present study, we set out to test the hypothesis that variation in gene expression both within and outside of this “pathway” might influence sensitivity to gemcitabine and AraC. Specifically, Affymetrix U133 Plus 2.0 GeneChip and cytotoxicity assays were performed to obtain basal mRNA expression and IC50 values for both drugs in 197 ethnically defined Human Variation Panel lymphoblastoid cell lines. Genes with a high degree of association with IC50 values were involved mainly in cell death, cancer, cell cycle, and nucleic acid metabolism pathways. We validated selected significant genes by performing real-time quantitative reverse transcription-PCR and selected two representative candidates, NT5C3 (within the pathway) and FKBP5 (outside of the pathway), for functional validation. Those studies showed that down-regulation of NT5C3 and FKBP5 altered tumor cell sensitivity to both drugs. Our results suggest that cell-based model system studies, when combined with complementary functional characterization, may help to identify biomarkers for response to chemotherapy with these cytidine analogues. [Cancer Res 2008;68(17):7050–8]
Cancer patients show large individual variation in their response to chemotherapeutic agents. Gemcitabine (dFdC) and AraC, two cytidine analogues, have shown significant activity against a variety of tumors. We previously used expression data from a lymphoblastoid cell line-based model system to identify genes that might be important for the two drug cytotoxicity. In the present study, we used that same model system to perform a genome-wide association (GWA) study to test the hypothesis that common genetic variation might influence both gene expression and response to the two drugs. Specifically, genome-wide single nucleotide polymorphisms (SNPs) and mRNA expression data were obtained using the Illumina 550K® HumanHap550 SNP Chip and Affymetrix U133 Plus 2.0 GeneChip, respectively, for 174 ethnically-defined “Human Variation Panel” lymphoblastoid cell lines. Gemcitabine and AraC cytotoxicity assays were performed to obtain IC50 values for the cell lines. We then performed GWA studies with SNPs, gene expression and IC50 of these two drugs. This approach identified SNPs that were associated with gemcitabine or AraC IC50 values and with the expression regulation for 29 genes or 30 genes, respectively. One SNP in IQGAP2 (rs3797418) was significantly associated with variation in both the expression of multiple genes and gemcitabine and AraC IC50. A second SNP in TGM3 (rs6082527) was also significantly associated with multiple gene expression and gemcitabine IC50. To confirm the association results, we performed siRNA knock down of selected genes with expression that was associated with rs3797418 and rs6082527 in tumor cell and the knock down altered gemcitabine or AraC sensitivity, confirming our association study results. These results suggest that the application of GWA approaches using cell-based model systems, when combined with complementary functional validation, can provide insights into mechanisms responsible for variation in cytidine analogue response.