Abstract
Yang Zhang1,2,*, Xiaoai Zhang1,3,*, Hongxing Zhang1, Yun Zhai1, Zhifu Wang1,2, Peiyao Li1, Lixia Yu1, Xia Xia1, Ying Zhang1, Yixin Zeng4, Fuchu He1, Gangqiao Zhou1
1State Key Laboratory of Proteomics, Beijing Proteome Research Center, Beijing Institute of Radiation Medicine, Beijing, China
2Institute of Basic Medical Sciences, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
3State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Microbiology and Epidemiology, Beijing, China
4Department of Experimental Research, Sun Yat-sen University Cancer Center, Guangzhou, China
*These authors contributed equally to this work
Correspondence to:
Fuchu He, e-mail: [email protected]
Gangqiao Zhou, e-mail: [email protected]
Keywords: TERT, CLPTM1L, polymorphism, nasopharyngeal carcinoma
Received: April 28, 2015 Accepted: November 14, 2015 Published: November 26, 2015
ABSTRACT
Associations between single nucleotide polymorphisms (SNPs) at 5p15 (TERT-CLPTM1L) and multiple cancer types have been reported. We examined whether polymorphisms in the TERT-CLPTM1L locus were related to the risk of developing nasopharyngeal carcinoma (NPC) among Chinese populations. In the first stage, 26 tag SNPs were genotyped in a Guangxi population (855 patients and 1036 controls). In the second stage, the SNPs, which showed significant association, were further genotyped in a Guangdong population (997 patients and 972 controls). Functional analyses were conducted to verify the biological relevance of the associated polymorphism. In the 1st stage, four SNPs (rs2736098, rs2735845, rs402710, and rs401681) were significantly associated with the risk of developing NPC. After the 2nd stage validation, rs2735845 and rs401681 were independently associated with the risk of developing NPC in the additive model (rs2735845, OR = 1.19, 95% CI = 1.04–1.37, P = 0.011; rs401681, OR = 0.85, 95% CI = 0.74–0.99, P = 0.034). Furthermore, we observed higher CLPTM1L messenger RNA levels in fetal mesenchymal stem cells from the rs2735845 G allele carriers compared with that from non-carriers. In addition, using an immunohistochemistry assay, we observed higher TERT and CLPTM1L levels in NPC tissues compared with that in non-cancerous nasopharyngeal tissues. Our findings suggest that polymorphisms in the TERT-CLPTM1L locus may play a role in mediating the susceptibility to NPC in Chinese populations.
INTRODUCTION
Nasopharyngeal carcinoma (NPC) is an epithelial malignancy of the head and neck. The incidence of NPC is higher in Southeast Asia and Africa but lower among Caucasians in North America and Europe. In southern China, the incidence rate is approximately 15–50 per 100,000 person-years, which is 100-fold higher than that in the Western world [1]. NPC is a complex disease caused by an interaction among Epstein-Barr virus infection, environmental risk factors and host genes in a multi-step process of carcinogenesis [2]. The identification of susceptibility genes contributing to NPC may help clarify the pathogenesis of carcinogenesis and improve the prevention and treatment of this malignancy.
The human 5p15.33 locus contains two well-known genes, the telomerase reverse transcriptase (TERT) and the cleft lip and palate transmembrane 1-like (CLPTM1L) genes, which have been implicated in carcinogenesis. TERT, as the rate-limiting catalytic subunit of the telomerase enzyme, plays an essential role in the maintenance of telomere DNA length, chromosomal stability, and cellular immortality [3]. TERT is expressed in the vast majority of human malignant celllines and tumors but not in the corresponding benign tissues [4]. With regard to NPC, elevated telomerase activity and TERT expression levels were detected in tumor tissues compared with normal nasopharyngeal epithelium and tissues [5, 6]. Furthermore, for the clinicopathological features in NPC patients, the telomerase activity was observed more frequently in advanced clinical stage and lymph node metastasis [7, 8]. The other gene in this region, CLPTM1L, which is named for its similarity to a gene implicated in the susceptibility to cleft lip palate, was identified through screening for cisplatin-resistance-related genes. Although the function of this gene is less known, it was found to be upregulated in cisplatin-resistant ovarian tumor cell lines and to induce apoptosis in cisplatin-sensitive cells [9]. On the basis of the above functional relevance of the TERT and CLPTM1L genes in the pathogenesis of cancers, we hypothesize that TERT and CLPTM1L may be excellent biological candidate susceptibility genes for cancers.
Single nucleotide polymorphisms (SNPs) in the TERT-CLPTM1L locus on chromosome 5p15.33 have been linked to a spectrum of cancers. Multiple independent genome-wide association studies (GWAS) have shown that SNPs in the TERT-CLPTM1L locus were significantly associated with the risk of many different types of cancer, including basal cell carcinoma, glioma, lung, cervical, prostate, bladder, pancreatic, chronic lymphocytic leukemia, and testicular germ cell cancers [10–17]. Association studies on candidate gene strategy also reported that polymorphisms in the TERT-CLPTM1L locus were associated with serous ovarian, breast, endometrial, lung, and colorectal cancers [18–20]. The roles of the polymorphism in the TERT-CLPTM1L locus in NPC, however, have not been fully investigated. In the present study, we examined whether the polymorphism in the TERT-CLPTM1L locus have any bearing on the risk or severity of NPC in Chinese populations.
RESULTS
Genetic association between polymorphisms in the TERT-CLPTM1L locus and NPC
Among the 26 tag SNPs chosen first, four SNPs were eliminated because of low frequency (rs2853691), poor genotyping reliability (rs2736118), or departure from Hardy-Weinberg equilibrium (P < 0.05) (rs2853668 and rs10073340). We therefore selected the remaining 22 SNPs for the subsequent analyses.
The genotyping results for the 22 SNPs are shown in Table 1 and Figure 1. The observed genotype frequencies for the 22 polymorphisms conformed to Hardy–Weinberg equilibrium in the Guangxi population (all P > 0.05, data not shown). On the basis of logistic regression analysis with adjustment for age, sex, smoking status, alcohol use, and family history, significant associations with the susceptibility to NPC were observed with the 12 SNPs in the Guangxi population (Table 1). After multiple corrections, four SNPs (rs2736098, rs2735845, rs402710, and rs401681) were significantly associated with NPC risk (P < 0.0029 after SNP SpD correction, Meff Li value = 17, Table 1). We also investigated the association between the 22 SNPs and the severity of NPC (as measured by TNM staging system) in the Guangxi population. After multiple corrections, no association was observed between the 22 SNPs and the severity of NPC (Table S8).
Table 1: Associations between tag SNPs on chromosome 5p15.33 and risk of nasopharyngeal carcinoma in the Guangxi population
Abbreviations: OR, odds ratio; CI, confidence interval.
aFrequency of the genotype in the patients/frequency of the genotype in the controls.
bThe ORs and P values were adjusted for age, sex, smoking and drinking status, smoking level, nationality, and family history of cancer in the additive model.
*The P value remained significant after correction for multiple comparisons.
Figure 1: Regional plots for associations at the TERT-CLPTM1L locus on chromosome 5p15.33. (A) Significance of the each tag SNPs. Shown are the −log10 association P values of SNPs in the Guangxi population. The intensity of the red shading indicates the strength of the linkage disequilibrium (LD) with the index SNP (rs2735845). Also shown are the SNP built 36 coordinates in kilobases (kb), the recombination rates in centimorgans (cM) per megabase (Mb) (in blue) and the genes in the region (in green). (B) Pairwise LD between tag SNPs at the TERT-CLPTM1L locus. The value within each diamond represents the pairwise correlation between tag SNPs (measured as D’) defined by the upper left and the upper right sides of the diamond. The diamond without a number corresponds to D’ = 1. Shading represents the magnitude and significance of pairwise LD, with a red-to-white gradient reflecting higher to lower LD values.
To confirm the initial associations of these four SNPs with NPC, the Guangdong population, composed of 997 patients and 972 controls, was also genotyped. The observed genotype frequencies for the four polymorphisms conformed to Hardy-Weinberg equilibrium in the Guangdong population (all P > 0.05, data not shown). There was no association between the rs2736098 and rs402710 polymorphisms and the risk of NPC in the Guangdong population (data not shown). However, the rs2735845 and rs401681 polymorphisms were significantly associated with risk of NPC in the Guangdong population (Table 2). For the 2735845 polymorphism, the G allele presented significant association with increased risk of NPC in the additive model (OR = 1.19, 95% CI = 1.04–1.37, P = 0.011). For the rs401681 polymorphism, the T allele presented significant association with decreased risk of NPC in the additive model (OR = 0.85, 95% CI = 0.74–0.99, P = 0.034). When the Guangxi and Guangdong populations were combined, stronger associations were observed between both the rs2735845 and rs401681 polymorphisms and NPC risk in the additive model (for rs2735845: adjusted OR = 1.23, 95% CI = 1.12–1.35, P = 4.64 × 10−5; for rs401681: adjusted OR = 0.81, 95% CI = 0.74–0.90, P = 1.00 × 10−4) (Table 2). The effects of rs2735845 and rs401681 were maintained when corrected for each other by logistic regression (OR = 1.15, 95% CI = 1.03–1.29, P = 0.012 and OR = 0.88, 95% CI = 0.78–0.99, P = 0.03, respectively). No association more significant than that of rs2735845 and rs401681 was identified through imputation of other polymorphisms in the TERT-CLPTM1L locus based on HapMap CHB data (data not shown). Furthermore, the correlation between rs2735845 and rs401681 is not strong (D’ = 0.91, r2 = 0.24 in HapMap CHB samples, D’ = 0.89, r2 = 0.24 in our control data; Figure 1). Taken together, these data strongly suggest that rs2735845 and rs401681 are independently associated with NPC.
Table 2: The genotype frequencies of rs2735845 and rs401681 polymorphisms in patients and controls of the Guangxi and Guangdong population
aFrequency of the genotype in the patients/frequency of the c genotype in the controls.
bThe ORs and P values were calculated by logistic regression analysis adjusted for age, sex, smoking and drinking status, smoking level, nationality (only in Guangxi population), and family history of cancer.
The associations between the rs2735845 and rs401681 polymorphisms and the susceptibility to NPC in both the Guangxi and Guangdong populations were further examined with stratification by age, sex, smoking status, smoking level, drinking status and family of history. For the rs401681 polymorphism, although the susceptibility to NPC seemed to be more pronounced in some subgroups, these differences could be attributed to chance (all P > 0.05, test for homogeneity; data not shown). For the rs2735845 polymorphism, the general pattern of results was similar (data not shown).
Effects of the rs2735845 and rs401681 polymorphisms on TERT and CLPTM1L mRNA expression
We evaluated the association between genotypes of the rs2735845 and rs401681 polymorphisms and TERT and CLPTM1L mRNA levels in 12 fetal mesenchymal stem cells (fMSCs). As shown in Figure 2, the relative CLPTM1L mRNA level from fMSCs with rs2735845 G allele (GG + CG genotypes) was significantly higher than that from fMSCs with CC genotype (t test; P = 0.02). However, there was no association between CLPTM1L mRNA level and rs401681 genotypes in fMSCs (data not shown).
Figure 2: Quantitative real-time RT-PCR for CLPTM1L. The expression of CLPTM1L for three different genotypes of rs2735845 was measured in RNA from fetal mesenchymal stem cells from human umbilical cord obtained from a group of 12 normal full-term-delivery women. Normalization for mRNA quantity was performed with human GAPDH control primers for each sample, and the final abundance values were adjusted to yield an arbitrary value of 1 for the rs2735845 CC genotype. Columns include the mean from triplicate measurements; the vertical bars indicate the standard deviation. Compared to the CC carriers, the G allele carriers exhibited markedly elevated CLPTM1L transcription (P = 0.02; t test).
Greater TERT and CLPTM1L expression in NPC tissues compared with that in non-cancerous nasopharyngeal tissues
We investigated the protein expression of TERT and CLPTM1L in NPC tissues and non-cancerous nasopharyngeal tissues by immunohistochemical (IHC) assay (Figure 3 and Table S1 in Supplementary Materials). Significantly higher expression of TERT was observed in the NPC tissues compared to that in the non-cancer nasopharyngeal tissues (P < 0.001). There was also significantly higher expression of CLPTM1L in the NPC tissues compared to that in the non-cancer nasopharyngeal tissues (P < 0.001). However, there was no association between the genotypes of the rs2735845 and rs401681 polymorphisms and the expression of TERT and CLPTM1L in NPC tissues and non-cancerous nasopharyngeal tissues.
Figure 3: Protein expression of TERT and CLPTM1L by immunohistochemical staining in representative nasopharyngeal carcinoma tissues and non-tumor nasopharyngeal tissues. Images in the box (left, 100×) were enlarged and are shown in the right box (400×). The scale bar represents 50 μm.
DISCUSSION
In this study, we investigated the associations of 26 tag SNPs in the TERT-CLPTM1L locus on 5p15.33 with the risk of occurrence and progression of NPC in populations in southern China. We found that two SNPs, rs2735845 and rs401681, were significantly associated with the onset of NPC in two independent populations. Functional analysis showed that the rs2735845 G allele was associated with higher CLPTM1L mRNA levels when compared with the C allele. In addition, we also found that the expression of TERT and CLPTM1L was higher in the NPC tissues compared to the non-cancer nasopharyngeal tissues. Our data, together with recent evidence for an association of the polymorphisms in the TERT-CLPTM1L locus with the risk of lung and colorectal cancer [19, 20], suggest that the polymorphisms at this locus are potent genetic risk factors for the onset of NPC.
The genetic associations observed in this study are biologically plausible. There are two genes within the 5p15.33 TERT-CLPTM1L locus through which the association might plausibly be mediated: TERT and CLPTM1L. TERT translates the catalytic protein subunit of the telomerase complex, which elongates the telomere in cells. Telomere-driven genome instability is a widespread cause of genome instability in cancer and is thought to be the crucial event in the development of tumorigenesis [3]. It has been shown that small hairpin RNA can direct against TERT to inhibit cell viability by regulating telomerase activity and its related protein expression in NPC cells [21]. Moreover, overexpression of TERT has been detected in a variety of cancer cells, including NPC [6]. In this study, we also found that the expression of TERT was higher in the NPC tissues compared to the non-cancer nasopharyngeal tissues. Our results, together with those of the earlier studies, suggest that TERT plays an important role in NPC carcinogenesis. CLPTM1L encodes a predicted transmembrane protein, which plays an important role in carcinogenesis. CLPTM1L is upregulated in a cisplatin-resistant ovary cancer cell line and might be involved in the apoptotic response of cells under genotoxic stress caused by cisplatin [9]. CLPTM1L has been reported in the cellular response to genotoxic stress and cisplatin resistance [9]. In addition, over-expression of CLPTM1L has been observed in lung tumors [22]. In this study, we found that the expression of CLPTM1L was higher in the NPC tissues compared to the non-cancer nasopharyngeal tissues, suggesting that the CLPTM1Lhas been implicated in NPC carcinogenesis.
Chromosome 5p15.33 is one of the strongest and most consistent association signals for cancer risk [23]. Both GWAS and candidate gene association studies have identified common variants in 5p15.33 that reproducibly associate with several types of cancer. In 2008, McKay et al. carried out a GWAS of lung cancer, and two independent polymorphisms, rs402710 and rs2736100, were detected [12]. Another GWAS of lung cancer identified another associated polymorphism: rs401681 [11]. Subsequently, Rafnar et al. found that the rs401681 polymorphism was associated with the risk of basal cell, lung, urinary bladder, prostate and cervical cancers and cutaneous melanoma, and rs2736098, a synonymous coding SNP in exon 2 of the TERT gene, was associated with the risk of basal cell, lung, urinary bladder, and prostate cancers [10, 24]. Additionally, a meta-analysis of two GWAS of glioma demonstrated that the rs2736100 polymorphism was associated with the risk of glioma [16]. Petersen et al. conducted a GWAS of pancreatic cancer, and rs401681 was one of the most significantly associated polymorphisms [14]. GWAS of a testicular germ cell tumor identified two independently associated polymorphisms in the TERT-CLPTM1L locus, rs4635969 and rs2736100 [15]. Association studies of candidate gene strategies also reported that SNPs in TERT were associated with serous ovarian, breast, and endometrial cancers [25–27]. With regard to NPC, Fachiroh et al. and Ko et al. have reported that certain polymorphisms in the TERT gene were associated with NPC risk in Thailand [28, 29]. In this study, we found that two independent SNPs, rs2735845 and rs401681, were significantly associated with the risk of NPC. Our results, together with those of the earlier studies, suggest that the TERT-CLPTM1L locus may be a susceptibility region for NPC.
The molecular mechanism by which the rs401681 and rs2735845 polymorphisms confer cancer risk is not well characterized. The rs2735845 polymorphism, which resides in the intergenic region, was never reported before. The rs401681 polymorphism, which resides in an LD block that contains the CLPTM1L gene and the 5′ end of the TERT gene, is one of the most extensively studied polymorphisms in the TERT-CLPTM1L locus. Rafnar et al. reported that the rs401681 polymorphism may lead to an increase in the gradual shortening of telomeres over time, but the effect may only become apparent after a certain age [10]. In this study, we did not find any association between rs401681 genotypes and TERT and CLPTM1L mRNA levels. We found that the CLPTM1L mRNA levels transcribed by the rs2735845 G allele were higher than those transcribed by the CC genotype in fMSCs. However, the difference in the induction of the reporter gene expression between the rs2735845 two alleles is not statistically significant (data not shown). Further studies are warranted to elucidate the molecular mechanism of these two polymorphisms.
When reviewing the results of this study, several potential limits should be kept in mind. First, in this hospital-based study, our patients were enrolled from hospitals, and the controls were selected from the community population; thus, inherent selection bias cannot be completely excluded. However, through further adjustment and stratification in the data analyses, the potential confounding effect might have been minimized. Second, a tag SNP approach has been shown to improve the power of the study for common genetic variants; however, the modest effects from rare genetic variants may remain undetected. Deep re-sequencing of the TERT-CLPTM1L locus may provide further help to uncover additional associated variants and to facilitate selection of potential causal variants for further functional studies.
In conclusion, our study shows an association between the genetic variations in the TERT-CLPTM1L locus and the risk of occurrence of NPC in two Chinese subpopulations. These findings may provide support for the importance of the TERT-CLPTM1L locus in the pathogenesis of NPC. If confirmed by other studies, knowledge of genetic factors contributing to the pathogenesis of the NPC as presented here may have implications for the screening and treatment of this disorder.
MATERIALS AND METHODS
Study population
This study consisted of two populations of patients with NPC and controls residing in Guangxi and Guangdong provinces, both of which are located in southern China (see Supplementary Materials, Table S2). The Guangxi population, which was used for the 1st stage scanning analysis, consisted of 855 patients with NPC and 1036 controls. All subjects were unrelated ethnic Han Chinese and residents in Nanning city and the surrounding regions. The diagnosis of patients, the inclusion and exclusion criteria for patients and controls, and the definition of smokers and drinkers were described in detail in our previous studies [30]. Briefly, the patients with NPC were newly diagnosed and pathologically confirmed, and they were consecutively recruited at the Guangxi Cancer Hospital (Nanning, China) between September 2003 and January 2008. The response rate for case patients was 95%. Patients who received chemotherapy or radiotherapy before surgery or had another type of cancer were excluded from this study. The controls were randomly selected from a community cancer screening program for early detection of cancer conducted in the same regions during the same time period in which the NPC patients were collected. The response rate for control subjects was 91%. The selection criteria for the controls included no individual history of cancer. The Guangdong population, which was used for the 2nd stage replication analysis, consisted of 997 patients and 972 controls. The diagnosis of patients and the inclusion and exclusion criteria for patients and controls were described in detail previously [31]. Briefly, the patients with NPC were newly diagnosed and pathologically confirmed, and they were consecutively recruited from Sun Yat-sen University Cancer Center (SYSUCC, Guangzhou, China) between October 2005 and October 2007. Only the patients who were self-reported as Guangdong Chinese and who lived in Guangdong province were recruited. During the same period, healthy controls were recruited from physical examination centers of several large comprehensive hospitals in the local communities of Guangdong. The selection criteria for the controls included no individual history of cancer. The response rate for cases and control was 90%. For both populations, tumor staging was performed according to the tumor-node-metastasis (TNM) classification of the 1997 American Joint Committee on Cancer (AJCC) system. All TNM classifications were determined by senior pathologists of the hospital on the basis of the postoperative histopathologic examination.
Of the patients with NPC in the Guangxi population, 41 patients who had undergone resection before receiving any further treatment at Guangxi Cancer Hospital were selected, and primary NPC biopsies were collected from them (see Supplementary Materials, Table S3). All 41 tumor tissues were poorly differentiated squamous cell carcinoma (SCC). Histological non-cancerous nasopharyngeal epithelium tissues were collected from 13 of the 1036 control subjects in the Guangxi population (see Supplementary Materials, Table S3). All of the tissues were fixed in paraformaldehyde, embedded in paraffin wax, and prepared for IHC staining.
For TERT and CLPTM1L expression assays, human umbilical cord (UC) was obtained from a group of 12 normal full-term-delivery women. At recruitment, informed consent was obtained from all participants, and personal information about demographic factors and medical history and tobacco and alcohol use was collected via structured questionnaire. This study was performed with the approval of the Ethical Committee of Beijing Institute of Radiation Medicine and conducted according to the principles expressed in the Declaration of Helsinki.
Selection of SNPs
Using the tagger program implemented in Haploview 4.1 [32], tag SNPs across the 115-kb of TERT-CLPTM1L locus on 5p15.33 (chr5:1,301,873–1,417,242) were selected on the basis of linkage disequilibrium (LD) patterns observed in the population of Han Chinese in Beijing, China (CHB), and genotyped as part of the International HapMap Project (data Rel 24/phase II, http://www.hapmap.org) (see Supplementary Materials, Figure S1). Two well-studied functional SNPs in the TERT gene (rs2853669 and rs2735940) and six SNPs (rs2736100, rs2736098, rs4975616, rs402710, rs401681, rs31489) reported in the GWAS studied were directly included as tag SNPs using the forced inclusion option [10–12, 33, 34]. In addition, 18 additional tag SNPs (rs3803309, rs4983540, rs3809454, rs3803304, rs2494738, rs2494743, rs2494750, rs4983559, rs2498806, rs10139525, rs8010024 and rs2582519) were selected with a pairwise r2 of > 0.80 and a minor allele frequency of ≥ 0.01 on the basis of HapMap CHB data. Overall, 26 tag SNPs in the TERT-CLPTM1L locus were selected for further genotyping (see Supplementary Materials, Table S4).
Genotype analysis
In the stage I association study, 15 SNPs were genotyped using the MassArray System (Sequenom) with primers (see Supplementary Materials, Table S6) that were described previously [35]. All 15 SNPs had genotyping call rates of 95% or better. Nine SNPs were genotyped using the SNPstream ultra-high throughput genotyping system (Beckman Coulter, Fullerton, CA) according to the manufacturer’s instructions, as described previously [36]. The primer sequences are listed in Table S6 in Supplementary Materials 1. The s2736118 polymorphism was eliminated because of poor genotyping reliability. The remaining eight SNPs had genotyping call rates of 95% or better. The rs2736100 polymorphism was genotyped by TaqMan MGB assay (Applied Biosystems, USA) according to the manufacturer’s protocol, as described previously [30]. The primer and TaqMan probes used for the rs2736100 polymorphism were as follows: 5′-AGACGGGGAACAAAGGAGGA-3′ (forward primer), 5′-AGTTCTATCTCAGGCATCTTGACACC-3′ (reverse primer), 5′-FAM-CAAAGCTAAAGAAACAC-MGB-3′, 5′-HEX-CAAAGCTACAGAAACA-MGB-3′. The rs2853669 polymorphism was genotyped using polymerase chain reaction-based restriction fragment length polymorphism (PCR-RFLP) analysis. The primers and the restriction enzymes used are listed in Table S7 in Supplementary Materials 1. The PCR conditions were identical to those used for the SNP discovery described previously [37].
In the stage II association study, four SNPs (rs2736098, rs2735845, rs401710, and rs401681) were further genotyped in the Guangdong population. The genotyping data of the rs401710 and rs401681 polymorphisms were collected from Human610-Quad data published in Bei et al.’s study [31]. The rs2735845 and rs2736098 polymorphisms were genotyped by PCR-RFLP analysis (Table S6 in Supplementary Materials). The PCR conditions were identical to those for the rs2853669 polymorphism.
Genotyping was performed in a blind manner such that the performers did not know the subjects’ patient or control status. For quality control, a 15% masked, random sample of patients and controls was tested twice by different people, and all results exhibited 100% concordance. We also designed PCR-RFLP assays to validate two random SNPs to compare with the results of the Sequneom and SNPstream methods, respectively, in 5% of the patient/control samples; the reproducibility was 98.8%. In addition, genotypes identified by PCR-RFLP were confirmed by DNA sequencing.
Isolation and expansion of fMSCs from UC
UC (n = 12) was obtained from normal full-term-delivery women. Briefly, after disinfection with 75% ethanol, the fresh UC was transferred to PBS where excess blood on the surface was removed. After removal of blood vessels, 2 cm of UC was sliced into very small pieces followed by treatment with 0.1% collagenase II (Gibco, USA), containing antibiotic solution (100 U/mL penicillin and 100 μg/mL streptomycin; CSPC, China) in Dulbecco’s Modified Eagle Medium/F12 (DMEM/F12) (Gibco) overnight at 37°C. The cells were cultured in DMEM/F12 supplemented with 10% fetal bovine serum (Hyclone, China) in 5% CO2 in a 37°C incubator for 3–4 days to allow cells to adhere. The adherent cells formed colonies and grew rapidly, exhibiting spindle-shaped morphology. The medium was replaced every three days. When the cells reached to 60–80% confluence, they were harvested with trypsin-EDTA solution (Neuronbc, China) and subcultured at a density of 3–6×103 cells/cm2 [38]. Amplification of pluripotent stem cell markers, such as Oct-4 and Rex-1, was used to demonstrate primitive properties of the cultured cells, as described previously [39].
Real-time analysis of TERT and CLPTM1L messenger RNA (mRNA) expression
Total RNA was extracted from 12 fetal MSCs using Trizol reagent (Invitrogen, USA) according to the manufacturer’s protocol. Reverse transcription was performed using a reverse transcriptase system (Promega, USA) according to the manufacturer’s instructions. Quantitative real-time PCR for CLPTM1L and full-length TERT, with GAPDH as an internal reference gene, was performed using the Bio-Rad IQ5 (Bio-Rad, USA) based on the SYBR-Green method. The primers used for TERT were 5′-TGTACTTTGTCAAGGTGGATGTG-3′ and 5′-GTAAGGCTGGAGGTCTGTCAA-3′ [40], those used for CLPTM1L were 5′-TAAGAGCTGGTACTCCTGGT-3′ and 5′-GTCATCAATGAAGGTGTTGA-3′, and those used for GAPDH were 5′-CATGAGAAGTATGACAACAGCCT-3′ and 5′-ACTCCTTCCACGATACCAAAGT-3′. All reactions performed in a 15-μL reaction volume containing 7.5 μL of 2X SYBR Green Master mix (Bio-Rad, USA), 1 μL of each primer pair (6 μmol/L), and 1 μL of cDNA templates. Cycling conditions were 95°C for 2 minutes, followed by 40 cycles at 95°C for 10 seconds, 54°C for 15 seconds, and 72°C for 15 seconds. The data were analyzed using the 2−ΔΔCt method [41]. All analyses were performed in a blinded fashion with the laboratory persons unaware of the genotyping data.
IHC staining analysis
NPC tissues (n = 41) and non-tumor nasopharyngeal tissues (n = 13), which had been fixed with paraformaldehyde and embedded in paraffin wax, were analyzed for protein expression of TERT and CLPTM1L. The IHC assay was performed as previously described [42]. Briefly, two slides of each biopsy were stained with hematoxylin-eosin for routine histological evaluation. The slides were washed in xylene to remove the paraffin and then rehydrated through serial dilutions of alcohol, followed by washings with a solution of PBS (pH 7.2). All subsequent washes were buffered via the same protocol. The slides were then incubated with 3% H2O2 for 10 min to reduce non-specific staining. The treated sections were then placed in a citrate buffer (pH 6.0) and heated in a pressure cooker for two minutes. The samples were then incubated with a polyclonal anti-human CLPTM1L antibody (Sigma Aldrich, HPA014791) or a monoclonal rabbit anti-human TERT antibody (Lifespan, LS-C49516) overnight at 4°C. The conventional streptavidin peroxidase method (MaxVision™ HRP-Polymer IHC Kit; Maxim Co., China) was performed for signal development, and the cells were counter-stained with hematoxylin. Negative controls and positive controls were prepared at the same time. The slides were mounted with gum for examination and capture by the Olympus BX51 microscopic/Digital Camera System for study comparison. IHC signals were scored as previously described [42]. Slides were scored by two pathologists (Wu JH and Wang AC) who did not have knowledge of the genotype results of rs2735845 and rs401681 polymorphisms or of the patient outcomes.
Statistical analysis
Genotype and allele frequencies for polymorphisms were determined by gene counting. The fitness according to the Hardy-Weinberg equilibrium was tested using the χ2 test. Associations between polymorphisms and risk of NPC were estimated by use of logistic regression analyses. Odds ratios (ORs) and 95% confidence intervals (CIs) were used to measure the strength of association. To correct for multiple testing, we used spectral decomposition (SpD) of matrices of LD between SNPs, which adjusts for multiple testing while taking into account LD among the tested SNPs [43]. This method showed that P values of 0.0029 and below can be considered statistically significant after correction for multiple testing. These analyses were performed using SPSS software (version 17.0, SPSS Inc., Chicago, IL, USA). The pairwise LD calculation (Lewontin’s D and r2) and haplotype block construction were performed using the program Haploview 4.2 [32]. The haplotype analysis was performed using PLINK software (version 1.07) (http://pngu.mgh.harvard.edu/purcell/plink/) [44].
ACKNOWLEDGMENTS
We thank all of the tested individuals, their families, and the collaborating clinicians for their participation.
FUNDING
This work was supported in part by the China National Funds for Distinguished Young Scientists (No. 81125017) and the National Natural Science Foundation of China (No. 81222027, 30901707 and 30901231).
CONFLICTS OF INTEREST
The authors declare that they have no conflicts of interests.
REFERENCES
1. Wei KR, Yu YL, Yang YY, Ji MF, Yu BH, Liang ZH, Reng X. Epidemiological trends of nasopharyngeal carcinoma in China. Asian Pac J Cancer Prev. 2010; 11:29–32.
2. Chou J, Lin YC, Kim J, You L, Xu Z, He B, Jablons DM. Nasopharyngeal carcinoma—review of the molecular mechanisms of tumorigenesis. Head Neck. 2008; 30:946–963.
3. Murnane JP. Telomere loss as a mechanism for chromosome instability in human cancer. Cancer Res. 2010; 70:4255–4259.
4. Shay JW, Bacchetti S. A survey of telomerase activity in human cancer. Eur J Cancer. 1997; 33:787–791.
5. Wang X, Xiao J, Zhao S, Tian Y, Wang G. [Expression of telomerase subunits and its relationship with telomerase activity in nasopharyngeal carcinoma]. Zhonghua Yi Xue Za Zhi. 2001; 81:553–556.
6. Yeh YM, Pan YT, Wang TC. Cdc42/Rac1 participates in the control of telomerase activity in human nasopharyngeal cancer cells. Cancer Lett. 2005; 218:207–213.
7. Cheng RY, Yuen PW, Nicholls JM, Zheng Z, Wei W, Sham JS, Yang XH, Cao L, Huang DP, Tsao SW. Telomerase activation in nasopharyngeal carcinomas. Br J Cancer. 1998; 77:456–460.
8. Peng H, Wang X, Yang G, Jiang S, Zhao T. [The study on the role of telomerase activity in human nasopharyngeal carcinoma]. Zhonghua Er Bi Yan Hou Ke Za Zhi. 2000; 35:289–291.
9. Yamamoto K, Okamoto A, Isonishi S, Ochiai K, Ohtake Y. A novel gene, CRR9, which was up-regulated in CDDP-resistant ovarian tumor cell line, was associated with apoptosis. Biochem Biophys Res Commun. 2001; 280:1148–1154.
10. Landi MT, Chatterjee N, Yu K, Goldin LR, Goldstein AM, Rotunno M, Mirabello L, Jacobs K, Wheeler W, Yeager M, Bergen AW, Li Q, Consonni D, et al. A genome-wide association study of lung cancer identifies a region of chromosome 5p15 associated with risk for adenocarcinoma. Am J Hum Genet. 2009; 85:679–691.
11. 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, et al. Common 5p15.33 and 6p21.33 variants influence lung cancer risk. Nat Genet. 2008; 40:1407–1409.
12. 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, et al. Lung cancer susceptibility locus at 5p15.33. Nature genetics. 2008; 40:1404–1406.
13. Landi MT, Chatterjee N, Yu K, Goldin LR, Goldstein AM, Rotunno M, Mirabello L, Jacobs K, Wheeler W, Yeager M, Bergen AW, Li Q, Consonni D, et al. A genome-wide association study of lung cancer identifies a region of chromosome 5p15 associated with risk for adenocarcinoma. Am J Hum Genet. 2009; 85:679–691.
14. Petersen GM, Amundadottir L, Fuchs CS, Kraft P, Stolzenberg-Solomon RZ, Jacobs KB, Arslan AA, Bueno-de-Mesquita HB, Gallinger S, Gross M, Helzlsouer K, Holly EA, Jacobs EJ, et al. A genome-wide association study identifies pancreatic cancer susceptibility loci on chromosomes 13q22.1, 1q32.1 and 5p15.33. Nat Genet. 2010; 42:224–228.
15. Turnbull C, Rapley EA, Seal S, Pernet D, Renwick A, Hughes D, Ricketts M, Linger R, Nsengimana J, Deloukas P, Huddart RA, Bishop DT, Easton DF, et al. Variants near DMRT1, TERT and ATF7IP are associated with testicular germ cell cancer. Nature genetics. 2010; 42:604–607.
16. Shete S, Hosking FJ, Robertson LB, Dobbins SE, Sanson M, Malmer B, Simon M, Marie Y, Boisselier B, Delattre JY, Hoang-Xuan K, El Hallani S, Idbaih A, et al. Genome-wide association study identifies five susceptibility loci for glioma. Nature genetics. 2009; 41:899–904.
17. Speedy HE, Di Bernardo MC, Sava GP, Dyer MJ, Holroyd A, Wang Y, Sunter NJ, Mansouri L, Juliusson G, Smedby KE, Roos G, Jayne S, Majid A, et al. A genome-wide association study identifies multiple susceptibility loci for chronic lymphocytic leukemia. Nature genetics. 2014; 46:56–60.
18. Carvajal-Carmona LG, O’Mara TA, Painter JN, Lose FA, Dennis J, Michailidou K, Tyrer JP, Ahmed S, Ferguson K, Healey CS, Pooley K, Beesley J, Cheng T, et al. Candidate locus analysis of the TERT-CLPTM1L cancer risk region on chromosome 5p15 identifies multiple independent variants associated with endometrial cancer risk. Human genetics. 2015; 134:231–245.
19. Zhong R, Liu L, Zou L, Zhu Y, Chen W, Zhu B, Shen N, Rui R, Long L, Ke J, Lu X, Zhang T, Zhang Y, et al. Genetic variations in TERT-CLPTM1L locus are associated with risk of lung cancer in Chinese population. Mol Carcinog. 2013; 52:E118–126.
20. Kinnersley B, Migliorini G, Broderick P, Whiffin N, Dobbins SE, Casey G, Hopper J, Sieber O, Lipton L, Kerr DJ, Dunlop MG, Tomlinson IP, Houlston RS. The TERT variant rs2736100 is associated with colorectal cancer risk. Br J Cancer. 2012; 107:1001–1008.
21. Shen Y, Wang Y, Chen S, Xiao B, Su J, Tao Z. [The effect of shRNA targeting hTERT on telomerase and the expression of PCNA and Caspase-3 in nasopharyngeal carcinoma cells]. Lin Chung Er Bi Yan Hou Tou Jing Wai Ke Za Zhi. 2008; 22:411–415.
22. James MA, Wen W, Wang Y, Byers LA, Heymach JV, Coombes KR, Girard L, Minna J, You M. Functional characterization of CLPTM1L as a lung cancer risk candidate gene in the 5p15.33 locus. PLoS One. 2012; 7:e36116.
23. Ahmadiyeh N, Pomerantz MM, Grisanzio C, Herman P, Jia L, Almendro V, He HH, Brown M, Liu XS, Davis M, Caswell JL, Beckwith CA, Hills A, et al. 8q24 prostate, breast, and colon cancer risk loci show tissue-specific long-range interaction with MYC. Proc Natl Acad Sci USA. 2010; 107:9742–9746.
24. Stacey SN, Sulem P, Masson G, Gudjonsson SA, Thorleifsson G, Jakobsdottir M, Sigurdsson A, Gudbjartsson DF, Sigurgeirsson B, Benediktsdottir KR, Thorisdottir K, Ragnarsson R, Scherer D, et al. New common variants affecting susceptibility to basal cell carcinoma. Nat Genet. 2009; 41:909–914.
25. Johnatty SE, Beesley J, Chen X, Macgregor S, Duffy DL, Spurdle AB, deFazio A, Gava N, Webb PM, Rossing MA, Doherty JA, Goodman MT, Lurie G, et al. Evaluation of candidate stromal epithelial cross-talk genes identifies association between risk of serous ovarian cancer and TERT, a cancer susceptibility “hot-spot”. PLoS Genet. 2010; 6:e1001016.
26. Shen J, Gammon MD, Wu HC, Terry MB, Wang Q, Bradshaw PT, Teitelbaum SL, Neugut AI, Santella RM. Multiple genetic variants in telomere pathway genes and breast cancer risk. Cancer Epidemiol Biomarkers Prev. 2010; 19:219–228.
27. Prescott J, McGrath M, Lee IM, Buring JE, De Vivo I. Telomere length and genetic analyses in population-based studies of endometrial cancer risk. Cancer. 2010; 116:4275–4282.
28. Fachiroh J, Sangrajrang S, Johansson M, Renard H, Gaborieau V, Chabrier A, Chindavijak S, Brennan P, McKay JD. Tobacco consumption and genetic susceptibility to nasopharyngeal carcinoma (NPC) in Thailand. Cancer Causes Control. 2012; 23:1995–2002.
29. Yee Ko JM, Dai W, Wun Wong EH, Kwong D, Tong Ng W, Lee A, Kai Cheong Ngan R, Chung Yau C, Tung S, Li Lung M. Multigene pathway-based analyses identify nasopharyngeal carcinoma risk associations for cumulative adverse effects of TERT-CLPTM1L and DNA double-strand breaks repair. International journal of cancer Journal international du cancer. 2014; 135:1634–1645.
30. Ma F, Zhang H, Zhai Y, Huang W, Zhao C, Ou S, Zhou H, Yuan W, Wang Z, Wang H, Yue W, Yu L, Li P, et al. Functional polymorphism −31C/G in the promoter of BIRC5 gene and risk of nasopharyngeal carcinoma among chinese. PLoS One. 2011; 6:e16748.
31. Bei JX, Li Y, Jia WH, Feng BJ, Zhou G, Chen LZ, Feng QS, Low HQ, Zhang H, He F, Tai ES, Kang T, Liu ET, et al. A genome-wide association study of nasopharyngeal carcinoma identifies three new susceptibility loci. Nat Genet. 2010; 42:599–603.
32. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005; 21:263–265.
33. Hsu CP, Hsu NY, Lee LW, Ko JL. Ets2 binding site single nucleotide polymorphism at the hTERT gene promoter—effect on telomerase expression and telomere length maintenance in non-small cell lung cancer. Eur J Cancer. 2006; 42:1466–1474.
34. Matsubara Y, Murata M, Yoshida T, Watanabe K, Saito I, Miyaki K, Omae K, Ikeda Y. Telomere length of normal leukocytes is affected by a functional polymorphism of hTERT. Biochem Biophys Res Commun. 2006; 341:128–131.
35. Buetow KH, Edmonson M, MacDonald R, Clifford R, Yip P, Kelley J, Little DP, Strausberg R, Koester H, Cantor CR, Braun A. High-throughput development and characterization of a genomewide collection of gene-based single nucleotide polymorphism markers by chip-based matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Proc Natl Acad Sci USA. 2001; 98:581–584.
36. Zhu X, Wang Y, Zhang H, Liu X, Chen T, Yang R, Shi Y, Cao W, Li P, Ma Q, Zhai Y, He F, Zhou G, et al. Genetic variation of the human alpha-2-Heremans-Schmid glycoprotein (AHSG) gene associated with the risk of SARS-CoV infection. PLoS One. 2011; 6:e23730.
37. Zhou G, Zhai Y, Dong X, Li Y, Zhang X, Zhang R, Li S, Li X, He F, Wei H, Chen X, Yao Z, Shen Y, et al. Variants in TNFRSF5 locus and association analysis with Hepatitis B virus (HBV) infection. Hum Mutat. 2004; 23:99–100.
38. Liu H, Zhang J, Liu CY, Wang IJ, Sieber M, Chang J, Jester JV, Kao WW. Cell therapy of congenital corneal diseases with umbilical mesenchymal stem cells: lumican null mice. PLoS One. 2010; 5:e10707.
39. Jo CH, Kim OS, Park EY, Kim BJ, Lee JH, Kang SB, Han HS, Rhee SH, Yoon KS. Fetal mesenchymal stem cells derived from human umbilical cord sustain primitive characteristics during extensive expansion. Cell Tissue Res. 2008; 334:423–433.
40. Lincz LF, Mudge LM, Scorgie FE, Sakoff JA, Hamilton CS, Seldon M. Quantification of hTERT splice variants in melanoma by SYBR green real-time polymerase chain reaction indicates a negative regulatory role for the beta deletion variant. Neoplasia. 2008; 10:1131–1137.
41. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001; 25:402–408.
42. Zhang H, Zhai Y, Hu Z, Wu C, Qian J, Jia W, Ma F, Huang W, Yu L, Yue W, Wang Z, Li P, Zhang Y, et al. Genome-wide association study identifies 1p36.22 as a new susceptibility locus for hepatocellular carcinoma in chronic hepatitis B virus carriers. Nat Genet. 2010; 42:755–758.
43. Nyholt DR. A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other. Am J Hum Genet. 2004; 74:765–769.
44. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007; 81:559–575.
