IPD--the Immuno Polymorphism Database

25 years of the IPD-IMGT/HLA Database

Twenty-five years ago, in 1998, the HLA Informatics Group of the Anthony Nolan Research Institute released the IMGT/HLA Database. Since this time, this online resource has acted as the repository for the numerous variant sequences of HLA alleles named by the WHO Nomenclature Committee for Factors of the HLA System. The IPD-IMGT/HLA Database has provided a stable, highly accessible, user-friendly repository for this work. During this time, the technology underlying HLA typing has undergone significant changes. Next generation sequencing (NGS) has superseded previous methodologies of HLA typing and can generate large amounts of high-resolution sequencing data. This has resulted in a drastic increase in the number and complexity of sequences submitted to the database. The challenge for the IPD-IMGT/HLA Database has been to maintain the highest standards of curation, while supporting the core set of tools and functionality to our users with increased numbers of submissions and sequences. Traditional methods of accessing and presenting data have been challenged and new methods utilising new computing technologies have had to be developed to keep pace and support a shifting user demographic.

Four novel candidate causal variants for deficient homozygous haplotypes in Holstein cattle

Mendelian variants can determine both insemination success and neonatal survival and thus influence fertility and rearing success of cattle. We present 24 deficient homozygous haplotype regions in the Holstein population of Switzerland and provide an overview of the previously identified haplotypes in the global Holstein breed. This study encompasses massive genotyping, whole-genome sequencing (WGS) and phenotype association analyses. We performed haplotype screenings on almost 53 thousand genotyped animals including 114 k SNP data with two different approaches. We revealed significant haplotype associations to several survival, birth and fertility traits. Within haplotype regions, we mined WGS data of hundreds of bovine genomes for candidate causal variants, which were subsequently evaluated by using a custom genotyping array in several thousand breeding animals. With this approach, we confirmed the known deleterious SMC2:p.Phe1135Ser missense variant associated with Holstein haplotype (HH) 3. For two previously reported deficient homozygous haplotypes that show negative associations to female fertility traits, we propose candidate causative loss-of-function variants: the HH13-related KIR2DS1:p.Gln159* nonsense variant and the HH21-related NOTCH3:p.Cys44del deletion. In addition, we propose the RIOX1:p.Ala133_Glu142del deletion as well as the PCDH15:p.Leu867Val missense variant to explain the unexpected low number of homozygous haplotype carriers for HH25 and HH35, respectively. In conclusion, we demonstrate that with mining massive SNP data in combination with WGS data, we can map several haplotype regions and unravel novel recessive protein-changing variants segregating at frequencies of 1 to 5%. Our findings both confirm previously identified loci and expand the spectrum of undesired alleles impairing reproduction success in Holstein cattle, the world's most important dairy breed.

Laboratory Codes in Nomenclature and Scientific Communication (Advancing Organism Nomenclature in Scientific Communication to Improve Research Reporting and Reproducibility)

Laboratory registration codes, also known as laboratory codes or lab codes, are a key element in standardized laboratory animal and genetic nomenclature. As such they are critical to accurate scientific communication and to research reproducibility and integrity. The original committee on Mouse Genetic Nomenclature published nomenclature conventions for mice genetics in 1940, and then conventions for inbred strains in 1952. Unique designations were needed, and have been in use since the 1950s, for the sources of animals and substrains, for the laboratories that identified new alleles or mutations, and then for developers of transgenes and induced mutations. Current laboratory codes are typically a 2- to 4-letter acronym for an institution or an investigator. Unique codes are assigned from the International Laboratory Code Registry, which was developed and is maintained by ILAR in the National Academies (National Academies of Sciences Engineering and Medicine and previously National Academy of Sciences). As a resource for the global research community, the registry has been online since 1997. Since 2003 mouse and rat genetic and strain nomenclature rules have been reviewed and updated annually as a joint effort of the International Committee on Standardized Genetic Nomenclature for Mice and the Rat Genome and Nomenclature Committee. The current nomenclature conventions (particularly conventions for non-inbred animals) are applicable beyond rodents, although not widely adopted. Ongoing recognition, since at least the 1930s, of the research relevance of genetic backgrounds and origins of animals, and of spontaneous and induced genetic variants speaks to the need for broader application of standardized nomenclature for animals in research, particularly given the increasing numbers and complexities of genetically modified swine, nonhuman primates, fish, and other species.

High-throughput Interpretation of Killer-cell Immunoglobulin-like Receptor Short-read Sequencing Data with PING

The killer-cell immunoglobulin-like receptor (KIR) complex on chromosome 19 encodes receptors that modulate the activity of natural killer cells, and variation in these genes has been linked to infectious and autoimmune disease, as well as having bearing on pregnancy and transplant outcomes. The medical relevance and high variability of KIR genes makes short-read sequencing an attractive technology for interrogating the region, providing a high-throughput, high-fidelity sequencing method that is cost-effective. However, because this gene complex is characterized by extensive nucleotide polymorphism, structural variation including gene fusions and deletions, and a high level of homology between genes, its interrogation at high resolution has been thwarted by bioinformatic challenges, with most studies limited to examining presence or absence of specific genes. Here, we present the PING (Pushing Immunogenetics to the Next Generation) pipeline, which incorporates empirical data, novel alignment strategies and a custom alignment processing workflow to enable high-throughput KIR sequence analysis from short-read data. PING provides KIR gene copy number classification functionality for all KIR genes through use of a comprehensive alignment reference. The gene copy number determined per individual enables an innovative genotype determination workflow using genotype-matched references. Together, these methods address the challenges imposed by the structural complexity and overall homology of the KIR complex. To determine copy number and genotype determination accuracy, we applied PING to European and African validation cohorts and a synthetic dataset. PING demonstrated exceptional copy number determination performance across all datasets and robust genotype determination performance. Finally, an investigation into discordant genotypes for the synthetic dataset provides insight into misaligned reads, advancing our understanding in interpretation of short-read sequencing data in complex genomic regions. PING promises to support a new era of studies of KIR polymorphism, delivering high-resolution KIR genotypes that are highly accurate, enabling high-quality, high-throughput KIR genotyping for disease and population studies.

HLA Class I Molecules as Immune Checkpoints for NK Cell Alloreactivity and Anti-Viral Immunity in Kidney Transplantation

Natural killer (NK) cells are innate lymphocytes that can kill diseased- or virally-infected cells, mediate antibody dependent cytotoxicity and produce type I immune-associated cytokines upon activation. NK cells also contribute to the allo-immune response upon kidney transplantation either by promoting allograft rejection through lysis of cells of the transplanted organ or by promoting alloreactive T cells. In addition, they protect against viral infections upon transplantation which may be especially relevant in patients receiving high dose immune suppression. NK cell activation is tightly regulated through the integrated balance of signaling via inhibitory- and activating receptors. HLA class I molecules are critical regulators of NK cell activation through the interaction with inhibitory- as well as activating NK cell receptors, hence, HLA molecules act as critical immune checkpoints for NK cells. In the current review, we evaluate how NK cell alloreactivity and anti-viral immunity are regulated by NK cell receptors belonging to the KIR family and interacting with classical HLA class I molecules, or by NKG2A/C and LILRB1/KIR2DL4 engaging non-classical HLA-E or -G. In addition, we provide an overview of the methods to determine genetic variation in these receptors and their HLA ligands.

Estimation of German KIR Allele Group Haplotype Frequencies

The impact of the highly polymorphic Killer-cell immunoglobulin-like receptor (KIR) gene cluster on the outcome of hematopoietic stem cell transplantation (HCST) is subject of current research. To further understand the involvement of this gene family into Natural Killer (NK) cell-mediated graft-versus-leukemia reactions, knowledge of haplotype structures, and allelic linkage is of importance. In this analysis, we estimate population-specific KIR haplotype frequencies at allele group resolution in a cohort of n = 458 German families. We addressed the polymorphism of the KIR gene complex and phasing ambiguities by a combined approach. Haplotype inference within first-degree family relations allowed us to limit the number of possible diplotypes. Structural restriction to a pattern set of 92 previously described KIR copy number haplotypes further reduced ambiguities. KIR haplotype frequency estimation was finally accomplished by means of an expectation-maximization algorithm. Applying a resolution threshold of ½ n, we were able to identify a set of 551 KIR allele group haplotypes, representing 21 KIR copy number haplotypes. The haplotype frequencies allow studying linkage disequilibrium in two-locus as well as in multi-locus analyses. Our study reveals associations between KIR haplotype structures and allele group frequencies, thereby broadening our understanding of the KIR gene complex.

The IPD Project: a centralised resource for the study of polymorphism in genes of the immune system

The Immuno Polymorphism Database (IPD), https://www.ebi.ac.uk/ipd/, is a set of specialist databases that enable the study of polymorphic genes which function as part of the vertebrate immune system. The major focus is on the hyperpolymorphic major histocompatibility complex (MHC) genes and the killer-cell immunoglobulin-like receptor (KIR) genes, by providing the official repository and primary source of sequence data. Databases are centred around humans as well as animals important for food security, for companionship and as disease models. The IPD project works with specialist groups or nomenclature committees who provide and manually curate individual sections before they are submitted for online publication. To reflect the recent advance of allele sequencing technologies and the increasing demands of novel tools for the analysis of genomic variation, the IPD project is undergoing a progressive redesign and reorganisation. In this review, recent updates and future developments are discussed, with a focus on the core concepts to better future-proof the project.

Development of a novel monoclonal antibody that binds to most HLA-A allomorphs in a conformation-dependent yet peptide-promiscuous fashion

Specificity analyses of peptide binding to human leukocyte antigen (HLA)-A molecules have been hampered due to a lack of proper monoclonal antibodies (mAbs) for certain allomorphs, such as the prevalent HLA-A1 for Caucasians and HLA-A11 for Asians. We developed a mAb that recognizes a conformational epitope common to most HLA-A allomorphs. The mAb, named A-1, does not discriminate peptides by amino acid sequences, making it suitable for measuring peptide binding. A stabilization assay using TAP-deficient cell lines and A-1 was developed to investigate the specificity of peptide binding to HLA-A molecules. Regarding the evolution of HLA-A genes, the A-1 epitope has been conserved among most HLA-A allomorphs but was lost when the HLA-A gene diversified into the HLA-A*32, HLA-A*31, and HLA-A*33 lineages together with HLA-A*29 after bifurcating from the HLA-A*25 and HLA-A*26 branchs. The establishment of A-1 is expected to help researchers investigate the peptide repertoire and develop computational tools to identify cognate peptides. Since no HLA-A locus-specific mAb has been available, A-1 will also be useful for analyzing the locus-specific regulation of the HLA gene expression.

KIR2DL1 allele sequence extensions and discovery of 2DL1*0010102 and 2DL1*0010103 alleles by DNA sequencing

Full-length KIR2DL1 allele sequence extensions characterised by single molecule real-time (SMRT) DNA sequencing.

Vaccine and antibody production in plants: developments and computational tools

Plants as bioreactors have been widely used to express efficient vaccine antigens against viral, bacterial and protozoan infections. To date, many different plant-based expression systems have been analyzed, with a growing preference for transient expression systems. Antibody expression in diverse plant species for therapeutic applications is well known, and this review provides an overview of various aspects of plant-based biopharmaceutical production. Here, we highlight conventional and gene expression technologies in plants along with some illustrative examples. In addition, the portfolio of products that are being produced and how they relate to the success of this field are discussed. Stable and transient gene expression in plants, agrofiltration and virus infection vectors are also reviewed. Further, the present report draws attention to antibody epitope prediction using computational tools, one of the crucial steps of vaccine design. Finally, regulatory issues, biosafety and public perception of this technology are also discussed.

Extensive Alternative Splicing of KIR Transcripts

The killer-cell Ig-like receptors (KIR) form a multigene entity involved in modulating immune responses through interactions with MHC class I molecules. The complexity of the KIR cluster is reflected by, for instance, abundant levels of allelic polymorphism, gene copy number variation, and stochastic expression profiles. The current transcriptome study involving human and macaque families demonstrates that KIR family members are also subjected to differential levels of alternative splicing, and this seems to be gene dependent. Alternative splicing may result in the partial or complete skipping of exons, or the partial inclusion of introns, as documented at the transcription level. This post-transcriptional process can generate multiple isoforms from a single KIR gene, which diversifies the characteristics of the encoded proteins. For example, alternative splicing could modify ligand interactions, cellular localization, signaling properties, and the number of extracellular domains of the receptor. In humans, we observed abundant splicing for KIR2DL4, and to a lesser extent in the lineage III KIR genes. All experimentally documented splice events are substantiated by in silico splicing strength predictions. To a similar extent, alternative splicing is observed in rhesus macaques, a species that shares a close evolutionary relationship with humans. Splicing profiles of Mamu-KIR1D and Mamu-KIR2DL04 displayed a great diversity, whereas Mamu-KIR3DL20 (lineage V) is consistently spliced to generate a homolog of human KIR2DL5 (lineage I). The latter case represents an example of convergent evolution. Although just a single KIR splice event is shared between humans and macaques, the splicing mechanisms are similar, and the predicted consequences are comparable. In conclusion, alternative splicing adds an additional layer of complexity to the KIR gene system in primates, and results in a wide structural and functional variety of KIR receptors and its isoforms, which may play a role in health and disease.

Cost-effective and fast KIR gene-content genotyping by multiplex melting curve analysis

Killer cell immunoglobulin-like receptor (KIR) genes encode cell surface molecules that recognize HLA molecules and modulate the activity of natural killer (NK) cells. KIR genes exhibit presence and absence polymorphism, which generates a variety of gene-content haplotypes in worldwide populations. KIR gene-content variation is implicated in many diseases and is also important for placentation and transplantation. Because of the complexity of KIR polymorphism, variation in this family is still mostly studied at the gene-content level, even with the advent of next-generation sequencing (NGS) methods. Gene-content determination is generally expensive and/or time-consuming. To overcome these difficulties, we developed a method based on multiplex polymerase chain reaction with specific sequence primers (PCR-SSP) followed by melting curve analysis that allows cost-effective, precise and fast generation of results. Our method was 100% concordant with a gel-based method and 99.9% concordant with presence and absence determination by NGS. The limit of detection for accurate typing was 30 ng of DNA (0.42 μM) with 260/230 and 260/280 ratios as low as 0.19 and of 0.44. In addition, we developed a user-friendly Java-based computational application called killerPeak that interprets the raw data generated by Viia7 or QuantStudio 7 quantitative PCR machines and reliably exports the final genotyping results in spreadsheet file format. The combination of a reliable method that requires low amount of DNA with an automated interpretation of results allows scaling the KIR genotyping in large cohorts with reduced turnaround time.

Comparative MHC nomenclature: report from the ISAG/IUIS-VIC committee 2018

Significant progress has been made over the last decade in defining major histocompatibility complex (MHC) diversity at the nucleotide, allele, haplotype, diplotype, and population levels in many non-human species. Much of this progress has been driven by the increased availability and reduced costs associated with nucleotide sequencing technologies. This report provides an update on the activities of the comparative MHC nomenclature committee which is a standing committee of both the International Society for Animal Genetics (ISAG) and the International Union of Immunological Societies (IUIS) where it operates under the umbrella of the Veterinary Immunology Committee (VIC). A previous report from this committee in 2006 defined the role of the committee in providing guidance in the development of a standardized nomenclature for genes and alleles at MHC loci in non-human species. It described the establishment of the Immuno Polymorphism Database, IPD-MHC, which continues to provide public access to high quality MHC sequence data across a range of species. In this report, guidelines for the continued development of a universal MHC nomenclature framework are described, summarizing the continued development of each species section within the IPD-MHC project.

Bridging immunogenetics and immunoproteomics: Model positional scanning library analysis for Major Histocompatibility Complex class II DQ in Tursiops truncatus

The Major Histocompatibility Complex (MHC) is a critical element in mounting an effective immune response in vertebrates against invading pathogens. Studies of MHC in wildlife populations have typically focused on assessing diversity within the peptide binding regions (PBR) of the MHC class II (MHC II) family, especially the DQ receptor genes. Such metrics of diversity, however, are of limited use to health risk assessment since functional analyses (where changes in the PBR are correlated to recognition/pathologies of known pathogen proteins), are difficult to conduct in wildlife species. Here we describe a means to predict the binding preferences of MHC proteins: We have developed a model positional scanning library analysis (MPSLA) by harnessing the power of mixture based combinatorial libraries to probe the peptide landscapes of distinct MHC II DQ proteins. The algorithm provided by NNAlign was employed to predict the binding affinities of sets of peptides generated for DQ proteins. These binding affinities were then used to retroactively construct a model Positional Scanning Library screen. To test the utility of the approach, a model screen was compared to physical combinatorial screens for human MHC II DP. Model library screens were generated for DQ proteins derived from sequence data from bottlenose dolphins from the Indian River Lagoon (IRL) and the Atlantic coast of Florida, and compared to screens of DQ proteins from Genbank for dolphin and three other cetaceans. To explore the peptide binding landscape for DQ proteins from the IRL, combinations of the amino acids identified as active were compiled into peptide sequence lists that were used to mine databases for representation in known proteins. The frequency of which peptide sequences predicted to bind the MHC protein are found in proteins from pathogens associated with marine mammals was found to be significant (p values <0.0001). Through this analysis, genetic variation in MHC (classes I and II) can now be associated with the binding repertoires of the expressed MHC proteins and subsequently used to identify target pathogens. This approach may be eventually applied to evaluate individual population and species risk for outbreaks of emerging diseases.

IPD-MHC: nomenclature requirements for the non-human major histocompatibility complex in the next-generation sequencing era

The IPD-MHC Database is the official repository for non-human MHC sequences, overseen and supported by the Comparative MHC Nomenclature Committee, providing access to curated MHC data and associated analysis tools. To address the increasing amount and complexity of data being submitted, an entirely upgraded version of the IPD-MHC Database was recently released to maintain IPD-MHC as the central platform for the comparison of curated MHC data. As a consequence, a new level of nomenclature standardisation is required between the different species to enable data submission and to allow the unambiguous inter- and intra-species comparison of alleles. However, any changes must retain the flexibility demanded by the unique biology of different taxonomic groups. Here, we describe the rationale for a standardised nomenclature system and summarise the changes that have been driven by the requirements of implementing the IPD-MHC database. This modified nomenclature system is essential to maintain the current functionality of IPD-MHC and provide a scalable future-proof database organisation to fully exploit the bioinformatic tools used for analysis.

KIR3DL1 alleles and their epistatic interactions with human leukocyte antigen class I influence resistance and susceptibility to HIV-1 acquisition in the Pumwani sex worker cohort

OBJECTIVE

To determine the associations of KIR3DL1/S1(3DL1/S1) and its epistatic interactions with human leukocyte antigen class I (HLA-I) alleles with resistance and susceptibility to HIV-1.

DESIGN

Despite repeated exposure to HIV-1, a subset of women enrolled in the Pumwani sex worker cohort remain HIV uninfected. Previous studies have shown that specific HLA class I and II alleles were associated with this natural immunity. In this study, we investigated the association of 3DL1/S1 and its epistatic interactions with HLA-I, with resistance or susceptibility to HIV-1 acquisition.

METHODS

We used a sequence-based typing method to genotype 3DL1/S1 of 641 women in this cohort. The association of 3DL1/S1 and its epistatic interactions with HLA-I were analyzed using SPSS statistics software.

RESULTS

3DL1041 is enriched in the HIV-1-resistant women [P = 0.009, Pc = 0.0468, odds ratio (OR): 3.359, 95% confidence interval (CI): 1.39-8.32], whereas, 3DL1020 was associated with susceptibility to HIV-1 infection before correction for multiple comparisons (P = 0.029, Pc = 0.0858, OR: 0.316, 95%CI: 0.10-1.04). Epistatic interactions between several 3DL1 alleles and specific HLA-I alleles were observed. Among them the cocarriage of 3DL1041 with Bw4 (P = 1E - 05, Pc = 0.0015, OR: 13.33, 95%CI: 3.43-51.9), or Bw6 (P = 0.008, Pc = 0.272, OR: 3.92, 95%CI: 1.51-10.17), increased the odds of remaining HIV-1 uninfected. Further, 3DL1041+/Bw4+ women who entered the cohort HIV negative remained uninfected (P = 0.032, Pc = 0.0858). Cocarriage of 3DL101501 with C02 : 10 (P = 2.73E - 07, Pc = 7.0954E - 06), B15 : 03 (P = 3.21E - 04, Pc = 0.0042), A24 supertype (P = 8.89E - 04, Pc = 0.0077), or A23 : 01 (P = 0.0036, Pc = 0.0236) was associated with increased susceptibility to seroconversion.

CONCLUSION

The effects of interactions between 3DL1 and HLA-I alleles on resistance/susceptibility to HIV-1 infection suggest that innate immunity plays an important role in HIV-1 acquisition and should be studied and explored for HIV prevention.

Weaker HLA Footprints on HIV in the Unique and Highly Genetically Admixed Host Population of Mexico

HIV circumvents HLA class I-restricted CD8⁺ T-cell responses through selection of escape mutations that leave characteristic mutational “footprints,” also known as HLA-associated polymorphisms (HAPs), on HIV sequences at the population level. While many HLA footprints are universal across HIV subtypes and human populations, others can be region specific as a result of the unique immunogenetic background of each host population. Using a published probabilistic phylogenetically informed model, we compared HAPs in HIV Gag and Pol (PR-RT) in 1,612 subtype B-infected, antiretroviral treatment-naive individuals from Mexico and 1,641 individuals from Canada/United States. A total of 252 HLA class I allele subtypes were represented, including 140 observed in both cohorts, 67 unique to Mexico, and 45 unique to Canada/United States. At the predefined statistical threshold of a q value of <0.2, 358 HAPs (201 in Gag, 157 in PR-RT) were identified in Mexico, while 905 (534 in Gag and 371 in PR-RT) were identified in Canada/United States. HAPs identified in Mexico included both canonical HLA-associated escape pathways and novel associations, in particular with HLA alleles enriched in Amerindian and mestizo populations. Remarkably, HLA footprints on HIV in Mexico were not only fewer but also, on average, significantly weaker than those in Canada/United States, although some exceptions were noted. Moreover, exploratory analyses suggested that the weaker HLA footprint on HIV in Mexico may be due, at least in part, to weaker and/or less reproducible HLA-mediated immune pressures on HIV in this population. The implications of these differences for natural and vaccine-induced anti-HIV immunity merit further investigation.

IMPORTANCE HLA footprints on HIV identify viral regions under intense and consistent pressure by HLA-restricted immune responses and the common mutational pathways that HIV uses to evade them. In particular, HLA footprints can identify novel immunogenic regions and/or epitopes targeted by understudied HLA alleles; moreover, comparative analyses across immunogenetically distinct populations can illuminate the extent to which HIV immunogenic regions and escape pathways are shared versus population-specific pathways, information which can in turn inform the design of universal or geographically tailored HIV vaccines. We compared HLA-associated footprints on HIV in two immunogenetically distinct North American populations, those of Mexico and Canada/United States. We identify both shared and population-specific pathways of HIV adaptation but also make the surprising observation that HLA footprints on HIV in Mexico overall are fewer and weaker than those in Canada/United States, raising the possibility that HLA-restricted antiviral immune responses in Mexico are weaker, and/or escape pathways somewhat less consistent, than those in other populations.

Improved full-length killer cell immunoglobulin-like receptor transcript discovery in Mauritian cynomolgus macaques

Killer cell immunoglobulin-like receptors (KIRs) modulate disease progression of pathogens including HIV, malaria, and hepatitis C. Cynomolgus and rhesus macaques are widely used as nonhuman primate models to study human pathogens, and so, considerable effort has been put into characterizing their KIR genetics. However, previous studies have relied on cDNA cloning and Sanger sequencing that lack the throughput of current sequencing platforms. In this study, we present a high throughput, full-length allele discovery method utilizing Pacific Biosciences circular consensus sequencing (CCS). We also describe a new approach to Macaque Exome Sequencing (MES) and the development of the Rhexome1.0, an adapted target capture reagent that includes macaque-specific capture probe sets. By using sequence reads generated by whole genome sequencing (WGS) and MES to inform primer design, we were able to increase the sensitivity of KIR allele discovery. We demonstrate this increased sensitivity by defining nine novel alleles within a cohort of Mauritian cynomolgus macaques (MCM), a geographically isolated population with restricted KIR genetics that was thought to be completely characterized. Finally, we describe an approach to genotyping KIRs directly from sequence reads generated using WGS/MES reads. The findings presented here expand our understanding of KIR genetics in MCM by associating new genes with all eight KIR haplotypes and demonstrating the existence of at least one KIR3DS gene associated with every haplotype.

Killer Immunoglobulin-Like Receptor Profiles Are not Associated with Risk of Amoxicillin-Clavulanate-Induced Liver Injury in Spanish Patients

Natural killer cells are an integral part of the immune system and represent a large proportion of the lymphocyte population in the liver. The activity of these cells is regulated by various cell surface receptors, such as killer Ig-like receptors (KIR) that bind to human leukocyte antigen (HLA) class I ligands on the target cell. The composition of KIR receptors has been suggested to influence the development of specific diseases, in particularly autoimmune diseases, cancer and reproductive diseases. The role played in idiosyncratic drug-induced liver injury (DILI) is currently unknown. In this study, we examined KIR gene profiles and HLA class I polymorphisms in amoxicillin-clavulanate (AC) DILI patients in search for potential risk associations. One hundred and two AC DILI patients and 226 controls were genotyped for the presence or absence of 16 KIR loci, including the two pseudogenes 2DP1 and 3DP1. No significant differences were found in the distribution of individual KIRs between patients and controls, which were comparable to previously reported KIR data from ethnically similar cohorts. The 21.6 and 21.2% of the patients and controls, respectively, were homozygous haplotype A carriers, while 78.4 and 78.8%, respectively, contained at least one B haplotype (Bx). The genotypes translated into 27 (AC DILI) and 46 (controls) different gene profiles, with 19 being present in both groups. The most frequent Bx gene profile containing KIRs 2DS2, 2DL2, 2DL3, 2DP1, 2DL1, 3DL1, 2DS4, 3DL2, 3DL3, 2DL4, and 3PD1 was present in 16% of the DILI patients and 14% of the controls. The distribution of HLA class I epitopes did not differ significantly between AC DILI patients and controls. The most frequent receptor-ligand combinations in the DILI patients were 2DL3 + epitope C1 (67%) and 3DL1 + Bw4 motif (67%), while 2DL1 + epitope C2 (69%) and 3DL1 + Bw4 motif (69%) predominated in the controls. This is to our knowledge the first analysis of KIR receptor-HLA ligand associations in DILI, although our findings do not support evidence of these genetic variations playing a major role in AC DILI development.

Expression of bovine non-classical major histocompatibility complex class I proteins in mouse P815 and human K562 cells

Major histocompatibility complex class I (MHC-I) proteins can be expressed as cell surface or secreted proteins. To investigate whether bovine non-classical MHC-I proteins are expressed as cell surface or secreted proteins, and to assess the reactivity pattern of monoclonal antibodies with non-classical MHC-I isoforms, we expressed the MHC proteins in murine P815 and human K562 (MHC-I deficient) cells. Following antibiotic selection, stably transfected cell lines were stained with H1A or W6/32 antibodies to detect expression of the MHC-I proteins by flow cytometry. Two non-classical proteins (BoLA-NC1*00501 and BoLA-NC3*00101) were expressed on the cell surface in both cell lines. Surprisingly, the BoLA-NC4*00201 protein was expressed on the cell membrane of human K562 but not mouse P815 cells. Two non-classical proteins (BoLA-NC1*00401, which lacks a transmembrane domain, and BoLA-NC2*00102) did not exhibit cell surface expression. Nevertheless, Western blot analyses demonstrated expression of the MHC-I heavy chain in all transfected cell lines. Ammonium-sulfate precipitation of proteins from culture supernatants showed that BoLA-NC1*00401 was secreted and that all surface expressed proteins where shed from the cell membrane by the transfected cells. Interestingly, the surface expressed MHC-I proteins were present in culture supernatants at a much higher concentration than BoLA-NC1*00401. This comprehensive study shows that bovine non-classical MHC-I proteins BoLA-NC1*00501, BoLA-NC3*00101, and BoLA-NC4*00201 are expressed as surface isoforms with the latter reaching the cell membrane only in K562 cells. Furthermore, it demonstrated that BoLA-NC1*00401 is a secreted isoform and that significant quantities of membrane associated MHC-I proteins can be shed from the cell membrane.

Association of KIR3DL1/S1 and HLA-Bw4 with CD4 T cell counts in HIV-infected Mexican mestizos

Certain genotypic combinations of killer-cell immunoglobulin-like receptors (KIR) and human leukocyte antigens (HLA) have been associated with favourable outcomes after exposure to human immunodeficiency virus in Caucasoid and African populations. Human immunodeficiency virus (HIV) infection is characterized by a rapid exhaustion of CD4 cells, which results in impaired cellular immunity. During this early phase of infection, it is thought that the natural killer (NK) cells represent the main effector arm of the host immune response to HIV. This study investigates whether KIR and HLA factors are associated to CD4 T cell numbers after HIV infection in Mexican mestizos as assessed at the time of initial medical evaluation and subsequent clinical follow-up. KIR and HLA-B gene carrier frequency differences were compared between groups of patients stratified by CD4 T cell numbers as assessed during their first medical evaluation (a point in time at which all patients were anti-retroviral therapy naïve). In addition, the influence that these genetic factors have on averaged historical CD4 cell counts in patients subjected to follow-up (mostly therapy-experienced) was also evaluated. Our results suggest a protective role for the HLA-Bw4 and KIR3D + Bw4 combination in both therapy-naïve and therapy-experienced patients. This report furthers our understanding on the way that immune genes modulate HIV disease progression in less-studied human populations such as the Mexican mestizos with a special focus on CD4 T cell number and behaviour.

2DL1, 2DL2 and 2DL3 all contribute to KIR phenotype variability on human NK cells

Natural killer (NK) cells are lymphocytes that function as part of the innate immune system. Their activity is controlled by a range of inhibitory and activating receptors, including the important killer-cell immunoglobulin-like receptors (KIR). The KIR are a multi-gene family of receptors that interact with the human leukocyte antigen (HLA) class I family of molecules and are characterised by extensive allelic polymorphism. Their expression on the cell surface of NK cells is highly variable, but the factors responsible for this variability are not yet clearly understood. In the current study, we investigated KIR expression in a healthy human cohort that we had previously characterised in depth at a genetic level, with KIR allele typing and HLA class I ligand genotypes available for all donors (n=198). Allelic polymorphism significantly affected the phenotypic expression of all KIR analysed, whereas HLA ligand background influenced the expression levels of 2DL1 and 2DL3. In particular, we found that although 2DL2 may influence 2DL1 expression, this appears to be owing to variation in 2DL1 copy number. Finally, the inhibitory receptor LILRB1 had higher expression levels in individuals with B/B KIR genotypes, suggesting a possible relationship between KIR and non-KIR receptors, which serves to balance NK cell activation potential.

A human-specific allelic group of the MHC DRB1 gene in primates

Background

Diversity among human leukocyte antigen (HLA) molecules has been maintained by host-pathogen coevolution over a long period of time. Reflecting this diversity, the HLA loci are the most polymorphic in the human genome. One characteristic of HLA diversity is long-term persistence of allelic lineages, which causes trans-species polymorphisms to be shared among closely related species. Modern humans have disseminated across the world after their exodus from Africa, while chimpanzees have remained in Africa since the speciation event between humans and chimpanzees. It is thought that modern humans have recently acquired resistance to novel pathogens outside Africa. In the present study, we investigated HLA alleles that could contribute to this local adaptation in humans and also studied the contribution of natural selection to human evolution by using molecular data.

Results

Phylogenetic analysis of HLA-DRB1 genes identified two major groups, HLA Groups A and B. Group A formed a monophyletic clade distinct from DRB1 alleles in other Catarrhini, suggesting that Group A is a human-specific allelic group. Our estimates of divergence time suggested that seven HLA-DRB1 Group A allelic lineages in humans have been maintained since before the speciation event between humans and chimpanzees, while chimpanzees possess only one DRB1 allelic lineage (Patr-DRB1*03), which is a sister group to Group A. Experimental data showed that some Group A alleles bound to peptides derived from human-specific pathogens. Of the Group A alleles, three exist at high frequencies in several local populations outside Africa.

Conclusions

HLA Group A alleles are likely to have been retained in human lineages for a long period of time and have not expanded since the divergence of humans and chimpanzees. On the other hand, most orthologs of HLA Group A alleles may have been lost in the chimpanzee due to differences in selective pressures. The presence of alleles with high frequency outside of Africa suggests these HLA molecules result from the local adaptations of humans. Our study helps elucidate the mechanism by which the human adaptive immune system has coevolved with pathogens over a long period of time.

Nonsynonymous substitution rate heterogeneity in the peptide-binding region among different HLA-DRB1 lineages in humans

An extraordinary diversity of amino acid sequences in the peptide-binding region (PBR) of human leukocyte antigen [HLA; human major histocompatibility complex (MHC)] molecules has been maintained by balancing selection. The process of accumulation of amino acid diversity in the PBR for six HLA genes (HLA-A, B, C, DRB1, DQB1, and DPB1) shows that the number of amino acid substitutions in the PBR among alleles does not linearly correlate with the divergence time of alleles at the six HLA loci. At these loci, some pairs of alleles show significantly less nonsynonymous substitutions at the PBR than expected from the divergence time. The same phenomenon was observed not only in the HLA but also in the rat MHC. To identify the cause for this, DRB1 sequences, a representative case of a typical nonlinear pattern of substitutions, were examined. When the amino acid substitutions in the PBR were placed with maximum parsimony on a maximum likelihood tree based on the non-PBR substitutions, heterogeneous rates of nonsynonymous substitutions in the PBR were observed on several branches. A computer simulation supported the hypothesis that allelic pairs with low PBR substitution rates were responsible for the stagnation of accumulation of PBR nonsynonymous substitutions. From these observations, we conclude that the nonsynonymous substitution rate at the PBR sites is not constant among the allelic lineages. The deceleration of the rate may be caused by the coexistence of certain pathogens for a substantially long time during HLA evolution.

The yin-yang of KIR3DL1/S1: molecular mechanisms and cellular function

Killer Immunoglobulin-like Receptors (KIR) are a family of receptors expressed on natural killer (NK) and T-cell subsets. KIR3DL1 is a highly polymorphic receptor that binds to groups of HLAA and HLA-B allotypes that express the Bw4 epitope. The variation in KIR3DL1 allotypes manifests at a number of levels. Most dramatically, a common allelic variant encodes an activating rather than an inhibitory receptor (KIR3DS1). In addition, sequence variants can affect both the frequency of expression within the NK cell population and the intensity of expression on a given cell. KIR3DL1 polymorphism also influences the interaction with HLA-Bw4 molecules, due to contacts with the HLA molecule itself and sensitivity to the presented peptide. A body of evidence from genetic association studies supports the biological significance not only of the interaction of KIR3DL1 with HLA-Bw4 but also the functional variation seen with different KIR3DL1 and HLA allotypes. In this review, we discuss our current understanding of KIR3DL1 function and our recent insights from the structure of the KIR3DL1 in complex with HLA. In addition, we will summarize our current understanding of KIR3DS1, including its ligand specificity and its role in immune responses.

External quality assessment in molecular immunohematology: the INSTAND proficiency test program

BACKGROUND

Genotyping for red blood cell (RBC), platelet (PLT), and granulocyte antigens is a new tool for clinical pathology, transfusion medicine services, and blood banks. Proficiency in laboratory tests can be established by external quality assessments (EQAs), which are required for clinical application in many health care systems. There are few EQAs for molecular immunohematology.

STUDY DESIGN AND METHODS

We analyzed the participation and pass rates in an EQA for RBC, PLT, and granulocyte antigens. This EQA was distributed by INSTAND, a large nonprofit provider of proficiency tests, twice per year since Fall 2006 as EQA Number 235 Immunohematology A (molecular diagnostic). The coordinators defined at the outset which alleles are mandatory for detection.

RESULTS

The number of participants steadily increased from 51 to 73 per proficiency by Fall 2012. More than 60 institutions utilized this EQA at least once a year. Approximately 80% of them participated in RBC, 68% in PLT, and 22% in granulocyte systems. With the exceptions of RHD (82%) and granulocytes (85%), pass rates exceeded 93%. While the pass rate increased for granulocyte and decreased for the ABO system, the pass rates for the other systems changed little over 6½ years.

CONCLUSIONS

The INSTAND proficiency test program was regularly used for EQA by many institutions, particularly in Central Europe. While the technical standards and pass rates in the participating laboratories were high, there has been little improvement in pass rates since 2006.

Polymorphisms and tissue expression of the feline leukocyte antigen class I loci FLAI-E, FLAI-H, and FLAI-K

Cytotoxic CD8+ T-cell immunosurveillance for intracellular pathogens, such as viruses, is controlled by classical major histocompatibility complex (MHC) class Ia molecules, and ideally, these antiviral T-cell populations are defined by the specific peptide and restricting MHC allele. Surprisingly, despite the utility of the cat in modeling human viral immunity, little is known about the feline leukocyte antigen class I complex (FLAI). Only a few coding sequences with uncertain locus origin and expression patterns have been reported. Of 19 class I genes, three loci--FLAI-E, FLAI-H, and FLAI-K--are predicted to encode classical molecules, and our objective was to evaluate their status by analyzing polymorphisms and tissue expression. Using locus-specific, PCR-based genotyping, we amplified 33 FLAI-E, FLAI-H, and FLAI-K alleles from 12 cats of various breeds, identifying, for the first time, alleles across three distinct loci in a feline species. Alleles shared the expected polymorphic and invariant sites in the α1/α2 domains, and full-length cDNA clones possessed all characteristic class Ia exons. Alleles could be assigned to a specific locus with reasonable confidence, although there was evidence of potentially confounding interlocus recombination between FLAI-E and FLAI-K. Only FLAI-E, FLAI-H, and FLAI-K origin alleles were amplified from cDNAs of multiple tissue types. We also defined hypervariable regions across these genes, which permitted the assignment of names to both novel and established alleles. As predicted, FLAI-E, FLAI-H, and FLAI-K fulfill the major criteria of class Ia genes. These data represent a necessary prerequisite for studying epitope-specific antiviral CD8+ T-cell responses in cats.

Comprehensive identification of high-frequency and co-occurring Mafa-DPA1, Mafa-DQA1, Mafa-DRA, and Mafa-DOA alleles in Vietnamese cynomolgus macaques

High-frequency alleles and/or co-occurring human leukocyte antigen alleles across loci appear to be more important than individual alleles as markers of disease risk and have clinical value as biomarkers for targeted screening or the development of new disease therapies. To better elucidate the major histocompatibility complex (MHC) background and to facilitate the experimental use of cynomolgus macaques, Mafa-DPA1, Mafa-DQA1, Mafa-DRA, and Mafa-DOA alleles were characterized, and their combinations were investigated in 30 Vietnamese macaques by gene cloning and sequencing. A total of 26 Mafa-DPA1, 18 Mafa-DQA1, 9 Mafa-DRA, and 15 Mafa-DOA alleles, including 7 high-frequency alleles, were identified in this study, respectively. In addition, 15 Mafa-DQA1, 17 Mafa-DPA1, 15 Mafa-DOA, and 2 Mafa-DRA alleles represented novel sequences that had not been documented in earlier studies. Our results also showed that the Vietnamese macaques might be valuable because no less than 30% of the test animals possessed Mafa-DRA*01:02:01 (90%), -DQA1*26:01:03 (37%), -DOA*01:02:07 (34%), and -DQA1*01:03:03 (30%). We previously reported that the combinations of MHC class II alleles, including the combination of DOA*01:02:07-DPA1*02:09 and DOA*01:02:07-DQA1*01:03:03, were detected in 17 and 14% of the animals, respectively. Interestingly, more than two Mafa-DQA1 and Mafa-DPA1 alleles were detected in one animal in this study, which suggested that they might be caused by a chromosomal duplication. If our findings can be validated by other studies, it will further enrich the number of known Mafa-DPA1 and Mafa-DQA1 polymorphisms. Our results identified the co-occurring MHC alleles across loci in a cohort of Vietnamese cynomolgus macaques, which emphasized the value of this species as a model for biomedical research.

Evaluating genome architecture of a complex region via generalized bipartite matching

With the remarkable development in inexpensive sequencing technologies and supporting computational tools, we have the promise of medicine being personalized by knowledge of the individual genome. Current technologies provide high throughput, but short reads. Reconstruction of the donor genome is based either on de novo assembly of the (short) reads, or on mapping donor reads to a standard reference. While such techniques demonstrate high success rates for inferring 'simple' genomic segments, they are confounded by segments with complex duplication patterns, including regions of direct medical relevance, like the HLA and the KIR regions.In this work, we address this problem with a method for assessing the quality of a predicted genome sequence for complex regions of the genome. This method combines two natural types of evidence: sequence similarity of the mapped reads to the predicted donor genome, and distribution of reads across the predicted genome. We define a new scoring function for read-to-genome matchings, which penalizes for sequence dissimilarities and deviations from expected read location distribution, and present an efficient algorithm for finding matchings that minimize the penalty. The algorithm is based on a formal problem, first defined in this paper, called Coverage Sensitive many-to-many min-cost bipartite Matching (CSM). This new problem variant generalizes the standard (one-to-one) weighted bipartite matching problem, and can be solved using network flows. The resulting Java-based tool, called SAGE (Scoring function for Assembled GEnomes), is freely available upon request. We demonstrate over simulated data that SAGE can be used to infer correct haplotypes of the highly repetitive KIR region on the Human chromosome 19.

Using databases and data mining in vaccinology

Throughout time functional immunology has accumulated vast amounts of quantitative and qualitative data relevant to the design and discovery of vaccines. Such data includes, but is not limited to, components of the host and pathogen genome (including antigens and virulence factors), T- and B-cell epitopes and other components of the antigen presentation pathway and allergens. In this review the authors discuss a range of databases that archive such data. Built on such information, increasingly sophisticated data mining techniques have developed that create predictive models of utilitarian value. With special reference to epitope data, the authors discuss the strengths and weaknesses of the available techniques and how they can aid computer-aided vaccine design deliver added value for vaccinology.

Identification of MhcMafa-DRB alleles in a cohort of cynomolgus macaques of Vietnamese origin

Cynomolgus macaques have been used widely to build a research model of infectious and chronic diseases, as well as in transplantation studies, where disease susceptibility and/or resistance are associated with the major histocompatibility complex (MHC). To better elucidate polymorphisms and genetic differences in the Mafa-DRB locus, and facilitate the experimental use of cynomolgus macaques, we used pool screening combined with cloning and direct sequencing of polymerase chain reaction products to characterize MhcMafa-DRB gene alleles in 153 Vietnamese cynomolgus macaques. We identified 30 Mafa-DRB alleles belonging to 17 allelic lineages, including four novel sequences that had not been documented in earlier reports. The highest frequency allele was Mafa-DRB*W27:04, which was present in 7 of 35 (20%) monkeys. The next most frequent alleles were Mafa-DRB*3:07 and Mafa-DRB*W7:01, which were detected in 5 of 35 (14.3%) and 4 of 35 (11.4%) of the monkeys, respectively. The high-frequency alleles in this Vietnamese population may be high priority targets for additional characterization of immune functions. Only the DRB1*03 and DRB1*10 lineages were also present in humans, whereas the remaining alleles were monkey-specific lineages. We found 25 variable sites by aligning the deduced amino acid sequences of 29 identified alleles. Evolutionary and population analyses based on these sequences showed that human, rhesus, and cynomolgus macaques share several Mhc-DRB lineages and the shared polymorphisms in the DRB region may be attributable to the existence of interbreeding between rhesus and cynomolgus macaques. This information will promote the understanding of MHC diversity and polymorphism in cynomolgus macaques and increase the value of this species as a model for biomedical research.

Comprehensive identification of high-frequency and co-occurring Mafa-B, Mafa-DQB1, and Mafa-DRB alleles in cynomolgus macaques of Vietnamese origin

High-frequency alleles and/or co-occurring human leukocyte antigen (HLA) alleles across loci appear to be more important than individual alleles, because they might be markers of disease risk that have clinical value as biomarkers for targeted screening or the development of new therapies. To better elucidate the major histocompatibility complex background and to facilitate the experimental use of cynomolgus macaques, Mafa-B, Mafa-DQB1, and Mafa-DRB alleles were characterized and their combinations were investigated from 30 macaques of Vietnamese origin by cloning and sequencing. A total of 48 Mafa-B, 22 Mafa-DQB1, and 42 Mafa-DRB alleles, were detected in this study, respectively. In addition, two Mafa-DQB1 and eight Mafa-DRB alleles represented novel sequences that had not been documented in earlier studies. Our results also showed that the macaque from Vietnam might be valuable because >30% of the test animals possessed Mafa-DRB*w304 (30%) and -DQB1*0616 (30%). We report that the combination of major histocompatibility complex (MHC) class I and II alleles, including the combination of DRB3*0403-DRB*w304, DRB1*1013-DRB*w304, and Mafa-B*007:01:01-DRB*w304, which was in 17%, 13%, and 13% of the animals, respectively. Interesting, more than two Mafa-DQB1 alleles detected in one animal in this study suggest that Mafa-DQB1, like Mafa-DRB, might be a duplication in the chromosome, which have ever been documented in cynomolgus monkeys but has not yet been observed in rhesus macaques or other primates. Our results for the high frequency of commonly co-occurring MHC alleles across loci in a cohort of the Vietnamese cynomolgus macaque emphasized the value of this species as a model for biomedical research.

Study of cynomolgus monkey (Macaca fascicularis) Mhc DRB gene polymorphism in four populations

The cynomolgus macaque (Macaca fascicularis) is currently used as an animal model in various fields of immunology especially in the development of innovative vaccines for the prevention and treatment of infectious diseases. The polymorphism of the major histocompatibility complex (MHC) influences the development of adaptive immune responses, and it is crucial to characterize the polymorphism of cynomolgus MHC genes. Among all macaque species, the cynomolgus macaque has the most diversified geographical area encompassing continental and insular populations. By the study of a large sample of animals from the Philippines (N = 359), we have characterized 20 DRB haplotypes. The DRB genotyping was performed by denaturing gradient gel electrophoresis (DGGE) sequencing of exon 2 and was confirmed by polymerase chain reaction-sequence-specific oligonucleotide. The DRB and DRA cDNA of 126 animals were characterized by cloning and sequencing. By means of DGGE sequencing, we characterized the polymorphism of genomic DRB exon 2 in three other cynomolgus macaque population samples (Java, Vietnam, and Mauritius), and we discuss about the origin of the founders of the Mauritian and the Filipino cynomolgus macaque populations.

Pregnancy close to the edge: an immunosuppressive infiltrate in the chorionic plate of placentas from uncomplicated egg cell donation

In pregnancies achieved after egg donation (ED) tolerance towards a completely allogeneic fetus is mediated by several complex immunoregulatory mechanisms, of which numerous aspects are still unknown. A distinct lesion not described previously in the literature, was repeatedly found in the chorionic plate in a substantial portion of placentas from ED pregnancies, but never in placentas from normal term pregnancies. The aim of this study was to assess its origin and its cellular composition. The relation between the lesion, the clinical and histological parameters were assessed. In addition we investigated the relation with the number of HLA-mismatches and KIR genotype of mother and child.In ten out of twenty-six (38.5%) placentas from ED pregnancies an inflammatory lesion was present in the chorionic plate. A significantly lower incidence of pre-eclampsia was found in the group with the lesion; 0% versus 45.5%. A significant relation was found between this lesion and the presence of intervillositis, chronic deciduitis, presence of plasma cells and fibrin deposition in the decidua. Fluorescent in situ hybridisation with X/Y-chromosome probes showed that the majority of cells present in the lesion are of maternal origin. The expression of the macrophage marker CD14+ and of the type 2 macrophage (M2) marker CD163+ was significantly higher in the lesion. The incidence of a fetal HLA-C2 genotype was significantly higher in cases with a lesion compared to the group without the lesion. In conclusion, a striking relationship was observed between the presence of a not previously described inflammatory lesion in the chorionic plate and the absence of pre-eclampsia in ED pregnancies. The lesion consists of mainly maternal cells with a higher expression of the macrophage marker CD14+ and the M2 marker CD163+. These findings suggest a protective immune mechanism which might contribute to the prevention of severe clinical complications like pre-eclampsia.

Killer-cell immunoglobulin-like receptor genotyping and HLA killer-cell immunoglobulin-like receptor-ligand identification by real-time polymerase chain reaction

The effector function of natural killer (NK) cells is modulated by surface expression of a range of killer-cell immunoglobulin-like receptors (KIRs) that interact with human leukocyte antigen (HLA) class I ligands. We describe the use of real-time polymerase chain reaction (PCR) assays that allow easy and quick detection of 16 KIR genes and the presence/absence of KIR-ligands based on allelic discrimination at codon 80 in the HLA-A/B Bw4 and HLA-C C1/C2 genes. These methods overcome the tedious and expensive nature of conventional KIR genotyping and HLA class I typing using sequence-specific primer (SSP) PCR, sequence-specific oligonucleotide (SSO) hybridization or sequence-based typing (SBT). Using these two cost-effective assays, we measured the frequencies of KIRs, KIR-ligands and KIR/KIR-ligand pairs in a cohort of Black women recruited in South Africa.

Enrichment of variations in KIR3DL1/S1 and KIR2DL2/L3 among H1N1/09 ICU patients: an exploratory study

Background

Infection by the pandemic influenza A (H1N1/09) virus resulted in significant pathology among specific ethnic groups worldwide. Natural Killer (NK) cells are important in early innate immune responses to viral infections. Activation of NK cells, in part, depend on killer-cell immunoglobulin-like receptors (KIR) and HLA class I ligand interactions. To study factors involved in NK cell dysfunction in overactive immune responses to H1N1 infection, KIR3DL1/S1 and KIR2DL2/L3 allotypes and cognate HLA ligands of H1N1/09 intensive-care unit (ICU) patients were determined.

Methodology and Findings

KIR3DL1/S1, KIR2DL2/L3, and HLA -B and -C of 51 H1N1/09 ICU patients and 105 H1N1-negative subjects (St. Theresa Point, Manitoba) were characterized. We detected an increase of 3DL1 ligand-negative pairs (3DL1/S1⁺ Bw6⁺ Bw4⁻), and a lack of 2DL1 HLA-C2 ligands, among ICU patients. They were also significantly enriched for 2DL2/L3 ligand-positive pairs (P<0.001, Pc<0.001; Odds Ratio:6.3158, CI95%:2.481–16.078). Relative to St. Theresa aboriginals (STh) and Venezuelan Amerindians (VA), allotypes enriched among aboriginal ICU patients (Ab) were: 2DL3 (Ab>VA, P = 0.024, Pc = 0.047; Odds Ratio:2.563, CI95%:1.109–5.923), 3DL1*00101 (Ab>VA, P<0.001, Pc<0.001), 3DL1*01502 (Ab>STh, P = 0.034, Pc = 0.268), and 3DL1*029 (Ab>STh, P  = 0.039, Pc = 0.301). Aboriginal patients ligand-positive for 3DL1/S1 and 2DL1 had the lowest probabilities of death (R_d) (R_d = 28%), compared to patients that were 3DL1/S1 ligand-negative (R_d = 52%) or carried 3DL1*029 (R_d = 52%). Relative to Caucasoids (CA), two allotypes were enriched among non-aboriginal ICU patients (NAb): 3DL1*00401 (NAb>CA, P<0.001, Pc<0.001) and 3DL1*01502 (CA<NAb, P = 0.012, Pc = 0.156). Non-aboriginal patients with ligands for all three KIRs (3DL1/S1, 2DL2/L3, and 2DL1) had the lowest probabilities of death (R_d = 36%), compared to subjects with 3DL1*01502 (R_d = 48%) and/or 3DL1*00401 (R_d = 58%).

Conclusions

Specific KIR3DL1/S1 allotypes, 3DL1/S1 and 2DL1 ligand-negative pairs, and 2DL2/L3 ligand-positive pairs were enriched among ICU patients. This suggests a possible association with NK cell dysfunction in patients with overactive immune responses to H1N1/09, leading to severe disease.

Expanded dog leukocyte antigen (DLA) single nucleotide polymorphism (SNP) genotyping reveals spurious class II associations

The dog leukocyte antigen (DLA) system contains many of the functional genes of the immune system, thereby making it a candidate region for involvement in immune-mediated disorders. A number of studies have identified associations between specific DLA class II haplotypes and canine immune hemolytic anemia, thyroiditis, immune polyarthritis, type I diabetes mellitus, hypoadrenocorticism, systemic lupus erythematosus-related disease complex, necrotizing meningoencephalitis (NME) and anal furunculosis. These studies have relied on sequencing approximately 300 bases of exon 2 of each of the DLA class II genes: DLA-DRB1, DLA-DQA1 and DLA-DQB1. In the present study, an association (odds ratio=4.29) was identified by this method between Weimaraner dogs with hypertrophic osteodystrophy (HOD) and DLA-DRB1∗01501. To fine map the association with HOD, a genotyping assay of 126 coding single nucleotide polymorphisms (SNPs) from across the entire DLA, spanning a region of 2.5 Mb (3,320,000-5,830,000) on CFA12, was developed and tested on Weimaraners with HOD, as well as two additional breeds with diseases associated with DLA class II: Nova Scotia duck tolling retrievers with hypoadrenocorticism and Pug dogs with NME. No significant associations were found between Weimaraners with HOD or Nova Scotia duck tolling retrievers with hypoadrenocorticism and SNPs spanning the DLA region. In contrast, significant associations were found with NME in Pug dogs, although the associated region extended beyond the class II genes. By including a larger number of genes from a larger genomic region, a SNP genotyping assay was generated that provides coverage of the extended DLA region and may be useful in identifying and fine mapping DLA associations in dogs.

Variable NK cell receptors exemplified by human KIR3DL1/S1

Variegated expression of variable NK cell receptors for polymorphic MHC class I broadens the range of an individual's NK cell response and the capacity for populations and species to survive disease epidemics and population bottlenecks. On evolutionary time scales, this component of immunity is exceptionally dynamic, unstable, and short-lived, being dependent on coevolution of ligands and receptors subject to varying, competing selection pressures. Consequently these systems of variable NK cell receptors are largely species specific and have recruited different classes of glycoprotein, even within the primate order of mammals. Such disparity helps to explain substantial differences in NK cell biology between humans and animal models, for which the population genetics is largely ignored. KIR3DL1/S1, which recognizes the Bw4 epitope of HLA-A and -B and is the most extensively studied of the variable NK cell receptors, exemplifies how variation in all possible parameters of function is recruited to diversify the human NK cell response.

Characterisation of Mhc class I and class II DRB polymorphism in red-bellied tamarins (Saguinus labiatus)

The infection of red-bellied tamarins (Saguinus labiatus) with GB virus B (GBV-B) is an important surrogate model of hepatitis C virus infection in man. To fully exploit the value of this model, we have characterised MHC class I G and class II DRB alleles in eight tamarins representing a cross-section of a UK breeding colony. The results indicated a high degree of classes I and II DRB allele sharing. Each animal transcribed three to four putative surface-expressed class I alleles and two to four class II DRB alleles. Most animals also transcribed at least one class I allele predicted to result in a C-terminal truncated protein. These results represent the first description of MHC polymorphism in this species and provide a foundation for characterisation of MHC diversity in breeding populations of red-bellied tamarins. The data will facilitate the identification of associations between MHC polymorphism and control of viral infections, and detailed dissection of cellular immune responses against GBV-B.

Progress of the comparative MHC Committee and a summary of the Comparative MHC Workshops held at the 32nd ISAG, Edinburgh and the 9th IVIS, Tokyo, 2010

The comparative MHC Nomenclature Committee was formed following a recommendation by the International Society for Animal Genetics (ISAG) at its meeting in Göttingen, Germany in 2002. The committee subsequently became affiliated to VIC, the Veterinary Immunology Committee of IUIS; the International Union of Immunological Societies, which is the umbrella organization for immunology societies worldwide. The committee currently fulfils three roles: (1) to establish guidelines for MHC nomenclature in species outside of human and mouse, (2) to facilitate the development of databases which provide a source of validated MHC sequences with official nomenclatures for associated species (immuno-polymorphism database IPD-MHC) and (3) to facilitate scientific discussion on comparative aspects of MHC biology. The committee currently comprises individuals representing a range of species nomenclature committees. This report summarises progress the committee has made in addressing its core objectives as well as the principal outputs from the Comparative MHC Workshops held at the 32nd International Society for Animal Genetics (ISAG) conference in Edinburgh, UK in July 2010 as well as in August 2010 at the 9th International Veterinary Immunology Symposium (IVIS) in Tokyo, Japan. Additional data from posters and presentations at related workshops focussing on aspects of MHC biology are also included where appropriate.

Major histocompatibility complex variation in insular populations of the Egyptian vulture: inferences about the roles of genetic drift and selection

Insular populations have attracted the attention of evolutionary biologists because of their morphological and ecological peculiarities with respect to their mainland counterparts. Founder effects and genetic drift are known to distribute neutral genetic variability in these demes. However, elucidating whether these evolutionary forces have also shaped adaptive variation is crucial to evaluate the real impact of reduced genetic variation in small populations. Genes of the major histocompatibility complex (MHC) are classical examples of evolutionarily relevant loci because of their well-known role in pathogen confrontation and clearance. In this study, we aim to disentangle the partial roles of genetic drift and natural selection in the spatial distribution of MHC variation in insular populations. To this end, we integrate the study of neutral (22 microsatellites and one mtDNA locus) and MHC class II variation in one mainland (Iberia) and two insular populations (Fuerteventura and Menorca) of the endangered Egyptian vulture (Neophron percnopterus). Overall, the distribution of the frequencies of individual MHC alleles (n=17 alleles from two class II B loci) does not significantly depart from neutral expectations, which indicates a prominent role for genetic drift over selection. However, our results point towards an interesting co-evolution of gene duplicates that maintains different pairs of divergent alleles in strong linkage disequilibrium on islands. We hypothesize that the co-evolution of genes may counteract the loss of genetic diversity in insular demes, maximize antigen recognition capabilities when gene diversity is reduced, and promote the co-segregation of the most efficient allele combinations to cope with local pathogen communities.

Genetic polymorphisms of killer cell immunoglobulin-like receptor 3DL2 in preeclampsia

AIMS

We evaluated if killer cell immunoglobulin-like receptor 3DL2 gene (KIR3DL2) polymorphisms are a key factor in the development of preeclampsia.

METHODS

In this case-control study, 105 pregnant women with PE (PE group) were enrolled. Their A52G in exon 3 and C32T in exon 9 polymorphisms of the KIR3DL2 genotypes were determined by polymorphism chain reaction-restriction fragment length polymorphism (PCR-RFLP) from venous blood samples and compared with the corresponding KIR3DL2 genotypes of 103 pregnant women with uncomplicated pregnancies (control group).

RESULTS

Carriers of the A allele in exon 3 of KIR3DL2 gene occurred less frequently in PE than in controls [P=0.001; odds ratio (OR)=2.65, range: 1.5-4.7]. No significant difference was found about allelic frequencies of KIR3DL2 gene C32T in exon 9 in women with preeclampsia as compared to controls. A significant difference between the two groups of genotypic frequencies of KIR3DL2 gene A52G in exon 3 and KIR3DL2 gene C32T in exon 9 polymorphisms was found (P=0.003 and P=0.000). There was no significant difference between genotypic or allelic frequencies in women with mild preeclampsia compared to sever preeclampsia.

CONCLUSIONS

Our results suggest that carriers of A allele in exon 3 have a decreased susceptibility to PE. It is likely that the presence of the CC genotype in exon 9 has a considerable effect on disease progression. The mutation of the two sites is not associated with the severity of preeclampsia.

Characterization of the major histocompatibility complex class II DOB, DPB1, and DQB1 alleles in cynomolgus macaques of Vietnamese origin

Major histocompatibility complex (MHC) molecules play an important role in the susceptibility and/or resistance to many diseases. To gain an insight into the MHC background and to facilitate the experimental use of cynomolgus macaques, the second exon of the MhcMafa-DOB, -DPB1, and -DQB1 genes from 143 cynomolgus macaques were characterized by cloning to sequencing. A total of 16 Mafa-DOB, 16 Mafa-DPB1, and 34 Mafa-DQB1 alleles were identified, which revealed limited, moderate, and marked allelic polymorphism at DOB, DPB1, and DQB1, respectively, in a cohort of cynomolgus macaques of Vietnamese origin. In addition, 16 Mafa-DOB, 5 Mafa-DPB1, and 8 Mafa-DQB1 alleles represented novel sequences that had not been reported in earlier studies. Almost of the sequences detected at the DOB and DQB1 locus in the present study belonged to DOB*01 (100%) and DQB1*06 (62%) lineages, respectively. Interestingly, four, three, and one high-frequency alleles were detected at Mafa-DOB, -DPB1, and -DQB1, respectively, in this monkeys. The alleles with the highest frequency among these monkeys were Mafa-DOB*010102, Mafa-DPB1*13, and Mafa-DQB1*0616, and these were found in 33 (25.6%) of 129 monkeys, 32 (31.37%) of 102 monkeys, and 30 (31%) of 143 monkeys, respectively. The high-frequency alleles may represent high priority targets for additional characterization of immune function. We also carried out evolutionary and population analyses using these sequences to reveal population-specific alleles. This information will not only promote the understanding of MHC diversity and polymorphism in the cynomolgus macaque but will also increase the value of this species as a model for biomedical research.

A single nomenclature and associated database for alleles at the major histocompatibility complex class II DRB1 locus of sheep

The development of standardised nomenclatures with associated databases containing reference sequences for alleles at polymorphic loci within the major histocompatibility complex (MHC) has been facilitated by the development of the immuno polymorphism database (IPD). Recently, included within IPD-MHC is information on allelic diversity within sheep species (IPD-MHC-OLA). Here, we present the first report of progress in populating the sheep IPD-MHC database with alleles at the class II MHC DRB1 locus. The sequence of 63 Ovar-DRB1 alleles within 24 allelic families is now held within the database, each meeting the minimum requirement of a complete second exon. These sequences are derived from a combination of genomic and cDNA-based approaches and represent the most extensive collection of validated alleles at the sheep DRB1 locus yet described. Although these 63 alleles probably represent only a fraction of the DRB1 allelic diversity in sheep species worldwide, we encourage the research community to use the official allelic nomenclature and to contribute allelic sequences to the database via its web-based submission tool. In time, the IPD-MHC-OLA resource will underpin population-based MHC genotyping studies and help to simplify meta-analyses of multi-source data from wild and domestic sheep populations.

An allelic typing method for 2DS4 variant used in study of haplotypes of killer cell immunoglobulin-like receptor gene

The KIR2DS4 variants differ in exon 5 and play a role in hematopoietic stem cells transplantation (HSCT). A sequence-based testing (SBT) and TOPO TA cloning system identifying and distinguishing alleles of the KIR2DS4 gene was established and applied to a total of 150 Chinese-Han individuals: 75 patients received T-cell-depleted HSCT and their unrelated donors. The majority (139) of the 150 samples (92.7%) were positive for KIR2DS4. Four of the nine known KIR2DS4 alleles, KIR2DS4 *00101, *003,*004, and *007, were identified. In the haplotype A/A group, a higher risk of acute graft-versus-host disease (aGVHD) was seen when the donor carried two full-length KIR2DS4 alleles (RR 9.0 [95% CI 1.2-66.9], P = 0.010). Our findings suggested that the expression of full-length 2DS4 (*001) in A/A group may contribute to a worse clinical outcome after URD-HSCT. These data would be beneficial for the selection of suitable donors.

Interactions of NK cell receptor KIR3DL1*004 with chaperones and conformation-specific antibody reveal a functional folded state as well as predominant intracellular retention

Variable interaction between the Bw4 epitope of HLA-B and the polymorphic KIR3DL1/S1 system of inhibitory and activating NK cell receptors diversifies the development, repertoire formation, and response of human NK cells. KIR3DL1*004, a common KIR3DL1 allotype, in combination with Bw4(+) HLA-B, slows progression of HIV infection to AIDS. Analysis in this study of KIR3DL1*004 membrane traffic in NK cells shows this allotype is largely misfolded but stably retained in the endoplasmic reticulum, where it binds to the chaperone calreticulin and does not induce the unfolded protein response. A small fraction of KIR3DL1*004 folds correctly and leaves the endoplasmic reticulum to be expressed on the surface of primary NK and transfected NKL cells, in a form that can be triggered to inhibit NK cell activation and secretion of IFN-γ. Consistent with this small proportion of correctly folded molecules, trace amounts of MHC class I coimmunoprecipitated with KIR3DL1*004. There was no indication of any extensive intracellular interaction between unfolded KIR3DL1*004 and cognate Bw4(+) HLA-B. A similarly limited interaction of Bw4 with KIR3DL1*002, when both were expressed by the same cell, was observed despite the efficient folding of KIR3DL1*002 and its abundance on the NK cell surface. Several positions of polymorphism modulate KIR3DL1 abundance at the cell surface, differences that do not necessarily correlate with the potency of allotype function. In this context, our results suggest the possibility that the effect of Bw4(+) HLA-B and KIR3DL1*004 in slowing progression to AIDS is mediated by interaction of Bw4(+) HLA-B with the small fraction of cell surface KIR3DL1*004.