Novel SNP improves differential survivability and mortality in non-small cell lung cancer patients
- Tzia Liang Mah1Email author,
- Xin Ning Adeline Yap†1,
- Vachiranee Limviphuvadh†3,
- Nanpu Li†1,
- Srinath Sridharan1,
- Vellaisemy Kuralmani1,
- Mengling Feng1,
- Natalia Liem2,
- Sharmila Adhikari3,
- Wei Peng Yong2,
- Ross A Soo2, 9,
- Sebastian Maurer-Stroh3, 4,
- Frank Eisenhaber3, 5, 6 and
- Joo Chuan Tong7, 8Email author
© Mah et al.; licensee BioMed Central Ltd. 2014
Published: 8 December 2014
Non-small cell lung cancer (NSCLC) is a major cause of cancer-related death worldwide due to poor patient prognosis and clinical outcome. Here, we studied the genetic variations underlying NSCLC pathogenesis based on their association to patient outcome after gemcitabine therapy.
Bioinformatics analysis was used to investigate possible effects of POLA2 G583R (POLA2+1747 GG/GA, dbSNP ID: rs487989) in terms of protein function. Using biostatistics, POLA2+1747 GG/GA (rs487989, POLA2 G583R) was identified as strongly associated with mortality rate and survival time among NSCLC patients. It was also shown that POLA2+1747 GG/GA is functionally significant for protein localization via green fluorescent protein (GFP)-tagging and confocal laser scanning microscopy analysis. The single nucleotide polymorphism (SNP) causes DNA polymerase alpha subunit B to localize in the cytoplasm instead of the nucleus. This inhibits DNA replication in cancer cells and confers a protective effect in individuals with this SNP.
The results suggest that POLA2+1747 GG/GA may be used as a prognostic biomarker of patient outcome in NSCLC pathogenesis.
Non-small cell lung cancer (NSCLC) is a leading cause of cancer mortality worldwide with over one million deaths annually . It accounts for 75% of lung cancer cases and consists of three major subtypes: adenocarcinoma, large-cell carcinoma, and squamous-cell carcinoma . Recent introduction of targeted therapy and increasing numbers of available chemotherapeutic regimens, such as platinums, taxanes and gemcitabine, do not effectively cure NSCLC patients, with varied response towards treatment and occurrence of drug toxicity [3, 4]. In addition, prognosis remains dismal in NSCLC patients albeit careful evaluation of clinico-pathological factors that determine patient response to therapy, such as tumor, nodes and metastasis (TNM) staging, performance status, gender and weight loss. The long-term survival rate is low with only 14% of patients surviving five years after diagnosis  and the risk for relapse is high.
Gemcitabine is a third generation chemotherapeutic agent that has shown activity in NSCLC. Preclinical studies have shown that the compound is a potent radiosensitizer, with response in stage III NSCLC . Gemcitabine can be administered as a single agent, or in platinum and non-platinum combination. The agent can also be combined with the chemotherapy drug pemetrexed, as well as the vascular endothelial growth factor (VEGF) inhibitor, for adenocarcinoma NSCLC. Due to its significant benefit and advantageous toxicity profile, gemcitabine has since evolved to become one of the most commonly used agents for lung cancer chemotherapy.
In recent years, much effort has been expended to identify genetic determinants in patient outcomes, so as to improve clinical treatment decisions and for the design of therapeutic agents. The epidermal growth factor receptor (EGFR) mutations, for instance, are common in patients with NSCLC , and are known to confer survival benefit and better clinical outcome when treated with EGFR tyrosine kinase inhibitors (TKIs) [8, 9]. To date, no known genetic variants have been reported, that could help determine the dose and clinical outcomes in NSCLC patients receiving gemcitabine chemotherapy.
Here, we studied the polymorphism of genes involved in gemcitabine transport, metabolism and activity, based on their association to patient outcome after gemcitabine therapy. We showed for the first time that the single nucleotide polymorphism (SNP) POLA2+1747 GG/GA (rs487989) is a key determinant of mortality and survival outcome in gemcitabine-treated NSCLC patients. The POLA2 gene encodes DNA polymerase alpha subunit B in humans, which is involved in the initiation of chromosomal DNA replication [10–12]. The SNP causes DNA polymerase alpha subunit B to localize in the cytoplasm instead of the nucleus. This inhibits DNA replication in cancer cells and confers a protective effect in individuals with this SNP. The results suggest that POLA2+1747 GG/GA (rs487989) may be used as a prognostic biomarker of patient outcome in NSCLC pathogenesis.
Results and discussion
Association of genotypes and the mortality of NSCLC patients after gemcitabine therapy
Association between the 21 SNP genotypes and mortality of 43 NSCLC patients receiving gemcitabine, based on the corresponding P value.
Significance Level (P value)
POLA2+1747=GG/GAimposes differential effects on mortality and survival time of NSCLC patients after gemcitabine therapy
Using conditional probability test, we studied the effects of this POLA2 variant on mortality of patients after gemcitabine therapy. We found that the probability of death (P value = 0.0128) for patients with wild-type GG genotype is significantly higher (89.19%) than those with GA variant (50%).
Association between POLA2+1747 GG/GA interaction pairs and the overall survival time of 43 NSCLC patients receiving gemcitabine, based on the corresponding p-value.
Genotype Interaction Pairs
POLA2+1747 GA together with SLC28A2+65 CC are observed to be associated with increased median survival time (Figure 1A). For POLA2+1747 GG/GA together with SLC28A2+65 CC, the median overall survival time of patients is 7.39 months and 13.18 months in patients with GG and GA genotypes, respectively (P value = 0.0004). Likewise, we also observed that the POLA2+1747 GA together with SLC28A2+225 CC are associated with increased median survival time (Figure 1B). For POLA2+1747 GG/GA together with SLC28A2+225=CC, the median overall survival time of patients with the wild-type GG genotype was 7.39 months and for GA variant is 13.17 months (P value = 0.0010) (Figure 1B).This result indicates that the non-synonymous POLA2+1747 GA SNP is an important interactor associated with increased survival time.
Computational prediction of functional effects of POLA2 G583R
Subcellular localization of wild type POLA2 and mutant POLA2 G583R proteins
In a recent study on genes involved in gemcitabine pharmacology in ethnic Asian populations, we reported on the use of a statistical approach to examine associations between genotypes and the outcome of NSCLC patients including response rate, time to progression, gemcitabine toxicity and overall survival. We have now extended the study to another aspect of NSCLC patient outcome that was not examined previously, i.e. mortality, and also shown here that the POLA2+174 GG/GA (rs487989) is strongly associated with mortality rate and survival time among NSCLC patients treated with gemcitabine. We have previously shown that this particular SNP by itself did not have a significant effect on survival time . Now, we found that its interaction with SLC28A2+65 CC and SLC28A2+225 CC led to an increase in the overall survival times of NSCLC patients.
This POLA2+1747 variant (rs487989) is not only present in the European and African populations , but is also prevalent in the Asian population among Chinese, Indians and Malays . This SNP encodes for a glycine to arginine amino acid change (G583R), where G is an ancestral allele, resulting inside chain polarity and charge reversal. Here, we showed, through biostatistics, that individuals with the ancestral allele G for POLA2 tend to have lower survival rates in NSCLC pathogenesis, compared to individuals with GA polymorphism. To unravel possible molecular mechanisms of functional effects of this mutation, we utilized multiple computational approaches based on evolutionary conservation, structural modelling and molecular dynamics simulations. Given its location in a surface loop of the structure and causing flexible rearrangements of this surface area, we hypothesized that it could disrupt protein interactions which may be important for subcellular localization. Indeed, we experimentally showed that this point mutation is functionally significant, leading to a change in localization that is likely to affect regulatory activity and induces better survival in NSCLC patients treated with gemcitabine. The wild type POLA2 that is known to facilitate nuclear DNA replication is predominantly found in the nucleus, whereas the mutant POLA2 G583R protein  that is strongly associated with better survival in NSCLC patients is mainly localized in the cytoplasm. DNA polymerase alpha subunit B is required for cell viability . By localizing in the cytoplasm, nuclear DNA polymerase alpha activity is inhibited. This confers a protective effect in NSCLC patients who possess the POLA2+1747 GG/GA SNP genotype, as the tumour DNA could not replicate. This inhibits tumour cell proliferation, and ultimately results in tumour cell death.
In summary, we established that the POLA2+1747 GG/GA (rs487989) is a genetic determinant of clinical outcomes in NSCLC patients receiving gemcitabine treatment. EGFR mutations are used for profiling NSCLC patients treated with EGFR tyrosine kinase inhibitors, and similarly, the findings in this article can become a stepping stone for the discovery of new options for gemcitabine-based therapy. Due to the lack of genetic variants that could help determine the dose and clinical outcomes in NSCLC patients receiving gemcitabine chemotherapy, such biomarkers would be useful for doctors in treating patients more efficiently to achieve satisfactory clinical outcome and better survival.
Materials and methods
Summary of the study population used in this work.
Total number of patients
Median time to progression
Range of age
Selection of gene variant loci
Here, we assessed 21 non-synonymous SNPs in 9 genes involved in gemcitabine transport, metabolism and activity . The SNPs are found in the respective gene variant loci, namely, CDA+79 (rs2072671), CDA+208 (rs60369023), CDA+435 (rs1048977), DCK+3122 (rs3775289), DCK+36791 (rs1803484), DCTD+315 (rs4742), POLA2+1747 (rs487989), RRM1-756(rs11030918), RRM1-269 (rs12806698), S28A1+419 (rs2277576), S28A1+565 (rs2290272), S28A1+709 (rs8187758), S28A1+1368 (rs2242048), S28A1+1528 (rs2242047), S28A1+1561 (rs2242046), S28A2+65 (rs61637002), S28A2+225 (rs1060896), S28A3+338 (rs10868138), TYMS-100 (rs34743033), TYMS-58 (rs2853542), TYMS+15705 (rs34489327).
Fisher's exact probability test was used to assess the relationship between each of the 21 SNPs and the mortality of 43 NSCLC patients based on the p-values between genotypes. Conditional probability of death given a geno-type of a SNP was used to characterize the differential effects on mortality. Chi-squared test was employed to confirm the significance (P value) of the difference between genotypes. Differences were considered statistically significant when the P value was less than 0.05. All statistical tests were two-sided. Kaplan-Meier method and log-rank test were used to compare overall survival time for interaction pairs. SPSS software version 14.0 (SPSS Inc., Chicago, IL) was used.
In order to investigate possible effects of POLA2 G583R (POLA2+1747 GG/GA, dbSNP ID: rs487989) in terms of protein function, we analysed the mutation with PolyPhen-2 version 2.2.2 using the rsid (rs487989) of POLA2 G583R as input, SNAP using the amino acid sequence of human POLA2 (RefSeq ID: NP_002680) as input (Bromberg et al., 2008) and SIFT using a curated alignment of orthologous sequences (Ng and Henikoff, 2001). Orthologous sequences of human POLA2 (RefSeq ID: NP_002680) were retrieved with the orthologue search in ANNOTATOR (Ooi et al, 2009). A multiple alignment was created using MAFFT with the L-INS-I algorithm (Katoh and Toh, 2008). After deleting sequence that had large gaps in Jalview (Waterhouse et al, 2009), we selected diverse organisms and 24 remaining sequences were compared.
Homology modelling of human POLA2 was performed with Modeller (Eswar et al, 2008) using the crystal structure of the carboxyl-terminal domain of yeast DNA polymerase alpha in complex with its B subunit (PDB:3FLO) (Klinge et al, 2009) as template. Then, the alignment containing 24 diverse orthologous sequences (mentioned above) was used to calculate the conservation on individual positions using the evolutionary trace algorithm (Lichtarge et al, 2003) and the level of residue conservation was mapped to its corresponding position in the model and visualized with YASARA (Krieger et al, 2004). Moreover, we used FoldX (Schymkowitz et al, 2005) with prior energy minimization using the RepairPDB function and 5 repetitions of the mutation stability change calculations to predict the SNP effect on protein structure stability. Lastly, we performed 5 wildtype and 5 mutant MD simulations over 10 ns in explicit water using the AMBER03 force field in YASARA (Krieger et al, 2004) following standard protocols to understand the effect of the SNP on protein structure flexibility.
Isolation of total RNA from HEK 293 cell culture
HEK 293 cells were lysed directly in a 10 cm culture dish using TRIZOL® Reagent (Invitrogen, Carls-bad, CA, USA). Total RNA was isolated and used for further experiments only if the RNA was found intact by running on 1% denaturing agarose gel.
Reverse transcription of total RNA and POLA2gene amplification by PCR
Using SuperScript™ III One-Step RT-PCR System with Platinum® Taq High Fidelity (Invitrogen), total RNA was reverse transcribed into complimentary DNA (cDNA), followed by amplification of the wild type POLA2 using forward primer with Kpn1 restriction site; 5'-AAGGTACCATGTCCGCATCCGCCCAGCA-3' and reverse primer with BamH1 restriction site; 5'-AAGGATCCGATCCTGACGACCTGCACAGCA-3'. Amplicon size was verified by running the PCR product and GeneRuler™ 100 base pairs DNA ladder on 0.8% agarose gel.
Cloning of POLA2amplicon and sequence verification
The POLA2 amplicon was cloned into pGEM®-T easy vector (Promega, USA), and subsequently transformed into Escherichia coli DH5α bacteria. The plasmid DNA was extracted and purified using QIAprep Spin Miniprep Kit (QIAGEN, Germany). Next, the concentration and purity of plasmid DNA was measured using NanoDrop (Thermo Fisher Scientific, USA). The resultant plasmid was digested with EcoRI and ran on 0.8% agarose gel to identify recombinant clones. Integrity of the wild type POLA2 constructs was verified by sequencing. Point mutation was introduced into POLA2 using the XL QuikChange Site-Directed Mutagenesis Kit, and POLA2+1747 GG/GA were verified by sequencing. Wild type POLA2 and POLA2+1747 GG/GA were then cloned into pEGFP-N3 (Clontech, USA) at Kpn1 and BamHI restriction sites.
Cell culture and transfection
HEK 293 cells were grown in Dulbecco's modification of Eagle's medium (DMEM) supplemented with 10% fetal bovine serum (FBS) and 1% PSG (penicillin/streptomycin/glutamine) on 6-well plate and maintained at 37°C and 5% CO2. HEK 293 cells were transiently transfected using Lipofectamine 2000 (Invitrogen), following the manufacturer's protocol. After 5 hours of transfection with POLA2-GFP constructs in Opti-MEM® I Reduced Serum Medium (Cat. No. 31985-062), 30% FBS was added into each well, and incubated at 37°C and 5% CO2 overnight. After 17-18 hours of transfection with POLA2-GFP constructs using Lipofectamine 2000 (Invitrogen), cells were lysed in 1× cell lysis buffer from Matchmaker™Chemiluminescent Co-IP System (Clontech) with 25× complete EDTA free protease inhibitor (Roche, Germany) and 100× Phenylmethylsulfonyl fluoride (Sigma-Aldrich®, USA).
Equivalent volumes (20 µl) of cell lysates were loaded onto 8% SDS-PAGE gels to resolve proteins. Then, proteins were transferred onto PVDF membrane and blocked using 5% non-fat dry milk for 1 hour to reduce non-specific binding and incubated overnight at 4°C with mouse anti-GFP (Roche) at 1:2500 dilution in 0.5% non-fat dry milk followed by Goat anti mouse IgG-HRP (Santa Cruz Biotechnology, US) at 1:10000 dilution in 0.5% non-fat dry milk for 1 hour at room temperature. Immunoblots were developed using Amersham™ ECL™ Prime Western Blotting Detection Reagent (GE healthcare, Sweden), following the manufacturer's protocol.
Confocal laser scanning microscopic analysis
After 17-18 hours of transfection with the POLA2-GFP constructs using Lipofectamine 2000 (Invitrogen), HEK293 cells grown on glass cover slips were fixed with 2% paraformaldehyde in PBS at room temperature. Slides were blocked at room temperature for 1 hour with 5% BSA in 0.1% Triton/PBS and then immunostained with mouse anti-GFP (Roche) at 1:100 dilution followed by the Alexa Fluor 488 donkey anti-mouseIgG (Invitrogen, Molecular Probes) (1:2000 dilution) at room temperature for an hour. Images were captured with Zeiss LSM Meta confocal inverted mi-croscope with a magnification of 63×.
This work, including funding for open access publication charges, was supported by the Agency for Science, Technology and Research (A*STAR) Joint Council Office (JCO) Grant JCOAG04_FG03_2009.
This article has been published as part of BMC Genomics Volume 15 Supplement 9, 2014: Thirteenth International Conference on Bioinformatics (InCoB2014): Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/15/S9.
- Jemal A, Siegel R, Ward E, Hao Y, Xu J, Murray T, Thun MJ: Cancer statistics, 2008. CA Cancer J Clin. 2008, 58: 71-96. 10.3322/CA.2007.0010.PubMedView ArticleGoogle Scholar
- Nacht M, Dracheva T, Gao Y, Fujii T, Chen Y, Player A, Akmaev V, Cook B, Dufault M, Zhang M, et al: Molecular characteristics of non-small cell lung cancer. Proc Natl Acad Sci USA. 2001, 98: 15203-15208. 10.1073/pnas.261414598.PubMedPubMed CentralView ArticleGoogle Scholar
- Shepherd FA: Chemotherapy for non-small cell lung cancer: have we reached a new plateau?. Semin Oncol. 1999, 26: 3-11.PubMedGoogle Scholar
- Tiseo M, Franciosi V, Grossi F, Ardizzoni A: Adjuvant chemotherapy for non-small cell lung cancer: ready for clinical practice?. Eur J Cancer. 2006, 42: 8-16. 10.1016/j.ejca.2005.08.031.PubMedView ArticleGoogle Scholar
- Spira A, Ettinger DS: Multidisciplinary management of lung cancer. N Engl J Med. 2004, 350: 379-392. 10.1056/NEJMra035536.PubMedView ArticleGoogle Scholar
- Le Chevalier T, Scagliotti G, Natale R, Danson S, Rosell R, Stahel R, Thomas P, Rudd RM, Vansteenkiste J, Thatcher N, et al: Efficacy of gemcitabine plus platinum chemotherapy compared with other platinum containing regimens in advanced non-small-cell lung cancer: a meta-analysis of survival outcomes. Lung Cancer. 2005, 47: 69-80. 10.1016/j.lungcan.2004.10.014.PubMedView ArticleGoogle Scholar
- Janne PA, Engelman JA, Johnson BE: Epidermal growth factor receptor mutations in non-small-cell lung cancer: implications for treatment and tumor biology. J Clin Oncol. 2005, 23: 3227-3234. 10.1200/JCO.2005.09.985.PubMedView ArticleGoogle Scholar
- Tsao MS, Sakurada A, Cutz JC, Zhu CQ, Kamel-Reid S, Squire J, Lorimer I, Zhang T, Liu N, Daneshmand M, et al: Erlotinib in lung cancer - molecular and clinical predictors of outcome. N Engl J Med. 2005, 353: 133-144. 10.1056/NEJMoa050736.PubMedView ArticleGoogle Scholar
- Zhu CQ, da Cunha Santos G, Ding K, Sakurada A, Cutz JC, Liu N, Zhang T, Marrano P, Whitehead M, Squire JA, et al: Role of KRAS and EGFR as biomarkers of response to erlotinib in National Cancer Institute of Canada Clinical Trials Group Study BR.21. J Clin Oncol. 2008, 26: 4268-4275. 10.1200/JCO.2007.14.8924.PubMedView ArticleGoogle Scholar
- Wang TS: Eukaryotic DNA polymerases. Annu Rev Biochem. 1991, 60: 513-552. 10.1146/annurev.bi.60.070191.002501.PubMedView ArticleGoogle Scholar
- Sugino A: Yeast DNA polymerases and their role at the replication fork. Trends Biochem Sci. 1995, 20: 319-323. 10.1016/S0968-0004(00)89059-3.PubMedView ArticleGoogle Scholar
- Foiani M, Lucchini G, Plevani P: The DNA polymerase alpha-primase complex couples DNA replication, cell-cycle progression and DNA-damage response. Trends Biochem Sci. 1997, 22: 424-427. 10.1016/S0968-0004(97)01109-2.PubMedView ArticleGoogle Scholar
- Mizuno T, Yamagishi K, Miyazawa H, Hanaoka F: Molecular architecture of the mouse DNA polymerase alpha-primase complex. Mol Cell Biol. 1999, 19: 7886-7896.PubMedPubMed CentralView ArticleGoogle Scholar
- Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, Kondrashov AS, Sunyaev SR: A method and server for predicting damaging missense mutations. Nature methods. 2010, 7: 248-249. 10.1038/nmeth0410-248.PubMedPubMed CentralView ArticleGoogle Scholar
- Soo RA, Wang LZ, Ng SS, Chong PY, Yong WP, Lee SC, Liu JJ, Choo TB, Tham LS, Lee HS, et al: Distribution of gemcitabine pathway genotypes in ethnic Asians and their association with outcome in non-small cell lung cancer patients. Lung Cancer. 2009, 63: 121-127. 10.1016/j.lungcan.2008.04.010.PubMedView ArticleGoogle Scholar
- Fukunaga AK, Marsh S, Murry DJ, Hurley TD, McLeod HL: Identification and analysis of single-nucleotide polymorphisms in the gemcitabine pharmacologic pathway. Pharmacogenomics J. 2004, 4: 307-314. 10.1038/sj.tpj.6500259.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.