Gemcitabine and arabinosylcytosin pharmacogenomics: genome-wide association and drug response biomarkers - PubMed (original) (raw)
Gemcitabine and arabinosylcytosin pharmacogenomics: genome-wide association and drug response biomarkers
Liang Li et al. PLoS One. 2009.
Abstract
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(R) 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 IC(50) values for the cell lines. We then performed GWA studies with SNPs, gene expression and IC(50) of these two drugs. This approach identified SNPs that were associated with gemcitabine or AraC IC(50) 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 IC(50). 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.
Conflict of interest statement
Competing Interests: The authors have declared that no competing interests exist.
Figures
Figure 1. Genome-wide SNP association with drug IC50.
(A). Genome-wide SNP-gemcitabine IC50 association data. (B). Genome-wide SNP-AraC IC50 association data. Genome-wide SNPs were obtained for all 171 cell lines using Illumina 550K Beadchips. Genome-wide SNP association analysis was then performed with gemcitabine IC50 as the drug-response phenotype (n = 171). The y-axis is the −log10 (p-value) for the SNP. SNPs are plotted on the x-axis on the basis of their chromosomal locations. The “cut-off” p-value of 0.001 is highlighted with a red line.
Figure 2. Top candidate SNPs that overlapped between gemcitabine and AraC.
49 SNPs were common (p<0.001) among 442 unique SNPs associated with gemcitabine and 504 associated with AraC.
Figure 3. Strategy applied in these studies.
A genome-wide SNP association study was performed with a drug-related phenotype, gemcitabine or AraC cytotoxicity (IC50). SNPs associated with cytotoxicity (p<10−3) were used to perform an association study with 54,000 expression probesets to identify SNPs that were associated with both gene expression (p<10−6) and with gemcitabine or AraC IC50. Finally, an association study was performed with gene expression and gemcitabine or AraC IC50 to identify SNPs that might be associated with drug IC50 values through an influence on gene expression.
Figure 4. Association among SNPs in IQGAP2 and TGM3 with expression and gemcitabine or AraC cytotoxicity.
(A) SNP rs3797418 (A/C) association with VAV3 expression, gemcitabine IC50 values and VAV3 expression association with gemcitabine IC50 in individual ethnic groups and all of the cell lines. (B) SNP rs6082527 (A/G) association with GPM6A gene expression, gemcitabine IC50 values and GPM6A gene expression association with gemcitabine IC50 values. (C). SNP rs3797418 (A/C) association with VAV3 expression, AraC IC50 values and VAV3 expression association with AraC IC50 in individual ethnic groups and all of the cell lines. Each dot represents one sample. Genotypes for each SNP are plotted against gene expression levels (upper panels) as well as gemcitabine or AraC IC50 values (middle panels). In the lower panel, correlations were determined between gene expression levels and gemcitabine or AraC IC50 values.
Figure 5. Functional characterization of candidate genes with siRNA knock down of VAV3 (A) and GPM6A (B).
Knock down of VAV3 and GPM6A expression in human SU86 and Hup-T3 pancreatic cancer cell lines as well as lymphoblastoid cell lines showed increased resistance to gemcitabine and AraC after siRNA knock down as determined by MTS assays. Error bars represent SEM values for 3 independent experiments. Quantitative RT-PCR was performed to assess VAV3 and GPM6A gene expression levels after knock down with specific siRNAs. Results are expressed as % of control. Error bars represent SEM values for 3 independent experiments. * = p<0.05.
Figure 6. Linkage disequilibrium within ∼300 Kb surrounding the rs3797418 and rs6082527 SNPs.
Red indicates combinations where D′ = 1 and linkage of disequilibrium (LOD) ≥2; light red, combinations where D′<1 and LOD ≥2; red squares, combinations where D′ = 1 and LOD <2; white squares, combinations with D′<1 and LOD <2. SNPs are arranged in order from 5′ to 3′ in each gene, as shown in the gene structure above each plot.
Similar articles
- Copy number variation and cytidine analogue cytotoxicity: a genome-wide association approach.
Kalari KR, Hebbring SJ, Chai HS, Li L, Kocher JP, Wang L, Weinshilboum RM. Kalari KR, et al. BMC Genomics. 2010 Jun 4;11:357. doi: 10.1186/1471-2164-11-357. BMC Genomics. 2010. PMID: 20525348 Free PMC article. - Gemcitabine and cytosine arabinoside cytotoxicity: association with lymphoblastoid cell expression.
Li L, Fridley B, Kalari K, Jenkins G, Batzler A, Safgren S, Hildebrandt M, Ames M, Schaid D, Wang L. Li L, et al. Cancer Res. 2008 Sep 1;68(17):7050-8. doi: 10.1158/0008-5472.CAN-08-0405. Cancer Res. 2008. PMID: 18757419 Free PMC article. - Discovery of genetic biomarkers contributing to variation in drug response of cytidine analogues using human lymphoblastoid cell lines.
Li L, Fridley BL, Kalari K, Niu N, Jenkins G, Batzler A, Abo RP, Schaid D, Wang L. Li L, et al. BMC Genomics. 2014 Feb 1;15:93. doi: 10.1186/1471-2164-15-93. BMC Genomics. 2014. PMID: 24483146 Free PMC article. - Structural basis for topoisomerase I inhibition by nucleoside analogs.
Gmeiner WH, Yu S, Pon RT, Pourquier P, Pommier Y. Gmeiner WH, et al. Nucleosides Nucleotides Nucleic Acids. 2003 May-Aug;22(5-8):653-8. doi: 10.1081/NCN-120022604. Nucleosides Nucleotides Nucleic Acids. 2003. PMID: 14565246 Review. - Pharmacogenetics and pharmacoepigenetics of gemcitabine.
Candelaria M, de la Cruz-Hernández E, Pérez-Cárdenas E, Trejo-Becerril C, Gutiérrez-Hernández O, Dueñas-González A. Candelaria M, et al. Med Oncol. 2010 Dec;27(4):1133-43. doi: 10.1007/s12032-009-9349-y. Epub 2009 Nov 10. Med Oncol. 2010. PMID: 19902390 Review.
Cited by
- Response and Toxicity to Cytarabine Therapy in Leukemia and Lymphoma: From Dose Puzzle to Pharmacogenomic Biomarkers.
Di Francia R, Crisci S, De Monaco A, Cafiero C, Re A, Iaccarino G, De Filippi R, Frigeri F, Corazzelli G, Micera A, Pinto A. Di Francia R, et al. Cancers (Basel). 2021 Feb 25;13(5):966. doi: 10.3390/cancers13050966. Cancers (Basel). 2021. PMID: 33669053 Free PMC article. Review. - USP15 promotes the apoptosis of degenerative nucleus pulposus cells by suppressing the PI3K/AKT signalling pathway.
Yu B, Shen B, Ba Z, Liu Z, Yuan J, Zhao W, Wu D. Yu B, et al. J Cell Mol Med. 2020 Dec;24(23):13813-13823. doi: 10.1111/jcmm.15971. Epub 2020 Nov 1. J Cell Mol Med. 2020. PMID: 33135363 Free PMC article. - Association of Genetic Variants at TRPC6 With Chemotherapy-Related Heart Failure.
Norton N, Crook JE, Wang L, Olson JE, Kachergus JM, Serie DJ, Necela BM, Borgman PG, Advani PP, Ray JC, Landolfo C, Di Florio DN, Hill AR, Bruno KA, Fairweather D. Norton N, et al. Front Cardiovasc Med. 2020 Aug 13;7:142. doi: 10.3389/fcvm.2020.00142. eCollection 2020. Front Cardiovasc Med. 2020. PMID: 32903434 Free PMC article. - Functional genomics based on germline genome-wide association studies of endocrine therapy for breast cancer.
Zayas J, Qin S, Yu J, Ingle JN, Wang L. Zayas J, et al. Pharmacogenomics. 2020 Jun;21(9):615-625. doi: 10.2217/pgs-2019-0191. Epub 2020 Jun 16. Pharmacogenomics. 2020. PMID: 32539536 Free PMC article. - A Humanized Yeast Phenomic Model of Deoxycytidine Kinase to Predict Genetic Buffering of Nucleoside Analog Cytotoxicity.
Santos SM, Icyuz M, Pound I, William D, Guo J, McKinney BA, Niederweis M, Rodgers J, Hartman JL IV. Santos SM, et al. Genes (Basel). 2019 Sep 30;10(10):770. doi: 10.3390/genes10100770. Genes (Basel). 2019. PMID: 31575041 Free PMC article.
References
- Sachidanandam R, Weissman D, Schmidt SC, Kakol JM, Stein LD, et al. A map of human genome sequence variation containing 1.42 million single nucleotide polymorphisms. Nature. 2001;409:928–933. - PubMed
- Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, et al. The structure of haplotype blocks in the human genome. Science. 2002;296:2225–2229. - PubMed
Publication types
MeSH terms
Substances
Grants and funding
- U19 GM061388/GM/NIGMS NIH HHS/United States
- GM61388/GM/NIGMS NIH HHS/United States
- R01 GM028157/GM/NIGMS NIH HHS/United States
- U01 GM061388/GM/NIGMS NIH HHS/United States
- P50 CA102701/CA/NCI NIH HHS/United States
- P20 CA102701/CA/NCI NIH HHS/United States
- R01 CA138461/CA/NCI NIH HHS/United States
- CA102701/CA/NCI NIH HHS/United States
- R01 CA132780/CA/NCI NIH HHS/United States
- K22 CA130828/CA/NCI NIH HHS/United States
- R01 GM28157/GM/NIGMS NIH HHS/United States
LinkOut - more resources
Full Text Sources
Other Literature Sources