Neda Nemat-Gorgani

Genetic and Environmental Determinants of Human NK Cell Diversity Revealed by Mass Cytometry

Both genetics and environment contribute to human NK cell diversity. NK Cell Nature Versus Nurture Natural killer (NK) cells were first discovered because of their ability to kill tumor cells without any previous exposure. However, this population is actually quite heterogeneous: Different subgroups of NK cells express different combinations of activating and inhibiting receptors that govern their specificity. Now, Horowitz et al. use mass cytometry to examine NK cell diversity in humans. The authors examined 35 parameters simultaneously in 5 sets of monozygotic twins as well as 12 unrelated donors. They found up to 30,000 phenotypic NK cell populations in a given individual. What’s more, by comparing the twins versus unrelated donors, they determined that although genetics primarily determined inhibitory receptor expression, activating receptors were controlled by the environment. These data suggest that inhibitory receptors may contribute more to NK cell self-tolerance, whereas activating receptors may guide response to pathogens and tumors. Natural killer (NK) cells play critical roles in immune defense and reproduction, yet remain the most poorly understood major lymphocyte population. Because their activation is controlled by a variety of combinatorially expressed activating and inhibitory receptors, NK cell diversity and function are closely linked. To provide an unprecedented understanding of NK cell repertoire diversity, we used mass cytometry to simultaneously analyze 37 parameters, including 28 NK cell receptors, on peripheral blood NK cells from 5 sets of monozygotic twins and 12 unrelated donors of defined human leukocyte antigen (HLA) and killer cell immunoglobulin-like receptor (KIR) genotype. This analysis revealed a remarkable degree of NK cell diversity, with an estimated 6000 to 30,000 phenotypic populations within an individual and >100,000 phenotypes in the donor panel. Genetics largely determined inhibitory receptor expression, whereas activation receptor expression was heavily environmentally influenced. Therefore, NK cells may maintain self-tolerance through strictly regulated expression of inhibitory receptors while using adaptable expression patterns of activating and costimulatory receptors to respond to pathogens and tumors. These findings further suggest the possibility that discrete NK cell subpopulations could be harnessed for immunotherapeutic strategies in the settings of infection, reproduction, and transplantation.

Co-evolution of Human Leukocyte Antigen (HLA) Class I Ligands with Killer-Cell Immunoglobulin-Like Receptors (KIR) in a Genetically Diverse Population of Sub-Saharan Africans

Interactions between HLA class I molecules and killer-cell immunoglobulin-like receptors (KIR) control natural killer cell (NK) functions in immunity and reproduction. Encoded by genes on different chromosomes, these polymorphic ligands and receptors correlate highly with disease resistance and susceptibility. Although studied at low-resolution in many populations, high-resolution analysis of combinatorial diversity of HLA class I and KIR is limited to Asian and Amerindian populations with low genetic diversity. At the other end of the spectrum is the West African population investigated here: we studied 235 individuals, including 104 mother-child pairs, from the Ga-Adangbe of Ghana. This population has a rich diversity of 175 KIR variants forming 208 KIR haplotypes, and 81 HLA-A, -B and -C variants forming 190 HLA class I haplotypes. Each individual we studied has a unique compound genotype of HLA class I and KIR, forming 1–14 functional ligand-receptor interactions. Maintaining this exceptionally high polymorphism is balancing selection. The centromeric region of the KIR locus, encoding HLA-C receptors, is highly diverse whereas the telomeric region encoding Bw4-specific KIR3DL1, lacks diversity in Africans. Present in the Ga-Adangbe are high frequencies of Bw4-bearing HLA-B*53:01 and Bw4-lacking HLA-B*35:01, which otherwise are identical. Balancing selection at key residues maintains numerous HLA-B allotypes having and lacking Bw4, and also those of stronger and weaker interaction with LILRB1, a KIR-related receptor. Correspondingly, there is a balance at key residues of KIR3DL1 that modulate its level of cell-surface expression. Thus, capacity to interact with NK cells synergizes with peptide binding diversity to drive HLA-B allele frequency distribution. These features of KIR and HLA are consistent with ongoing co-evolution and selection imposed by a pathogen endemic to West Africa. Because of the prevalence of malaria in the Ga-Adangbe and previous associations of cerebral malaria with HLA-B*53:01 and KIR, Plasmodium falciparum is a candidate pathogen.

A KIR B centromeric region present in Africans but not Europeans protects pregnant women from pre-eclampsia.

Significance Pre-eclampsia is especially common in women of African ancestry and a major cause of maternal death. The killer-cell immunoglobulin-like receptor (KIR) genes that we analyzed are expressed by natural killer cells—immune cells that populate the uterus and are essential for successful pregnancy. KIR proteins bind HLA ligands on the implanting placental trophoblast cells. African and European women share similar risk associations for pre-eclampsia, but protection is associated with different KIR genes. African women are protected by a combination of KIR B haplotype genes that is present almost exclusively in Africans. This study emphasizes the importance of studying diseases in Africans, where the KIR/HLA genetic system is at its most diverse and maternal mortality rates are the highest in the world. In sub-Saharan Africans, maternal mortality is unacceptably high, with >400 deaths per 100,000 births compared with <10 deaths per 100,000 births in Europeans. One-third of the deaths are caused by pre-eclampsia, a syndrome arising from defective placentation. Controlling placentation are maternal natural killer (NK) cells that use killer-cell immunoglobulin-like receptor (KIR) to recognize the fetal HLA-C molecules on invading trophoblast. We analyzed genetic polymorphisms of maternal KIR and fetal HLA-C in 484 normal and 254 pre-eclamptic pregnancies at Mulago Hospital, Kampala, Uganda. The combination of maternal KIR AA genotypes and fetal HLA-C alleles encoding the C2 epitope associates with pre-eclampsia [P = 0.0318, odds ratio (OR) = 1.49]. The KIR genes associated with protection are located in centromeric KIR B regions that are unique to sub-Saharan African populations and contain the KIR2DS5 and KIR2DL1 genes (P = 0.0095, OR = 0.59). By contrast, telomeric KIR B genes protect Europeans against pre-eclampsia. Thus, different KIR B regions protect sub-Saharan Africans and Europeans from pre-eclampsia, whereas in both populations, the KIR AA genotype is a risk factor for the syndrome. These results emphasize the importance of undertaking genetic studies of pregnancy disorders in African populations with the potential to provide biological insights not available from studies restricted to European populations.

Defining KIR and HLA Class I Genotypes at Highest Resolution via High-Throughput Sequencing.

The physiological functions of natural killer (NK) cells in human immunity and reproduction depend upon diverse interactions between killer cell immunoglobulin-like receptors (KIRs) and their HLA class I ligands: HLA-A, HLA-B, and HLA-C. The genomic regions containing the KIR and HLA class I genes are unlinked, structurally complex, and highly polymorphic. They are also strongly associated with a wide spectrum of diseases, including infections, autoimmune disorders, cancers, and pregnancy disorders, as well as the efficacy of transplantation and other immunotherapies. To facilitate study of these extraordinary genes, we developed a method that captures, sequences, and analyzes the 13 KIR genes and HLA-A, HLA-B, and HLA-C from genomic DNA. We also devised a bioinformatics pipeline that attributes sequencing reads to specific KIR genes, determines copy number by read depth, and calls high-resolution genotypes for each KIR gene. We validated this method by using DNA from well-characterized cell lines, comparing it to established methods of HLA and KIR genotyping, and determining KIR genotypes from 1000 Genomes sequence data. This identified 116 previously uncharacterized KIR alleles, which were all demonstrated to be authentic by sequencing from source DNA via standard methods. Analysis of just two KIR genes showed that 22% of the 1000 Genomes individuals have a previously uncharacterized allele or a structural variant. The method we describe is suited to the large-scale analyses that are needed for characterizing human populations and defining the precise HLA and KIR factors associated with disease. The methods are applicable to other highly polymorphic genes.

Class I HLA haplotypes form two schools that educate NK cells in different ways

Individuals heterozygous for HLA haplotypes have increased NK cell diversity because of distinct education pathways. NK cell immunity gets schooled Population-level genetic analysis can give us clues as to which factors guide immune evolution. Now, Horowitz et al. have applied this analysis to polymorphisms in HLA that affect NK cell education. They find that the HLA-B haplotype −21M that delivers functional peptides to the conserved CD94/NKG2A receptor rarely encodes ligands for the more diverse killer cell immunoglobulin-like receptors (KIRs), in contrast to the haplotype −21T, which does not deliver functional peptides. Individuals homo- or heterozygous for −21M are more likely to have more diverse CD94/NKG2A+ NK cells, suggesting that these HLA haplotypes may have specialized to either the KIR or CD94/NKG2A school through complementary coevolution. Natural killer (NK) cells are lymphocytes that have vital functions in innate and adaptive immunity, as well as placental reproduction. Polymorphic human leukocyte antigen (HLA) class I educates NK cells through interactions with killer cell immunoglobulin-like receptors (KIRs) and by supplying peptides that bind HLA-E to form ligands for CD94/NKG2A receptors on NK cells. HLA-B dimorphism in the leader peptide modulates this latter function: −21methionine (−21M) delivers functional peptides, but −21threonine (−21T) does not. Genetic analysis of human populations worldwide showed that haplotypes with −21M HLA-B rarely encoded the KIR ligands Bw4+HLA-B and C2+HLA-C. Thus, there are two fundamental forms of HLA haplotype: one preferentially supplying CD94/NKG2A ligands and the other preferentially supplying KIR ligands. This −21 HLA-B dimorphism divides the human population into three groups: M/M, M/T, and T/T. Mass cytometry and assays of immune function demonstrated that M/M and M/T individuals have CD94/NKG2A+ NK cells that are better educated, phenotypically more diverse, and functionally more potent than those in T/T individuals. The KIR school of NK cell education evolved in the context of the much older CD94/NKG2A school, and this complementary coevolution may have facilitated the specialization of HLA haplotypes toward one school or the other.Natural killer (NK) cells are lymphocytes that have vital functions in innate and adaptive immunity, as well as placental reproduction. Polymorphic human leukocyte antigen (HLA) class I educates NK cells through interactions with killer cell immunoglobulin-like receptors (KIRs) and by supplying peptides that bind HLA-E to form ligands for CD94/NKG2A receptors on NK cells. HLA-B dimorphism in the leader peptide modulates this latter function: -21methionine (-21M) delivers functional peptides, but -21threonine (-21T) does not. Genetic analysis of human populations worldwide showed that haplotypes with -21M HLA-B rarely encoded the KIR ligands Bw4+HLA-B and C2+HLA-C. Thus, there are two fundamental forms of HLA haplotype: one preferentially supplying CD94/NKG2A ligands and the other preferentially supplying KIR ligands. This -21 HLA-B dimorphism divides the human population into three groups: M/M, M/T, and T/T. Mass cytometry and assays of immune function demonstrated that M/M and M/T individuals have CD94/NKG2A+ NK cells that are better educated, phenotypically more diverse, and functionally more potent than those in T/T individuals. The KIR school of NK cell education evolved in the context of the much older CD94/NKG2A school, and this complementary coevolution may have facilitated the specialization of HLA haplotypes toward one school or the other.

Polymorphic HLA-C Receptors Balance the Functional Characteristics of KIR Haplotypes

The human killer cell Ig-like receptor (KIR) locus comprises two groups of KIR haplotypes, termed A and B. These are present in all human populations but with different relative frequencies, suggesting they have different functional properties that underlie their balancing selection. We studied the genomic organization and functional properties of the alleles of the inhibitory and activating HLA-C receptors encoded by KIR haplotypes. Because every HLA-C allotype functions as a ligand for KIR, the interactions between KIR and HLA-C dominate the HLA class I–mediated regulation of human NK cells. The C2 epitope is recognized by inhibitory KIR2DL1 and activating KIR2DS1, whereas the C1 epitope is recognized by inhibitory KIR2DL2 and KIR2DL3. This study shows that the KIR2DL1, KIR2DS1, and KIR2DL2/3 alleles form distinctive phylogenetic clades that associate with specific KIR haplotypes. KIR A haplotypes are characterized by KIR2DL1 alleles that encode strong inhibitory C2 receptors and KIR2DL3 alleles encoding weak inhibitory C1 receptors. In striking contrast, KIR B haplotypes are characterized by KIR2DL1 alleles that encode weak inhibitory C2 receptors and KIR2DL2 alleles encoding strong inhibitory C1 receptors. The wide-ranging properties of KIR allotypes arise from substitutions throughout the KIR molecule. Such substitutions can influence cell surface expression, as well as the avidity and specificity for HLA-C ligands. Consistent with the crucial role of inhibitory HLA-C receptors in self-recognition, as well as NK cell education and response, most KIR haplotypes have both a functional C1 and C2 receptor, despite the considerable variation that occurs in ligand recognition and surface expression.

Regulation of Adaptive NK Cells and CD8 T Cells by HLA-C Correlates with Allogeneic Hematopoietic Cell Transplantation and with Cytomegalovirus Reactivation

Mass cytometry was used to investigate the effect of CMV reactivation on lymphocyte reconstitution in hematopoietic cell transplant patients. For eight transplant recipients (four CMV negative and four CMV positive), we studied PBMCs obtained 6 mo after unrelated donor hematopoietic cell transplantation (HCT). Forty cell-surface markers, distinguishing all major leukocyte populations in PBMC, were analyzed with mass cytometry. This group included 34 NK cell markers. Compared with healthy controls, transplant recipients had higher HLA-C expression on CD56−CD16+ NK cells, B cells, CD33bright myeloid cells, and CD4CD8 T cells. The increase in HLA-C expression was greater for CMV-positive HCT recipients than for CMV negative recipients. Present in CMV-positive HCT recipients, but not in CMV-negative HCT recipients or controls, is a population of killer cell Ig-like receptor (KIR)-expressing CD8 T cells not previously described. These CD8 T cells coexpress CD56, CD57, and NKG2C. The HCT recipients also have a population of CD57+NKG2A+ NK cells that preferentially express KIR2DL1. An inverse correlation was observed between the frequencies of CD57+NKG2C+ NK cells and CD57+NKG2A+ NK cells. Although CD57+NKG2A+ NK cells are less abundant in CMV-positive recipients, their phenotype is of a more activated cell than the CD57+NKG2A+ NK cells of controls and CMV-negative HCT recipients. These data demonstrate that HCT and CMV reactivation are associated with an increased expression of HLA-C. This could influence NK cell education during lymphocyte reconstitution. The increased inhibitory KIR expression by proliferating CMV-specific CD8 T cells suggests regulatory interactions between HLA-C and KIR might promote Graft-versus-Leukemia effects following transplantation.

Exome capture from saliva produces high quality genomic and metagenomic data

BackgroundTargeted capture of genomic regions reduces sequencing cost while generating higher coverage by allowing biomedical researchers to focus on specific loci of interest, such as exons. Targeted capture also has the potential to facilitate the generation of genomic data from DNA collected via saliva or buccal cells. DNA samples derived from these cell types tend to have a lower human DNA yield, may be degraded from age and/or have contamination from bacteria or other ambient oral microbiota. However, thousands of samples have been previously collected from these cell types, and saliva collection has the advantage that it is a non-invasive and appropriate for a wide variety of research.ResultsWe demonstrate successful enrichment and sequencing of 15 South African KhoeSan exomes and 2 full genomes with samples initially derived from saliva. The expanded exome dataset enables us to characterize genetic diversity free from ascertainment bias for multiple KhoeSan populations, including new exome data from six HGDP Namibian San, revealing substantial population structure across the Kalahari Desert region. Additionally, we discover and independently verify thirty-one previously unknown KIR alleles using methods we developed to accurately map and call the highly polymorphic HLA and KIR loci from exome capture data. Finally, we show that exome capture of saliva-derived DNA yields sufficient non-human sequences to characterize oral microbial communities, including detection of bacteria linked to oral disease (e.g. Prevotella melaninogenica). For comparison, two samples were sequenced using standard full genome library preparation without exome capture and we found no systematic bias of metagenomic information between exome-captured and non-captured data.ConclusionsDNA from human saliva samples, collected and extracted using standard procedures, can be used to successfully sequence high quality human exomes, and metagenomic data can be derived from non-human reads. We find that individuals from the Kalahari carry a higher oral pathogenic microbial load than samples surveyed in the Human Microbiome Project. Additionally, rare variants present in the exomes suggest strong population structure across different KhoeSan populations.

Loss and Gain of Natural Killer Cell Receptor Function in an African Hunter-Gatherer Population.

Modulating natural killer cell functions in human immunity and reproduction are diverse interactions between the killer cell immunoglobulin-like receptors (KIR) of Natural Killer (NK) cells and HLA class I ligands on the surface of tissue cells. Dominant interactions are between KIR2DL1 and the C2 epitope of HLA-C and between KIR2DL2/3 and the C1 epitope of HLA-C. KhoeSan hunter-gatherers of Southern Africa represent the earliest population divergence known and are the most genetically diverse indigenous people, qualities reflected in their KIR and HLA genes. Of the ten KhoeSan KIR2DL1 alleles, KIR2DL1*022 and KIR2DL1*026 likely originated in the KhoeSan, and later were transmitted at low frequency to the neighboring Zulus through gene flow. These alleles arose by point mutation from other KhoeSan KIR2DL1 alleles that are more widespread globally. Mutation of KIR2DL1*001 gave rise to KIR2DL1*022, causing loss of C2 recognition and gain of C1 recognition. This makes KIR2DL1*022 a more avid and specific C1 receptor than any KIR2DL2/3 allotype. Mutation of KIR2DL1*012 gave rise to KIR2DL1*026, causing premature termination of translation at the end of the transmembrane domain. This makes KIR2DL1*026 a membrane-associated receptor that lacks both a cytoplasmic tail and signaling function. At higher frequencies than their parental allotypes, the combined effect of the KhoeSan-specific KIR2DL1*022 and KIR2DL1*026 is to reduce the frequency of strong inhibitory C2 receptors and increase the frequency of strong inhibitory C1 receptors. Because interaction of KIR2DL1 with C2 is associated with risk of pregnancy disorder, these functional changes are potentially advantageous. Whereas all other KhoeSan KIR2DL1 alleles are present on a wide diversity of centromeric KIR haplotypes, KIR2DL1*026 is present on a single KIR haplotype and KIR2DL1*022 is present on two very similar haplotypes. The high linkage disequilibrium across their haplotypes is consistent with a recent emergence for these KIR2DL1 alleles that have distinctive functions.

Two alternate strategies for innate immunity to Epstein-Barr virus: One using NK cells and the other NK cells and γδ T cells

Most humans become infected with Epstein–Barr virus (EBV), which then persists for life. Infrequently, EBV infection causes infectious mononucleosis (IM) or Burkitt lymphoma (BL). Type I EBV infection, particularly type I BL, stimulates strong responses of innate immune cells. Humans respond to EBV in two alternative ways. Of 24 individuals studied, 13 made strong NK and &ggr;&dgr; T cell responses, whereas 11 made feeble &ggr;&dgr; T cell responses but stronger NK cell responses. The difference does not correlate with sex, HLA type, or previous exposure to EBV or cytomegalovirus. Cohorts of EBV+ children and pediatric IM patients include both group 1 individuals, with high numbers of &ggr;&dgr; T cells, and group 2 individuals, with low numbers. The even balance of groups 1 and 2 in the human population points to both forms of innate immune response to EBV having benefit for human survival. Correlating these distinctive responses with the progress of EBV infection might facilitate the management of EBV-mediated disease.