Three-dimensional structure of the human class II histocompatibility antigen HLA-DR1

Progression to type 1 diabetes in the DPT-1 and TN07 clinical trials is critically associated with specific residues in HLA-DQA1-B1 heterodimers

AIMS/HYPOTHESIS

The aim of this work was to explore molecular amino acids (AAs) and related structures of HLA-DQA1-DQB1 that underlie its contribution to the progression from stages 1 or 2 to stage 3 type 1 diabetes.

METHODS

Using high-resolution DQA1 and DQB1 genotypes from 1216 participants in the Diabetes Prevention Trial-Type 1 and the Diabetes Prevention Trial, we applied hierarchically organised haplotype association analysis (HOH) to decipher which AAs contributed to the associations of DQ with disease and their structural properties. HOH relied on the Cox regression to quantify the association of DQ with time-to-onset of type 1 diabetes.

RESULTS

By numerating all possible DQ heterodimers of α- and β-chains, we showed that the heterodimerisation increases genetic diversity at the cellular level from 43 empirically observed haplotypes to 186 possible heterodimers. Heterodimerisation turned several neutral haplotypes (DQ2.2, DQ2.3 and DQ4.4) to risk haplotypes (DQ2.2/2.3-DQ4.4 and DQ4.4-DQ2.2). HOH uncovered eight AAs on the α-chain (-16α, -13α, -6α, α22, α23, α44, α72, α157) and six AAs on the β-chain (-18β, β9, β13, β26, β57, β135) that contributed to the association of DQ with progression of type 1 diabetes. The specific AAs concerned the signal peptide (minus sign, possible linkage to expression levels), pockets 1, 4 and 9 in the antigen-binding groove of the α1β1 domain, and the putative homodimerisation of the αβ heterodimers.

CONCLUSIONS/INTERPRETATION

These results unveil the contribution made by DQ to type 1 diabetes progression at individual residues and related protein structures, shedding light on its immunological mechanisms and providing new leads for developing treatment strategies.

DATA AVAILABILITY

Clinical trial data and biospecimen samples are available through the National Institute of Diabetes and Digestive and Kidney Diseases Central Repository portal ( https://repository.niddk.nih.gov/studies ).

Permissible HLA mismatches in 9/10 unrelated donor pediatric stem cell transplants using HLA-EMMA: an EBMT Inborn Errors Working Party study

ABSTRACT

Allogeneic hematopoietic stem cell transplantation (HSCT) with mismatched unrelated donors (MMUD) is associated with inferior outcome compared with matched unrelated donors (MUDs). We aimed to identify permissible mismatches using HLA epitope mismatch algorithm, which determines permissibility by analyzing amino acid sequences, in a single-center cohort of 70 pediatric 9/10 MMUD HSCTs and 157 10/10 MUDs for comparison. Amino acid matching was evaluated for the whole HLA protein, the α-helices, and the β-sheets, in both host vs graft (HvG) and graft vs host (GvH) direction. Superior event-free survival (EFS) was found in 13 patients permissibly mismatched in the HvG direction (totalHvG, 92% vs 58% at 1 year; P = .009) and in 21 patients matched on the α-helices (αHvG, 90% vs 53%; P = .002). These rates were similar to EFS rates in patients with 10/10 MUDs (90% vs 80%; P = .60). EFS was not related to β-sheet amino acid matching, nor to matching in the GvH direction. Overall survival (OS) rates trended similarly to those of EFS for amino acid mismatches (totalHvG, 92% vs 74%; P = .075; αHvG, 90% vs 71%; P = .072). These findings were reproduced in an EBMT Registry inborn errors cohort of 271 pediatric 9/10 MMUD HSCTs and 929 10/10 MUD HSCTs, showing a significant effect of αHvG matching on both OS and EFS and similar OS and EFS between αHvG matched MMUDs and 10/10 MUDs. In summary, HvG amino acid matching on the α-helices identifies 9/10 MMUDs with permissible mismatches, which are correlated with favorable transplant outcomes similar to those of matched donors.

Energy landscapes of peptide-MHC binding

Molecules of the Major Histocompatibility Complex (MHC) present short protein fragments on the cell surface, an important step in T cell immune recognition. MHC-I molecules process peptides from intracellular proteins; MHC-II molecules act in antigen-presenting cells and present peptides derived from extracellular proteins. Here we show that the sequence-dependent energy landscapes of MHC-peptide binding encode class-specific nonlinearities (epistasis). MHC-I has a smooth landscape with global epistasis; the binding energy is a simple deformation of an underlying linear trait. This form of epistasis enhances the discrimination between strong-binding peptides. In contrast, MHC-II has a rugged landscape with idiosyncratic epistasis: binding depends on detailed amino acid combinations at multiple positions of the peptide sequence. The form of epistasis affects the learning of energy landscapes from training data. For MHC-I, a low-complexity problem, we derive a simple matrix model of binding energies that outperforms current models trained by machine learning. For MHC-II, higher complexity prevents learning by simple regression methods. Epistasis also affects the energy and fitness effects of mutations in antigen-derived peptides (epitopes). In MHC-I, large-effect mutations occur predominantly in anchor positions of strong-binding epitopes. In MHC-II, large effects depend on the background epitope sequence but are broadly distributed over the epitope, generating a bigger target for escape mutations due to loss of presentation. Together, our analysis shows how an energy landscape of protein-protein binding constrains the target of escape mutations from T cell immunity, linking the complexity of the molecular interactions to the dynamics of adaptive immune response.

GIHP: Graph convolutional neural network based interpretable pan-specific HLA-peptide binding affinity prediction

Accurately predicting the binding affinities between Human Leukocyte Antigen (HLA) molecules and peptides is a crucial step in understanding the adaptive immune response. This knowledge can have important implications for the development of effective vaccines and the design of targeted immunotherapies. Existing sequence-based methods are insufficient to capture the structure information. Besides, the current methods lack model interpretability, which hinder revealing the key binding amino acids between the two molecules. To address these limitations, we proposed an interpretable graph convolutional neural network (GCNN) based prediction method named GIHP. Considering the size differences between HLA and short peptides, GIHP represent HLA structure as amino acid-level graph while represent peptide SMILE string as atom-level graph. For interpretation, we design a novel visual explanation method, gradient weighted activation mapping (Grad-WAM), for identifying key binding residues. GIHP achieved better prediction accuracy than state-of-the-art methods across various datasets. According to current research findings, key HLA-peptide binding residues mutations directly impact immunotherapy efficacy. Therefore, we verified those highlighted key residues to see whether they can significantly distinguish immunotherapy patient groups. We have verified that the identified functional residues can successfully separate patient survival groups across breast, bladder, and pan-cancer datasets. Results demonstrate that GIHP improves the accuracy and interpretation capabilities of HLA-peptide prediction, and the findings of this study can be used to guide personalized cancer immunotherapy treatment. Codes and datasets are publicly accessible at: https://github.com/sdustSu/GIHP.

H19 influenza A virus exhibits species-specific MHC class II receptor usage

Avian influenza A virus (IAV) surveillance in Northern California, USA, revealed unique IAV hemagglutinin (HA) genome sequences in cloacal swabs from lesser scaups. We found two closely related HA sequences in the same duck species in 2010 and 2013. Phylogenetic analyses suggest that both sequences belong to the recently discovered H19 subtype, which thus far has remained uncharacterized. We demonstrate that H19 does not bind the canonical IAV receptor sialic acid (Sia). Instead, H19 binds to the major histocompatibility complex class II (MHC class II), which facilitates viral entry. Unlike the broad MHC class II specificity of H17 and H18 from bat IAV, H19 exhibits a species-specific MHC class II usage that suggests a limited host range and zoonotic potential. Using cell lines overexpressing MHC class II, we rescued recombinant H19 IAV. We solved the H19 crystal structure and identified residues within the putative Sia receptor binding site (RBS) that impede Sia-dependent entry.

Associations among MHC genes, latitude, and avian malaria infections in the rufous-collared sparrow (Zonotrichia capensis)

The major histocompatibility complex (MHC) is a genetic region in jawed vertebrates that contains key genes involved in the immune response. Associations between the MHC and avian malaria infections in wild birds have been observed and mainly explored in the Northern Hemisphere, while a general lack of information remains in the Southern Hemisphere. Here, we investigated the associations between the MHC genes and infections with Plasmodium and Haemoproteus blood parasites along a latitudinal gradient in South America. We sampled 93 rufous-collared sparrows (Zonotrichia capensis) individuals from four countries, Colombia, Ecuador, Peru, and Chile, and estimated MHC-I and MHC-II allele diversity. We detected between 1-4 (MHC-I) and 1-6 (MHC-II) amino acidic alleles per individual, with signs of positive selection. We obtained generalized additive mixed models to explore the associations between MHC-I and MHC-II diversity and latitude. We also explored the relationship between infection status and latitude/biome. We found a non-linear association between the MHC-II amino acidic allele diversity and latitude. Individuals from north Chile presented a lower MHC genetic diversity than those from other locations. We also found an association between deserts and xeric shrublands and a lower prevalence of Haemoproteus parasites. Our results support a lower MHC genetic in arid or semi-arid habitats in the region with the lower prevalence of Haemoproteus parasites.

The Art of Bioimmunogenomics (BIGs) 5.0 in CAR-T Cell Therapy for Lymphoma Management

Purpose

Lymphoma, the most predominant neoplastic disorder, is divided into Hodgkin and Non-Hodgkin Lymphoma classifications. Immunotherapeutic modalities have emerged as essential methodologies in combating lymphoid malignancies. Chimeric Antigen Receptor (CAR) T cells exhibit promising responses in chemotherapy-resistant B-cell non-Hodgkin lymphoma cases.

Methods

This comprehensive review delineates the advancement of CAR-T cell therapy as an immunotherapeutic instrument, the selection of lymphoma antigens for CAR-T cell targeting, and the conceptualization, synthesis, and deployment of CAR-T cells. Furthermore, it encompasses the advantages and disadvantages of CAR-T cell therapy and the prospective horizons of CAR-T cells from a computational research perspective. In order to improve the design and functionality of artificial CARs, there is a need for TCR recognition investigation, followed by the implementation of a quality surveillance methodology.

Results

Various lymphoma antigens are amenable to CAR-T cell targeting, such as CD19, CD20, CD22, CD30, the kappa light chain, and ROR1. A notable merit of CAR-T cell therapy is the augmentation of the immune system's capacity to generate tumoricidal activity in patients exhibiting chemotherapy-resistant lymphoma. Nevertheless, it also introduces manufacturing impediments that are laborious, technologically demanding, and financially burdensome. Physical, physicochemical, and physiological limitations further exacerbate the challenge of treating solid neoplasms with CAR-T cells.

Conclusion

While the efficacy and safety of CAR-T cell immunotherapy remain subjects of fervent investigation, the promise of this cutting-edge technology offers valuable insights for the future evolution of lymphoma treatment management approaches. Moreover, CAR-T cell therapies potentially benefit patients, motivating regulatory bodies to foster international collaboration.

Diversity of major histocompatibility complex of II B gene and mate choice in a monogamous and long-lived seabird, the Little Auk (Alle alle)

The major histocompatibility complex (MHC) plays a key role in the adaptive immune system of vertebrates, and is known to influence mate choice in many species. In birds, the MHC has been extensively examined but mainly in galliforms and passerines while other taxa that represent specific ecological and evolutionary life-histories, like seabirds, are underexamined. Here, we characterized diversity of MHC Class II B exon 2 in a colonial pelagic seabird, the Little Auk (or Dovekie Alle alle). We further examined whether MHC variation could be maintained through balancing selection and disassortative mating. We found high polymorphism at the genotyped MHC fragment, characterizing 99 distinct alleles across 140 individuals from three populations. The alleles frequencies exhibited a similar skewed distribution in both sexes, with the four most commonly occurring alleles representing approximately 35% of allelic variation. The results of a Bayesian site-by-site selection analysis suggest evidence of balancing selection and no direct evidence for MHC-dependent disassortative mating preferences in the Little Auk. The latter result might be attributed to the high overall polymorphism of the examined fragment, which itself may be maintained by the large population size of the species.

Experimental Structures of Antibody/MHC-I Complexes Reveal Details of Epitopes Overlooked by Computational Prediction

mAbs to MHC class I (MHC-I) molecules have proved to be crucial reagents for tissue typing and fundamental studies of immune recognition. To augment our understanding of epitopic sites seen by a set of anti-MHC-I mAb, we determined X-ray crystal structures of four complexes of anti-MHC-I Fabs bound to peptide/MHC-I/β2-microglobulin (pMHC-I). An anti-H2-Dd mAb, two anti-MHC-I α3 domain mAbs, and an anti-β2-microglobulin mAb bind pMHC-I at sites consistent with earlier mutational and functional experiments, and the structures explain allelomorph specificity. Comparison of the experimentally determined structures with computationally derived models using AlphaFold Multimer showed that although predictions of the individual pMHC-I heterodimers were quite acceptable, the computational models failed to properly identify the docking sites of the mAb on pMHC-I. The experimental and predicted structures provide insight into strengths and weaknesses of purely computational approaches and suggest areas that merit additional attention.

Cancer neoepitopes viewed through negative selection and peripheral tolerance: a new path to cancer vaccines

A proportion of somatic mutations in tumors create neoepitopes that can prime T cell responses that target the MHC I-neoepitope complexes on tumor cells, mediating tumor control or rejection. Despite the compelling centrality of neoepitopes to cancer immunity, we know remarkably little about what constitutes a neoepitope that can mediate tumor control in vivo and what distinguishes such a neoepitope from the vast majority of similar candidate neoepitopes that are inefficacious in vivo. Studies in mice as well as clinical trials have begun to reveal the unexpected paradoxes in this area. Because cancer neoepitopes straddle that ambiguous ground between self and non-self, some rules that are fundamental to immunology of frankly non-self antigens, such as viral or model antigens, do not appear to apply to neoepitopes. Because neoepitopes are so similar to self-epitopes, with only small changes that render them non-self, immune response to them is regulated at least partially the way immune response to self is regulated. Therefore, neoepitopes are viewed and understood here through the clarifying lens of negative thymic selection. Here, the emergent questions in the biology and clinical applications of neoepitopes are discussed critically and a mechanistic and testable framework that explains the complexity and translational potential of these wonderful antigens is proposed.

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.

PANDORA v2.0: Benchmarking peptide-MHC II models and software improvements

T-cell specificity to differentiate between self and non-self relies on T-cell receptor (TCR) recognition of peptides presented by the Major Histocompatibility Complex (MHC). Investigations into the three-dimensional (3D) structures of peptide:MHC (pMHC) complexes have provided valuable insights of MHC functions. Given the limited availability of experimental pMHC structures and considerable diversity of peptides and MHC alleles, it calls for the development of efficient and reliable computational approaches for modeling pMHC structures. Here we present an update of PANDORA and the systematic evaluation of its performance in modelling 3D structures of pMHC class II complexes (pMHC-II), which play a key role in the cancer immune response. PANDORA is a modelling software that can build low-energy models in a few minutes by restraining peptide residues inside the MHC-II binding groove. We benchmarked PANDORA on 136 experimentally determined pMHC-II structures covering 44 unique αβ chain pairs. Our pipeline achieves a median backbone Ligand-Root Mean Squared Deviation (L-RMSD) of 0.42 Å on the binding core and 0.88 Å on the whole peptide for the benchmark dataset. We incorporated software improvements to make PANDORA a pan-allele framework and improved the user interface and software quality. Its computational efficiency allows enriching the wealth of pMHC binding affinity and mass spectrometry data with 3D models. These models can be used as a starting point for molecular dynamics simulations or structure-boosted deep learning algorithms to identify MHC-binding peptides. PANDORA is available as a Python package through Conda or as a source installation at https://github.com/X-lab-3D/PANDORA.

Experimental structures of antibody/MHC-I complexes reveal details of epitopes overlooked by computational prediction

Monoclonal antibodies (mAb) to major histocompatibility complex class I (MHC-I) molecules have proved to be crucial reagents for tissue typing and fundamental studies of immune recognition. To augment our understanding of epitopic sites seen by a set of anti-MHC-I mAb, we determined X-ray crystal structures of four complexes of anti-MHC-I antigen-binding fragments (Fab) bound to peptide/MHC-I/β2m (pMHC-I). An anti-H2-Dd mAb, two anti-MHC-I α3 domain mAb, and an anti-β2-microglobulin (β2m) mAb bind pMHC-I at sites consistent with earlier mutational and functional experiments, and the structures explain allelomorph specificity. Comparison of the experimentally determined structures with computationally derived models using AlphaFold Multimer (AF-M) showed that although predictions of the individual pMHC-I heterodimers were quite acceptable, the computational models failed to properly identify the docking sites of the mAb on pMHC-I. The experimental and predicted structures provide insight into strengths and weaknesses of purely computational approaches and suggest areas that merit additional attention.

Development of a semi-automated MHC-associated peptide proteomics (MAPPs) method using streptavidin bead-based immunoaffinity capture and nano LC-MS/MS to support immunogenicity risk assessment in drug development

Major histocompatibility complex (MHC)-Associated Peptide Proteomics (MAPPs) is an ex vivo method used to assess the immunogenicity risk of biotherapeutics. MAPPs can identify potential T-cell epitopes within the biotherapeutic molecule. Using adalimumab treated human monocyte derived dendritic cells (DCs) and a pan anti-HLA-DR antibody (Ab), we systematically automated and optimized biotin/streptavidin (SA)-capture antibody coupling, lysate incubation with capture antibody, as well as the washing and elution steps of a MAPPs method using functionalized magnetic beads and a KingFisher Magnetic Particle processor. Automation of these steps, combined with capturing using biotinylated-Ab/SA magnetic beads rather than covalently bound antibody, improved reproducibility as measured by minimal inter-and intra-day variability, as well as minimal analyst-to-analyst variability. The semi-automated MAPPs workflow improved sensitivity, allowing for a lower number of cells per analysis. The method was assessed using five different biotherapeutics with varying immunogenicity rates ranging from 0.1 to 48% ADA incidence in the clinic. Biotherapeutics with ≥10%immunogenicity incidence consistently presented more peptides (1.8-28 fold) and clusters (10-21 fold) compared to those with <10% immunogenicity incidence. Our semi-automated MAPPs method provided two main advantages over a manual workflow- the robustness and reproducibility affords confidence in the epitopes identified from as few as 5 to 10 donors and the method workflow can be readily adapted to incorporate different capture Abs in addition to anti-HLA-DR. The incorporation of semi-automated MAPPs with biotinylated-Ab/SA bead-based capture in immunogenicity screening strategies allows the generation of more consistent and reliable data, helping to improve immunogenicity prediction capabilities in drug development. MHC associated peptide proteomics (MAPPs), Immunogenicity risk assessment, in vitro/ex vivo, biotherapeutics, Major Histocompatibility Complex Class II (MHC II), LC-MS, Immunoaffinity Capture, streptavidin magnetic beads.

Diversity and evolution of the MHC class II DRB gene in the Capra sibirica experienced a demographic fluctuation in China

The major histocompatibility complex (MHC) genes are the most polymorphic genes in vertebrates, and their proteins play a critical role in adaptive immunity for defense against a variety of pathogens. MHC diversity was lost in many species after experiencing a decline in size. To understand the variation and evolution of MHC genes in the Siberian ibex, Capra sibirica, which has undergone a population decline, we analyzed the variation of the second exon of MHC class II DRB genes in samples collected from five geographic localities in Xinjiang, China, that belong to three diverged mitochondrial clades. Consequently, we identified a total of 26 putative functional alleles (PFAs) with 260 bp in length from 43 individuals, and found one (for 27 individuals) to three (for 5 individuals) PFAs per individual, indicating the presence of one or two DRB loci per haploid genome. The Casi-DRB1*16 was the most frequently occurring PFA, Casi-DRB1*22 was found in only seven individuals, 14 PFAs occurred once, 7 PFAs twice, implying high frequency of rare PFAs. Interestingly, more than half (15) of the PFAs were specific to clade I, only two and three PFAs were specific to clades II and III, respectively. So, we assume that the polygamy and sexual segregation nature of this species likely contributed to the allelic diversity of DRB genes. Genetic diversity indices showed that PFAs of clade II were lower in nucleotide, amino acid, and supertype diversity compared to those of the other two clades. The pattern of allele sharing and FST values between the three clades was to some extent in agreement with the pattern observed in mitochondrial DNA divergence. In addition, recombination analyses revealed no evidence for significant signatures of recombination events. Alleles shared by clades III and the other two clades diverged 6 million years ago, and systematic neighbor grids showed Trans-species polymorphism. Together with the PAML and MEME analyses, the results indicated that the DRB gene in C. sibirica evolved under balancing and positive selection. However, by comparison, it can be clearly seen that different populations were under different selective pressures. Our results are valuable in understanding the diversity and evolution of the DRB gene in a mountain living C. sibirica and in making decisions on future long-term protection strategies.

Monomeric agonist peptide/MHCII complexes activate T-cells in an autonomous fashion

Molecular crowding of agonist peptide/MHC class II complexes (pMHCIIs) with structurally similar, yet per se non-stimulatory endogenous pMHCIIs is postulated to sensitize T-cells for the recognition of single antigens on the surface of dendritic cells and B-cells. When testing this premise with the use of advanced live cell microscopy, we observe pMHCIIs as monomeric, randomly distributed entities diffusing rapidly after entering the APC surface. Synaptic TCR engagement of highly abundant endogenous pMHCIIs is low or non-existent and affects neither TCR engagement of rare agonist pMHCII in early and advanced synapses nor agonist-induced TCR-proximal signaling. Our findings highlight the capacity of single freely diffusing agonist pMHCIIs to elicit the full T-cell response in an autonomous and peptide-specific fashion with consequences for adaptive immunity and immunotherapeutic approaches.

Ancient duplication, coevolution, and selection at the MHC class IIA and IIB genes of birds

Introduction

The Major Histocompatibility Complex (MHC) of vertebrates is a dynamically evolving multigene family primarily responsible for recognizing non-self peptide antigens and triggering a pathogen-specific adaptive immune response. In birds, the MHC was previously thought to evolve via concerted evolution with high degree of gene homogenization and the rapid loss of orthology. However, the discovery of two ancient avian MHC-IIB gene lineages (DAB1 and DAB2) originating before the radiation of extant birds indicated that despite the action of concerted evolution, orthology may be detectable for long evolutionary periods.

Methods

Here, we take advantage of rapidly accumulating digital genomic resources to search for the signal of an ancient duplication at the avian MHC-IIA genes, as well as to compare phylogenetic distribution and selection between MHC-IIA and IIB gene lineages.

Results

The analysis of MHC sequences from over 230 species representing ca. 70 bird families revealed the presence of two ancient MHC-IIA gene lineages (DAA1 and DAA2) and showed that their phylogenetic distribution matches exactly the distribution of DAB1 and DAB2 lineages, suggesting tight coevolution. The early post-duplication divergence of DAA1 and DAA2 was driven by positive selection fixing radical amino acid differences within the membrane-proximal domain and, most probably, being functionally related to the interactions between α2 and β2 chains of the MHC-II heterodimer. We detected no evidence for an overall (gene-wide) relaxation or intensification of selection at either DAA1/DAB1 or DAA2/DAB2, but codon-specific differences in selection signature were found at the peptide-binding sites between the two gene lineages, perhaps implying specialization to different pathogen regimes.

Discussion

Our results suggest that specific pairing of MHC-II α and β chains may have an adaptive significance, a conclusion that advances knowledge on the macroevolution of the avian MHC-II and opens exciting novel directions for future research.

Conserved biophysical compatibility among the highly variable germline-encoded regions shapes TCR-MHC interactions

T cells are critically important components of the adaptive immune system primarily responsible for identifying and responding to pathogenic challenges. This recognition of pathogens is driven by the interaction between membrane-bound T cell receptors (TCRs) and antigenic peptides presented on major histocompatibility complex (MHC) molecules. The formation of the TCR-peptide-MHC complex (TCR-pMHC) involves interactions among germline-encoded and hypervariable amino acids. Germline-encoded and hypervariable regions can form contacts critical for complex formation, but only interactions between germline-encoded contacts are likely to be shared across many of all the possible productive TCR-pMHC complexes. Despite this, experimental investigation of these interactions have focused on only a small fraction of the possible interaction space. To address this, we analyzed every possible germline-encoded TCR-MHC contact in humans, thereby generating the first comprehensive characterization of these largely antigen-independent interactions. Our computational analysis suggests that germline-encoded TCR-MHC interactions that are conserved at the sequence level are rare due to the high amino acid diversity of the TCR CDR1 and CDR2 loops, and that such conservation is unlikely to dominate the dynamic protein-protein binding interface. Instead, we propose that binding properties such as the docking orientation are defined by regions of biophysical compatibility between these loops and the MHC surface.

CD1 lipidomes reveal lipid-binding motifs and size-based antigen-display mechanisms

The CD1 system binds lipid antigens for display to T cells. Here, we solved lipidomes for the four human CD1 antigen-presenting molecules, providing a map of self-lipid display. Answering a basic question, the detection of >2,000 CD1-lipid complexes demonstrates broad presentation of self-sphingolipids and phospholipids. Whereas peptide antigens are chemically processed, many lipids are presented in an unaltered form. However, each type of CD1 protein differentially edits the self-lipidome to show distinct capture motifs based on lipid length and chemical composition, suggesting general antigen display mechanisms. For CD1a and CD1d, lipid size matches the CD1 cleft volume. CD1c cleft size is more variable, and CD1b is the outlier, where ligands and clefts show an extreme size mismatch that is explained by uniformly seating two small lipids in one cleft. Furthermore, the list of compounds that comprise the integrated CD1 lipidome supports the ongoing discovery of lipid blockers and antigens for T cells.

Genetic variation at innate and adaptive immune genes - contrasting patterns of differentiation and local adaptation in a wild gull

Immunogenetic variation in natural vertebrate populations is expected to respond to spatial and temporal fluctuations in pathogen assemblages. While spatial heterogeneity in pathogen-driven selection enhances local immunogenetic adaptations and population divergence, different immune genes may yield contrasting responses to the environment. Here, we investigated population differentiation at the key pathogen recognition genes of the innate and adaptive immune system in a colonial bird species, the black-headed gull Chroicocephalus ridibundus. We assessed genetic variation at three toll-like receptor (TLR) genes (innate immunity) and the major histocompatibility complex (MHC) class I and II genes (adaptive immunity) in gulls from seven colonies scattered across Poland. As expected, we found much greater polymorphism at the MHC than TLRs. Population differentiation at the MHC class II, but not MHC-I, was significantly stronger than at neutral microsatellite loci, suggesting local adaptation. This could reflect spatial variation in the composition of extracellular parasite communities (e.g., helminths), possibly driven by sharp differences in habitat structure between colonies. Despite contrasting patterns of population differentiation, both MHC classes showed similar regimes of diversifying selection. Some significant population differentiation was also observed at TLRs, suggesting that innate immune receptors may respond to fine-scale spatial variation in pathogen pressure, although this pattern could have been enhanced by drift. Our results suggested that local adaptation at the pathogen recognition immune genes can be maintained at relatively small or moderate spatial scales in species with high dispersal potential and they highlighted the complexity of immunogenetic responses of animals to heterogeneous environments.

Pathogenic selection promotes adaptive immune variations against serious bottlenecks in early invasions of bullfrogs

Adaptive genetic variations are key for understanding evolutionary processes influencing invasions. However, we have limited knowledge on how adaptive genetic diversity in invasive species responds to new pathogenic environments. Here, we compared variations in immune major histocompatibility complex (MHC) class-II β gene and neutral loci in relation to pathogenic chytrid fungus (Batrachochytrium dendrobatidis, Bd) infection across invasive and native populations of American bullfrog between China and United States (US). Chinese invasive populations show a 60% reduction in neutral cytb variations relative to US native populations, and there were similar MHC variation and functional diversity between them. One MHC allele private to China was under recent positive selection and associated with decreased Bd infection, partly explaining the lower Bd prevalence for Chinese populations than for native US populations. These results suggest that pathogen-mediated selection favors adaptive MHC variations and functional diversity maintenance against serious bottlenecks during the early invasions (within 15 generations) of bullfrogs.

A survey of the adaptive immune genes of the polka-dot batfish Ogcocephalus cubifrons

BACKGROUND

The anglerfish, belonging to the teleost order Lophiiformes, are a diverse and species-rich group of fish that are known to exhibit a number of unique morphological, reproductive and immunological adaptations. Work to date has identified the loss of specific adaptive immune components in two of the five Lophiiformes sub-orders (Lophioidei and Ceratioidei), while no anomalies have been identified to date in two other sub-orders, Antennaroidei and Chaunacoidei. The immunogenome of the fifth sub-order, Ogcocephaloidei has not yet been investigated, and we have therefore used whole genome shotgun sequencing, combined with RNA-seq, to survey the adaptive immune capabilities of the polka-dot batfish, O. cubifrons, as a representative of this as yet unexplored sub-order.

RESULTS

We find that the O. cubifrons genome encodes the core genes needed to mount adaptive T and B cell responses. These genes include those necessary for rearranging and editing antigen receptors, the antigen receptors themselves; as well as the co-receptors, signalling molecules, and antigen presenting molecules (both class I and class II) needed for B cell and T cell development and activation.

CONCLUSIONS

From an immune perspective, the polka-dot batfish has a canonical complement of adaptive immune genes, and does not exhibit any of the adaptive immune changes previously identified in monkfish and oceanic anglerfish.

HLA-II immunopeptidome profiling and deep learning reveal features of antigenicity to inform antigen discovery

CD4+ T cell responses are exquisitely antigen specific and directed toward peptide epitopes displayed by human leukocyte antigen class II (HLA-II) on antigen-presenting cells. Underrepresentation of diverse alleles in ligand databases and an incomplete understanding of factors affecting antigen presentation in vivo have limited progress in defining principles of peptide immunogenicity. Here, we employed monoallelic immunopeptidomics to identify 358,024 HLA-II binders, with a particular focus on HLA-DQ and HLA-DP. We uncovered peptide-binding patterns across a spectrum of binding affinities and enrichment of structural antigen features. These aspects underpinned the development of context-aware predictor of T cell antigens (CAPTAn), a deep learning model that predicts peptide antigens based on their affinity to HLA-II and full sequence of their source proteins. CAPTAn was instrumental in discovering prevalent T cell epitopes from bacteria in the human microbiome and a pan-variant epitope from SARS-CoV-2. Together CAPTAn and associated datasets present a resource for antigen discovery and the unraveling genetic associations of HLA alleles with immunopathologies.

Evolutionarily conserved amino acids in MHC-II mediate bat influenza A virus entry into human cells

The viral hemagglutinins of conventional influenza A viruses (IAVs) bind to sialylated glycans on host cell surfaces for attachment and subsequent infection. In contrast, hemagglutinins of bat-derived IAVs target major histocompatibility complex class II (MHC-II) for cell entry. MHC-II proteins from various vertebrate species can facilitate infection with the bat IAV H18N11. Yet, it has been difficult to biochemically determine the H18:MHC-II binding. Here, we followed a different approach and generated MHC-II chimeras from the human leukocyte antigen DR (HLA-DR), which supports H18-mediated entry, and the nonclassical MHC-II molecule HLA-DM, which does not. In this context, viral entry was supported only by a chimera containing the HLA-DR α1, α2, and β1 domains. Subsequent modeling of the H18:HLA-DR interaction identified the α2 domain as central for this interaction. Further mutational analyses revealed highly conserved amino acids within loop 4 (N149) and β-sheet 6 (V190) of the α2 domain as critical for virus entry. This suggests that conserved residues in the α1, α2, and β1 domains of MHC-II mediate H18-binding and virus propagation. The conservation of MHC-II amino acids, which are critical for H18N11 binding, may explain the broad species specificity of this virus.

Antigen presentation in vertebrates: Structural and functional aspects

Antigen presentation is a key process of the immune system and is responsible for the activation of T cells. The main characters are the major histocompatibility complex class I (MHC-I) and class II (MHC-II) molecules, and accessory proteins that act as chaperones for these glycoproteins. Current knowledge of this process and also the elucidation of the structural features of these proteins, has been extensively reviewed in humans. Unfortunately, this is not the case for non-human species, wherein the function and structural characteristic of the antigen presentation proteins is far from being understood. The majority of previous studies in non-human species, especially in teleost fish and lower vertebrates, are limited to the transcriptomic level, which leads to gaps in the knowledge about the functional process of antigen presentation in these species. This review summarizes what is known so far about antigen presentation pathways in vertebrates from a structural and functional perspective. The focus is not only on the MHC receptors, but also, on the forgotten characters of these pathways such as the proteins of the peptide loading complex, and the MHC-II chaperone invariant chain.

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.

Evolution of MHC IIB Diversity Across Cichlid Fish Radiations

The genes of the major histocompatibility complex (MHC) are among the most polymorphic genes in vertebrates and crucial for their adaptive immune response. These genes frequently show inconsistencies between allelic genealogies and species phylogenies. This phenomenon is thought to be the result of parasite-mediated balancing selection maintaining ancient alleles through speciation events (trans-species polymorphism [TSP]). However, allele similarities may also arise from postspeciation mechanisms, such as convergence or introgression. Here, we investigated the evolution of MHC class IIB diversity in the cichlid fish radiations across Africa and the Neotropics by a comprehensive review of available MHC IIB DNA sequence information. We explored what mechanism explains the MHC allele similarities found among cichlid radiations. Our results showed extensive allele similarity among cichlid fish across continents, likely due to TSP. Functionality at MHC was also shared among species of the different continents. The maintenance of MHC alleles for long evolutionary times and their shared functionality may imply that certain MHC variants are essential in immune adaptation, even in species that diverged millions of years ago and occupy different environments.

Loading dynamics of one SARS-CoV-2-derived peptide into MHC-II revealed by kinetic models

Major histocompatibility complex class II (MHC-II) plays an indispensable role in activating CD4+ T cell immune responses by presenting antigenic peptides on the cell surface for recognition by T cell receptors. The assembly of MHC-II and antigenic peptide is therefore a prerequisite for the antigen presentation. To date, however, the atomic-level mechanism underlying the peptide-loading dynamics for MHC-II is still elusive. Here, by constructing Markov state models based on extensive all-atom molecular dynamics simulations, we reveal the complete peptide-loading dynamics into MHC-II for one SARS-CoV-2 S-protein-derived antigenic peptide (235ITRFQTLLALHRSYL249). Our Markov state model identifies six metastable states (S1-S6) during the peptide-loading process and determines two dominant loading pathways. The peptide could potentially approach the antigen-binding groove via either its N- or C-terminus. Then, the consecutive insertion of several anchor residues into the binding pockets profoundly dictates the peptide-loading dynamics. Notably, the MHC-II αA52-E55 motif could guide the peptide loading into the antigen-binding groove via forming β-sheets conformation with the incoming peptide. The rate-limiting step, namely S5→S6, is mainly attributed to a considerable desolvation penalty triggered by the binding of the peptide C-terminus. Moreover, we further examined the conformational changes associated with the peptide exchange process catalyzed by the chaperon protein HLA-DM. A flipped-out conformation of MHC-II αW43 captured in S1-S3 is considered a critical anchor point for HLA-DM to modulate the structural dynamics. Our work provides deep structural insights into the key regulatory factors in MHC-II responsible for peptide recognition and guides future design for peptide vaccines against SARS-CoV-2.

Characterization and genetic diversity of MHC class II DRB genes in the Arabian camel (Camelus dromedarius)

This study investigated the MHC DRB genes in the Arabian camel (Camelus dromedarius). The results revealed the presence of - at least - two transcribed DRB-like genes in chromosome 20, designated MhcCadr-DRB1 and MhcCadr-DRB2. These genes are 155 Kb apart, have similar gene structure, and are transcribed in opposite directions. Compared to DRB1, the DRB2 locus contains a deletion of 12 nucleotides in the second exon (270 bp), exhibits lower transcript abundance, and is expressed as two splice variants differing by exon 2 skipping. This gene seems to be of minor functional relevance in the dromedary camel. Conversely, the DRB1 is thought to be the main gene in this species showing higher transcript abundance and polymorphism levels. A total of seven DRB1 exon 2 alleles were identified in the Tunisian dromedary camel resulting from 18 amino acid substitutions. Six full length alleles were characterized at the mRNA level. Although there is no clear evidence for balancing selection (i.e., heterozygote advantage), signals of weak historical positive selection acting on the DRB1 gene were detected, as indicated by the limited number of the sites being positively selected. This trend might be related to the low exposure to pathogens and to the demographic history of the species. Comparative analysis with Bactrian and wild camel genomes suggested occurrence of trans species polymorphism (TSP) in the Camelus genus. The results lay the foundation for the MHC DRB1 genetic diversity analysis in this genus since the developed genotyping protocols are fully applicable in the three Camelus species.

Extensive MHC class IIβ diversity across multiple loci in the small-spotted catshark (Scyliorhinus canicula)

The major histocompatibility complex (MHC) is a multigene family responsible for pathogen detection, and initiation of adaptive immune responses. Duplication, natural selection, recombination, and their resulting high functional genetic diversity spread across several duplicated loci are the main hallmarks of the MHC. Although these features were described in several jawed vertebrate lineages, a detailed MHC IIβ characterization at the population level is still lacking for chondrichthyans (chimaeras, rays and sharks), i.e. the most basal lineage to possess an MHC-based adaptive immune system. We used the small-spotted catshark (Scyliorhinus canicula, Carcharhiniformes) as a case-study species to characterize MHC IIβ diversity using complementary molecular tools, including publicly available genome and transcriptome datasets, and a newly developed high-throughput Illumina sequencing protocol. We identified three MHC IIβ loci within the same genomic region, all of which are expressed in different tissues. Genetic screening of the exon 2 in 41 individuals of S. canicula from a single population revealed high levels of sequence diversity, evidence for positive selection, and footprints of recombination. Moreover, the results also suggest the presence of copy number variation in MHC IIβ genes. Thus, the small-spotted catshark exhibits characteristics of functional MHC IIβ genes typically observed in other jawed vertebrates.

Meta-analysis of major histocompatibility complex (MHC) class IIA reveals polymorphism and positive selection in many vertebrate species

Pathogen-mediated selection and sexual selection are important drivers of evolution. Both processes are known to target genes of the major histocompatibility complex (MHC), a gene family encoding cell-surface proteins that display pathogen peptides to the immune system. The MHC is also a model for understanding processes such as gene duplication and trans-species allele sharing. The class II MHC protein is a heterodimer whose peptide-binding groove is encoded by an MHC-IIA gene and an MHC-IIB gene. However, our literature review found that class II MHC papers on infectious disease or sexual selection included IIA data only 18% and 9% of the time, respectively. To assess whether greater emphasis on MHC-IIA is warranted, we analysed MHC-IIA sequence data from 50 species of vertebrates (fish, amphibians, birds, mammals) to test for polymorphism and positive selection. We found that the number of MHC-IIA alleles within a species was often high, and covaried with sample size and number of MHC-IIA genes assayed. While MHC-IIA variability tended to be lower than that of MHC-IIB, the difference was only ~25%, with ~3 fewer IIA alleles than IIB. Furthermore, the unexpectedly high MHC-IIA variability showed clear signatures of positive selection in most species, and positive selection on MHC-IIA was stronger in fish than in other surveyed vertebrate groups. Our findings underscore that MHC-IIA can be an important target of selection. Future studies should therefore expand the characterization of MHC-IIA at both allelic and genomic scales, and incorporate MHC-IIA into models of fitness consequences of MHC variation.

Intronic primers reveal unexpectedly high major histocompatibility complex diversity in Antarctic fur seals

The major histocompatibility complex (MHC) is a group of genes comprising one of the most important components of the vertebrate immune system. Consequently, there has been much interest in characterising MHC variation and its relationship with fitness in a variety of species. Due to the exceptional polymorphism of MHC genes, careful PCR primer design is crucial for capturing all of the allelic variation present in a given species. We therefore developed intronic primers to amplify the full-length 267 bp protein-coding sequence of the MHC class II DQB exon 2 in the Antarctic fur seal. We then characterised patterns of MHC variation among mother-offspring pairs from two breeding colonies and detected 19 alleles among 771 clone sequences from 56 individuals. The distribution of alleles within and among individuals was consistent with a single-copy, classical DQB locus showing Mendelian inheritance. Amino acid similarity at the MHC was significantly associated with genome-wide relatedness, but no relationship was found between MHC heterozygosity and genome-wide heterozygosity. Finally, allelic diversity was several times higher than reported by a previous study based on partial exon sequences. This difference appears to be related to allele-specific amplification bias, implying that primer design can strongly impact the inference of MHC diversity.

Salamanders reveal novel trajectories of amphibian MHC evolution

Genes of the major histocompatibility complex (MHC) code for immune proteins that are crucial for pathogen recognition in vertebrates. MHC research in nonmodel taxa has long been hampered by its genomic complexity that makes the locus-specific genotyping challenging. The recent progress in sequencing and genotyping methodologies allows an extensive phylogenetic coverage in studies of MHC evolution. Here, we analyzed the peptide-binding region of MHC class I (MHC-I) in 30 species of salamanders from six families representative of Urodela phylogeny. This extensive dataset revealed an extreme diversity of MHC-I in salamanders, both in terms of sequence diversity (about 3000 variants) and architecture (2-22 gene copies per species). The signal of positive selection was moderate and consistent between both peptide-binding domains, but varied greatly between genera. Positions of positively selected sites mostly coincided with human peptide-binding sites, suggesting similar structural properties of MHC-I molecules across distant vertebrate lineages. Finally, we provided evidence for the common intraexonic recombination at MHC-I and for the role of life history traits in the processes of MHC-I expansion/contraction. Our study revealed novel evolutionary trajectories of amphibian MHC and it contributes to the understanding of the mechanisms that generated extraordinary MHC diversity throughout vertebrate evolution.

Habitat fragmentation in the Brazilian Atlantic Forest is associated with erosion of frog immunogenetic diversity and increased fungal infections

Habitat fragmentation and infectious diseases threaten wildlife globally, but the interactions of these threats are poorly understood. For instance, while habitat fragmentation can impact genetic diversity at neutral loci, the impacts on disease-relevant loci are less well-studied. We examined the effects of habitat fragmentation in Brazil's Atlantic Forest on amphibian genetic diversity at an immune locus related to antigen presentation and detection (MHC IIB Exon 2). We used a custom high-throughput assay to sequence a fragment of MHC IIB and quantified Batrachochytrium dendrobatidis (Bd) infections in six frog species in two Atlantic Forest regions. Habitat fragmentation was associated with genetic erosion at MHC IIB Exon 2. This erosion was most severe in forest specialists. Significant Bd infections were detected only in one Atlantic Forest region, potentially due to relatively higher elevation. In this region, forest specialists showed an increase in both Bd prevalence and infection loads in fragmented habitats. Reduced population-level MHC IIB diversity was associated with increased Bd infection risk. On the individual level, MHC IIB heterozygotes exhibited a trend toward reduced Bd infection risk, although this was marginally non-significant. Our results suggest that habitat fragmentation increases Bd infection susceptibility in amphibians, mediated at least in part through erosion of immunogenetic diversity. Our findings have implications for management of fragmented populations in the face of emerging infectious diseases.

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.

Novel insights on aquatic mammal MHC evolution: Evidence from manatee DQB diversity

The low diversity in marine mammal major histocompatibility complex (MHC) appears to support the hypothesis of reduced pathogen selective pressure in aquatic systems compared to terrestrial environments. However, the lack of characterization of the aquatic and evolutionarily distant Sirenia precludes drawing more generalized conclusions. Therefore, we aimed to characterize the MHC DQB diversity of two manatee species and compare it with those reported for marine mammals. Our results identified 12 and 6 alleles in T. inunguis and T. manatus, respectively. Alleles show high rates of nonsynonymous substitutions, suggesting loci are evolving under positive selection. Among aquatic mammals, Pinnipeda DQB had smaller numbers of alleles, higher synonymous substitution rate, and a dN/dS ratio closer to 1, suggesting it may be evolving under more relaxed selection compared to fully aquatic mammals. This contradicts one of the predictions of the hypothesis that aquatic environments impose reduced pathogen pressure to mammalian immune system. These results suggest that the unique evolutionary trajectories of mammalian MHC may impose challenges in drawing ecoevolutionary conclusions from comparisons across distant vertebrate lineages.

Intra-specific copy number variation of MHC class II genes in the Siamese fighting fish

Duplicates of genes for major histocompatibility complex (MHC) molecules can be subjected to selection independently and vary markedly in their evolutionary rates, sequence polymorphism, and functional roles. Therefore, without a thorough understanding of their copy number variation (CNV) in the genome, the MHC-dependent fitness consequences within a species could be misinterpreted. Studying the intra-specific CNV of this highly polymorphic gene, however, has long been hindered by the difficulties in assigning alleles to loci and the lack of high-quality genomic data. Here, using the high-quality genome of the Siamese fighting fish (Betta splendens), a model for mate choice studies, and the whole-genome sequencing (WGS) data of 17 Betta species, we achieved locus-specific amplification of their three classical MHC class II genes - DAB1, DAB2, and DAB3. By performing quantitative PCR and depth-of-coverage analysis using the WGS data, we revealed intra-specific CNV at the DAB3 locus. We identified individuals that had two allelic copies (i.e., heterozygous or homozygous) or one allele (i.e., hemizygous) and individuals without this gene. The CNV was due to the deletion of a 20-kb-long genomic region harboring both the DAA3 and DAB3 genes. We further showed that the three DAB genes were under different modes of selection, which also applies to their corresponding DAA genes that share similar pattern of polymorphism. Our study demonstrates a combined approach to study CNV within a species, which is crucial for the understanding of multigene family evolution and the fitness consequences of CNV.

Molecular characteristics, polymorphism and expression analysis of mhc Ⅱ in yellow catfish(pelteobagrus fulvidraco)responding to Flavobacterium columnare infection

The major histocompatibility complex (MHC) is an important component of the immune system of vertebrates, which plays a vital role in presenting extrinsic antigens. In this study, we cloned and characterized the mhc ⅡA and mhc ⅡB genes of yellow catfish Pelteobagrus fulvidraco. The open reading frames (ORFs) of mhc ⅡA and mhc ⅡB genes were 708 bp and 747bp in length, encoding 235 and 248 amino acids, respectively. The structure of mhc ⅡA and mhc ⅡB includes a signal peptide, an α1/β1 domain, an α2/β2 domain, a transmembrane region and a cytoplasmic region. Homologous identity analysis revealed that both mhc ⅡA and mhc ⅡB shared high protein sequence similarity with that of Chinese longsnout catfish Leiocassis longirostris. mhc ⅡA and mhc ⅡB showed similar expression patterns in different tissues, with the higher expression level in spleen, head kidney and gill and lower expression in liver, stomach, gall bladder and heart. The mRNA expression level of mhc ⅡA and mhc ⅡB in different embryonic development stages also showed the similar trends. The higher expression was detected from fertilized egg to 32 cell stage, low expression from multicellular period to 3 days post hatching (dph), and then the expression increased to a higher level from 4 dph to 14 dph. The mRNA expression levels of mhc ⅡA and mhc ⅡB were significantly up-regulated not only in the body kidney and spleen, but also in the midgut, hindgut, liver and gill after challenge of Flavobacterium columnare. The results suggest that Mhc Ⅱ plays an important role in the anti-infection process of yellow catfish P. fulvidraco.

Characterization, Selection, and Trans-Species Polymorphism in the MHC Class II of Heermann’s Gull (Charadriiformes)

The major histocompatibility complex (MHC) enables vertebrates to cope with pathogens and maintain healthy populations, thus making it a unique set of loci for addressing ecology and evolutionary biology questions. The aim of our study was to examine the variability of Heermann’s Gull MHC class II (MHCIIB) and compare these loci with other Charadriiformes. Fifty-nine MHCIIB haplotypes were recovered from sixty-eight Heermann’s Gulls by cloning, of them, twelve were identified as putative true alleles, forty-five as unique alleles, and two as pseudogenes. Intra and interspecific relationships indicated at least two loci in Heermann’s Gull MHCIIB and trans-species polymorphism among Charadriiformes (coinciding with the documented evidence of two ancient avian MHCIIB lineages, except in the Charadriidae family). Additionally, sites under diversifying selection revealed a better match with peptide-binding sites inferred in birds than those described in humans. Despite the negative anthropogenic activity reported on Isla Rasa, Heermann’s Gull showed MHCIIB variability consistent with population expansion, possibly due to a sudden growth following conservation efforts. Duplication must play an essential role in shaping Charadriiformes MHCIIB variability, buffering selective pressures through balancing selection. These findings suggest that MHC copy number and protected islands can contribute to seabird conservation.

HLA-DQ heterodimers in hematopoietic cell transplantation

HLA-DQ heterodimers increase the susceptibility to autoimmune diseases, but their role in hematopoietic cell transplantation is unknown. We tested the hypothesis that outcome after HLA-matched and HLA-DQ-mismatched hematopoietic cell transplantation is influenced by HLA-DQ heterodimers. Heterodimers were defined in 5164 HLA-matched and 520 HLA-DQ-mismatched patients and their transplant donors according to well-established crystallographic criteria. Group 1 (G1) heterodimers are any DQA1*02/03/04/05/06α paired with any DQB1*02/03/04β. Group 2 (G2) heterodimers are DQA1*01α paired with any DQB1*05/06β. Multivariable models identified significantly higher relapse risk in G1G2 and G2G2 compared with G1G1 HLA-matched patients with malignant disease; risk increased with an increasing number of G2 molecules. In HLA-DQ-mismatched transplantation for malignant diseases, matching or mismatching for G2 increased relapse risk. G2 lowered disease-free survival after both HLA-matched and HLA-DQ-mismatched transplantation. A paradigm based on HLA-DQ heterodimers provides a functional definition of the hematopoietic cell transplantation barrier and a means to lower risks for future patients.

Molecular Determinants Regulating the Plasticity of the MHC Class II Immunopeptidome

In the last few years, advancement in the analysis of the MHC class II (MHC-II) ligandome in several mouse and human haplotypes has increased our understanding of the molecular components that regulate the range and selection of the MHC-II presented peptides, from MHC class II molecule polymorphisms to the recognition of different conformers, functional differences in endosomal processing along the endocytic tract, and the interplay between the MHC class II chaperones DM and DO. The sum of all these variables contributes, qualitatively and quantitatively, to the composition of the MHC II ligandome, altogether ensuring that the immunopeptidome landscape is highly sensitive to any changes in the composition of the intra- and extracellular proteome for a comprehensive survey of the microenvironment for MHC II presentation to CD4 T cells.

Diversity of the MHC class II DRB gene in the wolverine (Carnivora: Mustelidae: Gulo gulo) in Finland

The wolverine (Gulo gulo) in Finland has undergone significant population declines in the past. Since major histocompatibility complex (MHC) genes encode proteins involved in pathogen recognition, the diversity of these genes provides insights into the immunological fitness of regional populations. We sequenced 862 amplicons (242 bp) of MHC class II DRB exon 2 from 32 Finnish wolverines and identified 11 functional alleles and three pseudogenes. A molecular phylogenetic analysis indicated trans-species polymorphism, and PAML and MEME analyses indicated positive selection, suggesting that the Finnish wolverine DRB genes have evolved under balancing and positive selection. In contrast to DRB gene analyses in other species, allele frequencies in the Finnish wolverines clearly indicated the existence of two regional subpopulations, congruent with previous studies based on neutral genetic markers. In the Finnish wolverine, rapid population declines in the past have promoted genetic drift, resulting in a lower genetic diversity of DRB loci, including fewer alleles and positively selected sites, than other mustelid species analyzed previously. Our data suggest that the MHC region in the Finnish wolverine population was likely affected by a recent bottleneck.

Long term patterns of association between MHC and helminth burdens in the bank vole support Red Queen dynamics

Major histocompatibility complex (MHC) genes encode proteins crucial for adaptive immunity of vertebrates. Negative frequency-dependent selection (NFDS), resulting from adaptation of parasites to common MHC types, has been hypothesized to maintain high, functionally relevant polymorphism of MHC, but demonstration of this relationship has remained elusive. In particular, differentiation of NFDS from fluctuating selection, resulting from changes in parasite communities in time and space (FS), has proved difficult in short-term studies. Here, we used temporal data, accumulated through long-term monitoring of helminths infecting bank voles (Myodes glareolus), to test specific predictions of NFDS on MHC class II. Data were collected in three, moderately genetically differentiated subpopulations in Poland, which were characterized by some stable spatiotemporal helminth communities but also events indicating introduction of new species and loss of others. We found a complex association between individual MHC diversity and species richness, where intermediate numbers of DRB supertypes correlated with lowest species richness, but the opposite was true for DQB supertypes - arguing against universal selection for immunogenetic optimality. We also showed that particular MHC supertypes explain a portion of the variance in prevalence and abundance of helminths, but this effect was subpopulation-specific, which is consistent with both NFDS and FS. Finally, in line with NFDS, we found that certain helminths that have recently colonized or spread in a given subpopulation, more frequently or intensely infected voles with MHC supertypes that have been common in the recent past. Overall, our results highlight complex spatial and temporal patterns of MHC-parasite associations, the latter being consistent with Red Queen coevolutionary dynamics.

Activation pathways that drive CD4+ T cells to break tolerance in autoimmune diseases. Immunol Rev 2022;307:161-190. [PMID: 35142369 PMCID: PMC9255211 DOI: 10.1111/imr.13071]

Autoimmune diseases are characterized by dysfunctional immune systems that misrecognize self as non-self and cause tissue destruction. Several cell types have been implicated in triggering and sustaining disease. Due to a strong association of major histocompatibility complex II (MHC-II) proteins with various autoimmune diseases, CD4+ T lymphocytes have been thoroughly investigated for their roles in dictating disease course. CD4+ T cell activation is a coordinated process that requires three distinct signals: Signal 1, which is mediated by antigen recognition on MHC-II molecules; Signal 2, which boosts signal 1 in a costimulatory manner; and Signal 3, which helps to differentiate the activated cells into functionally relevant subsets. These signals are disrupted during autoimmunity and prompt CD4+ T cells to break tolerance. Herein, we review our current understanding of how each of the three signals plays a role in three different autoimmune diseases and highlight the genetic polymorphisms that predispose individuals to autoimmunity. We also discuss the drawbacks of existing therapies and how they can be addressed to achieve lasting tolerance in patients.

Using de novo genome assembly and high-throughput sequencing to characterize the MHC region in a non-model bird, the Eurasian coot

Genes of the Major Histocompatibility Complex (MHC) form a key component of vertebrate adaptive immunity, as they code for molecules which bind antigens of intra- and extracellular pathogens (MHC class I and II, respectively) and present them to T cell receptors. In general, MHC genes are hyper-polymorphic and high MHC diversity is often maintained within natural populations (via balancing selection) and within individuals (via gene duplications). Because of its complex architecture with tandems of duplicated genes, characterization of MHC region in non-model vertebrate species still poses a major challenge. Here, we combined de novo genome assembly and high-throughput sequencing to characterize MHC polymorphism in a rallid bird species, the Eurasian coot Fulica atra. An analysis of genome assembly indicated high duplication rate at MHC-I, which was also supported by targeted sequencing of peptide-binding exons (at least five MHC-I loci genotyped). We found high allelic richness at both MHC-I and MHC-II, although signature of diversifying selection and recombination (gene conversion) was much stronger at MHC-II. Our results indicate that Eurasian coot retains extraordinary polymorphism at both MHC classes (when compared to other non-passerine bird species), although they may be subject to different evolutionary mechanism.

High resolution HLA-DRB1 analysis and shared molecular amino acid signature of DRβ1 molecules in Occult hepatitis B infection

AIM

To investigate the association of human leukocyte antigen (HLA)-DRB1 alleles and the variations of polymorphic amino acid changes in DRβ1 chain in Shaanxi Han population with Occult hepatitis B infection (OBI).

METHODS

High-resolution HLA-DRB1 genotyping was performed in 107 OBI carriers and 280 normal controls. Sequence information was used to assign which amino acids were encoded at all polymorphic positions. Three-dimensional modeling was performed to explore the effect of the key residues on the HLA-DRB1 molecule.

RESULTS

Strong susceptible association for allele DRB1*07:01 was observed in OBI carriers. The amino acid variation at HLA-DRβ1 molecule revealed susceptible associations for residues Gln^4β, Val^57β(P9), Ser^60β(P9) and Val^78β(P4), the amino acids Arg^4β, Asp^57β(P9), Tyr^60β(P9) and Tyr^78β(P4) showed protective associations.

CONCLUSION

Alleles DRB1*07:01 showed strong susceptible associations in OBI carriers. The amino acid variations in DRβ molecules revealed significant molecular markers for susceptibility and protection from OBI in Shaanxi Han population.