The MHC, disease and selection

Functional variation at an expressed MHC class IIβ locus associates with Ranavirus infection intensity in larval anuran populations

Infectious diseases are causing catastrophic losses to global biodiversity. Iridoviruses in the genus Ranavirus are among the leading causes of amphibian disease-related mortality. Polymorphisms in major histocompatibility complex (MHC) genes are significantly associated with variation in amphibian pathogen susceptibility. MHC genes encode two classes of polymorphic cell-surface molecules that can recognize and bind to diverse pathogen peptides. While MHC class I genes are the classic mediators of viral-acquired immunity, larval amphibians do not express them. Consequently, MHC class II gene diversity may be an important predictor of Ranavirus susceptibility in larval amphibians, the life stage most susceptible to Ranavirus. We surveyed natural populations of larval wood frogs (Rana sylvatica), which are highly susceptible to Ranavirus, across 17 ponds and 2 years in Maryland, USA. We sequenced the peptide-binding region of an expressed MHC class IIβ locus and assessed allelic and genetic diversity. We converted alleles to functional supertypes and determined if supertypes or alleles influenced host responses to Ranavirus. Among 381 sampled individuals, 26% were infected with Ranavirus. We recovered 20 unique MHC class IIβ alleles that fell into two deeply diverged clades and seven supertypes. MHC genotypes were associated with Ranavirus infection intensity, but not prevalence. Specifically, MHC heterozygotes and supertype ST1/ST7 had significantly lower Ranavirus infection intensity compared to homozygotes and other supertypes. We conclude that MHC class IIβ functional genetic variation is an important component of Ranavirus susceptibility. Identifying immunogenetic signatures linked to variation in disease susceptibility can inform mitigation strategies for combatting global amphibian declines.

Pathogen diversity drives the evolution of generalist MHC-II alleles in human populations

Central players of the adaptive immune system are the groups of proteins encoded in the major histocompatibility complex (MHC), which shape the immune response against pathogens and tolerance to self-peptides. The corresponding genomic region is of particular interest, as it harbors more disease associations than any other region in the human genome, including associations with infectious diseases, autoimmune disorders, cancers, and neuropsychiatric diseases. Certain MHC molecules can bind to a much wider range of epitopes than others, but the functional implication of such an elevated epitope-binding repertoire has remained largely unclear. It has been suggested that by recognizing more peptide segments, such promiscuous MHC molecules promote immune response against a broader range of pathogens. If so, the geographical distribution of MHC promiscuity level should be shaped by pathogen diversity. Three lines of evidence support the hypothesis. First, we found that in pathogen-rich geographical regions, humans are more likely to carry highly promiscuous MHC class II DRB1 alleles. Second, the switch between specialist and generalist antigen presentation has occurred repeatedly and in a rapid manner during human evolution. Third, molecular positions that define promiscuity level of MHC class II molecules are especially diverse and are under positive selection in human populations. Taken together, our work indicates that pathogen load maintains generalist adaptive immune recognition, with implications for medical genetics and epidemiology.

Whereas specialist major histocompatibility complex (MHC) molecules initiate immune response against only relatively few pathogens, generalists provide protection against a broad range. Accordingly, this study shows that the geographical distribution of generalist MHC alleles in human populations reflects exposure to diverse infectious diseases.

Variation in the human genome influences our susceptibility to infectious diseases, but the causal link between disease and underlying mutation often remains enigmatic. Major histocompatibility complex II (MHC class II) molecules shape both our immune response against pathogens and our tolerance of self-peptides. The genomic region that encodes MHC molecules is of particular interest, as it is home to more genetic disease associations than any other region in the human genome, including associations with infectious diseases, autoimmune disorders, cancers, and neuropsychiatric diseases. Here, we propose that MHC class II molecules can be categorized into two major types; specialists initiate effective immune response against only relatively few pathogens, while generalists provide protection against a broad range of pathogens. As support, we demonstrate that generalist MHC class II variants are more prevalent in human populations residing in pathogen-rich areas.

The MHC in the era of next-generation sequencing: Implications for bridging structure with function

The MHC continues to have the most disease-associations compared to other regions of the human genome, even in the genome-wide association study (GWAS) and single nucleotide polymorphism (SNP) era. Analysis of non-coding variation and their impact on the level of expression of HLA allotypes has shed new light on the potential mechanisms underlying HLA disease associations and alloreactivity in transplantation. Next-generation sequencing (NGS) technology has the capability of delineating the phase of variants in the HLA antigen-recognition site (ARS) with non-coding regulatory polymorphisms. These relationships are critical for understanding the qualitative and quantitative implications of HLA gene diversity. This article summarizes current understanding of non-coding region variation of HLA loci, the consequences of regulatory variation on HLA expression, the role for evolution in shaping lineage-specific expression, and the impact of HLA expression on disease susceptibility and transplantation outcomes. A role for phased sequencing methods for the MHC, and perspectives for future directions in basic and applied immunogenetic studies of the MHC are presented.

Allele-Level KIR Genotyping of More Than a Million Samples: Workflow, Algorithm, and Observations

The killer-cell immunoglobulin-like receptor (KIR) genes regulate natural killer cell activity, influencing predisposition to immune mediated disease, and affecting hematopoietic stem cell transplantation (HSCT) outcome. Owing to the complexity of the KIR locus, with extensive gene copy number variation (CNV) and allelic diversity, high-resolution characterization of KIR has so far been applied only to relatively small cohorts. Here, we present a comprehensive high-throughput KIR genotyping approach based on next generation sequencing. Through PCR amplification of specific exons, our approach delivers both copy numbers of the individual genes and allelic information for every KIR gene. Ten-fold replicate analysis of a set of 190 samples revealed a precision of 99.9%. Genotyping of an independent set of 360 samples resulted in an accuracy of more than 99% taking into account consistent copy number prediction. We applied the workflow to genotype 1.8 million stem cell donor registry samples. We report on the observed KIR allele diversity and relative abundance of alleles based on a subset of more than 300,000 samples. Furthermore, we identified more than 2,000 previously unreported KIR variants repeatedly in independent samples, underscoring the large diversity of the KIR region that awaits discovery. This cost-efficient high-resolution KIR genotyping approach is now applied to samples of volunteers registering as potential donors for HSCT. This will facilitate the utilization of KIR as additional selection criterion to improve unrelated donor stem cell transplantation outcome. In addition, the approach may serve studies requiring high-resolution KIR genotyping, like population genetics and disease association studies.

Allelic nomenclature for the duplicated MHC class II DQ genes in sheep

The principal MHC class II molecules involved in the presentation of peptides to the antigen specific receptors on CD4⁺ T cells genes in sheep are derived from DR and DQ genes. Allelic nomenclature systems for the DRB1 and its partner DRA loci are available for Ovid's; however, no official nomenclature is available for the DQ genes which creates ambiguity within the research community. Ovine MHC haplotypes include at least two pairs of DQA and DQB genes, termed DQA1, DQB1 and DQA2, DQB2 and both sets are polymorphic and both seem to be functional. In a number of haplotypes, the DQA1 locus appears to be absent (DQA1-null) and is replaced by a second locus termed DQA2-like. Here, we identify families of alleles based on sequence similarity and phylogenetic clustering which correspond to each of the DQA and DQB genes identified in previous genomic and transcript analyses of homozygous animals. Using such criteria to cluster sequences, we have named 82 full-length and partial cDNA transcripts derived from domestic sheep (Ovis aries) which correspond to alleles at the Ovar-DQA1, DQA2, DQA2-like, DQB1, DQB2 and DQB2-like genes and provide associated sequence resources available to the research community through the IPD-MHC Database. This sets the basis for naming and annotation of DQ genes within the ovine MHC and may be used as a template for DQ genes in other ruminant species which will ultimately support research in livestock infectious disease.

HLA genetics in bipolar disorder

OBJECTIVE

Bipolar Disorder (BD) is characterized by deregulated adaptive immune processes. Recent genome-wide association studies (GWAS) implicate the major histocompatibility complex (MHC) region in BD. The present study investigates the potential influence of variations in human leukocyte antigen (HLA) on BD risk and/or clinical presentations. This may have relevance to the dysregulated inflammatory processes commonly found in BD.

METHOD

DNAs from 475 BD patients and 195 healthy controls (HC) were genotyped for classical HLA class I and II loci.

RESULTS

We found that: (i) the HLA-A*02~B*44~DRB1*07 sub-haplotype is less prevalent in BD, vs. HC (pc = 2.4 × 10^-2 ); (ii) the 57.1 and the 8.1-derived ancestral haplotypes i.e. HLA-A*02~B*57~Cw*06~DRB1*07~DQB1*09 and HLA-A*02~B*08~Cw*07 are associated with rapid cycling (pc = 1.9 × 10^-3 and 1.05 × 10^-2 , respectively); (iii) the 8.1AH-derived HLA class II-DRB*03~HLA-DQB1*02 sub-haplotype is more frequent in BD patients with a history of suicidal behaviors (pc = 2.1 × 10^-2 ); and (iv) disease onset by an hypomanic episode or by psychotic symptoms are, respectively, more frequent in BD patients bearing the 7.1 AH-derived A*03~B*07~DRB1*15 sub-haplotype (pc = 8.5 × 10^-3 ) and the HLA-A*02~B*07~DRB1*15 sub-haplotype (pc = 4.0 × 10^-2 ).

CONCLUSION

Corroborating the established link between these HLA haplotypes/sub haplotypes and common immune disorders, our findings suggest possible HLA-mediated proinflammatory processes operating in BD.

Microbial modulation of the gut microbiome for treating autoimmune diseases

Many studies have shown the relationship between autoimmune diseases and the gut microbiome in humans: those with autoimmune conditions display gut microbiome dysbiosis. The big question that needs to be addressed is if restoring eubiosis of the gut microbiota can help suppress the autoimmune condition by activating various immune regulatory mechanisms. Inducing these self-healing mechanisms should prolong good health in affected individuals. Area covered: Here, we review the available clinical and preclinical studies that have used selective bacteria for modulating gut microbiota for treating autoimmune diseases. The potential bacterial candidates and their mechanism of action in treating autoimmune diseases will be discussed. We searched for genetically modified and potential probiotics for diseases and discuss the most likely candidates. Expert commentary: To achieve eubiosis, manipulation of the gut microbiota must occur in some form. Several approaches for modulating gut microbiota include prebiotic diets, antimicrobial interventions, fecal microbiota transplants, and selective probiotics. One novel approach showing promising results is the use of selective bacterial candidates to modulate microbial composition. Use of single microbe for treatment has an advantage as compared to multi-species as microbes grow at different rates and if needed, a single microbe is easy to target.

Fine mapping of MHC region in lung cancer highlights independent susceptibility loci by ethnicity

The basis for associations between lung cancer and major histocompatibility complex genes is not completely understood. Here the authors further consider genetic variation within the MHC region in lung cancer patients and identify independent associations within HLA genes that explain MHC lung cancer associations in Europeans and Asian populations.

Lung cancer has several genetic associations identified within the major histocompatibility complex (MHC); although the basis for these associations remains elusive. Here, we analyze MHC genetic variation among 26,044 lung cancer patients and 20,836 controls densely genotyped across the MHC, using the Illumina Illumina OncoArray or Illumina 660W SNP microarray. We impute sequence variation in classical HLA genes, fine-map MHC associations for lung cancer risk with major histologies and compare results between ethnicities. Independent and novel associations within HLA genes are identified in Europeans including amino acids in the HLA-B*0801 peptide binding groove and an independent HLA-DQB1*06 loci group. In Asians, associations are driven by two independent HLA allele sets that both increase risk in HLA-DQB1*0401 and HLA-DRB1*0701; the latter better represented by the amino acid Ala-104. These results implicate several HLA–tumor peptide interactions as the major MHC factor modulating lung cancer susceptibility.

Architecture of polymorphisms in the human genome reveals functionally important and positively selected variants in immune response and drug transporter genes

Background

Genetic polymorphisms can contribute to phenotypic differences amongst individuals, including disease risk and drug response. Characterization of genetic polymorphisms that modulate gene expression and/or protein function may facilitate the identification of the causal variants. Here, we present the architecture of genetic polymorphisms in the human genome focusing on those predicted to be potentially functional/under natural selection and the pathways that they reside.

Results

In the human genome, polymorphisms that directly affect protein sequences and potentially affect function are the most constrained variants with the lowest single-nucleotide variant (SNV) density, least population differentiation and most significant enrichment of rare alleles. SNVs which potentially alter various regulatory sites, e.g. splicing regulatory elements, are also generally under negative selection.

Interestingly, genes that regulate the expression of transcription/splicing factors and histones are conserved as a higher proportion of these genes is non-polymorphic, contain ultra-conserved elements (UCEs) and/or has no non-synonymous SNVs (nsSNVs)/coding INDELs. On the other hand, major histocompatibility complex (MHC) genes are the most polymorphic with SNVs potentially affecting the binding of transcription/splicing factors and microRNAs (miRNA) exhibiting recent positive selection (RPS). The drug transporter genes carry the most number of potentially deleterious nsSNVs and exhibit signatures of RPS and/or population differentiation. These observations suggest that genes that interact with the environment are highly polymorphic and targeted by RPS.

Conclusions

In conclusion, selective constraints are observed in coding regions, master regulator genes, and potentially functional SNVs. In contrast, genes that modulate response to the environment are highly polymorphic and under positive selection.

Electronic supplementary material

The online version of this article (10.1186/s40246-018-0175-1) contains supplementary material, which is available to authorized users.

Employing proteomics in the study of antigen presentation: an update

INTRODUCTION

Our immune system discriminates self from non-self by examining the peptide cargo of human leukocyte antigen (HLA) molecules displayed on the cell surface. Successful recognition of HLA-bound non-self peptides can induce T cell responses leading to, for example, the destruction of infected cells. Today, largely due to advances in technology, we have an unprecedented capability to identify the nature of these presented peptides and unravel the true complexity of antigen presentation. Areas covered: In addition to conventional linear peptides, HLA molecules also present post-translationally modified sequences comprising a wealth of chemical and structural modifications, including a novel class of noncontiguous spliced peptides. This review focuses on these emerging themes in antigen presentation and how mass spectrometry in particular has contributed to a new view of the antigenic landscape that is presented to the immune system. Expert Commentary: Advances in the sensitivity of mass spectrometers and use of hybrid fragmentation technologies will provide more information-rich spectra of HLA bound peptides leading to more definitive identification of T cell epitopes. Coupled with improvements in sample preparation and new informatics workflows, studies will access novel classes of peptide antigen and allow interrogation of rare and clinically relevant samples.

Single nucleotide polymorphisms within HLA region are associated with disease relapse for patients with unrelated cord blood transplantation

Disease relapse occurs in unrelated cord blood transplantation (CBT) even when the alleles of human leukocyte antigen (HLA) are fully matched between donor and recipient. This is similar to that observed in other types of hematopoietic stem cell transplantation. Fourteen single nucleotide polymorphisms (SNPs) within the HLA region have been reported previously by Petersdorf et al. and Piras et al. as transplantation determinants in unrelated hematopoietic cell transplantation. In this study, the genomic sequences within 500 base pairs upstream and downstream of the fourteen transplantation-related SNPs from 53 patients and their HLA-matched unrelated donors were analyzed for determining whether or not genetic variants, conferred by either recipient or donor SNP genotype or by recipient-donor SNP mismatching, were associated with the risk of relapse. Seven SNPs were associated with the risk of relapse in unrelated CBT. These included the donor genotype with the SNPs of rs2523675 and rs2518028 at the telomeric end of HCP5 gene, rs2071479 in the intron of the HLA-DOB gene, and rs2523958 in the MICD gene; and the recipient genotype with SNPs of rs9276982 in the HLA-DOA gene, and rs435766 and rs380924 in the MICD gene. As measured by pair-wise linkage disequilibrium (LD) with D′ as the parameter for normalized standard measurement of LD which compares the observed and expected frequencies of one haplotype comprised by alleles at different loci, rs2523675 had high LD with rs4713466 (D′ = 0.86) and rs2523676 (D′ = 0.91) in the HCP5 gene. The rs2518028 had no LD with all other SNPs except rs2523675 (D′ = 0.76). This study provides the basis for developing a method or algorithm for selecting better unrelated CBT candidate donors.

Peptide presentation by HLA-DQ molecules is associated with the development of immune tolerance

HLA class II proteins are important elements of human adaptive immune recognition and are associated with numerous infectious and immune-mediated diseases. These highly variable molecules can be classified into DP, DQ and DR groups. It has been proposed that in contrast with DP and DR, epitope binding by DQ variants rather results in immune tolerance. However, the pieces of evidence are limited and controversial. We found that DQ molecules bind more human epitopes than DR. Pathogen-associated epitopes bound by DQ molecules are more similar to human proteins than the ones bound by DR. Accordingly, DQ molecules bind epitopes of significantly different pathogen species. Moreover, the binding of autoimmunity-associated epitopes by DQ confers protection from autoimmune diseases. Our results suggest a special role of HLA-DQ in immune homeostasis and help to better understand the association of HLA molecules with infectious and autoimmune diseases.

HLA-class II haplotypes and Autism Spectrum Disorders

Infections and autoimmunity are associated with autism spectrum disorders (ASD), with both strongly influenced by the genetic regulation of the human leukocyte antigen (HLA) system. The relationship between ASD and the HLA genetic diversity requires further investigation. Using a case control design, the distribution of HLA class II-DRB1 and DQB1 alleles, genotypes and haplotypes were investigated in ASD patients, versus healthy controls (HC). ASD patients meeting DSM-IV TR criteria and HC (474 and 350 respectively) were genotyped at medium resolution using a Luminex-based SSO technology. Comparisons of genotypes, allele frequencies associated with a haplotype analysis were performed. Results indicate: (i) the HLA-DRB1 *11-DQB1*07 haplotype was more prevalent in ASD patients, versus HC (Pc = 0.001), partially replicating previous data and possibly linking to gastro-intestinal (GI)-related pro-inflammatory processes, given that this haplotype associates with pediatric celiac disorders; (ii) the HLA-DRB1 *17-DQB1*02 haplotype was higher in HC, versus ASD patients (Pc = 0.002), indicating that this is a protective haplotype. Using the Autism Diagnostic Interview to assess clinical dimensions, higher scores on social (Pc = 0.006) and non-verbal functioning (Pc = 0.004) associated with the DRB1 *11 DQB1*07 haplotype. Our results support HLA involvement in ASD, with possible relevance to GI and gut-brain axis dysregulation.

Evolution of DNAase I Hypersensitive Sites in MHC Regulatory Regions of Primates

It has been challenging to determine the disease-causing variant(s) for most major histocompatibility complex (MHC)-associated diseases. However, it is becoming increasingly clear that regulatory variation is pervasive and a fundamentally important mechanism governing phenotypic diversity and disease susceptibility. We gathered DNase I data from 136 human cells to characterize the regulatory landscape of the MHC region, including 4867 DNase I hypersensitive sites (DHSs). We identified thousands of regulatory elements that have been gained or lost in the human or chimpanzee genomes since their evolutionary divergence. We compared alignments of the DHS across six primates and found 149 DHSs with convincing evidence of positive and/or purifying selection. Of these DHSs, compared to neutral sequences, 24 evolved rapidly in the human lineage. We identified 15 instances of transcription-factor-binding motif gains, such as USF, MYC, MAX, MAFK, STAT1, PBX3, etc, and observed 16 GWAS (genome-wide association study) SNPs associated with diseases within these 24 DHSs using FIMO (Find Individual Motif Occurrences) and UCSC (University of California, Santa Cruz) ChIP-seq data. Combining eQTL and Hi-C data, our results indicated that there were five SNPs located in human gains motifs affecting the corresponding gene's expression, two of which closely matched DHS target genes. In addition, a significant SNP, rs7756521, at genome-wide significant level likely affects DDR expression and represents a causal genetic variant for HIV-1 control. These results indicated that species-specific motif gains or losses of rapidly evolving DHSs in the primate genomes might play a role during adaptation evolution and provided some new evidence for a potentially causal role for these GWAS SNPs.

Kidney Transplantation: The Challenge of Human Leukocyte Antigen and Its Therapeutic Strategies

Kidney transplantation remains the treatment of choice for end-stage renal failure. When the immune system of the recipient recognizes the transplanted kidney as a foreign object, graft rejection occurs. As part of the host immune defense mechanism, human leukocyte antigen (HLA) is a major challenge for graft rejection in transplantation therapy. The impact of HLA mismatches between the donor and the potential recipient prolongs the time for renal transplantation therapy, tethered to dialysis, latter reduces graft survival, and increases mortality. The formation of pretransplant alloantibodies against HLA class I and II molecules can be sensitized through exposures to blood transfusions, prior transplants, and pregnancy. These preformed HLA antibodies are associated with rejection in kidney transplantation. On the other hand, the development of de novo antibodies may increase the risk for acute and chronic rejections. Allograft rejection results from a complex interplay involving both the innate and the adaptive immune systems. Thus, further insights into the mechanisms of tissue rejection and the risk of HLA sensitization is crucial in developing new therapies that may blunt the immune system against transplanted organs. Therefore, the purpose of this review is to highlight facts about HLA and its sensitization, various mechanisms of allograft rejection, the current immunosuppressive approaches, and the directions for future therapy.

The Intergenic Recombinant HLA-B∗46:01 Has a Distinctive Peptidome that Includes KIR2DL3 Ligands

HLA-B^∗46:01 was formed by an intergenic mini-conversion, between HLA-B^∗15:01 and HLA-C^∗01:02, in Southeast Asia during the last 50,000 years, and it has since become the most common HLA-B allele in the region. A functional effect of the mini-conversion was introduction of the C1 epitope into HLA-B^∗46:01, making it an exceptional HLA-B allotype that is recognized by the C1-specific natural killer (NK) cell receptor KIR2DL3. High-resolution mass spectrometry showed that HLA-B^∗46:01 has a low-diversity peptidome that is distinct from those of its parents. A minority (21%) of HLA-B^∗46:01 peptides, with common C-terminal characteristics, form ligands for KIR2DL3. The HLA-B^∗46:01 peptidome is predicted to be enriched for peptide antigens derived from Mycobacterium leprae. Overall, the results indicate that the distinctive peptidome and functions of HLA-B^∗46:01 provide carriers with resistance to leprosy, which drove its rapid rise in frequency in Southeast Asia.

•

The interlocus recombinant HLA-B∗46:01 is found at high frequency in Southeast Asia

•

HLA-B∗46:01 has a low-diversity peptidome that is distinct from both its parents

•

A subset of HLA-B∗46:01 peptides provides ligands for the NK cell receptor KIR2DL3

•

The unique features of HLA-B∗46:01 correlate with protection against leprosy

Microsatellite markers for evaluating the diversity of the natural killer complex and major histocompatibility complex genomic regions in domestic horses

Genotyping microsatellite markers represents a standard, relatively easy, and inexpensive method of assessing genetic diversity of complex genomic regions in various animal species, such as the major histocompatibility complex (MHC) and/or natural killer cell receptor (NKR) genes. MHC-linked microsatellite markers have been identified and some of them were used for characterizing MHC polymorphism in various species, including horses. However, most of those were MHC class II markers, while MHC class I and III sub-regions were less well covered. No tools for studying genetic diversity of NKR complex genomic regions are available in horses. Therefore, the aims of this work were to establish a panel of markers suitable for analyzing genetic diversity of the natural killer complex (NKC), and to develop additional microsatellite markers of the MHC class I and class III genomic sub-regions in horses. Nine polymorphic microsatellite loci were newly identified in the equine NKC. Along with two previously reported microsatellites flanking this region, they constituted a panel of 11 loci allowing to characterize genetic variation in this functionally important part of the horse genome. Four newly described MHC class I/III-linked markers were added to 11 known microsatellites to establish a panel of 15 MHC markers with a better coverage of the class I and class III sub-regions. Major characteristics of the two panels produced on a group of 65 horses of 13 breeds and on five Przewalski's horses showed that they do reflect genetic variation within the horse species.

In silico prediction of potential vaccine candidates on capsid protein of human bocavirus 1

Human bocavirus 1 (HBoV1) is a newly identified parvovirus that causes serious respiratory infection among children across the globe. Aim of the present study was to predict immunogenic residues located on the VP2 protein of HBoV1 towards development of epitope based vaccines. Several computational tools were employed to predict epitopes (bothT and B cell restricted) with stringent regulation for the improvement of confidence. After meticulous analysis, the peptide "TTPWTYFNFNQY" was identified as potential candidate for development of preventive vaccine. Of note, the epitope "TTPWTYFNFNQY" was found to be recognized by fifteen different alleles belonging to seven HLA supertypes (A1, A3, A24, A26, B7, B58 and B62). Further, mutational variability analysis pointed that most of the amino acids were well conserved. Docking scores obtained from ClusPro and Autodock Vina for selected epitopes displayed energetically favorable and stable interaction of peptide-HLA-I complexes. The core peptide "LLYQMPFFL" was found to recognize by wide range of HLA class II allele recognition thereby qualified as candidate for therapeutic vaccine. Five distinct linear peptides (withT cell epitope superimposition) belonging to B cells were identified in the VP2 protein. Further attention on the enlisted epitopes may shed light on the path for development of diagnostic, therapeutic and preventive tools against HBoV1 infection. Additionally, the predicted epitopes may help us to address the original antigenic sin phenomena observed during consecutive HBoV2-4 infection.

Genetic factors are stressed variably by onset age-based sample selection in psoriasis: A hint from major histocompatibility complex region-based analysis

BACKGROUND

Large cohort-based genetic association studies have been established over a decade. However, for certain diseases, different results with respect to the genome-wide association study level have been obtained among studies, even for those conducted within the same ethnic groups. We hypothesized that onset age-based sample variables might have a great impact on the results.

METHODS

In the present study, we divided psoriasis patients into several subgroups according to the onset age bracket. We conducted genetic association analysis in the major histocompatibility complex (MHC) region of each patient subgroup with shared control subjects.

RESULTS

We found decreases in the numbers of susceptible variants in each subgroup analysis as the onset age increased in the longitudinal analysis. Meanwhile, the pairwise analysis showed that younger patients exhibited greater numbers of genetic risks in the MHC region compared to elder patients, regardless of whether the cut-off values were defined as 20 or 30 years old. Similar results were also found among 11-20-, 21-30- and 31-40-year-old groups. Furthermore, when combining the results of both the stepwise regression analysis and the HLA-C*06:02 conditioning analysis, different variants were found to be independently associated with each psoriasis subgroup.

CONCLUSIONS

Onset age-based sample variables influence the results of genetic association studies, at least in MHC region-based genetic analysis. We suggest that caution is required when selecting samples for genetic association studies to prevent confounders that might be a result of onset age.

Review of pharmacogenetics studies of L-asparaginase hypersensitivity in acute lymphoblastic leukemia points to variants in the GRIA1 gene

Acute lymphoblastic leukemia (ALL) is a major pediatric cancer in developed countries. Although treatment outcome has improved owing to advances in chemotherapy, there is still a group of patients who experience severe adverse events. L-Asparaginase is an effective antineoplastic agent used in chemotherapy of ALL. Despite its indisputable indication, hypersensitivity reactions are common. In those cases, discontinuation of treatment is usually needed and anti-asparaginase antibody production may also attenuate asparaginase activity, compromising its antileukemic effect. Till now, six pharmacogenetic studies have been performed in order to elucidate possible genetic predisposition for inter-individual differences in asparaginase hypersensitivity. In this review we have summarized the results of those studies which describe the involvement of four different genes, being polymorphisms in the glutamate receptor, ionotropic, AMPA 1 (GRIA1) the most frequently associated with asparaginase hypersensitivity. We also point to new approaches focusing on epigenetics that could be interesting for consideration in the near future.

Evolutionary dynamics of an expressed MHC class IIβ locus in the Ranidae (Anura) uncovered by genome walking and high-throughput amplicon sequencing

The Major Histocompatibility Complex (MHC) is a genomic region encoding immune loci that are important and frequently used markers in studies of adaptive genetic variation and disease resistance. Given the primary role of infectious diseases in contributing to global amphibian declines, we characterized the hypervariable exon 2 and flanking introns of the MHC Class IIβ chain for 17 species of frogs in the Ranidae, a speciose and cosmopolitan family facing widespread pathogen infections and declines. We find high levels of genetic variation concentrated in the Peptide Binding Region (PBR) of the exon. Ten codons are under positive selection, nine of which are located in the mammal-defined PBR. We hypothesize that the tenth codon (residue 21) is an amphibian-specific PBR site that may be important in disease resistance. Trans-species and trans-generic polymorphisms are evident from exon-based genealogies, and co-phylogenetic analyses between intron, exon and mitochondrial based reconstructions reveal incongruent topologies, likely due to different locus histories. We developed two sets of barcoded adapters that reliably amplify a single and likely functional locus in all screened species using both 454 and Illumina based sequencing methods. These primers provide a resource for multiplexing and directly sequencing hundreds of samples in a single sequencing run, avoiding the labour and chimeric sequences associated with cloning, and enabling MHC population genetic analyses. Although the primers are currently limited to the 17 species we tested, these sequences and protocols provide a useful genetic resource and can serve as a starting point for future disease, adaptation and conservation studies across a range of anuran taxa.

Structural and functional diversity arising from intra- and inter-haplotype combinations of duplicated DQA and B loci within the ovine MHC

In sheep, the A and B loci encoding the α and β chains of the classical class II MHC molecules are DRA and DRB and DQA and DQB. Previous analyses described the duplication of the DQA and DQB genes. The majority of haplotypes include DQA1 and DQA2 loci, however, in a number of haplotypes, DQA1 appears absent and these haplotypes have been described as DQA1 null. In these haplotypes, the DQA2 locus is found in combination with a second locus which appeared more closely related to DQA2 than DQA1, hence the description of this locus as DQA2-like. Here we combine our previous analysis of the DQA transcripts with an analysis of the associated DQB transcripts in ten haplotypes from MHC homozygous animals. This allows the potential for surface expression of different haplotype combinations of DQA and B genes and the functional significance of DQA2-like and its predicted DQB partner to be determined. Atypical DQB transcripts (DQB2-like) were identified in haplotypes classified as DQA1-null and conserved DQB2-like orthologues were identified in other Bovidae indicating trans-species conservation of the allelic lineage. Functional combinations detected by co-transfection of DQ1, DQ2 and DQ2-like genes demonstrates the potential for a wide range of DQ molecules derived from both intra- and inter-haplotype as well as inter-locus combinations. We provide evidence that DQA2-like and B2-like genes form an evolutionary conserved pair which generates structurally distinct class II molecules that are likely to present a distinct range of peptides to CD4⁺ T cells.

NFM Cross-Reactivity to MOG Does Not Expand a Critical Threshold Level of High-Affinity T Cells Necessary for Onset of Demyelinating Disease

Of interest to the etiology of demyelinating autoimmune disease is the potential to aberrantly activate CD4⁺ T cells due to cross-recognition of multiple self-epitopes such as has been suggested for myelin oligodendrocyte glycoprotein epitope 35-55 (MOG_35-55) and neurofilament medium protein epitope 15-35 (NFM_15-35). NFM_15-35 is immunogenic in C57BL/6 mice but fails to induce demyelinating disease by polyclonal T cells despite having the same TCR contact residues as MOG_35-55, a known encephalitogenic Ag. Despite reported cross-reactivity with MOG-specific T cells, the polyclonal response to NFM_15-35 did not expand threshold numbers of MOG_38-49 tetramer-positive T cells. Furthermore, NFM lacked functional synergy with MOG to promote experimental autoimmune encephalomyelitis because NFM-deficient synonymous with knockout mice developed an identical disease course to wild-type mice after challenge with MOG_35-55 Single-cell analysis of encephalitogenic T cells using the peptide:MHC monomer-based two-dimensional micropipette adhesion frequency assay confirmed that NFM was not a critical Ag driving demyelinating disease because NFM_18-30-specific T cells in the CNS were predominantly reactive to MOG_38-49 The absence of NFM contribution to disease allowed mapping of the amino acids required for encephalitogenicity and expansion of high-affinity, MOG-specific T cells that defined the polyclonal response. Alterations of N-terminal residues outside of the NFM_15-35 core nonamer promoted expansion of high-affinity, MOG_38-49 tetramer-positive T cells and promoted consistent experimental autoimmune encephalomyelitis induction, unlike mice challenged with NFM_15-35 Although NFM_15-35 is immunogenic and cross-reactive with MOG at the polyclonal level, it fails to expand a threshold level of encephalitogenic, high-affinity MOG-specific T cells.

Cheetahs have a stronger constitutive innate immunity than leopards

As a textbook case for the importance of genetics in conservation, absence of genetic variability at the major histocompatibility complex (MHC) is thought to endanger species viability, since it is considered crucial for pathogen resistance. An alternative view of the immune system inspired by life history theory posits that a strong response should evolve in other components of the immune system if there is little variation in the MHC. In contrast to the leopard (Panthera pardus), the cheetah (Acinonyx jubatus) has a relatively low genetic variability at the MHC, yet free-ranging cheetahs are healthy. By comparing the functional competence of the humoral immune system of both species in sympatric populations in Namibia, we demonstrate that cheetahs have a higher constitutive innate but lower induced innate and adaptive immunity than leopards. We conclude (1) immunocompetence of cheetahs is higher than previously thought; (2) studying both innate and adaptive components of immune systems will enrich conservation science.

Targeted Next-Generation Sequencing for Human Leukocyte Antigen Typing in a Clinical Laboratory: Metrics of Relevance and Considerations for Its Successful Implementation

CONTEXT

- Numerous feasibility studies to type human leukocyte antigens (HLAs) by next-generation sequencing (NGS) have led to the development of vendor-supported kits for HLA typing by NGS. Some clinical laboratories have introduced HLA-NGS, and many are investigating the introduction. Standards from accrediting agencies form the regulatory framework for introducing this test into clinical laboratories.

OBJECTIVES

- To provide an assessment of metrics and considerations relevant to the successful implementation of clinical HLA-NGS typing, and to provide as a reference a validated HLA-NGS protocol used clinically since December 2013 at the Children's Hospital of Philadelphia (Philadelphia, Pennsylvania).

DATA SOURCES

- The HLA-NGS has been performed on 2532 samples. The initial 1046 and all homozygous samples were also typed by an alternate method. The HLA-NGS demonstrated 99.7% concordance with the alternate method. Ambiguous results were most common at the DPB1 locus because of a lack of phasing between exons 2 and 3 or the unsequenced exon 1 (533 of 2954 alleles; 18.04%) and the DRB1 locus because of not sequencing exon 1 (75 of 3972 alleles; 1.89%). No ambiguities were detected among the other loci. Except for 2 false homozygous samples, all homozygous samples (1891) demonstrated concordance with the alternate method. The article is organized to address the critical elements in the preanalytic, analytic, and postanalytic phases of introducing this assay into the clinical laboratory.

CONCLUSIONS

- The results demonstrate that HLA typing by NGS is a highly accurate, reproducible, efficient method that provides more-complete sequencing information for the length of the HLA gene and can be the single methodology for HLA typing in clinical immunogenetics laboratories.

HLA-DQ 1 alleles associated with Epstein-Barr virus (EBV) infectivity and EBV gp42 binding to cells

Epstein-Barr virus (EBV) infects B cells and ~95% of adults are infected. EBV glycoprotein gp42 is essential for entry of virus into B cells. EBV gp42 binds to the β1 chain of HLA-DQ, -DR, and -DP on B cells, and uses these molecules for infection. To investigate if certain HLA-DQ alleles are associated with EBV seronegativity, we recruited ~3,300 healthy adult blood donors, identified 106 EBV-seronegative individuals, and randomly selected a control group of EBV-seropositive donors from the donor pool. A larger than expected proportion of EBV-seronegative subjects were HLA-DQ β1 *04/*05 and *06/*06, and to a lesser extent, *02/*03, compared with the control group, while a larger than expected portion of EBV-seropositive persons were HLA-DQ β1 *02/*02. We examined the ability of EBV gp42 to bind to different HLA-DQ molecules using human and mouse cells stably expressing these alleles. EBV gp42 bound less effectively to cells expressing HLA-DQ β1 *04/*05, *06/*06, or *03/*03 than to cells expressing HLA-DQ β1 *02/*02. These data are consistent with our observations of increased EBV seronegativity with DQ β1 *04/*05 or *06/*06 alleles. These findings emphasize the importance of a single genetic locus (HLA-DQ β1) to influence infectivity with EBV.

Molecular characterization of MHC class IIB genes of sympatric Neotropical cichlids

BACKGROUND

The Major Histocompatibility Complex (MHC) is a key component of the adaptive immune system of all vertebrates and consists of the most polymorphic genes known to date. Due to this complexity, however, MHC remains to be characterized in many species including any Neotropical cichlid fish. Neotropical crater lake cichlids are ideal models to study evolutionary processes as they display one of the most convincing examples of sympatric and repeated parallel radiation events within and among isolated crater lakes.

RESULTS

Here, we characterized the genes of MHC class IIB chain of the Midas cichlid species complex (Amphilophus cf. citrinellus) including fish from five lakes in Nicaragua. We designed 19 new specific primers anchored in a stepwise fashion in order to detect all alleles present. We obtained 866 genomic DNA (gDNA) sequences from thirteen individuals and 756 additional sequences from complementary DNA (cDNA) of seven of those individuals. We identified 69 distinct alleles with up to 25 alleles per individual. We also found considerable intron length variation and mismatches of alleles detected in cDNA and gDNA suggesting that some loci have undergone pseudogenization. Lastly, we created a model of protein structure homology for each allele and identified their key structural components.

CONCLUSIONS

Overall, the Midas cichlid has one of the most diverse repertoires of MHC class IIB genes known, which could serve as a powerful tool to elucidate the process of divergent radiations, colonization and speciation in sympatry.

High-Throughput Sequencing of the Major Histocompatibility Complex following Targeted Sequence Capture

The Human Major Histocompatibility Complex (MHC) is a highly polymorphic region full of immunoregulatory genes. The MHC codes for the human leukocyte antigens (HLA), proteins that present on the cellular surface and that are involved in self-non-self recognition. For matching donors and recipients for organ and stem-cell transplants it is important to know an individual's HLA haplotype determinable in this region. Now, as next-generation sequencing (NGS) platforms mature and become more and more accepted as a standard method, NGS applications have spread from research laboratories to the clinic, where they provide valid genetic insights. Here, we describe a cost-effective microarray-based sequence capture, enrichment, and NGS sequencing approach to characterize MHC haplotypes. Using this approach, ~4 MB of MHC sequence for four DNA samples (donor, recipient and the parents of the recipient) were sequenced in parallel in one NGS instrument run. We complemented this approach using microarray-based genome-wide SNP analysis. Taken together, the use of recently developed tools and protocols for sequence capture and massively parallel sequencing allows for detailed MHC analysis and donor-recipient matching.

Polymorphism at expressed DQ and DR loci in five common equine MHC haplotypes

The polymorphism of major histocompatibility complex (MHC) class II DQ and DR genes in five common equine leukocyte antigen (ELA) haplotypes was determined through sequencing of mRNA transcripts isolated from lymphocytes of eight ELA homozygous horses. Ten expressed MHC class II genes were detected in horses of the ELA-A3 haplotype carried by the donor horses of the equine bacterial artificial chromosome (BAC) library and the reference genome sequence: four DR genes and six DQ genes. The other four ELA haplotypes contained at least eight expressed polymorphic MHC class II loci. Next generation sequencing (NGS) of genomic DNA of these four MHC haplotypes revealed stop codons in the DQA3 gene in the ELA-A2, ELA-A5, and ELA-A9 haplotypes. Few NGS reads were obtained for the other MHC class II genes that were not amplified in these horses. The amino acid sequences across haplotypes contained locus-specific residues, and the locus clusters produced by phylogenetic analysis were well supported. The MHC class II alleles within the five tested haplotypes were largely non-overlapping between haplotypes. The complement of equine MHC class II DQ and DR genes appears to be well conserved between haplotypes, in contrast to the recently described variation in class I gene loci between equine MHC haplotypes. The identification of allelic series of equine MHC class II loci will aid comparative studies of mammalian MHC conservation and evolution and may also help to interpret associations between the equine MHC class II region and diseases of the horse.

Gene conversion of the major histocompatibility complex class I Caja-G in common marmosets (Callithrix jacchus)

Currently, the amount of sequenced and classified MHC class I genes of the common marmoset is limited, in spite of the wide use of this species as an animal model for biomedical research. In this study, 480 clones of MHC class I G locus (Caja‐G) cDNA sequences were obtained from 21 common marmosets. Up to 10 different alleles were detected in each common marmoset, leading to the assumption that the Caja‐G loci duplicated in the marmoset genome. In the investigated population, four alleles occurred more often, giving evidence for higher immunological advantage of these alleles. In contrast to the human non‐classical MHC class I genes, Caja‐G shows high rates of polymorphism at the relevant peptide‐binding sites, despite its phylogenetic relationship to the non‐classical HLA‐G. Our results provide information for better understanding of the immunological properties of the common marmoset and confirm the theory of a gene conversion of the Caja‐G due to its detected plasticity and the absence of any known HLA‐A equivalent.

sNebula, a network-based algorithm to predict binding between human leukocyte antigens and peptides

Understanding the binding between human leukocyte antigens (HLAs) and peptides is important to understand the functioning of the immune system. Since it is time-consuming and costly to measure the binding between large numbers of HLAs and peptides, computational methods including machine learning models and network approaches have been developed to predict HLA-peptide binding. However, there are several limitations for the existing methods. We developed a network-based algorithm called sNebula to address these limitations. We curated qualitative Class I HLA-peptide binding data and demonstrated the prediction performance of sNebula on this dataset using leave-one-out cross-validation and five-fold cross-validations. This algorithm can predict not only peptides of different lengths and different types of HLAs, but also the peptides or HLAs that have no existing binding data. We believe sNebula is an effective method to predict HLA-peptide binding and thus improve our understanding of the immune system.

Homology modelling of frequent HLA class-II alleles: A perspective to improve prediction of HLA binding peptide and understand the HLA associated disease susceptibility

Human leukocyte antigen (HLA) plays significant role via the regulation of immune system and contribute in the progression and protection of many diseases. HLA molecules bind and present peptides to T- cell receptors which generate the immune response. HLA peptide interaction and molecular function of HLA molecule is the key to predict peptide binding and understanding its role in different diseases. The availability of accurate three dimensional (3D) structures is the initial step towards this direction. In the present work, homology modelling of important and frequent HLA-DRB1 alleles (07:01, 11:01 and 09:01) was done and acceptable models were generated. These modelled alleles were further refined and cross validated by using several methods including Ramachandran plot, Z-score, ERRAT analysis and root mean square deviation (RMSD) calculations. It is known that numbers of allelic variants are related to the susceptibility or protection of various infectious diseases. Difference in amino acid sequences and structures of alleles were also studied to understand the association of HLA with disease susceptibility and protection. Susceptible alleles showed more amino acid variations than protective alleles in three selected diseases caused by different pathogens. Amino acid variations at binding site were found to be more than other part of alleles. RMSD values were also higher at variable positions within binding site. Higher RMSD values indicate that mutations occurring at peptide binding site alter protein structure more than rest of the protein. Hence, these findings and modelled structures can be used to design HLA-DRB1 binding peptides to overcome low prediction accuracy of HLA class II binding peptides. Furthermore, it may help to understand the allele specific molecular mechanisms involved in susceptibility/resistance against pathogenic diseases.

The next generation of target capture technologies - large DNA fragment enrichment and sequencing determines regional genomic variation of high complexity

BACKGROUND

The ability to capture and sequence large contiguous DNA fragments represents a significant advancement towards the comprehensive characterization of complex genomic regions. While emerging sequencing platforms are capable of producing several kilobases-long reads, the fragment sizes generated by current DNA target enrichment technologies remain a limiting factor, producing DNA fragments generally shorter than 1 kbp. The DNA enrichment methodology described herein, Region-Specific Extraction (RSE), produces DNA segments in excess of 20 kbp in length. Coupling this enrichment method to appropriate sequencing platforms will significantly enhance the ability to generate complete and accurate sequence characterization of any genomic region without the need for reference-based assembly.

RESULTS

RSE is a long-range DNA target capture methodology that relies on the specific hybridization of short (20-25 base) oligonucleotide primers to selected sequence motifs within the DNA target region. These capture primers are then enzymatically extended on the 3'-end, incorporating biotinylated nucleotides into the DNA. Streptavidin-coated beads are subsequently used to pull-down the original, long DNA template molecules via the newly synthesized, biotinylated DNA that is bound to them. We demonstrate the accuracy, simplicity and utility of the RSE method by capturing and sequencing a 4 Mbp stretch of the major histocompatibility complex (MHC). Our results show an average depth of coverage of 164X for the entire MHC. This depth of coverage contributes significantly to a 99.94 % total coverage of the targeted region and to an accuracy that is over 99.99 %.

CONCLUSIONS

RSE represents a cost-effective target enrichment method capable of producing sequencing templates in excess of 20 kbp in length. The utility of our method has been proven to generate superior coverage across the MHC as compared to other commercially available methodologies, with the added advantage of producing longer sequencing templates amenable to DNA sequencing on recently developed platforms. Although our demonstration of the method does not utilize these DNA sequencing platforms directly, our results indicate that the capture of long DNA fragments produce superior coverage of the targeted region.

Excess of Deleterious Mutations around HLA Genes Reveals Evolutionary Cost of Balancing Selection

Deleterious mutations are expected to evolve under negative selection and are usually purged from the population. However, deleterious alleles segregate in the human population and some disease-associated variants are maintained at considerable frequencies. Here, we test the hypothesis that balancing selection may counteract purifying selection in neighboring regions and thus maintain deleterious variants at higher frequency than expected from their detrimental fitness effect. We first show in realistic simulations that balancing selection reduces the density of polymorphic sites surrounding a locus under balancing selection, but at the same time markedly increases the population frequency of the remaining variants, including even substantially deleterious alleles. To test the predictions of our simulations empirically, we then use whole-exome sequencing data from 6,500 human individuals and focus on the most established example for balancing selection in the human genome, the major histocompatibility complex (MHC). Our analysis shows an elevated frequency of putatively deleterious coding variants in nonhuman leukocyte antigen (non-HLA) genes localized in the MHC region. The mean frequency of these variants declined with physical distance from the classical HLA genes, indicating dependency on genetic linkage. These results reveal an indirect cost of the genetic diversity maintained by balancing selection, which has hitherto been perceived as mostly advantageous, and have implications both for the evolution of recombination and also for the epidemiology of various MHC-associated diseases.

HLA class I molecular variation and peptide-binding properties suggest a model of joint divergent asymmetric selection

The main function of HLA class I molecules is to present pathogen-derived peptides to cytotoxic T lymphocytes. This function is assumed to drive the maintenance of an extraordinary amount of polymorphism at each HLA locus, providing an immune advantage to heterozygote individuals capable to present larger repertories of peptides than homozygotes. This seems contradictory, however, with a reduced diversity at individual HLA loci exhibited by some isolated populations. This study shows that the level of functional diversity predicted for the two HLA-A and HLA-B genes considered simultaneously is similar (almost invariant) between 46 human populations, even when a reduced diversity exists at each locus. We thus propose that HLA-A and HLA-B evolved through a model of joint divergent asymmetric selection conferring all populations an equivalent immune potential. The distinct pattern observed for HLA-C is explained by its functional evolution towards killer cell immunoglobulin-like receptor (KIR) activity regulation rather than peptide presentation.

Deep sequencing of the MHC region in the Chinese population contributes to studies of complex disease

The human major histocompatibility complex (MHC) region has been shown to be associated with numerous diseases. However, it remains a challenge to pinpoint the causal variants for these associations because of the extreme complexity of the region. We thus sequenced the entire 5-Mb MHC region in 20,635 individuals of Han Chinese ancestry (10,689 controls and 9,946 patients with psoriasis) and constructed a Han-MHC database that includes both variants and HLA gene typing results of high accuracy. We further identified multiple independent new susceptibility loci in HLA-C, HLA-B, HLA-DPB1 and BTNL2 and an intergenic variant, rs118179173, associated with psoriasis and confirmed the well-established risk allele HLA-C*06:02. We anticipate that our Han-MHC reference panel built by deep sequencing of a large number of samples will serve as a useful tool for investigating the role of the MHC region in a variety of diseases and thus advance understanding of the pathogenesis of these disorders.

High-risk HLA alleles for severe acute graft-versus-host disease and mortality in unrelated donor bone marrow transplantation

HLA molecules play an important role for immunoreactivity in allogeneic hematopoietic stem cell transplantation. To elucidate the effect of specific HLA alleles on acute graft-versus-host disease, we conducted a retrospective analysis using 6967 Japanese patients transplanted with T-cell-replete marrow from an unrelated donor. Using unbiased searches of patient and donor HLA alleles, patient and/or donor HLA-B*51:01 (patient: HR, 1.37,P<0.001; donor: HR, 1.35,P<0.001) and patient HLA-C*14:02 (HR, 1.35,P<0.001) were significantly associated with an increased risk of severe acute graft-versus-host disease. The finding that donor HLA-C*14:02 was not associated with severe acute graft-versus-host disease prompted us to elucidate the relation of these high-risk HLA alleles with patient and donor HLA-C allele mismatches. In comparison to HLA-C allele match, patient mismatched HLA-C*14:02 showed the highest risk of severe acute graft-versus-host disease (HR, 3.61,P<0.001) and transplant-related mortality (HR, 2.53,P<0.001) among all patient mismatched HLA-C alleles. Although patient HLA-C*14:02 and donor HLA-C*15:02 mismatch was usually KIR2DL-ligand mismatch in the graft-versus-host direction, the risk of patient mismatched HLA-C*14:02 for severe acute graft-versus-host disease was obvious regardless of KIR2DL-ligand matching. The effect of patient and/or donor HLA-B*51:01 on acute graft-versus-host disease was attributed not only to strong linkage disequilibrium of HLA-C*14:02 and -B*51:01, but also to the effect of HLA-B*51:01 itself. With regard to clinical implications, patient mismatched HLA-C*14:02 proved to be a potent risk factor for severe acute graft-versus-host disease and mortality, and should be considered a non-permissive HLA-C mismatch in donor selection for unrelated donor hematopoietic stem cell transplantation.

Intrahaplotypic Variants Differentiate Complex Linkage Disequilibrium within Human MHC Haplotypes

Distinct regions of long-range genetic fixation in the human MHC region, known as conserved extended haplotypes (CEHs), possess unique genomic characteristics and are strongly associated with numerous diseases. While CEHs appear to be homogeneous by SNP analysis, the nature of fine variations within their genomic structure is unknown. Using multiple, MHC-homozygous cell lines, we demonstrate extensive sequence conservation in two common Asian MHC haplotypes: A33-B58-DR3 and A2-B46-DR9. However, characterization of phase-resolved MHC haplotypes revealed unique intra-CEH patterns of variation and uncovered 127 single nucleotide variants (SNVs) which are missing from public databases. We further show that the strong linkage disequilibrium structure within the human MHC that typically confounds precise identification of genetic features can be resolved using intra-CEH variants, as evidenced by rs3129063 and rs448489, which affect expression of ZFP57, a gene important in methylation and epigenetic regulation. This study demonstrates an improved strategy that can be used towards genetic dissection of diseases.

Genetics of leprosy: Expected-and unexpected-developments and perspectives

A solid body of evidence produced over decades of intense research supports the hypothesis that leprosy phenotypes are largely dependent on the genetic characteristics of the host. The early evidence of a major gene effect controlling susceptibility to leprosy came from studies of familial aggregation, twins, and complex segregation analysis. Later, linkage and association analysis, first applied to the investigation of candidate genes and chromosomal regions and more recently, to genome-wide scans, have revealed several HLA and non-HLA gene variants as risk factors for leprosy phenotypes such as disease per se, its clinical forms, and leprosy reactions. In addition, powerful, hypothesis-free strategies such as genome-wide association studies have led to an exciting, unexpected development: Leprosy susceptibility genes seem to be shared with Crohn's and Parkinson's disease. Today, a major challenge is to find the exact variants causing the biological effect underlying the genetic associations. New technologies, such as Next Generation Sequencing-that allows, for the first time, the cost- and time-effective sequencing of a complete human genome-hold the promise to reveal such variants; thus, strategies can be developed to study the functional impact of these variants in the context of infection, hopefully leading to the development of new targets for leprosy treatment and prevention.

A comprehensive mapping of the structure and gene organisation in the sheep MHC class I region

Background

The major histocompatibility complex (MHC) is a chromosomal region that regulates immune responsiveness in vertebrates. This region is one of the most important for disease resistance because it has been associated with resistance or susceptibility to a wide variety of diseases and because the MHC often accounts for more of the variance than other loci. Selective breeding for disease resistance is becoming increasingly common in livestock industries, and it is important to determine how this will influence MHC polymorphism and resistance to diseases that are not targeted for selection. However, in sheep the order and sequence of the protein coding genes is controversial. Yet this information is needed to determine precisely how the MHC influences resistance and susceptibility to disease.

Methods

CHORI bacterial artificial chromosomes (BACs) known to contain sequences from the sheep MHC class I region were sub-cloned, and the clones partially sequenced. The resulting sequences were analysed and re-assembled to identify gene content and organisation within each BAC. The low resolution MHC class I physical map was then compared to the cattle reference genome, the Chinese Merino sheep MHC map published by Gao, et al. (2010) and the recently available sheep reference genome.

Results

Immune related class I genes are clustered into 3 blocks; beta, kappa and a novel block not previously identified in other organisms. The revised map is more similar to Bovidae maps than the previous sheep maps and also includes several genes previously not annotated in the Chinese Merino BAC assembly and others not currently annotated in the sheep reference chromosome 20. In particular, the organisation of nonclassical MHC class I genes is similar to that present in the cattle MHC. Sequence analysis and prediction of amino acid sequences of MHC class I classical and nonclassical genes was performed and it was observed that the map contained one classical and eight nonclassical genes together with three possible pseudogenes.

Conclusions

The comprehensive physical map of the sheep MHC class I region enhances our understanding of the genetic architecture of the class I MHC region in sheep and will facilitate future studies of MHC function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1992-4) contains supplementary material, which is available to authorized users.

Machine Learning Methods for Predicting HLA-Peptide Binding Activity

As major histocompatibility complexes in humans, the human leukocyte antigens (HLAs) have important functions to present antigen peptides onto T-cell receptors for immunological recognition and responses. Interpreting and predicting HLA–peptide binding are important to study T-cell epitopes, immune reactions, and the mechanisms of adverse drug reactions. We review different types of machine learning methods and tools that have been used for HLA–peptide binding prediction. We also summarize the descriptors based on which the HLA–peptide binding prediction models have been constructed and discuss the limitation and challenges of the current methods. Lastly, we give a future perspective on the HLA–peptide binding prediction method based on network analysis.

ATPase4A Autoreactivity and Its Association With Autoimmune Phenotypes in the Type 1 Diabetes Genetics Consortium Study

Autoantibodies targeting the H+/K+-ATPase proton pump of the gastric parietal cell (parietal cell antibodies [PCA]) are diagnostic of atrophic body gastritis (ABG) leading to pernicious anemia (PA). PCA, ABG, and PA occur in increased frequency in patients with type 1 diabetes and their relatives and are considered "minor" components of forms of autoimmune polyglandular syndrome (APS). A customized radioimmunoprecipitation assay was applied to 6,749 samples from the Type 1 Diabetes Genetics Consortium to measure ATP4A autoreactivity. Autoantibody prevalence was correlated with variants in HLA class II, PTPN22, and CTLA4 genes. With an ATP4A radioimmunoprecipitation assay, PCA were detected in sera from 20.9% of affected individuals. PCA prevalence increased with age and was greater in females (25.3%) than males (16.5%) and among Hispanics (36.3%) and blacks (26.2%) compared with non-Hispanic whites (20.8%) and Asians (16.7%). PCA and other organ-specific autoantibodies GAD65, IA-2, thyroid peroxidase (TPO), 21-hydroxylase (21-OH), and transglutaminase (TG) clustered within families with heritability estimates from 71 to 95%. PCA clustered with TPO, 21-OH, and persistent GAD65 autoantibodies but not with celiac (TG) or IA-2 autoantibodies. PCA-positive subjects showed an increased frequency of DRB1*0404, DPB1*0201, and PTPN22 R620W (rs2476601-T) and a decreased frequency of DRB1*0101, DPB1*0301, and CTLA4 CT60 (rs3087243-T). Genetic variants accounted for 4-5% of the heritable risk for PCA. The same alleles were associated with other autoantibody phenotypes in a consistent pattern. Whereas most of the heritable risk for PCA and other antibodies reflects genetic effects that are tissue specific, parietal cell autoimmunity is a major pathogenetic contributor in APS2.

Factors influencing reproduction in captive-bred cynomolgus monkeys (Macaca fascicularis) from Mauritius

The cynomolgus monkey is widely used in reproductive research. However, the effects on their reproductive parameters of infant and maternal factors such as birth order, sex of infants, twin births, maternal age and lactation status have not been fully examined. The aim of this retrospective study was to determine how such infant and maternal factors impact on infant birth weight, birth viability, neonatal loss and retained placenta in cynomolgus monkeys. The study was based on birth data from a cohort of 789 females over an eight-year period. Consistent with reports made in other macaque species, female offspring had lower birth weight compared with males. Birth weights of firstborn infants were lower compared with birth weights of higher birth order infants. Results from the logistic regression analysis showed that the risk of non-viable births was increased by advancing maternal age and retained placenta. As in other non-human primates, maternal age had predictive value for non-viable births in cynomolgus monkeys. The risk of neonatal loss decreased with advancing maternal age but was not affected by birth order. Firstborn offspring did not have an increased risk for neonatal loss, possibly from the practice of retaining mothers in their natal groups, which improved maternal skills in primiparous females. However, infant low birth weight and non-lactating females increased the risk of neonatal loss, and the delivery of low birth weight infants was associated with retained placenta. The results from this study can be useful for scientists conducting reproductive studies and for colony managers in maximizing fertility and infant survival of cynomolgus monkeys.