Abstract
Nanzheng Chen1,*, Xiaomei Yang2,*, Wen Guo3, Jiangtao You1, Qifei Wu1, Guangjian Zhang1, Haijun Li1, Donghong Geng4, Tianbo Jin5, Junke Fu1, Yong Zhang1
1Department of Thoracic Surgery, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’an 710061, Shaanxi, China
2Hospital 521 of China’s Ordnance Industry Group, Xi’an 710065, Shaanxi, China
3Inner Mongolia Medical University, Hohhot 010050, Inner Mongolia, China
4School of Continuing Education of Xi’an Jiaotong University, Xi’an 710061, Shaanxi, China
5School of Life Sciences, Northwest University, Xi’an 710069, Shaanxi, China
*These authors have contributed equally to this work
Correspondence to:
Junke Fu, email: [email protected]
Yong Zhang, email: [email protected]
Keywords: association study, lung cancer, telomere-related gene, single nucleotide polymorphisms (SNPs)
Received: August 26, 2016 Accepted: October 14, 2016 Published: December 10, 2016
ABSTRACT
Single nucleotide polymorphisms (SNPs) in the telomere-associated gene ACYP2 are associated with increased lung cancer risk. We explored the correlation between ACYP2 SNPs and lung cancer susceptibility in the Chinese Han population. A total of 554 lung cancer patients and 603 healthy controls were included in this study. Thirteen SNPs in ACYP2 were selected. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated using unconditional logistic regression analysis. Multivariate logistic regression analysis was used assess the correlation between SNPs and lung cancer. We found that rs1682111 was associated with increased lung cancer risk in the recessive model (crude, OR=1.50, 95%CI: 1.04-2.16, p=0.029; adjusted for age, OR=1.55, 95%CI: 1.04-2.30, p=0.029), as was rs11896604 in the codominant model (crude, OR=0.65, 95%CI: 0.33-1.28, p=0.045; adjusted for age, OR=0.74, 95%CI: 0.36-1.53, p=0.049) and over-dominant model (crude, OR=1.30, 95%CI: 1.02-1.66, p=0.032; adjusted for age, OR=1.37, 95%CI: 1.05-1.78, p=0.020). Finally, rs843720 was associated with increased lung cancer risk in the recessive model (crude, OR=1.43, 95%CI: 1.02-2.02, p=0.040; adjusted for age, OR=1.48, 95%CI: 1.02-2.15, p=0.040). Thus three SNPs in ACYP2 (rs1682111, rs11896604 and rs843720) associate with lung cancer in the Chinese Han population.
INTRODUCTION
Lung cancer, the leading cause of cancer-related death in most countries, has been increasing rapidly in China for the last three decades [1]. The risk of lung cancer is strongly associated with exposure to environmental carcinogens, in particular to cigarette smoke [2]. However, only approximately 10% of smokers ever develop lung cancer. This highlights the potential role of hereditary components in determining the risk of lung cancer [3]. Through genome-wide association studies (GWAS), a number of genes involved in lung carcinogenesis have been identified in populations of European descent. These include CHRNA5, CHRNA3 and CHRNA4 on chromosome 15q25.1 [4–6], (TERT-CLPTM1L) on 5p15.33 [7, 8], and (BAT3-MSH5) on 6p21.33 [7].
It was also learned from GWAS that inherited single nucleotide polymorphisms (SNPs) in telomere-related genes are associated with various malignancies [9]. Human telomeres, composed of a tandem hexanucleotide repeat (TTAGGG), are many kilobases long in newborns but shorten by an average of 20-40 base-pairs annually [10]. SNPs in the telomere-related gene ACYP2 are reproducibly associated with inter-individual variation in leukocyte telomere length (LTL) [11]. In addition, recent Mendelian randomization studies indicate that LTL-associated SNPs related to longer telomere lengths are associated with increased risk of adult glioma, melanoma and lung cancer [12].
As far as we know, no studies have investigated the correlations between SNPs in ACYP2 and lung cancer susceptibility. We therefore performed a case-control association study to determine whether any of 13 SNPs in the ACYP2 gene are associated with lung cancer susceptibility in the Chinese Han population.
RESULTS
Participant characteristics
A total of 554 lung cancer patients (416 males, 138 females; median age at diagnosis 48.24 ± 13.05 years) and 603 healthy subjects (469 males, 134 females; median age 58.18 ± 10.53 years) were enrolled in our study. There is significant difference in age between patients in the case and control groups (p<0.001), but no difference in gender (Table 1).
Table 1: Characteristics of the cases and controls
Variables | Case (N=554) | Control (N=603) | Total | p |
|---|---|---|---|---|
Gender, No. (%) | 0.282a | |||
Male | 416 (75.1%) | 469 (77.8%) | 885 | |
Female | 138 (24.9%) | 134 (22.2%) | 272 | |
Mean age ±SD | 48.24 ± 13.05 | 58.18 ± 10.53 | <0.001b |
a p was calculated using the chi-square test.
b p was calculated using Student’s t-test.
Association between ACYP2 polymorphisms and lung cancer risk
The minor allele of each SNP was assumed to be a risk factor. The minor allele frequencies (MAF) are listed in Table 2. As a risk factor, the minor allele of each SNP was compared with the wild-type allele. All of the tested SNPs were consistent with the Hardy-Weinberg equilibrium for the control population of this study (p>0.05). Comparing the differences in the frequency distributions of the alleles between the cases and controls using the χ2 test, we found no correlation between the loci and lung cancer susceptibility.
Table 2: Basic information on the candidate SNPs
Rs number | Nucleotide position | Allele A/B | MAF frequency | HW p | Allele model | p | ||
|---|---|---|---|---|---|---|---|---|
case | control | OR | 95%CI | |||||
rs6713088 | 54345469 | G/C | 0.404 | 0.393 | 0.089 | 1.048 | 0.887-1.237 | 0.585 |
rs12621038 | 54391113 | T/C | 0.454 | 0.465 | 0.368 | 0.956 | 0.812-1.126 | 0.591 |
rs1682111 | 54427979 | A/T | 0.347 | 0.309 | 0.924 | 1.189 | 1.000-1.415 | 0.051 |
rs843752 | 54446587 | G/T | 0.234 | 0.260 | 0.089 | 0.870 | 0.720-1.052 | 0.151 |
rs10439478 | 54459450 | C/A | 0.430 | 0.433 | 0.245 | 0.987 | 0.837-1.164 | 0.880 |
rs843645 | 54474664 | G/T | 0.225 | 0.250 | 0.065 | 0.874 | 0.721-1.059 | 0.169 |
rs11125529 | 54475866 | A/C | 0.195 | 0.187 | 0.503 | 1.054 | 0.857-1.298 | 0.617 |
rs12615793 | 54475914 | A/G | 0.211 | 0.200 | 0.702 | 1.072 | 0.876-1.312 | 0.499 |
rs843711 | 54479117 | T/C | 0.436 | 0.448 | 0.510 | 0.953 | 0.809-1.123 | 0.567 |
rs11896604 | 54479199 | G/C | 0.214 | 0.200 | 0.800 | 1.087 | 0.889-1.330 | 0.415 |
rs843706 | 54480369 | A/C | 0.432 | 0.452 | 0.459 | 0.925 | 0.784-1.091 | 0.355 |
rs17045754 | 54496757 | C/G | 0.196 | 0.187 | 1.000 | 1.062 | 0.863-1.307 | 0.570 |
rs843720 | 54510660 | G/T | 0.380 | 0.343 | 0.527 | 1.171 | 0.988-1.388 | 0.068 |
SNPs: Single nucleotide polymorphisms; MAF: Minor allele frequency; HWE: Hardy-Weinberg equilibrium; OR: Odds ratio; CI: Confidence interval; A: Minor alleles; B: Major alleles.
After adjusting for age, the models were further analyzed using unconditional logistic regression analysis of SNPs associated with lung cancer (Table 3). We found that the AA genotype of rs1682111 was associated with increased lung cancer risk under the recessive model (crude, OR=1.50, 95%CI: 1.04-2.16, p=0.029; adjusted for age, OR=1.55, 95%CI: 1.04-2.30, p=0.029). We also noticed that rs11896604 was associated with lung cancer risk in codominant model (“GG”, crude, OR=0.65, 95%CI: 0.33-1.28, p=0.045; adjusted for age, OR=0.74, 95%CI: 0.36-1.53, p=0.049) and over-dominant model (“CG”, crude, OR=1.30, 95%CI: 1.02-1.66, p=0.032; adjusted for age, OR=1.37, 95%CI: 1.05-1.78, p=0.020). In addition, the GG genotype of rs843720 increased the risk of lung cancer under the recessive model (crude, OR=1.43, 95%CI: 1.02-2.02, p=0.040; adjusted for age, OR=1.48, 95%CI: 1.02-2.15, p=0.040).
Table 3: Single loci associated with lung cancer (adjusted by age)
Finally, we also performed a Wald test using unconditional multivariate regression analysis to assess the associations between SNP haplotypes and lung cancer risk. However, no positive results were observed.
DISCUSSION
In this study, three SNPs (rs1682111, rs11896604 and rs843720) were found to be associated with the lung cancer susceptibility in the Chinese Han population. Rs1682111 was associated with increased lung cancer risk in the recessive model. An association between this locus and other diseases has not been previously reported. Although this is the first report that rs11896604 associates with lung cancer susceptibility in the codominant over-dominant model, it was previously suggested that rs11896604 associates with a decreased risk of breast cancer [13]. In addition, rs843720 associated with an increased the risk of lung cancer in the recessive model.
ACYP2, located on chromosome 2p16.2, encodes a small cytosolic enzyme acylphosphatase that catalyzes the dephosphorylation of phospho-enzyme intermediates of various membrane pumps, particularly the Ca2+/Mg2+-ATPase from sarcoplasmic reticulum of skeletal muscle [14]. The physiological function of ACYP2 is not yet completely clear, however. Vos et al. reported that the rs1872328 variant of ACYP2 was associated with cisplatin-induced ototoxicity in patients with osteosarcoma who did not receive potentially ototoxic cranial irradiation [15]. It was further suggested that severe oxaliplatin-induced chronic peripheral neurotoxicity was potentially associated with ACYP2 rs843748 [16]. A recent genome-wide meta-analysis showed that ACYP2 rs11125529 affects telomere length and coronary heart disease in the Chinese Han population [17]. In addition, Yongjun et al found that ACYP2 rs12615793 and rs11896604 may significantly decrease the risk of high-altitude pulmonary edema [18]. However, the mechanism by which ACYP2 contributes to lung cancer will need to be explored in future studies.
Although the present study has several strengths, it also has intrinsic limitations. First, we present no data on the relation between ACYP2 SNPs and the most influential lung cancer risk: smoking. Second, the size of our sample is relatively small. Finally, other factors such as BMI, bacterial and viral infection, marriage, and economic status were not taken into account. Those issues will be addressed in future studies, along with further characterization of the molecular mechanism underlying the association of ACYP2 SNPs and lung cancer and prospective studies to validate their clinical utility.
In sum, we have identified several novel associations between three SNPs in ACYP2 (rs1682111, rs11896604 and rs843720) and lung cancer in the Chinese Han population. Our findings suggest ACYP2 may be a useful marker that informs clinical decisions, and may shed light on new candidate genes and new ideas about the mechanism governing the occurrence of lung cancer.
MATERIALS AND METHODS
Study participants
We recruited 554 patients with lung cancer and 603 healthy controls for this study. The patients were treated at the First Affiliated Hospital of Xi’an Jiao Tong University between January 2014 and August 2016. All demographic and related clinical data, including residential region, age, ethnicity, and education status, were collected through a face-to-face questionnaire and review of medical records. Patients recently diagnosed with primary lung cancer (confirmed by histopathological analysis) were included. Patients with other types of cancers or who underwent radiotherapy or chemotherapy were excluded. The controls underwent annual health evaluations in the checkup centers affiliated with our institution. All control patients were in good health and had no history of cancer, and they had no blood relatives with lung cancer going back three generations. This research was performed in accordance with the Chinese Department of Health and Human Services regulations for the protection of human research subjects. We obtained informed consent from all of the participants, and the study protocols were approved by the Institutional Review Board of Xi’an Jiao Tong University.
SNP selection and genotyping
Thirteen SNPs in ACYP2 that had a MAF >5% in the HapMap Asian population were selected for the association analysis [13, 14]. Venous blood samples (5 mL) were collected from each study participant during a laboratory examination. Blood samples from patients were collected prior to radiation or chemotherapy. DNA was extracted from whole blood samples using a Gold Mag-Mini Whole Blood Genomic DNA Purification Kit (GoldMag Ltd, Xian, China) and stored at –80°C after centrifugation. The DNA concentration was measured using spectrometry (DU530 UV/VIS spectrophotometer, Beckman Instruments, Fullerton, CA, USA). Sequenom MassARRAY Assay Design 3.0 software (Sequenom, Inc, San Diego, CA, USA) was used to design the multiplexed SNP Mass EXTEND assay, and genotyping was performed using a Sequenom MassARRAY RS1000 (Sequenom, Inc.) according to the manufacturer’s protocol. SequenomTyper 4.0 Software™ (Sequenom, Inc.) was used to manage and analyze the data. The primers corresponding to each SNP are shown in Table 4. Based on these results, the following 13 SNPs were selected: rs6713088, rs12621038, rs1682111, rs843752, rs10439478, rs843645, rs11125529, rs12615793, rs843711, rs11896604, rs843706, rs17045754 and rs843720. Basic information on the SNPs is shown in Table 2.
Table 4: Primers used
SNP_ID | 1st-PCRP | 2nd-PCRP | UEP_SEQ |
|---|---|---|---|
rs6713088 | ACGTTGGATGACACACACAGACTCCTTCAC | ACGTTGGATGGTCACCAAAACACGTAATG | gaggcCAGAATGGTCCACTAGAGA |
rs12621038 | ACGTTGGATGATTGTGCTAGGCACTTTAGG | ACGTTGGATGGGCATAAGTTTTATTGCCTC | ccATTGCCTCAGCTAGACT |
rs1682111 | ACGTTGGATGGAATTGCTGGGTTATTTGGC | ACGTTGGATGGCCAGTGGGAATGCAAAATG | tgtcATGCAAAATGAAACAGACACTT |
rs843752 | ACGTTGGATGTCCTCTTTTCAGAAACCTGC | ACGTTGGATGGAGACAACATAATGGAGGTC | cGAGTTTGGGTTTGAGGT |
rs10439478 | ACGTTGGATGTAGCACAAGACCTACACTGG | ACGTTGGATGCTACACTCTCCAGAGGAATG | TTGCTGTTTTCCCAGAA |
rs843645 | ACGTTGGATGGAAATCTGAATACCACCTAC | ACGTTGGATGACAGTGCCTTTAGCAAGGTG | TCATAGGCACTACTGTATC |
rs11125529 | ACGTTGGATGCCGAAGAAAAGAAGATGAC | ACGTTGGATGGAGCTTAGTTGTTTACAGATG | AGAAAAGAAGATGACTAAAACAT |
rs12615793 | ACGTTGGATGATCTTGGCCCTTGAAGAA | ACGTTGGATGTTTGAGCTTAGTTGTTTAC | AAATTGAGTGACAAATATAAACTAC |
rs843711 | ACGTTGGATGGACAAAGGACCTTACAACTC | ACGTTGGATGTGCCTTGTGGGAATTAGAGC | gggaTCAGGGAACCAGTGCAAA |
rs11896604 | ACGTTGGATGAAGTCAGAATAGTGCTTAC | ACGTTGGATGTGTCTCTGACCTAGCATGTA | GTTAAGCTTGCAAGGAG |
rs843706 | ACGTTGGATGTGAAAGCCATAAATATTTTG | ACGTTGGATGTGAATAACTTGGTCTTATC | cACTTGGTCTTATCTGATGC |
rs17045754 | ACGTTGGATGCTGTAAAAGTTCTGGCATGG | ACGTTGGATGGAAATCAGGGATATTAGTGC | caggTATTCAGCTTCCTAGAGTTA |
rs843720 | ACGTTGGATGCTTCACAACACTCCTGTAAG | ACGTTGGATGAGTCAGAGCTAGACCTCTGG | ccccAATCTGTCTCAGGGTCTT |
Statistical analysis
We used Chi-squared tests (categorical variables) and Student’s t-tests (continuous variables) to assess the differences in the demographic characteristics between the cases and controls. The Hardy-Weinberg equilibrium of each SNP was assessed in order to compare the expected frequencies of the genotypes in the control groups. All minor alleles were regarded as risk alleles for lung cancer susceptibility. To evaluate associations between each SNP and lung cancer risk in the five models (codominant, dominant, recessive, over-dominant and log-additive), ORs and 95% CIs were calculated using unconditional logistic regression analysis after adjusting for age, and gender. Linkage disequilibrium analysis and SNP haplotypes were analyzed using the Haploview software package (version 4.2) and the SHEsi software platform (http://www.nhgg.org/analysis/). All statistical analyses were performed using the SPSS version 17.0 statistical package (SPSS, Chicago, IL, USA) and Microsoft Excel. All statistical tests were two-sided, and values of p<0.05 was considered significant.
ACKNOWLEDGMENTS
This work was supported by the Fundamental Research Funds for the Central University of China [No.08143004] and the Funds for Science and Technology Project of Shaanxi Province of China [No.2014K11-02-03-07]. The authors are also grateful to the patients and control individuals for their participation in the study. We thank the clinicians and hospital staff who contributed to sample and data collection for this study.
CONFLICTS OF INTEREST
The authors declare no conflicts of interest.
REFERENCES
1. Hu Z, Wu C, Shi Y, Guo H, Zhao X, Yin Z, Yang L, Dai J, Hu L, Tan W, Li Z, Deng Q, Wang J, Wu W, Jin G, Jiang Y, et al. A genome-wide association study identifies two new lung cancer susceptibility loci at 13q12.12 and 22q12.2 in Han Chinese. Nature genetics. 2011; 43:792-796.
2. Chaturvedi AK, Caporaso NE, Katki HA, Wong HL, Chatterjee N, Pine SR, Chanock SJ, Goedert JJ, Engels EA. C-reactive protein and risk of lung cancer. Journal of clinical oncology: official journal of the American Society of Clinical Oncology. 2010; 28:2719-2726.
3. Spinola M, Meyer P, Kammerer S, Falvella FS, Boettger MB, Hoyal CR, Pignatiello C, Fischer R, Roth RB, Pastorino U, Haeussinger K, Nelson MR, Dierkesmann R, Dragani TA, Braun A. Association of the PDCD5 locus with lung cancer risk and prognosis in smokers. Journal of clinical oncology: official journal of the American Society of Clinical Oncology. 2006; 24:1672-1678.
4. Hung RJ, McKay JD, Gaborieau V, Boffetta P, Hashibe M, Zaridze D, Mukeria A, Szeszenia-Dabrowska N, Lissowska J, Rudnai P, Fabianova E, Mates D, Bencko V, Foretova L, Janout V, Chen C, et al. A susceptibility locus for lung cancer maps to nicotinic acetylcholine receptor subunit genes on 15q25. Nature. 2008; 452:633-637.
5. Amos CI, Wu X, Broderick P, Gorlov IP, Gu J, Eisen T, Dong Q, Zhang Q, Gu X, Vijayakrishnan J, Sullivan K, Matakidou A, Wang Y, Mills G, Doheny K, Tsai YY, et al. Genome-wide association scan of tag SNPs identifies a susceptibility locus for lung cancer at 15q25.1. Nature genetics. 2008; 40:616-622.
6. Thorgeirsson TE, Geller F, Sulem P, Rafnar T, Wiste A, Magnusson KP, Manolescu A, Thorleifsson G, Stefansson H, Ingason A, Stacey SN, Bergthorsson JT, Thorlacius S, Gudmundsson J, Jonsson T, Jakobsdottir M, et al. A variant associated with nicotine dependence, lung cancer and peripheral arterial disease. Nature. 2008; 452:638-642.
7. Wang Y, Broderick P, Webb E, Wu X, Vijayakrishnan J, Matakidou A, Qureshi M, Dong Q, Gu X, Chen WV, Spitz MR, Eisen T, Amos CI, Houlston RS. Common 5p15.33 and 6p21.33 variants influence lung cancer risk. Nature genetics. 2008; 40:1407-1409.
8. McKay JD, Hung RJ, Gaborieau V, Boffetta P, Chabrier A, Byrnes G, Zaridze D, Mukeria A, Szeszenia-Dabrowska N, Lissowska J, Rudnai P, Fabianova E, Mates D, Bencko V, Foretova L, Janout V, et al. Lung cancer susceptibility locus at 5p15.33. Nature genetics. 2008; 40:1404-1406.
9. Walsh KM, Whitehead TP, de Smith AJ, Smirnov IV, Park M, Endicott AA, Francis SS, Codd V, Group ECT, Samani NJ, Metayer C, Wiemels JL. Common genetic variants associated with telomere length confer risk for neuroblastoma and other childhood cancers. Carcinogenesis. 2016; 37:576-582.
10. Daniali L, Benetos A, Susser E, Kark JD, Labat C, Kimura M, Desai K, Granick M, Aviv A. Telomeres shorten at equivalent rates in somatic tissues of adults. Nature communications. 2013; 4:1597.
11. Codd V, Nelson CP, Albrecht E, Mangino M, Deelen J, Buxton JL, Hottenga JJ, Fischer K, Esko T, Surakka I, Broer L, Nyholt DR, Mateo Leach I, Salo P, Hagg S, Matthews MK, et al. Identification of seven loci affecting mean telomere length and their association with disease. Nature genetics. 2013; 45:422-427, 427e421-422.
12. Zhang C, Doherty JA, Burgess S, Hung RJ, Lindstrom S, Kraft P, Gong J, Amos CI, Sellers TA, Monteiro AN, Chenevix-Trench G, Bickeboller H, Risch A, Brennan P, McKay JD, Houlston RS, et al. Genetic determinants of telomere length and risk of common cancers: a Mendelian randomization study. Human molecular genetics. 2015; 24:5356-5366.
13. Liu M, Li B, Guo W, Zhang X, Chen Z, Li J, Yan M, Chen C, Jin T. Association between single nucleotide polymorphisms in the TSPYL6 gene and breast cancer susceptibility in the Han Chinese population. Oncotarget. 2016; 7:54771-54781. doi: 10.18632/oncotarget.10754.
14. Degl'Innocenti D, Marzocchini R, Rosati F, Cellini E, Raugei G, Ramponi G. Acylphosphatase expression during macrophage differentiation and activation of U-937 cell line. Biochimie. 1999; 81:1031-1035.
15. Vos HI, Guchelaar HJ, Gelderblom H, de Bont ES, Kremer LC, Naber AM, Hakobjan MH, van der Graaf WT, Coenen MJ, te Loo DM. Replication of a genetic variant in ACYP2 associated with cisplatin-induced hearing loss in patients with osteosarcoma. Pharmacogenetics and genomics. 2016; 26:243-247.
16. Oguri T, Mitsuma A, Inada-Inoue M, Morita S, Shibata T, Shimokata T, Sugishita M, Nakayama G, Uehara K, Hasegawa Y, Ando Y. Genetic polymorphisms associated with oxaliplatin-induced peripheral neurotoxicity in Japanese patients with colorectal cancer. International journal of clinical pharmacology and therapeutics. 2013; 51:475-481.
17. Mathieson I, McVean G. Differential confounding of rare and common variants in spatially structured populations. Nature genetics. 2012; 44:243-246.
18. He Y, Chen J, Zhang X, Li X, Du J, He X, Zhang Z, Zhang Y, Kang L, Jin T, Yuan D. Telomere length related gene ACYP2 polymorphism is associated with the risk of HAPE in Chinese Han population. The journal of gene medicine. 2016.
