Following Transplantation for Acute Myelogenous Leukemia, Donor KIR Cen B02 Better Protects against Relapse than KIR Cen B01.

In the treatment of acute myelogenous leukemia with allogeneic hematopoietic cell transplantation, we previously demonstrated that there is a greater protection from relapse of leukemia when the hematopoietic cell transplantation donor has either the Cen B/B KIR genotype or a genotype having two or more KIR B gene segments. In those earlier analyses, KIR genotyping could only be assessed at the low resolution of gene presence or absence. To give the analysis greater depth, we developed high-resolution KIR sequence-based typing that defines all the KIR alleles and distinguishes the expressed alleles from those that are not expressed. We now describe and analyze high-resolution KIR genotypes for 890 donors of this human transplant cohort. Cen B01 and Cen B02 are the common CenB haplotypes, with Cen B02 having evolved from Cen B01 by deletion of the KIR2DL5, 2DS3/5, 2DP1, and 2DL1 genes. We observed a consistent trend for Cen B02 to provide stronger protection against relapse than Cen B01 This correlation indicates that protection depends on the donor having inhibitory KIR2DL2 and/or activating KIR2DS2, and is enhanced by the donor lacking inhibitory KIR2DL1, 2DL3, and 3DL1. High-resolution KIR typing has allowed us to compare the strength of the interactions between the recipient's HLA class I and the KIR expressed by the donor-derived NK cells and T cells, but no clinically significant interactions were observed. The trend observed between donor Cen B02 and reduced relapse of leukemia points to the value of studying ever larger transplant cohorts.

Adaptive Admixture of HLA Class I Allotypes Enhanced Genetically Determined Strength of Natural Killer Cells in East Asians.

Human natural killer (NK) cells are essential for controlling infection, cancer, and fetal development. NK cell functions are modulated by interactions between polymorphic inhibitory killer cell immunoglobulin-like receptors (KIR) and polymorphic HLA-A, -B, and -C ligands expressed on tissue cells. All HLA-C alleles encode a KIR ligand and contribute to reproduction and immunity. In contrast, only some HLA-A and -B alleles encode KIR ligands and they focus on immunity. By high-resolution analysis of KIR and HLA-A, -B, and -C genes, we show that the Chinese Southern Han (CHS) are significantly enriched for interactions between inhibitory KIR and HLA-A and -B. This enrichment has had substantial input through population admixture with neighboring populations, who contributed HLA class I haplotypes expressing the KIR ligands B*46:01 and B*58:01, which subsequently rose to high frequency by natural selection. Consequently, over 80% of Southern Han HLA haplotypes encode more than one KIR ligand. Complementing the high number of KIR ligands, the CHS KIR locus combines a high frequency of genes expressing potent inhibitory KIR, with a low frequency of those expressing activating KIR. The Southern Han centromeric KIR region encodes strong, conserved, inhibitory HLA-C-specific receptors, and the telomeric region provides a high number and diversity of inhibitory HLA-A and -B-specific receptors. In all these characteristics, the CHS represent other East Asians, whose NK cell repertoires are thus enhanced in quantity, diversity, and effector strength, likely augmenting resistance to endemic viral infections.

Killer Cell Immunoglobulin-like Receptor Variants Are Associated with Protection from Symptoms Associated with More Severe Course in Parkinson Disease.

Immune dysfunction plays a role in the development of Parkinson disease (PD). NK cells regulate immune functions and are modulated by killer cell immunoglobulin-like receptors (KIR). KIR are expressed on the surface of NK cells and interact with HLA class I ligands on the surface of all nucleated cells. We investigated KIR-allelic polymorphism to interrogate the role of NK cells in PD. We sequenced KIR genes from 1314 PD patients and 1978 controls using next-generation methods and identified KIR genotypes using custom bioinformatics. We examined associations of KIR with PD susceptibility and disease features, including age at disease onset and clinical symptoms. We identified two KIR3DL1 alleles encoding highly expressed inhibitory receptors associated with protection from PD clinical features in the presence of their cognate ligand: KIR3DL1*015/HLA-Bw4 from rigidity (p c = 0.02, odds ratio [OR] = 0.39, 95% confidence interval [CI] 0.23-0.69) and KIR3DL1*002/HLA-Bw4i from gait difficulties (p c = 0.05, OR = 0.62, 95% CI 0.44-0.88), as well as composite symptoms associated with more severe disease. We also developed a KIR3DL1/HLA interaction strength metric and found that weak KIR3DL1/HLA interactions were associated with rigidity (pc = 0.05, OR = 9.73, 95% CI 2.13-172.5). Highly expressed KIR3DL1 variants protect against more debilitating symptoms of PD, strongly implying a role of NK cells in PD progression and manifestation.

A specific amino acid motif of HLA-DRB1 mediates risk and interacts with smoking history in Parkinson's disease.

Parkinson's disease (PD) is a neurodegenerative disease in which genetic risk has been mapped to HLA, but precise allelic associations have been difficult to infer due to limitations in genotyping methodology. Mapping PD risk at highest possible resolution, we performed sequencing of 11 HLA genes in 1,597 PD cases and 1,606 controls. We found that susceptibility to PD can be explained by a specific combination of amino acids at positions 70-74 on the HLA-DRB1 molecule. Previously identified as the primary risk factor in rheumatoid arthritis and referred to as the "shared epitope" (SE), the residues Q/R-K/R-R-A-A at positions 70-74 in combination with valine at position 11 (11-V) is highly protective in PD, while risk is attributable to the identical epitope in the absence of 11-V. Notably, these effects are modified by history of cigarette smoking, with a strong protective effect mediated by a positive history of smoking in combination with the SE and 11-V (P = 10-4; odds ratio, 0.51; 95% confidence interval, 0.36-0.72) and risk attributable to never smoking in combination with the SE without 11-V (P = 0.01; odds ratio, 1.51; 95% confidence interval, 1.08-2.12). The association of specific combinations of amino acids that participate in critical peptide-binding pockets of the HLA class II molecule implicates antigen presentation in PD pathogenesis and provides further support for genetic control of neuroinflammation in disease. The interaction of HLA-DRB1 with smoking history in disease predisposition, along with predicted patterns of peptide binding to HLA, provide a molecular model that explains the unique epidemiology of smoking in PD.

Diversity of KIR, HLA Class I, and Their Interactions in Seven Populations of Sub-Saharan Africans.

HLA class I and KIR sequences were determined for Dogon, Fulani, and Baka populations of western Africa, Mbuti of central Africa, and Datooga, Iraqw, and Hadza of eastern Africa. Study of 162 individuals identified 134 HLA class I alleles (41 HLA-A, 60 HLA-B, and 33 HLA-C). Common to all populations are three HLA-C alleles (C1+C*07:01, C1+C*07:02, and C2+C*06:02) but no HLA-A or -B Unexpectedly, no novel HLA class I was identified in these previously unstudied and anthropologically distinctive populations. In contrast, of 227 KIR detected, 22 are present in all seven populations and 28 are novel. A high diversity of HLA A-C-B haplotypes was observed. In six populations, most haplotypes are represented just once. But in the Hadza, a majority of haplotypes occur more than once, with 2 having high frequencies and 10 having intermediate frequencies. The centromeric (cen) part of the KIR locus exhibits an even balance between cenA and cenB in all seven populations. The telomeric (tel) part has an even balance of telA to telB in East Africa, but this changes across the continent to where telB is vestigial in West Africa. All four KIR ligands (A3/11, Bw4, C1, and C2) are present in six of the populations. HLA haplotypes of the Iraqw and Hadza encode two KIR ligands, whereas the other populations have an even balance between haplotypes encoding one and two KIR ligands. Individuals in these African populations have a mean of 6.8-8.4 different interactions between KIR and HLA class I, compared with 2.9-6.5 for non-Africans.

Sequences of 95 human MHC haplotypes reveal extreme coding variation in genes other than highly polymorphic HLA class I and II.

The most polymorphic part of the human genome, the MHC, encodes over 160 proteins of diverse function. Half of them, including the HLA class I and II genes, are directly involved in immune responses. Consequently, the MHC region strongly associates with numerous diseases and clinical therapies. Notoriously, the MHC region has been intractable to high-throughput analysis at complete sequence resolution, and current reference haplotypes are inadequate for large-scale studies. To address these challenges, we developed a method that specifically captures and sequences the 4.8-Mbp MHC region from genomic DNA. For 95 MHC homozygous cell lines we assembled, de novo, a set of high-fidelity contigs and a sequence scaffold, representing a mean 98% of the target region. Included are six alternative MHC reference sequences of the human genome that we completed and refined. Characterization of the sequence and structural diversity of the MHC region shows the approach accurately determines the sequences of the highly polymorphic HLA class I and HLA class II genes and the complex structural diversity of complement factor C4A/C4B It has also uncovered extensive and unexpected diversity in other MHC genes; an example is MUC22, which encodes a lung mucin and exhibits more coding sequence alleles than any HLA class I or II gene studied here. More than 60% of the coding sequence alleles analyzed were previously uncharacterized. We have created a substantial database of robust reference MHC haplotype sequences that will enable future population scale studies of this complicated and clinically important region of the human genome.

Human NK Cells Downregulate Zap70 and Syk in Response to Prolonged Activation or DNA Damage.

The extent of NK cell activity during the innate immune response affects downstream immune functions and, ultimately, the outcome of infectious or malignant disease. However, the mechanisms that terminate human NK cell responses have yet to be defined. When activation receptors expressed on NK cell surfaces bind to ligands on diseased cells, they initiate a signal that is propagated by a number of intracellular kinases, including Zap70 and Syk, eventually leading to NK cell activation. We assayed Zap70 and Syk content in NK cells from healthy human donors and identified a subset of NK cells with unusually low levels of these two kinases. We found that this Zap70lowSyklow subset consisted of NK cells expressing a range of surface markers, including CD56hi and CD56low NK cells. Upon in vitro stimulation with target cells, Zap70lowSyklow NK cells failed to produce IFN-γ and lysed target cells at one third the capacity of Zap70hiSykhi NK cells. We determined two independent in vitro conditions that induce the Zap70lowSyklow phenotype in NK cells: continuous stimulation with activation beads and DNA damage. The expression of inhibitory receptors, including NKG2A and inhibitory killer Ig-like receptors (KIRs), was negatively correlated with the Zap70lowSyklow phenotype. Moreover, expression of multiple KIRs reduced the likelihood of Zap70 downregulation during continuous activation, regardless of whether NK cells had been educated through KIR-HLA interactions in vivo. Our findings show that human NK cells are able to terminate their functional activity without the aid of other immune cells through the downregulation of activation kinases.

Different Selected Mechanisms Attenuated the Inhibitory Interaction of KIR2DL1 with C2+ HLA-C in Two Indigenous Human Populations in Southern Africa.

The functions of human NK cells in defense against pathogens and placental development during reproduction are modulated by interactions of killer cell Ig-like receptors (KIRs) with HLA-A, -B and -C class I ligands. Both receptors and ligands are highly polymorphic and exhibit extensive differences between human populations. Indigenous to southern Africa are the KhoeSan, the most ancient group of modern human populations, who have highest genomic diversity worldwide. We studied two KhoeSan populations, the Nama pastoralists and the ≠Khomani San hunter-gatherers. Comprehensive next-generation sequence analysis of HLA-A, -B, and -C and all KIR genes identified 248 different KIR and 137 HLA class I, which assort into ∼200 haplotypes for each gene family. All 74 Nama and 78 ≠Khomani San studied have different genotypes. Numerous novel KIR alleles were identified, including three arising by intergenic recombination. On average, KhoeSan individuals have seven to eight pairs of interacting KIR and HLA class I ligands, the highest diversity and divergence of polymorphic NK cell receptors and ligands observed to date. In this context of high genetic diversity, both the Nama and the ≠Khomani San have an unusually conserved, centromeric KIR haplotype that has arisen to high frequency and is different in the two KhoeSan populations. Distinguishing these haplotypes are independent mutations in KIR2DL1, which both prevent KIR2DL1 from functioning as an inhibitory receptor for C2+ HLA-C. The relatively high frequency of C2+ HLA-C in the Nama and the ≠Khomani San appears to have led to natural selection against strong inhibitory C2-specific KIR.

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.

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.

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

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.

Hematopoietic stem cell transplantation: Improving alloreactive Bw4 donor selection by genotyping codon 86 of KIR3DL1/S1.

KIR3DL1 is a natural killer (NK) cell receptor that recognizes the Bw4 epitope of human leukocyte antigen (HLA) class I molecules. Following hematopoietic stem cell transplantation for patients lacking Bw4, KIR3DL1-expressing NK cells from Bw4-positive donors can be alloreactive and eliminate tumor cells. However, KIR3DL1 alleles having T instead of C at nucleotide 320 (encoding leucine 86 instead of serine 86) are not expressed on the cell surface. Thus, not all individuals testing positive for KIR3DL1 are optimal donors for Bw4-negative recipients. Therefore, we developed a method for genotyping codon 86, which was validated by its perfect correlation with NK cell phenotype for 100 donors of diverse KIR3DL1/S1 genotype. We typed 600 donors and found that ∼12.2% had the KIR3DL1 gene, but did not express cell-surface KIR3DL1. By contrast, high-expressing allotypes were identified when haplotypes from four families with duplicated KIR3DL1/S1 genes were characterized at high resolution. Identifying donors who have KIR3DL1 but lack cell-surface KIR3DL1 would refine donor selection. With this technique, the number of individuals identified who may not be optimal donors for Bw4-negative patients increases by threefold, when compared with standard methods. Taken together, we propose that allele typing of killer cell Ig-like receptor (KIR) polymorphisms should become a standard practice when selecting donors.

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, CD33(bright) 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.

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.

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.

Very long haplotype tracts characterized at high resolution from HLA homozygous cell lines.

The HLA region of chromosome 6 contains the most polymorphic genes in humans. Spanning ~5 Mbp the densely packed region encompasses approximately 175 expressed genes including the highly polymorphic HLA class I and II loci. Most of the other genes and functional elements are also polymorphic, and many of them are directly implicated in immune function or immune-related disease. For these reasons, this complex genomic region is subject to intense scrutiny by researchers with the common goal of aiding further understanding and diagnoses of multiple immune-related diseases and syndromes. To aid assay development and characterization of the classical loci, a panel of cell lines partially or fully homozygous for HLA class I and II was assembled over time by the International Histocompatibility Working Group (IHWG). Containing a minimum of 88 unique HLA haplotypes, we show that this panel represents a significant proportion of European HLA allelic and haplotype diversity (60-95 %). Using a high-density whole genome array that includes 13,331 HLA region SNPs, we analyzed 99 IHWG cells to map the coordinates of the homozygous tracts at a fine scale. The mean homozygous tract length within chromosome 6 from these individuals is 21 Mbp. Within HLA, the mean haplotype length is 4.3 Mbp, and 65 % of the cell lines were shown to be homozygous throughout the entire region. In addition, four cell lines are homozygous throughout the complex KIR region of chromosome 19 (~250 kbp). The data we describe will provide a valuable resource for characterizing haplotypes, designing and refining imputation algorithms and developing assay controls.

Exome capture from saliva produces high quality genomic and metagenomic data.

BACKGROUND: Targeted 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. RESULTS: We 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. CONCLUSIONS: DNA 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.

KIR diversity in Maori and Polynesians: populations in which HLA-B is not a significant KIR ligand.

HLA class I molecules and killer cell immunoglobulin-like receptors (KIR) form a diverse system of ligands and receptors that individualize human immune systems in ways that improve the survival of individuals and populations. Human settlement of Oceania by island-hopping East and Southeast Asian migrants started ~3,500 years ago. Subsequently, New Zealand was reached ~750 years ago by ancestral Maori. To examine how this history impacted KIR and HLA diversity, and their functional interaction, we defined at high resolution the allelic and haplotype diversity of the 13 expressed KIR genes in 49 Maori and 34 Polynesians. Eighty KIR variants, including four 'new' alleles, were defined, as were 35 centromeric and 22 telomeric KIR region haplotypes, which combine to give >50 full-length KIR haplotypes. Two new and divergent variant KIR form part of a telomeric KIR haplotype, which appears derived from Papua New Guinea and was probably obtained by the Asian migrants en route to Polynesia. Maori and Polynesian KIR are very similar, but differ significantly from African, European, Japanese, and Amerindian KIR. Maori and Polynesians have high KIR haplotype diversity with corresponding allotype diversity being maintained throughout the KIR locus. Within the population, each individual has a unique combination of HLA class I and KIR. Characterizing Maori and Polynesians is a paucity of HLA-B allotypes recognized by KIR. Compensating for this deficiency are high frequencies (>50 %) of HLA-A allotypes recognized by KIR. These HLA-A allotypes are ones that modern humans likely acquired from archaic humans at a much earlier time.

Sarcoma microenvironment cell states and ecosystems are associated with prognosis and predict response to immunotherapy.

Characterization of the diverse malignant and stromal cell states that make up soft tissue sarcomas and their correlation with patient outcomes has proven difficult using fixed clinical specimens. Here, we employed EcoTyper, a machine-learning framework, to identify the fundamental cell states and cellular ecosystems that make up sarcomas on a large scale using bulk transcriptomes with clinical annotations. We identified and validated 23 sarcoma-specific, transcriptionally defined cell states, many of which were highly prognostic of patient outcomes across independent datasets. We discovered three conserved cellular communities or ecotypes associated with underlying genomic alterations and distinct clinical outcomes. We show that one ecotype defined by tumor-associated macrophages and epithelial-like malignant cells predicts response to immune-checkpoint inhibition but not chemotherapy and validate our findings in an independent cohort. Our results may enable identification of patients with soft tissue sarcomas who could benefit from immunotherapy and help develop new therapeutic strategies.

Comparison of NK alloreactivity prediction models based on KIR-MHC interactions in haematopoietic stem cell transplantation.

The biological processes underlying NK cell alloreactivity in haematopoietic stem cell transplantation (HSCT) remain unclear. Many different models to predict NK alloreactivity through KIR and MHC genotyping exist, raising ambiguities in its utility and application for clinicians. We assessed 27 predictive models, broadly divided into six categories of alloreactivity prediction: ligand-ligand, receptor-ligand, educational, KIR haplotype-based, KIR matching and KIR allelic polymorphism. The models were applied to 78 NGS-typed donor/recipient pairs undergoing allogeneic HSCT in genoidentical (n=43) or haploidentical (n=35) matchings. Correlations between different predictive models differed widely, suggesting that the choice of the model in predicting NK alloreactivity matters. For example, two broadly used models, educational and receptor-ligand, led to opposing predictions especially in the genoidentical cohort. Correlations also depended on the matching fashion, suggesting that this parameter should also be taken into account in the choice of the scoring strategy. The number of centromeric B-motifs was the only model strongly correlated with the incidence of acute graft-versus-host disease in our set of patients in both the genoidentical and the haploidentical cohorts, suggesting that KIR-based alloreactivity, not MHC mismatches, are responsible for it. To our best knowledge, this paper is the first to experimentally compare NK alloreactivity prediction models within a cohort of genoidentical and haploidentical donor-recipient pairs. This study helps to resolve current discrepancies in KIR-based alloreactivity predictions and highlights the need for deeper consideration of the models used in clinical studies as well as in medical practice.