Associations between HLA Class I alleles and the prevalence of nasopharyngeal carcinoma (NPC) among Tunisians
Journal of Translational Medicine volume 5, Article number: 22 (2007)
The high prevalence of nasopharyngeal cancer (NPC) in Southern Asia and Mediterranean Northern Africa suggests genetic predisposition among other factors. While Human Leukocyte Antigen (HLA) haplotypes have been conclusively associated with NPC predisposition in Asians, Northern African Maghrebians have been less intensely studied. However, low resolution serological methods identified weak positive associations with HLA-B5, B13 and B18 and a negative with HLA-B14. Using sequence based typing (SBT), we performed a direct comparison of HLA class I frequencies in a cohort of 136 Tunisian patients with NPC matched for gender, age and geographical residence to 148 normal Tunisians. The bimodal age distribution of NPC in Maghrebians was also taken into account. HLA frequencies in normal Tunisians were also compared with those of Northern Moroccan Berbers (ME) to evaluate whether the Tunisian population in this study could be considered representative of other Maghrebian populations. HLA-B14 and -Cw08 were negatively associated with NPC (odd ratio = 0.09 and 0.18 respectively, Fisher p2-value = 0.0001 and = 0.003). Moreover, positive associations were observed for HLA-B-18, -B51 (split of -B5) and -B57 (p2-value < 0.025 in all) confirming previous findings in Maghrebs. The HLA-B14/Cw*08 haplotype frequency (HF) was 0.007 in NPC patients compared to 0.057 in both Tunisian (OR = 0.12; p2-value = 0.001) and Moroccan controls. This study confirms several previous associations noted by serologic typing between HLA class I alleles and the prevalence of NPC in Maghrebians populations. In addition, we identified a putative haplotype rare in Tunisian patients with NPC that may serve as a genetic marker for further susceptibility studies.
The prevalence of Epstein-Barr virus (EBV)-associated anaplastic nasopharyngeal carcinoma (NPC) might be partially dependent upon genetic background affecting predominantly certain ethnic groups1. Among them, some regions of Southern Asia suffer about a 100-fold higher prevalence compared with populations not at risk and this high prevalence is maintained upon immigration to Western Countries such as in the U.S. Chinese population Northern African populations (Algerians, Moroccans and Tunisians) referred to as Maghrebians also suffer increased prevalence of NPC which persists upon migration to low risk areas [1, 2] Although heterogeneity of viral strains, dietary, environmental and socio-economical variables have been implicated as predisposing factors, it is likely that genetic traits play a significant role at least in Asians [3, 4].
Genetic predisposition of NPC is paralleled by a strong associations with some HLA class I alleles which have been extensively studied in the Asian population. Of them, HLA-A2, -B14, -B46 and -B58 predispose to the disease while HLA-A11 exercises a protective effect in Asians [5–12] Moreover, a meta-analysis confirmed a consistent association between HLA-A2, A11, B14, B46 and NPC . With few exceptions [10, 11]., most studies employed serologic HLA typing methods that cannot discriminate alleles within large serologic families . Such discrimination, however, is critical since different ethnic groups display significant differences in frequency of HLA alleles within a serologic family with different capacity to bind and present epitopes that may be relevant to the immune defense against EBV-coded or other NPC-associated antigens. For instance, the HLA-A2 serologic family includes about 100 alleles  whose distribution is different among ethnic groups with HLA-A*0201 (~90%) and HLA-A*0205 (~5%) predominant among Caucasian and HLA-A*0203 (~20%) -A*0206 (~15%) and -A*0207 (~40%) in Asians [16–19] Interestingly, sequence based typing (SBT) in Taiwanese identified a positive association between HLA-A*0207 (but not other HLA-A2 subtypes) and NPC . This association was particularly strong when the observation was extended to the HLA-A*0207/B*4601 haplotype. Since, chromosome 6p, in particular 6p21-23, is characterized by unusually high gene density harboring gene clusters coding for several protein families besides HLA molecules , it remains unclear whether NPC/HLA associations reflect distinct antigen presentation potentials of various alleles or rather represent a marker for a susceptibility/protection locus in strong linkage disequilibrium with the HLA region .
Maghrebian populations have been less intensely studied at high-resolution SBT. Serological analyses reported positive associations between NPC and HLA: these included increased frequency of HLA-B5 in Algerians (38.2% vs 24.4%) , HLA-B13 in Tunisians (15.5% vs 4%)  and HLA-B18 in Moroccans, of which only the last was significant after correction for number of specificities tested (relative risk = 4.14) . In 1983, Heirat P et al.  reported in a cohort of patients with NPC from Algeria a lower frequency of HLA-Aw33 (3.9% vs 16.8%), -B14 (1.3% vs 16%) and -DR4 (13.2% vs 29.1%) in NPC compared with matched controls. After correction for the number of specificities tested, these differences were not statistically significant and were not pursued further. Subsequently, HLA-A*3301 and HLA-B*1402 were found to belong to the same ancestral haplotype . Further studies done in other Maghrebians did not reproduce these findings. Independently, Mokni-Baizig N et al.  noted that the HLA-A23 allele was absent in NPC patients compared to an 18% frequency in non tumor-bearing Tunisians. This is interesting because HLA-A23 is part of an African extended haplotype that includes a -B14 allele (HLA-A*2301/HLA-B*1403) . Moreover, Dardari R et al.  reported a significantly lower frequency of HLA-A9 serologic family alleles in NPC patients compared with controls. Since the HLA-A9 serologic family includes HLA-A23, this observation suggests that the reduced frequency of HLA-A9 reflects a lower frequency of the HLA-A*2301/HLA-B*1403 haplotype. Thus, previous studies in Maghrebians suggest that HLA-B14 and HLA-A alleles linked to the HLA-B14 serological family may be under-represented in patients with NPC compared with normal controls. Conversely, HLA-B5, -B13 and -B18 could be associated with increased risk of NPC.
Molecular typing provides enhanced discrimination for population studies and large data bases have accumulated in African populations including Maghrebians [25–29] These studies, however, did not include Tunisians (TU). Therefore, the aim of this study was to 1) establish a high resolution SBT data base in normal Tunisians matched for gender, age and geographical residence to a similar cohort of NPC patients; 2) compare HLA allele and haplotype frequencies between normal healthy Tunisians and the previously characterized Northern Moroccan Berber population (ME)  to provide a broader frame for the interpretation of the results from the Tunisian study; 3) test the validity of the reported decreased prevalence of HLA-B14 and increased prevalence of HLA-B5, -B13 and -B18 in patients with NPC based on SBT-derived, high resolution information; 4) Find novel alleles and haplotypes associated with NPC in Tunisian that may serve as putative genetic markers for further susceptibility studies.
Materials and methods
147 Patients with undifferentiated nasopharyngeal carcinoma were recruited from the Department of Radiation Oncology of Sousse Hospital, between 1991 and 2004. The patients with nasopharyngeal cancer had a mean age of 41.5 ± 16 years. The clinical stages ranged from II to IV (TNM classification, 1987). The diagnosis of cancer was confirmed by histopathology analyses. The histology was undifferentiated carcinoma (type III, WHO classification) in all cases [30, 31].
Patients were age and sex matched with 150 control subjects and selected from the same population living in the middle coast of Tunisia sharing the same environmental conditions. Control subjects having a mean age of 40 ± 12.5 years were healthy blood donors having no evidence of any personal or family history of cancer (or other serious illness). In order to respect the Hardy-Weinberg conditions highly stringent criteria were used for recruitment. To avoid any consanguinity bias, both patients and control groups were formed with unrelated individuals.
The control Tunisian population data were also compared to an available database of SBT typing performed in a population from North Morocco  to provide a frame of reference for this study in the context of other previously studied Maghrebian populations.
Specimen collection, preparation and HLA genotyping
Blood samples were obtained from all subjects consenting to the collection of blood,  through venipuncture and processed for DNA preparation. Extracted DNA was shipped to the Department of Transfusion Medicine, Clinical Center; National Institutes of Health where the genotyping was performed. SBT of HLA class I loci was performed as previously described . The primary PCR amplification reaction produced a 1.5 kb amplicon encompassing exon 1 through intron 4 of the HLA class I locus. All reagents necessary for primary amplification and sequencing were included in the HLA-A, HLA-B and HLA-C allele SEQR Sequenced Based Typing Kits (Abbott Diagnositics, Abbott Park, IL U.S.A.). After the primary amplification PCR products were purified from excess primers, dNTPs and genomic DNA using ExoSAP-IT (American Life Science, Cleveland, OH, U.S.A.). Each template was sequenced in the forward and reverse orientation for exon 2, exon 3 and exon 4 according to protocols supplied with the SBT kits. Excess dye terminators were removed from the sequencing products utilizing an ethanol precipitation method with absolute ethanol. The reaction products were reconstituted with 15 μl of Hi-Di™Formamide (PE Applied Biosystems/Perkin-Elmer, Foster City, CA, U.S.A.) and analyzed on the ABI Prism* 3700 DNA Analyzer with Dye Set file: Z and mobility file: DT3700POP6 [ET].
SBT information was used to compare our data base with recent reports on the prevalence of HLA alleles and haplotypes in African populations [25–29] The Hardy-Weinberg equilibrium was tested using the Guo and Thomson method , and gene diversity was estimated by maximum likelihood at each locus. Haplotype frequencies (HF) were estimated using an expectation-maximisation (EM) algorithm for multi-locus genotypic data when the gametic phase is not known employing the Arlequin population genetic software v3.01 [34, 35] The Ewens-Watterson homozygosity test [36, 37] was used to examine the presence of selective forces influencing allelic diversity at each locus (balancing selection or directional selection). The significance was calculated by using the Ewens-Watterson test and Slatkin exact test [38, 39]. HLA class I allele frequencies between Tunisians and Moroccans were compared using a χ2 analysis with Bonferroni correction. HLA class I (haplotipic) genetic diversity between the 2 populations was assessed by computing Fst genetic distance and the exact test of population differentiation using Arlequin. The most frequent 2 loci (Cw-B, A-Cw, A-B) and class I (A,-Cw,-B) haplotypes were described. Linkage disequilibrium (D) and relative D (D') between two alleles at two different loci and their level of significance (p) for 2 × 2 comparisons were calculated .
Comparisons between NPC and controls were performed first at the low resolution (two digit discrimination level) to compare our results with previously reported information based on low resolution serological typing. For this descriptive analysis data were simply compared using a 2 × 2 comparison without correction for number of tests. Further, SBT based comparisons between NPC and control populations were performed using allele frequencies and haplotypes calculated with the statistical tests discussed in the previous paragraph.
HLA frequencies among non-tumor bearing healthy Tunisians (TU)
To our knowledge no reports are available on SBT-defined HLA class I frequencies among healthy Tunisians. However, extensive information exists on other African/Maghrebian populations [25–29] We, therefore, first analyzed HLA class I frequencies in Tunisian in the context of available information in other Maghrebian populations. In particular, we performed a direct comparison with Piancatelli D. et al  SBT data on Moroccans.
Of the 150 samples obtained from normal control subjects only 147 HLA-A, 148 HLA-B and 147 HLA-Cw locus typings could be performed due to insufficient quantity or quality of the recovered genomic DNA in the remaining cases. Because of the recognized bimodal age distribution of NPC in Tunisia an effort was made to collect DNA representative of all age groups in the control population (Figure 1). The allele and haplotype frequencies for the HLA class I loci in the Tunisian population compared to Moroccan are reported in Tables 1, 2, 3. The Tunisian population was in Hardy-Weinberg equilibrium for all loci (HLA-A, p = 0.568 ± 0.000; HLA-B, p = 0.505 ± 0.000; HLA-Cw, p = 0.297 ± 0.000); the allelic distribution showed a very high level of gene diversity: HLA-A 0.94; HLA-B 0.96 and HLA-Cw 0.90. In the Moroccan population, these values were slightly lower for HLA-A (0.92) and HLA-B (0.95), and not different for HLA-Cw. Population comparisons showed that the two populations were not significantly different (Fst p value: 0.189) and no statistically significant differences were observed in allele frequencies between the two populations (Table 1).
By comparing in Tunisians the observed homozygosity F values with those expected under the Ewens-Watterson neutral model, a tendency toward balancing selection was observed, (observed F less than that expected under neutrality) although no significant difference in homozygosity was found (A locus: observed F: 0.068; expected F: 0.090; B locus: observed F: 0.048; expected F: 0.056; Cw locus: observed F: 0.101; expected F: 0.124).
Single locus analysis
A total of 35 HLA-A, 51 HLA-B and 24 HLA-Cw alleles were observed in Tunisians. The most frequently detected alleles in Tunisians and Moroccans were: HLA-A*020101(TU:17.3%, ME:17.8%), A*010101 (TU:8.8%, ME:13.7%), A*2301 (TU:7.8%, ME:8.2%), HLA-B*5001(TU:13.9%, ME:7.3%), B*440301(TU:7.1%, ME:11.6%), B*0801(TU:6.4%, ME:8.7%), HLA-Cw*0602(TU:21.8%, ME:18.3%), Cw*040101(TU:12.6%, ME:15.9%) and Cw*07 group (TU:19.3%, ME:18.4%).
Noteworthy was the presence of rare alleles in Tunisians, as classified according to the National Marrow Donor Program , such as B*5303 (originally identified in Hispanics but present also in Africans and Caucasians) and A*0285 (Caucasian) and uncommon alleles present in Africans, such as A*3010 (identified in Moroccans), A*0103, A*3004, A*3601, A*8001, B*1403, B*4202, B*5702, Cw*0210 (present in Moroccans), A*030103, B*4405, B*4427, B*5108, B*7301, Cw*1604 (present in Caucasians), A*0212, B*390602, B*5303 (present in Hispanics) and B*3920 (present in Blacks). Although not significant, a lower frequency of B*44, with a wider number of alleles (some of which uncommon), was found in Tunisians compared to Moroccans (12.8% vs. 21.7%, p = ns), with a frequency similar to South Africans (Zulu 12%, Cameroon 8.7%). The HLA-B*50 allele frequencies were high in Tunisians (15.3%) as in Moroccans (11.6%) and in accordance with previously reported information about Tunisians of mixed origin (10%) . Many allele groups showed a high variability (4–6 alleles, A*02, A*30, B*15, B*35, B*44). For instance, the HLA-B*15 group included the -B*1503 and -B*1510 alleles, that are present at high frequency in Sub-Saharan Africa. Overall, this data support the concept that Tunisians like other Maghrebian populations maintain a Northern African genetic identity which, however, includes traits derived from Sub-Saharan African  as well as Spanish/European influence . The higher level of differentiation present in Moroccans, expressed by a reduced number of alleles compared to TU, can be a genetic demonstration of a higher isolation of this ethnic group.
Extended haplotype analysis
One hundred and eighty different haplotypes were calculated in Tunisians. Data on haplotypes are shown in table 2 (2 loci) and table 3 (3 loci). The most common class I haplotypes were A*0201-B*5001-Cw*0602, shared by Tunisians (haplotype frequency; hf = 0.054) and Moroccans (hf = 0.040) and also present in Hispanics  (Table 3), and A*0201-B*4402-Cw*05, present in Hispanics, Africans and Caucasians [43, 42]. The A*2301-B*5001-Cw*0602 was not found in other populations, like the A*2301-B*5001 2-loci association. Many Cw-B haplotypes, in strong linkage disequilibrium, are common haplotypes found in Hispanic, Caucasians, African and African American (Black) populations, and have been detected in other Tunisians . Most haplotypes were shared between Tunisians and Moroccans (Table 2 and 3, in bold). Exceptions were represented by A*2301-B*5001and B4403-Cw0202, not found in other populations, and there were also not present (A*2301-B*5001) or present but in no significant linkage disequilibrium (B4403-Cw0202) in Moroccans.
HLA class I frequencies in Tunisian patients with NPC
Of the 147 samples obtained from Tunisian subjects with NPC only 135 HLA-A, 135 HLA-B and 136 HLA-Cw locus typings could be carried out due to insufficient quantity or quality of the recovered genomic DNA (Figure 1). The age distribution of the cohort of NPC patients in this study well reflected the recognized bimodal age distribution of NPC in Tunisia. Indeed, the study included 32 patients with NPC that were 25 year old or younger corresponding to a frequency of 23% identical to the frequency reported for this age group by Daoud J et al. .
To compare our results with previously reported studies on NPC in Maghrebian populations, we first analyzed allele frequencies at a low resolution (two digits) level that more closely approximate the serological nomenclature. Simple two by two comparison tables identified HLA-B18, -B51 and -B57 to be associated with a predisposition to develop NPC (Table 4). On the contrary, HLA-B14 and Cw-08 were associated with a decreased risk of developing NPC. The haplotype frequency was strongly reduced in patients with NPC. At the phenotype level, 2 of 136 patients with NPC carried the HLA Cw*0802-B*1402 haplotype (1.47%) while 11 of 148 (7.4%) Tunisian controls carried this haplotype (Fisher's test p2-value < 0.001). With the exception of HLA-Cw-08, these findings specifically reproduced associations that had been previously described in Maghrebian populations [22, 24] Thus, these results encouraged an in depth analysis of the HLA frequencies in patients with NPC from Tunisia at the high resolution level.
Extended haplotype analysis in Tunisian patients with NPC
The comparisons with the normal Tunisian and Moroccan populations were performed only for haplotypes found in Tunisians and/or patients with NPC (Table 5). Several associations were noted between NPC prevalence and HLA haplotypes. HLA-B*1801 was positively associated with the prevalence of NPC at the allele and haplotype level (specific NPC haplotype: A*020101-Cw*0701-B*180101, absent in Tunisians and Moroccans). The haplotype Cw*0501-B180101 was observed in Moroccans and, therefore, it may be less specific; HLA-B*570101 was also positively associated at the allele and haplotype level (specific NPC haplotype: A*010101-Cw*0602-B*570101, absent in Tunisians and Moroccans). Finally, HLA B*510101 was associated with an increased risk of developing NPC only at the allele but not the haplotype level.
SBT identified several sub-haplotypes including an HLA-B*14 allele (Cw*0802-B*1402, A*6802-B*1402 and A*1101-Cw*0802-B*1402) that appeared to be responsible for the negative relationship between HLA-B14 and the prevalence of NPC in Tunisians. However, the most striking negative association was observed between the HLA-B*1402/Cw*0802 haplotype and the prevalence of NPC (HF = 0.007 in patient with NPC compared to a HF of 0.037 in Tunisians and 0.040 in Moroccans). In particular, all the HLA-Cw*08 alleles were identified by SBT as HLA-Cw*0802 while HLA-B14, split this haplotype into 3 sub-haplotypes in which HLA-B*1402 occurred 11 times (HF = 0.037) and HLA-B*1401 and B*1403 three times each (none of them seen in patients with NPC).
There was a trend toward an age distribution bias with the HLA class I loci associated with NPC predisposition occurring more frequently in the younger age bracket (<25 year-old patients). In particular, HLA-Cw*0602 was present at a frequency of 44% (14 of 32) in 25 year old or younger NPC patients compared to 16% (6 of 37) in patients who had NPC diagnosed after the age of 50 (Fisher test p2-value = 0.02). Similarly the HLA-B*5701/Cw*0602 haplotype was present in 12% of cases below age 25, 6% of cases between 25 and 50 and only 3% of cases above age 50. This trend was, however, not significant possibly because of to the relatively small number of cases. Interestingly, HLA-Cw*0602 and its extended haplotype HLA-B*5701/Cw*0602 were particularly frequent among young patients. Of 11 patients age 16 or younger, 7 (64%) carried the HLA-Cw*0602 allele and 4 the extended haplotype (37%) compared to a frequency of 25% and 4% respectively in the older patients (Fisher's exact test p2-value = 0.01 and 0.003 respectively). Interestingly, the same haplotype was totally absent in the same age group and it occurred only in one case among the 148 control subjects studied. Thus, it is possible that this haplotype may include a predisposition factor that may bear a particular influence on the penetrance of NPC before other environmental factors could contribute to the onset of NPC in later age groups.
NPC is among common cancers the one whose prevalence is most likely affected by a genetic predisposition that can be uncovered through population studies. This seems to particularly apply to the Southern Asian population repeatedly shown to be affected at high incidence. Although viral and environmental factors play an important role in the ethio-pathogenesis of NPC, genetic factors may also predispose as documented by several family and/or case-control studies [45, 46]. The short arm of the human chromosome 6, is not the only genomic area that has been implicated in NPC predisposition as a large family study located a predisposition area in chromosome 3 . However, many case-control studies pointed at HLA class I as a genomic region strongly associated with the prevalence of NPC in Southern Chinese [5–12, 47–49] Such associations may bear two possible explanations; first the antigen presenting function of HLA molecules may be responsible for an altered effectiveness of the immune response against tumor associated antigens expressed by NPC in individuals bearing different HLA alleles . This logical explanation is, however, challenged by the lack of candidate epitopic determinants associated to particular HLA class I alleles that clearly affect the natural history of NPC. In addition, most HLA class I associations reported in Asian populations indirectly support this concept since they are primarily associated with increased risk of developing NPC. This suggests that lack of antigen presenting capacity may result in reduced protection against the growth of NPC. For instance, the HLA-A*0207 allele association reported by Hildesheim et al.  in Taiwanese Chinese was explained by suggesting a reduced antigen presentation potential of this allele compared with other HLA-A*02 alleles present in the same population such as HLA-A*0201 which are not associated with higher prevalence of NPC. Similarly, the higher prevalence of HLA-A*0207 in Asians compared to Caucasians  may also explain the higher prevalence of the disease in Asian populations. However, associations of HLA phenotypes with increase prevalence of NPC are difficult to attribute to a single allele considering the high frequency of heterozygosity and the co-dominant expression of the three HLA class I loci aimed at broadening the antigen presenting repertoire of each individual.
An alternative explanation for the association between HLA and the prevalence of NPC attributes a broader significance of HLA class I genes as markers for other predisposition factors genetically controlled within a locus in close linkage disequilibrium . For instance, even in Hildesheim's report, the extended haplotype HLA-A*0207/HLA-B*4601 was characterized by a higher odds ratio to develop NPC than the HLA-A*0207 allele alone suggesting a more complex genetic explanation for such association than the displacement by the HLA-A*0207 allele of other HLA class I alleles belonging to the HLA-A*02 family. In addition, our and others' experiences do not suggest strong differences in epitope binding and antigen presentation between HLA-A*0201 (the prototypic Caucasian HLA-A*02 allele) and HLA-A*0207. These alleles differ only by one amino acid substitution from a Tyrosine to a Cysteine in position 99 within the binding groove of the HLA molecule. This substitution does not seem to significantly affect binding and presentation of canonical HLA-A*0201-associated epitopes [50–52] as well as EBV-associated epitopes . Finally, lack of antigen presentation capacity by a specific HLA class I allele should be compensated by other alleles expressed by the same individual; therefore, favouring homozygosity of the unfavourable allele in affected individuals. However, most reports do not suggest higher homozygosity values in high prevalence populations suggesting lack of balancing selection. Thus, dominant genetic traits directly related to oncogenesis may play a more important role than HLA class I alleles which are co-dominantly expressed.
Maghrebs have been less exhaustively studied. However, serological analyses observed positive associations between NPC and HLA class I alleles. These include increased frequency of HLA-B5 in Algerians (38.2% vs 24.4%) , HLA-B13 in Tunisians (15.5% vs 4)  and HLA-B18 in Moroccans (relative risk = 4.14) . In addition, negative associations have been reported such as HLA-Aw33, -B14 in Algerians with NPC. However, after correction for the number of specificities tested, these differences were not statistically significant and were not pursued further. HLA-A*3301 and HLA-B*1402 have subsequently been recognized to belong to the same ancestral haplotype  suggesting that this haplotype may harbor more relevance to NPC than previously recognized. This impression was supported also by the identification of a negative association between HLA-A23 (known to belong to the extended HLA-A*2301/HLA-B*1403 haplotype ) and the prevalence of NPC . This observation converges toward the relevance of HLA-B14 related haplotypes which may bear protective effects against NPC in North Africans. Moreover, Dardari R et al.  reported a significantly lower frequency of the HLA-A9 serologic family alleles in NPC patients compared with controls. Since the latter study did not provide molecular resolution, we can only speculate that the reduced frequency of HLA-A9 reflected a lower frequency of the HLA-A*2301/HLA-B*1403 haplotype since HLA-A23 is part of the HLA-A9 serological family. Thus, previous studies in Maghrebians suggested that HLA-B14 and HLA-A alleles linked to the HLA-B14 serological family may be under-represented in patients with NPC compared with normal controls. Conversely, HLA-B5, -B13 and -B18 could be associated with increased risk of NPC.
Low resolution (two digits level) analysis of allele frequency in our study quite remarkably and specifically confirmed previous positive associations described between NPC and HLA class I alleles in Berber populations. Both HLA-B*18 and -B*51 (the latter representing a split of the serologically-defined HLA-B5 family) were confirmed by our study (Table 4). In addition, a novel association was observed between HLA-B-57 (belonging to the HLA-B17 serological family) and the prevalence of NPC. Only the positive association observed in Tunisians between HLA-B13 and NPC prevalence was not reproduced by this study . Negative associations were also identified by this study of which only HLA-B*14 (OR = 0.09, Fisher's test p-2 value = 0.0001) strongly reproduced previous findings while no association was noted with HLA-A*23 nor HLA-A*33. Interestingly, however, another strong negative association was noted between HLA-Cw*0802 and NPC (OR = 0.18, Fisher's test p-2 value = 0.003). This is of particular significance because both HLA-B*1402 and HLA-Cw*0802 belong to the same ancestral haplotype which was present not only in Tunisians but also Moroccan Berbers but practically absent in Tunisian patients with NPC (OR = 0.12, Fisher's test p-2 value = 0.001). This strong negative association suggests that this haplotype may bear a locus in strong linkage disequilibrium with these HLA alleles with dominant protective effects against the onset of NPC in North Africans such as a tumor suppressor gene with stronger stability in this phenotype. Since protective alleles have also been described in Asian populations (such as HLA-A11) , it could be postulated that some HLA alleles may bear immunological functions facilitating the immune surveillance against NPC. However, there is no known functional similarity between alleles described in Chinese and North Africans to support this hypothesis suggesting that other factors may bear on this predisposition. Indeed, in Asian populations, micro-satellite analyses within the short arm of Chromosome 6 identified associations between putative loci and NPC prevalence that are stronger than those with HLA class I alleles . Hu C-C et al.  using high-resolution micro satellite mapping recently described a carcinoma-susceptibility locus located in a segment of Chromosome 6p21.3 containing the class I loci. Similarly, Nh SP et al.  concluded that an HLA-related recessive mutation is responsible for the prevalence of NPC in Chinese. Finally, Tian W et al.  identified in the Southern Chinese Han population a positive association between a predisposition to develop NPC and exon 5 of the MICA-STR encoded within the HLA class I region. Moreover, a recent meta-analysis of chromosomal aberrations associated with NPC  identified a small region in Chromosome 6p21.2-p23 that is frequently amplified in NPC and contains various candidate genes that may be responsible for NPC . Thus, our findings suggest that there are haplotypes within or close to HLA class I region, particularly the HLA-B and -C loci that may bear genetic characteristics favoring or protecting against the onset of NPC in Tunisian populations at risk. This information may help identify individuals belonging to different ethnic groups whose genetic risk to develop NPC could be marked by HLA class I haplotypes. These prototypic phenotypes could be scanned with novel high throughput technologies covering genomic areas in high linkage disequilibrium with such haplotypes. Since the analysis of the Tunisian control population did not identify differences between Tunisians and Moroccans, considering that a common Berber substratum is shared by North African countries [59–63], it is possible that our findings could be relevant to the other Northern African populations justifying the extension of future studies to Moroccans and Algerians .
Some of the alleles discussed by this study have been associated with other pathologies. HLA-B*1403 is the predominant allele in patients affected by ankylosing spondylitis from West Africa, where HLA-B*27 is uncommon . HLA-B*51 and HLA-B*15 alleles have been associated with Behcet disease , B*44, B*51, B*57, B*15 with HPV infection and protection/susceptibility to in cervical carcinoma . The present study, may offer novel insight on the molecular genetics of this region that may influence disease, although obvious relationships between the information presented and it causality is still missing.
Khlat M: Cancer in Mediterranean migrants – based on studies in France and Australia. Cancer Causes Control. 1995, 6: 525-531. 10.1007/BF00054161.
Bouchardy C, Parkin DM, Wanner P, Khlat M: Cancer mortality among north African migrants in France. Int J Epidemiol. 1996, 25: 5-13. 10.1093/ije/25.1.5.
Feng BJ, Huang W, Shugart YY, Lee MK, Zhang F, Xia JC, Wang HY, Huang TB, Jian SW, Huang P: Genome-wide scan for familial nasopharyngeal carcinoma reveals evidence of linkage to chromosome 4. Nat Genet. 2002, 31: 395-399.
Chan AT, Teo PM, Huang DP: Pathogenesis and treatment of nasopharyngeal carcinoma. Semin Oncol. 2004, 31: 794-801. 10.1053/j.seminoncol.2004.09.008.
Jing J, Louie E, Henderson BE, Terasaki IP: Histocompatibility leukocyte antigen patterns in nasopharyngeal carcinoma cases from California. Natl Cancer Inst Monogr. 1977, 47: 156-
Cui SH, Lin Y: Apparent correlation between nasopharyngeal carcinoma and HLA phenotype. Zhong Hua Zhong Liu Zhi. 1982, 4: 249-253.
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. 10.1002/ijc.2910320206.
Ou B-X, Ruan H-Y, Fan Y: Study on association between HLA-A, B, C DR, and nasopharyngeal carcinoma in Guangzhou area. Ai Zheng. 1985, 4: 5-8.
Zhang JZ: Correlation between nasopharyngeal carcinoma (NPC) and HLA in Hunan Province. Zhong Hua Zhong Liu Zhi. 1986, 8: 170-172.
Hildesheim A, Apple RJ, Chen C-J, Wang SS, Cheng Y-J, Klitz W, Mack SJ, Chen IH, Hsu MM, Yang CS, Brinton LA, Levine PH, Erlich HA: Association of HLA class I and II alleles and extended haplotypes with nasopharyngeal carcinoma in Taiwan. J Natl Cancer Inst. 2002, 94: 1780-1789.
Pimtanothai N, Charoenwongse P, Mutirangura A, Hurley CK: Distribution of HLA-B alleles in nasopharyngeal carcinoma patients and normal controls in Thailand. Tissue Antigens. 2002, 59: 223-225. 10.1034/j.1399-0039.2002.590308.x.
Lu C-C, Chen J-C, Jin Y-T, Yang H-B, Chan S-H, Tsai S-T: Genetic susceptibility to nasopharyngeal carcinoma within the HLA-A locus in Taiwanese. Int J Cancer. 2003, 103: 745-751. 10.1002/ijc.10861.
Goldsmith DB, West TM, Morton R: HLA associations with nasopharyngeal carconoma in Southern Chinese: a meta-analysis. Clin Otolaryngol. 2002, 27: 61-67. 10.1046/j.0307-7772.2001.00529.x.
Krausa P, Browning MJ: A comprehensive PCR-SSP typing system for identification of HLA-A locus alleles. Tissue Antigens. 1996, 47: 237-244.
Marsh SG, Albert ED, Bodmer WF, Bontrop RE, Dupont B, Erlich HA, Geraghty DE, Hansen JA, Mach B, Mayr WR, Parham P, Petersdorf EW, Sasazuki T, Schreuder GM, Strominger JL, Svejgaard A, Terasaki PI: for factors of the HLA system, 2002. Eur J Immunogen. 2002, 29: 463-515. 10.1046/j.1365-2370.2002.00359.x.
Krausa P, Brywka M, Savage D, Hui KM, Bunce M, Ngai JL, Teo DL, Ong YW, Barouch D, Allsop C: Genetic polymorphism within HLA-A*02: significant allelic variation revealed in different populations. Tissue Antigens. 1995, 45: 223-231.
Browning M, Krausa P: Genetic diversity of HLA-A2: evolutionary and functional significance. Immunol Today. 1996, 17: 165-170. 10.1016/0167-5699(96)80614-1.
Simonis TB, Barracchini KC, Hackett JA, Maurer D, Siauw P, Tonai R: HLA-A*02 allele frequencies of different ethnic groups. Hum Immunol. 1997, 55: 99-
Chang YW, Hawkins BR: HLA class I and class II frequencies of a Hong Kong Chinese population based on bone marrow donor registry data. Human Immunol. 1997, 56: 125-135. 10.1016/S0198-8859(97)00108-0.
Clamp M, Andrews D, Barker D, Bevan P, Cameron G, Chen Y, Clark L, Cox T, Cuff J, Curwen V: Ensembl 2002: accommodating comparative genomics. Nucleic Acids Res. 2003, 31: 38-42. 10.1093/nar/gkg083.
Simons MJ: HLA and nasopharyngeal carcinoma: 30 years on. ASHI Quarterly. 2003, 27: 52-55.
Herait P, Tursz T, Guillard MY, Hanna K, Lipinski M, Micheau C, Sancho-Garnier H, Schwaab G, Cachin Y, Degos L: HLA-A, -B, and -DR antigens in North African patients with nasopharyngeal carcinoma. Tissue Antigens. 1983, 22: 335-341.
Mokni-Baizig N, Ayed K, Ayed FB, Ayed S, Sassi F, Ladgham A, Bel Hadj O, El May A: Association between HLA-A/-B antigens and -DRB1 alleles and nasopharyngeal carcinoma in Tunisia. Oncology. 2001, 61: 55-58. 10.1159/000055353.
Dardari R, Khyatti M, Jouhadi H, Benider A, Ettayebi H, Kahlain A, Mansouri A, El Gueddari B, Benslimane A: Study of human leukocyte antigen class I phenotypes in Moroccan patients with nasopharyngeal carcinoma. Int J Cancer. 2001, 92: 294-297. 10.1002/1097-0215(200102)9999:9999<::AID-IJC1177>3.0.CO;2-9.
Williams F, Meenagh A, Darke C, Acosta A, Daar AS, Gorodezky C, Hammond M, Nascimento E, Middleton D: Analysis of the distribution of HLA-B alleles in populations from five continents. Hum Immunol. 2001, 62: 645-650. 10.1016/S0198-8859(01)00247-6.
Ellis JM, Mack SJ, Leke RF, Quakyi I, Johnson AH, Hurley CK: Diversity is demonstrated in class I HLA-A and HLA-B alleles in Cameroon, Africa: description of HLA-A*0 *2612, *3006 and HLA-B*1403, *4016, *4703. Tissue Antigens. 3012, 56: 291-302. 10.1034/j.1399-0039.2000.560401.x.
Sanchez-Mazas A, Steiner QG, Grundschober C, Tiercy JM: The molecular determination of HLA-Cw alleles in the Mandenka (West Africa) reveals a close genetic relationship between Africans and Europeans. Tissue Antigens. 2000, 56: 303-312. 10.1034/j.1399-0039.2000.560402.x.
Piancatelli D, Canossi A, Aureli A, Oumhani K, Del Beato T, Di Rocco M, Liberatore G, Tessitore A, Witter K, El Aouad R, Adorno D: Human leukocyte antigen-A, -B, and -Cw polymorphism in a Berber population from North Morocco using sequence-based typing. Tissue Antigens. 2004, 63: 158-172. 10.1111/j.1399-0039.2004.00161.x.
Cao K, Moormann AM, Lyke KE, Masaberg C, Sumba OP, Doumbo OK, Koech D, Lancaster A, Nelson M, Meyer D, Single R, Hartzman RJ, Plowe CV, Kazura J, Mann DL, Sztein MB, Thomson G, Fernandez-Vina MA: Differentiation between African populations is evidenced by the diversity of alleles and haplotypes of HLA class I loci. Tissue Antigens. 2004, 63: 293-325. 10.1111/j.0001-2815.2004.00192.x.
Shanmugaratnam K, Tye CY, Goh EH, Chia KB: Etiological factors in nasopharyngeal carcinoma: a hospital-based, retrospective, case-control, questionnaire study. IARC Sci Publ. 1978, 199-212.
Shanmugaratnam K, Sobin LH: World Health Organization. International Histological Classification: Histological typing of tumors of the upper respiratory tract and ear. 1991, Berlin: Springer-Verlag, 2
Adams SD, Barracchini KC, Chen D, Robbins F, Wang L, Larsen P, Luhm R, Stroncek DF: Ambiguous allele combinations in HLA Class I and Class II sequence-based typing: when precise nucleotide sequencing leads to imprecise allele identification. J Transl Med. 2004, 2: 30-10.1186/1479-5876-2-30.
Guo SW, Thompson EA: Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics. 1992, 48: 361-372. 10.2307/2532296.
Arlequin Population Software. [http://cmpg.unibe.ch/software/arlequin3]
Excoffier L, Laval LG, Schneider S: Arlequin ver 3.0: An integrated software package for population genetics data analysis. Evol Bionformatics On Line. 2005, 1: 47-50.
Ewens WJ: The sampling theory of selectively neutral alleles. Theor Popul Biol. 1972, 3: 87-112. 10.1016/0040-5809(72)90035-4.
Watterson GA: The homozygosity test of neutrality. Genetics. 1978, 88: 405-415.
Slatkin M: An exact test for neutrality based on the Ewens sampling distribution. Genet Res. 1994, 64: 71-74.
Slatkin M: A correction to the exact test based on the Ewens sampling distribution. Genet Res. 1996, 68: 259-260.
Imanishi T, Azaka T, Kimura A, Tokunaga K, Gojobori T: Estimation of allele and haplotype frequencies for HLA and complement loci. HLA 1991. Edited by: Tsuji K, Aizawa M, Sesazuki T. 1992, Oxford: Oxford University Press
NPDP Bioinformatics. [http://bioinformatics.nmdp.org]
Allele frequencies in worldwide populations. 2007, [http://www.allelefrequencies.net]
Cao K, Hollenbach J, Shi X, Shi W, Chopek M, Fernandez-Vina MA: Analysis of the frequencies of HLA-A, B, and C alleles and haplotypes in the five major ethnic groups of the United States reveals high levels of diversity in these loci and contrasting distribution patterns in these populations. Hum Immunol. 2001, 62: 1009-1030. 10.1016/S0198-8859(01)00298-1.
Daoud J, Toumi N, Bouaziz M, Ghorbel A, Jlidi R, Drira MM, Frikha M: Nasopharyngeal carcinoma in childhood and adolescence: analysis of a series of 32 patients treated with combined chemotherapy and radiotherapy. Eur J Cancer. 2003, 39: 2349-2354. 10.1016/S0959-8049(03)00512-4.
Chan ATC, Teo PML, Johnson PJ: Nasopharyngeal carcinoma. Ann Oncol. 2002, 13: 1007-1015. 10.1093/annonc/mdf179.
McDermott AL, Dutt SN, Watkinson JC: The aetiology of nasopharyngeal carcinoma. Clin Otolaryngol. 2001, 26: 82-92. 10.1046/j.1365-2273.2001.00449.x.
Simons MJ, Day NE, Wee GB, Shanmugaratnam K, Ho HC, Wong SH, Ti TK, Yong NK, Darmalingam S, De-The G: Nasopharyngeal carcinoma V: immunogenetic studies of Southeast Asian ethnic groups with high and low risk for the tumor. Cancer Res. 1974, 34: 1192-1195.
Yuan BW: Correlation between nasopharyngeal carcinoma and HLA in Sichuan. Zhonghua Zhong Liu Za Zhi. 1988, 10: 263-266.
Pimtanothai N, Kangwanshiratada O, Charoenwongse P: Serological analysis of human leukocyte antigens-A and -B antigens in Thai patients with nasopharyngeal carcinoma. J Med Assoc Thai. 2003, 86: S237-S241.
del Guercio MF, Sidney J, Hermanson G, Perez C, Grey HM, Kubo RT, Sette A: Binding of a peptide antigen to multiple HLA alleles allows definition of an A2-like supertype. J Immunol. 1995, 154: 685-693.
Rivoltini L, Loftus DJ, Barracchini K, Arienti F, Mazzocchi A, Biddison WE, Salgaller ML, Appella E, Parmiani G, Marincola FM: Binding and presentation of peptides derived from melanoma antigens MART-1 and gp100 by HLA-A2 subtypes: implications for peptide-based immunotherapy. J Immunol. 1996, 156: 3882-3891.
Bettinotti M, Kim CJ, Lee K-H, Roden M, Cormier JN, Panelli M, Parker KK, Marincola FM: Stringent allele/epitope requirements for MART-1/Melan A immunodominance: implications for peptide-based immunotherapy. J Immunol. 1998, 161: 877-889.
Provenzano M, Selleri S, Jin P, Wang E, Werden R, Slezak S, Adams SD, Panelli MC, Leitman SF, Stroncek DF, Marincola FM: Comprehensive epitope mapping of the Epstein-Barr virus latent membrane protein-2 in normal, non tumor-bearing individuals. Cancer Immunol Immunother. 2006
Lu CC, Chen JC, Tsai ST, Jin YT, Tsai 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.
Hu SP, Day NE, Li DR, Luben RN, Cai KL, Ou-Yang T, Li B, Lu XZ, Ponder BA: Further evidence for an HLA-related recessive mutation in nasopharyngeal carcinoma among the Chinese. Br J Cancer. 2005, 92: 967-970. 10.1038/sj.bjc.6602347.
Tian W, Zeng XM, Li LX, Jin HK, Luo QZ, Wang F, Guo SS, Cao Y: Gender-specific associations between MICA-STR and nasopharyngeal carcinoma in a southern Chinese Han population. Immunogenetics. 2006, 58: 113-121. 10.1007/s00251-006-0093-6.
Li X, Wang E, Zhao YD, Ren JQ, Jin P, Yao KT, Marincola FM: Chromosomal imbalances in nasopharyngeal carcinoma: a meta-analysis of comparative genomic hybridization results. J Transl Med. 2006, 4: 4-10.1186/1479-5876-4-4.
Jalbout M, Bouaouina N, Gargouri J, Corbex M, Ben Ahmed S, Chouchane L: Polymorphism of the stress protein HSP70-2 gene is associated with the susceptibility to the nasopharyngeal carcinoma. Cancer Lett. 2003, 193: 75-81. 10.1016/S0304-3835(02)00697-3.
Bosch E, Calafell F, Perez-Lezaun A, Comas D, Mateu E, Bertranpetit J: Population history of north Africa: evidence from classical genetic markers. Hum Biol. 1997, 69: 295-311.
Chaabani H, Cox DW: Genetic characterization and origin of Tunisian Berbers. Hum Hered. 1988, 38: 308-316.
Gomez-Casado E, del Moral P, Martinez-Laso J, Garcia-Gomez A, Allende L, Silvera-Redondo C, Longas J, Gonzalez-Hevilla M, Kandil M, Zamora J, Arnaiz-Villena A: HLA genes in Arabic-speaking Moroccans: close relatedness to Berbers and Iberians. Tissue Antigens. 2000, 55: 239-249. 10.1034/j.1399-0039.2000.550307.x.
Lucotte G, Aouizerate A, Berriche S: Y-chromosome DNA haplotypes in North African populations. Hum Biol. 2000, 72: 473-480.
Myles S, Bouzekri N, Haverfield E, Cherkaoui M, Dugoujon JM, Ward R: Genetic evidence in support of a shared Eurasian-North African dairying origin. Hum Genet. 2005, 117: 34-42. 10.1007/s00439-005-1266-3.
Lopez-Larrea C, Mijiyawa M, Gonzalez S, Fernandez-Morera JL, Blanco-Gelaz MA, Martinez-Borra , Lopez-Vazquez A: Association of ankylosing spondylitis with HLA-B*1403 in a West African population. Arthritis Rheum. 2002, 46: 2968-2971. 10.1002/art.10584.
Choukri F, Chakib A, Himmich H, Marih L, Caillat-Zucman S: HLA-B phenotype modifies the course of Behcet's disease in Moroccan patients. Tissue Antigens. 2003, 61: 92-96. 10.1034/j.1399-0039.2003.610110.x.
Zehbe I, Mytilineos J, Wikstrom I, Henriksen R, Edler L, Tommasino M: Association between human papillomavirus 16 E6 variants and human leukocyte antigen class I polymorphism in cervical cancer of Swedish women. Hum Immunol. 2003, 64: 538-542. 10.1016/S0198-8859(03)00033-8.
The authors wish to thank Dr. Khadija Oumhani and Prof. Rajae El Aouad for their precious contribution with the study of the original Moroccan samples.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Li, X., Ghandri, N., Piancatelli, D. et al. Associations between HLA Class I alleles and the prevalence of nasopharyngeal carcinoma (NPC) among Tunisians. J Transl Med 5, 22 (2007). https://doi.org/10.1186/1479-5876-5-22