The spectrum of HLA-DQ and HLA-DR alleles, 2006: a listing correlating sequence and structure with function

HLA peptide-binding pocket diversity modulates immunological complications after cord blood transplant in acute leukaemia

Pocket motifs and their amino acid positions of HLA molecules are known to govern antigen presentation to effector cells. Our objective was to analyse their influence on the risk of graft-versus-host disease (GVHD) and relapse after umbilical cord blood transplant (UCBT). The transplant characteristics of 849 patients with acute leukaemia were obtained from the Eurocord/EBMT database. Higher acute (a) GVHD was associated with homozygosity of UCB HLA-C amino acid positions 77 and 80 (NN/KK) (p = 0.008). Severe aGVHD was associated with HLA-A pocket B YSAVMENVHY motif (p = 0.002) and NN and RR genotypes of the HLA-C amino acid positions 77 and 156 (p = 0.006 and p = 0.002). Such risk was also increased in case of recipient and UCB mismatches in P4 (p < 0.0001) and P9 (p = 0.003) pockets of HLA-DQB1 alleles. For chronic GVHD, the pocket B YYAVMEISNY motif of the HLA-B*15:01 allele and the absence of mismatch between recipient and UCB in the P6 pocket of HLA-DRB1 were associated with a lower risk (p = 0.0007 and p = 0.0004). In relapse, both UCB pocket B YFAVMENVHY belonging to HLA-A*32:01 and recipient pocket B YDSVGENYQY motif of the HLA-C*07:01 allele were associated with higher risk (p = 0.0026 and p = 0.015). We provide clues on HLA-mediated cellular interactions and their role in the development of GVHD and relapse.

Genetic determinants of IgG antibody response to COVID-19 vaccination

Human humoral immune responses to SARS-CoV-2 vaccines exhibit substantial inter-individual variability and have been linked to vaccine efficacy. To elucidate the underlying mechanism behind this variability, we conducted a genome-wide association study (GWAS) on the anti-spike IgG serostatus of UK Biobank participants who were previously uninfected by SARS-CoV-2 and had received either the first dose (n = 54,066) or the second dose (n = 46,232) of COVID-19 vaccines. Our analysis revealed significant genome-wide associations between the IgG antibody serostatus following the initial vaccine and human leukocyte antigen (HLA) class II alleles. Specifically, the HLA-DRB1∗13:02 allele (MAF = 4.0%, OR = 0.75, p = 2.34e-16) demonstrated the most statistically significant protective effect against IgG seronegativity. This protective effect was driven by an alteration from arginine (Arg) to glutamic acid (Glu) at position 71 on HLA-DRβ1 (p = 1.88e-25), leading to a change in the electrostatic potential of pocket 4 of the peptide binding groove. Notably, the impact of HLA alleles on IgG responses was cell type specific, and we observed a shared genetic predisposition between IgG status and susceptibility/severity of COVID-19. These results were replicated within independent cohorts where IgG serostatus was assayed by two different antibody serology tests. Our findings provide insights into the biological mechanism underlying individual variation in responses to COVID-19 vaccines and highlight the need to consider the influence of constitutive genetics when designing vaccination strategies for optimizing protection and control of infectious disease across diverse populations.

Association Study Identified HLA-DQA1 as a Novel Genetic Risk of Systemic Lupus Erythematosus-Associated Pulmonary Arterial Hypertension

OBJECTIVE

Pulmonary arterial hypertension (PAH) is a severe complication of systemic lupus erythematosus (SLE). However, the genetic signatures of SLE-associated PAH have not been well studied. We aimed to identify genetic variants implicated in SLE-associated PAH susceptibility within the major histocompatibility complex (MHC) region and assess the contribution to clinical outcomes.

METHODS

A total of 172 patients with SLE-associated PAH confirmed by right heart catheterization, 1,303 patients with SLE without PAH, and 9,906 healthy controls were included. Deep sequencing of the MHC region was performed to identify alleles, single-nucleotide polymorphisms, and amino acids. We compared patients with SLE-associated PAH with patients with SLE without PAH and healthy controls. Clinical association study was conducted to explore the contribution to phenotypes.

RESULTS

A total of 19,881 genetic variants were identified within the MHC region. HLA-DQA1*03:02 was identified as a novel genetic variant associated with SLE-associated PAH in the discovery cohort (P = 5.68 × 10-12 ) and authenticated in an independent replication cohort (P = 1.30 × 10-9 ). The strongest associated amino acid position was mapped to HLA-DQα1 in the region affecting MHC/peptide-CD4+ T cell receptor affinity and antigen binding. Clinical association study demonstrated that patients with SLE-associated PAH with HLA-DQA1*03:02 had significantly lower rates of target role achievement (P = 0.005) and survival (P = 0.04).

CONCLUSION

This study, based on the largest cohort of SLE-associated PAH, is the first to investigate how MHC region genetic variants contribute to SLE-associated PAH susceptibility. HLA-DQA1*03:02 is a novel genetic risk factor and a prognostic factor in SLE-associated PAH. Patients with SLE with this allele require regular monitoring and careful follow-up for early diagnosis and interventions for potential PAH.

Cross-population genetic analysis revealed genetic variation and selection in the Ovar-DRB1 gene of Indian sheep breeds

In sheep, MHC variability is studied widely to explore disease association. The aim of the current study was to explore the genetic diversity of Ovar-DRB diversity across sheep breeds of India. Here, Ovar-DRB1 locus was studied across 20 sheep breeds. DRB1 was amplified (301 bp) and sequenced using a PCR-sequence-based typing approach. Results revealed a high degree of heterozygosity across breeds (mean: 73.99%). Overall mean distance for DRB1 was highest in Sangamneri (0.18) and lowest in Madgyal sheep (0.10). There was a higher rate of transition, across breeds. Further, 39 alleles were isolated in different breeds, out of which 10 were new. To allow easy access and use of the immune-polymorphic database, an online database management system was launched (http://www.mhcdbms.in/). Nucleotide content across breeds for the DRB1 region revealed the richness of GC content (59.26%). Wu-Kabat index revealed vast genetic variation across peptide binding sites (PBS) of DRB1. Residues 6, 66, 69, 52, and 81, were polymorphic showing utility for antigen presentation. All breeds were under positive selection for DRB1 locus (dN > dS). Study revealed the importance of DRB locus diversity for beta chain specifically at PBS across sheep breeds of the Indian subcontinent and presented evidence of positive selection for DRB owing to its evolutionary significance.

Important denominator between autoimmune comorbidities: a review of class II HLA, autoimmune disease, and the gut

Human leukocyte antigen (HLA) genes are associated with more diseases than any other region of the genome. Highly polymorphic HLA genes produce variable haplotypes that are specifically correlated with pathogenically different autoimmunities. Despite differing etiologies, however, many autoimmune disorders share the same risk-associated HLA haplotypes often resulting in comorbidity. This shared risk remains an unanswered question in the field. Yet, several groups have revealed links between gut microbial community composition and autoimmune diseases. Autoimmunity is frequently associated with dysbiosis, resulting in loss of barrier function and permeability of tight junctions, which increases HLA class II expression levels and thus further influences the composition of the gut microbiome. However, autoimmune-risk-associated HLA haplotypes are connected to gut dysbiosis long before autoimmunity even begins. This review evaluates current research on the HLA-microbiome-autoimmunity triplex and proposes that pre-autoimmune bacterial dysbiosis in the gut is an important determinant between autoimmune comorbidities with systemic inflammation as a common denominator.

High functional allelic diversity and copy number in both MHC classes in the common buzzard

BACKGROUND

The major histocompatibility complex (MHC), which encodes molecules that recognize various pathogens and parasites and initiates the adaptive immune response in vertebrates, is renowned for its exceptional polymorphism and is a model of adaptive gene evolution. In birds, the number of MHC genes and sequence diversity varies greatly among taxa, believed due to evolutionary history and differential selection pressures. Earlier characterization studies and recent comparative studies suggest that non-passerine species have relatively few MHC gene copies compared to passerines. Additionally, comparative studies that have looked at partial MHC sequences have speculated that non-passerines have opposite patterns of selection on MHC class I (MHC-I) and class II (MHC-II) loci than passerines: namely, greater sequence diversity and signals of selection on MHC-II than MHC-I. However, new sequencing technology is revealing much greater MHC variation than previously expected while also facilitating full sequence variant detection directly from genomic data. Our study aims to take advantage of high-throughput sequencing methods to fully characterize both classes and domains of MHC of a non-passerine bird of prey, the common buzzard (Buteo buteo), to test predictions of MHC variation and differential selection on MHC classes.

RESULTS

Using genetic, genomic, and transcriptomic high-throughput sequencing data, we established common buzzards have at least three loci that produce functional alleles at both MHC classes. In total, we characterize 91 alleles from 113 common buzzard chicks for MHC-I exon 3 and 41 alleles from 125 chicks for MHC-IIB exon 2. Among these alleles, we found greater sequence polymorphism and stronger diversifying selection at MHC-IIB exon 2 than MHC-I exon 3, suggesting differential selection pressures on MHC classes. However, upon further investigation of the entire peptide-binding groove by including genomic data from MHC-I exon 2 and MHC-IIA exon 2, this turned out to be false. MHC-I exon 2 was as polymorphic as MHC-IIB exon 2 and MHC-IIA exon 2 was essentially invariant. Thus, comparisons between MHC-I and MHC-II that included both domains of the peptide-binding groove showed no differences in polymorphism nor diversifying selection between the classes. Nevertheless, selection analysis indicates balancing selection has been acting on common buzzard MHC and phylogenetic inference revealed that trans-species polymorphism is present between common buzzards and species separated for over 33 million years for class I and class II.

CONCLUSIONS

We characterize and confirm the functionality of unexpectedly high copy number and allelic diversity in both MHC classes of a bird of prey. While balancing selection is acting on both classes, there is no evidence of differential selection pressure on MHC classes in common buzzards and this result may hold more generally once more data for understudied MHC exons becomes available.

A genome-wide association study for allergen component sensitizations identifies allergen component-specific and allergen protein group-specific associations

Background

Allergic diseases are some of the most common diseases worldwide. Genome-wide association studies (GWASs) have been conducted to elucidate the genetic factors of allergic diseases. However, no GWASs for allergen component sensitization have been performed.

Objective

We sought to detect genetic variants associated with differences in immune responsiveness against allergen components.

Methods

The participants of the present study were recruited from the Tokyo Children's Health, Illness, and Development study, and allergen component-specific IgE level at age 9 years was measured by means of allergen microarray immunoassays. We performed GWASs for allergen component sensitization against each allergen (single allergen component sensitization, number of allergen components analyzed, n = 31), as well as against allergen protein families (allergen protein group sensitization, number of protein groups analyzed, n = 16).

Results

We performed GWAS on 564 participants of the Tokyo Children's Health, Illness, and Development study and found associations between Amb a 1 sensitization and the immunoglobulin heavy-chain variable gene on chromosome 14 and between Phl p 1 sensitization and the HLA class II region on chromosome 6 (P < 5.0 × 10-8). A GWAS-significant association was also observed between the HLA class II region and profilin sensitization (P < 5.0 × 10-8).

Conclusions

Our data provide the first demonstration of genetic risk for allergen component sensitization and show that this genetic risk is related to immune response genes including immunoglobulin heavy-chain variable gene and HLA.

Amino acid polymorphisms in human histocompatibility leukocyte antigen class II and proinsulin epitope have impacts on type 1 diabetes mellitus induced by immune-checkpoint inhibitors

Introduction

Immune-checkpoint inhibitors are effective in various advanced cancers. Type 1 diabetes mellitus induced by them (ICI-T1DM) is a serious complication requiring prompt insulin treatment, but the immunological mechanism behind it is unclear.

Methods

We examined amino acid polymorphisms in human histocompatibility leukocyte antigen (HLA) molecules and investigated proinsulin epitope binding affinities to HLA molecules.

Results and Discussion

Twelve patients with ICI-T1DM and 35 patients in a control group without ICI-T1DM were enrolled in the study. Allele and haplotype frequencies of HLA-DRB1*04:05, DQB1*04:01, and most importantly DPB1*05:01 were significantly increased in patients with ICI-T1DM. In addition, novel amino acid polymorphisms in HLA-DR (4 polymorphisms), in DQ (12 polymorphisms), and in DP molecules (9 polymorphisms) were identified. These amino acid polymorphisms might be associated with the development of ICI-T1DM. Moreover, novel human proinsulin epitope clusters in insulin A and B chains were discovered in silico and in vitro peptide binding assays to HLA-DP5. In conclusion, significant amino acid polymorphisms in HLA-class II molecules, and conformational alterations in the peptide-binding groove of the HLA-DP molecules were considered likely to influence the immunogenicity of proinsulin epitopes in ICI-T1DM. These amino acid polymorphisms and HLA-DP5 may be predictive genetic factors for ICI-T1DM.

West Nile virus-associated HLA-DRB1 alleles in the Greek population: A structural perspective

The HLA system plays a significant role via the regulation of the immune system and contributes to the progression and protection of many diseases. In our previous study, several HLA-DRB1 alleles were found to have a susceptible or protective role toward infection and neuroinvasion of West Nile Virus (WNV) in the Greek population. As expected, the majority of polymorphic positions are located in the peptide-binding region of the molecule. In the present work, the structure of these alleles was studied in silico, to examine the effect of polymorphism on the conformation of DRB1 proteins, with the aspect of WNV association. More specifically, molecular dynamics simulations were used for structural prediction of 23 available alleles. These modeled alleles were evaluated using root-mean-square deviation (RMSD) and root-mean-square fluctuation analysis. Low RMSD values indicate that different alleles have similar structures. Furthermore, low fluctuation was observed in the peptide-binding region between alleles with the higher and the lowest RMSD values. These findings indicate that probably variable residues do not affect the behavior of DRB1 alleles in WNV disease, by causing structural differences between them.

Development of a high-resolution mass-spectrometry-based method and software for human leukocyte antigen typing

Introduction

The human leukocyte antigen (HLA) system plays a critical role in the human immune system and is strongly associated with immune recognition and rejection in organ transplantation. HLA typing method has been extensively studied to increase the success rates of clinical organ transplantation. However, while polymerase chain reaction sequence-based typing (PCR-SBT) remains the gold standard, cis/trans ambiguity and nucleotide sequencing signal overlay during heterozygous typing present a problem. The high cost and low processing speed of Next Generation Sequencing (NGS) also render this approach inadequate for HLA typing.

Methods and materials

To address these limitations of the current HLA typing methods, we developed a novel typing technology based on nucleic acid mass spectrometry (MS) of HLA. Our method takes advantage of the high-resolution mass analysis function of MS and HLAMSTTs (HLA MS Typing Tags, some short fragment PCR amplification target products) with precise primer combinations.

Results

We correctly typed HLA by measuring the molecular weights of HLAMSTTs with single nucleotide polymorphisms (SNPs). In addition, we developed a supporting HLA MS typing software to design PCR primers, construct the MS database, and select the best-matching HLA typing results. With this new method, we typed 16 HLA-DQA1 samples, including 6 homozygotes and 10 heterozygotes. The MS typing results were validated by PCR-SBT.

Discussion

The MS HLA typing method is rapid, efficient, accurate, and readily applicable to typing of homozygous and heterozygous samples.

Association of HLA-DQ Heterodimer Residues -18β and β57 With Progression From Islet Autoimmunity to Diabetes in the Diabetes Prevention Trial-Type 1

OBJECTIVE

The purpose was to test the hypothesis that the HLA-DQαβ heterodimer structure is related to the progression of islet autoimmunity from asymptomatic to symptomatic type 1 diabetes (T1D).

RESEARCH DESIGN AND METHODS

Next-generation targeted sequencing was used to genotype HLA-DQA1-B1 class II genes in 670 subjects in the Diabetes Prevention Trial-Type 1 (DPT-1). Coding sequences were translated into DQ α- and β-chain amino acid residues and used in hierarchically organized haplotype (HOH) association analysis to identify motifs associated with diabetes onset.

RESULTS

The opposite diabetes risks were confirmed for HLA DQA1*03:01-B1*03:02 (hazard ratio [HR] 1.36; P = 2.01 ∗ 10-3) and DQA1*03:03-B1*03:01 (HR 0.62; P = 0.037). The HOH analysis uncovered residue -18β in the signal peptide and β57 in the β-chain to form six motifs. DQ*VA was associated with faster (HR 1.49; P = 6.36 ∗ 10-4) and DQ*AD with slower (HR 0.64; P = 0.020) progression to diabetes onset. VA/VA, representing DQA1*03:01-B1*03:02 (DQ8/8), had a greater HR of 1.98 (P = 2.80 ∗ 10-3). The DQ*VA motif was associated with both islet cell antibodies (P = 0.023) and insulin autoantibodies (IAAs) (P = 3.34 ∗ 10-3), while the DQ*AD motif was associated with a decreased IAA frequency (P = 0.015). Subjects with DQ*VA and DQ*AD experienced, respectively, increasing and decreasing trends of HbA1c levels throughout the follow-up.

CONCLUSIONS

HLA-DQ structural motifs appear to modulate progression from islet autoimmunity to diabetes among at-risk relatives with islet autoantibodies. Residue -18β within the signal peptide may be related to levels of protein synthesis and β57 to stability of the peptide-DQab trimolecular complex.

Allelic diversity at BoLA DRB3 locus and association with predisposition to clinical mastitis in indicus and crossbred cattle

Bovine lymphocyte antigen (BoLA) DRB3 locus in healthy and mastitis affected cattle has been genotyped by a polymerase chain reaction and restriction fragment length polymorphisms (PCR-RLFP) using RsaI restriction enzyme, followed by sequencing. In 130 farm animals, 25 BoLA DRB3 alleles have been detected by PCR-RFLP. Three distinct allelic patterns significantly associated with mastitis in Karan Fries crossbred and Sahiwal indicus cattle have been identified, whereas, four other allelic patterns were significantly high in frequency among healthy animals. Sequencing of RFLP genotypes revealed 25 and 47 alleles among healthy Sahiwal and Karan Fries, respectively, while 17 and 38 patterns observed in mastitis affected Sahiwal and Karan Fries animals, respectively. From Tajima's D-test of neutrality, it was concluded that alleles associated with mastitis were expanding in the population, whereas those of healthy were under contraction. Phylogenetic analysis carried out to delineate the evolutionary relationship of the farm and field animals at DRB3 locus, differentiating allelic patterns into six different clusters. Among the phylogenetic lineages, five patterns DRB3*028:01, DRB3*011:03, DRB3*031:01, DRB3*001:01 and DRB3*043:01, were previously reported, whereas one novel allelic variant was observed in indicus and crossbred cattle. This information will help in further exploring the association between BoLA-DRB3 genetic diversity and disease resistance in distinct cattle breeds, important in designing breeding strategies for increasing the distribution of favorable alleles.

Next Generation Sequencing Identifies the HLA-DQA1*03:03 Allele in the Type 1 Diabetes Risk-Associated HLA-DQ8 Serotype

The highest genetic type 1 diabetes risk is conferred by HLA class II haplotypes defined by alleles at the HLA-DR and -DQ loci. The combination of HLA-DQA1*03:01 and DQB1*03:02 alleles (summarized as ‘HLA-DQ8′) is reported to be among the two most prevalent HLA class II haplotypes in Caucasian type 1 diabetes patients. This classification is based on conventional genotyping of exon 2 of the DQ gene locus and excludes exon 3. In this study, HLA genotyping on the type 1 diabetes susceptibility loci HLA-DRB1, DQA1 and DQB1 was performed using a high-resolution next generation sequencing method. In addition to the routinely examined exon 2, exon 3 was also sequenced. Samples from 229 children with type 1 diabetes were included and compared to a cohort of 9,786 controls. In addition to previously described HLA-DQ haplotypes in type 1 diabetes patients, we found that as well as HLA-DQA1*03:01,HLA-DQA1*03:03 also contributed to HLA-DQ8. HLA-DQA1*03:03 differs from HLA-DQA1*03:01 by one nucleotide substitution in exon 3 at position 160, leading to a single amino acid replacement. DRB1*04:05 was exclusively associated with DQA1*03:03 whereas the DRB1*04:01 haplotype comprised either DQA1*03:01 or DQA1*03:03. Significantly increased type 1 diabetes risk was confirmed for all these haplotypes with only minor differences between DQA1*03:01 and DQA1*03:03 alleles. This study identified the HLA-DQA1*03:03 allele as an addition to the already known type 1 diabetes risk haplotypes, and can contribute to more precise HLA genotyping approaches.

Low neutral and immunogenetic diversity in northern fringe populations of the green toad Bufotes viridis: implications for conservation

Genetic variation is often lower at high latitudes, which may compromise the adaptability and hence survival of organisms. Here we show that genetic variability is negatively correlated with northern latitude in European green toads (Bufotes viridis). The result holds true for both putatively neutral microsatellite variation and supposedly adaptive MHC Class IIB variation. In particular, our findings have bearing on the conservation status of this species in Sweden, on the northern limit of its distribution where local populations are small and fragmented. These genetically impoverished populations are closely related to other populations found around the Baltic Sea basin. The low neutral and adaptive variation in these fringe populations compared to population at central ranges confirms a pattern shared across all other amphibians so far studied. In Sweden, the situation of green toads is of concern as the remaining populations may not have the evolutionary potential to cope with present and future environmental challenges.

Associations between specific IgE sensitization to 26 respiratory allergen molecules and HLA class II alleles in the EGEA cohort

BACKGROUND

Allergy, the most frequent immune disorder affecting 30% of the world's population, is the consequence of immunoglobin E (IgE) sensitization to allergens. Among the genetic factors suspected to be involved in allergy, the HLA class-II genomic region is a strong candidate.

OBJECTIVE

To assess the association between HLA class-II alleles and specific IgE (sIgE) sensitization to a large number of respiratory allergen molecules.

METHODS

The analysis relied on 927 participants of the EGEA cohort, including 497 asthmatics. The study focuses on 26 aeroallergens recognized by sIgE in at least 5% of the study population (determined with the MEDALL chip with sIgE ≥ 0.3 ISU) and 23 imputed HLA class-II alleles. For each sIgE sensitization and HLA class-II allele, we fitted a logistic regression model accounting for familial dependence and adjusted for gender, age, and genetic principal components. p-values were corrected for multiple comparisons (False Discovery Rate).

RESULTS

Most of the 19 statistically significant associations observed regard pollen allergens (mugwort Art v 1, olive tree Ole e 1, timothy grass Phl p 2, Phl p 5 and plantain Pla l 1), three were mold allergen (Alternaria Alt a 1), and a single one regards house dust mite allergen (Der p 7). No association was observed with pet allergens. The strongest associations were found with mugwort Art v 1 (OR = 5.42 (95%CI, 3.30; 8.88), 4.14 (2.65; 6.47), 3.16 (1.88; 5.31) with DQB1*05:01, DQA1*01:01 and DRB1*01:01, respectively).

CONCLUSION

Our results support the important role of HLA class-II alleles as immune response genes predisposing their carriers for sensitization to various major pollen allergens.

HLA class II peptide-binding-region analysis reveals funneling of polymorphism in action

BACKGROUND

HLA-class II proteins hold important roles in key physiological processes. The purpose of this study was to compile all class II alleles reported in human population and investigate patterns in pocket variants and their combinations, focusing on the peptide-binding region (PBR).

METHODS

For this purpose, all protein sequences of DPA1, DQA1, DPB1, DQB1 and DRB1 were selected and filtered, in order to have full PBR sequences. Proportional representation was used for pocket variants while population data were also used.

RESULTS

All pocket variants and PBR sequences were retrieved and analyzed based on the preference of amino acids and their properties in all pocket positions. The observed number of pocket variants combinations was much lower than the possible inferred, suggesting that PBR formation is under strict funneling. Also, although class II proteins are very polymorphic, in the majority of the reported alleles in all populations, a significantly less polymorphic pocket core was found.

CONCLUSIONS

Pocket variability of five HLA class II proteins was studied revealing favorable properties of each protein. The actual PBR sequences of HLA class II proteins appear to be governed by restrictions that lead to the establishment of only a fraction of the possible combinations and the polymorphism recorded is the result of intense funneling based on function.

Next-Generation HLA Sequence Analysis Uncovers Shared Risk Alleles Between Clinically Distinct Forms of Childhood Uveitis

Purpose

Classical alleles of the human leukocyte antigen (HLA) complex have been linked to specific entities of pediatric noninfectious uveitis, yet genetic predisposition encoded by the HLA super-locus across the patient population remains understudied.

Methods

We performed next-generation full-length sequencing of HLA-A, HLA-B, HLA-C, HLA-DPB1, HLA-DQB1, and HLA-DRB1 in 280 cases. Dense genotype data from 499 Dutch controls from Genome of the Netherlands were imputed using an HLA-specific reference panel (n = 5225 samples from European ancestry). Cases and controls were compared using logistic regression models adjusting for sex.

Results

In total, 179 common and rare alleles were detected. Considering all cases and controls, HLA-DQB1*04:02 and HLA-DRB1*08:01 were identified as the principal HLA association, which was mainly driven by 92 cases with juvenile idiopathic arthritis-associated uveitis (JIA-U). The HLA-DQB1*04:02-HLA-DRB1*08:01 haplotype was also the primary association for the phenotypically similar idiopathic chronic anterior uveitis without arthritis (CAU). Also, HLA-DQB1*05:03 was an independent risk allele for CAU, but not in JIA-U. Analysis of 185 cases with other forms of uveitis revealed HLA-wide associations (P < 2.79 × 10⁻⁴) for HLA-DRB1*01:02, HLA-DRB1*04:03, and HLA-DQB1*05:03, which could be primarily attributed to cases with panuveitis. Finally, amino acid substitution modeling revealed that aspartic acid at position 57 that distinguishes the risk allele HLA-DQB1*05:03 (for CAU and panuveitis) from nonrisk alleles, significantly increased the binding capacity of naturally presented ligands to HLA-DQ.

Conclusions

These results uncovered novel shared HLA associations among clinically distinct phenotypes of pediatric uveitis and highlight genetic predisposition affecting the antigen presentation pathway.

The KAG motif of HLA-DRB1 (β71, β74, β86) predicts seroconversion and development of type 1 diabetes

BACKGROUND

HLA-DR4, a common antigen of HLA-DRB1, has multiple subtypes that are strongly associated with risk of type 1 diabetes (T1D); however, some are risk neutral or resistant. The pathobiological mechanism of HLA-DR4 subtypes remains to be elucidated.

METHODS

We used a population-based case-control study of T1D (962 patients and 636 controls) to decipher genetic associations of HLA-DR4 subtypes and specific residues with susceptibility to T1D. Using a birth cohort of 7865 children with periodically measured islet autoantibodies (GADA, IAA or IA-2A), we proposed to validate discovered genetic associations with a totally different study design and time-to-seroconversions prior to clinical onset of T1D. A novel analytic strategy hierarchically organized the HLA-DRB1 alleles by sequence similarity and identified critical amino acid residues by minimizing local genomic architecture and higher-order interactions.

FINDINGS

Three amino acid residues of HLA-DRB1 (β71, β74, β86) were found to be predictive of T1D risk in the population-based study. The "KAG" motif, corresponding to HLA-DRB1×04:01, was most strongly associated with T1D risk ([O]dds [R]atio=3.64, p = 3.19 × 10-64). Three less frequent motifs ("EAV", OR = 2.55, p = 0.025; "RAG", OR = 1.93, p = 0.043; and "RAV", OR = 1.56, p = 0.003) were associated with T1D risk, while two motifs ("REG" and "REV") were equally protective (OR = 0.11, p = 4.23 × 10-4). In an independent birth cohort of HLA-DR3 and HLA-DR4 subjects, those having the "KAG" motif had increased risk for time-to-seroconversion (Hazard Ratio = 1.74, p = 6.51 × 10-14) after adjusting potential confounders.

INTERPRETATIONS

DNA sequence variation in HLA-DRB1 at positions β71, β74, and β86 are non-conservative (β74 A→E, β71 E vs K vs R and β86 G vs V). They result in substantial differences in peptide antigen anchor pocket preferences at p1, p4 and potentially neighboring regions such as pocket p7. Differential peptide antigen binding is likely to be affected. These sequence substitutions may account for most of the HLA-DR4 contribution to T1D risk as illustrated in two HLA-peptide model complexes of the T1D autoantigens preproinsulin and GAD65.

FUNDING

National Institute of Diabetes and Digestive and Kidney Diseases and the Swedish Child Diabetes Foundation and the Swedish Research Council.

Positive selection on the MHC class II DLA-DQA1 gene in golden jackals (Canis aureus) from their recent expansion range in Europe and its effect on their body mass index

Background

In Europe, golden jackals (Canis aureus) have been expanding their range out of the southern and southeastern Balkans towards central Europe continually since the 1960s. Here, we investigated the level of functional diversity at the MHC class II DLA-DQA1 exon 2 in golden jackal populations from Bulgaria, Serbia, and Hungary. Specifically, we tested for positive selection on and geographic variation at that locus due to adaptation to supposedly regionally varying pathogenic landscapes. To test for potential fitness effects of different protein variants on individual body condition, we used linear modeling of individual body mass indexes (bmi) and accounted for possible age, sex, geographical, and climatic effects. The latter approach was performed, however, only on Serbian individuals with appropriate data.

Results

Only three different DLA-DQA1 alleles were detected, all coding for different amino-acid sequences. The neutrality tests revealed no significant but positive values; there was no signal of spatial structuring and no deviation from the Hardy–Weinberg equilibrium across the studied range of expansion. However, we found a signal of trans-species polymorphism and significant test results for positive selection on three codons. Our information-theory based linear modeling results indicated an effect of ambient temperature on the occurrence of individual DLA-DQA1 genotypes in individuals from across the studied expansion range, independent from geographical position. Our linear modeling results of individual bmi values indicated that yearlings homozygous for DLA-DQA1*03001 reached values typical for adults contrary to yearlings carrying other genotypes (protein combinations). This suggested better growth rates and thus a possible fitness advantage of yearlings homozygous for DLA-DQA1*03001.

Conclusions

Our results indicate a demographic (stochastic) signal of reduced DLA-DQA1 exon 2 variation, in line with the documented historical demographic bottleneck. At the same time, however, allelic variation was also affected by positive selection and adaptation to varying ambient temperature, supposedly reflecting geographic variation in the pathogenic landscape. Moreover, an allele effect on body mass index values of yearlings suggested differential fitness associated with growth rates. Overall, a combination of a stochastic effect and positive selection has shaped and is still shaping the variation at the studied MHC locus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01856-z.

Nine residues in HLA-DQ molecules determine with susceptibility and resistance to type 1 diabetes among young children in Sweden

HLA-DQ molecules account over 50% genetic risk of type 1 diabetes (T1D), but little is known about associated residues. Through next generation targeted sequencing technology and deep learning of DQ residue sequences, the aim was to uncover critical residues and their motifs associated with T1D. Our analysis uncovered (αa1, α44, α157, α196) and (β9, β30, β57, β70, β135) on the HLA-DQ molecule. Their motifs captured all known susceptibility and resistant T1D associations. Three motifs, “DCAA-YSARD” (OR = 2.10, p = 1.96*10⁻²⁰), “DQAA-YYARD” (OR = 3.34, 2.69*10⁻⁷²) and “DQDA-YYARD” (OR = 3.71, 1.53*10⁻⁶) corresponding to DQ2.5 and DQ8.1 (the latter two motifs) associated with susceptibility. Ten motifs were significantly associated with resistance to T1D. Collectively, homozygous DQ risk motifs accounted for 43% of DQ-T1D risk, while homozygous DQ resistant motifs accounted for 25% protection to DQ-T1D risk. Of the identified nine residues five were within or near anchoring pockets of the antigenic peptide (α44, β9, β30, β57 and β70), one was the N-terminal of the alpha chain (αa1), one in the CD4-binding region (β135), one in the putative cognate TCR-induced αβ homodimerization process (α157), and one in the intra-membrane domain of the alpha chain (α196). Finding these critical residues should allow investigations of fundamental properties of host immunity that underlie tolerance to self and organ-specific autoimmunity.

Association of Bovine Leukemia Virus-Induced Lymphoma with BoLA-DRB3 Polymorphisms at DNA, Amino Acid, and Binding Pocket Property Levels

Bovine leukemia virus (BLV) causes enzootic bovine leucosis, a malignant B-cell lymphoma in cattle. The DNA sequence polymorphisms of bovine leukocyte antigen (BoLA)-DRB3 have exhibited a correlation with BLV-induced lymphoma in Holstein cows. However, the association may vary between different cattle breeds. Furthermore, little is known about the relationship between BLV-induced lymphoma and DRB3 at the amino acid and structural diversity levels. Here, we comprehensively analyzed the correlation between BLV-induced lymphoma and DRB3 at DNA, amino acid, and binding pocket property levels, using 106 BLV-infected asymptomatic and 227 BLV-induced lymphoma Japanese black cattle samples. DRB3*011:01 was identified as a resistance allele, whereas DRB3*005:02 and DRB3*016:01 were susceptibility alleles. Amino acid association studies showed that positions 9, 11, 13, 26, 30, 47, 57, 70, 71, 74, 78, and 86 were associated with lymphoma susceptibility. Structure and electrostatic charge modeling further indicated that binding pocket 9 of resistance DRB3 was positively charged. In contrast, alleles susceptible to lymphoma were neutrally charged. Altogether, this is the first association study of BoLA-DRB3 polymorphisms with BLV-induced lymphoma in Japanese black cattle. In addition, our results further contribute to understanding the mechanisms regarding how BoLA-DRB3 polymorphisms mediate susceptibility to BLV-induced lymphoma.

HLA-DRB1 Alleles Associated with Lower Leishmaniasis Susceptibility Share Common Amino Acid Polymorphisms and Epitope Binding Repertoires

Susceptibility for leishmaniasis is largely dependent on host genetic and immune factors. Despite the previously described association of human leukocyte antigen (HLA) gene cluster variants as genetic susceptibility factors for leishmaniasis, little is known regarding the mechanisms that underpin these associations. To better understand this underlying functionality, we first collected all known leishmaniasis-associated HLA variants in a thorough literature review. Next, we aligned and compared the protection- and risk-associated HLA-DRB1 allele sequences. This identified several amino acid polymorphisms that distinguish protection- from risk-associated HLA-DRB1 alleles. Subsequently, T cell epitope binding predictions were carried out across these alleles to map the impact of these polymorphisms on the epitope binding repertoires. For these predictions, we used epitopes derived from entire proteomes of multiple Leishmania species. Epitopes binding to protection-associated HLA-DRB1 alleles shared common binding core motifs, mapping to the identified HLA-DRB1 amino acid polymorphisms. These results strongly suggest that HLA polymorphism, resulting in differential antigen presentation, affects the association between HLA and leishmaniasis disease development. Finally, we established a valuable open-access resource of putative epitopes. A set of 14 HLA-unrestricted strong-binding epitopes, conserved across species, was prioritized for further epitope discovery in the search for novel subunit-based vaccines.

Depletion of MHC supertype during domestication can compromise immunocompetence

The major histocompatibility complex (MHC) plays an important role in infectious disease resistance. The presence of certain MHC alleles and functionally similar groups of MHC alleles (i.e., supertypes) has been associated with resistance to particular parasite species. Farmed and domesticated fish stocks are often depleted in their MHC alleles and supertype diversity, possibly as a consequence of artificial selection for desirable traits, inbreeding (loss of heterozygosity), genetic drift (loss of allelic diversity) and/or reduced parasite biodiversity. Here we quantify the effects of depletion of MHC class II genotype and supertype variation on resistance to the parasite Gyrodactylus turnbulli in guppies (Poecilia reticulata). Compared to the descendants of wild‐caught guppies, ornamental fish had a significantly reduced MHC variation (i.e., the numbers of MHC alleles and supertypes per individual, and per population). In addition, ornamental fish were significantly more susceptible to G. turnbulli infections, accumulating peak intensity 10 times higher than that of their wildtype counterparts. Four out of 13 supertypes were associated with a significantly reduced parasite load, and the presence of some supertypes had a dramatic effect on the intensity of infection. Remarkably, the ornamental and wildtype fish differed in the supertypes that were associated with parasite resistance. Analysis with a genetic algorithm showed that resistance‐conferring supertypes of the wildtype and ornamental fish shared two unique amino acids in the peptide‐binding region of the MHC that were not found in any other alleles. These data show that the supertype demarcation captures some, but not all, of the variation in the immune function of the alleles. This study highlights the importance of managing functional MHC diversity in livestock, and suggests there might be some immunological redundancy among MHC supertypes.

Cartilaginous fish class II genes reveal unprecedented old allelic lineages and confirm the late evolutionary emergence of DM

Cartilaginous fish (chimaeras, rays and sharks) are the most basal extant jawed vertebrates with an adaptive immune system based on the Major Histocompatibility Complex (MHC). Despite being a key taxon in the evolution of vertebrate adaptive immunity, no comprehensive characterization of MHC class II genes has been undertaken for the group. We performed extensive bioinformatic searches on a taxonomically diverse dataset of transcriptomes and genomes of cartilaginous fish targeting MHC class II sequences. Class IIα and IIβ sequences were retrieved from all taxa analyzed and showed typical features of classical class II genes. Phylogenetic trees of the immunoglobulin superfamily domain showed two divergent and remarkably ancient lineages of class II genes in Selachians (sharks), originating >350 million years ago. Close linkage of lineage-specific pairs of IIα and IIβ genes was found, confirming previous results, with genes from distinct lineages segregating as alleles. Nonclassical class II DM sequences were not retrieved from these data and classical class II sequences lacked the conserved residues shown to interact with DM molecules, supporting claims that the DM system arose only in the lobe-finned fish lineage leading to tetrapods. Based on our search methods, other divergent class II genes are unlikely in cartilaginous fish.

Enabling Routine MHC-II-Associated Peptide Proteomics for Risk Assessment of Drug-Induced Immunogenicity

Major histocompatibility complex-II (MHC-II)-Associated Peptide Proteomics (MAPPs) is a mass spectrometry-based approach to identify and relatively quantitate naturally processed and presented MHC-II-associated peptides that can potentially activate T cells and contribute to the immunogenicity of a drug. Acceptance of the MAPPs technology as an appropriate preclinical (and potentially clinical) immunogenicity risk assessment tool depends not only on its technical stability and robustness but also on the ability to compare results across experiments and donors. To this end, we developed a specialized MAPPs data processing pipeline, dataMAPPs, which presents complex mass spectrometric data sets in the form of heat maps (heatMAPPs), enabling rapid and convenient comparison between conditions and donors. A customized normalization procedure based on identified endogenous peptides standardizes signal intensities within and between donors and enables cross-experimental comparison. We evaluated the technical reproducibility of the MAPPs platform using tool compounds with respect to the most prominent experimental factors and found that the systematic biological differences across donors by far outweighed any technical source of variation. We illustrate the capability of the MAPPs platform to generate data that may be used for preclinical risk assessment of drug-induced immunogenicity and discuss its applicability in the clinics.

Diversity analysis at MHC class II DQA locus in buffalo (Bubalus bubalis) indicates extensive duplication and trans-species evolution

Variation at MHC Class II-DQA locus in riverine and swamp buffaloes (Bubu) has been explored in this study. Through sequencing of buffalo DQA, 48 nucleotide variants identified from 17 individuals, reporting 42 novel alleles, including one pseudogene. Individual animal displayed two to seven variants, suggesting the presence of more than two Bubu-DQA loci, as an evidence of extensive duplication. dN values were found to be higher than dS values at peptide binding sites, separately for riverine and swamp buffaloes, indicating locus being under positive selection. Evolutionary analysis revealed numerous trans-species polymorphism with alleles from water buffalo assigned to at least three different loci (Bubu-DQA1, DQA2, DQA3). Alleles of both the sub-species intermixed within the cluster, showing convergent evolution of MHC alleles in bovines. The results thus suggest that both riverine and swamp buffaloes share con-current arrangement of DQA region, comparable to cattle in terms of copy number and population polymorphism.

Motifs of Three HLA-DQ Amino Acid Residues (α44, β57, β135) Capture Full Association With the Risk of Type 1 Diabetes in DQ2 and DQ8 Children

HLA-DQA1 and -DQB1 are strongly associated with type 1 diabetes (T1D), and DQ8.1 and DQ2.5 are major risk haplotypes. Next-generation targeted sequencing of HLA-DQA1 and -DQB1 in Swedish newly diagnosed 1- to 18 year-old patients (n = 962) and control subjects (n = 636) was used to construct abbreviated DQ haplotypes, converted into amino acid (AA) residues, and assessed for their associations with T1D. A hierarchically organized haplotype (HOH) association analysis allowed 45 unique DQ haplotypes to be categorized into seven clusters. The DQ8/9 cluster included two DQ8.1 risk and the DQ9 resistant haplotypes, and the DQ2 cluster included the DQ2.5 risk and DQ2.2 resistant haplotypes. Within each cluster, HOH found residues α44Q (odds ratio [OR] 3.29, P = 2.38 * 10^-85) and β57A (OR 3.44, P = 3.80 * 10^-84) to be associated with T1D in the DQ8/9 cluster representing all ten residues (α22, α23, α44, α49, α51, α53, α54, α73, α184, β57) due to complete linkage disequilibrium (LD) of α44 with eight such residues. Within the DQ2 cluster and due to LD, HOH analysis found α44C and β135D to share the risk for T1D (OR 2.10, P = 1.96 * 10^-20). The motif "QAD" of α44, β57, and β135 captured the T1D risk association of DQ8.1 (OR 3.44, P = 3.80 * 10^-84), and the corresponding motif "CAD" captured the risk association of DQ2.5 (OR 2.10, P = 1.96 * 10^-20). Two risk associations were related to GAD65 autoantibody (GADA) and IA-2 autoantibody (IA-2A) but in opposite directions. CAD was positively associated with GADA (OR 1.56, P = 6.35 * 10^-8) but negatively with IA-2A (OR 0.59, P = 6.55 * 10^-11). QAD was negatively associated with GADA (OR 0.88; P = 3.70 * 10^-3) but positively with IA-2A (OR 1.64; P = 2.40 * 10^-14), despite a single difference at α44. The residues are found in and around anchor pockets 1 and 9, as potential T-cell receptor contacts, in the areas for CD4 binding and putative homodimer formation. The identification of three HLA-DQ AAs (α44, β57, β135) conferring T1D risk should sharpen functional and translational studies.

Depauperate major histocompatibility complex variation in the endangered reticulated flatwoods salamander (Ambystoma bishopi)

Reticulated flatwoods salamander (Ambystoma bishopi) populations began decreasing dramatically in the 1900s. Contemporary populations are small, isolated, and may be susceptible to inbreeding and reduced adaptive potential because of low genetic variation. Genetic variation at immune genes is especially important as it influences disease susceptibility and adaptation to emerging infectious pathogens, a central conservation concern for declining amphibians. We collected samples from across the extant range of this salamander to examine genetic variation at major histocompatibility complex (MHC) class Iα and IIβ exons as well as the mitochondrial control region. We screened tail or toe tissue for ranavirus, a pathogen associated with amphibian declines worldwide. Overall, we found low MHC variation when compared to other amphibian species and did not detect ranavirus at any site. MHC class Iα sequencing revealed only three alleles with a nucleotide diversity of 0.001, while MHC class IIβ had five alleles with a with nucleotide diversity of 0.004. However, unique variation still exists across this species' range with private alleles at three sites. Unlike MHC diversity, mitochondrial variation was comparable to levels estimated for other amphibians with nine haplotypes observed, including one haplotype shared across all sites. We hypothesize that a combination of a historic disease outbreak and a population bottleneck may have contributed to low MHC diversity while maintaining higher levels of mitochondrial DNA variation. Ultimately, MHC data indicated that the reticulated flatwoods salamander may be at an elevated risk from infectious diseases due to low levels of immunogenetic variation necessary to combat novel pathogens.

Molecular evolution of MHC class II DRB exon 2 in Japanese and Russian raccoon dogs, Nyctereutes procyonoides (Carnivora: Canidae)

Raccoon dogs, Nyctereutes procyonoides, are native to East Asia, but have been introduced into western Russia and eastern Europe. To determine allelic diversity and elucidate the evolution of major histocompatibility complex (MHC) class II genes in the raccoon dog, we analysed a 237-bp region of DRB exon 2 from 36 individuals of native and introduced populations from Japan and Russia. We detected 23 DRB alleles (Nypr-DRBs), 22 of which were novel. Some alleles were found across the species’ range, while others were geographically restricted. For both native and introduced populations, the ratio of non-synonymous to synonymous substitution rates for codons at predicted antigen-binding sites was significantly greater than 2, indicating that Nypr-DRBs have evolved under positive selection. Mixed effect model evolution analysis and an algorithm to detect recombination showed five positively selected codons and one recombination breakpoint, respectively. Overall, our results suggest that the diversity of MHC class II DRB in N. procyonoides was influenced and maintained by recombination, pathogen-driven positive selection, geographical barriers and the founder effect. A Bayesian phylogenetic tree revealed no evidence of trans-species polymorphism (TSP), but instead showed monophyly for the Nypr-DRB alleles within a larger clade of canid sequences. The lack of TSP may have been due to long-term divergence of raccoon dogs from other canids, or to their having encountered different sets of pathogens due to occupying a different ecological niche.

Ligandomes obtained from different HLA-class II-molecules are homologous for N- and C-terminal residues outside the peptide-binding cleft

Human CD4+ T lymphocytes play an important role in inducing potent immune responses. T cells are activated and stimulated by peptides presented in human leucocyte antigen (HLA)-class II molecules. These HLA-class II molecules typically present peptides of between 12 and 20 amino acids in length. The region that interacts with the HLA molecule, designated as the peptide-binding core, is highly conserved in the residues which anchor the peptide to the molecule. In addition, as these peptides are the product of proteolytic cleavages, certain conserved residues may be expected at the N- and C-termini outside the binding core. To study whether similar conserved residues are present in different cell types, potentially harbouring different proteolytic enzymes, the ligandomes of HLA-DRB1*03:01/HLA-DRB > 1 derived from two different cell types (dendritic cells and EBV-transformed B cells) were identified with mass spectrometry and the binding core and N- and C-terminal residues of a total of 16,568 peptides were analysed using the frequencies of the amino acids in the human proteome. Similar binding motifs were found as well as comparable conservations in the N- and C-terminal residues. Furthermore, the terminal conservations of these ligandomes were compared to the N- and C-terminal conservations of the ligandome acquired from dendritic cells homozygous for HLA-DRB1*04:01. Again, comparable conservations were evident with only minor differences. Taken together, these data show that there are conservations in the terminal residues of peptides, presumably the result of the activity of proteases involved in antigen processing.

Eleven Amino Acids of HLA-DRB1 and Fifteen Amino Acids of HLA-DRB3, 4, and 5 Include Potentially Causal Residues Responsible for the Risk of Childhood Type 1 Diabetes

Next-generation targeted sequencing of HLA-DRB1 and HLA-DRB3, -DRB4, and -DRB5 (abbreviated as DRB345) provides high resolution of functional variant positions to investigate their associations with type 1 diabetes risk and with autoantibodies against insulin (IAA), GAD65 (GADA), IA-2 (IA-2A), and ZnT8 (ZnT8A). To overcome exceptional DR sequence complexity as a result of high polymorphisms and extended linkage disequilibrium among the DR loci, we applied a novel recursive organizer (ROR) to discover disease-associated amino acid residues. ROR distills disease-associated DR sequences and identifies 11 residues of DRB1, sequences of which retain all significant associations observed by DR genes. Furthermore, all 11 residues locate under/adjoining the peptide-binding groove of DRB1, suggesting a plausible functional mechanism through peptide binding. The 15 residues of DRB345, located respectively in the β49-55 homodimerization patch and on the face of the molecule shown to interact with and bind to the accessory molecule CD4, retain their significant disease associations. Further ROR analysis of DR associations with autoantibodies finds that DRB1 residues significantly associated with ZnT8A and DRB345 residues with GADA. The strongest association is between four residues (β14, β25, β71, and β73) and IA-2A, in which the sequence ERKA confers a risk association (odds ratio 2.15, P = 10^-18), and another sequence, ERKG, confers a protective association (odds ratio 0.59, P = 10^-11), despite a difference of only one amino acid. Because motifs of identified residues capture potentially causal DR associations with type 1 diabetes, this list of residuals is expected to include corresponding causal residues in this study population.

Discriminative T cell recognition of cross-reactive islet-antigens is associated with HLA-DQ8 transdimer-mediated autoimmune diabetes

Human leukocyte antigen (HLA)-DQ8 transdimer (HLA-DQA1*0501/DQB1*0302) confers exceptionally high risk in autoimmune diabetes. However, little is known about HLA-DQ8 transdimer-restricted CD4 T cell recognition, an event crucial for triggering HLA-DQ8 transdimer-specific anti-islet immunity. Here, we report a high degree of epitope overlap and T cell promiscuity between susceptible HLA-DQ8 and HLA-DQ8 transdimer. Despite preservation of putative residues for T cell receptor (TCR) contact, stronger disease-associated responses to cross-reactive, immunodominant islet epitopes are elicited by HLA-DQ8 transdimer. Mutagenesis at the α chain of HLA-DQ8 transdimer in complex with the disease-relevant GAD65_250-266 peptide and in silico analysis reveal the DQ α52 residue located within the N-terminal edge of the peptide-binding cleft for the enhanced T cell reactivity, altering avidity and biophysical affinity between TCR and HLA-peptide complexes. Accordingly, a structurally promiscuous but nondegenerate TCR-HLA-peptide interface is pivotal for HLA-DQ8 transdimer-mediated autoimmune diabetes.

Risk HLA class II alleles and amino acid residues in myeloperoxidase-ANCA-associated vasculitis

A genome-wide association study (GWAS) indicated that myeloperoxidase-ANCA associated vasculitis (AAV) is associated with HLA-DQ. However, susceptibility alleles in these loci have been under-investigated. Here we genotyped 258 Chinese patients with myeloperoxidase-AAV and 597 healthy control individuals at HLA DRB1, DQA1, DQB1 and DPB1, and extracted the encoded amino acid sequences from the IMGT/HLA database. The replication cohort included 97 cases and 107 controls. T cell epitopes of myeloperoxidase were predicted and docked to the HLA molecules. We found DQA1∗0302 (odds ratio 2.34 (95% confidence interval 1.75-3.14)) and DQB1∗0303 (odds ratio 1.89 (1.45-2.48)) were risk alleles for myeloperoxidase-AAV. They are in overt linkage disequilibrium (r² 0.69) and the haplotype DQA1∗0302-DQB1∗0303 presents a significant risk (haplotype score 6.39) as well. Aspartate160 on the DQ α chain (odds ratio 2.06 (1.60-2.67)), encoded by DQA1∗0302, and isoleucine185 on the DQ β chain (odds ratio 1.73 (1.38-2.18)), encoded by DQB1∗0303, both located in the α2β2 domains, conferred significant risk for myeloperoxidase-AAV. Homologous modeling showed that DQα∗160D may confer susceptibility to myeloperoxidase-AAV by altering dimerization of the HLA molecules. Thus, more attention should be paid to the roles of amino acids in the α2β2 domains in addition to the α1β1 binding groove of HLA class II molecules.

Characterisation of major histocompatibility complex class IIa haplotypes in an island sheep population

The ovine MHC class IIa is known to consist of six to eight loci located in close proximity on chromosome 20, forming haplotypes that are typically inherited without recombination. Here, we characterise the class IIa haplotypes within the Soay sheep (Ovis aries) on St. Kilda to assess the diversity present within this unmanaged island population. We used a stepwise sequence-based genotyping strategy to identify alleles at seven polymorphic MHC class IIa loci in a sample of 118 Soay sheep from four cohorts spanning 15 years of the long-term study on St. Kilda. DRB1, the most polymorphic MHC class II locus, was characterised first in all 118 sheep and identified six alleles. Using DRB1 homozygous animals, the DQA (DQA1, DQA2 and DQA2-like) and DQB (DQB1, DQB2 and DQB2-like) loci were sequenced, revealing eight haplotypes. Both DQ1/DQ2 and DQ2/DQ2-like haplotype configurations were identified and a single haplotype carrying three DQB alleles. A test sample of 94 further individuals typed at the DRB1 and DQA loci found no exceptions to the eight identified haplotypes and a haplotype homozygosity of 21.3%. We found evidence of historic positive selection at DRB1, DQA and DQB. The limited variation at MHC class IIa loci in Soay sheep enabled haplotype characterisation but showed that no single locus could capture the full extent of the expressed variation in the region.

Islands and hybrid zones: combining the knowledge from "Natural Laboratories" to explain phylogeographic patterns of the European brown hare

Background

The aim of the study was to use hybrid populations as well as island populations of the European brown hare (Lepus europaeus) to explore the effect of evolutionary events, such as the post-deglaciation translocations, spontaneous and human-mediated, local adaptation and the genetic drift in the shaping of the phylogeographic patterns of the species. For this purpose, we used molecular markers, both nuclear and mitochondrial, that are indicative for local adaptation as well as neutral markers to elucidate the patterns of population differentiation based on geographic isolation and the clade of origin. To broaden our analysis, we included data from our previous studies concerning mainland populations, to explore the genetic differentiation in the base of the geographic origin (mainland/island) of the populations.

Results

Our results suggest that local adaptation shapes the differentiation in both genomes, favoring specific alleles in nuclear genes (e.g. DQA) or haplotypes in mtDNA (e.g. Control Region, CR). mtDNA variation was found to be in a higher level and was able to give a phylogeographic signal for the populations. Furthermore, the degree of variation was influenced not only by the geographic origin, but also by the clade of origin, since specific island populations of Anatolian origin showed a greater degree of variation compared to specific mainland populations of the European clade. Concerning the hybrid population, we confirmed the existence of both clades in the territory and we provided a possible explanation for the lack of introgression between the clades.

Conclusion

Our results indicate that the Quaternary’s climatic oscillations played a major role in the shaping of the phylogeographic patterns of the species, by isolating populations in the distinct refugia, where they adapted and differentiate in allopatry, leading to genome incompatibilities observed nowadays.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1354-y) contains supplementary material, which is available to authorized users.

DRB4*01:01 Has a Distinct Motif and Presents a Proinsulin Epitope That Is Recognized in Subjects with Type 1 Diabetes

DRB4*01:01 (DRB4) is a secondary HLA-DR product that is part of the high-risk DR4/DQ8 haplotype that is associated with type 1 diabetes (T1D). DRB4 shares considerable homology with HLA-DR4 alleles that predispose to autoimmunity, including DRB1*04:01 and DRB1*04:04. However, the DRB4 protein sequence includes distinct residues that would be expected to alter the characteristics of its binding pockets. To identify high-affinity peptides that are recognized in the context of DRB4, we used an HLA class II tetramer-based approach to identify epitopes within multiple viral Ags. We applied a similar approach to identify antigenic sequences within glutamic acid decarboxylase 65 and pre-proinsulin that are recognized in the context of DRB4. Seven sequences were immunogenic, eliciting high-affinity T cell responses in DRB4⁺ subjects. DRB1*04:01-restricted responses toward many of these peptides have been previously described, but responses to a novel pre-proinsulin 9-28 peptide were commonly observed in subjects with T1D. Furthermore, T cells that recognized this peptide in the context of DRB4 were present at significantly higher frequencies in patients with T1D than in healthy controls, implicating this as a disease-relevant specificity that may contribute to the breakdown of β cell tolerance in genetically susceptible individuals. We then deduced a DRB4 motif and confirmed its key features through structural modeling. This modeling suggested that the core epitope within the pre-proinsulin 9-28 peptide has a somewhat unusual binding motif, with tryptophan in the fourth binding pocket of DRB4, perhaps influencing the availability of this complex for T cell selection.

Amino acid signatures in the HLA class II peptide-binding region associated with protection/susceptibility to the severe West Nile Virus disease

The MHC class II region in humans is highly polymorphic. Each MHC molecule is formed by an α and a β chain, produced by different genes, creating an antigen-binding groove. In the groove there are several pockets into which antigens anchor and fit. The affinity of this fitting determines the recognition specificity of a given peptide. Here, based on our previous results about the association of MHC class II with the WNV disease, we examined the role of the binding pockets of HLA-DPA1, -DQA1 and-DRB1 in the severe form of the disease. In HLA-DQA1, variants in all pockets 1, 6 and 9 were found to be associated with either protection and/or susceptibility to neuroinvasion caused by WNV. Similarly, pockets 7, 9 and 10 in HLA-DRB1 were associated with severe disease. Protein modeling of these molecules revealed structural and functional differences among alleles with opposite roles concerning the development of the disease. Different amino acids in positions α52 and α66 (HLA-DQA1) significantly influenced the peptide binding while DYWLR/EFA combination (HLA-DRB1) was associated with neuronal damage. Further studies could help us understand the selectivity of pocket variants in order to create suitable peptides for an effective response.

Positive selection in coding regions and motif duplication in regulatory regions of bottlenose dolphin MHC class II genes

The vertebrate immune response is mediated through highly adaptive, quickly evolving cell surface receptors, the major histocompatibility complex (MHC). MHC molecules bind and present a diverse array of pathogenic molecules and trigger a cascade of defenses. Use of MHC variation as a marker for population health has also evolved quickly following advances in sequencing methods. We applied a combination of traditional and next generation sequencing methodology to characterize coding (peptide binding region) and regulatory (proximal promoter) sequence variation in MHC Class II DQA and DQB genes between estuarine and coastal populations of the bottlenose dolphin, Tursiops truncatus, an apex predator whose health status is indicative of anthropogenic impacts on the ecosystem. The coding regions had 10 alleles each at DQA and DQB; the promoters had 6 and 7 alleles at DQA and DQB, respectively with variation within key regulatory motifs. Positive selection was observed for the coding regions of both genes while both coding and promoter regions exhibited geographic differences in allele composition that likely indicates diversifying selection across habitats. Most notable was the discovery of a complete duplication of a 14-bp T-box motif in the DQA promoter. Four class II promoter regions (DQA, DQB, DRA, DRB) were characterized in species from four cetacean families (Delphinidae, Monodontidae, Lipotidae, and Physeteridae) and revealed substantial promoter structural diversity across this order. Peptide binding regions may not be the only source of adaptive potential within cetacean MHC for responding to pathogenic threats. These findings are the first analysis of cetacean MHC regulatory motifs, which may divulge unique immunogenetic strategies among cetaceans and reveal how MHC transcriptional control continues to evolve. The combined MHC regulatory and coding data provide new genetic context for distinct vulnerability profiles between coastal and estuarine populations, which are key concerns for health and risk management.

HLA class II alleles differing by a single amino acid associate with clinical phenotype and outcome in patients with primary membranous nephropathy

Genome-wide associations and HLA genotyping have revealed associations between HLA alleles and susceptibility to primary membranous nephropathy. However, associations with clinical phenotypes and kidney outcome are poorly defined. We previously identified DRB1*1501 and DRB1*0301 as independent risk alleles for primary membranous nephropathy. Here, we investigated HLA associations with demographic characteristics, anti-phospholipase A2 receptor (PLA2R) antibody, treatment response and kidney outcome after a median follow-up of 52 months in 258 patients. DRB1*0301, but not DRB1*1501, was associated with a significantly higher level of PLA2R antibody (odds ratio 1.58, 95% confidence interval 1.13-2.22). Although DRB1*1502, which differs from DRB1*1501 by a single amino acid, was not a risk allele for primary membranous nephropathy (odds ratio 1.01), it was associated with significantly lower estimated glomerular filtration rates both at baseline (1.79, 1.18-2.72) and at last follow-up (1.72, 1.17-2.53), a significantly worse renal outcome by Kaplan-Meier analysis and a significantly higher risk of end-stage renal disease by Cox regression analysis (hazard ratio 4.52, 1.22-16.74). Nevertheless, the absence of remission remained the only independent risk factor for end-stage renal disease by multivariate analysis. DRB1*1502 was also associated with a significantly higher median PLA2R antibody level [161.4 vs. 36.3 U/mL] and showed interaction with DRB1*0301 for this variable. Thus, HLA genes control PLA2R antibody production and primary membranous nephropathy severity and outcome. Additionally, DRB1*1502 behaves like a modifier gene with a strong predictor value when associated with HLA risk alleles. Other modifier genes need further investigations in larger cohorts.

Positive selection and climatic effects on MHC class II gene diversity in hares (Lepus capensis) from a steep ecological gradient

In natural populations, allelic diversity of the major histocompatibility complex (MHC) is commonly interpreted as resulting from positive selection in varying spatiotemporal pathogenic landscapes. Composite pathogenic landscape data are, however, rarely available. We studied the spatial distribution of allelic diversity at two MHC class II loci (DQA, DQB) in hares, Lepus capensis, along a steep ecological gradient in North Africa and tested the role of climatic parameters for the spatial distribution of DQA and DQB proteins. Climatic parameters were considered to reflect to some extent pathogenic landscape variation. We investigated historical and contemporary forces that have shaped the variability at both genes, and tested for differential selective pressure across the ecological gradient by comparing allelic variation at MHC and neutral loci. We found positive selection on both MHC loci and significantly decreasing diversity from North to South Tunisia. Our multinomial linear models revealed significant effects of geographical positions that were correlated with mean annual temperature and precipitation on the occurrence of protein variants, but no effects of co-occurring DQA or DQB proteins, respectively. Diversifying selection, recombination, adaptation to local pathogenic landscapes (supposedly reflected by climate parameters) and neutral demographic processes have shaped the observed MHC diversity and differentiation patterns.

Epitope-Binding Characteristics for Risk versus Protective DRB1 Alleles for Visceral Leishmaniasis

HLA-DRB1 is the major genetic risk factor for visceral leishmaniasis (VL). We used SNP2HLA to impute HLA-DRB1 alleles and SNPTEST to carry out association analyses in 889 human cases and 977 controls from India. NetMHCIIpan 2.1 was used to map epitopes and binding affinities across 49 Leishmania vaccine candidates, as well as across peptide epitopes captured from dendritic cells treated with crude Leishmania Ag and identified using mass spectrometry and alignment to amino acid sequences of a reference Leishmania genome. Cytokines were measured in peptide-stimulated whole blood from 26 cured VL cases and eight endemic healthy controls. HLA-DRB1*1501 and DRB1*1404/DRB1*1301 were the most significant protective and risk alleles, respectively, with specific residues at aa positions 11 and 13 unique to protective alleles. We observed greater peptide promiscuity in sequence motifs for 9-mer core epitopes predicted to bind to risk (*1404/*1301) compared with protective (*1501) DRB1 alleles. There was a higher frequency of basic amino acids in DRB1*1404/*1301-specific epitopes compared with hydrophobic and polar amino acids in DRB1*1501-specific epitopes at anchor residues pocket 4 and pocket 6, which interact with residues at DRB1 positions 11 and 13. Cured VL patients made variable, but robust, IFN-γ, TNF, and IL-10 responses to 20-mer peptides based on captured epitopes, with peptides based on DRB1*1501-captured epitopes resulting in a higher proportion (odds ratio 2.23, 95% confidence interval 1.17-4.25, p = 0.017) of patients with IFN-γ/IL-10 ratios > 2-fold compared with peptides based on DRB1*1301-captured epitopes. Our data provide insight into the molecular mechanisms underpinning the association of HLA-DRB1 alleles with risk versus protection in VL in humans.

Combined effect of glutamine at position 70 of HLA-DRB1 and alanine at position 57 of HLA-DQB1 in type 1 diabetes: An epitope analysis

The contribution of specific HLA Class II alleles in type 1 diabetes is determined by polymorphic amino acid epitopes that direct antigen binding therefore, along with conventional allele frequency analysis, epitope analysis can provide important insights into disease susceptibility. We analyzed the highly heterogeneous Cypriot population for the HLA class II loci of T1DM patients and controls and we report for the first time their allele frequencies. Within our patient cohort we identified a subgroup that did not carry the DRB1*03:01-DQA1*05:01-DQB1*02:01 and DRB1*04:xx-DQA1*03:01-DQB1*03:02 risk haplotypes but a novel recombinant one, DRB1*04:XX-DQA1*03:01-DQB1*02:01 designated DR4-DQ2.3. Through epitope analysis we identified established susceptibility (DQB1 A⁵⁷, DRB1 H¹³) and resistance (DQB1 D⁵⁷) residues as well as other novel susceptibility residues DRB1 Q⁷⁰, DQB1 L²⁶ and resistance residues DRB1 D⁷⁰, R⁷⁰ and DQB1 Y⁴⁷. Prevalence of susceptibility epitopes was higher in patients and was not exclusively a result of linkage disequilibrium. Residues DRB1 Q⁷⁰, DQB1 L²⁶ and A⁵⁷ and a 10 amino acid epitope of DQA1 were the most significant in discriminating risk alleles. An extended haplotype containing these epitopes was carried by 92% of our patient cohort. Sharing of susceptibility epitopes could also explain the absence of risk haplotypes in patients. Finally, many significantly associated epitopes were non-pocket residues suggesting that critical immune functions may exist spanning further from the binding pockets.

Signatures of historical selection on MHC reveal different selection patterns in the moor frog (Rana arvalis)

MHC genes are key components in disease resistance and an excellent system for studying selection acting on genetic variation in natural populations. Current patterns of variation in MHC genes are likely to be influenced by past and ongoing selection as well as demographic fluctuations in population size such as those imposed by post-glacial recolonization processes. Here, we investigated signatures of historical selection and demography on an MHC class II gene in 12 moor frog populations along a 1700-km latitudinal gradient. Sequences were obtained from 207 individuals and consecutively assigned into two different clusters (northern and southern clusters, respectively) in concordance with a previously described dual post-glacial colonization route. Selection analyses comparing the relative rates of non-synonymous to synonymous substitutions (dN/dS) suggested evidence of different selection patterns in the northern and the southern clusters, with divergent selection prevailing in the south but uniform positive selection predominating in the north. Also, models of codon evolution revealed considerable differences in the strength of selection: The southern cluster appeared to be under strong selection while the northern cluster showed moderate signs of selection. Our results indicate that the MHC alleles in the north diverged from southern MHC alleles as a result of differential selection patterns.

The Implications of DNA Methylation on Food Allergy

Food allergy is a major clinical and public health concern worldwide. The risk factors are well defined, however, the mechanisms by which they affect immune development remain largely unknown, and unfortunately the effective treatment or prevention of food allergy is still being researched. Recent studies show that the genes that are critical for the development of food allergy are regulated through DNA methylation. Environmental factors can affect host DNA methylation status and subsequently predispose people to food allergy. DNA methylation is therefore an important mediator of gene-environment interactions in food allergy and key to understanding the mechanisms underlying the allergic development. Indeed, the modification and identification of the methylation levels of specific genetic loci have gained increasing attention for therapeutic and diagnostic application in combating food allergy. In this review, we summarize and discuss the recent developments of DNA methylation in food allergy, including the pathogenesis, therapy, and diagnosis. This review will also summarize and discuss the environmental factors that affect DNA methylation levels in food allergy.

HLA class II and rheumatoid arthritis: the bumpy road of revelation

Rheumatoid arthritis (RA) is a chronic auto-immune disease primarily targeting the joints. Approximately 1% of the population is affected by RA, and despite the improvements in therapeutic interventions, elucidation of the disease pathogenesis is still in its infancy. RA patients can be subdivided on basis of the presence of autoantibodies, especially anti-citrullinated protein antibodies (ACPA). ACPA⁺ and ACPA^- disease most likely differ in aetiology, as different genetic and environmental risk factors are associated with these two disease entities. For ACPA⁺ RA disease, the genetic factors associating with disease mainly comprised of human leukocyte antigen (HLA) class II molecules. The predisposing HLA-DR alleles have been depicted as the 'HLA Shared Epitope (SE) alleles', as these alleles encode a similar sequence, the shared epitope sequence, within the beta chain of the HLA-DR molecule. In addition to the involvement of the HLA-SE alleles in the development of ACPA⁺ RA disease, other HLA-DR molecules have been shown to confer protection against this disease entity. The protective HLA molecules have, instead of the SE-motif, a different but shared sequence at the same location in the beta chain of HLA-DR molecules, consisting of the amino acid residues DERAA. The possible contributions of the predisposing and protective HLA molecules in association with ACPA-positive RA are discussed in this review.

Insights into the ancestral organisation of the mammalian MHC class II region from the genome of the pteropid bat, Pteropus alecto

BACKGROUND

Bats are an extremely successful group of mammals and possess a variety of unique characteristics, including their ability to co-exist with a diverse range of pathogens. The major histocompatibility complex (MHC) is the most gene dense and polymorphic region of the genome and MHC class II (MHC-II) molecules play a vital role in the presentation of antigens derived from extracellular pathogens and activation of the adaptive immune response. Characterisation of the MHC-II region of bats is crucial for understanding the evolution of the MHC and of the role of pathogens in shaping the immune system.

RESULTS

Here we describe the relatively contracted MHC-II region of the Australian black flying-fox (Pteropus alecto), providing the first detailed insight into the MHC-II region of any species of bat. Twelve MHC-II genes, including one locus (DRB2) located outside the class II region, were identified on a single scaffold in the bat genome. The presence of a class II locus outside the MHC-II region is atypical and provides evidence for an ancient class II duplication block. Two non-classical loci, DO and DM and two classical, DQ and DR loci, were identified in P. alecto. A putative classical, DPB pseudogene was also identified. The bat's antigen processing cluster, though contracted, remains highly conserved, thus supporting its importance in antigen presentation and disease resistance.

CONCLUSIONS

This detailed characterisation of the bat MHC-II region helps to fill a phylogenetic gap in the evolution of the mammalian class II region and is a stepping stone towards better understanding of the immune responses in bats to viral, bacterial, fungal and parasitic infections.