KIR diversity in three ethnic minority populations in China
- Qiongxiu Zhou†1,
- Jue Wang†1,
- Zhi He2,
- Xiaojuan Li3,
- Song Mao4,
- Shu Huang1,
- Guohui Bian1 and
- Feng Ma1Email author
© Zhou et al. 2015
Received: 12 January 2015
Accepted: 21 May 2015
Published: 11 July 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
Multiple PCR primers combination pattern
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
Comparison of genotypes, haplotypes and linkage groups in Tibetan, Kazakh and Uyghur populations
N = 224
N = 125
N = 130
Tibetan vs Kazakh
Tibetan vs Uyghur
Kazakh vs Uyghur
KIR gene frequency
Comparison of carrier frequency of KIR genes in three ethnic minority populations in China
N = 224
N = 125
N = 130
Tibetan vs Kazakh
Tibetan vs Uyghur
Kazak vs Uyghur
A haplotype associated
B haplotype associated
Genetic relationships between populations
Genetic relationships between three ethnic populations centromeric and telomeric clusters of KIR genes (CxTx)
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.
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.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ashouri E, Farjadian S, Reed EF, Ghaderi A, Rajalingam R (2009) KIR gene content diversity in four Iranian populations. Immunogenetics 61:483–492PubMed CentralPubMedView ArticleGoogle Scholar
- Schleinitz N, Hamidou M, Vély F (2006) NK cells: new insights on physiology and clinical implication in diseases. Rev Med Interne 27:465–472PubMedView ArticleGoogle Scholar
- Alter G, Heckerman D, Schneidewind A, Fadda L, Kadie CM (2011) HIV-1 adaptation to NK-cell-mediated immune pressure. Nature 476:96–100PubMed CentralPubMedView ArticleGoogle Scholar
- 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–4783PubMed CentralPubMedView ArticleGoogle Scholar
- Faure M, Long EO (2002) KIR2DL4 (CD158d), an NK cell-activating receptor with inhibitory potential. J Immunol 168:6208–6214PubMedView ArticleGoogle Scholar
- Parham P (2005) MHC class Imolecules and KIRs in human history, health and survival. Nat Rev Immunol 5:201–214PubMedView ArticleGoogle Scholar
- 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–1119PubMedView ArticleGoogle Scholar
- 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–323PubMed CentralPubMedView ArticleGoogle Scholar
- 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–508View ArticleGoogle Scholar
- 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–134PubMedView ArticleGoogle Scholar
- Maxwell LD, Williams F, Gilmore PA (2004) Investigation of killer cell immunoglobulin-like receptor gene diversity: II. KIR2DS4. Hum Immunol 65:613–621PubMedView ArticleGoogle Scholar
- 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–632PubMedView ArticleGoogle Scholar
- 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–1672View ArticleGoogle Scholar
- 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–495PubMed CentralPubMedView ArticleGoogle Scholar
- 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–129PubMedView ArticleGoogle Scholar
- 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–3028Google Scholar
- 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–332Google Scholar
- Zhu FM, Jiang K (2005) Zhejiang Han population of killer cell immunoglobulin-like receptor gene polymorphism. Chin J Exp Hematol 13:1109–1112Google Scholar
- 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–272Google Scholar
- 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–337PubMedView ArticleGoogle Scholar
- 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–154PubMedView ArticleGoogle Scholar
- 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–550PubMedView ArticleGoogle Scholar
- National Bureau of Statistics of China: Chinese Statistical Yearbook 2000. Beijing, China, Statistics Publishing Group 2001Google Scholar
- 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–1249PubMed CentralPubMedView ArticleGoogle Scholar
- 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–894PubMed CentralPubMedView ArticleGoogle Scholar
- 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–375Google Scholar
- Felsenstein J (1989) PHYLIP—Phylogeny Inference Package (Version 3.2). Cladistics 5:164–166Google Scholar
- Saitou N, Nei M (1987) The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 4:406–425PubMedGoogle Scholar