KIR diversity in three ethnic minority populations in China
Journal of Translational Medicine volume 13, Article number: 221 (2015)
Killer cell immunoglobulin-like receptors (KIRs) show extensive variation in genetic content and allelic polymorphi sms among different populations.
Materials and methods
We analyzed the distribution of KIR genes in the Tibetan ethnic minority of Lhasa city, the Uyghur and Kazakh ethnic minorities of Urumqi city populations in China. Genotyping of 16 KIR genes was tested in 479 randomly selected individuals using the multiple PCR-SSP method.
A total of 42 KIR genotypes were detected, of which, 29 were predicted to be AB genotypes, 12 were BB genotypes and one was AA genotypes. 27 KIR genotypes were identified in Kazakhs, 30 KIR genotypes were identified in Uyghurs and 20 KIR genotypes were identified in Tibetans. The predominant genotype 1(AA genotypes) occurred most frequently in Tibetans (52.7%, 118/224), Kazakhs (43.2%, 54/125) and Uyghurs (34.9%, 45/130). Not only the four framework genes were present in all individuals, but the pseudogene 2DP1 could also be detected in all Uyghur individuals. Tibetans were different from Kazakh and Uyghur groups in KIR genetic content and KIR allelic variation. Intriguingly, Tibetans (29.5%, 66/224) had lower frequencies of 2DS4-v when compared with Uyghurs (60.8%, 79/130) and Kazakh s (59.2%, 74/125). Uyghurs (25.4%, 33/130) displayed higher frequencies of Bx genotypes with C4Tx (absence of KIR3DS1-2DL5-2DS5-2DS1) than both Kazakhs (11.2%, 14/125) and Tibetans (3.6%, 8/224).
The study showed that profile of KIR genotypes in three ethnic minority populations in China displayed ethnic diversity. It could be valuable for enriching the ethnical information resources for KIR gene, as well as facilitating further research on KIR-related diseases.
Killer cell immunoglobulin-like receptors (KIRs) are a group of glycoproteins expressed on the surface of Natural killer (NK) cells [1, 2]. NK cells are important components of the innate immune system that take part in various immune responses to different infectious agents .
KIR genes are located on human chromosome 19q13.4. Up to now, 16 KIR gene loci have been identified, including two pseudogenes (KIR 2DP1 and 3DP1). Four framework KIR genes (2DL4, 3DL2, 3DL3 and 3DP1) exist in almost all populations . Fourteen functional KIR genes have been identified and confirmed, of which seven are inhibitory (KIR2DL1-3, 2DL5 and 3DL1-3) and six are activating (KIR2DS1-5 and 3DS1), with only KIR2DL4 showing both inhibitory and activating potential [1, 5, 21]. According to the numbers and types, KIR haplotypes are broadly classified into two groups: group A haplotypes have a fixed KIR gene (2DL1, 2DL3-4, 3DL1-3, 2DS4, 2DP1 and 3DP1), in contrast, group B haplotypes have variations in genetic content which is comprised of several genes and alleles that are not found in A haplotypes. Especially, KIR2DS1-3, 2DS5, 2DL2, 2DL5 and 3DS1 are associated only with group B haplotypes, thus B haplotypes generally encode more activating KIR receptors than A haplotypes which encode a single activating receptor (KIR2DS4). Both group A and B are present in all populations, but their frequencies vary considerably [1, 6]. Many individuals carry two copies of A haplotypes (AA genotypes). The remainders carry Bx genotypes, which are consisted of either one copy of A haplotype and one copy of B haplotype (AB genotypes) or two copies of B haplotypes (BB genotypes). Based on the presence and absence of two distinct gene clusters (C4, KIR2DS2-2DL2-2DS3-2DL5; T4, KIR3DS1-2DL5-2DS5-2DS1), Bx genotype carries are divided into four subsets: C4Tx (presence of C4 and absence of T4), CxT4 (absence of C4 and presence of T4), C4T4 (presence of both C4 and T4), and CxTx (absence of both C4 and T4) [1, 7, 8].
Usually, Group B haplotypes encode more activating KIR genes. In individuals with AA genotype, the alleles of KIR2DS4 can encode either functional (KIR2DS4-f) or non-functional (KIR2DS4-v) variants [9, 10]. KIR2DS4-v alleles differ from the KIR2DS4-f alleles in that the former have a 22-bp deletion in exon 5, which leads to a frame shift mutation and produces a premature stop codon, preventing the formation of a functional membrane-bound receptor domain . Therefore, individuals with AA genotypes which harbor KIR2DS4-v alleles will not have a membrane-bound KIR2DS4 protein, but rather a soluble form of KIR2DS4. KIR2DS4-f and 2DS4-v frequencies vary amongst different populations. Increasing attention has been paid to the role of KIR genetic content and allelic variation on infectious diseases such as Hepatitis C and HIV, with some studies dedicating a specific focus on KIR 2DS4 [9–14]. KIR gene diversities have been studied in many different geographical populations as previously reported [1, 9, 15–23]. There are 56 ethnic groups in China . Most of the studies on KIR gene diversity have been reported on Han populations in China [18–21]. A large-scale survey on transfusion-transmitted HIV-1/2 infection among Chinese blood donors conducted by our institute showed that the positive rate for HIV infection was higher in some ethnic minorities (including data accumulated at Urumqi Blood Center, Xinjiang) than Han majority donors . Expanding our understanding of the ethnic intermarriage and possibly random demographic factors could help us to determine the variation of KIR gene frequencies, which might be useful for future research on ethnicity-based diseases. Hence, in this study we chose the main minorities in Urumqi city of Northwest China, the Uyghur (comprised 10,069,346 persons , presented a typical mixture of Eastern and Western anthropometric traits ) and Kazakh (comprised 1,462,588 persons ), and the Tibetan (the major ethnic minority in Tibet, comprised 6,282,000 persons ), mostly living in Lhasa city of Southwest China, as the research objects. (Additional file 1: Figure S1).
Samples and DNA isolation
A total of 479 randomly selected donors were studied from three ethnic minority populations in China, including Tibetan (N = 224) volunteer blood donors in Lhasa (April 2011 to March 2012), Kazakh (N = 125) and Uyghur (N = 130) volunteer blood donors in Urumqi (March to August 2012). Informed consents of all samples used in this study were approved by the ethical Committee of Institute of Blood Transfusion. Genomic DNA was extracted from PBMC by TIANGEN blood DNA kit (Cat: #DP318-03), following the manufacturer’s instructions. The concentration of extracted DNA was adjusted to 50–100 ng/μL, and the O.D. 260/O.D. 280 ratio ranged from 1.6 to 1.9.
KIR PCR-SSP genotyping
KIR genes were typed for the presence or absence of the 14 KIR genes, including KIR2DL1-5, 2DS1-5, 3DL1-3 and the 2 pseudogenes (3DP1 and 2DP1), using multiplex PCR-SSP. In order to ensure accuracy, two different sets of primers for amplification (except 3DP1 and 2DS1) were used while HLA-DRB1 intron amplification primers were added to the third well as a positive control. Primer designs were adopted from those previously reported by Campbell et al. . The primer arrangement was referred in Table 1 PCR master-mix for the multiplex PCR-SSP which contained 3 μL of primer (10 μM) pairs, 1.5 μL of ddH2O, 0.5 μL of DNA, 5 μL of GoTaq green master mix (Promega, Cat: #M7123), for a total of 10 μL. The conditions for PCR cycles were as followed: five cycles at 94°C for 3 min, 94°C for 15 s, 65°C for 15 s, 72°C for 30 s; twenty-one cycles at 94°C for 15 s, 60°C for 15 s, 72°C for 30 s; four cycles at 94°C for 15 s, 55°C for 1 min, 72°C for 2 min; lastly, at 72°C for 7 min and stored at 4°C. PCR products were visualized under UV light after electrophoresis in 3% agarose gel and stained with ethidium bromide, using DL1000 DNA Marker (TIANGEN, Cat: #MD108).
Genotypic frequency = n/N (n represented numbers of positive genotypes and N was the total number of individuals tested). The carrier frequencies (CF) of the KIR genes were determined by dividing the number of positive typing reactions by the total number of typed individuals. The estimated gene frequencies were calculated using the formula GF = 1 − (1 − CF)1/2. The frequencies of A and B haplotypes were calculated using the following formula: group A = 2nAA + nAB/2N and group-B = 2nBB + nAB/2N, where nAA, nAB, and nBB were the numbers of AA, AB and BB genotypes. Differences in the carrier frequencies between the studied population and other populations previously published were assessed by the standard Chi square test (x2) using statistical software PEMS3.1 Medicine, and p < 0.05 was considered to be statistically significant. Subsets of KIR genetic content profiles and their frequencies in populations were carried out using statistical software Graph Pad Prism5. Principal components analysis (PCA) of 22 populations  was performed using an R package: psych software. Genetic distances were calculated by Nei’s method using phylip . On the basis of the Nei’s genetic distance, a dendrogram was constructed by the neighbor-joining method  and visualized with MEGA software in order to compare the frequencies of KIR genes in 22 populations , using bootstrap values are calculated from 100 trees.
KIR genotypes and haplotypes
A total of 479 individuals, 42 genotypes were found to carry a different number and combination of 16 KIR genes (Figure 1), of which, 29 were predicted to be AB genotypes, 12 were predicted to be BB genotypes and one was AA genotype. 27 KIR genotypes were identified in Kazakhs, 30 KIR genotypes were identified in Uyghurs and 20 KIR genotypes were identified in Tibetans. The predominant genotype 1 (AA genotypes) occurred most frequently in Tibetans (52.7%, 118/224), Kazakhs (43.2%, 54/125) and Uyghurs (34.9%, 45/130). Variations in genotypes and haplotypes were significant between the Tibetan and Kazakh/Uyghur groups (Table 2): 3 displayed significant statistical difference between Tibetan and Kazakh populations, 6 displayed significant statistical difference between Tibetan and Uyghur populations, and 2 displayed significant statistical difference between Uyghur and Kazakh populations.
KIR gene frequency
16 KIR genes were detected in Tibetans, Kazakhs and Uyghurs. The carrier frequency (CF) was listed in Table 3. Four framework genes (2DL4, 3DL2, 3DL3 and pseudogene 3DP1) were present in all individuals. KIR2DL1, 2DL3, 3DL1 and 2DP1 were detected in high frequencies (>90%). Among them, pseudogene 2DP1 can be detected in all individuals of Uyghur. KIR2DL2, 2DL5, 2DS1, 2DS2, 2DS3, 2DS5 and 3DS1 were in low frequencies. In the Uyghurs, the frequencies of KIR2DS4-f and KIR2DS4-v were found to be 62.3% (81/130) and 60.8% (79/130); in the Kazakhs, the frequencies of KIR2DS4-f and KIR2DS4-v were found to be 64.8% (81/125) and 59.2% (74/125); in the Tibetans, the frequencies of KIR2DS4-f and 2DS4-v were found to be 86.6% (194/224) and 29.5% (66/224). Four KIR genes displayed significant statistical difference between the Tibetan and Kazakh populations, eight KIR genes showed significant statistical difference between the Tibetan and Uyghur populations, while no difference was observed between the Kazakh and Uyghur populations (Table 3).
Genetic relationships between populations
The PCA was generated depending on the basis of KIR gene frequencies of the three groups with other previously studied populations (as shown in Additional file 2: Table S1), PC1 and PC2 of 22 populations were given in Additional file 3: Table S2. Different geographic clusters of Asian and European were noted on the PCA plot. The most divergent were the clusters of African, Senegalese and Indian, as they did not cluster in a single group (Figure 2). Tibetans was plotted within the Asian group, the Uyghur was plotted within the European group, while the Kazakh was mapped between Asian and European groups. The neighbor-joining dendrogram was constructed, using the allelic frequency data of 12 KIR loci (KIR2DL1-3, 2DL5, 3DL1, 2DS1-5, 3DS1 and 2DP1) to show a relationship between the three groups and other previously reported populations  (Additional file 2: Table S1). In this dendrogram, there were four main clusters (Figure 3): European, African, Asian and Indian cluster. The Uyghur population was located in the European cluster, the Tibetan was located in the Asian cluster, while the Kazakh was located between European and Asian clusters. Genetic distance was shown in Figure 3.
Genetic relationships between three ethnic populations centromeric and telomeric clusters of KIR genes (CxTx)
We found that the three groups carried all four subsets of Bx genotypes (CxT4, C4Tx, C4T4, and CxTx, Figure 1). 55 with the expression of C4Tx genotypes carried three to five activating KIR genes, among them, 49.1% (27/55) carried three activating KIRs and 45.3% (25/55) carried five activating KIRs. 108 with the expression of CxT4 genotypes carried three to six activating KIR genes, among them, 73.2% (79/108) carried four activating KIRs. 1 with the expression of C4T4 genotypes carried six activating KIRs. 73 individuals displaying CxTx carried one to five activating KIRs, among them, 50.7% (37/73) carried two activating KIRs, while 37.0% (27/73) carried four activating KIRs (Figure 1). The Uyghurs (25.4%, 33/130) comprised higher frequencies of Bx genotypes with C4Tx than the Kazakhs (11.2%, 14/125) and Tibetans (3.6%, 8/224, Figure 4).
KIR genotypes show extensive variation in populations according to different geographical regions and different ethnic groups, due to the presence or absence of different KIRs in individual haplotypes and allelic polymorphism of KIR genes . In this study, Tibetans were detected to display a lower level of genotypic diversity, with the predominant genotypes ID1 (AA genotypes)  and ID2. The cumulative frequency of the top five genotypes 1–4, 6 was 84.8% (190/224), which showed that KIR genotype distribution in Tibetan populations was relatively concentrated. Meantime, Uyghurs and Kazakhs displayed a higher genotypic diversity, with the predominant genotypes ID1. Intriguingly, Uyghur displayed consistency with Caucasians on the level of KIR allelic frequency , which comprised high frequencies of Bx genotypes (absence of KIR3DS1-2DL5-2DS5-2DS1), while Tibetans shared a lower frequencies consistent with Northeast Asian (Han Chinese , Korean  and Japanese ). Our data on Kazakhs showed that the KIR allelic frequency is between that of Caucasians and Northeast Asians. Many studies on KIR gene diversity have been reported on Han populations in China [18–21]. Expanding our understanding of ethnic intermarriage and possibly random demographic factors for ethnic minorities could help us to determine the variation of KIR gene frequencies, which might be useful for future research on ethnicity-based diseases.
Activating KIR genes were found to show greater variability in frequencies than inhibitory KIR genes, which was similar to those previously reported in other populations . KIR2DS4 is the only activating gene in AA genotypes, with the other 5 activating genes (KIR2DS1-3, 2DS5 and 3DS1) being absent [5, 21]. In this study, KIR2DS4-v was found to have a higher frequency in Uyghurs and Kazakhs, which was consistent with the frequency in Caucasians . Tibetans had a lower KIR2DS4-v frequency when compared with Kazakhs and Uyghurs, more similar to Asian populations (Han Chinese, Japanese and Koreans) [21–23].
The PCA and the generated dendrogram aligned with the origins of ancestry distributing across vast geographic regions. In our data, the Uyghur was located in the European cluster, showing a similarity to the origins of European Caucasian ancestry. This was consistent with previously report that Uyghur population, settled in Xinjiang, was a population presenting a typical admixture of Eastern and Western anthropometric traits . The Kazakh was between Eastern Asian and European Caucasian groups, demonstrating an array of mixed anthropological features of Europeans and Asians. This could help to further understand the origins of ancestry, the pattern of human migration, the intermarriage between ethnicities and the genetic distances between different populations or ethnic groups.
The KIR gene frequencies of the Tibetan and Uyghur populations were consistent with the geographic and ethnic distribution, while the Kazakhs showed unique genetic content and allelic polymorphisms. Better understanding the ethnic intermarriage and possibly random demographic factors could help determine the factors in the variation of KIR gene frequencies, which could be useful for future research on ethnicity-based diseases.
Ashouri E, Farjadian S, Reed EF, Ghaderi A, Rajalingam R (2009) KIR gene content diversity in four Iranian populations. Immunogenetics 61:483–492
Schleinitz N, Hamidou M, Vély F (2006) NK cells: new insights on physiology and clinical implication in diseases. Rev Med Interne 27:465–472
Alter G, Heckerman D, Schneidewind A, Fadda L, Kadie CM (2011) HIV-1 adaptation to NK-cell-mediated immune pressure. Nature 476:96–100
Wilson MJ, Torkar M, Haude A, Milne S, Jones T (2000) Plasticity in the organization and sequences of human KIR/ILT gene families. Proc Natl Acad Sci USA 97:4778–4783
Faure M, Long EO (2002) KIR2DL4 (CD158d), an NK cell-activating receptor with inhibitory potential. J Immunol 168:6208–6214
Parham P (2005) MHC class Imolecules and KIRs in human history, health and survival. Nat Rev Immunol 5:201–214
Single RM, Martin MP, Gao X, Meyer D, Yeager M, Kidd JR et al (2007) Global diversity and evidence for coevolution of KIR and HLA. Nat Genet 39:1114–1119
McQueen KL, Dorighi KM, Guethlein LA, Wong R, Sanjanwala B, Parham P (2007) Donor–recipient combinations of group A and B KIR haplotypes and HLA class I ligand affect the outcome of HLA-matched, sibling donor hematopoietic cell transplantation. Hum Immunol 68:309–323
Heather AH, Maria P, Glenda EG, Louise K, Caroline TT (2013) KIR2DS4 allelic variants: differential effects on in utero and intrapartum HIV-1 mother-to-child transmission. Clin Immunol 149:498–508
Middleton D, Gonzalez A, Gilmore PM (2007) Studies on the expression of the deleted KIR2DS4*003 gene product and distribution of KIR2DS4 deleted and nondeleted versions in different populations. Hum Immunol 68:128–134
Maxwell LD, Williams F, Gilmore PA (2004) Investigation of killer cell immunoglobulin-like receptor gene diversity: II. KIR2DS4. Hum Immunol 65:613–621
Bao X, Hou L, Sun A (2010) An allelic typing method for 2DS4 variant used in study of haplotypes of killer cell immunoglobulin-like receptor gene. Int J Lab Hematol 32:625–632
Gllez-Hawkins GM, Franck AE, Li X (2011) Expression of activating KIR2DS2 and KIR2DS4 genes after hematopoietic cell transplantation: relevance to cytomegalovirus infection. Biol Blood Marrow Transplant 17:1662–1672
Merino A, Malhotra R, Morton M, Mulenga J (2011) Impact of a functional KIR2DS4 allele on heterosexual HIV-1 transmission among discordant Zambian couples. J Infect Dis 203:487–495
The allele frequency net database. Available at: http://www.allelefrequencies.net/
Badwe R, Joshi N, Carrington M (2008) Comparison of the rapidly evolving KIR locus in Parsis and natives of India. Immunogenetics 60:121–129
Wang HD, Zhang FX, Shen CM (2012) Genetic polymorphism analysis of killer cell immunoglobulin-like receptor genes in the Chinese Uygur population. Mol Biol 39:3017–3028
Zhan XT, Wang S, Wang J (2012) Relationship between KIR and HLA ligands in Han population of Sichuan Marrow Donor Program. Chin J Blood Transfus 25:327–332
Zhu FM, Jiang K (2005) Zhejiang Han population of killer cell immunoglobulin-like receptor gene polymorphism. Chin J Exp Hematol 13:1109–1112
Su PC, Yang TH, Zou H (2008) Polymorph ism of killer cell immunoglobulin-like receptors gene in Yunnan Han population, China. Chin J Blood Transfus 21:268–272
Bao X, Wang M, Zhou H, Wu X, Yang L (2013) Characterization of killer cell immunoglobulin-like receptor (KIR) genotypes and haplotypes in Chinese Han population. Tissue Antigens 82:327–337
Whang DH, Park H, Yoon JA (2005) ParkMH: haplotype analysis of killer cell immunoglobulin-like receptor genes in 77 Korean families. Hum Immunol 66:146–154
Yawata M, Yawata N, McQueen KL, Cheng NW, Guethlein LA, Rajalingam R et al (2002) Predominance of group A KIR haplotypes in Japanese associated with diverse NK cell repertoires of KIR expression. Immunogenetics 54:543–550
National Bureau of Statistics of China: Chinese Statistical Yearbook 2000. Beijing, China, Statistics Publishing Group 2001
Wang JX, Liu J, Yao FZ, Wen GX, Li JL, Huang Y et al (2013) Prevalence, incidence, and residual risks for transfusion transmitted HIV-1/2 infection among Chinese blood donors. Transfusion 53(6):1240–1249
Tabulation on the 2010 population census of the Peoples Republic of China. Available at: http://www.stats.gov.cn/tjsj/pcsj/rkpc/6rp/indexch.htm
Xu SH, Huang W, Qian J, Jin L (2008) Analysis of genomic admixture in Uyghur and its implication in mapping strategy. Am J Hum Genet 82:883–894
Campbell KS, Martin MP, Carrington M (2010) KIR genotyping by multiplex PCR-SSP//Campbell, SK. Natural killer cell protocols:methods in molecular biology, vol 612. Humana Press, New York, pp 365–375
Felsenstein J (1989) PHYLIP—Phylogeny Inference Package (Version 3.2). Cladistics 5:164–166
Saitou N, Nei M (1987) The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 4:406–425
QZ and JW designed and performed the experiments, analyzed and interpreted the data, and drafted the manuscript; ZH and XL collected and informed consents of all samples used in this study; SM, SH and GB participated in the experiments and helped analyzed data; FM designed the experiments and wrote the paper. All authors read and approved the final manuscript.
We thank Dr. Bo Chen and Ms. Ka Yi Li at Institute of Blood Transfusion, CAMS and PUMC, for their generous help in proofreading the manuscript and providing constructive advice. This work was sponsored by National Natural Science Foundation of China (NSFC 81001341) to Jue Wang.
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Qiongxiu Zhou and Jue Wang are co-first authors
Figure S1. Map of China showing the city of three study populations. DNA samples of the Kazakh and Uyghur ethnic minority populations were collected from the Xinjiang autonomous region (Urumqi) of Northwest China, the Tibetan ethnic minority populations were collected from the Tibet autonomous region (Lhasa) of Southwest China.
Table S1. KIR diversity in three ethnic minority populations in China with other populations.
Table S2. PC1 and PC2 of 22 populations.
About this article
Cite this article
Zhou, Q., Wang, J., He, Z. et al. KIR diversity in three ethnic minority populations in China. J Transl Med 13, 221 (2015). https://doi.org/10.1186/s12967-015-0544-7