Further evidence for the existence of major susceptibility of nasopharyngeal carcinoma in the region near HLA-A locus in Southern Chinese
Journal of Translational Medicinevolume 10, Article number: 57 (2012)
Nasopharyngeal carcinoma (NPC) is a multi-factorial malignancy closely associated with environmental factors, genetic factors and Epstein-Barr virus infection. Human leukocyte antigen (HLA) complex, specially the region near HLA-A locus, was regarded as a major candidate region bearing NPC genetic susceptibility loci in many previous studies including two recent genome-wide association (GWA) studies. To provide further evidence for the NPC susceptibility in the region near HLA-A locus based on other previous studies, we carried out a two-stage hospital-based case control association study including 535 sporadic NPC patients and 525 cancer-free control subjects from Guangdong, a high prevalence area of NPC in China.
38 tag SNPs were initially selected by Heploview from the segment around HLA-A locus (from D6S211 to D6S510) and genotyped on GenomeLab SNPstream platform in 206 cases and 180 controls in the stage 1. Subsequently, the stage 1 significant SNPs and 17 additional SNPs were examined on another platform (Sequenom iPlex Assay) in another independent set of study population including 329 cases and 345 controls.
Totally eight SNPs from the segment from D6S211 to D6S510 within HLA complex were found to be significantly associated with NPC. Two of the most significant SNPs (rs9260734 and rs2517716) located near to HLA-A and HCG9 respectively were in strong LD with some other SNPs of this region reported by two previous GWA studies. Meanwhile, Meanwhile, novel independent susceptibility loci (rs9404952, Pcombined = 6.6 × 10-5, OR combined = 1.45) was found to be close to HLA-G.
Therefore, our present study supports that the segment from D6S211 to D6S510 in HLA complex region might contain NPC susceptibility loci which indeed needs to be fully investigated in the future.
Nasopharyngeal carcinoma(NPC)is a sequamous cell carcinoma that originates from the epithelial lining of the nasopharynx. The incidence of NPC is low in most populations around the world but common in the southern part of China and other Southeast Asia region, where the incidence can reach 20 to 50 per 100,000 individuals[1, 2], suggesting a strong association between specific genetic factors and NPC. Since 1974, many studies have indicated that some specific human leukocyte antigen (HLA) haplotypes and genes within the HLA complex are potentially associated with NPC [3–5]. Our previous meta-analysis  showed that populations in geographical areas at higher risk of developing NPC display HLA distribution patterns different and opposite from areas of low incidence. HLA complex is undoubtedly playing an important role in NPC predisposition, either by playing a functional role in modulating an innate and adaptive immune response against EBV, or as a marker of an unrelated predisposition locus in close linkage. Using high-resolution microsatellite mapping, Lu C.C. et al.  verified his previous study  and hypothesized the existence of a NPC susceptibility locus within a 132 Kb segment (from D6S211 to D6S510) containing the HLA-A locus in Taiwanese. Recently two independent research groups [8, 9] conducted GWA studies in Taiwanese and southern Chinese respectively and consistently strengthened this hypothesis. Therefore, to further tested the hypothesis and provide further evidence for the possible existence of major susceptibility of Nasopharyngeal carcinoma in the region near HLA-A locus, we carried out a two-stage case-control association study to focus on the investigation of 132 Kb segment of interest with moderate sample size (535 cases and 525 control) in a southern Chinese population.
Materials and methods
This study was approved in advance by the institutional review Board of Southern Medical University. A total of 535 NPC patients with pathology-based diagnoses and 525 healthy unrelated controls without family or personal history of cancers and other major illnesses such as inflammation, diabetes, SLE, rheumatoid arthritis (RA), etc, were recruited from Jiangmen Center Hospital and Nanfang Hospital, Guangdong Province, between January 2005 and July 2008. The control individuals were frequency matched to the NPC cases by age, gender, geographic region, and ethnicity. Both cases and controls are of Cantonese origin from Southern China. All participants were informed by the aim of the study before participating in the study and required to sign the written informed consents. Peripheral blood samples were collected from all the participants and genomic DNA was extracted from peripheral blood samples using the Tiangen™ Genomic DNA Kit (Tiangen, China) following to the manufacturer's instructions and stored at -80°C before test.
SNP selection, genotyping and quality control
The 132 kb segment between D6S211 and D6S510 (region from 29903 kb to 30050 kb) is our target region, from which 100 tag SNPs were selected based on Han Chinese (CHB) and Japanese (JPT) database in the HapMap Project  by using the Tagger program  to capture SNPs with a minimum minor allele frequency (MAF) of 0.05 and an r2 of 0.80 or greater.
Of 100 tag SNPs picked out, 38 SNPs with C/T allele, which may be available for facilitating further DNA methylation analysis in our another project, were designed and genotyped for 206 cases and 180 controls in the first stage, using GenomeLab SNPstream genotyping platform (Beckman Coulter, Fullerton, California, USA) as previously described. Subsequently, significant SNPs in the first stage, were designed and genotyped again in the second stage in another set of population composed of 329 cases and 345 controls, using the Sequenom iPLEX MassARRAY platform (Sequenom) as previously described. To further increase the density of genetic markers, 17 additional SNPs were selected from those 100 tag SNPs for genotyping as well in the second stage. In the first stage, a sample call rate of over 95% was observed with 100% matching for quality control samples and blind replicates. In the second stage, the DNA sample quality control (QC) threshold was set at 90%.
Deviations from the Hardy-Weinberg equilibrium (HWE) were firstly assessed using the Fisher's exact test or the Chi-square test. SNPs with HWE test p value > 10-3 could be underwent further calculation. Association analysis based on unconditional logistic regression was carried out by figuring the odds ratio (OR) and 95% confidence interval (95% CI) for tag SNPs. The significance level was set at 0.05. All calculations were done with SPSS 13.0 software (SPSS, Chicago, IL).
The linkage disequilibrium (LD) and haplotype block analysis between significant SNPs in this study and other published SNPs located in HLA region (Additional file 1: Table S1) were calculated by Haploview software 4.2 using Han Chinese (CHB) and Japanese(JPT) data presented in the HapMap Project http://www.hapmap.org. D' and r2 were used to describe LD relationship. The confidence interval method was adopted to set the haplotype block. This algorithm was taken from Gabriel et al. , in which 95% confidence bounds on D' were generated and each comparison was labeled "strong LD" (linkage disequilibrium), "inconclusive," or "strong recombination".
Demographic characteristics of both NPC cases and cancer-free controls were presented in Additional file 2: Table S2. Briefly, there was no significant difference in the distribution of age (P = 0.904) and gender (P = 0.874) between the cases and the controls. Among 535 NPC cases, 470 were Poorly differentiated SCC ( Squamous cell carcinoma), 58 Undifferentiated cancer, 7 Differentiated SCC. Both cases and controls are of Cantonese origin from Southern China.
A total of 38 tag SNPs included in the first stage were successfully genotyped. However, of them, five (rs1611650, rs9295822, rs9260759, rs1611442 and rs6457116) were found to have minor allele frequencies (MAFs) of less than 5%, and other two SNPs (Rs3202637 and rs2517681) were found to deviate from HWE. Therefore, the remaining 31 SNPs were further analyzed for genotype and allele frequencies.
As shown in Table 1 there were five SNPs (rs2517716, rs3869062, rs9260734, rs9260475 and rs9404952) showed significant differences in genotype and allele frequency adjusted for age and gender between NPC case and control group in the first stage. Although rs9260475 was unfortunately excluded as it was not in HWE, and rs3869062 showed a lower call rate (less than 90%) in the second stage, rs2517716, rs9260734 and rs9404952 still remained significant. Combination of these variant genotypes of stage 1 and stage 2 clearly showed that rs2517716 (OR combined = 1.62, P combined = 5.9 × 10-7), rs9260734 (OR combined = 1.67, P combined = 3.3 × 10-7), and rs9404952 (OR combined = 1.45, P combined = 6.6 × 10-5) were associated with a significantly increased risk of NPC.
Interestedly, rs9260734 has been previously reported to be associated with NPC risk in two GWA studies [8, 9], but it was using less tag SNPs in a relatively small sample population that we could still detect its association effect. Moreover, rs9260734 and rs2517716 were found to map to a 35 kb region containing HCG9 and HLA-A (Table 1 Figure 1). For further confirming our result, we picked out 17 additional SNPs from this region for subsequent analysis. Of them, five (rs1632882, rs2571400, rs2735085, rs1632902 and rs16896742) displayed statistical significance (Additional file 3: Table S3) as well, strongly indicating the possible existence of candidate NPC susceptibility in this region.
Furthermore, the linkage disequilibrium (LD) analysis was performed according to HapMap CHB and JPT data to measure LD for significant SNPs in our data set vs. significant SNPs in HLA complex from two recent GWA studies (Additional file 1: Table S1). Three LD blocks were formed separately (Figure 2). Some of significant SNPs were assigned into separate LD blocks together with some significant SNPs reported by Taiwan GWAS (rs2517713 and rs2975042) and the southern Chinese GWAS (rs2860580). Rs2517716 was in high LD with rs2517713 in block 2(D' = 1, r2 = 0.7). Rs1632882, rs16896742 and rs2735085 of 17 additional SNPs and two GWAS rs2860580 and rs2975042 formed block 1 and 3. The LD coefficient D' was 1 and r2 was 0.631 between rs2860580 and rs1632882. The LD coefficients (D') ranged from 0.952 to 1 and r2 value ranged from 0.223 to 0.56 between rs2975042, rs16896742 and rs2735085. Therefore, although some of our significant loci were inconsistent with previous GWA studies, our LD analysis displayed strong LD relationships between them.
Rs9404952 as a novel SNP was located at another region downstream from HLA-G which has never been reported to be associated with NPC before. Accordingly, Rs9404952 did not merge with any blocks and the D' and r2 between it and other SNPs were very low (D' from 0.006 to 0.676, r2 from 0 to 0.271), suggesting that it might be an independent susceptibility factor.
Applying a two stage case-control analysis in southern Chinese, we demonstrated that multiple loci within the 132 Kb segment (from D6S211 to D6S510) were associated with NPC susceptibility in consistence with previous findings [6–9, 15]; rs2517716 near HLA-A and rs9260734 near HCG9 were found significant, and a novel significant SNP, rs9404952, was close to the HLA-G gene suggesting HLA-G might be a new susceptible locus for NPC, and needs to be further researched.
The HLA-A gene is a HLA class I gene, which is involved in the presentation of foreign peptide antigens including EBV-derived peptides. HCG9 belongs to a multigene family (HCG) and is believed to be non-coding, but its function has not been determined. HLA-A and HCG9 were reportedly associated with the development of EBV-associated Hodgkin lymphoma  and infectious mononucleosis (IM) after primary EBV infection . The association of HLA-A and HCG9 polymorphisms with NPC was reported in numerous studies including two GWA studies [6–9]. In addition, a recent association study was carried out by Li Y  validating Taiwan's GWAS and presenting three NPC susceptible SNPs, of which two located at HLA-A and HCG9 loci. In the current study, we found several NPC-associated risk SNPs near HLA-A and HCG9 locus as well. Interestedly, rs2517716 as one of the most significant SNPs was closed to the HLA-A locus and genetically linked to rs2517713 reported significant in the Taiwan group(r2 = 0.7 and D = 1)based on Han Chinese (CHB) and Japanese (JPT) ancestries from HapMap data. Rs9260734 near HCG9 had been reported by both GWA studies in Taiwanese and Southern Chinese respectively [8, 9], and also confirmed in the present study though it seems that its association signal might originate from the effects of HLA-A according to a previous research .
Moreover, in addition of the HLA-A and HCG9, our data for the first time showed that HLA-G might be involved in NPC susceptibility in southern Chinese. HLA-G belongs to the HLA non-classical class I heavy chain paralogues and seems to play a crucial role in immuno-tolerance mechanisms . More recently, some studies have suggested that HLA-G plays an important role through the immune system in tumor development and progression [19–23]. It possesses the ability to suppress immune cell functions, such as NK- and CTL-mediated cytolysis and T-cell proliferative progress. HLA-G expression has been identified in several tumors such as breast cancer, glioblastoma, classical Hodgkin's lymphoma, and renal and lung cancers [19–23]. Ghandri N et al.  examined HLA-G gene polymorphisms in NPC for the first time using the restriction fragment length polymorphism-polymerase chain reaction (RFLP-PCR) and the amplification refractory mutation system-polymerase chain reaction (ARMS-PCR) method in Tunisian and found three nonsynonymous SNPs in codons not to associate with NPC susceptibility but implicate in NPC progression and survival. In the present study, rs9404952 located downstream of HLA-G gene showed a significant association with NPC though a less significant difference was found in the stage 1 (P = 0.017) possibly due to its small sample size. Hence, HLA-G might be a previously unreported candidate susceptibility locus for NPC among southern Chinese.
It is noted that previous GWA studies could not find the association of rs9404952 with NPC after scanning the whole human genome. The possible reasons are inferred as follow. First, due to its relatively larger P value (e.g., P > 10-5) in the preliminary stage of GWA studies, this SNP might be omitted from bioinformatics analysis. Second, such deviation might be resulted from genetic heterogeneity between populations. Furthermore, perhaps this SNP was not successfully genotyped by the genome-wide SNP-typing arrays due to some unclear reasons.
In this study, we utilized tag SNPs according to Hapmap database to target the 132 kb segment between D6S211 and D6S510, and identified several loci significantly associated with NPC. In contrast to the GWA studies that examined large sample sets and more than thousands of genetic markers, our case-control analysis employed less tag SNPs in relatively small samples. Interestedly, we could also detect significant polymorphisms in HLA-A and HCG9, and even identified a new candidate susceptibility locus near HLA-G. Our study supported that the some HLA loci indeed have strong genetic effect for NPC.
Therefore, further investigations should be done in the future to expand the search for NPC susceptibility loci within or close to the HLA region, to explore functional loci depended on more advanced means, such as targeted exomic sequencing, and to further investigate the association between HLA-G locus and NPC risk.
Curado MP, Edwards B, Shin HR, Storm H, Ferlay J, Heanue M, Boyle P: Cancer incidence in five continents. Lyon, IARC Sci Publ. 2007, Electronic publication, IX:
Hildesheim A, Wang C-P: Genetic predisposition factors and nasopharyngeal carcinoma risk: A review of epidemiological association studies, 2000-2011. Rosetta Stone for NPC: Genetics, viral infection, and other environmental factors. Semin Cancer Biol. 2012, doi:10.1016/j.semcancer.2012.01.007
Li X, Fasano R, Wang E, Yao KT, Marincola FM: HLA associations with nasopharyngeal carcinoma. Curr Mol Med. 2009, 9: 751-765. 10.2174/156652409788970698.
Tang M, Zeng Y, Poisson A, Marti D, Guan L, Zheng Y, Deng H, Liao J, Guo X, Sun S, Nelson G, de The G, Winkler CA, O'Brien SJ, Carrington M, Gao X: Haplotype-dependent HLA susceptibility to nasopharyngeal carcinoma in a Southern Chinese population. Genes Immun. 2010, 11 (4): 334-342. 10.1038/gene.2009.109.
Cho WC: Recent progress in genetic variants associated with cancer and their implications in diagnostics development. Expert Rev Mol Diagn. 2010, 10 (6): 699-703. 10.1586/erm.10.64.
Lu CC, Chen JC, Tsai ST, Jin YT, Tasi JC, Chan SH, Su IJ: Nasopharyngeal carcinoma-susceptibility locus is localized to a 132 kb segment containing HLA-A using high-resolution microsatellite mapping. Int J Cancer. 2005, 115: 742-746. 10.1002/ijc.20946.
Lu CC, Chen JC, Jin YT, Yang HB, Chan SH, Tsai ST: Genetic susceptibility to nasopharyngeal carcinoma within the HLA-A locus in Taiwanese. Int J Cancer. 2003, 103: 745-751. 10.1002/ijc.10861.
Tse KP, Su WH, Chang KP, Tsang NM, Yu CJ, Tang P, See LC, Hsueh C, Yang ML, Hao SP, Li HY, Wang MH, Liao LP, Chen LC, Lin SR, Jorgensen TJ, Chang YS, Shugart YY: Genome-wide association study reveals multiple nasopharyngeal carcinoma-associated loci within the HLA region at chromosome 6p21.3. Am J Hum Genet. 2009, 85: 194-203. 10.1016/j.ajhg.2009.07.007.
Bei JX, Li Y, Jia WH, Feng BJ, Zhou GQ, Chen LZ, Feng QS, Low HQ, Zhang HX, He FC, Tai ES, Kang TB, Liu ET, Liu JJ, Zeng YX: A genome-wide association study of nasopharyngeal carcinoma identifies three new susceptibility loci. Nat Genet. 2010, 42: 599-603. 10.1038/ng.601.
The International HapMap Consortium: The International HapMap Project. Nature. 2003, 426 (6968): 789-796. 10.1038/nature02168.
de Bakker PIW, McVean G, Sabeti PC, Miretti MM, Green T, Marchini J, Ke X, Monsuur AJ, Whittaker P, Delgado M, Morrison J, Richardson A, Walsh EC, Gao XJ, Galver L, Hart J, Hafler DA, Pericak-Vance M, Todd JA, Daly MJ, Trowsdale J, Wijmenga C, Vyse TJ, Beck S, Murray SS, Carrington M, Gregory S, Deloukas P, Rioux JD: A high-resolution HLA and SNP haplotype map for disease association studies in the extended human MHC. Nat Genet. 2006, 38 (10): 1166-1172. 10.1038/ng1885.
Bell PA, Chaturvedi S, Gelfand CA, Huang CY, Kochersperger M, Kopla R, Modica F, Pohl M, Varde S, Zhao R, Zhao X, Boyce-Jacino MT, Yassen A: SNPstream UHT: ultra-high throughput SNP genotyping for pharmacogenomics and drug discovery. Biotechniques. 2002, 70-2. 74, 76-7, Suppl
Li X, Yang XX, Hu NY, Sun JZ, Li FX, Li M: A risk-associated single nucleotide polymorphism of SMAD7 is common to colorectal, gastric, and lung cancers in a Han Chinese population. Mol Biol Rep. 2010, 38 (8): 5093-5097.
Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, Liu-Cordero SN, Rotimi C, Adeyemo A, Cooper R, Ward R, Lander ES, Daly MJ, Altshuler D: The structure of haplotype blocks in the human genome. Science. 2002, 296: 2225-2229. 10.1126/science.1069424.
Li Y, Fu L, Wong AM, Fan YH, Li MX, Bei JX, Jia WH, Zeng YX, Chan D, Cheung KM, Sham P, Chua D, Guan XY, Song YQ: Identification of genes with allelic imbalance on 6p associated with nasopharyngeal carcinoma in southern Chinese. PLoS One. 2011, 6 (1): e14562-10.1371/journal.pone.0014562.
Niens M, van den Berg A, Diepstra A, Nolte IM, van der Steege G, Gallagher A, Taylor GM, Jarrett RF, Poppema S, te Meerman GJ: The human leukocyte antigen class I region is associated with EBV-positive Hodgkin's lymphoma: HLA-A and HLA complex group 9 are putative candidate genes. Cancer Epidemiol Biomarkers Prev. 2006, 15: 2280-2284. 10.1158/1055-9965.EPI-06-0476.
McAulay KA, Higgins CD, Macsween KF, Lake A, Jarrett RF, Robertson FL, Williams H, Crawford DH: HLA class I polymorphisms are associated with development of infectious mononucleosis upon primary EBV infection. J Clin Invest. 2007, 117: 3042-3048. 10.1172/JCI32377.
Sheu J, Shih IM: HLA-G and immune evasion in cancer cells. J Formos Med Assoc. 2010, 109 (4): 248-257. 10.1016/S0929-6646(10)60050-2.
Chen HX, Lin A, Shen CJ, Zhen R, Chen BG, Zhang X: Up-regulation of human leukocyte antigen-G expression and its clinical significance in ductal breast cancer. Hum Immunol. 2010, 71: 892-898. 10.1016/j.humimm.2010.06.009.
Wiendl H, Mitsdoerffer M, Hofmeister V, Wischhusen J, Bornemann A, Meyermann R: A functional role of HLA-G expression in human gliomas: An alternative strategy of immune escape. J Immunol. 2002, 168: 4772-4780.
Diepstra A, Poppema S, Boot M, Visser L, Nolte IM, Niens M, te Meerman GJ, van den Berg A: HLA-G protein expression as a potential immune escape mechanism in classical Hodgkin's lymphoma. Tissue Antigens. 2008, 71: 219-6. 10.1111/j.1399-0039.2008.01005.x.
Ibrahim EC, Guerra N, Lacombe MJ, Angevin E, Chouaib S, Carosella ED, Caignard A, Paul P: Tumor specific up-regulation of the nonclassical class I HLA-G antigen expression in renal carcinoma. Cancer Res. 2001, 61: 6838-6845.
Urosevic M, Kurrer MO, Kamarashev J, Mueller B, Weder W, Burg G: Human leukocyte antigen G up-regulation in lung cancer associates with high-grade histology, human leukocyte antigen class I loss and interleukin-10 production. Am J Pathol. 2001, 159: 817-824. 10.1016/S0002-9440(10)61756-7.
Ghandri N, Gabbouj S, Farhat K, Bouaouina N, Abdelaziz H, Nouri A, Chouchane L, Hassen E: Association of HLA-G polymorphisms with nasopharyngeal carcinoma risk and clinical outcome. Hum Immunol. 2011, 72: 150-158. 10.1016/j.humimm.2010.10.006.
This project was funded by grants from National Natural Science Foundation of China (No. 30371535, No. 30973292).
The authors declare that they have no competing interests.
XL, KY and ML designed the study and modified manuscript as well as participated in data analysis and interpretation. MZ, XL and HC carried out most experimental studies and drafted the manuscript. WF and XY were in charge of the data assembly and statistic analysis. HZ collected blood sample materials and clinical data. GZ and GW mainly established and maintained sample database including clinical data. All authors read and approved the final manuscript.
Manli Zhao, Hongbing Cai, Xin Li, Hang Zheng, Xuexi Yang contributed equally to this work.