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.

[1]  Mary Carrington,et al.  The impact of host genetic variation on infection with HIV-1 , 2015, Nature Immunology.

[2]  T. Wetter,et al.  Using the miraEST assembler for reliable and automated mRNA transcript assembly and SNP detection in sequenced ESTs. , 2004, Genome research.

[3]  P. Parham,et al.  Nature of polymorphism in HLA-A, -B, and -C molecules. , 1988, Proceedings of the National Academy of Sciences of the United States of America.

[4]  J. Chiaroni,et al.  RH diversity in Mali: characterization of a new haplotype RHD*DIVa/RHCE*ceTI(D2) , 2015, Transfusion.

[5]  M. Carrington,et al.  The killer immunoglobulin-like receptor gene cluster: tuning the genome for defense. , 2006, Annual review of genomics and human genetics.

[6]  P. Parham,et al.  Variable NK cell receptors and their MHC class I ligands in immunity, reproduction and human evolution , 2013, Nature Reviews Immunology.

[7]  K. Campbell,et al.  Structure/function of human killer cell immunoglobulin‐like receptors: lessons from polymorphisms, evolution, crystal structures and mutations , 2011, Immunology.

[8]  Sue Povey,et al.  Gene map of the extended human MHC , 2004, Nature Reviews Genetics.

[9]  Philip M. Kim,et al.  Paired-End Mapping Reveals Extensive Structural Variation in the Human Genome , 2007, Science.

[10]  T. Beddoe,et al.  Killer cell immunoglobulin-like receptor 3DL1-mediated recognition of human leukocyte antigen B , 2011, Nature.

[11]  S. Khakoo,et al.  NK cells: tuned by peptide? , 2015, Immunological reviews.

[12]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[13]  R. Zinkernagel,et al.  Restriction of in vitro T cell-mediated cytotoxicity in lymphocytic choriomeningitis within a syngeneic or semiallogeneic system , 1974, Nature.

[14]  Erik H. Rozemuller,et al.  13th IHWS Shared Resources Joint Report. IHWG Cell and Gene Bank and reference cell panels. , 2006 .

[15]  Jean L. Chang,et al.  Initial sequence of the chimpanzee genome and comparison with the human genome , 2005, Nature.

[16]  Peter Parham,et al.  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 , 2013, PLoS genetics.

[17]  M. Nei,et al.  Pattern of nucleotide substitution at major histocompatibility complex class I loci reveals overdominant selection , 1988, Nature.

[18]  J. Oksenberg,et al.  The killer immunoglobulin-like receptor KIR3DL1 combination with HLA-Bw4 is protective against multiple sclerosis in African Americans , 2016, Genes and Immunity.

[19]  K. Tadokoro,et al.  Large-scale comparative mapping of the MHC class I region of predominant haplotypes in Japanese , 1997, Immunogenetics.

[20]  M. Carrington,et al.  Single haplotype analysis demonstrates rapid evolution of the killer immunoglobulin-like receptor (KIR) loci in primates. , 2005, Genome research.

[21]  Sge Marsh,et al.  KIR gene characterisation of HLA homozygous cell lines. In Immunobiology of the Human MHC. , 2006 .

[22]  Ronald W. Davis,et al.  Unusual selection on the KIR3DL1/S1 natural killer cell receptor in Africans , 2007, Nature Genetics.

[23]  R. Nielsen,et al.  A Model-Based Approach for Identifying Signatures of Ancient Balancing Selection in Genetic Data , 2014, PLoS genetics.

[24]  Michael C. Schatz,et al.  Oxford Nanopore Sequencing, Hybrid Error Correction, and de novo Assembly of a Eukaryotic Genome , 2015 .

[25]  R. Kaslow,et al.  Conserved extended haplotypes of the major histocompatibility complex: further characterization , 2006, Genes and Immunity.

[26]  Gil McVean,et al.  Improved genome inference in the MHC using a population reference graph , 2014, Nature Genetics.

[27]  Jeffrey C. Boyington,et al.  Crystal structure of an NK cell immunoglobulin-like receptor in complex with its class I MHC ligand , 2000, Nature.

[28]  S. Beck,et al.  Plasticity in the organization and sequences of human KIR/ILT gene families. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[29]  H A Erlich,et al.  Analysis of enzymatically amplified beta-globin and HLA-DQ alpha DNA with allele-specific oligonucleotide probes. , 1986, Nature.

[30]  Kenny Q. Ye,et al.  Sensitive and accurate detection of copy number variants using read depth of coverage. , 2009, Genome research.

[31]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[32]  C. Vilches,et al.  Genotyping of human killer-cell immunoglobulin-like receptor genes by polymerase chain reaction with sequence-specific primers: an update. , 2002, Tissue antigens.

[33]  Heather G. Shilling,et al.  Allelic Polymorphism Synergizes with Variable Gene Content to Individualize Human KIR Genotype1 , 2002, The Journal of Immunology.

[34]  Ronald W Davis,et al.  Meiotic recombination generates rich diversity in NK cell receptor genes, alleles, and haplotypes. , 2009, Genome research.

[35]  P. Parham,et al.  Human-Specific Evolution and Adaptation Led to Major Qualitative Differences in the Variable Receptors of Human and Chimpanzee Natural Killer Cells , 2010, PLoS genetics.

[36]  Andrew R. Jones,et al.  Allele frequency net 2015 update: new features for HLA epitopes, KIR and disease and HLA adverse drug reaction associations , 2014, Nucleic Acids Res..

[37]  R. Vaughan,et al.  A multi-laboratory characterization of the KIR genotypes of 10th International Histocompatibility Workshop cell lines. , 2003, Human immunology.

[38]  C. Hurley,et al.  The characteristics of allelic polymorphism in killer-immunoglobulin-like receptor framework genes in African Americans , 2011, Immunogenetics.

[39]  Michael J. Bevan,et al.  CD8+ T Cells: Foot Soldiers of the Immune System , 2011, Immunity.

[40]  P. Parham,et al.  KIR diversity in Māori and Polynesians: populations in which HLA-B is not a significant KIR ligand , 2014, Immunogenetics.

[41]  Henry A. Erlich,et al.  Analysis of enzymatically amplified β-globin and HLA-DQα DNA with allele-specific oligonucleotide probes , 1986, Nature.

[42]  K. Shianna,et al.  Copy Number Variation of KIR Genes Influences HIV-1 Control , 2011, PLoS biology.

[43]  Christopher R. Gignoux,et al.  Exome capture from saliva produces high quality genomic and metagenomic data , 2014, BMC Genomics.

[44]  Richard Durbin,et al.  Extending reference assembly models , 2015, Genome Biology.

[45]  W. Klitz,et al.  Common and well-documented HLA alleles: 2012 update to the CWD catalogue. , 2013, Tissue antigens.

[46]  P Parham,et al.  Human diversity in killer cell inhibitory receptor genes. , 1997, Immunity.

[47]  M. A. Saper,et al.  The foreign antigen binding site and T cell recognition regions of class I histocompatibility antigens , 1987, Nature.

[48]  B. Dupont,et al.  Description of the Reference Panel of B-Lymphoblastoid Cell Lines for Factors of the HLA System: The B-Cell Line Panel Designed for the Tenth International Histocompatibility Workshop , 1989 .

[49]  Miles Parkes,et al.  Genetic insights into common pathways and complex relationships among immune-mediated diseases , 2013, Nature Reviews Genetics.

[50]  Loren Gragert,et al.  The Shaping of Modern Human Immune Systems by Multiregional Admixture with Archaic Humans , 2011, Science.

[51]  M. Colonna,et al.  Cloning of immunoglobulin-superfamily members associated with HLA-C and HLA-B recognition by human natural killer cells. , 1995, Science.

[52]  Eric O Long,et al.  Molecular clones of the p58 NK cell receptor reveal immunoglobulin-related molecules with diversity in both the extra- and intracellular domains. , 1995, Immunity.

[53]  Sara Goodwin,et al.  Oxford Nanopore sequencing, hybrid error correction, and de novo assembly of a eukaryotic genome , 2015, bioRxiv.

[54]  M. Caligiuri,et al.  Innate or Adaptive Immunity? The Example of Natural Killer Cells , 2011, Science.

[55]  Kenny Q. Ye,et al.  An integrated map of genetic variation from 1,092 human genomes , 2012, Nature.

[56]  Daniel S. Himmelstein,et al.  Association of HLA Genetic Risk Burden With Disease Phenotypes in Multiple Sclerosis. , 2016, JAMA neurology.

[57]  Amalio Telenti,et al.  Innate partnership of HLA-B and KIR3DL1 subtypes against HIV-1 , 2007, Nature Genetics.

[58]  D. Middleton,et al.  The silent KIR3DP1 gene (CD158c) is transcribed and might encode a secreted receptor in a minority of humans, in whom the KIR3DP1, KIR2DL4 and KIR3DL1/KIR3DS1 genes are duplicated , 2005, European journal of immunology.

[59]  Sge Marsh,et al.  The International Histocompatibility Workshop cell panel. , 1997 .

[60]  M. Carrington,et al.  Cutting Edge: Expansion of the KIR Locus by Unequal Crossing Over 1 , 2003, The Journal of Immunology.

[61]  James K. Bonfield,et al.  Genome analysis Advance Access publication May 30, 2010 Gap5—editing , 2010 .

[62]  M. Carrington,et al.  Inhibitory KIR3DL1 alleles are associated with psoriasis , 2016, The British journal of dermatology.

[63]  A. Tosti,et al.  Haploidentical hematopoietic transplantation from KIR ligand-mismatched donors with activating KIRs reduces nonrelapse mortality. , 2015, Blood.

[64]  P. Parham,et al.  Co-evolution of KIR2DL3 with HLA-C in a human population retaining minimal essential diversity of KIR and HLA class I ligands , 2009, Proceedings of the National Academy of Sciences.

[65]  R. Single,et al.  Killer cell immunoglobulin-like receptor (KIR) gene content variation in the HGDP-CEPH populations , 2012, Immunogenetics.

[66]  Paolo Rocco,et al.  Good laboratory practice for clinical next-generation sequencing informatics pipelines , 2015 .

[67]  S G E Marsh,et al.  Genotype List String: a grammar for describing HLA and KIR genotyping results in a text string , 2013, Tissue antigens.

[68]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[69]  Soumya Raychaudhuri,et al.  Interrogating the major histocompatibility complex with high-throughput genomics. , 2012, Human molecular genetics.

[70]  H. Muller.  The American Journal of Human Genetics Vol . 2 No . 2 June 1950 Our Load of Mutations 1 , 2006 .

[71]  J. Knight,et al.  Major histocompatibility complex genomics and human disease. , 2013, Annual review of genomics and human genetics.

[72]  Jyothi Jayaraman,et al.  Copy number variation leads to considerable diversity for B but not A haplotypes of the human KIR genes encoding NK cell receptors , 2012, Genome research.

[73]  Peter Parham,et al.  Different Patterns of Evolution in the Centromeric and Telomeric Regions of Group A and B Haplotypes of the Human Killer Cell Ig-Like Receptor Locus , 2010, PloS one.

[74]  P. Parham,et al.  Domain Shuffling Has Been the Main Mechanism Forming New Hominoid Killer Cell Ig-Like Receptors1 , 2004, The Journal of Immunology.

[75]  Peter Parham,et al.  Some human KIR haplotypes contain two KIR2DL5 genes: KIR2DL5A and KIR2DL5B , 2002, Immunogenetics.

[76]  Bi Zhou,et al.  Gene copy-number variation and associated polymorphisms of complement component C4 in human systemic lupus erythematosus (SLE): low copy number is a risk factor for and high copy number is a protective factor against SLE susceptibility in European Americans. , 2007, American journal of human genetics.

[77]  Gil McVean,et al.  Imputation of KIR Types from SNP Variation Data , 2015, American journal of human genetics.

[78]  M. Ronaghi,et al.  Very long haplotype tracts characterized at high resolution from HLA homozygous cell lines , 2015, Immunogenetics.

[79]  James Robinson,et al.  The IPD and IMGT/HLA database: allele variant databases , 2014, Nucleic Acids Res..

[80]  Mary B. Stewart,et al.  Allele-Level Haplotype Frequencies and Pairwise Linkage Disequilibrium for 14 KIR Loci in 506 European-American Individuals , 2012, PloS one.

[81]  P. Parham,et al.  Roles for HLA and KIR polymorphisms in natural killer cell repertoire selection and modulation of effector function , 2006, The Journal of experimental medicine.