Abstract

We investigated polymorphisms of the human leukocyte antigen (HLA) class I (A, B, and C) loci of a Han population (n, 239) from the Yunnan province, Southwest China, using high-resolution polymerase chain reaction-Luminex (PCR-Luminex) typing. We combined the HLA data from this study with the KIR genotypes from a previous study of this Han population to analyze the combination of KIR/HLA ligands. A total of 27 HLA-A, 54 HLA-B, and 31 HLA-C alleles were found in this population. The frequencies of 11:01, 24:02, 40:01, 46:01, 01:02, 03:04, and 07:02 were all > 10%. The following haplotypes were common, with frequencies > 5%: 1 A-B (02:07-46:01), 2 A-C (02:07-01:02, and 11:01-07:02), 4 C-B (13:01-03:04, 40:01-07:02, 46:01-01:02 and 58:01-03:02), and 1 A-C-B (02:07-01:02-46:01). Analysis of KIR3D and their ligands HLA-A3/A11 and HLA-Bw4 showed that the frequencies of 3DL2+-A3/A11+ and 3DL2+-A3/A11 were 0.527 and 0.473, and the frequencies of 3DL1+-Bw4+, 3DL1+-Bw4, 3DL1-Bw4+, and 3DL1-Bw4 were 0.552, 0.397, 0.038, and 0.013, respectively. The results of KIR/HLA-C combination analysis showed that all individuals had at least one inhibitory or activating KIR/HLA-C pair, and one KIR/HLA-C pair was the most frequent (157/239), followed by two pairs (46/239), three pairs (33/239), and no pairs (3/239). Comparison of KIR gene and HLA gene and their pair frequency between Yunnan Han and the isolated Han (FYDH) who also lived in Yunnan province showed no significant difference () in KIR frequencies, but significant differences () for some HLA allele frequencies. In addition, there was no significant difference () between the two populations for KIR/HLA pairs.

1. Introduction

The major histocompatibility complex (MHC) genomic region or human leukocyte antigen (HLA) super-locus at 6p21 encodes at least 132 protein-coding genes, including six classical transplantation HLA genes with multiple roles in immune regulation [1]. The high degree of diversity of alleles and haplotypes makes the classical HLA loci particularly valuable for investigating the origins and migration of human populations and for disease association studies.

Killer cell immunoglobulin-like receptors (KIRs), which are found on natural killer (NK) cells, have a central role in immune responses [2]. The KIR genes are located on chromosome 19q13.4 within the leukocyte receptor complex, which is clustered in one of the most variable regions of the human genome and shows extensive polymorphism in terms of both gene content and sequence polymorphism [2, 3].

HLA class I molecules and the KIR serve as ligand and receptor. Certain combinations of KIRs and HLA molecules can transmit activating or inhibitory signals to regulate the function of NK cells [4]. For HLA-A molecules, several studies have reported that KIR3DL2 interacted with HLA-A3 and -A11 allotypes [57]. For HLA-B molecules, the KIR3DL1 interacts with HLA-B Bw4 allotypes [6, 8, 9], and no KIRs have been found to bind the HLA-B Bw6 allotypes with high affinity [6, 10]. For HLA-C molecules, on the basis of a dimorphism (asparagine or lysine) at position 80 in the alpha helix, HLA-C alleles can be placed into one of two categories, C1 and C2, which are recognized differentially by the KIR receptors. NK cells that express KIR2DL2 and KIR2DL3 interact with HLA-C1 allotypes, whereas KIR2DL1 and KIR2DS1 interact with HLA-C2 allotypes [11, 12].

Chinese populations consist of 56 ethnic groups and Han is the largest population (~93%). The Han were derived from the Huaxia tribe living in the area of the Yellow River, who then migrated to the area of the Yangtze River and have now spread throughout the country.

On the basis of HLA-A, -B, and -DRB1 allele frequencies, our earlier work showed that the Han population can be divided into Southern and Northern Han [13]. We previously reported the distribution of HLA-A, -B, and -DRB1 genes and haplotypes at low resolution in a Han population (, 129) from the Yunnan province (YNH) and found that it had characteristics of both northern and southern Chinese groups based on HLA data [13]. In current study, we investigated polymorphisms of HLA class I (A, B, and C) loci of a Han population from Yunnan (, 239) using high-resolution polymerase chain reaction-Luminex (PCR-Luminex) typing to get the high-resolution HLA data of YNH. In addition, we combined the HLA data with KIR genotypes of the YNH to analyze the combination of KIR/HLA pairs and looked for any difference between the YNH and an isolated Han subpopulation [1416] that lived in the Fenyandong (FYDH) region in Yunnan province.

2. Materials and Methods

2.1. Participants

A total of 239 unrelated individuals, which were from the Yunnan Province of Southwest of China (YNH), were included in this study and the geographic origin, nationality, and pedigree (unrelated through at least three generations) of each individual were ascertained before sampling. All individuals were healthy and provided informed consent.

2.2. HLA Typing

Blood samples were collected and genomic DNA was extracted from peripheral lymphocytes using the QIAamp Blood Kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. DNA-based HLA typing was done with the Luminex multianalyte profiling system (xMAP) using the WAKFlow HLA typing kit (Wakunaga, Hiroshima, Japan).

2.3. KIR Typing

The results of the KIR typing in the Yunnan Han population are published [16] and here we have used the KIR results to analyze their relationship with HLA-A3/A11, HLA-Bw4, and the HLA-C1 and HLA-C2 groups that we obtained in this study.

2.4. Statistical Analysis

Allele and inferred haplotype frequencies (deduced from the phenotype) were calculated by the expectation-maximization algorithm in the Pypop software on the basis of the genotyping results [17, 18]. The HLA loci observed and expected heterozygosity were calculated by using MolKin v3.0 software [19]. For each locus, the Hardy-Weinberg equilibrium was assessed using the Guo and Thompson method [20] and the Ewens-Watterson homozygosity test of neutrality was done [21, 22]. Linkage disequilibrium (LD) was analyzed using the Pypop software [17, 18]. The HLA-A3/A11 and HLA-Bw4 groups were used to analyze KIR/HLA combinations. The HLA-C1 and HLA-C2 groups were also used to analyze KIR/HLA combinations. The presence of two inhibitory (2DL2/2DL3+HLA-C1 and 2DL1+HLA-C2) and one activating (2DS1+HLA-C2) receptor-ligand pairs was evaluated for each individual. The differences in significance between the HLA/KIR pair frequencies were determined by a contingency test. Statistical significance was defined at the 5% level.

3. Results

3.1. Hardy-Weinberg’s Equilibrium Observed Heterozygosity and Expected Heterozygosity and Ewens-Watterson Test of Neutrality for HLA-A, -B, and -C Loci

Table 1 showed the HLA-A, -B, and -C loci in Hardy-Weinberg’s equilibrium () and the results for the Evens-Watterson test of neutrality for HLA-A, -B, and -C loci in YNH. Data on HLA-A, -B, and -C did not reject the neutral model (). The observed and expected heterozygosity of HLA-A, -B, and -C loci were shown in Table 1.

3.2. HLA-A, -B, and -C Allele Frequencies

In all, 27 HLA-A, 54 HLA-B, and 31 HLA-C high resolution alleles were observed (Table 2). The three most common alleles at (1) HLA-A locus were 0.268 for HLA-11:01, 0.190 for HLA-24:02, and 0.094 for HLA-02:07; (2) HLA-B locus were 0.144 for HLA-46:01, 0.117 for HLA-40:01, and 0.084 for HLA-13:01; and (3) HLA-C locus were 0.203 for HLA-01:02, 0.153 for HLA-07:02, and 0.117 for HLA-03:04.

Compared with HLA alleles of an isolated Han population (FYDH) in Yunnan province [14], the frequencies of some alleles were significantly different () between these two Han populations. On the HLA-A locus, the frequencies of 02:01 and 02:07 were 0.086 and 0.094, respectively, in YNH and 0.010 and 0.183, respectively, in FYDH. On the HLA-B locus, 40:01 was present in YNH at frequencies of >10% (0.117) but were observed only at frequencies of 0.045 in FYDH. By contrast, the common allele of HLA-15:02 in FYDH (0.124) was present at a low frequency in YNH (0.044). At the HLA-C locus, HLA-07:02 was more common in YNH compared to FYDH (0.153 versus 0.094), whereas HLA-08:01 was more common in FYDH compared to YNH (0.149 versus 0.040).

3.3. The Linkage Disequilibrium of HLA-A, -B, and -C Loci and the Haplotype Frequencies

The likelihood ratio test of linkage disequilibrium found that all pairwise associations (Table 3) were statistically significant (). Therefore, we estimated the HLA haplotype frequencies bearing 2-3 loci and show the 2-3 loci haplotype frequencies of >1% in Table 4. The most common A-B haplotype was 02:07-46:01 (0.073), the most common A-C haplotypes were 02:07-01:02 (0.068) and 11:01-07:02 (0.070), and the most common C-B haplotypes were 13:01-03:04 (0.071), 40:01-07:02 (0.052), 46:01-01:02 (0.132), and 58:01-03:02 (0.054). The most common A-C-B haplotype was 02:07-01:02-46:01 (0.067).

The difference between FYDH and YNH for the HLA was obvious also at the haplotype level [14]. 11:01-08:01-15:02 was a common haplotype in FYDH (0.124), but the frequency of this haplotype was only 0.014 in YNH. Similarly, the frequency of 02:07-01:02-46:01 was 0.126 in FYDH but 0.067 in YNH, and the frequency of 33:03-03:02-58:01 was 0.074 in FYDH but 0.042 in YNH.

3.4. Two KIR3D and Their Ligands HLA-A3 and -A11

The distribution of two KIR3D genes and their HLA-A3/A11 ligands is given in Table 5. In all, the frequencies of 3DL2+-A3/A11+ and 3DL2+-A3/A11+ were 0.527 and 0.473 in YNH. No significant differences () were detected in the statistical comparison between the YNH and the FYDH.

3.5. Two KIR3D and Their Ligands HLA-Bw4

The distribution of two KIR3D genes and their HLA-Bw4 ligands is given in Table 6. In all, the frequencies of 3DL1+-Bw4+, 3DL1+-Bw4, 3DL1-Bw4+, and 3DL1-Bw4 were 0.552, 0.397, 0.038, and 0.013, respectively in YNH. There was no statistical difference () in the frequencies between the YNH and FYDH populations.

3.6. Combinations of KIR/HLA-C

The HLA and KIR genes segregate on two different chromosomes but have a functional correlation owing to the KIR/HLA-C ligand-receptor pairs. Using our KIR data from earlier studies [15], we analyzed two inhibitory (2DL2/2DL3+HLA-C1 and 2DL1+HLA-C2) and one activating (2DS1+HLA-C2) receptor-ligand pairs (Table 7). First, we analyzed inhibitory and activating KIR/HLA-C pairs in YNH. We observed that all individuals had at least one inhibitory or activating KIR/HLA-C pair. One KIR/HLA-C pair was the most frequent (157/239), followed by two pairs (46/239), three pairs (33/239), and no pairs (3/239). Second, we differentiated between inhibitory and activating KIR/HLA-C pairs. The frequencies for different numbers of inhibitory KIR/HLA-C pairs in YHN were 75/239 for two inhibitory KIR/HLA-C pairs (KIR2DL2/2DL3+HLA-C1, KIR2DL1+HLA-C2), 161/239 for one inhibitory KIR/HLA-C pair (152 KIR2DL2/2DL3+HLA-C1 and 9 KIR2DL1+HLA-C2), and three for no pairs (3/239). For the activating KIR/HLA-C pairs, 37/239 were for one KIR/HLA-C pair (35 KIR2DS1+HLA-C2) and 202/239 were for no pairs. There was no statistical difference () in the frequencies of KIR/HLA-C ligand-receptor pairs between the YNH and FYDH populations.

4. Discussion

Earlier, we reported the distribution of HLA-A, -B, and -DRB1 genes and haplotypes at low resolution in a Yunnan Han (YNH) population and found that it had characteristics of both northern and southern Chinese groups, although they live in the southwest of China [13]. In the present study, we genotyped the HLA-A, HLA-B, and HLA-C genes at high resolution in YNH and compared their HLA alleles, haplotypes, and combinations of KIR/HLA to an isolated Han subpopulation [14–16] that lived in the Fenyandong (FYDH) region of the Yunnan province.

The FYDH subpopulation became geographically isolated and lived in the FYDH regionof Yunnan province (Southwest China) after the group had emigrated from other places more than 200 years ago [23]. Since then, the FYDH did not usually communicate with other groups, as it took at least one day to walk to the nearest town. Our previous studies showed that the FYDH has a particular distribution of HLA alleles and KIR genotypes, and there is a difference between the YNH and FYDH populations in the distribution of HLA-A and -B alleles at low resolution [14]. In the present study, a significant difference was observed between the YNH and FYDH for the high resolution, HLA allele frequencies. Firstly, the number of alleles was less in the FYDH than the YNH population. Secondly, the HLA alleles of 11:01 (0.317), 02:07 (0.183), 15:02 (0.124), 46:01 (0.198), 08:01 (0.149), and 07:02 (0.094) were at a substantially higher frequency in the FYDH subpopulation than in the YNH (Table 2) or in other Chinese Han populations. Thirdly, the frequencies of common haplotypes such as 02:07- 01:02-46:01 (0.126 vs 0.067) and 11:01-08:01-15:02 (0.124 vs 0.014) were higher in the FYDH than in the YNH, respectively. Therefore, the high frequencies of some of these alleles in the FYDH could be a reflection of genetic drift, founder effects, and/or population bottlenecks [24, 25], suggesting that genetic drift and/or selection and subsequent geographic isolation has played a more important role in distribution of HLA alleles in FYDH [14] than the YNH.

By contrast, comparison of KIR frequencies between YNH and FYDH showed that the two populations were more similar for the KIR genes than the HLA genes [15, 16]. Firstly, the 15 KIR genes frequencies in YNH did not differ significantly compared to FYDH except for KIR2DS5 (). Secondly, more than 40 KIR genotypes showed no significant difference between YNH and FYDH, except for genotype 2 and genotype 33, respectively, and the KIR AA genotype was the most frequent in YNH and FYDH (0.478 and 0.570, resp.). Thirdly, the similar ratio of A and B haplotype frequencies in YNH and FYDH was observed (2.39 : 1 vs 2.91 : 1, resp.) [15, 16]. The KIR genes, genotypes, and haplotypes in Han populations were the same as those in Asian populations but were different from those in African and Caucasian populations [2630]. Moreover, the significant difference of two activating genes (KIR2DS4 and 2DS3) and one inhibitory gene (KIR2DL2) that was observed between northern and southern Han did not show any significant difference between the YNH and FYDH in the present study.

When the HLA/KIR combination was compared, there was still no significant difference between these two Han populations. Individuals with KIR3DL2 genes coupled with HLA-A3/A11 were predominant in YNH and FYDH with similar frequencies. Individuals with KIR3DL1 genes coupled with HLA-Bw4 were predominant in YNH and FYDH with similar frequencies. Individuals that did not have the KIR3DL1 genes coupled with HLA-Bw4 were also common in the YNH and the FYDH. When the HLA-C1/C2 ligands were considered, the frequencies of the individuals with HLA-C1 homozygosity coupled to the KIR2DLl and KIR2DL2/2DL3 genes (but lacking the KIR2DSl gene) were predominant in the YNH and the FYDH with similar frequencies around 50% (data no shown). The distribution of KIR2D-HLA-C1/C2 showed no difference between YNH and FYDH (data not shown). These data also showed the existence of a preponderant inhibitory KIR/HLA signaling pathway in these populations.

The extensive diversity of the HLA and KIR genes and the central role of their interaction in immune responses to different pathogens are expected to favor the coevolution of these loci as suggested in several disease, population, and comparative genetic studies across primate species [11, 31, 32]. However, in the present study, the main difference was observed for HLA but not for KIR and HLA/KIR combinations. We assumed this difference might be affected for two reasons. Firstly, the HLA genes were genotyped for SNPs at a high resolution level, allowing allelic differences to be more easily detected. Secondly, the KIR genes were detected only at the level of their presence or absence in each individual and the KIR SNPs were not genotyped. The difference between the KIR gene alleles of the YNH and FYDH might exist more obviously at the SNP level. In addition, the KIR gene cluster is located in a rapidly evolving genomic region where the activating KIR genes evolved from inhibitory KIR genes and were short lived in comparison to genes encoding the inhibitory KIR genes [33]. Our previous study as well as other population studies have suggested that the geographic distribution of populations was more closely related to activating KIR genes than to inhibitory KIR genes [16, 34]. The only difference in the KIR gene cluster observed in the present study was for the KIR2DS5. According to our epidemic data, the FYDH has been threatened with infectious diseases such as hepatitis A virus (HAV) infection with the antibody to HAV detected in more than 90% of the population. Thus, we deduced that after FYDH isolation, the activating KIR genes have evolved more rapidly than the inhibitory KIR genes in response to epidemics.

In conclusion, the current study determined the distribution of KIR/HLA gene combinations in the YNH. In combination with earlier published data, we demonstrated that the YNH population has a particular distribution of HLA alleles, which is different compared to FYDH population. However, the KIR genotypes and the combinations of KIR/HLA were similar in both Han populations. Our results will help us to understand the role of KIR/HLA in autoimmune diseases and viral infections. Epidemiology of the YNH population and the distributions of other immune system genes need to be investigated further.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contribution

The contributions of Yunsong Shen and Danfeng Cao are equal.

Acknowledgments

This work was supported by grants from National Natural Science Foundation of China (81160197 and 31270030), Peking Union Medical College Talented Youth Project, and Yunnan Provincial Science and Technology Department (2013FZ131 and 2012FB064).