PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Genes Immun. Author manuscript; available in PMC 2013 August 8.
Published in final edited form as:
PMCID: PMC3737777
NIHMSID: NIHMS493525

Haplotype-dependent HLA susceptibility to nasopharyngeal carcinoma in a Southern Chinese population

Abstract

We have conducted a comprehensive case–control study of a nasopharyngeal carcinoma (NPC) population cohort from Guangxi Province of Southern China, a region with one of the highest NPC incidences on record. A total of 1407 individuals including NPC patients, healthy controls, and their adult children were examined for the human leukocyte antigen (HLA) association, which is so far the largest NPC cohort reported for such studies. Stratified analysis performed in this study clearly demonstrated that while NPC protection is associated with independent HLA alleles, most NPC susceptibility is strictly associated with HLA haplotypes. Our study also detected for the first time that A*0206, a unique A2 subtype to South and Southeast Asia is also associated with a high risk for NPC. HLA-A*0206, HLA-B*3802 alleles plus the A*0207–B*4601 and A*3303–B*5801 haplotypes conferred high risk for NPC showing a combined odds ratio (OR) of 2.6 (P<0.0001). HLA alleles that associate with low risk for NPC include HLA-A*1101, B*27, and B*55 with a combined OR of 0.42 (P<0.0001). The overall high frequency of NPC-susceptible HLA factors in the Guangxi population is likely to have contributed to the high-NPC incidence in this region.

Keywords: HLA, nasopharyngeal carcinoma, haplotype, stratified analysis

Introduction

Nasopharyngeal carcinoma (NPC) is an epithelial malignancy caused by a combination of Epstein–Barr virus (EBV) infection and environmental factors.13 Genetic predisposition also has a function in NPC susceptibility causing the observed familial aggregation4 and distinct racial and geographical distribution of the disease incidence.58 Southern China and Southeast Asia have a disproportionally high incidence of NPC. In these regions, NPC is one of the most common cancers,9,10 but in Caucasians, it is a rare disease.11

Among the host genetic markers that have been associated with NPC, the class I human leukocyte antigen (HLA) genes have shown a strong and consistent association with disease risk. As it was first reported by Simons et al.12 in a Singaporean Chinese cohort in the mid 1970s, the HLA association with NPC has been widely detected in different racial groups even though the exact HLA factors (alleles and haplotypes) associated with the disease sometimes vary among racial groups due to population-dependent HLA distributions. Populations of Southern China and Southeast Asia share high-NPC incidence as well as common HLA characteristics. Mainland Chinese,9,13 Taiwanese,14 and Singaporean Chinese15 NPC cohorts also show similar HLA associations, where HLA-A*11 and B*13 seem protective against NPC, and A*02 (A*0207), A*33, B*46, and B*58 associate with susceptibility to this disease (Figure 1). Other HLA loci including HLA-C and the class II loci, DR and DQ, have not shown independent associations with NPC even though certain class II alleles might be part of the extended HLA haplotypes associated with the disease as described by Hildesheim et al.14

Figure 1
HLA associations with NPC in Chinese population.

Association of HLA polymorphism with human disease may indicate direct involvement of the HLA molecule in the disease pathogenesis. To date, however, there has been little in vitro or in vivo evidence that the NPC-associated HLA alleles affect differentially NPC-related EBV replication on pathogenesis. The role of HLA-A and -B alleles or haplotypes that include combinations of these alleles in NPC pathogenesis remains unknown. It is possible that the associated HLA class I alleles are simply marking by linkage disequilibrium (LD) the true NPC-causing gene(s). Large-scale cohort studies exploring combinations of associated HLA alleles may provide useful insights into the influence of HLA on NPC.

An early indicator of NPC development is the occurrence of immunoglobulin (Ig) A antibodies to EBV capsid antigens (EBV-IgA/VCA).1618 Even though >95% of adults in the general population of all ethnic groups are healthy carriers of EBV, <2.5% are EBV-IgA/VCA antibody positive. In comparison, >95% of all NPC patients are EBV-IgA/VCA antibody positive.10 If HLA diversity is indeed directly responsible for the individually varied NPC risks, it is plausible that the development of the EBV-IgA/VCA antibodies in EBV-positive individuals may also be affected by the HLA polymorphism.

Here, we have conducted a case–control study of an NPC cohort recruited from Guangxi Province in Southern China where the NPC incidence is as high as 25–50 cases per 10000 individuals. DNA-based high-resolution HLA typing was performed on a total of 1407 individuals including NPC cases, matched controls, and offspring of the study subjects. This large study has allowed a comprehensive stratification of NPC-associated HLA factors and provided novel insights into the nature of the HLA association with the disease.

Results

HLA typing was informative for 356 NPC patients, 287 NPC free EBV-IgA/VCA antibody positive healthy individuals, and 342 NPC free EBV-IgA/VCA antibody negative healthy individuals. Comparative analyses between the two healthy groups failed to detect any significant deviation in the frequency distribution of HLA alleles and haplotypes (Supplementary Tables 1 and 2), indicating that HLA polymorphism does not affect the occurrence of the EBV-IgA/VCA antibody. Therefore, in subsequent analyses, the two NPC-free groups were combined as the control group (N = 629) for NPC cases.

HLA alleles showing a significant difference in frequency distribution between cases and controls are listed in Table 1, and full analyses are presented in Supplementary Tables 3 and 4. For the reason of most of former NPC, HLA studies were based on serology typing, for better understanding the HLA influence in NPC, both HLA allotype and genotype were present in this table. For the HLA-A locus, 31 four-digit alleles were detected, 12 of which had an allele frequency >1% in either the case or control group. Five alleles, A*0206, A*0207, A*1101, A*3303, and A*7401/7402, showed a significant difference in frequency distribution between the case and control groups. After correction, however, only two alleles, A*1101 and A*3303, remained significant. Of the two detected A*11 alleles, A*1101 showed a reduced presence in patients compared with controls (P<0.0001), whereas the less common A*1102 showed no difference between these groups, but due to its low frequency, it need to be confirmed in even large study cohort. Four common A*02 subtypes, A*0201, A*0203, A*0206, and A*0207, were detected. Two of alleles, A*0206 and A*0207, showed an elevated frequency in patients. A*3303 also associated with increased risk of NPC (P = 0.0004).

Table 1
Gene frequencies (%) of the two and four-digit HLA-A, -B, and -C alleles detected in NPC patients (N = 356) and controls (N = 629)

Fifty-five HLA-B alleles were detected, but only 14 had a frequency >1%. Among the 14 B alleles, seven showed a significantly different distribution between cases and controls, including B*1301, B*2704, B*3802, B*4001, B*5502, B*5601, and B*5801. After P-value correction, however, only the decrease of B*5502 in the patient group remained significant (P corrected = 0.003). Three B*27 subtypes, B*2704, B*2705, and B*2706 were detected only in the control group with a combined frequency of 1.59% (P = 0.0014). B*3802, B*4001, and B*5801 were each observed more frequently among the case group but the significance disappeared after correction.

Twenty-four alleles were detected at the HLA-C locus, 13 of which had a frequency >1%. Cw*0403, Cw*1202, and Cw*1203 were observed at a lower frequency in cases and Cw*0302 was the only allele showing an elevated frequency in cases. After correction, however, none of the four-digit HLA-C alleles remained significant.

HLA typing was informative for 422 children of the study subjects, which enabled us to directly determine HLA haplotypes in 179 patients and 379 controls. For the remaining 177 patients and 250 controls, HLA haplotypes were assigned by population-based estimation methods. HLA-A/B haplotype frequencies were calculated on the basis of both methods of haplotype assignments. Stratification analyses were performed to determine the nature of the observed HLA association (Tables 2 and and3).3). Table 2 compares the frequency distribution of all HLA-A/B (where the HLA-A allele showed individual association with NPC) and -B/A (where the HLA-B allele showed individual association with NPC) haplotypes related to individual NPC-associated alleles between cases and controls to determine whether the observed allele association with NPC might be due to particular haplotypes or to an individual allele. Four of the five potentially NPC-related HLA-A alleles, A*0206, A*0207, A*1101, and A*3303, were each observed on HLA-A/B haplotypes with frequencies of >1% and were therefore included in Table 2A. Both the dominant A*0206–B*1502 haplotype (P = 0.0278) and all other A*0206-associated haplotypes combined (P = 0.0071) showed elevated frequencies in the patient group. A*0207 is predominantly associated with B*4601, and this was the only A*0207 haplotype associating with susceptibility to NPC (P = 0.0095). Similarly, A*3303–B*5801 was the only A*3303 haplotype associating with NPC susceptibility (P = 0.0003). The protective A*1101 allele was found on seven HLA-A/B haplotypes with frequencies >1% and these haplotypes showed inconsistent associations with NPC. A*1101–B*1301 (P = 0.0029), A*1101–B*4601 (P = 0.0198), and A*1101–B*others (P = 0.0009) showed reduced frequencies in the NPC group, whereas A*1101–B*3802 haplotype (P = 0.0228) showed an elevated frequency in the patient group.

Table 2
HLA-A/B and B/A haplotype frequencies (%) in NPC patients (N = 356) and controls (N = 629)a
Table 3
Effect of allelic combinatioNS of NPC-associated HLA alleles on the disease development

Five of the seven NPC-associated HLA-B alleles composed haplotypes with frequencies >1%. The B*5801–A*3303 haplotype showed an elevated frequency in NPC patients, whereas the other B*5801 haplotypes did not. The other NPC-associated B alleles were all found on at least two common HLA-B/A haplotypes. Both B*1301 and B*3802 were significantly more common in the control group only when A*1101 was also present on these haplotypes (P = 0.0029 and P = 0.0228, respectively). B*4001 was associated with four haplotypes with a frequency >1%. All haplotypes containing this allele were observed more frequently in the patients (B*4001–A*0203 and B*4001–A*1102 reaching significance; P = 0.0167 and P = 0.0138, respectively), except for the A*1101-associated B*4001 haplotype, which was more common in the controls albeit not significantly.

Table 3 stratifies six common NPC-associated HLA-A/B allelic combinations to examine independent as well as collective effects of HLA alleles on NPC risk analyzed in Table 2. For each tested HLA-A/B phenotypic combination, the study subjects were divided into three groups: those having both of the relevant HLA-A and HLA-B alleles and those that had one, but not the other. The protective A*1101 and the susceptible B*4001 allele combination is included first in Table 3 to show the collective effect of these two offsetting NPC-associated alleles. The protective B*5502 is in LD with A*0203, and stratification for this allele combination showed that the protective effect of B*5502 could be detected in the presence (odds ratio (OR) = 0.18, P = 0.0018) or absence (OR = 0.11, P = 0.0088) of A*0203. A*1101, which is commonly linked to the NPC-related B alleles B*1301 and B*4001, is the most common protective allele in the Guangxi NPC cohort. Stratified analysis showed that the A*1101 protection was stronger when B*1301 was also present (OR = 0.5, P = 0.0010; without B*1301; OR = 0.73, P = 0.0280). However, B*1301 had no effect (OR = 0.9, P>0.05) in the absence of A*1101. Remarkably, B*4001 actually associated strongly with risk of developing NPC if A*1101 was missing (OR = 2.56, P<0.0001), but in the presence of A*1101, the susceptibility effect of B*4001 was completely abrogated (OR = 0.66, P>0.05). The susceptibility effect of A*0206 seems to be largely independent of HLA-B, though the A*0206 and B*1502 combination was somewhat stronger (OR = 3.63, P = 0.0065) than that in the presence of other HLA-B alleles (OR = 1.99, P = 0.0184). Both A*0207 (Table 1) and B*460114,19 were thought to be independent high-risk factors for NPC in Southern Chinese and Southeast Asia. The two alleles are in strong LD in the Guangxi cohort, as in most Southern Chinese populations, and the A*0207–B*4601 haplotype associated with susceptibility (Table 3). However, A*0207 showed no effect in the absence of B*4601, whereas B*4601 had a weak protective effect in the absence of A*0207 (OR = 0.61, P = 0.0325). Finally, A*3303 and B*5801 forms the most common A–B haplotype in the Guangxi cohort and both alleles were associated with an elevated NPC risk (Table 1). Stratification of the A*3303 and B*5801 combination showed that the NPC effect was more strongly associated with the presence of both alleles (OR = 1.79, P = 0.0001) as one allele without the other had no effect.

Table 4 summarizes the HLA alleles and the allele combinations associated with the risk of developing NPC. HLA protection against NPC mainly involves independent allele influence except that the A*1101 protection may be overridden by the presence of the high-risk B*3802. The presence of one or more of the six protective factors has a frequency of 58.98% within the control group and 34.83% within the NPC patients. Together, these protective factors delivered a combined OR of 0.37 (P<0.0001). NPC susceptibility, on the other hand, associates most strongly with certain allele combinations, as opposed to single alleles, with the exception of the A*0206 effect. Individual positive for any of the five high-risk factors accounted for 50.08 and 73.03% of the controls and cases, respectively. These high-risk factors showed individual ORs between 2.05 and 3.56 and a collective OR of 2.70 (P<0.0001).

Table 4
Protective and susceptible HLA alleles, haplotypes, and genotypes for the development of NPC

On the basis of analyses shown, HLA haplotypes were classified as NPC protective (P), susceptible (S), and neutral (N), generating six genotypes. Figure 2 shows the ORs of the six genotypes for NPC development. Interestingly, neither the protective nor the susceptible genotypes seemed dominant over the other. The exact values are presented in Supplementary Table 5. The ORs of the remaining genotypes were distributed in a manner that was expected given that the protective and susceptible genotypes were defined in this same cohort. The validity of this scheme will be particularly interesting to test in an independent cohort.

Figure 2
ORs of HLA genotypes for NPC development.

The influence of HLA genotypes on NPC risk was analyzed by comparing all individual genotypes in HLA-A and -B loci separately and the compound genotypes of HLA-A and -B. Two HLA-A, six HLA-B genotypes and six HLA-A-B compound genotypes were significantly associated with either elevated NPC risk or protection (Supplementary Table 6). These genotypes all involve at least one allele showing a higher risk or protection for the disease.

Discussion

Southern China has a disproportionally high incidence of NPC compared with other parts of the world. The study population from Wuzhou City of Guangxi Province in Southern China holds, perhaps, the highest recorded NPC incidence.9,10 Apart from unique environmental factors (mainly traditional diet), our data support a role for genetic predisposition contributing to the high-disease incidence. The results from this study confirm and extend previously reported HLA and NPC associations in Southern Chinese populations,1214,19 such as the protective effect associated with HLA-A*1101 and susceptibility effects associated with B*0207, A*3303, and B*5801. This study provides further insights into the nature of the HLA association with NPC, such as the observation that the dominant A*11 protection can be entirely attributed to the major subtype A*1101. The other A*11 subtype detected in this population, A*1102, did not show any protective effect even though the two A*11 subtypes differ by only a single amino acid in position 19 outside of the peptide-binding groove.

Our stratified analyses demonstrated that a proportion of the HLA-associated NPC susceptibility is likely to be haplotype-dependent in this cohort (Table 3). In earlier cohort studies from Southern China13 and Taiwan,14 the A*0207–B*4601 allele combination was consistently associated with a high-NPC risk, but there are contradictory reports on whether the two alleles have independent effects. In our study cohort, two major high-risk HLA-A/B allele combinations, A*0207–B*4601 and A*3303–B*5801, were identified. Patients with either of these two allele combinations accounted for half of the NPC cohort and susceptibility conferred by these HLA-A/B combinations was strictly dependent on the presence of both alleles. The alleles comprising both pairs are in strong LD, forming the two most common HLA-A/B haplotypes in the Guangxi population.

Haplotype-dependent disease associations may indicate one of the two possibilities: (1) the two alleles on the haplotype are behaving in an epistatic manner to reach functional synergy or (2) the disease-associated HLA haplotype is tracking an unidentified disease locus present on that specific haplotype. An example of the former is that the structural stability of the HLA-DQ αβ dimer is affected by the type of DQA and DQB alleles carried on HLA haplotypes, which in turn influences the risk of insulin-dependent diabetes mellitus (IDDM).20 Functional synergy between HLA-A and -B molecules and any potential clinical relevance such synergy may have, however, are yet to be established. In terms of the alternative HLA marker model, there are no data that unequivocally support or refute the hypothesis in NPC studies, despite efforts to identify non-HLA disease genes in the extended MHC region.21,22 Still, direct involvement of HLA molecules in NPC pathogenesis, including that involving the causative EBV infection, also lacks functional data support. Interpretations of the present results tend to support the locus tracking hypothesis in that there is evidence for (1) a lack of HLA association with the occurrence of IgA antibodies to EBV-IgA/VCA, a well-known precursor of NPC, (2) haplotype-dependent but allele-independent disease susceptibility, and (3) the lack of association of A*1102 despite its identical peptide-binding structure shared with the well-documented protective A*1101 allele.14 Furthermore, distinct population-specific HLA alleles associate with susceptibility to NPC, which is unexpected if HLA is directly involved in disease pathogenesis.

In contrast to the haplotype-dependent but allele-independent susceptibility exhibited by the two major high-risk haplotypes A*0207–B*4601 and A*3303–B*5801, the dominant protective effect conferred by A*1101 did not show a clear haplotypic dependence. In fact, most A*1101-linked A/B haplotypes showed a trend of under-representation in the case group except the one linked with the susceptible B*3802. A*1101 seems to protect against other viral infections as well. For example, A*1101 may confer protection against AIDS through restriction of immunodominant HIV peptides,23,24 but the structural and functional basis for its direct involvement in NPC pathogenesis has not been thoroughly explored.

The detected allele and haplotype associations with NPC seem to be independent of HLA genotypes as genotypes showing significant ORs are essentially those with least two copies of NPC-associated alleles or haplotypes. Also, no particular genotypes showed stronger association than individual alleles or haplotypes.

Another noteworthy HLA association observed in the Guangxi NPC cohort is an opportune ‘complete’ protection of B*27 against disease development. B*27 was detected with a typical frequency in the control group (1.59%), but was completely absent in the patient group (OR = 0.04, P = 0.001). Earlier cohort studies from Southern China and Taiwan have also detected B*27 as a low-risk allele though the protection was never ‘complete.’14,19,25 B*27 is well known for its high-risk association with the inflammatory autoimmune disease ankylosing spondylitis and it confers protection against AIDS progression apparently through its peptide-binding properties for immunodominant HIV epitopes.26 Perhaps B*27 is also involved in EBV-related pathogenesis, impacting the risk of NPC development. Given the low frequency of B*27 in the study population, however, the absence of B*27 in the patient group should be confirmed in larger replication cohort studies.

The overall HLA profile of a population is likely to influence the incidence of HLA-associated diseases. An example is the correlation between the rates of IDDM incidence and the frequency distribution of the IDDM-associated HLA-DQ and DR alleles in world populations.27 The NPC incidence in a population may also correlate positively with the frequencies of NPC-susceptible HLA alleles and haplotypes of the population. The two dominant susceptible HLA haplotypes, A*0207–B*4601 and A*3303–B*5801, are the two most common HLA-A/B haplotypes in Guangxi and together they account for 20% of the total haplotype frequency of the study Wuzhou population, which is one of the highest in Far East Asian populations. Overall, the susceptible HLA factors were found in about half of the individuals in the Wuzhou population, which is also one of the highest. Therefore, the unusually high presence of NPC-associated HLA alleles and haplotypes may partly explain the disproportionally high-NPC incidence in Southern China and Southeast Asia.

Materials and methods

Study cohorts

NPC cases and controls were recruited from Wuzhou City and Cangwu County of Guangxi Province.10 All study subjects were of Han ethnic origin. Informed consent was obtained from all study participants. An effort was made to enroll triads consisting of a proband (either NPC patient or NPC-free but EBV-IgA/VCA antibody positive), an unaffected spouse and an adult child. The case group included 356 unrelated patients with biopsy-confirmed NPC. The mean age was 50.1 years (range 19–80), 95.5% of them being EBV-IgA/VCA antibody positive. Two groups of control were the case’s spouse or geographically matched residents who were NPC free at the time of study enrollment. An antibody to EBV capsid antigen (EBV-IgA/VCA) were confirmed by serologic testing at the time of study enrollment. One group was positive (N = 287) and the other negative (N = 342) for the EBV-IgA/VCA antibody. The mean age was 45.7 and 46.9, respectively, for the antibody positive and negative groups. The controls were matched to the cases by age, ethnicity, and geographic residence (Table 5). In addition, 422 adult children of the study subjects in the case and control groups were recruited to allow elucidation of HLA haplotypes, but they were excluded in all other analyses.

Table 5
Gender and age of NPC patients and EBV-IgA/VCA positive and negative controls

Detection of EBV-IgA/VCA antibody

The presence of the EBV-IgA/VCA antibody was detected using the immunoperoxidase assay as described earlier.28 EBV positive B95-8 cells fixed on slides were incubated with multiple dilutions of the testing serum followed by incubation with antihuman IgA horseradish peroxidase and staining with diaminobenzidine. Testing sera with a staining titer of 1:10 or higher dilution were considered positive for the EBV-IgA/VCA antibody.

HLA typing

HLA class I alleles were characterized using a PCR-SSOP (sequence-specific oligonucleotide probe) typing protocol developed by the 13th International histocompatibility Workshop.29 Briefly, the gene fragment spanning exon 2, intron 2, and exon 3 was amplified using locus-specific primers for HLA-A, -B, and -C separately. The PCR products were immobilized on nylon membranes and hybridized with a panel of P32-labeled oligonucleotide (19 mers) matching all known sequence variations of the HLA genes. Typing results were interpreted by SSOP hybridization patterns based on sequences of known HLA alleles. Typing ambiguities were resolved by sequencing exons 2 and 3 completely. For sequencing analysis, the PCR product of HLA-A, -B, or -C was used as the template for the sequencing reaction. For each of the HLA genes, two sequencing reactions (one for exon 2 and one for exon 3) were performed using exon-specific sequencing primers. The sequencing analysis was performed using the ABI Big Dye Terminator Cycle Sequencing Kit and ABI3730xl DNA analyzer (Applied Biosystems, Foster City, CA, USA). HLA alleles were assigned on the basis of the sequence database of known alleles with the help of the ASSIGN software developed by Conexio Genomics (Conexio Genomics, Perth City, WA, Australia). Ambiguous heterozygous genotypes were resolved by additional PCR and sequencing procedures using allele-specific PCR primers to selectively amplify only one of the two alleles.

Statistical analyses

HLA allele frequencies were calculated based on observed genotypes, and HLA-A and -B haplotype frequencies were estimated using one of the two methods: (1) unambiguous assignment based on familial segregation for the cases and controls using recruited spouse and child genotype data or (2) indirect assignment based on maximum likelihood estimation for the study subjects without recruited family members. For the latter, the haplotypic analysis was performed using the BLOCKHEAD genetic analysis software developed by George Nelson in the Laboratory of Genomic Diversity, National Cancer Institute.

The effect of HLA alleles on the development of NPC and EBV-IgA/VCA antibody was evaluated by computing ORs and 95% confidence intervals as well as exact P-values using the FREQ procedure of the SAS 9.1 software (The SAS Institute, NC, USA). A correction of 0.5 was applied on every cell of the 2 × 2 table that contains a zero. The analyses were performed at four-digit and two-digit resolution levels separately. P-value was calculated by χ2-test for each allele and was corrected by multiplying the number of all detected alleles. Significance was considered at P<0.05. Both uncorrected and corrected P-values were presented in our tables. Stratified analyses were applied to evaluate the effect of haplotypes of HLA-A and -B that contained alleles showing individual significant associations with NPC risk.

Supplementary Material

table 1-6

Acknowledgments

This project has been funded in whole or in part with federal funds from the National Cancer Institute, National Institutes of Health, under Contract No. HHSN261200800001E. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the US Government. This Research was supported in part by the Intramural Research Program of the NIH, National Cancer Institute, Center for Cancer Research.

Footnotes

Conflict of interest

The authors declare no conflict of interest

Supplementary Information accompanies the paper on Genes and Immunity website (http://www.nature.com/gene)

References

1. Brown TM, Heath CW, Lang RM, Lee SK, Whalley BW. Nasopharyngeal cancer in Bermuda. Cancer. 1976;37:1464–1468. [PubMed]
2. Yu MC, Garabrant DH, Huang TB, Henderson BE. Occupational and other non-dietary risk factors for nasopharyngeal carcinoma in Guangzhou, China. Int J Cancer. 1990;45:1033–1039. [PubMed]
3. Guo X, Johnson RC, Deng H, Liao J, Guan L, Nelson GW, et al. Evaluation of non-viral risk factors for nasopharyngeal carcinoma in a high-risk population of Southern China. Int J Cancer. 2009;124:2942–2947. [PubMed]
4. Coffin CM, Rich SS, Dehner LP. Familial aggregation of nasopharyngeal carcinoma and other malignancies. A clinicopathologic description. Cancer. 1991;68:1323–1328. [PubMed]
5. Jeannel D, Hubert A, de Vathaire F, Ellouz R, Camoun M, Ben Salem M, et al. Diet, living conditions and nasopharyngeal carcinoma in Tunisia–a case-control study. Int J Cancer. 1990;46:421–425. [PubMed]
6. Laramore GE, Clubb B, Quick C, Amer MH, Ali M, Greer W, et al. Nasopharyngeal carcinoma in Saudi Arabia: a retrospective study of 166 cases treated with curative intent. Int J Radiat Oncol Biol Phys. 1988;15:1119–1127. [PubMed]
7. Johansen LV, Mestre M, Overgaard J. Carcinoma of the nasopharynx: analysis of treatment results in 167 consecutively admitted patients. Head Neck. 1992;14:200–207. [PubMed]
8. Lee AW, Foo W, Mang O, Sze WM, Chappell R, Lau WH, et al. Changing epidemiology of nasopharyngeal carcinoma in Hong Kong over a 20-year period (1980–99): an encouraging reduction in both incidence and mortality. Int J Cancer. 2003;103:680–685. [PubMed]
9. Lu SJ, Day NE, Degos L, Lepage V, Wang PC, Chan SH, et al. Linkage of a nasopharyngeal carcinoma susceptibility locus to the HLA region. Nature. 1990;346:470–471. [PubMed]
10. Guo XC, Scott K, Liu Y, Dean M, David V, Nelson GW, et al. Genetic factors leading to chronic Epstein-Barr virus infection and nasopharyngeal carcinoma in South East China: study design, methods and feasibility. Hum Genomics. 2006;2:365–375. [PMC free article] [PubMed]
11. Burt RD, Vaughan TL, McKnight B, Davis S, Beckmann AM, Smith AG, et al. Associations between human leukocyte antigen type and nasopharyngeal carcinoma in Caucasians in the United States. Cancer Epidemiol Biomarkers Prev. 1996;5:879–887. [PubMed]
12. Simons MJ, Wee GB, Chan SH, Shanmugaratnam K, Day NE, de-The G. Immunogenetic aspects of nasopharyngeal carcinoma (NPC) III. HL-a type as a genetic marker of NPC predisposition to test the hypothesis that Epstein-Barr virus is an etiological factor in NPC. IARC Sci Publ. 1975:249–258. [PubMed]
13. Hu SP, Day NE, Li DR, Luben RN, Cai KL, Ou-Yang T, et al. Further evidence for an HLA-related recessive mutation in nasopharyngeal carcinoma among the Chinese. Br J Cancer. 2005;92:967–970. [PMC free article] [PubMed]
14. Hildesheim A, Apple RJ, Chen CJ, Wang SS, Cheng YJ, Klitz W, et al. Association of HLA class I and II alleles and extended haplotypes with nasopharyngeal carcinoma in Taiwan. J Natl Cancer Inst. 2002;94:1780–1789. [PubMed]
15. Chan SH, Day NE, Kunaratnam N, Chia KB, Simons MJ. HLA and nasopharyngeal carcinoma in Chinese—a further study. Int J Cancer. 1983;32:171–176. [PubMed]
16. Henle G, Henle W. Epstein-Barr virus-specific IgA serum antibodies as an outstanding feature of nasopharyngeal carcinoma. Int J Cancer. 1976;17:1–7. [PubMed]
17. Zeng Y, Zhang LG, Li HY, Jan MG, Zhang Q, Wu YC, et al. Serological mass survey for early detection of nasopharyngeal carcinoma in Wuzhou City, China. Int J Cancer. 1982;29:139–141. [PubMed]
18. Zeng Y, Zhang LG, Wu YC, Huang YS, Huang NQ, Li JY, et al. Prospective studies on nasopharyngeal carcinoma in Epstein-Barr virus IgA/VCA antibody-positive persons in Wuzhou City, China. Int J Cancer. 1985;36:545–547. [PubMed]
19. Wu SB, Hwang SJ, Chang AS, Hsieh T, Hsu MM, Hsieh RP, et al. Human leukocyte antigen (HLA) frequency among patients with nasopharyngeal carcinoma in Taiwan. Anticancer Res. 1989;9:1649–1653. [PubMed]
20. Kwok WW, Schwarz D, Nepom BS, Hock RA, Thurtle PS, Nepom GT. HLA-DQ molecules form alpha-beta heterodimers of mixed allotype. J Immunol. 1988;141:3123–3127. [PubMed]
21. Lu CC, Chen JC, Tsai ST, Jin YT, Tsai JC, Chan SH, et al. 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. [PubMed]
22. Ooi EE, Ren EC, Chan SH. Association between microsatellites within the human MHC and nasopharyngeal carcinoma. Int J Cancer. 1997;74:229–232. [PubMed]
23. Beyrer C, Artenstein AW, Rugpao S, Stephens H, VanCott TC, Robb ML, et al. Epidemiologic and biologic characterization of a cohort of human immunodeficiency virus type 1 highly exposed, persistently seronegative female sex workers in northern Thailand. Chiang Mai HEPS Working Group. J Infect Dis. 1999;179:59–67. [PubMed]
24. Fukada K, Tomiyama H, Chujoh Y, Miwa K, Kaneko Y, Oka S, et al. HLA-A*1101-restricted cytotoxic T lymphocyte recognition for a novel epitope derived from the HIV-1 Env protein. Aids. 1999;13:2597–2599. [PubMed]
25. Yu MC, Huang TB, Henderson BE. Diet and nasopharyngeal carcinoma: a case-control study in Guangzhou, China. Int J Cancer. 1989;43:1077–1082. [PubMed]
26. Gao X, Bashirova A, Iversen AK, Phair J, Goedert JJ, Buchbinder S, et al. AIDS restriction HLA allotypes target distinct intervals of HIV-1 pathogenesis. Nat Med. 2005;11:1290–1292. [PubMed]
27. Bao MZ, Wang JX, Dorman JS, Trucco M. HLA-DQ beta non-ASP-57 allele and incidence of diabetes in China and the USA. Lancet. 1989;2:497–498. [PubMed]
28. Cevenini R, Donati M, Rumpianesi F, Moroni A, Paolucci P. An immunoperoxidase assay for the detection of specific IgA antibody in Epstein-Barr virus infections. J Clin Pathol. 1984;37:440–443. [PMC free article] [PubMed]
29. Gao X, Nelson GW, Karacki P, Martin MP, Phair J, Kaslow R, et al. Effect of a single amino acid change in MHC class I molecules on the rate of progression to AIDS. N Engl J Med. 2001;344:1668–1675. [PubMed]